Next Article in Journal
An Assessment of the Recent Evolution of the Streamflow in a Near-Natural System: A Case Study in the Headwaters of the Mero Basin (Galicia, Spain)
Next Article in Special Issue
Water Budget Analysis Considering Surface Water–Groundwater Interactions in the Exploitation of Seasonally Varying Agricultural Groundwater
Previous Article in Journal
Investigating Tradeoffs between Agricultural Development and Environmental Flows under Climate Change in the Stung Chinit Watershed, Cambodia
Previous Article in Special Issue
Meeting SDG6 in the Kingdom of Tonga: The Mismatch between National and Local Sustainable Development Planning for Water Supply
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Institute of Life Sciences, Scuola Superiore Sant’Anna, 56127 Pisa, Italy
2
Helmholtz Centre for Environmental Research GmbH—UFZ, Permoserstraße 15, 04318 Leipzig, Germany
3
Genio Civile Toscana Nord, Regione Toscana, 55100 Lucca, Italy
*
Author to whom correspondence should be addressed.
Submission received: 31 October 2020 / Revised: 28 November 2020 / Accepted: 5 December 2020 / Published: 8 December 2020
(This article belongs to the Special Issue Integrated Surface Water and Groundwater Analysis)

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 Mm3/year, which is 85% of the total amount of water abstracted in a year (~16 Mm3). 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.

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,3,4]. In this sense managed aquifer recharge (MAR) schemes may provide reliable solutions in order to guarantee safe and continuous supply of water [5,6,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,11,12,13,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,18,19,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 Mm3/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−2 m2/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.
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 m3/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. 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 m3/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 m3/s during extreme rainfall events down to few m3/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.
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.

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,32,33,34,35,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).
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.
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 2H/1H and 18O/16O are expressed as δ notation [δ = ((Rsample/Rstandard) − 1) × 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 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 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:
(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 riverbed profile prior to the weir construction;
(b)
Simulation of the aquifer river exchange prior to any riverbed modification considering aquifer exploitation for 0.350 m3/s;
(c)
Simulation of the aquifer river exchange considering the Sant’Alessio weir in operation and aquifer exploitation for 0.350 m3/s;
(d)
Simulation of the aquifer river exchange considering the Sant’Alessio weir in operation and aquifer exploitation for 0.430 m3/s.

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.
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 × 10−3 m/s and 3.3 × 10−4 m/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.
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 (Figure 1 and Figure 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.
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 m3/s. Construction of the first set of wells and their operation at an overall abstraction rate of 0.350 m3/s would have caused an estimated recharge of about 0.490 m3/s, hence inducing in the aquifer around 0.340 m3/s. We then estimated in about 0.030 m3/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 Mm3. Finally, increasing groundwater withdrawal up to 0.430 m3/s (from initial 0.350 m3/s) induces an additional recharge of about 0.09 m3/s.
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.
As far as water hydrochemistry, according to the Piper classification diagrams (Figure 11 bottom), the samples mainly exhibit a Na-K-HCO3 hydrochemical facies with a clear trend towards the Ca-Mg-SO4 facies, perfectly aligned on a mixing line from the point 221 to the Serchio river (see Figure 5 and Figure 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 Cl- concentrations are more uniform in the groundwater.
The δ18O and δ2H values of water are very uniform in most of the samples and they have isotopic values ranging between −8.2‰ and −4.9‰ for δ18O, and −31.7‰ and −49‰ for δ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 Figure 7 and Figure 11) is the least depleted (in light isotopes, i.e., 16O and 1H), 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 δ18O and δ2H 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.
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).
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.

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 Mm3/year, which is 85% of the total amount of water abstracted in a year (~16 Mm3). 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

  1. WHO; UNICEF. Progress on Drinking Water, Sanitation and Hygiene: Joint Monitoring Programme 2017 Update and SDG Baselines. Available online: http://www.who.int/water_sanitation_health/publications/jmp-2017/en/ (accessed on 22 October 2020).
  2. Meffe, R.; de Bustamante, I. Emerging organic contaminants in surface water and groundwater: A first overview of the situation in Italy. Sci. Total Environ. 2014, 481, 280–295. [Google Scholar] [CrossRef] [PubMed]
  3. Paldor, A.; Shalev, E.; Katz, O.; Aharonov, E. Dynamics of saltwater intrusion and submarine groundwater discharge in confined coastal aquifers: A case study in northern Israel. Hydrogeol. J. 2019, 27, 1611–1625. [Google Scholar] [CrossRef]
  4. Wanner, P. Plastic in agricultural soils—A global risk for groundwater systems and drinking water supplies?—A review. Chemosphere 2021, 264, 128453. [Google Scholar] [CrossRef] [PubMed]
  5. Dillon, P.; Stuyfzand, P.; Grischek, T. Sixty years of global progress in managed aquifer recharge. Hydrogeol. J. 2019, 27, 1–30. [Google Scholar] [CrossRef] [Green Version]
  6. Dillon, P.; Fernández Escalante, E.; Megdal, S.B.; Massmann, G. Managed Aquifer Recharge for Water Resilience. Water 2020, 12, 1846. [Google Scholar] [CrossRef]
  7. Fernandez Escalante, E.; Henao Casas, J.D.; Vidal Medeiros, A.M.; San Sebastián Sauto, J. Regulations and guidelines on water quality requirements for Managed Aquifer Recharge. International comparison. Acque Sotter. Ital. J. Groundw. 2020, 9. [Google Scholar] [CrossRef]
  8. Dillon, P. Future management of aquifer recharge. Hydrogeol. J. 2005, 13, 313–316. [Google Scholar] [CrossRef]
  9. Stefan, C.; Ansems, N. Web-based global inventory of managed aquifer recharge applications. Sustain. Water Resour. Manag. 2018, 4, 153–162. [Google Scholar] [CrossRef] [Green Version]
  10. Umar, D.A.; Ramli, M.F.; Aris, A.Z.; Sulaiman, W.N.A.; Kura, N.U.; Tukur, A.I. An overview assessment of the effectiveness and global popularity of some methods used in measuring riverbank filtration. J. Hydrol. 2017, 550, 497–515. [Google Scholar] [CrossRef]
  11. Kühn, W.; Müller, U. Riverbank filtration. JAWWA 2000, 92, 60–69. [Google Scholar] [CrossRef]
  12. Hiscock, K.M.; Grischek, T. Attenuation of groundwater pollution by bank filtration. J. Hydrol. 2002, 266, 139–144. [Google Scholar] [CrossRef] [Green Version]
  13. Hunt, H.; Schubert, J.; Ray, C. Conceptual Design of Riverbank Filtration Systems. In Riverbank Filtration: Improving Source Water Quality; Springer: Berlin/Heidelberg, Germany, 2002; pp. 19–27. [Google Scholar]
  14. Ray, C. Riverbank Filtration Concepts and Applicability to Desert Environments. In Riverbank Filtration for Water Security in Desert Countries; Springer: Berlin/Heidelberg, Germany, 2011; pp. 1–4. [Google Scholar]
  15. Shankar, V.; Eckert, P.; Ojha, C. Transient three-dimensional modeling of riverbank filtration at Grind well field, Germany. Hydrogeol. J. 2009, 17, 321–326. [Google Scholar] [CrossRef]
  16. Ray, C.; Grischek, T.; Schubert, J.; Wang, Z.; Speth, T.F. A perspective of riverbank filtration. JAWWA 2002, 94, 149–160. [Google Scholar] [CrossRef]
  17. Hoppe-Jones, C.; Oldham, G.; Drewes, J.E. Attenuation of total organic carbon and unregulated trace organic chemicals in U.S. riverbank filtration systems. Water Res. 2010, 44, 4643–4659. [Google Scholar] [CrossRef] [PubMed]
  18. Massmann, G.; Dünnbier, U.; Heberer, T.; Taute, T. Chemosphere Behaviour and redox sensitivity of pharmaceutical residues during bank filtration–Investigation of residues of phenazone-type analgesics. Chemosphere 2008, 71, 1476–1485. [Google Scholar] [CrossRef] [PubMed]
  19. Sprenger, C.; Lorenzen, G.; Hülshoff, I. Vulnerability of bank filtration systems to climate change. Sci. Total Environ. 2011, 409, 655–663. [Google Scholar] [CrossRef]
  20. Ahmed, A.K.A.; Marhaba, T.F. Review on river bank filtration as an in situ water treatment process. Clean Technol. Environ. Policy 2017, 19, 349–359. [Google Scholar] [CrossRef]
  21. Barthel, R.; Banzhaf, S. Groundwater and Surface Water Interaction at the Regional-scale—A Review with Focus on Regional Integrated Models. Water Resour. Manag. 2016, 30, 1–32. [Google Scholar] [CrossRef] [Green Version]
  22. Cook, P.G. Quantifying river gain and loss at regional scales (Review). J. Hydrol. 2015, 531, 749–758. [Google Scholar] [CrossRef]
  23. Sottani, A.; Vielmo, A. Groundwater conservation and monitoring activities in the middle Brenta River plain (Veneto Region, Northern Italy): Preliminary results about aquifer recharge. Acque Sotter. Ital. J. Groundw. 2014, 3. [Google Scholar] [CrossRef]
  24. Wang, W.-S.; Oswald, S.E.; Gräff, T.; Lensing, H.-J.; Liu, T.; Strasser, D.; Munz, M. Impact of river reconstruction on groundwater flow during bank filtration assessed by transient three-dimensional modelling of flow and heat transport. Hydrogeol. J. 2020, 28, 723–743. [Google Scholar] [CrossRef] [Green Version]
  25. Barbagli, A.; Jensen, B.N.; Raza, M.; Schueth, C.; Rossetto, R. Assessment of soil buffer capacity on nutrients and pharmaceuticals in nature-based solution applications. Environ. Sci. Pollut. Res. 2019, 26, 759–774. [Google Scholar] [CrossRef] [PubMed]
  26. Rossetto, R.; Barbagli, A.; Borsi, I.; Mazzanti, G.; Vienken, T.; Bonari, E. Site investigation and design of the monitoring system at the Sant’Alessio Induced RiverBank Filtration plant (Lucca, Italy). Rend. Online Soc. Geol. Ital. 2015, 35, 248–251. [Google Scholar] [CrossRef]
  27. Studio Nolledi. Indagini Idrogeologiche Finalizzate All’individuazione Dell’area di Salvaguardia del Campo Pozzi di sant’Alessio—Lucca Art. 21 del D.Lgs. 152/99. Prima Fase Individuazione e Caratterizzazione Dell’area di Rispetto. 2003; (Comune di Lucca Unpublished Report); Available online: file:///C:/Users/mdpi/AppData/Local/Temp/D39.pdf (accessed on 7 December 2020).
  28. Rossetto, R.; Barbagli, A.; De Filippis, G.; Marchina, C.; Mazzanti, G.; De Caterini, A. Case Study 19: The Serchio River Bank Filtration for Drinking Water Supply in Sant’Alessio area of Lucca, Italy. In Managing Aquifer Recharge: A Showcase for Resilience and Sustainability; Zheng, Y., Ross, A., Villholth, K.G., Dillon, P., Eds.; A Unesco-Iah-Gripp Publication, 2002; (to be Advised); Available online: https://recharge.iah.org/unesco-exemplary-mar-projects-booklet (accessed on 8 December 2020).
  29. Rossetto, R.; Bockelmann-Evans, B. Modellazione numerica del flusso e del trasporto di soluti ai fini dell’investigazione dei processi di trasporto dell’erbicida terbutilazina nel sistema acquifero della pianura di S. Alessio (Lucca). G. Geol. Appl. 2007, 5, 29–44. [Google Scholar]
  30. EPA. Expedited Site Assessment Tools for underground Storage Tank Sites. A Guide for Regulators; U.S. Government Printing Office: Washington, DC, USA, 1997.
  31. Butler, J.J. Hydrogeological Methods for estimation of Spatial Variations in Hydraulic Conductivity. In Hydrogeophysics; Rubin, Y., Hubbard, S.S., Eds.; Springer: Berlin/Heidelberg, Germany, 2005; pp. 23–58. [Google Scholar]
  32. Dietrich, P.; Leven, C. Direct Push-Technologies. In Groundwater Geophysics; Kirsch, R., Ed.; Springer: Berlin/Heidelberg, Germany, 2006; pp. 321–340. [Google Scholar]
  33. McCall, W.; Nielsen, D.M.; Farrington, S.P.; Christy, T.M. Use of Direct-Push Technologies in Environmental Site Characterization and Ground-Water Monitoring. In Practical Handbook of Environmental Site Characterization and Ground-Water Monitoring; Nielsen, D.M., Ed.; CRC Press Taylor and Francis Group: Boca Raton, FL, USA, 2006; pp. 345–471. [Google Scholar]
  34. Liu, G.S.; Butler, J.J.; Reboulet, E.; Knobbe, S. Hydraulic conductivity profiling with direct push methods. Grundwasser 2012, 17, 19–29. [Google Scholar] [CrossRef]
  35. Vienken, T.; Leven, C.; Dietrich, P. Use of CPT and other direct push methods for (hydro-) stratigraphic aquifer characterization—A field study. Can. Geotech. J. 2012, 49, 197–206. [Google Scholar] [CrossRef]
  36. Vienken, T.; Kreck, M.; Hausmann, J.; Werban, U.; Dietrich, P. Innovative strategies for high resolution site characterization: Application to a flood plain. Ital. J. Groundw. 2014, 138, 7–14. [Google Scholar] [CrossRef] [Green Version]
  37. Dietrich, P.; Butler, J.J.; Faiss, K. A rapid method for hydraulic profiling in unconsolidated formations. Ground Water 2008, 46, 323–328. [Google Scholar] [CrossRef]
  38. Sellwood, S.M.; Healey, J.M.; Birk, S.; Butler, J.J. Direct-push hydrostratigraphic profiling: Coupling electrical logging and slug tests. Ground Water 2005, 43, 19–29. [Google Scholar] [CrossRef]
  39. Vienken, T.; Huber, E.; Kreck, M.; Huggenberger, P.; Dietrich, P. How to chase a tracer—Combining conventional salt tracer testing and direct push electrical conductivity profiling for enhanced aquifer characterization. Adv. Water Resour. 2017, 99, 60–66. [Google Scholar] [CrossRef]
  40. Marchina, C.; Zuecco, G.; Chiogna, G.; Bianchini, G.; Carturan, L.; Comiti, F.; Engel, M.; Natali, C.; Borga, M.; Penna, D. Alternative methods to determine the δ2H-δ18O relationship: An application to different water types. J. Hydrol. 2020, 587, 124951. [Google Scholar] [CrossRef]
  41. De Filippis, G.; Borsi, I.; Foglia, L.; Cannata, M.; Velasco Mansilla, V.; Vasquez-Suñe, E.; Ghetta, M.; Rossetto, R. Software tools for sustainable water resources management: The GIS-integrated FREEWAT platform. Rend. Online Soc. Geol. 2017, 42, 59–61. [Google Scholar] [CrossRef] [Green Version]
  42. Foglia, L.; Borsi, I.; Mehl, S.; De Filippis, G.; Cannata, M.; Vasquez-Sune, E.; Criollo, R.; Rossetto, R. Freewat, A Free and Open Source, GIS-Integrated, Hydrological Modeling Platform. Groundwater 2018, 56, 521–523. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Harbaugh, A.W. MODFLOW-2005, The U.S. Geological Survey Modular Ground-Water Model—The Ground-Water Flow Process; Techniques and Methods 6-A16; U.S. Geological Survey: Reston, VA, USA, 2005. [CrossRef] [Green Version]
  44. Rossetto, R.; Borsi, I.; Schifani, C.; Bonari, E.; Mogorovich, P.; Primicerio, M. SID&GRID: Integrating hydrological modeling in GIS environment. Rend. Online Soc. Geol. Ital. 2013, 24, 282–283. [Google Scholar]
  45. Rossetto, R.; Borsi, I.; Foglia, L. Freewat: Free and open source software tools for WATer resource management. Rend. Online Soc. Geol. Ital. 2015, 35, 252–255. [Google Scholar] [CrossRef] [Green Version]
  46. Criollo, R.; Velasco, V.; Nardi, A.; Vries, L.M.; Riera, C.; Scheiber, L.; Jurado, A.; Brouyère, S.; Pujades, E.; Rossetto, R.; et al. AkvaGIS: An open source tool for water quantity and quality management. Comput. Geosci. 2020, 127, 123–132. [Google Scholar] [CrossRef]
  47. Cannata, M.; Neumann, J.; Rossetto, R. Open source GIS platform for water resource modelling: FREEWAT approach in the Lugano Lake. Spat. Inf. Res. 2018, 26, 241–251. [Google Scholar] [CrossRef] [Green Version]
  48. De Filippis, G.; Pouliaris, C.; Kahuda, D.; Vasile, T.A.; Manea, V.A.; Zaun, F.; Panteleit, B.; Dadaser-Celik, F.; Positano, P.; Nannucci, M.S.; et al. Spatial Data Management and Numerical Modelling: Demonstrating the Application of the QGIS-Integrated FREEWAT Platform at 13 Case Studies for Tackling Groundwater Resource Management. Water 2020, 12, 41. [Google Scholar] [CrossRef] [Green Version]
  49. Joodavi, A.; Izady, A.; Karbasi Maroof, M.T.; Majidi, M.; Rossetto, R. Deriving optimal operational policies for off-stream man-made reservoir considering conjunctive use of surface- and groundwater at the Bar dam reservoir (Iran). J. Hydrol. Reg. Stud. 2020, 31, 100725. [Google Scholar] [CrossRef]
  50. La Vigna, F.; Hill, M.C.; Rossetto, R.; Mazza, R. Parameterization, sensitivity analysis, and inversion: An investigation using groundwater modeling of the surface-mined Tivoli-Guidonia basin (Metropolitan City of Rome, Italy). Hydrogeol. J. 2016, 24, 1423–1441. [Google Scholar] [CrossRef]
  51. Rossetto, R.; Barbagli, A.; De Filippis, G.; Marchina, C.; Di Bartolo, S.; Bonari, E.; Sabbatini, T.; Triana, F.; Tozzini, C.; Ercoli, L.; et al. Deliverable 8.4. Report on the Induced RiverBank Filtration MAR Plant at Sant’Alessio (Lucca, Italy). FP7 MARSOL (Demonstrating Managed Aquifer Recharge as a Solution to Water Scarcity and Drought). Available online: http://marsol.eu/ (accessed on 28 November 2020).
  52. Butler, J.J., Jr.; Garnett, E.J.; Healey, J.M. Analysis of Slug Tests in Formations of High Hydraulic Conductivity. Groundwater 2003, 41, 620–631. [Google Scholar] [CrossRef] [PubMed]
  53. Barbagli, A. Analysis of Water-Soil Interaction in Drainage Water Phyto-Treatment and in Aquifer Recharge Schemes. Ph.D. Thesis, Scuola Superiore Sant’Anna, Pisa, Italy. Available online: https://dta.santannapisa.it (accessed on 28 November 2020).
  54. Rozanski, K.; Araguás-Araguás, L.; Gonfiantini, R. Isotopic Patterns in Modern Global Precipitation, in Climate Change in Continental Isotopic Records; Swart, P.K., Lohmann, K.C., Mckenzie, J., Savin, S., Eds.; American Geophysical Union: Washington, DC, USA, 1993. [Google Scholar] [CrossRef]
  55. Longinelli, A.; Selmo, E. Isotopic composition of precipitation in Italy: A first overall map. J. Hydrol. 2003, 270, 75–88. [Google Scholar]
  56. La Ruffa, G.; Panichi, C. Caratterizzazione Chimico-Isotopica Delle Acque Fluviali: Il Caso Del Fiume Arno; Istituti Editoriali Poligrafici Internazionali: Pisa, Italy, 2000; p. 101. [Google Scholar]
  57. Rodríguez-Escales, P.; Canelles, A.; Sanchez-Vila, X.; Folch, A.; Kurtzman, D.; Rossetto, R.; Fernández-Escalante, E.; Lobo-Ferreira, J.P.; Sapiano, M.; San-Sebastián, J.; et al. A risk assessment methodology to evaluate the risk failure of managed aquifer recharge in the Mediterranean Basin. Hydrol. Earth Syst. Sci. 2018, 22, 3213–3227. [Google Scholar] [CrossRef] [Green Version]
  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).
Figure 1. Layout of the Serchio River riverbank filtration scheme in Lucca, Italy. Coordinates are in WGS84-UTM zone 32N (EPSGr: 32632).
Figure 1. Layout of the Serchio River riverbank filtration scheme in Lucca, Italy. Coordinates are in WGS84-UTM zone 32N (EPSGr: 32632).
Hydrology 07 00096 g001
Figure 2. The Sant’Alessio weir on the Serchio River.
Figure 2. The Sant’Alessio weir on the Serchio River.
Hydrology 07 00096 g002
Figure 3. Simplified stratigraphic cross section (top) and detailed stratigraphic logs around the Serchio River RBF area (bottom).
Figure 3. Simplified stratigraphic cross section (top) and detailed stratigraphic logs around the Serchio River RBF area (bottom).
Hydrology 07 00096 g003
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.
Hydrology 07 00096 g004
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 07 00096 g005
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.
Hydrology 07 00096 g006
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.
Hydrology 07 00096 g007
Figure 8. Direct Push profiling results (DPIL and Electrical Conductivity measurements).
Figure 8. Direct Push profiling results (DPIL and Electrical Conductivity measurements).
Hydrology 07 00096 g008
Figure 9. Simulated hydraulic head at the end of May 2016.
Figure 9. Simulated hydraulic head at the end of May 2016.
Hydrology 07 00096 g009
Figure 10. Simulated difference in groundwater head between pre- and post-operam of the Sant’Alessio weir.
Figure 10. Simulated difference in groundwater head between pre- and post-operam of the Sant’Alessio weir.
Hydrology 07 00096 g010
Figure 11. Piper diagram (bottom), Cl 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 vs Br- binary plot (top right) and δ18O vs δ2H binary plot (top left) of the samples collected during the dry season of 2015.
Hydrology 07 00096 g011
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.
Hydrology 07 00096 g012
Table 1. Simulated water budget for the Sant’Alessio aquifer for the hydrologic year October 2015 to September 2016.
Table 1. Simulated water budget for the Sant’Alessio aquifer for the hydrologic year October 2015 to September 2016.
Inflow TermsOutflow Terms
Cumulative Volume (m3)% Over the TotalCumulative Volume (m3)% Over the Total
Storage846Negligible938Negligible
Inflow from the Monte S.Quirico—Carignano hills207,0890.315,640,34126.4
Rainfall infiltration673,2871.1658Negligible
River leakage58,154,74498.238,439,23264.9
Southern boundary of the domain160,2140.35,109,2728.6
TOTAL59,196,180100.059,190,432100.0
Table 2. Simulated changes in river recharge from pristine conditions to completion of the Serchio River RBF managed aquifer recharge scheme.
Table 2. Simulated changes in river recharge from pristine conditions to completion of the Serchio River RBF managed aquifer recharge scheme.
No-Weir/No-WellsNo-Weir/Average Pumping (0.350 m3/s)Change in Recharge Rate (m3/s)
Net Aquifer Recharge (m3/s)0.1510.4880.337
No-Weir/Average Pumping (0.350 m3/s)Weir/Average Pumping (0.350 m3/s)Change in Recharge Rate (m3/s)
Net Aquifer Recharge (m3/s)0.4880.5210.033
Weir/Average Pumping (0.350 m3/s)Weir/High Pumping (0.430 m3/s)Change in Recharge Rate (m3/s)
Net Aquifer Recharge (m3/s)0.5210.6090.088
Table 3. Average values during dry and wet season of the presented analytes.
Table 3. Average values during dry and wet season of the presented analytes.
Dry SeasonWet Season
Sampling PointNa+
mg/L
K+
mg/L
Ca2+
mg/L
Mg2+
mg/L
HCO3
mg/L
Cl
mg/L
SO42−
mg/L
Br
mg/L
Li
µg/L
δ18O
δD
Na+
mg/L
K+
mg/L
Ca2+
mg/L
Mg2+
mg/L
HCO3
mg/L
Cl
mg/L
SO42−
mg/L
Br
mg/L
Li
µg/L
δ18O
δD
9313.972.2371.429.4723416.8355.790.067.33--16.192.09104.9011.3419129.04103.310.167.63--
9414.662.8679.8110.4122815.1047.700.0710.26-------------
22113.353.8052.674.5924920.8514.670.110.77−4.9−31.714.453.8569.845.4227022.6710.480.140.74−5.0−30.6
C_215.671.6068.329.5416426.1597.050.1211.66--19.051.6283.5510.5315027.0098.880.1313.30−6.7−46.6
C_519.301.5569.459.4415531.54101.770.1112.29−7.9−47.617.821.5574.939.0714427.38103.870.1313.13−6.4−43.7
MAR_1015.162.1863.3010.0932024.4291.970.087.93−6.3−43.321.392.81117.2713.8422428.9380.760.1010.32−6.8−42.3
P_117.731.6571.889.7316229.58116.050.1413.21--22.631.4286.1812.0114042.87129.960.2815.22--
PR_01a17.771.9874.069.2416233.49101.050.1613.25−8.2−48.116.621.3779.3910.2612027.80102.040.1812.31−7.2−44.0
PR_01b15.842.3661.098.408717.3060.800.079.29−7.4−45.7-----------
PR_01c16.061.6573.148.6015919.9562.980.0712.22−7.8−47.320.401.3783.2910.5317430.66110.710.2112.02−7.0−44.0
PR_02a13.761.4859.218.5220920.2562.650.1310.91−8.0−47.119.771.4880.8710.7915532.91126.160.2014.78−7.2−43.7
PR_02b15.771.7067.909.0813325.6999.800.1213.95−8.0−48.221.021.7878.7511.6312733.23125.070.3413.89--
PR_03b13.802.3468.657.6114425.5067.200.1611.51−7.4−46.5-----------
PR_04a14.701.7457.407.6618121.6769.990.119.98−7.6−46.918.971.6479.0410.3714332.99119.660.179.99--
PR_04b14.411.5461.048.5115924.0572.620.1211.11−7.6−49.020.141.6583.5510.7815331.09122.800.1413.02−7.2−46.1
PR_04c14.751.3659.658.0512327.2481.390.1711.18−7.5−46.7-----------
PR_05a13.891.3452.228.0313520.7258.030.0710.10--19.431.5177.1810.3823131.89121.680.1812.20−7.5−46.6
PR_05b14.031.2459.178.9117818.8868.470.058.08--19.981.9987.2910.9716331.12118.680.3313.42--
PR_06a11.111.1848.917.0116025.3175.400.108.86--19.541.9176.6810.8513331.73125.620.1715.10−7.1−44.6
PR_06b11.401.8652.517.6812923.5083.930.0811.40-------------
PR_0720.952.6731.0616.4418045.4570.11<LOD1.72--24.923.56137.7520.0226439.4043.000.092.99−7.3−43.4
PR_0817.101.5758.0410.0419024.2659.84<LOD7.60--17.382.50110.0611.9919632.95100.600.189.20−6.8−42.4
PR_0913.493.4655.249.8819016.5659.66<LOD4.71--16.563.84102.9611.0921427.0395.950.167.50−7.1−43.0
Freddanella Ditch40.026.3962.356.77294112.2157.110.115.43--47.044.4592.458.4320594.4241.950.116.14--
Serchio River20.921.5964.959.3514429.0697.600.1911.27−8.2−48.115.151.3558.698.1015323.0367.270.1112.30−6.5−44.2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rossetto, R.; Barbagli, A.; De Filippis, G.; Marchina, C.; Vienken, T.; Mazzanti, G. Importance of the Induced Recharge Term in Riverbank Filtration: Hydrodynamics, Hydrochemical, and Numerical Modelling Investigations. Hydrology 2020, 7, 96. https://0-doi-org.brum.beds.ac.uk/10.3390/hydrology7040096

AMA Style

Rossetto R, Barbagli A, De Filippis G, Marchina C, Vienken T, Mazzanti G. Importance of the Induced Recharge Term in Riverbank Filtration: Hydrodynamics, Hydrochemical, and Numerical Modelling Investigations. Hydrology. 2020; 7(4):96. https://0-doi-org.brum.beds.ac.uk/10.3390/hydrology7040096

Chicago/Turabian Style

Rossetto, Rudy, Alessio Barbagli, Giovanna De Filippis, Chiara Marchina, Thomas Vienken, and Giorgio Mazzanti. 2020. "Importance of the Induced Recharge Term in Riverbank Filtration: Hydrodynamics, Hydrochemical, and Numerical Modelling Investigations" Hydrology 7, no. 4: 96. https://0-doi-org.brum.beds.ac.uk/10.3390/hydrology7040096

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop