The five stable isotope compositions of Fig Tree barites: Implications on sulfur cycle in ca. 3.2 Ga oceans

The discovery of ³³S anomalies in Archean sedimentary rocks has established that the early Earth before ∼2.2 Ga (billion years ago) had a very different sulfur cycle than today. The origin of the anomalies and the nature of early sulfur cycle are, however, poorly known and debated. In this study, we analyzed the total sulfur and oxygen isotope compositions, the δ¹⁸O, Δ¹⁷O, δ³⁴S, Δ³³S, and Δ³⁶S, for the >3.2 Ga Fig Tree barite deposits from the Barberton Greenstone Belt, South Africa. The goal is to address two questions: (1) was Archean barite sulfate a mixture of ³³S-anomalous sulfate of photolysis origin and ³³S-normal sulfate of other origins? (2) did the underlying photochemical reactions that generated the observed ³³S anomalies for sulfide and sulfate also generate ¹⁷O anomalies for sulfate?We developed a new method in which pure barite sulfate is extracted for oxygen and sulfur isotope measurements from a mixture of barite sands, cherts, and other oxygen-bearing silicates. The isotope data reveal that (1) there is no distinct ¹⁷O anomaly for Fig Tree barite, with an average Δ¹⁷O value the same as that of the bulk Earth (−0.02 ± 0.07‰, N = 49); and (2) the average δ¹⁸O value is +10.6 ± 1.1‰, close to that of the modern seawater sulfate value (+9.3‰). Evidence from petrography and from the δ¹⁸O of barites and co-existing cherts suggest minimum overprinting of later metamorphism on the sulfate’s oxygen isotope composition. Assuming no other processes (e.g., biological) independently induced oxygen isotope exchange between sulfate and water, the lack of reasonable correlation between the δ¹⁸O and Δ³³S or between the δ³⁴S and Δ³³S suggests two mutually exclusive scenarios: (1) An overwhelming majority of the sulfate in the Archean ocean was of photolysis origin, or (2) The early Archean sulfate was a mixture of ³³S-normal sulfates and a small portion (<5%?) of ³³S-anomalous sulfate of photolysis origin from the atmosphere. Scenario 1 requires that sulfate of photolysis origin must have had only small ³³S or ³⁶S anomalies and no ¹⁷O anomaly. Scenario 2 requires that the photolysis sulfate have had highly negative δ³⁴S and Δ³³S values, recommending future theoretical and experimental work to look into photochemical processes that generate sulfate in Quadrant I and sulfide in Quadrant III in a δ³⁴S (X)-Δ³³S (Y) Cartesian plane. A total sulfur and oxygen isotope analysis has provided constraints on the underlying chemical reactions that produced the observed sulfate isotope signature as well as the accompanying atmospheric, oceanic, and biological conditions.     

Fractionation of multiple sulfur isotopes during phototrophic oxidation of sulfide and elemental sulfur by a green sulfur bacterium

We present multiple sulfur isotope measurements of sulfur compounds associated with the oxidation of H₂S and S⁰ by the anoxygenic phototrophic S-oxidizing bacterium Chlorobium tepidum. Discrimination between ³⁴S and ³²S was +1.8 ± 0.5‰ during the oxidation of H₂S to S⁰, and −1.9 ± 0.8‰ during the oxidation of S⁰ to , consistent with previous studies. The accompanying Δ³³S and Δ³⁶S values of sulfide, elemental sulfur, and sulfate formed during these experiments were very small, less than 0.1‰ for Δ³³S and 0.9‰ for Δ³⁶S, supporting mass conservation principles. Examination of these isotope effects within a framework of the metabolic pathways for S oxidation suggests that the observed effects are due to the flow of sulfur through the metabolisms, rather than abiotic equilibrium isotope exchange alone, as previously suggested. The metabolic network comparison also indicates that these metabolisms work to express some isotope effects (between sulfide, polysulfides, and elemental sulfur in the periplasm) and suppress others (kinetic isotope effects related to pathways for oxidation of sulfide to sulfate via the same enzymes involved in sulfate reduction acting in reverse). Additionally, utilizing fractionation factors for phototrophic S oxidation calculated from our experiments and for other oxidation processes calculated from the literature (chemotrophic and inorganic S oxidation), we constructed a set of ecosystem-scale sulfur isotope box models to examine the isotopic consequences of including sulfide oxidation pathways in a model system. These models demonstrate how the small δ³⁴S effects associated with S oxidation combined with large δ³⁴S effects associated with sulfate reduction (by SRP) and sulfur disproportionation (by SDP) can produce large (and measurable) effects in the Δ³³S of sulfur reservoirs. Specifically, redistribution of material along the pathways for sulfide oxidation diminishes the net isotope effect of SRP and SDP, and can mask the isotopic signal for sulfur disproportionation if significant recycling of S intermediates occurs. We show that the different sulfide oxidation processes produce different isotopic fields for identical proportions of oxidation, and discuss the ecological implications of these results to interpreting minor S isotope patterns in modern systems and in the geologic record.     

The Paleoarchean Northern Mundo Novo Greenstone Belt,São Francisco Craton: Geochemistry,U–Pb–Hf–O in zircon and pyrite δ34S-Δ33S-Δ36S signatures

Greenstone belts contain several clues about the evolutionary history of primitive Earth. Here, we describe the volcano-sedimentary rock association exposed along the eastern margin of the Gavião Block, named the Northern Mundo Novo Greenstone Belt (N-MNGB), and present data collected with different techniques, including U–Pb–Hf–O isotopes of zircon and multiple sulfur isotopes (³²S, ³³S, ³⁴S, and ³⁶S) of pyrite from this supracrustal sequence. A pillowed metabasalt situated in the upper section of the N-MNGB is 3337 ± 25 Ma old and has zircon with ε_Hf(t) =  ?2.47 to ?1.40, Hf model ages between 3.75 Ga and 3.82 Ga, and δ¹⁸O = +3.6‰ to +7.3‰. These isotopic data, together with compiled whole-rock trace element data, suggest that the mafic metavolcanic rocks formed in a subduction-related setting, likely a back-arc basin juxtaposed to a continental arc. In this context, the magma interacted with older Eoarchean crustal components from the Gavião Block. Detrital zircons from the overlying quartzites of the Jacobina Group are sourced from Paleoarchean rocks, in accordance with previous studies, yielding a maximum depositional age of 3353 ± 22 Ma. These detrital zircons have ε_Hf(t) =  ?5.40 to ?0.84, Hf model ages between 3.66 Ga and 4.30 Ga, and δ¹⁸O = +4.8‰ to +6.4‰. The pyrite multiple sulfur isotope investigation of the 3.3 Ga supracrustal rocks from the N-MNGB enabled a further understanding of Paleoarchean sulfur cycling. The samples have diverse isotopic compositions that indicate sulfur sourced from distinct reservoirs. Significantly, they preserve the signal of the anoxic Archean atmosphere, expressed by MIF-S signatures (Δ³³S between ?1.3‰ to +1.4‰) and a Δ³⁶S/Δ³³S slope of ?0.81 that is indistinguishable from the so-called Archean array. A BIF sample has a magmatic origin of sulfur, as indicated by the limited δ³⁴S range (0 to +2‰), Δ³³S ~ 0‰, and Δ³⁶S ~ 0‰. A carbonaceous schist shows positive δ³⁴S (2.1‰–3.5‰) and elevated Δ³³S (1.2‰–1.4‰) values, with corresponding negative Δ³⁶S between ?1.2‰ to ?0.2‰, which resemble the isotopic composition of Archean black shales and suggest a source from the photolytic reduction of elemental sulfur. The pillowed metabasalt displays heterogeneous δ³⁴S, Δ³³S, and Δ³⁶S signatures that reflect assimilation of both magmatic sulfur and photolytic sulfate during hydrothermal seafloor alteration. Lastly, pyrite in a massive sulfide lens is isotopically similar to barite of several Paleoarchean deposits worldwide, which might indicate mass dependent sulfur processing from a global and well-mixed sulfate reservoir at this time.     

Characterization of calcium isotopes in natural and synthetic barite

The mineral barite (BaSO₄) accommodates calcium in its crystal lattice, providing an archive of Ca-isotopes in the highly stable sulfate mineral. Holocene marine (pelagic) barite samples from the major ocean basins are isotopically indistinguishable from each other (δ^44/40Ca = −2.01 ± 0.15‰) but are different from hydrothermal and cold seep barite samples (δ^44/40Ca = −4.13 to −2.72‰). Laboratory precipitated (synthetic) barite samples are more depleted in the heavy Ca-isotopes than pelagic marine barite and span a range of Ca-isotope compositions, Δ^44/40Ca = −3.42 to −2.40‰. Temperature, saturation state, , and aCa²⁺/aBa²⁺ each influence the fractionation of Ca-isotopes in synthetic barite; however, the fractionation in marine barite samples is not strongly related to any measured environmental parameter. First-principles lattice dynamical modeling predicts that at equilibrium Ca-substituted barite will have much lower ⁴⁴Ca/⁴⁰Ca than calcite, by −9‰ at 0 °C and −8‰ at 25 °C. Based on this model, none of the measured barite samples appear to be in isotopic equilibrium with their parent solutions, although as predicted they do record lower δ^44/40Ca values than seawater and calcite. Kinetic fractionation processes therefore most likely control the extent of isotopic fractionation exhibited in barite. Potential fractionation mechanisms include factors influencing Ca²⁺ substitution for Ba²⁺ in barite (e.g. ionic strength and trace element concentration of the solution, competing complexation reactions, precipitation or growth rate, temperature, pressure, and saturation state) as well as nucleation and crystal growth rates. These factors should be considered when investigating controls on isotopic fractionation of Ca²⁺ and other elements in inorganic and biogenic minerals.     

Application of sulphur isotope ratios to examine weaning patterns and freshwater fish consumption in Roman Oxfordshire, UK

This study investigates the application of sulphur isotope ratios (δ³⁴S) in combination with carbon (δ¹³C) and nitrogen (δ¹⁵N) ratios to understand the influence of environmental sulphur on the isotopic composition of archaeological human and faunal remains from Roman era sites in Oxfordshire, UK. Humans (n = 83), terrestrial animals (n = 11), and freshwater fish (n = 5) were analysed for their isotope values from four locations in the Thames River Valley, and a broad range of δ³⁴S values were found. The δ³⁴S values from the terrestrial animals were highly variable (−13.6‰ to +0.5‰), but the δ³⁴S values of the fish were clustered and ³⁴S-depleted (−20.9‰ to −17.3‰). The results of the faunal remains suggest that riverine sulphur influenced the terrestrial sulphur isotopic signatures. Terrestrial animals were possibly raised on the floodplains of the River Thames, where highly ³⁴S-depleted sulphur influenced the soil. The humans show the largest range of δ³⁴S values (−18.8‰ to +9.6‰) from any archaeological context to date. No differences in δ³⁴S values were found between the males (−7.8 ± 6.0‰) and females (−5.3 ± 6.8‰), but the females had a linear correlation (R² = 0.71; p < 0.0001) between their δ¹⁵N and δ³⁴S compositions. These δ³⁴S results suggest a greater dietary variability for the inhabitants of Roman Oxfordshire than previously thought, with some individuals eating solely terrestrial protein resources and others showing a diet almost exclusively based on freshwater protein such as fish. Such large dietary variability was not visible by analysing only the carbon and nitrogen isotope ratios, and this research represents the largest and most detailed application of δ³⁴S analysis to examine dietary practices (including breastfeeding and weaning patterns) during the Romano-British Period.     

S-33 constraints on the seawater sulfate contribution in modern seafloor hydrothermal vent sulfides

Sulfide sulfur in mid-oceanic ridge hydrothermal vents is derived from leaching of basaltic-sulfide and seawater-derived sulfate that is reduced during high temperature water rock interaction. Conventional sulfur isotope studies, however, are inconclusive about the mass-balance between the two sources because ³⁴S/³²S ratios of vent fluid H₂S and chimney sulfide minerals may reflect not only the mixing ratio but also isotope exchange between sulfate and sulfide. Here, we show that high-precision analysis of S-33 can provide a unique constraint because isotope mixing and isotope exchange result in different Δ³³S (≡δ³³S-0.515 δ³⁴S) values of up to 0.04‰ even if δ³⁴S values are identical. Detection of such small Δ³³S differences is technically feasible by using the SF₆ dual-inlet mass-spectrometry protocol that has been improved to achieve a precision as good as 0.006‰ (2σ).Sulfide minerals (marcasite, pyrite, chalcopyrite, and sphalerite) and vent H₂S collected from four active seafloor hydrothermal vent sites, East Pacific Rise (EPR) 9-10°N, 13°N, and 21°S and Mid-Atlantic Ridge (MAR) 37°N yield Δ³³S values ranging from −0.002 to 0.033 and δ³⁴S from −0.5‰ to 5.3‰. The combined δ³⁴S and Δ³³S systematics reveal that 73 to 89% of vent sulfides are derived from leaching from basaltic sulfide and only 11 to 27% from seawater-derived sulfate. Pyrite from EPR 13°N and marcasite from MAR 37°N are in isotope disequilibrium not only in δ³⁴S but also in Δ³³S with respect to associated sphalerite and chalcopyrite, suggesting non-equilibrium sulfur isotope exchange between seawater sulfate and sulfide during pyrite precipitation. Seafloor hydrothermal vent sulfides are characterized by low Δ³³S values compared with biogenic sulfides, suggesting little or no contribution of sulfide from microbial sulfate reduction into hydrothermal sulfides at sediment-free mid-oceanic ridge systems. We conclude that ³³S is an effective new tracer for interplay among seawater, oceanic crust and microbes in subseafloor hydrothermal sulfur cycles.     

Rare sulfur and triple oxygen isotope geochemistry of volcanogenic sulfate aerosols

We present analyses of stable isotopic ratios ¹⁷O/¹⁶O, ¹⁸O/¹⁶O, ³⁴S/³²S, and ³³S/³²S, ³⁶S/³²S in sulfate leached from volcanic ash of a series of well known, large and small volcanic eruptions. We consider eruptions of Mt. St. Helens (Washington, 1980, ∼1 km³), Mt. Spurr (Alaska, 1953, <1 km³), Gjalp (Iceland, 1996, 1998, <1 km³), Pinatubo (Phillipines, 1991, 10 km³), Bishop tuff (Long Valley, California, 0.76 Ma, 750 km³), Lower Bandelier tuff (Toledo Caldera, New Mexico, 1.61 Ma, 600 km³), and Lava Creek and Huckleberry Ridge tuffs (Yellowstone, Wyoming, 0.64 Ma, 1000 km³ and 2.04 Ma 2500 km³, respectively). This list covers much of the diversity of sizes and the character of silicic volcanic eruptions. Particular emphasis is paid to the Lava Creek tuff for which we present wide geographic sample coverage.This global dataset spans a significant range in δ³⁴S, δ¹⁸O, and Δ¹⁷O of sulfate (29‰, 30‰, and 3.3‰, respectively) with oxygen isotopes recording mass-independent (Δ¹⁷O > 0.2‰) and sulfur isotopes exhibiting mass-dependent behavior. Products of large eruptions account for most of‘ these isotopic ranges. Sulfate with Δ¹⁷O > 0.2‰ is present as 1-10 μm gypsum crystals on distal ash particles and records the isotopic signature of stratospheric photochemical reactions. Sediments that embed ash layers do not contain sulfate or contain little sulfate with Δ¹⁷O near 0‰, suggesting that the observed sulfate in ash is of volcanic origin.Mass-dependent fractionation of sulfur isotopic ratios suggests that sulfate-forming reactions did not involve photolysis of SO₂, like that inferred for pre-2.3 Ga sulfates from Archean sediments or Antarctic ice-core sulfate associated with few dated eruptions. Even though the sulfate sulfur isotopic compositions reflect mass-dependent processes, the products of caldera-forming eruptions display a large δ³⁴S range and exhibit fractionation relationships that do not follow the expected equilibrium slopes of 0.515 and 1.90 for ³³S/³²S vs. ³⁴S/³²S and ³⁶S/³²S vs. ³⁴S/³²S, respectively. The data presented here are consistent with modification of a chemical mass-dependent fractionation of sulfur isotopes in the volcanic plume by either a kinetic gas phase reaction of volcanic SO₂ with OH and/or a Rayleigh processes involving a residual Rayleigh reactant—volcanic SO₂ gas, rather than a Rayleigh product. These results may also imply at least two removal pathways for SO₂ in volcanic plumes.Above-zero Δ¹⁷O values and their positive correlation with δ¹⁸O in sulfate can be explained by oxidation by high-δ¹⁸O and high-Δ¹⁷O compounds such as ozone and radicals such as OH that result from ozone break down. Large caldera-forming eruptions have the highest Δ¹⁷O values, and the largest range of δ¹⁸O, which can be explained by stratospheric reaction with ozone-derived OH radicals. These results suggest that massive eruptions are capable of causing a temporary depletion of the ozone layer. Such depletion may be many times that of the measured 3-8% depletion following 1991 Pinatubo eruption, if the amount of sulfur dioxide released scales with the amount of ozone depletion.     

Constraining atmospheric oxygen and seawater sulfate concentrations during Paleoproterozoic glaciation: In situ sulfur three-isotope microanalysis of pyrite from the Turee Creek Group, Western Australia

Previous efforts to constrain the timing of Paleoproterozoic atmospheric oxygenation have documented the disappearance of large, mass-independent sulfur isotope fractionation and an increase in mass-dependent sulfur isotope fractionation associated with multiple glaciations. At least one of these glacial events is preserved in diamictites of the ∼2.4 Ga Meteorite Bore Member of the Kungarra Formation, Turee Creek Group, Western Australia. Outcrop exposures of this unit show the transition from the Boolgeeda Iron Formation of the upper Hamersley Group into clastic, glaciomarine sedimentary rocks of the Turee Creek Group. Here we report in situ multiple sulfur isotope and elemental abundance measurements of sedimentary pyrite at high spatial resolution, as well as the occurrence of detrital pyrite in the Meteorite Bore Member. The 15.3‰ range of Δ³³S in one sample containing detrital pyrite (−3.6‰ to 11.7‰) is larger than previously reported worldwide, and there is evidence for mass-independent sulfur isotope fractionation in authigenic pyrite throughout the section (Δ³³S from −0.8‰ to 1.0‰). The 90‰ range in δ³⁴S observed (−45.5‰ to 46.4‰) strongly suggests microbial sulfate reduction under non-sulfate limiting conditions, indicating significant oxidative weathering of sulfides on the continents. Multiple generations of pyrite are preserved, typically represented by primary cores with low δ³⁴S (<−20‰) overgrown by euhedral rims with higher δ³⁴S (4-7‰) and enrichments in As, Ni, and Co. The preservation of extremely sharp sulfur isotope gradients (30‰/<4 μm) implies limited sulfur diffusion and provides time and temperature constraints on the metamorphic history of the Meteorite Bore Member. Together, these results suggest that the Meteorite Bore Member was deposited during the final stages of the “Great Oxidation Event,” when pO₂ first became sufficiently high to permit pervasive oxidative weathering of continental sulfides, yet remained low enough to permit the production and preservation of mass-independent sulfur isotope fractionation.     

The elemental and isotopic composition of sulfur and nitrogen in Chinese coals

Coal combustion is an important atmospheric pollution source in most Chinese cities, so systematic studies on sulfur and nitrogen in Chinese coals are needed. The sulfur contents in Chinese coals average 0.9 ± 1.0%, indicating that most Chinese coals are low in sulfur. A nearly constant mean δ³⁴S value is observed in low sulfur (TS < 1) Chinese coals of different ages (D, P₁, T₃ and J₃). High sulfur Chinese coals (OS > 0.8%), often found at late Carboniferous (C₃) and late Permian (P₂) in southern China, had two main sulfur sources (original plant sulfur and secondary sulfur). The wide variety of δ³⁴S values of Chinese coals (−15‰ to +50‰) is a result of a complex sulfur origin. The δ¹⁵N values of Chinese coals ranged from −6‰ to +4‰, showing a lack of correlation with coal ages, whereas nitrogen contents are higher in Paleozoic coals than in Mesozoic coals. This may be related to their original precursor plant species: high nitrogen pteridophytes for the Paleozoic coals and low nitrogen gymnosperms for the Mesozoic coals. Different to δ³⁴S values, Chinese coals showed higher δ¹⁵N values in marine environments than in freshwater environments.     

Abundance of mass 47 CO2 in urban air, car exhaust, and human breath

Atmospheric carbon dioxide is widely studied using records of CO₂ mixing ratio, δ¹³C and δ¹⁸O. However, the number and variability of sources and sinks prevents these alone from uniquely defining the budget. Carbon dioxide having a mass of 47 u (principally ¹³C¹⁸O¹⁶O) provides an additional constraint. In particular, the mass 47 anomaly (Δ₄₇) can distinguish between CO₂ produced by high temperature combustion processes vs. low temperature respiratory processes. Δ₄₇ is defined as the abundance of mass 47 isotopologues in excess of that expected for a random distribution of isotopes, where random distribution means that the abundance of an isotopologue is the product of abundances of the isotopes it is composed of and is calculated based on the measured ¹³C and ¹⁸O values. In this study, we estimate the δ¹³C (vs. VPDB), δ¹⁸O (vs. VSMOW), δ47, and Δ₄₇ values of CO₂ from car exhaust and from human breath, by constructing ‘Keeling plots’ using samples that are mixtures of ambient air and CO₂ from these sources. δ47 is defined as , where is the R⁴⁷ value for a hypothetical CO₂ whose δ¹³C_VPDB = 0, δ¹⁸O_VSMOW = 0, and Δ₄₇ = 0. Ambient air in Pasadena, CA, where this study was conducted, varied in [CO₂] from 383 to 404 μmol mol⁻¹, in δ¹³C and δ¹⁸O from −9.2 to −10.2‰ and from 40.6 to 41.9‰, respectively, in δ47 from 32.5 to 33.9‰, and in Δ₄₇ from 0.73 to 0.96‰. Air sampled at varying distances from a car exhaust pipe was enriched in a combustion source having a composition, as determined by a ‘Keeling plot’ intercept, of −24.4 ± 0.2‰ for δ¹³C (similar to the δ¹³C of local gasoline), δ¹⁸O of 29.9 ± 0.4‰, δ47 of 6.6 ± 0.6‰, and Δ₄₇ of 0.41 ± 0.03‰. Both δ¹⁸O and Δ₄₇ values of the car exhaust end-member are consistent with that expected for thermodynamic equilibrium at∼200 °C between CO₂ and water generated by combustion of gasoline-air mixtures. Samples of CO₂ from human breath were found to have δ¹³C and δ¹⁸O values broadly similar to those of car exhaust-air mixtures, −22.3 ± 0.2 and 34.3 ± 0.3‰, respectively, and δ47 of 13.4 ± 0.4‰. Δ₄₇ in human breath was 0.76  ± 0.03‰, similar to that of ambient Pasadena air and higher than that of the car exhaust signature.     

Sulfur isotope signatures in gypsiferous sediments of the Estancia and Tularosa Basins as indicators of sulfate sources, hydrological processes, and microbial activity

In order to reconstruct paleo-environmental conditions for the saline playa lakes of the Rio Grande Rift, we investigated sediment sulfate sources using sulfur isotope compositions of dissolved ions in modern surface water, groundwater, and precipitated in the form of gypsum sediments deposited during the Pleistocene and Holocene in the Tularosa and Estancia Basins. The major sulfate sources are Lower and Middle Permian marine evaporites (δ³⁴S of 10.9-14.4‰), but the diverse physiography of the Tularosa Basin led to a complex drainage system which contributed sulfates from various sources depending on the climate at the time of sedimentation. As inferred from sulfur isotope mass balance constraints, weathering of sulfides of magmatic/hydrothermal and sedimentary origin associated with climate oscillations during Last Glacial Maximum contributed about 35-50% of the sulfates and led to deposition of gypsum with δ³⁴S values of −1.2‰ to 2.2‰ which are substantially lower than Permian evaporates. In the Estancia Basin, microbial sulfate reduction appears to overprint sulfur isotopic signatures that might elucidate past groundwater flows. A Rayleigh distillation model indicates that about 3-18% of sulfates from an inorganic groundwater pool (δ³⁴S of 12.6-13.8‰) have been metabolized by bacteria and preserved as partially to fully reduced sulfur-bearing minerals species (elemental sulfur, monosulfides, disulfides) with distinctly negative δ³⁴S values (−42.3‰ to −20.3‰) compared to co-existing gypsum (−3.8‰ to 22.4‰). For the Tularosa Basin microbial sulfate reduction had negligible effect on δ³⁴S value of the gypsiferous sediments most likely because of higher annual temperatures (15-33 °C) and lower organic carbon content (median 0.09%) in those sediments leading to more efficient oxidation of H₂S and/or smaller rates of sulfate reduction compared to the saline playas of the Estancia Basin (5-28 °C; median 0.46% of organic carbon).The White Sands region of the Tularosa Basin is frequently posited as a hydrothermal analogue for Mars. High temperatures of groundwater (33.3 °C) and high δ¹⁸O(H₂O) values (1.1‰) in White Sands, however, are controlled predominantly by seasonal evaporation rather than the modern influx of hydrothermal fluids. Nevertheless, it is possible that some of the geochemical processes in White Sands, such as sulfide weathering during climate oscillations and upwelling of highly mineralized waters, might be considered as valid terrestrial analogues for the sulfate cycle in places such as Meridiani Planum on Mars.     

The oxygen-isotope composition of chondrules and isolated forsterite and olivine grains from the Tagish Lake carbonaceous chondrite

The oxygen-isotope compositions (obtained by laser fluorination) of hand-picked separates of isolated forsterite, isolated olivine and chondrules from the Tagish Lake carbonaceous chondrite describe a line (δ¹⁷O = 0.95 * δ¹⁸O − 3.24; R² = 0.99) similar to the trend known for chondrules from other carbonaceous chondrites. The isolated forsterite grains (Fo_99.6-99.8; δ¹⁸O = −7.2‰ to −5.5‰; δ¹⁷O = −9.6‰ to −8.2‰) are more ¹⁶O-rich than the isolated olivine grains (Fo_39.6-86.8; δ¹⁸O = 3.1‰ to 5.1‰; δ¹⁷O = −0.3‰ to 2.2‰), and have chemical and isotopic characteristics typical of refractory forsterite. Chondrules contain olivine (Fo_97.2-99.8) with oxygen-isotope compositions (δ¹⁸O = −5.2‰ to 5.9‰; δ¹⁷O = −8.1‰ to 1.2‰) that overlap those of isolated forsterite and isolated olivine. An inverse relationship exists between the Δ¹⁷O values and Fo contents of Tagish Lake isolated forsterite and chondrules; the chondrules likely underwent greater exchange with ¹⁶O-poor nebular gases than the forsterite. The oxygen-isotope compositions of the isolated olivine grains describe a trend with a steeper slope (1.1 ± 0.1, R² = 0.94) than the carbonaceous chondrite anhydrous mineral line (CCAM_slope = 0.95). The isolated olivine may have crystallized from an evolving melt that exchanged with ¹⁶O-poor gases of somewhat different composition than those which affected the chondrules and isolated forsterite. The primordial components of the Tagish Lake meteorite formed under conditions similar to other carbonaceous chondrite meteorite groups, especially CMs. Its alteration history has its closest affinities to CI carbonaceous chondrites.     

Sulfur isotope fractionation during the May 2003 eruption of Anatahan volcano, Mariana Islands: Implications for sulfur sources and plume processes

Sulfur isotope compositions of pumice and adsorbed volatiles on ash from the first historical eruption of Anatahan volcano (Mariana arc) are presented in order to constrain the sources of sulfur erupted during the period 10-21 May, 2003. The isotopic composition of S extracted from erupted pumice has a narrow range, from δ³⁴S_V-CDT +2.6‰ to +3.2‰, while the composition of sulfur adsorbed onto ash has a larger range (+2.8‰ to +5.3‰). Fractionation modeling for closed and open system scenarios suggests that degassing of SO₂ raised the δ³⁴S_V-CDT value of S dissolved in the melt from an initial composition of between +1.6‰ and +2.6‰ for closed-system degassing, or between −0.5‰ and +1.5‰ for open-system degassing, however closed-system degassing is the preferred model. The calculated values for the initial composition of the magma represent a MORB-like (δ³⁴S_V-CDT ∼ 0‰) mantle source with limited contamination by subducted seawater sulfate (δ³⁴S_V-CDT +21‰). Modeling also suggests that the δ³⁴S_V-CDT value of SO₂ gas in closed-system equilibrium with the degassed magma was between +0.9‰ and +2.5‰. The δ³⁴S_V-CDT value of sulfate adsorbed onto ash in the eruption plume (+2.8‰ to +5.1‰) is consistent with sulfate formation by oxidation of magmatic SO₂ in the eruption column. The sulfur isotope composition of sulfate adsorbed to ash changes from lower δ³⁴S values for ash erupted early in the eruption to higher δ³⁴S values for ash erupted later in the eruption. We interpret the temporal/stratigraphic change in sulfate isotopic composition to primarily reflect a change in the isotopic composition of magmatic SO₂ released from the progressively degassing magma and is attributed to the expulsion of an accumulated gas phase at the beginning of the eruption. More efficient oxidation of magmatic SO₂ gas to sulfate in the early water-rich eruption plume probably contributed to the change in S isotope compositions observed in the ash leachates.     

Sulfur cycling in a stratified euxinic lake with moderately high sulfate: Constraints from quadruple S isotopes

We present a 3-year study of concentrations and sulfur isotope values (δ³⁴S, Δ³³S, and Δ³⁶S) of sulfur compounds in the water column of Fayetteville Green Lake (NY, USA), a stratified (meromictic) euxinic lake with moderately high sulfate concentrations (12-16 mM). We utilize our results along with numerical models (including transport within the lake) to identify and quantify the major biological and abiotic processes contributing to sulfur cycling in the system. The isotope values of sulfide and zero-valent sulfur across the redox-interface (chemocline) change seasonally in response to changes in sulfide oxidation processes. In the fall, sulfide oxidation occurs primarily via abiotic reaction with oxygen, as reflected by an increase in sulfide δ³⁴S at the redox interface. Interestingly, S isotope values for zero-valent sulfur sampled at this time still reflect production and recycling by phototrophic S-oxidation. In the spring, sulfide S isotope values suggest an increased input from phototrophic oxidation, consistent with a more pronounced phototroph population at the chemocline. This trend is associated with smaller fractionations between sulfide and zero-valent sulfur, suggesting a metabolic rate control on fractionation similar to that for sulfate reduction. Comparison of our data with previous studies indicates that the S isotope values of sulfate and sulfide in the deep waters are remarkably stable over long periods of time, with consistently large fractionations of up to 58‰ in δ³⁴S. Models of the δ³⁴S and Δ³³S trends in the deep waters (considering mass transport via diffusion and advection along with biological processes) require that these fractionations are a consequence of sulfur compound disproportionation at and below the redox interface in addition to large fractionations during sulfate reduction. The large fractionations during sulfate reduction appear to be a consequence of the high sulfate concentrations and the distribution of organic matter in the water column. The occurrence of disproportionation in the lake is supported by profiles of intermediate sulfur compounds and by lake microbiology, but is not evident from the δ³⁴S trends alone. These results illustrate the utility of including minor S isotopes in sulfur isotope studies to unravel complex sulfur cycling in natural systems.     

Oxygen isotope systematics of chondrules in the Allende CV3 chondrite: High precision ion microprobe studies

The oxygen three-isotope systematics of 36 chondrules from the Allende CV3 chondrite are reported using high precision secondary ion mass spectrometer (CAMECA IMS-1280). Twenty-six chondrules have shown internally homogenous Δ¹⁷O values among olivine, pyroxene, and spinel within a single chondrule. The average Δ¹⁷O values of 19 FeO-poor chondrules (13 porphyritic chondrules, 2 barred olivine chondrules, and 4 chondrule fragments) show a peak at −5.3 ± 0.6‰ (2SD). Another 5 porphyritic chondrules including both FeO-poor and FeO-rich ones show average Δ¹⁷O values between −3‰ and −2‰, and 2 other FeO-poor barred olivine chondrules show average Δ¹⁷O values of −3.6‰ and 0‰. These results are similar to those for Acfer 094 chondrules, showing bimodal Δ¹⁷O values at −5‰ and −2‰. Nine porphyritic chondrules contain olivine grains with heterogeneous Δ¹⁷O values as low as −18‰, indicating that they are relict olivine grains and some of them were derived from precursors related to refractory inclusions. However, most relict olivine grains show oxygen isotope ratios that overlap with those in homogeneous chondrules. The Δ¹⁷O values of four barred olivine chondrules range from −5‰ to 0‰, indicating that not all BO chondrules plot near the terrestrial fractionation line as suggested by previous bulk chondrule analyses. Based on these data, we suggest the presence of multiple oxygen isotope reservoirs in local dust-rich protoplanetary disk, from which the CV3 parent asteroid formed.A compilation of 225 olivine and low-Ca pyroxene isotopic data from 36 chondrules analyzed in the present study lie between carbonaceous chondrite anhydrous mineral (CCAM) and Young and Russell lines. These data define a correlation line of δ¹⁷O = (0.982 ± 0.019) × δ¹⁸O − (2.91 ± 0.10), which is similar to those defined by chondrules in CV3 chondrites and Acfer 094 in previous studies. Plagioclase analyses in two chondrules plot slightly below the CCAM line with Δ¹⁷O values of −2.6‰, which might be the result of oxygen isotope exchange between chondrule mesostasis and aqueous fluid in the CV parent body.     

Implications of conservation of mass effects on mass-dependent isotope fractionations: Influence of network structure on sulfur isotope phase space of dissimilatory sulfate reduction

This paper presents worked solutions for the fractionations of all four stable sulfur isotopes (³²S, ³³S, ³⁴S, and ³⁶S) in several models of the sulfate reduction metabolism. We describe methods for obtaining solutions and how the predictions made by these solutions define different compositional fields (phase space) that can be used to gain new insights into sulfur metabolisms, specifically with respect to understanding the structure of and fractionations associated with the network of reactions that describe the transformations of sulfur within the cell. We show how this treatment can be used to evaluate data from experiments with dissimilatory sulfate reducers and to suggest that the expression of fractionations by the metabolic process is largely limited by the fraction of sulfate that is lost from the cell, and that the variation in observed fractionations reflects differences in the proportion of sulfur intermediates that are reoxidized to sulfate. This analysis provides a line of support for this assertion that depends only on the sulfur isotopic fractionations between sulfate and sulfide. This analysis also indicates that internal fractionations are consistent with a relationship given by ³³α = (³⁴α)^θ where α is the fractionation factor (e.g., ³³α_a−b = (³³ S/³²S)_a/(³³S/³²S)_b and ³⁴α_a−b = (³⁴ S/³²S)_a/(³⁴S/³²S)_b), and where θ is restricted to a value between 0.515 and 0.514. This finding is consistent with a control on isotopic fractionation effects within the cell that is rooted in the different partition coefficients (energetics) for the different isotopologs.     

Non-biological fractionation of stable Ca isotopes in soils of the Atacama Desert, Chile

We measured Ca stable isotope ratios (δ^44/40Ca) in an ancient (2 My), hyperarid soil where the primary source of mobile Ca is atmospheric deposition. Most of the Ca in the upper meter of this soil (3.5 kmol m⁻²) is present as sulfates (2.5 kmol m⁻²), and to a lesser extent carbonates (0.4 kmol m⁻²). In aqueous extracts of variably hydrated calcium sulfate minerals, δ^44/40Ca_E values (vs. bulk Earth) increase with depth (1.4 m) from a minimum of −1.91‰ to a maximum of +0.59‰. The trend in carbonate-δ^44/40Ca in the top six horizons resembles that of sulfate-δ^44/40Ca, but with values 0.1-0.6‰ higher. The range of observed Ca isotope values in this soil is about half that of δ^44/40Ca values observed on Earth. Linear correlation among δ^44/40Ca, δ³⁴S and δ¹⁸O values indicates either (a) a simultaneous change in atmospheric input values for all three elements over time, or (b) isotopic fractionation of all three elements during downward transport. We present evidence that the latter is the primary cause of the isotopic variation that we observe. Sulfate-δ³⁴S values are positively correlated with sulfate-δ¹⁸O values (R² = 0.78) and negatively correlated with sulfate δ^44/40Ca_E values (R² = 0.70). If constant fractionation and conservation of mass with downward transport are assumed, these relationships indicate a δ^44/40Ca fractionation factor of −0.4‰ in CaSO₄. The overall depth trend in Ca isotopes is reproduced by a model of isotopic fractionation during downward Ca transport that considers small and infrequent but regularly recurring rainfall events. Near surface low Ca isotope values are reproduced by a Rayleigh model derived from measured Ca concentrations and the Ca fractionation factor predicted by the relationship with S isotopes. This indicates that the primary mechanism of stable isotope fractionation in CaSO₄ is incremental and effectively irreversible removal of an isotopically enriched dissolved phase by downward transport during small rainfall events.     

Oxygen isotope compositions of chondrules in CR chondrites

We report in situ ion microprobe analyses of oxygen isotopic compositions of olivine, low-Ca pyroxene, high-Ca pyroxene, anorthitic plagioclase, glassy mesostasis, and spinel in five aluminum-rich chondrules and nine ferromagnesian chondrules from the CR carbonaceous chondrites EET92042, GRA95229, and MAC87320. Ferromagnesian chondrules are isotopically homogeneous within ±2‰ in Δ¹⁷O; the interchondrule variations in Δ¹⁷O range from 0 to −5‰. Small oxygen isotopic heterogeneities found in two ferromagnesian chondrules are due to the presence of relict olivine grains. In contrast, two out of five aluminum-rich chondrules are isotopically heterogeneous with Δ¹⁷O values ranging from −6 to −15‰ and from −2 to −11‰, respectively. This isotopic heterogeneity is due to the presence of ¹⁶O-enriched spinel and anorthite (Δ¹⁷O = −10 to −15‰), which are relict phases of Ca,Al-rich inclusions (CAIs) incorporated into chondrule precursors and incompletely melted during chondrule formation. These observations and the high abundance of relict CAIs in the aluminum-rich chondrules suggest a close genetic relationship between these objects: aluminum-rich chondrules formed by melting of spinel-anorthite-pyroxene CAIs mixed with ferromagnesian precursors compositionally similar to magnesium-rich (Type I) chondrules. The aluminum-rich chondrules without relict CAIs have oxygen isotopic compositions (Δ¹⁷O = −2 to −8‰) similar to those of ferromagnesian chondrules. In contrast to the aluminum-rich chondrules from ordinary chondrites, those from CRs plot on a three-oxygen isotope diagram along the carbonaceous chondrite anhydrous mineral line and form a continuum with amoeboid olivine aggregates and CAIs from CRs. We conclude that oxygen isotope compositions of chondrules resulted from two processes: homogenization of isotopically heterogeneous materials during chondrule melting and oxygen isotopic exchange between chondrule melt and ¹⁶O-poor nebular gas.     

Stable isotope and petrologic evidence for open-system degassing during the climactic and pre-climactic eruptions of Mt. Mazama, Crater Lake, Oregon

Evaluation of the extent of volatile element recycling in convergent margin volcanism requires delineating likely source(s) of magmatic volatiles through stable isotopic characterization of sulfur, hydrogen and oxygen in erupted tephra with appropriate assessment of modification by degassing. The climactic eruption of Mt. Mazama ejected approximately 50 km³ of rhyodacitic magma into the atmosphere and resulted in formation of a 10-km diameter caldera now occupied by Crater Lake, Oregon (lat. 43°N, long. 122°W). Isotopic compositions of whole-rocks, matrix glasses and minerals from Mt. Mazama climactic, pre-climactic and postcaldera tephra were determined to identify the likely source(s) of H₂O and S. Integration of stable isotopic data with petrologic data from melt inclusions has allowed for estimation of pre-eruptive dissolved volatile concentrations and placed constraints on the extent, conditions and style of degassing.Sulfur isotope analyses of climactic rhyodacitic whole rocks yield δ³⁴S values of 2.8-14.8‰ with corresponding matrix glass values of 2.4-13.2‰. δ³⁴S tends to increase with stratigraphic height through climactic eruptive units, consistent with open-system degassing. Dissolved sulfur concentrations in melt inclusions (MIs) from pre-climactic and climactic rhyodacitic pumices varies from 80 to 330 ppm, with highest concentrations in inclusions with 4.8-5.2 wt% H₂O (by FTIR). Up to 50% of the initial S may have been lost through pre-eruptive degassing at depths of 4-5 km. Ion microprobe analyses of pyrrhotite in climactic rhyodacitic tephra and andesitic scoria indicate a range in δ³⁴S from −0.4‰ to 5.8‰ and from −0.1‰ to 3.5‰, respectively. Initial δ³⁴S values of rhyodacitic and andesitic magmas were likely near the mantle value of 0‰. Hydrogen isotope (δD) and total H₂O analyses of rhyodacitic obsidian (and vitrophyre) from the climactic fall deposit yielded values οf −103 to −53‰ and 0.23-1.74 wt%, respectively. Values of δD and wt% H₂O of obsidian decrease towards the top of the fall deposit. Samples with depleted δD, and mantle δ¹⁸O values, have elevated δ³⁴S values consistent with open-system degassing. These results imply that more mantle-derived sulfur is degassed to the Earth’s atmosphere/hydrosphere through convergent margin volcanism than previously attributed. Magmatic degassing can modify initial isotopic compositions of sulfur by >14‰ (to δ³⁴S values of 14‰ or more here) and hydrogen isotopic compositions by 90‰ (to δD values of −127‰ in this case).     

Stable isotope fractionation during bacterial sulfate reduction is controlled by reoxidation of intermediates

Bacterial sulfate reduction is one of the most important respiration processes in anoxic habitats and is often assessed by analyzing the results of stable isotope fractionation. However, stable isotope fractionation is supposed to be influenced by the reduction rate and other parameters, such as temperature. We studied here the mechanistic basics of observed differences in stable isotope fractionation during bacterial sulfate reduction. Batch experiments with four sulfate-reducing strains (Desulfovibrio desulfuricans, Desulfobacca acetoxidans, Desulfonatronovibrio hydrogenovorans, and strain TRM1) were performed. These microorganisms metabolize different carbon sources (lactate, acetate, formate, and toluene) and showed broad variations in their sulfur isotope enrichment factors. We performed a series of experiments on isotope exchange of ¹⁸O between residual sulfate and ambient water. Batch experiments were conducted with ¹⁸O-enriched (δ¹⁸O_water = +700‰) and depleted water (δ¹⁸O_water = −40‰), respectively, and the stable ¹⁸O isotope shift in the residual sulfate was followed. For Desulfovibrio desulfuricans and Desulfonatronovibrio hydrogenovorans, which are both characterized by low sulfur isotope fractionation (ε_S > −13.2‰), δ¹⁸O values in the remaining sulfate increased by only 50‰ during growth when ¹⁸O-enriched water was used for the growth medium. In contrast, with Desulfobacca acetoxidans and strain TRM1 (ε_S < −22.7‰) the residual sulfate showed an increase of the sulfate δ¹⁸O close to the values of the enriched water of +700‰. In the experiments with δ¹⁸O-depleted water, the oxygen isotope values in the residual sulfate stayed fairly constant for strains Desulfovibrio desulfuricans, Desulfobacca acetoxidans and Desulfonatronovibrio hydrogenovorans. However, strain TRM1, which exhibits the lowest sulfur isotope fractionation factor (ε_S < −38.7‰) showed slightly decreasing δ¹⁸O values.Our results give strong evidence that the oxygen atoms of sulfate exchange with water during sulfate reduction. However, this neither takes place in the sulfate itself nor during formation of APS (adenosine-5′-phosphosulfate), but rather in intermediates of the sulfate reduction pathway. These may in turn be partially reoxidized to form sulfate. This reoxidation leads to an incorporation of oxygen from water into the “recycled” sulfate changing the overall ¹⁸O isotopic composition of the remaining sulfate fraction. Our study shows that such incorporation of ¹⁸O is correlated with the stable isotope enrichment factor for sulfur measured during sulfate reduction. The reoxidation of intermediates of the sulfate reduction pathway does also strongly influence the sulfur stable isotope enrichment factor. This aforesaid reoxidation is probably dependent on the metabolic conversion of the substrate and therefore also influences the stable isotope fractionation factor indirectly in a rate dependent manner. However, this effect is only indirect. The sulfur isotope enrichment factors for the kinetic reactions themselves are probably not rate dependent.