**Importance of the Induced Recharge Term in Riverbank Filtration: Hydrodynamics, Hydrochemical, and Numerical Modelling Investigations**

#### **Rudy Rossetto 1,\* , Alessio Barbagli <sup>1</sup> , Giovanna De Filippis <sup>1</sup> , Chiara Marchina <sup>1</sup> , Thomas Vienken <sup>2</sup> and Giorgio Mazzanti <sup>3</sup>**


Received: 31 October 2020; Accepted: 5 December 2020; Published: 8 December 2020

**Abstract:** While ensuring adequate drinking water supply is increasingly being a worldwide challenging need, managed aquifer recharge (MAR) schemes may provide reliable solutions in order to guarantee safe and continuous supply of water. This is particularly true in riverbank filtration (RBF) schemes. Several studies aimed at addressing the treatment capabilities of such schemes, but induced aquifer recharge hydrodynamics from surface water bodies caused by pumping wells is seldom analysed and quantified. In this study, after presenting a detailed description of the Serchio River RBF site, we used a multidisciplinary approach entailing hydrodynamics, hydrochemical, and numerical modelling methods in order to evaluate the change in recharge from the Serchio river to the aquifer due to the building of the RBF infrastructures along the Serchio river (Lucca, Italy). In this way, we estimated the increase in aquifer recharge and the ratio of bank filtrate to ambient groundwater abstracted at such RBF scheme. Results highlight that in present conditions the main source of the RBF pumping wells is the Serchio River water and that the groundwater at the Sant'Alessio plain is mainly characterized by mixing between precipitation occurring in the higher part of the plain and the River water. Based on chemical mixing, a precautionary amount of abstracted Serchio River water is estimated to be on average 13.6 Mm<sup>3</sup> /year, which is 85% of the total amount of water abstracted in a year (~16 Mm<sup>3</sup> ). RBF is a worldwide recognized MAR technique for supplying good quality and reliable amount of water. As in several cases and countries the induced recharge component is not duly acknowledged, the authors suggest including the term "induced" in the definition of this type of MAR technique (to become then IRBF). Thus, clear reference may be made to the fact that the bank filtration is not completely due to natural recharge, as in many cases of surface water/groundwater interactions, but it may be partly/almost all human-made.

**Keywords:** drinking water supply; water supply scheme; surface water/groundwater interactions; managed aquifer recharge; induced riverbank filtration; groundwater resource management

#### **1. Introduction**

Ensuring adequate drinking water supply is one of the most pressing needs for our societies [1]. However, finding reliable sources is more and more a difficult task because of the widespread deterioration of surface water resources (and the related costs for treatment) and groundwater

deterioration and overexploitation [2–4]. In this sense managed aquifer recharge (MAR) schemes may provide reliable solutions in order to guarantee safe and continuous supply of water [5–7]. Managed aquifer recharge consists in the intentional recharge of aquifers using excess water in wet periods, or non-conventional water sources while at the same time assuring adequate protection of the environment and human-health [8]. Several techniques (spreading methods, Aquifer Storage and Recovery, Soil Aquifer Treatment systems, etc.) have been devised and widely adopted at global scale in the last 60 years [5,9]. Among these, riverbank filtration (RBF) has been in use since about two centuries, being a popular way to tap surface water ensuring at the same time adequate water treatment [10].

RBF is a technique in which the bed and bank of a river serve as treatment zone for the induced river water [10–14], when pumping wells are placed adjacently to a surface water body with continuous and adequate discharge or storage (such as rivers, lakes or basins). As such pumped groundwater in RBF schemes is constituted by a mix of ambient groundwater and induced surface water recharge. Several Authors (i.e., [11,15,16]) well-stated that while groundwater quality depends on the land use in the catchment area, bank filtrate quality is a function of river water quality and the efficiency of the purification processes during RBF. Moreover, other studies discuss the treatment capabilities of the RBF MAR systems in dealing with common surface water contaminants, but also with organics and emerging contaminants [17–20]. At the same time, while many examples exist in literature in evaluating surface- and groundwater interactions, at different spatial scales [21] with a number of methods [22], still induced recharge hydrodynamics is seldom analysed and quantified. Sottani and Vielmo (2014) [23] evaluated the extent and the rate of the recharge effects in groundwater due to weir realization in the middle Brenta River plain (Italy). Wei-shi et al. (2020) [24] evaluated the impact of a river reach restoration on the groundwater flow on bank filtration by using a transient flow and heat transport numerical model. Shankar et al. (2009) [15] quantified the contributing ratio of bank filtrate to ambient groundwater by means of a transient numerical model simulation at Grind well field (Germany).

The assessment of the induced recharge and of the ratio of bank filtrate to ambient groundwater is then important as it may provide relevant insights on the risk of groundwater contamination at the RBF wells. As such, in this study, we applied a multidisciplinary approach entailing hydrodynamics, hydrochemical and numerical modelling methods in order to evaluate the change in recharge from the Serchio river (Lucca, Italy) to the aquifer due to the building of the RBF infrastructures along the river. In this way, we estimated the increase in aquifer recharge and the ratio of bank filtrate to ambient groundwater abstracted at such RBF scheme.

#### **2. Materials and Methods**

The Serchio Riverbank Filtration scheme in Sant'Alessio (Lucca, Italy) supplies an average volume of 16 Mm<sup>3</sup> /year of drinking water (data from years 2016–2019; source GEAL SpA, local water utility) to about 300,000 citizens of Lucca, Pisa, and Livorno. This plant can be considered an exemplary site for MAR RBF systems. As a blue infrastructure, it provides ecosystem services such as "water storage" and "water quality improvement" [25].

The RBF scheme is set along the Serchio River (Figure 1). Groundwater is tapped by Holocene coarse sand and gravel aquifer overlain by a silty surficial cover. Water is pumped by twelve vertical wells (four at the Sant'Alessio pumping station and eight of the Pisa-Livorno water supply scheme), located between 30 m and 100 m from the river reach, inducing river water into the high yield (10−<sup>2</sup> m<sup>2</sup> /s transmissivity on average [26]) sand and gravel aquifer. The aquifer storage is then increased by the presence of one weir downstream the well field (Figure 2) to raise the river head of about 1.5 m, therefore raising the saturated part of the aquifer along the river reach.

**Figure 1.** Layout of the Serchio River riverbank filtration scheme in Lucca, Italy. Coordinates are in WGS84-UTM zone 32N (EPSGr: 32632). m gravel and medium sand in fine matrix, 9.5–12 m gravel in a silty-sandy matrix, 12–15.50 m gravel and sand, 15.50 - 17.3 m silty compacted clay.

**Figure 2.** The Sant'Alessio weir on the Serchio River. **Figure 2.** The Sant'Alessio weir on the Serchio River.

3 *2.2. Hydrodynamic Investigations*  Investigations of the aquifer properties were performed during different campaigns. The first campaign was performed to characterize the aquifer structure in vertical and horizontal extent around the most western RBF well (Figure 1), with a special focus placed on the characterization of subsurface heterogeneity. The successive field campaigns aimed at characterizing the flow field in the vicinity of this pumping well. Direct Push technology, which describes an approach of pushing, hammering, and/or vibrating small-diameter hollow steel rods into the ground [30], was employed Since the 1960s, the Sant'Alessio area along the Serchio River was deemed suitable for groundwater abstraction for drinking water purposes and, initially, four vertical wells were set in operation in 1967. They were drilled about 100 m away from the Serchio River to supply the north western households of the Lucca town with about 0.1 m<sup>3</sup> /s [27]. Following hydrogeological investigations, in the 1980s a highly yielding sand and gravel aquifer capable of supplying the towns of Pisa and Livorno, about 20 km and 40 km away, respectively, was identified. At the end of the 1980s, an Expert Commission evaluated the potential impacts of increasing abstractions and of the building of a weir to increase aquifer storage.

allows vertical differentiation of units with different hydraulic properties, measured as relative hydraulic conductivity (KDPIL), a parameter that can be linked to absolute hydraulic conductivity [35]. DPIL measurements were performed in vertical steps of 0.4 m; see [37] for detailed information on the tool and data analysis routines. In contrast, Direct Push electrical conductivity was employed for vertical high resolution identification of clay containing, confining layers or lenses as an increase in subsurface electrical conductivity can be related to an increase of clay mineral content under non-saline conditions; see [33,38] for further information. Thirteen DPIL profiles as well as 10 Direct Push electrical conductivity logs were collected during the initial site investigation campaign.

Several tracer tests were additionally performed to characterize the flow field in the vicinity of the pumping well. Tracer test included heat tracer as well as salt tracer testing; the latter was combined with Direct Push electrical conductivity logging for enhanced in-situ monitoring,

From the hydrodynamic point of view, three main flow directions are converging to the RBF plant (Figure 5): (i) a flow from the alluvial plain northern boundary (San Quirico-Carignano hills) given by the local meteoric recharge infiltrating through low-permeability deposits (from north to south flow) and then infiltrating in the aquifer; (ii) from east to west and from south to north flows. During the year, no seasonal changes are affecting these main flows directions given by the Serchio

to end of the Direct Push rod string [31–36].

following the approach of Vienken et al. (2017) [39].

during the site investigation to install groundwater monitoring wells and to collect vertical high resolution subsurface information. The latter are obtained using specialized probes that are attached During the following ten years, the full MAR scheme was completed with the building of the river weir in Sant'Alessio, and it reached a total of about 0.500 m<sup>3</sup> /s supply capacity. More information on the historical development and the socio-economic importance of the Serchio river MAR scheme may be found in Rossetto et al. 2020a [28].

### *2.1. Site Description*

The Sant'Alessio plain surface geology is characterized by unconsolidated silty to sandy sediments covering a sandy-gravel aquifer (Figure 3). The aquifer is limited at the bottom by a level of silty clays of lacustrine origin, which outcrops as stiff clays in the northern part of the area. From the Serchio river (south) towards the San Quirico-Carignano hills (north) the silty-sandy cover of the aquifer becomes thinner, going from 6.5–7 m to less than 1 m (Figure 3 Bottom). According to Rossetto and Bockelmann-Evans (2007) [29], also the aquifer thickness follows the same trend, starting from about 30 m in the Serchio riverbed and pinching-out towards the San Quirico-Carignano hills, where the silty-sandy cover directly meets the clay sediments at the aquifer bottom (Figure 3 Top). The Serchio river shows variable discharge ranging from more than 1000 m<sup>3</sup> /s during extreme rainfall events down to few m<sup>3</sup> /s during the summer season in extreme dry years. In the Serchio riverbed the aquifer is covered by coarse and clean gravels outcrops, as well as sand bars.

**Figure 3.** Simplified stratigraphic cross section (top) and detailed stratigraphic logs around the Serchio River RBF area (bottom).

2

Outcrops along the Serchio river show a strong contrast in the sedimentological composition of the river deposits which may lead to changes in hydraulic properties over several orders of magnitude on short vertical distances. The vertical heterogeneity of the aquifer was further investigated by a 17 m deep core drilled directly in the riverbed during a low-flow period (December 2015; Figure 4). According to the drilled core, the river bed stratigraphy was then reconstructed as follows (from the top of the river bed): 0–2 m gravels and sand, 2–5 m sand with gravel in a weakly silty matrix, 5–9.5 m gravel and medium sand in fine matrix, 9.5–12 m gravel in a silty-sandy matrix, 12–15.50 m gravel and sand, 15.50–17.3 m silty compacted clay. *Hydrology* **2020**, *7*, x FOR PEER REVIEW 6 of 21


**Figure 4.** Stratigraphic log from the top of the Serchio river riverbed down to about 17 m. **Figure 4.** Stratigraphic log from the top of the Serchio river riverbed down to about 17 m.

#### *2.2. Hydrodynamic Investigations*

Investigations of the aquifer properties were performed during different campaigns. The first campaign was performed to characterize the aquifer structure in vertical and horizontal extent around the most western RBF well (Figure 1), with a special focus placed on the characterization of subsurface heterogeneity. The successive field campaigns aimed at characterizing the flow field in the vicinity of this pumping well. Direct Push technology, which describes an approach of pushing, hammering, and/or vibrating small-diameter hollow steel rods into the ground [30], was employed during the site investigation to install groundwater monitoring wells and to collect vertical high resolution subsurface information. The latter are obtained using specialized probes that are attached to end of the Direct Push rod string [31–36].

In this study, the Direct Push Injection Logger (DPIL) was used together with Direct Push electrical conductivity logging. The DPIL is used for vertical in-situ hydrostratigraphic profiling as it allows vertical differentiation of units with different hydraulic properties, measured as relative hydraulic conductivity (KDPIL), a parameter that can be linked to absolute hydraulic conductivity [35]. DPIL measurements were performed in vertical steps of 0.4 m; see [37] for detailed information on the tool and data analysis routines. In contrast, Direct Push electrical conductivity was employed for vertical high resolution identification of clay containing, confining layers or lenses as an increase in subsurface electrical conductivity can be related to an increase of clay mineral content under non-saline conditions; see [33,38] for further information. Thirteen DPIL profiles as well as 10 Direct Push electrical conductivity logs were collected during the initial site investigation campaign.

Several tracer tests were additionally performed to characterize the flow field in the vicinity of the pumping well. Tracer test included heat tracer as well as salt tracer testing; the latter was combined with Direct Push electrical conductivity logging for enhanced in-situ monitoring, following the approach of Vienken et al. (2017) [39].

From the hydrodynamic point of view, three main flow directions are converging to the RBF plant (Figure 5): (i) a flow from the alluvial plain northern boundary (San Quirico-Carignano hills) given by the local meteoric recharge infiltrating through low-permeability deposits (from north to south flow) and then infiltrating in the aquifer; (ii) from east to west and from south to north flows. During the year, no seasonal changes are affecting these main flows directions given by the Serchio River filtration into the aquifer. The presence of vertical upward groundwater flow was tested by means of multi-level piezometers, but no relevant vertical gradient was highlighted. The piezometric level at the well field shows a direct relationship with the Serchio river stage, with the range between the maximum and minimum groundwater head no larger than 2 m between 2011 and 2016 (Figure 6).

River discharge measurements were performed in May 2015 in the Serchio river at two cross-sections corresponding to the 2 weirs (Figure 1), measuring the flow speed with an acoustic flow meter (acoustic digital current meter, OTT Messtechnik GmbH) every 50 cm along three verticals.

#### *2.3. Hydrochemical Characterisation*

Surface water and groundwater sampling campaigns for hydrochemical characterization were repeated both in dry (2 campaigns from May 2015 to October 2015 and 3 from May 2016 to October 2016) and in wet (2 campaigns from November 2015 to April 2016) periods. Data were collected from 25 different sampling points (including multilevel piezometer clusters at point PR\_01-02-03-04-05-06) to assess hydrochemical seasonal variability. Groundwater samples were mostly collected using a low-flow pump after purging at least three well volumes and reaching stable values of the chemical-physical parameters monitored during sampling (Temperature, Electric Conductivity, Redox Potential, and Dissolved Oxygen). Surface water samples were collected 30 cm below the water surface, where the flow was sufficiently high to guarantee a representative sample. Sampling campaigns for isotopic characterization were performed on a selected sub sample of the monitoring points (14 in total) and on 2 different periods: dry season (May 2015) and wet season (January 2016).

*Hydrology* **2020**, *7*, x FOR PEER REVIEW 7 of 21

**Figure 5.** Monitoring network (groundwater head and hydrochemistry) in relation to the main groundwater flow directions. Piezometric head data from September 2003 (dry season; data from Studio Nolledi 2003). **Figure 5.** Monitoring network (groundwater head and hydrochemistry) in relation to the main groundwater flow directions. Piezometric head data from September 2003 (dry season; data from Studio Nolledi 2003). *Hydrology* **2020**, *7*, x FOR PEER REVIEW 8 of 21

concentrated Suprapur HNO3) was taken for the determination of cations and trace elements. After sampling all the aliquots were stored at 4°C until analyzed. Cations and trace elements were measured by inductively coupled plasma mass spectrometry (ICP-MS) using a Thermo-Scientific X Series instrument (Thermo, Bremen, Germany). The samples were previously diluted 1:10 by deionized Milli-Q water (resistivity of ca. 18.2 MΩ×cm) and a **Figure 6.** Monthly groundwater head data at the Serchio River RBF plant and Serchio river stage at the Ponte di Monte San Quirico monitoring station from May 2011 to December 2016 (data courtesy of Servizio Idrologico Regionale—Regione Toscana). The locations of the piezometer and the Serchio river monitoring station are marked in Figure 4. **Figure 6.** Monthly groundwater head data at the Serchio River RBF plant and Serchio river stage at the Ponte di Monte San Quirico monitoring station from May 2011 to December 2016 (data courtesy of Servizio Idrologico Regionale—Regione Toscana). The locations of the piezometer and the Serchio river monitoring station are marked in Figure 4.

known amount of Re and Rh have been introduced as internal standard; in each analytical session. The analysis of samples was verified with that of the reference materials EU-L-1 and ES-L1 provided by SCP-Science (www.scpscience.com). Anions were determined by ion chromatography using a Dionex ICS-1000 (Thermo, Bremen, Germany), calibrated by different dilutions of the Dionex "7-ion standard". The coherence of chemical data was verified by checking the ionic balance, as the sum of cation (expressed in meq/L) approaches that of anions with a relative error less than 5%. Hydrogen and oxygen isotope ratios were determined using the CRDS Los Gatos LWIA 24-d isotopic analyzer (Los Gatos Research, California, USA). The isotopic ratios of 2H/1H and 18O/16O are expressed as δ notation [δ = Water samples were filtered by 0.45 µm filters (Minisart® NML syringe cellulose acetate filters) and separated in different aliquots. 50 mL PE bottles were taken for the analysis of anions and oxygen/hydrogen isotopes in each site. An additional 50 mL PE bottle (acidified with 0.5 mL of concentrated Suprapur HNO3) was taken for the determination of cations and trace elements. After sampling all the aliquots were stored at 4 ◦C until analyzed.

software has been applied at several case studies; details may be found in, i.e., [48,49].

A groundwater flow numerical model was built by means of the FREEWAT software [41,42] using the MODFLOW-2005 code [43]. FREEWAT, as its predecessor SID&GRID [44], is a free and open source platform, QGIS-integrated, for planning and management of ground- and interactions with surface–water resources. It is the main result of the HORIZON 2020 FREEWAT project (FREE and open source software tools for WATer resource management [45]). FREEWAT provides tools for archiving, analyzing and processing the groundwater data and information (e.g., sensor and field data on groundwater levels and hydrochemical data) along with post-processing tools for visualization and analysis of the model results [46,47]. MODFLOW is one of the codes most widely used in the world for numerical simulation of groundwater flow in aquifers. It solves the groundwater flow equation in three dimensions using a finite difference scheme. The FREEWAT

We simulated an active domain extending over 4.2 km2 discretized using square cells 100 m2 wide and two model layers, a silty-sandy superficial cover and the sandy-gravelly aquifer. The conceptual model for the simulated domain is shown in Figure 7. The northern boundary of the active domain, corresponding to the outcrop of impervious clayey sediments of the Monte S. Quirico-Carignano hills, was defined according to the geological setting of the area (Figure 3). Hydraulic connection with the Serchio river and exploitation of the aquifer through the Sant'Alessio well field were simulated. The model has been implemented to simulate two years, from 20 November 2014 until 30 November 2016. The simulation started with a steady-state stress period

*2.4. Groundwater Flow Numerical Modelling* 

((Rsample/Rstandard) − 1) × 1000] with respect to the Vienna Standard Mean Ocean Water

Cations and trace elements were measured by inductively coupled plasma mass spectrometry (ICP-MS) using a Thermo-Scientific X Series instrument (Thermo, Bremen, Germany). The samples were previously diluted 1:10 by deionized Milli-Q water (resistivity of ca. 18.2 MΩ×cm) and a known amount of Re and Rh have been introduced as internal standard; in each analytical session. The analysis of samples was verified with that of the reference materials EU-L-1 and ES-L1 provided by SCP-Science (www.scpscience.com).

Anions were determined by ion chromatography using a Dionex ICS-1000 (Thermo, Bremen, Germany), calibrated by different dilutions of the Dionex "7-ion standard". The coherence of chemical data was verified by checking the ionic balance, as the sum of cation (expressed in meq/L) approaches that of anions with a relative error less than 5%. Hydrogen and oxygen isotope ratios were determined using the CRDS Los Gatos LWIA 24-d isotopic analyzer (Los Gatos Research, California, USA). The isotopic ratios of <sup>2</sup>H/ <sup>1</sup>H and <sup>18</sup>O/ <sup>16</sup>O are expressed as <sup>δ</sup> notation [<sup>δ</sup> <sup>=</sup> ((Rsample/Rstandard) <sup>−</sup> 1) <sup>×</sup> 1000] with respect to the Vienna Standard Mean Ocean Water (VSMOW) international standard. Four bracketing standards were run throughout the analytical sessions, as indicated in Marchina et al. (2020) [40].

#### *2.4. Groundwater Flow Numerical Modelling*

A groundwater flow numerical model was built by means of the FREEWAT software [41,42] using the MODFLOW-2005 code [43]. FREEWAT, as its predecessor SID&GRID [44], is a free and open source platform, QGIS-integrated, for planning and management of ground- and interactions with surface–water resources. It is the main result of the HORIZON 2020 FREEWAT project (FREE and open source software tools for WATer resource management [45]). FREEWAT provides tools for archiving, analyzing and processing the groundwater data and information (e.g., sensor and field data on groundwater levels and hydrochemical data) along with post-processing tools for visualization and analysis of the model results [46,47]. MODFLOW is one of the codes most widely used in the world for numerical simulation of groundwater flow in aquifers. It solves the groundwater flow equation in three dimensions using a finite difference scheme. The FREEWAT software has been applied at several case studies; details may be found in, i.e., [48,49].

We simulated an active domain extending over 4.2 km<sup>2</sup> discretized using square cells 100 m<sup>2</sup> wide and two model layers, a silty-sandy superficial cover and the sandy-gravelly aquifer. The conceptual model for the simulated domain is shown in Figure 7. The northern boundary of the active domain, corresponding to the outcrop of impervious clayey sediments of the Monte S. Quirico-Carignano hills, was defined according to the geological setting of the area (Figure 3). Hydraulic connection with the Serchio river and exploitation of the aquifer through the Sant'Alessio well field were simulated. The model has been implemented to simulate two years, from 20 November 2014 until 30 November 2016. The simulation started with a steady-state stress period lasting 11 days, and then with 24 transient monthly stress periods. Calibration was run using a mixed trial and error and automatic parameter estimation approach [50]. Details and more information on the model implementation may be found in Rossetto et al. (2017) [51].

The calibrated model was then used to simulate the evolution of the RBF scheme in order to evaluate the change in recharge from the Serchio river to the aquifer. The following scenarios were simulated:


simulated:

lasting 11 days, and then with 24 transient monthly stress periods. Calibration was run using a mixed trial and error and automatic parameter estimation approach [50]. Details and more

(a) Simulation of the aquifer river exchange prior to any riverbed modification or aquifer exploitation in order to evaluate pristine conditions: to this aim, we linearly reconstructed the

(b) Simulation of the aquifer river exchange prior to any riverbed modification considering aquifer

(c) Simulation of the aquifer river exchange considering the Sant'Alessio weir in operation and

(d) Simulation of the aquifer river exchange considering the Sant'Alessio weir in operation and

The calibrated model was then used to simulate the evolution of the RBF scheme in order to evaluate the change in recharge from the Serchio river to the aquifer. The following scenarios were

information on the model implementation may be found in Rossetto et al. (2017) [51].

riverbed profile prior to the weir construction;

exploitation for 0.350 m3/s;

aquifer exploitation for 0.350 m3/s;

**Figure 7.** Conceptual model for the simulated domain; main hydrological elements are shown. **Figure 7.** Conceptual model for the simulated domain; main hydrological elements are shown.

#### **3. Result and discussion 3. Result and Discussion**

In order to answer the research questions, we jointly used the results of the performed analysis to provide data for the definition of the conceptual model and for the implementation of the numerical modelling (i.e., data coming from the hydrodynamic investigations) or to cross-validate the results of each approach. In order to answer the research questions, we jointly used the results of the performed analysis to provide data for the definition of the conceptual model and for the implementation of the numerical modelling (i.e., data coming from the hydrodynamic investigations) or to cross-validate the results of each approach.

Results from the hydrodynamics investigations informed the implementation of the numerical model and provided additional insights on the groundwater flow field. DPIL measurements revealed a strong variation of DPIL values on short vertical distances within the aquifer, while measured electrical conductivity values showed only minor variation, ranging between 10 mS/m and 20 mS/m. These are typical values for clean sands and gravels. Figure 8 depicts Direct Push Results from the hydrodynamics investigations informed the implementation of the numerical model and provided additional insights on the groundwater flow field. DPIL measurements revealed a strong variation of DPIL values on short vertical distances within the aquifer, while measured electrical conductivity values showed only minor variation, ranging between 10 mS/m and 20 mS/m. These are typical values for clean sands and gravels. Figure 8 depicts Direct Push logging results at an investigation point located approximately 90 m north-east of the pumping well C5 (Figure 5). Based on the results, the identification of hydrostratigraphic units or prevailing layers within the aquifer itself was not possible. Hence, Direct Push based pneumatic slug testing (see [52]) was performed at four locations in proximity of the pumping well C5 at different depth intervals between 3.6 and 10 m below ground surface to obtain further information on the distribution of hydraulic conductivity. Measured values ranged between 5.6 <sup>×</sup> <sup>10</sup>−<sup>3</sup> <sup>m</sup>/s and 3.3 <sup>×</sup> <sup>10</sup>−<sup>4</sup> <sup>m</sup>/s.

Thanks to the tracer tests a maximum bulk groundwater flow velocity of 26 m/d was estimated based on the first arrival times of the sodium-chloride tracer test breakthrough curve measurement. Despite the use of vertical high-resolution exploration, it was not possible to resolve the aquifer heterogeneity and to obtain in-depth understanding of the hydrodynamics using field investigations alone. This clearly highlights the necessity to approach experimental complex questions in groundwater hydrology by combining sophisticated site investigation and advanced modelling.

logging results at an investigation point located approximately 90 m north-east of the pumping well C5 (Figure 5). Based on the results, the identification of hydrostratigraphic units or prevailing layers within the aquifer itself was not possible. Hence, Direct Push based pneumatic slug testing (see [52]) was performed at four locations in proximity of the pumping well C5 at different depth intervals between 3.6 and 10 m below ground surface to obtain further information on the distribution of

Thanks to the tracer tests a maximum bulk groundwater flow velocity of 26 m/d was estimated based on the first arrival times of the sodium-chloride tracer test breakthrough curve measurement. Despite the use of vertical high-resolution exploration, it was not possible to resolve the aquifer heterogeneity and to obtain in-depth understanding of the hydrodynamics using field investigations

hydraulic conductivity. Measured values ranged between 5.6 × 10−3 m/s and 3.3 × 10−4 m/s.

groundwater hydrology by combining sophisticated site investigation and advanced modelling.

**Figure 8.** Direct Push profiling results (DPIL and Electrical Conductivity measurements). **Figure 8.** Direct Push profiling results (DPIL and Electrical Conductivity measurements).

The simulated water budget for one hydrological year (October 2015 to September 2016) is presented in Table 1. As far as the inflow terms, the Serchio river recharge constitutes the bulk of aquifer recharge in the Sant'Alessio plain with an overall yearly volume of about 58 Mm3, in line with the loss of river water estimated by means of flow measurements (1.4 m3/s). This recharge occurs within the reach Ponte di Monte San Quirico—Sant'Alessio weir (Figures 1 and 9). A smaller inflow term is constituted by meteoric recharge and inflow from the adjoining hills. Outflows from the systems are dominated by RBF pumping wells abstraction (about 16 Mm3), outflow from the southern boundary (about 8 Mm3), and drainage of the aquifer operated by the river south of the Sant'Alessio weir. The effect of the Freddanella ditch in draining the aquifer is negligible. The simulated head contour map (Figure 9) clearly shows this situation, presenting also a large drainage axis at about the middle of the plain. The simulated water budget for one hydrological year (October 2015 to September 2016) is presented in Table 1. As far as the inflow terms, the Serchio river recharge constitutes the bulk of aquifer recharge in the Sant'Alessio plain with an overall yearly volume of about 58 Mm<sup>3</sup> , in line with the loss of river water estimated by means of flow measurements (1.4 m<sup>3</sup> /s). This recharge occurs within the reach Ponte di Monte San Quirico—Sant'Alessio weir (Figures 1 and 9). A smaller inflow term is constituted by meteoric recharge and inflow from the adjoining hills. Outflows from the systems are dominated by RBF pumping wells abstraction (about 16 Mm<sup>3</sup> ), outflow from the southern boundary (about 8 Mm<sup>3</sup> ), and drainage of the aquifer operated by the river south of the Sant'Alessio weir. The effect of the Freddanella ditch in draining the aquifer is negligible. The simulated head contour map (Figure 9) clearly shows this situation, presenting also a large drainage axis at about the middle of the plain.

**Table 1.** Simulated water budget for the Sant'Alessio aquifer for the hydrologic year October 2015 to September 2016.


**Figure 9.** Simulated hydraulic head at the end of May 2016.

Table 2 shows the results of the analyses run in order to evaluate the change in recharge from the Serchio river to the aquifer induced during the evolution of the RBF scheme. In pristine conditions (riverbed not modified and wells not built), simulation results show that the natural recharge of the aquifer was of about 0.150 m<sup>3</sup> /s. Construction of the first set of wells and their operation at an overall abstraction rate of 0.350 m<sup>3</sup> /s would have caused an estimated recharge of about 0.490 m<sup>3</sup> /s, hence inducing in the aquifer around 0.340 m<sup>3</sup> /s. We then estimated in about 0.030 m<sup>3</sup> /s the additional recharge provided as result of the construction in the 1990s of the Sant'Alessio weir. This recharge then creates an additional aquifer storage of about 1 Mm<sup>3</sup> . Finally, increasing groundwater withdrawal up to 0.430 m<sup>3</sup> /s (from initial 0.350 m<sup>3</sup> /s) induces an additional recharge of about 0.09 m<sup>3</sup> /s.


**Table 2.** Simulated changes in river recharge from pristine conditions to completion of the Serchio River RBF managed aquifer recharge scheme.

4 These results show that the increase in recharge induced by pumping is larger than that caused by the weir construction. The latter brought a limited, yet important, change in storage. Figure 10 shows the groundwater head increase in the aquifer after the construction of the Sant'Alessio weir compared to the previous condition, with only the well field in operation. At the RBF scheme, the river/aquifer connection ensures, thanks to Serchio river discharge, stable recharge of the aquifer induced by the pumping wells and by the presence of the river weir. Hence, this also guarantees a limited drawdown of the water table.

**Figure 10.** Simulated difference in groundwater head between pre- and post-operam of the Sant'Alessio weir.

3 As far as water hydrochemistry, according to the Piper classification diagrams (Figure 11 bottom), the samples mainly exhibit a Na-K-HCO<sup>3</sup> hydrochemical facies with a clear trend towards the Ca-Mg-SO<sup>4</sup> facies, perfectly aligned on a mixing line from the point 221 to the Serchio river (see Figures 5 and 7). The only exception to this trend is the sample at the "Freddanella" ditch, a canal draining the Sant'Alessio plain and also collecting untreated wastewater, hence to be considered as an end-member for potential groundwater pollution in the conceptual model [53]. This canal did not show any hydrochemical facies predominance during the dry season 2015.

The mixing process between surface water and ambient groundwater is also confirmed by the binary plots of conservative elements such as Cl− and Br− (Figure 11, top right). However, as for the geochemical facies, the mixing process of conservative elements is more pronounced during the dry season, when the chloride content of the Serchio river is not affected by dilution due to rainfall, and, at the same time, the Clconcentrations are more uniform in the groundwater.

The δ <sup>18</sup>O and δ <sup>2</sup>H values of water are very uniform in most of the samples and they have isotopic values ranging between −8.2‰ and −4.9‰ for δ <sup>18</sup>O, and <sup>−</sup>31.7‰ and <sup>−</sup>49‰ for <sup>δ</sup>D. These values are plotted in Figure 11 (top left), together with the Global Meteoric Water Line (GMWL) calculated by Rozansky (1993) [54] and the Local Meteoric Water Lines calculated by Longinelli and Selmo (2003) [55] and La Ruffa and Panichi (2000) [56]. The results evidence a trend of precipitation typical for the west coast of Italy. This data is also confirmed by an average deuterium excess of 14.6‰, respect the GMWL, which is in the typical range for the northern Tuscan area (10.8–14.6‰; Longinelli and Selmo, 2003) [55].

Isotopic data describes a clear mixing process between the Serchio river water and the groundwater. The isotopic composition of the local recharge (point 221, in Figures 7 and 11) is the least depleted (in light isotopes, i.e., <sup>16</sup>O and <sup>1</sup>H), while the Serchio river reveals an isotopic signature similar to those of the groundwater pumped at the RBF wells. The groundwater in the Sant'Alessio plain and at the RBF plant shows a clear signature characterized by the Serchio river. This River originates in the Apuan Alps catchment (with maximum elevation of about 2000 m amsl, and average elevation of 717 m amsl) and therefore shows values more depleted (therefore richer in light isotopes) than the ones

of the local recharge (originating in the San Quirico-Carignano hills, average elevation ~40 m amsl). In particular, the point sampled in the Sant'Alessio plain (MAR\_10) is placed almost at the middle of the mixing line (Figure 11 left). The mixing processes highlighted by the δ <sup>18</sup>O and δ <sup>2</sup>H can be detected during the dry and, also, the wet season. However, during the wet period, the Serchio river shows an enrichment in heavy isotopes due to the impact of local precipitation and the influence of tributaries from lower sub-catchments in the upstream part of the basin and, at the same time, its influence on the groundwater of the Sant'Alessio plain (also at MAR\_10) is higher due to its higher stage and discharge. *Hydrology* **2020**, *7*, x FOR PEER REVIEW 14 of 21

**Figure 11.** Piper diagram (**bottom**), Cl<sup>−</sup> vs Br binary plot (**top right**) and δ18O vs δ2H binary plot (**top left**) of the samples collected during the dry season of 2015. **Figure 11.** Piper diagram (**bottom**), Cl<sup>−</sup> vs Br- binary plot (**top right**) and δ <sup>18</sup>O vs δ <sup>2</sup>H binary plot (**top left**) of the samples collected during the dry season of 2015.

These results confirm that the main sources of the RBF pumping wells is the Serchio river water and that the groundwater at the Sant'Alessio plain is mainly characterized by mixing between precipitation occurring in the higher part of the plain and the river water. As an example, the percentage of Serchio river water at the RBF wells ranges from nearly 100% (during the wet season) to around 80% (in the dry period), while for the point MAR\_10 sample fraction of about 50% and 70% of Serchio river water may be estimated for the dry and wet season, respectively (Figure 11). Based on chemical mixing, a precautionary amount of abstracted Serchio river water is estimated to be on average 13.6 Mm3/year, this is 85% of the total amount of water abstracted in a year (~16 Mm3). As such, only the 15% (2.4 Mm3) of abstracted groundwater comes from ambient groundwater. Therefore, as suggested both by the hydrodynamic and the hydrochemical data, the Serchio river and the local recharge from the northern part of the aquifer (groundwater entering the aquifer along the northern hilly border) can be considered end-member terms of the mixing system. Based on this assumption, different seasonal mixing fractions were computed using the distribution of lithium as conservative tracer. Lithium is the trace element showing the widest range of values among the end-members in the Sant'Alessio area (with high concentration in the Serchio river and values below the detection limit in the point 221; Table 3). The mixing fractions showed the extension of the Serchio river recharge within the Sant'Alessio aquifer and its importance at the RBF plant (Figure 12).

*Hydrology* **2020**, *7*, 96


**Table 3.** Average values during dry and wet season of the presented analytes.

*Hydrology* **2020**, *7*, x FOR PEER REVIEW 15 of 21

**Figure 12.** Interpolation showing the Serchio river influence on the Sant'Alessio aquifer and the Serchio River RBF, using lithium as tracer. **Figure 12.** Interpolation showing the Serchio river influence on the Sant'Alessio aquifer and the Serchio River RBF, using lithium as tracer.

These results confirm that the main sources of the RBF pumping wells is the Serchio river water and that the groundwater at the Sant'Alessio plain is mainly characterized by mixing between precipitation occurring in the higher part of the plain and the river water. As an example, the percentage of Serchio river water at the RBF wells ranges from nearly 100% (during the wet season) to around 80% (in the dry period), while for the point MAR\_10 sample fraction of about 50% and 70% of Serchio river water may be estimated for the dry and wet season, respectively (Figure 11). Based on chemical mixing, a precautionary amount of abstracted Serchio river water is estimated to be on average 13.6 Mm<sup>3</sup> /year, this is 85% of the total amount of water abstracted in a year (~16 Mm<sup>3</sup> ). As such, only the 15% (2.4 Mm<sup>3</sup> ) of abstracted groundwater comes from ambient groundwater.

#### **4. Conclusions**

In this study, we presented a detailed description of the Serchio River RBF site, and we used a multidisciplinary approach entailing hydrodynamics, hydrochemical, and numerical modelling methods in order to evaluate the change in recharge from the Serchio river to the aquifer due to the building of the RBF infrastructures along the river, and the ratio of bank filtrate to ambient groundwater abstracted at such RBF scheme.

The use of Direct Push technology in combination with traditional site investigation approaches was feasible to provide critical information for reliable characterization of subsurface properties and to inform first the conceptual model and then the implementation of the numerical model. Thereby, the use of DPIL proved successful to produce quantitative information about vertical distribution of

hydraulic conductivities in a high K environment. Direct Push slug testing as well as tracer testing was successfully used for the in-situ characterization of hydraulic conductivity variations and to capture high flow velocities at the river/aquifer interface in the vicinity of the most western pumping well that were in the order of m/day.

Our results show that when discussing surface water/groundwater interactions, the Serchio river water constitutes the main source of water for the pumping wells at the RBF scheme. By means of numerical modelling of groundwater flow we simulated the increase in induced recharge, starting from natural conditions, up to present conditions with the Sant'Alessio weir and the vertical wells in full operation. The amount of induced recharge caused by the pumping wells is by far larger than that related to the construction of the Sant'Alessio weir raising the riverbed of about 1.5 m. The weir is however important as it contributes to guarantee a small drawdown in the well-field area. Hydrochemical investigations confirmed and specified the results of the numerical modelling activities. Groundwater at the Sant'Alessio plain is mainly characterized by mixing between precipitation occurring in the higher part of the plain and the river water. This mix seasonally varies depending on the Serchio river discharge, with a marked decrease of the induced component during the Serchio low-flow periods. As an example, the percentage of abstracted Serchio river water at the RBF wells ranges from nearly 100% (in the wet season) to around 80% (in the dry season). Based on chemical mixing, a precautionary amount of abstracted Serchio river water is estimated to be on average 13.6 Mm<sup>3</sup> /year, which is 85% of the total amount of water abstracted in a year (~16 Mm<sup>3</sup> ). Finally, we wish to stress that the assessment of the induced recharge and of the ratio of bank filtrate to ambient groundwater is very important at RBF sites as it may provide relevant insights on the risk of groundwater contamination. This is especially true for the Serchio River RBF as it was already affected by surface water related pollution at the beginning of the 2000 [29]. Disregarding the interactions between surface and groundwater brought to search the cause of contamination in the adjoining agricultural areas and to issue inappropriate local regulations for farming activities [29]. All these issues were discussed in Rodríguez-Escales et al. (2018) [57] in comparison with other MAR sites across the Mediterranean basin: the risk of failure of the Serchio River RBF scheme was defined low. Nowadays, the Serchio River RBF scheme is a reliable site for continuous water supply.

RBF is a worldwide recognized MAR technique for supplying good quality and reliable amounts of water. However, in many cases and countries (including Italy) government authorities do not appear to acknowledge that groundwater extraction is sustained by induced recharge through the river bank and it is not entirely from natural recharge/ambient groundwater. Additionally, several studies take this connection for granted and often only focus on bank filtration treatment capabilities disregarding the amount of recharge induced by the pumping wells. This "induced" component is even not taken into account in the definition of this MAR type in the literature (i.e., [58]). In fact, this remains one of the few MAR schemes where reference to the recharge mechanism is not adequately made in the definition, conversely to, i.e., spreading methods.

Because of this, we suggest to include the term "induced" to complete the definition of this type of MAR technique, hence defining it "*Induced RiverBank Filtration*". This way we make reference to the fact that the bank filtration is not completely due to natural recharge, as in several cases of natural surface water/groundwater interactions, but it may be partly/almost all human-made.

**Author Contributions:** Conceptualization, R.R. and T.V.; methodology, R.R., T.V., A.B., C.M.; software, R.R. and G.D.F.; validation, R.R. and T.V.; formal analysis, R.R., A.B., T.V., C.M., G.D.F.; investigation, R.R., A.B., T.V., C.M., G.D.F.; data curation, R.R., A.B., T.V., G.D.F.; writing—original draft preparation, R.R. and A.B.; writing—review and editing, R.R., A.B., T.V., C.M., G.D.F.; visualization, A.B., C.M. and G.D.F.; supervision, R.R.; project administration, R.R. and G.M.; funding acquisition, R.R. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the European Union, Grant Number 619120.

**Acknowledgments:** The work presented here was performed during the EU co-funded project FP7 MARSOL (Grant Agreement No. 619120). This research also exploited results from the H2020 FREEWAT project, funded by the European Union within the Horizon 2020 research and innovation programme (Grant number 642224). Alessio Barbagli (now Department of Physics and Earth Sciences, University of Ferrara, Italy) took part in the research activities during his three-year PhD in Agrobiosciences at the Institute of Life Sciences, Scuola Superiore Sant'Anna. Giovanna De Filippis (now AECOM URS, Milano, Italy) and Chiara Marchina (now Department of Land, Environment, Agriculture and Forestry, University of Padova, Italy) took part to the research activities during a post-doc scholarship at the Institute of Life Sciences, Scuola Superiore Sant'Anna. Thomas Vienken (now also Weihenstephan-Triesdorf University of Applied Sciences, TUM Campus Straubing for Biotechnology and Sustainability, Straubing, Germany) performed these activities while at UFZ. The authors thanks GEAL S.p.a. for the technical support provided and granting access to the well-field during the research activities.

**Conflicts of Interest:** The authors declare no conflict of interest.

### **References**


58. Natural Resource Management Ministerial Council, Environment Protection and Heritage Council National Health and Medical Research Council. Australian Guidelines for Water Recycling, Managing Health and Environmental Risks, Vol 2C: Managed Aquifer Recharge. Available online: http: //webarchive.nla.gov.au/gov/20130904195601/http://www.environment.gov.au/water/publications/quality/ water-recycling-guidelines-mar-24.html (accessed on 29 September 2020).

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

*Article*
