A low‐dimensional hillslope‐based catchment model for layered groundwater flow

Despite the strong interaction between surface and subsurface waters, groundwater flow representation is often oversimplified in hydrological models. For instance, the interplay between local or shallow aquifers and deeper regional‐scale aquifers is typically neglected. In this work, a novel hillslope‐based catchment model for the simulation of combined shallow and deep groundwater flow is presented. The model consists of the hillslope‐storage Boussinesq (hsB) model representing shallow groundwater flow and an analytic element (AE) model representing deep regional groundwater flow. The component models are iteratively coupled via a leakage term based on Darcy's law, representing delayed recharge to the regional aquifer through a low conductivity layer. Simulations on synthetic single hillslopes and on a two‐hillslope open‐book catchment are presented, and the results are compared against a benchmark three‐dimensional Richards equation model. The impact of hydraulic conductivity, hillslope plan geometry (uniform, convergent, divergent), and hillslope inclination (0.2%, 5%, and 30%) under drainage and recharge conditions are examined. On the single hillslopes, good matches for heads, hydrographs, and exchange fluxes are generally obtained, with the most significant differences in outflows and heads observed for the 30% slope and for hillslopes with convergent geometry. On the open‐book catchment, cumulative outflows are overestimated by 1–4%. Heads in the confined and unconfined aquifers are adequately reproduced throughout the catchment, whereas exchange fluxes are found to be very sensitive to the hillslope drainable porosity. The new model is highly efficient computationally compared to the benchmark model. The coupled hsB/AE model represents an alternative to commonly used groundwater flow representations in hydrological models, of particular appeal when surface–subsurface exchanges, local aquifer–regional aquifer interactions, and low flows play a key role in a watershed's dynamics. Copyright © 2011 John Wiley & Sons, Ltd.     

Seepage response along an alluvial valley in a semi‐arid catchment in north‐central California

Over a period of 12 months, soil moisture content and potential was monitored in an annual‐grass‐dominated 20 ha catchment in order to determine flow paths leading to exfiltration at the catchment outlet. Water was found to enter the catchment valley either through flow originating in the slopes or through surface infiltration during rainfall events. Although subsurface flow from the slopes to the catchment outlet occurred throughout the year, surface recharge was restricted to a few events during the wet season. In the deeper saturated profile of the valley, flow was directed upwards along the valley edges and gradually became horizontal towards the central axis of the valley. During the peak of the rainfall season, horizontal flow close to the catchment outlet intercepted the gradually sloping surface, resulting in exfiltration. Plants influenced the hydrology of the catchment by removing moisture from the root zone during spring and early summer, resulting in evapotranspiration losses from the vadose zone. Heterogeneities within the valley soil were evident as variable‐permeability layers that resulted in a seasonally confined water table within the valley. This investigation shows that the vadose zone plays an important role in redistributing surface recharge and emphasizes the importance of accounting for effective moisture in low‐yielding catchments with ephemeral surface runoff. Copyright © 2002 John Wiley & Sons, Ltd.     

Spatial and temporal analysis of groundwater recharge with application to sampling design

A numerical experiment of flow in variably saturated porous media was performed in order to evaluate the spatial and temporal distribution of the groundwater recharge at the phreatic surface for a shallow aquifer as a function of the input rainfall process and soil heterogeneity. The study focused on the groundwater recharge which resulted from the percolation of the excess rainfall for a 90-days period of an actual precipitation record. Groundwater recharge was defined as the water flux across the moving phreatic surface. The observed spatial non-uniformity of the groundwater recharge was caused by soil heterogeneity and is particularly pronounced during the stage of recharge peak (substantial percolation stage). During that stage the recharge is associated with preferential flow paths defined as soil zones of locally higher hydraulic conductivity. For the periods of low percolation intensity the groundwater recharge was exhibiting more uniform spatial characteristics. The temporal distribution of the recharge was found to be a function of the frequency and intensity of the rainfall events. Application of sampling design demonstrates the joint influence of the spatial and temporal recharge variability on the cost-effective monitoring of groundwater potentiometric surfaces.     

Studying the Flow Dynamics of a Karst Aquifer System with an Equivalent Porous Medium Model

The modeling of groundwater flow in karst aquifers is a challenge due to the extreme heterogeneity of its hydraulic parameters and the duality in their discharge behavior, that is, rapid response of highly conductive karst conduits and delayed drainage of the low‐permeability fractured matrix after recharge events. There are a number of different modeling approaches for the simulation of the karst groundwater dynamics, applicable to different aquifer as well as modeling problem types, ranging from continuum models to double continuum models to discrete and hybrid models. This study presents the application of an equivalent porous model approach (EPM, single continuum model) to construct a steady‐state numerical flow model for an important karst aquifer, that is, the Western Mountain Aquifer Basin (WMAB), shared by Israel and the West‐Bank, using MODFLOW2000. The WMAB was used as a catchment since it is a well‐constrained catchment with well‐defined recharge and discharge components and therefore allows a control on the modeling approach, a very rare opportunity for karst aquifer modeling. The model demonstrates the applicability of equivalent porous medium models for the simulation of karst systems, despite their large contrast in hydraulic conductivities. As long as the simulated saturated volume is large enough to average out the local influence of karst conduits and as long as transport velocities are not an issue, EPM models excellently simulate the observed head distribution. The model serves as a starting basis that will be used as a reference for developing a long‐term dynamic model for the WMAB, starting from the pre‐development period (i.e., 1940s) up to date.     

Effects of meteorological and hydrogeological factors on gross recharge percentage at unconfined sandy aquifers with an equatorial climate

A thorough understanding of rainfall recharge processes and their controlling factors is essential for management of groundwater systems. This study investigates the effects of various meteorological and hydrogeological factors on the gross recharge percentages, the rainfall–recharge relationships and the recharge threshold values for unconfined sandy aquifers under an equatorial climate. Among the meteorological factors investigated, rainfall intensity was found to have the most significant impact on the gross recharge rate. The effects of potential evaporation rate, relative humidity and air temperature on the gross recharge percentage were significant when the vadose zone thickness is larger than 2·5 m. The recharge threshold values were found to depend strongly on the vadose zone thickness. The rainfall–recharge relationships could generally be well defined by a normal–log relationship. The rainfall–recharge relationships derived here are applicable to yield estimates of gross recharge percentages for unconfined sandy aquifers under an equatorial climate, using rainfall intensity and vadose zone thickness as input variables. In this study, a theory was developed and validated to provide physical explanations for the observations, based on the residence time of the percolated rainwater within the vadose zone. Among the soil hydraulic parameters tested, porosity and saturated hydraulic conductivity were found to have the most pronounced effects on the gross recharge percentage. Utilizing the sensitivity results and the theory derived, an approach was developed for extending the application of the derived rainfall–recharge relationships to other sand textures. The approach was found to be capable of producing rough and fast estimations of gross recharge percentage for other sand textures. Copyright © 2007 John Wiley & Sons, Ltd.     

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.     

Effect of hydrogeomorphic setting on calcareous fen hydrology

Calcareous fens are species‐rich peatlands that are dependent on minerotrophic water sources for wetland functioning, with current conceptual models suggesting the water source is ubiquitously groundwater upwelling. By quantifying the water balance and subsurface water flow paths and fluxes over 3 growing seasons for calcareous fens in 3 different hydrogeomorphic settings (Riparian, Trough, and Basin), we show evidence that challenges this conceptual model. The Riparian Fen received an order of magnitude more water inputs than the Trough or Basin Fens and was dominated by stream recharge inputs and groundwater outputs. Precipitation and evaporation dominated the water balance of the Trough Fen whereas only the Basin Fen received sizeable groundwater inputs. Indeed, subsurface water fluxes were low at all fens due to weak hydraulic gradients and low saturated hydraulic conductivity in some areas of each wetland, though variations in growing season precipitation led to subsurface flow reversals in all 3 fens. Our results demonstrate the importance of understanding landscape position, or hydrogeomorphic setting, on calcareous fen hydrology for improving conservation, management, and restoration efforts of these important ecosystems.     

Physical controls on septic leachate movement in the vadose zone at the hillslope scale,Putnam County,New York,USA

The fate and transport of contaminants in the vicinity of septic fields remains poorly understood in many hydrogeomorphological environments. We report hydrometric data from an intensive hillslope‐scale experiment conducted between 29 August and 11 November 1998 at a residential leach field in New York State. The objective of our study was to characterize water flux within the vadose zone, understand the physical controls on the flux, and predict how this ultimately will affect subsurface water quality. Soil‐water flux was calculated using matric potential measurements from a network of 25 tensiometer nests, each nest consisting of three tensiometers installed to depths of 10, 50 and 130 cm. Unsaturated hydraulic conductivity curves were derived at each depth from field‐determined time‐domain reflectometry–tensiometry moisture‐release curves and borehole permeametry measurements. Flownets indicated that a strong upward flux of soil water occurred between rainstorms. Following the onset of (typically convective) rainfall, low near‐surface matric potentials were rapidly converted to near‐saturated and saturated conditions, promoting steep vertical gradients through the near‐surface horizons of the hillslope. Lateral hydraulic gradients were typically 10 times smaller than the vertical gradients. Resultant flow vectors showed that the flux was predominantly vertical through the vadose zone, and that the flux response to precipitation was short‐lived. The flux response was controlled primarily by the shape of the unsaturated hydraulic conductivity curves, which indicated a rapid loss of conductivity below saturation. Thus, soil water had a very high residence time in the vadose zone. The absence of rapid wetting at 130 cm and the delayed and small phreatic zone response to rainfall indicated that water movement through macropores did not occur on this hillslope. These results are consistent with a Cl tracing experiment, which demonstrated that the tracer was retained in the vadose zone for several months after injection to the system. Copyright © 2002 John Wiley & Sons, Ltd.     

A method for estimating spatially variable seepage and hydraulic conductivity in channels with very mild slopes

Infiltration along ephemeral channels plays an important role in groundwater recharge in arid regions. A model is presented for estimating spatial variability of seepage due to streambed heterogeneity along channels based on measurements of streamflow‐front velocities in initially dry channels. The diffusion‐wave approximation to the Saint‐Venant equations, coupled with Philip's equation for infiltration, is connected to the groundwater model MODFLOW and is calibrated by adjusting the saturated hydraulic conductivity of the channel bed. The model is applied to portions of two large water delivery canals, which serve as proxies for natural ephemeral streams. Estimated seepage rates compare well with previously published values. Possible sources of error stem from uncertainty in Manning's roughness coefficients, soil hydraulic properties and channel geometry. Model performance would be most improved through more frequent longitudinal estimates of channel geometry and thalweg elevation, and with measurements of stream stage over time to constrain wave timing and shape. This model is a potentially valuable tool for estimating spatial variability in longitudinal seepage along intermittent and ephemeral channels over a wide range of bed slopes and the influence of seepage rates on groundwater levels. Copyright © 2012 John Wiley & Sons, Ltd.     

A scale‐invariant property of the water retention curve in weakly heterogeneous vadose zones

Abrupt changes of hydraulic properties in a vadose zone are modelled within a stochastic framework, which regards the saturated conductivity and parameters related to the distribution of soil pores as stationary, log‐normally distributed, random space functions. As a consequence, flow variables become random fields, and we aim at deriving an effective Richards equation. To obtain the latter, we adopt a perturbation expansion truncated at the first order (weakly heterogeneous media), which leads to the effective hydraulic conductivity and water retention curves. Overall, the effective properties are scale dependent. However, within the proposed framework, we demonstrate that the inflection point of the laboratory scale retention curve is not affected by the heterogeneity of the vadose zone. Finally, to illustrate the quantitative implications of our results, we consider a monitoring experiment at field scale, and we show how our approach leads to an effective water retention curve, which differs significantly from that which would be obtained without accounting for the above scale‐invariance property.     

Spatial and temporal infiltration dynamics during managed aquifer recharge

Natural groundwater recharge is inherently difficult to quantify and predict, largely because it comprises a series of processes that are spatially distributed and temporally variable. Infiltration ponds used for managed aquifer recharge (MAR) provide an opportunity to quantify recharge processes across multiple scales under semi-controlled conditions. We instrumented a 3-ha MAR infiltration pond to measure and compare infiltration patterns determined using whole-pond and point-specific methods. Whole-pond infiltration was determined by closing a transient water budget (accounting for inputs, outputs, and changes in storage), whereas point-specific infiltration rates were determined using heat as a tracer and time series analysis at eight locations in the base of the pond. Whole-pond infiltration, normalized for wetted area, rose rapidly to more than 1.0 m/d at the start of MAR operations (increasing as pond stage rose), was sustained at high rates for the next 40 d, and then decreased to less than 0.1 m/d by the end of the recharge season. Point-specific infiltration rates indicated high spatial and temporal variability, with the mean of measured values generally being lower than rates indicated by whole-pond calculations. Colocated measurements of head gradients within saturated soils below the pond were combined with infiltration rates to calculate soil hydraulic conductivity. Observations indicate a brief period of increasing saturated hydraulic conductivity, followed by a decrease of one to two orders of magnitude during the next 50 to 75 d. Locations indicating the most rapid infiltration shifted laterally during MAR operation, and we suggest that infiltration may function as a "variable source area" processes, conceptually similar to catchment runoff.     

An approach to catchment‐scale groundwater nitrate risk assessment from diffuse agricultural sources: a case study in the Upper Bann,Northern Ireland

DRASTIC has drawbacks in groundwater risk assessment that are important in guiding activities to prevention agricultural diffuse groundwater pollution. This paper presents an improved and GIS‐based D‐DRASTIC approach for groundwater nitrate risk assessment from diffuse agricultural sources based on DRASTIC. D‐DRASTIC considers the risk concept, nitrate loading, pollutant transport with runoff, depth to water, net recharge, aquifer media, soil media, topography, impact of the vadose zone media, and the hydraulic conductivity of the aquifer. D‐DRASTIC was developed within an ArcGIS environment and applied to the Upper Bann Catchment, Northern Ireland as a case study. D‐DRASTIC shows that ‘very high’ and ‘high’ zones of groundwater nitrate risk occupy 5% and 11% of the case study area, respectively. When considering groundwater pollution sources and pathways, the results using D‐DRASTIC are helpful in guiding the activities of groundwater pollution prevention at the catchment scale in the context of better implementation of the EU Water Framework Directive. Copyright © 2008 John Wiley & Sons, Ltd.     

The Effect of Modeled Recharge Distribution on Simulated Groundwater Availability and Capture

Simulating groundwater flow in basin‐fill aquifers of the semiarid southwestern United States commonly requires decisions about how to distribute aquifer recharge. Precipitation can recharge basin‐fill aquifers by direct infiltration and transport through faults and fractures in the high‐elevation areas, by flowing overland through high‐elevation areas to infiltrate at basin‐fill margins along mountain fronts, by flowing overland to infiltrate along ephemeral channels that often traverse basins in the area, or by some combination of these processes. The importance of accurately simulating recharge distributions is a current topic of discussion among hydrologists and water managers in the region, but no comparative study has been performed to analyze the effects of different recharge distributions on groundwater simulations. This study investigates the importance of the distribution of aquifer recharge in simulating regional groundwater flow in basin‐fill aquifers by calibrating a groundwater‐flow model to four different recharge distributions, all with the same total amount of recharge. Similarities are seen in results from steady‐state models for optimized hydraulic conductivity values, fit of simulated to observed hydraulic heads, and composite scaled sensitivities of conductivity parameter zones. Transient simulations with hypothetical storage properties and pumping rates produce similar capture rates and storage change results, but differences are noted in the rate of drawdown at some well locations owing to the differences in optimized hydraulic conductivity. Depending on whether the purpose of the groundwater model is to simulate changes in groundwater levels or changes in storage and capture, the distribution of aquifer recharge may or may not be of primary importance.     

Investigation of the thermal regime and subsurface properties of a tidally affected,variably saturated streambed

Temperature and moisture content in the variably saturated subsurface are two of the most important physical parameters that govern a wide variety of geochemical and ecological processes. An understanding of thermal and hydraulic processes and properties of transient vadose zones is therefore fundamental in the evaluation of such processes. Here, an investigation of the thermal regime and subsurface properties of a tidally affected, variably saturated streambed is presented. Field and laboratory measurements, as well as a forward numerical model, are jointly employed in the investigation. Temperature, soil moisture, surface level, and water level data were recorded in a transect perpendicular to a tidally driven stream. Frequency‐domain analysis of the subsurface temperature measurements revealed the rapid decay of the tidal temperature driver within the top ～30 cm of sediment. Several techniques were used to evaluate subsurface thermal and hydraulic properties, including thermal conductivity and the soil water retention curve. These properties were used to constrain a forward numerical model that included coupled treatment of relevant variable saturation thermal and hydraulic physics. Even though the investigated vadose zone is intermittent and relatively shallow ( 20 cm), the results illustrate how error can be introduced into heat‐transport calculations if unsaturated conditions are not taken into account.     

A numerical investigation of bedrock groundwater recharge and exfiltration on soil mantled hillslopes

Here we use Richards Equation models of variably saturated soil and bedrock groundwater flow to investigate first-order patterns of the coupling between soil and bedrock flow systems. We utilize a Monte Carlo sensitivity analysis to identify important hillslope parameters controlling bedrock recharge and then model the transient response of bedrock and soil flow to seasonal precipitation. Our results suggest that hillslopes can be divided into three conceptual zones of groundwater interaction, (a) the zone of lateral unsaturated soil moisture accumulation (upper portion of hillslope), (b) the zone of soil saturation and bedrock recharge (middle of hillslope) and (c) the zone of saturated-soil lateral flow and bedrock groundwater exfiltration (bottom of hillslope). Zones of groundwater interaction expand upslope during periods of precipitation and drain downslope during dry periods. The amount of water partitioned to the bedrock groundwater system a can be predicted by the ratio of bedrock to soil saturated hydraulic conductivity across a variety of hillslope configurations. Our modelled processes are qualitatively consistent with observations of shallow subsurface saturation and groundwater fluctuation on hillslopes studied in our two experimental watersheds and support a conceptual model of tightly coupled shallow and deep subsurface circulation where groundwater recharge and discharge continuously stores and releases water from longer residence time storage.     

Hydraulic subsurface measurements and hydrodynamic modelling as indicators for groundwater flow systems in the Rotondo granite,Central Alps (Switzerland)

Regional groundwater flow in high mountainous terrain is governed by a multitude of factors such as geology, topography, recharge conditions, structural elements such as fracturation and regional fault zones as well as man‐made underground structures. By means of a numerical groundwater flow model, we consider the impact of deep underground tunnels and of an idealized major fault zone on the groundwater flow systems within the fractured Rotondo granite. The position of the free groundwater table as response to the above subsurface structures and, in particular, with regard to the influence of spatial distributed groundwater recharge rates is addressed. The model results show significant unsaturated zones below the mountain ridges in the study area with a thickness of up to several hundred metres. The subsurface galleries are shown to have a strong effect on the head distribution in the model domain, causing locally a reversal of natural head gradients. With respect to the position of the catchment areas to the tunnel and the corresponding type of recharge source for the tunnel inflows (i.e. glaciers or recent precipitation), as well as water table elevation, the influence of spatial distributed recharge rates is compared to uniform recharge rates. Water table elevations below the well exposed high‐relief mountain ridges are observed to be more sensitive to changes in groundwater recharge rates and permeability than below ridges with less topographic relief. In the conceptual framework of the numerical simulations, the model fault zone has less influence on the groundwater table position, but more importantly acts as fast flow path for recharge from glaciated areas towards the subsurface galleries. This is in agreement with a previous study, where the imprint of glacial recharge was observed in the environmental isotope composition of groundwater sampled in the subsurface galleries. Copyright © 2012 John Wiley & Sons, Ltd.     

Characterizing the heterogeneity of karst critical zone and its hydrological function: An integrated approach

Spatial heterogeneity in the subsurface of karst environments is high, as evidenced by the multiphase porosity of carbonate rocks and complex landform features that result in marked variability of hydrological processes in space and time. This includes complex exchange of various flows (e.g., fast conduit flows and slow fracture flows) in different locations. Here, we integrate various “state‐of‐the‐art” methods to understand the structure and function of this poorly constrained critical zone environment. Geophysical, hydrometric, and tracer tools are used to characterize the hydrological functions of the cockpit karst critical zone in the small catchment of Chenqi, Guizhou Province, China. Geophysical surveys, using electrical resistivity tomography (ERT), inferred the spatial heterogeneity of permeability in the epikarst and underlying aquifer. Water tables in depression wells in valley bottom areas, as well as discharge from springs on steeper hillslopes and at the catchment outlet, showed different hydrodynamic responses to storm event rainwater recharge and hillslope flows. Tracer studies using water temperatures and stable water isotopes (δD and δ¹⁸O) could be used alongside insights into aquifer permeability from ERT surveys to explain site‐ and depth‐dependent variability in the groundwater response in terms of the degree to which “new” water from storm rainfall recharges and mixes with “old” pre‐event water in karst aquifers. This integrated approach reveals spatial structure in the karst critical zone and provides a conceptual framework of hydrological functions across spatial and temporal scales.     

Lagged rejuvenation of groundwater indicates internal flow structures and hydrological connectivity

Large proportions of rainwater and snowmelt infiltrate into the subsurface before contributing to stream flow and stream water quality. Subsurface flow dynamics steer the transport and transformation of contaminants, carbon, weathering products and other biogeochemistry. The distribution of groundwater ages with depth is a key feature of these flow dynamics. Predicting these ages are a strong test of hypotheses about subsurface structures and time-varying processes. Chlorofluorocarbon (CFC)-based groundwater ages revealed an unexpected groundwater age stratification in a 0.47 km² forested catchment called Svartberget in northern Sweden. An overall groundwater age stratification, representative for the Svartberget site, was derived by measuring CFCs from nine different wells with depths of 2–18 m close to the stream network. Immediately below the water table, CFC-based groundwater ages of already 30 years that increased with depth were found. Using complementary groundwater flow models, we could reproduce the observed groundwater age stratification and show that the 30 year lag in rejuvenation comes from return flow of groundwater at a subsurface discharge zone that evolves along the interface between two soil types. By comparing the observed groundwater age stratification with a simple analytical approximation, we show that the observed lag in rejuvenation can be a powerful indicator of the extent and structure of the subsurface discharge zone, while the vertical gradient of the age-depth-relationship can still be used as a proxy of the overall aquifer recharge even when sampled in the discharge zone. The single age stratification profile measured in the discharge zone, close to the aquifer outlet, can reveal the main structure of the groundwater flow pattern from recharge to discharge. This groundwater flow pattern provides information on the participation of groundwater in the hydrological cycle and indicates the lower boundary of hydrological connectivity.     

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.