Spatial distribution and temporal trends of rainfall and erosivity in the Eastern Africa region

Soil erosion by water is one of the main environmental concerns in the drought‐prone Eastern Africa region. Understanding factors such as rainfall and erosivity is therefore of utmost importance for soil erosion risk assessment and soil and water conservation planning. In this study, we evaluated the spatial distribution and temporal trends of rainfall and erosivity for the Eastern Africa region during the period 1981–2016. The precipitation concentration index, seasonality index, and modified Fournier index have been analysed using 5 × 5‐km resolution multisource rainfall product (Climate Hazards Group InfraRed Precipitation with Stations). The mean annual rainfall of the region was 810 mm ranging from less than 300 mm in the lowland areas to over 1,200 mm in the highlands being influenced by orography of the Eastern Africa region. The precipitation concentration index and seasonality index revealed a spatial pattern of rainfall seasonality dependent on latitude, with a more pronounced seasonality as we go far from the equator. The modified Fournier index showed high spatial variability with about 55% of the region subject to high to very high rainfall erosivity. The mean annual R‐factor in the study region was calculated at 3,246 ± 1,895 MJ mm ha^?1 h^?1 yr^?1, implying a potentially high water erosion risk in the region. Moreover, both increasing and decreasing trends of annual rainfall and erosivity were observed but spatial variability of these trends was high. This study offers useful information for better soil erosion prediction as well as can support policy development to achieve sustainable regional environmental planning and management of soil and water resources.     

Drop size distribution and kinetic energy load of rainfall events in the highlands of the Central Rift Valley,Ethiopia

Abstract

Knowledge of rainfall characteristics is important for estimating soil erosion in arid areas. We determined basic rainfall characteristics (raindrop size distribution, intensity and kinetic energy), evaluated the erosivity of rainfall events, and established a relationship between rainfall intensity I and volume-specific kinetic energy KE_vol for the Central Rift Valley area of the Ethiopian highlands. We collected raindrops on dyed filter paper and calculated KE_vol and erosivity values for each rainfall event. For most rainfall intensities the median volume drop diameter (D₅₀) was higher than expected, or reported in most studies. Rainfall intensity in the region was not high, with 8% of rain events exceeding 30 mm h^-1. We calculated soil erosion from storm energy and maximum 30-min intensity for soils of different erodibility under conditions of fallow (unprotected soil), steep slope (about 9%) and no cover and management practice on the surface, and determined that 3 MJ mm ha^-1 h^-1 is the threshold erosivity, while erosivity of >7 MJ mm ha^-1 h^-1 could cause substantial erosion in all soil types in the area.

Editor Z.W. Kundzewicz; Associate Editor Q. Zhang     

Assessing changes in rainfall erosivity in Sicily during the twentieth century

Changes in rainfall erosivity are an expected consequence of climate change. Long‐term series of the single storm erosion index, EI, may be analysed to detect trends in rainfall erosivity. An indirect approach has to be applied for estimating EI, given that long series of rainfall intensities are seldom available. In this paper, a method for estimating EI from the corresponding rainfall amount, h_e, was developed for Sicily. This method was then applied at 17 Sicilian locations, representative of different climatic zones of the region, to generate a long series (i.e. from 1916 to 1999 in most cases) of EI values. Linear and step (step located at 1970) trends in annual and seasonal erosivity were detected by both classical approaches (Mann–Kendall test, Wilcoxon‐Mann‐Whitney rank‐sum test) and a new empirical approach (quantile approach, QA), based on the determination of the erosivity values corresponding to selected probability levels. A power relationship between EI and h_e with a space‐ and time‐variable scale factor and a time‐variable process parameter yielded the most accurate predictions of EI. However, a simpler model, using a time‐variable scale factor and a constant process parameter, yielded reasonably accurate EI estimates. Annual erosivity did not increase in Sicily during the twentieth century. At the most, it decreased at a few locations (three of the 17 considered locations). Significant trends were observed more frequently for winter erosivity (six locations) than for summer erosivity (two locations), suggesting that the erosive storms of winter determined the occasional occurrence of a negative trend in annual erosivity. In general, the QA compared reasonably well with more classical approaches. The QA appears promising since step trends for different return periods may be detected but efforts are needed to statistically formalize the proposed approach. Copyright © 2007 John Wiley & Sons, Ltd.     

Spatial and temporal variations in rainfall erosivity during 1960–2005 in the Yangtze River basin

Water resources and soil erosion are the most important environmental concerns in the Yangtze River basin, where soil erosion and sediment yield are closely related to rainfall erosivity. The present study explores the spatial and temporal changing patterns of the rainfall erosivity in the Yangtze River basin of China during 1960–2005 at annual, seasonal and monthly scales. The Mann–Kendall test is employed to detect the trends during 1960–2005, and the T test is applied to investigate possible changes between 1991–2005 and 1960–1990. Meanwhile the Rescaled Range Analysis is used for exploring future trend of rainfall erosivity. Moreover the continuous wavelet transform technique is using studying the periodicity of the rainfall erosivity. The results show that: (1) The Yangtze River basin is an area characterized by uneven spatial distribution of rainfall erosivity in China, with the annual average rainfall erosivity range from 131.21 to 16842 MJ mm ha^?1 h^?1. (2) Although the directions of trends in annual rainfall erosivity at most stations are upward, only 22 stations have significant trends at the 90 % confidence level, and these stations are mainly located in the Jinshajiang River basin and Boyang Lake basin. Winter and summer are the seasons showing strong upward trends. For the monthly series, significant increasing trends are mainly found during January, June and July. (3) Generally speaking, the results detected by the T test are quite consistent with those detected by the Mann–Kendall test. (4) The rainfall erosivity of Yangtze River basin during winter and summer will maintain a detected significant increasing trend in the near future, which may bring greater risks to soil erosion. (5) The annual and seasonal erosivity of Yangtze River basin all have one significant periodicity of 2–4 years.     

Spatial soil loss risk and reservoir siltation in semi-arid Tunisia

Abstract

Soil erosion vulnerability and extreme rainfall characteristics over the Mediterranean semi-arid region of Tunisia are crucial input for estimation of siltation rate in artificial reservoirs. A comprehensive high-resolution database on erosive rainfall, together with siltation records for 28 small reservoirs, were analysed for this region, the Tunisian Dorsal (the easternmost part of the Atlas Mountains). The general life-span of these reservoirs is only about 14 years. Depending on the soil degradation in the different catchments, the corresponding reservoirs display a wide range of soil erosion rates. The average soil loss was 14.5 t ha^?1 year^?1 but some catchments display values of up to 36.4 t ha^?1 year^?1. The maximum 15-min duration rainfall intensity was used to determine the spatial distribution of rainfall erosivity. The northwestern parts of the Tunisian Dorsal display the most extreme rainfall erosivity. Spatial erosion patterns are to some extent similar; however, they vary greatly according to their location in the “soil degradation cycle”. This cycle determines the soil particle delivery potential of the catchment. In general, the northwestern parts of the Dorsal display modest soil erosion patterns due to the already severely degraded soil structure. Here, the soil surface is often the original bedrock. However, the greatest soil erosion occurs in the mid-eastern parts of the Dorsal, which represents the “degradation front”. The latter corresponds to the area with highest erosion, which is continuously progressing westward in the Dorsal. The large variation between the erosive rainfall events and the annual soil loss rates was explained by two important factors. The first relates to the soil degradation cycle. The second factor corresponds to the degradation front with the highest soil loss rates. At present this front is located at 300 m altitude and appears to be moving along an 80-km westward path starting from the east coast. A better understanding of the above can be used to better manage soils and soil covers in the Tunisian Dorsal area and, eventually, to decrease the soil erosion and reservoir siltation risk.

Citation Jebari, S., Berndtsson, R., Bahri, A. & Boufaroua, M. (2010) Spatial soil loss risk and reservoir siltation in semi-arid Tunisia. Hydrol. Sci. J. 55(1), 121–137.     

EROSION RATES ON UPLAND AREAS IN THE MIDSOUTH USA

lINTR0DUCTIONSoilerosi0n0fslopingagriculturallandoccursduringandimmediatelyafterrainfalleventswherethesoilsurfaceisn0tadequatelyprotectedfromtheerosiveforces0frainfallandrunoffThebestaltemativefortheconservationplarinerist0minimizesoilerosionasmuchasp0ssible.Managementofsoilsurfaceconditionsinvolvesexperienceandjudgmeniinselectingmanagementpracticesthatutilizecr0presiduesandminimizedisturbancebytillage.Concentfatedrun0ffs0metimescausessevereerosi0n,particularlyonunprotectedsteepslopes.Ril…     

Estimating rainfall erosivity by aggregated drop size distributions

Rainfall erosivity is defined as the potential of the rain to cause erosion, and it can be represented by rainfall kinetic power. At first in this paper, the raindrop size distributions (DSD) measured by an optical disdrometer located at Palermo in the period June 2006–March 2014 and aggregated for intensity classes, are presented. Then an analysis of raindrop size characteristics is carried out, and the reliability of Ulbrich's distribution, using both the maximum likelihood and momentum estimate parameter methods, is tested. The raindrop size measurements are used to determine the experimental rainfall kinetic power values, which are compared with the ones calculated by a theoretically deduced relationship. This analysis demonstrates that the kinetic power is strictly related to the median volume diameter of DSD. Finally, the reliability of the simplest Marshall and Palmer exponential DSD for estimating the rainfall kinetic power is demonstrated. Copyright © 2015 John Wiley & Sons, Ltd.     

Reliability of rainfall kinetic power–intensity relationships

The rainfall erosivity plays a fundamental role in water soil erosion processes and it can be expressed by its kinetic power. At first in this paper, the raindrop‐size distributions measured, in the period June 2006–March 2014, by an optical disdrometer installed at the Department of Agricultural and Forestry Sciences of University of Palermo are aggregated into rainfall intensity classes, having different ranges, and the measured kinetic power values are determined. Measured kinetic power values are initially used for testing the applicability of the kinetic power‐rainfall intensity relationships proposed by Wischmeier and Smith ( 1978 ), used in Universal Soil Loss Equation (USLE), Brown and Foster ( 1987 ) (RUSLE), and McGregor et al. ( 1995 ) (RUSLE2). Then, the reliability of a theoretical relationship for estimating the kinetic power by rainfall intensity and median volume diameter is verified. Finally, using the literature available datasets, corresponding to measurements carried out by different techniques and in different geographical sites, the analysis demonstrated that the rainfall intensity is not sufficient to determine the rainfall kinetic power. On the contrary, the theoretically deduced relationship allows to reproduce adequately the kinetic power of all available datasets, demonstrating that the knowledge of both rainfall intensity and median volume diameter allows a reliable estimate of the rainfall erosivity.     

Drop size distribution and kinetic energy load of rainstorms in Hong Kong

Data on drop size distribution and kinetic energy load of rainstorms are basic for rainfall erosivity indices. A simple and relatively inexpensive instrument was used to asses the instantaneous intensity and kinetic energy load of rainstorms in Hong Kong. Both the drop size and the instantaneous kinetic energy load of rainfall in Hong Kong are greater than in temperate and subtropical climates. The high kinetic energy results from the large size and greater number of raindrops falling per unit time. A high correlation between the kinetic energy of rainfall and the amount of rainfall allows for a convenient estimate of the energy load of storms from the amount of rainfall. Of more significance to the erosion process is the determination that about 74% of the total annual rainfall is erosive, containing about three‐quarters of the total annual energy load of the rains. The variability of rainfall parameters within a rainfall and from storm to storm is shown. The energy–intensity relationship, seasonal and annual distributions of rainfall erosivity are presented. Copyright © 2000 John Wiley & Sons, Ltd.     

Use of remote sensing,GIS and C++ for soil erosion assessment in the Shakkar River basin,India

Soil is an essential resource for human livelihoods. Soil erosion is now a global environmental crisis that threatens the natural environment and agriculture. This study aimed to assess the annual rate of soil erosion using distributed information for topography, land use and soil, with a remote sensing (RS) and geographical information system (GIS) approach and comparison of simulated with observed sediment loss. The Shakkar River basin, situated in the Narsinghpur and Chhindwara districts of Madhya Pradesh, India, was selected for this study. The universal soil loss equation (USLE) with RS and GIS was used to predict the spatial distribution of soil erosion occurring in the study area on a grid-cell basis. Thematic maps of rainfall erosivity factor (R), soil erodibility factor (K), topographic factor (LS), crop/cover management factor (C), and conservation/support practice factor (P) were prepared using annual rainfall data, soil map, digital elevation model (DEM) and an executable C++ program, and a satellite image of the study area in the GIS environment. The annual rate of soil erosion was estimated for a 15-year period (1992–2006) and was found to vary between 6.45 and 13.74 t ha^?1 year^?1, with an average annual rate of 9.84 t ha^?1 year^?1. The percentage deviation between simulated and observed values varies between 2.68% and 18.73%, with a coefficient of determination (R²) of 0.874.     

Exploring the effect of different time resolutions to calculate the rainfall erosivity factor R in Calabria,southern Italy

The rainfall erosivity factor R of the Universal Soil Loss Equation is a good indicator of the potential of a storm to erode soil, as it quantifies the raindrop impact effect on the soil based on storm intensity. The R‐factor is defined as the average annual value of rainfall erosion index, EI, calculated by cumulating the EI values obtained for individual storms for at least 22 years. By definition, calculation of EI is based on rainfall measurements at short time intervals over which the intensity is essentially constant, i.e. using so‐called breakpoint data. Because of the scarcity of breakpoint rainfall data, many authors have used different time resolutions (Δt = 5, 10, 15, 30, and 60 min) to deduce EI in different areas of the world. This procedure affects the real value of EI because it is strongly dependent on Δt. In this contribution, after a general overview of similar studies carried out in different countries, the relationship between EI and Δt is explored in Calabria, southern Italy. The use of 17 139 storm events collected from 65 rainfall stations allowed the calculation of EI for different time intervals ranging from 5 to 60 min. The overall results confirm that calculation of EI is dependent on time resolution and a conversion factor able to provide its value for the required Δt is necessary. Based on these results, a parametric equation that gives EI as a function of Δt is proposed, and a regional map of the scale parameter a that represents the conversion factor for converting fixed‐interval values of (EI₃₀)_Δt to values of (EI₃₀)₁₅ is provided in order to calculate R anywhere in the region using rainfall data of 60 min. Copyright © 2015 John Wiley & Sons, Ltd.     

Number and volume raindrop size distributions in Taiwan

Knowledge of rainfall characteristics is very important for the accurate estimation of rainfall kinetic energy and prediction of soil erosion. In this study, a reliable and efficient data collection and analysis system was developed to analyse the natural raindrop data collected in subtropical Taiwan. Both raindrop size distributions by number and volume were carefully analysed. The seasonal variations of the rainfall erosivity factor R, which is an index of the erosive potential of rainfall and a function of rainfall kinetic energy, was also discussed. An isoerodent map of Taiwan was also developed based on the rainfall data recorded by 158 automated rainfall‐measuring stations within 26 years. Copyright © 2007 John Wiley & Sons, Ltd.     

Spatially varied soil erosion under different climates

Abstract

The spatial variability of the factors of the universal soil loss equation is examined on the mediterranean basin of Conca de Tremp covering 43.1 km² in Spain. The evaluation of the rainfall erosivity R and the soil erodibility K is relatively straightforward and spatially-averaged values of these parameters can be applied to the entire basin. Conversely, the spatial variability of annual soil erosion losses on large basins depends primarily on the factors L, S and C describing topographic, vegetation and land use parameters. A grid size analysis of soil erosion losses from the Conca de Tremp basin under mediterranean climatic conditions in Spain shows excellent agreement with the earlier results on the Chaudière basin in Canada. It is concluded for both basins that unbiassed estimates of soil erosion losses are obtained for grid sizes less than about 0.125 km². The analysis of the Conca de Tremp basin validates the use of the grid size factor proposed by Julien & Frenette (1987). It is also found that the grid size factor primarily depends on the average slope gradient which decreases with increasing grid size or drainage area. On the other hand, the grid size factor does not depend on the spatial variability of the factors R, K, L and C.     

Incorporating uncertainty analysis into a regional IDF formula

A rainfall intensity–duration–frequency (IDF) relationship was generated by pooling annual maximum rainfall series from 14 recording rain gauges in southern Taiwan. Dimensionless frequency curves, plotted by the growth curve method, can be well fitted by regression equations for a duration ranging from 10 mins to 24 hours. As the parameters in regression equations have a good statistical relationship with average annual rainfall, a generalized regional IDF formula was then formulated. The formula, based on average annual rainfall as an index, can be easily applied to non-recording rain gauges. This paper further applies the mean value first-order second moment (MFOSM) method to estimate the uncertainty of the proposed regional IDF formula. From a stochastic viewpoint, the generalized regional IDF formula can accurately simulate the IDF relationship developed using frequency analysis (EV1) at individual stations. The method can provide both rainfall intensity and variance isohyetal maps for various rainfall durations and return periods over the study area. © 1998 John Wiley & Sons, Ltd.     

Dynamic controls on wind erosion and dust generation on west‐central Free State agricultural land,South Africa

The west‐central part of South Africa's Free State Province falls within the transition zone between South Africa's sub‐humid, temperate grasslands to the east, and the semi‐arid Karoo and arid Kalahari to the south and west, respectively. The area is characterized by low rainfall (typically 500 mm or less) with high variability, but environmental conditions allow widespread dryland commercial agriculture (maize, sunflowers and stock farming). However, human activity promotes wind erosion and the area is susceptible to dust emissions. This study is the first to quantify the degree of wind erosion on the agricultural soils in the region under prevailing winter to spring climatic conditions and land management practices. Using arrays of cup anemometers, dust deposition traps and saltation impact sensors (Safires), measurements were made of the key erosivity and erodibility drivers that control the degree of wind erosion. Results demonstrate that significant quantities of dust are mobilized, particularly during the months of September and October. Thresholds of wind erosion are shown to respond particularly closely to changes in surface and aerodynamic roughness (z₀) with the amount of collected dust correlating well with measures of wind erosivity that weight the impact of higher wind speeds. Given the importance of surface roughness in controlling erosion thresholds, results show that the opportunity exists for well designed farming practices to control wind erosion. However, it is likely that climatically driven environmental change will impact on some of the identified controls on erosion (wind power, moisture availability) with the result that the wind erosion hazard is likely to increase within this marginal environment. Copyright © 2010 John Wiley & Sons, Ltd.     

The effect of wind averaging time on wind erosivity estimation

The Wind Erosion Prediction System (WEPS) and Revised Wind Erosion Equation (RWEQ) are widely used for estimating wind‐induced soil erosion at a field scale. Wind is the principal erosion driver in the two models. Wind erosivity, which describes the capacity of wind to cause soil erosion, is defined as erosive wind power density (WPD) in WEPS, and wind value (W) in RWEQ. In this study, the daily average WPD (AWPD) and the daily average W (Wf) were chosen to investigate the effect of averaging time on wind erosivity estimation based on observed wind data. We compare the daily AWPD and Wf calculated from 1, 5, 10, 15, 30, and 60 minute average wind speed data. The results of comparisons indicate that averaging wind speed can significantly influence estimates of wind erosivity. Compared with the daily AWPD and Wf calculated from one minute average wind speed data, all daily AWPD and Wf values calculated from 5, 10, 15, 30, and 60 minute averaged wind speeds tend to be significantly lower than values calculated from one minute values. In general, longer averaging times tend to produce smaller values of daily AWPD or Wf, which may lead to an under‐estimation of wind erosion. Further studies are needed to extend and apply the findings obtained in this study to actual wind erosion predictions. Copyright © 2012 John Wiley & Sons, Ltd.     

Orders of magnitude increase in soil erosion associated with land use change from native to cultivated vegetation in a Brazilian savannah environment

The Brazilian savanna (cerrado) is a large and important economic and environmental region that is experiencing significant loss of its natural landscapes due to pressures of food and energy production, which in turn has caused large increases in soil erosion. However the magnitude of the soil erosion increases in this region is not well understood, in part because scientific studies of surface runoff and soil erosion are scarce or nonexistent in the cerrado as well as in other savannahs of the world. To understand the effects of deforestation we assessed natural rainfall‐driven rates of runoff and soil erosion on an undisturbed tropical woodland classified as ‘cerrado sensu stricto denso’ and bare soil. Results were evaluated and quantified in the context of the cover and management factor (C‐factor) of the Universal Soil Loss Equation (USLE). Replicated data on precipitation, runoff, and soil loss on plots (5 × 20 m) under undisturbed cerrado and bare soil were collected for 77 erosive storms that occurred over 3 years (2012 through 2014). C‐factor was computed annually using values of rainfall erosivity and soil loss rate. We found an average runoff coefficient of ~20% for the plots under bare soil and less than 1% under undisturbed cerrado. The mean annual soil losses in the plots under bare soil and cerrado were 12.4 t ha^‐1 yr^‐1 and 0.1 t ha^‐1 yr^‐1, respectively. The erosivity‐weighted C‐factor for the undisturbed cerrado was 0.013. Surface runoff, soil loss and C‐factor were greatest in the summer and fall. Our results suggest that shifts in land use from the native to cultivated vegetation result in orders of magnitude increases in soil loss rates. These results provide benchmark values that will be useful to evaluate past and future land use changes using soil erosion models and have significance for undisturbed savanna regions worldwide. Copyright © 2015 John Wiley & Sons, Ltd.     

Introducing the recharge–discharge relationship to evaluate the recharge coefficient of karstic aquifers: thematic and indexing approaches

ABSTRACT

Karst aquifers and springs are important with respect to their potential for supplying drinking water in regions suffering from water scarcity in Iran. Accordingly, it is essential to determine the recharge potential of the catchment and the regions with higher obtainability potential. This study provides a road map for the Sheshpeer catchment in southern Iran. A recharge potential (RP) map was produced from which a recharge index (RI) was computed for several selected springs in the catchment. Furthermore, the unit discharge (q) – defined as the average annual discharge for a given catchment area and unit rainfall depth for each spring – was calculated. The plot of q versus RI for the springs showed a linear positive relationship between the two variables (R ² = 0.9). Applying the trend equation of this plot to the whole Sheshpeer karstic catchment reveals that its long-term recharge coefficient is 0.74.     

A new procedure to estimate the RUSLE EI30 index, based on monthly rainfall data and applied to the Algarve region, Portugal

One of the best indicators of the potential erosion risks is the rainfall–runoff erosivity factor (R) of the revised universal soil loss equation (RUSLE). Frequently, however, there is not enough data available to compute the R value, and other parameters, such as the modified Fournier index (F_mod), are used instead. But RUSLE is less effective if only the alternative procedures exist. One of the major discrepancies between R and the alternative parameters is time resolution: individual storms are used to calculate R while monthly averages over the year are used to calculate F_mod.

In this study, a multiple linear regression (r²=0.89) involving monthly EI₃₀, monthly rainfall for days with ≥10.0 mm and monthly number of days with rainfall ≥10.0 mm, for the Algarve region, is presented. Twenty-seven years of monthly rainfall erosivity values were computed for the 32 standard daily-read raingauge stations of the Algarve region.