Water flux and sediment transport within a forested landscape: the role of connectivity,subsurface flow,and slope length scale on transport mechanism

The connectivity and upscaling of overland runoff and sediment transport are important issues in hillslope hydrology to identify water flux and sediment transport within landscape. These processes are highly variable in time and space with regard to their interactions with vegetation and soil surface conditions. The generation of overland runoff and its spatial connectivity were examined along a slope to determine the variations in the transport mechanism of runoff and soil particles by rain splash and overland runoff. Field experiments were conducted by erosion plots on a steep hillslope at lengths of 5, 10, and 15 m. The overland runoff connectivity and flow transport distance decreased with the slope length, while spatial variability of infiltration increased significantly with the slope length. Observation of subsurface flow revealed that surface soil and litter layer could have important role in water transport. However, the surface soil water content and water flux transport along the slope was highly variable for different storm events; the variability was related to the complexity of the system, mainly by way of the initial wetness conditions and infiltration characteristics. Only net rain‐splashed soil was measurable, but examination of the water flux, overland runoff and sediment transport connectivity, characteristics of sheetwash, and the variability in spatial infiltration indicated an increase in the contribution of the rain splash transport mechanism along the slope. Copyright © 2013 John Wiley & Sons, Ltd.     

Field studies of sediment transport by overland flow

Field studies on sandy soils of the Cottenham Series in mid-Bedfordshire show that the mean annual rate of sediment transport by overland flow on an 11° mid-slope is 98 g cm¹. The feasibility of using sediment transport equations to predict erosion by overland flow on a storm basis is examined by comparing the observed values of sediment yield with values predicted by four sediment transport equations and a regression equation which relates soil loss to runoff energy and rainfall energy. An expression combining Engelund's sediment transport capacity equation and the Manning equation for flow velocity, as modified by Savat for disturbed flow, best reflects field conditions. Although there is a significant correlation (r = 0.69; n 30) between the observed and predicted values using this expression, the coefficient of determination is too low for predictive purposes. Reasons for this are presented.     

Soil erosion and pollutant transport during rainfall-runoff processes

The pollutant from land surface applied to agricultural chemicals is one of the major sources of contamination in water bodies. The pollutant transport within a watershed is profoundly influenced by the rainfall-runoff processes, especially the associated upland erosion and sediment transport processes because most of pollutant can be dissolved into water or attached to the soil particles. A set of soil experiments in laboratory was conducted in this paper to investigate the impacts of upland erosion and sediment transport on pollutant loads. The soil utilized for the experiments was the silty sand collected from Loess Plateau, China; and ammonium bicarbonate was applied on the soil surface as the pollutant source. Runoff discharge, soil loss, and ammonia- and nitrate-nitrogen concentrations were measured to establish the relationships which can help the numerical model to predict the pollutant losses coupled with upland soil erosion during the rain-fall-runoff processes. The experimental results indicate the ammonia-nitrogen concentration in runoff reaches the peak at the initial stage of the overland flow generation, and quickly decreases and approaches to the steady state. The ammonia-nitrogen transported by the soil loss also makes contributions to the nitrogen loss; and its amount mainly depends on the soil transport rate. The ammonia-nitrogen dissolved in overland flow is dominant due to the strong aqueous solution of ammonium bicarbonate during the first storm right after its application.     

Interstorm surface preparation and sediment detachment by vehicle traffic on unpaved mountain roads

Road survey and field rainfall simulation experiments have shown that the erodibility of a road surface is dynamic. In the absence of extreme runoff events, dynamic erodibility results from the generation and removal of easily entrained surface material by human road surface maintenance activities, vehicular detachment and overland flow events. Maintenance activities introduce easily transportable material to the road surface where it can be entrained by overland flow. Traffic in dry conditions detaches material that is quickly removed during subsequent overland flow events. The pre‐storm erodibility of a road is therefore largely a function of maintenance and vehicle activity since the last overland flow event. During rainstorms, vehicle passes increase sediment production by detaching/redistributing surface material and creating efficient overland flow pathways for sediment transport. However, if incision of tracks by overland flow does not occur, post‐pass sediment transport quickly returns to pre‐pass rates. Field rainfall simulation data suggest that sediment transport resulting from during‐storm vehicle passes is greatly influenced by the presence of existing loose material, which again is a function of prior road usage and maintenance activities. Incorporation of vehicular effects into physically based road erosion models may be possible by parameterizing both during‐storm and inter‐storm changes in the supply of loose surface material as changes in surface erodibility. Copyright © 2001 John Wiley & Sons, Ltd.     

Downslope soil detachment–transport on steep slopes via rain splash

This study developed a one‐dimensional model of downslope rain splash transport based on field experiments and previous studies. The developed model considers soil detachment processes, ground cover, probability densities, and the effect of overland run‐off in preventing detachment. Field monitoring was conducted to observe precipitation run‐off, ground cover, and sediment production on steep hillslopes. Field‐observed data were used to develop the splash detachment rate equation, probability densities for splash transport, and the maximum splash transport distance. Observed and estimated splash transport showed overall agreement, with some differences for small storm events or events with relatively low intensity, probably caused by variation of overland run‐off depth and connectivity as well as differences in soil surface cohesion at various degrees of wetness. Our model can provide insights on the interactions among rainfall intensity, soil surface condition, soil wetness, and splash transport on forested hillslopes. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling water erosion in the Sahel: application of a physically based soil erosion model in a gentle sloping environment

Water is a major limiting factor in arid and semi‐arid agriculture. In the Sahelian zone of Africa, it is not always the limited amount of annual rainfall that constrains crop production, but rather the proportion of rainfall that enters the root zone and becomes plant‐available soil moisture. Maximizing the rain‐use efficiency and therefore limiting overland flow is an important issue for farmers. The objectives of this research were to model the processes of infiltration, runoff and subsequent erosion in a Sahelian environment and to study the spatial distribution of overland flow and soil erosion. The wide variety of existing water erosion models are not developed for the Sahel and so do not include the unique Sahelian processes. The topography of the Sahelian agricultural lands in northern Burkina Faso is such that field slopes are generally low (0–5°) and overland flow mostly occurs in the form of sheet flow, which may transport large amounts of fine, nutrient‐rich particles despite its low sediment transport capacity. Furthermore, pool formation in a field limits overland flow and causes resettlement of sediment resulting in the development of a surface crust. The EUROSEM model was rewritten in the dynamic modelling code of PCRaster and extended to account for the pool formation and crust development. The modelling results were calibrated with field data from the 2001 rainy season in the Katacheri catchment in northern Burkina Faso. It is concluded that the modified version of EUROSEM for the Sahel is a fully dynamic erosion model, able to simulate infiltration, runoff routing, pool formation, sediment transport, and erosion and deposition by inter‐rill processes over the land surface in individual storms at the scale of both runoff plots and fields. A good agreement is obtained between simulated and measured amounts of runoff and sediment discharge. Incorporating crust development during the event may enhance model performance, since the process has a large influence on infiltration capacity and sediment detachment in the Sahel. Copyright © 2005 John Wiley & Sons, Ltd.     

An analytical approach to ascertain saturation‐excess versus infiltration‐excess overland flow in urban and reference landscapes

Uncontrolled overland flow drives flooding, erosion, and contaminant transport, with the severity of these outcomes often amplified in urban areas. In pervious media such as urban soils, overland flow is initiated via either infiltration‐excess (where precipitation rate exceeds infiltration capacity) or saturation‐excess (when precipitation volume exceeds soil profile storage) mechanisms. These processes call for different management strategies, making it important for municipalities to discern between them. In this study, we derived a generalized one‐dimensional model that distinguishes between infiltration‐excess overland flow (IEOF) and saturation‐excess overland flow (SEOF) using Green–Ampt infiltration concepts. Next, we applied this model to estimate overland flow generation from pervious areas in 11 U.S. cities. We used rainfall forcing that represented low‐ and high‐intensity events and compared responses among measured urban versus predevelopment reference soil hydraulic properties. The derivation showed that the propensity for IEOF versus SEOF is related to the equivalence between two nondimensional ratios: (a) precipitation rate to depth‐weighted hydraulic conductivity and (b) depth of soil profile restrictive layer to soil capillary potential. Across all cities, reference soil profiles were associated with greater IEOF for the high‐intensity set of storms, and urbanized soil profiles tended towards production of SEOF during the lower intensity set of storms. Urban soils produced more cumulative overland flow as a fraction of cumulative precipitation than did reference soils, particularly under conditions associated with SEOF. These results will assist cities in identifying the type and extent of interventions needed to manage storm water produced from pervious areas.     

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.     

Evaluating overland flow sediment transport capacity

An equation for evaluating the sediment transport capacity of overland flow is a necessary part of a physically based soil erosion model describing sediment detachment and transport as distributed processes. At first, for the hydraulic conditions of small-scale and large-scale roughness, the sediment transport capacity relationship used in the WEPP model is calibrated by Yalin and Govers' equation. The analysis shows that the transport coefficient K_t depends on the Shields parameter, Y, according to a semi-logarithmic (Yalin) or a linear (Govers) equation. The reliability of the semi-logarithmic equation is verified by Smart's, and Aziz and Scott's experimental data. Then the Low's formula, whose applicability is also proved by Smart's, and Aziz and Scott's data, is transformed as a stream power equation in which a stream power coefficient, K_SP, depending on Shields parameter, slope, sediment and water-specific weight, appears. A relationship between transport capacity and effective stream power is also proposed. Finally, the influence of rainfall on sediment transport capacity and the prediction of critical shear stress corresponding to overland flow are examined. © 1998 John Wiley & Sons, Ltd.     

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.     

On predicting upland erosion losses from rainfall depth

Point rainfall triggers the complex processes of overland flow and surface erosion. The probability density functions of rainfall duration and intensity are coupled with a physically based dynamic formulation of rainfall-runoff-sediment transport relationships for upland areas. When considering a single storm, rainfall depth alone is a poor predictor of sediment transport because of the dispersion introduced by the effect of rainfall intensity. On a long terms basis, however, the total amount of rainfall can be used to predict total erosion losses.     

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.     

A MATHEMATICAL QUASI-2D OVERLAND FLOW MODEL

1 INTRODUCTIONThe prediction of future impacts on terrestrial ecosystems by atmospheric, climatic and land-usechanges is the aim of watershed management. Meeting these requirements scientists, managers and policymakers try to achieve the sustainable management of the vitally important resources of watersheds due toan integrated ecosystem approach at the catchment scale. As composite landscapes often have a highdegree of contingency between its elements, the transport over these landscape s…     

Incremental distributed modelling investigation in a small agricultural catchment: 2. Erosion and phosphorus transport

Distributed physically based erosion and phosphorus (P) transport models, run by the overland flow model described in Taskinen and Bruen (2006. Hydrological Processes 20 : this issue), are described. In the erosion model, the additional components to the basic model were the outflow of the particles by infiltration and a new model component, i.e. deposition when rainfall stops. Two ways of calculating the shielding factor due to the flow depth were compared. The P transport model had both dissolved P (DP) and particulate P (PP) components. The processes included in the DP model were desorption from the soil surface, advection, storage in the overland flow and infiltration. The PP model accounted for advection, storage in the flow, infiltration, detachment from the soil surface by flow and rainfall and deposition both when transport capacity of suspended solids (SS) is exceeded and when rainfall ceases. When the models were developed and validated in small agricultural fields of cohesive soil types in southern Finland, comparisons were made between corresponding processes and the significance of added components were estimated in order to find out whether increased model complexity improves the model performance. The sedigraphs were found to follow the dynamics of rainfall, emphasizing the importance of the rainfall splash component. The basic model was too slow to react to changes in rainfall and flow rates, but infiltration and deposition that acts during the cessation in rainfall improved the model significantly by enabling the modelled SS to fall sharply enough. The shielding effect of flow depth from the splash detachment was found to play a significant role. Transport capacity should also be included in erosion models when they are applied to cohesive soils. In this study, the Yalin method worked well. A strong correlation was obvious between the measured SS and total P concentrations, indicating that the main form of P in runoff is PP. This emphasizes the importance of a good sediment transport model in P transport modelling. The submodel used for DP desorption from the soil surface produced plausible results without any calibration. Copyright © 2006 John Wiley & Sons, Ltd.     

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

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

Field investigation and modelling of coupled stream discharge and shallow water‐table dynamics in a small Mediterranean catchment (Sardinia)

In the semi‐arid Mediterranean environment, the rainfall–runoff relationships are complex because of the markedly irregular patterns in rainfall, the seasonal mismatch between evaporation and rainfall, and the spatial heterogeneity in landscape properties. Watersheds often display considerable non‐linear threshold behavior, which still make runoff generation an open research question. Our objectives in this context were: to identify the primary processes of runoff generation in a small natural catchment; to test whether a physically based model, which takes into consideration only the primary processes, is able to predict spatially distributed water‐table and stream discharge dynamics; and to use the hydrological model to increase our understanding of runoff generation mechanisms. The observed seasonal dynamics of soil moisture, water‐table depth, and stream discharge indicated that Hortonian overland‐flow was negligible and the main mechanism of runoff generation was saturated subsurface‐flow. This gives rise to base‐flow, controls the formation of the saturated areas, and contributes to storm‐flow together with saturation overland‐flow. The distributed model, with a 1D scheme for the kinematic surface‐flow, a 2D sub‐horizontal scheme for the saturated subsurface‐flow, and ignoring the unsaturated flow, performed efficiently in years when runoff volume was high and medium, although there was a smoothing effect on the observed water‐table. In dry years, small errors greatly reduced the efficiency of the model. The hydrological model has allowed to relate the runoff generation mechanisms with the land‐use. The forested hillslopes, where the calibrated soil conductivity was high, were never saturated, except at the foot of the slopes, where exfiltration of saturated subsurface‐flow contributed to storm‐flow. Saturation overland‐flow was only found near the streams, except when there were storm‐flow peaks, when it also occurred on hillslopes used for pasture, where soil conductivity was low. The bedrock–soil percolation, simulated by a threshold mechanism, further increased the non‐linearity of the rainfall–runoff processes. Copyright © 2013 John Wiley & Sons, Ltd.     

MODELING EPHEMERAL GULLY EROSION FOR CONSERVATION PLANNING

1INTRODUCTIONEphemeral gully erosion,which is caused by concentrated flow within cultivated farm fields,is distinct from rill erosion.Ephemeral gully erosion is also distinct from gully erosion in permanent,deep,incised channels,formed by headcuts moving upstream.Ephemeral gully erosion is often overlooked.It is not estimated with rill-interrill erosion prediction technology such as the Revised Universal Soil Loss Equation(Renard et al.,1997),and it is often not measured in field survey…     

Estimation of soil erosion and sediment yield using GIS

Abstract

A Geographical Information System (GIS) based method is proposed and demonstrated for the identification of sediment source areas and the prediction of storm sediment yield from catchments. Data from the Nagwa and Karso catchments in Bihar (India) have been used. The Integrated Land and Water Information System (ILWIS) GIS package has been used for carrying out geographic analyses. An Earth Resources Data Analysis System (ERDAS) Imagine image processor has been used for the digital analysis of satellite data for deriving the land cover and soil characteristics of the catchments. The catchments were discretized into hydrologically homogeneous grid cells to capture the catchment heterogeneity. The cells thus formed were then differentiated into cells of overland flow regions and cells of channel flow regions based on the magnitude of their flow accumulation areas. The gross soil erosion in each cell was calculated using the Universal Soil Loss Equation (USLE) by carefully determining its various parameters. The concept of sediment delivery ratio (SDR) was used for determination of the total sediment yield of each catchment during isolated storm events.