Multicycle Simulation of Strike-Slip Earthquake Rupture for Use in Near-Source Ground-Motion Simulations

Abstract

Realistic dynamic rupture modeling validated by observed earthquakes is necessary for estimating parameters that are poorly resolved by seismic source inversion, such as stress drop, rupture velocity, and slip rate function. Source inversions using forward dynamic modeling are increasingly used to obtain earthquake rupture models. In this study, to generate a large number of physically self-consistent rupture models, rupture process of which is consistent with the spatiotemporal heterogeneity of stress produced by previous earthquakes on the same fault, we use multicycle simulations under the rate and state (RS) friction law. We adopt a one-way coupling from multicycle simulations to dynamic rupture simulations; the quasidynamic solver QDYN is used to nucleate the seismic events and the spectral element dynamic solver SPECFEM3D to resolve their rupture process. To simulate realistic seismicity, with a wide range of magnitudes and irregular recurrence, several realizations of 2D-correlated heterogeneous random distributions of characteristic weakening distance (Dc) in RS friction are tested. Other important parameters are the normal stress, which controls the stress drop and rupture velocity during an earthquake, and the maximum value of Dc, which controls rupture velocity but not stress drop. We perform a parametric study on a vertical planar fault and generate a set of a hundred spontaneous rupture models in a wide magnitude range (Mw 5.5–7.4). We validate the rupture models by comparison of source scaling, ground motion (GM), and surface slip properties to observations. We compare the source-scaling relations between rupture area, average slip, and seismic moment of the modeled events with empirical ones derived from source inversions. Near-fault GMs are computed from the source models. Their peak ground velocities and peak ground accelerations agree well with the ground-motion prediction equation values. We also obtain good agreement of the surface fault displacements with observed values.

A numerical study of 3D topographic site effects considering wavefield incident direction and geomorphometric parameters

“…However, with the growth in computational capabilities, the development of codes capable of handling complex topographies, and more high-quality near-surface geologic data available, it is becoming easier to directly model groundmotion amplification due to topography. Some studies that consider surface topography has obtained simulation results that are in good agreement with observed ground motions (Magnoni et al, 2013;Galvez et al, 2021;Wang et al, 2021).…”

Rupture Dynamics of Cascading Earthquakes in a Multiscale Fracture Network

“…Dynamic rupture modeling across different fault or stress‐heterogeneity scales may use variable characteristic slip distances (Ando & Yamashita, 2007; Bizzarri & Cocco, 2003; Gallovič & Valentová, 2023; Galvez et al., 2021; Ide & Aochi, 2005; Ulrich et al., 2022). The critical nucleation size required to initiate spontaneous or runaway rupture (i.e., ruptures across the entire main fault (Galis et al., 2019)), scale with characteristic slip distance L under rate‐and‐state friction (Rubin & Ampuero, 2005), or with critical slip distance D c under linear slip‐weakening friction (Andrews, 1976; Day, 1982; Galis et al., 2015).…”

A quantitative comparison and validation of finite-fault models: The 2011 Tohoku-Oki earthquake

“…We note that our analysis does not take highly complex rupture propagation effects into account such as, for example, multiple slip-episodes inferred from multiple timewindow slip inversion or in dynamic rupture scenario simulations informed from local strong ground motions (Galvez et al, 2016(Galvez et al, , 2020.…”