Earth's sinking surface.

China's major cities show considerable subsidence from human activities.

Disappearing cities on US coasts.

The sea level along the US coastlines is projected to rise by 0.25-0.3 m by 2050, increasing the probability of more destructive flooding and inundation in major cities1-3. However, these impacts may be exacerbated by coastal subsidence-the sinking of coastal land areas4-a factor that is often underrepresented in coastal-management policies and long-term urban planning2,5. In this study, we combine high-resolution vertical land motion (that is, raising or lowering of land) and elevation datasets with projections of sea-level rise to quantify the potential inundated areas in 32 major US coastal cities. Here we show that, even when considering the current coastal-defence structures, further land area of between 1,006 and 1,389 km2 is threatened by relative sea-level rise by 2050, posing a threat to a population of 55,000-273,000 people and 31,000-171,000 properties. Our analysis shows that not accounting for spatially variable land subsidence within the cities may lead to inaccurate projections of expected exposure. These potential consequences show the scale of the adaptation challenge, which is not appreciated in most US coastal cities.

Slowly but surely: Exposure of communities and infrastructure to subsidence on the US east coast.

Coastal communities are vulnerable to multihazards, which are exacerbated by land subsidence. On the US east coast, the high density of population and assets amplifies the region's exposure to coastal hazards. We utilized measurements of vertical land motion rates obtained from analysis of radar datasets to evaluate the subsidence-hazard exposure to population, assets, and infrastructure systems/facilities along the US east coast. Here, we show that 2,000 to 74,000 km2 land area, 1.2 to 14 million people, 476,000 to 6.3 million properties, and >50% of infrastructures in major cities such as New York, Baltimore, and Norfolk are exposed to subsidence rates between 1 and 2 mm per year. Additionally, our analysis indicates a notable trend: as subsidence rates increase, the extent of area exposed to these hazards correspondingly decreases. Our analysis has far-reaching implications for community and infrastructure resilience planning, emphasizing the need for a targeted approach in transitioning from reactive to proactive hazard mitigation strategies in the era of climate change.

Hidden vulnerability of US Atlantic coast to sea-level rise due to vertical land motion.

The vulnerability of coastal environments to sea-level rise varies spatially, particularly due to local land subsidence. However, high-resolution observations and models of coastal subsidence are scarce, hindering an accurate vulnerability assessment. We use satellite data from 2007 to 2020 to create high-resolution map of subsidence rate at mm-level accuracy for different land covers along the ~3,500 km long US Atlantic coast. Here, we show that subsidence rate exceeding 3 mm per year affects most coastal areas, including wetlands, forests, agricultural areas, and developed regions. Coastal marshes represent the dominant land cover type along the US Atlantic coast and are particularly vulnerable to subsidence. We estimate that 58 to 100% of coastal marshes are losing elevation relative to sea level and show that previous studies substantially underestimate marsh vulnerability by not fully accounting for subsidence.

Widespread deep seismicity in the Delaware Basin, Texas, is mainly driven by shallow wastewater injection.

Industrial activity away from plate boundaries can induce earthquakes and has evolved into a global issue. Much of the induced seismicity in the United States' midcontinent is attributed to a direct pressure increase from deep wastewater disposal. This mechanism is not applicable where deep basement faults are hydraulically isolated from shallow injection aquifers, leading to a debate about the mechanisms for induced seismicity. Here, we compile industrial, seismic, geodetic, and geological data within the Delaware Basin, western Texas, and calculate stress and pressure changes at seismogenic depth using a coupled poroelastic model. We show that the widespread deep seismicity is mainly driven by shallow wastewater injection through the transmission of poroelastic stresses assuming that unfractured shales are hydraulic barriers over decadal time scales. A zone of seismic quiescence to the north, where injection-induced stress changes would promote seismicity, suggests a regional tectonic control on the occurrence of induced earthquakes. Comparing the poroelastic responses from injection and extraction operations, we find that the basement stress is most sensitive to shallow reservoir hydrogeological parameters, particularly hydraulic diffusivity. These results demonstrate that intraplate seismicity can be caused by shallow human activities that poroelastically perturb stresses at hydraulically isolated seismogenic depths, with impacts on seismicity that are preconditioned by regional tectonics.

Subsidence-Derived Volumetric Strain Models for Mapping Extensional Fissures and Constraining Rock Mechanical Properties in the San Joaquin Valley, California.

Large-scale subsidence due to aquifer-overdraft is an ongoing hazard in the San Joaquin Valley. Subsidence continues to cause damage to infrastructure and increases the risk of extensional fissures.Here, we use InSAR-derived vertical land motion (VLM) to model the volumetric strain rate due to groundwater storage change during the 2007-2010 drought in the San Joaquin Valley, Central Valley, California. We then use this volumetric strain rate model to calculate surface tensile stress in order to predict regions that are at the highest risk for hazardous tensile surface fissures. We find a maximum volumetric strain rate of -232 microstrain/yr at a depth of 0 to 200 m in Tulare and Kings County, California. The highest risk of tensile fissure development occurs at the periphery of the largest subsiding zones, particularly in Tulare County and Merced County. Finally, we assume that subsidence is likely due to aquifer pressure change, which is calculated using groundwater level changes observed at 300 wells during this drought. We combine pressure data from selected wells with our volumetric strain maps to estimate the quasi-static bulk modulus, K, a poroelastic parameter applicable when pressure change within the aquifer is inducing volumetric strain. This parameter is reflective of a slow deformation process and is one to two orders of magnitude lower than typical values for the bulk modulus found using seismic velocity data. The results of this study highlight the importance of large-scale, high-resolution VLM measurements in evaluating aquifer system dynamics, hazards associated with overdraft, and in estimating important poroelastic parameters.

Understanding of Contemporary Regional Sea-Level Change and the Implications for the Future.

Global sea level provides an important indicator of the state of the warming climate, but changes in regional sea level are most relevant for coastal communities around the world. With improvements to the sea-level observing system, the knowledge of regional sea-level change has advanced dramatically in recent years. Satellite measurements coupled with in situ observations have allowed for comprehensive study and improved understanding of the diverse set of drivers that lead to variations in sea level in space and time. Despite the advances, gaps in the understanding of contemporary sea-level change remain and inhibit the ability to predict how the relevant processes may lead to future change. These gaps arise in part due to the complexity of the linkages between the drivers of sea-level change. Here we review the individual processes which lead to sea-level change and then describe how they combine and vary regionally. The intent of the paper is to provide an overview of the current state of understanding of the processes that cause regional sea-level change and to identify and discuss limitations and uncertainty in our understanding of these processes. Areas where the lack of understanding or gaps in knowledge inhibit the ability to provide the needed information for comprehensive planning efforts are of particular focus. Finally, a goal of this paper is to highlight the role of the expanded sea-level observation network-particularly as related to satellite observations-in the improved scientific understanding of the contributors to regional sea-level change.

Tracking California's sinking coast from space: Implications for relative sea-level rise.

Coastal vertical land motion affects projections of sea-level rise, and subsidence exacerbates flooding hazards. Along the ~1350-km California coastline, records of high-resolution vertical land motion rates are scarce due to sparse instrumentation, and hazards to coastal communities are underestimated. Here, we considered a ~100-km-wide swath of land along California's coast and performed a multitemporal interferometric synthetic aperture radar (InSAR) analysis of large datasets, obtaining estimates of vertical land motion rates for California's entire coast at ~100-m dimensions-a ~1000-fold resolution improvement to the previous record. We estimate between 4.3 million and 8.7 million people in California's coastal communities, including 460,000 to 805,000 in San Francisco, 8000 to 2,300,00 in Los Angeles, and 2,000,000 to 2,300,000 in San Diego, are exposed to subsidence. The unprecedented detail and submillimeter accuracy resolved in our vertical land motion dataset can transform the analysis of natural and anthropogenic changes in relative sea-level and associated hazards.

Pore-pressure diffusion, enhanced by poroelastic stresses, controls induced seismicity in Oklahoma.

Induced seismicity linked to geothermal resource exploitation, hydraulic fracturing, and wastewater disposal is evolving into a global issue because of the increasing energy demand. Moderate to large induced earthquakes, causing widespread hazards, are often related to fluid injection into deep permeable formations that are hydraulically connected to the underlying crystalline basement. Using injection data combined with a physics-based linear poroelastic model and rate-and-state friction law, we compute the changes in crustal stress and seismicity rate in Oklahoma. This model can be used to assess earthquake potential on specific fault segments. The regional magnitude-time distribution of the observed magnitude (M) 3+ earthquakes during 2008-2017 is reproducible and is the same for the 2 optimal, conjugate fault orientations suggested for Oklahoma. At the regional scale, the timing of predicted seismicity rate, as opposed to its pattern and amplitude, is insensitive to hydrogeological and nucleation parameters in Oklahoma. Poroelastic stress changes alone have a small effect on the seismic hazard. However, their addition to pore-pressure changes can increase the seismicity rate by 6-fold and 2-fold for central and western Oklahoma, respectively. The injection-rate reduction in 2016 mitigates the exceedance probability of M5.0 by 22% in western Oklahoma, while that of central Oklahoma remains unchanged. A hypothetical injection shut-in in April 2017 causes the earthquake probability to approach its background level by ∼2025. We conclude that stress perturbation on prestressed faults due to pore-pressure diffusion, enhanced by poroelastic effects, is the primary driver of the induced earthquakes in Oklahoma.

Comment on "Short-lived pause in Central California subsidence after heavy winter precipitation of 2017" by K. D. Murray and R. B. Lohman.

In a study by Murray and Lohman (M&L), the authors suggest that remote sensing data are useful for monitoring land subsidence due to aquifer system compaction. We agree. To infer aquifer dynamics, we provide a more detailed and joint analysis of deformation and groundwater data. Investigating well data in the Tulare Basin, we find that groundwater levels stabilized before 2015 and show that M&L's observed continued subsidence through July 2016 is likely caused by the delayed compaction of the aquitard. Our analysis suggests the observed 2017 transient uplift is not due to recharge of the aquifer system after heavy winter rainfall because it requires an unrealistic vertical hydraulic gradient nearly five orders of magnitude larger than that typical of Tulare Basin. We find that, regardless of the amount of rainfall, transient annual uplifts of ~3 cm occur in May to June. Using an elastic skeletal storage coefficient of 5 × 10-3, we link this ground uplift to annual groundwater level changes.

Groundwater Loss and Aquifer System Compaction in San Joaquin Valley During 2012-2015 Drought.

California's millennium drought of 2012-2015 severely impacted the Central Valley aquifer system and caused permanent loss of groundwater and aquifer storage capacity. To quantify these impacts within the southern San Joaquin Valley, we analyze various complementary measurements, including gravity changes from Gravity Recovery and Climate Experiment (GRACE) satellites; vertical land motion from Global Positioning System, interferometric synthetic aperture radar, and extensometer; and groundwater level records. The interferometric data set acquired by the Sentinel-1 satellite only spans the period January 2015 and October 2017, while the other data sets span the entire drought period. Using GRACE observations, we find an average groundwater loss of 6.1 ± 2.3 km3/year as a lower bound estimate for the San Joaquin Valley, amounting to a total volume of 24.2 ± 9.3 km3 lost during the period October 2011 to September 2015. This is consistent with the total volume of 29.25 ± 8.7 km3, estimated using only Global Positioning System deformation data. Our results highlight the advantage of using vertical land motion data to evaluate groundwater loss and thus fill the gaps between GRACE and GRACE-Follow-On missions and complement their estimates. We further determine that 0.4-3.25% of the aquifer system storage capacity is permanently lost during this drought period. Comparing groundwater level and vertical land motion data following September 2015, we determine an equilibration time of 0.5-1.5 years for groundwater levels within aquitard and aquifer units, during which residual compaction of aquitard and land subsidence continues beyond the drought period. We suggest that such studies can advance the knowledge of evolving groundwater resources, enabling managers and decision makers to better assess water demand and supply during and in-between drought periods.

Sustained Groundwater Loss in California's Central Valley Exacerbated by Intense Drought Periods.

The accelerated rate of decline in groundwater levels across California's Central Valley results from overdrafting and low rates of natural recharge and is exacerbated by droughts. The lack of observations with an adequate spatiotemporal resolution to constrain the evolution of groundwater resources poses severe challenges to water management efforts. Here we present SAR interferometric measurements of high-resolution vertical land motion across the valley, revealing multiscale patterns of aquifer hydrogeological properties and groundwater storage change. Investigating the depletion and degradation of the aquifer-system during 2007-2010, when the entire valley experienced a severe drought, we find that ~2% of total aquifer-system storage was permanently lost, owing to irreversible compaction of the system. Over this period, the seasonal groundwater storage change amplitude of 10.11 ± 2.5 km3 modulates a long-term groundwater storage decline of 21.32 ± 7.2 km3. Estimates for subbasins show more complex patterns, most likely associated with local hydrogeology, recharge, demand, and underground flow. Presented measurements of aquifer-system compaction provide a more complete understanding of groundwater dynamics and can potentially be used to improve water security.

Episodic creep events on the San Andreas Fault caused by pore-pressure variations.

Recent seismic and geodetic observations indicate that interseismic creep rate varies in both time and space. The spatial extent of creep pinpoints locked asperities, while its temporary accelerations, known as slow-slip events, may trigger earthquakes. Although the conditions promoting fault creep are well-studied, the mechanisms for initiating episodic slow-slip events are enigmatic. Here we investigate surface deformation measured by radar interferometry along the central San Andreas Fault between 2003 and 2010 to constrain the temporal evolution of creep. We show that slow-slip events are ensembles of localized creep bursts that aseismically rupture isolated fault compartments. Using a rate and state friction model, we show that effective normal stress is temporally variable on the fault, and support this using seismic observations. We propose that, compaction-driven elevated pore fluid pressure in hydraulically isolated fault zone and subsequent frictional dilation cause the observed slow slip episodes. We further suggest that the 2004 Mw6 Parkfield earthquake might have been triggered by a slow-slip event, which increased the Coulomb failure stress by up to 0.45 bar per year. This implies that while creeping segments are suggested to act as seismic rupture barriers, slow-slip events on these zones might promote seismicity on adjacent locked segments.

Global climate change and local land subsidence exacerbate inundation risk to the San Francisco Bay Area.

The current global projections of future sea level rise are the basis for developing inundation hazard maps. However, contributions from spatially variable coastal subsidence have generally not been considered in these projections. We use synthetic aperture radar interferometric measurements and global navigation satellite system data to show subsidence rates of less than 2 mm/year along most of the coastal areas along San Francisco Bay. However, rates exceed 10 mm/year in some areas underlain by compacting artificial landfill and Holocene mud deposits. The maps estimating 100-year inundation hazards solely based on the projection of sea level rise from various emission scenarios underestimate the area at risk of flooding by 3.7 to 90.9%, compared with revised maps that account for the contribution of local land subsidence. Given ongoing land subsidence, we project that an area of 125 to 429 km2 will be vulnerable to inundation, as opposed to 51 to 413 km2 considering sea level rise alone.

Surface uplift and time-dependent seismic hazard due to fluid injection in eastern Texas.

Observations that unequivocally link seismicity and wastewater injection are scarce. Here we show that wastewater injection in eastern Texas causes uplift, detectable in radar interferometric data up to >8 kilometers from the wells. Using measurements of uplift, reported injection data, and a poroelastic model, we computed the crustal strain and pore pressure. We infer that an increase of >1 megapascal in pore pressure in rocks with low compressibility triggers earthquakes, including the 4.8-moment magnitude event that occurred on 17 May 2012, the largest earthquake recorded in eastern Texas. Seismic activity increased even while injection rates declined, owing to diffusion of pore pressure from earlier periods with higher injection rates. Induced seismicity potential is suppressed where tight confining formations prevent pore pressure from propagating into crystalline basement rocks.