地球物理学报  2021, Vol. 64 Issue (12): 4548-4561   PDF    
基于InSAR地表形变约束的玛多MS7.4地震孕育发生机理数值模拟研究
祝爱玉1,2,3, 王永哲1,2,3, 李永华1,2,3, 张东宁1,2     
1. 中国地震局地球物理研究所, 北京 100081;
2. 北京白家疃地球科学国家野外科学观测研究站, 北京 100095;
3. 中国地震局震源物理重点实验室, 北京 100095
摘要:2021年5月22日在青海玛多发生MS7.4地震,此次地震造成了严重的地表破坏.本文通过合成孔径雷达干涉测量(InSAR)技术得到了Sentinel-1升、降轨同震地表形变、东西方向和垂直方向的同震形变场.并根据东西向地表的错动,勾画出了地震破裂面分段结构.以InSAR地表形变观测为约束,构建了能反映玛多地震破裂特点的三维黏弹性有限元模型,模拟得到了玛多地震的三维同震形变及应力场变化和演化特征,讨论了玛多地震的成因及玛多地震的发生对周边断层的影响.结果显示:玛多地震为典型的左旋走滑型地震,其中东部错动量相对较大;玛多地震的余震都分布在库仑应力升高的影区,说明玛多地震序列使得断层上积累的应变得到了充分的释放,未来余震的发生会朝西南和东北方向发展.玛多地震的发生形成了多个应力强加载区,推测玛多地震的发生大大增加了部分临近断裂带滑动的可能.2010年玉树7.3级地震引起的库仑应力变化在玛多地震主破裂面上的投影为正数,推测玉树地震对玛多地震具有促进作用.研究结果可以为大震后的余震分布判断提供定量参考依据.
关键词: 玛多地震      同震位移      库仑应力变化      InSAR形变观测      数值模拟     
Numerical simulation on the mechanism of the Madoi, Qinghai MS7.4 earthquake constrained by InSAR deformation
ZHU AiYu1,2,3, WANG YongZhe1,2,3, LI YongHua1,2,3, ZHANG DongNing1,2     
1. Institute of Geophysics, China Earthquake Administration, Beijing 100081, China;
2. Beijing Baijiatuan Earth Science National Observation and Research Station, Beijing 100095, China;
3. KeyLaboratory of Earthquake Source Physics, China Earthquake Administration, Beijing 100095, China
Abstract: The Madoi MS7.4 earthquake occurred in Qinghai Province on 22 May 2021, which caused serious damage on the surface. In this work, using the interferometric synthetic aperture radar (InSAR) technology, coseismic surface deformation in ascending and descending directions of Sentinel-1, east-west and vertical directions were obtained, and the segmented structure of earthquake rupture surface was delineated. Constrained by the observed surface deformation from InSAR, a 3D viscoelastic finite element model considering the rupture characteristics of the Madoi earthquake was constructed. Then the 3D coseismic deformation and stress field of the Madoi earthquake were simulated. Results show that the Madoi earthquake is a typical left-lateral strike-slip faulting, in which the fault slip in the east is relatively large. The aftershocks of the Madoi earthquake are all distributed in the area where the Coulomb stress variation is negative, which indicates that the Madoi earthquake has fully released the accumulated strain on the fault and the occurrence of aftershocks in the future will spread towards the southwest and northeast. The occurrence of the Madoi earthquake has formed several stress enhancement areas, which indicates that the Madoi event has greatly increased the possibility of sliding of some other fault zones. The projection of Coulomb stress variation caused by the 2010 Yushu M7.3 earthquake on the main rupture surface of the Madoi earthquake is positive, which permits to speculate that the Yushu earthquake may have promoted the Madoi shock. This study can provide a quantitative reference for the judgment of aftershocks and possible surface damage following big earthquakes.
Keywords: Madoi earthquake    Coseismic displacement    Coulomb stress variation    InSAR deformation observation    Numerical simulation    
0 引言

根据中国地震台网中心测定, 2021年5月22日2时4分在青海省果洛州玛多县发生7.4级地震, 震中位于北纬34.59°, 东经98.34°, 震源深度17 km. 震后的地质科考结果表明, 玛多MS7.4地震发震构造西段发现大型地表破裂带.该破裂带总体走向NW230°—270°, 由两条NWW向的剪切破裂左阶羽列组成.破裂带由一系列的挤压鼓包和张裂隙相间排列, 张裂隙从十几厘米到2~3 m不等, 裂缝通常呈右阶雁列状排列, 并显示明显垂直位错, 垂直位错可达约30 cm, 并可清晰看到水平擦痕, 显示为具有左旋走滑兼正断位移的同震断层.

大地测量是获取地震同震地表形变的重要手段, 并为地震发生机理的研究提供重要基础数据(袁霜等, 2020; 方进等, 2019; 申文豪等, 2019; 陈威等, 2018; 季灵运等, 2017; Xu, 2017; Song et al., 2019; Jónsson et al., 2002; Funning et al., 2005). 同震地表形变观测中, 以全球定位系统(GPS)技术和合成孔径雷达干涉测量(InSAR)技术应用最广泛.GPS观测虽然得到的形变数据精度较高, 但其空间分辨率取决于观测站的分布和数量, 玛多地震所在区域的自然环境恶劣, 用于进行地表形变观测的GPS观测站非常少(朱亚戈等, 2021), 难以准确反映本次地震的地表详细形变特征.InSAR技术由于可以获得米级乃至分米级空间分辨率且覆盖几十公里范围内的形变信息, 使得该技术在1993年首次应用于强震地表形变观测(Massonnet et al., 1993)以来得到广泛应用(Wang et al., 2019; Wen et al., 2021; Yang et al., 2021; Xu et al., 2016).该技术还具有全天时、全天侯、受天气影响小等特点, 可以在地表植被覆盖稀疏的环境条件下保持高相干性, 获取雷达视线方向1~2 cm精度的形变(孙建宝等, 2007), 尽管观测精度不如GPS技术, 但对于强震引起的地表形变观测已经足够.此次地震区域地表植被稀疏, 非常适合利用InSAR技术获取此次地震的同震地表形变场.采用InSAR观测玛多地震的同震形变场, 不仅能深入的研究地震造成的地表影响, 同时能为数值模拟提供有力的约束.

玛多地震是2008年汶川地震之后中国发生的震级最高的一次强震, 该地震发生在巴颜喀拉地块内部的昆仑山口—江错断裂上(王未来等, 2021).其北部和南部的东昆仑断裂带和甘孜—玉树断裂带, 为青藏高原东北缘的巨型走滑断裂带(邓起东等, 2002; 刘光勋, 1996; Fu and Awata, 2005; Kirby et al., 2007; Karplus et al., 2011; 李陈侠等, 2011; 张培震等, 2008).根据历史地震记录, 1900年以来, 玛多地震震中周边300 km范围内共发生7级以上地震5次, 依次为1937年托索湖7.5级地震(李陈侠等, 2006)、1947年达日7.7级地震(刘雷等, 2021; 梁明剑等, 2020)、1963年阿兰湖7.0级地震、1990年共和7.0级地震(李旭和陈运泰, 1996)和2010年玉树7.3级地震(张勇等, 2010; 孙鑫喆, 2012).强震导致的库仑应力变化和后续地震的发生在时间、空间上具有一定的相关性, 断层带上库仑应力的增加对后续中强地震的发生具有显著促进作用(Toda et al., 2011a; Xiong et al., 2010; 石耀霖和曹建玲, 2010; Harris, 1998).为此开展历史强震对玛多地震的影响研究, 对于分析玛多地震的孕育机理具有重要的意义.

地震精定位时空序列显示玛多地震呈现出典型的主震-余震型特征, 余震序列总长度约170 km.主震位于余震区中间, 呈双侧破裂特征(王未来等, 2021).强震的同震位错会导致邻近区域应力状态的改变(Toda et al., 2008; 单斌等, 2009), 库仑应力变化大于0.01 MPa就会触发地震(Toda et al., 2011b; Ma et al., 2005), 为此, 通过研究玛多地震序列导致的应力场及库仑应力场变化, 可为确定本次地震对其周边断裂带的影响, 进而预判余震趋势具有重要的意义(Li et al., 2021; Liu et al., 2014; Stein and Okal, 2011; Luo and Liu, 2010; 严珍珍等, 2009; 朱守彪等, 2008; Rundle et al., 2006)(图 1).

图 1 玛多地震周边构造 Fig. 1 Geological structure around Madoi earthquake

本文采用合成孔径雷达干涉测量(InSAR)技术得到地表形变场及地震破裂面分段结构, 并建立能反映玛多地震分段破裂面的黏弹性有限元模型, 以InSAR观测的地表形变、该地区深部结构、岩石圈流变特性、GPS观测给出的块体运动为约束, 模拟计算了玛多地震引起的三维同震位移及同震应力变化, 讨论了玛多地震的发生对区域应力场及周边断层的影响; 以玛多破裂面为接收断层模拟了周边5次强震的同震库仑应力变化及震后松弛, 讨论了周边强震对玛多地震的触发作用.

1 InSAR同震地表形变获取

地震发生后, 为全面了解地震的地表形变特征, 我们收集了欧空局Sentinel-1A/1B两颗卫星升、降轨的合成孔径雷达影像, 组成两个干涉对.这两幅干涉对的时间跨度均为6 d, 即震前观测时间为2021年5月20日和震后2021年5月26日.美国Scripps海洋研究所Sandwell教授团队开发的开源GMTSAR软件(Sandwell et al., 2011)作为本文干涉处理的软件, 采用的处理策略为二轨差分方法, 利用水平分辨率3弧秒的SRTM地表高程数据(Farr et al., 2007)作为差分基础, 从而去除原始干涉图中的地形相位贡献.为保证获取的地表形变的可靠性, 在使用SNAPHU算法进行差分相位解缠时, 我们只对相干性高于0.1的区域进行了解缠, 而忽略了低于0.1的区域形变信息.最终, 将解缠后的雷达坐标系下的结果地理编码到WGS-84坐标系下, 得到此次地震的升、降轨雷达视线方向同震形变场.结果显示, 此次地震引起的升轨方向最大抬升量和沉降量分别为: 1.398 m和1.482 m, 而降轨则为: 1.291 m和1.506 m.为提高建模时的计算效率, 我们对升、降轨同震形变场进行了重采样.采样的方法是每3 km取一个观测数据点, 同时, 不考虑相干性小于0.3的观测数据点.用于有限元模型的地表形变约束数据, 如图 2所示.

图 2 玛多地震地表InSAR同震形变观测 (a) 升轨方向; (b) 降轨方向. Fig. 2 The observed surface deformation of Madoi earthquake by InSAR (a) Ascending direction; (b) Descending direction.

为确定断层所在位置以及对地震发生机制进行定性分析(杨亚夫等, 2016), 本文利用升、降轨雷达视线方向形变测量结果, 根据SAR图像的成像几何获得了东西方向和垂直方向的地表形变结果.具体的, 利用式(1)所表征的地表三维形变与InSAR雷达视线向形变的关系, 考虑南北方向形变对于雷达视线方向贡献最小, InSAR技术对南北向形变的识别也最不敏感(王永哲等, 2012).因此, 在忽略雷达视线方向中的南北方向形变分量情况下, 利用升、降轨视线方向形变解算出东西方向和垂直方向的同震形变场(图 3).

图 3 玛多地震的地表东西向(a)和垂直向(b)的InSAR同震形变场 Fig. 3 The coseismic surface deformation in the (a) East-West and (b) Vertical directions of the Madoi earthquake observed by InSAR

从图中可以看出, 地表的向西运动和向东运动存在明显的分界线, 可推断此次地震破裂出地表, 形成了较长的破裂迹线; 破裂迹线两侧存在明显的东西反向运动, 这与左旋走滑的机制非常吻合.另外, 垂直向运动相对于东西向运动来讲量级较小, 可以推断此次地震伴有少量的两侧相对升降运动.根据东西向地表的错动, 勾画了断层在地表的破裂迹线, 将整个断层共分成6段用于数值模拟的分析.考虑到东西向和垂直向形变的误差, 本文在建模时, 未将东西向和垂直向形变作为约束进行拟合.

(1)

其中, dlos为雷达视线方向形变, dedu分别为东西向和垂直向形变, θα分别为雷达入射角和飞行方向角.

2 有限元数值模拟 2.1 数值模型

InSAR形变观测给出了玛多地震的破裂面分为6段(图 3), 本文在建立几何模型时充分考虑玛多地震破裂面分段, 建立了能反映分段破裂面及周边断层(邓起东等, 2002)的三维黏弹性有限元模型, 模型范围为95°E—102.5°E, 31.5°N—37°N, 深度为100 km(图 4).研究区包含2条大型的走滑断裂带, 分别为东昆仑断裂带及甘孜—玉树断裂带,玛多地震的发震构造即昆仑山口—江错断裂, 以及一些较小的构造, 例如: 玛多—甘德断裂、达日断裂、鄂拉山断裂、日月山断裂、共和盆地的南缘和北缘断裂等.模型为XYZ坐标系, X轴朝东、Y轴朝北、Z轴朝上.

图 4 三维有限元模型 (a) 几何模型; (b) 弹性模量; (c) 黏滞系数. Fig. 4 Three dimensional finite element model (a) Geometric model; (b) Elastic modulus; (c) Viscosity coefficient.

采用最新的联合反演获得的青藏高原东北缘及邻区的高分辨率P、S波速及重力反演得到的绝对密度结果(夏思茹等, 2021; Shen et al., 2016; Stolk et al., 2013), 其中P、S波速及密度数据的分辨率约为0.5°×0.5°.计算泊松比ν=和杨氏模量E=vp2ρ(1+ν)(1-2ν)/(1-v), 其中vpvs分别为P波、S波波速, ρ为密度.黏滞系数采用孙玉军等(2013)给出的中国大陆及邻区岩石圈三维流变结构计算结果, 其分辨率约为20 km×20 km. 本文采用插值方法将弹性模型、泊松比及黏滞系数等参数赋值给模型中的每个单元如图 4b, 4c所示.由于断层的强度一般比周边地块低, 本文将断层的弹性模量和黏滞系数分别设置为10 GPa、10-19Pa·s.和以往采用的分块模型相比, 本模型能充分的反映黏弹性本构模型的材料非均匀性.网格尺寸在断裂附近控制在3 km内, 其他区域控制在10 km内.

Wang和Shen(2020)给出的中国大陆及其邻区地壳GPS速度场(相对于欧亚大陆)作为模型边界条件和检验模拟结果的依据.采用GPS的观测速度值插值到模型边界给出模型四周边界的水平速度值.该速度值乘以计算时间步长作为有限元模型的水平位移约束条件, 且假定模型四周的水平约束不随深度变化, 即假设从地表到100 km深度保持一致, 垂直方向位移可以保持自由.三维模型底面所有节点垂直方向(z方向)上位移约束为0, 水平向可以自由运动.模型上表面为自由表面, 即法向应力和剪应力均为零.

采用多物理场耦合有限元软件COMSOL进行模拟计算, 时间步长500年(依据模型给定弹性和流变材料参数计算出的数百年至数千年的黏弹性松弛时间), 根据时间步长逐步施加边界位移约束.经过了40个时间步, 即20000年的加载作用, 以最后得到的研究区稳定地壳构造应力场作为研究地区的背景构造应力场.通过模拟计算的地表形变和现今GPS观测形变对比分析(图 5), 可以看出计算值和测量值从方向和大小来看都比较相近, 说明了计算模型的可靠性.

图 5 地壳形变速度的计算值和测量值的对比(图中红色箭头为模拟结果; 蓝色箭头为Wang和Shen(2020)发表的相对欧亚板块的GPS观测速度结果) Fig. 5 Comparison of calculated and observed horizontal deformation rate (The red arrows are the simulated results; The blue arrows are GPS observations (Wang and Shen, 2020))
2.2 同震位移模拟结果

基于黏弹性三维有限元模型, 构建符合实际的背景应力场.通过降低玛多地震分段破裂区域的弹性模量, 采用试错法确定破裂段1, 2, 3, 4的弹性模量降低为7.8 GPa, 破裂段5, 6的弹性模量降低为5.8 GPa, 以地表的InSAR形变观测为约束, 从而确定玛多地震引起的位移场和应力场的变化.由于本文是将玛多地震的主震和震后大约8 d的时间内余震作为一个整体进行求解, 因此本文的应力场变化能反映主震和大部分余震形成的应力释放.

为了将模拟的地表形变与Sentinel-1升、降轨InSAR同震地表形变比较, 采用公式(1)将三维地表形变投影到雷达视线方向(升、降轨)上, 并与实际观测的形变(图 2a, 6b)相减, 得到观测值与模拟值之间的残差如图 6a(升轨)和6c(降轨)所示.根据残差数据分布的直方图(图 6b, 6d)可以看出升轨的残差超过80%的值小于0.1 m, 95%的残差不超过0.2 m; 降轨的残差分布85%的值不超过0.2 m. 在进行模型构建时, 本文将整个断裂简单地分成了6段.然而, 此次玛多地震实际的地表破裂比较复杂, 并非简单的直线分段, 这会造成本来位于上盘的观测点在模拟时按照下盘来模拟, 下盘的观测点按照上盘来模拟, 这必然会造成模拟值与观测值存在较大的残差.从残差图中可以看出, 尽管有少数几个观测点的残差值较大, 但这些残差较大的点对于此次地震的整个形变场模式影响不大, 模拟的升、降轨形变场与实际观测值具有较好地一致性.根据以上观测值和模拟值的对比分析, 可以看出模拟得到的地表形变分布特征和InSAR观测的地表形变基本吻合.

图 6 模拟值和InSAR观测的地表形变残差 (a) 升轨方向观测值与模拟值的残差; (b) 升轨方向残差直方图; (c) 降轨方向观测值与模拟值的残差; (d) 降轨方向残差直方图. Fig. 6 Surface deformation residuals of simulated results and InSAR observations (a) The residuals of the observed and simulated values in the ascending direction; (b) The histogram of the residual error in the ascending direction; (c) The residual error between the observed value and the simulated value in the descending direction; (d) Histogram of residual error in the descending direction.

本文模拟得出的三维同震位移(图 7)结果显示, 东西向分量最大, 其最大和最小值分别为1.1 m和-1.38 m; 其次是南北向分量, 其最大和最小值分别为0.7 m和-0.74 m; 垂向分量相对较小, 最大为0.3 m.玛多地震总的同震位移(图 7d)呈现出“四象限”分布, 是走滑断层错动的典型图像.同震位移结果还显示, 此次地震发震断层的北侧, 位移方向为北西向; 在发震断层的南侧, 位移方向为东南向; 发震断层的南侧位移要略大于北侧位移, 比较大的错动主要集中在东侧和西侧, 中间部位的错动相对较小.

图 7 三维同震位移模拟结果 (a) 东西向分量; (b) 南北向分量; (c) 垂直向分量; (d) 总位移. Fig. 7 The simulation results of 3D coseismic displacement (a) East-West component; (b) North-South component; (c) Vertical component; (d) Total displacement.
2.3 同震应力场变化模拟结果

以地表同震位移为约束, 本文给出了10 km深度的同震应力变化和库仑应力变化(图 8).库仑应力变化ΔCFS为

(2)

图 8 同震应力变化 (a) 正应力变化; (b) 剪应力变化; (c) 库仑应力变化; (d) 库仑应力变化和震后3个小时的余震分布; (e) 库仑应力变化和震后8天内的余震分布. Fig. 8 Coseismic stress variation (a) Normal stress variation; (b) Shear stress variation; (c) Coulomb stress variation; (d) Coulomb stress variation and aftershock distribution in 3 hours after the Madoi earthquake; (e) Coulomb stress variation and aftershock distribution within 8 days after the Madoi earthquake.

其中Δτ为断层面上的剪切力变化; Δσn为断层面上的正应力变化, 压应力为负数; μ为断层面上的有效摩擦系数, 一般取值范围为0.2~0.8, 本文取值为0.4(Stein et al., 1992; King et al., 1994; Freed, 2005; 石耀霖和曹建玲, 2010).有效摩擦系数的取值会在一定程度上改变库仑应力的大小, 但不会影响库仑应力正负的改变.

图 8a显示在整个破裂面上是以压应力为主, 在破裂面的西南向和北东向有拉张分量, 并呈现不规则的四象限分布.图 8b显示破裂面的北侧和南侧剪应力增强, 而东侧和西侧剪应力减低, 因此地震呈现左旋走滑趋势.以上正应力和剪应力的分布特征都揭示了玛多地震为典型的走滑型地震的特征.由于本文考虑了区域弹性模量等参数的非均匀性, 以及区域断层分布的影响, 而应力场分布受区域构造和深部结构的影响, 因此本文结果更能反映区域应力场分布特征.

根据玛多地震的震源机制解(郭祥云, http://www.cea-igp.ac.cn/kydt/278249.Html; 万永革等, https://mp.weixin.qq.com/s/60hFBoaBXa-c_Q9BVCNYHw), 其破裂面的整体走向为101°, 倾向为87°, 滑动角为-7°, 采用公式(2)计算得到10 km处的同震库仑应力变化(图 8c, 8d, 8e).图 8d给出了震后3个小时的地震精定位结果, 可以看出余震主要分布在西部和东部, 而中间部位有一段地震空区.图 8e给出了震后8天的余震精定位, 随着时间的推移, 余震数量越来越多, 可以看出余震向西部和东部继续扩展, 东部的余震逐渐转向为北东向, 中间部位的地震空区也在逐渐变小.对比库仑应力变化可以看出, 余震都分布在库仑应力影区, 说明玛多地震已经充分释放了应变,其和库仑应力影区有很好的一致性, 同时中间的地震空区和库仑应力增强区域也有很好的对应关系.

从库仑应力变化(图 8c)看出, 玛多地震发震断层两侧的震后库仑应力变化为负数.说明玛多地震的发生, 使得发震断层南北两侧区域应力都得到了充分的释放.而库仑应力的增加区域是发震断层的西南和东北端, 推测未来的余震的发生主要向西南和东北方向发展.

3 讨论 3.1 区域构造应力对玛多地震的影响

根据三维应力和应变模拟计算结果给出了未考虑历史强烈地震活动时研究区最大剪应力分布模拟结果(图 9), 其中最大剪应力为(σ1-σ3)/2, σ1为最大主压应力, σ3为最小主压应.图 9中的黑色五角星表示自1900年以来发生在研究区域M7级以上地震(表 1).首先从不同深度的最大剪应力分布可以看出, 10 km处的剪应力最大.根据地震精定位的研究结果(王未来等, 2021), 玛多地震的主震深度为10 km左右, 这与我们的计算结果一致.其次, 区域剪应力高值区的分布, 主要集中在东昆仑断裂带西部、昆仑山口—江错断裂、甘孜—玉树断裂的西部, 这些剪应力高值区和历史的强震有很好的对应关系, 尤其是1937年托索湖地震、1963年阿兰湖地震、2010年玉树地震和2021年玛多地震等.同时, 高剪应力区域和弹性模量低值区域相对应(图 4b), 即在周边块体的挤压作用下, 较低的弹性模量更容易累积应变, 因而使得断层更加容易滑动.此次玛多7.4级地震明显处于应力增强区域, 因此地震的发生和周边块体的挤压作用及地处低弹性模量区域密不可分.

图 9 不同深度的最大剪应力 (a) 5 km处; (b) 10 km处; (c) 15 km处; (d) 20 km处. Fig. 9 Maximum shear stress at different depths (a) 5 km; (b) 10 km; (c) 15 km; (d) 20 km.
表 1 区域1900年以来研究区7级以上地震目录和破裂方向数据 Table 1 Catalogue and fracture direction of earthquakes with magnitude 7 and above in the study area since 1900
3.2 周边历史强震对2021年玛多MS7.4地震的触发作用

本文根据历史地震目录统计了玛多地震震中周边300 km范围内, 1900年以来发生7级以上地震5次(表 1), 其中1990年共和地震是逆冲型地震, 其他4次强震发生在巴颜喀拉地块的边界或者内部, 都是走滑型地震, 而1937年的托索湖地震破裂尺度达到240 km, 影响范围极广.

本文计算了以上5次7级以上地震事件投影到玛多地震破裂面的同震库仑应力变化(图 10).根据同震库仑应力的影响可以看出(图 10a), 1937年托索湖7.5级地震、1963年阿兰湖7.0级地震和1990年共和7.0级地震在玛多地震主破裂面上的投影都为负数, 说明其对玛多地震破裂面有卸载的作用, 使得玛多地震的发生推迟.2010年玉树7.3级地震的发生在玛多地震主破裂面的库仑应力投影为正数, 说明玉树地震对玛多地震具有促进作用.

图 10 1900年以来区域5次强震投影在玛多地震破裂面上的库仑应力变化 (a) 同震库仑应力变化; (b) 考虑构造应力的库仑应力变化. Fig. 10 Coulomb stress variation of 5 strong earthquakes projected on the fracture surface of Madoi earthquake (a) Coseismic Coulomb stress variation; (b) Coulomb stress variation considering tectonic stress.

其次, 我们考虑了震后的黏弹松弛效应, 计算了以上5次强震在玛多地震主破裂面的库仑应力变化(如图 10b), 可以看出随着震后的松弛效应, 玛多主破裂面的库仑应力变化在块体边界的作用力下, 应变在不断的累积, 随着区域的地震的发生, 例如1963年的阿兰湖地震、1990年的共和地震使得玛多地震发震断层的应力得到了部分释放, 从而推迟了玛多地震的发生, 而2010年的玉树地震使得玛多地震发震断层的应力突然增加, 从而加速了玛多地震的发生.

3.3 玛多地震序列对周边断层的影响

由于以上的同震库仑应力的计算是以玛多地震发震断层为接收断层, 因此只能反映和发震断层相同走向和倾向的周边断层的库仑应力影响, 而不能反映其他周边断层的影响.为了进一步研究玛多地震对周边主要断层的库仑应力影响, 本文给出了以周边主要断层为接收断层的库仑应力变化(图 11所示), 可以看出玛多地震序列的发生, 使得南北向的应力得到了充分的释放, 导致东昆仑断裂带的西部、玛多—甘得断裂带的中部偏西部位、达日断裂带的西部以及甘孜—玉树断裂带的东部区域的库仑应力变化为负; 库仑应力变化明显增强的区域为昆仑山口—江错断裂的北部、玛多—甘德断裂的西部及东部部分区域、东昆仑断裂带中部、共和盆地南缘断裂和北缘断裂, 以上区域的库仑应力变化明显大于0.01 MPa, 说明玛多地震序列的发生对以上断裂带有很大影响, 大大的增加了上述断裂带滑动的可能; 甘孜—玉树断裂带西部及日月山断裂带的库仑应力变化为正, 说明玛多地震会加快以上断裂带滑动.

图 11 周边主要断裂带的同震库仑应力变化 Fig. 11 Variation of coseismic Coulomb stress of main faults around the Madoi earthquake
4 结论

在玛多地震发生的区域, Sentinel-1影像数据具有较好的覆盖, 本文利用较短时间基线的四景雷达影像数据, 通过干涉处理获得了玛多地震的升、降轨InSAR同震形变场, 得到了两个雷达视线方向的同震形变场.通过忽略南北向形变, 解算了东西方向和垂直方向的同震形变场.根据东西向地表的错动, 勾画了断层在地表的迹线, 将破裂面共分成6段用于数值模拟的分析.东西方向的地表形变表明, 此次地震发生了较大的东西向运动, 破裂迹线南侧向东、北侧向西运动, 且两侧的相对运动明显, 垂直向形变量级较小, 沿地表破裂迹线无明显地表错动.由此可知, 此次地震为左旋走滑为主的地震, 这与地震学所得结果一致.

本文模拟计算得到了玛多地震周边区域三维的地表形变, 通过与InSAR观测结果进行比较, 证实了模拟结果的可靠性.在此基础上, 给出了三维地表同震形变场, 可以看出玛多地震引起了较大的地表变形, 其中东边部分区域的地表形变较大, 错动量达到2.5 m, 中间有部分区域形变相对较小, 同时也可以看出玛多地震为典型的左旋走滑型地震, 模拟结果与震源机制及地质考察的结果一致.

本文以Sentinel-1升、降轨InSAR同震地表形变等多个结果作为约束, 给出了玛多地震引起的应力场变化, 结果显示玛多地震的余震都分布在库仑应力为负的区域, 说明玛多地震序列释放了发震断层上积累的应力; 而发震断层的西南和东北端库仑应力有所增加, 推测未来余震的发生可能会朝西南和东北方向发展.

为了探讨历史地震对玛多地震的影响, 本文模拟了周边5次7级以上强震对玛多地震断层面的影响.结果显示周边地块的加载使得玛多地震断层面的应变不断的累积, 2010年玉树地震对玛多地震的发生有很大的促进作用.本文也给出了玛多地震的发生对周边断层的影响, 模拟结果显示, 玛多地震造成了昆仑山口—江错断裂的西部、东昆仑断裂带中部、共和盆地南缘断裂、北缘断裂等多个应力强加载区(库仑应力变化明显大于0.01 MPa), 推测玛多地震的发生将大大增加上述断裂带滑动的可能; 玛多地震同时造成了甘孜—玉树断裂带西部、日月山段带等多个弱加载区(库仑应力>0), 由此推测玛多地震可能会对上述断裂带有一定的促进作用.

同时, 本文也存在一定的不足之处.本文采用软弱夹层反映断层的属性, 利用减低地震破裂区域的弹性模量的方式模拟地震引起的应力场的变化.因此本文的模拟结果主要是在InSAR观测形变的约束下, 反映玛多地震的孕育发震原因及玛多地震的发生可能会对周边断层的影响, 对于精细的模拟玛多地震的滑动特征, 接下来我们需要考虑具体的地震破裂模型.

致谢  本文采用的最新的波速和密度数据是中国地震局地球物理研究所石磊博士提供, 黏滞系数参数是国土资源部孙玉军博士提供.欧洲空间局(ESA)免费发布Sentinel-1A/1B SAR数据.在此表示感谢.同时感谢匿名审稿专家的建议.
References
Chen W, Qiao W X, Liu G, et al. 2018. Study on the coseismic slip model and Coulomb stress of the 2017 Jiuzhaigou MS7.0 earthquake constrained by GNSS and InSAR measurements. Chinese Journal of Geophysics (in Chinese), 61(5): 2122-2132. DOI:10.6038/cjg2018L0613
Deng Q D, Zhang P Z, Ran Y K, et al. 2002. Basic characteristics of active tectonics of China. Science in China Series D: Earth Sciences, 46(4): 356-372.
Fang J, Xu C J, Wen Y M, et al. 2019. Coseismic slip distribution of 2015 Gorkha (Nepal) MW7.8 earthquake determined using the Helmert variance component estimation. Chinese Journal of Geophysics (in Chinese), 62(3): 923-939. DOI:10.6038/cjg2019M0476
Farr T G, Rosen P A, Caro E, et al. 2007. The shuttle radar topography mission. Reviews of Geophysics, 45(2): RG2004. DOI:10.1029/2005RG000183
Freed A M. 2005. Earthquake triggering by static, dynamic, and postseismic stress transfer. Annual Review of Earth and Planetary Sciences, 33(1): 335-367. DOI:10.1146/annurev.earth.33.092203.122505
Fu B H, Awata Y, Du J G, et al. 2005. Complex geometry and segmentation of the surface rupture associated with the 14 November 2001 great Kunlun earthquake, northern Tibet, China. Tectonophysics, 407(1-2): 43-63. DOI:10.1016/j.tecto.2005.07.002
Funning G J, Parsons B, Wright T J, et al. 2005. Fielding, E. Surface displacements and source parameters of the 2003 Bam (Iran) earthquake from Envisat advanced synthetic aperture radar imagery. Journal of Geophysical Research, 110(B9): B09406. DOI:10.1029/2004JB003338
Harris R A. 1998. Introduction to special section: Stress triggers, stress shadows, and implications for seismic hazard. Journal of Geophysical Research: Solid Earth, 103(B10): 24347-24358. DOI:10.1029/98JB01576
Jónsson S, Zebker H, Segall P, et al. 2002. Fault slip distribution of the 1999 MW7.1 hector mine, california, earthquake, estimated from satellite radar and GPS measurements. Bulletin of the Seismological Society of America, 92(4): 1377-1389. DOI:10.1785/0120000922
Ji L Y, Liu C J, Xu J, et al. 2017. InSAR observation and inversion of the seismogenic fault for the 2017 Jiuzhaigou MS7.0 earthquake in China. Chinese Journal of Geophysics (in Chinese), 60(10): 4069-4082. DOI:10.6038/cjg20171032
Karplus M S, Zhao W, Klemperer S L, et al. 2011. Injection of Tibetan crust beneath the south Qaidam Basin: evidence from INDEPTH Ⅳ wide-angle seismic data. Journal of Geophysical Research: Solid Earth, 116(B7): B07301. DOI:10.1029/2010JB007911
King G C P, Stein R S, Lin J. 1994. Static stress changes and the triggering of earthquakes. Bulletin of the Seismological Society of America, 84(3): 935-953.
Kirby E, Harkins N, Wang E Q, et al. 2007. Slip rate gradients along the eastern Kunlun fault. Tectonics, 26(2): TC2010. DOI:10.1029/2006TC002033
Li C X, Dai H G, Chen Y M, et al. 2006. The discussion on some problems of the MS7.5 Tuosuohu Lake earthquake in 1937. Seismology and Geology (in Chinese), 28(1): 12-21. DOI:10.3969/j.issn.0253-4967.2006.01.002
Li C X, Xu X W, Wen X Z, et al. 2011. Rupture segmentation and slip partitioning of the mid-eastern part of the Kunlun Fault, north Tibetan Plateau. Science China Earth Sciences, 54(11): 1730-1745. DOI:10.1007/s11430-011-4239-5
Li P E, Liao L, Feng J Z, et al. 2019. Numerical simulation of relationship between stress evolution and strong earthquakes around the Bayan Har block since 1900. Chinese Journal of Geophysics (in Chinese), 62(11): 4170-4188. DOI:10.6038/cjg2019M0660
Li Q, Yao X Z, Yang S M, et al. 2012. A precise velocity field of tectonic deformation in China as inferred from intensive GPS observations. Science China Earth Sciences, 55(5): 695-698. DOI:10.1007/s11430-012-4412-5
Li Y J, Huang L Y, Ding R, et al. 2021. Coulombstress changes associated with the M7.3 Madoi earthquake and implications for seismic hazards. Natural Hazard Research, 1(2): 95-101. DOI:10.1016/j.nhres.2021.06.003
Liang M J, Yang Y, Du F, et al. 2020. Late quaternary activity of the central segment of the Dari fault and restudy of the surface rupture zone of the 1947 M73/4 Dari earthquake, Qinghai province. Seismology and Geology (in Chinese), 42(3): 703-714.
Liu G X. 1996. Eastern Kunlun active fault zone and its seismic activity. Earthquake Research in China (in Chinese), 12(2): 119-126.
Liu L, Li Y J, Zhu L Y, et al. 2021. Dari M7.7 earthquake on stress evolution along the boundary fault of the Bayan Har block: insights from numerical simulation. Chinese Journal of Geophysics (in Chinese), 64(7): 2221-2231. DOI:10.6038/cjg2021P0194
Liu M, Luo G, Wang H. 2014. The 2013 Lushan earthquake in China tests hazard assessments. Seismological Research Letters, 85(1): 40-43. DOI:10.1785/0220130117
Luo G, Liu M. 2010. Stress evolution and fault interactions before and after the 2008 Great Wenchuan Earthquake. Tectonophysics, 491(1-4): 127-140. DOI:10.1016/j.tecto.2009.12.019
Ma K F, Chan C H, Stein R S. 2005. Response of seismicity to Coulomb stress triggers and shadows of the 1999 MW=7.6 Chi-Chi, Taiwan, earthquake. Journal of Geophysical Research: Solid Earth, 110(B5): B05S19. DOI:10.1029/2004JB003389
Massonnet D, Rossi M, Carmona C, et al. 1993. The displacement field of the Landers earthquake mapped by radar interferometry. Nature, 364(6433): 138-142. DOI:10.1038/364138a0
Rundle J B, Rundle P B, Donnellan A, et al. 2006. Stress transfer in earthquakes, hazard estimation and ensemble forecasting: Inferences from numerical simulations. Tectonophysics, 413(1-2): 109-125. DOI:10.1016/j.tecto.2005.10.031
Sandwell D, Mellors R, Tong X P, et al. 2011. Open radar interferometry software for mapping surface deformation. Eos, Transactions American Geophysical Union, 92(28): 234. DOI:10.1029/2011EO280002
Shan B, Xiong X, Zheng Y, et al. 2009. Stress changes on major faults caused by MW7.9 Wenchuan earthquake, May 12, 2008. Science in China Series D: Earth Sciences, 52(5): 593-601. DOI:10.1007/s11430-009-0060-9
Shen W H, Li Y S, Jiao Q S, et al. 2019. Joint inversion of strong motion and InSAR/GPS data for fault slip distribution of the Jiuzhaigou 7.0 earthquake and its application in seismology. Chinese Journal of Geophysics (in Chinese), 62(1): 115-129. DOI:10.6038/cjg2019L0541
Shen W S, Ritzwoller M H, Kang D, et al. 2016. A seismic reference model for the crust and uppermost mantle beneath China from surface wave dispersion. Geophysical Journal International, 206(2): 954-979. DOI:10.1093/gji/ggw175
Shi Y L, Cao J L. 2010. Some aspects in static stress change calculation-case study on Wenchuan earthquake. Chinese Journal of Geophysics (in Chinese), 53(1): 102-110. DOI:10.3969/j.issn.0001-5733.2010.01.011
Song X G, Zhang Y F, Shan X J, et al. 2019. Geodetic observations of the 2018 MW7.5 Sulawesi earthquake and its implications for the kinematics of the Palu Fault. Geophysical Research Letters, 46(8): 4212-4220. DOI:10.1029/2019gl082045
Stein R S, King G C P, Lin J. 1992. Change in failure stress on the southern san andreas fault system caused by the 1992 Magnitude=7.4 Landers earthquake. Science, 258(5086): 1328-1332. DOI:10.1126/science.258.5086.1328
Stein S, Okal E A. 2011. The size of the 2011 Tohoku earthquake need not have been a surprise. Eos, Transactions American Geophysical Union, 92(27): 227-228.
Stolk W, Kaban M, Beekman F, et al. 2013. High resolution regional crustal models from irregularly distributed data: Application to Asia and adjacent areas. Tectonophysics, 602: 55-68. DOI:10.1016/j.tecto.2013.01.022
Sun J B, Xu X W, Shen Z K, et al. 2007. Parameter inversion of the 1997 Mani earthquake from INSAR co-seismic deformation field based on linear elastic dislocation model-Ⅰ. uniform slip inversion. Chinese Journal of Geophysics (in Chinese), 50(4): 1097-1110.
Sun XZ, Xu X W, Chen L C, et al. 2012. Surface rupture features of the 2010 Yushu earthquake and its tectonic implication. Chinese Journal of Geophysics (in Chinese), 55(1): 155-170. DOI:10.6038/j.issn.0001-5733.2012.01.015
Sun Y J, Dong S W, Fan T Y, et al. 2013. 3D rheological structure of the continental lithosphere beneath China and adjacent regions. Chinese Journal of Geophysics (in Chinese), 56(9): 2936-2946. DOI:10.6038/cjg20130908
Toda S, Lin J, Meghraoui M, et al. 2008. May 12 2008 M=7.9 Wenchuan, China, earthquake calculated to increase failure stress and seismicity rate on three major fault systems. Geophysical Research Letters, 35(17): L17305. DOI:10.1029/2008GL034903
Toda S, Lin J, Stein R S. 2011a. MW9.0 off the Pacific coast of Tohoku Earthquake to test the Coulomb stress triggering hypothesis and to calculate faults brought closer to failure. Earth Planets and Space, 63: 39. DOI:10.5047/eps.2011.05.010
Toda S, Stein R S, Lin J. 2011b. Widespread seismicity excitation throughout central Japan following the 2011 M=9.0 Tohoku earthquake and its interpretation by Coulomb stress transfer. Geophysical Research Letters, 38(7): L00G03. DOI:10.1029/2011GL047834
Wang M, Shen Z K. 2020. Present-day crustal deformation of continental China derived from GPS and its tectonic implications. Journal of Geophysical Research: Solid Earth, 125(2): e2019JB018774. DOI:10.1029/2019JB018774
Wang W L, Fang L H, Wu J P, et al. 2021. Aftershock sequence relocation of the 2021 MS7.4 Maduo Earthquake, Qinghai, China. Science China Earth Sciences, 64(8): 1371-1380. DOI:10.1007/s11430-021-9803-3
Wang Y Z, Feng W P, Chen K, et al. 2019. Source Characteristics of the 28 September 2018 MW7.4 Palu, Indonesia, earthquake derived from the advanced land observation satellite 2 data. Remote Sensing, 11(17): 1999. DOI:10.3390/rs11171999
Wang Y Z, Li Z W, Zhu J J, et al. 2012. Coseismic three-dimensional deformation of L'Aquila earthquake derived from multi-platform DInSAR data. Geomatics and Information Science of Wuhan University, 37(007): 859-863.
Wen Y M, Xiao Z H, He P, et al. 2021. Source characteristics of the 2020 MW7.4 Oaxaca, Mexico, earthquake estimated from GPS, InSAR, and teleseismic waveforms. Seismological Research Letters, 92(3): 1900-1912. DOI:10.1785/0220200313
Xia S R, Shi L, Li Y H, et al. 2021. Velocity structures of the crust and uppermost mantle beneath the northeastern margin of Tibetan plateau revealed by double-difference tomography. Chinese Journal of Geophysics (in Chinese), 64(9): 3194-3206. DOI:10.6038/cjg2021O0514
Xiong X, Shan B, Zheng Y, et al. 2010. Stress transfer and its implication for earthquake hazard on the Kunlun Fault, Tibet. Tectonophysics, 482(1-4): 216-225. DOI:10.1016/j.tecto.2009.07.020
Xu W B, Bürgmann R, Li Z W. 2016. An improved geodetic source model for the 1999 MW6.3 Chamoli earthquake, India. Geophysical Journal International, 205(1): 236-242. DOI:10.1093/gji/ggw016
Xu W B. 2017. Finite-fault slip model of the 2016 MW7.5 Chiloé earthquake, southern Chile, estimated from Sentinel-1 data. Geophysical Research Letters, 44(10): 4774-4780. DOI:10.1002/2017gl073560
Yan Z Z, Zhang H, Yang C C, et al. 2009. Spectral element analysis on the characteristics of seismic wave propagation triggered by Wenchuan MS8.0 earthquake. Science in China Series D: Earth Sciences, 52(6): 764-773. DOI:10.1007/s11430-009-0078-z
Yang Y F, Zhu J J, Wang Y Z, et al. 2017. 3D displacement field of 2016 Kaohsiung MS6.7 earthquake from D-InSAR and along-track interferometry with ascending and descending sentinel-1A images. Journal of Geodesy and Geodynamics (in Chinese), 37(4): 339-343.
Yang Y H, Hu J C, Chen Q, et al. 2021. Shallow slip of blind fault associated with the 2019 MS6.0 Changning earthquake in fold-and-thrust belt in salt mines of Southeast Sichuan, China. Geophysical Journal International, 224(2): 909-922. DOI:10.1093/gji/ggaa488
Yuan S, He P, Wen Y M, et al. 2020. Integrated InSAR and strain tensor to estimate three-dimensional coseismic displacements associated with the 2016 MW7.0 Kumamoto earthquake. Chinese Journal of Geophysics (in Chinese), 63(4): 1340-1356. DOI:10.6038/cjg2020N0308
Zhang P Z, Xu X W, Wen X Z, et al. 2008. Slip rates and recurrence intervals of the longmen Shan active fault zone, and tectonic implications for the mechanism of the May 12 Wenchuan earthquake, 2008, Sichuan, China. Chinese Journal of Geophysics (in Chinese), 51(4): 1066-1073.
Zhang Y, Xu L S, Chen Y T. 2010. Fast inversion of rupture process for 14 April 2010 Yushu, Qinghai, earthquake. Acta Seismologica Sinica (in Chinese), 32(3): 361-365.
Zhu S B, Zhang P Z. 2009. A study on the dynamical mechanisms of the Wenchuan MS8.0 earthquake, 2008. Chinese Journal of Geophysics (in Chinese), 52(2): 418-427.
Zhu Y G, Diao F Q, Fu Y C, et al. 2021. Slip rate of the seismogenic fault of the 2021 Madoi earthquake in western China inferred from GPS observations. Science China Earth Science, 64(8): 1363-1370. DOI:10.1007/s11430-021-9808-0
陈威, 乔学军, 刘刚, 等. 2018. 基于GNSS与InSAR约束的九寨沟MS7.0地震滑动模型及其库仑应力研究. 地球物理学报, 61(5): 2122-2132. DOI:10.6038/cjg2018L0613
邓起东, 张培震, 冉勇康, 等. 2002. 中国活动构造基本特征. 中国科学(D辑), 32(12): 1020-1030.
方进, 许才军, 温扬茂, 等. 2019. 基于方差分量估计的2015年尼泊尔MW7.8地震同震滑动分. 地球物理学报, 62(3): 923-939. DOI:10.6038/cjg2019M0476
季灵运, 刘传金, 徐晶, 等. 2017. 九寨沟MS7.0地震的InSAR观测及发震构造分析. 地球物理学报, 60(10): 4069-4082. DOI:10.6038/cjg20171032
李陈侠, 戴华光, 陈永明, 等. 2006. 对1937年托索湖7.5级地震若干问题的探讨. 地震地质, 28(1): 12-21. DOI:10.3969/j.issn.0253-4967.2006.01.002
李陈侠, 徐锡伟, 闻学泽, 等. 2011. 东昆仑断裂带中东部地震破裂分段性与走滑运动分解作用. 中国科学: 地球科学, 41(9): 1295-1310.
李平恩, 廖力, 奉建州, 等. 2019. 1900年以来巴颜喀拉块体应力演化与周缘强震关系的数值模拟研究. 地球物理学报, 62(11): 4170-4188. DOI:10.6038/cjg2019M0660
李旭, 陈运泰. 1996. 用长周期地震体波波形资料反演1990年青海共和地震的震源过程. 地震学报, 18(3): 279-286.
梁明剑, 杨耀, 杜方, 等. 2020. 青海达日断裂中段晚第四纪活动性与1947年M73/4地震地表破裂带再研究. 地震地质, 42(3): 703-714. DOI:10.3969/j.issn.0253-4967.2020.03.011
刘光勋. 1996. 东昆仑活动断裂带及其强震活动. 中国地震, 12(2): 119-126.
刘雷, 李玉江, 朱良玉, 等. 2021. 1947年达日M7.7地震对巴颜喀拉块体边界断裂应力影响的数值模拟. 地球物理学报, 64(7): 2221-2231. DOI:10.6038/cjg2021P0194
单斌, 熊熊, 郑勇, 等. 2009. 2008年5月12日MW7.9汶川地震导致的周边断层应力变化. 中国科学D辑: 地球科学, 39(5): 537-545.
申文豪, 李永生, 焦其松, 等. 2019. 联合强震记录和InSAR/GPS结果的四川九寨沟7.0级地震震源滑动分布反演及其地震学应用. 地球物理学报, 62(1): 115-129. DOI:10.6038/cjg2019L0541
石耀霖, 曹建玲. 2010. 库仑应力计算及应用过程中若干问题的讨论——以汶川地震为例. 地球物理学报, 53(1): 102-110. DOI:10.3969/j.issn.0001-5733.2010.01.011
孙建宝, 徐锡伟, 沈正康, 等. 2007. 基于线弹性位错模型及干涉雷达同震形变场反演1997年玛尼MW7.5级地震参数-Ⅰ.均匀滑动反演. 地球物理学报, 50(4): 1097-1110. DOI:10.3321/j.issn:0001-5733.2007.04.017
孙鑫喆, 徐锡伟, 陈立春, 等. 2012. 2010年玉树地震地表破裂带典型破裂样式及其构造意义. 地球物理学报, 55(1): 155-170. DOI:10.6038/j.issn.0001-5733.2012.01.015
孙玉军, 董树文, 范桃园, 等. 2013. 中国大陆及邻区岩石圈三维流变结构. 地球物理学报, 56(9): 2936-2946. DOI:10.6038/cjg20130908
王未来, 房立华, 吴建平, 等. 2021. 2021年青海玛多MS7.4地震序列精定位研究. 中国科学: 地球科学, 51(7): 1193-1202. DOI:10.1360/SSTe-2021-0149
王永哲, 李志伟, 朱建军, 等. 2012. 融合多平台DInSAR数据解算拉奎拉地震三维同震形变场. 武汉大学学报: 信息科学版, 37(007): 859-863.
夏思茹, 石磊, 李永华, 等. 2021. 青藏高原东北缘地壳及上地幔顶部速度结构研究. 地球物理学报, 64(9): 3194-3206. DOI:10.6038/cjg2021O0514
严珍珍, 张怀, 杨长春, 等. 2009. 汶川大地震地震波传播的谱元法数值模拟研究. 中国科学D辑: 地球科学, 39(4): 393-402.
杨亚夫, 朱建军, 王永哲, 等. 2017. 利用Sentinel-1A数据、D-InSAR和沿轨干涉技术获取2016年高雄MS6.7地震三维形变场. 大地测量与地球动力学, 37(4): 339-343.
袁霜, 何平, 温扬茂, 等. 2020. 综合InSAR和应变张量估计2016年MW7.0熊本地震同震三维形变场. 地球物理学报, 63(4): 1340-1356. DOI:10.6038/cjg2020N0308
张培震, 徐锡伟, 闻学泽, 等. 2008. 2008年汶川8.0级地震发震断裂的滑动速率、复发周期和构造成因. 地球物理学报, 51(4): 1066-1073. DOI:10.3321/j.issn:0001-5733.2008.04.015
张勇, 许力生, 陈运泰. 2010. 2010年4月14日青海玉树地震破裂过程快速反演. 地震学报, 32(3): 361-365. DOI:10.3969/j.issn.0253-3782.2010.03.011
朱守彪, 张培震. 2009. 2008年汶川MS8.0地震发生过程的动力学机制研究. 地球物理学报, 52(2): 418-427.
朱亚戈, 刁法启, 付誉超, 等. 2021. 基于GPS资料约束的2021年玛多地震发震断层的滑动速率. 中国科学: 地球科学, 51(10): 1788-1795. DOI:10.1360/SSTe-2021-0156