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     

Simulation of the size distribution and erosivity of raindrops and throughfall drops

This paper presents a model that simulates the size distribution and erosivity of raindrops and throughfall drops. It utilizes existing models of rainfall drop size distribution and fall velocity and combines them with newly collated evidence of throughfall drop size distributions. A sensitivity analysis reveals that the model is sensitive to parameters that are easily measured or estimated: rainfall intensity, the mean volume drop diameter of the intercepted throughfall, canopy cover, and canopy height. The results of the model may be used at two levels. Firstly, to calculate specifically the size and fall velocity of individual drops, parameters that are needed in studies examining the response of soil surfaces to forces applied by rainfall. Secondly, to produce erosivity indices, based on rainfall intensity but which take account of the effects of a vegetation canopy. The paper shows that while the kinetic energy of rainfall (E(0), J mm^?1 m^?2) may be calculated from an equation of the familiar form: the kinetic energy of throughfall under any canopy may be calculated by combining this equation with another that relates the energy of drops under a 100 per cent canopy cover (E(100)) and the canopy height: .     

Interception,drop size distributions and rainfall kinetic energy in four colombian forest ecosystems

Interception losses, rain and throughfall drop size spectra and kinetic energy were studied in four relatively undisturbed tropical forest ecosystems along a transect across the Central Andean Cordillera of Colombia at altitudes between 3000 and 1000 m above sea level. Interception amounts ranged from 11 to more than 20 per cent of the total rainfall and fell within the normal range of interception figures observed in natural tropical forests. Drop size spectra were established using the filter paper method; the drop size distributions of the open field rainfall were unimodal while the throughfall had bimodal distributions, with a higher percentage of the volume of rain falling as large drops. Disturbance of the natural forests, for example by logging activities or cattle grazing, will further increase the throughfall kinetic energy and may lead to higher splash erosion rates inside the forests than in the open field. The kinetic energy of the throughfall was higher than that of the open field rainfall (20-70 per cent), even after correcting for interception losses (4-30 per cent). Splash-cup experiments, conducted both in the field and in the laboratory, indicated that the kinetic energy is a good index of rainfall erosivity. Inside the forests the amounts of sand splashed from the splash-cups was, after correction for interception losses, 2-16 per cent higher than outside the forests.     

EROSION RATES ON UPLAND AREAS IN THE MIDSOUTH USA

lINTR0DUCTIONSoilerosi0n0fslopingagriculturallandoccursduringandimmediatelyafterrainfalleventswherethesoilsurfaceisn0tadequatelyprotectedfromtheerosiveforces0frainfallandrunoffThebestaltemativefortheconservationplarinerist0minimizesoilerosionasmuchasp0ssible.Managementofsoilsurfaceconditionsinvolvesexperienceandjudgmeniinselectingmanagementpracticesthatutilizecr0presiduesandminimizedisturbancebytillage.Concentfatedrun0ffs0metimescausessevereerosi0n,particularlyonunprotectedsteepslopes.Ril…     

Comparison of three types of laser optical disdrometers under natural rainfall conditions

ABSTRACT

Optical disdrometers can be used to estimate rainfall erosivity; however, the relative accuracy of different disdrometers is unclear. This study compared three types of optical laser-based disdrometers to quantify differences in measured rainfall characteristics and to develop correction factors for kinetic energy (KE). Two identical PWS100 (Campbell Scientific), one Laser Precipitation Monitor (Thies Clima) and a first-generation Parsivel (OTT) were collocated with a weighing rain gauge (OTT Pluvio²) at a site in Austria. All disdrometers underestimated total rainfall compared to the rain gauge with relative biases from 2% to 29%. Differences in drop size distribution and velocity resulted in different KE estimates. By applying a linear regression to the KE–intensity relationship of each disdrometer, a correction factor for KE between the disdrometers was developed. This factor ranged from 1.15 to 1.36 and allowed comparison of KE between different disdrometer types despite differences in measured drop size and velocity.     

Modelling raindrop impact and splash erosion processes within a spatial cell: a stochastic approach

A new approach is proposed to simulate splash erosion on local soil surfaces. Without the effect of wind and other raindrops, the impact of free‐falling raindrops was considered as an independent event from the stochastic viewpoint. The erosivity of a single raindrop depending on its kinetic energy was computed by an empirical relationship in which the kinetic energy was expressed as a power function of the equivalent diameter of the raindrop. An empirical linear function combining the kinetic energy and soil shear strength was used to estimate the impacted amount of soil particles by a single raindrop. Considering an ideal local soil surface with size of 1 m × 1 m, the expected number of received free‐falling raindrops with different diameters per unit time was described by the combination of the raindrop size distribution function and the terminal velocity of raindrops. The total splash amount was seen as the sum of the impact amount by all raindrops in the rainfall event. The total splash amount per unit time was subdivided into three different components, including net splash amount, single impact amount and re‐detachment amount. The re‐detachment amount was obtained by a spatial geometric probability derived using the Poisson function in which overlapped impacted areas were considered. The net splash amount was defined as the mass of soil particles collected outside the splash dish. It was estimated by another spatial geometric probability in which the average splashed distance related to the median grain size of soil and effects of other impacted soil particles and other free‐falling raindrops were considered. Splash experiments in artificial rainfall were carried out to validate the availability and accuracy of the model. Our simulated results suggested that the net splash amount and re‐detachment amount were small parts of the total splash amount. Their proportions were 0·15% and 2·6%, respectively. The comparison of simulated data with measured data showed that this model could be applied to simulate the soil‐splash process successfully and needed information of the rainfall intensity and original soil properties including initial bulk intensity, water content, median grain size and some empirical constants related to the soil surface shear strength, the raindrop size distribution function and the average splashed distance. Copyright © 2007 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.     

Impacts of annual precipitation extremes on soil and nutrient losses in vineyards of NE Spain

The objective of this research was to characterise annual precipitation extremes in a Mediterranean vineyard region. The number of exceptional events (P > 95th percentile) and annual extreme events (P > 99th percentile), as well as their strength, erosive character and return period were analysed for 2000–2004. The erosive character was evaluated according to the R‐factor (kinetic energy × maximum intensity in 30‐min periods). Soil and nutrient losses caused by these events were evaluated by combining field sampling and a hydrological model to estimate total runoff in a vineyard plot. The results show a clear increase in the number of very wet days and extreme events (P > 95th percentile), which represented up to 88% of annual rainfall. The severity of the extreme events (TS = precipitation event P > 99th percentile) reached values higher than 50 mm almost every year. These values were far exceeded in 2000, when one extraordinary event recorded 50% of the annual rainfall, with TS of 189 mm, about 80% of total rainfall being lost as runoff. Annual erosivity was driven not only by extreme events, but also by short events of less depth but high intensity. During some of the years analysed, rainfall erosivity was two or three times the average in the area. Most soil and nutrient losses occurred in a small number of events: one or two events every year were responsible for more than 75% of the annual soil and nutrient losses on average. Antecedent soil moisture conditions, runoff rates, and events with a return period higher than two years were responsible for the higher erosion rates. Apart from an exceptional event recorded in 2000, which produced more than 200 Mg ha^?1 soil losses, annual soil losses up to 25 Mg ha^?1 were recorded, which are much higher than the soil loss tolerance. Copyright © 2008 John Wiley & Sons, Ltd.     

Assessment of design rainfall uncertainty for hydrologic engineering applications in Hong Kong

Hydrologic engineering designs and analyses often require the specification of design storm which involves rainfall amount, duration and hyetograph. In practice, the determination of design rainfall in hydrologic engineering applications involves the frequency analysis of extreme rainfalls of different durations and the establishment of rainfall hyetograph for the design event under consideration. Sampling errors exist in the estimation of rainfall depth (or intensity) quantiles from frequency analysis, which will be transmitted in the process of determining the design rainfall hyetograph. This paper presents a practical methodological framework based on the bootstrap resampling scheme to assess the uncertainty features associated with the magnitude of estimated rainfall depth/intensity quantiles and the corresponding design hyetographs. The procedure is implemented to quantify uncertainty of design rainfall hyetograph following the Stormwater Drainage Manual of Hong Kong involving the use of rainfall intensity–duration–frequency (IDF) model. Of particular interesting is that the bootstrap resampling scheme implemented herein is modified to handle unequal record period of annual maximum rainfall data series of different durations and to account for their intrinsic correlations. According to the adopted rainfall IDF model, the design rainfall hyetograph is a function of the IDF model coefficients. Due to the correlation among rainfall quantiles of different durations, the IDF coefficients are found to be strongly related in a nonlinear fashion which should not be ignored in the establishment of the design hyetographs.     

Measuring drop size distribution and kinetic energy of rainfall using a force transducer

The relatively high cost of commercially available raindrop spectrometers and disdrometers has inhibited detailed and intensive research on drop size distribution, kinetic energy and momentum of rainfall which are important for understanding and modelling soil erosion caused by raindrop detachment. In this study, an approach to find the drop size distribution, momentum and kinetic energy of rainfall using a relatively inexpensive device that uses a piezoelectric force transducer for sensing raindrop impact response is introduced. The instrument continuously and automatically records, on a time‐scale, the amplitude of electrical pulses produced by the impact of raindrops on the surface of the transducer. The size distribution of the raindrops and their respective kinetic energy are calculated by analysing the number and amplitude of pulses recorded, and from the measured volume of total rainfall using a calibration curve. Simultaneous measurements of the instrument, a rain gauge and a dye‐stain method were used to assess the performance of the instrument. Test results from natural and simulated rainfalls are presented. Copyright © 2000 John Wiley & Sons, Ltd.     

The effect of a new zealand beech forest canopy on the kinetic energy of water drops and on surface erosion

Rain and throughfall drops were sampled during rain events in a New Zealand beech forest and the frequency distributions of drop mass and kinetic energy calculated. The kinetic energy of throughfall under the canopy was always greater than that of rainfall in the open, notwithstanding interception losses. During a typical rain event in which 51 mm fell in 36 h, the total kinetic energy of throughfail was 1.5 times greater than that of rainfall, and the mean amount of sand splashed from sample cups was 3.1 times greater under the canopy than in the open. It appears that where mineral soil is exposed at the surface, by animal trampling or burrowing for example, rates of soil detachment by splash under a forest canopy will probably exceed those in the open.     

Stochastic generation of hourly rainstorm events

Occurrence of rainstorm events can be characterized by the number of events, storm duration, rainfall depth, inter-event time and temporal variation of rainfall within a rainstorm event. This paper presents a Monte-Carlo based stochastic hourly rainfall generation model considering correlated non-normal random rainstorm characteristics, as well as dependence of various rainstorm patterns on rainfall depth, duration, and season. The proposed model was verified by comparing the derived rainfall depth–duration–frequency relations from the simulated rainfall sequences with those from observed annual maximum rainfalls based on the hourly rainfall data at the Hong Kong Observatory over the period of 1884–1990. Through numerical experiments, the proposed model was found to be capable of capturing the essential statistical features of rainstorm characteristics and those of annual extreme rainstorm events according to the available data.     

Application of piezoelectric transducers in simulated rainfall erosivity assessment

Abstract

Rainfall simulators have often been used to mimic natural rainfall for studies of various land-surface and water interaction processes. The characteristics of the simulated rainfall are the main indicators used to judge the performance of the rainfall simulators. The aim of this study is to investigate the potential of piezoelectric transducers for measuring and evaluating a dripper-type simulated rainfall drop-size distribution (DSD) and kinetic energy (KE). The directly measured KE was significantly correlated with the estimated KE using the drop-size distribution and empirical rain drop fall velocity relationships. This result emphasizes the potential use of the piezoelectric sensor to directly measure and evaluate rainfall kinetic energy. Also, the relationship between rainfall intensity and KE showed good patterns of agreement between simulated rainfall and natural rainfall.

Citation Abd Elbasit, M. A. M., Yasuda, H. & Salmi, A. (2011) Application of piezoelectric transducers in simulated rainfall erosivity assessment. Hydrol. Sci. J. 56(1), 187–194.     

The relationship between rainfall inputs and flood generation in south–east Spain

Rainfall and flood data are relatively sparse in semi‐arid areas; hence there have been relatively few investigations into the relationships between rainfall inputs and flood generation in these environments. Previous work has shown that flood properties are influenced by a combination of precipitation characteristics including amount, intensity, duration and spatial distribution. Therefore floods may be produced by high intensity, short duration storms, or longer duration, low intensity rainfall. Most of this research has been undertaken in small catchments in either hyper‐arid or relatively high rainfall Mediterranean climates. This paper presents results from a 6 year data record in south‐east Spain from research conducted in two basins, the Rambla Nogalte (171 km²) and the Rambla de Torrealvilla (200 km²). Data cover an area of approximately 500 km² and an annual average rainfall of 300 mm. At coarse temporal resolutions gauges spread over large areas record similar patterns of rainfall, although spells of rain show much more complexity; pulses of rain within storms can vary considerably in total rainfall, intensity and duration over the same area. The analysis for south‐east Spain shows that most storms occur over a period of less than 24 h, but that the number of rainfall events declines as the duration exceeds 8 h. This is at odds with data on floods for the study area suggesting that they are produced by storms lasting longer than 18 h. However, one flood event was produced by a very short (15 min) storm with high intensity rainfall. Most floods tended to occur in May/June or September, which coincides with wetter months of the year (September, October, December and May). Floods are also more highly related to the total rainfall occurring in a spell of rain, than to intensity. The complexity of storm rainfall increases with the storm total, which makes it difficult to generalize on the importance of rainfall intensity for flood generation. Copyright © 2007 John Wiley & Sons, Ltd.     

Characterization of rainfall generated by dripper‐type rainfall simulator using piezoelectric transducers and its impact on splash soil erosion

Rainfall erosivity represents the primary driver for particle detachment in splash soil erosion. Several raindrop erosivity indices have been developed in order to quantify the potential of rainfall to cause soil erosion. Different types of rainfall simulators have been used to relate rainfall characteristics to soil detachment. However, rainfall produced by different rainfall simulators has different characteristics, specifically different relationships between rainfall intensity and rainfall erosivity. For this reason, the effect of rainfall characteristics produced by a dripper‐type rainfall simulator on splash soil erosion (D_s) has been investigated. The simulated rainfall kinetic energy (KE) and drop size distribution (DSD) were measured using piezoelectric transducers, modified from the Vaisala RAINCAP^® rain sensor. The soil splash was evaluated under various simulated rainfall intensities ranging from 10 to 100 mm h^?1 using the splash‐cup method. The simulated rainfall intensity (I) and kinetic energy relationship (I–KE) was found to be different from natural rainfall. The simulated rainfall intensity and splash soil erosion relationship (I–D_s) also followed this same trend. The I–KE relationship was found to follow the natural rainfall trend until the rainfall intensity reached 30 mm h^?1 and above this limit the KE started to decrease. This emphasizes the importance of the I–KE relationship in determining the I‐D_s relationship, which can differ from one rainfall simulator to another. D_s was found to be highly correlated with KE (r = 0·85, P < 0·001), when data produced by the rainfall intensity ranged from 10 to 100 mm h^?1. However, when the threshold rainfall intensity (30 mm h^?1) was considered, the correlation coefficient further improved (r = 0·89, P = 0·001). Accordingly, to improve the soil splash estimation of simulated rainfall under various rainfall intensities the I–KE characterization relationship for rainfall simulators has to be taken into account. Copyright © 2010 John Wiley & Sons, Ltd.     

Interrill erosion on cultivated Greek soils: modelling sediment delivery

For interrill erosion, raindrop‐induced detachment and transport of sediment by rainfall‐disturbed sheet flow are the predominant processes, while detachment by sheet flow and transport by raindrop impact are negligible. In general, interrill subprocesses are inter‐actively affected by rainfall, soil and surface properties. The objective of this work was to study the relationships among interrill runoff and sediment loss and some selected para‐meters, for cultivated soils in central Greece, and also the development of a formula for predicting single storm sediment delivery. Runoff and soil loss measurement field experiments have been conducted for a 3·5‐year period, under natural storms. The soils studied were developed on Tertiary calcareous materials and Quaternary alluvial deposits and were textured from sandy loam to clay. The second group of soils showed greater susceptibility to sealing and erosion than the first group. Single storm sediment loss was mainly affected by rain and runoff erosivity, being significantly correlated with rain kinetic energy (r = 0·64***), its maximum 30‐minute intensity (r = 0·64***) and runoff amount (r = 0·56***). Runoff had the greatest correlation with rain kinetic energy (r = 0·64***). A complementary effect on soil loss was detected between rain kinetic energy and its maximum 30‐minute intensity. The same was true for rain kinetic energy and topsoil aggregate instability, on surface seal formation and thus on infiltration characteristics and overland flow rate. Empirical analysis showed that the following formula can be used for the successful prediction of sediment delivery (D_i): D_i = 0·638βEI₃₀tan(θ) (R² = 0·893***), where β is a topsoil aggregate instability index, E the rain kinetic energy, I₃₀ the maximum 30‐minute rain intensity and θ the slope angle. It describes soil erodibility using a topsoil aggregate instability index, which can be determined easily by a simple laboratory technique, and runoff through the product of this index and rain kinetic energy. Copyright © 2006 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.     

Use of the Connectivity of Runoff Model (CRUM) to investigate the influence of storm characteristics on runoff generation and connectivity in semi‐arid areas

Much attention has been given to the surface controls on the generation and transmission of runoff in semi‐arid areas. However, the surface controls form only one part of the system; hence, it is important to consider the effect that the characteristics of the storm event have on the generation of runoff and the transmission of flow across the slope. The impact of storm characteristics has been investigated using the Connectivity of Runoff Model (CRUM). This is a distributed, dynamic hydrology model that considers the hydrological processes relevant to semi‐arid environments at the temporal scale of a single storm event. The key storm characteristics that have been investigated are the storm duration, rainfall intensity, rainfall variability and temporal structure. This has been achieved through the use of a series of defined storm hydrographs and stochastic rainfall. Results show that the temporal fragmentation of high‐intensity rainfall is important for determining the travel distances of overland flow and, hence, the amount of runoff that leaves the slope as discharge. If the high‐intensity rainfall is fragmented, then the runoff infiltrates a short distance downslope. Longer periods of high‐intensity rainfall allow the runoff to travel further and, hence, become discharge. Therefore, storms with similar amounts of high‐intensity rainfall can produce very different amounts of discharge depending on the storm characteristics. The response of the hydrological system to changes in the rainfall characteristics can be explained using a four‐stage model of the runoff generation process. These stages are: (1) all water infiltrating, (2) the surface depression store filling or emptying without runoff occurring, (3) the generation and transmission of runoff and (4) the transmission of runoff without new runoff being generated. The storm event will move the system between the four stages and the nature of the rainfall required to move between the stages is determined by the surface characteristics. This research shows the importance of the variable‐intensity rainfall when modelling semi‐arid runoff generation. The amount of discharge may be greater or less than the amount that would have been produced if constant rainfall intensity is used in the model. Copyright © 2006 John Wiley & Sons, Ltd.