Triggering relationship in strong earthquake sequence around the Bayan Har block and its tectonic significance based on Coulomb stress changes and seismicity
-
摘要: 自本世纪以来，青藏高原巴颜喀拉块体周缘密集地发生了7次MW≥6.5强震事件，包括伤亡惨重的2008年汶川MW7.9地震和2017年九寨沟MW6.5地震。本文根据这7次强震事件中先前地震对后续地震的库仑应力改变（dCFS）的计算结果，结合基于ETAS模型得到的背景地震活动性的变化结果，研究了强震间的触发关系，以试图解释发生在环巴颜喀拉块体的几次强震的发震机理并探讨其构造意义。结果表明，各个强震的库仑应力变化与其造成的背景地震活动性变化呈正相关，并且在这7次强震中，汶川地震对芦山地震有显著的触发作用，并造成芦山地区背景地震活动性的显著提高。同时汶川地震对九寨沟地震的发生具有一定的延迟作用，造成九寨沟地区背景地震活动性降低。除此之外，其它地震之间均无明显的触发/延迟作用或显著的背景地震活动性变化。这表明该强震序列的孕震机制主要是巴颜喀拉块体东南向持续挤压的构造运动，推断巴颜喀拉块体目前仍处于构造运动活跃期，因此包括巴颜喀拉块体周缘在内的我国西南地区未来的强震危险性值得持续关注。Abstract: Seven major earthquakes of MW≥6.5 have occurred around the Bayan Har block in the Tibetan Plateau since this century, including the devastated 2008 Wenchuan earthquake and the recent 2017 Jiuzhaigou earthquake. Here we investigate the triggering effects among those major earthquakes by calculating the Coulomb stress changes and analyzing the background seismicity rate changes, and furthermore to understand the tectonic activity of the Bayan Har block. Our results show that, the Coulomb stress changes positively correlate with the background seismicity rate changes. Also, the results show that the Wenchuan earthquake triggered the Lushan earthquake and delayed the Jiuzhaigou earthquake. Besides, there is no significant triggering/delay effect among other major earthquakes. These findings imply that the remarkable earthquake sequence results from the southeastward movement of the Bayan Har block, and the Bayan Har block is still in the highly active seismic period. Furthermore, the seismic hazards of the Bayan Har block and its adjacent faults in the southwestern China should be paid more attention.
-
引言
自1964年电离层异常扰动在美国阿拉斯加大地震震前首次被发现（Leonard，Barnes，1965），地震-电离层耦合效应逐渐得到各国研究人员的广泛关注。对于地震电离层耦合机理的大量研究，异常电场理论能较好地解释震前电离层异常扰动现象。该理论认为，地震孕育过程中由地表岩石破裂引起的一系列复杂的物理化学反应，产生了大气电场异常，进而叠加作用于电离层，引起了电离层电子浓度、电子温度等参量的显著变化（Pulinets et al，2000；Sorokin et al，2001；Kuo et al，2011）。近年来，以全球导航卫星系统（Global Navigation Satellite System，简写为GNSS）连续观测和电离层垂直观测为代表的地基电离层探测技术已被广泛应用于地震前兆监测和研究中（张学民等，2016；吴健等，2016；何宇飞等，2020）。基于该技术，研究人员可获取电离层垂直总电子含量（vertical total electron content，简写为VTEC）、电离层F2层临界频率（foF2）等电离层形态特征的重要指标。其中，电离层VTEC可通过GNSS台站观测数据反演获取，foF2可通过电离层垂测站观测获取。大量研究结果表明，MS5.0以上强震发生前数天到数小时内，上述参量可能会出现不同程度的异常变化（张学民等，2014；王壮凯等，2020），而多参量同步分析可为地震-电离层异常扰动判识提供佐证，对于强震频发地区的震情研判具有重要意义。
川滇地区位于青藏高原东南缘，包括川滇地块、巴颜喀拉地块南部、华南地块西部等二级地块。长期以来，在印度板块的推挤及华南地块的阻挡作用下，川滇地区的构造变形活动强烈，发育了规模、产状和活动速率各异的断裂，致使强震频发。由于地震短临预报的困难性，地震电离层监测预报模式逐渐成为川滇地区地震监测预报的新方向。目前，地震-电离层现象的研究主要处于震例研究和统计研究阶段（闫相相等，2014），区域电离层实时监测系统开发实例较少，电离层多参量异常关联及时空演化分析不够深入，快速获取电离层多参量异常扰动信息及其临震监测预报的应用还有待加强。
鉴于此，本文设计开发川滇地区的电离层多参量异常监测系统，旨在实现各类观测数据的自动化处理、电离层多参量图形及数据结果的自动化产出以及川滇地区电离层多参量异常变化的实时监测，为该区域地震电离层异常扰动判识提供证据。同时，本文拟通过2019年6月17日四川长宁MS6.0地震分析测试该系统图形产出功能的可靠性，以进一步验证川滇地区地震电离层监测预报的应用效能。
1. 系统概述
为了实现川滇地区电离层异常扰动变化的实时监测，本文设计开发川滇地区电离层多参量异常监测系统，快速产出全球和中国区域VTEC异常分布、空间参数、站点VTEC和foF2时序曲线，同时实现网页方式下图形及数据结果的查询、浏览、对比和下载等功能。基于该系统可开展示范性地震电离层监测工作，分析地震发生前后孕震区附近电离层多参量的异常扰动变化，尝试探索地震电离层现象的内在联系，为该地区地震监测预报提供电离层前兆异常的判定依据。
1.1 数据处理
1） 全球和中国区域VTEC异常分布。在中国区域按照均匀分布原则选取陆态网络的70个GNSS基准站。基于该基准站的观测数据，利用球谐函数建立大规模高精度格网电离层模型（林剑等，2009；王泽民等，2016），利用Bernese软件解算处理（Dach et al，2015）获取中国区域电离层VTEC格网数据（经纬度分辨率为1°×1°，时间分辨率为2 h），并结合欧洲轨道确定中心（Center for Orbit Determination in Europe，缩写为CODE）提供的全球电离层格网数据（经纬度分辨率为5°×2.5°，时间分辨率为1 h）自动监测全球和中国区域VTEC异常分布特征变化。
2） 空间参数、站点VTEC和foF2时序曲线。采用川滇地区GNSS基准站观测数据，通过单站点VTEC解算（熊波，2012），获取站点上空的VTEC数据（时间分辨率为15 min）。结合川滇地区垂测站（乐山、道孚、普洱、腾冲）的foF2观测数据（时间分辨率为15 min）、空间环境参数（F10.7，Kp，Dst，AE），自动监测站点VTEC、站点foF2以及其它空间参数的时序曲线特征变化。川滇地区GNSS、垂测台站的分布如图1所示。
![]()
图 1 川滇地区GNSS和垂测台站分布图Figure 1. Distribution of GNSS and vertical ionosonde stations in Sichuan-Yunnan region1.2 异常检验
选用滑动四分位距法提取电离层VTEC和foF2的异常信息，基于前15天相同时刻的观测值求取相应四分位距（IQR）及中位数（M），确定异常检验倍数（N），建立上（下）边界阈值M±N×IQR，观测值超出上（下）边界即视为异常，高于上边界为正异常，低于下边界为负异常。然后，从二维空间分布获取全球和中国区域VTEC异常分布，所用公式如下：
$$ \Delta {\rm{TEC}}\text{\text{＝}}\left\{\begin{array}{l}K{\text{－}}{L}_{1}\qquad K {\text{＞}} {L}_{1}{\text{，}}\\ K{\text{－}}{L}_{2}\qquad K {\text{＜}} {L}_{2}{\text{，}}\\ 0\;\qquad\qquad{L}_{2}{\text{≤}} K{\text{≤}} {L}_{1}{\text{，}}\end{array} \right.$$ (1) 式中，ΔTEC为VTEC异常扰动变化量，L1和L2分别为VTEC背景上、下边界阈值，K为VTEC观测值。站点VTEC和foF2时序曲线的异常变化获取方法同上。
2. 系统设计
2.1 结构设计
为满足川滇地区电离层多参量异常扰动监测的总体需求，结合数据获取、数据解算、图形及数据结果展示等需求，采用CS＋BS架构将该系统分为电离层VTEC异常分布产出子系统（简写为“异常分布子系统”）、GNSS站点VTEC及垂测站foF2曲线产出子系统（简写为“时序曲线子系统”）和电离层多参量异常监测平台（简写为“监测平台”）等三部分。其中：异常分布子系统用于监测全球和中国区域的电离层VTEC异常变化，实时产出全球和中国区域VTEC异常分布；时序曲线子系统用于监测川滇地区GNSS单站点VTEC及垂测站foF2的异常时序变化，实时产出空间参数、站点VTEC和foF2时序曲线；监测平台将各类图形及数据结果以网页方式自动发布，方便用户按照权限查询、浏览、对比和下载。两个子系统布设在不同服务器上，独立运行，方便维护管理。系统总体结构如图2所示。
![]()
图 2 川滇地区电离层多参量异常监测系统的总体结构图Figure 2. Overall structure diagram of the ionospheric multi-parameter anomaly monitoring system in Sichuan-Yunnan region2.2 功能设计
基于系统结构设计、业务需求及日常监测工作的要求，系统功能模块划分如图3所示。可见：异常分布子系统分为参数设置、自动产出、手动产出、数据服务、任务列表等五部分，主要实现全球和中国区域VTEC异常分布图形结果产出，同时提供数据转换、数据导出、数据下载等功能；时序曲线子系统分为参数设置、自动产出、手动产出、数据服务、任务列表等五部分，主要实现站点VTEC和foF2、空间参数时序曲线图形结果的产出，同时提供数据导出、下载等功能。监测平台分为系统管理、图形展示、图形对比、数据下载等四部分，主要实现VTEC异常分布，空间参数、站点VTEC和foF2时序曲线查询、浏览，同时提供电离层多参量数据分析结果的下载功能。
![]()
图 3 川滇地区电离层多参量异常监 测系统功能图Figure 3. The function diagram of the ionospheric multi-parameter anomaly monitoring system in Sichuan-Yunnan region2.3 业务流程
异常分布子系统。该子系统可分别通过自动和手动模式以FTP方式下载各类数据，自动调用Bernese数据处理软件，基于球谐函数建立中国区域电离层格网模型，解算中国区域电离层VTEC格网数据。结合CODE全球电离层格网数据和空间环境参数（F10.7，Kp，Dst，AE），提取全球和中国区域VTEC异常信息，利用ArcGIS平台分别产出全球和中国区域VTEC异常分布以及格网点VTEC时序曲线，自动产出结果延迟1—2天，具体业务流程如图4所示。
时序曲线子系统。该子系统通过自动和手动模式以FTP方式下载各类数据，自动调用GPS-TEC数据处理软件，获取单站点VTEC数据。结合川滇地区垂测站（乐山、道孚、普洱、腾冲） foF2观测数据、空间环境参数（F10.7，Kp，Dst，AE），提取站点VTEC和foF2异常信息，并利用Matlab平台分别产出站点VTEC和foF2、空间参数的时序曲线，自动产出结果延迟1—2天，具体业务流程如图5所示。
监测平台。该平台实时访问异常分布子系统和时序曲线子系统，将各类图形、数据产出结果通过http自动发布，根据用户权限为用户提供查询、浏览、对比、下载等差异化功能。
3. 系统实现和震例分析
3.1 系统实现
该系统性能要求较高，需要实现各类数据的自动化下载、多个第三方软件的自动化运行和监测、各类图形及数据结果的自动化生成与发布。为实现系统运行稳定性、人机交互性及功能扩展性等多方面要求，具体实现内容如下：
1） 异常分布子系统。基于Windows 2008/2012操作系统，采用C/S （Client-Server，即客户端/服务器模式）架构，将面向对象的编程语言C#作为后台开发语言，MySQL作为后台数据库管理系统，ArcGIS作为图形结果产出平台，通过分层技术体系结构实现系统运行的稳定性和交互性。子系统主界面分为功能菜单、用户操作及图形展示等三部分，其中：功能菜单位于主界面左上方，包含系统配置、手动产出、时序曲线、数据转换、数据导出、数据下载、任务列表、站点管理等功能；用户操作位于功能菜单左下方，包含查询条件和查询结果两部分；图形展示位于用户操作界面的右侧，根据查询条件和查询结果查看电离层VTEC异常分布情况，全球VTEC异常分布间隔一小时，中国区域VTEC异常分布间隔两小时（图6a）。
![]()
图 6 系统主界面（a） 异常分布子系统；（b） 时序曲线子系统；（c） 监测平台Figure 6. Main system interfaces（a） Anomaly distribution subsystem；（b） Time series curve subsystem；（c） Monitoring platform2） 时序曲线子系统。基于Windows 2008/2012操作系统，采用C/S架构，将C#作为后台开发语言，MySQL作为后台数据库管理系统，Matlab作为图形结果产出平台，通过分层技术体系结构实现系统运行的稳定性和交互性。子系统主界面分为功能菜单、用户操作及图形展示等三部分，其中：功能菜单位于主界面左上方，包含手动产出、数据导出、数据下载、任务列表、站点设置、系统配置等功能；用户操作位于功能菜单左下方，包含查询条件和查询结果两部分；图形展示位于用户操作界面的右侧，根据查询条件和查询结果分别展示站点VTEC和foF2、空间参数时序曲线（图6b）。
3） 监测平台。基于Windows操作系统，采用B/S （Browser-Server，即浏览器/服务器模式）架构，将C#作为后台开发语言。该平台设置用户权限，通过网页方式登录。功能菜单位于主界面左侧，包含系统管理、VTEC异常分布、空间参数曲线、站点VTEC曲线、站点foF2曲线、图形对比和数据下载等功能。选择功能菜单，打开用户操作界面，根据查询条件和查询结果分别展示VTEC异常分布及空间参数、站点VTEC和foF2的时序曲线。同时，平台提供不同时期的电离层多参量图形结果对比展示功能，方便用户针对电离层前兆异常分析判断（图6c）。
3.2 震例分析
为测试系统图形产出功能的可靠性，验证川滇地区地震电离层监测预报的可行性，本文选择2019年6月17日四川长宁MS6.0地震前后的观测数据，利用该系统手动产出全球和中国区域VTEC异常分布及站点VTEC、foF2和空间参数的时序曲线，综合分析长宁地震震前电离层多参量异常扰动变化特征。鉴于篇幅有限，本文仅展示个别图形结果（图7，8）。
![]()
图 7 异常分布子系统关于长宁地震前2019年6月14日12:00 UT全球（a）和中国区域（b）VTEC异常分布图形结果，图中红色圆点代表震中Figure 7. Global （a） and China regional （b） VTEC anomaly distribution graphic results of the anomaly distribution subsystem at 12:00 UT on June 14，2019 before Changning earthquake，where red dots represent the epicenter从全球和中国区域VTEC异常分布（图7）可以看出，6月14日12:00 UT （北京时间20:00），震中附近区域均出现显著的电离层VTEC正异常，地震电离层耦合的局部性、集中性和显著性特征明显。同时，以个别站点为例，站点VTEC （图8a）和站点foF2 （图8b）的时序曲线显示6月14日四川宁南GNSS站点VTEC出现正异常，云南腾冲垂测站foF2同步出现正异常，该同步异常现象与VTEC异常分布结果基本吻合。结合空间参数时序曲线变化（图8c）可知，6月14日空间环境参数F10.7为68 sfu，Kp最大值为2.3，Dst最小值为−12 nT，AE最大值为338 nT，均未超过参数阈值，可以排除太阳活动、地磁变化对电离层异常扰动的影响。同时，长宁地震发生前后该异常区域并未发生其它强震事件。综上，2019年6月14日震中附近区域电离层多参量异常扰动满足地震电离层异常判识标准，可能是地震前兆异常，与长宁MS6.0地震具有一定的相关性。
![]()
图 8 时序曲线子系统图形结果（2019年6月1—20日）（a） 四川宁南站点VTEC时序曲线；（b） 云南腾冲站点foF2时序曲线；（c） 空间参数时序曲线，图中红色直线为参数阈值Figure 8. Graphic results of time series curve subsystem （June 1 to June 20，2019）（a） VTEC time series curves of Ningnan station in Sichuan Province；（b） foF2 time series curves of Tengchong station in Yunnan Province；（c） Time series curves of spatial parameter，where all red lines are parameter thresholds4. 讨论与结论
本文主要介绍了川滇地区电离层多参量异常监测系统的结构设计、功能设计和业务流程设计，并给出了系统实现结果。该系统实现了全球和中国区域VTEC、站点VTEC和foF2异常变化的实时监测，有助于开展示范性地震-电离层异常监测工作，目前该系统已成功应用于四川省地震局的监测预报日常工作，图形产出结果已参与地震监测预报会商。
该系统虽然较大地提高了川滇地区电离层异常变化监测能力，但对电离层参量的监测手段较为单一，电离层异常对比分析功能较弱。该系统下一步仍需不断补充完善，基于地基和空基电离层探测技术开展天地一体化联合观测，并结合临震电离层多参量同步分析为各参量异常指标提供量化信息，为地震监测任务提供更便捷的技术支撑，为深入研究地震孕育期、临震、主震及余震中的电离层异常扰动机理提供便捷的数据支撑。
CODE提供了全球电离层格网数据，中国地震台网中心和中国地震局地震预测研究所提供了GNSS台站和电离层垂测台站的观测数据，作者在此一并表示衷心的感谢。
-
![]()
图 1 巴颜喀拉块体构造背景及7次强震事件（活动断层数据引自张培震等，2003）
Figure 1. The tectonic background of the Bayan Har block and the seven major earthquakes（The active faults data is from Zhang et al，2003 ）
![]()
图 2 强震之间的库仑应力改变随时间的变化
（a） 汶川地震对芦山地震的库仑应力改变；（b） 汶川地震对九寨沟地震的库仑应力改变；（c） 芦山地震对九寨沟地震的库仑应力改变
Figure 2. Coulomb stress changes dCFS of strong earthquakes versus time
（a） The dCFS induced by the Wenchuan earthquake on the Lushan earthquake；（b） The dCFS induced by the Wenchuan earthquake on the Jiuzhaigou earthquake；（c） The dCFS induced by the Lushan earthquake on the Jiuzhaigou earthquake
![]()
图 3 2008年于田地震对2014年于田地震的库仑应力改变量
Figure 3. The dCFS induced by the 2008 Yutian earthquake on the 2014 Yutian earthquake
![]()
图 4 青藏高原北部地震活动性空间分布及各研究区空间范围
Figure 4. Seismicity of the northern Tibetan Plateau and the spatial region of the studied area
![]()
图 5 累积地震次数和背景地震概率
（a） 可可西里；（b） 于田；（c） 汶川；（d） 玉树；（e） 芦山；（f） 九寨沟
Figure 5. The cumulative total earthquakes and the cumulative background earthquake probabilities
（a） Hoh Xil；（b） Yutian；（c） Wenchuan；（d） Yushu；（e） Lushan；（f） Jiuzhaigou
表 1 巴颜喀拉块体周缘7次强震的震源机制解
Table 1 The focal mechanism solutions of seven major earthquakes occurred around the Bayan Har block
序号 发震时间
年-月-日东经/° 北纬/° MW 走向/° 倾角/° 滑动角/° 深度/km 1 2001−11−14 92.91 35.80 7.8 94 51 −12 15.0 2 2008−03−21 81.37 35.43 7.1 203 52 −74 12.0 3 2008−05−12 104.10 31.44 7.9 231 35 138 12.8 4 2010−04−14 96.79 33.05 6.9 300 88 23 15.7 5 2013−04−20 103.12 30.22 6.6 212 42 100 21.9 6 2014−02−12 82.57 36.22 6.9 242 86 −5 18.3 7 2017−08−08 103.90 33.21 6.5 150 78 −13 14.9 注：数据参考Ekström等 （2012）。 
下载: 导出CSV
表 2 巴颜喀拉强震序列在其后续地震发震位置所引起的库仑应力变化
Table 2 The dCFS results induced by the strong earthquake sequence occurred around the Bayan Har block on following earthquakes
2001年可可西里地震 2010年玉树地震 同震应力改变
/MPa震后应力改变
/MPa同震应力改变
/MPa震后应力改变
/MPa于田地震 7.00×10−5 1.30×10−4 7.28×10−7 8.00×10−7 汶川地震 1.14×10−5 1.20×10−5 −1.32×10−4 −2.50×10−4 玉树地震 −7.15×10−4 −1.35×10−3 注：震后应力改变指同震和震后黏弹性松弛共同导致的结果。
下载: 导出CSV
-
陈运泰, 杨智娴, 张勇, 刘超. 2013. 从汶川地震到芦山地震[J]. 中国科学: 地球科学, 43（6）: 1064-1072. Chen Y T, Yang Z X, Zhang Y, Liu C. 2013. From 2008 Wenchuan earthquake to 2013 Lushan earthquake[J]. Scientia Sinica （Terrae）, 43（6）: 1064-1072(in Chinese).
程佳, 刘杰, 甘卫军, 余怀忠. 2011. 1997年以来巴颜喀拉块体周缘强震之间的黏弹性触发研究[J]. 地球物理学报, 54（8）: 1997-2010. Cheng J, Liu J, Gan W J, Yu H Z. 2011. Coulomb stress interaction among strong earthquakes around the Bayan Har block since the Manyi earthquake in 1997[J]. Chinese Journal of Geophysics, 54（8）: 1997-2010（in Chinese）.
刁桂苓, 王晓山, 高国英, 聂晓红, 冯向东. 2010. 以震源机制类型划分汶川、玉树地震构造块体归属[J]. 地球物理学报, 53（8）: 1778-1783. Diao G L, Wang X S, Gao G Y, Nie X H, Feng X D. 2010. Tectonic block attribution of Wenchuan and Yushu earthquakes distinguished by focal mechanism type[J]. Chinese Journal of Geophysics, 53（5）: 1778-1783（in Chinese）.
杜方, 龙锋, 阮祥, 易桂喜, 宫悦, 赵敏, 张致伟, 乔慧珍, 汪智, 吴江. 2013. 四川芦山7.0级地震及其与汶川8.0级地震的关系[J]. 地球物理学报, 56（5）: 1772-1783. Du F, Long F, Ruan X, Yi G X, Gong Y, Zhao M, Zhang Z W, Qiao H Z, Wang Z, Wu J. 2013. The M7.0 Lushan earthquake and the relationship with the M8.0 wenchuan earthquake in sichuan, China[J]. Chinese Journal of Geophysics, 56（5）: 1772-1783（in Chinese）.
蒋长胜, 庄建仓. 2010. 基于时-空ETAS模型给出的川滇地区背景地震活动和强震潜在危险区[J]. 地球物理学报, 53（2）: 305-317. Jiang C S, Zhuang J C. 2010. Evaluation of background seismicity and potential source zones of strong earthquakes in the Sichuan-Yunan region base on the space-time ETAS model[J]. Chinese Journal of Geophysics, 53（2）: 305-317（in Chinese）.
单斌, 熊熊, 郑勇, 刁法启. 2009. 2008年5月12日MW7.9汶川地震导致的周边断层应力变化[J]. 中国科学: D辑, 39（5）: 537-545. Shan B, Xiong X, Zheng Y, Diao F Q. 2009. Stress changes on major faults caused by Mw7.9 Wenchuan earthquake, May 12, 2008[J]. Science in China: Series DEarth, 52（5）: 593-601(in Chinese). doi: 10.1007/s11430-009-0060-9.
单斌, 郑勇, 刘成利, 谢祖军, 孔军. 2017. 2017年M7.0级九寨沟地震同震库仑应力变化及其与2008年汶川地震的关系[J]. 中国科学: 地球科学, 47（11）: 1329-1338. doi: 10.1360/N072017-00268. Shan B, Zheng Y, Liu C L, Xie Z J, Kong J. 2017. Coseismic coulomb failure stress changes caused by the 2017 M7.0 Jiuzhaigou earthquake, and its relationship with the 2008 Wenchuan earthquake[J]. Science in China Earth Sciences, 60（12）: 2181-2189(in Chinese). doi: 10.1007/s11430-017-9125-2.
王辉, 曹建玲, 洪顺英, 徐岳仁, 荆凤. 2016. 2008年和2014年2次新疆于田M7地震之间的黏弹性应力转移[J]. 地震地质, 38（3）: 646-659. Wang H, Cao J L, Hong S Y, Xu Y R, Jing F. 2016. Viscoelastic stress transfer between 2008 and 2014 Yutian M7 earthquakes, Xinjiang[J]. Seismology and Geology, 38（3）: 646-659（in Chinese）.
王为民, 郝金来, 姚振兴. 2013. 2013年4月20日四川芦山地震震源破裂过程反演初步结果[J]. 地球物理学报, 56（4）: 1412-1417. doi: 10.6038/cjg20130436. Wang W M, Hao J L, Yao Z X. 2013. Preliminary result for rupture process of Apr. 20, 2013, Lushan Earthquake, Sichuan, China[J]. Chinese Journal of Geophysics, 56（4）: 1412-1417（in Chinese）. doi: 10.6038/cjg20130436.
闻学泽, 杜方, 张培震, 龙锋. 2011. 巴颜喀拉块体北和东边界大地震序列的关联性与2008年汶川地震[J]. 地球物理学报, 54（3）: 706-716. Wen X Z, Du F, Zhang P Z, Long F. 2011. Correlation of major earthquake sequences on the northern and eastern boundaries of the Bayan Har block, and its relation to the 2008 Wenchuan earthquake[J]. Chinese Journal of Geophysics, 54（3）: 706-716（in Chinese）.
徐锡伟, 陈桂华, 王启欣, 陈立春, 任治坤, 许冲, 魏占玉, 鲁人齐, 谭锡斌, 董绍鹏, 石峰. 2017. 九寨沟地震发震断层属性及青藏高原东南缘现今应变状态讨论[J]. 地球物理学报, 60（10）: 4018-4026. Xu X W, Chen G H, Wang Q X, Chen L C, Ren Z K, Xu C, Wei Z Y, Lu R Q, Tan X B, Dong S P, Shi F. 2017. Discussion on seismogenic structure of Jiuzhaigou earthquake and its implication for current strain state in the southeastern Qinghai-Tibet Plateau[J]. Chinese Journal of Geophysics, 60（10）: 4018-4026（in Chinese）.
张国宏, 屈春燕, 单新建, 张桂芳, 宋小刚, 汪荣江, 李振洪, 胡植庆. 2011. 2008年MS7.1于田地震InSAR同震形变场及其震源滑动反演[J]. 地球物理学报, 54（11）: 2753-2760. Zhang G H, Qu C Y, Shan X J, Zhang G F, Song X G, Wang R J, Li Z H, Hu Z Q. 2011. The Coseismic InSAR measurements of 2008 Yutian earthquake and its inversion for source parameters[J]. Chinese Journal of Geophysics, 2011, 54（11）: 2753-2760（in Chinese）.
张培震, 邓起东, 张国民, 马瑾, 甘卫军, 闵伟, 毛凤英, 王琪. 2003. 中国大陆的强震活动与活动地块[J]. 中国科学（D辑）, 33（增刊1）: 12-20. doi: 10.3321/j.issn:1006-9267.2003.z1.002. Zhang P Z, Deng Q D, Zhang G M, Ma J, Gan W J, Min W, Mao F Y, Wang Q. 2003. Active tectonic blocks and strong earthquakes in the continent of China[J]. Science in China Series D: Earth Sciences, 46（S2）: 13-24(in Chinese).
赵立波, 赵连锋, 谢小碧, 曹俊兴, 姚振兴. 2016. 2014年2月12日新疆于田Mw7.0地震源区静态库仑应力变化和地震活动率[J]. 地球物理学报, 59（10）: 3732-3743. Zhao L B, Zhao L F, Xie X B, Cao J X, Yao Z X. 2016. Static Coulomb stress changes and seismicity rate in the source region of the 12 February, 2014 Mw7.0 Yutian earthquake in Xinjiang, China[J]. Chinese Journal of Geophysics, 59（10）: 3732-3743（in Chinese）.
Aron A, Hardebeck J L. 2009. Seismicity rate changes along the central california coast due to stress changes from the 2003 M 6.5 San Simeon and 2004 M 6.0 Parkfield earthquakes[J]. Bull Seismol Soc Am, 99（4）: 2280-2292.
Beeler N M, Simpson R W, Hickman S H, Lockner D A. 2000. Pore fluid pressure, apparent friction, and Coulomb failure[J]. J Geophys Res, 105（B11）: 25533-25542. doi: 10.1029/2000jb900119.
Cocco M, Hainzl S, Catalli F, Enescu B, Lombardi A M, Woessner J. 2010. Sensitivity study of forecasted aftershock seismicity based on Coulomb stress calculation and rate- and state-dependent frictional response[J]. J Geophy Res, 115（B5）: B05307. doi: 10.1029/2009JB006838.
Dieterich J H. 1979. Modeling of rock friction: 1. Experimental results and constitutive equations[J]. J Geophys Res, 84（B5）: 2161-2168.
Dieterich J. 1994. A constitutive law for rate of earthquake production and its application to earthquake clustering[J]. J Geophys Res, 99（B2）: 2601-2618.
Ekström G, Nettles M, Dziewonski A M. 2012. The global CMT project 2004−2010: Centroid-moment tensors for 13 017 earthquakes[J], Phys Earth Planet Inter, 200/201: 1-9. doi: 10.1016/j.pepi.2012.04.002.
Elliott J R, Walters R J, England P C, Jackson J A, Li Z, Parsons B. 2010. Extension on the Tibetan plateau: Recent normal faulting measured by InSAR and body wave seismology[J]. Geophys J Int, 183（2）: 503-535.
Enescu B, Mori J, Miyazawa M. 2007. Quantifying early aftershock activity of the 2004 mid Niigata Prefecture earthquake （Mw6.6）[J]. J Geophys Res, 112（B4）: B04310.
Freed A M, Lin J. 2001. Delayed triggering of the 1999 Hector Mine earthquake by viscoelastic stress transfer[J]. Nature, 411（6834）: 180-183.
Harris R A. 1998. Introduction to special section: stress triggers, stress shadows, and implications for seismic hazard[J]. J Geophys Res, 103（B10）: 24347-24358.
Hearn E H. 2002. Dynamics of izmit earthquake postseismic deformation and loading of the duzce earthquake hypocenter[J]. Bull Seismol Soc Am, 92（1）: 172-193.
Hu J H, Fu L Y, Sun W J, Zhang Y. 2017. A study of the Coulomb stress and seismicity rate changes induced by the 2008 MW7.9 Wenchuan earthquake, SW China[J]. J Asian Earth Sci, 135: 303-319. doi: 10.1016/j.jseaes.2016.12.048.
Ishibe T, Satake K, Sakai S, Shimazaki K, Tsuruoka H, Yokota Y, Nakagawa S, Hirata N. 2015. Correlation between Coulomb stress imparted by the 2011 Tohoku-Oki earthquake and seismicity rate change in Kanto, Japan[J]. Geophys J Int, 201（1）: 112-134.
Ji C, Hayes G. 2008. Preliminary result of the May 12, 2008 MW7.9 eastern Sichuan, China earthquake[EB/OL]. [2013–08–01]. http://www.oalib.com/references/19014475. Jia K, Zhou S Y, Wang R. 2012. Stress interactions within the strong earthquake sequence from 2001 to 2010 in the bayankala block of Eastern Tibet[J]. Bull Seismol Soc Am, 102（5）: 2157-2164.
Jia K, Zhou S Y, Zhuang J C, Jiang C S. 2014. Possibility of the independence between the 2013 Lushan earthquake and the 2008 Wenchuan earthquake on Longmen shan fault, Sichuan, China[J]. Seismol Res Lett, 85（1）: 60-67. doi: 10.1785/0220130115.
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, 123. doi: 10.1002/2017JB015165.
Jiang Z S, Wang M, Wang Y Z, Wu Y Q, Che S, Shen Z K, Bürgmann R, Sun J B, Yang Y L, Liao H, Li Q. 2014. GPS constrained coseismic source and slip distribution of the 2013 Mw6.6 Lushan, China, earthquake and its tectonic implications[J]. Geophys Res Lett, 41（2）: 407-413.
King G C P, Stein R S, Lin J. 1994. Static stress changes and the triggering of earthquakes[J]. Bull Seismol Soc Am, 84（3）: 935-953.
King G C P, Cocco M. 2001. Fault interaction by elastic stress changes: new clues from earthquake sequences[J]. Adv Geophys, 44: 1-38. doi: 10.1016/S0065-2687(00)80006-0.
Lasserre C, Peltzer G, Crampé F, Klinger Y, Van Der Woerd J, Tapponnier P. 2005. Coseismic deformation of the 2001 MW=7.8 Kokoxili earthquake in Tibet, measured by synthetic aperture radar interferometry[J]. J Geophys Res, 110（B12）: B12408.
Li H B, Xue L, Brodsky E E, Mori J J, Fulton P M, Wang H, Kano Y, Yun K, Harris R N, Gong Z, Li C L, Si J L, Sun Z M, Pei J L, Zheng Y, Xu Z Q. 2015a. Long-term temperature records following the Mw 7.9 Wenchuan（China） earthquake are consistent with low friction[J]. Geology, 43（2）: 163-166. doi: 10.1130/g35515.1.
Li Y J, Chen L W, Liu S F, Yang S X, Yang X, Zhang G G. 2015b. Coseismic Coulomb stress changes caused by the Mw6.9 Yutian earthquake in 2014 and its correlation to the 2008 Mw7.2 Yutian earthquake[J]. J Asian Earth Sci, 105: 468-475.
Luo G, Liu M. 2010. Stress evolution and fault interactions before and after the 2008 Great Wenchuan earthquake[J]. Tectonophysics, 491（1/4）: 127-140.
Nalbant S S, McCloskey J. 2011. Stress evolution before and after the 2008 Wenchuan, China earthquake[J]. Earth Planet Sci Lett, 307（1/2）: 222-232. doi: 10.1016/j.jpgl.2011.04.039.
Ogata Y. 1988. Statistical models for earthquake occurrences and residual analysis for point processes[J]. J Am Stat Assoc, 83（401）: 9-27.
Ogata Y. 1998. Space-time point-process models for earthquake occurrences[J]. Ann Inst Stat Math, 50（2）: 379-402. doi: 10.1023/A:1003403601725.
Ogata Y, Zhuang J C. 2006. Space–time ETAS models and an improved extension[J]. Tectonophysics, 413（1/2）: 13-23. doi: 10.1016/j.tecto.2005.10.016.
Parsons T, Ji C, Kirby E. 2008. Stress changes from the 2008 Wenchuan earthquake and increased hazard in the Sichuan basin[J]. Nature, 454（7203）: 509-510. doi: 10.1038/nature07177.
Peng Y, Zhou S, Zhuang J, Shi J. 2012. An approach to detect the abnormal seismicity increase in Southwestern China triggered co-seismically by 2004 Sumatra MW9.2 earthquake[J]. Geophys J Int, 189（3）: 1734-1740.
Perfettini H, Stein R S, Simpson R, Cocco M. 1999. Stress transfer by the 1988-1989 M5.3 and 5.4 Lake Elsman foreshocks to the Loma Prieta fault: unclamping at the site of peak mainshock slip[J]. J Geophys Res, 104（B9）: 20169-20182. doi: 10.1029/1999jb900092.
Pollitz F, Vergnolle M, Calais E. 2003. Fault interaction and stress triggering of twentieth century earthquakes in Mongolia[J]. J Geophys Res, 108（B10）: 2503. doi: 10.1029/2002jb002375.
Reasenberg P A, Simpson R W. 1992. Response of regional seismicity to the static stress change produced by the loma prieta earthquake[J]. Science, 255（5052）: 1687-1690.
Robinson R. 2005. Stress Interactions within the Tangshan, China, Earthquake Sequence of 1976[J]. Bull Seismol Soc Am, 95（6）: 2501-2505. doi: 10.1785/0120050091.
Ruina A. 1983. Slip instability and state variable friction laws[J]. J Geophys Res, 88（B12）: 10359-10370.
Steacy S, Gomberg J, Cocco M. 2005. Introduction to special section: stress transfer, earthquake triggering, and time-dependent seismic hazard[J]. J Geophys Res, 110（B5）: B05S01. doi: 10.1029/2005JB003692.
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[J]. Science, 258（5086）: 1328-1332.
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 M7.9 Wenchuan, China, earthquake calculated to increase failure stress and seismicity rate on three major fault systems[J]. Geophys Res Lett, 35（17）: L17305. doi: 10.1029/2008GL034903.
Toda S, Enescu B. 2011. Rate/state Coulomb stress transfer model for the CSEP Japan seismicity forecast[J]. Earth Planets Space, 63（3）: 171-185.
Utsu T, Ogata Y, Ritsuko S, Matsu'ura. 1995. The centenary of the omori formula for a decay law of aftershock activity[J]. J Phys Earth, 43（1）: 1-33.
Wan Y G, Shen Z K. 2010. Static Coulomb stress changes on faults caused by the 2008 Mw 7.9 Wenchuan, China earthquake[J]. Tectonophysics, 491（1/4）: 105-118. doi: 10.1016/j.tecto.2010.03.017.
Wang J J, Xu C J, Freymueller J T, Li Z H. 2017. Probing Coulomb stress triggering effects for a MW＞6.0 earthquake sequence from 1997 to 2014 along the periphery of the Bayan Har block on the Tibetan Plateau[J]. Tectonophysics, 694: 249-267. doi: 10.1016/j.tecto.2016.11.009.
Wang Q, Qiao X J, Lan Q G, Freymueller J, Yang S M, Xu C J, Yang Y L, You X Z, Tan K, Chen G. 2011. Rupture of deep faults in the 2008 Wenchuan earthquake and uplift of the Longmen Shan[J]. Nat Geosci, 4（9）: 634-640.
Wang R J, 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.
Wang Y Z, Wang F, Wang M, Shen Z K, Wan Y G. 2014. Coulomb stress change and evolution induced by the 2008 Wenchuan earthquake and its delayed triggering of the 2013 Mw 6.6 Lushan earthquake[J]. Seismol Res Lett, 85（1）: 52-59.
Xiong X, Shan B, Zheng Y, Wang R J. 2010. Stress transfer and its implication for earthquake hazard on the Kunlun Fault, Tibet[J]. Tectonophysics, 482（1/4）: 216-225. doi: 10.1016/j.tecto.2009.07.020.
Xu X W, Wen X Z, Han Z J, Chen G H, Li C Y, Zheng W J, Zhnag S, Ren Z Q, Xu C, Tan X B, Wei Z Y, Wang M M, Ran J J, He Z T, Liang M J. 2013. Lushan MS7.0 earthquake: a blind reserve-fault event[J]. Chinese Science Bulletin, 58（28/29）: 3437-3443.
Zhang Y, Wang R J, Yun T C, Xu L S, Du F, Jin M P, Tu H W, Dahm T. 2014. Kinematic rupture model and hypocenter relocation of the 2013 Mw 6.6 Lushan earthquake constrained by strong-motion and teleseismic data[J]. Seismol Res Lett, 85（1）: 15-22.
Zhuang J C, Ogata Y, Vere-Jones D. 2002. Stochastic declustering of space-time earthquake occurrences[J]. J Am Stat Assoc, 97（458）: 369-380. doi: 10.1198/016214502760046925.
Zhuang J C, Chang C P, Ogata Y, Chen Y I. 2005. A study on the background and clustering seismicity in the Taiwan region by using point process models[J]. J Geophys Res, 110（B5）: B05S18. doi: 10.1029/2004jb003157.
-
期刊类型引用(0)
其他类型引用(1)
计量
- 文章访问数: 1995
- HTML全文浏览量: 1116
- PDF下载量: 140
- 被引次数: 1
