2016年阿克陶MS6.7地震震源复杂性与烈度

2016年11月25日在我国新疆克孜勒苏州阿克陶县发生M_S6.7地震（阿克陶M_S6.7地震）.我们收集国内外地震资料,对主震及4级以上余震进行了重新定位和震源机制反演,对434次余震进行了双差定位,对主震震源过程进行了反演确定和复杂性分析,并基于反演确定的有限动态源模型估计了此次地震的烈度分布.结果表明：这次地震发生在当地一个近乎东西向展布的小型盆地内,很可能由一条新断层或隐伏断层的活动所致.发震断层近乎直立,近东西向展布,总体上表现为右旋走滑.破裂首先向西扩展,紧接着向东,随后向东西两个方向同时扩展,然后西侧破裂首先停止,东侧破裂继续,最后破裂在东侧停止,整个过程持续~20 s,释放地震矩1.08×10¹⁹N·m,相当于M_W6.6.破裂过程最终形成两个位错高值区,分别位于初始破裂点的东西两侧,西侧高值区规模较小,东侧区规模较大.根据烈度估计,烈度椭圆长轴方向与主震破裂方向以及余震展布方向一致,最大烈度约为IX度,主要集中在震中以东很小的区域,VIII度区呈纺锤形,分布于震中东西两侧,V至VII度区呈椭圆形,总体上东侧烈度大于西侧.     

国家测震台网中心强震动数据处理系统的设计与实现

为满足国家测震台网海量强震动观测数据的处理时效性、格式标准化、产品丰富度等需求,开发了兼容多类强震动观测站点且具备数据快速汇集、处理及归档等功能的强震动数据处理系统。该系统提供地震波形人机交互数据预处理界面,对加速度记录进行预处理,进一步分析预处理后的加速度事件波形数据,计算得到地震动各项参数,包括峰值加速度(PGA)、峰值速度(PGV)、峰值位移(PGD)、仪器烈度、持时、傅氏谱、反应谱和三联谱等,可以导出地震元数据、地震记录波形,对各类数据进行归档存储。该系统具有平台统一性、功能集成性、数据完备性等特点,有效提升了日常数据处理和管理能力,能在地震应急、震害评估和科学研究中发挥实效。     

基于随机振动模型的强地面运动数值模拟及其在广西灵山M6(3/4)地震中的应用

1936年4月1日广西灵山县平山镇东南罗阳山附近发生M6(3/4)地震,该地震是华南大陆自有地震记载以来发生的最大地震。本研究收集整理了灵山M6(3/4)地震的地质资料、活动断层探测最新成果等,选取适当的研究区域,利用随机振动有限断层模型计算区域内网格点的峰值加速度(PGA)及峰值速度(PGV)等参数,并且加入浅层横波速度结构V_s³⁰对地震动参数的影响,最终得到此次地震的地震动分布并分析了地震动特征。本研究将模拟结果与野外调查烈度数据和地震动衰减关系进行对比,结果显示模拟结果与调查烈度值和衰减关系在整体衰减特征、极震区的分布等方面均符合较好,模拟结果可为该地区未来地震危害性评估提供依据。本研究所使用的方法流程亦可应用于本地区地震烈度速报,为震后应急救援及决策指挥提供帮助。     

吉林省四平市伊通满族自治县公主岭市交界4.3级地震强震动记录和分析

北京时间2009年03月20日14时48分,在吉林省四平市伊通满族自治县、公主岭市交界(北纬43.4°,东经124.9°)发生4.3级地震,震源深度约19公里.辽宁省强震动台网有6个强震台站记录到了地面运动,记录的最大加速度为57.68gal.本文对此次地震的强震记录进行了初步分析.     

吉林伊通M4.3地震强震记录分析

2009年吉林伊通M4.3地震,辽宁省及吉林省强震动台网有8个强震台站记录到地面运动,记录的最大加速度为30.8 cm·s-2,初步展现了强震动台网建设的功效.     

水库诱发地震的地震动特征

利用新丰江、美国奥洛维尔和蒙蒂赛洛水库地震极近场地面运动资料，研究了水库诱发地震的地震动峰值加速度衰减。结果表明，水库地震的峰值加速度高，衰减快。利用国内水库地震资料研究了水库地震烈度。水库地震烈度与同级天然地震烈度相比，一般是烈度高，衰减慢。水库地震的地震动反应谱频率高，地震动持续时间短     

Simulation of strong earthquake characteristics of a scenario earthquake (MS7.5) based on the enlightenment of 2022 MS6.9 earthquake in Menyuan

The Menyuan area is an important transportation hub in the Hexi Corridor. The Menyuan M_S6.9 earthquake that occurred on January 8, 2022 had a major impact on the local infrastructure and transportation of this region. Due to the high possibility of similar strong earthquakes occurring in this area in the future, preliminary assessment of the seismic intensity characteristics of destructive earthquakes in this region is essential for effective disaster control. This paper uses the empirical Green′s function (EGF) method as a numerical simulation tool to predict the ground motion intensity of Datong Autonomous County under the action of the scenario earthquake (M_S7.5). Seismic records of aftershocks of the 2016 Menyuan M_S6.4 earthquake were used as Green’s functions for this simulation. The uncertainties associated with various source parameters were considered, and 36 possible earthquake scenarios were simulated to obtain 72 sets of horizontal ground motions in Datong County. The obtained peak ground acceleration (PGA) vs. time histories of the horizontal ground motion were screened using the attenuation relationships provided by the fifth-edition of China's Seismic Ground Motion Parameter Zoning Map and the NGA-West2 dataset. Ultimately, 32 possible acceleration-time histories were selected for further analysis. The screened PGA values ranged from 78.8 to 153 cm/s². The uncertainty associated with the initial rupture point was found to greatly affect the results of the earthquake simulation. The average acceleration spectrum of the selected acceleration-time history exceeded the expected spectrum of a intermediate earthquake, which means that buildings in Datong County might sustain some damage should the scenario earthquake occur. This research can provide reliable ground motion input for urban earthquake damage simulation and seismic design in Datong County. Growing the dataset of small earthquakes recorded in this region will facilitate the large-scale simulation of ground motions under different earthquake scenarios.     

Strong motion envelope modelling of the source of the Chamoli earthquake of March 28, 1999 in the Garhwal Himalaya, India

Garhwal Himalaya has been rocked by two major earthquakes in the span of just eight years, viz. Uttarkashi earthquake of 20th Oct, 1991 and Chamoli earthquake of 28th March, 1999. Chamoli earthquake of March 28, 1999 was recorded at 11 different stations of a strong motion array installed in the epicentral region. The maximum peak ground acceleration (353 cm/s²) was recorded at an accelerograph located at Gopeshwar. The data from eleven stations has been used for comparison with the simulated acceleration envelopes due to a model of the rupture responsible for this earthquake. For simulation of acceleration envelope the method of Midorikawa (1993) has been modified for its applicability to Himalayan region. This method has earlier been used by Joshi and Patel (1997) and Joshi (1999) for the studyof Uttarkashi earthquake of 20th Oct, 1991. The same method has been used for study of Chamoli earthquake. Layered earth crust has been introduced in place of homogeneous one in this method. The model of rupture is placed at a depth of 12 km below the Munsiari thrust for modelling Chamoli earthquake. Peak ground acceleration was calculated from simulated acceleration envelope using layered as well as homogeneous earth crust. For the rupture placed in a layered crust model peak ground acceleration of order 312 cm/s² was simulated at Gopeshwar which is quite close to actually recorded value. The comparison of peak ground acceleration values in terms of root mean square error at eleven stations suggests that the root mean square error is reduced by inclusion of a layered earth crust in place of homogeneous earth crust.     

A PRELIMINARY STUDY ON THE CORRELATIVITY OF SEISMIC HAZARD BETWEEN BEIJING AREA AND XIONG'AN NEW AREA

Xiong'an New Area is established on April 1, 2017 and some non-capital functions from Beijing would be transferred to this new area. As a political, economic, and cultural center of China, Beijing has a highly dense population, buildings, and transportation. With the rapid development of the urban economy, the population and assets exposed to dangerous areas with M ≥ 7 earthquakes have accumulated in an exponential manner, leading to a continuous surge in seismic risk in Beijing. Studying on the correlativity of seismic hazard between Beijing area and Xiong'an New Area is of great significance to judge whether Xiong'an can play a role in the dispersal of seismic risk of Beijing. Using Monte Carlo method to simulate synthetic earthquake sequences, and for each simulated earthquake, the peak ground acceleration data sets on each site of Beijing and Xiong'an can be calculated through the attenuation relationship. Based on the statistical analysis of the ground motion peak acceleration data sets, this paper holds that the correlativity of the ground motion between Beijing area and Xiong'an New Area is not high; the probability of Beijing and Xiong'an suffering at the same time from an exceeded fortification level of ground motion effect is very low; the probability of Beijing and Xiong'an suffering at the same time from a rare ground motion effect is extremely low. Through defibering population, assets, and setting up a remote backup of earthquake emergency equipment and supplies, Xiong'an New Area can disperse the high seismic risk of Beijing to some extent.     

甘肃岷县漳县6.6级地震强震动观测记录与初步分析

2013年7月22日甘肃省定西市岷县漳县交界(东经104.2°,北纬34.5°)发生M6.6地震。甘肃强震动台网在该地区覆盖良好,获得了丰富的主震加速度记录。本文收集整理了此次地震中各强震动台站获得的加速度记录资料并进行了基本处理;经统计分析绘出了峰值加速度分布图。     

Modelling of rupture planes for peak ground accelerations and its application to the isoseismal map of MMI scale in Indian region

The rupture plane for an earthquake has been modelledby using the semi empirical technique of Midorikawa(1993). This technique estimates ground accelerationby modelling the rupture process during an earthquake.Modifications in this technique have been made for itsapplication to the Indian region. This has been tested forthe Uttarkashi earthquake of 20th Oct, 1991, India, whichwas well recorded at thirteen stations of installedstrong motion array in this region. After testingseveral possible rupture models, a final model has beenselected and peak ground acceleration due to thismodel is simulated at thirteen different stations.Dependency of methodology on model parameters, e.g.dip and mode of rupture propagation have also beenstudied in detail.Using this technique synthetic isoseismal maps wereprepared by converting peak ground acceleration intoMMI scale. Dependency of rupture models on syntheticisoseismals has also been studied in detail. Usingthis method, peak ground acceleration for the Laturearthquake of Sept 30, 1993 has been obtained atvarious places within meisoseismal area. Synthetic andfield intensity was compared at various well-knownsites. Since the region was not covered by anyinstrumental array during Latur earthquake, thesimulated peak ground accelerations are expected toserve basis of design criteria in this region.     

Broadband ground motion simulation using a hybrid approach of the May 21, 2021 M7.4 earthquake in Maduo,Qinghai, China

In this study, the broadband ground motions of the 2021 M7.4 Maduo earthquake were simulated to overcome the scarcity of ground motion recordings and the low resolution of macroseismic intensity map in sparsely populated high-altitude regions. The simulation was conducted with a hybrid methodology, combining a stochastic high-frequency simulation with a low-frequency ground motion simulation, from the regional 1-D velocity structure model and the Wang WM et al. (2022) source rupture model, respectively. We found that the three-component waveforms simulated for specific stations matched the waveforms recorded at those stations, in terms of amplitude, duration, and frequency content. The validation results demonstrate the ability of the hybrid simulation method to reproduce the main characteristics of the observed ground motions for the 2021 Maduo earthquake over a broad frequency range. Our simulations suggest that the official map of macroseismic intensity tends to overestimate shaking by one intensity unit. Comparisons of simulations with empirical ground motion models indicate generally good consistency between the simulated and empirically predicted intensity measures. The high-frequency components of ground motions were found to be more prominent, while the low-frequency components were not, which is unexpected for large earthquakes. Our simulations provide valuable insight into the effects of source complexity on the level and variability of the resulting ground motions. The acceleration and velocity time histories and corresponding response spectra were provided for selected representative sites where no records were available. The simulated results have important implications for evaluating the performance of engineering structures in the epicentral regions of this earthquake and for estimating seismic hazards in the Tibetan regions where no strong ground motion records are available for large earthquakes.     

Stochastic finite-fault simulation of ground motion from the August 11, 2012, <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> 6.4 Ahar earthquake,northwestern Iran

In this study, the 11 August 2012 M _w 6.4 Ahar earthquake is investigated using the ground motion simulation based on the stochastic finite-fault model. The earthquake occurred in northwestern Iran and causing extensive damage in the city of Ahar and surrounding areas. A network consisting of 58 acceleration stations recorded the earthquake within 8–217 km of the epicenter. Strong ground motion records from six significant well-recorded stations close to the epicenter have been simulated. These stations are installed in areas which experienced significant structural damage and humanity loss during the earthquake. The simulation is carried out using the dynamic corner frequency model of rupture propagation by extended fault simulation program (EXSIM). For this purpose, the propagation features of shear-wave including \( {Q}_s \) value, kappa value \( {k}_0 \), and soil amplification coefficients at each site are required. The kappa values are obtained from the slope of smoothed amplitude of Fourier spectra of acceleration at higher frequencies. The determined kappa values for vertical and horizontal components are 0.02 and 0.05 s, respectively. Furthermore, an anelastic attenuation parameter is derived from energy decay of a seismic wave by using continuous wavelet transform (CWT) for each station. The average frequency-dependent relation estimated for the region is \( Q=\left(122\pm 38\right){f}^{\left(1.40\pm 0.16\right)}. \) Moreover, the horizontal to vertical spectral ratio \( H/V \) is applied to estimate the site effects at stations. Spectral analysis of the data indicates that the best match between the observed and simulated spectra occurs for an average stress drop of 70 bars. Finally, the simulated and observed results are compared with pseudo acceleration spectra and peak ground motions. The comparison of time series spectra shows good agreement between the observed and the simulated waveforms at frequencies of engineering interest.     

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.     

2016年新疆阿克陶地震Sentinel-1 InSAR同震形变特征

2016年11月25日新疆克孜勒苏州阿克陶县发生M_W6.6地震。 本文利用合成孔径雷达差分干涉测量技术, 对Sentinel-1卫星获取的升、 降轨雷达数据进行了处理, 提取了该次地震的同震形变场, 并结合形变场特征与震源机制解, 采用梯度下降法反演发震断层的滑动分布。 结果表明, 升、 降轨LOS向同震形变场在发震断层两侧具有明显不同的形变特征, 主要形变区域分布在断层两侧, 升轨LOS向形变量可达-8.2 cm与11.2 cm, 降轨LOS向形变量可达-21.4 cm与13.1 cm; 反演的升、 降轨干涉形变场与InSAR测量值之间的残差得到有效控制, 大部分的残差介于±5 cm之间; 断层滑动分布主要集中于沿断层面深约2~18 km处, 最大滑动量位于沿断层面深约7 km处可达0.96 m; 平均滑动角约182.29°, 最大滑动处的滑动角约197.13°, 两个滑动分布中心的滑动角均接近180°, 表明阿克陶地震为一典型的右旋走滑破裂性事件; 当剪切模量取32 Gpa时, 反演的发震断层地震矩M₀可达9.75×10¹⁸, 相当于矩震级M_W6.60, 与地震波形反演结果一致。     

GROUND MOTION CHARACTERISTICS OF GAOYOU-BAOYING M4.9 EARTHQUAKE

An M4.9 earthquake occurred at the junction of Gaoyou and Baoying on July 20, 2012. In this paper, 43 sets of strong motion records of the main shock are analyzed. With these data, we analyzed the characteristics of the peak ground motion value, attenuation relation, duration and acceleration response spectrum. We draw the peak acceleration contour map of the region near the epicenter. The contour line is smooth and the trend of long axis is northwest-southeast. Distribution of peak acceleration of the observed records is basically consistent with the real intensity distribution. Compared with the predicted result based on the seismic attenuation relation proposed by Yu Yanxiang and Wang Suyunon for eastern China and the Fifth-generation ground motion zonation map, the horizontal PGA and PGV of Gaoyou-Baoying earthquake are higher than the predicted results that are based on the model of Fifth-generation ground motion zonation map, while the PGV is similar with the predicted results which are based on Yu Yanxiang and Wang Suyun's model. We regressively analyzed the spatial-temporal change curves of the two types of relative ground motion durations. Compared with the predicted results proposed by Bommer et al. (2009) based on the NGA strong motion records, the durations of all the three components of this earthquake are higher. 10 typical recordings' acceleration response spectra with 5% damping are calculated, their peak periods are around 0.1~0.3s. The acceleration response spectrum of the station 32BYT, which has the largest amplitude, is considerably larger than the Chinese code design spectra, while it becomes notably smaller when the period is larger than 0.4s. Compared with the horizontal bedrock acceleration response spectrum predicted by the attenuation relationship for the eastern part of China, the observed response spectrum shape is similar with the predicted ones, while almost all the observed response spectrum values (except station 32YCT)are smaller than the predicted bedrock acceleration response spectrum. These phenomena suggest that this earthquake has a weak impact on the seismic fortification standards in this area. Using H/V single-station spectral ratio method, amplitude and site amplification effect of the two typical stations are calculated, and the results show the H/V values are obviously larger than that of ground microtremor. This suggests that the site of the station has obvious amplification effect on ground motion.     

2016年新西兰MS8.0地震震源及强地面运动特征分析

2016年11月13日新西兰发生M_S8.0地震,该地震造成2人死亡,20多人受伤,仅10余栋房屋严重受损,致灾特征较轻. 本文针对新西兰M_S8.0地震的震源参数特征,对该地震的震源过程进行了分析. 结果表明,该地震的地震波辐射能量和视应力均偏低,震源破裂过程为应力上调模式,发震断层破裂相对充分,余震相对丰富.结合新一代衰减关系(NGA)的分析结果显示,该地震的强地面运动峰值加速度(PGA)观测记录整体偏低,与基于震源参数对该地震PGA的理论估算结果较为一致. 进一步选取10 km范围内6个台站的强震观测记录进行比较研究,初步认为影响该地震近断层区域PGA高值的因素主要是断层破裂面上最大滑动集中区的位置,而非主震的初始破裂位置. 综合研究表明,新西兰M_S8.0地震的强地面运动主要受到应力上调模式和最大滑动集中区位置的影响,致灾特征较轻则可能源于峰值加速度偏低和地表破坏较大地区的人口相对稀少.      

Strong ground motion simulation of the 2003 Bam, Iran, earthquake using the empirical Green’s function method

The 2003 Bam, Iran, earthquake caused catastrophic damage to the city of Bam and neighboring villages. Given its magnitude (M _w ) of 6.5, the damage was remarkably large. Large-amplitude ground motions were recorded at the Bam accelerograph station in the center of Bam city by the Building and Housing Research Center (BHRC) of Iran. We simulated the Bam earthquake acceleration records at three BHRC strong-motion stations—Bam, Abaraq, and Mohammad-Abad—by the empirical Green’s function method. Three aftershocks were used as empirical Green’s functions. The frequency range of the empirical Green’s function simulations was 0.5–10 Hz. The size of the strong motion generation area of the mainshock was estimated to be 11 km in length by 7 km in width. To estimate the parameters of the strong motion generation area, we used 1D and 2D velocity structures across the fault and a combined source model. The empirical Green’s function method using a combination of aftershocks produced a source model that reproduced ground motions with the best fit to the observed waveforms. This may be attributed to the existence of two distinct rupture mechanisms in the strong motion generation area. We found that the rupture starting point for which the simulated waveforms best fit the observed ones was near the center of the strong motion generation area, which reproduced near-source ground motions in a broadband frequency range. The estimated strong motion generation area could explain the observed damaging ground motion at the Bam station. This suggests that estimating the source characteristics of the Bam earthquake is very important in understanding the causes of the earthquake damage.     

盆地构造地震地面运动速度和加速度分布特征的数值模拟

本文用错格实数傅里叶变换的拟谱法的数值模拟方法分析了地震波在冲积扇、盆地等不均匀地震构造体区域的传播过程和地面运动分布. 结果表明, 地震波由岩石区进入盆地结构后,在盆地内上下多次反射振荡,对地面建筑物可能形成多次连续的振动和破坏,仅有极少量地震波能量返回岩石区域中,这是防灾研究中值得注意的地面运动特征;地震波在盆地边界地质构造条件下,形成的地震波体波与次生面波动的叠加干涉形成了大振幅的地面运动,它可能导致建筑物的极大破坏;破坏峰值的空间位置可能远离岩石和盆地沉积层的边界或者地震断层的位置.     

Strong ground motion from the November 12, 2017, <Emphasis Type="Bold">M</Emphasis> 7.3 Kermanshah earthquake in western Iran

In this paper, we analyzed the strong ground motion from the November 12, 2017, Kermanshah earthquake in western Iran with moment magnitude (M) of 7.3. Nonlinear and linear amplification of ground motion amplitudes were observed at stations with soft soil condition at hypocentral distances below and above 100 km, respectively. Observation of large ground motion amplitudes dominated with long-period pulses on the strike-normal component of the velocity time series suggests a right-lateral component of movement and propagation of rupture towards southeast. Comparison of the horizontal peak ground acceleration (PGA) from the M 7.3 earthquake with global PGA values showed a similar decay in ground motion amplitudes, although it seems that PGA from the M 7.3 Kermanshah earthquake is higher than global values for NEHRP site class B. We also found that the bracketed duration (D_b) was higher in the velocity domain than in the acceleration domain for the same modified Mercalli intensity (MMI) threshold. For example, D_b reached ~?30 s at the maximum PGA while it was ~?50 s at the maximum peak ground velocity above the threshold of MMI?=?5. Although the standard design spectrum from Iranian Code of Practice for Seismic Resistant Design of Buildings (standard No. 2800) seems to include appropriate values for the design of structures with fundamental period of 1 s and higher, it is underestimated for near-field ground motions at lower periods.