Effect of the time-weakening friction law during the nucleation process

1.   Introduction

The earthquake nucleation process, which is difficult to observe directly in natural earthquakes, is an essential question in dynamic rupture problems. Numerical simulations are usually preferred to investigate the problem. During the spontaneous rupture process, slip starts when the shear traction reaches the friction strength at each point on a given fault following a given constitution law. Both slip-weakening friction law and time-weakening friction law are widely used in rupture dynamics (Andrews 1976; Bizzarri 2010). However, the effect of time-weakening friction law in controlling rupture behaviors during the nucleation process is seldom studied.

Starting from the pioneering work of Ida (1972), and Andrews (1976), the slip-weakening friction law is one of the most adopted constitutive laws in spontaneous rupture simulations (Harris et al. 2009). Although conceptually simple, the corporation of the slip-weakening law and numerical codes can predict the main characteristics of dynamic rupture scenarios on fault systems (Harris et al. 1991; Kase and Kuge 2001; Hu et al. 2016a). The slip-weakening behavior has also been noticed using the rate- and state-dependent friction law (Cocco and Bizzarri 2002). Compared to the rate- and state-dependent friction law (Dieterich 1979; Ruina 1983), it is easier to accomplish in numerical codes, although a high shear stress patch has to be located on the fault to trigger rupture in the slip-weakening friction law.

Thus, the slip-weakening friction law cannot be used to simulate repeated ruptures on a planar fault unless there is external time variable loading inserted (Bizzarri 2010). It is usually assumed that an asperity exists on the fault plane to trigger the rupture, in which the initial shear stress is larger than the friction strength. Under this assumption rupture starts simultaneously from all the points in the asperity. There exists another nucleation method referred as the time-weakening friction law (Bizzarri 2010), in which rupture starts from one point in the asperity, and propagates with a given rupture velocity toward the boundary of the nucleation patch. The starting point is treated as the hypocenter. Although Bizzarri (2010) compared the rupture speed dependence on the two nucleation methods, it is still unclear about the factors influencing the time-weakening friction law and equivalence of these two constitution laws in controlling rupture behaviors.

In this paper, we compared dynamic rupture behaviors on a single fault plane in a homogeneous half-space following the time-weakening and slip-weakening friction laws, respectively. The characteristic time and rupture speed in the asperity are two key factors in controlling the rupture fronts following the time-weakening friction law. To have a similar rupture front contour on a planar fault, rupture following the slip-weakening law should have a smaller size of nucleation asperity compared to the result of the time-weakening law. The distributions of peak ground motions are also studied following the two friction laws.

2.   Two different nucleation methods

The slip-weakening (SW) friction law was first proposed by Ida (1972) and Andrews (1976). Its simplicity in use and accuracy in predicting dynamic rupture behaviors made it widely used in rupture dynamics (Andrews 1985; Madariaga et al. 1998; Zhang and Chen 2006a, b; Hu et al. 2014, 2016b; Xu et al. 2015). As shown in Eq. (1), the fault strength decreases linearly to the residual strength as slip increases from zero to the characteristic slip distance,

where u is the magnitude of the fault slip, τ_u is the fault strength, τ_f is the residual strength, and d₀ is the characteristic slip distance. However, such a constitution law does not contain any nucleation process. An asperity with a stress perturbation which makes the initial shear stress larger than the friction strength is needed to trigger the rupture. Thus, the fault traction is:

where τ₀ is the initial shear stress, Δτ_nuc is the additional shear stress inserted in the asperity to start the rupture. It should be noticed that the time-weakening (TW) friction law here is only used during the nucleation process. Rupture outside the nucleation patch also follows the SW law. The TW friction is specified as (Bizzarri 2010):

where τ^(SW) is given in Eq. (1) and τ^(TW) is given by,

where μ_s and μ_f are static and dynamic frictional coefficients, respectively. t₀ is the characteristic time. τ_n means the normal stress. And

where v_force is the rupture speed in the asperity, and t_force is the time at which rupture propagating at v_force reaches the point (y, z) from the hypocenter (y^H, z^H). Thus, the fault strength in the asperity decreases linearly as time increases according to the TW law, rather than the slip dependent in the SW law. Slip starts when fault strength decreases to the initial shear stress. Different from the nucleation strategy of the SW law, the initial shear stress using the TW law here is homogeneous through the whole fault plane.

3.   Models and methods

We construct a right-lateral strike-slip single fault in a homogeneous, isotropic, elastic half-space. The P wave speed is 6.0 km/s, and the S wave speed is 3.464 km/s. The density is 2670 kg/m³. As shown in Fig. 1, the fault is 30 km long, 15 km wide. The hypocenter is located at (y^H, z^H) = (7.5 km, -7.5 km). The stress and frictional parameters are shown in Table 1. As stated previously, rupture starts simultaneously from all the points in the nucleation asperity following the SW law, but it is triggered from the hypocenter and propagates at a given velocity toward the boundary of nucleation asperity following the TW law. For the SW law, rupture is simultaneously triggered by the higher initial shear stress τ₀+Δτ_nuc (Table 1) in the nucleation asperity. For the two constitution laws, we assume a circle nucleation asperity with a radius of R_nuc. The dynamic rupture process of the single fault is calculated using the finite difference method developed by Zhang et al. (2014). The spatial grid is set to be 100 m, and the time interval is 0.01 s.

Figure  1.  Illustration of a right-lateral strike-slip fault model in a half-space

DownLoad: Full-Size Img PowerPoint

Table  1.  Stress and frictional parameters

 | Show Table

DownLoad: CSV

4.   Rupture behavior controlled by the TW law

To understand the effects of the characteristic time t₀ in the TW law, three rupture fronts with different characteristic time t₀ of 0.3, 0.5 and 0.7 s are presented in Fig. 2. Rupture speed in the nucleation patch v_force is all set to be v_S here, where v_S is the S wave speed. A prominent feature is rupture with a smaller characteristic time propagates earlier. As shown in Eq. (4), a smaller t₀ means the fault strength in the asperity decreases faster to the residual strength. For all the three t₀ values, free-surface-induced supershear ruptures are all observed. Shortly after 7 s, rupture fronts encounter the end of the fault plane, 30 km along the strike direction. Figure 2 also shows the horizontal shear stresses at the hypocenter for the three cases. Different from the SW law, the initial shear stress at the hypocenter is 71.3 MPa, the same as that outside of the nucleation patch, which is far from the friction strength. As time increases, the friction strength decreases linearly. Slip starts when the friction strength decreases to the initial shear stress. And all the shear stress decreases to the residual stress at the characteristic time t₀. With a larger t₀, the shear stress decreases slower.

Figure  2.  (left) Rupture front contours of three different characteristic time t₀ following the TW law. (right) Shear stress versus time at the hypocenter for the three cases

DownLoad: Full-Size Img PowerPoint

Different from the effect of characteristic time, a larger v_force means a smaller t_force at the same point on a given fault (Eq. 5). Thus, at the same time, the fault strength in the asperity decreases faster with a larger v_force, which permits the rupture to propagate earlier on the fault plane. We compared the space-time plots of horizontal slip velocity along the fault trace (z = -7.5 km) and the rupture front contours following the TW and SW friction laws, respectively (Fig. 3). Three different values of v_force, including 0.75v_s, 1.73v_s (which equals to the P wave velocity) and ∞ are considered. Although the upper limit of rupture velocity is P wave velocity, the infinity rupture velocity here considered is only used to study the differences of the rupture behaviors caused by the friction laws themselves, because all the points in the nucleation area following the SW friction law start to rupture simultaneously. The characteristic time t₀ is set to be 0.5 s in the three cases following the TW law. However, three different values R_nuc of 1.3, 1.5 and 3.0 km are considered in the SW law. The plot also gives the S wave and P wave arrivals along the fault trace (white dashes).

Figure  3.  (left) Space-time plot horizontal slip velocity along the fault trace (z = -7.5 km). The dashed lines depict wave arrivals along the fault trace from the nucleation patch with three different values of v_force and R_nuc following the TW and SW friction laws, respectively. (right) Rupture front contours comparisons between two different friction laws

DownLoad: Full-Size Img PowerPoint

Two main features are noticed in Fig. 3. First, the rupture velocity in the nucleation patch is identical to the v_force. As shown in the space-time slip rate distribution in Fig. 3, when v_force is 0.75v_S, rupture propagates at sub-Rayleigh speed in the nucleation patch, which increases to P wave speed when v_force increases to 1.73v_S. If v_force is set to ∞, all the points in the asperity start to rupture simultaneously, which coincides with the nucleation process of the SW law.

Second, to have similar rupture fronts following the TW laws, rupture following the SW law should have a smaller nucleation radius R_nuc. As shown in Fig. 3, rupture following the TW law with a 3 km R_nuc and v_force of 0.75v_S shows similar characteristics of slip rate on the given fault trace with that following the SW law with R_nuc of 1.3 km. The free-surface-induced supershear has not reached the fault trace in the middle depth at 8 s, which is the same in the slip rate snaps shown in Fig. 3. When the radius of asperity R_nuc increases to 1.5 km in the SW law, the corresponding rupture following the TW law requires v_force of P wave velocity. Nearly at 8 s, the free-surface-induced supershear rupture can be observed in both slip rate snaps. If R_nuc increases to 3 km in the SW law, the rupture front following the TW law propagates slower than that of the SW law, even if the v_force is set to ∞. As noticed in Fig. 2, the shear stress of the TW law in the asperity decreases from 71.3 to 63.0 MPa, instead of 81.6 to 63.0 MPa of the SW law. Thus, the higher slip rate in the asperity is induced by the larger stress drop of the SW law.

An interesting feature is for the case with v_force of 1.73v_S, the rupture front travels slightly farther for the TW law at 1 s; however, the rupture front for the SW law travels goes beyond it after 1 s. Such differences may be caused by the transition process from TW to SW law at the boundary of the nucleation asperity (Eq. 3). Different from the TW law, stresses at the points just outside of the nucleation asperity have to rise from the initial shear stresses to the static friction strengths of the SW law, which takes some time. This helps to understand the difference between rupture fronts for the TW case with v_force of ∞ and the SW case with R_nuc of 3 km.

Another main factor that causes the smaller nucleation size of the SW law with respect to the TW law is the larger peak slip rate in the nucleation asperity. Figure 4 shows the slip rate versus time of the point in the center of the nucleation asperity of the different cases shown in Fig. 3. As the v_force increases in the TW law, the peak slip rate in the hypocenter also increases; however, even if the v_force is set to ∞, which is the same in the SW law, the peak slip rate is smaller than that of the SW law with a 3.0 km nucleation radius. For all the cases shown in Fig. 4, although their rupture fronts are similar, all peak slip rates of the SW law are larger than those of the TW law, which is caused by the larger stress drop in the SW law. A larger peak slip rate represents a larger fracture energy, which permits a rupture to propagate faster with the same nucleation size.

Figure  4.  Slip rate distribution of the single point located at the hypocenter with different cases shown in Fig. 3

DownLoad: Full-Size Img PowerPoint

5.   Peak ground velocity distribution

To further understand the differences between the TW and SW friction laws, we compared the peak ground velocity distributions on the free surface induced by the rupture following the TW friction law with a v_force of 0.75v_S and a R_nuc of 3 km. As shown in Fig. 5, the PGVx is the fault-perpendicular component, and the PGVy means the fault-parallel component. The single fault trace, 30 km long, starts from 0 km in the strike direction on the free surface. Strong rupture directivity can be noticed in the fault-perpendicular component. The point which has the maximum value of PGVx lies on the fault trace close to the end of the fault plane in the strike direction, in front of which a patch with strong PGVx can be observed. For the fault-parallel component PGVy, the largest value can be found in an area close to the nucleation patch.

Figure  5.  (left) Distribution of peak ground velocities on the free surface of rupture following the TW law with a v_force of 0.75v_S and a R_nuc of 3 km. (right) Perturbation of peak ground velocities (results of the SW law with a R_nuc of 3 km minus the result of Fig. 5 (left))

DownLoad: Full-Size Img PowerPoint

The distributions of peak ground velocity induced by the rupture with the same radius of asperity R_nuc of 3 km following the SW law are also calculated. Figure 5 also shows the perturbation of PGVx and PGVy between the two cases with the same R_nuc. For both the fault-perpendicular and fault-parallel components, the peak ground velocities are larger for the SW case. Compared to the results of the TW law with a v_force of 0.75v_S, the results of the SW case show stronger PGVx at the area outside of the nucleation area, which represent stronger rupture directivity. Moreover, the main difference of PGVy can be noticed close to the nucleation area near the fault trace.

6.   Conclusions

We have performed dynamic rupture simulations using the TW and SW friction laws, respectively, on a planar fault in a homogeneous half-space. Different from the SW law, two factors of the TW law, including the characteristic time t₀ and the rupture speed in the nucleation asperity v_force, control the main characteristics of rupture behaviors on a planar fault. Rupture with a smaller t₀ or a larger v_force propagates faster on a single fault plane following the TW law.

Rupture fronts following the TW law look similar with those of the SW law with a smaller nucleation asperity. There are two main reasons that cause the phenomenon. The first one is the different nucleation strategies. Rupture of the SW law starts simultaneously in the nucleation patch, while it is triggered in a single point (the hypocenter) in the TW friction law, and propagates at a given speed toward the boundary of the nucleation asperity. According to our simulations, the rupture following the TW law with a v_force of 0.75v_S and a R_nuc of 3 km has similar rupture front and peak ground velocity distribution with the rupture following the SW law with a R_nuc of 1.3 km. The other factor is that peak slip rates of the TW law are smaller than those of the SW law in the nucleation asperity. Even if the v_force is set to ∞, which is the same as that of the SW law, the rupture front also propagates slower than that of the SW law because of the lower slip rate value in the nucleation asperity following the TW friction law.

With the same radius of nucleation asperity, the peak ground velocity distribution near the fault trace caused by the rupture following the TW friction law shows similar characteristics with that of the SW friction law with the same nucleation radius, except the lower effect of rupture directivity and smaller value near the nucleation area in the fault-parallel component. Our work can provide useful insight to understanding the earthquake nucleation process.

Acknowledgements

This work is supported by the National Natural Science Foundation of China (Nos. 41504039, 41474037 and 41274053). All data used in this study are derived from numerical simulations.

Open Access

This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.