Malgré leur fort développement en assainissement domestique urbain, les bioréacteurs à membranes (BaM) pâtissent de phénomènes de colmatage, induisant des coûts énergétiques et de maintenance importants. L’aération séquencée des membranes par des grosses bulles est l’une des stratégies pour limiter le colmatage ; son impact a fait l’objet de plusieurs études mais reste globalement mal compris, notamment du fait de la complexité de ces systèmes, multiphasiques et opaques. L’hydrodynamique des réacteurs reste mal caractérisée en présence de boues. Pour apporter des éléments de compréhension aux mécanismes de limitation du colmatage par injection d’air, un pilote de filtration membranaire semi-industriel (2 m3, 3 sous modules fibres creuses Puron®) a été conçu, dimensionné et installé sur l’unité de traitement des jus (TDJ) de la station d’épuration Seine Aval (SIAAP), afin d’être alimenté en boues biologiques dans des conditions réelles de fonctionnement. Le suivi des paramètres opératoires du pilote et de ses performances de filtration sur une période de 5 mois avait un double objectif : (i) mieux caractériser la dispersion du gaz pour différentes conditions de fonctionnement (paramètres de l’aération, concentration en boues de l’alimentation), (ii) hiérarchiser les facteurs qui limitent le colmatage des membranes. Il s’agit in fine de proposer des stratégies d’aération adaptées et efficaces pour limiter le colmatage. Afin de caractériser la dispersion du gaz dans le réacteur, une méthodologie innovante basée sur la tomographie de résistivité électrique (ERT) a été adaptée au pilote. Les conditions d’utilisation de l’ERT (nombre d’électrodes de mesure, séquence de quadripôles) ont été sélectionnées à travers une étude numérique, de même que les paramètres d’inversion nécessaires pour reconstituer la cartographie des résistivités à partir des mesures expérimentales. Cette étude numérique poussée, réalisée sous COMSOL, a permis de conclure à l’intérêt de la méthode pour représenter la distribution des phases dans la géométrie considérée. L’ERT a donc été appliquée au pilote alimenté en boues, pour différentes conditions d’aération. Le jeu de données de filtration a par ailleurs été analysé par logique floue, à l’aide du logiciel FisPro. Les arbres de décision obtenus, en analysant les résultats de manière globale et en les regroupant par conditions opératoires similaires, ont mis en évidence l’impact prépondérant des variables suivantes sur la dérive de perméabilité observée (comprise entre - 9 et 2 LMH/bar) : la différence de DCO entre le surnageant des boues et le perméat (DDCO) traduisant une phase colloïdale complexe, et la concentration en matières en suspension (MES), ayant toutes deux un impact négatif sur les performances de filtration. Une augmentation du débit d’air conduirait à une limitation de la dérive de perméabilité, sauf lorsque la variable DCO est élevée (> à 500 mg/L), cette hypothèse restant cependant à vérifier sur une base de données plus conséquente. Le modèle ainsi obtenu par logique floue permet de mieux simuler les évolutions de perméabilité que les modèles obtenus par régression linéaire multivariée (erreurs de 0,61 et de 0,70 respectivement), et ce malgré une incertitude relative importante sur la mesure de perméabilité (jusqu’à 16 %). Ces résultats sont cohérents avec la dispersion du gaz observée par ERT : son homogénéité dépend de la concentration en MES et du débit d’air injecté. A forte concentration en MES (6 – 10 g/L), des zones préférentielles de passage des bulles ont été observées, en particulier à faible débit d’air, expliquant ainsi un colmatage plus important. L’utilisation nouvelle dans ce contexte de ces techniques, ERT et logique floue, donne des résultats qui confortent l’intérêt d’adapter l’aération (débit, séquençage) aux caractéristiques des boues notamment leurs concentrations, et qui permettent d’envisager des stratégies de contrôle de ces paramètres

Située à l'interface entre les eaux de surface et les eaux souterraines, la zone hyporhéique (ZH) est depuis maintenant plusieurs décennies considérée comme une zone cruciale pour la préservation des milieux aquatiques. Elle constitue souvent un indicateur fiable de la bonne qualité des eaux et une niche écologique primordiale pour de nombreuses espèces. Mais elle est aussi le lieu d'interaction entre deux masses d'eau de signature différente, ce qui conduit à la formation d'un milieu extrêmement fragile et siège d'un grand nombre de réactions biogéochimiques. L'objectif principal de cette thèse est de parvenir à une meilleure compréhension de la dynamique des échanges au sein de la zone hyporhéique. L'approche de cette problématique s'est faite sous un aspect innovant en couplant une démarche hydrogéologique "classique" à l'aide de mesures hydrodynamiques et géochimiques, et l'utilisation de la tomographie de résistivité électrique (ERT). Plusieurs campagnes de terrain ont été menées sur la rivière Essonne, choisie comme lieu d'expérimentation. Différents outils de prélèvement et/ou de mesure ont été mis en place et un grand nombre de mesures à différentes périodes de l'année ont été réalisées. Des expériences assez techniques et innovantes de suivi d'un abaissement et relèvement de barrage, ainsi qu'un traçage artificiel au sel ont pu être effectués grâce à la collaboration avec le syndicat chargé de la gestion et l'aménagement d'une partie du réseau hydrographique de l'Essonne (SIARCE). En parallèle avec cette étude expérimentale, une maquette numérique 3D de la zone d'étude a été réalisée à l'aide du logiciel HydroGeoSphere. Des tests de sensibilité ont permis d'identifier les paramètres hydrodynamiques les plus importants et de quantifier leur impact sur la formation et l'évolution de la zone hyporhéique. Finalement, les premières simulations des expériences menées sur le terrain ont permis de confronter l'approche expérimentale et l'approche théorique. / The Hyporheic Zone (HZ) is located at the interface between surface water and groundwater. For several decades it is considered as a hotspot for the development of a rich aquatic environment in rivers. Its system is often considered as a reliable indicator for water quality and a primary ecological niche for many species. From a hydrological point of view, it is also the place of interaction between two distinct water bodies with different geochemical signatures. This place of mixing forms a very fragile equilibrium where many biogeochemical reactions can occur. The objective of this thesis is to reach a better understanding of mixing and water fluxes in a dynamic context within the hyporheic zone. An innovative method was used by coupling a "classic" hydrogeological approach with hydrodynamic and geochemical measurements with Electrical Resistivity Tomography (ERT). Several field campaigns were done on the Essonne river as experimental site. A large number of measurements were done at various periods of the year and field equipment for water sampling and measurements were installed during these three years. Technical and innovative experiments were conducted such as a dam lowering and rising and an artificial salt tracer test in collaboration with the federation in charge of organization and management of the Essonne network. Finally, a 3D-model of the studied area was built with the HydroGeosphere software. The main hydrodynamic parameters have been tested in order to understand their impact and their variation in a static or dynamic environment on the hyporheic system and its development. In addition, field experiments were reproduced to compare the experimental and theoretical approach.

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.

Malgré leur fort développement en assainissement domestique urbain, les bioréacteurs à membranes (BaM) pâtissent de phénomènes de colmatage, induisant des coûts énergétiques et de maintenance importants. L’aération séquencée des membranes par des grosses bulles est l’une des stratégies pour limiter le colmatage ; son impact a fait l’objet de plusieurs études mais reste globalement mal compris, notamment du fait de la complexité de ces systèmes, multiphasiques et opaques. L’hydrodynamique des réacteurs reste mal caractérisée en présence de boues. Pour apporter des éléments de compréhension aux mécanismes de limitation du colmatage par injection d’air, un pilote de filtration membranaire semi-industriel (2 m3, 3 sous modules fibres creuses Puron®) a été conçu, dimensionné et installé sur l’unité de traitement des jus (TDJ) de la station d’épuration Seine Aval (SIAAP), afin d’être alimenté en boues biologiques dans des conditions réelles de fonctionnement. Le suivi des paramètres opératoires du pilote et de ses performances de filtration sur une période de 5 mois avait un double objectif : (i) mieux caractériser la dispersion du gaz pour différentes conditions de fonctionnement (paramètres de l’aération, concentration en boues de l’alimentation), (ii) hiérarchiser les facteurs qui limitent le colmatage des membranes. Il s’agit in fine de proposer des stratégies d’aération adaptées et efficaces pour limiter le colmatage. Afin de caractériser la dispersion du gaz dans le réacteur, une méthodologie innovante basée sur la tomographie de résistivité électrique (ERT) a été adaptée au pilote. Les conditions d’utilisation de l’ERT (nombre d’électrodes de mesure, séquence de quadripôles) ont été sélectionnées à travers une étude numérique, de même que les paramètres d’inversion nécessaires pour reconstituer la cartographie des résistivités à partir des mesures expérimentales. Cette étude numérique poussée, réalisée sous COMSOL, a permis de conclure à l’intérêt de la méthode pour représenter la distribution des phases dans la géométrie considérée. L’ERT a donc été appliquée au pilote alimenté en boues, pour différentes conditions d’aération. Le jeu de données de filtration a par ailleurs été analysé par logique floue, à l’aide du logiciel FisPro. Les arbres de décision obtenus, en analysant les résultats de manière globale et en les regroupant par conditions opératoires similaires, ont mis en évidence l’impact prépondérant des variables suivantes sur la dérive de perméabilité observée (comprise entre - 9 et 2 LMH/bar) : la différence de DCO entre le surnageant des boues et le perméat (DDCO) traduisant une phase colloïdale complexe, et la concentration en matières en suspension (MES), ayant toutes deux un impact négatif sur les performances de filtration. Une augmentation du débit d’air conduirait à une limitation de la dérive de perméabilité, sauf lorsque la variable DCO est élevée (> à 500 mg/L), cette hypothèse restant cependant à vérifier sur une base de données plus conséquente. Le modèle ainsi obtenu par logique floue permet de mieux simuler les évolutions de perméabilité que les modèles obtenus par régression linéaire multivariée (erreurs de 0,61 et de 0,70 respectivement), et ce malgré une incertitude relative importante sur la mesure de perméabilité (jusqu’à 16 %). Ces résultats sont cohérents avec la dispersion du gaz observée par ERT : son homogénéité dépend de la concentration en MES et du débit d’air injecté. A forte concentration en MES (6 – 10 g/L), des zones préférentielles de passage des bulles ont été observées, en particulier à faible débit d’air, expliquant ainsi un colmatage plus important. L’utilisation nouvelle dans ce contexte de ces techniques, ERT et logique floue, donne des résultats qui confortent l’intérêt d’adapter l’aération (débit, séquençage) aux caractéristiques des boues notamment leurs concentrations, et qui permettent d’envisager des stratégies de contrôle de ces paramètres / Membrane bioreactors (MBRs) are widely used in the wastewater treatment sector. However, membrane fouling mitigation remains challenging, and leads to important maintenance and energy costs. Several strategies have been developed industrially to enhance MBR productivity, including coarse bubble sequenced aeration for fouling mitigation. The way such aeration participates in sludge hydrodynamic patterns is an important research topic. However, the methods currently used for hydrodynamic characterisation suffer from several drawbacks, mainly due to the system’s complexity (three phases with opaque deformable solids). More research is needed to characterize hydrodynamics in MBRs filled with activated sludge. A semi-industrial reactor (2.0 m3) was designed and equipped with three hollow fiber membrane modules (KMS Puron) to contribute to fill this knowledge gap. The reactor was continuously fed with activated sludge from one of the reactors of the Seine Aval (SIAAP) wastewater treatment plant. Operating parameters, activated sludge properties as well as process performances were monitored for five months with two objectives: (i) to characterize gas dispersion for different operating conditions (aeration strategy, activated sludge concentrations), (ii) to rank activated sludge properties and operating conditions, according to their impact on fouling mitigation. The ultimate goal is to propose adapted aeration strategies allowing sustainable filtration performance. In order to characterize the gas dispersion in a complex geometry such as MBR membrane tanks, an innovative approach based on electrical resistivity tomography was adapted to the pilot conditions. A numerical approach was used to define the experimental design in terms of electrode positions, quadripole sequences but also inversion parameters, used to reconstruct resistivity maps from experimental datasets. This numerical study, performed on COMSOL, demonstrated the ability of the ERT method to observe different patterns in the membrane zone. ERT was therefore applied to characterise bubble dispersion in the semi-industrial membrane bioreactor filled with activated sludge and operated using different aeration conditions. The dataset obtained from filtration monitoring was analysed using fuzzy logic. The resulting fuzzy decision trees, constructed from the whole dataset or from subsets clustering similar operating conditions, pointed out the impact on the permeability evolution of two main factors: the COD difference between activated sludge supernatant and permeate (COD), standing for complex colloidal phase, and Mixed Liquor Suspended Solid (MLSS) concentration. Both variables emphasizing fouling at high values. Also, an air flow rate augmentation seemed to mitigate fouling, except at high COD concentrations (> 500 mg/L), where its impact was reversed. This last observation needs further investigation, on a larger dataset. The model obtained from fuzzy logic allows a better simulation of permeability evolutions compared to multivariate linear regression (errors of 0.61 and 0.70, respectively), despite a relatively high measurement uncertainty of permeability (up to 16 %). Those results are consistent with the gas dispersion observed with ERT measurements: homogeneity of gas dispersion depends on MLSS concentration and air flow rate. At high concentrations of MLSS (6 g/L and 10 g/L), a high air flow rate is required for the bubbles to be distributed on the whole membrane zone and preferential flow paths have been observed, especially at low air flow rate. This gas dispersion heterogeneity explains worsened filtration performance. In this context, the new use of these technics, ERT and fuzzy logic, provided results that reinforce the interest to link aeration parameters to sludge properties, their MLSS concentration in particular, leading to consider control strategies for these operating parameters.

Ce travail a pour objectif d’étudier l’impact d’une installation de stockage de déchets nondangereux sur l’eau souterraine. L’ISDND d’Étueffont (90) dont le mode d’exploitation associedes prétraitements mécanique et biologique (broyage et maturation aérobie) comporte un casiernon sécurisé ainsi qu’un casier étanche. La caractérisation des lixiviats sur 21 ans (1989-2010)montre une dégradation rapide des déchets. Un état de stabilisation est atteint en fin de suivi. Lesprétraitements ont permis d’accélérer la dégradation des déchets et de réduire le potentiel polluantde l’ISDND. Les lixiviats sont traités in situ par lagunage naturel. Bien que cette technique soitpeu répandue pour le traitement des lixiviats, elle offre pourtant des abattements significatifs (75-90 %) de la matière organique, des matières en suspension, de l’azote et du phosphore. Lesprécédentes études ayant démontré une contamination de l’eau souterraine par les lixiviats, unedélimitation du panache de contamination a été entreprise en combinant des prospectionsgéophysiques du secteur (tomographie de résistivité électrique) et des analyses physico-Chimiquesde l’eau. Les résultats montrent une infiltration des lixiviats à la base du casier non sécurisé ainsiqu’au travers de la géomembrane du casier sécurisé. Une légère infiltration des eaux de lagunagea aussi été observée au niveau des lagunes. Le panache a une étendue limitée et se concentreessentiellement sous la base de l’ISDND. L’ensemble des résultats a permis de cibler les zonescontaminées et de définir les paramètres physico-Chimiques à surveiller dans le cadre du suivipost-Exploitation auquel est soumise l’ISDND jusqu’en 2031. / The purpose of this work was to evaluate the impact of a landfill on groundwater quality.Municipal solid waste from the Etueffont landfill (Belfort area, France) was mechanically andbiologically treated (shredding and aerobic composting) before being disposed of in an unlinedarea and in a lined cell. The leachate characterization performed over a 21-Year period (1989-2010) showed a rapid degradation of waste. A steady state was reached at the end of themonitoring. Waste pretreatment enhanced the rate of biodegradation and reduced the pollutionpotential of the landfill. The leachate produced was treated on site in a natural lagooning. Whilelagooning is not a widespread practice for leachate treatment, significant removal (75-90 %) oforganic matter, suspended solids, nitrogen and phosphorus was observed. In previous studies,groundwater contamination was highlighted around the landfill. As a consequence, geophysicalinvestigations of the area by electrical resistivity tomography and physico-Chemical analyses wereconducted in order to delineate the contamination plume. The results showed leachate infiltrationinto the soil below the unlined part of the landfill and leakage through the watertightgeomembrane of the lined cell. Seepage of wastewater into the soil below the lagoons was alsoobserved. The extent of contamination plume was limited and was mostly limited to the landfillboundaries. The results of this work allowed to determine the location of the contaminated areasand helped to choose the appropriate physico-Chemical analytes of the post-Closure groundwatermonitoring program, which will be applied to the landfill until the end of 2031.

Le territoire burkinabè est constitué à plus de 80% de formations de socle. Ces formations renferment des ressources en eau souterraine qui constitue la seule source en eau potable de certaines populations. Cependant, le quart des forages réalisés en vue d’exploiter ces eaux sont négatifs. Une bonne compréhension du fonctionnement hydrogéologique et hydrodynamique des aquifères de socle qui sont très hétérogènes peut aider à résoudre ce problème. L’objectif de la thèse est de proposer un modèle conceptuel hydrogéologique en milieu de socle fracturé soudano-sahélien en combinant différentes méthodes hydrogéophysiques et hydrogéologiques classiques. Le site expérimental de Sanon (environ 14 km2) a été choisi comme site d’étude en raison de sa géologie qui est représentative de certaines formations de socle de l’Afrique de l’Ouest. Les résultats indiquent que le profil d’altération de Sanon est cohérent avec d'autres modèles proposés ces dernières années, suggérant que la conductivité hydraulique des aquifères de socle est une conséquence des processus d'altération plutôt que de la fracturation tectonique. Concernant les propriétés hydrodynamiques, une corrélation linéaire positive d’environ 80% a été obtenue entre les transmissivités et les paramètres RMP. La présence d’un dôme piézométrique dans la vallée centrale a été observée en toute saison. La recharge se fait directement au niveau de ce dôme piézométrique (11 et 13% de la pluie annuelle) et l’alimentation des zones périphériques se ferait plutôt par des transferts latéraux. L’ensemble des informations géologiques et hydrogéologiques a permis de faire une simulation numérique qui reproduit l’écoulement observé. / More than 80% of Burkina Faso subsoil is composed of hard rocks. These hard rocks contain ground water resources which are the only drinking water for some people. Yet, the quarter of the boreholes drilled to exploit these waters are negative. A good understanding of the hydrogeological and hydrodynamic functioning of hard rock aquifers which are very heterogeneous can help to solve that problem. The objective of the thesis is to propose a hydrogeological conceptual model of hard rock aquifer under Sudano-Sahelian climate by combining different hydrogeophysical and hydrogeological methods. The Sanon experimental site (14 km² catchment area) was chosen as the study site because of its geology, which is representative of some hard rocks in West Africa. The results show that Sanon weathering profile model is consistent with other models proposed in recent years, suggesting that the hydraulic conductivity of hard rock aquifers is a consequence of weathering processes, rather than tectonic fracturing. Concerning hydrodynamic properties, a positive linear correlation of about 80% was obtained between the transmissivities and the MRS parameters. The presence of a piezometric dome in the central valley on all the piezometric maps produced in both low and high waters was observed. As for recharging, it is done directly in the zone of the piezometric dome (between 11 and 13% of annual rainfall) and the feeding of the peripheral zones would be rather by lateral transfers. The set of geological and hydrogeological data allowed the performance of a numerical simulation in steady state which reproduces well the observed flow.

Le bassin molassique de Carpentras est constitué de formations miocènes sablo-argileuses présentant de nombreuses variations de faciès. Il renferme un réservoir aquifère complexe aux caractéristiques hydrodynamiques intrinsèques médiocres (perméabilité comprise entre 10-4 et 10-7 m.s-1) mais dont la puissance est importante (600 m d'épaisseur). La qualité chimique de l'eau est excellente (naturellement potable) ; l'exploitation est facile et bon marché. Depuis les années 1960, le nombre de forages a augmenté rapidement et sans contrôle. En 2006, on compte des milliers d'ouvrages sur 600 km2. Dans un contexte de rationalisation des usages de l'eau, il est apparu nécessaire de faire un état des connaissances et d'actualiser celles-ci pour construire un schéma de fonctionnement du système aquifère.<br />Ce travail aborde la géologie, l'hydrodynamique et l'hydrochimie du bassin. Chaque élément fait l'objet d'une revue bibliographique critique et est approfondi par l'acquisition et l'analyse de données nouvelles. L'étude géophysique des bordures du bassin par tomographie de résistivité électrique permet de proposer un schéma de l'organisation des corps sédimentaires. Le suivi mensuel du niveau de l'eau dans 51 ouvrages pendant 20 mois révèle l'influence saisonnière très nette des pompages. Enfin deux campagnes de prélèvements et d'analyses des ions caractéristiques des eaux du bassin (100 éch.) ont été menées en 2004 et 2005. La comparaison des résultats avec ceux acquis en 1985 et 1996 permet une approche spatio-temporelle des concentrations (évolution de la contamination nitratée, temps de séjour).<br />On aboutit à un modèle hydrogéologique simple mais complet. La modélisation est alors utilisée comme un outil de synthèse des informations et des hypothèses de fonctionnement. Elle permet ainsi une visualisation générale du fonctionnement de l'aquifère et met en évidence les éléments qui nécessiteraient des investigations supplémentaires pour améliorer la compréhension.

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.

En contexte insulaire volcanique, les aquifères du domaine côtier présentent de très fortes hétérogénéités liées à leur structure, à la disparité de leurs recharges et à l'influence des conditions en aval. Les travaux envisagés ici consistent à élaborer une méthodologie adaptée à différentes configurations hydrogéologiques dans l'optique de mieux imager et comprendre l'évolution du biseau salé. L'acquisition des données est effectuée sur différents sites présentant des configurations hydrogéologiques variées sur le plan de la recharge, de la structure de l'aquifère et de la géomorphologie de la limite en aval. Un réseau d'observation est mis en place sur la côte ouest de l'île de La Réunion, dans la zone littorale allant de L'Ermitage à Saint-Leu. Cette zone est particulièrement intéressante dans le sens où elle présente les types de limites en aval les plus communes aux îles volcaniques à savoir, des côtes rocheuses, des plaines d'arrière-récif et des plaines d'accumulation fluviomarine. Sur La Réunion et la Grande Comore, le milieu volcanique est investigué au moyen de méthodes géophysiques et par une analyse de séries temporelles hydrogéologiques. Les outils géophysiques employés confirment la présence d’un milieu souterrain hétérogène, avec la présence sporadique de chenaux d'écoulement souterrain préférentiel susceptibles de contrôler la dynamique des intrusions marines. Les simulations numériques effectuées par la suite corroborent ces observations. L'analyse des signaux enregistrés indique une forte influence des forçages hydro-climatiques externes sur les nappes côtières. Elle permet aussi de distinguer deux familles de comportement selon le contexte géologique des ouvrages instrumentés. Les résultats de ces travaux mettent en avant une méthodologie de gestion des nappes côtières adaptée au milieu volcanique. Ils contribuent également à l'orientation des programmes de prospection de la ressource en eau souterraine dans ces contextes très complexes. / In volcanic islands context, the aquifers of coastal areas have very strong heterogeneities related to their structure, the disparity of their recharge and the influence of the downstream boundaries (rocky coasts, detritic cones...). The works considered here consist to develop a methodology adapted to different hydrogeological configurations in the aim to better understanding the development of saltwater intrusions. Time series recording and geophysical investigations are carried out on various sites, with different hydrogeological conditions, on the west coast of Reunion Island and throughout Grande Comore Island. The groundwater observatory tested in the pilot study area of Reunion was exported to three others experimental sites situated in Grande Comore Island (Indian Ocean). The geophysical investigations reveal a very important lateral heterogeneity along the shoreline due to the existence of 3D volcanic structures such as palaeovalleys filled by recent lava flows. In coastal areas, these structures exert an important role towards the marine intrusion dynamics, like it was confirmed by numerical simulations. Time series recording has allowed the hydrodynamic characterisation of the coastal groundwater and has enhanced the understanding of low and high frequency variations of the groundwater resources quality

La végétation arborée sur les digues ou les barrages en remblai est un facteur de fragilité pouvant favoriser l'apparition de mécanismes de détérioration par érosion. Définir la structure géométrique des systèmes racinaires ainsi que la nature des sols dans lesquels ils croissent, à partir de méthodes géophysiques non destructives, est nécessaire afin d'appréhender les conséquences de leur développement sur la sécurité d'un ouvrage hydraulique. Des expériences en laboratoire ont d'abord permis de déterminer les propriétés acoustiques et électriques intrinsèques d'échantillons racinaires menant à l'identification de signatures indispensables à la discrimination de l'anomalie liée à la racine sur le terrain. Le montage de supports expérimentaux adaptés nous a conduits à progressivement étudier des paramètres d'influences en conditions contrôlées. Des essais en conditions semi-contrôlées sur un dispositif d'arbres plantés dans un sol homogène ont permis d'évaluer la pertinence des différentes méthodologies d'acquisition, tels que l'utilisation de la polarisation provoquée en tomographie de résistivité complexe ou encore la géométrie des capteurs pour la tomographie acoustique. Des traitements innovants tels que l'analyse en ondelettes sont associés afin d'exploiter la richesse des enregistrements. Les résultats obtenus ont été validés par le relevé des positions réelles des racines. Enfin, des campagnes d'auscultation sur un ouvrage réel ont été réalisées et ont permis de mettre en évidence une variabilité spatiale du corps de digue associée à la présence d'arbres. Une méthodologie adaptée au diagnostic géophysique de la végétation sur les ouvrages a été mise en place. / Woody vegetation from earth dikes or dams is a fragility factor which can promote mechanisms of degradation such as erosion. An accurate assessment of root system structure, from geophysical non-destructive methods, of root position into the embankment (depth, extension), and a good knowledge of soil conditions are critical in order to anticipate the consequences of vegetation development for the hydraulic structure’s safety. Laboratory experiments allowed determining intrinsic acoustical and electrical root properties leading to identify relevant signatures and discriminate anomalies related to roots in the field. The establishment of adapted experimental devices led us progressively to assess different parameters (roots mass, water content. . . ) under controlled conditions. Experiments in semi-controlled conditions with trees planted into a homogenous soil, were conducted to assess the relevance of different methodologies, such as the use of temporal induced polarization in complex resistivity tomography or the geometry of sensors for acoustical tomography. Innovative data processing such as wavelet analysis were used to valorize the rich database. The results were validated by the determination of actual root position.Finally, field investigations into an embankment have been performed to highlight a spatial variability of dike structures associated with trees presence. A methodology adapted to the geophysical diagnostic of vegetation roots in embankments was developed.