Site-specific Ground Motion Research Articles

Seismic Vibration Control of Multiple-Degrees-of-Freedom Steel Structures Through Optimal Impact Mass Dampers

ABSTRACT A novel device based on the energy transfer concept, named impact mass damper (IMD), is presented for the mitigation of large structural vibrations due to horizontal seismic excitations. To achieve the optimal performance of the nonlinear device, coupled with different standard steel buildings, an optimization procedure is set both employing the design of experiments (DOEs) method and the Kriging surrogate modelling. Concerning the seismic input, site-specific ground motions and record-to-record variability are taken into account, in agreement with the new European standards. The validated performance of the IMD device demonstrates its benefits, in particular for short-period steel buildings.

Amplitude and duration hazard-consistent ground-motion selection for seismic risk assessment in Mexico City

The emphasis of seismic design regulations on applying nonlinear dynamic analyses (NDAs) promotes using accelerograms that characterize site-specific ground motions. Commonly, amplitude levels of such accelerograms are defined by a target spectrum that could be based on a uniform hazard spectrum (UHS), which is determined by a probabilistic seismic hazard analysis (PSHA) and represents a response spectrum with ordinates having an equal probability of being exceeded within a given return period, Tr\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${T}_{r}$$\\end{document}. Conversely, the definition of ground-motion duration levels is not yet properly defined in current regulations to select accelerograms. Thus, adhering to data handling as that for amplitude ground-motion parameters, this study motivates executing PSHAs to define hazard-consistent levels for the ground-motion duration. That is, accelerograms can be selected to match both amplitude and duration ground-motion levels associated with Tr\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${T}_{r}$$\\end{document}. Further, fragility functions conditional on Tr\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${T}_{r}$$\\end{document} that cover typical performance objectives can be developed using sets of hazard-consistent accelerograms to implement, e.g., multiple stripe analyses (MSAs). To demonstrate the importance of choosing fully hazard-consistent accelerograms to perform NDAs, this study includes the displacement- and energy-based seismic-response evaluation of a steel frame building located at different soil-profile sites in Mexico City. Sets of fully hazard-consistent accelerograms and solely amplitude-based hazard-consistent accelerograms were artificially generated per site for values of Tr\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${T}_{r}$$\\end{document} up to 5000 years. Results indicate that the probability of failure can be underestimated if the ground-motion duration is unvaried in MSAs, e.g., structural damage caused by 50-year return-period or higher events can be more noticeable when fully hazard-consistent accelerograms take place.

Seismic Scenario Losses Estimation Based on the Vulnerability of Pre-Code Masonry Buildings in the Metropolitan Area of Lisbon

ABSTRACT The metropolitan area of Lisbon is the region of Portugal with the highest seismic risk, given the coexistence of moderate-to-high hazard with high population density and building stock exposure. The present work addresses the seismic risk assessment of pre-code masonry buildings in this region, accounting various typological classes and conservation states in order to evaluate the consequences in terms of economic losses. The analyses are conducted through a seismic probabilistic approach, considering the site-specific ground motion through a non-stationary stochastic method. The results can serve as a useful guide for decision-making in developing regional seismic risk mitigation strategies.

Temporal variations of the ‘in-situ’ nonlinear behaviour of shallow sediments during the 2016 Kumamoto Earthquake sequence

SUMMARY Strong ground shaking has the potential to generate significant dynamic strains in shallow materials such as soils and sediments, thereby inducing nonlinear site response resulting in changes in near-surface materials. The nonlinear behaviour of these materials can be characterized by an increase in wave attenuation and a decrease in the resonant frequency of the soil; these effects are attributed to increased material damping and decreased seismic wave propagation velocity, respectively. This study investigates the ‘in-situ’ seismic velocity changes and the predominant ground motion frequency evolution during the 2016 Kumamoto earthquake sequence. This sequence includes two foreshocks (Mw 6 and Mw 6.2) followed by a mainshock (Mw 7.2) that occurred 24 hr after the last foreshock. We present the results of the seismic velocity evolution during these earthquakes for seismological records collected by the KiK-net (32 stations) and K-NET (88 stations) networks between 2002 and 2020. We analyse the impulse response and autocorrelation functions to investigate the nonlinear response in near-surface materials. By comparing the results of the impulse response and autocorrelation functions, we observe that a nonlinear response occurs in near-surface materials. We then quantify the velocity reductions that occur before, during, and after the mainshock using both approaches. This allows us to estimate the ‘in-situ’ shear modulus reduction for different site classes based on V$_{S30}$ values (V$_{S30}\lt 360$ m s−1, $360\lt $V$_{S30}\lt 760$ m s−1 and V$_{S30}\gt 760$ m s−1). We also establish the relationships between velocity changes, shear modulus reduction, variations in predominant ground motion frequencies and site characteristics (V$_{S30}$). The results of this analysis can be applied to site-specific ground motion modelling, site response analysis and the incorporation of nonlinear site terms into ground motion models.

Site-Specific Ground-Motion Waveform Generation Using a Conditional Generative Adversarial Network and Generalized Inversion Technique

ABSTRACT Accurate ground-motion simulations are essential for seismic hazard assessments and engineering practices. Herein, we propose a novel method combining conditional generative adversarial networks (cGANs) and the generalized inversion technique (GIT) to generate site-specific and variability-controlled strong-motion seismograms. The cGANs calculate synthetic seismogram without amplitude scales. The GIT is to separate the source, path, and site characteristics from the Fourier amplitude spectrum (FAS) of the observed seismograms. This method is applied to plate boundary earthquakes off the Pacific coast of Tohoku, Japan. It successfully generates a set of strong-motion seismograms at a given magnitude, distance, and observation station. The output waveforms reproduce the P and S waves as well as coda waves. We validate the method through a quantitative comparison with observed seismograms in terms of both time-domain duration and frequency-domain amplitude characteristics, using metrics of peak ground acceleration (PGA), peak ground velocity, FASs, response spectra, and waveform duration time. The validation results show that the variation in the PGA of the observed seismograms and the synthetic seismograms has a standard deviation of 0.643, and the duration of the seismograms has a standard deviation of 0.346, comparable to the standard deviations seen in the previous studies. Our approach offers high accuracy in stochastic finite-source modeling for a period of 1 s or shorter. The two features of the method, site-specificity and variability control, can contribute to further improvements in seismic hazard assessment by incorporating empirical information based on observed seismograms.

Regional and Site-Specific Ground Motion Model for Probabilistic Seismic Hazard Analysis in Taiwan: A Case Study of I-Lan

In this study, a probabilistic seismic hazard analysis (PSHA) using a traditional approach with a simple logic tree that handles a single-site sigma is developed to present a site-specific ground motion hazard. The analysis relies on the accelerometer data collected from the I-Lan Plain, a deep sedimentary basin in Northeastern Taiwan where hundreds of earthquake sequences (with [Formula: see text] = 4.0–7.6) were recorded by the Taiwan Strong Motion Instrument Program (TSMIP) network between 1992 and 2016. The performance of PSHA on soil and hard sites is evaluated from two representative instrumented sites (ILA048 and ILA025), using residual measures available from a ground motion model (GMM) developed by Phung et al. [ 2023 ] (referred to as Ph23), from which the uncertainties in site terms ([Formula: see text] and single-site sigma ([Formula: see text] are estimated. We address key conceptual issues in the current PSHA approach and introduce a new region- and site-specific PSHA approach in which (1) site-to-site variability ([Formula: see text] is estimated as a random variance in a mixed effects GMM regression and (2) the GMM site-specific single-site sigma ([Formula: see text] is replaced with a generic site-corrected aleatory variability ([Formula: see text]). Comparison of the region- and site-specific hazard curves from our method against the traditional method estimates at two well-recorded sites in the I-Lan region shows an approximate 50% difference in prediction ground motion values considering for 2% and 10% probability of exceedance in 50 years.

Site-specific stochastic ground motion model utilizing deterministic physics-informed simulations: A Bayesian approach

Limited availability of recorded ground motions poses a challenge for reliable probabilistic seismic-hazard analysis (PSHA), even in highly seismic regions like the Western United States. Stochastic ground motions are commonly employed to address this challenge. However, the stochastic ground motion models (GMMs) may not consistently generate ground motions compatible with the site hazard due to their calibration using global data, failing to capture site-specific characteristics adequately. In the absence of recorded motions, physics-informed simulations provide a viable alternative but are deterministic with limitations of their own that makes them challenging to support PSHA. This article introduces a Bayesian framework that combines prior knowledge from a stochastic GMM, calibrated with global data, with site-specific data obtained from deterministic physics-informed simulations. The proposed framework utilizes the Rezaeian–Der Kiureghian (2010) model as the stochastic GMM and incorporates site-specific data from the CyberShake 15.12 study. By updating the mean and variance of the predictive relationships, along with the marginal distribution of the model parameters, through Bayesian inference, this framework allows for the simulation of site-specific ground motions consistent with the site characteristics. The statistics of peak ground acceleration distributions, as well as both the median and variability of the elastic response spectra, obtained from the calibrated stochastic GMM, demonstrate consistency with those derived using GMMs based on the Next Generation Attenuation (NGA) database.

A comparison of ground motions predicted through one-dimensional site response analyses and three-dimensional wave propagation simulations at regional scales

One-dimensional (1D) site response analysis (SRA), which considers vertically propagating seismic waves from the bedrock to the surface, has been a common technique among geotechnical engineers to examine site-specific ground shaking. However, observations from past earthquakes and analytical studies indicate that idealizations ingrained in 1D SRA may be too severe to capture the ground truth, such as the omissions of spatial variability of soil properties, surface topography, and basin and directivity effects. Physics-based three-dimensional ground motion simulations (GMSs) can incorporate these factors and yield more reliable predictions. In this study, we utilize ground motions from 57 physics-based broadband (from 0 to 8–12 Hz) GMS for a region of Istanbul. A total of 2912 sites with experimentally measured soil profiles that are distributed over the 30 km-by-12.5 km area are also modeled as soil columns and analyzed through 1D SRA. The ground responses from 1D SRA and three-dimensional (3D) GMS are then compared for all 57 earthquake scenarios. These systematic comparisons are then used for examining model features that are correlated with variations in the ratios of various ground motion intensity measures (IMs) and for developing regression-based formulas that can be used for determining simple factors for the considered region to correctly scale (up or down) the site-specific ground motion intensities obtained from 1D SRA, including peak ground acceleration (PGA), peak ground velocity (PGV), and spectral acceleration ( Sa) values.

Utility of Stochastic Ground Motion Models in Conditional Spectrum-Based Selection Consistent with the Causal Earthquakes

ABSTRACT Conditional Spectrum (CS), characterized by a conditional mean and conditional standard deviation, serves as a target spectrum linking seismic hazard information to ground motion selection for seismic demand analysis. Current Ground Motion Selection and Modification (GMSM) methods aim to align with this target spectrum but face challenges due to the limited availability of recorded ground motions and potential modifications, leading to difficulties in ensuring selected ground motions accurately reflect site characteristics. Addressing this issue, the study presents an algorithm to select hazard-consistent CS-based records utilizing a stochastic ground motion model with a dual objective: (1) alleviate the need for scaling during the selection process, and (2) select records consistent with the target hazard and the contributing causal scenarios. To achieve both objectives, the study generates a database of hazard-targeted ground motions utilizing the kriging surrogate with scenarios sampled from the disaggregation matrix. Then, a postprocessing step explicitly considers causative parameters in the selection process. The algorithm’s potential to select site-specific ground motions is demonstrated using sites in the Western United States, providing insights into the computational cost and accuracy. Additionally, statistical comparisons are conducted to explore circumstances where similar hazard consistency can be achieved without the postprocessing step, reducing the overall computational cost. Finally, the selected ground motions are favorably compared to those selected through traditional conditional spectrum-based approaches based on their spectral shape and ground motion intensity measures.

Mean-to-median spectral acceleration ratios in Senior Seismic Hazard Analysis Committee Level 3 probabilistic seismic hazard analyses: An analysis of reported results

In many countries, seismic characterization of the site selected for a critical structure or industrial facility is required in terms of site-specific seismic ground motion hazard. For this purpose, a probabilistic seismic hazard analysis (PSHA), performed under the Senior Seismic Hazard Analysis Committee (SSHAC) protocol, is an extended practice for nuclear facilities. In the past decade, SSHAC Level 3 studies have been performed for sites in North America, Europe, Japan, Taiwan, and South Africa. When analyzing PSHA results, the mean-to-median spectral acceleration ratios given by the hazard curves can be interpreted as a measure of the degree of epistemic uncertainty associated with the results. In this article, results of 33 SSHAC Level 3 studies have been used to determine mean-to-median spectral acceleration ratios and the statistics of these ratios, as a function of spectral frequency and annual frequency of exceedance (AFE). The purpose was to develop a reference for the range of uncertainty that is typically captured in this kind of studies. It has been found that, for a given AFE, ratios corresponding to different sites are within a relatively small interval, especially for the spectral frequency band between 2.5 and 10 Hz, which is the band normally more relevant for the seismic design of nuclear installations. In this band, for 10−4 yr−1 AFE, a mean/median ratio of 1.40 would envelop practically all investigated sites.

Construction of Site-Specific Ground Motion Model Based on Multi-Output Gaussian Process

Construction of Site-Specific Ground Motion Model Based on Multi-Output Gaussian Process

Empirical Green’s Function Simulations Toward Site-Specific Ground Motion Prediction for Kopili Fault of NER India

Empirical Green’s Function Simulations Toward Site-Specific Ground Motion Prediction for Kopili Fault of NER India

Characteristics and selection of near-fault simulated earthquake ground motions for nonlinear analysis of buildings

Earthquake-induced ground shaking near rupturing faults is highly sensitive to the rupture characteristics, seismic wave propagation patterns and site conditions, and field recordings of near-fault shaking are relatively sparse. These challenges complicate the assessment of the seismic performance of near-fault structures. A common approach to representing near-fault ground motion in engineering analysis is to explicitly consider and select records with strong directivity pulses (pulse records). We use three-dimensional high-resolution physics-based earthquake simulations to test this approach in the context of scenario-based ground motion record selection, and to study the important characteristics of near-fault ground shaking. We highlight the deficiencies associated with classifying near-fault simulated records as “pulse” or “non-pulse,” based on the presence of a single dominating pulse in the velocity time history. We show that this approach is inadequate for characterizing near-fault shaking on soft soils which can be dominated by both forward rupture directivity and basin amplification effects. We conduct ground motion selection experiments for the analysis of near-fault structures with and without explicit classification of the pulse features in the records, and evaluate the bias in the predicted structural demands. We find that the maximum interstory drift demands on building structures imposed by unscaled site-specific simulated ground motion records selected based on relevant spectral shape features are not sensitive to the classification of records as pulse/non-pulse. Therefore, with regard to predicting the maximum interstory drifts in near-fault buildings, we do not find justification for the binary pulse classification of near-fault records.

Characterizing site-specific ground motion at great depth in a low seismicity region: challenges and perspectives for a nuclear waste repository project

Characterizing site-specific ground motion at great depth in a low seismicity region: challenges and perspectives for a nuclear waste repository project

Seismic response of a single-layer reticulated dome: Fault-to-structure simulation

Seismic response of large-span spatial structures is particularly susceptible to soil-structure interaction and spatial variability of ground motions. Yet existing seismic analysis fail to properly account for the foundation soil, and the ground motions are commonly determined based on empirical formulas. Physics-based 3D seismic simulation has recently been developed to provide site-specific ground motions, and the fault-to-structure simulation facilitates the incorporation of seismic wave field and soil-structure interaction into the analysis of engineering structures. This article innovatively implements the fault-to-structure simulation on a single-layer reticulated dome. Structural response of fault-to-structure model and two commonly used models is compared and discussed. Results show that the damping of soil reduces the stress and horizontal acceleration of the reticulated dome, while the soil movement leads to an increase in displacement. Noticeably, the seismic propagation characteristics increase vertical acceleration of the inner rings. Compared with fault-to-structure model, the multi-support excitation model and uniform excitation model reduce the vertical displacement of 29.1% and 50.1%, respectively, and reduce the vertical acceleration of the inner rings of 59.9% and 69.0%, respectively.

EESD special issue: AI and data‐driven methods in earthquake engineering – (Part 1)

EESD special issue: AI and data‐driven methods in earthquake engineering – (Part 1)

Frequency-dependent site amplification functions for key geological units in Iceland from a Bayesian hierarchical model for earthquake strong-motions

A fundamental modeling assumption in engineering seismology is that seismic ground motion is affected by source, propagation path and site effects. Contrary to the first two, comprehensive studies on site effects, the localized amplification effects on seismic ground motion due to the uppermost geological structure, have been limited in Iceland. As a result, site amplification has been modeled in ground motion models in a simplistic and qualitative way. In this study, we address this limitation by using the new Bayesian hybrid GMM proposed by Kowsari et al. (2020) to model the recorded ground motions at regional strong-motion stations in southwest Iceland, specifically the peak ground acceleration and pseudo-acceleration response spectra as a function of oscillator period (i.e., PGA and PSA at T=0.01-3.0 s), earthquake magnitude, distance and hypocentral depth, and at first, we omit a site-term. Then, we incorporate the GMM in a Bayesian Hierarchical Model (BHM) that allows prior distributions of model parameters to be used and we partition the model residuals into various terms associated with the source, propagation path and station terms. Through Bayesian inference of the model parameters the results show that the station term residuals reveal a frequency-dependent site amplification at the recording sites. By inspecting this systematic behaviour, we propose four groups of station terms by their characteristics. We then carefully inspect the station locations in terms of their associated geological unit and observe that the geological units and station term groups are qualitatively correlated. By combining the two, we propose the first frequency-dependent average site amplification functions associated with key geological units in Iceland based on age and formation. Provisionally we refer to them as hard rock, rock, (Holocene) lava rock, and (unspecified) soil site classes, respectively. We then update the GMM by incorporating these amplification functions so that it now enables the prediction of PGA and PSA at these four site classes. As a result, the new model standard deviations are considerably lower, with those associated with the event term being τ = 0.05, station terms φS2S = 0.10, intra-event φ = 0.17 and total σ = 0.18 in common logarithm units for PGA. The new site class-specific GMM now enables the site-specific ground motion for separate geological units which in turn improves site-specific seismic hazard assessment and informed site-specific engineering decision making in Iceland.

Towards seismic performance quantification of viscous damped steel structure: Site-specific hazard analysis and response prediction

With an emphasis on predictable performance, the paramount importance of performance-based earthquake engineering (PBEE) in quantifying earthquake risks and facilitating the better-informed design of the built environment to achieve earthquake resilience has been widely acknowledged. And the uses of damping systems, i.e., structures incorporated with supplemental damping devices, have been recognized as an effective structural development to increase the resilience of structures to earthquakes. This paper presents a uniform hazard spectrum (UHS) based site-specific ground motion selection procedure, to implement the PBEE framework to relate earthquake hazard to structural performance for structures with variable dynamic properties. This ground motion selection procedure accounts for the effect of variable structural dynamic properties on hazard characterization for structures with damping systems. A paradigm building site, on which the site-specific hazard from the probabilistic seismic hazard analysis (PSHA) is consistent with the risk-targeted hazard from the design maps of ASCE/SEI 7-10, is selected for the implementation of the procedure. Three suites of 40 ground motions representing various hazard intensities at the building site are selected for time history dynamic analysis of a prototype structure with nonlinear viscous damping system. Evaluated is the probability distribution of major engineering demand parameters (EDPs) including story drift, residual story drift, floor velocity, and floor acceleration. The results demonstrated that the major EDPs of a building structure with supplemental nonlinear viscous dampers at a paradigm site can be estimated using the lognormal probability distribution, and the proposed UHS-based site-specific ground motion selection procedure is critical to the performance assessment of structures with variable dynamic properties in the PBEE framework.

An enhanced ground motion selection algorithm for seismic safety assessment of structures

This paper proposes a ground motion selection algorithm (GSA) incorporating a multi-objective optimization scheme for selecting site-specific ground motion records for nonlinear time history analysis. The optimization algorithm utilized is a modified NSGA-II algorithm (MNSGA-Ⅱ) with a custom initialization procedure for ground motion selection that assists in finding records that efficiently match the mean and variance of conditional spectrum. The proposed GSA is applied to three case study regions and is compared with other GSAs such as the greedy optimization algorithm and the original NSGA-II algorithm to test its efficiency. The results show that the proposed GSA efficiently returns records optimally fitting the target spectrums' mean and variance. Moreover, the results show that the proposed optimization scheme converges faster and gives a more optimal results compared to other GSAs.

On the development of region and site-specific ground motion prediction model for the region of I-lan, Taiwan

this study presents the development of a region- and -site specific ground motion model (GMM) based on ground motion data obtained from the I-lan alluvial basin in northeastern Taiwan. The development was performed by calibrating the ergodic Phung et al. (2020) (Ph20) GMM developed for Taiwan. The source of strong ground motion data includes events from shallow crustal and subduction zones with magnitudes (MW) from 4.0 to 7.6, records from distances shorter than 200 km, and focal depths up to 150 km. The resulting regional model accounts for peculiarities that are not captured by the ergodic Ph20 GMMs, such as (1) the difference in source type effects, (2) VS30-scaling from the local site effects, and (3) enhancement of long period spectral ordinates due to the deep sedimentary basins. The standard deviations of the local model are much lower than those of the ergodic model, particularly the component of the standard deviation that corresponds to the site-to-site variability (ϕS2S), indicating that better site characterization can significantly reduce the overall uncertainty in ground motion estimation.