A one‐dimensional model for simulating armouring and erosion on hillslopes: 2. Long term erosion and armouring predictions for two contrasting mine spoils

This paper investigates the dynamics of soil armouring as a result of fluvial erosion for a non‐cohesive sandy gravel spoil from the Ranger Mine, Australia, and a cohesive silt loam spoil from the Northparkes Mine, Australia, using a model for hillslope soil armouring. These long term predictions concentrate on the temporal and spatial changes of the spoil grading and erosion over 100–200 years for the flat cap regions (1–2%) and steep batter edges (10–30%) typically encountered on waste rock dumps. The existence of a significant rock fragment fraction in the Ranger spoil means that it armours readily, while Northparkes does not. For Ranger the waste rock showed reductions in (1) cumulative erosion of up to 81% from that obtained by extrapolating the initial erosion rate out 100 years and (2) the erosion/year by more than 10‐fold. For Northparkes reductions were less marked, with the maximum reduction in erosion/year being 37% after 200 years. For Ranger the reductions were greatest and fastest for intermediate gradient hillslopes. For the steepest hillslopes the armouring decreased because the flow shear stresses were large enough to mobilize all material in the armour layer. Model uncertainty was assessed with probabilistic confidence limits demonstrating that these erodibility reductions were statistically significant. A commonly used hillslope erosion model (sediment flux = β₁ discharge m₁ slope n₁) was fitted to these predictions. The erodibility, β₁, and m₁ decreased with time, which was consistent with our physical intuition about armouring. At Ranger the parameter m₁ asymptoted to 1·5–1·6 while at Northparkes it asymptoted to 1·2–1·3. At Ranger transient spatial trends in armouring led to a short term (50–200 years in the future) reduction in n₁, to below zero under certain circumstances, recovering to an asymptote of about 0·5–1. At Northparkes n₁ asymptoted to about 0·6, with no negative transients predicted. The m₁ and n₁ parameters predicted for Ranger were shown to be consistent with field data from a 10‐year‐old armoured hillslope and consistent with published relationships between erodibility and rock content for natural hillslopes. Copyright © 2007 John Wiley & Sons, Ltd.     

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.     

A simulation model for unified interrill erosion and rill erosion on hillslopes

A mathematical model was developed for simulating runoff generation and soil erosion on hillslopes. The model is comprised of three modules: one for overland flow, one for soil infiltration, and one for soil erosion including rill erosion and interrill erosion. Rainfall and slope characteristics affecting soil erosion on hillslopes were analysed. The model results show that the slope length and gradient, time distribution rainfall, and distribution of rills have varying influence on soil erosion. Erosion rate increases nonlinearly with increase in the slope length; a long slope length leads to more serious erosion. The effect of the slope gradient on soil erosion can be both positive and negative. Thus, there exists a critical slope gradient for soil erosion, which is about 45° for the rate of erosion at the end of the slope and about 25° for the accumulated erosion. Copyright © 2005 John Wiley & Sons, Ltd.     

Impact of model simplifications on soil erosion predictions: application of the GLUE methodology to a distributed event‐based model at the hillslope scale

In this paper, we analyse how the performance and calibration of a distributed event‐based soil erosion model at the hillslope scale is affected by different simplifications on the parameterizations used to compute the production of suspended sediment by rainfall and runoff. Six modelling scenarios of different complexity are used to evaluate the temporal variability of the sedimentograph at the outlet of a 60 m long cultivated hillslope. The six scenarios are calibrated within the generalized likelihood uncertainty estimation framework in order to account for parameter uncertainty, and their performance is evaluated against experimental data registered during five storm events. The Nash–Sutcliffe efficiency, percent bias and coverage performance ratios show that the sedimentary response of the hillslope in terms of mass flux of eroded soil can be efficiently captured by a model structure including only two soil erodibility parameters, which control the rainfall and runoff production of suspended sediment. Increasing the number of parameters makes the calibration process more complex without increasing in a noticeable manner the predictive capability of the model. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

A transport‐distance approach to scaling erosion rates: 1. Background and model development

The process basis of existing soil‐erosion models is shown to be ill‐founded. The existing literature builds directly or indirectly on Bennett's (1974) paper, which provided a blueprint for integrated catchment‐scale erosion modelling. Whereas Bennett recognized the inherent assumptions of the approach suggested, subsequent readings of the paper have led to a less critical approach. Most notably, the assumption that sediment movement could be approximated by a continuity equation that related to transport in suspension has produced a series of submodels that assume that all movement occurs in suspension. For commonly occurring conditions on hillslopes, this case is demonstrably untrue both on theoretical grounds and from empirical observations. Elsewhere in the catchment system, it is only partially true, and the extent to which the assumption is reasonable varies both spatially and temporally. A second ground‐breaking paper – that of Foster and Meyer (1972) – was responsible for subsequent uncritical application of a first‐order approximation to deposition based on steady‐state analysis and again a weak empirical basis. We describe in this paper an alternative model (Mahleran – Model for Assessing Hillslope‐Landscape Erosion, Runoff And Nutrients) based upon particle‐travel distance that overcomes existing limitations by incorporating parameterizations of the different detachment and transport mechanisms that occur in water erosion in hillslopes and small catchments. In the second paper in the series, we consider the sensitivity and general behaviour of Mahleran , and test it in relation to data from a large rainfall‐simulation experiment. The third paper of the sequence evaluates the model using data from plots of different sizes in monitored rainfall events. From this evaluation, we consider the scaling characteristics of the current form of Mahleran and suggest that integrated modelling, laboratory and field approaches are required in order to advance the state of the art in soil‐erosion modelling. Copyright © 2008 John Wiley & Sons, Ltd.     

Quantifying rainfall controls on catchment‐scale landslide erosion in Taiwan

Landslide erosion is a dominant hillslope process and the main source of stream sediment in tropical, tectonically active mountain belts. In this study, we quantified landslide erosion triggered by 24 rainfall events from 2001 to 2009 in three mountainous watersheds in Taiwan and investigated relationships between landslide erosion and rainfall variables. The results show positive power‐law relations between landslide erosion and rainfall intensity and cumulative rainfall, with scaling exponents ranging from 2·94 to 5·03. Additionally, landslide erosion caused by Typhoon Morakot is of comparable magnitude to landslide erosion caused by the Chi‐Chi Earthquake (M_W = 7·6) or 22–24 years of basin‐averaged erosion. Comparison of the three watersheds indicates that deeper landslides that mobilize soil and bedrock are triggered by long‐duration rainfall, whereas shallow landslides are triggered by short‐duration rainfall. These results suggest that rainfall intensity and watershed characteristics are important controls on rainfall‐triggered landslide erosion and that severe typhoons, like high‐magnitude earthquakes, can generate high rates of landslide erosion in Taiwan. Copyright © 2012 John Wiley & Sons, Ltd.     

Implications of model uncertainty for the mapping of hillslope‐scale soil erosion predictions

An integrated modelling approach (MIRSED) which utilizes the process‐based soil erosion model WEPP (Water Erosion Prediction Project) is presented for the assessment of hillslope‐scale soil erosion at five sites throughout England and Wales. The methodology draws upon previous uncertainty analysis of the WEPP hillslope soil erosion model by the authors to qualify model results within an uncertainty framework. A method for incorporating model uncertainty from a range of sources is discussed as a first step towards using and learning from results produced through the GLUE (Generalized Likelihood Uncertainty Estimation) technique. Results are presented and compared to available observed data, which illustrate that levels of uncertainty are significant and must be taken into account if a meaningful understanding of output from models such as WEPP is to be achieved. Furthermore, the collection of quality, observed data is underlined for two reasons: as an essential tool in the development of soil erosion modelling and also to allow further constraint of model uncertainty. Copyright © 2001 John Wiley & Sons, Ltd.     

Amorphous silica mobilization by inter‐rill erosion: insights from rainfall experiments

Amorphous silica (ASi) carried in suspension by rivers is an important component in the global Si budget. Water erosion processes in cultivated catchments are likely to drive ASi delivery to the river system. However, no studies have investigated the controls on ASi mobilization by erosional processes in croplands. Rainfall experiments were performed on split fields (i.e. a part conventionally ploughed and a part under reduced tillage) to simulate ASi mobilization by inter‐rill erosion in croplands, and identify its dependency on soil, field and rainfall characteristics. The ASi content of the soil and the inter‐rill erosion rate were determined as the major controls on ASi mobilization. Variables such as tillage technique and crop type did not have a consistent direct or indirect effect. Inter‐rill erosion is clearly selective with respect to ASi, indicating association of ASi with the fine soil fraction and with soil organic carbon. Our experiments demonstrate that erosion increases due to human perturbation will increase the delivery of reactive Si to aquatic systems. We estimate that globally, c. 7% of all reactive Si that enters aquatic systems is derived from erosion of agricultural soils. Copyright © 2015 John Wiley & Sons, Ltd.     

Deforestation of water‐repellent soils in Galicia (NW Spain): effects on surface runoff and erosion under simulated rainfall

At three adjacent sites in steeply sloping woodland in Galicia (NW Spain), surface runoff and associated erosion under simulated rainfall (64 mm h^?1) were measured on five occasions between June 1998 and July 1999. Two of the three sites had recently been deforested and topsoil added, and one of these two had been sown with grass, which was germinating at the onset of the study. Deforestation greatly increased runoff and erosion rates, and the recovery of plant cover reduced erosion. All three soils were very hydrophobic due to high levels of poorly humified organic matter, which led to higher runoff rates than expected, especially during dry periods. However, great structural stability prevented there being a significant correlation between runoff rate and soil erosion. Copyright © 2003 John Wiley & Sons, Ltd.     

Using sediment travel distance to estimate medium‐term erosion rates: a 16‐year record

Fine grained (80 µm) magnetite was introduced onto a semi‐arid grassland hillslope in 1992, as part of a set of rainfall‐simulation experiments. Using measurements of magnetic susceptibility, the median distance travelled by these magnetite grains during subsequent natural runoff events in the 16‐year period up to 2008 was estimated. Coupling this estimate to direct measurements of sediment flux obtained during the rainfall‐simulation experiments has enabled estimation of the erosion rate over this period. The estimated average erosion rate of between 2·61 × 10^?2 and 4·36 × 10^?2 kg m^?1 year^?1, is equivalent to a rate of ground lowering between 0·020 and 0·033 mm year^?1. This estimate is consistent with (in the sense of being less than) an estimate of total sediment detachment over the same period. The rate of erosion measured using this travel‐distance approach is an order of magnitude less that obtained from a study based on ¹³⁷Cs in a nearby catchment, and compatible with the longevity of continents. Copyright © 2010 John Wiley & Sons, Ltd.     

A low‐dimensional model of bedrock weathering and lateral flow coevolution in hillslopes: 1. Hydraulic theory of reactive transport

This is the first of a two‐part paper exploring the coevolution of bedrock weathering and lateral flow in hillslopes using a simple low‐dimensional model based on hydraulic groundwater theory (also known as Dupuit or Boussinesq theory). Here, we examine the effect of lateral flow on the downward fluxes of water and solutes through perched groundwater at steady state. We derive analytical expressions describing the decline in the downward flux rate with depth. Using these, we obtain analytical expressions for water age in a number of cases. The results show that when the permeability field is homogeneous, the spatial structure of water age depends qualitatively on a single dimensionless number, Hi. This number captures the relative contributions to the lateral hydraulic potential gradient of the relief of the lower‐most impermeable boundary (which may be below the weathering front within permeable or incipiently weathered bedrock) and the water table. A “scaled lateral symmetry” exists when Hi is low: age varies primarily in the vertical dimension, and variations in the horizontal dimension x almost disappear when the vertical dimension z is expressed as a fraction z/H(x) of the laterally flowing system thickness H(x). Taking advantage of this symmetry, we show how the lateral dimension of the advection–diffusion‐reaction equation can be collapsed, yielding a 1‐D vertical equation in which the advective flux downward declines with depth. The equation holds even when the permeability field is not homogeneous, as long as the variations in permeability have the same scaled lateral symmetry structure. This new 1‐D approximation is used in the accompanying paper to extend chemical weathering models derived for 1‐D columns to hillslope domains.     

Rill development and soil erosion: a laboratory study of slope and rainfall intensity

A total of 15 rainfall simulation experiments were conducted in a 1 m by 2 m box varying slope (10, 20, 30%) and rainfall intensity (60, 90, 120 mm h^?1). The experiments were performed to study how rill networks initiate and evolve over time under controlled conditions with regard to the treatment variables considered, and to allow for input in a computer simulation model. Runoff and sediment yield samples were collected. Digital elevation models were calculated by means of photogrammetry for several time steps of most experiments. The soil used in the experiments was a basal till derived Cambisol typical for the Swiss Plateau. While significant differences were found for sediment yield, runoff did not vary significantly with treatment combinations. Increasing rainfall intensity had a larger effect on sediment yield than increasing slope. Rill density and energy expenditure decreased with time, suggesting that energy expenditure was a useful parameter to describe the emergence of rill network at the laboratory scale. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

Insight into sediment transport processes on saline rangeland hillslopes using three‐dimensional soil microtopography changes

In arid and semi‐arid rangeland environments, an accurate understanding of runoff generation and sediment transport processes is key to developing effective management actions and addressing ecosystem response to changes. Yet, many primary processes (namely sheet and splash and concentrated flow erosion, as well as deposition) are still poorly understood due to a historic lack of measurement techniques capable of parsing total soil loss into these primary processes. Current knowledge gaps can be addressed by combining traditional erosion and runoff measurement techniques with image‐based three‐dimensional (3D) soil surface reconstructions. In this study, data (hydrology, erosion and high‐resolution surface microtopography changes) from rainfall simulation experiments on 24 plots in saline rangelands communities of the Upper Colorado River Basin were used to improve understanding on various sediment transport processes. A series of surface change metrics were developed to quantify and characterize various erosion and transport processes (e.g. plot‐wide versus concentrated flow detachment and deposition) and were related to hydrology and biotic and abiotic land surface characteristics. In general, erosivity controlled detachment and transport processes while factors modulating surface roughness such as vegetation controlled deposition. The extent of the channel network was a positive function of slope, discharge and vegetation. Vegetation may deflect runoff in many flow paths but promoted deposition. From a management perspective, this study suggests that effective runoff soil and salt load reduction strategies should aim to promote deposition of transported sediments rather than reducing detachment which might not be feasible in these resource‐limited environments. Copyright © 2016 John Wiley & Sons, Ltd.     

A non‐equilibrium sediment transport model for rill erosion

Sediment transport in rill flows exhibits the characteristics of non‐equilibrium transport, and the sediment transport rate of rill flow gradually recovers along the flow direction by erosion. By employing the concept of partial equilibrium sediment transport from open channel hydraulics, a dynamic model of rill erosion on hillslopes was developed. In the model, a parameter, called the restoration coefficient of sediment transport capacity, was used to express the recovery process of sediment transport rate, which was analysed by dimensional analysis and determined from laboratory experimental data. The values of soil loss simulated by the model were in agreement with observed values. The model results showed that the length and gradient of the hillslope and rainfall intensity had different influences on rill erosion. Copyright © 2006 John Wiley & Sons, Ltd.     

Rill length and plot‐scale effects on the hydrogeomorphologic response of gravelly roadbeds

Although unpaved roads are well‐recognized as important sources of Hortonian overland flow and sediment in forested areas, their role in agriculturally‐active rural settings still lacks adequate documentation. In this study, we assessed the effect of micro‐catchment size, slope, and ground cover on runoff and sediment generation from graveled roadbeds servicing a rural area in southern Brazil. Fifteen replications based on 30‐min‐long simulated rainfall experiments were performed at constant rainfall intensities of 22–58 mm h^?1 on roadbeds with varying characteristics including ~3–7 m² micro‐catchment areas, 2–11° slopes, 2–9.7‐m‐long shallow rill features, and 30–100% gravel cover. The contributions of micro‐catchment size and rill length were the most important physical characteristics affecting runoff response and sediment production; both the size of the micro‐catchment and the length of the rills were inversely related to sediment loss and this contradicts most of the rill erosion literature. The effect of micro‐catchment size on runoff and sediment response suggests a potentially problematic spatial‐scale subjectivity of experimental plot results. The inverse relationship between rill length and sediment generation is interpreted here as related to the predominance of coarse fragments within rills, the inability of the shallow flows generated during the simulations to erode this sediment, and their role as zones of net sediment storage. Copyright © 2015 John Wiley & Sons, Ltd.     

SCALES: a large‐scale assessment model of soil erosion hazard in Basse‐Normandie (northern‐western France)

The cartography of erosion risk is mainly based on the development of models, which evaluate in a qualitative and quantitative manner the physical reproduction of the erosion processes (CORINE, EHU, INRA). These models are mainly semi‐quantitative but can be physically based and spatially distributed (the Pan‐European Soil Erosion Risk Assessment, PESERA). They are characterized by their simplicity and their applicability potential at large temporal and spatial scales. In developing our model SCALES (Spatialisation d'éChelle fine de l'ALéa Erosion des Sols/large‐scale assessment and mapping model of soil erosion hazard), we had in mind several objectives: (1) to map soil erosion at a regional scale with the guarantee of a large accuracy on the local level, (2) to envisage an applicability of the model in European oceanic areas, (3) to focus the erosion hazard estimation on the level of source areas (on‐site erosion), which are the agricultural parcels, (4) to take into account the weight of the temporality of agricultural practices (land‐use concept). Because of these objectives, the nature of variables, which characterize the erosion factors and because of its structure, SCALES differs from other models. Tested in Basse‐Normandie (Calvados 5500 km²) SCALES reveals a strong predisposition of the study area to the soil erosion which should require to be expressed in a wet year. Apart from an internal validation, we tried an intermediate one by comparing our results with those from INRA and PESERA. It appeared that these models under estimate medium erosion levels and differ in the spatial localization of areas with the highest erosion risks. SCALES underlines here the limitations in the use of pedo‐transfer functions and the interpolation of input data with a low resolution. One must not forget however that these models are mainly focused on an interregional comparative approach. Therefore the comparison of SCALES data with those of the INRA and PESERA models cannot result on a convincing validation of our model. For the moment the validation is based on the opinion of local experts, who agree with the qualitative indications delivered by our cartography. An external validation of SCALES is foreseen, which will be based on a thorough inventory of erosion signals in areas with different hazard levels. Copyright © 2010 John Wiley & Sons, Ltd.     

A one‐dimensional process‐based approach to study reservoir sediment dynamics during management operations

A comprehensive understanding of the dynamics of erosion and sedimentation in reservoirs under different management conditions is required to anticipate sedimentation issues and implement effective sediment management strategies. This paper describes a unique approach combining fluvial geomorphology tools and morphodynamic modeling for analyzing the sediment dynamics of an elongated hydropower reservoir subjected to management operations: the Génissiat Reservoir on the Rhône River. Functional sub‐reaches representative of the reservoir morphodynamics were delineated by adapting natural river segmentation methods to elongated reservoirs. The segmentation revealed the link between the spatial and temporal reservoir changes and the variability of longitudinal flow conditions during reservoir management operations. An innovative modeling strategy, incorporating the reservoir segmentation into two sediment transport codes, was implemented to simulate the dynamics of erosion and sedimentation at the reach scale during historic events. One code used a bedload approach, based on the Exner equation with a transport capacity formula, and the other used a suspended load approach based on the advection–dispersion equation. This strategy provided a fair quantification of the dynamics of erosion and sedimentation at the reach scale during different management operations. This study showed that the reservoir morphodynamics is controlled by bedload transport in upper reaches, graded suspended load transport of sand in middle reaches and suspended load transport of fine sediments in lower reaches. Eventually, it allowed a better understanding of the impact of dam management on sediment dynamics. Copyright © 2017 John Wiley & Sons, Ltd.     

Effects of rainfall and slope on runoff,soil erosion and rill development: an experimental study using two loess soils

Runoff generation and soil loss from slopes have been studied for decades, but the relationships among runoff, soil loss and rill development are still not well understood. In this paper, rainfall simulation experiments were conducted in two neighbouring plots (scale: 1 m by 5 m) with four varying slopes (17.6%, 26.8%, 36.4% and 46.6%) and two rainfall intensities (90 and 120 mm h^?1) using two loess soils. Data on rill development were extracted from the digital elevation models by means of photogrammetry. The effects of rainfall intensity and slope gradient on runoff, soil loss and rill development were different for the two soils. The runoff and soil loss from the Anthrosol surface were generally higher than those from the Calcaric Cambisol surface. Higher rainfall intensity produced less runoff and more sediment for almost each treatment. With increasing slope gradient, the values of cumulative runoff and soil loss peaked, except for the treatments with 90 mm h^?1 rainfall on the slopes with Anthrosol. With rainfall duration, runoff discharge decreased for Anthrosol and increased for Calcaric Cambisol for almost all the treatments. For both soils, sediment concentration was very high at the onset of rainfall and decreased quickly. Almost all the sediment concentrations increased on the 17.6% and 26.8% slopes and peaked on the 36.4% and 46.6% slopes. Sediment concentrations were higher on the Anthrosol slopes than on the Calcaric Cambisol slopes. At 90 mm h^?1 rainfall intensity, increasingly denser rills appeared on the Anthrosol slope as the slope gradient increased, while only steep slopes (36.4% and 46.6%) developed rills for the Calcaric Cambisol soil. The contributions of rill erosion ranged from 36% to 62% of the cumulative soil losses for Anthrosol, while the maximum contribution of rill erosion to the cumulative soil loss was only 37.9% for Calcaric Cambisol. Copyright © 2014 John Wiley & Sons, Ltd.