3D architecture and complex behavior along the simple central San Andreas fault.

The central San Andreas Fault (CSAF) exhibits a simple linear large-scale fault geometry, yet seismic and aseismic deformation features vary in a complex way along the fault. Here we investigate fault zone behaviors using geodetic observation, seismicity and microearthquake focal mechanisms. We employ an improved focal-mechanism characterization method using relative earthquake radiation patterns on 75,164 Ml ≥ 1 earthquakes along a 2-km-wide, 190-km-long segment of the CSAF, from 1984 to 2015. The data reveal the 3D fine-scale structure and interseismic kinematics of the CSAF. Our findings indicate that the first-order spatial variations in interseismic fault creep rate, creep direction, and the fault zone stress field can be explained by a simple fault coupling model. The inferred 3D mechanical properties of a mechanically weak and poorly coupled fault zone provide a unified understanding of the complex fine-scale kinematics, indicating distributed slip deficits facilitating small-to-moderate earthquakes, localized stress heterogeneities, and complex multi-scale ruptures along the fault. Through this detailed mapping, we aim to relate the fine-scale fault architecture to potential future faulting behavior along the CSAF.

A weak subducting slab at intermediate depths below northeast Japan.

Knowledge of the state of stress in subducting slabs is essential for understanding their mechanical behavior and the physical processes that generate earthquakes. Here, we develop a framework which uses a high-resolution focal mechanism catalog to determine the change in the position of the neutral plane before and after the M9 Tohoku-oki earthquake to determine that the deviatoric stress within the slab at intermediate depths must be very low (∼1 MPa). We show that by combining the static stress calculated from coseismic slip distributions with the stress orientations before and after the mainshock, we can determine the full deviatoric stress tensor within the subducting slab at intermediate depths. These results preclude earthquake source mechanisms that require large background driving stresses, favoring a mechanically weak subducting slab, thus providing quantitative constraints on the physical processes that generate intermediate-depth earthquakes.

Reliable earthquake precursors?

Global Positioning System (GPS) measurements suggest hours-long precursors to many large earthquakes.

Aseismic slip and recent ruptures of persistent asperities along the Alaska-Aleutian subduction zone.

The frictional properties and slip behaviors of subduction thrusts play a key role in seismic and tsunami hazard assessment, especially in weakly coupled "seismic gaps". Here, we rely on GPS observations in the Shumagin Gap of the Aleutian subduction zone to derive the slip distribution of the 2020 Mw 7.8 Simeonof Island, Alaska earthquake and of the subsequent afterslip during the first 87-day period. Our modeling results show that the mainshock ruptured at depths of ∼30-40 km beneath Simeonof Island. Kinematic and stress-driven models indicate that the afterslip occurred both updip and downdip of the mainshock rupture. Physically plausible locking models derived from interseismic GPS velocities suggest that the 2020 Simeonof and 2021 Mw 8.2 Chignik earthquakes ruptured persistent asperities on the subduction thrust. We infer that there are several additional persistent asperities at depths of 20-50 km west ∼157°W. However, it is still uncertain whether there are additional locked asperities at shallow depths because of the current lack of geodetic observations close to the trench.

Bridging earthquakes and mountain building in the Santa Cruz Mountains, CA.

Relative crustal motions along active faults generate earthquakes, and repeated earthquake cycles build mountain ranges over millions of years. However, the long-term summation of elastic, earthquake-related deformation cannot produce the deformation recorded within the rock record. Here, we provide an explanation for this discrepancy by showing that increases in strain facilitated by plastic deformation of Earth's crust during the earthquake cycle, in conjunction with isostatic deflection and erosion, transform relative fault motions that produce individual earthquakes to geologic deformations. We focus our study on the data-rich Santa Cruz Mountains, CA, USA and compare predicted and observed quantities for rock uplift, apatite (U-Th)/He thermochronology, topographic relief, 10Be-based erosion rates, and interseismic surface velocities. This approach reconciles these disparate records of mountain-building processes, allowing us to explicitly bridge decadal measures of deformation with that produced by millions of years of plate motion.

Synchronized and asynchronous modulation of seismicity by hydrological loading: A case study in Taiwan.

Delineation of physical factors that contribute to earthquake triggering is a challenging issue in seismology. We analyze hydrological modulation of seismicity in Taiwan using groundwater level data and GNSS time series. In western Taiwan, the seismicity rate reaches peak levels in February to April and drops to its lowest values in July to September, exhibiting a direct correlation with annual water unloading. The elastic hydrological load cycle may be the primary driving mechanism for the observed synchronized modulation of earthquakes, as also evidenced by deep earthquakes in eastern Taiwan. However, shallow earthquakes in eastern Taiwan (<18 km) are anticorrelated with water unloading, which is not well explained by either hydrological loading, fluid transport, or pore pressure changes and suggests other time-dependent processes. The moderate correlation between stacked monthly trends of large historic earthquakes and present-day seismicity implies a modestly higher seismic hazard during the time of low annual hydrological loading.

A revised position for the primary strand of the Pleistocene-Holocene San Andreas fault in southern California.

The San Andreas fault has the highest calculated time-dependent probability for large-magnitude earthquakes in southern California. However, where the fault is multistranded east of the Los Angeles metropolitan area, it has been uncertain which strand has the fastest slip rate and, therefore, which has the highest probability of a destructive earthquake. Reconstruction of offset Pleistocene-Holocene landforms dated using the uranium-thorium soil carbonate and beryllium-10 surface exposure techniques indicates slip rates of 24.1 ± 3 millimeter per year for the San Andreas fault, with 21.6 ± 2 and 2.5 ± 1 millimeters per year for the Mission Creek and Banning strands, respectively. These data establish the Mission Creek strand as the primary fault bounding the Pacific and North American plates at this latitude and imply that 6 to 9 meters of elastic strain has accumulated along the fault since the most recent surface-rupturing earthquake, highlighting the potential for large earthquakes along this strand.

What's down there? The structures, materials and environment of deep-seated slow slip and tremor.

Deep-seated slow slip and tremor (SST), including slow slip events, episodic tremor and slip, and low-frequency earthquakes, occur downdip of the seismogenic zone of numerous subduction megathrusts and plate boundary strike-slip faults. These events represent a fascinating and perplexing mode of fault failure that has greatly broadened our view of earthquake dynamics. In this contribution, we review constraints on SST deformation processes from both geophysical observations of active subduction zones and geological observations of exhumed field analogues. We first provide an overview of what has been learned about the environment, kinematics and dynamics of SST from geodetic and seismologic data. We then describe the materials, deformation mechanisms, and metamorphic and fluid pressure conditions that characterize exhumed rocks from SST source depths. Both the geophysical and geological records strongly suggest the importance of a fluid-rich and high fluid pressure habitat for the SST source region. Additionally, transient deformation features preserved in the rock record, involving combined frictional-viscous shear in regions of mixed lithology and near-lithostatic fluid pressures, may scale with the tremor component of SST. While several open questions remain, it is clear that improved constraints on the materials, environment, structure, and conditions of the plate interface from geophysical imaging and geologic observations will enhance model representations of the boundary conditions and geometry of the SST deformation process. This article is part of a discussion meeting issue 'Understanding earthquakes using the geological record'.

Four-dimensional surface motions of the Slumgullion landslide and quantification of hydrometeorological forcing.

Landslides modify the natural landscape and cause fatalities and property damage worldwide. Quantifying landslide dynamics is challenging due to the stochastic nature of the environment. With its large area of ~1 km2 and perennial motions at ~10-20 mm per day, the Slumgullion landslide in Colorado, USA, represents an ideal natural laboratory to better understand landslide behavior. Here, we use hybrid remote sensing data and methods to recover the four-dimensional surface motions during 2011-2018. We refine the boundaries of an area of ~0.35 km2 below the crest of the prehistoric landslide. We construct a mechanical framework to quantify the rheology, subsurface channel geometry, mass flow rate, and spatiotemporally dependent pore-water pressure feedback through a joint analysis of displacement and hydrometeorological measurements from ground, air and space. Our study demonstrates the importance of remotely characterizing often inaccessible, dangerous slopes to better understand landslides and other quasi-static mass fluxes in natural and industrial environments, which will ultimately help reduce associated hazards.

Slow slip events in the roots of the San Andreas fault.

Episodic tremor and accompanying slow slip are observed at the down-dip edge of subduction seismogenic zones. While tremors are the seismic signature of this phenomenon, they correspond to a small fraction of the moment released; thus, the associated fault slip can be quantified only by geodetic observations. On continental strike-slip faults, tremors have been observed in the roots of the Parkfield segment of the San Andreas fault. However, associated transient aseismic slip has never been detected. By making use of the timing of transient tremor activity and the dense Parkfield-area global positioning system network, we can detect deep slow slip events (SSEs) at 16-km depth on the Parkfield segment with an average moment equivalent to M w 4.90 ± 0.08. Characterization of transient SSEs below the Parkfield locked asperity, at the transition with the creeping section of the San Andreas fault, provides new constraints on the seismic cycle in this region.

A shift from drought to extreme rainfall drives a stable landslide to catastrophic failure.

The addition of water on or below the earth's surface generates changes in stress that can trigger both stable and unstable sliding of landslides and faults. While these sliding behaviours are well-described by commonly used mechanical models developed from laboratory testing (e.g., critical-state soil mechanics and rate-and-state friction), less is known about the field-scale environmental conditions or kinematic behaviours that occur during the transition from stable to unstable sliding. Here we use radar interferometry (InSAR) and a simple 1D hydrological model to characterize 8 years of stable sliding of the Mud Creek landslide, California, USA, prior to its rapid acceleration and catastrophic failure on May 20, 2017. Our results suggest a large increase in pore-fluid pressure occurred during a shift from historic drought to record rainfall that triggered a large increase in velocity and drove slip localization, overcoming the stabilizing mechanisms that had previously inhibited landslide acceleration. Given the predicted increase in precipitation extremes with a warming climate, we expect it to become more common for landslides to transition from stable to unstable motion, and therefore a better assessment of this destabilization process is required to prevent loss of life and infrastructure.

Seasonal modulation of deep slow-slip and earthquakes on the Main Himalayan Thrust.

The interaction between seasonally-induced non-tectonic and tectonic deformation along the Himalayan plate boundary remains debated. Here, we propose that tectonic deformation along this plate boundary can be significantly influenced by the deformation induced by the non-tectonic hydrological loading cycles. We explore seasonal mass oscillations by continental water storage in Southeast Asia and Himalayan arc region using continuous Global Positioning System measurements and satellite data from the Gravity Recovery and Climate Experiment. We suggest that the substantially higher transient displacements above the base of the seismogenic zone indicate a role of changes in aseismic slip rate on the deep megathrust that may be controlled by seasonal hydrological loading. We invoke modulation of aseismic slip on the megathrust down-dip of the seismogenic zone due to a fault resonance process induced by the seasonal stress changes. This process modulates mid-crustal ramp associated micro-seismicity and influences the timing of Central Himalayan earthquakes.

Tracking the weight of Hurricane Harvey's stormwater using GPS data.

On 26 August 2017, Hurricane Harvey struck the Gulf Coast as a category four cyclone depositing ~95 km3 of water, making it the wettest cyclone in U.S. HISTORY: Water left in Harvey's wake should cause elastic loading and subsidence of Earth's crust, and uplift as it drains into the ocean and evaporates. To track daily changes of transient water storage, we use Global Positioning System (GPS) measurements, finding a clear migration of subsidence (up to 21 mm) and horizontal motion (up to 4 mm) across the Gulf Coast, followed by gradual uplift over a 5-week period. Inversion of these data shows that a third of Harvey's total stormwater was captured on land (25.7 ± 3.0 km3), indicating that the rest drained rapidly into the ocean at a rate of 8.2 km3/day, with the remaining stored water gradually lost over the following 5 weeks at ~1 km3/day, primarily by evapotranspiration. These results indicate that GPS networks can remotely track the spatial extent and daily evolution of terrestrial water storage following transient, extreme precipitation events, with implications for improving operational flood forecasts and understanding the response of drainage systems to large influxes of water.

The rise, collapse, and compaction of Mt. Mantap from the 3 September 2017 North Korean nuclear test.

Surveillance of clandestine nuclear tests relies on a global seismic network, but the potential of spaceborne monitoring has been underexploited. We used satellite radar imagery to determine the complete surface displacement field of up to 3.5 meters of divergent horizontal motion with 0.5 meters of subsidence associated with North Korea's largest underground nuclear test. Combining insight from geodetic and seismological remote sensing, we found that the aftermath of the initial explosive deformation involved subsidence associated with subsurface collapse and aseismic compaction of the damaged rocks of the test site. The explosive yield from the nuclear detonation with best-fitting source parameters for 450-meter depth was 191 kilotonnes of TNT equivalent. Our results demonstrate the capability of spaceborne remote sensing to help characterize large underground nuclear tests.

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.

Seasonal water storage, stress modulation, and California seismicity.

Establishing what controls the timing of earthquakes is fundamental to understanding the nature of the earthquake cycle and critical to determining time-dependent earthquake hazard. Seasonal loading provides a natural laboratory to explore the crustal response to a quantifiable transient force. In California, water storage deforms the crust as snow and water accumulates during the wet winter months. We used 9 years of global positioning system (GPS) vertical deformation time series to constrain models of monthly hydrospheric loading and the resulting stress changes on fault planes of small earthquakes. The seasonal loading analysis reveals earthquakes occurring more frequently during stress conditions that favor earthquake rupture. We infer that California seismicity rates are modestly modulated by natural hydrological loading cycles.

Space geodetic monitoring of engineered structures: The ongoing destabilization of the Mosul dam, Iraq.

We present a detailed survey of the ongoing destabilization process of the Mosul dam. The dam is located on the Tigris river and is the biggest hydraulic structure in Iraq. From a geological point of view the dam foundation is poor due to a site geology formed by alternating strata of highly soluble materials including gypsum, anhydrite, marl and limestone. Here we present the first multi-sensor cumulative deformation map for the dam generated from space-based interferometric synthetic aperture radar measurements from the Italian constellation COSMO-SkyMed and the European sensor Sentinel-1a over the period 2014-2016 that we compare to an older dataset spanning 2004-2010 acquired with the European Envisat satellite. We found that deformation was rapid during 2004-2010, slowed in 2012-2014 and increased since August 2014 when grouting operations stopped due to the temporary capture of the dam by the self proclaimed Islamic State. We model the inferred deformation using a Markov chain Monte Carlo approach to solve for change in volume for simple tensile dislocations. Results from recent and historical geodetic datasets suggests that the volume dissolution rate remains constant when the equivalent volume of total concrete injected during re-grouting operations is included in the calculations.

Asthenosphere rheology inferred from observations of the 2012 Indian Ocean earthquake.

The concept of a weak asthenospheric layer underlying Earth's mobile tectonic plates is fundamental to our understanding of mantle convection and plate tectonics. However, little is known about the mechanical properties of the asthenosphere (the part of the upper mantle below the lithosphere) underlying the oceanic crust, which covers about 60 per cent of Earth's surface. Great earthquakes cause large coseismic crustal deformation in areas hundreds of kilometres away from and below the rupture area. Subsequent relaxation of the earthquake-induced stresses in the viscoelastic upper mantle leads to prolonged postseismic crustal deformation that may last several decades and can be recorded with geodetic methods. The observed postseismic deformation helps us to understand the rheological properties of the upper mantle, but so far such measurements have been limited to continental-plate boundary zones. Here we consider the postseismic deformation of the very large (moment magnitude 8.6) 2012 Indian Ocean earthquake to provide by far the most direct constraint on the structure of oceanic mantle rheology. In the first three years after the Indian Ocean earthquake, 37 continuous Global Navigation Satellite Systems stations in the region underwent horizontal northeastward displacements of up to 17 centimetres in a direction similar to that of the coseismic offsets. However, a few stations close to the rupture area that had experienced subsidence of up to about 4 centimetres during the earthquake rose by nearly 7 centimetres after the earthquake. Our three-dimensional viscoelastic finite-element models of the post-earthquake deformation show that a thin (30-200 kilometres), low-viscosity (having a steady-state Maxwell viscosity of (0.5-10) × 1018 pascal seconds) asthenospheric layer beneath the elastic oceanic lithosphere is required to produce the observed postseismic uplift.

Periodic slow slip triggers megathrust zone earthquakes in northeastern Japan.

Both aseismic and seismic slip accommodate relative motion across partially coupled plate-boundary faults. In northeastern Japan, aseismic slip occurs in the form of decelerating afterslip after large interplate earthquakes and as relatively steady slip on uncoupled areas of the subduction thrust. Here we report on a previously unrecognized quasi-periodic slow-slip behavior that is widespread in the megathrust zone. The repeat intervals of the slow slip range from 1 to 6 years and often coincide with or precede clusters of large [magnitude (M) ≥ 5] earthquakes, including the 2011 M 9 Tohoku-oki earthquake. These results suggest that inherently periodic slow-slip events result in periodic stress perturbations and modulate the occurrence time of larger earthquakes. The periodicity in the slow-slip rate has the potential to help refine time-dependent earthquake forecasts.