Search

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.

Fault-damage zones comprise multiscale fracture networks that may slip dynamically and interact with the main fault during earthquake rupture. Using 3D dynamic rupture simulations and scale-dependent fracture energy, we examine dynamic interactions of more than 800 intersecting multiscale fractures surrounding a listric fault, emulating a major fault and its damage zone. We investigate ten distinct orientations of maximum horizontal stress, probing the conditions necessary for sustained slip within the fracture network or activating the main fault. Additionally, we assess the feasibility of nucleating dynamic rupture earthquake cascades from a distant fracture and investigate the sensitivity of fracture network cascading rupture to the effective normal stress level. We model either pure cascades or main fault rupture with limited off-fault slip. We find that cascading ruptures within the fracture network are dynamically feasible under certain conditions, including: (i) the state-evolutional distance scales with fracture and fault size, (ii) favorable relative pre-stress of fractures within the ambient stress field, and (iii) close proximity of fractures. We find that cascading rupture within the fracture network discourages rupture on the main fault. Our simulations suggest that favorable relative pre-stress fractures within a fault damage zone may lead to cascading earthquake rupture reaching off-fault moment magnitudes up to $Mw \approx 5.6$, shadowing the main fault slip. Our findings offer fundamental insights into physical processes governing cascading earthquake dynamic rupture within multiscale fracture networks. Our results have implications for the seismic hazard of naturally activated fracture or fault networks and earthquakes induced in geo-energy exploitation activities.