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.     

Particulate sulfur species in the water column of the Cariaco Basin

The biogeochemistry of iron sulfide minerals in the water column of the Cariaco Basin was investigated in November 2007 (non-upwelling season) and May 2008 (upwelling season) as part of the on-going CARIACO (CArbon Retention In A Colored Ocean) time series project. The concentrations of particulate sulfur species, specifically acid volatile sulfur (AVS), greigite, pyrite, and particulate elemental sulfur, were determined at high resolution near the O₂/H₂S interface. In November 2007, AVS was low throughout the water column, with the highest concentration at the depth where sulfide was first detected (260 m) and with a second peak at 500 m. Greigite, pyrite, and particulate elemental sulfur showed distinct concentration maxima near the interface. In May 2008, AVS was not detected in the water column. Maxima for greigite, pyrite, and particulate elemental sulfur were again observed near the interface. We also studied the iron sulfide flux using sediment trap materials collected at the Cariaco station. Pyrite comprised 0.2-0.4% of the total particulate flux in the anoxic water column, with a flux of 0.5-1.6 mg S m⁻² d⁻¹.Consistent with the water column concentration profiles for iron sulfide minerals, the sulfur isotope composition of particulate sulfur found in deep anoxic traps was similar to that of dissolved sulfide near the O₂/H₂S interface. We conclude that pyrite is formed mainly within the redoxcline where sulfur cycling imparts a distinct isotopic signature compared to dissolved sulfide in the deep anoxic water. This conclusion is consistent with our previous study of sulfur species and chemoautotrophic production, which suggests that reaction of sulfide with reactive iron is an important pathway for sulfide oxidation and sulfur intermediate formation near the interface. Pyrite and elemental sulfur distributions favor a pathway of pyrite formation via the reaction of FeS with polysulfides or particulate elemental sulfur near the interface. A comparison of thermodynamic predictions with actual concentration profiles for iron sulfides leads us to argue that microbes may mediate this precipitation.     

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.     

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.     

Stable sulfur and oxygen isotope ratios of the Świętokrzyski National Park spring waters generated by natural and anthropogenic factors (south-central Poland)

This paper presents the results of an isotopic study of spring waters in ?wi?tokrzyski (Holy Cross Mountain) National Park (?NP), south-central Poland. The δ³⁴S_V-CDT and δ¹⁸O_V-SMOW of soluble sulfates (n = 40) varied from 0.5‰ to 18.1‰ and from 3.5‰ to 12.2‰, respectively. The average δ³⁴S values are closely similar to those of rainwater, soils and rocks (comprising scattered pyrite). This suggests that soluble sulfates in the springs originated from mixing of recent and historic deposition, sulfates derived from pyrite oxidation, and CS-mineralization in soils and debris. An additional anthropogenic sulfur input (inorganic fertilizer) occurs in the water of spring S-61 located in the ?wi?tokrzyski National Park buffer zone. The δ¹⁸O_V-SMOW of spring waters (n = 4) were in the range of −10.6‰ to −10.2‰ indicating that they are derived from vadose groundwater in ?NP. This was the first isotope study of spring waters in the national parks of Poland. It enabled the determination of sulfur pathways and discrimination between natural and anthropogenic sources of this element in a relatively pristine area.     

O and S isotopic composition of dissolved and attached oxidation products of pyrite by Acidithiobacillus ferrooxidans: Comparison with abiotic oxidations

The acidophilic iron-oxidizing bacterium, Acidithiobacillus ferrooxidans, plays a part in the pyrite oxidation process and has been widely studied in order to determine the kinetics of the reactions and the isotopic composition of dissolved product sulphates, but the details of the oxidation processes at the surface of pyrite are still poorly known. In this study, oxygen and sulphur isotopic compositions (δ¹⁸O and δ³⁴S) were analyzed for dissolved sulphates and water from experimental aerobic acidic (pH < 2) pyrite oxidation by A. ferrooxidans. The oxidation products attached to the pyrite surfaces were studied for their morphology (SEM), their chemistry (Raman spectroscopy) and for their δ¹⁸O (ion microprobe). They were compared to abiotically (Fe³⁺, H₂O₂, O₂) oxidized pyrite surface compounds in order to constrain the oxidation pathways and to look for the existence of potential biosignatures for this system.The pyrite dissolution evolved from non-stoichiometric (during the first days) to stoichiometric (with increasing time) resulting in dissolved sulphates having distinct δ¹⁸O (e.g. +11.0‰ and −2.0‰, respectively) and δ³⁴S (+4.5‰ and +2.8‰, respectively) values. The “oxidation layer” at the surface of pyrite is complex and made of iron oxides, sulphate, polysulphide, elemental sulphur and polythionates. Bio- and Fe³⁺-oxidation favour the development of monophased micrometric bumps made of hematite or sulphate while other abiotic oxidation processes result in more variable oxidation products. The δ¹⁸O of these oxidation products at the surface of oxidized pyrites are strongly variable (from ≈−40‰ to ≈+30‰) for all experiments.Isotopic fractionation between sulphates and pyrite, Δ³⁴S_SO4-pyrite, is equal to −1.3‰ and +0.4‰ for sulphates formed by stoichiometric and non-stoichiometric processes, respectively. These two values likely reflect either a S-S or a Fe-S bond breaking process. The Δ¹⁸O_SO4-H2O and Δ¹⁸O_SO4-O2 are estimated to be ≈+16‰ and ≈−25‰, respectively. These values are higher than previously published data and may reflect biological effects. The large δ¹⁸O heterogeneity measured at the surfaces of oxidized pyrites, whatever the oxidant, may be related (i) to the existence of local surface environments isolated from the solution in which the oxidation processes are different and (ii) to the stabilization at the pyrite surface of reaction intermediates that are not in isotopic equilibrium with the solution. Though the oxygen isotopic composition of surface oxidation products cannot be taken as a direct biosignature, the combined morphological, chemical and isotopic characterization of the surfaces of oxidized pyrites may furnish clues about a biological activity on a mineral surface.     

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.     

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.     

Oxygen isotope partitioning during oxidation of pyrite by H2O2 and its dependence on temperature

A detailed experimental study was conducted to investigate mechanisms of pyrite oxidation by determining product yields and oxygen isotopic fractionation during reactions between powdered pyrite (FeS₂) with aqueous hydrogen peroxide (H₂O₂). Sealed silica-tube experiments utilized aliquots of pyrite that were reacted with 0.2 M H₂O₂ for 7 to 14 days at 4 to 150 °C. No volatile sulfur species were detected in any experiment. The only gaseous product recovered was elemental oxygen inferred to result from decomposition of H₂O₂. Aqueous sulfate (S_aq) was the only sulfur product recovered from solution. Solid hydrated ferric iron sulfates (i.e., water-soluble sulfate fraction, S_ws) were recovered from all experiments. Ferric oxide (hematite) was detected only in high temperature experiments.Reactants were selected with large differences in initial δ¹⁸O values. The oxygen isotopic compositions of oxygen-bearing reactants and products were analyzed for each experiment. Subsequent isotopic mass-balances were used to identify sources of oxygen for reaction products and to implicate specific chemical reaction mechanisms. δ¹⁸O of water did not show detectable change during any experiment. δ¹⁸O of sulfate was similar for S_aq and S_ws and indicated that both H₂O and H₂O₂ were sources of oxygen in sulfate. Low-temperature experiments suggest that H₂O-derived oxygen was incorporated into sulfate via Fe³⁺ oxidation, whereas H₂O₂-derived oxygen was incorporated into sulfate via oxidation by hydroxyl radicals (HO). These two competing mechanisms for oxygen incorporation into sulfate express comparable influences at 25 °C. With increasing reaction temperatures from 4 to 100 °C, it appears that accelerated thermal decomposition and diminished residence time of H₂O₂ limit the oxygen transfer from H₂O₂ into sulfate and enhance the relative importance of H₂O-derived oxygen for incorporation into sulfate. Notably, at temperatures between 100 and 150 °C there is a reversal in the lower temperature trend resulting in dominance of H₂O₂-derived oxygen over H₂O-derived oxygen. At such high temperatures, complete thermal decomposition of H₂O₂ to water and molecular oxygen (O₂) occurs within minutes in mineral-blank experiments and suggests little possibility for direct oxidation of pyrite by H₂O₂ above 100 °C. We hypothesize that a Fe-O₂ mechanism is responsible for oxygenating pyrite to sulfate using O₂ from the preceding thermal decomposition of H₂O₂.     

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

The effect of early diagenesis on the Fe isotope compositions of porewaters and authigenic minerals in continental margin sediments

Iron isotope compositions in marine pore fluids and sedimentary solid phases were measured at two sites along the California continental margin, where isotope compositions range from δ⁵⁶Fe = −3.0‰ to +0.4‰. At one site near Monterey Canyon off central California, organic matter oxidation likely proceeds through a number of diagenetic pathways that include significant dissimilatory iron reduction (DIR) and bacterial sulfate reduction, whereas at our other site in the Santa Barbara basin DIR appears to be comparatively small, and production of sulfides (FeS and pyrite) was extensive. The largest range in Fe isotope compositions is observed for Fe(II)_aq in porewaters, which generally have the lowest δ⁵⁶Fe values (minimum: −3.0‰) near the sediment surface, and increase with burial depth. δ⁵⁶Fe values for FeS inferred from HCl extractions vary between ∼−0.4‰ and +0.4‰, but pyrite is similar at both stations, where an average δ⁵⁶Fe value of −0.8 ± 0.2‰ was measured. We interpret variations in dissolved Fe isotope compositions to be best explained by open-system behavior that involves extensive recycling of Feflux. This study is the first to examine Fe isotope variations in modern marine sediments, and the results show that Fe isotopes in the various reactive Fe pools undergo isotopic fractionation during early diagenesis. Importantly, processes dominated by sulfide formation produce high-δ⁵⁶Fe values for porewaters, whereas the opposite occurs when Fe(III)-oxides are present and DIR is a major pathway of organic carbon respiration. Because shelf pore fluids may carry a negative δ⁵⁶Fe signature it is possible that the Fe isotope composition of ocean water reflects a significant contribution of shelf-derived iron to the open ocean. Such a signature would be an important means for tracing iron sources to the ocean and water mass circulation.     

The behavior of molybdenum and its isotopes across the chemocline and in the sediments of sulfidic Lake Cadagno, Switzerland

Molybdenum (Mo) isotope studies in black shales can provide information about the redox evolution of the Earth’s oceans, provided the isotopic consequences of Mo burial into its major sinks are well understood. Previous applications of the Mo isotope paleo-ocean redox proxy assumed quantitative scavenging of Mo when buried into sulfidic sediments. This paper contains the first complete suite of Mo isotope fractionation observations in a sulfidic water column and sediment system, the meromictic Lake Cadagno, Switzerland, a small alpine lake with a pronounced oxygen-sulfide transition reaching up to H₂S ∼ 200 μM in the bottom waters (or about 300 μM total sulfide: ΣS²⁻ = H₂S + HS⁻ + S²⁻). We find that Mo behaves conservatively in the oxic zone and non-conservatively in the sulfidic zone, where dissolved Mo concentrations decrease from 14 nM to 2-8 nM across this transition. Dissolved Mo in the upper oxic waters has a δ⁹⁸Mo_oxic = 0.9 ± 0.1‰, which matches that of the riverine input, δ⁹⁸Mo_river = 0.9 ± 0.1‰. In the deeper sulfidic waters, a subaquatic source delivers Mo at 1.55 ± 0.1‰, but the dissolved Mo is even heavier at δ⁹⁸Mo_sulfidic = 1.8‰. Sediment traps in the sulfidic zone of the lake collect particles increasingly enriched in Mo with depth, with δ⁹⁸Mo values significantly fractionated at −0.8‰ to −1.2‰ both near the chemocline and in the deepest trap. Suspended particulates in the sulfidic waters carry lighter Mo than the ambient dissolved Mo pool by ∼0.3-1.5‰. Sedimentary Mo concentrations correlate with total organic carbon and yield Mo levels which are two orders of magnitude higher than typical crustal values found in rocks from the catchment area. Solid-phase Mo in the sediment shows a slightly positive δ⁹⁸Mo trend with depth, from δ⁹⁸Mo = 1.2‰ to 1.4‰ while the pore waters show dramatic enrichments of Mo (>2000 nM) with a relatively light isotope signature of δ⁹⁸Mo = 0.9-1.0‰.These data are explained if Mo is converted to particle-reactive oxythiomolybdates in the sulfidic waters and is fractionated during removal from solution onto particles. Isotope fractionation is expressed in the water column, despite the high sulfide concentrations, because the rate of Mo removal is fast compared to the slow reaction kinetics of thiomolybdate formation. However, elemental and isotopic mass balances show that Mo is indeed quantitatively removed to the lake sediments and thus the isotopic composition of the sediments reflects sources to the sulfidic water. This efficient Mo drawdown is expected to occur in settings where H₂S is very much in excess over Mo or in a restricted setting where the water renewal rate is slow compared to the Mo burial rate. We present a model for the Mo isotope fractionation in sulfidic systems associated with the slow reaction kinetics and conclude that quantitative removal will occur in highly sulfidic and restricted marine systems.     

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

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.     

Triple-oxygen-isotope determination of molecular oxygen incorporation in sulfate produced during abiotic pyrite oxidation (pH = 2-11)

Aqueous oxidation of sulfide minerals to sulfate is an integral part of the global sulfur and oxygen cycles. The current model for pyrite oxidation emphasizes the role of Fe²⁺-Fe³⁺ electron shuttling and repeated nucleophilic attack by water molecules on sulfur. Previous δ¹⁸O-labeled experiments show that a variable fraction (0-60%) of the oxygen in product sulfate is derived from dissolved O₂, the other potential oxidant. This indicates that nucleophilic attack cannot continue all the way to sulfate and that a sulfoxyanion of intermediate oxidation state is released into solution. The observed variability in O₂% may be due to the presence of competing oxidation pathways, variable experimental conditions (e.g. abiotic, biotic, or changing pH value), or uncertainties related to the multiple experiments needed to effectively use the δ¹⁸O label to differentiate sulfate-oxygen sources. To examine the role of O₂ and Fe³⁺ in determining the final incorporation of O₂ oxygen in sulfate produced during pyrite oxidation, we designed a set of aerated, abiotic, pH-buffered (pH = 2, 7, 9, 10, and 11), and triple-oxygen-isotope labeled solutions with and without Fe³⁺ addition. While abiotic and pH-buffered conditions help to eliminate variables, triple oxygen isotope labeling and Fe³⁺ addition help to determine the oxygen sources in sulfate and examine the role of Fe²⁺-Fe³⁺ electron shuttling during sulfide oxidation, respectively.Our results show that sulfate concentration increased linearly with time and the maximum concentration was achieved at pH 11. At pH 2, 7, and 9, sulfate production was slow but increased by 4× with the addition of Fe³⁺. Significant amounts of sulfite and thiosulfate were detected in pH ? 9 reactors, while concentrations were low or undetectable at pH 2 and 7. The triple oxygen isotope data show that at pH ? 9, product sulfate contained 21-24% air O₂ signal, similar to pH 2 with Fe³⁺ addition. Sulfate from the pH 2 reactor without Fe³⁺ addition and the pH 7 reactors all showed 28-29% O₂ signal. While the O₂% in final sulfate apparently clusters around 25%, the measurable deviations (>experimental error) from the 25% in many reaction conditions suggest that (1) O₂ does get incorporated into intermediate sulfoxyanions (thiosulfate and sulfite) and a fraction survives sulfite-water exchange (e.g. the pH 2 with no Fe³⁺ addition and both pH 7 reactors); and (2) direct O₂ oxidation dominates while Fe³⁺ shuttling is still competitive in the sulfite-sulfate step (e.g. the pH 9, 10, and 11 and the pH 2 reactor with Fe³⁺ addition). Overall, the final sulfate-oxygen source ratio is determined by (1) rate competitions between direct O₂ incorporation and Fe³⁺ shuttling during both the formation of sulfite from pyrite and from sulfite to final sulfate, and (2) rate competitions between sulfite and water oxygen exchange and the oxidation of sulfite to sulfate. Our results indicate that thiosulfate or sulfite is the intermediate species released into solution at all investigated pH and point to a set of dynamic and competing fractionation factors and rates, which control the oxygen isotope composition of sulfate derived from pyrite oxidation.     

Stable isotope geochemistry of ground and surface waters associated with undisturbed massive sulfide deposits; constraints on origin of waters and water-rock reactions

Stable isotopes (H, O, C) were determined for ground and surface waters collected from two relatively undisturbed massive sulfide deposits (Halfmile Lake and Restigouche) in the Bathurst Mining Camp (BMC), New Brunswick, Canada. Additional waters from active and inactive mines in the BMC were also collected. Oxygen and hydrogen isotopes of surface and shallow groundwaters from both the Halfmile Lake and Restigouche deposits are remarkably uniform (− 13 to − 14‰ and − 85 to − 95‰ for δ¹⁸O_VSMOW and δ²H_VSMOW, respectively). These values are lighter than predicted for northern New Brunswick and, combined with elevated deuterium excess values, suggest that recharge waters are dominated by winter precipitation, recharged during spring melting. Deeper groundwaters from the Restigouche deposit, and from active and inactive mines have heavier δ¹⁸O_VSMOW ratios (up to − 10.8‰) than shallow groundwaters suggesting recharge under warmer climate or mixing with Shield-type brines. Some of the co-variation in Cl concentrations and δ¹⁸O_VSMOW ratios can be explained by mixing between saline and shallow recharge water end-members. Carbon isotopic compositions of dissolved inorganic carbon (DIC) are variable, ranging from − 15 to − 5‰ δ¹³C_VPDB for most ground and surface waters. Much of the variation in the carbon isotopes is consistent with closed system groundwater evolution involving soil zone CO₂ and fracture zone carbonate minerals (calcite, dolomite and siderite; average = − 6.5‰ δ¹³C_VPDB). The DIC of saline Restigouche deposit groundwater is isotopically heavy (∼+ 12‰ δ¹³C_VPDB), indicating carbon isotopic fractionation from methanogenesis via CO₂ reduction, consistent with the lack of dissolved sulfate in these waters and the observation of CH₄-degassing during sampling.     

The role of carbonate ions in pyrite oxidation in aqueous systems

The mechanism of pyrite oxidation in carbonate-containing alkaline solutions at 80 °C was investigated with the help of rate experiments, thermodynamic modeling and diffuse reflectance infrared spectroscopy (DRIFTS). Pyrite oxidation rate increased with pH and was enhanced by addition of bicarbonate/carbonate ions. The carbonate effect was found to be limited to moderately alkaline conditions (pH 8-11). Metastable Eh-pH diagrams, at 25 °C, indicate that soluble iron-carbonate complexes (FeHCO₃⁻, FeCO₃⁰, Fe(CO₃)(OH)⁻ and FeCO₃²⁻) may coexist with pyrite in the pH range of 6-12.5. Above pH 11 and 13, the Fe(II) and Fe(III) hydroxocomplexes, respectively, become stable, even in the presence of carbonate/bicarbonate ions. Surface-bound carbonate complexes on iron were also identified with DRIFTS as products of pyrite oxidation in addition to iron oxyhydroxides and soluble sulfate species. The conditions under which thermodynamic and DRIFTS analyses indicate the presence of carbonate compounds also correspond to those in which the fastest rate of pyrite oxidation in carbonate solutions was observed. Following the Singer-Stumm model for pyrite oxidation in acidic solutions, it is assumed that Fe(III) is the preferred pyrite oxidant under alkaline conditions. We propose that carbonate ions facilitate the electron transfer from soluble iron(II)-carbonate to O₂, increase the iron solubility, and provide buffered, favorable alkaline conditions at the reaction front, which in turn favors the overall kinetics of pyrite oxidation. Therefore, the electron transfer from sulfur atoms to O₂ is facilitated by the formation of the cycle of Fe(II)-pyrite/Fe(III)-carbonate redox couple at the pyrite surface.     

Different types of methane monooxygenases produce similar carbon and hydrogen isotope fractionation patterns during methane oxidation

We determined the stable carbon and hydrogen isotope fractionation factors for methane oxidation under oxic conditions using strains with known degradation pathways. The aerobic oxidation of methane can be initiated by two different forms of enzymes known as methane monooxygenases (MMO). The expression of these enzymes is type-specific and dependent upon the adjusted copper concentration in the medium (or environment). In this study, the expression of either the soluble MMO or the particulate MMO was supported by adjusting the copper concentrations in the growth medium. Taxonomically different aerobic methanotrophic strains, mainly belonging to the alpha- and gamma- classes of Proteobacteria, produced methane isotope enrichment factors (ε_bulk) ranging from −14.8 to −27.9‰ for carbon, and from −110.0 to −231.5‰ for hydrogen. The ratio of hydrogen versus carbon discrimination (Λ = (α_H⁻¹ − 1)/(α_C⁻¹ − 1) ≈ Δ(δ²H)/Δ(δ¹³C)) were similar for all tested cultures, and are also identical to values calculated from previously published enrichment factors for aerobic and anaerobic methane degradation. In contrast, Λ-values for the abiotic oxidation of methane with OH radicals (this process is considered as the main removal process for methane from the atmosphere) were significantly higher than the values derived from biotic oxidation. Due to the low variability of microbial methane isotope fractionation patterns, we propose that combined carbon and hydrogen isotope fractionation analyses can be used to monitor and assess the occurrence of microbial methane oxidation in marine or terrestrial environments. However, it is not possible to distinguish distinct aerobic or anaerobic methane-oxidation pathways by this approach.     

Effects of changing solution chemistry on Fe/Fe isotope fractionation in aqueous Fe-Cl solutions

The range in ⁵⁶Fe/⁵⁴Fe isotopic compositions measured in naturally occurring iron-bearing species is greater than 5‰. Both theoretical modeling and experimental studies of equilibrium isotopic fractionation among iron-bearing species have shown that significant fractionations can be caused by differences in oxidation state (i.e., redox effects in the environment) as well as by bond partner and coordination number (i.e., nonredox effects due to speciation).To test the relative effects of redox vs. nonredox attributes on total Fe equilibrium isotopic fractionation, we measured changes, both experimentally and theoretically, in the isotopic composition of an Fe²⁺-Fe³⁺-Cl-H₂O solution as the chlorinity was varied. We made use of the unique solubility of FeCl₄⁻ in immiscible diethyl ether to create a separate spectator phase against which changes in the aqueous phase could be quantified. Our experiments showed a reduction in the redox isotopic fractionation between Fe²⁺- and Fe³⁺-bearing species from 3.4‰ at [Cl⁻] = 1.5 M to 2.4‰ at [Cl⁻] = 5.0 M, due to changes in speciation in the Fe-Cl solution. This experimental design was also used to demonstrate the attainment of isotopic equilibrium between the two phases, using a ⁵⁴Fe spike.To better understand speciation effects on redox fractionation, we created four new sets of ab initio models of the ferrous chloride complexes used in the experiments. These were combined with corresponding ab initio models for the ferric chloride complexes from previous work. At 20 °C, 1000 ln β (β = ⁵⁶Fe/⁵⁴Fe reduced partition function ratio relative to a dissociated Fe atom) values range from 6.39‰ to 5.42‰ for Fe(H₂O)₆²⁺, 5.98‰ to 5.34‰ for FeCl(H₂O)₅⁺, and 5.91‰ to 4.86‰ for FeCl₂(H₂O)₄, depending on the model. The theoretical models predict ferric-ferrous fractionation about half as large (depending on model) as the experimental results.Our results show (1) oxidation state is likely to be the dominant factor controlling equilibrium Fe isotope fractionation in solution and (2) nonredox attributes (such as ligands present in the aqueous solution, speciation and relative abundances, and ionic strength of the solution) can also have significant effects. Changes in the isotopic composition of an Fe-bearing solution will influence the resultant Fe isotopic signature of any precipitates.