L'influence d'un système hydrothermal au sein d'un édifice volcanique est un concept encore peu étudiée au regard de son importance pour la stabilité d'un volcan. Ce travail s'intéresse à la relation entre système hydrothermal et tectonique à travers l'étude des déformations engendrées dans un cône volcanique par le fluage de roches affectées par l'altération hydrothermale. La présence de ces argiles, ductiles à l'échelle de temps considérée, modifie profondément le comportement d'un édifice volcanique au cours de son évolution, pouvant conduire à la formation de grandes structures d'effondrement de type caldera. La confrontation de modèles analogiques avec la géologie et les données géophysiques nous permet de contraindre dans ce sens l'interprétation des structures d'effondrement observées par exemple sur Nuku Hiva (Polynésie française) et sur le Piton de la Fournaise (Ile de la Réunion). La modélisation analogique montre également que la formation du plancher plat d'une caldera peut être le résultat d'un processus syneffondrement, ne nécessitant pas de resurfaçage postérieur par érosion ou mise en place de nouveaux produits éruptifs. Les expériences nous permettent d'identifier la pente du cône et la diminution des contraintes latérales sur l'édifice comme les paramètres prépondérants à l'origine de la déformation. Les variations de ces paramètres sont déterminantes pour la morphologie de la caldera.

Un des enjeux de la volcanologie réside dans la compréhension du fonctionnement du système de stockage et de transferts magmatiques associés aux éruptions. Au Piton de La Fournaise, cette étude est ici abordée par une analyse croisée des déformations, couplée à des modélisations numériques, et aux données géochimiques de la période 1998-2007. Au sein de cette activité, nous définissons, depuis 2000, un caractère cyclique ; chaque cycle est défini par une séquence d'éruptions sommitales et latérales terminées par une éruption distale à océanite. Ces éruptions successives sont caractérisées par des dykes s'initiant à des niveaux de plus en plus profonds et des laves de plus en plus primitives. Ces cycles durant lesquels l'inflation du cône est continue, pourraient témoigner d'une réalimentation continue du système de stockage superficiel. Ces résultats permettront dans le futur une meilleure prédiction de la localisation des éruptions, selon leur position dans un cycle d'activité.

L'éruption d'avril 2007 du Piton de la Fournaise (Île de la Réunion, Océan Indien) a été marquée par les plus gros volumes de lave émis de ces deux derniers siècles, ainsi que par l'effondrement du cratère sommital (le Dolomieu) sur plus de 300 mètres de haut. Des données d'interférométrie radar (InSAR) montrent que les déplacements associés à cette éruption sont inhabituels pour deux raisons : (1) ils ont affecté l'ensemble de l'enclos Fouqué durant l'éruption, (2) deux motifs de déformation ont persisté plus d'un an après la fin de l'éruption. Le premier résulte d'une subsidence centripète du cône central et le second d'un glissement vers l'est du flanc est du volcan. Une méthode, appelée tomographie de déplacements, a été développée pour modéliser les déplacements, basée sur une discrétisation du sous-sol en sources unitaires, et sur la minimisation de deux fonctions coût. Elle permet de trouver une répartition compacte des variations de volume des sources unitaires permettant de reproduire au mieux les déplacements observés. En parallèle, la procédure NA-MBEM, une méthode de modélisation basée sur la combinaison d'un modèle numérique (MBEM) et d'une inversion de type Monte Carlo (NA), a été modifiée pour diminuer le temps de calcul nécessaire à l'obtention d'un bon modèle, et nous montrons que l'utilisation de données temporellement interpolées permet d'améliorer les résultats d'inversion. L'application des deux méthodes de modélisation aux données de déplacements de l'éruption d'avril 2007, montre que durant la période post-éruptive, la subsidence du cône central est provoquée par une source localisée de manière superficielle sous le cône central. Cette source est interprétée comme un système hydrothermal en déflation, dont le drainage aurait été amorcé par l'effondrement du Dolomieu. Les déplacements du flanc est seraient dus à deux sources différentes, toutes les deux superficielles et parallèles à la topographie. L'une est interprétée comme un réservoir temporaire en cours de vidange durant la fin de l' éruption, et l'autre comme un niveau de glissement sur lequel glisserait le flanc est. Enfin, nous proposons un modèle préliminaire des déplacements ayant eu lieu durant l'éruption, ainsi qu'un scénario de la succession de tous les événements à l'origine des déformations enregistrées entre le 30 mars 2007 et juin 2008.

Les volcans boucliers basaltiques se déforment sous leur propre poids, tout particulièrement en période d'éruption. A Hawaii, cette déformation co-éruptive est en majeure partie expliquée par un plan de décollement à la base des édifices, dont le glissement serait associé à la force de poussée des intrusions dans les " rift zones " (i.e. concentration de dykes subverticaux dans les zones de faiblesse de l'édifice). Cependant, cette explication n'est pas valable pour beaucoup d'autres volcans basaltiques, où l'existence d'un décollement est peu probable. Nous avons profité de l'incision intense du Piton des Neiges (La Réunion) par l'érosion pour observer la structure interne d'un volcan basaltique, et comprendre, par le terrain et la modélisation numérique, comment la déformation d'un tel édifice s'articule avec son activité magmatique. Notre étude structurale révèle que trois populations d'intrusions perpendiculaires (deux subverticales et une subhorizontale) coexistent, et que les principales zones de faiblesses sont des " sill zones " (i.e. concentrations d'intrusions subhorizontales), Parallèlement, notre étude des déformations, essentiellement cassantes, révèle que deux extensions perpendiculaires dominent dans l'édifice, avec l'apparition de régimes localement décrochants et même compressifs à proximité des sill zones. Les directions de déformations et les directions d'intrusions sont compatibles et coïncident aussi avec la direction d'écoulement des grandes avalanches de débris qui ont affectée le Piton des Neiges au cours de son histoire. A partir de ces données de terrain, nous proposons un modèle de déformation du Piton des Neiges où les intrusions dans l'édifice génèrent des cycles de permutations de contraintes. Le stade ultime de chaque cycle serait la mise en place de sills en régime compressif. En testant numériquement l'effet de telles injections sur la déformation d'un édifice, nous montrons que les sills peuvent être un facteur majeur d'instabilité, capable de conduire à des déstabilisations de flanc telles que les avalanches observées. Notre modèle conceptuel, déduit du terrain et quantifié par la modélisation, constitue ainsi une alternative au modèle de déformation hawaiien, applicable sur des édifices tels que le Piton de la Fournaise (La Réunion), Tenerife (Canaries) et Fernandina (Galapagos). Notre étude démontre enfin l'intérêt essentiel d'étudier les systèmes anciens pour mieux comprendre les volcans actifs.