地球物理学报  2015, Vol. 58 Issue (1): 184-193   PDF    
2014年2月12日Mw6.9于田地震震源破裂过程及对周围断层的应力影响
周云, 王卫民, 熊林, 何建坤    
中国科学院青藏高原研究所, 大陆碰撞与高原隆升重点实验室, 北京 100101
摘要:2014年2月12日,在新疆于田县发生了里氏7.3级地震.在该地震震中附近,前人研究证明发育了大量规模不同的活动断层(如康西瓦断裂与贡嘎错断裂等).根据地震触发理论,地震发生后因地壳同震变形会导致其周边不同性质断裂破裂应力发生变化,进而影响其地震的潜在危险性.本文利用地震远场波形记录,反演了该地震滑动模型.之后,根据弹性无限半空间位错理论,计算了该地震在近场范围内活动断裂上的同震应力变化.其目的在于讨论于田地震引起的附近断裂上的库仑应力变化以及这些活动断裂可能潜在的地震危险性.在地震发生后,从国际地震学联合会(IRIS)地震数据中心,下载了震中距离介于30°~90°的地震远场波形记录,为保证台站方位角分布均匀,从中挑选了27个不同方位角的高信噪比地震记录参与理论地震图的生成和波形反演过程.我们采用广义射线理论计算生成远场理论地震波形数据.每个子断层参数的反演则利用基于全局化反演的快速模拟退火反演方法.在有限断层反演过程中,我们采用了强调波形拟合的相关误差函数作为待反演的目标函数,拟合的断层参数使目标函数为最小.之后,根据弹性无限半空间位错理论,以库仑破裂准则为基础,结合反演得到的地震震源机制解和地震位错模型,计算该地震引起的近场断层面上库仑应力的变化.由远场波形计算结果可以看到,于田地震的震源深度为10 km,地震断层的倾角约71.9°,破裂面上最大的同震位移达到210 cm,以左旋走滑为主并具有正倾滑分量,地震能量主要在前15 s内释放.由此得到该地震的地震矩为2.91×1019N·m,地震震级为Mw6.9.于田地震引发的余震,大致分布在三个区域内:普鲁断裂北部、康西瓦断裂东部和贡嘎错断裂中部.弹性应力计算结果表明,于田地震导致阿尔金断裂西段、普鲁断裂中段、康西瓦断裂东段和贡嘎错断裂中段的静态库仑应力明显增加,其中以康西瓦断裂东段和贡嘎错断裂中段应力增量为最大,分别达到了0.05 MPa和0.04 MPa.大量研究证明,当地震所导致的库仑应力变化大于0.01 MPa时将具有明显的地震触发作用.根据本文结果,2014年于田Mw6.9地震使普鲁断裂、贡嘎错断裂和康西瓦断裂上的库仑应力增量均超过了触发阈值,具有被触发出地震的潜在危险.因此,在以后的地震学研究中,应加强对该三条断裂地震危险性的研究和监测.此外,近6年以来,研究区域发生了3次6级以上的地震.这些地震均沿着贡嘎错断裂,由南西向北东迁移,逐步靠近阿尔金断裂,并且逐渐由正倾滑型地震转变为走滑型地震.阿尔金断裂的走滑速率达到了9 mm·a-1,所以,尽管本次地震导致的阿尔金断裂库仑应力增量小于0.01 MPa,阿尔金的地震危险性也应该加强监测.
关键词2014年于田地震     同震破裂     库仑应力     地震危险性    
Rupture process of 12 February 2014,Yutian Mw6.9 earthquake and stress change on nearby faults
ZHOU Yun, WANG Wei-Min, XIONG Lin, HE Jian-Kun    
Key Laboratory of Continental Collision and Plateau Uplift, Institute of Tibetan Plateau Research, CAS, Beijing 100101, China
Abstract: On February 12 2014, an Ms7.3 earthquake hit Yutian, Xinjiang, the epicenter is at the intersection of Karakax fault and Gonggacuo fault. Earthquake triggering theory indicates that Coulomb stress on nearby faults will change because of crustal coseismic slip after earthquake, and it will affect earthquake potential risk. In this paper, we estimate focal mechanism and rupture process with far field seismic wave data, and calculate coseismic stress change on nearby faults around epicenter. The purpose is to discuss Coulomb stress change and seismic potential hazard on these faults caused by Yutian earthquake. After the earthquake, we download seismic far field wave data of which epicenter is between 30°to 90° from IRIS. We select 27 high SNR(signal to noise ratio) seismic record to make theoretical seismogram. We use generalized ray theory to get synthetic seismic wave map, and sub fault parameter inverse process is based on simulated annealing algorithms. In the way of changing fault parameter to fit actual and synthetic wave form, optimum solution of every sub fault is found. Based on fault rupture inversion model and infinite elastic half space theory, we calculate Coulomb stress change on nearby faults caused by the earthquake.
The inversion results show that source depth of earthquake is 10 km, dip angle is 71.9°, largest coseismic displacement is 210 cm. Seismic moment is 2.91×1019N·m, the main seismic energy is released in former 10 second. Aftershocks are mainly distributed in three regions: north Pulu fault, east Karakax fault and centre Gonggacuo fault. Stress increased significantly on western segment of Altyn fault, central part of Pulu fault, eastern segment of Karakax fault and central segment of Gonggacuo fault. Among them the largest stress change is 0.05 MPa on Karakax fault and 0.04 MPa on Gonggacuo fault. Former research shows that, Coulomb stress change caused by earthquake larger than 0.01 MPa will dramatically increase seismic risk on faults. In our research, stress change on Pulu fault, Gonggacuo fault and Karakax fault exceed triggering threshold, so all these three faults have seismic risk. In the near future seismic potential hazard on these faults should be closely monitored. In the past 6 years, 3 moderately strong earthquakes happen in the study area, epicenters migrate from Gonggacuo fault to Altyn fault, from southwest to northeast. Though stress change on Altyn fault is much lower, in consideration of 9 mm·a-1 slip rate, seismic risk on Altyn fault should cause enough attention.
Key words: Yutian earthquake in 2014     Coseismic rupture     Coulomb stress     Seismic hazard    
1 引言

2014年2月12日,在青藏高原北缘发生里氏7.3级的地震(图 1).该地震发生在高原平均海拔约5000 m的人口稀少地区,地震发生后震中近场地区仅10天内余震就达4000多次.

图 1 于田地震区域构造与地震活动性
黑色粗线条为主要断裂,黑色细线条为小断裂;绿色圆圈为1900年以来6级以上地震;两个蓝色震源机制为2008年M7.3和2012年 M6.2地震(Global CMT); 红色震源机制为我们的反演结果;小圆圈为2014年于田地震引发的余震(余震信息来自中国地震台网).图中断层信息万永革(2010).
Fig. 1 Tectonic structure and seismic activity in Yutian earthquake zone
Black and thin lines denotes major and small faults,respectively. Green circles mean M>6 earthquakes after 1900. Two blue focal mechanisms mean 2008 M7.3 and 2012 M6.2(Global CMT). Red focal mechanism is our result. Small circles represent aftershocks caused by 2014 Yutian earthquakes(from CSN). Faults information is from Wan et al., 2010.

地震应力触发理论认为,当一个地震发生后它将改变其周边活动断裂的应力状态(Deng and Sykes, 1997Lin and Stein, 2004单斌等,2009He and Peltzer, 2010),进而触发或抑制周边活动断裂潜在地震危险性.这一理论的提出,使人们意识到一些大地震之间可能存在内在联系,相互之间存在触发机制.于是,通过计算地震之后周边断层上应力加载卸载情况,可以评估地震潜在危险性,并在一定程度上预测地震可能发生区域.如1992年美国加州Mw7.3兰德斯地震,它使距震中约30 km处圣安德列斯断层上应力值增加了0.3 MPa,触发了3 h后Mw6.3级大熊地震(King et al., 1994).Stein等(1997)计算了土耳其North Anatolian断裂带10次M≥6.7级地震的同震应力变化,认为Izmit地区会有强震发生.1999年M7.4级Izmit地震发生,印证了Stein的预测.Izmit地震导致的断层应力变化,使距震中约110 km处Duzce-Hendek断层地震危险性急剧升高(Barka,1999),三个月后M7.1级Duzce地震发生.2004年12月M9.1级苏门答腊地震发生后,McCloskey等(2005)计算了库仑应力分布,认为North Sumatra断层短期内发生大地震的危险性很高.2005年3月,苏门答腊北部发生M8.6级地震.2008年Mw7.9级汶川地震发生后,距震中 约50 km的雅安断裂上,库仑应力增加了近0.01 MPa(Parsons et al., 2008),导致2013年的Mw6.6级芦山地震的提前发生(单斌等,2013刘成利等,2013).

青藏高原是我国活动断裂分布最广泛的地区之一(邓起东等,2002).由于印度板块不断向欧亚板块俯冲,青藏高原开始对大陆岩石圈挤压,受到塔里木块体的阻挡,青藏高原急剧隆升,并形成了许多复杂 的断裂构造.在2014年于田地震周边,分布有规模 巨大的阿尔金断裂,以及普鲁断裂和康西瓦断裂等. 这些断裂自第四纪以来的滑动速率均在几个mm·a-1到1 cm·a-1间,具有强烈的地震活动性(Li et al., 2012Mériaux et al., 2004Shen et al., 2001).本次地震是继2001年昆仑山口西8.1级地震、2008年于田7.3级地震之后,发生在青藏块体西部边缘的又一次7级以上地震(蒋廷臣等,2012洪顺英等,2010).近年来,在该地区的一些正倾滑活动断层上 已经发生了2008年MS7.3级地震和2012年MS6.4 级地震.历史地震记录表明,该地区具有强烈活动的阿尔金断裂、普鲁断裂和康西瓦断裂等在近百年内一直处于弹性应变的累积状态而没有强地震事件.2014年于田地震是否会影响这些活动断裂潜在的地震危险性?这是值得开展深入研究的课题.为此,本文利用远场体波地震资料反演了该地震的震源机制解和震源破裂过程.之后应用反演的地震滑动模型,根据弹性无限半空间位错理论计算了该地震在近场范围内活动断裂上的同震应力变化.其目的在于讨论于田地震引起的附近断裂上的库仑应力变化及其这些活动断裂可能潜在的地震危险性. 2 震源过程反演

使用有限断层反演方法(姚振兴和纪晨等,1997王卫民等,2008郝金来等,2011),反演了2014年于田地震的震源破裂过程.地震发生后,从国际地震学联合会(IRIS)地震数据中心下载了地震波形记录,采用张光伟等(2014)重定位研究得到的震中位置(82.56°E,36.04°N),来研究地震震源机制和震源破裂过程.由于发震区域海拔很高,缺乏近场资料,在计算中我们只用了远场波形记录.

我们利用远场P波反演获得的这次地震的震源机制解,震源深度10 km.反演得到节面1的走向 240°,倾角71.9°,滑动角-2.2°;节面2的走向330.6°,倾角87.9°,滑动角-177.8°(图 2).节面1的走向、倾角与实际贡嘎断裂最为接近,判定节面1为断层破裂面.破裂过程模型计算的理论地震图与观测波形资料进行拟合,波形相关性较好(图 3).依据点源模型确定的断层深度、走向和倾角,用一条沿走向长85 km,沿倾向宽21 km的断层模拟破裂面,将之剖分成5 km×3 km的119个子断层.使用P波垂向和SH波切向波形记录,反演得到断层面上滑动分布(图 4).

图 2 2014年M7.3于田地震震源机制解
采用下半球投影,同时给出了点源模型的P波垂向位移理论图(红线)与资料(黑线)的拟合情况.图形下方给出了两组节面解(左下,λ,δ,θ,h别表示断层滑动角、断层倾角、断层走向、震源深度)和点源模型的震源时间函数(右下).“+”(绿)和“-”(蓝)圆圈表示P波初动极性.
Fig. 2 Focal mechanism of 2014 M7.3 Yutian earthquake
Lower hemisphere projection is used here. The observed P wave records(black line) and synthetic seismogram(red line)based on the simple point model are compared. The parameters of two possible fault planes are listed also,with λ,δ,θ,h indicating the rake angle,dip angle,strike direction and source depth respectively. The obtained source time function are plotted. Green and blue circles indicate “+” and “-”,respectively.

图 3 远场台站分布(a)和有限断层模型的远场波形拟合(b)
(a)中黑色三角为P波记录的台站,蓝色圆圈为SH波记录的台站,红色圆圈为震中位置.(b)左列为P波垂向位移理论图(红线)与记录(黑线)的对比;右列为SH波切向位移理论图(红线)与记录(黑线)的对比;并给出每个记录的方位角(左上)、震中距(左下)、 振幅比(上)和台站名(右).
Fig. 3 Distribution of teleseismic stations(a),comparison of observed records with synthetic waveforms(b)
Station distribution(a),black triangles show stations of P wave,blue circles show stations of SH wave,red circle is epicenter. Comparison of observed(red) and synthetic waveforms(black)for P vertical components, and SH waveforms(b)are shown in the panel,the azimuth,epicentral distance,amplitude ratio and stations name are also indicated.

图 4 地震断层面上的滑动分布(c),图中分别给出了地震断层滑动分布投影(a)和断层面滑动矢量特征(b)Fig. 4 Inverted slip distribution on the fault(c),the projection of slip distribution on the ground(a) and slip vector character on fault(b)

图 5给出了在空间上于田地震破裂随时间的扩展过程.可以看到,整个破裂持续了18~20 s,在6~8 s时破裂达到了最大值.

结果表明,于田地震是以左旋走滑为主并具有正倾滑分量的地震,破裂在断层面上的分布比较集中.地震矩为2.91×1019N·m,地震震级为Mw=6.9,最大滑动210 cm.

于田地震发生后,不同研究机构迅速给出了各自的反演结果(表 1).可以看到,本研究得到的各个断层参数与USGS和Harvard CMT的结果都较为一致.

图 5 于田地震震源破裂过程(地表投影)Fig. 5 Rupture process snapshots of Yutian earthquake

表 1 不同机构的震源及断层参数Table 1 Source and fault parameter from different research group
3 同震应力计算

根据库仑破裂定律

Δσf=ΔτsΔσn

可以计算断层面上的库仑应力变化(King et al., 1994).式中Δσf为断层面上库仑应力改变量,Δσn为正应力变化,Δτs为剪应力变化,μ为摩擦系数.摩擦系数无法直接观测,一般选取经验值.改变μ值对应力变化值有一定影响,对应力的空间分布影响较小(张国宏等,2008).在本文的计算中,参照Stein(1992)King(1994)的值,取摩擦系数μ=0.4.

基于以上反演得到的震源机制解和地震位错模型(图 4),计算了于田地震产生的静态库仑应力变化.考虑到所研究断层面在不同深度,库仑应力变化量也不一样,而触发作用一般取决于最大库仑应力增量(Lin and Stein, 2004单斌等,2009).在计算时,我们以5 km为间隔,计算了0~20 km的库仑应力分布,选择其中的最大值.计算过程中,相关断层参数如表 2所示.

表 2 研究区域内的断层信息(万永革等,2010)Table 2 Fault information of earthquake zone

图 6给出了于田地震在邻近断层(阿尔金断裂、普鲁断裂、贡嘎错断裂、康西瓦断裂和龙木错—邦达错断裂)上的库仑应力变化和余震分布情况.余震主要分布在三个区域内:普鲁断裂北部(塔里木块体内)、康西瓦断裂东部和贡嘎错断裂中部.可以看出,绝大多数余震位于断层应力增强区周边.距离于田地震震中较近的断层,库仑应力变化很大,发生的余 震也多;而较远处断层上的库仑应力则很快衰减,余震一般也较少.

图 6 于田地震引起的周边区域(a)及主要断层上(b)的库仑应力变化
(a)于田地震在周边区域导致的应力变化,小圆圈表示余震(余震信息来自中国地震台网),灰色区域为断层破裂面在地表投影;(b)于田地震在周边断层上引起的应力变化.红色五角星表示2014年于田地震,两个蓝色震源球表示2008年M7.3和2012年M6.2.
Fig. 6 Coulomb stress change in nearby zone and major faults caused by Yutian earthquake
Stress change caused by Yutian earthquake in nearby area(a),small circles denote aftershocks(information is from CSN),the gray rectangular region is projection of fault rupture on surface. Stress change on major faults is showed in Fig.b,red star denotes 2014 Yutian earthquake,two blue focal mechanisms indicate 2008 M7.3 and 2012 M6.2.

阿尔金断裂位于震中东北方向,距离震中超过了100 km,库仑应力增量不超过5 kPa,低于一般 意义上的应力触发值0.01 MPa(Harris and Simpson, 1998). 然而,当应力增加低于0.01 MPa时,也可能触发地震(Ziv and Rubin, 2000).所以阿尔金断裂的地震危险性,还需要继续研究.余震在阿尔金断裂上分布很少.

普鲁断裂东部应力增强,最大值达到了0.015 MPa; 断层中部达到0.02 MPa以上,有助于断层发生破裂;断层西部应力变化不大.从余震分布来看,断层北部为一个余震密集区.相对于另外两个区域,该区域与主震震中位置较远,余震如此集中比较反常.然而,在计算过程中,我们发现,于田地震在周边区域,产生的应力增强区,有两个位于破裂断层北部和南部.发生在塔里木块体内的余震,80%以上发生在北部的应力增强区内,与应力分布有很好的相关性.所 以,该余震集中区域,可以用地震静态触发理论解释.

贡嘎错断裂由于靠近于田地震震中位置,应力变化很复杂,主要以增加为主.沿着断层走向,断层 北缘应力变化不大,之后迅速增大,达到了0.04 MPa;然后是应力影区,卸载了0.05 MPa;之后应力增加,达到0.02 MPa;最后在断层西南部,由于离震中较远,应力增量逐渐减小.断层面附近的余震一部分位于应力增强区内,在2012年Mw6.2地震的震中位置,可能与该地震相关;另一部分位于应力减弱区 内,可能与地震动态触发效应有关(Gonzalez-Huizar and Velasco, 2011).

康西瓦断裂东部库仑应力以增加为主,最大值超过了0.05 MPa,大大增强了断层上发生地震的可能性;断裂西部库仑应力变化不大.此处也有一个余震聚集区,位于库仑应力增加的区域.

龙木错—邦达错断裂由于距离震中较远,库仑应力变化不大,断层中部应力稍微增加,该断裂周边,基本上没有余震的发生.

4 结论与讨论

通过远场地震波形反演,结合弹性无限半空间位错理论的静态库仑应力计算,就2014年2月12日于田地震有以下主要结论:

(1)于田地震为震级Mw6.9,震源深度10 km,以左旋走滑为主并具有正倾滑分量的地震,破裂在断层面上的分布比较集中,最大滑动210 cm.地震能量主要在前15 s内释放.断层走向240°,倾角72°,地震矩为2.91×1019N·m.

(2)于田地震显著增加了普鲁断裂中部、康西瓦断裂东部、贡嘎错断裂北部和中部的应力积累,大大提升了这些断层上地震发生的概率.其中以康西瓦断裂东部和贡嘎错断裂中部应力增量为最大,分别达到了0.05 MPa和0.04 MPa.

(3)于田地震引发的余震,大致分布在三个区域内:普鲁断裂北部、康西瓦断裂东部和贡嘎错断裂中部.前两个区域均为应力增强区,与于田地震破裂所导致的静态库仑应力触发有关.

Stein等(1994)Chan和Stein(2009)认为,当断层上库仑应力增加超过0.01 MPa,就可能会触发地震.在本事例中,绝大多数断裂已经达到了这个值.2008年于田地震之后,万永革等(2010)研究发现,于田地震普遍增强了附近大断裂上的库仑应力值.两次地震增加的库仑应力,对周边断裂的地震 已经具有了触发意义.这些断裂大多是活动断裂,走 滑速率很大.康西瓦断裂的活动速率为7±3 mm·a-1,阿尔金的活动速率可以达到9 mm·a-1(Mériaux et al., 2004He et al., 2013).活动断裂的特征地震具有一定的旋回性(Shimazaki and Nakata, 1980Segall and Harris, 1987),自1924年M7.1级地震(图 1),至2007年的近百年时间里,该区域没有发生过7级以上的大地震,这些断裂一直处于弹性应变的累积状态.增加0.01~0.05 MPa的库仑应力,有利于使这些断裂的某些段落接近破裂邻近状态所需时间提前.因此于田地震后本区的潜在地震危险性研究要引起重视.图 6显示,有部分余震位于断层应力减弱区内,这可能有以下几个原因:(1)余震的震源机制尚不清楚.计算中我们选择主震破裂面的走向、倾向和倾角作为参数,如果考虑余震震源机制,结果可能会有变化.(2)动态触发的作用没有考虑.此例中我们只考虑静态应力的因素.(3)地壳中流体作用也会对余震分布有影响.所以,于田地震余震的分布,还有待深入研究.

上述研究表明,2014Mw6.9级地震是发生在阿尔金断裂西段走向发生明显变化的贡嘎错断裂上(图 6).这一空间关系与2001年昆仑山Mw7.8大地震有相似性,该地震首先在其西端太阳湖次级断裂上发生走滑破裂(Lasserre et al., 2005Robinson et al., 2006),之后以超剪切速度破裂沿东昆仑断裂向东扩展达300余公里.2014 Mw6.9走滑为主的地震至少使阿尔金大型走滑断裂在其100 km长的段落上增加了5 kPa的库仑应力.该应力显然小于前人认为的地震触发阈值(Stein et al., 1994),但数值模拟结果显示(He and Chéry,2008),阿尔金断裂的强度可能很小(摩擦系数小于0.1),因此2014 Mw6.9级于田地震后阿尔金断裂的地震潜在危险性必须引起足够的重视.

2008年以来,该区已经发生了2008 Mw7.3、2012 Mw6.2和2014 Mw6.9级地震.从空间上看,这些地震基本上是沿着贡嘎错断裂,从南西向北东迁移.前两个地震以正倾滑型为主,到了最后一个则是以走滑型为主(图 1).这些地震的空间距离在50~100 km之间,那么,根据地震触发理论,这些地震是否是一个早期地震导致的库仑应力变化后的连续触发?在触发过程中为什么会有这样的时间跨度?这是否意味着有黏-弹性应力转变(Savage,2000Johnson et al., 2005Jia et al., 2012)?这些问题有待于我们以后进一步深入研究.

参考文献
[1] Barka A. 1999. The 17 august 1999 Izmit earthquake.Science, 285(5435): 1858-1859.
[2] Chan C H, Stein R S. 2009. Stress evolution following the 1999 Chi-Chi, Taiwan, earthquake: consequences for afterslip, relaxation, aftershocks and departures from Omori decay. Geophysical Journal International, 177(1): 179-192.
[3] Deng J S, Sykes L R. 1997. Stress evolution in southern California and triggering of moderate-, small-, and micro-size earthquakes. Journal of Geophysical Research: Solid Earth (1978—2012), 102(B11): 24411-24435.
[4] Deng Q D, Zhang P Z, Ran Y K, et al. 2002. Basic feature of China active structures. Science in China (D Series) (in Chinese), 32(12): 1020-1030.
[5] Gonzalez-Huizar H, Velasco A A. 2011. Dynamic triggering: Stress modeling and a case study. Journal of Geophysical Research: Solid Earth (1978—2012), 116(B2), doi: 10.1029/2009JB007000.
[6] Hao J L, Wang W M, Yao Z X, 2011. Source process of the 2011 Mw9.0 Tohuko Japan earthquake. Science China Earth Science (in Chinese), 41(6): 745-749, doi: 10.1007/s11430-011-4241-y.
[7] Harris R A, Simpson R W. 1998. Suppression of large earthquakes by stress shadows: A comparison of Coulomb and rate-and-state failure. Journal of Geophysical Research: Solid Earth (1978—2012), 103(B10): 24439-24451.
[8] He J,ChéryJ. 2008.Slip rates of the AltynTagh, the Kunlun and the Karakorum faults (Tibet) from 3D mechanical modeling. Earth and Planetary Science Letters, 274(1-2): 50-58.
[9] He J,PeltzerG. 2010.Poroelastic triggering in the 9-22 January 2008 Nima-Gaize (Tibet) earthquake sequence. Geology, 38(10): 907-910.
[10] He J K, VernantP, ChéryJ, et al. 2013. Nailing down the slip rate of the AltynTagh fault. Geophysical Research Letters, 40(20): 5382-5386.
[11] Hong S Y, Shen X H, Shan X J, et al. The calculation and analysis of the co-seismic deformation field of Yutian Ms7.3 earthquake basing on the ascending and descending orbit ASAR data. Remote Sensing for Land & Resources (in Chinese), 2010,84(1):98-102.
[12] Jia K, Zhou S, Wang R. 2012. Stress interactions within the strong earthquake sequence from 2001 to 2010 in the Bayankala block of eastern Tibet. Bulletin of the Seismological Society of America, 102(5): 2157-2164.
[13] Jiang Tingchen, et al. Deformation analysis of Yutianearthquake based on ScanSARinterferometry. Journal of Geodesy and Geodynamics, 2012, (6):13-16
[14] Johnson K M, Segall P, Yu S B. 2005. A viscoelastic earthquake cycle model for Taiwan. Journal of Geophysical Research: Solid Earth (1978—2012), 110(B10), doi: 10.1029/2004JB003516.
[15] 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.
[16] Lasserre C, PeltzerG, CrampéF, et al. 2005.Coseismic deformation of the 2001 Mw=7.8 Kokoxiliearthquake in Tibet, measured by synthetic aperture radar interferometry. J. Geophys. Res., 110: B12408, doi: 10.1029/2004JB003500.
[17] Li H B, Van der Woerd J, Sun Z M, et al. 2012. Co-seismic and cumulative offsets of the recent earthquakes along the Karakax left-lateral strike-slip fault in western Tibet. Gondwana Research, 21(1): 64-87.
[18] Lin J, Stein R S. 2004. Stress triggering in thrust and subduction earthquakes and stress interaction between the southern San Andreas and nearby thrust and strike-slip faults. Journal of Geophysical Research: Solid Earth (1978—2012), 109(B2), doi: 10.1029/2003JB002607.
[19] Liu C L, Zheng Y, Ge C, et al. 2013.Rupture process of the M7.0 Lushan earthquake. Science in China (Ser D), 43(6): 1020-1026.
[20] MériauxA S, Ryerson F J, Tapponnier P, et al. 2004. Rapid slip along the central AltynTagh Fault: Morphochronologic evidence from Cherchen He and SulamuTagh. J. Geophys. Res, 109(B6): B06401, doi: 10.1029/2003JB002558.
[21] McCloskey J, Nalbant S S, Steacy S. 2005. Indonesian earthquake: Earthquake risk from co-seismic stress. Nature, 434(7031): 291-291.
[22] Parsons T, Ji C, Kirby E. 2008. Stress changes from the 2008 Wenchuan earthquake and increased hazard in the Sichuan basin. Nature, 454(7203): 509-510.
[23] Robinson D P, Brough C, Das S. 2006. The Mw7.8, 2001 Kunlunshan earthquake: Extreme rupture speed variabilityand effect of fault geometry. J. Geophys. Res., 111(B8): B08303, doi: 10.1029/2005JB004137.
[24] Savage J C. 2000.Viscoelastic-coupling model for the earthquake cycle driven from below. Journal of Geophysical Research: Solid Earth (1978—2012), 105(B11): 25525-25532.
[25] Segall P, Harris R. 1987. Earthquake deformation cycle on the San Andreas fault near Parkfield, California. Journal of Geophysical Research: Solid Earth (1978—2012), 92(B10): 10511-10525.
[26] Shan B, Xiong X, Zheng Y, et al. Stress changes on major faults caused by Mw7.9 Wenchuan earthquake, May 12, 2008. Sci. China Earth Sci., 2009, 52(5):593-601.
[27] Shan B, Xiong X, Zheng Y, et al. Stress changes on major faults caused by 2013 Lushan earthquake, and its relationship with 2008 Wenchuan earthquake. Science China (Earth Sciences), 2013, 56(7): 1169-1176, doi: 10. 1007/s11430-013-4642-1.
[28] Shen Z K, Wang M, Li Y X, et al. 2001. Crustal deformation along the AltynTagh fault system, western China, from GPS. Journal of Geophysical Research: Solid Earth (1978—2012), 106(B12): 30607-30621.
[29] Shimazaki K, Nakata T. 1980. Time-predictable recurrence model for large earthquakes. Geophysical Research Letters, 7(4): 279-282.
[30] 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.
[31] Stein R S, King G C P, Lin J. 1994. Stress triggering of the 1994 M=6.7 Northridge, California, earthquake by its predecessors.Science, 265(5177): 1432-1435.
[32] Stein R S, Barka A A, Dieterich J H. 1997. Progressive failure on the North Anatolian fault since 1939 by earthquake stress triggering. Geophysical Journal International, 128(3): 594-604.
[33] Wan Y G, Shen Z K, Sheng S Z, et al. 2010. The mechanical effects of the 2008 M7.3 Yutian, Xinjiang earthquake on the neighboring faults and itstectonic origin of normal faulting mechanism. ChineseJournal of Geophysics, 53(2): 280-289. doi:10.3969/j.issn.000i-5733.2010.02.006.
[34] Wang W M, Zhao L F, Li J, et al. 2008. Rupture process of the Ms8.0 Wenchuan earthquake of Sichuan, China. Chinese J. Geophys. (in Chinese), 51(5): 1403-1410.
[35] Yao Z X, Ji C.The inverse problem of finite fault study in time domain. Chinese J. Geophysics. (in Chinese), 1997,40(5):691-701.
[36] Zhang G H, Shan X J, Li W D. The Coulomb failure stress change associated with the Ms8.0 Wenchuan earthquake and the risk prediction of its surrounding faults. Seismicity and Geology (in Chinese), 2008, 30(4):935-944.
[37] Zhang G W, Lei J S, Sun C Q. 2013. Relocation of the 12 February 2014 Yutian, Xinjiang, mainshock (Ms7.3) and its aftershock sequence. Chinese J. Geophys. (in Chinese), 57(3):1012-1020, doi:10.6038/cjg20140330.
[38] Ziv A, Rubin A M. 2000. Static stress transfer and earthquake triggering: No lower threshold in sight?. Journal of Geophysical Research: Solid Earth (1978—2012), 105(B6): 13631-13642.
[39] 邓起东, 张培震, 冉勇康等. 2002. 中国活动构造基本特征. 中国科学(D辑), 32(12): 1020-1030.
[40] 郝金来, 王卫民, 姚振兴. 2011. 2011年3月11日日本东北沿海Mw9.0级地震震源破裂过程. 中国科学: 地球科学, 41(6): 745-749.
[41] 洪顺英, 申旭辉, 单新建等. 2010. 基于升降轨ASAR的于田Ms7.3级地震同震形变场信息提取与分析. 国土资源遥感, (4): 98-102.
[42] 蒋廷臣, 李陶, 刘经南. 2012. 基于宽幅SAR干涉测量的新疆于田7.3级地震形变分析. 大地测量与地球动力学, 32(6): 13-16.
[43] 刘成利, 郑勇, 葛粲等. 2013. 2013年芦山7.0级地震的动态破裂过程. 中国科学: 地球科学, 43(6): 1020-1026
[44] 单斌, 熊熊, 郑勇等. 2009. 2008年5月12日Mw7.9汶川地震导致的周边断层应力变化. 中国科学D辑: 地球科学, 39(5): 537-545.
[45] 单斌, 熊熊, 郑勇等. 2013. 2013年芦山地震导致的周边断层应力变化及其与2008年汶川地震的关系. 中国科学: 地球科学, 43(6): 1002-1009.
[46] 万永革, 沈正康, 盛书中等. 2010. 2008年新疆于田7.3级地震对周围断层的影响及其正断层机制的区域构造解释. 地球物理学报, 53(2): 280-289, doi: 10.3969/j.issn.0001-5733.2010.02.006.
[47] 王卫民, 赵连锋, 李娟等. 2008. 四川汶川8.0级地震震源过程. 地球物理学报, 51(5): 1403-1410, doi: 10.3321/j.issn:0001-5733.2008.05.013
[48] 姚振兴, 纪晨. 1997. 时间域内有限地震断层的反演问题. 地球物理学报, 40(5): 691-701.
[49] 张国宏, 单新建, 李卫东. 2008. 汶川Ms8.0地震库仑破裂应力变化及断层危险性初步研究. 地震地质, 30(4): 935-944.
[50] 张光伟, 雷建设, 孙长青. 2014. 2014年2月12日新疆于田Ms7.3级地震主震及余震序列重定位研究. 地球物理学报, 57(3): 1012-1020, doi: 10.6038/cjg20140330.