Impact of Disdrometer Types on Rainfall Erosivity Estimation

Soil erosion by water is affected by the rainfall erosivity, which controls the initial detachment and mobilization of soil particles. Rainfall erosivity is expressed through the rainfall intensity (I) and the rainfall kinetic energy (KE). KE–I relationships are an important tool for rainfall erosivity estimation, when direct measurement of KE is not possible. However, the rainfall erosivity estimation varies depending on the chosen KE–I relationship, as the development of KE–I relationships is affected by the measurement method, geographical rainfall patterns and data handling. This study investigated how the development of KE–I relationships and rainfall erosivity estimation is affected by the use of different disdrometer types. Rainfall data were collected in 1-min intervals from six optical disdrometers at three measurement sites in Austria, one site in Czech Republic and one site in New Zealand. The disdrometers included two disdrometers of each of the following types: the PWS100 Present Weather Sensor from Campbell Scientific, the Laser Precipitation Monitor from Thies Clima and the first generation Parsivel from OTT Hydromet. The fit of KE–I relationships from the literature varied among disdrometers and sites. Drop size and velocity distributions and developed KE–I relationships were device-specific and showed similarities for disdrometers of the same type across measurement sites. This hindered direct comparison of results from different types of disdrometers, even when placed at the same site. Thus, to discern spatial differences in rainfall characteristics the same type of measurement instrument should be used.


Introduction
Soil erosion is recognized as a major threat to soil resources and the environment globally [1]. The initial detachment of soil particles through splash erosion, where soil particles are disaggregated and mobilized, is influenced by the rainfall erosivity. Erosivity is dependent on the rainfall characteristics, such as drop size and fall velocity, and the rainfall intensity (I). Thus, rainfall kinetic energy (KE) has been applied as a measure of potential soil loss [2]. Direct measurement of KE has been labor-intensive and costly, so as a means to calculate the kinetic energy of rainfall, when data of drop size and velocity are not available, KE-I relationships have been developed. The empirical relationships have been established from studies measuring the drop size and fall velocity at a specific site. The equations Other commonly used universal KE-I relationships are the Wischmeier and Smith [10] and Brown and Foster [16] equations as these have been used in erosivity estimations and to predict soil loss in the Universal Soil Loss Equation (USLE) and its revised version (RUSLE). A modified version of the Brown and Foster equation, as proposed by McGregor et al. [21], is now used in the updated RUSLE2 model ( Table 1). The (R)USLE is the most widely used quantitative model for soil erosion risk assessment and includes several parameters, one of them being the rainfall erosivity. The rainfall erosivity is defined as the product of the total event KE and the 30-minute maximum rainfall intensity, as rainfall KE alone may not provide the best estimation of erosivity [10]. However, if KE is estimated by an equation based on rainfall intensity, then the whole calculation of KE and rainfall erosivity is actually based on rainfall intensity [22].
As rainfall intensity, unlike KE, is an extensively measured rainfall parameter, this makes KE-I relationships useful for spatial and temporal rainfall erosivity estimation. This has enabled the mapping of rainfall erosivity to show global and regional erosion risk areas [23][24][25][26] and projections of future rainfall erosivity due to altered rainfall patterns [27].
The extensive use of these KE-I relationships in soil erosion risk assessment demonstrates the importance of choosing the most suitable model, as erosivity estimates will vary according to the chosen KE-I relationship [19,28,29]. Angulo-Martínez et al. [19] stated that disdrometer data are therefore necessary to check which KE-I relationship works best for the specific site conditions. Disdrometers are automated rainfall measurement instruments, which are able to continuously measure the number, size and velocity of falling raindrops and are thus able to give direct measurements of the kinetic energy. Disdrometer data are useful in a wide range of fields due to their high temporal resolution and has been used for the development of KE-I relationships at several sites worldwide [5,6,13,14,18,[30][31][32][33]. However, despite the technological improvement, disdrometer data are still subject to measurement uncertainty. This uncertainty stems from both instrumental limitations (e.g., splashing, reduced measuring area, internal processing of drops) and environmental influences such as wind. Several studies have tried to quantify these measurement uncertainties and establish possible correction procedures [34][35][36][37][38]. Variations in rainfall measurements between disdrometers and Water 2020, 12, 963 3 of 18 between disdrometers and rain gauges have been observed [35,[39][40][41]. In addition, disdrometers have been shown to measure different drop size distributions even when located at the same site [38,42,43]. The drop size distribution varies naturally in time and space, but disdrometers situated next to each other showed different distributions, which would not be explained by this natural variation [39]. Therefore, the development of a new KE-I relationship based on disdrometer data, is likely to be dependent not only on the geographical rainfall pattern, but also on the type of disdrometer. The measurement uncertainty of disdrometer data presents a need to assess the effect of different disdrometer types on the estimation of rainfall erosivity.
In this study, we aimed to (i) investigate the drop size and velocity distribution measured by three types of disdrometers, (ii) develop new KE-I relationships for five sites based on disdrometer data, (iii) compare the performance of KE-I relationships from the literature on these data, and (iv) validate the site-specific KE-I relationships on rain gauge data, in order to better assess the uncertainties in rainfall erosivity estimation.

Description of Sites and Devices
Measurement of rainfall was carried out at three sites in eastern Austria: Mistelbach (MI), Petzenkirchen (PE) and Rauchenwarth (RA), as well as one site in Prague (PR), Czech Republic and one site in Christchurch (CH), New Zealand. The measurement site locations and mean annual precipitation and temperature can be seen in Table 2 [44 -48]. Figure 1 shows the set-up of the disdrometers at the measurement sites. Rainfall data were collected in 1-minute intervals using optical disdrometers with laser-based sensors, which measure precipitation by determining the size and velocity of raindrops when falling through the measuring area. The rainfall was measured in the period between 2014 to 2019, but with varying amounts of data for each site and disdrometer. The number of rainy minutes used for analysis for each site (after filtering) can be seen in Table 4. The first version of OTT Parsivel (Parsivel) from OTT Hydromet was installed in Petzenkirchen (Parsivel PE) and in Christchurch (Parsivel CH). It has a measuring area of 54 cm 2 and measures intensities from 0.001 to 1200 mm h −1 . Parsivel has 32 drop size and velocity classes, although the two smallest classes are not used as these drop sizes are outside the measurement range of the device.

Filtering of Data
The disdrometers are capable of discerning the type of precipitation into different classes according to the World Meteorological Organizations significant present and forecast weather code (WMO METAR/SPECI w'w' Table 4678) [49]. Based on this classification, only the minutes with rainfall and drizzle or a mixture of the two were selected for further analysis. The Parsivel CH did not output this classification of the precipitation. Thus, data measured during the winter months, where snow could be present, were manually removed before further filtering.
Drops deviating more than 50% from the value of the Atlas et al. [50] diameter-terminal fall velocity relationship were not included in the calculation of KE for establishing the KE-I relationships. A similar filtering was also applied in other studies [38,[51][52][53][54][55][56] to remove erroneously measured drops. Both Parsivel data sets were also cut off for rainfall intensities above 60 mm h −1 to eliminate outliers that were not removed by the filtering.

Development of Site-Specific Rainfall Kinetic Energy-Intensity Relationships
The kinetic energy was determined from the matrix of drop size and velocity classes given by each disdrometer. The rainfall kinetic energy per area per unit time, KE (J m −2 h −1 ), was calculated for each 1-min data retrieval as: where Ni,j is the number of detected raindrops in size class i and velocity class j, A is the measuring area of the disdrometer (m 2 ), Δt (1/60 h) is the interval of data collection, ρ is the density of water (g The first version of OTT Parsivel (Parsivel) from OTT Hydromet was installed in Petzenkirchen (Parsivel PE) and in Christchurch (Parsivel CH). It has a measuring area of 54 cm 2 and measures intensities from 0.001 to 1200 mm h −1 . Parsivel has 32 drop size and velocity classes, although the two smallest classes are not used as these drop sizes are outside the measurement range of the device.

Filtering of Data
The disdrometers are capable of discerning the type of precipitation into different classes according to the World Meteorological Organizations significant present and forecast weather code (WMO METAR/SPECI w'w' Table 4678) [49]. Based on this classification, only the minutes with rainfall and drizzle or a mixture of the two were selected for further analysis. The Parsivel CH did not output this classification of the precipitation. Thus, data measured during the winter months, where snow could be present, were manually removed before further filtering.
Drops deviating more than 50% from the value of the Atlas et al. [50] diameter-terminal fall velocity relationship were not included in the calculation of KE for establishing the KE-I relationships.
A similar filtering was also applied in other studies [38,[51][52][53][54][55][56] to remove erroneously measured drops. Both Parsivel data sets were also cut off for rainfall intensities above 60 mm h −1 to eliminate outliers that were not removed by the filtering.

Development of Site-Specific Rainfall Kinetic Energy-Intensity Relationships
The kinetic energy was determined from the matrix of drop size and velocity classes given by each disdrometer. The rainfall kinetic energy per area per unit time, KE (J m −2 h −1 ), was calculated for each 1-min data retrieval as: where N i,j is the number of detected raindrops in size class i and velocity class j, A is the measuring area of the disdrometer (m 2 ), ∆t (1/60 h) is the interval of data collection, ρ is the density of water (g cm −3 ), D i is the mean drop diameter (mm) of size class i and v j is the mean fall velocity (m s −1 ) of velocity class j.
The rainfall amount, R (mm), and intensity, I (mm h −1 ), were calculated from the drop size distribution by where A is the measuring area of the disdrometer (mm 2 ). For each site and disdrometer, an exponential KE-I relationship was established based on the calculated intensity and KE for each minute. The three-parameter exponential equation (Equation (4)), where a, b and c are empirical constants, has one more parameter than the commonly used linear, logarithmic or power law equations, but it has been recommended as a better fit to the KE-I relationship [17,57]. Linear, logarithmic and power law models suggest that there is no maximum limit to kinetic energy at increasing intensities. However, research has shown that an upper limit may exist [7]. This was also recognized by Wischmeier and Smith [10] as they adapted their equation to remain constant for intensities above 76 mm h −1 .
The parameter a in the exponential model describes an upper limit to the kinetic energy content of rainfall, that is the maximum kinetic energy the rainfall can contain per mm. The exponential model also results in a lower limit of KE at very low intensities, although this may not necessarily be a real phenomenon of rainfall KE [15].
To obtain the best fitting exponential KE-I relationship for each site and disdrometer, the data were fitted with a nonlinear regression by nonlinear least squares [58]. The equation parameters were found by successive approximations to the three-parameter exponential model using the nonlinear least squares function (nls) within the statistical language R [59].
To compare the measured KE from each disdrometer and site to that calculated from the KE-I relationships in Table 1, the disdrometer intensities found from Equation 3 were used in the equations. The logarithmic equation (WS) yielded negative KE values at an intensity below 0.043 mm h −1 . These negative values were not included in the summation of KE.

Validation of Site-Specific Rainfall Kinetic Energy-Intensity Relationships on Rain Gauge Data
The daily rainfall measured by the rain gauge at each site was compared to that measured by the disdrometers for days where both rainfall measurement devices measured more than 1 mm rainfall. A weighing principle rain gauge with a rainfall resolution of 0.01 mm was situated in PE. In MI and PR, the rain gauges were of the tipping bucket type with a rainfall resolution of 0.1 mm. The rain gauge in CH was also a tipping bucket with a rainfall resolution of 0.2 mm. Rain gauge data were not available for the Rauchenwarth site. The number of rainy days compared were 42, 45, 37, 26 and 65 for PWS MI, PWS PE, Thies PR, Parsivel PE and Parsivel CH, respectively.
For validation of the site-specific KE-I relationship established on the disdrometer data, the intensity measured by the rain gauges was used in the site-specific relationship to calculate KE. For the validation of rain gauge values (rain and KE), the unfiltered values from the disdrometers were used, as this would otherwise not correspond to the values measured by the rain gauges.
The percent bias (PBIAS), root mean squared error (RMSE) and the mean absolute error (MAE) were used to quantify the errors in KE estimation between disdrometers and literature KE-I relationships, and in rainfall measurements between disdrometers and rain gauges.

Drop Size and Velocity Distributions
The drop size and velocity distribution obtained by the devices were similar for disdrometers of the same type, despite being placed at different sites ( Figure 2). PWS MI and PWS PE both measured the largest number of drops in the class with a mean diameter of 1.1 mm and mean velocity of 4.4 m s −1 . They also both followed the terminal velocity line well and had few drops in classes deviating more than 50% from this line. Thies RA and Thies PR both registered the highest number of drops at a mean diameter of 0.3 mm and mean velocity of 1.2 m s −1 . Both Thies disdrometers recorded a large number of small drops with high velocities, although it was more pronounced at the RA site. This was also observed by Angulo-Martínez et al. [19,38] and Wilken et al. [28], who speculated it may be due to splashing and break up of drops on the instrument, or a larger occurrence of edge events due to the smaller width of the laser beam. The two Parsivel disdrometers did not record the largest number of drops in the same class. Parsivel PE measured the most drops at a diameter of 0.8 mm and 4.4 m s −1 , while Parsivel CH classified the most drops into the class with a diameter of 0.4 mm and a velocity of 2.6 m s −1 . This may be due to the difference in rainfall regime. However, despite their spatial separation, both Parsivel disdrometers showed similar distributions with an overestimation of the velocity of the majority of drops (up to around a diameter of 1 mm) compared to the terminal velocity line. This distribution was also observed in other studies [41,43,60].
The difference in raindrop distribution between disdrometers arise from differences in the instrument design, measuring area, drop size and velocity classification and internal algorithms. This hinders the direct comparison of rainfall measurements between disdrometers. The device-specific distribution is seen from the comparison of PWS PE and Parsivel PE, which were situated at the same site, but presented different drop size and velocity distributions. This shows that the disdrometer type affects rainfall measurement to such a degree that it makes it difficult to determine the actual geographical rainfall characteristics of the measurement site.
According to the disdrometer manuals, all three disdrometer types have some inbuilt filtering processes, which correct for erroneously measured drops. The amount of information is limited for Thies and Parsivel, where only the effect of possible edge events is mentioned [61,62]. For PWS100 more detail is given e.g., about the so-called unknown drops. These are drops, that for example have been measured as travelling too fast for their size [63]. The unknown drops are removed from the matrix of drop sizes and velocities, which could be a reason for the small number of drops deviating from the terminal fall velocity line. The well-defined drop size and velocity distribution of the PWS100 may suggest that this disdrometer provides better rainfall measurement, especially with regard to kinetic energy.
number of drops in the same class. Parsivel PE measured the most drops at a diameter of 0.8 mm and 4.4 m s −1 , while Parsivel CH classified the most drops into the class with a diameter of 0.4 mm and a velocity of 2.6 m s −1 . This may be due to the difference in rainfall regime. However, despite their spatial separation, both Parsivel disdrometers showed similar distributions with an overestimation of the velocity of the majority of drops (up to around a diameter of 1 mm) compared to the terminal velocity line. This distribution was also observed in other studies [41,43,60].  [50]. Drops in classes outside the red lines deviate more than 50% from the terminal velocity line.  [50]. Drops in classes outside the red lines deviate more than 50% from the terminal velocity line.

Filtering of Drops
The filtering of drops deviating more than 50% from the terminal velocity line by Atlas et al. [50] resulted in a loss of drops. The number of drops filtered out varied between disdrometers and sites. As the PWS MI and PWS PE disdrometers generally follow the terminal velocity line very well, only 3.6% and 4.3% of all drops were removed for PWS MI and PWS PE, respectively. For Thies RA, the filtering excluded 53.6% drops, and for Thies PR it excluded 25.6% drops. Parsivel PE and Parsivel CH had 58.8% and 37.7% of the measured drops removed, respectively. The filtering of drops corresponded to a decrease in measured rainfall amounts of 0.4% for PWS MI and 2.0% for PWS PE. Even though Thies had up to half of the drops filtered out, the small size of the filtered drops resulted in a reduction in rainfall of 12.1% for Thies RA and 5.3% for Thies PR. The Parsivel disdrometers had the highest decrease in measured rainfall due to filtering with 27.3% and 20.3% for Parsivel PE and Parsivel CH, respectively.
Filtering of disdrometer data based on the drop size and velocity deviating from the terminal velocity line has been applied in several other studies [38,[53][54][55][56]. The percentage of removed drops ranged from 12% to 24% in these studies. On the data measured by a Parsivel disdrometer, Jaffrain and Berne [52] employed a filter of 60% of the terminal drop velocity model, which lead to 25% of drops being removed. This corresponded to 3.5% of the total rainfall recorded. The effect of filtering on the mean drop size and velocity can be seen in Table 3. For PWS MI and PWS PE, there was no difference before and after filtering. For Thies and Parsivel, there was a slight increase in mean drop size as all the small, but fast drops had been filtered out. For Thies PR and Parsivel PE the mean velocity increased. The reason for filtering the raindrops is to remove erroneously measured drops, which could potentially affect the estimation of rainfall characteristics measured by the disdrometer, such as KE. The spurious drops can occur due to multiple drops falling at one time, splashing from the instrument or wind effects [51,52]. However, filtering could remove some real raindrops, but in order to reduce the number of outliers, a filtering can be necessary. This was especially the case for the Thies and Parsivel data, as seen in Figure 3, where the scattered nature of the KE-I relationship for each disdrometer before filtering is shown. Thies RA and Parsivel CH especially showed a high degree of scatter. Points of high intensity but low KE were present in these plots, which could indicate the presence of snow. Wind may also have affected the degree of scatter of the KE-I points, as the velocity distribution has been found to broaden under windy conditions due to sub-and super-terminal drops [64]. Johannsen et al. [43] observed the importance of the installation of a wind protection shield for the Thies, as the drop velocity of small drops was seen to increase for events without the wind protection shield installed.

Development of New Rainfall Kinetic Energy-Intensity Relationships
The filtering of the data generally removed outliers and led to less scatter of the plotted KE and I, although for the Parsivel disdrometers some scatter still persists (Figure 3). For PWS MI, PWS PE and Thies PR filtering did not affect the fitting of the KE-I relationships greatly. However, the scattered nature of the data from Thies RA and Parsivel CH made filtering necessary to get a good fit of the model. The filtering of the disdrometer data clearly showed that the parameters of the developed KE-I relationships are sensitive to changes in the data input. Therefore, the amount and quality of the data are important, and any filtering of the data should be clearly stated, as it can affect the outcome of the data analysis to a high degree.
After filtering, the best fitting exponential curve and its equation parameters were found for each disdrometer and site ( Figure 3 and Table 4). PWS MI recorded minute intensities up to 140 mm h −1 and PWS PE recorded intensities up to 180 mm h −1 . This was 2-4 times higher than the other disdrometers, which recorded maximum intensities between 34 mm h −1 and 58 mm h −1 . This may limit the reliability of the developed KE-I relationships on higher intensities. PWS MI and PWS PE also recorded the highest KE values. Water 2020, 12, x FOR PEER REVIEW 9 of 19 The parameter a in the exponential model describes the maximum kinetic energy per mm of rainfall. In the newly developed KE-I relationships, this upper limit to the kinetic energy content of rainfall was the lowest for the Thies disdrometers, whereas it was the highest for the Parsivel disdrometers. van Dijk et al. [15] found an average a parameter of 28.3 (±2.9) from several datasets from the literature. The a parameters in the KE-I relationships, developed on the PWS data, lie within the standard deviation of this average. The Thies a parameter is below this average, while the Parsivel a parameter is above this average. However, other KE-I relationships with higher a parameters have been established on rainfall data from impact and optical disdrometers [8,30,31]. Sanchez-Moreno et al. [13] found an a parameter of 35 with a Parsivel disdrometer in Cape Verde, which corresponds well with the a parameter values found for the Parsivel disdrometers in the present study.  The parameter a in the exponential model describes the maximum kinetic energy per mm of rainfall. In the newly developed KE-I relationships, this upper limit to the kinetic energy content of rainfall was the lowest for the Thies disdrometers, whereas it was the highest for the Parsivel disdrometers. van Dijk et al. [15] found an average a parameter of 28.3 (±2.9) from several datasets from the literature. The a parameters in the KE-I relationships, developed on the PWS data, lie within the standard deviation of this average. The Thies a parameter is below this average, while the Parsivel a parameter is above this average. However, other KE-I relationships with higher a parameters have been established on rainfall data from impact and optical disdrometers [8,30,31]. Sanchez-Moreno et al. [13] found an a parameter of 35 with a Parsivel disdrometer in Cape Verde, which corresponds well with the a parameter values found for the Parsivel disdrometers in the present study.
The KE per unit time (KE=J m −2 h −1 ) as used in the KE-I relationship here is related to the KE per unit rainfall (KE=J m −2 mm −1 ) through the rain intensity [65]. For clarity of comparison of the shape of the curves, the KE has been converted to KE per unit rainfall in Figure 4. It can be seen that the PWS MI curve reaches its maximum KE at lower intensities than the other equations. This behavior is determined by the c parameter of the equation, which defines the general shape of the curve [15]. As PWS MI has a higher value of c the curve reaches its max KE faster. The lowest value of c is found in the relationship for Parsivel CH, and it is clearly seen that the curve for this relationship only gradually reaches its max KE value at higher intensities.

Comparison with Rainfall Kinetic Energy-Intensity Relationships from the Literature
The KE sum of all analyzed rainfall minutes measured by the disdrometer, was compared to that found by inserting the disdrometer intensities into the four literature KE-I relationships ( Table  1). The results presented in Table 5 show that there was no one KE-I relationship that fit the best at all sites. For PWS MI, the lowest percent bias of the four literature KE-I relationships was from WS, which underestimated KE from the disdrometer by 8.3%. For PWS PE the lowest PBIAS resulted from

Comparison with Rainfall Kinetic Energy-Intensity Relationships from the Literature
The KE sum of all analyzed rainfall minutes measured by the disdrometer, was compared to that found by inserting the disdrometer intensities into the four literature KE-I relationships ( Table 1). The results presented in Table 5 show that there was no one KE-I relationship that fit the best at all sites. For PWS MI, the lowest percent bias of the four literature KE-I relationships was from WS, which underestimated KE from the disdrometer by 8.3%. For PWS PE the lowest PBIAS resulted from the VD equation with −6.6%. For Thies RA the WS had the lowest PBIAS of 5.8%. The BF equation gave the nearest fit to the Thies PR data with a PBIAS of −2.1%. Parsivel PE agreed the most with VD with a PBIAS of −4.6%, while Parsivel CH had the lowest PBIAS of 2.3% by MG. The site-specific relationships were calibrated on these data, so the error parameters show how well the calibration of the relationship was. Table 5. Comparison of the measured rainfall kinetic energy to that found from the calibration of the site-and device-specific rainfall kinetic energy-intensity relationships and the rainfall kinetic energy-intensity relationships Wischmeier   The only consistent trend was that BF always underestimated the total KE with PBIAS ranging from −2.1% to −30.2%. Several other studies based on disdrometer data also found that BF underestimates KE [5,8,18,32,33,66]. Nearing et al. [22] stated that BF should not be used for erosivity calculations, as it underestimates KE, and recommended using WS or MG, as they give the best results. This was also found in the present study, where together with VD these relationships gave good results for most datasets. Contrary to all other sites, BF was the best fitting equation to the Thies PR data measured on all error parameters. From Figure 4 it can be seen that the data points for Thies PR generally cluster at a lower KE value than the other sites, which causes a better fit with the BF curve. Angulo-Martínez et al. [19] collected data from a Thies disdrometer in Spain and also found that Brown and Foster performed quite well at estimating total KE for their data.
Petan et al. [5] measured KE at two sites in Slovenia, one equipped with a Parsivel and the other with a Thies disdrometer. They tested the agreement of the disdrometer KE with that found by the WS, BF and VD relationships and found similar percentage deviations for each relationship for both disdrometer types. However, they calculated drop velocity after Atlas et al. [50], so velocity differences between the devices had been eliminated.

Validation of Site-Specific Rainfall Kinetic Energy-Intensity Relationships on Rain Gauge Data
The new site-specific KE-I relationships developed from the PWS MI, PWS PE, Thies PR, Parsivel PE and Parsivel CH were tested by applying them to the intensities measured by rain gauges at the measurement sites and calculating daily KE and comparing it to the measured KE from the disdrometers ( Figure 5). For all sites but CH, the rain gauges generally measured higher total daily rainfall amounts and KE than the disdrometers. The overall lowest deviation between disdrometer and rain gauge was found for PWS PE with a PBIAS of 6.2%. The rain gauge in PE was of the weighing type and these have been found to have higher accuracy than tipping buckets [67,68]. The highest deviation in daily rainfall amount was observed between Parsivel PE and the same rain gauge with a PBIAS of 39.7%, which suggests degraded performance of the disdrometer.
The lowest overall deviation of KE measured by a disdrometer and KE estimated with rain gauge intensity from the developed site-specific KE-I relationship was found for Parsivel CH. The underestimation of total daily rainfall amount also resulted in an underestimation of total daily KE with a PBIAS of −11.7%. The comparison between Parsivel PE and the rain gauge gave the highest overall deviation of daily KE with an overestimation of 104.2% by the rain gauge. This was expected due to the low agreement in rainfall amount measurements between Parsivel PE and the rain gauge.
Even though PWS PE showed the best agreement with the rain gauge, the KE was overestimated by 28.5% by the rain gauge. This may indicate that the developed KE-I relationship is prone to overestimation of KE. It may also suggest that rainfall intensity measurement between the PWS PE and the rain gauge may deviate more within events than can be seen from the daily comparison. This was shown by Johannsen et al. [43], who found that the rain gauge often measured higher intensities than four collocated disdrometers at the minute-scale.
Fraile et al. [37] observed lower rainfall amounts from a disdrometer than rain gauges and attributed this to the use of a standard measuring area. By using the effective measuring area, which decreases as the drop size increases, they found, that the rainfall amount measured by the disdrometer increased and generally became closer to that of the rain gauge. Other reasons, such as minimum rainfall resolution, smaller sampling area of the disdrometers compared to the rain gauges [39], exclusion of the smallest drop classes (Parsivel, [61]) or difference in drop size classification [69,70] may also affect the rainfall measurements. In addition, insufficient calibration of the tipping buckets may also have played a role, especially at the MI site, where an older model, which seems to have overestimated rainfall, was installed. However, these effects have not been explicitly investigated in the present study.
From these results, it is seen that the validation of the site-specific KE-I relationships using rain gauge data is only as good as the agreement between the rain gauge and disdrometer rainfall measurement. Depending on which rainfall measurement method is used, the KE estimates will vary due to differences in measured rainfall intensity and amount.

Implications of Disdrometer Differences for the Estimation of Rainfall Erosivity
The implication of the disdrometers measuring rainfall differently is that results are devicespecific and therefore, the direct comparison of results from different disdrometer types is hindered. This can clearly be seen by the comparison of PWS PE and Parsivel PE, which were situated at the same site, but gave distinctive drop size and velocity distributions, which resulted in different KE-I relationships. This indicates that using one device instead of another may have a greater effect on the outcome of the developed KE-I relationship, than the local rainfall regime. For the formulation of the KE-I relationship, the PWS PE and Parsivel PE data did not cover exactly the same time period. Thus, a difference in measured rainfall events may explain some of the differences between the two devices. However, in Johannsen et al. [43], the KE-I relationships of the PWS100 and Parsivel disdrometers were established for the same rainfall events and were shown to differ due to instrumental differences.
In the present study, the actual measurement difference between the two disdrometers was quantified, by comparing the daily rainfall amount and KE for 27 rainy days, where both disdrometers measured more than 1 mm rainfall. The total rainfall amount measured by Parsivel PE

Implications of Disdrometer Differences for the Estimation of Rainfall Erosivity
The implication of the disdrometers measuring rainfall differently is that results are device-specific and therefore, the direct comparison of results from different disdrometer types is hindered. This can clearly be seen by the comparison of PWS PE and Parsivel PE, which were situated at the same site, but gave distinctive drop size and velocity distributions, which resulted in different KE-I relationships. This indicates that using one device instead of another may have a greater effect on the outcome of the developed KE-I relationship, than the local rainfall regime. For the formulation of the KE-I relationship, the PWS PE and Parsivel PE data did not cover exactly the same time period. Thus, a difference in measured rainfall events may explain some of the differences between the two devices.
However, in Johannsen et al. [43], the KE-I relationships of the PWS100 and Parsivel disdrometers were established for the same rainfall events and were shown to differ due to instrumental differences.
In the present study, the actual measurement difference between the two disdrometers was quantified, by comparing the daily rainfall amount and KE for 27 rainy days, where both disdrometers measured more than 1 mm rainfall. The total rainfall amount measured by Parsivel PE was 26.0% lower than that measured by PWS PE and the MAE in rainfall amount was 2.7 mm. The total accumulated daily KE was also underestimated by Parsivel PE compared to PWS PE with 37.9% and a MAE of 65.3 J m −2 . This deviation in measured rainfall amount and KE between devices placed at the same site shows the impact of the use of different disdrometers on the estimation of rainfall erosivity.
The sites Mistelbach and Rauchenwarth are within the same region of Austria and have similar average annual precipitation, and as such should not differ much in rainfall characteristics. However, due to the use of different devices, the drop size and velocity distributions, mean drop size and velocity and KE-I relationships found for these two sites were highly different. This shows the impact of using different devices under similar rainfall regimes.
On the other hand, data from the same type of disdrometer should be comparable, also when situated at different sites. As such, the results of PWS MI and PWS PE could enable comparison of the spatial differences in rainfall characteristics at the Mistelbach and Petzenkirchen sites. The comparison showed that mean drop size and velocity is larger in MI than in PE ( Table 2). The maximum KE per rainfall unit is slightly lower for PWS MI, but the rise of the curve of the KE-I relationship reaches its maximum KE at lower intensities than that of PWS PE. This means, that at lower intensities (until the intersect of the curves), the KE calculated from the site-specific KE-I relationships is higher in Mistelbach than in Petzenkirchen.
The difference in raindrop measurement by the different disdrometers is related to the design and internal processing of each disdrometer. Variations in instrument measuring area, rainfall resolution, number of classes of drop size and velocity and built-in correction processes of drop measurement all affect the final rainfall measurements. The device-specific differences could be corrected for by calculating drop fall velocity from a terminal fall velocity equation [5,32,71], applying filters or correction schemes to the drop size distribution [38,41] or by applying a correction factor for KE between disdrometers [43]. However, this will not make results fully comparable and might introduce other biases. Thus, to have comparable results from multiple measurement sites, the same type of disdrometer should be used.

Conclusions
This study demonstrated that the rainfall erosivity estimation is influenced by the use of different types of rainfall measurement devices. Depending on the disdrometer type, rainfall was classified in a certain way, leading to distinctive drop size and velocity distributions for each disdrometer. Disdrometers of the same type, however, showed similarities in drop size and velocity distributions across the measurement sites. The drop sizes and velocities measured by the PWS100 disdrometers generally followed the terminal velocity line well. A large number of small drops with high velocities was measured by both Thies disdrometers and an overestimation of the velocity of smaller drops was measured by the two Parsivel disdrometers. The broad drop size and velocity distribution led to a scattered nature of the KE-I data points for the Thies and Parsivel disdrometers, which required a filtering of the data in order to develop new KE-I relationships.
Exponential KE-I relationships were developed for each site and disdrometer based on the filtered disdrometer data. The found parameters of the exponential relationship were different for each type of disdrometer and site.
The best fit of KE-I relationships from the literature (WS, BF, MG, VD) varied among disdrometers and sites. One universal KE-I relationship for all the sites could not be recommended. The BF relationship consistently underestimated KE compared to that measured by the disdrometers.
Comparison with the literature KE-I relationships showed the importance of the site-specific disdrometer calibration as the percent bias of the relationships from literature was generally higher.
The site-specific KE-I relationships were successfully validated on rain gauge data, even though when rain gauge and disdrometer for any reason recorded different rainfall amount, the bias in KE estimation increased.
Results showed that, rainfall measurements by multiple types of disdrometers were associated with uncertainties stemming from instrumental differences. This led to variations in rainfall erosivity estimations between disdrometers even when placed at the same site. This was seen by the large percent bias in daily rainfall and KE measured by the PWS100 and Parsivel placed at the same site. Therefore, various disdrometers should not be used for comparison of the rainfall characteristics of different sites, as this will hinder the determination of spatial differences. If the use of the same type of disdrometer is not possible, or when comparing results from several studies, the possible effect of the specific disdrometer on the results must be recognized.
From our results, we cannot conclude which disdrometer type is better for rainfall measurements. Further research with methods such as artificial rainfall simulators with well-defined rainfall characteristics or the implementation of e.g., video disdrometers could help compare disdrometers to a certain reference and eliminate some of the uncertainties in rainfall erosivity estimation between instruments.