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.

An 40Ar/39Ar, Rb/Sr, and stable isotope study of micas in low-grade fold-and-thrust belt: An example from the Swiss Helvetic Alps

White micas in carbonate-rich tectonites and a few other rock types of large thrusts in the Swiss Helvetic fold-and-thrust belt have been analyzed by Ar-40/Ar-39 and Rb/Sr techniques to better constrain the timing of Alpine deformation for this region. Incremental Ar-40/Ar-39 heating experiments of 25 weakly metamorphosed (anchizone to low greenschist) samples yield plateau and staircase spectra. We interpret most of the staircase release spectra result from variable mixtures of syntectonic (neoformed) and detrital micas. The range in dates obtained within individual spectra depends primarily on the duration of mica nucleation and growth, and relative proportions of neoformed and detrital mica. Rb/Sr analyses of 12 samples yield dates of ca. 10-39 Ma (excluding one anomalously young sample). These dates are slightly younger than the Ar-40/Ar-39 total gas dates obtained for the same samples. The Rb/ Sr dates were calculated using initial Sr-87/Sr-86 ratios obtained from the carbonate-dominated host rocks, which are higher than normal Mesozoic carbonate values due to exchange with fluids of higher Sr-87/Sr-86 ratios (and lower O-18/O-16 ratios). Model dates calculated using Sr-87/Sr-86 values typical of Mesozoic marine carbonates more closely approximate the Ar-40/Ar-39 total gas dates for most of the samples. The similarities of Rb/Sr and Ar-40/Ar-39 total gas dates are consistent with limited amounts of detrital mica in the samples. The delta(18)O values range from 24-15%. (VSMOW) for 2-6 mum micas and 27-16parts per thousand for the carbonate host rocks. The carbonate values are significantly lower than their protolith values due to localized fluid-rock interaction and fluid flow along most thrust surfaces. Although most calcite-mica pairs are not in oxygen isotope equilibrium at temperatures of ca. 200-400 degreesC, their isotopic fractionations are indicative of either 1) partial exchange between the minerals and a common external fluid, or 2) growth or isotopic exchange of the mica with the carbonate after the carbonate had isotopically exchanged with an external fluid. The geological significance of these results is not easily or uniquely determined, and exemplifies the difficulties inherent in dating very fine-grained micas of highly deformed tectonites in low-grade metamorphic terranes. Two generalizations can be made regarding the dates obtained from the Helvetic thrusts: 1) samples from the two highest thrusts (Mt. Gond and Sublage) have all of their Ar-40/Ar-39 steps above 20 Ma, and 2) most samples from the deepest Helvetic thrusts have steps (often accounting for more than 80% of Ar-39 release) between 15 and 25 Ma. These dates are consistent with the order of thrusting in the foreland-imbricating system and increase proportions of neoformed to detrital mica in the more metamorphosed hinterland and deeply buried portions of the nappe pile. Individual thrusts accommodated the majority of their displacement during their initial incorporation into the foreland-imbricating system, and some thrusts remained active or were reactivated down to 15 Ma.

The role of a transcrustal shear zone in orogenic gold mineralization at the Aijanahalli Mine, Dharwar Craton, South India

The Ajjanahalli gold mine is spatially associated with a Late Archean craton-scale shear zone in the eastern Chitradurga greenstone belt of the Dharwar craton, India. Gold mineralization is hosted by an similar to100-m-wide antiform in a banded iron formation. Original magnetite and siderite are replaced by a peak metamorphic alteration assemblage of chlorite, stilpnomelane, minnesotaite, sericite, ankerite, arsenopyrite, pyrite, pyrrhotite, and gold at ca. 300degrees to 350degreesC. Elements enriched in the banded iron formation include Ca, Mg, C, S, An, As, Bi. Cu, Sb, Zn, Pb, Se, Ag, and Te, whereas in the wall rocks As, Cu, Zn, Bi, Ag, and An are only slightly enriched. Strontium correlates with CaO, MgO, CO2, and As, which indicates cogenetic formation of arsenopyrite and Mg-Ca carbonates. The greater extent of alteration in the Fe-rich banded iron formation layers than in the wall rock reflects the greater reactivity of the banded iron formation layers. The ore fluids, as interpreted from their isotopic composition (delta(18)O = 6.5-8.5parts per thousand; initial Sr-87/Sr-86 = 0.7068-0.7078), formed by metamorphic devolatilization of deeper levels of the Chitradurga greenstone belt. Arsenopyrite, chalcopyrite, and pyrrhotite have delta(34)S values within a narrow range between 2.1 and 2.7 per mil, consistent with a sulfur source in Chitradurga greenstone belt lithologies. Based on spatial and temporal relationships between mineralization, local structure development, and sinistral strike-slip deformation in the shear zone at the eastern contact of the Chitradurga greenstone belt, we suggest that the Ajjanahalli gold mineralization formed by fluid infiltration into a low strain area within the first-order structure. The ore fluids were transported along this shear zone into relatively shallow crustal levels during lateral terrane accretion and a change from thrust to transcurrent tectonics. Based on this model of fluid flow, exploration should focus on similar low strain areas or potentially connected higher order splays of the first-order shear zone.

Stable isotope geochemistry and formation mechanisms of quartz veins; Extreme paleoaltitudes of the Central Alps in the Neogene

Quartz veins ranging in size from less than 50 cm length and 5 cm width to greater than 10 m in length and 5 m in width are found throughout the Central Swiss Alps. In some cases, the veins are completely filled with milky quartz, while in others, sometimes spectacular void-filling quartz crystals are found. The style of vein filling and size is controlled by host rock composition and deformation history. Temperatures of vein formation, estimated using stable isotope thermometry and mineral equilibria, cover a range of 450 degrees C down to 150 degrees C. Vein formation started at 18 to 20 Ma and continued for over 10 My. The oxygen isotope values of quartz veins range from 10 to 20 permil, and in almost all cases are equal to those of the hosting lithology. The strongly rock-buffered veins imply a low fluid/rock ratio and minimal fluid flow. In order to explain massive, nearly morromineralic quartz formation without exceptionally large fluid fluxes, a mechanism of differential pressure and silica diffusion, combined with pressure solution, is proposed for early vein formation. Fluid inclusions and hydrous minerals in late-formed veins have extremely low delta D values, consistent with meteoric water infiltration. The change from rock-buffered, static fluid to infiltration from above can be explained in terms of changes in the large-scale deformation style occurring between 20 and 15 Ma. The rapid cooling of the Central Alps identified in previous studies may be explained in part, by infiltration of cold meteoric waters along fracture systems down to depths of 10 km or more. An average water flux of 0.15 cm 3 cm(-2)yr(-1) entering the rock and reemerging heated by 40 degrees C is sufficient to cool rock at 10 km depth by 100 degrees C in 5 million years. The very negative delta D values of < -130 permil for the late stage fluids are well below the annual average values measured in meteoric water in the region today. The low fossil delta D values indicate that the Central Alps were at a higher elevation in the Neogene. Such a conclusion is supported by an earlier work, where a paleoaltitude of 5000 meters was proposed on the basis of large erratic boulders found at low elevations far from their origin.

Devolatilization-induced pressure build-up: Implications for reaction front movement and breccia pipe formation

Generation of fluids during metamorphism can significantly influence the fluid overpressure, and thus the fluid flow in metamorphic terrains. There is currently a large focus on developing numerical reactive transport models, and with it follows the need for analytical solutions to ensure correct numerical implementation. In this study, we derive both analytical and numerical solutions to reaction-induced fluid overpressure, coupled to temperature and fluid flow out of the reacting front. All equations are derived from basic principles of conservation of mass, energy and momentum. We focus on contact metamorphism, where devolatilization reactions are particularly important owing to high thermal fluxes allowing large volumes of fluids to be rapidly generated. The analytical solutions reveal three key factors involved in the pressure build-up: (i) The efficiency of the devolatilizing reaction front (pressure build-up) relative to fluid flow (pressure relaxation), (ii) the reaction temperature relative to the available heat in the system and (iii) the feedback of overpressure on the reaction temperature as a function of the Clapeyron slope. Finally, we apply the model to two geological case scenarios. In the first case, we investigate the influence of fluid overpressure on the movement of the reaction front and show that it can slow down significantly and may even be terminated owing to increased effective reaction temperature. In the second case, the model is applied to constrain the conditions for fracturing and inferred breccia pipe formation in organic-rich shales owing to methane generation in the contact aureole.

Seismic attenuation and velocity dispersion in heterogeneous partially saturated porous rocks

Using a numerical approach, we explore wave-induced fluid flow effects in partially saturated porous rocks in which the gas-water saturation patterns are governed by mesoscopic heterogeneities associated with the dry frame properties. The link between the dry frame properties and the gas saturation is defined by the assumption of capillary pressure equilibrium, which in the presence of heterogeneity implies that neighbouring regions can exhibit different levels of saturation. To determine the equivalent attenuation and phase velocity of the synthetic rock samples considered in this study, we apply a numerical upscaling procedure, which permits to take into account mesoscopic heterogeneities associated with the dry frame properties as well as spatially continuous variations of the pore fluid properties. The multiscale nature of the fluid saturation is taken into account by locally computing the physical properties of an effective fluid, which are then used for the larger-scale simulations. We consider two sets of numerical experiments to analyse such effects in heterogeneous partially saturated porous media, where the saturation field is determined by variations in porosity and clay content, respectively. In both cases we also evaluate the seismic responses of corresponding binary, patchy-type saturation patterns. Our results indicate that significant attenuation and modest velocity dispersion effects take place in this kind of media for both binary patchy-type and spatially continuous gas saturation patterns and in particular in the presence of relatively small amounts of gas. The numerical experiments also show that the nature of the gas distribution patterns is a critical parameter controlling the seismic responses of these environments, since attenuation and velocity dispersion effects are much more significant and occur over a broader saturation range for binary patchy-type gas-water distributions. This analysis therefore suggests that the physical mechanisms governing partial saturation should be accounted for when analysing seismic data in a poroelastic framework. In this context, heterogeneities associated with the dry frame properties, which do not play important roles in wave-induced fluid flow processes per se, should be taken into account since they may determine the kind of gas distribution pattern taking place in the porous rock.

Metamorphism and kinetics in the Torres del Paine contact aureole

Contact aureoles provide an excellent geologic environment to study the mechanisms of metamorphic reactions in a natural system. The Torres del Paine (TP) intrusion is one of the most spectacular natural laboratories because of its excellent outcrop conditions. It formed in a period from 12.59 to 12.43 Ma and consists of three large granite and four smaller mafic batches. The oldest granite is on top, the youngest at the bottom of the granitic complex, and the granites overly the mafic laccolith. The TP intruded at a depth of 2-3 km into regional metamorphic anchizone to greenschist facies pelites, sandstones, and conglomerates of the Cerro Toro and Punta Barrosa formations. It formed a thin contact aureole of 150-400 m width. This thesis focuses on the reaction kinetics of the mineral cordierite in the contact aureole using quantitative textural analysis methods. First cordierite was formed from chlorite break¬down (zone I, ca. 480 °C, 750 bar). The second cordierite forming reaction was the muscovite break-down, which is accompanied by a modal decrease in biotite and the appearance of k- feldspar (zone II, 540-550 °C, 750 bar). Crystal sizes of the roundish, poikiloblastic cordierites were determined from microscope thin section images by manually marking each crystal. Images were then automatically processed with Matlab. The correction for the intersection probability of each crystal radius yields the crystal size distribution in the rock. Samples from zone I below the laccolith have the largest crystals (0.09 mm). Cordierites from zone II are smaller, with a maximum crystal radius of 0.057 mm. Rocks from zone II have a larger number of small cordierite crystals than rocks from zone I. A combination of these quantitative analysis with numerical modeling of nucleation and growth, is used to infer nucleation and growth parameters which are responsible for the observed mineral textures. For this, the temperature-time paths of the samples need to be known. The thermal history is complex because the main body of the intrusion was formed by several intrusive batches. The emplacement mechanism and duration of each batch can influence the thermal structure in the aureole. A possible subdivision of batches in smaller increments, so called pulses, will focus heat at the side of the intrusion. Focusing all pulses on one side increases the contact aureole size on that side, but decreases it on the other side. It forms a strongly asymmetric contact aureole. Detailed modeling shows that the relative thicknesses of the TP contact aureole above and below the intrusion (150 and 400 m) are best explained by a rapid emplacement of at least the oldest granite batch. Nevertheless, temperatures are significantly too low in all models, compared to observed mineral assemblages in the hornfelses. Hence, an other important thermal mechanisms needs to take place in the host rock. Clastic minerals in the immature sediments outside the contact aureole are hydrated due to small amounts of expelled fluids during contact metamorphism. This leads to a temperature increase of up to 50 °C. The origin of fluids can be traced by stable isotopes. Whole rock stable isotope data (6D and δ180) and chlorine concentrations in biotite document that the TP intrusion induced only very small amounts of fluid flow. Oxygen whole rock data show δ180 values between 9.0 and 10.0 %o within the first 5 m of the contact. Values increase to 13.0 - 15.0 %o further away from the intrusion. Whole rock 6D values display a more complex zoning. First, host rock values (-90 to -70 %o) smoothly decrease towards the contact by ca. 20 %o, up to a distance of ca. 150 m. This is followed by an increase of ca. 20 %o within the innermost 150 m of the aureole (-97.0 to -78 %o at the contact). The initial decrease in 6D values is interpreted to be due to Rayleigh fractionation accompanying the dehydration reactions forming cordierite, while the final increase reflects infiltration of water-rich fluids from the intrusion. An over-estimate on the quantity and the corresponding thermal effect yields a temperature increase of less than 30 °C. This suggests that fluid flow might have contributed only for a small amount to the thermal evolution of the system. A combination of the numerical growth model with the thermal model, including the hydration reaction enthalpies but neglecting fluid flow and incremental growth, can be used to numerically reproduce the observed cordierite textures in the contact aureole. This yields kinetic parameters which indicate fast cordierite crystallization before the thermal peak in the inner aureole, and continued reaction after the thermal peak in the outermost aureole. Only small temperature dependencies of the kinetic parameters seem to be needed to explain the obtained crystal size data. - Les auréoles de contact offrent un cadre géologique privilégié pour l'étude des mécanismes de réactions métamorphiques associés à la mise en place de magmas dans la croûte terrestre. Par ses conditions d'affleurements excellentes, l'intrusion de Torres del Paine représente un site exceptionnel pour améliorer nos connaissances de ces processus. La formation de cette intrusion composée de trois injections granitiques principales et de quatre injections mafiques de volume inférieur couvre une période allant de 12.50 à 12.43 Ma. Le plus vieux granite forme la partie sommitale de l'intrusion alors que l'injection la plus jeune s'observe à la base du complexe granitique; les granites recouvrent la partie mafique du laccolite. L'intrusion du Torres del Paine s'est mise en place a 2-3 km de profondeur dans un encaissant métamorphique. Cet encaissant est caractérisé par un métamorphisme régional de faciès anchizonal à schiste vert et est composé de pélites, de grès, et des conglomérats des formations du Cerro Toro et Punta Barrosa. La mise en place des différentes injections granitiques a généré une auréole de contact de 150-400 m d'épaisseur autour de l'intrusion. Cette thèse se concentre sur la cinétique de réaction associée à la formation de la cordiérite dans les auréoles de contact en utilisant des méthodes quantitatives d'analyses de texture. On observe plusieurs générations de cordiérite dans l'auréole de contact. La première cordiérite est formée par la décomposition de la chlorite (zone I, environ 480 °C, 750 bar), alors qu'une seconde génération de cordiérite est associée à la décomposition de la muscovite, laquelle est accompagnée par une diminution modale de la teneur en biotite et l'apparition de feldspath potassique (zone II, 540-550 °C, 750 bar). Les tailles des cristaux de cordiérites arrondies et blastic ont été déterminées en utilisant des images digitalisées des lames minces et en marquant individuellement chaque cristal. Les images sont ensuite traitées automatiquement à l'aide du programme Matlab. La correction de la probabilité d'intersection en fonction du rayon des cristaux permet de déterminer la distribution de la taille des cristaux dans la roche. Les échantillons de la zone I, en dessous du lacolite, sont caractérisés par de relativement grands cristaux (0.09 mm). Les cristaux de cordiérite de la zone II sont plus petits, avec un rayon maximal de 0.057 mm. Les roches de la zone II présentent un plus grand nombre de petits cristaux de cordiérite que les roches de la zone I. Une combinaison de ces analyses quantitatives avec un modèle numérique de nucléation et croissance a été utilisée pour déduire les paramètres de nucléation et croissance contrôlant les différentes textures minérales observées. Pour développer le modèle de nucléation et de croissance, il est nécessaire de connaître le chemin température - temps des échantillons. L'histoire thermique est complexe parce que l'intrusion est produite par plusieurs injections successives. En effet, le mécanisme d'emplace¬ment et la durée de chaque injection peuvent influencer la structure thermique dans l'auréole. Une subdivision des injections en plus petits incréments, appelés puises, permet de concentrer la chaleur dans les bords de l'intrusion. Une mise en place préférentielle de ces puises sur un côté de l'intrusion modifie l'apport thermique et influence la taille de l'auréole de contact produite, auréole qui devient asymétrique. Dans le cas de la première injection de granite, une modélisation détaillée montre que l'épaisseur relative de l'auréole de contact de Torres del Paine au-dessus et en dessous de l'intrusion (150 et 400 m) est mieux expliquée par un emplacement rapide du granite. Néanmoins, les températures calculées dans l'auréole de con¬tact sont trop basses pour que les modèles thermiques soient cohérants par rapport à la taille de cette auréole. Ainsi, un autre mecanisme exothermique est nécessaire pour permettre à la roche encais¬sante de produire les assemblages observés. L'observation des roches encaissantes entourant les granites montre que les minéraux clastiques dans les sédiments immatures au-dehors de l'auréole sont hydratés suite à la petite quantité de fluide expulsée durant le métamorphisme de contact et/ou la mise en place des granites. Les réactions d'hydratation peuvent permettre une augmentation de la température jusqu'à 50 °C. Afin de déterminer l'origine des fluides, une étude isotopique de roches de l'auréole de contact a été entreprise. Les isotopes stables d'oxygène et d'hydrogène sur la roche totale ainsi que la concentration en chlore dans la biotite indiquent que la mise en place des granites du Torres del Paine n'induit qu'une circulation de fluide limitée. Les données d'oxygène sur roche totale montrent des valeurs δ180 entre 9.0 et 10.0%o au sein des cinq premiers mètres du contact. Les valeurs augmentent jusqu'à 13.0 - 15.0 plus on s'éloigne de l'intrusion. Les valeurs 5D sur roche totale montrent une zonation plus complexe. Les valeurs de la roche encaissante (-90 à -70%o) diminuent progressivement d'environ 20%o depuis l'extérieur de l'auréole jusqu'à une distance d'environ 150 m du granite. Cette diminution est suivie par une augmentation d'environ 20%o au sein des 150 mètres les plus proches du contact (-97.0 à -78%o au contact). La diminution initiale des valeurs de 6D est interprétée comme la conséquence du fractionnement de Rayleigh qui accompagne les réactions de déshydratation formant la cordiérite, alors que l'augmentation finale reflète l'infiltration de fluide riche en eau venant de l'intrusion. A partir de ces résultats, le volume du fluide issu du granite ainsi que son effet thermique a pu être estimé. Ces résultats montrent que l'augmentation de température associée à ces fluides est limitée à un maximum de 30 °C. La contribution de ces fluides dans le bilan thermique est donc faible. Ces différents résultats nous ont permis de créer un modèle thermique associé à la for¬mation de l'auréole de contact qui intègre la mise en place rapide du granite et les réactions d'hydratation lors du métamorphisme. L'intégration de ce modèle thermique dans le modèle numérique de croissance minérale nous permet de calculer les textures des cordiérites. Cepen¬dant, ce modèle est dépendant de la vitesse de croissance et de nucléation de ces cordiérites. Nous avons obtenu ces paramètres en comparant les textures prédites par le modèle et les textures observées dans les roches de l'auréole de contact du Torres del Paine. Les paramètres cinétiques extraits du modèle optimisé indiquent une cristallisation rapide de la cordiérite avant le pic thermique dans la partie interne de l'auréole, et une réaction continue après le pic thermique dans la partie la plus externe de l'auréole. Seules de petites dépendances de température des paramètres de cinétique semblent être nécessaires pour expliquer les don¬nées obtenues sur la distribution des tailles de cristaux. Ces résultats apportent un éclairage nouveau sur la cinétique qui contrôle les réactions métamorphiques.

Detection of Optimal Models in Parameter Space with Support Vector Machines

The paper proposes an approach aimed at detecting optimal model parameter combinations to achieve the most representative description of uncertainty in the model performance. A classification problem is posed to find the regions of good fitting models according to the values of a cost function. Support Vector Machine (SVM) classification in the parameter space is applied to decide if a forward model simulation is to be computed for a particular generated model. SVM is particularly designed to tackle classification problems in high-dimensional space in a non-parametric and non-linear way. SVM decision boundaries determine the regions that are subject to the largest uncertainty in the cost function classification, and, therefore, provide guidelines for further iterative exploration of the model space. The proposed approach is illustrated by a synthetic example of fluid flow through porous media, which features highly variable response due to the parameter values' combination.

Structure et évolution thermique de la croûte supérieure des chaînes de subduction : approches thermologique et modélisation numérique, dans les North Cascades (U.S.A) et de la zone de faille Motagua (Guatemala)

Cette thèse cible l'étude de la structure thermique de la croûte supérieure (<10km) dans les arcs magmatiques continentaux, et son influence sur l'enregistrement thermochronologique de leur exhumation et de leur évolution topographique. Nous portons notre regard sur deux chaînes de montagne appartenant aux Cordillères Américaines : Les Cascades Nord (USA) et la zone de faille Motagua (Guatemala). L'approche utilisée est axée sur la thermochronologie (U-Th-Sm)/He sur apatite et zircon, couplée avec la modélisation numérique de la structure thermique de la croûte. Nous mettons en évidence la variabilité à la fois spatiale et temporelle du gradient géothermique, et attirons l'attention du lecteur sur l'importance de prendre en compte la multitude des processus géologiques perturbant la structure thermique dans les chaînes de type cordillère, c'est à dire formées lors de la subduction océanique sous un continent.Une nouvelle approche est ainsi développée pour étudier et contraindre la perturbation thermique autour des chambres magmatiques. Deux profiles âge-elevation (U-Th-Sm)/He sur apatite et zircon, ont été collectées 7 km au sud du batholithe de Chilliwack, Cascades Nord. Les résultats montrent une variabilité spatiale et temporelle du gradient géothermique lors de l'emplacement magmatique qui peut être contrainte et séparé de l'exhumation. Durant l'emplacement de l'intrusion, la perturbation thermique y atteint un état d'équilibre (-80-100 °C/km) qui est fonction du flux de magma et de ia distance à la source du magma, puis rejoint 40 °C/km à la fin du processus d'emplacement magmatique.Quelques nouvelles données (U-Th)/He, replacées dans une compilation des données existantes dans les Cascades Nord, indiquent une vitesse d'exhumation constante (-100 m/Ma) dans le temps et l'espace entre 35 Ma et 2 Ma, associée à un soulèvement uniforme de la chaîne contrôlé par l'emplacement de magma dans la croûte durant toute l'activité de l'arc. Par contre, après ~2 Ma, le versant humide de la chaîne est affecté par une accélération des taux d'exhumation, jusqu'à 3 km de croûte y sont érodés. Les glaciations ont un triple effet sur l'érosion de cette chaîne: (1) augmentation des vitesses d'érosion, d'exhumation et de soulèvement la où les précipitations sont suffisantes, (2) limitation de l'altitude contrôlé par la position de Γ Ε LA, (3) élargissement du versant humide et contraction du versant aride de la chaîne.Les modifications des réseaux de drainage sont des processus de surface souvent sous-estimés au profil d'événements climatiques ou tectoniques. Nous proposons une nouvelle approche couplant une analyse géomorphologique, des données thermochronologiques de basse température ((U-Th-Sm)/He sur apatite et zircon), et l'utilisation de modélisation numérique thermo-cinématique pour les mettre en évidence et les dater; nous testons cette approche sur la gorge de la Skagit river dans les North Cascades.De nouvelles données (U-Th)/He sur zircons, complétant les données existantes, montrent que le déplacement horizontal le long de la faille transformante continentale Motagua, la limite des plaques Caraïbe/Amérique du Nord, a juxtaposé un bloc froid, le bloc Maya (s.s.), contre un bloque chaud, le bloc Chortis (s.s.) originellement en position d'arc. En plus de donner des gammes d'âges thermochronologiques très différents des deux côtés de la faille, le déplacement horizontal rapide (~2 cm/a) a produit un fort échange thermique latéral, résultant en un réchauffement du côté froid et un refroidissement du côté chaud de la zone de faille de Motagua.Enfin des données (U-Th-Sm)/He sur apatite témoignent d'un refroidissement Oligocène enregistré uniquement dans la croûte supérieure de la bordure nord de la zone de faille Motagua. Nous tenterons ultérieurement de reproduire ce découplage vertical de la structure thermique par la modélisation de la formation d'un bassin transtensif et de circulation de fluides le long de la faille de Motagua. - This thesis focuses on the influence of the dynamic thermal structure of the upper crust (<10km) on the thermochronologic record of the exhumational and topographic history of magmatic continental arcs. Two mountain belts from the American Cordillera are studied: the North Cascades (USA) and the Motagua fault zone (Guatemala). I use a combined approach coupling apatite and zircon (U-Th-Sm}/He thermochronology and thermo- kinematic numerical modelling. This study highlights the temporal and spatial variability of the geothermal gradient and the importance to take into account the different geological processes that perturb the thermal structure of Cordilleran-type mountain belts (i.e. mountain belts related to oceanic subduction underneath a continent}.We integrate apatite and zircon (U-Th)/He data with numerical thermo-kinematic models to study the relative effects of magmatic and surface processes on the thermal evolution of the crust and cooling patterns in the Cenozoic North Cascades arc (Washington State, USA). Two age-elevation profiles that are located 7 km south of the well-studied Chiliiwack intrusions shows that spatial and temporal variability in geothermal gradients linked to magma emplacement can be contrained and separated from exhumation processes. During Chiliiwack batholith emplacement at -35-20 Ma, the geothermal gradient of the country rocks increased to a very high steady-state value (80-100°C/km), which is likely a function of magma flux and the distance from the magma source area. Including temporally varying geothermal gradients in the analysis allows quantifying the thermal perturbation around magmatic intrusions and retrieving a relatively simple denudation history from the data.The synthesis of new and previously published (U-Th)/He data reveals that denudation of the Northern Cascades is spatially and temporally constant at -100 m/Ma between ~32 and ~2 Ma, which likely reflects uplift due to magmatic crustal thickening since the initiation of the Cenozoic stage of the continental magmatic arc. In contrast, the humid flank of the North Cascades is affected by a ten-fold acceleration in exhumation rate at ~2 Ma, which we interpret as forced by the initiation of glaciations; around 3 km of crust have been eroded since that time. Glaciations have three distinct effects on the dynamics of this mountain range: (1) they increase erosion, exhumation and uplift rates where precipitation rates are sufficient to drive efficient glacial erosion; (2) they efficiently limit the elevation of the range; (3) they lead to widening of the humid flank and contraction of the arid flank of the belt.Drainage reorganizations constitute an important agent of landscape evolution that is often underestimated to the benefit of tectonic or climatic events. We propose a new method that integrates geomorphology, low-temperature thermochronometry (apatite and zircon {U-Th-Sm)/He), and 3D numerical thermal-kinematic modelling to detect and date drainage instability producing recent gorge incision, and apply this approach to the Skagit River Gorge, North Cascades.Two zircon (U-Th)/He age-elevation profiles sampled on both sides of the Motagua Fault Zone (MFZ), the boundary between the North American and the Caribbean plates, combined with published thermochronological data show that strike-slip displacement has juxtaposed the cold Maya block (s.s.) against the hot, arc derived, Chortis block (s.s ), producing different age patterns on both sides of the fault and short-wavelength lateral thermal exchange, resulting in recent heating of the cool side and cooling of the hot side of the MFZ.Finally, an apatite (U-Th-Sm)/He age-elevation profile records rapid cooling at -35 Ma localized only in the upper crust along the northern side of the Motagua fault zone. We will try to reproduce these data by modeling the thermal perturbation resulting from the formation of a transtensional basin and of fluid flow activity along a crustal- scale strike-slip fault.

Volatile transport during the crystallization of anatectic melts: Oxygen, boron and hydrogen stable isotope study on the metamorphic complex of Naxos, Greece

Crystallization of anatectic melts in high-temperature metamorphic terrains releases volatile-rich magmas that can be transported into adjacent lithologies. This study addresses the variations in the oxygen, boron and hydrogen isotopic composition of aplite-pegmatite dikes that formed during the crystallization of anatectic melts in regional high-temperature metamorphism on the island of Naxos, Greece, and propagated upward into the overlying sequences of metamorphic schist. The transport distance of these dikes was increased through a significant horizontal component of travel that was imposed by contemporaneous low-angle extensional shearing. Laser fluorination oxygen isotope analyses of quartz, tourmaline, garnet, and biotite mineral separates from the aplite-pegmatite dikes show a progressive rise in delta(18)O values with increasing distance from the core. Oxygen isotope fractionations among quartz, tourmaline, and garnet show temperature variations from > 700degreesC down to similar to400degreesC. This range is considered to reflect isotopic fractionation beginning with crystallization at high temperatures in water-undersaturated conditions and then evolving through lower temperature crystallization and retrograde sub-solidus exchange. Two processes are examined for the cause of the progressive increase in delta(18)O values: (1) heterogeneous delta(18)O sources and (2) fluid-rock exchange between the aplite/pegmatite magmas and their host rock. Although the former process cannot be ruled out, there is as yet no evidence in the exposed sequences on Naxos for the presence of a suitable high delta(18)O magma source. In contrast, a tendency for the delta(18)O of quartz in the aplite/pegmatite dikes to approach that of the quartz in the metamorphic rock suggests that fluid-rock exchange with the host rock may potentially be an important process. Advection of fluid into the magma is examined based on Darcian fluid flow into an initially water-undersaturated buoyantly propagating aplitic dike magma. It is shown that such advective flow could only account for part of the O-18-enrichment, unless it were amplified by repeated injection of magma pulses, fluid recycling, and deformation-assisted post-crystallization exchange. The mechanism is, however, adequate to account for hydrogen isotope equilibration between dike and host rock. In contrast, variations in the delta(11)B values of tourmalines suggest that B-11/B-10 fractionation during crystallization and/or magma degassing was the major control of boron geochemistry rather than fluid-rock interaction and that the boron isotopic system was decoupled from that of oxygen. Copyright (C) 2003 Elsevier Ltd.

Seismic attenuation in heterogeneous partially saturated rocks

We study wave-induced fluid flow effects in porous rocks partially saturated with gas and water, where the saturation patterns are governed by mesoscopic heterogeneities associated with the dry frame properties. The link between the dry frame properties and the gas saturation is defined by the assumption of capillary pressure equilibrium, which in the presence of heterogeneity implies that neighboring regions can exhibit different levels of saturation. In order to determine the equivalent attenuation and phase velocity of the synthetic rock samples considered in this study, we apply a numerical upscaling procedure, which permits to take into account mesoscopic heterogeneities associated with the dry frame properties as well as spatially continuous variations of the pore fluid properties. We consider numerical experiments to analyze such effects in heterogeneous partially saturated porous media, where the saturation field is determined by realistic variations in porosity. Our results indicate that the spatially continuous nature of gas saturation inherent to this study is a critical parameter controlling the seismic response of these environments, which in turn suggests that the physical mechanisms governing partial saturation should be accounted for when analyzing seismic data in a poro-elastic context.

Suprapineal recess: an alternate site for third ventriculostomy? Case report.

This 30-year-old woman presented with clinical symptoms and signs of intracranial hypertension and Parinaud syndrome secondary to ventriculoperitoneal shunt dysfunction. Magnetic resonance (MR) imaging revealed gross triventricular hydrocephalus with a large suprapineal recess due to aqueductal stenosis. Using an endoscopic approach, a ventriculostomy was performed within the floor of the dilated suprapineal recess. Following this procedure the patient experienced alleviation of all her neurological symptoms and signs. Postoperative MR imaging and cerebrospinal fluid flow studies demonstrated a functioning ventriculostomy. The anatomy of the suprapineal recess and its suitability for endoscopic ventriculostomy are discussed.

Geomodelling of a fluvial system with semi-supervised support vector regression

Fluvial deposits are a challenge for modelling flow in sub-surface reservoirs. Connectivity and continuity of permeable bodies have a major impact on fluid flow in porous media. Contemporary object-based and multipoint statistics methods face a problem of robust representation of connected structures. An alternative approach to model petrophysical properties is based on machine learning algorithm ? Support Vector Regression (SVR). Semi-supervised SVR is able to establish spatial connectivity taking into account the prior knowledge on natural similarities. SVR as a learning algorithm is robust to noise and captures dependencies from all available data. Semi-supervised SVR applied to a synthetic fluvial reservoir demonstrated robust results, which are well matched to the flow performance

The sandstone-hosted Beverley uranium Deposit, Lake Frome Basin, South Australia: Mineralogy, geochemistry, and a time-constrained model for its genesis

The sandstone-hosted Beverley uranium deposit is located in terrestrial sediments in the Lake Frome basin in the North Flinders Ranges, South Australia. The deposit is 13 km from the U-rich Mesoproterozoic basement of the Mount Painter inlier, which is being uplifted 100 to 200 m above the basin by neotectonic activity that probably initiated in the early Pliocene. The mineralization was deposited mainly in organic matter-poor Miocene lacustrine sands and partly in the underlying reductive strata comprising organic matter-rich clays and silts. The bulk of the mineralization consists of coffinite and/or uraninite nodules, growing around Co-rich pyrite with an S isotope composition (delta S-34 = 1.0 +/- 0.3 parts per thousand), suggestive of an early diagenetic lacustrine origin. In contrast, authigenic sulfides in the bulk of the sediments have a negative S isotope signature (delta S-34 ranges from -26.2 to -35.5 parts per thousand), indicative of an origin via bacterially mediated sulfate reduction. Minor amounts of Zn-bearing native copper and native lead also support the presence of specific, reducing microenvironments in the ore zone. Small amounts of carnotite are associated with the coffinite ore and also occur beneath a paleosoil horizon overlying the uranium deposit. Provenance studies suggest that the host Miocene sediments were derived from the reworking of Early Cretaceous glacial or glaciolacustrine sediments ultimately derived from Paleozoic terranes in eastern Australia. In contrast, the overlying Pliocene strata were in part derived from the Mesoproterozoic basement inlier. Mass-balance and geochemical data confirm that granites of the Mount Painter domain were the ultimate source of U and BEE at Beverley. U-Pb dating of coffinite and carnotite suggest that the U mineralization is Pliocene (6.7-3.4 Ma). The suitability of the Beverley deposit for efficient mining via in situ leaching, and hence its economic value, are determined by the nature of the hosting sand unit, which provides the permeability and low reactivity required for high fluid flow and low chemical consumption. These favorable sedimentologic and geometrical features result from a complex conjunction of factors, including deposition in lacustrine shore environment, reworking of angular sands of glacial origin, deep Pliocene weathering, and proximity to an active fault exposing extremely U rich rocks.

Emplacement of ultramafic rocks into the continental crust monitored by light and other trace elements: An example from the Geisspfad body (Swiss-Italian Alps)

In order to evaluate the influence of continental crustal rocks on trace element budgets of serpentinized peridotites incorporated into the continental crust, we have analyzed the chemical composition of whole rock samples and minerals of the Geisspfad ultramafic complex (Swiss-Italian Alps). This complex represents a relict oceanic succession composed of serpentinites, ophicarbonates and metabasic rocks, emplaced into crustal gneisses during Alpine collision. Following peak metamorphic amphibolite facies conditions, fluid flow modified some of the trace element contents of ophicarbonates and deformed serpentinites close to the contact with country rocks. The fluid originated from the surrounding continental crustal rocks as documented by the increase of Pb in the serpentinites, and by the strongly negative all) values (-112 parts per thousand) of some ultramafic rocks close to the contact with surrounding gneisses. Little or no modification of the fluid mobile elements Li, B or U was observed in the serpentinite. In-situ analysis of light elements of serpentinite minerals indicate redistribution of light elements coupled to changes of mineral modes towards the outer 100-150 m of the massif. In the centre of the massif, Li is preferentially concentrated in olivine, while Be and B are hosted by tremolite. In contrast, at the outer rim of the massif, Li and Be are preferentially incorporated into diopside, and B into antigorite. This redistribution of light elements among the different minerals is visible in the serpentinite, at a maximum distance of -100-150 m from the ophicarbonate-metabasite contact. Our results show that interaction of ultramafic rocks and crust-derived fluids can be easily detected by studies of Pb and partial derivative D in whole rocks. We argue that small ultramafic bodies potentially record an emplacement-related trace element signature, and that crustal light element values in ultramafic rocks are not necessarily derived from a subducting slab. (C) 2008 Elsevier B.V. All rights reserved.