Application of a watershed runoff model to North-east Pond River,Newfoundland: to study water balance and hydrological characteristics owing to atmospheric change

The hydrological sensitivities to long-term climate change of a watershed in Eastern Canada were analysed using a deterministic watershed runoff model developed to simulate watershed acidification. This model was modified to study atmospheric change effects in the watershed. Water balance modelling techniques, modified for assessing climate effects, were developed and tested for a watershed using atmospheric change scenarios from both state of the art general circulation models and a series of hypothetical scenarios. The model computed daily surface, inter- and groundwater flows from the watershed. The moisture, infiltration and recharge rate are also computed in the soil reservoirs. The thirty years of simulated data can be used to evaluate the effects of climatic change on soil moisture, recharge rate and surface and subsurface flow systems. The interaction between surface and subsurface water is discussed in relation to climate change. These hydrological results raise the possibility of major environmental and socioeconomic difficulties and have significant implications for future water resource planning and management. © 1997 John Wiley & Sons, Ltd.     

Vertical surface water–groundwater exchange processes within a headwater floodplain induced by experimental floods

Restoring hydrologic connectivity between channels and floodplains is common practice in stream and river restoration. Floodplain hydrology and hydrogeology impact stream hydraulics, ecology, biogeochemical processing, and pollutant removal, yet rigorous field evaluations of surface water–groundwater exchange within floodplains during overbank floods are rare. We conducted five sets of experimental floods to mimic floodplain reconnection by pumping stream water onto an existing floodplain swale. Floods were conducted throughout the year to capture seasonal variation and each involved two replicate floods on successive days to test the effect of varying antecedent moisture. Water levels and specific conductance were measured in surface water, soil, and groundwater within the floodplain, along with surface flow into and out of the floodplain. Vegetation density varied seasonally and controlled the volume of surface water storage on the floodplain. By contrast, antecedent moisture conditions controlled storage of water in floodplain soils, with drier antecedent moisture conditions leading to increased subsurface storage and slower flood wave propagation across the floodplain surface. The site experienced spatial heterogeneity in vertical connectivity between surface water and groundwater across the floodplain surface, where propagation of hydrostatic pressure, preferential flow, and bulk Darcy flow were all mechanisms that may have occurred during the five floods. Vertical connectivity also increased with time, suggesting higher frequency of floodplain inundation may increase surface water–groundwater exchange across the floodplain surface. Understanding the variability of floodplain impacts on water quality noted in the literature likely requires better accounting for seasonal variations in floodplain vegetation and antecedent moisture as well as heterogeneous exchange flow mechanisms. Copyright © 2016 John Wiley & Sons, Ltd.     

The role of overland flow and subsurface flow on the spatial distribution of soil moisture in the topsoil

Many investigations show relationships between topographical factors and the spatial distribution of soil moisture in catchments. However, few quantitative analyses have been carried out to elucidate the role of different hydrological processes in the spatial distribution of topsoil moisture in catchments. A spatially distributed rainfall—runoff model was used to investigate contributions of subsurface matric flow, macropore flow and surface runoff to the spatial distribution of soil moisture in a cultivated catchment. The model results show that lateral subsurface flow in the soil matrix or in macropores has a minor effect on the spatial distribution of soil moisture. Only when a perched groundwater table is maintained long enough, which is only possible if the subsurface is completely impermeable, may a spatial distribution in moisture content occur along the slope. Surface runoff, producing accumulations of soil moisture in flat flow paths of agricultural origin (field boundaries), was demonstrated to cause significant spatial variations in soil moisture within a short period after rainfall (<2 days). When significant amounts of surface runoff are produced, wetter moisture conditions will be generated at locations with larger upstream contributing areas. Copyright © 2001 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.     

A comparison of two physics-based numerical models for simulating surface water–groundwater interactions

Problems in hydrology and water management that involve both surface water and groundwater are best addressed with simulation models that can represent the interactions between these two flow regimes. In the current generation of coupled models, a variety of approaches is used to resolve surface–subsurface interactions and other key processes such as surface flow propagation. In this study we compare two physics-based numerical models that use a 3D Richards equation representation of subsurface flow. In one model, surface flow is represented by a fully 2D kinematic approximation to the Saint–Venant equations with a sheet flow conceptualization. In the second model, surface routing is performed via a quasi-2D diffusive formulation and surface runoff follows a rill flow conceptualization. The coupling between the land surface and the subsurface is handled via an explicit exchange term resolved by continuity principles in the first model (a fully-coupled approach) and by special treatment of atmospheric boundary conditions in the second (a sequential approach). Despite the significant differences in formulation between the two models, we found them to be in good agreement for the simulation experiments conducted. In these numerical tests, on a sloping plane and a tilted V-catchment, we examined saturation excess and infiltration excess runoff production under homogeneous and heterogeneous conditions, the dynamics of the return flow process, the differences in hydrologic response under rill flow and sheet flow parameterizations, and the effects of factors such as grid discretization, time step size, and slope angle. Low sensitivity to vertical discretization and time step size was found for the two models under saturation excess and homogeneous conditions. Larger sensitivity and differences in response were observed under infiltration excess and heterogeneous conditions, due to the different coupling approaches and spatial discretization schemes used in the two models. For these cases, the sensitivity to vertical and temporal resolution was greatest for processes such as reinfiltration and ponding, although the differences between the hydrographs of the two models decreased as mesh and step size were progressively refined. In return flow behavior, the models are in general agreement, with the largest discrepancies, during the recession phase, attributable to the different parameterizations of diffusion in the surface water propagation schemes. Our results also show that under equivalent parameterizations, the rill and sheet flow conceptualizations used in the two models produce very similar responses in terms of hydrograph shape and flow depth distribution.     

Modelling hydrological response to a fully‐monitored urban bioretention cell

Municipalities and agencies use green infrastructure to combat pollution and hydrological impacts (e.g., flooding) related to excess stormwater. Bioretention cells are one type of infiltration green infrastructure intervention that infiltrate and redistribute otherwise uncontrolled stormwater volume. However, the effects of these installations on the rest of the local water cycle is understudied; in particular, impacts on stormwater return flows and groundwater levels are not fully understood. In this study, full water cycle monitoring data were used to construct and calibrate a two‐dimensional Richards equation model (HYDRUS‐2D/3D) detailing hydrological implications of an unlined bioretention cell (Cleveland, Ohio) that accepts direct runoff from surrounding impervious surfaces. Using both preinstallation and postinstallation data, the model was used to (a) establish a mass balance to determine reduction in stormwater return flow, (b) evaluate green infrastructure effects on subsurface water dynamics, and (c) determine model sensitivity to measured soil properties. Comparisons of modelled versus observed data indicated that the model captured many hydrological aspects of the bioretention cell, including subsurface storage and transient groundwater mounding. Model outputs suggested that the bioretention cell reduced stormwater return flows into the local sewer collection system, though the extent of this benefit was attenuated during high inflow events that may have exhausted detention capacity. The model also demonstrated how, prior to bioretention cell installation, surface and subsurface hydrology were largely decoupled, whereas after installation, exfiltration from the bioretention cell activated a new groundwater dynamic. Still, the extent of groundwater mounding from the cell was limited in spatial extent and did not threaten other subsurface infrastructure. Finally, the sensitivity analysis demonstrated that the overall hydrological response was regulated by the hydraulics of the bioretention cell fill material, which controlled water entry into the system, and by the water retention parameters of the native soil, which controlled connectivity between the surface and groundwater.     

Comparison of two modeling approaches for groundwater–surface water interactions

An assessment of interactions between groundwater and surface water was carried out by applying two different modeling approaches to a small‐scale study area in the municipality of Havelock, Quebec. The first approach involved a commonly used sequential procedure that consists in determining the daily recharge rate using a quasi 2D infiltration model (HELP), applied in the next step as an imposed flux to a 3D finite‐element groundwater flow model. The flow model was calibrated under steady‐state and transient conditions against measured water levels. The second approach was based on a recently developed physically based, 3D fully coupled groundwater–surface water flow model (CATHY) applied to the entire flow domain in an integrated manner. Implementation, calibration, and results of the simulations for both approaches are presented and discussed. For equal annual precipitation (1038 mm/y) and evapotranspiration (556 mm/y), the second approach computed a recharge rate of 233 mm/y (8.9% higher than the first approach) and a net upward flow from the fractured aquifer (the first approach predicted a net downward flow to the rock). The simulated annual discharge was similar for the two approaches (9.6% difference). Both approaches were found to be useful in understanding the interactions between groundwater and surface water, although limitations are apparent in the sequential procedure's inability to account for surface–subsurface feedbacks, for instance near stream reaches where groundwater discharge is prevalent. The decoupled, two‐model approach provides disaggregated surface, vadose, and aquifer flows, and a simple aperçu at the different components of total discharge. The fully coupled model accounts for continuous water exchanges between the land surface, subsurface, and stream channel in a more complex manner, and produces a better match against observed data. Copyright © 2012 John Wiley & Sons, Ltd.     

The impact of subsurface conceptualization on land energy fluxes

There is a significant body of work demonstrating the importance of hydrologic control on land energy feedbacks. Yet, quantitative data on aquifer conductivity can be difficult to assemble. Furthermore, how subsurface uncertainty propagates into land-surface processes is not well understood. This study analyzes the impact of aquifer characterization on land energy fluxes, using a coupled hydrology–land-surface model. Four gridded subsurface conductivity fields are developed for the Upper Klamath basin using two data sources and different levels of imposed heterogeneity. Each model is forced with the same transient, observed meteorology for 3 years prior to the final year presented here. Results are analyzed to quantify the impact of subsurface heterogeneity on groundwater surface water interactions and spatial patterns in hydrologic variables. Analysis shows that heterogeneity does not fundamentally alter the connection between groundwater and land surface processes. However, differences between scenarios impact the extent and location of the critical zone.     

The subsurface–land surface–atmosphere connection under convective conditions

The dynamics of the free groundwater table influence land surface soil moisture and energy balance components, and are therefore also linked to atmospheric processes. In this study, the sensitivity of the atmosphere to groundwater table dynamics induced heterogeneity in land surface processes is examined under convective conditions. A fully coupled subsurface–land surface–atmosphere model is applied over a 150 km × 150 km study area located in Western Germany and ensemble simulations are performed over two convective precipitation events considering two separate model configurations based on groundwater table dynamics. Ensembles are generated by varying the model atmospheric initial conditions following the prescribed ensemble generation method by the German Weather Service in order to account for the intrinsic, internal atmospheric variability. The results demonstrate that especially under strong convective conditions, groundwater table dynamics affect atmospheric boundary layer height, convective available potential energy, and precipitation via the coupling with land surface soil moisture and energy fluxes. Thus, this study suggests that systematic uncertainties may be introduced to atmospheric simulations if groundwater table dynamics are neglected in the model.     

Role of bedrock groundwater in the rainfall–runoff process in a small headwater catchment underlain by volcanic rock

The role of bedrock groundwater in rainfall–runoff processes is poorly understood. Hydrometric, tracer and subsurface water potential observations were conducted to study the role of bedrock groundwater and subsurface flow in the rainfall–runoff process in a small headwater catchment in Shiranui, Kumamoto prefecture, south‐west Japan. The catchment bedrock consists of a strongly weathered, fractured andesite layer and a relatively fresh continuous layer. Major chemical constituents and stable isotopic ratios of δ¹⁸O and δD were analysed for spring water, rainwater, soil water and bedrock groundwater. Temporal and spatial variation in SiO₂ showed that stream flow under the base flow condition was maintained by bedrock groundwater. Time series of three components of the rainstorm hydrograph (rainwater, soil water and bedrock groundwater) separated by end member mixing analysis showed that each component fluctuated during rainstorm, and their patterns and magnitudes differed between events. During a typical mid‐magnitude storm event, a delayed secondary runoff peak with 1·0 l s⁻¹ was caused by increase in the bedrock groundwater component, whereas during a large rainstorm event the bedrock groundwater component increased to ≈ 2·5 l s⁻¹. This research shows that the contribution of bedrock groundwater and soil water depends strongly on the location of the groundwater table, i.e. whether or not it rises above the soil–bedrock interface. Copyright © 2010 John Wiley & Sons, Ltd.     

Snowmelt runoff and soil moisture dynamics on steep subalpine hillslopes

Snowmelt water supplies streamflow and growing season soil moisture in mountain regions, yet pathways of snowmelt water and their effects on moisture patterns are still largely unknown. This study examined how flow processes during snowmelt runoff affected spatial patterns of soil moisture on two steep sub‐alpine hillslope transects in Rocky Mountain National Park, CO, USA. The transects have northeast‐facing and east‐facing aspects, and both extend from high‐elevation bedrock outcrops down to streams in valley bottoms. Spatial patterns of both snow depth and near‐surface soil moisture were surveyed along these transects in the snowmelt and summer seasons of 2008–2010. To link these patterns to flow processes, soil moisture was measured continuously on both transects and compared with the timing of discharge in nearby streams. Results indicate that both slopes generated shallow lateral subsurface flow during snowmelt through near‐surface soil, colluvium and bedrock fractures. On the northeast‐facing transect, this shallow subsurface flow emerged through mid‐slope seepage zones, in some cases producing saturation overland flow, whereas the east‐facing slope had no seepage zones or overland flow. At the hillslope scale, earlier snowmelt timing on the east‐facing slope led to drier average soil moisture conditions than on the northeast‐facing slope, but within hillslopes, snow patterns had little relation to soil moisture patterns except in areas with persistent snow drifts. Results suggest that lateral flow and exfiltration processes are key controls on soil moisture spatial patterns in this steep sub‐alpine location. Copyright © 2014 John Wiley & Sons, Ltd.     

An analysis of mean transition time between flood and drought in the large river basin

 A soil moisture balance equation over large spatial regions is studied at seasonal and annual time scales for the Arkansas river basin. Interaction and feedback effects between land-surface and atmospheric moisture are studied in the parameterization for this basin. Due to the interaction between the land-surface and atmosphere at large scales, the surface hydrology of large land areas is susceptible to two distinct stable modes in the long-term probability density function: a dry and a wet state. In the soil moisture balance equation, stochastic fluctuations lead to separate preferred statistical stable states with transitions between these stable states induced by environmental fluctuations. On the basis of historical data, the soil moisture balance equation is calibrated for the Arkansas river basin. The transition times between the stable modes in the model are studied based on the stochastic representation of the physical processes and the calibrated model parameters. This study has implications for prediction of the transition times between stable modes or residence times, that is, the time the system spends in a given stable mode, since this would be equivalent to predicting the duration of droughts or wet conditions.     

Development and verification of a 3-D integrated surface water–groundwater model

Coupled modelling of surface and subsurface systems is a valuable tool for quantifying surface water–groundwater interactions. In the present paper, the 3-D non-steady state Navier–Stokes equations, after Reynolds averaging and with the assumption of a hydrostatic pressure distribution, are for the first time coupled to the 3-D saturated groundwater flow equations in an Integrated suRface watEr–grouNdwater modEl (IRENE). A finite-difference method is used for the solution of the governing equations of IRENE. A semi-implicit scheme is used for the discretisation of the surface water flow equations and a fully implicit scheme for the discretisation of the groundwater flow equations. The two sets of equations are coupled at the common interface of the surface water and groundwater bodies, where water exchange takes place, using Darcy’s law. A new approach is proposed for the solution of the coupled surface water and groundwater equations in a simultaneous manner, in such a fashion that gives computational efficiency at low computational cost. IRENE is verified against three analytical solutions of surface water–groundwater interaction, which are chosen so that different components of the model can be tested. The model closely reproduces the results of the analytical solutions and can therefore be used for analysing and predicting surface water–groundwater interactions in real-world cases.     

Evaluating the dual‐boundary forcing concept in subsurface–land surface interactions of the hydrological cycle

Subsurface and land surface processes (e.g. groundwater flow, evapotranspiration) of the hydrological cycle are connected via complex feedback mechanisms, which are difficult to analyze and quantify. In this study, the dual‐boundary forcing concept that reveals space–time coherence between groundwater dynamics and land surface processes is evaluated. The underlying hypothesis is that a simplified representation of groundwater dynamics may alter the variability of land surface processes, which may eventually affect the prognostic capability of a numerical model. A coupled subsurface–land surface model ParFlow.CLM is applied over the Rur catchment, Germany, and the mass and energy fluxes of the coupled water and energy cycles are simulated over three consecutive years considering three different lower boundary conditions (dynamic, constant, and free‐drainage) based on groundwater dynamics to substantiate the aforementioned hypothesis. Continuous wavelet transform technique is applied to analyze scale‐dependent variability of the simulated mass and energy fluxes. The results show clear differences in temporal variability of latent heat flux simulated by the model configurations with different lower boundary conditions at monthly to multi‐month time scales (~32–91 days) especially under soil moisture limited conditions. The results also suggest that temporal variability of latent heat flux is affected at even smaller time scales (~1–3 days) if a simple gravity drainage lower boundary condition is considered in the coupled model. This study demonstrates the importance of a physically consistent representation of groundwater dynamics in a numerical model, which may be important to consider in local weather prediction models and water resources assessments, e.g. drought prediction. Copyright © 2015 John Wiley & Sons, Ltd.     

IHMS—Integrated Hydrological Modelling System. Part 1. Hydrological processes and general structure

A newly Integrated Hydrological Modelling System (IHMS) has been developed to study the impact of changes in climate, land use and water management on groundwater and seawater intrusion (SWI) into coastal areas. The system represents the combination of three models, which can, if required, be run separately. It has been designed to assess the combined impact of climate, land use and groundwater abstraction changes on river, drainage and groundwater flows, groundwater levels and, where appropriate, SWI. The approach is interdisciplinary and reflects an integrated water management approach. The system comprises three packages: the Distributed Catchment Scale Model (DiCaSM), MODFLOW (96 and 2000) and SWI models. In addition to estimating all water balance components, DiCaSM, produces the recharge data that are used as input to the groundwater flow model of the US Geological Survey, MODFLOW. The latter subsequently generates the head distribution and groundwater flows that are used as input to the SWI model, SWI. Thus, any changes in land use, rainfall, water management, abstraction, etc. at the surface are first handled by DiCaSM, then by MODFLOW and finally by the SWI. The three models operate at different spatial and temporal scales and a facility (interface utilities between models) to aggregate/disaggregate input/output data to meet a desired spatial and temporal scale was developed allowing smooth and easy communication between the three models. As MODFLOW and SWI are published and in the public domain, this article focuses on DiCaSM, the newly developed unsaturated zone DiCaSM and equally important the interfacing utilities between the three models. DiCaSM simulates a number of hydrological processes: rainfall interception, evapotranspiration, surface runoff, infiltration, soil water movement in the root zone, plant water uptake, crop growth, stream flow and groundwater recharge. Input requirements include distributed data sets of rainfall, land use, soil types and digital terrain; climate data input can be either distributed or non‐distributed. The model produces distributed and time series output of all water balance components including potential evapotranspiration, actual evapotranspiration, rainfall interception, infiltration, plant water uptake, transpiration, soil water content, soil moisture (SM) deficit, groundwater recharge rate, stream flow and surface runoff. This article focuses on details of the hydrological processes and the various equations used in DiCaSM, as well as the nature of the interface to the MODFLOW and SWI models. Furthermore, the results of preliminary tests of DiCaSM are reported; these include tests related to the ability of the model to predict the SM content of surface and subsurface soil layers, as well as groundwater levels. The latter demonstrates how the groundwater recharge calculated from DiCaSM can be used as input into the groundwater model MODFLOW using aggregation and disaggregation algorithms (built into the interface utility). SWI has also been run successfully with hypothetical examples and was able to reproduce the results of some of the original examples of Bakker and Schaars ( 2005 ). In the subsequent articles, the results of applications to different catchments will be reported. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

Drought hotspot analysis and risk assessment using probabilistic drought monitoring and severity–duration–frequency analysis

Drought hotspot identification requires continuous drought monitoring and spatial risk assessment. The present study analysed drought events in the agriculture‐dominated mid‐Mahanadi River Basin in Odisha, India, using crop water stress as a drought indicator. This drought index incorporated different factors that affect crop water deficit such as the cropping pattern, soil characteristics, and surface soil moisture. The drought monitoring framework utilized a relevance vector machine model‐based classification that provided the uncertainty associated with drought categorization. Using the proposed framework, drought hotspots are identified in the study region and compared with indices based on precipitation and soil moisture. Further, a bivariate copula is employed to model the agricultural drought characteristics and develop the drought severity–duration–frequency (S–D–F) relationships. The drought hotspot maps and S–D–F curves are developed for different locations in the region. These provided useful information on the site‐specific drought patterns and the characteristics of the devastating droughts of 2002 and 2012, characterized by an average drought duration of 7 months at several locations. The site‐specific risk of short‐ and long‐term agricultural droughts are then investigated using the conditional copula. The results suggest that the conditional return periods and the S–D–F curves are valuable tools to assess the spatial variability of drought risk in the region.     

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.     

Subsurface lateral flow generation in aspen and conifer‐dominated hillslopes of a first order catchment in northern Utah

Mountain headwater catchments in the semi‐arid Intermountain West are important sources of surface water because these high elevations receive more precipitation than neighboring lowlands. This study examined subsurface runoff in two hillslopes, one aspen dominated, the other conifer dominated, adjacent to a first order stream in snow‐driven northern Utah. Snow accumulation, soil moisture, trenchflow and streamflow were examined in hillslopes and their adjacent stream. Snow water equivalents (SWEs) were greater under aspen stands compared to conifer, the difference increasing with higher annual precipitation. Semi‐variograms of shallow spatial soil moisture patterns and transects of continuous soil moisture showed no increase in soil moisture downslope, suggesting the absence of subsurface flow in shallow (～12 cm) soil layers of either vegetation type. However, a clear threshold relationship between soil moisture and streamflow indicated hillslope–stream connectivity, deeper within the soil profile. Subsurface flow was detected at ～50 cm depth, which was sustained for longer in the conifer hillslope. Soil profiles under the two vegetation types varied, with deep aspen soils having greater water storage capacity than shallow rocky conifer soils. Though SWEs were less under the conifers, the soil profile had less water storage capacity and produced more subsurface lateral flow during the spring snowmelt. Copyright © 2010 John Wiley & Sons, Ltd.