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.


  • NSF Public Access
  • Search Results
  • Spectral Characteristics of Hydraulic Fracturing-Induced Seismicity Can Distinguish between Activation of Faults and Fractures
Title: Spectral Characteristics of Hydraulic Fracturing-Induced Seismicity Can Distinguish between Activation of Faults and Fractures
Abstract Analysis of earthquake spectra can aid in understanding source characteristics like stress drop and rupture complexity. There is growing interest in probing the similarities and differences of fault rupture for natural and human-induced seismic events. Here, we analyze waveform data from a shallow, buried geophone array that recorded seismicity during a hydraulic fracturing operation near Fox Creek, Alberta. Starting from a quality-controlled catalog of 4000 events between magnitude 0 and 3.2, we estimate source-spectral corner frequencies using methods that account for the band-limited nature of the sensor response. The stress-drop values are found to be approximately self-similar, but with a slight magnitude dependence in which larger events have higher stress drop (∼10 MPa). Careful analysis of the relative corner frequencies shows that individual fault and fracture segments experienced systematic variations in relative corner frequency over time, indicating a possible change in the stress state. Clustering analysis of source spectra based on the relative proportion of high- and low-frequency content relative to the Brune model further shows that event complexity evolves over time. In addition, the faults produce earthquakes with systematically larger stress-drop values than the fractures. Combined, these results indicate that the features activated by hydraulic fracturing experience observable changes in source behavior over time and exhibit different properties depending on the orientation, scale, and fabric of the structural feature on which they occur.  more » « less
Award ID(s):
2231705
PAR ID:
10410962
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Seismological Society of America
Date Published:
Journal Name:
Seismological Research Letters
ISSN:
0895-0695
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; ; ; ; ; et al (, Bulletin of the Seismological Society of America)
    ABSTRACT We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress-drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude-spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress-drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two-week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. In this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite-fault modeling, ground-motion-based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between-method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress-drop values. 
    more » « less
  2. ; ; (, Bulletin of the Seismological Society of America)
    ABSTRACT Although the Brune source model describes earthquake moment release as a single pulse, it is widely used in studies of complex earthquakes with multiple episodes of high moment release (i.e., multiple subevents). In this study, we investigate how corner frequency estimates of earthquakes with multiple subevents are biased if they are based on the Brune source model. By assuming complex sources as a sum of multiple Brune sources, we analyze 1640 source time functions of Mw 5.5–8.0 earthquakes in the seismic source characteristic retrieved from deconvolving teleseismic body waves catalog to estimate the corner frequencies, onset times, and seismic moments of subevents. We identify more subevents for strike-slip earthquakes than dip-slip earthquakes, and the number of resolvable subevents increases with magnitude. We find that earthquake corner frequency correlates best with the corner frequency of the subevent with the highest moment release (i.e., the largest subsevent). This suggests that, when the Brune model is used, the estimated corner frequency and, therefore, the stress drop of a complex earthquake is determined primarily by the largest subevent rather than the total rupture area. Our results imply that, in addition to the simplified assumption of a radial rupture area with a constant rupture velocity, the stress variation of asperities, rather than the average stress change of the whole fault, contributes to the large variance of stress-drop estimates. 
    more » « less
  3. ; ; (, SSA Annual Meeting)
    The 2019 Ridgecrest, CA earthquake sequence has provided a unique opportunity and a rich dataset to understand earthquake source properties and near-fault structure. Using the high-quality seismic data provided by the SCEC Stress Drop Validation group, we first estimate the corner frequency of M2.0-4.5 earthquakes by applying the spectral ratio method based on empirical Green’s function (Liu et al., 2020). We relate corner frequency estimates to stress drops assuming the Brune source model and circular cracks. Our preliminary results show increasing median stress drops with magnitude for both P and S waves, from 1 MPa for M2.0 events to 10 MPa for M4.0 events, though the limited frequency bandwidth may cause underestimation for small events. The estimated moment magnitude is proportional to the catalog magnitude by a factor of 0.72, which is close to 0.74 estimated by Trugman (2020) for the Ridgecrest earthquake sequence. In the second part of the study, we examine the impact of fault zone structure on the azimuthal variation of the source spectra. Using kinematic simulations and observations of the 2003 Big Bear earthquake sequence, Huang et al. (2016) showed that fault damage zones can act as an effective wave guide and cause high-frequency wave amplification along directions close to fault strike. We use clusters of M1.5-3 earthquakes in the Ridgecrest region to further examine the azimuthal variation of the stacked source spectra and investigate if the near-source structure can affect our corner frequency estimates. We aim to develop robust methods that utilize high-quality seismic data to illuminate earthquake source processes and fault zone properties. 
    more » « less
  4. ; ; ; ; ; (, Bulletin of the Seismological Society of America)
    ABSTRACT Fault zones exhibit geometrical complexity and are often surrounded by multiscale fracture networks within their damage zones, potentially influencing rupture dynamics and near-field ground motions. In this study, we investigate the ground-motion characteristics of cascading ruptures across damage zone fracture networks of moderate-size earthquakes (Mw 5.5–6.0) using high-resolution 3D dynamic rupture simulations. Our models feature a listric normal fault surrounded by more than 800 fractures, emulating a major fault and its associated damage zone. We analyze three cases: a cascading rupture propagating within the fracture network (Mw 5.5), a non-cascading main-fault rupture with off-fault fracture slip (Mw 6.0), and a main-fault rupture without a fracture network (Mw 6.0). Cascading ruptures within the fracture network produce distinct ground-motion signatures with enriched high-frequency content, arising from simultaneous slip of multiple fractures and parts of the main fault, resembling source coda-wave-like signatures. This case shows elevated near-field characteristic frequency (fc) and stress drop, approximately an order of magnitude higher than the estimation directly on the fault of the dynamic rupture simulation. The inferred fc of the modeled vertical ground-motion components reflects the complexity of the radiation pattern and rupture directivity of fracture-network cascading earthquakes. We show that this is consistent with observations of strong azimuthal dependence of corner frequency in the 2009–2016 central Apennines, Italy, earthquake, sequence. Simulated ground motions from fracture-network cascading ruptures also show pronounced azimuthal variations in peak ground acceleration (PGA), peak ground velocity, and pseudospectral acceleration, with average PGA nearly double that of the non-cascading cases. Cascading ruptures radiate high-frequency seismic energy, yield nontypical ground-motion characteristics including coda-wave-like signatures, and may result in a significantly higher seismologically inferred stress drop and PGA. Such outcomes emphasize the critical role of fault-zone complexity in affecting rupture dynamics and seismic radiation and have important implications for physics-based seismic hazard assessment. 
    more » « less
  5. ; ; (, Bulletin of the Seismological Society of America)
    ABSTRACT Earthquake stress drop—a key parameter for describing the energetics of earthquake rupture—can be estimated in several different, but theoretically equivalent, ways. However, independent estimates for the same earthquakes sometimes differ significantly. We find that earthquake source complexity plays a significant role in why theoretically (for simple rupture models) equivalent methods produce different estimates. We apply time- and frequency-domain methods to estimate stress drops for real earthquakes in the SCARDEC (Seismic source ChAracteristics Retrieved from DEConvolving teleseismic body waves, Vallée and Douet, 2016) source time function (STF) database and analyze how rupture complexity drives stress-drop estimate discrepancies. Specifically, we identify two complexity metrics—Brune relative energy (BRE) and spectral decay—that parameterize an earthquake’s complexity relative to the standard Brune model and strongly correlate with the estimate discrepancies. We find that the observed systematic magnitude–stress-drop trends may reflect underlying changes in STF complexity, not necessarily trends in actual stress drop. Both the decay and BRE parameters vary systematically with magnitude, but whether this magnitude–complexity relationship is real remains unresolved. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email