Search USGSSearch

Geology topics

Justin T. Brandt

Publications and source records attributed to Justin T. Brandt.

At least 19 recordsLinked to original sources

Groundwater sustainability and land subsidence in California’s Central Valley

The Central Valley of California is one of the most prolific agricultural regions in the world. Agriculture is reliant on the conjunctive use of surface-water and groundwater. The lack of available surface-water and land-use changes have led to pumping-induced groundwater-level and storage declines, land subsidence, changes to streamflow and the environment, and the degradation of water quality. As a result, in part, the Sustainable Groundwater Management Act (SGMA) was developed. An examination of the components of SGMA and contextualizing regional model applications within the SGMA framework was undertaken to better understand and quantify many of the components of SGMA. Specifically, the U.S. Geological Survey (USGS) updated the Central Valley Hydrologic Model (CVHM) to assess hydrologic system responses to climatic variation, surface-water availability, land-use changes, and groundwater pumping. MODFLOW-OWHM has been enhanced to simulate the timing of land subsidence and attribute its inelastic and elastic portions. In addition to extending CVHM through 2019, the new version, CVHM2, includes several enhancements as follows: managed aquifer recharge (MAR), pumping with multi-aquifer wells, inflows from ungauged watersheds, and more detailed water-balance subregions, streamflow network, diversions, tile drains, land use, aquifer properties, and groundwater level and land subsidence observations. Combined with historical approximations, CVHM2 estimates approximately 158 km 3 of storage loss in the Central Valley from pre-development to 2019. About 15% of the total storage loss is permanent loss of storage from subsidence that has caused damage to infrastructure. Climate extremes will likely complicate the efforts of water managers to store more water in the ground. CVHM2 can provide data in the form of aggregated input datasets, simulate climatic variations and changes, land-use changes or water management scenarios, and resulting changes in groundwater levels, storage, and land subsidence to assist decision-makers in the conjunctive management of water supplies.

California

Hydrogeology, land-surface subsidence, and documentation of the Gulf Coast Land Subsidence and Groundwater-Flow (GULF) model, southeast Texas, 1897–2018

Executive Summary As a part of the Texas Water Development Board groundwater availability modeling program, the U.S. Geological Survey developed the Gulf Coast Land Subsidence and Groundwater-Flow model (hereinafter, the “GULF model”) and ensemble to simulate groundwater flow and land-surface subsidence in the northern part of the Gulf Coast aquifer system (the study area) in Texas from predevelopment (1897) through 2018. Since the publication of a previous groundwater model for the greater Houston area in 2012, there have been changes to the distribution of groundwater withdrawals and advances in modeling tools. To reflect these changes and to simulate more recent conditions, the GULF model was developed in cooperation with the Harris-Galveston and Fort Bend Subsidence Districts to provide an updated Groundwater Availability Model. Since the early 1900s, most of the groundwater withdrawals in the study area have been from three of the hydrogeologic units that compose the Gulf Coast aquifer system—the Chicot, Evangeline, and Jasper aquifers and, more recently, from the Catahoula confining unit. Withdrawals from these hydrogeologic units are used for municipal supply, commercial and industrial use, and irrigation purposes. Withdrawals of large quantities of groundwater in the greater Houston area have caused widespread groundwater-level declines in the Chicot, Evangeline, and Jasper aquifers of more than 300 feet (ft). Early development of the aquifer system, which began before 1900, resulted in nearly 50 percent of the eventual historical groundwater-level minimums having been reached as early as 1946 in some areas. These groundwater-level declines led to more than 9 ft of land-surface subsidence—historically in central and southeastern Harris County and Galveston County, but more recently in northern, northwestern, and western Harris County, Montgomery County, and northern Fort Bend County—from depressurization and compaction of clay and silt layers interbedded in the aquifer sediments. In a generalized conceptual model of the Gulf Coast aquifer system, water enters the groundwater system in topographically high outcrops of the hydrogeologic units in the northwestern part of the aquifer system. Groundwater that does not discharge to streams flows to intermediate and deep zones of the aquifer system southeastward of the outcrop areas where it is discharged by wells and by upward leakage in topographically low areas near the coast. The uppermost parts of the aquifer system, which include outcrop areas, are under water-table (unconfined) conditions where the groundwater is not confined under pressure. As depth increases in the aquifer system and interbedded clay and silt layers accumulate, water-table conditions evolve into confined conditions where the groundwater is under pressure. Groundwater flow and land-surface subsidence in the GULF model and ensemble were simulated by using MODFLOW 6 with the Skeletal Storage, Compaction, and Subsidence package. The model consists of six layers, one for each of the five hydrogeologic units in the northern part of the Gulf Coast aquifer system and a surficial top layer that includes part of each hydrogeologic unit. Transient groundwater flow was simulated during 1897–2018 by using a combination of multiyear, annual, and monthly stress periods. An initial steady-state stress period was configured to represent predevelopment mean annual inflows and outflows. The subsidence package used in the GULF model and ensemble uses a head-based subsidence formulation that simulates the delayed drainage response from clay and silt sediment to changes in groundwater levels. The GULF model and ensemble were history matched to groundwater-level observations at selected wells, land-surface subsidence at benchmarks, aquifer compaction at borehole extensometers, and vertical displacement from Global Positioning System stations. A Bayesian framework was used to represent uncertainty in modeled parameters and simulated outputs of interest. History matching and uncertainty quantification were performed by using a Monte Carlo approach enabled through iterative ensemble smoother software to produce an ensemble of models fit to historical data. The iterative ensemble smoother substantially reduced the computational demand of parameter estimation by approximating the first-order relation between model inputs and outputs, thereby allowing 183,207 adjustable parameters to be used for history matching at a relatively low computational and time cost. The history-matched parameter values are within the ranges of previously published values and agree with the current understanding of the spatial and temporal patterns of parameter uncertainty for the Gulf Coast aquifer system. A good agreement between the observed (or estimated) and simulated groundwater levels, land-surface subsidence, compaction, and vertical displacement was obtained across the modeled area based on qualitative and quantitative comparisons. Ensemble mean annual groundwater-flow rates to the Chicot, Evangeline, Jasper aquifers and Catahoula confining unit were 0.0–0.49 inch (in.), 0.09–0.33 in., 0.01–0.07 in., and 0.01–0.05 in., respectively. GULF model mean annual groundwater-flow rates to the Chicot, Evangeline, and Jasper aquifers and Catahoula confining unit were 0.31 in., 0.19 in., 0.03 in., and 0.03 in., respectively. The GULF-model-simulated recharge to the outcrop area was the largest inflow (75 percent), and recharge to other areas was 25 percent of the model inflow. The simulated outflows included (1) net surface-water/groundwater exchange with study area streams (50 percent), (2) groundwater use (49 percent), and (3) net surface-water/groundwater exchange with the Gulf of Mexico (1 percent). The sum of the simulated values of the outflows (1,041,973 acre-feet per year [acre-ft/yr]) and the elastic expansion of the fine-grained sediment and numerical solver error (339 acre-ft/yr) minus the inflows (654,172 acre-ft/yr) represents the reduction of storage from the Gulf Coast aquifer system (388,140 acre-ft/yr). Most of the storage depletion is caused by the long-term groundwater-level declines that have resulted primarily in inelastic compaction. The GULF model was used to estimate Jasper aquifer compaction at selected benchmarks in Montgomery County and northern Harris County, which are the primary locations of Jasper aquifer groundwater use. Simulated Jasper aquifer compaction in northern Harris County was between 0.2 and 0.5 ft, or between about 5 and 16 percent of simulated subsidence at the benchmark locations. Simulated Jasper aquifer compaction in Montgomery County was between 0.8 and 1.2 ft, or between about 33 and 57 percent of simulated subsidence at the benchmark locations.

Texas

Mapping structural control through analysis of land-surface deformation for the Rialto-Colton groundwater subbasin, San Bernardino County, California, 1992–2010

The locations of many faults in and near the Rialto-Colton groundwater subbasin are not precisely known because the spatial density of existing lithologic and hydrologic data used to infer the locations of faults can be sparse. The U.S. Geological Survey, in cooperation with the San Bernardino Valley Municipal Water District, analyzed structural control of groundwater flow in and near the Rialto-Colton groundwater subbasin using Interferometric Synthetic Aperture Radar (InSAR) methods. Faults commonly are barriers to groundwater flow, and the high spatial resolution of InSAR imagery can be used to infer the locations of buried faults where groundwater pumping occurs. InSAR results have revealed three areas in and near the Rialto-Colton groundwater subbasin where buried faults are interpreted as groundwater-flow barriers: the northwestern area about 3 miles northwest of the City of Rialto, the San Jacinto fault area west of the City of San Bernardino, and the southeastern area about 2 miles southeast of the City of Colton. The InSAR results were combined with knowledge gained from previous studies to better define the location and extent of faults acting as groundwater-flow barriers. New data about faults acting as groundwater-flow barriers can be incorporated into future conceptual and hydrologic models of the Rialto-Colton groundwater subbasin and provide water managers information to help effectively manage groundwater resources.

California

Detection and measurement of land-surface deformation, Pajaro Valley, Santa Cruz and Monterey counties, California, 2015–18

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.

California

Documentation of methods and inventory of irrigation information collected for the 2015 U.S. Geological Survey estimated use of water in the United States

The U.S. Geological Survey (USGS) National Water-Use Science Project strives to report water-use estimates using the best available information for the period of the estimates. The information available on water used for irrigation activities varies from State to State and in some areas from county to county within a State, which results in many information sources and methods being used to estimate water withdrawals and consumption for the Nation. The variety of estimation methods makes it difficult to compare information across States and makes it difficult to understand how different methods or data sources bias irrigation water-use estimates and trends over time. The sources of information and methods used by USGS Water Science Centers to estimate irrigation water use (the number of irrigated acres by irrigation system type, withdrawal values by water source type, and consumed-water values) for 2015 are compiled and described herein to assist with interpreting the water-use estimates. State-level summaries of information sources and methods are compiled in appendix 1, and a dataset of calendar-year, county-level estimates of actual evapotranspiration for the conterminous United States and Hawaii is provided in an associated USGS data release.

Scientific Investigations Report

Detection and measurement of land subsidence and uplift using Global Positioning System surveys and interferometric synthetic aperture radar, Coachella Valley, California, 2010–17

Groundwater has been a major source of agricultural, recreational, municipal, and domestic supply in the Coachella Valley of California since the early 1920s. Pumping of groundwater resulted in groundwater-level declines as large as 50 feet (ft) or 15 meters (m) by the late 1940s. Because of concerns that the declines could cause land subsidence, the Coachella Valley Water District (CVWD) and the U.S. Geological Survey (USGS) have cooperatively investigated subsidence in the Coachella Valley since 1996. Importation of Colorado River water to the southern Coachella Valley began in 1949, resulting in a reduction in groundwater pumping and a recovery of groundwater levels during the 1950s through the 1970s. Since the late 1970s, the demand for water in the valley increased to the point that groundwater levels again declined in response to increased pumping and, consequently, increased the potential for land subsidence caused by aquifer-system compaction. Several management actions to increase recharge or to reduce reliance on groundwater have been implemented since as early as 1973 to address overdraft in the Coachella Valley. The implementation of three particular projects has markedly improved groundwater conditions in some of the historically most overdrafted areas of the valley: (1) groundwater substitution with surface-water imports since 2006 using Colorado River water through the Mid-Valley Pipeline project, which was expanded through 2017; (2) budget-based, tiered rates since 2009; and (3) managed aquifer recharge at the Thomas E. Levy Groundwater Replenishment Facility since 2009. 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 Coachella Valley during 2010–17, updating 1993–2010 information presented in previous USGS reports. The GPS measurements taken at 24 geodetic monuments in August 2010 and September 2015 indicated that the land-surface elevation was stable at 17 monuments but changed at seven monuments during the 5-year period. Subsidence ranged from 0.17 to 0.43 ±0.09 ft (52 to 132 ±28 millimeters, or mm) at three monuments, and uplift ranged from 0.11 to 0.18 ±0.09 ft (33 to 54 ±28 mm) at four monuments between 2010 and 2015. At two of the monuments that subsided, the subsidence rates decreased between 2010 and 2015 from those computed between 2005 and 2010. Data prior to 2010 were not available for the third monument that subsided; thus, the 2010–15 subsidence rate could not be compared to an earlier period. At three of the monuments that uplifted between 2010 and 2015, data collected in 2005 and 2010 indicated stability. Data prior to 2010 were not available for the fourth monument that uplifted; thus, the 2010–15 uplift rate could not be compared to an earlier period. InSAR analyses for December 28, 2014–June 27, 2017, indicated that the land surface uplifted as much as about 0.20 ft (60 mm) near the Whitewater River Groundwater Replenishment Facility in the northern Coachella Valley and subsided as much as about 0.26 ft (80 mm) in the La Quinta area and less in Palm Desert, Indian Wells, and other localized areas in the southern Coachella Valley. These areas were identified as subsidence areas in previous reports covering periods during 1993–2010. The comparison of 2014–17 subsidence rates with those derived for 1995–2010 generally indicated a substantial slowing of subsidence, however. Analyses of deformation in the northern Coachella Valley were not included in the previous reports, so a comparison to deformation during the earlier period could not be made. Water levels in wells near the subsiding geodetic monuments, in and near the three subsiding areas shown by InSAR, and throughout the valley generally indicated seasonal fluctuations and longer-term stability or rising groundwater levels since about 2010. These results mark a reversal in trends of groundwater-level declines during the preceding decades. This trend reversal provides new insights into aquifer-system mechanics. Although many areas have stopped subsiding, and a few have even uplifted, the few areas that did subside during 2010–17—albeit at a slower rate—indicate a mixed aquifer-system response. Subsidence when groundwater levels are stable or recovering indicates that residual compaction may have occurred. At the same time, coarse-grained materials and thin aquitards may have expanded as groundwater levels recovered. The continued valley-wide stabilization and recovery of groundwater levels since 2010 likely is a result of various projects designed to increase recharge or to reduce reliance on groundwater.

California

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

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.

California

Mitigating land subsidence in the Coachella Valley, California, USA: An emerging success story

Groundwater has been a major source of agricultural, municipal, and domestic water supply since the early 1920s in the Coachella Valley, California, USA. Land subsidence, resulting from aquifer-system compaction and groundwater-level declines, has been a concern of the Coachella Valley Water District (CVWD) since the mid-1990s. As a result, the CVWD has implemented several projects to address groundwater overdraft that fall under three categories – groundwater substitution, conservation, and managed aquifer-recharge (MAR). The implementation of three projects in particular – replacing groundwater extraction with surface water from the Colorado River and recycled water (Mid-Valley Pipeline project), reducing water usage by tiered-rate costs, and increasing groundwater recharge at the Thomas E. Levy Groundwater Replenishment Facility – are potentially linked to markedly improved groundwater levels and subsidence conditions, including in some of the historically most overdrafted areas in the southern Coachella Valley. Groundwater-level and subsidence monitoring have tracked the effect these projects have had on the aquifer system. Prior to about 2010, water levels persistently declined, and some had reached historically low levels by 2010. Since about 2010, however, groundwater levels have stabilized or partially recovered, and subsidence has stopped or slowed substantially almost everywhere it previously had been observed; uplift was observed in some areas. Furthermore, results of Interferometric Synthetic Aperture Radar analyses for 1995–2017 indicate that as much as about 0.6 m of subsidence occurred; nearly all of which occurred prior to 2010. Continued monitoring of water levels and subsidence is necessary to inform the CVWD about future mitigation measures. The water management strategies implemented by the CVWD can inform managers of other overdrafted and subsidence-prone basins as they seek solutions to reduce overdraft and subsidence.

California

Evaluation of land subsidence and ground failures at Bicycle Basin, Fort Irwin National Training Center, California, 1992–2017

Groundwater has been pumped in the Bicycle Basin at Fort Irwin National Training Center since the 1960s, and the amount pumped has generally increased since the 1990s. After a large crack (approximately 0.5-kilometer long) formed at the surface of Bicycle Lake playa during 2005–06 in the area used as an aircraft runway, a monitoring study was initiated by the U.S. Geological Survey, in cooperation with the U.S. Army Fort Irwin National Training Center, to help determine the cause. The extent and effect of groundwaterlevel declines and land-surface deformation in Bicycle Basin were evaluated using a number of approaches, including water-level measurements from December 2007 to June 2017, land surveys across the playa area, interferometric synthetic aperture radar (InSAR) analyses, geophysical surveys of the playa area, and numerical experiments to test hypotheses about soil mechanical processes. A specific objective of this study was to evaluate the recent development of ground failures in the form of large, surface cracks that pose a hazard to aircraft operations on the Bicycle Lake playa airstrip. Another large crack and set of interconnected cracks formed at the surface of Bicycle Lake playa in 2013 following a period of inundation similar to that of 2005–06; however, the 2013 cracks formed a network rather than a single feature. Groundwater pumping resulted in more than 27 meters (90 feet) of water-level declines in wells north of the Bicycle Lake playa from 1990 to 2017, and InSAR interferograms indicated more than 400 millimeters (16 inches) of subsidence during 1993–2015. Subsidence rates calculated from InSAR interferograms were variable, temporally and spatially. Results of leveling surveys indicated differential subsidence between 2009 and 2016 along a transect across the playa to the area of maximum subsidence; there was less subsidence south of the 2005–06 crack than north of it. The steepest subsidence gradient for this time was in the area of maximum subsidence. Repeat tape-extensometer measurements from April 2009 to November 2017 across the main-runway crack indicated slight opening along part of the crack. A baseline lidar (light detection and ranging) survey of the main-runway crack was done in January 2009 to track the development of this feature. During the 95-months from January 2009 to December 2016, 0.5 meters (19.7 inches) of subsidence was observed next to the subsidence pit in the crack along the western end of the scanned area, presumably due to erosion and backfilling when the lake was inundated. Standing water in the crack along the eastern end of the scanned area prevented determining change for that part of the crack. Time-series, shaded-relief images show the progression of change in the crack from January 2009 through December 2016, with the crack closing and opening, depending on the time of most recent inundation relative to when the scans were taken. Volumetric changes in the size of the crack, calculated between the sequential lidar surveys, show that from 2009 to 2016, the cumulative volumetric change was 1.5 cubic meters (54 cubic feet). This volume is a minimum estimate because the crack volume below the water surface in the pits was not included in the calculations. Subsurface imaging of the 2005–06 main-runway crack by a series of electrical resistivity tomography surveys in 2008 and 2017 indicated that the crack could extend 5 meters (16 feet) in the subsurface. Electromagnetic induction surveys in 2008, 2014, and 2015 evaluated the technique for ground-failure monitoring and showed that the technique was effective at revealing anomalies correlated with the features of concern in the playa area, such as the 2005–06 crack, numerous ”healed” macropolygon features, and scattered sink-like depressions. Results from numerical experiments simulating watertable decline at depth indicated that the material deposits in the desaturating capillary fringe zone might not be able to transmit large enough stresses up through the overlying soil to cause cracks at the land surface. Results from simulations of desiccation in the presence of a regional tectonically induced stress field, however, tended to support the hypothesis that the combined processes could control the formation of giant desiccation macropolygons in certain areas, such as Bicycle Lake playa.

California

Land subsidence along the California Aqueduct in west-central San Joaquin Valley, California, 2003–10

Extensive groundwater withdrawal from the unconsolidated deposits in the San Joaquin Valley caused widespread aquifer-system compaction and resultant land subsidence from 1926 to 1970—locally exceeding 8.5 meters. The importation of surface water beginning in the early 1950s through the Delta-Mendota Canal and in the early 1970s through the California Aqueduct resulted in decreased groundwater pumping, recovery of water levels, and a reduced rate of compaction in some areas of the San Joaquin Valley. However, drought conditions during 1976–77, 1987–92, and drought conditions and operational reductions in surface-water deliveries during 2007–10 decreased surface-water availability, causing pumping to increase, water levels to decline, and renewed compaction. Land subsidence from this compaction has reduced freeboard and flow capacity of the California Aqueduct, Delta-Mendota Canal, and other canals that deliver irrigation water and transport floodwater. The U.S. Geological Survey, in cooperation with the California Department of Water Resources, assessed more recent land subsidence near a 145-kilometer reach of the California Aqueduct in the west-central part of the San Joaquin Valley as part of an effort to minimize future subsidence-related damages to the California Aqueduct. The location, magnitude, and stress regime of land-surface deformation during 2003–10 were determined by using data and analyses associated with extensometers, Global Positioning System surveys, Interferometric Synthetic Aperture Radar, spirit-leveling surveys, and groundwater wells. Comparison of continuous Global Positioning System, shallow-extensometer, and groundwater-level data indicated that most of the compaction in this area took place beneath the Corcoran Clay, the primary regional confining unit. The integration of measurements strengthens confidence in individual measurement methods and provides the information at spatial and temporal scales that water managers need to design and implement groundwater sustainability plans in compliance with California’s Sustainable Groundwater Management Act. Measurements of land-surface deformation during 2003–10 indicated that the parts of the California Aqueduct closest to the Coast Ranges in the west-central part of the San Joaquin Valley were fairly stable or minimally subsiding on an annual basis; some areas show seasonal periods of subsidence and uplift that resulted in little or no longer-term elevation loss. Many groundwater levels in these areas did not reach historical lows during 2003–10, indicating that deformation nearest the Coast Ranges was likely primarily elastic. Land-surface deformation measurements indicated that some parts of the California Aqueduct that traverse farther from the Coast Ranges toward the valley center subsided. Some parts of the California Aqueduct subsided locally, but generally the California Aqueduct is within part of a 12,000-square-kilometer area affected by 25 millimeters or more of subsidence during 2008–10, with maxima in Madera County, south of the town of El Nido near the San Joaquin River and the Eastside Bypass (540 millimeters), and in Tulare County, west of the town of Pixley (345 millimeters). Interferometric Synthetic Aperture Radar-derived subsidence maps for various periods during 2003–10 show that the area of maximum active subsidence (that is, the largest rates of subsidence) shifted from its historical (1926–70) location southwest of the town of Mendota to these areas nearer the valley center. Calculations indicated that the subsidence rate doubled in 2008 in parts of the study area. Water levels declined during 2007–10 in many shallow and deep wells in the most rapidly subsiding areas, where water levels in many deep wells reached their historical lows, indicating that subsidence measured during this period was largely inelastic. Continued groundwater-level and land-subsidence monitoring in the San Joaquin Valley is important because (1) operational- and drought-related reductions in surface-water deliveries since 1976 have resulted in increased groundwater pumping and associated water-level declines and land subsidence, (2) land use and associated pumping continue to change throughout the valley, and (3) subsidence management is stipulated in the Sustainable Groundwater Management Act. The availability of surface water remains uncertain; even during record-setting precipitation years, such as 2010–11, water deliveries fell short of requests and groundwater pumping was required to meet the irrigation demand. In some areas, the infrastructure is not available to supply surface water, and groundwater is the only source of water. Because of the expected continued demand for water and the limitations and uncertainty of surface-water supplies, groundwater pumping and associated land subsidence remains a concern. Spatially detailed information on land subsidence is needed to minimize future subsidence-related damages to the California Aqueduct and other infrastructure in the San Joaquin Valley, as well as alterations to natural resources such as stream gradients, water depths, and water temperatures. The integration of data on land-surface elevation, subsurface deformation, and water levels—particularly continuous measurements—enables the analysis of aquifer-system response to groundwater pumping, which in turn, enables estimation of the preconsolidation head and calculation of aquifer-system storage properties. This information can be used to improve numerical model simulations of groundwater flow and aquifer-system compaction and allow for consideration of land subsidence in the evaluation of water resource management alternatives and compliance with the Sustainable Groundwater Management Act.

California

Land subsidence and recovery in the Albuquerque Basin, New Mexico, 1993–2014

The Albuquerque Bernalillo County Water Utility Authority (ABCWUA) drinking water supply was almost exclusively sourced from groundwater from within the Albuquerque Basin before 2008. In 2008, the San Juan-Chama Drinking Water Project (SJCDWP) provided surface-water resources to augment the groundwater supply, allowing for a reduction in groundwater pumping in the Albuquerque Basin. In 2013, the U.S. Geological Survey, in cooperation with the ABCWUA, began a study to measure and compare aquifer-system and land-surface elevation change before and after the SJCDWP in 2008. Three methods of data collection with different temporal and spatial resolutions were used for this study: (1) aquifer-system compaction data collected continuously at a single extensometer from 1994 to 2013; (2) land-surface elevation change from Global Positioning System (GPS) surveys of a network of monuments collected in 1994–95, 2005, and 2014; and (3) spatially distributed Interferometric Synthetic Aperture Radar (InSAR) satellite data from 1993 to 2010. Collection of extensometer data allows for direct and continuous measurement of aquifer-system compaction at the extensometer location. The GPS surveys of a network of monuments allow for periodic measurements of land-surface elevation change at monument locations. Interferograms are limited in time by lifespan of the satellite, orbital pattern, and data quality but allow for measurement of gridded land-surface elevation change over the study area. Each of these methods was employed to provide a better understanding of aquifer-system compaction and land-surface elevation change for the Albuquerque Basin. Results do not show large magnitudes of subsidence in the Albuquerque Basin. High temporal-resolution but low spatial-resolution data measurements of aquifer-system compaction at the Albuquerque extensometer show elastic aquifer-system response to recovering groundwater levels. Results from the GPS survey of the network of monuments show inconsistent land-surface elevation changes over the Albuquerque Basin, likely because of the lack of significant change and the complexity of subsurface stratigraphy in addition to the spatial and temporal heterogeneity of groundwater withdrawals over the study period. Results from the InSAR analysis show areas of land-surface elevation increase after 2008, which could be attributed to elastic recovery of the aquifer system. The spatial extent to which elastic recovery of the aquifer system has resulted in recovery of land-surface elevation is limited to the in-situ measurements at the extensometer. Examination of spatially distributed InSAR data relative to limited spatial extent of the complex heterogeneity subsurface stratigraphy may explain some of the heterogeneity of land-surface elevation changes over this study period.

New Mexico

Land subsidence in the southwestern Mojave Desert, California, 1992–2009

Groundwater has been the primary source of domestic, agricultural, and municipal water supplies in the southwestern Mojave Desert, California, since the early 1900s. Increased demands on water supplies have caused groundwater-level declines of more than 100 feet (ft) in some areas of this desert between the 1950s and the 1990s (Stamos and others, 2001; Sneed and others, 2003). These water-level declines have caused the aquifer system to compact, resulting in land subsidence. Differential land subsidence (subsidence occurring at different rates across the landscape) can alter surface drainage routes and damage surface and subsurface infrastructure. For example, fissuring across State Route 247 at Lucerne Lake has required repairs as has pipeline infrastructure near Troy Lake. Land subsidence within the Mojave River and Morongo Groundwater Basins of the southwestern Mojave Desert has been evaluated using InSAR, ground-based measurements, geology, and analyses of water levels between 1992 and 2009 (years in which InSAR data were collected). The results of the analyses were published in three USGS reports— Sneed and others (2003), Stamos and others (2007), and Solt and Sneed (2014). Results from the latter two reports were integrated with results from other USGS/ MWA cooperative groundwater studies into the broader scoped USGS Mojave Groundwater Resources Web site ( http://ca.water.usgs.gov/ mojave/ ). This fact sheet combines the detailed analyses from the three subsidence reports, distills them into a longer-term context, and provides an assessment of options for future monitoring.

California

Water-resources and land-surface deformation evaluation studies at Fort Irwin National Training Center, Mojave Desert, California

The U.S. Army Fort Irwin National Training Center (NTC), in the Mojave Desert, obtains all of its potable water supply from three groundwater basins (Irwin, Langford, and Bicycle) within the NTC boundaries (fig. 1; California Department of Water Resources, 2003). Because of increasing water demands at the NTC, the U.S. Geological Survey (USGS), in cooperation with the U.S. Army, completed several studies to evaluate water resources in the developed and undeveloped groundwater basins underlying the NTC. In all of the developed basins, groundwater withdrawals exceed natural recharge, resulting in water-level declines. However, artificial recharge of treated wastewater has had some success in offsetting water-level declines in Irwin Basin. Additionally, localized water-quality changes have occurred in some parts of Irwin Basin as a result of human activities (i.e., wastewater disposal practices, landscape irrigation, and/or leaking pipes). As part of the multi-faceted NTC-wide studies, traditional datacollection methods were used and include lithological and geophysical logging at newly drilled boreholes, hydrologic data collection (i.e. water-level, water-quality, aquifer tests, wellbore flow). Because these data cover a small portion of the 1,177 square-mile (mi2 ) NTC, regional mapping, including geologic, gravity, aeromagnetic, and InSAR, also were done. In addition, ground and airborne electromagnetic surveys were completed and analyzed to provide more detailed subsurface information on a regional, base-wide scale. The traditional and regional ground and airborne data are being analyzed and will be used to help develop preliminary hydrogeologic framework and groundwater-flow models in all basins. This report is intended to provide an overview of recent water-resources and land-surface deformation studies at the NTC.

California

Hydrogeology, hydrologic effects of development, and simulation of groundwater flow in the Borrego Valley, San Diego County, California

Executive Summary The Borrego Valley is a small valley (110 square miles) in the northeastern part of San Diego County, California. Although the valley is about 60 miles northeast of city of San Diego, it is separated from the Pacific Ocean coast by the mountains to the west and is mostly within the boundaries of Anza-Borrego Desert State Park. From the time the basin was first settled, groundwater has been the only source of water to the valley. Groundwater is used for agricultural, recreational, and municipal purposes. Over time, groundwater withdrawal through pumping has exceeded the amount of water that has been replenished, causing groundwater-level declines of more than 100 feet in some parts of the basin. Continued pumping has resulted in an increase in pumping lifts, reduced well efficiency, dry wells, changes in water quality, and loss of natural groundwater discharge. As a result, the U.S. Geological Survey began a cooperative study of the Borrego Valley with the Borrego Water District (BWD) in 2009. The purpose of the study was to develop a greater understanding of the hydrogeology of the Borrego Valley Groundwater Basin (BVGB) and to provide tools to help evaluate the potential hydrologic effects of future development. The objectives of the study were to (1) improve the understanding of groundwater conditions and land subsidence, (2) incorporate this improved understanding into a model that would assist in the management of the groundwater resources in the Borrego Valley, and (3) use this model to test several management scenarios. This model provides the capability for the BWD and regional stakeholders to quantify the relative benefits of various options for increasing groundwater storage. The study focuses on the period 1945–2010, with scenarios 50 years into the future. This report documents and presents (1) an analysis of the conceptual model, (2) a description of the hydrologic features, (3) a compilation and analysis of water-quality data, (4) the measurement and analysis of land subsidence by using geophysical and remote sensing techniques, (5) the development and calibration of a two-dimensional borehole-groundwater-flow model to estimate aquifer hydraulic conductivities, (6) the development and calibration of a three-dimensional (3-D) integrated hydrologic flow model, (7) a water-availability analysis with respect to current climate variability and land use, and (8) potential future management scenarios. The integrated hydrologic model, referred to here as the “Borrego Valley Hydrologic Model” (BVHM), is a tool that can provide results with the accuracy needed for making water-management decisions, although potential future refinements and enhancements could further improve the level of spatial and temporal resolution and model accuracy. Because the model incorporates time-varying inflows and outflows, this tool can be used to evaluate the effects of temporal changes in recharge and pumping and to compare the relative effects of different water-management scenarios on the aquifer system. Overall, the development of the hydrogeologic and hydrologic models, data networks, and hydrologic analysis provides a basis for assessing surface and groundwater availability and potential water-resource management guidelines.

California

Water availability and land subsidence in the Central Valley, California, USA

The Central Valley in California (USA) covers about 52,000 km 2 and is one of the most productive agricultural regions in the world. This agriculture relies heavily on surface-water diversions and groundwater pumpage to meet irrigation water demand. Because the valley is semi-arid and surface-water availability varies substantially, agriculture relies heavily on local groundwater. In the southern two thirds of the valley, the San Joaquin Valley, historic and recent groundwater pumpage has caused significant and extensive drawdowns, aquifer-system compaction and subsidence. During recent drought periods (2007–2009 and 2012-present), groundwater pumping has increased owing to a combination of decreased surface-water availability and land-use changes. Declining groundwater levels, approaching or surpassing historical low levels, have caused accelerated and renewed compaction and subsidence that likely is mostly permanent. The subsidence has caused operational, maintenance, and construction-design problems for water-delivery and flood-control canals in the San Joaquin Valley. Planning for the effects of continued subsidence in the area is important for water agencies. As land use, managed aquifer recharge, and surface-water availability continue to vary, long-term groundwater-level and subsidence monitoring and modelling are critical to understanding the dynamics of historical and continued groundwater use resulting in additional water-level and groundwater storage declines, and associated subsidence. Modeling tools such as the Central Valley Hydrologic Model, can be used in the evaluation of management strategies to mitigate adverse impacts due to subsidence while also optimizing water availability. This knowledge will be critical for successful implementation of recent legislation aimed toward sustainable groundwater use.

California

Land subsidence in the San Joaquin Valley, California, USA, 2007-14

Rapid land subsidence was recently measured using multiple methods in two areas of the San Joaquin Valley (SJV): between Merced and Fresno (El Nido), and between Fresno and Bakersfield (Pixley). Recent land-use changes and diminished surface-water availability have led to increased groundwater pumping, groundwater-level declines, and land subsidence. Differential land subsidence has reduced the flow capacity of water-conveyance systems in these areas, exacerbating flood hazards and affecting the delivery of irrigation water. Vertical land-surface changes during 2007–2014 were determined by using Interferometric Synthetic Aperture Radar (InSAR), Continuous Global Positioning System (CGPS), and extensometer data. Results of the InSAR analysis indicate that about 7600 km 2 subsided 50–540 mm during 2008–2010; CGPS and extensometer data indicate that these rates continued or accelerated through December 2014. The maximum InSAR-measured rate of 270 mm yr −1 occurred in the El Nido area, and is among the largest rates ever measured in the SJV. In the Pixley area, the maximum InSAR-measured rate during 2008–2010 was 90 mm yr −1 . Groundwater was an important part of the water supply in both areas, and pumping increased when land use changed or when surface water was less available. This increased pumping caused groundwater-level declines to near or below historical lows during the drought periods 2007–2009 and 2012–present. Long-term groundwater-level and land-subsidence monitoring in the SJV is critical for understanding the interconnection of land use, groundwater levels, and subsidence, and evaluating management strategies that help mitigate subsidence hazards to infrastructure while optimizing water supplies.

California

Land subsidence, groundwater levels, and geology in the Coachella Valley, California, 1993-2010

Land subsidence associated with groundwater-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 groundwater resulted in water-level declines as much 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 groundwater 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 groundwater-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 during 1993–2010. The GPS measurements taken at 11 geodetic monuments in 1996 and in 2010 in the southern Coachella Valley indicated that the elevation of the land surface changed –136 to –23 millimeters (mm) ±54 mm (–0.45 to –0.08 feet (ft) ±0.18 ft) during the 14-year period. Changes at 6 of the 11 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 2010. GPS measurements taken at 17 geodetic monuments in 2005 and 2010 indicated that the elevation of the land surface changed –256 to +16 mm ±28 mm (–0.84 to +0.05 ft ±0.09 ft) during the 5-year period. Changes at 5 of the 17 monuments exceeded the maximum expected uncertainty of ±28 mm (±0.09 ft) at the 95-percent confidence level, indicating that subsidence occurred at these monuments between August 2005 and August 2010. At each of these five monuments, subsidence rates were about the same between 2005 and 2010 as between 2000 and 2005. InSAR measurements taken between June 27, 1995, and September 19, 2010, indicated that the land surface subsided from about 220 to 600 mm (0.72 to 1.97 ft) in three areas of the Coachella Valley: near Palm Desert, Indian Wells, and La Quinta. In Palm Desert, the average subsidence rates increased from about 39 millimeters per year (mm/yr), or 0.13 foot per year (ft/yr), during 1995–2000 to about 45 mm/yr (0.15 ft/yr) during 2003–10. In Indian Wells, average subsidence rates for two subsidence maxima were fairly steady at about 34 and 26 mm/yr (0.11 and 0.09 ft/yr) during both periods; for the third maxima, average subsidence rates increased from about 14 to 19 mm/yr (0.05 to 0.06 ft/yr) from the first to the second period. In La Quinta, average subsidence rates for five selected locations ranged from about 17 to 37 mm/yr (0.06 to 0.12 ft/yr) during 1995–2000; three of the locations had similar rates during 2003–mid-2009, while the other two locations had increased subsidence rates. Decreased subsidence rates were calculated throughout the La Quinta subsidence area during mid-2009–10, however, and uplift was observed during 2010 near the southern extent of this area. Water-level measurements taken at wells near the subsiding monuments and in the three subsiding areas shown by InSAR generally indicated that the water levels fluctuated seasonally and declined annually from the early 1990s, or earlier, to 2010; some water levels in 2010 were at the lowest levels in their recorded histories. An exception to annually declining water levels in and near subsiding areas was observed beginning in mid-2009 in the La Quinta subsidence area, where recovering water levels coincided with increased recharge operations at the Thomas E. Levy Recharge Facility; decreased pumpage also could cause groundwater levels to recover. Subsidence concomitant with declining water levels and land-surface uplift concomitant with recovering water levels indicate that aquifer-system compaction could 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 could be permanent.

California

Geology, water-quality, hydrology, and geomechanics of the Cuyama Valley groundwater basin, California, 2008--12

To assess the water resources of the Cuyama Valley groundwater basin in Santa Barbara County, California, a series of cooperative studies were undertaken by the U.S. Geological Survey and the Santa Barbara County Water Agency. Between 2008 and 2012, geologic, water-quality, hydrologic and geomechanical data were collected from selected sites throughout the Cuyama Valley groundwater basin. Geologic data were collected from three multiple-well groundwater monitoring sites and included lithologic descriptions of the drill cuttings, borehole geophysical logs, temperature logs, as well as bulk density and sonic velocity measurements of whole-core samples. Generalized lithologic characterization from the monitoring sites indicated the water-bearing units in the subsurface consist of unconsolidated to partly consolidated sand, gravel, silt, clay, and occasional cobbles within alluvial fan and stream deposits. Analysis of geophysical logs indicated alternating layers of finer- and coarser-grained material that range from less than 1 foot to more than 20 feet thick. On the basis of the geologic data collected, the principal water-bearing units beneath the monitoring-well sites were found to be composed of younger alluvium of Holocene age, older alluvium of Pleistocene age, and the Tertiary-Quaternary Morales Formation. At all three sites, the contact between the recent fill and younger alluvium is approximately 20 feet below land surface. Water-quality samples were collected from 12 monitoring wells, 27 domestic and supply wells, 2 springs, and 4 surface-water sites and were analyzed for a variety of constituents that differed by site, but, in general, included trace elements; nutrients; dissolved organic carbon; major and minor ions; silica; total dissolved solids; alkalinity; total arsenic and iron; arsenic, chromium, and iron species; and isotopic tracers, including the stable isotopes of hydrogen and oxygen, activities of tritium, and carbon-14 abundance. Of the 39 wells sampled, concentrations of total dissolved solids and sulfate from 38 and 37 well samples, respectively, were greater than the U.S. Environmental Protection Agency’s secondary maximum contaminant levels. Concentrations greater than the maximum contaminant levels for nitrate were observed in five wells and were observed for arsenic in four wells. Differences in the stable-isotopic values of hydrogen and oxygen among groundwater samples indicated that water does not move freely between different formations or between different zones within the Cuyama Valley. Variations in isotopic composition indicated that recharge is derived from several different sources. The age of the groundwater, expressed as time since recharge, was between 600 and 38,000 years before present. Detectable concentrations of tritium indicated that younger water, recharged since the early 1950s, is present in parts of the groundwater basin. Hydrologic data were collected from 12 monitoring wells, 56 domestic and supply wells, 3 surface-water sites, and 4 rainfall-gaging stations. Rainfall in the valley averaged about 8 inches annually, whereas the mountains to the south received between 12 and 19 inches. Stream discharge records showed seasonal variability in surface-water flows ranging from no-flow to over 1,500 cubic feet per second. During periods when inflow to the valley exceeds outflow, there is potential recharge from stream losses to the groundwater system Water-level records included manual quarterly depth-to-water measurements collected from 68 wells, time-series data collected from 20 of those wells, and historic water levels from 16 wells. Hydrographs of the manual measurements showed declining water levels in 16 wells, mostly in the South-Main zone, and rising water levels in 14 wells, mostly in the Southern Ventucopa Uplands. Time-series hydrographs showed daily, seasonal, and longer-term effects associated with local pumping. Water-level data from the multiple-well monitoring sites indicated seasonal fluctuations as great as 80 feet and water-level differences between aquifers as great as 40 feet during peak pumping season. Hydrographs from the multiple-well groundwater monitoring sites showed vertical hydraulic gradients were upward during the winter months and downward during the irrigation season. Historic hydrographs showed water-level declines in the Southern-Main, Western Basin, Caliente Northern-Main, and Southern Sierra Madre zone ranging from 1 to 7 feet per year. Hydrographs of wells in the Southern Ventucopa Uplands zone showed several years with marked increases in water levels that corresponded to increased precipitation in the Cuyama Valley. Investigation of hydraulic properties included hydraulic conductivity and transmissivity estimated from aquifer tests performed on 63 wells. Estimates of horizontal hydraulic conductivity ranged from about 1.5 to 28 feet per day and decreased with depth. The median estimated hydraulic conductivity for the older alluvium was about five times that estimated for the Morales Formation. Estimates of transmissivity ranged from 560 to 163,400 gallons per day per foot and decreased with depth. The median estimated transmissivity for the younger alluvium was about three times that estimated for the older alluvium. Geomechanical analysis included land-surface elevation changes at five continuously operating global positioning systems (GPS) and land-subsidence detection at five interferometric synthetic aperture radar (InSAR) reference points. Analysis of data collected from continuously operating GPS stations showed the mountains to the south and west moved upward about 1 millimeter (mm) annually, whereas the station in the center of the Southern-Main zone moved downward more than 7 mm annually, indicating subsidence. It is likely that this subsidence is inelastic (permanent) deformation and indicates reduced storage capacity in the aquifer sediments. Analysis of InSAR data showed local and regional changes that appeared to be dependent, in part, on the time span of the interferogram, seasonal variations in pumping, and tectonic uplift. Long-term InSAR time series showed a total maximum detected subsidence rate of approximately 12 mm per year at one location and approximately 8 mm per year at a second location, while short-term InSAR time series showed maximum subsidence of about 15 mm at one location and localized maximum uplift of about 10 mm at another location.

California