- Research article
- Open Access
Multiple sulfur isotopes fractionations associated with abiotic sulfur transformations in Yellowstone National Park geothermal springs
© Kamyshny et al.; licensee BioMed Central Ltd. 2014
- Received: 1 September 2013
- Accepted: 12 March 2014
- Published: 28 May 2014
The paper presents a quantification of main (hydrogen sulfide and sulfate), as well as of intermediate sulfur species (zero-valent sulfur (ZVS), thiosulfate, sulfite, thiocyanate) in the Yellowstone National Park (YNP) hydrothermal springs and pools. We combined these measurements with the measurements of quadruple sulfur isotope composition of sulfate, hydrogen sulfide and zero-valent sulfur. The main goal of this research is to understand multiple sulfur isotope fractionation in the system, which is dominated by complex, mostly abiotic, sulfur cycling.
Water samples from six springs and pools in the Yellowstone National Park were characterized by pH, chloride to sulfate ratios, sulfide and intermediate sulfur species concentrations. Concentrations of sulfate in pools indicate either oxidation of sulfide by mixing of deep parent water with shallow oxic water, or surface oxidation of sulfide with atmospheric oxygen. Thiosulfate concentrations are low (<6 μmol L-1) in the pools with low pH due to fast disproportionation of thiosulfate. In the pools with higher pH, the concentration of thiosulfate varies, depending on different geochemical pathways of thiosulfate formation. The δ34S values of sulfate in four systems were close to those calculated using a mixing line of the model based on dilution and boiling of a deep hot parent water body. In two pools δ34S values of sulfate varied significantly from the values calculated from this model. Sulfur isotope fractionation between ZVS and hydrogen sulfide was close to zero at pH < 4. At higher pH zero-valent sulfur is slightly heavier than hydrogen sulfide due to equilibration in the rhombic sulfur–polysulfide – hydrogen sulfide system. Triple sulfur isotope (32S, 33S, 34S) fractionation patterns in waters of hydrothermal pools are more consistent with redox processes involving intermediate sulfur species than with bacterial sulfate reduction. Small but resolved differences in ∆33S among species and between pools are observed.
The variation of sulfate isotopic composition, the origin of differences in isotopic composition of sulfide and zero–valent sulfur, as well as differences in ∆33S of sulfide and sulfate are likely due to a complex network of abiotic redox reactions, including disproportionation pathways.
- Isotope geochemistry
- Sulfur cycle
- Sulfide oxidation
- Sulfur disproportionation
- Multiple sulfur isotopes
- Yellowstone national park
- Hydrothermal springs
Yellowstone National Park (YNP) contains the highest density of accessible thermal features in the world, with a range of pH and sulfur chemistry conditions uniquely suited to investigating sulfur transformations. YNP is situated above a large magma chamber situated at a depth greater than 8 km and derived from a mantle plume interacting with the North American Plate that extends down 660 km into the mantle transition zone, tilting at a 60 WNW angle [1, 2]. The last significant eruption in Yellowstone was 0.6 Ma, with the last known volcanic activity occurring 70,000 years ago; hydrothermal activity is related to the convection of meteoric water driven by heat associated with the cooling magma chamber . The hydrothermal reservoirs in this area are contained in units of alternating ash flow tuff and rhyolitic volcanic deposits, with the largest volumes of water likely contained in the relatively more permeable rhyolitic flows [2–4]. Hydrothermal water discharge in Yellowstone National Park is structurally controlled by fracture intersections associated with both the caldera rim and broader tectonic activity in the area [2, 3, 5]. Hydrothermal water composition is influenced by a combination of magmatic source fluids and gases, water-rock interactions, and water-vapor phase separation [3, 6]. Significant gaseous sulfur is released from magmatic sources in Yellowstone as H2S, with relatively little input of SO2[7–9]. The remarkable range of chemistries characterizing YNP thermal springs is due to a combination of these processes affecting the deep hydrothermal reservoir and the mixing with relatively oxic shallow meteoric waters.
Various water types identified at YNP are ultimately produced from a single deep sulfide-rich parent water body. According to Truesdell et al.  the chloride and sulfate concentrations and temperature of this water body are 8.74 mmol L-1, 115 μmol L-1, and 360°C, respectively. Fournier  proposed values of 11.3 mmol L-1 for Cl- concentration and 335–340°C temperature. In this work we based our quantitative interpretation of the data on the model proposed by Truesdell et al. . This model is based on the observation that composition of most YNP thermal waters can be explained by steam loss during adiabatic cooling of mixtures of a single deep parent water body with shallow, cold waters. Water-rock interactions as well as dilution and boiling alter water composition. Mixing with oxygen-rich subsurface waters leads to oxidation of hydrogen sulfide to sulfuric acid, which was proposed to occur by both abiotic and microbial processes [9, 10].
In addition to sulfate, some pools, especially those with circumneutral pH, contain high concentrations of thiosulfate , the intermediate product of sulfide oxidation by oxygen . Zero-valent sulfur (ZVS) is another common constituent of hydrothermal springs and pools, being more common in acidic sulfate spring ( and references therein). A source of ZVS may be either oxidation of hydrogen sulfide [13, 14] or from buried solfataras .
Sulfur cycling in geothermal springs is known to be mostly abiotic [9, 15] though bacterial sulfate reduction occurs in sediments surrounding these springs. Microbial sulfate reduction rates in sediments of the Cinder Pool have been reported to be ≤2 nmol cm-3 day-1)  and 1 ± 0.27 nmol cm-3 day-1. Even higher rates of microbial sulfate reduction were reported in the same works for microbial mats in run-off streams. Most of the microorganisms detected in YNP hot springs, including Cinder Pool, rely on hydrogen, rather than on hydrogen sulfide oxidation , in agreement with thermodynamic considerations presented in the same publication.
Springs also vary in complexity of sulfur cycling. Some pools discharge in outflow channels, that become oxygenated and cool rapidly, other springs form deep pools, which retain water for long times, resulting in a more complex sulfur cycle. Cinder Pool is the best studied of the pools with a complex sulfur cycle in surface waters . The bottom of the Cinder Pool hosts a lens of molten sulfur. Hydrolysis of molten sulfur generates hydrogen sulfide, thiosulfate and sulfate. The surface of Cinder Pool is partially covered with the hollow “cinders”, which consist of elemental sulfur and pyrite. Oxidation of thiosulfate to tetrathionate with atmospheric oxygen is catalyzed by pyrite .
Quantitatively, zero-valent sulfur is one of the most important sulfide oxidation intermediates in the Yellowstone National Park springs and pools. Zero-valent sulfur may be present in natural aquatic systems as four main forms: 1) solid (most probably crystalline, rhombic sulfur), colloidal sulfur , sulfur in the form of dissolved polysulfides [19–23], and dissolved (mostly cycloocta-) sulfur [24–26]. Existing analytical techniques for analysis of sulfur speciation in natural aquatic systems [27–29] as well as in understanding thermodynamics of water-sulfur-hydrogen sulfide system [30, 31], allows estimation and in some cases, precise quantification of ZVS speciation.
Thiocyanate is another sulfur species of interest in hydrothermal systems. Thiocyanate was detected in geothermal springs . It was also found at the Red Sea Atlantic II brine in concentrations ranging from 23 to 40 μmol L-1. Reactions between abiotic hydrogen cyanide and reduced sulfur species were proposed as the source of thiocyanate in the brine. Such reactions are well studied under controlled conditions [29, 34–37] as well as in salt marsh sediments .
Analysis of quadruple sulfur isotope fractionation is an important novel tool, which was successfully applied to biogeochemical sulfur cycling in modern marine  and limnic [40–42] aquatic systems as well as for understanding of sulfur cycling in ancient oceans [43–47]. Small variations in mass-dependent fractionation of multiple sulfur isotopes were proposed to result from complex mass-flow of this element through enzymatic systems . A limited number of studies were performed on multiple sulfur isotope fractionation during sulfate reduction, sulfur disproportionation and sulfide oxidation. Multiple sulfur isotope fractionation during abiotic sulfur cycling in natural aquatic systems was not studied and thus was neglected in previous studies.
In this paper we present a combined study of concentrations of various sulfur species as well as of quadruple sulfur isotope composition of major (sulfate and sulfide) sulfur species and of elemental sulfur in the YNP hydrothermal springs and pools. This study was performed in order to understand multiple sulfur isotope fractionation in the system, which is dominated by a mostly abiotic oxidative part of the sulfur cycle, and to understand how isotopic signals are transformed in sediment. Such data should bear on understanding of sulfur cycling in modern aquatic systems and ancient ocean.
Temperature, pH and concentrations of sulfur species in YNP springs
CP, 0 m depth
CP, 4 m depth
CP, 8 m depth
CP, 12 m depth
CP, 16 m depth
5659 ± 40
520 ± 5
1276 ± 15
984 ± 8
32.8 ± 0.8
52.3 ± 1.1
121.1 ± 8.3
15.5 ± 1.0
30.7 ± 2.2
28.7 ± 5.4
26.2 ± 0.8
26.1 ± 1.1
22.8 ± 1.0
25.5 ± 0.7
461 ± 205
91.7 ± 8.6
5.52 ± 3.32
3005 ± 504
0.148 ± 0.063
14.8 ± 2.7
15.1 ± 3.0
13.3 ± 0.5
47.9 ± 49.1
27.2 ± 7.5
5.29 ± 0.19
2.94 ± 0.06
3.79 ± 0.36
559 ± 25
4.02 ± 2.34
127.8 ± 4.6
125.6 ± 1.8
125.2 ± 2.5
129.9 ± 3.9
126.1 ± 1.4
0.365 ± 0.090
0.584 ± 0.060
1.47 ± 0.08
2.57 ± 0.70
1.07 ± 0.60
4.49 ± 0.49
4.76 ± 0.07
3.30 ± 1.84
4.49 ± 0.38
4.16 ± 0.39
0.192 ± 0.013
0.164 ± 0.005
0.196 ± 0.013
0.240 ± 0.030
Frying Pan Spring (FPL) and a small spring situated between FPL and the road (FPS)
Isotopic composition of sulfur species
CP, 0 m depth
CP, 4 m depth
CP, 8 m depth
CP, 12 m depth
CP, 16 m depth
The δ34S values of ZVS (0.28–0.43‰) are slightly lower than those of sulfide. This can be interpreted as an absence of equilibration between sulfide and ZVS through polysulfide formation (Eq. 1, ) due to the low pH, at which polysulfide anions are not stable with respect to hydrogen sulfide and rhombic S8.
Concentrations of sulfite were 365 ± 90 nmol L-1 in FPL and 584 ± 60 nmol L-1 in FPS. Concentrations of thiosulfate were 5.29 ± 0.19 μmol L-1 in FPL and 2.94 ± 0.06 μmol L-1 in FPS. Low concentrations of sulfite and thiosulfate may be explained by fast oxidation of sulfite ) and low stability of thiosulfate under low pH conditions .
Thiocyanate was detected in both FPL and FPS in sub-micromolar concentrations, 192 ± 13 nmol L-1 and 164 ± 5 nmol L-1, respectively.
Sulfur speciation and isotopic composition in the surface sediments collected near FPS and FPL
[AVS], mmol kg-1 (wt% as FeS)
[S-Py], mmol kg-1 (wt% as FeS2), ‰
-0.983 ± 0.025
-0.028 ± 0.007
-0.298 ± 0.140
[S0]-tot, mmol kg-1 (wt% as S0), ‰
-0.048 ± 0.008
-0.030 ± 0.007
-0.371 ± 0.035
[AVS], mmol kg-1 (wt% as FeS), ‰
-1.053 ± 0.021
-0.128 ± 0.038
0.216 ± 0.289
[S-Py], mmol kg-1 (wt% as FeS2), ‰
0.290 ± 0.016
-0.040 ± 0.001
-0.219 ± 0.140
[S0]-tot, mmol kg-1 (wt% as S0), ‰
0.701 ± 0.067
-0.043 ± 0.015
-0.139 ± 0.053
The surface sediment collected near the FPL was richer in sulfur species (430 μmol kg-1 AVS, 227 mmol kg-1 of non-ZVS CRS, and 6.57 mol kg-1 of ZVS (>21 weight%) (Table 3). Non-ZVS CRS and ZVS were somewhat heavier than in the FPL, but differences between their δ34S values were less than 1‰. AVS in the sediment was lighter than non-ZVS CRS by 1.34‰, and lighter than ZVS by 1.75‰. This difference may result either from AVS oxidation or from bacterial sulfate reduction in the sediment. Microbial sulfate reduction is known to produce isotope fractionations characterized by the difference between δ34S values between sulfide and sulfate of –70‰–3‰ [52, 53]. Microbial hydrogen sulfide oxidation to sulfate produces fractionation with δ34S in the range of –2‰–3‰ [54 and references therein]. Microbial oxidation of sulfide to elemental sulfur produces zero-valent sulfur, which is heavier by 0-3‰ [54 and references therein], although the same or even slightly larger isotope fractionation may be produced by equilibration between sulfur and hydrogen sulfide due to formation of soluble polysulfide species .
Dragon Spring (DS)
DS is another acidic spring with pH of 3.4, but is very different in chemical composition from FPL and FPS (Table 1). The sulfate concentration in this spring is 1276 ± 15 μmol L-1, and the chloride to sulfate molar ratio is 8.65. This molar ratio requires that a source of at least 78% of water in the spring is from the parent water body (Figure 3). Shallow sulfide oxidation is, therefore, responsible for production of 85–89% of sulfate in the spring waters (Figure 4). Overall sulfide concentration in the pool is relatively high, 121.1 ± 8.3 μmol L-1 (Table 1).
In this spring the sulfate isotope composition is very similar to that calculated using the mixing line based approach of Truesdell et al. : The measured value of δ34S is 2.44‰ and the calculated value is 2.38‰ (Figure 2, Table 2). In this pool, the Δ33S values of sulfide and sulfate were significantly different (-0.030‰ and–0.063‰, respectively). The origin of this pattern is unclear. While it has been demonstrated that bacterial sulfate reduction can produce variations in Δ33S, it also is possible that abiotic disproportionation reactions may produce variations in Δ33S.
The total ZVS concentration in the spring is quite low, 5.52 ± 3.32 μmol L-1, and only <1 μmol L-1 of it is in non-solid, cyanide-reactive form. Up to 4.9 μmol L-1 of ZVS may be dissolved in the DS water, 1.0 μmol L-1 as polysulfide sulfur and 3.9 μmol L-1 as dissolved elemental sulfur. Dissolved and polysulfide zero-valent sulfur pools are apparently not in equilibrium with rhombic sulfur and hydrogen sulfide.
The δ34S values of ZVS were close to those of sulfide (0.374‰) (Table 2). This observation confirms an absence of sulfide-polysulfide-sulfur equilibration at low pH values. The Δ33S values of ZVS (-0.072‰) were lower than those of sulfide and more similar to Δ33S value of sulfate.
The sulfite concentration in the DS water was 1.47 ± 0.08 1.0 μmol L-1, and the thiosulfate concentration was 3.79 ± 0.36 μmol L-1. Sub-micromolar concentrations of thiocyanate were detected in DS waters as well (196 ± 13 nM).
In the surface sediments collected at the outflow of DS, both acid volatile AVS and non-ZVS CRS (pyrite) contents were too low for isotopic analysis (<400 μmol kg-1). Sedimentary ZVS was found to be the main sulfur pool (5.66 mol kg-1, 18.1 weight%), and its δ34S values were <1‰ lower than the values for ZVS in the water column.
Isotopic composition of sulfur species in the sediment core taken in the vicinity of DS
CRS (non ZVS)
-5.035 ± 0.018
-0.076 ± 0.031
-0.098 ± 0.133
-6.413 ± 0.017
0.012 ± 0.016
-0.461 ± 0.148
-1.509 ± 0.011
-0.046 ± 0.020
-0.383 ± 0.172
-4.172 ± 0.015
-0.025 ± 0.019
-0.294 ± 0.271
-5.352 ± 0.004
-0.031 ± 0.013
-0.539 ± 0.153
-2.029 ± 0.006
-0.027 ± 0.015
-0.411 ± 0.144
-1.603 ± 0.016
-0.026 ± 0.010
-0.276 ± 0.114
-2.749 ± 0.010
-0.021 ± 0.008
-0.333 ± 0.041
-0.255 ± 0.006
-0.042 ± 0.004
-0.061 ± 0.040
-2.356 ± 0.013
-0.062 ± 0.015
-0.283 ± 0.192
-2.051 ± 0.001
-0.029 ± 0.008
-0.224 ± 0.144
0.034 ± 0.006
-0.043 ± 0.014
-0.218 ± 0.069
Cinder Pool (CP)
These findings are in apparent contradiction with results of voltammetric measurements of sulfur species concentrations in CP, which show high variability of sulfide and thiosulfate concentrations with time. Possible explanations of this discrepancy are different spatial and temporal scales of sampling. Tens of milliliters of water were sampled for analysis of concentrations of sulfur species by chromatographic and photometric techniques. Sampling was performed with a peristaltic pump, and we cannot refute the possibility that tubing inlet shifted during the sampling. On the other hand voltammetric measurements were done in situ, with the sampling volume surrounding a 100 micron diameter working electrode surface being on the order of hundreds of microliters, and a full scan requiring between 3 and 5 seconds depending on instrument settings. Voltammetry thus records heterogeneity in the water column in ways that sampling for chromatography and spectroscopy did not.
In one of the samples from 12 m depth, the whole cinder was sampled together with the water sample, and the ZVS concentration measured by chloroform extraction was as high as 816 μmol L-1 (this sample was omitted from the calculations and graph). In the surface waters of CP, solid sulfur is the main ZVS pool (11.5 μmol L-1). Concentration of cyanide-reactive ZVS is only 3.3 μmol L-1. As much as 3.5 μmol L-1 of ZVS can be solubilized in pool waters as cyclooctasulfur and polysulfides (Figure 5). Thus solid sulfur in the pool is close to equilibrium with dissolved cyclooctasulfur.
Chloride and sulfate concentrations in Cinder Pool have shown that at least 87% of water in the spring is the parent water body (Figure 3) and 54-63% of sulfate is produced by shallow sulfide oxidation. Interestingly, that the δ34S values for sulfate (7.96–9.09‰) are 3.5–4.8‰ higher than the values calculated using a mixing line based on the approach of . The reason for this discrepancy may root in an additional fractionation due to the processes involving oxidation of sulfur intermediates (e.g. ZVS, thiosulfate, tetrathionate). Schoen and Rye  attributed fractionation between sulfate and sulfide of up to 3.9‰ to abiotic sulfide oxidation to elemental sulfur followed by biological sulfur to sulfate oxidation. Additional support for this explanation may be found in the sulfur isotope fractionation between ZVS and sulfide. In the Cinder Pool water column, δ34S values of ZVS are 0.22–2.04‰ higher than δ34S values of sulfide, indicating partial equilibration between these two species through polysulfide formation (Table 2). Relatively large differences between sulfide and sulfate Δ33S values (up to –0.051‰) (Figure 7b) are likely to be a result of complex sulfur cycling in the pool, involving multiple turnover steps for each sulfur atom . The direction of the Δ33S fractionation, opposite to that produced by bacterial sulfate reduction, but consistent with abiotic disproportionation reactions, supports an assumption that it is produced by a complex redox processes network, possibly not enhanced by microbial sulfur reduction in the water column .
Concentrations of sulfite in the pool were found to be 3.3–4.8 μmol L-1. Concentrations of thiosulfate, 125–130 μmol L-1, were higher than in the previously discussed pools, possibly due to a longer water residence time, more complicated sulfur cycle, and less acidic conditions. Concentration of thiocyanate in the pool was as high as 595 nmol L-1 (at 4 m depth).
In the “cinders”, the δ34S values for non-ZVS CRS (pyrite) and ZVS are 0.597 and 0.879‰, respectively. The values of pyrite and ZVS δ34S in the molten sulfur on the bottom of the pool are 0.540‰ and 0.783‰, respectively. These values are quite similar to the δ34S values of sulfide and ZVS of dispersed sulfur at the surface of the pool itself (0.228 and 0.750‰, respectively). The Δ33S values of non-ZVS CRS and ZVS of these “cinders” and the molten bottom of the pool (-0.036--0.062‰) were slightly higher than the values for sulfide and ZVS at the surface of the pool (-0.064‰ and–0.093‰, respectively).
Evening Primrose Spring (EPS)
The water column of Evening Primrose Spring was found to be slightly acidic (pH = 6.5). The values of chloride and sulfate concentrations in the EPS show that at least 85% of water in the spring is from the parent water body (Figure 3) and 77-81% of the sulfate is produced by shallow sulfide oxidation. In this spring, the sulfate isotope composition agrees well with a composition calculated from the mixing line reported in : the measured value of δ34S is 4.68‰ and the calculated value is 3.91‰. The Δ33S values of sulfide and sulfate do not significantly differ. ZVS and thiosulfate were present in this pool in extraordinary high concentrations, 3 mmol L-1 and 559 μmol L-1, respectively. More than 99% of ZVS is solid sulfur dispersed in the water column of the pool. Concentration of colloidal sulfur in the pool is in the range of 20.6–26.5 μmol L-1. Up to 5 μmol L-1 and 1 μmol L-1 of cyclooctasulfur ZVS and of polysulfide ZVS, respectively, may be present as a soluble form in equilibrium with rhombic elemental sulfur. In contrast to that, the concentration of sulfite was relatively low, only 2.6 μmol L-1. Thiocyanate was also detected in this system at 240 nmol L-1 concentration.
ZVS δ34S values were 1‰ higher than sulfide values, showing that partial equilibration in sulfide–polysulfide–ZVS system was reached. Polysulfides are readily formed at the circumneutral pH and elevated temperatures . On the other hand, only 0.4% of the total ZVS in the EPS passes through a 5 μm pore-size filter. Large sulfur particles have a relatively small surface area, which slows its reaction with sulfide, according to reaction (1). The Δ33S value of ZVS was by 0.029‰ lower than the Δ33S value of sulfide.
The crust from the side of the pool was found to contain 0.15, 128 and 4134 mmol kg-1 of AVS, non-ZVS CRS and ZVS, respectively. AVS δ34S values were 2.75‰ lower than the values for sulfide in the pool. The δ34S values non-ZVS CRS and ZVS in the crust, equal to–0.242 and–0.244‰, respectively, were similar to the values of ZVS dispersed in the pool water (-0.402‰). The Δ33S value of AVS (-0.059‰) was slightly lower than those of non-ZVS CRS and ZVS (-0.024 and–0.026‰, respectively).
Ojo Caliente Spring (OCS)
In the Ojo Caliente Spring with circumneutral pH (7.4) at least 71% of water in the spring is from the parent water body (Figure 3) and less than 24% of the sulfide is oxidized in the shallow aquifer. In this pool, similar to the Cinder Pool, δ34S of sulfate varies from the calculated model value, but it is lighter by 5.6‰ than the calculated values. The Δ33S value of sulfate is only slightly (by 0.024‰) higher than that of sulfide. Currently, we have no explanation for the discrepancy between the measured value of sulfate δ34S and the one predicted by the model.
ZVS concentrations in the pool (148 nmol L-1) were too low to allow speciation or isotope analysis. Concentrations of sulfite and thiosulfate in the OCS were 1.07 μmol L-1 and 4.02 μmol L-1, respectively. The highest concentration of thiocyanate, 791 nmol L-1, was detected in OCS, that is higher compared to ZVS concentrations and only slightly lower than concentrations of sulfite and thiosulfate.
The main fraction of sulfur in all springs was not reactive toward cyanide (>59% in all pools). This means that the most abundant form of zero-valent sulfur in the YNP hydrothermal springs is solid orthorhombic sulfur. We explain it by coagulation and crystallization of colloidal sulfur, which may be produced initially by oxidation of hydrogen sulfide. This explanation is consistent with extremely low fractions of cyanide-reactive sulfur in pools with long water residence times (<1% in FPL and EPS). It is as well consistent with fast (in hours) coagulation and crystallization of colloidal sulfur in the Wadden Sea tidal flat pools that has been recently reported by Kamyshny and Ferdelman . In Cinder Pool, a relatively high fraction of zero-valent sulfur (19 ± 10%) is non-crystalline (e.g. not reactive toward hydrogen cyanide). This may be explained by a) fast turnover of sulfur in the water column of the pool, and b) the fact that one of important zero-valent sulfur reservoir in CP is cinders, which were excluded from sampling for total sulfur concentrations.
The sulfite concentration in all springs was <5 μmol L-1 (Table 1). These results agree with data of Zinder and Brock (1977). Sulfite concentrations detected by Zinder and Brock (1977) in CP, EPS, and OCS are 1.9, 9.4 and 0.0 μmol L-1, respectively. In this study we detected sulfite in the same systems at 4.5 ± 0.5, 2.6 ± 0.7 and 1.1 ± 0.6 μmol L-1, respectively. There are a number of reasons for low sulfite concentration in YNP springs. First, hydrogen sulfide and not sulfur dioxide is a primary volcanic gas in YNP. Second, rates of sulfite oxidation by oxygen are relatively high, especially at low pH . Third, at circumneutral and basic pH, sulfite will react with polysulfides and sulfur to form thiosulfate (reaction opposite to Eq. 4).
The concentration of hydrogen cyanide in the springs was not measured in the spring waters. Possibly it should be very low due to fast conversion of cyanide to thiocyanate by reactions with ZVS species, polythionates and thiosulfate [29, 34–37].
Thus thiocyanate may serve as a proxy for hydrogen cyanide concentration in the parent water body. Concentrations of thiocyanate in the pools were in the range of 0.16-0.79 μmol L-1 (Table 1). Interestingly, the highest concentration of thiocyanate was measured in the OCS, where it was more than 5 times higher than concentration of ZVS. Reaction between ZVS and hydrogen cyanide may be a contributing factor for low ZVS content in this spring.
Measurements of hydrogen cyanide concentrations in the waters of springs, especially of OCS, are required for understanding of its impact on geochemistry of hydrothermal springs and thiocyanate formation mechanisms.
Negative Δ33S of sulfide relative to sulfate was found by Johnston et al.  in cultures of sulfur disproportionating bacteria and attributed to backflow in the metabolism, which illustrates one case where Δ33S is depleted in product sulfide. This case also involved large fractionations, which we do not see here. We propose that complex multi-cycle transformations of reduced sulfur species including disproportionation reactions could be responsible for negative Δ33S difference between sulfide and sulfate in these pools and springs. One way to do this would be through disproportionation reactions that yielded fractionated products from an intermediate oxidation state reactant.
The Δ36S values were suggested to be helpful in interpretation of biogeochemical sulfur cycling . In our work, measurement of Δ36S values of sulfur species does not add to understanding of biogeochemical processes due to large standard deviation (±0.2‰). On the other hand, we suggest that performing more precise Δ36S measurements in future will possibly enable extracting independent information regarding mechanisms of sulfur cycling in hydrothermal systems.
Implications for sulfur isotope signatures in modern aquatic systems
Microbial sulfate reduction leads to isotope fractionation with differences in δ34S values for sulfide and sulfate ranging between –70‰ and 3‰ [40, 52, 53]. Absolute values of sulfur isotope fractionation have been shown to decrease with a decrease in sulfate concentration at [SO42-] below approximately 2 mmol L-1, and values between –30‰ and 3‰ are diagnostic for microbial sulfate reduction at sulfate concentration below 500 μmol L-1. In the sampled systems, only in the OCS is [SO42-] < 500 μmol L-1. Contrary to what is expected in microbially driven process, in this pool we have a relatively high difference, –7.14‰, between sulfide and sulfate. In the pool with highest sulfate concentration (FPL, [SO42-] = 5659 μmol L-1) fractionation between sulfide and sulfate is represented by δ34S value of 2.28‰. Thus, the model of Truesdell et al.  provides a sound explanation of the observed δ34S values of sulfate and sulfide in waters of hydrothermal springs and pools, and no fractionations that are diagnostic to microbial sulfur transformations are present.
These observations put forward the question regarding possible ways to differentiate between isotope fractionation by chemical processes at high temperatures and microbial processes at low sulfate concentration. We suggest that taking into account fractionation of three sulfur isotopes (32S, 33S, and 34S) of sulfate and sulfide may help in answering this question. A growing body of data is available on fractionation of multiple sulfur isotopes by microbial sulfate reduction, sulfur disproportionation , and sulfide oxidation [48, 55, 65]. On the other hand, the data on multiple sulfur isotope fractionation by chemical processes involving sulfur species both under controlled conditions and in natural systems is lacking. One of the aims of this work is to shed light on the multiple sulfur isotope fingerprints of complex, mostly abiotic, sulfur cycling.
Three processes may be responsible for formation of low-Δ33S sulfide: 1) deep water body sulfate-sulfide equilibration; 2) reoxidation of sulfide by mixing with oxic shallow waters; 3) geochemical processes during translation of isotopic composition of dissolved sulfide to sedimentary pyrite. As our dataset is quite limited (Figure 13), we can suggest only preliminary conclusions regarding the effects of these processes. In order to explore the first two scenarios, we addressed a connection between a) difference in Δ33S values of sulfide and sulfate in the waters of springs and pools; and b) the fraction of sulfate, which is produced by oxidation of sulfide at shallow depths or at the surface (Figure 11). We can see that in waters both unaffected by the surface sulfide oxidation (OCS), and waters affected by a complex surface sulfur cycle (CP) a relatively large negative difference between Δ33S values of sulfide and sulfate is detected. Noteworthy in the pools that are the most affected by surface sulfide oxidation (FPL, EPS), a difference between Δ33S values of sulfide and sulfate is slightly positive or zero. Thus both equilibration between sulfate and sulfide in the deep water body as well as complex surface sulfur cycling may lead to negative difference between Δ33S values of sulfide and sulfate. Unfortunately, our data set is not extensive enough to draw more solid conclusions. Shifts in Δ33S values of sulfide and sulfate during signal translation from water column to sediment also does not provide an unequivocal trend. In FPS and EPS sediments, as well as in sediments near the DS (possibly affected by microbial sulfate reduction), we detected a positive shift in the difference between Δ33S values of sulfide and sulfate. The opposite trend was observed for the FPL.
Although our work clearly shows that hydrothermal systems affected by shallow sulfide oxidation may produce relatively large negative difference between Δ33S values of sulfide and sulfate, our data set is not comprehensive enough to explain mechanisms of formation of this signal.
Implications for sulfur isotope signatures in geologic record
Multiple sulfur isotope fractionation in Archean ocean was affected not only by biological and abiotic transformations in aquatic systems, but by various arrays of mass-independent isotope fractionations produced by photolysis of atmospheric SO2, which originated from volcanic activity [68–71]. A mass-independent isotopic signals were superimposed on mass-dependent isotopic fractionation by microbial sulfate reduction [45, 72, 73], and possibly by microbial sulfur disproportionation . Other isotopic signals were formed by abiotic and biological transformations involving sulfur sources with mass-dependent sulfur isotope patterns that were not affected by atmospheric chemistry . To make this story even more complicated, the existence of mass-independent fractionation of sulfur isotopes by chemical reactions involving amino-acids and sulfate at high temperatures in the absence of UV radiation was recently discovered [74, 75].
Our research provides a framework for evaluation of multiple sulfur isotope signals produced by a complex net of abiotic transformations in a hydrothermal system. We have not observed any mass-independent isotope fractionation between sulfate and sulfide at elevated temperatures during migration of sulfate- and sulfide rich waters to the surface discharge. On the other hand, Δ33S values, which were relatively far from zero for mass-dependent fractionation (e.g. -0.128‰ for AVS in the solid samples near FPL, –0.093‰ for zero-valent sulfur in the CP water column, and –0.064‰ for sulfide in the CP water column), were detected. A combination of sulfate-sulfide isotope equilibration at high temperature, phase separation, and oxidative part of the sulfur cycle were found to produce δ34S fractionation values between sulfate and sulfide of up to 9‰. This signal may be transformed to even higher fractionation in the solid phase (e.g. in the Evening Primrose Spring difference between δ34S values of sulfate and sulfide are 3.3‰ and 4.9‰ in the water column and rocks surrounding the spring, respectively). Thus sulfur isotope fractionation input formed by sulfur transformations in hydrothermal systems may affect interpretation of multiple sulfur isotope fractionation patterns in Archean samples, especially in samples with δ34S values close to zero.
Understanding of mass-dependent multiple sulfur isotope fractionation on abiotic sulfur cycling may be especially helpful for the explanation of multiple sulfur isotope fractionation in the sulfide-rich Proterozoic ocean. In the Proterozoic, the concentration of atmospheric oxygen was high enough to prevent formation of mass-independent fractionation of sulfur isotopes, thus allowing detection of small mass dependent multiple sulfur isotope effects produced by microbial and chemical transformations of sulfur species. Existence of an intensive re-oxidative part of the sulfur cycle in the Proterozoic was proposed by results from two sulfur isotopes  and multiple sulfur isotope fractionation [43, 77, 78]. In the works of Johnston and co-authors [43, 77, 78], it was shown that microbial sulfur reduction in combination with microbial sulfur disproportionation explains the measured triple isotope composition (32S, 33S, 34S) of sulfur species better than microbial sulfate reduction alone. In fact, even the combination of microbial sulfate reduction and microbial sulfur disproportionation failed to explain the composition of some sedimentary sulfide samples with excessively low Δ33S values [77, 78]. This set of multiple sulfur isotope fractionation data from Yellowstone Natural Park provides observation of fractionation patterns with Δ33Ssulfide significantly smaller than Δ33Ssulfate in the natural aquatic system.
We sampled six springs with sulfide concentrations >15 μmol L-1: I) Frying Pan Spring (FPL; UTM coordinates: 521959 Easting, 4955476 Northing); II) Small spring between Frying Pan Spring and the road (FPS; UTM coordinates: 521959 Easting, 4955476 Northing); III) Dragon Spring (DS; UTM coordinates: 522877 Easting, 4953208 Northing); IV) Evening Primrose Spring (EPS; UTM coordinates: 517696 Easting, 4948724 Northing); V) Ojo Caliente Spring (OCS; UTM coordinates: 512863 Easting, 4934398 Northing); VI) Cinder Pool (CP; UTM coordinates: 523035 Easting, 4953286 Northing). Depth profile containing 5 points as well as isotopic composition of molten sulfur from the bottom of the pool and floating sulfur cinders were sampled at Cinder Pool as well.
All hydrothermal springs and pools water was sampled with a peristaltic environmental pump (Masterflex E/S Variable Speed Water Pump, EnviroTech, CA) and immediately preserved with 1/50 v/v 200 g/L zinc acetate solution for analysis of concentrations of sulfate, sulfide, thiocyanate and ZVStot. Molten sulfur at the bottom of Cinder Pool was as well sampled by a peristaltic environmental pump.
Solid phase was sampled either manually (surface sediments) or by insertion of a tube, constructed by cutting out of the bottom of 60 mL syringe, into the sediment.
Quantitative analytical methods and calculations
Analyses of total sulfide concentrations (protonated and deprotonated hydrogen sulfide and acid-soluble metal sulfides) were performed by a spectrophotometric technique according to , and sulfate concentrations were measured by ion-chromatography. Sulfite and thiosulfate were analyzed by monobromobimane derivatization followed by analysis by HPLC with a fluorescence detector . Derivatization was performed on-site immediately after sampling and samples were immediately frozen and shipped for analyses on dry ice. Speciation of ZVS was measured according to scheme proposed by . Total zero-valent sulfur (ZVStot) was extracted with chloroform and analyzed by HPLC with a UV-visible detector . Sum of all zero-valent sulfur species except solid sulfur (e.g. polysulfidic zero-valent sulfur as well as dissolved and colloidal elemental sulfur, ZVScyan) was analyzed by cyanolysis . Cyanolysis was performed on-site immediately after sampling. Samples were stored at ambient conditions prior to analyses. Individual polysulfide concentrations were not measured as the pH in most pools was too low for application of fast single phase derivatization with methyl trifluoromethanesulfonate .ZVS pools presented in Figure 5 were defined and calculated in the following way:
[ZVSelemental,dissolved,max,calcd] – concentration of S0 in cyclooctasulfur in equilibrium with rhombic sulfur calculated according to thermodynamic constants from .
[ZVSPS,max,calcd] –concentration of S0 in polysulfides in equilibrium with hydrogen sulfide and rhombic sulfur calculated according to thermodynamic constants from . Each polysulfide chain is supposed to have composition S2-S0n-1.
Preparation of samples and analysis for quadruple sulfur isotope composition for sulfide, sulfate and zero-valent sulfur was performed according to .
Conversion of sulfur species to silver sulfide for isotope analysis
Water samples for sulfide isotope composition analyses were filtered through a 0.4 μm Nuclepore Track-Etched Membrane filter within 6 hours of sampling with 20% zinc acetate preservation and frozen on dry ice. Filters with precipitate were boiled for 3 hours with 20 mL 5 N HCl in order to distill the acid-volatile sulfide (AVS).
BaCl2 solution was added to the filtrate in order to precipitate BaSO4. BaSO4 samples were reduced to H2S by boiling for 3 hours with 25 mL of a reagent prepared from 500 mL 36% HCl, 320 mL 57% HI, and 156 mL 85% H3PO4.
Zero-valent sulfur was extracted from 5–10 L aqueous samples with CHCl3 (3 × 1/50 v/v). Chloroform was evaporated on a rotor evaporator to 20–30 ml volume. Sample was transferred to the reduction reactor and remaining chloroform was evaporated under gentle flow of nitrogen. The zero-valent sulfur was converted to H2S by boiling for 2 h in a solution of 10 mL 12 mol L-1 HCl, 40 mL 99.98% ethanol, and 20 mL CrCl2 solution filtered through 0.20 μm filter . CrCl2 solution was prepared by stirring of the mixture of 208 g CrCl3 · 6H2O, 120 g Zn and 400 mL 0.5 N HCl for two hours under a gentle flow of nitrogen.
Sediments and cinders were shaken overnight with methanol (c.a. 30 ml MeOH per 1 g of wet sediment) in order to extract zero-valent sulfur. Methanol was evaporated on a rotor evaporator; sulfur was extracted from methanol-water mixture with dichloromethane. Dichloromethane extract was dried with anhydrous calcium chloride and evaporated in the reduction reactor under gentle flow of nitrogen.
After extraction of zero-valent sulfur with methanol sediment was subjected to AVS distillation, followed by CrCl2 reduction . This reaction was performed by boiling of sediment for 3 hours with 20 mL of CrCl2 reagent. This analysis accounts for non-zero-valent chromium-reducible sulfur, e.g. pyrite sulfur.
Hydrogen sulfide produced in all reactions was trapped in 5% zinc acetate solution. A portion of solution was used for measurement of concentrations, and to the rest of solution, silver sulfide solution was added in order to convert zinc sulfide to silver sulfide. Each Ag2S sample was aged for at least 1 week and cleaned with sequential washes of 200 mL Milli-Q water, rinsed in 50 mL of 1 M NH4OH overnight, then washed with 150 mL Milli-Q water and dried at 55°C.
Quadruple sulfur isotope analysis
Conversion of silver sulfide to sulfur hexafluoride and analysis for quadruple sulfur isotope composition was performed according to . Silver sulfide was quantitatively converted to SF6 by reaction with 10-fold excess of F2 gas at ~250°C for ~8 hours in Ni bombs. SF6 was purified cryogenically by distillation in an ethanol slurry at -115ºC, and by gas chromatography on a 12’ molecular sieve 5 Å/Hasep Q column with thermal conductivity detector. The isotopic abundance of the purified SF6 was analyzed on a Finnigan MAT 253 dual inlet mass spectrometer at m/e-values of 127, 128, 129, and 131 (32SF5+, 33SF5+, 34SF5+, and 36SF5+). Analytical uncertainties on sulfur isotope measurements, estimated from long-term reproducibility of Ag2S fluorinations, are 0.14, 0.008, and 0.20 (1σ) for δ34S, Δ33S, and Δ36S, respectively. Typical standard deviations between analyses of δ34S, Δ33S and Δ36S were 0.02–0.05‰, 0.01-0.02‰, and 0.1–0.2‰, respectively.
Isotopic measurements are reported on a VCDT scale, assuming the composition of IAEA S-1 is δ34S = -0.30‰, Δ33S = 0.094‰, and Δ36S = -0.70‰.
Isotopic composition notations
where 3XRVCDT is the isotopic ratio of the starting sulfide relative to Vienna-Cañon Diabolo Troilite (VCDT).
We define the fractionation between two sulfur species, A and B, as δ34SA-δ34SB.
Funding for this research was provided by Marie Curie Outgoing International Fellowship SULFUTOPES number POIF-GA-2008-219586 (to A. K.), NASA Exobiology and Evolutionary Biology program, solicitation number NNH08ZDA001N-EXOB (to A. K., G. K. D. and J. F.). GKD additionally acknowledges funding support from NSF EAR CAREER grant 1304352 and NSF MRI grant 0922961. The Research Permit office at Yellowstone National Park (NPS permit YELL-05474) is gratefully acknowledged for arranging sampling permission at Cinder Pool. Alyssa Findlay (UVM) and Don Nuzzio (Analytical Instrument Systems) are gratefully acknowledged for field, lab, and instrumental assistance.
- Smith RB, Jordan M, Steinberger B, Puskas CM, Farrell J, Waite GP, Husen S, Chang WL, O’Connell R: Geodynamic of the Yellowstone hotspot and mantle plume: Seismic and GPS imaging, kinematics, and mantle flow. J Volcanol Geoth Res. 2009, 188: 26-56. 10.1016/j.jvolgeores.2009.08.020.View ArticleGoogle Scholar
- Fournier RO: Geochemistry and dynamics of the Yellowstone national park hydrothermal system. Geothermal Biology and Geochemistry in Yellowstone National Park. Proceeding of the Thermal Biology Institute workshop. Edited by: Inskeep WP, McDermott TR. 2005, Bozeman: Montana State University Publications, 3-30.Google Scholar
- Fournier RO: Geochemistry and dynamics of the Yellowstone National Park hydrothermal system. Annu Rev Earth Pl Sc. 1989, 17: 13-53. 10.1146/annurev.ea.17.050189.000305.View ArticleGoogle Scholar
- Morgan LA, Shanks WC: Influences of rhyolitic lava flows on hydrothermal processes in Yellowstone Lake and on the Yellowstone plateau. Geothermal Biology and Geochemistry in Yellowstone National Park. Proceeding of the Thermal Biology Institute workshop. Edited by: Inskeep WP, McDermott TR. 2005, Bozeman: Montana State University Publications, 31-52.Google Scholar
- Christiansen RL: Geology of Yellowstone National Park. The Quaternary and Pliocene Yellowstone Plateau Volcanic field of Wyoming, Idaho, and Montana. U.S. Geological Survey Professional Paper 729-G. 2001Google Scholar
- Nordstrom DK, McCleskey RB, Ball JW: Sulfur geochemistry of hydrothermal waters in Yellowstone National Park: IV Acid-sulfate waters. Appl Geochem. 2009, 24: 191-207. 10.1016/j.apgeochem.2008.11.019.View ArticleGoogle Scholar
- Werner C, Hurwitz S, Evans WC, Lowenstern JB, Bergfeld D, Heasler H, Jaworowski C, Hunt A: Volatile emissions and gas geochemistry of Hot Spring Basin, Yellowstone National Park, USA. J Volcanol Geoth Res. 2008, 178: 751-762. 10.1016/j.jvolgeores.2008.09.016.View ArticleGoogle Scholar
- Zinder S, Brock TD: Sulfur dioxide in geothermal waters and gases. Geochim Cosmochim Acta. 1977, 41: 73-79. 10.1016/0016-7037(77)90187-9.View ArticleGoogle Scholar
- Xu Y, Schoonen MAA, Nordstrom DK, Cunningham KM, Ball JW: Sulfur geochemistry of hydrothermal waters in Yellowstone National Park: I. The origin of thiosulfate in hot spring waters. Geochim Cosmochim Acta. 1998, 62: 3729-3743. 10.1016/S0016-7037(98)00269-5.View ArticleGoogle Scholar
- Truesdell AH, Rye RO, Whelan JF, Thompson JM: Sulfate chemical and isotopic patterns in thermal waters of Yellowstone Park Wyoming. USGS Open File Report. 1978, 78–701: 435-436.Google Scholar
- Fry B, Ruf W, Gest H, Hayes JM: Sulfur isotope effects associated with oxidation of sulfide by O2 in aqueous solution. Chem Geol. 1988, 73: 205-210.Google Scholar
- Zhang J-Z, Millero FJ: The products from the oxidation of H2S in seawater. Geochim Cosmochim Acta. 1993, 57: 1705-1718. 10.1016/0016-7037(93)90108-9.View ArticleGoogle Scholar
- Chen KY, Morris JC: Kinetics of oxidation of aqueous sulfide by O2. Environ Sci Technol. 1972, 6: 529-537. 10.1021/es60065a008.View ArticleGoogle Scholar
- Yao W, Millero FJ: Oxidation of hydrogen sulfide by hydrous Fe(III) oxides in seawater. Mar Chem. 1996, 52: 1-16. 10.1016/0304-4203(95)00072-0.View ArticleGoogle Scholar
- Xu Y, Schoonen MAA, Nordstrom DK, Cunningham KM, Ball JW: Sulfur geochemistry of hydrothermal waters in Yellowstone National Park, Wyoming, USA. II. Formation and decomposition of thiosulfate and polythionate in Cinder Pool. J Volcanol Geoth Res. 2000, 97: 407-423. 10.1016/S0377-0273(99)00173-0.View ArticleGoogle Scholar
- Fishbain S, Dillon JG, Gough HL, Stahl DA: Linkage of high rates of sulfate reduction in Yellowstone hot springs to unique sequence types in the dissimilatory sulfate respiration pathway. Appl Environ Microb. 2003, 69: 3663-3667. 10.1128/AEM.69.6.3663-3667.2003.View ArticleGoogle Scholar
- Roychoudhury AN: Sulfate respiration in extreme environments: A kinetic study. Geomicrobiol J. 2004, 21: 33-43. 10.1080/01490450490253446.View ArticleGoogle Scholar
- Spear JR, Walker JJ, McCollom TM, Pace NR: Hydrogen and bioenergetics in the Yellowstone geothermal ecosystem. Proc Natl Acad Sci U S A. 2005, 102: 2555-2560. 10.1073/pnas.0409574102.View ArticleGoogle Scholar
- Kamyshny A, Ferdelman TG: Dynamics of zero-valent sulfur species including polysulfides at seep sites on intertidal sand flats (Wadden Sea, North Sea). Mar Chem. 2010, 121: 17-26. 10.1016/j.marchem.2010.03.001.View ArticleGoogle Scholar
- Boulegue J, Lord CJ: Church TM Sulfur speciation and associated trace metals (Fe, Cu) in the pore waters of Great Marsh, Delaware. Geochim Cosmochim Acta. 1982, 46: 453-464. 10.1016/0016-7037(82)90236-8.View ArticleGoogle Scholar
- Kamyshny A, Zilberbrand M, Ekeltchik I, Voitsekovski T, Gun J, Lev O: Speciation of polysulfides and zero-valent sulfur in sulfide-rich water wells in southern and central Israel. Aquat Geochem. 2008, 14: 171-192. 10.1007/s10498-008-9031-6.View ArticleGoogle Scholar
- Holmkvist L, Kamyshny A, Vogt C, Vamvakopoulos K, Ferdelman TG, Jørgensen BB: Sulfate reduction below the sulfate-methane transition in Black Sea sediments. Deep-Sea Res Pt I. 2011, 58: 493-504.View ArticleGoogle Scholar
- Lichtschlag A, Kamyshny A, Ferdelman TG, de Beer D: Intermediate sulfur oxidation state compounds in the euxinic surface sediments of the Dvurechenskii mud volcano (Black Sea). Geochim Cosmochim Acta. 2013, 105: 130-145.View ArticleGoogle Scholar
- Boulegue J: Solubility of elemental sulfur in water at 298 K. Phosphorus and Sulfur. 1978, 5: 127-128.View ArticleGoogle Scholar
- Wang F, Tessier A: Zero-valent sulfur and metal speciation in sediment porewaters of freshwater lakes. Environ Sci Technol. 2009, 43: 7252-7257. 10.1021/es8034973.View ArticleGoogle Scholar
- Kamyshny A: Solubility of cyclooctasulfur in pure water and sea water at different temperatures. Geochim Cosmochim Acta. 2009, 73: 6022-6028. 10.1016/j.gca.2009.07.003.View ArticleGoogle Scholar
- Fahey RC, Newton GL: Determination of low-weight thiols using monobromobimane fluorescent labeling and high-performance liquid chromatography. Method Enzymol. 1987, 143: 85-96.View ArticleGoogle Scholar
- Kamyshny A, Ekeltchik I, Gun J, Lev O: Method for the determination of inorganic polysulfide distribution in aquatic systems. Anal Chem. 2006, 78: 2631-2639. 10.1021/ac051854a.View ArticleGoogle Scholar
- Kamyshny A: Improved cyanolysis protocol for detection of zero-valent sulfur in natural aquatic systems. Limnol Oceanogr: Meth. 2009, 7: 442-448.View ArticleGoogle Scholar
- Kamyshny A, Goifman A, Gun J, Rizkov D, Lev O: Equilibrium distribution of polysulfide ions in aqueous solutions at 25°C: a new approach for the study of polysulfides’ equilibria. Environ Sci Technol. 2004, 38: 6633-6644. 10.1021/es049514e.View ArticleGoogle Scholar
- Kamyshny A, Gun J, Rizkov D, Voitsekovski T, Lev O: Equilibrium distribution of polysulfide ions in aqueous solutions at different temperatures by rapid single phase derivatization. Environ Sci Technol. 2007, 41: 6633-6644. 10.1021/es072617n.View ArticleGoogle Scholar
- Mukhin L: Evolution of organic compounds in volcanic regions. Nature. 1974, 251: 50-51. 10.1038/251050a0.View ArticleGoogle Scholar
- Dowler MJ, Ingmanson DE: Thiocyanate in Red Sea brine and its implications. Nature. 1979, 279: 51-52. 10.1038/279051a0.View ArticleGoogle Scholar
- Kelly DP, Chambers LA, Trudinger PA: Cyanolysis and spectrophotometric estimation of trithionate in mixture with thiosulfate and tetrathionate. Anal Chem. 1969, 41: 898-901. 10.1021/ac60276a029.View ArticleGoogle Scholar
- Szekeres L: Analytical chemistry of sulfur acids. Talanta. 1974, 21: 1-44. 10.1016/0039-9140(74)80061-5.View ArticleGoogle Scholar
- Luthy RG, Bruce SG: Kinetics of reaction of cyanide and reduced sulfur species in aqueous solutions. Environ Sci Technol. 1979, 13: 1481-1487. 10.1021/es60160a016.View ArticleGoogle Scholar
- Kamyshny A, Borkenstein CG, Ferdelman TG: Protocol for quantitative detection of elemental sulfur and polysulfide zero-valent sulfur distribution in natural aquatic samples. Geostand Geoanal Res. 2009, 33: 415-435. 10.1111/j.1751-908X.2009.00907.x.View ArticleGoogle Scholar
- Kamyshny A, Oduro H, Mansaray ZF, Farquhar J: Hydrogen cyanide accumulation and transformations in non-polluted salt marsh sediments. Aquat Geochem. 2013, 19: 97-113. 10.1007/s10498-012-9180-5.View ArticleGoogle Scholar
- Li X, Gilhooly WP, Zerkle AL, Lyons TW, Farquhar J, Werne JP, Varela R, Scranton MI: Stable sulfur isotopes in water column of the Cariaco Basin. Geochim Cosmochim Acta. 2010, 74: 6764-6778. 10.1016/j.gca.2010.08.020.View ArticleGoogle Scholar
- Canfield DE, Farquhar J, Zerkle AL: High isotope fractionation during sulfate reduction in a low-sulfate euxinic ocean analog. Geology. 2010, 38: 415-418. 10.1130/G30723.1.View ArticleGoogle Scholar
- Zerkle AL, Kamyshny A, Kump LR, Farquhar J, Oduro H, Arthur MA: Sulfur cycling in a stratified euxinic lake with moderately high sulfate: Constraints from quadruple S isotopes. Geochim Cosmochim Acta. 2010, 74: 4953-4970. 10.1016/j.gca.2010.06.015.View ArticleGoogle Scholar
- Kamyshny A, Zerkle AL, Mansaray Z, Ciglenečki I, Bura-Nakić E, Farquhar J, Ferdelman TG: Biogeochemical sulfur cycle in water column of shallow stratified sea-water lake: Speciation and quadruple sulfur isotope composition. Mar Chem. 2011, 127: 144-154. 10.1016/j.marchem.2011.09.001.View ArticleGoogle Scholar
- Johnston DT, Wing BA, Farquhar J, Kaufman AJ, Strauss H, Lyons TW, Kah LC, Canfield DE: Active microbial sulfur disproportionation in the Mesoproterozoic. Science. 2005, 310: 1477-1479. 10.1126/science.1117824.View ArticleGoogle Scholar
- Philippot P, Van Zuilen M, Lepot K, Thomazo C, Farquhar J, Van Kranendonk MJ: Early Archaean microorganisms preferred elemental sulfur, not sulfate. Science. 2007, 317: 1534-1537. 10.1126/science.1145861.View ArticleGoogle Scholar
- Shen YN, Farquhar J, Masterson A, Kaufman AJ, Buick R: Evaluating the role of microbial sulfate reduction in the early Archean using quadruple isotopic systematic. Earth Planet Sc Lett. 2009, 279: 381-393.View ArticleGoogle Scholar
- Wacey D, McLoughlin N, Whitehouse MJ, Kilburn MR: Two coexisting sulfur metabolisms in ca. 3400 Ma sandstone. Geology. 2010, 38: 1115-1118. 10.1130/G31329.1.View ArticleGoogle Scholar
- Farquhar J, Cliff J, Zerkle AL, Kamyshny A, Poulton SW, Claire M, Adams D, Harms B: Pathways for Neoarchean pyrite formation constrained by mass-independent sulfur isotopes. Proc Natl Acad Sci U S A. 2013, Early On-Line EditionGoogle Scholar
- Farquhar J, Johnston DT, Wing BA: Implications of conservation of mass effects on mass-dependent isotope fractionations: Influence of network structure on sulfur isotope phase space of dissimilatory sulfate reduction. Geochim Cosmochim Acta. 2007, 71: 5862-5875. 10.1016/j.gca.2007.08.028.View ArticleGoogle Scholar
- Amrani A, Kamyshny A, Lev O, Aizenshtat Z: Sulfur stable isotope distribution of polysulfide anions in (NH4)2Sn aqueous solution. Inorg Chem. 2006, 45: 1427-1429. 10.1021/ic051748r.View ArticleGoogle Scholar
- Zhang J-Z, Millero FJ: The rate of sulfite oxidation in seawater. Geochim Cosmochim Acta. 1991, 55: 677-685. 10.1016/0016-7037(91)90333-Z.View ArticleGoogle Scholar
- Davis RE: Displacement reactions at the sulfur atom. I. An interpretation of the decomposition of acidified thiosulfate. J Am Chem Soc. 1958, 80: 3565-3569. 10.1021/ja01547a018.View ArticleGoogle Scholar
- Harrison AG, Thode HG: Mechanism of the bacterial reduction of sulphate from isotope fractionation studies. Trans Faraday Soc. 1958, 54: 84-92.View ArticleGoogle Scholar
- Sim MS, Bosak T, Ono S: Large sulfur isotope fractionation does not require disproportionation. Science. 2011, 333: 74-77. 10.1126/science.1205103.View ArticleGoogle Scholar
- Brunner B, Bernasconi SM: A revised isotope fractionation model for dissimilatory sulfate reduction in sulfate reducing bacteria. Geochim Cosmochim Acta. 2005, 69: 4759-4771. 10.1016/j.gca.2005.04.015.View ArticleGoogle Scholar
- Johnston DT, Farquhar J, Wing BA, Kaufman A, Canfield DE, Habicht KS: Multiple sulfur isotope fractionations in biological systems: A case study with sulfate reducers and sulfur disproportionators. Am J Sci. 2005, 305: 645-660. 10.2475/ajs.305.6-8.645.View ArticleGoogle Scholar
- Cline JD: Spectrophotometric determination of hydrogen sulfide in natural waters. Limnol Oceanogr. 1969, 14: 454-458. 10.4319/lo.1969.14.3.0454.View ArticleGoogle Scholar
- Schoen R, Rye RO: Sulfur isotope distribution in solfataras, Yellowstone National Park. Science. 1970, 170: 1082-1084. 10.1126/science.170.3962.1082.View ArticleGoogle Scholar
- Gartman A, Yücel M, Madison AS, Chu DW, Ma S, Janzen CP, Becker EL, Beinart RA, Girgius PR, Luther GW: Sulfide oxidation across diffuse flow zones of hydrothermal vents. Aquat Geochem. 2011, 17: 583-601. 10.1007/s10498-011-9136-1.View ArticleGoogle Scholar
- Luther GW, Findlay AJ, MacDonald DJ, Owings SM, Hanson TE, Beinart RA, Girgius PR: Thermodynamics and kinetics of sulfide oxidation by oxygen: a look at inorganically controlled reactions and biologically mediated processes in the environment. Front Microbiol. 2011, 2: 62-View ArticleGoogle Scholar
- Johnston F, McAmish L: A study of rates of sulfur production in acid thiosulfate solutions using S-35. J Colloid Interf Sci. 1973, 42: 112-119. 10.1016/0021-9797(73)90013-1.View ArticleGoogle Scholar
- Xu Y, Schoonen MAA: The stability of thiosulfate in the presence of pyrite in low-temperature aqueous solutions. Geochim Cosmochim Acta. 1995, 59: 4605-4622. 10.1016/0016-7037(95)00331-2.View ArticleGoogle Scholar
- Hoefs J: Stable Isotope Geochemistry. 2009, Berlin, Heidelberg: Springer-Verlag, 6Google Scholar
- Farquhar J, Johnston DT, Wing BA, Habicht KS, Canfield DE, Airieau S, Thiemens MH: Multiple sulphur isotopic interpretations of biosynthetic pathways: implications for biological signatures in the sulphur isotope record. Geobiology. 2003, 1: 27-36. 10.1046/j.1472-4669.2003.00007.x.View ArticleGoogle Scholar
- Zeikus JG, Dawson MA, Thompson TE, Ingvorsen K, Hatchikian EC: Microbial ecology of volcanic sulphidogenesis: Isolation and characterization of Thermosulfobacterium commune gen. nov. and sp. nov. J Gen Microbiol. 1983, 129: 1159-1169.Google Scholar
- Zerkle AL, Farquhar J, Johnston DT, Cox RP, Canfield DE: Fractionation of multiple sulfur isotopes during phototrophic oxidation of sulfide and elemental sulfur by a green sulfur bacterium. Geochim Cosmochim Acta. 2009, 73: 291-306. 10.1016/j.gca.2008.10.027.View ArticleGoogle Scholar
- Johnston DT, Farquhar J, Habicht K, Canfield DE: Sulphur isotopes and the search for life: strategies for identifying sulphur metabolisms in the rock record and beyond. Geobiology. 2008, 6: 425-435. 10.1111/j.1472-4669.2008.00171.x.View ArticleGoogle Scholar
- Habicht KS, Gade M, Thamdrup B, Berg P, Canfield DE: Calibration of sulfate levels in the Archean Ocean. Science. 2002, 298: 2372-2374. 10.1126/science.1078265.View ArticleGoogle Scholar
- Farquhar J, Bao H, Thiemens M: Atmospheric influence on Earth’s earliest sulfur cycle. Science. 2000, 289: 756-758. 10.1126/science.289.5480.756.View ArticleGoogle Scholar
- Farquhar J, Savarino J, Airieau S, Thiemens MH: Observation of wavelength-sensitive mass-independent sulfur isotope effects during SO2 photolysis: Implications for the early atmosphere. J Geophys Res. 2001, 106: 32829-32839. 10.1029/2000JE001437.View ArticleGoogle Scholar
- Lyons JR: Atmospherically-derived mass-independent sulfur isotope signatures, and incorporation into sediments. Chem Geol. 2009, 267: 164-174. 10.1016/j.chemgeo.2009.03.027.View ArticleGoogle Scholar
- Philippot P, Van Zuilen M, Rollion-Bard C: Variations in atmospheric sulphur chemistry on early Earth linked to volcanic activity. Nat Geosci. 2012, 5: 668-674. 10.1038/ngeo1534.View ArticleGoogle Scholar
- Shen YN, Buick R: The antiquity of sulfate reduction. Earth Sci Rev. 2004, 64: 243-272. 10.1016/S0012-8252(03)00054-0.View ArticleGoogle Scholar
- Ueno Y, Ono S, Rumble D, Maruyama S: Quadruple sulfur isotope analysis of ca. 3.5 Ga Dresser Formation: New evidence of microbial sulfate reduction in the early Archean. Geochim Cosmochim Acta. 2008, 72: 5675-5691. 10.1016/j.gca.2008.08.026.View ArticleGoogle Scholar
- Watanabe Y, Farquhar J, Ohmoto H: Anomalous fractionations of sulfur isotopes during thermochemical sulfate reduction. Science. 2009, 324: 370-373. 10.1126/science.1169289.View ArticleGoogle Scholar
- Oduro H, Harms B, Sintim HO, Kaufman AJ, Cody G, Farquhar J: Evidence of magnetic isotope effects during thermochemical sulfate reduction. Proc Natl Acad Sci U S A. 2011, 108: 17635-17638. 10.1073/pnas.1108112108.View ArticleGoogle Scholar
- Ries JB, Fike DA, Pratt LM, Lyons TW, Grotzinger JP: Superheavy pyrite (δ34Spyr > δ34SCAS) in the terminal Proterozoic Nama Group, southern Namibia: A consequence of low seawater sulfate at the dawn of animal life. Geology. 2009, 37: 743-746. 10.1130/G25775A.1.View ArticleGoogle Scholar
- Johnston DT, Poulton SW, Fralick PW, Wing BA, Canfield DE, Farquhar J: Evolution of the oceanic sulfur cycle at the end of the Paleoproterozoic. Geochim Cosmochim Acta. 2006, 70: 5723-5739. 10.1016/j.gca.2006.08.001.View ArticleGoogle Scholar
- Johnston DT, Farquhar J, Summons RE, Shen Y, Kaufman AJ, Masterson AL, Canfield DE: Sulfur isotope biogeochemistry of the Proterozoic McArthur Basin. Geochim Cosmochim Acta. 2008, 72: 4278-4290. 10.1016/j.gca.2008.06.004.View ArticleGoogle Scholar
- Zopfi J, Böttcher ME, Jørgensen BB: Biogeochemistry of sulfur and iron in Thioploca-colonized surface sediments in the upwelling area off central Chile. Geochim Cosmochim Acta. 2008, 72: 827-843. 10.1016/j.gca.2007.11.031.View ArticleGoogle Scholar
- Forrest J, Newman L: Ag-110 microgram sulfate analysis for short time resolution of ambient levels of sulfur aerosol. Anal Chem. 1977, 49: 1579-1584. 10.1021/ac50019a030.View ArticleGoogle Scholar
- Gröger J, Franke J, Hamer K, Schulz HD: Quantitative recovery of elemental sulfur and improved selectivity in a chromium-reducible sulfur distillation. Geostand Geoanal Res. 2009, 33: 17-27. 10.1111/j.1751-908X.2009.00922.x.View ArticleGoogle Scholar
- Fossing H, Jørgensen BB: Measurement of bacterial sulfate reduction in sediments: Evaluation of a single-step chromium reduction method. Biogeochemistry. 1989, 8: 205-222.View ArticleGoogle Scholar
- Hulston JR, Thode HG: Cosmic ray produced 36S and 33S in metallic phase of iron meteorites. J Geophys Res. 1965, 70: 4435-4442. 10.1029/JZ070i018p04435.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.