Assessment of Various Empirical Soil Loss Estimation Equations in Arid Regions ()
1. Introduction
Soil in most regions of the world is playing an important role in supporting the natural ecosystems (Singer & Warkentin, 1996). The main medium for plant growth is given by soil, which provides water, nutrients, and air. Many countries suffer from soil erosion which is considered a significant problem for some of them (Lal, 2001), particularly in developing countries throughout the world that rely on agriculture (Morgan et al., 1998). Soil erosion is also considered the source of the world’s economic and environmental problems (Lal, 1998). A decline in soil fertility and quality, loss of topsoil, soil pollution, sedimentation and pollution of water bodies, and an overall rise in soil degradation are some of the many issues associated with soil loss and erosion (Martin-Moreno et al., 2008).
Soil erosion is defined as a physical process that detaches soil particles from the surface of the ground and removes them due to wind, rainfall, water flow, tillage or irrigation, and other forces (Segura et al., 2014; Brady & Weil, 2010). The principal type of erosion is considered to be erosion by water. Soil erosion may not be visible on the ground surface when any amount of particles has been transported by raindrops, but soil erosion would be more evident when water forms rills and gullies flow concentrated (Kim & Julien, 2006).
The Revised Universal Soil Loss Equation (RUSLE) (Renard et al., 1997) is the most comprehensive work for estimating the soil loss. RUSLE is an empirical model and is one of the most common methods for soil erosion estimation (Renard et al., 2011). Thanks to its simplicity and the large database, it is a very popular model (Charlton, 2007). One of the most important factors in the RUSLE equation is the rainfall erosivity R-factor which is the focus of this study. The RUSLE was enhanced by revising the climatic factor, the seasonal soil erodibility factor, updating the gradient and slope length, and inventing a new technique for calculating the cover management factor (Renard et al., 1997). The RUSLE assumes that the sediment content of the flow controls detachment and deposition (Pitt, 2007). Erosion is limited by the flow’s carrying capacity but is not limited by the source (Pitt, 2007). When the sediment load has reached the carrying capacity of the flow, detachment can no longer take place (Pitt, 2007). A full historical review of the R-factor development is given in Nearing et al. (2017).
Numerous erosion methods were used to estimate the R-Factor and the resulting soil loss over the last few decades. Yin et al. (2017) and Benavidez et al. (2018) reviewed methodologies to estimate the R-Factor, some of the different worldwide studies that developed rainfall erosivity equations, the locations in which the study was conducted, as well as the other studies that used these equations. Numerous studies focused on Europe with a special concern on impact of climate change on rainfall erosivity (Panagos et al. 2015, 2016, 2017a). Moreover, several studies estimated the R-factor to determine the effect of rainfall on soil in many sites in the Arabian Peninsula, which is study area. In Saudi Arabia, Mallick et al. (2014) used the empirical equation developed by Arnoldus (1980) to estimate the soil loss in Abha mountainous watersheds in Assir Province, Saudi Arabia. Bahrawi et al. (2016) assessed the soil erosion in Wadi Yalamlam watershed to the Southeast of Jeddah city, Saudi Arabia, using the Barfield et al. (1981) equation with annual precipitation as input. Based on the equation developed by Fournier (1960) (which is the previous form of Arnoldus (1980) equation), Azaiez et al. (2020) provided the average soil loss due to rainfall in Wadi Mirabah, located west of the city of Abha. To the northeast of the study area, in Jordan, the equation developed by Eltaif et al. (2010) was used by Farhan et al. (2013) and Farhan & Nawaiseh (2015) to estimate soil erosion in Wadi Kufranja watershed in the northern highlands of Jordan and in Wadi Kerak watershed is in the southeast of the Dead Sea, Jordan, respectively. All used equations were originally developed to evaluate the R-factor equation outside the Arabian Peninsula and even not in the arid and hyper-arid regions.
Unfortunately, none of the studies mentioned above calibrated the equation used to calculate the R-factor in their study area. On the other hand, Panagos et al. (2017b) produced a global map of the R-factor for the entire world, where the stations located in our study region were only 15 stations in Kuwait, which represents, in terms of area, 5.27% of the Arabian Peninsula. The distribution of the R-factor, which is shown in Panagos et al. (2017b) study through the Arabian Peninsula, is thus mainly produced by interpolation. Therefore, the soil loss equations applied in the Arabian Peninsula were not calibrated, and the current research tries to fill this scientific gap which is considered the main objective for this paper.
2. Study Area Characteristics and Data Sources
2.1. Study Area
The study area is the Arabian Peninsula, located from 12˚N to 32˚N and from 36˚E to 61˚E. The Arabian Peninsula is one of the biggest deserts in the world. The area of the Arabian Peninsula is around 3.2 million km2, approximately the size of the Indian subcontinent. It is located in Southwest Asia and joins Asia with Africa. It is bordered by the Levant area (a region known as the Fertile Crescent) to the north, the Red Sea to the west and southwest, the Arabian Gulf to the northeast, and the Indian Ocean to the southeast. The length of the Red Sea boundary is approximately 1900 kilometers and the maximum width is 2500 kilometers from Yemen to Oman (Nijim et al., 2019).
2.2. Data Sources
The available sub-daily rainfall data obtained from the rainfall stations in the study area are obtained from 104 stations. Exactly 49 stations are available in Oman, one (1) station in Dubai and 54 stations in Saudi Arabia. The 104 stations of Oman, Dubai and Saudi Arabia provide sub-daily rainfall data. The study area and the used rainfall stations with available sub-daily rainfall data are shown in Figure 1.
3. Methodology
3.1. Evaluation of the EI30 Using Sub-Daily-Duration Rainfall Data
The RUSLE model equation is given below:
Figure 1. Study area and rainfall stations with available sub-daily data.
(1)
where:
A—Annual soil loss, in tons·ha−1·year−1·
R—Rainfall erosivity factor in MJ·mm/(ha·hr·year).
K—Soil erodibility factor, the erosion rate for a specific soil in continuous fallow condition on a 9% slope having a length of 22.1 m in ton·ha·hr/(MJ·mm·ha).
LS—Topographic factor that reflects the length of the slope and steepness of the slope. It is the soil loss ratio of a particular site to that from a unit site with the same soil and slope but with a length of 22.1 m.
C—Cover management factor, which is the protective covering of the canopy and organic material in direct contact with the soil. It is calculated as the ratio of soil loss from land cultivated under particular circumstances to the corresponding loss from tilled land under clean-tilled continuous fallow conditions (Renard et al., 1997).
PS—Support practice factor that reflects the activities of soil protection or other erosion control measures. It is measured as the soil loss ratio with a particular support practice to the corresponding loss with up and downslope plowing (Renard et al., 1997).
In general, the rainfall data have been recorded in the rainfall stations in the form of the depth of rainfall versus its duration. EI30 was evaluated according to Wischmeier & Smith (1978) method. The R-factor is obtained by applying in Equation (2).
(2)
where:
R—Rainfall erosivity factor;
E—total storm kinetic energy (MJ/ha);
I30—maximum 30 minutes rainfall intensity;
j—index for the number of years used to calculate the average;
k—index of the number of storms every year;
n—number of years to find the average;
m—number of storms every year.
3.2. Investigation of the Published Equations
Previously published equations estimating the R factor based on available rainfall data such as annual rainfall data (P) were investigated to determine their applicability in the Arabian Peninsula by comparing the outcome of the various equations against the calculated R-factor using sub-daily rainfall data. Table 1 shows the forms and parameters of these equations, which vary between fractional equation, linear equation, power equation and exponential equation as shown in Table 1.
3.3. Performance Criteria
The performance criteria were used to find out how similar the values and model simulations observed are. To assess model performance, correlation-based methods are widely used. To compare models, these mathematical criteria are used to calculate some kind of distance between the simulated values and the reference values (Waseem et al., 2017). There are several performance criteria to evaluate models. Three of these criteria are used in this study:
Table 1. Published equations using commonly available rainfall data as input.
1) the Root Mean Square Error (RMSE),
2) the Mean Absolute Percent Error (MAPE), and
3) the Mean Absolute Error (MAE).
The formulas to calculate these criteria are given hereafter in Equations (3) to (5).
(3)
(4)
(5)
where
is the “reference” rainfall erosivity value calculated from the sub-daily information;
is the estimated (predicted) value using any of the empirical equations;
is the average value for the exact series; and
N is the number of data points.
4. Results and Discussions
4.1. Estimating Rainfall Erosivity for the Available Stations with Sub-Daily Duration Data
Using the procedure of Wischmeier & Smith (1978) Equation (2) calculating R-factor per storm, adding them and then averaging for each year, the values of R-factor were obtained. Table 2 shows the R-factor calculated for the 104 stations in the study area. R-factor ranges from 53.96 (Mj·mm)/(ha·hr) to 2986.43 (Mj·mm)/(ha·hr) with a mean of 680.62 (Mj·mm)/(ha·hr) and a standard deviation of 576.47 (Mj·mm)/(ha·hr).
4.2. Assessment of the R-Factor Empirical Equations
Various equations estimate the R-factor based on the mean annual rainfall as input. In this study, seven (7) previously published equations were assessed against the exact R-factor estimated from sub-daily rainfall data. The scatter plots between the exact (reference) R-factor (which is the value of the R-factor obtained from sub-daily rainfall data) and the calculated R-factor by previously published equations are illustrated superimposed distributed on the 45˚ (1:1) line (Figure 2). The 1:1 line allows identifying how the reference values compare with the estimated ones by published equations. All tested equations except Roose (1976) significantly underestimated the R-factor. In fact, equations like Mikhailova et al. (1997) and Torri et al. (2006) even produced negative R-factor values. On the contrary, Roose (1976) equation produced overestimated values. Furthermore, the obtained RMSE and MAE values were almost as high as the average R-factor, with MAPE exceeding 100% (Table 3). Based on the above, the
Table 2. Reference rainfall erosivity values for the 104 stations.
Table 3. Results of assessment of published equations to evaluate the R-factor.
equation that gave the best results (although unsatisfactory) was that of Roose (1976), which is a simple linear equation. We tested a recalibration of this simple form and the obtained calibrated equation (without intercept as it was not statistically significant from zero) gave an RMSE of 280 (Mj·mm/(ha·hr)) and the MAPE dropped to 47.6% (Table 3).
5. Conclusion
Sub-daily rainfall data in 104 rainfall stations in the study area of the Arabian Peninsula are used to calculate the Erosivity R-Factor of the RUSLE equation. Using the procedure developed by the Soil Conservation Service in the late seventies, the values of the R-factor were calculated for the 104 stations and were used as reference values. The R-factor ranged from 53.96 (Mj·mm)/(ha·hr) to 2986.43 (Mj·mm)/(ha·hr) with a mean of 680.62 (Mj·mm)/(ha·hr) and a standard deviation of 576.47 (Mj·mm)/(ha·hr).
The results of seven empirical equations previously used (without calibration) in the study area were tested against the exact (reference) R-factor. All equations proved not suitable for application in the Arabian Peninsula, as an example of an arid region. All tested equations except one significantly underestimated the R-factor. In fact, some of the used equations produced negative R-factor values. Furthermore, the obtained RMSE and MAE values were almost as high as the average R-factor, with MAPE exceeding 100%. Therefore, we advise against applying these equations in the Arabian Peninsula because of their poor results. All of them require parameter recalibration. Recalibrating a simple linear regression equation gave an RMSE of 280 (Mj·mm/(ha·hr)) and the MAPE dropped to 47.6%, when compared to the reference R-factor values.
The current research opens a new starting point for extensive research on the different methods to predict the rainfall erosivity factor within areas with rainfall stations where only daily, monthly, or yearly precipitations are available, or regions where rainfall data is scarce. Other explanatory factors (such as elevation, wind or lee sides …) might be investigated to see if they can improve the outcomes. We also recommend investigating the daily rainfall, or some kind of characteristics derived from it, as an explanatory factor to estimate the R-factor in arid region.
Acknowledgements
The authors would like to thank the Ministry of Regional Municipalities and Water Resources (MRMWR), Oman, and the Ministry of Environment, Water and Agriculture (MEWA), Saudi Arabia, for providing the raw rainfall data for this study.