Delineating coastal groundwater discharge processes in a wetland area by means of electrical resistivity imaging, 224Ra and 222Rn

Coastal groundwater discharge (CGD) plays an important role in coastal hydrogeological systems as they are a water resource that needs to be managed, particularly in wetland areas. Despite its importance, identifying and monitoring CGD often presents physical and logistical constraints, restraining the application of more traditional submarine groundwater discharge surveying techniques. Here we investigate the capability of electrical resistivity imaging (ERI) in the Peníscola wetland (Mediterranean coast, Spain). ERI surveying made it possible to identify and delineate an ascending regional groundwater flow of thermal and Ra‐enriched groundwater converging with local flows and seawater intrusion. The continuous inputs of Ra‐rich groundwater have induced high activities of Ra isotopes and ²²²Rn into the marsh area, becoming among the highest previously reported in wetlands and coastal lagoons. Geoelectrical imaging enabled inferring focused upward discharging areas, leaking from the aquifer roof through a confining unit and culminating as spring pools nourishing the wetland system. Forward modelling over idealized subsurface configurations, borehole datasets, potentiometric records from standpipe piezometers, petrophysical analysis, and four natural and independent tracers (²²⁴Ra, ²²²Rn, temperature and salinity) permitted assessing the geoelectrical model and a derived hydrogeological pattern. The research highlights the potential of ERI to improve hydrogeological characterization of subsurface processes in complex contexts, with different converging flows. Additionally, a hydrogeological conceptual model for a groundwater‐fed coastal wetland was proposed, based on the integration of surveying datasets. Copyright © 2013 John Wiley & Sons, Ltd.     

Divergence and flow direction as indicators of subsurface heterogeneity and stage‐dependent storage in alluvial floodplains

Assuming homogeneity in alluvial aquifers is convenient, but limits our ability to accurately predict stream‐aquifer interactions. Research is needed on (i) identifying the presence of focused, as opposed to diffuse, groundwater discharge/recharge to streams and (ii) the magnitude and role of large‐scale bank and transient storage in alluvial floodplains relative to changes in stream stage. The objective of this research was to document and quantify the effect of stage‐dependent aquifer heterogeneity and bank storage relative to changes in stream stage using groundwater flow divergence and direction. Monitoring was performed in alluvial floodplains adjacent to the Barren Fork Creek and Honey Creek in northeastern Oklahoma. Based on results from subsurface electrical resistivity mapping, observation wells were installed in high and low electrical resistivity subsoils. Water levels in the wells were recorded real time using pressure transducers (August to October 2009). Divergence was used to quantify heterogeneity (i.e. variation in hydraulic conductivity, porosity, and/or aquifer thickness), and flow direction was used to assess the potential for large‐scale (100 m) bank or transient storage. Areas of localized heterogeneity appeared to act as divergence zones allowing stream water to quickly enter the groundwater system, or as flow convergence zones draining a large groundwater area. Maximum divergence or convergence occurred with maximum rates of change in flow rates or stream stage. Flow directions in the groundwater changed considerably between base and high flows, suggesting that the floodplains acted as large‐scale bank storage zones, rapidly storing and releasing water during passage of a storm hydrograph. During storm events at both sites, the average groundwater direction changed by at least 90° from the average groundwater direction during baseflow. Aquifer heterogeneity in floodplains yields hyporheic flows that are more responsive and spatially and temporally complex than would be expected compared to more common assumptions of homogeneity. Copyright © 2012 John Wiley & Sons, Ltd.     

Effect of variable‐density groundwater flow on nitrate flux to coastal waters

A two‐dimensional variable‐density groundwater flow and transport model was developed to provide a conceptual understanding of past and future conditions of nitrate (NO₃) transport and estimate groundwater nitrate flux to the Gulf of Mexico. Simulation results show that contaminant discharge to the coast decreases as the extent of saltwater intrusion increases. Other natural and/or artificial surface waters such as navigation channels may serve as major sinks for contaminant loading and act to alter expected transport pathways discharging contaminants to other areas. Concentrations of NO₃ in the saturated zone were estimated to range between 30 and 160 mg?L^?1 as NO₃. Relatively high hydraulic vertical gradients and mixing likely play a significant role in the transport processes, enhancing dilution and contaminant migration to depth. Residence times of NO₃ in the deeper aquifers vary from 100 (locally) to about 300 years through the investigated aquifer system. NO₃ mass fluxes from the shallow aquifers (0 to 5.7 × 10⁴ mg?m^?2?day^?1) were primarily directed towards the navigation channel, which intersects and captures a portion of the shallow groundwater flow/discharge. Direct NO₃ discharge to the sea (i.e. Gulf of Mexico) from the shallow aquifer was very low (0 to 9.0 × 10¹ mg · m^?2?day^?1) compared with discharge from the deeper aquifer system (0 to 8.2 × 10³ mg?m^?2?day^?1). Both model‐calibrated and radiocarbon tracer‐determined contaminant flux estimates reveal similar discharge trends, validating the use of the model for density‐dependent flow conditions. The modelling approach shows promise to evaluate contaminant and nutrient loading for similar coastal regions worldwide. Copyright © 2015 John Wiley & Sons, Ltd.     

Groundwater,salt and contaminant transport in a coastal aquifer system with a low-permeability layer in the intertidal zone

Low-permeability layer (LPL), formed by natural deposit or artificial reclamation and commonly found below the intertidal zone of coastal groundwater system, can retard the ingress of seawater and contaminants, and shorten the travel time of the land-sourced contaminant to the marine environment compared with a homogenous sandy coastal aquifer. However, there is limited understanding on how an intertidal LPL, a condition occurred in a coastal aquifer at Moreton Bay, Australia, influences the groundwater and contaminant transport across the shallow beach aquifer system. We characterized the aquifer hydrological parameters, monitored the in situ groundwater heads, and constructed a 2-D numerical model to analyses the cross-shore hydrological processes in this stratified system. The calibrated model suggests that in the lower aquifer, the inland-source fresh groundwater flowed horizontally towards the sea, upwelled along the freshwater–saltwater interface, and exited the aquifer at the shore below the LPL. Whereas in the upper aquifer, the tidally driven seawater circulation formed a barrier that prevented fresh groundwater from horizontal transport and discharge to the beach above the LPL, thereby directing its leakage to the lower aquifer. A contaminant represented by a conservative tracer was ‘released’ the upper aquifer in the model and results showed that the spreading extent of the contaminant plume, the maximum rate of contaminant discharge to the ocean, and its plume length decreased compared with a simulation case in a homogenous sandy aquifer. Sensitivity analysis was also conducted to investigate the characteristics of the LPL, including its continuity and hydraulic conductivity, which were found to vary along the beach at Moreton Bay. The result shows that with a lower hydraulic conductivity and continuous layer of LPL reduced the groundwater exchange and contaminant transport between upper and lower aquifer. The findings from the combined field and modelling investigations on the impact of an intertidal LPL on coastal aquifer systems highlight its significant implications to alter the groundwater and mass transport across the land–ocean interface.     

The impact of subsurface conceptualization on land energy fluxes

There is a significant body of work demonstrating the importance of hydrologic control on land energy feedbacks. Yet, quantitative data on aquifer conductivity can be difficult to assemble. Furthermore, how subsurface uncertainty propagates into land-surface processes is not well understood. This study analyzes the impact of aquifer characterization on land energy fluxes, using a coupled hydrology–land-surface model. Four gridded subsurface conductivity fields are developed for the Upper Klamath basin using two data sources and different levels of imposed heterogeneity. Each model is forced with the same transient, observed meteorology for 3 years prior to the final year presented here. Results are analyzed to quantify the impact of subsurface heterogeneity on groundwater surface water interactions and spatial patterns in hydrologic variables. Analysis shows that heterogeneity does not fundamentally alter the connection between groundwater and land surface processes. However, differences between scenarios impact the extent and location of the critical zone.     

Factors Influencing Variability in Groundwater Monitoring Data Sets

“Random” variability in groundwater monitoring data sets reduces the ability to identify long‐term concentration trends. This, in turn, increases the time and cost required to evaluate the effectiveness of natural attenuation and other groundwater remedies. To better understand the factors influencing variability in groundwater monitoring results, we have analyzed three large groundwater monitoring data sets. For the three data sets, the long‐term trend in contaminant concentration in each well accounted for an average of 30% to 40% of the overall variation in contaminant concentration. Understanding the causes of the remaining variability would support the development of improved groundwater monitoring methods. All three data sets show large differences in the temporal monitoring records between individual wells (e.g., coefficient of variation for monitoring results from individual wells ranges from 0.08 to 4.6) indicating that well and aquifer factors are more important contributors to variability than sample collection and analysis factors. However, the depth to groundwater (R² = 0.020) and distance between water level and screened interval (R² = 0.049) accounted for only a portion of the differences in variability between wells and other aquifer characteristics evaluated and were not correlated with the observed variability in monitoring results. Unidentified factors were apparently much more important contributors to variability than these factors. The monitoring data sets exhibited two distinct timescales for variability: Time‐independent variability that was apparent even when wells were re‐sampled within a few days and a long‐term variability likely associated with the long‐term concentration trend. The observation of time‐independent variability suggests that frequent monitoring of contaminated monitoring wells serves primarily to characterize sources of variability unrelated to the long‐term trend of primary interest.     

Determination of aquifer parameters based on measurements of tidal effects on a coastal aquifer near Beihai,China

Measurements of the tide and groundwater levels in coastal zones are of importance in determining the properties of coastal aquifers. The solution to a one‐dimensional unsteady groundwater flow model in a coastal confined aquifer with sinusoidal fluctuation of the tide shows that the tidal efficiency decreases exponentially with distance and the time lag increases linearly with distance from the coast. The aquifer property described by the ratio of storage coefficient to transmissivity is determined if the damping constant of the tidal efficiency or the slope of the time lag with distance are obtained on the basis of tidal measurements. Hourly observations of the tide and groundwater levels at 10 wells on the northern coast near Beihai, China show that with distance from the coast, tidal efficiency decreases roughly exponentially and the time lag increases roughly linearly. The estimated ratio of storage coefficient to transmissivity of the confined aquifer ranges from 1·169 × 10^?6 d m^?2 to 1·83 × 10^?7 d m^?2. For a given transmissivity of 750 m² d^?1, the storage coefficient of the aquifer is 8·7675 × 10^?4 with the tidal efficiency method and 1·3725 × 10^?4 with the time lag method. The damping constant of the tidal efficiency with distance can be defined as the tidal propagation coefficient. The value of the confined aquifer is determined as 0·0018892 m^?1. Copyright © 2007 John Wiley & Sons, Ltd.     

Characterizing the heterogeneity of karst critical zone and its hydrological function: An integrated approach

Spatial heterogeneity in the subsurface of karst environments is high, as evidenced by the multiphase porosity of carbonate rocks and complex landform features that result in marked variability of hydrological processes in space and time. This includes complex exchange of various flows (e.g., fast conduit flows and slow fracture flows) in different locations. Here, we integrate various “state‐of‐the‐art” methods to understand the structure and function of this poorly constrained critical zone environment. Geophysical, hydrometric, and tracer tools are used to characterize the hydrological functions of the cockpit karst critical zone in the small catchment of Chenqi, Guizhou Province, China. Geophysical surveys, using electrical resistivity tomography (ERT), inferred the spatial heterogeneity of permeability in the epikarst and underlying aquifer. Water tables in depression wells in valley bottom areas, as well as discharge from springs on steeper hillslopes and at the catchment outlet, showed different hydrodynamic responses to storm event rainwater recharge and hillslope flows. Tracer studies using water temperatures and stable water isotopes (δD and δ¹⁸O) could be used alongside insights into aquifer permeability from ERT surveys to explain site‐ and depth‐dependent variability in the groundwater response in terms of the degree to which “new” water from storm rainfall recharges and mixes with “old” pre‐event water in karst aquifers. This integrated approach reveals spatial structure in the karst critical zone and provides a conceptual framework of hydrological functions across spatial and temporal scales.     

Vertical Wellbore Flow Monitoring for Assessing Spatial and Temporal Flow Relationships with a Dynamic River Boundary

A useful tool for identifying the temporal and spatial ambient wellbore flow relationships near a dynamic river boundary is to monitor ambient vertical wellbore flow with an electromagnetic borehole flowmeter. This is important because the presence of the wellbore can result in significant mixing or exchange of groundwater vertically across the aquifer. Mixing or exchanging groundwater within the well-screen section can have significant impacts on the distribution of contaminants within the aquifer and adverse effects on the representativeness of groundwater samples collected from the monitoring well. Ambient monitoring data, collected from long screened wells at Hanford’s 300-Area Integrated Field Research Challenge site, located approximately 260 m from the Columbia River, demonstrate that vertical wellbore flow exhibits both a positive and inverse temporal relationship with periodic river-stage fluctuations that can change over short distances between wells. The spatial distribution of these vertical flows across the well field indicates two general regions of ambient wellbore flow behavior. The western region of the site is characterized by vertical flows that are positively related to river-stage fluctuations. In contrast, the eastern region of the site exhibits vertical flows that are inversely related to river-stage fluctuations. The cause of this opposite relationship is not completely understood; however, the positive relationships appear to be associated with high-energy Hanford formation flood deposits. These flood deposits have a well-defined northwest-southeast trend and are believed to coincide with a local paleochannel. The inverse relationships are attributed to an erosional, subsurface high in the Hanford/Ringold Formation contact between the site and the Columbia River. Under these complex hydrogeologic and hydrodynamic conditions, the behavior of ambient vertical wellbore flow in monitoring wells near a dynamic river boundary can have important implications for collecting groundwater-quality samples, for contributing to contaminant distribution within an aquifer system, and for implementing effective remediation strategies.     

Surface water–groundwater exchange dynamics in a tidal freshwater zone

In coastal rivers, tides can propagate for tens to hundreds of kilometres inland beyond the saltwater line. Yet the influence of tides on river–aquifer connectivity and solute transport in tidal freshwater zones (TFZs) is largely unknown. We estimate that along the TFZ of White Clay Creek (Delaware, USA), 11% of river water exchanges through tidal bank storage zones. Additional hyporheic processes such as flow through bedforms likely contribute even more exchange. The turnover length associated with tidal bank storage is 150 km, on the order of turnover lengths for all hyporheic exchange processes in non‐tidal rivers of similar size. Based on measurements at a transect of piezometers located 17 km from the coast, tides exchange 0.36 m³ of water across the banks and 0.86 m³ across the bed per unit river length. Exchange fluxes range from ?1.66 to 2.26 m day^?1 across the bank and ?0.84 to 1.88 m day^?1 across the bed. During rising tide, river water infiltrates into the riparian aquifer, and the downstream transport rate in the channel is low. During falling tide, stored groundwater is released to the river, and the downstream transport rate in the channel increases. Tidal bank storage zones may remove nutrients or other contaminants from river water and attenuate nutrient loads to coasts. Alternating expansion and contraction of aerobic zones in the riparian aquifer likely influence contaminant removal along flow paths. A clear need exists to understand contaminant removal and other ecosystem services in TFZs and adopt best management practices to promote these ecosystem services. Copyright © 2015 John Wiley & Sons, Ltd.     

Conceptualizing a multi-layered shingle aquifer model based on volcanic stratigraphy and water inflow to lava caves in Jeju Island,Korea

Volcanic aquifers supply a substantial portion of water resources in many parts of the world, including islands, and their productivity depends strongly on volcanic stratigraphy, which exhibits considerable heterogeneity. We investigated water inflow to lava tube caves formed from numerous basaltic lava flows in the northeastern coastal area of Jeju Island after storm events and monitored relative inflow rates monthly over 1 year to characterize groundwater flow processes in the upper parts of volcanic aquifers, and to evaluate the applicability of the previous hydrogeological models proposed for the island. Considerable water inflow arose shortly after storms from exposed palaeosol layers on the walls of the caves. The monthly monitoring results showed that wall inflow associated with these palaeosol layers is substantial. In both cases, discharge from ceiling drips was much less and more temporally variable compared to wall inflow discharge. Water flowing into the caves was rapidly drained through the floor at all monitoring sites. The lateral extent of the palaeosol layers was identified using drill core logs near the cave and outcrops in the coastal area. Based on these results, we inferred that multiple perched aquifers are formed by low-permeability palaeosol layers between lava flows, which are connected by vertical flows at discontinuities in the palaeosol layer, eventually reaching the basal aquifer. This study revealed the water inflow processes observed in lava tube caves constrained by palaeosol layers, and established a hydrogeological conceptual model incorporating multiple perched aquifers in both coastal and mountainous areas associated with extensive palaeosol layers formed during volcanic hiatuses. This finding would help elucidate recharge, groundwater flow, and contaminant transport processes in many volcanic aquifers that are not adequately represented by the previous models, and contribute to better management of groundwater in those areas.     

Tidal influence on BTEX biodegradation in sandy coastal aquifers

A numerical study was conducted to investigate the influence of tides on the fate of terrestrially derived BTEX discharging through an unconfined aquifer to coastal waters. Previous studies have revealed that tide-induced seawater circulations create an active salt–freshwater mixing zone in the near-shore aquifer and alter the specific subsurface pathway for contaminants discharging to the coastal environment. Here the coupled density-dependent flow and multi-species reactive transport code PHWAT was used to examine the impact of these tidal effects on the aerobic biodegradation of BTEX released in a coastal aquifer and its subsequent loading to coastal waters. Simulations indicated that tides significantly enhance BTEX attenuation in the near-shore aquifer. They also reduce the rate of chemical transfer from the aquifer to the ocean and exit concentrations at the beach face. For the base case consisting of toluene transport and biodegradation, 79% of toluene initially released in the aquifer was attenuated prior to discharge with tides present, compared to only 1.8% for the non-tidal case. The magnitude of tidal forcing relative to the fresh groundwater flow rate was shown to influence significantly the extent of biodegradation as it controls the intensity of salt–freshwater mixing, period of exposure of the contaminant to the mixing zone and rate of oxygen delivery to the aquifer. The oxygen available for biodegradation also depends on the rate at which oxygen is consumed by natural processes such as organic matter decomposition. While simulations conducted with heterogeneous conductivity fields highlighted the uncertainties associated with predicting contaminant loadings, the study revealed overall that BTEX may undergo significant attenuation in tidally influenced aquifers prior to discharge.     

Persistence of a Groundwater Contaminant Plume after Hydraulic Source Containment at a Chlorinated‐Solvent Contaminated Site

The objective of this study was to characterize the behavior of a groundwater contaminant (trichloroethene, TCE) plume after implementation of a source‐containment operation at a site in Arizona. The plume resides in a quasi‐three‐layer system comprising a sand/gravel unit bounded on the top and bottom by relatively thick silty clayey layers. The system was monitored for 60 months beginning at start‐up in 2007 to measure the change in contaminant concentrations within the plume, the change in plume area, the mass of the contaminant removed, and the integrated contaminant mass discharge (CMD). The concentrations of TCE in groundwater pumped from the plume extraction wells have declined significantly over the course of operation, as have concentrations for groundwater sampled from 40 monitoring wells located within the plume. The total CMD associated with operation of the plume extraction wells peaked at 0.23 kg/d, decreased significantly within 1 year, and thereafter began an asymptotic decline to a current value of approximately 0.03 kg/d. Despite an 87% reduction in contaminant mass and a comparable 87% reduction in CMD for the plume, the spatial area encompassed by the plume has decreased by only approximately 50%. This is much less than would be anticipated based on ideal flushing and mass‐removal behavior. Simulations produced with a simplified three‐dimensional (3D) numerical model matched reasonably well to the measured data. The results of the study suggest that permeability heterogeneity, back diffusion, hydraulic factors associated with the specific well field system, and residual discharge from the source zone are all contributing to the observed persistence of the plume, as well as the asymptotic behavior currently observed for mass removal and for the reduction in CMD.     

Suspended sediment discharge in subsurface flow from the head hollow of a small forested watershed,northern Japan

Surface flow and suspended sediment discharge from the head hollow of the Jozankei Experimental Watershed in Hokkaido, northern Japan, were measured to clarify the implications of subsurface hydrology for soil movement. Subsurface discharges during the extremely large storms of 1993 to 1994 were measured in a V-notch weir installed at a natural spring near the bottom of the head hollow, and shallow groundwater levels were observed in the wells excavated in the hollow. Sediment samples whose particle size range from 0·001 to 0·1 mm were manually and automatically collected at 15 to 60 min intervals, by use of 1 or 21 polyethylene bottles. Maximum concentration and flux of suspended sediment during the storms preceded the peak discharge of subsurface flow by several hours. Neither the changes in concentration (mg l⁻¹) nor flux (mg s⁻¹) of suspended sediment coincided with those in subsurface discharge (l s⁻¹). Furthermore, sediment concentration was poorly correlated with the rate of change in subsurface discharge (l s⁻²) during the rising limb of the hydrograph. Suspended sediment flux during the acceleratory limb, however, was closely correlated with the rate of change in subsurface discharge. The relationship between suspended sediment flux and rate of change in subsurface discharge were in inverse proportion to initial subsurface discharge before the storm runoff and they represented rare seasonal variation. Subsurface hydraulic erosion and transport of suspended sediment resulting from changes in rate of change in subsurface discharge actively occur during the acceleratory rising limb of the hydrograph. Accordingly, subsurface hydraulic erosion during the acceleratory rising limb of the hydrograph can be physically understood by analysing suspended sediment flux associated with rate of change in subsurface discharge and initial subsurface discharge. © 1997 John Wiley & Sons, Ltd.     

Lagged rejuvenation of groundwater indicates internal flow structures and hydrological connectivity

Large proportions of rainwater and snowmelt infiltrate into the subsurface before contributing to stream flow and stream water quality. Subsurface flow dynamics steer the transport and transformation of contaminants, carbon, weathering products and other biogeochemistry. The distribution of groundwater ages with depth is a key feature of these flow dynamics. Predicting these ages are a strong test of hypotheses about subsurface structures and time-varying processes. Chlorofluorocarbon (CFC)-based groundwater ages revealed an unexpected groundwater age stratification in a 0.47 km² forested catchment called Svartberget in northern Sweden. An overall groundwater age stratification, representative for the Svartberget site, was derived by measuring CFCs from nine different wells with depths of 2–18 m close to the stream network. Immediately below the water table, CFC-based groundwater ages of already 30 years that increased with depth were found. Using complementary groundwater flow models, we could reproduce the observed groundwater age stratification and show that the 30 year lag in rejuvenation comes from return flow of groundwater at a subsurface discharge zone that evolves along the interface between two soil types. By comparing the observed groundwater age stratification with a simple analytical approximation, we show that the observed lag in rejuvenation can be a powerful indicator of the extent and structure of the subsurface discharge zone, while the vertical gradient of the age-depth-relationship can still be used as a proxy of the overall aquifer recharge even when sampled in the discharge zone. The single age stratification profile measured in the discharge zone, close to the aquifer outlet, can reveal the main structure of the groundwater flow pattern from recharge to discharge. This groundwater flow pattern provides information on the participation of groundwater in the hydrological cycle and indicates the lower boundary of hydrological connectivity.     

A Geochemical Approach to Determine Sources and Movement of Saline Groundwater in a Coastal Aquifer

Geochemical evaluation of the sources and movement of saline groundwater in coastal aquifers can aid in the initial mapping of the subsurface when geological information is unavailable. Chloride concentrations of groundwater in a coastal aquifer near San Diego, California, range from about 57 to 39,400 mg/L. On the basis of relative proportions of major‐ions, the chemical composition is classified as Na‐Ca‐Cl‐SO₄, Na‐Cl, or Na‐Ca‐Cl type water. δ²H and δ¹⁸O values range from ?47.7‰ to ?12.8‰ and from ?7.0‰ to ?1.2‰, respectively. The isotopically depleted groundwater occurs in the deeper part of the coastal aquifer, and the isotopically enriched groundwater occurs in zones of sea water intrusion. ⁸⁷Sr/⁸⁶Sr ratios range from about 0.7050 to 0.7090, and differ between shallower and deeper flow paths in the coastal aquifer. ³H and ¹⁴C analyses indicate that most of the groundwater was recharged many thousands of years ago. The analysis of multiple chemical and isotopic tracers indicates that the sources and movement of saline groundwater in the San Diego coastal aquifer are dominated by: (1) recharge of local precipitation in relatively shallow parts of the flow system; (2) regional flow of recharge of higher‐elevation precipitation along deep flow paths that freshen a previously saline aquifer; and (3) intrusion of sea water that entered the aquifer primarily during premodern times. Two northwest‐to‐southeast trending sections show the spatial distribution of the different geochemical groups and suggest the subsurface in the coastal aquifer can be separated into two predominant hydrostratigraphic layers.     

Evaluation of long-term groundwater level data in regular monitoring wells,Barka, Sultanate of Oman

A detailed investigation was carried out to evaluate long-term groundwater level fluctuation in regular monitoring wells constructed by the Ministry of Water Resources in Barka, Sultanate of Oman. For this study, groundwater level data for 71 wells and rainfall data from six stations were collected from 1984 to 2003 and analysed. Based on long-term water level fluctuation, groundwater wells are classified into three groups. In group 1, water level shows a long-term cyclic trend without yearly fluctuation whereas in group 2 the water level declined continuously until 1995 followed by a constant water level. In group 3, water level decreases continuously throughout the study periods with rapid annual cyclic variation. Group 1 wells show high water-level fluctuations (5 to 10 m) and seem to be regulated by discharge (lateral flow) from this aquifer and recharge from the adjacent Jabal Akhdar mountainous region. Constant trend in water level after 1995 in group 2 wells illustrates the advancement of saline–fresh water interface to the inland due to heavy pumping which is justified by higher electrical conductivity and Cl/HCO₃ ratio. In group 3 wells, the water level dropped continuously due to overabstraction by agricultural farms and human settlements. In addition, wells existing near the recharge dams express the influences of recharge dams and rainfall, and exhibit high water-level fluctuations during heavy rainfall periods. The long-term regional variation indicates that water level drops continuously in the coastal and central parts of the study region. Linear regression analysis revealed that the decline in water level is 0·3–0·4 m year⁻¹ near the coastal and central parts of the study area and is almost constant in the remaining area. We conclude that the contribution of man-made activities on groundwater level is well compared with natural factors. Copyright © 2007 John Wiley & Sons, Ltd.     

Sulfur Hexafluoride and Potassium Bromide as Groundwater Tracers for Managed Aquifer Recharge

Sulfur hexafluoride (SF₆) is an established tracer for use in managed aquifer recharge projects. SF₆ exsolves from groundwater when it encounters trapped air according to Henry's law. This results in its retardation relative to groundwater flow, which can help determine porous media saturation and flow dynamics. SF₆ and the conservative, nonpartitioning tracer, bromide (Br⁻ added as KBr), were introduced to recharge water infiltrated into stacked glacial aquifers in Thurston County, Washington, providing the opportunity to observe SF₆ partitioning. Br⁻, which is assumed to travel at the same velocity as the groundwater, precedes SF₆ at most monitoring wells (MWs). Average groundwater velocity in the unconfined aquifer in the study area ranges from 3.9 to 40 m/d, except in the southwestern corner where it is slower. SF₆ in the shallow aquifer exhibits an average retardation factor of 2.5 ± 3.8, suggesting an air-to-water ratio on the order of 10⁻³ to 10⁻² in the pore space. Notable differences in tracer arrival times at adjacent wells indicate very heterogeneous conductivity. One MW exhibits double peaks in concentrations of both tracers with different degrees of retardation for the first and second peaks. This suggests multiple flowpaths to the well with variable saturation. The confining layer between the upper two aquifers appears to allow intermittent connection between aquifers but serves as an aquitard in most areas. This study demonstrates the utility of SF₆ partitioning for evaluating hydrologic conditions at prospective recharge sites.     

Localizing and quantifying groundwater-surface water interactions at different scales: A tracer approach at the River Moselle,Germany

Groundwater-surface water interactions (GSI) connect rivers and streams with riparian areas and the adjacent aquifer. Although these interactions exert a substantial control of quantity and quality of both groundwater and surface water, knowledge on GSI along rivers at the regional scale, particularly for inland waterways, is still limited. We investigated GSI along the river Moselle, an important federal inland waterway in Germany, by using radon and tritium to identify gaining (water flux from the aquifer to the surface water) and losing (water flux from the surface water to the aquifer) stream conditions, respectively. Gaining stream conditions were identified by continuously measuring radon along the river during boat surveys with a high spatial resolution (every 2 km) during intermediate (October 2020) and near low flow conditions (August/September 2021). The tritium concentrations in surface water and groundwater and the resulting tritium inventories were used to characterize losing stream conditions Monthly tritium inventories from 2017 to 2022 revealed a mean loss for the whole period of 20.3 % and a mean gain of 21.8%. Both were probably triggered by a combination of losing stream conditions and flood-induced mass transfer of water from the aquifer back into the river as well as discharge fluctuations. At the investigated site Lehmen there were direct indications of an influence of surface water due to elevated tritium concentrations in the groundwater (up to 13.3 Bq L⁻¹). Using radon mass balance modelling, good agreements of simulated versus measured radon data with respect to two groundwater end-member scenarios were obtained during intermediate flow (Spearman's ρ: 0.97 and 0.99; MAE: 10.1 and 3.4 Bq L⁻¹) and near low flow (Spearman's ρ: 0.97 and 0.99; MAE: 11 and 6.5 Bq L⁻¹). Considerable groundwater inflow was limited to the meander of Detzem, where cumulated groundwater inflow of about 19 m³ s⁻¹ (9.5% of total discharge) and 4.2 m³ s⁻¹ (3.8% of total discharge) was simulated during intermediate and near low flow, respectively. However, the groundwater inflow was relatively low compared to alpine streams, for example. The study will help to better identify and quantify GSI at the regional scale and provide methodological guidance for future studies focusing on inland waterways.