Detection and measurement of land subsidence and uplift using interferometric synthetic aperture radar, San Diego, California, USA, 2016–2018

International Association of Hydrological Sciences (IAHS) | April 22nd, 2020

Land subsidence associated with groundwater-level declines is stipulated as an “undesirable effect” in California’s Sustainable Groundwater Management Act (SGMA), a

Land subsidence associated with groundwater-level declines is stipulated as an “undesirable effect” in California’s Sustainable Groundwater Management Act (SGMA), and has been identified as a potential issue in San Diego, California, USA. The United States Geological Survey (USGS), the Sweetwater Authority, and the City of San Diego, undertook a cooperative study to better understand the hydromechanical response of the coastal aquifer system using Interferometric Synthetic Aperture Radar (InSAR) techniques. Three periods of interest were analyzed for this study that correspond to the periods before and after two substantial changes were made to the location and volume of pumpage: (1) April–August 2016 when groundwater levels and land surface elevation were relatively stable during normal pumping, (2) September 2016–May 2017 when groundwater levels recovered and the land surface uplifted during a period of substantially reduced pumping, (3) June 2017–October 2018 when groundwater levels declined and land subsidence occurred when pumpage resumed and expanded to new wells. Spatial and temporal characterization of the hydromechanical response to changes in pumpage is important for managing land subsidence. Further study using InSAR techniques, especially when combined with ground-based geodetic and monitoring-well networks, will provide water managers information to help effectively manage groundwater resources as stipulated in the SGMA.

Detection and Measurement of Land Subsidence Using Global Positioning System Surveying and Interferometric Synthetic Aperture Radar, Coachella Valley, California, 1996–2005

U.S. Geological Survey (USGS) | June 13th, 2013

Land subsidence associated with ground-water-level declines has been investigated by the U.S. Geological Survey in the Coachella Valley, California, since 1996. G

Land subsidence associated with ground-water-level declines has been investigated by the U.S. Geological Survey in the Coachella Valley, California, since 1996. Groundwater has been a major source of agricultural, municipal, and domestic supply in the valley since the early 1920s. Pumping of ground water resulted in water-level declines as large as 15 meters (50 feet) through the late 1940s. In 1949, the importation of Colorado River water to the southern Coachella Valley began, resulting in a reduction in ground-water pumping and a recovery of water levels during the 1950s through the 1970s. Since the late 1970s, demand for water in the valley has exceeded deliveries of imported surface water, resulting in increased pumping and associated ground-water-level declines and, consequently, an increase in the potential for land subsidence caused by aquifer-system compaction. Global Positioning System (GPS) surveying and interferometric synthetic aperture radar (InSAR) methods were used to determine the location, extent, and magnitude of the vertical land-surface changes in the southern Coachella Valley. GPS measurements made at 13 geodetic monuments in 1996 and in 2005 in the southern Coachella Valley indicate that the elevation of the land surface had a net decline of 124 to 9 mm ±54 mm (0.41 to 0.03 ft ±0.18 ft) during the 9-year period. Changes at 9 of the 13 monuments exceeded the maximum expected uncertainty of ±54 mm (±0.18 ft) at the 95-percent confidence level, indicating that subsidence occurred at these monuments between June 1996 and August 2005. GPS measurements made at 20 geodetic monuments in 2000 and in 2005 indicate that the elevation of the land surface changed –192 to +51 mm ±36 mm (–0.63 to +0.17 ft ±0.12 ft) during the 5-year period. Changes at 6 of the 20 monuments exceeded the maximum expected uncertainty of ±36 mm (±0.12 ft) at the 95-percent confidence level—subsidence occurred at five monuments and uplift occurred at one monument between August 2000 and August 2005. GPS measurements at two of the five subsiding monuments for which subsidence rates could be compared indicate that subsidence rates decreased during this period compared with subsidence rates before 2000. InSAR measurements made between May 7, 2003, and September 25, 2005, indicate that land subsidence, ranging from about 75 to 180 millimeters (0.25 to 0.59 foot), occurred in three areas of the Coachella Valley: near Palm Desert, Indian Wells, and La Quinta; the equivalent subsidence rates range from about 3 to more than 6 mm/month (0.01 to 0.02 ft/month). The subsiding areas near Palm Desert, Indian Wells, and La Quinta were previously identified using InSAR measurements for 1996–2000, which indicated that about 35 to 150 mm (0.11 to 0.49 ft) of subsidence occurred during the four-year period; the equivalent subsidence rates range from about 1 to 3 mm/month (0.003 to 0.01 ft/month). Comparison of the InSAR results indicates that subsidence rates have increased 2 to 4 times since 2000 in these three areas. Water-level measurements made at wells near the subsiding monuments and in the three subsiding areas generally indicated that the water levels fluctuated seasonally and declined annually between 1996 and 2005; some water levels in 2005 were at the lowest levels in their recorded histories. The coincident areas of subsidence and declining water levels suggest that aquifer-system compaction may be causing subsidence. If the stresses imposed by the historically lowest water levels exceeded the preconsolidation stress, the aquifer-system compaction and associated land subsidence may be permanent. Although the localized character of the subsidence signals is typical of the type of subsidence characteristically caused by localized ground-water pumping, the subsidence may also be related to tectonic activity in the valley.

Detection and Measurement of Land Subsidence Using Interferometric Synthetic Aperture Radar and Global Positioning System, San Bernardino County, Mojave Desert, California

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

Detection and Measurement of Land-Surface Deformation, Pajaro Valley, Santa Cruz and Monterey Counties, California, 2015–18

U.S. Geological Survey (USGS) | October 22nd, 2021

Land-surface deformation (subsidence) caused by groundwater withdrawal is identified as an undesirable result in the Pajaro Valley Water Management Agency’s Basin Manag

Land-surface deformation (subsidence) caused by groundwater withdrawal is identified as an undesirable result in the Pajaro Valley Water Management Agency’s Basin Management Plan and California’s Sustainable Groundwater Management Act. In Pajaro Valley, groundwater provides nearly 90 percent of the total water supply. To aid the development of sustainable groundwater management criteria, the U.S. Geological Survey, in cooperation with the Pajaro Valley Water Management Agency, performed an analysis of land-surface deformation (subsidence and uplift) in Pajaro Valley for 2015–18, using Interferometric Synthetic Aperture Radar and continuous Global Positioning System methods. Land-surface deformation results were then compared with subsurface geology and groundwater altitudes to better understand the hydromechanical response of the coastal aquifer system. The results indicate the land surface is generally stable with only small magnitudes (less than 1 inch) of seasonal land-surface deformation (subsidence in the summer and uplift in the winter) during 2015–18. During this time, the largest magnitude of land-surface deformation was less than 2 inches of subsidence and was localized in one area just north of the city limits of Watsonville, California. Groundwater altitudes during 2015–18 demonstrated seasonal variability and annual to multi-annual increases after reaching historical lows by the mid-1990s. The small magnitudes of land-surface deformation coupled with groundwater-altitude increases in most areas indicate that the subsidence likely is largely elastic and recoverable. The Corralitos-Pajaro Valley groundwater basin contains fine-grained (clay) sediments that have the potential for permanent aquifer-system compaction and resultant land subsidence. However, groundwater altitudes throughout the Pajaro Valley have increased above historical lows, and observed increases in groundwater altitudes coincided with changes in groundwater management activities. Observed relations between groundwater management activities and groundwater altitudes indicate that management of groundwater supplies could minimize the potential for permanent land-surface deformation in Pajaro Valley.

Detection of aquifer system compaction and land subsidence using interferometric synthetic aperture radar, Antelope Valley, Mojave Desert, California

U.S. Geological Survey (USGS) | April 15th, 1998

Interferometric synthetic aperture radar (InSAR) has great potential to detect and quantify land subsidence caused by aquifer system compaction. InSAR maps with h

Interferometric synthetic aperture radar (InSAR) has great potential to detect and quantify land subsidence caused by aquifer system compaction. InSAR maps with high spatial detail and resolution of range displacement (10 mm in change of land surface elevation) were developed for a groundwater basin (103 km2) in Antelope Valley, California, using radar data collected from the ERS-1 satellite. These data allow comprehensive comparison between recent (1993–1995) subsidence patterns and those detected historically (1926–1992) by more traditional methods. The changed subsidence patterns are generally compatible with recent shifts in land and water use. The InSAR-detected patterns are generally consistent with predictions based on a coupled model of groundwater flow and aquifer system compaction. The minor inconsistencies may reflect our imperfect knowledge of the distribution and properties of compressible sediments. When used in conjunction with coincident measurements of groundwater levels and other geologic information, InSAR data may be useful for constraining parameter estimates in simulations of aquifer system compaction.

Determination of Study Reporting Limits for Pesticide Constituent Data for the California Groundwater Ambient Monitoring and Assessment Program Priority Basin Project, 2004–2018—Part 1: National Water Quality Laboratory Schedules 2003, 2032, or 2033, and 2060

U.S. Geological Survey (USGS) | December 17th, 2019

The California Groundwater Ambient Monitoring and Assessment Program Priority Basin Project (GAMA-PBP) is a long-term cooperative project designed to assess the quality o

The California Groundwater Ambient Monitoring and Assessment Program Priority Basin Project (GAMA-PBP) is a long-term cooperative project designed to assess the quality of groundwater resources used for public and domestic drinking water supplies in the State of California, to monitor and evaluate changes to that quality, to investigate the human and natural factors controlling water quality, and to improve the availability of comprehensive groundwater quality data and information. Between May 18, 2004, and May 3, 2018, the GAMA-PBP collected 3001 groundwater samples for analysis of pesticide constituents by the U.S. Geological Survey (USGS) National Water Quality Laboratory (NWQL)(note that ‘pesticide constituents’ includes parent compounds and degradates). Of these samples, 2994 were analyzed for pesticide constituents on schedules 2003, 2032, or 2033 (65 to 84 constituents), and 840 were analyzed for pesticide constituents on schedule 2060 (58 constituents). The original dataset reported by the NWQL to the USGS National Water Information System (NWIS) database contained a total of 2,688 detections of 78 pesticide constituents and 253,825 non-detections. In this original dataset, 33 percent of the 3,001 samples analyzed had reported detections of one or more pesticide constituents. This report describes the GAMA-PBP data-quality objectives for pesticide data, the procedures used to establish study reporting limits, and use of those reporting limits to censor the data from the NWQL so that the final data published by the GAMA-PBP meet these data-quality objectives. The final GAMA-PBP dataset for samples collected from May 2004 to May 2018, after censoring, had a total of 1,632 detections of 37 pesticide constituents. In the final GAMA-PBP dataset, 25 percent of the 3,001 samples analyzed had detections of one or more pesticide constituents. The presence of pesticides in groundwater is commonly evaluated by calculating detection frequencies. Detection frequencies for pesticides are sensitive to detection limits and method performance for concentrations near those limits; therefore, the two primary data quality issues addressed in the GAMA-PBP data-quality objectives for pesticides are (1) establishing criteria for classifying data from the laboratory as detections or non-detections for the purpose of data reporting by the project and (2) accounting for changes in analytical methods or method performance over time. The GAMA-PBP addresses these issues by developing study reporting limits that are used as the boundary between detections and non-detections for the reporting of GAMA-PBP results. These reporting limits are defined from method detection limits (MDLs) provided by the NWQL, unless examination of results from laboratory set blanks (LSBs) and GAMA-PBP field blanks indicates that a higher concentration censoring limit is warranted. The GAMA-PBP selected the MDL as the primary choice for defining study reporting limits for consistency with U.S. Environmental Protection Agency (EPA) guidelines for reporting detections of pesticides and other organic constituents. A five-step procedure is used to develop study reporting limits and censor the GAMA-PBP dataset accordingly. The effect of the censoring at each step is described to provide information about the relative effect of each step on the overall censoring of the dataset. Steps 1 and 2 can be implemented at the time the data are received, whereas steps 3−5 require information accumulated over an extended period. • Step 1: Reject results that were most likely the result of specific contamination instances attributable to unusual field or laboratory conditions during sample collection or processing. Two such instances were identified, leading to rejection of 25 detections, which were assigned a data-quality indicator code of “Q” for “reviewed and rejected” in the NWIS database. • Step 2: Use the NWQL MDLs in effect at the time each sample was analyzed as the reporting limit. A total of 506 detections were censored on this basis. • Step 3: Use the maximum MDL established by the NWQL during July 2004–August 2018 (MDLmax) as the reporting limit. The rationale for using the MDLmax as the reporting limit is based primarily on the observation that the concentrations of MDLs generally increased over time. A total of 438 detections were censored on this basis. • Step 4: Use the LSBs to identify periods of greater potential laboratory contamination bias and define raised reporting limits to be used during those periods. These periods were defined by using a moving average detection frequency approach. For consistency with the NWQL procedures for defining raised reporting limits on the basis of detections in LSBs, the raised reporting limits were defined as equal to three times the highest concentration measured in an LSB during the period. A total of 25 detections in groundwater samples analyzed during periods of increased laboratory contamination bias were censored. • Step 5: Use the LSBs and field blanks to identify potential contamination bias from field or laboratory processes outside of the time periods identified in step 4. The NWQL protocols were used to define the MDLs from blanks analyzed outside of the periods identified in step 4. If an MDL defined from blanks was greater than the MDLmax, the MDL defined from blanks was used to censor the data. One constituent had a study reporting limit defined on this basis, and a total of 62 detections in groundwater samples were censored. As of 2019, the USGS NWIS database does not have the capability to store both the original value reported by the NWQL and the final value published by the GAMA-PBP that reflects application of the quality-control censoring described in this report. In the interim, while this capability is developed, the 1,031 results censored in steps 2−5 are blocked from public release in NWIS, and the GAMA-PBP has published the original and final values in a USGS data release accompanying this report. The entire GAMA-PBP final dataset for pesticide constituents on schedules 2003, 2032, or 2033, or on schedule 2060 is publicly available in that USGS data release, through the USGS GAMA-PBP public web portal, and through the California State Water Resources Control Board GAMA public groundwater information system.

Development and Application of a 1D Compaction Model to Understand 65 Years of Subsidence in the San Joaquin Valley

American Geophysical Union (AGU) | June 2nd, 2022

High rates of land subsidence, caused by groundwater overdraft, are resulting in millions of dollars of infrastructure damage in California's San Joaquin Valley (SJV). In

High rates of land subsidence, caused by groundwater overdraft, are resulting in millions of dollars of infrastructure damage in California's San Joaquin Valley (SJV). In recent years, the use of interferometric synthetic aperture radar (InSAR) has enabled us to substantially improve our understanding of this subsidence. However, only very occasionally have the InSAR data been integrated with a physical model of subsurface compaction. Here, we have used InSAR and other data to parameterize and calibrate a 1D compaction model. We applied our model to a study area in the SJV where we had access to the necessary information on hydraulic head to develop model inputs. Our model simulated subsidence in the three aquifer system layers over the period 1952–2017, and is the first 1D compaction model in the SJV to simulate multiple aquifer system layers from the 1950s to 2017. The results from our model suggest that previous studies have significantly underestimated the time constants governing the slow, residual compaction of subsurface clays. We suggest that residual compaction of clays is a process that continues for decades-to-centuries, indicating that to significantly reduce subsidence requires some recovery of head, not just a stabilization. We also show how compaction in the lower, confined aquifer has accounted for over 90% of subsidence in the past 20 years. Although our study area is small, our findings are likely representative of the subsiding regions of the SJV, and our methodology can be applied to unconsolidated aquifer systems exhibiting subsidence worldwide.

Development and Calibration of a Two-Dimensional Digital Model for the Analysis of the Ground-Water Flow System in the San Antonio Creek Valley, Santa Barbara County, California

U.S. Geological Survey (USGS) | August 1st, 1985

Development of a Groundwater-Simulation Model in the Los Angeles Coastal Plain, Los Angeles County, California

U.S. Geological Survey (USGS) | September 28th, 2021

The Los Angeles Coastal Plain (LACP) covers about 580 square miles and is the largest coastal plain of semiarid southern California. The LACP is heavily developed with mo

The Los Angeles Coastal Plain (LACP) covers about 580 square miles and is the largest coastal plain of semiarid southern California. The LACP is heavily developed with mostly residential, commercial, and industrial land uses that rely heavily on groundwater for water supply. In 2010, the LACP was home to about 14 percent of California’s population, or about 5.4 million residents. The LACP is also a major commercial and industrial hub with industries including manufacturing, aerospace, entertainment, and tourism. There has been a heavy reliance on groundwater from the LACP for many years. An average of 305,000 acre-feet per year (acre-ft/yr) of groundwater was used annually from the LACP from 1971 to 2015. The need to replenish the groundwater basins within the LACP was recognized as far back as the 1930s, when spreading grounds were first used to replenish groundwater basins and store water underground during times of water surplus to meet demands in times of shortage. Seawater intrusion resulting from freshwater pumping was first observed in the 1940s. As a result, injection of imported water through wells at what is now the West Coast Basin Barrier Project began on an experimental basis in 1951. Managed aquifer recharge from the spreading grounds and barrier wells is now a substantial component of the LACP’s groundwater supply. The average annual recharge from water spreading from 1971 to 2015 was about 120,000 acre-ft/yr, and the average annual injection into the barrier wells was about 33,000 acre-ft/yr. Other inflows include areal recharge, underflow from San Gabriel and San Fernando Valleys, and onshore flow from the ocean. The average annual recharge from these sources was 100,000 acre-feet (acre-ft) from 1971 to 2015. Additionally, cross-boundary flow from Orange County into the western Orange County subareas of the LACP was simulated as 48,000 acre-ft from 1971 to 2015. This study, conducted in cooperation with the Water Replenishment District of Southern California (WRD), involved an assessment of the historical and present status of groundwater resources in the LACP and the development of tools to better understand the groundwater system. These efforts were built upon results from previous studies and incorporate new information and developments in modeling capabilities to provide a more detailed analysis of the aquifer systems. This study includes a comprehensive compilation of geologic and hydrologic data (Chapter A), development of a chronostratigraphic model that provides a detailed description of the LACP aquifer systems (Chapter B), characterization of the groundwater hydrology of the LACP, including a down-hole analysis of grain size using lithologic and geophysical logs (Chapter C), and development and application of the Los Angeles Coastal Plain Groundwater-flow Model (LACPGM) to simulate past groundwater conditions, estimate groundwater-budget components and flow paths, and approximate future groundwater conditions under different scenarios (Chapter D).

DGMETA (Version 1): Dissolved Gas Modeling and Environmental Tracer Analysis Computer Program

U.S. Geological Survey (USGS) | January 8th, 2021

DGMETA (Dissolved Gas Modeling and Environmental Tracer Analysis) is a Microsoft Excel-based computer program that is used for modeling air-water equilibrium conditions f

DGMETA (Dissolved Gas Modeling and Environmental Tracer Analysis) is a Microsoft Excel-based computer program that is used for modeling air-water equilibrium conditions from measurements of dissolved gases and for computing concentrations of environmental tracers that rely on air-water equilibrium model results. DGMETA can solve for the temperature, salinity, excess air, fractionation of gases, or pressure/elevation of water when it is equilibrated with the atmosphere. Models are calibrated inversely using one or more measurements of dissolved gases such as helium, neon, argon, krypton, xenon, and nitrogen. Excess nitrogen gas, originating from denitrification or other sources, also can be included as a fitted parameter or as a separate calculation from the dissolved gas modeling results. DGMETA uses the air-water equilibrium models to separate measured concentrations of gases and isotopes of gases into components that are used for tracing water in the environment. DGMETA calculates atmospheric dry-air mole fractions (mixing ratios) for transient atmospheric gas tracers such as chlorofluorocarbons, sulfur hexafluoride, and bromotrifluoromethane (Halon-1301); and concentrations of tritiogenic helium-3 and radiogenic helium-4, which accumulate from the decay of tritium in water and the decay of uranium and thorium in rocks, respectively. Sample data can be graphed to identify applicable models of excess air, samples that contain excess nitrogen gas, or samples that have partially degassed, for example. Monte Carlo analysis of errors associated with dissolved gas equilibrium model results can be carried through computations of environmental tracer concentrations to provide robust estimates of error. In addition, graphical routines for separating helium sources using helium isotopes are included to refine estimates of tritiogenic helium-3 when terrigenic helium from mantle or crustal sources is present in samples. Environmental tracer concentrations and their errors computed from DGMETA can be used with other programs, such as TracerLPM (Jurgens and others, 2012), to determine groundwater ages and biogeochemical reaction rates. DGMETA also produces output files in a format that meets the U.S. Geological Survey open data requirements for documentation of model inputs and outputs. DGMETA is a versatile and adaptable program that allows users to add solubility data for new gases, modify the existing set of gas solubility data, modify the default set of gases used for modeling, choose calculations based on real (non-ideal) gas behavior, and select various concentration units for data entry and results to match laboratory reports and study objectives. DGMETA comes with a set of gases widely used in hydrology and oceanography and many gases include multiple solubilities from previous work. Seventeen dissolved gases are included in the default version of the program: noble gases (helium, neon, argon, krypton, and xenon), reactive gases (nitrogen, oxygen, methane, carbon dioxide, carbon monoxide, hydrogen, and nitrous oxide), and environmental tracers (chlorofluorocarbon-11, chlorofluorocarbon-12, chlorofluorocarbon-113, sulfur hexafluoride, and Halon-1301).