skip to main content

Title: Isotopic Compositions of Sulfides in Exhumed High‐Pressure Terranes: Implications for Sulfur Cycling in Subduction Zones

Subduction is a key component of Earth's long‐term sulfur cycle; however, the mechanisms that drive sulfur from subducting slabs remain elusive. Isotopes are a sensitive indicator of the speciation of sulfur in fluids, sulfide dissolution‐precipitation reactions, and inferring fluid sources. To investigate these processes, we report δ34S values determined by secondary ion mass spectroscopy in sulfides from a global suite of exhumed high‐pressure rocks. Sulfides are classified into two petrogenetic groups: (1) metamorphic, which represent closed‐system (re)crystallization from protolith‐inherited sulfur, and (2) metasomatic, which formed during open system processes, such as an influx of oxidized sulfur. The δ34S values for metamorphic sulfides tend to reflect their precursor compositions: −4.3 ‰ to +13.5 ‰ for metabasic rocks, and −32.4 ‰ to −11.0 ‰ for metasediments. Metasomatic sulfides exhibit a range of δ34S from −21.7 ‰ to +13.9 ‰. We suggest that sluggish sulfur self‐diffusion prevents isotopic fractionation during sulfide breakdown and that slab fluids inherit the isotopic composition of their source. We estimate a composition of −11 ‰ to +8 ‰ for slab fluids, a significantly smaller range than observed for metasomatic sulfides. Large fractionations during metasomatic sulfide precipitation from sulfate‐bearing fluids, and an evolving fluid composition during reactive transport may account for the entire ~36 ‰ range of metasomatic sulfide compositions. Thus, we suggest that sulfates are likely the dominant sulfur species in slab‐derived fluids.

more » « less
Award ID(s):
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geochemistry, Geophysics, Geosystems
Page Range / eLocation ID:
p. 3347-3374
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Outcrops of brittle faults are rare in Marie Byrd Land, West Antarctica, because fault damage zones commonly undergo enhanced erosion and form bedrock troughs occupied by glacier ice. Where exposures do exist, faults yield information about regional strain in the West Antarctic Rift System (WARS) and may host minerals that contain a record of the temperature and chemistry of fluids during regional-scale faulting. In MBL’s southern Ford Ranges, bordering Ross Sea, a distinctive fault array was sampled that hosts tourmaline and quartz, a mineral-pair that can provide temperature and composition of fault-associated fluids, using 18O. Host rocks are tourmaline-free. At three separate sites, fault surfaces display strongly aligned tourmaline, suggesting that mineralization occurred during tectonism. One site features highly polished, or mirrored, surfaces, a characteristic that may indicate tourmaline precipitation during seismic slip. The orientation and kinematics of the high angle faults are NNW-striking: normal-slip, and WNW-ESE striking: right-lateral strike-slip. The timing of mineralization is yet to be determined, but viable possibilities are that the faults formed during broad intracontinental extension during formation of Ross Embayment in the Cretaceous, or during development of deep, narrow basins beneath the RIS grounding zone, in the Neogene (newly detected, see Tankersley et al., this meeting). Once formed, tourmaline is resistant to chemical and isotopic re-equilibration, and therefore can retain a record of its conditions during formation. We used oxygen isotope compositions of tourmaline and quartz pairs to investigate temperatures, fluid-rock ratios, and fluid sources, with bearing on fault-localized flux of fluids and geothermal heat. Analyzed tourmaline and quartz were separated from the upper ~2mm of the fault surfaces, as well as quartz separated from host rock in the same hand samples. Tourmaline 18O ratios (n=4) fall within a range of +9.2 to +10.4 ± 0.1 ‰ VSMOW (average 9.7‰, StDev = 0.7). Paired quartz yield 18O values of +11.1 to +10.3 ± 0.1 ‰. Relative isotopic homogeneity between sites suggests similar fluid conditions were present across the region and supports field evidence for that the structures form a regional fault array. ∆Qtz-Trm values fall between 1.3 and 2.0, and 18O of quartz in faults closely resembles 18O of host rock quartz. We tentatively determine the water oxygen isotope ratio as greater than ~7.7 ‰. Plutonic-metamorphic associations in the immediate region, and comparisons with similar faults elsewhere (i.e. Isola d’Elba, Italy), suggest temperatures as high as 500°C for the fluids that circulated into the faults. The data are interpreted to show that brittle faults provided pathways for hot fluids derived from mid-crustal processes to make their way to shallow crustal depths. 18O values indicate magmatic and/or metamorphic fluid sources, with minor to no introduction of meteoric fluids. Tourmaline-quartz pairs did not attain equilibrium, likely due to tourmaline’s rapid crystallization. On-going investigation includes analysis of H and B isotopes in tourmaline, which will better characterize the relationship between fault-hosted and mid-crustal fluids. 
    more » « less
  2. The sulfur isotope composition of volcanic rocks in arcs can be difficult to constrain because significant fractionation can occur during degassing. Mafic and ultramafic cumulates represent the least degassed part of the magmatic arc system, thereby offering an opportunity to investigate undegassed sulfur in arcs. Recent work on high pressure metamorphic rocks has suggested that subducted materials can retain their original isotopic composition to sub-arc depths. In particular, extreme negative δ34S values can be retained in subducted sediments. The purpose of this project is to investigate to what extent these deep subduction zone processes are reflected in the sulfur isotope signature of arc magmas. In the Lesser Antilles arc, there is a gradual decrease in terrigenous sediment being subducted from south to north. An estimated ~15% subducted sediment in the south and ~2% in the north is reflected in the chemical and isotopic composition of the Lesser Antilles arc magmas. Sulfides in these magma- derived cumulates record the earliest stages of magma evolution and are a more faithful monitor of the sulfur isotopic composition of the magma source region in the mantle than erupting lavas. We hypothesize that the decrease in terrigenous sediment being subducted from the south to north will be reflected in the S isotopes in cumulate samples. Samples of mafic and ultramafic cumulates have been collected from fourteen islands across the Lesser Antilles arc. Primary rock types are olivine gabbro, amphibole gabbro, plagioclase gabbro, and olivine gabbronorite. Sulfide minerals include pyrite, chalcopyrite, and pyrrhotite, and typically occur as spherical blebs. Sulfides are found primarily as inclusions in clinopyroxene, amphibole, olivine, and plagioclase. Sulfides occur less frequently as inclusions in magnetite and within the matrix. Analyses of sulfur isotopes in cumulate sulfides are currently underway. The decrease in the amount sediment being subducted from south to north in the Lesser Antilles arc should result in δ34S values that increase from south to north (more sediment subducted = more negative δ34S values). 
    more » « less
  3. Abstract

    Basalts from the Samoan volcanoes sample contributions from all of the classical mantle endmembers, including extreme EM II and high3He/4He components, as well as dilute contributions from the HIMU, EM I, and DM components. Here, we present multiple sulfur isotope data on sulfide extracted from subaerial and submarine whole rocks (N = 16) associated with several Samoan volcanoes—Vailulu‘u, Malumalu, Malutut, Upolu, Savai‘i, and Tutuila—that sample the full range of geochemical heterogeneity at Samoa and upon exhaustive compilation of S‐isotope data for Samoan lavas, allow for an assessment of the S‐isotope compositions associated with the different mantle components sampled by the Samoan hotspot. We observe variable S concentrations (10–1,000 ppm) and δ34S values (−0.29‰ ± 0.30 to +4.84‰ ± 0.30, 2σ). The observed variable S concentrations are likely due to sulfide segregation and degassing processes. The range in δ34S reflects mixing between the mantle origin and recycled components, and isotope fractionations associated with degassing. The majority of samples reveal Δ33S within uncertainty of Δ33S = 0‰ ± 0.008. Important exceptions to this observation include: (a) a negative Δ33S (−0.018‰ ± 0.008, 2σ) from a rejuvenated basalt on Upolu island (associated with a diluted EM I component) and (b) previously documented small (but resolvable) Δ33S values (up to +0.027 ± 0.016) associated with the Vai Trend (associated with a diluted HIMU component). The variability we observed in Δ33S is interpreted to reflect contributions of sulfur of different origins and likely multiple crustal protoliths. Δ36S versus Δ33S relationships suggest all recycled S is of post‐Archean origin.

    more » « less
  4. Abstract

    Sulfur belongs among H2O, CO2, and Cl as one of the key volatiles in Earth’s chemical cycles. High oxygen fugacity, sulfur concentration, and δ34S values in volcanic arc rocks have been attributed to significant sulfate addition by slab fluids. However, sulfur speciation, flux, and isotope composition in slab-dehydrated fluids remain unclear. Here, we use high-pressure rocks and enclosed veins to provide direct constraints on subduction zone sulfur recycling for a typical oceanic lithosphere. Textural and thermodynamic evidence indicates the predominance of reduced sulfur species in slab fluids; those derived from metasediments, altered oceanic crust, and serpentinite have δ34S values of approximately −8‰, −1‰, and +8‰, respectively. Mass-balance calculations demonstrate that 6.4% (up to 20% maximum) of total subducted sulfur is released between 30–230 km depth, and the predominant sulfur loss takes place at 70–100 km with a net δ34S composition of −2.5 ± 3‰. We conclude that modest slab-to-wedge sulfur transport occurs, but that slab-derived fluids provide negligible sulfate to oxidize the sub-arc mantle and cannot deliver34S-enriched sulfur to produce the positive δ34S signature in arc settings. Most sulfur has negative δ34S and is subducted into the deep mantle, which could cause a long-term increase in the δ34S of Earth surface reservoirs.

    more » « less
  5. null (Ed.)
    Geologic processes at convergent plate margins control geochemical cycling, seismicity, and deep biosphere activity in subduction zones and suprasubduction zone lithosphere. International Ocean Discovery Program (IODP) Expedition 366 was designed to address the nature of these processes in the shallow to intermediate depth of the Mariana subduction channel. Although no technology is available to permit direct sampling of the subduction channel of an intraoceanic convergent margin at depths up to 18 km, the Mariana forearc region (between the trench and the active volcanic arc) provides a means to access this zone. Active conduits, resulting from fractures in the forearc, are prompted by along- and across-strike extension that allows slab-derived fluids and materials to ascend to the seafloor along associated faults, resulting in the formation of serpentinite mud volcanoes. Serpentinite mud volcanoes of the Mariana forearc are the largest mud volcanoes on Earth. Their positions adjacent to or atop fault scarps on the forearc are likely related to the regional extension and vertical tectonic deformation in the forearc. Serpentinite mudflows at these volcanoes include serpentinized forearc mantle clasts, crustal and subducted Pacific plate materials, a matrix of serpentinite muds, and deep-sourced formation fluid. Mud volcanism on the Mariana forearc occurs within 100 km of the trench, representing a range of depths and temperatures to the downgoing plate and the subduction channel. These processes have likely been active for tens of millions of years at this site and for billions of years on Earth. At least 10 active serpentinite mud volcanoes have been located in the Mariana forearc. Two of these mud volcanoes are Conical and South Chamorro Seamounts, which are the furthest from the Mariana Trench at 86 and 78 km, respectively. Both seamounts were cored during Ocean Drilling Program (ODP) Legs 125 and 195, respectively. Data from these two seamounts represent deeper, warmer examples of the continuum of slab-derived materials as the Pacific plate subducts, providing a snapshot of how slab subduction affects fluid release, the composition of ascending fluids, mantle hydration, and the metamorphic paragenesis of subducted oceanic lithosphere. Data from the study of these two mud volcanoes constrain the pressure, temperature, and composition of fluids and materials within the subduction channel at depths of about 18 to 19 km. Understanding such processes is necessary for elucidating factors that control seismicity in convergent margins, tectonic and magma genesis processes in the forearc and volcanic arc, fluid and material fluxes, and the nature and variability of environmental conditions that impact subseafloor microbial communities. Expedition 366 centered on data collection from cores recovered from three serpentinite mud volcanoes that define a continuum of subduction-channel processes defined by the two previously cored serpentinite mud volcanoes and the trench. Three serpentinite mud volcanoes (Yinazao, Fantangisña, and Asùt Tesoro) were chosen at distances 55 to 72 km from the Mariana Trench. Cores were recovered from active sites of eruption on their summit regions and on the flanks where ancient flows are overlain by more recent ones. Recovered materials show the effects of dynamic processes that are active at these sites, bringing a range of materials to the seafloor, including materials from the lithosphere of the Pacific plate and from subducted seamounts (including corals). Most of the recovered material consists of serpentinite mud containing lithic clasts, which are derived from the underlying forearc crust and mantle and the subducting Pacific plate. Cores from each of the three seamounts drilled during Expedition 366, as well as those from Legs 125 and 195, include material from the underlying Pacific plate. A thin cover of pelagic sediment was recovered at many Expedition 366 sites, and at Site U1498 we cored through serpentinite flows to the underlying pelagic sediment and volcanic ash deposits. Recovered serpentinites are largely uniform in major element composition, with serpentinized ultramafic rocks and serpentinite muds spanning a limited range in SiO2 , MgO, and Fe2 O3 compositions. However, variation in trace element composition reflects pore fluid composition, which differs as a function of the temperature and pressure of the underlying subduction channel. Dissolved gases H2 , CH4 , and C2 H6 are highest at the site furthest from the trench, which also has the most active fluid discharge of the Expedition 366 serpentinite mud volcanoes. These dissolved gases and their active discharge from depth likely support active microbial communities, which were the focus of in-depth subsampling and preservation for shore-based analytical and culturing procedures. The effects of fluid discharge were also registered in the porosity and GRA density data indicated by higher than expected values at some of the summit sites. These higher values are consistent with overpressured fluids that minimize compaction of serpentinite mud deposits. In contrast, flank sites have significantly greater decreases in porosity with depth, suggesting that processes in addition to compaction are required to achieve the observed data. Thermal measurements reveal higher heat flow values on the flanks (~31 mW/m2) than on the summits (~17 mW/m2) of the seamounts. The new 2G Enterprises superconducting rock magnetometer (liquid helium free) revealed relatively high values of both magnetization and bulk magnetic susceptibility of discrete samples related to ultramafic rocks, particularly in dunite. Magnetite, a product of serpentinization, and authigenic carbonates were observed in the mudflow matrix materials. In addition to coring operations, Expedition 366 focused on the deployment and remediation of borehole casings for future observatories and set the framework for in situ experimentation. Borehole work commenced at South Chamorro Seamount, where the original-style CORK was partially removed. Work then continued at each of the three summit sites following coring operations. Cased boreholes with at least three joints of screened casing were deployed, and a plug of cement was placed at the bottom of each hole. Water samples were collected from two of the three boreholes, revealing significant inputs of formation fluids. This suggests that each of the boreholes tapped a hydrologic zone, making these boreholes suitable for experimentation with the future deployment of a CORK-lite. An active education and outreach program connected with many classrooms on shore and with the general public through social media. 
    more » « less