Spatial patterns of sediment delivery to valley floors: sensitivity to sediment transport capacity and hillslope hydrology relations

There is considerable interest in large‐scale spatial patterns of sediment transport in catchments, and this topic is often approached using terrain‐based modelling. In such models topography influences the discharge of overland flow and its sediment transport capacity. The sediment transport capacity of overland flow is commonly expressed as a power function of slope and discharge (i.e. q_s=k₁q^βS^γ). The relationship between discharge and contributing area can also be expressed as a power function. Several reviews reveal a limited range of values for the two exponents β and γ. In this paper we examine the sensitivity of catchment‐scale patterns of sediment delivery to valley floors to a range of sediment transport capacity and hillslope hydrology parameterizations, using two catchments on the southern tablelands of New South Wales. The results indicate that, over the limited range of β and γ identified within the literature, sediment deliveries to valley floors across the two catchments are similar for all but one of five sediment transport capacity relationships. The patterns are dominated by the trend in slope through each catchment. The sensitivity to hillslope hydrology of predicted sediment delivery patterns is strong in the catchment with systematic variation in unit hillslope area, and weak in the catchment for which there are no systematic trends in unit hillslope area. We believe there is less experimental evidence to restrict choice of hillslope hydrology parameters than there is for sediment transport capacity. Copyright © 2001 John Wiley & Sons, Ltd.     

Runoff and sediment generation on bench‐terraced hillsides: measurements and up‐scaling of a field‐based model

Despite widespread bench‐terracing, stream sediment yields from agricultural hillsides in upland West Java remain high. We studied the causes of this lack of effect by combining measurements at different spatial scales using an erosion process model. Event runoff and sediment yield from two 4‐ha terraced hillside subcatchments were measured and field surveys of land use, bench‐terrace geometry and storage of sediment in the drainage network were conducted for two consecutive years. Runoff was 3·0–3·9% of rainfall and sediment yield was 11–30 t ha⁻¹ yr⁻¹ for different years, subcatchments and calculation techniques. Sediment storage changes in the subcatchment drainage network were less than 2 t ha⁻¹, whereas an additional 0·3–1·5 t ha⁻¹ was stored in the gully between the subcatchment flumes and the main stream. This suggests mean annual sediment delivery ratios of 86–125%, or 80–104% if this additional storage is included. The Terrace Erosion and Sediment Transport (TEST) model developed and validated for the studied environment was parameterized using erosion plot studies, land use surveys and digital terrain analysis to simulate runoff and sediment generation on the terraced hillsides. This resulted in over‐estimates of runoff and under‐estimates of runoff sediment concentration. Relatively poor model performance was attributed to sample bias in the six erosion plots used for model calibration and unaccounted covariance between important terrain attributes such as slope, infiltration capacity, soil conservation works and vegetation cover. Copyright © 2005 John Wiley & Sons, Ltd.     

Computing spatially distributed sediment delivery ratios: inferring functional sediment connectivity from repeat high‐resolution digital elevation models

High‐resolution digital elevation models (DEMs) from repeat LiDAR (light detection and ranging) or SfM (structure from motion) surveys have become an important tool in process geomorphology. The spatial pattern of negative and positive changes of surface elevation on raster DEMs of difference (DoD) can be interpreted in terms of geomorphic processes, and has been used for morphological budgeting. We show how the application of flow routing algorithms and flow accumulation opens new opportunities for the analysis of DoD. By accumulating the values of the DoD along downslope flowpaths delineated on a DEM, these algorithms lend themselves to computing the net balance, i.e. sediment yield (SY), for the contributing area of each cell. Doing the same for the negative subset of the DoD yields a minimum estimate of erosion (E) within the contributing area. The division of SY by E yields (a maximum estimate of) the sediment delivery ratio (SDR), that is the proportion of material eroded within the contributing area of each cell that has been exported from that area. The resulting SDR is a spatially distributed measure of functional sediment connectivity. In this letter, we develop the computationally simple approach by means of an example DoD from a lateral moraine section in the Upper Kaunertal Valley, Austrian Central Alps. We also discuss advantages, assumptions and limitations, and outline potential applications to connectivity research using field‐, laboratory‐, and model‐based DoD. Copyright © 2018 John Wiley & Sons, Ltd.     

Using sediment travel distance to estimate medium‐term erosion rates: a 16‐year record

Fine grained (80 µm) magnetite was introduced onto a semi‐arid grassland hillslope in 1992, as part of a set of rainfall‐simulation experiments. Using measurements of magnetic susceptibility, the median distance travelled by these magnetite grains during subsequent natural runoff events in the 16‐year period up to 2008 was estimated. Coupling this estimate to direct measurements of sediment flux obtained during the rainfall‐simulation experiments has enabled estimation of the erosion rate over this period. The estimated average erosion rate of between 2·61 × 10^?2 and 4·36 × 10^?2 kg m^?1 year^?1, is equivalent to a rate of ground lowering between 0·020 and 0·033 mm year^?1. This estimate is consistent with (in the sense of being less than) an estimate of total sediment detachment over the same period. The rate of erosion measured using this travel‐distance approach is an order of magnitude less that obtained from a study based on ¹³⁷Cs in a nearby catchment, and compatible with the longevity of continents. Copyright © 2010 John Wiley & Sons, Ltd.     

Meso‐scale catchment sediment budgets: combining field surveys and modeling in the Dragonja catchment,southwest Slovenia

In this paper, we present a methodology to construct a sediment budget for meso‐scale catchments. We combine extensive field surveys and expert knowledge of the catchment with a sediment delivery model. The meso‐scale Mediterranean drainage basin of the Dragonja (91 km²), southwest Slovenia, was chosen as case study area. During the field surveys, sheet wash was observed on sloping agricultural fields during numerous rainfall events, which was found to be the main source of sediment. With the sediment yield model WATEM/SEDEM the estimated net erosion on the hillslopes 4·1 t ha^–1 y^–1 (91% of inputs). The second source, bank erosion (4·2%; 0·25 t ha^–1 y^–1) was monitored during several years with erosion pins and photogrammetric techniques. The last source, channel incision, was derived from geomorphological mapping and lichenomery and provided 3·8% (0·17 t ha^–1 y^–1) of the sediment input. The river transports its suspended sediment mainly during high‐flow events (sampled with automated water samplers). About 27% (1·2 t ha^–1 y^–1) of the sediment delivered to the channel is deposited on floodplains and low terraces downstream (estimated with geomorphological mapping, coring and cesium‐137 measurements). The sediment transported as bedload disintegrates during transport to the outlet due to the softness of the bedrock material. As a result, the river carries no bedload when it reaches the sea. The results imply a build‐up of sediment in the valleys catchment. However, extreme flood events may flush large amounts of sediment stored in the lower parts of the system. Geomorphological evidence exists in the catchment that such high magnitude, low frequency events have happened in the past. Copyright © 2009 John Wiley & Sons, Ltd.     

三峡水库区间来沙量估算及水库排沙效果分析

三峡水库是控制和开发长江的重要工程,具有巨大的防洪、发电、航运等综合效益.水库排沙比是反映水库拦截泥沙程度的重要指标.针对目前三峡水库排沙比计算没有考虑三峡水库区间来沙,结果偏大的问题,依据三峡水库区间水文站2003-2016年实测水文资料,采用水文学法估算了三峡水库区间月、年入库输沙量,分析了三峡水库排沙效果及主要影响因素.结果表明：1）三峡水库区间年均来沙量约1775×10⁴ t,占总入库沙量的10.3%,近年来沙占比有所增大,其中2013-2016年来沙量占总入库沙量的26.9%;2）三峡水库年均排沙比为21.6%,其中围堰蓄水期为34.1%,初期蓄水期为17.0%,试验性蓄水期为15.4%;3）三峡水库采用"蓄清排浑"方式运行,主汛期7-9月的排沙比一般大于枯季,但2013年以后,出现了主汛期排沙比小于枯季的现象;4）入库细颗粒泥沙的排沙比大于粗颗粒泥沙,其中粒径d ≤ 0.062 mm的细颗粒泥沙排沙比为23.4%,0.062 mm0.125 mm的中粗沙排沙比分别为5.5%和11.1%;5）三峡水库排沙比汛期主要受V/Q影响,枯季主要受入库含沙量的影响;当V/Q约为170×10⁴ s时,水库排沙效果最差;分别建立了汛期和枯季排沙比经验计算式.本文的研究成果可为三峡水库水沙优化调度等提供参考.     

A new approach for linking event‐based upland sediment sources to downstream suspended sediment transport

In this study, we proposed a new approach for linking event sediment sources to downstream sediment transport in a watershed in central New York. This approach is based on a new concept of spatial scale, sub‐watershed area (SWA), defined as a sub‐watershed within which all eroded soils are transported out without deposition during a hydrological event. Using (rainfall) event data collected between July and November, 2007 from several SWAs of the studied watershed, we developed an empirical equation that has one independent variable, mean SWA slope. This equation was then used to determine event‐averaged unit soil erosion rate, Q_S/A, (in kg/km²/hr) for all SWAs in the studied watershed and calculate event‐averaged gross erosion E_ea (in kg/hr). The event gross erosion E_t (in kilograms) was subsequently computed as the product of E_ea and the mean event duration, T (in hours) determined using event hydrographs at the outlet of the studied watershed. Next, we developed two linear sediment rating curves (SRCs) for small and big events based on the event data obtained at the watershed outlet. These SRCs, together with T, allowed us to determine event sediment yield SY_e (in kilograms) for all events during the study period. By comparing E_t with SY_e, developing empirical equations (i) between E_t and SY_e and (ii) for event sediment delivery ratio, respectively, we revealed the event dynamic processes connecting sediment sources and downstream sediment transport. During small events, sediment transport in streams was at capacity and dominated by the deposition process, whereas during big events, it was below capacity and controlled by the erosion process. The key of applying this approach to other watersheds is establishing their empirical equations for Q_S/A and appropriately determining their numbers of SWAs. Copyright © 2011 John Wiley & Sons, Ltd.     

The impact of pro‐glacial lakes on downstream sediment delivery in Norway

Glacier recessions caused by climate change may uncover pro‐glacial lakes that form important sedimentation basins regulating the downstream sediment delivery. The impact of modern pro‐glacial lakes on fluvial sediment transport from three different Norwegian glaciers: Nigardsbreen, Engabreen and Tunsbergdalsbreen, and their long‐term development has been studied. All of these lakes developed in modern times in overdeepened bedrock basins. The recession of Nigardsbreen uncovered a 1.8 km long and on average 15 m deep pro‐glacial lake basin during 1937 to 1968. Since then the glacier front has been situated entirely on land, and the sediment input and output of the lake has been measured. The suspended sediment transport into and out of the lake averaged 11 730 t yr^?1 and 2340 t yr^?1 respectively. Thus, 20% remained in suspension at the outlet. The measured mean annual bedload supplied to the lake was 11 800 t yr^?1, giving a total transport of 23 530 t yr^?1 which corresponds to a specific sediment yield of 561 t km^?2 yr^?1. A 1.9 km long and up to 90 m deep pro‐glacial lake basin downstream from Engabreen glacier was uncovered during 1890 to 1944. The average suspended sediment load delivered from the glacier during the years 1970–1981 amounted to 12 375 t yr^?1and the transport out of the lake was 2021 t yr^?1, giving an average of 16% remaining in suspension. The mean annual bedload was 8000 t yr^?1, thus the total transport was 20 375 t yr^?1, giving a specific sediment yield of 566 t km^?2 yr^?1. For Tunsbergdalsbreen glacier, measurements in the early 1970s indicated that the suspended sediment transport was on average 44 000 t yr^?1. From 1987 to 1993 the recession of the glacier uncovered a small pro‐glacial lake, 0.3 km long and around 9 m deep. Downstream from this, the suspended sediment load measured in 2009 was 28 000 t yr^?1, indicating that as much as 64% remained in suspension. Flow velocity, grain size of sediment, and morphology of the lake are important factors controlling the sedimentation rate in the pro‐glacial lakes. A survey of the sub‐glacial morphology of Tunsbergdalsbreen revealed that there are several overdeepened basins beneath the glacier. The largest is 4 km long and 100 m deep. When the glacier melts back they will become lakes and act as sedimentation basins. Despite an expected increase in sediment yield from the glacier, little sediment will pass these lakes and downstream sediment delivery will be reduced markedly. Beneath Nigardsbreen there was only a small depression that may form a lake and the sediment delivery will not be significantly affected. © 2014 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Rainstorm-generated sediment yield model based on soil moisture proxies (SMP)

This study develops improved Soil Moisture Proxies (SMP) based suspended sediment yield (SMPSY) models corresponding to three antecedent moisture conditions (AMCs) (i.e., AMC-I-AMC-III) by coupling the improved initial abstraction (I_a-λ) model, the SMA procedure and the SMP concept for modelling the rainfall generated suspended sediment yield. The SMPSY models specifically incorporate a watershed storage index (S) model to accentuate the transformation from storm to storm and to avoid the sudden jumps in sediment yield computation. The workability of the SMPSY models is tested using a large dataset of rainfall and sediment yield (98 storm events) from twelve small watersheds and a comparison has been made with the existing MSY model. The goodness-of-fit (GOF) statistics is evaluated in terms of the Nash Sutcliffe efficiency (NSE), and error indices, i.e., root mean square error (RMSE), normalized root mean square error (nRMSE), standard error (SE), mean absolute error (MAE), and RMSE-observations standard deviation ratio (RSR). The NSE values vary from 74.31% to 96.57% and from 75.21% to 91.78%, respectively for the SPMSY and MSY model. The NSE statistics indicate that the SMPSY model has lower uncertainty in simulating sediment yield as compared to the MSY model. The error indices are lower for the SMPSY model than the MSY model for most of the watersheds. These results show that the SMPSY model has less uncertainty and performs better than the MSY model. A sensitivity analysis of the SMPSY model shows that the parameter β is most sensitive followed by parameter S, α and A. Overall, the results show that the characterization of soil moisture variability in terms of SMPs and incorporation of improved delivery ratio and runoff coefficient relationship improves the simulation of the erosion and sediment yield generation process.     

Evaluation of a dynamic multi‐class sediment transport model in a catchment under soil‐conservation agriculture

Soil erosion models are essential tools for the successful implementation of effective and adapted soil conservation measures on agricultural land. Therefore, models are needed that predict sediment delivery and quality, give a good spatial representation of erosion and deposition and allow us to account for various soil conservation measures. Here, we evaluate how well a modified version of the spatially distributed multi‐class sediment transport model (MCST) simulates the effectiveness of control measures for different event sizes. We use 8 year runoff and sediment delivery data from two small agricultural watersheds (0·7 and 3·7 ha) under optimized soil conservation. The modified MCST model successfully simulates surface runoff and sediment delivery from both watersheds; one of which was dominated by sheet and the other was partly affected by rill erosion. Moreover, first results of modelling enrichment of clay in sediment delivery are promising, showing the potential of MCST to model sediment enrichment and nutrient transport. In general, our results and those of an earlier modelling exercise in the Belgian Loess Belt indicate the potential of the MCST model to evaluate soil erosion and deposition under different agricultural land uses. As the model explicitly takes into account the dominant effects of soil‐conservation agriculture, it should be successfully applicable for soil‐conservation planning/evaluation in other environments. Copyright © 2008 John Wiley & Sons, Ltd.     

Bioturbation on a south‐east Australian hillslope: estimating contributions to soil flux

The contribution of bioturbation to downslope soil transport is significant in many situations, particularly in the context of soil formation, erosion and creep. This study explored the direct flux of soil caused by Aphaenogaster ant mounding, vertebrate scraping and tree‐throw on a wildfire‐affected hillslope in south‐east Australia. This included the development of methods previously applied to Californian gopher bioturbation, and an evaluation of methods for estimating the volume of soil displaced by tree‐throw events. All three bioturbation types resulted in a net downslope flux, but any influence of hillslope angle on flux rates appeared to be overshadowed by environmental controls over the spatial extent of bioturbation. As a result, the highest flux rates occurred on the footslope and lower slope. The overall contribution of vertebrate scraping (57.0 ± 89.4 g m^?1 yr^?1) exceeded that of ant mounding (36.4 ± 66.0 g m^?1 yr^?1), although mean rates were subject to considerable uncertainty. Tree‐throw events, which individually cause major disturbance, were limited in their importance by their scarcity relative to faunalturbation. However, tree‐throw might be the dominant mechanism of biotic soil flux on the mid‐slope provided that it occurs at a frequency of at least 2–3 events ha^?1 yr^?1. Although direct biotic soil flux appears to be geomorphologically significant on this hillslope, such transport processes are probably subordinate to other impacts of bioturbation at this site such as the enhancement of infiltration following wildfire. Copyright © 2011 John Wiley & Sons, Ltd.     

Tracking sediment redistribution in a small watershed: implications for agro‐landscape evolution

A new, multi‐tracer method is used to track erosion, translocation, and redeposition of sediment in a small watershed, thus allowing for the ?rst time a complete, spatially distributed, sediment balance to be made as a function of landscape position. A 0·68 ha watershed near Coshocton, Ohio, USA was divided into six morphological units, each tagged with one of six rare earth element oxides. Sediment translocation was evaluated by collecting run‐off and by spatially sampling the soil surface. Average measured erosion rate was 6·1 t ha^?1, but varied between 40·4 t ha^?1 loss from the lower channels to 24·1 t ha^?1 gain on the toeslope. With this technique it was possible for the ?rst time to itemize the sediment budget for landscape elements into three components: (1) the soil from the element that left the watershed with run‐off; (2) soil from the element that was redeposited on lower positions, with the spatial distribution of that deposition; and (3) soil originating from the upper positions and deposited on the element, with quanti?cation of relative source areas. The results are incongruous with the current morphology of the watershed, suggesting that diffusion‐type erosion must also play a major role in de?ning the evolution of this landscape. Copyright © 2004 John Wiley & Sons, Ltd.     

Distributed empirical algorithms to estimate catchment scale sediment connectivity and yield in a subtropical region

The intensity of soil loss and sediment delivery, representing hydrologic and geomorphic processes within a catchment, accelerates with rapid changes in land cover and rainfall events. An underlying component of sustainable management of water resources is an understanding of spatial and temporal variability and the adverse influences of regional parameters involved in generating sediment following widespread changes in land cover. A calibrated algorithm of soil loss coupled with a sediment delivery ratio (SDR) was applied in raster data layers to improve the capability of a combined model to estimate annual variability in sediment yields related to changes in vegetation cover identified by analyses of SPOT imagery. Four catchments in Kangaroo River State forest were assessed for annual changes in sediment yields. Two catchments were selectively logged in 2007, while the two other sites remained undisturbed. Results of SDR estimates indicated that only a small proportion of total eroded sediment from hillslopes is transported to catchment outlets. Larger SDR values were estimated in regions close to catchment outlets, and the SDR reduced sharply on hillslopes further than 200–300 m from these areas. Estimated sediment yield increased by up to 30% two years after land cover change (logging) in 2009 when more storm events were recorded, despite the moderate density of vegetation cover in 2009 having almost recovered to its initial pre‐logging (2005) condition. Rainfall had the most significant influence on streamflow and sediment delivery in all catchments, with steeply sloping areas contributing large amounts of sediment during moderate and high rainfall years in 2007 and 2009. It is concluded that the current scenario of single‐tree selection logging utilized in the study area is an acceptable and environmentally sound land management strategy for preservation of soil and water resources. Copyright © 2013 John Wiley & Sons, Ltd.     

Flume‐ and field‐based evaluation of a time‐integrated suspended sediment sampler for the analysis of sediment properties

Suspended sediment has been identified as a vector for nutrient and contaminant transport in the fluvial environment. A time‐integrated sampler (the Phillips sampler), which emerged over a decade ago as a cost‐effective tool for in situ suspended sediment collection, is increasingly being used to collect samples for the analysis of sediment properties such as particle size composition, and nutrient and contaminant concentrations. This study evaluates the sampler under both flume and field conditions for efficiency in the mass and grain size of the suspended sediment collected. The sampler was tested in a flume using both kaolinite and sediment samples (sieved to < 180 µm) collected from the Quesnel River, British Columbia, Canada. In the kaolinite trails, the sampler preferentially collected coarser grain sizes compared to the original sediment, probably due to finer sediment remaining in suspension and therefore passing through the sampler, and also possibly due to flocculation of the kaolinite upon introduction to the flume. Conversely, the sampler collected river sediment that was finer than the original sediment, probably due to some settling of coarser sediment observed at the bottom of the flume. Once allowance was made for these operational issues associated with the flume, maximum sediment mass efficiency for kaolinite and river sediment was 43% and 87%, respectively. Sediment collected by the time‐integrated sampler during field deployment and adjacent channel bed sediment were also compared. The sampler collected sediment with a representative grain size distribution. However, there were differences in the geochemical (arsenic and selenium) concentrations of channel bed sediment and sediment collected by the Phillips sampler which may be a function of differences in the behavior of geochemical elements associated with the two types of sediment. This work suggests that further research is needed to evaluate the role of the Phillips sampler in collecting sediment for contaminant and nutrient analysis. Copyright © 2014 John Wiley & Sons, Ltd.     

Flood‐based analysis of high‐magnitude sediment transport using a non‐parametric method

Upland erosion and the resulting reservoir siltation is a serious issue in the Isábena catchment (445 km² Central Spanish Pyrenees). During a three‐month period, water and sediment fluxes have been monitored at the catchment outlet (Capella), two adjacent subcatchments (Villacarli, 41 km²; Cabecera, 145 km²) and the elementary badland catchment Torrelaribera (8 ha). This paper presents the results of the monitoring, a method for the calculation of a sedigraph from intermittent measurements and the derived sediment yields at the monitored locations. The observed suspended sediment concentrations (SSCs) demonstrate the role of badlands as sediment sources: SSCs of up to 280 g l^?1 were encountered for Villacarli, which includes large badland areas. SSCs at the Cabecera catchment, with great areas of woodland, barely exceeded 30 g l^?1. SSCs directly at the sediment source (Torrelaribera) were comparable to those at Villacarli, suggesting a close connection within this subcatchment. At Capella, SSCs of up to 99 g l^?1 were observed. For all sites, SSC displayed only a loose correlation with discharge, inhibiting the application of a simple sediment rating curve. Instead, ancillary variables acting as driving forces or proxies for the processes (rainfall energy, cumulative discharge, rising/falling limb data) were included in a quantile regression forest model to explain the variability in SSC. The variables with most predictive power vary between the sites, suggesting the predominance of different processes. The subsequent flood‐based calculation of sediment yields attests high specific sediment yields for Torrelaribera and Villacarli (6277 and 1971 t km^?2) and medium to high yields for Cabecera and Capella (139 and 410 t km^?2) during the observation period. In all catchments, most of the sediment was exported during intense storms of late summer. Later flood events yield successively less sediment. Relating upland sediment production to yield at the outlet suggests considerable effects of sediment storage within the river channel. Copyright © 2008 John Wiley & Sons, Ltd.     

Estimation of suspended sediment loads and delivery in an incised upland headwater catchment,south‐eastern Australia

Historically upland headwater catchments in south‐eastern Australia have undergone extensive gully erosion that has removed large amounts of sediment to lowlands. Recent research suggests these upland areas may continue to dominate fine sediment loads in lowland rivers. Improved understanding of sediment transfer through upland headwater catchments may have implications for interpreting downstream sediment supply. In this study a nested catchment design was utilized to examine suspended sediment yields and delivery from a small tributary sub‐catchment (1·64 km²) to the study catchment outlet (53·5 km²). Monitoring of suspended sediment concentration and discharge was undertaken for a period of nearly two years and used to estimate suspended sediment loads. Estimated total suspended sediment exports over the period of monitoring were 24·16 t from the sub‐catchment and 550·3 t from the catchment, which are generally less than previous reported small catchment yields in south‐eastern Australia. The extent of sediment delivery was examined using between‐site ratios of specific sediment yield per unit area and incised channel length. Sediment delivery was high under average rainfall conditions, but seasonally dependent. Both suspended sediment yields and the extent of delivery peaked over spring months, supplemented by remobilization of sediment stored during summer months in the main catchment channel. The findings of this study suggest much of the suspended sediment exported from small incised upland sub‐catchments (1–2 km²) may be delivered to downstream reaches under average rainfall conditions, which, in conjunction with the findings of previous research supports the potential importance of contributions from these areas to suspended sediment loads in lowland rivers during high flow periods. Copyright © 2007 John Wiley & Sons, Ltd.     

Gypsies in the palace: experimentalist's view on the use of 3‐D physics‐based simulation of hillslope hydrological response

As a fundamental unit of the landscape, hillslopes are studied for their retention and release of water and nutrients across a wide range of ecosystems. The understanding of these near‐surface processes is relevant to issues of runoff generation, groundwater–surface water interactions, catchment export of nutrients, dissolved organic carbon, contaminants (e.g. mercury) and ultimately surface water health. We develop a 3‐D physics‐based representation of the Panola Mountain Research Watershed experimental hillslope using the TOUGH2 sub‐surface flow and transport simulator. A recent investigation of sub‐surface flow within this experimental hillslope has generated important knowledge of threshold rainfall‐runoff response and its relation to patterns of transient water table development. This work has identified components of the 3‐D sub‐surface, such as bedrock topography, that contribute to changing connectivity in saturated zones and the generation of sub‐surface stormflow. Here, we test the ability of a 3‐D hillslope model (both calibrated and uncalibrated) to simulate forested hillslope rainfall‐runoff response and internal transient sub‐surface stormflow dynamics. We also provide a transparent illustration of physics‐based model development, issues of parameterization, examples of model rejection and usefulness of data types (e.g. runoff, mean soil moisture and transient water table depth) to the model enterprise. Our simulations show the inability of an uncalibrated model based on laboratory and field characterization of soil properties and topography to successfully simulate the integrated hydrological response or the distributed water table within the soil profile. Although not an uncommon result, the failure of the field‐based characterized model to represent system behaviour is an important challenge that continues to vex scientists at many scales. We focus our attention particularly on examining the influence of bedrock permeability, soil anisotropy and drainable porosity on the development of patterns of transient groundwater and sub‐surface flow. Internal dynamics of transient water table development prove to be essential in determining appropriate model parameterization. Copyright © 2010 John Wiley & Sons, Ltd.