The phenomenon of seawater intrusion is becoming increasingly problematic, particularly in low-lying coastal regions and areas that rely heavily on aquifers for their freshwater supply. It is, therefore, vital to address the causes and consequences of this phenomenon in order to ensure the security of water resources and the sustainable use of water. The objective of this paper was twofold: firstly, to delineate zones with different salinization levels over time; secondly, to investigate the factors controlling seawater intrusion of the Nile Delta aquifer. Aquifer data were collected in Sharkia governorate, Egypt, over three historical periods of years: 1996, 2007, and 2018. The dataset used to create the linear model of coregionalization consisted of hydrogeological (water level), hydrodynamic (pH, EC, Na, Mg, K, Ca, HCO3, SO4), and auxiliary (distances from salt and freshwater sources) variables. Cokriging was applied to produce spatial thematic maps of the studied variables for the three years of the survey. In addition, factorial cokriging was applied to understand the processes beyond the change in the aquifer water quality and map the zones with similar characteristics. Results of mapping the first factor at long range over the three years indicated that there was an increase in seawater intrusion, especially in the northeastern part of the study area. The main cause of aquifer salinization over time was the depletion of the groundwater resource due to overexploitation. @en

The Gran Sasso carbonate aquifer is the largest and most productive in the Apennines. Its hydrogeological structure has been studied since the middle of the last century for the springs’ characterization for drinking purposes and for a motorway tunnel. Meanwhile, its hydrodynamic parametrization is less developed and has been limited to monitoring the discharge and chemical and isotopic parameters. Secondary porosity characterizes the aquifer, and an underlying impermeable marly complex represents the basal aquiclude. It might appear inappropriate to characterize the hydraulic properties via pumping tests, as their reliability has been proven in homogeneous and isotropic media. However, the high extent of the aquifer, the wells’ location, the scarcity of information available and the lack of alternatives has forced the estimation of hydrodynamic parameters as in porous aquifers and the experimental testing of the aquifer, especially in maximum pumping conditions, for a possible exploitation increase. Since aquifer testing was performed during the normal well field’s activities, it was not possible to perform typical tests. Therefore, the step-drawdown test was conducted by turning on an increasing number of wells over time and keeping the observation points fixed. As results, a mean hydraulic conductivity of 5 × 10−3 m/s and a mean transmissivity of 0.3 m2/s were established without interrupting the water supply; meanwhile, the influence radius and flow directions were also estimated.@en

Recently, new numerical methods have been applied to weather data for the estimation of water budget, especially when the lack of measured data is considerable. Geostatistics is one of the most powerful approaches when it comes to studying spatially relevant natural phenomena, as it considers the spatial correlation among measurements over a specific study area and provides the associate uncertainty. In this study, we tested the feasibility of using a geostatistical method to provide a reliable estimation of the water budget of the Foro river basin (Central Italy) by comparing the obtained results with those of a traditional yet robust method. The results obtained with the geostatistical approach proved to be in line with the ones from the traditional method. Additionally, it was possible to quantify the uncertainty associated with the discharge values, making the estimates more reliable than the ones obtained with the traditional approach. However, the yearly distribution of river discharge obtained using both methods appeared to be dissimilar to the measured ones. The surface water uses, as well as the regulatory effect of the carbonate and alluvial aquifer regime, may affect the river discharge variability over the year and then can account for similar discrepancies between the inflow and outflow water volumes.@en

Rain gauge spatial sparsity and temporal discontinuity of data represent one of the major issues for reliable recharge estimations. In the past decades, the use of ground-based microwave weather RaDAR has dramatically improved quantitative rainfall estimation by providing spatially continuous estimates of rainfall over an area of more than 400 km² every 10 minutes. Furthermore, weather RaDAR data have also proved relatively reliable in mountainous areas. These paramount features of RaDAR-derived precipitation data could improve the estimation of potential recharge of aquifers, which rely on geospatializations (e.g., Thiessen polygons) of rainfall data collected by a sparse rain gauge network which often shows lacking at high altitude (i.e., recharge areas), introducing additional uncertainty in the inflow volumes. Weather RaDAR rainfall estimation is also affected by various sources of error, which can be reduced by proper post-processing; however, uncertainties remain, especially for surface rain rate estimations. Despite the currently necessary complex numerical processing, the purpose of the study is to evaluate the use of the weather RaDAR data as an alternative or in addition to meteorological data. Based on the above considerations, the feasibility of using RaDAR-based precipitation data to estimate aquifer potential recharge and calculate a detailed water budget in the areas characterized by high elevations, such as the Majella massif in the central Apennines, has been evaluated. To address this objective, the water budget has been calculated in the 2017-2018 period using both RaDAR-based precipitation data and rain gauge data, as well as adopting different methods (i.e., Turc and Thornthwaite). Although intrinsically uncertain, the RaDAR-based precipitation data provided solid results, pointed out by comparing it with water budget obtained by rain gauge data, and especially with experimental literature data. This interdisciplinary work may pave the way for continuous monitoring of aquifer potential recharge at extremely high temporal and spatial resolution.

@en

The concept of natural background level (NBL) aims at distinguishing the natural and anthropogenic contributions to concentrations of specific contaminants, as groundwater management and protection tools. This is usually defined as a unique value at a regional scale, even when the hydrogeological and geochemical features of a certain territory are far from homogeneous. The concentration of target contaminants is affected by multiple hydrogeochemical processes. This is the case of arsenic in the Calabria region, where concentrations are definitely variable in groundwater. To overcome the limitation of a traditional approach and to include the intrinsic hydrogeological and geochemical heterogeneity into the definition of the natural contribution to As content in groundwater, an integrated probabilistic approach to the NBL assessment combining aquifer-based preselection criteria and multivariate non-parametric geostatistics was proposed. In detail, different NBL values were selected, based on the aquifer type and/or hydrogeochemical features. Then, these aquifer-based NBL values of arsenic were used in the Probability Kriging method to map the probability of exceedance and to provide contamination risk management tools. This multivariate geostatistical approach that takes advantage of the physico-chemical variables used in the aquifer-based NBL values definition allowed mapping the probability of exceedance of As in a physically-based way. The hydrogeochemical diversity of the study area and all the processes affecting As concentrations in the aquifers have been considered too. As a result, the obtained map was characterized by a short-range and long-range variability due to local hydrogeochemical anomalies and water-rock interaction and/or atmospheric precipitation. By this approach, the NBL exceedance probability maps proved to be less “noisy”, because the local hydrogeochemical conditions were filtered, and more capable of pointing out anthropogenic inputs or very anomalous natural contributions, which need to be investigated more in detail and properly managed.

@en

Accurate knowledge of the rain amount is a crucial driver in several hydrometeorological applications. This is especially true in complex orography territories, which are typically impervious, thus, leaving most mountain areas ungauged. Due to their spatial and temporal coverage, weather radars can potentially overcome such an issue. However, weather radar, if not accurately processed, can suffer from several limitations (e.g., beam blocking, altitude of the observation, path attenuation, and indirectness of the measurement) that can hamper the reliability of the rain estimates performed. In this study, a comparison between rain gauge and weather radar retrievals is performed in the target area of the Abruzzo region in Italy, which is characterized by a heterogeneous orography ranging from the seaside to Apennine ridge. Consequently, the Abruzzo region has an inhomogeneous distribution of the rain gauges, with station density decreasing with the altitude reaching approximately 1500 m a.s.l. Notwithstanding, pluviometric inflow spatial distribution shows a subregional dependency as a function of four climatic and altimetric factors: coastal, hilly, mountain, and inner plain areas (i.e., Marsica). Such areas are used in this analysis to characterize the radar retrieval vs. rain gauge amounts in each of those zones. Compared to previous studies on the topic, the analysis presented the importance of an accurate selection of the climatic and altimetric subregional areas where the radar vs. rain gauge comparison is undertaken. This aspect is not only of great importance to correct biases in radar retrieval in a more selective way, but it also paves the way for more accurate hydrometeorological applications (e.g., hydrological model initialization and quantification of aquifer recharge), which, in general, require the accurate knowledge of rain amounts upstream of a basin. To fill the gap caused by the uneven rain gauge distribution, ordinary Kriging (OK) was applied on a regional scale to obtain 2D maps of rainfall data, which were cumulated on a monthly and yearly basis. Weather radar data from the Italian mosaic were also considered, in terms of rain rate retrievals and cumulations performed on the same time frame used for rain gauges. The period considered for the analysis was two continuous years: 2017 and 2018. The output of the elaborations included raster maps for both radar and interpolated rain gauges, where each pixel contained a rainfall quantity. Although the results showed a general underestimation of the weather radar data, especially in mountain and Marsica areas, they were within the 95% confidence interval of the OK estimation. Our analysis highlighted that the average bias between radar and rain gauges, in terms of precipitation amounts, was a function of altitude and was almost constant in each of the selected areas. This achievement suggests that after a proper selection of homogeneous target areas, radar retrieval can be corrected using the denser network of rain gauges typically distributed at lower altitudes, and such correction can be extended at higher altitudes without loss of generality.@en

The hilly landscape of the Periadric area in Central Italy is characterized by mainly marly– clayey foredeep basin deposits (Plio–Pleistocene age). These lithotypes are generally considered aquicludes, if compared with the regional limestone aquifers of Apennines. However, a coarsening upward trend characterizes the upper portion of this stratigraphic sequence, with arenaceous deposits and even conglomerates on the top. From a geomorphological viewpoint, the areas with coarser outcrops show a flat shape and sub‐vertical slopes, like boundaries. At the base of these scarps, springs can be found at the interface between coarse and fine deposits, whereas these arenaceous bodies are actual aquifers. Until now, the hydrodynamics and hydrochemical features of this kind of aquifer have not been investigated deeply, because they have always been considered a worthy water resource. However, they could play a crucial role in integrated water management, especially to cope with climate changes and drought periods. Considering these, the main purpose of this study was to investigate from a hydrogeological point of view and to assess the groundwater quantity and quality. Five examples throughout the Abruzzo region were considered. For evaluation and comparisons between water resources, the water volume that infiltrates yearly at each squared kilometer of an aquifer (Mm³/y/km²) was applied. This value was calculated through three different approaches to provide a recharge estimation for this kind of aquifer that is as exhaustive and representative as possible. The results allowed us to characterize the hydrogeological model and to quantify the resources between 0.1 and 0.16 Mm³/y/km², to be suitable for multi–purpose utilization.

@en

The Empirical Rainfall Thresholds (ERTs) for shallow landslide initiation are commonly devised worldwide mostly to be implemented within landslide early warning systems. Nonetheless, since the pioneering works on ERTs in the 1980s, only meteorological variables - that are cumulated E or intensity I and duration D values of rainfalls that are likely to trigger landslides - have been used to predict landslide occurrence, even though they are characterized by a large uncertainty. Over time, many efforts have been devoted to constrain ERT to geo-morphological characters of the landslide locations but, since nowadays, they did not get to a sound new method to derive ERT and strengthen its ability to forecast future rainfall-induced landslide. In this study, local geo-morphological characters have been taken into account by means of the co-kriging technique to constrain the E and D mean values of a regional ERT and their confidence intervals. The study area, where the proposed method was trained, is the hilly side of the Abruzzo region (Italy). Here, 62 shallow landslides have been analyzed in the time span of 2013–2017 by collecting 62 (D,E) pairs related to the rainfalls that were likely to trigger them. The relevant geo-morphological features for the considered territory have been selected through the principal component analysis. Then, the Multi-Collocated Co-Kriging technique, through ISATIS Geovariances software, has been applied to derive the spatial variability structures of E and D values conditioned by the selected geo-morphological parameters. Therefore, threshold values of E and D and their confidence intervals have been calculated generating a new shape of regional ERT, consisting of maps of continuous estimated threshold values of (D,E) and confidence interval values suitable for being used in early warning systems for shallow landslide initiation.

@en

Understanding the groundwater flow in carbonate aquifers represents a challenging aspect in hydrogeology, especially when they have been struck by strong seismic events. It has been proved that large earthquakes change springs hydrodynamic behaviour showing transitory or long-lasting variations and making their management much more difficult. This is the case of Sibillini Massif (central Italy), which has been hit by the well-known 2016–2017 seismic period. This work aims to improve the knowledge of carbonate aquifers groundwater circulation and their possible changes in the hydrodynamic behaviour, during and after a series of strong seismic events. The goal has been achieved by comparing long-time tracer tests and transient time-series analysis, based on a sliding-window approach. This approach allowed investigating transient variations in the carbonate aquifers recharge system, highlighting the changes of relationships between the inflow contributions to the spring discharge in the area. As a result, the seismically triggered pore pressure distribution, and the hydraulic conductivity variations, because of the ground shaking and the fault systems activation, account for all the mid-and long-term modifications in the recharge system of Sibillini aquifers, respectively. These outcomes provide valuable insights to the knowledge of aquifer response under similar hydrogeological conditions, that are vital for water management.

@en

The incorrect wastewater management and the land use distribution lead to severe environmental problems, creating heavy eutrophication condition in surface-water; when surface-water/groundwater relationships exist, the organic matter transferred to the aquifer oxidizes and triggers redox processes (i.e. Terminal Electron Accepting Processes, TEAPs), that provoke severe groundwater quality modifications and complicate its exploitation and management. For this reason, the definition of the redox zonation within an aquifer is an effective tool for the identification of the contamination sources and for the conceptual model refinement, when remediation strategies need to be planned. Although the redox processes are dynamic reactions, the aquifer redox zonation is generally aimed to identify homogenous zones, characterized by a predominant TEAP. To overcome this methodological approach, the Multi-Collocated Factorial Kriging (MCFK) was applied to redox-related physico-chemical parameters, which allowed identifying their spatial relationships at different scales, transferring this method from precision agriculture and soil science to hydrogeochemistry. The selected study area is the San Pedro Sula aquifer (Honduras), a multi-layer alluvial aquifer characterized by well-known surface-water/groundwater interactions and heavy eutrophicated streams. Here, high concentrations of Mn and Fe were found in the aquifer. The MCFK results identified a short-range (2300 m) factor, highlighting a strong relation between Mn concentrations and anoxic conditions, due to the organic matter transfer from eutrophicated surface-water into the aquifer. Simultaneously, the relationship between Fe and turbidity is related to a fine Fe(III) oxi-hydroxide colloidal phase, developed when different redox conditions of groundwater mix up in the wells. The long-range (6000 m) factor points out that Fe is related to redox processes at a wider scale, especially in the northern San Pedro Sula alluvial plain. These results are supported by both the Principal Component Analysis and the hydrogeochemical numerical modeling. As a result, different TEAPs occur simultaneously in contaminated areas, acting at multiple scales.

@en

The wide regional aquifer system springs are one of the most important sources of freshwater, especially because of the ever-increasing global water demand. Thus, these springs must be protected and managed in a sustainable manner. To this purpose, a detailed knowledge about the corresponding aquifer systems hydrodynamic is vital. In this general framework, the Majella carbonate aquifer is one of the most important as well as heterogenous aquifer systems in Central Italy. It is characterized by high elevation, a heterogeneous hydraulic conductivity distribution and groundwater flow, and its recharge is mainly due to both rainfall and snowmelt. Its most important spring is the Verde spring, located in the eastern side, exploited for drinking and hydroelectric purposes. The main objectives of this research are (1) to get a deeper insight into the hydrogeological conceptual model of the Verde spring, combining multiparameter input and output time-series, and (2) to quantify each inflow effect on the Verde spring parameters by a methodological approach that is quite uncommon in literature: comparing univariate and bivariate statistical analyses of both raw and residual (non-systematic component) time-series. The results obtained highlighted that the snowmelt is the most important inflow of the Majella aquifer system. The snow cover melting leaks slowly large amount of water into the aquifer and creates smoothly significative modifications to the spring discharge, electrical conductivity and temperature, in terms of groundwater volume increase and dilution. Contrariwise, although the rainfall inflow volume and the corresponding spring parameters changes are very limited, the transient behavior of the rainfall inflow allowed identifying different recharge modes. These recharge modes depend on several flow paths in the unsaturated zone, characterized by different water volumes brought toward the saturated zone, sizes, hydraulic conductivities, and distances from the Verde spring.

@en

The water–rock interaction is discussed in this study for some Central Apennine aquifers and their relative springs, where the geological and hydrogeological setting is potentially responsible for hydrocarbon contamination. The contamination is related to the presence of limestone formations with high organic content that are connected to the genesis of hydrocarbons in the Central-Southern Apennines. Thanks to a multidisciplinary approach based on a seasonal monitoring of hydrogeological, hydrogeochemical, organic chemical, and isotopic variables, and to a detailed geological analysis, it was possible to demonstrate that the bituminous schists within the bituminous dolomite formation (a Triassic Formation presents in the Gran Sasso acquifer), the bituminous intercalations within the Bolognano Formation (an Olocenic calcareous Formation in the Majella acquifer), and the petroleum-saturated rocks of the Bolognano Formation (in the Morrone acquifer) are all able to leach hydrocarbons into groundwater. The results suggest that more detailed studies are required for areas where lithologies with fossil organic components are present. Insights should further investigate the interaction between groundwater and rocks in terms of organic compounds as well as inorganic compounds. In particular, the study also suggests that the supplementary quantification of hydrocarbon compounds in bituminous rock and the leaching tests are analyses that provide reliable results. From the normative point of view, the results of the study can be useful when dealing with hydrocarbon contamination resulting from anthropogenic activities within catchments where formations with high concentration of organic matter are present. In these cases, it will possible to assess the natural background concentrations and review the threshold values upwards.

@en

In this research, a multi-scenario numerical modeling was implemented to assess the effects of changes to abstraction patterns in the Sant'Angelo well-field (central Italy) and their implications on the aquifer hydrodynamic and the advective transport of contaminants. Once implemented and calibrated the steady-state numerical model by means of MODFLOW-2005, the well-field turning off scenario was modelled. In addition, the numerical results were analyzed by means of the post-processors ZONEBUDGET and MODPATH, to assess respectively the contribution of each hydrogeological feature to the total budget and the advective transport of contaminant particles. Comparing the two steady-state numerical models and the relative particle tracking analyses, the well-field turning off, although no longer acting as a hydraulic barrier, increased the residence time of contaminant particles and limited their mobility in the aquifer. Furthermore, the general decrease in groundwater abstractions also caused a higher increase in river flow, favoring contaminants' dilution in surface water.

@en

In this research, univariate and bivariate statistical methods were applied to rainfall, river and piezometric level datasets belonging to 24-year time series (1986-2009). These methods, which often are used to understand the effects of precipitation on rivers and karstic springs discharge, have been used to assess piezometric level response to rainfall and river level fluctuations in a porous aquifer. A rain gauge, a river level gauge and three wells, located in Central Italy along the lower Pescara River valley in correspondence of its important alluvial aquifer, provided the data. Statistical analysis has been used within a known hydrogeological framework, which has been refined by mean of a photo-interpretation and a GPS survey. Water-groundwater relationships were identified following the autocorrelation and cross-correlation analyses. Spectral analysis and mono-fractal features of time series were assessed to provide information on multi-year variability, data distributions, their fractal dimension and the distribution return time within the historical time series. The statistical-mathematical results were interpreted through fieldwork that identified distinct groundwater flowpaths within the aquifer and enabled the implementation of a conceptual model, improving the knowledge on water resources management tools.

@en

The present work aims to assess the origin of the boron contamination of groundwater that affect some areas of central-Adriatic Italy, through the use of boron and strontium isotopic analysis. The research results indicate that boron contamination is mainly related with salinization of groundwater attributable to the uprising of connate groundwater and mixing processes between more shallow fresh groundwater and deeper saline groundwater. As a result, the main boron contamination of groundwater is not related to anthropogenic inputs but rather to natural processes.

@en

A multidisciplinary approach has been used, integrating statistical multivariate analysis and hydrogeochemical modeling, in order to evaluate processes that cause high concentration of iron and manganese in drinking water wells exploiting the San Pedro Sula porous aquifer. The analyses carried out have allowed to individuate a close relationship between the heavy eutrophication conditions of surface water, due to uncorrect management of waste water and to land use impact, and redox processes diffusion in groundwater, that cause Fe and Mn mobilisation.

@en

Iron and manganese are two of the most common contaminants that exceed the threshold imposed by international and national legislation. When these contamination occurs in groundwater, the use of the water resource is forbidden for any purposes. Several studies investigated these two metals in groundwater, but research focused in the Central Adriatic area are still lacking. Thus, the objective of this study is to identify the origin of Fe and Mn contamination in groundwater and the hydrogeochemical processes that can enrich aquifers with these metals. This work is based on hydrogeochemical and multivariate statistical analysis of analytical results undertaken on soils and groundwater. Fe and Mn contamination are widespread in the alluvial aquifers, and their distribution is regulated by local conditions (i.e. long residence time, presence of peat or organic-rich fine sediments or anthropic pollution) that control redox processes in the aquifers and favour the mobilisation of these two metals in groundwater. The concentration of iron and manganese identified within soil indicates that the latter are a concrete source of the two metals. Anthropic impact on Fe and Mn contamination of groundwater is not related to agricultural activities, but on the contrary, the contribution of hydrocarbons (e.g. spills) is evident.

@en