Field investigation and modelling of coupled stream discharge and shallow water‐table dynamics in a small Mediterranean catchment (Sardinia)

In the semi‐arid Mediterranean environment, the rainfall–runoff relationships are complex because of the markedly irregular patterns in rainfall, the seasonal mismatch between evaporation and rainfall, and the spatial heterogeneity in landscape properties. Watersheds often display considerable non‐linear threshold behavior, which still make runoff generation an open research question. Our objectives in this context were: to identify the primary processes of runoff generation in a small natural catchment; to test whether a physically based model, which takes into consideration only the primary processes, is able to predict spatially distributed water‐table and stream discharge dynamics; and to use the hydrological model to increase our understanding of runoff generation mechanisms. The observed seasonal dynamics of soil moisture, water‐table depth, and stream discharge indicated that Hortonian overland‐flow was negligible and the main mechanism of runoff generation was saturated subsurface‐flow. This gives rise to base‐flow, controls the formation of the saturated areas, and contributes to storm‐flow together with saturation overland‐flow. The distributed model, with a 1D scheme for the kinematic surface‐flow, a 2D sub‐horizontal scheme for the saturated subsurface‐flow, and ignoring the unsaturated flow, performed efficiently in years when runoff volume was high and medium, although there was a smoothing effect on the observed water‐table. In dry years, small errors greatly reduced the efficiency of the model. The hydrological model has allowed to relate the runoff generation mechanisms with the land‐use. The forested hillslopes, where the calibrated soil conductivity was high, were never saturated, except at the foot of the slopes, where exfiltration of saturated subsurface‐flow contributed to storm‐flow. Saturation overland‐flow was only found near the streams, except when there were storm‐flow peaks, when it also occurred on hillslopes used for pasture, where soil conductivity was low. The bedrock–soil percolation, simulated by a threshold mechanism, further increased the non‐linearity of the rainfall–runoff processes. Copyright © 2013 John Wiley & Sons, Ltd.     

Assessing the performance of a physically based hydrological model using a proxy‐catchment approach in an agricultural environment

Physically based models are useful frameworks for testing intervention strategies designed to reduce elevated sediment loads in agricultural catchments. Evaluating the success of these strategies depends on model accuracy, generally established by a calibration and evaluation process. In this contribution, the physically based SHETRAN model was assessed in two similar U.K. agricultural catchments. The model was calibrated on the Blackwater catchment (18 km²) and evaluated in the adjacent Kit Brook catchment (22 km²) using 4 years of 15 min discharge and suspended sediment flux data. Model sensitivity to changes in single and multiple combinations of parameters and sensitivity to changes in digital elevation model resolution were assessed. Model flow performance was reasonably accurate with a Nash–Sutcliffe efficiency coefficient of 0.78 in Blackwater and 0.60 in Kit Brook. In terms of event prediction, the mean of the absolute percentage of difference (μAbs_diff) between measured and simulated flow volume (Qv), peak discharge (Qp), sediment yield (Sy), and peak sediment flux (Sp) showed larger values in Kit Brook (48% [Qv], 66% [Qp], 298% [Sy], and 438% [Sp]) compared with the Blackwater catchment (30% [Qv], 41% [Qp], 106% [Sy], and 86% [Sp]). Results indicate that SHETRAN can produce reasonable flow prediction but performs less well in estimation of sediment flux, despite reasonably similar hydrosedimentary behaviour between catchments. The sensitivity index showed flow volume sensitive to saturated hydraulic conductivity and peak discharge to the Strickler coefficient; sediment yield was sensitive to the overland flow erodibility coefficient and peak sediment flux to raindrop/leaf soil erodibility coefficient. The multiparameter sensitivity analysis showed that different combinations of parameters produced similar model responses. Model sensitivity to grid resolution presented similar flow volumes for different digital elevation model resolutions, whereas event peak and duration (for both flow and sediment flux) were highly sensitive to changes in grid size.     

Parameter estimation in semi‐distributed hydrological catchment modelling using a multi‐criteria objective function

Output generated by hydrologic simulation models is traditionally calibrated and validated using split‐samples of observed time series of total water flow, measured at the drainage outlet of the river basin. Although this approach might yield an optimal set of model parameters, capable of reproducing the total flow, it has been observed that the flow components making up the total flow are often poorly reproduced. Previous research suggests that notwithstanding the underlying physical processes are often poorly mimicked through calibration of a set of parameters hydrologic models most of the time acceptably estimates the total flow. The objective of this study was to calibrate and validate a computer‐based hydrologic model with respect to the total and slow flow. The quick flow component used in this study was taken as the difference between the total and slow flow. Model calibrations were pursued on the basis of comparing the simulated output with the observed total and slow flow using qualitative (graphical) assessments and quantitative (statistical) indicators. The study was conducted using the Soil and Water Assessment Tool (SWAT) model and a 10‐year historical record (1986–1995) of the daily flow components of the Grote Nete River basin (Belgium). The data of the period 1986–1989 were used for model calibration and data of the period 1990–1995 for model validation. The predicted daily average total flow matched the observed values with a Nash–Sutcliff coefficient of 0·67 during calibration and 0·66 during validation. The Nash–Sutcliff coefficient for slow flow was 0·72 during calibration and 0·61 during validation. Analysis of high and low flows indicated that the model is unbiased. A sensitivity analysis revealed that for the modelling of the daily total flow, accurate estimation of all 10 calibration parameters in the SWAT model is justified, while for the slow flow processes only 4 out of the set of 10 parameters were identified as most sensitive. Copyright © 2007 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.     

The use of agent based modelling techniques in hydrology: determining the spatial and temporal origin of channel flow in semi‐arid catchments

Information on the spatial and temporal origin of runoff entering the channel during a storm event would be valuable in understanding the physical dynamics of catchment hydrology; this knowledge could be used to help design flood defences and diffuse pollution mitigation strategies. The majority of distributed hydrological models give information on the amount of flow leaving a catchment and the pattern of fluxes within the catchment. However, these models do not give any precise information on the origin of runoff within the catchment. The spatial and temporal distribution of runoff sources is particularly intricate in semi‐arid catchments, where there are complex interactions between runoff generation, transmission and re‐infiltration over short temporal scales. Agents are software components that are capable of moving through and responding to their local environment. In this application, the agents trace the path taken by water through the catchment. They have information on their local environment and on the basis of this information make decisions on where to move. Within a given model iteration, the agents are able to stay in the current cell, infiltrate into the soil or flow into a neighbouring cell. The information on the current state of the hydrological environment is provided by the environment generator. In this application, the Connectivity of Runoff Model (CRUM) has been used to generate the environment. CRUM is a physically based, distributed, dynamic hydrology model, which considers the hydrological processes relevant for a semi‐arid environment at the temporal scale of a single storm event. During the storm event, agents are introduced into the model across the catchment; they trace the flows of water and store information on the flow pathways. Therefore, this modelling approach is capable of giving a novel picture of the temporal and spatial dynamics of flow generation and transmission during a storm event. This is possible by extracting the pathways taken by the agents at different time slices during the storm. The agent based modelling approach has been applied to two small catchments in South East Spain. The modelling approach showed that the two catchments responded differently to the same rainfall event due to the differences in the runoff generation and overland flow connectivity between the two catchments. The model also showed that the time of travel to the nearest flow concentration is extremely important for determining the connectivity of a point in the landscape with the catchment outflow. Copyright © 2007 John Wiley & Sons, Ltd.     

The effect of forest cover on peak flow and sediment discharge—an integrated field and modelling study in central–southern Chile

An integrated field and modelling study was carried out on the 35‐ha La Reina catchment, Chile, to test the hypothesis that the effect of forest cover on flood peaks becomes less important as the size of the hydrological event increases. Meteorological and discharge data were measured at the catchment before and after the pine plantation that covered 80% of the catchment area was logged. Analysis of the measured response of the catchment provides support for the hypothesis but is not conclusive. Therefore, modelling of the catchment using 1000 years of generated rainfall data representative of the current conditions was carried out for the forested and logged states. The simulations show that the absolute difference in discharge between the two cases remains approximately constant as the discharge increases: thus as a percentage of discharge it decreases. This relative convergence appears to become significant at return periods of greater than approximately 10 years. Tests with different hypothetical soil depths for the forested and logged catchments show an absolute convergence in discharge between the two cases for shallow soils and no convergence for deep soils. Sediment transport simulations show that forest cover provides a clear benefit in protecting the soil from erosion. Copyright © 2010 John Wiley & Sons, Ltd.     

Parameterization and multi‐criteria calibration of a distributed storm flow model applied to a Mediterranean agricultural catchment

The principal challenge in the parameterization of storm flow models for agricultural catchments with an artificial drainage network and fields with different degrees of tillage lies in the parsimonious definition of distributed model parameters in a way that reduces the number of calibration parameters to a justifiable minimum. This paper presents a comprehensive case study for the parameter estimation of a distributed storm flow model applied to an agricultural catchment (0.91 km²) in the Mediterranean region. Model parameterization was combined with procedures for multi‐criteria, multi‐storm calibration, where we automatically calibrated three parameters related to flow velocity and infiltration, and compared single and multi‐storm criteria that are based on discharge volume, peak flow, and the Nash–Sutcliffe coefficient. Multi‐storm calibration yielded a set of parameter values for the simulation batch with best multi‐storm overall performance, which are close to the median values in the pre‐calibration of individual storms. Our results suggest that flow velocities and proportionality of the channel infiltration rate do not vary significantly over the course of 11 years. Copyright © 2012 John Wiley & Sons, Ltd.     

Alteration of hydrological processes and streamflow with juniper (Juniperus virginiana) encroachment in a mesic grassland catchment

To understand the effect of woody plant encroachment on hydrological processes of mesic grasslands, we quantified infiltration capacity in situ, the temporal changes in soil water storage, and streamflow of a grassland catchment and a catchment heavily encroached by juniper (Juniperus virginiana, eastern redcedar) in previously cultivated, non‐karst substrate grasslands in north‐central Oklahoma for 3 years. The initial and steady‐state infiltration rates under the juniper canopy were nearly triple to that of the grassland catchment and were intermediate in the intercanopy spaces within the encroached catchment. Soil water content and soil water storage on the encroached catchment were generally lower than on the grassland catchment, especially when preceding the seasons of peak rainfall in spring and fall. Frequency and magnitude of streamflow events were reduced in the encroached catchment. Annual runoff coefficients for the encroached catchment averaged 2.1%, in contrast to 10.6% for the grassland catchment. Annual streamflow duration ranged from 80 to 250 h for the encroached catchment compared with 600 to 800 h for the grassland catchment. Our results showed that the encroachment of juniper into previously cultivated mesic grasslands fundamentally alters catchment hydrological function. Rapid transformation of mesic grassland to a woodland state with juniper encroachment, if not confined, has the potential to drastically reduce soil water, streamflow and flow duration of ephemeral streams in the Southern Great Plains. Copyright © 2013 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.     

A GIS‐based variable source area hydrology model

Effective control of nonpoint source pollution from contaminants transported by runoff requires information about the source areas of surface runoff. Variable source hydrology is widely recognized by hydrologists, yet few methods exist for identifying the saturated areas that generate most runoff in humid regions. The Soil Moisture Routing model is a daily water balance model that simulates the hydrology for watersheds with shallow sloping soils. The model combines elevation, soil, and land use data within the geographic information system GRASS, and predicts the spatial distribution of soil moisture, evapotranspiration, saturation‐excess overland flow (i.e., surface runoff), and interflow throughout a watershed. The model was applied to a 170 hectare watershed in the Catskills region of New York State and observed stream flow hydrographs and soil moisture measurements were compared to model predictions. Stream flow prediction during non‐winter periods generally agreed with measured flow resulting in an average r² of 0·73, a standard error of 0·01 m³/s, and an average Nash‐Sutcliffe efficiency R² of 0·62. Soil moisture predictions showed trends similar to observations with errors on the order of the standard error of measurements. The model results were most accurate for non‐winter conditions. The model is currently used for making management decisions for reducing non‐point source pollution from manure spread fields in the Catskill watersheds which supply New York City's drinking water. Copyright © 1999 John Wiley & Sons, Ltd.     

Anomalous behaviour of specific electrical conductivity at a karst spring induced by variable catchment boundaries: the case of the Podstenjšek spring,Slovenia

Anomalous behaviour of specific electrical conductivity (SEC) was observed at a karst spring in Slovenia during 26 high‐flow events in an 18‐month monitoring period. A conceptual model explaining this anomalous SEC variability is presented and reproduced by numerical modelling, and the practical relevance for source protection zoning is discussed. After storm rainfall, discharge increases rapidly, which is typical for karst springs. SEC displays a first maximum during the rising limb of the spring hydrograph, followed by a minimum indicating the arrival of freshly infiltrated water, often confirmed by increased levels of total organic carbon (TOC). The anomalous behaviour starts after this SEC minimum, when SEC rises again and remains elevated during the entire high‐flow period, typically 20–40 µS/cm above the baseflow value. This is explained by variable catchment boundaries: When the water level in the aquifer rises, the catchment expands, incorporating zones of groundwater with higher SEC, caused by higher unsaturated zone thickness and subtle lithologic changes. This conceptual model has been checked by numerical investigations. A generalized finite‐difference model including high‐conductivity cells representing the conduit network (“discrete‐continuum approach”) was set up to simulate the observed behaviour of the karst system. The model reproduces the shifting groundwater divide and the nearly simultaneous increase of discharge and SEC during high‐flow periods. The observed behaviour is relevant for groundwater source protection zoning, which requires reliable delineation of catchment areas. Anomalous behaviour of SEC can point to variable catchment boundaries that can be checked by tracer tests during different hydrologic conditions. Copyright © 2011 John Wiley & Sons, Ltd.     

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

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

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.     

Daily streamflow modelling and assessment based on the curve‐number technique

A cell‐based long‐term hydrological model (CELTHYM) that can be integrated with a geographical information system (GIS) was developed to predict continuous stream flow from small agricultural watersheds. The CELTHYM uses a cell‐by‐cell soil moisture balance approach. For surface runoff estimation, the curve number technique considering soil moisture on a daily basis was used, and release rate was used to estimate baseflow. Evapotranspiration was computed using the FAO modified Penman equation that considered land‐use‐based crop coefficients, soil moisture and the influence of topography on radiation. A rice paddy field water budget model was also adapted for the specific application of the model to East Asia. Model sensitivity analysis was conducted to obtain operational information about the model calibration parameters. The CELTHYM was calibrated and verified with measured runoff data from the WS#1 and WS#3 watersheds of the Seoul National University, Department of Agricultural Engineering, in Hwaseong County, Kyounggi Province, South Korea. The WS#1 watershed is comprised of about 35·4% rice paddy fields and 42·3% forest, whereas the WS#3 watershed is about 85·0% forest and 11·5% rice paddy fields. The CELTHYM was calibrated for the parameter release rate, K, and soil moisture storage coefficient, STC, and results were compared with the measured runoff data for 1986. The validation results for WS#1 considering all daily stream flow were poor with R², E² and RMSE having values of 0·40, ?6·63 and 9·69 (mm), respectively, but validation results for days without rainfall were statistically significant (R² = 0·66). Results for WS#3 showed good agreement with observed data for all days, and R², E² and RMSE were 0·92, 0·91 and 2·23 (mm), respectively, suggesting potential for CELTHYM application to other watersheds. The direct runoff and water balance components for watershed WS#1 with significant areas of paddy fields did not perform well, suggesting that additional study of these components is needed. Copyright © 2002 John Wiley & Sons, Ltd.     

Role of temporal resolution of meteorological inputs for process‐based snow modelling

Accurate snow accumulation and melt simulations are crucial for understanding and predicting hydrological dynamics in mountainous settings. As snow models require temporally varying meteorological inputs, time resolution of these inputs is likely to play an important role on the model accuracy. Because meteorological data at a fine temporal resolution (~1 hr) are generally not available in many snow‐dominated settings, it is important to evaluate the role of meteorological inputs temporal resolution on the performance of process‐based snow models. The objective of this work is to assess the loss in model accuracy with temporal resolution of meteorological inputs, for a range of climatic conditions and topographic elevations. To this end, a process‐based snow model was run using 1‐, 3‐, and 6‐hourly inputs for wet, average, and dry years over Boise River Basin (6,963 km²), which spans rain dominated (≤1,400 m), rain–snow transition (>1,400 and ≤1,900 m), snow dominated below tree line (>1,900 and ≤2,400 m), and above tree line (>2,400 m) elevations. The results show that sensitivity of the model accuracy to the inputs time step generally decreases with increasing elevation from rain dominated to snow dominated above tree line. Using longer than hourly inputs causes substantial underestimation of snow cover area (SCA) and snow water equivalent (SWE) in rain‐dominated and rain–snow transition elevations, due to the precipitation phase mischaracterization. In snow‐dominated elevations, the melt rate is underestimated due to errors in estimation of net snow cover energy input. In addition, the errors in SCA and SWE estimates generally decrease toward years with low snow mass, that is, dry years. The results indicate significant increases in errors in estimates of SCA and SWE as the temporal resolution of meteorological inputs becomes coarser than an hour. However, use of 3‐hourly inputs can provide accurate estimates at snow‐dominated elevations. The study underscores the need to record meteorological variables at an hourly time step for accurate process‐based snow modelling.     

A tracer‐based assessment of hydrological pathways at different spatial scales in a mesoscale Scottish catchment

Geochemically based hydrograph separation techniques were used in a preliminary assessment to infer how runoff processes change with landscape characteristics and spatial scale (1–233 km²) within a mesoscale catchment in upland Scotland. A two‐component end‐member mixing analysis (EMMA) used Gran alkalinity as an assumed conservative tracer. Analysis indicated that, at all scales investigated, acidic overland flow and shallow subsurface storm flows from the peaty soils covering the catchment headwaters dominated storm runoff generation. The estimated groundwater contribution to annual runoff varied from 30% in the smallest (ca 1 km²) peat‐dominated headwater catchment with limited groundwater storage, to >60% in larger catchments (>30 km²) with greater coverage of more freely draining soils and more extensive aquifers in alluvium and other drift. This simple approach offers a useful, integrated conceptualization of the hydrological functioning in a mesoscale catchment, which can be tested and further refined by focused modelling and process‐based research. However, even as it stands, the simple conceptualization of system behaviour will have significant utility as a tool for communicating hydrological issues in a range of planning and management decisions. Copyright © 2002 John Wiley & Sons, Ltd.     

Evaluating explicit and implicit routing for watershed hydro‐ecological models of forest hydrology at the small catchment scale

This paper explores the behaviour and sensitivity of a watershed model used for simulating lateral soil water redistribution and runoff production. In applications such as modelling the effects of land‐use change in small headwater catchments, interactions between soil moisture, runoff and ecological processes are important. Because climate, soil and canopy characteristics are spatially variable, both the pattern of soil moisture and the associated outflow must be represented in modelling these processes. This study compares implicit and explicit routing approaches to modelling the evolution of soil moisture pattern and spatially variable runoff production. It also addresses the implications of using different landscape partitioning strategies. This study presents the results of calibration and application of these different routing and landscape partitioning approaches on a 60 ha forested watershed in Western Oregon. For comparison, the different approaches are incorporated into a physically based hydro‐ecological model, RHESSys, and the resulting simulated soil moisture, runoff production and sensitivity to unbiased error are examined. Results illustrate that both routing approaches can be calibrated to achieve a reasonable fit between observed and modelled outflow. Calibrated values for effective watershed hydraulic conductivity are higher for the explicit routing approach, which illustrates differences between the two routing approaches in their representation of internal watershed dynamics. The explicit approach illustrates a seasonal shift in drainage organization from watershed to more local control as climate goes from a winter wet to a summer dry period. Assumptions used in the implicit approach maintain the same pattern of drainage organization throughout the season. The implicit approach is also more sensitive to random error in soil and topographic input information, particularly during wetter periods. Comparison between the two routing approaches illustrates the advantage of the explicit routing approach, although the loss of computational efficiency associated with the explicit routing approach is noted. To compare different strategies for partitioning the landscape, the use of a non‐grid‐based method of partitioning is introduced and shown to be comparable to grid‐based partitioning in terms of simulated soil moisture and runoff production. Copyright © 2001 John Wiley & Sons, Ltd.