Identifying and linking source areas of flow and P transport in dairy‐grazed headwater catchments,North Island,New Zealand

We examined the applicability of the critical‐source area (CSA) concept to the dairy‐grazed 192‐ha Upper Toenepi catchment and its 8·7‐ha Kiwitahi sub‐catchment, New Zealand. We evaluated if phosphorus (P) transport from land into stream is dominated by saturation‐excess (SE) and infiltration‐excess (IE) runoff during stormflow and by sub‐surface (<1·5 m depth) flows during baseflow. We measured stream flow and shallow groundwater levels, collected monthly stream, tile drain (TDA) and groundwater samples, and flow‐proportional stream samples from the Kiwitahi sub‐catchment, and determined their dissolved reactive phosphorus (DRP) and total phosphorus (TP) concentrations. In the Kiwitahi sub‐catchment, during storm events, IE contributions were significant. Contributions from SE appeared significant in the Upper Toenepi catchment. However, in both catchments, sub‐surface contributions dominated stormflow and baseflow periods. Absence of water table at the surface and the water table gradient towards the stream indicated that P transport during events was not limited to surface runoff. The dynamics of the groundwater table and the occurrence of SE areas were influenced by proximity to the stream and hillslope positions. Baseflow accounted for 42% of the annual flow in the Kiwitahi sub‐catchment, and contributed 37 and 52% to the DRP and TP loads, respectively. The P transport during baseflow appeared equally important as P losses from CSAs during stormflow. The close resemblance in P levels between groundwater and stream samples during baseflow demonstrates the importance of shallow groundwater for stream flow. In the Upper Toenepi catchment, contributions from effluent ponds (EFFs) dominated P loads. Management strategies should focus on controlling P release from EFFs, and on decreasing Olsen P concentrations in soil to minimize leaching of P via sub‐surface flow to streams. Research is needed to quantify the role of sub‐surface flow as well as to expand management strategies to minimize P transfers during stormflow and baseflow conditions. Copyright © 2010 John Wiley & Sons, Ltd.     

An analysis of the impact of spatial variability in rainfall on runoff and sediment predictions from a distributed model

This paper investigates the effect of introducing spatially varying rainfall fields to a hydrological model simulating runoff and erosion. Pairs of model simulations were run using either spatially uniform (i.e. spatially averaged) or spatially varying rainfall fields on a 500‐m grid. The hydrological model used was a simplified version of Thales which enabled runoff generation processes to be isolated from hillslope averaging processes. Both saturation excess and infiltration excess generation mechanisms were considered, as simplifications of actual hillslope processes. A 5‐year average recurrence interval synthetic rainfall event typical of temperate climates (Melbourne, Australia) was used. The erosion model was based on the WEPP interrill equation, modified to allow nonlinear terms relating the erosion rate to rainfall or runoff‐squared. The model results were extracted at different scales to investigate whether the effects of spatially varying rainfall were scale dependent. A series of statistical metrics were developed to assess the variability due to introducing the spatially varying rainfall field. At the catchment (approximately 150 km²) scale, it was found that particularly for saturation excess runoff, model predictions of runoff were insensitive to the spatial resolution of the rainfall data. Generally, erosion processes at smaller sub‐catchment scales, particularly when the sediment generation equation had non linearity, were more sensitive to spatial rainfall variability. Introducing runon infiltration reduced the total runoff and sediment yield at all scales, and this process was also most sensitive to the rainfall resolution. Copyright © 2011 John Wiley & Sons, Ltd.     

Mapping runoff generating areas using AGNPS-VSA model

ABSTRACT

In humid regions, surface runoff is often generated by saturation-excess runoff mechanisms from relatively small variable source areas (VSAs). However, the majority of the current hydrologic models are based on infiltration-excess mechanisms. In this study, the AGricultural Non-Point Source Pollution (AGNPS) model was used to integrate the VSA concept using topographic wetness index (TWI). Both the original and AGNPS-VSA models were evaluated for a small agricultural field in Ontario, Canada. The results indicate that the AGNPS-VSA model performed better than original model. The AGNPS-VSA model predicted that only the saturated portion of the field with higher TWI values produced runoff, whereas the original AGNPS model showed uniform hydrologic response from the entire field. The results of this study are important for accurately mapping the locations of VSAs. This new model could be a powerful tool in identifying critical source areas for applying targeted best management practices to minimize pollutant loads to receiving waters.     

Development of a soil moisture‐based distributed hydrologic model for determining hydrologically based critical source areas

A simple grid cell‐based distributed hydrologic model was developed to provide spatial information on hydrologic components for determining hydrologically based critical source areas. The model represents the critical process (soil moisture variation) to run‐off generation accounting for both local and global water balance. In this way, it simulates both infiltration excess run‐off and saturation excess run‐off. The model was tested by multisite and multivariable evaluation on the 50‐km² Little River Experimental Watershed I in Georgia, U.S. and 2 smaller nested subwatersheds. Water balance, hydrograph, and soil moisture were simulated and compared to observed data. For streamflow calibration, the daily Nash‐Sutcliffe coefficient was 0.78 at the watershed outlet and 0.56 and 0.75 at the 2 nested subwatersheds. For the validation period, the Nash‐Sutcliffe coefficients were 0.79 at the watershed outlet and 0.85 and 0.83 at the 2 subwatersheds. The per cent bias was less than 15% for all sites. For soil moisture, the model also predicted the rising and declining trends at 4 of the 5 measurement sites. The spatial distribution of surface run‐off simulated by the model was mainly controlled by local characteristics (precipitation, soil properties, and land cover) on dry days and by global watershed characteristics (relative position within the watershed and hydrologic connectivity) on wet days when saturation excess run‐off was simulated. The spatial details of run‐off generation and travel time along flow paths provided by the model are helpful for watershed managers to further identify critical source areas of non‐point source pollution and develop best management practices.     

Dynamics of stormflow generation—A hillslope‐scale field study in east‐central Pennsylvania,USA

A 40 m × 20 m mowed, grass hillslope adjacent to a headwater stream within a 26‐ha watershed in east‐central Pennsylvania, USA, was instrumented to identify and map the extent and dynamics of surface saturation (areas with the water table at the surface) and surface runoff source areas. Rainfall, stream flow and surface runoff from the hillslope were recorded at 5‐min intervals from 11 August to 22 November 1998, and 13 April to 12 November 1999. The dynamics of the water table (0 to 45 cm depth from the soil surface) and the occurrence of surface runoff source areas across the hillslope were recorded using specially designed subsurface saturation and surface runoff sensors, respectively. Detailed data analyses for two rainfall events that occurred in August (57·7 mm in 150 min) and September (83·6 mm in 1265 min) 1999, illustrated the spatial and temporal dynamics of surface saturation and surface runoff source areas. Temporal data analyses showed the necessity to measure the hillslope dynamics at time intervals comparable to that of rainfall measurements. Both infiltration excess surface runoff (runoff caused when rainfall intensity exceeds soil infiltration capacity) and saturation excess surface runoff (runoff caused when soil moisture storage capacity is exceeded) source areas were recorded during these rainfall events. The August rainfall event was primarily an infiltration excess surface runoff event, whereas the September rainfall event produced both infiltration excess and saturation excess surface runoff. Occurrence and disappearance of infiltration excess surface runoff source areas during the rainfall events appeared scattered across the hillslope. Analysis of surface saturation and surface runoff data showed that not all surface saturation areas produced surface runoff that reached the stream. Emergence of subsurface flow to the surface during the post‐rainfall periods appeared to be a major flow process dominating the hillslope after the August rainfall event. Copyright © 2002 John Wiley & Sons, Ltd.     

Infiltration,runoff and sediment production in blanket peat catchments: implications of field rainfall simulation experiments

Blanket peat covers the headwaters of many major European rivers. Runoff production in upland blanket peat catchments is flashy with large flood peaks and short lag times; there is minimal baseflow. Little is known about the exact processes of infiltration and runoff generation within these upland headwaters. This paper presents results from a set of rainfall simulation experiments performed on the blanket peat moorland of the North Pennines, UK. Rainfall was simulated at low intensities (3–12 mm h^?1), typical of natural rainfall, on bare and vegetated peat surfaces. Runoff response shows that infiltration rate increases with rainfall intensity; the use of low‐intensity rainfall therefore allows a more realistic evaluation of infiltration rates and flow processes than previous studies. Overland flow is shown to be common on both vegetated and bare peat surfaces although surface cover does exert some control. Most runoff is produced within the top few centimetres of the peat and runoff response decreases rapidly with depth. Little vertical percolation takes place to depths greater than 10 cm owing to the saturation of the peat mass. This study provides evidence that the quickflow response of upland blanket peat catchments is a result of saturation‐excess overland flow generation. Rainfall–runoff response from small plots varies with season. Following warm, dry weather, rainfall tends to infiltrate more readily into blanket peat, not just initially but to the extent that steady‐state surface runoff rates are reduced and more flow takes place within the peat, albeit at shallow depth. Sediment erosion from bare peat plots tends to be supply limited. Seasonal weather conditions may affect this in that after a warm, dry spell, surface desiccation allows sediment erosion to become transport limited. Copyright © 2002 John Wiley & Sons, Ltd.     

Development of a one-parameter variable source area runoff model for ungauged basins

This research develops a one-parameter model of saturated source area dynamics and the spatial distribution of soil moisture. The single required parameter is the maximum soil moisture deficit within the catchment. The concept behind the development of the model comes from the fact that the complexity of topographically-driven runoff generation can be reduced through the use of geomorphological scaling relations. The scaling formulation allows the prediction of the dynamics of saturated source areas as a function of basin-wide soil moisture state. This model offers a number of potential advantages. Firstly, the model parameter is independent of topographic index distribution and its associated scale effects. Secondly, it may be possible to measure this single parameter using field measurements or perhaps remote sensing, which gives the model significant potential for application in ungauged basins. Finally, the fact that this parameter is a physical characteristic of the basin, estimation of this parameter avoids regionalization and parameter transferability problems. The model is tested using rainfall–runoff data from the 10.4 ha experimental catchment known as Tarrawara in Australia, the 37 km² Town Creek catchment in U.S.A., and the 620 km² Balaphi and the 850 km² Likhu sub-catchments of the Koshi river in Nepal. In sub-catchments of Koshi river, the simulation results compare favorably against the calibrated TOPMODEL both in terms of direct runoff and the spatial distribution of soil moisture state. In the Tarrawara and Town Brook catchments, simulation results compare favorably against observed storm runoff using all observed data, without calibration.     

Portable irrigation system for studying hillslope and wetland runoff generation processes

The advantage of an irrigation system is that experiments of varying precipitation duration and intensity can be performed in a controlled situation. An effective portable irrigation system is described for use in experimental plot hillslope and wetland runoff studies. The system consists of four parts:(1) a water pump; (2) a tracer reservoir; (3) a chemical feed pump; and (4) distribution hoses. Relatively uniform water application is achieved by a series of manual valves controlling water flow from the main carrier hose to the distribution hoses. The irrigation system materials are inexpensive and installation and operation costs are minor. The system was used to study runoff generation from a small‐saturated area in a spring‐fed swamp for a range of precipitation intensities and durations. The irrigation system applied a maximum intensity of 14·0mm h^?1, for a maximum duration of 180 min, to a 190m² area. This range of application incorporated all storms up to a one in three year event. The variance of the tracer load was almost three times greater for natural (60%) than with the irrigation system (22%). The irrigation system reduced the uncertainty in the pre‐event water fraction (using a two‐component hydrograph separation) from 3·4% in natural events to only 0·7%. The irrigation system design, operation, calibration and cost are presented. Copyright © 2001 John Wiley & Sons, Ltd.     

The hydrology of inland valleys in the sub‐humid zone of West Africa: rainfall‐runoff processes in the M'bé experimental watershed

Inland valleys with wet lowlands are an important water source for farming communities in the sub‐humid zone of West Africa. An inland valley and surrounding contributing watershed area located in the sub‐humid zone near M'bé in central Côte d'Ivoire was instrumented to study surface runoff and base flow mechanisms. Four flumes at different distances down the main stream and more than 100 piezometers were installed. Measurements were taken during two rainfall seasons in 1998 and 1999. Under initial wet conditions, a typical single‐peak hydrograph was observed. Under low antecedent moisture conditions, however, runoff was characterized by a double‐peaked hydrograph. The first peak, which occurred during the storm, was caused by rain falling on the saturated valley bottom. The second peak was delayed by minutes to hours from the first peak and consisted of rain flowing via the subsurface of the hydromorphic zone that surrounds the valley bottom. The duration of the delay was a function of the water table depth in the hydromorphic zone before the storm. The volume of the second peak constituted the largest portion of the stream flow. Copyright © 2003 John Wiley & Sons, Ltd.     

The influence of geomorphological position and vegetation cover on the erosional and hydrological processes on a Mediterranean hillslope

Soil erosion and runoff rates are assumed to be highly dependent on slope position. However, little knowledge exists about the hydrogeomorphological processes at the pedon scale that support this idea. In order to assess the hydrological and erosional behaviour of soils at different slope positions, simulated rainfall experiments (55 mm was applied during one hour) were carried out on a south-facing slope with underlying limestone in south-east Spain. In the mean terms, the erosion rates (9 g m² hr⁻¹) and the runoff coefficients (12%) were very low at the scale of measurement (0·25 m²). The slope position does not affect erosion rates when the measurements are carried out under extreme dry conditions during summer. The low runoff rates found in summer under thunderstorms of high intensity (5 year return period) and the runon into surfaces with higher infiltration rates resulted in no detectable direct surface runoff (Hortonian) at the slope scale. This implies that erosion as a consequence of surface overland flow will only take place during events of high magnitude (55 mm hr⁻¹) and low frequency (>5 years). Vegetation is the most important factor determining the soil erosion and runoff rates within the slope. © 1998 John Wiley & Sons, Ltd.     

Predicting saturation‐excess runoff distribution with a lumped hillslope model: SWAT‐HS

Saturation‐excess runoff is the major runoff mechanism in humid well‐vegetated areas where infiltration rates often exceed rainfall intensity. Although the Soil and Water Assessment Tool (SWAT) is one of the most widely used models, it predicts runoff based mainly on soil and land use characteristics, and is implicitly an infiltration‐excess runoff type of model. Previous attempts to incorporate the saturation‐excess runoff mechanism in SWAT fell short due to the inability to distribute water from one hydrological response unit to another. This paper introduces a modified version of SWAT, referred to as SWAT‐Hillslope (SWAT‐HS). This modification improves the simulation of saturation‐excess runoff by redefining hydrological response units based on wetness classes and by introducing a surface aquifer with the ability to route interflow from “drier” to “wetter” wetness classes. Mathematically, the surface aquifer is a nonlinear reservoir that generates rapid subsurface stormflow as the water table in the surface aquifer rises. The SWAT‐HS model was tested in the Town Brook watershed in the upper reaches of the West Branch Delaware River in the Catskill region of New York, USA. SWAT‐HS predicted discharge well with a Nash‐Sutcliffe Efficiency of 0.68 and 0.87 for daily and monthly time steps. Compared to the original SWAT model, SWAT‐HS predicted less surface runoff and groundwater flow and more lateral flow. The saturated areas predicted by SWAT‐HS were concentrated in locations with a high topographic index and were in agreement with field observations. With the incorporation of topographic characteristics and the addition of the surface aquifer, SWAT‐HS improved streamflow simulation and gave a good representation of saturated areas on the dates that measurements were available. SWAT‐HS is expected to improve water quality model predictions where the location of the surface runoff matters.     

Improving the spatial representation of soil properties and hydrology using topographically derived initialization processes in the SWAT model

Topography exerts critical controls on many hydrologic, geomorphologic and biophysical processes. However, many watershed modelling systems use topographic data only to define basin boundaries and stream channels, neglecting opportunities to account for topographic controls on processes such as soil genesis, soil moisture distributions and hydrological response. Here, we demonstrate a method that uses topographic data to adjust spatial soil morphologic and hydrologic attributes: texture, depth to the C‐horizon, saturated conductivity, bulk density, porosity and the water capacities at field (33 kpa) and wilting point (1500 kpa) tensions. As a proof of concept and initial performance test, the values of the topographically adjusted soil parameters and those from the Soil Survey Geographic Database (SSURGO; available at 1 : 20 000 scale) were compared with measured soil pedon pit data in the Grasslands Soil and Water Research Lab watershed in Riesel, TX. The topographically adjusted soils were better correlated with the pit measurements than were the SSURGO values. We then incorporated the topographically adjusted soils into an initialization of the Soil and Water Assessment Tool model for 15 Riesel research watersheds to investigate how changes in soil properties influence modelled hydrological responses at the field scale. The results showed that the topographically adjusted soils produced better runoff predictions in 50% of the fields, with the SSURGO soils performing better in the remainder. In addition, the a priori adjusted soils result in fewer calibrated model parameters. These results indicate that adjusting soil properties based on topography can result in more accurate soil characterization and, in some cases, improve model performance. Copyright © 2016 John Wiley & Sons, Ltd.     

Quantifying aggregation and change in runoff source in accordance with catchment area increase in a forested headwater catchment

There has been a great deal of research interest regarding changes in flow path/runoff source with increases in catchment area. However, there have been very few quantitative studies taking subscale variability and convergence of flow path/runoff source into account, especially in relation to headwater catchments. This study was performed to elucidate how the contributions and discharge rates of subsurface water (water in the soil layer) and groundwater (water in fractured bedrock) aggregate and change with catchment area increase, and to elucidate whether the spatial variability of the discharge rate of groundwater determines the spatial variability of stream discharge or groundwater contribution. The study area was a 5‐km² forested headwater catchment in Japan. We measured stream discharge at 113 points and water chemistry at 159 points under base flow conditions. End‐member mixing analysis was used to separate stream water into subsurface water and groundwater. The contributions of both subsurface water and groundwater had large variability below 1 km². The contribution of subsurface water decreased markedly, while that of groundwater increased markedly, with increases in catchment area. The specific discharge of subsurface water showed a large degree of variability and decreased with catchment area below 0.1 km², becoming almost constant above 0.1 km². The specific discharge of groundwater showed large variability below 1 km² and increased with catchment area. These results indicated that the variabilities of stream discharge and groundwater contribution corresponded well with the variability of the discharge rate of groundwater. However, below 0.1 km², it was necessary to consider variations in the discharge rates of both subsurface water and groundwater. Copyright © 2016 John Wiley & Sons, Ltd.     

Rainfall intensity affects runoff responses in a semi-arid catchment

The source and hydrochemical makeup of a stream reflects the connectivity between rainfall, groundwater, the stream, and is reflected to water quantity and quality of the catchment. However, in a semi-arid, thick, loess covered catchment, temporal variation of stream source and event associated behaviours are lesser known. Thus, the isotopic and chemical hydrographs in a widely distributed, deep loess, semi-arid catchment of the northern Chinese Loess Plateau were characterized to determine the source and hydrochemical behaviours of the stream during intra-rainfall events. Rainfall and streamflow were sampled during six hydrologic events coupled with measurements of stream baseflow and groundwater. The deuterium isotope (²H), major ions (Cl⁻, SO₄²⁻, NO₃⁻, Ca²⁺, K⁺, Mg²⁺, and Na⁺) were evaluated in water samples obtained during rainfall events. Temporal variation of ²H and Cl⁻ measured in the groundwater and stream baseflow prior to rainfall was similar; however, the isotope compositions of the streamflow fluctuated significantly and responded quickly to rainfall events, likely due to an infiltration excess, overland dominated surface runoff during torrential rainfall events. Time source separation using ²H demonstrated greater than 72% on average, the stream composition was event water during torrential rainfall events, with the proportion increasing with rainfall intensity. Solutes concentrations in the stream had loglinear relationships with stream discharge, with an outling anomaly with an example of an intra-rainfall event on Oct. 24, 2015. Stream Cl⁻ behaved nonconservative during rainfall events, temporal variation of Cl⁻ indicated a flush and washout at the onset of small rainfall events, a dilution but still high concentration pattern in high discharge and old water dominated in regression flow period. This study indicates rainfall intensity affects runoff responses in a semi-arid catchment, and the stored water in the thick, loess covered areas was less connected with stream runoff. Solute transport may threaten water quality in the area, requiring further analysis of the performance of the eco-restoration project.     

Using a topographic index to distribute variable source area runoff predicted with the SCS curve‐number equation

Because the traditional Soil Conservation Service curve‐number (SCS‐CN) approach continues to be used ubiquitously in water quality models, new application methods are needed that are consistent with variable source area (VSA) hydrological processes in the landscape. We developed and tested a distributed approach for applying the traditional SCS‐CN equation to watersheds where VSA hydrology is a dominant process. Predicting the location of source areas is important for watershed planning because restricting potentially polluting activities from runoff source areas is fundamental to controlling non‐point‐source pollution. The method presented here used the traditional SCS‐CN approach to predict runoff volume and spatial extent of saturated areas and a topographic index, like that used in TOPMODEL, to distribute runoff source areas through watersheds. The resulting distributed CN–VSA method was applied to two subwatersheds of the Delaware basin in the Catskill Mountains region of New York State and one watershed in south‐eastern Australia to produce runoff‐probability maps. Observed saturated area locations in the watersheds agreed with the distributed CN–VSA method. Results showed good agreement with those obtained from the previously validated soil moisture routing (SMR) model. When compared with the traditional SCS‐CN method, the distributed CN–VSA method predicted a similar total volume of runoff, but vastly different locations of runoff generation. Thus, the distributed CN–VSA approach provides a physically based method that is simple enough to be incorporated into water quality models, and other tools that currently use the traditional SCS–CN method, while still adhering to the principles of VSA hydrology. Copyright © 2004 John Wiley & Sons, Ltd.     

Variations in runoff,sediment load,and their relationship for a major sediment source area of the Jialing River basin,southern China

Investigation of the variations in runoff, sediment load, and their dynamic relation is conducive to understanding hydrological regime changes and supporting channel regulation and fluvial management. This study is undertaken in the Xihanshui catchment, which is known for its high sediment-laden in the Jialing River of the Yangtze River basin, southern China, to evaluate the change characteristics of runoff, sediment load, and their relationship at multi-temporal scales from 1966 to 2016. The results showed that runoff changed significantly for more months, whereas the significant changes in monthly sediment load occurred from April to September. The contributions of runoff in summer and autumn and sediment load in summer to their annual value changes were greater. Annual runoff and sediment load in the Xihanshui catchment both exhibited significant decreasing trends (p < 0.05) with a significant mutation in 1993 (p < 0.05). The average annual runoff in the change period (1994–2016) decreased by 49.58% and annual sediment load displayed a substantial decline with a reduction of 77.77% in comparison with the reference period (1966–1993) due to climate change and intensive human activity. The power functions were satisfactory to describe annual and extreme monthly runoff–sediment relationships, whereas the monthly runoff–sediment relationship and extreme monthly sediment-runoff relationship were changeable. Spatially, annual runoff–sediment relationship alteration could be partly attributed to sediment load changes in the upstream area and runoff variations in the downstream region. Three quantitative methods revealed that the main driver for significant reductions of annual runoff and sediment load is the human activity dominated by soil and water conservation measures, while climate change only contributed 22.73%–38.99% (mean 32.07%) to the total runoff reduction and 3.39%–35.56% (mean 17.32%) to the total decrease in sediment load.     

Scale relationships in hillslope runoff and erosion

Eight runoff plots, located within a small catchment within the Walnut Gulch Experimental Watershed, southern Arizona, were constructed to test the argument that sediment yield (kg m^?2) decreases as plot length increases. The plots ranged in length from 2 m to 27·78 m. Runoff and sediment loss from these plots were obtained for ten natural storm events. The pattern of sediment yield from these plots conforms to the case in which sediment yield first increases as plot length increases, but then subsequently decreases. Data from the present experiment indicate that maximum sediment yield would occur from a plot 7 m long. Analysis of both runoff and sediment yield from the plots indicates that the relationship of sediment yield to plot length derives both from the limited travel distance of individual entrained particles and from a decline in runoff coefficient as plot length increases. Particle‐size analysis of eroded sediment confirms the role of travel distance in controlling sediment yield. Whether in response to the finite travel distance of entrained particles or the relationship of runoff coefficient to plot length, the experiment clearly demonstrates that the erosion rates for hillslopes and catchments cannot be simply extrapolated from plot measurements, and that alternative methods for estimating large‐area erosion rates are required. Copyright © 2006 John Wiley & Sons, Ltd.     

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.