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.

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.

Underground pumped storage hydroelectricity (UPSH) plants using open-pit or deep mines can be used in flat regions to store the excess of electricity produced during low-demand energy periods. It is essential to consider the interaction between UPSH plants and the surrounding geological media. There has been little work on the assessment of associated groundwater flow impacts. The impacts on groundwater flow are determined numerically using a simplified numerical model which is assumed to be representative of open-pit and deep mines. The main impact consists of oscillation of the piezometric head, and its magnitude depends on the characteristics of the aquifer/geological medium, the mine and the pumping and injection intervals. If an average piezometric head is considered, it drops at early times after the start of the UPSH plant activity and then recovers progressively. The most favorable hydrogeological conditions to minimize impacts are evaluated by comparing several scenarios. The impact magnitude will be lower in geological media with low hydraulic diffusivity; however, the parameter that plays the more important role is the volume of water stored in the mine. Its variation modifies considerably the groundwater flow impacts. Finally, the problem is studied analytically and some solutions are proposed to approximate the impacts, allowing a quick screening of favorable locations for future UPSH plants.