We present theoretical calculations for all three isotope ratios of sulfur (33S/32S, 34S/32S, 36S/32S) at the B3LYP/6-31+G(d,p) level of theory for aqueous sulfur compounds modeled in 30–40H2O clusters spanning the range of sulfur oxidation state (Sn, n = −2 to +6) for estimating equilibrium fractionation factors in aqueous systems. Computed 34β values based on major isotope (34S/32S) reduced partition function ratios (RPFRs) scale to a first order with sulfur oxidation state and coordination, where 34β generally increase with higher oxidation state and increasing coordination of the sulfur atom. Exponents defining mass dependent relationships based on β values (x/34κ = ln(xβ)/ln(34β), x = 33 or 36) conform to tight ranges over a wide range of temperature for all aqueous compounds (33/34κ ≈ 0.5148–0.5159, 36/34κ ≈ 1.89–1.90 from T ⩾ 0 °C). The exponents converge near a singular value for all compounds at the high temperature limit (33/34κT→∞ = 0.51587 ± 0.00003 and 36/34 κT→∞ = 1.8905 ± 0.0002; 1 s.d. of all computed compounds), and typically follow trends based on oxidation state and coordination similar to those seen in 34β values at lower temperatures. Theoretical equilibrium fractionation factors computed from these β-values are compared to experimental constraints for HSO3−T(aq)/SO2(g, aq), SO2(aq)/SO2(g), H2S(aq)/H2S(g), H2S(aq)/HS−(aq), SO42−(aq)/H2ST(aq), S2O32−(aq) (intramolecular), and S2O32−(aq)/H2ST(aq), and generally agree within a reasonable estimation of uncertainties. We make predictions of fractionation factors where other constraints are unavailable. Isotope partitioning of the isomers of protonated compounds in the sulfite and sulfoxylate systems depend strongly on whether protons are bound to either sulfur or oxygen atoms. The magnitude of the HSO3−T/SO32− major isotope (34S/32S) fractionation factor is predicted to increase with temperature from 0 to 70 °C due to the combined effects of the large magnitude (HS)O3−/SO32− fractionation factor (1000ln34α(HS)bisulfite-sulfite = 19.9‰, 25 °C) relative to the (HO)SO2−/SO32− fractionation factor (1000ln34α(HO)bisulfite–sulfite = −2.2‰, 25 °C), and the increased stability of the (HS)O3− isomer with increasing temperature. We argue that isomerization phenomenon should be considered in models of the sulfur cycle, including models that describe the overall sulfur isotope fractionations associated with microbial metabolism (e.g., microbial sulfate reduction).

We have investigated the quadruple sulfur isotopic composition of inorganic sulfur-bearing phases from 13 carbonaceous chondrites of CM type. Our samples include 4 falls and 9 Antarctic finds. We extracted sulfur from sulfides, sulfates, and elemental sulfur (S0) from all samples. On average, we recover a bulk sulfur (S) content of 2.11 ± 0.39 wt.% S (1σ). The recovered sulfate, S0 and sulfide contents represent 25 ± 12%, 10 ± 7% and 65 ± 15% of the bulk S, respectively (all 1σ). There is no evidence for differences in the bulk S content between falls and finds, and there is no correlation between the S speciation and the extent of aqueous alteration. We report ranges of Δ33S and Δ36S values in CMs that are significantly larger than previously observed. The largest variations are exhibited by S0, with Δ33S values ranging between −0.104 ± 0.012‰ and +0.256 ± 0.018‰ (2σ). The Δ36S/33S ratios of S0 are on average −3.1 ± 1.0 (2σ). Two CMs show distinct Δ36S/33S ratios, of +1.3 ± 0.1 and +0.9 ± 0.1. We suggest that these mass independent S isotopic compositions record H2S photodissociation in the nebula. The varying Δ36S/Δ33S ratios are interpreted to reflect photodissociation that occurred at different UV wavelengths. The preservation of these isotopic features requires that the S-bearing phases were heterogeneously accreted to the CM parent body. Non-zero Δ33S values are also preserved in sulfide and sulfate, and are positively correlated with S0 values. This indicates a genetic relationship between the S-bearing phases: We argue that sulfates were produced by the direct oxidation of S0 (not sulfide) in the parent body. We describe two types of models that, although imperfect, can explain the major features of the CM S isotope compositions, and can be tested in future studies. Sulfide and S0 could both be condensates from the nebula, as the residue and product, respectively, of incomplete H2S photodissociation by UV light (wavelength <150 nm). This idea requires that FeS formation and the S0 condensation co-occur. As an alternative, ice accretion to the CM parent body could allow the delivery of S-MIF in CMs. In that case, sulfides would have been the only S-bearing condensate in CM precursors, and S0 would have been derived from the oxidation of H2S trapped in ices, after its photodissociation at low temperature (<500 K) in the nebula. In our models, the observations of H2S UV photodissociation is required to occur at the disk surface, and allowed in nebular environments with canonical C/O ratios. Vertical motions in the disk would redistribute phases that condensed at high altitude to the midplane, where they accreted in the phases that make up the chondritic matrix.

Variability in the sulfur isotopic composition in sediments can reflect atmospheric, geologic and biological processes. Evidence for ancient fluvio-lacustrine environments at Gale crater on Mars and a lack of efficient crustal recycling mechanisms on the planet suggests a surface environment that was once warm enough to allow the presence of liquid water, at least for discrete periods of time, and implies a greenhouse effect that may have been influenced by sulfur-bearing volcanic gases. Here we report in situ analyses of the sulfur isotopic compositions of SO2 volatilized from ten sediment samples acquired by NASA's Curiosity rover along a 13 km traverse of Gale crater. We find large variations in sulfur isotopic composition that exceed those measured for Martian meteorites and show both depletion and enrichment in 34S. Measured values of Δ34S range from - 47 ± 14‰ to 28 ± 7‰, similar to the range typical of terrestrial environments. Although limited geochronological constraints on the stratigraphy traversed by Curiosity are available, we propose that the observed sulfur isotopic signatures at Gale crater can be explained by equilibrium fractionation between sulfate and sulfide in an impact-driven hydrothermal system and atmospheric processing of sulfur-bearing gases during transient warm periods.