| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|
|
|
|||||||
|
|||||||||
Laboratoire de Géologie URA CNRS 1716. Ecole Normale Supérieure, 24 rue Lhomond, 75231 Paris Cedex 05, France
Laboratoire de Sismologie, Université Denis Diderot de Paris 7
Institute for Crustal Studies UC Santa Barbara, Santa Barbara, California 93016-1100
Abstract
We propose a fourth-order staggered-grid finite-difference method to study dynamic faulting in three dimensions. The method uses an implementation of the boundary conditions on the fault that allows the use of general friction models including slip weakening and rate dependence. Because the staggered-grid method defines stresses and particle velocities at different grid points, we preserve symmetry by implementing a two-grid-row "thick" fault zone. Slip is computed between points located at the borders of the fault zone, while the two components of shear traction on the fault are forced to be symmetric inside the fault zone. We study the properties of the numerical method comparing our simulations with well-known properties of seismic ruptures in 3D. Among the properties that are well modeled by our method are full elastic-wave interactions, frictional instability, rupture initiation from a finite initial patch, spontaneous rupture growth at subsonic and supersonic speeds, as well as healing by either stopping phases or rate-dependent friction. We use this method for simulating spontaneous rupture propagation along an arbitrarily loaded planar fault starting from a localized asperity on circular and rectangular faults. The shape of the rupture front is close to elliptical and is systematically elongated in the inplane direction of traction drop. This elongation is due to the presence of a strong shear stress peak that moves ahead of the rupture in the in-plane direction. At high initial stresses the rupture front becomes unstable and jumps to super-shear speeds in the direction of in-plane shear. Another interesting effect is the development of relatively narrow rupture fronts due to the presence of rate-weakening friction. The solutions for the "thick fault" boundary conditions scale with the slip-weakening distance (D0) and are stable and reproducible for D0 greater than about 4 in terms of 2Tu/µ x 
x. Finally, a comparison of scalar and vector boundary conditions for the friction shows that slip is dominant along the direction of the prestress, with the largest deviations in slip-rate direction occurring near the rupture front and the edges of the fault.
This article has been cited by other articles:
|
|
 |
|
  S. M. Day, S. H. Gonzalez, R. Anooshehpoor, and J. N. Brune Scale-Model and Numerical Simulations of Near-Fault Seismic Directivity Bulletin of the Seismological Society of America, June 1, 2008; 98(3): 1186 - 1206. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
S. M. Corish, C. R. Bradley, and K. B. Olsen Assessment of a Nonlinear Dynamic Rupture Inversion Technique Applied to a Synthetic Earthquake Bulletin of the Seismological Society of America, June 1, 2007; 97(3): 901 - 914. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
L. F. F. Miguel and J. D. Riera A Constitutive Criterion for the Fault: Modified Velocity-Weakening Law Bulletin of the Seismological Society of America, June 1, 2007; 97(3): 915 - 925. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
L. A. Dalguer and S. M. Day Comparison of Fault Representation Methods in Finite Difference Simulations of Dynamic Rupture Bulletin of the Seismological Society of America, October 1, 2006; 96(5): 1764 - 1778. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
D. E. Smith, B. T. Aagaard, and T. H. Heaton Teleseismic Body Waves from Dynamically Rupturing Shallow Thrust Faults: Are They Opaque for Surface-Reflected Phases? Bulletin of the Seismological Society of America, June 1, 2005; 95(3): 800 - 817. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
L. R. Johnson and R. M. Nadeau Asperity Model of an Earthquake: Dynamic Problem Bulletin of the Seismological Society of America, February 1, 2005; 95(1): 75 - 108. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
B. T. Aagaard and T. H. Heaton Near-Source Ground Motions from Simulations of Sustained Intersonic and Supersonic Fault Ruptures Bulletin of the Seismological Society of America, December 1, 2004; 94(6): 2064 - 2078. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
The Importance of the Dynamic Source Effects on Strong Ground Motion during the 1999 Chi-Chi, Taiwan, Earthquake: Brief Interpretation of the Damage Distribution on Buildings Bulletin of the Seismological Society of America, October 1, 2004; 91(5): 1112 - 1127.
|
|
|
|
|
|
B. T. Aagaard, J. F. Hall, and T. H. Heaton Effects of Fault Dip and Slip Rake Angles on Near-Source Ground Motions: Why Rupture Directivity Was Minimal in the 1999 Chi-Chi, Taiwan, Earthquake Bulletin of the Seismological Society of America, February 1, 2004; 94(1): 155 - 170. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
H. Aochi and R. Madariaga The 1999 Izmit, Turkey, Earthquake: Nonplanar Fault Structure, Dynamic Rupture Process, and Strong Ground Motion Bulletin of the Seismological Society of America, June 1, 2003; 93(3): 1249 - 1266. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
Stress-Breakdown Time and Slip-Weakening Distance Inferred from Slip-Velocity Functions on Earthquake Faults Bulletin of the Seismological Society of America, February 1, 2003; 93(1): 264 - 282.
|
|
|
|
|
|
Effect of a Shallow Weak Zone on Fault Rupture: Numerical Simulation of Scale-Model Experiments Bulletin of the Seismological Society of America, December 1, 2002; 92(8): 3022 - 3041.
|
|
|
|
|
|
Dynamic Earthquake Ruptures in the Presence of Lithostatic Normal Stresses: Implications for Friction Models and Heat Production Bulletin of the Seismological Society of America, December 1, 2001; 91(6): 1765 - 1796.
|
|
|
|
|
|
D. J. Andrews Test of two methods for faulting in finite-difference calculations Bulletin of the Seismological Society of America, August 1, 1999; 89(4): 931 - 937. [Abstract] [PDF]
|
|
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
