skip to main content

Explore Research Products in the NSF-PAR

  • NSF Public Access
  • Search Results
  • Decadal‐Scale Aquifer Dynamics and Structural Complexities at a Municipal Wellfield Revealed by 25 Years of InSAR and Recent Groundwater Temperature Observations
Title: Decadal‐Scale Aquifer Dynamics and Structural Complexities at a Municipal Wellfield Revealed by 25 Years of InSAR and Recent Groundwater Temperature Observations
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.

  more » « less
Award ID(s):
1301346 1830172
NSF-PAR ID:
10372995
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Water Resources Research
Volume:
55
Issue:
12
ISSN:
0043-1397
Page Range / eLocation ID:
p. 10636-10656
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; ; ; ; ; et al ( , Seismological Research Letters)
    null (Ed.)
    Abstract The Mw 7.1 47 km deep earthquake that occurred on 30 November 2018 had deep societal impacts across southcentral Alaska and exhibited phenomena of broad scientific interest. We document observations that point to future directions of research and hazard mitigation. The rupture mechanism, aftershocks, and deformation of the mainshock are consistent with extension inside the Pacific plate near the down‐dip limit of flat‐slab subduction. Peak ground motions >25%g were observed across more than 8000  km2, though the most violent near‐fault shaking was avoided because the hypocenter was nearly 50 km below the surface. The ground motions show substantial variation, highlighting the influence of regional geology and near‐surface soil conditions. Aftershock activity was vigorous with roughly 300 felt events in the first six months, including two dozen aftershocks exceeding M 4.5. Broad subsidence of up to 5 cm across the region is consistent with the rupture mechanism. The passage of seismic waves and possibly the coseismic subsidence mobilized ground waters, resulting in temporary increases in stream flow. Although there were many failures of natural slopes and soils, the shaking was insufficient to reactivate many of the failures observed during the 1964 M 9.2 earthquake. This is explained by the much shorter duration of shaking as well as the lower amplitude long‐period motions in 2018. The majority of observed soil failures were in anthropogenically placed fill soils. Structural damage is attributed to both the failure of these emplaced soils as well as to the ground motion, which shows some spatial correlation to damage. However, the paucity of instrumental ground‐motion recordings outside of downtown Anchorage makes these comparisons challenging. The earthquake demonstrated the challenge of issuing tsunami warnings in complex coastal geographies and highlights the need for a targeted tsunami hazard evaluation of the region. The event also demonstrates the challenge of estimating the probabilistic hazard posed by intraslab earthquakes. 
    more » « less
  2. ; ; ; ( , Geosphere)
    Abstract Uturuncu volcano in southern Bolivia last erupted around 250 ka but is exhibiting signs of recent activity, including over 50 yr of surface uplift, elevated seismic activity, and fumarolic activity. We studied the spatial and temporal scales of surface deformation from 1992 to 2021 to better understand subsurface activity. We tracked Uturuncu’s recent deformation using interferometric synthetic aperture radar (InSAR) data and the global navigation satellite system (GNSS) station UTUR, located near Uturuncu’s summit. We observed a spatially coherent signal of uplift from 2014 to 2021 from Sentinel-1 A/B satellites that indicates the Altiplano-Puna magma body, located 19–24 km below ground level, and previously noted as the source of the large region of deformation, is still active. The ground is now uplifting at a rate of ~3 mm/yr compared to prior rates of ~10 mm/yr. We corroborated this waning uplift with in situ data from station UTUR. We combined the Sentinel-1 data with TerraSAR-X interferograms to constrain an ~25 km2 region of subsidence located 11 km SSW of Uturuncu, with a source depth of 2.1 km below ground level to an active period of ~2.5 yr with ~5 mm/yr subsidence. We developed a conceptual model that relates these varying depths and time scales of activity in a transcrustal magmatic system. We associate the surface uplift with pressurization from ascending gases and brines from magmatic reservoirs in the midcrust. We infer the existence of brine lenses in the shallow hydrothermal system based on low subsurface resistivity correlated with surface subsidence. 
    more » « less
  3. ; ; ; ; ( , Water Resources Research)
    Abstract

    Mountain‐front recharge (MFR), or all inflow to a basin‐fill aquifer with its source in the mountain block, is an important component of recharge to basin‐fill aquifer systems. Distinguishing and quantifying the surface from subsurface components of MFR is necessary for water resource planning and management, particularly as climate change may impact these components in distinct ways. This study tests the hypothesis that MFR components can be distinguished in long‐screened, basin‐fill production wells by (1) groundwater age and (2) the median elevation of recharge. We developed an MFR characterization approach by combining age distributions in six wells using tritium, krypton‐85, argon‐39, and radiocarbon, and median recharge elevations from noble gas thermometry combined with numerical experiments to determine recharge temperature lapse rates using flow and energy transport modeling. We found that groundwater age distributions provided valuable information for characterizing the dominant flow system behavior captured by the basin‐fill production wells. Tracers indicated the presence of old (i.e., no detectable tritium) water in a well completed in weathered bedrock located close to the mountain front. Two production wells exhibited age distributions of binary mixing between modern and a small fraction of old water, whereas the remaining wells captured predominantly modern flow paths. Noble gas thermometry provided important complementary information to the age distributions; however, assuming constant recharge temperature lapse rates produced improbable recharge elevations. Numerical experiments suggest that surface MFR, if derived from snowmelt, can locally suppress water table temperatures in the basin‐fill aquifer, with implications for recharge elevations estimated from noble gas thermometry.

     
    more » « less
  4. ; ; ; ; ; ; ( , Journal of Hydrometeorology)
    null (Ed.)
    Abstract Some of the most intense convective storms on Earth initiate near the Sierras de Córdoba mountain range in Argentina. The goal of the RELAMPAGO field campaign was to observe these intense convective storms and their associated impacts. The intense observation period (IOP) occurred during November–December 2018. The two goals of the hydrometeorological component of RELAMPAGO IOP were 1) to perform hydrological streamflow and meteorological observations in previously ungauged basins and 2) to build a hydrometeorological modeling system for hindcast and forecast applications. During the IOP, our team was able to construct the stage–discharge curves in three basins, as hydrological instrumentation and personnel were successfully deployed based on RELAMPAGO weather forecasts. We found that the flood response time in these river locations is typically between 5 and 6 h from the peak of the rain event. The satellite-observed rainfall product IMERG-Final showed a better representation of rain gauge–estimated precipitation, while IMERG-Early and IMERG-Late had significant positive bias. The modeling component focuses on the 48-h simulation of an extreme hydrometeorological event that occurred on 27 November 2018. Using the Weather Research and Forecasting (WRF) atmospheric model and its hydrologic component WRF-Hydro as an uncoupled hydrologic model, we developed a system for hindcast, deterministic forecast, and a 60-member ensemble forecast initialized with regional-scale atmospheric data assimilation. Critically, our results highlight that streamflow simulations using the ensemble forecasting with data assimilation provide realistic flash flood forecast in terms of timing and magnitude of the peak. Our findings from this work are being used by the water managers in the region. 
    more » « less
  5. ; ; ; ; ; ; ( , Geophysical Research Letters)
    Abstract

    In 2013 and 2018, earthquake swarms with a maximum moment magnitude of 4.5 occurred ~5 km from the northern section of the Dead Sea Transform Fault. Here we show that aquifer pressure data, interferometric synthetic aperture radar surface deformation time series, and seismic monitoring suggest that groundwater withdrawal triggered these earthquakes. Continuous groundwater extraction from several wells located ~10 km west of the swarms has accelerated since 2010 and resulted in a total decrease of ~50 m of the groundwater level at the time of the 2018 earthquake swarm. The withdrawal also corresponds to surface subsidence of ~10 mm/year based on repeat interferometric synthetic aperture radar measurements. The temporal correlation, extensive subsidence, anomalous swarm characteristics, and normal faulting orientation suggest a connection between the groundwater withdrawal and recent earthquakes. Poroelastic modeling demonstrates that pumping‐induced pore pressure decrease west of the earthquake could have caused significant dilatational stresses that led to normal faulting events outside the aquifer.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email