Using index and physically-based models to evaluate the intrinsic groundwater vulnerability to non-point source pollutants in an agricultural area in Sardinia (Italy)

This research aims at studying the intrinsic vulnerability of groundwater to diffuse environmental pollutants in the Muravera coastal agricultural area of Sardinia, Italy. The area faces contamination risks arising from agricultural practices, especially the use of fertilizers, pesticides, and various chemicals that can seep into the groundwater. The study examined the interplay among hydrological elements, including soil characteristics, groundwater depth, climate conditions, land use, and aquifer properties. To do that, the outcomes of FLOWS 1D physically-based agrohydrological model were analyzed in parallel with those of the overlay-and-index model SINTACS, in a sort of reciprocal benchmarking. By using FLOWS, water movement and solute transport in the unsaturated zone were simulated by, respectively, solving the Richard Equation (RE) and the Advection-Dispersion equation (ADE). As such, this model allowed to account for the role of soil hydraulic and hydro-dispersive properties variability in determining the travel times of a conservative solute through the soil profile to the groundwater. For FLOWS simulations, a complete dataset was used as input, including soil horizons, soil physical and hydraulic properties of 36 soil profiles, average annual depth to groundwater table at each soil profile (ranging from 1 to 50 meters), and climatic temporal series data on rainfall and evapotranspiration. Detailed analyses of travel times for the movement of 25, 50, 75, and 100% of the solute mass to reach groundwater were conducted, revealing that the depth to groundwater predominantly influences vulnerability. This result was coherent with SINTACS vulnerability map due to the large impact of the depth to groundwater on SINTACS analysis.