A model to couple overland flow and infiltration into macroporous vadose zone

Most vegetated land surfaces contain macropores that may have a significant effect on the rate of infiltration of water under ponded conditions on the ground surface. Owing to the small-scale variations of the land topography (microtopography), only portions of the land area may get ponded during the process of overland flow. As the macropores transmit water at much higher rates than the primary soil matrix, higher macropore activation in ponded areas produces larger effective infiltration rates into the soil. Therefore, overland flow and infiltration into the macroporous vadose zone are interrelated. Representing the microtopographic variation of the land surface by a simple sine wave function, a method was developed to relate the ponding area to the average ponding depth which was determined by overland flow. A numerical model coupling overland flow and infiltration into the macroporous vadose zone was developed. Overland flow was simulated using the St. Venant equations with the inertia terms neglected. A single macropore model was used to simulate the infiltration into the macroporous vadose zone. The interaction between overland flow and the infiltration into the macroporous vadose zone was analyzed for a hypothetical watershed. The sensitivity analysis revealed that the interaction of macropore flow and overland flow is significant. For the conditions tested, the macropore flow and the overland flow were found to be more sensitive to the macroporosity and less sensitive to the microtopographic surface variation.     

Matching hydrologic response to measured effective hydraulic conductivity

The objective of this study was to test the practicability of defining hydrologic response units as combinations of soil, land use and topography for modelling infiltration at the hillslope and catchment scales. In an experimental catchment in the East African Highlands (Kwalei, Tanzania), three methods of measuring infiltration were compared for their ability to capture the spatial variability of effective hydraulic conductivity: the constant head (CH) method; the tension infiltration (TI) method; and the mini‐rainfall simulation (RS) method. The three methods yielded different probability distributions of effective hydraulic conductivity and suggested different types of hydrologic response units. Independently from these measurements, the occurrence of infiltration‐excess overland flow was monitored over an area of 6 ha by means of overland flow detectors. The observed pattern of overland flow occurrence did not match any of the patterns suggested by the infiltration measurements. Instead, clusters of spots with overland flow were practically independent from field borders. Geostatistical analysis of the overland flow confirmed the absence of spatial correlation for distances over 40 m. The RS method yielded the pattern closest to the observations, probably because the method simulated better the processes that trigger infiltration‐excess overland flow, i.e. soil sealing and infiltration through macroporosity. The RS hydrologic response unit correlated significantly with observed overland flow frequency. However, the location of clusters and ‘hot spots’ of overland flow remained largely unexplained by land use, soil and topographic variables. It is concluded that using such landscape variables to define hydrologic units may create artificial boundaries that do no correspond to physical realities, especially if the stochastic component within hydrologic units is neglected. Copyright © 2005 John Wiley & Sons, Ltd.     

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

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

Topographic controls upon soil macropore flow

Macropores are important components of soil hydrology. The spatial distribution of macropore flow as a proportion of saturated hydraulic conductivity was tested on six humid–temperate slopes using transects of tension infiltrometer measurements. Automated water table and overland flow monitoring allowed the timing of, and differentiation between, saturation‐excess overland flow and infiltration‐excess overland flow occurrence on the slopes to be determined and related to tension‐infiltrometer measurements. Two slopes were covered with blanket peat, two with stagnohumic gleys and two with brown earth soils. None of the slopes had been disturbed by agricultural activity within the last 20 years. This controlled the potential for tillage impacts on macropores. The proportion of near‐surface macropore flow to saturated hydraulic conductivity was found to vary according to slope position. The spatial patterns were not the same for all hillslopes. On the four non‐peat slopes there was a relationship between locations of overland flow occurrence and reduced macroporosity. This relationship did not exist for the peat slopes investigated because they experienced overland flow across their whole slope surfaces. Nevertheless, they still had a distinctive spatial pattern of macropore flow according to slope position. For the other soils tested, parts of slopes that were susceptible to saturation‐excess overland flow (e.g. hilltoes or flat hilltops) tended to have least macropore flow. To a lesser extent, for the parts of slopes susceptible to infiltration‐excess overland flow, the proportion of macropore flow as a component of infiltration was also smaller compared with the rest of the slope. The roles of macropore creation and macropore infilling by sheet wash are discussed, and it is noted that the combination of these may result in distinctive topographically controlled spatial patterns of macropore flow. Copyright © 2008 John Wiley & Sons, Ltd.     

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

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

The errors in surface runoff prediction by neglecting the relationship between infiltration rate and overland flow depth

Many simplifications are used in modeling surface runoff over a uniform slope. A very common simplification is to determine the infiltration rate independent of the overland flow depth and to combine it afterward with the kinematic-wave equation to determine the overland flow depth. Another simplication is to replace the spatially variable infiltration rates along the slope i(x, t) due to the water depth variations h(x,t) with an infiltration rate that is determined at a certain location along the slope. The aim of this study is to evaluate the errors induced by these simplications on predicted infiltration rates, overland flow depths, and total runoff volume. The error analysis is accomplished by comparing a simplified model with a model where the interaction between the overland flow depth and infiltration rate is counted. In this model, the infiltration rate is assumed to vary along the slope with the overland flow depth, even for homogeneous soil profiles. The kinematic-wave equation with interactive infiltration rate, calculated along the slopy by Richard's equation, are then solved by a finite difference scheme for a 100-m-long uniform slope. In the first error analysis, we study the effect of combining an ‘exact’ and ‘approximate’ one-dimensional infiltration rate with the kinematic-wave equation for three different soil surface roughness coefficients. The terms ‘exact’ and ‘approximate’ stand for the solution of Richard's equation with and without using the overland flow depth in the boundary condition, respectively. The simulations showed that higher infiltration rates and lower overland flow depths are obtained during the rising stage of the hydrograph when overland flow depth is used in the upper boundary condition of the one-dimensional Richard's equation. During the recession period, the simplified model predicts lower infiltration rates and higher overland flow depths. The absolute relative errors between the ‘exact’ and ‘approximate’ solutions are positively correlated to the overland flow depths which increase with the soil surface roughness coefficient. For this error analysis, the relative errors in surface runoff volume per unit slope width throughout the storm are much smaller than the relative errors in momentary overland flow depths and discharges due to the alternate signs of the deviations along the rising and falling stages. In the second error analysis, when the spatially variable infiltration rate along the slope i(x, t) is replaced in the kinematic-wave equation by i(t), calculated at the slope outlet, the overland flow depth is underestimated during the rising stage of the hydrograph and overestimated during the falling stage. The deviations during the rising stage are much smaller than the deviations during the falling stage, but they are of a longer duration. This occurs because the solution with i(x, t) recognizes that part of the slope becomes dry after rainfall stops, while overland flow still exists with i(t) determined at the slope outlet. As obtained for the first error analysis, the relative errors in surface runoff volume per unit slope width are also much smaller than the relative errors in momentary overland flow depths and discharges. The relation between the errors in overland flow depth and discharge to different mathematical simplifications enables to evaluate whether certain simplifications are justified or more computational efforts should be used.     

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

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

LISEM: A SINGLE-EVENT,PHYSICALLY BASED HYDROLOGICAL AND SOIL EROSION MODEL FOR DRAINAGE BASINS. II: SENSITIVITY ANALYSIS,VALIDATION AND APPLICATION

A new hydrological and soil erosion model has been developed and tested: LISEM, the Limburg soil erosion model. The model uses physically based equations to describe interception, infiltration and soil water transport, storage in surface depressions, splash and flow detachment, transport capacity and overland and channel flow. From the validation results it is clear that, although the model has several advantages over other models, the results of LISEM 1.0 are far from perfect. Based on the sensitivity analysis and field observations, the main reasons for these differences seems to be the spatial and temporal variability of the soil hydraulic conductivity and the initial pressure head at the basin scale. Another reason for the differences between measured and simulated results is our lack or understanding of the theory of hydrological and soil erosion processes.     

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

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

Assessing plot-scale impacts of land use on overland flow generation in Central Panama

Land use in Panama has changed dramatically with ongoing deforestation and conversion to cropland and cattle pastures, potentially altering the soil properties that drive the hydrological processes of infiltration and overland flow. We compared plot-scale overland flow generation between hillslopes in forested and actively cattle-grazed watersheds in Central Panama. Soil physical and hydraulic properties, soil moisture and overland flow data were measured along hillslopes of each land-use type. Soil characteristics and rainfall data were input into a simple, 1-D representative model, HYDRUS-1D, to simulate overland flow that we used to make inferences about overland flow response at forest and pasture sites. Runoff ratios (overland flow/rainfall) were generally higher at the pasture site, although no overall trends were observed between rainfall characteristics and runoff ratios across the two land uses at the plot scale. Saturated hydraulic conductivity (K_s) and bulk density were different between the forest and pasture sites (p < 10⁻⁴). Simulating overland flow in HYDRUS-1D produced more outputs similar to the overland flow recorded at the pasture site than the forest site. Results from our study indicate that, at the plot scale, Hortonian overland flow is the main driver for overland flow generation at the pasture site during storms with high-rainfall totals. We infer that the combination of a leaf litter layer and the activation of shallow preferential flow paths resulting in shallow saturation-excess overland flow are likely the main drivers for plot scale overland flow generation at the forest site. Results from this study contribute to the broader understanding of the delivery of freshwater to streams, which will become increasingly important in the tropics considering freshwater resource scarcity and changing storm intensities.     

Parametric study of a physically-based,plot-scale hillslope hydrological model through virtual experiments

Abstract

A physically-based hillslope hydrological model with shallow overland flow and rapid subsurface stormflow components was developed and calibrated using field experiments conducted on a preferential path nested hillslope in northeast India. Virtual experiments were carried out to perform sensitivity analysis of the model using the automated parameter estimation (PEST) algorithm. Different physical parameters of the model were varied to study the resulting effects on overland flow and subsurface stormflow responses from the theoretical hillslopes. It was observed that topographical shapes had significant effects on overland flow hydrographs. The slope profiles, surface storage, relief, rainfall intensity and infiltration rates primarily controlled the overland flow response of the hillslopes. Prompt subsurface stormflow responses were mainly dominated by lateral preferential flow, as soil matrix flow rates were very slow. Rainfall intensity and soil macropore structures were the most influential parameters on subsurface stormflow. The number of connected soil macropores was a more sensitive parameter than the size of macropores. In hillslopes with highly active vertical and lateral preferential pathways, saturation excess overland flow was not evident. However, saturation excess overland flow was generated if the lateral macropores were disconnected. Under such conditions, rainfall intensity, duration and preferential flow rate governed the process of saturation excess overland flow generation from hillslopes.

Editor D. Koutsoyiannis; Associate editor C. Perrin     

Macroporosity and infiltration in blanket peat: the implications of tension disc infiltrometer measurements

Little is known about the processes of infiltration and water movement in the upper layers of blanket peat. A tension infiltrometer was used to measure hydraulic conductivity in a blanket peat in the North Pennines, England. Measurements were taken from the surface down to 20 cm in depth for peat under four different vegetation covers. It was found that macropore flow is a significant pathway for water in the upper layers of this soil type. It was also found that peat depth and surface vegetation cover were associated with macroporosity and saturated hydraulic conductivity. The proportion of macropore flow was found to be greater at 5 cm depth than at 0, 10 and 20 cm depth. Peat beneath a Sphagnum cover tends to be more permeable and a greater proportion of macropore flow can occur beneath this vegetation type. Functional macroporosity and matrix flow in the near‐surface layers of bare peat appear to have been affected by weathering processes. Comparision of results with rainfall records demonstrates that infiltration‐excess overland flow is unlikely to be a common runoff‐generating mechanism on blanket peat; rather, a saturation‐excess mechanism combined with percolation‐excess above much less permeable layers dominates the runoff response. 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.     

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.     

LISEM: A SINGLE-EVENT PHYSICALLY BASED HYDROLOGICAL AND SOIL EROSION MODEL FOR DRAINAGE BASINS. I: THEORY,INPUT AND OUTPUT

A new physically based hydrological and soil erosion model has been developed, which can be used for planning and conservation purposes: the LImburg Soil Erosion Model (LISEM). The LISEM model is one of the first examples of a physically based model that is completely incorporated in a raster Geographical Information System. This incorporation facilitates easy application in larger catchments, improves the user friendliness by avoiding conversion routines and allows remotely sensed data to be used. Processes incorporated in the model are rainfall, interception, surface storage in micro-depressions, infiltration and vertical movement of water in the soil, overland flow, channel flow, detachment by rainfall and throughfall, detachment by overland flow and transport capacity of the flow. Special attention has been given to the influence of tractor wheelings, small roads and surface sealing. Vertical movement of water in the soil is simulated using the Richard's equation. Optionally, the user can choose the Holtan or the Green–Ampt infiltration model. For the distribution flow routing, a four-point finite-difference solution of the kinematic wave is used together with Manning's equation.     

Simplified modelling of areal average infiltration at the hillslope scale

Two models for estimating expected areal‐average infiltration rate, ī, at the hillslope scale are presented. The first relies upon the condition of a negligible infiltration of surface water running downslope (run‐on process) into a previous heterogeneous soil. It is an adapted version of an earlier semi‐analytical model. The second incorporates the run‐on process and is based on a lumped approach that uses an effective saturated hydraulic conductivity. This latter was parameterized in terms of the main characteristics of rainfall and soil. Both the models were tested by comparison with the results carried out by Monte‐Carlo simulations over different soil types. It was found that the first model simulated ī with maximum errors in magnitude typically less than 10%. The second model provided similar errors in the total volume of overland flow, and the rising limb of the hydrograph experienced a distortion. Lastly, satisfactory results were obtained by comparing the model without run‐on with an empirical approach particularly accurate for fine‐textured soils. Copyright © 2002 John Wiley & Sons, Ltd.     

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.     

Numerical analysis of one-dimensional unsaturated flow in layered soils

The paper deals with numerical solutions to the Richards equation to simulate one-dimensional flow processes in the unsaturated zone of layered soil profiles. The equation is expressed in the pressure-based form and a finite-difference algorithm is developed for accurately estimating the values of the hydraulic conductivity between two neighboring nodes positioned in different soil layers, often referred to as the interlayer hydraulic conductivity. The algorithm is based upon flux conservation and continuity of pressure potential at the interface between two consecutive layers, and does not add significantly to simulation run time. The validity of the model is established for a number of test problems by comparing numerical results with the analytical solutions developed by Srivastava and Yeh²⁹ which hold for vertical infiltration towards the water table in a two-layer soil profile. The results show a significant reduction in relative mass balance errors when using the proposed model. Some specific insights into its numerical performance are also gained by comparisons with a numerical model in which the more common geometric averaging operator acts on the interlayer conductivities.