Physically based modelling of sheet erosion (detachment and deposition processes) in complex hillslopes

A one‐dimensional uncoupled model governed by this research is a physics‐based modelling of the rainfall‐runoff induced erosion process. The presented model is composed of three parts of a three‐dimensional (3D) hillslope geometry, a nonlinear storage (kinematic wave) model for hillslope hydrological response, and an unsteady physically based surface erosion model. The 3D hillslope geometry model allows describing of the hillslope morphology by defining their plan shape and profile curvature. By changing these two topographic parameters, nine basic hillslope types are derived. The modelling of hillslope hydrological response is based on a flow continuity equation as the relation of discharge and flow depth is passed on kinematic wave approximation. The erosion model is based on a mass conservation equation for unsteady flow. The model assumes that suspended sediment does not affect flow dynamics. The model also accounts for the effect of flow depth plus loose soil depth on soil detachment. The presented model was run for two different precipitations, slope content, and length, and results were plotted for sediment detachment/deposition rate. Based on the obtained results, in hillslopes with convex and straight profile curvatures, sediment detachment only occurred in the whole length of the hillslope. However, in concave ones, sediment detachment and deposition only occurred together in hillslope. The hillslopes with straight profiles and convergent plans have the highest rate of detachment. Also, results show that most detachment rates occur in convex profile curvatures, which are about 15 times more than in straight profiles. Copyright © 2015 John Wiley & Sons, Ltd.     

The Daily Erosion Project – daily estimates of water runoff,soil detachment,and erosion

Water runoff and sediment transport from agricultural uplands are substantial threats to water quality and sustained crop production. To improve soil and water resources, farmers, conservationists, and policy‐makers must understand how landforms, soil types, farming practices, and rainfall interact with water runoff and soil erosion processes. To that end, the Iowa Daily Erosion Project (IDEP) was designed and implemented in 2003 to inventory these factors across Iowa in the United States. IDEP utilized the Water Erosion Prediction Project (WEPP) soil erosion model along with radar‐derived precipitation data and government‐provided slope, soil, and management information to produce daily estimates of soil erosion and runoff at the township scale (93 km² [36 mi²]). Improved national databases and evolving remote sensing technology now permit the derivation of slope, soil, and field‐level management inputs for WEPP. These remotely sensed parameters, along with more detailed meteorological data, now drive daily WEPP hillslope soil erosion and water runoff estimates at the small watershed scale, approximately 90 km² (35 mi²), across sections of multiple Midwest states. The revisions constitute a substantial improvement as more realistic field conditions are reflected, more detailed weather data are utilized, hill slope sampling density is an order of magnitude greater, and results are aggregated based on surface hydrology enabling further watershed research and analysis. Considering these improvements and the expansion of the project beyond Iowa it was renamed the Daily Erosion Project (DEP). Statistical and comparative evaluations of soil erosion simulations indicate that the sampling density is adequate and the results are defendable. The modeling framework developed is readily adaptable to other regions given suitable inputs. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Equifinality and uncertainty in physically based soil erosion models: application of the GLUE methodology to WEPP–the Water Erosion Prediction Project–for sites in the UK and USA

Despite the wealth of soil erosion models available for the prediction of both runoff and soil loss at a variety of scales, little quantification is made of uncertainty and error associated with model output. This in part reflects the need to produce unequivocal or optimal results for the end user, which will often be an unrealistic goal. This paper presents a conceptually simple methodology, Generalized Likelihood Uncertainty Estimation (GLUE), for assessing the degree of uncertainty surrounding output from a physically based soil erosion model, the Water Erosion Prediction Project (WEPP). The ability not only to be explicit about model error but also to evaluate future improvements in parameter estimation, observed data or scientific understanding is demonstrated. This approach is applied to two sets of soil loss/runoff plot replicates, one in the UK and one in the USA. Although it is demonstrated that observations can be largely captured within uncertainty bounds, results indicate that these uncertainty bounds are often wide, reflecting the need to qualify results that derive from ‘optimum’ parameter sets, and to accept the concept of equifinality within soil erosion models. Attention is brought to the problem of under‐prediction of large events/over‐prediction of small events, as an area where model improvements could be made, specifically in the case of relatively dry years. Finally it is proposed that such a technique of model evaluation be employed more widely within the discipline so as to aid the interpretation and understanding of complex model output. Copyright © 2000 John Wiley & Sons, Ltd.     

Hillslope erosion measurement—a simple approach to a complex process

The measurement of hillslope erosion can be a difficult, costly and time‐consuming activity. Many techniques are available, ranging from using environmental tracers, to LiDAR. Erosion measurements using erosion pins are assessed and compared with regional scale erosion data, hillslope data obtained using ¹³⁷Cs and erosion modelling results. The pins produced erosion rates which are within the range determined using ¹³⁷Cs and model data but above that of regional denudation rates. Our findings demonstrate that inexpensive erosion pins can provide reliable data on hillslope erosion. © 2015 Commonwealth of Australia. Hydrological Processes © 2015 John Wiley & Sons Ltd.     

Capabilities and limitations of detailed hillslope hydrological modelling

Hillslope hydrological modelling is considered to be of great importance for the understanding and quantification of hydrological processes in hilly or mountainous landscapes. In recent years a few comprehensive hydrological models have been developed at the hillslope scale which have resulted in an advanced representation of hillslope hydrological processes (including their interactions), and in some operational applications, such as in runoff and erosion studies at the field scale or lateral flow simulation in environmental and geotechnical engineering. An overview of the objectives of hillslope hydrological modelling is given, followed by a brief introduction of an exemplary comprehensive hillslope model, which stimulates a series of hydrological processes such as interception, evapotranspiration, infiltration into the soil matrix and into macropores, lateral and vertical subsurface soil water flow both in the matrix and preferential flow paths, surface runoff and channel discharge. Several examples of this model are presented and discussed in order to determine the model's capabilities and limitations. Finally, conclusions about the limitations of detailed hillslope modelling are drawn and an outlook on the future prospects of hydrological models on the hillslope scale is given.The model presented performed reasonable calculations of Hortonian surface runoff and subsequent erosion processes, given detailed information of initial soil water content and soil hydraulic conditions. The vertical and lateral soil moisture dynamics were also represented quite well. However, the given examples of model applications show that quite detailed climatic and soil data are required to obtain satisfactory results. The limitations of detailed hillslope hydrological modelling arise from different points: difficulties in the representations of certain processes (e.g. surface crusting, unsaturated–saturated soil moisture flow, macropore flow), problems of small‐scale variability, a general scarcity of detailed soil data, incomplete process parametrization and problems with the interdependent linkage of several hillslopes and channel–hillslope interactions. Copyright © 1999 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.     

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.     

Predicting soil erosion hazard in Lattakia Governorate(W Syria)

The main objective of this study is to predict soil erosion in the Lattakia Governorate(WSyria)using the Water Erosion Prediction Project model(WEPP)and to compare the result with that of the RUSLE.Field survey and data collection were carried out,and 44 soil samples were analyzed.In addition,all the necessary input files were prepared for use in the WEPP model and RUSLE.Results show that more than of 80%of the locations studied experience slight to moderate erosion(less than 5 t/ha/y),whereas the rest of the locations experience severe soil erosion hazard.Moreover,the volume of runoff estimated by the WEPP model is in the range of 51e321 mm,and the R^2 between the simulated soil erosion and the predicted runoff reached 0.68.Interestingly,the R^2 between the WEPP model and RUSLE is 0.56,which indicates a good correlation between the two models.     

Effects of gradient,distance, curvature and aspect on steep burned and unburned hillslope soil erosion and deposition

Wildfire denudes vegetation and impacts chemical and physical soil properties, which can alter hillslope erosion rates. Post‐wildfire erosion can also contribute disproportionately to long‐term erosion rates and landscape evolution. Post‐fire hillslope erosion rates remain difficult to predict and document at the hillslope scale. Here we use ²¹⁰Pb_aex (lead‐210 mineral‐adsorbed excess) inventories to describe net sediment erosion on steep, convex hillslopes in three basins (unburned, moderately and severely burned) in mountainous central Idaho. We analyzed nearly 300 soil samples for ²¹⁰Pb_aex content with alpha spectrometry and related net sediment erosion to burn severity, aspect, gradient, curvature and distance from ridgetop. We also tested our data against models for advective, linear and non‐linear diffusive erosion. Statistically lower net soil losses on north‐ versus south‐facing unburned hillslopes suggest that greater vegetative cover and soil cohesion on north‐facing slopes decrease erosion. On burned hillslopes, erosion differences between aspects were less apparent and net erosion was more variable, indicating that vegetation influences erosion magnitude and fire drives erosion variability. We estimated net soil losses throughout the length of unburned hillslopes, including through a footslope transition to concave form. In contrast, on burned hillslopes, the subtle shift from convex to concave form was associated with deposition of a post‐fire erosion pulse. Such overall patterns of erosion and deposition are consistent with predictions from a non‐linear diffusion equation. This finding also suggests that concave sections of overall convex hillslopes affect post‐disturbance soil erosion and deposition. Despite these patterns, no strong relationships were evident between local net soil losses and gradient, curvature, distance from ridgetop, or erosion predicted with advection or diffusion equations. The observed relationship between gradient and erosion is therefore likely more complex or stochastic than often described theoretically, especially over relatively short timescales (60–100 years). Copyright © 2016 John Wiley & Sons, Ltd.     

An analysis of the impact of spatial variability in rainfall on runoff and sediment predictions from a distributed model

This paper investigates the effect of introducing spatially varying rainfall fields to a hydrological model simulating runoff and erosion. Pairs of model simulations were run using either spatially uniform (i.e. spatially averaged) or spatially varying rainfall fields on a 500‐m grid. The hydrological model used was a simplified version of Thales which enabled runoff generation processes to be isolated from hillslope averaging processes. Both saturation excess and infiltration excess generation mechanisms were considered, as simplifications of actual hillslope processes. A 5‐year average recurrence interval synthetic rainfall event typical of temperate climates (Melbourne, Australia) was used. The erosion model was based on the WEPP interrill equation, modified to allow nonlinear terms relating the erosion rate to rainfall or runoff‐squared. The model results were extracted at different scales to investigate whether the effects of spatially varying rainfall were scale dependent. A series of statistical metrics were developed to assess the variability due to introducing the spatially varying rainfall field. At the catchment (approximately 150 km²) scale, it was found that particularly for saturation excess runoff, model predictions of runoff were insensitive to the spatial resolution of the rainfall data. Generally, erosion processes at smaller sub‐catchment scales, particularly when the sediment generation equation had non linearity, were more sensitive to spatial rainfall variability. Introducing runon infiltration reduced the total runoff and sediment yield at all scales, and this process was also most sensitive to the rainfall resolution. Copyright © 2011 John Wiley & Sons, Ltd.     

The sensitivity of hillslope bedrock erosion to precipitation

Decoupling the impacts of climate and tectonics on hillslope erosion rates is a challenging problem. Hillslope erosion rates are well known to respond to changes in hillslope boundary conditions (e.g. channel incision rates) through their dependence on soil thickness, and precipitation is an important control on soil formation. Surprisingly though, compilations of hillslope denudation rates suggest little precipitation sensitivity. To isolate the effects of precipitation and boundary condition, we measured rates of soil production from bedrock and described soils on hillslopes along a semi‐arid to hyperarid precipitation gradient in northern Chile. In each climate zone, hillslopes with contrasting boundary conditions (actively incising channels versus non‐eroding landforms) were studied. Channel incision rates, which ultimately drive hillslope erosion, varied with precipitation rather than tectonic setting throughout the study area. These precipitation‐dependent incision rates are mirrored on the hillslopes, where erosion shifts from relatively fast and biologically‐driven to extremely slow and salt‐driven as precipitation decreases. Contrary to studies in humid regions, bedrock erosion rates increase with precipitation following a power law, from ～1 m Ma^?1 in the hyperarid region to ～40 m Ma^?1 in the semi‐arid region. The effect of boundary condition on soil thickness was observed in all climate zones (thicker soils on hillslopes with stable boundaries compared to hillslopes bounded by active channels), but the difference in bedrock erosion rates between the hillslopes within a climate region (slower erosion rates on hillslopes with stable boundaries) decreased as precipitation decreased. The biotic‐abiotic threshold also marks the precipitation rate below which bedrock erosion rates are no longer a function of soil thickness. Our work shows that hillslope processes become sensitive to precipitation as life disappears and the ability of the landscape to respond to tectonics decreases. Copyright © 2010 John Wiley & Sons, Ltd.     

Sensitivity analysis of EUROSEM using Monte Carlo simulation I: hydrological,soil and vegetation parameters

Knowledge about model uncertainty is essential for erosion modelling and provides important information when it comes to parameterizing models. In this paper a sensitivity analysis of the European soil erosion model (EUROSEM) is carried out using Monte Carlo simulation, suitable for complex non‐linear models, using time‐dependent driving variables. The analysis revealed some important characteristics of the model. The variability of the static output parameters was generally high, with the hydrologic parameters being the most important ones, especially saturated hydraulic conductivity and net capillary drive followed by the percentage basal area for the hydrological and vegetation parameters and detachability and cohesion for the soil erosion parameters. Overall, sensitivity to vegetation parameters was insignificant. The coefficient of variation for the sedigraph was higher than for the hydrograph, especially from the beginning of the rainstorm and up to the peak, and may explain difficulties encountered when trying to match simulated hydrographs and sedigraphs with observed ones. The findings from this Monte Carlo simulation calls for improved within‐storm modelling of erosion processes in EUROSEM. Information about model uncertainty will be incorporated in a new EUROSEM user interface. Copyright © 2000 John Wiley & Sons, Ltd.     

SOIL EROSION PROCESS RESEARCH AND ITS POTENTIAL IMPACT ON EROSION PREDICTION MODEL DEVELOPMENT

During the past 50 years, many research efforts have been invested in understanding soil erosion process and development of erosion prediction models at various scales. This paper briefly introduces the erosion process and prediction model development in the USA. Especially, this paper focuses on discussing potential impacts of the erosion process on erosion model development, and future directions of the soil erosion process research and process- based model development. 1 DEVELOPMENT O…     

Use of caesium‐137 data to evaluate SHETRAN simulated long‐term erosion patterns in arable lands

The caesium‐137 method of quantifying soil erosion is used to provide field data for validating the capability of the SHETRAN modelling system for predicting long‐term (30‐year) erosion rates and their spatial variability. Simulations were carried out for two arable farm sites (area 3–5 ha) in central England for which average annual erosion rates of 6·5 and 10·4 t ha^?1 year^?1 had already been determined using caesium‐137 measurements. These rates were compared with a range of simulated values representing the uncertainty in model output derived from uncertainty in the evaluation of model parameters. A successful validation was achieved in that the simulation range contained the measured rate at both sites, whereas the spatial variability was reproduced excellently at one site and partially at the other. The results indicate that, as the caesium‐137 technique measures the erosion caused by all the processes acting at a site, it is relevant to hydrologically based models such as SHETRAN only if erosion by wind, agricultural activities and other processes not represented in the model are insignificant. The results also indicate a need to reduce the uncertainty in model parameter evaluation. More generally, the caesium‐137 technique is shown to provide field data that improve the severity of the validation procedure (accounting for internal as well as outlet conditions) and that add spatial variability to magnitude as a condition for identifying unrealistic parameter sets when seeking to reduce simulation uncertainty. Copyright © 2004 John Wiley & Sons, Ltd.     

Runoff and water erosion modelling using WEPP on a Mediterranean cultivated catchment

Soil erosion due to water is a major environmental problem in many parts of the world. Most of Mediterranean countries are concerned because of their specific climate and soils sensitivity, but also because of the recent intensification of human activities and agricultural practices. Accurate estimation of soil water erosion for various land-use and climate scenarios is so an important key to define sustainable management policies. In the last decades, several studies have been carried out to build models suitable for quantifying soil erosion. Among these models, the Water Erosion Prediction Project (WEPP, Flanagan, D.C., Nearing, M.A., 1995. USDA-Water Erosion Prediction Project: Hillslope profile and watershed model documentation. NSERL Report 10, USDA-ARS National Soil Erosion Research Laboratory, West Lafayette, IN, USA.) is a physically based, distributed-parameter model that has been developed and mainly validated in USA. Only few studies have investigated its applicability to environmental conditions that differs from those where the model was developed. The aim of this work is to test the efficiency of WEPP model to predict soil erosion at catchment scale in a Mediterranean semi-arid area. Continuous simulations have been conducted between 1995 and 2002 on an cultivated experimental catchment located upstream from a hill reservoir (Kamech catchment, 2.45 km², Cap Bon, Tunisia) where runoff and soil erosion measurements are available at the outlet. Comparison between predictions and measurements shows significant differences. Processes related to seasonal effects (as cracking soils) are pointed out as a weakness of WEPP model for Mediterranean conditions.     

First‐year post‐fire erosion rates in Bitterroot National Forest,Montana

Accelerated runoff and erosion commonly occur following forest fires due to combustion of protective forest floor material, which results in bare soil being exposed to overland flow and raindrop impact, as well as water repellent soil conditions. After the 2000 Valley Complex Fires in the Bitterroot National Forest of west‐central Montana, four sets of six hillslope plots were established to measure first‐year post‐wildfire erosion rates on steep slopes (greater than 50%) that had burned with high severity. Silt fences were installed at the base of each plot to trap eroded sediment from a contributing area of 100 m². Rain gauges were installed to correlate rain event characteristics to the event sediment yield. After each sediment‐producing rain event, the collected sediment was removed from the silt fence and weighed on site, and a sub‐sample taken to determine dry weight, particle size distribution, organic matter content, and nutrient content of the eroded material. Rainfall intensity was the only significant factor in determining post‐fire erosion rates from individual storm events. Short duration, high intensity thunderstorms with a maximum 10‐min rainfall intensity of 75 mm h^?1 caused the highest erosion rates (greater than 20 t ha^?1). Long duration, low intensity rains produced little erosion (less than 0·01 t ha^?1). Total C and N in the collected sediment varied directly with the organic matter; because the collected sediment was mostly mineral soil, the C and N content was small. Minimal amounts of Mg, Ca, and K were detected in the eroded sediments. The mean annual erosion rate predicted by Disturbed WEPP (Water Erosion Prediction Project) was 15% less than the mean annual erosion rate measured, which is within the accuracy range of the model. Published in 2007 by John Wiley & Sons, Ltd.     

Simulation of surface runoff and sediment yield using the water erosion prediction project (WEPP) model: a study in Kaneli watershed,Himalaya, India / Simulation de ruissellement de surface et d'érosion à l'aide du modèle WEPP: cas du bassin versant de Kaneli,Himalaya, Inde

Abstract

The Water Erosion Prediction Project (WEPP) model was calibrated and evaluated for estimation of runoff and sediment yield in the data-scarce conditions of the Indian Himalaya. The inputs derived from remote sensing and geographic information system technologies were combined in the WEPP modelling system to simulate surface runoff and sediment yield from the hilly Kaneli watershed. The model parameters were calibrated using measured data on runoff volumes and sediment yield. The calibrated model was validated by producing the monthly runoff and sediment yield simulations and comparing them with data that were not used in calibration. The model was also used to make surface runoff and sediment yield simulations for each of the individual watershed elements, comprising 18 hillslopes and seven channels, and the detailed monthly results for each are presented. Although, no field data on hillslope runoff and sediment yield are currently available for the validation of distributed results produced by the model, the present investigation has demonstrated clearly the applicability of the WEPP model in predicting hydrological variables in a data-scarce situation.     

Evaluation of the LISEM soil erosion model in two catchments in the East African Highlands

Under increasing population pressure, soil erosion has become a threat in the East African Highlands, and erosion modelling can be useful to quantify this threat. To test its applicability for this region, the LISEM soil erosion model was applied to two small catchments, one in the Usumbara Mountains, Tanzania, and the other on the slopes of Mount Kenya. Input data for the model were collected in both catchments, as were data on runoff and erosion that were used for calibration and validation of the model. LISEM was first calibrated on catchment outlet data, and afterwards simulated spatial patterns of erosion were compared to available erosion data. The results showed that LISEM can, after calibration, give good discharge predictions for some events, but not for all. However, LISEM generally overpredicted soil loss from the catchments. Comparison with observed erosion patterns did not show overprediction, but according to the model, erosion was more widespread than was observed. There are several reasons for these discrepancies. First, it is difficult to obtain enough accurate data to run the model, such as accurate maps, rainfall data and soil and plant characteristics. Second, it is also difficult to obtain accurate data to evaluate the performance of the model, either for the catchment outlet or spatially, therefore observed erosion rates are also uncertain. Third, the model could not deal correctly with complex events, i.e. those having double rainfall peaks, and might also have difficulties with catchment characteristics such as soil type and the complexity of land use. Finally, LISEM could not deal with events in which throughflow or baseflow played a role, which was to be expected since those processes are not simulated by LISEM. Nevertheless, LISEM could be calibrated to give good discharge predictions for some events, and also gave reasonable results when compared to data obtained from erosion plots. Furthermore, only complex, distributed, storm‐based models such as LISEM can give spatial predictions for single storms. Therefore, it is concluded that if the aim is spatial prediction on an event basis, there is no alternative to complex erosion models such as LISEM, but if the aim is to predict average annual erosion, the data‐demanding, physically based LISEM erosion model may not be the most appropriate model. Copyright © 2005 John Wiley & Sons, Ltd.     

A one‐dimensional model for simulating armouring and erosion on hillslopes: 1. model development and event‐scale dynamics

This paper presents an erosion model, ARMOUR, which simulates time‐varying runoff, erosion, deposition and surface armour evolution down a hillslope either as a result of a single erosion event or as the cumulative impact of many events over periods up to decades. ARMOUR simulates sediment transport for both cohesive and non‐cohesive soil and dynamically differentiates between ‘transport‐limited’ and ‘source‐limited’ processes. A variety of feasible processes for entrainment of different size classes can be modelled and evaluated against data. The generalized likelihood of uncertainty estimation (GLUE) technique was used to calibrate and validate ARMOUR using data collected during rainfall simulator experiments at two contrasting sites: (1) non‐cohesive stony sediments at Ranger Uranium Mine, Northern Territory, Australia; and (2) cohesive silty sediments at Northparkes Gold Mine, NSW, Australia. The spatial and temporal variations of model predictions within the individual runoff events showed that some entrainment processes could not model the spikes in concentration and subsequent depletion, while the hiding model of Andrews and Parker best simulated the concentration trends for both calibrated and independent runoff events. ARMOUR also successfully captured the coarsening of the surface material, though small, over the duration of the rainfall simulator trials. This was driven by the depletion of the finest size class of the soil. For a constant discharge, ARMOUR simulated higher sediment flux at the start of the storm with the sediment flux and concentration diminishing with time. For natural rainfall a power law relationship between sediment flux and discharge was observed. The calibration exercise showed that sediment concentration and discharge alone are insufficient to calibrate all aspects of the physics, in particular the armour depth. This appears to be because the armouring during the short duration events is driven by depletion of the finest classes of the sediments (diameters less then 62·5 mm), which are not normally measured. Copyright © 2006 John Wiley & Sons, Ltd.