Size characteristics of sediment in interrill overland flow on a semiarid hillslope,Southern Arizona

This study examines the size characteristics of sediment removed from a semiarid hillslope by interrill overland flow. Rainfall simulation experiments were conducted on a runoff plot 18 m wide and 35 m long established on a piedmont hillslope in southern Arizona. The top of the plot coincided with the hillslope divide, and its outlet was located within a shallow rill. Samples of runoff were obtained from two cross-sections located in the interrill portion of the plot upslope of the rill and from a calibrated flume through which was directed interrill overland flow reaching the bottom of the plot. Analyses of sediment contained in these samples showed that sediment in interrill flow is finer than the matrix soil. The fineness of the interrill sediment compared to the matrix soil appears to be due to the inability of interrill overland flow to transport the coarser fraction of the sediment supplied to it by raindrop detachment. This finding implies that the rate of soil erosion in interrill areas is not. as is commonly supposed, limited by the rate at which raindrops can detach sediment but by the rate at which they detach sediment of a size that the overland flow is competent to transport. The relative fineness of sediment eroded from this hillslope is consistent with other evidence for the recent evolution of shrub-covered hillslopes in southern Arizona.     

Effects of rainfall patterns on runoff and rainfall-induced erosion

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

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

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

Impact of subsurface rock fragments on runoff and interrill soil loss from cultivated soils

Field and laboratory studies have indicated that rock fragments in the topsoil may have a large impact on soil properties, soil quality, hydraulic, hydrological and erosion processes. In most studies, the rock fragments investigated still remain visible at the soil surface and only properties of these visible rock fragments are used for predicting runoff and soil loss. However, there are indications that rock fragments completely incorporated in the topsoil could also significantly influence the percolation and water distribution in stony soils and therefore, also infiltration, runoff and soil loss rates. Therefore, in this study interrill laboratory experiments with simulated rainfall for 60 min were conducted to assess the influence of subsurface rock fragments incorporated in a disturbed silt loam soil at different depths below the soil surface (i.e. 0.001, 0.01, 0.05 and 0.10 m), on infiltration, surface runoff and interrill erosion processes for small and large rock fragment sizes (i.e. mean diameter 0.04 and 0.20 m, respectively). Although only small differences in infiltration rate and runoff volume are observed between the soil without rock fragments (control) and the one with subsurface rock fragments, considerable differences in total interrill soil loss are observed between the control treatment and both contrasting rock fragments sizes. This is explained by a rapid increase in soil moisture in the areas above the rock fragments and therefore a decrease in topsoil cohesion compared with the control soil profile. The observed differences in runoff volume and interrill soil loss between the control plots and those with subsurface rock fragments is largest after a cumulative rainfall (P_cum) of 11 mm and progressively decreases with increasing P_cum. The results highlight the impacts and complexity of subsurface rock fragments on the production of runoff volume and soil loss and requires their inclusion in process‐based runoff and erosion models. Copyright © 2011 John Wiley & Sons, Ltd.     

Impact of rainfall pattern on interrill erosion process

The impact of rainfall pattern on the interrill erosion process is not fully understood despite its importance. Systematic rainfall simulation experiments involving various rainfall intensities, stages, intensity sequences, and surface cover conditions were conducted in this study to investigate their effects on the interrill erosion process. Five rainfall patterns designed with the same total kinetic energy/precipitation (increasing, decreasing, rising–falling, falling–rising and constant patterns) were randomly delivered to a pre‐wet clay loam soil surface at a 10° slope gradient. Significant differences in soil losses were observed among the different rainfall patterns and stages, but there was no obvious difference in runoff. Kinetic energy flux (KE_r) was a governing factor for interrill erosion, and constant rainfall pattern (CST) produced nine times greater soil loss than runs with no KE_r. Varied‐intensity patterns had a profound effect on raindrop‐induced sediment transport processes; path analysis results indicated that said effect was complex, interactive and intensity‐dependent. Low hydraulic parameter thresholds further indicated that KE_r was the dominant factor in detaching soil particles, while overland flow mainly contributed to transporting the pre‐detached particles. This study not only sheds light on the mechanism of interrill sediment transport capacity and detachability, but also may provide a useful database for developing event‐based interrill erosion prediction models. Copyright © 2017 John Wiley & Sons, Ltd.     

Sensitivity analysis of EUROSEM using Monte Carlo simulation II: the effect of rills and rock fragments

A sensitivity analysis of the surface and catchment characteristics in the European soil erosion model (EUROSEM) was carried out with special emphasis on rills and rock fragment cover. The analysis focused on the use of Monte Carlo simulation but was supplemented by a simple sensitivity analysis where input variables were increased and decreased by 10%. The study showed that rock fragments have a significant effect upon the static output parameters of total runoff, peak flow rate, total soil loss and peak sediment discharge, but with a high coefficient of variation. The same applied to the average hydrographs and sedigraphs although the peak of the graphs was associated with a low coefficient of variation. On average, however, the model was able to simulate the effect of rock fragment cover quite well. The sensitivity analysis through the Monte Carlo simulation showed that the model is particularly sensitive to changes in parameters describing rills and the length of the plane when no rock fragments are simulated but that the model also is sensitive to changes in the fraction of non‐erodible material and interrill slope when rock fragments were embedded in the topsoil. For rock fragments resting on the surface, changes in parameter values did not affect model output significantly. The simple sensitivity analysis supported the findings from the Monte Carlo simulation and illustrates the importance when choosing input parameters to describe both rills and rock fragment cover when modelling with EUROSEM. Copyright © 2000 John Wiley & Sons, Ltd.     

An un-mixing model to study watershed erosion processes

An un-mixing model is formulated within a Bayesian Markov Chain Monte Carlo framework for use within land-use fingerprinting to study watershed erosion processes. The model has two new components: (1) An equation and erosion process parameter are used to weight tracer signatures from each erosion process within a land-use. (2) An extra tracer distribution and episodic erosion parameter are used to represent soil eroded throughout the sampling duration and thus include the episodic nature of erosion. To test specification of these new parameters, the un-mixing model is applied in the 15 km² Jerome Creek Watershed in the Palouse Region of Northwestern Idaho. Erosion processes include surface erosion upon mountain slopes due to logging in the forest land-use and rill/interrill erosion on cultivated slopes and headcut erosion in riparian floodplains of the agricultural land-use (winter wheat/peas rotation and hay pasture). Episodic erosion occurs for the event where the model is applied. A sensitivity analysis shows that the smallest Bayesian credible set results when the new parameters are specified using hydrologic data and process-based models. The un-mixing model predicts that 90% of the eroded-soil originated from the agricultural land-use and 10% originated from the forest land-use. A comparative study is performed that estimates 90.5% and 9.5% of eroded-soil originated from the agricultural and forest land-uses. Successful performance of the un-mixing model highlights future application as a standalone probabilistic tool to monitor watershed erosion processes that exhibit non-equilibrium conditions and provide calibration data for process-based watershed models.     

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.     

Plot‐scale measurement of soil erosion at the experimental area of Sparacia (southern Italy)

Obtaining good quality soil loss data from plots requires knowledge of the factors that affect natural and measurement data variability and of the erosion processes that occur on plots of different sizes. Data variability was investigated in southern Italy by collecting runoff and soil loss from four universal soil‐loss equation (USLE) plots of 176 m², 20 ‘large’ microplots (0·16 m²) and 40 ‘small’ microplots (0·04 m²). For the four most erosive events (event erosivity index, R_e ≥ 139 MJ mm ha^?1 h^?1), mean soil loss from the USLE plots was significantly correlated with R_e. Variability of soil loss measurements from microplots was five to ten times greater than that of runoff measurements. Doubling the linear size of the microplots reduced mean runoff and soil loss measurements by a factor of 2·6–2·8 and increased data variability. Using sieved soil instead of natural soil increased runoff and soil loss by a factor of 1·3–1·5. Interrill erosion was a minor part (0·1–7·1%) of rill plus interrill erosion. The developed analysis showed that the USLE scheme was usable to predict mean soil loss at plot scale in Mediterranean areas. A microplot of 0·04 m² could be used in practice to obtain field measurements of interrill soil erodibility in areas having steep slopes. Copyright © 2003 John Wiley & Sons, Ltd.     

MODELING OF INTERRILL EROSION IN THE LOESS PLATEAU OF CHINA

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

Selecting best mapping strategies for storm runoff modeling in a mountainous semi‐arid area

Accurate runoff and soil erosion modeling is constrained by data availability, particularly for physically based models such as OpenLISEM that are data demanding, as the processes are calculated on a cell‐by‐cell basis. The first decision when using such models is to select mapping units that best reflect the spatial variability of the soil and hydraulic properties in the catchment. In environments with limited data, available maps are usually generic, with large units that may lump together the values of the soil properties, affecting the spatial patterns of the predictions and output values in the outlet. Conversely, the output results may be equally acceptable, following the principle of equifinality. To studyhow the mapping method selected affects the model outputs, four types of input maps with different degrees of complexity were created: average values allocated to general soil map units (ASG1), average values allocated to detailed map units (ASG2), values interpolated by ordinary kriging (OK) and interpolated by kriging with external drift (KED). The study area was Ribeira Seca, a 90 km² catchment located in Santiago Island, Cape Verde (West Africa), a semi‐arid country subject to scarce but extreme rainfall during the short tropical summer monsoon. To evaluate the influence of rainfall on runoff and erosion, two storm events with different intensity and duration were considered. OK and KED inputs produced similar results, with the latter being closer to the observed hydrographs. The highest soil losses were obtained with KED (43 ton ha^{? 1} for the strongest event). To improve the results of soil loss predictions, higher accurate spatial information on the processes is needed; however, spatial information of input soil properties alone is not enough in complex landscapes. The results demonstrate the importance of selecting the appropriate mapping strategy to obtain reliable runoff and erosion estimates. Copyright © 2013 John Wiley & Sons, Ltd.     

Erodibility of soil and organic matter: independence of organic matter resistance to interrill erosion

The enrichment of organic matter in interrill sediment is well documented; however, the respective roles of soil organic matter (SOM) and interrill erosion processes for the enrichment are unclear. In this study, organic matter content of sediment generated on two silts with almost identical textures, but different organic matter contents and aggregations, was tested. Artificial rainfall was applied to the soils in wet, dry and crusted initial conditions to determine the effects of soil moisture and rainfall and drying history on organic matter enrichment in interrill sediment. While erosional response of the soils varied significantly, organic matter enrichment of sediment was not sensitive to initial soil conditions. However, enrichment was higher on the silt with a lower organic matter content and lower interrill erodibility. The results show that enrichment of organic matter in interrill sediment is not directly related to either SOM content or soil interrill erodibility, but is dominated by interrill erosion processes. As a consequence of the complex interaction between soil, organic matter and interrill erosion processes, erodibility of organic matter should be treated as a separate variable in erosion models. Further research on aggregate breakdown, in particular the content and fate of the organic matter in the soil fragments, is required. Copyright © 2007 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 GIS-BASED DISTRIBUTED SOIL EROSION MODEL： A CASE STUDY OF TYPICAL WATERSHED, SICHUAN BASIN

Based on the measuring data and Digital Elevation Data (DEM) in a typical watershed--Hemingguan Watershed, Nanbu County, Sichuan Province of China, a GIS-based distributed soil erosion model was developed particularly for the purple soil type. It takes 20 m × 20 m grid as calculating unit and operates at 10-minute time interval. The required input data to the model include DEM, soil, land use, and time-series of precipitation and evaporation loss. The model enables one to estimate runoff, erosion and sediment yield for each grid cell and route the flow along its flow path to the watershed outlet. Furthermore, the model is capable of calculating the total runoff; erosion and sediment yield for the entire watershed by recursion algorithm. The validation of the model demonstrated that it could quantitatively simulate the spatial distribution of hydrological variables in a watershed, such as runoff, vegetation entrapment, soil erosion, the degree of soil and water loss. Moreover, it can evaluate the effect of land use change on the runoff generation and soil erosion with an accuracy of 80% and 75% respectively. The application of this model to a neighboring watershed with similar conditions indicates that this distributed model could be extended to other similar regions in China.     

Predicting the impact of linear landscape elements on surface runoff,soil erosion,and sedimentation in the Wahnbach catchment,Germany

Model predictions concerning the endangerment of on‐site and off‐site damages due to runoff, soil erosion and sedimentation under alternative design and operation policies are of particular importance in recent catchment planning and management. By using the raster‐based model approach, linear landscape elements, such as streets and roads, and their impacts on flow paths are often neglected. Therefore, the aim of this study was to analyse the effects of linear landscape elements on patterns of soil erosion, sediment transport and sedimentation. To accomplish this, roads are considered while determining flow paths. Simulations in the well‐investigated catchment of the Wahnbach River (54 km²) in a low mountain range in Germany were carried out using a combination of different models for hydrology and soil erosion. Although the study focuses on the catchment scale of the Wahnbach River, detailed investigations concerning the sub‐catchment scale (21 ha) were also conducted. The simulation results show that these spatial structures mainly affect the pattern of soil erosion and sedimentation. On the sub‐catchment scale, improved identification of active zones for sediment dynamic becomes possible. On the catchment scale, the predicted runoff is about 20% higher, and sediment outputs were four times larger than predicted when roads were considered. Soil erosion increases by 37% whereas sedimentation is reduced by 29%. The model improvement could not be evaluated on the catchment scale because of the high variability and heterogeneity of land use and soils, but road impacts could be explained by simulations on the sub‐catchment scale. It can be concluded that runoff concentration due to rerouted flow paths leads to lower non‐concentrated and higher concentric‐linear surface runoff. Thus, infiltration losses decline and surface runoff and soil erosion increase because sedimentation is reduced. Further, runoff concentration can cause soil erosion hot spots. In the model concept used in this study, buffering of runoff and sediments on the upslope side of roads and in local depressions adjacent to roads cannot be simulated. Flow paths will only be rerouted because of road impacts, but the temporal ponding of water is not simulated. Therefore, the drastic increase of predicted sediment output due to road impact does not seem to be reliable. However, results indicate that the consideration of roads when determining flow paths enabled more detailed simulations of surface runoff, soil erosion and sedimentation. Thus, progress in model‐based decision‐making support for river catchment planning and management can be achieved. Copyright © 2011 John Wiley & Sons, Ltd.     

RUNOFF GENERATION CHARACTERISTICS IN TYPICAL EROSION REGIONS ON THE LOESS PLATEAU

1 mTRODUCT1ONThe Loess Plateau in the central and westem China with an area of 430,(X)0 lQn2 is well known for itsancient cultUre and serious soil and water loss. About 287,(X)0 km2 of area in this region has an annualerosion be greatC than 1(XX) tlkIn2. The eroded soil totally amounts to more than 2.2 X l0' tons, amongwhich about l.6X 10' tons, on average, is delivered into the Yellow kiver mdeng the river one of thehighest sediment-laden rivers in the word. The serious soil erosion …     

Evaluation of an interrill soil erosion model using laboratory catchment data

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

Fire effects on rangeland hydrology and erosion in a steep sagebrush‐dominated landscape

Post‐fire runoff and erosion from wildlands has been well researched, but few studies have researched the degree of control exerted by fire on rangeland hydrology and erosion processes. Furthermore, the spatial continuity and temporal persistence of wildfire impacts on rangeland hydrology and erosion are not well understood. Small‐plot rainfall and concentrated flow simulations were applied to unburned and severely burned hillslopes to determine the spatial continuity and persistence of fire‐induced impacts on runoff and erosion by interrill and rill processes on steep sagebrush‐dominated sites. Runoff and erosion were measured immediately following and each of 3 years post‐wildfire. Spatial and temporal variability in post‐fire hydrologic and erosional responses were compared with runoff and erosion measured under unburned conditions. Results from interrill simulations indicate fire‐induced impacts were predominantly on coppice microsites and that fire influenced interrill sediment yield more than runoff. Interrill runoff was nearly unchanged by burning, but 3‐year cumulative interrill sediment yield on burned hillslopes (50 g m^?2) was twice that of unburned hillslopes (25 g m^?2). The greatest impact of fire was on the dynamics of runoff once overland flow began. Reduced ground cover on burned hillslopes allowed overland flow to concentrate into rills. The 3‐year cumulative runoff from concentrated flow simulations on burned hillslopes (298 l) was nearly 20 times that measured on unburned hillslopes (16 l). The 3‐year cumulative sediment yield from concentrated flow on burned and unburned hillslopes was 20 400 g m^?2 and 6 g m^?2 respectively. Fire effects on runoff generation and sediment were greatly reduced, but remained, 3 years post‐fire. The results indicate that the impacts of fire on runoff and erosion from severely burned steep sagebrush landscapes vary significantly by microsite and process, exhibiting seasonal fluctuation in degree, and that fire‐induced increases in runoff and erosion may require more than 3 years to return to background levels. Published in 2008 by John Wiley & Sons, Ltd.     

Interrill erosion in the sloping lands of northern Laos subjected to shifting cultivation

In this study our main objective was to quantify water interrill erosion in the sloping lands of Southeast Asia, one of the most bio‐geochemically active regions of the world. Investigations were performed on a typical hillslope of Northern Laos subjected to slash and burn agriculture practiced as shifting cultivation. Situations with different periods of the shifting cultivation cycle (secondary forest, upland rice cultivation following a four‐year fallow period and three‐year continuous upland rice cultivation) and soil orders (Ultisols, Alfisols, Inceptisols) were selected. One metre square micro‐plots were installed to quantify the soil material removed by either detachment of entire soil aggregate or aggregate destruction, and the detached material transported by thin sheet flow, the main mechanisms of interrill erosion. In addition, laboratory tests were carried out to quantify the aggregate destruction in the process of water erosion by slaking, dispersion and mechanical breakdown. The average runoff coefficient (R) evaluated throughout the 2002 rainy season was 30·1 per cent and the interrill erosion was 1413 g m^?2 yr^?1 for sediments and 68 g C m^?2 yr^?1 for soil organic carbon, which was relatively high. Among the mechanisms of interrill water erosion, aggregate destruction was low and mostly caused by mechanical breakdown due to raindrops, thus leading to the conclusion that detachment and further transport by the shallow runoff of macro‐aggregates predominates. R ranged from 23·1 to 35·8 per cent. It decreased with the proportion of mosses on the soil surface and soil surface coverage, and increased with increasing proportion of structural crust, thus confirming previous results. Water erosion varied from 621 to 2433 g m^?2 yr^?1 for sediments and from 31 to 146 g C m^?2 yr^?1 for soil organic carbon, and significantly increased with increasing clay content of the surface horizon, probably due to the formation of easily detachable and transportable sand‐size aggregates, and proportion of macro‐aggregates not embedded in the soil matrix and prone to transport. In addition, water erosion decreased with increasing proportion of structural crusts, probably due to their higher hardness, and when cultivation follows a fallow period rather than after a long period of cultivation due to the greater occurrence of algae on the soil surface, which affords physical protection and greater aggregate stability through binding and gluing. This study based on simultaneous field and laboratory investigations allowed successful identification and quantification of the main erosion mechanisms and controlling factors of interrill erosion, which will give arguments to further set up optimal strategies for sustainable use of the sloping lands of Southeast Asia. Copyright © 2006 John Wiley & Sons, Ltd.     

EVALUATION OF EROSION PRODUCTIVITY IMPACT CALCULATOR （EPIC） MODEL FOR MIDDLE MOUNTAIN REGION OF NEPAL

This study verifies the applicability of EPIC model for an erosion plot (61 .2 m~2) and an uplandterraced watershed (72 ha) using a total of 94 rainfall events over a study period of two years. Inorder to analyze the effect of storm size on runoff and soil loss processes, rainfall events aredivided into three groups: small (<25mm), moderate (25--50mm) and large (>50mm). Resultsindicate that the model could predict reasonably well the runoff and soil loss from the erosion plotand the watershed for the moderate and large rainfall events. However, the runoff and soil lossprediction for the small rainfall events is found to be poor. On annual basis, both surface runoff andsoil loss predictions match well the observations. In ligh of the importance of the moderate andlarge rainfall events in producing most of the annual runoff and soil loss in the study area, the EPICmodel is applied to assess the impacts of erosion on agricultural productivity and to evaluatemanagement practices to protect watersheds in the