Using CPTu profiles for subsoil characterisation, transformation equations must be used to obtain the hydro-mechanical properties for structures and infrastructure designing. Additionally, the uncertainty and the spatial variability of measured parameters must be taken into account for a reliable geotechnical design. In this work, we used a Stochastic Simulation approach to define reliable 3D models of two geotechnical designing variables for granular soils (friction angle–ϕ’ and the Darcy permeability coefficient–k) from tip resistance (qc), sleeve friction (fs), and pore pressure (u2) profiles. The selected method – the Sequential Gaussian Co-Simulation (SGCS) – provided reliable optimized 3D models of the spatial distribution of the variables of interest and allowed quantifying the propagation of the estimation uncertainty associated with the raw measurement models through the transformation equations. Overestimation (OE) and Underestimation (UE) percentages for a confidence interval of 68% were calculated throughout the 3D model: granular soils showed a larger uncertainty than fine soils concerning the measured variables (qc, fs, and u2). In granular soils, the measured variable uncertainty varies up to 100% but the derived variables show different behavior: ϕ’ shows UE and OE less than 25% while k reaches 100%. These differences in the propagated uncertainties depend on the transformation equations and the measured variable dependence. @en

Biomarkers such as functional gene mRNA (transcripts) and proteins (enzymes) provide direct proof of metabolic regulation during the reductive dechlorination (RD) of chlorinated ethenes (CEs). Yet, current models to simulate their spatiotemporal variability are not flexible enough to mimic the homologous behavior of RDase functional genes. To this end, we developed new enzyme-based kinetics to model the concentrations of CEs together with the transcript and enzyme levels during RD. First, the model was calibrated to existing microcosm data on RD of cis-DCE. The model mirrored the tceA and vcrA gene expression and the production of their enzymes in Dehalococcoides spp. Considering tceA and vcrA as homologous instead of nonhomologous improved fitting of the mRNA time series. Second, CEs and biomarker patterns were explored as a proof of concept under groundwater flow conditions, considering degraders occurring in immobile and mobile states. Under both microcosm and flow conditions, biomarker-rate relationships were nonlinear hysteretic because tceA and vcrA acted as homologous genes. The mobile biomarkers additionally undergo advective-dispersive transport, which increases the nonlinearity and makes the observed patterns even more challenging to interpret. The model offers a thorough mechanistic description of RD while also allowing simulation of spatiotemporal dynamic patterns of various key biomarkers in aquifers.@en

In the past, soil-layer delineation methods can usually only take a single type of input data, e.g., soil-type data at boreholes. However, this does not fit in the geotechnical engineering practice where multiple types of data are usually available during site investigation (e.g., borehole data and cone penetration test data are both available). This paper adopts a novel data-driven method for soil-layer delineation that accommodates multiple types of site investigation data. The basic idea is to include liquid limit (LL), plasticity index (PI), and fines content (FC) into the soil parameters of analysis. According to the Unified Soil Classification System (USCS), the information of (LL, PI, FC) can be used to determine whether the soil is sand, silt, or clay. As a result, the conditional random field simulation results for (LL, PI, FC) can be used to delineate sand, silt, and clay layers. If extra soil parameters (such as cone penetration test results) are incorporated, the novel method can accommodate multiple types of site investigation data. A real example of the Fucino Basin in Italy is adopted to demonstrate the application of the novel data-driven soil-delineation method.@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 2016–2017 seismic events that struck central Italy led the Government to carry out a project to produce the third level Seismic Microzonation studies in 138 municipalities. These activities have involved many experts in different disciplines such as geology, geomorphology, geophysics, seismology and geotechnical engineering. This project represented the first opportunity to perform nationally coordinated third level Seismic Microzonation studies over a wide area in a quite short time (6 months). It provided the chance to improve methodological procedures, to test the reliability of methods and models for site response analyses and to produce a huge amount of validated data. This paper focuses on the contribution of geological disciplines and concerns: (a) the definition of the main “morphostructural domains” of the Central-Northern Apennines; (b) the creation of an archive of all the lithostratigraphic units occurring in the study area with their conversion into engineering-geological units and their distribution in the different morphostructural domains; (c) the construction of the reference geological and geotechnical models, which are essential to classify the territory into seismically homogeneous microzones and to perform the successive 1D and 2D numerical analyses of the local site response. The geophysical dataset acquired for the study allowed a first statistical characterization of the Vs values typical of the engineering-geological units identified in this study. Some examples of the recurrent geological and geotechnical models are shown to explain the complexity and variety of the geological and geomorphological features of the investigated area and to highlight the different seismostratigraphic behavior of rocks and cover terrains. The analysis of third level Seismic Microzonation data made it possible to identify recurrent subsoil models and to note the main stratigraphic and morphological control-factors of the ground motion modification in the different morphostructural domains.

@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 Cone Penetration Test (CPT) measures enable to recognize vertical lithological sequence at each investigated point. From the tip resistance q_c and sleeve resistance f_s profiles, the Soil Behavior Type index I_SBT has been calculated, in order to identify the lithotypes alongside depth. The present study focuses on the combination of different variables to provide a lithological and mechanical subsoil characterization. The main objectives of the paper are: (1) to model the 3D spatial variability structure of the soil lithotypes and mechanical properties using q_c and f_s profiles; 2) to evaluate the uncertainties of the estimates for designing purposes. 182 CPTs were collected in a 900 km² area (corresponding to a subsoil volume of about 12 km³) located in the study site in the Bologna province (Italy). The study area is made up of fine-grained soils, silt and clay mixtures that are intercalated at different depths by sandy and gravelly soils. These variations of each soil fraction affect the engineering properties of these alluvial deposits. For 3D modeling, two geostatistical methods, Ordinary Kriging (stationary method), and Intrinsic Random Function theory (non-stationary method) have been used. The results show that the non-stationary method allows to obtain more reliable q_c and f_s predicted values. The final stochastic mechanical and lithological model enables engineers and geologists to detect the emerging of fan and paleochannels bodies where mean resistance values can abruptly change in terms of bearing capacity, liquefaction potential and static and dynamic settlement occurrence.

@en

The definition of hydrogeological conceptual models sometimes suffers on uncertainties depending on the available dataset. The availability and the elaboration of long-term hydrogeological dataset allow ensuring a better definition of the aquifer dynamics and features. One of the main features, which can benefice on long-term datasets, is the recharge quantification. According to the literature, several methods exist to calculate the recharge and to get a deeper insight in recharge mode assessment. Most of them allow the quantification of infiltrating water by means of meteo-climatic data and land use information. In this research, a combined approach on hydrogeological time-series is proposed. The Cross-Correlation among the rainfall and the piezometric levels was applied to point out the response time of each aquifer layer. In addition, the Water Table Fluctuation method was used to indirectly quantify the recharge. The study area is the Roccamonfina Volcano caldera (Italy), where lateral inflows are absent, and the almost total natural groundwater discharge occurs in the Savone delle Ferriere stream. Here, a monitoring scheme of the major hydrogeological features was implemented. The local rainfall, the groundwater levels in 2 wells tapping the different aquifer layers, and the river stage of the Savone delle Ferriere stream were measured daily during the 2016–2017 period. The analyses defined different recharge modes for each aquifer layer. The aquifer recharge rates, calculated with different approaches, are comparable to each other and in line with the mean groundwater discharge through the Savone delle Ferriere streambed springs.

@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

This research deals with some hydrogeochemical and hydrodynamic features that affect redox processes in groundwater, presenting three different case studies and methodological approaches. This infor-mation is vital when characterizing contaminated sites, planning monitoring, selecting remedial techniques, and assessing health risks. In the San Pedro Sula site case study (Honduras), new insights and an advanced method for the aquifer redox zonation were provided. The Multi-Collocated Factorial Kriging, a multivariate geostatistical technique borrowed from agriculture and soil science, identified simultaneous redox processes, acting at different scales and mainly due to eutrophicated surface-water/groundwater interaction: at short-range, Mn oxi-hydroxide reductive dissolution and Fe oxi-hydroxide colloidal phase precipitation; at long-range, Fe mobilization. The obtained results are supported by the Principal Component Analysis and hydrogeochemical numerical modeling. The interaction among different contamination events in the VR site allowed to get a deeper insight into the competition for organic substrate among different redox processes. To this purpose, a 1-D Advective-Reactive-Dispersive transport numerical model was implemented, and its results highlighted that the reductive dechlorination is inhibited by the inorganic Terminal Electron Accepting Processes, enriching groundwater in highly toxic by-products of the reductive dechlorination (i.e. Vinyl Chloride). As a matter of fact, the bacteria using inorganic Terminal Electron Acceptors to oxidize the organic substrate are more efficient than the halo-respiring bacteria, that favor reductive dechlorination. The contaminated site of Bussi Officine case study allowed to clarify the hydrodynamic effect on the redox processes and the anoxic condition distribution. The comparison between dissolved chlorinated solvents and the inorganic Terminal Electron Acceptors distribution in the aquifer pointed out that dispersion inhibits all the redox processes because it dilutes the organic substrate concentration in groundwater. On the other hand, the lacustrine deposits rich organic matter (i.e. peat) can be considered as "chemical reactors", that releases by back-diffusion toxic by-product of the reductive dichlorination in groundwater.

@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

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 m³/s, totaling more than 0.1 km³ 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.

@en