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.     

Mineralogic and sulfur isotopic effects accompanying oxidation of pyrite in millimolar solutions of hydrogen peroxide at temperatures from 4 to 150 °C

Oxidation of pyrite by hydrogen peroxide (H₂O₂) at millimolar levels has been studied from 4 to 150 °C in order to evaluate isotopic effects potentially associated with radiolytic oxidation of pyrite. Gaseous, aqueous, and solid phases were collected and measured following sealed-tube experiments that lasted from 1 to 14 days. The dominant gaseous product was molecular oxygen. No volatile sulfur species were recovered from any experiment. Sulfate was the only aqueous sulfur species detected in solution, with sulfite and thiosulfate below the detection limits. X-ray diffraction patterns and images from scanning electron microscopy reveal solid residues composed primarily of hydrated ferric iron sulfates and sporadic ferric-ferrous iron sulfates. Hematite was detected only in solid residue produced during high temperature experiments. Elemental sulfur and/or polysulfides are inferred to be form on reacting pyrite surface based on extraction with organic solvents. Pyrite oxidation by H₂O₂ increases in rate with increasing H₂O₂concentration, pyrite surface area, and temperature. Rates measured in sealed-tube experiments at 25°C, for H₂O₂ concentration of 2 × 10⁻³ M are 8.8 × 10⁻⁹ M/m²/sec, which are higher than previous estimates. A combination of reactive oxygen species from H₂O₂ decomposition products and reactive iron species from pyrite dissolution is inferred to aggressively oxidize the receding pyrite surface. Competing oxidants with temperature-dependent oxidation efficiencies results in multiple reaction mechanisms for different temperatures and surface conditions. Sulfur isotope values of remaining pyrite were unchanged during the experiments, but showed distinct enrichment of ³⁴S in produced sulfate and depletion in elemental sulfur. The Δ_{sulfate-pyrite} and Δ_{elemental sulfur-pyrite} was +0.5 to +1.5‰ and was −0.2 to −1‰, respectively. Isotope data from high-temperature experiments indicate an additional ³⁴S-depleted sulfur fraction, with up to 4‰ depletion of ³⁴S, in the hematite. Sulfur isotope trends were not influenced by H₂O₂ concentration, temperature, or reaction time. Results of this study indicate that radiolytically produced oxidants, such as hydrogen peroxide and hydroxyl radicals, could efficiently oxidize pyrite in an otherwise oxygen-limited environment. Although H₂O₂ is generally regarded as being of minor geochemical significance on Earth, the H₂O₂ molecule plays a pivotal role in Martian atmospheric and soil chemistry. Additional experimental and field studies are needed to characterize sulfur and oxygen isotope systematics during radiolytical oxidation of metallic sulfides and elemental sulfur.     

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.     

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.     

Molecular H2O in armenite, BaCa2Al6Si9O30·2H2O, and epididymite, Na2Be2Si6O15·H2O: Heat capacity, entropy and local-bonding behavior of confined H2O in microporous silicates

Armenite, ideal formula BaCa₂Al₆Si₉O₃₀·2H₂O, and its dehydrated analog BaCa₂Al₆Si₉O₃₀ and epididymite, ideal formula Na₂Be₂Si₆O₁₅·H₂O, and its dehydrated analog Na₂Be₂Si₆O₁₅ were studied by low-temperature relaxation calorimetry between 5 and 300 K to determine the heat capacity, C_p, behavior of their confined H₂O. Differential thermal analysis and thermogravimetry measurements, FTIR spectroscopy, electron microprobe analysis and powder Rietveld refinements were undertaken to characterize the phases and the local environment around the H₂O molecule.The determined structural formula for armenite is Ba_0.88(0.01)Ca_1.99(0.02)Na_0.04(0.01)Al_5.89(0.03)Si_9.12(0.02)O₃₀·2H₂O and for epididymite Na_1.88(0.03)K_0.05(0.004)Na_0.01(0.004)Be_2.02(0.008)Si_6.00(0.01)O₁₅·H₂O. The infrared (IR) spectra give information on the nature of the H₂O molecules in the natural phases via their H₂O stretching and bending vibrations, which in the case of epididymite only could be assigned. The powder X-ray diffraction data show that armenite and its dehydrated analog have similar structures, whereas in the case of epididymite there are structural differences between the natural and dehydrated phases. This is also reflected in the lattice IR mode behavior, as observed for the natural phases and the H₂O-free phases. The standard entropy at 298 K for armenite is S° = 795.7 ± 6.2 J/mol K and its dehydrated analog is S° = 737.0 ± 6.2 J/mol K. For epididymite S° = 425.7 ± 4.1 J/mol K was obtained and its dehydrated analog has S° = 372.5 ± 5.0 J/mol K. The heat capacity and entropy of dehydration at 298 K are Δ = 3.4 J/mol K and ΔS^rxn = 319.1 J/mol K and Δ = −14.3 J/mol K and ΔS^rxn = 135.7 J/mol K for armenite and epididymite, respectively. The H₂O molecules in both phases appear to be ordered. They are held in place via an ion-dipole interaction between the H₂O molecule and a Ca cation in the case of armenite and a Na cation in epididymite and through hydrogen-bonding between the H₂O molecule and oxygen atoms of the respective silicate frameworks. Of the three different H₂O phases ice, liquid water and steam, the C_p behavior of confined H₂O in both armenite and epididymite is most similar to that of ice, but there are differences between the two silicates and from the C_p behavior of ice. Hydrogen-bonding behavior and its relation to the entropy of confined H₂O at 298 K is analyzed for various microporous silicates.The entropy of confined H₂O at 298 K in various silicates increases approximately linearly with increasing average wavenumber of the OH-stretching vibrations. The interpretation is that decreased hydrogen-bonding strength between a H₂O molecule and the silicate framework, as well as weak ion-dipole interactions, results in increased entropy of H₂O. This results in increased amplitudes of external H₂O vibrations, especially translations of the molecule, and they contribute strongly to the entropy of confined H₂O at T < 298 K.     

Kinetics of H2-O2-H2O redox equilibria and formation of metastable H2O2 under low temperature hydrothermal conditions

Hydrothermal experiments were conducted to evaluate the kinetics of H_2(aq) oxidation in the homogeneous H₂-O₂-H₂O system at conditions reflecting subsurface/near-seafloor hydrothermal environments (55-250 °C and 242-497 bar). The kinetics of the water-forming reaction that controls the fundamental equilibrium between dissolved H_2(aq) and O_2(aq), are expected to impose significant constraints on the redox gradients that develop when mixing occurs between oxygenated seawater and high-temperature anoxic vent fluid at near-seafloor conditions. Experimental data indicate that, indeed, the kinetics of H_2(aq)-O_2(aq) equilibrium become slower with decreasing temperature, allowing excess H_2(aq) to remain in solution. Sluggish reaction rates of H_2(aq) oxidation suggest that active microbial populations in near-seafloor and subsurface environments could potentially utilize both H_2(aq) and O_2(aq), even at temperatures lower than 40 °C due to H_2(aq) persistence in the seawater/vent fluid mixtures. For these H₂-O₂ disequilibrium conditions, redox gradients along the seawater/hydrothermal fluid mixing interface are not sharp and microbially-mediated H_2(aq) oxidation coupled with a lack of other electron acceptors (e.g. nitrate) could provide an important energy source available at low-temperature diffuse flow vent sites.More importantly, when H_2(aq)-O_2(aq) disequilibrium conditions apply, formation of metastable hydrogen peroxide is observed. The yield of H₂O_2(aq) synthesis appears to be enhanced under conditions of elevated H_2(aq)/O_2(aq) molar ratios that correspond to abundant H_2(aq) concentrations. Formation of metastable H₂O₂ is expected to affect the distribution of dissolved organic carbon (DOC) owing to the existence of an additional strong oxidizing agent. Oxidation of magnetite and/or Fe⁺⁺ by hydrogen peroxide could also induce formation of metastable hydroxyl radicals (•OH) through Fenton-type reactions, further broadening the implications of hydrogen peroxide in hydrothermal environments.     

Pyrite (FeS2) oxidation: A sub-micron synchrotron investigation of the initial steps

Pyrite is an environmentally significant mineral being the major contributor to acid rock drainage. Synchrotron based SPEM (scanning photoelectron microscopy) and micro-XPS (X-ray photoelectron spectroscopy) have been used to characterise fresh and oxidised pyrite (FeS₂) with a view to understanding the initial oxidation steps that take place during natural weathering processes. Localised regions of the pyrite surface containing Fe species of reduced coordination have been found to play a critical role. Such sites not only initiate the oxidation process but also facilitate the formation of highly reactive hydroxyl radical species, which then lead the S oxidation process.Four different S species are found to be present on fresh fractured pyrite surfaces: S₂²⁻_(bulk) (4-fold coordination), S₂²⁻_(surface) (3-fold coordination), S²⁻ and S⁰/S_n²⁻ (metal deficient sulfide and polysulfide respectively). These species were found to be heterogeneously distributed on the fractured pyrite surface. Both O₂ and H₂O gases are needed for effective oxidation of the pyrite surface. The process is initiated when O₂ dissociatively and H₂O molecularly adsorb onto the surface Fe sites where high dangling bond densities exist. H₂O may then dissociate to produce OH radicals. The adsorption of these species leads to the formation of Fe-oxy species prior to the formation of sulfoxy species. Evidence suggests that Fe-O bonds form prior to Fe-OH bonds. S oxidation occurs through interactions of OH radicals formed at the Fe sites, with formation of SO₄²⁻ occurring via S₂O₃²⁻/SO₃²⁻ intermediates. The pyrite oxidation process is electrochemical in nature and was found to occur in patches, where site specific adsorption of O₂ and H₂O has occurred. Fe and S oxidation was found to occur within the same area of oxidation probably in atomic scale proximity. Furthermore, the O in SO₄²⁻ arises largely from H₂O; however, depending on the surface history, SO₄²⁻ formed early in the oxidation process may also contain O from O₂.     

Kinetic nitrogen isotope fractionation associated with thermal decomposition of NH3: Experimental results and potential applications to trace the origin of N2 in natural gas and hydrothermal systems

Ammonia (NH₃) is the major intermediate phase in the pathway of nitrogen (N) transfer from the fixed N phases (e.g., in crustal material) to free N₂ (e.g., in natural gas reservoirs and volcanic gases). Yet the N isotopic behavior during these N-cycling processes remains poorly known. In an attempt to contribute to the understanding of N cycling using N isotopes, we carried out laboratory experiments to investigate the N isotopic effect associated with thermal decomposition of ammonia (2NH₃ → N₂ + 3H₂). Pure NH₃ (with initial δ¹⁵N_NH3 of ∼ −2‰, relative to air standard) was sealed into quartz tubes and thermally decomposed at 600, 700 or 800 °C from 2 hours to 500 days. With the progress of the reaction, the δ¹⁵N of the remaining NH₃ and the accumulated N₂ increased from −2 to +35‰ and from −20 to −2‰, respectively. The differences of the N-isotope fractionations at the three temperatures are not significant. Modeling using the Rayleigh distillation model yielded similar kinetic N-isotope fractionation factors (α_N2-NH3) of 0.983 ± 0.002 for 600, 700 and 800 °C. Applied to geological settings, this significant isotope discrimination (∼17‰) associated with partial decomposition of NH₃/NH₄⁺ from crustal sources (δ¹⁵N_average ∼ +6.3‰) can produce mantle-like (i.e. ∼ −5‰) or even lower δ¹⁵N values of N₂. This may explain the large variation of δ¹⁵N (−20 to +30‰) of N₂ in natural gas reservoirs. It can also possibly explain the extreme ¹⁵N-depletion of N₂ in some volcanic gases. This possibility has to be carefully considered when using N isotopes to trace geological N cycling across subduction zones by analysis of volcanic N₂.     

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.     

Interaction between sulphide and H2O in silicate melts

Reaction between dissolved water and sulphide was experimentally investigated in soda-lime-silicate (NCS) and sodium trisilicate (NS3) melts at temperatures from 1000 to 1200 °C and pressures of 100 or 200 MPa in internally heated gas pressure vessels. Diffusion couple experiments were conducted at water-undersaturated conditions with one half of the couple being doped with sulphide (added as FeS or Na₂S; 1500-2000 ppm S by weight) and the other with H₂O (∼3.0 wt.%). Additionally, two experiments were performed using a dry NCS glass cylinder and a free H₂O fluid. Here, the melt was water-saturated at least at the melt/fluid interface. Profiling by electron microprobe (sulphur) and infrared microscopy (H₂O) demonstrate that H₂O diffusion in the melts is faster by 1.5-2.3 orders of magnitude than sulphur diffusion and, hence, H₂O can be considered as a rapidly diffusing oxidant while sulphur is quasi immobile in these experiments.In Raman spectra a band at 2576 cm⁻¹ appears in the sulphide - H₂O transition zone which is attributed to fundamental S-H stretching vibrations. Formation of new IR absorption bands at 5025 cm⁻¹ (on expense of the combination band of molecular H₂O at 5225 cm⁻¹) and at 3400 cm⁻¹ was observed at the front of the in-diffusing water in the sulphide bearing melt. The appearance and intensity of these two IR bands is correlated with systematic changes in S K-edge XANES spectra. A pre-edge excitation at 2466.5 eV grows with increasing H₂O concentration while the sulphide peak at 2474.0 eV decreases in intensity relative to the peak at 2477.0 eV and the feature at 2472.3 eV becomes more pronounced (all energies are relative to the sulphate excitation, calibrated to 2482.5 eV). The observations by Raman, IR and XANES spectroscopy indicate a well coordinated S²⁻ - H₂O complex which was probably formed in the glasses during cooling at the glass transition. No oxidation of sulphide was observed in any of the diffusion couple experiments. On the contrary, XANES spectra from experiments conducted with a free H₂O fluid show complete transformation of sulphide to sulphate near the melt surface and coexistence of sulphate and sulphide in the center of the melt. This can be explained by a lower H₂O activity in the diffusion couple experiments or by the need of a sink for hydrogen (e.g., a fluid which can dissolve high concentration of hydrogen) to promote oxidation of sulphide by H₂O via the reaction S²⁻ + 4H₂O = SO₄²⁻ + 4H₂. Sulphite could not be detected in any of the XANES spectra implying that this species, if it exists in the melt, it is a subordinate or transient species only.     

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.     

Iron isotope fractionation between pyrite (FeS2), hematite (Fe2O3) and siderite (FeCO3): A first-principles density functional theory study

In addition to equilibrium isotopic fractionation factors experimentally derived, theoretical predictions are needed for interpreting isotopic compositions measured on natural samples because they allow exploring more easily a broader range of temperature and composition. For iron isotopes, only aqueous species were studied by first-principles methods and the combination of these data with those obtained by different methods for minerals leads to discrepancies between theoretical and experimental isotopic fractionation factors. In this paper, equilibrium iron isotope fractionation factors for the common minerals pyrite, hematite, and siderite were determined as a function of temperature, using first-principles methods based on the density functional theory (DFT). In these minerals belonging to the sulfide, oxide and carbonate class, iron is present under two different oxidation states and is involved in contrasted types of interatomic bonds. Equilibrium fractionation factors calculated between hematite and siderite compare well with the one estimated from experimental data (ln α⁵⁷Fe/⁵⁴Fe = 4.59 ± 0.30‰ and 5.46 ± 0.63‰ at 20 °C for theoretical and experimental data, respectively) while those for Fe(III)_aq-hematite and Fe(II)_aq-siderite are significantly higher that experimental values. This suggests that the absolute values of the reduced partition functions (β-factors) of aqueous species are not accurate enough to be combined with those calculated for minerals. When compared to previous predictions derived from Mössbauer or INRXS data [Polyakov V. B., Clayton R. N., Horita J. and Mineev S. D. (2007) Equilibrium iron isotope fractionation factors of minerals: reevaluation from the data of nuclear inelastic resonant X-ray scattering and Mössbauer spectroscopy. Geochim. Cosmochim. Acta71, 3833-3846], our iron β-factors are in good agreement for siderite and hematite while a discrepancy is observed for pyrite. However, the detailed investigation of the structural, electronic and vibrational properties of pyrite as well as the study of sulfur isotope fractionation between pyrite and two other sulfides (sphalerite and galena) indicate that DFT-derived β-factors of pyrite are as accurate as for hematite and siderite. We thus suggest that experimental vibrational density of states of pyrite should be re-examined.     

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

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

Stable sulfur isotope partitioning during simulated petroleum formation as determined by hydrous pyrolysis of Ghareb Limestone, Israel

Hydrous pyrolysis experiments at 200 to 365°C were carried out on a thermally immature organic-rich limestone containing Type-IIS kerogen from the Ghareb Limestone in North Negev, Israel. This work focuses on the thermal behavior of both organic and inorganic sulfur species and the partitioning of their stable sulfur isotopes among organic and inorganic phases generated during hydrous pyrolyses. Most of the sulfur in the rock (85%) is organic sulfur. The most dominant sulfur transformation is cleavage of organic-bound sulfur to form H₂S_(gas). Up to 70% of this organic sulfur is released as H₂S_(gas) that is isotopically lighter than the sulfur in the kerogen. Organic sulfur is enriched by up to 2‰ in ³⁴S during thermal maturation compared with the initial δ³⁴S values. The δ³⁴S values of the three main organic fractions (kerogen, bitumen and expelled oil) are within 1‰ of one another. No thermochemical sulfate reduction or sulfate formation was observed during the experiments. The early released sulfur reacted with available iron to form secondary pyrite and is the most ³⁴S depleted phase, which is 21‰ lighter than the bulk organic sulfur. The large isotopic fractionation for the early formed H₂S is a result of the system not being in equilibrium. As partial pressure of H₂S_(gas) increases, retro reactions with the organic sulfur in the closed system may cause isotope exchange and isotopic homogenization. Part of the δ³⁴S-enriched secondary pyrite decomposes above 300°C resulting in a corresponding decrease in the δ³⁴S of the remaining pyrite. These results are relevant to interpreting thermal maturation processes and their effect on kerogen-oil-H₂S-pyrite correlations. In particular, the use of pyrite-kerogen δ³⁴S relations in reconstructing diagenetic conditions of thermally mature rocks is questionable because formation of secondary pyrite during thermal maturation can mask the isotopic signature and quantity of the original diagenetic pyrite. The main transformations of kerogen to bitumen and bitumen to oil can be recorded by using both sulfur content and δ³⁴S of each phase including the H₂S_(gas). H₂S generated in association with oil should be isotopically lighter or similar to oil. It is concluded that small isotopic differentiation obtained between organic and inorganic sulfur species suggests closed-system conditions. Conversely, open-system conditions may cause significant isotopic discrimination between the oil and its source kerogen. The magnitude of this discrimination is suggested to be highly dependent on the availability of iron in a source rock resulting in secondary formation of pyrite.     

The role of sulfur in chemical weathering and atmospheric CO2 fluxes: Evidence from major ions, δCDIC, and δSSO4 in rivers of the Canadian Cordillera

Water samples from the Fraser, Skeena and Nass River basins of the Canadian Cordillera were analyzed for dissolved major element concentrations (HCO₃⁻, SO₄²⁻, Cl⁻, Ca²⁺, Mg²⁺, K⁺, Na⁺), δ¹³C of dissolved inorganic carbon (δ¹³C_DIC), and δ³⁴S of dissolved sulfate (δ³⁴S_SO4) to quantify chemical weathering rates and exchanges of CO₂ between the atmosphere, hydrosphere, and lithosphere. Weathering rates of silicates and carbonates were determined from major element mass balance. Combining the major element mass balance with δ³⁴S_SO4 (−8.9 to 14.1‰_CDT) indicates sulfide oxidation (sulfuric acid production) and subsequent weathering of carbonate and to a lesser degree silicate minerals are important processes in the study area. We determine that on average, 81% of the riverine sulfate can be attributed to sulfide oxidation in the Cordilleran rivers, and that 25% of the total weathering cation flux can be attributed to carbonate and silicate dissolution by sulfuric acid. This result is validated by δ¹³C_DIC values (−9.8 to −3.7‰ _VPDB) which represents a mixture of DIC produced by the following weathering pathways: (i) carbonate dissolution by carbonic acid (−8.25‰) > (ii) silicate dissolution by carbonic acid (−17‰) ≈ (iii) carbonate dissolution by sulfuric acid derived from the oxidation of sulfides (coupled sulfide-carbonate weathering) (+0.5‰).δ³⁴S_SO4 is negatively correlated with δ¹³C_DIC in the Cordilleran rivers, which further supports the hypothesis that sulfuric acid produced by sulfide oxidation is primarily neutralized by carbonates, and that sulfide-carbonate weathering impacts the δ¹³C_DIC of rivers. The negative correlation between δ³⁴S_SO4 and δ¹³C_DIC is not observed in the Ottawa and St. Lawrence River basins. This suggests other factors such as landscape age (governed by tectonic uplift) and bedrock geology are important controls on regional sulfide oxidation rates, and therefore also on the magnitude of sulfide-carbonate weathering—i.e., it is more significant in tectonically active areas.Calculated DIC fluxes due to Ca and Mg silicate weathering by carbonic acid (38.3 × 10³ mol C · km⁻² · yr⁻¹) are similar in magnitude to DIC fluxes due to sulfide-carbonate weathering (18.5 × 10³ mol C · km⁻² · yr⁻¹). While Ca and Mg silicate weathering facilitates a transfer of atmospheric CO₂ to carbonate rocks, sulfide-carbonate weathering can liberate CO₂ from carbonate rocks to the atmosphere when sulfide oxidation exceeds sulfide deposition. This implies that in the Canadian Cordillera, sulfide-carbonate weathering can offset up to 48% of the current CO₂ drawdown by silicate weathering in the region.     

Structure of H2O-saturated peralkaline aluminosilicate melt and coexisting aluminosilicate-saturated aqueous fluid determined in-situ to 800 °C and ∼800 MPa

The structure of H₂O-saturated silicate melts and of silicate-saturated aqueous solutions, as well as that of supercritical silicate-rich aqueous liquids, has been characterized in-situ while the sample was at high temperature (to 800 °C) and pressure (up to 796 MPa). Structural information was obtained with confocal microRaman and with FTIR spectroscopy. Two Al-bearing glasses compositionally along the join Na₂O•4SiO₂-Na₂O•4(NaAl)O₂-H₂O (5 and 10 mol% Al₂O₃, denoted NA5 and NA10) were used as starting materials. Fluids and melts were examined along pressure-temperature trajectories of isochores of H₂O at nominal densities (from PVT properties of pure H₂O) of 0.85 g/cm³ (NA10 experiments) and 0.86 g/cm³ (NA5 experiments) with the aluminosilicate + H₂O sample contained in an externally-heated, Ir-gasketed hydrothermal diamond anvil cell.Molecular H₂O (H₂O°) and OH groups that form bonds with cations exist in all three phases. The OH/H₂O° ratio is positively correlated with temperature and pressure (and, therefore, fugacity of H₂O, f_H2O) with (OH/H₂O°)^melt > (OH/H₂O°)^fluid at all pressures and temperatures. Structural units of Q³, Q², Q¹, and Q⁰ type occur together in fluids, in melts, and, when outside the two-phase melt + fluid boundary, in single-phase liquids. The abundance of Q⁰ and Q¹ increases and Q² and Q³ decrease with f_H2O. Therefore, the NBO/T (nonbridging oxygen per tetrahedrally coordination cations), of melt is a positive function of f_H2O. The NBO/T of silicate in coexisting aqueous fluid, although greater than in melt, is less sensitive to f_H2O.The melt structural data are used to describe relationships between activity of H₂O and melting phase relations of silicate systems at high pressure and temperature. The data were also combined with available partial molar configurational heat capacity of Qⁿ-species in melts to illustrate how these quantities can be employed to estimate relationships between heat capacity of melts and their H₂O content.     

Grain-scale iron isotopic distribution of pyrite from Precambrian shallow marine carbonate revealed by a femtosecond laser ablation multicollector ICP-MS technique: Possible proxy for the redox state of ancient seawater

The redox state of Precambrian shallow seas has been linked with material cycle and evolution of the photosynthesis-based ecosystem. Iron is a redox-sensitive element and exists as a soluble Fe(II) species or insoluble Fe(III) species on Earth’s surface. Previous studies have shown that the iron isotopic ratio of marine sedimentary minerals is useful for understanding the ocean redox state, although the redox state of the Archean shallow sea is poorly known. This is partly because the conventional bulk isotope analytical technique has often been used, wherein the iron isotopic record may be dampened by the presence of isotopically different iron-bearing minerals within the same sample. Here we report a microscale iron isotopic ratio of individual pyrite grains in shallow marine stromatolitic carbonates over geological time using a newly developed, near-infrared femtosecond laser ablation multicollector ICP-MS technique (NIR-fs-LA-MC-ICP-MS).We have determined that the grain-scale iron isotopic distribution of pyrite from coeval samples shows a bimodal (2.7 and 2.3 Ga) or unimodal pattern (2.9, 2.6, and 0.7 Ga). In particular, pyrite from the 2.7 Ga Fortescue Group shows a unique bimodal distribution with highly positive (+1.0‰ defined as Type 1) and negative δ⁵⁶Fe values (−1.8‰ defined as Type 2). Type 1 and 2 pyrites occasionally occur within different siliceous layers in the same rock specimen. Layer-scale iron isotopic heterogeneity indicates that the iron isotopic ratios of the two types of pyrite are not homogenized by diagenesis after deposition. Some cubic pyrites have a core with a positive δ⁵⁶Fe value (1‰) and a rim with a crustal δ⁵⁶Fe value (0‰). The observed isotopic zoning suggests that the positive δ⁵⁶Fe value is a primary signature at the time of stromatolite formation, while secondary pyrite precipitated during diagenesis.The positive δ⁵⁶Fe value of Type 1 and the large iron isotopic difference between Type 1 and 2 (2.8‰.) suggest partial Fe(II) oxidation in the 2.7-Ga shallow sea, i.e., pyritization of ⁵⁶Fe-enriched ferric oxyhydroxide (Type 1) and ⁵⁶Fe depleted Fe²⁺aq in seawater (Type 2). Type 2 pyrite was probably not produced by microbial iron redox cycling during diagenesis because this scenario requires a higher abundance of pyrite with δ⁵⁶Fe of 0‰ than of −1.8‰. Consequently, the degree of Fe(II) oxidation in the 2.7-Ga shallow sea can be estimated by a Fe²⁺aq steady-state model. The model calculation shows that half the Fe²⁺aq influx was oxidized in the seawater. This implies that O₂ produced by photosynthesis would have been completely consumed by oxidation of the Fe²⁺aq influx. Grain-scale iron isotopic distribution of pyrite could be a useful index for reconstructing the redox state of the Archean shallow sea.     

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.     

Oxygen isotope exchange between H2O and CO2 at elevated CO2 pressures: Implications for monitoring of geological CO2 storage

Traditionally, the application of stable isotopes in Carbon Capture and Storage (CCS) projects has focused on δ¹³C values of CO₂ to trace the migration of injected CO₂ in the subsurface. More recently the use of δ¹⁸O values of both CO₂ and reservoir fluids has been proposed as a method for quantifying in situ CO₂ reservoir saturations due to O isotope exchange between CO₂ and H₂O and subsequent changes in δ¹⁸O_H2O values in the presence of high concentrations of CO₂. To verify that O isotope exchange between CO₂ and H₂O reaches equilibrium within days, and that δ¹⁸O_H2O values indeed change predictably due to the presence of CO₂, a laboratory study was conducted during which the isotope composition of H₂O, CO₂, and dissolved inorganic C (DIC) was determined at representative reservoir conditions (50 °C and up to 19 MPa) and varying CO₂ pressures. Conditions typical for the Pembina Cardium CO₂ Monitoring Pilot in Alberta (Canada) were chosen for the experiments. Results obtained showed that δ¹⁸O values of CO₂ were on average 36.4 ± 2.2‰ (1σ, n = 15) higher than those of water at all pressures up to and including reservoir pressure (19 MPa), in excellent agreement with the theoretically predicted isotope enrichment factor of 35.5‰ for the experimental temperatures of 50 °C. By using ¹⁸O enriched water for the experiments it was demonstrated that changes in the δ¹⁸O values of water were predictably related to the fraction of O in the system sourced from CO₂ in excellent agreement with theoretical predictions. Since the fraction of O sourced from CO₂ is related to the total volumetric saturation of CO₂ and water as a fraction of the total volume of the system, it is concluded that changes in δ¹⁸O values of reservoir fluids can be used to calculate reservoir saturations of CO₂ in CCS settings given that the δ¹⁸O values of CO₂ and water are sufficiently distinct.     

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.