How to Promote Earthquake Ruptures: Different Nucleation Strategies in a Dynamic Model with Slip-Weakening Friction

Abstract

The introduction of the linear slip-weakening friction law permits the solution of the elastodynamic equation for a rupture that develops on a fault by removing the singularity in the components of stress tensor, thereby ensuring a finite energy flux at the crack tip. With this governing model, largely used by seismologists, it is possible to simulate a single earthquake event; but, in the absence of remote tectonic loading, it requires the introduction of an artificial procedure to initiate the rupture (i.e., to reach the failure stress point). In this article, by studying the dynamic rupture propagation and the solutions on the fault and on the free surface, I systematically compare three conceptually and algorithmically different nucleation strategies widely adopted in the literature: the imposition of an initially constant rupture speed, the introduction of a shear stress asperity, and the perturbation to the initial particle velocity field. My results show that, contrary to supershear ruptures, which tend to forget their origins, subshear ruptures are quite sensitive to the adopted nucleation procedure, which can bias the runaway rupture. I confirm that the most gradual transition from imposed nucleation and spontaneous propagation is obtained by initially forcing the rupture to expand at a properly chosen, constant speed (0.75 times the Rayleigh speed). I also numerically demonstrate that a valid alternative to this strategy is an appropriately smoothed, elliptical shear stress asperity. Moreover, I evaluate the optimal size of the nucleation patch where the procedure is applied; the simulations indicate that its size has to equal the critical distance of Day (1982) in the case of supershear ruptures and to exceed it in the case of subshear ruptures.

The transition of dynamic rupture styles in elastic media under velocity‐weakening friction

“…Hence, dynamic rupture simulations often adopt artificial nucleation procedures. Several studies have revealed that details of rupture initiation can significantly affect the subsequent rupture propagation [Festa and Vilotte, 2006;Shi et al, 2008;Lu et al, 2009;Bizzarri, 2010]. Here we assess the role of nucleation conditions in determining the rupture style by studying systematically the influence of the size of nucleation.…”

Landers 1992 "reloaded": Integrative dynamic earthquake rupture modeling

“…Within this zone, the friction coefficient is gradually reduced from its static to its dynamic value over a specified time of 0.5 s [Bizzarri , 2010]. Outside this zone, forced rupture is smoothly overtaken by spontaneous rupture.…”

“…Interestingly, our numerical tests reveal that fully elastic simulations can partially emulate the reference simulation when increasing the critical slip distance D c . In this case, moment-release rate and rupture transfer dynamics are preserved, but exhibit slower rupture speeds and longer delays when transferring to adjacent segments due to an increased critical size [e.g., Ampuero et al, 2002;Bizzarri , 2010;Galis et al, 2014] to initiate self-sustained rupture by dynamic triggering.…”

The mechanics of lubricated faults: Insights from 3‐D numerical models

“…The governing model assumed here (equation (10)) does not allow for a description of the hardening stage of the fault, unlike other friction models [e.g., Ruina, 1983;Sone and Shimamoto, 2009]. As thoroughly discussed in Bizzarri [2010a], to simulate the ongoing rupture propagation it is necessary to initially impose the nucleation, which is obtained by initially forcing the rupture to enlarge at a constant rupture speed (which equals to 2.4 km/s in our simulations). Also the breakdown stress drop is imposed in the nucleation patch; Dt b = t u À t f = 57.24 MPa.…”