Field and computer simulation experiments on the formation of desert pavement

A series of rainfall simulation experiments was carried out at the Walnut Gulch Experimental Watershed, Tombstone, Arizona (31° 43′N, 110° 41′W), to observe the speed at which desert pavement surfaces could be re-established following disturbance. The results of these experiments, which consisted of repeated, 5 min rainfall events, demonstrate that pavements can reform within 10 events, which is compatible with observations of the recovery of surfaces under natural rainfall on an annual cycle. A model for the development of pavements by raindrop erosion processes had previously shown the importance of these processes. The rainfall simulation experiments were used to test the general applicability of this model. The model was able to reproduce the general characteristics of the regenerated surfaces and the timing of their development. However, details of the particle size fractions produced were less well simulated by the model. Testing of the sensitivity of the model to the sediment transport parameters suggests that this problem is not related to the soil characteristics, but is more likely to be an indication of a poor understanding of all the feedbacks operating in the raindrop erosion processes. Copyright © 1999 John Wiley & Sons, Ltd.     

Unsteady soil erosion due to rainfall impact: a model of sediment sorting on the hillslope

A new method is presented for predicting sediment sorting associated with soil erosion by raindrop impact for non-equilibrium conditions. The form of soil erosion considered is that which results from raindrop impact in the presence of shallow overland flow itself where the flow is not capable of eroding sediment. The method specifically considers early time runoff and erosion when sediment leaving an eroding area is generally finer and thus may have a higher potential for transport of sorbed pollutants. The new mechanism described is the formation of a deposited layer on the soil surface, which is shown to lead to sediment sorting during an erosion event. The deposited layer is taken to have two roles in this process: to temporarily store sediment on the surface between successive trajectories, and to shield the underlying soil from erosive stresses. Equations describing the dynamics of the suspended sediment mixture and the deposited layer are developed. By integrating these equations over the length of eroding land element and over the duration of the erosion event, an event-based solution is proposed which predicts total sediment sorting over the event. This solution is shown to be consistent with experimentally observed trends in enrichment of fine sediment. Predictions using this approach are found to only partly explain measured enrichment for sets of experimental data for two quite different soils, but to be in poor agreement for an aridsol of dispersive character. It is concluded that the formation of the deposited layer is a significant mechanism in the enrichment of fine sediment and associated sorbed pollutants, but that processes in the dispersive soil are not as well described by the theory presented.     

Desert pavement dynamics: numerical modeling and field‐based calibration

Desert pavements are widely used as a relative surface‐dating tool because they are progressively better developed on surfaces ranging from thousands to hundreds of thousands of years in age. Recent work, however, has highlighted the dynamic nature of pavements and undermined their use as surface‐age indicators. Quade (2001) proposed that latest Pleistocene vegetation advances destroyed all Mojave Desert pavements above 400 m elevation, making all such pavements Holocene in age. In an effort to reconcile young‐pavement evidence with their widespread use as Pleistocene surface‐age indicators, we developed a numerical model based on the classic conceptual model in which pavements co‐evolve with their underlying eolian epipedons over millennial timescales. In this co‐evolutionary process, fine‐grained eolian deposition and A_v‐horizon development within the eolian epipedon promotes surface clast motion and pavement development, enhancing the eolian‐sediment‐trapping ability of the pavement in a positive feedback. Model results illustrate the multi‐scale nature of pavement dynamics: pavements may require tens of thousands of years to fully develop from a newly abandoned alluvial surface, but may heal over timescales of decades to centuries if a mature eolian epipedon is present. As such, there is no inconsistency between rapid pavement healing and a Pleistocene age for the underlying alluvial surface. To calibrate the model, we conducted surficial geologic mapping and pavement‐sedimentological analysis on two desert piedmonts. Our study areas include both proximal and distal fan environments, illustrating the role of parent‐material texture in controlling the mode of pavement formation. Using available geochronology, our work provides a rigorous calibration of pavement formation rates in our study areas and provides evidence supporting the use of pavements as local relative surface‐age indicators over Holocene to late Pleistocene timescales. Copyright © 2007 John Wiley & Sons, Ltd.     

Dating desert pavements – First results from a challenging environmental archive

Desert pavements are widespread landforms of arid environments. They consist of a monolayer of clasts at the surface, associated with an underlying unit of eolian fines. In some situations, buried desert pavements can be observed, which is interpreted as a change in the environmental conditions. Therefore, it is believed that desert pavements represent important paleoenvironmental sediment archives, especially for arid environments, where natural archives of past environments are rare. To better understand the formation process of desert pavements and to enable the paleoenvironmental interpretation of these valuable sediment archives, reliable chronologies are of crucial importance. Thus, OSL dating was applied to samples from well-developed desert pavements in two different study areas, the Cima Volcanic Field, eastern Mojave Desert, USA, and the desert of northeastern Badia, Jordan. To test the suitability of the sediments for OSL dating, the luminescence characteristics of the fine- and coarse-grain quartz fraction are described and compared. Finally, first OSL ages are presented.     

Raindrop‐induced saltation and the enrichment of sediment discharged from sheet and interrill erosion areas

Sheet and interrill erosion areas are sources of soil material rich in nutrients and pollutants. The loss of soil, nutrients and other chemicals from these areas is a matter of concern both in terms of maintaining soil productivity and the health of offsite environments. Many experiments on rainfall erosion have shown enrichment of fine material, nutrients and other chemicals in the sediment discharged for sheet and interrill erosion areas, but often these results were obtained over short periods of time. A qualitative mechanistic model of raindrop‐induced saltation is used to illustrate how this transport mechanism influences the composition of sediment discharged by rain‐impacted flow. Initially, fast moving particles are enriched in the sediment discharge but, over time, during a rainfall event, slower moving particles become more represented. Raindrop‐induced saltation promotes the storage of material on the soil surface with a coarser composition than the original soil. Winnowing of material from this storage by the development of flow‐driven saltation during high‐intensity events can modify the composition of the sediment discharged later by raindrop‐induced saltation. Given stable soil particles, the composition of the sediment discharged at the steady state is the same as the original soil. Enrichment is a non‐steady‐state phenomenon and failure to recognize the transient nature of enrichment may lead to inappropriate interpretation of the implications of the results from short‐term experiments. Copyright © 2011 John Wiley & Sons, Ltd.     

Evaluation of an interrill soil erosion model using laboratory catchment data

Physically based soil erosion simulation models require input parameters of soil detachment and sediment transport owing to the action and interactions of both raindrops and overland flow. A simple interrill soil water transport model is applied to a laboratory catchment to investigate the application of raindrop detachment and transport in interrill areas explicitly. A controlled laboratory rainfall simulation study with slope length simulation by flow addition was used to assess the raindrop detachment and transport of detached soil by overland flow in interrill areas. Artificial rainfall of moderate to high intensity was used to simulate intense rain storms. However, experiments were restricted to conditions where rilling and channelling did not occur and where overland flow covered most of the surface. A simple equation with a rainfall intensity term for raindrop detachment, and a simple sediment transport equation with unit discharge and a slope term were found to be applicable to the situation where clear water is added at the upper end of a small plot to simulate increased slope length. The proposed generic relationships can be used to predict raindrop detachment and the sediment transport capacity of interrill flow and can therefore contribute to the development of physically‐based erosion models. Copyright © 1999 John Wiley & Sons, Ltd.     

A physically-based model of the dispersion of splash droplets ejected from a water drop impact

A model has been developed which predicts the dispersion of splash droplets produced by the impact of a water drop on a sloping soil surface. Experimental results of the ejection velocities and ejection angles of the splash droplets are generalized to a planar slope and the resulting splash distances are calculated taking into account the effects of air resistance. The predictions are presented in terms of the numbers of splash droplets from the impact point to surrounding squares arranged in a grid on the slope. The model explains many experimentally observed features of raindrop splash in terms of the mechanics of the processes involved and can make predictions of the effects of slope, wind, raindrop size, and soil properties on droplet dispersion. The component of the raindrop velocity parallel to the surface of the slope is identified as the main factor determining the degree and the direction of the asymmetry in the splash droplet dispersion. By combining the model with a theory of the entrapment of soil in the splash droplets it is possible to extend it to predict the dispersion of soil particles by raindrop impact, which is the basis of a model of soil erosion by rainsplash.     

Physically based modelling of sediment generation and transport under a large rainfall simulator

A series of large rainfall simulator experiments was conducted in 2002 and 2003 on a small plot located in an experimental catchment in the North Island of New Zealand. These experiments measured both runoff and sediment transport under carefully controlled conditions. A physically based hydrological modelling system (SHETRAN) was then applied to reproduce the observed hydrographs and sedigraphs. SHETRAN uses physically based equations to represent flow and sediment transport, and two erodibility coefficients to model detachment of soil particles by raindrop erosion and overland flow erosion. The rate of raindrop erosion also depended on the amount of bare ground under the simulator; this was estimated before each experiment. These erodibility coefficients were calibrated systematically for summer and winter experiments separately, and lower values were obtained for the summer experiments. Earlier studies using small rainfall simulators in the vicinity of the plot also found the soil to be less erodible in summer and autumn. Limited validation of model parameters was carried out using results from a series of autumn experiments. The modelled suspended sediment load was also sensitive to parameters controlling the generation of runoff from the rainfall simulator plot; therefore, we found that accurate runoff predictions were important for the sediment predictions, especially from the experiments where the pasture cover was good and overland flow erosion was the dominant mechanism. The rainfall simulator experiments showed that the mass of suspended sediment increased post‐grazing, and according to the model this was due to raindrop detachment. The results indicated that grazing cattle or sheep on steeply sloping hill‐country paddocks should be carefully managed, especially in winter, to limit the transport of suspended sediment into watercourses. Copyright © 2006 John Wiley & Sons, Ltd.     

The impact of slope length on the discharge of sediment by rain impact induced saltation and suspension

Simulations using a mechanistic model of raindrop driven erosion in rain‐impacted flow were performed with particles travelling by suspension, raindrop induced saltation and flow driven saltation. Results generated by both a high intensity storm, and a less intense one, indicate that, because of the effect of flow depth on the delivery of raindrop energy to the bed, there is a decline in sediment concentration, and hence soil loss per unit area, with slope length when particles are transported by raindrop induced saltation. However, that decline is reversed when the critical velocities that lead to flow driven saltation are episodically exceeded during an event. The simulations were performed on smooth surfaces and a single drop size but the general relationships are likely to apply for rain made up of a wide range of drop size. Although runoff is not always produced uniformly, as a general rule, flow velocities increase with slope length so that, typically, the distance particles travel before being discharged during an event increase with slope length. The effect of slope length on soil loss per unit area is often considered to vary with slope length to a power greater than zero and less that 1·0. The simulations show that effect of slope length on sediment discharge is highly dependent on the variations in runoff response resulting from variations in rainfall duration‐intensity‐infiltration conditions rather than plot length per se. Consequently, predicting soil loss per unit area using slope length with positive powers close to zero when sheet erosion occurs may not be as effective as commonly expected. Erosion by rain‐impacted flow is a complex process and that complexity needs to be considered when analysing the results of experiments associated with rain‐impacted flow under both natural and artificial conditions. Copyright © 2009 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.     

Numerical modeling of gravitational erosion in rill systems

A self-organizing model was developed for simulating rill erosion process on slopes with particular attention to the role of gravitational erosion.For a complete simulation circle,processes such as precipitation,infiltration,runoff,scouring,gravitational erosion and elevation variation were fully considered.Precipitation time(or runoff time) was regarded as iteration benchmark in the model.To specify the contribution of gravitational erosion to the process of rill formation and development,a gravitational erosion module was inserted into the model.Gravitational erosion in rill development was regarded as a Gaussian random process.A model was calibrated by our experimental data,and further validated satisfactorily with 22 runs of experimental results from different investigators. Systematic comparison was made between sediment yields with and without consideration of gravitational erosion module.It was demonstrated that the model could reasonably simulate the rill erosion process under a variety of slope gradients,rainfall intensities and soil conditions upon the gravitational erosion being considered.However,the role of gravitational erosion on sediment yields in rill systems varies significantly under different conditions,although it is of the utmost importance in steeper slopes.The process of gravitational erosion in rill development was studied by a newly-defined parameter a>,which is defined as the volume ratio of gravitational erosion over hydraulic-related erosion.The gravitational contribution to the total erosion could be over 50%for the rill systems with higher rainfall intensity and steeper slopes.     

A model of the redistribution of disaggregated soil particles by rainsplash

A model of the dispersion of splash droplets from a single raindrop impact on a sloping soil surface is combined with a theory of the entrainment of mineral particles from a disaggregated mixture in splash droplets to obtain a model of the dispersion of such particles by a raindrop impact. Stochastic modelling techniques extend this further to a model of the spatial redistribution of soil on a plot after a period of rainfall. Since the model is probabilistic and physically based it enables the incorporation of further advances in the understanding of splash erosion at all stages and can simulate the effect of the stochastic nature of rainfall and soil properties on the process. Several different situations are simulated. These include the movement of marked soil particles from point sources and the spatial patterns of erosion on a sloping plot. The model can also simulate the differential erosion of different soil particle size fractions.     

Modelling water erosion in the Sahel: application of a physically based soil erosion model in a gentle sloping environment

Water is a major limiting factor in arid and semi‐arid agriculture. In the Sahelian zone of Africa, it is not always the limited amount of annual rainfall that constrains crop production, but rather the proportion of rainfall that enters the root zone and becomes plant‐available soil moisture. Maximizing the rain‐use efficiency and therefore limiting overland flow is an important issue for farmers. The objectives of this research were to model the processes of infiltration, runoff and subsequent erosion in a Sahelian environment and to study the spatial distribution of overland flow and soil erosion. The wide variety of existing water erosion models are not developed for the Sahel and so do not include the unique Sahelian processes. The topography of the Sahelian agricultural lands in northern Burkina Faso is such that field slopes are generally low (0–5°) and overland flow mostly occurs in the form of sheet flow, which may transport large amounts of fine, nutrient‐rich particles despite its low sediment transport capacity. Furthermore, pool formation in a field limits overland flow and causes resettlement of sediment resulting in the development of a surface crust. The EUROSEM model was rewritten in the dynamic modelling code of PCRaster and extended to account for the pool formation and crust development. The modelling results were calibrated with field data from the 2001 rainy season in the Katacheri catchment in northern Burkina Faso. It is concluded that the modified version of EUROSEM for the Sahel is a fully dynamic erosion model, able to simulate infiltration, runoff routing, pool formation, sediment transport, and erosion and deposition by inter‐rill processes over the land surface in individual storms at the scale of both runoff plots and fields. A good agreement is obtained between simulated and measured amounts of runoff and sediment discharge. Incorporating crust development during the event may enhance model performance, since the process has a large influence on infiltration capacity and sediment detachment in the Sahel. Copyright © 2005 John Wiley & Sons, Ltd.     

MODELING OF INTERRILL EROSION IN THE LOESS PLATEAU OF CHINA

IINTRODUCTIONTheinterrillerosiononafieldplotisaffectedbythekineticenergyoftherainfall,wind,topographyfactors,propertiesofsoilandthecanopy.Theinterrillerosionoccursasthefirstdropimpactsthehillslopes.Theinterrillerosionoccursinallkindsofrainfallandtheamountofthesplasherosion,whichisthemainpartofinterrillerosion,canaccountforagreatpanofthetotalerosionamountinaheavystorm(Baner1990,Glymph1957,QianandWan1986,Zhou1981).Therefore,itisveryimportanttorevealthemechanismtoestimatetheamountofinterri…     

Vegetation controls on small‐scale runoff and erosion dynamics in a degrading dryland environment

This paper investigates the controls of vegetation on runoff and erosion dynamics in the dryland environment of Jornada, New Mexico, USA. As the American southwest has seen significant shifts in the dominant vegetation species in the past 150 years, an understanding of the vegetation effects on hydrological and erosional processes is vital for understanding and managing environmental change. Small‐scale rainfall simulations were carried out to identify the hydrological and erosional processes resulting from the grassland and shrubland vegetation species. Results obtained using tree‐regression analysis suggested that the primary vegetation control on runoff and erosion is the shrub type and canopy density, which directly affects the local microtopographic gradient of mounds beneath the shrubs. Significant interactions and feedbacks were found to occur among the local mound gradient, crust cover, soil aggregate stability and antecedent soil moisture between the different vegetation species for both the runoff and erosion responses. Although some of the shrub species were found to produce higher sediment yields than the grass species, the distinguishing feature of the grassland was the significantly higher enrichment in the fine sediment fraction compared to all other surface cover types. This enrichment in fines has important implications for nutrient movement in such environments. Copyright © 2009 John Wiley & Sons, Ltd.     

Modelling wind‐blown sediment transport around single vegetation elements

Wind erosion is an important soil erosion and hence a soil degradation problem in the Sahelian zone of West Africa. Potentially, the characteristic dryland vegetation with scattered trees and shrubs can provide for soil erosion protection from wind erosion, but so far adequate quantification of vegetation impacts is lacking. The aim of this study was to develop a model of wind‐blown soil erosion and sediment transport around a single shrub‐type vegetation element. Starting with the selection of a suitable transport equation from four possible sediment transport equations, the effects of a single vegetation element on wind speed were parameterized. The modified wind speed was then applied to a sediment transport equation to model the change in sediment mass flux around a shrub. The model was tested with field data on wind speed and sediment transport measured around isolated shrubs in a farmer's field in the north of Burkina Faso. The simple empirical equation of Radok (Journal of Glaciology 19 : 123–129, 1977) performed best in modelling soil erosion and sediment transport, both for the entire event duration and for each minute within an event. Universal values for the empirical constants in the sediment transport equation could not be obtained because of the large variability in soil and roughness characteristics. The pattern of wind speed, soil erosion and sediment transport behind a shrub and on either side of it was modelled. The wind speed changed in the lee of the vegetation element depending on its porosity, height and downwind position. Wind speed was recovered to the upstream speed at a downwind distance of 7·5 times the height of the shrub. The variability in wind direction created a ‘rotating’ area of influence around the shrub. Compared to field measurements the model predicted an 8% larger reduction in sediment transport in the lee of the vegetation element, and a 22% larger increase beside the vegetation element. Copyright © 2011 John Wiley & Sons, Ltd.     

Soil detachment and transport on field‐ and laboratory‐scale interrill areas: erosion processes and the size‐selectivity of eroded sediment

Field‐ and laboratory‐scale rainfall simulation experiments were carried out in an investigation of the temporal variability of erosion processes on interrill areas, and the effects of such variation upon sediment size characteristics. Poorly aggregated sandy soils from the semi‐arid environment of Senegal, West Africa, were used on both a 40 m² field plot and a 0·25 m² laboratory plot; rainfall intensity for all experiments was 70 mm h^?1 with a duration of 1 to 2 hours. Time‐series measurements were made of the quantity and the size distribution of eroded material: these permitted an estimate of the changing temporal balance between the main erosion processes (splash and wash). Results from both spatial scales showed a similar temporal pattern of runoff generation and sediment concentration. For both spatial scales, the dominant erosional process was detachment by raindrops; this resulted in a dynamic evolution of the soil surface under raindrop impact, with the rapid formation of a sieving crust followed by an erosion crust. However, a clear difference was observed between the two scales regarding the size of particles detached by both splash and wash. While all measured values were lower than the mean weight diameter (MWD) value of the original soil (mean 0·32 mm), demonstrating the size‐selective nature of wash and splash processes, the MWD values of washed and splashed particles at the field scale ranged from 0·08 to 0·16 mm and from 0·12 to 0·30 mm respectively, whereas the MWD values of washed and splashed particles at the laboratory scale ranged from 0·13 to 0·29 mm and from 0·21 to 0·32 mm respectively. Thus only at the field scale were the soil particles detached by splash notably coarser than those transported by wash. This suggests a transport‐limited erosion process at the field scale. Differences were also observed between the dynamics of the soil loss by wash at the two scales, since results showed wider scatter in the field compared to the laboratory experiments. This scatter is probably related to the change in soil surface characteristics due to the size‐selectivity of the erosion processes at this spatial scale. Copyright © 2006 John Wiley & Sons, Ltd.     

Effects of rainfall patterns on runoff and rainfall-induced erosion

Rainfall-induced erosion involves the detachment of soil particles by raindrop impact and their transport by the combined action of the shallow surface runoff and raindrop impact.Although temporal variation in rainfall intensity(pattern)during natural rainstorms is a common phenomenon,the available information is inadequate to understand its effects on runoff and rainfall-induced erosion processes.To address this issue,four simulated rainfall patterns(constant,increasing,decreasing,and increasing-decreasing)with the same total kinetic energy were designed.Two soil types(sandy and sandy loam)were subjected to simulated rainfall using 15 cm×30 cm long detachment trays under infiltration conditions.For each simulation,runoff and sediment concentration were sampled at regular intervals.No obvious difference was observed in runoff across the two soil types,but there were significant differences in soil losses among the different rainfall patterns and stages.For varying-intensity rainfall patterns,the dominant sediment transport mechanism was not only influenced by raindrop detachment but also was affected by raindrop-induced shallow flow transport.Moreover,the efficiency of equations that predict the interrill erosion rate increased when the integrated raindrop impact and surface runoff rate were applied.Although the processes of interrill erosion are complex,the findings in this study may provide useful insight for developing models that predict the effects of rainfall pattern on runoff and erosion.     

Modeling of the effect of flow depth on sediment discharged by rain‐impacted flows from sheet and interrill erosion areas: a review

Sediment, nutrients and pollutants discharged from sheet and interrill erosion areas by rain‐impacted flows may influence water quality in streams and rivers. The depth of water on the soil surface influences the capacity of raindrop impacts to detach soil material underlying rain‐impacted flows, and a number of so‐called process‐based and mechanistic models erroneously use equations on the basis of the effect of water depth on splash erosion to account for this effect. Also, a number of these models require complex mathematical solutions to make them operate and can only predict sediment composition and discharges well if many of their parameters are calibrated specifically to the situations where they are being applied. Experiments with rain‐impacted flows, where flow depth and velocity over eroding surfaces have been controlled, have been reported in the literature and provide more appropriate equations to account for the drop size – flow depth interactions that affect detachment and transport of particles in rain‐impacted flows. There is a need to develop modeling approaches that rely on relevant data obtained under well‐controlled flow conditions where flow depths and velocities are known. Copyright © 2012 John Wiley & Sons, Ltd.