Riccardo Pinardi

Few hydrogeological studies have focused on possible per- and poly-fluoroalkyl substance (PFAS) contamination in groundwater with particular attention to the role of hydraulic interconnections and to the interdigitations present between shallow and deep aquifer layers in heterogeneous alluvial systems. In general, deeper groundwater is considered chemically safer and less impacted by contamination, especially in multilayer aquifers characterized by low permeability apparently confining horizons. Therefore, this research analyzed PFAS in groundwater at depths ranging from 20 to 120 m below ground level, combining stratigraphic, hydrogeological, and chemical data with GIS mapping to identify industrial activities potentially contributing to PFAS contamination using the cross-checking methodology. During the second survey, the monitoring network was extended along a hydrogeological transect, including two springs located upstream and downstream of the deep wells, to assess PFAS concentration in shallow groundwater and the possible transfer along the groundwater flow path. The intra-site comparative analysis reveals, for the same sampling locations, a differentiation in the PFAS profiles detected across the two monitoring campaigns, indicating a temporal evolution in the chemical composition. Furthermore, chemical results show the presence of PFAS exclusively in deep monitoring wells, confirming a spatially heterogeneous distribution within the aquifer system. These results highlight both the temporal and spatial evolution of PFAS concentration, suggesting a complex contaminant migration pathway along preferential gravel and sand horizons in deeper aquifer layers. The conceptual hydrogeological model confirmed hydraulic interconnections among aquifer layers and identified zones of higher vulnerability to contamination. The analysis of possible PFAS migration pathways at the basin scale raised some questions about the influence of wells features and management practices on PFAS distribution in shallow and deep groundwater. The findings of this research contribute to environmental sustainability, providing initial insights for measuring and managing the presence and pathways of PFAS in deep alluvial aquifers.

Chlorinated organic compounds are Persistent Organic Pollutants (POPs) with high environmental continuity. These chemicals possess the ability to permeate into the environment across both unsaturated and saturated zones. This study examines the potential impact of perchloroethylene (PCE) releases in aquifer systems consisting of layers with varying permeability. The numerical simulations utilized the CactusHydro numerical code that employs a high-resolution shock-capturing flux conservative method to solve the non-linear partial differential equations of a three-phase immiscible fluid flow and study the migration of PCE into variably saturated zones to predict the spatial and temporal distribution of free products across multilayered aquifers. Two scenarios were simulated to predict how a first low-permeable horizon would affect the downward migration of PCE in its pure phase. The numerical simulations show that the multilayered aquifer system is vulnerable to an in-depth study of the PCE migration paths in the pure phase. However, very long times (about a year) are required for the DNAPL to be able to cross the semi-permeable horizons. The results in both scenarios show the quantity of mass (in kg and percentage concerning the initial one) that is left along the multilayered aquifer during the migration and after several years.

The significance of the interconnection between water and energy, known as the water–energy (WE) nexus, is highly regarded in scientific publications. This study used a narrative review method to analyze the existing WE nexus studies performed before 2024 in 26 European countries. The aim of this study is to provide a comprehensive analysis of the existing WE nexus to identify research gaps and to report a conceptual overview of energy consumption related to groundwater use phases, ranging from the tapping to distribution. This information is valuable as a guideline for any future estimates in this field. The results indicate that the WE nexus in 26 European countries comprises a variety of topics, including the water supply system, wastewater treatment, hydropower, desalination, and biofuel production. Most of the focus has been on fossil fuel production, while water supply and desalination were considered rarely. Italy and Portugal had the largest WE nexus. It is highlighted that there have been no studies on the WE nexus focusing on the groundwater supply system that consider the conceptual hydrological model or hydrodynamic processes. In this work, a view of these aspects was provided by taking into account different hydrogeological and hydraulic scenarios that may affect the amount of energy required for groundwater exploitation. Most scientific publications have focused on quantitative analysis. In the future, it will be necessary for WE nexus models to place a greater emphasis on governance and the implications of the WE nexus approach.

In recent years, there has been a surge in interest concerning emerging contaminants, also known as contaminants of emerging concern (CECs), due to their presence in environmental matrices. Despite lacking regulation, these chemicals pose potential health and environmental safety risks. Disodium EDTA, a widely utilized chelating agent, has raised concerns regarding its environmental impact. The present work aimed to verify the presence of Disodium EDTA at the exit of eight wastewater treatment plants discharging into some losing streams flowing within a large alluvial aquifer. Conducted in the Province of Parma (Northern Italy), the research employs a multidisciplinary approach, incorporating geological, hydrogeological, chemical, and microbial community analyses. Following a territorial analysis to assess industries in the region, through the use of ATECO codes (a classification system for economic activities), the study investigated the concentration of Disodium EDTA in effluents from eight diverse wastewater treatment plants, noting that all discharges originate from an activated sludge treatment plant, released into surface water courses feeding the alluvial aquifer. Results revealed detectable levels of Disodium EDTA in all samples, indicating its persistence post-treatment. Concentrations ranged from 80 to 980 µg/L, highlighting the need for further research on its environmental fate and potential mitigation strategies. Additionally, the microbial communities naturally occurring in shallow groundwater were analyzed from a hydrogeological perspective. The widespread presence of a bacterial community predominantly composed of aerobic bacteria further confirmed that the studied aquifer is diffusely unconfined or semi-confined and/or diffusely fed by surface water sources. Furthermore, the presence of fecal bacteria served as a marker of diffuse leakage from sewage networks, which contain pre-treated wastewater. Although concentrations of Disodium EDTA above the instrumental quantification limit have not been found in groundwater to date, this research highlights the significant vulnerability of aquifers to Disodium EDTA. It reveals the critical link between surface waters, which receive treated wastewaters impacted by Disodium EDTA, and groundwater, emphasizing how this connection can expose aquifers to potential contamination. At this stage of the research, dilution of wastewaters in surface- and groundwater, as well as hydrodynamic dispersion within the alluvial aquifer, seem to be the main factors influencing the decrease in Disodium EDTA concentration in the subsurface below the actual quantification limit. Consequently, there is a pressing need to enhance methodologies to lower the instrumental quantification limit within aqueous matrices. In a broader context, urgent measures are needed to address the risk of diffuse transport of CECs contaminants like Disodium EDTA and safeguard the integrity of surface and groundwater resources, which are essential for sustaining ecosystems and human health.

The occurrence of emerging pollutants (EPs) such as pharmaceuticals and personal care products (PPCPs) has raised serious concerns about the possible adverse effects on ecosystem integrity and human health. Wastewater treatment facilities appear to be the main sources of PPCPs released in aquatic environments. This research examines the effectiveness of groundwater microbial community activities to remove phenoxyethanol (Phy-Et), currently exploited as a preservative in many cosmetic formulations at a maximum concentration of 1% but which has shown, at higher levels of exposure, adverse systemic effects on animals. Mesocosm experiments were carried out for 28 days using two different concentrations of the substance (5.2 mg/L and 27.4 mg/L). The main results obtained through chemical and microbiological investigations revealed a significant Phy-Et reduction (≈100% when added at a concentration of 5.2 mg/L and ≈84% when added at a concentration of 27.4 mg/L), demonstrating that some autochthonous microorganisms in the analyzed samples played a “key role” in removing this compound, despite its proven antimicrobial activity. Nevertheless, the decrease in the “natural attenuation” efficacy (≈16%) when using higher concentrations of the chemical suggests the existence of a “dose-dependent effect” of Phy-Et on the process of biodegradation. Biomolecular investigations carried out through next-generation sequencing (NGS) revealed (i) the presence of a significant fraction of hidden microbial diversity to unravel, (ii) variations of the composition and species abundance of the groundwater microbial communities induced by Phy-Et, and (iii) a biodiversity reduction trend correlated to the increase of Phy-Et concentrations. Overall, the preliminary information obtained from the experiments carried out at the laboratory scale appears encouraging, although it reflects only partially the complexity of the phenomena that occur in natural environments and influences their “auto-purification capability”. Accordingly, this research paves the way for more in-depth investigations to develop appropriate tools and protocols to evaluate the occurrence and fate of Phy-Et in nature and assess the impact of its release and the effects of long-term exposure (even at low concentrations) on ecosystems and health.

Dense nonaqueous phase liquids (DNAPLs) are known to be denser than water and immiscible with other fluids. Once released into the environment, they migrate downward through the variably saturated zone, causing severe damage. For this reason, it is essential to properly develop a rapid response strategy, including predictions of contaminant migration trajectories from numerical simulations modeling. This paper presents a series of simulations of free-product DNAPL extraction by means of a purpose-designed pumping well. The objective is to minimize the environmental impact caused by DNAPL release in the subsurface, estimating the recoverable free-product DNAPL, depending on the hydraulic properties of the aquifer medium, and estimating the leaving residual DNAPL that could act as a long-term pollution source. Coupling the numerical simulations to the bacterial community characterization (through biomolecular analyses), it was verified that (i) the DNAPL recovery (mainly PCE at the study site) through a pumping well would be almost complete and (ii) the application of other remediation techniques (such as bioremediation) would not be necessary to remove the pollution source because (iii) a natural attenuation process is provided by the autochthonous bacterial community, which is characterized by genera (such as Dechloromonas, Rhodoferax, and Desulfurivibrio) that have metabolic pathways capable of favoring the degradation of chlorinated compounds.

Chlorinated organic compounds are widespread aquifer contaminants. They are known to be dense non-aqueous phase liquids (DNAPLs). Therefore, they are denser than water and immiscible with other fluids. Their migration into the environment in variably saturated zones can cause severe damage. For this reason, optimizing those actions that minimize the negative impact of these compounds in the subsurface is essential. This paper presented a numerical model simulating the free-product DNAPL migration and extraction through a purpose-designed pumping well in a potential emergency scenario. The numerical simulations were performed using CactusHydro, a numerical code that uses a high-resolution shock-capturing flux conservative method to resolve the non-linear coupled partial differential equations of a three-phase immiscible fluid flow recently proposed in the literature, including the contaminant extraction at the base of the aquifer. We investigated the temporal (and spatial) evolution of its migration in the Parma (Northern Italy) porous alluvial aquifer following the saturation contour profiles of the three-phase fluid flow in variably saturated zones. The results indicated that this numerical approach can simulate the contaminant migration in the subsurface and the pumping of the free-product from a well screened at the base of the aquifer system. Moreover, the simulation showed the possibility of recovering about two-thirds of the free-product, in agreement with the scientific literature.

Hydrogeological maps must synthesize scientific knowledge about the hydraulic features and the hydrogeological behavior of a specific area, and, at the same time, they must meet the expectations of land planners and administrators. Thus, hydrogeological maps can be fully effective when they are purpose-designed, especially in complex interconnected systems. In this case study, purpose-designed graphical solutions emphasize all the hydraulic interconnections that play significant roles in recharging the multilayered alluvial aquifer, where the majority of wells have been drilled for human purposes, artificial channels are used for agricultural purposes, and the shallow groundwater feeds protected groundwater-dependent ecosystems. The hydrogeological map was then designed to be the synthesis of three different and hydraulically interconnected main contexts: (i) the alluvial aquifer, (ii) the hydrographic basin of the Taro losing river, and (iii) those hard-rock aquifers whose springs feed the same river. The main hydrogeological map was integrated with two smaller sketches and one hydrogeological profile. One small map was drawn from a modeling perspective because it facilitates visualization of the alluvial aquifer bottom and the “no-flow boundaries.” The other small sketch shows the artificial channel network that emphasizes the hydraulic connection between water courses and groundwater within the alluvial aquifer. The hydrogeological profile was reconstructed to be able to (i) show the main heterogeneities within the aquifer system (both layered and discontinuous), (ii) visualize the coexistence of shallower and deeper groundwater, (iii) emphasize the hydraulic interconnections between subsystems, and (iv) suggest the coexistence of groundwater pathways with different mean residence times.

The contamination impact and the migration of the contaminant into the surrounding environment due to the presence of a spilled oil pipeline will cause significant damage to the natural ecosystem. For this reason, developing a rapid response strategy that might include accurate predictions of oil migration trajectories from numerical simulation modeling is decisive. This paper uses a three-dimensional model based on a high-resolution shock-capturing conservative method to resolve the nonlinear governing partial differential equations of the migration of a spilled light nonaqueous liquid oil contaminant in a variably saturated zone employed to investigate the migration of the oil pipeline leakage with great accuracy. The effects of the oil type density, gasoline, and diesel oil, the unsaturated zone depth, its saturation, the hydraulic gradient, and the pressure oil pipeline are investigated through the temporal evolution of the contaminant migration following the saturation profiles of the three-phase fluid flow in the variably saturated zone. The calculation results indicate that the leaking oil’s pressure is the parameter that significantly affects the contaminants’ arrival time at the groundwater table. Additionally, the water saturation of the unsaturated zone influences the arrival time, as the water saturation increases at a fixed depth. The unsaturated zone depth significantly influences the contaminant migration in the unsaturated zone. At the same time, the oil density and the hydraulic gradient have limited effects on the contaminant migration in the variably saturated zone.

Hydrogeological maps must synthesize the scientific knowledge about the hydraulic features and the hydrogeological behavior of a specific area, and, at the same time, they must meet the expectations of land planners and administrators. Thus, especially in complex interconnected systems, hydrogeological maps can be fully effective when they are purpose-designed. In these contexts, usual graphical approaches in creating hydrogeological maps could result uncomplete and/or ineffective. Therefore, in case of complex systems and/or specific management/protection aims, new and purpose-designed graphical solutions must be applied to enhance the maps effectiveness. In the case study, these solutions show and emphasize all the hydraulic interconnections playing significant roles in recharging the multilayered alluvial aquifer, where the majority of wells have been drilled for drinking/industrial/agricultural purposes, artificial channels are used for agricultural purposes, and the shallow groundwater feeds protected groundwater-dependent-ecosystems. The hydrogeological map was then designed to be the synthesis of three different and hydraulically interconnected main contexts: (i) the main heterogeneous alluvial aquifer (the main target of the purpose-designed map), (ii) the hydrographic basin of a losing river that feeds the main alluvial aquifer, and (iii) those hard-rock aquifers (mainly turbiditic and ophiolitic) whose springs feed the same river.