The impact of groundwater–surface water interactions on the water balance of a mesoscale lowland river catchment in northeastern Germany

The glacially formed northeastern German lowlands are characterized by extensive floodplains, often interrupted by relatively steep moraine hills. The hydrological cycle of this area is governed by the tight interaction of surface water dynamics and the corresponding directly connected shallow groundwater aquifer. Runoff generation processes, as well as the extent and spatial distribution of the interaction between surface water and groundwater, are controlled by floodplain topography and by surface water dynamics. A modelling approach based on extensive experimental analyses is presented that describes the specific water balance of lowland areas, including the interactions of groundwater and surface water, as well as reflecting the important role of time‐variable shallow groundwater stages for runoff generation in floodplains. In the first part, experimental investigations of floodplain hydrological characteristics lead to a qualitative understanding of the water balance processes and to the development of a conceptual model of the water balance and groundwater dynamics of the study area. Thereby model requirements which allow for an adequate simulation of the floodplain hydrology, considering also interactions between groundwater and surface water have been characterized. Based on these analyses, the Integrated Modelling of Water Balance and Nutrient Dynamics (IWAN) approach has been developed. This consists of coupling the surface runoff generation and soil water routines of the deterministic, spatially distributed hydrological model WASIM‐ETH‐I with the three‐dimensional finite‐difference‐based numerical groundwater model MODFLOW and Processing MODFLOW. The model was applied successfully to a mesoscale subcatchment of the Havel River in northeast Germany. It was calibrated for two small catchments (1·4 and 25 km²), where the importance of the interaction processes between groundwater and surface waters and the sensitivity of several controlling parameters could be quantified. Validation results are satisfying for different years for the entire 198 km² catchment. The model approach was further successfully tested for specific events. The experimental area is a typical example of a floodplain‐dominated landscape. It was demonstrated that the lateral flow processes and the interactions between groundwater and surface water have a major importance for the water balance and periodically superimposed on the vertical runoff generation. Copyright © 2006 John Wiley & Sons, Ltd.     

Salinograph trends as indicators of the recession characteristics of stream components

Compared to hydrograph recession analysis, which is widely applied in engineering hydrology, the quantitative assessment of stream salinity with time (i.e. the salinograph) has received significantly less attention. In particular, while in many previous hydrological studies an inverse relationship between hydrograph and salinograph responses is apparent, the concept of salinity accession (the inversely related salinity counterpart to hydrograph recession) has not been introduced nor quantitatively evaluated in previous literature. In this study, we conduct a mathematical analysis of salinograph accession, and determine new quantitative relationships between salinity accession and hydrograph recession parameters. An equation is formulated that reproduces the general trend in salinity accession. A salinity accession parameter k_c is then introduced and is shown to be the ratio of direct runoff to total stream flow recession parameters: k_r/k. The groundwater recession parameter k_g was estimated using a simple and rapid method that uses both salinograph and hydrograph data. Salinity accession type‐curves illustrate that under certain conditions, the relative steepness of individual salinographs is dependent upon the ratio of groundwater salinity to direct runoff salinity: C_g/C_r. The salinity accession algorithms are applied to two contrasting field settings: Scott Creek, South Australia and Sandy Creek, northern Queensland, Australia. It was found that k_g > k during periods of obvious stream flow recession, for the events analysed. Salinograph accession behaviour was fairly similar for both sites, despite contrasting environments. Using assumed end‐member salinities for groundwater and direct runoff based upon field observations, the behaviour of k_c from the Scott Creek site was approximately reproduced by varying the initial groundwater to runoff flow ratio: Q_g0/Q_r0, within reasonable parameter ranges. The use of salinograph information when used in addition to standard hydrograph analyses provided useful information on recession characteristics of stream components. Copyright © 2008 John Wiley & Sons, Ltd.     

Mining lakes as groundwater‐dominated hydrological systems: assessment of the water balance of Mining Lake Plessa 117 (Lusatia,Germany) using stable isotopes

In the present study, the stable isotopes δ¹⁸O and δ²H were used for assessment of the water balance in a heterogeneously structured catchment area in the Lusatian Lignite Mining District, in particular, for estimation of the annual groundwater inflow and outflow (I_GW and O_GW) of Mining Lake Plessa 117. The application of stable isotopes was possible since the water exchange in the catchment area had reached steady‐state conditions after the abandonment of mining activities in 1968 and the filling of the voids and aquifers by re‐rising groundwater in the years thereafter. Diverging slopes of the Evaporation Line and the Global Meteoric Water Line manifested as evaporation from the lake catchment area. The calculated isotope water balance was compared with the commonly used surface water balance, which is unable to differentiate between I_GW and O_GW, and with a local groundwater model. The groundwater model calculated an I_GW of about 811 000 m³ yr^?1 and an O_GW close to zero, whereas the isotope water balance showed fluxes of about 914 000 and 140 000 m³ yr^?1, respectively. Considering the contribution of the groundwater inflow to the total annual input into the lake (ΔI_T) and the mean residence time (τ), where the groundwater model and the isotope water balance calculated 42 and 47% for ΔI_T and 4·3 and 3·9 years for τ, respectively, it was shown that both water balance calculation methods led to comparable results despite the differences in I_GW and O_GW. Copyright © 2008 John Wiley & Sons, Ltd.     

Solute dynamics across the stream‐to‐riparian continuum under different flood waves

To enhance the understanding of solute dynamics within the stream‐to‐riparian continuum during flood event‐driven water fluctuation (i.e., flood wave), a variable saturated groundwater flow and solute transport model were developed and calibrated against in situ measurements of the Inbuk stream, Korea, where seasonal flooding prevails. The solute dynamics were further investigated for flood waves (varying by amplitude [A], duration [T], roundness [r], and skewness [t_p]) that were parameterised by real‐time stream stage fluctuations. We found that the solute transferred faster and farther in the riparian zone, especially within the phreatic zone, above which in the variable saturated zone the concentration required a significantly longer time, particularly at higher altitudes, to return to the initial state. By comparison, solute transferred shallowly in the streambed where the solute plume exhibited an exponential growth trend from the centre to the bank. The dynamic changes of solute flux and mass along the stream–aquifer interface and stream concentration were linked to the shape of flood wave. As the flood wave became higher (A↗), wider (T↗), rounder (r↘), and less skewed (t_p↗), the maximum solute storage in aquifer increased. Maximum stream concentration (C_str?max) not only presented a positive linear relationship with A or t_p but also showed a negative logarithmic trend with increasing T or r. The sensitivity of C_{str_max} to A was approximately two times that of t_p, and between these values, the r was slightly more sensitive than T. C_str?max linearly increased as hydraulic conductivity increased and logarithmically increased as longitudinal dispersivity increased. The former relationship was more sensitive than the latter.     

Drip irrigation enhances shallow groundwater contribution to crop water consumption in an arid area

Shallow groundwater plays a key role in agro‐hydrological processes of arid areas. Groundwater often supplies a necessary part of the water requirement of crops and surrounding native vegetation, such as groundwater‐dependent ecosystems. However, the impact of water‐saving irrigation on cropland water balance, such as the contribution of shallow groundwater to field evapotranspiration, requires further investigation. Increased understanding of quantitative evaluation of field‐scale water productivity under different irrigation methods aids policy and decision‐making. In this study, high‐resolution water table depth and soil water content in field maize were monitored under conditions of flood irrigation (FI) and drip irrigation (DI), respectively. Groundwater evapotranspiration (ET_g) was estimated by the combination of the water table fluctuation method and an empirical groundwater–soil–atmosphere continuum model. The results indicate that daily ET_g at different growth stages varies under the two irrigation methods. Between two consecutive irrigation events of the FI site, daily ET_g rate increases from zero to greater than that of the DI site. Maize under DI steadily consumes more groundwater than FI, accounting for 16.4% and 14.5% of ET_a, respectively. Overall, FI recharges groundwater, whereas DI extracts water from shallow groundwater. The yield under DI increases compared with that under FI, with less ET_a (526 mm) compared with FI (578 mm), and irrigation water productivity improves from 3.51 kg m^?3 (FI) to 4.58 kg m^?3 (DI) through reducing deep drainage and soil evaporation by DI. These results highlight the critical role of irrigation method and groundwater on crop water consumption and productivity. This study provides important information to aid the development of agricultural irrigation schemes in arid areas with shallow groundwater.     

Coupled surface water–groundwater model and its application in the arid Shiyang River basin,China

The planning and management of water resources in the Shiyang River basin, China require a tool for assessing the impact of groundwater and stream use on water supply reliabilities and improving many environment‐related problems such as soil desertification induced by recent water‐related human activities. A coupled model, integrating rule‐based lumped surface water model and distributed three‐dimensional groundwater flow model, has been established to investigate surface water and groundwater management scenarios that may be designed to restore the deteriorated ecological environment of the downstream portion of the Shiyang River basin. More than 66% of the water level among 24 observation wells have simulation error less than 1·0 m. The overall trend of the temporal changes of simulated and observed surface runoff at the Caiqi gauging station remains almost the same. The calibration was considered satisfactory. Initial frameworks for water allocation, including agricultural water‐saving projects, water diversion within the basin and inter‐basin water transfer, reducing agricultural irrigation area and surface water use instead of groundwater exploitation at the downstream were figured out that would provide a rational use of water resources throughout the whole basin. Sixteen scenarios were modelled to find out the most appropriate management strategies. The results showed that in the two selected management options, the groundwater budget at the Minqin basin was about 1·4 × 10⁸ m³/a and the ecological environment would be improved significantly, but the deficit existed at the Wuwei basin and the number was about 0·8 × 10⁸ m³/a. Water demand for domestic, industry and urban green area would be met in the next 30 years, but the water shortage for meeting the demand of agricultural water use in the Shiyang River basin was about 2·2 × 10⁸ m³/a. It is suggested that more inter‐basin water transfer should be required to obtain sustainable water resource use in the Shiyang River basin. Copyright © 2009 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.     

Dynamic river–groundwater exchange in the presence of a saline,semi‐confined aquifer

Understanding groundwater–surface water exchange in river banks is crucial for effective water management and a range of scientific disciplines. While there has been much research on bank storage, many studies assume idealized aquifer systems. This paper presents a field‐based study of the Tambo Catchment (southeast Australia) where the Tambo River interacts with both an unconfined aquifer containing relatively young and fresh groundwater (<500 μS/cm and <100 years old) and a semi‐confined artesian aquifer containing old and saline groundwater (electrical conductivity > 2500 μS/cm and >10 000 years old). Continuous groundwater elevation and electrical conductivity monitoring within the different aquifers and the river suggest that the degree of mixing between the two aquifers and the river varies significantly in response to changing hydrological conditions. Numerical modelling using MODFLOW and the solute transport package MT3DMS indicates that saline water in the river bank moves away from the river during flooding as hydraulic gradients reverse. This water then returns during flood recession as baseflow hydraulic gradients are re‐established. Modelling also indicates that the concentration of a simulated conservative groundwater solute can increase for up to ~34 days at distances of 20 and 40 m from the river in response to flood events approximately 10 m in height. For the same flood event, simulated solute concentrations within 10 m of the river increase for only ~15 days as the infiltrating low‐salinity river water drives groundwater dilution. Average groundwater fluxes to the river stretch estimated using Darcy's law were 7 m³/m/day compared with 26 and 3 m³/m/day for the same periods via mass balance using Radon (²²²Rn) and chloride (Cl), respectively. The study shows that by coupling numerical modelling with continuous groundwater–surface water monitoring, the transient nature of bank storage can be evaluated, leading to a better understanding of the hydrological system and better interpretation of hydrochemical data. Copyright © 2015 John Wiley & Sons, Ltd.     

Groundwater contamination risks from conservative point source pollutants in a future climate

ABSTRACT

The groundwater contamination risk in future climates was investigated at three locations in Sweden. Solute transport penetration depths were simulated using the HYDRUS-1D model using historical data and an ensemble of climate projections including two global climate models (GCMs), three emission scenarios and one regional climate model. Most projections indicated increasing precipitation and evapotranspiration until mid-century with a further increase at end-century. Results showed both increasing and decreasing groundwater contamination risks depending on emission scenario and GCM. Generally, the groundwater contamination risk is likely to be unchanged until mid-century, but higher at the end of the century. Soil and site specific relationships between Δ(P – PET) (i.e. change in the difference between precipitation, P, and potential evapotranspiration, PET) and changes in solute transport depths were determined. Using this, changes in solute transport depths for other climate projections can be assessed.     

Coupled groundwater flow and heat transport simulation for estimating transient aquifer–stream exchange at the lowland River Spree (Germany)

Subsurface flow and heat transport near Freienbrink, NE Germany, was simulated in order to study groundwater–surface water exchange between a floodplains aquifer and a section of the lowland River Spree and an adjacent oxbow. Groundwater exfiltration was the dominant process, and only fast surface water level rises resulted in temporary infiltration into the aquifer. The main groundwater flow paths are identified based on a 3D groundwater flow model. To estimate mass fluxes across the aquifer–surface water interfaces, a 2D flow and heat transport modelling approach along a transect of 12 piezometers was performed. Results of steady‐state and transient water level simulations show an overall high accuracy with a Spearman coefficient ρ = 0.9996 and root mean square error (RMSE) = 0.008 m. Based on small groundwater flow velocities of about 10^?7 to 10^?6 ms^?1, mean groundwater exfiltration rates of 233 l m^?2 d^?1 are calculated. Short periods of surface water infiltration into the aquifer do not exceed 10 days, and the infiltration rates are in the same range. The heat transport was modelled with slightly less accuracy (ρ = 0.8359 and RMSE = 0.34 °C). In contrast to the predominant groundwater exfiltration, surface water temperatures determine the calculated temperatures in the upper aquifer below both surface water bodies down to 10 m during the whole simulation period. These findings emphasize prevailing of heat conduction over advection in the upper aquifer zones, which seems to be typical for lowland streams with sandy aquifer materials and low hydraulic gradients. Moreover, this study shows the potential of coupled numerical flow and heat transport modelling to understand groundwater–surface water exchange processes in detail. Copyright © 2013 John Wiley & Sons, Ltd.     

Performance of the White method for estimating groundwater evapotranspiration under conditions of deep and fluctuating groundwater

The White method is a simple but the most frequently applied approach to estimate groundwater evapotranspiration (ET_g) from groundwater level diurnal signals. Because of a lack of direct measurements of ET_g, it is difficult to evaluate the performance of the White method, particularly in field environments with variable groundwater fluctuations. A 2‐year field observation in a hyper‐arid riparian tamarisk (Tamarix spp.) stand with deep groundwater depth in the lower Tarim River basin of China was conducted to measure the surface evapotranspiration (ET_s) and groundwater table. The performance of the White method and the influences of the variable groundwater fluctuations on the determinations of the specific yield (S_y) and recharge rate of groundwater (r) in the White method were investigated. The results showed that the readily available S_y determined by Loheide's method was feasible but must be finely determined based on the soil textures in the layers in which the groundwater level fluctuated. A constant S_y value for a defined porous medium could be assumed regardless of the discharge or recharge processes of groundwater. The time span of 0000 h to 0600 h for r estimation for the White method worked best than other time spans. A 2‐day moving average of r values further enhance ET_g estimation. Slight effects of environmental or anthropogenic disturbances on the diurnal fluctuations of groundwater level did not influence the ET_g estimations by the White method. Our results provide valuable references to the application of the White method for estimating daily ET_g in desert riparian forests with deep groundwater depth. Copyright © 2015 John Wiley & Sons, Ltd.     

The effect of transient storage on nitrate uptake lengths in streams: an inter‐site comparison

Surface water–groundwater interaction in the hyporheic zone may enhance biogeochemical cycling in streams, and it has been hypothesized that streams exchanging more water with the hyporheic zone should have more rapid nitrate utilization. We used simultaneous conservative solute and nitrate addition tracer tests to measure transient storage (which includes hyporheic exchange and in‐stream storage) and the rate of nitrate uptake along three reaches within the Red Canyon Creek watershed, Wyoming. We calibrated a one‐dimensional transport model, incorporating transient storage (OTIS‐P), to the conservative solute breakthrough curves and used the results to determine the degree of transient storage in each reach. The nitrate uptake length was quantified from the exponential decrease in nitrate concentration with distance during the tracer tests. Nitrate uptake along the most downstream reach of Red Canyon Creek was rapid (turnover time K^?1_c = 32 min), compared with nitrate uptake reported in other studies (K^?1_c = 12 to 551 min), but other sites within the watershed showed little nitrate retention or loss. The uptake length S_w‐NO^?₃ for the most downstream reach was 500 m and the mass transfer coefficient V_f‐NO^?₃ was 6·3 m min^?1. Results from 15 other nitrate‐addition tracer tests were used to create a regression model relating transient storage and measures of stream flow to nitrate uptake length. The model, which includes specific discharge and transient storage area, explains almost half the variability in nitrate uptake length (adjusted R² = 0·44) and is most effective for comparing sites with very different stream characteristics. Although large differences in specific discharge and storage zone area explain inter‐site differences in nitrate uptake, other unmeasured variables, such as available organic carbon and microbial community composition, are likely important for predicting differences in nitrate uptake between sites with similar specific discharge rates and storage zone areas, such as when making intra‐site comparisons. Copyright © 2007 John Wiley & Sons, Ltd.     

A framework for sourcing of evaporation between saturated and unsaturated zone in bare soil condition

ABSTRACT

Sourcing subsurface evaporation (E_ss) into groundwater (E_g) and unsaturated zone (E_u) components has received little scientific attention so far, despite its importance in water management and agriculture. We propose a novel sourcing framework, with its implementation in dedicated post-processing software called SOURCE (used along with the HYDRUS1D model), to study evaporation sourcing dynamics, define quantitatively “shallow” and “deep” water table conditions and test the applicability of water table fluctuation (WTF) and “bucket” methods for estimation of E_g and E_u separately.

For the “shallow” and “deep” water table we propose E_g?>?0.95E_ss and E_g = 0 criteria, respectively. Assessment of the WTF method allowed sourcing of very small fluxes otherwise neglected by standard hydrological methods. Sourcing with SOURCE software was more accurate than the standard “bucket” method mainly because of greater flexibility in spatio-temporal discretization. This study emphasized the dry condition relevance of groundwater evaporation which should be analysed by applying coupled flow of heat, vapour and liquid water.

Editor D. Koutsoyiannis; Associate editor S. Kanae     

Hydrochemical characteristics and water quality of the Mušnica River catchment,Bosnia and Herzegovina

Abstract

The spatial position of the Mu?nica River catchment in the Dinaric Karst results in specific conditions of the drainage, distribution and chemical characteristics of the surface water and groundwater. This catchment is mountainous, located between 930 and 1985 m above sea level. The chemical characteristics of the waters are strongly dependent on the local lithological composition and the precipitation regime. The dominant cations are Ca²⁺ and Mg²⁺, and anions HCO₃ ^- and SO₄ ^2-. The surface water is characterized by low and medium mineralization, with resultant soft, medium-hard and very-hard total hardness levels. The groundwater is more mineralized than the surface water. Gacko Town and its associated thermal power plant, Gacko 1, which are located in the southern part of the Mu?nica River catchment, have no modern facilities for sewage and wastewater transport or treatment. In order to prevent further water quality deterioration and to enable its sustainable use, it is recommended to ensure adequate protection measures in the catchment area.

Editor Z.W. Kundzewicz; Associate editor S. Faye

Citation Banjak, D. and Nikoli?, J., 2012. Hydrochemical characteristics and water quality of the Mu?nica River catchment, Bosnia and Herzegovina. Hydrological Sciences Journal, 57 (3), 562–575.     

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.     

Quantitative guidance for efficient vertical flow measurements at the sediment–water interface using temperature–depth profiles

Upward discharge to surface water bodies can be quantified using analytical models based on temperature–depth (T-z) profiles. The use of sediment T-z profiles is attractive as discharge estimates can be obtained using point-in-time data that are collected inexpensively and rapidly. Previous studies have identified that T-z methods can only be applied at times of the year when there is significant difference between the streambed–water interface and deeper sediment temperatures (e.g., winter and summer). However, surface water temperatures also vary diurnally, and the influence of these variations on discharge estimates from T-z methods is poorly understood. For this study, synthetic T-z profiles were generated numerically using measured streambed interface temperature data to assess the influence of diurnal temperature variations on discharge estimation and provide insight into the suitable application of T-z methods. Results show that the time of day of data collection can have a substantial influence on vertical flux estimates using T-z methods. For low groundwater discharge fluxes (e.g., 0.1 m d⁻¹), daily transience in streambed temperatures led to relatively large errors in estimated flow magnitude and direction. For higher discharge fluxes (1.5 m d⁻¹), the influence of transient streambed temperatures on discharge estimates was strongly reduced. Discharge estimates from point-in-time T-z profiles were most accurate when the uppermost point in the T-z profile was near the bed interface daily mean (two time periods daily). Where temperature time series data are available, daily averaged T-z profiles can produce accurate discharge estimates across a wide range of discharge rates. Seasonality in shallow groundwater temperature generally had a negligible influence on vertical flow estimates. These findings can be used to plan field campaigns and provide guidance on the optimal application of T-z methods to quantify vertical groundwater discharge to surface water bodies.     

Estimating agricultural deep drainage lag times to groundwater: application to Antelope Valley,California, USA

Estimates of groundwater volumes available in semiarid regions that rely on water balance calculations require the determination of both surface to groundwater lag times and volumes from irrigation or rainfall initiated recharge. Subsurface geologic material hydraulic properties (e.g. hydraulic conductivities, water retention functions) necessary for unsaturated flow modelling are rarely available as are the instrumented field tests that might determine such lag times. Here we develop a simple two‐parameter (specific yield, S_y, and pore‐size distribution index, λ), one‐dimensional unsaturated flow model from simplifications of the Richards equation (using the Brooks‐Corey relationships) to determine lag times from agricultural deep drainage associated with the irrigation of alfalfa hay and various row crops in the Antelope Valley of California, USA. Model‐predicted lag times to depths of 85 m bgs (below ground surface) were similar to that measured in a 2‐year ponded recharge field trial, slightly overestimating that measured by approximately 15% (0.51 vs 0.44 years). Lag time estimates were most sensitive to estimated deep percolation rates and roughly equally sensitive to the model hydraulic parameters. Generally, as subsurface material textures coarsen towards larger S_y and λ values for all S_y >10%, lag times progressively increase; however, at S_y <10%, lag times decrease substantially suggesting that particular combinations of S_y and λ values that may be associated with similarly textured materials can result in the prediction of different lag times for S_y approximately 10%. Overall, lag times of 1–3 years to a depth of 69 m bgs were estimated from deep drainage of agricultural irrigation across a variety of irrigation schedules and subsurface materials. Copyright © 2012 John Wiley & Sons, Ltd.     

First‐order exchange coefficient coupling for simulating surface water–groundwater interactions: parameter sensitivity and consistency with a physics‐based approach

Distributed hydrologic models capable of simulating fully‐coupled surface water and groundwater flow are increasingly used to examine problems in the hydrologic sciences. Several techniques are currently available to couple the surface and subsurface; the two most frequently employed approaches are first‐order exchange coefficients (a.k.a., the surface conductance method) and enforced continuity of pressure and flux at the surface‐subsurface boundary condition. The effort reported here examines the parameter sensitivity of simulated hydrologic response for the first‐order exchange coefficients at a well‐characterized field site using the fully coupled Integrated Hydrology Model (InHM). This investigation demonstrates that the first‐order exchange coefficients can be selected such that the simulated hydrologic response is insensitive to the parameter choice, while simulation time is considerably reduced. Alternatively, the ability to choose a first‐order exchange coefficient that intentionally decouples the surface and subsurface facilitates concept‐development simulations to examine real‐world situations where the surface‐subsurface exchange is impaired. While the parameters comprising the first‐order exchange coefficient cannot be directly estimated or measured, the insensitivity of the simulated flow system to these parameters (when chosen appropriately) combined with the ability to mimic actual physical processes suggests that the first‐order exchange coefficient approach can be consistent with a physics‐based framework. Copyright © 2009 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.     

Stable isotopic and geochemical identification of groundwater evolution and recharge sources in the arid Shule River Basin of Northwestern China

Stable isotopic (δD_VSMOW and δ¹⁸O_VSMOW) and geochemical signatures were employed to constrain the geochemical evolution and sources of groundwater recharge in the arid Shule River Basin, Northwestern China, where extensive groundwater extraction occurs for agricultural and domestic supply. Springs in the mountain front of the Qilian Mountains, the Yumen‐Tashi groundwater (YTG), and the Guazhou groundwater (GZG) were Ca‐HCO₃, Ca‐Mg‐HCO₃‐SO₄ and Na‐Mg‐SO₄‐Cl type waters, respectively. Total dissolved solids (TDS) and major ion (Mg²⁺, Na⁺, Ca²⁺, K⁺, SO₄^2?, Cl^? and NO₃^?) concentrations of groundwater gradually increase from the mountain front to the lower reaches of the Guazhou Basin. Geochemical evolution in groundwater was possibly due to a combination of mineral dissolution, mixing processes and evapotranspiration along groundwater flow paths. The isotopic and geochemical variations in melt water, springs, river water, YTG and GZG, together with the end‐member mixing analysis (EMMA) indicate that the springs in the mountain front mainly originate from precipitation, the infiltration of melt water and river in the upper reaches; the lateral groundwater from the mountain front and river water in the middle reaches are probably effective recharge sources for the YTG, while contribution of precipitation to YTG is extremely limited; the GZG is mainly recharged by lateral groundwater flow from the Yumen‐Tashi Basin and irrigation return flow. The general characteristics of groundwater in the Shule River Basin have been initially identified, and the results should facilitate integrated management of groundwater and surface water resources in the study area. Copyright © 2015 John Wiley & Sons, Ltd.