Le stockage est la méthode de traitement des déchets non dangereux la plus communément utilisée dans le monde entier car elle est un moyen simple et économique pour leur élimination. Malgré une volonté nationale et européenne pour réduire le stockage de déchets biodégradables, une quantité non négligeable est encore enfouie entrainant la mise en place de modes de gestion spécifique. Autrefois exploitées comme de simples fosses de remplissage, les installations de stockage de déchets non dangereux (ISDND) sont aujourd’hui des ouvrages complexes dont l’objectif est de réduire l’impact environnemental et de valoriser énergétiquement le biogaz produit par la biodégradation de la matière organique. Afin d’accélérer les processus de biodégradation un casier de déchets peut être exploité en mode bioréacteur en réinjectant des lixiviats pour augmenter la teneur en eau des déchets. Le défi des années 2000 était de dimensionner et valider les systèmes de réinjection de lixiviat pour garantir une distribution optimale de la teneur en eau dans les massifs de déchets. Aujourd’hui, la question est de de suivre l’évolution de la biodégradation en tous points d’un massif pour notamment comprendre l’effet de ce mode de gestion. Les méthodes géophysiques en plus d’être non-destructives et spatialisantes sont utilisées depuis des années sur les ISDND et pourraient être sensibles à l’évolution de la biodégradation d’un massif de déchets, comme cela a été démontré pour la biodégradation d’autres milieux. Ainsi la problématique de cette thèse est d’évaluer la capacité de certaines méthodes géophysiques pour suivre l’évolution des paramètres bio-physico-chimiques d’un massif de déchets au cours de sa biodégradation. Un premier travail bibliographique a permis d’identifier quatre méthodes électriques parmi les méthodes géophysiques disponibles pour répondre à cette question :• La résistivité électrique• Le potentiel spontané• La polarisation provoquée• La polarisation provoquée spectraleAprès cet état de l’art, le travail de thèse a été séparé en trois parties. La première a été consacrée à la mise en place du suivi de ces quatre méthodes à l’échelle du laboratoire dans des conditions contrôlées, la seconde à analyser le suivi géophysique long terme sur le site industriel de la SAS Les Champs Jouault et la dernière a évaluée les observations à ces deux échelles. Enfin, la conclusion présente le potentiel de la méthode de mesure du potentiel provoquée comme la plus pertinente pour suivre l’évolution de la biodégradation d’un déchet non dangereux au cours du temps et aborde son utilisation dans un cadre industriel. / Storage is the most commonly used waste treatment method in the world because it is a simple and economical way to dispose of solid waste. Despite a national and European desire to reduce the storage of biodegradable waste, a significant amount is still buried, leading to the implementation of specific management methods. Formerly exploited as mere filling pits, MSWL are today complex structures whose objective is to reduce the environmental impact and energetically valorize the biogas produced by the biodegradation of organic matter. In order to accelerate biodegradation processes, a waste cell can be operated in bioreactor mode by reinjecting leachates to increase the water content of the waste. The challenge of the 2000s was to size and validate leachate re-injection systems to ensure optimal distribution of water content in the waste mass. Today, the question is to monitor evolution of the biodegradation in all points of a waste mass in particular to understand the effect of this management mode. Geophysical methods in addition to being non-destructive and spatializing have been used for years on MSWLs and could be sensitive to the evolution of a waste mass biodegradation, as has been demonstrated for the biodegradation of others environments. Thus the problematic of this thesis is to evaluate the capacity of certain geophysical methods to monitor the evolution of the bio-physicochemical parameters of a waste mass during its biodegradation. A first bibliographic work identified four electrical methods among the geophysical methods available to answer this question:• Electrical resistivity• Self potential• Time domain induced polarization• Spectral induced polarizationAfter this state of the art, the thesis work was separated into three parts. The first one was devoted to the implementation of the monitoring of these four methods at the laboratory scale under controlled conditions, the second to analyze the long-term geophysical monitoring at the industrial site of SAS Les Champs Jouault and the last one to evaluate the observations at these two scales. Finally, the conclusion presents the potential of time domain induced polarization method as the most relevant to monitor the evolution of a waste mass biodegradation over time and discusses its use in an industrial setting.

La caractérisation d'un aquifère est classiquement réalisée à partir d'études géologiques, piézométriques, géochimiques et de pompages d'essai. Ces techniques d'étude présentent certaines limites, qui peuvent néanmoins être réduites par l'utilisation de méthodes complémentaires comme certaines méthodes géophysiques. L'objectif de cette thèse est de quantifier l'amélioration apportée par les informations géophysiques dans la caractérisation hydrogéologique des aquifères. Deux contextes géologiquessont sélectionnés pour les difficultés qu'ils posent à la caractérisation hydrogéologique, l'importance des questions sociétales qui s'y posent et leur pertinence dans le cadre des interventions humanitaires : les aquifères de socle cristallin et les lentilles d'eau douce des cordons sableux côtiers. Les résultats d'une étude numérique et leur validation sur deux sites expérimentaux en Inde montrent que les informations complémentaires apportées par la technique de panneau électrique et la méthode de sondage par résonance magnétique des protons (RMP) permettent de préciser l'interprétation de pompages d'essai en contexte de socle. L'incertitude sur la transmissivité est réduite en moyenne de 39 % et celle sur le coefficient d'emmagasinement de 34 %. Cette approche hydrogéophysique permet ainsi de mieux décrire le milieu souterrain, et de proposer des modèles conceptuels plus complets. En milieu côtier sableux dans le Sud de l'Inde, l'utilisation combinée d'observations piézométriques, de la loi d'Archie et des méthodes de sondage électromagnétiques en domaine temporel (TDEM) et RMP permet de quantifier un volume d'eau douce disponible de 510 litres en moyenne par mètre carré de surface au sol avant la mousson dans la partie inhabitée (contre 420 L/m2 dans la partie habitée), et 670 L/m2 après la mousson (contre 450 L/m2 dans la partie habitée). Ces volumes auraient été surestimés de près de 40 % si une approche basée uniquement sur le TDEM et la loi d'Archie avait été utilisée. Le suivi temporel des mesures géophysiques permet également d'estimer la recharge nette de l'année hydrologique 2008 à 140 mm sur ce cordon (pour 2030 mm de précipitations). De plus, ce cordon sableux présente un champ géomagnétique hétérogène, qui rend la mesure de signaux de Free induction decay, (FID, classiquement enregistrés en sondage RMP) difficile. Une nouvelle procédure de mesure et d'interprétation basée sur des signaux d'écho de spin a donc été testée et validée expérimentalement. La caractérisation du milieu avec ce nouveau protocole permet de préciser la teneur en eau du milieu, sous-estimée de 26 % par le protocole FID. Ce résultat ouvre de nouvelles perspectives dans la recherche d'une relation quantifiée entre la teneur en eau RMP et les différentes porosités hydrogéologiques.

Cette étude a pour objet d’apporter des éléments de compréhension quant à l’amélioration de la résistance à la corrosion des alliages de ZnMg(Al), en comparaison du Zn, en milieu atmosphérique chloré. Le cheminement de réflexion part de considérations fondamentales et tend vers des conditions réelles d’utilisation. La première étape a consisté à synthétiser et caractériser des films de ZnO dopé Mg et notamment l’évolution de la structure cristalline, la structure électronique, la résistivité, (…) avec le dopage. La pertinence de ces paramètres a ensuite été évaluée et discutée au regard de la stabilité de ces films en solution. Dans un second temps, les produits de corrosion formés, en laboratoire, en présence ou non de chlorures, sur des phases pures de Zn et ZnMg, ont été caractérisés. Les processus physico-chimiques liés à la formation de ces produits ont ensuite été discutés, afin de mettre en lumière le rôle du Mg dans l’amélioration de la résistance à la corrosion. La dernière étape, s’est attachée à caractériser des produits de corrosion formés après 6 mois en milieu naturel, en atmosphère marine, sur des phases pures de Zn et ZnMg, et des revêtements industriel de type ZnMgAl. Les résultats obtenus ont permis de mettre en évidence une meilleure stabilité en solution des films de Zn0.84Mg0.16O en comparaison du ZnO, en corrélation avec une présence accrue de liaisons hydroxyles, une augmentation de la résistivité et de la fonction de travail. Ces trois paramètres sont apparus également pertinents, sur les produits de corrosion formés en laboratoire et naturellement en présence de Mg et de Mg, Al, et corrélés à l’amélioration de la résistance à la corrosion. / The aim of this study was to provide some understanding about the improvement of the corrosion resistance of ZnMg and ZnMgAl alloys in comparison to Zn, in chloride content atmospheric environment, from fundamental considerations gradually towards the actual use of these materials. The first approach consisted in the synthesis and physicochemical characterization of films of magnesium doped zinc oxide and particularly the evolution of parameters such as the crystal structure, the electronic structure (band gap), the resistivity, (…) with doping. The relevance of these parameters was then evaluated and discussed regarding to the stability of these films in a solution of pH 8.4. The second stage focused on the characterization of corrosion products formed in the laboratory with or without chloride, on pure Zn and ZnMg. The physico-chemical processes related to the formation of these products were then discussed to highlight the role of Mg in the improvement of the corrosion resistance of ZnMg alloys in comparison with the Zn. The third and final stage consisted in the characterization of corrosion products formed in marine atmosphere on pure Zn and ZnMg, as well as on industrial coating-type ZnMgAl, after 6 months of exposure. The results obtained in this study showed a better stability in solution of solid solution films of Zn0.84Mg0.16O in comparison with the ZnO, correlated with the increase of hydroxyl bonds, associated with an increase of resistivity and work function. These three parameters appeared also relevant regarding the corrosion products formed naturally in the laboratory and in the presence of Mg and Mg-Al, and correlated with improvement of the corrosion resistance.

Au cours de cette thèse, nous étudions la compétition entre les états onde de densité de charge (ODC) et supraconducteur dans le composé multi-bandes à deux chaines TTF[Ni(dmit)2]2. Nous avons réalisé des mesures de résistivité, de pouvoir thermoélectrique et de conductivité thermique sous pression hydrostatique jusqu’aux basses températures. A basse pression, deux ondes de densité de charge détectées par mesures de rayons X et observées par nos mesures de résistivité transverse sont associées aux bandes 1D LUMO et HOMOI de la chaine Ni(dmit)2. A 12kbar, la fusion de ces deux instabilités est associée à l’emboîtement des bandes LUMO avec les bandes HOMOI à travers le point Γ de la première zone de Brillouin. A 18kbar et sur un intervalle de 5-6kbar, on observe un pic de commensurabilité attribué au vecteur d’emboîtement commensurable 2kF = 1/3b*.La supraconductivité est observée à partir de 2kbar et jusqu’à 22kbar avec une température critique 0.6K qui augmente avec la pression et est corrélée à la variation des températures de transition ODC. La supraconductivité est associée à la bande 2D HOMOII du Ni(dmit)2. Les mesures de champ magnétique critique ont permis de donner des informations sur l’évolution de la texture de coexistence métal- ODC en fonction de la pression. Elles mettent aussi en évidence une supraconductivité non conventionnelle avec des nœuds dans le gap.Ces résultats nous ont permis de revisiter le diagramme de phase Température-Pression de ce composé qui s'est révélé beaucoup plus riche que ce qui avait été rapporté jusqu’à maintenant. / In this thesis, we studied the competition between charge density wave (CDW) and superconductivity in the 1D multiband compound TTF[Ni(dmit)2]2. Resistivity, thermoelectric power and thermal conductivity measurements have been performed under pressure and down to very low temperatures.At low pressure, two charge density waves, detected by X ray measurements and observed in our transverse resistivity measurements, are associated to the 1D bands LUMO and HOMO of the Ni(dmit)2 chains. At 12kbar, the fusion of those two CDWs is associated to the nesting of the LUMO bands with the HOMO ones through the Γ point of the Brillouin zone. At 18 and around 5-6kbar, we observe a commensurability pic due to the commensurate wave vector 2kF=1/3b*.From 2 to 22kbar, we observe a superconductive transition at a critical temperature 0.6K which grow with pressure and is correlated to critical temperature of CDW transition. The superconductivity is associated to the 2D band HOMOII of the Ni(dmit)2 chains. Measurements of the upper critical field gave information about the evolution with pressure of metal-charge density wave regions. They bring out the unconventional nature of the superconductivity and the presence of node in the gap.Thanks to those results, we obtained a more rich Temperature – Pressure phase diagram of the TTF[Ni(dmit)2]2 than the existing one.

Les résultats présentés dans ce document portent sur les anisotropies sismique et électrique d'un massif calcaire fracturé dans des galeries parallèles souterraines. Les fractures étant orientées verticalement, leurs propriétés sont approximées par isotropie transverse horizontale (HTI). Plusieurs méthodes d'inversion ont été étudié pour traiter les données sismiques (temps d’arrivée des ondes P et S); et pour traiter les données de résistivité électrique mesuré dans une configuration pôle-pôle (valeur de potentiel mesuré).Pour traiter les données sismiques (quatre campagnes), les méthodes étudiées dans ce document sont: une tomographie isotrope, une approximation de cosinus, une inversion anisotrope basée sur l'algorithme Monte-Carlo pour les paramètres de la matrice de rigidité (de l'isotropie transverse horizontale) et une tomographie anisotrope pour un milieu avec une isotropie transverse inclinée. Toutes les méthodes conduisent à la conclusion qu'il y a effectivement une anisotropie présente dans le massif rocheux et confirment la direction de la vélocité maximale parallèle à la direction de la fracturation. Une forte anisotropie de 15 % est présente dans la zone étudiée.Pour traiter les données électriques (deux campagnes), les méthodes étudiées sont une approximation de cosinus et une inversion anisotrope basée sur l'algorithme Monte-Carlo pour les paramètres de résistivité transversale ρT, de résistivité longitudinale ρL et de l'angle d'orientation du modèle par rapport au système de référence. La présence à 180° de deux maxima de résistivité apparente par rapport à l'azimut de mesure ne peuvent être modélisés par des approches conventionnelles. Une modélisation conceptuelle d'un maillage des fils conducteurs dans un milieu non conducteur a été développé et montre des résultats prometteurs, avec le succès de la modélisation des deux maxima.Des mesures répétées montrent de légères variations de résistivité apparentes et des variations des ondes P. Les variations des ondes S ne sont pas prononcés, ce qui reflète un changement de la saturation en eau. La porosité de 10-15% est estimé à partir des mesures sismiques. / Les résultats présentés dans ce document portent sur les anisotropies sismique et électrique d'un massif calcaire fracturé dans des galeries parallèles souterraines. Les fractures étant orientées verticalement, leurs propriétés sont approximées par isotropie transverse horizontale (HTI). Plusieurs méthodes d'inversion ont été étudié pour traiter les données sismiques (temps d’arrivée des ondes P et S); et pour traiter les données de résistivité électrique mesuré dans une configuration pôle-pôle (valeur de potentiel mesuré).Pour traiter les données sismiques (quatre campagnes), les méthodes étudiées dans ce document sont: une tomographie isotrope, une approximation de cosinus, une inversion anisotrope basée sur l'algorithme Monte-Carlo pour les paramètres de la matrice de rigidité (de l'isotropie transverse horizontale) et une tomographie anisotrope pour un milieu avec une isotropie transverse inclinée. Toutes les méthodes conduisent à la conclusion qu'il y a effectivement une anisotropie présente dans le massif rocheux et confirment la direction de la vélocité maximale parallèle à la direction de la fracturation. Une forte anisotropie de 15 % est présente dans la zone étudiée.Pour traiter les données électriques (deux campagnes), les méthodes étudiées sont une approximation de cosinus et une inversion anisotrope basée sur l'algorithme Monte-Carlo pour les paramètres de résistivité transversale ρT, de résistivité longitudinale ρL et de l'angle d'orientation du modèle par rapport au système de référence. La présence à 180° de deux maxima de résistivité apparente par rapport à l'azimut de mesure ne peuvent être modélisés par des approches conventionnelles. Une modélisation conceptuelle d'un maillage des fils conducteurs dans un milieu non conducteur a été développé et montre des résultats prometteurs, avec le succès de la modélisation des deux maxima.Des mesures répétées montrent de légères variations de résistivité apparentes et des variations des ondes P. Les variations des ondes S ne sont pas prononcés, ce qui reflète un changement de la saturation en eau. La porosité de 10-15% est estimé à partir des mesures sismiques.

Les bassins versants anthropisés sont soumis à des forçages qui évoluent rapidement et déséquilibrent les écosystèmes. Les écosystèmes côtiers et en particulier lagunaires sont très sensibles à la qualité et à la quantité de la source biogène d’origine continentale. La lagune d’Arcachon (Sud Ouest de la France) est alimentée par un bassin versant de 3.000 km2 dont l’occupation des sols est dominée par la forêt depuis plus de 150 ans. La maïsiculture s’est fortement développée depuis les années 1970 et un centre de stockage de déchets ménagers et industriels (CSDU) a été aménagé à moins de 3 km du littoral. Cette thèse se propose d’évaluer la dynamique des éléments biogènes dans le continuum fluvio-estuarien du bassin versant. L’analyse de l’occupation des sols par SIG, une mission d’observation de deux ans, et des études ponctuelles de crues et du continuum fluvial ont permis d’évaluer précisément le flux biogène d’origine continentale et de construire les bases d’un modèle de transport : les taux d’exportation dépendent de l’occupation des sols avec 45, 93, 2.850 et 61 kg N-Nitrate km-2 an-1 calculés pour les parcelles forestières, déforestées, agricoles et urbaines respectivement. L’exportation a lieu préférentiellement en tout début de crue mais le lessivage est rapidement limité. En revanche, les eaux souterraines peuvent contribuer significativement au flux annuel pour les bassins versants les plus agricoles. On a pu démontrer qu’une réduction de 10 % des apports d’engrais azotés réduirait de plus de 50 % le flux d’azote vers la lagune. Une étude approfondie intégrant hydrologie, biogéochimie et géophysique a permis de comprendre les mécanismes de dispersion et d’atténuation d’un panache de pollution issu du CSDU. Celui-ci s’écoule vers la lagune avec l’aquifère superficiel à une vitesse de 23 m an-1 pour l’ammonium et environ 50 m an-1 pour le chlorure mettant en évidence une forte atténuation de l’ammonium dans l’aquifère par sorption et nitrification. Actuellement, le panache n’a pas atteint le littoral et seul le drainage par la rivière permet une exportation rapide jusqu’à la lagune. Ce drainage favorise une très forte nitrification dans la zone hyporhéique et dans la rivière, ce qui limite le flux d’ammonium vers la lagune. Dans ces conditions, le temps de résidence de l’ammonium dans l’aquifère est entre 7 et 18 ans. L’intégration des processus biogéochimiques estuariens dans la dynamique des flux biogènes vers la lagune est une nouveauté. L’étude des éléments biogènes et de la matière organique particulaire permet de démontrer que les estuaires de la lagune sont autotrophes à l’échelle de l’année et permettent un abattement d’environ 30 % du flux de nitrate continental. Saisonnièrement toutefois, les processus de reminéralisation peuvent dominer. Les estuaires ont peu d’impact sur la biogéochimie du phosphore contrairement au sédiment de la lagune. / The anthropized watersheds shift quickly and may disturb the ecosystems. Coastal ecosystems and coastal lagoons are highly sensitive to the amount and the quality of biogenic material originating from continents. The Arcachon lagoon (France South West) raises continental waters originating from a watershed of 3,000 km2 . Pine forestry cover the main surface area of the watershed since 150 years. Cornfields are appeared since 1970 and quickly increase. A solid waste disposal has been created near the littoral. The purpose of this PhD is to assess the global cycle of biogenic materials through the riverine and estuarine continuum. Land-use analyze by GIS, two years of monitoring and some studies of flood events and river continuum have permit to assess the flux of biogenic material and to build a transport model. The export rates depend of land-use and we assess the annual export rate to 45, 93, 2,850, and 61 kg N-Nitrate km-2 year -1 in forested parcel, deforested parcel, agricultural parcel and urban parcel respectively. The export is more efficiency at the beginning of flood events. However the groundwater may supply high amount of nitrogen in highly agricultural watershed. We have shown that a decrease of 10 % in fertilizer input to farming areas would probably reduce the annual flux of more than 50 %. A study based on hydrology, biogeochemistry and geophysics permit to understand the behavior of a contaminated plume originating from the waste disposal. The ammonium plume flows towards the lagoon at a velocity of 23 m year -1 whereas the chloride plume should flows at 50 m year -1. It highlights an efficiency attenuation of ammonium dispersion due to sorption in soil and nitrification. The plume needs about 100 years to reach the lagoon and nowadays the only way of dispersion is the drainage of contaminated groundwater by the river. This way permits an efficiency nitrification in the hyporheic zone and limits the flux of ammonium. Considering that the input of ammonium from the landfill is balanced by the present day output via the river, the residence time of ammonium in the aquifer is between 7 and 18 years. The integration of estuarine biogeochemistry has been realized for the first time. The survey of biogenic materials and organic matter emphasizes that estuaries are autotrophic at annual scale and decrease the nitrogen continental flux of about 30 %. But in autumn, the mineralization and nutrient recycling is observed. In the Arcachon lagoon, the estuaries have a negligible impact on phosphorus biogeochemistry contrary to the sediment.

Ce travail de thèse se compose de deux parties : tout d’abord nous avons caractérisé sur le plan morphologique des motifs de TiOx réalisés par anodisation localisée par AFM (LAO) dans des couches pleines plaques de 5 nm de titane. Nous avons étudié l’influence de la tension d’oxydation, de la vitesse de balayage de la pointe AFM, de l’humidité relative de l’environnement, du mode AFM (contact ou intermittent) et du type de pointe. Les motifs les plus fins atteignent une largeur à mi-hauteur de 21 nm pour 2.2 nm de hauteur, ils sont obtenus avec une pointe PtSi utilisée en mode intermittent sous une tension de polarisation de -7V, une vitesse de balayage de 0.4 µm.s-1 et dans un environnement comportant une humidité de 43%. La deuxième partie de notre travail a été consacrée à l’élaboration et à la caractérisation de jonctions planaires MIM Ti/TiOx/Ti. Ces jonctions sont des motifs TLM de titane comportant une ligne transverse de TiOx réalisée par LAO. Lorsque les jonctions sont stressées électriquement sous air, une transformation morphologique irréversible se produit pour une densité de courant et un champ électrique atteignant de l’ordre de 7.1010 A.m-2 et 3.107 V.m-1 respectivement. Des analyses chimiques et structurales basées sur la microscopie électronique à transmission ont montré que la ligne initiale de TiOx amorphe s’était considérablement élargie et est constituée d’une zone de TiOx cristallin. Cette transformation peut être évitée en appliquant le stress électrique sous vide. Enfin des mesures électriques en température ont permis d’élucider les mécanismes de conduction : émission Schottky sous vide et conduction ionique sous air. / This work is divided in two specific parts: first of all we caracterized oxide patterns made by local anodic oxidation using an AFM on 5nm titanium wafers. We caracterized the morphology of the patterns. We studied the influence of several parameters such as oxidation voltage, writing speed of the AFM tip, relative humidity of the environment, AFM modes (contact or tapping)and the type of the tips we used. Most thinnest pattern we made reaches a full width at half maximum of 21nm with a 2.2nm height. It was obtained using a PtSi coating tip in tapping mode with an oxidation voltage of -7V, a writing speed of 0.4 um/s and a relative humidity of 43%. The second part of our work was dedicated to the realization and the characterization of planar MIM junction Ti/TiOx/Ti. These junctions are TLM patterns with a TiOx line cross-ways over the microwire of the TLM pattern. When the junctions are stressed electrically under ambient atmosphere, an irreversible morphological transformation is happenning for a current density and an electric field of 7.1010 A/m² and 3.107 V/m respectively. Chemical and structural analysis based on transmission electronic microscopy have shown that the initial amorphous TiOx junction have grown importantly with an area of crystalline TiOx. This transformation can be avoided by applying the electric stress under vacuum. Finally, electrical measurements in temperature highlighted the transport mecanisms within the junction: Schottky emission under vacuum and ionic conduction under ambient atmosphere.

Le potentiel rédox et la conductivitée du fluide d'un panache de contamination sont deux paramètres clefs pour la surveillance des sites contaminées. L'objectif de ce travail de thèse est de montrer que les méthodes de résistivité et de potentiel spontané apportent des informations spatialisées sur ces deux paramètres. Ces méthodes permettent ainsi de minimiser et d'optimiser l'implantation de piézomètres de contrôle. L'application de la méthode de résistivité électrique sur le site d'Entressen en Provence (C.E.T. de Classe 2) a permis d'obtenir une carte de la conductivité du fluide de l'aquifère, définissant ainsi l'extension du panache de contamination. Les signaux de potentiel spontané mesurés sur les sites contaminés sont associés à l'écoulement de l'eau dans le sous-sol (phénomène d'électrofiltration) et à la présence de contaminants organiques (phénomène que nous appelons "électro-rédox"). Les études antérieures avaient mis en évidence l'existence d'anomalies de potentiel spontané aux abords de sites contaminés, sans pour autant parvenir à expliquer physiquement ce phénomène. Nous avons donc développé un modèle de (bio)-géobatterie naturelle, dans lequel les bactéries jouent le rôle de transfert d'électrons entre la zone réduite du panache et la zone oxydée. Le courant électrique qui en résulte est à l'origine des anomalies de potentiel électrique mesurées à la surface du sol. Pour valider ce modèle, nous avons réalisé des expériences de terrain (Entressen en Provence et Trecate en Italie) ainsi que des expériences en cuve en présence de bactéries sulfato-réductrices. Les résultats de terrain et de laboratoire sont en accord avec le modèle, montrant un couplage linéaire entre potentiel spontané et potentiel rédox.

Ce travail de thèse porte sur l'étude de la structure et du fonctionnement des systèmes hydrothermaux volcaniques. Nous nous sommes attachés à l'étude géophysique multi-méthodes de deux systèmes hydrothermaux : Waimangu, en Nouvelle-Zélande, et la Solfatara, en Italie. Le premier en raison de sa dynamique, et parce qu'il apparait comme une zone laboratoire purement hydrothermale, sans influence magmatique; le second pour l'observation et la caractérisation géophysique d'un panache hydrothermal accessible à la surface. Dans un premier temps, nous nous sommes intéressés à la reconnaissance de la structure de Waimangu à grande échelle, et à la caractérisation des principales composantes du système (lacs de cratères, zones fumeroliennes et geysers). Cette étude a permis d'associer les méthodes de potentiel spontané et de résistivité électrique avec des mesures ponctuelles de température et de flux de CO2. Nous avons ensuite ciblé notre étude sur des laboratoires à petite-échelle présents sur Waimangu : Iodine vent, Inferno Crater Lake et le site du Old Geyser. Ces systèmes différent les uns des autres par leur taille, mais aussi par leurs caractéristiques physico-chimiques, qui jouent un rôle determinant dans la sensibilité des méthodes géo-électriques. Sur le site de Iodine Pool, l'application de techniques récentes d'imagerie acoustique sous-marine a permis la localisation précise des sources de bruit hydrothermal. Le suivi par polarisation spontanée d'un évent intermittent a montré une cyclicité dans les signaux, qui peut être expliquée par une variation de hauteur piézométrique. L'application conjointe de l'imagerie acoustique et de l'imagerie électrique sur le site du Old Geyser a mis en évidence la complémentarité des deux techniques dans la localisation d'une structure hydrothermale qui correspond à la fois à une source acoustique et électrique. Dans un troisième temps, nous avons étudié la signature géo-électrique des mouvements de fluides qui accompagnent l'activité cyclique du Lac Inferno. Le suivi temporel de résistivité électrique montre que les variations du lac et les mouvements de fluide observés sont caractérisés par des variations importantes de résistivité dans le sous-sol qui pourraient correspondre à des changements de phase du fluide en profondeur, en accord avec des modélisations analogiques. La campagne réalisée à la Solfatare a permis d'imager le panache hydrothermal par les méthodes électriques (tomographie de résistivité en trois dimensions), acoustiques et thermiques. Nous présentons les résultats obtenus avec la tomographie de résistivité électrique et les mesures de température, qui ont permis de déterminer avec précision les limites spatiales du panache hydrothermal.

La caractérisation des circulations thermo-convectives dans les milieux naturels est particulièrement étudiée pour ses applications industrielles, notamment en géothermie. Cette caractérisation correspond à deux enjeux liés à (i) la détection des remontées hydrothermales et (ii) la quantification de l'énergie évacuée par le système actif. La méthodologie proposée par cette thèse comporte ainsi deux volets. Le premier volet vise la localisation des systèmes actifs par caractérisation géophysique de sub-surface. Les techniques de prospection mises en œuvre sont la tomographie de résistivité électrique (TRE), la mesure in situ de température et la mesure de polarisation spontanée (PS). L'objet étudié dans cette thèse est la remontée hydrothermale située au niveau de la limite supérieure de l'effondrement de Rina Grande, au sommet du volcan Stromboli. L'approche géophysique de sub-surface comprend une cartographie spatiale (à l’échelle métrique) en TRE, PS et température et un suivi temporel des signaux PS et thermique. La cartographie spatiale nous renseigne sur l'extension spatiale du système hydrothermal (TRE) et le caractère actif des circulations (PS et température). Le modèle de résistivité électrique a été amélioré en tenant compte de la topographie, issue d'un Modèle Numérique de Terrain généré par photogrammétrie. Le traitement des données de surveillance montre qu'une relation linéaire lie le signal thermique et le signal PS, une relation qui ne peut être purement thermo-électrique. Le jeu de données temporel met l'accent sur la dynamique du système convectif et notamment la probable dépendance de la vigueur du système convectif avec les variations saisonnières de température atmosphérique. Le second volet vise à quantifier l'énergie évacuée par des systèmes thermo-convectifs à l'aide d'une approche numérique. L'approche géophysique met en évidence le fait que les zones endommagées telles que des limites structurales favorisent la circulation des fluides hydrothermaux. Ceci s'explique par le fait que fractures, failles ouvertes, zones d'altération etc. changent à petite échelle la perméabilité du milieu. L'approche numérique est donc vouée à la simulation 3D de l'écoulement et du transfert de chaleur dans des milieux poreux fracturés. Le travail s'est basé sur un code d'écoulement en milieux poreux fracturé pré-existant. Ce code a été adapté afin de résoudre le problème thermique. L'étude théorique vise à quantifier l'influence des paramètres du réseau de fractures sur l'énergie dissipée. L'influence de la densité de fractures, de la transmissivité de fractures et l'anisotropie du réseau de fractures sont évaluées. L'étude compare également les flux de chaleur dégagés par des modèles où les fractures sont insérées de manière discrète avec des modèles homogènes et de même propriétés macroscopiques (approche effective). Les résultats montrent que la validité de l'approche effective est fortement dépendante des paramètres du réseau de fracture. / The characterization of thermo-convective circulations in natural environments is particularly studied especially for geothermal applications. This characterization corresponds to two issues related to (i) the detection of the hydrothermal flows(ii) the quantification of the energy discharged by the active system. The methodology proposed by this thesis thus comprises two parts. The first component aims at locating active systems with geophysical methods. The prospecting techniques used are electrical resistivity tomography (ERT), in situ temperature measurement and self potential (SP) measurements. The object studied in this thesis is the hydrothermal system located at the upper limit of the Rina Grande sector collapse, at the top of the Stromboli volcano. The sub-surface geophysical approach includes spatial (metric) mapping in ERT, SP and temperature, and monitoring of SP and thermal signals. Spatial mapping informs us about the spatial extension of the hydrothermal system (ERT) and the active character of the circulations (SP and temperature). The model of electrical resistivity has been improved taking into account the topography, resulting from a Digital Elevation Model generated by photogrammetry. The treatment of the monitoring data shows that a linear relationship links the thermal and the PS signals, a relationship that can not be entirely due to thermoelectric effect. The temporal dataset focuses on the dynamics of the convective system and in particular the probable dependence of the vigor of the convective system on seasonal variations in atmospheric temperature.The second part aims at quantifying the energy released by thermo-convective systemsusing a numerical approach. The geophysical approach highlights the fact that damaged areas such as structural boundaries favor the circulation of hydrothermal fluids. This is due to the fact that fractures, open faults, alteration zones, etc. change the permeability of the medium on a small scale. The numerical approach is therefore dedicated to the 3D simulation of the flow and heat transfer in fractured porous media. The work is based on a pre-existing code in fracturing porous media, able to solve the flow problem. This code has been adapted to solve the thermal problem also. The theoretical study aims to quantify the influence of fracture network parameters on the released energy. The influence of fracture density, fracture transmissivity and fracture network anisotropy are evaluated. The study also compares the heat fluxes generated by models where the fractures are discretely inserted with homogeneous models with the same macroscopic properties (effective approach). The results show that the validity of the effective approach is highly dependent on the parameters of the fracture network.