skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Detecting Repeating Earthquakes on the San Andreas Fault with Unsupervised Machine Learning of Spectrograms
Abstract Repeating earthquakes—sequences of colocated, quasi-periodic earthquakes of similar size—are widespread along California’s San Andreas fault (SAF) system. Catalogs of repeating earthquakes are vital for studying earthquake source processes, fault properties, and improving seismic hazard models. Here, we introduce an unsupervised machine learning-based method for detecting repeating earthquake sequences (RES) to expand existing RES catalogs or to perform initial, exploratory searches. We implement the “SpecUFEx” algorithm (Holtzman et al., 2018) to reduce earthquake spectrograms into low-dimensional, characteristic fingerprints, and apply hierarchical clustering to group similar fingerprints together independent of location, allowing for a global search for potential RES throughout the data set. We then relocate the potential RES and subject them to the same detection criteria as Waldhauser and Schaff (2021). We apply our method to ∼4000 small (ML 0–3.5) earthquakes located on a 10 km long segment of the creeping SAF and double the number of detected RES, allowing for greater spatial coverage of slip-rate estimations at seismogenic depths. Our method is novel in its ability to detect RES independent of initial locations and is complimentary to existing cross-correlation-based methods, leading to more complete RES catalogs and a better understanding of slip rates at depth.  more » « less
Award ID(s):
2103741
PAR ID:
10552418
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
SSA
Date Published:
Journal Name:
The Seismic Record
Volume:
3
Issue:
4
ISSN:
2694-4006
Page Range / eLocation ID:
376 to 384
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; (, Journal of Geophysical Research: Solid Earth)
    Abstract The Calaveras Fault (CF) branches from the San Andreas Fault (SAF) near San Benito, extending sub‐parallel to the SAF for about 50 km with only 2–6 km separation and diverging northeastward. Both the SAF and CF are partially coupled, exhibit spatially variable aseismic creep and have hosted moderate to large earthquakes in recent decades. Understanding how slip partitions among the main fault strands of the SAF system and establishing their degree of coupling is crucial for seismic hazard evaluation. We perform a timeseries analysis using more than 5 years of Sentinel‐1 data covering the Bay Area (May 2015–October 2020), specifically targeting the spatiotemporal variations of creep rates around the SAF‐CF junction. We derive the surface creep rates from cross‐fault InSAR timeseries differences along the SAF and CF including adjacent Sargent and Quien Sabe Faults. We show that the variable creep rates (0–20 mm/yr) at the SAF‐CF junction are to first order controlled by the angle between the fault strike and the background stress orientation. We further examine the spatiotemporal variation of creep rates along the SAF and CF and find a multi‐annual coupling increase during 2016–2018 the subparallel sections of both faults, with the CF coupling change lagging behind the SAF by 3–6 months. Similar temporal variations are also observed in bothb‐values inferred from declustered seismicity and aseismic slip rates inferred from characteristic repeating earthquakes. 
    more » « less
  2. ; ; ; (, Seismological research letters)
    Allison Bent, Editor-in-Chief (Ed.)
    Seismic data are generally scrutinized for repeating earthquakes (REs) to evaluate slip rates, changes in the mechanical properties of a fault zone, and accelerating nucleation processes in foreshock and aftershock sequences. They are also used to study velocity changes in the medium, earthquake physics and prediction, and for constraining creep rate models at depth. For a robust detection of repeaters, multiple constraints and different parameter configurations related to waveform similarity have been proposed to measure cross‐correlation values at a local seismic network and evaluate the location of overlapping sources. In this work, we developed a Python code to identify REs (FINDRES), inspired by previous literature, which combines both seismic waveform similarity and differential S‐P travel time measured at each seismic station. A cross‐spectral method is applied to evaluate precise differential arrival travel times between earthquake pairs, allowing a subsample precision and increasing the capacity to resolve an overlapping common source radius. FINDRES is versatile and works with and without P‐ and S‐wave phase pickings, and has been validated using synthetic and real data, and provides reliable results. It would contribute to the implementation of open‐source Python packages in seismology, supporting the activities of researchers and the reproducibility of scientific results. 
    more » « less
  3. ; ; ; ; ; (, Nature Communications)
    Abstract Earthquakes occur in clusters or sequences that arise from complex triggering mechanisms, but direct measurement of the slow subsurface slip responsible for delayed triggering is rarely possible. We investigate the origins of complexity and its relationship to heterogeneity using an experimental fault with two dominant seismic asperities. The fault is composed of quartz powder, a material common to natural faults, sandwiched between 760 mm long polymer blocks that deform the way 10 meters of rock would behave. We observe periodic repeating earthquakes that transition into aperiodic and complex sequences of fast and slow events. Neighboring earthquakes communicate via migrating slow slip, which resembles creep fronts observed in numerical simulations and on tectonic faults. Utilizing both local stress measurements and numerical simulations, we observe that the speed and strength of creep fronts are highly sensitive to fault stress levels left behind by previous earthquakes, and may serve as on-fault stress meters. 
    more » « less
  4. ; ; ; (, Nature Communications)
    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
  5. ; (, Journal of Geophysical Research: Solid Earth)
    Abstract Many natural faults are believed to consist of velocity weakening (VW) patches surrounded by velocity strengthening (VS) sections. Numerical studies routinely employ this framework to study earthquake sequences including repeating earthquakes. In this laboratory study, we made a VW asperity, of lengthL, from a bare Poly(methyl methacrylate) PMMA frictional interface and coated the surrounding interface with Teflon to make VS fault sections. Behavior of this isolated asperity was studied as a function ofL(ranging from 100 to 400 mm) and the critical nucleation length, , which is inversely proportional to the applied normal stress (2–16 MPa). Consistent with recent numerical simulations, we observed aseismic slip for  < 2, periodic slip for 2 <  < 6, and non‐periodic slip for 10 < . Furthermore, we compared the experiments whereLwas contained by VS material to standard stick‐slip events whereLwas bounded by free surfaces (i.e.,L = the total sample length). The free surface case produced ∼10 times larger slip during stick‐slip events compared to the contained fault ruptures, even with identical . This disparity highlights how standard, complete‐rupture stick‐slip events differ from contained events expected in nature, due to both the free surface conditions and the heterogeneous normal stress along the fault near the free ends, as confirmed by Digital Image Correlation analysis. This study not only introduces the Teflon coating experimental technique for containing laboratory earthquake ruptures, but also highlights the utility of as a predictive parameter for earthquake behavior. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email