Sulphate, Sulphide and the Sulphur Cycle

Sulphur is an essential nutrient for vegetation. It is also a major element of seawater and marine sediments. Its four main oxidation states ranging from +VI to –II make it both an electron acceptor and donor redox reactions. Major forms of sulphur in the subsurface include sulphate and sulphide minerals, dissolved sulphate (SO42–), dissolved sulphide (HS) and hydrogen sulphide gas (H2S). Organic sulphur is a component of organic compounds such as humic substances, kerogen and hydrocarbons. Its oxidation and recycling in soils produces the "terrestrial" sulphates found in semi-arid regions. Atmospheric sulphur sources include natural and "technogenic" or industrial SO2, particulate sulphur and aerosols of marine sulphate. The movement of sulphur through these various reservoirs in soils and the hydrosphere constitutes the sulphur cycle.

Sulphur compounds from these various sources participate in the geochemical evolution of groundwater. They also contribute to groundwater salinization. Some of the applications of sulphur isotope geochemistry include the cycling of sulphur in agricultural watersheds, the origin of salinity in coastal aquifers or sedimentary strata, groundwater contamination by landfill leachate plumes, acid mine drainage, and dating sulphate-reducing groundwaters.

Sulphur-34 is highly fractionated between sulphur compounds due to biological cycling. Similarly, the 18O content of sulphate is an important tool to trace the sulphur cycle. A summary of the ranges for 34S in natural materials is shown in Fig. 6-1. Meteorites and magmatic sulphur are close to the standard Cañon Diablo Troilite (CDT). Values exceeding +20‰ are found in association with evaporites and limestones. Negative d34S values are typical of diagenetic environments where reduced sulphur compounds are formed (Krouse, 1980). The most common reaction product is pyrite, which is present in many shales or other organic-rich sedimentary rocks and is formed by bacteria reducing seawater sulphate in marine sediments.

 

 

Fig. 6-1 Ranges in d34S contents of sulphur and sulphur compounds in different materials and environments (modified from Krouse 1980).

 

Marine sulphate

The sulphate concentration in groundwater can readily exceed 1 g/L due to the high solubility of gypsum. Gypsum [CaSO4·2H2O] and its unhydrated polymorph, anhydrite [CaSO4], are principal constituents of marine evaporites in sedimentary strata and can be major contributors to groundwater sulphate. Gypsum also accumulates in the soil of arid regions and in sabkha environments. In addition to dissolution of evaporites, marine sulphate in groundwaters comes from mixing with seawater in coastal aquifers and seawater-derived sulphate that accumulates in soils by evaporation. Sulphur isotopes and geochemistry can be used to distinguish between these sources. Other tools include ionic ratios, mineral saturation indices, carbonate isotope geochemistry, and 37Cl to trace chloride salinity. Table 6-1 gives the ionic composition of seawater, along with values for isotopic composition.
 

The solubility of gypsum in pure water is governed by the equilibrium equation:
 

CaSO4·2H2O ® Ca2+ + SO42– + 2H2O

and Kgyp = [Ca2+][SO42-] = 10–4.36

As equal moles of SO42– and Ca2+ are produced, the activity of SO42– in groundwater which has dissolved gypsum to saturation will be:
  [SO42–] = (10–4.36)½ = 10–2.1 = 0.0079 mole/L = 763 mg/L. This simple calculation overlooks the activity coefficient, gSO4, (mSO4 = aSO4/gSO4) which would be in the order of 0.7 (recall the influence of ionic strength I on g from Chapter 5), and so concentrations would exceed 1000 mg/L. In more saline groundwaters, gSO4 increases with ion complexation, allowing even higher concentrations in solution.
 

Table 6-1 The major ion geochemistry and isotopes of seawater (from Drever, 1988; Longinelli, 1989; Tan, 1989; Bassett, 1990; Chan, 1992)

 

         
        Species mg/L mmole/L Isotopic composition
        Cl 19,350  546 d37Cl = 0‰ SMOC
        SO42– 2710  28.3 d34S = 21‰ CDT
        " " " d18OSO4 = 9.5‰ VSMOW
        HCO3 142  2.33 d13C = –1 to +2‰ VPDB
        Br 67 0.84  
        F 1.3  0.068  
        NO3 0.005–2  0.8-300 · 10–4  
        PO43– 0.001–0.05  1-50 · 10–5 d18OPO4 = 19.7‰ VSMOW
        DOC 0.3–2 as C 0.02-0.2 d13C ~ –25 to –30‰ VPDB
        Na+ 10,760 468  
        Mg2+ 1290  53.1  
        Ca2+ 411  10.3  
        K+ 399  10.2  
        Sr2+ 8 0.091 87Sr/86Sr = 0.70924 ± 0.00003
        SiO2 0.5–10  < 0.35  
        B3+ 4.5  0.42 d11B = 40.4 ± 1.7‰ NBS 951
        Li+ 0.18  0.026 d6Li = –32.3 ± 0.5‰ LSVEC
         
Gypsum and anhydrite are not restricted to major evaporite sequences in the geological record. Minor sulphate nodules and laminations on bedding planes can occur throughout limestone and dolomite sequences, and can be overlooked during geological mapping. The dissolution of sulphate minerals increases permeability, and so can enhance concentrations in groundwater.

In pure water, the Ca2+/SO42– molar ratio equal to 1 would distinguish gypsum dissolution from other sources of sulphate salinity such as seawater (0.36) (Table 6-2). Calcite precipitation also affects this ratio. Groundwaters typically move towards calcite saturation within the soil or along the flow path. Dissolution of sulphate then forces calcite precipitation according to the common ion effect, resulting in a disproportionate increase in SO42–:

Ca2+ + HCO3 + CaSO4·2H2O ® Ca2+ + SO42– + CaCO3 + H+ + 2H2O

The dissolution of gypsum or anhydrite occurs without measurable isotope effects, and so the isotope contents of SO42– can be used as a tracer for sulphate origin. The sulphate of modern seawater has a very homogeneous and well-defined isotopic composition (Table 6-2):
 

d34SSO4 = 21‰ CDT

d18OSO4 = 9.5‰ VSMOW

 

 

Fig. 6-2 The 34S and 18O composition of marine sulphate through geologic time (modified from Claypool et al., 1980, with data from Fritz et al., 1988).

 
This composition was not constant in the past. Pronounced variations exist throughout geologic time (Claypool et al., 1980; Fig. 6-2). These variations are found in all major marine evaporite deposits and were most likely controlled by major inputs or removal of sulphide from the oceanic reservoirs during changes in tectonic activity and weathering rates. Simple removal of sulphate (increase in evaporite formation) would not be accompanied by such dramatic isotope effects. The d18O content of seawater sulphate has been more stable over geologic time, and is controlled largely by the sulphide weathering reactions that contribute sulphate to seawater.

Comparison of the isotope contents of sulphate in groundwater with the appropriate geological period in Fig. 6-2 can distinguish a geogenic source from other sources of sulphate salinity.