Correlation of JMA instrumental seismic intensity with strong motion parameters

The Japan Meteorological Agency (JMA) seismic intensity (I_JMA) has been used as a measure of strong shaking for many years in Japan, and it necessitates to know the correlation between the JMA seismic intensity and other strong motion indices, e.g. Peak Ground Acceleration (PGA), Peak Ground Velocity (PGV), and Spectrum Intensity (SI). In this study, two strong motion data sets were selected; in which, the first set consists of 879 three‐components non‐liquefied records selected from 13 major earthquake events that occurred in Japan, the United States, and Taiwan, and the second set consists of 17 liquefied records selected from 7 major earthquake events that occurred in Japan and the United States. The JMA seismic intensity and other ground motion indices were calculated using the selected data sets. The relationships between the JMA seismic intensity and PGA, PGV, and SI were then derived performing a two‐stage linear regression analysis. Results from the analysis show that the JMA instrumental seismic intensity shows higher correlation with SI than PGA or PGV, and it shows the highest correlation with the parameters such as the combination of PGA and SI or the product of PGA and SI. The obtained relationships are also compared with the ones obtained by other studies, which may be useful for the disaster management agencies in Japan and deployment of new SI‐sensors that monitor both PGA and SI. Copyright © 2002 John Wiley & Sons, Ltd.     

Freestanding block overturning fragilities: Numerical simulation and experimental validation

The overturning fragilities of symmetric and asymmetric freestanding blocks, ranging in height from 0.54 to 3.6 m and with height‐to‐width ratios ranging from 2.1 to 6.6, are determined numerically. A probabilistic formulation regularizes the overturning responses when exposed to earthquake‐like random‐vibration waveforms. The peak amplitude of the forcing excitation (peak ground acceleration or PGA) is parameterized as a function of the block size, block shape, overturning probability, and either the PGA normalized peak ground velocity (PGV/PGA), spectral acceleration at 1 s (S_a(1)/PGA), or spectral acceleration at 2 s (S_a(2)/PGA). These later intensity measures are correlated with the duration of the predominant acceleration pulse. The overturning fragilities are compared with shake table experiments using blocks ranging in height from ～0.2 to 1.2 m and with height‐to‐width ratios ranging from ～2 to 10. Excitations utilized in the shake table experiments include recordings of the 1979 Imperial Valley, 1985 Michoacan, 1999 Duzce, 1999 Chi‐Chi, and 2002 Denali Earthquakes along with synthetic waveforms. The overturning fragilities accurately represent the overturning responses of blocks with simple basal contact conditions. Objects with multiple rocking points, such as precariously balanced rocks, are more fragile than predicted. Nondestructive tilting tests are used to account for blocks with complex basal contact conditions, demonstrating that these blocks overturn similarly to more slender blocks with simple contact conditions. Copyright © 2008 John Wiley & Sons, Ltd.     

Design Earthquakes in the U.K.

Three sites in the UK are taken, representative of low, medium and high hazard levels (by UK standards). For each site, the hazard value at 10⁻⁴ annual probability is computed using a generic seismic source model, and a variety of ground motion parameters: peak ground acceleration (PGA), spectral acceleration at 10 Hz and 1 Hz, and intensity. Disaggregation is used to determine the nature of the earthquakes most likely to generate these hazard values. It is found (as might be expected) that the populations are quite different according to which ground motion parameter is used. When PGA is used, the result is a rather flat magnitude distribution with a tendency to low magnitude events (\le 4.5 ML) which are probably not really hazardous. Hazard-consistent scenario earthquakes computed using intensity are found to be in the range 5.8–5.9 ML, which is more in accord with the type of earthquake that one expects to be a worst-case event in the UK. This revised version was published online in July 2006 with corrections to the Cover Date.     

Determination of Q S, spreading factor and S-wave acceleration by inversion of seismic intensity data

In this paper, on the basis of intensity data from 85 earthquakes occurred in China, the method of resolving overdetermined equations by using the damped least squares method is applied to inversing for the mean value of 1 s. S wave quality factor, the spreading factor and the S-wave acceleration at the foci of earthquakes in 7 regions of China. The relation between S-wave acceleration at the foci of earthquakes and magnitude is discussed. As an example, 2-DQ _S distribution in Taiwan Province is obtained by inversion. It is found that there exists some corresponding relation between this distribution and Bouguer gravitational anomaly and seismicity. Preliminary analysis and discussion on results of the inversion, and estimation of error inQ _S are made. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 202–211, 1991. The English version of this paper is improved by Professor Yushou Xie.     

The space and time distribution characteristics of the shear stress field for the sequence of the Wuding earthquake

Follow Chen and Duda's model of spectral fall-off of (3, the dependence of peak parameters of ground motion, peak displacement dm, peak velocity vm and peak acceleration am, upon the environment stress (0-values are studied using near source seismic digital recordings for the sequence of the Wuding, Yunnan, M = 6.5 earthquake, in which, as a new thought, the peak parameters are assumed to be related to the medium Q-value. Three formulae for estimating the environment stress (0-values by the peak parameters of three types of ground motions are derived. Using these formulae, the environment stress (0-values are calculated for the sequence of the Wuding earthquake. The result show that (0-values calculated by the three formulae are constant largely, the averages of (0 are in the range of 5.0～35 MPa for most earthquakes. It belongs to the high-stress earthquakes sequence: the high-stress values are restricted to the relatively small area closely near to the epicenter of the main shock. The fine distribution structure for the contours of the environment stress (0-values is related closely to the strong aftershocks. The analysis of spatial and temporal feature of (0-values suggests that the earthquakes sequence in a rupture process generated at the specific intersection zone of seismo-tectonics under high-stress background.     

Seismic fragility of traction elevators

Second‐generation performance‐based earthquake engineering (PBEE‐2) requires a library of component fragility functions to estimate probabilistic damage to a wide variety of building components. The present work draws on a large body of (mostly) post‐earthquake reconnaissance and (some) post‐earthquake survey observations of traction elevators to create fragility functions useful in PBEE‐2. Two surveys provide detailed observations of 115 representative elevators at 12 hospitals shaken in the 1989 Loma Prieta and 1994 Northridge earthquakes and selected without regard to or foreknowledge of damage. Of these, 55 failed and 60 did not. Approximately half were installed after an important code change of 1972, so one can distinguish the performance of pre‐1973 and post‐1973 elevator construction. They experienced a range of strong motion: 22 with peak ground acceleration (PGA) < 0.25 g, 93 with 0.25 g < PGA < 0.85 g. The hospitals had elevator failure rates as low as 0% and as high as 100%. A third survey describes damage qualitatively for six sites with PGA ≤ 0.25 and per‐site failure rates of 0% to perhaps 30%. Fragility functions are offered where the damage state is the loss of functionality of the elevator. The elevators in these surveys exhibit a median capacity of PGA ≈ 0.35 g with a logarithmic standard deviation of 0.40. Capacity is modestly sensitive to whether the elevator was installed before or after 1973. Using building‐specific intensity measures such as S_a(T₁) does not improve the fragility functions. Copyright © 2015 John Wiley & Sons, Ltd.     

Collapse assessment of moment frame buildings,considering vertical ground shaking

While many cases of structural damage in past earthquakes have been attributed to strong vertical ground shaking, our understanding of vertical seismic load effects and their influence on collapse mechanisms of buildings is limited. This study quantifies ground motion parameters that are capable of predicting trends in building collapse because of vertical shaking, identifies the types of buildings that are most likely affected by strong vertical ground motions, and investigates the relationship between element level responses and structural collapse under multi‐directional shaking. To do so, two sets of incremental dynamic analyses (IDA) are run on five nonlinear building models of varying height, geometry, and design era. The first IDA is run using the horizontal component alone; the second IDA applies the vertical and horizontal motions simultaneously. When ground motion parameters are considered independently, acceleration‐based measures of the vertical shaking best predict trends in building collapse associated with vertical shaking. When multiple parameters are considered, Housner intensity (SI), computed as a ratio between vertical and horizontal components of a record (SI_V/SI_H), predicts the significance of vertical shaking for collapse. The building with extensive structural cantilevered members is the most influenced by vertical ground shaking, but all frame structures (with either flexural and shear critical columns) are impacted. In addition, the load effect from vertical ground motions is found to be significantly larger than the nominal value used in US building design. Copyright © 2016 John Wiley & Sons, Ltd.     

谱强度SI和地震动峰值加速度PGA参数在高铁列车地震预警中的应用

为了降低没有破坏性的高频小震对高铁地震报警的干扰,引入谱强度SI参数作为报警参数,分析其在高铁地震报警中的适用性。研究发现低频振动对SI的影响比高频振动大,对于高铁线路附近的高频小震,SI参数能有效的排除,但对于破坏性不大的远震大震可能会引发报警,从而产生误报现象,影响高铁列车运行效率。为了降低SI在破坏性小的远震大震中的误报率,本文引入PGA与SI作为报警参数。最后给出了不同时速下基于PGA与SI联合报警的阈值供参考。     

A study of residuals for strong ground motions in Adapazari basin,NW Turkey,by ground motion prediction equations (GMPEs)

The city of Adapazarı — located in the Marmara Region of northwest Turkey — is situated on a deep sedimentary basin and was the city most heavily damaged by the strong ground motion of the 17 August 1999 Kocaeli earthquake (moment magnitude M_w = 7.4). This study determines site amplifications of the attenuation relationships for shallow earthquakes in the Adapazarı basin by using the previous ground motion prediction equations (GMPEs) and the traditional spectral ratio method. The site amplifications are determined empirically by averaging the residuals between the observed and predicted peak ground acceleration (PGA) and spectral acceleration (SA) values for various periods. Residuals are significantly correlated with the known characteristics of geological units. A new attenuation model has also been developed for 5% damped spectral acceleration to determine the dependence of strong ground motions on frequency.     

EFFECT OF MODELLING ASSUMPTIONS AND INPUT MOTION CHARACTERISTICS ON SEISMIC DESIGN PARAMETERS OF RC BRIDGE PIERS

Modelling assumptions, boundary and loading conditions have a significant effect on analytical assessment of ductility supply and demand measures for RC bridges, a structural form which had suffered extensively in recent earthquakes. In recognition of the important role played by analysis in advancing seismic design of bridges, this paper is concerned with assessing the effect of model characteristics and earthquake strong-motion selection on analytical action and deformation seismic design parameters. This is of particular significance when viewed in the light of the large capital investment and problems with the satisfaction of dynamic similitude encountered in physical testing of piers and pier-deck assemblies. The models studied range between simple fixed-base cantilever and inclusion of both soil and deck effects, represented by assemblies of springs in translational and rotational degrees of freedom. Moreover, two sets of earthquake records are used in dynamic analysis, each comprising six records covering low, intermediate and high a/v, where a and v are the peak ground acceleration and velocity, respectively. The two sets differ in the scaling procedure employed to bring them to a common level of severity; the first set is obtained by direct acceleration scaling whilst the second utilizes the concept of velocity spectral intensity. The results from static and dynamic analysis, using advanced material characterization and solution procedures, are assessed and discussed. Subject to the limitations of the study, outlined in the paper, the results indicate that the inclusion of deck stiffness and/or soil representation is essential to avail of accurate seismic response parameters. However, the effect of variations in soil stiffness and/or deck torsional rigidity applied in the analysis is rather small, compared to the inclusion/exclusion of the model feature. Moreover, it is also observed that using acceleration scaling leads to much larger scatter in the results than when velocity spectral intensity scaling is used. Finally, the results from two particular earthquakes, Friuli and El Centro, highlight the peril of using a small number of records selected without due consideration to the relationship between their wave form, predominant periods and spectral shapes on the one hand and the response periods of the structure on the other.     

Development of damage probability matrices based on Greek earthquake damage data

A comprehensive study is presented for empirical seismic vulnerability assessment of typical structural types, representative of the building stock of Southern Europe, based on a large set of damage statistics. The observational database was obtained from post-earthquake surveys carried out in the area struck by the September 7, 1999 Athens earthquake. After analysis of the collected observational data, a unified damage database has been created which comprises 180,945 damaged buildings from/after the near-field area of the earthquake. The damaged buildings are classified in specific structural types, according to the materials, seismic codes and construction techniques in Southern Europe. The seismic demand is described in terms of both the regional macroseismic intensity and the ratio αg/ao, where αg is the maximum peak ground acceleration (PGA) of the earthquake event and ao is the unique value PGA that characterizes each municipality shown on the Greek hazard map. The relative and cumulative frequencies of the different damage states for each structural type and each intensity level are computed in terms of damage ratio. Damage probability matrices (DPMs) and vulnerability curves are obtained for specific structural types. A comparison analysis is fulfilled between the produced and the existing vulnerability models.     

Empirical model for estimating Fourier amplitude spectra of ground acceleration in Taiwan region

A collection of ground‐motion recordings (1070 acceleration records) of moderate (5.1⩽M_L⩽6.5) earthquakes obtained during the execution of the Taiwan Strong Motion Instrumentation Program (TSMIP) since 1991 was used to study source scaling model and attenuation relations for a wide range of earthquake magnitudes and distances and to verify the models developed recently for the Taiwan region. The results of the analysis reveal that the acceleration spectra of the most significant part of the records, starting from S‐wave arrival, can be modelled accurately using the Brune's ω‐squared source model with magnitude‐dependent stress parameter Δσ, that should be determined using the recently proposed regional relationships between magnitude (M_L) and seismic moment (M₀) and between M₀ and Δσ. The anelastic attenuation Q of spectral amplitudes with distance may be described as Q=225 ƒ^1.1 both for deep (depth more than 35 km) and shallow earthquakes. The source scaling and attenuation models allow a satisfactory prediction of the peak ground acceleration for magnitudes 5.1⩽M⩽6.5 and distances up to about 200 km in the Taiwan region, and may be useful for seismic hazard assessment. Copyright © 2000 John Wiley & Sons, Ltd.     

Estimation of Strong Ground Motions for 2001 Bhuj (M w 7.6), India Earthquake

The strong ground motions for the 2001 Bhuj (M _w 7.6) India earthquake have been estimated on hard rock and B/C boundary (NEHRP) levels using a recently modified version of stochastic finite fault modeling based on dynamic corner frequency (Motazedian and Atkinson in Bull Seismol Soc Am 95, 995–1010 2005). Incorporation of dynamic corner frequency removes the limitations of earlier stochastic methods. Simulations were carried out at 13 sites in Gujarat where structural response recorder (SRR) recordings are available. In addition, accelerograms were simulated at the B/C boundary at a large number of points distributed on a grid. The corresponding response spectra have also been estimated. The values of peak ground accelerations and spectral accelerations at three periods (0.4, 0.75 and 1.25 s) are presented in the form of contour maps. The maximum value of peak ground acceleration (PGA) in the center of meizoseismal zone is 550 cm/s². The response spectral acceleration in same zone is 900 cm/s² (T = 0.4 s), 600 cm/s² (T = 0.75 s) and 300 cm/s² (T = 1.25 s). The innermost PGA contour is on the fault plane. A comparison of the PGA values obtained at 13 sites in this study with those obtained in earlier studies on the same sites, but employing different methods, show that the present PGA values are comparable at most of the sites. The rate of decay of PGA values is fast at short distances as compared to that at longer distances. The PGA values obtained here put some constraints on the expected values from a similar earthquake in the region. A synthetic intensity map has been prepared from the estimated values of PGA using an empirical relation. A comparison with the reported intensity map of the earthquake shows the synthetic MMI values, as expected, are lower by 1 unit compared to reported intensity map. The contour map of PGA along with the contour maps of spectral acceleration at various periods permit the assessment of damage potential to various categories of houses and other structures. Such information will be quite important in planning of mitigation and disaster management programs in the region.     

Prediction of spatially distributed seismic demands in specific structures: Structural response to loss estimation

A companion paper has investigated the effects of intensity measure (IM) selection in the prediction of spatially distributed response in a multi‐degree‐of‐freedom structure. This paper extends from structural response prediction to performance assessment metrics such as probability of structural collapse; probability of exceeding a specified level of demand or direct repair cost; and the distribution of direct repair loss for a given level of ground motion. In addition, a method is proposed to account for the effect of varying seismological properties of ground motions on seismic demand that does not require different ground motion records to be used for each intensity level. Results illustrate that the conventional IM, spectral displacement at the first mode, S_de(T₁), produces higher risk estimates than alternative velocity‐based IM's, namely spectrum intensity, SI, and peak ground velocity, PGV, because of its high uncertainty in ground motion prediction and poor efficiency in predicting peak acceleration demands. Copyright © 2009 John Wiley & Sons, Ltd.     

Evidence of nonlinear site response in horizontal-to-vertical spectral ratio from near-field earthquakes

The issue is addressed as to whether the horizontal-to-vertical spectral ratio (HVSR) method is sensitive to the amplitude of ground motion from near-field earthquakes. Twenty-one three-component accelerograms from two closely located similar soil sites in the town of Lefkas are used. The recordings represent 17 earthquakes covering a wide range of magnitudes, epicentral distances and azimuths. Peak horizontal accelerations (PGA) and velocities (PGV) lie in the ranges 20–540 cm/s² and 1.4–55.2 cm/s. For each HVS ratio, the site's fundamental-resonance frequency, f_res, is determined visually. Linear correlation analysis shows that f_res is strongly (negatively) correlated to PGA and PGV (r between −0.7 and −0.8); no correlation is found with resonance amplitude or epicentral distance. We show that the observed correlation is attributable to soil nonlinearity and indicate how weak-motion estimates of f_res can be corrected for use in assessing site response during strong shaking.     

Evidence of nonlinear site response in HVSR from SMART1 (Taiwan) data

Deamplification of strong motion and the increase of the effective period of soil deposits are typical nonlinear effects; we seek them in SMART1-array data by applying the horizontal-to-vertical spectral ratio (HVSR) technique. The recordings, from four soil and one rock stations, represent 23 earthquakes (M_L 4.9–7.0); PGA varies between 20–260 cm/s². For each station, mean HVSR curves are calculated for two PGA ranges: <75 cm/s² and >100 cm/s² (weak and strong motion). At the soil stations, the “weak” (linear) and “strong” (nonlinear) responses are significantly different. Below 1–1.8 Hz, the nonlinear response exceeds the linear one. Above 2 Hz, the nonlinear response drops below the linear one and above 4–6 Hz below unity (deamplification). From 10 to 16 Hz, the two responses converge. One soil site shows significant negative correlation between resonance frequency and ground acceleration. Such behaviour agrees with other empirical studies and theoretical predictions. Our results imply that the HVSR technique is sensitive to ground-motion intensity and can be used to detect and study nonlinear site response.     

Determination ofQ S,spreading factor and S-wave acceleration by inversion of seismic intensity data

In this paper, on the basis of intensity data from 85 earthquakes occurred in China, the method of resolving overdetermined equations by using the damped least squares method is applied to inversing for the mean value of 1 s. S wave quality factor, the spreading factor and the S-wave acceleration at the foci of earthquakes in 7 regions of China. The relation between S-wave acceleration at the foci of earthquakes and magnitude is discussed. As an example, 2-DQ _S distribution in Taiwan Province is obtained by inversion. It is found that there exists some corresponding relation between this distribution and Bouguer gravitational anomaly and seismicity. Preliminary analysis and discussion on results of the inversion, and estimation of error inQ _S are made.     

Site Response Evaluation Using the H/V Ratio at the Yan-Liau Station in Hualien, Taiwan

—?Four days after the December 13, 1990 Hualien, Taiwan earthquake (M _L = 6.9), a temporary array of fifteen triaxial digital accelerographs was deployed in the epicenter area to monitor aftershocks. Approximately 600 earthquakes triggered this array during the three-month deployment period. The Yan-Liau station (S63) alone recorded 162 events. Most of the accelerograms at S63 exhibit resonance. We have estimated site responses at the Yan-Liau station using both the single-station spectral ratio (or H/V ratio) method and the traditional spectral ratio method that compares ratios at a soil site with those at a reference hard rock site. Based on site response analyses of S waves and coda waves of ground motion recordings, both types of waves show that the H/V ratio provide a good estimate at the resonant frequency although the site amplification factor is overestimated. In addition, the study of microtremor is also a good alternate for estimating the site predominant frequency. While the ground acceleration (or PGA) gradually increases, the resonant frequency shifts to lower frequencies.     

Housner谱烈度的一种快速近似算法

Housner谱烈度及修正谱烈度作为基于加速度记录时程直接得到的地震动强度参数,与建筑结构破坏及地震宏观烈度存在较高的相关性,是可靠的地震仪器烈度物理参数指标。然而,相对于地面加速度峰值、地面速度峰值等地震动峰值参数,三分量加速度记录对应的谱烈度计算过程较为复杂,耗时相对较长,影响了利用谱烈度确定地震仪器烈度的时效性。基于对强震动加速度记录的统计分析,本文提出了谱烈度的快速近似算法,仅计算4个方向上的谱烈度值,采用其中3点作圆即可获得水平面内谱烈度迹线的近似最大值,使计算速度提高了45倍,且保持了谱烈度作为地震仪器烈度物理指标的精度。利用在汶川M_S 8.0地震等386次M_S > 3.0地震中获取的2701组强震动加速度记录,经可靠性检验,结果表明所提出的Housner谱烈度快速近似算法的计算误差在±4.5%以内,可以同时满足地震仪器烈度速报的可靠性和时效性需求。     

Unified modeling of flow liquefaction and cyclic mobility

Flow liquefaction and cyclic mobility are two different phenomena that are encountered during strong earthquakes. Flow liquefaction is associated with the contractive behavior of loose granular materials and cyclic mobility is associated with the dilative response of both loose and dense granular materials at low confining stresses. These two types of response pertain to the same material and therefore should be modeled in a consistent manner. Whether a soil in a given state exhibits contractive or dilative behavior is dependent on its dilatancy, d=d_v^p/d_q^p. A form of the dilatancy d=d(η,ψ,C) is proposed, where η=q/p is the stress ratio, and ψ and C denote the internal state variables and the intrinsic properties, respectively. It has been shown that such a state-dependent dilatancy is effective in describing both the contractive and dilative behavior. This allows the soil behavior associated with flow liquefaction and cyclic mobility during earthquakes to be modeled in a unified way.