Deciphering polymetamorphism by garnet growth modelling : an application to the Austroalpine crystalline basement east of the Tauern Window

This thesis is based upon studies conducted between June 2003 and June 2007 at the Institute of Mineralogy and Petrography at Basel University. It was developed under the supervision of Rainer Abart (FU Berlin) and Christian de Capitani (Basel University) and essentially benefited from the guidance of Ralf Schuster (Austrian Geological Survey) in terms of the regional geological framework. The thesis is divided into six chapters, four of which (Chapter 2 to 5) can be regarded as individual manuscripts: Chapter 1 (‘Introduction’) briefly presents the general background and the aim of the thesis. Chapter 2 (‘Characterization of polymetamorphism in the Austroalpine basement east of the Tauern Window using garnet isopleth thermobarometry’) focuses on the description of the method ‘garnet isopleth thermobarometry’ and its application to polymetamorphic garnet porphyroblasts from the Austroalpine basement. This method allows to estimate P-T conditions for the formation of pre-existing equilibrium assemblages, employing bulkrock chemistries and garnet as the only relict mineral phase. This study is published in the Journal of Metamorphic Geology. Chapter 3 (‘THERIA G: A software program to numerically model prograde garnet growth’) contains the description of the software program THERIA G, which was developed throughout the course of the preparation of this thesis. THERIA G extends earlier work to simulate the formation of garnet in a given volume of rock to garnet populations. Chapter 3 clearly demonstrates that the chemical composition of garnet porphyroblasts may be influenced by chemical fractionation between garnet and rock matrix during growth, intracrystalline diffusion in garnet, and the nucleation history. This study is under review in Contributions to Mineralogy and Petrology. Chapter 4 [‘Prograde garnet growth along complex P-T-t paths: Results from numerical experiments on polyphase garnet from the W¨olz Complex (Austroalpine basement)’] presents an application of THERIA G. Due to the application of THERIA G, essential parts of the P-T-t paths of two different metamorphic events can be obtained from a single polyphase garnet porphyroblast of the Austroalpine basement. This study is under review in Contributions to Mineralogy and Petrology. Chapter 5 [‘Coupling garnet growth forward modelling with monazite geochronology: An application to the Rappold Complex (Austroalpine crystalline basement)’] links thermobarometrical information obtained from garnet growth modelling with microprobe monazite ages. As a result, a possible geothermal history is suggested, which indicates the complex metamorphic evolution of the Austroalpine crystalline basement. This study is to be submitted for publication in the Journal of Metamorphic Geology. Chapter 6 (‘Overall conclusion’) presents the main summery and conclusion of the entire thesis. The Appendix, which contains supplementary material to all the chapters, and the enclosed CD, at which the remaining data are stored, complete the thesis

THERIA_G: a software program to numerically model prograde garnet growth

We present the software program THERIA_G, which allows for numerical simulation of garnet growth in a given volume of rock along any pressure-temperature-time (P-T-t) path. THERIA_G assumes thermodynamic equilibrium between the garnet rim and the rock matrix during growth and accounts for component fractionation associated with garnet formation as well as for intracrystalline diffusion within garnet. In addition, THERIA_G keeps track of changes in the equilibrium phase relations, which occur during garnet growth along the specified P-T-t trajectory. This is accomplished by the combination of two major modules: a Gibbs free energy minimization routine is used to calculate equilibrium phase relations including the volume and composition of successive garnet growth increments as P and T and the effective bulk rock composition change. With the second module intragranular multi-component diffusion is modelled for spherical garnet geometry. THERIA_G allows to simulate the formation of an entire garnet population, the nucleation and growth history of which is specified via the garnet crystal size frequency distribution. Garnet growth simulations with THERIA_G produce compositional profiles for the garnet porphyroblasts of each size class of a population and full information on equilibrium phase assemblages for any point along the specified P-T-t trajectory. The results of garnet growth simulation can be used to infer the P-T-t path of metamorphism from the chemical zoning of garnet porphyroblasts. With a hypothetical example of garnet growth in a pelitic rock we demonstrate that it is essential for the interpretation of the chemical zoning of garnet to account for the combined effects of the thermodynamic conditions of garnet growth, the nucleation history and intracrystalline diffusio

Prograde garnet growth along complex P-T-t paths: results from numerical experiments on polyphase garnet from the Wölz Complex (Austroalpine basement)

Garnet in metapelites from the Wölz Complex of the Austroalpine crystalline basement east of the Tauern Window characteristically consists of two growth phases, which preserve a comprehensive record of the geothermal history during polymetamorphism. From numerical modelling of garnet formation, detailed information on the pressure-temperature-time (P-T-t) evolution during prograde metamorphism is obtained. In that respect, the combined influences of chemical fractionation associated with garnet growth, modification of the original growth zoning through intragranular diffusion and the nucleation history on the chemical zoning of garnet as P and T change during growth are considered. The concentric chemical zoning observed in garnet and the homogenous rock matrix, which is devoid of chemical segregation, render the simulation of garnet growth through successive equilibrium states reliable. Whereas the first growth phase of garnet was formed at isobaric conditions of ∼3.8kbar at low heating/cooling rates, the second growth phase grew along a Barrovian P-T path marked with a thermal peak of ∼625°C at ∼10kbar and a maximum in P of ∼10.4kbar at ∼610°C. For the heating rate during the growth of the second phase of garnet, average rates faster than 50°CMa−1 are obtained. From geochronological investigations the first growth phase of garnet from the Wölz Complex pertains to the Permian metamorphic event. The second growth phase grew in the course of Eo-Alpine metamorphism during the Cretaceou

Integration of natural data within a numerical model of ablative subduction: A possible interpretation for the Alpine dynamics of the Austroalpine crust

A numerical modelling approach is used to validate the physical and ge- ological reliability of the ablative subduction mechanism during Alpine con- vergence in order to interpret the tectonic and metamorphic evolution of an inner portion of the Alpine belt: the Austroalpine Domain. The model pre- dictions and the natural data for the Austroalpine of the Western Alps agree very well in terms of P-T peak conditions, relative chronology of peak and exhumation events, P-T-t paths, thermal gradients and the tectonic evolu- tion of the continental rocks. These findings suggest that a pre-collisional evolution of this domain, with the burial of the continental rocks (induced by ablative subduction of the overriding Adria plate) and their exhumation (driven by an upwelling flow generated in a hydrated mantle wedge) could be a valid mechanism that reproduces the actual tectono-metamorphic config- uration of this part of the Alps. There is less agreement between the model predictions and the natural data for the Austroalpine of the Central-Eastern Alps. Based on the natural data available in the literature, a critical discus- sion of the other proposed mechanisms is presented, and additional geological factors that should be considered within the numerical model are suggested to improve the fitting to the numerical results; these factors include varia- tions in the continental and/or oceanic thickness, variation of the subduction rate and/or slab dip, the initial thermal state of the passive margin, the oc- currence of continental collision and an oblique convergence.Comment: 11 Figures and 3 Tabe

Argon behaviour in an inverted Barrovian sequence, Sikkim Himalaya: the consequences of temperature and timescale on ⁴⁰Ar/³⁹Ar mica geochronology

40Ar/39Ar dating of metamorphic rocks sometimes yields complicated datasets which are difficult to interpret in terms of timescales of the metamorphic cycle. Single-grain fusion and step-heating data were obtained for rocks sampled through a major thrust-sense shear zone (the Main Central Thrust) and the associated inverted metamorphic zone in the Sikkim region of the eastern Himalaya. This transect provides a natural laboratory to explore factors influencing apparent 40Ar/39Ar ages in similar lithologies at a variety of metamorphic pressure and temperature (P–T) conditions. The 40Ar/39Ar dataset records progressively younger apparent age populations and a decrease in within-sample dispersion with increasing temperature through the sequence. The white mica populations span ~ 2–9 Ma within each sample in the structurally lower levels (garnet grade) but only ~ 0–3 Ma at structurally higher levels (kyanite-sillimanite grade). Mean white mica single-grain fusion population ages vary from 16.2 ± 3.9 Ma (2σ) to 13.2 ± 1.3 Ma (2σ) from lowest to highest levels. White mica step-heating data from the same samples yields plateau ages from 14.27 ± 0.13 Ma to 12.96 ± 0.05 Ma. Biotite yield older apparent age populations with mean single-grain fusion dates varying from 74.7 ± 11.8 Ma (2σ) at the lowest structural levels to 18.6 ± 4.7 Ma (2σ) at the highest structural levels; the step-heating plateaux are commonly disturbed. Temperatures > 600 °C at pressures of 0.4–0.8 GPa sustained over > 5 Ma, appear to be required for white mica and biotite ages to be consistent with diffusive, open-system cooling. At lower temperatures, and/or over shorter metamorphic timescales, more 40Ar is retained than results from simple diffusion models suggest. Diffusion modelling of Ar in white mica from the highest structural levels suggests that the high-temperature rocks cooled at a rate of ~ 50–80 °C Ma− 1, consistent with rapid thrusting, extrusion and exhumation along the Main Central Thrust during the mid-Miocene

Pervasive Eclogitization Due to Brittle Deformation and Rehydration of Subducted Basement: Effects on Continental Recycling?

The buoyancy of continental crust opposes its subduction to mantle depths, except where mineral reactions substantially increase rock density. Sluggish kinetics limit such densification, especially in dry rocks, unless deformation and hydrous fluids intervene. Here we document how hydrous fluids in the subduction channel invaded lower crustal granulites at 50–60 km depth through a dense network of probably seismically induced fractures. We combine analyses of textures and mineral composition with thermodynamic modeling to reconstruct repeated stages of interaction, with pulses of high‐pressure (HP) fluid at 650–670°C, rehydrating the initially dry rocks to micaschists. SIMS oxygen isotopic data of quartz indicate fluids of crustal composition. HP growth rims in allanite and zircon show uniform U‐Th‐Pb ages of ∼65 Ma and indicate that hydration occurred during subduction, at eclogite facies conditions. Based on this case study in the Sesia Zone (Western Italian Alps), we conclude that continental crust, and in particular deep basement fragments, during subduction can behave as substantial fluid sinks, not sources. Density modeling indicates a bifurcation in continental recycling: Chiefly mafic crust, once it is eclogitized to >60%, are prone to end up in a subduction graveyard, such as is tomographically evident beneath the Alps at ∼550 km depth. By contrast, dominantly felsic HP fragments and mafic granulites remain positively buoyant and tend be incorporated into an orogen and be exhumed with it. Felsic and intermediate lithotypes remain positively buoyant even where deformation and fluid percolation allowed them to equilibrate at HP

Developing an inverted Barrovian sequence; insights from monazite petrochronology

In the Himalayan region of Sikkim, the well-developed inverted metamorphic sequence of the Main Central Thrust (MCT) zone is folded, thus exposing several transects through the structure that reached similar metamorphic grades at different times. In-situ LA-ICP-MS U–Th–Pb monazite ages, linked to pressure–temperature conditions via trace-element reaction fingerprints, allow key aspects of the evolution of the thrust zone to be understood for the first time. The ages show that peak metamorphic conditions were reached earliest in the structurally highest part of the inverted metamorphic sequence, in the Greater Himalayan Sequence (GHS) in the hanging wall of the MCT. Monazite in this unit grew over a prolonged period between ~37 and 16 Ma in the southerly leading-edge of the thrust zone and between ~37 and 14.5 Ma in the northern rear-edge of the thrust zone, at peak metamorphic conditions of ~790 ◦C and 10 kbar. Monazite ages in Lesser Himalayan Sequence (LHS) footwall rocks show that identical metamorphic conditions were reached ~4–6 Ma apart along the ~60 km separating samples along the MCT transport direction. Upper LHS footwall rocks reached peak metamorphic conditions of ~655 ◦C and 9 kbar between ~21 and 16 Ma in the more southerly-exposed transect and ~14.5–12 Ma in the northern transect. Similarly, lower LHS footwall rocks reached peak metamorphic conditions of ~580 ◦C and 8.5 kbar at ~16 Ma in the south, and 9–10 Ma in the north. In the southern transect, the timing of partial melting in the GHS hanging wall (~23–19.5 Ma) overlaps with the timing of prograde metamorphism (~21 Ma) in the LHS footwall, confirming that the hanging wall may have provided the heat necessary for the metamorphism of the footwall. Overall, the data provide robust evidence for progressively downwards-penetrating deformation and accretion of original LHS footwall material to the GHS hanging wall over a period of ~5 Ma. These processes appear to have occurred several times during the prolonged ductile evolution of the thrust. The preserved inverted metamorphic sequence therefore documents the formation of sequential ‘paleothrusts’ through time, cutting down from the original locus of MCT movement at the LHS–GHS protolith boundary and forming at successively lower pressure and temperature conditions. The petrochronologic methods applied here constrain a complex temporal and thermal deformation history, and demonstrate that inverted metamorphic sequences can preserve a rich record of the duration of progressive ductile thrusting

Permian high-temperature metamorphism in the Western Alps (NW Italy)

During the late Palaeozoic, lithospheric thinning in part of the Alpine realm caused high-temperature low-to-medium pressure metamorphism and partial melting in the lower crust. Permian metamorphism and magmatism has extensively been recorded and dated in the Central, Eastern, and Southern Alps. However, Permian metamorphic ages in the Western Alps so far are constrained by very few and sparsely distributed data. The present study fills this gap. We present U/Pb ages of metamorphic zircon from several Adria-derived continental units now situated in the Western Alps, defining a range between 286 and 266 Ma. Trace element thermometry yields temperatures of 580-890°C from Ti-in-zircon and 630-850°C from Zr-in-rutile for Permian metamorphic rims. These temperature estimates, together with preserved mineral assemblages (garnet-prismatic sillimanite-biotite-plagioclase-quartz-K-feldspar-rutile), define pervasive upper-amphibolite to granulite facies conditions for Permian metamorphism. U/Pb ages from this study are similar to Permian ages reported for the Ivrea Zone in the Southern Alps and Austroalpine units in the Central and Eastern Alps. Regional comparison across the former Adriatic and European margin reveals a complex pattern of ages reported from late Palaeozoic magmatic and metamorphic rocks (and relics thereof): two late Variscan age groups (~330 and ~300 Ma) are followed seamlessly by a broad range of Permian ages (300-250 Ma). The former are associated with late-orogenic collapse; in samples from this study these are weakly represented. Clearly, dominant is the Permian group, which is related to crustal thinning, hinting to a possible initiation of continental rifting along a passive margin

Characterisation of a garnet population from the Sikkim Himalaya: insights into the rates and mechanisms of porphyroblast crystallisation

The compositional zoning of a garnet population contained within a garnet-grade metapelitic schist from the Lesser Himalayan Sequence of Sikkim (India) provides insight into the rates and kinetic controls of metamorphism, and the extent of chemical equilibration during porphyroblast crystallisation in the sample. Compositional profiles across centrally sectioned garnet crystals representative of the observed crystal size distribution indicate a strong correlation between garnet crystal size and core composition with respect to major end-member components. Systematic steepening of compositional gradients observed from large to small grains is interpreted to reflect a progressive decrease in the growth rate of relatively late-nucleated garnet as a result of an increase in interfacial energies during progressive crystallisation. Numerical simulation of garnet nucleation and growth using an equilibrium approach accounting for chemical fractionation associated with garnet crystallisation reproduces both the observed crystal size distribution and the chemical zoning of the entire garnet population. Simulation of multicomponent intracrystalline diffusion within the population indicates rapid heating along the pressure–temperature path, in excess of 100 ∘C Myr- 1. Radial garnet growth is correspondingly rapid, with minimum rates of 1.4 mm Myr- 1. As a consequence of such rapid crystallisation, the sample analysed in this study provides a close to primary record of the integrated history of garnet nucleation and growth. Our model suggests that nucleation of garnet occurred continuously between incipient garnet crystallisation at ∼ 520 ∘C, 4.5 kbar and peak metamorphic conditions at ∼ 565 ∘C, 5.6 kbar. The good fit between the observed and predicted garnet growth zoning suggests that the departure from equilibrium associated with garnet nucleation and growth was negligible, despite the particularly fast rates of metamorphic heating. Consequently, rates of major element diffusion in the intergranular medium during garnet crystallisation are interpreted to have been correspondingly rapid. It is, therefore, possible to simulate the prograde metamorphic history of our sample as a succession of equilibrium states of a chemical system modified by chemical fractionation associated with garnet crystallisation