Crustal stress field in the focal area of MS7.0 Jiuzhaigou,Sichuan,earthquake by shear wave splitting
-
摘要: 本文讨论了2017年8月8日九寨沟MS7.0 (MW6.5)地震之后震源区的地震活动性,并针对震后九寨沟地区架设的4个流动地震台站记录到的地震波形数据进行了剪切波分裂参数的计算,分析了震源区地壳应力场的特征性变化。九寨沟MS7.0地震的地震活动性表现为:8月8日主震后震源区地震活动明显增强,直至8月22日之后才有所下降,余震震级整体也呈降低趋势,震级多集中在MS2.0以下。震源区4个流动台站的剪切波分裂结果表明:快波偏振方向近似为东西向,揭示了九寨沟地区的区域应力场方向为近似东西向;区域内慢波延迟时间的平均值为5.04 ms/km,慢波延迟时间在九寨沟MS7.0地震之后降低,反映了随着余震的发生,地壳应力不断调整,应力持续释放,直至8月底地壳应力场逐渐趋于稳定状态。
-
关键词:
- 地震活动性 /
- 剪切波分裂 /
- 地壳应力 /
- 九寨沟MS7.0地震
Abstract: This paper discussed the seismic activity in the focal area after Jiuzhaigou MS7.0 (MW6.5) earthquake on August 8, 2017, and calculated the shear wave splitting parameters by using the seismic waveform data from four temporary seismic stations in Jiuzhaigou area, then analyzed the characteristic change of crustal stress in the focal area. After the MS7.0 Jiuzhaigou earthquake, the seismic activity was significantly enhanced in the focal area on August 8, 2017. Until after 22 August 2017, the seismic activity dropped slightly; correspondingly, the magnitude of the aftershocks also became smaller, concentrating below MS2.0. The results of the polarization directions of fast shear wave at four seismic stations in the focal area show that the polarization directions are approximately east-west, which reveals the direction of the regional stress field in Jiuzhaigou area is approximately east-west. Meanwhile, the average value of the time delay of slow shear wave in the studied area is 5.04 ms/km, and after Jiuzhaigou earthquake the time delay of slow shear wave decreased, which reflects the release and adjustment of the crustal stress with the aftershocks occurrence. And thus the stress field in the crust tends to be stable at the end of August. -
引言
地震波的震相在地震学发展及地球结构研究中起着重要的作用。随着数值计算方法的发展和观测技术水平的提高,逐渐可以利用地震波形中尽可能多的信号来研究地球结构(崇加军,2013)。地震波是源自震源和震中区的弹性波,震相是地震波显示在地震记录图上的信号(赵荣国,1999)。震相分析对深入研究地球的内部构造和探查地震活动规律具有极为重要的意义,地震震相识别也是省级地震分析和编目工作的一项常规工作。地震震相识别中,人工识别的成分占较大比例,而且由于震源机制和传播介质的复杂性,同一种震相在不同地震中的形态均不相同(杨配新等,2004),因此常见的震相Pn,Pg,Sn,Sg较易识别,但非常见的震相识别难度较大。而康拉德面是不连续面,地震波速度在深度变化的过程中均匀变化,因此Pb震相不易被发现,难以得到地震分析人员的关注。
1923年,奥地利研究人员康拉德发现,大陆地壳内部还存在这样一个界面(Conrad,1925):该界面以上,岩石密度较低,为硅铝层,由中酸性岩石组成;该界面以下,岩石密度较高,为硅镁层(铁镁层),由基性岩石组成。该界面被称为康拉德(Conrad)不连续面,简称康氏面或康拉德面(刘瑞丰等,2014)。此后许多地震学研究人员在大陆地壳内不同区域和不同尺度上也探测到了该界面(Berry,Fuchs,1973;胡德昭等,1989;郭杰等,2013),并对其进行了相关研究。胡德昭等(1989)在中国东南部地壳内探查到康拉德面;1970—1994年期间,苏联的超深钻井项目研究显示康拉德界面上的地震波速度跃变非常大(Richard,1989;Pavlenkova,1993)。但康拉德界面是否普遍存在仍有争议。相关研究认为,康拉德界面与研究区的地壳结构及其演化历史密切相关,在海洋地壳内或者有时在构造活动比较强烈的大陆地壳内,康拉德界面可能是缺失的(Litak,Brown,1989;Richard,1989;Pavlenkova,1993)。近期研究已在有些地质构造活动比较强烈的区域记录到的地震波形中发现有来自康拉德界面的强振幅震相(刘赛君等,2011;郭杰等,2013),且焦煜媛等(2017)已在青藏高原东北缘找到康拉德面存在的直接地震学证据。郭杰等(2013)在位于豫北东的濮阳市地震台的观测中发现了康拉德界面的反射波震相。
关于海南地区Pb震相的研究,也有一些相关的研究成果。范玉兰等(1990)在研究华南地区近震走时表时,描述到海口台、文昌台、定安台、那大台、琼中台等海南省台站在震中距大于151 km时,可记录到清晰的Pb震相;刘赛君等(2011)在海南岛西南海域地壳剖面海陆联合探测研究中提到Pb震相在50—130 km偏移距范围内可被追踪到。
这样看来,有些地区不仅有康拉德面存在的证据,而且前人已记录到海南地区存在Pb震相,因此有必要进一步研究海南地区的Pb震相,故本文一方面以云南省地震台网和广东省地震台网记录到的Pb震相为基础,采用折合走时分析,将其结果作为参考数据;另一方面重点分析海南省地震台网记录到的地震事件,并使用PTD方法(朱元清等,1990)重新定位地震事件的震源深度,以Pb震相理论到时结合实际波形,标注Pb震相,进而拟合标注的震相速度(林建民等,2008)。 此外,综合前人的研究成果(李志雄等,2008;刘赛君等,2011;黄海波等,2012;Kumar et al,2016 )来分析Pb震相的物理特性,最后确定Pb震相。通过识别Pb震相,探查海南地区康拉德面的深度变化,对该地区一维速度模型的建立、地震定位精度的提高以及地震活动、地震定位深度等研究均具有重要意义。
1. 资料和方法
1.1 资料选取
本文研究区域为海南地区(106°E—114°E,17°N—23°N),选取2009年1月1日至2016年8月25日期间发生在研究区域内的地震,并按照以下标准严格筛选地震事件:① 至少被10个台站记录到;② 地震震级ML≥2.0;③ 地震事件含有Pn震相。最终统计筛选得到Pg震相1 549个,Pn震相588个,无Pb震相。由于研究区域的数据中不包含Pb震相,故以2009年至2016年云南省、广东省地震台网的数据为参考资料,按照上述标准,最终云南地区的震相筛选结果为:Pg震相49 521个,Pb震相3 406个,Pn震相10 781个;广东地区的震相筛选结果为:Pg震相5 038个,Pb震相23个,Pn震相2 240个。
1.2 方法
1.2.1 理论走时
由于近震波的传播路程短,受地球曲率影响小,因此在研究近震问题时,通常将地球表面及各层界面看作水平界面(刘瑞丰等,2014)。Pb震相是在康拉德面上的纵波性首波,其传播路径如图1所示。由该图所示的Pb震相走时路径,依据斯奈尔(Snell)定理和地震波走时方程,可求得Pb的理论走时为
${t_{{\rm{Pb}}}} {\text{=}} \frac{\varDelta }{{{v_2}}} {\text{+}} \left( {2{H_1} - h} \right)\frac{{\cos i}}{{{v_1}}},$
(1) 式中,tPb为Pb震相走时,Δ为震中距,H1为康拉德面深度,h为震源深度,v1为一维速度模型中第一层的Pg波平均速度,v2为第二层的Pb波速度,i为Pb波入射角。
1.2.2 折合走时
折合走时与介质厚度(康拉德界面或者莫霍面深度)、震源深度和波速相关(王莉婵等,2016),在同一震源深度的情况下,理论折合走时应为常数值b,震相走时可理解为横向走时和折合走时,因此将式(1)变换为
${t_{\rm redu}} {\text{=}} t - \frac{\varDelta }{v},$
(2) 式中tredu为折合走时,t为理论震相或观测震相走时,v为波速。
1.2.3 PTD定位方法
PTD定位方法是采用初至波为Pg直达波的台站到时和初至波为Pn首波的台站到时,经转换后的到时差来确定地震的震源深度(朱元清等,1990)。测定震源深度的分辨率为深度每改变5 km,到时差改变0.7 s,换句话说,初动到时测量误差每增加0.1 s,深度误差则增加0.7 km (朱元清等,1997)。
本文可辨别出Pb震相的地震共计52次,由于原震源深度定位偏浅,8 km左右的深度占大部分,因此重新测定震源深度,便于下一步震相理论到时的计算。图2给出了使用PTD定位前后的震源深度比较,可见使用PTD方法测定深度可减小初动到时测量误差引起的深度误差,在深度测定方面有较大的优势。
1.2.4 模型选取
数据处理涉及理论值计算,速度模型是必不可少的部分,虽然目前大多采用IASP91模型,但由于本文涉及云南省(2015模型)、广东省(华南模型)和海南省(2015模型)的数据(朱元清等,2017),考虑到一维速度模型与到时数据之间的一致性(Wang,2014),本文采用各省的一维速度模型作为各自省份参与计算的理论到时,各省的速度模型列于表1。
表 1 云南、广东和海南各省份的地壳速度模型表Table 1. The crustal velocity model of Yunnan,Guangdong and Hainan Provinces省份 v1/(km·s−1) v2/(km·s−1) v3/(km·s−1) H1/km H/km 云南 6.01 6.60 7.89 20 41 广东 6.00 6.87 7.96 22 33 海南 6.00 6.84 7.97 21 30 注:v1和v2分别表示第一、二层的平均速度,v3表示莫霍面的平均速度,H1表示第一层的平均厚度,H表示莫霍面的平均厚度。 2. Pb震相识别
2.1 云南、广东Pb震相分析
采用折合走时分析云南、广东两省份的数据,结果显示:云南省的Pb震相数据基本处于理论震相范围内,标注的Pb震相接近该地区的康拉德面;广东省的Pb震相数据几乎全在理论震相范围内,标注的Pb震相在康拉德面上。这表明结果良好,与预期结果相符,故可用此方法分析海南的Pb震相数据。
2.2 海南Pb震相识别
震相识别的步骤为:① 对资料进行预处理,筛选出含有Pn震相的地震事件,共计281次;② 对281次地震事件使用PTD法重新定位其震源深度;③ 计算震相理论到时,结合人工分析识别出Pb震相;④ 对识别出的Pb震相进行速度拟合;⑤ 分析和佐证识别的Pb震相结果。
以近几年海南地区发生的影响较大的两次地震为例:事件1为2012年11月5日19时51分海南万宁ML4.1地震,那大台记录的波形见图3左;事件2为2014年7月28日10时10分海南儋州ML3.3地震,万宁台记录的波形见图3右。
那大台的震中距为169 km,PTD重新定位深度为10 km,计算其Pg,Pb,Pn的理论走时分别为28.21,27.26,26.06 s,在实际波形识别的相应结果分别为28.54,27.33,26.33 s。依据波形特性、振幅、周期等因素,Pb震相处于Pn震相与Pg震相之间,其振幅稍大于Pn,但小于Pg振幅,符合绕射波的动力学特征。实际人工识别Pb震相走时t2与理论走时t1相差无几,人工识别的Pb震相标注见图3。
以同样的方法,对万宁台记录到的事件2进行PTD定位计算,结果显示重定位深度为13 km,Pg,Pb,Pn震相的理论走时分别为21.60,21.18,20.71 s,实际波形识别的相应结果分别为21.61,21.22,21.04 s。
依据上述震相识别步骤,在281次地震事件中,从52次地震事件能识别出Pb震相57个,含有Pb震相的地震分布如图4所示。通过分析识别结果可知,能记录到Pb震相的台站共计16个,占全部台站个数(24)的67%,且分布均匀,统计台站记录到Pb震相的震中距处于66.75—234 km范围内,这与海南岛西南海域地壳剖面海陆联合探测研究中追踪到的Pb震相的震中距结果(刘赛君等,2011)基本相近。
采用速度线性拟合(图5)和折合走时(图6)两种方式分析Pb震相的速度,以佐证识别的震相是否与其物理意义相符,图6中折合走时理论值是基于海南模型计算而得。从图5可知,Pb速度为6.47 km/s,在对数据进行折合走时分析的过程中,扰动第二层的速度为6.65 km/s时,能达到图6的效果,这说明实际震相处于理论值范围内。海南地区在区域构造的伸展作用下,其地壳厚度相对正常大陆型地壳较薄,具有西南厚、东北薄的特点(黄海波等,2012),康拉德面作为上下地壳的分界面,穿过界面的波速由6.0 km/s激增至6.4 km/s (刘赛君等,2011)。
使用折合走时分析全部Pb震相,结果表明,实际震相基本处于理论震相范围内,尽管有些震相的实际走时与理论走时有偏差,这可能与一维模型单一、不能有效地反映局部地壳内部结构的复杂性和不均匀性有关。综合上述速度拟合分析,识别的Pb震相符合海南地区地震波在康拉德面滑行的物理意义。
3. 海南地区康拉德面
康拉德面是地球内部的次级不连续面,其深度介于10—40 km之间,在陆壳内的平均深度约为20 km (殷伟伟等,2017)。考虑到海南地区的地壳较薄,选取10—30 km作为反演海南地区的康拉德面深度Hi。当Hi以一定步长从10 km依次增加至30 km时,将记录到Pb震相的台站震中距、地震震源深度代入式(1),可得到理论走时t1,t1与实测走时t2之差的绝对值的最小值,即Pb震相的走时残差,为Y=min|t1−t2|,则57次地震的平均走时残差Si为
${S_i} = \frac{1}{{57}}\sum\limits_{i = 1}^{57} {{Y_i}} .$
(3) 为了较好地反映康拉德面的深度和速度情况,通过计算康拉德面的速度和深度变化,可得到平均走时残差Si随深度和速度变化的分布。取Pb波速度为6.40—6.84 km/s,以0.1 km/s的步长计算所有走时残差;深度范围取10—30 km,以1 km的步长计算。计算结果表明速度为6.60—6.80 km/s、深度为19—22 km时残差最小。基于该结果,按上述方法重复,进一步细算速度和深度变化,速度以0.02 km/s为步长,深度以0.5 km为步长,计算结果如图7所示。结果表明,海南地区的Pb波速度为6.60—6.72 km/s、深度为19—21.5 km较为合理。
4. 讨论与结论
本文利用海南省记录到的地震事件,识别了Pb震相并反演得到海南地区的康拉德面深度。结果显示,海南地区的Pb波速度介于6.60—6.72 km/s之间、康拉德面在19—21.5 km左右较合理。海南地区Pb波速度与前人的研究结果(范玉兰等,1990;刘赛君等,2011)基本吻合,首次在海南地区得出康拉德面的合理范围。由于初至Pb震相难识别,缺乏实际震相支持,今后会持续关注这方面的震相数据。
本文结果能对地震分析识别Pb震相起到辅助作用,可为建立海南地区的地壳速度模型提供参考资料。但是,由于本文研究的大部分震源位于上地壳,研究结果具有一定的局限性,震源位于下地壳时震相的识别有待进一步研究。
上海市地震局的朱元清研究员和海南省地震局的李志雄研究员、张慧高工对本研究给予了指导,江苏省地震局廖发军高工分别为本文提供了PTD软件,审稿专家提出了修改意见,作者在此一并表示衷心的感谢!
-
图 3 L5112台站记录的地震数据剪切波分裂分析
所记录事件的发震时间为2017年8月10日20时12分59秒;震中位于(33.31°N,103.74°N),震源深度为7 km;震级为MS1.2;震中距为6.5 km。 (a) 三分量地震波行记录,两条红色竖虚线是截取的用于质点偏振图分析的剪切波段;(b) 质点偏振图,箭头标示了快慢波的质点运动方向,红点标记了快慢波的到达
Figure 3. The shear wave splitting analyses of the earthquakes recorded by station L5112
The earthquake with MS1.2 occurred at 20:12:59 on 10 August 2017,the epicenter is located at (33.31°N,103.74°E) with focal depth 7 km and epicentral distance 6.5 km。 (a) Three-component seismic waveforms,where the two red vertical dashed lines delineate the scope of shear wave shown in polarization figure;(b) Particle trail where the arrows are their particle motion directions of fast and slow waves,and red dots mark their arrivals
表 1 九寨沟地区4个流动台站的剪切波分裂参数
Table 1 Parameters of shear wave splitting at the four portable seismic stations in Jiuzhaigou area
台站 窗口内波形个数 平均偏振方向/° 偏振方向误差/° 平均慢波延时
/(ms·km−1)慢波延时误差
/(ms·km−1)L5110 4 104 6.9 4.23 1.63 L5111 217 97 26.7 5.35 2.35 L5112 304 99 23.8 5.08 2.18 L5113 58 82 19.5 4.76 2.14 -
常利军,丁志峰,王椿镛. 2015. 2013年芦山MS7.0地震震源区横波分裂的变化特征[J]. 中国科学:地球科学,45(2):161–168 Chang L J,Ding Z F,Wang C Y. 2014. Variations of shear wave splitting in the 2013 Lushan MS7.0 earthquake region[J]. Science China Earth Sciences,57(9):2045–2052 doi: 10.1007/s11430-014-4866-8
邓起东,高翔,陈桂华,杨虎. 2010. 青藏高原昆仑—汶川地震系列与巴颜喀喇断块的最新活动[J]. 地学前缘,17(5):163–178 Deng Q D,Gao X,Chen G H,Yang H. 2010. Recent tectonic activity of Bayankala fault-block and the Kunlun-Wenchuan earthquake series of the Tibetan Plateau[J]. Earth Science Frontiers,17(5):163–178 (in Chinese)
邓起东,程绍平,马冀,杜鹏. 2014. 青藏高原地震活动特征及当前地震活动形势[J]. 地球物理学报,57(7):2025–2042 Deng Q D,Cheng S P,Ma J,Du P. 2014. Seismic activities and earthquake potential in the Tibetan Plateau[J]. Chinese Journal of Geophysics,57(7):2025–2042 (in Chinese)
高原,吴晶. 2008. 利用剪切波各向异性推断地壳主压应力场: 以首都圈地区为例[J]. 科学通报,53(23):2933–2939 doi: 10.3321/j.issn:0023-074X.2008.23.015 Gao Y,Wu J. 2008. Compressive stress field in the crust deduced from shear-wave anisotropy: An example in capital area of China[J]. Chinese Science Bulletin,53(18):2840–2848
高原,吴晶,易桂喜,石玉涛. 2010. 从壳幔地震各向异性初探华北地区壳幔耦合关系[J]. 科学通报,55(29):2837–2843 Gao Y,Wu J,Yi G X,Shi Y T. 2010. Crust-mantle coupling in North China: Preliminary analysis from seismic anisotropy[J]. Chinese Science Bulletin,55(31):3599–3605 doi: 10.1007/s11434-010-4135-y
黄禄渊,张贝,程惠红,罗纲,瞿武林,石耀霖. 2017. 2015年智利Illapel MW8.3地震同震效应及其对南美大陆地震危险性影响[J]. 地球物理学报,60(1):163–173 Huang L Y,Zhang B,Cheng H H,Luo G,Qu W L,Shi Y L. 2017. The co-seismic effect of the 2015 Chile Illapel MW8.3 earthquake and its effect on seismic risk of South America[J]. Chinese Journal of Geophysics,60(1):163–173 (in Chinese)
赖院根,刘启元,陈九辉,刘洁,李顺成,郭飙,黄志斌. 2006. 首都圈地区横波分裂与地壳应力场特征[J]. 地球物理学报,49(1):189–196 doi: 10.3321/j.issn:0001-5733.2006.01.025 Lai Y G,Liu Q Y,Chen J H,Liu J,Li S C,Guo B,Huang Z B. 2006. Shear wave splitting and the features of the crustal stress field in the Capital Circle[J]. Chinese Journal of Geophysics,49(1):189–196 (in Chinese)
刘莎,吴朋,杨建思,苏金蓉. 2014. 吉林省前郭地区地震各向异性的初步探讨[J]. 地球物理学报,57(7):2088–2098 Liu S,Wu P,Yang J S,Su J R. 2014. Preliminary study of seismic anisotropy in Qianguo area,Jilin Province[J]. Chinese Journal of Geophysics,57(7):2088–2098 (in Chinese)
缪淼,朱守彪. 2012. 俯冲带上特大地震静态库仑应力变化对后续余震触发效果的研究[J]. 地球物理学报,55(9):2982–2993 Miao M,Zhu S B. 2012. A study of the impact of static Coulomb stress changes of megathrust earthquakes along subduction zone on the following aftershocks[J]. Chinese Journal of Geophysics,55(9):2982–2993 (in Chinese)
庞亚瑾,程惠红,张怀,石耀霖. 2017. 巴颜喀拉块体东缘形变及九寨沟地震孕震环境数值分析[J]. 地球物理学报,60(10):4046–4055 doi: 10.6038/cjg20171030 Pang Y J,Cheng H H,Zhang H,Shi Y L. 2017. Numerical modeling of crustal deformation in the eastern margin of the Bayan Har block and analysis of seismogenic environment of the 2017 Jiuzhaigou earthquake[J]. Chinese Journal of Geophysics,60(10):4046–4055 (in Chinese)
任俊杰,徐锡伟,张世民,罗毅,梁欧博,赵俊香. 2017. 东昆仑断裂带东端的构造转换与2017年九寨沟MS7.0地震孕震机制[J]. 地球物理学报,60(10):4027–4045 doi: 10.6038/cjg20171029 Ren J J,Xu X W,Zhang S M,Luo Y,Liang O B,Zhao J X. 2017. Tectonic transformation at the eastern termination of the Eastern Kunlun fault zone and seismogenic mechanism of the 8 August 2017 Jiuzhaigou MS7.0 earthquake[J]. Chinese Jour-nal of Geophysics,60(10):4027–4045 (in Chinese)
单新建,屈春燕,龚文瑜,赵德政,张迎峰,张国宏,宋小刚,刘云华,张桂芳. 2017. 2017年8月8日四川九寨沟7.0级地震InSAR同震形变场及断层滑动分布反演[J]. 地球物理学报,60(12):4527–4536 doi: 10.6038/cjg20171201 Shan X J,Qu C Y,Gong W Y,Zhao D Z,Zhang Y F,Zhang G H,Song X G,Liu Y H,Zhang G F. 2017. Coseismic deformation field of the Jiuzhaigou MS7.0 earthquake from Sentinel-1A InSAR data and fault slip inversion[J]. Chinese Journal of Geophysics,60(12):4527–4536 (in Chinese)
闻学泽,杜方,张培震,龙锋. 2011. 巴颜喀拉块体北和东边界大地震序列的关联性与2008年汶川地震[J]. 地球物理学报,54(3):706–716 doi: 10.3969/j.issn.0001-5733.2011.03.010 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)
谢祖军,郑勇,姚华建,房立华,张勇,刘成利,王毛毛,单斌,张会平,任俊杰,季灵运,宋美琴. 2018. 2017年九寨沟MS7.0地震震源性质及发震构造初步分析[J]. 中国科学: 地球科学,48(1):79–92 Xie Z J,Zheng Y,Yao H J,Fang L H,Zhang Y,Liu C L,Wang M M,Shan B,Zhang H P,Ren J J,Ji L Y,Song M Q. 2018. Preliminary analysis on the source properties and seismogenic structure of the 2017 MS7.0 Jiuzhaigou earthquake[J]. Science China Earth Sciences,61(3):339–352 doi: 10.1007/s11430-017-9161-y
徐纪人,赵志新. 2006. 中国岩石圈应力场与构造运动区域特征[J]. 中国地质,33(4):782–792 doi: 10.3969/j.issn.1000-3657.2006.04.008 Xu J R,Zhao Z X. 2006. Regional characteristics of the lithospheric stress field and tectonic motions in China and its adjacent areas[J]. Geology in China,33(4):782–792 (in Chinese)
徐锡伟,闻学泽,陈桂华,于贵华. 2008. 巴颜喀拉地块东部龙日坝断裂带的发现及其大地构造意义[J]. 中国科学: D辑,38(5):529–542 Xu X W,Wen X Z,Chen G H,Yu G H. 2008. Discovery of the Longriba fault zone in eastern Bayan Har block,China and its tectonic implication[J]. Science in China:Series D,51(9):1209–1223 doi: 10.1007/s11430-008-0097-1
徐锡伟,陈桂华,王启欣,陈立春,任治坤,许冲,魏占玉,鲁人齐,谭锡斌,董绍鹏,石峰. 2017. 九寨沟地震发震断层属性及青藏高原东南缘现今应变状态讨论[J]. 地球物理学报,60(10):4018–4026 doi: 10.6038/cjg20171028 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)
许忠淮. 2001. 东亚地区现今构造应力图的编制[J]. 地震学报,23(5):492–501 doi: 10.3321/j.issn:0253-3782.2001.05.005 Xu Z H. 2001. A present-day tectonic stress map for eastern Asia region[J]. Acta Seismologica Sinica,23(5):492–501 (in Chinese)
杨宜海,范军,花茜,高见,王朝亮,周鲁,赵韬. 2017. 近震全波形反演2017年九寨沟M7.0地震序列震源机制解[J]. 地球物理学报,60(10):4098–4104 doi: 10.6038/cjg20171034 Yang Y H,Fan J,Hua Q,Gao J,Wang C L,Zhou L,Zhao T. 2017. Inversion for the focal mechanisms of the 2017 Jiuzhaigou M7.0 earthquake sequence using near-field full waveforms[J]. Chinese Journal of Geophysics,60(10):4098–4104 (in Chinese)
张军龙,任金卫,付俊东,胡朝忠,熊仁伟,陈长云,杨攀新. 2012. 东昆仑断裂带东部塔藏断裂地震地表破裂特征及其构造意义[J]. 地震,32(1):1–16 doi: 10.3969/j.issn.0253-3782.2012.01.001 Zhang J L,Ren J W,Fu J D,Hu C Z,Xiong R W,Chen C Y,Yang P X. 2012. Earthquake rupture features and tectonic signifi-cance of the Tazang fault in the eastern part of the east Kunlun fault zones[J]. Earthquake,32(1):1–16 (in Chinese)
张培震,邓起东,张国民,马瑾,甘卫军,闵伟,毛凤英,王琪. 2003. 中国大陆的强震活动与活动地块[J]. 中国科学: D辑,33(增刊1):12–20 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,46(S2):13–24
张旭,冯万鹏,许力生,李春来. 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)
赵博,高原,马延路. 2018. 利用面波振幅谱确定四川九寨沟M7.0地震震源深度[J]. 科学通报,63(13):1223–1234 Zhao B,Gao Y,Ma Y L. 2018. Using surface wave amplitude spectra to estimate the source depth of Sichuan Jiuzhaigou M7.0 earthquake[J]. Chinese Science Bulletin,63(13):1223–1234 (in Chinese) doi: 10.1360/N972017-01150
中国地震局震害防御司. 1999. 中国近代地震目录(公元1912—1990年MS≥4.7)[M]. 北京: 中国科学技术出版社: 1–637. Department of Earthquake Disaster Prevention, China Earthquake Administration. 1999. Recent Earthquake Catalog of China (1912—1990, MS≥4.7)[M]. Beijing: Chinese Sciences and Technology Press: 1–637 (in Chinese).
周荣军,李勇,Densmore A L,Ellis M A,何玉林,王凤林,黎小刚. 2006. 青藏高原东缘活动构造[J]. 矿物岩石,26(2):40–51 doi: 10.3969/j.issn.1001-6872.2006.02.007 Zhou R J,Li Y,Densmore A L,Ellis M A,He Y L,Wang F L,Li X G. 2006. Active tectonics of the eastern margin of the Tibet Plateau[J]. Journal of Mineralogy and Petrology,26(2):40–51 (in Chinese)
Burchfiel B C,Chen Z L,Liu Y,Royden L H. 1995. Tectonics of the Longmen Shan and adjacent regions,central China[J]. Int Geol Rev,37(8):661–735 doi: 10.1080/00206819509465424
Crampin S. 1999. Calculable fluid-rock interactions[J]. J Geol Soc,156(3):501–514 doi: 10.1144/gsjgs.156.3.0501
Crampin S,Peacock S. 2005. A review of shear wave splitting in the compliant crack-critical anisotropic Earth[J]. Wave Motion,41(1):59–77 doi: 10.1016/j.wavemoti.2004.05.006
Crampin S,Gao Y. 2006. A review of techniques for measuring shear-wave splitting above small earthquakes[J]. Phys Earth Planet Inter,159(1/2):1–14
Gan W J,Zhang P Z,Shen Z K,Niu Z J,Wang M,Wan Y G,Zhou D M,Cheng J. 2007. Present-day crustal motion within the Tibetan Plateau inferred from GPS measurements[J]. J Geophys Res,112(B8):B08416 doi: 10.1029/2005JB004120
Gao Y,Crampin S. 2004. Observations of stress relaxation before earthquakes[J]. Geophys J Int,157(2):578–582 doi: 10.1111/gji.2004.157.issue-2
Gao Y,Crampin S. 2008. Shear-wave splitting and earthquake forecasting[J]. Terra Nova,20(6):440–448 doi: 10.1111/ter.2008.20.issue-6
Gao Y,Wu J,Fukao Y,Shi Y T,Zhu A L. 2011. Shear wave splitting in the crust in North China:Stress,faults and tectonic impli-cations[J]. Geophys J Int,187(2):642–654 doi: 10.1111/gji.2011.187.issue-2
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[J]. J Geophys Res,109(B2):B02303
Liu S,Crampin S,Luckett R,Yang J S. 2014. Changes in shear wave splitting before the 2010 Eyjafjallajökull eruption in Iceland[J]. Geophys J Int,199(1):102–112 doi: 10.1093/gji/ggu202
Stein R S. 2005. Earthquake conversations[J]. Scientific American,15:82–89 doi: 10.1038/scientificamerican0705-82sp
USGS. 2017. M6.5: 36 km WSW of Yongle, China[EB/OL]. [2018-02-12]. https://earthquake.usgs.gov/earthquakes/eventpage/us2000a5x1/executive.
Yang Y H,Liang C T,Li Z Q,Su J R,Zhou L,He F J. 2017. Stress distribution near the seismic gap between Wenchuan and Lushan earthquakes[J]. Pure Appl Geophys,174(6):2257–2267 doi: 10.1007/s00024-016-1360-6
Zatsepin S V,Crampin S. 1997. Modelling the compliance of crustal rock: I. Response of shear-wave splitting to differential stress[J]. Geophys J Int,129(3):477–494