Geochronology of metasomatic events

In order to date any geological event, suitable mineral geochronometers that record that and only that event must be identified and analyzed. In the case of metasomatism, recrystallisation is a key process that controls both the petrology and the isotopic record of minerals. It can occur both in the form of complete neocrystallisation (e.g. in a vein) and in the form of pseudomorphism, whereby dissolution/reprecipitation at the submicroscopic scale plays a central role. Recrystallisation may be complete or not, raising the possibility that relicts of a pre-metasomatic assemblage may be preserved. Because recrystallisation is energetically less costly at almost any temperature than diffusion, and because radiogenic isotopes (except 4He) never diffuse faster than major elements forming the mineral structure, there is a strong causal link between petrographic relicts and isotopic inheritance (as demonstrated for zircon, monazite, titanite, amphibole, K-feldspar, biotite, and muscovite). Metasomatic assemblages commonly contain such mixtures between relicts and newly formed phases, whose geochronology is slightly more complex than that of simple, ideal systems, but can be managed by techniques that have become routine in the last decade and which are described in this chapter. Because of its crucial role in controlling the isotope systematics, the petrogenesis of a mineral needs to be understood in extreme detail, especially using microchemical analyses and micro-imaging techniques, before mineral ages can be correctly interpreted. As the occurrence of recrystallization is limited by the availability of water, minerals act as “geohygrometers” that allow constraints to be placed on the nature and age of fluid circulation episodes, especially metasomatic events.

Ar diffusion in K-feldspar: present and absent

The mechanisms of Ar release from K-feldspar samples in laboratory experiments and during their geological history are assessed here. Modern petrology clearly established that the chemical and isotopic record of minerals is normally dominated by aqueous recrystallization. The laboratory critique is trickier, which explains why so many conflicting approaches have been able to survive long past their expiration date. Current models are evaluated for self-consistency; especially Arrhenian non-linearity leads to paradoxes. The models’ testable geological predictions suggest that temperature-based downslope extrapolations often overestimate observed geological Ar mobility substantially. An updated interpretation is based on the unrelatedness of geological behaviour to laboratory experiments. The isotopic record of K-feldspar in geological samples is not a unique function of temperature, as recrystallisation promoted by aqueous fluids is the predominant mechanism controlling isotope transport. K-feldspar should therefore be viewed as a hygrochronometer. Laboratory degassing proceeds from structural rearrangements and phase transitions such as are observed in situ at high temperature in Na and Pb feldspars. These effects violate the mathematics of an inert Fick’s Law matrix and preclude downslope extrapolation. The similar upward-concave, non-linear shapes of Arrhenius trajectories of many silicates, hydrous and anhydrous, are likely common manifestations of structural rearrangements in silicate structures.

A re-examination of petrogenesis and 40Ar/39Ar systematics in the the Chain of Ponds K-feldspar: "diffusion domain" archetype versus polyphase hygrochronology

K-feldspar (Kfs) from the Chain of Ponds Pluton (CPP) is the archetypal reference material, on which thermochronological modeling of Ar diffusion in discrete “domains” was founded. We re-examine the CPP Kfs using cathodoluminescence and back-scattered electron imaging, transmission electron microscopy, and electron probe microanalysis. 40Ar/39Ar stepwise heating experiments on different sieve fractions, and on handpicked and unpicked aliquots, are compared. Our results reproduce the staircase-shaped age spectrum and the Arrhenius trajectory of the literature sample, confirming that samples collected from the same locality have an identical Ar isotope record. Even the most pristine-looking Kfs from the CPP contains successive generations of secondary, metasomatic/retrograde mineral replacements that post-date magmatic crystallization. These chemically and chronologically distinct phases are responsible for its staircase-shaped age spectra, which are modified by handpicking. While genuine within-grain diffusion gradients are not ruled out by these data, this study demonstrates that the most important control on staircase-shaped age spectra is the simultaneous presence of heterochemical, diachronous post-magmatic mineral growth. At least five distinct mineral species were identified in the Kfs separate, three of which can be traced to external fluids interacting with the CPP in a chemically open system. Sieve fractions have size-shifted Arrhenius trajectories, negating the existence of the smallest “diffusion domains”. Heterochemical phases also play an important role in producing non-linear trajectories. In vacuo degassing rates recovered from Arrhenius plots are neither related to true Fick’s Law diffusion nor to the staircase shape of the age spectra. The CPP Kfs used to define the "diffusion domain" model demonstrates the predominance of metasomatic alteration by hydrothermal fluids and recrystallization in establishing the natural Ar distribution amongst different coexisting phases that gives rise to the staircase-shaped age spectrum. Microbeam imaging of textures is as essential for 40Ar-39Ar hygrochronology as it is for U-Pb geochronology.

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.

Probing the Allende meteorite with a miniature laser-ablation mass analyser for space application

We measured the elemental composition on a sample of Allende meteorite with a miniature laser ablation mass spectrometer. This Laser Mass Spectrometer (LMS) has been designed and built at the University of Bern in the Department of Space Research and Planetary Sciences with the objective of using such an instrument on a space mission. Utilising the meteorite Allende as the test sample in this study, it is demonstrated that the instrument allows the in situ determination of the elemental composition and thus mineralogy and petrology of untreated rocky samples, particularly on planetary surfaces. In total, 138 measurements of elemental compositions have been carried out on an Allende sample. The mass spectrometric data are evaluated and correlated with an optical image. It is demonstrated that by illustrating the measured elements in the form of mineralogical maps, LMS can serve as an element imaging instrument with a very high spatial resolution of µm scale. The detailed analysis also includes a mineralogical evaluation and an investigation of the volatile element content of Allende. All findings are in good agreement with published data and underline the high sensitivity, accuracy and capability of LMS as a mass analyser for space exploration.

Melting of metasomatized peridotite at 4–6 GPa and up to 1200 °C: an experimental approach

The phase assemblages and compositions in a K-bearing lherzolite + H2O system are determined between 4 and 6 GPa and 850–1200 °C, and the melting reactions occurring at subarc depth in subduction zones are constrained. Experiments were performed on a rocking multi-anvil apparatus. The experiments had around 16 wt% water content, and hydrous melt or aqueous fluid was segregated and trapped in a diamond aggregate layer. The compositions of the aqueous fluid and hydrous melt phases were measured using the cryogenic LA-ICP-MS technique. The residual lherzolite consists of olivine, orthopyroxene, clinopyroxene, and garnet, while diamond (C) is assumed to be inert. Hydrous and alkali-rich minerals were absent from the run products due to preferred dissolution of K2O (and Na2O) to the aqueous fluid/hydrous melt phases. The role of phlogopite in melting relations is, thus, controlled by the water content in the system: at the water content of around 16 wt% used here, phlogopite is unstable and thus does not participate in melting reactions. The water-saturated solidus, i.e., the first appearance of hydrous melt in the K–lherzolite composition, is located between 900 and 1000 °C at 4 GPa and between 1000 and 1100 °C at 5 and 6 GPa. Compositional jumps between hydrous melt and aqueous fluid at the solidus include a significant increase in the total dissolved solids load. All melts/fluids are peralkaline and calcium-rich. The melting reactions at the solidus are peritectic, as olivine, clinopyroxene, garnet, and H2O are consumed to generate hydrous melt plus orthopyroxene. Our fluid/melt compositional data demonstrate that the water-saturated hybrid peridotite solidus lies above 1000 °C at depths greater than 150 km and that the second critical endpoint is not reached at 6 GPa for a K2O–Na2O–CaO–FeO–MgO–Al2O3–SiO2–H2O–Cr2O3(–TiO2) peridotite composition.

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.