A thermodynamic model for di-trioctahedral chlorite from experimental and natural data in the system MgO-FeO-Al2O3-SiO2-H2O. Applications to P-T sections and geothermometry

We present a new thermodynamic activity-composition model for di-trioctahedral chlorite in the system FeO–MgO–Al2O3–SiO2–H2O that is based on the Holland–Powell internally consistent thermodynamic data set. The model is formulated in terms of four linearly independent end-members, which are amesite, clinochlore, daphnite and sudoite. These account for the most important crystal-chemical substitutions in chlorite, the Fe–Mg, Tschermak and di-trioctahedral substitution. The ideal part of end-member activities is modeled with a mixing-on-site formalism, and non-ideality is described by a macroscopic symmetric (regular) formalism. The symmetric interaction parameters were calibrated using a set of 271 published chlorite analyses for which robust independent temperature estimates are available. In addition, adjustment of the standard state thermodynamic properties of sudoite was required to accurately reproduce experimental brackets involving sudoite. This new model was tested by calculating representative P–T sections for metasediments at low temperatures (<400 °C), in particular sudoite and chlorite bearing metapelites from Crete. Comparison between the calculated mineral assemblages and field data shows that the new model is able to predict the coexistence of chlorite and sudoite at low metamorphic temperatures. The predicted lower limit of the chloritoid stability field is also in better agreement with petrological observations. For practical applications to metamorphic and hydrothermal environments, two new semi-empirical chlorite geothermometers named Chl(1) and Chl(2) were calibrated based on the chlorite + quartz + water equilibrium (2 clinochlore + 3 sudoite = 4 amesite + 4 H2O + 7 quartz). The Chl(1) thermometer requires knowledge of the (Fe3+/ΣFe) ratio in chlorite and predicts correct temperatures for a range of redox conditions. The Chl(2) geothermometer which assumes that all iron in chlorite is ferrous has been applied to partially recrystallized detrital chlorite from the Zone houillère in the French Western Alps.

Multiple metamorphic stages within an eclogite-facies terrane (Sesia Zone, Western Alps) revealed by Th-U-Pb petrochronology

Convergent plate margins typically experience a transition from subduction to collision dynamics as massive continental blocks enter the subduction channel. Studies of high-pressure rocks indicate that tectonic fragments are rapidly exhumed from eclogite facies to midcrustal levels, but the details of such dynamics are controversial.To understand the dynamics of a subduction channel we report the results of a petrochronological study from the central Sesia Zone, a key element of the internalWestern Alps.This comprises two polymetamorphic basement complexes (Eclogitic Micaschist Complex and Gneiss Minuti Complex) and a thin, dismembered cover sequence (Scalaro Unit) associated with pre-Alpine metagabbros and metasediments (Bonze Unit). Structurally controlled samples from three of these units (Eclogitic Micaschist Complex and Scalaro-Bonze Units) yield unequivocal petrological and geochronological evidence of two distinct high-pressure stages. Ages (U-Th-Pb) of growth zones in accessory allanite and zircon, combined with inclusion and textural relationships, can be tied to the multi-stage evolution of single samples.Two independent tectono-metamorphic ‘slices’ showing a coherent metamorphic evolution during a given time interval have been recognized: the Fondo slice (which includes Scalaro and Bonze rocks) and the Druer slice (belonging to the Eclogitic Micaschist Complex).The new data indicate separate stages of deformation at eclogite-facies conditions for each recognized independent kilometer-sized tectono-metamorphic slice, between ~85 and 60 Ma, with evidence of intermittent decompression (∆P~0.5 GPa) within only the Fondo slice. The evolution path of the Druer slice indicates a different P-T-time evolution with prolonged eclogite-facies metamorphism between ~85 and 75Ma. Our approach, combining structural, petrological and geochronological techniques, yields field-based constraints on the duration and rates of dynamics within a subduction channel.

Oxygen isotope record of oceanic and high-pressure metasomatism: a P-T-time-fluid path for the Monviso eclogites (Italy)

Fluids are considered a fundamental agent for chemical exchanges between different rock types in the subduction system. Constraints on the sources and pathways of subduction fluids thus provide crucial information to reconstruct subduction processes. The Monviso ophiolitic sequence is composed of mafic, ultramafic and minor sediments that have been subducted to ~80 km depth. In this sequence, both localized fluid flow and channelized fluids along major shear zones have been documented. We investigate the timing and source of the fluids that affected the dominant mafic rocks using microscale U-Pb dating of zircon and oxygen isotope analysis of mineral zones (garnet, zircon and antigorite) in high pressure rocks with variable degree of metasomatic modification. In mafic eclogites, Jurassic zircon cores are the only mineralogical relicts of the protolith gabbros and retain δ18O values of 4.5–6 ‰, typical of mantle melts. Garnet and metamorphic zircon that grew during prograde to peak metamorphism display low δ18O values between 0.2 and 3.8 ‰, which are likely inherited from high-temperature alteration of the protolith on the sea floor. This is corroborated by δ18O values of 3.0 and 3.6 ‰ in antigorite from surrounding serpentinites. In metasomatised eclogites within the Lower Shear Zone, garnet rim formed at the metamorphic peak shows a shift to higher δ18O up to 6‰. The age of zircons in high-pressure veins and metasomatised eclogites constrains the timing of fluid flow at high pressure at around 45–46 Ma. Although the oxygen data do not contradict previous reports of interaction with serpentinite-derived fluids, the shift to isotopically heavier oxygen compositions requires contribution from sediment-derived fluids. The scarcity of metasediments in the Monviso sequence suggests that such fluids were concentrated and fluxed along the Lower Shear Zone in a sufficient amount to modify the oxygen composition of the eclogitic minerals.

UHP Metamorphism Documented in Ti-chondrodite- and Ti-clinohumite-bearing Serpentinized Ultramafic Rocks from Chinese Southwestern Tianshan

The southwestern Tianshan (China) metamorphic belt records high-pressure (HP) to ultrahigh-pressure (UHP) conditions corresponding to a cold oceanic subduction-zone setting. Serpentinites enclosing retrogressed eclogite and rodingite occur as lenses within metapelites in the UHP unit, which also hosts coesite-bearing eclogites. Based on the petrology and petrography of these serpentinites, five events are recognized: (1) formation of a wehrlite–harzburgite–dunite association in the mantle; (2) retrograde metamorphism and partial hydration during exhumation of the mantle rocks close to the seafloor; (3) oceanic metamorphism leading to the first serpentinization and rodingitization; (4) UHP metamorphism during subduction; (5) retrograde metamorphism during exhumation together with a second serpentinization. The peak metamorphic mineral assemblage of the serpentinized wehrlite comprises Ti-chondrodite + olivine + antigorite + chlorite + magnetite + brucite. A computed pseudosection for this serpentinized wehrlite shows that the Al content in antigorite is mostly sensititive to temperature but can also be used to constrain pressure. The average XAl = 0·204 ± 0·026 of antigorite (XAl = Al (a.p.f.u.)/8, where Al is in atoms per formula unit for a structural formula M48T34O85(OH)62, and M and T are octahedral and tetrahedral sites, respectively) included in Ti-chondrodite and average XAl = 0·203 ± 0·019 of antigorite in the matrix result in a well-constrained peak metamorphic temperature of 510–530°C. Peak pressures are less precisely constrained at 37 ± 7 kbar. The Tianshan serpentinites thus record UHP metamorphic conditions and represent the deepest subducted serpentinites discovered so far. The retrograde evolution occurs within the stability field of brucite + antigorite + olivine + chlorite and formation of Ti-clinohumite at the expense of Ti-chondrodite has been observed, suggesting isothermal decompression. The resulting P–T path is in excellent agreement with the metamorphic evolution of country rocks, indicating that the UHP unit in Tianshan was subducted and exhumed as a coherent block. To refine the metamorphic path of the ultramafic rocks, we have investigated the stability fields of Ti-chondrodite and Ti-clinohumite using piston-cylinder experiments. A total of 11 experiments were conducted at 25–55 kbar and 600–750°C in a F-free natural system. Combined with previous experiments and information from natural rocks we constructed a petrogenetic grid for the stability of Ti-chondrodite and Ti-clinohumite in F-free peridotite compositions. The formation of Ti-chondrodite in serpentinites requires a minimum pressure of about 26 kbar, whereas in Ti-rich systems it can form at considerably lower pressures. A key finding is that at UHP conditions, F-free Ti-chondrodite or Ti-clinohumite breaks down in the presence of orthopyroxene between 700 and 750°C, at temperatures that are significantly lower than those of the terminal breakdown reactions of these humite minerals. These breakdown reactions are an additional source of fluid during prograde subduction of serpentinites.

Dating prograde fluid pulses during subduction by in situ U–Pb and oxygen isotope analysis

Keywords High-pressure fluids · Whiteschists · U–Pb dating · Oxygen isotopes · Ion microprobe · Metasomatism Introduction The subduction of crustal material to mantle depths and its chemical modification during burial and exhumation contribute to element recycling in the mantle and the formation of new crust through arc magmatism. Crustal rocks that Abstract The Dora-Maira whiteschists derive from metasomatically altered granites that experienced ultrahighpressure metamorphism at ~750 °C and 40 kbar during the Alpine orogeny. In order to investigate the P–T–time– fluid evolution of the whiteschists, we obtained U–Pb ages from zircon and monazite and combined those with trace element composition and oxygen isotopes of the accessory minerals and coexisting garnet. Zircon cores are the only remnants of the granitic protolith and still preserve a Permian age, magmatic trace element compositions and δ18O of ~10 ‰. Thermodynamic modelling of Si-rich and Si-poor whiteschist compositions shows that there are two main fluid pulses during prograde subduction between 20 and 40 kbar. In Si-poor samples, the breakdown of chlorite to garnet + fluid occurs at ~22 kbar. A first zircon rim directly overgrowing the cores has inclusions of prograde phlogopite and HREE-enriched patterns indicating zircon growth at the onset of garnet formation. A second main fluid pulse is documented close to peak metamorphic conditions in both Si-rich and Si-poor whiteschist when talc + kyanite react to garnet + coesite + fluid. A second metamorphic overgrowth on zircon with HREE depletion was observed in the Si-poor whiteschists, whereas a single metamorphic overgrowth capturing phengite and talc inclusions was observed in the Si-rich whiteschists. Garnet rims, zircon rims and monazite are in chemical and isotopic equilibrium for oxygen, demonstrating that they all formed at peak metamorphism at 35 Ma as constrained by the age of monazite (34.7 ± 0.4 Ma) and zircon rims (35.1 ± 0.8 Ma). The prograde zircon rim in Si-poor whiteschists has an age that is within error indistinguishable from the age of peak metamorphic conditions, consistent with a minimum rate of subduction of 2 cm/year for the Dora-Maira unit. Oxygen isotope values for zircon rims, monazite and garnet are equal within error at 6.4 ± 0.4 ‰, which is in line with closed-system equilibrium fractionation during prograde to peak temperatures. The resulting equilibrium Δ18Ozircon-monazite at 700 ± 20 °C is 0.1 ± 0.7 ‰. The in situ oxygen isotope data argue against an externally derived input of fluids into the whiteschists. Instead, fluidassisted zircon and monazite recrystallisation can be linked to internal dehydration reactions during prograde subduction. We propose that the major metasomatic event affecting the granite protolith was related to hydrothermal seafloor alteration post-dating Jurassic rifting, well before the onset of Alpine subduction.

Constraints on the thermal evolution of the Adriatic margin during Jurassic continental break-up: U–Pb dating of rutile from the Ivrea–Verbano Zone, Italy

The Ivrea–Verbano Zone (IVZ), northern Italy, exposes an attenuated section through the Permian lower crust that records high-temperature metamorphism under lower crustal conditions and a protracted history of extension and exhumation associated partly with the Jurassic opening of the Alpine Tethys ocean. This study presents SHRIMP U–Pb geochronology of rutile from seven granulite facies metapelites from the base of the IVZ, collected from locations spanning ~35 km along the strike of Paleozoic fabrics. Rutile crystallised during Permian high-temperature metamorphism and anatexis, yet all samples give Jurassic rutile U–Pb ages that record cooling through 650–550 °C. Rutile age distributions are dominated by a peak at ~160 Ma, with a subordinate peak at ~175 Ma. Both ~160 and ~175 Ma age populations show excellent agreement between samples, indicating that the two distinctive cooling stages they record were synchronous on a regional scale. The ~175 Ma population is interpreted to record cooling in the footwall of rift-related faults and shear zones, for which widespread activity in the Lower Jurassic has been documented along the western margin of the Adriatic plate. The ~160 Ma age population postdates the activity of all known rift-related structures within the Adriatic margin, but coincides with extensive gabbroic magmatism and exhumation of sub-continental mantle to the floor of the Alpine Tethys, west of the Ivrea Zone. We propose that this ~160 Ma early post-rift age population records regional cooling following episodic heating of the distal Adriatic margin, likely related to extreme lithospheric thinning and associated advection of the asthenosphere to shallow levels. The partial preservation of the ~175 Ma age cluster suggests that the post-rift (~160 Ma) heating pulse was of short duration. The regional consistency of the data presented here, which is in contrast to many other thermochronometers in the IVZ, demonstrates the value of the rutile U–Pb technique for probing the thermal evolution of high-grade metamorphic terrains. In the IVZ, a significant decoupling between Zr-in-rutile temperatures and U–Pb ages of rutile is observed, with the two systems recording events ~120 Ma apart.

A systematic evaluation of the Zr-in-rutile thermometer in ultra-high temperature (UHT) rocks

The Zr-in-rutile geothermometer is potentially a widely applicable tool to estimate peak metamorphic temperatures in rocks from diverse geological settings. In order to evaluate its usefulness and reliability to record and preserve high temperatures in granulite facies rocks, rutile from UHT rocks was investigated to assess different mechanisms of Zr (re-)distribution following cooling from high temperature. Granulite facies paragneisses from the lowermost part of the Ivrea Zone, Italy, incorporated as thin sheets into the extensive basaltic body of the Mafic Complex were selected for this study. The results show that Zr-in-rutile thermometry, if properly applied, is well suited to identify and study UHT terranes as it preserves a record of temperatures up to 1190 °C, although the thermometer is susceptible to partial post-peak metamorphic resetting by Zr diffusion. Texturally homogeneous rutile grains preserve Zr concentrations corresponding to temperatures of prograde rutile growth. Diverse rutile textures and relationships between some rutile host grains and included or adjacent Zr-bearing phases bear testimony to varying mechanisms of partial redistribution and resetting of Zr in rutile during cooling and link Zr-in-rutile temperatures to different steps of the metamorphic evolution. Rutile grains that equilibrated their Zr concentrations at temperatures above 1070 °C (i.e. 1.1 wt% Zr) could not retain all Zr in the rutile structure during cooling and exsolved baddeleyite (ZrO2). By subsequent reaction of baddeleyite exsolution lamellae with SiO2, zircon needles formed before the system finally closed at 650–700 °C without significant net loss of Zr from the whole host rutile grain. By reintegration of zircon exsolution needles, peak metamorphic temperatures of up to 1190 °C are derived for the studied rocks, which demonstrates the suitability of this solution thermometer to record UHT conditions and also confirms the extraordinary geological setting of the lowermost part of the Ivrea Zone.

The Interplay between Melting, Refertilization and Carbonatite Metasomatism in Off-Cratonic Lithospheric Mantle under Zealandia: an Integrated Major, Trace and Platinum Group Element Study

It is widely accepted that stabilization of the continental crust requires the presence of sub-continental lithospheric mantle. However, the degree of melt depletion required to stabilize the lithosphere and whether widespread refertilization is a significant process remain unresolved. Here, major and trace element, including platinum group elements (PGE), characterization of 40 mantle xenoliths from 13 localities is used to constrain the melt depletion, refertilization and metasomatic history of lithospheric mantle underneath the micro-continent Zealandia. Our previously published Re–Os isotopic data for a subset of these xenoliths indicate Phanerozoic to Paleoproterozoic ages and, reinterpreted with the new major and trace element data presented here, demonstrate that a large volume (>2 million km3) of lithospheric mantle with an age of 1·99 ± 0·21 Ga is present below the much younger crust of Zealandia. A peritectic melting model using moderately incompatible trace elements (e.g. Yb) in bulk-rocks demonstrates that these peridotites experienced a significant range of degrees of partial melting, between 3 and 28%. During subsolidus equilibration clinopyroxene gains significant rare earth elements (REE), which then leads to the underestimation of the degree of partial melting by ≤12% in fertile xenoliths. A new approach taking into account the effects of subsolidus re-equilibration on clinopyroxene composition effectively removes discrepancies in the calculated degree of melting and provides consistent estimates of between 4 and 29%. The estimated amount of melting is independent of the Re–Os model ages of the samples. The PGE patterns record simple melt depletion histories and the retention of primary base metal sulfides in the majority of the xenoliths. A rapid decrease in Pt/IrN observed at c. 1·0 wt % Al2O3 is a direct result of the exhaustion of sulfide in the mantle residue at c. 20–25% partial melting and the inability of Pt to form a stable alloy phase. Major elements preserve evidence for refertilization by a basaltic component that resulted in the formation of secondary clinopyroxene and low-forsterite olivine. The majority of xenoliths show the effects of cryptic metasomatic overprinting, ranging from minor to strong light REE enrichments in bulk-rocks (La/YbN = 0·16–15·9). Metasomatism is heterogeneous, with samples varying from those with weak REE enrichment and notable positive Sr and U–Th anomalies and negative Nb–Ta anomalies in clinopyroxene to those that have extremely high concentrations of REE, Th–U and Nb. Chemical compositions are consistent with a carbonatitic component contributing to the metasomatism of the lithosphere under Zealandia. Notably, the intense metasomatism of the samples did not affect the PGE budget of the peridotites as this was controlled by residual sulfides.