Magnitude-Scaling Rate in Ground-Motion Prediction Equations for Response Spectra from Large, Shallow Crustal Earthquakes

Abstract We have evaluated the magnitude-scaling rates (MSRs, the rate of increase in the predicted spectrum with increasing moment magnitude) of five modern ground-motion prediction equations (GMPEs) for response spectra, including four Next Generation Attenuation (NGA) models and a model developed for data from Japan. We have found that MSRs for crustal earthquakes with a moment magnitude over 7 vary significantly among the five models, by a factor of 2–3 for some cases. The variation of MSRs among the four NGA models is alarmingly large, considering that they were derived from largely the same dataset and used the same site parameters. We have selected 641 strong-motion records from shallow crustal earthquakes with a moment magnitude over 7 from the NGA dataset and 69 from the 2008 Wenchuan, China, earthquake with a moment magnitude of 7.9, all within a short distance of 200 km from the fault rupture plane. To illustrate the extreme extent of magnitude scaling, we have fitted an attenuation model without a magnitude term to this dataset, based on an observation that an increase in one moment magnitude unit is related to an increase in fault length by a factor of 5–10, a significant increase in shaking duration, which cannot be fully accounted for by response spectra, and limited or no increase in ground-motion amplitude. The statistical analyses of the results indeed suggest that a zero magnitude scaling at spectral periods over 0.6 s may be reasonable for our dataset, while the required apparent magnitude scaling at short periods may be due to other factors such as stress drop, a parameter that is not used in any of the models considered. We will provide some plausible explanations for the possible zero magnitude scaling that can be considered as the lower limit of the uncertain magnitude-scaling rate.

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

[2]  P. Somerville,et al.  The influence of critical Moho Reflections on strong ground motions recorded in San Francisco and Oakland during the 1989 Loma Prieta Earthquake , 1990 .

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

[4]  Y. Tsai,et al.  Large Effects of Moho Reflections (SmS) on Peak Ground Motion in Northwestern Taiwan , 2009 .

[5]  B. Bradley,et al.  THE MW 7.6 DUSKY SOUND EARTHQUAKE OF 2009: PRELIMINARY REPORT , 2010 .

[6]  Xiao Wen An,et al.  A Preliminary Study on the Near-Source Strong-Motion Characteristics of the Great 2008 Wenchuan Earthquake in China , 2010 .

[7]  ChengHu Zhou,et al.  Spatio-temporal rupture process of the 2008 great Wenchuan earthquake , 2009 .

[8]  A Comparison of Recorded Response Spectra from the 2008 Wenchuan, China, Earthquake with Modern Ground-Motion Prediction Models , 2010 .

[9]  J. E. Glynn,et al.  Numerical Recipes: The Art of Scientific Computing , 1989 .

[10]  Jonathan P. Stewart,et al.  Comparison of Ground Motions from Hybrid Simulations to NGA Prediction Equations , 2011 .

[11]  John X. Zhao Geometric Spreading Functions and Modeling of Volcanic Zones for Strong-Motion Attenuation Models Derived from Records in Japan , 2010 .

[12]  A. Pitarka,et al.  Differences in Earthquake Source and Ground Motion Characteristics between Surface and Buried Crustal Earthquakes , 2006 .

[13]  H. Thio,et al.  Attenuation Relations of Strong Ground Motion in Japan Using Site Classification Based on Predominant Period , 2006 .

[14]  Jian Zhang,et al.  An Empirical Site-Classification Method for Strong-Motion Stations in Japan Using h/v Response Spectral Ratio , 2006 .

[15]  Julian J. Bommer,et al.  The Acquisition of Source, Path, and Site Effects from Microearthquake Recordings Using Q Tomography: Application to the United Kingdom , 2008 .

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

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