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.     

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.     

Changes in overland flow and infiltration after a rangeland fire in a Mediterranean scrubland

Changes in overland flow and infiltration after a wildfire (summer 1989) in a typical Mediterranean scrubland were measured during the winters of 1990, 1991, 1992 and 1995 by means of simulated rainfall. Infiltration increases gradually from 1990 (sixth months after the forest fire) to 1995 (five and a half years after the forest fire). Overland flow decreases from 45% of rainfall after the forest fire to less than 6% five and a half years later. The reduction in overland flow was greatest in the first two years after the fire because of the quick recovery of vegetation. The steady-state infiltration capacity increased every year after the fire. Runoff and infiltration changes are mainly determined by the gradual recovery of vegetation. © 1998 John Wiley & Sons, Ltd.     

Robert E. Horton's perceptual model of infiltration processes

Robert E. Horton is best known as the originator of the infiltration excess overland flow concept for storm hydrograph analysis and prediction, which, in conjunction with the unit hydrograph concept, provided the foundation for engineering hydrology for several decades. Although these concepts, at least in their simplest form, have been largely superseded, a study of Horton's archived scientific papers reveals that his perceptual model of infiltration processes and appreciation of scale problems in modelling were far more sophisticated and complete than normally presented in hydrological texts. His understanding of surface controls on infiltration remain relevant today. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

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.     

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.     

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     

A new modeling approach for simulating microtopography-dominated,discontinuous overland flow on infiltrating surfaces

Realistic modeling of discontinuous overland flow on irregular topographic surfaces has been proven to be a challenge. This study is aimed to develop a new modeling framework to simulate the discontinuous puddle-to-puddle (P2P) overland flow dynamics for infiltrating surfaces with various microtopographic characteristics. In the P2P model, puddles were integrated in a well-delineated, cascaded drainage system to facilitate explicit simulation of their dynamic behaviors and interactions. Overland flow and infiltration were respectively simulated by using the diffusion wave model and a modified Green–Ampt model for the DEM-derived flow drainage network that consisted of a series of puddle-based units (PBUs). The P2P model was tested by using a series of data from laboratory overland flow experiments for various microtopography, soil, and rainfall conditions. The modeling results indicated that the hierarchical relationships and microtopographic properties of puddles significantly affected their connectivity, filling–spilling dynamics, and the associated threshold flow. Surface microtopography and rainfall characteristics also exhibited strong influences on the spatio-temporal distributions of infiltration rates, runoff fluxes, and unsaturated flow. The model tests demonstrated its applicability in simulating microtopography-dominated overland flow on infiltrating surfaces.     

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.     

A sensitivity analysis of Hortonian flow

We present a sensitivity analysis for infiltration excess (Hortonian) overland flow based on a classic laboratory experiment by Smith and Woolhiser [Smith RE, Woolhiser DA. Overland flow on an infiltrating surface. Water Resour Res 1971;7(4):899–913]. The model components of the compartment approach are comprised of a diffusive wave approximation to the Saint–Venant equations for overland flow, a Richards model for flow in the variably saturated zone, and an interface coupling concept that combines the two components. In the coupling scheme a hydraulic interface is introduced to allow the definition of an exchange flux between the surface and the unsaturated zone. The effects of friction processes, soil capillarity, hydraulic interface, and vertical soil discretization on both infiltration and runoff prediction are investigated in detail. The corresponding sensitivity analysis is conducted using a small-perturbation method. As a result the importance of the hydraulic processes and related parameters are evaluated for the coupled hydrosystem.     

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.     

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.     

A distributed TOPMODEL for modelling impacts of land‐cover change on river flow in upland peatland catchments

There is global concern about headwater management and associated impacts on river flow. In many wet temperate zones peatlands can be found covering headwater catchments. In the UK there is major concern about how environmental change, driven by human interventions, has altered the surface cover of headwater blanket peatlands. However, the impact of such land‐cover changes on river flow is poorly understood. In particular, there is poor understanding of the impacts of different spatial configurations of bare peat or well‐vegetated, restored peat on river flow peaks in upland catchments. In this paper, a physically based, distributed and continuous catchment hydrological model was developed to explore such impacts. The original TOPMODEL, with its process representation being suitable for blanket peat catchments, was utilized as a prototype acting as the basis for the new model. The equations were downscaled from the catchment level to the cell level. The runoff produced by each cell is divided into subsurface flow and saturation‐excess overland flow before an overland flow calculation takes place. A new overland flow module with a set of detailed stochastic algorithms representing overland flow routing and re‐infiltration mechanisms was created to simulate saturation‐excess overland flow movement. The new model was tested in the Trout Beck catchment of the North Pennines of England and found to work well in this catchment. The influence of land cover on surface roughness could be explicitly represented in the model and the model was found to be sensitive to land cover. Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

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.     

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.     

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.     

Mechanisms responsible for streamflow generation on a small,salt-affected and deeply weathered hillslope

This paper considers the contributions of overland flow, throughflow and deep seepage to the generation of streamflow in a salt-affected, deeply weathered landscape. Runoff mechanisms on a small hillslope in south-western Australia were dependent on the extent and development of variable source areas. In winter, streamflow generation was controlled by returnflow, saturation overland flow and throughflow. In summer, post-ponding, infiltration-excess and saturation overland flow dominated. The extent of the variable source area and the magnitude of streamflow were due to antecedent soil moisture, rainfall and slope morphology. Concave hillslope sections accumulated soil moisture due to both saturated and unsaturated lateral flow processes. Throughflow provided the mechanism and vehicle for solute movement from the groundwater discharge area to the stream. However, discharge from the deep aquifer was the primary mechanism responsible for soil salinity and maintaining the core of the variable source area. Estimates of throughflow which only take account of soil-water movement and disregard returnflow, will underestimate the magnitude of throughflow.