Assessing the Transport of Pharmaceutical Compounds in a Layered Aquifer Discharging to a Stream

A groundwater plume containing high concentrations of pharmaceutical compounds, mainly sulfonamides, barbiturates, and ethyl urethane, in addition to chlorinated ethenes and benzene was investigated. The contamination originating from a former pharmaceutical industry discharges into a multilayered aquifer system and a downgradient stream. In this study, geological and hydrogeological data were integrated into a numerical flow model to examine identified trends using statistical approaches, including principal component analysis and hierarchal cluster analysis. A joint interpretation of the groundwater flow paths and contaminant concentrations in the different compartments (i.e., groundwater and hyporheic zone) provided insight on the transport processes of the different contaminant plumes to the stream. The analysis of historical groundwater concentrations of pharmaceutical compounds at the site suggested these compounds are slowly degrading. The pharmaceutical compounds migrate in both a deep semiconfined aquifer, as well as in the shallow unconfined aquifer, and enter the stream along a 2-km stretch. This contrasted with the chlorinated ethenes, which mainly discharge to the stream as a focused plume from the unconfined aquifer. The integrated approach developed here, combining groundwater flow modeling and statistical analyses of the contaminant concentration data collected in groundwater and the hyporheic zone, lead to an improved understanding of the observed distribution of contaminants in the unconfined and semiconfined aquifers, and thus to their discharge to the stream. This approach is particularly relevant for large and long-lasting contaminant sources and plumes, such as abandoned landfills and industrial production sites, where field investigations may be very expensive.     

The use of regional and alluvial groundwater by riparian trees in the wet-dry tropics of northern Australia

Riparian trees play a critical role in the ecological function of rivers, yet are threatened by anthropogenic change to the hydrological cycle. Identifying the sources of water used by riparian trees can inform sustainable water policy. We used isotopic analysis complemented by measurements of plant water relations to assess water sources for riparian trees at two sites with contrasting hydrogeological processes; one with an alluvial aquifer overlaying an aquitard, and one where fault-induced preferential pathways in the aquitard allowed the flow of deeper, older groundwater from a regional aquifer to the alluvium. At both sites, plant water potential, stomatal conductance, and plant water isotope composition in the xylem sap of riparian trees were collected from two landscape positions, the riverbank and floodplain. We used a Bayesian mixing model (MixSIAR) to assess differences in the proportion of water sources for sites and landscape positions. We found that xylem water isotope values differed between the two sites in line with their hydrogeological characteristics, with trees at the regional aquifer site using water sourced from the regional groundwater and trees at the site with only an alluvial aquifer present using a mixture of water sources, with no dominant source identified. Higher plant predawn water potential values at the regional site indicated greater water availability and support the inference that plants were using more groundwater at the regional site compared to the alluvial site. Trees closer to the river had higher isotope values, indicative of surficial water sources i.e. shallow soil water and river water. Our findings show that the water sources used by riparian trees reflect local hydrogeology and resource availability. Water managers should identify and protect plant water sources to ensure maintenance of riparian trees.     

Numerical modelling of stream–aquifer interaction: Quantifying the impact of transient streambed permeability and aquifer heterogeneity

Stream–aquifer interaction plays a vital role in the water cycle, and a proper study of this interaction is needed for understanding groundwater recharge, contaminants migration, and for managing surface water and groundwater resources. A model‐based investigation of a field experiment in a riparian zone of the Schwarzbach river, a tributary of the Rhine River in Germany, was conducted to understand stream–aquifer interaction under alternative gaining and losing streamflow conditions. An equivalent streambed permeability, estimated by inverting aquifer responses to flood waves, shows that streambed permeability increased during infiltration of stream water to aquifer and decreased during exfiltration. Aquifer permeability realizations generated by multiple‐point geostatistics exhibit a high degree of heterogeneity and anisotropy. A coupled surface water groundwater flow model was developed incorporating the time‐varying streambed permeability and heterogeneous aquifer permeability realizations. The model was able to reproduce varying pressure heads at two observation wells near the stream over a period of 55 days. A Monte Carlo analysis was also carried out to simulate groundwater flow, its age distribution, and the release of a hypothetical wastewater plume into the aquifer from the stream. Results of this uncertainty analysis suggest (a) stream–aquifer exchange flux during the infiltration periods was constrained by aquifer permeability; (b) during exfiltration, this flux was constrained by the reduced streambed permeability; (c) the effect of temporally variable streambed permeability and aquifer heterogeneity were found important to improve the accurate capture of the uncertainty; and (d) probabilistic infiltration paths in the aquifer reveal that such pathways and the associated prediction of the extent of the contaminant plume are highly dependent on aquifer heterogeneity.     

Transport of pharmaceutically active compounds in saturated laboratory columns

Occurrences of pharmaceutically active compounds in surface water and sewage water have been widely reported. Investigations show the presence of several classes of pharmaceuticals such as antirheumatics (e.g., diclofenac), analgesics (e.g., propyphenazone), and blood lipid regulators (clofibric acid), even in ground water. Compared to their occurrences in surface water, however, the reported incidences of drugs in ground water are much rarer. This may be due to the input, but also to transport processes and degradation in the aquifer. In field studies investigating ground water sampled at a bank infiltration site at Lake Tegel, Berlin, Germany, clofibric acid was found at concentrations up to 290 ng/L, and propyphenazone up to 250 ng/L, whereas concentrations of diclofenac were around the detection limit. The aim of this study was to investigate the ground water transport behavior of the pharmaceuticals clofibric acid, propyphenazone, and diclofenac with a laboratory soil column experiment. Results show that clofibric acid exhibits no degradation and almost no retardation (Rf = 1.1). Diclofenac (Rf = 2.0) and propyphenazone (Rf = 1.6) are retarded, whereas significant degradation was not observed for both pharmaceuticals under the prevailing conditions in the soil column. We conclude that the concentration distribution of the pharmaceuticals at the bank filtration site at Lake Tegel is controlled by sorption, desorption, and input variation, rather than by degradation.     

Anthropogenic Impacts on a Bedrock Aquifer at the Village Scale

This study focuses on assessing groundwater potability in a highly complex and heterogeneous fractured bedrock aquifer having variable overburden cover. Eight monitoring wells were installed in a privately serviced lakeside village, and groundwater was routinely sampled over a 2‐year timeframe for concentration analysis of nitrate, fecal indicator bacteria, stable isotopes, and a total of 41 pharmaceutical compounds. While pollutant concentrations remained low throughout the study, the presence of fecal indicator bacteria and pharmaceuticals was noted at least once (but not always consistently) in most sampling intervals. An interpretation based on the integration of chemical, bacterial, and site characterization datasets suggests that: (1) the fracture network is complex and heterogeneous with limited vertical connectivity; (2) existing pathways are sufficient for the quick and widespread migration of surface contaminants to depth; (3) anthropogenic contaminants from both septic systems and agriculture are likely sourced in the surrounding uplands where overburden is thin; and (4) fecal contamination, as observed over the long term, is ubiquitous at the village scale. Groundwater quality is continually changing in this hydrogeologic environment and the determination of potability on the larger scale is not likely to be adequately captured with infrequent domestic well sampling (i.e., voluntary annual sampling by homeowners).     

Effectiveness of Monitored Natural Attenuation (MNA) as a Groundwater Remedy for Arsenic in Phosphatic Wastes

A study was conducted to evaluate monitored natural attenuation (MNA) as a remedy for arsenic in groundwater at a former phosphate mining and manufacturing facility. The mineralogy, speciation, and lability of arsenic in phosphatic wastes present in soils were characterized using sequential extraction procedures, leaching experiments, batch adsorption tests, and microchemical speciation analysis. A PHREEQC-based reactive transport model was also parameterized using these laboratory results, and it was used to evaluate the importance of identified attenuation mechanisms on arsenic concentrations along a vertical flow path from a shallow, alluvial aquifer to the underlying Floridan aquifer. Arsenic was found to occur in several chemical forms in phosphatic wastes, including unstable sulfide minerals, adsorbed surface complexes, and relatively insoluble phosphate and oxide minerals. Most arsenic was associated with stable minerals. The reactive transport model predicted that historical leaching of solid-phase waste materials in soils would not have generated enough arsenic to explain the concentrations observed in downgradient groundwater; instead, the source of arsenic to groundwater was likely acidic and saline process water that infiltrated though unlined ponds and ditches during historical manufacturing operations. A key factor affecting the long-term effectiveness of natural attenuation of arsenic in groundwater is the occurrence and stability of iron oxyhydroxides in aquifer sediments. According to laboratory and reactive transport model results, sufficient levels were found to be present at the site to effectively limit arsenic migration at concentrations exceeding drinking water standards in the future in the Floridan aquifer. This study presents the geochemical evaluations that are needed to satisfy EPA guidelines on determining whether or not MNA is an acceptable remedy for a site. It specifically details the characterization and modeling that were used to demonstrate effectiveness at a site where MNA was ultimately selected as the remedy for arsenic in groundwater.     

A geomorphology‐based integrated stream–aquifer interaction model for semi‐gauged catchments

Many of the existing stream–aquifer interaction models available in the literature are very complex with limited applicability in semi‐gauged and ungauged catchments. In this study, to estimate the influent and effluent subsurface water fluxes under limited geo‐hydrometeorological data availability conditions, a simple stream–aquifer interaction model, namely, the variable parameter McCarthy–Muskingum (VPMM) hillslope‐storage Boussinesq (hsB) model, has been developed. This novel model couples the VPMM streamflow transport with the hsB groundwater flow transport modules in online mode. In this integrated model, the surface water–groundwater flux exchange process is modelled by the Darcian approach with the variable hydraulic heads between the river stage and groundwater table accounting for the rainfall forcing. Considering the exchange fluxes in the hyporheic zone and lateral overland flow contribution, this approach is field tested in a typical 48‐km stretch of the Brahmani River in eastern India to simulate the streamflow and its depth with the minimum Nash–Sutcliffe efficiency of 94% and 88%; the maximum root mean square error of 134 m³/s and 0.35 m; and the minimum index of agreement of 98% and 97%, respectively. This modelling approach could be very well utilized in data‐scarce world‐river basins to estimate the stream–aquifer exchange flux due to rainfall forcings.     

Modelling the impact of extensive irrigation on the groundwater resources

Drastic groundwater resource depletion due to excessive extraction for irrigation is a major concern in many parts of India. In this study, an attempt was made to simulate the groundwater scenario of the catchment using ArcSWAT. Due to the restriction on the maximum initial storage, the deep aquifer component in ArcSWAT was found to be insufficient to represent the excessive groundwater depletion scenario. Hence, a separate water balance model was used for simulating the deep aquifer water table. This approach is demonstrated through a case study for the Malaprabha catchment in India. Multi‐site rainfall data was used to represent the spatial variation in the catchment climatology. Model parameters were calibrated using observed monthly stream flow data. Groundwater table simulation was validated using the qualitative information available from the field. The stream flow was found to be well simulated in the model. The simulated groundwater table fluctuation is also matching reasonably well with the field observations. From the model simulations, deep aquifer water table fluctuation was found very severe in the semi‐arid lower parts of the catchment, with some areas showing around 60 m depletion over a period of eight years. Copyright © 2012 John Wiley & Sons, Ltd.     

Groundwater contributions of flow and nitrogen in a headwater agricultural watershed

Nonpoint sources of nitrogen (N) and other nutrients are a major source of water pollution within the Chesapeake Bay watershed and other basins around the world. Human activities associated with agricultural practices can account for a large percentage of N loadings delivered to streams and rivers. This work aims to improve understanding of N transport from groundwater to surface waters, quantifying the principal hydrological processes driving water and N fluxes into and out of a headwater agricultural stream reach. The study site is a 175-m stream reach in a heavily cultivated 40-ha watershed in east-central Pennsylvania. This subwatershed is underlain by fractured shale bedrock, and receives most of its baseflow from groundwater, either by diffuse matrix discharge through the streambed or by localized discharge through riparian seeps. Samples of stream, seep, and shallow groundwater were collected approximately monthly under steady hydrologic conditions in 2017. Calculated matrix flow from hydraulic head and conductivity measurements paired with differential stream gauging was used to solve for the riparian seep flux using a mass balance approach. Riparian seep fluxes ranged from 45 to 217 m³/d, transporting 0.6–4.2 kg N d⁻¹ of nitrate-N from the fractured bedrock aquifer to the stream. Hydrochemical data suggest that the stream is mainly disconnected from the underlying aquifer and that seeps supply essentially all water and N to the system. Seeps are likely sourced with N in nearby agricultural fields and accelerated through the system with shorter residence times than shallow groundwater. Water isotope data reinforced this notion. This study underscores the importance of agriculture as a source of N to ground and surface waters. Identifying source areas that are causing groundwater enrichment of N and seep areas where N discharges to streams is beneficial for developing N pollution mitigation strategies and implementing management practices that aim to reduce nutrient loads to the Chesapeake Bay.     

Reducing monitoring gaps at the aquifer–river interface by modelling groundwater–surface water exchange flow patterns

This study investigates spatial patterns and temporal dynamics of aquifer–river exchange flow at a reach of the River Leith, UK. Observations of sub‐channel vertical hydraulic gradients at the field site indicate the dominance of groundwater up‐welling into the river and the absence of groundwater recharge from surface water. However, observed hydraulic heads do not provide information on potential surface water infiltration into the top 0–15 cm of the streambed as these depths are not covered by the existing experimental infrastructure. In order to evaluate whether surface water infiltration is likely to occur outside the ‘window of detection’, i.e. the shallow streambed, a numerical groundwater model is used to simulate hydrological exchanges between the aquifer and the river. Transient simulations of the successfully validated model (Nash and Sutcliff efficiency of 0·91) suggest that surface water infiltration is marginal and that the possibility of significant volumes of surface water infiltrating into non‐monitored shallow streambed sediments can be excluded for the simulation period. Furthermore, the simulation results show that with increasing head differences between river and aquifer towards the end of the simulation period, the impact of streambed topography and hydraulic conductivity on spatial patterns of exchange flow rates decreases. A set of peak flow scenarios with altered groundwater‐surface water head gradients is simulated in order to quantify the potential for surface water infiltration during characteristic winter flow conditions following the observation period. The results indicate that, particularly at the beginning of peak flow conditions, head gradients are likely to cause substantial increase in surface water infiltration into the streambed. The study highlights the potential for the improvement of process understanding of hyporheic exchange flow patterns at the stream reach scale by simulating aquifer‐river exchange fluxes with a standard numerical groundwater model and a simple but robust model structure and parameterization. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

Isolating the impacts of anthropogenic water use within the hydrological regime of north India

The effects of anthropogenic water use play a significant role in determining the hydrological cycle of north India. This paper explores anthropogenic impacts within the region's hydrological regime by explicitly including observed human water use behaviour, irrigation infrastructure and the natural environment in the CHANSE (Coupled Human And Natural Systems Environment) socio-hydrological modelling framework. The model is constrained by observed qualitative and quantitative information collected in the study area, along with climate and socio-economic variables from additional sources. Four separate scenarios, including business as usual (BAU, representing observed irrigation practices), groundwater irrigation only (where the influence of the canal network is removed), canal irrigation only (where all irrigation water is supplied by diverted surface water) and rainfed only (where all human interventions are removed) are used. Under BAU conditions the modelling framework closely matched observed groundwater levels. Following the removal of the canal network, which forces farmers to rely completely on groundwater for irrigation, water levels decrease, while under a canal-only scenario flooding occurs. Under the rainfed-only scenario, groundwater levels similar to current business-as-usual conditions are observed, despite much larger volumes of recharge and discharge entering and leaving the system under BAU practices. While groundwater abstraction alone may lead to aquifer depletion, the conjunctive use of surface and groundwater resources, which includes unintended contributions of canal leakage, create conditions similar to those where no human interventions are present. Here, the importance of suitable water management practices, in maintaining sustainable water resources, is shown. This may include augmenting groundwater resources through managed aquifer recharge and reducing the impacts on aquifer resources through occasional canal water use where possible. The importance of optimal water management practices that highlight trade-offs between environmental impact and human wellbeing are shown, providing useful information for policy makers, water managers and users. © 2019 John Wiley & Sons, Ltd.     

Characterising groundwater–surface water interactions in idealised ephemeral stream systems

Transmission losses from the beds of ephemeral streams are thought to be a widespread mechanism of groundwater recharge in arid and semi-arid regions and support a range of dryland hydro-ecology. Dryland areas cover ~40% of the Earth's land surface and groundwater resources are often the main source of freshwater. It is commonly assumed that where an unsaturated zone exists beneath a stream, the interaction between surface water and groundwater is unidirectional and that groundwater does not exert a significant feedback on transmission losses. To test this assumption, we conducted a series of numerical model experiments using idealised two-dimensional channel-transects to assess the sensitivity and degree of interaction between surface and groundwater for typical dryland ephemeral stream geometries, hydraulic properties and flow regimes. We broaden the use of the term ‘stream–aquifer interactions’ to refer not just to fluxes and water exchange but also to include the ways in which the stream and aquifer have a hydraulic effect on one another. Our results indicate that deep water tables, less frequent streamflow events and/or highly permeable sediments tend to result in limited bi-directional hydraulic interaction between the stream and the underlying groundwater which, in turn, results in high amounts of infiltration. With shallower initial depth to the water table, higher streamflow frequency and/or lower bed permeability, greater ‘negative’ hydraulic feedback from the groundwater occurs which in turn results in lower amounts of infiltration. Streambed losses eventually reach a constant rate as initial water table depths increase, but only at depths of 10s of metres in some of the cases studied. Our results highlight that bi-directional stream–aquifer hydraulic interactions in ephemeral streams may be more widespread than is commonly assumed. We conclude that groundwater and surface water should be considered as connected systems for water resource management unless there is clear evidence to the contrary.     

Coupled flow and salinity transport modelling in semi-arid environments: The Shashe River Valley,Botswana

Numerical groundwater modelling is used as the base for sound aquifer system analysis and water resources assessment. In many cases, particularly in semi-arid and arid regions, groundwater flow is intricately linked to salinity transport. A case in point is the Shashe River Valley in Botswana. A freshwater aquifer located around an ephemeral stream is depleted by the combined effect of transpiration and pumping. Quantitative system analysis reveals that the amount of water taken by transpiration is far more than the quantities pumped for water supply. Furthermore, the salinity distribution in and around Shashe River Valley as well as its temporal dynamics can be satisfactorily reproduced if the transpiration is modelled as a function of groundwater salinity. The location and dynamics of the saltwater–freshwater interface are highly sensitive to the parameterization of evaporative and transpirative salt enrichment. An existing numerical code for coupled flow/transport simulations (SEAWAT) was adapted to this situation. Model results were checked against a large set of field data including water levels, water chemistry, isotope data and ground and airborne geophysical data. The resulting groundwater model was able to reproduce the long-term development of the freshwater lens located in Shashe River Valley as well as the decline in piezometric heads observed over the last decade. Furthermore, the old age of the saline water surrounding the central freshwater lens could be explained.     

Innovative Contaminant Mass Flux Monitoring in an Aquifer Subject to Tidal Effects

Exposure from groundwater contamination to aquatic receptors residing in receiving surface water is dependent upon the rate of contaminated groundwater discharge. Characterization of groundwater fluxes is challenging, especially in coastal environments where tidal fluctuations result in transient groundwater flows towards these receptors. This can also be further complicated by the high spatial heterogeneity of subsurface deposits enhanced by anthropogenic influences such as the mixing of natural sediments and backfill materials, the presence of subsurface built structures such as sheet pile walls or even occurrence of other sources of contaminant discharge. In this study, the finite volume point dilution method (FVPDM) was successfully used to characterize highly transient groundwater flows and contaminant mass fluxes within a coastal groundwater flow system influenced by marked tides. FVPDM tests were undertaken continuously for more than 48 h at six groundwater monitoring wells, in order to evaluate groundwater flow dynamics during several tide cycles. Contaminant concentrations were measured simultaneously which allowed calculating contaminant mass fluxes. The study highlighted the importance of the aquifer heterogeneity, with groundwater fluxes ranging from 10⁻⁷ to 10⁻³ m/s. Groundwater flux monitoring enabled a significant refinement of the conceptual site model, including the fact that inversion of groundwater fluxes was not observed at high tide. Results indicated that contaminant mass fluxes were particularly higher at a specific monitoring well, by more than three orders of magnitude, than at other wells of the investigated aquifer. This study provided crucial information for optimizing further field investigations and risk mitigation measures.     

Coupled groundwater flow and heat transport simulation for estimating transient aquifer–stream exchange at the lowland River Spree (Germany)

Subsurface flow and heat transport near Freienbrink, NE Germany, was simulated in order to study groundwater–surface water exchange between a floodplains aquifer and a section of the lowland River Spree and an adjacent oxbow. Groundwater exfiltration was the dominant process, and only fast surface water level rises resulted in temporary infiltration into the aquifer. The main groundwater flow paths are identified based on a 3D groundwater flow model. To estimate mass fluxes across the aquifer–surface water interfaces, a 2D flow and heat transport modelling approach along a transect of 12 piezometers was performed. Results of steady‐state and transient water level simulations show an overall high accuracy with a Spearman coefficient ρ = 0.9996 and root mean square error (RMSE) = 0.008 m. Based on small groundwater flow velocities of about 10^?7 to 10^?6 ms^?1, mean groundwater exfiltration rates of 233 l m^?2 d^?1 are calculated. Short periods of surface water infiltration into the aquifer do not exceed 10 days, and the infiltration rates are in the same range. The heat transport was modelled with slightly less accuracy (ρ = 0.8359 and RMSE = 0.34 °C). In contrast to the predominant groundwater exfiltration, surface water temperatures determine the calculated temperatures in the upper aquifer below both surface water bodies down to 10 m during the whole simulation period. These findings emphasize prevailing of heat conduction over advection in the upper aquifer zones, which seems to be typical for lowland streams with sandy aquifer materials and low hydraulic gradients. Moreover, this study shows the potential of coupled numerical flow and heat transport modelling to understand groundwater–surface water exchange processes in detail. Copyright © 2013 John Wiley & Sons, Ltd.     

Three-dimensional flow patterns at the river–aquifer interface — a case study at the Danube

The three-dimensional groundwater flow patterns in a gravel bar at the Danube east of Vienna were investigated and are discussed in this paper. The observed groundwater level gradients are highly dynamic and respond very quickly to changes in the river water levels. A variably saturated groundwater model was calibrated to the data to describe the complex dynamics of flow in the gravel bar. The model results suggest that short-term (6–48 h) fluctuations of river water levels cause variations in the exchange flow rates from − 35 l/s to 82 l/s. The highest rates occur during brief infiltration after rapidly rising river water levels. Simulations of different scenarios indicate that riverbank clogging will decrease the exchange fluxes by up to 80%, while clogging of both riverbank and riverbed essentially stops the flow exchange. The groundwater model is also used to simulate the transport of a conservative tracer. The variation of river water levels over time is shown to increase the extent of the active river–aquifer mixing zone in the gravel bar. These dynamic factors significantly enhance the dilution of conservative tracer concentrations in this zone.     

Spatial delineation of riparian groundwater within alluvium deposit of mountainous region using Laplace equation

The interaction between surface and groundwater plays a key role in a riparian ecosystem while the size of riparian groundwater has not been typically incorporated into hydrological modelling systems. An extensive geophysical survey composed of 25 individual DC electrical resistivity profiles was conducted at the Blair–Wallis site in Wyoming. The observed resistivity images show a near‐surface aquifer interpreted as the saturated alluvium deposit along the channel, rather than the geological bedrock. Based on the electrical resistivity images, it can be inferred that only the near‐surface portion of the groundwater actively interacts with the stream flow in the mountainous and hilly watershed. This study attempted the spatial extrapolation of the measured riparian aquifer depths by means of fitting functions based on the surface topography. The analysis indicated that the boundary of the riparian aquifer well corresponds to the topographical inflexion point of the hill slope profile. It was also demonstrated that the extent of alluvium deposit, where the area of riparian aquifer is indicated, can be delineated using the slope and curvature maps in the geographic information system. Then, the parabolic and biharmonic functions were tested for the groundwater depth estimation using the developed alluvium deposit map. The proposed methodology was effective if geological diffusion processes by wind and water dominated the topography. The spatial map of the active aquifer will be useful in hydrological drought analysis because it is considered to be a main source of baseflow during dry seasons.     

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.