The origin and isotopic composition of dissolved sulfide in groundwater from carbonate aquifers in Florida and Texas

The δ³⁴S values of dissolved sulfide and the sulfur isotope fractionations between dissolved sulfide and sulfate species in Floridan ground water generally correlate with dissolved sulfate concentrations which are related to flow patterns and residence time within the aquifer. The dissolved sulfide derives from the slow in situ biogenic reduction of sulfate dissolved from sedimentary gypsum in the aquifer. In areas where the water is oldest, the dissolved sulfide has apparently attained isotopic equilibrium with the dissolved sulfate (Δ³⁴S = 65 per mil) at the temperature (28°C) of the system. This approach to equilibrium reflects an extremely slow reduction rate of the dissolved sulfate by bacteria; this slow rate probably results from very low concentrations of organic matter in the aquifer.In the reducing part of the Edwards aquifer, Texas, there is a general down-gradient increase in both dissolved sulfide and sulfate concentrations, but neither the δ³⁴S values of sulfide nor the sulfide-sulfate isotope fractionation correlates with the ground-water flow pattern. The dissolved sulfide species appear to be derived primarily from biogenic reduction of sulfate ions whose source is gypsum dissolution although upgradient diffusion of H₂S gas from deeper oil field brines may be important in places. The sulfur isotope fractionation for sulfide-sulfate (about 38 per mil) is similar to that observed for modern oceanic sediments and probably reflects moderate sulfate reduction in the reducing part of the aquifer owing to the higher temperature and significant amount of organic matter present; contributions of isotopically heavy H₂S from oil field brines are also possible.     

Characterization of chromophoric dissolved organic matter (CDOM) in the Baltic Sea by excitation emission matrix fluorescence spectroscopy

Chromophoric dissolved organic matter (CDOM) is the major light absorber in the Baltic Sea. In this study, excitation emission matrix (EEM) fluorescence spectra and UV–visible absorption spectra of CDOM are reported as a function of salinity. Samples from different locations and over different seasons were collected during four cruises in 2002 and 2003 in the Baltic Sea in both Pomeranian Bay and the Gulf of Gdansk. Absorption by CDOM decreased with increased distance from the riverine source and reached a relatively stable absorption background in the open sea. Regression analysis showed that fluorescence intensity was linearly related to absorption by CDOM at 375 nm and a_CDOM(375) absorption coefficients were inversely related to salinity. Analysis of CDOM-EEM spectra indicated that a change in composition of CDOM occurred along the salinity gradient in the Baltic Sea. Analysis of percent contribution of respective fluorophore groups to the total intensity of EEM spectra indicated that the fluorescence peaks associated with terrestrial humic components of the CDOM and total integrated fluorescence decreased with decreasing CDOM absorption. In contrast, the protein-like fraction of CDOM decreased to a lesser degree than the others. Analysis of the percent contribution of fluorescence peak intensities to the total fluorescence along the salinity gradient showed that the contribution of protein-like fluorophores increased from 2.6% to 5.1% in the high-salinity region of the transect. Fluorescence and absorption changes observed in the Baltic Sea were similar to those observed in similar transects that have been sampled elsewhere, e.g. in European estuaries, Gulf of Mexico, Mid-Atlantic Bight and the Cape Fear River plume in the South Atlantic Bight, although the changes in the Baltic Sea occurred over a much smaller salinity gradient.     

A method of alternating characteristics with application to advection-dominated environmental systems

We present a numerical method for solving a class of systems of partial differential equations (PDEs) that arises in modeling environmental processes undergoing advection and biogeochemical reactions. The salient feature of these PDEs is that all partial derivatives appear in linear expressions. As a result, the system can be viewed as a set of ordinary differential equations (ODEs), albeit each one along a different characteristic. The method then consists of alternating between equations and integrating each one step-wise along its own characteristic, thus creating a customized grid on which solutions are computed. Since the solutions of such PDEs are generally smoother along their characteristics, the method offers the potential of using larger time steps while maintaining accuracy and reducing numerical dispersion. The advantages in efficiency and accuracy of the proposed method are demonstrated in two illustrative examples that simulate depth-resolved reactive transport and soil carbon cycling.