Search

The oxidation state of sulfur in slab fluids is controversial, with both dominantly oxidized and reduced species proposed. Here we use in situ X-ray absorption spectroscopy analysis of sulfur-in-apatite to monitor changes in the oxidation state of sulfur during high-P metasomatism by slab fluids in the subduction channel. Our samples include a 73 cm continuous transect of reaction zones between a metagabbroic eclogite block and serpentinite matrix from a mélange zone on the island of Syros, Greece. The block core consists of garnet, omphacite, phengite, paragonite, epidote-clinozoisite, and rutile. In this region, apatite is only observed as elongate inclusions in omphacite cores. From the core outwards micas are increasingly replaced by epidote-clinozoisite, garnets are smaller and more frequent, pyrite + bornite is observed as inclusions in recrystallized omphacite, and apatite is increasingly abundant in the matrix and inclusions in garnet. A major transition at 48 cm separates an assemblage of Ca-Na amphibole, omphacite, chlorite, pyrite, and apatite from the inner garnet-bearing eclogite assemblages. Omphacite disappears from the assemblage at ~56 cm and amphibole compositions sharply transition to tremolite at 59 cm. Finally, the assemblage tremolite + talc + pyrite is observed after ~70 cm.Apatites in the eclogite assemblages exclusively display S6+ peaks in their absorption spectra. This includes apatite inclusions in omphacite in the least altered lithology, as well as matrix apatite and isolated apatite inclusions in garnet in the outermost metasomatized eclogite zone. In the intermediate pyrite-rich (~1–5 vol %) amphibole + omphacite + chlorite zone, apatite displays a strong S1- absorption peak in most grains, with rare analyses showing mixed S1- and S6+. Finally, apatite in the outermost tremolite-bearing assemblages only displays a S6+ peak. The pyrite-rich zone at 48 cm occurs at the initial interface between the serpentinite matrix and eclogite block, characterized by a dramatic decrease in Na content and Mg#. Our data suggest that reduction of S6+ in infiltrating fluids to S1- in pyrite became focused as Fe diffused across the steep Mg# gradient, resulting in pyrite precipitation. In contrast, S reduction in the Mg-rich tremolite-dominant portions of the transect was limited by a lack of Fe, resulting in low modes of pyrite and fluid buffered S6+ in apatite. Finally, S6+-bearing apatite is also observed in reaction zone lithologies from elsewhere on Syros, suggesting our observations are not isolated.Two important conclusions are drawn from these data and observations: (1) In the case of Syros, slab fluids at eclogite-facies conditions carried oxidized S6+, and (2) The interaction of these fluids with eclogites composed of ferrous-Fe silicates resulted in extensive sulfide precipitation.

MinPlot is a MATLAB®-based mineral formula recalculation and compositional plotting program for electron microprobe analyses (EPMA). The program offers recalculation and structural formula assignment for 15 different mineral groups: Garnet, pyroxene, olivine, amphibole, feldspar, mica, staurolite, cordierite, chlorite, chloritoid, talc, epidote, titanite, spinel, and sulfides. MinPlot is a fast and easy to use command line program and requires no prior computer programming knowledge. Percent mass fractions of oxides are loaded from datafiles and the user answers simple prompts to select mineral type, normalization scheme, and plotting options. Recalculated mineral formulas are automatically saved as output files and plots may be further manually customized by the user prior to saving. MinPlot can perform thousands of calculations in seconds and the modular nature of the program makes it simple to add new calculation routines in future releases. Finally, MinPlot utilizes simple matrix algebra and linear programing steps, and is thus anticipated to be forwards compatible with future releases of MATLAB®. Combined, these features make MinPlot a powerful and useful program for the processing of EPMA data.

Sulfur is a key element in the subduction zone-volcanic arc system; however, the mechanism(s) that recycle sulfur from the slab into the overlying volcanic arc are debated. Here we summarize recent advances in quantifying this component of the deep sulfur cycle. First, primary metamorphic or inherited sulfides in oceanic-type eclogites are only rarely observed as inclusions and are typically absent from the rock matrix. Additionally, sulfides are relatively common in rocks metasomatized at the slab-mantle interface by slab-derived fluids during exhumation. Combined, these two observations suggest that sulfur loss from subducted mafic crust is relatively efficient. Thermodynamic modeling in Perple_X using the Holland and Powell (2011) database combined with the Deep Earth Water model suggests that the efficiency and speciation of sulfur loss varies depending on the degree of seafloor alteration prior to subduction and the geothermal gradient of the slab. In relatively cold subduction zones, such as Honshu, slab-fluids derived from subducted mafic crust are predicted to exhibit elevated concentrations of HSO4-, SO42-, HSO3-, and CaSO4(aq), whereas hot subduction zones, such as Cascadia, are predicted to produce slab fluids enriched in HS- and H2S at lower pressures. The oxidation of sulfur expelled from subducted pyrite is balanced by the reduction of Fe3+ to Fe2+, consistent with the low Fe3+/SFe of exhumed eclogites relative to blueschists and altered oceanic crust. Where oxidized S-bearing fluids are produced, they are anticipated to interact with more reduced rocks at the slab-mantle interface and within the mantle wedge, resulting in sulfide precipitation and significant isotopic fractionation. The δ34S values of slab fluids are estimated to fall between -11 and +8 ‰. Rayleigh fractionation during progressive fluid-rock interaction results in fractionations of tens of per mil as oxidized species are depleted and sulfides are precipitated, resulting in δ34S values of sulfides that easily span the -21.7 to +13.9 ‰ range observed in metasomatic sulfides in exhumed high-pressure rocks. However, in subduction zones where reduced species prevail, the S isotopic signature of slab fluids is expected to reflect their source and will exhibit a narrower range in δ34S values. As a result, the δ34S values measured in arc magmas may not always be a reliable indicator of the contribution of different components of the slab, such as sediments vs. AOC. Additionally, the impact of S recycling on the oxygen fugacity of arc magmas is expected to vary both spatially and temporally throughout Earth history.