Les Champs Phlégréens, situés dans la métropole napolitaine (Italie du sud), forment l’une des plus grandes structures volcaniques au monde. Depuis 1950, ce complexe volcanique manifeste un regain d’activité, qui s’est amplifié au cours de la dernière décennie. Cette accélération s’exprime au travers d’une intensification de la sismicité, de la déformation du sol ainsi qu’une extension de la zone de dégazage. L’ensemble des récentes études s’accorde à dire que le système s’achemine actuellement vers un point critique, sans toutefois pouvoir préciser quand et où pourrait avoir lieu une éventuelle éruption. Cette difficulté à prédire l’état réel du système est principalement associée à la présence d’un système hydrothermal relativement développé. Aux Champs Phlégréens, il est en effet difficile de déconvoluer les signaux provenant du forçage magmatique de ceux résultant de la réponse hydrothermale. L’objectif de cette thèse est donc d’améliorer les connaissances actuelles du système hydrothermal superficiel du volcan de la Solfatara, lieu où se concentre actuellement la reprise d’activité. Pour cela, une approche multidisciplinaire a été menée en deux phases : l’imagerie géophysique du volcan puis la modélisation de son système hydrothermal.La tomographie haute-résolution de résistivité électrique 3-D du cratère a permis de reconnaître les principales formations géologiques et leurs connexions avec les structures et écoulements hydrothermaux. L’interprétation du modèle de résistivité électrique a été réalisée grâce à un ensemble de mesures superficielles complémentaires : flux de CO2, température, potentiel spontané, capacité d’échange cationique et pH du sol. Deux panaches à dominante liquide ont été identifiés : la mare de boue de la Fangaia et la fumerole de Pisciarelli. À la Fangaia, une étude conjointe des modèles de résistivité électrique et de vitesses du sous-sol (obtenues par l’INGV) établit la présence de forts gradients, à la frontière entre panache hydrothermal et zone de dégazage diffus. Au niveau du principal secteur fumerolien, le modèle de résistivité électrique et la localisation des sources acoustiques révèlent clairement l’anatomie d’une zone fumerolienne. Deux conduits séparés, saturés en gaz, alimentent les fumeroles de Bocca Grande et de Bocca Nuova, depuis un même réservoir de gaz situé à ~50 mètres de profondeur. L’intense dégazage diffus produit à proximité de ces fumeroles occasionne la condensation de vapeur. Le modèle de résistivité électrique met en évidence la circulation souterraine de cet important volume d’eau, canalisée à l’intérieur d’une zone fracturée.En utilisant l’ensemble de ces informations structurelles, un modèle thermodynamique des écoulements multiphasiques de la principale zone fumerolienne a été réalisé. Ce modèle reproduit fidèlement les observables des fumeroles : température, flux et rapport CO2/H2O. Il valide l’imagerie géophysique et confirme l’interaction entre la circulation d’eau de condensation et l’un des conduits fumeroliens. Ainsi, cette simulation explique, pour la première fois par un effet d’interaction superficiel, les différentes signatures géochimiques des deux fumeroles : Bocca Nuova et Bocca Grande. L’approche multidisciplinaire, employée dans cette thèse, constitue une nouvelle étape vers une meilleure connaissance des interactions hydrothermales. Celles-ci doivent être prise en compte dans l’objectif de réaliser des modélisations dynamiques précises permettant d’appréhender in fine l’état réel du système volcanique. / The Campi Flegrei caldera is located in the metropolitan area of Naples (Italy), and it is one of the largest volcanic systems on Earth. Since 1950, this volcanic complex shows significant unrest, which accelerated over the last decade with a rise in the seismic activity, ground deformation, and the extent of the degassing area. Recent studies indicate that the volcanic system is potentially moving toward a critical state, although their authors remain unable to point out when and where a possible eruption could take place. The difficulty of predicting the real volcanic state is here mainly related to the hydrothermal system. Indeed, at the Campi Flegrei, it is difficult to separate the magmatic input signal from the hydrothermal response. Hence, the aim of this thesis is to improve our knowledge on the shallow hydrothermal system of the Solfatara volcano, where most of the renewal activity takes place. A multidisciplinary approach has been performed in two steps: first a geophysical imagery of the volcano and second the modeling of its hydrothermal system.The 3-D electrical resistivity tomography of the crater allows to recognize the main geological units, and their connection with hydrothermal fluid flow features. The interpretation of the resistivity model has been realized thanks to numerous soil complementary measurements: CO2 flux, temperature, self-potential, Cation Exchange Capacity and pH. We identify two liquid-dominated plumes: the Fangaia mud pool and the Pisciarelli fumarole. In the Fangaia area, the comparison between electrical resistivity and velocity models reveals strong gradients related to a sharp transition at the border between the hydrothermal plume and the high diffuse degassing region. Combining electrical resistivity model with hydrothermal tremor sources localization reveal the anatomy of the main fumarolic area. Two separated conduits, gas-saturated, feed the two fumaroles Bocca Grande and Bocca Nuova. These conduits originate from the same gas reservoir located 60 m below the surface. The intense degassing activity, produced in the vicinity of fumaroles, creates large amounts of vapor condensation. The resistivity model reveals this condensate circulation, within a fractured area.All these results are incorporated into a multiphase flow model of the main fumarolic area. The simulation accurately reproduces the fumaroles observables: temperature, flux and CO2/H2O ratio. The model validates the geophysical imagery and confirms the interaction between Bocca Nuova fumarolic conduit and the condensate flow. Hence, this simulation explains for the first time the distinct geochemical signature of the two fumaroles due to a shallow water-interaction. The multidisciplinary approach performed in this thesis constitutes a new step toward a better understanding of hydrothermal interactions. Those phenomena have to be taken into account in order to perform dynamic modelling, and thus apprehend the real state of the volcanic system.

Bei oberflächennahen geothermischen Anwendungen kommt es zu Temperaturveränderungen im Untergrund, um Wärme zu gewinnen oder zu speichern. Die dadurch verursachten thermischen Auswirkungen können die thermischen, hydraulischen und chemischen Bedingungen des Untergrunds verändern und sich auf die Systemleistung auswirken. Die Situation wird noch komplizierter, wenn die Grundwasserleiter unterschiedlich gesättigt sind und eine Verunreinigung des Untergrunds vorliegt. In dieser Arbeit werden vollständig gekoppelte Wärme- und Stofftransportmodelle angewandt, um die Wechselwirkung zwischen oberflächennahen geothermischen Anwendungen und der unterirdischen Umgebung besser zu charakterisieren. Es werden drei verschiedene Szenarien untersucht. Zunächst werden die thermischen Auswirkungen und die Nachhaltigkeit einer intensiven oberflächennahen geothermischen Nutzung in einem Wohngebiet über einen Zeitraum von 24 Jahren bewertet. Überwachungsdaten des Standorts werden in das Modell integriert, und die Bedeutung standortspezifscher Kenntnisse für die Planung wird hervorgehoben. Zweitens wird der gekoppelte Feuchtigkeits- und Wärmetransport eines Erdwärmespeichersystems (BTES) untersucht. Es wird eine Sensitivitätsanalyse durchgeführt und die Wärmeentzugsefzienz zwischen verschiedenen Szenarien verglichen. Bei Austrocknung oder starker Gasphasenkonvektion werden signifkante Änderungen im Wärmetransportverhalten erwartet. Drittens wird ein nichtisothermes Drei-Komponenten-Zweiphasenströmungsmodell abgeleitet und mit der Finite-Elemente-Methode implementiert. Die Validierung des numerischen Modells bestätigt seine Fähigkeit, die gekoppelte Strömung sowie den Wärme- und Stofftransport in einem ungespannten Grundwasserleiter zu simulieren, der einem BTES-Betrieb ausgesetzt ist. Insbesondere wird das Potenzial für die thermisch verstärkte Verflüchtigung von Trichlorethylen (TCE) in wässriger Phase bewertet. Da durch die Wärmezufuhr eine Auftriebsströmung induziert wird, nimmt die Verringerung der Schadstoffmasse erheblich zu und erreicht nach fünf Jahren mehr als 70 %. Die in dieser Arbeit gewonnenen Erkenntnisse können durch die Verbesserung ihrer ökologischen und ökonomischen Leistungen zu einer breiteren Akzeptanz von Technologien der flachen Geothermie beitragen. / Shallow geothermal applications induce temperature changes in the subsurface for heat extraction or storage purposes. Their induced thermal impacts may alter the thermal, hydraulic, and chemical conditions of the subsurface and feedback on the system performance. The situation is further complicated in variably saturated aquifers and when subsurface contamination exists. In this work, fully coupled heat and mass transport models are applied to improve the characterization of interaction between shallow geothermal applications and the subsurface environment. Three different scenarios are investigated. First, the thermal impact and sustainability of intensive shallow geothermal exploitation in a residential area is evaluated over a 24-year period. Monitoring data from the site is integrated into the model and the importance of site-specifc knowledge for planning is highlighted. Second, the coupled moisture and heat transfer of a soil borehole thermal energy storage (BTES) system is explored. A sensitivity analysis is performed and the heat extraction efciency is compared among various scenarios. Signifcant changes in the heat transport behavior are expected when drying out or strong gas phase convection occurs. Third, a non-isothermal three-component two-phase flow model is derived and implemented with the fnite element method. Validation of the numerical model confrms its ability to simulate the coupled flow, heat and mass transport in an unconfned aquifer subject to BTES operation. In particular, the potential for thermally-enhanced volatilization of aqueousphase trichloroethylene (TCE) is assessed. As buoyant flow is induced due to heat injection, reduction of the contaminant mass grows considerably, reaching more than 70% after fve years. The fndings obtained from this thesis can contribute to a wider adoption of shallow geothermal technologies through the enhancement of their environmental and economical services.