Eric M. Thompson

Regional Correlations of VS30 and Velocities Averaged Over Depths Less Than and Greater Than 30 Meters

Using velocity profiles from sites in Japan, California, Turkey, and Europe, we find that the time-averaged shear-wave velocity to 30 m (VS30), used as a proxy for site amplification in recent ground-motion prediction equations (GMPEs) and building codes, is strongly correlated with average velocities to depths less than 30 m (VSz, with z being the averaging depth). The correlations for sites in Japan (corresponding to the KiK-net network) show that VS30 is systematically larger for ag ivenVSz than for profiles from the other regions. The difference largely results from the placement of the KiK-net station locations on rock and rocklike sites, whereas stations in the other regions are generally placed in urban areas underlain by sediments. Using the KiK-net velocity profiles, we provide equations relating VS30 to VSz for z ranging from 5 to 29 m in 1-m increments. These equations (and those for California velocity profiles given in Boore, 2004b) can be used to estimate VS30 from VSz for sites in which velocity profiles do not extend to 30 m. The scatter of the residuals decreases with depth, but, even for an averaging depth of 5 m, a variation in logVS30 of � 1 standard deviation maps into less than a 20% uncertainty in ground motions given by recent GMPEs at short periods. The sensitivity of the ground motions to VS30 uncertainty is considerably larger at long periods (but is less than a factor of 1.2 for averaging depths greater than about 20 m). We also find that VS30 is correlated with VSz for z as great as 400 m for sites of the KiK-net network, providing some justification for using VS30 as a site-response variable for predicting ground motions at periods for which the wavelengths far exceed 30 m. Online Material: Estimates of VS30 at K-NET stations.

Impediments to predicting site response: Seismic property estimation and modeling simplifications

Abstract We compare estimates of the empirical transfer function (ETF) to the plane SH -wave theoretical transfer function (TTF) within a laterally constant medium for invasive and noninvasive estimates of the seismic shear-wave slownesses at 13 Kiban-Kyoshin network stations throughout Japan. The difference between the ETF and either of the TTFs is substantially larger than the difference between the two TTFs computed from different estimates of the seismic properties. We show that the plane SH -wave TTF through a laterally homogeneous medium at vertical incidence inadequately models observed amplifications at most sites for both slowness estimates, obtained via downhole measurements and the spectral analysis of surface waves. Strategies to improve the predictions can be separated into two broad categories: improving the measurement of soil properties and improving the theory that maps the 1D soil profile onto spectral amplification. Using an example site where the 1D plane SH -wave formulation poorly predicts the ETF, we find a more satisfactory fit to the ETF by modeling the full wavefield and incorporating spatially correlated variability of the seismic properties. We conclude that our ability to model the observed site response transfer function is limited largely by the assumptions of the theoretical formulation rather than the uncertainty of the soil property estimates.

Impediments to predicting site response

Abstract We compare estimates of the empirical transfer function (ETF) to the plane SH -wave theoretical transfer function (TTF) within a laterally constant medium for invasive and noninvasive estimates of the seismic shear-wave slownesses at 13 Kiban-Kyoshin network stations throughout Japan. The difference between the ETF and either of the TTFs is substantially larger than the difference between the two TTFs computed from different estimates of the seismic properties. We show that the plane SH -wave TTF through a laterally homogeneous medium at vertical incidence inadequately models observed amplifications at most sites for both slowness estimates, obtained via downhole measurements and the spectral analysis of surface waves. Strategies to improve the predictions can be separated into two broad categories: improving the measurement of soil properties and improving the theory that maps the 1D soil profile onto spectral amplification. Using an example site where the 1D plane SH -wave formulation poorly predicts the ETF, we find a more satisfactory fit to the ETF by modeling the full wavefield and incorporating spatially correlated variability of the seismic properties. We conclude that our ability to model the observed site response transfer function is limited largely by the assumptions of the theoretical formulation rather than the uncertainty of the soil property estimates.

Path Durations for Use in the Stochastic‐Method Simulation of Ground Motions

The stochastic method of ground-motion simulation assumes that the en- ergy in a target spectrum is spread over a duration DT. DT is generally decomposed into the duration due to source effects (DS) and to path effects (DP). For the most commonly used source, seismological theory directly relates DS to the source corner frequency, accounting for the magnitude scaling of DT. In contrast, DP is related to propagation effects that are more difficult to represent by analytic equations based on the physics of the process. We are primarily motivated to revisit DT because the func- tion currently employed by many implementations of the stochastic method for active tectonic regions underpredicts observed durations, leading to an overprediction of ground motions for a given target spectrum. Further, there is some inconsistency in the literature regarding which empirical duration corresponds to DT. Thus, we begin by clarifying the relationship between empirical durations and DT as used in the first authors implementation of the stochastic method, and then we develop a new DP relationship. The new DP function gives significantly longer durations than in the previous DP function, but the relative contribution of DP to DT still diminishes with increasing magnitude. Thus, this correction is more important for small events or subfaults of larger events modeled with the stochastic finite-fault method.

Geotechnical Effects of the 2015 Magnitude 7.8 Gorkha, Nepal, Earthquake and Aftershocks

This article summarizes the geotechnical effects of the 25 April 2015 M 7.8 Gorkha, Nepal, earthquake and aftershocks, as documented by a reconnaissance team that undertook a broad engineering and scientific assessment of the damage and collected perishable data for future analysis. Brief descriptions are provided of ground shaking, surface fault rupture, landsliding, soil failure, and infrastructure performance. The goal of this reconnaissance effort, led by Geotechnical Extreme Events Reconnaissance, is to learn from earthquakes and mitigate hazards in future earthquakes.

Revisions to Some Parameters Used in Stochastic-Method Simulations of Ground Motion

Abstract The stochastic method of ground‐motion simulation specifies the amplitude spectrum as a function of magnitude ( M ) and distance ( R ). The manner in which the amplitude spectrum varies with M and R depends on physical‐based parameters that are often constrained by recorded motions for a particular region (e.g., stress parameter, geometrical spreading, quality factor, and crustal amplifications), which we refer to as the seismological model. The remaining ingredient for the stochastic method is the ground‐motion duration. Although the duration obviously affects the character of the ground motion in the time domain, it also significantly affects the response of a single‐degree‐of‐freedom oscillator. Recently published updates to the stochastic method include a new generalized double‐corner‐frequency source model, a new finite‐fault correction, a new parameterization of duration, and a new duration model for active crustal regions. In this article, we augment these updates with a new crustal amplification model and a new duration model for stable continental regions. Random‐vibration theory (RVT) provides a computationally efficient method to compute the peak oscillator response directly from the ground‐motion amplitude spectrum and duration. Because the correction factor used to account for the nonstationarity of the ground motion depends on the ground‐motion amplitude spectrum and duration, we also present new RVT correction factors for both active and stable regions. Online Material: Files of coefficients for evaluating distance ( D rms ), time‐domain–to–random‐vibration ratios, and SMSIM parameter files.

Empirical improvements for estimating earthquake response spectra with random‐vibration theory

The stochastic method of ground-motion simulation is often used in combination with the random-vibration theory to directly compute ground-motion in- tensity measures, thereby bypassing the more computationally intensive time-domain simulations. Key to the application of random-vibration theory to simulate response spectra is determining the duration (Drms) used in computing the root-mean-square oscillator response. Boore and Joyner (1984) originally proposed an equation for Drms, which was improved upon by Liu and Pezeshk (1999). Though these equations are both substantial improvements over using the duration of the ground-motion ex- citation for Drms, we document systematic differences between the ground-motion intensity measures derived from the random-vibration and time-domain methods for both of these Drms equations. These differences are generally less than 10% for most magnitudes, distances, and periods of engineering interest. Given the systematic nature of the differences, however, we feel that improved equations are warranted. We empirically derive new equations from time-domain simulations for eastern and western North America seismological models. The new equations improve the random-vibration simulations over a wide range of magnitudes, distances, and oscil- lator periods. Online Material: SMSIM parameter files, tables of coefficients and model parameters, and shaded contour plots of TD/RV ratios for two WNA models.

Geotechnical Reconnaissance of the 2002 Denali Fault, Alaska, Earthquake

The 2002 M7.9 Denali fault earthquake resulted in 340 km of ruptures along three separate faults, causing widespread liquefaction in the fluvial deposits of the alpine valleys of the Alaska Range and eastern lowlands of the Tanana River. Areas affected by liquefaction are largely confined to Holocene alluvial deposits, man-made embankments, and backfills. Liquefaction damage, sparse surrounding the fault rupture in the western region, was abundant and severe on the eastern rivers: the Robertson, Slana, Tok, Chisana, Nabesna and Tanana Rivers. Synthetic seismograms from a kinematic source model suggest that the eastern region of the rupture zone had elevated strong-motion levels due to rupture directivity, supporting observations of elevated geotechnical damage. We use augered soil samples and shear-wave velocity profiles made with a portable apparatus for the spectral analysis of surface waves (SASW) to characterize soil properties and stiffness at liquefaction sites and three trans-Alaska pipeline pump station accelerometer locations.

A Geospatial Liquefaction Model for Rapid Response and Loss Estimation

We describe an approach to model liquefaction extent that focuses on identifying broadly available geospatial variables (e.g., derived from digital elevation models) and earthquake-specific parameters (e.g., peak ground acceleration, PGA). A key step is database development: We focus on the 1995 Kobe and 2010–2011 Christchurch earthquakes because the presence/absence of liquefaction has been mapped so that the database is unbiased with respect to the areal extent of liquefaction. We derive two liquefaction models with explanatory variables that include PGA, shear-wave velocity, compound topographic index, and a newly defined normalized distance parameter (distance to coast divided by the sum of distance to coast and distance to the basin inland edge). To check the portability/reliability of these models, we apply them to the 2010 Haiti earthquake. We conclude that these models provide first-order approximations of the extent of liquefaction, appropriate for use in rapid response, loss estimation, and simulations.

A Flatfile of Ground Motion Intensity Measurements from Induced Earthquakes in Oklahoma and Kansas

We have produced a uniformly processed database of orientation-independent (RotD50, RotD100) ground motion intensity measurements containing peak horizontal ground motions (accelerations and velocities) and 5-percent-damped pseudospectral accelerations (0.1–10 s) from more than 3,800 M ≥ 3 earthquakes in Oklahoma and Kansas that occurred between January 2009 and December 2016. Ground motion time series were collected from regional, national, and temporary seismic arrays out to 500 km. We relocated the majority of the earthquake hypocenters using a multiple-event relocation algorithm to produce a set of near uniformly processed hypocentral locations. Ground motion processing followed standard methods, with the primary objective of reducing the effects of noise on the measurements. Regional wave propagation features and the high seismicity rate required careful selection of signal windows to ensure that we captured the entire ground motion record and that contaminating signals from extraneous earthquakes did not contribute to the database. Processing was carried out with an automated scheme and resulted in a database comprising more than 174,000 records (https://dx.doi.org/10.5066/F73B5X8N). We anticipate that these results will be useful for improved understanding of earthquake ground motions and for seismic hazard applications.