water
Article Karst Recharge Areas Identified by Combined Application of Isotopes and Hydrogeological Budget
Silvia Iacurto , Gerardo Grelle, Francesco Maria De Filippi and Giuseppe Sappa *
Department of Civil, Constructional and Environmental Engineering (DICEA), Sapienza University of Rome, 00184 Rome, Italy; [email protected] (S.I.); [email protected] (G.G.); francescomaria.defi[email protected] (F.M.D.F.) * Correspondence: [email protected]; Tel.: +39-0644-585-010 or +39-3452-808-882
Abstract: The identification of recharge areas in karst aquifers allows us to perform sustainable management of these groundwater resources. Stable isotopes (δ18O and δ2H) have been largely used to provide information about recharge elevation in many mountainous regions. In this paper, an improved version of a recent “isotope-driven model”, for the identification of recharge areas, was applied to Capodacqua di Spigno Spring (south of the Latium region). The model upgrade consists of a preliminary check procedure to estimate the degree of influence of the rainfall’s isotopic variability on the spring water. This additional procedure gives us an indication of the reliability of the model and its applicability conditions. Moreover, the dataset of the spring was updated to analyze the degree of reliability of the isotope-driven model. The purpose of this study was to combine the previously mentioned isotope-driven model with hydrogeological tools. A quantitative study of the basin, based on the estimation of the average monthly infiltration volume, was performed by using
the inverse hydrogeological water budget. In this way, the qualitative model for the recharge areas’ estimation was validated by a quantitative hydrogeological tool. Both models show that, for karst
Citation: Iacurto, S.; Grelle, G.; De mountain basins, the recharge areas decrease as the average recharge elevations increase, including Filippi, F.M.; Sappa, G. Karst areas at high altitudes. Recharge Areas Identified by Combined Application of Isotopes Keywords: recharge area; karst aquifer; stable isotope; inverse hydrogeological water budget; and Hydrogeological Budget. Water isotopic tracers; spring discharge; rainfall isotopic variability 2021, 13, 1965. https://doi.org/ 10.3390/w13141965
Academic Editor: Juraj Parajka 1. Introduction Karst groundwater is a primary source of freshwater in many regions around the Received: 19 May 2021 world. About 25% of the global population is partly or entirely supplied by freshwater Accepted: 13 July 2021 Published: 17 July 2021 from karst. Most of the largest springs in the world are located in karst areas [1]. It was estimated that the percentage of karst water consumers in 2016 was 9.2% of the world’s
Publisher’s Note: MDPI stays neutral population [2]. Indeed, for a large part of the Mediterranean region, karst springs play an with regard to jurisdictional claims in essential role for water supply. Karst springs provide fresh and high-quality water, which published maps and institutional affil- has been an important resource for human development in this region since antiquity [3]. iations. Specifically, in Italy, the karst carbonate aquifers of the central Apennines represent the largest groundwater reservoir [3,4]. Therefore, for the protection and sustainable management of groundwater resources, it is critical to study the recharge–discharge mechanism of a karst aquifer and determine the recharge areas [5]. These aquifers are very complex systems, as each one has its own Copyright: © 2021 by the authors. Licensee MDPI, Basel, Switzerland. distinctive characteristics [6,7]. The high spatial heterogeneity of a karst aquifer, often This article is an open access article characterized by poor available data, requires specific investigation techniques such as distributed under the terms and environmental tracers. At the common temperature of groundwater, there is no sensitive conditions of the Creative Commons isotope exchange between the groundwater and minerals. Therefore, stable isotopes such as 2 18 Attribution (CC BY) license (https:// H and O can act as conservative tracers [8]. In this framework, the isotopic tracers play a creativecommons.org/licenses/by/ key role. As the basic conceptual model, it is possible to claim that the isotopic composition 4.0/). and concentration values are affected by a number of environmental variables, which
Water 2021, 13, 1965. https://doi.org/10.3390/w13141965 https://www.mdpi.com/journal/water Water 2021, 13, 1965 2 of 24
affect evapotranspiration output. Starting from Dansgaard (1964) [9], many attempts were undertaken to include the relationship of the isotopic composition of monthly precipitation with the latitude, altitude and amount of precipitation, the distance from the coastline, and the air temperature [10]. Some of these effects have local variability and were represented by the Local Meteoric Water Line (LMWL) [11,12]; in this way, light changing of the 2H to 18O ratio also identifies a specific composition marker for an assigned meteorological region. In addition to this regional behavior, other local effects contribute to change in the rainfall’s isotopic concentration. Among these ones, the altitude–topographic effect seems to have a relevant role in mountainous areas. This effect is related to the evapotranspiration of air that induces an enrichment of the heavy-isotopic content during precipitation, up to the time the rainwater reaches the soil surface. The enrichment is strictly due to the air temperature gradient and the travel time from the clouds to the topographic surface. Therefore, in mountainous environments, the isotopic concentration of the rainwater impacting the soil is strictly linked to the altitude at which it happens [12]. Stable isotopes may be used as conservative tracers to understand water dynamics, e.g., groundwater recharge areas, its source and movement [13,14]. Water isotopes are a natural part of the water molecule, and so these can be regarded as the most ideal tracer. Combined studies on natural tracers, such as isotopes, and spring discharge may provide integrated knowledge regarding the key characteristics of the karst groundwater system [15]. Oxygen-18 and deuterium have been largely used to determine groundwater recharge areas, with their elevations, transit times and water flow paths [16–20]. Considering the local vertical isotopic gradient, the isotopic values of groundwater can be used to identify the average recharge area elevation of aquifers [17]. This average elevation can be estimated starting from the δ18O and δ2H values of the water sampled, at a known altitude, in a pluviometric station located near the study area. Identifying a specific elevation at which groundwater has infiltrated is a technique frequently applied in isotope hydrology [10,17]. This, however, assumes that the spring sample is the product of rainwater infiltration at a single elevation, which is a strong simplification. In fact, the elevation identified by the isotope is always an average recharge elevation. Based on this consideration, an isotope-driven model has been set up [21], which estimates, for any discharge value of a spring, the percentage of the topographic area involved in the aquifer recharge and its spatial distribution. Specifically, the model identifies the main recharge areas according to the oxygen and hydrogen isotopes’ (δ18O and δ2H) rates of variation, which are detected in groundwater samples taken at the spring. Considering this framework, the identification of recharge areas has a relevant role for a proper water governance, with the aim of safeguarding spring water quality.
Objectives The present study focuses on the identification of recharge areas for karst springs by the combined application of stable isotopes (δ18O and δ2H) and hydrogeological budget. The latter is performed using the empirical method for estimating potential evapotranspi- ration (PET) proposed by Thornthwaite (1948) [22], applied to the inverse hydrogeological budget method (IHBM) [23–25]. Specifically, in the present study, the previously men- tioned isotope-driven model [21] has been improved. The model extension consists of a preliminary check procedure to estimate the influence degree of seasonality combined with amount effect of rainfall (input) on spring water (output). This additional procedure aims to give an overview of the favorable conditions to the applicability of the model. Moreover, in this study, the updating of the previous dataset, regarding Capodacqua di Spigno Spring (in the southern part of the Latium region), made it possible to verify the consistency and soundness of the model. Water 2021, 13, 1965 3 of 24
2. Study Area Characterization 2.1. Geological and Hydrogeological Setting of the Study Area Capodacqua di Spigno Spring, one of the main groundwater outlets in the southeast Latium region, is sited in the province of Latina, about 4 km from the Tyrrhenian Sea [21]. This spring, located at the bottom of Petrella Mountain, has been known since Roman times. In 72 BC, the emperor Vespasian built an 11 km long aqueduct that collected spring waters to supply the city of Minturno (LT). The current aqueduct of Capodacqua di Spigno Spring feeds the larger area including Spigno Saturnia, Formia and Esperia municipalities, in Latina province [26] (Figure1). The aquifer’s altitude ranges from about 100 m a.s.l. to about 1500 m a.s.l., with an average value of 880 m a.s.l; the limits of the aquifer are assumed considering the watershed topographic pattern. The altitude of the spring, instead, is of 35 m a.s.l. This spring is located in the karst coastal region of the western Aurunci Mountains, which, together with the Lepini and the Ausoni Mountains, belong to the Pre-Apennines of Latium and form the carbonatic platform of the Volsci Ridge, separated from the Apennine ridge by the Latina Valley [21]. The spring water comes out from the permeable limestone of La Civita Mountain and Castello Mountain and flows above the Upper Miocene clays, at the lowest point of the limestone–clay contact [27]. The water table map described by previous studies [28] shows that the main direction of the aquifer is generally south–south east. Differently, there are also important local flows along the two faults, which delimit the carbonate series outcropping, especially the one to the east, which seems to represent the main water conduit towards the spring. More detailed geological information on the study area can be derived from previous studies [28]. The stratigraphic succession of the soils, from the bottom upwards, is as follows:. • Carbonate Series (CS1): limestone—dolomite series that forms the western Aurunci Mountains (stratigraphic alternation Palaeocene–Upper Lias; sedimentary marine fa- cies); • Kaolinitic Clays (KC): consisting of a polychrome clay–silty blanket (marginal marine siltstone facies); • Clutch Breach (Br); • Miocene Series (MS): formed by the grouping of two hydrogeologically similar forma- tions (marine to continental–marine transitional facies); • Deposits of the lower Pliocene (p) (continental–marine facies); • Quaternary deposits (q): river and lake deposits and groundwater debris (continen- tal facies). The Carbonate Series (CS1), present in the area of the Capodacqua di Spigno basin, is characterized by limestones and dolomites (Palaeocene—Upper Lias) alternations. As found in the Geological Map of Italy (scale 1:100,000, sheet no. 171 “Gaeta”) [29] referring to the area under study, identified by the Civita Mountain, the aforementioned Carbon- ate Series (CS1) is divided into two formations: a younger Carbonate Series formed in the Cretaceous period and an older one formed during the Jurassic period. Therefore, in the geological stratigraphic succession, the Carbonate Series (CS1), Jurassic carbon- ate Series (JCS) and Cretaceous Carbonate Series (CCS) were taken into consideration. Figure1 depicts a sketch of the structural setting of the Volsci chain. This scheme has been modified starting from that proposed by Baldi et al. in 2005 [28]. The new interpretation of the geology and structural characteristics of the area came from the study of both the Geological Map of Italy (scale 1:100,000, sheet no. 171 “Gaeta”) [29] and the geophysical survey (tomography) performed in the area [28]. In this regard, the geological section AA’ is shown in Figure2.
2.2. Climate Framework Regarding the climate assessment of the study area, it has been evaluated on the basis of monthly rainfall time series collected in the “Esperia” rainfall gauge, the closest to the recharge area of the spring, in the northern part of the hydrogeological basin. However, Water 2021, 13, 1965 4 of 24
rainfall values from 1994 to 2011 are not available for this pluviometric station. Therefore, the rainfall behavior for the longest time series available, i.e., from 1959 to 1994, was studied. The mean annual precipitation (MAP) of the time series ranges from 811 to 2103 mm/year in the period 1959–1994. The rainfall seasonal trend (Figure3) shows that the most important contribution is due to the rainiest season (from November to February), with a trend inversion in the last eight years (2012–2020), where a deficit has Water 2021, 13, 1965 4 of 26 been registered, mostly due to a reduced contribution of rainfall in January and a positive increase in February and March.
FigureFigure 1. 1.Structural Structural map map and and main main groundwater groundwater flow flow direction direction of Capodacquaof Capodacqua di Spignodi Spigno Spring Spring (up ); excerpt(up); excerpt of the geologicalof the geological map and map positioning and positionin of theg geologicalof the geological section section AA’ represented AA’ represented in Figure in 3 (downFigure)—based 3 (down)—based on Baldi et on al., Baldi 2005 et [28 al.,]. 2005 [28]. Water 2021,, 13,, 19651965 5 of 26 24 Water 2021, 13, 1965 5 of 26
Figure 2. Geological sectionsection AA’AA’ (trace (trace in in Figure Figure1) modified1) modified from from the studythe study of the of geologicalthe geological map ofmap Italy of [Italy29] and [29] geological and geo- logicalmapFigure edited map2. Geological by edited Baldi by et section al.,Baldi 2005 etAA’ [28al., ]:(trace 2005 the image in[28]: Figure the below image1) depictsmodified belo thew from structuraldepicts the thestudy pattern structural of the by geoelectric geological pattern by tomographymap geoelectric of Italy interpretation. [29]tomography and geo- logical map edited by Baldi et al., 2005 [28]: the image below depicts the structural pattern by geoelectric tomography interpretation.Legend: JCS is Legend: Jurassic CarbonateJCS is Jurassic Series; Carb CCSonate is CretaceousSeries; CCS Carbonate is Cretaceous Series; Carbonate KC is Kaolinitic Series; KC Clays; is Kaolinitic Br is Clutch Clays; Breach; Br is Clutchinterpretation. Breach; MSLegend: is Miocene JCS is SeriJurassices and Carb q isonate Quaternary Series; CCSdeposits. is Cretaceous Carbonate Series; KC is Kaolinitic Clays; Br is MS is Miocene Series and q is Quaternary deposits. Clutch Breach; MS is Miocene Series and q is Quaternary deposits.
Figure 3. Mean Monthly Precipitation (MMP) seasonal trend in the study area, based on [26]. FigureFigure 3.3. MeanMean MonthlyMonthly PrecipitationPrecipitation (MMP)(MMP) seasonalseasonal trendtrend inin thethe studystudy area,area, basedbased onon [[26].26].
Water 2021, 13, 1965 6 of 24
3. Materials and Methods 3.1. Data Collection In addition to the 14 samples already analyzed, whose results have been previously published [21], the isotope geochemistry laboratory of the University of Parma (Italy) tested 10 more groundwater samples, which were taken always from Capodacqua di Spigno Spring. This means that we considered 70% more samples. The IRMS (Isotope- 18 ratio mass spectrometry) continuous low-equilibration method, with CO2 for δ O and H2 for δ2H, was applied. The results of all samples collected at the spring, shown in Table1, are expressed as per mil deviations from the internationally accepted standard VSMOW (Vienna Standard Mean Ocean Water). The “¦Ä” notation is a way to express the relative abundances of oxygen (16O, 18O) and hydrogen (1H, 2H) isotopes in the water molecules [30]. The analytical error is ±0.2‰ for δ18O and ±1‰ for δ2H.
Table 1. Isotope composition of spring water (δ18O, δ2H), flow rates exploited by the local water supply (Qe), excess flow rates measured in the river (Qm) and reconstruction of the discharges supplied by Capodacqua di Spigno Spring (Qs), (new data in bold).
δ18O δ2H Date Qe (l/s) Qm (l/s) Qs (l/s) (‰ VSMOW) (‰ VSMOW) 26 April 2018 −7.22 −40.09 520 915 1435 25 May 2018 −7.24 −40.35 506 780 1286 27 June 2018 −7.29 −41.02 537 490 1027 30 July 2018 −7.32 −41.07 543 251 794 29 August 2018 −7.27 −40.94 554 159 713 1 October 2018 −7.37 −41.20 489 168 657 8 November 2018 −6.91 −37.92 454 1657 2111 18 December 2018 −6.78 −37.57 486 2164 2650 18 January 2019 −7.14 −39.80 503 829 1333 19 February 2019 −7.23 −40.29 504 1050 1555 20 March 2019 −7.26 −40.58 454 837 1291 30 April 2019 −7.23 −40.09 433 558 991 29 May 2019 −6.74 −37.15 435 2500 2935 26 June 2019 −7.37 −41.22 458 671 1129 25 July 2019 −7.36 −40.54 524 337 861 29 August 2019 −7.31 −40.47 534 102 636 27 September 2019 −7.40 −40.73 487 102 589 30 October 2019 −7.49 −41.60 466 52 518 27 November 2019 −6.43 −36.05 449 2242 2690 23 December 2019 −6.95 −39.00 425 3000 3425 28 January 2020 −6.76 −37.63 443 1440 1883 27 February 2020 −7.17 −39.80 450 784 1234 28 April 2020 −7.22 −39.86 417 595 1012 25 May 2020 −7.30 −40.54 421 550 971
During the same time period the Department of Civil, Building and Environmental Engineering (DICEA) team took spring water samples and carried out discharge measure- ments. Specifically, the authors made measurements along the river downward of the spring to determinate the excess flow rate of the spring, which has to be added to the discharge, collected by the local water supply agency Acqualatina Spa, in order to obtain the total amount of the spring outflow, Qs, as represented in Table1. Measurements were conducted by the application of traditional current-meter method [21]. The stream flow velocity was measured through a SEBA horizontal axis current-meter F1, according to EN ISO 748:2007 requirements. This activity consisted mainly of measuring the excess flow rate of the spring in a discharge channel [26]. The current meter measures velocity at a one point, so it requires determination of the mean velocity in each of the selected verticals. When possible, to measure water velocity, the “Three-Point Method” described by the “National Environmental Monitoring Standards (NEMS)” was used [31]. Water 2021, 13, 1965 7 of 24
The morphometric parameters were defined by using the digital elevation model (DEM) discretized into finite square elements (FSE), with the size of 100 × 100 m, using the open source software Q-GIS. In this way the elevation classes distribution of Capodacqua di Spigno basin was extracted from the resampled DEM. The value of the monthly supply volume of the aquifer was calculated through the inverse hydrogeological water budget method application [23–25]. For the average monthly groundwater recharge estimation (i.e., the effective infiltration), data referring to four meteorological stations near the study area, for a time series of 40 years (1959–1999) for rainfall and 6 years (2012–2018) for temperature, were used. The thermo-pluviometric Water 2021, 13, 1965 stations are Gaeta, Itri, Esperia and Santi Cosma e Damiano, which cover an area of about8 of 26
220 km2 (Figure4)[32].
FigureFigure 4. 4. Thermo-pluviometricThermo-pluviometric stations stations of of the the study study area area (latitude (latitude and and longitude longitude are are in in decimal decimal de- de- grees).grees).
3.2.3.2. Average Average Recharge Recharge Areas Areas Elevation (I-Elevation) InIn a a given hydrogeological basin, basin, groundwater sources sources may be fingerprintedfingerprinted using isotopes,isotopes, such such as as 2HH and 18O,O, which which depend depend on infiltration infiltration altitude of rainfalls. Higher infiltrationinfiltration altitude altitude will will have have precipitations, precipitations, which which are are depleted depleted in 2 inH 2andH and 18O.18 TheO. Thealtitude alti- effecttude effectis used is usedin hydrogeological in hydrogeological studies studies because because it allows it allows researchers researchers to determine to determine at whichat which altitude altitude recharge recharge areas areas are arelocated. located. However, However, some some studies studies on this on thistopic topic have have em- phasizedemphasized that that the thevalue value assumed assumed by bythe thegradie gradientnt of isotopic of isotopic content content varies varies according according to meteo-climaticto meteo-climatic characteristics, characteristics, too [33]. too In [33 Italy,]. In the Italy, average the average values of values vertical of isotopic vertical gra- iso- 18 18 dienttopic are gradient given areby givenGrad byO =Grad −0.24‰/100O = −0.24 m ‰/of elevation100 m of for elevation δ18O [12] for andδ O[ Grad12] H= and 2 2 −1.8‰/100Grad H = −m1.8 of‰/ elevation100 m of for elevation δ2H [12]. for Therefore,δ H[12]. in Therefore, this study, in thisthe study,average the elevation average 18 valueselevation of the values recharge of the area recharge (I-eleva areation) (I-elevation) were estimated were starting estimated from starting δ18O and from δ2δH valuesO and 2 detectedδ H values for detectedSabaudia for Station Sabaudia (17 m Station a.s.l.) (17[11]. m These a.s.l.) values [11]. These are computed values are by computed solving the by followingsolving the formulas: following formulas: 18δ O +(Grad18 O∙∆ )= δ18 O δ OSab + Grad O·∆q = δ O Spring;; (1) δ H + Grad H∙∆ )= δ H ; (2) 2 2 2 δ HSab + (Grad H·∆q) = δ Hspr; (2) where: 2 18 δ H and δ O are the δ H and δ O values detected for Sabaudia Station; 2 18 δ H and δ O are the δ H and δ O values detected for the spring. Considering that ∆q = difference between the elevation value of the mean recharge area and the altitude of the Sabaudia station, it results that: |δ O −δ O | = + ; (3) Grad O