Solute transport by surface runoff from low‐angle slopes: theory and application

The removal of chemicals in solution by overland flow from agricultural land has the potential to be a significant source of chemical loss where chemicals are applied to the soil surface, as in zero tillage and surface‐mulched farming systems. Currently, we lack detailed understanding of the transfer mechanism between the soil solution and overland flow, particularly under field conditions. A model of solute transfer from soil solution to overland flow was developed. The model is based on the hypothesis that a solute is initially distributed uniformly throughout the soil pore space in a thin layer at the soil surface. A fundamental assumption of the model is that at the time runoff commences, any solute at the soil surface that could be transported into the soil with the infiltrating water will already have been convected away from the area of potential exchange. Solute remaining at the soil surface is therefore not subject to further infiltration and may be approximated as a layer of tracer on a plane impermeable surface. The model fitted experimental data very well in all but one trial. The model in its present form focuses on the exchange of solute between the soil solution and surface water after the commencement of runoff. Future model development requires the relationship between the mass transfer parameters of the model and the time to runoff to be defined. This would enable the model to be used for extrapolation beyond the specific experimental results of this study. The close agreement between experimental results and model simulations shows that the simple transfer equation proposed in this study has promise for estimating solute loss to surface runoff. Copyright © 2000 John Wiley & Sons, Ltd.     

The generation and redistribution of overland flow on a massive oxic soil in a eucalypt woodland within the semi-arid tropics of North Australia

There is a dearth of knowledge on the runoff processes of eucalypt woodland communities in the semi-arid tropics of Australia. The work was undertaken on a 100 m transect of a 0·8 degree hillslope typical of the ‘smooth plainlands’ of central-north Queensland. This paper introduces a new experimental design for measuring overland flow in such areas by way of a cascade system of unbounded runoff plots which allow the inputs and outputs between troughs to be calculated. Most storms generate overland flow. Time to overland flow ranges between 1 and 18 min where rain intensities are above 10mm hr⁻¹ and when the average detention storage of 3·6 mm is exceeded. The bare soil surfaces within the scattered grass understory control the runoff generation process through the temporal variability of field saturated hydraulic conductivity. The study demonstrated that overland flow is mainly redistributed over the freely-draining oxic soil. Some areas export more overland flow than they gain from upslope (runoff), others gain more overland flow than they export (runon). Over the study period only 2 per cent of total rain is transferred out of this 100 m transect as overland flow due to the short duration of storms, the relatively high soil permeability, and the low slope angle. The remainder adds to the large soil water store or deep drainage. The variability of runoff–runon over these ‘smooth plainlands’ highlights how results from bounded plots would be misleading in such areas.     

Multiple modes of storm runoff generation in a North Carolina coastal plain watershed

The results of field measurements conducted in a small (19·37 ha) agricultural watershed on the North Carolina coastal plain during the summer of 1996 are presented. The objective of the study was to develop a more complete understanding of basin response in the region with respect to stormflow generation and, in particular, to identify the processes that determine storm runoff and the conditions under which such processes occur. Twenty‐four storm events were monitored, including two tropical storm systems and two hurricanes. The data demonstrate considerable spatial and temporal heterogeneity in runoff generation within the watershed. Surface flowpaths, in the form of Hortonian overland flow and saturation overland flow, were found to be the dominant runoff processes during the storm events measured. The hillslope flowpaths had the same response time as the basin streamflow, but significantly shorter time of rise and lag times. The importance of Hortonian flow in a basin with sandy, permeable soils, as well as the rapid stormflow response in a low‐relief area with a humid climate, was contrary to expectations. This, coupled with the contingency of runoff response, suggests that it may be difficult to generalize about runoff generation mechanisms in broad terms, and that a synoptic approach may be more appropriate. Copyright © 2006 John Wiley & Sons, Ltd.     

Assessing the impacts of land cover and climate on runoff and sediment yield of a river basin

ABSTRACT

This study investigated the impacts of changes in land cover and climate on runoff and sediment yield in a river basin in India. Land Change Modeler was used to derive the future land cover and its changes using the Sankey diagram approach. The future climatic parameters were derived from five general circulation models for two emission scenarios with representative concentration pathways (RCPs) 4.5 and 8.5. The land cover and climate change impacts on runoff and sediment yield were estimated using SWAT model. The results show important changes in land cover and indicate that urban and agricultural areas strongly influence the runoff and sediment yield. Among the land cover and climate change impacts, climate has more predominant (70%–95%) impact. Runoff and sediment yield are likely to decrease in both RCP scenarios in the future period. The impacts of land cover changes are more prominent on sediment yield than runoff.     

Annual distribution of river runoff with estimated contribution of winter low-water season

Proportions of the overland, delayed, and groundwater runoff are given for different phases of the hydrological regimes of rivers under different landscape–climate conditions in Russia. The concept of runoff hydrological regime is in agreement with the formation features of the annual and seasonal runoff on small watersheds under current climate conditions. The increase in winter runoff is due to the impulse discharge regime of soil water and groundwater, which receive additional recharge during frequent thaws and at weak freezing of the aeration zone. Variations in the values of runoff, total precipitation, and mean temperature were shown to be both synchronous and periodic.     

MEASURING AND MODELLING OF SOIL WATER DYNAMICS AND RUNOFF GENERATION IN AN AGRICULTURAL LOESSIAL HILLSLOPE

Surface runoff may be generated when the rainfall intensity exceeds the infiltration capacity, or when the soil profile is saturated with water. Indications exist that both types of overland flow may occur in hilly agricultural loess regions. Here, for a loessial hillslope under maize in the southern part of The Netherlands, it was shown, with pressure head and runoff measurements, that Hortonian overland flow occurs during typical summer rain events. Surface runoff was initiated after saturation of the top 5–10 cm of the soil. Deeper in the soil, unsaturated conditions prevailed while runoff took place. Peak runoff discharges at the outlet of the subcatchment occurred a few minutes after peak rainfall intensities were measured. It appeared that SWMS_2D, a two-dimensional water flow model, was capable in simulating observed pressure head changes and runoff. Simulated potential runoff for the transect studied was higher by a magnitude of three than the measured areal average. This indicates effects of surface ponding, and the probable location of this particular transect in a region with high runoff production.     

Simulating the effects of spatial configurations of agricultural ditch drainage networks on surface runoff from agricultural catchments

The study of runoff is a crucial issue because it is closely related to flooding, water quality and erosion. In cultivated catchments, agricultural ditch drainage networks are known to influence runoff. As anthropogenic elements, agricultural ditch drainage networks can therefore be altered to better manage surface runoff in cultivated catchments. However, the relationship between the spatial configuration, i.e. the density and the topology, of agricultural ditch drainage networks and surface runoff in cultivated catchments is not understood. We studied this relationship by using a random network simulator that was coupled to a distributed hydrological model. The simulations explored a large variety of spatial configurations corresponding to a thousand stochastic agricultural ditch drainage networks on a 6.4 km² Mediterranean cultivated catchment. Next, several distributed hydrological functions were used to compute water flow paths and runoff for each simulation. The results showed that (i) denser networks increased the drained volume and the peak discharge and decreased hillslopes runoff, (ii) greater network density did not affect the surface runoff any further above a given network density, (iii) the correlation between network density and runoff was weaker for small subcatchments (< 2 km²) where the variability in the drained area that resulted from changes in agricultural ditch drainage networks increased the variability of runoff and (iv) the actual agricultural ditch drainage network appeared to be well optimized for managing runoff as compared with the simulated networks. Finally, our results highlighted the role of agricultural ditch drainage networks in intercepting and decreasing overland flow on hillslopes and increasing runoff in drainage networks. Copyright © 2011 John Wiley & Sons, Ltd.     

Hydrological simulation of a conservation bench terrace system in a subhumid climate

Abstract

A finite element model to simulate runoff and soil erosion from agricultural lands has been developed. The sequential solutions of the governing differential equations were found: Richards' equation with a sink term for infiltration and soil water dynamics under cropped conditions; St Venant equation with kinematic wave approximation for overland and channel flow; and sediment continuity equation, for soil erosion. The model developed earlier has been improved to simulate erosion/deposition in impoundments and predicted and observed soil loss values were in reasonably good agreement when the model was tested for a conservation bench terrace (CBT) system. The finite element model was extensively applied to study the hydrological behaviour of a CBT system vis-à-vis the conventional system of sloping borders. The model estimates runoff and soil loss reasonably well, under varying conditions of rainfall and at different crop growth stages. The probable reasons for discrepancies between observation and simulation are reported and discussed. Sensitivity analysis was carried out to study the effect of various hydrological, soil and topographical parameters, such as ratio of contributing to receiving areas, weir length, depth of impoundment, slope of contributing area, etc. on the flow behaviour in a CBT system.     

RUNOFF GENERATION AND SOIL EROSION IN SMALL AGRICULTURAL CATCHMENTS WITH LOESS-DERIVED SOILS

Over a two-year period, rainfall, runoff and sediment output were measured in six small agricultural catchments (3–10 ha) in south Limburg (The Netherlands). These measurements were needed for validation of an erosion model for South Limburg (LISEM). In this paper, results of the measurements are presented and processes that determine surface runoff and sediment yield during winter and summer rainfall are identified. Before the start of the measurement programme, surface slaking and crust formation on the erodible loess soils were thought to be the main cause of overland flow and soil erosion in South Limburg. This was the starting point for soil conservation measures in the area. The measurement results discussed in this paper show that in some catchments much runoff occurred in winter and that soil moisture storage capacity may be just as important for runoff generation as infiltration capacity. Therefore, when modelling soil erosion and optimizing erosion control measures for South Limburg, runoff generation through Hortonian as well as through saturation overland flow must be considered.     

Modelling of water quality processes in lakes and reservoirs. II

ABSTRACT

The snowmelt runoff process from small basins is discussed. A differentiation is made between overland flow in the snowpack and groundwater flow induced by infiltrating meltwater. The effect of variations of the snowmelt intensity on streamflow is studied. It is shown that the runoff is high from the first day of snowmelt runoff if the streamflow is caused by overland flow, and that there are pronounced peaks every day, which almost correspond with the snowmelt intensity during daytime. Streamflow originating from groundwater, on the other hand, increases continuously during the snowmelt and shows only small daily peaks in the flow. Simultaneous overland and groundwater flow are also discussed. Observed runoff hydrographs from small basins are analysed in some detail. For the open fields studied the runoff shows the typical character of overland flow. For a rather large forested area the surface runoff also constitutes an important part of the runoff, but the groundwater baseflow is considerable.     

Do not only connect: a model of infiltration‐excess overland flow based on simulation

The paper focusses on connectivity in the context of infiltration‐excess overland flow and its integrated response as slope‐base overland flow hydrographs. Overland flow is simulated on a sloping surface with some minor topographic expression and spatially differing infiltration rates. In each cell of a 128 × 128 grid, water from upslope is combined with incident rainfall to generate local overland flow, which is stochastically routed downslope, partitioning the flow between downslope neighbours. Simulations show the evolution of connectivity during simple storms. As a first approximation, total storm runoff is similar everywhere, discharge increasing proportionally with drainage area. Moderate differences in plan topography appear to have only a second‐order impact on hydrograph form and runoff amount. Total storm response is expressed as total runoff, runoff coefficient or total volume infiltrated; each plotted against total storm rainfall, and allowing variations in average gradient, overland flow roughness, infiltration rate and storm duration. A one‐parameter algebraic expression is proposed that fits simulation results for total runoff, has appropriate asymptotic behaviour and responds rationally to the variables tested. Slope length is seen to influence connectivity, expressed as a scale distance that increases with storm magnitude and can be explicitly incorporated into the expression to indicate runoff response to simple events as a function of storm size, storm duration, slope length and gradient. The model has also been applied to a 10‐year rainfall record, using both hourly and daily time steps, and the implications explored for coarser scale models. Initial trails incorporating erosion continuously update topography and suggest that successive storms produce an initial increase in erosion as rilling develops, while runoff totals are only slightly modified. Other factors not yet considered include the dynamics of soil crusting and vegetation growth. Copyright © 2014 John Wiley & Sons, Ltd.     

Runoff generation in zero-order precambrian shield catchments: The stormflow response of a heterogeneous landscape

Hydrometric and isotopic (oxygen-18) observations were used to delineate the runoff processes operating in several headwater catchments on the Precambrian Shield of Canada. The catchments comprise patches of conifer forest situated on thin soils among areas of lichen-covered granitic bedrock. Horton overland flow occurred from the lichen-bedrock areas in all precipitation events that exceeded 4–6 mm. Runoff from the forest stands occurred mainly as subsurface stormflow, but in some instances saturation overland flow was observed. The occurrence of saturation overland flow was controlled by the topography of the bedrock beneath the forest soils. The area contributing runoff and the pathway by which water was conveyed to the catchment outflow switched from the open lichen-bedrock areas producing overland flow on the rising limb of the storm hydrograph to the forest stands contributing subsurface stormflow on the recession limb of the hydrograph. The areal extent and position of the landscape units in the basin were important to the rate and magnitude of stormflow production. Runoff was generated from the catchments only during and immediately after snowmelt and/or rainfall events. The catchments were dry and/or frozen for about 70% of the year.     

Effect of climate change on overland flow generation: a case study in central Germany

The impact of global climate change on runoff components, especially on the type of overland flow, is of utmost significance. High‐resolution temporal rainfall plays an important role in determining the hydrological response of quick runoff components. However, hydrological climate change scenario analyses with high temporal resolution are rare. This study investigates the impact of climate change on discharge peak events generated by rainfall, snowmelt, and soil‐frost induced runoff using high‐resolution hydrological modelling. The study area is Schäfertal catchment (1.44 km²) in the lower Harz Mountains in central Germany. The WaSiM‐ETH hydrological model is used to investigate the rainfall response of runoff components under near future (2021–2050) and far‐distant future (2071–2100) climatic conditions. Disaggregated daily climate variables of WETTREG2010 SRES scenario A1B are used on a temporal resolution of 10 min. Hydrological model parameter optimization and uncertainty analysis was conducted using the Differential Evolution Adaptive Metropolis (DREAM_(ZS)) uncertainty tool. The scenario results show that total runoff and interflow will increase by 3.8% and 3.5% in the near future and decrease by 32.85% and 31% in the far‐distant future compared to the baseline scenario. In contrast, overland flow and the number and size of peak runoff will decrease moderately for the near future and drastically for the far‐distant future compared to the baseline scenario. We found the strongest decrease for soil‐frost induced discharge peaks at 79.6% in the near future and at 98.2% in the far‐distant future scenario. It can be concluded that high‐resolution hydrological modelling can provide detailed predictions of future hydrological regimes and discharge peak events of the catchment. Copyright © 2014 John Wiley & Sons, Ltd.     

Understanding plot‐scale hydrology of Lesser Himalayan watershed—A field study and HYDRUS‐2D modelling approach

Soil moisture dynamics have a significant effect on overland flow generation. Catchment aspect is one of the major controlling factors of overland flow and soil moisture behaviour. A few experimental studies have been carried out in the uneven topography of the Himalayas. This study presents plot‐scale experiments using portable rainfall simulator at an altitude of 1,230 m above mean sea level and modelling of overland flow using observed datasets. Two plots were selected in 2 different aspects of Aglar watershed of Lesser Himalaya; the agro‐forested (AF) plot was positioned at the north aspect whereas the degraded (DE) plot was located at the south aspect of the hillslope. HS flumes and rain gauges were installed to measure the runoff at the outlet of the plot and the rainfall depth during rainfall simulation experiments. Moreover, 10 soil moisture sensors were installed at upslope and downslope locations of both the plots at 5, 15, 25, 35, and 45 cm depth from ground level to capture the soil moisture dynamics. The tests were conducted at intensities of 79.8 and 75 mm/hr in AF plot and 82.2 and 72 mm/hr in the DE plot during Test 1 and Test 2, respectively. The observed data indicate the presence of reinfiltration process only in the AF plot. The high water holding capacity and the presence of reinfiltration process results in less runoff volume in the AF plot compared with the DE plot. The Hortonian overland flow mechanism was found to be the dominant overland flow mechanism as only a few layers of top soil get saturated during all of the rainfall–runoff experiments. The runoff, rainfall, and soil moisture data were subsequently used to calibrate the parameters of HYDRUS‐2D overland flow module to simulate the runoff hydrograph and soil moisture. The components of hydrograph were evaluated in terms of peak discharge, runoff volume and time of concentration, the results were found to be within the satisfactory range. The goodness of fit of simulated hydrographs were more than 0.85 and 0.95 for AF and DE plot, respectively. The model produced satisfactory simulation results of soil moisture for all of the rainfall–runoff experiments. The HYDRUS‐2D overland flow module was found promising to simulate the runoff hydrograph and soil moisture in plot‐scale research.     

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.     

Response of snow cover and runoff to climate change in high Alpine catchments of Eastern Switzerland

A model study on the impact of climate change on snow cover and runoff has been conducted for the Swiss Canton of Graubünden. The model Alpine3D has been forced with the data from 35 Automatic Weather Stations in order to investigate snow and runoff dynamics for the current climate. The data set has then been modified to reflect climate change as predicted for the 2021–2050 and 2070–2095 periods from an ensemble of regional climate models.The predicted changes in snow cover will be moderate for 2021–2050 and become drastic in the second half of the century. Towards the end of the century the snow cover changes will roughly be equivalent to an elevation shift of 800 m. Seasonal snow water equivalents will decrease by one to two thirds and snow seasons will be shortened by five to nine weeks in 2095.Small, higher elevation catchments will show more winter runoff, earlier spring melt peaks and reduced summer runoff. Where glacierized areas exist, the transitional increase in glacier melt will initially offset losses from snow melt. Larger catchments, which reach lower elevations will show much smaller changes since they are already dominated by summer precipitation.     

Hydrological response types for gypsiferous soils in a semi‐arid region during nine years of continuous record

A 9‐year rainfall, sediment yield and runoff generation record from four experimental plots has been studied. Plots are located in the central Ebro Basin over smooth hillslope developed over gypsum and marl Miocene deposits. The hydrological response of these areas is a function of soil properties, final infiltration capacity and permeability of soils and rainfall characteristics, such as intensity and amount. Results show that there are two types of hydrological response in these areas. First is Hortonian like, which takes place during wet periods and it is responsible for the main part of total sediment yield. Second is like saturation excess overland flow, and it appears after long period. The presence of either type is controlled by a double threshold, starting from when runoff is significant. Copyright © 2002 John Wiley & Sons, Ltd.     

THE RELATIVE CONTRIBUTION OF SOIL TILLAGE AND OVERLAND FLOW EROSION TO SOIL REDISTRIBUTION ON AGRICULTURAL LAND

This study uses evidence for the long-term (35 years) pattern of soil redistribution within two agricultural fields in the UK to identify the relative importance of tillage and overland flow erosion. Spatially distributed long-term total soil redistribution data for the fields (Dalicott Farm and Rufford Forest Farm) were obtained using the caesium-137 (¹³⁷Cs) technique. These data were compared with predicted patterns of soil redistribution. Recent studies have demonstrated that the redistribution of soil by tillage may be described as a diffusive process. A two-component model was, therefore, developed which accounts for soil redistribution by both overland flow and diffusive processes. Comparison of the predicted patterns of overland flow erosion alone with the observed (¹³⁷Cs-derived) data indicated a poor agreement (r² = 0.17 and 0.11). In contrast, a good agreement exists between the predicted pattern of diffusive redistribution and the observed data (r² = 0.43 and 0.41). These results give a clear indication that diffusive processes are dominant in soil redistribution within these fields. Possible diffusive processes include splash erosion, soil creep and tillage. However, the magnitude of the diffusion coefficients for the optimum predicted pattern (c. 350–400 kg m⁻¹ a⁻¹) demonstrates that tillage is the only process capable of explaining the very significant soil redistribution which is indicated by the ¹³⁷Cs data. Consideration is given to the implications of these results for both soil erosion prediction and landscape interpretation.     

A kinematic‐wave‐based distributed watershed model using FEM,GIS and remotely sensed data

For the appropriate management of water resources in a watershed, it is essential to calculate the time distribution of runoff for the given rainfall event. In this paper, a kinematic‐wave‐based distributed watershed model using finite element method (FEM), geographical information systems (GIS) and remote‐sensing‐based approach is presented for the runoff simulation of small watersheds. The kinematic wave equations are solved using FEM for overland and channel flow to generate runoff at the outlet of the watershed concerned. The interception loss is calculated by an empirical model based on leaf area index (LAI). The Green‐Ampt Mein Larson (GAML) model is used for the estimation of infiltration. Remotely sensed data has been used to extract land use (LU)/land cover (LC). GIS have been used to prepare finite element grid and input files such as Manning's roughness and slope. The developed overland flow model has been checked with an analytical solution for a hypothetical watershed. The model has been applied to a gauged watershed and an ungauged watershed. From the results, it is seen that the model is able to simulate the hydrographs reasonably well. A sensitivity analysis of the model is carried out with the calibrated infiltration parameters, overland flow Manning's roughness, channel flow Manning's roughness, time step and grid size. The present model is useful in predicting the hydrograph in small, ungauged watersheds. Copyright © 2007 John Wiley & Sons, Ltd.