Many irrigated agricultural areas seek to prolong the lifetime of their groundwater resources by reducing pumping. However, it is unclear how lagged responses, such as reduced groundwater recharge caused by more efficient irrigation, may impact the long‐term effectiveness of conservation initiatives. Here, we use a variably saturated, simplified surrogate groundwater model to: (a) analyze aquifer responses to pumping reductions, (b) quantify time lags between reductions and groundwater level responses, and (c) identify the physical controls on lagged responses. We explore a range of plausible model parameters for an area of the High Plains aquifer (USA) where stakeholder‐driven conservation has slowed groundwater depletion. We identify two types of lagged responses that reduce the long‐term effectiveness of groundwater conservation, recharge‐dominated and lateral‐flow‐dominated, with vertical hydraulic conductivity (
- Award ID(s):
- 2019561
- NSF-PAR ID:
- 10252185
- Date Published:
- Journal Name:
- Frontiers in Water
- Volume:
- 3
- ISSN:
- 2624-9375
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract K Z) the major controlling variable. When highK Zallows percolation to reach the aquifer, more efficient irrigation reduces groundwater recharge. By contrast, when lowK Zimpedes vertical flow, short term changes in recharge are negligible, but pumping reductions alter the lateral flow between the groundwater conservation area and the surrounding regions (lateral‐flow‐dominated response). For the modeled area, we found that a pumping reduction of 30% resulted in median usable lifetime extensions of 20 or 25 years, depending on the dominant lagged response mechanism (recharge‐ vs. lateral‐flow‐dominated). These estimates are far shorter than estimates that do not account for lagged responses. Results indicate that conservation‐based pumping reductions can extend aquifer lifetimes, but lagged responses can create a sizable difference between the initially perceived and actual long‐term effectiveness. -
more » « less
Abstract Across South Asia, millions of villagers have reduced their exposure to high‐arsenic (As) groundwater by switching to low‐As wells. Isotopic tracers and flow modeling are used in this study to understand the groundwater flow system of a semi‐confined aquifer of Pleistocene (>10 kyr) age in Bangladesh that is generally low in As but has been perturbed by massive pumping at a distance of about 25 km for the municipal water supply of Dhaka. A 10‐ to 15‐m‐thick clay aquitard caps much of the intermediate aquifer (>40‐ to 90‐m depth) in the 3‐km2study area, with some interruptions by younger channel sand deposits indicative of river scouring. Hydraulic heads in the intermediate aquifer below the clay‐capped areas are 1–2 m lower than in the high‐As shallow aquifer above the clay layer. In contrast, similar heads in the shallow and intermediate aquifer are observed where the clay layer is missing. The head distribution suggests a pattern of downward flow through interruptions in the aquitard and lateral advection from the sandy areas to the confined portion of the aquifer. The interpreted flow system is consistent with3H‐3He ages, stable isotope data, and groundwater flow modeling. Lateral flow could explain an association of elevated As with high methane concentrations within layers of gray sand below certain clay‐capped portions of the Pleistocene aquifer. An influx of dissolved organic carbon from the clay layer itself leading to a reduction of initially orange sands has also likely contributed to the rise of As.
-
more » « less
Abstract Over the past 35 years the Buckman wellfield near Santa Fe, New Mexico, experienced production well drawdowns in excess of 180 m, resulting in ground subsidence and surface cracks. Increased reliance on surface water diversions since 2011 has reduced pumping and yielded water level recovery. To characterize the impact of wellfield management decisions on the aquifer system, we reconstruct the surface deformation history through the European Remote Sensing Satellite, Advanced Land Observing Satellite, and Sentinel‐1 Interferometric Synthetic Aperture Radar (InSAR) time series analysis during episodes of drawdown (1993–2000), recovery (2007–2010), and modern management (2015–2018) in discontinuous observations over a 25‐year period. The observed deformation generally reflects changes in hydraulic head. However, at times during the wellfield recovery, the deformation signal is complex, with patterns of uplift and subsidence suggesting a compartmentalized aquifer system. Recent records of locally high geothermal gradients and an overall warming of the system (~0.5°C during the water level recovery) obtained from repeat temperature measurements between 2013 and 2018 constrain a conceptual model of convective heat transfer that requires a vertical permeable zone near an observed fault. To reproduce observed temperature patterns at monitoring wells, high basal heat flow and convective cooling associated with downward flow of water from cool shallow aquifers during the drawdown period is necessary. The fault, however, appears to die out southward or may be locally permeable, as conceptual cross‐sectional hydrologic modeling reproduces the surface deformation without such a structure. Our work demonstrates the importance of incorporating well‐constrained stratigraphy and structure when modeling near‐surface deformation induced by, for instance, groundwater production.
-
Abstract Non-perennial rivers and streams make up over half the global river network and are becoming more widespread. Transitions from perennial to non-perennial flow are a threshold-type change that can lead to alternative stable states in aquatic ecosystems, but it is unknown whether streamflow itself is stable in either wet (flowing) or dry (no-flow) conditions. Here, we investigated drivers and feedbacks associated with regime shifts between wet and dry conditions in an intermittent reach of the Arkansas River (USA) over the past 23 years. Multiple lines of evidence suggested that these regimes represent alternative stable states, including (a) significant jumps in discharge time series that were not accompanied by jumps in flow drivers such as precipitation and groundwater pumping; (b) a multi-modal state distribution with 92% of months experiencing no-flow conditions for <10% or >90% of days, despite unimodal distributions of precipitation and pumping; and (c) a hysteretic relationship between climate and flow state. Groundwater levels appear to be the primary control over the hydrological regime, as groundwater levels in the alluvial aquifer were higher than the stream stage during wet regimes and lower than the streambed during dry regimes. Groundwater level variation, in turn, was driven by processes occurring at both the regional scale (surface water inflows from upstream, groundwater pumping) and the reach scale (stream–aquifer exchange, diffuse recharge through the soil column). Historical regime shifts were associated with diverse pressures including network disconnection caused by upstream water use, increased flow stability potentially associated with reservoir operations, and anomalous wet and dry climate conditions. In sum, stabilizing feedbacks among upstream inflows, stream–aquifer interactions, climate, vegetation, and pumping appear to create alternative wet and dry stable states at this site. These stabilizing feedbacks suggest that widespread observed shifts from perennial to non-perennial flow will be difficult to reverse.more » « less
-
In this study, a stress-dependent groundwater model, MODFLOW-SD, has been developed and coupled with the nonlinear subsidence model, NDIS, to predict vertical deformation occurring in basins with highly compressible deposits. The MODFLOW-SD is a modified version of MODFLOW (the USGS Modular Three-Dimensional Groundwater Flow Model) with two new packages, NONK and NONS, to update hydraulic conductivity and skeletal specific storage due to change in effective stress. The NDIS package was developed based on Darcy–Gersevanov Law and bulk flux to model land subsidence. Results of sample simulations run for a conceptual model showed that hydraulic heads calculated by MODFLOW significantly overestimated for confining units and slightly underestimated for aquifer ones. Moreover, it showed that applied stress due to pumping changed initially homogeneous layers to be heterogeneous ones. Comparison of vertical deformations calculated by NDIS andMODFLOW-SUB showed that neglecting horizontal strain and stress-dependency of aquifer parameters can overestimate future subsidence. Furthermore, compared to the SUB (Subsidence and Aquifer-System Compaction) package, NDIS is more likely to provide a more accurate compaction model for a complex aquifer system with vertically variable compression (Cc), recompression (Cr), and hydraulic conductivity change (Ck) indices.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3389/frwa.2021.623476
