Quantification of Uncertainties for Conducting Partially Non-ergodic Probabilistic Seismic Hazard Analysis

Abstract

Estimating local site effects and modifying the uncertainty in ground motion predictions are two indispensable parts of partially non-ergodic site-specific PSHA. Local site effects can be estimated using site response simulations or recorded ground motions at the site. When such predictions are available, the aleatory variability of ground motions used in PSHA can be changed to the single station sigma value. However, in these cases, the epistemic uncertainty in predicting site effects must be incorporated into the hazard analyses. This research focuses on the challenges specific to conducting partially non-ergodic site-specific PSHA using recorded ground motions or site response analysis. The main challenge in estimating local site effects using recorded data is whether ground motions collected in a relatively short time can be used to estimate site effects for long return period events. We first develop a database for recorded ground motions at the KiK-net array to investigate this question and use this database to develop a predictive model for the Fourier Amplitude Spectra of ground motions. The ground motion model (GMM) residuals are used to investigate the stability of site terms across different tectonic regimes. We observe that empirical site terms are stable across different tectonic regimes. This observation allows the use of ground motions from any tectonic regime (whether they belong to the tectonic regime that controls the hazard or not) to estimate local site effects. Moreover, in Fourier amplitude, site effects are not dependent on event magnitude and source to site distance; therefore, estimates of site effects from low magnitude events can be easily extrapolated to larger events. The Fourier amplitude GMM developed in this study adds to the library of Fourier amplitude models to be used in future partially non-ergodic site-specific PSHAs. In practice, one of the most common tools for simulating wave propagation is 1-D site response analysis. Two central assumptions in 1-D site response analysis are that the soil profile is comprised of horizontal soil layers of infinite extent and that the vertically propagating SH-waves control the horizontal component of ground motion. SH-waves tend to propagate vertically near the surface because as earthquake waves hit softer layers traveling from the source to the site, they refract until the path becomes steeply inclined. The validity of both assumptions in 1-D site response depends on the geological setting at the site and the geology between the earthquake source and the site, raising the question of which sites are suitable for 1-D site response analysis and what the model error in 1-D site response analysis is. We use the GMM developed for FAS to estimate observed and empirical site terms. The empirical site effects are then compared with the theoretical site effects to determine whether sites are amenable to 1-D site response analyses, and to quantify the model error in the analyses.