1. Introduction
Anthropogenic acceleration of soil erosion processes, due to tillage, nonconservative agricultural practices, overgrazing, and deforestation, has a strong impact on the soil ecosystem, and it affects the environmental quality [
1].
Soil erosion due to rainfall is a process characterized by the detachment of soil particles attributable to rainfall impact on the soil surface and the transport of the detached particles by rainfall and runoff.
Rainfall erosivity, which is defined as the capability of rainfall to determine soil particle detachment, is the variable used to characterize the erosive agent; it is necessary for quantifying soil erosion processes and can be expressed by the kinetic energy per unit time, called kinetic power
Pn [
2,
3,
4,
5]. The kinetic power of a given rainfall can be obtained by adding the contribution of each raindrop constituting the precipitation [
6], and it varies with the size, falling velocity, and number of impacting drops on soil surface area in unit time. In other words, the rainfall kinetic power is calculated by relating the drop-size distribution (DSD) to rainfall intensity and raindrop terminal velocity [
7,
8].
Several empirical equations estimating kinetic power by rainfall intensity,
I (mm h
−1), characterized by different mathematical shapes were proposed by different authors [
9,
10,
11,
12]. The most common empirical equation for the estimate of
Pn (J m
−2 h
−1) that considers its dependence on only rainfall intensity
I is the one presented by Wischmeier and Smith [
11] and applied in the Universal Soil Loss Equation (USLE):
where
It is the intensity threshold that Wischmeier and Smith [
11] assumed equal to 76 mm h
−1. Equation (1) demonstrates that the ratio
Pn/
I (J m
−2 mm
−1), which is the kinetic energy per unit volume of rainfall, grows for rainfall intensities ≤
It and becomes constant for rainfall intensities >
It. Wischmeier and Smith [
11] explained this threshold by considering that the median volume diameter,
D0, which is the diameter dividing the DSD into two parts of the same volume, does not grow with rainfall intensities > 76 mm h
−1.
Kinnell [
13] described the
Pn-
I relationship as follows:
in which
a,
b, and
c are parameters to be estimated.
Kinnell [
14] reported an
a value of 29 J m
−2 mm
−1 and considered
b and
c as site-specific parameters [
15]. As suggested by Brown and Foster [
16], Equation (2) is used in the RUSLE (Revised Universal Soil Loss Equation) model with
a = 29 J m
−2 mm
−1,
b = 0.72, and
c = 0.05 h mm
−1. For the RUSLE2 equation, Foster [
17] used Equation (2) with
a = 29 J m
−2 mm
−1,
b = 0.72, and
c = 0.082 h mm
−1 [
18].
Instead of using an empirical
Pn-
I relationship, the kinetic power of the precipitation can be obtained by adding the contribution of each raindrop constituting precipitation, if its DSD and the relationship between raindrop terminal velocity and drop diameter are known. The DSD establishing the number of drops
N(
D)
dD, having a drop diameter between
D and
D + dD, that hit a unit horizontal area in the unit time is useful to develop hydrological studies [
6,
19]. Ulbrich [
20] proposed, as DSD, a Gamma distribution that is characterized by wide applicability in several climatic conditions. This theoretical DSD is a generalization of the exponential distribution by Marshall and Palmer [
21], and it has the following form [
8]:
In which
μ (−),
Λ (cm
−1), and
No (m
−2 s
−1 cm
−1−μ) are the three parameters of the Gamma distribution and
N(D) dD is expressed as m
−2 s
−1.
Using the available experimental measurements of raindrop terminal velocity
V(
D) (m s
−1), Ferro [
22] deduced the following equation:
in which
Vn = 9.5 m s
−1,
an = 6 cm
−1, and
D is the drop diameter expressed as cm. Equation (4) was obtained using experimental measurements of raindrop terminal velocity carried out by Laws [
23], Gunn and Kinzer [
24], Blanchard [
25], Beard [
26], Epema and Riezebos [
27], and Jayawardena and Rezaur [
28] using drop diameters in the range 0.01–0.6 cm; it is characterized by a good agreement with the available measurements [
8,
22]. As the velocity
V(
D) increases for
D < 0.56 cm until it becomes constant, Equation (4) is also able to follow the physical boundary for
D greater than 0.55 cm (the drops are unstable and break up) [
8].
If both the DSD and the terminal velocity
V(
D) are known, the kinetic power
Pn (J s
−1 m
−2) can be obtained as follows [
15]:
where
ρ is the water density (kg m
−3), and
V(
D) is expressed as m s
−1 and
D as cm.
Carollo and Ferro [
6], coupling Equations (3)–(5), obtained the following expression of the kinetic power
Pn (J s
−1 m
−2):
in which
I (mm h
−1) is the rainfall intensity expressed as [
5,
15]
in which Γ is the gamma function.
Carollo et al. [
3] proved that, if the exponential distribution by Marshall and Palmer [
21] is referred to the unit area and time, then it becomes coincident with a Gamma distribution [
20] with
μ = 0.67 [
19], and Equation (6) can be rewritten as
The following relationship between the median volume diameter
D0 (cm) and the
Λ parameter [
4,
19,
20],
can be introduced into Equation (8) for deducing the following equation for estimating the kinetic power
Pn (J s
−1 m
−2):
>Equation (10) shows that
Pn/
I depends only on
D0. Hence, Equation (10) establishes that the whole DSD can be represented by
D0 when the rainfall kinetic power is calculated. According to Carollo et al. [
3,
5] and Serio et al. [
8], when the rainfall intensity
I >
It,
D0 and, according to Equation (10), the
Pn/
I ratio are constant. Thereby, Equation (10) theoretically confirms the hypothesis by Wischmeier and Smith [
11] that, for
I >
It, the
Pn/
I ratio becomes constant because
D0 has a constant value.
Firstly, in this paper, the DSDs recorded by the same optical disdrometer at the Sparacia experimental area from March 2017 to December 2019 and the Palermo experimental area from June 2006 to March 2014 are presented. The DSDs, aggregated into rainfall intensity classes characterized by different amplitudes, are used to test the reliability of Ulbrich’s distribution (Equation (3)) and to estimate its parameters μ and Λ. Then, the measured kinetic power values are used to assess the applicability of Equations (1), (2) and (10). The comparison between the observed relations of rainfall kinetic power vs. rainfall intensity allows for assessing if the well-known literature result, concerning the variation with the location of this type of relationship, is also valid for two sites located within the same region at such a short distance of about 100 km. The use of the same disdrometer in two different sites and the prolonged investigation period (11 years) make the developed analysis very robust and lend originality to the present investigation relative to previous studies on rainfall kinetic power performed elsewhere.
2. Materials and Methods
The DSDs were recorded using optical disdrometers (model ODM 70 made by Eigenbrodt) placed in two experimental Sicilian sites (Palermo and Sparacia) (
Figure 1).
The Sparacia experimental area (37.637300° N, 13.760671° E) is situated in western Sicily (in the Sicilian hinterland, away from the coast), southern Italy, approximately 100 km south of Palermo (
Figure 2a) at 415 m a.s.l. The climate is of Mediterranean semiarid type and is characterized by a 700-mm mean annual rainfall. The Palermo (northwestern Sicily, on the coast) experimental station is situated at the Department of Agricultural, Food and Forest Sciences of the University of Palermo (38.108171° N, 13.351861° E), at 40 m a.s.l. (
Figure 2b). The climate is Mediterranean temperate, with dry and hot summers and mild and rainy winters. The mean temperatures vary between 15°C and 22°C. Rainfalls mainly occur in autumn and winter and have a mean annual value equal to 654 mm.
For each minute of the rainfall event, the disdrometer (
Figure 3) measures the drop diameter in the range 0.035–0.6 cm. Every drop is separately measured and included in classes having about a 0.005-cm width. The disdrometer shares the diameter range into 128 classes and detects the number of drops that belong to a given class for every raining minute. The passage of the drop in the control volume between two diodes causes a light damping, which allows for measuring the drop diameter (
Figure 3). Further details are reported in Carollo et al. [
3,
4].
The disdrometer recorded 544 rainy days from June 2006 to April 2014 at the Palermo experimental area, and 318 rainy days at Sparacia from March 2017 to December 2019. For each registered DSD, rainfall intensity was calculated by summing the volumes of the single droplets under the hypothesis of a spherical shape.
4. Discussion
For a given rainfall intensity,
Figure 5 and
Figure 6 demonstrate that the two sites are characterized by similar values of
D25 and
D50 and, especially for
I < 100 mm h
−1, very different values of the characteristic diameters
D75 and
D0 that, according to Serio et al. [
29], correspond to even higher percentiles. Accordingly, these figures demonstrate that different DSDs, especially for smaller rainfall intensities, are registered in the two sites. For the highest values of
I, the differences between the DSDs registered in the two sites are less evident.
In particular,
Figure 6, confirming the similarity among the DSDs at the two locations for the highest rainfall intensity values, shows comparable values of the
D0 diameters for
I > 100 mm h
−1. In other words, the spatial variability within the Sicilian region affects the DSD up to a threshold intensity value of 100 mm h
−1, while for the highest
I values, this spatial variability of DSD is not clearly detected.
As detected for the relationship
D0-
I,
Figure 8 highlights that, in both sites, the kinetic energy per unit volume of rainfall,
Pn/
I, increases with rainfall intensity until a threshold value of rainfall intensity,
It, is reached. This threshold
It is equal to 40 mm h
−1 for the Palermo site, while it becomes 100 mm h
−1 for Sparacia. The comparison between
Figure 6 and
Figure 8 highlights that above the threshold
It, the median volume diameter
D0 assumes a constant value that can be considered very close for both sites. This result agrees with that obtained by Carollo et al. [
3,
5] and Serio et al. [
8], which established that when the rainfall intensity
I >
It,
D0 and, according to Equation (10), the
Pn/
I ratio are constant. This result also confirms the hypothesis of Wischmeier and Smith [
11] that
Pn/
I does not continue to increase for the highest intensity values (>
It) since
D0 is constant, too. According to Carollo et al. [
4], this result can be justified by considering that the unstable large drops break into small drops (disaggregation), and this circumstance, determining an increase of the number of drops, favors an effect of the collapse of drops. For low rainfall intensities (
I <
It), the drop collapse effect prevails on the disaggregation phenomenon, and the rainfall kinetic power increases more rapidly than rainfall intensity. Conversely, for high intensity (
I >
It), an equilibrium condition between disaggregation and collapse is reached; as a consequence, the DSD and
D0 become almost invariable, and the kinetic power is proportional to
I.
For rainfall intensity values
I >
It,
Pn/
I assumes a quasi constant value (27 J m
−2 mm
−1), which is very close to the asymptotic value obtained by Equation (1) (28.3 J m
−2 mm
−1) or Equation (2) (29 J m
−2 mm
−1). In any case,
Figure 8 clearly demonstrates that the relation
Pn/
I-
I is site specific.
For
I <
It and
I >
It,
Figure 11 shows the comparison between the parameters of Ulbrich’s distribution obtained by the momentum method using Equations (11) and (12). It is evident that, for both cases, the two parameters are related to each other even if this relation is described by different equations for
I <
It and
I >
It. According to Equation (6), this circumstance implies that the kinetic power per unit volume of rainfall can be determined by a single parameter, and this parameter could be the median volume diameter as established by Equation (10) and suggested by
Figure 10. Moreover, from a statistical point of view,
Figure 11 suggests that the DSDs can be described by two parameters alone, estimated, for example, by
I and
D0, but these estimating relationships differ between
I <
It and
I >
It. This result does not agree with that by Chandrasekar and Bringi [
30], who did not find any explanation concerning the relation between
μ and
Λ.
Figure 10 allows us to state that the hypotheses on which Equation (10) is based hold for the Palermo dataset, while for Sparacia measurements, Equation (10) can be applied only for
D0 > 0.17 cm, which, according to
Figure 6, corresponds to
I > 50 mm h
−1. However, for the Sparacia dataset, the (
D0,
Pn/
I) pairs are well aligned along a single curve whose equation is different from Equation (10). In other words, for the Sparacia site, a relationship between
Pn/
I and the only predictive variable
D0 can be identified, but the hypotheses adopted for deducing Equation (10) are no longer valid. In particular, the applicability of the Marshall and Palmer [
21] distribution is no longer verified, especially for the lower values of
D0.
Finally, all the above considerations highlight that the disdrometric characteristics of rainfall are site specific, especially for rainfall intensities less than 100 mm h
−1. The knowledge of rainfall intensity alone does not allow one to determine these characteristics or, consequently, to accurately estimate its erosivity, as clearly highlighted by
Figure 9.
5. Conclusions
Rainfall erosivity, which can be represented by its kinetic power, is a very important variable for defining erosion processes. In this investigation, the DSDs measured by the same optical disdrometer in two Sicilian experimental sites (Palermo and Sparacia) were compared. For a given value of rainfall intensity, the DSDs registered at Palermo are different from the DSDs measured at Sparacia even if the two sites are located in a similar geographical context.
The measured DSDs were used to deduce the kinetic power by coupling them with falling velocity estimated by the relationship proposed by Ferro [
22]. For a given value of rainfall intensity, the kinetic power per unit volume of rainfall registered in the two sites is very different and can be considered roughly coincident only for rainfall intensity higher than 100 mm h
−1.
This investigation clearly showed that the relation Pn/I-I is site specific, and a relationship such as Equation (1) or (2), which uses the rainfall intensity as a single predictive variable, is not applicable everywhere. In other words, the knowledge of rainfall intensity does not allow us to accurately estimate erosivity.
Finally, the relationship proposed by Carollo et al. [
3] (Equation (10)) is reliable for the Palermo dataset, while for the Sparacia dataset, it is reliable only for
D0 greater than 0.17 cm. For
D0 < 0.17 cm, the use of Equation (10) implies an appreciable systematic underestimation of
Pn/
I. This result suggests that the hypotheses on which Equation (10) is based are valid for the Palermo dataset, while they are valid only for
D0 > 0.17 cm for the Sparacia dataset.
Further investigations are necessary to verify the presented results in other climatic contexts.