Potential evapotranspiration (ET0) is an indicator of great interest for water budget analysis and the agricultural sector. Therefore, the purpose of this study was to make the calculation reliable even if only the temperature data were present. In this research, the ET0 was initially calculated for a limited number of weather stations (12) using the Penman–Monteith method. In some cases, the simplified Penman–Monteith formula was adopted, while in others, as in the case of mountain weather stations, the complete formula was employed to consider the differences in vegetation, deduced from satellite surveys. Subsequently, the ET0 was calculated with the Hargreaves–Samani (HS) formula, calibrating the Hargreaves coefficient, through the spatialization of ET0, by the geostatistical method. The results showed a high reliability of the HS method in comparison with simplified PM (PM) method, and complete Penman–Monteith (cPM) method, with a minimum calibration of the empirical Hargreaves coefficient. In particular, a very good correlation between the results obtained in the mountain environment with the uncalibrated HS method and the cPM method was also observed in this area, while PM showed discordant and much higher results than ET0 compared with the other methods. It follows that this procedure allowed a more accurate estimate of potential evapotranspiration with a view to territory management, both in terms of water resources and the irrigation needs of the vegetation.

The extreme rainfall events of recent years in central Italy are producing an increase in hydrogeological risk, with disastrous flooding in terms of human lives and economic losses, as well as triggering landslide phenomena in correspondence with these events. A correct prediction of 100-year return levels could encourage better land planning, sizing works correctly according to the expected extreme events and managing emergencies more consciously through real-time alerts. In the recent period, it has been observed that the return levels predicted by the main forecasting methods for extreme rainfall events have turned out to be lower than observed within a few years. In this context, a model widely used in the literature, the generalised extreme value (GEV) with the “block maxima” approach, was used to assess the dependence of this model on the length of the collected precipitation time series and the possible addition of years with extreme events of great intensity. A total of 131 rainfall time series were collected from the Adriatic slope in central Italy comparing two periods: one characterised by 70 years of observations (1951–2020), the other by only 30 years (1991–2020). At the same time, a decision was made to analyse what the effect might be—in terms of the 100-year return level—of introducing an additional extreme event to the 1991–2020 historical series, in this case an event that actually occurred in the area on 15 September 2022. The results obtained were rather surprising, with a clear indication that the values of the 100-year return level calculated by GEV vary according to the length of the historical series examined. In particular, the shorter time series 1991–2020 provided higher return level values than those obtained from the 1951–2020 period; furthermore, the addition of the extreme event of 2022 generated even higher return level values. It follows that, as shown by the extreme precipitation events that have occurred in recent years, it is more appropriate to consider a rather short period because the ongoing climate change does not allow true estimates to be obtained using longer time series, which are preferred in the scientific literature, or possibly questioning the real reliability of the GEV model.

The Italian territory is subject to a high level of hydrogeological instability that periodically results in the loss of lives, buildings and productive activities. Therefore, the recognition of areas susceptible to hydrogeological instability is the basis for preparing countermeasures. In this context, landslide susceptibility in the mid-Adriatic slope was analyzed using a statistical method, the weight of evidence (WoE), which uses information from several independent sources to provide sufficient evidence to predict possible system developments. Only flows, slides, debris flows and mud flows were considered, with a total of 14,927 landslides obtained from the IFFI (Inventory of Franous Phenomena in Italy) database. Seven climatic–environmental factors were used for mapping landslide susceptibility in the study area: slope, aspect, extreme precipitation, normalized difference vegetation index (NDVI), CORINE land cover (CLC), and topographic wetness index (TWI). The introduction of these factors into the model resulted in rasters that allowed calculation by GIS-type software of a susceptibility map. The result was validated by the ROC curve method, using a group of landslides, equal to 20% of the total, not used in the modeling. The performance of the model, i.e., the ability to predict the presence or absence of a landslide movement correctly, was 0.75, indicating a moderately accurate model, which nevertheless appears innovative for two reasons: the first is that it analyzes an inhomogeneous area of more than 9000 km², which is very large compared to similar analyses, and the second reason is the causal factors used, which have high weights for some classes despite the heterogeneity of the area. This research has enabled the simultaneous introduction of unconventional factors for landslide susceptibility analysis, which, however, could be successfully used at larger scales in the future.