Rainfall Parameters A ﬀ ecting Splash Erosion under Natural Conditions

Featured Application: Practical and cost-e ﬀ ective splash cup method for splash erosion measurements in ﬁeld and laboratory conditions. Abstract: The interaction between rainfall erosivity parameters and splash erosion is crucial for describing the soil erosion process; however, it is rarely investigated under natural rainfall conditions. In this study, we conducted splash erosion experiments under natural rainfall on three sites in Central Europe. The main goal was to obtain the relationship between splash erosion of the bare soil in seedbed condition and commonly used rainfall erosivity parameters (kinetic energy, intensity, and rainfall erosivity ( EI 30 )). All sites were equipped with a rain gauge and an optical laser disdrometer where the splash erosion was measured, with modiﬁed Morgan splash cups. In order to investigate which parameter best describes the splash erosion process for all sites, a regression analysis was performed. In total, 80 splash erosion events were evaluated. Splash erosion can be described as a linear function of total kinetic energy and a non-linear function of EI 30. However, the use of the total kinetic energy led to underestimation of the splash erosion rates for highly intensive rainfalls. Therefore, better results were obtained when using average rainfall intensity as the splash erosion predictor or the kinetic energy divided by the rainfall duration. Minor di ﬀ erences between the replicates during splash erosion measurements indicate that the modiﬁed Morgan splash cup provides a good tool for soil erosion assessment.


Introduction
Soil erosion by water is the most common soil degradation process globally, and in arable cropland it is consistently higher than soil formation [1]. Detailed knowledge of the processes that control erosion on arable croplands contributes to better application of soil management techniques that minimize and control soil erosion risk [2]. Splash erosion starts with the raindrop impact on the soil surface, which represents the first stage in soil erosion by water [3]. The detached soil particles transported by raindrop impact are deposited on the near-distance soil surface or are transported further by surface runoff if the infiltration capacity of the soil is reached [4]. The major splash erosion driver is the erosivity of rainfall, which can be expressed by parameters like rainfall intensity [5][6][7], kinetic energy (KE) [8][9][10], rainfall erosivity (EI 30 ) [11], or raindrop momentum [12]. Apart from the rainfall properties, the detachment of soil particles also depends on soil physico-chemical characteristics, such as infiltration capacity [13,14], initial water content [15,16], the ability to form stable aggregates and crusts [17,18], organic matter content [19], texture, cohesion, porosity, capacity of ionic interchange, and clay content [20].
Splash erosion measurements on a small scale are usually done through splash cups or splash containers [21][22][23][24][25][26]. Most of the splash erosion studies have been conducted in the laboratories using rainfall simulators, where by controlling the raindrop size and fall height, the KE of raindrops is adjusted [12,[27][28][29][30]. However, rainfall simulators often do not reproduce the same rainfall drop size and velocity distribution characteristics as in nature [31]. As the velocity of raindrops is controlled by the height at which the nozzles are located, due to space or design limitations, sometimes raindrops cannot reach the terminal velocity of natural raindrops [32]. By applying different pressure at the nozzle, the raindrop velocity can be adjusted regarding to the raindrop size; however, large drops are unlikely to reach their terminal velocity, and consequently the KE, of natural rainfall [33]. Nevertheless, laboratory experiments improve the consistency of the results by minimizing the effects of the various uncontrolled factors that are present in the field [34], and also allow experiments to be repeated.
Splash erosion experiments under natural rainfall investigate the relationship between rainfall erosivity and splash detachment [3]. Morgan [23] observed the splash erosion under natural rainfall for 100 consecutive days, comparing four different soil textures. The KE of rainfall was calculated from the 10 min rainfall intensity values, using the formula from Hudson [35]. Splash erosion of the bare soil was significantly correlated with KE. Govers [36] collected data at 21 sites in Belgium using circular splash cups. He found that the product of rainfall KE and drop circumference are better at expressing the rainfall erosivity compared to KE and intensity, or when the 0.75 power of rainfall intensity is used. However, a detailed drop size distribution (DSD) was not available at the time, and the fall velocities of raindrops were based on data by Laws [37]. Splash erosion under natural conditions is primarily affected by rainfall DSD. The ability of raindrop impact to cause splash erosion (rainfall erosivity) is mainly dependent on drop size and drop fall velocity [12]. Direct measurements of raindrop size and velocity provide precise information about the erosivity of rainstorms-namely, KE. When the raindrop size and velocity is not directly measured, the rainfall KE is estimated from the experimentally based equations between rainfall intensity and KE from other studies. Theoretically obtained rainfall KE could underestimate or overestimate the real KE [38][39][40]. Furthermore, DSD obtained from other studies can significantly vary depending on rain type and geographical location [41]. With the development of optical laser techniques (disdrometer), the continuous and direct measurement of raindrop size and velocity has become easily available to assess rainfall KE.
A recent study with splash erosion measurements under natural rainfall, using the splash cup technique and rainfall monitoring with a disdrometer, was performed by Fernández-Raga et al. [42]. They used a funnel and cup installed directly in the field for splash erosion measurements, and found a good correlation between splash erosion and rainfall KE; however, their findings were based on only nine sampling periods. Angulo-Martínez et al. [43] conducted a study in Spain where the splash erosion of three soil types was measured with Morgan splash cups [23]. A significant relationship was found between splash erosion and the rainfall erosivity index EI 30 , and high variabilities between the replicates indicated the heterogeneity in splash erosion spatial distribution. According to the results reported from these studies, there are still many uncertainties concerning the changes in surface condition and spatial distribution of splash erosion.
The studies investigating splash erosion under natural rainfall are limited to local conditions. Consequently, monitoring of the rainfall characteristics on higher temporal and spatial resolution is crucial for describing the dominant rainfall parameters on splash erosion related to a specific location. Apart from the field studies of Fernández-Raga et al. [42] and Angulo-Martínez et al. [43], there are very few experiments that include both the monitoring of splash erosion and rainfall characteristics, including DSD, in the same location. Considering the local influences and lack of the data sets on rainfall DSD, it is difficult to define the role of splash in soil erosion process and predict it relative to local conditions. Furthermore, Bauer [44] pointed out that many rain events in Central Europe do not generate overland flow, but splash erosion is initiated already from the first drop impact, which emphasizes the importance of this soil degradation process. Lack of knowledge about the effect of erosive rainfall events on splash detachment in the agriculturally active Central European area was the main motivation for the present study.
This study presents the results from the splash erosion measurements collected during three consecutive summer seasons at three sites in Central Europe. Together with splash erosion, rainfall parameters, including rainfall intensity and KE, were monitored at the sites, with the aim of analyzing performance of the most common rainfall erosivity parameters (KE, intensity, and rainfall erosivity (EI 30 )), in order to predict splash erosion under natural rainfall.

Study Sites
The monitoring of rainfall parameters and splash erosion measurements was located at two sites in Austria and one in the Czech Republic ( Figure 1) [48] in Lower Austria, most of the erosive storms occurred during the summer period. Panagos et al. [49] also confirmed the highest rainfall erosivity during summer in the Central European region. Therefore, the selected monitoring period was during late spring and summer from 2017 to 2019.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 3 of 17 data sets on rainfall DSD, it is difficult to define the role of splash in soil erosion process and predict it relative to local conditions. Furthermore, Bauer [44] pointed out that many rain events in Central Europe do not generate overland flow, but splash erosion is initiated already from the first drop impact, which emphasizes the importance of this soil degradation process. Lack of knowledge about the effect of erosive rainfall events on splash detachment in the agriculturally active Central European area was the main motivation for the present study. This study presents the results from the splash erosion measurements collected during three consecutive summer seasons at three sites in Central Europe. Together with splash erosion, rainfall parameters, including rainfall intensity and KE, were monitored at the sites, with the aim of analyzing performance of the most common rainfall erosivity parameters (KE, intensity, and rainfall erosivity (EI30)), in order to predict splash erosion under natural rainfall.

Study Sites
The monitoring of rainfall parameters and splash erosion measurements was located at two sites in Austria and one in the Czech Republic ( Figure 1). The Austrian sites, Petzenkirchen (48°9'17" N, 15°14'46" E) and Mistelbach (48°34'59" N, 16°35'15" E), are situated in the region of Lower Austria, where most of the country agricultural activity takes place. The Czech site was located in Prague (50°6'17" N, 14°23'15" E). The long-term, average annual precipitation for the three sites is 902, 537 and 459 mm for Petzenkirchen, Mistelbach, and Prague, respectively [45][46][47]. According to data measured by Klik and Truman [48] in Lower Austria, most of the erosive storms occurred during the summer period. Panagos et al. [49] also confirmed the highest rainfall erosivity during summer in the Central European region. Therefore, the selected monitoring period was during late spring and summer from 2017 to 2019.

Investigated Soils
Soil samples were taken from two locations in Austria close to the experimental sites, Zwerbach (ZW) ( 20.0" E). The locations of three sites are marked in Figure 1. The samples were collected from agricultural land in the first 10 cm during April 2017, after seedbed preparation. Soil was dried and sieved through a 10 mm sieve and distributed to the three experimental sites. Particle size distribution was determined with a combined wet sieving and sedimentation method, as defined in the Austrian Norm for soil physical analysis [50,51]. Accordingly, soil textures were determined using the Austrian soil texture triangle [52]. Soil pH was obtained with electrometric pH meter. Total organic carbon was measured according to the Austrian Norm for the determination of soil organic carbon [53], and the aggregate stability of soils was determined with the modified Kemper and Rosenau method [54]. The physical and chemical properties of the soils obtained in the laboratory analysis are listed in Table 1.

Splash Erosion Measurements
Splash erosion was measured with the splash cup technique proposed by Morgan [23]. The splash cup was produced from a light polypropylene material, with an inner diameter of 10.3 cm and a standing height of 6 cm. On the bottom of the splash cup, holes were drilled to ensure water drainage through the soil; however, two fine meshes (500 and 1000 µm) were placed on the bottom to prevent soil loss trough the holes (Figure 2a). Air-dried and sieved soil (<10 mm) was filled in three layers up to 1 cm below the splash cup edge, to prevent the overflow of soil on the surface during the high intensity rainfall. While levelling it continuously using a long needle, each soil layer was slightly compacted to reach similar conditions (bulk density) as in the field. Aggregates at the top layer surface were randomly distributed to achieve heterogeneous arrangement. For each soil, more or less the same mass was filled into the splash cups, to keep soil density (in seedbed condition) constant within the replicates ( Figure 2b). Major differences in the soil structure of the samples prepared for experiment and the original soil in seedbed condition are the aggregate size and their arrangement, due to sieving and sample preparation. However, soil bulk density and porosity were within the same range of the soil in the seedbed condition in the field. Splash cups filled with soils were placed in the middle of a splash collector, with standing height of 30 cm and diameter of 47 cm (Figure 2c). The splash collector had an outlet, ensuring the drainage of rainfall water with splashed soil into collectors placed underneath ( Figure 2d). The water collector underneath was completely closed, ensuring that only splashed soil was trapped into the collector. Splash cups were installed in the field in such a way that the soil surface was 1 m above ground level, which corresponds to the same height as the disdrometer (described in next chapter). A detailed splash cup and splash collector design can be found in Zumr et al. [55].
diameter of 47 cm (Figure 2c). The splash collector had an outlet, ensuring the drainage of rainfall water with splashed soil into collectors placed underneath (Figure 2d). The water collector underneath was completely closed, ensuring that only splashed soil was trapped into the collector. Splash cups were installed in the field in such a way that the soil surface was 1 m above ground level, which corresponds to the same height as the disdrometer (described in next chapter). A detailed splash cup and splash collector design can be found in Zumr et al. [55].  Splash cups were exchanged, and splashed sediment was collected after each rainfall event, with accumulated precipitation of 5 mm. The events with precipitation lower than 5 mm were considered as erosive if the accumulated rainfall of 2.5 mm was reached after 15 min. This threshold was adopted based on preliminary data analysis and recommended by some European authors [44,56]. The intensity threshold of 12.7 mm (~0.5 in), reported in the Revised Universal Soil Erosion Equation (RUSLE) [11], is too high for European conditions. A rainfall event was defined as the rain period separated between proceeding and succeeding rainfall by 6 h or more. The splashed soil particles on the rim of the splash collector were completely washed off and drained into the water collector with rainfall and sediment. In the laboratory, splashed particles were filtered from water and oven-dried at 40 • C. The average mass of the splashed soil per each rainfall event was calculated from the three or four replicates, depending on the study site.

Rainfall Parameters
At all sites, rainfall data were collected with a rain gauge (tipping bucket or balance principle). To obtain the KE of rainfall, each site was equipped with an optical laser disdrometer, which measures raindrop size and velocity distribution in one-minute intervals. The PWS100 Present Weather Sensors from Campbell Scientific (PWS100) were installed in Mistelbach and Petzenkirchen, and the Laser Precipitation Monitor from Thies Clima (Thies) in Prague. The devices differ in measurement principle, sampling area, and drop size and velocity distribution classes. The PWS100 sampling area was 40 cm 2 , in which the drops are categorized in 34 size and 34 velocity classes. The Thies had a sampling area of 44.1 cm 2 , with 22 drop sizes and 20 velocity classes. Both disdrometer types differentiate raindrops from hail, ice pellets, and graupel. A detailed description of the disdrometer and rainfall monitoring set-ups at the three sites is given in Johannsen et al. [57].
The rainfall kinetic energy per area, KE (J m −2 ), was calculated for the diameter class i and velocity class j that are provided by the disdrometer, as follows: where N i,j is number of detected raindrops of a certain size class i and velocity class j; A is the sampling area of the disdrometer (m 2 ); ρ is density of water (g cm −3 ); D i is mean drop diameter (mm) of size class i; and v j is mean fall velocity (m s −1 ) of velocity class j. The mass of the raindrop is calculated assuming a spherical drop shape. Total KE is the sum of kinetic energies for each drop size and velocity, multiplied by the number of drops in the corresponding classes. Cumulative KE of a single rainfall event, KE sum (J m −2 ), was defined as Appl. Sci. 2020, 10, 4103 6 of 17 where KE i represents the i-th minute of the rainfall event. The KE per rainfall duration was calculated as where the KE h is kinetic energy per hour (J m −2 h −1 ) and T is total duration of a rainfall event (h), measured from the beginning of rainfall. The results obtained in the recent study by Johannsen et al. [57] showed that the differences in rainfall measurements from different disdrometer types, including PWS100 and Thies, influenced the interpretation of rainfall KE. That study concluded that the Thies disdrometer measured higher numbers of smaller drops; therefore, it underestimated the rainfall KE. To ensure comparable results obtained from the two disdrometer types used in this study, we applied the suggested correction factor by Johannsen et al. [57] for the Thies disdrometer. The correction is the ratio of the slope factors obtained from the KE and the intensity relationships for the PWS100 and Thies disdrometers. Accordingly, the final KE for the Thies disdrometer is calculated as in Equation (2) and multiplied with a correction factor of 1.36.
For the periods when the disdrometer data were not available due to disruptions or errors, the KE was calculated based on KE and intensity relationships developed for the three sites, according to the data investigated in Johannsen et al. [40]. The KE is presented as the function of the intensity, which was available from the rain gauge.
The average rainfall intensity I av (mm h −1 ) of single rainfall event considered the amount of rainfall over a time period, and is expressed as where P is the cumulative rainfall precipitation (mm) over rainfall duration T (h). The rainfall erosivity EI 30 (MJ mm ha −1 h −1 ) [11] for soil erosion estimation was also considered in analysis as an important rainfall parameter affecting splash erosion. The erosivity of a single rainfall event has been defined by the following expression: where I 30 is the maximum 30 min rainfall intensity of one rainfall event (mm h −1 ). Rainfall erosivity density expresses the mean rainfall erosivity per rainfall unit [49]. The mean monthly erosivity density MED (MJ ha −1 h −1 ) is expressed as: where n is the number of rainfall events (i) recorded during one month. The total splash erosion rate per one measurement period was calculated as the total amount of the splashed material from the test area per unit area, which is calculated as follows: where S is average splash erosion rate (g m −2 ), m s is the mass of splashed material, and A s is the area of the splash cup (0.0084 m 2 ). Splash erosion S h (g m −2 h −1 ) per rainfall event duration (T) was calculated as

Data Analysis
The coefficient of determination (R 2 ) and the root mean square error (RMSE) were used to validate the efficiency of the influencing rainfall erosivity parameters as the predictors of splash erosion. Pearson correlation analysis was used to estimate the relationship between the splash erosion rates and the sand, silt, or clay content of the three soils. The Kruskal-Wallis [58] test by ranks was used to determinate the differences in splash rates between the three soils. During the measuring period, the Mistelbach site had highest average monthly rainfall, KE, EI30, and erosivity density, as well as the largest mean drop diameter (Table 2). Petzenkirchen was the site with the noticeably lowest average monthly EI30 and erosivity density. The Prague site had lowest monthly precipitation among the three sites. However, higher average monthly EI30 and erosivity density compared to the Petzenkirchen site were caused by intensive storms recorded in June and August 2018 (Figure 3). It should be noted that calculated mean drop diameter could be affected by During the measuring period, the Mistelbach site had highest average monthly rainfall, KE, EI 30 , and erosivity density, as well as the largest mean drop diameter (Table 2). Petzenkirchen was the site with the noticeably lowest average monthly EI 30 and erosivity density. The Prague site had lowest monthly precipitation among the three sites. However, higher average monthly EI 30 and erosivity density compared to the Petzenkirchen site were caused by intensive storms recorded in June and August 2018 (Figure 3). It should be noted that calculated mean drop diameter could be affected by differences in measured drop size distributions between the disdrometer types used in the study [40].

Splash Erosion as the Function of Total Kinetic Energy
During the measuring period, 99 splash erosion records of the three soils were obtained from the investigated sites. After the data evaluation, a total of 80 records that had complete rainfall and splash erosion measurements were selected for further data analysis.
For measured KE sum , the mean splash erosion rates ranged between 4 and 2503 g m −2 for ZW, 5 to 1972 g m −2 for MI, and 12 to 2508 g m −2 for BK soils (Figure 4). The variabilities between the splash erosion replicates were larger for the measurements of KE sum above 780 J m −2 .

Splash Erosion as the Function of Total Kinetic Energy
During the measuring period, 99 splash erosion records of the three soils were obtained from the investigated sites. After the data evaluation, a total of 80 records that had complete rainfall and splash erosion measurements were selected for further data analysis.
For measured KEsum, the mean splash erosion rates ranged between 4 and 2503 g m −2 for ZW, 5 to 1972 g m −2 for MI, and 12 to 2508 g m −2 for BK soils (Figure 4). The variabilities between the splash erosion replicates were larger for the measurements of KEsum above 780 J m −2 . Significant (p < 0.05) linear regression was obtained for the three soils, as the splash erosion for each soil was positively correlated to KEsum with R 2 of 0.52, 0.50, and 0.45 for the BK, ZW, and MI soils, respectively (Table 3). According to the linear model, BK soil yielded the highest splash erosion, followed by ZW and MI soils. RMSE indicated a high deviation of measured values and values predicted by linear relationships.   Significant (p < 0.05) linear regression was obtained for the three soils, as the splash erosion for each soil was positively correlated to KE sum with R 2 of 0.52, 0.50, and 0.45 for the BK, ZW, and MI soils, respectively (Table 3). According to the linear model, BK soil yielded the highest splash erosion, followed by ZW and MI soils. RMSE indicated a high deviation of measured values and values predicted by linear relationships. Table 3. Main outputs of the regression analysis between splash erosion S (g m −2 ) for Zwerbach (ZW), Mistelbach (MI), and Bílkovice (BK) soils and rainfall event cumulative kinetic energy KE sum (J m −2 ). R 2 is the determination coefficient of the regression model, and RMSE (g m −2 ) indicates the root mean squared error.

Impact of Rainfall Intensity on Splash Erosion
The points deviating above the fitted linear regression line for the BK soil were mostly measured in Mistelbach (Figure 4). Considering the site-specific differences ( Table 2), it was necessary to distinguish the erosive events with higher rainfall intensities from the ones with low intensities having the same KE sum . The example in Figure 5 shows two rainfall events with similar cumulative precipitation and KE sum measured in Mistelbach (Event 1) and Prague (Event 2) with the corresponding splash erosion. The rainfall Event 1 was recorded on June 6, 2018, with a duration of 3.5 h. The rainfall Event 2 was recorded on August 8, 2018, with a duration of 43 min. Cumulative rainfall reached 19 and 20 mm for Event 1 and Event 2, respectively. Corresponding KE sum values were 467 J m −2 for Event 1 and 446 J m −2 for Event 2. However, during Event 2, up to 86% higher splash erosion rates were measured. Therefore, to compare the effect of kinetic energies characterized by different rainfall intensities on splash erosion, the KE sum was divided by rainfall duration (T) (defined in Equation (3)).
Appl. Sci. 2020, 10, x FOR PEER REVIEW 9 of 17 characterized by different rainfall intensities on splash erosion, the KEsum was divided by rainfall duration (T) (defined in Equation (3)).

Regression Analysis of Rainfall Erosivity Parameters with Splash Erosion
Previous results have indicated that the KEsum could not represent the realistic ability of rainfall to produce splash erosion. Since rainfall intensity plays an important role in rainfall erosivity,

Regression Analysis of Rainfall Erosivity Parameters with Splash Erosion
Previous results have indicated that the KE sum could not represent the realistic ability of rainfall to produce splash erosion. Since rainfall intensity plays an important role in rainfall erosivity, together with I av and EI 30 the results of regression analysis with KE h are presented in the following section. KE h represents the average energy load per rainfall duration.

Splash Erosion and Kinetic Energy Per Rainfall Duration
The relationship between KE h and splash erosion resulted in a non-linear (power) regression function, with an R 2 of 0.75 for MI soil and 0.76 for ZW and BK soil (Figure 6a, Table 4). The results indicated a less scattered distribution compared to results for KE sum (Figure 4). However, several observations in Figure 6a present a difference of 97% between the lowest and the highest splash erosion rate for a similar range of KE h (220-280 J m −2 h −1 ). The splash erosion measurement with KE h of 220 J m −2 h −1 (E 1 on Figure 6a) was the result of a rainfall event with short duration. Peak intensity of 42 mm h −1 was reached at the beginning of a rainfall event, with a short duration of 3 min. This might indicate that the peak intensity duration was too short to produce higher splash erosion rates.  Excluding these observations from the regression analysis, the resulting relationship between splash erosion and KEh was linear, with R 2 values of 0.91 for ZW and MI soil, and 0.90 for BK soil (Figure 6b; Table 4). The difference of 62% between the RMSEs for linear and nonlinear models indicates better performance of the splash erosion-KEh relationship without extreme observations. According to RMSE, correlation coefficients, and slope coefficients of regression equations, fewer minor differences between the soils were noted for the linear than for non-linear regression.    Excluding these observations from the regression analysis, the resulting relationship between splash erosion and KE h was linear, with R 2 values of 0.91 for ZW and MI soil, and 0.90 for BK soil (Figure 6b; Table 4). The difference of 62% between the RMSEs for linear and nonlinear models indicates better performance of the splash erosion-KE h relationship without extreme observations. According to RMSE, correlation coefficients, and slope coefficients of regression equations, fewer minor differences between the soils were noted for the linear than for non-linear regression.

Splash Erosion and Mean Rainfall Intensity
Splash erosion plotted against the I av (calculated according to Equation (4)), resulted in linear relationships for the three soils (Figure 7). Compared to previous results with KE h , most of the data were grouped in the range of low intensities up to 10 mm h −1 , and the highest I av corresponds to a KE h of 650 J m −2 h −1 . The three extreme measurements from the previous example with KE h (E 1-3 ) were more linearly distributed with increasing I av ; however, they still deviate from the regression lines (Figure 7a). The regression analysis without extreme observations did not greatly affect the linear regression for the three soils; nevertheless, R 2 increased for all soils, with the highest increase from 0.81 to 0.93 calculated for the ZW soil (Table 5). Table 5. Main outputs of the regression analysis between mean splash erosion Sh (g m −2 h −1 ) for Zwerbach (ZW), Mistelbach (MI), and Bílkovice (BK) soil, and mean rainfall intensity Iav (mm h −1 ). The R 2 is the determination coefficient of the regression model, and RMSE (g m −2 h −1 ) indicates the root mean squared error. E1, E2, and E3 denote extreme observations. The relationship between splash erosion and EI30 resulted in a non-linear power function (Figure 8a), with R 2 values of 0.60 for ZW, 0.64 for MI, and 0.65 for BK soil, including the extreme observations ( Table 6). Most of the observations were grouped up to an EI30 value of 250 (MJ ha mm −1 h −1 ), with high variations in splash erosion (from 5 to 1586 g m −2 ) between the single observations. Minor differences in R 2 and RMSE between the analysis with and without extreme observations The regression analysis without extreme observations did not greatly affect the linear regression for the three soils; nevertheless, R 2 increased for all soils, with the highest increase from 0.81 to 0.93 calculated for the ZW soil (Table 5). Table 5. Main outputs of the regression analysis between mean splash erosion S h (g m −2 h −1 ) for Zwerbach (ZW), Mistelbach (MI), and Bílkovice (BK) soil, and mean rainfall intensity I av (mm h −1 ). The R 2 is the determination coefficient of the regression model, and RMSE (g m −2 h −1 ) indicates the root mean squared error. E 1 , E 2 , and E 3 denote extreme observations.

Splash Erosion and Rainfall Erosivity (EI 30 )
The relationship between splash erosion and EI 30 resulted in a non-linear power function (Figure 8a), with R 2 values of 0.60 for ZW, 0.64 for MI, and 0.65 for BK soil, including the extreme observations ( Table 6). Most of the observations were grouped up to an EI 30 value of 250 (MJ ha mm −1 h −1 ), with high variations in splash erosion (from 5 to 1586 g m −2 ) between the single observations. Minor differences in R 2 and RMSE between the analysis with and without extreme observations indicated better correlation of extreme splash erosion observations to EI 30 ( Table 6). The highest splash erosion rates were indicated for BK soil, similar to the results obtained with KE sum seen in Figure 4.

Discussion
Comparable studies to our splash erosion experiments under natural rainfall were made in Portugal by Fernández-Raga et al. [42], and in Spain by Angulo-Martínez et al. [43]. The Thies disdrometer was used in both studies to directly assess rainfall KE. The splash erosion rates measured by Fernández-Raga et al. [42] were between 2.3 and 100 g m −2 . In the same range of total KE measured at our sites, splash erosion for loamy sand soil, which was most similar to the texture from the study in Portugal, was between 12 and 2508 g m −2 . However, the Portugal study was based on only nine splash erosion records, during which low rainfall intensities characterized by small raindrops (<0.55 mm) were measured. According to findings by Bubenzer and Jones [59], smaller drops produce significantly less splash erosion than larger ones, even for the same amount of KE. This would explain the lower splash erosion rates compared to our measurements, where more erosive rainfall events with larger mean drop sizes (>0.6 mm) were measured. Furthermore, differences in splash erosion measuring principles could also play a role when comparing results. In Portugal, the cup was placed directly on the soil bed, and splashed particles were collected from the surrounding soil. We prepared the soil samples and measured the particles splashed into the collector surrounding the soil. Fernández-Raga et al. [42] Table 6. Main outputs of the regression analysis between mean splash erosion S (g m −2 ) for Zwerbach (ZW), Mistelbach (MI), and Bílkovice (BK) soil and rainfall erosivity factor EI 30 (MJ ha mm −1 h −1 ). R 2 is the determination coefficient of the regression model, and RMSE (g m −2 ) indicates the root mean squared error. E 1 , E 2 , and E 3 denote extreme observations.

Discussion
Comparable studies to our splash erosion experiments under natural rainfall were made in Portugal by Fernández-Raga et al. [42], and in Spain by Angulo-Martínez et al. [43]. The Thies disdrometer was used in both studies to directly assess rainfall KE. The splash erosion rates measured by Fernández-Raga et al. [42] were between 2.3 and 100 g m −2 . In the same range of total KE measured at our sites, splash erosion for loamy sand soil, which was most similar to the texture from the study in Portugal, was between 12 and 2508 g m −2 . However, the Portugal study was based on only nine splash erosion records, during which low rainfall intensities characterized by small raindrops (<0.55 mm) were measured. According to findings by Bubenzer and Jones [59], smaller drops produce significantly less splash erosion than larger ones, even for the same amount of KE. This would explain the lower splash erosion rates compared to our measurements, where more erosive rainfall events with larger mean drop sizes (>0.6 mm) were measured. Furthermore, differences in splash erosion measuring principles could also play a role when comparing results. In Portugal, the cup was placed directly on the soil bed, and splashed particles were collected from the surrounding soil. We prepared the soil samples and measured the particles splashed into the collector surrounding the soil. Fernández-Raga et al. [42] described splash erosion as the linear function of total KE, with the R 2 being 0.51 and 0.69 for different drop size and intensities thresholds used. That corresponds with our observations for the loamy sand soil, where an R 2 of 0.52 was obtained.
The study from Angulo-Martínez et al. [43] was more comparable to ours, considering that the Morgan splash cups were used for the splash erosion measurements. However, the samples were kept undisturbed during the whole monitoring period, whereas our samples were exchanged after each rainfall event. Splash erosion was measured for three soils with silt, sandy loam, and clay loam textures. The authors suggested EI 30 as a controlling factor for splash erosion where no differences in detached rates between the soils were reported. Comparable splash erosion rates from our analysis with EI 30 were found up to 200 MJ ha mm −1 h −1 ; however, with increasing EI 30 , our splash erosion rates increased up to 2500 g m −2 , whereas the rates from Spain remained constant with an average rate of 337 g m −2 . The fact that the samples were exchanged between measurements may contribute to higher rates obtained for soils in our study, which was in seedbed conditions.
There is still no general agreement on which rainfall parameters define splash erosion [12]. Parameters dependent on raindrop size and fall velocity, such as rainfall KE, momentum, intensity, or a combination of these, are commonly used to describe the raindrop impact on splash detachment. According to our analysis, KE sum could not (Figure 4) explain the variabilities between splash erosion rates obtained for the same amount of KE. The reason for that lies in different rainfall intensities between the rainfall events, where high-intensity rainfall produced more splash erosion than low-intensity rainfall ( Figure 5). In the field study by Govers [36], it was also concluded that the use of KE as an estimate of the rainfall detachment power leads to an underestimation of the relative impact of events with high intensities. From the strong linear relationship between splash erosion and rainfall intensity obtained in our study (Figure 7, Table 5), it can be stated that the splash erosion was more related to rainfall intensity than to other analyzed parameters (KE sum , KE h and EI 30 ). Nevertheless, we found good agreement between the splash erosion and KE sum divided by rainfall duration (T) ( Figure 6, Table 4). This indicates that KE can also be used as the parameter to predict splash erosion, even when events with different rainfall intensities are analyzed, but its erosive impact has to be expressed through rainfall duration.
The detailed information about drop size distribution allowed us to discern the differences in rainfall characteristics between the study sites. The differences in the splash rates between the rainfall events for the same range of KE also contributed to the differing drop size distribution. This was also noticed for the splash rates measured at the Mistelbach site, which is characterized as the site with the highest average raindrop diameter. Another example of this is the extreme event (E 2 ) reported in the results (Figure 6), where the high splash erosion rates were affected by the large drop size measured for this event. Bubenzer and Jones [59] found that rainfall with larger drops produce more detachment than rainfall with smaller drops, for rainfall having the same total KE. Recently, Fu et al. [60] also reported the gradual increase of splash erosion rates with increasing raindrop diameter. Detailed information about raindrop size distribution plays an important role for splash erosion studies like ours, where the direct measurements of the parameters is needed to describe the factors affecting the splash erosion process.
The soil's physical characteristics (texture, soil moisture, organic matter, structure, infiltration capacity, etc.) play an important role in understanding the soil detachment by raindrop splash [2]. Splash erosion of the three soils was positively correlated to the sand content, and significantly (p < 0.05) negatively correlated to clay content. For this reason, cumulative splash erosion rates were highest for the BK soil with highest sand content, and significantly different (p < 0.05) from ZW soil with the highest clay content (Table 1). Equivalent results were reported in a recent study by Zambon et al. [61], using simulated rainfall on same soils. The high splash detachability of soils with dominating sand content was also confirmed in experiments by Salles et al. [3], Cheng et al. [19], and Xiao et al. [62]. However, the results of the regression and correlation analysis between splash erosion per rainfall duration and KE h , as well as I av , indicates small differences between the three soils. Other soil properties, such as soil moisture, also have a significant impact on splash erosion [17,63,64]. The results reported by Zambon et al. [61] show that lower splash erosion rates are related to high initial soil water content, followed by surface ponding and changes in saturated hydraulic conductivity induced by surface crusting under high rainfall intensities. Although it was not possible to monitor the changes in soil moisture and surface conditions in the field, these effects probably contributed to results obtained in this field study, especially for the extreme observation E 3 . During this observation, two rainfall sub-events were recorded. The second major rainfall sub-event (with total rainfall of 40 mm) occurred 48 h later. Therefore, low splash erosion rates could be related to the long drying period between the first and second rainfall-sub event, resulting in increasing soil surface resistance against the raindrop impact [65]. A more detailed study, including the temporal monitoring of surface changes and soil moisture properties, would possibly contribute to clarifying the complex interaction between soil properties and rainfall controlling the splash erosion process.
Apart from the differences in the rainfall characteristic and soil properties, the experimental design for splash erosion assessment plays an important role when comparing the results from different studies. Recently, a study was published by Fernández-Raga et al. [25] that compared different devices for splash erosion measurements, where the results were strongly affected by the measurement device. This was also visible when comparing our results to the above-described studies. However, low standard deviations between the replicates for each soil obtained in our study confirm that the modified version of the Morgan splash cup provided reliable results for splash erosion measurements.

Conclusions
This study analyzed the effect of rainfall parameters on splash erosion under natural rainfall on three sites in Central Europe. Based on the results obtained during the three seasons of measurements, we conclude that splash erosion is more dependent on rainfall intensity than on total kinetic energy and rainfall erosivity (EI 30 ). Still, the kinetic energy of a natural rainfall event can be used as suitable erosivity parameter when dividing it by the rainfall duration. Monitoring of rainfall properties (intensity and kinetic energy) is important to discern the spatial and temporal differences in rainfall characteristics, which influence splash erosion. The dynamic changes in soil moisture, infiltration capacity, and surface roughness affected by weather conditions can lead to uncertainties in the evaluation of splash erosion in the field. Minor differences between the replicates during splash erosion measurements indicate that the modified Morgan splash cup provides a good tool for soil erosion assessment.