Next Article in Journal
Assessing the Role of Extreme Mediterranean Events on Coastal River Outlet Dynamics
Next Article in Special Issue
Spatiotemporal Analysis of Groundwater Storage Changes, Controlling Factors, and Management Options over the Transboundary Indus Basin
Previous Article in Journal
Problems of Local Flooding in Functional Urban Areas in Poland
Previous Article in Special Issue
Hybrid Fuzzy Multi-Criteria Analysis for Selecting Discrete Preferable Groundwater Recharge Sites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Boundary Conditions Dynamics on Groundwater Budget in the Campania Region (Italy)

1
Department of Material, Environmental Sciences and Urban Planning, Polytechnic University of Marche, Via Brecce Bianche 12, 60131 Ancona, Italy
2
Department of Environmental, Biological and Pharmaceutical Sciences and Technologies, Second University of Naples, Via Vivaldi 43, 81100 Caserta, Italy
3
Jealott’s Hill International Research Centre, Syngenta Ltd., Warfield, Bracknell RG42 6EY, UK
*
Author to whom correspondence should be addressed.
Water 2022, 14(16), 2462; https://doi.org/10.3390/w14162462
Submission received: 14 June 2022 / Revised: 6 August 2022 / Accepted: 8 August 2022 / Published: 9 August 2022

Abstract

:
Groundwater budgets and fluxes are affected by human activities and climate change. Numerical models are cost-effective tools to investigate the different components of the hydrologic cycle. In this study, a groundwater flow model of the unconfined aquifers of the Campania region (Italy) has been developed and calibrated in Processing Modflow 11, resulting in an accurate assessment of groundwater fluxes and their trends over fifteen years (2000–2015). The model was implemented using a high-resolution grid to capture small hydrogeological features such as wells and rivers and informed by time variable datasets used as boundary conditions (i.e., river and sea levels, aquifer recharge, evapotranspiration, and discharge from adjacent systems). Good calibration and validation performances were achieved for piezometric heads (R2 = 0.958). A set of scenarios was developed using constant boundary conditions (i.e., constant sea-level BC, uniform extinction depth BC), and the outputs were compared, quantitively assessing differences in groundwater fluxes. Simulations pointed out that using time series to inform boundary conditions in the model does not always result in a significant change in the computed fluxes. Overall, non-uniform extinction depth was the most influential condition, while both rivers and sea level conditions barely affected groundwater budgets. In addition, results highlighted the need for an accurate estimation of spatiotemporal variations of both recharge and evapotranspiration, due to their strong seasonal variability and their massive contribution to the hydrogeological cycle. Finally, a marked increase of evapotranspiration fluxes controlled by interannual variability of precipitation and atmospheric temperatures has been quantified over the modelled period.

1. Introduction

The Mediterranean region, where exceptional levels of biodiversity occur [1], has been identified as one of the most prominent “Hot-Spots” in actual and future climate change (CC) projections [2]. Climatic disturbance in terms of increasing temperatures [3,4], global hydrological cycle disturbance [5], and sea-level rise [6,7,8] could greatly affect both groundwater systems and groundwater-dependent ecosystems. Indeed, groundwater inputs support or compensate for surface water inputs and play a vital role in maintaining healthy ecosystems [9]. Since CC has a significant effect on the hydrological cycle [10], is mandatory for water managers to plan appropriate measures to avoid negative potential impacts on water resources [11]. Most of the studies, published so far focused on the above-ground components of the hydrological cycle, analysing both historical and projected changes, while the scientific knowledge on the subsurface components, such as recharge, groundwater levels, aquifer fluxes, and groundwater quality have been only recently emerging in literature [12,13,14]. The main tools that the research community has developed to face these issues are observation networks, which allow to record the hydraulic parameters, and mathematical/physical models which simulate their evolution in time and space. Indeed, the implementation of numerical groundwater flow models is crucial to investigating groundwater-related issues by evaluating the trends of groundwater fluxes over a simulation time and thus quantifying the ongoing CC effects on the subsurface hydrological cycle [15]. Since its release in 1988, MODFLOW has become the worldwide standard for groundwater modelling because of its flexible modular structure, complete coverage of hydrogeological processes, and free availability [16]. MODFLOW can also be integrated with a geographic information system (GIS) to provide a visual environment for groundwater resources evaluation and management [17]. Application of groundwater flow models to large domains started a long time ago, in 1978, with the Regional Aquifer System Analysis (RASA) program of the U.S. Geological Survey [18], but only a few studies, so far, aimed to numerically evaluate changes in model outputs derived from constant or variable boundary conditions (BCs).
Coelho et al. [19] developed a study focusing on the influence of different BCs on numerical simulations of an unconfined aquifer in Brazil to define the most suitable BC, without considering scenarios with constant or variable BCs. Bastani and Harter [20] developed a three-dimensional transient flow and nitrate-transport model, comparing reference results to temporally upscaled models and considering annual-averaged flow and steady-state stresses across different scenarios. However, that study only focused on the concentrations of the nitrate discharging points without quantitively addressing changes in the groundwater budget. Moreover, employing proper spatial and temporal scales for both flow and transport simulations remains one of the most important challenges in numerical modelling [21,22,23]. One of the main tasks in transient groundwater modelling is the need for large amounts of hydrogeological data for both model implementation and calibration/validation; those constraints have often limited the use of transient numerical flow models to well-characterized field sites, while transient large-scale models are still rare [24]. Despite this, transient groundwater flow modelling has been recognized as one of the hydrological research fields extremely relevant for scientific progress [25]. In recent years, open datasets on hydrology and hydrogeology have been released by many research groups and international organizations, but their widespread use in numerical models still faces problems such as the lack of universally accepted guidelines or insufficient resolution [26]. To explore the new opportunities (and limitations) offered by open datasets to represent variable BCs in large-scale groundwater modelling, the unconfined aquifers of the Campania region (Italy) were selected as representative of potentially vulnerable to CC groundwater bodies. In fact, newly released open datasets on hydrology and hydrogeology can theoretically increase the model’s capability to quantify the variations of groundwater fluxes induced by spatially and temporally variable BCs. Over the last few years, several studies have focused on developing conceptual and numerical models of the Campania region (Italy) to tackle a series of groundwater-related issues [27,28,29,30], but without investigating the impact of BCs variability on the model outputs. This study aimed to fill this gap, assessing the influence of different implementation of BCs (e.g., constant versus variable sea-level BC, uniform versus variable extinction depth BC, etc.) over the computed groundwater flow budget from 2000 to 2015. This kind of numerical analysis would also help in reducing both human and computational efforts in similar settings, by inferring what datasets are worth being simplified (e.g., averaged over time and space) and which need to be treated as spatially variable and transient over time.

2. Materials and Methods

2.1. The Study Area

The Campania region, located in the southern part of the Italian Peninsula, covers an area of about km2. The region embraces a wide part of the Apennines chain, including composite geological and geomorphological units. The western slope of the Apennines is characterized by graben-like coastal depressions [31] and the plains alongside the coastal strip result from the aggradation of such structural depressions, whose bed is several thousand meters lower than the Apennine massifs [32]. The lithological and structural complexity of the region controls the hydrogeological fluxes.
The domain extends from the Tyrrhenian Sea to approximately 50 km inland, delimited in the South by the Cilento mountain range (Figure 1). The study area includes different alluvial plains (e.g., Garigliano, Volturno, Sarno, and Sele plains), characterized by debris and alluvial deposits. An important part is also occupied by limestone (e.g., M. Massico and Lattari Mts.) and volcanic units caused by the aggregation of volcanic material released during historical eruptions in the area (i.e., Vesuvio, Phlegraean Fields, and Roccamonfina volcanos are located within the model domain). In addition, the Volturno plain is characterized by a vast lowland area, subjected to land reclamation since the seventeenth century, which promoted the development of agriculture and farming along with strong coastal urbanization [33]. For toponyms used in the text refer to Figure S1.

2.2. Field Characterization and Boundary Conditions

Before the implementation of the groundwater flow model, all available datasets have been collected to define the conceptual model of the study area and the BCs (Table 1).
The regional SRTM Digital Terrain Model with a spatial resolution of 20 × 20 m cells was used to retrieve topographic information and represent morphology. The data files contain the elevation in a digital format, averaged on the model grid. Vegetation, buildings, and other cultural features were removed digitally, resulting in a proper terrain model.
The hydraulic conductivity of lithological units was obtained from permeability tests, by mining data from 91 borehole reports available on the online database of Instituto Superiore per la Protezione e Ricerca Ambientale [34]. All the data about pumping well rates, effective radius, screen length, and groundwater drawdown were retrieved from the same reports. The depth of the unconfined aquifer for each borehole was collected along with locations (latitude, longitude, and UTM coordinates). Other information on the aquifer’s depth has been retrieved from previous studies carried out in the Campania region, such as logs on the Volturno river plain [35] or aquifer information in Cuma [36,37], Bagnoli, Giuliano, and San Valentino Torio.
Measurments on sea levels were downloaded from the open dataset Servizio Mareografico Nazionale (SMN). SMN exploits a network of 36 measuring stations uniformly distributed over the national territory and located mainly inside the port facilities. The hydrometric level is monitored since 2010 using a new microwave (radar) level meter with millimetric precision, coupled with a second level meter with a floater. Each level is referred to the absolute altimetric network of Istituto Geografico Militare (IGM). To gain a better understanding of the sea level behaviour across the entire region, Napoli, Salerno, Gaeta, and Palinuro stations were considered. However, since SMN open dataset does not cover the entire period of simulation, the datasets of the Copernicus Marine Environment Monitoring Service (CMEMS) were also exploited for sea-level trends at the same locations. CMEMS altimetric satellite data were elaborated together with the SMN ones allowing to cover the entire period of study.
The Moderate Resolution Imaging Spectroradiometer (MODIS) products have been deployed to obtain the trends of evapotranspiration in the study area for the period 2000–2015. Data were downloaded using the AppEEARS [38] interface and the MYD16A3GF model [39] to obtain spatially interpolated values. The algorithm used in MODIS is based on the Penman-Monteith equation, which includes daily meteorological data along with remotely sensed variables, such as vegetation dynamics, albedo, and land cover [40]. Land use and soil texture maps (Table 1) were overlapped via GIS to derive a two-dimensional matrix representing the evapotranspiration extinction depth [41].
Time series of groundwater recharge has been extracted from the outputs of PCR-GLOBWB v2.0 (Utrecht University, Utrecht, The Netherland) [42], a grid-based global hydrology and water resources model. Global hydrological models (GHMs) are tools able to quantify the global terrestrial water cycle, simulating the distributed hydrological response to weather and climate variations at high resolution.
Hydrometric levels were downloaded from the Centro Funzionale Multirischi (CFM), a Civil Protection Agency of the Campania region, to represent the variation of the river stage over the simulation period. CFM provides a list of all hydrometric stations of the region recording data with a daily frequency. The closest stations to the model’s boundaries were considered representative of the baseflow entering the domain from the fractured aquifers. Hydraulic head fluctuation in the fractured aquifers was modified according to the average baseflow variations over the simulation time.

2.3. Numerical Modelling

To estimate the groundwater flow budget at a regional scale, assessing the influence of variable BCs, all the data were elaborated in a GIS environment and processed via Processing Modflow 11. The flow simulation was carried out using the MODFLOW-2005 code [43]. The equation that stands at the base of MODFLOW-2005 is the partial-differential equation that describes the three-dimensional movement of groundwater of constant density through a porous matrix:
x K xx h x + y K yy h y + z K zz h z + W = S s h t
where Kxx, Kyy and Kzz are values of hydraulic conductivity along x, y and z coordinate axes (LT−1); h is the potentiometric head (L); W is a volumetric flux per unit volume representing sources and/or sinks of water (T−1); Ss is the specific storage of the porous material (L−1); and t is time (T). The standard version of MODFLOW-2005 treating the fractured units as equivalent porous media and not complex ones such as the Conduit Flow Package [44] was used given the lack of information on the location and dimensions of such zones. Thus, for the principle of parsimony in numerical modelling, a simpler representation was here adopted [45].
The model domain consisted of 656 rows, 335 columns, and 3 layers with variable width and thickness. The top of the grid represents the local ground surface, ranging from −3.49 m to 1642 m above sea level (a.s.l.). The bottom of the grid represents the bottom of the unconfined aquifers and ranges from −83 m below sea level (b.s.l.) to 144 m a.s.l.
Initially, a steady-state simulation has been run and the model has been calibrated via PEST [46] against an observed piezometric map of the Campania region for the year 2000 [47]. Given that the piezometric map was available only as a vector map, it was sampled along the contour lines to create 500 observation points that were compared versus calculated ones with an estimated error of ±1.0 m. Tyrrhenian Sea levels have been initially set to 0 m as constant head cells. Values of hydraulic conductivity, corresponding to boreholes in the same cluster, have been averaged to get an average value for each geological formation in the study area. In this way, the model domain was subdivided into polygons whose extension corresponds to the extension of the different geological formations reported in the geological map. Figure 2 highlights the heterogeneity of the Campania region, where the horizontal hydraulic conductivity values range from 1 × 10−6 m/s in the coastal part of the Volturno plain, where clay layers are present, to 6.6 × 10−4 m/s in the Vesuvio area, where highly fractured volcanic rocks are present. The value of the vertical hydraulic conductivity was set equal to the value of the horizontal hydraulic conductivity assigned to that cell in fractured units, while it was set as a tenth of the horizontal hydraulic conductivity in sedimentary units, to account for internal layering. It must be stressed, however, that local scale heterogeneities such as intensively faulted zones, small lenses, etc. have not been explicitly incorporated in this large-scale model.
For the steady-state simulation, evapotranspiration rates for the year 2000 were averaged from the MODIS dataset and interpolated on the numerical grid; while the recharge rates were averaged from the PCR-GLOBWB v2.0 for the same period and interpolated on the numerical grid. The Processing Modflow polyline input method in the River Package was employed to simulate, from North to South, Garigliano, Savone, Agnena, Volturno, Regi Lagni, Sarno, Picentino, Tusciano, Sele, and Calore rivers, drawing their actual pathways and defining river parameters. The riverbed conductance per unit length Crivl (LT−1), the river stage hriv (L), and the elevation of the riverbed bottom Briv (L) have been assigned to the most upstream and downstream vertices of each river polyline and then applied to the model cells employing a linear interpolation. River stages were assumed close to the assigned starting head value and kept constant for the entire steady-state simulation. The flow rate Qriv (L3T−1) between each river reach and the model cell is calculated by the following equations:
Q riv =   C riv · h riv h ,         i f   h > B riv   C riv · ( h riv B riv ) ,   i f   h < B riv  
where h (L) is the calculated hydraulic head in the cell containing the reach.
The wells analysed from the online ISPRA database were entered into the model through the Well Package. Industrial wells were set with continuous pumping rates, while the agricultural wells were set with discontinuous pumping rates, active only during the spring and summer seasons. Since the specified wells penetrate more than one model layer, the total pumping rate has not been distributed over the penetrated layers, but it has been assigned to the lowest cell of the well, belonging to the third model layer. The Drain Package was used to simulate the reclamation channels in the Marina di Licola area and also to simulate the effect of the underground metro lines in Naples, the mean height of each feature a.s.l. was derived from the topographical online database, while the cell conductance Cd (L2T−1) was set to a very high value (0.01) to ensure optimal groundwater withdrawal from the aquifer. The location of the depth of the different metro lines was sourced from Russo et al. [48]. The discharge rate Qd (L3T−1) from the model cell to the drain is calculated by:
Q d = C d · h h d ,     i f   h > h d 0 ,     i f   h h d  
The horizontal hydraulic conductivity was estimated by inverse modelling using the numerical code PEST, assigning a parameter number to it. PEST searches optimum parameter values for which the sum of squared deviations between model-calculated and observed values from the piezometric map of the region is reduced to a minimum. Such comparison has been performed using a scatter diagram of the simulated and observed data (Figure 3).
To investigate possible changes in the groundwater system associated with temporal variations of evapotranspiration and recharge patterns, a transient simulation was carried out. Sixty-four stress periods of seasonal length (3 months), subdivided into three time steps of 1 month each, were used. As recommended by Harbaugh [43], the starting hydraulic heads for the transient simulation have been set equal to the head values calculated at the end of the steady-state simulation. Table 2 shows the specific storage (Ss) and specific yield (Sy) values, depending on the geologic material, together with the riverbed conductance (Criv) gained after the calibration procedure.
To consider the variability of the evapotranspiration over time, the MODIS data have been spatialized using the kriging technique and inserted into the model with seasonal frequency. The values of recharge rate from the global numerical hydrologic model published by de Graaf et al. [49] have been processed in python to obtain spatially interpolated values for the period 2000–2015, which have been assigned to the model with seasonal frequency. The calibrated model was implemented with the before-mentioned transient datasets accounting for the seasonal changes in hydrometric levels, Tyrrhenian Sea levels, and fractured aquifer seasonal head fluctuations as variable BCs. The Evapotranspiration Package was also integrated with spatially variable extinction depths dependent on the land use and soil texture maps.
The validation of the transient state model was performed by exploiting a piezometric map of the Campania region for the year 2013, which was derived from a secondary elaboration of the global groundwater depth dataset [50] to match modelled and measured hydraulic heads. Besides, selected monitoring wells (Figure 1) with piezometric trend data were compared versus model results.
Both the model calibration and validation were assessed by calculating and comparing the value of the coefficient of determination (R2) and the mean absolute error (MAE).
The Processing Modflow ZoneBudget routine [43] can produce water budgets for user-defined zones within the model domain. Zones are defined by assigning zone numbers to model cells, allowing for to calculation of the volumetric water budget in a different part of the domain, accounting for evapotranspiration, fractured aquifers, rivers, and the Tyrrhenian Sea component separately. Each trend calculated by the model, considering all the transient input datasets as variable BCs (transient model), was compared with the respective trend calculated by an upscaled model considering a uniform extinction depth, a constant head at the fractured aquifers, and constant sea and river levels (basic model). Moreover, to investigate which datasets are worth being averaged, a set of scenarios have been developed applying to the transient model each of the averaged BCs individually:
  • constant head at the fractured aquifers BC (scenario S1);
  • uniform extinction depth BC (scenario S2);
  • constant Tyrrhenian Sea level BC (scenario S3);
  • constant river levels BC (scenario S4).
In addition, to numerically quantify the impact of the averaged conditions, the 15-year average flux of each groundwater budget component has been estimated and differences between basic, transient, S1, S2, S3, and S4 simulations were graphically compared.
Finally, since the pivotal contribution to the groundwater budget of both recharge and evapotranspiration, an additional recharge and evapotranspiration rates dataset has been investigated. Alternatively, to the PCR-GLOBWB v.2.0 published by Sutanudjaja et al. and the MODIS products, the BIGBANG (Bilancio Idrogeologico GIS BAsed a scala Nazionale su Griglia regolare) [51] distributed model was employed, and results compared.

3. Results and Discussion

3.1. Flow Model Performance

Given the extremely challenging regional scale of the domain (3672 km2), reasonably good calibration and validation were achieved for hydraulic heads, resulting in an R2 of 0.958 (Figure 3). MAE values of 3.56 m for calibration and 2.87 m for validation, were judged acceptable in the wide range of the simulated piezometric levels (>100 m). Observations of groundwater levels from selected monitoring wells were used to verify the ability of the model to describe water table fluctuations in the unconfined aquifers (Figure 3). Unfortunately, a widespread monitoring network of wells is not available in the region, and this is a limitation for the model performance. Indeed, monitoring groundwater levels at regional scale may enable a better assessment of groundwater depletion, as well as inferring changes in flow directions over time, as recently pointed out by Brakkee et al. [52]. Moreover, piezometric data are missing in some karst and fractured areas of the domain (Massico mountain, Vesuvio volcano, and Monti Lattari) consequently model reliability is lower on those areas than in alluvial zones. By means of the Processing Modflow ZoneBudget routine [43] a defined zone number was assigned to those areas and their importance in relation to the overall model domain was assessed in terms of water budget component percentages. As expected, groundwater recharge played a very important role in karst environments, representing 12.7% of the overall domain, while evapotranspiration was just 1.5% with respect to the total, given that in mountain areas the groundwater depth from the surface is higher than the extinction depth. The storage capacity turns out to be consistent, whereas all the other budget components resulted to be negligible. The general high permeability grade due to fracturing leads to a high groundwater recharge rate [53] that varies between 48% and 78% of the mean annual effective precipitation [54,55], depending on the relative abundance of limestone and dolomites lithologies.
The validated transient model allowed analysing and quantifying relevant stresses in different portions of the domain at different times. Figure 4 shows an example of both the three-dimensional and planar view of the calculated hydraulic heads map at the end of the simulation. It pointed out regions of the domain characterized by high hydraulic gradient, especially the northern part in correspondence of the Roccamonfina Volcano, which acted as a groundwater mounding area, as well as the Vesuvio, Phlegraean Fields, and Monti Lattari. Hydraulic heads resulted to be very low in the Volturno river plain, following the topographic gradient. Near the coast, the water level was very shallow, due to reclamation channels located in the Marina di Licola area (lower Volturno plain).
In addition, the model produced accurate estimates of the groundwater budget with a mean error of less than 0.3% among in and out components. This allowed focusing on the single water budget components, whose fluxes are difficult to estimate without a physically based model even though they are crucial for better management of the hydrogeological basins.

3.2. Water Fluxes Trends and Comparative Analysis

The fluxes trends have been estimated accounting for the inflow from the fractured aquifers, the evapotranspiration, the river-groundwater interactions, and the Tyrrhenian Sea component, separately. To highlight the influence of the time variable datasets, the basic model calculated trends (red line), the transient model calculated trends (blue line), and the scenarios (S1, S2, etc.) calculated trends (dashed green line) associated with the budget component have been plotted together.
Figure 5a displays the inflow component from the fractured aquifers. The transient model calculated trend shows a quite higher flux with a wider seasonal variation throughout the simulation time, with respect to the basic and S1 model.
The evapotranspiration component represents the most important outflow contribution. Figure 5b clearly shows its seasonal variability with increasing fluxes over the warm seasons, reaching values up to 20 m3/s. Moreover, since the extinction depth affects the evapotranspiration fluxes, the S2 trend perfectly reflects the basic model, underlying the importance of considering a non-uniform extinction depth over the domain to account for different land cover and different soil types that directly drive the extinction depth [41].
The seawater intrusion component (Figure 5c) is characterized by a very low flux and can be considered negligible with respect to submarine groundwater discharge. However, importing a transient sea-level dataset leads the seawater intrusion flux of the transient model to display positive peaks due to the seasonal oscillation of the Tyrrhenian Sea. This could be further investigated at the local scale by performing a variable-density flow and transport simulation. There are no evident changes between the basic, transient, and S3 scenarios in submarine groundwater discharge, which shows a mean flux over the simulation time between 5 and 6 m3/s (Figure 5d).
Both river inflow and outflow (Figure 5e,f), which represent the losing and gaining flow rate from and to the river, show no significant changes between basic, transient, and S4 simulations. Therefore, considering the effort to import all required data in the models, a constant river level dataset averaged over the simulation time should be considered for large-scale simulations.
To numerically quantify the impact of time variable input datasets on groundwater budget, the 15-year average flux of each component has been estimated and differences between basic, transient, S1, S2, S3, and S4 simulations were represented in Figure 6. A difference of ~1.67 m3/s was found in the mean evapotranspiration between the basic and the transient simulations. Almost the same difference (~1.70 m3/s) was found between S2 and the transient model. This confirmed that extinction depth is the most influencing parameter, requiring a precise spatiotemporal estimation to obtain accurate predictions of the groundwater budget.
The difference in the mean inflow from the fractured aquifers component (~0.65 m3/s) between the basic and transient models did not seem to be entirely due to the associated transient inflow parameters, but also accounted for the variable river level dataset.
Even if the difference in the averaged seawater intrusion component was relatively high (almost doubled from the basic to transient model), it should be considered negligible in absolute terms (0.02 m3/s). The mean submarine groundwater discharge showed no appreciable changes between the scenarios, both in relative and absolute terms. Effectively, that component seems to be more affected by the non-uniform extinction depth dataset (0.16 m3/s difference) than by the transient Tyrrhenian Sea levels (0.06 m3/s difference), suggesting that sea level could be kept constant over the simulation time.
The same analysis can be conducted considering the mean groundwater flux towards the river (river outflow). The difference of 0.98 m3/s between the basic and the transient simulations was mostly due to the non-uniform extinction depth rather than to the variable river level dataset: the value rose to 1.02 m3/s between the S2 and the transient model, while it dropped to 0.05 m3/s between the S4 and the transient model. Considering also that the mean river flux towards the unconfined aquifers (river inflow) slightly differs among the various scenarios, accounting for a variable river level dataset does not seem to be crucial.
The comparative analysis pointed out the fundamental role of the extinction depth in shifting the allocation of groundwater fluxes among different compartments, altering the contribution of each component with respect to the total budget.
Moreover, since the evaluation of each groundwater contribution underlines the pivotal role of both the recharge (with a mean flux of 18.1 m3/s) and evapotranspiration components, the additional BIGBANG dataset has been investigated. The derived recharge and evapotranspiration trends have been calculated and compared with the transient model results (Figure 7).
The recharge flow rates strongly differ among the two models, while the evapotranspiration fluxes were more similar. The linear temporal trends show an increase in the recharge flux over time for both the PCR-GLOBWB v2.0 and the BIGBANG products, reflecting the positive trends of the hydraulic heads observed in the monitoring wells P1, P2, P3, and P4 (Figure 3). The process of groundwater recharge is driven by a decadal variability of the North Atlantic Oscillation and its robust correlation with the rainfall patterns [56]. Both the calculated evapotranspiration fluxes, controlled by interannual variability of precipitation and air temperature, showed increasing trends with a small positive shift (~0.10 m3/s per year), highlighting that CC impact on groundwater resources is confirmed by the two independent databased here employed. Furthermore, because of the overestimation of the recharge flux due to topographic and geological heterogeneities of the Campania region, running the transient model using the BIGBANG products slightly decrease the performance parameters, with an R2 = 0.948 and an MAE = 3.42 m, producing flooded areas during the peak recharge seasons that were not physically based.
Finally, it must be stressed that all other scenarios showed very little variation of the model performance indicators with an R2 and an MAE that never exceeded 0.947 and 2.91 m, respectively. This is due to the low sensitivity of piezometric heads in capturing the changes induced by evapotranspiration at the basin scale [57], while to solve this complex issue, monitoring wells with continuous heads dataloggers should be employed coupled with baseflow measurements [58].

4. Conclusions

Open online datasets, along with historical observations and literature review were used to untangle the complexity of regional hydrogeological settings, supporting the implementation of a robust numerical model to describe groundwater fluxes for the unconfined aquifers of the Campania region. The quantification of the total flux of each component in the budget highlighted the pivotal role of recharge and evapotranspiration rates, with a mean flux of 18.1 and 12.3 m3/s, respectively. Their variable patterns throughout the simulation period affected the whole hydrogeological budget, pointing out the need for their accurate estimation over time, to account for their high seasonal variability and climatic trends. Recharge can reach a value of above 60 m3/s during wet periods, while evapotranspiration variability ranges from 5 to maximum values up to 20 m3/s over the warm seasons. Time variable input datasets had very different effects on the components of the groundwater budget, controlling water fluxes between different compartments. Simulations pointed out that both fractured aquifers head fluctuations and non-uniform extinction depth datasets used as BCs should not be averaged. Considering transient data of extinction depth had a massive effect on the budget results, not only influencing evapotranspiration rates directly but also affecting all the other components of the groundwater budget, indirectly. Conversely, the inclusion of other datasets, namely variable river and sea levels, barely impacted the groundwater budget, not making it worth the extra effort. Therefore, the inclusion of transient BCs in a numerical model should be based on expert judgment, balancing accuracy and complexity and considering the human and computational efforts required. Finally, the results have underlined that continuous groundwater level and baseflow estimates are crucial to minimize model uncertainties and increase the reliability of groundwater budgets at the basin scale. Therefore, such measurements need to be included in an ideal groundwater monitoring network at a regional scale.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w14162462/s1, Figure S1: The domain is limited by the Tyrrhenian Sea (light blue) and by the Apennines and the Cilento mountain range (light green). The area includes the Garigliano, Volturno, Sarno, and Sele plains (yellow), the volcanic units Vesuvio, Phlegraean Fields, and Roccamonfina (white) and the limestone units Mount Massico, Mount Maggiore, and Mounts Lattari (white). More in detail: the Garigliano plain hosts the Garigliano river, the Volturno plain hosts the Savone, Agnena, Volturno, and Regi Lagni rivers, the Sarno plain hosts the Sarno rivers, and the Sele plain hosts the Picentino, Tusciano, Sele and Calore rivers.

Author Contributions

Conceptualization, M.M. and N.C.; methodology, N.C., M.M. and M.G.; software, M.G. and N.C.; validation, M.M. and N.C.; formal analysis, M.G.; investigation, M.G.; data curation, M.G., F.R. and G.B.; writing—original draft preparation, M.G.; writing—review and editing, N.C., M.M. and F.R.; visualization, M.G. and G.B.; supervision, N.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data will be made available by the authors upon request.

Acknowledgments

Simone Tavoni is acknowledged for helping to gather the raw data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Blondel, J.; Aronson, J. Biology and Wildlife of the Mediterranean Region, 1st ed.; Oxford University Press: Oxford, UK, 1999; p. 328. [Google Scholar]
  2. Lionello, P.; Scarascia, L. The relation between climate change in the Mediterranean region and global warming. Reg. Environ. Chang. 2018, 18, 1481–1493. [Google Scholar] [CrossRef]
  3. Cramer, W.; Guiot, J.; Fader, M.; Garrabou, J.; Gattuso, J.P.; Iglesias, A.; Xoplaki, E. Climate change and interconnected risks to sustainable development in the Mediterranean. Nat. Clim. Chang. 2018, 8, 972–980. [Google Scholar] [CrossRef] [Green Version]
  4. Hallegatte, S. Strategies to adapt to an uncertain climate change. Global. Environ. Chang. 2009, 19, 240–247. [Google Scholar] [CrossRef]
  5. IPCC. Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; Masson-Delmotte, V.P., Zhai, A., Pirani, S.L., Connors, C., Péan, S., Berger, N., Caud, Y., Chen, L., Goldfarb, M.I., Gomis, M., et al., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2021; in press. [Google Scholar] [CrossRef]
  6. FitzGerald, D.M.; Fenster, M.S.; Argow, B.A.; Buynevich, I.V. Coastal impacts due to sea-level rise. Earth Planet. Sci. 2008, 36, 601–647. [Google Scholar] [CrossRef] [Green Version]
  7. Carrasco, A.R.; Ferreira, Ó.; Roelvink, D. Coastal lagoons and rising sea level: A review. Earth Sci. 2016, 154, 356–368. [Google Scholar] [CrossRef] [Green Version]
  8. Benjamin, J.; Rovere, A.; Fontana, A.; Furlani, S.; Vacchi, M.; Inglis, R.H.; Gehrels, R. Late Quaternary sea-level changes and early human societies in the central and eastern Mediterranean Basin: An interdisciplinary review. Quat. Int. 2017, 449, 29–57. [Google Scholar] [CrossRef] [Green Version]
  9. Erostate, M.; Huneau, F.; Garel, E.; Ghiotti, S.; Vystavna, Y.; Garrido, M.; Pasqualini, V. Groundwater dependent ecosystems in coastal Mediterranean regions: Characterization, challenges and management for their protection. Water Res. 2020, 172, 115461. [Google Scholar] [CrossRef]
  10. Famiglietti, J.S. The global groundwater crisis. Nat. Clim. Chang. 2014, 4, 945–948. [Google Scholar] [CrossRef]
  11. Erturk, A.; Ekdal, A.; Gurel, M.; Karakaya, N.; Guzel, C.; Gonenc, E. Evaluating the impact of climate change on groundwater resources in a small Mediterranean watershed. Sci. Tot. Environ. 2014, 499, 437–447. [Google Scholar] [CrossRef] [PubMed]
  12. Butscher, C.; Huggenberger, P. Modeling the temporal variability of karst groundwater vulnerability, with implications for climate change. Environ. Sci. Technol. 2009, 43, 1665–1669. [Google Scholar] [CrossRef] [Green Version]
  13. Oude Essink, G.H.P.; Van Baaren, E.S.; De Louw, P.G. Effects of climate change on coastal groundwater systems: A modeling study in the Netherlands. Water Resour. Res. 2010, 46. [Google Scholar] [CrossRef]
  14. Mastrocicco, M.; Gervasio, M.P.; Busico, G.; Colombani, N. Natural and anthropogenic factors driving groundwater resources salinization for agricultural use in the Campania plains (Southern Italy). Sci. Total Environ. 2021, 758, 144033. [Google Scholar] [CrossRef]
  15. Allen, D.M.; Mackie, D.C.; Wei, M. Groundwater and climate change: A sensitivity analysis for the Grand Forks aquifer, southern British Columbia, Canada. Hydrogeol. J. 2004, 12, 270–290. [Google Scholar] [CrossRef]
  16. Zhou, Y.; Li, W. A review of regional groundwater flow modeling. Geosci. Front. 2011, 2, 205–214. [Google Scholar] [CrossRef] [Green Version]
  17. Aghlmand, R.; Abbasi, A. Application of MODFLOW with boundary conditions analyses based on limited available observations: A case study of Birjand plain in East Iran. Water 2019, 11, 1904. [Google Scholar] [CrossRef] [Green Version]
  18. Sun, R.J.; Johnston, R.H. Regional Aquifer-System Analysis Program of the US Geological Survey, 1978–1992; US Government Printing Office: Denver, CO, USA, 1994; p. 126. [Google Scholar] [CrossRef]
  19. Coelho, C.D.; Faria, A.C.S.; Marques, E.A.G. Comparative analysis of different boundary conditions and their influence on numerical hydrogeological modeling of Palmital watershed, southeast Brazil. J. Hydrol. 2017, 12, 210–219. [Google Scholar] [CrossRef]
  20. Bastani, M.; Harter, T. Effects of upscaling temporal resolution of groundwater flow and transport boundary conditions on the performance of nitrate-transport models at the regional management scale. Hydrogeol. J. 2020, 28, 1299–1322. [Google Scholar] [CrossRef] [Green Version]
  21. Cushman, J.H.; Bennethum, L.S.; Hu, B.X. A primer on upscaling tools for porous media. Adv. Water Resour. 2002, 25, 1043–1067. [Google Scholar] [CrossRef]
  22. Dentz, M.; Le Borgne, T.; Englert, A.; Bijeljic, B. Mixing, spreading and reaction in heterogeneous media: A brief review. J. Contam. Hydrol. 2011, 120, 1–17. [Google Scholar] [CrossRef] [PubMed]
  23. Scheibe, T.D.; Murphy, E.M.; Chen, X.; Rice, A.K.; Carroll, K.C.; Palmer, B.J.; Wood, B.D. An analysis platform for multiscale hydrogeologic modeling with emphasis on hybrid multiscale methods. Groundwater 2015, 53, 38–56. [Google Scholar] [CrossRef] [PubMed]
  24. Hellwig, J.; de Graaf, I.E.M.; Weiler, M.; Stahl, K. Large-scale assessment of delayed groundwater responses to drought. Water Resour. Res. 2020, 56, e2019WR025441. [Google Scholar] [CrossRef] [Green Version]
  25. Fatichi, S.; Vivoni, E.R.; Ogden, F.L.; Ivanov, V.Y.; Mirus, B.; Gochis, D.; Downer, C.W.; Camporese, M.; Davison, J.H.; Ebel, B.A.; et al. An overview of current applications, challenges, and future trends in distributed process-based models in hydrology. J. Hydrol. 2016, 537, 45–60. [Google Scholar] [CrossRef] [Green Version]
  26. Rau, G.C.; Cuthbert, M.O.; Post, V.E.A.; Daniel Schweizer, D.; Acworth, R.I.; Andersen, M.S.; Blum, P.; Carrara, E.; Rasmussen, T.C.; Ge, S. Future-proofing hydrogeology by revising groundwater monitoring practice. Hydrogeol. J. 2020, 28, 2963–2969. [Google Scholar] [CrossRef]
  27. De Vita, P.; Fusco, F.; Tufano, R.; Cusano, D. Seasonal and event-based hydrological and slope stability modeling of pyroclastic fall deposits covering slopes in Campania (Southern Italy). Water 2018, 10, 1140. [Google Scholar] [CrossRef] [Green Version]
  28. Corniello, A.; Ducci, D.; Ruggieri, G.; Iorio, M. Complex groundwater flow circulation in a carbonate aquifer: Mount Massico (Campania Region, Southern Italy). Synergistic hydrogeological understanding. J. Geochem. Explor. 2018, 190, 253–264. [Google Scholar] [CrossRef]
  29. Mastrocicco, M.; Busico, G.; Colombani, N.; Vigliotti, M.; Ruberti, D. Modelling actual and future seawater intrusion in the Variconi coastal wetland (Italy) due to climate and landscape changes. Water 2019, 11, 1502. [Google Scholar] [CrossRef] [Green Version]
  30. Lancia, M.; Petitta, M.; Zheng, C.; Saroli, M. Hydrogeological insights and modelling for sustainable use of a stressed carbonate aquifer in the Mediterranean area: From passive withdrawals to active management. J. Hydrol. 2020, 32, 100749. [Google Scholar] [CrossRef]
  31. Brancaccio, L.; Cinque, A.; Romano, P.; Rosskopf, C.; Russo, F.; Santangelo, N.; Santo, A. Geomorphology and neotectonic evolution of a sector of the Tyrrhenian flank of the southern Apennines (Region of Naples, Italy). Z. Geomorphol. 1991, 82, 47–58. [Google Scholar]
  32. Pappone, G.; Alberico, I.; Amato, V.; Aucelli, P.P.C.; Di Paola, G. Recent evolution and present-day conditions of the Campanian Coastal plains (South Italy): The case history of Sele River Coastal plain. WIT Trans. Ecol. Envir. 2011, 149, 15–27. [Google Scholar] [CrossRef]
  33. Ruberti, D.; Vigliotti, M. Land use and landscape pattern changes driven by land reclamation in a coastal area: The case of Volturno delta plain, Campania region, southern Italy. Environ. Earth Sci. 2017, 76, 694. [Google Scholar] [CrossRef]
  34. ISPRA. National Geological Database. 2021. Available online: http://sgi2.isprambiente.it/viewersgi2/ (accessed on 10 October 2021).
  35. Corniello, A.; Ducci, D. Hydrogeochemical characterization of the main aquifer of the “Litorale Domizio-Agro Aversano NIPS” (Campania—Southern Italy). J. Geochem. Explor. 2014, 137, 1–10. [Google Scholar] [CrossRef]
  36. Allocca, V.; Coda, S.; De Vita, P.; Di Rienzo, B.; Ferrara, L.; Giarra, A.; Mangoni, O.; Stellato, L.; Trifuoggi, M.; Arienzo, M. Hydrogeological and hydrogeochemical study of a volcanic-sedimentary coastal aquifer in the archaeological site of Cumae (Phlegraean Fields, southern Italy). J. Geochem. Explor. 2018, 185, 105–115. [Google Scholar] [CrossRef]
  37. Stellato, L.; Coda, S.; Arienzo, M.; De Vita, P.; Di Rienzo, B.; D’Onofrio, A.; Ferrara, L.; Marzaioli, F.; Trifuoggi, M.; Allocca, V. Natural and anthropogenic groundwater contamination in a coastal volcanic-sedimentary aquifer: The case of the archaeological site of Cumae (Phlegraean Fields, Southern Italy). Water 2020, 12, 3463. [Google Scholar] [CrossRef]
  38. AppEEARS Team. Application for Extracting and Exploring Analysis Ready Samples (AppEEARS); Ver. 2.54.1. NASA EOSDIS Land Processes Distributed Active Archive Center (LP DAAC); USGS/Earth Resources Observation and Science (EROS) Center: Sioux Falls, SD, USA, 2020. Available online: https://lpdaacsvc.cr.usgs.gov/appeears (accessed on 10 October 2021).
  39. Running, S.; Mu, Q.; Zhao, M.; Moreno, A. MOD16A2GF MODIS/Terra Net Evapotranspiration Gap-Filled 8-Day L4 Global 500m SIN Grid V006. NASA EOSDIS Land Processes DAAC. 2019. Available online: https://doi.org/10.5067/MODIS/MOD16A2GF.006 (accessed on 10 October 2021).
  40. Mu, Q.; Zhao, M.; Running, S.W. Improvements to a MODIS global terrestrial evapotranspiration algorithm. Rem. Sens Environ. 2011, 115, 1781–1800. [Google Scholar] [CrossRef]
  41. Shah, N.; Nachabe, M.; Ross, M. Extinction depth and evapotranspiration from ground water under selected land covers. Groundwater 2007, 45, 329–338. [Google Scholar] [CrossRef]
  42. Sutanudjaja, E.H.; van Beek, R.; Wanders, N.; Wada, Y.; Bosmans, J.H.C.; Drost, N.; van der Ent, R.J.; de Graaf, I.E.M.; Hoch, J.M.; de Jong, K.; et al. PCR-GLOBWB 2: A 5 arcmin global hydrological and water resources model. Geosci. Model Dev. 2018, 11, 2429–2453. [Google Scholar] [CrossRef] [Green Version]
  43. Harbaugh, A.W. MODFLOW-2005, the US Geological Survey Modular Ground-Water Model: The Ground-Water Flow Process; U.S. Geological Survey Techniques and Methods 6-A16; US Department of the Interior; US Geological Survey: Reston, VA, USA, 2005; p. 253. [Google Scholar]
  44. Reimann, T.; Hill, M.E. MODFLOW-CFP: A New Conduit Flow Process for MODFLOW-2005. Groundwater 2009, 47, 321–325. [Google Scholar] [CrossRef]
  45. Hill, M.C. The practical use of simplicity in developing ground water models. Groundwater 2006, 44, 775–781. [Google Scholar] [CrossRef]
  46. Doherty, J. PEST-Model-Independent Parameter Estimation. Version 12. Watermark Computing. Australia. 2010. Available online: http://www.pesthomepage.org/ (accessed on 10 October 2021).
  47. De Vita, P.; Allocca, V.; Celico, F.; Fabbroncino, S.; Cesaria, M.; Monacelli, G.; Musilli, I.; Piscopo, V.; Scalise, A.R.; Summa, G.; et al. Hydrogeology of continental southern Italy. J. Maps 2018, 14, 230–241. [Google Scholar] [CrossRef] [Green Version]
  48. Russo, G.; Viggiani, C.; Viggiani, G. Geotechnical design and construction issues for Lines 1 and 6 of the Naples Underground/Geotechnische Planung und Ausführung der U-Bahnlinien 1 und 6 in Neapel. Geomech. Tunnelbau 2012, 5, 300–311. [Google Scholar] [CrossRef]
  49. de Graaf, I.E.; van Beek, R.L.; Gleeson, T.; Moosdorf, N.; Schmitz, O.; Sutanudjaja, E.H.; Bierkens, M.F. A global-scale two-layer transient groundwater model: Development and application to groundwater depletion. Adv. Water Resour. 2017, 102, 53–67. [Google Scholar] [CrossRef]
  50. Fan, Y.; Li, H.; Miguez-Macho, G. Global patterns of groundwater table depth. Science 2013, 339, 940–943. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Braca, G.; Bussettini, M.; Ducci, D.; Lastoria, B.; Mariani, S. Evaluation of national and regional groundwater resources under climate change scenarios using a GIS-based water budget procedure. Rend. Fis. Acc. Lincei 2019, 30, 109–123. [Google Scholar] [CrossRef]
  52. Brakkee, E.; Van Huijgevoort, M.H.; Bartholomeus, R.P. Improved understanding of regional groundwater drought development through time series modelling: The 2018–2019 drought in the Netherlands. Hydrol. Earth Syst. Sci. 2022, 26, 551–569. [Google Scholar] [CrossRef]
  53. Ruggieri, G.; Allocca, V.; Borfecchia, F.; Cusano, D.; Marsiglia, P.; De Vita, P. Testing Evapotranspiration Estimates Based on MODIS Satellite Data in the Assessment of the Groundwater Recharge of Karst Aquifers in Southern Italy. Water 2021, 13, 118. [Google Scholar] [CrossRef]
  54. Allocca, V.; Manna, F.; De Vita, P. Estimating annual groundwater recharge coefficient for karst aquifers of the southern Apennines (Italy). Hydrol. Earth Syst. Sci. 2014, 18, 803–817. [Google Scholar] [CrossRef] [Green Version]
  55. Allocca, V.; De Vita, P.; Manna, F.; Nimmo, J.R. Groundwater recharge assessment at local and episodic scale in a soil mantled perched karst aquifer in southern Italy. J. Hydrol. 2015, 529, 843–853. [Google Scholar] [CrossRef]
  56. De Vita, P.; Allocca, V.; Manna, F.; Fabbrocino, S. Coupled decadal variability of the North Atlantic Oscillation, regional rainfall and karst spring discharges in the Campania region (southern Italy). Hydrol. Earth Syst. Sci. 2012, 16, 1389–1399. [Google Scholar] [CrossRef] [Green Version]
  57. Colombani, N.; Gaiolini, M.; Busico, G.; Postacchini, M. Quantifying the impact of evapotranspiration at the aquifer scale via groundwater modelling and MODIS data. Water 2021, 13, 950. [Google Scholar] [CrossRef]
  58. Lehr, C.; Lischeid, G. Efficient screening of groundwater head monitoring data for anthropogenic effects and measurement errors. Hydrol. Earth Syst. Sci. 2020, 24, 501–513. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Model location with initial heads and boundary conditions (BCs).
Figure 1. Model location with initial heads and boundary conditions (BCs).
Water 14 02462 g001
Figure 2. Distribution of the lithological units and their horizontal hydraulic conductivity values in layer 1 of the model domain.
Figure 2. Distribution of the lithological units and their horizontal hydraulic conductivity values in layer 1 of the model domain.
Water 14 02462 g002
Figure 3. Head Scatter Diagrams of the calibration (a panel) and validation processes (b panel), and comparison between calculated and observed heads at piezometers P1, P2, P3 and P4 (cf panels), represented in Figure 1.
Figure 3. Head Scatter Diagrams of the calibration (a panel) and validation processes (b panel), and comparison between calculated and observed heads at piezometers P1, P2, P3 and P4 (cf panels), represented in Figure 1.
Water 14 02462 g003
Figure 4. Three-dimensional sketch (a panel) and plan view (b panel) of the model calculated hydraulic heads map at the end of the simulation time.
Figure 4. Three-dimensional sketch (a panel) and plan view (b panel) of the model calculated hydraulic heads map at the end of the simulation time.
Water 14 02462 g004
Figure 5. Comparison of the fluxes’ trends over the simulation time among the basic (red line), the transient (blue line) and the S1, S2, S3, and S4 (green dashed lines) simulations: (a) Inflow from the fractured aquifer, (b) Evapotranspiration, (c,d) Seawater intrusion and Submarine groundwater discharge, (e,f) River inflow and outflow.
Figure 5. Comparison of the fluxes’ trends over the simulation time among the basic (red line), the transient (blue line) and the S1, S2, S3, and S4 (green dashed lines) simulations: (a) Inflow from the fractured aquifer, (b) Evapotranspiration, (c,d) Seawater intrusion and Submarine groundwater discharge, (e,f) River inflow and outflow.
Water 14 02462 g005
Figure 6. Comparison of the average flow rates and storage over the simulation time between: Basic, Transient, S1, S2, S3, and S4 simulations.
Figure 6. Comparison of the average flow rates and storage over the simulation time between: Basic, Transient, S1, S2, S3, and S4 simulations.
Water 14 02462 g006
Figure 7. Left panel: comparison between the recharge rates obtained by PCR-GLOBWB v2.0 and BIGBANG, with their linear trends (solid lines); Right panel: comparison between the evapotranspiration rates obtain by MODIS products and BIGBANG, with their linear trends (solid lines). The coloured areas between the dashed lines represent the 95% confidence interval.
Figure 7. Left panel: comparison between the recharge rates obtained by PCR-GLOBWB v2.0 and BIGBANG, with their linear trends (solid lines); Right panel: comparison between the evapotranspiration rates obtain by MODIS products and BIGBANG, with their linear trends (solid lines). The coloured areas between the dashed lines represent the 95% confidence interval.
Water 14 02462 g007
Table 1. National and worldwide open online datasets used to set up the model.
Table 1. National and worldwide open online datasets used to set up the model.
DatasetParameterRepositoryURL
SoilGrain size distributionISRIChttps://maps.isric.org/ accessed on 17 July 2022
GeologyAquifer thicknessISPRAhttp://sgi2.isprambiente.it/viewersgi2 accessed on 17 July 2022
HydrogeologyHydraulic conductivityISPRAhttp://sgi2.isprambiente.it/viewersgi2 accessed on 17 July 2022
CORINE Land CoverLand useCOPERNICUShttps://land.copernicus.eu/pan-european/corine-land-cover accessed on 17 July 2022
SEALEVELEURPHYL4MY008068Sea levelsCMEMShttps://doi.org/10.48670/moi-00141 accessed on 13 June 2022
Hydrometric station dependentRivers stagesCFMhttp://centrofunzionale.regione.campania.it/#/pages/sensori/archivio-idrometrici accessed on 17 July 2022
MYD16A3GFEvapotranspirationMODIShttps://doi.org/10.5067/MODIS/MOD16A2GF.006 accessed on 10 May 2022
PCRGLOBWB 2RechargeUU-HYDROhttps://github.com/UU-Hydro/PCR-GLOBWB_model accessed on 12 June 2022
BIGBANG40EvapotranspirationISPRAhttps://groupware.sinanet.isprambiente.it/bigbang-data/library/bigbang40/grids/it_evapotraspirazione_ae accessed on 17 July 2022
BIGBANG40RechargeISPRAhttps://groupware.sinanet.isprambiente.it/bigbang-data/library/bigbang40/grids/it_ricarica_degli_acquiferi_gr accessed on 17 July 2022
Table 2. Ss and Sy depending on the geologic material (upper panel), and Criv gained after the calibration procedure (lower panel).
Table 2. Ss and Sy depending on the geologic material (upper panel), and Criv gained after the calibration procedure (lower panel).
MaterialSs (m−1)Sy (-)
Clayey and calcareous clayey units0.0010.2
Debris, alluvial, fluviolacusters and fluvioglacial deposits1 × 10−40.33
Lavas, ignibrites, pyroclastics1 × 10−50.21
Organogenic, neritic, and platform limestone1 × 10−50.14
RiverCriv (m2/s)
Gariglianofrom 0.021 to 0.235
Savonefrom 0.017 to 0.058
Agnenafrom 0.001 to 0.033
Volturnofrom 0.009 to 0.044
Regi Lagnifrom 0.007 to 0.072
Sarnofrom 0.007 to 0.040
Picentinofrom 0.009 to 0.011
Tuscianofrom 0.011 to 0.019
Selefrom 0.014 to 0.116
Calorefrom 0.060 to 0.086
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gaiolini, M.; Colombani, N.; Busico, G.; Rama, F.; Mastrocicco, M. Impact of Boundary Conditions Dynamics on Groundwater Budget in the Campania Region (Italy). Water 2022, 14, 2462. https://doi.org/10.3390/w14162462

AMA Style

Gaiolini M, Colombani N, Busico G, Rama F, Mastrocicco M. Impact of Boundary Conditions Dynamics on Groundwater Budget in the Campania Region (Italy). Water. 2022; 14(16):2462. https://doi.org/10.3390/w14162462

Chicago/Turabian Style

Gaiolini, Mattia, Nicolò Colombani, Gianluigi Busico, Fabrizio Rama, and Micòl Mastrocicco. 2022. "Impact of Boundary Conditions Dynamics on Groundwater Budget in the Campania Region (Italy)" Water 14, no. 16: 2462. https://doi.org/10.3390/w14162462

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