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.     

Connectivity of post-fire runoff and sediment from nested hillslopes and watersheds

Wildfire increases the potential connectivity of runoff and sediment throughout watersheds due to greater bare soil, runoff and erosion as compared to pre-fire conditions. This research examines the connectivity of post-fire runoff and sediment from hillslopes (< 1.5 ha; n = 31) and catchments (< 1000 ha; n = 10) within two watersheds (< 1500 ha) burned by the 2012 High Park Fire in northcentral Colorado, USA. Our objectives were to: (1) identify sources and quantify magnitudes of post-fire runoff and erosion at nested hillslopes and watersheds for two rain storms with varied duration, intensity and antecedent precipitation; and (2) assess the factors affecting the magnitude and connectivity of runoff and sediment across spatial scales for these two rain storms. The two summer storms that are the focus of this research occurred during the third summer after burning. The first storm had low intensity rainfall over 11 hours (return interval <1–2 years), whereas the second event had high intensity rainfall over 1 hour (return interval <1–10 years). The lower intensity storm was preceded by high antecedent rainfall and led to low hillslope sediment yields and channel incision at most locations, whereas the high intensity storm led to infiltration-excess overland flow, high sediment yields, in-stream sediment deposition and channel substrate fining. For both storms, hillslope-to-stream sediment delivery ratios and area-normalised cross-sectional channel change increased with the percent of catchment that burned at high severity. For the high intensity storm, hillslope-to-stream sediment delivery ratios decreased with unconfined channel length (%). The findings quantify post-fire connectivity and sediment delivery from hillslopes and streams, and highlight how different types of storms can cause varying magnitues and spatial patterns of sediment transport and deposition from hillslopes through stream channel networks.     

Identification of soil erosion and fluvial sediment problems

A computer model has been used to estimate soil loss and sediment yield from irregular field-size units of small watersheds. Input to the model includes spring data (i.e. relating to February through May) for the independent variables of the Universal Soil Loss Equation, and for factors such as surface roughness, an index of overland runoff, and proximity to the stream. Output from the model includes maps of seasonal estimates of potential soil losses, field sediment delivery ratios, and expected sediment yields. On the basis of selected erosion and sediment yield tolerances, the output information has been analysed to identify watershed areas which (1) exhibit both erosion and sediment yield problems; (2) exhibit only erosion problems; (3) exhibit only sediment yield problems; and (4) exhibit neither erosion nor sediment yield problems. The percentage of the watershed area in each category and the percentage of the watershed soil loss and sediment loads contributed by each category are also identified. Application of the procedure for planning remedial control programs for five watersheds is discussed.     

An integral-balance nonlinear model to simulate changes in soil moisture,groundwater and surface runoff dynamics at the hillslope scale

We present a system of ordinary differential equations (ODEs) capable of reproducing simultaneously the aggregated behavior of changes in water storage in the hillslope surface, the unsaturated and the saturated soil layers and the channel that drains the hillslope. The system of equations can be viewed as a two-state integral-balance model for soil moisture and groundwater dynamics. Development of the model was motivated by the need for landscape representation through hillslopes and channels organized following stream drainage network topology. Such a representation, with the basic discretization unit of a hillslope, allows ODEs-based simulation of the water transport in a basin. This, in turn, admits the use of highly efficient numerical solvers that enable space–time scaling studies. The goal of this paper is to investigate whether a nonlinear ODE system can effectively replicate observations of water storage in the unsaturated and saturated layers of the soil. Our first finding is that a previously proposed ODE hillslope model, based on readily available data, is capable of reproducing streamflow fluctuations but fails to reproduce the interactions between the surface and subsurface components at the hillslope scale. However, the more complex ODE model that we present in this paper achieves this goal. In our model, fluxes in the soil are described using a Taylor expansion of the underlying storage flux relationship. We tested the model using data collected in the Shale Hills watershed, a 7.9-ha forested site in central Pennsylvania, during an artificial drainage experiment in August 1974 where soil moisture in the unsaturated zone, groundwater dynamics and surface runoff were monitored. The ODE model can be used as an alternative to spatially explicit hillslope models, based on systems of partial differential equations, which require more computational power to resolve fluxes at the hillslope scale. Therefore, it is appropriate to be coupled to runoff routing models to investigate the effect of runoff and its uncertainty propagation across scales. However, this improved performance comes at the expense of introducing two additional parameters that have no obvious physical interpretation. We discuss the implications of this for hydrologic studies across scales.     

WATERSHED SEDIMENT YIELD AND RANGELAND HEALTH

lINTRoDUCTIoNDifferencesintheprevailinglanduseandmanagementofaridandsemiaridareasaredeterminedinlargepartbyclimate.AridareasgenerallyreceivetoolittleprecipitationtosupportdrylandagricultureordomesticlivestockgrazingalthoughtheyaregrazedbywildIife,andattimes,bydomesticlivestock.Incontrast,insemiaridareasadequatemoistureisusuallyavaiIableatsometimeduringtheyeartoproduceforageforlivestockandwildlife,andtherearesomeyearswhendrylandcropproductionissuccessful.However,bothclimatesarecharacterize…     

Plot and small‐watershed scale runoff from two reclaimed surface‐mined watersheds in Alberta

Infiltration excess overland flow has been identified as the dominant flow pathway in recently reclaimed surface mined watersheds as a result of compaction and sorting during the reclamation procedure. Therefore, there could be a fairly direct relationship between runoff generated from the hillslopes to that measured at the watershed outlet. A 3‐year study was initiated in 1993 to determine how well surface runoff at a watershed scale could be predicted from 1‐m² runoff frames placed on hillslopes in two reclaimed surface‐mined watersheds in central Alberta. Runoff from the hillslope frames suggests outlet discharge should be high from the 3\4‐ha Sandy Subsoil Watershed and much less for the 9\8‐ha West Watershed, but the opposite occurred. Most of the hillslope runoff from the Sandy Subsoil Watershed infiltrated once it reached the channel and depression storage played an insignificant role in determining runoff. In contrast, most of the runoff from the West Watershed originated from rain falling directly on the saturated channel (depression storage) or near‐channel saturated areas, rather than the hillslopes. Neither watershed runoff magnitude nor timing could be predicted from the same parameters for hillslope runoff frames for either reclaimed watershed. Copyright © 2000 John Wiley & Sons, Ltd.     

Spatial and temporal analysis of hillslope–channel coupling and implications for the longitudinal profile in a dryland basin

The long‐term evolution of channel longitudinal profiles within drainage basins is partly determined by the relative balance of hillslope sediment supply to channels and the evacuation of channel sediment. However, the lack of theoretical understanding of the physical processes of hillslope–channel coupling makes it challenging to determine whether hillslope sediment supply or channel sediment evacuation dominates over different timescales and how this balance affects bed elevation locally along the longitudinal profile. In this paper, we develop a framework for inferring the relative dominance of hillslope sediment supply to the channel versus channel sediment evacuation, over a range of temporal and spatial scales. The framework combines distinct local flow distributions on hillslopes and in the channel with surface grain‐size distributions. We use these to compute local hydraulic stresses at various hillslope‐channel coupling locations within the Walnut Gulch Experimental Watershed (WGEW) in southeast Arizona, USA. These stresses are then assessed as a local net balance of geomorphic work between hillslopes and channel for a range of flow conditions generalizing decadal historical records. Our analysis reveals that, although the magnitude of hydraulic stress in the channel is consistently higher than that on hillslopes, the product of stress magnitude and frequency results in a close balance between hillslope supply and channel evacuation for the most frequent flows. Only at less frequent, high‐magnitude flows do channel hydraulic stresses exceed those on hillslopes, and channel evacuation dominates the net balance. This result suggests that WGEW exists mostly (~50% of the time) in an equilibrium condition of sediment balance between hillslopes and channels, which helps to explain the observed straight longitudinal profile. We illustrate how this balance can be upset by climate changes that differentially affect relative flow regimes on slopes and in channels. Such changes can push the long profile into a convex or concave condition. © 2018 The Authors. Earth Surface Processes and Landforms published by 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.     

A numerical investigation of bedrock groundwater recharge and exfiltration on soil mantled hillslopes

Here we use Richards Equation models of variably saturated soil and bedrock groundwater flow to investigate first-order patterns of the coupling between soil and bedrock flow systems. We utilize a Monte Carlo sensitivity analysis to identify important hillslope parameters controlling bedrock recharge and then model the transient response of bedrock and soil flow to seasonal precipitation. Our results suggest that hillslopes can be divided into three conceptual zones of groundwater interaction, (a) the zone of lateral unsaturated soil moisture accumulation (upper portion of hillslope), (b) the zone of soil saturation and bedrock recharge (middle of hillslope) and (c) the zone of saturated-soil lateral flow and bedrock groundwater exfiltration (bottom of hillslope). Zones of groundwater interaction expand upslope during periods of precipitation and drain downslope during dry periods. The amount of water partitioned to the bedrock groundwater system a can be predicted by the ratio of bedrock to soil saturated hydraulic conductivity across a variety of hillslope configurations. Our modelled processes are qualitatively consistent with observations of shallow subsurface saturation and groundwater fluctuation on hillslopes studied in our two experimental watersheds and support a conceptual model of tightly coupled shallow and deep subsurface circulation where groundwater recharge and discharge continuously stores and releases water from longer residence time storage.     

Physically-based distributed soil erosion and sediment yield model (DREAM) for simulating individual storm events

Abstract

A relatively simple process-oriented, physically-based distributed (PBD) hydrological model, the distributed runoff and erosion assessment model (DREAM), is described, and a validation study conducted in the semi-forested watershed of Pathri Rao, in the Garhwal Himalayas, India, is reported. DREAM takes account of watershed heterogeneity as reflected by land use, soil type, topography and rainfall, measured in the field or estimated through remote sensing, and generates estimates of runoff and sediment yield in spatial and temporal domains. The model is based on simultaneous solution of flow dynamics, based on kinematic wave theory, followed by solution of soil erosion dynamics. As the storm rainfall proceeds, the process of overland flow generation is dependent on the interception storage and infiltration rates. The components of the soil erosion model have been modified to provide better prediction of sediment flow rates and sediment yields. The validation study conducted to test the performance of the model in simulating soil erosion and sediment yield during different storm events monitored in the study watershed showed that the model outputs are satisfactory. Details of a sensitivity analysis, model calibration and the statistical evaluation of the results obtained are also presented and discussed. It is noteworthy that the distributed nature of the model combined with the use of geographical information system (GIS) techniques permits the computation and representation of the spatial distribution of sediment yield for simulated storm events, and a map of the spatial distribution of sediment yield for a simulated storm event is presented to highlight this capability.

Citation Ramsankaran, R., Kothyari, U.C., Ghosh, S.K., Malcherek, A., and Murugesan, K., 2013. Physically-based distributed soil erosion and sediment yield model (DREAM) for simulating individual storm events. Hydrological Sciences Journal, 58 (4), 872–891.     

Evaluation of the WEPP model and digital elevation grid size,for simulation of streamflow and sediment yield in a heterogeneous catchment

The paper presents the result of an application of the GeoWEPP model in a heterogeneous semi‐agricultural catchment located in the northern Italian Apennines mountain range. The objectives were: (a) to evaluate the GeoWEPP model in a heterogeneous catchment in a Mediterranean climate and (b) to examine the effect of digital elevation model grid size on hydrological and sediment yield simulations. The catchment is characterized by large heterogeneity in geology, soil type, vegetation cover and topography. In addition, 10% of its area is occupied by calanchi (badlands), characterized by steep, bare soil and accentuated erosion. Experimental streamflow data were compared with those simulated by GeoWEPP for a period of eight years and the results were evaluated by means of statistical indices, with the analysis of the flow duration curve. Simulated sediment yields were compared with experimental data for one year. The streamflow cumulative annual results were satisfactory with NSE oscillating between 0.40 and 0.83 and RMSE between 1.1 and 2.9 mm. Also, the performance of the model with daily streamflow data was positive (NSE = 0.68 and RMSE = 1.9 mm). The flow duration curve indicated that GeoWEPP could represent the experimental streamflow for fluxes over 1 mm d^?1. The model performance for simulation of sediment yield was satisfactory with both digital elevation models of different grid sizes (NSE = 0.84 and 0.87). Indeed, the sensitivity analysis tests of the model showed that there was no statistically significant improvement in the accuracy of the digital elevation model between 10 and 2 m resolution. These results were confirmed for both streamflow as well as sediment yield. Additional sensitivity analysis of other model parameters performed on the entire catchment and badlands hillslopes showed that bedrock hydraulic conductivity primarily affected the model in both settings. Copyright © 2014 John Wiley & Sons, Ltd.     

Rainfall–runoff responses on Arctic hillslopes underlain by continuous permafrost,North Slope,Alaska, USA

The Arctic hydrologic cycle is intensifying, as evidenced by increased rates of precipitation, evapotranspiration, and riverine discharge. However, the controls on water fluxes from terrestrial to aquatic systems in upland Arctic landscapes are poorly understood. Upland landscapes account for one third of the Arctic land surface and are often drained by zero‐order geomorphic flowpath features called water tracks. Previous work in the region attributed rapid runoff response at larger stream orders to water tracks, but models suggest water tracks are hydrologically disconnected from the surrounding hillslope. To better understand the role of water tracks in upland landscapes, we investigated the surface and subsurface hydrologic responses of 6 water tracks and their hillslope watersheds to natural patterns of rainfall, soil thaw, and drainage. Between storms, both water track discharge and the water table in the hillslope watersheds exhibited diel fluctuations that, when lagged by 5 hr, were temporally correlated with peak evapotranspiration rate. Water track soils remained saturated for more of the summer season than soils in their surrounding hillslope watersheds. When rainfall occurred, the subsurface response was nearly instantaneous, but the water tracks took significantly longer than the hillslopes to respond to rainfall, and longer than the responses previously observed in nearby larger order Arctic streams. There was also evidence for antecedent soil water storage conditions controlling the magnitude of runoff response. Based on these observations, we used a broken stick model to test the hypothesis that runoff production in response to individual storms was primarily controlled by rainfall amount and antecedent water storage conditions near the water track outlet. We found that the relative importance of the two factors varied by site, and that water tracks with similar watershed geometries and at similar landscape positions had similar rainfall–runoff model relationships. Thus, the response of terrestrial water fluxes in the upland Arctic to climate change depends on the non‐linear interactions between rainfall patterns and subsurface water storage capacity on hillslopes. Predicting these interactions across the landscape remains an important challenge.     

How long is a hillslope?

Hillslope length is a fundamental attribute of landscapes, intrinsically linked to drainage density, landslide hazard, biogeochemical cycling and hillslope sediment transport. Existing methods to estimate catchment average hillslope lengths include inversion of drainage density or identification of a break in slope–area scaling, where the hillslope domain transitions into the fluvial domain. Here we implement a technique which models flow from point sources on hilltops across pixels in a digital elevation model (DEM), based on flow directions calculated using pixel aspect, until reaching the channel network, defined using recently developed channel extraction algorithms. Through comparisons between these measurement techniques, we show that estimating hillslope length from plots of topographic slope versus drainage area, or by inverting measures of drainage density, systematically underestimates hillslope length. In addition, hillslope lengths estimated by slope–area scaling breaks show large variations between catchments of similar morphology and area. We then use hillslope length–relief structure of landscapes to explore nature of sediment flux operating on a landscape. Distinct topographic forms are predicted for end‐member sediment flux laws which constrain sediment transport on hillslopes as being linearly or nonlinearly dependent on hillslope gradient. Because our method extracts hillslope profiles originating from every ridgetop pixel in a DEM, we show that the resulting population of hillslope length–relief measurements can be used to differentiate between linear and nonlinear sediment transport laws in soil mantled landscapes. We find that across a broad range of sites across the continental United States, topography is consistent with a sediment flux law in which transport is nonlinearly proportional to topographic gradient. © 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Variation in the sediment deposition behind check‐dams under different soil erosion conditions on the Loess Plateau,China

To maintain a reasonable sediment regulation system in the middle reaches of the Yellow River, it is critical to determine the variation in sediment deposition behind check‐dams for different soil erosion conditions. Sediment samples were collected by using a drilling machine in the Fangta watershed of the loess hilly–gully region and the Manhonggou watershed of the weathered sandstone hilly–gully (pisha) region. On the basis of the check‐dam capacity curves, the soil bulk densities and the couplet thickness in these two small watersheds, the sediment yields were deduced at the watershed scale. The annual average sediment deposition rate in the Manhonggou watershed (702.0 mm/(km²·a)) from 1976 to 2009 was much higher than that in the Fangta watershed (171.6 mm/(km²·a)) from 1975 to 2013. The soil particle size distributions in these two small watersheds were generally centred on the silt and sand fractions, which were 42.4% and 50.7% in the Fangta watershed and 60.6% and 32.9% in the Manhonggou watershed, respectively. The annual sediment deposition yield exhibited a decreasing trend; the transition years were 1991 in the Fangta watershed and 1996 in the Manhonggou watershed (P < 0.05). In contrast, the annual average sediment deposition yield was much higher in the Manhonggou watershed (14011.1 t/(km²·a)) than in the Fangta watershed (3149.6 t/(km²·a)). In addition, the rainfalls that induced sediment deposition at the check‐dams were greater than 30 mm in the Fangta watershed and 20 mm in the Manhonggou watershed. The rainfall was not the main reason for the difference in the sediment yield between the two small watersheds. The conversion of farmland to forestland or grassland was the main reason for the decrease in the soil erosion in the Fangta watershed, while the weathered sandstone and bare land were the main factors driving the high sediment yield in the Manhonggou watershed. Knowledge of the sediment deposition process of check‐dams and the variation in the catchment sediment yield under different soil erosion conditions can serve as a basis for the implementation of improved soil erosion and sediment control strategies, particularly in semi‐arid hilly–gully regions. Copyright © 2018 John Wiley & Sons, Ltd.     

Base flow separation for soil erosion simulation in a granitic forested headwater catchment using a process-based model,GeoWEPP

Distributed erosion models, which simulate the physical processes of water flow and soil erosion, are effective for predicting soil erosion in forested catchments. Although subsurface flow through multiple pathways is dominant for runoff generation in forested headwater catchments, the process-based erosion model, Geo-spatial interface for Water Erosion Prediction Project(Geo WEPP), does not have an adequate subsurface component for the simulation of hillslope water flow. In the current study, t...     

Use of fallout tracers 7Be, 210Pb and 137Cs to distinguish the form of sub‐surface soil erosion delivering sediment to rivers in large catchments

Fallout radionuclides (FRNs) ¹³⁷Cs and ²¹⁰Pb are well established as tracers of surface and sub‐surface soil erosion contributing sediment to river systems. However, without additional information, it has not been possible to distinguish sub‐surface soil erosion sources. Here, we use the FRN ⁷Be (half‐life 53 days) in combination with ¹³⁷Cs and excess ²¹⁰Pb to trace the form of erosion contributing sediment in three large river catchments in eastern Australia; the Logan River (area 3700 km²), Bowen River (9400 km²) and Mitchell River (4700 km²). We show that the combination of ¹³⁷Cs, excess ²¹⁰Pb and ⁷Be can discriminate horizontally aligned sub‐surface erosion sources (rilled and scalded hillslopes and the floors of incised drainage lines and gully ‘badland’ areas) from vertical erosion sources (channel banks and gully walls). Specifically, sub‐surface sources of sediment eroded during high rainfall and high river flow events have been distinguished by the ability of rainfall‐derived ⁷Be to label horizontal soil surfaces, but not vertical. Our results indicate that in the two northern catchments, erosion of horizontal sub‐surface soil sources contributed almost as much fine river sediment as vertical channel banks, and several times the contribution of hillslope topsoils. This result improves on source discrimination provided previously and indicates that in some areas erosion of hillslope soils may contribute significantly to sediment yield, but not as topsoil loss. We find that in north‐eastern Australia, scalded areas on hillslopes and incising drainage lines may be sediment sources of comparable importance to vertical channel banks. Previous studies have used the combination of ¹³⁷Cs, excess ²¹⁰Pb and ⁷Be to estimate soils losses at the hillslope scale. Here, we show that with timely and judicious sampling of soil and sediment during and immediately after high flow events ⁷Be measurements can augment fallout ¹³⁷Cs and ²¹⁰Pb to provide important erosion source information over large catchments. Copyright © 2013 John Wiley & Sons, Ltd.     

Sediment delivery from hillslopes and the Universal Soil Loss Equation: some perceptions and misconceptions

The Universal Soil Loss Equation (USLE) or the revised USLE (RUSLE) are often used together with sediment delivery ratios in order to predict sediment delivery from hillslopes. In using sediment delivery ratios for this purpose, it is assumed that the sediment delivery ratio for a given hillslope does not vary with the amount of erosion occurring in the upslope area. This assumption is false. There is a perception that hillslope erosion is calculated on the basis that hillslopes are, in effect, simply divided into 22·1 m long segments. This perception fails to recognize the fact the inclusion of the 22·1 m length in the calculation has no physical significance but simply produces a value of 1·0 for the slope length factor when slopes have a length equal to that of the unit plot. There is a perception that the slope length factor is inappropriate because not all the dislodged sediment is discharged. This perception fails to recognize that the USLE and the RUSLE actually predict sediment yield from planar surfaces, not the total amount of soil material dislocated and removed some distance by erosion within an area. The application of the USLE/RUSLE to hillslopes also needs to take into account the fact that runoff may not be generated uniformly over that hillslope. This can be achieved by an equation for the slope length factor that takes account of spatial variations in upslope runoff on soil loss from a segment or grid cell. Several alternatives to the USLE event erosivity index have been proposed in order to predict event erosion better than can be achieved using the EI₃₀ index. Most ignore the consequences of changing the event erosivity index on the values for the soil, crop and soil conservation protection factors because there is a misconception that these factors are independent of one another. Copyright © 2007 John Wiley & Sons, Ltd.     

ESTIMATION OF SOIL EROSION IN A RESERVOIR WATERSHED USING ^137CS FALLOUT RADIONUCLIDE

Sedimentation from soil erosion is a critical reservoir watershed management issue. Due to the difficulty of field investigations, empirical formulas are commonly used to estimate the soil erosion rate. However, these estimations are often far from accurate. An effective alternative to estimating soil erosion is to analyze the spatial variation of 137Cs inventory in the soil. 137Cs can be adsorbed by the soil and is widely assumed to change its distribution only when disturbed by rainfall and human activities. Thus, 137Cs distributed in soils can be a useful environmental tracer to estimate soil erosion. In this study, the net soil loss estimate is 108,346 t/yr and the gross erosion and net erosion rates are 10.1 and 9 t/ha yr respectively. The sediment delivery ratio is therefore estimated to be 0.9 based on the two erosion rates. Because of the steep hillsides in the watershed, only 10% of the sediment yield stayed in the deposition sites and 90% was transported to the river as the sediment output. Soil erosion estimates from spatial variations of the 137Cs activity in the Baishi river watershed showed satisfactory accuracy when compared to sediment yield data. Using soil 137Cs concentrations is therefore a feasible method for estimating soil loss or deposition in Taiwan. Data sampling, analysis and result of this approach are given in this paper.     

PHYSICAL PROCESS BASED SOIL EROSION MODEL IN A SMALL WATERSHED IN THE HILLY LOESS REGION

ＰＨＹＳＩＣＡＬＰＲＯＣＥＳＳＢＡＳＥＤＳＯＩＬＥＲＯＳＩＯＮＭＯＤＥＬＩＮＡＳＭＡＬＬＷＡＴＥＲＳＨＥＤＩＮＴＨＥＨＩＬＬＹＬＯＥＳＳＲＥＧＩＯＮ１ＣＡＩＱｉａｎｇｇｕｏ２ＡＢＳＴＲＡＣＴＡｐｈｙｓｉｃａｌｐｒｏｃｅｓｂａｓｅｄｐｅｒｓｔｏｒｍ...     

The design of post‐mining landscapes using geomorphic principles

Nature can provide analogues for post‐mining landscapes in terms of landscape stability and also in terms of the rehabilitated structure ‘blending in’ with the surrounding undisturbed landscape. In soil‐mantled landscapes, hillslopes typically have a characteristic pro?le that has a convex upper hillslope pro?le with a concave pro?le lower down the slope. In this paper hillslope characteristic form is derived using the area–slope relationship from pre‐mining topography at two sites in Western Australia. Using this relationship, concave hillslope pro?les are constructed and compared to linear hillslopes in terms of sediment loss using the SIBERIA erosion model. It is found that concave hillslopes can reduce sediment loss by up to ?ve times that of linear slopes. Concave slopes can therefore provide an alternative method for the construction of post‐mining landscapes. An understanding of landscape geomorphological properties and the use of erosion models can greatly assist in the design of post‐mining landscapes. Copyright © 2003 John Wiley & Sons, Ltd.