Highly oxidising fluids generated during serpentinite breakdown in subduction zones.

Subduction zones facilitate chemical exchanges between Earth's deep interior and volcanism that affects habitability of the surface environment. Lavas erupted at subduction zones are oxidized and release volatile species. These features may reflect a modification of the oxidation state of the sub-arc mantle by hydrous, oxidizing sulfate and/or carbonate-bearing fluids derived from subducting slabs. But the reason that the fluids are oxidizing has been unclear. Here we use theoretical chemical mass transfer calculations to predict the redox state of fluids generated during serpentinite dehydration. Specifically, the breakdown of antigorite to olivine, enstatite, and chlorite generates fluids with high oxygen fugacities, close to the hematite-magnetite buffer, that can contain significant amounts of sulfate. The migration of these fluids from the slab to the mantle wedge could therefore provide the oxidized source for the genesis of primary arc magmas that release gases to the atmosphere during volcanism. Our results also show that the evolution of oxygen fugacity in serpentinite during subduction is sensitive to the amount of sulfides and potentially metal alloys in bulk rock, possibly producing redox heterogeneities in subducting slabs.

Introduction

During subduction, the increase of pressure and temperature conditions in the subducting plate results in hydrous mineral breakdown and the release of volatile-rich fluids. The oxidized or reduced nature of the released fluids is controlled by mineral local equilibrium changes and can be monitored through the oxygen fugacity (fO2). Although large variabilities of fO2 were shown to occur in subduction zones, especially in melange zones constituting the plate/mantle interface[1], previous petrological and geochemical studies have shown that the redox state of the subducted mafic crust does not change during prograde metamorphism and that its redox budget remains relatively constant[2]. Therefore, other sources of fluids must be considered to explain the oxidized nature of arc magmas[3-5].

Serpentinized mantle peridotites can constitute a major source of water at depth in subduction zones[6]. Those rocks comprise a significant part of subducting oceanic lithosphere hydrated near the seafloor at slow or ultra-slow spreading ridges[7, 8] or slab-bending related faults[9], and are also present as part of the slab-wedge interface or mantle wedge that is percolated by aqueous fluids emanating from the slab during subduction[10, 11]. Serpentinites are thus ubiquitous in subduction zones. Furthermore, several studies have suggested that the breakdown of these rocks is the most likely source of fluid-mobile elements (e.g. B), volatiles and halogens enriched in arc magmas (e.g. [12-17]). The dehydration of serpentinized mantle results in a net decrease in bulk rock volatile concentrations[6, 16, 17] and Fe3+/∑e[18, 19]. In addition, observations of hematite-magnetite assemblages in dehydrated serpentinites[20, 21] suggest a high fO2, from one to five log units above the Quartz-Fayalite-Magnetite (QFM) oxygen buffer, during antigorite breakdown in subduction zones which could be compatible with the release of oxidized S and C in slab-derived fluids. If correct, this model would have major implications for the cycling of volatiles between the Earth’s surface and interior, volcanic degassing (e.g. SOX, COX species emitted from volcanoes), and the occurrence of metal sulfide ore deposits[22, 23]. Although several studies have indicated some evidence of variations of fO2 in antigorite-bearing serpentinites that can be related to the extent of initial serpentinization[24, 25], there is at this time no consensus on the evolution of fO2 or the redox state of serpentinite-derived fluids during antigorite breakdown in subduction zones.

Recent advances in theoretical and experimental aqueous geochemistry resulting in the Deep Earth Water (DEW) model now enable the calculation of equilibrium constants involving minerals and aqueous ions, metal complexes, and organics to 6.0 GPa and 1,200 °C[26-33]. In the present study, we model the dehydration of antigorite-bearing serpentinite to investigate its potential to liberate highly oxidizing fluids in subduction zones. The model uses a conceptual scenario in which an initial fluid chemistry was set by reaction of water with an antigorite-bearing serpentinite at 630 °C before undergoing an increase of temperature to 660 °C at constant pressure. The initial model assemblage is composed of antigorite, olivine (XMg [Mg/(Fe + Mg)] = 0.885), clinochlore, magnetite, and tremolite in agreement with field observations from Padron-Navarta et al.[34]. This assemblage sets the initial fO2 which ranges from near QFM at 500 °C and 2.0 GPa to QFM + 4.2 at 650 °C and 2.0 GPa (Supp. Info Fig. S1). At the lowest temperature in this range the fO2 near QFM is in agreement with recent work on natural samples[24]. However, at 630 °C the fO2 is several units above QFM just before the breakdown of antigorite. In this study, we explore the further dramatic fO2 changes on heating through the breakdown of antigorite.

Modelling serpentinite dehydration during subduction

Figure 1 displays the evolution of the fO2, mineral assemblages, and fluid composition during antigorite breakdown at 2 GPa (modelling carried out at 4 GPa displays a similar evolution as presented in Supp. Info Fig. S2). In agreement with natural[21, 34] observations, the reaction progress leads to the successive growth of olivine, clinochlore, and orthopyroxene and the progressive decrease of the amount of magnetite (Fig. 1b,d). During the first part of the reaction path (logξ < −4), the crystallization of olivine and chlorite is accompanied by a progressive increases of fO2 up to four log units above QFM oxygen buffer (Fig. 1a). At these conditions, hematite becomes in equilibrium with magnetite and buffers the fO2. The reaction progress is then accompanied with a progressive decrease of the amount of magnetite and the appearance of orthopyroxene (+/− tremolite) in the rock (Fig. 1d). A second increase of fO2 appears later (logξ > −1.4) and corresponds to the disappearance of magnetite (Fig. 1a,b).

Figure 1

Predicted evolution of fO2, mineral reactants, aqueous phase composition, and mineral products during antigorite breakdown at 2 GPa in a sulfur-free model. (a) Evolution of temperature and fO2 in the system. QFM and magnetite-hematite buffers are reported in dashed lines. (b) Abundances of key minerals involved.Hematite appears at the beginning of the reaction and remains stable until the disappearance of antigorite and magnetite. (c) Concentrations of major species in the aqueous fluid. (d) Minerals produced during antigorite breakdown. The x axis represents the logarithm of the reaction progress variable ξ, which is equal to the number of moles of each reactant mineral destroyed during the reaction progress.

Few variations in pH and the composition of the released fluids (Fig. 1c) are observed during the reaction progress. The pH is alkaline and varies slightly around 5 (neutral pH is about 3.4 at these conditions; Supp. Info Fig. S3), while the aqueous species are dominated by Na+ (aq), Cl− (aq) and Si(OH)4(aq) and SiO(OH)3 − (aq), with minor amounts of NaHCO3(aq), Al(OH)4 − (aq), Ca(OH)+ (aq), Mg(OH)+ (aq), and FeCl2(aq). The presence of equilibrated hematite-magnetite assemblages in our models is compatible with recent observations in meta-peridotites from the Cerro Del Almirez massif[21] or partly dehydrated serpentinites from Western Alps meta-ophiolites[20]. It therefore suggests that the dehydration of serpentinites in subduction zones can take place at high fO2, close to the hematite-magnetite buffer. At such P-T-fO2 -pH conditions, volatile and redox-sensitive elements, such as sulfur or carbon, are expected to be mobilized under in their oxidized form (e.g. CO2 0 aq) or SO4 2− (aq)) rather than reduced species (e.g. CH4(aq) or HS− (aq); Fig. 1a).

The importance of sulfides

Although the progressive increase in fO2 during prograde metamorphism is suggestive of oxidized fluids release during antigorite breakdown, previous studies have shown that the presence of Ni-Fe alloys and/or sulfides can buffer the fO2 down to four or five log units below the QFM buffer in serpentinites[35, 36]. In addition, several studies have shown that sulfides can persist[17] or even be formed[24, 25] during serpentinite devolatilization. In order to test the impact of sulfur-bearing phases during the progressive devolatilization of serpentinites in subduction zones, we ran a suite of models containing various amounts of pyrrhotite (0.001, 0.01 and 0.1 moles), which is commonly observed in abyssal and orogenic serpentinites[17]. The most pronounced difference between sulfide-free and sulfide-bearing models is the evolution of fO2 during the reaction progress (Fig. 2a). The presence of a relatively small amount of pyrrhotite (0.001 to 0.01 moles of pyrrhotite) is accompanied with a decrease of fO2 from three to two log units above QFM buffer (logξ < −3 or −2, Fig. 2a). This stage is associated with the total destruction of the pyrrhotite and an increase of sulfur concentration in the serpentinite-derived fluids. Interestingly, even if hematite is absent from the system at this stage, sulfate bearing species (HSO4 − (aq), SO4 2− (aq), CaSO4 (aq), MgSO4 (aq)) are dominating the fluid composition relative to sulfide bearing species (H2S and HS−, Fig. 2b). This suggests that even at modestly oxidizing fO2 values, sulfates are continuously released in serpentinite-derived fluids. As the reaction progresses (from logξ of −3 to −1.4), the destruction of sulfide and then magnetite is accompanied by an increase of fO2 and the appearance of hematite (Fig. 2b). In these conditions, sulfide species are at negligible concentrations in the fluid phase which is mainly composed of HSO4 − (aq), SO4 2− (aq), CaSO4 (aq) and MgSO4 (aq). For high sulfide concentrations (0.1 mole of pyrrhotite), aqueous sulphide species remain at significant concentrations in the fluid during the whole reaction progress and become more abundant than sulphate species above logξ of about −1.5 (Fig. 2c). These species buffer the fO2 at three units above QFM buffer but below the Hematite-Magnetite buffer, therefore preventing hematite crystallization (Fig. 2a). However, it should be noted that, even if the presence of these high amounts of pyrrhotite buffer the fO2 at relatively low values, the amount of sulfur and sulfate dissolved in the fluid is significantly higher in those models (Fig. 2b and c).

Figure 2

Predicted evolution of fO2 and aqueous phase composition in the presence of pyrrhotite during antigorite breakdown at 2 GPa. (a) Evolution of fO2 for different amounts of pyrrhotite (No pyrrhotite, 0.001, 0.01 and 0.1 moles). (b,c) Number of moles of sulfur dissolved into the fluid per kg of water with the evolution of major aqueous sulfur species. In b and c, the amount of sulfur present in the initial assemblage varies from 0.001 to 0.1 moles of pyrrhotite.

Evolution of oxygen fugacity and implications for the nature of slab derived fluids

Iron, carbon and sulfur are the main redox sensitive elements that can be transferred by fluids from the slab to the mantle wedge, having thus the potential to modify the redox state of the source region of arc magmas above subduction zones[4, 5, 37]. In a previous Fe isotope study on subducted meta-ophiolites, it was shown that Fe can be mobile in slab-derived fluids in the form of Fe(II)-Cl2 and/or Fe(II)-SOX [20]. In good agreement with these observations, thermodynamic calculations predict that during serpentinite dehydration, the magnetite modal amount significantly decreases (Fig. 1a,b) while the amount of Fe release in the fluid remains very low (Fig. 1c) and is dominated by the FeCl2 species. We note that the complexation of Fe(II) to sulfate is not considered in the DEW model and even though Fe(III)-complexes are included, they are completely unimportant at the relatively high temperatures and low chloride concentrations of the present study. In other words, with our present state of knowledge of aqueous Fe-complexes, there are none known that could contribute significantly to the transport of Fe(III). It should also be noted that chlorine concentrations in fluid inclusions hosted by metamorphic olivines range up to 2 wt% NaCl. Such chlorine concentrations could potentially lead to a higher mobility of Fe(II) in serpentinites-derived fluids[14, 15].

Carbon is mainly transported via CO2,aq under the conditions we investigated in subduction zones. Sources of this C could be in sediments and ophicarbonates (up to 3 wt%; e.g. refs 38 and 39). Although there are few constraints on carbon mobility during serpentinite dehydration, eclogitic serpentinites and meta-peridotites from Almirez massif display similar low concentrations (between 200 and 800 ppm) and preserve δ13C and C consistent with mixing at low pressure between reduced carbon and seawater derived carbonates[40]. It thus suggests that carbon is not the main oxidizing agent release by slab serpentinite dehydration at high pressure. Instead, the sulfur concentrations of subducted ultramafic rocks are highly variable (from about 50 ppm to 2000 ppm) reflecting a mobility of sulfur in serpentinite-derived fluids during subduction[40, 41]. Indeed the presence of sulfur-bearing daughter phases precipitating after fluid entrapment in dehydrated serpentinite shows a significant mobility of sulfur in serpentinite derived fluids during subduction[15]. In agreement with these studies, sulfur-bearing models reveal the total destruction of sulfur-bearing phases during the first stages of serpentinite dehydration and the presence of sulfate in serpentinite-derived fluids. The ratio of total sulfate over total sulfur in serpentinite derived fluids increases with the fO2 (Fig. 2b and c).

The fO2 evolution of serpentinites during antigorite breakdown is clearly highly influenced by pre-subduction characteristics of the rocks, specifically the initial assemblage (e.g. modal amount of sulfide, Fig. 2a). At mid-oceanic ridges, serpentinites formed in the deep part of the lithosphere or nearby metagabbroic intrusions are likely to contain a high amount of mantle sulfides[17, 42] and other reduced species such as awaruite[25]. In contrast, serpentinites sampled at the seafloor level are formed at high water-rock ratio and can contain a large amount of sulfate[14, 42]. Such large-scale heterogeneities are likely to be preserved during subduction[43] resulting in a variety of mineralogical and fO2 evolution pathways in ultramafic rocks during oceanic lithosphere dehydration. Indeed, the existence of subducted serpentinites equilibrated at variable fO2, from −3 to + 5[21, 25], in different subduction settings is in good agreement with this scenario. It thus suggests that in the deep and/or less serpentinized parts of the subducted lithosphere, the presence of reduced assemblages (e.g. awaruite, sulfide) can buffer the fO2 to relatively low values, until their eventual destruction during antigorite breakdown, while toward to the top of the slab and/or nearby slab-bending related fractures the ultramafic rocks will preferentially crystallize assemblages equilibrated at high fO2. As a consequence, the fO2 of the slab is likely to be heterogeneous, reflecting different pre-subduction mineralogy. However, in highly serpentinized peridotites, regardless of the considered models, sulfur seems to be highly mobile in fluids during antigorite breakdown (Fig. 2b,c). Interestingly, although sulfide can buffer the fO2 to relatively modest elevations above QFM, sulfate-bearing species (HSO4 − (aq), SO4 2− (aq), CaSO4 (aq), MgSO4 (aq)) always represent a significant proportion of the sulfur-bearing phases carried by the fluid. The migration of these fluids from the slab to the slab-mantle interface or mantle wedge can therefore enhance the oxidation of the mantle wedge Fe2+ to Fe3+ in response to the reduction of slab fluid SO4 2− to S2−. Such processes may provide an oxidized mantle source region for primary arc magmas. In this way, the origin of the oxidized signature carried by fluids from the slab may be an inevitable consequence of the breakdown of antigorite during subduction.

Supplementary Figures