skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Sulfur_X: A Model of Sulfur Degassing During Magma Ascent
Key Points Sulfur_X is a new open‐source magma degassing model that accurately predicts the volatile and redox evolution of ascending arc magmas Sulfur_X shows that sulfur can start degassing in the lower crust or near‐surface, depending on the initial S 6+ /ΣS and H 2 O in the melt The vapor compositions predicted by Sulfur_X can be used to interpret the CO 2 /S T ratios in high‐ T volcanic gases, an eruption precursor  more » « less
Award ID(s):
1933773 2017814
PAR ID:
10459039
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Geochemistry, Geophysics, Geosystems
Volume:
24
Issue:
4
ISSN:
1525-2027
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The mainly deep-submarine Ontong Java Plateau (OJP) is the result of the largest outpouring of lava in the geologic record. Volcanic events of this magnitude can have dramatic environmental impacts due to volatile emissions. We report new S measurements in naturally glassy, olivine-hosted melt inclusions and pillow basalt glasses from the OJP. We combined these data with previous S measurements in OJP glasses to quantify S degassing in a suite of OJP glasses. Comparison with an S degassing model suggests OJP lavas that erupted at depths ~>1500 m did not degas S; OJP lavas that erupted at depths ~<1500 m degassed up to ~40% initial S, but these lavas likely made up a small fraction of OJP lavas. This result suggests that despite its large volume compared to continental large igneous provinces (LIPs), OJP lavas emitted less S, potentially contributing to its muted environmental impact. The OJP may provide a framework for the temporal evolution of S degassing at oceanic LIPs, with early eruptions at great water depths releasing limited to no S, and later eruptions at shallow water depths releasing larger, but still limited amounts of S. This framework may also have implications for continental LIP magmas, which may release significant amounts of CO2 but limited amounts of S during intrusive activity, with magmatic S emissions only becoming important during extrusive phases. 
    more » « less
  2. null (Ed.)
    During the formation of terrestrial planets, volatile loss may occur through nebular processing, planetesimal differentiation, and planetary accretion. We investigate iron meteorites as an archive of volatile loss during planetesimal processing. The carbon contents of the parent bodies of magmatic iron meteorites are reconstructed by thermodynamic modeling. Calculated solid/molten alloy partitioning of C increases greatly with liquid S concentration, and inferred parent body C concentrations range from 0.0004 to 0.11 wt%. Parent bodies fall into two compositional clusters characterized by cores with medium and low C/S. Both of these require significant planetesimal degassing, as metamorphic devolatilization on chondrite-like precursors is insufficient to account for their C depletions. Planetesimal core formation models, ranging from closed-system extraction to degassing of a wholly molten body, show that significant open-system silicate melting and volatile loss are required to match medium and low C/S parent body core compositions. Greater depletion in C relative to S is the hallmark of silicate degassing, indicating that parent body core compositions record processes that affect composite silicate/iron planetesimals. Degassing of bare cores stripped of their silicate mantles would deplete S with negligible C loss and could not account for inferred parent body core compositions. Devolatilization during small-body differentiation is thus a key process in shaping the volatile inventory of terrestrial planets derived from planetesimals and planetary embryos. 
    more » « less
  3. Porphyry Cu ± Mo ± Au and iron oxide-apatite (IOA) deposits rarely occur in spatial and temporal proximity in Phanerozoic arc-related settings, and the formation of these mineral deposit types in an evolving arc setting remains poorly understood. Specifically, the roles of magma composition and the tectonic regime remain the subject of some debate. Here, we systematically estimated the P-T-fO2 conditions and H2O-S-Cl contents for dioritic to granodioritic source magmas for porphyry and skarn Cu ± Au (150–135 Ma) and IOA deposits (~130 Ma) that formed in transpressional and transtensional settings in the Middle-Lower Yangtze River metallogenic belt, China. Our estimates show that, compared to IOA deposits, the porphyry- and skarn-related magmas were relatively felsic, cooler, and more hydrous. These geochemical features are consistent with the tectonic transition from subduction to slab rollback of the paleo-Pacific plate in the East Asia continental margin at <135 Ma and concomitant crustal extension and steepening of the regional geothermal gradient. Apatite data reveal that the silicate melts associated with the porphyry and skarn Cu ± Au and IOA deposits had comparable predegassed S concentrations (~0.13 ± 0.06 wt % vs. ~0.16 ± 0.09 wt % on average), but that IOA-related melts contained higher predegassed Cl/H2O ratios (~0.11 ± 0.03 vs. ~0.04 ± 0.03 for porphyry- and skarn-related magmas) that decreased by one order of magnitude after magmatic degassing. Magmatic fO2 estimated using zircon and amphibole, reported in log units relative to the fayalite-magnetite-quartz (FMQ) redox buffer, gradually increased during cooling of the porphyry- and skarn-related magmas (ΔFMQ +0.7 to +2.5) at 950° to 800°C and decreased to ΔFMQ +1 at 700°C owing to fractionation of Fe2+-rich minerals and subsequent S degassing, respectively. In contrast, the magmatic fO2 values for the IOA-related source magmas varied significantly from ΔFMQ –1.5 to ΔFMQ +2.5 but generally show an increasing trend with cooling from 970° to 700°C that probably resulted from variable degrees of evaporite assimilation, fractionation of Fe2+-rich minerals, and Cl degassing. These results are consistent with the hypothesis that Cl enrichment of the IOA-related source magmas played a determinant role in their formation. We propose that the porphyry and skarn Cu ± Au deposits in the Middle-Lower Yangtze River metallogenic belt formed in a transpressional setting in response to paleo-Pacific flat-slab subduction that favored storage and evolution of S-rich hydrous ore-forming magmas at variable crustal levels. A subsequent extensional setting formed due to slab rollback, leading to rapid degassing of Cl-rich IOA-related magmas. For the latter scenario, assimilation of evaporite by mafic to intermediate magmas would lead to an enrichment of Cl in the predegassed magmas and subsequent exsolution of hypersaline magmatic-hydrothermal fluid enriched in Fe as FeCl2. This Fe-rich ore fluid efficiently transported Fe to the apical parts of the magma bodies and overlying extensional normal faults where IOA mineralization was localized. The concomitant loss of S, H2O, and Cu with Cl by volcanic outgassing may have inhibited sulfide mineralization at lower temperatures. 
    more » « less
  4. Let $$f\in C^{2}(\mathbb{T}^{2})$$ have mean value 0 and consider $$\begin{eqnarray}\sup _{\unicode[STIX]{x1D6FE}\,\text{closed geodesic}}\frac{1}{|\unicode[STIX]{x1D6FE}|}\biggl|\int _{\unicode[STIX]{x1D6FE}}f\,d{\mathcal{H}}^{1}\biggr|,\end{eqnarray}$$ where $$\unicode[STIX]{x1D6FE}$$ ranges over all closed geodesics $$\unicode[STIX]{x1D6FE}:\mathbb{S}^{1}\rightarrow \mathbb{T}^{2}$$ and $$|\unicode[STIX]{x1D6FE}|$$ denotes its length. We prove that this supremum is always attained. Moreover, we can bound the length of the geodesic $$\unicode[STIX]{x1D6FE}$$ attaining the supremum in terms of the smoothness of the function: for all $$s\geq 2$$ , $$\begin{eqnarray}|\unicode[STIX]{x1D6FE}|^{s}{\lesssim}_{s}\biggl(\max _{|\unicode[STIX]{x1D6FC}|=s}\Vert \unicode[STIX]{x2202}_{\unicode[STIX]{x1D6FC}}f\Vert _{L^{1}(\mathbb{T}^{2})}\biggr)\Vert \unicode[STIX]{x1D6FB}f\Vert _{L^{2}}\Vert f\Vert _{L^{2}}^{-2}.\end{eqnarray}$$ 
    more » « less
  5. The π-hole above the plane of the X 2 T′Y molecule (T′ = Si, Ge, Sn; X = F, Cl, H; Y = O, S) was allowed to interact with the TH hydride of TH(CH 3 ) 3 (T = Si, Ge, Sn). The resulting TH⋯T′ tetrel bond is quite strong, with interaction energies exceeding 30 kcal mol −1 . F 2 T′O engages in the strongest such bonds, as compared to F 2 T′S, Cl 2 T′O, or Cl 2 T′S. The bond weakens as T′ grows larger as in Si > Ge > Sn, despite the opposite trend in the depth of the π-hole. The reverse pattern of stronger tetrel bond with larger T is observed for the Lewis base TH(CH 3 ) 3 , even though the minimum in the electrostatic potential around the H is nearly independent of T. The TH⋯T′ arrangement is nonlinear which can be understood on the basis of the positions of the extrema in the molecular electrostatic potentials of the monomers. The tetrel bond is weakened when H 2 O forms an O⋯T′ tetrel bond with the second π-hole of F 2 T′O, and strengthened if H 2 O participates in an OH⋯O H-bond. 
    more » « less