Trace element cycling in a subterranean estuary: Part 2. Geochemistry of the pore water

Submarine groundwater discharge (SGD) is an important source of dissolved elements to the ocean, yet little is known regarding the chemical reactions that control their flux from sandy coastal aquifers. The net flux of elements from SGD to the coastal ocean is dependent on biogeochemical reactions in the groundwater-seawater mixing zone, recently termed the “subterranean estuary.” This paper is the second in a two part series on the biogeochemistry of the Waquoit Bay coastal aquifer/subterranean estuary. The first paper addressed the biogeochemistry of Fe, Mn, P, Ba, U, and Th from the perspective of the sediment composition of cores Charette et al. [Charette, M.A., Sholkovitz, E.R., Hansell, C.M., 2005. Trace element cycling in a subterranean estuary: Part 1. Geochemistry of the permeable sediments. Geochim. Cosmochim. Acta, 69, 2095-2109]. This paper uses pore water data from the subterranean estuary, along with Bay surface water data, to establish a more detailed view into the estuarine chemistry and the chemical diagenesis of Fe, Mn, U, Ba and Sr in coastal aquifers. Nine high-resolution pore water (groundwater) profiles were collected from the head of the Bay during July 2002. There were non-conservative additions of both Ba and Sr in the salinity transition zone of the subterranean estuary. However, the extent of Sr release was significantly less than that of its alkaline earth neighbor Ba. Pore water Ba concentrations approached 3000 nM compared with 25-50 nM in the surface waters of the Bay; the pore water Sr-salinity distribution suggests a 26% elevation in the amount of Sr added to the subterranean estuary. The release of dissolved Ba to the mixing zone of surface estuaries is frequently attributed to an ion-exchange process whereby seawater cations react with Ba from river suspended clay mineral particles at low to intermediate salinity. Results presented here suggest that reductive dissolution of Mn oxides, in conjunction with changes in salinity, may also be an important process in maintaining high concentrations of Ba in the pore water of subterranean estuaries. In contrast, pore water U was significantly depleted in the subterranean estuary, a result of SGD-driven circulation of seawater through reducing permeable sediments. This finding is supported by surface water concentrations of U in the Bay, which were significantly depleted in U compared with adjacent coastal waters. Using a global estimate of SGD, we calculate U removal in subterranean estuaries at 20 × 10⁶ mol U y⁻¹, which is the same order of magnitude as the other major U sinks for the ocean. Our results suggest a need to revisit and reevaluate the oceanic budgets for elements that are likely influenced by SGD-associated processes.     

Iron isotope systematics in estuaries: The case of North River, Massachusetts (USA)

Recent studies have suggested that rivers may present an isotopically light Fe source to the oceans. Since the input of dissolved iron from river water is generally controlled by flocculation processes that occur during estuarine mixing, it is important to investigate potential fractionation of Fe-isotopes during this process. In this study, we investigate the influence of the flocculation of Fe-rich colloids on the iron isotope composition of pristine estuarine waters and suspended particles. The samples were collected along a salinity gradient from the fresh water to the ocean in the North River estuary (MA, USA). Estuarine samples were filtered at 0.22 μm and the iron isotope composition of the two fractions (dissolved and particles) were analyzed using high-resolution MC-ICP-MS after chemical purification. Dissolved iron results show positive δ⁵⁶Fe values (with an average of 0.43 ± 0.04‰) relative to the IRMM-14 standard and do not display any relationships with salinity or with percentage of colloid flocculation. The iron isotopic composition of the particles suspended in fresh water is characterized by more negative δ⁵⁶Fe values than for dissolved Fe and correlate with the percentage of Fe flocculation. Particulate δ⁵⁶Fe values vary from −0.09‰ at no flocculation to ∼0.1‰ at the flocculation maximum, which reflect mixing effects between river-borne particles, lithogenic particles derived from coastal seawaters and newly precipitated colloids. Since the process of flocculation produces minimal Fe-isotope fractionation in the dissolved Fe pool, we suggest that the pristine iron isotope composition of fresh water is preserved during estuarine mixing and that the value of the global riverine source into the ocean can be identified from the fresh water values. However, this study also suggests that δ⁵⁶Fe composition of rivers can also be characterized by more positive δ⁵⁶Fe values (up to 0.3‰) relative to the crust than previously reported. In order to improve our current understanding of the oceanic iron isotope cycling, further work is now required to determine the processes controlling the fractionation of Fe-isotopes during continental run-off.     

Application of iron and zinc isotopes to track the sources and mechanisms of metal loading in a mountain watershed

Here the hydrogeochemical constraints of a tracer dilution study are combined with Fe and Zn isotopic measurements to pinpoint metal loading sources and attenuation mechanisms in an alpine watershed impacted by acid mine drainage. In the tested mountain catchment, δ⁵⁶Fe and δ⁶⁶Zn isotopic signatures of filtered stream water samples varied by ∼3.5‰ and 0.4‰, respectively. The inherent differences in the aqueous geochemistry of Fe and Zn provided complimentary isotopic information. For example, variations in δ⁵⁶Fe were linked to redox and precipitation reactions occurring in the stream, while changes in δ⁶⁶Zn were indicative of conservative mixing of different Zn sources. Fen environments contributed distinctively light dissolved Fe (<−2.0‰) and isotopically heavy suspended Fe precipitates to the watershed, while Zn from the fen was isotopically heavy (>+0.4‰). Acidic drainage from mine wastes contributed heavier dissolved Fe (∼+0.5‰) and lighter Zn (∼+0.2‰) isotopes relative to the fen. Upwelling of Fe-rich groundwater near the mouth of the catchment was the major source of Fe (δ⁵⁶Fe ∼ 0‰) leaving the watershed in surface flow, while runoff from mining wastes was the major source of Zn. The results suggest that given a strong framework for interpretation, Fe and Zn isotopes are useful tools for identifying and tracking metal sources and attenuation mechanisms in mountain watersheds.     

Paleo-redox boundaries in fractured granite

At the Earth’s surface, Fe(II) often oxidises and forms insoluble Fe(III)-(oxyhydr)oxides, whose particle size and structure depend on solution composition and temperature during formation and afterwards. Bacterial processes and exposure to reducing environments reduces them again, releasing dissolved iron to the groundwater. During such cycling, the Fe isotopes fractionate to an extent that is expected to depend on temperature. In this study, we report on the use of Fe-oxides as paleo-redox indicators, using their structure, morphology and Fe-composition as a clue for formation conditions. In samples taken from ∼120 m drill cores in granite from SE Sweden, X-ray amorphous, superparamagnetic, nanometre-sized Fe-oxides are confined to fractures of the upper ∼50 m, whereas well-crystalline Fe-oxides, with particle sizes typical for soils, occur down to ∼110 m. We also identified hematite with a particle size of 100 nm, similar to hematite of hydrothermal origin. The Fe isotope composition of the fine-grained Fe-oxides (−1‰ < δ⁵⁶Fe < 1‰, IRMM-14 referenced) scatter significantly compared to the distribution previously observed for hydrothermal material (−0.26‰ < δ⁵⁶Fe < 0.12‰) and they are dominantly heavier than Fe-bearing silicates from fractures (−0.56‰ < δ⁵⁶Fe < −0.35‰). This is consistent with formation by low-temperature weathering, where Fe-silicates dissolve, Fe(II) oxidises and Fe(III)-oxides precipitate. The X-ray amorphous, nanometre-sized nature of near-surface Fe-oxides suggests recent formation. The deeper situated, well-crystalline Fe-oxides are more mature and we interpret that they record earlier oxidising events. They exist in fractures that are not significantly altered, indicating formation during periods of oxidation. Our results show that oxygenated water may reach depths of ∼110 m in fractured granite. The absence of natural, low-temperature Fe-oxides from deeper drill cores suggests that oxygenated waters do not readily penetrate beyond about 100 m and suggests that radioactive waste repositories located at a depth of ∼500 m should be well-protected from oxygenated waters.     

The effect of early diagenesis on the Fe isotope compositions of porewaters and authigenic minerals in continental margin sediments

Iron isotope compositions in marine pore fluids and sedimentary solid phases were measured at two sites along the California continental margin, where isotope compositions range from δ⁵⁶Fe = −3.0‰ to +0.4‰. At one site near Monterey Canyon off central California, organic matter oxidation likely proceeds through a number of diagenetic pathways that include significant dissimilatory iron reduction (DIR) and bacterial sulfate reduction, whereas at our other site in the Santa Barbara basin DIR appears to be comparatively small, and production of sulfides (FeS and pyrite) was extensive. The largest range in Fe isotope compositions is observed for Fe(II)_aq in porewaters, which generally have the lowest δ⁵⁶Fe values (minimum: −3.0‰) near the sediment surface, and increase with burial depth. δ⁵⁶Fe values for FeS inferred from HCl extractions vary between ∼−0.4‰ and +0.4‰, but pyrite is similar at both stations, where an average δ⁵⁶Fe value of −0.8 ± 0.2‰ was measured. We interpret variations in dissolved Fe isotope compositions to be best explained by open-system behavior that involves extensive recycling of Feflux. This study is the first to examine Fe isotope variations in modern marine sediments, and the results show that Fe isotopes in the various reactive Fe pools undergo isotopic fractionation during early diagenesis. Importantly, processes dominated by sulfide formation produce high-δ⁵⁶Fe values for porewaters, whereas the opposite occurs when Fe(III)-oxides are present and DIR is a major pathway of organic carbon respiration. Because shelf pore fluids may carry a negative δ⁵⁶Fe signature it is possible that the Fe isotope composition of ocean water reflects a significant contribution of shelf-derived iron to the open ocean. Such a signature would be an important means for tracing iron sources to the ocean and water mass circulation.     

In situ iron isotope ratio determination using UV-femtosecond laser ablation with application to hydrothermal ore formation processes

The feasibility of in situ stable Fe isotope ratio measurements using UV-femtosecond laser ablation connected to a multiple-collector inductively coupled plasma mass spectrometer (MC-ICP-MS) has been investigated. Different types of matrices, independently determined by solution MC-ICP-MS after chromatographic separation of Fe, have been analysed by laser ablation using the isotopically certified iron reference material IRMM-014 as the bracketing standard. The samples have been pure iron metal (JM Puratronic), Fe-meteorites (North Chile, Glenormiston and Toluca), the meteorite phases kamacite and taenite in Toluca and Fe-sulphides. Furthermore, Fe isotope ratios from hydrothermal hematite, siderite and goethite from an old mining area in the Schwarzwald, Germany, and of magnetite from the metamorphic Biwabik iron formation have been determined. The results show that a precision of better than 0.1‰ (2 sigma) can be achieved with laser ablation and that all the results obtained agree with those determined by solution ICP to better than 0.1‰. This precision and accuracy is achievable in both raster and spot ablation mode. A matrix-matched bracketing standard is not required , and all these materials can be measured accurately against a metal standard. The hydrothermal minerals show significant Fe isotope zonations. In some samples the range of δ⁵⁶Fe in a single aggregate encompasses the entire spectrum of ratios found by bulk solution analyses in multiple samples distributed over the whole mining district. For example, isotopic zonations found in secondary fibrous hematites show a continuous change in the δ⁵⁶Fe values from −0.5‰ in the core to −1.8‰ in the rim. Primary hydrothermal siderite shows the reverse pattern with lighter values in the core than in the rim. While the siderite is thought to record primary fluid histories, the hematite pattern is interpreted as a reworked isotopic signature generated by oxic dissolution of primary zoned siderite and immediate close range re-precipitation of the oxidized Fe. Abrupt changes are documented for secondary goethite showing a distinct overgrowth that is 0.4‰ lighter than the core of the grain. If indeed Fe isotopes in secondary minerals from hydrothermal ore deposits record the initial isotopic signatures of their precursor minerals, and these in turn record hydrothermal fluid histories, then the tools are in place for a detailed reconstruction of the deposit‘s genesis. We expect similar observations from other Fe-rich deposits formed at intermediate and low-temperatures (e.g. banded iron formations). Laser ablation now provides us with the spatial resolution that adds a further dimension to our interpretation of stable Fe-isotope fractionation.     

Submarine groundwater discharge is an important net source of light and middle REEs to coastal waters of the Indian River Lagoon, Florida, USA

Porewater (i.e., groundwater) samples were collected from multi-level piezometers across the freshwater-saltwater seepage face within the Indian River Lagoon subterranean estuary along Florida’s (USA) Atlantic coast for analysis of the rare earth elements (REE). Surface water samples for REE analysis were also collected from the water column of the Indian River Lagoon as well as two local rivers (Eau Gallie River, Crane Creek) that flow into the lagoon within the study area. Concentrations of REEs in porewaters from the subterranean estuary are 10-100 times higher than typical seawater values (e.g., Nd ranges from 217 to 2409 pmol kg⁻¹), with submarine groundwater discharge (SGD) at the freshwater-saltwater seepage face exhibiting the highest REE concentrations. The elevated REE concentrations for SGD at the seepage face are too high to be the result of simple, binary mixing between a seawater end-member and local terrestrial SGD. Instead, the high REE concentrations indicate that geochemical reactions occurring within the subterranean estuary contribute substantially to the REE cycle. A simple mass balance model is used to investigate the cycling of REEs in the Indian River Lagoon and its underlying subterranean estuary. Mass balance modeling reveals that the Indian River Lagoon is approximately at steady-state with respect to the REE fluxes into and out of the lagoon. However, the subterranean estuary is not at steady-state with respect to the REE fluxes. Specifically, the model suggests that the SGD Nd flux, for example, exported from the subterranean estuary to the overlying lagoon waters exceeds the combined input to the subterranean estuary from terrestrial SGD and recirculating marine SGD by, on average, ∼100 mmol day⁻¹. The mass balance model also reveals that the subterranean estuary is a net source of light REEs (LREE) and middle REEs (MREE) to the overlying lagoon waters, but acts as a sink for the heavy REEs (HREE). Geochemical modeling and statistical analysis further suggests that this fractionation occurs, in part, due to the coupling between REE cycling and iron redox cycling within the Indian River Lagoon subterranean estuary. The net SGD flux of Nd to the Indian River Lagoon is ∼7-fold larger than the local effective river flux to these coastal waters. This previously unrecognized source of Nd to the coastal ocean could conceivably be important to the global oceanic Nd budget, and help to resolve the oceanic “Nd paradox” by accounting for a substantial fraction of the hypothesized missing Nd flux to the ocean.     

Trace element cycling in a subterranean estuary: Part 1. Geochemistry of the permeable sediments

Subterranean estuaries are characterized by the mixing of terrestrially derived groundwater and seawater in a coastal aquifer. Subterranean estuaries, like their river water-seawater counterparts on the surface of the earth, represent a major, but less visible, hydrological and geochemical interface between the continents and the ocean. This article is the first in a two-part series on the biogeochemistry of the subterranean estuary at the head of Waquoit Bay (Cape Cod, MA, USA). The pore-water distributions of salinity, Fe and Mn establish the salt and redox framework of this subterranean estuary. The biogeochemistry of Fe, Mn, P, Ba, U and Th will be addressed from the perspective of the sediment composition. A second article will focus on the groundwater and pore-water chemistries of Fe, Mn, U and Ba.Three sediment cores were collected from the head of Waquoit Bay where the coastal aquifer consists of permeable sandy sediment. A selective dissolution method was used to measure the concentrations of P, Ba, U and Th that are associated with “amorphous (hydr)oxides of iron and manganese” and “crystalline Fe and Mn (hydr)oxides.” The deeper sections of the cores are characterized by large amounts of iron (hydr)oxides that are precipitated onto organic C-poor quartz sand from high-salinity pore waters rich in dissolved ferrous iron. Unlike Fe (hydr)oxides, which increase with depth, the Mn (hydr)oxides display midcore maxima. This type of vertical stratification is consistent with redox-controlled diagenesis in which Mn (hydr)oxides are formed at shallower depths than iron (hydr)oxides. P and Th are enriched in the deep sections of the cores, consistent with their well-documented affinity for Fe (hydr)oxides. In contrast, the downcore distribution of Ba, especially in core 3, more closely tracks the concentration of Mn (hydr)oxides. Even though Mn (hydr)oxides are 200-300 times less abundant than Fe (hydr)oxides in the cores, Mn (hydr)oxides are known to have an affinity for Ba which is many orders of magnitude greater than iron (hydr)oxides. Hence, the downcore distribution of Ba in Fe (hydr)oxide rich sediments is most probably controlled by the presence of Mn (hydr)oxides. U is enriched in the upper zones of the cores, consistent with the formation of highly reducing near-surface sediments in the intertidal zone at the head of the Bay. Hence, the recirculation of seawater through this type of subterranean estuary, coupled with the abiotic and/or biotic reduction of soluble U(VI) to insoluble U(IV), leads to the sediments acting as a oceanic net sink of U. These results highlight the importance of permeable sediments as hosts to a wide range of biogeochemical reactions, which may be impacting geochemical budgets on scales ranging from coastal aquifers to the continental shelf.     

Iron isotope fractionation during hydrothermal ore deposition and alteration

Iron isotopes fractionate during hydrothermal processes. Therefore, the Fe isotope composition of ore-forming minerals characterizes either iron sources or fluid histories. The former potentially serves to distinguish between sedimentary, magmatic or metamorphic iron sources, and the latter allows the reconstruction of precipitation and redox processes. These processes take place during ore formation or alteration. The aim of this contribution is to investigate the suitability of this new isotope method as a probe of ore-related processes. For this purpose 51 samples of iron ores and iron mineral separates from the Schwarzwald region, southwest Germany, were analyzed for their iron isotope composition using multicollector ICP-MS. Further, the ore-forming and ore-altering processes were quantitatively modeled using reaction path calculations. The Schwarzwald mining district hosts mineralizations that formed discontinuously over almost 300 Ma of hydrothermal activity. Primary hematite, siderite and sulfides formed from mixing of meteoric fluids with deeper crustal brines. Later, these minerals were partly dissolved and oxidized, and secondary hematite, goethite and iron arsenates were precipitated. Two types of alteration products formed: (1) primary and high-temperature secondary Fe minerals formed between 120 and 300 °C, and (2) low-temperature secondary Fe minerals formed under supergene conditions (<100 °C). Measured iron isotope compositions are variable and cover a range in δ⁵⁶Fe between −2.3‰ and +1.3‰. Primary hematite (δ⁵⁶Fe: −0.5‰ to +0.5‰) precipitated by mixing oxidizing surface waters with a hydrothermal fluid that contained moderately light Fe (δ⁵⁶Fe: −0.5‰) leached from the crystalline basement. Occasional input of CO₂-rich waters resulted in precipitation of isotopically light siderite (δ⁵⁶Fe: −1.4 to −0.7‰). The difference between hematite and siderite is compatible with published Fe isotope fractionation factors. The observed range in isotopic compositions can be accounted for by variable fractions of Fe precipitating from the fluid. Therefore, both fluid processes and mass balance can be inferred from Fe isotopes. Supergene weathering of siderite by oxidizing surface waters led to replacement of isotopically light primary siderite by similarly light secondary hematite and goethite, respectively. Because this replacement entails quantitative transfer of iron from precursor mineral to product, no significant isotope fractionation is produced. Hence, Fe isotopes potentially serve to identify precursors in ore alteration products. Goethites from oolitic sedimentary iron ores were also analyzed. Their compositional range appears to indicate oxidative precipitation from relatively uniform Fe dissolved in coastal water. This comprehensive iron isotope study illustrates the potential of the new technique in deciphering ore formation and alteration processes. Isotope ratios are strongly dependent on and highly characteristic of fluid and precipitation histories. Therefore, they are less suitable to provide information on Fe sources. However, it will be possible to unravel the physico-chemical processes leading to the formation, dissolution and redeposition of ores in great detail.     

Hydrogen isotopes in individual alkenones from the Chesapeake Bay estuary

Hydrogen isotope ratios of individual alkenones from haptophyte algae were measured in suspended particles and surface sediment from the Chesapeake Bay (CB) estuary, eastern USA, in order to determine their relationship to water δD values and salinity. δD values of four alkenones (MeC_37:2, MeC_37:3, EtC_38:2, EtC_38:3) from particles and sediments were between −165‰ and −221‰ and increased linearly (R² = 0.7-0.9) with water δD values from the head to the mouth of the Bay. Individual alkenones were depleted in deuterium by 156-188‰ relative to water. The MeC₃₇ alkenones were consistently enriched by ∼12‰ relative to the EtC₃₈ alkenones, and the di-unsaturated alkenones of both varieties were consistently enriched by ∼20‰ relative to the tri-unsaturated alkenones. All of the increase in alkenone δD values could be accounted for by the water δD increase. Consequently, no net change in alkenone-water D/H fractionation occurred as a result of the salinity increase from 10 to 29. This observation is at odds with results from culture studies with alkenone-producing marine coccolithophorids, and from two field studies, one with a dinoflagellate sterol in the CB, and one with a wide variety of lipids in saline ponds on Christmas Island, that indicate a decline in D/H fractionation with increasing salinity. Why D/H fractionation in alkenones in the CB showed no dependence on salinity, while D/H fractionation in CB dinsoterol decreased by 1‰ per unit increase in salinity remains to be determined. Two hypotheses we consider to be valid are that (i) the assemblage of alkenone-producing haptophytes changes along the Bay and each species has a different sensitivity to salinity, such that no apparent trend in α_{alkenone-water} occurs along the salinity gradient, and (ii) greater osmoregulation capacity in coastal haptophytes may result in a diminished sensitivity of alkenone-water D/H fractionation to salinity changes.     

Carbon, oxygen, and strontium isotope geochemistry of carbonate rocks of the upper Miocene Kudankulam Formation, southern India: Implications for paleoenvironment and diagenesis

The carbon, oxygen, and strontium isotope compositions of carbonate rocks from the upper Miocene Kudankulam Formation, southern India, were measured to understand palaeoenvironment and carbonate diagenesis of this formation. Both carbon and oxygen isotope ratios of various carbonate phases including whole rocks, ooids, molluscan mold-fill and sparry pore-fill calcite cements are depleted in ¹⁸O and ¹³C compared to those of contemporaneous seawater, indicating that the Kudankulam carbonates underwent extensive meteoric diagenesis. Based on δ¹³C and δ¹⁸O values for sparry calcite cements (pore-fill and molluscan mold-fill) formed in the meteoric diagenetic realm (δ¹³C from −7.8‰ to −6.0‰ and −9.0‰ to −7.0‰; δ¹⁸O from −9.2‰ to −6.5‰ and −9.4‰ to −2.6‰, respectively), it is interpreted that the diagenetic system was open and was proximal to the vadose water recharge zone. The negative δ¹⁸O values of various carbonate components (about −9.4‰ to −4.1‰ for whole rocks; about −8.4‰ to −2.6‰ for ooids) suggest that during the late Miocene the paleoclimate of the study area was humid, unlike today, probably due to the intense Indian monsoon system. The carbon isotope compositions (−7.9‰ to −3.6‰ for whole rocks; −4.9‰ to −1.5‰ for ooids) are consistent with the interpretation that the paleo-ecosystem comprised a significant proportion of C₄ type plants, supporting a scenario of expansion of C₄ plants during the late Miocene in the Indian subcontinent as far south as the southern tip of India. The ⁸⁷Sr/⁸⁶Sr ratios of the Kudankulam carbonates (0.70920 to 0.72130) are much greater than those of the contemporaneous or modern seawater (between 0.7089 and 0.7091) and show a general decrease up-sequence. Such high Sr isotope ratios indicate significant radiogenic ⁸⁷Sr influx to the system from the Archean rocks exposed in the drainage area, implying that the deep-seated Archean rocks were already exposed in southern India by the late Miocene.     

pH-Dependent iron oxide precipitation in a subterranean estuary

Iron-oxide-coated sediment particles in subterranean estuaries can act as a geochemical barrier (“iron curtain”) for various chemical species in groundwater (e.g. phosphate), thus limiting their discharge to coastal waters. Little is known about the factors controlling this Fe-oxide precipitation. Here, we implement a simple reaction network in a 1D reactive transport model (RTM), to investigate the effect of O₂ and pH gradients along a flow-line in the subterranean estuary of Waquoit Bay (Cape Cod, Massachusetts) on oxidative precipitation of Fe(II) and subsequent PO₄ sorption. Results show that the observed O₂ gradient is not the main factor controlling precipitation and that it is the pH gradient at the mixing zone of freshwater (pH 5.5) and seawater (pH 7.9) near the beach face that causes a  7-fold increase in the rate of oxidative precipitation of Fe(II) at  15 m. Thus, the pH gradient determines the location and magnitude of the observed iron oxide accumulation and the subsequent removal of PO₄ in this subterranean estuary.     

Deciphering interaction of regional aquifers in Southern Tunisia using hydrochemistry and isotopic tools

Groundwater is the most important source of water supply in southern Tunisia. Previous hydrogeologic and isotopic studies carried out in this region revealed the existence of two major aquifer systems: the “Complex Terminal” (CT) and the “Continental Intercalaire” (CI). Turonian carbonates constitute one of the major aquifer levels of the CT multilayered aquifer. It extends over most of southern Tunisia, and its hydrodynamic regime is largely influenced by tectonics, lithology and recharge conditions. Forty-eight groundwater samples from the CI and Turonian aquifers were collected between January and April 2004 for chemical and isotopic analyses. Hydrochemistry and isotopic tools were combined to get an insight into the processes controlling chemical composition of groundwater and wide-scale interaction of these two aquifer systems. Analysis of the dissolved constituents revealed that several processes control the observed chemical composition: (i) incongruent dissolution of carbonate minerals, (ii) dissolution of evaporitic minerals, and (iii) cation exchange. Dissolution alone cannot account for the observed high supersaturation states of groundwater with respect to calcite and dolomite. The observed supersaturation is most probably linked to geogenic CO₂ entering water-bearing horizons of the CT and CI aquifers via deep tectonic faults and discontinuities and subsequent degassing in the exploitation wells. Presence of geogenic CO₂ in the investigated region was confirmed by C isotope data of the DIC reservoir. The radiocarbon content of the Turonian samples varied between 9.5 and 43 pmc. For CI samples generally lower values were recorded, between 3.8 and 22.5 pmc. Stable isotope composition of Turonian groundwater samples varied from −8.3 to −5.3‰ for δ¹⁸O and from −60 to −25‰ for δ²H. The corresponding ranges of δ values for the Continental Intercalaire samples were from −8.9‰ to −6.9‰ for δ¹⁸O and from −68.2‰ to −45.7‰ for δ²H. Stable isotope composition of groundwater representing CT and CI aquifers provide strong evidence for regional interaction between both systems.     

Stable silicon isotopes of groundwater, feldspars, and clay coatings in the Navajo Sandstone aquifer, Black Mesa, Arizona, USA

We present some of the first analyses of the stable isotopic composition of dissolved silicon (Si) in groundwater. The groundwater samples were from the Navajo Sandstone aquifer at Black Mesa, Arizona, USA, and the Si isotope composition of detrital feldspars and secondary clay coatings in the aquifer were also analyzed. Silicon isotope compositions were measured using high-resolution multi-collector inductively coupled mass spectrometry (HR-MC-ICP-MS) (Nu1700 & NuPlasma HR). The quartz dominated bulk rock and feldspar separates have similar δ³⁰Si of −0.09 ± 0.04‰ and −0.15 ± 0.04‰ (±95% SEM), respectively, and clay separates are isotopically lighter by up to 0.4‰ compared to the feldspars. From isotopic mass-balance considerations, co-existing aqueous fluids should have δ³⁰Si values heavier than the primary silicates. Positive δ³⁰Si values were only found in the shallow aquifer, where Si isotopes are most likely fractionated during the dissolution of feldspars and subsequent formation of clay minerals. However, δ³⁰Si decreases along the flow path from 0.56‰ to −1.42‰, representing the most negative dissolved Si isotope composition so far found for natural waters. We speculate that the enrichment in ²⁸Si is due to dissolution of partly secondary clay minerals and low-temperature silcretes in the Navajo Sandstone. The discovery of the large range and systematic shifts of δ³⁰Si values along a groundwater flow path illustrates the potential utility of stable Si isotopes for deciphering the Si cycling in sedimentary basins, tracing fluid flow, and evaluating global Si cycle.     

Subsurface processes at the lucky strike hydrothermal field, Mid-Atlantic ridge: evidence from sulfur, selenium, and iron isotopes

At Lucky Strike near the Azores Triple Junction, the seafloor setting of the hydrothermal field in a caldera system with abundant low-permeability layers of cemented breccia, provides a unique opportunity to study the influence of subsurface geological conditions on the hydrothermal fluid evolution. Coupled analyses of S isotopes performed in conjunction with Se and Fe isotopes have been applied for the first time to the study of seafloor hydrothermal systems. These data provide a tool for resolving the different abiotic and potential biotic near-surface hydrothermal reactions. The δ³⁴S (between 1.5‰ and 4.6‰) and Se values (between 213 and 1640 ppm) of chalcopyrite suggest a high temperature end-member hydrothermal fluid with a dual source of sulfur: sulfur that was leached from basaltic rocks, and sulfur derived from the reduction of seawater sulfate. In contrast, pyrite and marcasite generally have lower δ³⁴S within the range of magmatic values (0 ± 1‰) and are characterized by low concentrations of Se (<50 ppm). For ⁸²Se/⁷⁶Se ratios, the δ⁸²Se values range from basaltic values of near −1.5‰ to −7‰. The large range and highly negative values of hydrothermal deposits observed cannot be explained by simple mixing between Se leached from igneous rock and Se derived from seawater. We interpret the Se isotope signature to be a result of leaching and mixing of a fractionated Se source located beneath hydrothermal chimneys in the hydrothermal fluid. At Lucky Strike we consider two sources for S and Se: (1) the “end-member” hydrothermal fluid with basaltic Se isotopic values (−1.5‰) and typical S isotope hydrothermal values of 1.5‰; (2) a fractionated source hosted in subsurface environment with negative δ³⁴S values, probably from bacterial reduction of seawater sulfate and negative δ⁸²Se values possibly derived from inorganic reduction of Se oxyanions. Fluid trapped in the subsurface environment is conductively cooled and has restricted mixing and provide favorable conditions for subsurface microbial activity which is potentially recorded by S isotopes. Fe isotope systematic reveals that Se-rich high temperature samples have δ⁵⁷Fe values close to basaltic values (∼0‰) whereas Se-depleted samples precipitated at medium to low temperature are systematically lighter (δ⁵⁷Fe values between −1 to −3‰). An important implication of our finding is that light Fe isotope composition down to −3.2‰ may be explained entirely by abiotic fractionation, in which a reservoir effect during sulfide precipitation was able to produce highly fractionated compositions.     

Paleowaters in Silurian-Devonian carbonate aquifers: Geochemical evolution of groundwater in the Great Lakes region since the Late Pleistocene

Changes in the climatic conditions during the Late Quaternary and Holocene greatly impacted the hydrology and geochemical evolution of groundwaters in the Great Lakes region. Increased hydraulic gradients from melting of kilometer-thick Pleistocene ice sheets reorganized regional-scale groundwater flow in Paleozoic aquifers in underlying intracratonic basins. Here, we present new elemental and isotopic analyses of 134 groundwaters from Silurian-Devonian carbonate and overlying glacial drift aquifers, along the margins of the Illinois and Michigan basins, to evaluate the paleohydrology, age distribution, and geochemical evolution of confined aquifer systems. This study significantly extends the spatial coverage of previously published groundwaters in carbonate and drift aquifers across the Midcontinent region, and extends into deeper portions of the Illinois and Michigan basins, focused on the freshwater-saline water mixing zones. In addition, the hydrogeochemical data from Silurian-Devonian aquifers were integrated with deeper basinal fluids, and brines in Upper Devonian black shales and underlying Cambrian-Ordovician aquifers to reveal a regionally extensive recharge system of Pleistocene-age waters in glaciated sedimentary basins. Elemental and isotope geochemistry of confined groundwaters in Silurian-Devonian carbonate and glacial drift aquifers show that they have been extensively altered by incongruent dissolution of carbonate minerals, dissolution of halite and anhydrite, cation exchange, microbial processes, and mixing with basinal brines. Carbon isotope values of dissolved inorganic carbon (DIC) range from −10 to −2‰, ⁸⁷Sr/⁸⁶Sr ratios range from 0.7080 to 0.7090, and δ³⁴S-SO₄ values range from +10 to 30‰. A few waters have elevated δ¹³C_DIC values (>15‰) from microbial methanogenesis in adjacent organic-rich Upper Devonian shales. Radiocarbon ages and δ¹⁸O and δD values of confined groundwaters indicate they originated as subglacial recharge beneath the Laurentide Ice Sheet (14-50 ka BP, −15 to −13‰ δ¹⁸O). These paleowaters are isolated from shallow flow systems in overlying glacial drift aquifers by lake-bed clays and/or shales. The presence of isotopically depleted waters in Paleozoic aquifers at relatively shallow depths illustrates the importance of continental glaciation on regional-scale groundwater flow. Modern groundwater flow in the Great Lakes region is primarily restricted to shallow unconfined glacial drift aquifers. Recharge waters in Silurian-Devonian and unconfined drift aquifers have δ¹⁸O values within the range of Holocene precipitation: −11 to −8‰ and −7 to −4.5‰ for northern Michigan and northern Indiana/Ohio, respectively. Carbon and Sr isotope systematics indicate shallow groundwaters evolved through congruent dissolution of carbonate minerals under open and closed system conditions (δ¹³C_DIC = −14.7 to−11.1‰ and ⁸⁷Sr/⁸⁶Sr = 0.7080-0.7103). The distinct elemental and isotope geochemistry of Pleistocene- versus Holocene-age waters further confirms that surficial flow systems are out of contact with the deeper basinal-scale flow systems. These results provide improved understanding of the effects of past climate change on groundwater flow and geochemical processes, which are important for determining the sustainability of present-day water resources and stability of saline fluids in sedimentary basins.     

A chlorine isotope study of DSDP/ODP serpentinized ultramafic rocks: Insights into the serpentinization process

Eight DSDP/ODP cores were analyzed for major ion concentrations and δ³⁷Cl values of water-soluble chloride (δ³⁷Cl_WSC) and structurally bound chloride (δ³⁷Cl_SBC) in serpentinized ultramafic rocks. This diverse set of cores spans a wide range in age, temperature of serpentinization, tectonic setting, and geographic location of drilled serpentinized oceanic crust. Three of the cores were sampled at closely spaced intervals to investigate downhole variation in Cl concentration and chlorine isotope composition.The average total Cl content of all 86 samples is 0.26 ± 0.16 wt.% (0.19 ± 0.10 wt.% as water-soluble Cl (X_WSC) and 0.09 ± 0.09 wt.% as structurally bound Cl (X_SBC)). Structurally bound Cl concentration nearly doubles with depth in all cores; there is no consistent trend in water-soluble Cl content among the cores. Chlorine isotope fractionation between the structurally bound Cl⁻ site and the water-soluble Cl⁻ site varies from − 1.08‰ to + 1.16‰, averaging to + 0.21‰. Samples with negative fractionations may be related to reequilibration of the water-soluble chloride with seawater post-serpentinite formation. Six of the cores have positive bulk δ³⁷Cl values (+ 0.05‰ to + 0.36‰); the other two cores (173-1068A (Leg-Hole) and 84-570) have negative bulk δ³⁷Cl values (− 1.26‰ and − 0.54‰). The cores with negative δ³⁷Cl values also have variable Cl⁻ / SO₄² ratios, in contrast to all other cores. The isotopically positive cores (153-920D and 147-895E) show no isotopic variation with depth; the isotopically negative core (173-1068A) decreases by ∼1‰ with depth for both the water-soluble and structurally bound Cl fractions.Non-zero bulk δ³⁷Cl values indicate Cl in serpentinites was incorporated during original hydration and is not an artifact of seawater infiltration during drilling. Cores with positive δ³⁷Cl values are most likely explained by open system fractionation during hydrothermal alteration, with preferential incorporation of ³⁷Cl from seawater into the serpentinite and loss of residual light Cl back to the ocean. Fluid / rock ratios were probably low as evidenced by the presence of water-soluble salts. The two isotopically negative cores are characterized by a thick overlying sedimentary package that was in place prior to serpentinization. We believe the low δ³⁷Cl values of these cores are a result of hydration of ultramafic rock by infiltrating aqueous pore fluids from the overlying sediments. The resulting serpentinites inherit the characteristic negative δ³⁷Cl values of the pore waters. Chlorine stable isotopes can be used to identify the source of the serpentinizing fluid and ultimately discern chemical and tectonic processes involved in serpentinization.     

Isotopic constraints on ice age fluids in active geothermal systems: Reykjanes, Iceland

The Reykjanes geothermal system is located on the landward extension of the Mid-Atlantic Ridge in southwest Iceland, and provides an on-land proxy to high-temperature hydrothermal systems of oceanic spreading centers. Previous studies of elemental composition and salinity have shown that Reykjanes geothermal fluids are likely hydrothermally modified seawater. However, δD values of these fluids are as low as −23‰, which is indicative of a meteoric water component. Here we constrain the origin of Reykjanes hydrothermal solutions by analysis of hydrogen and oxygen isotope compositions of hydrothermal epidote from geothermal drillholes at depths between 1 and 3 km. δD_EPIDOTE values from wells RN-8, -9, -10 and -17 collectively range from −60 to −78‰, and δ¹⁸O_EPIDOTE in these wells are between −3.0 and 2.3‰. The δD values of epidote generally increase along a NE trend through the geothermal field, whereas δ¹⁸O values generally decrease, suggesting a southwest to northeast migration of the geothermal upflow zone with time that is consistent with present-day temperatures and observed hydrothermal mineral zones. For comparative analysis, the meteoric-water dominated Nesjavellir and Krafla geothermal systems, which have a δD_FLUID of ∼ −79‰ and −89‰, respectively, show δD_EPIDOTE values of −115‰ and −125‰. In contrast, δD_EPIDOTE from the mixed meteoric-seawater Svartsengi geothermal system is −68‰; comparable to δD_EPIDOTE from well RN-10 at Reykjanes.Stable isotope compositions of geothermal fluids in isotopic equilibrium with the epidotes at Reykjanes are computed using published temperature dependent hydrogen and oxygen isotope fractionation curves for epidote-water, measured isotope composition of the epidotes and temperatures approximated from the boiling point curve with depth. Calculated δD and δ¹⁸O of geothermal fluids are less than 0‰, suggesting that fluids of meteoric or glacial origin are a significant component of the geothermal solutions. Additionally, δD_FLUID values in equilibrium with geothermal epidote are lower than those of modern-day fluids, whereas calculated δ¹⁸O_FLUID values are within range of the observed fluid isotope composition. We propose that modern δD_EPIDOTE and δD_FLUID values are the result of diffusional exchange between hydrous alteration minerals that precipitated from glacially-derived fluids early in the evolution of the Reykjanes system and modern seawater-derived geothermal fluids. A simplified model of isotope exchange in the Reykjanes geothermal system, in which the average starting δD_ROCK value is −125‰ and the water to rock mass ratio is 0.25, predicts a δD_FLUID composition within 1‰ of average measured values. This model resolves the discrepancy between fluid salinity and isotope composition of Reykjanes geothermal fluids, explains the observed disequilibrium between modern fluids and hydrothermal epidote, and suggests that rock-fluid interaction is the dominant control over the evolution of fluid isotope composition in the hydrothermal system.     

Oxygen isotopic studies of the interaction between xenoliths and mafic magma, Voisey’s Bay Intrusion, Labrador, Canada

Sulfide mineralization in the Voisey’s Bay Intrusion, Labrador, Canada, is closely associated with country rock xenoliths that have extensively reacted with basaltic magma. In order to better understand the processes that control the assimilation of country rocks by mafic magma, a detailed study of oxygen isotope systematics related to magma-country rock interaction in the Voisey’s Bay area was undertaken. Protracted interaction of the xenoliths with magma produced refractory mineral assemblages in the xenoliths (2-10 cm in diameter) composed of Ca-rich plagioclase, corundum, hercynite, and minor magnetite. Overgrowth rims of plagioclase and biotite that surround most xenoliths separate the restites from the enclosing igneous matrix. The δ¹⁸O values of minerals from regionally metamorphosed pelitic and quartzofeldspathic protoliths are: plagioclase (8.7-12.3‰), orthoclase (9.5-9.8‰), biotite (5.2-8.7‰), garnet (8.3-10.8‰), pyroxene (8.0-10.1‰), and quartz (9.6-14.0). The δ¹⁸O values of minerals from the hornfels in the contact aureole of the intrusion are consistent with modeling which indicates that as a result of essentially closed system contact metamorphism oxygen isotope values should differ only slightly from those of the protoliths. Hercynite, plagioclase, and corundum separates from the xenoliths have δ¹⁸O values that vary from 2.9‰ to 10.5‰, 5.6‰ to 10.9‰, and 2.0‰ to 6.8‰, respectively. Although a siliceous ¹⁸O-enriched melt has been lost from the xenoliths, corundum, and feldspar δ¹⁸O values are significantly lower than expected through melt loss alone. The relatively low δ¹⁸O values of minerals from the xenoliths may be a function of incomplete isotopic exchange with surrounding mafic magma which had a δ¹⁸O value of ∼5.5‰ to 6.0‰. The high-¹⁸O melt that was released from the xenoliths is partially recorded in the plagioclase overgrowth on the margin of the xenoliths (δ¹⁸O values from 6.2‰ to 10.7‰), and in hercynite that replaced corundum. However, mass balance calculations indicate that a portion of the partial melt must have been transferred to magma that was moving through the conduit system. δ¹⁸O and δD values of biotite surrounding the plagioclase overgrowth range from 5.0‰ to 6.2‰ and −58‰ to −80‰, respectively. These data suggest that the outermost rim associated with many xenoliths has closely approached isotopic equilibrium with uncontaminated mafic magma. The current gabbroic to troctolitic matrix of the xenoliths shows no evidence for contamination by the high-¹⁸O partial melt from the xenoliths. The feldspar and biotite overgrowths on the xenoliths that formed after the motion of the xenoliths relative to the magma had stopped prevented further isotopic exchange between the xenoliths and final magma. The minerals within the xenoliths are not in oxygen isotopic equilibrium with each other, due in part to rapid thermal equilibration, partial melting, and partial exchange with flow through magma.