Groundwater is among the most important freshwater resources. Worldwide, aquifers are experiencing an increasing threat of pollution from urbanization, industrial development, agricultural activities and mining enterprise. Thus, practical actions, strategies and solutions to protect groundwater from these anthropogenic sources are widely required. The most efficient tool, which helps supporting land use planning, while protecting groundwater from contamination, is represented by groundwater vulnerability assessment. Over the years, several methods assessing groundwater vulnerability have been developed: overlay and index methods, statistical and process-based methods. All methods are means to synthesize complex hydrogeological information into a unique document, which is a groundwater vulnerability map, useable by planners, decision and policy makers, geoscientists and the public. Although it is not possible to identify an approach which could be the best one for all situations, the final product should always be scientific defensible, meaningful and reliable. Nevertheless, various methods may produce very different results at any given site. Thus, reasons for similarities and differences need to be deeply investigated. This study demonstrates the reliability and flexibility of a spatial statistical method to assess groundwater vulnerability to contamination at a regional scale. The Lombardy Plain case study is particularly interesting for its long history of groundwater monitoring (quality and quantity), availability of hydrogeological data, and combined presence of various anthropogenic sources of contamination. Recent updates of the regional water protection plan have raised the necessity of realizing more flexible, reliable and accurate groundwater vulnerability maps. A comparison of groundwater vulnerability maps obtained through different approaches and developed in a time span of several years has demonstrated the relevance of the continuous scientific progress, recognizing strengths and weaknesses of each research.

Contamination by chlorinated solvents is typically associated with point sources, which are able to release high concentrations and to generate well defined plumes. Nevertheless, in urban settings (especially in functional urban areas—FUAs), multiple-point sources are frequently present, consisting of a series of unidentifiable small sources clustered within large areas, generating a diffuse, anthropogenic contamination. This situation results in the coexistence of single plumes with higher contaminant concentrations, and larger areas where the concentration is lower but still higher than the maximum admissible concentration limits. This paper proposes a methodology devised to cope with the diffuse contamination by chlorinated solvents within shallow aquifers due to multiple-point sources in FUAs. The approach is based on a Bayesian model that helps to spatially evaluate the likelihood of having active multiple-point sources, and to relate their impact on the shallow aquifer to the hydrogeological features of the area. Moreover, the approach allows testing of the efficiency of the monitoring network to properly characterize the contamination in the aquifer. The consistency of the results of the analysis was also checked for the Milan FUA (Italy) by a comparison to a previous study, performed through an inverse numerical modelling approach within a Monte Carlo statistical framework to identify the areas with the highest likelihood to host potential multiple-point sources.

Groundwater resources in mountain areas are strategically important to maintain adequate water supply for domestic uses, farming, industrial activities, and energy production, also considering the expected growing demand due to ongoing climate changes. Within this framework, the objective of the study is to develop a regional approach, compliant with the European requirements of the Water Framework Directive 2000/60/EC and Groundwater Directive 2006/118/EC, that could support public agencies and water companies to efficiently manage and protect the available water resources in mountainous environments. The proposed approach identifies and delineates groundwater bodies by coupling a 3D hydro-stratigraphic model with the definition of the water budget and water hydrochemical fingerprints in a geologically complex Alpine environment in Northern Italy. Sixteen groundwater bodies (GWBs) have been identified all over the 10.290 km² area, showing an average storage capacity of more than 500 Mm³ y^-1 (about 3% of the average total inflow from precipitation and snowmelt), with differences up to four times between GWBs mainly constituted of carbonate rocks and those prevalently composed of crystalline or terrigenous rocks. Groundwater quality in the study domain is generally excellent, with few exceptions due to geogenic (i.e., natural) or anthropogenic sources of contamination. The results of this study show the advantages of coupling 3D hydro-stratigraphic modelling combined with meteorological, hydrological and hydrogeological information, which consist in: i) identifying the most Strategic Storage Reservoir both in terms of quality and storage capacity; ii) evaluating the present ground- and surface water availability; iii) detecting areas of specific interest for implementing groundwater monitoring networks; iv) recognising recharge areas of the most relevant springs, to implement protection strategies of the resource.