Crystal chemistry and hydrogen bonding of rustumite Ca10(Si2O7)2(SiO4)(OH)2Cl2 with variable OH, Cl, F

Three samples of the skarn mineral rustumite Ca10(Si2O7)2(SiO4)(OH)2Cl2, space group C2/c, a ≈7.6, b ≈ 18.5, c ≈ 15.5 Å, β ≈ 104°, with variable OH, Cl, F content were investigated by electron microprobe, single-crystal X-ray structure refinements, and Raman spectroscopy. “Rust1LCl” is a low chlorine rustumite Ca10(Si2O7)2(SiO4)(OH1.88F0.12)(Cl1.28,OH0.72) from skarns associated with the Rize batholith near Ikizedere, Turkey. “Rust2F” is a F-bearing rustumite Ca10(Si2O7)2(SiO4)(OH1.13F0.87) (Cl1 96OH0.04) from xenoliths in ignimbrites of the Upper Chegem Caldera, Northern Caucasus, Russia. “Rust3LClF” represents a low-Cl, F-bearing rustumite Ca10(Si2O7)2(SiO4)0.87(H4O4)0.13(OH1.01F0.99) (Cl1.00 OH1.00) from altered merwinite skarns of the Birkhin massif, Baikal Lake area, Eastern Siberia, Russia. Rustumite from Birkhin massif is characterized by a significant hydrogarnet-like or fluorine substitution at the apices of the orthosilicate group, leading to specific atomic displacements. The crystal structures including hydrogen positions have been refined from single-crystal X-ray data to R1 = 0.0205 (Rust1_LCl), R1 = 0.0295 (Rust2_F), and R1 = 0.0243 (Rust3_LCl_F), respectively. Depletion in Cl and replacement by OH is associated with smaller unit-cell dimensions. The substitution of OH by F leads to shorter hydrogen bonds O-H⋯F instead of O-H⋯OH. Raman spectra for all samples have been measured and confirm slight strengthening of the hydrogen bonds with uptake of F.This study discusses the complex crystal chemistry of the skarn mineral rustumite and may provide a wider understanding of the chemical reactions related to contact metamorphism of limestones.

Bias associated with detrital zircon geochronology and thermochronometry

Detrital studies that utilize zircon U–Pb geochronology and fission-track (FT) thermochronometry are subject to a range of potential sources of bias that should be properly evaluated and minimized. Some of them are common to any single-grain mineral analysis (e.g., variable bedrock mineral fertility, hydraulic sorting during transport, selective grain loss during sample processing), whereas others are intrinsic to zircon, and are related to radiation damage and age discordance. In this article, we quantify the impact of intrinsic bias on detrital studies thanks to the analysis of modern detritus shed from the European Alps, and illustrate the general implications on geological interpretations. We show that detrital zircon U–Pb age distributions based on statistically robust datasets are highly reproducible and representative of the parent bedrock ages in the catchment. Arbitrary or selective removal of discordant grain ages can be minimized by using the Kolmogorov–Smirnov test to identify an appropriate cutoff level. Loss of metamict (α-damaged) zircon has a minor impact on data representativeness, and is mainly controlled by regional metamorphism rather than by mechanical abrasion during river transport. Zircon FT grain-age distributions were found to have poor reproducibility, although age spectra are consistent with bedrock data. However, unlike the U–Pb datasets, U-rich zircon grains (> 1000 ppm) are systematically missed, and undatable grains may exceed 50%. We identify two major sources of distribution bias specific to zircon FT datasets: (i) sediment sources dominated by U-rich zircon grains are markedly underrepresented in the detrital record, because such grains often have uncountable high densities of fission tracks (“U concentration bias”); (ii) sediment sources that shed zircon grains with high levels of α-damage are underrepresented, because these grains are lost during chemical etching for FT revelation (“etching bias”). In the case of multimethod dating on the same grains (e.g., FT and U–Pb double dating), bias affecting detrital zircon FT dating propagates to the entire dataset. These effects may not impact on exhumation-rate studies that utilize the youngest grain ages (i.e., lag-time approach). However, they represent a limiting factor for conventional provenance studies, and generally preclude application of zircon FT dating to sediment budget calculations.

Neotethys closure history of Anatolia: insights from 40Ar–39Ar geochronology and P–T estimation in high pressure metasedimentary rocks

The multiple high-pressure (HP), low-temperature (LT) metamorphic units of Western and Central Anatolia offer a great opportunity to investigate the subduction- and continental accretion-related evolution of the eastern limb of the long-lived Aegean subduction system. Recent reports of the HP–LT index mineral Fe-Mg-carpholite in three metasedimentary units of the Gondwana-derived Anatolide–Tauride continental block (namely the Afyon Zone, the Ören Unit and the southern Menderes Massif) suggest a more complicated scenario than the single-continental accretion model generally put forward in previous studies. This study presents the first isotopic dates (white mica 40Ar–39Ar geochronology), and where possible are combined with P–T estimates (chlorite thermometry, phengite barometry, multi-equilibrium thermobarometry), on carpholite-bearing rocks from these three HP–LT metasedimentary units. It is shown that, in the Afyon Zone, carpholite-bearing assemblages were retrogressed through greenschist-facies conditions at c. 67–62 Ma. Early retrograde stages in the Ören Unit are dated to 63–59 Ma. In the Kurudere–Nebiler Unit (HP Mesozoic cover of the southern Menderes Massif), HP retrograde stages are dated to c. 45 Ma, and post-collisional cooling to c. 26 Ma. These new results support that the Ören Unit represents the westernmost continuation of the Afyon Zone, whereas the Kurudere–Nebiler Unit correlates with the Cycladic Blueschist Unit of the Aegean Domain. In Western Anatolia, three successive HP–LT metamorphic belts thus formed: the northernmost Tavşanlı Zone (c. 88–82 Ma), the Ören–Afyon Zone (between 70 and 65 Ma), and the Kurudere–Nebiler Unit (c. 52–45 Ma). The southward younging trend of the HP–LT metamorphism from the upper and internal to the deeper and more external structural units, as in the Aegean Domain, points to the persistence of subduction in Western Anatolia between 93–90 and c. 35 Ma. After the accretion of the Menderes–Tauride terrane, in Eocene times, subduction stopped, leading to continental collision and associated Barrovian-type metamorphism. Because, by contrast, the Aegean subduction did remain active due to slab roll-back and trench migration, the eastern limb (below Southwestern Anatolia) of the Hellenic slab was dramatically curved and consequently teared. It therefore is suggested that the possibility for subduction to continue after the accretion of buoyant (e.g. continental) terranes probably depends much on palaeogeography.

REE and Hf distribution among mineral phases in the CV-CK clan : a way to explain present-day Hf isotopic variations in chondrites

Chondrites are among the most primitive objects in the Solar System and constitute the main building blocks of telluric planets. Among the radiochronometers currently used for dating geological events, Sm–Nd and Lu–Hf are both composed of refractory, lithophile element. They are thought to behave similarly as the parent elements (Sm and Lu) are generally less incompatible than the daughter elements (Nd and Hf) during geological processes. As such, their respective average isotopic compositions for the solar system should be well defined by the average of chondrites, called Chondritic Uniform Reservoir (CHUR). However, while the Sm–Nd isotopic system shows an actual spread of less than 4% in the average chondritic record, the Lu–Hf system shows a larger variation range of 28% [Bouvier A., Vervoort J. D. and Patchett P. J. (2008) The Lu–Hf and Sm–Nd isotopic composition of CHUR: Constraints from unequilibrated chondrites and implications for the bulk composition of terrestrial planets. Earth Planet. Sci. Lett.273, 48–57]. To better understand the contrast between Sm–Nd and Lu–Hf systems, the REE and Hf distribution among mineral phases during metamorphism of Karoonda (CK) and Vigarano-type (CV) carbonaceous chondrites has been examined. Mineral modes were determined from elemental mapping on a set of five CK chondrites (from types 3–6) and one CV3 chondrite. Trace-element patterns are obtained for the first time in all the chondrite-forming minerals of a given class (CK chondrites) as well as one CV3 sample. This study reveals that REE are distributed among both phosphates and silicates. Only 30–50% of Sm and Nd are stored in phosphates (at least in chondrites types 3–5); as such, they are not mobilized during early stages of metamorphism. The remaining fraction of Sm and Nd is distributed among the same mineral phases; these elements are therefore not decoupled during metamorphism. Of the whole-rock total of Lu, the fraction held in phosphate decreases significantly as the degree of metamorphism increases (30% for types 3 and 4, less than 5% in type 6). In contrast to Lu, Hf is mainly hosted by silicates with little contribution from phosphates throughout the CK metamorphic sequence. A significant part of Sm and Nd are stored in phosphates in types 3–5, and these elements behave similarly during CK chondrite metamorphism. That explains the robustness of the Sm/Nd ratios in chondrites through metamorphism, and the slight discrepancies observed in the present-day isotopic Nd values in chondrites. On the contrary, Lu and Hf are borne by several different minerals and consequently they are redistributed during metamorphism–induced recrystallization. The Lu/Hf ratios are therefore significantly disturbed during chondrites metamorphism, leading to the high discrepancies observed in present-day Hf isotopic values in chondrites.

Metamorphism, Transient Mid-Crustal Rheology, Strain Localization and the Exhumation of High-Grade Metamorphic Rocks

We present a series of three-dimensional numerical models investigating the effects of metamorphic strengthening and weakening on the geodynamic evolution of convergent orogens that are constrained by observations from an exposed mid-crustal section in the New England Appalachians. The natural mid-crustal section records evidence for spatially and temporally variable mid-crustal strength as a function of metamorphic grade during prograde polymetamorphism. Our models address changes in strain rate partitioning and topographic uplift as a function of strengthening/weakening in the middle crust, as well as the resultant changes in deformation kinematics and potential exhumation patterns of high-grade metamorphic rock. Results suggest that strengthening leads to strain rate partitioning around the zone and suppressed topographic uplift rates whereas weakening leads to strain rate partitioning into the zone and enhanced topographic uplift rates. Deformation kinematics recorded in the orogen are also affected by strengthening/weakening, with complete reversals in shear sense occurring as a function of strengthening/weakening without changes in plate boundary kinematics.

Dating polygenetic metamorphic assemblages along a transect across the Western Alps

Multichronometric analyses were performed on samples from a transect in the French-Italian Western Alps crossing nappes derived from the Briançonnais terrane and the Piemonte-Liguria Ocean, in an endeavour to constrain the high-pressure (HP) metamorphism and the retrogression history. 12 samples of white mica were analysed by 39Ar-40Ar stepwise heating, complemented by 2 samples from the Monte Rosa 100 km to the NE and also attributed to the Briançonnais terrane. One Sm-Nd and three Lu-Hf garnet ages from eclogites were also obtained. White mica ages decrease from ca. 300 Ma in the westernmost samples (Zone Houillère), reaching ca. 300 °C during Alpine metamorphism, to < 48 Ma in the internal units to the East, which reached ca. 500 °C during Alpine orogeny. The conventional “thermochronological” interpretation postulates Cretaceous Eo-Alpine HP metamorphism and younger “cooling ages” in the higher-temperature samples. However, Eocene Lu-Hf and Sm-Nd ages from the same samples cannot be interpreted as post-metamorphic cooling ages, which makes a Cretaceous eclogitization untenable. The age date from this transect require instead to replace conventional “thermochronology” by an approach combining age dating with detailed geochemical, petrological and microstructural investigations. Petrology reveals important mineralogical differences along the transect. Samples from the Zone Houillère mostly contain detrital mica. White mica with Si > 6.45 atoms per formula unit becomes more abundant eastward. Across the whole traverse, HP phengitic mica forms the D1 foliation. Syn-D2 mica is Si-poorer and associated with nappe stacking, exhumation, and hydrous retrogression under greenschist facies conditions. D1 phengite is very often corroded, overgrown or intergrown by syn-D2 muscovite. Most importantly, syn-D2 recrystallization is not limited to S2 schistosity domains; microchemical fingerprinting shows that it also can form pseudomorphs after crystals that could be mistaken to have formed during D1 based on microstructural arguments alone. Thereby the Cl concentration in white mica is a useful discriminator, since D2 retrogression was associated with a less saline fluid than eclogitization. Once the petrological stage is set, geochronology is straightforward. All samples contain mixtures of detrital, syn-D1 and syn-D2 mica, and retrogression phases (D3) in greatly varying proportions according to local pressure-temperature-fluid activity-deformation conditions. The correlation of age vs. Cl/K clearly identifies 47 ± 1 Ma as the age of formation of syn-D1 mica along the entire transect, including the Monte Rosa nappe samples. The inferred age of the greenschist-facies low-Si syn-D2 mica generation ranges within 39-43 Ma, with local variations. Coexistence of D1 and D2 ages, and the constancy of non-reset D1 ages along the entire transect, are strong evidence that the D1 white mica ages are very close to formation ages. Volume diffusion of Ar in white mica (activation energy E = 250 kJ/mol; pressure-adjusted diffusion coefficient D’0 < 0.03 cm2 s-1) has a subordinate effect on mineral ages compared to both prograde and retrograde recrystallization in most samples. Eocene Lu-Hf and Sm-Nd garnet ages are prograde and predate the HP peak.

Ediacaran 2,500-km-long synchronous deep continental subduction in the West Gondwana Orogen

The deeply eroded West Gondwana Orogen is a major continental collision zone that exposes numerous occurrences of deeply subducted rocks, such as eclogites. The position of these eclogites marks the suture zone between colliding cratons, and the age of metamorphism constrains the transition from subduction-dominated tectonics to continental collision and mountain building. Here we investigate the metamorphic conditions and age of high-pressure and ultrahigh-pressure eclogites from Mali, Togo and NE-Brazil and demonstrate that continental subduction occurred within 20 million years over at least a 2,500-km-long section of the orogen during the Ediacaran. We consider this to be the earliest evidence of large-scale deep-continental subduction and consequent appearance of Himalayan-scale mountains in the geological record. The rise and subsequent erosion of such mountains in the Late Ediacaran is perfectly timed to deliver sediments and nutrients that are thought to have been necessary for the subsequent evolution of sustainable life on Earth.

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.

The mobility of Nb in rutile-saturated NaCl- and NaF-bearing aqueous fluids from 1-6.5 GPa and 300-800 °C

Rutile (TiO2) is an important host phase for high field strength elements (HFSE) such as Nb in metamorphic and subduction zone environments. The observed depletion of Nb in arc rocks is often explained by the hypothesis that rutile sequesters HFSE in the subducted slab and overlying sediment, and is chemically inert with respect to aqueous fluids evolved during prograde metamorphism in the forearc to subarc environment. However, field observations of exhumed terranes, and experimental studies, indicate that HFSE may be soluble in complex aqueous fluids at high pressure (i.e., >0.5 GPa) and moderate to high temperature (i.e., >300 degrees C). In this study, we investigated experimentally the mobility of Nb in NaCl- and NaF-bearing aqueous fluids in equilibrium with Nb-bearing rutile at pressure-temperature conditions applicable to fluid evolution in arc environments. Niobium concentrations in aqueous fluid at rutile saturation were measured directly by using a hydrothermal diamond-anvil cell (HDAC) and synchrotron X-ray fluorescence (SXRF) at 2.1 to 6.5 GPa and 300-500 degrees C, and indirectly by performing mass loss experiments in a piston-cylinder (PC) apparatus at similar to 1 GPa and 700-800 degrees C. The concentration of Nb in a 10 wt% NaCl aqueous fluid increases from 6 to 11 mu g/g as temperature increases from 300 to 500 degrees C, over a pressure range from 2.1 to 2.8 GPa, consistent with a positive temperature dependence. The concentration of Nb in a 20 wt% NaCl aqueous fluid varies from 55 to 150 mu g/g at 300 to 500 degrees C, over a pressure range from 1.8 to 6.4 GPa; however, there is no discernible temperature or pressure dependence. The Nb concentration in a 4 wt% NaF-bearing aqueous fluid increases from 180 to 910 mu g/g as temperature increases from 300 to 500 degrees C over the pressure range 2.1 to 6.5 GPa. The data for the F-bearing fluid indicate that the Nb content of the fluid exhibits a dependence on temperature between 300 and 500 degrees C at >= 2 GPa, but there is no observed dependence on pressure. Together, the data demonstrate that the hydrothermal mobility of Nb is strongly controlled by the composition of the fluid, consistent with published data for Ti. At all experimental conditions, however, the concentration of Nb in the fluid is always lower than coexisting rutile, consistent with a role for rutile in moderating the Nb budget of arc rocks.

Fluid-related inclusions in Alpine high-pressure peridotite reveal trace element recycling during subduction-zone dehydration of serpentinized mantle (Cima di Gagnone, Swiss Alps).

Serpentinites release at sub-arc depths volatiles and several fluid-mobile trace elements found in arc magmas. Constraining element uptake in these rocks and defining the trace element composition of fluids released upon serpentinite dehydration can improve our understanding of mass transfer across subduction zones and to volcanic arcs. The eclogite-facies garnet metaperidotite and chlorite harzburgite bodies embedded in paragneiss of the subduction melange from Cima di Gagnone derive from serpentinized peridotite protoliths and are unique examples of ultramafic rocks that experienced subduction metasomatism and devolatilization. In these rocks, metamorphic olivine and garnet trap polyphase inclusions representing the fluid released during high-pressure breakdown of antigorite and chlorite. Combining major element mapping and laser-ablation ICP-MS bulk inclusion analysis, we characterize the mineral content of polyphase inclusions and quantify the fluid composition. Silicates, Cl-bearing phases, sulphides, carbonates, and oxides document post-entrapment mineral growth in the inclusions starting immediately after fluid entrapment. Compositional data reveal the presence of two different fluid types. The first (type A) records a fluid prominently enriched in fluid-mobile elements, with Cl, Cs, Pb, As, Sb concentrations up to 10(3) PM (primitive mantle), similar to 10(2) PM Tit Ba, while Rb, B, Sr, Li, U concentrations are of the order of 10(1) PM, and alkalis are similar to 2 PM. The second fluid (type B) has considerably lower fluid-mobile element enrichments, but its enrichment patterns are comparable to type A fluid. Our data reveal multistage fluid uptake in these peridotite bodies, including selective element enrichment during seafloor alteration, followed by fluid-rock interaction along with subduction metamorphism in the plate interface melange. Here, infiltration of sediment-equilibrated fluid produced significant enrichment of the serpentinites in As, Sb, B, Pb, an enriched trace element pattern that was then transferred to the fluid released at greater depth upon serpentine dehydration (type A fluid). The type B fluid hosted by garnet may record the composition of the chlorite breakdown fluid released at even greater depth. The Gagnone study-case demonstrates that serpentinized peridotites acquire water and fluid-mobile elements during ocean floor hydration and through exchange with sediment-equilibrated fluids in the early subduction stages. Subsequent antigorite devolatilization at subarc depths delivers aqueous fluids to the mantle wedge that can be prominently enriched in sediment-derived components, potentially triggering arc magmatism without the need of concomitant dehydration/melting of metasediments or altered oceanic crust.

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.

Provenance analysis using Raman spectroscopy of carbonaceous material: A case study in the Southern Alps of New Zealand

Detrital provenance analyses in orogenic settings, in which sediments are collected at the outlet of a catchment, have become an important tool to estimate how erosion varies in space and time. Here we present how Raman Spectroscopy on Carbonaceous Material (RSCM) can be used for provenance analysis. RSCM provides an estimate of the peak temperature (RSCM-T) experienced during metamorphism. We show that we can infer modern erosion patterns in a catchment by combining new measurements on detrital sands with previously acquired bedrock data. We focus on the Whataroa catchment in the Southern Alps of New Zealand and exploit the metamorphic gradient that runs parallel to the main drainage direction. To account for potential sampling biases, we also quantify abrasion properties using flume experiments and measure the total organic carbon content in the bedrock that produced the collected sands. Finally, we integrate these parameters into a mass-conservative model. Our results first demonstrate that RSCM-T can be used for detrital studies. Second, we find that spatial variations in tracer concentration and erosion have a first-order control on the RSCM-T distributions, even though our flume experiments reveal that weak lithologies produce substantially more fine particles than do more durable lithologies. This result implies that sand specimens are good proxies for mapping spatial variations in erosion when the bedrock concentration of the target mineral is quantified. The modeling suggests that highest present-day erosion rates (in Whataroa catchment) are not situated at the range front but around 10 km into the mountain belt.

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.

Contrasting PT paths within the Barchi-Kol UHP terrain (Kokchetav Complex): Implications for subduction and exhumation of continental crust

The Barchi-Kol terrain is a classic locality of ultrahigh-pressure (UHP) metamorphism within the Kokchetav metamorphic belt. We provide a detailed and systematic characterization of four metasedimentary samples using dominant mineral assemblages, mineral inclusions in zircon and monazite, garnet zonation with respect to major and trace elements, and Zr-in-rutile and Ti-in-zircon temperatures. A typical diamond-bearing gneiss records peak conditions of 49 ± 4 kbar and 950–1000 °C. Near isothermal decompression of this rock resulted in the breakdown of phengite associated with a pervasive recrystallization of the rock. The same terrain also contains mica schists that experienced peak conditions close to those of the diamond-bearing rocks, but they were exhumed along a cooler path where phengite remained stable. In these rocks, major and trace element zoning in garnet has been completely equilibrated. A layered gneiss was metamorphosed at UHP conditions in the coesite field, but did not reach diamond-facies conditions (peak conditions: 30 kbar and 800–900 °C). In this sample, garnet records retrograde zonation in major elements and also retains prograde zoning in trace elements. A garnet-kyanite-micaschist that reached significantly lower pressures (24 ± 2 kbar, 710 ± 20 °C) contains garnet with major and trace element zoning. The diverse garnet zoning in samples that experienced different metamorphic conditions allows to establish that diffusional equilibration of rare earth element in garnet likely occurs at ~900–950 °C. Different metamorphic conditions in the four investigated samples are also documented in zircon trace element zonation and mineral inclusions in zircon and monazite. U-Pb geochronology of metamorphic zircon and monazite domains demonstrates that prograde (528–521 Ma), peak (528–522 Ma), and peak to retrograde metamorphism (503–532 Ma) occurred over a relatively short time interval that is indistinguishable from metamorphism of other UHP rocks within the Kokchetav metamorphic belt. Therefore, the assembly of rocks with contrasting P-T trajectories must have occurred in a single subduction-exhumation cycle, providing a snapshot of the thermal structure of a subducted continental margin prior to collision. The rocks were initially buried along a low geothermal gradient. At 20–25 kbar they underwent near isobaric heating of 200 °C, which was followed by continued burial along a low geothermal gradient. Such a step-wise geotherm is in good agreement with predictions from subduction zone thermal models.