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.     

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

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

Assessment of soil erosion in a monsoon-dominated mountain river basin in India using RUSLE-SDR and AHP

The long-term average annual soil loss (A) and sediment yield (SY) in a tropical monsoon-dominated river basin in the southern Western Ghats, India (Muthirapuzha River Basin, MRB; area: 271.75 km²), were predicted by coupling the Revised Universal Soil Loss Equation (RUSLE) and sediment delivery ratio (SDR) models. Moreover, the study also delineated soil erosion risk zones based on the soil erosion potential index (SEPI) using the analytic hierarchy process (AHP) technique. Mean A of the basin is 14.36 t ha^?1 year^?1, while mean SY is only 3.65 t ha^?1 year^?1. Although the land use/land cover types with human interference show relatively lower A compared to natural vegetation, their higher SDR values reflect the significance of anthropogenic activities in accelerated soil erosion. The soil erosion risk in the MRB is strongly controlled by slope, land use/land cover and relative relief, compared to geomorphology, drainage density, stream frequency and lineament frequency.     

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.     

The contribution of aeolian processes to fluvial sediment yield from a desert watershed in the Ordos Plateau,China

In arid zones, many active aeolian dunes terminate at ephemeral and perennial desert rivers. The desert rivers show very high rates of sediment transport that cause deleterious downstream effects on the river system and ecology. High sediment loading has been attributed to severe water erosion of sparsely covered watersheds during infrequent but heavy rainfall. Although aeolian erosion is known to lead to high rates of wind‐blown sand transport, direct confirmation of whether the aeolian processes accelerate or inhibit fluvial sediment loss is lacking. Here, we show that an aeolian‐fluvial cycling process is responsible for the high rate of suspended sediment transport in a Sudalaer ephemeral desert channel in the Ordos Plateau of China. Frequent aeolian processes, but low frequency (once every 3 years on average) flooding, occur in this region. Wind‐blown saltating grains appeared to be unable to cross the desert channel because of interruption of channel‐induced recirculating air flow, and therefore tended to settle in the channel during the windy seasons, leading to channel narrowing. During flooding, this narrowed channel was found to yield a threefold increase in suspended sediment loading and a 3.4‐fold increase in the weight percentage of the 0.08–0.2 mm sediment fraction on 18 July 2012. Loss of stored aeolian sand due to channel erosion accounted for about half of the total sediment yield in this watershed. These findings show that aeolian processes play an essential role in accelerating the sediment yield from a watershed characterized by aeolian‐fluvial interplay and also suggest that the drier the region and the greater the aeolian process, the more the aeolian process contributes to fluvial sediment yield. Copyright © 2013 John Wiley & Sons, Ltd.     

Evaluating the effect of land use changes on soil erosion and sediment yield using a grid‐based distributed modelling approach

Processes of soil erosion and sediment transport are strongly influenced by land use changes so the modelling of land use changes is important with respect to the simulation of soil degradation and its on‐site and off‐site consequences. The reliability of simulation results from erosion models is circumscribed by considerable spatial variation in many parameters. However, most of the currently widely used erosion models at the mesoscale are semidistributed, which leads to difficulties in incorporating a high degree of spatial information, especially land use information, so that the effects of land use changes on soil erosion have hitherto not been investigated in detail using these models. In this article, a grid‐based distributed erosion and sediment transport model is introduced, which simulates the spatial pattern of erosion and deposition rates and sediment transport processes in river channels. In this model, land use affects soil erosion through altering soil loss and influencing sediment delivery. Simulated soil erosion for events recorded in 1989 and 1996 in the Lushi basin in China was analyzed by comparing it with historical land use maps. The results indicated that even relatively minor land use changes had a significant effect on regional soil erosion rates and sediment transport to rivers. The average erosion rate increased from 1989 to 1996, after the transformation of forest to farmland. The results of the study suggest that the proposed soil erosion model can be applied in similar river basins. Copyright © 2011 John Wiley & Sons, Ltd.     

SCALES: a large‐scale assessment model of soil erosion hazard in Basse‐Normandie (northern‐western France)

The cartography of erosion risk is mainly based on the development of models, which evaluate in a qualitative and quantitative manner the physical reproduction of the erosion processes (CORINE, EHU, INRA). These models are mainly semi‐quantitative but can be physically based and spatially distributed (the Pan‐European Soil Erosion Risk Assessment, PESERA). They are characterized by their simplicity and their applicability potential at large temporal and spatial scales. In developing our model SCALES (Spatialisation d'éChelle fine de l'ALéa Erosion des Sols/large‐scale assessment and mapping model of soil erosion hazard), we had in mind several objectives: (1) to map soil erosion at a regional scale with the guarantee of a large accuracy on the local level, (2) to envisage an applicability of the model in European oceanic areas, (3) to focus the erosion hazard estimation on the level of source areas (on‐site erosion), which are the agricultural parcels, (4) to take into account the weight of the temporality of agricultural practices (land‐use concept). Because of these objectives, the nature of variables, which characterize the erosion factors and because of its structure, SCALES differs from other models. Tested in Basse‐Normandie (Calvados 5500 km²) SCALES reveals a strong predisposition of the study area to the soil erosion which should require to be expressed in a wet year. Apart from an internal validation, we tried an intermediate one by comparing our results with those from INRA and PESERA. It appeared that these models under estimate medium erosion levels and differ in the spatial localization of areas with the highest erosion risks. SCALES underlines here the limitations in the use of pedo‐transfer functions and the interpolation of input data with a low resolution. One must not forget however that these models are mainly focused on an interregional comparative approach. Therefore the comparison of SCALES data with those of the INRA and PESERA models cannot result on a convincing validation of our model. For the moment the validation is based on the opinion of local experts, who agree with the qualitative indications delivered by our cartography. An external validation of SCALES is foreseen, which will be based on a thorough inventory of erosion signals in areas with different hazard levels. Copyright © 2010 John Wiley & Sons, Ltd.     

Determination of runoff and sediment yield from a small watershed in sub‐humid subtropics using the HSPF model

The Hydrologic Simulation Programme‐Fortran (HSPF), a hydrologic and water quality computer model, was employed for simulating runoff and sediment yield during the monsoon months (June–October) from a small watershed situated in a sub‐humid subtropical region of India. The model was calibrated using measured runoff and sediment yield data for the monsoon months of 1996 and was validated for the monsoon months of 2000 and 2001. During the calibration period, daily‐calibrated runoff had a Nash‐Sutcliffe efficiency (E_NS) value of 0·68 and during the validation period it ranged from 0·44 to 0·67. For daily sediment yield E_NS was 0·71 for the calibration period and it ranged from 0·68 to 0·90 for the validation period. Sensitivity analysis was performed to assess the impact of important watershed characteristics. The model parameters obtained in this study could serve as reference values for model application in similar climatic regions, with practical implications in watershed planning and management and designing best management practices. Copyright © 2007 John Wiley & Sons, Ltd.     

Modeling the sediment yield and the impact of vegetated filters using an event‐based soil erosion model—a case study of a small Canadian watershed

This paper presents the first application of the event‐based MHYDAS‐Erosion model to a small agricultural watershed under temperate climate conditions (Quebec, Canada). Simulation results based on observed and synthetic rainfall events revealed a bimodal behaviour of sediment yield. During high‐intensity rainfall events, most of the sediments reaching the watershed outlet originate from cropland. Meanwhile, during low‐intensity events, most of the sediments come from the drainage network. Furthermore, simulation results show that implementation of 5‐m and 20‐m wide vegetated filters throughout the watershed or at the edge of the most problematic fields (4% of the total fields) could reduce soil loss by 52% and 31%, respectively. The modeling framework could be used for the design and location of beneficial management practices such as grass strips and riparian zones Copyright © 2016 John Wiley & Sons, Ltd.     

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

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

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.     

Using GIS and a digital elevation model to assess the effectiveness of variable grade flow diversion terraces in reducing soil erosion in northwestern New Brunswick,Canada

Flow diversion terraces (FDT) are commonly used beneficial management practice (BMP) for soil conservation on sloped terrain susceptible to water erosion. A simple GIS‐based soil erosion model was designed to assess the effectiveness of the FDT system under different climatic, topographic, and soil conditions at a sub‐basin level. The model was used to estimate the soil conservation support practice factor (P‐factor), which inherently considered two major outcomes with its implementation, namely (1) reduced slope length, and (2) sediment deposition in terraced channels. A benchmark site, the agriculture‐dominated watershed in northwestern New Brunswick (NB), was selected to test the performance of the model and estimated P‐factors. The estimated P‐factors ranged from 0·38–1·0 for soil conservation planning objectives and ranged from 0·001 to 0·45 in sediment yield calculations for water‐quality assessment. The model estimated that the average annual sediment yield was 773 kg ha^?1 yr ^?1 compared with a measured value of 641 kg ha^?1 yr^?1. The P‐factors estimated in this study were comparable with predicted values obtained with the revised universal soil loss equation (RUSLE2). The P‐factors from this study have the potential to be directly used as input in hydrological models, such as the soil and water assessment tool (SWAT), or in soil conservation planning where only conventional digital elevation models (DEMs) are available. Copyright © 2009 John Wiley & Sons, Ltd.     

Evaluating the use of a gridded climate surface for modelling groundwater recharge in a semi‐arid region (Okanagan Basin,Canada)

Spatially distributed groundwater recharge was simulated for a segment of a semi‐arid valley using three different treatments of meteorological input data and potential evapotranspiration (PET). For the same area, timeframe, land cover characteristics and soil properties, groundwater recharge was estimate using (i) single‐station climate data with monthly PET calculated by the Thornthwaite method; (ii) single‐station climate data with daily PET calculated by the Penman–Monteith method; and (iii) daily gridded climate data with spatially distributed PET calculated using the Penman–Monteith method. For each treatment, the magnitude and distribution of actual evapotranspiration (AET) for summer months compared well with those estimated for a 5‐year crop study, suggesting that the near‐surface hydrological processes were replicated and that subsequent groundwater recharge rates are realistic. However, for winter months, calculated AET was near zero when using the Thornthwaite PET method. Mean annual groundwater recharge varied from ～3·2 to 10·0 mm when PET was calculated by the Thornthwaite method, and from ～1·8 to 7·5 mm when PET was calculated by the Penman–Monteith method. Comparisons of bivariate plots of seasonal recharge rates estimated from single‐station versus gridded surface climate reveal that there is greater variability between the different methods for spring months, which is the season of greatest recharge. Furthermore, these seasonal differences are shown to provide different results when compared to the depth to water table, which could lead to different results of evaporative extinction depth. These findings illustrate potential consequences of using different approaches for representing spatial meteorological input data, which could provide conflicting predictions when modelling the influence of climate change on groundwater recharge. Copyright © 2010 John Wiley & Sons, Ltd.     

Season‐based rainfall–runoff modelling using the probability‐distributed model (PDM) for large basins in southeastern Brazil

Southeastern Brazil is characterized by seasonal rainfall variability. This can have a great social, economic, and environmental impact due to both excessive and deficient water availability. During 2014 and 2015, the region experienced one of the most severe droughts since 1960. The resulting water crisis has seriously affected water supply to the metropolitan region of São Paulo and hydroelectric power generation throughout the entire country. This research considered the upstream basins of the southeastern Brazilian reservoirs Cantareira (2,279 km²; water supply) and Emborcação (29,076 km²), Três Marias (51,576 km²), Furnas (52,197 km²), and Mascarenhas (71,649 km²; hydropower) for hydrological modelling. It made the first attempt at configuring a season‐based probability‐distributed model (PDM‐CEMADEN) for simulating different hydrological processes during wet and dry seasons. The model successfully reproduced the intra‐annual and interannual variability of the upstream inflows during 1985–2015. The performance of the model was very satisfactory not only during the wet, dry, and transitional seasons separately but also during the whole period. The best performance was obtained for the upstream basin of Furnas, as it had the highest quality daily precipitation and potential evapotranspiration data. The Nash–Sutcliffe efficiency and logarithmic Nash–Sutcliffe efficiency were 0.92 and 0.93 for the calibration period 1984–2001, 0.87 and 0.88 for the validation period 2001–2010, and 0.93 and 0.90 for the validation period 2010–2015, respectively. Results indicated that during the wet season, the upstream basins have a larger capacity and variation of soil water storage, a larger soil water conductivity, and quicker surface water flow than during the dry season. The added complexity of configuring a season‐based PDM‐CEMADEN relative to the traditional model is well justified by its capacity to better reproduce initial conditions for hydrological forecasting and prediction. The PDM‐CEMADEN is a simple, efficient, and easy‐to‐use model, and it will facilitate early decision making and implement adaptation measures relating to disaster prevention for reservoirs with large‐sized upstream basins.     

The influence of preferential flow on hillslope hydrology in a semi‐arid watershed (in the Spanish Dehesas)

Preferential flow is known to influence hillslope hydrology in many areas around the world. Most research on preferential flow has been performed in temperate regions. Preferential infiltration has also been found in semi‐arid regions, but its impact on the hydrology of these regions is poorly known. The aim of this study is to describe and quantify the influence of preferential flow on the hillslope hydrology from small scale (infiltration) to large scale (subsurface stormflow) in a semi‐arid Dehesa landscape. Precipitation, soil moisture content, piezometric water level and discharge data were used to analyse the hydrological functioning of a catchment in Spain. Variability of soil moisture content during the transition from dry to wet season (September to November) within horizontal soil layers leads to the conclusion that there is preferential infiltration into the soils. When the rainfall intensity is high, a water level rapidly builds up in the piezometer pipes in the area, sometimes even reaching soil surface. This water level also drops back to bedrock within a few hours (under dry catchment conditions) to days (under wet catchment conditions). As the soil matrix is not necessarily wet while this water layer is built up, it is thought to be a transient water table in large connected pores which drain partly to the matrix, partly fill up bedrock irregularities and partly drain through subsurface flow to the channels. When the soil matrix becomes wetter the loss of water from macropores to the matrix and bedrock decreases and subsurface stormflow increases. It may be concluded that the hillslope hydrological system consists of a fine matrix domain and a macropore domain, which have their own flow characteristics but which also interact, depending on the soil matrix and macropore moisture contents. The macropore flow can result in subsurface flow, ranging from 13% contribution to total discharge for a large event of high intensity rainfall or high discharge to 80% of total discharge for a small event with low intensity rainfall or low discharge. During large events the fraction of subsurface stormflow in the discharge is suppressed by the large amount of surface runoff. Copyright © 2008 John Wiley & Sons, Ltd.