Search

Rill erosion is a major issue on a global scale, and predicting the presence, position, and development of erosive forms on hillslopes is a significant challenge for the scientific community. Several plot-scale investigations confirmed the reliability of the terrestrial photogrammetric (TP) technique for studying rill erosion and the reliability of a method for extracting the rill network from Digital Surface Models (DSMs) and measuring the corresponding volume. In this paper, for an intense erosive event that occurred at the Sparacia experimental area (Sicily, Southern Italy), TP surveys of three plots, with different length and steepness, incised by rills, were performed to reconstruct the DSMs. For each plot, the rill network was extracted from the DSMs, and the non-contributing network was distinguished from the contributing one, from which the soil loss and the consequent eroded volumes V were determined. The specific aims were to (i) establish the effect of plot steepness on rill depths and some morphometric characteristics of the drainage rill network; (ii) test and calibrate the relationship between V and the total rill length L, using all rill measurements available in the literature and those obtained in this study; and (iii) modify the V–L relationship by including climate forcing and assessing the related performance. The rill depths, h, the drainage frequency, and drainage density of the rill networks detected in the three plots were compared. The analysis demonstrated that h and the morphometric parameters of the contributing rill network increase with plot steepness s. In particular, the mean depth increases from 2.79 to 4.85 cm for slope increasing from 14.9 to 26%. Moreover, the drainage frequency of the contributing rill network varies from 0.16 m−2 for s = 14.9% to 0.47 m−2 for s = 26%, while the drainage density of the contributing rill network varies from 0.92 m−1 for s = 14.9% to 2.1 m−1 for s = 26%. Finally, using the data available in the literature and those obtained in this investigation, an empirical relationship between V and the total rill length L was firstly tested and then rearranged considering the event rainfall erosivity Re. Including Re in the rearranged equation guaranteed the best performance in V estimation.

The rainfall kinetic energy, which affects soil erosion processes, can be calculated by the drop-size distribution (DSD) and falling velocity. This study presents the outcomes derived by the DSDs recorded with the same optical disdrometer in two experimental areas, located in Sicily (southern Italy). Specifically, the DSDs were recorded from March 2017 to December 2019 at Sparacia and from June 2006 to April 2014 at Palermo. The aims of this paper are both to compare the DSDs for the two sites and to evaluate the applicability of Gamma theoretical distribution. Moreover, the relations of rainfall kinetic power vs. rainfall intensity are assessed. Differences in DSDs, especially for rainfall intensity I lower than 100 mm h−1, are detected, while they are less evident for the highest I values. The kinetic power per unit volume of rainfall varies between the two sites and can be considered roughly coincident only for I higher than 100 mm h−1. The observed kinetic power-rainfall intensity relations are not well described by the equations by Wischmeier and Smith and Brown and Foster. Finally, the reliability of a theoretical relationship for estimating the kinetic power by I and the median volume diameter D0 is positively tested for the Palermo dataset, while it is found to be reliable only for D0 greater than 0.17 cm for the Sparacia dataset.