Evaluation of Bias Associated with Capture Maps Derived from Nonlinear Groundwater Flow Models

The impact of groundwater withdrawal on surface water is a concern of water users and water managers, particularly in the arid western United States. Capture maps are useful tools to spatially assess the impact of groundwater pumping on water sources (e.g., streamflow depletion) and are being used more frequently for conjunctive management of surface water and groundwater. Capture maps have been derived using linear groundwater flow models and rely on the principle of superposition to demonstrate the effects of pumping in various locations on resources of interest. However, nonlinear models are often necessary to simulate head‐dependent boundary conditions and unconfined aquifers. Capture maps developed using nonlinear models with the principle of superposition may over‐ or underestimate capture magnitude and spatial extent. This paper presents new methods for generating capture difference maps, which assess spatial effects of model nonlinearity on capture fraction sensitivity to pumping rate, and for calculating the bias associated with capture maps. The sensitivity of capture map bias to selected parameters related to model design and conceptualization for the arid western United States is explored. This study finds that the simulation of stream continuity, pumping rates, stream incision, well proximity to capture sources, aquifer hydraulic conductivity, and groundwater evapotranspiration extinction depth substantially affect capture map bias. Capture difference maps demonstrate that regions with large capture fraction differences are indicative of greater potential capture map bias. Understanding both spatial and temporal bias in capture maps derived from nonlinear groundwater flow models improves their utility and defensibility as conjunctive‐use management tools.     

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.     

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

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

Analytical Solutions for Flow Fields near Drain-and-Gate Reactive Barriers

Permeable reactive barriers (PRBs) are a popular technology for passive contaminant remediation in aquifers through installation of reactive materials in the pathway of a plume. Of fundamental importance are the degree of remediation inside the reactor (residence time) and the portion of groundwater intercepted by a PRB (capture width). Based on a two-dimensional conformal mapping approach (previously used in related work), the latter is studied in the present work for drain-and-gate (DG) PRBs, which may possess a collector and a distributor drain (“full” configuration) or a collector drain only (“simple” configuration). Inherent assumptions are a homogeneous unbounded aquifer with a uniform far field, in which highly permeable drains establish constant head boundaries. Solutions for aquifer flow fields in terms of the complex potential are derived, illustrated, and analyzed for doubly symmetric DG configurations and arbitrary reactor hydraulic resistance as well as ambient groundwater flow direction. A series of practitioner-friendly charts for capture width is given to assist in PRB design and optimization without requiring complex mathematics. DG PRBs are identified as more susceptible to flow divergence around the reactor than configurations using impermeable side structures (e.g., funnel-and-gate), and deployment of impermeable walls on drains is seen to mitigate this problem under certain circumstances.     

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.     

Water Scarcity: Moving Beyond Indexes to Innovative Institutions

Water scarcity is a media darling often times described as a trigger of conflict in arid regions, a by‐product of human influences ranging from desertification to climate change, or a combination of natural‐ and human‐induced changes in the water cycle. A multitude of indexes have been developed over the past 20 years to define water scarcity to map the “problem” and guide international donor investment. Few indexes include groundwater within the metrics of “scarcity.” Institutional communication contributes to the recognition of local or regional water scarcity. However, evaluations that neglect groundwater resources may incorrectly define conditions as scarce. In cases where there is a perception of scarcity, the incorporation of groundwater and related storage in aquifers, political willpower, new policy tools, and niche diplomacy often results in a revised status, either reducing or even eliminating the moniker locally. Imaginative conceptualization and innovative uses of aquifers are increasingly used to overcome water scarcity.     

Temporal sampling and role of flux measurements for subsurface heterogeneous characterization in groundwater basins using hydraulic tomography

Accurate characterization of heterogeneity in groundwater basins is crucial to the sustainable management of groundwater resources. This study explores the temporal sampling issues and the role of flux measurements in the characterization of heterogeneity in groundwater basins using numerical experiments. The experiments involve a digital basin imitating the groundwater basin of the North China Plain (NCP), where the groundwater exploitation reduction program is ongoing. Using the experiments, we champion that the reduction program could collect groundwater level information induced by operational variations of existing pumping wells at different locations in the basin. Such a dataset could serve as a basin-scale hydraulic tomography (HT) to characterize the basin-scale heterogeneity cost-effectively. Both steady-state and transient-state inversion experiments demonstrate the advantage of HT surveys in characterizing basin-scale heterogeneity over conventional pumping tests at fixed well locations. Additionally, head data at the early, intermediate, and late time from well hydrographs should be selected for the HT analysis to maximize HT's power and save computational costs. When accurate geological zones are incorporated in prior information, flux measurements significantly improve parameter estimates based on conventional pumping tests. However, their effects are less noticeable for long-term HT surveys in such basin-scale aquifers without fissures or fractures. This basin-scale tomographic survey example serves a guide for field data collection and optimization of the analysis of future basin-scale HT.     

Tracking Groundwater Discharge to a Large River using Tracers and Geophysics

Few studies have investigated large reaches of rivers in which multiple sources of groundwater are responsible for maintaining baseflow. This paper builds upon previous work undertaken along the Fitzroy River, one of the largest perennial river systems in north‐western Australia. Synoptic regional‐scale sampling of both river water and groundwater for a suite of environmental tracers (⁴He, ⁸⁷Sr/⁸⁶Sr, ²²²Rn and major ions), and subsequent modeling of tracer behavior in the river, has enabled definition and quantification of groundwater input from at least three different sources. We show unambiguous evidence of both shallow “local” groundwater, possibly recharged to alluvial aquifers beneath the adjacent floodplain during recent high‐flow events, and old “regional” groundwater introduced via artesian flow from deep confined aquifers. We also invoke hyporheic exchange and either bank return flow or parafluvial flow to account for background ²²²Rn activities and anomalous chloride trends along river reaches where there is no evidence of the local or regional groundwater inputs. Vertical conductivity sections acquired through an airborne electromagnetic (AEM) survey provide insights to the architecture of the aquifers associated with these sources and general groundwater quality characteristics. These data indicate fresh groundwater from about 300 m below ground preferentially discharging to the river, at locations consistent with those inferred from tracer data. The results demonstrate how sampling rivers for multiple environmental tracers of different types—including stable and radioactive isotopes, dissolved gases and major ions—can significantly improve conceptualization of groundwater—surface water interaction processes, particularly when coupled with geophysical techniques in complex hydrogeological settings.     

Groundwater modelling and marine intrusion in the semi‐arid systems (Cap‐Bon,Tunisia)

This article studies the effect of drought and pumping discharge on groundwater supplies and marine intrusion in the Korba aquifer (Cap‐Bon peninsula, Tunisia). The Groundwater Modelling System has been used to model the groundwater flow and to simulate the seawater intrusion. The calibration is based on the groundwater levels in the steady state from 1963, and in the transient state from the groundwater levels from 1963 to 2005. The main objective is to quantify the components of the groundwater mass balance and to estimate the hydraulic conductivity distribution. The impact of pumping discharge on the groundwater level evolution has been examined by two pumping scenarios P1 (no. 8420) and P2 (no. 8862) wells. The hydrodynamic modelling shows the increasing drawdowns after 14 years of pumping: 4 m in P1 well and about 5 m in P2 well below sea level. The drawdowns are accompanied by the inverse hydraulic gradient. The numerical model was used to discuss the management of the groundwater resources of Cap‐Bon. As the population continues to grow and the demand for groundwater pumping intensifies beyond the 1963 level, it can be expected that the actual extent of seawater intrusion in the future would be more severe than the model prediction. Better strategies for groundwater development and management will be necessary to protect the freshwater aquifers to the marine intrusion. Copyright © 2011 John Wiley & Sons, Ltd.     

Solving Groundwater Flow Problems with Time Series Analysis: You May Not Even Need Another Model

Time series analysis is a data-driven approach to analyze time series of heads measured in an observation well. Time series models are commonly much simpler and give much better fits than regular groundwater models. Time series analysis with response functions gives insight into why heads vary, while such insight is difficult to gain with black box models out of the artificial intelligence world. An important application is to quantify the contributions to the head variation of different stresses on the aquifer, such as rainfall and evaporation, pumping, and surface water levels. Time series analysis may be applied to answer many groundwater questions without the need for a regular groundwater model, such as what is the drawdown caused by a pumping station? Or, how long will it take before groundwater levels recover after a period of drought? Even when a regular groundwater model is needed to solve a groundwater problem, time series analysis can be of great value. It can be used to clean up the data, identify the major stresses on the aquifer, determine the most important processes that affect flow in the aquifer, and give an indication of the fit that can be expected. In addition, it can be used to determine calibration targets for steady-state models, and it can provide several alternative calibration methods for transient models. In summary, the overarching message of this paper is that it would be wise to do time series analysis for any application that uses measured groundwater heads.     

GIS integration of hydrogeological and geoelectrical data for groundwater potential modeling in the western part of greater Kushtia district of Bangladesh

Geographical Information System (GIS) has been used in this paper to delineate groundwater resources potential in the western part of greater Kushtia district of Bangladesh, where urgent attention for augmentation of irrigation water supply is required. Thematic maps of transmissivity, hydraulic conductivity, specific yield, net recharge, aquifer thickness, surface water bodies, aquifer resistivity, overburden aquitard thickness and its resistivity have been prepared and assigned weight according to their relative importance using Analytical Hierarchical Process for the preparation of groundwater potential model. Since the values within each thematic map vary significantly, they are classified into various ranges or types and assigned ratings. Finally, the thematic maps are integrated using GIS to prepare the groundwater potential map for the study area in terms of Ground Water Potential Index (GWPI). The evolved map indicates that 22.51% of the study area have GWPI more than 0.70 and therefore, have excellent prospective for exploitation. About 69.12% of the area with GWPI ranging from 0.50 to 0.70 is also quite promising for groundwater abstraction, while the rest 8.37% area having GWPI below 0.50 indicates moderate potential. The obtained map of groundwater potential is found in good agreement with the yields of available pumping test data.     

Conceptualizing a multi-layered shingle aquifer model based on volcanic stratigraphy and water inflow to lava caves in Jeju Island,Korea

Volcanic aquifers supply a substantial portion of water resources in many parts of the world, including islands, and their productivity depends strongly on volcanic stratigraphy, which exhibits considerable heterogeneity. We investigated water inflow to lava tube caves formed from numerous basaltic lava flows in the northeastern coastal area of Jeju Island after storm events and monitored relative inflow rates monthly over 1 year to characterize groundwater flow processes in the upper parts of volcanic aquifers, and to evaluate the applicability of the previous hydrogeological models proposed for the island. Considerable water inflow arose shortly after storms from exposed palaeosol layers on the walls of the caves. The monthly monitoring results showed that wall inflow associated with these palaeosol layers is substantial. In both cases, discharge from ceiling drips was much less and more temporally variable compared to wall inflow discharge. Water flowing into the caves was rapidly drained through the floor at all monitoring sites. The lateral extent of the palaeosol layers was identified using drill core logs near the cave and outcrops in the coastal area. Based on these results, we inferred that multiple perched aquifers are formed by low-permeability palaeosol layers between lava flows, which are connected by vertical flows at discontinuities in the palaeosol layer, eventually reaching the basal aquifer. This study revealed the water inflow processes observed in lava tube caves constrained by palaeosol layers, and established a hydrogeological conceptual model incorporating multiple perched aquifers in both coastal and mountainous areas associated with extensive palaeosol layers formed during volcanic hiatuses. This finding would help elucidate recharge, groundwater flow, and contaminant transport processes in many volcanic aquifers that are not adequately represented by the previous models, and contribute to better management of groundwater in those areas.     

Assessment of groundwater resources of Quaternary aquifer system in concordance with avoidable negative impact on geoenvironment,south-eastern part of Lithuania

This paper focuses upon south-eastern Lithuanian Quaternary aquifer system groundwater resources formation modelling. Groundwater model calibration has been performed for a pre-development and transient flow conditions. The results demonstrate that there is an intense interaction between groundwater and surface water bodies which form groundwater resources and runoff. For Quaternary cover the majority of unconfined groundwater outflows to surface water streams, the remaining part discharges through the confined interglacial/interstadial aquifers and lateral outflows across the boundaries. Groundwater prognostic exploitable resources can be obtained without a significant negative impact on the geoenvironment. The main sources of exploitable resources formation are increase in groundwater recharge and lateral inflow (40.7%), decreased outflow via streams and lateral boundaries (41.2%) and drawing up additional flow from the hydraulically connected streams (14.5%).     

Contribution of baseflow nitrate export to non-point source pollution

As a common pollutant of nitrogen in groundwater, nitrate contamination has become a major concern worldwide. Baseflow, one of the dominant hydrological pathways for nitrate migration to streamflow, has been confirmed as a leading nitrate source for stream water where groundwater or subsurface flow contaminated heavily by nitrate. That is, sufficient improvements of water quality may not be attained without proper management for baseflow, even if non-point sources (NPS) pollutants discharged through surface runoff are being well managed. This article reviews the primary nitrate sources, the main factors affecting its transport, and the methodologies for baseflow nitrate estimation, to give some recommendations for future works, including: (1) giving sufficient consideration for the effects of climatological, morphological, and geological factors on baseflow recessions to obtain more reliable and accurate baseflow separation; (2) trying to solve calibration and validation problems for baseflow loads determining in storm flow period; (3) developing a simple and convenient algorithm with certain physics that can be used to separate baseflow NPS pollution from the total directly in different regions, for a reliable estimation of baseflow NPS pollution at larger scale (e.g., national scale); (4) improving groundwater quality simulation module of existing NPS pollution models to have a better simulation for biogeochemical processes in shallow aquifers; (5) taking integrated measures of “source control”, “process interception” and “end remediation” to prevent and control NPS nitrate pollution effectively, not just only the strict control of nutrients loss from surface runoff.     

The unusual and large drawdown response of buried-valley aquifers to pumping

The buried-valley aquifers that are common in the glacial deposits of the northern hemisphere are a typical case of the strip aquifers that occur in many parts of the world. Pumping from a narrow strip aquifer leads to much greater drawdown and much more distant drawdown effects then would occur in a sheet aquifer with a similar transmissivity and storage coefficient. Widely used theories for radial flow to wells, such as the Theis equation, are not appropriate for narrow strip aquifers. Previously published theory for flow to wells in semiconfined strip aquifers is reviewed and a practical format of the type curves for pumping-test analysis is described. The drawdown response of strip aquifers to pumping tests is distinctive, especially for observation wells near the pumped well. A case study is presented, based on extensive pumping test experience for the Estevan Valley Aquifer in southern Saskatchewan, Canada. Evaluation of groundwater resources in such buried-valley aquifers needs to take into account the unusually large drawdowns in response to pumping.     

Groundwater recharge estimation in arid hardrock‐alluvium aquifers using combined water‐table fluctuation and groundwater balance approaches

This paper proposes an approach to estimate groundwater recharge using an optimization‐based water‐table fluctuation method combined with a groundwater balance model in an arid hardrock‐alluvium region, located at the Oman–United Arab Emirates border. We introduce an “effective hardrock thickness” term to identify the percentage of the considered hardrock thickness in which effective groundwater flow takes place. The proposed method is based upon a Thiessen polygon zoning approach. The method includes subpolygons to represent specific geologic units and to enhance the confidence of the estimated groundwater recharge. Two linear and 1 nonlinear submodels were developed to evaluate the model components for the calibration (October 1996 to September 2008) and validation (October 2008 to September 2013) periods. Long‐term annual groundwater recharge from rainfall and return flow over the model domain are estimated as 24.62 and 5.71 Mm³, respectively, while the effective groundwater flow circulation is found to occur in the upper 7% of the known hardrock thickness (42 m), confirming conclusions of previous field studies. Considering a total difference in groundwater levels between eastern and western points of the study area of the order of 220 m and a 12‐year monthly calibration period, a weighted root mean squared error in predicted groundwater elevation of 2.75 m is considered quite reasonable for the study area characterized by remarkable geological and hydrogeological diversity. The proposed approach provides an efficient and robust method to estimate groundwater recharge in regions with a complex geological setting in which interaction between fractured and porous media cannot be easily assessed.     

Hydraulic Conductivity Calibration of Logging NMR in a Granite Aquifer,Laramie Range,Wyoming

In granite aquifers, fractures can provide both storage volume and conduits for groundwater. Characterization of fracture hydraulic conductivity (K) in such aquifers is important for predicting flow rate and calibrating models. Nuclear magnetic resonance (NMR) well logging is a method to quickly obtain near-borehole hydraulic conductivity (i.e., K_NMR) at high-vertical resolution. On the other hand, FLUTe flexible liner technology can produce a K profile at comparable resolution but requires a fluid driving force between borehole and formation. For three boreholes completed in a fractured granite, we jointly interpreted logging NMR data and FLUTe K estimates to calibrate an empirical equation for translating borehole NMR data to K estimates. For over 90% of the depth intervals investigated from these boreholes, the estimated K_NMR are within one order of magnitude of K_FLUTe. The empirical parameters obtained from calibrating the NMR data suggest that “intermediate diffusion” and/or “slow diffusion” during the NMR relaxation time may occur in the flowing fractures when hydraulic aperture are sufficiently large. For each borehole, “intermediate diffusion” dominates the relaxation time, therefore assuming “fast diffusion” in the interpretation of NMR data from fractured rock may lead to inaccurate K_NMR estimates. We also compare calibrations using inexpensive slug tests that suggest reliable K_NMR estimates for fractured rock may be achieved using limited calibration against borehole hydraulic measurements.     

Viruses as Groundwater Tracers: Using Ecohydrology to Characterize Short Travel Times in Aquifers

Viruses are attractive tracers of short (<3 year) travel times in aquifers because they have unique genetic signatures, are detectable in trace quantities, and are mobile in groundwater. Virus “snaphots” result from infection and disappearance in a population over time; therefore, the virus snapshot shed in the fecal wastes of an infected population at a specific point in time can serve as a marker for tracking virus and groundwater movement. The virus tracing approach and an example application are described to illustrate their ability to characterize travel times in high‐groundwater velocity settings, and provide insight unavailable from standard hydrogeologic approaches. Although characterization of preferential flowpaths does not usually characterize the majority of other travel times occurring in the groundwater system (e.g., center of plume mass; tail of the breakthrough curve), virus approaches can trace very short times of transport, and thus can fill an important gap in our current hydrogeology toolbox.     

Extending the Capture Map Concept to Estimate Discrete and Risk-Based Streamflow Depletion Potential

A popular and contemporary use of numerical groundwater models is to estimate the discrete relation between groundwater extraction and surface-water/groundwater exchange. Previously, the concept of a “capture map” has been put forward as a means to effectively summarize this relation for decision-making consumption. While capture maps have enjoyed success in the environmental simulation industry, they are deterministic, ignoring uncertainty in the underlying model. Furthermore, capture maps are not typically calculated in a manner that facilitates analysis of varying combinations of extraction locations and/or reaches. That is, they are typically constructed with focus on a single reach or group of reaches. The former of these limitations is important for conveying risk to decision makers and stakeholders, while the latter is important for decision-making support related to surface-water management, where future foci may include reaches that were not the focus of the original capture analysis. Herein, we use the concept of a response matrix to generalize the theory of the capture-map approach to estimate spatially discrete streamflow depletion potential. We also use first-order, second-moment uncertainty estimation techniques with the concept of “risk shifting” to place capture maps and streamflow depletion potential in a stochastic, risk-based framework. Our approach is demonstrated for an integrated groundwater/surface-water model of the lower San Antonio River, Texas, USA.     

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.