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.     

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.     

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.     

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.     

Oxygen and sulfur isotope systematics of sulfate produced by bacterial and abiotic oxidation of pyrite

To better understand reaction pathways of pyrite oxidation and biogeochemical controls on δ¹⁸O and δ³⁴S values of the generated sulfate in acid mine drainage (AMD) and other natural environments, we conducted a series of pyrite oxidation experiments in the laboratory. Our biological and abiotic experiments were conducted under aerobic conditions by using O₂ as an oxidizing agent and under anaerobic conditions by using dissolved Fe(III)_aq as an oxidant with varying δ¹⁸O_H2O values in the presence and absence of Acidithiobacillus ferrooxidans. In addition, aerobic biological experiments were designed as short- and long-term experiments where the final pH was controlled at ∼2.7 and 2.2, respectively. Due to the slower kinetics of abiotic sulfide oxidation, the aerobic abiotic experiments were only conducted as long term with a final pH of ∼2.7. The δ³⁴S_SO4 values from both the biological and abiotic anaerobic experiments indicated a small but significant sulfur isotope fractionation (∼−0.7‰) in contrast to no significant fractionation observed from any of the aerobic experiments. Relative percentages of the incorporation of water-derived oxygen and dissolved oxygen (O₂) to sulfate were estimated, in addition to the oxygen isotope fractionation between sulfate and water, and dissolved oxygen. As expected, during the biological and abiotic anaerobic experiments all of the sulfate oxygen was derived from water. The percentage incorporation of water-derived oxygen into sulfate during the oxidation experiments by O₂ varied with longer incubation and lower pH, but not due to the presence or absence of bacteria. These percentages were estimated as 85%, 92% and 87% from the short-term biological, long-term biological and abiotic control experiments, respectively. An oxygen isotope fractionation effect between sulfate and water (ε¹⁸O_SO4-H2O) of ∼3.5‰ was determined for the anaerobic (biological and abiotic) experiments. This measured value was then used to estimate the oxygen isotope fractionation effects between sulfate and dissolved oxygen in the aerobic experiments which were −10.0‰, −10.8‰, and −9.8‰ for the short-term biological, long-term biological and abiotic control experiments, respectively. Based on the similarity between δ¹⁸O_SO4 values in the biological and abiotic experiments, it is suggested that δ¹⁸O_SO4 values cannot be used to distinguish biological and abiotic mechanisms of pyrite oxidation. The results presented here suggest that Fe(III)_aq is the primary oxidant for pyrite at pH < 3, even in the presence of dissolved oxygen, and that the main oxygen source of sulfate is water-oxygen under both aerobic and anaerobic conditions.     

Oxygen isotope fractionation between apatite-bound carbonate and water determined from controlled experiments with synthetic apatites precipitated at 10-37 °C

The oxygen isotope fractionation between the structural carbonate of inorganically precipitated hydroxyapatite (HAP) and water was determined in the range 10-37 °C. Values of 1000 ln α() are linearly correlated with inverse temperature (K) according to the following equation: 1000 ln α() = 25.19 (±0.53)·T⁻¹ − 56.47 (±1.81) (R² = 0.998). This fractionation equation has a slightly steeper slope than those already established between calcite and water ( [O’Neil et al., 1969] and [Kim and O’Neil, 1997]) even though measured fractionations are of comparable amplitude in the temperature range of these experimental studies. It is consequently observed that the oxygen isotope fractionation between apatite carbonate and phosphate increases from about 7.5‰ up to 9.1‰ with decreasing temperature from 37 °C to 10 °C. A compilation of δ¹⁸O values of both phosphate and carbonate from modern mammal teeth and bones confirms that both variables are linearly correlated, despite a significant scattering up to 3.5‰, with a slope close to 1 and an intercept corresponding to a 1000 ln α() value of 8.1‰. This apparent fractionation factor is slightly higher or close to the fractionation factor expected to be in the range 7-8‰ at the body temperature of mammals.     

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.     

An integrated sulfur isotope model for Namibian shelf sediments

In this study the sulfur cycle in the organic-rich mud belt underlying the highly productive upwelling waters of the Namibian shelf is quantified using a 1D reaction-transport model. The model calculates vertical concentration and reaction rate profiles in the top 500 cm of sediment which are compared to a comprehensive dataset which includes carbon, sulfur, nitrogen and iron compounds as well as sulfate reduction (SR) rates and stable sulfur isotopes (³²S, ³⁴S). The sulfur dynamics in the well-mixed surface sediments are strongly influenced by the activity of the large sulfur bacteria Thiomargaritanamibiensis which oxidize sulfide (H₂S) to sulfate () using sea water nitrate () as the terminal electron acceptor. Microbial sulfide oxidation (SOx) is highly efficient, and the model predicts intense cycling between and H₂S driven by coupled SR and SOx at rates exceeding 6.0 mol S m⁻² y⁻¹. More than 96% of the SR is supported by SOx, and only 2-3% of the pool diffuses directly into the sediment from the sea water. A fraction of the produced by Thiomargarita is drawn down deeper into the sediment where it is used to oxidize methane anaerobically, thus preventing high methane concentrations close to the sediment surface. Only a small fraction of total H₂S production is trapped as sedimentary sulfide, mainly pyrite (FeS₂) and organic sulfur (S_org) (∼0.3 wt.%), with a sulfur burial efficiency which is amongst the lowest values reported for marine sediments (<1%). Yet, despite intense SR, FeS₂ and S_org show an isotope composition of ∼5 ‰ at 500 cm depth. These heavy values were simulated by assuming that a fraction of the solid phase sulfur exchanges isotopes with the dissolved sulfide pool. An enrichment in H₂S of ³⁴S towards the sediment-water interface suggests that Thiomargarita preferentially remove H₂³²S from the pore water. A fractionation of 20-30‰ was estimated for SOx (ε_SOx) with the model, along with a maximum fractionation for SR (ε_SR-max) of 100‰. These values are far higher than previous laboratory-based estimates for these processes. Mass balance calculations indicate negligible disproportionation of autochthonous elemental sulfur; an explanation routinely cited in the literature to account for the large fractionations in SR. Instead, the model indicates that repeated multi-stepped sulfide oxidation and intracellular disproportionation by Thiomargarita could, in principle, allow the measured isotope data to be simulated using much lower fractionations for ε_SOx (5‰) and ε_SR (78‰).     

Kinetic oxygen isotope effects during dissimilatory sulfate reduction: A combined theoretical and experimental approach

Kinetic isotope effects related to the breaking of chemical bonds drive sulfur isotope fractionation during dissimilatory sulfate reduction (DSR), whereas oxygen isotope fractionation during DSR is dominated by exchange between intercellular sulfur intermediates and water. We use a simplified biochemical model for DSR to explore how a kinetic oxygen isotope effect may be expressed. We then explore these relationships in light of evolving sulfur and oxygen isotope compositions (δ³⁴S_SO4 and δ¹⁸O_SO4) during batch culture growth of twelve strains of sulfate-reducing bacteria. Cultured under conditions to optimize growth and with identical δ¹⁸O_H2O and initial δ¹⁸O_SO4, all strains show ³⁴S enrichment, whereas only six strains show significant ¹⁸O enrichment. The remaining six show no (or minimal) change in δ¹⁸O_SO4 over the growth of the bacteria. We use these experimental and theoretical results to address three questions: (i) which sulfur intermediates exchange oxygen isotopes with water, (ii) what is the kinetic oxygen isotope effect related to the reduction of adenosine phosphosulfate (APS) to sulfite (SO₃²⁻), (iii) does a kinetic oxygen isotope effect impact the apparent oxygen isotope equilibrium values? We conclude that oxygen isotope exchange between water and a sulfur intermediate likely occurs downstream of APS and that our data constrain the kinetic oxygen isotope fractionation for the reduction of APS to sulfite to be smaller than 4‰. This small oxygen isotope effect impacts the apparent oxygen isotope equilibrium as controlled by the extent to which APS reduction is rate-limiting.     

Isotopic fractionations associated with phosphoric acid digestion of carbonate minerals: Insights from first-principles theoretical modeling and clumped isotope measurements

Phosphoric acid digestion has been used for oxygen- and carbon-isotope analysis of carbonate minerals since 1950, and was recently established as a method for carbonate ‘clumped isotope’ analysis. The CO₂ recovered from this reaction has an oxygen isotope composition substantially different from reactant carbonate, by an amount that varies with temperature of reaction and carbonate chemistry. Here, we present a theoretical model of the kinetic isotope effects associated with phosphoric acid digestion of carbonates, based on structural arguments that the key step in the reaction is disproportionation of H₂CO₃ reaction intermediary. We test that model against previous experimental constraints on the magnitudes and temperature dependences of these oxygen isotope fractionations, and against new experimental determinations of the fractionation of ¹³C-¹⁸O-containing isotopologues (‘clumped’ isotopic species). Our model predicts that the isotope fractionations associated with phosphoric acid digestion of carbonates at 25 °C are 10.72‰, 0.220‰, 0.137‰, 0.593‰ for, respectively, ¹⁸O/¹⁶O ratios (1000 lnα^∗) and three indices that measure proportions of multiply-substituted isotopologues . We also predict that oxygen isotope fractionations follow the mass dependence exponent, λ of 0.5281 (where ). These predictions compare favorably to independent experimental constraints for phosphoric acid digestion of calcite, including our new data for fractionations of ¹³C-¹⁸O bonds (the measured change in Δ₄₇ = 0.23‰) during phosphoric acid digestion of calcite at 25 °C.We have also attempted to evaluate the effect of carbonate cation compositions on phosphoric acid digestion fractionations using cluster models in which disproportionating H₂CO₃ interacts with adjacent cations. These models underestimate the magnitude of isotope fractionations and so must be regarded as unsucsessful, but do reproduce the general trend of variations and temperature dependences of oxygen isotope acid digestion fractionations among different carbonate minerals. We suggest these results present a useful starting point for future, more sophisticated models of the reacting carbonate/acid interface. Examinations of these theoretical predictions and available experimental data suggest cation radius is the most important factor governing the variations of isotope fractionation among different carbonate minerals. We predict a negative correlation between acid digestion fractionation of oxygen isotopes and of ¹³C-¹⁸O doubly-substituted isotopologues, and use this relationship to estimate the acid digestion fractionation of for different carbonate minerals. Combined with previous theoretical evaluations of ¹³C-¹⁸O clumping effects in carbonate minerals, this enables us to predict the temperature calibration relationship for different carbonate clumped isotope thermometers (witherite, calcite, aragonite, dolomite and magnesite), and to compare these predictions with available experimental determinations. The success of our models in capturing several of the features of isotope fractionation during acid digestion supports our hypothesis that phosphoric acid digestion of carbonate minerals involves disproportionation of transition state structures containing H₂CO₃.     

A revised isotope fractionation model for dissimilatory sulfate reduction in sulfate reducing bacteria

Sulfur isotope fractionation during dissimilatory sulfate reduction has been conceptually described by the widely accepted Rees model as related to the stepwise reduction of sulfate to sulfide within the cells of bacteria. The magnitude of isotope fractionation is determined by the interplay between different reduction steps in a chain of reactions. Here we present a revision of Rees’ model for bacterial sulfate reduction that includes revised fractionation factors for the sulfite-sulfide step and incorporates new forward and reverse steps in the reduction of sulfite to sulfide, as well as exchange of sulfide between the cell and ambient water. With this model we show that in contrast to the Rees model, isotope fractionations well in excess of −46‰ are possible. Therefore, some of the large sulfur isotope fractionations observed in nature can be explained without the need of alternate pathways involving the oxidative sulfur cycle. We use this model to predict that large fractionations should occur under hypersulfidic conditions and where electron acceptor concentrations are limiting.     

Seasonal δS variations in two high elevation snow pits measured by S-S double spike thermal ionization mass spectrometry

δ³⁴S and sulfate concentrations were determined in snow pit samples using a thermal ionization mass spectrometric technique capable of 0.2‰ accuracy and requires ≈5 μg (0.16 μmol) natural S. The technique utilizes a ³³S-³⁶S double spike for instrumental mass fractionation correction, and has been applied to snow pit samples collected from the Inilchek Glacier, Kyrgyzstan and from Summit, Greenland. These δ³⁴S determinations provide the first high-resolution seasonal data for these sites, and are used to estimate seasonal sulfate sources. Deuterium (δD) and oxygen (δ¹⁸O) isotope data show that the Inilchek and Summit snow pit samples represent precipitation over ≈20 months.The δ³⁴S values for the Inilchek ranged from +2.6 ± 0.4‰ to +7.6 ± 0.4‰ on sample sizes ranging from 0.3 to 1.8 μmol S. δ³⁴S values for Greenland ranged from +3.6 ± 0.7‰ to +13.3 ± 5‰ for sample sizes ranging from 0.05 to 0.29 μmol S. The concentration ranged from 92.6 ± 0.4 to 1049 ± 4 ng/g for the Inilchek and 18 ± 9 to 93 ± 6 ng/g for the Greenland snow pit. Anthropogenic sulfate dominates throughout the sampled time interval for both sites based on mass balance considerations. Additionally, both sites exhibit a seasonal signature in both δ³⁴S and concentration. The thermal ionization mass spectrometric technique has three advantages compared to gas source isotopic methods: (1) sample size requirements of this technique are 10-fold less permitting access to the higher resolution S isotope record of low concentration snow and ice, (2) the double spike technique permits δ³⁴S and S concentration to be determined simultaneously, and (3) the double spike is an internal standard.     

Iron isotope fractionation and atom exchange during sorption of ferrous iron to mineral surfaces

The application of stable Fe isotopes as a tracer of the biogeochemical Fe cycle necessitates a mechanistic knowledge of natural fractionation processes. We studied the equilibrium Fe isotope fractionation upon sorption of Fe(II) to aluminum oxide (γ-Al₂O₃), goethite (α-FeOOH), quartz (α-SiO₂), and goethite-loaded quartz in batch experiments, and performed continuous-flow column experiments to study the extent of equilibrium and kinetic Fe isotope fractionation during reactive transport of Fe(II) through pure and goethite-loaded quartz sand. In addition, batch and column experiments were used to quantify the coupled electron transfer-atom exchange between dissolved Fe(II) (Fe(II)_aq) and structural Fe(III) of goethite. All experiments were conducted under strictly anoxic conditions at pH 7.2 in 20 mM MOPS (3-(N-morpholino)-propanesulfonic acid) buffer and 23 °C. Iron isotope ratios were measured by high-resolution MC-ICP-MS. Isotope data were analyzed with isotope fractionation models. In batch systems, we observed significant Fe isotope fractionation upon equilibrium sorption of Fe(II) to all sorbents tested, except for aluminum oxide. The equilibrium enrichment factor, , of the Fe(II)_sorb-Fe(II)_aq couple was 0.85 ± 0.10‰ (±2σ) for quartz and 0.85 ± 0.08‰ (±2σ) for goethite-loaded quartz. In the goethite system, the sorption-induced isotope fractionation was superimposed by atom exchange, leading to a δ^56/54Fe shift in solution towards the isotopic composition of the goethite. Without consideration of atom exchange, the equilibrium enrichment factor was 2.01 ± 0.08‰ (±2σ), but decreased to 0.73 ± 0.24‰ (±2σ) when atom exchange was taken into account. The amount of structural Fe in goethite that equilibrated isotopically with Fe(II)_aq via atom exchange was equivalent to one atomic Fe layer of the mineral surface (∼3% of goethite-Fe). Column experiments showed significant Fe isotope fractionation with δ^56/54Fe(II)_aq spanning a range of 1.00‰ and 1.65‰ for pure and goethite-loaded quartz, respectively. Reactive transport of Fe(II) under non-steady state conditions led to complex, non-monotonous Fe isotope trends that could be explained by a combination of kinetic and equilibrium isotope enrichment factors. Our results demonstrate that in abiotic anoxic systems with near-neutral pH, sorption of Fe(II) to mineral surfaces, even to supposedly non-reactive minerals such as quartz, induces significant Fe isotope fractionation. Therefore we expect Fe isotope signatures in natural systems with changing concentration gradients of Fe(II)_aq to be affected by sorption.     

Biogeochemistry of sulfur and iron in Thioploca-colonized surface sediments in the upwelling area off central chile

The biogeochemistry of sedimentary sulfur was investigated on the continental shelf off central Chile at water depths between 24 and 88 m under partial influence of an oxygen minimum zone. Dissolved and solid iron and sulfur species, including the sulfur intermediates sulfite, thiosulfate, and elemental sulfur, were analyzed at high resolution in the top 20 cm. All stations were characterized by high rates of sulfate reduction, but only the sediments within the Bay of Concepción contained dissolved sulfide. Due to advection and/or in-situ reoxidation of sulfide, dissolved sulfate was close to bottom water values. Whereas the concentrations of sulfite and thiosulfate were mostly in the submicromolar range, elemental sulfur was by far the dominant sulfur intermediate. Although the large nitrate- and sulfur-storing bacteria Thioploca were abundant, the major part of S⁰ was located extracellularly. The distribution of sulfur species and dissolved iron suggests the reaction of sulfide with FeOOH as an important pathway for sulfide oxidation and sulfur intermediate formation. This is in agreement with the sulfur isotope composition of co-existing elemental sulfur and iron monosulfides. In the Bay of Concepción, sulfur isotope data suggest that pyrite formation proceeds via the reaction of FeS with polysulfides or H₂S. At the shelf stations, on the other hand, pyrite was significantly depleted in ³⁴S relative to its potential precursors FeS and S⁰. Isotope mass balance considerations suggest further that pyritization at depth includes light sulfide, potentially originating from bacterial sulfur disproportionation. The δ³⁴S-values of pyrite down to −38‰ vs. V-CDT are among the lightest found in organic-rich marine sediments. Seasonal variations in the sulfur isotope composition of dissolved sulfate indicated a dynamic non-steady-state sulfur cycle in the surface sediments. The ¹⁸O content of porewater sulfate increased with depth at all sites compared to the bottom water composition due to intracellular isotope exchange reactions during microbial sulfur transformations.     

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)

Microbially mediated sulfate reduction affects the isotopic composition of dissolved and solid sulfur species in marine sediments. Experiments and field data show that the composition is also modified in the presence of sulfate-reducing microorganisms. This has been attributed either to a kinetic isotope effect during the reduction of sulfate to sulfite, cell-internal exchange reactions between enzymatically-activated sulfate (APS), and/or sulfite with cytoplasmic water. The isotopic fingerprint of these processes may be further modified by the cell-external reoxidation of sulfide to elemental sulfur, and the subsequent disproportionation to sulfide and sulfate or by the oxidation of sulfite to sulfate. Here we report values from interstitial water samples of ODP Leg 182 (Site 1130) and provide the mathematical framework to describe the oxygen isotope fractionation of sulfate during microbial sulfate reduction. We show that a purely kinetic model is unable to explain our data, and that the data are well explained by a model using oxygen isotope exchange reactions. We propose that the oxygen isotope exchange occurs between APS and cytoplasmic water, and/or between sulfite and adenosine monophosphate (AMP) during APS formation. Model calculations show that cell external reoxidation of reduced sulfur species would require up to 3000 mol/m³ of an oxidant at ODP Site 1130, which is incompatible with the sediment geochemical data. In addition, we show that the volumetric fluxes required to explain the observed data are on average 14 times higher than the volumetric sulfate reduction rates (SRR) obtained from inverse modeling of the porewater data. The ratio between the gross sulfate flux into the microbes and the net sulfate flux through the microbes is depth invariant, and independent of sulfide concentrations. This suggests that both fluxes are controlled by cell density and that cell-specific sulfate reduction rates remain constant with depth.     

Experimental investigation of the effects of temperature and ionic strength on Mo isotope fractionation during adsorption to manganese oxides

The Mo stable isotope system is being applied to study changes in ocean redox. Such applications implicitly assume that Mo isotope fractionation in aqueous systems is relatively insensitive to frequently changing environmental variables such as temperature (T) and ionic strength (I). A major driver of fractionation is the adsorption of Mo to Mn oxyhydroxide surfaces [Barling J. and Anbar A. D. (2004) Molybdenum isotope fractionation during adsorption by manganese oxides. Earth Planet. Sci. Lett.217(3-4), 315-329]. Here, we report the results of experiments that determine the extent to which Mo isotope fractionation during adsorption of Mo to the Mn oxyhydroxide mineral birnessite is sensitive to T and I. The results are compared to new predictions from quantum chemical computations. We measured fractionation from 1 to 50 °C at I = 0.1 m and found that Δ^97/95Mo_{dissolved-adsorbed} varies from 1.9‰ to 1.6‰ over this temperature range. Experiments were also performed at 25 °C in synthetic seawater (I = 0.7); fractionation at this condition was the same within analytical error as in low ionic strength experiments. These findings confirm that the Mo isotope fractionation during adsorption to Mn oxyhydroxides is relatively insensitive to variations and T and I over environmentally relevant ranges. To relate these findings to potential mechanisms of Mo isotope fractionation, we also report results for density functional theory computations of the fractionation between and various possible structures of molybdic acid as a function of temperature. Because no plausible species fractionates from with a magnitude matching the experiments, we are left with three possibilities to explain the fractionation: (1) solvation effects on the vibrational frequencies of aqueous species considered thus far are significant, such that our calculations in vacuo yield inaccurate fractionations; (2) a trace aqueous species not yet considered fractionates from and then adsorbs to birnessite; or (3) a surface complex not present in solution forms on birnessite in which Mo is not tetrahedrally coordinated. Our findings help validate assumptions underlying paleoceanographic applications of the Mo isotope system and also lead us closer to understanding the mechanism of isotope fractionation during adsorption of Mo to Mn oxyhydroxides.     

Modeling sulfur isotope fractionation and differential diffusion during sulfate reduction in sediments of the Cariaco Basin

Sulfur isotope composition (δ³⁴S) profiles in sediment pore waters often show an offset between sulfate and sulfide much greater in magnitude than S isotope fractionations observed in pure cultures. A number of workers have invoked an additional reaction, microbial disproportionation of sulfur intermediates, to explain the offset between experimental and natural systems. Here, we present an alternative explanation based on modeling of pore water sulfate and sulfide concentrations and stable isotope data from the Cariaco Basin (ODP Leg 165, Site 1002B). The use of unique diffusion coefficients for and , based on their unequal molecular masses, resulted in an increase in the computed fractionation by almost 10‰, when compared to the common assumption of equal diffusion coefficients for the two species. These small differences in diffusion coefficients yield calculated isotopic offsets between coeval sediment pore water sulfate and sulfide without disproportionation (up to 53.4‰) that exceed the largest fractionations observed in experimental cultures. Furthermore, the diffusion of sulfide within sediment pore waters leads to values that are even greater than those predicted by our model for sulfate reduction with unique diffusion coefficients. These diffusive effects on the sulfur isotope composition of pore water sulfate and sulfide can impact our interpretations of geologic records of sulfate and sulfide minerals, and should be considered in future studies.     

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.     

Modeling the effects of bond environment on equilibrium iron isotope fractionation in ferric aquo-chloro complexes

Four or five sets of ab initio models, including Unrestricted Hartree Fock (UHF) and hybrid Density Functional Theory (DFT) are calculated for each species in a series of aqueous ferric aquo-chloro complexes: , , , FeCl₃(H₂O)₃, FeCl₃(H₂O)₂, , FeCl₅H₂O²⁻, , ) in order to determine the relative isotopic fractionation among the complexes, to compare the results of different models for the same complexes, to examine factors that influence the magnitude of the isotopic fractionation, and to compare bond-partner-driven fractionation with redox-driven fractionation.Relative to , all models show a nearly linear decrease in ⁵⁶Fe/⁵⁴Fe as the number of Cl⁻ ions per Fe³⁺ ion increases, with slopes of −0.8‰ to −1.0‰ per Cl⁻ at 20 °C. At 20 °C, 1000 ln β (β = ⁵⁶Fe/⁵⁴Fe reduced partition function ratio relative to a dissociated Fe atom) values range from 8.93‰ to 9.73‰ for , 8.04-9.12‰ for , 7.61-8.73‰ for , 7.14-8.25‰ for , and 3.09-4.41‰ for . The fractionation between and ranges from 1.5‰ to 2.6‰, depending on the model; this is comparable in magnitude to fractionation effects due to Fe³⁺/Fe²⁺ redox reactions. β values from the UHF models are consistently higher than those from the hybrid DFT models.Isotopic fractionation is shown to be sensitive to differences in ligand bond stiffness (above), coordination number, bond length, and the frequency of the asymmetric Fe-X stretching vibrational mode, as predicted by previous theoretical studies. Complexes with smaller coordination numbers have higher 1000 ln β (7.46‰, 5.25‰, and 3.48‰ for , ,, respectively, from the B3LYP/6-31G(d) model). Species with the same number of chlorides but fewer waters also show the effect of coordination number on 1000 ln β: (7.46‰ vs. 7.05‰ for FeCl₃(H₂O)₂ vs. FeCl₃(H₂O)₃ and 5.25‰ vs. 4.94‰ for vs. FeCl₅H₂O²⁻ with the B3LYP/6-31G(d) model). As more Fe-Cl bonds substitute for Fe-OH₂ bonds (with a resulting decrease in β), the lengths of the Fe-Cl bonds and the Fe-O bonds increase.Preliminary modeling of shows an Fe³⁺/Fe²⁺ fractionation of 3.2‰ for the B3LYP/6-31G(d) model, in agreement with previous studies. The addition of an explicit outer hydration sphere of 12 H₂O molecules to models of improves agreement with measured vibrational frequencies and bond lengths; 1000 ln β increases by 0.8-1.0‰. An additional hydration sphere around increases 1000 ln β by only 0.1‰.Isotopic fractionations predicted for this simple system imply that ligands present in an aqueous iron environment are potentially important drivers of fractionation, and suggest that significant fractionation effects are likely in other aqueous systems containing sulfides or organic ligands. Fractionation effects due to both speciation and redox must be considered when interpreting iron isotope fractionations in the geological record.     

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.