Capturing epistemic uncertainty in site response

2020 ◽  
pp. 875529302097097
Author(s):  
Adrian Rodriguez-Marek ◽  
Julian J Bommer ◽  
Robert R Youngs ◽  
Maria J Crespo ◽  
Peter J Stafford ◽  
...  
Keyword(s):  
Ground Motion ◽  
Site Response ◽  
Epistemic Uncertainty ◽  
Hazard Analysis ◽  
Discrete Distribution ◽  
Seismic Source ◽  
Site Specific ◽  
Underlying Distribution ◽  
Response Characteristics ◽  
Motion Models

The incorporation of local amplification factors (AFs) determined through site response analyses has become standard practice in site-specific probabilistic seismic hazard analysis (PSHA). Another indispensable feature of the current state of practice in site-specific PSHA is the identification and quantification of all epistemic uncertainties that influence the final hazard estimates. Consequently, logic trees are constructed not only for seismic source characteristics and ground-motion models (GMMs) but also for the site AFs, the latter generally characterized by branches for alternative shear-wave velocity ( VS) profiles. However, in the same way that branch weights on alternative GMMs can give rise to unintentionally narrow distributions of predicted ground-motion amplitudes, the distribution of AFs obtained from a small number of weighted VS profiles will often be quite narrow at some oscillator frequencies. We propose an alternative approach to capturing epistemic uncertainty in site response in order to avoid such unintentionally constricted distributions of AFs using more complete logic trees for site response analyses. Nodes are included for all the factors that influence the calculated AFs, which may include shallow VS profiles, deeper VS profiles, depth of impedance contrasts, low-strain soil damping, and choice of modulus reduction and damping curves. Site response analyses are then executed for all branch combinations to generate a large number of frequency-dependent AFs. Finally, these are re-sampled as a discrete distribution with enough branches to capture the underlying distribution of AFs. While this approach improves the representation of epistemic uncertainty in the dynamic site response characteristics, modeling uncertainty in the AFs is not automatically captured in this way, for which reason it is also proposed that a minimum level of epistemic uncertainty should be imposed on the final distribution.

scholarly journals Seismic hazard maps of Peshawar District for various return periods

2020 ◽  
Vol 20 (6) ◽  
pp. 1639-1661
Author(s):  
Khalid Mahmood ◽  
Naveed Ahmad ◽  
Usman Khan ◽  
Qaiser Iqbal
Keyword(s):  
Seismic Hazard ◽  
Ground Motion ◽  
Target Location ◽  
Epistemic Uncertainty ◽  
Hazard Analysis ◽  
Source Parameters ◽  
Seismic Source ◽  
Return Periods ◽  
Hazard Maps ◽  
Seismic Hazard Maps

Abstract. Probabilistic seismic hazard analysis of Peshawar District has been performed for a grid size of 0.01∘. The seismic sources for the target location are defined as the area polygon with uniform seismicity. The earthquake catalogue was developed based on the earthquake data obtained from different worldwide seismological networks and historical records. The earthquake events obtained at different magnitude scales were converted into moment magnitude using indigenous catalogue-specific regression relationships. The homogenized catalogue was subdivided into shallow crustal and deep-subduction-zone earthquake events. The seismic source parameters were obtained using the bounded Gutenberg–Richter recurrence law. Seismic hazard maps were prepared for peak horizontal acceleration at bedrock level using different ground motion attenuation relationships. The study revealed the selection of an appropriate ground motion prediction equation is crucial for defining the seismic hazard of Peshawar District. The inclusion of deep subduction earthquakes does not add significantly to the seismic hazard for design base ground motions. The seismic hazard map developed for shallow crustal earthquakes, including also the epistemic uncertainty, was in close agreement with the map given in the Building Code of Pakistan Seismic Provisions (2007) for a return period of 475 years on bedrock. The seismic hazard maps for other return periods i.e., 50, 100, 250, 475 and 2500 years, are also presented.

Residual analysis of large strong-motion flatfiles as a tool for detecting data error and anomalies

2021 ◽  
Author(s):  
Claudia Mascandola ◽  
Giovanni Lanzano ◽  
Francesca Pacor
Keyword(s):  
Ground Motion ◽  
Site Response ◽  
Epistemic Uncertainty ◽  
Strong Motion ◽  
Consistency Check ◽  
Mixed Effect ◽  
Motion Models ◽  
Tree Classifier ◽  
The Common ◽  
Ground Motion Models

<p>The rapid increase of seismic waveforms, due to the increment of seismic stations and continuous real-time streaming to data centres, leads to the need for automatic procedures aimed at supporting data processing and data quality control. In this study, we propose a semi-automatic procedure for the consistency check of large strong-motion datasets, classifying the anomalies observed on the residuals analysis and identifying the possible causes.</p><p>The data collected in the strong-motion databases are usually arranged as parametric tables (called flatfiles), used to disseminate the Intensity Measures (IMs) and the associated metadata of the processed waveforms. This is the current practice for the ITalian ACcelerometric Archive (ITACA, D’Amico et al., 2020) and Engineering Strong Motion (ESM; Lanzano et al. 2019a) databases. The adopted criteria for flatfile compilation are designed to collect IMs and related metadata in a uniform, updated, and traceable way, with the aim of providing datasets useful to develop Ground Motion Models (GMMs) for Probabilistic Seismic Hazard Assessment (PSHA) and engineering applications. Therefore, the consistency check of the flatfiles is a crucial task to improve the quality of the products provided by the waveform services.</p><p>The proposed procedure is based on the residual distributions obtained from ad-hoc ground motion prediction equations for the ordinates of the 5% damped acceleration response spectra. In this study, we focus on the active shallow crust events in ITACA, considering the ITA18 ground motion model (Lanzano et al., 2019b) as a reference for Italy. The total residuals, computed as logarithm difference between observations and predictions, are decomposed in between-event, between-station and event-and-station corrected residuals by applying a mixed-effect regression (Bates et al., 2015). This is the common practice for the (partial) removal of the ergodic assumption in empirical GMMs (e.g., Stafford 2014), where the contribution of the systematic corrective effects of event and station on aleatory variability are identified and shifted to the epistemic uncertainty. Afterward, the proposed procedure is applied to raise a warning in case of anomalous residual values. Warnings are provided when the normalized residuals exceed a certain threshold, in three ranges of periods (i.e., 0.01-0.15 s, 0.15-1 s, 1-5 s). The causes of warnings may be several and may concern the event, the site, the waveform, or a combination of them. Among the possible sources of anomalous trends, the more common are: preliminary or inaccurate event localization or magnitude, wrong soil category assigned based on proxies, misleading tectonic regime assigned to the earthquake, and fault directivity that may cause strong-ground motion amplification in certain directions. Warnings may also raise for peculiarities in the site-response (e.g., large amplifications/de-amplifications at certain frequency-bands) and to the occurrence of near-source effects in the waveforms (see Pacor et al., 2018). Based on the raised warnings, a decision tree classifier is developed to identify the common anomaly sources and to support the consistency check of the semi-automatic procedure.</p><p>This study may help to enhance the waveform services and related products, besides reducing the variability of ground motion models and guiding decisions for site characterization studies and network maintenance.</p>

Benefits of site-specific hazard analyses for seismic design in New Zealand

2015 ◽  
Vol 48 (2) ◽  
pp. 92-99 ◽  
Author(s):  
Brendon A. Bradley
Keyword(s):  
New Zealand ◽  
Seismic Design ◽  
Site Response ◽  
Hazard Analysis ◽  
Seismic Source ◽  
Design Guidelines ◽  
Assessment Process ◽  
Site Specific ◽  
Comparative Examination

This paper summarizes the role site-specific seismic hazard analyses can play in seismic design and assessment in New Zealand. The additional insights and potential improvements in the seismic design and assessment process through a better understanding of the ground motion hazard are examined through a comparative examination with prescriptive design guidelines. Benefits include the utilization of state-of-the-art knowledge, improved representation of site response, reduced conservatism, and the determination of dominant seismic source properties, among others. The paper concludes with a discussion of these relative benefits so that the efficacy of site-specific hazard analysis for a particular project can be better judged by the engineer.

NGA-subduction global ground motion models with regional adjustment factors

2021 ◽  
pp. 875529302110348
Author(s):  
Grace A Parker ◽  
Jonathan P Stewart ◽  
David M Boore ◽  
Gail M Atkinson ◽  
Behzad Hassani
Keyword(s):  
Ground Motion ◽  
Site Response ◽  
Epistemic Uncertainty ◽  
Peak Ground Velocity ◽  
Earthquake Ground Motion ◽  
Motion Models ◽  
Magnitude Scaling ◽  
Aleatory Variability ◽  
Single Station ◽  
Ground Motion Models

We develop semi-empirical ground motion models (GMMs) for peak ground acceleration, peak ground velocity, and 5%-damped pseudo-spectral accelerations for periods from 0.01 to 10 s, for the median orientation-independent horizontal component of subduction earthquake ground motion. The GMMs are applicable to interface and intraslab subduction earthquakes in Japan, Taiwan, Mexico, Central America, South America, Alaska, the Aleutian Islands, and Cascadia. The GMMs are developed using a combination of data inspection, data regression with respect to physics-informed functions, ground-motion simulations, and geometrical constraints for certain model components. The GMMs capture observed differences in source and path effects for interface and intraslab events, conditioned on moment magnitude, rupture distance, and hypocentral depth. Site effect and aleatory variability models are shared between event types. Regionalized GMM components include the model constant (that controls ground motion amplitude), anelastic attenuation, magnitude-scaling break point, linear site response, and sediment depth terms. We develop models for the aleatory between-event variability [Formula: see text], within-event variability [Formula: see text], single-station within-event variability [Formula: see text], and site-to-site variability [Formula: see text]. Ergodic analyses should use the median GMM and aleatory variability computed using the between-event and within-event variability models. An analysis incorporating non-ergodic site response should use the median GMM at the reference shear-wave velocity condition, a site-specific site response model, and aleatory variability computed using the between-event and single-station within-event variability models. Epistemic uncertainty in the median model is represented by standard deviations on the regional model constants, which facilitates scaled-backbone representations of model uncertainty in hazard analyses.

Non-Ergodic Site Response in Seismic Hazard Analysis

2017 ◽  
Vol 33 (4) ◽  
pp. 1385-1414 ◽  
Author(s):  
Jonathan P. Stewart ◽  
Kioumars Afshari ◽  
Christine A. Goulet
Keyword(s):  
Standard Deviation ◽  
Seismic Hazard ◽  
Ground Motion ◽  
Site Response ◽  
Hazard Analysis ◽  
Specific Site ◽  
Uniform Hazard Spectra ◽  
Motion Models ◽  
Site Term ◽  
Site Variability

Probabilistic seismic hazard analyses are usually performed with semi-empirical ground motion models (GMMs) following the ergodic assumption whereby average source, path, and site effects from global databases apply for a specific site of interest. Site-specific site response is likely to differ from the global average conditional on site parameters used in GMMs (typically V S30 and basin depth). Non-ergodic site response can be evaluated using on-site ground motion recordings and/or one-dimensional wave propagation analyses, and allows site-to-site variability to be removed from the within-event standard deviation. Relative to ergodic, non-ergodic hazard analyses often reduce ground motions at long return periods. We describe procedures for replacing the site term in GMMs with a non-ergodic nonlinear mean over its appropriate range of periods (returning to the ergodic mean outside that range). We also present procedures for computing non-ergodic standard deviation by removing site-to-site variability while considering effects of soil nonlinearity. We illustrate application of these procedures, and their effect on hazard curves and uniform hazard spectra, as implemented in OpenSHA.

scholarly journals A ground motion logic tree for seismic hazard analysis in the stable cratonic region of Europe: regionalisation, model selection and development of a scaled backbone approach

2020 ◽  
Vol 18 (14) ◽  
pp. 6119-6148
Author(s):  
Graeme Weatherill ◽  
Fabrice Cotton
Keyword(s):  
United States ◽  
Seismic Hazard ◽  
Ground Motion ◽  
Epistemic Uncertainty ◽  
Hazard Analysis ◽  
Seismic Hazard Analysis ◽  
Eastern United States ◽  
Logic Tree ◽  
Short Period ◽  
Northeastern Europe

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.

A Suite of Alternative Ground-Motion Models (GMMs) for Israel

2021 ◽  
Author(s):  
Soumya Kanti Maiti ◽  
Gony Yagoda-Biran ◽  
Ronnie Kamai
Keyword(s):  
Seismic Hazard ◽  
Ground Motion ◽  
Empirical Data ◽  
Epistemic Uncertainty ◽  
Ground Motions ◽  
Plate Boundary ◽  
Strong Motion ◽  
Engineering Applications ◽  
Motion Models ◽  
Ground Motion Models

ABSTRACT Models for estimating earthquake ground motions are a key component in seismic hazard analysis. In data-rich regions, these models are mostly empirical, relying on the ever-increasing ground-motion databases. However, in areas in which strong-motion data are scarce, other approaches for ground-motion estimates are sought, including, but not limited to, the use of simulations to replace empirical data. In Israel, despite a clear seismic hazard posed by the active plate boundary on its eastern border, the instrumental record is sparse and poor, leading to the use of global models for hazard estimation in the building code and all other engineering applications. In this study, we develop a suite of alternative ground-motion models for Israel, based on an empirical database from Israel as well as on four data-calibrated synthetic databases. Two host models are used to constrain model behavior, such that the epistemic uncertainty is captured and characterized. Despite the lack of empirical data at large magnitudes and short distances, constraints based on the host models or on the physical grounds provided by simulations ensure these models are appropriate for engineering applications. The models presented herein are cast in terms of the Fourier amplitude spectra, which is a linear, physical representation of ground motions. The models are suitable for shallow crustal earthquakes; they include an estimate of the median and the aleatory variability, and are applicable in the magnitude range of 3–8 and distance range of 1–300 km.

Ground-Motion Prediction Model Based on Neural Networks to Extract Site Properties from Observational Records

2021 ◽  
Author(s):  
Tomohisa Okazaki ◽  
Nobuyuki Morikawa ◽  
Asako Iwaki ◽  
Hiroyuki Fujiwara ◽  
Tomoharu Iwata ◽  
...  
Keyword(s):  
Ground Motion ◽  
Ground Acceleration ◽  
Site Specific ◽  
Motion Data ◽  
Motion Models ◽  
Proposed Model ◽  
Single Station ◽  
Input Variables ◽  
Ground Condition ◽  
Ground Motion Models

ABSTRACT Choosing the method for inputting site conditions is critical in reducing the uncertainty of empirical ground-motion models (GMMs). We apply a neural network (NN) to construct a GMM of peak ground acceleration that extracts site properties from ground-motion data instead of referring to ground condition variables given for each site. A key structure of the model is one-hot representations of the site ID, that is, specifying the collection site of each ground-motion record by preparing input variables corresponding to all observation sites. This representation makes the best use of the flexibility of NN to obtain site-specific properties while avoiding overfitting at sites where a small number of strong motions have been recorded. The proposed model exhibits accurate and robust estimations among several compared models in different aspects, including data-poor sites and strong motions from large earthquakes. This model is expected to derive a single-station sigma that evaluates the residual uncertainty under the specification of estimation sites. The proposed NN structure of one-hot representations would serve as a standard ingredient for constructing site-specific GMMs in general regions.

NGA-Subduction research program

2021 ◽  
pp. 875529302110560
Author(s):  
Yousef Bozorgnia ◽  
Norman A Abrahamson ◽  
Sean K Ahdi ◽  
Timothy D Ancheta ◽  
Linda Al Atik ◽  
...  
Keyword(s):  
Ground Motion ◽  
Research Program ◽  
Epistemic Uncertainty ◽  
Specific Model ◽  
Peak Ground Velocity ◽  
Ground Acceleration ◽  
Motion Data ◽  
Motion Models ◽  
Aleatory Variability ◽  
Adjustment Factors

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.

Sign in / Sign up

Export Citation Format

Share Document