Simulated effect of a forest road on near‐surface hydrologic response: redux

In the work reported here the comprehensive physics‐based Integrated Hydrology Model (InHM) was employed to conduct both three‐ and two‐dimensional (3D and 2D) hydrologic‐response simulations for the small upland catchment known as C3 (located within the H. J. Andrews Experimental Forest in Oregon). Results from the 3D simulations for the steep unchannelled C3 (i) identify subsurface stormflow as the dominant hydrologic‐response mechanism and (ii) show the effect of the down‐gradient forest road on both the surface and subsurface flow systems. Comparison of the 3D results with the 2D results clearly illustrates the importance of convergent subsurface flow (e.g. greater pore‐water pressures in the hollow of the catchment for the 3D scenario). A simple infinite‐slope model, driven by subsurface pore‐water pressures generated from the 3D and 2D hydrologic‐response simulations, was employed to estimate slope stability along the long‐profile of the C3 hollow axis. As expected, the likelihood of slope failure is underestimated for the lower pore pressures from the 2D hydrologic‐response simulation compared, in a relative sense, to the higher pore pressures from the 3D hydrologic response simulation. The effort reported herein provides a firm quantitative foundation for generalizing the effects that forest roads can have on near‐surface hydrologic response and slope stability at the catchment scale. Copyright © 2006 John Wiley & Sons, Ltd.     

Rapid simulated hydrologic response within the variably saturated near surface

Column and field experiments have shown that the hydrologic response to increases in rainfall rates can be more rapid than expected from simple estimates. Physics‐based hydrologic response simulation, with the Integrated Hydrology Model (InHM), is used here to investigate rapid hydrologic response, within the variably saturated near surface, to temporal variations in applied flux at the surface boundary. The factors controlling the speed of wetting front propagation are discussed within the Darcy–Buckingham conceptual framework, including kinematic wave approximations. The Coos Bay boundary‐value problem is employed to examine simulated discharge, pressure head, and saturation responses to a large increase in applied surface flux. The results presented here suggest that physics‐based simulations are capable of representing rapid hydrologic response within the variably saturated near surface. The new InHM simulations indicate that the temporal discretization and measurement precision needed to capture the rapid subsurface response to a spike increase in surface flux, necessary for both data‐based analyses and evaluation of physics‐based models, are smaller than the capabilities of the instrumentation deployed at the Coos Bay experimental catchment. Copyright © 2007 John Wiley & Sons, Ltd.     

Understanding the water balance paradox in the Athabasca River Basin,Canada

This study demonstrates the importance of the including and appropriately parameterizing peatlands and forestlands for basin‐scale integrated surface–subsurface models in the northern boreal forest, with particular emphasis on the Athabasca River Basin (ARB). With a long‐term water balance approach to the ARB, we investigate reasons why downstream mean annual stream flow rates are consistently higher than upstream, despite the subhumid water deficit conditions in the downstream regimes. A high‐resolution 3D variably saturated subsurface and surface water flow and evapotranspiration model of the ARB is constructed based on the bedrock and surficial geology and the spatial distribution of peatlands and their corresponding eco‐regions. Historical climate data were used to drive the model for calibration against 40‐year long‐term average surface flow and groundwater observations during the historic instrumental period. The simulation results demonstrate that at the basin‐scale, peatlands and forestlands can have a strong influence on the surface–subsurface hydrologic systems. In particular, peatlands in the midstream and downstream regimes of the ARB increase the water availability to the surface–subsurface water systems by reducing water loss through evapotranspiration. Based on the comparison of forestland evapotranspiration between observation and simulation, the overall spatial average evapotranspiration in downstream forestlands is larger than that in peatlands and thus the water contribution to the stream flow in downstream areas is relatively minor. Therefore, appropriate representation of peatlands and forestlands within the basin‐scale hydrologic model is critical to reproduce the water balance of the ARB.     

Effects of hillslope topography on hydrological responses in a weathered granite mountain,Japan: comparison of the runoff response between the valley‐head and the side slope

To evaluate the effects of hillslope topography on storm runoff in a weathered granite mountain, discharge rate, soil pore water pressures, and water chemistry were observed on two types of hillslope: a valley‐head (a concave hillslope) and a side slope (a planar hillslope). Hydrological responses on the valley‐head and side slope reflected their respective topographic characteristics and varied with the rainfall magnitude. During small rainfall events (<35 mm), runoff from the side slope occurred rapidly relative to the valley‐head. The valley‐head showed little response in storm runoff. As rainfall amounts increased (35–60 mm), the valley‐head yielded a higher flow relative to the side slope. For large rainfall events (>60 mm), runoff from both hillslopes increased with rainfall, although that from the valley‐head was larger than that from the side slope. The differences in the runoff responses were caused by differences in the roles of lower‐slope soils and the convergence of the hillslope. During small rainfall events, the side slope could store little water; in contrast, all rainwater could be stored in the soils at the valley‐head hollow. As the amount of rainfall increased, the subsurface saturated area of the valley‐head extended from the bottom to the upper portion of the slope, with the contributions of transient groundwater via lateral preferential flowpaths due to the high concentration of subsurface water. Conversely, saturated subsurface flow did not contribute to runoff responses, and the subsurface saturated area at the side slope did not extend to the upper slope for the same storm size. During large rainfall events, expansion of the subsurface saturated area was observed in both hillslopes. Thus, differences in the concentration of subsurface water, reflecting hillslope topography, may create differences in the extension of the subsurface saturated area, as well as variability in runoff responses. Copyright © 2007 John Wiley & Sons, Ltd.     

Physically‐based modelling of double‐peak discharge responses at Slapton Wood catchment

Heavy winter rainfall produces double‐peak hydrographs at the Slapton Wood catchment, Devon, UK. The first peak is saturation‐excess overland flow in the hillslope hollows and the second (i.e. the delayed peak) is subsurface stormflow. The physically‐based spatially‐distributed model SHETRAN is used to try to improve the understanding of the processes that cause the double peaks. A three‐stage (multi‐scale) approach to calibration is used: (1) water balance validation for vertical one‐dimensional flow at arable, grassland and woodland plots; (2) two‐dimensional flow for cross‐sections cutting across the stream valley; and (3) three‐dimensional flow in the full catchment. The main data are for rainfall, stream discharge, evaporation, soil water potential and phreatic surface level. At each scale there was successful comparison with measured responses, using as far as possible parameter values from measurements. There was some calibration but all calibrated values at one scale were used at a larger scale. A large proportion of the subsurface runoff enters the stream from three dry valleys (hillslope hollows), and previous studies have suggested convergence of the water in the three large hollows as being the major mechanism for the production of the delayed peaks. The SHETRAN modelling suggests that the hillslopes that drain directly into the stream are also involved in producing the delayed discharges. The model shows how in the summer most of the catchment is hydraulically disconnected from the stream. In the autumn the catchment eventually ‘wets up’ and shallow subsurface flows are produced, with water deflected laterally along the soil‐bedrock interface producing the delayed peak in the stream hydrograph. Copyright © 2007 John Wiley & Sons, Ltd.     

Runoff generation from logged and burnt convergent hillslopes: rainfall simulation and modelling

This paper reports results from field experiments and hydrological modelling on the dynamics of runoff generation in highly convergent parts of the landscape in a logged and burnt eucalypt forest in south‐eastern Victoria, Australia. Large‐scale rainfall simulation experiments were conducted to explore runoff generating mechanisms from harvested areas, and to assess the effectiveness of standard water quality protective measures, here a disturbed filter strip, in preventing accession of sediment to near‐stream areas. We then examined the likely effects of varying antecedent moisture conditions on surface and subsurface runoff generating mechanisms. Very small volumes of surface runoff were generated only at very high rainfall intensity rates that exceeded a 100 year recurrence interval event during the simulated experiments. There was little or no identifiable impact of either compaction from logging operations or fire‐induced hydrophobicity on surface infiltration or generation of surface runoff. Measured soil hydraulic properties and soil depths explained the paucity of surface runoff, and the dominance of subsurface storm flow as the prime runoff generating mechanism. Deep lateral subsurface flow was observed from the cut‐face of a fire access track and into a streamhead downslope of the experimental plots. Water balance modelling using Topog_Dynamic indicated the conditions under which saturated overland flow in this environment could be generated are rare, but that care should be taken in siting of roads and tracks in lower parts of convergent landscapes. Copyright © 2004 John Wiley & Sons, Ltd.     

Characterizing hydrological processes in a semiarid rangeland watershed: A hydrogeophysical approach

The complex ecohydrological processes of rangelands can be studied through the framework of ecological sites (ESs) or hillslope‐scale soil–vegetation complexes. High‐quality hydrologic field investigations are needed to quantitatively link ES characteristics to hydrologic function. Geophysical tools are useful in this context because they provide valuable information about the subsurface at appropriate spatial scales. We conducted 20 field experiments in which we deployed time‐lapse electrical resistivity tomography (ERT), variable intensity rainfall simulation, ground‐penetrating radar (GPR), and seismic refraction, on hillslope plots at five different ESs within the Upper Crow Creek Watershed in south‐east Wyoming. Surface runoff was measured using a precalibrated flume. Infiltration data from the rainfall simulations, coupled with site‐specific resistivity–water content relationships and ERT datasets, were used to spatially and temporally track the progression of the wetting front. First‐order constraints on subsurface structure were made at each ES using the geophysical methods. Sites ranged from infiltrating 100% of applied rainfall to infiltrating less than 60%. Analysis of covariance results indicated significant differences in the rate of wetting front progression, ranging from 0.346 m min^?1/2 for sites with a subsurface dominated by saprolitic material to 0.156 m min^?1/2 for sites with a well‐developed soil profile. There was broad agreement in subsurface structure between the geophysical methods with GPR typically providing the most detail. Joint interpretation of the geophysics showed that subsurface features such as soil layer thickness and the location of subsurface obstructions such as granite corestones and material boundaries had a large effect on the rate of infiltration and subsurface flow processes. These features identified through the geophysics varied significantly by ES. By linking surface hydrologic information from the rainfall simulations with subsurface information provided by the geophysics, we can characterize the ES‐specific hydrologic response. Both surface and subsurface flow processes differed among sites and are directly linked to measured characteristics.     

Impact of road‐generated storm runoff on a small catchment response

The impact of road‐generated runoff on the hydrological response of a zero‐order basin was monitored for a sequence of 24 storm events. The study was conducted in a zero‐order basin (C1; 0·5ha) with an unpaved mountain road; an adjacent unroaded zero‐order basin (C2; 0·2 ha) with similar topography and lithology was used to evaluate the hydrological behaviour of the affected zero‐order basin prior to construction of the road. The impact of the road at the zero‐order basin scale was highly dependent on the antecedent soil‐moisture conditions, total storm precipitation, and to some extent rainfall intensity. At the beginning of the monitoring period, during dry antecedent conditions, road runoff contributed 50% of the total runoff and 70% of the peak flow from the affected catchment (C1). The response from the unroaded catchment was almost insignificant during dry antecedent conditions. As soil moisture increased, the road exerted less influence on the total runoff from the roaded catchment. For very wet conditions, the influence of road‐generated runoff on total outflow from the roaded catchment diminished to only 5·4%. Both catchments, roaded and unroaded, produced equivalent amount of outflow during very wet antecedent conditions on a unit area basis. The lag time between the rainfall and runoff peaks observed in the unroaded catchment during the monitoring period ranged from 0 to 4 h depending on the amount of precipitation and antecedent conditions, owing mainly to much slower subsurface flow pathways in the unroaded zero‐order basin. In contrast, the lag time in the roaded zero‐order basin was virtually nil during all storms. Copyright © 2009 John Wiley & Sons, Ltd.     

Impacts of forest conversion and agriculture practices on water pathways in Southern Brazil

Land‐use/cover change (LUCC), and more specifically deforestation and multidecadal agriculture, is one of the various controlling factors of water fluxes at the hillslope or catchment scale. We investigated the impact of LUCC on water pathways and stream stormflow generation processes in a subtropical region in southern Brazil. We monitored, sampled and analysed stream water, pore water, subsurface water, and rainwater for dissolved silicon concentration (DSi) and ¹⁸O/¹⁶O (δ¹⁸O) signature to identify contributing sources to the streamflow under forest and under agriculture. Both forested and agricultural catchments were highly responsive to rainfall events in terms of discharge and shallow groundwater level. DSi versus δ¹⁸O scatter plots indicated that for both land‐use types, two run‐off components contributed to the stream discharge. The presence of a dense macropore network, combined with the presence of a compact and impeding B‐horizon, led to rapid subsurface flow in the forested catchment. In the agricultural catchment, the rapid response to rainfall was mostly due to surface run‐off. A 2‐component isotopic hydrograph separation indicated a larger contribution of rainfall water to run‐off during rainfall event in the agricultural catchments. We attributed this higher contribution to a decrease in topsoil hydraulic conductivity associated with agricultural practices. The chemical signature of the old water component in the forested catchment was very similar to that of the shallow groundwater and the pore soil water: It is therefore likely that the shallow groundwater was the main source of old water. This is not the case in the agricultural catchments where the old water component had a much higher DSi concentration than the shallow groundwater and the soil pore water. As the agricultural catchments were larger, this may to some extent simply be a scale effect. However, the higher water yields under agriculture and the high DSi concentration observed in the old water under agriculture suggest a significant contribution of deep groundwater to catchment run‐off under agriculture, suggesting that LUCC may have significant effects on weathering rates and patterns.     

A physically based distributed subsurface–surface flow dynamics model for forested mountainous catchments

This study was designed to develop a physically based hydrological model to describe the hydrological processes within forested mountainous river basins. The model describes the relationships between hydrological fluxes and catchment characteristics that are influenced by topography and land cover. Hydrological processes representative of temperate basins in steep terrain that are incorporated in the model include intercepted rainfall, evaporation, transpiration, infiltration into macropores, partitioning between preferential flow and soil matrix flow, percolation, capillary rise, surface flow (saturation‐excess and return flow), subsurface flow (preferential subsurface flow and baseflow) and spatial water‐table dynamics. The soil–vegetation–atmosphere transfer scheme used was the single‐layer Penman–Monteith model, although a two‐layer model was also provided. The catchment characteristics include topography (elevation, topographic indices), slope and contributing area, where a digital elevation model provided flow direction on the steepest gradient flow path. The hydrological fluxes and catchment characteristics are modelled based on the variable source‐area concept, which defines the dynamics of the watershed response. Flow generated on land for each sub‐basin is routed to the river channel by a kinematic wave model. In the river channel, the combined flows from sub‐basins are routed by the Muskingum–Cunge model to the river outlet; these comprise inputs to the river downstream. The model was applied to the Hikimi river basin in Japan. Spatial decadal values of the normalized difference vegetation index and leaf area index were used for the yearly simulations. Results were satisfactory, as indicated by model efficiency criteria, and analysis showed that the rainfall input is not representative of the orographic lifting induced rainfall in the mountainous Hikimi river basin. Also, a simple representation of the effects of preferential flow within the soil matrix flow has a slight significance for soil moisture status, but is insignificant for river flow estimations. Copyright © 2005 John Wiley & Sons, Ltd.     

Modelling long‐term contaminant migration in a catchment at fine spatial and temporal scales using the UP system

A method has been developed to simulate the long‐term migration of radionuclides in the near‐surface of a river catchment, following their release from a deep underground repository for radioactive waste. Previous (30‐year) simulations, conducted using the SHETRAN physically based modelling system, showed that long‐term (many decades) simulations are required to allow the system to reach steady state. Physically based, distributed models, such as SHETRAN, tend to be too computationally expensive for this task. Traditional lumped catchment‐scale models, on the other hand, do not give sufficiently detailed spatially distributed results. An intermediate approach to modelling has therefore been developed which allows flow and transport processes to be simulated with the spatial resolution normally associated with distributed models, whilst being computationally efficient.The approach involves constructing a lumped model in which the catchment is represented by a number of conceptual water storage compartments. The flow rates to and from these compartments are prescribed by functions that summarize the results from physically based distributed models run for a range of characteristic flow regimes. The physically based models used were, SHETRAN for the subsurface compartments, a particle tracking model for overland flow and an analytical model for channel routing. One important advantage of the method used in constructing the lumped model is that it makes down scaling possible, in the sense that fine‐scale information on the distributed hydrological regime, as simulated by the physically based distributed models, can be inferred from the variables in the lumped model that describe the hydrology at the catchment scale. A 250‐year flow simulation has been run and the down scaling process used to infer a 250‐year time‐series of three‐dimensional velocity fields for the subsurface of the catchment. This series was then used to drive a particle tracking simulation of contaminant migration. The concentration and spatial distribution of contaminants simulated by this model for the first 30 years were in close agreement with SHETRAN results. The remaining 220 years highlighted the fact that some of the most important transport pathways to the surface carry contaminants only very slowly so both the magnitude and spatial distribution of concentration in surface soils are not apparent over the shorter SHETRAN simulations. Copyright © 1999 John Wiley & Sons, Ltd.     

The effects of topography and forest management on water storage in catchments in south‐central Chile

Our work analyses the intra‐annual variability of the volume of water stored in 15 forested headwater catchments from south‐central Chile, aiming at understanding how forest management, hydrology, and climate influence the dynamic components of catchment storage. Thus, we address the following questions: (a) How does the annual water storage vary in catchments located in diverse hydroclimatic conditions and subject to variable forest management? (b) Which natural (i.e., hydrologic regime and physiographic setting) and anthropogenic factors explain the variance in water storage? Results show that the annual catchment storage increases at the beginning of each hydrological year in direct response to increases in rainfall. The maximum water storage ranges from 666 to 1,272 mm in these catchments. The catchments with Pinus or Eucalyptus spp. cover store less water than the catchments with mixed forest species cover. Forest cover (biomass volume, plantation density, and percentage of plantation and age) has the primary control on dynamic storage in all catchments. These results indicate that forest management may alter the catchment water storage.     

Thresholds for run‐off generation in a drained closed depression

Hydrological threshold behaviour has been observed across hillslopes and catchments with varying characteristics. Few studies, however, have evaluated rainfall–run‐off response in areas dominated by agricultural land use and artificial subsurface drainage. Hydrograph analysis was used to identify distinct hydrological events over a 9‐year period and examine rainfall characteristics, dynamic water storage, and surface and subsurface run‐off generation in a drained and farmed closed depression in north‐eastern Indiana, USA. Results showed that both surface flow and subsurface tile flow displayed a threshold relationship with the sum of rainfall amount and soil moisture deficit (SMD). Neither surface flow nor subsurface tile flow was observed unless rainfall amount exceeded the SMD. Timing of subsurface tile flow relative to soil moisture response on the shoulder slope of the depression indicated that the formation and drainage of perched water tables on depression hillslopes were likely the main mechanism that produced subsurface connectivity. Surface flow generation was delayed compared with subsurface tile flow during rainfall events due to differences in soil water storage along depression hillslopes and run‐off generation mechanisms. These findings highlight the substantial impact of subsurface tile drainage on the hydrology of closed depressions; the bottom of the depression, the wettest area prior to drainage installation, becomes the driest part of the depression after installation of subsurface drainage. Rapid connectivity of localized subsurface saturation zones during rainfall events is also greatly enhanced because of subsurface drainage. Thus, less fill is required to generate substantial spill. Understanding hydrologic processes in drained and farmed closed depressions is a critical first step in developing improved water and nutrient management strategies in this landscape.     

Homogenization of spatial patterns of hydrologic response in artificially drained agricultural catchments

Anthropogenic modifications to the landscape, with agricultural activities being a primary driver, have resulted in significant alterations to the hydrologic cycle. Artificial drainage, including surface and subsurface drainage (tile drains), is one of the most extensive manipulations in agricultural landscapes and thus is expected to provide a distinct signature of anthropogenic modification. This study adopts a data synthesis approach in an effort to characterize the signature of artificial subsurface drainage. Daily discharge data from 24 basins across the state of Iowa, which encapsulate a range of anthropogenic modifications, are assessed using a variety of flow metrics. Results indicate that the presence of artificial subsurface drainage leads to a homogenization of landscape hydrologic response. Non‐tiled watersheds exhibit a decrease in the area‐normalized peak discharge and an increase in the baseflow ratio (baseflow/streamflow) with increases in the spatial scale, while scale invariance is apparent in tiled basins. Within‐basin variability in hydrograph recession coefficients also appears to decrease with increases in the proportion of the catchment that is artificially drained. Finally, the differences between tiled and non‐tiled landscapes disappear at scales greater than approximately 2200 km², indicating that this may be a threshold scale for studying the effects of tile drainage. This decrease in within‐basin variability and the scale invariance of hydrologic metrics in artificially drained watersheds are attributed to the creation of a bypass flow hydrologic pathway that bypasses the complexity of the catchment travel paths. Spatial homogeneity in responses implies that it may be possible to develop more parsimonious hydrologic models for these regions. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

Aspect control of water movement on hillslopes near the rain–snow transition of the Colorado Front Range

In the Colorado Front Range, forested catchments near the rain–snow transition are likely to experience changes in snowmelt delivery and subsurface water transport with climate warming and associated shifts in precipitation patterns. Snowpack dynamics are strongly affected by aspect: Lodgepole pine forested north‐facing slopes develop a seasonal snowpack, whereas Ponderosa pine‐dotted south‐facing slopes experience intermittent snow accumulation throughout winter and spring. We tested the degree to which these contrasting water input patterns cause different near‐surface hydrologic response on north‐facing and south‐facing hillslopes during the snowmelt period. During spring snowmelt, we applied lithium bromide (LiBr) tracer to instrumented plots along a north–south catchment transect. Bromide broke through immediately at 10‐ and 30‐cm depths on the north‐facing slope and was transported out of soil waters within 40 days. On the south‐facing slope, Br^? was transported to significant depths only during spring storms and remained above the detection limit throughout the study. Modelling of unsaturated zone hydrologic response using Hydrus‐1D corroborated these aspect‐driven differences in subsurface transport. Our multiple lines of evidence suggest that north‐facing slopes are dominated by connected flow through the soil matrix, whereas south‐facing slope soils experience brief periods of rapid vertical transport following snowmelt events and are drier overall than north‐facing slopes. These differences in hydrologic response were largely a function of energy‐driven differences in water supply, emphasizing the importance of aspect and climate forcing when considering contributions of water and solutes to streamflow in catchments near the snow line. Copyright © 2012 John Wiley & Sons, Ltd.     

Impacts of model initialization on an integrated surface water–groundwater model

Integrated hydrologic models characterize catchment responses by coupling the subsurface flow with land surface processes. One of the major areas of uncertainty in such models is the specification of the initial condition and its influence on subsequent simulations. A key challenge in model initialization is that it requires spatially distributed information on model states, groundwater levels and soil moisture, even when such data are not routinely available. Here, the impact of uncertainty in initial condition was explored across a 208 km² catchment in Denmark using the ParFlow.CLM model. The initialization impact was assessed under two meteorological conditions (wet vs dry) using five depth to water table and soil moisture distributions obtained from various equilibrium states (thermal, root zone, discharge, saturated and unsaturated zone equilibrium) during the model spin‐up. Each of these equilibrium states correspond to varying computation times to achieve stability in a particular aspect of the system state. Results identified particular sensitivity in modelled recharge and stream flow to the different initializations, but reduced sensitivity in modelled energy fluxes. Analysis also suggests that to simulate a year that is wetter than the spin‐up period, an initialization based on discharge equilibrium is adequate to capture the direction and magnitude of surface water–groundwater exchanges. For a drier or hydrologically similar year to the spin‐up period, an initialization based on groundwater equilibrium is required. Variability of monthly subsurface storage changes and discharge bias at the scale of a hydrological event show that the initialization impacts do not diminish as the simulations progress, highlighting the importance of robust and accurate initialization in capturing surface water–groundwater dynamics. Copyright © 2015 John Wiley & Sons, Ltd.     

Catchment‐scale contribution of forest roads to stream exports of sediment,phosphorus and nitrogen

The relative contribution of forest roads to total catchment exports of suspended sediment, phosphorus, and nitrogen was estimated for a 13 451 ha forested catchment in southeastern Australia. Instrumentation was installed for 1 year to quantify total in‐stream exports of suspended sediment, phosphorus, and nitrogen. In addition, a total of 101 road–stream crossings were mapped and characterized in detail within the catchment to identify the properties of the road section where the road network and the stream network intersect. Sediment and nutrient generation rates from different forest road types within the catchment were quantified using permanent instrumentation and rainfall simulation. Sediment and nutrient generation rates, mapped stream crossing information, traffic data and annual rainfall data were used to estimate annual loads of sediment, phosphorus, and nitrogen from each stream crossing in the catchment. The annual sum of these loads was compared with the measured total catchment exports to estimate the proportional contribution of loads from roads within the catchment. The results indicated that 3·15 ha of near‐stream unsealed road surface with an average slope of 8·4% delivered an estimated 50 t of the 1142 t of total suspended sediment exported from the catchment, or about 4·4% of the total sediment load from the forest. Stream discharge over this period was 69 573 Ml. The unsealed road network delivered an estimated maximum of 22 kg of the 1244 kg of total phosphorus from the catchment, or less than 1·8% of the total load from the forest. The average sediment and phosphorous load per crossing was estimated at 0·5 t (standard deviation 1·0 t) and 0·22 kg (standard deviation 0·30 kg) respectively. The lower proportional contribution of total phosphorus resulted from a low ratio of total phosphorus to total suspended sediment for the road‐derived sediment. The unsealed road network delivered approximately 33 kg of the 20 163 kg of total nitrogen, about 0·16% of the total load of nitrogen from the forest. The data indicate that, in this catchment, improvement of stream crossings would yield only small benefits in terms of net catchment exports of total suspended sediment and total phosphorus, and no benefit in terms of total nitrogen. These results are for a catchment with minimal road‐related mass movement, and extrapolation of these findings to the broader forested estate requires further research. Copyright © 2007 John Wiley & Sons, Ltd.     

Simulated response of an intermittent stream to rainfall frequency patterns

The response of intermittent catchments to rainfall is complex and difficult to model. This study uses the spatially distributed CATchment HYdrology (CATHY) model to explore how the frequency of daily rainfall (λ) can affect the hydrologic regime of intermittent catchments. After a multi-objective calibration and validation of CATHY against experimental measurements of streamflow and groundwater levels in a catchment used as a pasture, the role of λ in affecting streamflow characteristics was explored using different scenarios. With different values of λ for the dry and wet periods of the year, CATHY showed that a series of frequent rainfall events was often associated with incipient streamflow, independent of the season. Activation of streamflow during the wet season was related to multiple factors and was not often associated with the shallow groundwater levels near the outlet of the catchment. The interplay between rainfall depth and intensity acted as the most important factor for the generation of streamflow. Using the difference between accumulated rainfall and evapotranspiration as a measure of wetness, saturated subsurface flow mechanism generated streamflow in simulations with wetness at least three times larger than mean wetness of other simulations. Although groundwater uprise near the outlet did not effectively contribute to streamflow in the initial days of flow, it strongly correlated with the magnitude of the runoff coefficient. Values of λ close or equal to the maximum value in the wet season can sustain the connectivity between groundwater and streamflow in the riparian zone. This connectivity increases the catchment wetness, which consequently results in an increase of the generated streamflow. Our study showed that rainfall regimes characterized by different λ were able to identify distinct flow regimes typical of either intermittent, ephemeral, or nonflowing catchments. Decrease of λ in the wet season is likely associated with a reduction of streamflow, with a shift of flow regime from intermittent to ephemeral or no-flow.