Source and high-frequency decay parameters for the Azores region for stochastic finite-fault ground motion simulations

Strong ground motion prediction based on finite-fault simulation requires the identification of the fault (strike, dip, length and width), source kinematics parameters (stress drop, rupture velocity and slip distribution), regional crustal properties (geometrical spreading, anelastic structure, and upper crustal amplification and attenuation parameters) and the determination of amplification effects due to the local site geology. The general purpose of this study is to understand source and attenuation properties in the Azores, by the determination of stress drop, quality factor and kappa, through records obtained by the Portuguese digital seismic and accelerometer network. Source Spectra were obtained, for each record, after correcting observed spectra from geometrical spreading and anelastic attenuation effect: quality factor was estimated based on coda decay in the time domain and the kappa parameter was estimated by fitting the high-frequency decay of the acceleration spectrum with a straight line in a log-linear scale. Mean stress drop value was obtained considering that ω-squared model for the source spectra prevails. Parameters kappa, k, and quality factor, Q, have been estimated to be k = (0.075 ± 0.02)s and \(Q(f) = (76 \pm 11)f^{0.69 \pm 0.09}\), respectively. A mean value of stress drop was estimated to be around 90–130 bars for the earthquakes in analysis. The knowledge of source and path parameters, in association with other ground motion parameters, allows improving ground-motion estimates for the Azores and, consequently, will lead to more accurate seismic hazard assessment for the Azores and better characterization of seismic scenarios.     

Empirical analysis of path effects on prediction equations of pseudo‐velocity response spectra in northern Japan

We study path effects on prediction equations of pseudo‐velocity response spectra (natural period of 0.1–5.0 s) in northern Japan, where heterogeneous attenuation structure exists. The path effects have been examined by comparing the regression analysis results for two different prediction equations. The first equation consists of a single term of anelastic attenuation conventionally. The second equation consists of two terms of anelastic attenuation in consideration of the heterogeneous attenuation structure. In the second equation, we divide a source‐to‐site distance into two distances at the attenuation boundary beneath the volcanic front. The boundary is considered to separate the relatively high Q fore‐arc side mantle wedge (FAMW) from the low Q back‐arc side mantle wedge (BAMW). Strong motion records (hypocentral distances less than 300 km) from interplate and intraslab events with Mw 5.1–7.3 are used. Regression analysis results show that the standard errors are significantly reduced by the second prediction equation at short periods (0.1–0.5 s), whereas the difference in standard errors from both prediction equations is negligible at intermediate and long periods. The Qs values (quality factor for S‐wave) converted from two anelastic attenuation coefficients for the second prediction equation are remarkably similar to the path‐averaged Qs values for the FAMW and BAMW by other studies using spectral inversion method. From these findings, we conclude that the path effects on the prediction equation of pseudo‐velocity response spectra are satisfactorily accomplished by the second prediction equation. Copyright © 2009 John Wiley & Sons, Ltd.     

Simulation of synthetic ground motions for specified earthquake and site characteristics

A method for generating a suite of synthetic ground motion time‐histories for specified earthquake and site characteristics defining a design scenario is presented. The method employs a parameterized stochastic model that is based on a modulated, filtered white‐noise process. The model parameters characterize the evolving intensity, predominant frequency, and bandwidth of the acceleration time‐history, and can be identified by matching the statistics of the model to the statistics of a target‐recorded accelerogram. Sample ‘observations’ of the parameters are obtained by fitting the model to a subset of the NGA database for far‐field strong ground motion records on firm ground. Using this sample, predictive equations are developed for the model parameters in terms of the faulting mechanism, earthquake magnitude, source‐to‐site distance, and the site shear‐wave velocity. For any specified set of these earthquake and site characteristics, sets of the model parameters are generated, which are in turn used in the stochastic model to generate the ensemble of synthetic ground motions. The resulting synthetic acceleration as well as corresponding velocity and displacement time‐histories capture the main features of real earthquake ground motions, including the intensity, duration, spectral content, and peak values. Furthermore, the statistics of their resulting elastic response spectra closely agree with both the median and the variability of response spectra of recorded ground motions, as reflected in the existing prediction equations based on the NGA database. The proposed method can be used in seismic design and analysis in conjunction with or instead of recorded ground motions. Copyright © 2010 John Wiley & Sons, Ltd.     

Non‐structure‐specific intensity measure parameters and characteristic period of near‐fault ground motions

This paper focuses on the effects of long‐period pulse of near‐fault ground motions on the structural damage potential. Two sets of near‐fault ground motion records from Chi‐Chi, Taiwan earthquake and Northridge earthquake with and without distinct pulse are selected as the input, and the correlation analysis between 30 non‐structure‐specific intensity measure parameters and maximum inelastic displacements and energy responses (input energy and hysteretic energy) of bilinear single degree of freedom systems are conducted. Based on the frequency characteristic of near‐fault ground motions with remarkable long‐period components, two intensity indices are proposed, namely, the improved effective peak acceleration (IEPA) and improved effective peak velocity (IEPV). In addition a new characteristic period of these ground motions is defined based on IEPA and IEPV. Numerical results illustrate that the intensity measure parameters related to ground acceleration present the best correlation with the seismic responses for rigid systems; the velocity‐related and displacement‐related parameters are better for medium‐frequency systems and flexible systems, respectively. The correlation curves of near‐fault ground motions with velocity pulse differ from those of ground motions without pulse. Moreover, the improved parameters IEPA and IEPV of near‐fault impulsive ground motions enhance the performance of intensity measure of corresponding conventional parameters, i.e. EPA and EPV. The new characteristic period based on IEPA and IEPV can better reflect the frequency content of near‐fault ground motions. Copyright © 2009 John Wiley & Sons, Ltd.     

An improved energy‐based approach for selecting pulse‐like ground motions

This study proposes an improved energy‐based approach for quantitative classification of velocity‐pulse‐like ground motions. The pulse amplitude is determined, in its value and in time location, by the amplitude of the half‐cycle pulse having the largest seismic energy. After conducting statistical analyses, a newly‐determined threshold level for selecting pulse‐like ground motions is derived; and then what followed is a comparison analysis of three pulse‐detecting schemes, one using the wavelet analysis, the other two using the energy concept. It is believed that other than providing a useful way of classifying pulse‐like ground motions for structural demand analysis, knowledge of this work could also benefit the development of the ground motion prediction equations accounting for pulse effects, and further to aid the probabilistic seismic hazard analysis in a near‐fault environment. Copyright © 2016 John Wiley & Sons, Ltd.     

Response spectral attenuation relationships for Singapore and the Malay Peninsula due to distant Sumatran‐fault earthquakes

As part of the effort to assess the seismic hazards of Singapore and the Malay Peninsula, representative ground motion prediction models have to be established. Seven existing attenuation relationships developed for shallow crustal earthquakes in stable continent and active tectonic regions are examined, and they are found to consistently over‐predict the ground motions of Sumatran‐fault earthquakes recently recorded in Singapore. This may be attributed to the differences in the regional crustal structures and distance ranges considered. Since the number of recorded ground motions in the region is very limited, a new set of attenuation relationships is derived based on synthetic seismograms. The uncertainties in rupture parameters, such as stress drop, focal depth, dip and rake angles, are defined according to the regional geological and tectonic settings as well as the ruptures of previous earthquakes. Ground motions are simulated for earthquakes with Mw ranging from 4.0 to 8.0, within a distance range from 174 to 1379km. Besides magnitude and distance, source‐to‐station azimuth is found to influence the amplitudes of the ground motions simulated. Thus, the azimuth is taken as an independent variable in the derived ground motion attenuation relationships. The Sumatran‐fault segments that have the potential to generate a specified level of response spectral accelerations in Singapore and Kuala Lumpur are identified based on the newly derived ground motion models. Copyright © 2003 John Wiley & Sons, Ltd.     

Ground motion characteristics of the Chi‐Chi earthquake of 21 September 1999

The purpose of this paper is to investigate the ground motion characteristics of the Chi‐Chi earthquake (21 September 1999) as well as the interpretation of structural damage due to this earthquake. Over 300 strong motion records were collected from the strong motion network of Taiwan for this earthquake. A lot of near‐field ground motion data were collected. They provide valuable information on the study of ground motion characteristics of pulse‐like near‐field ground motions as well as fault displacement. This study includes: attenuation of ground motion both in PGA and spectral amplitude, principal direction, elastic and inelastic response analysis of a SDOF system subjected to near‐field ground motion collected from this event. The distribution of spectral acceleration and spectral velocity along the Chelungpu fault is discussed. Based on the mode decomposition method the intrinsic mode function of ground acceleration of this earthquake is examined. A long‐period wave with large amplitude was observed in most of the near‐source ground acceleration. The seismic demand from the recorded near‐field ground motion is also investigated with an evaluation of seismic design criteria of Taiwan Building Code. Copyright © 2000 John Wiley & Sons, Ltd.     

A simple ground‐motion prediction model for cumulative absolute velocity and model validation

Cumulative absolute velocity (CAV) is an important ground motion intensity measure used in seismic hazard analysis. Based on the Next Generation Attenuation strong motion database, a simple ground‐motion prediction equation is proposed for the geometric mean of as‐recorded horizontal components of CAVs using mixed regression analysis. The proposed model employs only four parameters and has a simple functional form. Validation tests are conducted to compare the proposed model with the recently developed Campbell–Bozorgnia (CB10) model using subsets of the strong motion database, as well as several recent earthquakes that are not used in developing the model. It is found that the predictive capability of the proposed model is comparable with the CB10 model, which employs a complex functional form and more parameters. The study also corroborates previous findings that CAV has higher predictability than other intensity measures such as the peak ground acceleration. The high predictability of CAV warrants the use of the proposed simple model as an alternative in seismic hazard analysis. Copyright © 2012 John Wiley & Sons, Ltd.     

Establishing absorbed energy spectra—an attenuation approach

For energy‐based seismic design, energy demand in the form of absorbed energy spectra was established by an attenuation relationship. The absorbed energy is proposed for evaluating the energy demand in an inelastic system because the absorbed energy is directly related to the pseudo‐velocity in the elastic case. Based on a total of 273 ground motion records from 15 significant earthquakes in California, an attenuation relationship of the absorbed energy was established from a two‐stage non‐linear regression analysis. This relationship was established for a given earthquake magnitude, source‐to‐site distance, site class, and ductility factor. A similar expression for the normalized absorbed energy was also developed. This study showed that the absorbed energy for near‐field ground motions can be significantly larger than that predicted by the attenuation relationship for normal ground motions. Copyright © 2000 John Wiley & Sons, Ltd.     

Ground-motion prediction equations for the intermediate depth Vrancea (Romania) earthquakes

We present the regional ground-motion prediction equations for peak ground acceleration (PGA), peak ground velocity (PGV), pseudo-spectral acceleration (PSA), and seismic intensity (MSK scale) for the Vrancea intermediate depth earthquakes (SE-Carpathians) and territory of Romania. The prediction equations were constructed using the stochastic technique on the basis of the regional Fourier amplitude spectrum (FAS) source scaling and attenuation models and the generalised site amplification functions. Values of considered ground motion parameters are given as the functions of earthquake magnitude, depth and epicentral distance. The developed ground-motion models were tested and calibrated using the available data from the large Vrancea earthquakes. We suggest to use the presented equations for the rapid estimation of seismic effect after strong earthquakes (Shakemap generation) and seismic hazard assessment, both deterministic and probabilistic approaches.     

Direct use of peak ground motion parameters for the estimation of inelastic displacement ratio of SDOF systems subjected to repeated far fault ground motions

This study is aimed at developing statistical equations to estimate the inelastic displacement ratio of singledegree- of-freedom systems subjected to far fault repeated earthquakes. In the study, peak ground motion parameters are used to define the scatter of the original data. The ratio of peak ground acceleration to peak ground velocity, and peak ground velocity of the ground motion records and structural parameters such as period of vibration and lateral strength ratio are used in the proposed equations. For the development of the equations, nonlinear time history analyses of single-degree-offreedom systems are conducted. Then, the results are used in a multivariate regression procedure. The equations are verified by comparing the estimated results with the calculated results. The average error and coefficient of variation of the proposed equations are presented. The analyses results revealed that the direct use of peak ground motion parameters for the estimation of inelastic displacement ratio significantly reduced the scatter in the original data and yielded accurate results. From the comparative results it is also observed that results obtained using equations specific to peak ground velocity or peak ground acceleration to peak ground velocity ratio are similar.     

Sliding fragility of block‐type non‐structural components. Part 2: Restrained components

This paper focuses on seismic vulnerability assessment of restrained block‐type non‐structural components under sliding response on the basis of seismic inputs specified by current seismic codes. The general representation of restrained equipment considered in this study consists of a rigid block restrained by four post‐tensioned, symmetrically arranged cables. Two sliding‐related failure modes are considered: restraint breakage and excessive absolute acceleration. Fragility analysis is proposed as an appropriate tool to evaluate these failure modes. Sample fragility curves developed through Monte‐Carlo simulations show that the restraint breakage limit state is sensitive to the parameters of the equation of motion. For instance, fragility estimates obtained without taking into account vertical base accelerations can be significantly unconservative for relatively large values of the coefficient of friction. In contrast, the excessive absolute acceleration limit state exhibits little sensitivity to the parameters of the equation of motion. Peak absolute acceleration response is almost always equal to or greater than the horizontal peak base acceleration. Representative results suggest that reasonable response estimates for blocks located at stories other than the ground in multistorey buildings can in general be obtained by simply scaling the ground acceleration to the peak acceleration at the corresponding storey. Copyright © 2002 John Wiley & Sons, Ltd.     

Ground motion parameters of Shillong plateau: One of the most seismically active zones of northeastern India

Strong ground motion parameters for Shillong plateau of northeastern India are examined. Empirical relations are obtained for main parameters of ground motions as a function of earthquake magnitude, fault type, source depth, velocity characterization of medium and distance. Correlation between ground motion parameters and characteristics of seismogenic zones are established. A new attenuation relation for peak ground acceleration is developed, which predicts higher expected PGA in the region. Parameters of ...     

Semi‐empirical model for site effects on acceleration time histories at soft‐soil sites. Part 2: calibration

A previously developed simplified model of ground motion amplification is applied to the simulation of acceleration time histories at several soft‐soil sites in the Valley of Mexico, on the basis of the corresponding records on firm ground. The main objective is to assess the ability of the model to reproduce characteristics such as effective duration, frequency content and instantaneous intensity. The model is based on the identification of a number of parameters that characterize the complex firm‐ground to soft‐soil transfer function, and on the adjustment of these parameters in order to account for non‐linear soil behavior. Once the adjusted model parameters are introduced, the statistical properties of the simulated and the recorded ground motions agree reasonably well. For the sites and for the seismic events considered in this study, it is concluded that non‐linear soil behavior may have a significant effect on the amplification of ground motion. The non‐linear soil behavior significantly affects the effective ground motion duration for the components with the higher intensities, but it does not have any noticeable influence on the lengthening of the dominant ground period. Copyright © 2004 John Wiley & Sons, Ltd.     

Wavelet‐based non‐stationary seismic rocking response of flexibly supported tanks

The non‐stationary rocking response of liquid storage tanks under seismic base excitations including soil interaction has been developed based on the wavelet domain random vibration theory. The ground motion has been characterized through statistical functionals of wavelet coefficients of the ground acceleration history. The tank–liquid–foundation system is modelled as a multi‐degree‐of‐freedom (MDOF) system with both lateral and rocking motions of vibration of the foundation. The impulsive and convective modes of vibration of the liquid in the tank have been considered. The wavelet domain coupled dynamic equations are formulated and then solved to get the expressions of instantaneous power spectral density function (PSDF) in terms of functionals of input wavelet coefficients. The moments of the instantaneous PSDF are used to obtain the stochastic responses of the tank in the form of coefficients of hydrodynamic pressure, base shear and overturning base moment for the largest expected peak responses. Parametric variations are carried out to study the effects of various governing parameters like height of liquid in the tank, height–radius ratio of the tank, ratio of total liquid mass to mass of foundation, and shear wave velocity in the soil medium, on the responses of the tank. Copyright © 2003 John Wiley & Sons, Ltd.     

Copula‐based joint probability function for PGA and CAV: a case study from Taiwan

This study aims to develop a joint probability function of peak ground acceleration (PGA) and cumulative absolute velocity (CAV) for the strong ground motion data from Taiwan. First, a total of 40,385 earthquake time histories are collected from the Taiwan Strong Motion Instrumentation Program. Then, the copula approach is introduced and applied to model the joint probability distribution of PGA and CAV. Finally, the correlation results using the PGA‐CAV empirical data and the normalized residuals are compared. The results indicate that there exists a strong positive correlation between PGA and CAV. For both the PGA and CAV empirical data and the normalized residuals, the multivariate lognormal distribution composed of two lognormal marginal distributions and the Gaussian copula provides adequate characterization of the PGA‐CAV joint distribution observed in Taiwan. This finding demonstrates the validity of the conventional two‐step approach for developing empirical ground motion prediction equations (GMPEs) of multiple ground motion parameters from the copula viewpoint. Copyright © 2016 John Wiley & Sons, Ltd.     

Stochastic procedure for the simulation of synthetic main shock‐aftershock ground motion sequences

According to the current seismic codes, structures are designed to resist the first damaging earthquake during their service life. However, after a strong main shock, a structure may still face damaging aftershocks. The main shock‐aftershock sequence may result in major damage and eventually the collapse of a structure. Current studies on seismic hazard mainly focus on the modeling and simulation of main shocks. This paper proposes a 3‐step procedure to generate main shock‐aftershock sequences of pairs of horizontal components of a ground motion at a site of interest. The first step generates ground motions for the main shock using either a source‐based or site‐based model. The second step generates sequences of aftershocks' magnitudes, locations, and times of occurrence using either a fault‐based or seismicity‐based model. The third step simulates pairs of ground motion components using a new empirical model proposed in this paper. We develop prediction equations for the controlling parameters of a ground motion model, where the predictors are the site condition and the aftershock characteristics from the second step. The coefficients in the prediction equations and the correlation between the model parameters (of the 2 horizontal components of 1 record and of several records in 1 sequence) are estimated using a database of aftershock accelerograms. A backward stepwise deletion method is used to simplify the initial candidate prediction equations and avoid overfitting the data. The procedure, based on easily identifiable engineering parameters, is a useful tool to incorporate effects of aftershocks into seismic analysis and design.     

Non‐stationary seismic response of tanks with soil interaction by wavelets

A wavelet‐based random vibration theory has been developed for the non‐stationary seismic response of liquid storage tanks including soil interaction. The ground motion process has been characterized via estimates of statistical functionals of wavelet coefficients obtained from a single time history of ground accelerations. The tank–liquid–soil system has been modelled as a two‐degree‐of‐freedom (2‐DOF) system. The wavelet domain equations have been formulated and the wavelet coefficients of the required response state are obtained by solving two linear simultaneous algebraic equations. The explicit expression for the instantaneous power spectral density function (PSDF) in terms of the functionals of the input wavelet coefficients has been obtained. The moments of this PSDF are used to estimate the expected pseudo‐spectral acceleration (PSA) response of the tank. Parametric variations are carried out to study the effects of tank height, foundation natural frequency, shear wave velocity of soil and ratio of the mass of tank (including liquid) to the mass of foundation on the PSA responses of tanks. Copyright © 2001 John Wiley & Sons, Ltd.     

Tall building performance‐based seismic design using SCEC broadband platform site‐specific ground motion simulations

The scarcity of strong ground motion records presents a challenge for making reliable performance assessments of tall buildings whose seismic design is controlled by large‐magnitude and close‐distance earthquakes. This challenge can be addressed using broadband ground‐motion simulation methods to generate records with site‐specific characteristics of large‐magnitude events. In this paper, simulated site‐specific earthquake seismograms, developed through a related project that was organized through the Southern California Earthquake Center (SCEC) Ground Motion Simulation Validation (GMSV) Technical Activity Group, are used for nonlinear response history analyses of two archetype tall buildings for sites in San Francisco, Los Angeles, and San Bernardino. The SCEC GMSV team created the seismograms using the Broadband Platform (BBP) simulations for five site‐specific earthquake scenarios. The two buildings are evaluated using nonlinear dynamic analyses under comparable record suites selected from the simulated BBP catalog and recorded motions from the NGA‐West database. The collapse risks and structural response demands (maximum story drift ratio, peak floor acceleration, and maximum story shear) under the BBP and NGA suites are compared. In general, this study finds that use of the BBP simulations resolves concerns about estimation biases in structural response analysis which are caused by ground motion scaling, unrealistic spectral shapes, and overconservative spectral variations. While there are remaining concerns that strong coherence in some kinematic fault rupture models may lead to an overestimation of velocity pulse effects in the BBP simulations, the simulations are shown to generally yield realistic pulse‐like features of near‐fault ground motion records.