Groundwater–surface water interactions in a large semi‐arid floodplain: implications for salinity management

Flow regulation and water diversion for irrigation have considerably impacted the exchange of surface water between the Murray River and its floodplains. However, the way in which river regulation has impacted groundwater–surface water interactions is not completely understood, especially in regards to the salinization and accompanying vegetation dieback currently occurring in many of the floodplains. Groundwater–surface water interactions were studied over a 2 year period in the riparian area of a large floodplain (Hattah–Kulkyne, Victoria) using a combination of piezometric surface monitoring and environmental tracers (Cl⁻, δ²H, and δ¹⁸O). Despite being located in a local and regional groundwater discharge zone, the Murray River is a losing stream under low flow conditions at Hattah–Kulkyne. The discharge zone for local groundwater, regional groundwater and bank recharge is in the floodplain within ∼1 km of the river and is probably driven by high rates of transpiration by the riparian Eucalyptus camaldulensis woodland. Environmental tracers data suggest that the origin of groundwater is principally bank recharge in the riparian zone and a combination of diffuse rainfall recharge and localized floodwater recharge elsewhere in the floodplain. Although the Murray River was losing under low flows, bank discharge occurred during some flood recession periods. The way in which the water table responded to changes in river level was a function of the type of stream bank present, with point bars providing a better connection to the alluvial aquifer than the more common clay‐lined banks. Understanding the spatial variability in the hydraulic connection with the river channel and in vertical recharge following inundations will be critical to design effective salinity remediation strategies for large semi‐arid floodplains. Copyright © 2005 John Wiley & Sons, Ltd.     

A tracer‐based simulation approach to quantify seasonal dynamics of surface‐groundwater interactions in the Pantanal wetland

The Pantanal wetland is one of the least explored regions of South America. It is characterized by an outstanding flora and fauna adapted to a seasonal flood pulse controlled by a dry and a wet season within each year. The resulting inundation covers in average an area of approximately 150 000 km² and is seen as the most important driver for ecological integrity. Evaporation from the large floodplain is supposed to influence the climate of the whole continent. The regional groundwater is connected to the surface water and plays an important role for the characteristic flooding regime by regulating the wetland's water table. The water balance assessment of the wetland and the internal water exchange between surface and groundwater is therefore of high relevance for the conservation of the Pantanal biodiversity. Despite of its importance, water balance studies including groundwater–surface water interactions based on field data are rarely undertaken. This is mainly due to the remoteness and difficulty in accessing this area, which results in lack of data. In our study, we developed a new tracer‐based model to simulate the spatio–temporal surface and subsurface fluxes for a range of water bodies. The model was able to simulate these fluxes considering a dynamic simulation of inflow and outflow using a newly collected 2‐year dataset of water levels, stable water isotopes and chloride collected from several water bodies in the northern Pantanal region. Quantitative differences between water bodies according to their location in the floodplain were determined by the flooding regime and connectivity as well as site‐specific characteristics, such as hydraulic conductivity and water depth. Our model simulated water balance fluxes with a Nash–Sutcliffe efficiency of 0.61, whereas it simulated stable water isotopic compositions better than chloride. We present the first study based on field data for the Pantanal, which is able to quantify water balances fluxes. Because their representation in global climate and land cover products is insufficient, our simulation results are valuable for validating large‐scale models. Copyright © 2016 John Wiley & Sons, Ltd.     

Field assessment of surface water–groundwater connectivity in a semi‐arid river basin (Murray–Darling,Australia)

In semi‐arid and arid river basins, understanding the connectivity between rivers and alluvial aquifers is one of the key challenges for the management of groundwater resources. The type of connection present (gaining, losing‐connected, transitional and losing‐disconnected) was assessed at 12 sites along six Murray–Darling Basin river reaches. The assessments were made by measuring the hydraulic head in the riparian zone near the rivers to evaluate if the water tables intersected the riverbeds and by measuring fluid pressure (ψ) in the riverbeds. The rationale for the latter was that ψ will always be greater than or equal to zero under connected conditions (either losing or gaining) and always lesser than or equal to zero under losing‐disconnected conditions. A mixture of losing‐disconnected, losing‐connected and gaining conditions was found among the 12 sites. The losing‐disconnected sites all had a riverbed with a lower hydraulic conductivity than the underlying aquifer, usually in the form of a silty clay or clay unit 0.5–2 m in thickness. The riparian water tables were 6 to 25 m below riverbed level at the losing‐disconnected sites but never lower than 1 m below riverbed level at the losing‐connected ones. The contrast in water table depth between connected and disconnected sites was attributed to the conditions at the time of the study, when a severe regional drought had generated a widespread decline in regional water tables. This decline was apparently compensated near losing‐connected rivers by increased infiltration rates, while the decline could not be compensated at the losing‐disconnected rivers because the infiltration rates were already maximal there. Copyright © 2012 John Wiley & Sons, Ltd.     

Quantifying groundwater–surface water interactions in a proglacial valley,Cordillera Blanca,Peru

A myriad of downstream communities and industries rely on streams fed by both groundwater discharge and glacier meltwater draining the Cordillera Blanca, Northern Peruvian Andes, which contains the highest density of glaciers in the tropics. During the dry season, approximately half the discharge in the region's proglacial streams comes from groundwater. However, because of the remote and difficult access to the region, there are few field methods that are effective at the reach scale to identify the spatial distribution of groundwater discharge. An energy balance model, Rhodamine WT dye tracing, and high‐definition kite‐borne imagery were used to determine gross and net groundwater inputs to a 4‐km reach of the Quilcay River in Huascaran National Park, Peru. The HFLUX computer programme ( http://hydrology.syr.edu/hflux.html ) was used to simulate the Quilcay River's energy balance using stream temperature observations, meteorological measurements, and kite‐borne areal photography. Inference from the model indicates 29% of stream discharge at the reach outlet was contributed by groundwater discharge over the study section. Rhodamine WT dye tracing results, coupled with the energy balance, show that approximately 49% of stream water is exchanged (no net gain) with the subsurface as gross gains and losses. The results suggest that gross gains from groundwater are largest in a moraine subreach but because of large gross losses, net gains are larger in the meadow subreaches. These insights into pathways of groundwater–surface water interaction can be applied to improve hydrological modelling in proglacial catchments throughout South America. Copyright © 2016 John Wiley & Sons, Ltd.     

Groundwater–surface water interactions in a lowland watershed: source contribution to stream flow

The lower coastal plain of the Southeast USA is undergoing rapid urbanisation as a result of population growth. Land use change has been shown to affect watershed hydrology by altering stream flow and, ultimately, impairing water quality and ecologic health. However, because few long‐term studies have focused on groundwater–surface water interactions in lowland watersheds, it is difficult to establish what the effect of development might be in the coastal plain region. The objective of this study was to use an innovative improvement to end‐member mixing analysis (EMMA) to identify time sequences of hydrologic processes affecting storm flow. Hydrologic and major ion chemical data from groundwater, soil water, precipitation and stream sites were collected over a 2‐year period at a watershed located in USDA Forest Service's Santee Experimental Forest near Charleston, South Carolina, USA. Stream flow was ephemeral and highly dependent on evapotranspiration rates and rainfall amount and intensity. Hydrograph separation for a series of storm events using EMMA allowed us to identify precipitation, riparian groundwater and streambed groundwater as main sources to stream flow, although source contribution varied as a function of antecedent soil moisture condition. Precipitation, as runoff, dominated stream flow during all storm events while riparian and streambed groundwater contributions varied and were mainly dependent on antecedent soil moisture condition. Sensitivity analyses examined the influence of 10% and 50% increases in analyte concentration on EMMA calculations and found that contribution estimates were very sensitive to changes in chemistry. This study has implications on the type of methodology used in traditional forms of EMMA research, particularly in the recognition and use of median end‐member water chemistry in hydrograph separation techniques. Potential effects of urban development on important hydrologic processes (groundwater recharge, interflow, runoff, etc.) that influence stream flow in these lowland watersheds were qualitatively examined. Copyright © 2011 John Wiley & Sons, Ltd.     

Life of a bifurcation in a gravel‐bed braided river

Channel bifurcation is a key element in braided rivers, determining the water and sediment distribution and hence controlling the morphological evolution. Recent theoretical and experimental findings, as well as field observations, showed that bifurcations in gravel‐bed braided rivers are often asymmetrical and highly unstable. In this paper field data are presented on a bifurcation in the Tagliamento River, northeast Italy. The planform configuration of the bifurcation and its temporal evolution was monitored by an automatic digital camera during a series of seven floods with different magnitudes. This remote sensing technique allowed a high temporal resolution (pictures were acquired every hour) that was proved to be essential in a highly dynamic system as the one considered here. Digitized maps of the channels provided information on the location of the bifurcation, the width of the anabranches, the angle between them, along with the occurrence and migration of sediment bars. Data were acquired at two different water levels, giving the possibility to compare low and high flow conditions. The monitored bifurcation is largely unstable and shows sudden changes in the water distribution, mainly driven by the bar migrating in the upstream channel and entering the distributaries. A relationship between width asymmetry and flood magnitude was observed, confirming previous analyses. Moreover, recent theoretical findings were applied, in order to test the possibility to estimate general trends in bifurcation evolution. The analysis pointed out the relevance of a correct assessment of the characteristic temporal scales, as the bifurcation evolves on a timescale similar to that of bar migration and flood duration. Understanding the interactions between these processes is therefore crucial in order to increase the ability to model and predict the morphological evolution of a braided network. Copyright © 2012 John Wiley & Sons, Ltd.     

Groundwater–surface water interactions,nutrient fluxes and ecological response in river corridors: Translating science into effective environmental management

A summary is provided of the second in a series of Integrated Science Initiative workshops supported by the UNESCO International Hydrological Programme. The workshop brought together hydrologists, ecologists, biogeochemists, hydrogeologists and natural resource managers to discuss the processes that occur in hyporheic and riparian ecotones. The principal objectives were to share new ideas on the importance of biogeochemical processes that affect nutrients at the groundwater–surface water interface, to understand the impact of nutrient flux on stream (principally hyporheic) ecology, and to identify the management strategies for river corridors to mitigate the effects of nutrients applied to land and discharged via groundwater into rivers. The workshop concluded that: (1) more interdisciplinary research and environmental management practices are needed to better understand, predict and manage processes at the interface of environmental compartments; (2) the goal of environmental regulations to improve ecological health requires a holistic approach integrating our understanding of the ecological, hydrological, biogeochemical and physical processes; (3) upscaling spatially and temporally variable processes remains difficult and may hinder translation of research at micro‐scales (molecular to grain size) into macro‐scale (reach to catchment) decision‐making; (4) scientists need to better communicate existing research to river managers, while smanagers must better communicate policy and regulatory‐driven science requirements to researchers. Existing models, such as those that simulate stream–hyporheic exchange, are not widely known and rarely used by environmental managers. Copyright © 2007 John Wiley & Sons, Ltd.     

A tracer‐based assessment of hydrological pathways at different spatial scales in a mesoscale Scottish catchment

Geochemically based hydrograph separation techniques were used in a preliminary assessment to infer how runoff processes change with landscape characteristics and spatial scale (1–233 km²) within a mesoscale catchment in upland Scotland. A two‐component end‐member mixing analysis (EMMA) used Gran alkalinity as an assumed conservative tracer. Analysis indicated that, at all scales investigated, acidic overland flow and shallow subsurface storm flows from the peaty soils covering the catchment headwaters dominated storm runoff generation. The estimated groundwater contribution to annual runoff varied from 30% in the smallest (ca 1 km²) peat‐dominated headwater catchment with limited groundwater storage, to >60% in larger catchments (>30 km²) with greater coverage of more freely draining soils and more extensive aquifers in alluvium and other drift. This simple approach offers a useful, integrated conceptualization of the hydrological functioning in a mesoscale catchment, which can be tested and further refined by focused modelling and process‐based research. However, even as it stands, the simple conceptualization of system behaviour will have significant utility as a tool for communicating hydrological issues in a range of planning and management decisions. Copyright © 2002 John Wiley & Sons, Ltd.     

Effect of temporal resolution of water level and temperature inputs on numerical simulation of groundwater–surface water flux exchange in a heavily modified urban river

Groundwater interacts with surface water features nearly in all types of landscapes. Understanding these interactions has practical consequences on the quantity and quality of water in either system, because the depletion or contamination of one of the systems will eventually affect the other one. Many studies have shown that the use of heat as natural tracer in conjunction with water level measurements is an effective method for estimating water flow (fluxes) between groundwater and surface water. A number of studies have explored the effects of spatial and temporal variability of groundwater–surface water flux exchanges using temperature and water level measurements; however, the effect of temporal resolution of water level and temperature data on estimating flux remains unexplored. Therefore, this study investigated the effect of temporal resolution of input data on temporal variation of groundwater–surface water flux exchanges. To this end, we calibrated a variably saturated two‐dimensional groundwater flow and heat transport model (VS2DH) at hourly and daily time scales using temperatures measured at multiple depths below the riverbed of the Zenne River, located at a well‐known Belgian brownfield site. Results of the study showed that the computed water flux through the streambed ranged between ?32 mm/day and +25 mm/day using the hourly model and from ?10 mm/day to ?37 mm/day using the daily model. The hourly model resulted in detecting reversal of flow direction inducing short‐term surface water flow into the streambed. However, such events were not captured if daily temperature and water level measurements were used as input. These findings have important implications for understanding contaminant mass flux and their attenuation in the mixing zone of groundwater and surface water. Copyright © 2012 John Wiley & Sons, Ltd.     

Heterogeneity in ground water–surface water interactions in the hyporheic zone of a salmonid spawning stream

Spatial and temporal variability in ground water–surface water interactions in the hyporheic zone of a salmonid spawning stream was investigated. Four locations in a 150‐m reach of the stream were studied using hydrometric and hydrochemical tracing techniques. A high degree of hydrological connectivity between the riparian hillslope and the stream channel was indicated at two locations, where hydrochemical changes and hydraulic gradients indicated that the hyporheic zone was dominated by upwelling ground water. The chemistry of ground water reflected relatively long residence times and reducing conditions with high levels of alkalinity and conductivity, low dissolved oxygen (DO) and nitrate. At the other locations, connectivity was less evident and, at most times, the hyporheic zone was dominated by downwelling stream water characterized by high DO, low alkalinity and conductivity. Substantial variability in hyporheic chemistry was evident at fine (<10 m) spatial scales and changed rapidly over the course of hydrological events. The nature of the hydrochemical response varied among locations depending on the strength of local ground water influence. It is suggested that greater emphasis on spatial and temporal heterogeneity in ground water–surface water interactions in the hyporheic zone is necessary for a consideration of hydrochemical effects on many aspects of stream ecology. For example, the survival of salmonid eggs in hyporheic gravels varied considerably among the locations studied and was shown to be associated with variation in interstitial chemistry. River restoration schemes and watershed management strategies based only on the surface expression of catchment characteristics risk excluding consideration of potentially critical subsurface processes. Copyright © 2002 John Wiley & Sons, Ltd.     

River–aquifer interactions in a semi‐arid environment stressed by groundwater abstraction

Rivers and aquifers are, in many cases, a connected resource and as such the interactions between them need to be understood and quantified for the resource to be managed appropriately. The objective of this paper is to advance the understanding of river–aquifer interactions processes in semi‐arid environments stressed by groundwater abstraction. This is performed using data from a specific catchment where records of precipitation, evapotranspiration, river flow, groundwater levels and groundwater abstraction are analysed using basic statistics, hydrograph analysis and a simple mathematical model to determine the processes causing the spatial and temporal changes in river–aquifer interactions. This combined approach provides a novel but simple methodology to analyse river–aquifer interactions, which can be applied to catchments worldwide. The analysis revealed that the groundwater levels have declined (~ 3 m) since the onset of groundwater abstraction. The decline is predominantly due to the abstraction rather than climatic changes (r = 0.84 for the relationship between groundwater abstraction and groundwater levels; r = 0.92 for the relationship between decline in groundwater levels and magnitude of seasonal drawdown). It is then demonstrated that, since the onset of abstraction, the river has changed from being gaining to losing during low‐flow periods, defined as periods with flow less than 0.5, 1.0 or 1.5 GL/day (1 GL/day = 1 × 10⁶ m³/day). If defined as < 1.0 GL/day, low‐flow periods constitute approximately 65% of the river flows; the periods where the river is losing at low‐flow conditions are thus significant. Importantly, there was a significant delay (> 10 years) between the onset of groundwater abstraction and the changeover from gaining to losing conditions. Finally, a relationship between the groundwater gradient towards the river and the river flow at low‐flow is demonstrated. The results have important implications for water management as well as water ecology and quality. Copyright © 2012 John Wiley & Sons, Ltd.     

Modelling coastal ground‐ and surface‐water interactions using an integrated approach

Beach water table fluctuations have an impact on the transport of beach sediments and the exchange of solute and mass between coastal aquifer and nearby water bodies. Details are given of the refinement of a dynamically integrated ground‐ and surface‐water model, and its application to study ground‐ and surface‐water interactions in coastal regions. The depth‐integrated shallow‐water equations are used to represent the surface‐water flow, and the extended Darcy's equation is used to represent the groundwater flow, with a hydrostatic pressure distribution being assumed to apply for both these two types of flows. At the intertidal region, the model has two layers, with the surface‐water layer being located on the top of the groundwater layer. The governing equations for these two types of flows are discretized in a similar manner and they are combined to give one set of linear algebraic equations that can be solved efficiently. The model is used to predict water level distributions across sloping beaches, where the water table in the aquifer may or may not decouple from the free water surface. Five cases are used to test the model for simulating beach water table fluctuations induced by tides, with the model predictions being compared with existing analytical solutions and laboratory and field data published in the literature. The numerical model results show that the integrated model is capable of simulating the combined ground‐ and surface‐water flows in coastal areas. Detailed analysis is undertaken to investigate the capability of the model. Copyright © 2009 John Wiley & Sons, Ltd.     

From streambed temperature measurements to spatial‐temporal flux quantification: using the LPML method to study groundwater–surface water interaction

Knowledge on groundwater–surface water interaction and especially on exchange fluxes between streams and aquifers is an important prerequisite for the study of transport and fate of contaminants and nutrients in the hyporheic zone. One possibility to quantify groundwater–surface water exchange fluxes is by using heat as an environmlental tracer. Modern field equipment including multilevel temperature sticks and the novel open‐source analysis tool LPML make this technique ever more attractive. The recently developed LPML method solves the one‐dimensional fluid flow and heat transport equation by combining a local polynomial method with a maximum likelihood estimator. In this study, we apply the LPML method on field data to quantify the spatial and temporal variability of vertical fluxes and their uncertainties from temperature–time series measured in a Belgian lowland stream. Over several months, temperature data were collected with multilevel temperature sticks at the streambed top and at six depths for a small stream section. Long‐term estimates show a range from gaining fluxes of ?291 mm day^?1 to loosing fluxes of 12 mm day^?1; average seasonal fluxes ranged from ?138 mm day^?1 in winter to ?16 mm day^?1 in summer. With our analyses, we could determine a high spatial and temporal variability of vertical exchange fluxes for the investigated stream section. Such spatial and temporal variability should be taken into account in biogeochemical cycling of carbon, nutrients and metals and in fate analysis of contaminant plumes. In general, the stream section was gaining during most of the observation period. Two short‐term high stream stage events, seemingly caused by blockage of the stream outlet, led to a change in flow direction from gaining to losing conditions. We also found more discharge occurring at the outer stream bank than at the inner one indicating a local flow‐through system. With the conducted analyses, we were able to advance our understanding of the regional groundwater flow system. Copyright © 2015 John Wiley & Sons, Ltd.     

Scales and magnitude of hyporheic,river–aquifer and bank storage exchange fluxes

Many studies have investigated the exchange processes that occur between rivers and groundwater systems and have successfully quantified the water fluxes involved. Specifically, these exchange processes include hyporheic exchange, river–aquifer exchange (groundwater discharge and river loss) and bank storage exchange. Remarkably, there are relatively few examples of field studies where more than one exchange process is quantified, and as a consequence, the relationships between them are not well understood. To compare the relative magnitudes of these common exchange processes, we have collected data from 54 studies that have quantified one or more of these exchange flux types. Each flux value is plotted against river discharge at the time of measurement to allow the different exchange flux types to be compared. We show that there are positive relationships between the magnitude of each exchange flux type and increasing river discharge across the different studies. For every one order of magnitude increase in river discharge, the hyporheic, river–aquifer and bank storage exchange fluxes increase by factors of 2.7, 2.9 and 2.5, respectively. On average, hyporheic exchange fluxes are almost an order of magnitude greater than river–aquifer exchange fluxes, which are, in turn, approximately four times greater than bank storage exchange fluxes for the same river discharge. Unless measurement approaches that can distinguish between different types of exchange flux are used, there is potential for hyporheic exchange fluxes to be misinterpreted as river–aquifer exchange fluxes, with possible implications for water resource management decisions. Copyright © 2015 John Wiley & Sons, Ltd.     

Using new mass balance methods to estimate gross surface water and groundwater exchange with naturally occurring tracer 222Rn in data poor regions: a case study in northwest China

Given that the concentration of ²²²Rn in groundwater is much higher than that in surface water and that its radioactive half‐life (3.83 d) is short, ²²²Rn is an effective tracer of groundwater–surface water interactions. In this study, a new mass balance method is presented, which can be used to estimate specific groundwater–surface water interactions within a river reach. Three possible situations of interaction between groundwater and surface water are considered, and equations based on the mass conservation of ²²²Rn are formulated for judging specific groundwater–surface water interaction processes and for calculating water flux. A case study was conducted for the Nalenggele River, Northwest China, to demonstrate the usefulness of this method. Samples of river water and groundwater containing ²²²Rn were collected from the study area to estimate the interactions between groundwater and surface water. The amount of water exchanged during these interactions was estimated and the results show that transformations between groundwater and surface water are frequent along the stream. The ²²²Rn mass balance method is highly sensitive for studying such interactions, even in areas for which conventional hydrologic data are sparse. Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

Groundwater and surface water connectivity within the recharge area of Guarani aquifer system during El Niño 2014–2016

Recharge areas of the Guarani Aquifer System (GAS) are particularly sensitive and vulnerable to climate variability; therefore, the understanding of infiltration mechanisms for aquifer recharge and surface run‐off generation represent a relevant issue for water resources management in the southeastern portion of the Brazilian territory, particularly in the Jacaré‐Pepira River watershed. The main purpose of this study is to understand the interactions between precipitation, surface water, and groundwater using stable isotopes during the strong 2014–2016 El Niño Southern Oscillation event. The large variation in the isotopic composition of precipitation (from ?9.26‰ to +0.02‰ for δ¹⁸O and from ?63.3‰ to +17.6‰ for δ²H), mainly associated with regional climatic features, was not reflected in the isotopic composition of surface water (from ?7.84‰ to ?5.83‰ for δ¹⁸O and from ?49.7‰ to +33.6‰ for δ²H), mainly due to the monthly sampling frequency, and groundwater (from ?7.04‰ to ?7.76‰ for δ¹⁸O and from ?49.5‰ to ?44.7‰ for δ²H), which exhibited less variation throughout the year. However, variations in deuterium excess (d‐excess) in groundwater and surface water suggest the occurrence of strong secondary evaporation during the infiltration process, corresponding with groundwater level recovery. Similar isotopic composition in groundwater and surface water, as well as the same temporal variations in d‐excess and line‐conditioned excess denote the strong connectivity between these two reservoirs during baseflow recession periods. Isotopic mass balance modelling and hydrograph separation estimate that the groundwater contribution varied between 70% and 80%, however, during peak flows, the isotopic mass balance tends to overestimate the groundwater contribution when compared with the other hydrograph separation methods. Our findings indicate that the application of isotopic mass balance methods for ungauged rivers draining large groundwater reservoirs, such as the GAS outcrop, could provide a powerful tool for hydrological studies in the future, helping in the identification of flow contributions to river discharge draining these areas.     

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.     

Studies of a regulated dryland river: surface–groundwater interactions

Unlike rivers in humid regions, dryland rivers typically exhibit reduced flow in the downstream direction as a result of transmission losses, which include seepage of streamflow into the aquifer, evaporation, and transpiration. However, much remains to be learned about the nature of the exchange between surface water and groundwater in these landscapes, especially in terms of spatial and temporal variability. Our study focused on streambank seepage and groundwater flow in the alluvial aquifer, specifically on answering questions such as: Is there seasonal variability in seepage losses? Is seepage permanently lost? Can losses be reduced by killing riparian vegetation? To better understand the magnitude, variability, and fate of streambank seepage, we assessed river stages, groundwater hydraulic gradients, and groundwater flow paths at two sites along a reach of the Pecos River, a dryland perennial river in West Texas. We found that along this reach the river was losing water to the aquifer even under low‐flow conditions; but seepage was controlled by a number of different mechanisms. Seepage increased not only during high‐flow events but also when the groundwater level was declining owing to long periods of no irrigation release. Tamarix (saltcedar) control did not affect hydraulic gradients nor reduce streambank seepage and given that this reach of the Pecos River is a losing one, streamflow will not be enhanced by controlling saltcedar. These findings can be used to improve basic conceptual models of dryland river systems and to predict hydrologic responses to changes in the timing and magnitude of streamflows and to riparian vegetation management. Copyright © 2012 John Wiley & Sons, Ltd.