Identification and quantification of groundwater and surface-water interactions provide important scientific insights for managing groundwater and surface-water conjunctively. This is especially relevant in semi-arid areas where groundwater is often the main source to feed river discharge and to maintain groundwater dependent ecosystems. Multiple field measurements were taken in the semi-arid Bulang sub-catchment, part of the Hailiutu River basin in northwest China, to identify and quantify groundwater and surface-water interactions. Measurements of groundwater levels and stream stages for a 1-year investigation period indicate continuous groundwater discharge to the river. Temperature measurements of stream water, streambed deposits at different depths, and groundwater confirm the upward flow of groundwater to the stream during all seasons. Results of a tracer-based hydrograph separation exercise reveal that, even during heavy rainfall events, groundwater contributes much more to the increased stream discharge than direct surface runoff. Spatially distributed groundwater seepage along the stream was estimated using mass balance equations with electrical conductivity measurements during a constant salt injection experiment. Calculated groundwater seepage rates showed surprisingly large spatial variations for a relatively homogeneous sandy aquifer.

Concentrations of dissolved organic matter (DOM) and ultraviolet/visible light absorbance decrease systematically as groundwater moves through the unsaturated zones overlying aquifers and along flowpaths within aquifers. These changes occur over distances of tens of meters (m) implying rapid removal kinetics of the chromophoric DOM that imparts color to groundwater. A one-compartment input-output model was used to derive a differential equation describing the removal of DOM from the dissolved phase due to the combined effects of biodegradation and sorption. The general solution to the equation was parameterized using a 2-year record of dissolved organic carbon (DOC) concentration changes in groundwater at a long-term observation well. Estimated rates of DOC loss were rapid and ranged from 0.093 to 0.21 micromoles per liter per day (μM d⁻¹), and rate constants for DOC removal ranged from 0.0021 to 0.011 per day (d⁻¹). Applying these removal rate constants to an advective-dispersion model illustrates substantial depletion of DOC over flow-path distances of 200 m or less and in timeframes of 2 years or less. These results explain the low to moderate DOC concentrations (20–75 μM; 0.26–1 mg L⁻¹) and ultraviolet absorption coefficient values (a ₂₅₄ < 5 m⁻¹) observed in groundwater produced from 59 wells tapping eight different aquifer systems of the United States. The nearly uniform optical clarity of groundwater, therefore, results from similarly rapid DOM-removal kinetics exhibited by geologically and hydrologically dissimilar aquifers.

Isotope hydrological studies carried out since 1987 in Haiti have proven to be effective in improving conceptual models of the principal aquifers in this Caribbean island nation. However, many of these studies have been limited and are mostly unpublished. It is, therefore, important to share the available datasets with scientists and water managers to support further studies and guide decision-making. Eleven datasets were compiled in a unique database for continued advancement and application. This dataset now contains 248 δ¹⁸O and δ²H results from water samples collected at 190 locations throughout the country. δ¹⁸O and δ²H values range, respectively, from –3.8 to –8.95‰ and –67.4 to –16.8‰ for precipitation (n = 5), 0.9 to 2.76‰ and 8.8 to 18.4‰ for lakes (n = 4), –2.61 to 4.98‰ and –26.63 to 28.2‰ for rivers (n = 11), –5.11 to –2.6‰ and –29.8 to –6.9‰ for springs (n = 82) excluding hot springs, and –5.1 to –1.37‰ and –29.4 to –2.9‰ for boreholes (n = 142). The overall presentation of the available data is complemented with an overview of the studies in which stable isotope data acquisition was performed. The global dataset is presented under three primary contexts: (1) elaboration of a local isotope meteoric water line of equation δ²H = 7.2δ¹⁸O + 8.8 supported with data predominately from the Ouest jurisdictional department, (2) spatial variability analysis of the dataset, and (3) definition of the altitude gradient (per 100 m) using freshwater springs of the Plaine du Cul-de-Sac/Massif de la Selle area of –6.67‰ for δ²H and –0.09‰ for δ¹⁸O.

When implementing remediation programs to mitigate diffuse-source contamination of aquifers, tools are required to anticipate if the measures are sufficient to meet groundwater quality objectives and, if so, in what time frame. Transfer function methods are an attractive approach, as they are easier to implement than numerical groundwater models. However, transfer function approaches as commonly applied in environmental tracer studies are limited to a homogenous input of solute across the catchment area and a unique transfer compartment. The objective of this study was to develop and test an original approach suitable for the transfer of spatially varying inputs across multiple compartments (e.g. unsaturated and saturated zone). The method makes use of a double convolution equation accounting for transfer across two compartments separately. The modified transfer function approach was applied to the Wohlenschwil aquifer (Switzerland), using a formulation of the exponential model of solute transfer for application to subareas of aquifer catchments. A minimum of information was required: (1) delimitation of the capture zone of the outlet of interest; (2) spatial distribution of historical and future pollution input within the capture zone; (3) contribution of each subarea of the recharge zone to the flow at the outlet; (4) transfer functions of the pollutant in the aquifer. A good fit to historical nitrate concentrations at the pumping well was obtained. This suggests that the modified transfer function approach is suitable to explore the effect of environmental projects on groundwater concentration trends, especially at an early screening stage.

Well ageing is mostly caused by mechanical and biogeochemical clogging processes, which affect the gravel pack, screen slots and casing. Clogging deposits increase head losses due to a constriction of the hydraulically effective area. For this study, clogging is mimicked by systematically reducing the gravel pack porosity, the screen open area and the nominal inner casing diameter. Groundwater flow velocity strongly increases close to the well, inducing inertial and turbulent flow components. Therefore, gravel pack head losses were calculated using the Forchheimer-Engelund equation, in conjunction with the Kozeny-Carman equation, which relates gravel pack porosity and hydraulic conductivity. Screen losses were assessed using the Orifice equation and turbulent casing losses with the Darcy-Weisbach equation. For the settings chosen here, a dramatic increase of head losses occurs when the clogging has reduced the effective porosity in the gravel pack by ~65%, the open area of the screen by ≥98%, and the casing diameter by ~50%. Since the latter two conditions are rarely reached in actual wells, the clogging of the gravel pack is the decisive parameter that controls well ageing. Regular monitoring of the well yield is therefore needed, since processes in the gravel pack are difficult to track directly. Unlike the deposits on the casing and in the screen slots, obstructions in the gravel pack are much more difficult to remove.

Past discussions around water-resources management and development in the River Nile basin disregard groundwater resources from the equation. There is an increasing interest around factoring the groundwater resources as an integral part of the Nile Basin water resources. This is hampered by knowledge gap regarding the groundwater resources dynamics (recharge, storage, flow, quality, surface-water/groundwater interaction) at basin scale. This report provides a comprehensive analysis of the state of surface-water/groundwater interaction from the headwater to the Nile Delta region. Piezometric and isotopic (δ¹⁸O, δ²H) evidence reveal that the Nile changes from a gaining stream in the headwater regions to mostly a loosing stream in the arid lowlands of Sudan and Egypt. Specific zones of Nile water leakage to the adjacent aquifers is mapped using the two sources of evidence. Up to 50% of the surface-water flow in the equatorial region of the Nile comes from groundwater as base flow. The evidence also shows that the natural direction and rate of surface-water/groundwater interaction is largely perturbed by human activities (diversion, dam construction) particularly downstream of the Aswan High Dam in Egypt. The decrease in discharge of the Nile River along its course is attributed to leakage to the aquifers as well as to evaporative water loss from the river channel. The surface-water/groundwater interaction occurring along the Nile River and its sensitivity to infrastructure development calls for management strategies that account groundwater as an integral part of the Nile Basin resources.

Sharp-interface (or interface) flow models with Dupuit-Forchheimer approximation are widely used to assess, to first order, an aquifer’s vulnerability to seawater intrusion (SWI) and to evaluate sustainable management options for coastal groundwater resources at the screening level. Recognising that interface flow models overestimate SWI, corrections have been proposed to account for the neglected mixing and also for the outflow through a finite gap. These corrections, however, were introduced in the context of specific studies and may not be generally applicable as proposed. The interface model is revisited, placing its corrections in the context of variable-density flow (VDF) theory, by expressing them in terms of the dimensionless parameters governing VDF in schematised (aspect ratio = thickness/length) homogeneous confined coastal aquifers: the coupling parameter (α), a Péclet number (Pe), and the dispersivities ratio (rα). Interfaces are compared to the 50%-salinity lines of VDF numerical solutions and regression equations are developed for estimating the outflow gap and for correcting the length of the interface (terminating with a blunted edge); the dispersion correction, which modifies the interface curvature, is restated with a variable exponent. The corrections for dispersion and for the interface length appear to be the most effective; an outflow gap is important only at small α values (strong advection relative to vertical flow due to density differences). These concepts are applied successfully to calculate the interface position in the lowermost confined sub-unit of the Coastal Plain aquifer of Israel, as an estimate of SWI.

The interactions between groundwater and surface water have been significantly affected by human activities in the semi-arid Hailiutu catchment, northwest China. Several methods were used to investigate the spatial and temporal interactions between groundwater and surface water. Isotopic and chemical analyses of water samples determined that groundwater discharges to the Hailiutu River, and mass balance equations were employed to estimate groundwater seepage rates along the river using chemical profiles. The hydrograph separation method was used to estimate temporal variations of groundwater discharges to the river. A numerical groundwater model was constructed to simulate groundwater discharges along the river and to analyze effects of water use in the catchment. The simulated seepage rates along the river compare reasonably well with the seepage estimates derived from a chemical profile in 2012. The impacts of human activities (river-water diversion and groundwater abstraction) on the river discharge were analyzed by calculating the differences between the simulated natural groundwater discharge and the measured river discharge. Water use associated with the Hailiutu River increased from 1986 to 1991, reached its highest level from 1992 to 2000, and decreased from 2001 onwards. The reduction of river discharge might have negative impacts on the riparian ecosystem and the water availability for downstream users. The interactions between groundwater and surface water as well as the consequences of human activities should be taken into account when implementing sustainable water resources management in the Hailiutu catchment.

VDFST-CFP (variable-density flow and solute transport–conduit flow process) is a density-dependent discrete-continuum numerical model for simulating seawater intrusion in a dual-permeability coastal karst aquifer. A previous study (Xu and Hu 2017) simulates variable-density flow only in a single conduit, and studies the parameter sensitivities only in the horizontal case (2D domain as horizontal section) by the VDFST-CFP model. This paper focuses on the density-dependent vertical case (2D domain as vertical section) with two major improvements: 1) when implementing double-conduit networks in the domain, simulated intruded plumes in the porous medium are extended in the double-conduit scenario, compared to the single-conduit system; 2) by quantifying micro-textures on the conduit wall by the Goudar-Sonnad equation and considering macro-structures as local head loss. Sensitivity analysis shows that medium hydraulic conductivity, conduit diameter and effective porosity are important parameters for simulating seawater intrusion in the discrete-continuum system. On the other hand, rougher micro-structures and additional macro-structure components on the conduit wall would reduce the distance of seawater intrusion to the conduit system, but, rarely affect salinity distribution in the matrix. Compared to the equivalent mean roughness height, the new method (with more detailed description of structure) simulates seawater intrusion slightly landward in the conduit system. The macro-structure measured by local head loss is more reasonable but needs further study on conduit flow. Xu and Hu (2017) Development of a discrete-continuum VDFST-CFP numerical model for simulating seawater intrusion to a coastal karst aquifer with a conduit system. Water Resources Research: 53, 688-711.

The Thiem equation of radial groundwater flow to a well is more than 100 years old and is still commonly used. Here, deviations caused by some of its simplifications are quantified by comparing the analytical to a numerical model that allows the implementation of more complex geometries. The assumption of horizontal flow in the Thiem equation, which necessitates uniform inflow over the entire screen length of the fully penetrating well, was found to cause deviations from actual pumping wells where the pump is placed above the screen, resulting in non-uniform inflow and additional drawdown. The same applies to partially penetrating wells, where inflow peaks and additional drawdown occur, especially when the well is screened in the lower part of the aquifer. The use of the Thiem equation in the near-field of a well should thus be restricted to situations where the screen inflow is relatively uniformly distributed, e.g. when it covers large portions of the aquifer thickness. The presence of a gravel pack and a background gradient, on the other hand, are of limited importance.