Ground-motion modeling of Hayward fault scenario earthquakes, part II: Simulation of long-period and broadband ground motions

Abstract We simulate long-period ( T >1.0–2.0 s) and broadband ( T >0.1 s) ground motions for 39 scenario earthquakes ( M w  6.7–7.2) involving the Hayward, Calaveras, and Rodgers Creek faults. For rupture on the Hayward fault, we consider the effects of creep on coseismic slip using two different approaches, both of which reduce the ground motions, compared with neglecting the influence of creep. Nevertheless, the scenario earthquakes generate strong shaking throughout the San Francisco Bay area, with about 50% of the urban area experiencing modified Mercalli intensity VII or greater for the magnitude 7.0 scenario events. Long-period simulations of the 2007 M w  4.18 Oakland earthquake and the 2007 M w  5.45 Alum Rock earthquake show that the U.S. Geological Survey’s Bay Area Velocity Model version 08.3.0 permits simulation of the amplitude and duration of shaking throughout the San Francisco Bay area for Hayward fault earthquakes, with the greatest accuracy in the Santa Clara Valley (San Jose area). The ground motions for the suite of scenarios exhibit a strong sensitivity to the rupture length (or magnitude), hypocenter (or rupture directivity), and slip distribution. The ground motions display a much weaker sensitivity to the rise time and rupture speed. Peak velocities, peak accelerations, and spectral accelerations from the synthetic broadband ground motions are, on average, slightly higher than the Next Generation Attenuation (NGA) ground-motion prediction equations. We attribute much of this difference to the seismic velocity structure in the San Francisco Bay area and how the NGA models account for basin amplification; the NGA relations may underpredict amplification in shallow sedimentary basins. The simulations also suggest that the Spudich and Chiou (2008) directivity corrections to the NGA relations could be improved by increasing the areal extent of rupture directivity with period.

[1]  D. Boore.  Stochastic simulation of high-frequency ground motions based on seismological models of the radiated spectra , 1983 .

[2]  N. Abrahamson,et al.  Summary of the Abrahamson & Silva NGA Ground-Motion Relations , 2008 .

[3]  Walter H. F. Smith,et al.  New, improved version of generic mapping tools released , 1998 .

[4]  Pengcheng Liu,et al.  Prediction of Broadband Ground-Motion Time Histories: Hybrid Low/High- Frequency Method with Correlated Random Source Parameters , 2006 .

[5]  N. Anders Petersson,et al.  Broadband Waveform Modeling of Moderate Earthquakes in the San Francisco Bay Area and Preliminary Assessment of the USGS 3D Seismic Velocity Model , 2008 .

[6]  Shuo Ma,et al.  Modeling of the Perfectly Matched Layer Absorbing Boundaries and Intrinsic Attenuation in Explicit Finite-Element Methods , 2006 .

[7]  Howard Bundock,et al.  Modified Mercalli Intensity Maps for the 1868 Hayward Earthquake Plotted in ShakeMap Format , 2008 .

[8]  Allen H. Olson,et al.  Finite faults and inverse theory with applications to the 1979 Imperial Valley earthquake , 1982 .

[9]  Stefan Nilsson,et al.  Stable Difference Approximations for the Elastic Wave Equation in Second Order Formulation , 2007, SIAM J. Numer. Anal..

[10]  E. Brabb,et al.  Compilation of 59 sonic and density logs from 51 oil test wells in the San Francisco Bay area, California , 1997 .

[11]  Douglas S. Dreger,et al.  Seismic remote sensing for the earthquake source process and near‐source strong shaking: A case study of the October 16, 1999 Hector Mine earthquake , 2000 .

[12]  N. Anders Petersson,et al.  Ground-Motion Modeling of the 1906 San Francisco Earthquake, Part II: Ground-Motion Estimates for the 1906 Earthquake and Scenario Events , 2007 .

[13]  J. Bardet,et al.  Kocaeli, Turkey, earthquake of August 17, 1999 reconnaissance report , 2000 .

[14]  A. Pitarka,et al.  Broadband Ground-Motion Simulation Using a Hybrid Approach , 2010 .

[15]  J. Bernard Minster,et al.  Numerical simulation of attenuated wavefields using a Padé approximant method , 1984 .

[16]  D. Appelö,et al.  A stable finite difference method for the elastic wave equation on complex geometries with free surfaces , 2007 .

[17]  N. Abrahamson,et al.  Characterizing Crustal Earthquake Slip Models for the Prediction of Strong Ground Motion , 1999 .

[18]  Arthur Frankel,et al.  Simulating strong motions of large earthquakes using recordings of small earthquakes: the Loma Prieta mainshock as a test case , 1995, Bulletin of the Seismological Society of America.

[19]  Badie Rowshandel.  Directivity Correction for the Next Generation Attenuation (NGA) Relations , 2010 .

[20]  B. Chiou,et al.  Directivity in NGA Earthquake Ground Motions: Analysis Using Isochrone Theory , 2008 .

[21]  P. Segall,et al.  Slip in the 1868 Hayward earthquake from the analysis of historical triangulation data , 1996 .

[22]  BrianS-J. Chiou,et al.  An NGA Model for the Average Horizontal Component of Peak Ground Motion and Response Spectra , 2008 .

[23]  Robert W. Graves,et al.  Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences , 1996, Bulletin of the Seismological Society of America.

[24]  William Hung Kan Lee,et al.  International handbook of earthquake and engineering seismology , 2002 .

[25]  Gregory C. Beroza,et al.  Testing Community Velocity Models for Southern California Using the Ambient Seismic Field , 2008 .

[26]  Timothy E. Dawson,et al.  Uniform California Earthquake Rupture Forecast, Version 2 (UCERF 2) , 2009 .

[27]  Jonathan P. Stewart,et al.  Broadband simulations for Mw 7.8 southern San Andreas earthquakes: Ground motion sensitivity to rupture speed , 2008 .

[28]  W. Ellsworth,et al.  Near-Field Ground Motion of the 2002 Denali Fault, Alaska, Earthquake Recorded at Pump Station 10 , 2004 .

[29]  T. Brocher.  Compressional and Shear-Wave Velocity versus Depth Relations for Common Rock Types in Northern California , 2008 .

[30]  T. Iwata,et al.  Rupture Process of the 1999 Kocaeli, Turkey, Earthquake Estimated from Strong-Motion Waveforms , 2002 .

[31]  Ares J. Rosakis,et al.  How fast is rupture during an earthquake? New insights from the 1999 Turkey Earthquakes , 2001 .

[32]  G. Atkinson,et al.  Ground-Motion Prediction Equations for the Average Horizontal Component of PGA, PGV, and 5%-Damped PSA at Spectral Periods between 0.01 s and 10.0 s , 2008 .

[33]  Barbara Romanowicz,et al.  Toward real-time estimation of regional moment tensors , 1996, Bulletin of the Seismological Society of America.

[34]  W. Mooney,et al.  Data report for seismic refraction surveys conducted from 1980 to 1982 in the Livermore Valley and the Santa Cruz Mountains, California , 1999 .

[35]  T. Brocher.  A Regional View of Urban Sedimentary Basins in Northern California Based on Oil Industry Compressional-Wave Velocity and Density Logs , 2005 .

[36]  A. Pitarka,et al.  BROADBAND TIME HISTORY SIMULATION USING A HYBRID APPROACH , 2002 .

[37]  William A. Bryant,et al.  A Site-Conditions Map for California Based on Geology and Shear-Wave Velocity , 2000 .

[38]  David M. Boore,et al.  Site amplifications for generic rock sites , 1997, Bulletin of the Seismological Society of America.

[39]  N. Abrahamson,et al.  Nonlinear Site Amplification Factors for Constraining the NGA Models , 2008 .

[40]  Michel Bouchon,et al.  Observation of Long Supershear Rupture During the Magnitude 8.1 Kunlunshan Earthquake , 2003, Science.

[41]  D. Wald,et al.  Best practices for using macroseismic intensity and ground motion to intensity conversion equations for hazard and loss models , 2010 .

[42]  Ralph J. Archuleta,et al.  A faulting model for the 1979 Imperial Valley earthquake , 1984 .

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

[44]  C. Thurber,et al.  Three-dimensional P wave velocity model for the San Francisco Bay region, California , 2007 .

[45]  John G. Anderson,et al.  A MODEL FOR THE SHAPE OF THE FOURIER AMPLITUDE SPECTRUM OF ACCELERATION AT HIGH FREQUENCIES , 1984 .

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

[47]  David Carver,et al.  Modeling and Validation of a 3D Velocity Structure for the Santa Clara Valley, California, for Seismic-Wave Simulations , 2006 .

[48]  Ralph J. Archuleta,et al.  Evidence for a Supershear Transient during the 2002 Denali Fault Earthquake , 2004 .

[49]  R. Graymer.  Geologic Map and Map Database of the Oakland Metropolitan Area, Alameda, Contra Costa, and San Francisco Counties, California , 2000 .

[50]  Mustafa Aktar,et al.  Seismic imaging of the 1999 Izmit (Turkey) Rupture inferred from the near‐fault recordings , 2000 .

[51]  M. Bouchon,et al.  The Aftershock Signature of Supershear Earthquakes , 2008, Science.

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

[53]  C. W. Stover,et al.  Seismicity of the United States, 1568-1989 (revised) , 1993 .

[54]  Thomas H. Heaton,et al.  Characterization of Near-Source Ground Motions with Earthquake Simulations , 2001 .

[55]  S. Hartzell,et al.  Simulated ground motion in Santa Clara Valley, California, and vicinity from M≥6.7 scenario earthquakes , 2008 .

[56]  Thomas M. Brocher,et al.  Ground-Motion Modeling of the 1906 San Francisco Earthquake, Part I: Validation Using the 1989 Loma Prieta Earthquake , 2008 .

[57]  W. Bakun.  Seismic activity of the San Francisco Bay region , 1999, Bulletin of the Seismological Society of America.

[58]  Robert W. Graves,et al.  Ground-motion modeling of Hayward fault scenario earthquakes, part I: Construction of the suite of scenarios , 2010 .

[59]  Thomas H. Heaton,et al.  Near-Source Ground Motions from Simulations of Sustained Intersonic and Supersonic Fault Ruptures , 2004 .

[60]  P. Spudich,et al.  Direct observation of rupture propagation during the 1979 Imperial Valley earthquake using a short baseline accelerometer array , 1984 .

[61]  Shawn Larsen,et al.  Calculation of broadband time histories of ground motion: Comparison of methods and validation using strong-ground motion from the 1994 Northridge earthquake , 1999, Bulletin of the Seismological Society of America.

[62]  Kim B. Olsen,et al.  Model for Basin Effects on Long-Period Response Spectra in Southern California , 2008 .

[63]  Shawn Larsen,et al.  Moderate Earthquake Ground-Motion Validation in the San Francisco Bay Area , 2010 .

[64]  N. Abrahamson,et al.  Modification of Empirical Strong Ground Motion Attenuation Relations to Include the Amplitude and Duration Effects of Rupture Directivity , 1997 .

[65]  Joakim O. Blanch,et al.  Viscoelastic finite-difference modeling , 1994 .