Neo-deterministic seismic hazard scenarios for India—a preventive tool for disaster mitigation

Current computational resources and physical knowledge of the seismic wave generation and propagation processes allow for reliable numerical and analytical models of waveform generation and propagation. From the simulation of ground motion, it is easy to extract the desired earthquake hazard parameters. Accordingly, a scenario-based approach to seismic hazard assessment has been developed, namely the neo-deterministic seismic hazard assessment (NDSHA), which allows for a wide range of possible seismic sources to be used in the definition of reliable scenarios by means of realistic waveforms modelling. Such reliable and comprehensive characterization of expected earthquake ground motion is essential to improve building codes, particularly for the protection of critical infrastructures and for land use planning. Parvez et al. (Geophys J Int 155:489–508, 2003) published the first ever neo-deterministic seismic hazard map of India by computing synthetic seismograms with input data set consisting of structural models, seismogenic zones, focal mechanisms and earthquake catalogues. As described in Panza et al. (Adv Geophys 53:93–165, 2012), the NDSHA methodology evolved with respect to the original formulation used by Parvez et al. (Geophys J Int 155:489–508, 2003): the computer codes were improved to better fit the need of producing realistic ground shaking maps and ground shaking scenarios, at different scale levels, exploiting the most significant pertinent progresses in data acquisition and modelling. Accordingly, the present study supplies a revised NDSHA map for India. The seismic hazard, expressed in terms of maximum displacement (Dmax), maximum velocity (Vmax) and design ground acceleration (DGA), has been extracted from the synthetic signals and mapped on a regular grid over the studied territory.     

1998年宁蒗5.3、5.2、6.2级地震的短临预报

１９９８年１０月２日宁蒗、盐源交界地带发生了５．３级地震。接着在原震区附近又相继发生了１０月２７日５．２级,１１月１９日５．０级和６．２级、１９９９年１月３日５．０级系列地震,组成了一次震群事件。四川省地震局的专、群地震队伍对前４次地震作了很好的短临预报,当地政府及时采取了有效的防震减灾措施,取得了良好的社会经济效益,得到了中国地震局和四川省地震局的表彰。     

1997年伽师地震临震预报的经验与启示

１９９７年１￣４月新疆伽师发生罕见的强震群活动。几个月内在十分狭小的区域内发生７次６级强震,是新疆几十年来从未遇到过的。概述了伽师地震３次准确和较为准确的临震预报;总结了地点已知的后续强震临震预报的科学和实践问题;论述了震例总结与现实预报的差异,以及地震孕育的复杂性和地震预报的现实做法。     

Consistency of ground-motion predictions from the past four decades

Due to the limited observational datasets available for the derivation of ground-motion prediction equations (GMPEs) there is always epistemic uncertainty in the estimated median ground motion. Because of the increasing quality and quantity of strong-motion datasets it would be expected that the epistemic uncertainty in ground-motion prediction (related to lack of knowledge and data) is decreasing. In this study the predicted median ground motions from over 200 GMPEs for various scenarios are plotted against date of publication to examine whether the scatter in the predictions (a measure of epistemic uncertainty) is decreasing with time. It is found that there are still considerable differences in predicted ground motions from the various GMPEs and that the variation between estimates is not reducing although the ground motion estimated by averaging median predictions is roughly constant. For western North America predictions for moderate earthquakes have show a high level of consistency since the 1980s as do, but to a lesser extent, predictions for moderate earthquakes in Europe, the Mediterranean and the Middle East. A good match is observed between the predictions from GMPEs and the median ground motions based on observations from similar scenarios. Variations in median ground motion predictions for stable continental regions and subduction zones from different GMPEs are large, even for moderate earthquakes. The large scatter in predictions of the median ground motion shows that epistemic uncertainty in ground-motion prediction is still large and that it is vital that this is accounted for in seismic hazard assessments.     

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.     

山西地区6次MS≥4.5地震前地震活动异常总结

系统梳理了2000年以来山西地区6次M_S≥4.5地震前地震活动异常,结果表明,地震空区/平静、地震条带、显著地震/震群、大同地震窗“开窗”活动等异常在地震发生前具有一定普遍性,且异常基本围绕在震中及附近地区分布,特别是在地震平静/空区、地震条带等异常发展后期出现的显著地震/震群活动,对未来地震发生的地点和时间具有较好的预测意义。异常持续时间与发震间隔统计表明：异常多出现在主震发生前6个月以内,显著地震/震群、大同地震窗“开窗”对未来主震的发生具有短临预测意义。此外,随着区域应力水平的不断增强,在特定敏感地区会发生成组极微震密集活动,监视跟踪这些有别于正常活动背景的极微震活动,对地震短临预测具有一定意义。     

用博弈评分法评估分析中国年度地震危险区的预测效能

自20世纪70年代至今,中国地震局持续召开年度全国地震趋势会商会,对来年中国大陆地区的地震趋势作出综合判断,按照"东部5级、西部6级以上"的目标划定年度地震危险区.因为年度危险区的形状不规则、预测震级不同以及每年发生地震数量有限,难以对其预测效能进行科学、客观的评价,也不利于通过有真正信息增益危险区来积累预测经验.本研究采用新近发展的"博弈评分"(gambling score)方法评估了年度地震危险区预测结果的显著性.博弈评分是利用"参考模型"计算每个危险区预测成功与否的参考概率作为"奖励"标准,根据实际地震对应情况进行加分或减分.采用泊松模型和古登堡/里克特定律作为参考模型,对1990-2003年期间的年度地震危险区的预测效能进行了评估和分析.结果表明,虽然年度预测效能之间的差异较大,但是年度地震危险区预测结果明显优于非均匀泊松模型.这一方面表明年度地震危险区的圈定含有一定的地震前兆信息,得分高的危险区的判定依据将为提高地震预测准确率提供有益信息并积累有效的经验,另一方面也揭示了年度会商结果在前兆观测资料和地震活动背景空间分布知识之间整合的不足,即年度地震危险区的圈定在技术上仍有提高的余地.     

From Earthquake Prediction Research to Time-Variable Seismic Hazard Assessment Applications

The first part of the paper defines the terms and classifications common in earthquake prediction research and applications. This is followed by short reviews of major earthquake prediction programs initiated since World War II in several countries, for example the former USSR, China, Japan, the United States, and several European countries. It outlines the underlying expectations, concepts, and hypotheses, introduces the technologies and methodologies applied and some of the results obtained, which include both partial successes and failures. Emphasis is laid on discussing the scientific reasons why earthquake prediction research is so difficult and demanding and why the prospects are still so vague, at least as far as short-term and imminent predictions are concerned. However, classical probabilistic seismic hazard assessments, widely applied during the last few decades, have also clearly revealed their limitations. In their simple form, they are time-independent earthquake rupture forecasts based on the assumption of stable long-term recurrence of earthquakes in the seismotectonic areas under consideration. Therefore, during the last decade, earthquake prediction research and pilot applications have focused mainly on the development and rigorous testing of long and medium-term rupture forecast models in which event probabilities are conditioned by the occurrence of previous earthquakes, and on their integration into neo-deterministic approaches for improved time-variable seismic hazard assessment. The latter uses stress-renewal models that are calibrated for variations in the earthquake cycle as assessed on the basis of historical, paleoseismic, and other data, often complemented by multi-scale seismicity models, the use of pattern-recognition algorithms, and site-dependent strong-motion scenario modeling. International partnerships and a global infrastructure for comparative testing have recently been developed, for example the Collaboratory for the Study of Earthquake Predictability (CSEP) with test regions in California, Italy, Japan, New Zealand, and the Western Pacific. Algorithms and data bases are operated in a permanently learning and upgrading mode. Future perspectives and research requirements and the feasibility and possible problems encountered with the implementation of earthquake predictions in practice are briefly summarized.     

重庆地震台地倾斜观测资料分析

分析了重庆地震台地倾斜观测资料,结合重庆辖区内10多年来的地震活动,提取地倾斜观测异常,为地震预测提供前兆参考。资料分析认为:对于发生在震中距为80~120km的5.0级以下地震,地倾斜观测地震异常少;5.0级以上的地震,地倾斜观测有异常记录,地震中期异常往往表现为震前几个月倾斜量的较大幅度变化,短期异常多数表现为震前差分值的低值变化(在均值附近变化,幅值较小,持续时间20d左右),部分表现为差分值震前突跳,临震异常主要是震前几天短周期波动变化。     

中国大陆地震震例异常统计与分析

以1986～1999年中国大陆地震震例为基础, 对中国大陆及其不同分区的测震和前兆异常进行了统计与分析。 研究结果认为: 地震异常数目、 异常站点和项次百分比随着震级的增大而增多; 震级高异常出现的时间早, 异常范围大; 中期异常多于短期异常, 短期异常多于临震异常, 目前还缺少较多有效的观测到短临异常的方法。 测震中的频度、 空区、 b值、 缺震、 条带、 地震窗、 前兆震(群)、 地震活动分布、 应变释放测项和前兆中的水化学、 水物理、 地形变、 地电测项是出现较多的异常测项。 不同的区域异常出现率高的测震和前兆测项不同, 显示出区域差异。     

对1989年1991年山西地区3次5级以上地震预报的回顾

1989、1991年，山西发生了大同－阳高6.1、忻州5.1和大同－阳高5.8级3次5级以上地震（不含余震），有对大同－阳高2次地震普提出明确的中期预报意见。山西省地震局在中期预报的基础上，对1991年1月29日忻州5.1级、3月26日大同－阳高5.8级地震分别做出了一定程度的短期或临震预报，本文对这3次地震的预报过程做了简要的回顾。     

Retrospective Prediction of Macroseismic Intensities Using Strong Ground Motion Simulation: The Case of the 1978 Thessaloniki (Greece) Earthquake (M6.5)

The densely populated city of Thessaloniki (Northern Greece) is situated in~the vicinity of active seismic faults, capable of producing moderate to strong earthquakes. The city has been severely affected by such events several times during the last 15 centuries. The most recent event occurred on 20 June 1978 (M6.5) in the Mygdonian graben, with an epicentral distance of about 30 km, causing extended damage in the city, with macroseismic intensities between MSK V+ and VIII+. The majority of buildings affected by the earthquake were of reinforced-concrete typology, typical to many southern European metropolitan areas. The source properties of the normal-faulting causative event and the source-to-city propagation path are well known from previous studies. The soil structure under the metropolitan area of Thessaloniki is assigned NEHRP categories B, C, D on the basis of geotechnical and geologic information and single-station ambient-noise measurements. A finite source model and various rupture scenarios of the June 1978 earthquake are used to perform forward stochastic modeling of strong ground motion in terms of peak ground and spectral acceleration. Rock motion is assessed under the city and it is transferred to the surface in accordance with the respective soil category. A GIS tool is employed to compare the estimated strong-motion parameters with the observed detailed damage pattern induced by the 1978 earthquake. For selected natural periods, a satisfactory correlation is established between macroseismic intensity and peak ground and spectral acceleration, thus encouraging the application of stochastic modeling for generating realistic ground-shaking scenarios in metropolitan areas.     

1997年西藏玛尼7.5级地震2001年青海昆仑山口8.1级地震震前中长期预测

作者在1996年底对中国大陆2005年前6级以上强震危险区（13处）作出了预测，其中对应1997年玛尼7.5级地震、2001年昆仑山口8.1级地震，取得了很好的预报效果。本文分时段：1996年底前、1997年11月前、2001年11月前，分别总结了预测依据及震情跟踪分析。加上对同期（1997－2001年）其他强震预测的实际对应，表明地震活动图像分析方法在分析预报研究中是可以有所作为的。     

云南地区强震的短临中强前震序列特征初探

对1965￣2003年云南地区强震前3个月内发生过的短临中强前震进行整理,将其按距离主震发生时间分为直接前震、临震前震与短期前震三类,然后通过地震序列动态分析地震活动,发现各类前震有以下显著特征:(1)直接前震:3、4级地震时间间隔临震时骤减;将地震序列从时间上分成三等分后,第二时段的地震频次没有明显衰减。(2)临震前震:震后第一天地震的频次较高而第二天相对过低,同时同一台站P波初动符号基本上一致或振幅比稳定在某一范围。(3)短期前震:震后10天前震频度衰减至零;20-30天左右地震频度骤增,同时h、p值也反应出类似的特征,30天h值迅速减小接近1或小于1,而b值基本稳定且均低于区域背景值。     

大震区的长期活动

本文通过对我国内陆地区历史上的一些7级以上大震区地震活动的时、空、强分析后,发现许多大震区时过几十年、几百年乃至近千年,都还有不少的小震活动,有时竟然出现了中强地震,甚至还有与主震震级相当的复发性大地震发生。这里根据我国内陆地区17次8级以上地震和37次7级以上地震区的统计分析后,提出了大震区的长期活动概念。这个概念得到了六十年代中期全国陆续设立区域台网后,特别是1970年以来对大震区的观测事实的支持。 本文除简单地讨论了长期活动的物理意义外,进一步用一些具体事实说明了这种长期活动在趋势预报工作中所具有的现实意义。      

长岛地区分维结构变化与区域地震活动关系及短临前兆特征

本文研究了长岛附近地区地震分维结构及月频度趋势性变化与渤海、山东及邻区大范围中强地震活动的关系,得出:大区中强以上地震活跃阶段,研究区存在分维结构,小震活动频度出现趋势性增强;中强以上地震活跃期间分维值较低,活跃阶段将要结束时分维值较高,进入地震平静阶段分维结构不存在。还研究了长岛小震月频度变化的短临前兆特征,得出长岛区地震活动对本区Ms4级以上地震、200公里范围内4.5级以上地震及主震为4.0级的震群、600公里范围内5.5级以上地震,长岛月频度一般在震前1至5个月内会出现8或14次以上的异常现象。反映了长岛小震活动对监视渤海及邻区大范围中强以上地震活动具有一定的窗口作用。     

Recipe for Predicting Strong Ground Motion from Crustal Earthquake Scenarios

We developed a recipe for predicting strong ground motions based on a characterization of the source model for future crustal earthquakes. From recent developments of waveform inversion of strong motion data used to estimate the rupture process, we have inferred that strong ground motion is primarily related to the slip heterogeneity inside the source rather than average slip in the entire rupture area. Asperities are characterized as regions that have large slip relative to the average slip on the rupture area. The asperity areas, as well as the total rupture area, scale with seismic moment. We determined that the areas of strong motion generation approximately coincide with the asperity areas. Based on the scaling relationships, the deductive source model for the prediction of strong ground motions is characterized by three kinds of parameters: outer, inner, and extra fault parameters. The outer fault parameters are defined as entire rupture area and total seismic moment. The inner fault parameters are defined as slip heterogeneity inside the source, area of asperities, and stress drop on each asperity based on the multiple-asperity model. The pattern of rupture nucleation and termination are the extra fault parameters that are related to geomorphology of active faults. We have examined the validity of the earthquake sources constructed by our recipe by comparing simulated and observed ground motions from recent inland crustal earthquakes, such as the 1995 Kobe and 2005 Fukuoka earthquakes.     

九寨沟MS7.0地震强地震动模拟及漳扎镇地震动强度预测

破坏性地震强度预测可用于工程领域抗震设防以及地震危险性分析评估,是防震减灾中一项很重要的基础工作.为了再现九寨沟地震的地震动强度,评估缺失强震记录的九寨章扎台站的地震动强度,本文用经验格林函数法对九寨沟地震进行了数值模拟.选取了震源周边地震动峰值加速度超过10 Gal的10个强震台站进行模拟.因未得到九寨沟地震的余震,初次尝试将汶川地震和定西地震的余震作为格林函数模拟九寨沟地震.模拟结果整体上可以反映各台站地震动的强度特征,尤其是地震动高频成份拟合较好.模拟值的地震动峰值加速度、时程数据、反应谱等与观测值拟合较好.预测结果显示漳扎镇的地震动峰值加速度值约为180~200 Gal.预测结果也表明在缺少大震的余震记录时,经验格林函数法使用其他大震的余震同样可以再现目标地震的强度特征.本研究也为经验格林函数方法在缺乏小震记录地区的使用积累了经验.最后总结了格林函数的选取标准,为经验格林函数方法来预测未来强震动特征积累了经验.     

四川地区强震前弱震空区的变化特征

本文使用四川台网地震目录全部资料,适当参考邻省地震目录,对于四川省从1970年下半年以来发生的15组5级以上地震(另有5组缺乏完整资料无法讨论),作 R-t 图、D-t 图寻找时间空区,并作相应的平面空区图.通过对比分析,发现6组6级以上和4组5级地震前震中附近出现空区现象.根据形成空区的变化过程可分为先扩(大)后缩(小)型和逐渐扩大的开放型两种类型.平面空区的位置和形状往往与构造展布有关.有5组5级震例,或因震前震中区小震不断,或因震中及其附近地区小震太少,看不出明显的空区.本文利用出现空区的10个震例,总结出空区参数与未来地震的震级的经验关系:M=5.50 lg Rmax-2.331.32 M=2.98 lgT+2.941.06 M=2.52lg(RmaxT)-0.510.96文中将上述结果与大竹政和的结果作了对比.空区不同变化形态可用金森博雄概括的突起点模式加以解释.      

Northridge, California, earthquake of 1994: density of pipe breaks and surface strains

Empirical scaling equations are presented which relate the average number of water pipe breaks per km², , with the peak strain in the soil or intensity of shaking at the site. These equations are based on contour maps of peak surface strain evaluated from strong motion recordings, and observations of intensity of ground shaking and damage following the Northridge, California, earthquake of 17 January 1994. Histograms for the number of pipe breaks per km², n, are presented for several ranges of values of the horizontal peak strain and for several values of the site intensity. The observed distribution of pipe breaks is also used to speculate on possible more detailed geographical distribution of near surface strains in the San Fernando Valley and in central Los Angeles. The results can be used to predict number of pipe breaks in the San Fernando Valley and in Los Angeles, for a scenario earthquake or in probabilistic seismic hazard calculations, considering all possible scenarios that contribute to the hazard and the likelihood of their occurrence during specified exposure. Such predictions will be useful for emergency response planning (as the functioning of the water supply is critical for sanitation and in fighting fires caused by earthquakes), to estimate strains during future and past earthquakes in areas where no strong motion was recorded and in defining design guidelines for pipelines and other structures and structural systems sensitive to strains in the ground.