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.     

Topographic site effects of Xishan Park ridge in Zigong city,Sichuan considering epicentral distance

Ground motion amplification induced by topography plays a vital role in engineering seismology. A topographic array of 8 accelerographs has been operating along the ridge in Xishan Park since 2007. The topographic site effects in Zigong city are studied based on the strong motion data of 2008 M_s 8.0 Wenchuan earthquake (the epicentral distance?=?225 km) and 2019 M_s 5.2 Zizhong earthquake (the epicentral distance?=?29 km). We compare the peak ground acceleration (PGA) of the two earthquakes and find that the PGA of Station 7#, which locates on a relatively steep slope, is amplified by 4.41 times comparing with the reference station in Zizhong earthquake, while this value is only 1.62 in Wenchuan earthquake. Fourier amplitude spectrum shows that the high frequency content of Zizhong earthquake is more abundant because of its smaller epicentral distance. By using the standard spectral ratio (SSR) method, we conclude that the amplification occurs because high-frequency ground motion is likely to resonate at small-scale features. Finally, the 3D numerical simulations are used to verify these conclusions. Our work indicates that more sophisticated numerical models need to be established for more accurate topographic site effects quantification. In addition, the influence of nearby topographic features should be considered when selecting reference stations.

     

Characterizing and quantifying earthquake-induced site nonlinearity

This study proposes a recording-based approach to characterize and quantify earthquake-induced site nonlinearity, exemplified as soil nonlinearity and/or liquefaction. Alternative to Fourier spectral analysis (FSA), the paper introduces time–frequency analysis of earthquake ground motion recordings with the aid of so-called Hilbert-Huang transform (HHT), and offers justification for the HHT in addressing the nonlinear features shown in the recordings. With the use of the 2001 Nisqually earthquake recordings, this study shows that the proposed approach is effective in characterizing site nonlinearity and quantifying the influences in seismic ground responses. Major results from this study are listed below.

1. HHT-based site amplification is defined as the ratio of marginal Hilbert amplitude spectra, similar to the Fourier-based one that is the ratio of Fourier amplitude spectra. The HHT-based site amplification can be used in effectively quantifying site nonlinearity in terms of frequency downshift in the low-frequency range and amplification-reduction factor in intermediate-frequency range in comparison with the Fourier-based one.

2. Instantaneous damping, and Hilbert and marginal Hilbert damping spectra are defined in ways similar to instantaneous frequency, and Hilbert and marginal Hilbert amplitude spectra, respectively. Consequently, the HHT-based site damping is found as the difference of marginal Hilbert damping spectra, which can be used as an alternative, complementary index to measure the influences of site nonlinearity in seismic ground responses.

The response law of far-field seismic ground motion of the Wenchuan earthquake and its damaging mechanism in the Loess Plateau

A series of housing collapses and other serious damage was caused by the 2008 Wenchuan M_S 8.0 earthquake in the seismic intensity Ⅵ areas of the Loess Plateau, which is hundreds of kilometers away from the epicenter, and which showed a remarkable seismic intensity anomaly. The seismic disasters are closely related to the seismic response characteristics of the site, therefore, the systematic study of the far-field seismic response law of the Wenchuan earthquake in the Loess Plateau is of great significance to prevent the far-field disaster of great earthquake. In this paper, the seismic acceleration records of several bedrock stations and loess stations from the seismogenic fault of the Wenchuan earthquake to the Loess Plateau were collected, and the attenuation law of ground motion along the propagation path and the characteristics of seismic response on the loess site are studied, and the mechanism of amplification effect of ground motion is analyzed based on the dynamic feature parameters of the loess site obtained through the HVSR method. Taking a typical loess site of thick deposit as the prototype, a series of shaking table tests of dynamic response of loess site models with different thicknesses were carried out. Amplification effect, spectral characteristics of acceleration in model sites were analyzed under the action of a far-field seismic wave of the Wenchuan earthquake. The results show that seismic attenuation on the propagation path along the NE strike of the seismogenic fault to the Loess Plateau is slower than that in other directions, and the predominant period range of ground motion on bedrock site of the Loess Plateau presents broadband characteristics. Because the natural periods of loess sites with thick deposits are within the predominant period range of bedrock input wave, loess sites appear significant amplification effect of ground motion, the horizontal acceleration of ground motion exceeds 0.1 ?g, the seismic intensity reaches 7°. The thicker the loess deposit is, the more significant the change of spectral characteristics of ground motion on loess sites, and the narrower the predominant period range of ground motion becomes, and the closer it is to the natural period of loess sites. Therefore, for some old houses on thick loess sites, the poor seismic performance and strong seismic response eventually led to their collapses and damages because their natural periods are very close to the predominant period of ground motion of the Wenchuan earthquake on thick loess sites; For these damaged high-rise buildings, the resonance effect might be the main reason for their damages because their natural periods are included in the predominant period range of ground motion of the Wenchuan earthquake on thick loess sites.These research results would provide a basis for seismic disasters prediction and evaluation and seismic design of construction engineering in the Loess Plateau.     

海域不同类别场地地震动参数变化规律研究

基于海域场地分类标准,选取南海海域实测钻孔作为计算模型,同时人工构造部分钻孔计算模型,对126条不同特性地震动输入下5种类别场地计算模型开展土层反应分析计算,分析不同地震动输入下不同海域场地峰值加速度和特征周期变化规律.结果表明,场地类别和地震动输入强度显著影响峰值加速度放大系数和特征周期,场地土越软,地震动输入强度越...     

黄土斜坡动力响应特征分析

斜坡动力响应特征与斜坡形态密切相关,若入射地震波主频接近斜坡卓越频率就会放大斜坡动力响应,甚至造成斜坡失稳。汶川地震对远离震中的黄土地区造成了较为严重的破坏,局部场地震害和地震动放大效应显著。选取汶川地震典型黄土斜坡场地,利用地形台阵流动观测和数值模拟计算相结合的方法,系统开展强震动作用下黄土斜坡场地动力响应特征研究。结果表明:坡顶卓越频率最小,其PGA放大系数甚至达到坡底的1.98,这种现象可能与斜坡高差和入射波波长之比密切相关,比值0.2时坡顶放大效应达到最大。随斜坡坡度增加,放大效应增强,坡顶反应谱卓越周期放大系数可达5,说明斜坡地形对强震地面运动有显著影响。数值计算结果与实际强震观测基本吻合,其结果对黄土地区建设工程抗震设防具有重要的科学与实际意义。     

Evaluation of ASCE 7 equations for designing acceleration‐sensitive nonstructural components using data from instrumented buildings

This study uses instrumented buildings and models of code‐based designed buildings to validate the results of previous studies that highlighted the need to revise the ASCE 7 F_p equation for designing nonstructural components (NSCs) through utilizing oversimplified linear and nonlinear models. The evaluation of floor response spectra of a large number of instrumented buildings illustrates that, unlike the ASCE 7 approach, the in‐structure and the component amplification factors are a function of the ratio of NSC period to the supporting building modal periods, the ground motion intensity, and the NSC location. It is also shown that the recorded ground motions at the base of instrumented buildings in most cases are significantly lower than design earthquake (DE) ground motions. Because ASCE 7 is meant to provide demands at a DE level, for a more reliable evaluation of the F_p equation, 2 representative archetype buildings are designed based on the ASCE 7‐16 seismic provisions and exposed to various ground motion intensity levels (including those consistent with the ones experienced by instrumented buildings and the DE). Simulation results of the archetype buildings, consistent with previous numerical studies, illustrate the tendency of the ASCE 7 in‐structure amplification factor, [1 + 2(z/h)] , to significantly overestimate demands at all floor levels and the ASCE 7 limit of to in many cases underestimate the calculated NSC amplification factors. Furthermore, the product of these 2 amplification factors (that represents the normalized peak NSC acceleration) in some cases exceeds the ASCE 7 equation by a factor up to 1.50.     

Case study of strong ground motion variations across cut slope

We evaluate the influence of topography on motions recorded at the base and crest of an approximate 3H:1V, 20 m single-faced slope. The motions were recorded during the 1983 Coalinga earthquake mainshock and two aftershocks. Mainshock peak accelerations at the crest and base transverse to the slope face were 0.59 and 0.38 g, respectively. The spectral amplification of crest motion occurred across T≈0–2 s. Differences between the crest/base motions are postulated to result principally from soil-structure interaction (base instrument is in a structure), variations in local ground response, and topography. Transfer functions quantifying soil-structure interaction (SSI) effects are evaluated and the base motion is modified at short periods to correct it to an equivalent free-field motion. The different levels of ground response at the crest and base are identified based on location-specific measurements of soil shear wave velocities. Differences between crest/base motions not accounted for by SSI or differential ground response are attributed to topographic effects. By these means, topographic spectral amplification (i.e. amplification relative to level ground conditions) is estimated to be about 1.2 at the crest and about 0.85–0.9 at the base across the period range T≈0.4–1.0 s.     

Site amplification in lava rock on soft sediments

Seismicity and volcanic activity in Iceland are related to the Mid-Atlantic plate boundary that crosses the island. Due to volcanic activity, different sea levels through the ages and glacial drift, the geology in Iceland is quite complex at many sites. In June 2000, two major earthquakes of magnitude 6.6 (M_w) and 6.5 (M_w) occurred in southern Iceland. Ground motion from these main shocks and a number of aftershocks were recorded at the Icelandic strong motion network operated by the University of Iceland. At one of the instrumented sites considerable amplification was recorded on lava-rock overlying alluvial deposits. This site fits poorly in the soil classification systems of most earthquake codes. In this paper the recorded data is analyzed with different methods and the results are compared with the results of a one-dimensional site response analysis. The different methods produce the same characteristic of soil amplification at the site. The findings of this study have important implications for design criteria for sites with similar geology.     

基于小田原试验场地数据的不同基岩输入模型对预测地震动特征的影响研究

随着强震台网的密布及观测记录的增加,为研究各类局部场地地震反应预测模型的合理性提供了有效的参考依据,也使利用强震记录及场地条件研究地震动特征成为可能。选取场地地质参数资料和地震记录数据齐全的日本小田原（Ashigara Valley）盲测试验场地,通过对比不同地震动输入方式及场地反应分析模型,研究地震动特征,分析现有模型的优劣。基于1990年8月5日M5.1强震事件的地表基岩记录和地下基岩地震记录,采用地下台强震记录直接输入、地表基岩台强震记录减半为基底地震动输入、地表基岩台强震记录反演为基底地震动输入作为3种基岩地震动输入。基于局部场地条件分别建立一维等效线性模型、二维黏弹性模型及二维时域等效线性化模型等工程中常用的场地数值分析模型,进行局部场地地震反应分析,预测该盲测场地的地表地震动特征,并与对应的实测强震记录结果进行对比,分析不同基岩地震动输入方式对预测地震动特征及地表土层反应谱特征的影响,重点分析地震动输入、土体非线性、场地横向不均匀性及几何与非线性特征共同作用等因素对地表地震动特征的影响,以期为地表地震动的合理预测提供参考。     

Seismotectonic and liquefaction studies of an industrial site in Northern India

Foundation soil of the proposed fertiliser complex in Northern India is examined for its susceptibility to liquefaction during an earthquake. Information on geotectonic set up and earthquake occurrences in the region around the site is used for defining the earthquake parameters of the ground motion. The effective peak ground acceleration for the site is estimated to be of the order of 0.15 g. Laboratory tests were carried out on soil samples obtained from the site on a horizontal vibration table. The test results were used in determining the possibility of liquefaction employing the methodology developed at the Department of Earthquake Engineering, University of Roorkee. About 10 m of a thick soil layer below the top 1.5 m stiff clay is likely to liquefy. Remedial measures used at the site to counter the possibility of liquefaction are mentioned.     

Estimation of Site Response in Kachchh,Gujarat, India,Region Using H/V Spectral Ratios of Aftershocks of the 2001 <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> 7.7 Bhuj Earthquake

Site response in the aftershock zone of 2001 Bhuj M_w 7.7 earthquake has been studied using the H/V spectral ratio method using 454 aftershocks (M_w 2.5–4.7) recorded at twelve three-component digital strong motion and eight three-component digital seismograph sites. The mean amplification factor obtained for soft sediment sites (Quaternary/Tertiary) varies from 0.75–6.03 times for 1–3 Hz and 0.49–3.27 times for 3–10 Hz. The mean amplification factors obtained for hard sediment sites (hard Jurassic/Mesozoic sediments) range from 0.32–3.24 times for 1–3 Hz and 0.37–2.18 times for 310 Hz. The upper bounds of the larger mean amplification factors for 1–3 Hz are found to be of the order of 3.13–6.03 at Chopadwa, Vadawa, Kavada, Vondh, Adhoi, Jahwarnagar and Gadhada, whereas, the upper bounds of the higher mean amplification factors at 3–10 Hz are estimated to be of the order of 2.00–3.27° at Tapar, Chopadwa, Adhoi, Jahwarnagar, Gandhidham and Khingarpur. The site response estimated at Bhuj suggests a typical hard-rock site behavior. Preliminary site response maps for 1–3 Hz and 310 Hz frequency ranges have been prepared for the area extending from 23–23.85 °N and 69.65–70.85°E. These frequency ranges are considered on the basis of the fact that the natural frequencies of multi-story buildings (3 to 10 floor) range between 1–3 Hz, while the natural frequencies for 1 to 3 story buildings vary from 3–10 Hz. The 1–3 Hz map delineates two distinct zones of maximum site amplification (>3 times): one lying in the NW quadrant of the study area covering Jahwarnagar, Kavada and Gadadha and the other in the SE quadrant of the study area with a peak of 6.03 at Chopadwa covering an area of 70 km × 50 km. While the 3–10 Hz map shows more than 2 times site amplification value over the entire study area except, NE quadrant, two patches in the southwest corner covering Bhuj and Anjar, and one patch at the center covering Vondh, Manfara and Sikara. The zones for large site amplification values (∼3 times) are found at Tapar, Chopadwa, Adhoi and Chobari. The estimated site response values show a good correlation with the distribution of geological formations as well as observed ground deformation in the epicentral zone.     

Seismic microzonation studies for the city of Ragusa (Italy)

The seismic history of the city of Ragusa (Italy), the geotechnical characterisation of the subsoil and the site response analysis should be correctly evaluated for the definition of the Seismic Geotechnical Hazard of the city of Ragusa, through geo-settled seismic microzoning maps. Basing on the seismic history of the city of Ragusa, the following earthquake scenarios have been considered: the “Val di Noto” earthquake of January 11, 1693 (with intensity X–XI on MCS scale, magnitude M_W=7.41 and epicentral distance of about 53 km); the “Etna” earthquake of February 20, 1818 (with intensity IX on MCS scale, magnitude M_W=6.23 and epicentral distance of about 64 km); the Vizzini earthquake of April 13, 1895 (with intensity I=VII–VIII on MCS scale, magnitude M_W=5.86 and epicentral distance of about 26 km); the “Modica” earthquake of January 23, 1980 (with intensity I=V–VI on MCS scale, magnitude M_W=4.58 and epicentral distance of about 10 km); the “Sicilian” earthquake of December 13, 1990 (with intensity I=VII on MCS scale, magnitude M_W=5.64 and epicentral distance of about 50 km). Geotechnical characterisation has been performed by in situ and laboratory tests, with the definition of shear wave velocity profiles in the upper 30 m of soil. Soil response analyses have been evaluated for about 120 borings location by some non-linear 1-D models. Finally the seismic microzonation of the city of Ragusa has been obtained in terms of maps with different peak ground acceleration at the surface; shaking maps for the central area of the city of Ragusa were generated via GIS for the earthquake scenarios.     

Evidence for fault-related directionality and localized site effects from strong motion recordings of the 2003 Boumerdes (Algeria) earthquake: Consequences on damage distribution and the Algerian seismic code

The Algiers–Boumerdes region has been struck by a destructive magnitude 6.8 (M_w) earthquake on May 21, 2003. The study presented in this paper is based on main shock strong motions from 13 stations of the Algerian accelerograph network. A maximum 0.58g peak ground acceleration (PGA) has been recorded at 20 km from the epicenter, only about 150 m away from a PGA of 0.34g, with both a central frequency around 5 Hz, explained by a strong very localized site effect, confirmed by receiver function technique results showing peaks at 5 Hz with amplitudes changing by a factor of 2. Soil amplifications are also evidenced at stations located in the quaternary Mitidja basin, explaining the higher PGA values recorded at these stations than at stations located on firm soil at similar distances from the epicenter. A fault-related directionality effect observed on the strong motion records and confirmed by the study of the seismic movement anisotropy, in agreement with the N65 fault plan direction, explains the SW–NE orientation of the main damage zone. In the near field, strong motions present a high-frequency content starting at 3 Hz with a central frequency around 8 Hz, while in the far field their central frequency is around 3 Hz, explaining the high level of damage in the 3- to 4-story buildings in the epicentral zone. The design spectra overestimate the recorded mean response spectra, and its high corner frequency is less than the recorded one, leading to a re-examination of the seismic design code that should definitively integrate site-related coefficient, to account for the up to now neglected site amplification, as well as a re-modeling of the actual design spectra. Finally, both the proposed Algerian attenuation law and the worldwide laws usually used in Algeria underestimate the recorded accelerations of the 6.8 (M_w) Boumerdes earthquake, clearly showing that it is not possible to extrapolate the proposed Algerian law to major earthquakes.     

The 1986 Dharamsala earthquake of Himachal Himalaya – estimates of source parameters,average intrinsic attenuation and site amplification functions

The 26th April 1986 Dharamsala earthquake (m_b 5.5) occurred in the Kangra region of Himachal Himalaya, which lies in the rupture zone of great Kangra earthquake of 1905. This was the first moderate sized earthquake to be recorded at a few sites of the strong ground motion array in the NW Himalaya. The accelerograms of this earthquake have been used to estimate its source parameters, site amplification functions and to estimate the effective shear wave attenuation factor Q_β in the frontal region of Himachal Himalaya. A double couple fault plane solution for the earthquake has been obtained based on the spectra of the transverse component of the accelerograms. The estimated values of the source parameters are seismic moment: 2.1×10²⁴ dyne-cm, static stress drop (Δσ): 36 bars, source radius (r): 2.8 km and moment magnitude (M_w): 5.4. The estimated average values of effective shear wave attenuation factor Q_β for various sites are in the range of 125 to 300 with an overall spatial average of 239. The influence of local site effects on the observed PGA values have been examined on the basis of site amplification functions.     

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.     

山地与平原地形对地震动影响的对比分析

山地和平原地形在破坏性地震中对震害和地震动有很大影响,不同地形对地震动有不同的地形效应。为研究山地和平原地形对地震动的影响规律,基于集集地震后续M L6.6余震的强震观测记录,对位于不同地形且与震中位于同一直线上的四个强震台记录处理计算,从波形特征、频谱分析、加速度反应谱计算、与规范谱比较等方面进行对比分析,完成长周期地震动反应谱分布的计算。研究表明,山地地形对地震动的高频成分有地形放大效应,沉积平原对地震动的低频长周期部分有放大作用,并且长周期反应谱的峰值区域始终位于沉积平原中部。平原上强震记录产生的长周期放大系数谱已超出规范谱,沉积平原上长周期建筑的抗震设防应引起重视。     

Computation of ground motion amplification in Kolkata megacity (India) using finite-difference method for seismic microzonation

This paper presents the ground motion amplification scenario along with fundamental frequency (F ₀) of sedimentary deposit for the seismic microzonation of Kolkata City, situated on the world’s largest delta island with very soft soil deposit. A 4th order accurate SH-wave viscoelastic finite-difference algorithm is used for computation of response of 1D model for each borehole location. Different maps, such as for F ₀, amplification at F ₀, average spectral amplification (ASA) in the different frequency bandwidth of earthquake engineering interest are developed for a variety of end-users communities. The obtained ASA of the order of 3–6 at most of the borehole locations in a frequency range of 0.25–10.0 Hz reveals that Kolkata City may suffer severe damage even during a moderate earthquake. Further, unexpected severe damage to collapse of multi-storey buildings may occur in localities near Hoogly River and Salt Lake area due to double resonance effects during distant large earthquakes.     

Earthquake ground motion modeling of induced seismicity in the Groningen gas field

A physics‐based numerical approach is used to characterize earthquake ground motion due to induced seismicity in the Groningen gas field and to improve empirical ground motion models for seismic hazard and risk assessment. To this end, a large‐scale (20 km × 20 km) heterogeneous 3D seismic wave propagation model for the Groningen area is constructed, based on the significant bulk of available geological, geophysical, geotechnical, and seismological data. Results of physics‐based numerical simulations are validated against the ground motion recordings of the January 8, 2018, M_L 3.4 Zeerijp earthquake. Taking advantage of suitable models of slip time functions at the seismic source and of the detailed geophysical model, the numerical simulations are found to reproduce accurately the observed features of ground motions at epicentral distances less than 10 km, in a broad frequency range, up to about 8 Hz. A sensitivity analysis is also addressed to discuss the impact of 3D underground geological features, the stochastic variability of seismic velocities and the frequency dependence of the quality factor. Amongst others, results point out some key features related to 3D seismic wave propagation, such as the magnitude and distance dependence of site amplification functions, that may be relevant to the improvement of the empirical models for earthquake ground motion prediction.