The water and fluid-mobile element cycles during serpentinite subduction. A review

The key role of serpentinites in the global cycles of volatiles, halogens and fluid-mobile elements in oceans and in subduction zones is now ascertained by many studies quantifying their element budgets and the composition of fluids they release during subduction. Geochemical tracers (e.g. B, As, Sb; stable B and radiogenic Sr and Pb isotopes) have also been employed to trace the provenance of serpentinites (slab or forearc mantle?) accreted to the plate interface of fossil subduction zones. In turn, this helps defining the tectonic processes, seismicity and mass transfer attending rock burial and exhumation within subduction zones. The results suggest that the sole use of geochemical data is insufficient to track the origin of subduction-zone serpentinites and the timing of serpentinization, whether oceanic or subduction-related. Integrated multidisciplinary studies of ophiolitic serpentinite show that pristine, oceanic, geochemical imprints (e.g. high 11B, marine Sr isotopes, low As + Sb) become reset towards more radiogenic Sr, lower 11B, and higher As + Sb via metasomatic exchange with crust-derived fluids during subduction accretion to the plate interface.The dehydration fluids released by serpentinite dehydration at various subduction stages and still preserved in these rocks as inclusions, carry significant amounts of halogens and fluid-mobile elements. The key compositional similarities of antigorite-breakdown fluids from different localities (Betic Cordillera, Spain; Central Alps, Switzerland) indicate that rocks record comparable subduction processes. We individuate the fluid-mediated exchange with sedimentary and/or crustal reservoirs during subduction as the key mechanism for geochemical hybridization of serpentinite. The antigorite dehydration fluids produced by hybrid serpentinites have high Cs, Rb, Ba, B, Pb, As, Sb and Li overlapping those of the arc lavas and representing the mixed serpentinite\u2013sediment (crustal) component released to arcs. This helps discriminating the mass transfer processes responsible for supra-subduction mantle metasomatism and arc magmatism. The studied plate-interface hybrid serpentinites are also proxies of forearc mantle metasomatized by slab fluids. Based on the above observations, we propose that the mass transfer from slabs to plate interface and/or forearc mantle and the subsequent down-drag of this altered mantle to subarc depths potentially is a major process operating in subduction zones.The nominally anhydrous olivine, orhopyroxene, clinopyroxene and garnet produced by serpentinite dehydration host appreciable amounts of halogens and fluid-mobile elements that can be recycled in the deep mantle beyond arcs. Involvement of de-serpentinized residues in lower mantle metasomatism begins to be increasingly recognized by studies of ocean island basalts (OIB) and of B-bearing blue diamonds and by the isotopic serpentinite compositions presented here

Petrology and Geochemistry of Serpentinites Associated with the Ultra-High Pressure Lago di Cignana Unit (ItalianWestern Alps)

In the Western Alps, the ophiolitic Zermatt–Saas Zone (ZSZ) and the Lago di Cignana Unit (LCU) record oceanic lithosphere subduction to high (540°C, 2·3GPa) and ultra-high pressure (600°C, 3·2GPa), respectively. The top of the Zermatt–Saas Zone in contact with the Lago di Cignana Unit consists of olivine þ Ti-clinohumite-bearing serpentinites (the Cignana serpentinite) hosting olivine þ Ti-clinohumite veins and dykelets of olivine þ Ti-chondrodite þ Ti-clinohumite. The composition of this serpentinite reveals a refertilized oceanic mantle peridotite protolith that became subsequently enriched in fluid-mobile elements (FME) during oceanic serpentinization. The olivine þ Ti-clinohumite veins in the Cignana serpentinite display Rare Earth Element (REE) and FME compositions quite similar to the host-rock, which suggests closed-system dehydration of this serpentinite during subduction. The Ti-chondrodite-bearing dykelets are richer in REE and FME than the host-rock and the dehydration olivine þ Ti-clinohumite veins: their Nd composition points to a mafic protolith, successively overprinted by oceanic metasomatism and by subduction zone recrystallization. These dykelets are comparable in composition to eclogites within the ultra-high pressure LCU that derive from subducted oceanic mafic crust. Different from the LCU, serpentinites from the core domains of the ZSZ display REE compositions indicating a depleted mantle protolith. The oceanic serpentinization of these rocks led to an increase in FME and to seawater-like Sr isotope compositions. The serpentinites sampled at increasing distance from the ultra-high pressure LCU reveal different mantle protoliths, still preserve an oceanic geochemical imprint and contain mafic dykelets affected by oceanic metasomatism. The subduction zone history of these rocks thus occurred under relatively closed system conditions, the only possible change during subduction being an enrichment in As and Sb recorded by the serpentinites closer to the crustal LCU. The ZSZ and Cignana serpentinites thus likely evolved in a slab setting and were weakly exposed to interaction with slab-derived fluids characteristic of plate interface settings. Our data suggest two possible scenarios for the evolution of the studied ZSZ and Cignana serpentinites. They are either part of a coherent ophiolite unit whose initial lithospheric mantle was variably affected by depletion and re-fertilization processes, or they belong to separate tectonic slices derived from two different oceanic mantle sections. In the Cignana serpentinite atop the ZSZ, the presence of Ti-chondrodite dykelets similar in composition to the LCU eclogites suggests these two domains were closely associated in the oceanic lithosphere and shared the same evolution to ultra-high pressure conditions during Alpine subduction

Commensurate Growth of Magnetite Microinclusions in Olivine under Mantle Conditions

Magnetite-bearing multiphase solid inclusions hosted in metamorphic olivine have been interpreted as final products of the trapping of the aqueous fluid produced by the subduction-zone dehydration of former serpentinites. We provide here a careful analysis performed by microfocus single-crystal X-ray diffraction of inclusions found in harzburgites from the Almirez Complex (Bétic Cordillera, Spain) to determine the occurrence of preferential crystallographic orientation relationships between the olivine host and the magnetite inclusion. The results demonstrate that the magnetite–olivine interface selectively displays parallelism between crystallographic planes (111) and (100) and between crystallographic directions ⟨110⟩ and ⟨011⟩, respectively. This evidence points to a clear epitaxial growth of magnetite on olivine. The calculation of the geometrical misfit between the two lattices in contact as a function of their relative azimuthal orientation shows that, under the aforementioned reciprocal orientation, a perfect commensurism is achieved; i.e., all of the nodes of the magnetite lattice coincide with nodes of the olivine lattice. This particular relationship must be interpreted as a unique occurrence, playing a fundamental role in favoring the heterogeneous nucleation of magnetite on olivine

Raman elastic geobarometry for anisotropic mineral inclusions

Elastic geobarometry for host-inclusion systems can provide new constraints to assess the pressure and temperature conditions attained during metamorphism. Current experimental approaches and theory are developed only for crystals immersed in a hydrostatic stress field, whereas inclusions experience deviatoric stress. We have developed a method to determine the strains in quartz inclusions from Raman spectroscopy using the concept of the phonon-mode Gr\ufcneisen tensor. We used ab initio Hartree-Fock/Density Functional Theory to calculate the wavenumbers of the Raman-active modes as a function of different strain conditions. Least-squares fits of the phonon-wavenumber shifts against strains have been used to obtain the components of the mode Gr\ufcneisen tensor of quartz (\u2060\u3b3m1 and \u3b3m3\u2060) that can be used to calculate the strains in inclusions directly from the measured Raman shifts. The concept is demonstrated with the example of a natural quartz inclusion in eclogitic garnet from Mir kimberlite and has been validated against direct X-ray diffraction measurement of the strains in the same inclusio

Raman elastic geobarometry for anisotropic mineral inclusions

Elastic geobarometry for host-inclusion systems can provide new constraints to assess the pressure and temperature conditions attained during metamorphism. Current experimental approaches and theory are developed only for crystals immersed in a hydrostatic stress field, whereas inclusions experience deviatoric stress. We have developed a method to determine the strains in quartz inclusions from Raman spectroscopy using the concept of the phonon-mode Grüneisen tensor. We used ab initio Hartree-Fock/Density Functional Theory to calculate the wavenumbers of the Raman-active modes as a function of different strain conditions. Least-squares fits of the phonon-wavenumber shifts against strains have been used to obtain the components of the mode Grüneisen tensor of quartz (γm1 and γm3) that can be used to calculate the strains in inclusions directly from the measured Raman shifts. The concept is demonstrated with the example of a natural quartz inclusion in eclogitic garnet from Mir kimberlite and has been validated against direct X-ray diffraction measurement of the strains in the same inclusion

Fluid‐mediated mass transfer between mafic and ultramafic rocks in subduction zones

© The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Codillo, E., Klein, F., Dragovic, B., Marschall, H., Baxter, E., Scambelluri, M., & Schwarzenbach‬, E. Fluid‐mediated mass transfer between mafic and ultramafic rocks in subduction zones. Geochemistry Geophysics Geosystems, 23, (2022): e2021GC010206, https://doi.org/10.1029/2021gc010206.Metasomatic reaction zones between mafic and ultramafic rocks exhumed from subduction zones provide a window into mass-transfer processes at high pressure. However, accurate interpretation of the rock record requires distinguishing high-pressure metasomatic processes from inherited oceanic signatures prior to subduction. We integrated constraints from bulk-rock geochemical compositions and petrophysical properties, mineral chemistry, and thermodynamic modeling to understand the formation of reaction zones between juxtaposed metagabbro and serpentinite as exemplified by the Voltri Massif (Ligurian Alps, Italy). Distinct zones of variably metasomatized metagabbro are dominated by chlorite, amphibole, clinopyroxene, epidote, rutile, ilmenite, and titanite between serpentinite and eclogitic metagabbro. Whereas the precursor serpentinite and oxide gabbro formed and were likely already in contact in an oceanic setting, the reaction zones formed by diffusional Mg-metasomatism between the two rocks from prograde to peak, to retrograde conditions in a subduction zone. Metasomatism of mafic rocks by Mg-rich fluids that previously equilibrated with serpentinite could be widespread along the subduction interface, within the subducted slab, and the mantle wedge. Furthermore, the models predict that talc formation by Si-metasomatism of serpentinite in subduction zones is limited by pressure-dependent increase in the silica activity buffered by the serpentine-talc equilibrium. Elevated activities of aqueous Ca and Al species would also favor the formation of chlorite and garnet. Accordingly, unusual conditions or processes would be required to stabilize abundant talc at high P-T conditions. Alternatively, a different set of mineral assemblages, such as serpentine- or chlorite-rich rocks, may be controlling the coupling-decoupling transition of the plate interface.M. Scambelluri acknowledges the Italian Ministry of Research MUR for granting the PRIN project n. 2017ZE49E7. This research was funded by NSF-OISE (Office of International Science & Engineering, Petrology & Geochemistry) PIRE, Award #1545903, and the WHOI Ocean Ventures Fund

Fossil subduction recorded by quartz from the coesite stability field

Metamorphic rocks are the records of plate tectonic processes whose reconstruction relies on correct estimates of the pressures and temperatures (P-T) experienced by these rocks through time. Unlike chemical geothermobarometry, elastic geobarometry does not rely on chemical equilibrium between minerals, so it has the potential to provide information on overstepping of reaction boundaries and to identify other examples of non-equilibrium behavior in rocks. Here we introduce a method that exploits the anisotropy in elastic properties of minerals to determine the unique P and T of entrapment from a single inclusion in a mineral host. We apply it to preserved quartz inclusions in garnet from eclogite xenoliths hosted in Yakutian kimberlites (Russia). Our results demonstrate that quartz trapped in garnet can be preserved when the rock reaches the stability field of coesite (the high-pressure and hightemperature polymorph of quartz) at 3 GPa and 850 \ub0C. This supports a metamorphic origin for these xenoliths and sheds light on the mechanisms of craton accretion from a subducted crustal protolith. Furthermore, we show that interpreting P and T conditions reached by a rock from the simple phase identification of key inclusion minerals can be misleading