Coulomb stress evolution and seismic hazards along the Qilian-Haiyuan fault zone
-
摘要: 祁连—海原断裂带是青藏高原东北缘重要的活动断裂带,调节着青藏高原北东向推挤作用和阿拉善地块的东西向运动。已有地震地质和数值模型结果显示,祁连—海原断裂带目前存在几个强震破裂空段且其上应力积累显著、断层闭锁程度高,2022年1月8日门源MS6.9地震即发生在祁连—海原断裂带西段的断层高闭锁、应力积累显著的破裂空段。为进一步认识祁连—海原断裂带未来的强震危险性,本文收集整理了青藏高原北部的强震破裂模型,并基于分层黏弹性流变模型计算了青藏高原北部1900年以来的强震对祁连—海原断裂带的库仑应力加载。结果显示,祁连—海原断裂带西段的木里—江仓断裂和托莱山断裂以及中段的金强河—老虎山断裂应力增强显著,最大库仑应力加载可达1 MPa以上。显著的强震库仑应力加载、强震破裂空段与已有数值模型给出的高应力积累区域和断层高闭锁区域吻合,这表明祁连—海原断裂带西段的木里—江仓断裂和托莱山断裂以及中段的金强河—老虎山断裂未来地震危险性高,亟需进一步关注。Abstract: The Qilian-Haiyuan fault zone (QHF) is one of the most active faults in the northeastern margin of the Tibetan Plateau, which regulates the northeastward spreading of the Tibetan Plateau and the eastward movement of Alxa block. Previous studies indicate that there are several seismic gaps with high fault locking and high stress accumulation rate along the QHF. The Menyuan MS6.9 earthquake on January 8, 2022 ruptured the western segment of the QHF with high fault locking and high stress accumulation. To further understanding the seismic hazards of the QHF, we investigated the Coulomb stress evolution along the QHF based on a layered viscoelastic lithosphere model and the updated coseismic rupture models of the strong earthquakes in the northern Tibetan Plateau since 1900 based on the geological investigation. Our results show there are two Coulomb stress loading segments along the QHF, they are the Muli-Jiangcang fault and Tuolaishan fault in the west, and the Jinqianghe-Laohushan fault in central of the QHF, with the maximum Coulomb stress loading as large as more than 1 MPa. Moreover, the two segments with high Coulomb stress loading are also consistent with the seismic gap with high fault locking and high interseismic stress accumulation rate, manifesting that the two segments are at high seismic risk that deserves more attention in the further researches.
-
引言
青藏高原东北缘是印度板块俯冲、青藏高原向北东向推挤扩展的最前缘,区域内发育有多条大型活动构造,祁连—海原断裂带即是其中一条以走滑为主的大型活动断裂带。该断裂带向西可延伸至哈拉湖附近的木里—江仓断裂,向东经过门源、古浪、景泰、海原、固原,直至六盘山断裂,全长约980 km (图1)。地震活动资料表明,祁连—海原断裂带强震活跃,仅1900年以来就发生了1920年海原M8½地震和1927年古浪M8地震,这两次M8地震共引起祁连—海原断裂带长约350 km的地表破裂(国家地震局兰州地震研究所,宁夏回族自治区地震局,1980;Guo et al,2020;Han et al,2021)。根据地震地质研究,祁连—海原断裂带还存在三个强震危险区,分别为包含祁连—海原断裂带西段的木里—江仓断裂和拖莱山断裂在内的祁连山中段危险区、包含祁连—海原断裂带中段的金强河—毛毛山—老虎山断裂在内的天祝空区和包含六盘山断裂在内的六盘山南—西秦岭东空区(M7专项工作组,2012)。基于大地测量的断层闭锁反演(Li et al,2017)和基于动力学模拟(石富强等,2018)的断层应力状态分析结果均显示,这些区域断层闭锁强,应力积累速率高,未来强震危险性不容忽视。
图 1 青藏高原北部构造环境和强震活动沿祁连—海原断裂带(F3)的彩色散点为数值模拟给出的断层剪应力积累速率(石富强等,2018)。F1:阿尔金断裂带;F2:祁连山断裂带;F3:祁连—海原断裂带;F4:东昆仑断裂带;F5:柴达木盆地北缘断裂;F6:鄂拉山断裂;F7:日月山断裂;F8:狼山山前断裂;F9:六盘山断裂;F10:西秦岭北缘断裂;F11:甘孜—玉树断裂;F12:青川—平武断裂Figure 1. The tectonic setting and the strong earthquake ruptures of the north Tibetan PlateauThe colored dots are the maximum shear stress rates of the Qilian-Haiyuan fault zone based on the finite element simulations (Shi et al,2018)。F1:Altyn fault zone;F2:Qilianshan fault zone;F3:Qilian-Haiyuan fault zone;F4:East Kunlun fault zone;F5:Northern marginal fault of Qaidam basin;F6:Elashan fault;F7:Riyueshan fault;F8:Langshan piedmont fault;F9:Liupanshan fault;F10:Northern marginal fault of west Qinling;F11:Garze-Yushu fault;F12:Qingchuan-Pingwu fault据中国地震台网中心(2022)记录,2022年1月8日在祁连—海原断裂带青海门源发生的MS6.9地震是继1927年古浪M8地震后,祁连—海原断裂带发生的震级最大的一次地震,该地震造成祁连—海原断裂带长约22 km的地表破裂(袁道阳,2022;Yang et al,2022)。此次门源MS6.9地震同时也是继1990年共和M7.0地震后,南北地震带北段发生的最为显著的一次地震,可能预示着青藏高原东北缘地震活动水平增强。此外,青藏高原东北缘及其周边区域1900年以来还发生了1932年昌马M7.6地震和1954年腾格里和山丹M7以上地震,区域外围发生了2001年昆仑山口西M8.1和2008年汶川M8.0两次特大地震以及1963年阿兰胡M7.0、2017年九寨沟M7.0和2021年玛多M7.4等多次M7以上强震。针对这些强震之间的应力触发关系和青藏高原东北缘的地震危险性,研究人员开展了大量的研究工作(傅征祥等,2001;万永革等,2007;陈为涛等,2013;Xiao,He,2015;孙云强,罗纲,2018;张瑞等,2021),指出青藏高原东北缘的历史强震之间存在显著的应力触发关系,暗示区域活动构造之间存在动力学关联,并据此对该区域的地震危险性进行了讨论分析。
然而随着地震地质研究的不断深入,对于一些历史强震的破裂有了更深入的认识。以1927年古浪M8地震为例,以往研究由于对该地震的同震位错和破裂长度的认识还不够充分,所给出的同震位错模型(傅征祥等,2001;万永革等,2007;Xiao,He,2015)与最新的地震地质研究结果(Guo et al,2020)存在一定差异。此外,2017年九寨沟地震、2021年玛多地震和2022年门源M6.9地震的发生,势必也会引起祁连—海原断裂带的应力变化。为进一步认识祁连—海原断裂带当前应力累积水平和地震危险性,本文拟调研整理该断裂带上历史地震同震位错的最新研究成果,并基于分层黏弹性模型计算青藏高原东北缘及其周边区域1900年以来M7以上强震对祁连—海原断裂带的应力加载,以期讨论祁连—海原断裂带未来的强震危险性。
1. 理论方法及模型参数
一次强震的发生,可以改变周边区域岩石圈的应力状态,进而触发区域地震。基于库仑破裂准则,经过一系列的假定,可以得到我们常见的库仑应力表达(King et al,1994) ΔCFS=Δτ+μ′Δσn,式中:Δτ和△σn分别为指定断层面上的剪应力变化和正应力变化;μ′为断层有效摩擦系数,与断层介质物性、流体侵入等相关,表征断层阻碍两侧地块相对运动的能力,也有研究人员认为其与断层类型(Ali et al,2008)和滑动速率(Parsons,Dreger,2000)相关。He等(2013)基于GPS观测的青藏高原东北缘地壳变形机制模拟研究显示海原断裂带的有效摩擦系数可低至0.05,本文认为该研究基于现有观测资料约束给出的断层有效摩擦系数更接近于真实的断层力学状态。综合考虑青藏高原北部的断层类型和He等(2013)的结果,本文取μ′=0.1。由于M6地震影响的空间范围有限(Verdecchia et al,2018),本文主要考虑青藏高原北部及其周边区域发生的M7以上地震,震源参数详见表1。计算程序采用汪荣江教授开发的PSGRN/PSCMP (Wang et al,2006),库仑应力计算深度参考区域震源平均深度,取10 km。
表 1 青藏高原北部及周边1900年以来的强震同震位错模型Table 1. Coseismic rupture models of the strong earthquakes in northern Tibetan Plateau since 1900发震日期
年−月−日地名 M 震中位置 走向
/°倾向
/°滑动角
/°破裂长度
/km破裂宽度
/kmSS
/mDS
/m来源 东经/° 北纬/° 1920−12−16 海原 8½ 104.10 37.04 110 88 14 23 30 4.39 −1.09 ① 104.50 36.90 112.5 88 14 64 30 6.31 −1.57 104.95 36.74 310.6 88 14 31 30 4.99 −1.24 105.58 36.52 112.6 88 14 39 30 6.21 −1.55 105.96 36.16 148.4 88 14 74 30 6.93 −1.73 1927−05−23 古浪 8.0 101.50 37.69 139 80 0 11 20 3 0 ② 101.58 37.61 122 80 0 11 20 3.3 0 101.69 37.56 116 80 0 12 20 3.6 0 101.81 37.51 122 80 0 11 20 4 0 101.91 37.46 102 80 0 16 20 6 0 102.09 37.43 99 80 0 16 20 7 0 102.27 37.41 92 80 0 14 20 7 0 102.42 37.41 95 80 0 18 20 3 0 102.63 37.39 90 80 0 14 20 2 0 102.25 38.07 139 40 90 42 20 0 −1.50 1932−12−25 昌马 7.6 96.70 39.70 115 79 30 116 20 2.34 −1.35 ① 1936−02−07 康乐 6.8 103.40 35.40 270 70 10 36 15 0.73 −0.13 ③ 1937−01−07 阿兰湖 7.5 97.60 35.50 110 70 15 208 20 3.96 −1.06 ④ 1947−03−17 达日 7.7 99.50 33.30 135 60 60 150 20 2.00 −3.46 ④ 1954−02−11 山丹 7.3 101.30 39.00 290 45 35 90 20 1.45 −1.02 ① 1954−07−31 腾格里 7.0 104.18 38.80 153 83 171 58 15 −0.94 −0.15 ① 1963−04−19 阿兰湖 7.0 97.00 35.70 277 80 −10 68 15 1.16 0.2 ① 1973−07−14 玛尼 7.0 86.48 35.18 81 60 −35 59 15 0.78 0.55 ① 1976−08−16 松潘 7.2 104.08 32.60 165 63 40 30 15 1.10 0.90 ④ 1976−08−23 平武 7.2 104.30 32.50 155 65 40 22 15 1.10 0.90 ④ 1990−04−26 共和 7.0 100.33 36.06 346 78 128 59 15 −0.41 −0.52 ① 1997−11−08 玛尼 7.5 87.33 35.07 76 90 −5 170 20 5 0.44 ⑤ 2001−11−14 昆仑山 8.1 90.54 35.95 99 90 5 346 20 4 −0.35 ⑤ 2008−05−12 汶川* 7.9 103.32 31.00 − − − − − − − ⑥ 2016−01−21 门源* 6.4 101.60 37.70 − − − − − − − ⑦ 2017−08−08 九寨沟* 7.0 103.82 33.20 − − − − − − − ⑧ 2021−05−21 玛多* 7.4 98.30 34.60 − − − − − − − ⑨ 2022−01−08 门源* 6.9 101.26 37.77 − − − − − − − ⑦ 注,*表示这些地震的相关震源参数采用有限断层模型反演给出;SS为走向方向滑动量,左旋为正;DS为倾向方向滑动量,正断为正。最后一列来源:① 万永革等(2007); ② Guo et al (2020) ;③ 梅秀苹等(2012);④Shan et al (2015) ;⑤ 沈正康(2003);⑥Shen et al (2009) ;⑦ 李振洪等(2022);⑧ 张旭等(2017);⑨ USGS (2021)。岩石圈介质物性参数中的弹性参数参考邵志刚等(2008),详见表2。震后形变研究(Broerse et al,2015)表明,岩石圈流变性质对区域地壳形变影响显著,且表现出时间尺度差异性(Huang et al,2014)。本文采用能够协调短期变形和长期变形的伯格斯(Burgers)体模型,其中表征短期变形的开尔文(Kelvin)体黏滞系数ηk参考邵志刚等(2008)关于昆仑山口西MS8.1地震震后变形模拟研究,表征长期变形的麦克斯韦尔(Maxwell)体黏滞系数ηm参考He等(2013)的震间应力场模拟结果,具体见表2。
表 2 青藏高原东北缘岩石圈介质模型参数Table 2. The model parameters of the lithosphere structure in the northeastern margin of the Tibetan Plateau分层 厚度/km vP/(km·s−1) vS/(km·s−1) 密度
/(kg·m−3)ηk/(1018 Pa·s) ηm/(1019 Pa·s) 沉积层 10 5.9 3.41 2500 弹 性 上地壳 10 6.175 3.57 2700 低速层 12 5.85 3.38 2600 中地壳 10 6.4 3.70 3000 6.30 1.00 下地壳 20 6.8 3.93 3100 6.30 1.00 上地幔 − 8.1 4.68 3350 0 10.00 注:ηk为开尔文体黏滞系数,ηm为麦克斯韦尔体黏滞系数。 2. 计算结果分析
2.1 祁连—海原断裂带强震应力触发
前人已对祁连—海原断裂带强震应力触发开展了大量的研究(傅征祥等,2001;万永革等,2007;Xiao,He,2015)。本文结合地震地质现场调查研究成果更新1927年古浪M8地震的同震破裂模型(表1),在此基础上结合2022年门源M6.9地震进一步讨论祁连—海原断裂带的应力触发。结果显示:1920年海原M8½ 地震对1927年古浪M8地震破裂面的同震和震后松弛累积库仑应力加载约为0.01—0.2 MPa (图2a);由于空间距离较远,1920年海原M8½ 地震对2022年门源M6.9地震的累积库仑应力加载约为0.01 MPa (图2b)。地震地质调查发现,2022年门源M6.9地震发生在1927年古浪M8地震破裂区的北端,且二者的地表破裂有一定重合(韩竹军等,2022)。此外,本文计算结果显示1927年古浪M8地震对2022年门源M6.9地震的应力加载最为显著,同震库仑应力加载约0.3 MPa,考虑到近百年的岩石圈流变松弛,同震和震后松弛累积库仑应力加载可达0.75 MPa (图2c,图4a)。
图 2 祁连—海原断裂带三次显著强震之间的应力触发(a) 1920年海原地震对1927年古浪地震的库仑应力加载;(b) 1920年海原地震对2022年门源地震的库仑应力加载;(c) 1927年古浪地震对2022年门源地震的库仑应力加载;(d) 1920年海原地震和1927年古浪地震对2022年门源地震的库仑应力加载Figure 2. Stress interaction among the three strong earthquakes along the Qilian-Haiyuan fault zoneThe magenta lines and beach balls are the current earthquake ruptures and the related focal mechanisms,and the light blue lines and beach balls express the receive faults and the related focal mechanisms. (a,b) The cumulated Coulomb stress changes associated with 1920 Haiyuan earthquake just before 1927 Gulang earthquake and 2022 Menyuan earthquake;(c) The cumulated Coulomb stress changes associated with the 1927 Gulang earthquake just before 2022 Menyuan earthquake;(d) The joint Coulomb stress interaction on 2022 Menyuan earthquake associated with 1920 Haiyuan and 1927 Gulang earthquakes图 4 (a) 青藏高原北部历史强震对2022年门源M6.9地震的累积库仑应力加载;(b) 扣除1920年海原M8.5地震、1927年古浪M8.0地震和1954年山丹M7.3地震应力影响后,2022年门源M6.9地震断层面的库仑应力累积变化Figure 4. (a) The temporal evolution of the Coulomb stress on the rupture plane of the 2022 Menyuan M6.9 earthquake associated with the strong earthquakes in Table 1;(b) Same as Fig. (a),but without the stress loading associated with the 1920 Haiyuan M8.5,1927 Gulang M8.0 and 1954 Shandan M7.3 earthquakes2.2 历史强震对门源M6.9地震的应力加载
以2022年门源M6.9地震震源机制为接收断层,利用表2中的区域岩石圈分层流变模型计算表1所列区域历史强震对2022年门源M6.9地震的应力加载,结果显示1920年海原M8½ 地震、1927年古浪M8地震和1954年山丹M7.3地震对2022年门源M6.9地震的库仑应力加载最为显著(图2,图4a)。扣除这三次地震的应力影响,给出祁连—海原断裂带周边历史强震对2022年门源M6.9地震的库仑应力加载(图3,图4b)。结果显示:发生在祁连—海原断裂带周边的历史强震对2022年门源M6.9地震的同震库仑应力加载相对较低,约为0.002 MPa;与之形成鲜明对比的是,数十年至百年的岩石圈流变松弛,可造成2022年门源M6.9地震震源区约0.006 MPa的库仑应力加载,同震和震后松弛库仑应力加载可达0.008 MPa,接近常见的触发阈值0.01 MPa,这表明祁连—海原断裂带周边历史强震对2022年门源地震的发生同样具有一定的应力触发作用。从应力加载的时序图可见,祁连—海原断裂带周边历史强震对2022年门源地震的应力加载主要来自于1932年昌马M7.6地震、1954年腾格里M7.0地震、1990年共和地震和2008年汶川地震;而1947年达日M7.7地震、1963年阿兰湖M7.0地震和2001年昆仑山口西M8.1地震对2022年门源地震表现为库仑应力卸载,但卸载作用不明显,累积库仑应力卸载约为0.001 MPa。
2.3 祁连—海原断裂带累积库仑应力演化
将祁连—海原断裂带(图1)离散为10 km的小段,根据每段首尾点坐标计算每个断层段的走向,通过最近点插值将地震地质研究给出的不同断层倾角和滑动角赋值到10 km的离散断层段,给出计算祁连—海原断裂带库仑应力的接收断层信息(图5)。基于表1给出的青藏高原北部1900年以来的M7以上强震和表2构建的岩石圈分层流变结构,计算出祁连—海原断裂带1900年以来周边历史强震同震和震后松弛加载的库仑应力累积变化(图5)。结果显示,受1920年海原M8½地震、1927年古浪M8地震和1932年昌马M7.6地震的影响,祁连—海原断裂带目前还存在两个显著的库仑应力加载区域,即祁连—海原断裂带中段的金强河—老虎山断裂和西段的木里—江仓断裂和托莱山断裂,累积库仑应力加载已达到1.0 MPa以上。这两个区域是地震地质研究给出的大震破裂空段(M7专项工作组,2012),同时也是动力学模拟(石富强等,2018)给出的高应力累积区域(图1)和大地测量给出的断层高闭锁区域(Li et al,2017),未来地震危险性值得高度关注。
3. 参数不确定分析
3.1 有效摩擦系数影响
有效摩擦系数是库仑应力计算的不确定性参数,通常采用不同的有效摩擦系数取值讨论库仑应力计算结果的可靠性(Steacy et al,2004;Wan,Shen,2010;Wang et al,2014;石富强等,2020)。本文以老虎山断裂西端的库仑应力时序变化为对象,讨论不同有效摩擦系数取值对库仑应力计算结果的影响,结果如图6所示。可以看到,对于给定的三个有效摩擦系数取值,老虎山断裂西段的库仑应力变化趋势一致,近百年的流变松弛引起的库仑应力变化最大差值约为0.05 MPa。
实际上,有效摩擦系数对库仑应力计算结果的影响源于断层面正应力的变化。当正应力和孔隙应力的变化相较于剪切应力变化较小时,有效摩擦系数的选取一般不会影响库仑应力的计算结果(Pollitz et al,2006);反之,当剪切应力变化相对于正应力和孔隙应力的变化较小时,基于不同有效摩擦系数计算得到的库仑应力之间存在显著的差异(Jia et al,2018,2021)。祁连—海原断裂带乃至整个青藏高原北部的历史强震,特别是M7.5以上特大地震,基本以走滑为主且走向近乎一致(图1),这使得库仑应力计算结果主要以接收断层面的剪应力变化为主,正应力变化相对较小。当断层面剪应力变化绝对值远大于正应力变化时,有效摩擦系数在库仑应力计算中发挥的作用有限,对计算结果影响不大。
3.2 古浪地震位错模型对模拟结果的影响
在本文引言中,我们简单介绍了不同研究给出的1927年古浪M8地震同震破裂模型存在一定的差异。为进一步理解同震破裂模型的选取对模拟结果的影响,本文分别基于Xiao和He (2015)、傅征祥等(2001)、万永革等(2007)和Guo等(2020)给出的1927年古浪地震同震位错模型,计算1927年古浪地震对周边断层的库仑应力影响,结果如图7所示。可见:对于鄂拉山断裂、日月山断裂、云雾山断裂以及黄河—灵武断裂等距离古浪地震震源较远的活动断裂,不同的古浪地震同震位错模型引起的库仑应力变化均以加载为主且量级相当。但对于震源附近的断裂,基于古浪地震的不同位错模型所得的断层库仑应力计算结果存在差异。对于托莱山断裂,基于万永革等(2007)和Guo等(2020)的古浪地震破裂模型得到的断层库仑应力变化表现为加载(图7c,d),而基于Xiao和He (2015)和傅征祥等(2001)的破裂模型表现为库仑应力卸载(图7a,b);对于金强河断裂,结果刚好相反,基于Xiao和He (2015)和傅征祥等(2001)的古浪地震破裂模型表现为库仑应力加载(图7a,b),基于万永革等(2007)和Guo等(2020)的结果表现为库仑应力卸载(图7c,d),与托莱山断裂的情形刚好相反;对于毛毛山断裂和老虎山断裂,基于Xiao和He (2015)和万永革等(2007)的古浪地震破裂模型的结果表现为库仑应力卸载(图7a,c),而基于Guo等(2020)的结果表现为库仑应力加载。此外,对于祁连山断裂、香山—天景山断裂、六盘山断裂和海原断裂,基于不同模型的库仑应力计算结果也存在不同程度的差异。
图 7 古浪地震位错模型对周边主要断裂的库仑应力影响F1:祁连山北缘断裂;F2:香山—天景山断裂;F3:鄂拉山断裂;F4:日月山断裂;F5:海原断裂;F6:云雾山断裂;F7:龙首山断裂;F8:皇城—双塔断裂;F9:金强河断裂;F10:毛毛山断裂;F11:黄河—灵武断裂;F12:木里—江仓断裂;F13:拖莱山断裂;F14:冷龙岭断裂;F15:老虎山断裂;F16:六盘山断裂(a) 基于Xiao和He (2015)的位错模型;(b) 基于傅征祥等(2001)的位错模型;(c) 基于万永革等(2007)的位错模型;(d) 基于Guo等(2020)的位错模型Figure 7. The cumulated Coulomb stress changes on the main active faults associated with different coseismic models of the 1927 Gulang earthquakeF1:Qilianshan northern marginal fault;F2:Xiangshan-Tianjingshan fault;F3:Elashan fault;F4:Riyueshan fault;F5:Haiyuan fault;F6:Yunwushan fault;F7:Longshoushan fault ;F8:Huangcheng-Shuangta fault;F9:Jinqianghe fault;F10:Maomaoshan fault;F11:Huanghe-Lingwu fault;F12:Muli-Jiangcang fault;F13:Tuolaishan fault;F14:Lenglongling fault;F15:Laohushan fault;F16:Liupanshan fault(a) The coseismic rupture model from Xiao and He (2015);(b) The coseismic rupture model from Fu et al (2001); (c) The coseismic rupture model from Wan et al (2007);(d) The coseismic rupture model from Guo et al (2020)总之,Guo等(2020)的研究是在前人地震地质调查工作的基础上,通过进一步的现场调查测量和测年分析而得到的结果,其给出的古浪地震破裂模型更接近于真实,且计算所得的祁连—海原断裂带库仑应力加载区域、断层高闭锁(Li et al,2017)与M7专项工作组(2012)圈定的强震破裂空段一致,因此本文采用Guo等(2020)给出的古浪地震破裂模型。
3.3 1540年冷龙岭断裂M7地震和1888老虎山断裂M7地震对模拟结果的影响
地震地质调查研究(M7专项工作组,2012)指出,祁连—海原断裂带冷龙岭段和老虎山段分别于1540年和1888年发生各发生一次M7地震。由于这两次M7地震的地震地质研究程度较低,尚无详细的震源参数研究结果。本文根据M7专项工作组(2012)给出的破裂范围给出了这两次M7地震的破裂面长度,破裂面宽度假设为20 km,破裂面走向、倾角、滑动角参考图5给出的结果,根据经验回归关系(Wells,Coppersmith,1994)所得平均位错量1.3 m计算了这两次M7地震对祁连—海原断裂带应力的影响(图8a),并将其叠加到图5给出的1900年以来青藏高原北部M7地震对祁连—海原断裂带的应力影响(图8b)。结果显示,这两次地震可以引起祁连—海原断裂带西段的托莱山断裂、中段的金强河—毛毛山断裂和老虎山断裂部分段落约0.1 MPa库仑应力增强。综合考虑1540年和1888年两次M7地震以及表1给出的1900年以来青藏高原北部强震的共同作用,祁连—海原断裂带中段的金强河—老虎山断裂和西段的木里—江仓断裂和托莱山断裂库仑应力增强显著,与基于图5得到的认识一致。
图 8 祁连—海原断裂带库仑应力变化(a) 1540年和1888年两次M7.0地震对祁连—海原断裂带的库仑应力影响;(b) 1540年M7.0地震、1888年M7.0地震和表1中的历史强震引起的祁连—海原断裂带的累积库仑应力变化Figure 8. The Coulomb stress changes of the Qilian-Haiyuan fault zone(a) The Coulomb stress changes associated with the 1540 M7 and 1888 M7 earthquakes;(b) The cumulated Coulomb stress changes of Qilian-Haiyuan fault zone associated with the 1540 M7 and 1888 M7 earthquakes as well as the strong earthquakes in Table 14. 讨论与结论
本文收集了青藏高原北部历史强震同震破裂模型,并结合最新地震地质研究成果对部分历史强震的同震位错模型进行了修正。在此基础上,基于分层黏弹性流变模型计算了祁连—海原断裂带几次强震之间的相互作用以及青藏高原北部1900年以来的强震对该断裂带的库仑应力加载。为保证分析结果的可靠性,本文还补充分析了有效摩擦系数、古浪地震同震破裂模型差异以及祁连—海原断裂带两次离逝时间较近的M7地震对模拟计算结果的影响。
有效摩擦系数的选取是库仑应力计算过程中一个重要的参数,Ali等(2008)曾建议走滑断层的有效摩擦系数取0.2,正断层的有效摩擦系数取0.6,逆冲断层的有效摩擦系数取0.8。本文的参数不确定计算结果显示,对于走向相近的走滑型地震而言,库仑应力计算结果主要受断层剪应力变化的控制,有效摩擦系数对库仑应力计算结果影响不显著,这一认识与万永革等(2007)的认识一致。He等(2013)基于大地测量结果的祁连—海原断裂带动力学模拟结果显示,祁连—海原断裂带断层的有效摩擦系数可低至0.05,本文认为这一结果是对祁连—海原断裂带有效摩擦系数较为可靠的量化表征。因此,本文在计算过程中取有效摩擦系数为0.1,以确保计算分析结果的可靠性。
研究表明,地震活动与断层应力变化密切相关,即在库仑应力增强的地方,地震活动性显著提高,反之亦然(Harris,1998;Stein,1999;Toda et al,2008;Kroll et al,2017)。强震相互作用计算结果显示,1920年海原M8½地震对1927年古浪M8地震的累积库仑应力加载最大可达0.2 MPa,1927年古浪M8地震对2022年门源M6.9地震的累积库仑应力加载可达0.75 MPa,均远大于通常所用的触发阈值0.01 MPa,表明祁连—海原断裂带强震之间的库仑应力触发作用显著。该过程与同为大型走滑断裂的土耳其北部北安托利亚(Anatolian)断层强震触发过程相似(Stein et al,1997)。祁连—海原断裂带周边强震对2022年门源M6.9地震的累积库仑应力加载接近0.01 MPa,且其随时间的演化过程与祁连—海原断裂带M5以上地震活动时序过程相似(图9),进一步表明青藏高原东北缘各活动构造之间存在动力学相关性(Xiao,He,2015),是一个相互作用的动力学系统,2022年门源M6.9地震在孕育过程中也受到了祁连—海原断裂周边强震活动的应力触发作用。
祁连海原断裂带库仑应力的时空演化过程显示:青藏高原北部历史强震对祁连—海原断裂带西段的木里—江仓断裂和托莱山断裂以及中段的金强河—老虎山断裂的应力增强显著,木里—江仓断裂、托莱山断裂和金强河—老虎山断裂的同震库仑应力累积分别为0.000 7—0.013 3 MPa,0.002 3—1.608 MPa和0.094—5.06 MPa,相应各段的平均同震库仑应力加载分别为0.006 45 MPa,0.367 MPa和1.125 MPa。考虑岩石圈的流变松弛作用,木里—江仓断裂、托莱山断裂和金强河—老虎山断裂的最大累积库仑应力变化分别可达0.064 MPa,2.02 MPa和5.75 MPa,各段的平均累积库仑应力分别可达0.036 MPa,0.568 MPa和1.572 MPa,达到常用的触发阈值0.01 MPa。显著的强震库仑应力加载与地震地质研究给出的强震破裂空段(M7专项工作组,2012)、数值模拟给出的高应力积累区域(石富强等,2018)和断层高闭锁区域(Li et al,2017)吻合,表明祁连—海原断裂带西段的木里—江仓断裂和托莱山断裂以及中段的金强河—老虎山断裂未来的地震危险性高。考虑到青藏高原东北缘各活动构造之间的动力学相关性(Xiao,He,2015),这些区域未来的强震危险性不容忽视。
长安大学李振洪教授提供了2016年和2022年两次门源地震的破裂模型,德国地学研究中心(GFZ)汪荣江教授提供了库仑应力计算程序PSGRN/PSCMP,三位评审专家对本文提出了宝贵的意见和建议,作者在此一并表示感谢。
-
图 7 古浪地震位错模型对周边主要断裂的库仑应力影响
F1:祁连山北缘断裂;F2:香山—天景山断裂;F3:鄂拉山断裂;F4:日月山断裂;F5:海原断裂;F6:云雾山断裂;F7:龙首山断裂;F8:皇城—双塔断裂;F9:金强河断裂;F10:毛毛山断裂;F11:黄河—灵武断裂;F12:木里—江仓断裂;F13:拖莱山断裂;F14:冷龙岭断裂;F15:老虎山断裂;F16:六盘山断裂(a) 基于Xiao和He (2015)的位错模型;(b) 基于傅征祥等(2001)的位错模型;(c) 基于万永革等(2007)的位错模型;(d) 基于Guo等(2020)的位错模型
Figure 7. The cumulated Coulomb stress changes on the main active faults associated with different coseismic models of the 1927 Gulang earthquake
F1:Qilianshan northern marginal fault;F2:Xiangshan-Tianjingshan fault;F3:Elashan fault;F4:Riyueshan fault;F5:Haiyuan fault;F6:Yunwushan fault;F7:Longshoushan fault ;F8:Huangcheng-Shuangta fault;F9:Jinqianghe fault;F10:Maomaoshan fault;F11:Huanghe-Lingwu fault;F12:Muli-Jiangcang fault;F13:Tuolaishan fault;F14:Lenglongling fault;F15:Laohushan fault;F16:Liupanshan fault(a) The coseismic rupture model from Xiao and He (2015);(b) The coseismic rupture model from Fu et al (2001); (c) The coseismic rupture model from Wan et al (2007);(d) The coseismic rupture model from Guo et al (2020)
图 1 青藏高原北部构造环境和强震活动
沿祁连—海原断裂带(F3)的彩色散点为数值模拟给出的断层剪应力积累速率(石富强等,2018)。F1:阿尔金断裂带;F2:祁连山断裂带;F3:祁连—海原断裂带;F4:东昆仑断裂带;F5:柴达木盆地北缘断裂;F6:鄂拉山断裂;F7:日月山断裂;F8:狼山山前断裂;F9:六盘山断裂;F10:西秦岭北缘断裂;F11:甘孜—玉树断裂;F12:青川—平武断裂
Figure 1. The tectonic setting and the strong earthquake ruptures of the north Tibetan Plateau
The colored dots are the maximum shear stress rates of the Qilian-Haiyuan fault zone based on the finite element simulations (Shi et al,2018)。F1:Altyn fault zone;F2:Qilianshan fault zone;F3:Qilian-Haiyuan fault zone;F4:East Kunlun fault zone;F5:Northern marginal fault of Qaidam basin;F6:Elashan fault;F7:Riyueshan fault;F8:Langshan piedmont fault;F9:Liupanshan fault;F10:Northern marginal fault of west Qinling;F11:Garze-Yushu fault;F12:Qingchuan-Pingwu fault
图 2 祁连—海原断裂带三次显著强震之间的应力触发
(a) 1920年海原地震对1927年古浪地震的库仑应力加载;(b) 1920年海原地震对2022年门源地震的库仑应力加载;(c) 1927年古浪地震对2022年门源地震的库仑应力加载;(d) 1920年海原地震和1927年古浪地震对2022年门源地震的库仑应力加载
Figure 2. Stress interaction among the three strong earthquakes along the Qilian-Haiyuan fault zone
The magenta lines and beach balls are the current earthquake ruptures and the related focal mechanisms,and the light blue lines and beach balls express the receive faults and the related focal mechanisms. (a,b) The cumulated Coulomb stress changes associated with 1920 Haiyuan earthquake just before 1927 Gulang earthquake and 2022 Menyuan earthquake;(c) The cumulated Coulomb stress changes associated with the 1927 Gulang earthquake just before 2022 Menyuan earthquake;(d) The joint Coulomb stress interaction on 2022 Menyuan earthquake associated with 1920 Haiyuan and 1927 Gulang earthquakes
图 4 (a) 青藏高原北部历史强震对2022年门源M6.9地震的累积库仑应力加载;(b) 扣除1920年海原M8.5地震、1927年古浪M8.0地震和1954年山丹M7.3地震应力影响后,2022年门源M6.9地震断层面的库仑应力累积变化
Figure 4. (a) The temporal evolution of the Coulomb stress on the rupture plane of the 2022 Menyuan M6.9 earthquake associated with the strong earthquakes in Table 1;(b) Same as Fig. (a),but without the stress loading associated with the 1920 Haiyuan M8.5,1927 Gulang M8.0 and 1954 Shandan M7.3 earthquakes
图 8 祁连—海原断裂带库仑应力变化
(a) 1540年和1888年两次M7.0地震对祁连—海原断裂带的库仑应力影响;(b) 1540年M7.0地震、1888年M7.0地震和表1中的历史强震引起的祁连—海原断裂带的累积库仑应力变化
Figure 8. The Coulomb stress changes of the Qilian-Haiyuan fault zone
(a) The Coulomb stress changes associated with the 1540 M7 and 1888 M7 earthquakes;(b) The cumulated Coulomb stress changes of Qilian-Haiyuan fault zone associated with the 1540 M7 and 1888 M7 earthquakes as well as the strong earthquakes in Table 1
表 1 青藏高原北部及周边1900年以来的强震同震位错模型
Table 1 Coseismic rupture models of the strong earthquakes in northern Tibetan Plateau since 1900
发震日期
年−月−日地名 M 震中位置 走向
/°倾向
/°滑动角
/°破裂长度
/km破裂宽度
/kmSS
/mDS
/m来源 东经/° 北纬/° 1920−12−16 海原 8½ 104.10 37.04 110 88 14 23 30 4.39 −1.09 ① 104.50 36.90 112.5 88 14 64 30 6.31 −1.57 104.95 36.74 310.6 88 14 31 30 4.99 −1.24 105.58 36.52 112.6 88 14 39 30 6.21 −1.55 105.96 36.16 148.4 88 14 74 30 6.93 −1.73 1927−05−23 古浪 8.0 101.50 37.69 139 80 0 11 20 3 0 ② 101.58 37.61 122 80 0 11 20 3.3 0 101.69 37.56 116 80 0 12 20 3.6 0 101.81 37.51 122 80 0 11 20 4 0 101.91 37.46 102 80 0 16 20 6 0 102.09 37.43 99 80 0 16 20 7 0 102.27 37.41 92 80 0 14 20 7 0 102.42 37.41 95 80 0 18 20 3 0 102.63 37.39 90 80 0 14 20 2 0 102.25 38.07 139 40 90 42 20 0 −1.50 1932−12−25 昌马 7.6 96.70 39.70 115 79 30 116 20 2.34 −1.35 ① 1936−02−07 康乐 6.8 103.40 35.40 270 70 10 36 15 0.73 −0.13 ③ 1937−01−07 阿兰湖 7.5 97.60 35.50 110 70 15 208 20 3.96 −1.06 ④ 1947−03−17 达日 7.7 99.50 33.30 135 60 60 150 20 2.00 −3.46 ④ 1954−02−11 山丹 7.3 101.30 39.00 290 45 35 90 20 1.45 −1.02 ① 1954−07−31 腾格里 7.0 104.18 38.80 153 83 171 58 15 −0.94 −0.15 ① 1963−04−19 阿兰湖 7.0 97.00 35.70 277 80 −10 68 15 1.16 0.2 ① 1973−07−14 玛尼 7.0 86.48 35.18 81 60 −35 59 15 0.78 0.55 ① 1976−08−16 松潘 7.2 104.08 32.60 165 63 40 30 15 1.10 0.90 ④ 1976−08−23 平武 7.2 104.30 32.50 155 65 40 22 15 1.10 0.90 ④ 1990−04−26 共和 7.0 100.33 36.06 346 78 128 59 15 −0.41 −0.52 ① 1997−11−08 玛尼 7.5 87.33 35.07 76 90 −5 170 20 5 0.44 ⑤ 2001−11−14 昆仑山 8.1 90.54 35.95 99 90 5 346 20 4 −0.35 ⑤ 2008−05−12 汶川* 7.9 103.32 31.00 − − − − − − − ⑥ 2016−01−21 门源* 6.4 101.60 37.70 − − − − − − − ⑦ 2017−08−08 九寨沟* 7.0 103.82 33.20 − − − − − − − ⑧ 2021−05−21 玛多* 7.4 98.30 34.60 − − − − − − − ⑨ 2022−01−08 门源* 6.9 101.26 37.77 − − − − − − − ⑦ 注,*表示这些地震的相关震源参数采用有限断层模型反演给出;SS为走向方向滑动量,左旋为正;DS为倾向方向滑动量,正断为正。最后一列来源:① 万永革等(2007); ② Guo et al (2020) ;③ 梅秀苹等(2012);④Shan et al (2015) ;⑤ 沈正康(2003);⑥Shen et al (2009) ;⑦ 李振洪等(2022);⑧ 张旭等(2017);⑨ USGS (2021)。表 2 青藏高原东北缘岩石圈介质模型参数
Table 2 The model parameters of the lithosphere structure in the northeastern margin of the Tibetan Plateau
分层 厚度/km vP/(km·s−1) vS/(km·s−1) 密度
/(kg·m−3)ηk/(1018 Pa·s) ηm/(1019 Pa·s) 沉积层 10 5.9 3.41 2500 弹 性 上地壳 10 6.175 3.57 2700 低速层 12 5.85 3.38 2600 中地壳 10 6.4 3.70 3000 6.30 1.00 下地壳 20 6.8 3.93 3100 6.30 1.00 上地幔 − 8.1 4.68 3350 0 10.00 注:ηk为开尔文体黏滞系数,ηm为麦克斯韦尔体黏滞系数。 -
陈为涛,甘卫军,万永革,肖根如,孟令超. 2013. 青藏高原东北缘4个强震重点监视区库仑破裂应力的近百年变化和危险性分析[J]. 吉林大学学报(地球科学版),43(2):494–505. Chen W T,Gan W J,Wan Y G,Xiao G R,Meng L C. 2013. Coulomb failure stress and earthquake risk analysis of several four earthquake concerned zone on northeastern Tibetan Plateau in nearly hundred years[J]. Journal of Jilin University (Earth Science Edition)
,43(2):494–505 (in Chinese). 傅征祥,刘桂萍,陈棋福. 2001. 青藏高原北缘海原、古浪、昌马大地震间相互作用的动力学分析[J]. 地震地质,23(1):35–42. doi: 10.3969/j.issn.0253-4967.2001.01.004 Fu Z X,Liu G P,Chen Q F. 2001. Dynamic analysis on interaction between the Haiyuan-Gulang-Changma great earthquake in the north boundary of the Tibetan Plateau[J]. Seismology and Geology,23(1):35–42 (in Chinese).
国家地震局兰州地震研究所, 宁夏回族自治区地震局. 1980. 1920年海原大地震[M]. 北京: 地震出版社: 118. Lanzhou Institute of Seismology, State Seismological Bureau, Seismological Bureau of Ningxia Hui Autonomous Region. 1980. The 1920 Haiyuan Earthquake[M]. Beijing: Seismological Press: 118 (in Chinese).
韩竹军, 牛鹏飞, 李科长, 吕丽星. 2022. 2022年1月8日青海门源6.9级地震的一些初步认识[EB/OL]. [2022-01-18]. https://www.eq-igl.ac.cn/zhxw/info/2022/36632.html. Han Z J, Niu P F, Li K Z, Lü L X. 2022. Preliminary acknowledgements on the Jan. 8, 2022, Menyuan MS6.9 earthquake[EB/OL]. [2022-01-18]. https://www.eq-igl.ac.cn/zhxw/info/2022/36632.html (in Chinese).
李振洪, 韩炳权, 刘振江, 张苗苗, 余琛, 陈博, 刘海辉, 杜静, 张双成, 朱武, 张勤, 彭建兵. 2022. InSAR数据约束下的2016年和2022年青海门源地震震源参数及其滑动分布[J]. 武汉大学学报(信息科学版), doi: 10.13203/j.whugis20220037. Li Z H, Han B Q, Liu Z J, Zhang M M, Yu C, Chen B, Liu H H, Du J, Zhang S C, Zhu W, Zhang Q, Peng J B. Source parameters and slip distributions of the 2016 and 2022 Menyuan, Qinghai earthquakes constrained by InSAR observations[J]. Geomatics and Information Science of Wuhan University, doi: 10.13203/j.whugis20220037 (in Chinese).
M7专项工作组. 2012. 中国大陆大地震中-长期危险性研究[M]. 北京: 地震出版社: 196. Working Group of M7. 2012. Study on the Mid-to-Long Term Potential of Large Earthquakes on the Chinese Continent[M]. Beijing: Seismological Press: 196 (in Chinese).
梅秀苹,邵志刚,张浪平,冯建刚,代炜. 2012. 南北地震带北段强震破裂空段的地震危险性研究[J]. 地震学报,34(4):509–525. doi: 10.3969/j.issn.0253-3782.2012.04.008 Mei X P,Shao Z G,Zhang L P,Feng J G,Dai W. 2012. Study on potential earthquake risk of unbroken active faults in the northern segment of the North-South Seismic Zone[J]. Acta Seismologica Sinica,34(4):509–525 (in Chinese).
邵志刚,傅容珊,薛霆虓,黄建华. 2008. 昆仑山MS8.1级地震震后变形场数值模拟与成因机理探讨[J]. 地球物理学报,51(3):805–816. doi: 10.3321/j.issn:0001-5733.2008.03.021 Shao Z G,Fu R S,Xue T X,Huang J H. 2008. The numerical simulation and discussion on mechanism of postseismic deformation after Kunlun MS8.1 earthquake[J]. Chinese Journal of Geophysics,51(3):805–816 (in Chinese).
沈正康. 2003. 东昆仑活动断裂带大地震之间的黏弹性应力触发研究[J]. 地球物理学报,46(6):786–795. doi: 10.3321/j.issn:0001-5733.2003.06.010 Shen Z K. 2003. Viscoelastic triggering among large earthquakes along the east Kunlun fault system[J]. Chinese Journal of Geophysics,46(6):786–795 (in Chinese).
石富强,邵志刚,占伟,丁晓光,朱琳,李玉江. 2018. 青藏高原东北缘活动断裂剪切模量及应力状态数值模拟[J]. 地球物理学报,61(9):3651–3663. doi: 10.6038/cjg2018L0631 Shi F Q,Shao Z G,Zhan W,Ding X G,Zhu L,Li Y J. 2018. Numerical modeling of the shear modulus and stress state of active faults in the northeastern margin of the Tibetan Plateau[J]. Chinese Journal of Geophysics,61(9):3651–3663 (in Chinese).
石富强,张辉,邵志刚,徐晶,邵辉成,李玉江. 2020. 华北地区库仑应力演化与强震活动关系[J]. 地球物理学报,63(9):3338–3354. doi: 10.6038/cjg2020O0094 Shi F Q,Zhang H,Shao Z G,Xu J,Shao H C,Li Y J. 2020. Coulomb stress evolution and stress interaction among strong earthquakes in North China[J]. Chinese Journal of Geophysics,63(9):3338–3354 (in Chinese).
孙云强,罗纲. 2018. 青藏高原东北缘地震时空迁移的有限元数值模拟[J]. 地球物理学报,61(6):2246–2264. doi: 10.6038/cjg2018L0401 Sun Y Q,Luo G. 2018. Spatial-temporal migration of earthquakes in the northeastern Tibetan Plateau:Insights from a finite element model[J]. Chinese Journal of Geophysics,61(6):2246–2264 (in Chinese).
万永革,沈正康,曾跃华,盛书中. 2007. 青藏高原东北部的库仑应力积累演化对大地震发生的影响[J]. 地震学报,29(2):115–129. doi: 10.3321/j.issn:0253-3782.2007.02.001 Wan Y G,Shen Z K,Zeng Y H,Sheng S Z. 2007. Evolution of cumulative Coulomb failure stress in northeastern Qinghai-Xizang (Tibetan) Plateau and its effect on large earthquake occurrence[J]. Acta Seismologica Sinica,29(2):115–129 (in Chinese).
袁道阳. 2022. 科考简讯| 第二次青藏高原综合科学考察研究任务九专题: 对2022年1月8日青海门源6.9级地震应急科学考察取得重要阶段性成果[EB/OL]. [2022-01-14]. https://mp.weixin.qq.com/s/ArKFPnu_lby6tRcH69XLOQ. Yuan D Y. 2022. Special topic of the ninth research task of the secondary comprehensive scientific investigation on the Tibetan Plateau: Achievement of the emergency scientific investigation of the Jan. 8, 2022 Menyuan M6.9 earthquake[EB/OL]. [2022-01-14]. https://mp.weixin.qq.com/s/ArKFPnu_lby6tRcH69XLOQ (in Chinese).
张瑞,张竹琪,郑德文,刘兴旺,雷启云,邵延秀. 2021. 鄂尔多斯活动地块西缘强震间库仑应力作用[J]. 地球物理学报,64(10):3576–3599. doi: 10.6038/cjg2021P0008 Zhang R,Zhang Z Q,Zheng D W,Liu X W,Lei Q Y,Shao Y X. 2021. Coulomb stress transfer of strong earthquakes within tectonic belts near western Ordos block[J]. Chinese Journal of Geophysics,64(10):3576–3599 (in Chinese).
张旭,冯万鹏,许力生,李春来. 2017. 2017年九寨沟MS7.0级地震震源过程反演与烈度估计[J]. 地球物理学报,60(10):4105–4116. doi: 10.6038/cjg20171035 Zhang X,Feng W P,Xu L S,Li C L. 2017. The source-process inversion and the intensity estimation of the 2017 MS7.0 Jiuzhaigou earthquake[J]. Chinese Journal of Geophysics,60(10):4105–4116 (in Chinese).
中国地震台网中心. 2022. 1月8日1时45分在青海海北州门源县发生6.9级地震[EB/OL]. [2022-01-08]. https://www.cenc.ac.cn/cenc/dzxx/396391/index.html. China Earthquake Networks Center. 2022. An earthquake with MS6.9 occurred in Menyuan County, Haibei Prefecture, Qinghai Province, at 1:45 on January 8[EB/OL]. [2022-01-08]. https://www.cenc.ac.cn/cenc/dzxx/396391/index.html (in Chinese).
Ali S T,Freed A M,Calais E,Manaker D M,Mccann W R. 2008. Coulomb stress evolution in northeastern Caribbean over the past 250 years due to coseismic,postseismic and interseismic deformation[J]. Geophys J Int,174(3):904–918. doi: 10.1111/j.1365-246X.2008.03634.x
Broerse T,Riva R,Simons W,Govers R,Vermeersen B. 2015. Postseismic GRACE and GPS observations indicate a rheology contrast above and below the Sumatra slab[J]. J Geophys Res:Solid Earth,120(7):5343–5361. doi: 10.1002/2015JB011951
Guo P,Han Z,Gao F,Zhu C,Gai H. 2020. A new tectonic model for the 1927 M8.0 Gulang earthquake on the NE Tibetan Plateau[J]. Tectonics,39(9):e2020T–e6064T.
Han L,Liu Z J,Yao W,Shao Y,Yuan Z,Wang Y. 2021. Coseismic slip gradient at the western terminus of the 1920 Haiyuan MW7.9 earthquake[J]. J Struct Geol,152:104442. doi: 10.1016/j.jsg.2021.104442
Harris R A. 1998. Introduction to special section:Stress triggers,stress shadows,and implications for seismic hazard[J]. J Geophys Res,103(10):24347–24358.
He J,Lu S,Wang W. 2013. Three-dimensional mechanical modeling of the GPS velocity field around the northeastern Tibetan Plateau and surrounding regions[J]. Tectonophysics,584:257–266. doi: 10.1016/j.tecto.2012.03.025
Huang M,Bürgmann R,Freed A M. 2014. Probing the lithospheric rheology across the eastern margin of the Tibetan Plateau[J]. Earth Planet Sci Lett,396:88–96. doi: 10.1016/j.jpgl.2014.04.003
Jia K,Zhou S Y,Zhuang J C,Jiang C S,Guo Y C,Gao Z H,Gao S S. 2018. Did the 2008 MW7.9 Wenchuan earthquake trigger the occurrence of the 2017 MW6.5 Jiuzhaigou earthquake in Sichuan,China?[J]. J Geophys Res:Solid Earth,123(4):2965–2983. doi: 10.1002/2017JB015165
Jia K,Zhou S Y,Zhuang J C,Jiang C S. 2021. Stress transfer along the western boundary of the Bayan Har block on the Tibet Plateau from the 2008 to 2020 Yutian earthquake sequence in China[J]. Geophys Res Lett,48(15):e2021G–e94125G.
King G C,Stein R S,Lin J. 1994. Static stress changes and the triggering of earthquakes[J]. Bull Seismol Soc Am,84(3):935–953.
Kroll K A,Richards‐Dinger K B,Dieterich J H,Cochran E S. 2017. Delayed seismicity rate changes controlled by static stress transfer[J]. J Geophys Res:Solid Earth,122(10):7951–7965. doi: 10.1002/2017JB014227
Li Y C,Shan X J,Qu C Y,Zhang Y F,Song X G,Jang Y,Zhang G H,Nocquet J M,Gong W Y,Gan W J,Wang C S. 2017. Elastic block and strain modeling of GPS data around the Haiyuan-Liupanshan fault,northeastern Tibetan Plateau[J]. J Asian Earth Sci,150:87–97. doi: 10.1016/j.jseaes.2017.10.010
Parsons T,Dreger D S. 2000. Static‐stress impact of the 1992 Landers earthquake sequence on nucleation and slip at the site of the 1999 M=7.1 Hector Mine earthquake,southern California[J]. Geophys Res Lett,27(13):1949–1952. doi: 10.1029/1999GL011272
Pollitz F F, Banerjee P, Bürgmann R, Hashimoto M, Choosakul N. 2006. Stress changes along the Sunda trench following the 26 December 2004 Sumatra‐Andaman and 28 March 2005 Nias earthquakes[J]. Geophys Res Lett, 33(6): L06309. doi: 10.1029/2005GL024558.
Shan B,Xiong X,Wang R J,Zheng Y,Yadav R B S. 2015. Stress evolution and seismic hazard on the Maqin-Maqu segment of east Kunlun fault zone from co-,post- and interseismic stress changes[J]. Geophys J Int,200(1):244–253. doi: 10.1093/gji/ggu395
Shen Z K,Sun J B,Zhang P Z,Wan Y G,Wang M,Bürgmann R,Zeng Y H,Gan W J,Liao H,Wang Q L. 2009. Slip maxima at fault junctions and rupturing of barriers during the 2008 Wenchuan earthquake[J]. Nat Geosci,2(10):718–724. doi: 10.1038/ngeo636
Steacy S, Marsan D, Nalbant S S, Mccloskey J. 2004. Sensitivity of static stress calculations to the earthquake slip distribution[J]. J Geophys Res: Solid Earth, 109(B4): B04303.
Stein R S,Barka A A,Dieterich J H. 1997. Progressive failure on the North Anatolian fault since 1939 by earthquake stress triggering[J]. Geophys J Int,128(3):594–604. doi: 10.1111/j.1365-246X.1997.tb05321.x
Stein R S. 1999. The role of stress transfer in earthquake occurrence[J]. Nature,402(6762):605–609. doi: 10.1038/45144
Toda S, Lin J, Meghraoui M, Stein R S. 2008. 12 May 2008 M=7.9 Wenchuan, China, earthquake calculated to increase failure stress and seismicity rate on three major fault systems[J]. Geophys Res Lett, 35(17), doi: 10.1029/2008GL034903.
USGS. 2021. M7.3: Southern Qinghai, China[EB/OL]. [2021-05-21]. https://earthquake.usgs.gov/earthquakes/eventpage/us7000e54r/executive.
Verdecchia A,Pace B,Visini F,Scotti O,Peruzza L,Benedetti L. 2018. The role of viscoelastic stress transfer in long‐term earthquake cascades:Insights after the central Italy 2016–2017 seismic sequence[J]. Tectonics,37(10):3411–3428. doi: 10.1029/2018TC005110
Wan Y G,Shen Z K. 2010. Static Coulomb stress changes on faults caused by the 2008 MW7.9 Wenchuan,China earthquake[J]. Tectonophysics,491(1/2/3/4):105–118. doi: 10.1016/j.tecto.2010.03.017
Wang J J,Xu C J,Freymueller J T,Li Z H,Shen W B. 2014. Sensitivity of Coulomb stress change to the parameters of the Coulomb failure model:A case study using the 2008 MW7.9 Wenchuan earthquake[J]. J Geophys Res:Solid Earth,119(4):3371–3392. doi: 10.1002/2012JB009860
Wang R,Lorenzo-Martín F,Roth F. 2006. PSGRN/PSCMP:A new code for calculating co- and post-seismic deformation,geoid and gravity changes based on the viscoelastic-gravitational dislocation theory[J]. Comput Geosci,32(4):527–541. doi: 10.1016/j.cageo.2005.08.006
Wells D L,Coppersmith K J. 1994. New empirical relationships among magnitude,rupture length,rupture width,rupture area,and surface displacement[J]. Bull Seismol Soc Am,84(4):974–1002.
Xiao J,He J. 2015. 3D finite-element modeling of earthquake interaction and stress accumulation on main active faults around the northeastern Tibetan Plateau edge in the past ~100 years[J]. Bull Seismol Soc Am,105(5):2724–2735. doi: 10.1785/0120140342
Yang H F, Wang D, Guo R M, Xie M Y, Zang Y, Wang Y, Yao Q, Cheng C, An Y R, Zhang Y Y. 2022. Rapid report of the 8 January 2022 MS6.9 Menyuan earthquake, Qinghai, China[J]. Earthquake Research Advances, 100113.
-
期刊类型引用(15)
1. 夏新宇,高原. 青藏高原东北缘构造特征与分层各向异性. 地震学报. 2025(01): 1-20 . 本站查看
2. 姚生海,盖海龙,殷翔,李鑫,刘炜. 青海门源M6.9级地震地表破裂特征及区域地震活动趋势分析. 地质通报. 2024(Z1): 340-349 . 百度学术
3. 蒲颂文,闻鑫,周志伟,汪汉胜,江利明,李晨程. 基于Sentinel-1的门源Ms 6.9三维同震形变提取与分析. 地理空间信息. 2024(06): 109-112+117 . 百度学术
4. 徐伟,程理,李光涛,陈文龙. 马雅雪山断裂金强河段晚第四纪断错地貌及滑动速率确定. 震灾防御技术. 2024(03): 457-467 . 百度学术
5. 曾锐,蒋亚楠,燕翱翔,程燕,罗袆沅. InSAR约束的青海门源M_w6.7级地震同震断层滑动反演与三维形变场模拟. 地质科技通报. 2024(06): 212-225 . 百度学术
6. 程燕,蒋亚楠,侯中健,曾锐,罗袆沅. 2016年和2022年青海门源强震活动的InSAR形变观测与区域强震危险性分析. 地质力学学报. 2024(06): 965-977 . 百度学术
7. 袁道阳,谢虹,苏瑞欢,李智敏,文亚猛,司国军,薛善余,陈干,刘炳旭,梁淑敏,彭慧,段磊,魏拾其. 2022年1月8日青海门源M_S6.9地震地表破裂带特征与发震机制. 地球物理学报. 2023(01): 229-244 . 百度学术
8. 朱俊文,蔡艺萌,张波,姚赟胜. 昌马M_S7.6地震三维同震形变场数值模拟及库仑应力触发研究. 地震工程学报. 2023(03): 598-605 . 百度学术
9. 龚炜程,孙云强,邱鑫鹏,王甘娇,李敏娟,阮云凯. 青藏高原东北缘主要断裂带的库仑应力演化及其强震危险性. 地震工程学报. 2023(03): 585-597 . 百度学术
10. 汤大委,葛伟鹏,袁道阳,郑文俊,胡才博. 青藏高原北部历史强震对2022年门源M_S6.9地震及后续地震库仑应力触发作用. 地球物理学报. 2023(07): 2772-2795 . 百度学术
11. Sixin Zhang,Peng Jia,Bowen Hou,Ming Hao. Vertical deformation before and after the 2022 Menyuan Ms6.9 earthquake and analysis of earthquake precursors. Geodesy and Geodynamics. 2023(06): 582-588 . 必应学术
12. 李媛,杨周胜,庞亚瑾,梁洪宝,刘峡. 2022年门源M_S6.9地震前断层活动及应力状态的数值模拟. 地震地质. 2023(06): 1286-1308 . 百度学术
13. 韩立波. 2022年青海门源M_S6.9地震震源机制解. 地震科学进展. 2022(02): 49-54 . 百度学术
14. 刘璐,刘兴旺,张波,何文贵,朱俊文,蔡艺萌. 利用高分影像识别门源M_S6.9地震地表破裂带. 地震工程学报. 2022(02): 440-449 . 百度学术
15. 苑争一,赵静,牛安福. 2022年1月8日青海门源M_S6.9地震震前变形特征分析. 中国地震. 2022(03): 389-398 . 百度学术
其他类型引用(6)