Interrill soil erosion processes and their interaction on low slopes

Soil erosion by water is mostly the result of rainfall‐driven and runoff‐driven processes taking place simultaneously during a storm event. However, the effect of interaction between these two erosion processes has received limited attention. Most laboratory experiments indicate that the rate of erosion in a rain‐impacted flow is greater than for un‐impacted flows of similar depth and velocity; however, negative interaction between the two processes has also been reported. There is no provision for any such interaction in any of the current erosion models. This paper reports on the results of a number of exact experiments on three soil types carried out in the flume of Griffith University's large rainfall simulator to study interaction between rain and runoff processes. The results show that interaction is generally positive under approximately steady state condition and there is very limited sign of negative interaction reported by others. Results provide strong evidence that raindrops continuously peel fine sediment from larger stable aggregates. This mechanism could be the reason for positive interaction during simultaneous rainfall and flow driven erosion in well aggregated soils as a result of increased fine particles in the eroded sediment. Strong positive interaction between rain and runoff erosion also occurs for medium to large aggregates. This strongly suggests that mechanisms that are not well understood are operational. It is quite possible that particle movement can be stimulated by rolling or creeping in a size‐selective manner. Indeed, such additional mechanisms may well be largely responsible for the positive interaction observed between rain and surface flow. Copyright © 2006 John Wiley & Sons, Ltd.     

Simulations demonstrating interaction between coarse and fine sediment loads in rain‐impacted flow

Rain‐impacted flows dominate sheet and interrill erosion and are important in eroding soil rich in nutrients and other chemicals which may have deleterious effects on water quality. Erosion in rain‐impacted flow is associated with raindrop detachment followed by transport either by the combination of flow velocity and raindrop impact (raindrop‐induced flow transport, RIFT) or the inherent capacity of the flow to transport detached material. Coarse particles tend to be transported by RIFT, while fine particles tend to be transported without any assistance from raindrop impact. Because the transport process associated with coarse particles is not 100 per cent efficient, it generates a layer of loose particles on the soil surface and this layer protects the underlying soil from detachment. Simulations were performed by modelling the uplift and downstream movement of both fine and coarse particles detached from the soil surface by individual raindrop impacts starting with a surface where no loose material was present. The simulations produced a flush of fine material followed by a decline in the discharge of fine material as the amount of loose material built up on the bed. The decline in the discharge of fine material was accompanied by an increase in the discharge of coarse material. The relative amounts of coarse and fine material discharged in the flow varied with flow velocity and cohesion in the surface of the soil matrix. The results indicate that the discharge of various sized sediments is highly dependent on local soil, rain and flow conditions and that extrapolating the results from one situation to another may not be appropriate. 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.     

Particle travel distances and bed and sediment compositions associated with rain‐impacted flows

Laboratory experiments were performed with rain of uniform drop size (2·7 mm, 5·1 mm) impacting flows over non‐cohesive beds of uniform sized sand (0·11–0·9 mm) and coal (0·2–0·9 mm) particles with flow velocities (20 mm s^?1, 40 mm s^?1) that were insufficient for the flow to entrain the particles without the aid of raindrop impact. Measurement of particle travel distance under rain made up of 2·7 mm drops confirmed a theoretical relationship between settling velocity and the distance particles travel after being disturbed by drop impact. Although, in theory, a relationship between settling velocity and particle travel distance exists, settling velocity by itself was unable to account for the effect of changes in both particle size and density on sediment discharge from beds of uniform non‐cohesive material. Particle density was also a factor. Further study of how particle characteristics influence sediment discharge will aid modelling of the impact of the soil in process‐based models of erosion by rain‐impacted flow. Copyright © 2001 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.     

Bed scour by debris flows: experimental investigation of effects of debris‐flow composition

Debris flows can grow greatly in size by entrainment of bed material, enhancing their runout and hazardous impact. Here, we experimentally investigate the effects of debris‐flow composition on the amount and spatial patterns of bed scour and erosion downstream of a fixed to erodible bed transition. The experimental debris flows were observed to entrain bed particles both grain by grain and en masse, and the majority of entrainment was observed to occur during passage of the flow front. The spatial bed scour patterns are highly variable, but large‐scale patterns are largely similar over 22.5–35° channel slopes for debris flows of similar composition. Scour depth is generally largest slightly downstream of the fixed to erodible bed transition, except for clay‐rich debris flows, which cause a relatively uniform scour pattern. The spatial variability in the scour depth decreases with increasing water, gravel (= grain size) and clay fraction. Basal scour depth increases with channel slope, flow velocity, flow depth, discharge and shear stress in our experiments, whereas there is no correlation with grain collisional stress. The strongest correlation is between basal scour and shear stress and discharge. There are substantial differences in the scour caused by different types of debris flows. In general, mean and maximum scour depths become larger with increasing water fraction and grain size, and decrease with increasing clay content. However, the erodibility of coarse‐grained experimental debris flows (gravel fraction = 0.64) is similar on a wide range of channel slopes, flow depths, flow velocities, discharges and shear stresses. This probably relates to the relatively large influence of grain‐collisional stress to the total bed stress in these flows (30–50%). The relative effect of grain‐collisional stress is low in the other experimental debris flows (<5%), causing erosion to be largely controlled by basal shear stress. Copyright © 2016 John Wiley & Sons, Ltd.     

Meltwater erosion process of frozen soil as affected by thawed depth under concentrated flow in high altitude and cold regions

Changes in thawed depth of frozen soil caused by diurnal and seasonal temperature fluctuations are commonly found in high altitude and latitude regions of the world. These changes significantly influence hydrologic and erosion processes. Experimental data are necessary to improve the understanding and modeling of the phenomenon. Laboratory experiments were conducted in Beijing to assess the impacts of thawed soil depth, slope gradient, and flow rate on soil erosion by concentrated meltwater flow over an underlying frozen soil layer. Soil samples from watershed were filled in flumes, saturated before being frozen. After the soil was completely frozen, flumes were taken out of storage to thaw the frozen soil from top to the designed depths. Meltwater flow was simulated using a tank filled with water and icecubes at approximately 0°C. The erosion experiments involved four thawed soil depths of 1, 2, 5, and 10 cm; three slope gradients of 5°, 10°, and 15°; and three flow rates of 1, 2, and 4 l/min; and seven rill lengths of 0.5, 1, 2, 3, 4, 5, and 6 m. Sediment‐laden water samples were collected at the lower end of the flume for determination of sediment concentration. The results showed that sediment concentration increased exponentially with rill length to approach a maximum value. The sediment concentrations were closely correlated with thawed soil depth, flow rate, and slope gradient. Shallower thawed depths delivered more sediments than deeper thawed depths. Slope gradient was the primary factor responsible for severe erosion. The effect of flow rate on sediment concentration which decreased with increasing slope gradient, was not as significant as that of slope gradient. Results from these experiments are useful for understanding the effect of thawed soil depth on erosion process in thawed soils subject to freezing and for estimating erosion model parameters. Copyright © 2017 John Wiley & Sons, Ltd.     

Modified MMF (Morgan–Morgan–Finney) model for evaluating effects of crops and vegetation cover on soil erosion

Modifications are made to the revised Morgan–Morgan–Finney erosion prediction model to enable the effects of vegetation cover to be expressed through measurable plant parameters. Given the potential role of vegetation in controlling water pollution by trapping clay particles in the landscape, changes are also made to the way the model deals with sediment deposition and to allow the model to incorporate particle‐size selectivity in the processes of erosion, transport and deposition. Vegetation effects are described in relation to percentage canopy cover, percentage ground cover, plant height, effective hydrological depth, density of plant stems and stem diameter. Deposition is modelled through a particle fall number, which takes account of particle settling velocity, flow velocity, flow depth and slope length. The detachment, transport and deposition of soil particles are simulated separately for clay, silt and sand. Average linear sensitivity analysis shows that the revised model behaves rationally. For bare soil conditions soil loss predictions are most sensitive to changes in rainfall and soil parameters, but with a vegetation cover plant parameters become more important than soil parameters. Tests with the model using field measurements under a range of slope, soil and crop covers from Bedfordshire and Cambridgeshire, UK, give good predictions of mean annual soil loss. Regression analysis of predicted against observed values yields an intercept value close to zero and a line slope close to 1·0, with a coefficient of efficiency of 0·81 over a range of values from zero to 38·6 t ha^?1. Copyright © 2007 John Wiley & Sons, Ltd.     

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.     

Effects of water quality on infiltration,runoff and interrill erosion processes during simulated rainfall

This paper discusses the effects of water quality on the hydrological and erosion response of non‐saline, non‐sodic soils during simulated rain experiments. It is well known that rain water quality affects the behaviour of saline soils. In particular, rain simulation experiments cannot be run using tap water if realistic values of infiltration rates and soil erosion are to be found. This paper reports on similar effects for non‐saline, non‐sodic soils. Two soils – a well‐aggregated clay‐rich soil developed on marine silty clay deposits and a soil developed on silt loam – were selected and subjected to a series of simulated rainstorms using demineralized water and tap water. The experiments were conducted in two different laboratories in order to obtain results independent of the tap water quality or the rainfall simulator characteristics. The results indicate that time‐to‐ponding is largely delayed by solute‐rich water (tap water). When tap water is used, infiltration rates are significantly overestimated, i.e. by more than 100 per cent. Interrill erosion rates increase by a factor of 2·5–3 when demineralized water is used. The silty clay soil was more affected by the water quality than the silt loam soil, with respect to infiltration and runoff production. Regarding interrill erosion rates, the two tested soils were similarly affected by the water quality. Therefore, it can be concluded that rainfall simulation experiments with non‐dispersive soils (e.g. non‐saline, non‐sodic) must also be conducted using water with very low electrical conductivity (i.e. less than 30–50 µS cm⁻¹), close to that of distilled water. The use of tap water certainly hampers comparisons and the relative ranking of the hydrological and erosion response of different soils, while parameter values, such as final infiltration rate or time‐to‐ponding, cannot be extrapolated and extended to natural situations. Therefore, the majority of hydrological and erosion models and parameter values measured during rainfall simulations in the past should be used with caution for all types of soils. Copyright © 2001 John Wiley & Sons, Ltd.     

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.     

Field measurements of interrill erosion under ­different slopes and plot sizes

Despite numerous studies, the effect of slope on interrill erosion is not clearly established. Several interactions exist between erosion parameters that are not taken into account under experimental laboratory measurements and results need to be validated in the field. The influence of slope steepness (2 to 8 per cent) on soil loss for a crusted interrill area and the detachment and transport processes involved in the interaction between slope, rain characteristics and plot size were investigated. Sediment discharge and runoff rates were measured in bounded plots (1 m2 and 10 m2) under natural and simulated rainfall, allowing the analysis of a combination of detachment and transport processes at various scales in the field. Runoff rate increased from 20 to 90 per cent with increasing slope and rain intensity for both plot sizes, whereas sediment concentration increased from 2 to 6 g l−1 with increasing slope only for the 10 m2 plots. At the 1 m2 scale, erosion was transport‐limited due to the reduced rain‐impacted flow. Interactions between slope angle and rain intensity were observed for detachment and transport processes in interrill erosion. Results show the importance of an adapted experimental set‐up to get reference data for interrill erosion model development and validation. Copyright © 2000 John Wiley & Sons, Ltd.     

The increase of rainfall erosivity and initial soil erosion processes due to rainfall acidification

The drastic growth of population in highly industrialized urban areas, as well as fossil fuel use, is increasing levels of airborne pollutants and enhancing acid rain. In rapidly developing countries such as Iran, the occurrence of acid rain has also increased. Acid rain is a driving factor of erosion due to the destructive effects on biota and aggregate stability; however, little is known about its impact on specific rates of erosion at the pedon scale. Thus, the present study aimed to investigate the effect of acid rain at pH levels of 5.25, 4.25, and 3.75 for rainfall intensities of 40, 60, and 80 mm h^?1 on initial soil erosion processes under dry and saturated soil conditions using rainfall simulations. The results were compared using a two‐way ANOVA and Duncan tests and showed that initial soil erosion rates with acidic rain and non‐acidic rain under dry soil conditions were significantly different. The highest levels of soil particle loss due to splash effects in all rainfall intensities were observed with the most acidic rain (pH = 3.75), reaching maximum values of 16 g m^?2 min^?1. The lowest levels of particle losses were observed in the control plot where non‐acidic rain was used, with values ranging from 3.8 to 8.1 g m^?2 min^?1. Similarly, under saturated soil conditions, the lowest level of soil particle loss was observed in the control plot, and the highest peaks of soil loss were observed for the most acidic rains (pH = 3.75 and pH = 4.25), reaching maximum average values of 40 g m^?2 min^?1. However, for saturated soils with acidic water but with non‐acidic rain, the highest soil particle loss was observed for the control plot for all the rainfall intensities. In conclusion, acidic rain has a negative impact on soils, which can be more intense with a concomitant increase in rainfall intensity. Rapid solutions, therefore, need to be found to reduce the emission of pollutants into the air, otherwise, rainfall erosivity may drastically increase.     

One-dimensional morphodynamic model for retrogressive erosion based on a sediment entrainment theory at high flow velocity

Retrogressive erosion is a high-speed erosion process that usually occurs during the rapid release of stored water in reservoirs built on sandy rivers.Retrogressive erosion has been utilized in the practice of reservoir sedimentation control,but accurate prediction of the bed deformation process by numerical models has rarely been reported.The current study presents a one-dimensional morphodynamic model for simulating the evolution process of retrogressive erosion induced by high-velocity flows on steep slopes.The governing equations apply a Cartesian coordinate system with a vertically oriented z axis.The bed surface gradient and friction terms in the flow equations include correction factors to take account of the effects of high slope on flow movement.The net vertical sediment flux term in the sediment transport and bed deformation equations is calculated using an equation of erosion velocity.Particularly,this equation is based on an empirical relation between the sediment entrainment rate and the Shields parameter in contrast to the traditional sediment transport capacity,and the critical Shields parameter is modified by taking into account the permeability of the sediment layer and the stability of particles on a slope.The feedback of scoured sediment on the flow movement is considered by additional terms in the governing equations.Flume experiments of retrogressive erosion in literature were simulated to validate the model.The temporal variations of the longitudinal profiles of the free surface and channel bed and the sediment transport rate were well predicted.The algorithm calculating sediment entrainment in the proposed model also was validated for an experiment measuring entrainment rate from the literature.More importantly,it was found that the morphodynamic model using the sediment transport capacity equation predicts the trend of cumulative erosion contrary to the measurements,while results of the proposed model can follow a similar trend with the observed data in the retrogressive erosion process.     

Terrestrial laser scanning soil surfaces: a field methodology to examine soil surface roughness and overland flow hydraulics

Conventional roughness–resistance relationships developed for pipe and open‐channel flows cannot accurately describe shallow overland flows over natural rough surfaces. This paper develops a new field methodology combining terrestrial laser scanning (TLS) and overland flow simulation to provide a high‐resolution dataset of surface roughness and overland flow hydraulics as simulated on natural bare soil surfaces. This method permits a close examination of the factors controlling flow velocity and a re‐evaluation of the relationship between surface roughness and flow resistance. The aggregate effect of flow dynamics, infiltration and depression storage on retarding the passage of water over a surface is important where runoff‐generating areas are distant from well‐defined channels. Experiments to separate these effects show that this ‘effective resistance’ is dominated by surface roughness. Eight measurements of surface roughness are found to be related to flow resistance: standard deviation of elevations, inundation ratio, pit density (measured both perpendicular and parallel to the flow direction), slope, median depth, skewness of the depth distribution and frontal area. Hillslope position is found to affect the significant roughness measures. In contrast, infiltration rate has little effect on the velocity of water fronts advancing over the soil surfaces examined here and the effect of depression storage is limited. Overland flow resistance is depth dependent where complex microtopographic structures are progressively inundated. Copyright © 2010 John Wiley & Sons, Ltd.     

MODELING HEADCUT DEVELOPMENT AND MIGRATION IN UPLAND CONCENTRATED FLOWS

On hillslopes and agricultural fields, discrete areas of intense, localized soil erosion commonly take place in the form of migrating headcuts. These erosional features significantly increase soil loss and landscape degradation, yet the unsteady, transient, and migratory habits of headcuts complicate their phenomenological and erosional characterization. Here a unique experimental facility was constructed to examine actively migrating headcuts typical of upland concentrated flows. Essential components of the facility include a deep soil cavity with external drainage, rainfall simulator, capacity for overland flow, and a video recording technique for data collection. Results from these experiments show that： （1） after a short period of adjustment, headcut migration attained a steady-state condition, where the rate of migration, scour hole geometry, and sediment discharge remain constant with time; （2） boundary conditions of higher rates of overland flow, steeper bed slopes, and larger initial headcut heights produced systematically larger scour holes with higher rates of soil erosion; and （3） during migration, the turbulent flow structure within the scour hole remained unchanged, consisting of an overfall nappe at the brink transitioning into a reattached wall jet with two recirculation eddies within the plunge pool. The systematic behavior of headcut development and migration enabled the application of modified jet impingement theory to predict with good success the characteristics of the impinging jet, the depth of maximum scour, the rate of headcut migration, and the rate of sediment erosion. These laboratory data and the analytical formulation can be used in conjunction with soil erosion prediction technology to improve the management of agricultural areas impacted by headcut development and ephemeral gully erosion.     

Sub‐grid scale parameterization of hillslope runoff and erosion processes for catchment‐scale models of semi‐arid landscapes

The processes of hillslope runoff and erosion are typically represented at coarse spatial resolution in catchment‐scale models due to computational limitations. Such representation typically fails to incorporate the important effects of topographic heterogeneity on runoff generation, overland flow, and soil erosion. These limitations currently undermine the application of distributed catchment models to understand the importance of thresholds and connectivity on hillslope and catchment‐scale runoff and erosion, particularly in semi‐arid environments. This paper presents a method for incorporating high‐resolution topographic data to improve sub‐grid scale parameterization of hillslope overland flow and erosion models. Results derived from simulations conducted using a kinematic wave overland flow model at 0.5 m spatial resolution are used to parameterize the depth–discharge relationship in the overland flow model when applied at 16 m resolution. The high‐resolution simulations are also used to derive a more realistic parameterization of excess flow shear stress for use in the 16 m resolution erosion model. Incorporating the sub‐grid scale parameterization in the coarse‐resolution model (16 m) leads to improved predictions of overland flow and erosion when evaluated using results derived from high‐resolution (0.5 m) model simulations. The improvement in performance is observed for a range of event magnitudes and is most notable for erosion estimates due to the non‐linear dependency between the rates of erosion and overland flow. Copyright © 2013 John Wiley & Sons, Ltd.     

Sediment yield modelling for small agricultural catchments: land‐cover parameterization based on remote sensing data analysis

Vegetation and soil properties and their associated changes through time and space affect the various stages of soil erosion. The island of Ishigaki in Okinawa Prefecture, Japan is of particular concern because of the propensity of the red‐soil‐dominated watersheds in the area to contribute substantial sediment discharge to adjacent coastal areas. This paper discusses the application of remote sensing techniques in the retrieval of vegetation and soil parameters necessary for the distributed soil‐loss modelling in small agricultural catchments and analyses the variation in erosional patterns and sediment distribution during rainfall events using numerical solutions of overland flow simulations and sediment continuity equations. To account for the spatial as well as temporal variability of selected parameters of the soil‐loss equations, a method is proposed to account for the variability of associated vegetation cover based on their spectral characteristics as captured by remotely sensed data. To allow for complete spatial integration, modelling the movement of sediment is accomplished under a loose‐coupled GIS computational framework. This study lends a theoretical support and empirical evidence to the role of vegetation as a potential agent for soil erosion control. Copyright © 2003 John Wiley & Sons, Ltd.     

Bistable plant–soil dynamics and biogenic controls on the soil production function

Soil formation results from opposite processes of bedrock weathering and erosion, whose balance may be altered by natural events and human activities, resulting in reduced soil depth and function. The impacts of vegetation on soil production and erosion and the feedbacks between soil formation and vegetation growth are only beginning to be explored quantitatively. Since plants require suitable soil environments, disturbed soil states may support less vegetation, leading to a downward spiral of increased erosion and decline in ecosystem function. We explore these feedbacks with a minimal model of the soil–plant system described by two coupled nonlinear differential equations, which include key feedbacks, such as plant‐driven soil production and erosion inhibition. We show that sufficiently strong positive plant–soil feedback can lead to a ‘humped’ soil production function, a necessary condition for soil depth bistability when erosion is assumed to vary monotonically with vegetation biomass. In bistable plant–soil systems, the sustainable soil condition engineered by plants is only accessible above a threshold vegetation biomass and occurs in environments where the high potential rate of erosion exerts a strong control on soil production and erosion. Vegetation removal for agriculture reduces the stabilizing effect of vegetation and lowers the system resilience, thereby increasing the likelihood of transition to a degraded soil state. Copyright © 2016 John Wiley & Sons, Ltd.