Ground Motion Characteristics of Cascading Earthquakes in a Multiscale Fracture Network

TL;DR

This study uses high-resolution simulations to analyze how cascading ruptures in complex fracture networks influence ground motion characteristics, revealing significant effects on seismic signals and hazard assessments.

Contribution

It introduces detailed 3D dynamic rupture models of cascading earthquakes in multiscale fracture networks, highlighting their impact on ground motion signatures and seismic hazard evaluation.

Findings

Cascading ruptures produce high-frequency, source coda-wave-like signatures.

Elevated characteristic frequency and stress drop in cascading events.

Pronounced azimuthal variations in ground motion parameters.

Abstract

Fault zones exhibit geometrical complexity and are often surrounded by multiscale fracture networks within their damage zones, influencing rupture dynamics and near-field ground motions. We investigate the ground-motion characteristics of cascading ruptures across damage zone fracture networks of moderate-sized earthquakes using high-resolution 3D dynamic rupture simulations. Our models feature a listric fault surrounded by over 800 fractures, emulating a major fault and its associated damage zone. We analyze three cases: a cascading rupture propagating within the fracture network, a non-cascading main-fault rupture with off-fault fracture slip, and a main-fault rupture without a fracture network. Cascading ruptures within the fracture network produce distinct ground-motion signatures with 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 components reflects the complexity of the radiation pattern and rupture directivity of 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 cascading ruptures also show pronounced azimuthal variations in peak ground acceleration (PGA), peak ground velocity, and pseudo-spectral acceleration, with average PGA nearly double that of the non-cascading cases. 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.