Comparison of stormflow responses of surface‐mined and forested watersheds in the Appalachian Mountains,USA

The results of a hydrological analysis that was conducted as part of a larger, multifaceted, collaborative effort to quantify ecosystem functions in watersheds subjected to land‐use and land‐cover change are presented. The primary goal of the study was to determine whether a small watershed in the Appalachian region (USA) that was recently subjected to surface mining and reclamation practices produces stormflow responses to rain events that are different from those produced by a nearby reference watershed covered by young, second‐growth forest. Water balances indicated that runoff yields did not vary significantly between the two watersheds on an annual basis. Statistically significant differences (p?0·05) in runoff responses were observed on an event basis, however, with the mined/reclaimed watershed producing, on average (a) higher storm runoff coefficients (2·5×), (b) greater total storm runoff (3×), and (c) higher peak hourly runoff rates (2×) when compared with the reference watershed. Results of a unit hydrograph analysis also showed, unexpectedly, that the modelled unit responses of the two watersheds to effective rainfall pulses were similar, despite the noted differences in land cover. Differences in stormflow responses were thus largely explained by dramatic reductions in cumulative rates of rainfall abstraction (measured using infiltrometers) attributable to soil compaction during land reclamation. Additional field hydrological measurements on other mined watersheds will be needed to generalize our results, as well as to understand and predict the cumulative hydrological impacts of widespread surface mining in larger watersheds and river basins. Copyright © 2006 John Wiley & Sons, Ltd.     

A simplified model for estimating field‐scale surface runoff hydrographs

The problem of obtaining field‐scale surface response to rainfall events is complicated by the spatial variability of infiltration characteristics of the soil and rainfall. In this paper, we develop and test a simplified model for generating surface runoff over fields with spatial variation in both rainfall rate and saturated hydraulic conductivities. The model is able to represent the effects of local variation in infiltration, as well as the run‐on effect that controls infiltration of excess water from saturated upstream areas. The effective rainfall excess is routed to the slope outlet using a simplified solution of the kinematic wave approximation. Model results are compared to averaged hydrographs from numerically‐intensive Monte–Carlo simulations for observed and design rainfall events and soil patterns that are typical of Central Italy. The simplified model is found to yield satisfactory results at a relatively small computational expense. A proposal to include a simple channel routing scheme is also presented as a prelude to extend this conceptualization to watershed scales. Copyright © 2006 John Wiley & Sons, Ltd.     

Vegetation controls on small‐scale runoff and erosion dynamics in a degrading dryland environment

This paper investigates the controls of vegetation on runoff and erosion dynamics in the dryland environment of Jornada, New Mexico, USA. As the American southwest has seen significant shifts in the dominant vegetation species in the past 150 years, an understanding of the vegetation effects on hydrological and erosional processes is vital for understanding and managing environmental change. Small‐scale rainfall simulations were carried out to identify the hydrological and erosional processes resulting from the grassland and shrubland vegetation species. Results obtained using tree‐regression analysis suggested that the primary vegetation control on runoff and erosion is the shrub type and canopy density, which directly affects the local microtopographic gradient of mounds beneath the shrubs. Significant interactions and feedbacks were found to occur among the local mound gradient, crust cover, soil aggregate stability and antecedent soil moisture between the different vegetation species for both the runoff and erosion responses. Although some of the shrub species were found to produce higher sediment yields than the grass species, the distinguishing feature of the grassland was the significantly higher enrichment in the fine sediment fraction compared to all other surface cover types. This enrichment in fines has important implications for nutrient movement in such environments. Copyright © 2009 John Wiley & Sons, Ltd.     

Effects of forest roads on the hydrological response of a small‐scale mountain watershed in Greece

The effects of land use changes on the ecology and hydrology of natural watersheds have long been debated. However, less attention has been given to the hydrological effects of forest roads. Although less studied, several researchers have claimed that streamflow changes related to forest roads can cause a persistent and pervasive effect on hillslope hydrology and the functioning of the channel system. The main potential direct effects of forest roads on natural watersheds hydrologic response are runoff production on roads surfaces due to reduced infiltration rates, interruption of subsurface flow by road cutslopes and rapid transfer of the produced runoff to the stream network through roadside ditches. The aforementioned effects may significantly modify the total volume and timing of the hillslope flow to the stream network. This study uses detailed field data, spatial data, hydro‐meteorological records, as well as numerical simulation to investigate the effects of forest roads on the hydrological response of a small‐scale mountain experimental watershed, which is situated in the east side of Penteli Mountain, Attica, Greece. The results of this study highlight the possible effects of forest roads on the watersheds hydrological response that may significantly influence direct runoff depths and peak flow rates. It is demonstrated that these effects can be very important in permeable watersheds and that more emphasis should be given on the impact of roads on the watersheds hydrological response. Copyright © 2014 John Wiley & Sons, Ltd.     

Estimation of surface runoff and water‐covered area during filling of surface microrelief depressions

During the filling of surface microrelief depressions the precipitation excess (precipitation minus infiltration and interception) is divided between surface storage and runoff, i.e. runoff starts before the surface depressions are filled. Information on the division of precipitation excess is needed for modelling surface runoff during the filling of surface depressions. Furthermore, information on the surface of the area covered with water is needed for calculating infiltration of water stored in soil surface depressions. Thirty‐two soil surface microreliefs were determined in Danish erosion study plots. The slope was c. 10% for all plots. Data were treated initially by removing the slope, after which 20 ‘artificial’ slopes (1–20%) were introduced producing 640 new data sets. Runoff during filling of the microrelief storage was calculated for each of the 640 data sets using a model developed for calculating surface storage and runoff from grid elevation measurements. Runoff started immediately after the first addition of water for all data sets. On a field scale, however, runoff has to travel some distance as overland flow and storage in smaller and larger depressions below the runoff initiation point must be taken into consideration. The runoff increases by intermittent steps. Whenever a depression starts to overflow to the border of the plot, the runoff jumps accordingly. In spite of the jumps, the distribution between surface storage and runoff was closely related to the quotient between precipitation excess and depression storage capacity. Surface area covered with water was exponentially related to the amount of water stored in surface depressions. Models for calculating surface storage and runoff from grid elevation measurements are cumbersome and require time‐consuming measurements of the soil surface microrelief. Therefore, estimation from roughness indices requiring fewer measurements is desirable. New improved equations for such estimations are suggested. Copyright © 2000 John Wiley & Sons, Ltd.     

Measurement of surface water runoff from plots of two different sizes

Intensities and amounts of water infiltration and runoff on sloping land are governed by the rainfall pattern and soil hydraulic conductivity, as well as by the microtopography and soil surface conditions. These components are closely interrelated and occur simultaneously, and their particular contribution may change during a rainfall event, or their effects may vary at different field scales. The scale effect on the process of infiltration/runoff was studied under natural field and rainfall conditions for two plot sizes: small plots of 0·25 m² and large plots of 50 m². The measurements were carried out in the central region of Chile in a piedmont most recently used as natural pastureland. Three blocks, each having one large plot and five small plots, were established. Cumulative rainfall and runoff quantities were sampled every 5 min. Significant variations in runoff responses to rainfall rates were found for the two plot sizes. On average, large plots yielded only 40% of runoff quantities produced on small plots per unit area. This difference between plot sizes was observed even during periods of continuous runoff. Copyright © 2002 John Wiley & Sons, Ltd.     

Infiltration and runoff generation processes in fire‐affected soils

Post‐wildfire runoff was investigated by combining field measurements and modelling of infiltration into fire‐affected soils to predict time‐to‐start of runoff and peak runoff rate at the plot scale (1 m²). Time series of soil‐water content, rainfall and runoff were measured on a hillslope burned by the 2010 Fourmile Canyon Fire west of Boulder, Colorado during cyclonic and convective rainstorms in the spring and summer of 2011. Some of the field measurements and measured soil physical properties were used to calibrate a one‐dimensional post‐wildfire numerical model, which was then used as a ‘virtual instrument’ to provide estimates of the saturated hydraulic conductivity and high‐resolution (1 mm) estimates of the soil‐water profile and water fluxes within the unsaturated zone. Field and model estimates of the wetting‐front depth indicated that post‐wildfire infiltration was on average confined to shallow depths less than 30 mm. Model estimates of the effective saturated hydraulic conductivity, K_s, near the soil surface ranged from 0.1 to 5.2 mm h^?1. Because of the relatively small values of K_s, the time‐to‐start of runoff (measured from the start of rainfall), t_p, was found to depend only on the initial soil‐water saturation deficit (predicted by the model) and a measured characteristic of the rainfall profile (referred to as the average rainfall acceleration, equal to the initial rate of change in rainfall intensity). An analytical model was developed from the combined results and explained 92–97% of the variance of t_p, and the numerical infiltration model explained 74–91% of the variance of the peak runoff rates. These results are from one burned site, but they strongly suggest that t_p in fire‐affected soils (which often have low values of K_s) is probably controlled more by the storm profile and the initial soil‐water saturation deficit than by soil hydraulic properties. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

Comparison of two stream gauging systems for measuring runoff and sediment yield for a semi‐arid watershed

Our ability to understand erosion processes in semi‐arid ecosystems depends on establishing relationships between rainfall and runoff. This requires collection of extensive and accurate hydrologic and sediment data sets. A supercritical flume with a total load traversing slot sediment sampler used on several sites at the Walnut Gulch Experimental Watershed (WGEW) near Tombstone, AZ has proven to be a reliable way to measure flow and sediment discharge from small watersheds. However, it requires installation of a costly structure that is only suitable for relatively small flows. A more commonly used method based on ease of installation and expense is the pump sampler. One example of this is a set of instrumentation developed by the Australian Commonwealth Scientific and Industrial Research Organization (CSIRO), in which the pump sediment sampler is part of an in‐channel, fully automated system for measuring water velocity, depth, turbidity and collecting runoff samples. A 3.7 ha arid watershed at WGEW was instrumented with both systems and hydrologic and sediment data were collected and compared during a 2 year period. Total sediment yield for the entire period measured by the CSIRO pump sampler (11.6 t ha^‐1) was similar to that by traversing slot sampler (11.5 t ha^‐1). The pump sampler accurately estimated the amount of fine (< 0.5 mm) sediment fractions exported, but consistently underestimated the coarse (>0.5 mm) sediment fractions. Median sediment diameter of samples collected by traversing slot and pump sampler were 0.32 and 0.22 mm, respectively. This study outlines the benefits and limitations of the pump sampler based system for monitoring sediment concentration and yield in high‐energy headwater catchments, and makes recommendations for improvement of its performance. Copyright © 2012 John Wiley & Sons, Ltd.     

A rainfall‐runoff model parameterized from GIS and runoff data

A guiding principle in hydrological modelling should be to keep the number of calibration parameters to a minimum. A reduced number of parameters to be calibrated, while maintaining the accuracy and detail required by modern hydrological models, will reduce parameter and model structure uncertainty and improve model diagnostics. In this study, the dynamics of runoff are derived from the distribution of distances from points in the catchments to the nearest stream. This distribution is unique for each catchment and can be determined from a geographical information system. The distribution of distances, will, when a celerity of (subsurface) flow is introduced, provide a distribution of travel times, or a unit hydrograph (UH). For spatially varying levels of saturation deficit, we have different celerities and, hence, different UHs. Runoff is derived from the superposition of the different UHs. This study shows how celerities can be estimated if we assume that recession events represent the combined UHs for different levels of saturation deficit. A new soil moisture routine which estimates saturated and unsaturated volumes of subsurface water and with only one parameter to calibrate is included in the new model. The performance of the new model is compared with that of the Swedish HBV model and is found to perform equally well for eight Norwegian catchments although the number of parameters to be calibrated in the module concerning soil moisture and runoff dynamics is reduced from seven in the HBV model to one in the new model. It is also shown that the new model has a more realistic representation of the subsurface hydrology. Copyright © 2013 John Wiley & Sons, Ltd.     

Quantifying and simulating stormwater runoff in watersheds

We developed the Stormwater Runoff Modeling System (SWARM) based on curve number and unit hydrograph methods of the U.S. Department of Agriculture, Natural Resources Conservation Service. SWARM models single events, targets watersheds fitting easily within hydrologic units with 12‐digit codes, and has been calibrated for low‐gradient topography of the Southeast coastal plain. We established protocols; made changes related to peak rate factors, travel time formulas, curve numbers, and the initial abstraction ratio; and then tested the output with multi‐site validation using U.S. Geological Survey measurements of discharge and rainfall. Validation results from both undeveloped and developed watersheds support the robustness of our system in quantifying and simulating runoff: rainfall to runoff differences between measured and simulated volumes ranged from 3 to 11%; r² for hydrograph curves ranged from 0.82 to 0.98. SWARM can be a useful tool for scientific research and for coastal resource management and decision making in the Southeast coastal plain specifically and also may be applied to other areas by recalibrating parameters and modifying calculation templates. Copyright © 2012 John Wiley & Sons, Ltd.     

Delineating runoff processes and critical runoff source areas in a pasture hillslope of the Ozark Highlands

The identification of runoff contributing areas would provide the ideal focal points for water quality monitoring and Best Management Practice (BMP) implementation. The objective of this study was to use a field‐scale approach to delineate critical runoff source areas and to determine the runoff mechanisms in a pasture hillslope of the Ozark Highlands in the USA. Three adjacent hillslope plots located at the Savoy Experimental Watershed, north‐west Arkansas, were bermed to isolate runoff. Each plot was equipped with paired subsurface saturation and surface runoff sensors, shallow groundwater wells, H‐flumes and rain gauges to quantify runoff mechanisms and rainfall characteristics at continuous 5‐minute intervals. The spatial extent of runoff source areas was determined by incorporating sensor data into a geographic information‐based system and performing geostatistical computations (inverse distance weighting method). Results indicate that both infiltration excess runoff and saturation excess runoff mechanisms occur to varying extents (0–58% for infiltration excess and 0–26% for saturation excess) across the plots. Rainfall events that occurred 1–5 January 2005 are used to illustrate the spatial and temporal dynamics of the critical runoff source areas. The methodology presented can serve as a framework upon which critical runoff source areas can be identified and managed for water quality protection in other watersheds. Copyright © 2008 John Wiley & Sons, Ltd.     

Connectivity of post-fire runoff and sediment from nested hillslopes and watersheds

Wildfire increases the potential connectivity of runoff and sediment throughout watersheds due to greater bare soil, runoff and erosion as compared to pre-fire conditions. This research examines the connectivity of post-fire runoff and sediment from hillslopes (< 1.5 ha; n = 31) and catchments (< 1000 ha; n = 10) within two watersheds (< 1500 ha) burned by the 2012 High Park Fire in northcentral Colorado, USA. Our objectives were to: (1) identify sources and quantify magnitudes of post-fire runoff and erosion at nested hillslopes and watersheds for two rain storms with varied duration, intensity and antecedent precipitation; and (2) assess the factors affecting the magnitude and connectivity of runoff and sediment across spatial scales for these two rain storms. The two summer storms that are the focus of this research occurred during the third summer after burning. The first storm had low intensity rainfall over 11 hours (return interval <1–2 years), whereas the second event had high intensity rainfall over 1 hour (return interval <1–10 years). The lower intensity storm was preceded by high antecedent rainfall and led to low hillslope sediment yields and channel incision at most locations, whereas the high intensity storm led to infiltration-excess overland flow, high sediment yields, in-stream sediment deposition and channel substrate fining. For both storms, hillslope-to-stream sediment delivery ratios and area-normalised cross-sectional channel change increased with the percent of catchment that burned at high severity. For the high intensity storm, hillslope-to-stream sediment delivery ratios decreased with unconfined channel length (%). The findings quantify post-fire connectivity and sediment delivery from hillslopes and streams, and highlight how different types of storms can cause varying magnitues and spatial patterns of sediment transport and deposition from hillslopes through stream channel networks.     

Landscape influence on small‐scale water temperature variations in a moorland catchment

We monitored temperatures in stream water, groundwater and riparian wetland surface water over 18 months in a 3.2‐km² moorland catchment in the Scottish Highlands. The stream occupies a glaciated valley, aligned east–west. It has three main headwater tributaries with a large north facing catchment, a south facing catchment and the smallest east facing headwater. The lower catchment sampling locations begin after the convalescence of all three headwaters. Much of the stream network is fringed by riparian peatlands. Stream temperatures are mainly regulated by energy exchanges at the air–water interface. However, they are also influenced by inflows from the saturated riparian zone, where surface water source areas are strongly connected with the stream network. Consequently, the spatial distribution of stream temperatures exhibits limited variability. Nevertheless, there are significant summer differences between the headwaters, despite their close proximity to each other. This is consistent with aspect (and incident radiation), given the south and east facing headwaters having higher temperatures. The largest, north‐facing sub‐catchment shows lower summer diurnal temperature variability, suggesting that lower radiation inputs dampen temperature extremes. Whilst stream water temperature regimes in the lower catchment exhibit little change along a 1‐km reach, they are similar to those in the largest headwater; probably reflecting size and comparable catchment aspect and hydrological flow paths. Our results suggest that different parts of the channel network and its connected wetlands have contrasting sensitivity to higher summer temperatures. This may be important in land management strategies designed to mitigate the impacts of projected climatic warming. Copyright © 2015 John Wiley & Sons, Ltd.     

Determination of runoff and sediment yield from a small watershed in sub‐humid subtropics using the HSPF model

The Hydrologic Simulation Programme‐Fortran (HSPF), a hydrologic and water quality computer model, was employed for simulating runoff and sediment yield during the monsoon months (June–October) from a small watershed situated in a sub‐humid subtropical region of India. The model was calibrated using measured runoff and sediment yield data for the monsoon months of 1996 and was validated for the monsoon months of 2000 and 2001. During the calibration period, daily‐calibrated runoff had a Nash‐Sutcliffe efficiency (E_NS) value of 0·68 and during the validation period it ranged from 0·44 to 0·67. For daily sediment yield E_NS was 0·71 for the calibration period and it ranged from 0·68 to 0·90 for the validation period. Sensitivity analysis was performed to assess the impact of important watershed characteristics. The model parameters obtained in this study could serve as reference values for model application in similar climatic regions, with practical implications in watershed planning and management and designing best management practices. Copyright © 2007 John Wiley & Sons, Ltd.     

Lumped surface and sub‐surface runoff for erosion modeling within a small hilly watershed in northern Vietnam

Developing models to predict on‐site soil erosion and off‐site sediment transport at the agricultural watershed scale represent an on‐going challenge in research today. This study attempts to simulate the daily discharge and sediment loss using a distributed model that combines surface and sub‐surface runoffs in a small hilly watershed (< 1 km²). The semi‐quantitative model, Predict and Localize Erosion and Runoff (PLER), integrates the Manning–Strickler equation to simulate runoff and the Griffith University Erosion System Template equation to simulate soil detachment, sediment storage and soil loss based on a map resolution of 30 m × 30 m and over a daily time interval. By using a basic input data set and only two calibration coefficients based, respectively, on water velocity and soil detachment, the PLER model is easily applicable to different agricultural scenarios. The results indicate appropriate model performance and a high correlation between measured and predicted data with both Nash–Sutcliffe efficiency (Ef) and correlation coefficient (r²) having values > 0.9. With the simple input data needs, PLER model is a useful tool for daily runoff and soil erosion modeling in small hilly watersheds in humid tropical areas. Copyright © 2013 John Wiley & Sons, Ltd.     

Geospatially based distributed rainfall‐runoff modelling for simulation of internal and outlet responses in a semi‐forested lower Himalayan watershed

This study presents a Geographic Information System (GIS)‐based distributed rainfall‐runoff model for simulating surface flows in small to large watersheds during isolated storm events. The model takes into account the amount of interception storage to be filled using a modified Merriam ( 1960 ) approach before estimating infiltration by the Smith and Parlange ( 1978 ) method. The mechanics of overland and channel flow are modelled by the kinematic wave approximation of the Saint Venant equations which are then numerically solved by the weighted four‐point implicit finite difference method. In this modelling the watershed was discretized into overland planes and channels using the algorithms proposed by Garbrecht and Martz ( 1999 ). The model code was first validated by comparing the model output with an analytical solution for a hypothetical plane. Then the model was tested in a medium‐sized semi‐forested watershed of Pathri Rao located in the Shivalik ranges of the Garhwal Himalayas, India. Initially, a local sensitivity analysis was performed to identify the parameters to which the model outputs like runoff volume, peak flow and time to peak flow are sensitive. Before going for model validation, calibration was performed using the Ordered‐Physics‐based Parameter Adjustment (OPPA) method. The proposed Physically Based Distributed (PBD) model was then evaluated both at the watershed outlet as well as at the internal gauging station, making this study a first of its kind in Indian watersheds. The results of performance evaluation indicate that the model has simulated the runoff hydrographs reasonably well within the watershed as well as at the watershed outlet with the same set of calibrated parameters. The model also simulates, realistically, the temporal variation of the spatial distribution of runoff over the watershed and the same has been illustrated graphically. Copyright © 2011 John Wiley & Sons, Ltd.     

Comparison of simple runoff models used in Korea for small watersheds

This study analysed runoff characteristics for three small watersheds in the Su‐Young River basin in Korea using the storage function method, the linear reservoir cascade model, and the discrete linear input–output model. The models performed well for the watersheds. The runoff analysis of the observed flood events showed that the storage function method was the most accurate of the three, and was followed by the linear reservoir cascade model and the discrete linear input–output model. Copyright © 1999 John Wiley & Sons, Ltd.     

Prescribed‐fire effects on rill and interrill runoff and erosion in a mountainous sagebrush landscape

Changing fire regimes and prescribed‐fire use in invasive species management on rangelands require improved understanding of fire effects on runoff and erosion from steeply sloping sagebrush‐steppe. Small (0·5 m²) and large (32·5 m²) plot rainfall simulations (85 mm h^–1, 1 h) and concentrated flow methodologies were employed immediately following burning and 1 and 2 years post‐fire to investigate infiltration, runoff and erosion from interrill (rainsplash, sheetwash) and rill (concentrated flow) processes on unburned and burned areas of a steeply sloped sagebrush site on coarse‐textured soils. Soil water repellency and vegetation were assessed to infer relationships in soil and vegetation factors that influence runoff and erosion. Runoff and erosion from rainfall simulations and concentrated flow experiments increased immediately following burning. Runoff returned to near pre‐burn levels and sediment yield was greatly reduced with ground cover recovery to 40 per cent 1 year post‐fire. Erosion remained above pre‐burn levels on large rainfall simulation and concentrated flow plots until ground cover reached 60 per cent two growing seasons post‐fire. The greatest impact of the fire was the threefold reduction of ground cover. Removal of vegetation and ground cover and the influence of pre‐existing strong soil‐water repellency increased the spatial continuity of overland flow, reduced runoff and sediment filtering effects of vegetation and ground cover, and facilitated increased velocity and transport capacity of overland flow. Small plot rainfall simulations suggest ground cover recovery to 40 per cent probably protected the site from low‐return‐interval storms, large plot rainfall and concentrated flow experiments indicate the site remained susceptible to elevated erosion rates during high‐intensity or long duration events until ground cover levels reached 60 per cent. The data demonstrate that the persistence of fire effects on steeply‐sloped, sandy sagebrush sites depends on the time period required for ground cover to recover to near 60 per cent and on the strength and persistence of ‘background’ or fire‐induced soil water repellency. Published in 2009 by John Wiley & Sons, Ltd.     

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

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

Simulating pan‐Arctic runoff with a macro‐scale terrestrial water balance model

A terrestrial hydrological model, developed to simulate the high‐latitude water cycle, is described, along with comparisons with observed data across the pan‐Arctic drainage basin. Gridded fields of plant rooting depth, soil characteristics (texture, organic content), vegetation, and daily time series of precipitation and air temperature provide the primary inputs used to derive simulated runoff at a grid resolution of 25 km across the pan‐Arctic. The pan‐Arctic water balance model (P/WBM) includes a simple scheme for simulating daily changes in soil frozen and liquid water amounts, with the thaw–freeze model (TFM) driven by air temperature, modelled soil moisture content, and physiographic data. Climate time series (precipitation and air temperature) are from the National Centers for Environmental Prediction (NCEP) reanalysis project for the period 1980–2001. P/WBM‐generated maximum summer active‐layer thickness estimates differ from a set of observed data by an average of 12 cm at 27 sites in Alaska, with many of the differences within the variability (1σ) seen in field samples. Simulated long‐term annual runoffs are in the range 100 to 400 mm year^?1. The highest runoffs are found across northeastern Canada, southern Alaska, and Norway, and lower estimates are noted along the highest latitudes of the terrestrial Arctic in North America and Asia. Good agreement exists between simulated and observed long‐term seasonal (winter, spring, summer–fall) runoff to the ten Arctic sea basins (r = 0·84). Model water budgets are most sensitive to changes in precipitation and air temperature, whereas less affect is noted when other model parameters are altered. Increasing daily precipitation by 25% amplifies annual runoff by 50 to 80% for the largest Arctic drainage basins. Ignoring soil ice by eliminating the TFM sub‐model leads to runoffs that are 7 to 27% lower than the control run. The results of these model sensitivity experiments, along with other uncertainties in both observed validation data and model inputs, emphasize the need to develop improved spatial data sets of key geophysical quantities (particularly climate time series) to estimate terrestrial Arctic hydrological budgets better. Copyright © 2003 John Wiley & Sons, Ltd.