Derivation of consistent hard rock (1000 < VS < 3000 m/s) GMPEs from surface and down-hole recordings: analysis of KiK-net data

A key component in seismic hazard assessment is the estimation of ground motion for hard rock sites, either for applications to installations built on this site category, or as an input motion for site response computation. Empirical ground motion prediction equations (GMPEs) are the traditional basis for estimating ground motion while VS30 is the basis to account for site conditions. As current GMPEs are poorly constrained for VS30 larger than 1000 m/s, the presently used approach for estimating hazard on hard rock sites consists of “host-to-target” adjustment techniques based on VS30 and κ0 values. The present study investigates alternative methods on the basis of a KiK-net dataset corresponding to stiff and rocky sites with 500 < VS30 < 1350 m/s. The existence of sensor pairs (one at the surface and one in depth) and the availability of P- and S-wave velocity profiles allow deriving two “virtual” datasets associated to outcropping hard rock sites with VS in the range [1000, 3000] m/s with two independent corrections: 1/down-hole recordings modified from within motion to outcropping motion with a depth correction factor, 2/surface recordings deconvolved from their specific site response derived through 1D simulation. GMPEs with simple functional forms are then developed, including a VS30 site term. They lead to consistent and robust hard-rock motion estimates, which prove to be significantly lower than host-to-target adjustment predictions. The difference can reach a factor up to 3–4 beyond 5 Hz for very hard-rock, but decreases for decreasing frequency until vanishing below 2 Hz.

[1]  P. Renault,et al.  Importance and Impact of Host-to-Target Conversions for Ground Motion Prediction Equations in PSHA , 2012 .

[2]  Adrian Rodriguez-Marek,et al.  Erratum to: Site effect assessment using KiK-net data: Part 1. A simple correction procedure for surface/downhole spectral ratios , 2011, Bulletin of Earthquake Engineering.

[3]  Susan E. Hough,et al.  High-frequency spectra observed at Anza, California: Implications for Q structure , 1988 .

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

[5]  Frank Scherbaum,et al.  On the Discrepancy of Recent European Ground-Motion Observations and Predictions from Empirical Models: Analysis of KiK-net Accelerometric Data and Point-Sources Stochastic Simulations , 2008 .

[6]  B. Kennett,et al.  Reflections, rays, and reverberations , 1974, Bulletin of the Seismological Society of America.

[7]  Norman A. Abrahamson,et al.  Empirical Estimation of High‐Frequency Ground Motion on Hard Rock , 2016 .

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

[9]  Stéphane Drouet,et al.  Analysis of the Origins of κ (Kappa) to Compute Hard Rock to Rock Adjustment Factors for GMPEs , 2011 .

[10]  Fabrice Hollender,et al.  Earthquake Ground Motion in the Mygdonian Basin, Greece: The E2VP Verification and Validation of 3D Numerical Simulation up to 4 Hz , 2015 .

[11]  Pierre-Yves Bard,et al.  The seismic response of two-dimensional sedimentary deposits with large vertical velocity gradients , 1986 .

[12]  R. Borcherdt Effects of local geology on ground motion near San Francisco Bay , 1970 .

[13]  Wei Li,et al.  Site Amplification and Attenuation via Downhole Array Seismogram Inversion: A Comparative Study of the 2003 Miyagi-Oki Aftershock Sequence , 2008 .

[14]  John Douglas,et al.  A κ Model for Mainland France , 2010 .

[15]  Aurore Laurendeau,et al.  Rock and Stiff‐Soil Site Amplification: Dependency on VS30 and Kappa (κ0) , 2013 .

[16]  N. Lam,et al.  Near-surface attenuation modelling based on rock shear-wave velocity profile , 2006 .

[17]  F. D. Martin,et al.  Source‐Related Variability of Site Response in the Mygdonian Basin (Greece) from Accelerometric Recordings and 3D Numerical Simulations , 2017 .

[18]  P. Bard,et al.  Horizontal to vertical spectral ratio and geological conditions: an analysis of strong motion data from Greece and Taiwan (SMART-1) , 1995 .

[19]  Pierre-Yves Bard,et al.  Selecting time windows of seismic phases and noise for engineering seismology applications: a versatile methodology and algorithm , 2018, Bulletin of Earthquake Engineering.

[20]  Y. Fukushima,et al.  Site-dependent attenuation relations of seismic motion parameters at depth using borehole data , 1995, Bulletin of the Seismological Society of America.

[21]  D. Boore The Uses and Limitations of the Square‐Root‐Impedance Method for Computing Site Amplification , 2013 .

[22]  Julian J. Bommer,et al.  Why Do Modern Probabilistic Seismic-Hazard Analyses Often Lead to Increased Hazard Estimates? , 2006 .

[23]  Jamison H. Steidl,et al.  What is a reference site? , 1996, Bulletin of the Seismological Society of America.

[24]  Robert E. Kayen,et al.  A taxonomy of site response complexity , 2012 .

[25]  Frank Scherbaum,et al.  On the Relationship between Fourier and Response Spectra: Implications for the Adjustment of Empirical Ground‐Motion Prediction Equations (GMPEs) , 2016 .

[26]  J. Bommer,et al.  Application of Single-Station Sigma and Site-Response Characterization in a Probabilistic Seismic-Hazard Analysis for a New Nuclear Site , 2014 .

[27]  K. Campbell PREDICTION OF STRONG GROUND MOTION USING THE HYBRID EMPIRICAL METHOD AND ITS USE IN THE DEVELOPMENT OF GROUND-MOTION (ATTENUATION) RELATIONS IN EASTERN NORTH AMERICA , 2003 .

[28]  Kim B. Olsen,et al.  Estimation of Q for Long-Period (>2 sec) Waves in the Los Angeles Basin , 2003 .

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

[30]  Julian J. Bommer,et al.  Criteria for Selecting and Adjusting Ground-Motion Models for Specific Target Regions: Application to Central Europe and Rock Sites , 2006 .

[31]  P. Bard,et al.  Site effect assessment using KiK-net data: Part 1. A simple correction procedure for surface/downhole spectral ratios , 2011, Bulletin of Earthquake Engineering.

[32]  P. Bard,et al.  Site Effect Study in Urban Area: Experimental Results in Grenoble (France) , 2001 .

[33]  F. Scherbaum,et al.  Development of a Response Spectral Ground‐Motion Prediction Equation (GMPE) for Seismic‐Hazard Analysis from Empirical Fourier Spectral and Duration Models , 2015 .

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

[35]  W. Silva,et al.  NGA-West2 Database , 2014 .

[36]  Shigeo Kinoshita,et al.  Kyoshin Net (K-NET) , 1998 .

[37]  Tadahiro Kishida,et al.  Semi-automated procedure for windowing time series and computing Fourier amplitude spectra (FAS) for the NGA-West2 database , 2016 .

[38]  Norman A. Abrahamson,et al.  Taxonomy of κ: A Review of Definitions and Estimation Approaches Targeted to Applications , 2014 .

[39]  C. Beauval,et al.  SITE-SPECIFIC PSHA: COMBINED EFFECTS OF SINGLE STATION SIGMA, HOST-TO-TARGET ADJUSTMENTS AND NON-LINEAR BEHAVIOR , 2017 .

[40]  L. Luzi,et al.  Reference database for seismic ground-motion in Europe (RESORCE) , 2014, Bulletin of Earthquake Engineering.

[41]  T. Ohmachi,et al.  Ground Motion Characteristics Estimated from Spectral Ratio between Horizontal and Verticcl Components of Mietremors. , 1997 .

[42]  W. B. Joyner,et al.  The effect of Quaternary alluvium on strong ground motion in the Coyote Lake, California, earthquake of 1979 , 1981 .

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

[44]  T. Brocher Empirical relations between elastic wavespeeds and density in the Earth's crust , 2005 .

[45]  Linda Al Atik,et al.  Kappa ( ) Scaling of Ground-Motion Prediction Equations Using an Inverse Random Vibration Theory Approach , 2014 .

[46]  Hiroyuki Fujiwara,et al.  Recent Progress of Seismic Observation Networks in Japan , 2004 .

[47]  D. Fäh,et al.  The contribution of scattering to near-surface attenuation , 2017, Journal of Seismology.

[48]  Adrian Rodriguez-Marek,et al.  Analysis of Single-Station Standard Deviation Using the KiK-net Data , 2011 .

[49]  Stéphane Drouet,et al.  vS30, κ, regional attenuation and Mw from accelerograms: application to magnitude 3–5 French earthquakes , 2010 .

[50]  Francisco J. Chávez-García,et al.  Site effect evaluation using spectral ratios with only one station , 1993, Bulletin of the Seismological Society of America.

[51]  Benjamin Edwards,et al.  Reference S‐Wave Velocity Profile and Attenuation Models for Ground‐Motion Prediction Equations: Application to Japan , 2013 .

[52]  Benjamin Edwards,et al.  Attenuation of seismic shear wave energy in Switzerland , 2011 .

[53]  Jamison H. Steidl,et al.  Borehole Response Studies at the Garner Valley Downhole Array, Southern California , 2002 .

[54]  N. Abrahamson,et al.  Understanding the physics of kappa (κ): Insights from a downhole array , 2015 .

[55]  Dino Bindi,et al.  Influence of Soil-Layer Properties on k Evaluation , 2004 .

[56]  Pierre-Yves Bard,et al.  Defining a Standard Rock Site: Propositions Based on the KiK-net Database , 2010 .

[57]  David M. Boore,et al.  Simulation of Ground Motion Using the Stochastic Method , 2003 .