skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Friday, April 12 until 2:00 AM ET on Saturday, April 13 due to maintenance. We apologize for the inconvenience.

Title: Structural Controls Over the 2019 Ridgecrest Earthquake Sequence Investigated by High‐Fidelity Elastic Models of 3D Velocity Structures

We develop finite element models of the coseismic displacement field accounting for the 3D elastic structures surrounding the epicentral area of the 2019 Ridgecrest earthquake sequence containing two major events of Mw7.1 and Mw6.4. The coseismic slip distribution is inferred from the surface displacement field recorded by interferometric synthetic aperture radar. The rupture dip geometry is further optimized using a novel nonlinear‐crossover‐linear inversion approach. It is found that accounting for elastic heterogeneity and fault along‐strike curvilinearity improves the fit to the observed displacement field and yields a more accurate estimate of geodetic moment and Coulomb stress changes. We observe spatial correlations among the locations of aftershocks and patches of high slip, and rock anomalous elastic properties, suggesting that the shallow crust's elastic structures possibly controlled the Ridgecrest earthquake sequence. Most of the coseismic slip with a peak slip of 7.4 m at 3.6 km depth occurred above a zone of reducedS‐wave velocity and significant post‐Mw7.1 afterslip. This implies that viscous materials or fluid presence might have contributed to the low rupture velocity of the mainshock. Moreover, the zone of high slip on the northwest‐trending fault segment is laterally bounded by two aftershock clusters, whose location is characterized by intermediate rock rigidity. Notably, some minor orthogonal faults consistently end above a subsurface rigid body. Overall, these observations of structural controls improve our understandings of the seismogenesis within incipient fault systems.

more » « less
Award ID(s):
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    On 4 and 6 July 2019, two large strike‐slip earthquakes withW‐phase moment magnitudesMWW6.5 (foreshock) andMWW7.1 (mainshock) struck the Eastern California Shear Zone, northeast of Ridgecrest. The faulting geometry and kinematic coseismic slip distribution of both events are determined by jointly inverting seismological and geodetic observations guided by aftershock and surface rupture locations. The foreshock ruptured two orthogonal faults with a prominent L‐shaped geometry with maximum slip of ~1.1 m on the NE‐SW segment. The mainshock faulting extended NW‐SE along several primary fault segments that straddle the foreshock slip. The surface rupture and slip model indicate mostly near‐horizontal strike‐slip motion with maximum slip of ~3.7 m, but there is a localized vertical dip‐slip motion. Both the foreshock and mainshock ruptures terminate in regions of complex surface offsets. High aftershock productivity and low rupture velocity may be the result of rupture of a relatively immature fault system.

    more » « less
  2. null (Ed.)
    ABSTRACT The 2019 Ridgecrest earthquake sequence culminated in the largest seismic event in California since the 1999 Mw 7.1 Hector Mine earthquake. Here, we combine geodetic and seismic data to study the rupture process of both the 4 July Mw 6.4 foreshock and the 6 July Mw 7.1 mainshock. The results show that the Mw 6.4 foreshock rupture started on a northwest-striking right-lateral fault, and then continued on a southwest-striking fault with mainly left-lateral slip. Although most moment release during the Mw 6.4 foreshock was along the southwest-striking fault, slip on the northwest-striking fault seems to have played a more important role in triggering the Mw 7.1 mainshock that happened ∼34  hr later. Rupture of the Mw 7.1 mainshock was characterized by dominantly right-lateral slip on a series of overall northwest-striking fault strands, including the one that had already been activated during the nucleation of the Mw 6.4 foreshock. The maximum slip of the 2019 Ridgecrest earthquake was ∼5  m, located at a depth range of 3–8 km near the Mw 7.1 epicenter, corresponding to a shallow slip deficit of ∼20%–30%. Both the foreshock and mainshock had a relatively low-rupture velocity of ∼2  km/s, which is possibly related to the geometric complexity and immaturity of the eastern California shear zone faults. The 2019 Ridgecrest earthquake produced significant stress perturbations on nearby fault networks, especially along the Garlock fault segment immediately southwest of the 2019 Ridgecrest rupture, in which the coulomb stress increase was up to ∼0.5  MPa. Despite the good coverage of both geodetic and seismic observations, published coseismic slip models of the 2019 Ridgecrest earthquake sequence show large variations, which highlight the uncertainty of routinely performed earthquake rupture inversions and their interpretation for underlying rupture processes. 
    more » « less
  3. ABSTRACT The July 2019 Ridgecrest, California, earthquake sequence involved two large events—the M 6.4 foreshock and the M 7.1 mainshock that ruptured a system of intersecting strike-slip faults. We present analysis of space geodetic observations including Synthetic Aperture Radar and Global Navigation Satellite System data, geological field mapping, and seismicity to constrain the subsurface rupture geometry and slip distribution. The data render a complex pattern of faulting with a number of subparallel as well as cross-cutting fault strands that exhibit variations in both strike and dip angles, including a “flower structure” formed by shallow splay faults. Slip inversions are performed using both homogeneous and layered elastic half-space models informed by the local seismic tomography data. The inferred slip distribution suggests a moderate amount of the shallow coseismic slip deficit. The peak moment release occurred in the depth interval of 3–4 km, consistent with results from previous studies of major strike-slip earthquakes, and the depth distribution of seismicity in California. We use the derived slip models to investigate stress transfer and possible triggering relationships between the M 7.1 mainshock and the M 6.4 foreshock, as well as other moderate events that occurred in the vicinity of the M 7.1 hypocenter. Triggering is discouraged for the average strike of the M 7.1 rupture (320°) but encouraged for the initial orientation of the mainshock rupture suggested by the first-motion data (340°). This lends support to a scenario according to which the earthquake rupture nucleated on a small fault that was more optimally oriented with respect to the regional stress and subsequently propagated along the less-favorably oriented pre-existing faults, possibly facilitated by dynamic weakening. The nucleation site of the mainshock experienced positive dynamic Coulomb stress changes that are much larger than the static stress changes, yet the former failed to initiate rupture. 
    more » « less
  4. Abstract

    Fault zone structures at many scales largely dictate earthquake ruptures and are controlled by the geologic setting and slip history. Characterizations of these structures at diverse scales inform better understandings of earthquake hazards and earthquake phenomenology. However, characterizing fault zones at sub‐kilometer scales has historically been challenging, and these challenges are exacerbated in urban areas, where locating and characterizing faults is critical for hazard assessment. We present a new procedure for characterizing fault zones at sub‐kilometer scales using distributed acoustic sensing (DAS). This technique involves the backprojection of the DAS‐measured scattered wavefield generated by natural earthquakes. This framework provides a measure of the strength of scattering along a DAS array and thus constrains the positions and properties of local scatterers. The high spatial sampling of DAS arrays makes possible the resolution of these scatterers at the scale of tens of meters over distances of kilometers. We test this methodology using a DAS array in Ridgecrest, CA which recorded much of the 2019 Mw7.1 Ridgecrest earthquake aftershock sequence. We show that peaks in scattering along the DAS array are spatially correlated with mapped faults in the region and that the strength of scattering is frequency‐dependent. We present a model of these scatterers as shallow, low‐velocity zones that is consistent with how we may expect faults to perturb the local velocity structure. We show that the fault zone geometry can be constrained by comparing our observations with synthetic tests.

    more » « less
  5. Abstract

    We analyze seismograms recorded by four arrays (B1–B4) with 100 m station spacing and apertures of 4–8 km that cross the surface rupture of the 2019 Mw 7.1 Ridgecrest earthquake. The arrays extend from B1 in the northwest to B4 in the southeast of the surface rupture. Delay times betweenPwave arrivals associated with ∼1,200 local earthquakes and four teleseismic events are used to estimate local velocity variations beneath the arrays. Both teleseismic and localPwaves travel faster on the northeast than the southwest side of the fault beneath arrays B1 and B4, but the velocity contrast is less reliably resolved at arrays B2 and B3. We identify several 1–2 km wide low‐velocity zones with much slower inner cores that amplifySwaveforms, inferred as damage zones, beneath each array. The damage zones at arrays B2 and B4 also generate fault‐zone head and trapped waves. An automated detector, based on peak ground velocities and durations of high‐amplitude waves, identifies candidate fault‐zone trapped waves (FZTWs) in a localized zone for ∼600 earthquakes at array B4. Synthetic waveform modeling of averaged FZTWs, generated by ∼30 events with high‐quality signals, indicates that the trapping structure at array B4 has a width of ∼300 m, depth of 3–5 km,Swave velocity reduction of ∼20% with respect to the surrounding rock,Q‐value of ∼30, andSwave velocity contrast of ∼4% across the fault (faster on the northeast side). The results show complex fault‐zone internal structures (velocity contrasts and low‐velocity zones) that vary along fault strike.

    more » « less