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.     

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.     

Br/Cl ratios and O, H, C, and B isotopic constraints on the origin of saline waters from eastern Canada

Saline groundwaters were recovered from undisturbed (Restigouche deposit) and active (Brunswick #12 mine) Zn-Pb volcanogenic massive sulfide deposits in the Bathurst Mining Camp (BMC), northern New Brunswick, Canada. These groundwaters, along with fresh to brackish meteoric ground and surface waters from the BMC, have been analyzed to determine their major, trace element and stable isotopic (O, H, C, and B) compositions. Saline groundwaters (total dissolved solids = 22-45 g/L) are characterized by relatively high Na/Ca ratios compared to brines from the Canadian Shield and low Na/Cl_molar and δ¹¹B isotopic compositions (−2.5‰ to 11.1‰) compared to seawater. Although saline waters from the Canadian Shield commonly have oxygen and hydrogen isotopic compositions that plot to the left of the global meteoric water line, those from the BMC fall close to the water line. Fracture and vein carbonate minerals at the Restigouche deposit have restricted carbon isotopic compositions of around −5‰ to −6‰. The carbon isotopic compositions of the saline waters at the Restigouche deposit (+12‰ δ¹³C_DIC) are the result of fractionation of dissolved inorganic carbon by methanogenesis. We suggest that, unlike previous models for shield brines, the composition of saline waters in the BMC is best explained by prolonged water-rock reaction, with no requirement of precursor seawater. We suggest that elevated Br/Cl ratios of saline waters compared to seawater may be explained by differential uptake of Br and Cl during groundwater evolution through water-rock reaction.     

Chlorine stable isotope studies of old groundwater,southwestern Great Artesian Basin,Australia

Stable Cl isotope ratios (³⁷Cl/³⁵Cl) were measured in groundwater samples from the southwestern flow system of the Great Artesian Basin, Australia to gain a better understanding of the Cl⁻ sources and transport mechanisms. δ³⁷Cl values range from 0‰ to −2.5‰ (SMOC), and are inversely correlated with Cl⁻ concentration along the inferred flow direction. The Cl isotopic compositions, in conjunction with other geochemical parameters, suggest that Cl⁻ in groundwaters is not derived from salt dissolution. Mixing of the recharge water with saline groundwater cannot explain the relationship between δ³⁷Cl and Cl⁻ concentration measured. Marine aerosols deposited via rainfall and subsequent evapotranspiration appear to be responsible for the Cl⁻ concentrations observed in wells that are close to the recharge area, and in groundwaters sampled along the southern transect. δ³⁷Cl values measured in the leachate of the Bulldog shale suggest that the aquitard is the subsurface source of Cl⁻ for the majority of groundwater samples studied. Diffusion is likely the mechanism through which Cl⁻ is transported from the pore water of the Bulldog shale to the aquifer. However, a more detailed study of the aquitard rocks is required to verify this hypothesis.     

Arsenic mobility in groundwater/surface water systems in carbonate-rich Pleistocene glacial drift aquifers (Michigan)

Within the Lower Peninsula of Michigan, groundwaters from the Marshall Formation (Mississippian) contain As derived from As-rich pyrites, often exceeding the World Heath Organization drinking water limit of 10 μg/L. Many Michigan watersheds, established on top of Pleistocene glacial drift derived from erosion of the underlying Marshall Formation, also have waters with elevated As. The Huron River watershed in southeastern Lower Michigan is a well characterized hydrogeochemical system of glacial drift deposits, proximate to the Marshall Fm. subcrop, which hosts carbonate-rich groundwaters, streams, and wetlands (fens), and well-developed soil profiles. Aqueous and solid phase geochemistry was determined for soils, soil waters, surface waters (streams and fens) and groundwaters from glacial drift aquifers to better understand the hydrogeologic and chemical controls on As mobility. Soil profiles established on the glacial drift exhibit enrichment in both Fe and As in the oxyhydroxide-rich zone of accumulation. The amounts of Fe and As present as oxyhydroxides are comparable to those reported from bulk Marshall Fm. core samples by previous workers. However, the As host in core samples is largely unaltered pyrite and arsenopyrite. This suggests that the transformation of Fe sulfides to Fe oxyhydroxides largely retains As and Fe at the oxidative weathering site. Groundwaters have the highest As values of all the waters sampled, and many were at or above the World Health limit. Most groundwaters are anaerobic, within the zones of Fe³⁺ and As(V) reduction. Although reduction of Fe(III) oxyhydroxides is the probable source of As, there is no correlation between As and Fe concentrations. The As/Fe mole ratios in drift groundwaters are about an order of magnitude greater than those in soil profiles, suggesting that As is more mobile than Fe. This is consistent with the dominance of As(III) in these groundwaters and with the partitioning of Fe²⁺ into carbonate cements. Soil waters have very low As and Fe contents, consistent with the stability of oxyhydroxides under oxidizing vadose conditions. When CO₂ charged groundwaters discharge in streams and fens, dissolved As is effectively removed by adsorption onto Fe-oxides or carbonate marls. Although Fe does not display conservative behavior with As in groundwaters, a strong positive correlation exists between As and Sr concentrations. As water–rock interactions proceed, the As/Fe and Sr/Ca ratios would be expected to increase because both As and Sr behave as incompatible elements. Comparisons with groundwater chemistries from other drift-hosted aquifers proximate to the Marshall sandstone are consistent with these relations. Thus, the Sr content of carbonate-rich groundwaters may provide useful constraints on the occurrence, origin and evolution of dissolved As in such systems.     

Magnesium isotope systematics of the lithologically varied Moselle river basin, France

Magnesium and strontium isotope signatures were determined during different seasons for the main rivers of the Moselle basin, northeastern France. This small basin is remarkable for its well-constrained and varied lithology on a small distance scale, and this is reflected in river water Sr isotope compositions. Upstream, where the Moselle River drains silicate rocks of the Vosges mountains, waters are characterized by relatively high ⁸⁷Sr/⁸⁶Sr ratios (0.7128-0.7174). In contrast, downstream of the city of Epinal where the Moselle River flows through carbonates and evaporites of the Lorraine plateau, ⁸⁷Sr/⁸⁶Sr ratios are lower, down to 0.70824.Magnesium in river waters draining silicates is systematically depleted in heavy isotopes (δ²⁶Mg values range from −1.2 to −0.7‰) relative to the value presently estimated for the continental crust and a local diorite (−0.5‰). In comparison, δ²⁶Mg values measured in soil samples are higher (∼0.0‰). This suggests that Mg isotope fractionation occurs during mineral leaching and/or formation of secondary clay minerals. On the Lorraine plateau, tributaries draining marls, carbonates and evaporites are characterized by low Ca/Mg (1.5-3.2) and low Ca/Sr (80-400) when compared to local carbonate rocks (Ca/Mg = 29-59; Ca/Sr = 370-2200), similar to other rivers draining carbonates. The most likely cause of the Mg and Sr excesses in these rivers is early thermodynamic saturation of groundwater with calcite relative to magnesite and strontianite as groundwater chemistry progressively evolves in the aquifer. δ²⁶Mg of the dissolved phases of tributaries draining mainly carbonates and evaporites are relatively low and constant throughout the year (from −1.4‰ to −1.6‰ and from −1.2‰ to −1.4‰, respectively), within the range defined for the underlying rocks. Downstream of Epinal, the compositions of the Moselle River samples in a δ²⁶Mg vs. ⁸⁷Sr/⁸⁶Sr diagram can be explained by mixing curves between silicate, carbonate and evaporite waters, with a significant contribution from the Vosgian silicate lithologies (>70%). Temporal co-variation between δ²⁶Mg and ⁸⁷Sr/⁸⁶Sr for the Moselle River throughout year is also observed, and is consistent with a higher contribution from the Vosges mountains in winter, in terms of runoff and dissolved element flux. Overall, this study shows that Mg isotopes measured in waters, rocks and soils, coupled with other tracers such as Sr isotopes, could be used to better constrain riverine Mg sources, particularly if analytical uncertainties in Mg isotope measurements can be improved in order to perform more precise quantifications.     

Iron isotope fractionation in subterranean estuaries

Dissolved Fe concentrations in subterranean estuaries, like their river-seawater counterparts, are strongly controlled by non-conservative behavior during mixing of groundwater and seawater in coastal aquifers. Previous studies at a subterranean estuary of Waquoit Bay on Cape Cod, USA demonstrate extensive precipitation of groundwater-borne dissolved ferrous iron and subsequent accumulation of iron oxides onto subsurface sands. Waquoit Bay is thus an excellent natural laboratory to assess the mechanisms of Fe-isotope fractionation in redox-stratified environments and determine potential Fe-isotope signatures of groundwater sources to coastal seawater. Here, we report Fe isotope compositions of iron-coated sands and porewaters beneath the intertidal zone of Waquoit Bay. The distribution of pore water Fe shows two distinct sources of Fe: one residing in the upward rising plume of Fe-rich groundwater and the second in the salt-wedge zone of pore water. The groundwater source has high Fe(II) concentration consistent with anoxic conditions and yield δ⁵⁶Fe values between 0.3 and −1.3‰. In contrast, sediment porewaters occurring in the mixing zone of the subterranean estuary have very low δ⁵⁶Fe values down to −5‰. These low δ⁵⁶Fe values reflect Fe-redox cycling and result from the preferential retention of heavy Fe-isotopes onto newly formed Fe-oxyhydroxides. Analysis of Fe-oxides precipitated onto subsurface sands in two cores from the subterranean estuary revealed strong δ⁵⁶Fe and Fe concentration gradients over less than 2m, yielding an overall range of δ⁵⁶Fe values between −2 and 1.5‰. The relationship between Fe concentration and δ⁵⁶Fe of Fe-rich sands can be modeled by the progressive precipitation of Fe-oxides along fluid flow through the subterranean estuary. These results demonstrate that large-scale Fe isotope fractionation (up to 5‰) can occur in subterranean estuaries, which could lead to coastal seawater characterized by very low δ⁵⁶Fe values relative to river values.     

Mg isotope constraints on soil pore-fluid chemistry: Evidence from Santa Cruz, California

Mg isotope ratios (²⁶Mg/²⁴Mg) are reported in soil pore-fluids, rain and seawater, grass and smectite from a 90 kyr old soil, developed on an uplifted marine terrace from Santa Cruz, California. Rain water has an invariant ²⁶Mg/²⁴Mg ratio (expressed as δ²⁶Mg) at −0.79 ± 0.05‰, identical to seawater δ²⁶Mg. Detrital smectite (from the base of the soil profile, and therefore unweathered) has a δ²⁶Mg value of 0.11‰, potentially enriched in ²⁶Mg by up to 0.3‰ compared to the bulk silicate Earth Mg isotope composition (although within the range of all terrestrial silicates). The soil pore-waters show a continuous profile with depth for δ²⁶Mg, ranging from −0.99‰ near the surface to −0.43‰ at the base of the profile. Shallow pore-waters (<1 m) have δ²⁶Mg values that are similar to, or slightly lower than the rain waters. This implies that the degree of biological cycling of Mg in the pore-waters is relatively small and is quantified as <32%, calculated using the average Mg isotope enrichment factor between grass and rain (δ²⁶Mg_grass-δ²⁶Mg_rain) of 0.21‰. The deep pore-waters (1-15 m deep) have δ²⁶Mg values that are intermediate between the smectite and rain, ranging from −0.76‰ to −0.43‰, and show a similar trend with depth compared to Sr isotope ratios. The similarity between Sr and Mg isotope ratios confirms that the Mg in the pore-waters can be explained by a mixture between rain and smectite derived Mg, despite the fact that Mg and Sr concentrations may be buffered by the exchangeable reservoir. However, whilst Sr isotope ratios in the pore-waters span almost the complete range between mineral and rain inputs, Mg isotopes compositions are much closer to the rain inputs. If Mg and Sr isotope ratios are controlled uniquely by a mixture, the data can be used to estimate the mineral weathering inputs to the pore-waters, by correcting for the rain inputs. This isotopic correction is compared to the commonly used chloride correction for precipitation inputs. A consistent interpretation is only possible if Mg isotope ratios are fractionated either by the precipitation of a secondary Mg bearing phase, not detected by conventional methods, or selective leaching of ²⁴Mg from smectite. There is therefore dual control on the Mg isotopic composition of the pore-waters, mixing of two inputs with distinct isotopic compositions, modified by fractionation. The data provide (1) further evidence for Mg isotope fractionation at the surface of the Earth and (2) the first field evidence of Mg isotope fractionation during uptake by natural plants. The coherent behaviour of Mg isotope ratios in soil environments is encouraging for the development of Mg isotope ratios as a quantitative tracer of both weathering inputs of Mg to waters, and the physicochemical processes that cycle Mg, a major cation linked to the carbon cycle, during continental weathering.     

Geothermal waters from the Taupo Volcanic Zone,New Zealand: Li,B and Sr isotopes characterization

Chemical and isotopic data for 23 geothermal water samples collected in New Zealand within the Taupo Volcanic Zone (TVZ) are reported. Major and trace elements including Li, B and Sr and their isotopic compositions (δ⁷Li, δ¹¹B, ⁸⁷Sr/⁸⁶Sr) were determined in high temperature geothermal waters collected from deep boreholes in different geothermal fields (Ohaaki, Wairakei, Mokai, Kawerau and Rotokawa geothermal systems). Lithium concentrations are high (from 4.5 to 19.9 mg/L) and Li isotopic compositions (δ⁷Li) are homogeneous, ranging between −0.5‰ and +1.4‰. In particular, it is noteworthy that, except for the samples from the Kawerau geothermal field having slightly higher δ⁷Li values (+1.4%), the other geothermal waters have a near constant δ⁷Li signature around a mean value of 0‰ ± 0.6 (2σ, n = 21). Boron concentrations are also high and relatively homogeneous for the geothermal samples, falling between 17.5 and 82.1 mg/L. Boron isotopic compositions (δ¹¹B) are all negative, and display a range between −6.7‰ and −1.9‰. These B isotope compositions are in agreement with those of the Ngawha geothermal field in New Zealand. Lithium and B isotope signatures are in a good agreement with a fluid signature mainly derived from water/rock interaction involving magmatic rocks with no evidence of seawater input. On the other hand, Sr concentrations are lower and more heterogeneous and fall between 2 and 165 μg/L. The ⁸⁷Sr/⁸⁶Sr ratios range from 0.70549 to 0.70961. These Sr isotope compositions overlap those of the Rotorua geothermal field in New Zealand, confirming that some geothermal waters (with more radiogenic Sr) have interacted with bedrocks from the metasedimentary basement. Each of these isotope systems on their own reveals important information about particular aspects of either water source or water/rock interaction processes, but, considered together, provide a more integrated understanding of the geothermal systems from the TVZ in New Zealand.     

Application of carbon isotopes to detect seepage out of coalbed natural gas produced water impoundments

Coalbed natural gas (CBNG) production from coal bed aquifers requires large volumes of produced water to be pumped from the subsurface. The produced water ranges from high quality that meets state and federal drinking water standards to low quality due to increased salinity and/or sodicity. The Powder River Basin of northeastern Wyoming is a major coalbed natural gas producing region, where water quality generally decreases moving from the southeastern portion of the basin towards the center. Most produced water in Wyoming is disposed into impoundments and other surface drainages, where it may infiltrate into shallow groundwater. Groundwater degradation caused by infiltration of CBNG produced water holding impoundments into arid, soluble salt-rich soils is an issue of immense importance because groundwater is a major source for stock water, irrigation, and drinking water for many small communities in these areas. This study examines the potential of using stable C isotope signatures of dissolved inorganic C (δ¹³C_DIC) to track the fate of CBNG produced water after it is discharged into the impoundments. Other geochemical proxies like the major cations and major anions were used in conjunction with field water quality measurements to understand the geochemical differences between CBNG produced waters and ambient waters in the study area. Samples were collected from the CBNG discharge outfalls, produced water holding impoundments, and monitoring wells from different parts of the Powder River Basin and analyzed for δ¹³C_DIC. The CBNG produced waters from outfalls and impoundments have positive δ¹³C_DIC values that fall within the range of +12‰ to +22‰, distinct from the ambient regional surface and groundwaters with δ¹³C_DIC values ranging from −10‰ to −14‰. The results from the study demonstrate that these contrasting δ¹³C_DIC signatures can be used to trace seepage out of CBNG produced water impoundments into shallow groundwaters.     

Carbon isotope systematics of the Cambrian–Vendian aquifer system in the northern Baltic Basin: Implications to the age and evolution of groundwater

Groundwater in the Cambrian–Vendian aquifer system has a strongly depleted stable isotope composition (δ¹⁸O values of about −22‰) and a low radiocarbon concentration, which suggests that the water is of glacial origin from the last Ice Age. The aim of this paper was to elucidate the timing of infiltration of glacial waters and to understand the geochemical evolution of this groundwater. The composition of the dissolved inorganic C (DIC) in Cambrian–Vendian groundwater is influenced by complex reactions and isotope exchange processes between water, organic materials and rock matrix. The δ¹³C composition of dissolved inorganic C in Cambrian–Vendian water also indicates a bacterial modification of the isotope system. The corrected radiocarbon ages of groundwater are between 14,000 and 27,000 radiocarbon years, which is coeval with the advance of the Weichselian Glacier in the area.     

The sources and budget for dissolved sulfate in a fractured carbonate aquifer,southern Sacramento Mountains,New Mexico,USA

Climate change in the SW USA is likely to involve drier conditions and higher surface temperatures. In order to better understand the evolution of water chemistry and the sources of aqueous SO₄ in these semi-arid settings, chemical and S isotope compositions were determined of springs, groundwater, and bedrock associated with a Permian fractured carbonate aquifer located in the southern Sacramento Mountains, New Mexico, USA. The results suggest that the evolution of water chemistry in the semi-arid carbonate aquifer is mainly controlled by dedolomitization of bedrock, which was magnified by increasing temperature and increasing dissolution of gypsum/anhydrite along the groundwater flow path. The δ³⁴S of dissolved SO₄ in spring and groundwater samples varied from +9.0‰ to +12.8‰, reflecting the mixing of SO₄ from the dissolution of Permian gypsum/anhydrite (+12.3‰ to +13.4‰) and oxidation of sulfide minerals (−24.5‰ to −4.2‰). According to S isotope mass balance constraints, the contribution of sulfide-derived SO₄ was considerable in the High Mountain recharge areas, accounting for up to ∼10% of the total SO₄ load. However, sulfide weathering decreased in importance in the lower reaches of the watershed. A smaller SO₄ input of ∼2–4% was contributed by atmospheric wet deposition. This study implies that the δ³⁴S variation of SO₄ in semi-arid environments can be complex, but that S isotopes can be used to distinguish among the different sources of weathering. Here it was found that H₂SO₄ dissolution due to sulfide oxidation contributes up to 5% of the total carbonate weathering budget, while most of the SO₄ is released from bedrock sources during dedolomitization.     

The origin of thermal waters in the northeastern part of the Eger Rift,Czech Republic

An investigation of the thermal waters in the Ústí nad Labem area in the northeastern part of the Eger Rift has been carried out, with the principal objective of determining their origin. Waters from geothermal reservoirs in the aquifers of the Bohemian Cretaceous Basin (BCB) from depths of 240 to 616 m are exploited here. For comparison, thermal waters of the adjacent Teplice Spa area were also incorporated into the study. Results based on water chemistry and isotopes indicate mixing of groundwater from aquifers of the BCB with groundwater derived from underlying crystalline rocks of the Erzgebirge Mts. Unlike thermal waters in Dě?ín, which are of Ca–HCO₃ type, there are two types of thermal waters in Ústí nad Labem, Na–HCO₃–Cl–SO₄ type with high TDS values and Na–Ca–HCO₃–SO₄ type with low TDS values. Carbon isotope data, speciation calculations, and inverse geochemical modeling suggest a significant input of endogenous CO₂ at Ústí nad Labem in the case of high TDS groundwaters. Besides CO₂ input, both silicate dissolution and cation exchange coupled with dissolution of carbonates may explain the origin of high TDS thermal waters equally well. This is a consequence of similar δ¹³C and ¹⁴C values in endogenous CO₂ and carbonates (both sources have ¹⁴C of 0 pmc, endogenous CO₂ δ¹³C around −3‰, carbonates in the range from −5‰ to +3‰ V-PDB). The source of Cl⁻ seems to be relict brine formed in Tertiary lakes, which infiltrated into the deep rift zone and is being flushed out. The difference between high and low TDS groundwaters in Ústí nad Labem is caused by location of the high mineralization groundwater wells in CO₂ emanation centers linked to channel-like conduits. This results in high dissolution rates of minerals and in different δ¹³C(DIC) and ¹⁴C(DIC) fingerprints. A combined δ³⁴S and δ¹⁸O study of dissolved SO₄ indicates multiple SO₄ sources, involving SO₄ from relict brines and oxidation of H₂S. The study clearly demonstrates potential problems encountered at sites with multiple sources of C, where several evolutionary groundwater scenarios are possible.     

The use of O,H, B,Sr and S isotopes for tracing the origin of dissolved boron in groundwater in Central Macedonia,Greece

The groundwater B concentration in Mesozoic karst, Neogene and alluvial aquifers in the West part of Chalkidiki province in Central Macedonia, Greece reaches 6.45 mg L⁻¹, which exceeds the limit of 1 mg L⁻¹, set by the European Union for drinking water. The high B contents have been detected in this area, not only near the shoreline, where seawater intrusion occurs, but also in the inland part of the basin. Multi isotope (²H, ¹⁸O, ³⁴S, ¹⁸O_(SO4), ¹¹B, ⁸⁷Sr/⁸⁶Sr) data from borehole and thermal water springs allow identification of the possible B sources. The B dissolved in groundwater in the Chalkidiki area is mainly geogenic. The low δ¹¹B values, 0–1‰, similar to those of thermal fluids from continental geothermal fields, and the low Cl/B ratio compared to seawater both indicate a geothermal origin for B and reflect deep circulation and interaction with igneous rocks. The ⁸⁷Sr/⁸⁶Sr ratio also indicates that the deep-aquifer granodiorite is the predominant rock source of Sr, while the shallow limestone unit has negligible effects on the dissolved Sr budget in these thermal karst waters which O and H isotopes show to be of meteoric origin. The main source of high B in borehole water is mainly mixing with B-rich geothermal water. The mixing between geothermal water and water from the Neogene aquifer is also reflected by isotopic contents of SO₄.     

Pore water evolution in oilfield sandstones: constraints from oxygen isotope microanalyses of quartz cement

Oxygen isotope microanalyses of authigenic quartz, in combination with temperatures of quartz precipitation constrained by fluid inclusion microthermometry and burial history modelling, are employed to trace the origin and evolution of pore waters in three distinct reservoirs of the Brae Formation in the Miller and Kingfisher Fields (North Sea). Oxygen isotope ratios of quartz cements were measured in situ in nine sandstone thin sections with a Cameca ims-4f ion microprobe. In conjunction with quartz cement paragenesis in the reservoirs, constrained from textural and cathodoluminescence (CL) microscopy studies, pore water evolution was reconstructed from the time of deposition of the sandstones in the Upper Jurassic until the present.CL photomicrographs of quartz overgrowths in the Brae Formation sandstones show three cement zones (A, B and C) which can be related to different oxygen isotope compositions: (1) the earliest, and thinnest, zone A (homogeneous CL pattern with probable δ¹⁸O values between +23‰ and +26‰—direct measurements were not possible) precipitated in the sandstones at temperatures <60 °C; (2) the second zone B (complex CL pattern and directly measured δ¹⁸O values between +15‰ and +18‰) precipitated in the sandstones most likely between 70 and 90 °C; (3) the third zone C (homogeneous CL pattern and directly measured δ¹⁸O values between +16‰ and +22‰) precipitated in the sandstones most likely at temperatures >90 °C. Calculated oxygen isotope compositions of pore waters show that zone A quartz cements, and enclosing concretionary calcite, precipitated from a meteoric-type fluid (∼−7‰) during shallow burial (<1.5 km). Zone B quartz cements precipitated from fluids which evolved in composition from a meteoric-type fluid (δ¹⁸O −7‰) to a more ¹⁸O-enriched fluid (δ¹⁸O −4‰) as burial continued to ∼3.0 km. Data from zone C quartz cements are consistent with further fluid evolution from δ¹⁸O −4‰ to basinal-type fluids with δ¹⁸O similar to the present-day formation water oxygen isotope composition (+0.6‰ at 4.0 km burial). A similar pore water evolution can be derived for all three reservoirs studied, indicating that hydrogeologic evolution was similar across sandstones of the whole Brae Formation.The quartz cement zones observed in the Brae Formation sandstones, and the pore water history derived for the area studied, is analogous to published petrographic and pore water evolution data from the nearby Brent Group reservoirs and from reservoirs located in the Haltenbanken area on the Atlantic margin offshore Norway. Considering quartz cement is a major porosity-occluding phase in many reservoir sandstones, and because pore waters both dissolve quartz and carry the dissolved silica to cementation sites, the data presented are valuable for improving the understanding and prediction of reservoir quality development in sandstones globally.     

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.     

Note on the isotopic geochemistry of fossil-lacustrine tufas in carbonate plateau—A study from Dungul region (SW Egypt)

Lithological, chemical, and stable isotope data are used to characterize lacustrine tufas dating back to pre-late Miocene and later unknown times, capping different surfaces of a Tertiary carbonate (Sinn el-Kedab) plateau in Dungul region in the currently hyperarid south-western Egypt. These deposits are composed mostly of calcium carbonate, some magnesium carbonate and clastic particles plus minor amounts of organic matter. They have a wide range of (Mg/Ca)_molar ratios, from 0.03 to 0.3. The bulk-tufa carbonate has characteristic isotope compositions: (δ¹³C_mean = −2.49 ± 0.99‰; δ¹⁸O_mean = −9.43 ± 1.40‰). The δ¹³C values are consistent with a small input from C4 vegetation or thinner soils in the recharge area of the tufa-depositing systems. The δ¹⁸O values are typical of fresh water carbonates. Covariation between δ¹³C and δ¹⁸O values probably is a reflection of climatic conditions such as aridity. The tufas studied are isotopically similar to the underlying diagenetic marine chalks, marls and limestones (δ¹³C_mean = −2.06 ± 0.84‰; δ¹⁸O_mean = −10.06 ± 1.39‰). The similarity has been attributed to common meteoric water signatures. This raises large uncertainties in using tufas (Mg/Ca)_molar, δ¹³C and δ¹⁸O records as proxies of paleoclimatic change and suggests that intrinsic compositional differences in material sources within the plateau may mask climatic changes in the records.     

Pleistocene recharge to midcontinent basins: effects on salinity structure and microbial gas generation

The hydrogeochemistry of saline-meteoric water interface zones in sedimentary basins is important in constraining the fluid migration history, chemical evolution of basinal brines, and physical stability of saline formation waters during episodes of freshwater recharge. This is especially germane for interior cratonic basins, such as the Michigan and Illinois basins. Although there are large differences in formation water salinity and hydrostratigraphy in these basins, both are relatively quiescent tectonically and have experienced repeated cycles of glaciation during the Pleistocene. Exploration for unconventional microbial gas deposits, which began in the upper Devonian-age Antrim Shale at the northern margin of the Michigan Basin, has recently extended into the age-equivalent New Albany Shale of the neighboring Illinois Basin, providing access to heretofore unavailable fluid samples. These reveal an extensive regional recharge system that has profoundly changed the salinity structure and induced significant biogeochemical modification of formation water elemental and isotope geochemistry.New-formation water and gas samples were obtained from Devonian-Mississippian strata in the Illinois Basin. These included exploration wells in the New Albany Shale, an organic-rich black shale of upper Devonian age, and formation waters from over- and underlying regional aquifer systems (Siluro-Devonian and Mississippian age). The hydrostratigraphic relations of major aquifers and aquitards along the eastern margin of the Illinois Basin critically influenced fluid migration into the New Albany Shale. The New Albany Shale formation water chemistry indicates significant invasion of meteoric water, with δD values as low as −46.05‰, into the shale. The carbon stable isotope system (δ¹³C values as high as 29.4‰), coupled with δ¹⁸O, δD, and alkalinity of formation waters (alkalinity ≤24.08 meq/kg), identifies the presence of microbial gas associated with meteoric recharge. Regional geochemical patterns identify the underlying Siluro-Devonian carbonate aquifer system as the major conduit for freshwater recharge into the fractured New Albany Shale reservoirs. Recharge from overlying Mississippian carbonates is only significant in the southernmost portion of the basin margin where carbonates directly overlie the New Albany Shale.Recharge of dilute waters (Cl⁻ <1000 mM) into the Siluro-Devonian section has suppressed formation water salinity to depths as great as 1 km across the entire eastern Illinois Basin margin. Taken together with salinity and stable isotope patterns in age-equivalent Michigan Basin formation waters, they suggest a regional impact of recharge of δ¹⁸O- and δD-depleted fluids related to Pleistocene glaciation. Devonian black shales at both basin margins have been affected by recharge and produced significant volumes of microbial methane. This recharge is also manifested in different salinity gradients in the two basins because of their large differences in original formation water salinity. Given the relatively quiet tectonic history and subdued current topography in the midcontinent region, it is likely that repeated cycles of glacial meltwater invasion across this region have induced a strong disequilibrium pattern in fluid salinity and produced a unique class of unconventional shale-hosted gas deposits.