Aquifer Mechanical Properties and Decelerated Compaction in Tucson, Arizona

In recent decades, high groundwater extraction rates, often coincident with periods of severe drought, result in the widespread decline of water levels. Overexploitation of aquifers also causes land subsidence, which poses a severe threat to infrastructure. Tucson, Arizona experiences land subsidence coupled with the depletion of groundwater, a critical water resource for the desert city. To understand the spatiotemporal evolution of land subsidence and its implications for aquifer properties, we examine long time series of surface deformation and head levels. Measurements at extensometer stations indicate rapid compaction of fine‐grained material up to 8.5 mm/yr from 1990 to 2005, which results in permanent storage volume losses up to 4.1%. The analysis of densely populated sets of interferograms generated from Envisat and RadarSAT C band acquisitions yields multitemporal maps of surface deformation at unprecedented resolution. These maps reveal that subsidence significantly slows by the late 2000s, corresponding with the implementation of artificial recharge efforts. Subsequent to groundwater level recovery, we observe a brief 6.6 year interval of residual compaction, suggesting a high vertical hydraulic conductivity, which is then shown to be up to 9.8 × 10−4 m/d. We also estimate the average elastic and inelastic skeletal storage coefficients for the aquifer system to be 3.78 × 10−3 and 6.01 × 10−3, respectively. Interferometric synthetic aperture radar shows deformation nearly ceases by 2015, likely reducing hazards associated with Earth fissuring and infrastructure damage. This study highlights successful outcomes of water management and conservation plans that preserve existing groundwater reserves and increase artificial recharge.

[1]  E. Fielding,et al.  Applicability of Sentinel‐1 Terrain Observation by Progressive Scans multitemporal interferometry for monitoring slow ground motions in the San Francisco Bay Area , 2017 .

[2]  T. Farr,et al.  Estimating the permanent loss of groundwater storage in the southern San Joaquin Valley, California , 2017 .

[3]  Willem A. Schreüder,et al.  Confined aquifer head measurements and storage properties in the San Luis Valley, Colorado, from spaceborne InSAR observations , 2016 .

[4]  Kristine Uhlman,et al.  Enhancing drought resilience with conjunctive use and managed aquifer recharge in California and Arizona , 2016 .

[5]  B. Conway Land subsidence and earth fissures in south-central and southern Arizona, USA , 2016, Hydrogeology Journal.

[6]  Di Long,et al.  Hydrologic implications of GRACE satellite data in the Colorado River Basin , 2015 .

[7]  Sharon B. Megdal,et al.  How a Drought‐Resilient Water Delivery System Rose Out of the Desert: The Case of Tucson Water , 2015 .

[8]  Zhong Lu,et al.  Ground subsidence in Tucson, Arizona, monitored by time-series analysis using multi-sensor InSAR datasets from 1993 to 2011 , 2015 .

[9]  Manoochehr Shirzaei,et al.  Spatiotemporal characterization of land subsidence and uplift in Phoenix using InSAR time series and wavelet transforms , 2015 .

[10]  Matthew Rodell,et al.  Groundwater depletion during drought threatens future water security of the Colorado River Basin , 2014, Geophysical research letters.

[11]  E. Fielding,et al.  Predictability of hydraulic head changes and characterization of aquifer‐system and fault properties from InSAR‐derived ground deformation , 2014 .

[12]  W. Cunningham,et al.  Groundwater Technical Procedures of the U.S. Geological Survey , 2014 .

[13]  Peter K. Kitanidis,et al.  Estimating temporal changes in hydraulic head using InSAR data in the San Luis Valley, Colorado , 2014 .

[14]  Felix W. Landerer,et al.  Seasonal variation in total water storage in California inferred from GPS observations of vertical land motion , 2014 .

[15]  G. Preisig,et al.  Regional Flow and Deformation Analysis of Basin‐Fill Aquifer Systems Using Stress‐Dependent Parameters , 2014, Ground water.

[16]  Manoochehr Shirzaei,et al.  Time‐dependent model of creep on the Hayward fault from joint inversion of 18 years of InSAR and surface creep data , 2013 .

[17]  Manoochehr Shirzaei,et al.  A Wavelet-Based Multitemporal DInSAR Algorithm for Monitoring Ground Surface Motion , 2013, IEEE Geoscience and Remote Sensing Letters.

[18]  D. Argus Uncertainty in the velocity between the mass center and surface of Earth , 2012 .

[19]  Manoochehr Shirzaei,et al.  Topography correlated atmospheric delay correction in radar interferometry using wavelet transforms , 2012 .

[20]  T. Burbey,et al.  Review: Regional land subsidence accompanying groundwater extraction , 2011 .

[21]  Thomas R. Walter,et al.  Estimating the Effect of Satellite Orbital Error Using Wavelet-Based Robust Regression Applied to InSAR Deformation Data , 2011, IEEE Transactions on Geoscience and Remote Sensing.

[22]  Stanley A. Leake,et al.  Trends in groundwater levels in wells in the active management areas of Arizona, USA , 2010 .

[23]  Empirical Modeling of Solar Radiation Pressure Forces Affecting GPS Satellites , 2010 .

[24]  Thomas R. Walter,et al.  Randomly iterated search and statistical competency as powerful inversion tools for deformation source modeling: Application to volcano interferometric synthetic aperture radar data , 2009 .

[25]  C. Faunt,et al.  Groundwater availability of the Central Valley Aquifer, California , 2009 .

[26]  Timothy H. Dixon,et al.  Space geodetic imaging of rapid ground subsidence in Mexico City , 2008 .

[27]  Yi Liu,et al.  Inverse procedure for calibrating parameters that control land subsidence caused by subsurface fluid withdrawal: 1. Methods , 2008 .

[28]  A. Ferretti,et al.  Permanent scatterer InSAR reveals seasonal and long‐term aquifer‐system response to groundwater pumping and artificial recharge , 2008 .

[29]  Jyr‐Ching Hu,et al.  Thirty-year land elevation change from subsidence to uplift following the termination of groundwater pumping and its geological implications in the Metropolitan Taipei Basin, Northern Taiwan , 2007 .

[30]  Peter Steigenberger,et al.  Generation of a consistent absolute phase-center correction model for GPS receiver and satellite antennas , 2007 .

[31]  Carl E. Anderson,et al.  Land Subsidence and Aquifer-System Compaction in the Tucson Active Management Area, South-Central Arizona, 1987-2005 , 2007 .

[32]  S. Megdal Arizona's recharge and recovery programs , 2006 .

[33]  Devin L. Galloway,et al.  Impacts of land subsidence caused by withdrawal of underground fluids in the United States , 2005 .

[34]  K. Jacobs,et al.  Managing for sustainability in an arid climate: lessons learned from 20 years of groundwater management in Arizona, USA , 2004 .

[35]  Michael T. Pavelko Estimates of hydraulic properties from a one-dimensional numerical model of vertical aquifer-system deformation, Lorenzi site, Las Vegas, Nevada , 2004 .

[36]  Howard A. Zebker,et al.  Inverse modeling of interbed storage parameters using land subsidence observations, Antelope Valley, California , 2003 .

[37]  W. Kustas,et al.  Comparison of energy balance modeling schemes using microwave‐derived soil moisture and radiometric surface temperature , 2003 .

[38]  S. A. Leake,et al.  MODFLOW-2000 Ground-Water Model?User Guide to the Subsidence and Aquifer-System Compaction (SUB) Package , 2003 .

[39]  M. Palecki,et al.  THE DROUGHT MONITOR , 2002 .

[40]  Barbara J. Morehouse,et al.  Sensitivity of urban water resources in Phoenix, Tucson, and Sierra Vista, Arizona, to severe drought , 2002 .

[41]  Howard A. Zebker,et al.  Seasonal subsidence and rebound in Las Vegas Valley, Nevada, observed by Synthetic Aperture Radar Interferometry , 2001 .

[42]  T. Burbey Storage Coefficient Revisited: Is Purely Vertical Strain a Good Assumption? , 2001, Ground water.

[43]  T. Burbey Stress‐Strain Analyses for Aquifer‐System Characterization , 2001 .

[44]  Michelle Sneed Hydraulic and mechanical properties affecting ground-water flow and aquifer-system compaction, San Joaquin Valley, California , 2001 .

[45]  D. Pool,et al.  Aquifer compaction and ground-water levels in south-central Arizona , 2000 .

[46]  H. Zebker,et al.  Sensing the ups and downs of Las Vegas: InSAR reveals structural control of land subsidence and aquifer-system deformation , 1999 .

[47]  Riccardo Lanari,et al.  Synthetic Aperture Radar Processing , 1999 .

[48]  C. Torrence,et al.  A Practical Guide to Wavelet Analysis. , 1998 .

[49]  J. Zumberge,et al.  Precise point positioning for the efficient and robust analysis of GPS data from large networks , 1997 .

[50]  T. W. Anderson Summary of the Southwest Alluvial Basins, Regional Aquifer-System Analysis, South-Central Arizona and Parts of Adjacent States , 1995 .

[51]  E. Frind,et al.  Hydraulic response of highly compressible aquitards during consolidation , 1991 .

[52]  Herbert F. Wang,et al.  Specific Storage as a Poroelastic Coefficient , 1990 .

[53]  R. Hanson AQUIFER-SYSTEM COMPACTION, TUCSON BASIN AND AVRA VALLEY, ARIZONA , 1988 .

[54]  S. R. Anderson,et al.  Land-subsidence measurements and aquifer-compaction monitoring in Tucson basin and Avra Valley, Arizona , 1988 .

[55]  S. R. Anderson Potential for aquifer compaction, land subsidence, and earth fissures in the Tucson basin, Pima County, Arizona , 1987 .

[56]  David E. Prudic,et al.  Ground-water flow in the Central Valley, California , 1985 .

[57]  J. F. Poland Guidebook to studies of land subsidence due to ground-water withdrawal , 1984 .

[58]  Edward M. Mikhail,et al.  Observations And Least Squares , 1983 .

[59]  Donald C. Helm,et al.  One‐dimensional simulation of aquifer system compaction near Pixley, California: 2. Stress‐Dependent Parameters , 1976 .

[60]  J. F. Poland,et al.  LAND SUBSIDENCE DUE TO WITHDRAWAL OF FLUIDS , 1969 .

[61]  C. E. Jacob On the flow of water in an elastic artesian aquifer , 1940 .

[62]  C. V. Theis The relation between the lowering of the Piezometric surface and the rate and duration of discharge of a well using ground‐water storage , 1935 .