Some relationships between lithology,basin form and hydrology: a case study from the Thames basin,UK

The role of lithology in influencing basin form and function is explored empirically by investigating correlations between a range of catchment variables, where the spatial unit of analysis is not surface catchments but lithologically coherent groundwater units. Using the Thames basin, UK, as a case study, nine groundwater units have been identified. Values for 11 hydrological and geomorphological variables, including rainfall, drainage density, Baseflow Index, aquifer porosity, storage coefficient and log‐hydraulic conductivity, aquifer and drainage elevation, river incision, and hypsometric integral have been estimated for each of the groundwater units in the basin, and Pearson correlation coefficients calculated for all pairs of variables. Seven of the correlation coefficients are found to be significant at a confidence level of > 99%. Negative correlations between drainage density and log aquifer hydraulic conductivity, and between drainage density and river incision, and positive correlations between log‐hydraulic conductivity and river incision, log‐hydraulic conductivity and Baseflow Index, and between Baseflow Index and river incision are inferred to have consistent causal explanations. For example, incision of rivers into aquifers leads to relative increases in hydraulic gradients in the vicinity of rivers which, in turn, promotes the development of secondary porosity increasing both aquifer hydraulic conductivity and, hence, Baseflow Index. The implication of this interpretation is that the geomorphological evolution of basins is intimately linked to the evolution of hydraulic conductivity of the underlying aquifers. This is consistent with, and supports the notion of a coupled complexly evolving surface water‐groundwater system. Copyright © 2011 John Wiley & Sons, Ltd.     

Simulation of groundwater flow in a crystalline rock aquifer system in Southern Ghana – An evaluation of the effects of increased groundwater abstraction on the aquifers using a transient groundwater flow model

Monitored groundwater level data, well logs, and aquifer data as well as the relevant surface hydrological data were used to conceptualise the hydrogeological system of the Densu Basin in Southern Ghana. The objective was to numerically derive the hydraulic conductivity field for better characterization of the aquifer system and for simulating the effects of increasing groundwater abstraction on the aquifer system in the basin. The hydraulic conductivity field has been generated in this study through model calibration. This study finds that hydraulic conductivity ranges between a low of 2 m/d in the middle sections of the basin and about 40 m/d in the south. Clear differences in the underlying geology have been indicated in the distribution of aquifer hydraulic conductivities. This is in consonance with the general assertion that the hydrogeological properties of the aquifers in the crystalline basement terrains are controlled by the degree of fracturing and/or weathering of the country rock. The transient model suggest aquifer specific storage values to range between 6.0 × 10^?5 m^?1 and 2.1 × 10^?4 m^?1 which are within acceptable range of values normally quoted for similar lithologies in the literature. There is an apparent subtle decrease in groundwater recharge from about 13% of the annual precipitation in 2005 to about 10.3% of the precipitation in 2008. The transient model was used to simulate responses of the system to annual increment of groundwater abstraction by 20% at the 2008 recharge rates for the period 2009 – 2024. The results suggest that the system will not be able to sustain this level of abstraction as it would lead to a basin wide drawdown in the hydraulic head by 4 m by the end of the prediction period. It further suggests a safe annual increment in groundwater abstraction by 5% under business as usual recharge conditions. Identification and protection of groundwater recharge areas in the basin are recommended in order to safeguard the integrity of the resource under the scenario of increased abstraction for commercial activities in the basin. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

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.     

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

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

Distinguishing Mountain Front and Mountain Block Recharge in an Intermontane Basin of the Himalayan Region

Intermontane basin aquifers worldwide, particularly in the Himalayan region, are recharged largely by the adjoining mountains. Recharge in these basins can occur either by water infiltrating from streams near mountain fronts (MFs) as mountain front recharge (MFR) or by sub-surface mountain block infiltration as mountain block recharge (MBR). MFR and MBR recharge are challenging to distinguish and are least quantified, considering the lack of extensive understanding of the hydrological processes in the mountains. This study used oxygen and hydrogen isotopes (δ¹⁸O and δ²H), electrical conductivity (EC) data, hydraulic head, and groundwater level data to differentiate MFR and MBR. Groundwater level data provide information about the groundwater-surface water interactions and groundwater flow directions, whereas isotopes and EC data are used to distinguish and quantify different recharge sources. The present methodology is tested in an intermontane basin of the Himalayan region. The results suggest that karst springs (KS) and deep groundwater (DGW) recharge are dominated by snowmelt (47% ± 10% and 46% ± 9%) as MBR from adjacent mountains, insignificantly affected by evaporation. The hydraulic head data and isotopes indicate Quaternary shallow groundwater (SGW) aquifer system recharge as MFR of local meteoric water with significant evaporation. The results indicate several flow paths in the aquifer system, a local flow for KS, intermediate flow for SGW, and regional flow for DGW. The findings will significantly impact water resource management in the area and provide vital baseline knowledge for sustainable groundwater management in other Himalayan intermontane basins.     

Assessing hydrogeologic controls on dynamic groundwater storage using long‐term instrumental records of water table levels

This study analyzes a long‐term regional compilation of water table response to climate variability based on 124 long‐term groundwater wells distributed across New England, USA, screened in a variety of geologic materials. The New England region of the USA is located in a humid‐temperature climate underlain by low‐storage‐fractured metamorphic and crystalline bedrock dissected by north–south trending valleys filled with glacial and post‐glacial valley fill sediments. Uplands are covered by thin glacial till that comprises more than 60% of the total area. Annual and multi‐annual responses of the water table to climate variability are assessed to understand how local hydraulic properties and hydrogeologic setting (located in recharge/discharge region) of the aquifer influence the hydrologic sensitivity of the aquifer system to climate variability. This study documents that upland aquifer systems dominated by thin deposits of surface till comprise ~70% of the active and dynamic storage of the region. Total aquifer storage changes of +5 to ?7 km³ occur over the region during the study interval. The storage response is dominated by thin and low permeability surficial till aquifer that fills and drains on a multi‐annual basis and serves as the main mechanism to deliver water to valley fill aquifers and underlying bedrock aquifers. Whereas the till aquifer system is traditionally neglected as an important storage reservoir, this study highlights the importance of a process‐based understanding of how different landscape hydrogeologic units contribute to the overall hydrologic response of a region.     

Hydrogeological interactions between fault zones and alluvial aquifers in regional flow systems

The hydrological influence of fault zones in tectonic areas is usually difficult to depict from field data. Numerical simulation allows representation of such flow systems and an estimation of flow lines and rates. This paper reports on simulations of the groundwater flow in a range‐and‐basin area affected by a regional fault zone, which may drain or recharge an overlaying alluvial aquifer. Different hydraulic conductivity values for the range rocks, the fault‐zone, and the sedimentary infill of the basin are considered, as well as different fault‐zone widths and boundary conditions. Results show that upward and downward fluxes develop in the upper part of the fault zone, controlled by the action of the alluvial aquifer, influencing the recharge of the sedimentary basin. This paper shows the hydrological efficiency of fault zones as preferential flow; it also analyses the constraints that determine groundwater recharge to the surrounding basins. These results contribute to the understanding of hydrogeological dynamics in tectonic areas. Copyright © 2008 John Wiley & Sons, Ltd.     

Dissolution in a variably confined carbonate platform: effects of allogenic runoff,hydraulic damming of groundwater inputs,and surface–groundwater exchange at the basin scale

In variably confined carbonate platforms, impermeable confining units collect rainfall over large areas and deliver runoff to rivers or conduits in unconfined portions of platforms. Runoff can increase river stage or conduit heads in unconfined portions of platforms faster than local infiltration of rainfall can increase groundwater heads, causing hydraulic gradients between rivers, conduits and the aquifer to reverse. Gradient reversals cause flood waters to flow from rivers and conduits into the aquifer where they can dissolve limestone. Previous work on impacts of gradient reversals on dissolution has primarily emphasized individual caves and little research has been conducted at basin scales. To address this gap in knowledge, we used legacy data to assess how a gradient of aquifer confinement across the Suwannee River Basin, north‐central Florida affected locations, magnitudes and processes of dissolution during 2005–2007, a period with extreme ranges of discharge. During intense rain events, runoff from the confining unit increased river stage above groundwater heads in unconfined portions of the platform, hydraulically damming inputs of groundwater along a 200 km reach of river. Hydraulic damming allowed allogenic runoff with SI_CAL < ?4 to fill the entire river channel and flow into the aquifer via reversing springs. Storage of runoff in the aquifer decreased peak river discharges downstream and contributed to dissolution within the aquifer. Temporary storage of allogenic runoff in karst aquifers represents hyporheic exchange at a scale that is larger than found in streams flowing over non‐karst aquifers because conduits in karst aquifers extend the area available for exchange beyond river beds deep into aquifers. Post‐depositional porosity in variably confined carbonate platforms should thus be enhanced along rivers that originate on confining units. This distribution should be considered in models of porosity distribution used to manage water and hydrocarbon resources in carbonate rocks. Copyright © 2013 John Wiley & Sons, Ltd.     

Groundwater flow and implications for groundwater contamination north of Prewitt,New Mexico,U.S.A.

In the southern San Juan Basin, New Mexico, strata of Permian and younger age dip gently toward the center of the basin. Most previous investigators believed that recharge to these strata occurred by precipitation on the outcrops and groundwater flowed downdip to the north and northeast. Recent water-level measurements in an undeveloped part of the basin near Prewitt, New Mexico, show that groundwater at shallow depths in alluvium and bedrock flows southward, opposite to the dip direction, and toward a major ephemeral drainage in a strike valley. North of this area, groundwater in deep bedrock aquifers does appear to flow northward. This information suggests that there are two groundwater circulation patterns; a shallow one controlled by topography and a deeper one controlled by geologic structure.Significant amounts of recharge to sandstone aquifers by infiltration through outcrops is unlikely due to the near-vertical exposures on cliffs, the gentle dip of the strata, and small annual precipitation. Numerical model results suggest that recharge to bedrock aquifers may be from downward leakage via aquitards over large areas and leakage from narrow alluvial aquifers in the subcrop area. The recharge mechanism is controlled by the hydraulic conductivity of the strata.As the flow path is controlled by hydraulic conductivity contrasts, geologic structure, and topography, contamination movement from surface impoundments is likely to be difficult to predict without a thorough hydrogeological site investigation.     

Performance Assessments of a Novel Well Design for Reducing Exposure to Bedrock‐Derived Arsenic

Arsenic in groundwater is a serious problem in New England, particularly for domestic well owners drawing water from bedrock aquifers. The overlying glacial aquifer generally has waters with low arsenic concentrations but is less used because of frequent loss of well water during dry periods and the vulnerability to surface‐sourced bacterial contamination. An alternative, novel design for shallow wells in glacial aquifers is intended to draw water primarily from unconsolidated glacial deposits, while being resistant to drought conditions and surface contamination. Its use could greatly reduce exposure to arsenic through drinking water for domestic use. Hypothetical numerical models were used to investigate the potential hydraulic performance of the new well design in reducing arsenic exposure. The aquifer system was divided into two parts, an upper section representing the glacial sediments and a lower section representing the bedrock. The location of the well, recharge conditions, and hydraulic properties were systematically varied in a series of simulations and the potential for arsenic contamination was quantified by analyzing groundwater flow paths to the well. The greatest risk of arsenic contamination occurred when the hydraulic conductivity of the bedrock aquifer was high, or where there was upward flow from the bedrock aquifer because of the position of the well in the flow system.     

Identification of Groundwater Parameters Using an Adaptative Multiscale Method

The identification of groundwater parameters in heterogeneous systems is a major challenge in groundwater modeling. Flexible parameterization methods are needed to assess the complexity of the spatial distributions of these parameters in real aquifers. In this article, we introduce an adaptative parameterization to identify the distribution of hydraulic conductivity within the large‐scale (4400 km²) Upper Rhine aquifer. The method is based on adaptative multiscale triangulation (AMT) coupled with an inverse problem procedure that identifies the parameters' distributions by reducing the error between measured and simulated heads. The AMT method has the advantage of combining both zonation and interpolation approaches. The AMT method uses area‐based interpolation rather than an interpolation based on stochastic features. The method is applied to a standard 2D groundwater model that takes into account the interactions between the aquifer and surface water bodies, groundwater recharge, and pumping wells. The simulation period covers 204 months, from January 1986 to December 2002. Recordings at 109 piezometers are used for model calibration. The simulated heads are globally quite accurate and reproduce the main dynamics of the system. The local hydraulic conductivities resulting from the AMT method agree qualitatively with existing local experimental observations across the Rhine aquifer.     

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.     

Modification of the SWAT model to simulate regional groundwater flow using a multicell aquifer

The soil and water assessment tool (SWAT) has been widely used and thoroughly tested in many places in the world. The application of the SWAT model has pointed out that 2 of the major weaknesses of SWAT are related to the nonspatial reference of the hydrologic response unit concept and to the simplified groundwater concept, which contribute to its low performance in baseflow simulation and its inability to simulate regional groundwater flow. This study modified the groundwater module of SWAT to overcome the above limitations. The modified groundwater module has 2 aquifers. The local aquifer, which is the shallow aquifer in the original SWAT, represents a local groundwater flow system. The regional aquifer, which replaces the deep aquifer of the original SWAT, represents intermediate and regional groundwater flow systems. Groundwater recharge is partitioned into local and regional aquifer recharges. The regional aquifer is represented by a multicell aquifer (MCA) model. The regional aquifer is discretized into cells using the Thiessen polygon method, where centres of the cells are locations of groundwater observation wells. Groundwater flow between cells is modelled using Darcy's law. Return flow from cell to stream is conceptualized using a non‐linear storage–discharge relationship. The SWAT model with the modified aquifer module, the so‐called SWAT‐MCA, was tested in 2 basins (Wipperau and Neetze) with porous aquifers in a lowland area in Lower Saxony, Germany. Results from the Wipperau basin show that the SWAT‐MCA model is able (a) to simulate baseflow in a lowland area (where baseflow is a dominant source of streamflow) better than the original model and (b) to simulate regional groundwater flow, shown by the simulated groundwater levels in cells, quite well.     

Steady state flow with hydraulic conductivity change around large diameter wells

Darcian flow law in aquifers assumes that the aquifer hydraulic conductivity is constant and the groundwater movement is due only to the piezometric level changes through hydraulic gradient. In practice, after the well development the aquifer just around the well has comparatively larger hydraulic conductivity and gradient. Patchy aquifer solutions in the literature consider sudden hydraulic conductivity changes with distance for the steady state flow. The change of transmissivity is demonstrated by the application of slope‐matching procedure to actual field data. It is the main purpose of this paper to derive simple analytical expressions for aquifer parameter evaluations with steadily decreasing hydraulic conductivity around the well. Spatial nonlinear hydraulic conductivity changes around a large‐diameter well within the depression cone of a confined aquifer are considered as exponentially decreasing functions of the radial distance. Copyright © 2010 John Wiley & Sons, Ltd.     

Is there a superior conceptual groundwater model structure for baseflow simulation?

Previous work has shown that streamflow response during baseflow conditions is a function of storage, but also that this functional relationship varies among seasons and catchments. Traditionally, hydrological models incorporate conceptual groundwater models consisting of linear or non‐linear storage–outflow functions. Identification of the right model structure and model parameterization however is challenging. The aim of this paper is to systematically test different model structures in a set of catchments where different aquifer types govern baseflow generation processes. Nine different two‐parameter conceptual groundwater models are applied with multi‐objective calibration to transform two different groundwater recharge series derived from a soil‐atmosphere‐vegetation transfer model into baseflow separated from streamflow data. The relative performance differences of the model structures allow to systematically improve the understanding of baseflow generation processes and to identify most appropriate model structures for different aquifer types. We found more versatile and more aquifer‐specific optimal model structures and elucidate the role of interflow, flow paths, recharge regimes and partially contributing storages. Aquifer‐specific recommendations of storage models were found for fractured and karstic aquifers, whereas large storage capacities blur the identification of superior model structures for complex and porous aquifers. A model performance matrix is presented, which highlights the joint effects of different recharge inputs, calibration criteria, model structures and aquifer types. The matrix is a guidance to improve groundwater model structures towards their representation of the dominant baseflow generation processes of specific aquifer types. Copyright © 2014 John Wiley & Sons, Ltd.     

Stream water infiltration, bank storage, and storage zone changes due to stream-stage fluctuations

During a flood period, stream-stage increases induce infiltration of stream water into an aquifer; subsequent declines in stream stage cause a reverse motion of the infiltrated water. This paper presents the results of the water exchange rate between a stream and aquifer, the storage volume of the infiltrated stream water in the surrounding aquifer (bank storage), and the storage zone. The storage zone is the part of aquifer where groundwater is replaced by stream water during the flood. MODFLOW was used to simulate stream–aquifer interactions and to quantify rates of stream infiltration and return flow. MODPATH was used to trace the pathlines of the infiltrated stream water and to determine the size of the storage zone. Simulations were focused on the analyses of the effects of the stream-stage fluctuation, aquifer properties, the hydraulic conductivity of streambed sediments, regional hydraulic gradients, and recharge and evapotranspiration (ET) rates on stream–aquifer interactions. Generally, for a given stream–aquifer system, larger flow rates result from larger stream-stage fluctuations; larger storage volumes and storage zones are produced by larger and longer-lasting fluctuations. For a given stream-stage hydrograph, a lower-permeable streambed, an aquitard, or an anisotropic aquifer of low vertical hydraulic conductivity can significantly reduce the rate of infiltration and limit the size of the storage zone. The bank storage solely caused by the stage fluctuation differs slightly between gaining and losing streams. Short-term rainfall recharge and ET loss in the shallow groundwater slightly influence on the flow rate, but their effects on bank storage in a larger area for a longer period can be considerable.     

Influence of Groundwater Hydraulic Gradient on Bank Storage Metrics

The hydraulic gradient between aquifers and rivers is one of the most variable properties in a river/aquifer system. Detailed process understanding of bank storage under hydraulic gradients is obtained from a two‐dimensional numerical model of a variably saturated aquifer slice perpendicular to a river. Exchange between the river and the aquifer occurs first at the interface with the unsaturated zone. The proportion of total water exchanged through the river bank compared to the river bed is a function of aquifer hydraulic conductivity, partial penetration, and hydraulic gradient. Total exchange may be estimated to within 50% using existing analytical solutions provided that unsaturated zone processes do not strongly influence exchange. Model‐calculated bank storage is at a maximum when no hydraulic gradient is present and increases as the hydraulic conductivity increases. However, in the presence of a hydraulic gradient, the largest exchange flux or distance of penetration does not necessarily correspond to the highest hydraulic conductivity, as high hydraulic conductivity increases the components of exchange both into and out of an aquifer. Flood wave characteristics do not influence ambient groundwater discharge, and so in large floods, hydraulic gradients must be high to reduce the volume of bank storage. Practical measurement of bank storage metrics is problematic due to the limitations of available measurement technologies and the nested processes of exchange that occur at the river‐aquifer interface. Proxies, such as time series concentration data in rivers and groundwater, require further development to be representative and quantitative.     

Investigating Groundwater Flow Between Edwards and Trinity Aquifers in Central Texas

Understanding the nature of communication between aquifers can be challenging when using traditional physical and geochemical groundwater sampling approaches. This study uses two multiport wells completed within Edwards and Trinity aquifers in central Texas to determine the degree of groundwater inter‐flow between adjacent aquifers. Potentiometric surfaces, hydraulic conductivities, and groundwater major ion concentrations and Sr isotope values were measured from multiple zones within three hydrostratigraphic units (Edwards and Upper and Middle Trinity aquifers). Physical and geochemical data from the multiport wells were combined with historical measurements of groundwater levels and geochemical compositions from the region to characterize groundwater flow and identify controls on the geochemical compositions of the Edwards and Trinity aquifers. Our results suggest that vertical groundwater flow between Edwards and Middle Trinity aquifers is likely limited by low permeability, evaporite‐rich units within the Upper and Middle Trinity. Potentiometric surface levels in both aquifers vary with changes in wet vs. dry conditions, indicating that recharge to both aquifers occurs through distinct recharge areas. Geochemical compositions in the Edwards, Upper, and Middle Trinity aquifers are distinct and likely reflect groundwater interaction with different lithologies (e.g., carbonates, evaporites, and siliceous sediments) as opposed to mixing of groundwater between the aquifers. These results have implications for the management of these aquifers as they indicate that, under current conditions, pumping of either aquifer will likely not induce vertical cross‐formational flow between the aquifers. Inter‐flow between the Trinity and the Edwards aquifers, however, should be reevaluated as pumping patterns and hydrogeologic conditions change.     

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.