On the 4th and 5th of March 2005, about 100 rainfall-induced landslides occurred along volcanic slopes of Camaldoli Hill in Naples, Italy. These started as soil slips in the upper substratum of incoherent and welded volcaniclastic deposits, then evolved downslope according to debris avalanche and debris flow mechanisms. This specific case of slope instability on complex volcaniclastic deposits remains poorly characterized and understood, although similar shallow landsliding phenomena have largely been studied in other peri-volcanic areas of the Campania region underlain by carbonate bedrock. Considering the landslide hazard in this urbanized area, this study focused on quantitatively advancing the understanding of the predisposing factors and hydrological conditions contributing to the initial landslide triggering. Borehole drilling, trial pits, dynamic penetrometer tests, topographic surveys, and infiltration tests were conducted on a slope sector of Camaldoli Hill to develop a geological framework model. Undisturbed soil samples were collected for laboratory testing to further characterize hydraulic and geotechnical properties of the soil units identified. In situ soil pressure head monitoring probes were also installed. A numerical model of two-dimensional variably saturated subsurface water flow was parameterized for the monitored hillslope using field and laboratory data. Based on the observed soil pressure head dynamics, the model was calibrated by adjusting the evapotranspiration parameters. This physically based hydrologic model was combined with an infinite-slope stability analysis to reconstruct the critical unsaturated/saturated conditions leading to slope failure. This coupled hydromechanical numerical model was then used to determine intensity–duration (I-D) thresholds for landslide initiation over a range of plausible rainfall intensities and topographic slope angles for the region. The proposed approach can be conceived as a practicable method for defining a warning criterion in urbanized areas threatened by rainfall-induced shallow landslides, given the unavailability of a consistent inventory of past landslide events that prevents a rigorous empirical analysis.

In many Italian regions, and particularly in southern Italy, karst aquifers are the main sources of drinking water and play a crucial role in the socio-economic development of the territory. Hence, estimating the groundwater recharge of these aquifers is a fundamental task for the proper management of water resources, while also considering the impacts of climate changes. In the southern Apennines, the assessment of hydrological parameters that is needed for the estimation of groundwater recharge is a challenging issue, especially for the spatial and temporal inhomogeneity of networks of rain and air temperature stations, as well as the variable geomorphological features and land use across mountainous karst areas. In such a framework, the integration of terrestrial and remotely sensed data is a promising approach to limit these uncertainties. In this research, estimations of actual evapotranspiration and groundwater recharge using remotely sensed data gathered by the Moderate Resolution Imaging Spectrometer (MODIS) satellite in the period 2000–2014 are shown for karst aquifers of the southern Apennines. To assess the uncertainties affecting conventional methods based on empirical formulas, the values estimated by the MODIS dataset were compared with those calculated by Coutagne, Turc, and Thornthwaite classical empirical formulas, which were based on the recordings of meteorological stations. The annual rainfall time series of 266 rain gauges and 150 air temperature stations, recorded using meteorological networks managed by public agencies in the period 2000–2014, were considered for reconstructing the regional distributed models of actual evapotranspiration (AET) and groundwater recharge. Considering the MODIS AET, the mean annual groundwater recharge for karst aquifers was estimated to be about 448 mm·year⁻¹. In contrast, using the Turc, Coutagne, and Thornthwaite methods, it was estimated as being 494, 533, and 437 mm·year⁻¹, respectively. The obtained results open a new methodological perspective for the assessment of the groundwater recharge of karst aquifers at the regional and mean annual scales, allowing for limiting uncertainties and taking into account a spatial resolution greater than that of the existing meteorological networks. Among the most relevant results obtained via the comparison of classical approaches used for estimating evapotranspiration is the good matching of the actual evapotranspiration estimated using MODIS data with the potential evapotranspiration estimated using the Thornthwaite formula. This result was considered linked to the availability of soil moisture for the evapotranspiration demand due to the relevant precipitation in the area, the general occurrence of soils covering karst aquifers, and the dense vegetation.

Shallow aquifers of coastal and internal alluvial plains of developed countries are commonly characterized by the challenging management of groundwater resources due to the intense agricultural and industrial activities that determine a high risk of groundwater contamination. Among the principal origins of pollution in these areas are agricultural practices based on the amendment of soils by nitrate fertilizers, which have been recognized as one of the most severe environmental emergencies for which specific policies and regulations have been issued (e.g., EU Directive 2006/118/EC). In such a framework, the results of research aimed at assessing the specific vulnerability of shallow alluvial aquifers to nitrate fertilizer pollutants by coupled process-based and empirical approaches are here proposed. The research focused on assessing the specific vulnerability to nitrate pollution of a shallow alluvial aquifer of the Campania region (southern Italy), which was selected due to its representativeness to other recurrent hydrogeological settings occurring in alluvial plains of the region and worldwide. In this area, 1D hydro-stratigraphic models of the unsaturated zone were reconstructed and applied for simulating the transport of nitrate pollutants at the water table and estimating the associated travel times. Numerical modeling was carried out by the finite differences VS2TDI code and considered a 10-year time series of rainfall and evapotranspiration as well as typical local farming practices of nitrate fertilizer input. Results of the travel time calculated for the 1D hydro-stratigraphic models considered and at different depths were recognized as a proxy to assess the specific vulnerability to nitrate fertilizer pollution. Among the principal outcomes is an empirical multiple correlation between the travel time of the nitrate fertilizer pollutant, water table depth, and equivalent saturated hydraulic conductivity of the unsaturated zone or hydraulic resistance, which was used to assess the travel time at the distributed scale over the whole area studied as well as the related specific vulnerability. Given such results, the coupled process-based and empirical approach is proposed as generally applicable for assessing and mapping groundwater vulnerability in shallow aquifers, for which detailed stratigraphic and piezometric data are available.