USGS ScienceSearch

Geology topics

Norman A. Abrahamson

Publications and source records attributed to Norman A. Abrahamson.

16 recordsLinked to original sources

Why do seismic hazard models worldwide appear to overpredict historical intensity observations?

Probabilistic seismic hazard assessments (PSHAs) provide the scientific basis for building codes to reduce damage from earthquakes. Despite their substantial impact, little is known about how well PSHA predicts actual shaking. Recent PSHA for California, Japan, Italy, Nepal, and France appear to consistently overpredict historically observed earthquake shaking intensities. Numerical simulations show that observed shaking is equally likely to be above or below predictions. This result from independently developed models and datasets in different countries and tectonic settings indicates possible systematic bias in the hazard models, the observations, or both. Analysis of possible causes shows that much of the discrepancy is due to a subtle and rarely considered issue: the conversion equations used in comparing the models—which forecast shaking as peak ground acceleration or velocity—and observations—parameterizations of qualitative shaking reports. Historical shaking reports fill a crucial data gap, but more research is warranted on how qualitative observations relate to instrumental shaking measures for earthquakes.

Science Advances

Uncertainties in intensity-based earthquake magnitude estimates

Estimating the magnitude of historical earthquakes is crucial for assessing seismic hazard. Magnitudes of early‐instrumental earthquakes can be inferred using a combination of instrumental records, field observations, and the observed distribution of shaking intensity determined from macroseismic observations. For earthquakes before 1900, shaking intensity distributions often provide the only information to constrain earthquake magnitude. Considerable effort has been made to develop methods to estimate the magnitude of moderate‐to‐large historical earthquakes using shaking intensities derived from macroseismic data. In this study, we consider earthquakes in California with known instrumental magnitudes to explore uncertainties in estimating the magnitude of historical earthquakes from intensity information alone. We use three California‐specific intensity prediction equations (IPEs) and an IPE based on a global ground‐motion model (GMM) to determine optimum intensity‐based magnitudes for 33 moderate‐to‐large California earthquakes between 1979 and 2021. Intensity‐based magnitudes are close to instrumental magnitudes on average. However, intensity‐based magnitudes for individual events differ by as much as 2.2 magnitude units from instrumental magnitudes. This result reflects the weak dependence of ground motions and shaking intensities on moment magnitude and their strong dependence on stress drop. Considering the intensity distributions of the 1906 San Francisco and 1989 Loma Prieta earthquakes, we show that information that could constrain rupture length is discarded when considering only the 2D decay of intensity with distance. We also show that ground‐motion intensity conversion equations used in a GMM‐based approach may cause a systematic overestimation of large historical earthquake magnitudes. This study underscores both the reducible and potentially irreducible uncertainties associated with using intensity data to estimate magnitudes of historical earthquakes using IPEs and highlights the value of using additional information to constrain rupture dimensions. Using intensity observations alone, moment magnitude uncertainties are typically on the order of a full unit.

California

Panel review of Ground Motion Characterization Model in 2023 NSHM

The 2023 National Seismic Hazard Model (NSHM; Petersen et al., 2023) has two major components – a seismic source characterization (SSC) model and a ground motion characterization (GMC) model. The US Geological Survey (USGS) established separate panels to review and provide input on these two models. Both panels are advisory, meaning that they provide input on technical issues for consideration by the USGS NSHM team, but they do not have decision making authority. Here, we report on the activities and recommendations of the Ground Motion Characterization Panel, made up of the authors of this review. Final modeling decisions are presented in separate USGS documents, including Petersen et al., (2023) and Moschetti et al., (2023). Where modeling decisions depart from our recommendations, the rationale is explained in those publications.

Report

Estimation of site terms in ground-motion models for California using horizontal-to-vertical spectral ratios from microtremor

The horizontal‐to‐vertical spectral ratios from microtremor (mHVSR) data obtained at 196 seismic stations in California are used to evaluate three alternative microtremor‐based proxies for site amplification for use in ground‐motion models (GMMs): the site fundamental period ( ⁠ f 0 "> f 0 ⁠ ), the period‐dependent amplitude of the mHVSR( T ), and the normalized amplitude of the mHVSR( T ). The alternative parameters are evaluated for the sites with and without measurements of V S 30 "> V S 30 ⁠ . If a V S 30 "> V S 30 measurement is not available for a site, then f 0 "> f 0 has the highest correlation with the site amplification for short periods ( T <1 s) and the normalized amplitude of the mHVSR( T ) has the highest correlation for long periods ( T ≥1 s). If a measurement of the V S 30 "> V S 30 is available for a site, then the normalized amplitude of the mHVSR( T ) has the highest correlation for the site amplification not explained by V S 30 "> V S 30 for all periods. For both cases, the correlations are strongest at the longer periods as mHVSR( T ) measurements excel at providing valuable information for sites with long‐period amplification due to the deeper velocity structure. In particular, for sites with a V S 30 "> V S 30 measurement, the normalized mHVSR( T ) amplitude provides more information about the long‐period site terms than the basin depth currently used in GMMs. Empirical models of the median and standard deviation of the site terms based on the normalized mHVSR( T ) curves are developed for the two cases. These models can be used directly in the ASK14 GMM to modify the median and aleatory standard deviation or they can be used to estimate the site‐specific site term in the context of a partially nonergodic GMM. Including the mHVSR( T ) measurement can have a significant effect on estimates of the ground motion at a site: the range 5%–95% on the observed HVSR( T ) values corresponds to factors of 0.6–1.6 for the median spectral acceleration for periods between 0.5 and 4 s.

California

A study on the effect of site response on California seismic hazard map assessment

Prior studies have repeatedly shown that probabilistic seismic hazard maps from several different countries predict higher shaking than that observed. Previous map assessments have not, however, considered the influence of site response on hazard. Seismologists have long acknowledged the influence of near-surface geology, in particular low-impedance sediment layers, on earthquake ground-motion at frequencies of engineering concern. Although the overall effects of site response are complex, modern ground-motion models (GMMs) account for site effects using terms based on V S30 , the time-averaged shear-wave velocity in the upper 30 m of the Earth’s surface. In this study, we consider general implications of incorporating site terms from modern GMMs using site-specific V S30 as a proxy in probabilistic seismic hazard maps for California. At the long periods (1–5 s) that affect tall buildings, site terms amplify the mapped hazard by factors of 1–3 at many sites relative to maps calculated for the standard reference soft-rock site condition, V S30 = 760 m/s. However, at the short periods of ground-motion that are the main contributors to peak ground acceleration (PGA) and thus affect smaller structures, only negligible effects occur due to nonlinear deamplification of strong ground-motion at high frequencies. Nonlinear deamplification increases as the shaking level increases. For very strong shaking, deamplification can overcome the linear amplification, yielding net deamplification. We explore the implications of these results for the evaluation of hazard maps. Because site effects do not change the maps appreciably at short periods, we can exclude site response as an explanation for why the maps overpredict historically observed shaking as captured by the California Historical Intensity Mapping Project (CHIMP) dataset. The results are expected to be generalizable to regions that are comparable to California in terms of structure and seismicity rates. In low-to-moderate-seismicity regions where the hazard reflects weaker shaking, nonlinear site response is expected to be less important for the hazard.

California

Comparisons of the NGA-Subduction ground motion models

In this article, ground-motion models (GMMs) for subduction earthquakes recently developed as part of the Next Generation Attenuation-Subduction (NGA-Sub) project are compared. The four models presented in this comparison study are documented in their respective articles submitted along with this article. Each of these four models is based on the analysis of the large NGA-Sub database. Three of the four current models are developed for a global version as well as separate regionalized models. The fourth model was developed based on earthquakes only from Japan, and as such is applicable only for Japan. As part of this comparison study, a general discussion on the parameterization of the four models and the regionalization of the three models is provided. The specific strengths and or weaknesses or the technical decisions and justifications of any one model are not part of this comparison. A selected suite of deterministic attenuation curves and spectra are presented for the models along with a selected suite of currently used subduction models. A limited number of comparisons are presented in this article with a larger number of comparisons and the digital values provided in the electronic attachment. In addition to these scenario calculation comparisons, the results from a standard probabilistic seismic hazard analysis (PSHA) for two sites located in the Pacific Northwest Region in the state of Washington are presented. These calculations highlight the potential impact of using the new GMMs. Based on the comparisons presented here, a general understanding of these new GMMs can be obtained with the expectation that the implementation of a specific seismic hazard study should incorporate similar and additional comparisons and sensitivity studies pertinent to the site of interest.

Earthquake Spectra

NGA-Subduction research program

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.

Earthquake Spectra

Selection of random vibration theory procedures for the NGA-East project and ground-motion modeling

Traditional ground-motion models (GMMs) are used to compute pseudo-spectral acceleration (PSA) from future earthquakes and are generally developed by regression of PSA using a physics-based functional form. PSA is a relatively simple metric that correlates well with the response of several engineering systems and is a metric commonly used in engineering evaluations; however, characteristics of the PSA calculation make application of scaling factors dependent on the frequency content of the input motion, complicating the development and adaptability of GMMs. By comparison, Fourier amplitude spectrum (FAS) represents ground-motion amplitudes that are completely independent from the amplitudes at other frequencies, making them an attractive alternative for GMM development. Random vibration theory (RVT) predicts the peak response of motion in the time domain based on the FAS and a duration, and thus can be used to relate FAS to PSA. Using RVT to compute the expected peak response in the time domain for given FAS therefore presents a significant advantage that is gaining traction in the GMM field. This article provides recommended RVT procedures relevant to GMM development, which were developed for the Next Generation Attenuation (NGA)-East project. In addition, an orientation-independent FAS metric—called the effective amplitude spectrum (EAS)—is developed for use in conjunction with RVT to preserve the mean power of the corresponding two horizontal components considered in traditional PSA-based modeling (i.e., RotD50). The EAS uses a standardized smoothing approach to provide a practical representation of the FAS for ground-motion modeling, while minimizing the impact on the four RVT properties ( zeroth moment, m0 ; bandwidth parameter, δ ; frequency of zero crossings, fz ; and frequency of extrema, fe ). Although the recommendations were originally developed for NGA-East, they and the methodology they are based on can be adapted to become portable to other GMM and engineering problems requiring the computation of PSA from FAS.

Earthquake Spectra

Earthquake stress drop and Arias Intensity

We determine earthquake stress drops directly from the Arias intensity database of NGA-West2. Arias intensity (Arias, 1970) is an engineering measure proportional to the integral of the absolute value of acceleration squared, over the significant duration of the signal. As such, it is closely related to root-mean-square acceleration, and can readily be connected to earthquake stress drop (Hanks and McGuire, 1981). Arias intensity records out to 100 km yield stable stress drops for moderate-to-large magnitude earthquakes, M6.5+; for smaller events ~M4.5 – 6.5, only closer-in records yield stable results. For the 116 events considered, stress drops are about 35% larger for Class 1 mainshocks than for traditional on-fault Class 2 aftershocks, and smaller for those aftershocks close to the main fault plane. Aftershock stress drops show large variability, however, implying that on average they re-rupture weakened patches, but can also rupture intact rock or high-stress asperities. We observe an increase of stress drop with earthquake depth similar to that of other studies but do not find any significant faulting mechanism dependence. The variability of the Arias intensity-based stress drop is lower than that of eGf-based stress drops from Baltay et al. (2010, 2011), and nearly on par with variability seen in ground-motion prediction equations. The Arias intensity stress drop is a novel and promising method to estimate stress drop without the need for path and site corrections, and yields further insight into the connection between source physics and ground-motion.

Journal of Geophysical Research B: Solid Earth

A generalization of the double-corner-frequency source spectral model and its use in the SCEC BBP validation exercise

The stochastic method of simulating ground motions requires the specification of the shape and scaling with magnitude of the source spectrum. The spectral models commonly used are either single-corner-frequency or double-corner-frequency models, but the latter have no flexibility to vary the high-frequency spectral levels for a specified seismic moment. Two generalized double-corner-frequency ω 2 source spectral models are introduced, one in which two spectra are multiplied together, and another where they are added. Both models have a low-frequency dependence controlled by the seismic moment, and a high-frequency spectral level controlled by the seismic moment and a stress parameter. A wide range of spectral shapes can be obtained from these generalized spectral models, which makes them suitable for inversions of data to obtain spectral models that can be used in ground-motion simulations in situations where adequate data are not available for purely empirical determinations of ground motions, as in stable continental regions. As an example of the use of the generalized source spectral models, data from up to 40 stations from seven events, plus response spectra at two distances and two magnitudes from recent ground-motion prediction equations, were inverted to obtain the parameters controlling the spectral shapes, as well as a finite-fault factor that is used in point-source, stochastic-method simulations of ground motion. The fits to the data are comparable to or even better than those from finite-fault simulations, even for sites close to large earthquakes.

Bulletin of the Seismological Society of America

NGA-West2 Research Project

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&ndash;7.9; updating NGA ground-motion prediction equations (GMPEs) for the &ldquo;average&rdquo; 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.

Earthquake Spectra

Adding fling effects to processed ground‐motion time histories

Fling is the engineering term for the effects of the permanent tectonic offset, caused by a rupturing fault in the recorded ground motions near the fault. It is expressed by a one‐sided pulse in ground velocity and a nonzero final displacement at the end of shaking. Standard processing of earthquake time histories removes some of the fling effects that may be required for engineering applications. A method to parameterize the fling‐step time history and to superimpose it onto traditionally processed time histories has been developed by Abrahamson (2002) . In this paper, we first present an update to the Abrahamson (2002) fling‐step models, in which the fling step is parameterized as a single cycle of a sine wave. Parametric models are presented for the sine‐wave amplitude ( D site ) and period ( T f ). The expressions for D site and T f are derived from an extensive set of finite‐fault simulations conducted on the Southern California Earthquake Center broadband platform (see Data and Resources). The simulations were run with the Graves and Pitarka (2010) hybrid simulation method and included strike‐slip and reverse scenarios for magnitudes of 6.0–8.2 and dips of 30 through 90. Next, an improved approach for developing design ground motions with fling effects is presented, which deals with the problem of double‐counting intermediate period components that were not removed by the standard ground‐motion processing. Finally, the results are validated against a set of 84 empirical recordings containing fling.

Bulletin of the Seismological Society of America

Strong ground motions generated by earthquakes on creeping faults

A tenet of earthquake science is that faults are locked in position until they abruptly slip during the sudden strain-relieving events that are earthquakes. Whereas it is expected that locked faults when they finally do slip will produce noticeable ground shaking, what is uncertain is how the ground shakes during earthquakes on creeping faults. Creeping faults are rare throughout much of the Earth's continental crust, but there is a group of them in the San Andreas fault system. Here we evaluate the strongest ground motions from the largest well-recorded earthquakes on creeping faults. We find that the peak ground motions generated by the creeping fault earthquakes are similar to the peak ground motions generated by earthquakes on locked faults. Our findings imply that buildings near creeping faults need to be designed to withstand the same level of shaking as those constructed near locked faults.

Geophysical Research Letters

Engineering uses of physics-based ground motion simulations

This paper summarizes validation methodologies focused on enabling ground motion simulations to be used with confidence in engineering applications such as seismic hazard analysis and dynmaic analysis of structural and geotechnical systems. Numberical simullation of ground motion from large erthquakes, utilizing physics-based models of earthquake rupture and wave propagation, is an area of active research in the earth science community. Refinement and validatoin of these models require collaboration between earthquake scientists and engineering users, and testing/rating methodolgies for simulated ground motions to be used with confidence in engineering applications. This paper provides an introduction to this field and an overview of current research activities being coordinated by the Souther California Earthquake Center (SCEC). These activities are related both to advancing the science and computational infrastructure needed to produce ground motion simulations, as well as to engineering validation procedures. Current research areas and anticipated future achievements are also discussed.

Conference Paper

Extreme ground motions and Yucca Mountain

Yucca Mountain is the designated site of the underground repository for the United States' high-level radioactive waste (HLW), consisting of commercial and military spent nuclear fuel, HLW derived from reprocessing of uranium and plutonium, surplus plutonium, and other nuclear-weapons materials. Yucca Mountain straddles the western boundary of the Nevada Test Site, where the United States has tested nuclear devices since the 1950s, and is situated in an arid, remote, and thinly populated region of Nevada, ~100 miles northwest of Las Vegas. Yucca Mountain was originally considered as a potential underground repository of HLW because of its thick units of unsaturated rocks, with the repository horizon being not only ~300 m above the water table but also ~300 m below the Yucca Mountain crest. The fundamental rationale for a geologic (underground) repository for HLW is to securely isolate these materials from the environment and its inhabitants to the greatest extent possible and for very long periods of time. Given the present climate conditions and what is known about the current hydrologic system and conditions around and in the mountain itself, one would anticipate that the rates of infiltration, corrosion, and transport would be very low—except for the possibility that repository integrity might be compromised by low-probability disruptive events, which include earthquakes, strong ground motion, and (or) a repository-piercing volcanic intrusion/eruption. Extreme ground motions (ExGM), as we use the phrase in this report, refer to the extremely large amplitudes of earthquake ground motion that arise at extremely low probabilities of exceedance (hazard). They first came to our attention when the 1998 probabilistic seismic hazard analysis for Yucca Mountain was extended to a hazard level of 10 -8 /yr (a 10-4/yr probability for a 10 4 -year repository “lifetime”). The primary purpose of this report is to summarize the principal results of the ExGM research program as they have developed over the past 5 years; what follows will be focused on Yucca Mountain, but not restricted to it.

Nevada

Representation of bidirectional ground motions for design spectra in building codes

The 2009 NEHRP Provisions modified the definition of horizontal ground motion from the geometric mean of spectral accelerations for two components to the peak response of a single lumped mass oscillator regardless of direction. These maximum-direction (MD) ground motions operate under the assumption that the dynamic properties of the structure (e.g., stiffness, strength) are identical in all directions. This assumption may be true for some in-plan symmetric structures, however, the response of most structures is dominated by modes of vibration along specific axes (e.g., longitudinal and transverse axes in a building), and often the dynamic properties (especially stiffness) along those axes are distinct. In order to achieve structural designs consistent with the collapse risk level given in the NEHRP documents, we argue that design spectra should be compatible with expected levels of ground motion along those principal response axes. The use of MD ground motions effectively assumes that the azimuth of maximum ground motion coincides with the directions of principal structural response. Because this is unlikely, design ground motions have lower probability of occurrence than intended, with significant societal costs. We recommend adjustments to make design ground motions compatible with target risk levels.

Earthquake Spectra