Quadruple sulfur isotope analysis of ca. 3.5 Ga Dresser Formation: New evidence for microbial sulfate reduction in the early Archean

Multiple sulfur isotope system is a powerful new tracer for atmospheric, volcanic, and biological influences on sulfur cycles in the anoxic early Earth. Here, we report high-precision quadruple sulfur isotope analyses (³²S/³³S/³⁴S/³⁶S) of barite, pyrite in barite, and sulfides in related hydrothermal and igneous rocks occurring in the ca. 3.5 Ga Dresser Formation, Western Australia. Our results indicate that observed isotopic variations are mainly controlled by mixing of mass-dependently (MD) and non-mass-dependently fractionated (non-MD) sulfur reservoirs. Based on the quadruple sulfur isotope systematics (δ³⁴S-Δ³³S-Δ³⁶S) for these minerals, four end-member sulfur reservoirs have been recognized: (1) non-MD sulfate (δ³⁴S = −5 ± 2‰; Δ³³S = −3 ± 1‰); (2) MD sulfate (δ³⁴S = +10 ± 3‰); (3) non-MD sulfur (δ³⁴S > +6‰; Δ³³S > +4‰); and (4) igneous MD sulfur (δ³⁴S = Δ³³S = 0‰). The first and third components show a clear non-MD signatures, thus probably represent sulfate and sulfur aerosol inputs. The MD sulfate component (2) is enriched in ³⁴S (+10 ± 3‰) and may have originated from microbial and/or abiotic disproportionation of volcanic S or SO₂. Our results reconfirm that the Dresser barites contain small amounts of pyrite depleted in ³⁴S by 15-22‰ relative to the host barite. These barite-pyrite pairs exhibit a mass-dependent relationship of δ³³S/δ³⁴S with slope less than 0.512, which is consistent with that expected for microbial sulfate reduction and is significantly different from that of equilibrium fractionation (0.515). The barite-pyrite pairs also show up to 1‰ difference in Δ³⁶S values and steep Δ³⁶S/Δ³³S slopes, which deviate from the main Archean array (Δ³⁶S/Δ³³S = −0.9) and are comparable to isotope effects exhibited by sulfate reducing microbes (Δ³⁶S/Δ³³S = −5 to −11). These new lines of evidence support the existence of sulfate reducers at ca. 3.5 Ga, whereas microbial sulfur disproportionation may have been more limited than recently suggested.     

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.     

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).     

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 and carbon isotope records from 1700 to 800 Ma carbonates of the Jixian section, northern China: Implications for secular isotope variations in Proterozoic seawater and relationships to global supercontinental events

This study is a comprehensive, stable isotope survey of the marine carbonate-dominated, upper Paleo- to lower Neoproterozoic stratigraphy of Jixian County, China. Carbonate-associated sulfate (CAS) was extracted and measured for δ³⁴S_CAS using the same samples analyzed for δ¹³C_carbonate. This integrated proxy approach is a step towards a more comprehensive picture of secular variation in the composition of Proterozoic seawater. We specifically sampled marine carbonate intervals from the lower section of the Chuanlinggou Formation, Changcheng Group (ca. 1700 Ma) to the top of the Jingeryu Formation, Qingbaikou Group (ca. 800 Ma). δ¹³C_carbonate values are mostly negative in the upper Paleoproterozoic Changcheng Group, with an ascending trend from −3‰ to 0‰. We observed variation of approximately 0 ± 1‰ in the Mesoproterozoic Jixian Group, and positive values of +2 ± 2‰ characterize the lower Neoproterozoic Qingbaikou Group. Stratigraphic variations in δ³⁴S_CAS are more remarkable in their ranges and magnitudes, including conspicuously high values exceeding +30‰ in the three intervals at ca. 1700 Ma, 1300-1100 Ma, and 1000-900 Ma. In the Changcheng Group, δ³⁴S_CAS values are typically higher than +25‰, with only a few values of less than +15‰. In contrast, most of the data spanning from the Mesoproterozoic Tieling Formation of the Jixian Group to the lower Neoproterozoic Jingeryu Formation of the Qingbaikou Group are highly variable between +10‰ and +25‰, with some values exceeding +25‰.In the late Paleoproterozoic (1700-1600 Ma), a >10‰ decrease in δ³⁴S_CAS and ∼3‰ increase in δ¹³C_carbonate are coincident with, and likely related to, the breakup of Columbia, a supercontinent that predated Rodinia. Carbon and sulfur isotope data from the Mesoproterozoic, when global tectonic activity was comparatively weaker, fall mostly in the ranges of +15 ± 10‰ and 0 ± 1‰, respectively, but fluctuations of >20‰ for δ³⁴S_CAS and >3‰ for the δ¹³C_carbonate at ca. 1450-1400 Ma may reflect subduction and large-scale magmatic activity in island arcs marking the end of Columbia breakup. From the late Mesoproterozoic (ca. 1300-1100 Ma) to the early Neoproterozoic (ca. 800 Ma), the δ¹³C and δ³⁴S of seawater increased gradually with increasing variability. Most impressive areδ³⁴S_CAS values that exceed +30‰ in two intervals at ca. 1300-1100 Ma and ca. 1000-900 Ma, which may reflect the assembly and early breakup of Rodinia. Although gaps in the record remain, and studies of even higher resolution are warranted, our results suggest that changes in paleoceanographic conditions linked to global tectonics strongly influenced the biogeochemical cycles of C and S. Furthermore, periods of the Proterozoic previously noted for their isotopic invariability show clear isotopic expressions of this tectonic activity.     

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.     

The triple isotopic composition of oxygen in leaf water

The isotopic composition of atmospheric O₂ depends on the rates of oxygen cycling in photosynthesis, respiration, photochemical reactions in the stratosphere and on δ¹⁷O and δ¹⁸O of ocean and leaf water. While most of the factors affecting δ¹⁷O and δ¹⁸O of air O₂ have been studied extensively in recent years, δ¹⁷O of leaf water—the substrate for all terrestrial photosynthesis—remained unknown. In order to understand the isotopic composition of atmospheric O₂ at present and in fossil air in ice cores, we studied leaf water in field experiments in Israel and in a European survey. We measured the difference in δ¹⁷O and δ¹⁸O between stem and leaf water, which is the result of isotope enrichment during transpiration. We calculated the slopes of the lines linking the isotopic compositions of stem and leaf water. The obtained slopes in ln(δ¹⁷O + 1) vs. ln(δ¹⁸O + 1) plots are characterized by very high precision (∼0.001) despite of relatively large differences between duplicates in both δ¹⁷O and δ¹⁸O (0.02-0.05‰). This is so because the errors in δ¹⁸O and δ¹⁷O are mass-dependent. The slope of the leaf transpiration process varied between 0.5111 ± 0.0013 and 0.5204 ± 0.0005, which is considerably smaller than the slope linking liquid water and vapor at equilibrium (0.529). We further found that the slope of the transpiration process decreases with atmospheric relative humidity (h) as 0.522-0.008 × h, for h in the range 0.3-1. This slope is neither influenced by the plant species, nor by the environmental conditions where plants grow nor does it show strong variations along long leaves.     

Calibration of the calcite-water oxygen-isotope geothermometer at Devils Hole, Nevada, a natural laboratory

The δ¹⁸O of ground water (−13.54 ± 0.05 ‰) and inorganically precipitated Holocene vein calcite (+14.56 ± 0.03 ‰) from Devils Hole cave #2 in southcentral Nevada yield an oxygen isotopic fractionation factor between calcite and water at 33.7 °C of 1.02849 ± 0.00013 (1000 ln α_{calcite-water} = 28.09 ± 0.13). Using the commonly accepted value of ∂(α_{calcite-water})/∂T of −0.00020 K⁻¹, this corresponds to a 1000 ln α_{calcite-water} value at 25 °C of 29.80, which differs substantially from the current accepted value of 28.3. Use of previously published oxygen isotopic fractionation factors would yield a calcite precipitation temperature in Devils Hole that is 8 °C lower than the measured ground water temperature. Alternatively, previously published fractionation factors would yield a δ¹⁸O of water, from which the calcite precipitated, that is too negative by 1.5 ‰ using a temperature of 33.7 °C. Several lines of evidence indicate that the geochemical environment of Devils Hole has been remarkably constant for at least 10 ka. Accordingly, a re-evaluation of calcite-water oxygen isotopic fractionation factor may be in order.Assuming the Devils Hole oxygen isotopic value of α_{calcite-water} represents thermodynamic equilibrium, many marine carbonates are precipitated with a δ¹⁸O value that is too low, apparently due to a kinetic isotopic fractionation that preferentially enriches ¹⁶O in the solid carbonate over ¹⁸O, feigning oxygen isotopic equilibrium.     

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.     

A paired sulfate-pyrite δS approach to understanding the evolution of the Ediacaran-Cambrian sulfur cycle

An anomalous enrichment in marine sulfate δ³⁴S_SO4 is preserved in globally-distributed latest Ediacaran-early Cambrian strata. The proximity of this anomaly to the Ediacaran-Cambrian boundary and the associated evolutionary radiation has invited speculation that the two are causally related. Here we present a high-resolution record of paired sulfate (δ³⁴S_SO4) and pyrite (δ³⁴S_pyr) from sediments spanning ca. 547-540 million years ago (Ma) from the Ara Group of the Huqf Supergroup, Sultanate of Oman. We observe an increase in δ³⁴S_SO4 from ∼20‰ to ∼42‰, beginning at ca. 550 Ma and continuing at least through ca. 540 Ma. There is a concomitant increase in δ³⁴S_pyr over this interval from ∼ −15‰ to 10‰. This globally correlative enrichment, here termed the Ara anomaly, constitutes a major perturbation to the sulfur cycle. The absolute values of δ³⁴S_pyr reported here and in equivalent sections around the world, require the isotopic composition of material entering the ocean (δ³⁴S_in) to be significantly more enriched than modern (∼3‰) values, likely in excess of 12‰ during the late Ediacaran-early Cambrian. Against this background of elevated δ³⁴S_in, the Ara anomaly is explained not by increased fractionation between sulfate and pyrite (Δδ³⁴S), but by an increase in pyrite burial (f_pyr), most likely driven by enhanced primary production and sequestration of organic carbon, consistent with earlier reports of elevated organic carbon burial and widespread phosphorite deposition.     

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.     

Sulfur and oxygen isotope study of sulfate reduction in experiments with natural populations from Fællestrand, Denmark

This study investigates the sulfur and oxygen isotope fractionations of dissimilatory sulfate reduction and works to reconcile the relationships between the oxygen and sulfur isotopic and elemental systems. We report results of experiments with natural populations of sulfate-reducing bacteria using sediment and seawater from a marine lagoon at Fællestrand on the northern shore of the island of Fyn, Denmark. The experiments yielded relatively large magnitude sulfur isotope fractionations for dissimilatory sulfate reduction (up to approximately 45‰ for ³⁴S/³²S) with higher δ¹⁸O accompanying higher δ³⁴S, similar to that observed in previous studies. The seawater used in the experiments was spiked by addition of ¹⁷O-labeled water and the ¹⁷O content of residual sulfate was found to depend on the fraction of sulfate reduced in the experiments. The ¹⁷O data provides evidence for recycling of sulfur from metabolic intermediates and for an ¹⁸O/¹⁶O fractionation of ∼25-30‰ for dissimilatory sulfate reduction. The close correlation between the ¹⁷O data and the sulfur isotope data suggests that isotopic exchange between cell water and external water (reactor water) was rapid under experimental conditions. The molar ratio of oxygen exchange to sulfate reduction was found to be about 2.5. This value is slightly lower than observed in studies of natural ecosystems [e.g., Wortmann U. G., Chernyavsky B., Bernasconi S. M., Brunner B., Böttcher M. E. and Swart P. K. (2007) Oxygen isotope biogeochemistry of pore water sulfate in the deep biosphere: dominance of isotope exchange reactions with ambient water during microbial sulfate reduction (ODP Site 1130). Geochim. Cosmochim. Acta71, 4221-4232]. Using recent models of sulfur isotope fractionations we find that our combined sulfur and oxygen isotopic data places constraints on the proportion of sulfate recycled to the medium (78-96%), the proportion of sulfur intermediate sulfite that was recycled by way of APS to sulfate and released back to the external sulfate pool (∼70%), and also that a fraction of the sulfur intermediates between sulfite and sulfide were recycled to sulfate. These parameters can be constrained because of the independent information provided by δ¹⁸O, δ³⁴S, δ¹⁷O labels, and Δ³³S.     

Oxygen isotope heterogeneity and disequilibria of olivine crystals in large volume Holocene basalts from Iceland: Evidence for magmatic digestion and erosion of Pleistocene hyaloclastites

This work considers petrogenesis of the largest Holocene basaltic fissure eruptions of Iceland, which are also the largest in the world: Laki (1783-84 AD, 15 km³), Eldgjá (934 AD, 18 km³), Veidivötn (900, 1480 AD, multiple eruptions, >2 km³), Núpahraun (ca. 4000 BP, >1 km³) and Thjórsárhraun (ca 8000 BP, >20 km³). We present oxygen isotope laser fluorination analyses of 55 individual and bulk olivine crystals, coexisting individual and bulk plagioclase phenocrysts, and their host basaltic glasses with average precision of better than 0.1‰ (1SD). We also report O isotope analyses of cores and rims of 61 olivine crystals by SIMS with average precision on single spots of 0.24‰ (1SD) in 13 samples coupled with electron microprobe data for major and trace elements in these olivines. Within each individual sample, we have found that basaltic glass is relatively homogeneous with respect to oxygen isotopes, plagioclase phenocrysts exhibit crystal to crystal variability, while individual olivines span from the values in equilibrium with the low-δ¹⁸O matrix glass to those being three permil higher in δ¹⁸O than the equilibrium. Olivine cores with maximum value of 5.2‰ are found in many of these basalts and suggest that the initial magma was equilibrated with normal-δ¹⁸O mantle. No olivines or their intracrystalline domains are found with bulk or spot value higher than those found in MORB olivines. The δ¹⁸O variability of 0.3-3‰ exists for olivine grains from different lavas, and variable core-to-rim oxygen isotopic zoning is present in selected olivine grains. Many olivines in the same sample are not zoned, while a few grains are zoned with respect to oxygen isotopes and exhibit small core-to-core variations in Fe-Mg, Ni, Mn, Ca. Grains that are zoned in both Mg# and δ¹⁸O exhibit positive correlation of these two parameters. Electron microprobe analysis shows that most olivines equilibrated with the transporting melt, and thin Fe-richer rim is present around many grains, regardless of the degree of olivine-melt oxygen isotope disequilibrium.The preservation of isotopic and compositional zoning in selected grains, and subtle to severe Δ¹⁸O (melt-olivine) and Δ¹⁸O (plagioclase-olivine) disequilibria suggests rather short crystal residence times of years to centuries. Synglacially-altered upper crustal, tufaceous hyaloclastites of Pleistocene age serve as a viable source for low-δ¹⁸O values in Holocene basalts through assimilation, mechanical and thermal erosion, and devolatilization of stoped blocks. Cumulates formed in response to cooling during assimilation, and xenocrysts derived from hyaloclastites, contribute to the diverse δ¹⁸O crystalline cargo. The magma plumbing systems under each fissure are likely to include a network of interconnected dikes and sills with high magma flow rates that contribute to the efficacy of magmatic erosion of large quantities (10-60% mass) of hyaloclastites required by isotopic mass balance.Olivine diversity and the pervasive lack of phenocryst-melt oxygen isotopic equilibrium suggest that a common approach of analyzing bulk olivine for oxygen isotopes, as a proxy for the basaltic melt or to infer mantle δ¹⁸O value, needs to proceed with caution. The best approach is to analyze olivine crystals individually and demonstrate their equilibrium with matrix.     

Oxygen- and magnesium-isotope compositions of calcium-aluminum-rich inclusions from Rumuruti (R) chondrites

We report oxygen- and magnesium-isotope compositions of Ca,Al-rich inclusions (CAIs) from several Rumuruti (R) chondrites measured in situ using a Cameca ims-1280 ion microprobe. On a three-isotope oxygen diagram, δ¹⁷O vs. δ¹⁸O, compositions of individual minerals in most R CAIs analyzed fall along a slope-1 line. Based on the variations of Δ¹⁷O values (Δ¹⁷O = δ¹⁷O − 0.52 × δ¹⁸O) within individual inclusions, the R CAIs are divided into (i) ¹⁶O-rich (Δ¹⁷O ∼ −23-26‰), (ii) uniformly ¹⁶O-depleted (Δ¹⁷O ∼ −2‰), and (iii) isotopically heterogeneous (Δ¹⁷O ranges from −25‰ to +5‰). One of the hibonite-rich CAIs, H030/L, has an intermediate Δ¹⁷O value of −12‰ and a highly fractionated composition (δ¹⁸O ∼ +47‰). We infer that like most CAIs in other chondrite groups, the R CAIs formed in an ¹⁶O-rich gaseous reservoir. The uniformly ¹⁶O-depleted and isotopically heterogeneous CAIs subsequently experienced oxygen-isotope exchange during remelting in an ¹⁶O-depleted nebular gas, possibly during R chondrite chondrule formation, and/or during fluid-assisted thermal metamorphism on the R chondrite parent asteroid.Three hibonite-bearing CAIs and one spinel-plagioclase-rich inclusion were analyzed for magnesium-isotope compositions. The CAI with the highly fractionated oxygen isotopes, H030/L, shows a resolvable excess of ²⁶Mg (²⁶Mg^∗) corresponding to an initial ²⁶Al/²⁷Al ratio of ∼7 × 10⁻⁷. Three other CAIs show no resolvable excess of ²⁶Mg (²⁶Mg^∗). The absence of ²⁶Mg^∗ in the spinel-plagioclase-rich CAI from a metamorphosed R chondrite NWA 753 (R3.9) could have resulted from metamorphic resetting. Two other hibonite-bearing CAIs occur in the R chondrites (NWA 1476 and NWA 2446), which appear to have experienced only minor degrees of thermal metamorphism. These inclusions could have formed from precursors with lower than canonical ²⁶Al/²⁷Al ratio.     

Copper isotope fractionation in acid mine drainage

We measured the Cu isotopic composition of primary minerals and stream water affected by acid mine drainage in a mineralized watershed (Colorado, USA). The δ⁶⁵Cu values (based on ⁶⁵Cu/⁶³Cu) of enargite (δ⁶⁵Cu = −0.01 ± 0.10‰; 2σ) and chalcopyrite (δ⁶⁵Cu = 0.16 ± 0.10‰) are within the range of reported values for terrestrial primary Cu sulfides (−1‰ < δ⁶⁵Cu < 1‰). These mineral samples show lower δ⁶⁵Cu values than stream waters (1.38‰ ? δ⁶⁵Cu ? 1.69‰). The average isotopic fractionation (Δ_aq-min = δ⁶⁵Cu_aq − δ⁶⁵Cu_min, where the latter is measured on mineral samples from the field system), equals 1.43 ± 0.14‰ and 1.60 ± 0.14‰ for chalcopyrite and enargite, respectively. To interpret this field survey, we leached chalcopyrite and enargite in batch experiments and found that, as in the field, the leachate is enriched in ⁶⁵Cu relative to chalcopyrite (1.37 ± 0.14‰) and enargite (0.98 ± 0.14‰) when microorganisms are absent. Leaching of minerals in the presence of Acidithiobacillus ferrooxidans results in smaller average fractionation in the opposite direction for chalcopyrite (Δ_aq-min^o=-0.57±0.14‰, where min^o refers to the starting mineral) and no apparent fractionation for enargite (Δ_aq-min^o=0.14±0.14‰). Abiotic fractionation is attributed to preferential oxidation of ⁶⁵Cu⁺ at the interface of the isotopically homogeneous mineral and the surface oxidized layer, followed by solubilization. When microorganisms are present, the abiotic fractionation is most likely not seen due to preferential association of ⁶⁵Cu_aq with A. ferrooxidans cells and related precipitates. In the biotic experiments, Cu was observed under TEM to occur in precipitates around bacteria and in intracellular polyphosphate granules. Thus, the values of δ⁶⁵Cu in the field and laboratory systems are presumably determined by the balance of Cu released abiotically and Cu that interacts with cells and related precipitates. Such isotopic signatures resulting from Cu sulfide dissolution should be useful for acid mine drainage remediation and ore prospecting purposes.     

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 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.     

Oxygen isotopic compositions of chondrules from the metal-rich chondrites Isheyevo (CH/CBb), MAC 02675 (CBb) and QUE 94627 (CBb)

It has been recently suggested that (1) CH chondrites and the CB_b/CH-like chondrite Isheyevo contain two populations of chondrules formed by different processes: (i) magnesian non-porphyritic (cryptocrystalline and barred) chondrules, which are similar to those in the CB chondrites and formed in an impact-generated plume of melt and gas resulted from large-scale asteroidal collision, and (ii) porphyritic chondrules formed by melting of solid precursors in the solar nebula. (2) Porphyritic chondrules in Isheyevo and CH chondrites are different from porphyritic chondrules in other carbonaceous chondrites ( [Krot et al., 2005], [Krot et al., 2008a] and [Krot et al., 2008b]). In order to test these hypotheses, we measured in situ oxygen isotopic compositions of porphyritic (magnesian, Type I and ferroan, Type II) and non-porphyritic (magnesian and ferroan cryptocrystalline) chondrules from Isheyevo and CB_b chondrites MAC 02675 and QUE 94627, paired with QUE 94611, using a Cameca ims-1280 ion microprobe.On a three-isotope oxygen diagram (δ¹⁷O vs. δ¹⁸O), compositions of chondrules measured follow approximately slope-1 line. Data for 19 magnesian cryptocrystalline chondrules from Isheyevo, 24 magnesian cryptocrystalline chondrules and 6 magnesian cryptocrystalline silicate inclusions inside chemically-zoned Fe,Ni-metal condensates from CB_b chondrites have nearly identical compositions: Δ¹⁷O = −2.2 ± 0.9‰, −2.3 ± 0.6‰ and −2.2 ± 1.0‰ (2σ), respectively. These observations and isotopically light magnesium compositions of cryptocrystalline magnesian chondrules in CB_b chondrites (Gounelle et al., 2007) are consistent with their single-stage origin, possibly as gas-melt condensates in an impact-generated plume. In contrast, Δ¹⁷O values for 11 Type I and 9 Type II chondrules from Isheyevo range from −5‰ to +4‰ and from −17‰ to +3‰, respectively. In contrast to typical chondrules from carbonaceous chondrites, seven out of 11 Type I chondrules from Isheyevo plot above the terrestrial fractionation line. We conclude that (i) porphyritic chondrules in Isheyevo belong to a unique population of objects, suggesting formation either in a different nebular region or at a different time than chondrules from other carbonaceous chondrites; (ii) Isheyevo, CB and CH chondrites are genetically related meteorites: they contain non-porphyritic chondrules produced during the same highly-energetic event, probably large-scale asteroidal collision; (iii) the differences in mineralogy, petrography, chemical and whole-rock oxygen isotopic compositions between CH and CB chondrites are due to various proportions of the nebular and the impact-produced materials.     

Effect of electron donors on the fractionation of sulfur isotopes by a marine Desulfovibrio sp.

Sulfur isotope effects produced by microbial dissimilatory sulfate reduction are used to reconstruct the coupled cycling of carbon and sulfur through geologic time, to constrain the evolution of sulfur-based metabolisms, and to track the oxygenation of Earth’s surface. In this study, we investigate how the coupling of carbon and sulfur metabolisms in batch and continuous cultures of a recently isolated marine sulfate reducing bacterium DMSS-1, a Desulfovibrio sp., influences the fractionation of sulfur isotopes.DMSS-1 grown in batch culture on seven different electron donors (ethanol, glycerol, fructose, glucose, lactate, malate and pyruvate) fractionates ³⁴S/³²S ratio from 6‰ to 44‰, demonstrating that the fractionations by an actively growing culture of a single incomplete oxidizing sulfate reducing microbe can span almost the entire range of previously reported values in defined cultures. The magnitude of isotope effect correlates well with cell specific sulfate reduction rates (from 0.7 to 26.1 fmol/cell/day). DMSS-1 grown on lactate in continuous culture produces a larger isotope effect (21-37‰) than the lactate-grown batch culture (6‰), indicating that the isotope effect also depends on the supply rate of the electron donor and microbial growth rate. The largest isotope effect in continuous culture is accompanied by measurable changes in cell length and cellular yield that suggest starvation. The use of multiple sulfur isotopes in the model of metabolic fluxes of sulfur shows that the loss of sulfate from the cell and the intracellular reoxidation of reduced sulfur species contribute to the increase in isotope effects in a correlated manner. Isotope fractionations produced during sulfate reduction in the pure culture of DMSS-1 expand the previously reported range of triple sulfur isotope effects (³²S, ³³S, and ³⁴S) by marine sulfate reducing bacteria, implying that microbial sulfur disproportionation may have a smaller ³³S isotopic fingerprint than previously thought.     

Influence of the enzyme dissimilatory sulfite reductase on stable isotope fractionation during sulfate reduction

The stable isotopes of sulfate are often used as a tool to assess bacterial sulfate reduction on the macro scale. However, the mechanisms of stable isotope fractionation of sulfur and oxygen at the enzymatic level are not yet fully understood. In batch experiments with water enriched in ¹⁸O we investigated the effect of different nitrite concentrations on sulfur isotope fractionation by Desulfovibrio desulfuricans.With increasing nitrite concentrations, we found sulfur isotope enrichment factors ranging from −11.2 ± 1.8‰ to −22.5 ± 3.2‰. Furthermore, the δ¹⁸O values in the remaining sulfate increased from approximately 50-120‰ when ¹⁸O-enriched water was supplied. Since ¹⁸O-exchange with ambient water does not take place in sulfate, but rather in intermediates of the sulfate reduction pathway (e.g. ), we suggest that nitrite affects the steady-state concentration and the extent of reoxidation of the metabolic intermediate sulfite to sulfate during sulfate reduction. Given that nitrite is known to inhibit the production of the enzyme dissimilatory sulfite reductase, our results suggest that the activity of the dissimilatory sulfite reductase regulates the kinetic isotope fractionation of sulfur and oxygen during bacterial sulfate reduction. Our novel results also imply that isotope fractionation during bacterial sulfate reduction strongly depends on the cell internal enzymatic regulation rather than on the physico-chemical features of the individual enzymes.