Resolving rupture processes of great earthquakes: Reviews and perspective from fast response to joint inversion

Resolving rupture processes of great earthquakes has fundamental importance to the study of earthquake physics,rupture dynamics, fault zone structure, and evolving processes. It also plays an essential role in earthquake hazard estimation,emergency response and seismic hazard mitigation. This paper reviews the major progress of the earthquake rupture process studies in the last decades, with an emphasize on the research directions of the department geophysics of Peking University including real-time response, back-projection techniques, geodetic data analysis, joint inversion and inversion in complex earth medium. We discussed the advantages and limitations of tradition methods; proposed a systematic and integrated approach from fast-response to detailed study. We also raised perspectives of using source models for ground motion prediction and the possibility of full-dynamic inversion.     

Intensifying Earthquake Prediction in the Light of Action of International Decade for Natural Disaster Reduction

This paper is a review on earthquake prediction and forecast research,progress in earthquake prediction work and pre-estimation of earthquake hazard degree in China in recent years.It indicates that China is the first country,the government of which has promoted and organized the state administration department on reduction of seismic hazards and ensured the socialization of earthquake prediction and forecast in the world.A program of earthquake prevention and protection and hazard reduction based on the results of research on earthquake occurrence regularities and prediction of earthquake preparation trend has been completed,and hence the socialization of results of earthquake prediction and forecast research can be expected to be in practice.The practical seismological works in last 20 years indicate that the earthquakes are not considered to be unpredictable,but it is a challenge remaining to be accepted.We are willing to cooperate with all friends who are engaged in earthquake prediction and forecas     

Potential rupture surface model and its application on probabilistic seismic hazard analysis

Potential sources are simplified as point sources or linear sources in current probabilistic seismic hazard analysis (PSHA) methods. Focus size of large earthquakes is considerable, and fault rupture attitudes may have great influ-ence upon the seismic hazard of a site which is near the source. Under this circumstance, it is unreasonable to use the simplified potential source models in the PSHA, so a potential rupture surface model is proposed in this paper. Adopting this model, we analyze the seismic hazard near the Chelungpu fault that generated the Chi-Chi (Jiji) earthquake with magnitude 7.6 and the following conclusions are reached. 1 This model is reasonable on the base of focal mechanism, especially for sites near potential earthquakes with large magnitude; 2 The attitudes of poten-tial rupture surfaces have great influence on the results of probabilistic seismic hazard analysis and seismic zoning.     

Prediction and verification of earthquakes induced by the Xiluodu hydropower station reservoir

Research has been conducted on reservoir-induced earthquakes in China since the Xinfengjiang reservoir-induced earthquakes in the 1960s. Regulations now require the risk of reservoir-induced earthquakes to be evaluated in the pre-research stage of all hydropower projects. Although nearly 40 cases of reservoir-induced earthquakes have been reported in China, analyses comparing the changes in seismic activity following reservoir impoundment with predictions are rare. In this study, we compared seismic activities observed in the reservoir area before and after the impoundment of the Xiluodu hydropower station in terms of the spatial distribution, frequency, and focal depths of the earthquakes, and clarified the correlation between their frequency/timing and reservoir level after impoundment. We then concluded that the seismic events in the head region were karst-type earthquakes, while those in the second segment of the reservoir were tectonic earthquakes. The spatial distribution of the earthquake epicenters and the seismic intensities validated some of the results for the reservoir-induced seismic risk assessment for the Xiluodu hydropower station, indicating that the proposed earthquake triggers and predictive models are reasonable. This study can provide a valuable reference for investigating the mechanism (s) of reservoir-induced earthquakes, revising reservoir-induced earthquake hazard assessment codes, and predicting the hazard zones of reservoir-induced seismicity under similar conditions.     

Simulations of Seismic Hazard for the Pacific Northwest of the United States from Earthquakes Associated with the Cascadia Subduction Zone

-- We investigate the impact of different rupture and attenuation models for the Cascadia subduction zone by simulating seismic hazard models for the Pacific Northwest of the U.S. at 2% probability of exceedance in 50 years. We calculate the sensitivity of hazard (probabilistic ground motions) to the source parameters and the attenuation relations for both intraslab and interface earthquakes and present these in the framework of the standard USGS hazard model that includes crustal earthquakes. Our results indicate that allowing the deep intraslab earthquakes to occur anywhere along the subduction zone increases the peak ground acceleration hazard near Portland, Oregon by about 20%. Alternative attenuation relations for deep earthquakes can result in ground motions that differ by a factor of two. The hazard uncertainty for the plate interface and intraslab earthquakes is analyzed through a Monte-Carlo logic tree approach and indicates a seismic hazard exceeding 1 g (0.2 s spectral acceleration) consistent with the U.S. National Seismic Hazard Maps in western Washington, Oregon, and California and an overall coefficient of variation that ranges from 0.1 to 0.4. Sensitivity studies indicate that the paleoseismic chronology and the magnitude of great plate interface earthquakes contribute significantly to the hazard uncertainty estimates for this region. Paleoseismic data indicate that the mean earthquake recurrence interval for great earthquakes is about 500 years and that it has been 300 years since the last great earthquake. We calculate the probability of such a great earthquake along the Cascadia plate interface to be about 14% when considering a time-dependent model and about 10% when considering a time-independent Poisson model during the next 50-year interval.     

Rupture model of the 2013 M_W 6.6 Lushan (China) earthquake constrained by a new GPS data set and its effects on potential seismic hazard

Vertical records are critically important when determining the rupture model of an earthquake, especially a thrust earthquake. Due to the relatively low fitness level of near-field vertical displacements, the precision of previous rupture models is relatively low, and the seismic hazard evaluated thereafter should be further updated. In this study, we applied three-component displacement records from GPS stations in and around the source region of the 2013 MW6.6 Lushan earthquake to re-investigate the rupture model.To improve the resolution of the rupture model, records from both continuous and campaign GPS stations were gathered, and secular deformations of the GPS movements were removed from the records of the campaign stations to ensure their reliability. The rupture model was derived by the steepest descent method(SDM), which is based on a layered velocity structure. The peak slip value was about 0.75 m, with a seismic moment release of 9.89 × 10~(18) N·m, which was equivalent to an M_W6.6 event. The inferred fault geometry coincided well with the aftershock distribution of the Lushan earthquake. Unlike previous rupture models, a secondary slip asperity existed at a shallow depth and even touched the ground surface. Based on the distribution of the co-seismic ruptures of the Lushan and Wenchuan earthquakes, post-seismic relaxation of the Wenchuan earthquake, and tectonic loading process, we proposed that the seismic hazard is quite high and still needs special attention in the seismic gap between the two earthquakes.     

走滑型断层地表永久位移预测模型研究

论述了断层地表永久位移概率地震危险性分析方法的发展过程和应用现状,指出了断层地表破裂预测模型研究的意义。根据我国最近的活动断层地表破裂震后调查结果,基于走滑型断层地表破裂数据,拟合了地表破裂预测模型。采用地震活动性模型和地表破裂预测模型,计算了则木河断裂地表迹线上的永久位移危险性曲线,对预测模型在工程场地地震危险性分析中的应用,提出了建议。     

2008年汶川大地震单侧破裂过程的动力学机制研究

2008年汶川大地震的破裂过程极其发杂,向东北方向的破裂距离长达300 km,而向西南方向的破裂长度很小,呈现出单侧破裂的主要特征.尽管汶川地震破裂呈单侧传播的现象引起许多地震学家的关注,但其物理机制至今还不是十分清楚.本文利用有限单元计算方法,模拟了汶川地震的破裂过程.模型中根据龙门山断裂带两侧（东南侧为四川盆地,西北侧为川西高原）实际的地震波速度来确定模型的介质物性参数,利用目前观测的应力环境来选定初始应力条件.模拟结果表明：破裂在汶川地震的震中处成核后,先向断层两侧自发传播,但向东北方向的传播距离明显大于向西南方向;断层面上的正应力在东北方向（破裂的正方向）随着传播距离的增大而不断减小,位错速率随着破裂的传播距离而越来越大,其脉冲变得越来越尖锐,即产生了Weertman脉冲.研究结果显示：由于这种脉冲的出现,破裂在正方向上（东北方向）能够自己放大、自己愈合、自行维持,摩擦热极小,所以破裂能够沿着东北方向一直传播,直到应力场方位发生变化,不利于破裂时才最后终止.但在西南方向,破裂过程中断层面上的正应力增大,阻碍破裂继续扩展.最后就出现了汶川地震中破裂朝东北方向单侧优势传播的基本格局.模拟结果还表明：若断层面两侧介质均匀,则破裂向两侧是对称传播,且破裂距离很短,因此这种情况无法产生像汶川大地震那样的特大地震.因此,文中的模拟结果表明龙门山断裂带两侧的物性差异是造成汶川大地震单侧传播的决定性因素.断层两侧物性差异（bimaterial contrast）影响断层破裂过程的研究对于深入认识地震动力学过程、地震灾害预测及评估等有重要的科学意义.     

地震动衰减模型在玛多Mw7.3级地震影响场快速评估中的应用

2021年5月22日青海省玛多县发生M_w7.3级地震。震后,根据初步估计的断层走向和破裂长度,基于YU15地震动衰减模型和三种NGA-West2（Next Generation Attenuation-West2）地震动衰减模型快速产出地震区震动图及理论烈度图。在获得强震记录和地表破裂长度信息后,对预测结果进行修正。通过比较理论烈度与调查烈度,并结合震动图分布形态以及衰减模型在2016年新疆呼图壁M_w6.0地震中的应用情况对四种地震动衰减模型的适用性进行了分析。结果表明:在台网稀疏地区,基于地震动衰减模型可在震后快速获得地震动分布,并产出具有应用价值的地震影响场;NGA-West2模型在断层破裂较长的大震中表现优于YU15模型,而在中强地震中后者适用性更强;近实时强震动记录可用来检验模型的适用性并对预测结果进行修正;断层破裂尺度、震源机制和破裂过程等信息的准确估计可有效提高地震影响场预测精度。     

设定地震的分类及其确定方法

讨论了震害预测中的设定地震,将其分为检验型设定地震、预测型设定地震和重现型设定地震,并论述了３种设定地震的确定原则和方法．作者还进一步论述了３种设定地震在不同震害预测方法中的应用     

我国地震预报研究近十年的发展与展望

简要回顾了我国地震预报的发展过程，并着重对最近10年（尤其是90年代以来）的地震预报研究进展作了较系统的阐述．主要分为：（1）大地震震源区地壳细结构和介质物性的探测研究；（2）地震和前兆数字化观测技术的研制；（3）地震短临预报方法及其判据、指标的研究；（4）地震序列类型和震后趋势早期应用判断方法及其指标的研究；（5）华北北部地震短临预报方法、前兆综合解释模型与短临预报综合判别系统的研究；（6）地震预报新技术、新方法的探索；（7）地震中长期预报研究及我国近期（10年左右尺度）地震重点监视防御区的确定；（8）板内地震孕育发生的物理模型及前兆场物理解释的探索研究．最后还对我国地震预报的未来发展进行了展望．     

近年我国震后趋势判定与后续强震预测研究进展

详细论述了近年来我国在震后趋势判定与后续强震预测研究中取得的重要进展：（１）对近年发生的重大社会影响的一系列６级以上强震的震后趋势后判断迅速,正确,取得了突出的社会减灾实效;（２）对中强地震发生后后续强震的短临预报取得明显进展,曾较成功地预报１９９５年７月１２日孟连西７．３级强震以及１９９７年４月伽师４次６级以上强震;（３）在国外大震现场的国际合作了开创了震后趋势判定的成功范例;（４）后续强震预测     

论发震构造特性在潜在震源区参数确定中的应用

发震构造特性是潜在震源区划分及其地震年发生率确定的重要依据。潜在震源区除了反映“未来具有发生破坏性地震的地区”的内涵外,还应反映高震级档地震具有相似复发特征的涵义。由于在地震活动性参数统计单元内,有一些具有不同本底地震的活动构造块体,为更好地反映地震活动的空间不均匀性,考虑潜在震源区的三级划分是有必要的。通过分析潜在震源区内高震级档地震的复发特征,计算预测时段内潜在震源区的高震级档地震的发震概率,采用预测时段内概率等效转换获得地震年平均发生率的方法,有助于在中国地震危险性分析框架内考虑潜在震源区的强震复发特性。另外,文中还对潜在震源区内特征地震次级震级档频度不足的特性和发震构造上强震非均匀性在地震危险性分析中的应用问题进行了探讨     

设定地震及其烈度影响判别

设定地震常用于震害预测、地震小区划和重大工程选址,烈度衰减模型反映了地震引起的地面震动及其影响的强弱程度分布。设定地震包括确定性和非确定性设定两种方法,确定性方法基于构造或历史地震,非确定性方法是基于概率危险性方法,用于估计区域或城市未来可能遭遇的地震危险。缺失等震线或震害记载不详的历史地震和概率设定地震都不能确切地反映地震破坏影响,借助于烈度衰减关系模型和GIS,可直观地判别其影响分布情况,便于设定地震的取舍。     

地震潮汐触发

地震潮汐触发是指由潮汐引起的应力变化而触发地震的现象,对地震潮汐触发的研究可能为断层错动引发地震的条件提供有价值的线索,一直受到国内外研究者的关注,本文对相关研究方法和研究结果进行了简要介绍和评述。在早期,地震潮汐触发研究主要是在地震活动中寻找潮汐周期。随着研究的深入,运用统计检验方法对地震活动与潮汐之间的相关性进行严格检验。时至今日,地震是否被潮汐触发的问题仍然悬而未决。不过,有研究结果表明,在强震发生前数年,震源区地震活动往往受到潮汐触发,如果这个结果被证实,将为地震预测研究提供新的线索。      

基于潜在破裂面源模型的南北地震带南段地震区划研究

南北地震带南段大震活动频繁。已有的研究结果表明,大震近场范围场点的地震危险性与地震破裂面产状及其尺度密切相关。因此,在南北地震带南段需要考虑潜在震源三维特征进行地震危险性分析和地震区划研究。本文在充分搜集大震发震构造资料的基础上,在南北地震带南段构建了考虑震源尺度和产状的潜在震源模型,改进了地震危险性概率分析方法,进而对该地区进行地震区划研究。结果表明,考虑潜在震源三维特征的地震危险性分析结果可以有效地反映南北地震带南段发震构造的产状和尺寸特征,提高地震区划结果的合理性。     

The long-term earthquake prediction for the Kuril–Kamchatka island arc for the April 2016 through March 2021 period,its modification and application; the Kuril–Kamchatka seismicity before and after the May 24, 2013, M 8.3 deep-focus earthquake in the Sea of Okhotsk

This paper discusses results from ongoing research on long-term earthquake prediction for the Kuril–Kamchatka island arc based on the concepts of seismic gaps and the seismic cycle. We developed a forecast for the next 5 years (April 2016 through March 2021) for all segments of the earthquake-generating zone along the Kuril–Kamchatka arc. The 20 segments of the arc were analyzed to develop forecasts of the appropriate phases of the seismic cycle, a normalized parameter of the rate of small earthquakes (A₁₀), the magnitudes of moderate earthquakes that are expected with probabilities of 0.8, 0.5, and 0.15, the maximum expected magnitudes, and the probabilities of great (M ≥ 7.7) earthquakes. We discuss the seismic process in the Kuril–Kamchatka earthquake-generating zone before and after the deep-focus May 24, 2013 M 8.3 earthquake in the Sea of Okhotsk. The results corroborate the high seismic hazard in the area of Petropavlovsk-Kamchatskii and the urgent need for continuing with and expanding the ongoing work of seismic retrofitting and seismic safety enhancement. We quote practical results from applications of the method during 30 years.     

中国海域及邻区地震时间分布特征研究

地震时间分布特征研究是进行地震预测和地震危险性分析的重要基础.以中国海域统一地震目录为基础资料,以指数分布模型、伽马分布模型、威布尔分布模型、对数正态分布模型以及布朗过程时间分布(BPT)模型为目标模型,采用极大似然法估算模型参数.根据赤池信息准则(AIC)、贝叶斯信息准则(BIC)以及K-S检验结果确定能够描述海域地...     

Development of the ElarmS methodology for earthquake early warning: Realtime application in California and offline testing in Japan

In July 2009, the California Integrated Seismic Network concluded a three-year study of earthquake early warning systems in California. Three algorithms were expanded and examined during the study. Here we discuss the history, methodology, and performance of one of the algorithms, ElarmS. Earthquake Alarm Systems, or ElarmS, uses peak displacement and maximum predominant frequency of the P-wave to detect earthquakes and quantify their hazard in the seconds after rupture begins. ElarmS was developed for Northern and Southern California, and now processes waveforms in realtime from 603 seismic sensors across the state. We outline the methodology as currently implemented, present several example events from different regions of California, and summarize the performance in terms of false and missed alarms. ElarmS was also tested offline with a dataset of 84 large magnitude earthquakes from Japan. The results from the Japan dataset were used to create a statistical error model for the algorithm. The model can be used to provide realtime uncertainty estimates at any stage in processing. In August 2009 the CISN embarked on a second three-year study of earthquake early warning. As part of this ongoing research, we identify the technological and methodological challenges facing ElarmS. Telemetry latencies and false alarm rates are two key opportunities for improvement.     

夏垫断裂MW≥7.5地震动的预测

针对夏垫断裂开展了M_W≥7.5地震动预测研究。首先基于全破裂模式设定震源（使其尽可能涵盖夏垫断裂的未知信息）模拟得到夏垫断裂发生M_W≥7.5地震时研究区域内的地面地震动场,进而依据分位数筛选出各场点的地震动空间分布,讨论了包含不确定震源下的加速度峰值和速度峰值的分布特征,结果显示当夏垫断裂发生M_W7.9地震时,通州城区、北京中心城区均会发生强烈的运动。之后对比讨论了仿真震源下M_W7.5地震所引起的地面运动场的空间变化,结果显示对于同等震级而言,两种震源的模拟结果可以相互印证。