Recharge Areas and Geochemical Evolution of Groundwater in an Alluvial Aquifer System in the Sultanate of Oman

A regional hydrogeochemical model was developed to evaluate the geochemical evolution of different groundwaters in an alluvial aquifer system in the Interior of Oman. In combination with environmental isotopes the model is able to extract qualitative and quantitative information about recharge, groundwater flow paths and hydraulic connections between different aquifers. The main source of water to the alluvial aquifer along the flow paths ofWadi Abyadh andWadi M’uaydin in the piedmont is groundwater from the high-altitude areas of the Jabal Akhdar and local infiltration along the wadi channels. In contrast, the piedmont alluvial aquifer alongWadi Halfayn is primarily replenished by lateral recharge from the ophiolite foothills to the east besides smaller contributions from the Jabal Akhdar and local infiltration. Further down gradient in the Southern Alluvial Plain aquifer a significant source of recharge is direct infiltration of rain and surface runoff, originating from a moisture source that approaches Oman from the south. The model shows that the main geochemical evolution of the alluvial groundwaters occurs along the flow path from the piedmont to the Southern Alluvial Plain, where dedolomitization is responsible for the observed changes in the chemical and carbon isotope composition in these waters.

Time-lapse pressure tomography for characterizing CO2 plume evolution in a deep saline aquifer

A time-lapse pressure tomography inversion approach is applied to characterize the CO2 plume development in a virtual deep saline aquifer. Deep CO2 injection leads to flow properties of the mixed-phase, which vary depending on the CO2 saturation. Analogous to the crossed ray paths of a seismic tomographic experiment, pressure tomography creates streamline patterns by injecting brine prior to CO2 injection or by injecting small amounts of CO2 into the two-phase (brine and CO2) system at different depths. In a first step, the introduced pressure responses at observation locations are utilized for a computationally rapid and efficient eikonal equation based inversion to reconstruct the heterogeneity of the subsurface with diffusivity (D) tomograms. Information about the plume shape can be derived by comparing D-tomograms of the aquifer at different times. In a second step, the aquifer is subdivided into two zones of constant values of hydraulic conductivity (K) and specific storage (Ss) through a clustering approach. For the CO2 plume, mixed-phase K and Ss values are estimated by minimizing the difference between calculated and “true” pressure responses using a single-phase flow simulator to reduce the computing complexity. Finally, the estimated flow property is converted to gas saturation by a single-phase proxy, which represents an integrated value of the plume. This novel approach is tested first with a doublet well configuration, and it reveals a great potential of pressure tomography based concepts for characterizing and monitoring deep aquifers, as well as the evolution of a CO2 plume. Still, field-testing will be required for better assessing the applicability of this approach.

(234)U/(238)U activity ratio disequilibrium technique for studying uranium mobility in the Opalinus Clay at Mont Terri, Switzerland

Mobility of naturally occurring U-238 and U-234 radionuclides was studied in a low permeability, reducing claystone formation (Opalinus Clay) near its contact with an overlying oxidising aquifer (Dogger Limestones) at Mont Terri, Switzerland. Our data point to a limited redistribution of U in some of the studied samples. Observed centimetre-scale U mobility is explained by slow diffusive transport of U-234 in the pore waters of the Opalinus Clay driven by spatially variable in situ supply (by alpha-recoil) of U-234 from the rock matrix. Metre-scale mobility is interpreted as a result of infiltration of meteoric water into the overlying aquifer which developed gradients of U concentration across the two rock formations. This triggered a slow in-diffusion of U with (U-234/U-238) > 1 into the Opalinus Clay as attested by a clear-cut pattern of decreasing bulk rock (U-234/U-238) inwards the Opalinus Clay, away from the Dogger Limestones.

Groundwater age distributions at a public drinking water supply well field derived from multiple age tracers (⁸⁵Kr, ³H/³He, and ³⁹Ar)

Groundwater age is a key aspect of production well vulnerability. Public drinking water supply wells typically have long screens and are expected to produce a mixture of groundwater ages. The groundwater age distributions of seven production wells of the Holten well field (Netherlands) were estimated from tritium-helium (3H/3He), krypton-85 (85Kr), and argon-39 (39Ar), using a new application of a discrete age distribution model and existing mathematical models, by minimizing the uncertainty-weighted squared differences of modeled and measured tracer concentrations. The observed tracer concentrations fitted well to a 4-bin discrete age distribution model or a dispersion model with a fraction of old groundwater. Our results show that more than 75 of the water pumped by four shallow production wells has a groundwater age of less than 20 years and these wells are very vulnerable to recent surface contamination. More than 50 of the water pumped by three deep production wells is older than 60 years. 3H/3He samples from short screened monitoring wells surrounding the well field constrained the age stratification in the aquifer. The discrepancy between the age stratification with depth and the groundwater age distribution of the production wells showed that the well field preferentially pumps from the shallow part of the aquifer. The discrete groundwater age distribution model appears to be a suitable approach in settings where the shape of the age distribution cannot be assumed to follow a simple mathematical model, such as a production well field where wells compete for capture area.

Dating groundwater in the Bohemian Cretaceous Basin: Understanding tracer variations in the subsurface

The northern section of the Bohemian Cretaceous Basin has been the site of intensive U exploitation with harmful impacts on groundwater quality. The understanding of groundwater flow and age distribution is crucial for the prediction of the future dispersion and impact of the contamination. State of the art tracer methods (3H, 3He, 4He, 85Kr, 39Ar and 14C) were, therefore, used to obtain insights to ageing and mixing processes of groundwater along a north–south flow line in the centre of the two most important aquifers of Cenomanian and middle Turonian age. Dating of groundwater is particularly complex in this area as: (i) groundwater in the Cenomanian aquifer is locally affected by fluxes of geogenic and biogenic gases (e.g. CO2, CH4, He) and by fossil brines in basement rocks rich in Cl and SO4; (ii) a thick unsaturated zone overlays the Turonian aquifer; (iii) a periglacial climate and permafrost conditions prevailed during the Last Glacial Maximum (LGM), and iv) the wells are mostly screened over large depth intervals. Large disagreements in 85Kr and 3H/3He ages indicate that processes other than ageing have affected the tracer data in the Turonian aquifer. Mixing with older waters (>50 a) was confirmed by 39Ar activities. An inverse modelling approach, which included time lags for tracer transport throughout the unsaturated zone and degassing of 3He, was used to estimate the age of groundwater. Best fits between model and field results were obtained for mean residence times varying from modern up to a few hundred years. The presence of modern water in this aquifer is correlated with the occurrence of elevated pollution (e.g. nitrates). An increase of reactive geochemical indicators (e.g. Na) and radiogenic 4He, and a decrease in 14C along the flow direction confirmed groundwater ageing in the deeper confined Cenomanian aquifer. Radiocarbon ages varied from a few hundred years to more than 20 ka. Initial 14C activity for radiocarbon dating was calibrated by means of 39Ar measurements. The 14C age of a sample recharged during the LGM was further confirmed by depleted stable isotope signatures and near freezing point noble gas temperature. Radiogenic 4He accumulated in groundwater with concentrations increasing linearly with 14C ages. This enabled the use of 4He to validate the dating range of 14C and extend it to other parts of this aquifer. In the proximity of faults, 39Ar in excess of modern concentrations and 14C dead CO2 sources, elevated 3He/4He ratios and volcanic activity in Oligocene to Quaternary demonstrate the influence of gas of deeper origin and impeded the application of 4He, 39Ar and 14C for groundwater dating.

Assessing the Environmental Hazard of Using Seawater for Ore Processing at the Lasail Mine Site in the Sultanate of Oman

The Lasail mining area (Sultanate of Oman) was contaminated by acid mine drainage during the exploitation and processing of local and imported copper ore and the subsequent deposition of sulphide-bearing waste material into an unsealed tailings dump. In this arid environment, the use of seawater in the initial stages of ore processing caused saline contamination of the fresh groundwater downstream of the tailings dump. After detection of the contamination in the 1980s, different source-controlled remediation activities were conducted including a seepage water collection system and, in 2005, surface sealing of the tailings dump using an HDPE-liner to prevent further infiltration of meteoric water. We have been assessing the benefits of the remediation actions undertaken so far. We present chemical and isotopic (δ18O, δ 2H, 3H) groundwater data from a long-term survey (8–16 years) of the Wadi Suq aquifer along a 28 km profile from the tailings dump to the Gulf of Oman. Over this period, most metal concentrations in the Wadi Suq groundwater decreased below detection limits. In addition, in the first boreholes downstream of the tailings pond, the salinity contamination has decreased by 30 % since 2005. This decrease appears to be related to the surface coverage of the tailings pond, which reduces flushing of the tailings by the sporadic, but commonly heavy, precipitation events. Despite generally low metal concentrations and the decreased salinity, groundwater quality still does not meet the WHO drinking water guidelines in more than 90 % of the Wadi Suq aquifer area. The observations show that under arid conditions, use of seawater for ore processing or any other industrial activity has the potential to contaminate aquifers for decades.

Natural Tracers in recent groundwaters from different Alpine aquifers

Groundwater with underground residence times between days and a few years have been investigated over more than 20 years from 487 remote sites located in different aquifer types in the Alpine belt. Analysis of the data reveals that groundwaters evolved in crystalline, evaporite, carbonate, molasse, and flysch aquifers can be clearly distinguished based on their major and trace element composition and degree of mineralisation. A further subdivision can be made even within one aquifer type based on the trace element compositions, which are characteristic for the lithologic environment. Major and trace element concentrations can be quantitatively described by interaction of the groundwater with the aquifer- specific mineralogy along the flow path. Because all investigated sites show minimal anthropogenic influences, the observed concentration ranges represent the natural background concentrations and can thus serve as a “geo-reference” for recent groundwaters from these five aquifer types. This “geo-reference” is particularly useful for the identification of groundwater contamination. It further shows that drinking water standards can be grossly exceeded for critical elements by purely natural processes

Investigation of young water inflow in karst aquifers using SF6–CFC–3H/He–85Kr–39Ar and stable isotope components

Karst aquifers are known for their wide distribution of water transfer velocities. From this observation, a multiple geochemical tracer approach seems to be particularly well suited to provide a significant assessment of groundwater flows, but the choice of adapted tracers is essential. In this study, several common tracers in karst aquifers such as physicochemical parameters, major ions, stable isotopes, and d13C to more specific tracers such as dating tracers – 14C, 3H, 3H–3He, CFC-12, SF6 and 85Kr, and 39Ar – were used, in a fractured karstic carbonated aquifer located in Burgundy (France). The information carried by each tracer and the best sampling strategy are compared on the basis of geochemical monitoring done during several recharge events and over longer time periods (months to years). This study’s results demonstrate that at the seasonal and recharge event time scale, the variability of concentrations is low for most tracers due to the broad spectrum of groundwater mixings. The tracers used traditionally for the study of karst aquifers, i.e., physicochemical parameters and major ions, efficiently describe hydrological processes such as the direct and differed recharge, but require being monitored at short time steps during recharge events to be maximized. From stable isotopes, tritium, and Cl� contents, the proportion of the fast direct recharge by the largest porosity was estimated using a binary mixing model. The use of tracers such as CFC-12, SF6, and 85Kr in karst aquifers provides additional information, notably an estimation of apparent age, but they require good preliminary knowledge of the karst system to interpret the results suitably. The CFC-12 and SF6 methods efficiently determine the apparent age of baseflow, but it is preferable to sample the groundwater during the recharge event. Furthermore, these methods are based on different assumptions such as regional enrichment in atmospheric SF6, excess air, and flow models among others. 85Kr and 39Ar concentrations can potentially provide a more direct estimation of groundwater residence time. Conversely, the 3H–3He method is inefficient in the karst aquifer for dating due to 3He degassing.

Continental degassing of 4He by surficial discharge of deep groundwater

Radiogenic He is produced by the decay of uranium and thorium in the Earth’s mantle and crust. From here, it is degassed to the atmosphere and eventually escapes to space. Assuming that all of the 4He produced is degassed, about 70% of the total He degassed from Earth comes from the continental crust. However, the outgoing flux of crustal He has not been directly measured at the Earth’s surface and the migration pathways are poorly understood. Here we present measurements of helium isotopes and the long-lived cosmogenic radio-isotope Kr in the deep, continental-scale Guarani aquifer in Brazil and show that crustal He reaches the atmosphere primarily by the surficial discharge of deep groundwater. We estimate that He in Guarani groundwater discharge accounts for about 20% of the assumed global flux from continental crust, and that other large aquifers may account for about 33%. Old groundwater ages suggest that He in the Guarani aquifer accumulates over half- to one-millionyear timescales. We conclude that He degassing from the continents is regulated by groundwater discharge, rather than episodic tectonic events, and suggest that the assumed steady state between crustal production and degassing of He, and its resulting atmospheric residence time, should be re-examined.