Hazard characterization for alternative intensity measures using the total probability theorem

Since their inception in the 1980s, simplified procedures for the analysis of liquefaction hazards have typically characterized seismic loading using a combination of peak ground acceleration and earthquake magnitude. However, more recent studies suggest that certain evolutionary intensity measures (IMs) such as Arias intensity or cumulative absolute velocity may be more efficient and sufficient predictors of liquefaction triggering and its consequences. Despite this advantage, widespread hazard characterizations for evolutionary IMs are not yet feasible due to a relatively incomplete representation of the ground motion models (GMMs) needed for probabilistic seismic hazard analysis (PSHA). Without widely available hazard curves for evolutionary IMs, current design codes often rely on spectral targets for ground motion selection and scaling, which are shown in this study to indirectly result in low precision of evolutionary IMs often associated with liquefaction hazards. This study presents a method to calculate hazard curves for arbitrary intensity measures, such as evolutionary IMs for liquefaction hazard analyses, without requiring an existing GMM. The method involves the conversion of a known IM hazard curve into an alternative IM hazard curve using the total probability theorem. The effectiveness of the method is illustrated by comparing hazard curves calculated using the total probability theorem to the results of a PSHA to demonstrate that the proposed method does not result in additional uncertainty under idealized conditions and provides a range of possible hazard values under most practical conditions. The total probability theorem method can be utilized by practitioners and researchers to select ground motion time series that target alternative IMs for liquefaction hazard analyses or other geotechnical applications. This method also allows researchers to investigate the efficiency, sufficiency, and predictability of new, alternative IMs without necessarily requiring GMMs.

The power of the little ones: Computed and observed aftershock hazard in Central Italy

Seismic hazard in Central Italy due to the 2016–2017 seismic sequence is modeled using a standard probabilistic aftershock seismic hazard model. Two key features of the model are the consideration of time-dependent aftershock occurrence, modeled by stacking Omori decay curves associated with the three largest ( Mw > 5.5) events, and the incorporation of geologic information by modeling the locations of expected seismicity along realistic fault surfaces. The computed seismic hazard at Amatrice indicates higher hazard values compared to those computed using a conventional time-independent hazard analysis. We then compare the computed hazard curves against empirical hazard curves constructed for 12 individual recording stations in terms of peak ground acceleration, each with at least 35 (and up to 231) recordings. At eight sites, the observed exceedances fall within one standard deviation of the expected mean, while at the remaining sites, the observed exceedances fall outside this range indicating a poorer match. The soil sites are among the stations with the poorest match, suggesting that site effects may not be accurately modeled with the current approach.

Risk-targeted hazard maps for Spain

Many studies have demonstrated that the design of structures in a region through the uniform hazard principle does not guarantee a uniform collapse risk. Even in regions with similar Peak Ground Accelerations (PGAs) corresponding to the same mean return period, the seismic risk in terms of collapse probability will be significantly different mainly due to the shape of the hazard curves as well as uncertainties in structural capacities. In this paper, risk-targeted hazard mapping is being explored in peninsular Spain using a recently updated seismic hazard map. Since risk targeting involves multiple input parameters such as the model parameters of fragility curves, their variability was considered through their probability distribution as observed in reinforced concrete (RC) moment frame buildings, representing the most common building typology in Spain. The influence of the variation of these parameters on the risk results were investigated, and different assumptions for estimating the model parameters of fragility curves are illustrated. These assumptions were included in a fixed (generic) fragility curve or building-site-specific fragility curves. Different acceptable damage states (i.e., collapse and yielding) were considered concerning Spain’s seismicity level. Finally, the maps for risk-targeted design ground motions and risk coefficients are presented. It is outlined that the employment of risk-targeted analysis leads to the modifications for existing design ground motions due to the different shape of the hazard curves across Spain and considering the uncertainty of structural capacity. Moreover, it is found that using the building- and site-specific fragility curves could result in a more uniform seismic risk across the country.

Probabilistic demand assessment of bridges considering the effect of corrosion and seismic loading over time

This paper proposes an approach to calculate demand hazard curves considering the effect of both corrosion and seismic loadings over time. The corrosion is defined as the reduction of the cross-sectional area in the reinforced bars of concrete, induced by chloride ions. Three corrosion phases are considered: starting time of corrosion, cracking, and evolution time. Seismic loads are characterized as a stochastic Poisson process. Uncertainties related to the randomness of geometric properties, mechanical properties, and seismic loadings are considered. The approach is illustrated in a continuous bridge designed to comply with a drift of 0.002. The structure is located in Acapulco, Guerrero, Mexico. Fragility curves and demand hazard curves are obtained at 0, 45, 57, 75, 100, and 125 years, based on the global drift. The effect of both corrosion and seismic loadings over time increase the annual rate of demand up to 308% between 0 years (without damage) and 125 years after the bridge construction.

Simulation-Based Seismic Hazard Assessment Using Monte-Carlo Earthquake Catalogs: Application to CyberShake

ABSTRACT The use of numerical simulations in probabilistic seismic hazard analysis (PSHA) has achieved a promising level of reliability in recent years. One example is the CyberShake project, which incorporates physics-based 3D ground-motion simulations within seismic hazard calculations. Nonetheless, considerable computational time and resources are required due to the significant processing requirements imposed by source-based models on one hand, and the large number of seismic sources and possible rupture variations on the other. This article proposes to use a less computationally demanding simulation-based PSHA framework for CyberShake. The framework can accurately represent the seismic hazard at a site, by only considering a subset of all the possible earthquake scenarios, based on a Monte-Carlo simulation procedure that generates earthquake catalogs having a specified duration. In this case, ground motions need only be simulated for the scenarios selected in the earthquake catalog, and hazard calculations are limited to this subset of scenarios. To validate the method and evaluate its accuracy in the CyberShake platform, the proposed framework is applied to three sites in southern California, and hazard calculations are performed for earthquake catalogs with different lengths. The resulting hazard curves are then benchmarked against those obtained by considering the entire set of earthquake scenarios and simulations, as done in CyberShake. Both approaches yield similar estimates of the hazard curves for elastic pseudospectral accelerations and inelastic demands, with errors that depend on the length of the Monte-Carlo catalog. With 200,000 yr catalogs, the errors are consistently smaller than 5% at the 2% probability of exceedance in 50 yr hazard level, using only ∼3% of the entire set of simulations. Both approaches also produce similar disaggregation patterns. The results demonstrate the potential of the proposed approach in a simulation-based PSHA platform like CyberShake and as a ground-motion selection tool for seismic demand analyses.

Simulation-Based Seismic Hazard Assessment Using Monte-Carlo Earthquake Catalogs: Application to CyberShake

ABSTRACT The use of numerical simulations in probabilistic seismic hazard analysis (PSHA) has achieved a promising level of reliability in recent years. One example is the CyberShake project, which incorporates physics-based 3D ground-motion simulations within seismic hazard calculations. Nonetheless, considerable computational time and resources are required due to the significant processing requirements imposed by source-based models on one hand, and the large number of seismic sources and possible rupture variations on the other. This article proposes to use a less computationally demanding simulation-based PSHA framework for CyberShake. The framework can accurately represent the seismic hazard at a site, by only considering a subset of all the possible earthquake scenarios, based on a Monte-Carlo simulation procedure that generates earthquake catalogs having a specified duration. In this case, ground motions need only be simulated for the scenarios selected in the earthquake catalog, and hazard calculations are limited to this subset of scenarios. To validate the method and evaluate its accuracy in the CyberShake platform, the proposed framework is applied to three sites in southern California, and hazard calculations are performed for earthquake catalogs with different lengths. The resulting hazard curves are then benchmarked against those obtained by considering the entire set of earthquake scenarios and simulations, as done in CyberShake. Both approaches yield similar estimates of the hazard curves for elastic pseudospectral accelerations and inelastic demands, with errors that depend on the length of the Monte-Carlo catalog. With 200,000 yr catalogs, the errors are consistently smaller than 5% at the 2% probability of exceedance in 50 yr hazard level, using only ∼3% of the entire set of simulations. Both approaches also produce similar disaggregation patterns. The results demonstrate the potential of the proposed approach in a simulation-based PSHA platform like CyberShake and as a ground-motion selection tool for seismic demand analyses.

The Making of the NEAM Tsunami Hazard Model 2018 (NEAMTHM18)

The NEAM Tsunami Hazard Model 2018 (NEAMTHM18) is a probabilistic hazard model for tsunamis generated by earthquakes. It covers the coastlines of the North-eastern Atlantic, the Mediterranean, and connected seas (NEAM). NEAMTHM18 was designed as a three-phase project. The first two phases were dedicated to the model development and hazard calculations, following a formalized decision-making process based on a multiple-expert protocol. The third phase was dedicated to documentation and dissemination. The hazard assessment workflow was structured in Steps and Levels. There are four Steps: Step-1) probabilistic earthquake model; Step-2) tsunami generation and modeling in deep water; Step-3) shoaling and inundation; Step-4) hazard aggregation and uncertainty quantification. Each Step includes a different number of Levels. Level-0 always describes the input data; the other Levels describe the intermediate results needed to proceed from one Step to another. Alternative datasets and models were considered in the implementation. The epistemic hazard uncertainty was quantified through an ensemble modeling technique accounting for alternative models’ weights and yielding a distribution of hazard curves represented by the mean and various percentiles. Hazard curves were calculated at 2,343 Points of Interest (POI) distributed at an average spacing of ∼20 km. Precalculated probability maps for five maximum inundation heights (MIH) and hazard intensity maps for five average return periods (ARP) were produced from hazard curves. In the entire NEAM Region, MIHs of several meters are rare but not impossible. Considering a 2% probability of exceedance in 50 years (ARP≈2,475 years), the POIs with MIH >5 m are fewer than 1% and are all in the Mediterranean on Libya, Egypt, Cyprus, and Greece coasts. In the North-East Atlantic, POIs with MIH >3 m are on the coasts of Mauritania and Gulf of Cadiz. Overall, 30% of the POIs have MIH >1 m. NEAMTHM18 results and documentation are available through the TSUMAPS-NEAM project website (http://www.tsumaps-neam.eu/), featuring an interactive web mapper. Although the NEAMTHM18 cannot substitute in-depth analyses at local scales, it represents the first action to start local and more detailed hazard and risk assessments and contributes to designing evacuation maps for tsunami early warning.

The 2018 update of the US National Seismic Hazard Model: Additional period and site class data

As part of the update of the 2018 National Seismic Hazard Model (NSHM) for the conterminous United States (CONUS), new ground motion and site effect models for the central and eastern United States were incorporated, as well as basin depths from local seismic velocity models in four western US (WUS) urban areas. These additions allow us, for the first time, to calculate probabilistic seismic hazard curves for an expanded set of spectral periods (0.01 to 10 s) and site classes (VS30 = 150 to 1500 m/s) for the CONUS, as well as account for amplification of long-period ground motions in deep sedimentary basins in the Los Angeles, San Francisco Bay, Seattle, and Salt Lake City areas. Two sets of 2018 NSHM hazard data (hazard curves and uniform-hazard ground motions) are available: (1) 0.05°-latitude-by-0.05°-longitude gridded data for the CONUS and (2) higher resolution 0.01°-latitude-by-0.01°-longitude gridded data for the four WUS basins. Both sets of data contain basin effects in the WUS deep sedimentary basins. Uniform-hazard ground motion data are interpolated for 2, 5, and 10% probability of exceedance in 50 years from the hazard curves. The gridded data for the hazard curves and uniform-hazard ground motions, for all periods and site classes, are available for download at the U.S. Geological Survey ScienceBase Catalog ( https://doi.org/10.5066/P9RQMREV ). The design ground motions derived from the hazard curves have been accepted by the Building Seismic Safety Council for adoption in the 2020 National Earthquake Hazard Reduction Program Recommended Seismic Provisions.