Search

We describe source and path attributes of the Next Generation of Ground-Motion Attenuation for Subduction zones (NGA-Sub) project relational database. The database contains 991 earthquakes between 1937 and 2016 that meet quality assurance standards and have assigned event types (mostly interface or intraslab). Data curation emphasized large magnitudes: 73% of events have M > 5 and 13 events have M > 8. Event attributes in the database can be broadly divided into those related to the moment tensor, those related to rupture surface geometry (referred to as finite-fault parameters), and categorizations by type (e.g. interface, intraslab) and sequence classification (mainshock, aftershock). Earthquake information and moment tensor parameters compiled from earthquake catalogs include event date and origin time, hypocenter location, seismic moment, and the strike, dip, and rake angles of nodal planes. Finite-fault parameters include along-strike length, down-dip width, and depth to top-of-rupture for one or more rectangles used to parameterize rupture surface geometry. These are modified from models in literature where available, and otherwise are simulated using procedures customized for subduction earthquakes. Rupture distance and other metrics are computed using finite-fault representations of sources, and (as applicable) are partitioned into forearc and backarc components. Forearc and backarc regions are defined based on volcano locations, and are assigned to events and sites.

This article summarizes the Next Generation Attenuation (NGA) Subduction (NGA-Sub) project, a major research program to develop a database and ground motion models (GMMs) for subduction regions. A comprehensive database of subduction earthquakes recorded worldwide was developed. The database includes a total of 214,020 individual records from 1,880 subduction events, which is by far the largest database of all the NGA programs. As part of the NGA-Sub program, four GMMs were developed. Three of them are global subduction GMMs with adjustment factors for up to seven worldwide regions: Alaska, Cascadia, Central America and Mexico, Japan, New Zealand, South America, and Taiwan. The fourth GMM is a new Japan-specific model. The GMMs provide median predictions, and the associated aleatory variability, of RotD50 horizontal components of peak ground acceleration, peak ground velocity, and 5%-damped pseudo-spectral acceleration (PSA) at oscillator periods ranging from 0.01 to 10 s. Three GMMs also quantified “within-model” epistemic uncertainty of the median prediction, which is important in regions with sparse ground motion data, such as Cascadia. In addition, a damping scaling model was developed to scale the predicted 5%-damped PSA of horizontal components to other damping ratios ranging from 0.5% to 30%. The NGA-Sub flatfile, which was used for the development of the NGA-Sub GMMs, and the NGA-Sub GMMs coded on various software platforms, have been posted for public use.

The Next-Generation Attenuation for subduction zone regions project (NGA-Sub) has developed data resources and ground motion models for global subduction zone regions. Here we describe the NGA-Sub database. To optimize the efficiency of data storage, access, and updating, data resources for the NGA-Sub project are organized into a relational database consisting of 20 tables containing data, metadata, and computed quantities (e.g. intensity measures, distances). A database schema relates fields in tables to each other through a series of primary and foreign keys. Model developers and other users mostly interact with the data through a flatfile generated as a time-stamped output of the database. We describe the structure of the relational database, the ground motions compiled for the project, and the means by which the data can be accessed. The database contains 71,340 three-component records from 1880 earthquakes from seven global subduction zone regions: Alaska, Central America and Mexico, Cascadia, Japan, New Zealand, South America, and Taiwan. These data were processed on a component-specific basis to minimize noise effects in the data and remove baseline drifts. Provided ground motion intensity measures include peak acceleration, peak velocity, and 5%-damped pseudo-spectral accelerations for a range of oscillator periods.

The Sacramento-San Joaquin Delta in central California is particularly susceptible to damage in a large earthquake due to the vulnerability of the levees that protect cities, farms, and infrastructure. The Delta is located adjacent to the seismically active San Andreas fault system and is also subject to strong ground shaking from numerous other seismic sources in central California, including faults within the Delta. We performed a probabilistic seismic hazard analysis (PSHA) to provide seismic design ground motions for the Metropolitan Water District (MWD) Emergency Freshwater Pathway. We have evaluated the appropriateness of the Next Generation of Attenuation (NGA)-West2 ground motion models for use in our analyses of the Delta, evaluated shear-wave velocity ( VS) data in the vicinity of the Pathway, and performed site response analyses. The latter was performed to compute the probabilistic hazard at the top of the peat at five sites along the Pathway. The sixth site was located outside the Delta and on firm soil. The probabilistic hazard for the six sites and for a range of return periods of engineering relevance were computed in the PSHA. For a return period of 2500 years, the peak horizontal ground acceleration (PGA) values on peat ranged from 0.40 g to 0.53 g. The seismic sources that control the hazard at these sites vary as a function of return period and spectral frequency, but in general, the closer the sites are to faults within the San Andreas fault system, the higher the hazard.

The damaging 4 May 2018 M 6.9 Kalapana earthquake and its aftershocks have provided the largest suite of strong motion records ever produced for an earthquake sequence in Hawaii exceeding the number of records obtained in the deep 2006 M 6.7 Kiholo Bay earthquake. These records provided the best opportunity to understand the processes of strong ground shaking in Hawaii from shallow crustal (20 km) earthquakes. To provide inputs into the stochastic model, we performed an inversion to estimate kappa, stress drops, Ro, and Q(f) using the shallow crustal earthquake database. The GMM is valid from M 4.0 to 8.0 and at Joyner–Boore (RJB) distances up to 400 km. Models were developed for eight VS30 (time-averaged shear-wave velocity in the top 30 m) values corresponding to the National Earthquake Hazards Reduction Program (NEHRP) site bins: A (1500 m/s), B (1080 m/s), B/C (760 m/s), C (530 m/s), C/D (365 m/s), D (260 m/s), D/E (185 m/s), and E (150 m/s). The GMM is for PGA, peak horizontal ground velocity (PGV), and 5%-damped pseudo-spectral acceleration (SA) at 26 periods from 0.01 to 10 s. In addition, we updated our GMM for deep earthquakes (>20 km) to include the same NEHRP site bins using the same approach for the crustal earthquake GMM.

This article documents the earthquake ground motion database developed for the NGA-East Project, initiated as part of the Next Generation Attenuation (NGA) research program and led by the Pacific Earthquake Engineering Research Center (PEER). The project was focused on developing a ground motion characterization model (GMC) model for horizontal ground motions for the large region referred to as Central and Eastern North America (CENA). The CENA region covers most of the U.S. and Canada, from the Rocky Mountains to the Atlantic Ocean and is characterized tectonically as a stable continental region (SCR). The ground-motion database includes the two- and three-component ground-motion recordings from numerous selected events relevant to CENA ( M > 2.5, with distances up to 3500 km) that have been recorded since 1976. The final database contains over 27,000 time series from 82 earthquakes and 1271 recording stations. The ground motion database includes uniformly processed time series, 5% damped pseudo-spectral acceleration (PSA) median-component ordinates for 429 periods ranging from 0.01 to 10 s, duration and Arias intensity in 5% increments, and Fourier amplitude spectra for different time windows. Ground motions and metadata for source, path, and site conditions were subjected to quality checks by topical working groups and the ground-motion model (GMM) developers. The NGA-East database constitutes the largest database of processed recorded ground motions in SRCs and is publicly available from the PEER ground-motion database website.

This article describes the work undertaken within the Next Generation Attenuation (NGA)-East project with the aim of estimating κ0 (the site-specific component of the high-frequency decay parameter, κ) for rock sites in Central and Eastern North America (CENA), using the project’s shallow crustal dataset. We introduce a methodology to address the numerous challenges in CENA: a large dataset in a low-seismicity stable continental region, with poor magnitude and distance coverage, undesirable recording sensor characteristics (low sampling rates leading to poor high-frequency resolution), high uncertainty in the regional stress drop, and lack of site-specific velocity characterization. We use two band-limited κ estimation approaches, the acceleration and displacement spectrum (AS and DS), applied above and below the source corner frequency ( fc), respectively. For band-limited approaches, the key requirement is an estimate of fc, which—apart from the event magnitude readily available in the flatfile—also heavily depends on the highly uncertain stress drop. By considering lower and upper bounds on regional stress drop, we propose a new method to quickly and automatically screen such very large datasets to identify all possible recordings for which band-limited κ approaches can be used. Combining them produces better-quantify estimates of κ and its epistemic uncertainties for this challenging dataset. The mean κ0 values combining the two methods are 13 ± 23 ms for horizontal ground motion.

The Carquinez Bridge geotechnical arrays are operated by the California Strong Motion Instrumentation Program (CSMIP) and recorded a peak ground acceleration (PGA) of approximately 1.0g at ground surface during the 2014 South Napa earthquake. The recorded PGA was significantly larger than those at the nearby surface sites. This study considers surface and downhole recordings from the additional 28 earthquakes recorded at the same arrays to understand the effects of wave propagation and site response at these arrays. Several site response analyses are performed to understand soil nonlinearity using the observed ground accelerations during the 2014 South Napa sequence. Apparent shear wave velocities are calculated from downhole records, which show clear reduction as ground motion intensity increases. Empirical transfer functions (EFTs) are also calculated in which the resonance frequencies became lower during strong shaking during the 2014 South Napa main shock. The in-situ critical damping ratio appears to be frequency dependent in the soft clay deposits. Lower damping at frequencies greater than about 5 Hz may have contributed to the observed PGA at Array #1 during the main shock.

Homogeneous estimates of moment magnitude (M) for small-to-moderate-magnitude (SMM) earthquakes are important to assess regional ground motion compared to global models. However, an estimate of M was available for only 6% of the shallow crustal earthquakes in Iran from 1973 to 2013 in the database developed by the Pacific Earthquake Engineering Research Center. Therefore, this study evaluates the existing conversion models, from ML, Ms, and mb, to M, for M > 4.0 using the Iranian data set. Correlations between the residuals of these conversions were computed. A weighted estimate for M is presented that combines these conversion models using these correlations based on multivariate statistics. The conversion to M was improved with smaller standard deviations, especially when MLand mbwere both available. The inverse-variance–weighted approach, which is common in practice, underestimates the standard deviations because it neglects the positive correlations between these variables.

Until recently, no ground motion prediction model was available for deep ( >20 km) Hawaiian earthquakes, including the 2006 M6.7 Kiholo Bay earthquake. We developed such a model based on the stochastic point-source model. Strong motion data from the 2006 event and 15 other deep Hawaiian earthquakes of M3.3 to M6.2 were inverted using a nonlinear least-squares inversion of Fourier amplitude spectra to estimate stress drops for input into the stochastic modeling and for the few larger events (M ≥ 5.0), to calibrate the ground motion prediction model. The ground motion model is valid for M3.5 to M7.5 over the Joyner-Boore ( RJB) distance range of 20 km to 400 km and are for 5%-damped horizontal spectral acceleration at 27 periods from PGA (0.01 s) to 10.0 s. The shallow site condition assumed for the model is soil and weathered basalt with a mean VS30of 428 m/s.

The NGA-West2 project is a large multidisciplinary, multi-year research program on the Next Generation Attenuation (NGA) models for shallow crustal earthquakes in active tectonic regions. The research project has been coordinated by the Pacific Earthquake Engineering Research Center (PEER), with extensive technical interactions among many individuals and organizations. NGA-West2 addresses several key issues in ground-motion seismic hazard, including updating the NGA database for a magnitude range of 3.0–7.9; updating NGA ground-motion prediction equations (GMPEs) for the “average” horizontal component; scaling response spectra for damping values other than 5%; quantifying the effects of directivity and directionality for horizontal ground motion; resolving discrepancies between the NGA and the National Earthquake Hazards Reduction Program (NEHRP) site amplification factors; analysis of epistemic uncertainty for NGA GMPEs; and developing GMPEs for vertical ground motion. This paper presents an overview of the NGA-West2 research program and its subprojects.

The NGA-West2 site database (SDB) contains information on site condition and instrument housing for 4,147 strong-motion stations with recordings in the project flatfile. The stations are from active tectonic regions, mainly in California, Japan, Taiwan, China, and the Mediterranean area. The principal site parameter is the time-averaged shear wave velocity in the upper 30 m ( V S30), which we characterize using measurements where available (2,013 stations) and proxy-based relationships otherwise. We also provide basin depths from published models for 2,761 sites mostly in California and Japan. We improved the documentation and consistency of site descriptors used as proxies for V S30 estimation (surface geology, ground slope, and geotechnical or geomorphic categories) and analyzed proxy performance relative to V S30 values from measurements. We present protocols for V S30 estimation from proxies that emphasize methods minimizing bias and dispersion relative to data. For each site, we provide the preferred V S30 and its dispersion.

The NGA-West2 project database expands on its predecessor to include worldwide ground motion data recorded from shallow crustal earthquakes in active tectonic regimes post-2000 and a set of small-to-moderate-magnitude earthquakes in California between 1998 and 2011. The database includes 21,336 (mostly) three-component records from 599 events. The parameter space covered by the database is M 3.0 to M 7.9, closest distance of 0.05 to 1,533 km, and site time-averaged shear-wave velocity in the top 30 m of VS30= 94 m/s to 2,100 m/s (although data becomes sparse for distances >400 km and VS30> 1,200 m/s or

A key component of the NGA research project was the development of a strong-motion database with improved quality and content that could be used for ground-motion research as well as for engineering practice. Development of the NGA database was executed through the Lifelines program of the PEER Center with contributions from several research organizations and many individuals in the engineering and seismological communities. Currently, the data set consists of 3551 publicly available multi-component records from 173 shallow crustal earthquakes, ranging in magnitude from 4.2 to 7.9. Each acceleration time series has been corrected and filtered, and pseudo absolute spectral acceleration at multiple damping levels has been computed for each of the 3 components of the acceleration time series. The lowest limit of usable spectral frequency was determined based on the type of filter and the filter corner frequency. For NGA model development, the two horizontal acceleration components were further rotated to form the orientationindependent measure of horizontal ground motion (GMRotI50). In addition to the ground-motion parameters, a large and comprehensive list of metadata characterizing the recording conditions of each record was also developed. NGA data have been systematically checked and reviewed by experts and NGA developers. DOI: 10.1193/1.2894831

The 1984 Halls Valley (Morgan Hill, California) earthquake had a complex seismic source. Velocities of the major seismic phases measured from continuous broadband seismograms at Berkeley Seismographic Station (BKS) and Richmond Field Station (RFS) show unambiguously that the earthquake is predominantly a double event with the second source hypocenter located approximately 17 km southeast of the mainshock hypocenter given by Bolt, Uhrhammer and Darragh (1985). The southeasterly fault rupture of the first source and the location of the focus of the second source have critical implications for the observed spatial variation of the recorded accelerograms. Of particular engineering interest, the high frequency 1.29g pulse of horizontal ground acceleration measured at Coyote Lake dam can be explained primarily as due to the second source and constructive interference of the principal S waves from the two sources.