Evaluation of models for Fourier amplitude spectra for the Taiwan region

In the Taiwan region, the empirical spectral models for estimating ground-motion parameters were obtained recently on the basis of recordings of small to moderate (5.0≤M_L≤6.5) earthquakes. A large collection of acceleration records from the M_L=7.3 Chi-Chi earthquake (21 September, 1999) makes it possible to test the applicability of the established relationships in the case of larger events. The comparison of ground-motion parameters (Fourier amplitude spectra, peak accelerations and response spectra), which were calculated using the models, and the observed data demonstrates that the models could provide an accurate prediction for the case of the Chi-Chi earthquake and the largest aftershocks. However, there are some peculiarities in the ground-motion frequency content and attenuation that, most probably, are caused by the features of the rupture process of the large shallow earthquake source.     

Predictive ground motion scaling in Switzerland: Best estimates and uncertainties

Spectral ground motion (1 to 15 Hz) as a function of distance is modeled for events spanning 3.0 <M_w ≤ 7.0 in Switzerland. The parameters required to simulate ground motion with a stochastic approach are inverted from 2958 horizontal and vertical component waveforms of small to moderate size events (2.0 ≤ M_{L} ≤ 5.2) in the distance range 10 to 300 km recorded on hard rock sites. Using a Monte Carlo simulation, we establish a significantly different amplification of about a factor of 1.9 between the Alpine Foreland and the Alps. To assess the trade-off between the free parameters of our stochastic model and their influence on the predictive ground motion relationship, we perform a grid search over the five-dimensional solution space. The uncertainties are separated into epistemic and aleatory parts; the main epistemic uncertainty is attributed to the lack of data forM > 5. To constrain the viable models at large magnitudes, results from worldwide scaling studies are evaluated in light of the Swiss data. The model that explains best the low observed stress drops at small magnitudes (Δσ ≅ 3 bar) yet matches observed intensities of historical earthquakes assumes a stress drop increasing with moment asM₀^0.25. For three sites in Switzerland we evaluate the sensitivity of the epistemic uncertainty by computing probabilistic hazard curves. Our model offers the most comprehensive and detailed study of spectral ground motion for Switzerland to date.     

Q and high-frequency strong motion spectra

Empirical scaling equations for Fourier amplitude spectra of strong ground motion are used to describe A₀ and τ in the assumed (high-frequency) shape of strong motion amplitudes: FS(_φ) = A₀e-^πτφ. The res of computed A₀ and τ with other related estimates of spectral amplitudes; (2) smooth decay of strong motion spectral amplitudes up to φ = 25 Hz, without an abrupt low-pass filtering of high frequecies; and (3) good agreement with other estimates of the regionally specific attenuation of high-frequncy seismic waves.As the recorded strong earthquake shaking in the western United States typically samples only the shallow (10 km) and local (100km) characteristics of wave attenuation, the processed strong motion accelerograms can be used as the most direct means of describing the nature of the high-frequency attenuation of the entire strong motion signal for use in earthquake engineering applications. Seismological body wave, Lg and coda wave estimates of Q sample different volumes of the crust surrounding the station, and involve different paths of the waves. These differences must be carefully documented and understood before the results can be used in earthquake engineering characterization of strong motion amplitudes.     

An earthquake scenario for the microzonation of Sofia and the vulnerability of structures designed by use of the Eurocodes

The study of the site effects and the microzonation of a part of the metropolitan Sofia, based on the modelling of seismic ground motion along three cross-sections are performed. Realistic synthetic strong motion waveforms are computed for scenario earthquakes (M=7) applying a hybrid modelling method, based on the modal summation technique and finite differences scheme. The synthesized ground motion time histories are source and site specific. The site amplification is determined in terms of response spectra ratio (RSR). A suite of time histories and quantities of earthquake engineering interest are provided. The results of this study constitute a “database” that describes the ground shaking of the urban area. A case study of experiment-based assessment of vulnerability of a cast-in-situ single storey, industrial, reinforced concrete frame, designed according to Eurocodes 2 and 8 is presented. The main characteristics of damage index and storey drift are discussed for the purposes of microzonation.     

Canyon topography effects on ground motion at Feitsui damsite

The free-field accelerograms along Feitsui Canyon are analyzed and modeled by a numerical scheme to study the effect of canyon topography. Since six strong-motion accelerometers (SC1–SC6) were deployed along the Feitsui Canyon in 1991; there are 14 earthquakes (4.9≤M_L≤6.6) recorded by these stations until June 1996, but only five triggered all six stations. The maximum PGA value is 68.6 cm s⁻² recorded at station SC1. According to the present data, the effect of the dam on the ground motions at canyon stations can be negligible. The amplitude of ground motion on the slopes of the canyon is bigger than that at its trough. The integral equation method is applied to a two dimensional model of Feitsui Canyon to study the effects of the canyon topography. We choose the ground motion of SC3 or SC4 station at the trough of the canyon as the input motion for the model, which is then used to predict the ground motion at the other five stations. Apart from the earthquake close to the damsite, the simple model can reproduce the observed accelerations at all frequencies below 4 Hz. Overall, the numerical method can well predict the ground motion along the canyon, although the high-frequency simulation is underestimated.     

Near‐fault ground motions,and the response of elastic and inelastic single‐degree‐of‐freedom (SDOF) systems

In order to investigate the response of structures to near‐fault seismic excitations, the ground motion input should be properly characterized and parameterized in terms of simple, yet accurate and reliable, mathematical models whose input parameters have a clear physical interpretation and scale, to the extent possible, with earthquake magnitude. Such a mathematical model for the representation of the coherent (long‐period) ground motion components has been proposed by the authors in a previous study and is being exploited in this article for the investigation of the elastic and inelastic response of the single‐degree‐of‐freedom (SDOF) system to near‐fault seismic excitations. A parametric analysis of the dynamic response of the SDOF system as a function of the input parameters of the mathematical model is performed to gain insight regarding the near‐fault ground motion characteristics that significantly affect the elastic and inelastic structural performance. A parameter of the mathematical representation of near‐fault motions, referred to as ‘pulse duration’ (T_P), emerges as a key parameter of the problem under investigation. Specifically, T_P is employed to normalize the elastic and inelastic response spectra of actual near‐fault strong ground motion records. Such normalization makes feasible the specification of design spectra and reduction factors appropriate for near‐fault ground motions. The ‘pulse duration’ (T_P) is related to an important parameter of the rupture process referred to as ‘rise time’ (τ) which is controlled by the dimension of the sub‐events that compose the mainshock. Copyright © 2004 John Wiley & Sons, Ltd.     

Ground Motion Zoning of Santiago de Cuba: An Approach by SH Waves Modelling

— The expected ground motion in Santiago de Cuba basin from earthquakes which occurred in the Oriente fault zone is studied. Synthetic SH-waves seismograms have been calculated along four profiles in the basin by the hybrid approach (modal summation for the path source-profile and finite differences for the profile) for a maximum frequency of 1 Hz. The response spectra ratio (RSR) has been determined in 49 sites, distributed along all considered profiles with a spacing of 900 m. The corresponding RSR versus frequency curves have been classified using a logical-combinatorial algorithm. The results of the classification, in combination with the uppermost geological setting (geotechnical information and geological geometry of the subsoil) are used for the seismic zoning of the city. Three different main zones are identified, and a small sector characterized by major resonance effects, due to the particular structural conditions. Each zone is characterized in terms of its expected ground motion parameters for the most probable strong earthquake (M_S=7), and for the maximum possible (M_S=8).     

2021年青海玛多7.4级地震强震动记录及特征分析

2021年5月22日,青海省果洛州玛多县发生7.4级地震,中国强震动观测网络在主震中捕获16组强震动数据。对48条三分向加速度记录进行基线校正、滤波等常规处理,计算相应的地震动参数,发现位于断层破裂前向位置的63DAW台NS向记录的地震动速度波形具有长周期分量丰富的特征。分析6个典型台站的单自由度加速度反应谱,并与我国建筑抗震设计谱比较,分析此次地震的频谱特性。将实际观测到的PGA、PGV和S_a(T=0.1 s、T=1.0 s、T=2.0 s、T=5.0 s)与国内广泛使用的几种地震动预测模型对比,研究此次地震的影响场。通过分析S_a-S_d曲线,探讨此次地震靠近断层区域地面运动大位移与桥梁落梁震害间可能存在的关系。     

Analysis of the spatial variation of seismic waves and ground movements from smart-1 array data

Specially designed arrays of strong motion seismographs located near earthquake sources are required for engineering studies of the near-source properties and the spatial variation of seismic waves. The SMART-1 array in Taiwan provides good records for this type of study. Careful study of the observed strong motion data permits the identification of wave types, directions and apparent wave velocities. In this paper, a principal direction ratio R (f,α) is defined; this indicates the principal direction of the motion (along a nearly straight line) within the range 0 < R < 1. Vertical motion of the ground is also included in this study. Orbit spectrum analysis is used to verify the identification of wave directions and wave types. The spatial variation of seismic waves along the principal direction is studied. From frequency-domain analysis, mathematical models of the spatial variation of ground displacement are developed using a wave-number spectrum and the cross-spectral density function between two spatial coordinates; these models in turn can provide two alternative models for the random vibration analysis of extensive structures subject to multiple point seismic excitation. The SMART-1 array data gathered during the January 29, 1981 earthquake also are used to demonstrate calculation of the ground strains and differential movements of the array site. From time-domain analysis, the spatial variation of seismic waves is defined for ground motion along the identified principal direction. The time variation of evolutionary spectra characterized by frequency-dependent parameters is used for this formulation. The SMART-1 array data again form the basis for discussion of the spatial variation of model parameters.     

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.     

Intraplate seismicity in the western Bohemian Massif (central Europe): A possible correlation with a paleoplate junction

Locations of the Eger Rift, Cheb Basin, Quaternary volcanoes, crustal earthquake swarms and exhalation centers of CO₂ and ³He of mantle origin correlate with the tectonic fabric of the mantle lithosphere modelled from seismic anisotropy. We suggest that positions of the seismic and volcanic phenomena, as well as of the Cenozoic sedimentary basins, correlate with a “triple junction” of three mantle lithospheres distinguished by different orientations of their tectonic fabric consistent within each unit. The three mantle domains most probably belong to the originally separated microcontinents – the Saxothuringian, Teplá-Barrandian and Moldanubian – assembled during the Variscan orogeny. Cenozoic extension reactivated the junction and locally thinned the crust and mantle lithosphere. The rigid part of the crust, characterized by the presence of earthquake foci, decoupled near the junction from the mantle probably during the Variscan. The boundaries (transitions) of three mantle domains provided open pathways for Quaternary volcanism and the ascent of ³He- and CO₂-rich fluids released from the asthenosphere. The deepest earthquakes, interpreted as an upper limit of the brittle–ductile transition in the crust, are shallower above the junction of the mantle blocks (at about 12 km) than above the more stable Saxothuringian mantle lithosphere (at about 20 km), probably due to a higher heat flow and presence of fluids.     

Simulating time-histories and pseudo-spectral accelerations from the 1992 Cairo earthquake at the proposed El-Fayoum New City Site,Egypt

El-Fayoum New City represents one of the new urban settlements that are recently erected all over Egypt. Because seismic recordings are not available, I used the stochastic method to simulate the largest damaging earthquake from the closest seismic source to the proposed area of the city. To verify the method and its computed results in Egypt, a study termed “method verification” was performed. I found that the October 12, 1992, earthquake (M _b = 5.8) that occurred southwest of Cairo in the vicinity of the Dahshour region, at the coordinates 29.77°N, 31.07°E, is a significant earthquake to the city. The parameters of the path from the hypocenter of the event to the city were taken into consideration. To determine the site parameters, a shallow seismic refraction survey was carried out in the studied area. Accordingly, I simulated time-histories and pseudo-spectral accelerations from the October 12, 1992, earthquake at the location of seismic profiles. Finally, it is demonstrated that the site is characterized by high ground motion amplification factors, producing a high ground motion acceleration value.     

Long period microtremors, microseisms and earthquake damage: Northridge, CA, earthquake of 17 January 1994

Three studies of site amplification factors, based on the recorded aftershocks, and one study based on strong motion data, are compared one with another and with the observed distribution of damage from the Northridge, CA, earthquake of 17 January 1994 (M_L=6.4). In the epicentral area, when the peak ground velocities are larger than v_m≈15 cm/s, nonlinear response of soil begins to distort the amplification factors determined from small amplitude (linear) wave motion. Moving into the area of near-field and strong ground motion (v_m>30 cm/s), the site response becomes progressively more affected by the nonlinear soil response. Based on the published results, it is concluded that site amplification factors determined from small amplitude waves (aftershocks, small earthquakes, coda waves) and their transfer-function representation may be useful for small and distant earthquake motions, where soils and structures respond to earthquake waves in a linear manner. However in San Fernando Valley, during the Northridge earthquake, the observed distribution of damage did not correlate with site amplification determined from spectra of recorded weak motions. Mapping geographical distribution of site amplification using other than very strong motion data, therefore appears to be of little use for seismic hazard analyses.     

基于烈度指标的核电厂地震预警停堆参数研究

核电厂在遭遇超设计基准地震时需要考虑安全停堆,以避免造成核泄漏事故。采用地震动的峰值加速度参数(PGA)作为判别停堆的参数具有一定的局限性,PGA参数不能反映地震动的频谱和持时特征,因而可能引起不必要的停堆。针对此问题,基于我国本土大量实际强震记录,遴选出7种典型的工程相关地震动参数,进而基于地震烈度指标,分析发现地震动的标准累计绝对速度参数能更好地表征地震对核电厂的整体潜在破坏能力,适合作为判别停堆的参数,然后提出两种确定预警参数阈值的方法,最终建议考虑我国强震数据特征的核电厂判别停堆的预警参数阈值,为我国核电厂的地震安全停堆参数的确定提供了参考。     

Peak and integral seismic parameters of L��Aquila 2009 ground motions: observed versus code provision values

The M _w 6.3 L’Aquila earthquake of April 6, 2009 hit a wide area of the Abruzzo region (Central Italy). The epicentre of the main shock was very close to the urban centre of L’Aquila, the regional capital, with an epicentral distance less than 10 km. It was the strongest earthquake ever recorded in Italy which has provided ground motion recordings from accelerometric stations located in close proximity to the epicentre. Because of this, several remarkable results can be achieved by analysing the strong motion recorded signals in terms of peak (PGA, PGV and PGD) and integral (Housner Intensity, I _H) seismic parameters. Additionally, an alternative time-domain representation of recorded signals has been used to furnish a rapid comparison of traces recorded at different stations and along different directions. Some comparisons between the response spectra derived from the recordings and the elastic demand spectra provided in the new seismic Italian code have also been performed. PGA recorded values are very high and generally higher than code values for seismic actions with return period T _R = 475 years. In some cases, this also happens for seismic actions with T _R up to 2,475 years. With regard to I _H, recorded values are generally higher for T _R = 475 years, whilst they are remarkably lower for T _R = 2,475 years. Accurate analyses have been carried out in the article to better understand the above differences and their significance and implications.     

2018年陕西宁强5.3级地震强地面运动特征及局部场地效应分析

2018年9月12日陕西省宁强县发生5.3级地震,中国数字强震动台网的39个专业台站在此次地震中触发。文章中通过处理捕获的117条三分向加速度记录,给出近场台站的地震动参数,绘制震中附近区域峰值加速度等值线图,其长轴呈西南-东北方向展布。采用实际观测数据与几种常用地震动衰减关系对比,发现霍俊荣衰减预测模型能更好地反映此次地震的影响场。将振幅最大的51GYD台的反应谱与我国抗震设计反应谱比较,采用最小二乘法拟合出不同震中距5个台站各周期谱加速度衰减特性,总结出此次地震的反应谱基本特征。运用H/V谱比法对51GYD土层台和62ZM台阵进行局部场地地震反应分析,研究覆盖土层对地震动的放大作用,及局部地形对峰值加速度和峰值速度的影响过程。     

Real-time mapping of earthquake-induced landslides

Rigid sliding block analysis is a common analytical procedure used to predict the potential for earthquake-induced landslides for natural slopes. Currently, predictive models provide the expected level of displacement as a function of the characteristics of the slope (e.g., geometry, strength, yield acceleration) and the characteristics of earthquake shaking (e.g., peak ground acceleration, peak ground velocity). These predictive models are used for developing seismic landslide hazard maps which identify zones with risk of earthquake-induced landslides. Alternatively, these models can be combined with Shakemaps to generate “near-real-time” Slidemaps which could be used, among others, as a tool in disaster management. Shakemaps (a publicly available free service of the United States Geological Survey, USGS) provide near-real-time ground motion conditions during the time of an earthquake event. The ground motion parameters provided by a Shakemap are very useful for the development of Slidemaps. By providing ground motion parameters from an actual earthquake event, Shakemaps also serve as a tool to decouple the uncertainty of the ground motion in sliding displacements prediction. Campania region in Italy is studied for assessing the applicability of using Shakemaps for regional landslide-risk assessment. This region is selected based on the availability of soil shear strength parameters and the proximity to the 1980 Irpina (M _w  = 6.9) Earthquake.     

2016年8月24日意大利中部MW6.2地震强地面运动特征分析

对意大利国家强震台网在2016年8月24日获得的其中部拉齐奥大区阿库莫利市发生的M_W6.2地震强震动三分向记录进行处理和分析。完成原始数据基线校正、滤波等基本数据处理,回归此次地震动幅值衰减规律,发现其整体与ITA08及BA08的衰减趋势一致,但远场实际值低于预测值,不同场地条件下的衰减特性与ITA10一致,近震源幅值较大,且方向性明显;计算并回归分析几种持时,与全球经验预测方程均基本吻合;比较4个幅值较大的近震源台站的反应谱,发现其明显高于欧洲抗震设计规范中的设计反应谱。结合此次震害特点,该地区在实际建设中仍需提高抗震设防能力,以确保安全性等级。     

Seismic vulnerability assessment to slight damage based on experimental modal parameters

The aim of this paper is to adjust behaviour models for each class of structure for vulnerability assessment by using ambient vibration. A simple model based on frequencies, mode shapes and damping, taken from ambient vibrations, allows computation of the response of the structures and comparison of inter‐storey drifts with the limits found in the literature for the slight damage grade, considered here as the limit of elastic behaviour. Two complete methodologies for building fragility curves are proposed: (1) using a multi‐degree of freedom system including higher modes and full seismic ground‐motion and (2) using a single‐degree of freedom model considering the fundamental mode f₀ of the structure and ground‐motion displacement response spectra S_D(f₀). These two methods were applied to the city of Grenoble, where 60 buildings were studied. Fragility curves for slight damage were derived for the various masonry and reinforced concrete classes of buildings. A site‐specific earthquake scenario, taking into account local site conditions, was considered, corresponding to an M_L = 5.5 earthquake at a distance of 15 km. The results show the benefits of using experimental models to reduce variability of the slight damage fragility curve. Moreover, by introducing the experimental modal model of the buildings, it is possible to improve seismic risk assessment at an overall scale (the city) or a local scale (the building) for the first damage grade (slight damage). This level of damage, of great interest for moderate seismic‐prone regions, may contribute to the seismic loss assessment. Copyright © 2011 John Wiley & Sons, Ltd.