Search

Changes in groundwater flow in response to strong earthquakes are widely described in many tectonic environments. For example, a post-seismic discharge variation is often attributed to an increase of bulk permeability due to co-seismic fracturing and/or to a change in the role of faults in acting as conduits/barrier to groundwater flow.We take as an example the fractured aquifer of the Mts. Sibillini carbonate massif, in Central Italy, which were affected by a strong and prolonged extensional seismic sequence in 2016-17. The sequence was characterized by an M=6.5 event (mainshock), an M=6 event, an M=5.9 event, up to 60 M>4 events and several M>5 earthquakes. The strongest events caused rupturing of the topographic surface for a cumulative length in the order of 30 km and an important portion of aftershocks occurred at depths where groundwater is stored.As a response to the seismic sequence, the main NNW-directed groundwater flow was diverted to the west and a discharge deficit was observed at the foot-wall of the activated fault system with a relevant discharge increase, accompanied by geochemical variations, at the fault system hanging-wall.By integrating geo-structural reconstructions, seismological and ground deformation data, artificial tracer tests results and a 4-years discharge and geochemical monitoring campaign data, we show that the observed groundwater variations are due to a combination of permeability increase along the activated fault systems and hydraulic conductivity increase of the hanging-wall block due to fracturing, extension and subsidence, which determined a fast aquifers emptying. Seismicity temporarily triggered a change of the pre-existing predominant along-faults-strike NNW-SSE oriented regional flow to a west-directed flow, perpendicular to faults strike. We discuss the position of the aquifer with respect to the activated faults and how this affected the observed phenomena. REFERENCECambi, C., Mirabella, F., Petitta, M., Banzato, F., Beddini, G., Cardellini, C., ... & Valigi, D. (2022). Reaction of the carbonate Sibillini Mountains Basal aquifer (Central Italy) to the extensional 2016–2017 seismic sequence. Scientific Reports, 12(1), 1-13. DOI: 10.1038/s41598-022-26681-2

Understanding hydrological flow pathways and spatio-temporal origin of surface and subsurface water in catchments with highly fractured geology is particularly challenging. In this work, we relied on the integration of hydrometric measurements with stable oxygen and hydrogen isotope data in two adjacent catchments in the Sibillini Mountains National Park, Central Italy, to better understand the drivers of the catchment hydrological response and the spatio-temporal origin of stream and spring waters. The Ussita catchment is 44 km2 and its highest elevation is 2204 m a.s.l. The Nera catchment is 110 km2 and its highest peak is 2233 m a.s.l. The two rivers merge at the Visso Village, at 615 m a.s.l.The area is characterised by heavily fissured and fractured calcareous rocks that foster the occurrence of several springs, some of them of karst origin. Both catchments host a dense hydrometerological network. The experimental apparatus is completed with one piezometer, soil moisture probes at five locations, lysimeters at two depths and four throughfall plots under beeches and oaks. Monthly samples for isotopic analysis are being collected since fall 2020 from precipitation at three different elevations and four locations, the streams at different sections, and four springs in the Nera catchment only.Preliminary results show a distinct hydrological behaviour in the annual streamflow regimes: the Ussita stream slightly reacts only to the largest storms and during intense snowmelt periods, whereas the Nera stream has a very damped response during all the year, revealing a clear buffer effect of the large subsurface reservoir, facilitated by the highly fractured nature of the geological setting. As expected, there is an elevation and seasonal effect in the isotopic composition of precipitation, although the seasonal effect is partly masked by the exceptionally high temperatures occurred in fall 2021. However, the time series of isotope data in stream water show a damped signal and very low seasonal variability in both streams, matching the observed low variability of streamflow. Only the Ussita catchment shows some more enriched outliers likely reflecting runoff response during large storm events. Interestingly, stream and spring samples from both catchments lie along but also above and below the Local Meteoric Water Line, suggesting that the sampled spring and stream water was either originated i) from precipitation fell, infiltrated, and stored well before the collection of the precipitation samples, and released; ii) and/or from areas outside the topographic catchments, and therefore not adequately characterized by the isotopic signal of sampled precipitation. The isotopic composition of the streams and springs is statistically the same, revealing that spring groundwater is the main component of stream runoff. Moreover, the isotope signature of both springs and streams is much closer to that of winter precipitation rather than summer precipitation indicating a major role of winter precipitation in recharging the catchments, consistently with the precipitation seasonal regime. On-going work is assessing the spatial difference in the isotopic composition and quantifying the temporal origin of stream and spring water of the two catchments.

Over the years the scientific community underlined several problems related to the use of isotopic hydrology techniques in areas characterized by a complex orography, usually occurring in mountain areas, or in those cases where the hydrological setting is complicated by contacts between different aquifers (Nanni et al., 2013).In this study, an innovative isotopic model, able to identify the most probable recharge area for several springs exploited for drinking purposes, has been developed and applied to the Nera catchment in the Sibillini Mountain National Park (central Italy). The isotopic investigation consists of a preliminary definition of a new δ18O - elevation relationship, considering the morphological and meteorological heterogeneities within the area and their possible influences on the precipitation isotope values (e.g., shaded areas, snow drift effect, etc.). Second, an advanced δ18O distribution model, supported by statistical and GIS-based procedures, has been implemented by clipping the precipitation δ18O values (depicted from the δ18O – elevation relationship) over an upstream area for each analyzed spring. The new isotopic modeling approach can be conveniently applied if the infiltration rate of the meteoric water is fast enough to avoid fractionation processes that may alter the isotopic signal of the precipitation input within the aquifer, and if peculiar meteoric recharge phenomena, altering the springs' isotopic signal, are treated as outliers.This research highlights if the most used isotopic approach based on the determination of groundwater recharge areas starting from δ18O - elevation gradient (Jeelani et al., 2010; Jasechko, 2019) applied to a selected spring isotopic data agrees with the hydrogeological setting of the spring recharge area which is often complicated by the topography and the contacts between different aquifers both for stratigraphic and tectonic reasons.The ultimate goal of this study is to quantify the aquifers recharge under the impact of drought to improve the water resources management operations in the area.

In the framework of landslides, the hydrogeological features play an essential role in slope stability, governing water movement and thus resulting in modification of the effective stress in the soil. In this framework, the hydrogeological conceptualization of landslide areas and the identification of groundwater origin are key points to developing risk mitigation measures. In fact, groundwater recharge cannot always be attributed to local precipitation alone. Mixing processes between water derived from local infiltration and deep water upflow along tectonic lineaments or anthropogenic water can affect the groundwater balance on a local scale. This study aims to define the potential groundwater origin of one of the highest risk urban landslides in central Italy and to define a hydrogeological conceptual model by exploiting its existing drainage system network. This research is based on a multiple-techniques approach based on hydrological water balance, artificial tracer tests during recharge period, seasonal monitoring of the water stable-isotope content, hydro-chemical survey during low-flow periods, and analysis of the piezometric level fluctuation. All these analyses are coupled with a detailed reconstruction of the geology of the area depicted from boreholes and drill holes. Two groundwater bodies have been evidenced from the study. The shallower one is located in the landslide unstable zone and is hydraulically connected to a deeper groundwater body hosted in the underlying bedrock. Results highlighted that the local rainfall regime could not fully explain the hydro-chemical facies. Local water contributions to the landslide area coming from leakage of the urban sewerage system have been evidenced, excluding deep groundwater upflow from the fault system.

Placed at the center of the Mediterranean, the Apennines chains provide a critical water supply for people living in the Italian Peninsula. Yet, the quantification of water resources in this region is challenging given that the different components of the water cycles (i.e., snowmelt, evapotranspiration, and subsurface water recharge) are highly variable in space and time due to the specificity of the climate, the reforestation trend, and the complex landscapes and geology. In this study, we investigated the challenging hydrological river regime of a complex carbonate basin with significant external (and partially karst) groundwater contribution – the Upper Nera basin – affected by recent important seismic sequences. When dealing with such type of basins, the generic approach to delineate the basin boundaries based on the geomorphology of the area can lead to questionable results potentially yielding significant water budget imbalances. Therefore, both (hydro)geological and hydrological features have to be considered for understanding the challenging hydrological behavior of these basins. Here we proposed a specific analysis of precipitation-runoff time series corroborated with hydro-geological survey to obtain information on the basin response time and the real contribution area of the basin. Then, we integrated this information within the structure of a hydrological model (i.e., the Geoframe modelling system) by using a parsimonious type of approach to improve the description of the different components of the water balance with a specific focus on external groundwater recharge. We validated the model against in situ discharge observations and with remote sensing information of evapotranspiration and snow. Since these two variables (snow cover and evapotranspiration) play the crucial role in correct estimation of water balance in the basin, in this study we additionally focused on that. We show that the model (tested with several hydrological signatures and a new conceptual evaluation based on an empirical probability function) performs relatively well in reproducing the different water balance components (including remotely sensed evapotranspiration and snow) and that the upper river basin is significantly impacted by carbonate rock river discharge (i.e., up to 85 % in proportion to the total discharge for some stations) coming from outside the geomorphological boundary of the basin. Yet, the groundwater recharge effects on the river, gradually attenuates at the outlet of the basin (Visso station).

Using artificial tracer tests, this study confirms the presence of a single basal aquifer feeding the springs in the wide and complex hydrogeological boundary of the Sibillini Mountains of central Italy. The tracer was introduced into the sinkhole of the Castelluccio di Norcia plain. The tracer test results, observed at the studied springs, highlight the changes induced by the 2016/2017 earthquake in the water circulation of the aquifer system. In particular, the seismic events increased the hydraulic conductivity of the basal aquifer, with a consequent increase in the springs’ flow rates at the western hydrogeological boundary and a decrease in the flow rates at the eastern Adriatic hydrogeological boundary. This phenomenon is in accordance with the hydro-structural framework of the area. The study also investigated the relation between groundwater circulation and tracer behaviour in the springs during the pre- and post-earthquake periods. The tracer test results led to the formulation of hypotheses about water circulation of the area. The trend of the tracer breakthrough curves demonstrates that the upper portion of the basal aquifer is characterised by fast water circulation similar to that in the epiphreatic area of karst aquifers, while a slow circulation due to fissures with interconnected drains occurs in the deepest portion of the aquifers. The obtained results highlight how tectonics and karst can affect the hydrogeological setting of the Apennine carbonate chain; in particular, the seismicity of central Italy may alter groundwater circulation for a long period of time when great magnitude earthquakes occur.

The analysis of the water budget in the Upper Nera River basin, typical karst and fissured rocks catchment located in the Apennines chains in central Italy, has been performed to simulate snow, Evapotranspiration (ET), and runoff through different components of the GEOframe system. During this study, we showed that using an unsupervised approach for extracting the basin boundary could provide significant issues in the correct estimation of water budget components. To overcome this problem, both hydrogeological and hydrological information -obtained through a new type of time-series analyses and recent geological surveys- have been considered to estimate the contribution area and time response of the karst discharge. According to the mentioned information and benefiting the flexibility of the GEOframe-NewAge modeling system, a conceptual reservoir with a 30-day time response, derived from the time series analysis, has been added to estimate the karst river discharge contributed to Nera. The model, evaluated by different signatures (including mean daily discharge, high flow, low flow, low flow duration frequency, and flow duration curve slope and a new empirical probability function) has been shown to reproduce the water fluxes of the hydrological cycle in the basin relatively well (KGE values equal to 0.61, 0.80, and 0.71 in different sections, respectively). The karst discharge flux has a significant effect on the water budget of the basin especially in the upstream part (Castelsantangelo section) and this effect decreased through the river downstream to the outlet of Visso. We showed that 85% of the total discharge at Castelsantangelo station comes from outside of the geomorphological boundary of the basin. According to the water balance analysis, the maximum karst flux that happened in 2014 could be mainly caused by the maximum precipitation that happened in 2013 over the basin. To further cross-validation of the model performance, MODIS ET and Sentinel-1 snow depth products were used. The comparison of remote-sensed MODIS ET and GEOframe ET shows a systematic difference, with higher values of MODIS ET than our model estimations. As well, the spatial correlation of snow cover retrieved from Sentinel-1 snow depth and GEOframe Snow Water Equivalent has been examined and a good correlation has been reported especially for Castelsantangelo. The values of Sentinel-1 were also verified through some in-situ snow depth data.

As the rise in pore water pressure is one of the main factors in triggering landslides, the understanding of groundwater processes taking place at the hillslope scale is a crucial issue in slope stability analysis. However, identifying flow-paths travelled by water molecules from their infiltration is still a complex task. Hydrochemistry is recognized as a powerful tool that can help to gain useful hydrogeological information and has gradually become increasingly used in addition to conventional study methods. This manuscript presents a comprehensive geochemical investigation consisting of leaching tests and quantitative mineralogical analyses on soil samples, chemical analyses on groundwater samples and modelling. Our results highlighted the usefulness of coupling, even in hydrogeological studies focusing on landslides, geochemical surveys on both water and the soil matrix to constrain the interactions between host-rocks and groundwater. Moreover, it demonstrated that kinetic-based geochemical models, if properly calibrated on leaching tests, can provide valuable information on groundwater dynamics, allowing us to elucidate water-mixing processes beneath the soil surface.

Surface water, groundwater, phytosanitary products, soil, isotope., Flowpath 2017, National meeting on hydrogeology 2017

Dye tracer tests are a good technique useful to determine the hydrogeological parameters of aquifers, especially in alluvial plains, where the heterogeneity of the lithology is a key factor leading the groundwater circulation. In the last decade, new tracers have been developed and introduced in the hydrological and geological context, including DNA tracers, that are used in the same way as the traditionaltracers. In this study, some experiments were performed in a porous aquifer with the aim of comparing the behavior of different tracers for the evaluation of aquifer’s parameters, and of testing the suitability of a DNA tracer in this geological context. Single well and multiwell tests were performed during the winter and spring 2017, involving one piezometer and 4 wells. The tracer test results highlighted fasterzones of flow in the alluvial plain, being located at different depths, whose levels are characterized by gravel and sand layers. A maximum groundwater velocity of around 170 cm/day was revealed by the single well tests, whereas the main direction of flow was along the watercourse axis, as resulted by the multi well tracer test.The overall results obtained confirmed that fluorescent tracers are a very valuable and reliable tool for the determination of groundwater velocity, groundwater flow direction and dispersion parameters, even in cases of multi-layered aquifers, whereas DNA tracers are very helpful in the determination of faster circuits and vertical flow even in multi-layered aquifers.

This Editorial paper sums up the contents of the Special Issue named “Applying Artificial and Environmental Tracing Techniques in Hydrogeology” [...]

The interaction between fluids and tectonic structures such as fault systems is a much-discussed issue. Many scientific works are aimed at understanding what the role of fault systems in the displacement of deep fluids is, by investigating the interaction between the upper mantle, the lower crustal portion and the upraising of gasses carried by liquids. Many other scientific works try to explore the interaction between the recharge processes, i.e., precipitation, and the fault zones, aiming to recognize the function of the abovementioned structures and their capability to direct groundwater flow towards preferential drainage areas. Understanding the role of faults in the recharge processes of punctual and linear springs, meant as gaining streams, is a key point in hydrogeology, as it is known that faults can act either as flow barriers or as preferential flow paths. In this work an investigation of a fault system located in the Nera River catchment (Italy), based on geo-structural investigations, tracer tests, geochemical and isotopic recharge modelling, allows to identify the role of the normal fault system before and after the 2016–2017 central Italy seismic sequence (Mmax = 6.5). The outcome was achieved by an integrated approach consisting of a structural geology field work, combined with GIS-based analysis, and of a hydrogeological investigation based on artificial tracer tests and geochemical and isotopic analyses.

Nowadays, the balance between incoming precipitation and stream or spring discharge is a challenging aspect in many scientific disciplines related to water management. In this regard, although advances in the methodologies for water balance calculation concerning each component of the water cycle have been achieved, the Thornthwaite–Mather method remains one of the most used, especially for hydrogeological purposes. In fact, in contrast to physical-based models, which require many input parameters, the Thornthwaite–Mather method is a simple, empirical, data-driven procedure in which the error associated with its use is smaller than that associated with the measurement of input data. The disadvantage of this method is that elaboration times can be excessively long if a classical MS Excel file is used for a large amount of data. Although many authors have attempted to automatize the procedure using simple algorithms or graphical user interfaces, some bugs have been detected. For these reasons, we propose a WebApp for monthly water balance calculation, called WaterbalANce. WaterbalANce was written in Python and is driven by a serverless computing approach. Two respective European watersheds are selected and presented to demonstrate the application of this method.

This study wants to give a contribution to the semi-automatic evaluation of rock mass discontinuities, orientation and spacing, as important parameters used in Engineering. In complex and inaccessible study areas, a traditional geological survey is hard to conduct, therefore, remote sensing techniques have proven to be a very useful tool for discontinuity analysis. However, critical expert judgment is necessary to make reliable analyses. For this reason, the open-source Python tool named DCS (Discontinuities Classification and Spacing) was developed to manage point cloud data. The tool is written in Python and is based on semi-supervised clustering. By this approach the users can: (a) estimate the number of discontinuity sets (here referred to as “clusters”) using the Error Sum of Squares (SSE) method and the K-means algorithm; (b) evaluate step by step the quality of the classification visualizing the stereonet and the scatterplot of dip vs. dip direction from the clustering; (c) supervise the clustering procedure through a manual initialization of centroids; (d) calculate the normal spacing. In contrast to other algorithms available in the literature, the DCS method does not require complex parameters as inputs for the classification and permits the users to supervise the procedure at each step. The DCS approach was tested on the steep coastal cliff of Ancona town (Italy), called the Cardeto–Passetto cliff, which is characterized by a complex fracturing and is largely affected by rockfall phenomena. The results of discontinuity orientation were validated with the field survey and compared with the ones of the FACETS plug-in of CloudCompare. In addition, the algorithm was tested and validated on regular surfaces of an anthropic wall located at the bottom of the cliff. Eventually, a kinematic analysis of rock slope stability was performed, discussing the advantages and limitations of the methods considered and making fundamental considerations on their use.

A large tank (1.4 m x 4.0 m x 1.3 m) filled with medium-coarse sand was employed to measure evaporation rates from shallow groundwater at controlled laboratory conditions, to determine drivers and mechanisms. To monitor the groundwater level drawdown 12 piezometers were installed in a semi regular grid and equipped with high precision water level, temperature, and electrical conductivity (EC) probes. In each piezometer, 6 micro sampling ports were installed every 10 cm to capture vertical salinity gradients. Moreover, the soil water content, temperature and EC were measured in the unsaturated zone using TDR probes placed at 5, 20 and 40 cm depth. The monitoring started in February 2020 and lasted for 4 months until the groundwater drawdown became residual. To model the groundwater heads, temperature, and salinity variations SEAWAT 4.0 was employed. The calibrated model was then used to obtain the unknown parameters, such as: maximum evaporation rates (1.5-4.4 mm/d), extinction depth (0.90 m), mineral dissolution (5.0e-9 g/d) and evaporation concentration (0.35 g/L). Despite the drawdown was uniformly distributed, the increase of groundwater salinity was rather uneven, while the temperature increase mimicked the atmospheric temperature increase. The initial groundwater salinity and the small changes in the evaporation rate controlled the evapoconcentration process in groundwater, while the effective porosity was the most sensitive parameter. This study demonstrates that shallow groundwater evaporation from sandy soils can produce homogeneous water table drawdown but appreciable differences in the distribution of groundwater salinity.

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&ndash, 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.

Groundwater is the most used water resource around the world, but due to population growth and climate change the alluvial lowland aquifers are often polluted and over-exploited. Thus, more and more frequently water managers need to shift their attention to mountain regions to identify groundwater resources for drinking purposes. This study presents a monitoring and modelling approach that allowed to quantify the inflow from the &ldquo, Montagna dei Fiori&rdquo, fractured aquifer to the Castellano stream. Continuous monitoring of flow discharge and temperature during an entire hydrological year (2018&ndash, 2019) at two monitoring stations along the stream allowed to discriminate between the baseflow (on average, 0.891 m3/s) and the run-off (on average, 0.148 m3/s) components. A hydrogeological basin-wide numerical flow model (using MODFLOW-2005) was set up using information from hydrogeological and geomechanical surveys. The model was calibrated using the daily baseflow observations made in the Castellano stream (R2 = 0.75). The calibrated model allowed to quantify groundwater/surface water interactions. After an automated sensitivity analysis (using MODFLOW-2000), the recharge was found to be the most uncertain parameter, followed by the hydraulic conductivity zonation. This methodology could be applied in other mountain regions where groundwater monitoring networks are usually lacking to improve water resources management.

The dynamic of groundwater systems feeding several springs of the Sibillini Mountains was deeply affected by nine Mw 5.0÷6.5 seismic events occurred in central Italy starting from August 2016. The strongest shock occurred on October 30th 2016 about 5 km NNE of Norcia Town, 9 km below the surface, as a result of upper crust normal faulting on the nearly 30 km-long Mt Vettore - Mt Bove fault system, a NW-SE trending, SW-dipping fault system outcropping on the western slope of Mt Vettore, the highest peak of Sibillini Mountains. Soon after this event, a general increase of springs and rivers discharge and groundwater levels was observed both in the Visso and Norcia areas, west of the Sibillini Mountains. In the Visso area the hydrogeological changes due to the seismic sequence exhausted in the 2019, while nowadays both discharges and groundwater levels are still higher than before in the Norcia area. Discharge data of the main springs located east, south-east of the Sibillini Mountains were analysed to verify whether the general increase observed on the western side was associated to a decrease on the eastern and southern-east area. The results show that the springs located on the eastern side and southern-east side of Mt Vettore experienced a significant long-term discharge decrease. In this preliminary work, the analysis of the historical discharge series of the Pescara di Arquata spring (SE of Mt Vettore), and its relationship with the Standard Precipitation Index (SPI) shows that the very low discharge values recorded during the post-seismic period are not associated with SPI as low as documented in the past for similar discharges. Moreover, the stable isotopic composition of Pescara di Arquata water during the post-earthquake period is slightly different from that measured before the seismic events; this suggests that a lower amount of water having more enriched isotopic δ18O content reaches the spring after the seismic sequence. These aspects seem to indicate that groundwater circulation in the southern-east area of Sibillini Mountains has been affected by the 2016-2017 seismic sequence.

We compared five regression approaches, namely, ordinary least squares, major axis, reduced major axis, robust, and Prais–Winsten to estimate δ18O–δ2H relationships in four water types (precipitation, surface water, groundwater collected in wells from lowlands, and groundwater from low-yield springs) from the northern Italian Apennines. Differences in terms of slopes and intercepts of the different regressions were quantified and investigated by means of univariate, bivariate, and multivariate statistical analyses. We found that magnitudes of such differences were significant for water types surface water and groundwater (both in the case of wells and springs), and were related to robustness of regressions (i.e., standard deviations of the estimates and sensitiveness to outliers). With reference to surface water, we found the young water fraction was significant in inducing changes of slopes and intercepts, leading us to suppose a certain role of kinetic fractionation processes as well (i.e., modification of former water isotopes from both snow cover in the upper part of the catchments and precipitation linked to pre-infiltrative evaporation and evapotranspiration processes). As final remarks, due to the usefulness of δ18O–δ2H relationships in hydrological and hydrogeological studies, we provide some recommendations that should be followed when assessing the abovementioned water types from the northern Italian Apennines.

Deep-seated landslides are complex systems. In many cases, multidisciplinary studies are necessary to unravel the key hydrological features that can influence their evolution in space and time. The deep-seated Berceto landslide, in the northern Apennines of Italy, has been investigated in order to define the origin and geochemical evolution of groundwater (GW), to identify the slope system hydrological boundary, and to highlight the GW flow paths, transit time and transfer modalities inside the landslide body. This research is based on a multidisciplinary approach that involves monitoring GW levels, obtaining analyses of water chemistry and stable and unstable isotopes (δ18O-δ2H, 3H, 87Sr/86Sr), performing soil leaching tests, geochemical modelling (PHREEQC), and principal component analysis (PCA). The results of δ18O-δ2H and 87Sr/86Sr analyses show that the source of GW recharge in the Berceto landslide is local rainwater, and external contributions from a local stream can be excluded. In the landslide body, two GW hydrotypes (Ca-HCO3 and Na-HCO3) are identified, and the results of PHREEQC and PCA confirm that the chemical features of the GW depend on water–rock interaction processes occurring inside the landslide. The 3H content suggests a recent origin for GW and appears to highlight mixing between shallow and deep GW aliquots. The 3H content and GW levels data confirm that shallow GW is mainly controlled by a mass transfer mechanism. The 3H analyses with GW levels also indicate that only deep GW is controlled by a pressure transfer mechanism, and this mechanism is likely the main influence on the landslide kinematics.

The recharge settings of aquifers (as those of the carbonate ridges in central Apennines) is strongly dependent on hydrogeological and geo-structural complexity; in fact, geochemical data, hydrographs interpretation, hydrogeological and geo-structural surveys, tracer and isotope data are usually not so developed and quite difficult to put into relation each other. In this paper, the importance of relating isotope techniques to hydrogeological survey on the recharge area evaluation is demonstrated by studying the system of the Mt Conero limestone ridge, located in central Italy close to the Adriatic sea. A closed system, at a pilot scale, with sufficiently well defined outcropping areas of aquifer, which is offering a contribution to the aquifer recharge area evaluation. In this study a hydrogeological investigation was performed to identify the infiltration areas in the aquifer, analysing the entire geological formation divided into several members, each of them showing different hydrogeological behaviour. In addition, periodic sampling of groundwater and of precipitation at different elevations were carried out to study the isotopic composition over a period of about 4 years. Weighted mean and arithmetic mean of isotopes were used in the gradient calculation and compared giving actually different results. The mean recharge altitude of the aquifer, as evaluated by weighted mean of isotopes, is about 100 m higher than the average elevation of the aquifer outcrop; moreover, it was proven that the vast majority of the recharge is due to the more permeable members of the aquifer (about 20% of the outcropping rocks, located at higher elevation). The recharge mechanism of the aquifer is therefore depending both on the isotope values of the local precipitation and on the distinct infiltration conditions existing in the lithotypes.

The Soil and Water Assessment Tool (SWAT) is a physical model designed to predict the hydrological processes that could characterize natural and anthropized watersheds. The model can be forced using input data of climate prediction models, soil characteristics and land use scenarios to forecast their effect on hydrological processes. In this study, the SWAT model has been applied in the Aspio basin, a small watershed, highly anthropized and characterized by a short runoff generation. Three simulations setup, named SL1, SL2 and SL3, were investigated using different soil resolution to identify the best model performance. An increase of space requirement and calibration time has been registered in conjunction with the increasing soil resolution. Among all simulations, SL1 has been chosen as the best one in describing watershed streamflow, despite it was characterized by the lower soil resolution. A map of susceptibility to runoff for the entire basin was so created reclassifying the runoff amount of four years in five classes of susceptibility, from very low to very high. Eleven sub-basins, coinciding with the main urban settlements, were identified as highly susceptible to runoff generation. Considering future climate predictions, a slight increase of runoff has been forecasted during summer and autumn. The map of susceptibility successfully identified as highly prone to runoff those sub-basins where extreme flood events were yet recorded in the past, remarking the reliability of the proposed assessment and suggesting that this methodology could represent a useful tool in flood managing plan.

The design and construction of engineering works such as road infrastructures, port facilities and waste disposal plants require that important environmental issues be addressed, in order to guarantee sustainability and facilitate acceptability by the population. With these objectives, much of the scientific work produced by the geological and geotechnical research group of the Universita Politecnica delle Marche in the past 50 years has involved aspects such as understanding the effects of human activities on the environment in order to prevent or minimize the risks of permanent damage, exploration of solutions that optimize the use of natural resources, implementation of technologies that promote the use of improved natural materials for construction, application of design procedures that are based on high quality geological, geotechnical and geo-environmental models of the construction sites. This paper will present, via the illustration of some practical examples, the contribution that such research activities have provided on the solution of some important engineering problems in the past decades.

Several prior studies investigated the use of stable isotopes of water in hydrogeological applications, most on a local scale and often involving the isotopic gradient (evaluated by exploiting the so-called altitude effect), calculated on the basis of rainwater isotopes. A few times, this gradient has been obtained using the stable isotopic contents of low-yield springs in a limited time series. Despite the fact that this method has been recognized by the hydrogeological community, marked differences have been observed with respect to the mean stable isotopes content of groundwater and rainwater. The present investigation compares the stable isotopic signatures of 23 low-yield springs discharging along two transects from the Tyrrhenian sea to the Po Plain of Italy, evaluates the different isotopic gradients and assesses their distribution in relation to some climatic and topographic conditions. Stable isotopes of water show that groundwater in the study area is recharged by precipitation and that the precipitation regime in the eastern portion of the study area is strongly controlled by a shadow effect caused by the Alps chain on the air masses from central Europe. Stable isotopes (in particular the &delta, 18O and deuterium excess (d-excess) contents together with the obtained isotopic gradients) allow us to identify in the study area an opposite oriented orographic effect and a different provenance of the air masses. When the windward slope is located on the Tyrrhenian side, the precipitation shows a predominant oceanic origin, when the windward slope moves to the Adriatic side, the precipitation is characterized by a continental origin. The main results of this study confirm the usefulness of low-yield springs and the need for a highly detailed survey-scale hydrological investigation in the mountainous context.

Co-seismic changes in groundwater regime are often observed after moderate to strong earthquakes. The 24 August 2016 Mw 6.0 extensional Amatrice earthquake, which was the first event of a long-lasting seismic sequence, including the 30 October 2016 Mw 6.5 Norcia event, triggered a significant discharge alteration to the Pescara di Arquata spring, located in the Umbria-Marche Apennines (Northern Apennines, Central Italy) and exploited for drinking purposes. During the first five months after the first mainshock, an extra flow of about 30% was recorded, while both water chemistry and temperature did not show significant changes. Thereafter, the spring discharge decreased significantly, and at the end of 2019 it was still lower than normal. The Standardized Precipitation Index (SPI) indicates that these low mean monthly discharge values are not related to particularly dry conditions. The increase in post-seismic depletion coefficients indicates that the aquifer empties faster than it did during the inter-seismic period. The observed transient increase and subsequent decrease of discharge are consistent with a transient, earthquake-related increase in hydraulic conductivity.

A seismic sequence in central Italy from August 2016 to January 2017 affected groundwater dynamics in fractured carbonate aquifers. Changes in spring discharge, water-table position, and streamflow were recorded for several months following nine Mw 5.0–6.5 seismic events. Data from 22 measurement sites, located within 100 km of the epicentral zones, were analyzed. The intensity of the induced changes were correlated with seismic magnitude and distance to epicenters. The additional post-seismic discharge from rivers and springs was found to be higher than 9 m3/s, totaling more than 0.1 km3 of groundwater release over 6 months. This huge and unexpected contribution increased streamflow in narrow mountainous valleys to previously unmeasured peak values. Analogously to the L’Aquila 2009 post-earthquake phenomenon, these hydrogeological changes might reflect an increase of bulk hydraulic conductivity at the aquifer scale, which would increase hydraulic heads in the discharge zones and lower them in some recharge areas. The observed changes may also be partly due to other mechanisms, such as shaking and/or squeezing effects related to intense subsidence in the core of the affected area, where effects had maximum extent, or breaching of hydraulic barriers.

Tracing techniques are commonly used to investigate groundwater quality and dynamics, as well as to measure the hydrogeological parameters of aquifers. The last decade has seen a growing interest in environmentally friendly tracers, including single-stranded DNA molecules. In this study, an electrolytic tracer and a synthetic DNA tracer are comparatively evaluated in laboratory scale tests to assess their potential application in field studies aimed at investigating groundwater environments. A real-time quantitative Polymerase Chain Reaction assay was developed and optimized to detect and quantify the DNA tracer, while tracer column tests were performed to investigate the DNA tracer behavior and to compare it to the electrolytic tracer. The results show that the DNA tracer has an almost pure convective flow, while the KCl tracer experiences dispersive behavior. The tracing method proposed can be applied in hydrogeological field studies involving calcareous fractured rock systems, with the DNA tracer particularly suitable in tracing karst systems, which are often characterized by several conduits of flow. To test the DNA tracer in operation, a preliminary test was conducted in the field.

The investigated area, located in the inner part of the Marche region (central Italy) and belonging to the carbonate Umbria- Marche ridges in the central Apennines, is characterised by very complex geo-structural setting and widespread karst phenomena that make difficult the definition of the relation among the aquifers basing only on the hydrogeological survey. Hence, the presence of different flowpaths among aquifers of the Umbria-Marche hydrostratigraphic sequence and of tectonic contacts among the different structures is verified using tracer tests. In particular, the tests showed that the Calcare Massiccio and the Maiolica aquifers are connected under certain tectonic conditions. A new tracer given by a single stranded DNA molecule and traditional fluorescent dyes have been injected into the Montelago sinkhole in different periods (during the recharge and during the discharge) and recovered in several points along the expected hydrogeological basin, using either manual and automatic sampling. Fluorescent traps were positioned in creeks, rivers and springs. The DNA molecule is useful to trace surface water and groundwater, is detectable even at very low concentrations, no significant change in water density and viscosity can be observed and its use is not dangerous for the environment. The results stress the suitability of DNA as hydrogeological tracer, capable to identify connections among aquifers and study different flowpaths even in high flow conditions when traditional tracers are more and more diluted. Moreover, fluorescein tracer allowed for the transport parameter determination, giving mean velocities ranging from 100 to 3000 m/day and mean residence time from some tens to hundreds of hours, and determining the aquifer volumes.

The purpose of this paper is to highlight the hydrogeological contact among aquifers in karst and fissured systems and study the different flowpaths present in said aquifers; given the complex hydrogeological and tectonics settings, detailed surveys in the Umbria-Marche limestone ridges and tracer field tests were used to achieve this goal. DNA and fluorescent tracers were injected into a sinkhole and thereafter recovered at different points, allowing both for the identification of contact among aquifers, as well as the main and secondary directions of groundwater. The tracer tests have shown that the majority of the groundwater flow is directed towards the axis of the tectonic structures, emerging in more depressed areas. In addition, they have confirmed the influence of karstification on the recharge, which allows for a more detailed characterisation of water circulation in the unsaturated zone of the Maiolica and the Calcare Massiccio limestone aquifers thus offering confirmation of the hydrogeological survey results, and providing additional information that was previously unknown. Furthermore, the DNA tracer provided good results in the field especially in the karst systems and fissured rocks, thus confirming its appropriacy for flowpath investigations and vulnerability analyses of springs. In particular, DNA is suitable for very large amounts of groundwater owing to its very low detection limit. Further investigations and studies are, however, required to assess its validity also for the evaluation of hydrogeological parameters.

Discharge measurements in natural watercourses are performed in order to determine the value of the surface outflow of a basin, its temporal variability, and the outflow characteristics. The methods conventionally used for these measurements utilize an immersed current meter in different points of a river section, which acquires the mean flow velocity. Using this measurement, the discharge can be calculated. Some experimental problems arise, however, when there is a very high discharge. An important method, valid in such cases, is the artificial tracing method. In particular, the use of chemical tracers for small watercourses is very convenient because they are low cost, easily handled, low impact and provide satisfactory results. In the past, radioactive tracers such as tritium have been used in large rivers, while fluorescent tracers have been commonly exploited in the USA and now elsewhere. However, if the water is turbid, the suspended sediments may easily absorb some tracers. In this paper, the preliminary results of a comparison between current meter and artificial tracer measurements are reported. In particular, field tests in a small tributary have been performed, in order to investigate the behaviour of different tracers.

The town of Orvieto, located on the Rock of the same name, is an example of ''vulnerable town''; problems of slope instability connected with the lithological and morphological characteristics of the Rock have been thor- oughly examined and discussed during previous research studies. Hydrogeochemical data about groundwater recharging the springs present in the area were never taken into account. Pollution of the springs is well known but still occurs for unclear reasons. The aim of this work is there- fore to present the results of a hydrogeological and geo- chemical investigation of all the springs along the slopes of Orvieto hill and at the foot of the tuffaceous Rock, to characterize the groundwater flow paths and to suggest a possible source of contamination. The research study was carried out during three hydrogeochemical surveys in the years 1998-1999, 2003-2004, and 2007-2008.

This study examines the use of distributed power plants coupled to desalination systems to achieve combined production of electrical energy and of water for drinking and similar purposes in Marches, an Adriatic region in central Italy. To reduce Marches’ large electricity deficit, the regional government has drawn up an environmental and energy plan promoting distributed generation through power plants built right where the electricity is consumed. At the same time, the region’s abundant water resources are increasingly threatened by aquifer pollution, due particularly to fertilizers and to the phenomenon of saltwater intrusion near the coast. Distributed electrical energy generation combined with water production can help meet the goal of reducing the electrical deficit while at the same time improving ground-water quality and reducing the hydrogeological risk. The solution is also efficient, since the more the useful heat recovered by the energy system, the higher its global efficiency. Furthermore, cou...

The study of the hydrologic characters of a water course permits the correct management of the corresponding basin and a greater control over the water resources of the whole basin; therefore, a suitable planning and maintenance of the necessary interventions along the water course, especially in proximity of the outlet to sea, becomes necessary. An evaluation of the solid transport allows an estimation of the erosion to which the basin is subjected as a result of the river flow, and further helps to prevent hydrologic disasters in the possible risk zones. Among the various experimental techniques in use for measuring the suspended-solid transport, nuclear methods have been preferably used in this research, which are based on monitoring the concentration of the suspended sediments. The suspended-solid concentration is detected by the attenuation of radioactivity emitted by a source of 241Am dipped in the water. This attenuation, due to the presence of the sediments transported in great amounts during events of flood is measured using a scintillation detector made up of a crystal of NaI(Tl). With appropriate calibration curves built both in the laboratory and in the field, it is possible to trace the amount of suspended-solid transport in a certain river section that is located in the proximity of the river outlet. This methodology, applied to different equipped basins in Italy and Africa, is particularly useful for small and medium water courses (similar to those of the Apennine ranges in Italy), allowing an assessment of the erosion in the whole watershed. In this note, the techniques used are introduced in detail, with particular attention to the instrument calibration, and the numerical results obtained for some basins in the Marche region (Italy) are compared with some empirical formulae used in previous reports for the calculation of erosion.

12 springs from the northern Apennines of Italy were studied by means of a comprehensive hydrogeological investigation to unravel recharge processes taking place in a highly fractured slab of flysch rock hosting the corresponding aquifers. Several campaigns were carried out during the period 2005–2008 to gather springs’ discharge together with electrical conductivity and temperature data. Water samples were collected and allowed the determination of the major ions (Ca2+, Mg2+, Na+, K+, SO4 2−, HCO3 −, Cl−) and the isotopic content (δ18O). Additionally, groundwater discharge from 3 selected springs was measured continuously using pressure-transducers. Over the same period, reference δ18O values for rainfall and snow-melt water were collected monthly by a rain-gauge. Results highlighted an aquifer-like behavior of this slab of flysch rock. The springs are characterized by a Ca–HCO3 hydro-facies and an increase of groundwater mineralization has been noticed moving toward the base of the slab; the mean discharges are between 0.3 and 0.7 l·s−1 and suffered a marked variability during the year; the total maximum yield is about 155,000 m3 while the total discharge volume is around 320,000 m3·y−1. By processing the δ18O isotope values from springs (mean annual values ranged from −9.67 and −10.42 ‰) and by combining them with rainfalls datasets, results show that the main aliquot of recharge occurs in the winter—spring months and it is principally related to the solid precipitations (snow-melt).

The purpose of this paper is to study floods, groundwater pollution and landslides in a small watershed characterized by high hydrogeological hazard, high population density and widespread industrial and commercial settlements. Measurements were performed to study hydrogeological hazard in the watershed during 2011, 2012 and 2013. First the relation between groundwater and surface water was estimated; second, the rainfall-runoff relationship was assessed. Some rainfall events were analysed to investigate the watershed behaviour in different periods over the hydrologic year. The analysis stressed the importance of the rainfall events, soil moisture and hydrological conditions. From this study it’s apparent that the natural hazard (i.e. floods, landslides and pollution hazard) can be better individuated and predicted in a small basin rather than in a wider one. Furthermore, these studies can be easily implemented in other basins with similar features in all the northern Mediterranean area.