Devin L. Galloway

Sensing the ups and downs of Las Vegas: InSAR reveals structural control of land subsidence and aquifer-system deformation

Land subsidence in Las Vegas, Nevada, United States, between April 1992 and December 1997 was measured using spaceborne interferometric synthetic aperture radar. The detailed deformation maps clearly show that the spatial extent of subsidence is controlled by geologic structures (faults) and sediment composition (clay thickness). The maximum detected subsidence during the 5.75 yr period is 19 cm. Comparison with leveling data indicates that the subsidence rates declined during the past decade as a result of rising ground-water levels brought about by a net reduction in ground-water extraction. Temporal analysis also detects seasonal subsidence and uplift patterns, which provide information about the elastic and inelastic properties of the aquifer system and their spatial variability.

DETECTION OF AQUIFER SYSTEM COMPACTION AND LAND SUBSIDENCE USING INTERFEROMETRIC SYNTHETIC APERTURE RADAR, ANTELOPE VALLEY, MOJAVE DESERT, CALIFORNIA

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.

Seasonal subsidence and rebound in Las Vegas Valley, Nevada, observed by Synthetic Aperture Radar Interferometry

Analyses of areal variations in the subsidence and rebound occurring over stressed aquifer systems, in conjunction with measurements of the hydraulic head fluctuations causing these displacements, can yield valuable information about the compressibility and storage properties of the aquifer system. Historically, stress-strain relationships have been derived from paired extensometer/piezometer installations, which provide only point source data. Because of the general unavailability of spatially detailed deformation data, areal stress-strain relations and their variability are not commonly considered in constraining conceptual and numerical models of aquifer systems. Interferometric synthetic aperture radar (InSAR) techniques can map ground displacements at a spatial scale of tens of meters over 100 km wide swaths. InSAR has been used previously to characterize larger magnitude, generally permanent aquifer system compaction and land subsidence at yearly and longer timescales, caused by sustained drawdown of groundwater levels that produces intergranular stresses consistently greater than the maximum historical stress. We present InSAR measurements of the typically small-magnitude, generally recoverable deformations of the Las Vegas Valley aquifer system occurring at seasonal timescales. From these we derive estimates of the elastic storage coefficient for the aquifer system at several locations in Las Vegas Valley. These high-resolution measurements offer great potential for future investigations into the mechanics of aquifer systems and the spatial heterogeneity of aquifer system structure and material properties as well as for monitoring ongoing aquifer system compaction and land subsidence.

Water-level changes induced by local and distant earthquakes at Long Valley caldera, California

Distant as well as local earthquakes have induced groundwater-level changes persisting for days to weeks at Long Valley caldera, California. Four wells open to formations as deep as 300 m have responded to 16 earthquakes, and responses to two earthquakes in the 3-km-deep Long Valley Exploratory Well (LVEW) show that these changes are not limited to weathered or unconsolidated near-surface rocks. All five wells exhibit water-level variations in response to earth tides, indicating they can be used as low-resolution strainmeters. Earthquakes induce gradual water-level changes that increase in amplitude for as long as 30 days, then return more slowly to pre-earthquake levels. The gradual water-level changes are always drops at wells LKT, LVEW, and CH-10B, and always rises at well CW-3. At a dilatometer just outside the caldera, earthquake-induced strain responses consist of either a step followed by a contractional strain-rate increase, or a transient contractional signal that reaches a maximum in about seven days and then returns toward the pre-earthquake value. The sizes of the gradual water-level changes generally increase with earthquake magnitude and decrease with hypocentral distance. Local earthquakes in Long Valley produce coseismic water-level steps; otherwise the responses to local earthquakes and distant earthquakes are indistinguishable. In particular, water-level and strain changes in Long Valley following the 1992 M7.3 Landers earthquake, 450 km distant, closely resemble those initiated by a M4.9 local earthquake on November 22, 1997, during a seismic swarm with features indicative of fluid involvement. At the LKT well, many of the response time histories are identical for 20 days after each earthquake, and can be matched by a theoretical solution giving the pore pressure as a function of time due to diffusion of a nearby, instantaneous, pressure drop. Such pressure drops could be produced by accelerated inflation of the resurgent dome by amounts too small to be detected by the two-color electronic distance-measuring network. Opening-mode displacement in the south moat, inferred to have followed a M4.9 earthquake on November 22, 1997, could also create extensional strain on the dome and lead to water-level changes similar to those following dome inflation. Contractional strain that could account for earthquake-induced water-level rises at the CW-3 well is inconsistent with geodetic observations. We instead attribute these water-level rises to diffusion of elevated fluid pressure localized in the south moat thermal aquifer. For hydraulic diffusivities appropriate to the upper few hundred meters at Long Valley, an influx of material at temperatures of 300°C can thermally generate pressure of 6 m of water or more, an order of magnitude larger than needed to account for the CW-3 water-level rises. If magma or hot aqueous fluid rises to within 1 km of the surface in the eastern part of the south moat, then hydraulic diffusivities are high enough to allow fluid pressure to propagate to CW-3 on the time scale observed. The data indicate that seismic waves from large distant earthquakes can stimulate upward movement of fluid in the hydrothermal system at Long Valley.

Inferences on the hydrothermal system beneath the resurgent dome in Long Valley Caldera, east-central California, USA, from recent pumping tests and geochemical sampling

Abstract Quaternary volcanic unrest has provided heat for episodic hydrothermal circulation in the Long Valley caldera, including the present-day hydrothermal system, which has been active over the past 40 kyr. The most recent period of crustal unrest in this region of east-central California began around 1980 and has included periods of intense seismicity and ground deformation. Uplift totaling more than 0.7 m has been centered on the caldera’s resurgent dome, and is best modeled by a near-vertical ellipsoidal source centered at depths of 6–7 km. Modeling of both deformation and microgravity data now suggests that (1) there are two inflation sources beneath the caldera, a shallower source 7–10 km beneath the resurgent dome and a deeper source ∼15 km beneath the caldera’s south moat and (2) the shallower source may contain components of magmatic brine and gas. The Long Valley Exploration Well (LVEW), completed in 1998 on the resurgent dome, penetrates to a depth of 3 km directly above this shallower source, but bottoms in a zone of 100°C fluid with zero vertical thermal gradient. Although these results preclude extrapolations of temperatures at depths below 3 km, other information obtained from flow tests and fluid sampling at this well indicates the presence of magmatic volatiles and fault-related permeability within the metamorphic basement rocks underlying the volcanic fill. In this paper, we present recently acquired data from LVEW and compare them with information from other drill holes and thermal springs in Long Valley to delineate the likely flow paths and fluid system properties under the resurgent dome. Additional information from mineralogical assemblages in core obtained from fracture zones in LVEW documents a previous period of more vigorous and energetic fluid circulation beneath the resurgent dome. Although this system apparently died off as a result of mineral deposition and cooling (and/or deepening) of magmatic heat sources, flow testing and tidal analyses of LVEW water level data show that relatively high permeability and strain sensitivity still exist in the steeply dipping principal fracture zone penetrated at a depth of 2.6 km. The hydraulic properties of this zone would allow a pressure change induced at distances of several kilometers below the well to be observable within a matter of days. This indicates that continuous fluid pressure monitoring in the well could provide direct evidence of future intrusions of magma or high-temperature fluids at depths of 5–7 km.

Assessing Groundwater Depletion and Dynamics Using GRACE and InSAR: Potential and Limitations

In the last decade, remote sensing of the temporal variation of ground level and gravity has improved our understanding of groundwater dynamics and storage. Mass changes are measured by GRACE (Gravity Recovery and Climate Experiment) satellites, whereas ground deformation is measured by processing synthetic aperture radar satellites data using the InSAR (Interferometry of Synthetic Aperture Radar) techniques. Both methods are complementary and offer different sensitivities to aquifer system processes. GRACE is sensitive to mass changes over large spatial scales (more than 100,000 km2 ). As such, it fails in providing groundwater storage change estimates at local or regional scales relevant to most aquifer systems, and at which most groundwater management schemes are applied. However, InSAR measures ground displacement due to aquifer response to fluid-pressure changes. InSAR applications to groundwater depletion assessments are limited to aquifer systems susceptible to measurable deformation. Furthermore, the inversion of InSAR-derived displacement maps into volume of depleted groundwater storage (both reversible and largely irreversible) is confounded by vertical and horizontal variability of sediment compressibility. During the last decade, both techniques have shown increasing interest in the scientific community to complement available in situ observations where they are insufficient. In this review, we present the theoretical and conceptual bases of each method, and present idealized scenarios to highlight the potential benefits and challenges of combining these techniques to remotely assess groundwater storage changes and other aspects of the dynamics of aquifer systems.

Coastal subsidence and relative sea level rise

Subsurface fluid-pressure declines caused by pumping of groundwater or hydrocarbons can lead to aquifer-system compaction and consequent land subsidence. This subsidence can be rapid, as much as 30 cm per year in some instances, and large, totaling more than 13 m in extreme examples. Thus anthropogenic subsidence may be the dominant contributor to relative sea-level rise in coastal environments where subsurface fluids are heavily exploited. Maximum observed rates of human-induced subsidence greatly exceed the rates of natural subsidence of unconsolidated sediments (∼0.1–1cmyr �1 ) and the estimated rates of ongoing global sea-level rise (∼0.3 cm yr �1 ).

Hydraulic Characterization of Overpressured Tuffs in Central Yucca Flat, Nevada Test Site, Nye County, Nevada

A sequence of buried, bedded, air-fall tuffs has been used extensively as a host medium for underground nuclear tests detonated in the central part of Yucca Flat at the Nevada Test Site. Water levels within these bedded tuffs have been elevated hundreds of meters in areas where underground nuclear tests were detonated below the water table. Changes in the ground-water levels within these tuffs and changes in the rate and distribution of land-surface subsidence above these tuffs indicate that pore-fluid pressures have been slowly depressurizing since the cessation of nuclear testing in 1992. Declines in ground-water levels concurrent with regional land subsidence are explained by poroelastic deformation accompanying ground-water flow as fluids pressurized by underground nuclear detonations drain from the host tuffs into the overlying water table and underlying regional carbonate aquifer. A hydraulic conductivity of about 3 x 10-6 m/d and a specific storage of 9 x 10-6 m-1 are estimated using ground-water flow models. Cross-sectional and three-dimensional ground-water flow models were calibrated to measured water levels and to land-subsidence rates measured using Interferometric Synthetic Aperture Radar. Model results are consistent and indicate that about 2 million m3 of ground water flowed from the tuffs to the carbonate rock as a result of pressurization caused by underground nuclear testing. The annual rate of inflow into the carbonate rock averaged about 0.008 m/yr between 1962 and 2005, and declined from 0.005 m/yr in 2005 to 0.0005 m/yr by 2300.

Preface: Land subsidence processes

Both natural and anthropogenic land subsidence are global phenomena caused by a variety of factors, many of which are related to hydrogeologic processes. Common natural subsidence processes include consolidation related to sediment loading, tectonics, volcanism, and dissolution of relatively soluble carbonate and evaporite minerals. Some natural subsidence processes are directly influenced by human activities related to land and water use and by climatic variability. The development of water resources to support human habitation and cultivation for agriculture typically results in the use and diversion of available surface-water supplies and a reliance on groundwater supplies. These practices can alter the natural hydrologic system in ways that amplify natural subsidence processes or create new anthropogenic subsidence. For example, engineered diversion of runoff can focus recharge in areas susceptible to mineral dissolution which can lead to sinkholes or other collapse features in the karst landscape, or engineered drainage of wetlands or saturated organic soils can cause oxidation and consolidation of the soils. Anthropogenic groundwater abstraction from susceptible (generally, unconsolidated alluvial, fluvial and lacustrine sediments) aquifer systems to support water use principally for agriculture, municipal-industrial and energy development typically can lead to local and regional groundwater storage depletion and accompanying aquifersystem compaction and land subsidence related to increases in effective stresses caused by groundwater-level declines. Note that the terms ‘compaction’, commonly used by geologists, and ‘consolidation’, commonly used in soil mechanics, are used interchangeably in this preface and theme issue. Climate variation in terms of global warming, whether natural or anthropogenic, can indirectly cause glacial isostatic adjustments (uplift and subsidence) of the Earth’s crust related to melting of ice sheets, or can thaw permafrost with subsequent loss of ice volume and drainage of shallow groundwater leading to mechanical and even oxidation mediated subsidence. Climate variation may result in either reductions (droughts) or enhancements (wet periods) of precipitation, surface runoff and groundwater recharge. These reductions can cause subsidence owing to lowered groundwater levels contributing to aquifersystem compaction and to oxidation and consolidation of organic soils; enhancements can cause subsidence owing to increased dissolution of karst minerals and reduced mechanical support for pre-existing karst features. The 14 articles (13 papers and one essay) constituting this theme issue address several of the subsidence processes where human use of natural resources and climate variability have combined to create critical anthropogenic land subsidence problems: oxidation and consolidation of organic soils (three articles), dissolution and collapse of carbonate and evaporite rocks (karst) (three articles), thawing permafrost (thermokarst) (one article) and aquifer-system compaction (seven articles). Each of these subsidence processes is intricately related to the Published in the theme issue BLand Subsidence Processes^

Time-varying land subsidence detected by radar altimetry: California, Taiwan and north China

Contemporary applications of radar altimetry include sea-level rise, ocean circulation, marine gravity, and icesheet elevation change. Unlike InSAR and GNSS, which are widely used to map surface deformation, altimetry is neither reliant on highly temporally-correlated ground features nor as limited by the available spatial coverage, and can provide long-term temporal subsidence monitoring capability. Here we use multi-mission radar altimetry with an approximately 23 year data-span to quantify land subsidence in cropland areas. Subsidence rates from TOPEX/POSEIDON, JASON-1, ENVISAT, and JASON-2 during 1992–2015 show time-varying trends with respect to displacement over time in California’s San Joaquin Valley and central Taiwan, possibly related to changes in land use, climatic conditions (drought) and regulatory measures affecting groundwater use. Near Hanford, California, subsidence rates reach 18 cm yr−1 with a cumulative subsidence of 206 cm, which potentially could adversely affect operations of the planned California High-Speed Rail. The maximum subsidence rate in central Taiwan is 8 cm yr−1. Radar altimetry also reveals time-varying subsidence in the North China Plain consistent with the declines of groundwater storage and existing water infrastructure detected by the Gravity Recovery And Climate Experiment (GRACE) satellites, with rates reaching 20 cm yr−1 and cumulative subsidence as much as 155 cm.