skip to main content


Title: Slow‐Growing and Extended‐Duration Seismicity Swarms: Reactivating Joints or Foliations in the Cahuilla Valley Pluton, Central Peninsular Ranges, Southern California
Abstract

Three prolific earthquake swarms and numerous smaller ones have occurred since 1980 in the Mesozoic igneous plutonic rocks of the Perris block of the Peninsular Ranges, Southern California. The major swarms occurred in 1980–1981, 1983–1984, and 2016–2018, with the latest swarm still ongoing. These swarms have no clear mainshock, with the largest events ofML3.6,ML3.7, andMw4.4. Each successive swarm had larger cumulative seismic moment release with about 314 and 411 events ofM ≥ 1.5, while the third swarm has produced about 451 events ofM ≥ 1.5 (as of September 2018). The concurrent strike‐slip faulting occurred on north to northwest striking planes but with no orthogonal northeast trending seismicity alignments. These shallow swarms are probably driven by intrablock Pacific‐North America plate boundary stress loading of the two bounding major late Quaternary strike‐slip faults, the Elsinore and San Jacinto faults. The state of stress within the Cahuilla Valley pluton has a ~40° angle between the maximum principal stress and the average trend of the swarms, suggesting that migrating pore fluid pressures aid in the formation and growth of zones of weakness. These swarms, which last more than 600 days each, exhibit clear bilateral spatial migration for distances of up to ~7–8 km and reach their full length in about 20 months. The slow spatial‐temporal development of the swarms corresponds to a fluid diffusivity of 0.006 to 0.01 m2/s, consistent with very low permeability rocks as expected for this block. There is no geodetic or other evidence for a slow slip event driving the swarms.

 
more » « less
NSF-PAR ID:
10453259
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
124
Issue:
4
ISSN:
2169-9313
Page Range / eLocation ID:
p. 3933-3949
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    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
  2. Abstract

    We conduct a detailed study of the foreshock sequence preceding the 2010Mw6.7 Yushu, Qinghai earthquake in the Tibetan plateau by examining continuous waveforms recorded at a seismic station near the mainshock rupture zone. By using a deep learning phase picker—EQTransformer and a matched‐filter technique, we identify 120 foreshocks with magnitude ranging from −0.7 to 1.6, starting with aMw4.6 foreshock approximately 2 hr before theMw6.7 Yushu mainshock. Our analyses show that the foreshock sequence follows a typical Omori's law decay with ap‐value of 0.73 and the Gutenberg‐Richer frequency‐magnitudeb‐value of 0.66. We do not find any evidence of accelerating events leading up to the Yushu mainshock. Hence, they could be considered as aftershocks of theMw4.6 earthquake. We further invert for the focal mechanisms and rupture directions for both the largest foreshock and the mainshock. TheMw4.6 foreshock likely occurred on a NE‐SW trending fault conjugating to the NW‐SE trending fault of the mainshock. Coulomb stress analysis suggests theMw4.6 foreshock induces negative stress on the mainshock source area. These observations do not support either the pre‐slip or the cascade triggering model for foreshock generation. The occurrence of the foreshock, mainshock and large aftershocks appear to be modulated by the Earth's tidal forces, likely reflecting the role of high pore‐fluid pressures. Our observations, together with other recent studies, suggest that extensional step‐overs and conjugate faults along major strike‐slip faults play an important role in generating short‐term foreshock sequences.

     
    more » « less
  3. Abstract

    Fault-zone fluids control effective normal stress and fault strength. While most earthquake models assume a fixed pore fluid pressure distribution, geologists have documented fault valving behavior, that is, cyclic changes in pressure and unsteady fluid migration along faults. Here we quantify fault valving through 2-D antiplane shear simulations of earthquake sequences on a strike-slip fault with rate-and-state friction, upward Darcy flow along a permeable fault zone, and permeability evolution. Fluid overpressure develops during the interseismic period, when healing/sealing reduces fault permeability, and is released after earthquakes enhance permeability. Coupling between fluid flow, permeability and pressure evolution, and slip produces fluid-driven aseismic slip near the base of the seismogenic zone and earthquake swarms within the seismogenic zone, as ascending fluids pressurize and weaken the fault. This model might explain observations of late interseismic fault unlocking, slow slip and creep transients, swarm seismicity, and rapid pressure/stress transmission in induced seismicity sequences.

     
    more » « less
  4. Cooling ages of tectonic blocks between the Yakutat microplate and the Fairweather transform boundary fault reveal exhumation due to strike-slip faulting and subsequent collision into this tectonic corner. The Yakutat and Boundary faults are splay faults that define tectonic panels with bounding faults that have evidence of both reverse and strike-slip motion, and they are parallel to the northern end of the Fairweather fault. Uplift and exhumation simultaneous with strike-slip motion have been significant since the late Miocene. The blocks are part of an actively deforming tectonic corner, as indicated by the ~14–1.5 m of coseismic uplift from the M 8.1 Yakutat Bay earthquake of 1899 and 4 m of strike-slip motion in the M 7.9 Lituya Bay earthquake in 1958 along the Fairweather fault. New apatite (U-Th-Sm)/He (AHe) and zircon (U-Th)/He (ZHe) data reveal that the Boundary block and the Russell Fiord block have different cooling histories since the Miocene, and thus the Boundary fault that separates them is an important tectonic boundary. Upper Cretaceous to Paleocene flysch of the Russell Fiord block experienced a thermal event at 50 Ma, then a relatively long period of burial until the late Miocene when initial exhumation resulted in ZHe ages between 7 and 3 Ma, and then very rapid exhumation in the last 1–1.5 m.y. Exhumation of the Russell Fiord block was accommodated by reverse faulting along the Yakutat fault and the newly proposed Calahonda fault, which is parallel to the Yakutat fault. The Eocene schist of Nunatak Fiord and 54–53 Ma Mount Stamy and Mount Draper granites in the Boundary block have AHe and ZHe cooling ages that indicate distinct and very rapid cooling between ca. 5 Ma and ca. 2 Ma. Rocks of the Chugach Metamorphic Complex to the northeast of the Fairweather fault and in the fault zone were brought up from 10–12 km at extremely high rates (>5 km/m.y.) since ca. 3 Ma, which implies a significant component of dip-slip motion along the Fairweather fault. The adjacent rocks of the Boundary block were exhumed with similar rates and from similar depths during the early Pliocene, when they may have been located 220–250 km farther south near Baranof Island. The profound and significant exhumation of the three tectonic blocks in the last 5 m.y. has probably been driven by uplift and erosional exhumation due to contraction as rocks collide into this tectonic corner. The documented spatial and temporal pattern of exhumation is in agreement with the southward shift of focused exhumation at the St. Elias syntaxial corner and the southeast propagation of the fold-and thrust belt. 
    more » « less
  5. Abstract

    Seamounts are found at many subduction zones and act as seafloor heterogeneities that affect slip behavior on megathrusts. At the Hikurangi subduction zone offshore the North Island, New Zealand, seamounts have been identified on the incoming Pacific plate and below the accretionary prism, but there is little concrete evidence for seamounts subducted beyond the present‐day coastline. Using a high‐resolution, adjoint tomography‐derived velocity model of the North Island, we identify two high‐velocity anomalies below the East Coast and an intraslab low‐velocity zone up‐dip of one of these anomalies. We interpret the high‐velocity anomalies as previously unidentified, deeply subducted seamounts, and the low‐velocity zone as fluid in the subducting slab. The seamounts are inferred to be 10–30 km wide and on the plate interface at 12–15 km depth. Resolution analysis using point spread functions confirms that these are well‐resolved features. The locations of the two seamounts coincide with bathymetric features whose geometries are consistent with those predicted from analog experiments and numerical simulations of seamount subduction. The spatial characteristics of seismicity and slow slip events near the inferred seamounts agree well with previous numerical modeling predictions of the effects of seamount subduction on megathrust stress and slip. Anomalous geophysical signatures, magnetic anomalies, and swarm seismicity have also been observed previously at one or both seamount locations. We propose that permanent fracturing of the northern Hikurangi upper plate by repeated seamount subduction may be responsible for the dichotomous slow slip behavior observed geodetically, and partly responsible for along‐strike variations in plate coupling on the Hikurangi subduction interface.

     
    more » « less