Linking frequency of rainstorms,runoff generation and sediment transport across hyperarid talus-pediment slopes

Documenting hillslope response to hydroclimatic forcing is crucial to our understanding of landscape evolution. The evolution of talus-pediment sequences (talus flatirons) in arid areas was often linked to climatic cycles, although the physical processes that may account for such a link remain obscure. Our approach is to integrate field measurements, remote sensing of rainfall and modeling to link between storm frequency, runoff, erosion and sediment transport. We present a quantitative hydrometeorological analysis of rainstorms, their geomorphic impact and their potential role in the evolution of hyperarid talus-pediment slopes in the Negev desert, Israel. Rainstorm properties were defined based on intensity–duration–frequency curves and using a rainfall simulator, artificial rainstorms were executed in the field. Then, the obtained measured experimental results were up-scaled to the entire slope length using a fully distributed hydrological model. In addition, natural storms and their hydro-geomorphic impacts were monitored using X-band radar and time-lapse cameras. These integrated analyses constrain the rainfall threshold for local runoff generation at rain intensity of 14 to 22 mm h^-1 for a duration of five minutes and provide a high-resolution characterization of small-scale runoff-generating rain cells. The current frequency of such runoff-producing rainstorms is ~1–3 per year. However, extending this local value into the full extent of hillslope runoff indicates that it occurs only under rainstorms with ≥ 100-years return interval, or 1% annual exceedance probability. Sheetwash efficiency rises with downslope distance; beyond a threshold distance of ~100 m, runoff during rainstorms with such annual exceedance probability are capable of transporting surface clasts. The erosion efficiency of these discrete rare events highlights their potential importance in shaping the landscape of arid regions. Our results support the hypothesis that a shift in the properties and frequency of extreme events can trigger significant geomorphic transitions in areas that remained hyperarid during the entire Quaternary. © 2020 John Wiley & Sons, Ltd.     

Suspended sediment discharge in subsurface flow from the head hollow of a small forested watershed,northern Japan

Surface flow and suspended sediment discharge from the head hollow of the Jozankei Experimental Watershed in Hokkaido, northern Japan, were measured to clarify the implications of subsurface hydrology for soil movement. Subsurface discharges during the extremely large storms of 1993 to 1994 were measured in a V-notch weir installed at a natural spring near the bottom of the head hollow, and shallow groundwater levels were observed in the wells excavated in the hollow. Sediment samples whose particle size range from 0·001 to 0·1 mm were manually and automatically collected at 15 to 60 min intervals, by use of 1 or 21 polyethylene bottles. Maximum concentration and flux of suspended sediment during the storms preceded the peak discharge of subsurface flow by several hours. Neither the changes in concentration (mg l⁻¹) nor flux (mg s⁻¹) of suspended sediment coincided with those in subsurface discharge (l s⁻¹). Furthermore, sediment concentration was poorly correlated with the rate of change in subsurface discharge (l s⁻²) during the rising limb of the hydrograph. Suspended sediment flux during the acceleratory limb, however, was closely correlated with the rate of change in subsurface discharge. The relationship between suspended sediment flux and rate of change in subsurface discharge were in inverse proportion to initial subsurface discharge before the storm runoff and they represented rare seasonal variation. Subsurface hydraulic erosion and transport of suspended sediment resulting from changes in rate of change in subsurface discharge actively occur during the acceleratory rising limb of the hydrograph. Accordingly, subsurface hydraulic erosion during the acceleratory rising limb of the hydrograph can be physically understood by analysing suspended sediment flux associated with rate of change in subsurface discharge and initial subsurface discharge. © 1997 John Wiley & Sons, Ltd.     

Climate controls on coupled processes of chemical weathering,bioturbation, and sediment transport across hillslopes

Most hillslope studies examining the interplay between climate and earth surface processes tend to be biased towards eroding parts of landscapes. This limitation makes it difficult to assess how entire upland landscapes, which are mosaics of eroding and depositional areas, evolve physio‐chemically as a function of climate. Here we combine new soil geochemical data and published ¹⁰Be‐derived soil production rates to estimate variations in chemical weathering across two eroding‐to‐depositional hillslopes spanning a climate gradient in southeastern Australia. At the warmer and wetter Nunnock River (NR) site, rates of total soil (–3 to –14 g m^‐2 yr^‐1; negative sign indicates mass loss) and saprolite (–18 to –32 g m^‐2 yr^‐1) chemical weathering are uniform across the hillslope transect. Alternatively, the drier hillslope at Frog's Hollow (FH) is characterized by contrasting weathering patterns in eroding soils (–30 to –53 g m^‐2 yr^‐1) vs. depositional soils (+91 g m^‐2 yr^‐1; positive sign indicates mass addition). This difference partly reflects mineral grain size sorting as a result of upslope bioturbation coupled with water‐driven soil erosion, as well as greater vegetative productivity in moister depositional soils. Both of these processes are magnified in the drier climate. The data reveal the importance of linking the erosion–deposition continuum in hillslope weathering studies in order to fully capture the coupled roles of biota and erosion in driving the physical and chemical evolution of hillslopes. Our findings also highlight the potential limitations of applying current weathering models to landscapes where particle‐sorting erosion processes are active. Copyright © 2018 John Wiley & Sons, Ltd.     

Soil detachment and transport on field‐ and laboratory‐scale interrill areas: erosion processes and the size‐selectivity of eroded sediment

Field‐ and laboratory‐scale rainfall simulation experiments were carried out in an investigation of the temporal variability of erosion processes on interrill areas, and the effects of such variation upon sediment size characteristics. Poorly aggregated sandy soils from the semi‐arid environment of Senegal, West Africa, were used on both a 40 m² field plot and a 0·25 m² laboratory plot; rainfall intensity for all experiments was 70 mm h^?1 with a duration of 1 to 2 hours. Time‐series measurements were made of the quantity and the size distribution of eroded material: these permitted an estimate of the changing temporal balance between the main erosion processes (splash and wash). Results from both spatial scales showed a similar temporal pattern of runoff generation and sediment concentration. For both spatial scales, the dominant erosional process was detachment by raindrops; this resulted in a dynamic evolution of the soil surface under raindrop impact, with the rapid formation of a sieving crust followed by an erosion crust. However, a clear difference was observed between the two scales regarding the size of particles detached by both splash and wash. While all measured values were lower than the mean weight diameter (MWD) value of the original soil (mean 0·32 mm), demonstrating the size‐selective nature of wash and splash processes, the MWD values of washed and splashed particles at the field scale ranged from 0·08 to 0·16 mm and from 0·12 to 0·30 mm respectively, whereas the MWD values of washed and splashed particles at the laboratory scale ranged from 0·13 to 0·29 mm and from 0·21 to 0·32 mm respectively. Thus only at the field scale were the soil particles detached by splash notably coarser than those transported by wash. This suggests a transport‐limited erosion process at the field scale. Differences were also observed between the dynamics of the soil loss by wash at the two scales, since results showed wider scatter in the field compared to the laboratory experiments. This scatter is probably related to the change in soil surface characteristics due to the size‐selectivity of the erosion processes at this spatial scale. Copyright © 2006 John Wiley & Sons, Ltd.     

Effects of hillslope topography on hydrological responses in a weathered granite mountain,Japan: comparison of the runoff response between the valley‐head and the side slope

To evaluate the effects of hillslope topography on storm runoff in a weathered granite mountain, discharge rate, soil pore water pressures, and water chemistry were observed on two types of hillslope: a valley‐head (a concave hillslope) and a side slope (a planar hillslope). Hydrological responses on the valley‐head and side slope reflected their respective topographic characteristics and varied with the rainfall magnitude. During small rainfall events (<35 mm), runoff from the side slope occurred rapidly relative to the valley‐head. The valley‐head showed little response in storm runoff. As rainfall amounts increased (35–60 mm), the valley‐head yielded a higher flow relative to the side slope. For large rainfall events (>60 mm), runoff from both hillslopes increased with rainfall, although that from the valley‐head was larger than that from the side slope. The differences in the runoff responses were caused by differences in the roles of lower‐slope soils and the convergence of the hillslope. During small rainfall events, the side slope could store little water; in contrast, all rainwater could be stored in the soils at the valley‐head hollow. As the amount of rainfall increased, the subsurface saturated area of the valley‐head extended from the bottom to the upper portion of the slope, with the contributions of transient groundwater via lateral preferential flowpaths due to the high concentration of subsurface water. Conversely, saturated subsurface flow did not contribute to runoff responses, and the subsurface saturated area at the side slope did not extend to the upper slope for the same storm size. During large rainfall events, expansion of the subsurface saturated area was observed in both hillslopes. Thus, differences in the concentration of subsurface water, reflecting hillslope topography, may create differences in the extension of the subsurface saturated area, as well as variability in runoff responses. Copyright © 2007 John Wiley & Sons, Ltd.     

On the Panday and Huyakorn surface–subsurface hydrology test case: analysis of internal flow dynamics

A number of previous studies using models of integrated surface‐subsurface hydrology have adopted the Panday and Huyakorn (P&H) tilted V‐catchment test case (Panday S, Huyakorn PS. 2004. A fully coupled spatially distributed model for evaluating surface/subsurface flow. Advances in Water Resources 27: 361–382) to show inter‐code comparability. The P&H test case is used to evaluate models that simulate a broad range of hydrological processes, and yet only the catchment outflow hydrograph has been presented as verification of the consistency between codes. Therefore, a more comprehensive evaluation of the surface‐subsurface hydrology of the P&H case is needed. This study explores the internal catchment functioning of the P&H case, using the popular catchment simulator MODHMS. The processes leading to streamflow generation in the model are illustrated, including separation of overland flow (OLF) and groundwater discharge to the stream. The results identify non‐physical flow processes due to the problem set‐up, and modifications to the P&H case are suggested that include changes to stream roughness and incision of the stream channel to overcome these shortcomings. A modified P&H case produced more plausible transfers between OLF and the stream, and an increased groundwater discharge to the stream (6·5% of streamflow in the modified case compared to 0·5% in the original case). Despite changes to internal flow processes, near‐identical outflow hydrographs were obtained, showing the importance of considering and comparing internal flow processes when using surface‐subsurface hydrology test cases to evaluate integrated hydrological simulators. Copyright © 2011 John Wiley & Sons, Ltd.     

The impact of slope length on the discharge of sediment by rain impact induced saltation and suspension

Simulations using a mechanistic model of raindrop driven erosion in rain‐impacted flow were performed with particles travelling by suspension, raindrop induced saltation and flow driven saltation. Results generated by both a high intensity storm, and a less intense one, indicate that, because of the effect of flow depth on the delivery of raindrop energy to the bed, there is a decline in sediment concentration, and hence soil loss per unit area, with slope length when particles are transported by raindrop induced saltation. However, that decline is reversed when the critical velocities that lead to flow driven saltation are episodically exceeded during an event. The simulations were performed on smooth surfaces and a single drop size but the general relationships are likely to apply for rain made up of a wide range of drop size. Although runoff is not always produced uniformly, as a general rule, flow velocities increase with slope length so that, typically, the distance particles travel before being discharged during an event increase with slope length. The effect of slope length on soil loss per unit area is often considered to vary with slope length to a power greater than zero and less that 1·0. The simulations show that effect of slope length on sediment discharge is highly dependent on the variations in runoff response resulting from variations in rainfall duration‐intensity‐infiltration conditions rather than plot length per se. Consequently, predicting soil loss per unit area using slope length with positive powers close to zero when sheet erosion occurs may not be as effective as commonly expected. Erosion by rain‐impacted flow is a complex process and that complexity needs to be considered when analysing the results of experiments associated with rain‐impacted flow under both natural and artificial conditions. Copyright © 2009 John Wiley & Sons, Ltd.     

The role of overland flow and subsurface flow on the spatial distribution of soil moisture in the topsoil

Many investigations show relationships between topographical factors and the spatial distribution of soil moisture in catchments. However, few quantitative analyses have been carried out to elucidate the role of different hydrological processes in the spatial distribution of topsoil moisture in catchments. A spatially distributed rainfall—runoff model was used to investigate contributions of subsurface matric flow, macropore flow and surface runoff to the spatial distribution of soil moisture in a cultivated catchment. The model results show that lateral subsurface flow in the soil matrix or in macropores has a minor effect on the spatial distribution of soil moisture. Only when a perched groundwater table is maintained long enough, which is only possible if the subsurface is completely impermeable, may a spatial distribution in moisture content occur along the slope. Surface runoff, producing accumulations of soil moisture in flat flow paths of agricultural origin (field boundaries), was demonstrated to cause significant spatial variations in soil moisture within a short period after rainfall (<2 days). When significant amounts of surface runoff are produced, wetter moisture conditions will be generated at locations with larger upstream contributing areas. Copyright © 2001 John Wiley & Sons, Ltd.     

Assessing the impact of the hydraulic properties of a crusted soil on overland flow modelling at the field scale

Soil surface crusts are widely reported to favour Hortonian runoff, but are not explicitly represented in most rainfall‐runoff models. The aim of this paper is to assess the impact of soil surface crusts on infiltration and runoff modelling at two spatial scales, i.e. the local scale and the plot scale. At the local scale, two separate single ring infiltration experiments are undertaken. The first is performed on the undisturbed soil, whereas the second is done after removal of the soil surface crust. The HYDRUS 2D two‐dimensional vertical infiltration model is then used in an inverse modelling approach, first to estimate the soil hydraulic properties of the crust and the subsoil, and then the effective hydraulic properties of the soil represented as a single uniform layer. The results show that the crust hydraulic conductivity is 10 times lower than that of the subsoil, thus illustrating the limiting role the crust has on infiltration. Moving up to the plot scale, a rainfall‐runoff model coupling the Richards equation to a transfer function is used to simulate Hortonian overland flow hydrographs. The previously calculated hydraulic properties are used, and a comparison is undertaken between a single‐layer and a double‐layer representation of the crusted soil. The results of the rainfall‐runoff model show that the soil hydraulic properties calculated at the local scale give acceptable results when used to model runoff at the plot scale directly, without any numerical calibration. Also, at the plot scale, no clear improvement of the results can be seen when using a double‐layer representation of the soil in comparison with a single homogeneous layer. This is due to the hydrological characteristics of Hortonian runoff, which is triggered by a rainfall intensity exceeding the saturated hydraulic conductivity of the soil surface. Consequently, the rainfall‐runoff model is more sensitive to rainfall than to the subsoil's hydrodynamic properties. Therefore, the use of a double‐layer soil model to represent runoff on a crusted soil does not seem necessary, as the increase of precision in the soil discretization is not justified by a better performance of the model. Copyright © 2005 John Wiley & Sons, Ltd.     

The importance of surface controls on overland flow connectivity in semi‐arid environments: results from a numerical experimental approach

In semi‐arid environments, the characteristics of the land surface determine how rainfall is transformed into surface runoff and influences how this runoff moves from the hillslopes into river channels. Whether or not water reaches the river channel is determined by the hydrological connectivity. This paper uses a numerical experiment‐based approach to systematically assess the effects of slope length, gradient, flow path convergence, infiltration rates and vegetation patterns on the generation and connectivity of runoff. The experiments were performed with the Connectivity of Runoff Model, 2D version distributed, physically based, hydrological model. The experiments presented are set within a semi‐arid environment, characteristic of south‐eastern Spain, which is subject to low frequency high rainfall intensity storm events. As a result, the dominant hydrological processes are infiltration excess runoff generation and surface flow dynamics. The results from the modelling experiments demonstrate that three surface factors are important in determining the form of the discharge hydrograph: the slope length, the slope gradient and the infiltration characteristics at the hillslope‐channel connection. These factors are all related to the time required for generated runoff to reach an efficient flow channel, because once in this channel, the transmission losses significantly decrease. Because these factors are distributed across the landscape, they have a fundamental role in controlling the landscape hydrological response to storm events. Copyright © 2013 John Wiley & Sons, Ltd.     

Response of soil erosion and sediment sorting to the transport mechanism on a steep rocky slope

A deeper understanding of the sediment characteristics associated with rock fragment content can improve our knowledge of the erosional processes and transport mechanisms of sediments on steep rocky slopes. This research used simulated rainfall experiments lasting for 1 h at a rate of 90 mm h⁻¹ and employed 5 × 1 × 0.4 m parallel troughs filled with purple soils with different rock fragment volumetric contents (0, 5, 10, 20, 30 and 40%) on a 15° slope gradient. For each simulated event, runoff and sediment were sampled at 1- and 3-min intervals, respectively, to study, in detail, the temporal changes in the size distributions of the eroded sediments. The results show that sediment concentrations, soil erosion rates and soil loss ratios significantly decreased as rock fragment content increased for rock fragment contents from 0 to 40% in purple soils. During the transportation process, clay particles often formed aggregates and were then transported as larger particles. Silt particles were more likely to be transported as primary particles with a low degree of sediment aggregation. Sand-sized particles, which constituted a greater proportion of the original soil than the eroded sediments, were formed from other fine particles and transported as aggregates rather than as primary particles. Suspension-saltation, which mainly transports fine particles of 0.02–0.05 mm and coarse particles larger than 0.5 mm in size, was the most important transport mechanism on steep rocky slopes. The results of this study can help to explain the inherent laws of erosional processes on steep rocky slopes and can provide a foundation for improving physical models of soil erosion. © 2019 John Wiley & Sons, Ltd.     

The influence of landscape connectivity and landslide dynamics upon channel adjustments and sediment flux in the Liwu Basin,Taiwan

Catchment‐scale analyzes of spatial and temporal variability in landscape connectivity are critical considerations in appraisals of landscape evolution and disaster mitigation in tectonically active mountain belts such as Taiwan. This study uses historical aerial photographs, flow discharge and seismic data to analyze landslide changes and channel adjustments over a 30 year period in the Liwu Basin. Recurrent earthquakes and typhoon events trigger frequent landslide activity, channel adjustment and sediment reworking in this system. Spatial variability in magnitude–frequency relations of hillslope‐valley floor (lateral) and upstream–downstream (longitudinal) connectivity during the study period are shown to reflect annual reworking in source and accumulation zones, while partly‐confined valleys in the mid‐catchment area trap sediment behind landslide‐induced dams that are formed and breached on an approximately decadal basis. This promotes partial longitudinal connectivity in these areas. Landscape responses to disturbance events were especially pronounced following combinations of seismic and typhoon events prior to the 1998 and 2005 images. Although single high magnitude events and series of moderate events affect patterns of landscape connectivity in the Liwu Basin, residence times for sediment storage are very short in this highly‐connected river system, where confined valley settings extend virtually to the coast. Copyright © 2014 John Wiley & Sons, Ltd.     

Rill characteristics and sediment transport as a function of slope length during a storm event on loess soil

On the basis of detailed rill surveys carried out on bare plots of different lengths at slopes of 12 per cent, basic rill parameters were derived. Rill width and maximum depth increased with plot length, whereas rill amount and cross‐sectional area, expressed per unit length, remained similar. On smaller plots, all rills were connected in a continuous transport system reaching the plot outlet, whilst on larger plots (10 and 20 m long) part of the rills ended with a deposition areas inside the plots. Amounts of erosion, calculated from rill volume and soil bulk density, were compared with soil loss measured at the plot outlets. On plots 10 and 20 m long, erosion estimated from volume of all rills was larger than measured soil loss. The latter was larger than erosion estimated from volume of contributing rills. To identify contributing soil loss area on these plots, two methods were applied: (i) ratio of total soil loss to maximum soil loss per unit area, and (ii) partition of plot area according to the ratio of contributing to total rill volume. Both methods resulted in similar areas of 21·8–23·5 m² for the plot 10 m long and 31·2 m² for the plot 20 m long. Identification of contributing areas enabled rill (5·9 kg m^?2) and interrill (2·6 kg m^?2) erosion rate to be calculated, the latter being very close to the value predicted from the Universal Soil Loss Equation. Although rill and interrill rates seemed to be similar on all plots, their ratio increased slightly with plot length. Application of this ratio to compute slope length factor of the Revised Universal Soil Loss Equation resulted in similar values to those predicted with the model. The achieved balance of soil loss suggested that all the sediment measured at the plot outlet originated from contributing rills and associated contributing rill areas. The results confirmed the utility of different plot lengths as a research tool for analysing the dynamic response of soil to rainfall–runoff. Copyright © 2005 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.     

The role of run-on for overland flow and the characteristics of runoff generation in the Loess Plateau,China

Abstract

The Loess Plateau in China is overlain by deep and loose soil. As in other semi-arid regions, convective precipitation produces storms, typically of short duration, relatively high intensity and limited areal extent. Infiltration excess (Hortonian mechanism) of precipitation is conventionally assumed to be more prominent than saturation excess (Dunne mechanism) for storm runoff generation. This assumption is true at a point during the storm. However, the runoff generation mechanism is altered when the runoff is conditioned by a lateral redistribution movement of water, i.e. run-on, as the spatial scale increases. In the Loess Plateau, the effects of run-on may be significant, because of the deep and loose surface soil layer. In this study, the role of run-on for overland flow in the Upper Wei River basin, located in the Loess Plateau, is evaluated by means of a simple numerical model at the hillslope scale. The results show that almost all the Hortonian overland flow infiltrates into the soil along the flat hillslope and dry gully before it reaches the river channel. Most of the runoff is generated from the saturated soil near the river channel and from the subsurface. The run-on process takes much longer than the infiltration, facilitating rainfall–runoff modelling at a daily time step. A hydrological model is employed to investigate the characteristics of runoff generation in the Upper Wei River basin. The analysis shows that the subsurface flow contribution to total streamflow is more than 53% from October to March, while the overland flow contribution exceeds 72% from April to September.

Editor D. Koutsoyiannis; Associate editor Dawen Yang

Citation Liu, D.F., Tian, F.Q., Hu, H.C., and Hu, H.P., 2012. The role of run-on for overland flow and the characteristics of runoff generation in the Loess Plateau, China. Hydrological Sciences Journal, 57 (6), 1107–1117.     

The effect of antecedent moisture conditions on sediment and phosphorus loss during overland flow: Mahantango Creek catchment,Pennsylvania, USA

The loss of P in overland flow from most cultivated soils is controlled by erosion, and in‐turn soil moisture. We evaluated the effect of soil moisture on erosion and P transport in overland flow by applying rainfall (7 cm h^?1) to packed soil boxes (1 m long and 0·15 m wide) and field plots (1 and 10 m long by 1 m wide) of silt loams in a central Pennsylvania (USA) catchment. Flow from packed soil boxes took longer to initiate as antecedent soil moisture decreased from field capacity (2 min) to air dried (8 to 9 min). Even in the more complex field plots (i.e. soil heterogeneity and topography), the wetter site (1 by 10 m plot; 70% field capacity) produced flow more quickly (3 min) and in greater volume (439 L) than the drier site (1 by 10 m plot; 40% field capacity, 15 min, and 214 L, respectively). However, less suspended sediment was transported from wetter soil boxes (1·6 to 2·5 g L^?1) and field plots (0·9 g L^?1) than drier boxes (2·9 to 4·2 g L^?1) and plots (1·2 g L^?1). Differences are attributed to their potential for soil aggregate breakdown, slaking and dispersion, which contribute to surface soil sealing and crusting, as dry soils are subject to rapid wetting (by rainfall). During flow, selective erosion and antecedent moisture conditions affected P transport. At field capacity, DRP and PP transport varied little during overland flow. Whereas P transport from previously dry soil decreased rapidly after the initiation of flow (6 to 1·5 mg TP L^?1), owing to the greater slaking and dispersion of P‐rich particles into flow at the beginning than end of the flow event. These results indicate that soil moisture fluctuations greatly effect erosion and P transport potential and that management to decrease the potential for loss should consider practices such as conservation tillage and cover crops, particularly on areas where high soil P and erosion coincide. Copyright © 2002 John Wiley & Sons, Ltd.     

The effect of soil type,meteorological forcing and slope gradient on the simulation of internal erosion processes at the local scale

Numerical simulation experiments of water erosion at the local scale (20 × 5 m) using a process‐based model [Plot Soil Erosion Model_2D (PSEM_2D)] were carried out to test the effects of various environmental factors (soil type, meteorological forcing and slope gradient) on the runoff and erosion response and to determine the dominant processes that control the sediment yield at various slope lengths. The selected environmental factors corresponded to conditions for which the model had been fully tested beforehand. The use of a Green and Ampt model for infiltration explained the dominant role played by rainfall intensity in the runoff response. Sediment yield at the outlet of the simulated area was correlated positively with rainfall intensity and slope gradient, but was less sensitive to soil type. The relationship between sediment yield (soil loss per unit area) and slope length was greatly influenced by all environmental factors, but there was a general tendency towards higher sediment yield when the slope was longer. Contribution of rainfall erosion to gross erosion was dominant for all surfaces with slope lengths ranging from 4 to 20 m. The highest sediment yields corresponded to cases where flow erosion was activated. An increase in slope gradient resulted in flow detachment starting upstream. Sediment exported at the outlet of the simulated area came predominantly from the zone located near the outlet. The microrelief helped in the development of a rill network that controlled both the ratio between rainfall and flow erosion and the relationship between sediment yield and slope length. Copyright © 2010 John Wiley & Sons, Ltd.     

Impact of grid resolution on the integrated and distributed response of a coupled surface–subsurface hydrological model for the des Anglais catchment,Quebec

Digital elevation models (DEMs) at different resolutions (180, 360, and 720 m) are used to examine the impact of different levels of landscape representation on the hydrological response of a 690‐km² catchment in southern Quebec. Frequency distributions of local slope, plan curvature, and drainage area are calculated for each grid size resolution. This landscape analysis reveals that DEM grid size significantly affects computed topographic attributes, which in turn explains some of the differences in the hydrological simulations. The simulations that are then carried out, using a coupled, process‐based model of surface and subsurface flow, examine the effects of grid size on both the integrated response of the catchment (discharge at the main outlet and at two internal points) and the distributed response (water table depth, surface saturation, and soil water storage). The results indicate that discharge volumes increase as the DEM is coarsened, and that coarser DEMs are also wetter overall in terms of water table depth and soil water storage. The reasons for these trends include an increase in the total drainage area of the catchment for larger DEM cell sizes, due to aggregation effects at the boundary cells of the catchment, and to a decrease in local slope and plan curvature variations, which in turn limits the capacity of the watershed to transmit water downslope and laterally. The results obtained also show that grid resolution effects are less pronounced during dry periods when soil moisture dynamics are mostly controlled by vertical fluxes of evaporation and percolation. Copyright © 2010 John Wiley & Sons, Ltd.     

Analysis of the contributions of topographic,soil, and vegetation features on the spatial distributions of surface soil moisture in a steep natural forested headwater catchment

Surface soil moisture has been extensively studied for various land uses and landforms. Although many studies have reported potential factors that control surface soil moisture over space or time, the findings have not always been consistent, indicating a need for identification of the main factors. This study focused on the static controls of topographic, soil, and vegetation features on surface soil moisture in a steep natural forested headwater catchment consisting of three hillslope units of a gully area, side slope, and valley‐head slope. Using a simple correlation analysis to investigate the effects of the static factors on surface soil moisture at depths of 0–20 cm at 470 points in 13 surveys, we addressed the characteristics of surface soil moisture and its main controlling factors. The results indicated that the mean of surface soil moisture was in the decreasing order of gully area > valley‐head slope > side slope. The relationship between the mean and standard deviation of surface soil moisture showed a convex‐upward shape in the headwater catchment, a negative curvilinear shape in the gully area, and positive curvilinear shapes at the side and valley‐head slopes. At the headwater catchment and valley‐head slope, positive contributions of soil porosity and negative contributions of slope gradient and saturated hydraulic conductivity were the main controlling factors of surface soil moisture under wetter conditions, whereas positive contributions of topographic wetness index and negative contributions of vegetation density were the main controlling factors of surface soil moisture under drier conditions. At the side slope underlain by fractured bedrocks, only saturated hydraulic conductivity and vegetation density were observed to be the controlling factors. Surface soil moisture in the gully area was mainly affected by runoff rather than were static features. Thus, using hillslope units is effective for approximately estimating the hydrological behaviours of surface moisture on a larger scale, whereas dependency between the main static factors and moisture conditions is helpful for estimating the spatial distributions of surface moisture on a smaller scale.     

Modelling the impacts of subsurface drainage on surface runoff and sediment yield in the Le Sueur Watershed,Minnesota, USA

Abstract

Tile drainage influences infiltration and surface runoff and is thus an important factor in the erosion process. Tile drainage reduces surface runoff, but questions abound on its influence on sediment transport through its dense network and into the stream network. The impact of subsurface tiling on upland erosion rates in the Le Sueur River watershed, USA, was assessed using the Water Erosion Prediction Project (WEPP) model. Six different scenarios of tile drainage with varying drainage coefficient and management type (no till and autumn mulch-till) were evaluated. The mean annual surface runoff depth, soil loss rate and sediment delivery ratio (SDR) for croplands, based on a 30-year simulation for the watershed with untiled autumn mulch-till (Scenario 1), were estimated to be 83.5 mm, 0.27 kg/m² and 86.7%, respectively; on no-till management systems (Scenario 4), the respective results were 72.3 mm, 0.06 kg/m² and 88.2%. Tile drains reduced surface runoff, soil loss and SDR estimates for Scenario 1 by, on average, 14.5, 8.1 and 7.9%, respectively; and for Scenario 4 by an estimated 31.5, 22.1 and 20.2%, respectively. The impact of tile drains on surface runoff, soil loss and SDR was greater under the no-till management system than under the autumn mulch-till management system. Comparison of WEPP outputs with those of the Soil Water Assessment Tool (SWAT) showed differences between the two methods.

Editor Z.W. Kundzewicz

Citation Maalim, F.K. and Melesse. A.M., 2013. Modelling impacts of subsurface drainage on surface runoff and sediment yield in the Le Sueur Watershed in Minnesota, USA. Hydrological Sciences Journal, 58 (3), 570–586.