Factors controlling seasonal groundwater and solute flux from snow‐dominated basins

Critical zone influences on hydrologic partitioning, subsurface flow paths and reactions along these flow paths dictate the timing and magnitude of groundwater and solute flux to streams. To isolate first‐order controls on seasonal streamflow generation within highly heterogeneous, snow‐dominated basins of the Colorado River, we employ a multivariate statistical approach of end‐member mixing analysis using a suite of daily chemical and isotopic observations. Mixing models are developed across 11 nested basins (0.4 to 85 km²) spanning a gradient of climatological, physical, and geological characteristics. Hydrograph separation using rain, snow, and groundwater as end‐members indicates that seasonal contributions of groundwater to streams is significant. Mean annual groundwater flux ranges from 12% to 33% whereas maximum groundwater contributions of 17% to 50% occur during baseflow. The direct relationship between snow water equivalent and groundwater flux to streams is scale dependent with a trend toward self‐similarity when basins exceed 5.5 km². We find groundwater recharge increases in basins of high relief and within the upper subalpine where maximum snow accumulation is coincident with reduced conifer cover and lower canopy densities. The mixing model developed for the furthest downstream site did not transfer to upstream basins. The resulting error in predicted stream concentrations points toward weathering reactions as a function of source rock and seasonal shifts in flow path. Additionally, the potential for microbial sulfate reduction in floodplain sediments along a low‐gradient, meandering portion of the river is sufficient to modify hillslope contributions and alter mixing ratios in the analysis. Soil flushing in response to snowmelt is not included as an end‐member but is identified as an important mechanism for release of solutes from these mountainous watersheds. End‐member mixing analysis used in combination with high‐frequency observations reveals important aspects of catchment hydrodynamics across scale.     

Automated calibration applied to watershed‐scale flow simulations

Complexity in simulating the hydrological response in large watersheds over long times has prompted a significant need for procedures for automatic calibration. Such a procedure is implemented in the basin‐scale hydrological model (BSHM), a physically based distributed parameter watershed model. BSHM simulates the most important basin‐scale hydrological processes, such as overland flow, groundwater flow and stream–aquifer interaction in watersheds. Here, the emphasis is on estimating the groundwater parameters with water levels in wells and groundwater baseflows selected as the calibration targets. The best set of parameters is selected from within plausible ranges of parameters by adjusting the values of hydraulic conductivity, storativity, groundwater recharge and stream bed permeability. The baseflow is determined from stream flow hydrographs by using an empirical scheme validated using a chemical approach to hydrograph separation. Field studies determined that the specific conductance for components of the composite hydrograph were sufficiently unique to make the chemical approach feasible. The method was applied to the Big Darby Creek Watershed, Ohio. The parameter set selected for the groundwater system provides a good fit with the estimated baseflow and observed water well data. Copyright © 1999 John Wiley & Sons, Ltd.     

Numerical studies on the influences of the South‐to‐North Water Transfer Project on groundwater level changes in the Beijing Plain,China

Groundwater is an important component of the water supply, and overexploitation has triggered many problems in the Beijing Plain. The South‐to‐North Water Transfer Project has been proposed as a promising solution to alleviate these problems. Evaluation of different scenarios of groundwater management after the implementation of the South‐to‐North Water Transfer Project is necessarily required. In this study, a numerical model of groundwater flow was established using FEFLOW software and was well calibrated by parameter optimization and groundwater withdrawal inversion in the Beijing Plain. Sixteen scenarios that considered groundwater exploitation, artificial recharge, and precipitation were designed to simulate the groundwater dynamics after 11 years of the project. The results showed that the groundwater level in the study area would recover to various degrees due to the reductions of groundwater withdrawal and the increments of infiltration; additionally, it was concluded that groundwater was significantly affected by precipitation. Generally, in the designed scenarios, the groundwater‐level increment in the upper streams of the model area was higher than that in the lower streams. The groundwater level would obviously increase from artificial recharge in the immediate and adjacent areas. In addition, modes of reducing exploitation had no significant influence on the change in groundwater level during the 11‐year study period. The developed model offers a reliable and effective way to improve groundwater management.     

Effects of Linking a Soil‐Water‐Balance Model with a Groundwater‐Flow Model

A previously published regional groundwater‐flow model in north‐central Nebraska was sequentially linked with the recently developed soil‐water‐balance (SWB) model to analyze effects to groundwater‐flow model parameters and calibration results. The linked models provided a more detailed spatial and temporal distribution of simulated recharge based on hydrologic processes, improvement of simulated groundwater‐level changes and base flows at specific sites in agricultural areas, and a physically based assessment of the relative magnitude of recharge for grassland, nonirrigated cropland, and irrigated cropland areas. Root‐mean‐squared (RMS) differences between the simulated and estimated or measured target values for the previously published model and linked models were relatively similar and did not improve for all types of calibration targets. However, without any adjustment to the SWB‐generated recharge, the RMS difference between simulated and estimated base‐flow target values for the groundwater‐flow model was slightly smaller than for the previously published model, possibly indicating that the volume of recharge simulated by the SWB code was closer to actual hydrogeologic conditions than the previously published model provided. Groundwater‐level and base‐flow hydrographs showed that temporal patterns of simulated groundwater levels and base flows were more accurate for the linked models than for the previously published model at several sites, particularly in agricultural areas.     

Transient Recharge Estimability Through Field‐Scale Groundwater Model Calibration

The estimation of recharge through groundwater model calibration is hampered by the nonuniqueness of recharge and aquifer parameter values. It has been shown recently that the estimability of spatially distributed recharge through calibration of steady‐state models for practical situations (i.e., real‐world, field‐scale aquifer settings) is limited by the need for excessive amounts of hydraulic‐parameter and groundwater‐level data. However, the extent to which temporal recharge variability can be informed through transient model calibration, which involves larger water‐level datasets, but requires the additional consideration of storage parameters, is presently unknown for practical situations. In this study, time‐varying recharge estimates, inferred through calibration of a field‐scale highly parameterized groundwater model, are systematically investigated subject to changes in (1) the degree to which hydraulic parameters including hydraulic conductivity (K) and specific yield (S_y) are constrained, (2) the number of water‐level calibration targets, and (3) the temporal resolution (up to monthly time steps) at which recharge is estimated. The analysis involves the use of a synthetic reality (a reference model) based on a groundwater model of Uley South Basin, South Australia. Identifiability statistics are used to evaluate the ability of recharge and hydraulic parameters to be estimated uniquely. Results show that reasonable estimates of monthly recharge (<30% recharge root‐mean‐squared error) require a considerable amount of transient water‐level data, and that the spatial distribution of K is known. Joint estimation of recharge, S_y and K, however, precludes reasonable inference of recharge and hydraulic parameter values. We conclude that the estimation of temporal recharge variability through calibration may be impractical for real‐world settings.     

Semi‐Arid Aquifer Responses to Forest Restoration Treatments and Climate Change

The purpose of this study was to develop an interpretive groundwater‐flow model to assess the impacts that planned forest restoration treatments and anticipated climate change will have on large regional, deep (>400 m), semi‐arid aquifers. Simulations were conducted to examine how tree basal area reductions impact groundwater recharge from historic conditions to 2099. Novel spatial analyses were conducted to determine areas and rates of potential increases in groundwater recharge. Changes in recharge were applied to the model by identifying zones of basal area reduction from planned forest restoration treatments and applying recharge‐change factors to these zones. Over a 10‐year period of forest restoration treatment, a 2.8% increase in recharge to one adjacent groundwater basin (the Verde Valley sub‐basin) was estimated, compared to conditions that existed from 2000 to 2005. However, this increase in recharge was assumed to quickly decline after treatment due to regrowth of vegetation and forest underbrush and their associated increased evapotranspiration. Furthermore, simulated increases in groundwater recharge were masked by decreases in water levels, stream baseflow, and groundwater storage resulting from surface water diversions and groundwater pumping. These results indicate that there is an imbalance between water supply and demand in this regional, semi‐arid aquifer. Current water management practices may not be sustainable into the far future and comprehensive action should be taken to minimize this water budget imbalance.     

Application of an integrated surface water‐groundwater model to multi‐aquifers modeling in Choushui River alluvial fan,Taiwan

This research proposes a combination of SWAT and MODFLOW, MD‐SWAT‐MODFLOW, to address the multi‐aquifers condition in Choushui River alluvial fan, Taiwan. The natural recharge and unidentified pumping/recharge are separately estimated. The model identifies the monthly pumping/recharge rates in multi‐aquifers so that the daily streamflow can be simulated correctly. A multi‐aquifers condition means a subsurface formation composed of at least the unconfined aquifer, the confined aquifer, and an in‐between aquitard. In such a case, the variation of groundwater level is related to pumping/recharge activities in vertically adjacent aquifer and the river‐aquifer interaction. Both factors in turn affect the streamflow performance. Results show that MD‐SWAT‐MODFLOW performs better than SWAT alone in terms of simulated streamflow, especially during low flow period, when pumping/recharge rates are properly estimated. A sensitivity analysis of individual parameter suggests that the vertical leakance may be the most sensitive among all investigated MODFLOW parameters in terms of the estimated pumping/recharge among aquifers, and the Latin‐Hypercube‐One‐factor‐At‐a‐Time sensitivity analysis indicates that the hydraulic conductivity of channel is the most sensitive to the model performance. It also points out the necessity to simultaneously estimate pumping/recharge rates in multi‐aquifers. The estimated net pumping rate can be treated as a lower bound of the actual local pumping rate. As a whole, the model provides the spatio‐temporal groundwater use, which gives the authorities insights to manage groundwater resources. Copyright © 2012 John Wiley & Sons, Ltd.     

Modelling two‐dimensional steady‐state groundwater flow and flow sensitivity to boundary conditions in blanket peat complexes

This study used a two‐dimensional steady‐state finite‐element groundwater flow model to simulate groundwater flow in two Newfoundland blanket peat complexes and to examine flow system sensitivity to changes in water table recharge and aquifer properties. The modelling results were examined within the context of peat‐forming processes in the two complexes. Modelled flow compared favourably with observed flow. The sensitivity analyses suggested that more highly decomposed bog peat along bog margins probably has/had a positive impact on net peat accumulation within bog interiors. Peat with lower hydraulic conductivity along bog margins effectively impedes lateral drainage, localizes water table drawdown to extreme bog margins, and elevates water tables along bog interiors. Peat formation and elevated water tables in adjacent poor fens/laggs currently rely on placic and ortstein horizons impeding vertical drainage and water flow inputs from adjacent bogs. Modest reductions in atmospheric recharge were found to govern bog‐flow‐system geometries in a way that would adversely affect paludification processes in adjacent fens/laggs. Copyright © 2004 John Wiley & Sons, Ltd.     

A geomorphology‐based integrated stream–aquifer interaction model for semi‐gauged catchments

Many of the existing stream–aquifer interaction models available in the literature are very complex with limited applicability in semi‐gauged and ungauged catchments. In this study, to estimate the influent and effluent subsurface water fluxes under limited geo‐hydrometeorological data availability conditions, a simple stream–aquifer interaction model, namely, the variable parameter McCarthy–Muskingum (VPMM) hillslope‐storage Boussinesq (hsB) model, has been developed. This novel model couples the VPMM streamflow transport with the hsB groundwater flow transport modules in online mode. In this integrated model, the surface water–groundwater flux exchange process is modelled by the Darcian approach with the variable hydraulic heads between the river stage and groundwater table accounting for the rainfall forcing. Considering the exchange fluxes in the hyporheic zone and lateral overland flow contribution, this approach is field tested in a typical 48‐km stretch of the Brahmani River in eastern India to simulate the streamflow and its depth with the minimum Nash–Sutcliffe efficiency of 94% and 88%; the maximum root mean square error of 134 m³/s and 0.35 m; and the minimum index of agreement of 98% and 97%, respectively. This modelling approach could be very well utilized in data‐scarce world‐river basins to estimate the stream–aquifer exchange flux due to rainfall forcings.     

A large weighing lysimeter for evapotranspiration and soil‐water–groundwater exchange studies

A large weighing lysimeter was installed at Yucheng Comprehensive Experimental Station, north China, for evapotranspiration and soil‐water–groundwater exchange studies. Features of the lysimeter include the following: (i) mass resolution equivalent to 0·016 mm of water to accurately and simultaneously determine hourly evapotranspiration, surface evaporation and groundwater recharge; (ii) a surface area of 3·14 m² and a soil profile depth of 5·0 m to permit normal plant development, soil‐water extraction, soil‐water–groundwater exchanges, and fluctuations of groundwater level; (iii) a special supply–drainage system to simulate field conditions of groundwater within the lysimeter; (iv) a soil mass of about 30 Mg, including both unsaturated and saturated loam. The soil consists mainly of mealy sand and light loam. Monitoring the vegetated lysimeter during the growing period of winter wheat, from October 1998 through to June 1999, indicated that during the period groundwater evaporation contributed 16·6% of total evapotranspiration for a water‐table depth from 1·6 m to 2·4 m below ground surface. Too much irrigation reduced the amount of upward water flow from the groundwater table, and caused deep percolation to the groundwater. Data from neutron probe and tensiometers suggest that soil‐water‐content profiles and soil‐water‐potential profiles were strongly affected by shallow groundwater. Copyright © 2000 John Wiley & Sons, Ltd.     

Thermal effects of lateral supra‐permafrost water flow around a thermokarst lake on the Qinghai–Tibet Plateau

Both the inflow and outflow of supra‐permafrost water to lakes play important roles in the hydrologic process of thermokarst lakes. The accompanying thermal effects on the adjacent permafrost are required for assessing their influences on the development of thermokarst lakes. For these purposes, the lake water level, temperature dynamics, and supra‐permafrost water flow of a lake were monitored on the Qinghai‐Tibet Plateau. In addition, the spatial and temporal variation of the active layer thickness and permafrost distribution around the lake were investigated by combining ground penetrating radar, electrical resistivity tomography, and borehole temperature monitoring. The results revealed that the yearly unfrozen supra‐permafrost water flow around the lake lasted approximately 5 months. The temperature and water level measurements during this period indicate that the lake water was recharged by relatively colder supra‐permafrost water from the north‐western lakeshore and was discharged through the eastern lakeshore. This process, accompanied by heat exchange with the underlying permafrost, might cause a directional difference of the active layer thickness and permafrost characteristics around the lake. Specifically, the active layer thickness variation was minimal, and the ice‐rich permafrost was found adjacent to the lakeshore along the recharge groundwater pathways, whereas a deeper active layer and ice‐poor permafrost were observed close to the lakeshore from which the warm lake water was discharged. This study suggests that the lateral flow of warm lake water can be a major driver for the rapid expansion of thermokarst lakes and provides clues for evaluating the relationships between the thermokarst expansion process and climate warming.     

Isotopic evidence for widespread cold‐season‐biased groundwater recharge and young streamflow across central Canada

Transformations of precipitation into groundwater and streamflow are fundamental hydrological processes, critical to irrigated agriculture, hydroelectric power generation, and ecosystem health. Our understanding of the timing of groundwater recharge and streamflow generation remains incomplete, limiting our ability to predict fresh water, nutrient, and contaminant fluxes, especially in large basins. Here, we analyze thousands of rain, snow, groundwater, and streamflow δ¹⁸O and δ²H values in the Nelson River basin, which covers 1.2 million km² of central Canada. We show that the fraction of precipitation that recharges aquifers is ~1.3–5 times higher for precipitation falling during cold months with subzero mean monthly temperatures than for precipitation falling during warmer months. The near‐ubiquity of cold‐season‐biased groundwater recharge implies that changes to winter water balances may have disproportionate impacts on annual groundwater recharge rates. We also show that young streamflow—defined as precipitation that enters a river in less than ~2.3 months—comprises ~27% of annual streamflow but varies widely among tributaries in the Nelson River basin (1–59%). Young streamflow fractions are lower in steep catchments and higher in flatter catchments such as the transboundary Red River basin. Our findings imply that flat, lower permeability, heavily tiled landscapes favor more rapid transmission of precipitation into rivers, possibly mobilizing excess soluble fertilizers and exacerbating eutrophication events in Lake Winnipeg.     

Assessing regional‐scale spatio‐temporal patterns of groundwater–surface water interactions using a coupled SWAT‐MODFLOW model

Interaction between groundwater and surface water in watersheds has significant impacts on water management and water rights, nutrient loading from aquifers to streams, and in‐stream flow requirements for aquatic species. Of particular importance are the spatial patterns of these interactions. This study explores the spatio‐temporal patterns of groundwater discharge to a river system in a semi‐arid region, with methods applied to the Sprague River Watershed (4100 km²) within the Upper Klamath Basin in Oregon, USA. Patterns of groundwater–surface water interaction are explored throughout the watershed during the 1970–2003 time period using a coupled SWAT‐MODFLOW model tested against streamflow, groundwater level and field‐estimated reach‐specific groundwater discharge rates. Daily time steps and coupling are used, with groundwater discharge rates calculated for each model computational point along the stream. Model results also are averaged by month and by year to determine seasonal and decadal trends in groundwater discharge rates. Results show high spatial variability in groundwater discharge, with several locations showing no groundwater/surface water interaction. Average annual groundwater discharge is 20.5 m³/s, with maximum and minimum rates occurring in September–October and March–April, respectively. Annual average rates increase by approximately 0.02 m³/s per year over the 34‐year period, negligible compared with the average annual rate, although 70% of the stream network experiences an increase in groundwater discharge rate between 1970 and 2003. Results can assist with water management, identifying potential locations of heavy nutrient mass loading from the aquifer to streams and ecological assessment and planning focused on locations of high groundwater discharge. Copyright © 2016 John Wiley & Sons, Ltd.     

Impact of complex aquifer geometry on groundwater storage in high‐elevation meadows of the Sierra Nevada Mountains,CA

Recent research has indicated that Sierra Nevada meadows are hydrologically more complex than previously considered. Improved understanding of the effects of aquifer parameters and climate change on water resources in and downstream of meadows is critically needed to effectively manage mountain meadows for ecosystem services and watershed contributions. This research investigates the roles of bedrock geometry, saturated hydraulic conductivity, and meadow gradient in affecting groundwater storage dynamics and surface‐water outflows in site‐scale high‐elevation meadows. Under current and projected lower snowpack conditions, we modeled groundwater flow in representative high‐elevation meadows considering 2 conceptual aquifer thickness models: uniform and variable thickness. Spatially, variable aquifer thicknesses interpreted from bedrock depths (0–28 m) were identified from a high‐resolution ground‐penetrating radar survey conducted at Tuolumne Meadows, CA. Our interpreted bedrock surface indicated several buried U‐shaped valleys including a buried ridge that separates 2 U‐shaped valleys. Groundwater flow simulations show that an increase in meadow gradient and hydraulic conductivity led to a decrease in seasonal storage and an increase in surface‐water outflow. However, models with varying bedrock geometries change the magnitude and timing of these processes. Uniform thickness models overestimated storage at the model edges and resulted in higher projected volumes of water being released to streams earlier than previously observed.     

Evaluating the use of a gridded climate surface for modelling groundwater recharge in a semi‐arid region (Okanagan Basin,Canada)

Spatially distributed groundwater recharge was simulated for a segment of a semi‐arid valley using three different treatments of meteorological input data and potential evapotranspiration (PET). For the same area, timeframe, land cover characteristics and soil properties, groundwater recharge was estimate using (i) single‐station climate data with monthly PET calculated by the Thornthwaite method; (ii) single‐station climate data with daily PET calculated by the Penman–Monteith method; and (iii) daily gridded climate data with spatially distributed PET calculated using the Penman–Monteith method. For each treatment, the magnitude and distribution of actual evapotranspiration (AET) for summer months compared well with those estimated for a 5‐year crop study, suggesting that the near‐surface hydrological processes were replicated and that subsequent groundwater recharge rates are realistic. However, for winter months, calculated AET was near zero when using the Thornthwaite PET method. Mean annual groundwater recharge varied from ～3·2 to 10·0 mm when PET was calculated by the Thornthwaite method, and from ～1·8 to 7·5 mm when PET was calculated by the Penman–Monteith method. Comparisons of bivariate plots of seasonal recharge rates estimated from single‐station versus gridded surface climate reveal that there is greater variability between the different methods for spring months, which is the season of greatest recharge. Furthermore, these seasonal differences are shown to provide different results when compared to the depth to water table, which could lead to different results of evaporative extinction depth. These findings illustrate potential consequences of using different approaches for representing spatial meteorological input data, which could provide conflicting predictions when modelling the influence of climate change on groundwater recharge. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

Water use strategies of two co‐occurring tree species in a semi‐arid karst environment

The flow of precipitation from the surface through to groundwater in karst systems is a complex process involving storage in the unsaturated zone and diffuse and preferential recharge pathways. The processes associated with this behaviour are not well understood, despite the prevalence of karst aquifers being used as freshwater supplies. As a result, uncertainty regarding the ecohydrological processes in this geological setting remains large. In response to the need to better understand the impact of woody vegetation on groundwater recharge, annual evapotranspiration (ET) rates and tree water sources were measured for two years above a shallow, fresh karst aquifer. Water use strategies of the co‐occurring Eucalyptus diversifolia subsp. diversifolia Bonpl. and Allocasuarina verticillata (Lam.) L. Johnson were investigated using a monthly water balance approach, in conjunction with measurement of the stable isotopes of water, leaf water potentials and soil matric potentials. The results suggest that it is unlikely groundwater resources are required to sustain tree transpiration, despite its shallow proximity to the soil surface, and that similarities exist between ET losses and the estimated long‐term average rainfall for this area. Irrespective of stand and morphological differences, E. diversifolia and A. verticillata ET rates showed remarkable convergence, demonstrating the ability of these co‐occurring species to maximise their use of the available precipitation, which avoids the requirement to differentiate between these species when estimating ET at a landscape scale. We conclude that the water holding capacity of porous geological substrates, such as those associated with karst systems, will play an important role in equilibrating annual rainfall variability and should be considered when assessing ecohydrological links associated with karst systems. Copyright © 2013 John Wiley & Sons, Ltd.     

Divergence and flow direction as indicators of subsurface heterogeneity and stage‐dependent storage in alluvial floodplains

Assuming homogeneity in alluvial aquifers is convenient, but limits our ability to accurately predict stream‐aquifer interactions. Research is needed on (i) identifying the presence of focused, as opposed to diffuse, groundwater discharge/recharge to streams and (ii) the magnitude and role of large‐scale bank and transient storage in alluvial floodplains relative to changes in stream stage. The objective of this research was to document and quantify the effect of stage‐dependent aquifer heterogeneity and bank storage relative to changes in stream stage using groundwater flow divergence and direction. Monitoring was performed in alluvial floodplains adjacent to the Barren Fork Creek and Honey Creek in northeastern Oklahoma. Based on results from subsurface electrical resistivity mapping, observation wells were installed in high and low electrical resistivity subsoils. Water levels in the wells were recorded real time using pressure transducers (August to October 2009). Divergence was used to quantify heterogeneity (i.e. variation in hydraulic conductivity, porosity, and/or aquifer thickness), and flow direction was used to assess the potential for large‐scale (100 m) bank or transient storage. Areas of localized heterogeneity appeared to act as divergence zones allowing stream water to quickly enter the groundwater system, or as flow convergence zones draining a large groundwater area. Maximum divergence or convergence occurred with maximum rates of change in flow rates or stream stage. Flow directions in the groundwater changed considerably between base and high flows, suggesting that the floodplains acted as large‐scale bank storage zones, rapidly storing and releasing water during passage of a storm hydrograph. During storm events at both sites, the average groundwater direction changed by at least 90° from the average groundwater direction during baseflow. Aquifer heterogeneity in floodplains yields hyporheic flows that are more responsive and spatially and temporally complex than would be expected compared to more common assumptions of homogeneity. Copyright © 2012 John Wiley & Sons, Ltd.     

Groundwater Recharge Rate and Zone Structure Estimation Using PSOLVER Algorithm

The quantification of groundwater recharge is an important but challenging task in groundwater flow modeling because recharge varies spatially and temporally. The goal of this study is to present an innovative methodology to estimate groundwater recharge rates and zone structures for regional groundwater flow models. Here, the unknown recharge field is partitioned into a number of zones using Voronoi Tessellation (VT). The identified zone structure with the recharge rates is associated through a simulation‐optimization model that couples MODFLOW‐2000 and the hybrid PSOLVER optimization algorithm. Applicability of this procedure is tested on a previously developed groundwater flow model of the Tahtal? Watershed. Successive zone structure solutions are obtained in an additive manner and penalty functions are used in the procedure to obtain realistic and plausible solutions. One of these functions constrains the optimization by forcing the sum of recharge rates for the grid cells that coincide with the Tahtal? Watershed area to be equal to the areal recharge rate determined in the previous modeling by a separate precipitation‐runoff model. As a result, a six‐zone structure is selected as the best zone structure that represents the areal recharge distribution. Comparison to results of a previous model for the same study area reveals that the proposed procedure significantly improves model performance with respect to calibration statistics. The proposed identification procedure can be thought of as an effective way to determine the recharge zone structure for groundwater flow models, in particular for situations where tangible information about groundwater recharge distribution does not exist.     

A comparison of groundwater recharge estimation methods in a semi‐arid,coastal avocado and citrus orchard (Ventura County,California)

We assess the relative merits of application of the most commonly used field methods (soil‐water balance (SWB), chloride mass balance (CMB) and soil moisture monitoring (NP)) to determine recharge rates in micro‐irrigated and non‐irrigated areas of a semi‐arid coastal orchard located in a relatively complex geological environment. Application of the CMB method to estimate recharge rates was difficult owing to the unusually high, variable soil‐water chloride concentrations. In addition, contrary to that expected, the chloride concentration distribution at depths below the root zone in the non‐irrigated soil profiles was greater than that in the irrigated profiles. The CMB method severely underestimated recharge rates in the non‐irrigated areas when compared with the other methods, although the CMB method estimated recharge rates for the irrigated areas, that were similar to those from the other methods, ranging from 42 to 141 mm/year. The SWB method, constructed for a 15‐year period, provided insight into the recharge process being driven by winter rains rather than summer irrigation and indicated an average rate of 75 mm/year and 164 mm/year for the 1984 – 98 and 1996 – 98 periods, respectively. Assuming similar soil‐water holding capacity, these recharge rates applied to both irrigated and non‐irrigated areas. Use of the long period of record was important because it encompassed both drought and heavy rainfall years. Successful application of the SWB method, however, required considerable additional field measurements of orchard ET_c, soil‐water holding capacity and estimation of rainfall interception – runoff losses. Continuous soil moisture monitoring (NP) was necessary to identify both daily and seasonal seepage processes to corroborate the other recharge estimates. Measured recharge rates during the 1996 – 1998 period in both the orchards and non‐irrigated site averaged 180 mm/year. The pattern of soil profile drying during the summer irrigation season, followed by progressive wetting during the winter rainy season was observed in both irrigated and non‐irrigated soil profiles, confirming that groundwater recharge was rainfall driven and that micro‐irrigation did not ‘predispose’ the soil profile to excess rainfall recharge. The ability to make this recharge assessment, however, depended on making multiple field measurements associated with all three methods, suggesting that any one should not be used alone. Copyright © 2000 John Wiley & Sons, Ltd.