Geochemistry and 40Ar/39Ar geochronology of pseudotachylyte associated with UHP whiteschists from the Dora Maira massif, Italy

In situ UV-Iaser ablation Ar-40/(39) Ar geochronological and geochemical data, together with rock and mineral compositional data, have been determined from pseudotachylyte and surrounding mylonitic gneiss associated with the UHP whiteschists of the Dora Maira Massif, Italy. Several generations of fresh pseudotachylyte occur as irregular veins up to a few cur thick both parallel and at high angles to the foliation. Whole rock XRF data collected from representative lithologies of mylonitic gneiss are uniformly consistent with a mildly alkalic granitic protolith. Minimal compositional variation is observed between the pseudotachylyte and its surrounding mylonitic gneiss. The pseudotachylyte contains newly crystallized grains of biotite and K-feldspar in a matrix of glass with partially fused grains of quartz, zircon, apatite, and titanite. Electron microprobe analyses of the glass show significant compositional variation that is probably strongly influenced by micrometer-scale changes in mineralogy. UV-Iaser ablation ICP-MS traverses across the mylonitic gneiss-pseudotachylyte contact are consistent with cataclastic communition of REE carriers such as epidote, monazite, allanite, zircon, and apatite before melting as an efficient mechanism of REE homogenization in the pseudotachylyte. The 40Ar/39Ar data from one band of pseudotachylyte indicate formation at 20.1 +/- 0.5 Ma, when the mylonitic gneisses were already in a near surface position. The variable effects of top-to-the-west shear deformation within outcrops of the coesite-bearing unit are reflected in localized zones of protomylonite, cataclasite, ultracataclasite, and pseudotachylyte. Preservation of several generations of pseudotachylyte suggests that seismic events may have played a significant role in triggering late unroofing of the UHP rocks. It is speculated that deeper crustal seismic events potentially played a role in the unroofing of the UHP rocks at earlier stages in their exhumation history. (c) 2005 Elsevier B.V. All rights reserved.

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.

Large calcite and bulk-rock volume loss in metacarbonate xenoliths from the Qu,rigut massif (French Pyrenees)

Chemical mass transfer was quantified in a metacarbonate xenolith enclosed within the granodiorite of the Qu,rigut massif (Pyrenees, France). Mass balance calculations suggest a strong decrease of CaO, SrO and CO(2) contents (up to -90%), correlated with a decrease of modal calcite content as the contact is approached. Most other chemical elements behave immobile during metasomatism. They are therefore passively enriched. Only a small increase of SiO(2), Al(2)O(3) and Fe(2)O(3) contents occurs in the immediate vicinity of the contact. Hence, in this study, skarn formation is characterized by the lack of large chemical element influx from the granitoid protolith. A large decrease of the initial carbonate volume (up to -86%) resulted from a combination of decarbonation reactions and loss of CaO and CO(2). The resulting volume change has potentially important consequences for the interpretation of stable isotope profiles: the isotope alteration could have occured over greater distances than those observed today.

Metallogenic model of the Trepca Pb-Zn-Ag skarn deposit, Kosovo: Evidence from fluid inclusions, Rare Earth Elements, and stable isotope data

The Trepca Pb-Zn-Ag skarn deposit (29 Mt of ore at 3.45% Pb, 2.30% Zn, and 80 g/t Ag) is located in the Kopaonik block of the western Vardar zone, Kosovo. The mineralization, hosted by recrystallized limestone of Upper Triassic age, was structurally and lithologically controlled. Ore deposition is spatially and temporally related with the postcollisional magmatism of Oligocene age (23-26 Ma). The deposit was formed during two distinct mineralization stages: an early prograde closed-system and a later retrograde open-system stage. The prograde mineralization consisting mainly of pyroxenes (Hd(54-100)Jo(0-45)Di(0-45)) resulted from the interaction of magmatic fluids associated with Oligocene (23-26 Ma) postcollisional magmatism. Whereas there is no direct contact between magmatic rocks and the mineralization, the deposit is classified as a distal Pb-Zn-Ag skarn. Abundant pyroxene reflects low oxygen fugacity (<10(-31) bar) and anhydrous environment. Fluid inclusion data and mineral assemblage limit the prograde stage within a temperature range between 390 degrees and 475 degrees C. Formation pressure is estimated below 900 bars. Isotopic composition of aqueous fluid, inclusions hosted by hedenbergite (delta D = -108 to -130 parts per thousand; delta O-18 = 7.5-8.0 parts per thousand), Mn-enriched mineralogy and high REE content of the host carbonates at the contact with the skarn mineralization suggest that a magmatic fluid was modified during its infiltration through the country rocks. The retrograde mineral assemblage comprises ilvaite, magnetite, arsenopyrite, pyrrhotite, marcasite, pyrite, quartz, and various carbonates. Increases in oxygen and sulfur fugacities, as well as a hydrous character of mineralization, require an open-system model. The opening of the system is related to phreatomagmatic explosion and formation of the breccia. Arsenopyrite geothermometer limits the retrograde stage within the temperature range between 350 degrees and 380 degrees C and sulfur fugacity between 10(-8.8) and 10(-7.2) bars. The principal ore minerals, galena, sphalerite, pyrite, and minor chalcopyrite, were deposited from a moderately saline Ca-Na chloride fluid at around 350 degrees C. According to the isotopic composition of fluid inclusions hosted by sphalerite (delta D = -55 to -74 parts per thousand; delta O-18 = -9.6 to -13.6 parts per thousand), the fluid responsible for ore deposition was dominantly meteoric in origin. The delta S-31 values of the sulfides spanning between -5.5 and +10 parts per thousand point to a magmatic origin of sulfur. Ore deposition appears to have been largely contemporaneous with the retrograde stage of the skarn development. Postore stage accompanied the precipitation of significant amount of carbonates including the travertine deposits at the deposit surface. Mineralogical composition of travertine varies from calcite to siderite and all carbonates contain significant amounts of Mn. Decreased formation temperature and depletion in the REE content point to an influence of pH-neutralized cold ground water and dying magmatic system.

Middle and late Cenomanian oceanic anoxic events in shallow and deeper shelf environments of western Morocco

The response of shallow-water sequences to oceanic anoxic event 2 and mid-Cenomanian events 1a and 1b was investigated along the west African margin of Morocco north of Agadir (Azazoul) and correlated with the deep-water sequence of the Tarfaya Basin (Mohammed Beach) based on biostratigraphy, mineralogy, phosphorus and stable isotopes. In the deeper Mohammed Beach section results show double peaks in delta 13C(org) for mid-Cenomanian events 1a and 1b (Rotalipora reicheli biozone, lower CC10a biozone), the characteristic oceanic anoxic event 2 delta 13C excursion (Rotalipora cushmani extinction, top of CC10a biozone) and laminated (anoxic) black shale. In the shallow environment north of Agadir, a fluctuating sea-level associated with dysoxic, brackish and mesotrophic conditions prevailed during the middle to late Cenomanian, as indicated by oyster biostromes, nannofossils, planktonic and benthonic foraminiferal assemblages. Anoxic conditions characteristic of oceanic anoxic event 2 (for example, laminated black shales) did not reach into shallow-water environments until the maximum transgression of the early Turonian. Climate conditions decoupled along the western margin of Morocco between mid-Cenomanian event 1b and the Cenomanian-Turonian boundary, as also observed in eastern Tethys. North of Agadir alternating humid and dry seasonal conditions prevailed, whereas in the Tarfaya Basin the climate was dry and seasonal. This climatic decoupling can be attributed to variations in the Intertropical Convergence Zone and in the intensity of the north-east trade winds in tropical areas.

Alpine and pre-Alpine magmatism in the root-zone of the western Central Alps

The highest grade of metamorphism and associated structural elements in orogenic belts may be inherited from earlier orogenic events. We illustrate this point using magmatic and metamorphic rocks from the southern steep belt of the Lepontine Gneiss Dome (Central Alps). The U-Pb zircon ages from an anatectic granite at Verampio and migmatites at Corcapolo and Lavertezzo yield 280-290 Ma, i.e., Hercynian ages. These ages indicate that the highest grade of metamorphism in several crystalline nappes of the Lepontine Gneiss Dome is pre-Alpine. Alpine metamorphism reached sufficiently high grade to reset the Rb-Sr and K-Ar systematics of mica and amphibole, but generally did not result in crustal melting, except in the steep belt to the north of the Insubric Line, where numerous 29 to 26 Ma old pegmatites and aplites had intruded syn- and post-kinematically into gneisses of the ductile Simplon Shear Zone. The emplacement age of these pegmatites gives a minimum estimate for the age of the Alpine metamorphic peak in the Monte Rosa nappe. The U-Pb titanite ages of 33 to 31 Ma from felsic porphyritic veins represent a minimum-age estimate for Alpine metamorphism in the Sesia Zone. A porphyric vein emplaced at 448 +/- 5 Ma (U-Pb monazite) demonstrates that there existed a consolidated Caledonian basement in the Sesia Zone.

Properties of debris flow deposits and source materials compared: implications for debris flow characterization

In the Alps, debris flow deposits generally contain < 5% clay-size particles, and the role of the surface-charged < 2 mu m particles is often neglected, although these particles may have a significant impact on the rheological properties of the interstitial fluid. The objective of this study was to compare debris flow deposits and parent materials from two neighbouring catchments of the Swiss Alps, with special emphasis on the colloidal constituents. The catchments are small in area (4 km(2)), 2.5 km long, similar in morphology, but different in geology. The average slopes are 35-40%. The catchments were monitored for debris flow events and mapped for surface aspect and erosion activity. Debris flow deposits and parent materials were sampled, the clay and silt fractions extracted and the bulk density, < 2 mm fraction bulk density, particle size distribution, chemical composition, cation exchange capacity (CEC) and mineralogy analysed. The results show that the deposits are similar to the parent screes in terms of chemical composition, but differ in terms of: (i) particle size distribution; and (ii) mineralogy, reactivity and density of the < 2 mm fraction. In this fraction, compared with the parent materials the deposits show dense materials enriched in coarse monocrystalline particles, of which the smallest and more reactive particles were leached. The results suggest that deposit samples should not be considered as representative of source or flow materials, particularly with respect to their physical properties.

Alkali-calcic and alkaline post-orogenic (PO) granite magmatism: petrologic constraints and geodynamic settings

The end of an orogenic Wilson cycle corresponds to amalgamation of terranes into a Pangaea and is marked by widespread magmatism dominated by granitoids. The post-collision event starts with magmatic processes still influenced by subducted crustal materials. The dominantly calc-alkaline suites show a shift from normal to high-K to very high-K associations. Source regions are composed of depleted and later enriched orogenic subcontinental lithospheric mantle, affected by dehydration melting and generating more and more K- and LILE-rich magmas. In the vicinity of intra-crustal magma chambers, anatexis by incongruent melting of hydrous minerals may generate peraluminous granitoids bearing mafic enclaves. The post-collision event ends with emplacement of bimodal post-orogenic (PO) suites along transcurrent fault zones. Two suites are defined, (i) the alkali-calcic monzonite-monzogranite-syenogranite-alkali feldspar granite association characterised by [biotite + plagioclase] fractionation and moderate [LILE + HFSE] enrichments and (ii) the alkaline monzonite-syenite-alkali feldspar granite association characterised by [amphibole + alkali feldspar] fractionation and displaying two evolutionary trends, one peralkaline with sodic mafic mineralogy and higher enrichments in HFSE than in LILE, and the other aluminous biotite-bearing marked by HFSE depletion relative to LILE due to accessory mineral precipitation. Alkali-calcic and alkaline suites differ essentially in the amounts of water present within intra-crustal magma chambers, promoting crystallisation of various mineral assemblages. The ultimate enriched and not depleted mantle source is identical for the two PO suites. The more primitive LILE and HFSE-rich source rapidly replaces the older orogenic mantle source during lithosphere delamination and becomes progressively the thermal boundary layer of the new lithosphere. Present rock compositions are a mixture of major mantle contribution and various crustal components carried by F-rich aqueous fluids circulating within convective cells created around magma chambers. In favourable areas, PO suites pre-date a new orogenic Wilson cycle. (C) 1998 Elsevier Science B.V. All rights reserved.

The roles of flux- and decompression melting and their respective fractionation lines for continental crust formation: Evidence from the Kohistan arc

Delamination and foundering of the lower continental crust (LCC) into the mantle is part of the crust-forming mechanism. However, knowledge of the composition and mineralogy of the preserved or delaminated LCC over geological timescales remains scarce. We provide a synopsis of recent research within the Kohistan arc (Pakistan) and demonstrate that hydrous and less hydrous liquid lines of descent related to flux assisted and decompression mantle melting, respectively, produce compositionally different lower crustal rocks. The argument refers to two lower crustal sections exposed in Kohistan, the older Southern Plutonic Complex (SPC) and the younger Chilas Complex. The SPC typifies a hydrous, high-pressure fractionation sequence of olivine-pyroxenes-garnet-Fe/Ti-oxide-amphibole-plagioclase. The Chilas Complex illustrates a less hydrous fractionation sequence of olivine-clinopyroxene-orthopyroxene-plagioclase-amphibole. Despite the similarity of the Chilas Complex rocks to proposed lower crust compositions, the less hydrous fractionation results in unrealistically small volumes of silica-rich rocks, precluding the Chilas Complex gabbros to represent the magmatic complement to the upper crust. The composition of the SPC lower crust differs markedly from bulk lower crust estimates, but is complementary to silica-rich rocks exposed along this section and in the Kohistan batholith. These observations inspire a composite model for the formation of continental crust (CC) where the negatively buoyant delaminated and the buoyant preserved lower continental crusts (LCC) differ in genesis, mineralogy, and composition. We propose that the upper, non-sedimentary subsequent removal of the complementary, negatively buoyant garnet-pyroxene-amphibole-plagioclase-rich cumulates. In contrast, the LCC, which is buoyant and preserved over geological timescales, is formed by less hydrous parental mantle melts. We suggest that the bulk continental crust composition is related to mixing of these petrologically not directly related end members. Published by Elsevier B.V.

Geochemical and H-O-Sr-Nd isotope evidence for magmatic processes and meteoric-water interactions in the basal complex of La Gomera, Canary Islands

The plutonic rocks of the Basal Complex of La Gomera, Canary Islands, Spain, were studied by means of major and trace element contents and by H-O-Sr-Nd isotope compositions in order to distinguish primary magmatic characteristics and late-stage alteration products. Deciphering the effects of alteration allowed us to determine primary, plume-related compositions that indicated D- and (18)O-depletion relative to normal upper mantle, supporting the conclusions of earlier studies on the plutonic rocks of Fuerteventura and La Palma. Late-stage alteration took place during the formation of the intrusive series induced by interaction with meteoric water. Inferred isotopic compositions of the meteoric water indicate that the water infiltrated into the rock edifice at a height of about 1500 m above sea level, suggesting the existence of a subaerial volcano which was active during the intrusive activity and that it has been either distroyed or remain buried by later volcanic and landslide events.

Origin of siderite mineralisation in Petrova and Trgovska Gora Mts., NW Dinarides

The Petrova and Trgovska Gora Mts. (Gora=Mountain) are Variscan basement units incorporated into the northwestern Dinarides during the Alpine orogeny. They host numerous siderite-quartz-polysulphide, siderite-chalcopyrite, siderite-galena and barite veins, as well as stratabound hydrothermal-replacement ankerite bodies within carbonates in non-metamorphosed, flysch-like Permo-Carboniferous sequences. The deposits have been mined for Cu, Pb, Ag and Fe ores since Medieval times. Fluid inclusion studies of quartz from siderite-polysulphide-quartz and barite veins of both regions have shown the presence of primary aqueous NaCl-CaCl(2)+/- MgCl(2)-H(2)O +/- CO(2) inclusions. The quartz-sulphide stage of both regions show variable salinities; 2.7-26.2 wt% NaCl eq. for the Trgovska Gora region and 3.4-23.4 wt% NaCl eq. for the Petrova gora region, and similar homogenisation temperatures (100-230A degrees C). Finally, barite is precipitated from low salinity-low temperature solutions (3.7-15.8 wt % NaCl equ. and 115-145A degrees C). P-t conditions estimated via isochore construction yield formation temperatures between 180-250A degrees C for the quartz-sulphide stage and 160-180A degrees C for the barite stage, using a maximum lithostatic pressure of 1 kbar (cc. 3 km of overburden). The sulphur isotope composition of barite from both deposits indicates the involvement of Permian seawater in ore fluids. This is supported by the elevated bromium content of the fluid inclusion leachates (120-660 ppm in quartz, 420-960 ppm in barite) with respect to the seawater, indicating evaporated seawater as the major portion of the ore-forming fluids. Variable sulphur isotope compositions of galena, pyrite and chalcopyrite, between -3.2 and +2.7aEuro degrees, are interpreted as a product of incomplete thermal reduction of the Permian marine sulphate mixed with organically- and pyrite-bound sulphur from the host sedimentary rocks. Ore-forming fluids are interpreted as deep-circulating fluids derived primarily from evaporated Permian seawater and later modified by interaction with the Variscan basement rocks. (40)Ar/(39)Ar data of the detrital mica from the host rocks yielded the Variscan age overprinted by an Early Permian tectonothermal event dated at 266-274 Ma. These ages are interpreted as those reflecting hydrothermal activity correlated with an incipient intracontinental rifting in the Tethyan domain. Nevertheless, 75 Ma recorded at a fine-grained sericite sample from the alteration zone is interpreted as a result of later resetting of white mica during Campanian opening/closure of the Sava back arc in the neighbouring Sava suture zone (Ustaszewski et al. 2008).

The 500 Ma-old Thyon metagranite: A new A-type granite occurrence in the western Penninic Alps (Wallis, Switzerland)

The Thyon metagranite is located in the frontal part of the Siviez-Mischabel Nappe, in the western Penninic Alps. It is intrusive in a polymetamorphic banded volcanic complex as leucocratic concordant sills with pseudoaplitic rims. A distinct metamorphic schistosity is defined by dark-green Fe-rich biotite. Abundant mesoperthites, chess-board albite and low microcline are presumably related to magmatic stages and/or greenschist-facies metamorphic retrogression. Major, trace element and REE geochemistry, zircon typology, Y and Nb-bearing accessory minerals such as fergusonite and euxenite, all point to a metaluminous to peraluminous alkaline A-type granite. High-precision U-Pb zircon dating yielded a sub-concordant age of 500 +3/-4 Ma. The Thyon metagranite is the third record of a Cambro-Ordovician alkaline magmatic activity in the Alps. As A-type granitic magmatism is common in post-orogenic to anorogenic extensional tectonic regime, the Thyon intrusion could mark the transition between the Cadomian and the Caledonian orogenies.

Melt variability in percolated peridotite: an experimental study applied to reactive migration of tholeiitic basalt in the upper mantle

Melt-rock reaction in the upper mantle is recorded in a variety of ultramafic rocks and is an important process in modifying melt composition on its way from the source region towards the surface. This experimental study evaluates the compositional variability of tholeiitic basalts upon reaction with depleted peridotite at uppermost-mantle conditions. Infiltration-reaction processes are simulated by employing a three-layered set-up: primitive basaltic powder ('melt layer') is overlain by a 'peridotite layer' and a layer of vitreous carbon spheres ('melt trap'). Melt from the melt layer is forced to move through the peridotite layer into the melt trap. Experiments were conducted at 0.65 and 0.8 GPa in the temperature range 1,170-1,290 degrees C. In this P-T range, representing conditions encountered in the transition zone (thermal boundary layer) between the asthenosphere and the lithosphere underneath oceanic spreading centres, the melt is subjected to fractionation, and the peridotite is partially melting (T (s) similar to 1,260 degrees C). The effect of reaction between melt and peridotite on the melt composition was investigated across each experimental charge. Quenched melts in the peridotite layers display larger compositional variations than melt layer glasses. A difference between glasses in the melt and peridotite layer becomes more important at decreasing temperature through a combination of enrichment in incompatible elements in the melt layer and less efficient diffusive equilibration in the melt phase. At 1,290A degrees C, preferential dissolution of pyroxenes enriches the melt in silica and dilutes it in incompatible elements. Moreover, liquids become increasingly enriched in Cr(2)O(3) at higher temperatures due to the dissolution of spinel. Silica contents of liquids decrease at 1,260 degrees C, whereas incompatible elements start to concentrate in the melt due to increasing levels of crystallization. At the lowest temperatures investigated, increasing alkali contents cause silica to increase as a consequence of reactive fractionation. Pervasive percolation of tholeiitic basalt through an upper-mantle thermal boundary layer can thus impose a high-Si 'low-pressure' signature on MORB. This could explain opx + plag enrichment in shallow plagioclase peridotites and prolonged formation of olivine gabbros.

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.