Visco-elastic stress triggering model of Tangshan earthquake sequence

We calculated the Coulomb failure stress change generated by the 1976 Tangshan earthquake that is projected onto the fault planes and slip directions of large subsequent aftershocks.Results of previous studies on the seismic fail-ure distribution,crustal velocity and viscosity structures of the Tangshan earthquake are used as model constraints.Effects of the local pore fluid pressure and impact of soft medium near the fault are also considered.Our result shows that the subsequent Luanxian and Ninghe earthquakes occurred in the regions with a positive Coulomb fail-ure stress produced by the Tangshan earthquake.To study the triggering effect of the Tangshan,Luanxian,and Ninghe earthquakes on the follow-up small earthquakes,we first evaluate the possible focal mechanisms of small earthquakes according to the regional stress field and co-seismic slip distributions derived from previous studies,assuming the amplitude of regional tectonic stress as 10 MPa.By projecting the stress changes generated by the above three earthquakes onto the possible fault planes and slip directions of small earthquakes,we find that the "butterfly" distribution pattern of increased Coulomb failure stress is consistent with the spatial distribution of follow-up earthquakes,and 95% of the aftershocks occurred in regions where Coulomb failure stresses increase,indicating that the former large earthquakes modulated occurrences of follow-up earthquakes in the Tangshan earthquake sequence.This result has some significance in rapid assessment of aftershock hazard after a large earthquake.If detailed failure distribution,seismogenic fault in the focal area and their slip features can be rapidly determined after a large earthquake,our algorithm can be used to predict the locations of large aftershocks.     

The environment shear stress field for the 1976 Tangshan earthquake sequence

（陈培善，肖磊，白彤霞，王溪莉）Ｔｈｅｅｎｖｉｒｏｎｍｅｎｔｓｈｅａｒｓｔｒｅｓｓｆｉｅｌｄｆｏｒｔｈｅ１９７６Ｔａｎｇｓｈａｎｅａｒｔｈｑｕａｋｅｓｅｑｕｅｎｃｅ￥Ｐｅｉ－ＳｈａｎＣＨＥＮ；ＬｅｉＸＩＡＯ；Ｔｏｎｇ－ＸｉａＢＡＩａｎｄＸｉ－ＬｉＷＡ...     

2010年MS7.1级玉树地震同震库仑应力变化以及对2011年MS5.2级囊谦地震的影响

2010年4月14日青海省玉树藏族自治州发生M_S7.1级地震.和传统的板内地震相比,玉树M_S7.1级地震的余震具有数量少、震级大的特点.研究玉树地震主震与余震之间的关系,对于我们了解余震的发震机理具有十分重要的参考价值.本文利用弹性位错理论和分层岩石圈模型,计算玉树地震引起的同震及震后黏弹松弛应力场变化,讨论M_S7.1级玉树地震对余震分布的影响以及与2011年囊谦M_S5.2级地震之间的触发关系.结果显示,玉树地震导致了四处明显的库仑应力增强的扇区,2010年4月13日至6月17日的870次M_L>1.0级余震主要分布于主震破裂面附近区域以及破裂面东北端的应力增强扇区.分析玉树地震对余震分布的影响时,有效摩擦系数以及计算深度的选取对计算结果的影响较小,是否考虑区域构造应力场的影响较大.考虑区域构造应力场时,占总数86.7%的余震位于库仑应力增强区,地震应力触发理论较好地解释了余震的分布.选取囊谦地震震源机制解的两个节面作为库仑应力计算中的接收断层参数,并且考虑不同黏滞系数下的玉树地震同震及震后黏弹松弛效应,模型计算结果均表明囊谦地震位于玉树地震所导致应力影区,仅依靠地震的静态、震后黏弹松弛应力触发理论,无法解释囊谦地震的发生,说明该次地震可能是一次独立的事件.     

Comparison and analysis of the stress field in the source region of Tangshan and Lancang earthquake sequences

ＣｏｍｐａｒｉｓｏｎａｎｄａｎａｌｙｓｉｓｏｆｔｈｅｓｔｒｅｓｓｆｉｅｌｄｉｎｔｈｅｓｏｕｒｃｅｒｅｇｉｏｎｏｆＴａｎｇｓｈａｎａｎｄＬａｎｃａｎｇｅａｒｔｈｑｕａｋｅｓｅｑｕｅｎｃｅｓＧｕｉ－ＬｉｎｇＤＩＡＯ；（刁桂苓）Ｌｉ－ＭｉｎＹＵ；（于利民）...     

唐山余震活动的分期分区特征

本文研宄唐山余震活动时空分布特征,提出唐山余震活动可以分为四个时期和三个余震活动亚区。文中讨论了这一分布特征与震源区地质构造、深部结构以及主震前后地面形变特征的关系。最后指出,主震后应力场的重新分布以及膨胀-流体流动理论可以作为唐山余震成因的一种物理解释。     

1976年唐山强震群震后库仑应力演化及其与2020年古冶5.1级地震的关系

采用能够综合协调长期变形和震后短期变形的Burgers流变模型,模拟了1976年唐山强震群引起的震后形变场以及同震和震后库仑应力变化。结果显示:1976年唐山强震群中主震的两个破裂面以及滦县和宁河两次强余震均对2020年古冶5.1级地震表现为库仑应力加载。岩石圈粘弹性松弛效应引起的库仑应力变化显示,震后15年前后,库仑应力演化状态呈现显著的差异性:震后15年内,库仑应力变化剧烈;而震后15年后,库仑应力呈现缓慢的稳定增加状态。该过程与唐山强震群余震区地震活动过程相似,可能暗示1976年唐山强震群余震区应力调整过程已基本稳定。     

1976年唐山MS7.8地震余震序列持续时间及对地震危险性分析的意义

1976年7月28日唐山MS7.8大地震对唐山及其周边地区造成了重大的人员伤亡和财产损失. 主震之后约15小时滦县又发生了MS7.1地震; 同年11月15日宁河也发生了MS6.9地震. 唐山MS7.8主震后的余震一直持续至今,使该区域至今保持了与主震前相比具有较高的地震活动性．如何估计余震的持续时间,并进一步将余震从主震目录中去除,一直是地震学中所关注的问题．该文通过对数线性回归和理论计算,从不同角度求取并讨论了1976年唐山MS7.8大地震的余震持续时间．结果表明,由对数线性回归计算得到的余震持续时间约为80 a．而基于Dieterich的余震触发理论所得到的余震持续时间则与区域剪应力变化率有关．区域剪应力变化率可有几种不同方法求得: ① 根据剪应力变化率和静态应力降Delta;tau;e及地震回复周期tr之间的关系求取应力变化率,该方法所得到的余震持续时间约为70——100 a;② Ziv和Rubin对Dieterich的方法进行了修正,给出了通过远场加载速率和断层宽度求取应力变化率, 该方法得到的余震持续时间约为80 a;③ 由背景场地震活动性求取远场剪应力速率, 可以得到该地区二维分布式的余震持续时间,此方法得到的研究区域内余震持续时间为130——160 a．综上,唐山地区余震持续时间约为70——140 a,据此, 该地区现今所发生的地震仍为MS7.8唐山地震所触发的余震．      

关于Coulomb应力变化/扰动作用下的Dieterich余震触发机制的广义解

基于单自由度弹簧-滑块模型,滑移速率和状态相依赖的摩擦关系可用于对断层内部地震成核和断层失稳过程的定量化描述.Dieterich(1994)余震触发理论模型首次给出中强震后区域应力场受到静态应力扰动后所导致的区域地震活动性的时空变化特征,近期研究也表明Dieterich理论模型可进一步推广至依赖时间的地震预测模型的建立.本文从简单直观的断层群体化概念模型出发,推导出了包括静态剪应力和正应力扰动作用下广义Dieterich解.同Dieterich经典解相比,广义Coulomb应力变化:ΔCFF^G=Δτ－(μ₀－α)Δσ取代了Dieterich方程中原有的剪切应力扰动Δτ.从而表明余震发生率R同作用于断层上的正应力的变化(扰动)有着密切的相关性.进一步,我们讨论了传统Coulomb应力变化(扰动)模型在地震预测过程中可能存在的局限性.以1976年M_S7.8唐山大地震的主余震序列为例,采用本文中得到的结果,并结合视时窗分段方法,拟合了该地区地震活动性的时间演化过程.结果表明,除Coulomb应力变化(扰动)的影响外,主震前后加载于断层上的剪应力速率变化可对早期余震发生率产生很大影响.     

Numerical simulation of dynamic Coulomb stress changes induced by M 6.5 earthquake in Wuding, Yunnan and its relationship with aftershocks

Based on the discrete wavenumber method, we calculate the fields of dynamic Coulomb rupture stress changes and static stress changes caused by M6.5 earthquake in Wuding, and study their relationship with the subsequent aftershocks. The results show that the spatial distribution patterns of the positive region of dynamic stress peak value and static stress peak value are similarly asymmetric, which are basically identical with distribution features of aftershock. The dynamic stress peak value and the static stress in the positive region are more than 0.1 MPa and 0.01 MPa of the triggering threshold, respectively, which indicates that the dynamic and static stresses are helpful for the occurrence of aftershock. This suggests that both influences of dynamic and static stresses should be considered other than only either of them when studying aftershock triggering in near field. Foundation item: National Natural Science Foundation of China (40364001) and Joint Seismological Science Foundation of China (605014).     

地震的"序列归属"问题与ETAS模型——以唐山序列为例

历史上发生过强震地区的余震活动可能持续较长时间,而余震序列在何时可被看作正常的"背景地震活动",即"序列归属"问题在地球动力学和地震物理中有重要意义.时-空"传染型余震序列"(ETAS)模型可分离"背景"地震和"丛集"地震,并用概率形式表示作为相应事件的可能性,为考察此问题提供了可能.本文以1976年唐山M_S7.8地震序列为例,对唐山地区1970年以来的M_L4.0以上地震进行了时-空ETAS模型拟合,并以2010年以来发生的3次M_S4.0以上地震为例讨论了它们的"序列归属"问题.研究结果显示,3次M_S4.0以上地震的背景地震概率分别为0.72、0.88和0.76,表明它们作为1976年唐山M_S7.8的余震的可能性较低,更可能为背景地震.     

Numerical simulations of deformation and movement of blocks within North China in response to 1976 Tangshan earthquake

Using methods of discontinuous deformation analysis and finite element (DDA+FEM), this paper simulates dynamic processes of the Tangshan earthquake of 1976, which occurred in the northern North China where its internal blocks apparently interacted. Studies focus upon both the movement and deformation of the blocks, in particular, the Ordos block, and variations of stress states on the boundary faults. The Tangshan earthquake was composed of three events: slipping motions of NNE-striking major fault, NE-striking fault near the northeastern end of the NNE-striking fault, and NW-striking fault on the southeastern side of the NNE-striking fault. Compared with previous studies, our model yields a result that is more agreeable with the configuration of aftershock distributions. A number of data are presented, such as the principle stress field during the earthquake, contours of the maximum shear stress, the strike-slip deformation between blocks near the earthquake focus, time-dependent variations of slips of earthquake-triggered faulting, the maximum slip distance, and stress drops. These results are in accord with the earthquake source mechanism, basic parameters from earthquake wave study, macro-isoseismic line, observed horizontal displacement vectors, etc. The Tangshan earthquake exerted different influences on the adjacent blocks and boundary faults between them, thus resulting in differential movement and deformation. The Ordos block seems to have experienced the small-scale counterclockwise rotation and deformation, but its northeast part, bounded on the east by the Taihangshan and on the north by the Yanshan and Yinshan belts, underwent relatively stronger deformation. The Tangshan earthquake also changed the stress state of boundary faults of the North China, leading to an increase in shear stress and a decrease in normal stress in the NW-trending Zhangjiakou-Penglai fault through Tangshan City and the northern border faults of the Ordos block, and therefore raises the potential risk of earthquake occurrence. This result is supported by the facts that a series of Ms≥ 6 earthquakes took place at the northern margin of the Ordos block after the Tangshan earthquake.     

Temporal and spatial variation characteristics of present stress field in the northern area of North China

From focal mechanism solutions of the earthquakes in the northern part of North China during the period of 2002～2006,the directions of principal stress axes in 4 stress sub-districts are analyzed using a grid test method.The characteristic of present crustal stress directions is discussed.Based on this result and on the focal mechanism solutions calculated for some events in the period of 1977～1998,in combination with some other study results,the temporal variation of present crustal stress directions in the northern part of North China is investigated.The re-sults confirm that the direction of crustal principal stress in some regions had somewhat rotated after the 1976 Tangshan M7.8 earthquake.The mean P axes of the focal mechanism solutions rotated clockwise not only in Tangshan sub-district,but also in Beijing and Xingtai sub-districts after the Tangshan earthquake.In Beijing and Xingtai sub-districts the orientations of principal stress axes in the period of 2002～2006 are consistent with that before the Tangshan earthquake,implying that the stress orientations has rotated back to the state before the Tang-shan earthquake in these two sub-districts.The directions of the mean P axes are nearly E-W in Tangshan sub-dis-trict since the M7.8 earthquake.The present stress field in the sub-district northwest to Beijing,or in the western part of the Zhangjiakou-Bohai fault zone,is relatively stable during the time period concerned in this study.Because of the limitation of data,this paper only states a possible variation of stress field in the northern part of North China in the recent decades.     

岩石圈结构对大地震震后形变的影响——以1976年唐山大地震和2001年昆仑山大地震为例

大地震发生之后通常会诱发一系列的余震序列,对比1976年M_S7.8唐山大地震和2001年M_S8.1昆仑山大地震周边区域的地震事件可以看出,唐山大地震余震活动时间要明显长于昆仑山大地震余震活动时间.余震序列往往与震后形变密切相关,而影响震后形变的因素不仅与地震发震断层和震级有关,同时与岩石圈的结构有关.考虑到唐山大地震的发震区华北地块和昆仑山大地震的发震区青藏高原有着较大的岩石圈结构差异,本文采用PSGRN/PSCMP软件计算了岩石圈分层模型的大地震同震和震后形变,分析了地壳弹性模量、弹性厚度以及黏滞性系数对同震和震后形变的影响,进而讨论了影响唐山地震和昆仑山地震余震序列差异的原因.计算结果显示,震后形变会在黏弹性效应的作用下逐渐调整,震后形变的持续时间与地壳弹性模量、地壳弹性厚度和下地壳黏滞性系数有关.上地壳和下地壳弹性模量越大,震后形变达到稳定值的时间越短,弹性模量对震后形变稳定值影响很小.地壳弹性厚度越大,震后形变达到稳定值的时间越短,当断层面底端深度小于地壳弹性厚度时,地壳弹性厚度的增加会引起震后形变稳定值的减小;下地壳厚度对震后形变达到稳定值的时间和稳定值基本无影响.下地壳黏滞性系数越大,震后形变达到稳定值的时间越长,反之亦然.结合唐山地震区的华北地块和昆仑山地震的青藏高原深部结构发现,两者之间的上地壳弹性模型差别不大,唐山地震区地壳弹性厚度略大于昆仑山地震区,但昆仑山地震区下地壳黏滞性系数明显低于唐山地震区.这些因素均决定了昆仑山地震的震后形变持续时间短（余震时间序列短）而唐山地震的震后形变持续时间长（余震时间序列长）.由此可见,岩石圈结构差异可能是导致唐山地震和昆仑山地震余震序列差异的主要因素之一.     

The Longnan MS5.5 Earthquake on September 12, 2008: A Very Shallow Event Probably Triggered by the Wenchuan Earthquake

The 2008 Wenchuan earthquake has a significant impact on the seismicity of nearby regions. The Longnan earthquake which occurred on September 12,2008 in Gansu Province was out of the aftershock zone. Reliable source parameters are essential for understanding the seismogenic process of this earthquake. Therefore,three approaches are adopted to study the source parameters of this event. The focal mechanism is obtained with the g CAP method that takes non-Double-Couple(non-DC)component into account. The two fault planes are NP1:150°/45°/81° and NP2:342°/45°/98°,while the non-DC component is about 53%. The focal depth is 1. 6 km,which indicates the Longnan earthquake is a shallow event. Furthermore,this result is also in good agreement with results obtained with two other approaches:amplitude spectra of Rayleigh wave and surface displacement from In SAR measurement. To analyze the cause of the event,coulomb failure stress change caused by the Wenchuan earthquake on the Longnan earthquake fault plane is calculated. The result shows that coulomb stress change is 30 k Pa around the Longnan earthquake hypocenter,which exceeds the typical triggering threshold of 10 k Pa. The research indicates that the Wenchuan earthquake probably promote the happening of the Longnan earthquake.     

2020年7月12日河北唐山5.1级地震总结

系统梳理2020年7月12日唐山古冶区5.1级地震前出现的地震活动和地球物理观测等异常,结果如下：①地震活动：震前存在地震平静、小震高频、中等地震成组等中短期异常;②地球物理观测：震前主要为电磁和流体异常,短期异常主要分布在震中100 km范围内;③震后回溯：北京及周边震前存在地震发生率异常,利用其他综合方法主要识别出年尺度异常。此次唐山5.1级地震发生在1976年唐山7.8级地震余震区,而7.8级地震序列地震活动呈非均匀衰减特征。此次唐山5.1级地震震源机制为走滑型破裂,截至7月31日,序列b值为0.70、h值为1.8,序列参数基本正常。综合分析认为,此次地震发生前,地震活动中短期异常较显著,地球物理观测异常占比明显偏低;震前对该区域出现的中期和短期异常有所察觉,但短期异常的预测强度偏低,震后总结时按震前异常预测的地震强度为4—5级,接近实际发生地震的强度。     

Three dimensional viscoelastic simulation on dynamic evolution of stress field in North China induced by the 1966 Xingtai earthquake

Introduction Among the mechanisms of earthquake, the essentials may be that the earthquake is a dynamic process of energy accumulation, reaching critical status and failure in focal region. Some nonlin-ear mechanical studies (Bak, Tang, 1989; Ito, Matsuzaki, 1990) show that earthquakes as a self-organized critical phenomenon may be induced by relative small stress disturbance. WANG, et al (1980) already carried out the numerical simulation on earthquake migration in the beginning of 1980s. …     

Comparative Study of Changes in Stress-drop of the Jinggu MS6.6 and Ludian MS6.5 Earthquake Sequences

The earthquake stress-drop values of two sequences were accurately calculated after taking away the effects due to regional earthquake anelastic attenuation and station site response, using waveform data and seismic phase data of sequences of the Jinggu M_S6.6, and Ludian M_S6.5 earthquakes in Yunnan. These results show that the stress drop with magnitude increases within the scope of this study of magnitude. After eliminating the influence of the magnitude, the average value of stress-drop in the Jinggu sequence is higher than that of the Ludian sequence at the same magnitude range. This may be related to the stress state in different regions. In terms of the changes of time and space of stress-drop, before M_S5.8 strong aftershock, the stress-drop is "slowing down-turning up-keeping a high value" after the mainshock, meanwhile, almost all of the abnormally high stress drop value is distributed around the M_S5.8 strong aftershock, showing that the stress environment in the region was increasing after the mainshock. And after the M_S5.9 strong aftershock, stress-drop rapidly declines to a relatively stable state, meanwhile, the high value of stress-drop is distributed around the strong aftershock, showing that the regional tectonic stress gets more fully release, its stress environment begins to rapidly decrease. For the Ludian sequence without a strong aftershock occurring, the average value of stress drop is lower than that of the Jinggu earthquake sequence at the same magnitude range, while at the same time, the stress-drop of the aftershock sequence almost hasn''t changed much. In the time after the mainshock, combined with the release characteristics of the main energy, the stress in the region is excessively released, the subsequent stress in the region gradually returns to normal. This may be the reason why the activity of Ludian aftershocks significantly was weaker and subsequently there were no strong aftershocks occurred.     

1966年邢台地震引起的华北地区应力场动态演化过程的三维粘弹性模拟

利用三维粘弹性有限元模型，首先研究了1966年3月22日邢台7.2级地震所引起的库仑破裂应力的震时变化，以及百年时间尺度流变效应的动态演化图象及其特征，然后模拟了自该地震开始，经过1976年7月28日唐山7.8级地震，直到1980年的15年间，两次7级以上强地震对华北地区应力场的震时扰动作用及其引起的１年时间尺度的库仑破裂应力的动态变化速率，探讨了一次强地震对潜在的下一次强地震可能的加速触发作用．      

2013年芦山MS7.0地震序列参数的早期 特征: 传染型余震序列模型计算结果

为考察2013年4月20日芦山M_S7.0地震震后序列参数的早期特征, 利用“传染型余震序列”(ETAS)模型和最大似然法进行了参数估计. 设定截止震级M_c=M_L2.0, 拟合时段为震后0.31—24.12天, 计算获得α=1.89, p=1.22, 同时利用最大似然法估计获得b=0.72. 与中国大陆地区其它中强震的余震序列参数的比较表明, 芦山M_S7.0地震序列参数表现为触发次级余震的能力较弱和序列衰减速率较快的特征, 反映出余震区相对较高的应力水平. 为检测结果的稳定性, 设定不同的截止震级M_c以及不同的拟合截止时间, 分别进行参数拟合和参数标准差估计. 结果表明, M_c的选取对α值影响明显, 对p值影响则较小. 此外, 震后10天内获得的参数拟合结果随时间变化较为明显, 而其后各参数变化总体较为平稳.      

2011年日本9.0级大地震的应力触发作用初步研究

选用Hayes和Guangfu Shao等给出的震源破裂模型(2011年3月网上公布)、哈佛CMT目录和日本F-net目录给出的余震目录,使用Coulomb3.2软件,对2011年3月11日发生在日本本州东海岸附近海域的M_W9.0地震序列间的静态库仑应力触发关系进行了初步研究。结果表明:1)3月9日发生的M_W7.2前震破裂分布产生的库仑应力对随后发生的M_W9.0主震存在触发作用;2)M_W9.0主震对余震的触发统计结果表明,选用不同的主震模型、余震目录、等效摩擦系数和震源机制解中不同的节面作为接收断层面时,会得到不同的触发统计结果,该研究得到的主震对余震的触发率最小为56.8%,最大为75.3%;3)利用震源机制解计算库仑应力时,理论上震源机制解的2个节面上的剪切应力是相同的,但在实际计算中可能由于2个节面的非正交或震源机制解结果的取整,导致2个节面上计算出的剪切应力不同,但一般差异很小。由于节面的选取对接收断层面上的剪切应力有一定影响,而对接收断层面上正应力的影响较大,因而会影响到库仑应力的计算结果,因此讨论某个具体余震是否被触发或统计余震触发率时,对接收断层面的选取应加以注意。