Toward Physics-Based Nonergodic PSHA: A Prototype Fully Deterministic Seismic Hazard Model for Southern California

Author(s):  
Kevin R. Milner ◽  
Bruce E. Shaw ◽  
Christine A. Goulet ◽  
Keith B. Richards-Dinger ◽  
Scott Callaghan ◽  
...  

ABSTRACT We present a nonergodic framework for probabilistic seismic-hazard analysis (PSHA) that is constructed entirely of deterministic, physical models. The use of deterministic ground-motion simulations in PSHA calculations is not new (e.g., CyberShake), but prior studies relied on kinematic rupture generators to extend empirical earthquake rupture forecasts. Fully dynamic models, which simulate rupture nucleation and propagation of static and dynamic stresses, are still computationally intractable for the large simulation domains and many seismic cycles required to perform PSHA. Instead, we employ the Rate-State earthquake simulator (RSQSim) to efficiently simulate hundreds of thousands of years of M≥6.5 earthquake sequences on the California fault system. RSQSim produces full slip-time histories for each rupture, which, unlike kinematic models, emerge from frictional properties, fault geometry, and stress transfer; all intrinsic variability is deterministic. We use these slip-time histories directly as input to a 3D wave-propagation code within the CyberShake platform to obtain simulated Fmax=0.5  Hz ground motions. The resulting 3 s spectral acceleration ground motions closely match empirical ground-motion model (GMM) estimates of median and variability of shaking. When computed over a range of sources and sites, the variability is similar to that of ergodic GMMs. Variability is reduced for individual pairs of sources and sites that repeatedly sample a single path, which is expected for a nonergodic model. This results in increased exceedance probabilities for certain characteristic ground motions for a source–site pair, while decreasing probabilities at the extreme tails of the ergodic GMM predictions. We present these comparisons and preliminary fully deterministic physics-based RSQSim–CyberShake hazard curves, as well as a new technique for estimating within- and between-event variability through simulation.

1999 ◽  
Vol 89 (2) ◽  
pp. 501-520 ◽  
Author(s):  
Paolo Bazzurro ◽  
C. Allin Cornell

Abstract Probabilistic seismic hazard analysis (PSHA) integrates over all potential earthquake occurrences and ground motions to estimate the mean frequency of exceedance of any given spectral acceleration at the site. For improved communication and insights, it is becoming common practice to display the relative contributions to that hazard from the range of values of magnitude, M, distance, R, and epsilon, ɛ, the number of standard deviations from the median ground motion as predicted by an attenuation equation. The proposed disaggregation procedures, while conceptually similar, differ in several important points that are often not reported by the researchers and not appreciated by the users. We discuss here such issues, for example, definition of the probability distribution to be disaggregated, different disaggregation techniques, disaggregation of R versus ln R, and the effects of different binning strategies on the results. Misconception of these details may lead to unintended interpretations of the relative contributions to hazard. Finally, we propose to improve the disaggregation process by displaying hazard contributions in terms of not R, but latitude, longitude, as well as M and ɛ. This permits a display directly on a typical map of the faults of the surrounding area and hence enables one to identify hazard-dominating scenario events and to associate them with one or more specific faults, rather than a given distance. This information makes it possible to account for other seismic source characteristics, such as rupture mechanism and near-source effects, during selection of scenario-based ground-motion time histories for structural analysis.


1987 ◽  
Vol 77 (4) ◽  
pp. 1110-1126
Author(s):  
Anne S. Kiremidjian ◽  
Shigeru Suzuki

Abstract A stochastic model is presented for estimating probabilities of exceeding site ground motions due to temporally dependent earthquake events. The model reflects the hypothesized dependence of the size of large earthquake events on the time of occurrence of the last major earthquake. An empirical attenuation relationship is used to describe the ground motion at a site originating from a well-defined fault system. The application of the model to the Middle America Trench is discussed. The seismic hazard potential in Mexico City is computed in terms of probabilities of exceeding peak ground acceleration levels. The results indicate that consideration of the seismic gap is important for estimating the seismic hazard at a site. It is also observed that site hazard estimates are greatly dependent on the specific attenuation relationship used. The need for other approaches of ground motion estimation is recognized.


2019 ◽  
Vol 35 (2) ◽  
pp. 759-786 ◽  
Author(s):  
Karim Tarbali ◽  
Brendon A. Bradley ◽  
Jack W. Baker

This paper focuses on the selection of ground motions for seismic response analysis in the near-fault region, where directivity effects are significant. An approach is presented to consider forward directivity velocity pulse effects in seismic hazard analysis without separate hazard calculations for ‘pulse-like’ and ‘non-pulse-like’ ground motions, resulting in a single target hazard (at the site of interest) for ground motion selection. The ability of ground motion selection methods to appropriately select records that exhibit pulse-like ground motions in the near-fault region is then examined. Applications for scenario and probabilistic seismic hazard analysis cases are examined through the computation of conditional seismic demand distributions and the seismic demand hazard. It is shown that ground motion selection based on an appropriate set of intensity measures (IMs) will lead to ground motion ensembles with an appropriate representation of the directivity-included target hazard in terms of IMs, which are themselves affected by directivity pulse effects. This alleviates the need to specify the proportion of pulse-like motions and their pulse periods a priori as strict criteria for ground motion selection.


1998 ◽  
Vol 25 (2) ◽  
pp. 305-318 ◽  
Author(s):  
Gail M Atkinson ◽  
Igor A Beresnev

Ground-motion time histories which are compatible with the uniform hazard spectra (UHS) provided by the new national seismic hazard maps of the Geological Survey of Canada (GSC) are simulated. Time histories are simulated for the following cities: Halifax, La Malbaie, Québec, Montreal, Ottawa, Toronto, Prince George, Tofino, Vancouver, and Victoria. The target UHS for the time history simulations are the GSC 5% damped horizontal-component spectra for "firm ground" (Class B) sites for an annual probability of 1/500. The Canadian Council on Earthquake Engineering is currently considering the adoption of these maps as the seismological basis for the earthquake design requirements for future editions of the National Building Code of Canada. It is therefore useful to have compatible time histories for these spectra, in order that dynamic analysis methods requiring the use of time histories can be employed. The simulated records provide a realistic representation of ground motion for the earthquake magnitudes and distances that contribute most strongly to hazard at the selected cities and probability level. For each selected city, two horizontal components are generated for a moderate earthquake nearby, and two horizontal components are generated for a larger earthquake farther away. These records match the short- and long-period ends of the target UHS, respectively. These simulations for local and regional crustal earthquakes are based on a point-source stochastic simulation procedure. For cities in British Columbia, records are also simulated for a scenario M8.5 earthquake on the Cascadia subduction zone, using a stochastic finite-fault simulation model. Four different rupture scenarios are considered. The ground motions for this scenario event are not associated with a specific probability level, but current information suggests that their probability of occurrence is comparable to that of the 1/500 UHS (the probabilistic analyses performed for the national hazard maps do not explicitly include the Cascadia subduction event). Thus it would be reasonable to conduct engineering analyses for cities in British Columbia using both the simulated crustal-event motions and the simulated Cascadia-event motions for the Cascadia event. The time histories simulated for this study are available free of charge to all interested parties.Key words: compatible time-histories, seismic hazard, ground motions.


2007 ◽  
Vol 23 (4) ◽  
pp. 867-891 ◽  
Author(s):  
Polsak Tothong ◽  
C. Allin Cornell ◽  
J. W. Baker

Probabilistic seismic hazard analysis (PSHA) is widely used to estimate the ground motion intensity that should be considered when assessing a structure's performance. Disaggregation of PSHA is often used to identify representative ground motions in terms of magnitude and distance for structural analysis. Forward directivity–induced velocity pulses, which may occur in near-fault (or near-source) motions, are known to cause relatively severe elastic and inelastic response in structures of certain periods. Here, the principles of PSHA are extended to incorporate the possible occurrence of a velocity pulse in a near-fault ground motion. For each magnitude and site-source geometry, the probability of occurrence of a pulse is considered along with the probability distribution of the pulse period given that a pulse does occur. A near-source “narrowband” attenuation law modification to predict ground motion spectral acceleration ( Sa) amplitude that takes advantage of this additional pulse period information is utilized. Further, disaggregation results provide the probability that a given level of ground motion intensity is caused by a pulse-like ground motion, as well as the conditional probability distribution of the pulse period associated with that ground motion. These extensions improve the accuracy of PSHA for sites located near faults, as well as provide a rational basis for selecting appropriate near-fault ground motions to be used in the dynamic analyses of a structure.


2021 ◽  
Author(s):  
Otilio Rojas ◽  
Juan Esteban Rodriguez ◽  
Josep de la Puente ◽  
Scott Callaghan ◽  
Claudia Abril ◽  
...  

<p>Traditional Probabilistic Seismic Hazard Analysis (PSHA) estimates the level of earthquake ground shaking that is expected to be exceeded with a given recurrence time on the basis of  historical earthquake catalogues and empirical and time-independent Ground Motion Prediction Equations (GMPEs). The smooth nature of GMPEs usually disregards some well known drivers of ground motion characteristics associated with fault rupture processes, in particular in the near-fault region, complex source-site propagation of seismic waves, and sedimentary basin response. Modern physics-based earthquake simulations can consider all these effects, but require a large set of input parameters for which constraints may often be scarce. However, with the aid of high-performance computing (HPC) infrastructures  the parameter space may be sampled in an efficient and scalable manner allowing for a large suite of site-specific ground motion simulations that approach the center, body and range of expected ground motions. </p><p>CyberShake is a HPC platform designed to undertake physics-based PSHA from a large suite of earthquake simulations. These simulations are based on seismic reciprocity, rendering PSHA computationally tractable for hundreds of thousands potential earthquakes. For each site of interest, multiple kinematic rupture scenarios, derived by varying slip distributions and hypocenter location across the pre-defined fault system, are generated from an input Earthquake Forecast Model (EFM). Each event is simulated to determine ground motion intensities, which are synthesized into hazard results. CyberShake has been developed by the Southern California Earthquake Center, and used so far to assess seismic hazard in California. This work focuses on the CyberShake migration to the seismic region of South Iceland (63.5°- 64.5°N, 20°-22°W) where the largely sinistral East-West transform motion across the tectonic margin is taken up by a complex array of near-vertical and parallel North-South oriented dextral transform faults in the South Iceland Seismic Zone (SISZ) and the Reykjanes Peninsula Oblique Rift (RPOR). Here, we describe the main steps of migrating CyberShake to the SISZ and RPOR, starting by setting up a relational input database describing potential causative faults and rupture characteristics, and key sites of interest. To simulate our EFM, we use the open source code SHERIFS, a logic-tree method that converts the slip rates of complex fault systems to the corresponding annual seismicity rate. The fault slip rates are taken from a new 3D physics-based fault model for the SISZ-RPOR transform fault system. To validate model and simulation parameters, two validation steps using key CyberShake modeling tools have been carried out. First, we perform simulations of historical earthquakes and compare the synthetics with recorded ground motions and results from other forward simulations. Second, we adjust the rupture kinematics to make slip distributions more representative of SISZ-type earthquakes by comparing with static slip distributions of past significant earthquakes. Finally, we run CyberShake and compare key parameters of the synthetic ground motions with new GMPEs available for the study region. The successful migration and use of CyberShake in South Iceland is the first step of a full-scale physics-based PSHA in the region, and showcases the implementation of CyberShake in new regions.</p>


2006 ◽  
Vol 4 (1) ◽  
pp. 41-50
Author(s):  
Borko Bulajic ◽  
Miodrag Manic

This paper presents a discussion regarding the most common approaches to the deterministic seismic hazard analysis, as well as their relation with the probabilistic hazard analysis. Different methodologies for estimation of the strong earthquake ground motion at a site of interest on the territory of the Republic of Serbia are also discussed. When generation of the synthetic ground motion time histories on the territory of the Republic of Serbia is concerned, a method developed by Trifunac and his associates is suggested having in mind that this approach uses only those input parameters that can be easily and accurately defined while at the same time being able to model all properties of strong earthquake ground motion that are presently known as well as to consider the probabilistic nature of earthquake occurrence.


2021 ◽  
pp. 875529302098197
Author(s):  
Jack W Baker ◽  
Sanaz Rezaeian ◽  
Christine A Goulet ◽  
Nicolas Luco ◽  
Ganyu Teng

This manuscript describes a subset of CyberShake numerically simulated ground motions that were selected and vetted for use in engineering response-history analyses. Ground motions were selected that have seismological properties and response spectra representative of conditions in the Los Angeles area, based on disaggregation of seismic hazard. Ground motions were selected from millions of available time series and were reviewed to confirm their suitability for response-history analysis. The processes used to select the time series, the characteristics of the resulting data, and the provided documentation are described in this article. The resulting data and documentation are available electronically.


2020 ◽  
Vol 18 (14) ◽  
pp. 6119-6148
Author(s):  
Graeme Weatherill ◽  
Fabrice Cotton

Abstract Regions of low seismicity present a particular challenge for probabilistic seismic hazard analysis when identifying suitable ground motion models (GMMs) and quantifying their epistemic uncertainty. The 2020 European Seismic Hazard Model adopts a scaled backbone approach to characterise this uncertainty for shallow seismicity in Europe, incorporating region-to-region source and attenuation variability based on European strong motion data. This approach, however, may not be suited to stable cratonic region of northeastern Europe (encompassing Finland, Sweden and the Baltic countries), where exploration of various global geophysical datasets reveals that its crustal properties are distinctly different from the rest of Europe, and are instead more closely represented by those of the Central and Eastern United States. Building upon the suite of models developed by the recent NGA East project, we construct a new scaled backbone ground motion model and calibrate its corresponding epistemic uncertainties. The resulting logic tree is shown to provide comparable hazard outcomes to the epistemic uncertainty modelling strategy adopted for the Eastern United States, despite the different approaches taken. Comparison with previous GMM selections for northeastern Europe, however, highlights key differences in short period accelerations resulting from new assumptions regarding the characteristics of the reference rock and its influence on site amplification.


Sign in / Sign up

Export Citation Format

Share Document