Solving the dynamic rupture problem with different numerical approaches and constitutive laws

Summary We study the dynamic initiation, propagation and arrest of a 2-D in-plane shear rupture by solving the elastodynamic equation by using both a boundary integral equation method and a finite difference approach. For both methods we adopt different constitutive laws: a slip-weakening (SW) law, with constant weakening rate, and rate- and state-dependent friction laws (Dieterich–Ruina). Our numerical procedures allow the use of heterogeneous distributions of constitutive parameters along the fault for both formulations. We first compare the two solution methods with an SW law, emphasizing the required stability conditions to achieve a good resolution of the cohesive zone and to avoid artificial complexity in the solutions. Our modelling results show that the two methods provide very similar time histories of dynamic source parameters. We point out that, if a careful control of resolution and stability is performed, the two methods yield identical solutions. We have also compared the rupture evolution resulting from an SW and a rate- and state-dependent friction law. This comparison shows that despite the different constitutive formulations, a similar behaviour is simulated during the rupture propagation and arrest. We also observe a crack tip bifurcation and a jump in rupture velocity (approaching the P-wave speed) with the Dieterich–Ruina (DR) law. The rupture arrest at a barrier (high strength zone) and the barrier-healing mechanism are also reproduced by this law. However, this constitutive formulation allows the simulation of a more general and complex variety of rupture behaviours. By assuming different heterogeneous distributions of the initial constitutive parameters, we are able to model a barrier-healing as well as a self-healing process. This result suggests that if the heterogeneity of the constitutive parameters is taken into account, the different healing mechanisms can be simulated. We also study the nucleation phase duration Tn, defined as the time necessary for the crack to reach the half-length lc. We compare the Tn values resulting from distinct simulations calculated using different constitutive laws and different sets of constitutive parameters. Our results confirm that the DR law provides a different description of the nucleation process than the SW law adopted in this study. We emphasize that the DR law yields a complete description of the rupture process, which includes the most prominent features of SW.

[1]  Keiiti Aki,et al.  A numerical study of two-dimensional spontaneous rupture propagation , 1977 .

[2]  D. J. Andrews,et al.  Bulletin of the Seismological Society of America , 1985 .

[3]  Raul Madariaga,et al.  Dynamic faulting under rate-dependent friction , 1994 .

[4]  A. Ruina,et al.  Slip motion and stability of a single degree of freedom elastic system with rate and state dependent friction , 1984 .

[5]  Gutuan Zheng,et al.  Self-healing slip pulse on a frictional surface , 1995 .

[6]  J. Rice,et al.  Slip complexity in earthquake fault models. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[7]  Kim B. Olsen,et al.  Three-Dimensional Dynamic Simulation of the 1992 Landers Earthquake , 1997 .

[8]  P. Bodin,et al.  Earthquake triggering by transient and static deformations , 1998 .

[9]  T. Heaton.  Evidence for and implications of self-healing pulses of slip in earthquake rupture , 1990 .

[10]  B. Shibazaki,et al.  Slip-dependent friction law and nucleation processes in earthquake rupture , 1992 .

[11]  R. Simpson,et al.  Stress transfer by the 1988–1989 M = 5.3 and 5.4 Lake Elsman foreshocks to the Loma Prieta fault: Unclamping at the site of peak mainshock slip , 1999 .

[12]  J. Rice.  Spatio‐temporal complexity of slip on a fault , 1993 .

[13]  P. Okubo.  Dynamic rupture modeling with laboratory‐derived constitutive relations , 1989 .

[14]  R. Madariaga,et al.  Rupture dynamics of a planar fault in a 3D elastic medium: Rate- and slip-weakening friction , 1998, Bulletin of the Seismological Society of America.

[15]  Steven M. Day,et al.  Effects of a low-velocity zone on a dynamic rupture , 1997, Bulletin of the Seismological Society of America.

[16]  J. Dieterich,et al.  Implications of fault constitutive properties for earthquake prediction. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[17]  A. Papageorgiou,et al.  A specific barrier model for the quantitative description of inhomogeneous faulting and the prediction of strong ground motion. I. Description of the model , 1983 .

[18]  C. Marone,et al.  Earthquake nucleation on model faults with rate‐ and state‐dependent friction: Effects of inertia , 1996 .

[19]  Kim B. Olsen,et al.  Influence of friction and fault geometry on earthquake rupture , 2000 .

[20]  Yoshiaki Ida,et al.  Cohesive force across the tip of a longitudinal‐shear crack and Griffith's specific surface energy , 1972 .

[21]  P. Mirwald,et al.  The Curie Point and the α-γ transition of iron to 53 Kbar — A reexamination , 1979 .

[22]  D. J. Andrews,et al.  Rupture velocity of plane strain shear cracks , 1976 .

[23]  J. Dieterich.  Earthquake nucleation on faults with rate-and state-dependent strength , 1992 .

[24]  M. E. Hanson,et al.  Difference equations for two-dimensional elastic flow , 1968 .

[25]  J. Rice,et al.  Conditions under which velocity-weakening friction allows a self-healing versus a cracklike mode of rupture , 1998, Bulletin of the Seismological Society of America.

[26]  Teruo Yamashita,et al.  A cohesive zone model for dynamic shear faulting based on experimentally inferred constitutive relation and strong motion source parameters , 1989 .

[27]  Dynamic fault rupture and stress recovery processes in continental crust under depth-dependent shear strength and frictional parameters , 1992 .

[28]  B. Shibazaki,et al.  Transition process from nucleation to high-speed rupture propagation: scaling from stick-slip experiments to natural earthquakes , 2002 .

[29]  P. Spudich,et al.  What Can Strong-Motion Data Tell Us about Slip-Weakening Fault-Friction Laws? , 2000 .

[30]  Nonuniformity of the constitutive law parameters for shear rupture and quasistatic nucleation to dynamic rupture: a physical model of earthquake generation processes. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[31]  Steven M. Day,et al.  Three-dimensional finite difference simulation of fault dynamics: Rectangular faults with fixed rupture velocity , 1982 .

[32]  Terry E. Tullis,et al.  The roles of time and displacement in the evolution effect in rock friction , 1994 .

[33]  Breaking of a single asperity: Rupture process and seismic radiation , 1983 .

[34]  Y. Umeda.  The bright spot of an earthquake , 1992 .

[35]  J. Boatwright,et al.  Frictional constraints on crustal faulting , 1996 .

[36]  G. I. Barenblatt.  Concerning equilibrium cracks forming during brittle fracture. The stability of isolated cracks. Relationships with energetic theories , 1959 .

[37]  A. Ruina,et al.  Stability of Steady Frictional Slipping , 1983 .

[38]  J. Dieterich.  A constitutive law for rate of earthquake production and its application to earthquake clustering , 1994 .

[39]  Carlson,et al.  Mechanical model of an earthquake fault. , 1989, Physical review. A, General physics.

[40]  M. Ohnaka,et al.  Scaling of the shear rupture process from nucleation to dynamic propagation: Implications of geometric irregularity of the rupturing surfaces , 1999 .

[41]  M. Campillo,et al.  Influence of the shape of the friction law and fault finiteness on the duration of initiation , 1999 .

[42]  R. Madariaga,et al.  Modeling dynamic rupture in a 3D earthquake fault model , 1998, Bulletin of the Seismological Society of America.

[43]  Kojiro Irikura,et al.  Some characteristics of the stress field of the 1995 Hyogo‐ken Nanbu (Kobe) earthquake , 1998 .

[44]  John R. Rice,et al.  Dynamic motion of a single degree of freedom system following a rate and state dependent friction law , 1986 .

[45]  G. Spada,et al.  Spherical versus flat models of coseismic and postseismic deformations , 1999 .

[46]  D. J. Andrews,et al.  Properties and implications of dynamic rupture along a material interface , 1998, Bulletin of the Seismological Society of America.

[47]  John R. Rice,et al.  Crustal Earthquake Instability in Relation to the Depth Variation of Frictional Slip Properties , 1986 .

[48]  J. Rice,et al.  Stability of quasi-static slip in a single degree of freedom elastic system with rate and state dependent friction , 1999 .

[49]  Chris Marone,et al.  The depth of seismic faulting and the upper transition from stable to unstable slip regimes , 1988 .

[50]  M. F. Linker,et al.  Effects of variable normal stress on rock friction: Observations and constitutive equations , 1992 .

[51]  Y. Iio.  Observations of the slow initial phase generated by microearthquakes: Implications for earthquake nucleation and propagation , 1995 .

[52]  D. Andrews.  Dynamic growth of mixed-mode shear cracks , 1994, Bulletin of the Seismological Society of America.

[53]  Keiiti Aki,et al.  Fault plane with barriers: A versatile earthquake model , 1977 .

[54]  D. J. Andrews,et al.  A numerical study of tectonic stress release by underground explosions , 1973 .

[55]  Gregory C. Beroza,et al.  Properties of the seismic nucleation phase , 1996 .

[56]  C. Marone.  LABORATORY-DERIVED FRICTION LAWS AND THEIR APPLICATION TO SEISMIC FAULTING , 1998 .

[57]  A. Ruina.  Slip instability and state variable friction laws , 1983 .

[58]  W. Ellsworth,et al.  Seismic Evidence for an Earthquake Nucleation Phase , 1995, Science.

[59]  M. Takeo,et al.  Determination of constitutive relations of fault slip based on seismic wave analysis , 1997 .

[60]  Michel Campillo,et al.  Initiation of in-plane shear instability under slip-dependent friction , 1997, Bulletin of the Seismological Society of America.

[61]  M. Ohnaka.  Critical Size of the Nucleation Zone of Earthquake Rupture Inferred from Immediate Foreshock Activity , 1993 .

[62]  Yehuda Ben-Zion,et al.  Wrinkle-like slip pulse on a fault between different , 1997 .

[63]  D. J. Andrews,et al.  Rupture propagation with finite stress in antiplane strain , 1976 .