Arsenic is a common trace element in groundwater and its fate and transport are controlled by combination of (i) natural processes, including redox conditions, salinity and pH, (ii) sedimentary and geochemical environment, and (iii) anthropogenic influences such as groundwater extraction, managed aquifer recharge (MAR), Aquifer Thermal Energy Storage (ATES), and pollution. We investigated the relative influence of these processes by presenting 10 cases from The Netherlands. Our review showed that the primary controlling factor for arsenic mobility in natural coastal dune systems is the redox state of groundwater, with concentrations between 2 and 10 μg/L. Strongly reduced greensands (containing glauconite and, more importantly, associated minerals) exhibited elevated As concentrations, with concentrations up to 40 μg/L. Groundwater systems modified by MAR or those that are influenced by nitrate pollution (NP) showed elevated As concentrations (20–110 μg/L), as a result of either pyrite oxidation (MAR, NP) or reductive dissolution of iron(hydr)oxides (MAR). Increasing temperature at ATES systems may cause mobilization of As at temperatures beyond 25 °C. The highest As concentrations were observed at sites where muddy sediments were recently deposited in surface water bodies (200–820 μg/L), for example in dammed Rhine River tributaries and sand pit lakes south of the city of Amsterdam. The reduction of arsenate to arsenite and competitive desorption during intrusion of polluted water also form important As mobilizing processes. The data for the Netherlands show that high CH₄ and NH₄ concentrations may form a risk indicator of elevated As levels in some fresh groundwater systems.

Well clogging was studied at an aquifer storage transfer and recovery (ASTR) site used to secure freshwater supply for a flower bulb farm. Tile drainage water (TDW) was collected from a 10-ha parcel, stored in a sandy brackish coastal aquifer via well injection in wet periods, and reused during dry periods. This ASTR application has been susceptible to clogging, as the TDW composition largely exceeded most clogging mitigation guidelines. TDW pretreatment by sand filtration did not cause substantial clogging at a smaller ASR site (2 ha) at the same farm. In the current (10 ha) system, sand filtration was substituted by 40-μm disc filters to lower costs (by 10,000–30,000 Euro) and reduce space (by 50–100 m²). This measure treated TDW insufficiently and injection wells rapidly clogged. Chemical, biological, and physical clogging occurred, as observed from elemental, organic carbon, 16S rRNA, and grain-size distribution analyses of the clogging material. Physical clogging by particles was the main cause, based on the strong relation between injected turbidity load and normalized well injectivity. Periodical backflushing of injection wells improved operation, although the disc filters clogged when the turbidity increased (up to 165 NTU) during a severe rainfall event (44 mm in 3 days). Automated periodical backflushing, together with regulating the maximum turbidity (<20 NTU) of the TDW, protected ASTR operation, but reduced the injected TDW volume by ~20–25%. The studied clogging-prevention measures collectively are only viable as an alternative for sand filtration when the injected volume remains sufficient to secure the farmer’s needs for irrigation.

Larger well diameters allow higher groundwater abstraction rates. But particularly for the construction of wells at greater depth, it may be more cost-efficient to only expand the borehole in the target aquifer. However, current drilling techniques for unconsolidated formations are limited by their expansion factors (<2) and diameters (<1000 mm). Therefore, we developed a new technique aiming to expand borehole diameters at depth in a controlled manner using a low-pressure water jet perpendicular to the drilling direction and extendable by means of a pivoting arm. During a first field test, the borehole diameter was expanded 2.6-fold from 600 to 1570 mm at a depth of 53.5 to 68 m and equipped with a well screen to create an expanded diameter gravel well (EDGW). In keeping with the larger diameter, the volume flux per m screen length was two times higher than conventional 860 mm diameter wells at the site in the subsequent 3 year production period. Although borehole clogging was slower on a volumetric basis and similar when normalized to borehole wall area, rehabilitation of particle clogging at the borehole wall was more challenging due to the thickness of the gravel pack. While jetting the entire borehole wall before backfilling holds promise to remove filter cake and thus limit particle clogging, we found that a second borehole (expanded 4.1-fold to 2460 mm) collapsed during jetting. Overall, the EDGW technique has potential to make the use of deeper unconsolidated aquifers economically (more) feasible, although further understanding of the borehole stability and rehabilitation is required to assess its wider applicability.

An aquifer storage transfer and recovery (ASTR) system was studied in which tile drainage water (TDW) was injected with relatively high NO₃ (about 14 mg/L) concentrations originating from fertilizers. Here we present the evolution of denitrification kinetics at 6 different depths in the aquifer before, and during ASTR operation. First-order denitrification rate constants increased over time before and during the first days of ASTR operation, likely due to microbial adaptation of the native bacterial community and/or bioaugmentation of the aquifer by denitrifying bacteria present in injected TDW. Push-pull tests were performed in the native aquifer before ASTR operation. Obtained first-order denitrification rate constants were negligible (0.00–0.03 d⁻¹) at the start, but increased to 0.17–0.83 d⁻¹ after a lag-phase of about 6 days. During the first days of ASTR operation in autumn 2019, the arrival of injected TDW was studied at 2.5 m distance from the injection well. First-order denitrification rate constants increased again over time (maximum >1 d⁻¹). Three storage periods without injection were monitored in winter 2019, fall 2020, and spring 2021 during ASTR operation. First-order rate constants ranged between 0.12 and 0.61 d⁻¹. Denitrification coupled to pyrite oxidation occurred at all depths, but other oxidation processes were indicated as well. NO₃ concentration trends resembled Monod kinetics but were fitted also to a first-order decay rate model to facilitate comparison. Rate constants during the storage periods were substantially lower than during injection, probably due to a reduction in the exchange rate between aquifer solid phases and injected water during the stagnant conditions. Denitrification rate constants deviated maximally a factor 5 over time and depth for all in-situ measurement approaches after the lag-phase. The combination of these in-situ approaches enabled to obtain more detailed insights in the evolution of denitrification kinetics during AS(T)R.

Degradation of 7 common pesticides (bentazon, boscalid, chloridazon, fluopyram, flutolanil, imidacloprid, and methoxyfenozide) and 2 metabolites of chloridazon (desphenyl-chloridazon, and methyl-desphenyl-chloridazon) was studied in an anoxic and brackish sandy aquifer before and during Aquifer Storage Transfer and Recovery (ASTR) operation. Fresh tile drainage water was injected and stored for later re-use as irrigation water. We hypothesized that electron acceptors (O₂, NO₃), dissolved organic carbon (∼24.7 mg/L), nutrients (NO₃: ∼14.1 mg/L, NH₄: ∼0.13 mg/L, PO₄: ∼5.2 mg/L), and biodegrading bacteria in tile drainage water could stimulate degradation of the pesticides and metabolites (ranging between 0.013 and 10.8 μg/L) introduced in the aquifer. Pesticide degradation was studied at 6 depths in the aquifer using push-pull tests lasting ±18 days before the onset of ASTR operation. Degradation was too limited to quantify and/or could not be assessed because of the potential occurrence of pesticide retardation. Utilizing push-pull tests to obtain degradation constants should only be considered in future studies for non-retarding pesticides with relative low half-lives (here <20 days). During ASTR operation, pesticide degradation was studied at the same depths during 3 storage periods equally spread over 1.5 years of ASTR operation. Overall, trends of degradation were observed, although with relatively high half-lives of at least 53 days. Microbial adaptation of the aquifer and/or bioaugmentation by the injected biodegrading bacteria did not result in enhanced degradation during consecutive storage periods. Operational monitoring data over longer periods and distances yielded half-lives of at least 141 days. The slow degradation mostly agrees with previous studies. The injected tile drainage water composition did therefore not notably stimulate pesticide degradation. The relatively persistent behavior of the studied pesticides/metabolites implies that ASTR abstracted water will have generally high pesticide concentrations, and non-abstracted water may form a contamination risk for the surrounding native brackish groundwater.

A field injection experiment was performed in an anoxic sandy aquifer over 6 days to assess sorption characteristics of 7 commonly applied pesticides in agriculture and 2 frequently detected metabolites. Pesticide use changed considerably in the last decades, and there is insufficient knowledge of the fate of currently used pesticides in aquifers. Injected water arrival was monitored at 6 depth intervals of 1 m ranging from 11.4 to 32.2 m-below surface level with varying organic carbon contents (0.057–0.91%d.w.) to examine intra-aquifer variations in sorption. Observed pesticide concentrations were fit using a non-linear least squares routine to an advection-dispersion equation, from which retardation factors (R) were obtained. Pesticide degradation did not significantly influence the simulated R during the experiment. We observed that bentazon and cycloxydim were most mobile with R < 1.1 at all depths. Desphenyl chloridazon, methyl desphenyl chloridazon, and imidacloprid were, on average, less mobile, with maximum R of 1.5. Boscalid, chloridazon, fluopyram, and flutolanil showed a larger range of R, and R > 2.0 were observed in the shallowest part of the aquifer. Largest R were observed at the top of the aquifer and decreased with depth. K_oc values varied similarly, which indicates that sorption is not only influenced by sedimentary organic matter (SOM) content but also by its sorption reactivity. Obtained sorption parameters were substantially lower than reported in a widely used pesticide sorption database, which suggests that sorption parameters are influenced by methodological differences and variations in the sorption reactivity of SOM. The large intra-aquifer variations in pesticide sorption highlights that aquifer heterogeneity should be considered in groundwater risk assessments.

Various hydrogeochemical processes can modify the quality of river water during riverbank filtration (RBF). Identifying the subsurface processes responsible for the bank-filtered water quality is challenging, but essential for predicting water quality changes and determining the necessity of post-treatment. However, no systematic approach for this has been proposed yet. In this study, the subsurface hydrogeochemical processes that caused the high concentrations of total iron (Fe) and sulfate (SO₄²⁻) in the bank-filtered water were investigated at a pilot-scale RBF site in South Korea. For this purpose, water quality variations were monitored in both the extraction well and the adjacent river over five months. The volumetric mixing ratio between the river water and the native groundwater in the RBF well was calculated to understand the effect of mixing on the quality of water from the well and to assess the potential contribution of subsurface reactions to water quality changes. To identify the subsurface processes responsible for the evolution of Fe and SO₄²⁻ during RBF, an inverse modeling based on the chemical mass balance was conducted using the water quality data and the calculated volumetric mixing ratio. The modeling results suggest that pyrite oxidation by abundant O₂ present in an unsaturated zone could be a primary process explaining the evolution of total Fe and SO₄²⁻ during RBF at the study site. The presence of pyrite in the aquifer was indirectly supported by iron sulfate hydroxide (Fe(SO₄)(OH)) detected in oxidized aquifer sediments.

Arsenic (As) is a highly toxic element which naturally occurs in drinking water. In spite of substantial evidence on the association between many illnesses and chronic consumption of As, there is still a considerable uncertainty about the health risks due to low As concentrations in drinking water. In the Netherlands, drinking water companies aim to supply water with As concentration of <1 μg/L – a water quality goal which is tenfold more stringent than the current WHO guideline. This paper provides (i) an account on the assessed lung cancer risk for the Dutch population due to pertinent low-level As in drinking water and cost-comparison between health care provision and As removal from water, (ii) an overview of As occurrence and mobility in drinking water sources and water treatment systems in the Netherlands and (iii) insights into As removal methods that have been employed or under investigation to achieve As reduction to <1 µg/L at Dutch water treatment plants. Lowering of the average As concentration to <1μg/L in the Netherlands is shown to result in an annual benefit of 7.2–14 M€. This study has a global significance for setting drinking water As limits and provision of safe drinking water.

As part of an integrated water-cycle management strategy, City West Water (CWW) is conducting research to develop an aquifer storage recovery (ASR) scheme utilizing recycled water. In this contribution, we address the risk of well clogging based on two ASR bore pilots, each with intensive monitoring. Well clogging is a critical aspect of the strategy due to a projected high injection rate, a high clogging potential of recycled water, and a small diameter injection borehole. Microscopic and geochemical analysis of suspended solids in the injectant and backflushed water, demonstrate a significant contribution of diatoms, algae and colloidal or precipitating Fe(OH)₃, Al(OH)₃ and MnO₂. CWW is, therefore, testing additional prefiltration that includes a 20 μm spin Klin disc and 1-5 μm bag filter operating in series. In this paper, we present optimized methods to (i) detect the contribution of the injectant and aquifer particles to total suspended solids in backflushed water by hydrogeochemical analysis; and (ii) predict and reduce the risk of physical and biological clogging, by combination of the membrane filter index (MFI) method of Buik and Willemsen, a modification of the total suspended solids method of Bichara and an amendment of the exponential bacterial growth method of Huisman and Olsthoorn.

Soil passage of (pretreated) surface water to remove pathogenic microorganisms is a highly efficient process under oxic conditions, reducing microorganism concentrations about 8 log¹⁰ within tens of meters. However, under anoxic conditions, it has been shown that removal of microorganisms can be limited very much. Setback distances for adequate protection of natural groundwater may, therefore, be too short if anoxic conditions apply. Because removal of microorganisms under suboxic conditions is unknown, this research investigated removal of bacteriophage MS2 and PRD1 by soil passage under suboxic conditions at field scale. At the field location (dune area), one injection well and six monitoring wells were installed at different depths along three suboxic flow lines, where oxygen concentrations ranged from 0.4 to 1.7 mg/l and nitrate concentrations ranged from 13 to 16 mg/L. PRD1 and MS2 were injected directly at the corresponding depths and their removal in each flow line was determined. The highest bacteriophage removal was observed in the top layer, with about 9 log removal of MS2, and 7 log removal of PRD1 after 16 meters of aquifer transport. Less removal was observed at 12 m below surface, probably due to a higher groundwater velocity in this coarser grained layer. MS2 was removed more effectively than PRD1 under all conditions. Due to short travel times, inactivation of the phages was limited and the reported log removal was mainly associated with attachment of phages to the aquifer matrix. This study shows that attachment of MS2 and PRD1 is similar for oxic and suboxic sandy aquifers, and, therefore, setback distances used for sandy aquifers under oxic and suboxic conditions provide a similar level of safety. Sticking efficiency and the attachment rate coefficient, as measures for virus attachment, were evaluated as a function of the physico-chemical conditions.

In the coastal dunes of the Western Netherlands, managed aquifer recharge (MAR) is applied for drinking water supply since 1957. The MAR systems belong to the Aquifer Transfer Recovery (ATR) type, because recharge and recovery are operated without interruption. This makes these systems very vulnerable to intake interruptions, which are expected to increase in frequency and duration due to climate change. Such interruptions are problematic, because: (i) groundwater recovery from dunes needs to continue to supply fresh drinking water to the Western Netherlands; (ii) risks of salt water intrusion are high, and (iii) MAR bordering wet dune slacks with an EU Natura 2000 status cannot survive for long without MAR. In this paper, effects of intake stops are discussed and quantified. The hydrological effects consist of the decline of water tables, disappearance of flow-through dune lakes, reservoir depletion, salt water intrusion, disruption of rainwater lenses, and entrapped air hampering a rapid refill of the groundwater reservoir. Water quality effects include changes in (i) redox environment of the flushed aquifer, impacting the behavior of nutrients, calcium, sulfate and organic micro-pollutants, and (ii) the mixing ratio of water types. The main ecological impacts comprise the dying of organisms in recharge ponds and dune lakes, and a decline of biodiversity. Effects of very long intake interruptions (years) are predicted via historical observations during the long overexploitation period (1900–1957) prior to MAR. A closed form analytical solution for safe yield of a semiconfined aquifer is proposed, together with a related upconing risk index. Both also apply to the pumping from any fresh water lens without MAR. Some mitigation strategies are discussed, such as a dual intake, raising the storage capacity, earlier mud removal, and accelerated refilling of the reservoir. A magnitude scale for intake stops (MIS) is proposed.

To date, there has been no agreement on the best way to simulate saltwater intrusion (SWI) in karst aquifers. An equivalent porous medium (EPM) is usually assumed without justification of its applicability. In this paper, SWI in a poorly karstified aquifer in Lebanon is simulated in various ways and compared to measurements. Time series analysis of rainfall and aquifer response is recommended to decide whether quickflow through conduits can be safely ignored. This aids in justifying the selection of the exemplified EPM model. To examine the improvement of SWI representation when discrete features (DFs) are embedded in the model domain, the results of a coupled discrete-continuum (CDC) approach (a hybrid EPM-DF approach) are compared to the EPM model. The two approaches yielded reasonable patterns of hydraulic head and groundwater salinity, which seem trustworthy enough for management purposes. The CDC model also reproduced some local anomalous chloride patterns, being more adaptable with respect to the measurements. It improved the overall accuracy of salinity predictions at wells and better represented the fresh–brackish water interface. Therefore, the CDC approach can be beneficial in modeling SWI in poorly karstified aquifers, and should be compared with the results of the EPM method to decide whether the differences in the outcome at local scale warrant its (more complicated) application. The simulation utilized the SEAWAT code since it is density dependent and public domain, and it enjoys widespread application. Including DFs necessitated manual handling because the selected code has no built-in option for such features.

This study demonstrates groundwater quality differences between a limestone and a dolomitic limestone, (sub)oxic coastal aquifer in the Eastern Mediterranean (Lebanon), with and without ongoing moderate salinization since the last decades. For this purpose, 8 major and 50 trace elements (TEs) were analyzed in 80 water and 65 rock samples, and interpreted with a quad-fold approach utilizing: (1) nonparametric statistical tests, (2) concentration deviations from ideal conservative freshwater–seawater mixing lines, (3) a new parameter called Mixing Enrichment Factor to assess the mobility of chemical constituents under salinizing conditions, and (4) 1-D dual porosity flow path modeling with PHREEQC. Dissolution/precipitation of Ca_xMg_ySr_zCO₃ and cation exchange were the main disclosed hydrogeochemical processes besides minor signs of organic matter oxidation. In the dolomitic limestone aquifer, less carbonate dissolved as compared to the limestone aquifer, partly because of lower pCO₂ in addition to seawater inflow triggering Mg-calcite precipitation by cation exchange. Saltwater intrusion led to mobilization of As, Ba, Cu, Ni, Rb, Sr and U in both aquifers, sometimes likely by cation exchange (e.g. Ba and Sr). Some of these TEs (notably Cu and Ni) recorded higher concentrations in the dolomitic limestone regardless of salinization. Other elements like Al, Be, Ce, Cr, Nb, Pb, V, Y and Zr revealed no or a low mobilization tendency. The concentration of all TEs in groundwater remained below drinking water limits notwithstanding moderate salinization. This classifies carbonate rocks as a weak geogenic source of TEs, whereas encroaching seawater appears to be a more important source.

Coastal aquifers and the deeper subsurface are increasingly exploited. The accompanying perforation of the subsurface for those purposes has increased the risk of short-circuiting of originally separated aquifers. This study shows how this short-circuiting negatively impacts the freshwater recovery efficiency (RE) during aquifer storage and recovery (ASR) in coastal aquifers. ASR was applied in a shallow saltwater aquifer overlying a deeper, confined saltwater aquifer, which was targeted for seasonal aquifer thermal energy storage (ATES). Although both aquifers were considered properly separated (i.e., a continuous clay layer prevented rapid groundwater flow between both aquifers), intrusion of deeper saltwater into the shallower aquifer quickly terminated the freshwater recovery. The presumable pathway was a nearby ATES borehole. This finding was supported by field measurements, hydrochemical analyses, and variable-density solute transport modeling (SEAWAT version 4; Langevin et al., 2007). The potentially rapid short-circuiting during storage and recovery can reduce the RE of ASR to null. When limited mixing with ambient groundwater is allowed, a linear RE decrease by short-circuiting with increasing distance from the ASR well within the radius of the injected ASR bubble was observed. Interception of deep short-circuiting water can mitigate the observed RE decrease, although complete compensation of the RE decrease will generally be unattainable. Brackish water upconing from the underlying aquitard towards the shallow recovery wells of the ASR system with multiple partially penetrating wells (MPPW-ASR) was observed. This leakage may lead to a lower recovery efficiency than based on current ASR performance estimations.

Well field Heel, in the south east of the Netherlands, consists of a row of wells drilled in an anoxic pyrite-containing aquifer alongside a former gravel pit, which now serves as a recharge basin, where water is actively aerated. All wells are seriously affected by chemical (screen slot) and/or mechanical (well bore) clogging. The objective of this study is to explain this combined occurrence. A combination of chemical, hydraulic and well-maintenance data indicate three groundwater quality types: (1) oxic basin water, (2) anoxic iron-containing basin water after oxidation of the traversed aquifer, and (3) deeply anoxic native groundwater. Wells abstracting a mixture of oxic basin water and anoxic basin water and/or native groundwater experience chemical well clogging, whereas wells abstracting (only or partly) native groundwater are vulnerable to mechanical well clogging. In the end, after oxic basin water has completely oxidized the traversed the aquifer, only two groundwater quality types will be present. Wells abstracting only oxic basin water will show no clogging, and wells abstracting a mixture of native groundwater and oxic basin water will experience chemical and possibly also mechanical well clogging. In this reasoning, the sequence in abstracted groundwater quality types coincides with a sequence in well clogging: from mechanical to chemical to no clogging. As well field Heel is situated in sloping terrain, the interplay between regional hydraulic gradient and different water qualities results in one-sided chemical clogging in the upper part of the well screen during abstraction, and in the lower part during the resting phase.

Managed Aquifer Recharge (MAR) is a promising method of increasing water availability in water stressed areas by subsurface infiltration and storage, to overcome periods of drought, and to stabilize or even reverse salinization of coastal aquifers. Moreover, MAR could be a key technique in making alternative water resources available, such as reuse of communal effluents for agriculture, industry and even indirect potable reuse. As exemplified by the papers in this Special Issue, consideration of water quality plays a major role in developing the full potential for MAR application, ranging from the improvement of water quality to operational issues (e.g., well clogging) or sustainability concerns (e.g., infiltration of treated waste water). With the application of MAR expanding into a wider range of conditions, from deserts to urban and coastal areas, and purposes, from large scale strategic storage of desalinated water and the reuse of waste water, the importance of these considerations are on the rise. Addressing these appropriately will contribute to a greater understanding, operational reliability and acceptance of MAR applications, and lead to a range of engineered MAR systems that help increase their effectiveness to help secure the availability of water at the desired quality for the future.

Artificial recharge of aquifers can be performed for various purposes and under varying hydrogeological conditions. We present an overview of deep-well recharge applications which have taken place in the Netherlands over the last two decades. We present the purpose of each application, the issues which had to be resolved, the preventive measures which were taken to improve performance and the lessons learned from each experience. Examples are given of applications which aimed at the storage of water for drinking and other purposes such as irrigation, achieving environmental goals and disposal of wastewater. Applications aiming at drinking water production usually faced issues related to the quality of the abstracted water not meeting drinking water standards with respect to various elements, such as iron, manganese and arsenic. Storage of water in brackish aquifers was complicated by buoyancy effects making part of the recharged water irrecoverable. Recharge of water with the purpose of recovering declined groundwater tables and fighting seawater intrusion was hindered by clogging of the injection well while the disposal of wastewater was limited to aquifers of lower groundwater quality.