Dynamic Rupture Models for the Southern San Andreas Fault

Abstract Dynamic rupture, and resultant ground motions up to 0.25 Hz, are simulated for an M w  7.6 earthquake on the southern San Andreas fault. Spontaneous rupture is modeled with slip-weakening friction, and 3D viscoelastic wave solutions are computed with a support-operator numerical method. The initial traction model is derived from inversions of the M w  7.3 1992 Landers strong ground-motion records, and borrows heavily from that used for the TeraShake2 simulations by Olsen et al. (2008). Heterogeneity in the traction model leads to focusing of the rupture front, and the focusing produces cases of supershear rupture velocity in asperities (areas of high initial traction), as well as cases of high peak slip rate and cohesive zone contraction in antiasperities. Separate solutions are computed for version 3.0 and 4.0, respectively, of the Southern California Earthquake Center Community Velocity Model (SCEC-CVM). We also compare the case of a flat ground surface (a common simplification made for finite-difference simulations) to the case of the ground surface conformed to regional topography. The overall distribution of simulated ground motion intensity is consistent with that derived from the empirical model of Campbell and Bozorgnia (2008), in the sense that the bulk of simulated pseudospectral velocity (PSV) values are within the 68% confidence intervals of the empirical model. Simulated PSVs corresponding to low probability in the empirical model are principally associated with basin wave-guide and directivity effects. An important example, first identified by the TeraShake1 simulations (Olsen et al. , 2006), is the stronger than expected ground motions at the site of Montebello due to a basin wave-guide effect. We find that this effect is lessened for version 4.0 of the SCEC-CVM, relative to version 3.0, due to a shallower model for the Chino basin.

[1]  Terry E. Tullis,et al.  Self-healing slip pulses in dynamic rupture models due to velocity-dependent strength , 1996, Bulletin of the Seismological Society of America.

[2]  N. Beeler,et al.  Constitutive relationships and physical basis of fault strength due to flash heating , 2008 .

[3]  Morgan T. Page,et al.  Effects of Large-Scale Surface Topography on Ground Motions, as Demonstrated by a Study of the San Gabriel Mountains, Los Angeles, California , 2007 .

[4]  Philip J. Maechling,et al.  TeraShake2: Spontaneous Rupture Simulations of Mw 7.7 Earthquakes on the Southern San Andreas Fault , 2008 .

[5]  L. Seeber,et al.  The San Andreas Fault system through the Transverse Ranges as illuminated by earthquakes , 1995 .

[6]  Thomas C. Hanks,et al.  A Bilinear Source-Scaling Model for M-log A Observations of Continental Earthquakes , 2002 .

[7]  David D. Jackson,et al.  Seismic hazards in southern California: probable earthquakes, 1994 to 2024 , 1996 .

[8]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[9]  Jean-Bernard Minster,et al.  A support-operator method for 3-D rupture dynamics , 2009 .

[10]  Christopher R. Bradley,et al.  Memory-Efficient Simulation of Anelastic Wave Propagation , 2001 .

[11]  Kim B. Olsen,et al.  On the implementation of perfectly matched layers in a three‐dimensional fourth‐order velocity‐stress finite difference scheme , 2003 .

[12]  J. Carlson,et al.  Rupture Pulse Characterization: Self-Healing, Self-Similar, Expanding Solutions in a Continuum Model of Fault Dynamics , 2000 .

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

[14]  J. Bérenger.  Three-Dimensional Perfectly Matched Layer for the Absorption of Electromagnetic Waves , 1996 .

[15]  P. Maechling,et al.  Strong shaking in Los Angeles expected from southern San Andreas earthquake , 2006 .

[16]  J. D. Eshelby.  Uniformly Moving Dislocations , 1949 .

[17]  D. Dreger,et al.  Inverse Kinematic and Forward Dynamic Models of the 2002 Denali Fault Earthquake, Alaska , 2004 .

[18]  Robert W. Graves,et al.  The SCEC Southern California Reference Three-Dimensional Seismic Velocity Model Version 2 , 2000 .

[19]  Arthur Frankel,et al.  Three-dimensional simulations of ground motions in the San Bernardino Valley, California, for hypothetical earthquakes on the San Andreas fault , 1993, Bulletin of the Seismological Society of America.

[20]  J. Rice,et al.  The growth of slip surfaces in the progressive failure of over-consolidated clay , 1973, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[21]  S. Harmsen,et al.  Documentation for the 2002 update of the national seismic hazard maps , 2002 .

[22]  K. Campbell,et al.  NGA Ground Motion Model for the Geometric Mean Horizontal Component of PGA, PGV, PGD and 5% Damped Linear Elastic Response Spectra for Periods Ranging from 0.01 to 10 s , 2008 .

[23]  Robert W. Graves,et al.  Three-dimensional finite-difference modeling of the San Andreas fault: Source parameterization and ground-motion levels , 1998, Bulletin of the Seismological Society of America.

[24]  S. Day,et al.  Comparison of finite difference and boundary integral solutions to three‐dimensional spontaneous rupture , 2005 .

[25]  D. J. Andrews,et al.  Rupture models with dynamically determined breakdown displacement , 2004 .

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

[27]  R. Madariaga,et al.  Dynamic modeling of the 1992 Landers earthquake , 2001 .

[28]  D. Wells,et al.  New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement , 1994, Bulletin of the Seismological Society of America.

[29]  Toshihiko Shimamoto,et al.  High‐velocity frictional properties of gabbro , 1997 .

[30]  P. Somerville.  REVIEW OF MAGNITUDE - RUPTURE AREA SCALING OF CRUSTAL EARTHQUAKES , 2007 .

[31]  David J. Wald,et al.  Dynamic stress changes during earthquake rupture , 1998, Bulletin of the Seismological Society of America.

[32]  Geoffrey Palarz Ely,et al.  A support-operator method for viscoelastic wave modelling in 3-D heterogeneous media , 2008 .

[33]  Steven M. Day,et al.  Efficient simulation of constant Q using coarse-grained memory variables , 1998, Bulletin of the Seismological Society of America.

[34]  Monica D. Kohler,et al.  Mantle Heterogeneities and the SCEC Reference Three-Dimensional Seismic Velocity Model Version 3 , 2003 .

[35]  Glenn P. Biasi,et al.  Wrightwood and the earthquake cycle: What a long recurrence record tells us about how faults work , 2004 .

[36]  E. Dunham,et al.  A Supershear Transition Mechanism for Cracks , 2003, Science.

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

[38]  Thomas C. Hanks,et al.  M-logA Observations for Recent Large Earthquakes , 2008 .

[39]  J. Rice.  Heating and weakening of faults during earthquake slip , 2006 .

[40]  Gregory C. Beroza,et al.  Short slip duration in dynamic rupture in the presence of heterogeneous fault properties , 1996 .

[41]  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.