Nucleation and growth of metamorphic minerals in contact aureoles

THESIS ABSTRACT Nucleation and growth of metamorphic minerals are the consequence of changing P-T-X-conditions. The thesis presented here focuses on processes governing nucleation and growth of minerals in contact metamorphic environments using a combination of geochemical analytics (chemical-, isotope-, and trace element composition), statistical treatments of spatial data, and numerical models. It is shown, that a combination of textural modeling and stable isotope analysis allows a distinction between several possible reaction paths for olivine growth in a siliceous dolomite contact aureole. It is suggested that olivine forms directly from dolomite and quartz. The formation of olivine from this metastable reaction implies metamorphic crystallization far from equilibrium. As a major consequence, the spatial distribution of metamorphic mineral assemblages in a contact aureole cannot be interpreted as a proxy for the temporal evolution of a single rock specimen, because each rock undergoes a different reaction path, depending on temperature, heating rate, and fluid-infiltration rate. A detailed calcite-dolomite thermometry study was initiated on multiple scales ranging from aureole scale to the size of individual crystals. Quantitative forward models were developed to evaluate the effect of growth zoning, volume diffusion and the formation of submicroscopic exsolution lamellae (<1 µm) on the measured Mg-distribution in individual calcite crystals and compare the modeling results to field data. This study concludes that Mg-distributions in calcite grains of the Ubehebe Peak contact aureole are the consequence of rapid crystal growth in combination with diffusion and exsolution. The crystallization history of a rock is recorded in the chemical composition, the size and the distribution of its minerals. Near the Cima Uzza summit, located in the southern Adamello massif (Italy), contact metamorphic brucite bearing dolomite marbles are exposed as xenoliths surrounded by mafic intrusive rocks. Brucite is formed retrograde pseudomorphing spherical periclase crystals. Crystal size distributions (CSD's) of brucite pseudomorphs are presented for two profiles and combined with geochemistry data and petrological information. Textural analyses are combined with geochemistry data in a qualitative model that describes the formation periclase. As a major outcome, this expands the potential use of CSD's to systems of mineral formation driven by fluid-infiltration. RESUME DE LA THESE La nucléation et la croissance des minéraux métamorphiques sont la conséquence de changements des conditions de pression, température et composition chimique du système (PT-X). Cette thèse s'intéresse aux processus gouvernant la nucléation et la croissance des minéraux au cours d'un épisode de métamorphisme de contact, en utilisant la géochimie analytique (composition chimique, isotopique et en éléments traces), le traitement statistique des données spatiales et la modélisation numérique. Il est montré que la combinaison d'un modèle textural avec des analyses en isotopes stables permet de distinguer plusieurs chemins de réactions possibles conduisant à la croissance de l'olivine dans une auréole de contact riche en Silice et dolomite. Il est suggéré que l'olivine se forme directement à partir de la dolomie et du quartz. Cette réaction métastable de formation de l'olivine implique une cristallisation métamorphique loin de l'équilibre. La principale conséquence est que la distribution spatiale des assemblages de minéraux métamorphiques dans une auréole de contact ne peut pas être considérée comme un témoin de l'évolution temporelle d'un type de roche donné, puisque chaque type de roche suit différents chemins de réactions, en fonction de la température, la vitesse de réchauffement et le taux d'infiltration du fluide. Une étude thermométrique calcite-dolomite détaillée a été réalisée à diverses échelles, depuis l'échelle de l'auréole de contact jusqu'à l'échelle du cristal. Des modèles numériques quantitatifs ont été développés pour évaluer l'effet des zonations de croissance, de la diffusion volumique et de la formation de lamelles d'exsolution submicroscopiques (<1µm) sur la distribution du magnésium mesuré dans des cristaux de calcite individuels. Les résultats de ce modèle ont été comparés ä des échantillons naturels. Cette étude montre que la distribution du Mg dans les grains de calcite de l'auréole de contact de l'Ubehebe Peak (USA) résulte d'une croissance cristalline rapide, associée aux processus de diffusion et d'exsolution. L'histoire de cristallisation d'une roche est enregistrée dans la composition chimique, la taille et la distribution de ses minéraux. Près du sommet Cima Uzza situé au sud du massif d'Adamello (Italie), des marbres dolomitiques à brucite du métamorphisme de contact forment des xénolithes dans une intrusion mafique. La brucite constitue des pseudomorphes rétrogrades du périclase. Les distributions de taille des cristaux (CSD) des pseudomorphes de brucite sont présentées pour deux profiles et sont combinées aux données géochimiques et pétrologiques. Les analyses textorales sont combinées aux données géochimiques dans un modèle qualitatif qui décrit la formation du périclase. Ceci élargit l'utilisation potentielle de la C5D aux systèmes de formation de minéraux controlés par les infiltrations fluides. THESIS ABSTRACT (GENERAL PUBLIC) Rock textures are essentially the result of a complex interaction of nucleation, growth and deformation as a function of changing physical conditions such as pressure and temperature. Igneous and metamorphic textures are especially attractive to study the different mechanisms of texture formation since most of the parameters like pressure-temperature-paths are quite well known for a variety of geological settings. The fact that textures are supposed to record the crystallization history of a rock traditionally allowed them to be used for geothermobarometry or dating. During the last decades the focus of metamorphic petrology changed from a static point of view, i.e. the representation of a texture as one single point in the petrogenetic grid towards a more dynamic view, where multiple metamorphic processes govern the texture formation, including non-equilibrium processes. This thesis tries to advance our understanding on the processes governing nucleation and growth of minerals in contact metamorphic environments and their dynamic interplay by using a combination of geochemical analyses (chemical-, isotope-, and trace element composition), statistical treatments of spatial data and numerical models. In a first part the thesis describes the formation of metamorphic olivine porphyroblast in the Ubehebe Peak contact aureole (USA). It is shown that not the commonly assumed succession of equilibrium reactions along a T-t-path formed the textures present in the rocks today, but rather the presence of a meta-stable reaction is responsible for forming the olivine porphyroblast. Consequently, the spatial distribution of metamorphic minerals within a contact aureole can no longer be regarded as a proxy for the temporal evolution of a single rock sample. Metamorphic peak temperatures for samples of the Ubehebe Peak contact aureole were determined using calcite-dolomite. This geothermometer is based on the temperature-dependent exchange of Mg between calcite and dolomite. The purpose of the second part of this thesis was to explain the interfering systematic scatter of measured Mg-content on different scales and thus to clarify the interpretation of metamorphic temperatures recorded in carbonates. Numerical quantitative forward models are used to evaluate the effect of several processes on the distribution of magnesium in individual calcite crystals and the modeling results were then compared to measured field. Information about the crystallization history is not only recorded in the chemical composition of grains, like isotope composition or mineral zoning. Crystal size distributions (CSD's) provide essential information about the complex interaction of nucleation and growth of minerals. CSD's of brucite pseudomorphs formed retrograde after periclase of the southern Adamello massif (Italy) are presented. A combination of the textural 3D-information with geochemistry data is then used to evaluate reaction kinetics and to constrain the actual reaction mechanism for the formation of periclase. The reaction is shown to be the consequence of the infiltration of a limited amount of a fluid phase at high temperatures. The composition of this fluid phase is in large disequilibrium with the rest of the rock resulting in very fast reaction rates. RESUME DE LA THESE POUR LE GRAND PUBLIC: La texture d'une roche résulte de l'interaction complexe entre les processus de nucléation, croissance et déformation, en fonction des variations de conditions physiques telles que la pression et la température. Les textures ignées et métamorphiques présentent un intérêt particulier pour l'étude des différents mécanismes à l'origine de ces textures, puisque la plupart des paramètres comme les chemin pression-température sont relativement bien contraints dans la plupart des environnements géologiques. Le fait que les textures soient supposées enregistrer l'histoire de cristallisation des roches permet leur utilisation pour la datation et la géothermobarométrie. Durant les dernières décennies, la recherche en pétrologie métamorphique a évolué depuis une visualisation statique, c'est-à-dire qu'une texture donnée correspondait à un point unique de la grille pétrogénétique, jusqu'à une visualisation plus dynamique, où les multiples processus métamorphiques qui gouvernent 1a formation d'une texture incluent des processus hors équilibre. Cette thèse a pour but d'améliorer les connaissances actuelles sur les processus gouvernant la nucléation et la croissance des minéraux lors d'épisodes de métamorphisme de contact et l'interaction dynamique existant entre nucléation et croissance. Pour cela, les analyses géochimiques (compositions chimiques en éléments majeurs et traces et composition isotopique), le traitement statistique des données spatiales et la modélisation numérique ont été combinés. Dans la première partie, cette thèse décrit la formation de porphyroblastes d'olivine métamorphique dans l'auréole de contact de l'Ubehebe Peak (USA). Il est montré que la succession généralement admise des réactions d'équilibre le long d'un chemin T-t ne peut pas expliquer les textures présentes dans les roches aujourd'hui. Cette thèse montre qu'il s'agirait plutôt d'une réaction métastable qui soit responsable de la formation des porphyroblastes d'olivine. En conséquence, la distribution spatiale des minéraux métamorphiques dans l'auréole de contact ne peut plus être interprétée comme le témoin de l'évolution temporelle d'un échantillon unique de roche. Les pics de température des échantillons de l'auréole de contact de l'Ubehebe Peak ont été déterminés grâce au géothermomètre calcite-dolomite. Celui-ci est basé sur l'échange du magnésium entre la calcite et la dolomite, qui est fonction de la température. Le but de la deuxième partie de cette thèse est d'expliquer la dispersion systématique de la composition en magnésium à différentes échelles, et ainsi d'améliorer l'interprétation des températures du métamorphisme enregistrées dans les carbonates. Des modèles numériques quantitatifs ont permis d'évaluer le rôle de différents processus sur la distribution du magnésium dans des cristaux de calcite individuels. Les résultats des modèles ont été comparés aux échantillons naturels. La composition chimique des grains, comme la composition isotopique ou la zonation minérale, n'est pas le seul témoin de l'histoire de la cristallisation. La distribution de la taille des cristaux (CSD) fournit des informations essentielles sur les interactions entre nucléation et croissance des minéraux. La CSD des pseudomorphes de brucite retrograde formés après le périclase dans le sud du massif Adamello (Italie) est présentée dans la troisième partie. La combinaison entre les données textorales en trois dimensions et les données géochimiques a permis d'évaluer les cinétiques de réaction et de contraindre les mécanismes conduisant à la formation du périclase. Cette réaction est présentée comme étant la conséquence de l'infiltration d'une quantité limitée d'une phase fluide à haute température. La composition de cette phase fluide est en grand déséquilibre avec le reste de la roche, ce qui permet des cinétiques de réactions très rapides.

Characterization of functional role of α-ketoglutarate receptor in the kidney : role of the intrinsic circadian timing system in renin expressing cells

Le rein joue un rôle essential dans le maintien de l'homéostasie des fluides extracellulaires (FEC) et la pression artérielle. L'objectif de notre groupe est d'identifier de nouveaux mécanismes impliqués dans le contrôle de l'homéostasie des FEC et de la pression artérielle par le rein. Projet 1) Caractérisation du rôle fonctionnel du récepteur à l'a-cétogluatarate Oxgrl dans le rein Oxgrl est le récepteur spécifique de l'a-cétogluatarate, une moléule intermédiaire du cycle de l'acide citrique, filtrée par le rein et réabsorbée ou secrétée au niveau des tubules proximaux. Le rôle fonctionnel de ces deux récepteurs reste inconnu. Nos résultats montrent qu'Oxgrl est localisé au niveau des cellules intercalaires du tube collecteur (CCD). Des souris (Oxgrr/_) montrent une diminution du pH urinaire ,une augmentation de la concentration de l'acide urinaire titrable et une augmentation des niveaux d'a-cétoglutarate. Le traitement au Na-bicarbonate provoque une augmentation plus prononcée de l'alcalose métabolique chez les souris Oxgrl"7"' accompagnée d'une augmentation de la concentration de bicarbonate et une diminution du niveau de chlore plasmatique. En parallèle, des études de microperfusion ont montré que a-cétoglutarate stimule la réabsorption éléctroneutre de NaCl dans le CCD des souris de type sauvage mais pas des souris Oxgrl"7". En résumé, ces résultats montrent que l'a-cétoglutarate joue un rôle de molécule messagère du tubule proximal jusqu'au tube collecteur au niveau du rein et qu'Oxgrl pourrait être impliqué dans la régulation de l'échange Cl/bicarbonate et la réabsorption du NaCl dans les cellules intercalées. Projet 2) Rôle du système circadien dans les cellules productrices de rénine. Le système chronologique circadien est un mécanisme moléculaire ubiquitaire qui permet à l'organisme de coordonner ses fonctions principales en fonction du temps géophysique. Comme l'activité de la rénine plasmatique montre une rythmicité circadienne nette chez l'homme et la souris ; dans ce projet, nous avons abordé la question à savoir dans quelle mesure le système circadien est impliqué dans cette variabilité circadienne. Pour cela, le gène Bmall, élément principal de l'horloge moléculaire, a été perturbé dans les cellules granulaires productrices de rénine par le système Cre-LoxP. Nos résultats montrent que les souris Renld- Cre/Bmalllox/lox (cKO) présentent de faibles taux d'ARNm de Reni, altèrent la dynamique d'expression de la protéine rénine, mais il le niveau de concentration plasmatique de la rénine reste le même. Cependant, les souris cKO montrent une réduction significative de la concentration plasmatique de l'aldostérone. Nos analyses de l'urine récupérée dans des intervalles de temps de 24 et 1 heure montrent une augmentation du volume urinaire, une tendance à une hypercalciurie, ainsi qu'une altération de la dynamique d'excrétion urinaire de sodium chez les souris cKO. Plusieurs gènes impliqués dans la production/sécrétion de la rénine et dans le contrôle de la fonction rénale montrent une altération de l'expression circadienne d'ARNm. Par ailleurs, les souris cKO montrent une baisse significative de la pression artérielle. Nos résultats suggèrent que l'horloge intrinsèque des cellules productrices de la rénine joue un rôle important dans le control des FEC et l'homéostasie de la pression artérielle via régulation de la fonction rénale.

Partial melting of basic dykes in the contact aureoles of two shallow-level, gabbro-pyroxenite intrusions, Fuerteventura

THESIS ABSTRACT Low-pressure anatexis of basic dykes gave rise to unusual, zebra-like migmatites, in the contact metamorphic aureoles of two layered gabbro-pyroxenite intrusions, PXl and PX2, in the root zone of an ocean island, Fuerteventura Basal Complex (Canary Islands). This thesis focuses on the understanding of processes attributing to the partial melting and formation of these migmatites, characterised by a dense network of closely spaced, millimetre-wide leucocratic segregations with perfectly preserved igneous textures. The presence of fluids are required to decrease the solidus of basic igneous lithologies, to allow partial melting in such aloes-pressure (1-2 kb) environment. An oxygen isotope study was thus carried out on dykes inside and beyond the PX2 aureole, in order to decipher the nature and origin of such fluids. Low or negative δ18O values were obtained for whole rocks and mineral-separates, decreasing towards the contact, with the intrusion itself retaining fairly high values. This trend has been attributed to the advection of meteoric water during magma emplacement, with increasing fluid/rock ratios (higher dyke intensities towards the intrusion acting as fluid-pathways) and higher temperatures promoting increasing exchange during recrystallisation. A comparison of whole rock and mineral major- and trace- element data allowed the redistribution of elements .between different mineral phases and generations, during contact metamorphism and partial melting to be assessed. Certain trace-elements, e.g. Zr, Hf, Y, and REEs, were internally redistributed during contact metamorphic recrystallisation, causing- the enrichment of neocrystallised diopsides compared to relict phenocrysts. This has been assigned to the liberation of trace elements on the breakdown of primary minerals, kaersutite and sphene, on entering the thermal aureole. Major and trace element compositions of minerals in migmatite melanosomes and leucosomes are almost identical, pointing to a syn- or post- solidus reequilibration on cooling of the migmatite terrain. The mineralogical, textural and geochemical evolution of dykes in a contact metamorphic aureole, is recorded around an apophysis of the PX1 intrusion, where there is evidence of incipient partial melting. Hydrothermal mineral pseudomorphs in the outer parts of the aureole are progressively replaced by dry mineral assemblages, with increasingly recrystallised diopside and evidence of partial melting -the extent of which varies from one lithology to another. The appearance of more mafic lithologies towards the intrusion, with lower whole rock SiO2 and mobile element abundances, e.g. Rb, Cs, K, has been explained by the migration and accumulation of feldspathic material into leucosomes outside the samples. A micro-structural study of leucosomes and leucocratic pods, with the aid of high-resolution X-ray computed micro-tomography (HRXµCT), allowing the visualization and quantification of shapes and orientations, was carried out in order to better understand the processes of melt segregation in the PX1 aureole. Leucocratic pods, representing former amygdales, are considered as natural strain ellipsoids. Their short axes are oriented perpendicular to leucosome planes, which sub-parallel the intrusive contact. Leucosomes thus effectively represent foliation planes. This implies that the direction of maximum shortening, during migmatisation, was perpendicular to the orientation of leucosomes, contradicting earlier models that suggest leucosomes represent tension veins. RESUME DE LA THESE Un phénomène rare de fusion partielle de filons basiques à basse pression a été étudié dans les auréoles de contact de deux intrusions litées de gabbro-pyroxénite, PX1 et PX2, localisées dans le soubassement de l'île volcanique de Fuerteventura aux Canaries. Cette anatexie a engendré des migmatites finement zébrées d'aspect très inhabituel, dont les processus de formation ont été étudiés dans le présent travail. Ces roches sont caractérisées par un réseau dense de veinules leucocrates d'épaisseur millimétrique, dont les textures ignées sont parfaitement préservées. La fusion partielle de roches basiques à basse pression (1-2 kbar) requiert la présence d'eau afin d'abaisser le solidus du système à des températures géologiquement réalistes. Une étude comparative des isotopes de l'oxygène a ainsi été menée sur des filons respectivement affectés et non affectés par le métamorphisme de contact, afin de confirmer la présence de ces fluides, de déterminer l'importance de leur interaction avec les roches et leur origine. Des valeurs de δ180 basses ou négatives ont été mesurées sur roche totale et minéraux séparés, décroissantes en direction du contact, alors que l'intrusion elle-même a conservé des valeurs élevées. Ce gradient a été attribué à l'advection d'eau météorique durant la mise en place du magma, les températures les plus élevées favorisant d'autant plus la circulation des fluides et les échanges isotopiques durant la recristallisation des roches. Cette recristallisation engendré une redistribution chimique complète des éléments entre les différentes générations de minéraux résultant du métamorphisme de contact et de l'anatexie, mise en évidence par microanalyse. Certains éléments traces comme Zr, Hf, Y et les REE ont été concentrés dans le diopside néoformé consécutivement à la déstabilisation de minéraux primaires riches en ces éléments comme la kaersutite ou le sphène. Les compositions en éléments majeurs et traces des minéraux des mélanosomes et leucosomes des migmatites sont pratiquement identiques, indiquant une rééquilibration syn- à postsolidus lors du refroidissement de l'auréole de contact. La transformation progressive des filons basiques au niveau de leur minéralogie, textures et composition chimique a pu être observée en détail à l'approche du contact d'une apophyse de l'intrusion PX1. La paragenèse magmatique initiale n'est jamais préservée, les faciès les plus distants du contact étant constitués d'un assemblage pseudomorphique hydrothermal. Ce dernier est progressivement remplacé par des assemblages anhydres incluant du diopside néoformé, puis apparaissent les premiers signes de fusion partielle, dont l'importance varie fortement d'une lithologie à l'autre. L'apparition de faciès plus basiques en direction du contact, avec des teneurs réduites en SiO2 et en éléments incompatibles tels Rb, Cs, K, a été attribuée à l'échappement de leucosomes feldspathiques hors du système. Une étude microstructurale de la distribution spatiale du matériel leucocrate au sein des migmatites par microtomographie X de haute résolution (HRXµCT) a été menée pour mieux comprendre les processus de ségrégation des liquides dans l'auréole de PX1. De petites entités ovoïdes, représentant d'anciennes structures amygdalaires au sein des filons, ont été considérées comme des ellipsoïdes marqueurs de la déformation finie. Leur petit axe est orienté perpendiculairement aux plans définis par les leucosomes, eux-mêmes subparallèles au contact intrusif. Les leucosomes matérialisent donc des plans de clivage. Ainsi, la direction de raccourcissement maximum durant la fusion partielle était perpendiculaire à l'orientation des leucosomes, contrairement à ce qui a été dit dans de précédentes publications, qui suggéraient que les leucosomes représentaient des veines de tension. RESUME DE LA THESE (POUR LE GRAND PUBLIC) L'observation directe du soubassement d'une île volcanique est une occasion rare, accessible dans le «complexe de base » de l'île canarienne de Fuerteventura. Ce dernier a enregistré divers phénomènes magmatiques, métamorphiques et de fusion partielle induits par l'intrusion répétée de magmas alimentant des appareils volcaniques sus jacents, sous forme de petits plutons, essaims de filons et complexes annulaires de gabbros alcalins, pyroxénites, syénites et carbonatites. Dans ce contexte de flux de chaleur élevé, des filons basiques ont subi une fusion partielle au contact de deux intrusions de gabbro-pyroxénite, un phénomène extrêmement rare à une profondeur aussi réduite, estimée à quelque 3-6 km. Les produits de cette fusion partielle sont des liquides très riches en feldspath, concentrés en un réseau dense de veinules blanches (leucosomes) au sein du matériau résiduel sombre non fondu (mélanosome) pour former ce qu'on appelle des migmatites. Outre les aspects pétrologiques liés à la formation de ces migmatites, l'intérêt majeur du phénomène réside dans le fait qu'il puisse représenter la source des magmas évolués parfois observés sur les îles océaniques. A des pressions aussi faibles que dans le soubassement de Fuerteventura (1-2 Kbar), la présence de fluides abondants est nécessaire pour abaisser la température de début de fusion des roches (solidus) à des valeurs géologiquement réalistes. Des expériences ont montré que même en présence de plusieurs %-poids d'eau, une température de 1000°C était encore nécessaire pour obtenir une proportion de liquide équivalente à celle observée sur le terrain, soit 25%. Or les magmas alcalins des îles océaniques, bien qu'hydratés, n'en contiennent de loin pas autant, ce qui implique une source d'eau externe. Une étude isotopique de l'oxygène a été entreprise afin de tester cette hypothèse. Les valeurs obtenues en 5180 sont basses ou négatives et indiquent l'influence d'eau d'origine météorique. Cette eau de pluie se serait infiltrée le long des filons depuis la surface du volcan et les aurait complètement hydrothermalisés en profondeur (situation encore visible à l'extérieur de l'auréole de contact), leur permettant ainsi de stocker l'eau nécessaire à leur fusion partielle ultérieure. L'interaction entre eau de pluie et filons a été d'autant plus importante que ces derniers étaient proches du contact avec l'intrusion, ce qui suggère que la circulation de ces eaux et leur interaction avec les roches a été favorisée par la chaleur fournie par l'intrusion elle-même. Un autre aspect de ce travail s'est focalisé sur la redistribution des éléments traces au sein des minéraux des filons basiques durant le métamorphisme de contact et la fusion partielle. Ainsi, le pyroxène de seconde génération est-il sensiblement enrichi en traces telles Zr, Hf, Y et les terres rares, par rapport au pyroxène magmatique originel, en relation avec la déstabilisation de minéraux primaires riches en ces éléments tels le sphène et la kaersutite. Cependant, les compositions en éléments majeurs et traces des minéraux recristallisés des migmatites sont pratiquement les mêmes dans les leucosomes et les mélanosomes, suggérant une rééquilibration chimique complète durant le refroidissement de ces lithologies. Si certaines migmatites se sont comportées en système chimiquement fermé (hormis l'eau météorique), d'autres filons ont manifestement perdu une partie de leurs leucosomes, ainsi qu'en témoigne leur composition progressivement appauvrie en silice et autres éléments incompatibles mobiles, tels K et Rb à l'approche du contact de l'intrusion. Parallèlement à cette évolution chimique, les paragenèses hydrothermales distantes du contact sont progressivement remplacées par des paragenèses anhydres, puis par l'apparition des premiers leucosomes, tandis que les textures magmatiques initiales sont complètement effacées au profit d'une combinaison de textures magmatique dans les leucosomes et en mosaïque dans les mélanosomes. Enfin, la distribution spatiale des liquides de fusion partielle a été étudiée par microtomographie X de haute résolution, dans des filons contenant des entités ovoïdes leucocrates, sans doute d'anciennes amygdales à zéolites. Ces dernières ont été considérées comme des ellipsoïdes de la déformation finie. L'orientation de leur petit axe, perpendiculaire au plan défini par les veinules de leucosomes, indique que ces derniers représentent des plans de clivage perpendiculaires à la direction de raccourcissement maximum. Ainsi, la ségrégation des liquides de fusion partielle se serait faite dans les plans de compression et non dans des plans de dilatation, contrairement à ce que laisserait penser le sens commun.

Sur une conception nouvelle des forces intérieures dans un fluide en mouvement / par M. Stanislas Zaremba,...

Collection : Mémorial des sciences mathématiques ; fasc. 82

Time-lapse and probabilistic inversion strategies for plane-wave electromagnetic data.

L'utilisation efficace des systèmes géothermaux, la séquestration du CO2 pour limiter le changement climatique et la prévention de l'intrusion d'eau salée dans les aquifères costaux ne sont que quelques exemples qui démontrent notre besoin en technologies nouvelles pour suivre l'évolution des processus souterrains à partir de la surface. Un défi majeur est d'assurer la caractérisation et l'optimisation des performances de ces technologies à différentes échelles spatiales et temporelles. Les méthodes électromagnétiques (EM) d'ondes planes sont sensibles à la conductivité électrique du sous-sol et, par conséquent, à la conductivité électrique des fluides saturant la roche, à la présence de fractures connectées, à la température et aux matériaux géologiques. Ces méthodes sont régies par des équations valides sur de larges gammes de fréquences, permettant détudier de manières analogues des processus allant de quelques mètres sous la surface jusqu'à plusieurs kilomètres de profondeur. Néanmoins, ces méthodes sont soumises à une perte de résolution avec la profondeur à cause des propriétés diffusives du champ électromagnétique. Pour cette raison, l'estimation des modèles du sous-sol par ces méthodes doit prendre en compte des informations a priori afin de contraindre les modèles autant que possible et de permettre la quantification des incertitudes de ces modèles de façon appropriée. Dans la présente thèse, je développe des approches permettant la caractérisation statique et dynamique du sous-sol à l'aide d'ondes EM planes. Dans une première partie, je présente une approche déterministe permettant de réaliser des inversions répétées dans le temps (time-lapse) de données d'ondes EM planes en deux dimensions. Cette stratégie est basée sur l'incorporation dans l'algorithme d'informations a priori en fonction des changements du modèle de conductivité électrique attendus. Ceci est réalisé en intégrant une régularisation stochastique et des contraintes flexibles par rapport à la gamme des changements attendus en utilisant les multiplicateurs de Lagrange. J'utilise des normes différentes de la norme l2 pour contraindre la structure du modèle et obtenir des transitions abruptes entre les régions du model qui subissent des changements dans le temps et celles qui n'en subissent pas. Aussi, j'incorpore une stratégie afin d'éliminer les erreurs systématiques de données time-lapse. Ce travail a mis en évidence l'amélioration de la caractérisation des changements temporels par rapport aux approches classiques qui réalisent des inversions indépendantes à chaque pas de temps et comparent les modèles. Dans la seconde partie de cette thèse, j'adopte un formalisme bayésien et je teste la possibilité de quantifier les incertitudes sur les paramètres du modèle dans l'inversion d'ondes EM planes. Pour ce faire, je présente une stratégie d'inversion probabiliste basée sur des pixels à deux dimensions pour des inversions de données d'ondes EM planes et de tomographies de résistivité électrique (ERT) séparées et jointes. Je compare les incertitudes des paramètres du modèle en considérant différents types d'information a priori sur la structure du modèle et différentes fonctions de vraisemblance pour décrire les erreurs sur les données. Les résultats indiquent que la régularisation du modèle est nécessaire lorsqu'on a à faire à un large nombre de paramètres car cela permet d'accélérer la convergence des chaînes et d'obtenir des modèles plus réalistes. Cependent, ces contraintes mènent à des incertitudes d'estimations plus faibles, ce qui implique des distributions a posteriori qui ne contiennent pas le vrai modèledans les régions ou` la méthode présente une sensibilité limitée. Cette situation peut être améliorée en combinant des méthodes d'ondes EM planes avec d'autres méthodes complémentaires telles que l'ERT. De plus, je montre que le poids de régularisation des paramètres et l'écart-type des erreurs sur les données peuvent être retrouvés par une inversion probabiliste. Finalement, j'évalue la possibilité de caractériser une distribution tridimensionnelle d'un panache de traceur salin injecté dans le sous-sol en réalisant une inversion probabiliste time-lapse tridimensionnelle d'ondes EM planes. Etant donné que les inversions probabilistes sont très coûteuses en temps de calcul lorsque l'espace des paramètres présente une grande dimension, je propose une stratégie de réduction du modèle ou` les coefficients de décomposition des moments de Legendre du panache de traceur injecté ainsi que sa position sont estimés. Pour ce faire, un modèle de résistivité de base est nécessaire. Il peut être obtenu avant l'expérience time-lapse. Un test synthétique montre que la méthodologie marche bien quand le modèle de résistivité de base est caractérisé correctement. Cette méthodologie est aussi appliquée à un test de trac¸age par injection d'une solution saline et d'acides réalisé dans un système géothermal en Australie, puis comparée à une inversion time-lapse tridimensionnelle réalisée selon une approche déterministe. L'inversion probabiliste permet de mieux contraindre le panache du traceur salin gr^ace à la grande quantité d'informations a priori incluse dans l'algorithme. Néanmoins, les changements de conductivités nécessaires pour expliquer les changements observés dans les données sont plus grands que ce qu'expliquent notre connaissance actuelle des phénomenès physiques. Ce problème peut être lié à la qualité limitée du modèle de résistivité de base utilisé, indiquant ainsi que des efforts plus grands devront être fournis dans le futur pour obtenir des modèles de base de bonne qualité avant de réaliser des expériences dynamiques. Les études décrites dans cette thèse montrent que les méthodes d'ondes EM planes sont très utiles pour caractériser et suivre les variations temporelles du sous-sol sur de larges échelles. Les présentes approches améliorent l'évaluation des modèles obtenus, autant en termes d'incorporation d'informations a priori, qu'en termes de quantification d'incertitudes a posteriori. De plus, les stratégies développées peuvent être appliquées à d'autres méthodes géophysiques, et offrent une grande flexibilité pour l'incorporation d'informations additionnelles lorsqu'elles sont disponibles. -- The efficient use of geothermal systems, the sequestration of CO2 to mitigate climate change, and the prevention of seawater intrusion in coastal aquifers are only some examples that demonstrate the need for novel technologies to monitor subsurface processes from the surface. A main challenge is to assure optimal performance of such technologies at different temporal and spatial scales. Plane-wave electromagnetic (EM) methods are sensitive to subsurface electrical conductivity and consequently to fluid conductivity, fracture connectivity, temperature, and rock mineralogy. These methods have governing equations that are the same over a large range of frequencies, thus allowing to study in an analogous manner processes on scales ranging from few meters close to the surface down to several hundreds of kilometers depth. Unfortunately, they suffer from a significant resolution loss with depth due to the diffusive nature of the electromagnetic fields. Therefore, estimations of subsurface models that use these methods should incorporate a priori information to better constrain the models, and provide appropriate measures of model uncertainty. During my thesis, I have developed approaches to improve the static and dynamic characterization of the subsurface with plane-wave EM methods. In the first part of this thesis, I present a two-dimensional deterministic approach to perform time-lapse inversion of plane-wave EM data. The strategy is based on the incorporation of prior information into the inversion algorithm regarding the expected temporal changes in electrical conductivity. This is done by incorporating a flexible stochastic regularization and constraints regarding the expected ranges of the changes by using Lagrange multipliers. I use non-l2 norms to penalize the model update in order to obtain sharp transitions between regions that experience temporal changes and regions that do not. I also incorporate a time-lapse differencing strategy to remove systematic errors in the time-lapse inversion. This work presents improvements in the characterization of temporal changes with respect to the classical approach of performing separate inversions and computing differences between the models. In the second part of this thesis, I adopt a Bayesian framework and use Markov chain Monte Carlo (MCMC) simulations to quantify model parameter uncertainty in plane-wave EM inversion. For this purpose, I present a two-dimensional pixel-based probabilistic inversion strategy for separate and joint inversions of plane-wave EM and electrical resistivity tomography (ERT) data. I compare the uncertainties of the model parameters when considering different types of prior information on the model structure and different likelihood functions to describe the data errors. The results indicate that model regularization is necessary when dealing with a large number of model parameters because it helps to accelerate the convergence of the chains and leads to more realistic models. These constraints also lead to smaller uncertainty estimates, which imply posterior distributions that do not include the true underlying model in regions where the method has limited sensitivity. This situation can be improved by combining planewave EM methods with complimentary geophysical methods such as ERT. In addition, I show that an appropriate regularization weight and the standard deviation of the data errors can be retrieved by the MCMC inversion. Finally, I evaluate the possibility of characterizing the three-dimensional distribution of an injected water plume by performing three-dimensional time-lapse MCMC inversion of planewave EM data. Since MCMC inversion involves a significant computational burden in high parameter dimensions, I propose a model reduction strategy where the coefficients of a Legendre moment decomposition of the injected water plume and its location are estimated. For this purpose, a base resistivity model is needed which is obtained prior to the time-lapse experiment. A synthetic test shows that the methodology works well when the base resistivity model is correctly characterized. The methodology is also applied to an injection experiment performed in a geothermal system in Australia, and compared to a three-dimensional time-lapse inversion performed within a deterministic framework. The MCMC inversion better constrains the water plumes due to the larger amount of prior information that is included in the algorithm. The conductivity changes needed to explain the time-lapse data are much larger than what is physically possible based on present day understandings. This issue may be related to the base resistivity model used, therefore indicating that more efforts should be given to obtain high-quality base models prior to dynamic experiments. The studies described herein give clear evidence that plane-wave EM methods are useful to characterize and monitor the subsurface at a wide range of scales. The presented approaches contribute to an improved appraisal of the obtained models, both in terms of the incorporation of prior information in the algorithms and the posterior uncertainty quantification. In addition, the developed strategies can be applied to other geophysical methods, and offer great flexibility to incorporate additional information when available.

Estudi microtermomètric dels filons de quars de les Esquerdes de Rojà, Massís del Canigó, Pirineu Oriental

Els materials pre-hercinians i el granit de Costabona estan travessats per un conjunt de filons de quars de fins a 15 km de longitud, coneguts com a 'Esquerdes de Roja'. Aquests filons estan localment afectats per un sistema de bandes milonítiques i per un altre posterior de fractures. A partir de l'anilisi microtermometrica de les inclusions fluides primàries en quars, ha estat obtinguda una temperatura de formació de l'ordre de + 280 C per la major part del quars, i una temperatura de + 250 C per als monocristalls que omplen esquerdes i cavitats. Els filons s'haurien format a partir de solucions mineralitzants que contenien ~ a + C, a++, K+ i, en menor proporció, ~ g + + .fins una salinitat total de l'ordre del 20/@ en pes de la Na CI.

Fluid-rock-strain feedbacks in extensional shear zones

La présence de fluide météorique synchrone à l'activité du détachement (Farmin, 2003 ; Mulch et al., 2007 ; Gébelin et al., 2011), implique que les zones de cisaillement sont des systèmes ouverts avec des cellules de convections à l'échelle crustale et un intense gradient géothermique au sein du détachement (Morrison et Anderson, 1998, Gottardi et al., 2011). De plus, les réactions métamorphiques liées à des infiltrations fluides dans les zones de cisaillement extensionnel peuvent influencer les paramètres rhéologiques du système (White and Knipe, 1978), et impliquer la localisation de la déformation dans la croûte. Dans ce manuscrit, deux zones de cisaillement infiltrées par des fluides météoriques sont étudiées, l'une étant largement quartzitique, et l'autre de nature granitique ; les relations entre déformation, fluides, et roches s'appuient sur des approches structurales, microstructurales, chimiques et isotopiques. L'étude du détachement du Columbia river (WA, USA) met en évidence que la déformation mylonitique se développe en un million d'années. La phase de cisaillement principal s'effectue à 365± 30°C d'après les compositions isotopiques en oxygène du quartz et de la muscovite. Ces minéraux atteignent l'équilibre isotopique lors de leur recristallisation dynamique contemporaine à la déformation. La zone de cisaillement enregistre une baisse de température, remplaçant le mécanisme de glissement par dislocation par celui de dissolution- précipitation dans les derniers stades de l'activité du détachement. La dynamique de circulation fluide bascule d'une circulation pervasive à chenalisée, ce qui engendre localement la rupture des équilibres d'échange isotopiques. La zone de cisaillement de Bitterroot (MT, USA) présente une zone mylonitique de 600m d'épaisseur, progressant des protomylonites aux ultramylonites. L'intensité de la localisation de la déformation se reflète directement sur l'hydratation des feldspaths, réaction métamorphique majeure dite de « rock softening ». Une étude sur roche totale indique des transferts de masse latéraux au sein des mylonites, et d'importantes pertes de volume dans les ultramylonites. La composition isotopique en hydrogène des phyllosilicates met en évidence la présence (1) d'une source magmatique/métamorphique originelle, caractérisée par les granodiorites ayant conservé leur foliation magmatique, jusqu'aux protomylonites, et (2) une source météorique qui tamponne les valeurs des phyllosilicates des fabriques mylonitiques jusqu'aux veines de quartz non-déformées. Les compositions isotopiques en oxygène des minéraux illustrent le tamponnement de la composition du fluide météorique par l'encaissant. Ce phénomène cesse lors du processus de chloritisation de la biotite, puisque les valeurs des chlorites sont extrêmement négatives (-10 per mil). La thermométrie isotopique indique une température d'équilibre isotopique de la granodiorite entre 600-500°C, entre 500-300°C dans les mylonites, et entre 300 et 200°C dans les fabriques cassantes (cataclasites et veines de quartz). Basé sur les résultats issus de ce travail, nous proposons un modèle général d'interactions fluide-roches-déformation dans les zones de détachements infiltrées par des fluides météoriques. Les zones de détachements évoluent rapidement (en quelques millions d'années) au travers de la transition fragile-ductile ; celle-ci étant partiellement contrôlée par l'effet thermique des circulations de fluide météoriques. Les systèmes de détachements sont des lieux où la déformation et les circulations fluides sont couplées ; évoluant rapidement vers une localisation de la déformation, et de ce fait, une exhumation efficace. - The presence of meteoric fluids synchronous with the activity of extensional detachment zones (Famin, 2004; Mulch et al., 2007; Gébelin et al., 2011) implies that extensional systems involve fluid convection at a crustal scale, which results in high geothermal gradients within active detachment zones (Morrison and Anderson, 1998, Gottardi et al., 2011). In addition, the metamorphic reactions related to fluid infiltration in extensional shear zones can influence the rheology of the system (White and Knipe, 1978) and ultimately how strain localizes in the crust. In this thesis, two shear zones that were permeated by meteoric fluids are studied, one quartzite-dominated, and the other of granitic composition; the relations between strain, fluid, and evolving rock composition are addressed using structural, microstructural, and chemical/isotopic measurements. The study of the Columbia River detachment that bounds the Kettle core complex (Washington, USA) demonstrates that the mylonitic fabrics in the 100 m thick quartzite- dominated detachment footwall developed within one million years. The main shearing stage occurred at 365 ± 30°C when oxygen isotopes of quartz and muscovite equilibrated owing to coeval deformation and dynamic recrystallization of these minerals. The detachment shear zone records a decrease in temperature, and dislocation creep during detachment shearing gave way to dissolution-precipitation and fracturing in the later stages of detachment activity. Fluid flow switched from pervasive to channelized, leading to isotopic disequilibrium between different minerals. The Bitterroot shear zone detachment (Montana, USA) developed a 600 m thick mylonite zone, with well-developed transitions from protomylonite to ultramylonite. The localization of deformation relates directly to the intensity of feldspar hydration, a major rock- softening metamorphic reaction. Bulk-rock analyses of the mylonitic series indicate lateral mass transfer in the mylonite (no volume change), and significant volume loss in ultramylonite. The hydrogen isotope composition of phyllosilicates shows (1) the presence of an initial magmatic/metamorphic source characterized by the granodiorite in which a magmatic, and gneissic (protomylonite) foliation developed, and (2) a meteoric source that buffers the values of phyllosilicates in mylonite, ultramylonite, cataclasite, and deformed and undeformed quartz veins. The mineral oxygen isotope compositions were buffered by the host-rock compositions until chloritization of biotite started; the chlorite oxygen isotope values are negative (-10 per mil). Isotope thermometry indicates a temperature of isotopic equilibrium of the granodiorite between 600-500°C, between 500-300°C in the mylonite, and between 300 and 200°C for brittle fabrics (cataclasite and quartz veins). Results from this work suggest a general model for fluid-rock-strain feedbacks in detachment systems that are permeated by meteoric fluids. Phyllosilicates have preserved in their hydrogen isotope values evidence for the interaction between rock and meteoric fluids during mylonite development. Fluid flow generates mass transfer along the tectonic anisotropy, and mylonites do not undergo significant volume change, except locally in ultramylonite zones. Hydration of detachment shear zones attends mechanical grain size reduction and enhances strain softening and localization. Self-exhuming detachment shear zones evolve rapidly (a few million years) through the transition from ductile to brittle, which is partly controlled by the thermal effect of circulating surface fluids. Detachment systems are zones in the crust where strain and fluid flow are coupled; these systems. evolve rapidly toward strain localization and therefore efficient exhumation.

Hydrogen isotope fractionations between amphiboles, micas, and fluids in alkaline igneous intrusions

RÉSUMÉ DE LA THÈSE Les teneurs des amphiboles en éléments majeurs et en isotopes stables ont été analysées dans plusieurs complexes ignés alcalins et hyperalcalins, dans le but de déterminer l'importance des variations de composition des minéraux pour le fractionnement isotopique de l'hydrogène dans un système naturel minéral-magma-fluide. Cette étude se concentre principalement sur les syénites néphéliniques de complexes intrusifs alcalins bien connus mais à chimie variable, dont les amphiboles, ainsi que d'autres silicates hydratés tels que micas et eudialytes, lorsque cela était possible, ont été séparés. L'intérêt principal s'est porté sur le complexe alcalin d'Ilímaussaq de la Province du Gardar, au Sud du Groenland. Dans une optique de comparaison, nous avons collecté et analysé d'autres échantillons provenant du complexe de Tugtutôq (Sud Groenland), des complexes de Khibina et Lovozero (Péninsule de Kola, Russie), du Mont St-Hilaire et du Mont Royal (Canada) et de 6 autres du nord-ouest de la Namibie (Cape Cross, Okenyenya, Messum, Etaneno, Kalkfeld,et Okorusu). Les compositions isotopiques de l'hydrogène des amphiboles des ces différentes zones présentent de grandes variations (-227 à -700/00), ce qui est atypique pour des magmas d'origine mantellique. Les valeurs comprises entre -80 et -400/00 indiquent une provenance du manteau. Ces larges variations de compositions ainsi que l'extrême appauvrissement en isotope lourd de l'hydrogène (D), en comparaison avec d'autres roches ignées, semblent être propres.aux roches alcalines et hyperalcalines de ce type, ce qui indiquerait un processus commun. Les différents complexes alcalins choisis présentent un large intervalle de composition chimique des amphiboles. La caractérisation des amphiboles par microscopie électronique et par spectroscopie Mössbauer contribuent à observer le contrôle du Fe sur le fractionnement des isotopes de l'hydrogène. En effet, cela a mis en évidence un contrôle du Fe sur le fractionnement et même, dans le cas du complexe hyperalcalin d'Ilímaussaq, une relation entre le rapport Fei3+/FeT et les variations du rapport D/H. Les complexes étudiés diffèrent de par leur index agpaïtique (Na+K/Al) et également de par leur contenu en fer. Les plus hautes valeurs en Fe (27-35 wt%) et en éléments alcalins dans les amphiboles, ainsi que les teneurs de D/H les plus basses et leur grande variation, sont celles du complexe d'Ilímaussaq. Les amphiboles de la Péninsule de Kola et du Canada sont similaires, mais toutefois moins appauvries en D. En ce qui concerne les amphiboles des complexes du NO de la Namibie, elles présentent des compositions isotopiques de l'hydrogène magmatiques normales (-73 à -100 0/00), contiennent moins de Fe (15-17 wt%) et sont fortement enrichies en Ca et moins en Na. Dans ce cas, l'alcalinité est moins importante en comparaison des autres complexes étudiés. En dehors des teneurs en éléments alcalins des amphiboles, l'alcalinité des fluides s'avère également un facteur important, ce qui est cohérent avec certaines suggestions à partir de systèmes expérimentaux. Afin de mieux contraindre ce facteur, des expériences d'échanges hydrothermaux entre les amphiboles et les fluides de salinité différente ont été effectuées en simulant des conditions naturelles. L'approximation d'amphiboles naturelles de complexes ignés alcalins, couplée aux expériences d'échange, aide à préciser les facteurs contrôlant le fractionnement des isotopes de l'hydrogène dans les roches alcalines. Les valeurs extrêmement basses de 3D des amphiboles de ces complexes alcalins peuvent être dues à une combinaison de différents facteurs, telles qu'une haute alcalinité, une haute teneur en Fe et une faible profondeur d'intrusion. Les grandes variations ainsi que les faibles valeurs de SD des amphiboles étudiées peuvent résulter d'un processus magmatique interne et il est peu probable que de l'eau météorique soit impliquée et/ou que le dégazage magmatique ait joué un rôle. THESIS ABSTRACT Major element and stable isotope compositions of amphiboles were analyzed from a number of alkaline and peralkaline igneous complexes in order to determine the importance of compositional variations in minerals to hydrogen isotope fractionations in natural mineral-melt-fluid systems. The thesis mainly focuses on nepheline syenites of well-studied, but chemically variable alkaline intrusive rocks, from which amphiboles and, if possible, other hydrous silicates such as micas and eudialytes were separated. The system of primary interest was the alkaline Ilímaussaq Complex of the Gardar Province of South Greenland. For the purpose of comparison additional samples were collected and examined from the Tugtutôq Complex (South Greenland), the Khibina and Lovozero Complexes (Kola Peninsular, Russia), Mount St-Hilaire and Mount Royal (Canada) and six further complexes from NW Namibia (Cape Cross, Okenyenya, Messum, Etaneno, Kalkfeld, and Okorusu). The hydrogen isotope compositions of amphiboles from the localities studied differ greatly, which is atypical for amphiboles from mantle, range between - 227 and - 700/00 (latter compatible with a simple mantle origin). As this wide range in compositions and the extreme depletion in the heavy hydrogen isotope (D) content relative to other igneous rocks appear to be unique to alkaline to peralkaline rocks of this type, a common process is indicated. The different alkaline complexes chosen cover a wide range of amphibole chemical compositions. Detailed chemical characterization of amphiboles by electron microprobe and Mössbauer spectroscopy analyses helped to constrain the control of Fe on the H-isotope fractionations. Complete characterization of the chemical compositions of the amphiboles support Fe-control on fractionations and at least for the peralkaline Ilímaussaq complex a relationship between Fe3+/FeT ratios and variations in D/H. The studied complexes differ in their agpaitic index (Na+K/Al) and also in their Fe-content. The most iron (27-35 wt. %) and alkaline element rich amphiboles, with the lowermost D/H ratio, as well with very wide range, are the ones from Ilímaussaq complex. Similar, but less D depleted amphiboles are from the Kola Peninsula and the Canadian localities. The complexes described from NW Namibia have amphiboles with normal magmatic hydrogen isotope composition (-730/00 to -1000/00), and have less Fe-content (15-17 wt. %), and are more Ca-and less Na-rich. In this case alkalinity is not that important in comparison to the other studied complexes. Beside the alkaline element contents in the amphiboles, the alkalinity of the fluids has been found to be an important factor, in conjunction with earlier suggestions from experimental systems. To further constrain this factor, hydrothermal exchange experiments between amphiboles and fluids of different salinity simulating natural conditions were performed. The approach of examining natural amphiboles from alkaline igneous complexes in parallel to performing exchange experiments - helped to further constrain the factors controlling the H-isotope fractionations in alkaline rocks. The observed changes between the hydrogen and oxygen isotope compositions of amphiboles and fluids before and after the experiments suggest that another phase was produced during the experiments, which influenced the final hydrogen isotope composition of the system. This presumably hydrous phase has also influenced the Fe3 +/Fe2+ ratio of the amphiboles, which became more oxidized. The extremely low SD values of amphiboles in these alkaline complexes may be due to a combination of different factors such as high alkalinity, high Fe-content, and shallow intrusion depths. This wide range and the low SD values of the amphiboles studied might be a result of internal, magmatic processes and it is unlikely that meteoric water was involved and/or magmatic degassing played an important role. RÉSUMÉ DE LA THÈSE (pour le grand public) Fractionnement isotopique de l'hydrogène entre amphiboles, micas et fluides dans des intrusions alcalines Zsófia Wáczek Directeur de thèse, Prof. Torsten W. Vennemann Institut de Minéralogie et Géochimie, Université de Lausanne Les roches alcalines et celles qui leurs sont associées sont des sources importantes de nombreux minéraux et minerais, tels l'apatite, le niobium, le diamant et autres pierres précieuses. Cette étude se concentre sur des complexes alcalins localisés dans le sud du Groenland, au Canada, dans la péninsule de Kola en Russie et au nord-ouest de la Namibie. Ces complexes sont composés de roches ayant cristallisé à partir de magmas et de fluides très enrichis en alcalins. Cet enrichissement permet la précipitation de minéraux inhabituels riches en potassium et/ou sodium, telles les amphiboles sodiques, également enrichies en fer. Les amphiboles étudiées ont des compositions calciques, sodi-calciques et sodiques, qui reflètent leurs différents environnements de formation. Des études précédentes ont révélé une large gamme de rapports isotopiques de l'hydrogène dans les amphiboles de roches hyperalcalines, dont certains extrêmement bas. Cette variation importante est très intrigante, sachant que des valeurs entre -40 et -800/00 correspondent à des silicates ignés hydratés et non altérés, alors que des valeurs descendant jusqu'a -1500/00 nécessiteraient une altération par de l'eau météorique et/ou une contamination par les roches environnantes ou des sédiments riches en matière organique. Dans lé cas précis du complexe d'Ilímaussaq (sud du Groenland), aucune de ces explications n'a pu être démontrée et des valeurs encore plus faibles ont été trouvées. Le complexe d'Ilímaussaq présente des valeurs de rapport isotopique de l'hydrogène entre -227 et -500/00 dans les amphiboles. Une origine mantellique permet d'expliquer les valeurs élevées, mais d'autres processus doivent entrer en jeu pour engendrer les valeurs les plus négatives. C'est à l'identification de ces processus que nous nous sommes attachés dans ce travail. Les grandes variations observées dans les teneurs en fer et dans le rapport Fe3+/FeT des roches et des minéraux de ces complexes sont corrélées avec d'autres paramètres chimiques, tels que la composition isotopique de l'hydrogène dans les amphiboles. Nous avons dès lors abordé les questions suivantes: quelle est la relation entre la teneur en fer des amphiboles et leur composition isotopique? Que nous apprennent les changements de la teneur en fer et les changements dans le rapport Fe3+/FeT sur les processus pétrologiques dans ces roches? Pour répondre à ces questions, nous avons analysé les compositions isotopiques de l'oxygène et de l'hydrogène dans les amphiboles et d'autres silicates hydratés. La composition chimique et le rapport Fe3+/FeT des amphiboles ont également été déterminés. Des expériences hydrothermales simulant des conditions naturelles ont été entreprises afin de mieux comprendre les processus de fractionnement isotopiques dans ces systèmes très alcalins. Nos conclusions sont les suivantes: (1) Les valeurs extrêmement faibles ainsi que les larges variations des rapports isotopiques de l'hydrogène des amphiboles de ces complexes alcalins sont dues à une combinaison de facteurs tels que la forte alcalinité, la haute teneur en fer et la profondeur très faible de l'intrusion. (2) Ces valeurs sont probablement le résultat de processus magmatiques internes. (3) Il est peu probable que les eaux météoriques et/ou le dégazage magmatique aient joué un rôle lors de la formation de ces amphiboles. (4) Certaines corrélations, en accord avec les études précédentes, ont pu être trouvées au niveau des concentrations en fer. (5) Dans le cas du complexe d'Ilímaussaq exclusivement, une relation a été trouvée entre le rapport Fe3+/FeT et la composition isotopique de l'hydrogène des amphiboles.

Elemental and isotopic zoning in natural alpine quartz

THESIS ABSTRACT : Stable isotope geochemistry is used to help resolve a large number of geological questions. In order to do this, it is essential to understand the different mechanisms that govern isotopic fractionation processes between different phases and to identify the conditions required to reach equilibrium fractionation. However, at low temperatures, these processes are poorly constrained and many factors can induce differential partitioning of the isotopes between sectors of a mineral species and the fluid during mineral growth. This can result in so-called 'sector zoning' of a mineral species. The aim of this thesis is to evaluate the occurrence of sector zoning of the oxygen isotopes and trace elements in natural α-quartz crystals and to identify the reasons for such zoning. The implications for the fluid-mineral interactions are studied in the context of the Alpine metamorphism. The approach chosen has focused on examining the crystal structure, cathodoluminescence appearance (CL), and on relating elemental (e.g. Li, Na, Al, P, K, Ca, Ge, Ti, Fe) to stable oxygen isotope compositions between and along different growth sectors. Low temperature quartz samples were selected from Alpine veins in different localities, where growth conditions have already been well constrained. The mineralogy as well as the isotopic compositions of the host rocks were also investigated, in order to interpret the variations obtained between the different growth stages in the framework of fluid-rock interaction during Alpine metamorphism. Depending on the growth conditions, most of the studied quartz is strongly zoned in CL, and it reveals corresponding zonations in the trace element content (e.g. growth zoning). Aluminium, substituting for Si in the lattice, was found in concentrations up to 1000's ppma, and its distribution is strongly related to Li and H and to a lesser extent, to Ge. Elemental sector zoning is evident from the distribution of these three elements since they exhibit differences in their respective concentrations between faces for distinct growth zones, with prismatic faces having the lowest Al contents. Quartz from veins in magmatic rocks, for example, tend to have lower Al concentrations and similar concentrations of Li and Ti suggesting also a contribution of these elements from the host rock. The relationship between Al and Li is still correlated. Only Alpine crystals grown at higher temperatures (~400°C) without any CL zoning feature are free of these impurities and do not show such zoning characteristics. Differences in the δ18O values were measured between different faces principally in the AIenriched growth zones or stages. These results were confirmed by the means of two different methods (in situ/non in situ). However, it was determined that the Al concentrations do not affect significantly oxygen isotope fractionations at 300°C. The results altogether suggest that the presence of sector zoning in quartz crystals is real, but not universal, and henceforth should be taken into consideration for any use of these systems. The occurrence of disequilibrium partitioning has been enhanced and is possibly related to kinetic processes as well as structural effects that do not affect similarly trace element incorporation and isotopic fractionation. In situ measurements also revealed fine scale δ18O zonations along growth paths that are useful to constrain fluid-rock interactions during Alpine metamorphism. Variations in the δ18O values present along growth vectors indicate changes in the fluid composition and origin. Association with oxygen isotope composition of the host rock allows for the deduction of interactions between rocks, veins and consequently fluids, as well as fluid regimes. RESUME DE LA THESE : A basses températures, (i.e. <400°C) les différents mécanismes qui régissent le fractionnement isotopique ainsi que les conditions nécessaires pour établir un état d'équilibre sont peu connus et nombre de paramètres peuvent entraîner un partitionnement chimique différentiel entre différents secteurs d'un minéral et le fluide en contact. Ainsi, ce travail de thèse a pour but d'évaluer la possible présence de zonages sectoriels en isotopes de l'oxygène mais aussi en éléments traces dans des cristaux naturels de quartz-α de basses températures, ainsi que les raisons d'un tel phénomène et enfin ses implications sur les interactions fluide-roche, principalement dans le cadre du métamorphisme Alpin. La structure et l'apparence en cathodoluminescence (CL) des échantillons ont été caractérisées avant de retracer en détail les compositions en élément traces (Li, Na, Al, P, K, Ca, Ge, Ti, Fe) et en isotopes de l'oxygène, le long et entre différents secteurs. Les échantillons de quartz sélectionnés proviennent majoritairement de veines Alpine de différentes localités, où les conditions de croissance ont été déjà bien caractérisées. Les compositions minéralogiques et isotopiques de la roche encaissante ont aussi été examinées, pour contraindre les variations obtenues dans un contexte Alpin. Selon leurs conditions de croissance, la plupart des cristaux étudiés sont fortement zonés, ce qui est souligné par un zonage des concentrations en éléments traces (e.g. zonage de croissance). L'Aluminium, qui peut se substituer à la Silice dans le réseau cristallin, a été retrouvé jusqu'en très grandes concentrations dans certaines zones (plusieurs milliers de ppma). De plus, la distribution en Al est fortement liée à celles de Li et H, ainsi que dans une moindre mesure à Ge. La présence de zonage sectoriel est évidente au niveau de ces éléments qui montrent de larges différences de concentrations entre différentes faces pour une même zone de croissance, avec les concentrations les plus basses retrouvées dans les faces prismatiques. Les quartz de veines situées dans des roches magmatiques par exemple possèdent des concentrations en Li et Ti de même ordre de grandeur, confirmant le rôle de la composition de la roche encaissante. La relation Li/Al est toujours fortement présente, mais ce rapport est fonction de la face mesurée. Seuls les cristaux Alpins de plus hautes températures (400°C) ne possédant pas de zones en CL ne présentent aucune de ces caractéristiques. Des différences dans les valeurs de δ18O de zones identiques enrichies en Al ont clairement été mesurées entre les différentes faces r, z, et m, mais aussi au sein d'une même seule zone, indiquant que le fractionnement a probablement eu lieu en déséquilibre. Il a été déterminé que la présence d'Al dans ces teneurs n'avait qu'un faible effet sur le fractionnement isotopique de l'oxygène. L'utilisation de deux méthodes différentes a permis d'obtenir des résultats in situ et non in situ concordants. La comparaison des résultats obtenus permet de démontrer que le zonage sectoriel est bien présent dans certains cristaux de quartz, et dépend des conditions de formation. La présence d'un partitionnement différentiel des éléments traces peut être due à des effets cinétiques aussi bien que structuraux, alors que le zonage sectoriel des isotopes de l'oxygène aurait d'autres origines. Il est alors évident que la possibilité de zonage sectoriel doit être désormais pris en considération avant toute interprétations de données isotopiques de cristaux zonés. Les mesures in situ ont de plus permis de distinguer de fines variations des valeurs δ18O au cours de la croissance, qui peuvent aider à retracer la circulations des fluides dans les Alpes durant cette période. En association avec les compositions des roches encaissantes, ii est possible de déduire les interactions entre roches, veines, et par conséquent fluides, au cours de différentes étapes. RESUME GRAND PUBLIC : La géochimie des isotopes stables a pris beaucoup d'importance depuis ces dernières années pour aider à résoudre nombre de questions géologiques, en se basant sur les caractéristiques du fractionnement isotopiques pour différents systèmes. Il est donc nécessaire d'avoir une connaissance approfondie des mécanismes qui s'appliquent au fractionnement isotopique entre les minéraux et les fluides à partir desquels ils se forment. Ces mécanismes ont été bien approchés par différents types de calibrations pour des systèmes à hautes températures, cependant cela n'est pas aussi évident pour les systèmes à des températures inférieures à 400-500°C. Ce travail de thèse a pour but d'aider à la description et la compréhension des phénomènes qui peuvent affecter le fractionnement isotopique à basses températures, ainsi que leurs implications, à partir de l'étude de cristaux de quartz. Le choix des échantillons s'est porté sur des cristaux naturels formés à des températures inférieures ou égales à 400°C, provenant majoritairement de fissures hydrothermales Alpines dont les conditions de formation ont déjà été déterminées. L'étude des cristaux Alpin permet de plus de replacer les résultats obtenus dans le contexte du métamorphisme Alpin au cours du Miocène (21-13 Ma). Après examen de la structure et de la morphologie des cristaux, et leur caractérisation par cathodoluminescence (CL), des analyses chimiques détaillées sur les éléments en traces pouvant entrer dans le réseau cristallin du quartz comme impuretés (i.e. Li, Na, Al, P, K, Ca, Ge, Ti), et des isotopes stables de l'oxygène, ont été menées. En fonction des conditions de croissance, la plupart des cristaux présentent des zonations, qui peuvent être facilement reliées à la distribution des éléments traces analysés par microsonde électronique, sonde ionique (SIMS) et LA-ICPMS. De fortes concentrations d'Aluminium (plusieurs milliers de parties par million atomique) ont pu être observées dans les zones les plus externes des cristaux. De plus, les concentrations en Al et en Li sont toujours corrélées; la présence d'Hydrogène déduite à partir d'analyses par FTIR suit cette même tendance. Les différentes faces des cristaux présentent des concentrations distinctes d'Al, Li et H pour des mêmes zones de croissance, avec par exemple les concentrations les plus faibles dans les zones des faces prismatiques. Cela implique la présence d'un zonage sectoriel, qui a déjà été observé principalement dans des carbonates mais jamais décrit auparavant pour des quartz. Seuls les cristaux alpins homogènes en CL dont la croissance s'est faite à plus haute température (400°C) ne présentent aucune de ces caractéristiques. Par analogie avec le zonage sectoriel en Al, élément qui se substitue au Si dans le réseau cristallin du quartz, il est possible de penser qu'un zonage sectoriel pourrait aussi s'appliquer aux isotopes de l'oxygène. Des précédentes études avaient en effet émis cette hypothèse. Nos résultats ont été obtenus à partir d'analyses à la fois in- situ par SIMS, et par extraction assistée par laser-CO2 sur des parties de quartz soigneusement séparées, et sont en accord entre les deux méthodes. Un zonage sectoriel est en effet bien présent pour les cristaux alpins, mais principalement au niveau des zones très riches en Aluminium. Cependant, il a été déterminé que la présence d'Al dans ces teneurs avait un effet plus que minimal sur le fractionnement isotopique de l'oxygène. Des différences importantes ont été observées entre les faces r & z mais aussi au sein d'une même et seule zone, indiquant que le fractionnement a pu avoir lieu en déséquilibre, ce qui est aussi visible au niveau des valeurs totalement opposées entre faces pour la dernière phase de croissance de certains cristaux. Ainsi l'association de ces résultats laisse suggérer que la présence d'un zonage sectoriel peut être liée à différents paramètres tels que le taux de croissance ou la structure de surface du cristal, mais qui n'affectent pas de la même façon l'incorporation des éléments traces et le fractionnement isotopique. La possibilité d'un zonage sectoriel est importante à prendre en compte lors de toute interprétation de données isotopiques. Les analyses des isotopes de l'oxygène effectuées par SIMS ont aussi permis de distinguer des variations importantes à petite échelle au cours de la croissance. Des mesures faites par laser CO2 sur certaines roches encaissantes, ont permis distinguer plusieurs étapes dans la croissance des minéraux et de déduire le rôle de l'encaissant et le type de fluide. En association avec de précédentes études, il a été ainsi possible de mieux contraindre la formation de ces cristaux dans le contexte alpin et la circulation de fluide au cours du métamorphisme alpin durant le Miocène.

Estudi microtermomètric dels filons de quars de les Esquerdes de Rojà, Massís del Canigó, Pirineu Oriental

Els materials pre-hercinians i el granit de Costabona estan travessats per un conjunt de filons de quars de fins a 15 km de longitud, coneguts com a 'Esquerdes de Roja'. Aquests filons estan localment afectats per un sistema de bandes milonítiques i per un altre posterior de fractures. A partir de l'anilisi microtermometrica de les inclusions fluides primàries en quars, ha estat obtinguda una temperatura de formació de l'ordre de + 280 C per la major part del quars, i una temperatura de + 250 C per als monocristalls que omplen esquerdes i cavitats. Els filons s'haurien format a partir de solucions mineralitzants que contenien ~ a + C, a++, K+ i, en menor proporció, ~ g + + .fins una salinitat total de l'ordre del 20/@ en pes de la Na CI.

L'Ouverture émotionnelle: une nouvelle approche du vécu et du traitement émotionnels

L'Ouverture émotionnelle (OE) est une nouvelle approche théorique et pratique qui s'inscrit dans le domaine du vécu et du traitement des émotions. Le modèle de l'OE se base sur l'évaluation et l'analyse de cinq dimensions élémentaires ancrées dans les théories modernes des émotions (la représentation cognitive des états émotionnels, la perception des indicateurs émotionnels corporels internes et externes, l'expression et la communication des émotions et la régulation des émotions). Ces éléments permettent de caractériser l'affectivité des personnes et de réfléchir ainsi à une prise en charge la plus adaptée possible. Cet ouvrage se veut un véritable guide de la pratique de l'Ouverture émotionnelle tant pour les professionnels de la santé mentale (psychologues, psychothérapeutes...) que pour les étudiants en psychologie. Après une présentation des bases conceptuelles de l'OE et de ses aspects développementaux, les auteurs exposent de façon détaillée les résultats empiriques du modèle dans différents domaines de la psychologie clinique et de la psychologie de la santé (troubles de la personnalité, troubles de la dépendance, troubles anxieux, phobies, troubles somatoformes ou burnout, satisfaction dans le couple, bien-être individuel). Enfin ils fournissent des pistes et conseils pour l'application pratique du modèle de l'OE en intervention psychologique et en psychothérapie. En annexe, le lecteur trouvera l'ensemble des instruments (questionnaires et instruments d'auto-observation), sur lesquels repose le modèle de l'OE et qui sont inédits à ce jour.

Fluid-rock interaction during late stages of the Alpine exhumation

Résumé de la thèseLa fracturation des roches au cours de phases compressives ou extensives est un souvent évoquée pour expliquer la circulation de fluide au sein des roches cristallines. Dans le cadre de cette thèse, la circulation des fluides lors de l'exhumation tardive des Alpes a été étudiée en utilisant deux approches différentes: analyses structurales de la déformation fragile d'une part et analyses géochimiques des roches et des minéraux (isotopes stables, datations U/Pb, thermochronologie (U-Th)/He) d'autre part. Cette approche combinée a permis de mieux comprendre l'interaction existante entre les fluides métamorphiques et les fluides météoriques, ainsi que leur interaction avec les roches encaissantes. Le travail a été effectué dans la zone Pennique du Valais suisse.La première partie était focalisée sur la déformation fragile, le but étant de définir les différents types de déformations existantes et de déterminer l'âge relatif des différentes familles de failles. Dans la région d'étude, quatre domaines ont été distingués. Chacun d'eux comportent deux types de structures fragiles, certaines sont minéralisées alors que d'autre non. Au sein de chaque domaine, la direction principale des structures minéralisées correspond à l'orientation des accidents tectoniques majeurs de la région (Aosta- Ranzola Line au Sud, Rhône Line au Nord et Simplon Fault Zone à l'Est), alors que les structures non- minéralisées montrent des orientations plus variables. Ainsi, le premier type de structure est interprété comme résultant d'une dislocation tectonique alors que le deuxième type de structure résulterait d'une dislocation gravitaire locale. Il n'est néanmoins pas possible de classer chronologiquement la formation de ces deux types de structure ni d'attribuer un âge relatif aux changements d'orientation des contraintes majeures.La deuxième étude a été effectuée dans la région de la zone de faille du Simplon. Dans cette zone, la composition isotopique des minéraux ayant cristallisé à l'intérieur des fractures tardives permet de distinguer différents types de circulation de fluide. Les valeurs δ180 du quartz de la roche encaissante ainsi que ceux des veines tardives du bloque inférieur de la faille sont comparables. Ces valeurs indiquent un rééquilibrage et un tamponnage isotopique des fluides tardifs au contact de la roche encaissante lors de la fracturation de cette dernière et de la cristallisation des veines tardives. La même situation est observée dans la partie nord du bloque supérieur ainsi que dans sa partie sud. Ceci n'est néanmoins pas le cas pour la partie centrale du bloque supérieur où les valeurs isotopiques des minéraux dans les veines tardives sont approximativement 3 %o plus basses (avec des valeurs extrêmes négatifs), indiquant une contribution d'eau météorique aux fluides circulant dans les veines. Ces données suggèrent qu'une infiltration d'eau météorique a pu avoir lieu dans le bloque supérieur, où la fracturation des roches est plus intensive car le déplacement relatif le long de la faille y fut plus important, et la température maximale du métamorphisme plus basse. La troisième contribution traite de la géo-thermochronologie de la zone de contact entre la klippe de la Dent Blanche et la nappe de Tsaté. De petits zircons euhédraux ont été trouvés dans un plan de faille minéralisé (parallèle à la Faille du Rhône, voir première partie de l'étude), riche en hématite et quartz, de la zone d'étude. Les analyses U/Pb donnent des âges radiométriques autour de 270 - 280 Ma aux zircons extraits de la minéralisation ainsi que ceux extraits de la roche encaissante, ce qui correspond à l'âge de la nappe de la Dent Blanche et non celui de la nappe du Tsaté qui est elle-même classiquement interprétée comme une ophiolite Jurassique de l'Océan Liguro-Piémontais. Ces données suggèrent que les zircons contenus dans la veine ont été hérités de la roche encaissante. Les résultats (U-Th)/He indiquent un âge de refroidissement différent pour la roche encaissante (25.5 ± 2.0 Ma) que celui de la minéralisation (17.7 ±1.4 Ma). Le thermomètre isotopique quartz-hématite indique une température d'équilibre, et donc de mise en place de la minéralisation, d'environ 170 °C, température très proche de la température de -180 °C de fermeture du zircon pour le système (U-Th)/He. Ceci suggère que l'âge de refroidissement des zircons de la minéralisation correspond aussi à l'âge de formation de la faille.Thesis abstractFluid circulation in fractured rocks is a common process in geology, and it is generally the consequence of faulting and fracturing during both tectonic compression and extension. This thesis is focused on fluid circulation during late stages of the Alpine exhumation. After a structural analysis of the late brittle deformation of the studied samples, several analytical methods (stable isotope investigations, U/Pb radiometric dating, (U-Th)/He thermochronology) have been applied to understand the interaction of metamorphic and meteoric fluids with one another as well as with the host rock. This thesis is articulated around three study directions. All studies were conducted in the Penninic Zone of the Valais, Switzerland. The first study deals with late, brittle deformation and focuses on the different deformation styles and on the relative age of the different families of fractures. In order to do this, late brittle structures observed in four different domains have been subdivided as a function of the existence (or not) and type of mineralization. Comparisons between mineralized and non-mineralized strike directions for all four domains show that mineralized structures follow the strike orientation of major tectonic movements indicated in the Penninic Zone of the Valais (Aosta-Ranzola Line to the S, Rhône Line to the Ν and Simplon Fault Zone to the E), whereas non-mineralized fractures have a more variable strike orientation. This difference could be interpreted as indicative of tectonic-related faulting (mineralized structures) vs. local, collapse-related faulting (non-mineralized fractures), but it is not strong enough to indicate a relative age of the late brittle structures, and/or a change in the orientation of the strain field in post-Miocene times. The second studied area is focused on the Simplon Fault Zone (SFZ). Stable isotope analyses of minerals filling these late fractures indicate that there are two different fluid circulation systems in the footwall and hanging wall of the SFZ. In the footwall, δ180 values of quartz from both the host rock and the late veins range from +10 %o to +12 %o. This is consistent with buffering of circulating fluids by the host rock during fracturing and vein precipitation. In the hanging wall, δΙ80 values for quartz crystals from the host rock and the late veins are similar in both the northern and southern parts of the detachment that are both affected by the same degree of metamorphism (greenschist to the Ν and amphibolite to the S). This is not the case in the central part of the SFZ, where there is a jump from amphibolite facies in the footwall to greenschist facies in the hanging wall. δ,80 values for quartz from the hanging wall late veins are approximately 3.0 %o lower (down to negative values in some cases) than the values observed in the footwall These data suggest that infiltration of meteoric water may have occurred in the most fractured parts of the hanging wall, where relative displacement on the SFZ was the greatest and the peak temperature lower. In the less fractured footwall the δ180 values reflect a host rock-buffered system.The third study is focused on geo-thermochronology at the contact between the Dent Blanche nappe and the Tsaté nappe where small, euhedral zircons were found in a hematite- and quartz-rich mineralization on a late normal fault plane parallel to the Rhône Line (see first part of the study). U/Pb analysis indicates that the zircons - both in the late mineralization and in the host rock - have absolute radiometric ages clustering around 270 - 280 Ma, which is the accepted age for intrusive rocks from the Austroalpine Dent Blanche units but not for the Tsaté nappe. The latter is classically interpreted as an ophiolitic remnant of the Jurassic Liguro-Piemontais Ocean. U/Pb analyses suggest that zircons in late mineralization are all inherited from the host rock; however, results of (U-Th)/He analyses indicate that cooling ages for the host rocks are different to the cooling ages for the zircons in late mineralization. Indeed, the calculated cooling age for the Arolla gneiss is 25.5 ± 2.0 Ma, whilst the cooling age for the associated mineralized fault plane is 17.7 ±1.4 Ma. Oxygen stable isotope fractionation between quartz and hematite in the same late mineralization corresponds to temperatures of about 170 °C. The proximity of the calculated emplacement temperature for the mineralization and the lower accepted closure temperature for zircon in the (U-Th)/He system (-180 °C) imply that the age of 17.7 ± 1.4 Ma can also be interpreted as the formation age of this late brittle fault.Résumé grand publicLa circulation des fluides dans les roches fracturées est typique de nombreux processus géologiques, et très souvent est la conséquence de la fracturation des roches. Cette thèse aborde la question de la circulation des fluides pendant les dernières phases du soulèvement des Alpes. Après une analyse structurale de la fracturation directement sur le terrain, plusieurs méthodes géochimiques ont été appliquées pour comprendre l'interaction entre les différents fluides circulants, et avec leur propre roche mère. L'étude, concentrée sur trois directions principales, a été conduite dans la zone Pennique du Valais suisse. La première partie traite de la déformation cassante dans le secteur cité. L'analyse détaillée des fractures a permis de les subdiviser en structures minéralisées et non-minéralisées, sur quatre domaines différents. La comparaison entre les directions des structures minéralisées et non-minéralisées a permis de montrer que les premières suivent l'orientation des accidents tectoniques majeurs de la région, alors que les structures non- minéralisées ont une orientation plus variable. Cette différence pourrait être interprétée comme indication d'une dislocation tectonique (structures minéralisées) contre une dislocation gravitaire locale (structures non-minéralisées), mais elle n'est pas assez forte pour indiquer un âge relatif des structures tardives et/ou un changement de l'orientation des contraintes après -20 Ma vers le présent.A partir de ces observations, la deuxième étude est concentrée dans la région de la faille du Simplon. Les analyses géochimiques sur les minéraux remplissant les structures tardives indiquent qu'il y a deux différents systèmes de circulation des fluides dans les deux parties (toit et mur) de la faille. Dans le mur, les valeurs isotopiques des minéraux cristallisés à partir d'un fluide tardif sont les mêmes de ceux de la roche mère, donc il y a eu rééquilibration chimique entre fluide et roche pendant la fracturation de cette dernière et la précipitation des minéraux. Dans le toit, les valeurs isotopiques dans la roche mère et dans les minéraux des veines tardives sont comparables dans les parties Ν et S de la faille, où les roches du toit et du mur ont atteint une température maximale - pendant phase prograde de la formation des Alpes - comparable. Au contraire, dans la partie centrale, où le mur a atteint des températures maximales plus élevées par rapport au toit, les valeurs géochimiques des minéralisations tardives du toit sont parfois plus basses que les valeurs observées dans le mur. Ces données suggèrent que l'infiltration de l'eau de surface aurait pu se produire dans la partie plus fracturée du toit, où le déplacement relatif le long de la faille était majeur et les températures maximales mineures. Au contraire, les données géochimiques du mur de la partie centrale indiquent un système isotopique équilibré par la roche mère.La troisième partie de ce travail se base sur l'étude géochimique intégrée des isotopes stables d'Oxygène et radioactifs du Plomb, Uranium, Thorium et Hélium, auprès d'une faille normale minéralisée et des roches de la région à cheval entre deux nappes, la nappe de la Dent Blanche et la nappe de Tsaté. Ici, des petits zircons ont été trouvés dans la minéralisation citée, riche en hématite et quartz. L'analyse radiométrique Uranium/Plomb a montré que les zircons dans la minéralisation et dans les roches autour ont des âges comparables (autour 280 Ma). Cela signifie que les zircons dans la minéralisation tardive ont été hérités de la roche mère pendant la fracturation et la circulation des fluides tardives. De l'autre coté, les résultats des analyses Uranium-Thorium/Hélium indiquent que les âges de refroidissement pour les roches mères sont différents comparés aux âges de refroidissement pour les zircons dans la minéralisation tardive: ces derniers sont plus jeunes d'environ 8 Ma (autour 25 Ma et autour 17 Ma respectivement). Les analyses des isotopes de l'oxygène sur quartz et hématite dans la même minéralisation donnent une température de mise en place de cette dernière d'environ 170° C. La température de fermeture du système chimique des zircons dans le système (Uranium-Thorium)/Hélium est d'environ 180 °C: la proximité de ces deux températures implique que l'âge de refroidissement de la minéralisation tardive peut également être interprété comme âge de formation de la faille.

Pore-scale modeling of two-phase flow instabilities in porous media

Les problèmes d'écoulements multiphasiques en média poreux sont d'un grand intérêt pour de nombreuses applications scientifiques et techniques ; comme la séquestration de C02, l'extraction de pétrole et la dépollution des aquifères. La complexité intrinsèque des systèmes multiphasiques et l'hétérogénéité des formations géologiques sur des échelles multiples représentent un challenge majeur pour comprendre et modéliser les déplacements immiscibles dans les milieux poreux. Les descriptions à l'échelle supérieure basées sur la généralisation de l'équation de Darcy sont largement utilisées, mais ces méthodes sont sujettes à limitations pour les écoulements présentant de l'hystérèse. Les avancées récentes en terme de performances computationnelles et le développement de méthodes précises pour caractériser l'espace interstitiel ainsi que la distribution des phases ont favorisé l'utilisation de modèles qui permettent une résolution fine à l'échelle du pore. Ces modèles offrent un aperçu des caractéristiques de l'écoulement qui ne peuvent pas être facilement observées en laboratoire et peuvent être utilisé pour expliquer la différence entre les processus physiques et les modèles à l'échelle macroscopique existants. L'objet premier de la thèse se porte sur la simulation numérique directe : les équations de Navier-Stokes sont résolues dans l'espace interstitiel et la méthode du volume de fluide (VOF) est employée pour suivre l'évolution de l'interface. Dans VOF, la distribution des phases est décrite par une fonction fluide pour l'ensemble du domaine et des conditions aux bords particulières permettent la prise en compte des propriétés de mouillage du milieu poreux. Dans la première partie de la thèse, nous simulons le drainage dans une cellule Hele-Shaw 2D avec des obstacles cylindriques. Nous montrons que l'approche proposée est applicable même pour des ratios de densité et de viscosité très importants et permet de modéliser la transition entre déplacement stable et digitation visqueuse. Nous intéressons ensuite à l'interprétation de la pression capillaire à l'échelle macroscopique. Nous montrons que les techniques basées sur la moyenne spatiale de la pression présentent plusieurs limitations et sont imprécises en présence d'effets visqueux et de piégeage. Au contraire, une définition basée sur l'énergie permet de séparer les contributions capillaires des effets visqueux. La seconde partie de la thèse est consacrée à l'investigation des effets d'inertie associés aux reconfigurations irréversibles du ménisque causé par l'interface des instabilités. Comme prototype pour ces phénomènes, nous étudions d'abord la dynamique d'un ménisque dans un pore angulaire. Nous montrons que, dans un réseau de pores cubiques, les sauts et reconfigurations sont si fréquents que les effets d'inertie mènent à différentes configurations des fluides. A cause de la non-linéarité du problème, la distribution des fluides influence le travail des forces de pression, qui, à son tour, provoque une chute de pression dans la loi de Darcy. Cela suggère que ces phénomènes devraient être pris en compte lorsque que l'on décrit l'écoulement multiphasique en média poreux à l'échelle macroscopique. La dernière partie de la thèse s'attache à démontrer la validité de notre approche par une comparaison avec des expériences en laboratoire : un drainage instable dans un milieu poreux quasi 2D (une cellule Hele-Shaw avec des obstacles cylindriques). Plusieurs simulations sont tournées sous différentes conditions aux bords et en utilisant différents modèles (modèle intégré 2D et modèle 3D) afin de comparer certaines quantités macroscopiques avec les observations au laboratoire correspondantes. Malgré le challenge de modéliser des déplacements instables, où, par définition, de petites perturbations peuvent grandir sans fin, notre approche numérique apporte de résultats satisfaisants pour tous les cas étudiés. - Problems involving multiphase flow in porous media are of great interest in many scientific and engineering applications including Carbon Capture and Storage, oil recovery and groundwater remediation. The intrinsic complexity of multiphase systems and the multi scale heterogeneity of geological formations represent the major challenges to understand and model immiscible displacement in porous media. Upscaled descriptions based on generalization of Darcy's law are widely used, but they are subject to several limitations for flow that exhibit hysteric and history- dependent behaviors. Recent advances in high performance computing and the development of accurate methods to characterize pore space and phase distribution have fostered the use of models that allow sub-pore resolution. These models provide an insight on flow characteristics that cannot be easily achieved by laboratory experiments and can be used to explain the gap between physical processes and existing macro-scale models. We focus on direct numerical simulations: we solve the Navier-Stokes equations for mass and momentum conservation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. In the VOF the distribution of the phases is described by a fluid function (whole-domain formulation) and special boundary conditions account for the wetting properties of the porous medium. In the first part of this thesis we simulate drainage in a 2-D Hele-Shaw cell filled with cylindrical obstacles. We show that the proposed approach can handle very large density and viscosity ratios and it is able to model the transition from stable displacement to viscous fingering. We then focus on the interpretation of the macroscopic capillary pressure showing that pressure average techniques are subject to several limitations and they are not accurate in presence of viscous effects and trapping. On the contrary an energy-based definition allows separating viscous and capillary contributions. In the second part of the thesis we investigate inertia effects associated with abrupt and irreversible reconfigurations of the menisci caused by interface instabilities. As a prototype of these phenomena we first consider the dynamics of a meniscus in an angular pore. We show that in a network of cubic pores, jumps and reconfigurations are so frequent that inertia effects lead to different fluid configurations. Due to the non-linearity of the problem, the distribution of the fluids influences the work done by pressure forces, which is in turn related to the pressure drop in Darcy's law. This suggests that these phenomena should be taken into account when upscaling multiphase flow in porous media. The last part of the thesis is devoted to proving the accuracy of the numerical approach by validation with experiments of unstable primary drainage in a quasi-2D porous medium (i.e., Hele-Shaw cell filled with cylindrical obstacles). We perform simulations under different boundary conditions and using different models (2-D integrated and full 3-D) and we compare several macroscopic quantities with the corresponding experiment. Despite the intrinsic challenges of modeling unstable displacement, where by definition small perturbations can grow without bounds, the numerical method gives satisfactory results for all the cases studied.

Textural controls on induced polarization and permeability in granular media

The spatial resolution visualized with hydrological models and the conceptualized images of subsurface hydrological processes often exceed resolution of the data collected with classical instrumentation at the field scale. In recent years it was possible to increasingly diminish the inherent gap to information from point like field data through the application of hydrogeophysical methods at field-scale. With regards to all common geophysical exploration techniques, electric and electromagnetic methods have arguably to greatest sensitivity to hydrologically relevant parameters. Of particular interest in this context are induced polarisation (IP) measurements, which essentially constrain the capacity of a probed subsurface region to store an electrical charge. In the absence of metallic conductors the IP- response is largely driven by current conduction along the grain surfaces. This offers the perspective to link such measurements to the characteristics of the solid-fluid-interface and thus, at least in unconsolidated sediments, should allow for first-order estimates of the permeability structure.¦While the IP-effect is well explored through laboratory experiments and in part verified through field data for clay-rich environments, the applicability of IP-based characterizations to clay-poor aquifers is not clear. For example, polarization mechanisms like membrane polarization are not applicable in the rather wide pore-systems of clay free sands, and the direct transposition of Schwarz' theory relating polarization of spheres to the relaxation mechanism of polarized cells to complex natural sediments yields ambiguous results.¦In order to improve our understanding of the structural origins of IP-signals in such environments as well as their correlation with pertinent hydrological parameters, various laboratory measurements have been conducted. We consider saturated quartz samples with a grain size spectrum varying from fine sand to fine gravel, that is grain diameters between 0,09 and 5,6 mm, as well as corresponding pertinent mixtures which can be regarded as proxies for widespread alluvial deposits. The pore space characteristics are altered by changing (i) the grain size spectra, (ii) the degree of compaction, and (iii) the level of sorting. We then examined how these changes affect the SIP response, the hydraulic conductivity, and the specific surface area of the considered samples, while keeping any electrochemical variability during the measurements as small as possible. The results do not follow simple assumptions on relationships to single parameters such as grain size. It was found that the complexity of natural occurring media is not yet sufficiently represented when modelling IP. At the same time simple correlation to permeability was found to be strong and consistent. Hence, adaptations with the aim of better representing the geo-structure of natural porous media were applied to the simplified model space used in Schwarz' IP-effect-theory. The resulting semi- empiric relationship was found to more accurately predict the IP-effect and its relation to the parameters grain size and permeability. If combined with recent findings about the effect of pore fluid electrochemistry together with advanced complex resistivity tomography, these results will allow us to picture diverse aspects of the subsurface with relative certainty. Within the framework of single measurement campaigns, hydrologiste can than collect data with information about the geo-structure and geo-chemistry of the subsurface. However, additional research efforts will be necessary to further improve the understanding of the physical origins of IP-effect and minimize the potential for false interpretations.¦-¦Dans l'étude des processus et caractéristiques hydrologiques des subsurfaces, la résolution spatiale donnée par les modèles hydrologiques dépasse souvent la résolution des données du terrain récoltées avec des méthodes classiques d'hydrologie. Récemment il est possible de réduire de plus en plus cet divergence spatiale entre modèles numériques et données du terrain par l'utilisation de méthodes géophysiques, notamment celles géoélectriques. Parmi les méthodes électriques, la polarisation provoquée (PP) permet de représenter la capacité des roches poreuses et des sols à stocker une charge électrique. En l'absence des métaux dans le sous-sol, cet effet est largement influencé par des caractéristiques de surface des matériaux. En conséquence les mesures PP offrent une information des interfaces entre solides et fluides dans les matériaux poreux que nous pouvons lier à la perméabilité également dirigée par ces mêmes paramètres. L'effet de la polarisation provoquée à été étudié dans différentes études de laboratoire, ainsi que sur le terrain. A cause d'une faible capacité de polarisation des matériaux sableux, comparé aux argiles, leur caractérisation par l'effet-PP reste difficile a interpréter d'une manière cohérente pour les environnements hétérogènes.¦Pour améliorer les connaissances sur l'importance de la structure du sous-sol sableux envers l'effet PP et des paramètres hydrologiques, nous avons fait des mesures de laboratoire variées. En détail, nous avons considéré des échantillons sableux de quartz avec des distributions de taille de grain entre sables fins et graviers fins, en diamètre cela fait entre 0,09 et 5,6 mm. Les caractéristiques de l'espace poreux sont changées en modifiant (i) la distribution de taille des grains, (ii) le degré de compaction, et (iii) le niveau d'hétérogénéité dans la distribution de taille de grains. En suite nous étudions comment ces changements influencent l'effet-PP, la perméabilité et la surface spécifique des échantillons. Les paramètres électrochimiques sont gardés à un minimum pendant les mesures. Les résultats ne montrent pas de relation simple entre les paramètres pétro-physiques comme par exemples la taille des grains. La complexité des media naturels n'est pas encore suffisamment représenté par les modèles des processus PP. Néanmoins, la simple corrélation entre effet PP et perméabilité est fort et consistant. En conséquence la théorie de Schwarz sur l'effet-PP a été adapté de manière semi-empirique pour mieux pouvoir estimer la relation entre les résultats de l'effet-PP et les paramètres taille de graines et perméabilité. Nos résultats concernant l'influence de la texture des matériaux et celles de l'effet de l'électrochimie des fluides dans les pores, permettront de visualiser des divers aspects du sous-sol. Avec des telles mesures géo-électriques, les hydrologues peuvent collectionner des données contenant des informations sur la structure et la chimie des fluides des sous-sols. Néanmoins, plus de recherches sur les origines physiques de l'effet-PP sont nécessaires afin de minimiser le risque potentiel d'une mauvaise interprétation des données.

« Représentations de la maltraitance infantile d'un groupe d'enseignantes d'une école primaire du canton de Vaud, leur rôle dans ce domaine et leurs connaissances sur les procédures de signalement. »

Le rapport international le plus récent concernant la maltraitance infantile date de 2006 : il s'agit du Rapport mondial sur la violence contre les enfants, du Secrétaire général des Nations Unies (1). La définition retenue pour la maltraitance infantile s'inspire de celle du Rapport mondial sur la violence et la santé, de l'OMS en 2002 (2) : «La menace ou l'utilisation intentionnelle de la force physique ou du pouvoir contre un enfant par un individu ou un groupe qui entraîne ou risque fortement de causer un préjudice à la santé, à la survie, au développement ou à la dignité de l'enfant.». Il existe différentes formes de maltraitance : - la maltraitance physique (brutalités, coups, blessures, brûlures, etc.)  la maltraitance psychologique (insultes, humiliation, isolement, terroriser l'enfant, etc.) - la maltraitance sexuelle (exhibitionnisme, attouchements, relations sexuelles, etc.) - les négligences (manque d'attention et de soins) Dans la majorité des cas, plusieurs formes de maltraitances sont présentes chez un enfant victime de mauvais traitements ; elles se chevauchent (3). L'Observatoire national de l'Action Sociale Décentralisée (ODAS) a réalisé une classification des enfants à protéger, les définitions sont les suivantes (4): L'enfant maltraité est « celui qui est victime de violences physiques, d'abus sexuels, de cruauté mentale, de négligences lourdes ayant des conséquences sur son développement physique et psychologique. » L'enfant en risque est « celui qui connaît des conditions d'existence qui risquent de mettre en danger sa santé, sa sécurité, sa moralité, son éducation ou son entretien, mais qui n'est pas pour autant maltraité. » L'enfant en souffrance est « un enfant aimé et soigné mais qui souffre des conditions d'existences qui fragilisent ou menacent son développement et son épanouissement personnel. » En Suisse, peu de données sont disponibles concernant la prévalence de la maltraitance étant donné la difficulté à récolter des données. Selon l'Office Fédéral de la Statistique suisse, les résultats d'une étude de 2004 montre une diminution des châtiments corporels par rapport à une étude semblable réalisée 12 ans auparavant (5). Cependant, la maltraitance infantile est un problème de santé publique du fait de la gravité de ses conséquences sur la santé physique, mentale et sociale de l'individu et de son retentissement sur la communauté ainsi que de sa fréquence estimée dans la population suisse. Elle a des effets néfastes sur la santé de l'enfant par mortalité directe ou morbidité directe ou indirecte et représente également un facteur de risque pour la santé physique et mentale, le développement et les perspectives de réalisation personnelle du jeune adulte et de l'adulte (6). On sait aujourd'hui que le nombre de cas de maltraitance signalés en Suisse est en augmentation. Ceci démontre que la maltraitance est un phénomène courant. Cependant, les professionnels ne pensent pas MF / Travail de master en médecine / 2011-2012 3 que le phénomène de la maltraitance infantile soit en augmentation, mais que les cas de maltraitance sont mieux repérés, que les professionnels s'occupant d'enfants sont plus sensibles à cette problématique et qu'il y a donc davantage de signalements (7). La prévention de la maltraitance est nécessaire et possible. Des interventions ont établi leur efficacité et il a été démontré que plus l'intervention est précoce, plus elle a de chances de réussite (2). C'est la raison pour laquelle il est important de repérer les cas de maltraitance précocement afin de pouvoir intervenir, aider les familles et garantir la protection de l'enfant. Des mesures de prévention ont été mises en place au niveau international, comme au niveau fédéral, pour assurer la reconnaissance et la prise en charge de l'enfant victime de maltraitance. Au niveau international, la Convention internationale des droits de l'enfant a été adoptée par l'Assemblée Générale en 1989 (8). Elle reconnaît l'enfant comme personne indépendante ayant des droits propres. Cette convention est divisée en quatre parties comportant : les principes directeurs (la non-discrimination, viser les meilleurs intérêts pour l'enfant, le droit de vivre, de survivre et de se développer, le droit de participation), les droits de survie et de développement (le droit à avoir les ressources, les compétences et les contributions nécessaires pour pouvoir survivre et pouvoir profiter d'un développement complet), les droits de protection (de toutes les formes de maltraitance envers les enfants, négligences, exploitation et cruauté), les droits de participation (la liberté d'expression de leurs opinions, de parler de sujets qui concernent leur vie sociale, économique, religieuse, culturelle ou politique et d'être écouté, la liberté d'information et la liberté d'association). Les stratégies de prévention de la maltraitance infantile visent à réduire les causes sous- jacentes et les facteurs de risque, tout en renforçant les facteurs de protection, de façon à prévenir de nouveaux cas (9). Elles comprennent : les stratégies sociétales et communautaires (mise en place de réformes juridiques et des droits de la personne humaine, instauration des politiques sociales et économiques favorables, correction des normes sociales et culturelles, réduction des inégalités économiques, réduction du facteur de risque environnemental, formation des professionnels), les stratégies relationnelles (formation parentale et des adultes s'occupant d'enfants), les stratégies individuelles (apprendre aux enfants à reconnaître et à éviter les situations de violence potentielle). En plus des mesures structurelles mises en place par les états (scolarisation obligatoire, dispositif légal, service de protection des enfants et des jeunes, services de santé spécialisés, etc.), des associations de lutte contre la maltraitance infantile existent et jouent également un rôle important dans la prévention. Par exemple, la Fondation Suisse pour la Protection de l'Enfant s'emploie à analyser les causes de la violence envers les MF / Travail de master en médecine / 2011-2012 4 enfants et à les combattre, à protéger les enfants contre la violence physique, psychologique, sexuelle et structurelle ainsi que contre la négligence par le biais d'un travail de prévention ciblé à l'échelle nationale. Elle vise également à apprendre aux enfants comment se protéger eux-mêmes et demander de l'aide, à sensibiliser les adultes qui les entourent au fait que les enfants ont une personnalité propre et qu'ils ont le droit d'être protégés et encouragés et à demander au niveau politique que l'on mette en place des structures adaptées aux enfants (10).