Sulfur Isotope Fractionation Research Articles

Sulfur isotope anomalies in coal combustion: Applications to the present and early Earth environments

The observation of mass-independent sulfur isotope fractionations (S-MIF) in Archean-Paleoproterozoic rocks has been instrumental in constraining oxygen levels on early Earth. The S-MIF effect, experimentally demonstrated to result from photochemical reactions, has now been observed in coal combustion, expanding our understanding of this phenomenon. Our study reveals that the negative Δ33S anomalies produced by coal combustion are consistent with similar anomalies observed in present-day sulfate aerosols in Beijing, China, and the black crust sulfates formed on building stones, monument walls, and sculptures in Europe that contribute to carbonate stone deterioration and cultural heritage damage. This finding provides independent evidence for a critical role of atmospheric sulfate from coal combustion in maintaining isotopic balance and offers an effective method for tracing sulfate aerosol sources. These insights are vital for developing more effective regulatory policies to control air pollution and protect public health. Given that coal energy production remains a significant issue in climate science, accurately mapping the global distribution of its by-products is imperative.

Effect of Cd2+ and Cu2+ on Desulfovibrio vulgaris strain ATCC 7757: Insights from sulfur isotope fractionation

Microbial sulfate reduction (MSR) is the predominant microbial metabolic process involved in sulfur cycling in metal-sulfide mining areas. In this study, we investigated the effect of Cd2+ and Cu2+ on MSR driven by Desulfovibrio vulgaris strain ATCC 7757, determined the corresponding sulfur isotope compositions to establish a fundamental relationship between toxicity and isotope fractionation. Cd2+ did not cause significant cellular irrecoverable damage until concentrations exceeded 3.5 ppm, and the bacteria would enhance the transfer of SO42- to APS to ensure the bacterial growth. The repair of DNA damage caused by Cd2+ was a gradually stimulated process, which resulted in a tolerance period. Cu2+ inhibited MSR in a dose-dependent manner at concentrations below 1 ppm, and the bacterial growth significantly inhibited. When exposed to 1 ppm Cu2+ for 48 h, ATCC 7757 might enhance the substrates transport to maintain its culturability through up-regulating genes encoding H+-transporting ATPase, ABC transporters and Protein export. Cd2+ had minimal impact on δ34Ssulfate, with an enrichment factor (34ε) ranging from −11.7 ‰ to −12.3 ‰. Conversely, Cu2+ significantly affected the δ34Ssulfate (34ε = −11.1 ‰ to −15.4 ‰), resulting in a substantial decrease in δ34Ssulfide, with the lowest value of δ34Ssulfide recorded as −18.8 ‰ (-11.5 ‰ for control group). The study demonstrated that sulfur isotope fractionation could indicate different toxicity mechanisms of heavy metal on ATCC 7757, aiding in the identification of environmental contamination caused by heavy metals.

Sulfur isotopic fractionation during hydrolysis of carbonyl sulfide

Carbonyl Sulfide (OCS) is the most abundant sulfur-containing gas in the atmosphere, and it is used as a proxy for terrestrial gross primary productivity (GPP). Oceans are the major source of OCS to the atmosphere, produced by photochemical and “dark” reactions. Hydrolysis to H2S and CO2 is the major removal process of OCS from the ocean's surface. Measuring the sulfur isotope values (δ34S) and the isotopic fractionation (ε) associated with these major OCS sources and sinks could decrease the uncertainties in its fluxes. In the current study, we aim to determine the ε during the hydrolysis process of OCS (εh). We used a purge and trap system coupled to a GC/MC-ICPMS to measure δ34S values during hydrolysis under different temperatures (4–40 °C), salinities (0.2–40 g/L), and pH (4–9), representing various natural environmental conditions. In addition, we use the quantum chemical method to calculate the equilibrium εh and compare it to the empirical results. Our results for the low salinity (S =0.2 g/L; pH 8.0) water show a temperature dependency of the εh from −3.9 ‰ ± 0.2 ‰ (4 °C,) to −2.2 ± 0.6 ‰ (40 °C). The higher fractionation at low temperatures has implication for ice-core data interpretation. However, in natural seawater at 4°C and 22 °C (S = 40 g/L, pH 8.2) there was no such temperature dependency and the εh averaged −2.6 ± 0.3 ‰. Thus, it seems that salinity cancels the temperature effect close to the freezing temperature of water. Varying the pH between 4 and 9 (at 22 °C) did not result in any εh trend. Ab-initio calculations suggest that OCS hydrolysis is not controlled by equilibrium. The εh values we report will aid in quantifying the impact of OCS's hydrolysis on the observable sulfur isotopic signature of OCS in oceanic and in freshwater environments. This in turn will facilitate more accurate mass-balance calculations for the OCS budget from the ocean to the atmosphere.

A new Re-Os age constraint informs the dynamics of the Great Oxidation Event

Abstract The early Paleoproterozoic (ca. 2.5–2.2 Ga) represents a critical juncture in Earth history, marking the inception of an oxygenated atmosphere while bearing witness to potentially multiple widespread and severe glaciations. Deciphering the nature of this glacial epoch and its connection with atmospheric oxygenation has, however, proven difficult, hindered by a reliance on disputed stratigraphic correlations given the paucity of direct radiometric age constraints. Nowhere is this more acute than within the South African Transvaal Supergroup: Here, while the loss of oxygen-sensitive mass-independent sulfur isotope fractionation (S-MIF) has been reported from both the Duitschland and Rooihoogte formations, divided opinion surrounding the time-equivalence of these units has prompted authors to argue for vastly different oxygenation trajectories. Addressing this debate, we present a depositional Re-Os age (2443 ± 33 Ma) from diamictite samples preserved in drillcore of the upper Duitschland Formation. The 100-million-year separation between the Duitschland Formation and its previously presumed equivalent reveals at least two isolated disappearances of S-MIF, requiring that the Great Oxidation Event was dynamic and proceeded via discrete oxygenation episodes whose structure remains incompletely understood. Importantly, our revised framework aligns the lower Duitschland diamictite with the low-latitude glacigenic Makganyene Formation, supporting hypotheses of widespread regional, and potentially global, early Paleoproterozoic glaciation.

Oceanic and Sedimentary Microbial Sulfur Cycling Controlled by Local Organic Matter Flux During the Ediacaran Shuram Excursion in the Three Gorges Area, South China.

The increased difference in the sulfur isotopic compositions of sedimentary sulfate (carbonate-associated sulfate: CAS) and sulfide (chromium-reducible sulfur: CRS) during the Ediacaran Shuram excursion is attributed to increased oceanic sulfate concentration in association with the oxidation of the global ocean and atmosphere. However, recent studies on the isotopic composition of pyrites have revealed that CRS in sediments has diverse origins of pyrites. These pyrites are formed either in the water column/shallow sediments, where the system is open with respect to sulfate, or in deep sediments, where the system is closed with respect to sulfate. The δ34S value of sulfate in the open system is equal to that of seawater; on the contrary, the δ34S value of sulfate in the closed system is higher than that of seawater. Therefore, obtaining the isotopic composition of pyrites formed in an open system, which most likely retain microbial sulfur isotope fractionation, is essential to reconstruct the paleo-oceanic sulfur cycle. In this study, we carried out multiple sulfur isotope analyses of CRS and mechanically separated pyrite grains (>100 μm) using a fluorination method, in addition to secondary ion mass spectrometry (SIMS) analyses of insitu δ34S values of pyrite grains in drill core samples of Member 3 of the Ediacaran Doushantuo Formation in the Three Gorges area, South China. The isotope fractionation of microbial sulfate reduction (MSR) in the limestone layers of the upper part of Member 3 was calculated to be 34ε = 55.7‰ and 33λ = 0.5129 from the δ34S and Δ33S' values of medium-sized pyrite grains ranging from 100 to 300 μm and the average δ34S and Δ33S' values of CAS. Model calculations revealed that the influence of sulfur disproportionation on the δ34S values of these medium-sized pyrite grains was insignificant. In contrast, within the dolostone layers of the middle part of Member 3, isotope fractionation was determined to be 34ε = 47.5‰. The 34ε value in the middle part of Member 3 was calculated from the average δ34S values of the rim of medium-sized pyrite grains and the average δ34S values of CAS. This observation revealed an increase in microbial sulfur isotope fractionation during the Shuram excursion at the drill core site. Furthermore, our investigation revealed correlations between δ34SCRS values and CRS concentrations and between CRS and TOC concentrations, implying that organic matter load to sediments controlled the δ34SCRS values rather than oceanic sulfate concentrations. However, these CRS and TOC concentrations are local parameters that can change only at the kilometer scale with local redox conditions and the intensity of primary production. Therefore, the decreasing δ34SCRS values likely resulted from local redox conditions and not from a global increase in the oceanic sulfate concentration.

Tracing sulfate sources in a tropical agricultural catchment with a stable isotope Bayesian mixing model

Sulfate (SO42−) is an essential anion in drinking water and a vital macronutrient for plant growth. However, elevated sulfate levels can impact ecosystem or human health and could be an important indicator of acid rock drainage or pollution. Therefore, monitoring SO42− sources and transport is important for water quality assessments. This study focused on exploring the sources and transformations of SO42− as well as estimating the proportional contribution of the potential SO42− pollutant sources to groundwater and surface water in a tropical river basin, the Densu River Basin. The study used major ions combined with stable sulfur and oxygen isotope compositions and a Bayesian isotope mixing model, MixSIAR. The major ion characteristics indicate that SO42− concentrations remain stable throughout the rainy and dry seasons but originate from diverse sources. The multi-isotope model (δ34SSO4, δ18OSO4) identified four potential SO42− sources: detergent, precipitation, sewage, and sulfate fertilizer. However, the δ34SSO4 and δ18OSO4 values of the fertilizer source signatures overlapped with those of precipitation and sewage. Nevertheless, the contributions from each source were disentangled using the MixSIAR model, which revealed sewage as the most dominant SO42− pollutant in the Densu Basin, accounting for ~47 % of sulfate in groundwater and ~ 56 % of sulfate in surface water. Sulfate fertilizer (~33 %) was the second most important source after sewage for groundwater, while detergent (~23 %) was the second most important source for surface water. The redox processes of bacterial sulfate reduction and sulfide oxidation were determined to have a minimal impact on the sulfur isotope fractionation within the basin. This study highlights the benefits of combining major ions, sulfur isotopes and the MixSIAR model for identifying sources of sulfate. This approach accounts for uncertainties in source contributions which allows for more robust and reliable apportionment of sulfate sources. The study emphasizes the need for effective waste management and pollution control measures to protect water quality and provides vital guidelines on how to partition sulfate sources on a large catchment scale and evidence for making pollution management decisions on water resources.

Evidence of both molecular cloud and fluid chemistry in Ryugu regolith.

The sulfur chemistry of (162173) Ryugu particles can be a powerful tracer of molecular cloud chemistry and small body processes, but it has not been well explored. We report identification of organosulfurs and a sulfate grain in two Ryugu particles, A0070 and A0093. The sulfate grain shows oxygen isotope ratios (δ17O = -11.0 ± 4.3 per mil, δ18O = -7.8 ± 2.3 per mil) that are akin to silicates in Ryugu but exhibit mass-independent sulfur isotopic fractionation (Δ33S = +5 ± 2 per mil). A methionine-like coating on the sulfate grain is isotopically anomalous (δ15N = +62 ± 2 per mil). Both the sulfate and organosulfurs can simultaneously form and survive during aqueous alteration within Ryugu's parent body, under reduced conditions, low temperature, and a pH >7 in the presence of N-rich organic molecules. This work extends the heliocentric zone where anomalous sulfur, formed by selective photodissociation of H2S gas in the molecular cloud, is found.

Mass-independent fractionation of oxygen and sulfur isotopes

Mass-independent fractionation of oxygen and sulfur isotopes

Oxidized Sulfur Species in Slab Fluids as a Source of Enriched Sulfur Isotope Signatures in Arcs

AbstractRecycling of oxidized sulfur from subducting slabs to the mantle wedge provides simultaneous explanations for the elevated oxygen fugacity (fO2) in subduction zones, their high hydrothermal and magmatic sulfur outputs, and the enriched sulfur isotopic signatures (i.e., δ34S > 0‰) of these outputs. However, a quantitative understanding of the abundance and speciation of sulfur in slab fluids consistent with high pressure experiments is lacking. Here we analyze published experimental data for anhydrite solubility in H2O‐NaCl solutions to calibrate a high‐pressure aqueous speciation model of sulfur within the framework of the deep earth water model. We characterize aqueous complexes, required to account for the high experimental anhydrite solubilities. We then use this framework to predict the speciation and solubility of sulfur in chemically complex fluids in equilibrium with model subducting mafic and ultramafic lithologies, from 2 to 3 GPa and 400 to 800°C at log fO2 from FMQ‐2 to FMQ+4. We show that sulfate complexes of calcium and sodium markedly enhance the stability of sulfate in moderately oxidized fluids in equilibrium with pyrite at fO2 conditions of FMQ+1 to +2, causing large sulfur isotope fractionations up to 10‰ in the fluid relative to the slab. Such fluids could impart oxidized, sulfur‐rich and high δ34S signatures to the mantle wedge that are ultimately transferred to arc magmas, without the need to invoke 34S‐rich subducted lithologies.

High-accuracy measurement of 36SF5 + signal using an ultrahigh-resolution isotope ratio mass spectrometer.

The Δ36S standard deviation measured in a conventional isotope ratio mass spectrometer such as MAT 253 is ca 0.1‰ to 0.3‰. At this precision, it is difficult to resolve the origin of non-mass-dependent sulfur isotope fractionation in tropospheric sulfate aerosol and in Martian meteorites or small deviations from the canonical mass-dependent fractionation laws. Interfering ions with m/z at 131 of 36SF5 + are suggested by the community as the cause of the poor precision, but the exact ion species has not been identified or confirmed. Here we examined the potential interfering ions by using a Thermo Scientific ultrahigh-resolution isotope ratio mass spectrometer to measure SF6 working gas and SF6 gases converted from IAEA-S1/2/3 Ag2S reference materials. We found that there are two resolvable peaks to the right of the 36SF5 + peak when a new filament was installed, which are 186WF4 2+ followed by 12C3F5 +. However, only the 12C3F5 + interference peak was observed after more than three days of filament use. 12C3F5 + is generated inside the instrument during the ionization process. Avoiding the interfering signals, we were able to achieve a Δ36S standard deviation of 0.046‰ (n = 8) for SF6 zero-enrichment and 0.069‰ (n = 8) for overall measurement start from silver sulfide IAEA-S1. Aging the filament with SF6 gas can avoid the interference of 186WF4 2+. Minimizing the presence of carbon-bearing compounds and avoiding the interfering signals of 12C3F5 + from 36SF5 +, we can improve Δ36S measurement accuracy and precision, which helps to open new territories for research using quadruple sulfur isotope composition.

Sedimentary and stratigraphic architecture of the Duitschland and Rooihoogte formations (Palaeoproterozoic, South Africa): implications for tempo of the Great Oxidation Event

Abstract Determining the tempo and causality of key palaeoclimate events recorded in sedimentary strata depends on high-resolution numerical ages and well-constrained stratigraphic correlations at regional and global scale. This requirement is not necessarily met in Precambrian strata due to poorer age resolution, limited preservation, and secondary overprints. A good example includes the Palaeoproterozoic Rooihoogte and Duitschland formations in South Africa, which document the disappearance of mass-independent fractionation of sulphur isotopes (MIF-S) and contain glacial diamictites at their base. They are thus key records of Earth’s surface oxygenation during the Great Oxidation Event (GOE). However, previous studies have either correlated these units, resulting in a unidirectional oxidation trend; or have regarded them as successive strata, causing an interpretation of oscillating oxidation. This study uses extensive outcrop and new core material to investigate correlation between these two units, and to establish depositional models. Results show that key stratigraphic markers can be traced around the entire Kaapvaal craton both in outcrop and the subsurface. In particular, the basal Bevets breccia and the top Duitschland breccia are here re-interpreted as two separate units that are present at the base and top of both formations, supporting correlation of the formations. Consequently, the base of both formations records a major craton-wide event of uplift and karstification, leading to carbonate dissolution and chert brecciation. Erosion of older rocks from across the craton also delivered material for the basal glacial diamictite. The majority of mixed siliciclastic-carbonate sediments were deposited on a storm- and/or delta-influenced shelf. Depositional packages in both formations reflect post-glacial relative sea level rise, followed by progradation of a deltaic, storm or shoreline depositional system. There is a relatively short-lived depositional hiatus to overlying shales of the Timeball Hill Formation. Both Rooihoogte and Duitschland formations thus record only a single glacial event at their base, and a unidirectional trend of oxidation.

Atmospheric oxygenation at the onset of Earth’s Great Oxidation forced enhanced marine anoxia

Abstract Capturing the loss of mass-independent sulphur isotope fractionation (MIF-S), the correlative South African Duitschland and Rooihoogte formations are widely held to bear the isotopic fingerprint of the first atmospheric oxygenation at the onset of the so-called Great Oxidation Event (GOE). Surprisingly, however, while the multiple sulphur isotope systematics of these formations remain central to our understanding of the GOE, until now, comparatively little work has been done to elucidate the repercussions within the marine realm. Here we present chemostratigraphic records from four drill cores covering a large area of the Transvaal Basin, transcending these crucial units and continuing into the overlying Timeball Hill Formation (TBH), that document the immediate, yet counterintuitive, marine response to atmospheric oxygenation. Specifically, irrespective of the interpretative framework employed, our basin-wide redox-sensitive trace element data document an environmental change from oxic/suboxic conditions within the lower and middle parts of the Duitschland and Rooihoogte formations to suboxic/anoxic conditions within their upper reaches. Interestingly, in concert with a ~35‰ negative δ34S excursion that implicates increased sulphate availability and bacterial sulphate reduction, δ98/95Mo3134+0.25 values increase by ~1.0 to 1.5‰. Combining these observations with increased Fe/Mn ratios, elevated total sulphur and carbon contents and a trend towards lower δ13Corg values imply a shift toward less oxygenated conditions across the Transvaal Basin. The combined observations in the mentioned parameters expose a geobiological feedback-driven causality between the earliest oxygenation of the atmosphere and decreased redox potentials of medium to deep marine environments, at least within the Transvaal Basin.

Using sulfur isotopes to constrain the sources of sulfate in PM2.5 during the winter in Jiaozuo City

Sulfate (SO42−) is one of the main driving forces behind the explosive growth of atmospheric fine particulate matter (PM2.5), and it scatters sunlight, contributing to atmospheric cooling. Investigating the sources and transformation of SO42− is crucial for understanding PM2.5 pollution and its implications for climate change. Despite a continuous decline in PM2.5 concentrations in China since 2015, it remains unclear whether there have been changes in the sources and transformation processes of SO42−. Furthermore, previous studies that utilized isotope techniques for source identification have often overlooked the impact of isotope fractionation. This study focuses on atmospheric PM2.5 in Jiaozuo based on samples collected from December 2017 to February 2018. The analysis encompassed both the water-soluble ions and sulfur isotope compositions of SO42− (δ34S–SO42−). By employing a Bayesian isotope mixing model, quantitative analysis of SO42− sources in winter atmospheric PM2.5 in the study area was conducted. The findings revealed that the δ34S–SO42− values of winter PM2.5 ranged from −0.1‰ to 9.8‰, with an arithmetic mean of 4.9 ± 1.5‰ (n = 74). After considering sulfur isotope fractionation effects and the impact of primary SO42−, the calculated results suggested that on average, the contributions of primary SO42−, coal combustion, vehicle exhaust, biomass burning, and oil combustion to the SO42− in winter PM2.5 were 24.9 ± 15.7%, 41.2 ± 2.6%, 27.9 ± 5.3%, 4.8 ± 3.7%, and 1.2 ± 0.7%, respectively. Compared with the clean period, the contribution of the primary SO42− decreased significantly during haze episodes, while the contribution of the secondary SO42− increased. The contribution of coal combustion to the secondary SO42− increased notably by 21.9%, while the contributions from other secondary sulfate remained relatively stable. The lower δ34S values of the atmospheric PM2.5 in Jiaozuo affirm the significant impact of the anthropogenic SO2 input on the sulfur isotope composition of the atmospheric sulfate. Analysis confirmed that the decrease in the contribution of coal combustion to atmospheric SO42− in the study area, indicating the positive effects of regional air pollution control measures. Nevertheless, continued attention should be paid to the influences of coal combustion and vehicle exhaust on regional air pollution.

Direct prediction of isotopic properties from molecular dynamics trajectories: Application to sulfide, sulfate and sulfur radical ions in hydrothermal fluids

In hydrothermal fluids, disulfur (S2•−) and trisulfur (S3•−) radical anions have been observed to coexist with the major hydrogen sulfide and sulfate species. These radical ions have potentially important effects on the solubility, transport and fractionation of metals and sulfur, with consequences for ore deposit formation and, more generally, for geochemical cycles of metals and volatiles. It is therefore essential to know the intrinsic isotopic properties of these important sulfur species in order to use sulfur isotopes for tracing different geological processes. Here, the theoretical equilibrium isotopic properties of the disulfur and trisulfur radical ions are computed and compared to the hydrogen sulfide (H2S) and sulfate ion (SO42−), using, for the first time, a first-principles molecular dynamics (MD) approach. The isotopic properties are calculated directly from molecular dynamics trajectories using the vibrational density of states and the atomic kinetic energy, and then compared to the more established method based on sampling of several snapshots. This comparison allowed us to validate the new modelling method and to assess its advantages and limitations. The predicted equilibrium isotope fractionation in terms of 34S/32S between S2•− and S3•− is small, i.e. <1‰, with a slight enrichment in the heavier isotope for S3•−, over the temperature range 200–500 °C. Both radical ions are slightly depleted in the heavier isotope, by 1 to 2‰, relative to aqueous H2S. Our results help tuning sulfur isotope fractionation models used for tracing the origin and evolution of hydrothermal fluids. Our method opens large perspectives for using the rapidly growing body of MD simulation data in geosciences on structure and stability of aqueous complexes to assess in parallel element isotope fractionations from MD-generated trajectories.

Sulfur Isotope Geochemistry of Ice‐Wedges in Yakutia, East Siberia

ABSTRACTSulfur, with its highly varying stable isotope ratio and involvement in numerous biogeochemical processes, is one of the most widely used elements as an isotopic paleoenvironmental proxy, yet the sulfur isotope ratios of ice‐wedges and their insoluble fraction remain unexplored. This study first presents the sulfur isotopic compositions of soluble sulfate, particulate organic matter (POM), and lithic particles recovered from East Siberian ice‐wedges. Soluble sulfate, primarily representing atmospheric sulfate deposition, shows comparable sulfur isotope ranges in Zyryanka and Batagay, while in Central Yakutia, ice‐wedge sulfate is more enriched in 34S, consistent with the orogenic and cratonic terranes in East Siberia. Given the wedge growth during the inland winter, it is likely that sulfate aerosols were derived mainly from erosion and weathering of regional basement rocks rather than from sea salt spray or biogenic emissions. Within individual ice‐wedges, however, the sulfur isotopic composition of soluble sulfate varies by as much as 7‰, possibly reflecting changes in the relative contributions of sulfur‐isotopically distinct source regions. Beyond the origin of sulfate, greater sulfur isotope fractionations between POM and sulfate during the last glaciation suggest that sulfate may have been anaerobically reduced to sulfide and vice versa in the adjacent root zone.

Energy flux couples sulfur isotope fractionation to proteomic and metabolite profiles in Desulfovibrio vulgaris.

Microbial sulfate reduction is central to the global carbon cycle and the redox evolution of Earth's surface. Tracking the activity of sulfate reducing microorganisms over space and time relies on a nuanced understanding of stable sulfur isotope fractionation in the context of the biochemical machinery of the metabolism. Here, we link the magnitude of stable sulfur isotopic fractionation to proteomic and metabolite profiles under different cellular energetic regimes. When energy availability is limited, cell-specific sulfate respiration rates and net sulfur isotope fractionation inversely covary. Beyond net S isotope fractionation values, we also quantified shifts in protein expression, abundances and isotopic composition of intracellular S metabolites, and lipid structures and lipid/water H isotope fractionation values. These coupled approaches reveal which protein abundances shift directly as a function of energy flux, those that vary minimally, and those that may vary independent of energy flux and likely do not contribute to shifts in S-isotope fractionation. By coupling the bulk S-isotope observations with quantitative proteomics, we provide novel constraints for metabolic isotope models. Together, these results lay the foundation for more predictive metabolic fractionation models, alongside interpretations of environmental sulfur and sulfate reducer lipid-H isotope data.

Procedures from samples to sulfur isotopic data: A review.

Sulfur isotopes have been widely used to solve some key scientific questions, especially in the last two decades with advanced instruments and analytical schemes. Different sulfur speciation and multiple isotopes analyzed in laboratories worldwide and in situ microanalysis have also been reported in many articles. However, methods of sampling to measurements are multifarious, and occasionally some inaccuracies are present in published papers. Vague methods may mislead newcomers to the field, puzzle readers, or lead to incorrect data-based correlations. We have reviewed multiple methods on sulfur isotopic analyses from the perspectives of sampling, laboratory work, and instrumental analysis in order to help reduce operational inhomogeneity and ensure the fidelity of sulfur isotopic data. We do not deem our proposed solutions as the ultimate standard methods but as a lead-in to the overall introduction and summary of the current methods used. It has been shown that external contamination and transformation of different sulfur species should be avoided during the sampling, pretreatment, storage, and chemical treatment processes. Conversion rates and sulfur isotopic fractionations during sulfur extraction, purification, and conversion processes must be verified by researchers using standard or known samples. The unification of absence of isotopic fractionation is needed during all steps, and long-term monitoring of standard samples is recommended. This review compiles more details on different methods in sampling, laboratory operation, and measurement of sulfur isotopes, which is beneficial for researchers' better practice in laboratories. Microanalyses and molecular studies are the frontier techniques that compare the bulk sample with the elemental analysis/continuous flow-gas source stable isotope ratio mass spectrometry method, but the latter is widely used. The development of sulfur isotopic measurements will lead to the innovation in scientific issues with sulfur proxies.

Sulfur isotope and trace element geochemistry of sulfides from the unique Yaojialing Zn-Au-Cu deposit, Lower Yangtze River Metallogenic Belt, China: Implications for ore-forming process and exploration

The Yaojialing Zn-Au‑Cu deposit is a porphyry-skarn-epithermal vein-type compound deposit that consists of epithermal vein-type lead‑zinc‑silver orebodies, skarn Zn-Cu-Au orebodies, and porphyry CuAu orebodies from shallow level to depth. We decipher the source and metallogenic mechanism by studying the trace element and S isotopic compositions of sulfides (pyrite, sphalerite, chalcopyrite, and galena) from three types of orebodies. The collected pyrite can be divided into two types. PyI coexists with chalcopyrite, and PyII coexists with sphalerite or galena. In addition, PyI was further divided into PyIa (collected from porphyry copper bodies) and PyIb (collected from skarn copper bodies). PyII can be further divided into PyIIa (collected from skarn-type lead‑zinc ore body) and PyIIb (collected from vein lead‑zinc ore body in strata). Sphalerite and galena from skarn-type PbZn ore bodies and shallow vein-type PbZn ore bodies are named SpI and GnI, and SpII and GnII, respectively.Pyrite and sphalerite trace element thermometers revealed that the temperature of the ore-forming fluids decreased from 500 to 600 °C in the porphyry ore body to 300– 360 °C in the skarn ore body and then to 240– 280 °C in the shallow vein ore body. The decrease in the Co and Ni contents of pyrite collected from deep to shallow depths may indicate that meteoric water precipitated in the late ore-forming hydrothermal system. Rapid crystallization and variations in the physicochemical states (such as temperature, pH, fO2, and geochemical composition) of the fluids resulted in an obvious oscillating zone of euhedral PyI pyrite particles. It also affects the solubility of trace metal elements and leads to the selective entry of these elements into pyrite. PyII also shows obvious zonation characteristics rather than oscillating zonation, indicating that the growth rate of pyrite is relatively slow. Moreover, the relationships between the Au and As contents in the two types of pyrite are different. PyI is coupled, while PyII shows decoupling. We believe that this phenomenon is due to the high content of As in ore-forming fluid systems, which may inhibit the absorption of Au on the surface of pyrite rather than the decoupling of Au and As caused by rapid crystallization, based on the rapid increase in As content and non-oscillating zonation in PyII. According to the analysis of the sulfur isotopes of various sulfides, the solid evidences suggest that the ore-forming materials of each type of orebody in the Yaojialing deposit were mainly derived from magmatic-hydrothermal processes and that the fractionation of sulfur isotopes was the result of variations in physicochemical conditions caused by magma-hydrothermal evolution rather than the addition of foreign sulfur sources.

Using Io's Sulfur Isotope Cycle to Understand the History of Tidal Heating

AbstractStable isotope fractionation of sulfur offers a window into Io's tidal heating history, which is difficult to constrain because Io's dynamic atmosphere and high resurfacing rates leave it with a young surface. We constructed a numerical model to describe the fluxes in Io's sulfur cycle using literature constraints on rates and isotopic fractionations of relevant processes. Combining our numerical model with measurements of the 34S/32S ratio in Io's atmosphere, we constrain the rates for the processes that move sulfur between reservoirs and model the evolution of sulfur isotopes over time. Gravitational stratification of SO2 in the upper atmosphere, leading to a decrease in 34S/32S with increasing altitude, is the main cause of sulfur isotopic fractionation associated with loss to space. Efficient recycling of the atmospheric escape residue into the interior is required to explain the 34S/32S enrichment magnitude measured in the modern atmosphere. We hypothesize this recycling occurs by SO2 surface frost burial and SO2 reaction with crustal rocks, which founder into the mantle and/or mix with mantle‐derived magmas as they ascend. Therefore, we predict that magmatic SO2 plumes vented from the mantle to the atmosphere will have lower 34S/32S than the ambient atmosphere, yet are still significantly enriched compared to solar‐system average sulfur. Observations of atmospheric variations in 34S/32S with time and/or location could reveal the average mantle melting rate and hence whether the current tidal heating rate is anomalous compared to Io's long‐term average. Our modeling suggests that tides have heated Io for &gt;1.6 Gyr if Io today is representative of past Io.

Atmospheric ozone destruction and the end-Permian crisis: Evidence from multiple sulfur isotopes

The end-Permian mass extinction (EPME) is considered the largest biotic crisis of the Phanerozoic. To explain the worldwide destruction and mutation of land plants, previous work has emphasized the role of enhanced UV irradiation linked to volcanism-induced disruption of the ozone shield. However, direct evidence for a link between volcanic SO2 emission and stratospheric ozone deterioration is missing. Previous bulk analysis (i.e., SF6 method) found no Mass-Independent Fractionation of sulfur isotopes (MIF-S, noted Δ33S). However, possible grain scale sulfur isotope anomaly can be easily overlooked during bulk analysis. To detect possible sulfur isotope anomaly in grain scale, we applied in-situ multiple sulfur isotope analysis of the sulfide using secondary ion mass spectrometry (SIMS). We found a marked positive shift (higher than +0.30‰ and up to +0.94‰) just before the end-Permian extinction in the Meishan section. Modelling shows that these positive Δ33S anomalies cannot be explained by Mass Dependent Fractionation (MDF) processes alone, rather by UV-induced photolysis of volcanic SO2. The formation and preservation of MIF-S anomalies in the atmosphere and subsequent transport to the oceans require intense and prolonged irradiation of the Earth's surface with solar UV, thus fingerprinting severe disruption of the atmospheric ozone balance at the onset of Earth's largest mass extinction.