An Earthquake Classification Scheme Adapted for Japan Determined by the Goodness of Fit for Ground‐Motion Prediction Equations

Abstract Establishing a set of ground‐motion prediction equations (GMPEs) for Japan requires earthquake source categories in the dataset. Earthquakes are typically divided into three groups: shallow crustal events that occur in the Earth’s crust, subduction interface events that occur at the interface between the crust or mantle and the subducting plate, and the subduction slab events that occur within the subducting plate. In the present study, we compared the hypocentral locations published in the catalogs of the International Seismological Centre/Engdahl–van der Hilst–Buland (ISC‐EHB; Engdahl et al. , 1998), the Japan Meteorological Agency (JMA), and the National Earthquake Information Center (NEIC) of the U.S. Geological Survey (USGS). The hypocentral location for the same earthquake varies significantly from one catalog to another. We used the subduction interface model from the USGS, Slab 1.0, to help guide the classification. We designed four classification schemes using locations from these three catalogs. We then fitted a set of random effects models to the strong‐motion dataset from these earthquakes to assess the merits of the classification schemes. Our results showed that using ISC‐EHB locations for events before 2005, and then using the preference order of catalogs as (1) JMA locations with high precision levels, (2) ISC‐EHB, and (3) NEIC (excluding the events with a fixed depth) for events since 2005, together with some conditions for subduction interface events, produced the best GMPEs in terms of the maximum log likelihood. We also found that having a separate group for the earthquakes above the subduction interface, but with a depth over 25 km, improved the goodness of fit of the GMPEs. Online Material: Figures of azimuthal‐dependent distance shift of epicenters between catalogs, and between‐event and within‐event residuals.

[1]  N. Umino,et al.  Offshore double-planed shallow seismic zone in the NE Japan forearc region revealed by sP depth phases recorded by regional networks , 2009 .

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

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

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

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

[6]  H. Akaike A new look at the statistical model identification , 1974 .

[7]  E. Engdahl,et al.  Global teleseismic earthquake relocation with improved travel times and procedures for depth determination , 1998, Bulletin of the Seismological Society of America.

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

[9]  Tomoki Watanabe,et al.  Aftershock observation of the 2003 Tokachi-oki earthquake by using dense ocean bottom seismometer network , 2004 .

[10]  Hiroyuki Fujiwara,et al.  A New Attenuation Relation for Strong Ground Motion in Japan Based on Recorded Data , 2006 .

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

[12]  Magnitude‐Scaling Rate in Ground‐Motion Prediction Equations for Response Spectra from Large Subduction Interface Earthquakes in Japan , 2012 .

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

[14]  Hua Xu,et al.  A Comparison of VS30 and Site Period as Site‐Effect Parameters in Response Spectral Ground‐Motion Prediction Equations , 2013 .

[15]  N. A. Abrahamson,et al.  A stable algorithm for regression analyses using the random effects model , 1992, Bulletin of the Seismological Society of America.

[16]  Ying-bin Zhang,et al.  Nonlinear Site Models Derived from 1D Analyses for Ground‐Motion Prediction Equations Using Site Class as the Site Parameter , 2015 .

[17]  John X. Zhao,et al.  Magnitude-Scaling Rate in Ground-Motion Prediction Equations for Response Spectra from Large, Shallow Crustal Earthquakes , 2011 .

[18]  Norman A. Abrahamson,et al.  New Zealand acceleration response spectrum attenuation relations for crustal and subduction zone earthquakes , 2006 .

[19]  David J. Wald,et al.  Slab1.0: A three‐dimensional model of global subduction zone geometries , 2012 .

[20]  W. Silva,et al.  Strong Ground Motion Attenuation Relationships for Subduction Zone Earthquakes , 1997 .

[21]  Y. Fukushima,et al.  Reply to T. Masuda and M. Ohtake's “Comment on ‘A new attenuation relation for peak horizontal acceleration of strong earthquake ground motion in Japan’” , 1990, Bulletin of the Seismological Society of America.

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