Earthquake-induced landslide displacement attenuation models and application in probabilistic seismic landslide displacement analysis

A landslide displacement (DLL) attenuation model has been developed using spectral intensity and a ratio of critical acceleration coefficient to ground acceleration coefficient. In the development of the model,a New Zealand earthquake record data set with magnitudes ranging from 5.0 to 7.2 within a source distance of 175 km is used. The model can be used to carry out deterministic landslide displacement analysis,and readily extended to carry out probabilistic seismic landslide displacement analysis. DLL attenuation models have also been developed by using earthquake source terms,such as magnitude and source distance,that account for the effects of earthquake faulttype,source type,and site conditions. Sensitivity analyses show that the predicted DLL values from the new models are close to those from the Romeo model that was developed from an Italian earthquake record data set. The proposed models are also applied to an analysis of landslide displacements in the Wenchuan earthquake,and a comparison between the predicted and the observed results shows that the proposed models are reliable,and can be confidently used in mapping landslide potential.     

The Umbria–Marche case: some suggestions for the Italian seismic zonation

Probabilistic seismic hazard maps, in terms of spectral acceleration and uniform hazard response spectra at given sites, considering local soil conditions, represent a much more complete estimate of the seismic hazard than the traditional maps in terms of peak ground acceleration or macroseismic intensity. This is particularly true when the requests of urban planners and engineers have to be met. The present analysis shows how some hazard parameters, such as the effective peak acceleration and the spectrum intensity, can well synthesise the overall information available from traditional probabilistic studies, but also suggests that soil condition is a first-order ingredient for effective seismic hazard mapping at national level. Three Italian towns, damaged by the 1997 Umbria–Marche earthquake sequence, are considered as example to demonstrate that: (1) soil condition dependent uniform hazard spectra well approximate actual spectra recorded during some events of the seismic sequence; (2) for these localities, the design spectrum of the present Italian seismic code does not seem adequate.

These considerations have induced the Italian scientific community to propose an updating of the national seismic zonation on the basis of several hazard parameters, that are described in this paper.     

Ground motion prediction equations for horizontal and vertical components of acceleration in Northern Iran

Recent studies have shown that the vertical component of ground motion can be quite destructive on a variety of structural systems. Development of response spectrum for design of buildings subjected to vertical component of earthquake needs ground motion prediction equations (GMPEs). The existing GMPEs for northern Iranian plateau are proposed for the horizontal component of earthquake, and there is not any specified GMPE for the vertical component of earthquake in this region. Determination of GMPEs is mostly based on regression analyses on earthquake parameters such as magnitude, site class, distance, and spectral amplitudes. In this study, 325 three-component records of 55 earthquakes with magnitude ranging from M _w 4.1 to M _w 7.3 are used for estimation on the regression coefficients. Records with distances less than 300 km are selected for analyses in the database. The regression analyses on earthquake parameters results in determination of GMPEs for peak ground acceleration and spectral acceleration for both horizontal and vertical components of the ground motion. The correlation between the models for vertical and horizontal GMPEs is studied in details. These models are later compared with some other available GMPEs. According to the result of this investigation, the proposed GMPEs are in agreement with the other relationships that were developed based on the local and regional data.     

A Bayesian Approach to Seismic Hazard Estmation: Maximum Values of Magnitudes and Peak Ground Accelerations

Earthquake Researeh in Ch一na461 .METHODLet R be some value measured or estimated as a sequenee ina’‘Past”time interval(一丁,O)(I)万‘月,=(RI,…,R。),R,之R。,R=nlaX l二f匕11(RI,…,R,,) Values(l)eould have an arbitrary Physieal nature.BelowweshalleonsiderEq.(l)asearthquakemagnitudes in a given seismic aetive region or logarithms of seismie Peak ground aeeelerations at习given site.Ro isa而nimum eutoff value;it 15 defined by Possibilities of registration systems or wasehosen as the …     

Spectral parameters of ground motion in different regions: comparison of empirical models

Fourier-amplitude spectrum is one of the most important parameters describing earthquake ground motion, and it is widely used for strong ground motion prediction and seismic hazard estimation. The relationships between Fourier-acceleration spectra, earthquake magnitude and distance were analysed for different seismic regions (the Caucasus and Taiwan island) on the basis of ground motion recordings of small to moderate (3.5≤M_L≤6.5) earthquakes. It has been found that the acceleration spectra of the most significant part of the records, starting from S-wave arrival, can be modelled accurately by the Brune's “ω-squared” point-source model. Parameters of the model are found to be region-dependent. Peak ground accelerations and response spectra for condition of rock sites were calculated using stochastic simulation technique and obtained models of source spectra. The modelled ground-motion parameters are compared with those predicted by recent empirical attenuation relationship for California.     

基岩地震动参数与震级和距离的关系

关于强地震对结构的作用,应当考虑随时间变化的地面运动的全过程。 这个过程的统计特征是由震级、中间介质及地表土层的特性决定的。 因为现在已经有了从基岩运动推算地表运动的可靠方法,所以搞清基岩运动的特征并由此对基岩运动的过程作出估计,显得特别重要。 本文对美国27个基岩台上的强震记录的反应普以及几个最大值与震级和震源距离的关系进行了统计分析,找到它们的关系。从这些关系出发,对指定的震级和震源距离的基岩,可以估计它的运动特征,并人工造成一组运动过程。      

砂土震陷分析中多维地震荷载等效循环周数计算

借助有限元数值方法模拟不同砂土试样在多向和单向地震荷载输入条件下的动三轴试验,选取覆盖大震、中震、近场、远场及不同土质条件的155组多向地震输入时程对不同砂土试样分别进行单向和多向加载,建立考虑多向地震荷载作用的等效循环周数计算方法,并研究震级、震中距和砂土特性对等效循环周数比的影响。研究表明,震级和震中距对等效循环周数比影响不明显,砂土特性对等效循环周数比的影响具有主导作用,砂土相对密度为45%、60%、80%和100%时,其对应的等效循环周数比均值分别为1.58、1.75、1.93和2.08,相对密度越大,比值越大,结合该比值,可以较好地应用文中提出方法进行多维地震荷载等效循环周数计算。     

Acceleration Time Histories for a Scenario Earthquake in Moscow at Sites with Different Soil Conditions

According to general seismic zoning maps, Moscow is in an area with a seismic intensity of 5, in which the maximum seismic effect is expected from remote deep-focal earthquakes in the Vrancea zone (Eastern Carpathians, Romania). In our previous studies, an earthquake with a hypocenter at a depth of 80–150 km in the Vrancea zone, a moment magnitude of M_w = 8.0, and a drop in stress of Δσ = 325 bar was used as a scenario earthquake for Moscow. A series of model acceleration time histories for ground vibrations was calculated for this earthquake for the reference local conditions of the Moskva seismic station (Moscow, Pyzhevskii per. 3). In this paper, these acceleration time histories are used to calculate the acceleration time histories and estimate the ground vibration parameters for an scenario earthquake at other sites on the territory of Moscow for which information on soil conditions is available. Since the epicentral distance is large (~1300 km), it can be assumed that changes in the shape and spectral content of the acceleration time histories at different sites in Moscow are only caused by different local conditions.     

Stochastic earthquake source model for ground motion simulation

In the analysis and design of important structures with relatively long life spans, there is a need to generate strong motion data for possible large events. The source of an earthquake is characterized by the spatial distribution of slip on the fault plane. For future events, this is unknown. In this paper, a stochastic earthquake source model is developed to address this issue. Here, 1D and 2D stochastic models for slip distribution developed by Lavallée et al.(2006) are used. The random field associated with the slip distribution is heavy-tailed stable distribution which can be used for large events. Using 236 past rupture models, the spectral scaling parameter and the four stable or Levy's parameters against empirical relationship for known quantities like magnitude or fault length are developed. The model is validated with data from 411 stations of 1999 Chi-Chi earthquake. The simulated response spectrum showed good agreement to actual data. Further the proposed model is used to generate ground motion for the 1993 Killari Earthquake where strong motion data is not available. The simulated mean peak ground velocity was in turn related to the intensity(MSK) and compared against values in the literature.     

An energy-based envelope function for the stochastic simulation of earthquake accelerograms

An energy-based envelope function is developed for use in the stochastic simulation of earthquake ground motion. The envelope function is directly related to the Arias intensity of the ground motion as well to the manner in which this Arias intensity is built-up over time. It is shown that this build-up, represented by a Husid plot, can be very well modelled using a simple lognormal distribution. The proposed envelope makes use of parameters that are commonly available in seismic design situations, either following a deterministic scenario-type analysis or following a more comprehensive probabilistic seismic hazard analysis (PSHA), either in terms of Arias intensity or the more common spectral acceleration. The shape parameters of the envelope function are estimated following the calculation of the analytic envelopes for a large number of records from PEER Next Generation of Attenuation (NGA) database. The envelope may also be used to predict the distribution of peak ground acceleration values corresponding to an earthquake scenario. The distribution thus obtained is remarkably consistent with those of the recent NGA models.     

硬场地实测记录下几种土层地震反应分析程序可靠性对比

以日本KiK-net强震观测台网中硬场地井下记录为样本,对传统等效线性化方法LSSRLI-1（频域）、SHAKE2000,时域非线性方法DEEPSOIL和频域一致等效线性化方法SOILQUAKE等几种计算程序在硬土场地地震反应分析中的可靠性进行对比检验。检验工况包括KiK-net井下台网中地表峰值加速度不小于0.05g的水平硬场地的总计344台次的加速度记录,涉及5个台站,土层厚度6~50m,地表峰值加速度范围0.050~0.805g。结果表明:认为硬土场地以往计算方法能够体现土层放大的认识是片面的,在中强地震动情况下,现有流行方法计算出的地表响应会偏小,强地震动下会严重偏小,会给工程抗震设计提供偏于危险的输入;硬场地中烈度8度以下（地表PGA在0.19g以下）,SHAKE2000和DEEPSOIL、LSSRLI-1（频域）计算结果与实际记录差距可以接受,但烈度8度以上,计算出地表响应均较实际记录明显偏小,且随地震动强度增加差距急剧增大;硬场地中,无论何种烈度,SOILQUAKE16计算的地表加速度响应与实测相当,可体现出土层放大作用。     

Strong motion characteristics of the M_w 6.6 Lushan earthquake,Sichuan,China——an insight into the spatial difference of a typical thrust fault earthquake

Near-field strong ground motions are useful for engineering seismology studies and seismic design, but dense observation networks of damaging earthquakes are still rare. In this study, based on the strong-motion data from the M w 6.6 Lushan earthquake, the ground motion parameters in different spatial regions are systematically analyzed, and the contributions from different effects, like the hanging-wall effect, directivity effect, and attenuation effect are separated to the extent possible. Different engineering parameters from the observed ground motions are compared with the local design response spectra and a new attenuation relation of Western China. General results indicate that the high frequency ground motion, like the peak ground acceleration, on two sides of the fault plane is sensitive to the hanging-wall effect, whereas the low frequency ground motion, like the long period spectral acceleration, in the rupture propagation direction is affected by the directivity effect. Moreover, although the M w 6.6 Lushan earthquake is not a large magnitude event, the spatial difference of ground motion is still obvious; thus, for a thrust faulting earthquake, in addition to the hanging effect, the directivity effect should also be considered.     

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.     

基于人工神经网络的仪器地震烈度预测模型研究

地震发生后,针对能够快速预测震中附近的烈度分布情况的问题,首先对632次地震触发的台站进行筛选,对2 231个台站触发后20 s内有效的7个地震动参数以及震级和震源距的信息进行提取,并利用人工神经网络对所选数据样本进行训练,建立三种有效的预测模型。研究结果显示模型一所选的输入参数为7个,不利用震源参数,在预测中有着较好的时效性,从第1 s到20 s,预测的平均烈度差值逐渐减小到0.45;模型二所选的输入参数为8个,利用了震源距信息,可以用于烈度级别的预测,预测的平均烈度差值逐渐减小到0.36;模型三所选的输入参数为9个,预测结果较好,可用于震后烈度场的实时预测,平均烈度差值逐渐减小到0.31。利用提出的3种模型对两次地震事件进行烈度预测,预测烈度差值取整后分别有95%和76%以上在1以内,有着较好的结果,可以用于地震预警当中。     

Arma models for earthquake ground motions

This paper outlines the use of discrete, autoregressive/moving-average (ARMA) models for identification and estimation of parameters in models derived from analysis of uniformly digitized earthquake ground motion acceleration data. Such models are of equal generality as compared to continuous-time models and have a number of significant advantages for purposes of digital analysis and simulation. The structure of ARMA models is briefly described, their relation to continuous models noted, and results of their application to a number of recorded accelerograms summarized.     

Investigating the role of source mechanism,surface topography,and attenuation on the observed PGA pattern in May 28, 2004, Mw 6.2 Baladeh earthquake (Iran)

In this paper, we use seismic waveform simulation to investigate the influence of source mechanism complexity, surface topography, and quality factor on the observed peak ground motions in May 28, 2004, moment magnitude (Mw) 6.2 Baladeh earthquake. The observed peak ground acceleration (PGA) pattern in this event, which is the biggest earthquake to hit the Central Alborz Mountains of Iran in modern instrumental era, is irregular in some respects. First, the observed PGA contours are elongated toward north-west and, second, the maximum observed PGA value of 1049 cm/s² on the horizontal component of Hasan Keyf station 50 km away from the epicenter is quite high and irregular for an earthquake of this magnitude, at such long distance. In this study, we employ the spectral element method, implemented in SPECFEM3D software package to simulate the 3D wave propagation from several source models in the area. Our results suggest directivity effect is the main cause of the anomalous observations in this earthquake and could account for the elongation of PGA contours and also the anomalous maximum PGA value observed at Hasan Keyf strong motion station. We show that the surface topography has minor effect on the observed peak ground acceleration and the resulting PGA maps. Also by finding the bounds of seismic quality factor effect on the peak ground acceleration values, we show that this factor could not account for the elongation of iso-acceleration contours in the north-west direction.     

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.     

利用加速度或速度记录快速测定矩震级的一种新方法

根据峰值加速度、速度、位移对构造环境剪应力场和地震矩的依赖关系,导出了直接计算矩震级的公式．用北京怀来小台网（速度记录）、唐山强震观测台网（加速度记录）和北京地区台网（速度记录）的地震记录资料,检验了这些公式的可用性．结果表明,利用加速度或速度记录测定的矩震级与常规方法（根据位移谱高度）测定的矩震级十分接近．证明了这些公式的正确性,从而拓宽了用于测定震级的资料类型,即不仅限于位移记录,速度和加速度记录都可用于测定震级．     

相对局部区域震源参数随机不确定性经验关系研究

在进行未来破坏性地震的强地面运动数值模拟时,震源参数选取的准确性对地震动预测的结果影响很大。震源参数的确定存在很多不确定性因素,既包含随机的不确定性因素,又包含认知的不确定性因素。本文在大量地震事件及文献调研的基础上,运用统计学方法对具备随机不确定性特征的震源参数进行统计研究,以震源参数经验公式的形态建立解释其随机性和不确定性的数学模型。为了研究局部地区震源参数的定标关系特征,获得更加适用于局部地震密集区域,尤其是包含中国大陆地区在内的局部区域的震源参数的经验关系,本文从GCMT地震目录中选取了1 700多个M_W≥5.5的地震事件,运用统计学方法研究地震密集地区的震源参数经验关系,包括震级、地震矩、破裂面积等,增加了相对较大的局部范围内凹凸体的地震样本数量,从统计学角度计算更加适合局部区域的震源参数的经验关系。统计结果表明:局部区域震例获得的震源参数的经验关系与不限区域震例获得的经验关系存在差异,尤其是涉及到断层破裂面积、凹凸体相关参数时差异较大,局部区域内震例获得的震源参数的经验关系将更具有代表性。应用本文获得的相对局部区域的经验公式计算未来破坏性地震的强地面运动所需的震源参数时,获得的地震动预测结果将更能体现目标区域真实的地震动特征,进而提高地震动预测结果的可靠性。      

Strong motion characteristics of the <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> 6.6 Lushan earthquake,Sichuan, China — an insight into the spatial difference of a typical thrust fault earthquake

Near-field strong ground motions are useful for engineering seismology studies and seismic design, but dense observation networks of damaging earthquakes are still rare. In this study, based on the strong-motion data from the M w 6.6 Lushan earthquake, the ground motion parameters in different spatial regions are systematically analyzed, and the contributions from different effects, like the hanging-wall effect, directivity effect, and attenuation effect are separated to the extent possible. Different engineering parameters from the observed ground motions are compared with the local design response spectra and a new attenuation relation of Western China. General results indicate that the high frequency ground motion, like the peak ground acceleration, on two sides of the fault plane is sensitive to the hanging-wall effect, whereas the low frequency ground motion, like the long period spectral acceleration, in the rupture propagation direction is affected by the directivity effect. Moreover, although the M w 6.6 Lushan earthquake is not a large magnitude event, the spatial difference of ground motion is still obvious; thus, for a thrust faulting earthquake, in addition to the hanging effect, the directivity effect should also be considered.