In soil erosion estimation models, the variable with the greatest impact is rainfall erosivity (RE), which is the measurement of precipitation energy and its potential capacity to cause erosion, and erosivity density (ED), which relates RE to precipitation. The RE requires high temporal resolution records for its estimation. However, due to the limited observed information and the increasing availability of rainfall estimates based on remote sensing, recent research has shown the usefulness of using observed-corrected satellite data for RE estimation. This study evaluates the performance of a new gridded dataset of RE and ED in Peru (PISCO_reed) by merging data from the IMERG v06 product, through a new calibration approach with hourly records of automatic weather stations, during the period of 2000-2020. By using this method, a correlation of 0.7 was found between the PISCO\_reed and RE obtained by the observed data. An average annual RE for Peru of 4831 MJ·mm·ha-1·h-1 was estimated with a general increase towards the lowland Amazon regions and high values are found on the north-coast Pacific area of Peru. The spatial identification of the most risk areas of erosion, was carried out through a relationship between the ED and rainfall. Both erosivity data sets will allow us to expand our fundamental understanding and quantify soil erosion with greater precision.

Floods are a common natural disaster whose severity in terms of duration, water resource volume, peak and accumulated rainfall-based damage are likely to differ significantly for different geographical regions. In this paper, we first propose a novel hourly flood index (SWRI24−hr−S), derived from normalising the existing 24-hourly water resources index (WRI24−hr−S) in literature, to monitor flash flood risks on an hourly scale. The proposed SWRI24−hr−S is adopted to identify a flood situation and derive its characteristics, such as the duration (D), volume (V) and peak (Q). The comprehensive result analysis establishes the practical utility of SWRI24−hr−S in identifying flood situations at seven study sites in Fiji between 2014 and 2018 and deriving their characteristics (i.e., D, V, and Q). Secondly, this study develops a vine copula-probabilistic risk analysis system that models the joint distribution of flood characteristics (i.e., D, V, and Q) to extract their joint exceedance probability for the seven study sites in Fiji, enabling probabilistic flood risk assessment. The vine copula approach, particularly suited to Fiji’s study sites, introduces a novel probabilistic framework for flood risk assessment. The results show moderate differences in the spatial patterns of joint exceedance probability of flood characteristics in different combination scenarios generated by the proposed vine copula approach. In the worst-case scenario, the probability of any flood event occurring where the flood volume, peak and duration are likely to exceed the 95th-quantile value (representing an extreme flood event) is found to be less than 5% for all study sites. The proposed hourly flood index and the vine copula approach can be feasible and cost-effective tools for flood risk monitoring and assessment. The methodologies proposed in this study can be applied to other data-scarce regions where only rainfall data is available, offering crucial information for flood risk monitoring, assessment, and the development of effective mitigation strategies.

The study of precipitation amount at different time scales constitutes an important issue in climate research and risk assessment. When dealing with daily totals, a frequent but sometimes underestimated problem is at what time the observation day begins. The choice of different starting times may lead to incompatibility between stations and incorrect identification of extreme events. In this work, the problem of temporal misalignment between precipitation datasets characterized by different starting time of the observation day is analyzed. The most widely used adjustment methods (1day and uniform shift) and two methods based on reanalysis (NOAA and ERA5) are evaluated in terms of temporal alignment, precipitation statistics and percentile distributions. As test series, the precipitation amount collected from 9 a.m. local time (09 LT) on the previous day to 09 LT on the target day (9-9 datasets) of the Padua and nearby stations in the period 1993-2022 have been selected. Results show that the reanalysis-based methods, in particular ERA5, outperform the others in temporal alignment, regardless the station. But, for the periods in which reanalysis data are not available, “1day” method, which shifts the daily amount back one calendar day, and “unif” method, which distributes uniformly the daily total from a 2-day moving window surrounding the target date, can be considered valid alternatives. On the other side, concerning the precipitation statistics, the reanalysis-based methods are not the best option, as they increase the precipitation frequency and reduce the mean value over wet days, NOAA much more than ERA5. Nevertheless, the uniform method provides a larger deviation from the original daily series. The use of the series of a station nearby the target one, which is mandatory in case of missing data, gives similar or better results than applying any adjustment method to the 9-9 series. General conclusions can hardly be drawn as they depend on the method and station. For the Padua dataset, the analysis was repeated at monthly and seasonal resolution. In general, the adjustment series show the most relevant changes in the precipitation statistics in summer and less temporal alignment with the original series in summer and autumn, the two seasons mainly affected by heavy rains in Padua. Finally, the percentiles distribution, analyzed for all the methods and stations, indicates that any adjustment method underestimates the percentile values, except ERA5. Only Legnaro, the station most correlated with Padua, gives results like ERA5.