An Advanced Estimation Algorithm for Ground‐Motion Models with Spatial Correlation

Ground‐motion prediction equations (GMPEs), also called ground‐motion models and attenuation relationships, are empirical models widely used in probabilistic seismic hazard analysis (PSHA). They estimate the conditional distribution of ground shaking at a site given an earthquake of a certain magnitude occurring at a nearby location. In the past decade, the increasing interest in assessing earthquake risk and resilience of spatially distributed portfolios of buildings and infrastructure has motivated the modeling of ground‐motion spatial correlation. This introduces further challenges for researchers to develop statistically rigorous and computationally efficient algorithms to perform ground‐motion model estimation with spatial correlation. To this goal, we introduce a one‐stage ground‐motion estimation algorithm, called the scoring estimation approach, to fit ground‐motion models with spatial correlation. The scoring estimation approach is introduced theoretically and numerically, and it is proven to have desirable properties on convergence and computation. It is a statistically robust method, producing consistent and statistically efficient estimators of inter‐ and intraevent variances and parameters in spatial correlation functions. The performance of the scoring estimation approach is assessed through a comparison with the multistage algorithm proposed by Jayaram and Baker (2010) in a simulation‐based application. The results of the simulation study show that the proposed scoring estimation approach presents comparable or higher accuracy in estimating ground‐motion model parameters, especially when the spatial correlation becomes smoother. The simulation study also shows that ground‐motion models with spatial correlation built via the scoring estimation approach can be used for reliable ground‐shaking intensity predictions. The performance of the scoring estimation approach is further discussed under the ignorance of spatial correlation, and we find that neglecting spatial correlation in ground‐motion models may result in overestimation of interevent variance and underestimation of intraevent variance and thus inaccurate predictions.

[1]  R. Puglia,et al.  Pan-European ground-motion prediction equations for the average horizontal component of PGA, PGV, and 5 %-damped PSA at spectral periods up to 3.0 s using the RESORCE dataset , 2013, Bulletin of Earthquake Engineering.

[2]  Richard F. Gunst,et al.  Applied Regression Analysis , 1999, Technometrics.

[3]  Norman A. Abrahamson,et al.  Summary of the ASK14 Ground Motion Relation for Active Crustal Regions , 2014 .

[4]  Beata Bylina,et al.  The influence of a matrix condition number on iterative methods' convergence , 2011, 2011 Federated Conference on Computer Science and Information Systems (FedCSIS).

[5]  Rodolfo Puglia,et al.  Erratum to: Pan-European ground-motion prediction equations for the average horizontal component of PGA, PGV, and 5 %-damped PSA at spectral periods up to 3.0 s using the RESORCE dataset , 2014, Bulletin of Earthquake Engineering.

[6]  D. Wald,et al.  Topographic Slope as a Proxy for Seismic Site-Conditions (VS30) and Amplification Around the Globe , 2007 .

[7]  Jeffrey M. Woodbridge.  Econometric Analysis of Cross Section and Panel Data , 2002 .

[8]  K. Campbell,et al.  NGA-West2 Ground Motion Model for the Average Horizontal Components of PGA, PGV, and 5% Damped Linear Acceleration Response Spectra , 2014 .

[9]  Yun Wang,et al.  Application of EM Algorithm in Statistics Natural Language Processing , 2013 .

[10]  Simona Esposito,et al.  Spatial Correlation of Spectral Acceleration in European Data , 2012 .

[11]  Steven M. Lalonde,et al.  A First Course in Multivariate Statistics , 1997, Technometrics.

[12]  David R. Brillinger,et al.  Further analysis of the Joyner-Boore attenuation data , 1985 .

[13]  Mario Ordaz,et al.  Multivariate Bayesian Regression Analysis Applied to Ground-Motion Prediction Equations, Part 2: Numerical Example with Actual Data , 2010 .

[14]  Dale L. Zimmerman,et al.  Classical Geostatistical Methods , 2010 .

[15]  Jack W. Baker,et al.  Spatial and Spectral Interpolation of Ground-Motion Intensity Measure Observations , 2018 .

[16]  Gail M. Atkinson,et al.  Probabilistic Characterization of Spatially Correlated Response Spectra for Earthquakes in Japan , 2009 .

[17]  G. Atkinson,et al.  Intraevent Spatial Correlation of Ground-Motion Parameters Using SK-net Data , 2010 .

[18]  Robert R. Youngs,et al.  Update of the Chiou and Youngs NGA Model for the Average Horizontal Component of Peak Ground Motion and Response Spectra , 2014 .

[19]  T. Subba Rao,et al.  Statistics for Spatial Data, Revised Edition, by Noel Cressie. Published by Wiley Classics Library, John Wiley, 2015. Total number of pages: 928. ISBN: 978-1-119-11518-2 , 2016 .

[20]  J. Baker,et al.  Correlation model for spatially distributed ground‐motion intensities , 2009 .

[21]  Jonathan P. Stewart,et al.  NGA-West2 Equations for Predicting PGA, PGV, and 5% Damped PSA for Shallow Crustal Earthquakes , 2014 .

[22]  I. M. Idriss.  An NGA-West2 Empirical Model for Estimating the Horizontal Spectral Values Generated by Shallow Crustal Earthquakes , 2014 .

[23]  Simona Esposito,et al.  PGA and PGV Spatial Correlation Models Based on European Multievent Datasets , 2011 .

[24]  Tsuyoshi Takada,et al.  Macrospatial Correlation Model of Seismic Ground Motions , 2005 .

[25]  W. B. Joyner,et al.  Methods for regression analysis of strong-motion data , 1993, Bulletin of the Seismological Society of America.

[26]  Kaare Brandt Petersen,et al.  The Matrix Cookbook , 2006 .

[27]  Han Ping Hong,et al.  Effect of spatial correlation on estimated ground motion prediction equations , 2009 .

[28]  Brent Kerby.  Semivariogram estimation: Asymptotic theory and applications , 2016 .

[29]  Mike Rees,et al.  5. Statistics for Spatial Data , 1993 .

[30]  Friedemann Wenzel,et al.  Uncertainty and Spatial Correlation of Earthquake Ground Motion in Taiwan , 2010 .

[31]  Carl E. Rasmussen,et al.  Gaussian processes for machine learning , 2005, Adaptive computation and machine learning.

[32]  Yi-Hau Chen,et al.  A New Method for Estimation of the Attenuation Relationship with Variance Components , 2002 .

[33]  Siam Rfview,et al.  CONVERGENCE CONDITIONS FOR ASCENT METHODS , 2016 .

[34]  Roger Woodard,et al.  Interpolation of Spatial Data: Some Theory for Kriging , 1999, Technometrics.

[35]  Noel A Cressie,et al.  On asymptotic distribution and asymptotic efficiency of least squares estimators of spatial variogram parameters , 2002 .

[36]  David R. Brillinger,et al.  An exploratory analysis of the Joyner-Boore attenuation data , 1984 .

[37]  Robert Haining,et al.  Statistics for spatial data: by Noel Cressie, 1991, John Wiley & Sons, New York, 900 p., ISBN 0-471-84336-9, US $89.95 , 1993 .

[38]  Eugene Demidenko,et al.  Mixed Models: Theory and Applications with R , 2013 .

[39]  N. Laird,et al.  Maximum likelihood computations with repeated measures: application of the EM algorithm , 1987 .

[40]  D. Bates,et al.  Nonlinear mixed effects models for repeated measures data. , 1990, Biometrics.

[41]  P. Wolfe.  Convergence Conditions for Ascent Methods. II: Some Corrections , 1971 .

[42]  Mario Ordaz,et al.  Multivariate Bayesian Regression Analysis Applied to Ground-Motion Prediction Equations, Part 1: Theory and Synthetic Example , 2010 .

[43]  D. Wald,et al.  Review Article Topographic Slope as a Proxy for Seismic Site Conditions and Amplification , 2007 .

[44]  Jack W. Baker,et al.  Considering Spatial Correlation in Mixed-Effects Regression and the Impact on Ground-Motion Models , 2010 .

[45]  Philip E. Gill,et al.  Practical optimization , 1981 .

[46]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[47]  Christophe Couvreur,et al.  The EM Algorithm: A Guided Tour , 1997 .

[48]  Gene H. Golub,et al.  Matrix computations , 1983 .

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

[50]  H. Hong,et al.  Spatial correlation of peak ground motions and response spectra , 2008 .

[51]  D. Bates,et al.  Mixed-Effects Models in S and S-PLUS , 2001 .

[52]  J. Bommer,et al.  Empirical Equations for the Prediction of PGA, PGV, and Spectral Accelerations in Europe, the Mediterranean Region, and the Middle East , 2010 .

[53]  Rory A. Fisher,et al.  Theory of Statistical Estimation , 1925, Mathematical Proceedings of the Cambridge Philosophical Society.

[54]  Andrej Pázman,et al.  Nonlinear Regression , 2019, Handbook of Regression Analysis With Applications in R.

[55]  H. Crowley,et al.  Seismic Hazard Assessment (2003–2009) for the Italian Building Code , 2011 .

[56]  P. Wolfe.  Convergence Conditions for Ascent Methods. II , 1969 .