Search

Mature strike‐slip faults are usually surrounded by a narrow zone of damaged rocks characterized by low seismic wave velocities. Observations of earthquakes along such faults indicate that seismicity is highly concentrated within this fault damage zone. However, the long‐term influence of the fault damage zone on complete earthquake cycles, that is, years to centuries, is not well understood. We simulate aseismic slip and dynamic earthquake rupture on a vertical strike‐slip fault surrounded by a fault damage zone for a thousand‐year timescale using fault zone material properties and geometries motivated by observations along major strike‐slip faults. The fault damage zone is approximated asan elastic layer with lower shear wave velocity than the surrounding rock. We find that dynamic wave reflections, whose characteristics are strongly dependent on the width and the rigidity contrast of the fault damage zone, have a prominent effect on the stressing history of the fault. The presence of elastic damage can partially explain the variability in the earthquake sizes and hypocenter locations along a single fault, which vary with fault damage zone depth, width and rigidity contrast from the host rock. The depth extent of the fault damage zone has a pronounced effect on the earthquake hypocenter locations, and shallower fault damage zones favor shallower hypocenters with a bimodal distribution of seismicity along depth. Our findings also suggest significant effects on the hypocenter distribution when the fault damage zone penetrates to the nucleation sites of earthquakes, likely being influenced by both lithological (material) and rheological (frictional) boundaries. Plain Language Summary: Large strike‐slip earthquakes tend to create a zone of fractured network surrounding the main fault. This zone, referred to as a fault damage zone, becomes highly localized as the fault matures, with a width of few hundred meters. The influence of this fault damage zone on earthquake characteristics remains elusive since we do not have enough long‐term observations along a single fault. We use numerical simulations to examine the behavior of earthquake nucleation and rupture dynamics on a fault surrounded by a damage zone over a thousand‐year timescale. Our simulations reveal that the reflection of seismic waves from the fault damage zone boundaries leads to complexity in earthquake sequences, such as variability in earthquake locations and sizes. We also show that a shallowfault damage zone produces shallower earthquakes with the earthquake depths centered around two locations (bimodal), as opposed to a deep fault damage zone with the earthquake depths centered around a single location (unimodal). Our study suggests that imaging the geometry and physical properties of fault damage zones could potentially give us clues about depths of future earthquakes and improve earthquake probabilistic hazard assessment. Key Points: Fully dynamic earthquake cycle simulations show persistent heterogeneous stress distribution generated by fault zone wavesFaults surrounded by low‐velocity damage zones lead to more complexities in earthquakelocation, size, and slip patternsBoth lithology and rheology influence the depth distribution of seismicity, with shallow fault damage zones exhibiting bimodal distribution [ABSTRACT FROM AUTHOR]

There is a strong need to model potential rupture behaviors for the next Cascadia megathrust earthquake. However, there exists significant uncertainty regarding the extent of downdip rupture and rupture speed. To address this problem, we study how the transition region (i.e., the gap), which separates the locked from slow‐slip regions, influences coseismic rupture propagation using 2‐D dynamic rupture simulations governed by a slip‐weakening friction law. We show that rupture propagation through the gap is strongly controlled by the amount of accumulated tectonic initial shear stress and gap friction level. A large amplitude negative dynamic stress drop is needed to arrest downdip rupture. We also observe downdip supershear rupture when the gradient in effective normal stress from the locked to slow‐slip regions is dramatic. Our results justify kinematic rupture models that extend below the gap and suggests the possibility of high‐frequency energy radiation during the next Cascadia megathrust earthquake. Plain Language Summary: How large, deep, and damaging a future earthquake will be depends on factors such as energy release that must be constrained by precise observations of previous earthquakes in the same area. But such data are rarely available. Instead, computer models of earthquakes guided by the laws of physics can provide us with estimates of potential ground shaking for a future event. In our study, we design two‐dimensional earthquake simulations for the Cascadia fault below the northwestern United States coast and test different hypotheses for how stress may be accumulating at depth along this fault. Our models focus on a portion of the fault referred to as the "gap." The gap physically separates a shallow region that slips during large earthquakes from a deeper region that experiences intermittent slip between large earthquakes. A gap region similar to that in Cascadia is also found in Japan, Mexico, and around other active faults worldwide. We find that our simulated rupture is able to extend to deeper regions at faster speeds given the current understanding of stress levels and earthquake fault friction in the gap. While this work represents only a first step toward understanding how stresses and friction influence how the Cascadia fault might slip, it lays the foundation for modeling more complex physics that can help scientists better predict shaking from seismic waves. Key Points: We examine dynamic source effects on along‐dip rupture propagation for a Cascadia megathrust earthquakeSimulated earthquake rupture is able to penetrate through the transition zone and reach the deeper slow‐slip regionOur results underscore the potential for a deeper downdip rupture and faster rupture speed than previously assumed in kinematic models [ABSTRACT FROM AUTHOR]

Geological and geophysical observations reveal along‐strike fault zone heterogeneity on major strike‐slip faults, which can play a significant role in earthquake rupture propagation and termination. I present 2‐D dynamic rupture simulations to demonstrate rupture characteristics in such heterogeneous fault zone structure. The modeled rupture is nucleated in a damaged fault zone and propagates on a preexisting fault toward the zone of intact rocks. There is an intermediate range of nucleation lengths that only allow rupture to spontaneously propagate in the damaged fault zone but not in a homogeneous medium given the same stresses and frictional parameters. Rupture with an intermediate nucleation length tends to stop when it reaches the zone of intact rocks for uniform fault stress conditions, especially when the rupture propagation distance in the damaged fault zone is relatively short and when the damaged fault zone is relatively narrow or smooth in the fault‐normal direction. Pronounced small‐scale heterogeneity within the damaged fault zone also contributes to such early rupture termination. In asymmetric fault zones bisected by a bimaterial fault, rupture moving in the direction of slip of faster rocks tends to terminate under the same conditions as in symmetric fault zones, whereas rupture moving in the direction of slip of slower rocks can penetrate into the zone of intact rocks. A sufficiently large asperity located at the edge of the zone of intact rocks also allows break‐through rupture. The results suggest that the along‐strike fault zone heterogeneity can play a critical role in seismicity distribution. Plain Language Summary: Natural faults are surrounded by a zone of deformed rocks to accommodate strain localization. Such fault zone is not continuous along the fault, but rather includes segments of relatively damaged rocks adjacent to segments of relatively intact rocks. By simulating the dynamic interactions between fault stress, friction, and fault zone heterogeneities during the earthquake rupture process, I show that rupture is more likely to spontaneously propagate inside the damaged fault zone and stop when it reaches the relatively intact zone for uniform fault stress conditions. This phenomenon is less pronounced when the damaged fault zone becomes wider, sharper, and more damaged, indicating a higher likelihood of having large earthquakes that can penetrate into the relatively intact zone on more mature faults. The results suggest that a priori knowledge of the fault zone heterogeneity is critical for understanding the spatial distribution of earthquakes and the likelihood of having large earthquakes. Key Points: Rupture nucleated in a damaged fault zone tends to terminate when it propagates along strike to an intact zone for uniform fault stressesRupture tends to penetrate into the relatively intact zone when the damaged fault zone becomes wider, sharper, and more damagedAn asperity at the edge of the relatively intact zone can facilitate rupture when its size is comparable to the nucleation half‐length [ABSTRACT FROM AUTHOR]