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
  • Dynamic Rupture Models of the 2016 ML 5.8 Gyeongju, South Korea, Earthquake, Constrained by a Kinematic Rupture Model and Seismic Waveform Data
Title: Dynamic Rupture Models of the 2016 ML 5.8 Gyeongju, South Korea, Earthquake, Constrained by a Kinematic Rupture Model and Seismic Waveform Data
ABSTRACT The ML 5.8 earthquake that jolted Gyeongju in southeastern Korea in 2016 was the country’s largest inland event since instrumental seismic monitoring began in 1978. We developed dynamic rupture models of the Gyeongju event constrained by near-source ground-motion data using full 3D spontaneous dynamic rupture modeling with the slip-weakening friction law. Based on our results, we propose two simple dynamic rupture models with constant strength excess (SE) and slip-weakening distance (Dc) that produce near-source ground-motion waveforms compatible with recorded ones in the low-frequency band. Both dynamic models exhibit relatively large stress-drop values, consistent with previous estimates constrained by source spectrum analyses. The fracture energy estimates were also larger than those predicted by a scaling relationship with the seismic moment. The dynamic features constrained in this study by spontaneous rupture modeling and waveform comparison may help understand the source and ground-motion characteristics of future large events in southeastern Korea and thus the seismic hazard of the region.  more » « less
Award ID(s):
2147340
PAR ID:
10511279
Author(s) / Creator(s):
;
Publisher / Repository:
Seismological Society of America
Date Published:
Journal Name:
Bulletin of the Seismological Society of America
Volume:
114
Issue:
2
ISSN:
0037-1106
Page Range / eLocation ID:
710 to 725
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; (, Journal of Geophysical Research: Solid Earth)
    Abstract The Húsavík‐Flatey Fault Zone (HFFZ) is the largest strike‐slip fault in Iceland and poses a high seismic risk to coastal communities. To investigate physics‐based constraints on earthquake hazards, we construct three fault system models of varying geometric complexity and model 79 3‐D multi‐fault dynamic rupture scenarios in the HFFZ. By assuming a simple regional prestress and varying hypocenter locations, we analyze the rupture dynamics, fault interactions, and the associated ground motions up to 2.5 Hz. All models account for regional seismotectonics, topo‐bathymetry, 3‐D subsurface velocity, viscoelastic attenuation, and off‐fault plasticity, and we explore the effect of fault roughness. The rupture scenarios obey earthquake scaling relations and predict magnitudes comparable to those of historical events. We show how fault system geometry and segmentation, hypocenter location, and prestress can affect the potential for rupture cascading, leading to varying slip distributions across different portions of the fault system. Our earthquake scenarios yield spatially heterogeneous near‐field ground motions modulated by geometric complexities, topography, and rupture directivity, particularly in the near‐field. The average ground motion attenuation characteristics of dynamic rupture scenarios of comparable magnitudes and mean stress drop are independent of variations in source complexity, magnitude‐consistent and in good agreement with the latest regional empirical ground motion models. However, physics‐based ground motion variability changes considerably with fault‐distance and increases for unilateral compared to bilateral ruptures. Systematic variations in physics‐based near‐fault ground motions provide important insights into the mechanics and potential earthquake hazard of large strike‐slip fault systems, such as the HFFZ. 
    more » « less
  2. ; ; ; ; (, Geophysical Research Letters)
    Abstract Advances in physics‐based earthquake simulations, utilizing high‐performance computing, have been exploited to better understand the generation and characteristics of the high‐frequency seismic wavefield. However, direct comparison to ground motion observations of a specific earthquake is challenging. We here propose a new approach to simulate data‐fused broadband ground motion synthetics using 3D dynamic rupture modeling of the 2016Mw6.2 Amatrice, Italy earthquake. We augment a smooth, best‐fitting model from Bayesian dynamic rupture source inversion of strong‐motion data (<1 Hz) with fractal fault roughness, frictional heterogeneities, viscoelastic attenuation, and topography. The required consistency to match long periods allows us to quantify the role of small‐scale dynamic source heterogeneities, such as the 3D roughness drag, from observational broadband seismic waveforms. We demonstrate that 3D data‐constrained fully dynamic rupture synthetics show good agreement with various observed ground‐motion metrics up to ∼5 Hz and are an important avenue toward non‐ergodic, physics‐based seismic hazard assessment. 
    more » « less
  3. ; ; (, Journal of Geophysical Research: Solid Earth)
    Abstract Determining conditions for earthquake slip on faults is a key goal of fault mechanics highly relevant to seismic hazard. Previous studies have demonstrated that enhanced dynamic weakening (EDW) can lead to dynamic rupture of faults with much lower shear stress than required for rupture nucleation. We study the stress conditions before earthquake ruptures of different sizes that spontaneously evolve in numerical simulations of earthquake sequences on rate‐and‐state faults with EDW due to thermal pressurization of pore fluids. We find that average shear stress right before dynamic rupture (aka shear prestress) systematically varies with the rupture size. The smallest ruptures have prestress comparable to the local shear stress required for nucleation. Larger ruptures weaken the fault more, propagate over increasingly under‐stressed areas due to dynamic stress concentration, and result in progressively lower average prestress over the entire rupture. The effect is more significant in fault models with more efficient EDW. We find that, as a result, fault models with more efficient weakening produce fewer small events and result in systematically lower b‐values of the frequency‐magnitude event distributions. The findings (a) illustrate that large earthquakes can occur on faults that appear not to be critically stressed compared to stresses required for slip nucleation; (b) highlight the importance of finite‐fault modeling in relating the local friction behavior determined in the lab to the field scale; and (c) suggest that paucity of small events or seismic quiescence may be the observational indication of mature faults that operate under low shear stress due to EDW. 
    more » « less
  4. ; (, Seismological Research Letters)
    Abstract Dynamic rupture simulations generate synthetic waveforms that account for nonlinear source and path complexity. Here, we analyze millions of spatially dense waveforms from 3D dynamic rupture simulations in a novel way to illuminate the spectral fingerprints of earthquake physics. We define a Brune-type equivalent near-field corner frequency (fc) to analyze the spatial variability of ground-motion spectra and unravel their link to source complexity. We first investigate a simple 3D strike-slip setup, including an asperity and a barrier, and illustrate basic relations between source properties and fc variations. Next, we analyze >13,000,000 synthetic near-field strong-motion waveforms generated in three high-resolution dynamic rupture simulations of real earthquakes, the 2019 Mw 7.1 Ridgecrest mainshock, the Mw 6.4 Searles Valley foreshock, and the 1992 Mw 7.3 Landers earthquake. All scenarios consider 3D fault geometries, topography, off-fault plasticity, viscoelastic attenuation, and 3D velocity structure and resolve frequencies up to 1–2 Hz. Our analysis reveals pronounced and localized patterns of elevated fc, specifically in the vertical components. We validate such fc variability with observed near-fault spectra. Using isochrone analysis, we identify the complex dynamic mechanisms that explain rays of elevated fc and cause unexpectedly impulsive, localized, vertical ground motions. Although the high vertical frequencies are also associated with path effects, rupture directivity, and coalescence of multiple rupture fronts, we show that they are dominantly caused by rake-rotated surface-breaking rupture fronts that decelerate due to fault heterogeneities or geometric complexity. Our findings highlight the potential of spatially dense ground-motion observations to further our understanding of earthquake physics directly from near-field data. Observed near-field fc variability may inform on directivity, surface rupture, and slip segmentation. Physics-based models can identify “what to look for,” for example, in the potentially vast amount of near-field large array or distributed acoustic sensing data. 
    more » « less
  5. ; ; ; ; ; ; ; ; (, Nature Communications)
    Abstract Two major earthquakes (MW7.8 and MW7.7) ruptured left-lateral strike-slip faults of the East Anatolian Fault Zone (EAFZ) on February 6, 2023, causing >59,000 fatalities and ~$119B in damage in southeastern Türkiye and northwestern Syria. Here we derived kinematic rupture models for the two events by inverting extensive seismic and geodetic observations using complex 5-6 segment fault models constrained by satellite observations and relocated aftershocks. The larger event nucleated on a splay fault, and then propagated bilaterally ~350 km along the main EAFZ strand. The rupture speed varied from 2.5-4.5 km/s, and peak slip was ~8.1 m. 9-h later, the second event ruptured ~160 km along the curved northern EAFZ strand, with early bilateral supershear rupture velocity (>4 km/s) followed by a slower rupture speed (~3 km/s). Coulomb Failure stress increase imparted by the first event indicates plausible triggering of the doublet aftershock, along with loading of neighboring faults. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email