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.     

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.     

Evaluating post-wildfire logging-slash cover treatment to reduce hillslope erosion after salvage logging using ground measurements and remote sensing

Continuing long and extensive wildfire seasons in the Western US emphasize the need for better understanding of wildfire impacts including post-fire management scenarios. Advancements in our understanding of post-fire hillslope erosion and watershed response such as flooding, sediment yield, and debris flows have recently received considerable attention. The potential impacts of removing dead trees, called salvage logging, has been studied, however the use of remotely sensed imagery after salvage logging to evaluate spatial patterns and recovery is novel. The 2015 North Star Fire provided an opportunity to evaluate hillslope erosion reduction using two field experiments and coincidental remotely sensed imagery over 3 years. Simulated rill experiments with four flow rates were used to quantify hillslope erosion on skidder trails with and without added logging slash compared with a burned-only control. Seven replicated hillslope silt fence plots with the same treatments were also evaluated for natural rainfall events. WorldView-2 satellite imagery was used to relate ground cover and erodible bare soil between the two experiments using multi-temporal Normalized Differenced Vegetation Index (NDVI) values. Results indicate that the skid trails produced significantly more sediment (0.70 g s⁻¹) than either the slash treated skid trail (0.34 g s⁻¹) or controls (0.04 g s⁻¹) with the simulated rill experiment. Similarly, under natural rainfall conditions sediment yield from hillslope silt fence plots was significantly greater for the skid trail (3.42 Mg ha⁻¹) than either the slash treated skid trail (0.18 Mg ha⁻¹) or controls (0 Mg ha⁻¹). An NDVI value of 0.32 on all plots over all years corresponded to a ground cover of about 60% which is an established threshold for erosion reduction. Significant relationships between NDVI, ground cover, and sediment values suggest that NDVI may help managers evaluate ground cover and erosion potential remotely after disturbances such as a wildfire or salvage logging.     

Application of the WEPP model to determine sources of run‐off and sediment in a forested watershed

This study investigates critical run‐off and sediment production sources in a forested Kasilian watershed located in northern Iran. The Water Erosion Prediction Project (WEPP) watershed model was set up to simulate the run‐off and sediment yields. WEPP was calibrated and validated against measured rainfall–run‐off–sediment data. Results showed that simulated run‐off and sediment yields of the watershed were in agreement with the measured data for the calibration and validation periods. While low and medium values of run‐off and sediment yields were adequately simulated by the WEPP model, high run‐off and sediment yield values were underestimated. Performance of the model was evaluated as very good and satisfactory during the calibration and validation stages, respectively. Total soil erosion and sediment load of the study watershed during the study period were determined to be 10 108 t yr^?1 and 8735 t yr^?1, respectively. The northern areas of the watershed with dry farming were identified as the critical erosion prone zones. To prioritize the subwatersheds based on their contribution to the run‐off and sediment production at the watershed's main outlet, unit response approach (URA) was applied. In this regard, subwatersheds close to the main outlet were found to have the highest contribution to sediment yield of the whole watershed. Results indicated that depending on the objective of land and water conservation practices, particularly, for controlling sediment yield at the main outlet, critical areas for implementing the best management practices may be identified through conjunctive application of WEPP and URA. Copyright © 2014 John Wiley & Sons, Ltd.     

MATHEMATICAL MODEL OF OVERLAND FLOW AND MECHANISM OF SOIL CONSERVATION FOR FORESTED STEEP HILLSLOPE－－（II）Mechanism of Soil Conservation Over Steep Hillslope

ＭＡＴＨＥＭＡＴＩＣＡＬＭＯＤＥＬＯＦＯＶＥＲＬＡＮＤＦＬＯＷＡＮＤＭＥＣＨＡＮＩＳＭＯＦＳＯＩＬＣＯＮＳＥＲＶＡＴＩＯＮＦＯＲＦＯＲＥＳＴＥＤＳＴＥＥＰＨＩＬＬＳＬＯＰＥ－－（ＩＩ）ＭｅｃｈａｎｉｓｍｏｆＳｏｉｌＣｏｎｓｅｒｖａｔｉｏｎＯｖｅｒＳｔ...     

Assessing dominant factors affecting soil erosion using a portable rainfall simulator

Investigating the causes of soil erosion is difficult in natural conditions owing to the presence of other factors. Without simplifying the experimental conditions, studying soil behavior with its numerous parameters while considering factors such as vegetation cover, topography, and rainfall is difficult and in most conditions impossible. The application of simulation approaches is therefore necessary to simplify the prototype. In this research, the effects of physical soil factors such as texture and antecedent soil moisture, along with land slope and vegetation cover were evaluated in the Taleghan watershed, Iran, using a rainfall simulator and soil erosion plots. For this purpose, a 89 × 120 cm rainfall simulator producing 24.5 and 32 mm/h rainfall intensities of 30 min duration, as a common condition of the study area, was used at 144 locations over soil erosion plots with dimensions of 95 × 125 cm. Plots had slope classes of 12-20 and 20-30 %, different soil textures, different antecedent soil moistures, and medium to poor vegetation cover conditions. It was found that for 24.5 and 32 mm/h rainfall intensities, the sediment yield had high correlations of-0.771 and -0.796 with vegetation cover and slight correlations of 0.045 and 0.029 with land slope respectively. Regression equations for predicting the sediment yield were also developed for different conditions.     

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.     

Spatial discretization of large watersheds and its influence on the estimation of hillslope sediment yield

The combined use of water erosion models and geographic information systems has facilitated soil loss estimation at the watershed scale. Tools such as the Geo‐spatial interface for the Water Erosion Prediction Project (GeoWEPP) model provide a convenient spatially distributed soil loss estimate but require discretization to identify hillslopes and channels. In GeoWEPP, the TOpographic PArameteriZation (TOPAZ) model is used as an automated procedure to extract a watershed boundary, hillslopes and channels from a digital elevation model (DEM). Previous studies in small watersheds have shown that the size of the hillslopes and the channel distribution affect the model estimates, but in large watersheds, the effects on the soil loss estimates have yet to be tested. Therefore, the objective of this study was to evaluate the effect of discretization on the hillslope sediment yield estimates using GeoWEPP in two large watersheds (>10 km²). The watersheds were selected and discretized varying the TOPAZ parameters [critical source area (CSA) and minimum source channel length (MSCL)] in a 30‐m resolution digital elevation model. The drainage networks built with TOPAZ were compared with each other using the drainage density index. The results showed that the discretization affected hillslope sediment yield estimates and their spatial distribution more than the total runoff. The drainage density index and the hillslope sediment yield were proportional but inversely related; thus, soil loss estimates were highly affected by the spatial discretization. As a result of this analysis, a method to choose the CSA and MSCL values that generates the greatest fraction of hillslopes having profile lengths less than 200 m was developed. This slope length condition is particularly crucial when using the WEPP and GeoWEPP models, in order for them to produce realistic estimates of sheet and rill erosion. Finally, and as a result of this analysis, a more reliable method was developed for selecting the TOPAZ channel network parameters (CSA and MSCL). Copyright © 2015 John Wiley & Sons, Ltd.     

The effect of wildfire on runoff and erosion in native Eucalyptus forest

Wildfires raise concerns over the risk of accelerated erosion as a result of increased overland flow and decreased protection of the soil by litter and ground vegetation cover. We investigated these issues following the 1994 fires that burnt large areas of native Eucalyptus forest surrounding Sydney, Australia. A review of previous studies identifies the fire and rainfall conditions that are likely to lead to increased runoff and accelerated erosion. We then compare runoff and erosion between burnt and unburnt sites for 10 months after the 1994 fires. At the scale of hillslope plots, the 1994 fire increased runoff by enhancing soil hydrophobicity, and greatly increased sediment transport, mainly through the reduced ground cover, which lowered substantially the threshold for initial sediment movement. However, both runoff and sediment transport were very localized, resulting in little runoff or sediment yield after the fire at the hillslope catchment scale. We identify that after moderately intense fires, rainfall events of greater than one year recurrence interval are required to generate substantial runoff and sediment yield. Such events did not occur during the monitoring period. Past work shows that mild burns have little effect on erosion, and it is only after the most extreme fires that erosion is produced from small, frequent storms. © 1998 John Wiley & Sons, Ltd.     

Application of the WEPP model for prioritization and evaluation of best management practices in an Indian watershed

The pre‐calibrated and validated physically based watershed model, water erosion prediction project (WEPP) was used as a modelling tool for the identification of critical watersheds and evaluation of best management practices for a small hilly watershed (Karso) of India. The land use/cover of the study area was generated using IRS‐1C LISS‐III (linear imaging self scanner) satellite data. The watershed and sub‐watershed boundaries, drainage, slope and soil map of the study area were generated using ARC/INFO geographic information system (GIS). The WEPP model was finally applied to the Karso watershed which lies within Damodar Barakar catchment of India to identify the critical sub‐watersheds on the basis of their simulated average annual sediment yields. Priorities were fixed on the basis of ranks assigned to each critical sub‐watershed based on the susceptibility to erosion. The sub‐watershed having the highest sediment yield was assigned a priority number 1, the next highest value was assigned a priority number 2, and so on. Subsequently, the model was used for evaluating the effectiveness of best management practices (crop and tillage) for conservation of soil for all the sub‐watersheds. On the basis of this study, it is realized that cash crops like soyabean should be encouraged in the upland portion of the sub‐watersheds, and the existing tillage practice (country plough/mould board plough) may be replaced by a field cultivation system for conservation of soil and water in the sub‐watersheds. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

Hydrological modelling of a small watershed using generated rainfall in the soil and water assessment tool model

An adequately tested soil and water assessment tool (SWAT) model was applied to the runoff and sediment yield of a small agricultural watershed in eastern India using generated rainfall. The capability of the model for generating rainfall was evaluated for a period of 18 years (1981–1998). The watershed and subwatershed boundaries, drainage networks, slope, soil series and texture maps were generated using a geographical information system (GIS). A supervised classification method was used for land‐use/cover classification from satellite imageries. Model simulated monthly rainfall for the period of 18 years was compared with observations. Simulated monthly rainfall, runoff and sediment yield values for the monsoon season of 8 years (1991–1998) were also compared with their observed values. In general monthly average rainfall predicted by the model was in close agreement with the observed monthly average values. Also, simulated monthly average values of surface runoff and sediment yield using generated rainfall compared well with observed values during the monsoon season of the years 1991–1998. Results of this study revealed that the SWAT model can generate monthly average rainfall satisfactorily and thereby can produce monthly average values of surface runoff and sediment yield close to the observed values. Therefore, it can be concluded that the SWAT model could be used for developing a multiple year management plan for the critical erosion prone areas of a small watershed. Copyright © 2004 John Wiley & Sons, Ltd.     

Model study of the relationship between sediment yield and river basin area

The SHETRAN physically based, spatially distributed model is used to investigate the scaling relationship linking specific sediment yield to river basin area, for two contrasting topographies of upland and more homogeneous terrain and as a function of sediment source, land use and rainfall distribution. Modelling enables the effects of the controls to be examined on a systematic basis, while avoiding the difficulties associated with the use of field data (which include limited data, lack of measurements for nested basins and inability to isolate the effects of individual controls). Conventionally sediment yield is held to decrease as basin area increases, as the river network becomes more remote from the headwater sediment sources (an inverse relationship). However, recent studies have reported the opposite variation, depending on the river basin characteristics. The simulation results are consistent with these studies. If the sediment is supplied solely from hillslope erosion (no channel bank erosion) then, with uniform land use, sediment yield either decreases or is constant as area increases. The downstream decrease is accentuated if rainfall (and thence erosion) is higher in the headwaters than at lower elevations. Introducing a non‐uniform land use (e.g. forest at higher elevations, wheat at lower elevations) can reverse the trend, so that sediment yield increases downstream. If the sediment is supplied solely from bank erosion (no hillslope erosion), the sediment yield increases downstream for all conditions. The sediment yield/basin area relationship can thus be inverse or direct, depending on basin characteristics. There still remains, therefore, considerable scope for defining a universal scaling law for sediment yield. Copyright © 2006 John Wiley & Sons, Ltd.     

The changing dynamics of rill erosion on sloping farmland during the different growth stages of a maize crop

Rill erosion is a serious concern in the hilly region of China with purple soil, and maize is extensively cultivated in this region. Evaluations of the dynamic mechanisms of rill erosion in sloping farmland areas are particularly important during the maize growing season to determine whether rill erosion can occur. A new ridge tillage (RT) system was designed using local agricultural methods in China. Twelve artificial rainfall experiments were conducted in three 1 × 2 m experimental plots with a slope of 15°, which is a typical slope in the study area. The rainfall intensities were designated as 1.0, 1.5, and 2.0 mm min^?1. The rainfall experiments were performed in the field to determine the characteristics of run‐off and sediment transport related to rill erosion processes during different stages of maize growth and to analyse how hydraulic parameters and the sediment yield of the rill erosion process are related. The results showed that rill flow patterns were mainly classified as subcritical transition flow during all the growth stages of maize. The effects of hydrodynamic parameters on the sediment yield were ordered as follows: Reynolds number > stream power > Froude number > shear stress. The total sediment yield varied by stage as follows: seedling stage > jointing stage > mature stage > tasseling stage. The sediment yield and run‐off rate exhibited a linear relationship that was well described at the hillslope scale. To initiate soil loss in sloping farmland areas with purple soil during the maize growing season, the critical hydrodynamic shear stress and stream power must be at least 46.505 Pa and 1.541 N m^?1 s^?1, respectively.     

Processes of accelerated pluvial erosion on desert hillslopes modified by vehicular traffic

Accelerated pluvial erosion on hillslopes modified by off-road vehicles (ORVs) is analysed using results from 50 rainfall simulation experiments conducted in the Mojave Desert, California. Sediment yield from 1 m² hillslope plots subjected to intense, 20-minute rainfalls is typically increased 10 to 20-fold following ORV use. Salient effects of vehicle traffic, which reduce infiltration, increase runoff sediment transport efficiency, and enhance gully formation, are further studied by combining simple theoretical relations with experimental data. This analysis helps identify factors controlling erosion on natural desert hillslopes, as well as those used by ORVs. Erosion of natural or vehicle-used desert surfaces is heavily influenced by runoff hydraulics. Calculated Darcy-Weisbach friction factors decrease by an average of 13-fold following vehicular slope modification, whereas runoff Reynolds numbers increase by an average of 5 1/2-fold. The capacity of overland flow to transport sediment is related to runoff power and its degree of localization, which usually increase considerably following ORV activity; however, the ability of overland flow to move large grains (competency) is related to a combination of factors not always systematically influenced by ORV use. Kinematic runoff routing, which is used to extrapolate experimental results to longer slope lengths, leads to the suggestion that the hydraulic roughness of desert hillslopes strongly influences their erosional behaviour.     

The role of vegetation patterns in structuring runoff and sediment fluxes in drylands

The dynamics of vegetation‐driven spatial heterogeneity (VDSH) and its function in structuring runoff and sediment fluxes have received increased attention from both geomorphological and ecological perspectives, particularly in arid regions with sparse vegetation cover. This paper reviews the recent findings in this area obtained from field evidence and numerical simulation experiments, and outlines their implications for soil erosion assessment. VDSH is often observed at two scales, individual plant clumps and stands of clumps. At the patch scale, the local outcomes of vegetated patches on soil erodibility and hydraulic soil properties are well established. They involve greater water storage capacity as well as increased organic carbon and nutrient inputs. These effects operate together with an enhanced capacity for the interception of water and windborne resources, and an increased biological activity that accelerates breakdown of plant litter and nutrient turnover rates. This suite of relationships, which often involve positive feedback mechanisms, creates vegetated patches that are increasingly different from nearby bare ground areas. By this way a mosaic builds up with bare ground and vegetated patches coupled together, respectively, as sources and sinks of water, sediments and nutrients. At the stand scale within‐storm temporal variability of rainfall intensity controls reinfiltration of overland flow and its decay with slope length. At moderate rainfall intensity, this factor interacts with the spatial structure of VDSH and the mechanism of overland flow generation. Reinfiltration is greater in small‐grained VDSH and topsoil saturation excess overland flow. Available information shows that VDSH structures of sources and sinks of water and sediments evolve dynamically with hillslope fluxes and tune their spatial configurations to them. Rainfall simulation experiments in large plots show that coarsening VDSH leads to significantly greater erosion rates even under heavy rainfall intensity because of the flow concentration and its velocity increase. Copyright © 2005 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…     

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.     

Sediment transport,redistribution and storage on logged forest hillslopes in south‐eastern Australia

Post‐logging changes in catchment sediment yield have traditionally been attributed to increases in hillslope erosion and delivery rates as a result of forest harvesting activities. Linking hillslope erosion to catchment yield in forestry environments remains difficult, however, primarily because of the scarcity of data on the nature of hillslope sediment storage and delivery processes. A large rainfall simulator (350 m²) was used to apply rainstorms to a logged hillslope containing a snig track (skid trail) and a general logging or harvesting area (GHA) on 10 forest compartments in south‐eastern Australia. The experiments confirmed that the compacted, disturbed surfaces, such as roads and tracks, are the dominant sources of sediment in forestry areas. Sediment transport rates were limited by available sediment supply on both the snig track and the GHA, introducing important implications for the modelling of these surfaces using sediment transport capacity theories. Sediment delivery from the snig track to the adjacent GHA, via a cross‐bank (drainage diversion), was strongly influenced by the percentage fine fraction in the eroded sediment. Preferential deposition of coarse aggregates was measured at erosion control structures and on the adjacent GHA. Over 50% of fine‐grained material were deposited on the hillslope over a relatively short, flow path length of <5 m, highlighting the effectiveness of runoff diversion as a practice in reducing sediment flux. The transfer of water and sediment from disturbed to less disturbed parts of the landscape, and the associated potential for sediment storage, needs to be considered as part of any catchment impact assessment. Copyright © 1999 John Wiley & Sons, Ltd.     

Soil and water losses from new citrus orchards growing on sloped soils in the western Mediterranean basin

Ten representative research sites were selected in eastern Spain to assess soil erosion rates and processes in new citrus orchards on sloping soils. The experimental plots were located at representatives sites on limestone, in areas with 498 to 715 mm year^?1 mean annual rainfall, north‐facing slopes, herbicide treated, and new (less than 3 years old) plantations. Ten rainfall simulation experiments (1 h at 55 mm h^?1 on 0·25 m² plots) were carried out at each of the 10 selected study sites to determine the interill soil erosion and runoff rates. The 100 rainfall simulation tests (10 × 10 m) showed that ponding and runoff occurred in all the plots, and quickly: 121 and 195 s, respectively, following rainfall initiation. Runoff discharge was one third of the rainfall, and sediment concentration reached 10·4 g L^?1. The soil erosion rates were 2·4 Mg ha^?1 h^?1 under 5‐year return period rainfall thunderstorms. These are among the highest soil erosion rates measured in the western Mediterranean basin, similar to badland, mine spoil and road embankment land surfaces. The positive relationship between runoff discharge and sediment concentration (r² = 0·83) shows that the sediment availability is very high. Soil erosion rates on new citrus orchards growing on sloped soils are neither tolerable nor sustainable. Copyright © 2009 John Wiley & Sons, Ltd.