Modeling Land Subsidence Using InSAR and Airborne Electromagnetic Data

American Geophysical Union (AGU) | March 13th, 2019

Land subsidence as a result of groundwater overpumping in the San Joaquin Valley, California, is associated with the loss of groundwater storage and aquifer contaminatio

Land subsidence as a result of groundwater overpumping in the San Joaquin Valley, California, is associated with the loss of groundwater storage and aquifer contamination. Although the physical processes governing land subsidence are well understood, building predictive models of subsidence is challenging because so much subsurface information is required to do so accurately. For the first time, we integrate airborne electromagnetic data, representing the subsurface, with subsidence data, mapped by interferometric synthetic aperture radar (InSAR), to model deformation. By combining both data sets, we are able to solve for hydrologic and geophysical properties of the subsurface to effectively model the complex spatiotemporal process of deformation. The resulting model reveals that roughly 3 m of subsidence has occurred at one location of our study area from 1990 to 2018. This model also allows us to predict subsidence more accurately under future hydrologic scenarios, which is needed to develop plans for sustainable groundwater use.

One‐dimensional simulation of aquifer system compaction near Pixley, California: 1. Constant parameters

American Geophysical Union (AGU) | June 12th, 1975

Aggregate one‐dimensional compaction (consolidation) of a series of aquitards in a compacting aquifer system has been simulated through use of a finite difference repre

Aggregate one‐dimensional compaction (consolidation) of a series of aquitards in a compacting aquifer system has been simulated through use of a finite difference representation of the vertical stress distribution within an idealized aquitard. Among the parameters affecting the simulated compaction are two storage coefficients (compressibility values), one for recoverable and the other for nonrecoverable compression. These two storage coefficients introduce a transient heterogeneity within an aquitard that is generally ignored by hydrologists. A computer program with two sets of constant coefficients calculates the daily deformation due to observed changes in applied stress near Pixley, California. Although water levels fluctuate annually, no long‐term water level decline occurred near Pixley between January 1, 1959, and February 4, 1971. During this period, 3.19 ft (0.972 m) of compaction was observed. The net difference between simulated and observed compaction on February 4, 1971, was 1.3% of the observed value. Maximum deviation occurred in mid‐1964 and equaled 7% of the observed compaction.

Potential Impacts of Future Geological Storage of CO2 on the Groundwater Resources in California’s Central Valley: Southern San Joaquin Basin Oil and Gas Production Analog for Geologic Carbon Storage

California Energy Commission (CEC) | May 1st, 2013

This study investigated the potential impacts of the injection and storage of large quantities of carbon dioxide on shallow water resources in California’s southern San

This study investigated the potential impacts of the injection and storage of large quantities of carbon dioxide on shallow water resources in California’s southern San Joaquin Basin, including the possibility that if such injection was properly planned it might help mitigate the land subsidence impacts related to over-pumping of groundwater. The project used simulation studies to investigate whether the basin-scale pressure changes and brine displacement caused by future carbon dioxide storage in the deep sediments could impact the groundwater-surface water systems in the area. Two simulation models were developed. The Deep Basin Reservoir Simulation Model was developed for accurately simulating the multiphase processes of carbon dioxide and brine flow occurring in deep sequestration reservoirs. The Shallow Water Resources Impact Simulation Model was developed for accurately predicting the impact of deep carbon dioxide sequestration on a stressed water system, including evaluation of carbon dioxide-related impacts in the context of the long-term groundwater and subsidence trends occurring in the area as a result of other water stresses such as pumping and irrigation. The Deep Basin Reservoir Simulation Model showed that the caprock and overlying geological units had sufficient barrier attributes to prevent pressure propagation from the injection layer to the shallow aquifers. The Shallow Water Resources Impact Simulation Model simulation results showed that deep carbon dioxide injection would create only minor head change and surface uplift in the shallow model layers. Overall the simulation results suggested that even large carbon dioxide storage projects should have minimal impact on shallow groundwater and surface water systems if thick layers of shales and sands separated the deep CO2 injection layers from shallow aquifers. Researchers believed this study justified devoting resources to evaluating the potential impact from localized pathways or from induced seismicity to more thoroughly investigate the feasibility of carbon dioxide storage in the southern San Joaquin Basin. 

Probabilistic Subsidence Forecast Model for the California Aqueduct Subsidence Program, San Joaquin Valley, California: Revision 1 Design Report

California Department of Water Resources (DWR) | October 4th, 2024

This report documents the development of a probabilistic subsidence forecast model for simulating a plausible range of future land-surface altitude conditions along the C

This report documents the development of a probabilistic subsidence forecast model for simulating a plausible range of future land-surface altitude conditions along the California Aqueduct (Aqueduct) and San Luis Canal in the San Joaquin Valley (SJV), with an emphasis on areas of localized subsidence (i.e., “subsidence bowls”). The model forecasts will be used to inform long-term planning and on-going analyses of potential investments needed to provide a suitable level of performance for the Aqueduct. The forecast model is primarily based on an empirical relationship between historical subsidence rates and annual water deliveries from the Central Valley Project (CVP) and State Water Project (SWP). A key assumption of the model is that the rate of groundwater overdraft, which contributes to loss of aquifer storage and permanent land subsidence, is correlated with CVP and SWP deliveries, and specifically with higher groundwater storage loss during severe drought years. Another key assumption is that the same rates of subsidence will continue even after groundwater levels decline and extraction may be coming from deeper parts of the aquifer system. If the geology and/or aquifer properties change dramatically in response, then these new conditions may have an unmodeled impact on future rates. 

Progress Report: Subsidence in California, March 2015 – September 2016

Jet Propulsion Laboratory (JPL) | February 7th, 2017

Subsidence caused by groundwater pumping in the Central Valley has been a problem for decades. Over the last few years, interferometric synthetic aperture radar (InSAR) h

Subsidence caused by groundwater pumping in the Central Valley has been a problem for decades. Over the last few years, interferometric synthetic aperture radar (InSAR) has been used from satellites and aircraft to produce maps of subsidence with sensitivity of fractions of an inch. For this study, we have obtained and analyzed data from the European Space Agency’s satellite-borne Sentinel-1A from the period March 2015 – September 2016 and the NASA airborne UAVSAR for the period March 2015 – June 2016, and produced maps of total subsidence from the two data sets. These data add to the earlier data processed from the Japanese PALSAR for 2006 – 2010, Canadian Radarsat-2 for the period May 2014 – January 2015, and UAVSAR for July 2013 - March 2015, for which subsidence measurements were reported previously (Farr et al., 2015). Here we also present results for the South-Central coast of California including Ventura, Oxnard, Santa Barbara and north to the San Joaquin Valley as well as the Santa Clara Valley from colleagues who have processed Sentinel-1A data covering March 2015 – March 2016. As multiple scenes were acquired during these periods, we can also produce time histories of subsidence at selected locations and transects showing how subsidence varies both spatially and temporally.

Pumping and Ground-Water Storage Depletion in Las Vegas Valley, Nevada, 1955-1974

U.S. Geological Survey (USGS) | July 1st, 1976

Quantification of record-breaking subsidence in California’s San Joaquin Valley

Nature Portfolio (Springer Nature) | November 19th, 2024

Rapid subsidence over oil fields measured by SAR interferometry: Geophysical Research Letters

American Geophysical Union (AGU) | September 1st, 1998

The Lost Hills and Belridge oilfields are in the San Joaquin Valley, California. The major oil reservoir is high porosity and low permeability diatomite. Extraction of la

The Lost Hills and Belridge oilfields are in the San Joaquin Valley, California. The major oil reservoir is high porosity and low permeability diatomite. Extraction of large volumes from shallow depths causes reduction in pore pressure and subsequent compaction, forming a surface sub- sidence bowl. We measure this subsidence from space using interferometric analysis of SAR (Synthetic Aperture Radar) data collected by the European Space Agency Remote Sensing Satellites (ERS-1 and ERS-2). Maximum subsidence rates are as high as 40 mm in 35 days or > 400 mm/yr, measured from interferograms with time separations ranging from one day to 26 months. The 8- and 26-month inter-ferograms contain areas where the subsidence gradient exceeds the measurement possible with ERS SAR, but shows increased detail in areas of less rapid subsidence. Synoptic mapping of subsidence distribution from satellite data powerfully complements ground-based techniques, permits measurements where access is difficult, and aids identification of underlying causes

Record of Decision Friant-Kern Canal Middle Reach Capacity Correction Project

U.S. Bureau of Reclamation (USBR) | November 4th, 2020

The United States Department of the Interior, Bureau of Reclamation (Reclamation), as the Federal lead agency under the National Environmental Policy Act (NEPA), and Fria

The United States Department of the Interior, Bureau of Reclamation (Reclamation), as the Federal lead agency under the National Environmental Policy Act (NEPA), and Friant Water Authority (FWA) as the State lead agency under the California Environmental Quality Act (CEQA), prepared the Friant-Kern Canal Middle Reach Capacity Correction Project (Project) Environmental Impact Statement/Environmental Impact Report (EIS/EIR) to assess impacts of the Project. The Project is being proposed to restore the capacity of a 33-mile reach of the Friant-Kern Canal (FKC) to its original design and constructed conveyance capacity. Click here for a directory with links to the federal environmental documents for the Friant-Kern Canal Middle Reach Capacity Correction Project.

Regional land subsidence accompanying groundwater extraction

Hydrogeology Journal (Springer) | August 26th, 2011

The extraction of groundwater can generate land subsidence by causing the compaction of susceptible aquifer systems, typically unconsolidated alluvial or basin-fill aquif

The extraction of groundwater can generate land subsidence by causing the compaction of susceptible aquifer systems, typically unconsolidated alluvial or basin-fill aquifer systems comprising aquifers and aquitards. Various ground-based and remotely sensed methods are used to measure and map subsidence. Many areas of subsidence caused by groundwater pumping have been identified and monitored, and corrective measures to slow or halt subsidence have been devised. Two principal means are used to mitigate subsidence caused by groundwater withdrawal—reduction of groundwater withdrawal, and artificial recharge. Analysis and simulation of aquifer-system compaction follow from the basic relations between head, stress, compressibility, and groundwater flow and are addressed primarily using two approaches—one based on conventional groundwater flow theory and one based on linear poroelasticity theory. Research and development to improve the assessment and analysis of aquifer-system compaction, the accompanying subsidence and potential ground ruptures are needed in the topic areas of the hydromechanical behavior of aquitards, the role of horizontal deformation, the application of differential synthetic aperture radar interferometry, and the regional-scale simulation of coupled groundwater flow and aquifer-system deformation to support resource management and hazard mitigation measures.