Sulfur geochemistry of hydrothermal waters in Yellowstone National Park: IV Acid–sulfate waters

Many waters sampled in Yellowstone National Park, both high-temperature (30–94 °C) and low-temperature (0–30 °C), are acid–sulfate type with pH values of 1–5. Sulfuric acid is the dominant component, especially as pH values decrease below 3, and it forms from the oxidation of elemental S whose origin is H₂S in hot gases derived from boiling of hydrothermal waters at depth. Four determinations of pH were obtained: (1) field pH at field temperature, (2) laboratory pH at laboratory temperature, (3) pH based on acidity titration, and (4) pH based on charge imbalance (at both laboratory and field temperatures). Laboratory pH, charge imbalance pH (at laboratory temperature), and acidity pH were in close agreement for pH < 2.7. Field pH measurements were predominantly used because the charge imbalance was <±10%. When the charge imbalance was generally >±10%, a selection process was used to compare acidity, laboratory, and charge balance pH to arrive at the best estimate. Differences between laboratory and field pH can be explained based on Fe oxidation, H₂S or S₂O₃ oxidation, CO₂ degassing, and the temperature-dependence of pK₂ for H₂SO₄. Charge imbalances are shown to be dependent on a speciation model for pH values <3. The highest SO₄ concentrations, in the thousands of mg/L, result from evaporative concentration at elevated temperatures as shown by the consistently high δ¹⁸O values (−10‰ to −3‰) and a δD vs. δ¹⁸O slope of 3, reflecting kinetic fractionation. Low SO₄ concentrations (<100 mg/L) for thermal waters (>350 mg/L Cl) decrease as the Cl⁻ concentration increases from boiling which appears inconsistent with the hypothesis of H₂S oxidation as a source of hydrothermal SO₄. This trend is consistent with the alternate hypothesis of anhydrite solubility equilibrium. Acid–sulfate water analyses are occasionally high in As, Hg, and NH₃ concentrations but in contrast to acid mine waters they are low to below detection in Cu, Zn, Cd, and Pb concentrations. Even concentrations of SO₄, Fe, and Al are much lower in thermal waters than acid mine waters of the same pH. This difference in water chemistry may explain why certain species of fly larvae live comfortably in Yellowstone’s acid waters but have not been observed in acid rock drainage of the same pH.     

The oxidation state of sulfur in synthetic and natural glasses determined by X-ray absorption spectroscopy

Sulfur K-edge X-ray absorption near edge structure (XANES) spectra were recorded for experimental glasses of various compositions prepared at different oxygen fugacities (fO₂) in one-atmosphere gas-mixing experiments at 1400 °C. This sample preparation method only results in measurable S concentrations under either relatively reduced (log fO₂ < −9) or oxidised (log fO₂ > −2) conditions. The XANES spectra of the reduced samples are characterised by an absorption edge crest at 2476.4 eV, typical of S²⁻. In addition, spectra of Fe-bearing compositions exhibit a pronounced absorption edge shoulder. Spectra for all the Fe-free samples are essentially identical, as are the spectra for the Fe-bearing compositions, despite significant compositional variability within each group. The presence of a sulfide phase, such as might exsolve on cooling, can be inferred from a pre-edge feature at 2470.5 eV.The XANES spectra of the oxidised samples are characterised by an intense transition at 2482.1 eV, typical of the sulfate anion SO₄²⁻. Sulfite (SO₃²⁻) has negligible solubility in silicate melts at low pressures. The previous identification of sulfite species in natural glass samples is attributed to an artefact of the analysis (photoreduction of S⁶⁺). S⁴⁺ does, however, occur unambiguously with S⁶⁺ in Fe-free and Fe-poor compositions prepared in equilibrium with CaSO₄ at 4-16 kbar, and when buffered with Re/ReO₂ at 10 kbar. Solubility of S⁴⁺ thus requires partial pressures of SO₂ considerably in excess of 1 bar. A number of experiments were undertaken in an attempt to access intermediate fO₂s more applicable to terrestrial volcanism. Although these were largely unsuccessful, S²⁻ and S⁶⁺ were found to coexist in some samples that were not in equilibrium with the imposed fO₂.The XANES spectra of natural olivine-hosted melt inclusions and submarine glasses representative of basalts at, or close to, sulfide saturation show mainly dissolved S²⁻, but with minor sulfate, and additionally a peak at 2469.5 eV, which, although presumably due to immiscible sulfide, is 1 eV lower than that typical of FeS. These sulfate and sulfide-related peaks disappear with homogenisation of the inclusions by heating to 1200 °C followed by rapid quenching, suggesting that both these features are a result of cooling under natural conditions. The presence of small amounts of sulfate in otherwise reduced basaltic magmas may be explained by the electron exchange reaction: S²⁻ + 8Fe³⁺ = S⁶⁺ + 8Fe²⁺, which is expected to proceed strongly to the right with decreasing temperature. This reaction would explain why S²⁻ and S⁶⁺ are frequently found together despite the very limited fO₂ range over which they are thermodynamically predicted to coexist. The S XANES spectra of water-rich, highly oxidised, basaltic inclusions hosted in olivine from Etna and Stromboli confirm that nearly all S is dissolved as sulfate, explaining their relatively high S contents.     

Irrigation produces elevated arsenic in the underlying groundwater of a semi-arid basin in Southwestern Idaho

The shallow aquifer beneath the Western Snake River Plain (Idaho, USA) exhibits widespread elevated arsenic concentrations (up to 120 μg L⁻¹). While semi-arid, crop irrigation has increased annual recharge to the aquifer from approximately 1 cm prior to a current rate of >50 cm year⁻¹. The highest aqueous arsenic concentrations are found in proximity to the water table (all values >50 μg L⁻¹ within 50 m) and concentrations decline with depth. Despite strong vertical redox stratification within the aquifer, spatial distribution of aqueous species indicates that redox processes are not primary drivers of arsenic mobilization. Arsenic release and transport occur under oxidizing conditions; groundwater wells containing dissolved arsenic at >50 μg L⁻¹ exhibit elevated concentrations of O₂ (average 4 mg L⁻¹) and NO₃ (average 8 mg L⁻¹) and low concentrations of dissolved Fe (<20 μg L⁻¹). Sequential extractions and spectroscopic analysis of surficial soils and sediments indicate solid phase arsenic is primarily arsenate and is present at elevated concentrations (4–45 mg kg⁻¹, average: 17 mg kg⁻¹) relative to global sedimentary abundances. The highest concentrations of easily mobilized arsenic (up to 7 mg kg⁻¹) are associated with surficial soils and sediments visibly stained with iron oxides. Batch leaching experiments on these materials using irrigation waters produce pore water arsenic concentrations approximating those observed in the shallow aquifer (up to 152 μg L⁻¹). While As:Cl aqueous phase relationships suggest minor evaporative enrichment, this appears to be a relic of the pre-irrigation environment. Collectively, these data indicate that infiltrating irrigation waters leach arsenic from surficial sediments to the underlying aquifer.     

Sedimentary iron geochemistry in acidic waterways associated with coastal lowland acid sulfate soils

We examined the solubility, mineralogy and geochemical transformations of sedimentary Fe in waterways associated with coastal lowland acid sulfate soils (CLASS). The waterways contained acidic (pH 3.26-3.54), Fe^III-rich (27-138 μM) surface water with low molar Cl:SO₄ ratios (0.086-5.73). The surficial benthic sediments had high concentrations of oxalate-extractable Fe(III) due to schwertmannite precipitation (kinetically favoured by 28-30% of aqueous surface water Fe being present as the Fe^III species). Subsurface sediments contained abundant pore-water HCO₃ (6-20 mM) and were reducing (Eh < −100 mV) with pH 6.0-6.5. The development of reducing conditions caused reductive dissolution of buried schwertmannite and goethite (formed via in situ transformation of schwertmannite). As a consequence, pore-water Fe^II concentrations were high (>2 mM) and were constrained by precipitation-dissolution of siderite. The near-neutral, reducing conditions also promoted SO₄-reduction and the formation of acid-volatile sulfide (AVS). The results show, for the first time for CLASS-associated waterways, that sedimentary AVS consisted mainly of disordered mackinawite. In the presence of abundant pore-water Fe^II, precipitation-dissolution of disordered mackinawite maintained very low (i.e. <0.1 μM) S^−II concentrations. Such low concentrations of S^−II caused slow rates for conversion of disordered mackinawite to pyrite, thereby resulting in relatively low concentrations of pyrite (<300 μmol g⁻¹ as Fe) compared to disordered mackinawite (up to 590 μmol g⁻¹ as Fe). This study shows that interactions between schwertmannite, goethite, siderite, disordered mackinawite and pyrite control the geochemical behaviour of sedimentary Fe in CLASS-associated waterways.     

Pyrite dissolution in acidic media

Oxidation of pyrite in aqueous solutions in contact with air (oxygen 20%) was studied at 25°C using short-term batch experiments. Fe²⁺ and SO₄²⁻ were the only dissolved Fe and S species detected in these solutions. After a short period, R = [S]_tot/[Fe]_tot stabilized from 1.25 at pH = 1.5 to 1.6 at pH = 3. These R values were found to be consistent with previously published measurements (as calculated from the raw published data). This corresponds to a nonstoichiometric dissolution (R < 2) resulting from a deficit in aqueous sulfur. Thermodynamics indicate that S(−I) oxidation can only produce S_(s)⁰ and SO₄²⁻ under these equilibrium conditions. However, Pourbaix diagrams assuming the absence of SO₄²⁻ indicate that S₂O₃²⁻ and S₄O₆²⁻ can appear in these conditions. Using these species the simplest expected oxidation mechanism is

     

Landscape variations in stream water SO4 and δSSO4 in a boreal stream network

Despite reduced anthropogenic deposition during the last decades, deposition sulphate may still play an important role in the biogeochemical cycles of S and many catchments may act as net sources of S that may remain for several decades. The aim of this study is to elucidate the temporal and spatial dynamics of both SO₄²⁻ and δ³⁴S_SO4 in stream water from catchments with varying percentage of wetland and forest coverage and to determine their relative importance for catchment losses of S. Stream water samples were collected from 15 subcatchments ranging in size from 3 to 6780 ha, in a boreal stream network, northern Sweden. In forested catchments (<2% wetland cover) S-SO₄²⁻ concentrations in stream water averaged 1.7 mg L⁻¹ whereas in wetland dominated catchments (>30% wetland cover) the concentrations averaged 0.3 mg L⁻¹. A significant negative relationship was observed between S-SO₄²⁻ and percentage wetland coverage (r² = 0.77, p < 0.001) and the annual export of stream water SO₄²⁻ and wetland coverage (r² = 0.76, p < 0.001). The percentage forest coverage was on the other hand positively related to stream water SO₄²⁻ concentrations and the annual export of stream water SO₄²⁻ (r² = 0.77 and r² = 0.79, respectively). The annual average δ³⁴S_SO4 value in wetland dominated streams was +7.6‰ and in streams of forested catchments +6.7‰. At spring flood the δ³⁴S_SO4 values decreased in all streams by 1‰ to 5‰. The δ³⁴S_SO4 values in all streams were higher than the δ³⁴S_SO4 value of +4.7‰ in precipitation (snow). The export of S ranged from 0.5 kg S ha⁻¹ yr⁻¹ (wetland headwater stream) to 3.8 kg S ha⁻¹ yr⁻¹ (forested headwater stream). With an average S deposition in open field of 1.3 kg S ha⁻¹ yr⁻¹ (2002-2006) the mass balance results in a net export of S from all catchments, except in catchments with >30% wetland. The high temporal and spatial resolution of this study demonstrates that the reducing environments of wetlands play a key role for the biogeochemistry of S in boreal landscapes and are net sinks of S. Forested areas, on the other hand were net sources of S.     

Solubility of cyclooctasulfur in pure water and sea water at different temperatures

The solubility of cyclooctasulfur in water and sea water at various temperatures in the range between 4 and 80 °C was determined. Cyclooctasulfur in equilibrium with rhombic sulfur reacted with hot acidic aqueous potassium cyanide to form thiocyanate anion which was measured by anion chromatography. Sulfur solubility in pure water was found to increase with temperature by more than 78 times: from 6.1 nM S₈ at 4 °C to 478 nM S₈ at 80 °C. The following thermodynamic values for solubilisation of S₈ in water were calculated from the experimental data: K° = 3.01 ± 1.04 × 10⁻⁸, ΔG_r° = 42.93 ± 0.73 kJ mol⁻¹, ΔH_r° = 47.4 ± 3.6 kJmol⁻¹, ΔS_r° = 15.0 ± 11.7 J mol⁻¹ K⁻¹). Solubility of cyclooctasulfur in sea water was found to be 61 ± 13% of the solubility in pure water regardless of the temperature.     

Reduced sulfur species in a stratified seawater lake (Rogoznica Lake, Croatia); seasonal variations and argument for organic carriers of reactive sulfur

Over a period of a year, Hg⁰-reactive, total reduced sulfur species (RSS_T), as well as a non-volatile fraction that cannot be gas-stripped at pH ∼2 (RSS_NV), have been measured by voltammetry in a stratified, saline lake. In the hypolimnion, RSS_T is dominated by unusually high (up to 5 mM) dissolved divalent sulfur (S^−II), present as H₂S + HS⁻ and as inorganic polysulfides (H_xS_n^x−2). Less abundant RSS_NV is attributed to dissolved zero-valent sulfur (S⁰) in inorganic polysulfides. Assuming negligible contribution of organic S⁰ species in the hypolimnion, the equilibrium distribution of polysulfide ions is calculated; S₅²⁻ is found to predominate. In the epilimnion, all RSS_T consists of RSS_NV within analytical uncertainty. Through spring and summer, RSS_T and RSS_NV display little vertical or seasonal variation, but they increase dramatically when stratification breaks down in autumn. Based on decay rate, RSS during mixing events is attributed to dissolved S₈ from oxidation of sulfide and decomposition of inorganic polysulfides. This hypothesis quantitatively predicts precipitation of elemental sulfur in a year when colloidal sulfur was observed and predicts no precipitation in a year when it was not observed. Except during mixing events, the entire water column is undersaturated with respect to both rhombic sulfur and biologic sulfur, and the limited variations of RSS exclude hydrophobic and volatile aqueous S₈ as a major species. During such periods, RSS (typically 8 nM) may be associated with organic carbon, perhaps as adsorbed S₈ or as covalently bound polysulfanes or polysulfides. The hypolimnion is viewed as a zero-valent sulfur reactor that creates S⁰-containing, dissolved organic macromolecules during stable stratification periods. Some are sufficiently degradation-resistant and hydrophilic to be dispersed throughout the lake during mixing events, subsequently giving rise to ∼10⁻⁸ M RSS in the oxic water column. Voltammetrically determined RSS in oxic natural waters has often been described as “sulfide” or “metal complexed sulfide”, implying an oxidation state of S^−II; we argue that RSS in oxic Rogoznica Lake waters is mainly S⁰.     

Arsenic, iron and sulfur co-diagenesis in lake sediments

Profiles of porewater pH and dissolved As, Fe, Mn, sulfate, total sulfide (ΣS^−II), total zero-valent sulfur (ΣS⁰), organic carbon and major ion concentrations, as well as those of solid As, acid-volatile sulfide (AVS), total S, Fe, Mn, Al, organic C, ²¹⁰Pb and ¹³⁷Cs were determined in the sediment of four lakes spanning a range of redox and geochemical conditions. An inverse modeling approach, based on a one-dimensional transport-reaction equation assuming steady-state, was applied to the porewater As profiles and used to constrain the net rates of reactions involving As (). The model defines depth intervals where As is either released to (positive ) or removed from (negative ) the porewaters.At two of the sites, whose bottom water were oxygenated at sampling time, a production zone ( = 12 × 10⁻¹⁸ mol cm⁻³ s⁻¹-71 × 10⁻¹⁸ mol cm⁻³ s⁻¹) is inferred a few cm below the sediment-water interface, coincident with sharp porewater As and Fe peaks that indicate an intense coupled recycling of As and Fe. This process is confirmed by solid As and Fe maxima just below the sediment surface. In these two lakes a zone of As consumption ( = −5 × 10⁻¹⁸ mol cm⁻³ s⁻¹ to −53 × 10⁻¹⁸ mol cm⁻³ s⁻¹), attributed to the slow adsorption of As to authigenic Fe oxyhydroxides, occurs just above the production zone. A second-order rate constant of 0.12 ± 0.03 cm³ mol⁻¹ s⁻¹ is estimated for this adsorption reaction.Such features in the porewater and solid profiles were absent from the two other lakes that develop a seasonally anoxic hypolimnion. Thermodynamic calculations indicate that the porewaters of the four lakes, when sulfidic (i.e., ΣS^−II ? 0.1 μM), were undersaturated with respect to all known solid As sulfides; the calculation also predicts the presence of As^V oxythioanions in the sulfidic waters, as suggested by a recent study. In the sulfidic waters, the removal of As ( = −1 × 10⁻¹⁸ mol cm⁻³ s⁻¹ to −23 × 10⁻¹⁸ mol cm⁻³ s⁻¹) consistently occurred when saturation, with respect to FeS_(s), was reached and when As^V oxythioanions were predicted to be significant components of total dissolved As. This finding has potential implications for As transport in other anoxic waters and should be tested in a wider variety of natural environments.     

The role of sulfate groups in controlling CaCO3 polymorphism

The nucleation and growth of CaCO₃ phases from aqueous solutions with SO₄²⁻:CO₃²⁻ ratios from 0 to 1.62 and a pH of ∼10.9 were studied experimentally in batch reactors at 25 °C. The mineralogy, morphology and composition of the precipitates were characterized by X-ray diffraction, Fourier transform infrared spectroscopy, scanning electron microscopy and microanalyses. The solids recovered after short reaction times (5 min to 1 h) consisted of a mixture of calcite and vaterite, with a S content that linearly correlates with the SO₄²⁻:CO₃²⁻ ratio in the aqueous solution. The solvent-mediated transformation of vaterite to calcite subsequently occurred. After 24 h of equilibration, calcite was the only phase present in the precipitate for aqueous solutions with SO₄²⁻:CO₃²⁻ ? 1. For SO₄²⁻:CO₃²⁻ > 1, vaterite persisted as a major phase for a longer time (>250 h for SO₄²⁻:CO₃²⁻ = 1.62). To study the role of sulfate in stabilizing vaterite, we performed a molecular simulation of the substitution of sulfate for carbonate groups into the crystal structure of vaterite, aragonite and calcite. The results obtained show that the incorporation of small amounts (<3 mole%) of sulfate is energetically favorable in the vaterite structure, unfavorable in calcite and very unfavorable in aragonite. The computer modeling provided thermodynamic information, which, combined with kinetic arguments, allowed us to put forward a plausible explanation for the observed crystallization behavior.     

Schwertmannite transformation to goethite via the Fe(II) pathway: Reaction rates and implications for iron-sulfide formation

Schwertmannite (Fe₈O₈(OH)₆SO₄) is a common Fe(III)-oxyhydroxysulfate mineral in acid-sulfate systems, where its formation and fate strongly influence water quality. The present study examines transformation of schwertmannite to goethite (FeOOH), as catalyzed by interactions with Fe(II) in anoxic aquatic environments. This study also evaluates the role of the Fe(II) pathway in influencing the formation of iron-sulfide minerals in such environments. At pH > 5, the rates of Fe(II)-catalyzed schwertmannite transformation were several orders of magnitude faster than transformation in the absence of Fe(II). Complete transformation of schwertmannite occurred within only 3-5 h at pH > 6 and Fe(II)_(aq) ? 5 mmol L⁻¹. Model calculations indicate that the Fe(II)-catalyzed transformation of schwertmannite to goethite greatly decreases the reactivity of the Fe(III) pool, thereby favoring SO₄-reduction and facilitating the formation of iron-sulfide minerals (particularly mackinawite, tetragonal FeS). Examination of in situ sediment geochemistry in an acid-sulfate system revealed that the rapid Fe(II)-catalyzed transformation was consistent with an abrupt shift from an acidic Fe(III)-reducing regime with abundant schwertmannite near the sediment surface, to a near-neutral mackinawite-forming regime where goethite was dominant. This study demonstrates that the Fe(II) pathway exerts a major influence on schwertmannite transformation and iron-sulfide formation in anoxic acid-sulfate systems. These findings have important implications for understanding acidity dynamics and trace element mobility in such systems.     

An isopiestic study of aqueous NaBr and KBr at 50 °C: Chemical equilibrium model of solution behavior and solubility in the NaBr-H2O, KBr-H2O and Na-K-Br-H2O systems to high concentration and temperature

The isopiestic method has been used to determine the osmotic coefficients of the binary solutions NaBr-H₂O (from 0.745 to 5.953 mol kg⁻¹) and KBr-H₂O (from 0.741 to 5.683 mol kg⁻¹) at the temperature t = 50 °C. Sodium chloride solutions have been used as isopiestic reference standards. The isopiestic results obtained have been combined with all other experimental thermodynamic quantities available in literature (osmotic coefficients, water activities, bromide mineral’s solubilities) to construct a chemical model that calculates solute and solvent activities and solid-liquid equilibria in the NaBr-H₂O, KBr-H₂O and Na-K-Br-H₂O systems from dilute to high solution concentration within the 0-300 °C temperature range. The Harvie and Weare [Harvie C., and Weare J. (1980) The prediction of mineral solubilities in naturalwaters: the Na-K-Mg-Ca-Cl-SO₄-H₂O system from zero to high concentration at 25 °C. Geochim. Cosmochim. Acta44, 981-997] solubility modeling approach, incorporating their implementation of the concentration-dependent specific interaction equations of Pitzer [Pitzer K. (1973) Thermodynamics of electrolytes. I. Theoretical basis and general equations. J. Phys. Chem.77, 268-277] is employed. The model for binary systems is validated by comparing activity coefficient predictions with those given in literature, and not used in the parameterization process. Limitations of the mixed solutions model due to data insufficiencies are discussed. This model expands the variable temperature sodium-potassium model of Greenberg and Moller [Greenberg J., and Moller N. (1989) The prediction of mineral solubilities in natural waters: a chemical equilibrium model for the Na-K-Ca-Cl-SO₄-H₂O system to high concentration from 0 to 250 °C. Geochim. Cosmochim. Acta53, 2503-2518] by evaluating Br⁻ pure electrolyte and mixing solution parameters and the chemical potentials of three bromide solid phases: NaBr-2H₂O (cr), NaBr (cr) and KBr (cr).     

Carbon dioxide degassing and thermal energy release in the Monte Amiata volcanic-geothermal area (Italy)

The quaternary volcanic complex of Mount Amiata is located in southern Tuscany (Italy) and represents the most recent manifestation of the Tuscan Magmatic Province. The region is characterised by a large thermal anomaly and by the presence of numerous CO₂-rich gas emissions and geothermal features, mainly located at the periphery of the volcanic complex. Two geothermal systems are located, at increasing depths, in the carbonate and metamorphic formations beneath the volcanic complex. The shallow volcanic aquifer is separated from the deep geothermal systems by a low permeability unit (Ligurian Unit). A measured CO₂ discharge through soils of 1.8 × 10⁹ mol a⁻¹ shows that large amounts of CO₂ move from the deep reservoir to the surface. A large range in δ¹³C_TDIC (−21.07 to +3.65) characterises the waters circulating in the aquifers of the region and the mass and isotopic balance of TDIC allows distinguishing a discharge of 0.3 × 10⁹ mol a⁻¹ of deeply sourced CO₂ in spring waters. The total natural CO₂ discharge (2.1 × 10⁹ mol a⁻¹) is slightly less than minimum CO₂ output estimated by an indirect method (2.8 × 10⁹ mol a⁻¹), but present-day release of 5.8 × 10⁹ mol a⁻¹ CO₂ from deep geothermal wells may have reduced natural CO₂ discharge. The heat transported by groundwater, computed considering the increase in temperature from the infiltration area to the discharge from springs, is of the same order of magnitude, or higher, than the regional conductive heat flow (>200 mW m⁻²) and reaches extremely high values (up to 2700 mW m⁻²) in the north-eastern part of the study area. Heat transfer occurs mainly by conductive heating in the volcanic aquifer and by uprising gas and vapor along fault zones and in those areas where low permeability cover is lacking. The comparison of CO₂ flux, heat flow and geological setting shows that near surface geology and hydrogeological setting play a central role in determining CO₂ degassing and heat transfer patterns.     

Complexation of cadmium to sulfur and oxygen functional groups in an organic soil

Cadmium (Cd) is a toxic trace element and due to human activities soils and waters are contaminated by Cd both on a local and global scale. It is widely accepted that chemical interactions with functional groups of natural organic matter (NOM) is vital for the bioavailability and mobility of trace elements. In this study the binding strength of cadmium (Cd) to soil organic matter (SOM) was determined in an organic (49% organic C) soil as a function of reaction time, pH and Cd concentration. In experiments conducted at native Cd concentrations in soil (0.23 μg g⁻¹ dry soil), halides (Cl, Br) were used as competing ligands to functional groups in SOM. The concentration of Cd in the aqueous phase was determined by isotope-dilution (ID) inductively-coupled-plasma-mass-spectrometry (ICP-MS), and the activity of Cd²⁺ was calculated from the well-established Cd-halide constants. At higher Cd loading (500-54,000 μg g⁻¹), the Cd²⁺ activity was directly determined by an ion-selective electrode (ISE). On the basis of results from extended X-ray absorption fine structure (EXAFS) spectroscopy, a model with one thiolate group (RS⁻) was used to describe the complexation (Cd²⁺ + RS⁻ ? CdSR⁺; log K_CdSR) at native Cd concentrations. The concentration of thiols (RSH; 0.047 mol kg⁻¹ C) was independently determined by X-ray absorption near-edge structure (XANES) spectroscopy. Log K_CdSR values of 11.2-11.6 (pK_a for RSH = 9.96), determined in the pH range 3.1-4.6, compare favorably with stability constants for the association between Cd and well-defined thiolates like glutathione. In the concentration range 500-54,000 μg Cd g⁻¹, a model consisting of one thiolate and one carboxylate (RCOO⁻) gave the best fit to data, indicating an increasing role for RCOOH groups as RSH groups become saturated. The determined log K_CdOOCR of 3.2 (Cd²⁺ +  RCOO⁻ ? CdOOCR⁺; log K_CdOOCR; pK_a for RCOOH = 4.5) is in accordance with stability constants determined for the association between Cd and well-defined carboxylates. Given a concentration of reduced sulfur groups of 0.2% or higher in NOM, we conclude that the complexation to organic RSH groups may control the speciation of Cd in soils, and most likely also in surface waters, with a total concentration less than 5 mg Cd g⁻¹ organic C.     

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.     

The universal ratio of boron to chlorinity for the North Pacific and North Atlantic oceans

We report seawater boron concentration (mg kg⁻¹) and chlorinity (‰) values measured in seawater samples (n = 139) collected at various depths in the North Pacific and North Atlantic oceans and the East/Japan Sea (located in the western temperate North Pacific). Our results indicate that variations in seawater boron concentration are strongly coupled to variations in chlorinity (and salinity), yielding a mean boron to chlorinity ratio of 0.2414 ± 0.0009 mg kg⁻¹ ‰⁻¹ (boron to salinity ratio = 0.1336 ± 0.0005 mg kg⁻¹ ‰⁻¹). This ratio was surprisingly universal throughout the water column in the three marine basins and across widely different ocean surface regimes, but differs from the generally accepted ratio of 0.232 ± 0.005 mg kg⁻¹ ‰⁻¹ determined by Uppström (1974), which was based on only 20 measurements at four sites in the tropical Pacific Ocean. In converting total alkalinity to carbonate alkalinity (and vice versa) for thermodynamic calculations, the difference between these two ratios leads to a difference of 5 μmol kg⁻¹ in estimates for ocean surface waters, where the contribution of borate to total alkalinity is typically greatest. We suggest the use of the new boron to chlorinity ratio for predicting seawater boron concentrations using chlorinity (or salinity) data.     

Characterization and cycling of atmospheric mercury along the central US Gulf Coast

Concentrations of atmospheric Hg species, elemental Hg (Hg^°), reactive gaseous Hg (RGM), and fine particulate Hg (Hg-PM_2.5) were measured at a coastal site near Weeks Bay, Alabama from April to August, 2005 and January to May, 2006. Mean concentrations of the species were 1.6 ± 0.3 ng m⁻³, 4.0 ± 7.5 pg m⁻³ and 2.7 ± 3.4 pg m⁻³, respectively. A strong diel pattern was observed for RGM (midday maximum concentrations were up to 92.7 pg m⁻³), but not for Hg^° or Hg-PM_2.5. Elevated RGM concentrations (>25 pg m⁻³) in April and May of 2005 correlated with elevated average daytime O₃ concentrations (>55 ppbv) and high light intensity (>500 W m⁻²). These conditions generally corresponded with mixed continental-Gulf and exclusively continental air mass trajectories. Generally lower, but still elevated, RGM peaks observed in August, 2005 and January–March, 2006 correlated significantly (p < 0.05) with peaks in SO₂ concentration and corresponded to periods of high light intensity and lower average daytime O₃ concentrations. During these times air masses were dominated by trajectories that originated over the continent. Elevated RGM concentrations likely resulted from photochemical oxidation of Hg^° by atmospheric oxidants. This process may have been enhanced in and by the near-shore environment relative to inland sites. The marine boundary layer itself was not found to be a significant source of RGM.     

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.     

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.