Relocation of the aftershock sequence of Milin MS6.9 earthquake in 2007 and spatio-temporal variation characteristics of b-value
-
摘要: 本文采用双差定位方法对2017年西藏米林MS6.9地震的余震序列进行重定位,获得了较准确的余震空间分布。在此基础上,计算了米林地震震源区的b值空间分布,并结合前人资料研究了米林地震前后的b值变化。重定位结果显示:米林地震余震序列呈NW−SE向分布,主要分布在雅鲁藏布江大拐弯缝合带和比鲁构造岩片的北东边界,震源深度集中于3—20 km范围内。b值时间分布显示:米林地震震前震源区的b值降低,表明震前震源区存在较强的应力积累;地震后震源区的b值呈现跳跃式上升,反映出震源区应力释放,同时随时间的推移,b值逐渐趋于稳定。米林地震震源区的b值在0.52—1.35范围内变化,15 km深度以上b值在东、西方向上存在差异,15 km以下在东北、南西方向上存在差异,由此推测不同的b值分布与地下的结构特征有关,浅层的b值变化与震源区的破裂程度相关,深部变化则反映了不同构造单元的岩性差异。Abstract: The aftershock sequence of the Milin MS6.9 earthquake in Tibet in 2017 are relocated by double-difference algorithm in this paper. Based on the relocated results, the spatial distribution of b-value in the source region of Milin earthquake is obtained. Combining historical b-value data from previous researches, the temporal distribution of b-value in the studied area is also obtained. The relocation results show that the aftershock sequence of the Milin earthquake are distributed in NE−SE direction along the NE boundary of Indus-Yarlung Tsangpo big-turned suture and Bilu tectonic slice, most aftershocks are concentrated in 3−20 km depth. The temporal distribution of b-value shows before Milin earthquake b-value is lower than that in the year of 2008, indicating that the stress accumulation in the studied area is strong before Milin earthquake. After Milin earthquake, there was a step-rise of b-value, showing the releases of stress, and b-value tended to stable. b-value ranges in 0.52−1.35 in the source region of Milin earthquake, the b-value variation above the depth of 15 km is distributed along the east-west direction, and below 15 km along the NE−SW direction. Therefore, it is deduced that different b-value distribution represents different structure characteristics, i.e., the b-value variation of the shallow area is related to the degree of rupture in the source region, and the b-value variation of deep area reflects the lithology difference of different tectonic units.
-
引言
南极洲位于地球最南端,是冈瓦纳古陆的中心区域(Craddock,1982),以横贯南极山脉为界分为东南极和西南极,南极半岛在西南极的北部。西南极是在冈瓦纳裂解中形成的活跃造山带,南极半岛处于古太平洋板块、南极板块和南美板块的交会部位,南极半岛的最北端(含南设得兰群岛)是南极大陆遗留的最后一段冈瓦纳俯冲陆缘(Kaminuma,2006;何高文等,2018),也是南极大陆现存的唯一一套“沟—弧—盆”体系(吴水根,吕文正,1988)。南设得兰群岛由欺骗岛、利文斯顿岛、格林尼治岛、罗伯特岛、纳尔逊岛、乔治王岛和布里奇曼岛等一系列岛屿组成,它和南极半岛原本为一个整体,是叠加在古太平洋板块边缘活动带上的岩浆弧(刘小汉等,1991;姚伯初等,1995)。自中生代起,大洋板块持续向南极半岛下方俯冲,尤其在大洋板块的扩张中心阿鲁克脊的俯冲期间,南极半岛的火山岛弧逐渐形成。至少在4 Ma之前,阿鲁克脊与南极大陆架边缘的相互作用导致了布兰斯菲尔德海峡的打开和南设得兰群岛与南极半岛的分离(Barker,1982;陈圣源等,1997;王改云等,2013;Li et al,2021)。
南设得兰海沟是现今菲尼克斯板块向南极半岛下俯冲的边界。有人认为俯冲作用在约4 Ma就结束了(Barker,1982),但也有人认为俯冲仍在继续,只是俯冲的速率减小了(Pelayo,Wiens,1989;Robertson Maurice et al,2003)。海洋热流值测量、重震联合反演和沉积物分析的研究结果均显示布兰斯菲尔德海峡中存在较为频繁的火山和岩浆活动(Lawver et al,1996;何高文等,2018;马龙,邢健,2020),且受到东西两侧英雄断裂带和沙克尔顿断裂带的控制。1969年和1970年在海峡西部的欺骗岛还发生过两次火山喷发,后续经历了多次的地震动荡期(Jiménez‑Morales et al,2022)。Kaminuma (1995,2006)根据国际地震中心(International Seismological Centre,缩写为ISC)公布的地震目录发现,尽管大陆整体活动性较弱,但火山活动较强的欺骗岛地区是一个地震相对活跃的区域。Robertson Maurice等(2003)利用短期海底地震仪观测发现在海峡中段也存在与海底火山和岩浆活动密切相关的地震活动。关于菲尼克斯板块向南设得兰群岛下的俯冲作用是否仍然进行着?在布兰斯菲尔德海峡,与裂谷作用导致的岩浆和火山活动相关的地震活动有何特征?该区域地震分布特征如何?这些科学问题仍需进一步研究。
南设得兰群岛及周边地区受俯冲、裂谷、岩浆和火山等活动影响,构造复杂,是开展地震学研究的理想场所。但南极大陆极端的气候条件限制了地震观测,相对于其它大陆,南极大陆的地震学研究程度较低。由于全球地震台网(Global Seismographic Network,缩写为GSN)的地震台站分布稀疏,距离南设得兰群岛最近的为南极半岛的PMSA台,距离也在500 km以上。全球地震目录中该区域的地震事件多在M4.3以上,说明该区域中小地震的监测能力不足,从而也影响了对该区域详细地震活动特征的分析。
中国南极长城站(Chinese Antarctic Great Wall Station,缩写为CCZ)地震台位于南设得兰群岛的乔治王岛上(图1),于1986年中国第二次南极科学考察期间建立并开始观测(林云芳,贾国强,1989)。长城站地震台早期采用模拟地震仪和图纸记录地震波形,受观测条件和仪器性能的限制,早期观测的数据制约了科学研究工作。中途由于南极科考任务的调整等原因,长城站地震观测一度停止,直至2007年第24次南极科学考察在原址架设了宽频带数字地震仪才重新恢复了观测任务。但是站区建设对台站观测造成了干扰,中国地震局地球物理研究所在2010年第27次南极科学考察期间为地震台重新选址,避开了干扰源,重建了长城站地震台,恢复了地震观测(常利军等,2012),观测质量显著提高,不仅记录到清晰的远震,还记录到发生在长城站附近的近震和疑似冰震波形(常利军,2013)。
图 1 研究区构造背景图中绿色三角形是中国南极长城站(CCZ)地震台的位置,地震数据来自美国国家地质调查局地震信息中心提供的1990—2017年的美国国家地震监测台网系统(Advanced National Seismic System,缩写为ANSS)综合地震目录(USGS,Earthquake Hazards Program,2017),震源机制解来自全球质心矩张量项目提供的1990—2017年的全球质心矩张量(Centroid-moment tensors,缩写为CMT)目录(Erik et al,2013),断层数据来自德克萨斯大学地球物理研究所的UTIG数据库(Institute of Geophysics,University of Texas,2015)和Dziak等(2010)。右上角插图为斯科舍板块及其周边的地震与台站分布,其中蓝色三角形是全球地震台网(GSN)的台站,蓝色方框代表主图的位置Figure 1. Tectonic setting of studied areaThe green triangle is the location of Chinese Great Wall Seismic Station (CCZ),Antarctic. The seismic data is from the ANSS Comprehensive Earthquake Catalog between 1990 and 2017,provided by USGS National Earthquake Information Center (USGS,Earthquake Hazards Program,2017). The focal mechanism solution is from CMT catalog between 1990 and 2017 provided by the Global Centroid-Moment-Tensor Project (Erik et al,2013). The fault data is from the University of Texas Institute for Geophysics database (Institute of Geophysics,University of Texas,2015) and Dziak et al (2010). The insert in the upper right is the distribution of earthquakes and seismic stations in and around the Scotia Plate,the blue triangles are Global Seismographic Network Stations,and the blue rectangle represents the location of the main figure本文对长城站地震台的三分量近震记录进行震相拾取、振幅测量和交错网格搜索定位,确定了长城站地区的地震活动特征,进一步讨论了与板块碰撞引起的俯冲活动和裂谷作用相关的孕震动力学机制等问题。
1. 资料和方法
1.1 数据及处理
1.1.1 台站和数据
重新选址和改建后的地震台位于长城站南部的八达岭山上(58.96°W,62.22°S),高程为41 m。台站布设了CMG-3ESP宽频带地震计,频带为50 Hz—60 s,采用REFTEK130数字采集器,采样率设为50 Hz,GPS实时授时(常利军等,2012)。本研究收集了CCZ台2015—2017年的三分量连续地震波形资料。但由于南极地区极端恶劣天气的影响导致仪器部分时段工作不正常,2016年1月16日至2月29日和2016年11月25日至2017年1月14日两个时段的波形记录异常,无法使用。
1.1.2 数据处理
数据处理包括震相拾取和振幅测量。在震相拾取过程中,为了准确地识别近震事件,根据近震记录的频率范围对三分量波形数据进行1—20 Hz的带通滤波;然后通过目视方法识别地震信号并拾取近震事件的P波和S波的到时;最后计算P波和S波的到时差$ {t}_{\mathrm{S}-\mathrm{P}} $。其中信噪比较差的地震事件波形(图2c)只用作地震数量统计,不进行定位。挑选具有清晰Pg波和Sg波的事件进行振幅测量(图2a),获得三分量Pg波初动信号前半个周期的最大振幅(图2b)和两水平分量的Sg波最大振幅,用于反方位角、入射角和震级的计算。垂直向记录uU向上为正,东西向记录uE向东为正,南北向记录uN向北为正。当uU<0时,地震射线背向震中运动,此时uE和uN需取相反数。入射角$ \overline{{i}_{\mathrm{P}}} $的计算公式如下:
图 2 中国南极长城站地震台的近震波形记录(a) 可定位地震事件的三分量波形记录;(b) 可定位地震事件的P波初动信号三分量记录;(c) 不可定位地震事件的三分量波形记录;(d) 不可定位地震事件的P波初动信号三分量记录Figure 2. Near earthquake waveforms recorded at CCZ Station(a) Three-component waveform recordings of a locatable seismic event;(b) Three-component recordings of P wave initial motion signal of a locatable seismic event;(c) Three-component waveform recordings of an unlocatable seismic event; (d) Three-component recordings of P wave initial motion signal of an unlocatable seismic event$$ \overline{{i}_{\mathrm{P}}}=\mathrm{atan}\left(\frac{\sqrt{{{u}_{{\mathrm{E}}}^{2}} + {{u}_{{\mathrm{N}}}^{2}}}}{\left|{u}_{{\mathrm{U}}}\right|}\right) . $$ (1) 根据《地震学导论》(万永革,2016)4.2.4节,由P波地震记录计算的入射角$ \overline{{i}_{\mathrm{P}}} $是视入射角,需用下式
$$ \mathrm{sin}{i}_{\mathrm{P}}=\frac{{v}_{\mathrm{P}}}{{v}_{\mathrm{S}}}\mathrm{sin}\frac{\overline{{i}_{\mathrm{P}}}}{2} $$ (2) 计算真入射角$ {i}_{\mathrm{P}} $,式中vP和vS分别是P波速度和S波速度。
1.2 方法
本研究使用交错网格搜索法对长城站地震台的连续波形数据进行单台定位,所谓交错网格,即先建立一套网格进行单台搜索定位,然后将网格进行空间平移,基于平移后的网格采用同样的数据进行定位,获取不同网格下的定位结果。单套网格如果划分太细会导致搜索获得的定位结果稳定性差,如果太粗则会造成定位精度低,而交错网格搜索法可以有效克服单套网格的缺点。通过多套网格在空间上的交叉,提高地震定位的空间采样率,保证定位结果的稳定性和精度。
1.2.1 构建速度模型
速度模型是地震定位最重要的参数之一,本研究使用的长城站地区速度模型(表1)参考了波兰在1979—1991年进行西南极深地震测深实验获得的南设得兰群岛及其周边地区的深部结构特征(Janik et al,2006),以及中国第28次南极科学考察期间,常利军(2013)在法尔兹半岛进行的近地表浅层结构探测的实验结果。
表 1 中国南极长城站地区速度模型Table 1. The velocity model of the Chinese Antarctic Great Wall Station area深度/km vP/(km·s−1) vS/(km·s−1) 0 4.10 2.20 2 5.70 3.10 9 6.65 3.59 20 7.28 3.94 30 8.05 4.43 50 8.30 4.52 1.2.2 交错网格搜索定位
交错网格搜索定位方法的具体操作包括如下几个步骤:
1) 单套网格的构建。根据实测数据走时信息选择研究区的网格划分范围和网格间距,以台站位置为中心建立三维网格,水平方向为经度-纬度网格(单位为度),深度方向以震中距-深度构建网格(单位为km)。根据研究区位置以及由到时差$ {t}_{\mathrm{S}-\mathrm{P}} $和虚波速度[vPvS/(vP-vS)=8 km/s]估算的可定位地震震中距范围设置网格尺寸,考虑到长城站位于约62°S的高纬度地区,这里1个经度对应的陆地距离大约是1个纬度的1/2,因此将经度和纬度范围分别设置为2°和1°,在其内部以0.02°×0.01°为单元划分小格。深度网格根据可定位地震的震中距范围将震中距和深度范围分别设为200 km和160 km,内部小格尺寸对应经纬度网格设为1 km×1 km。
2) 构建走时表。计算水平网格点坐标(lon’,lat’)与台站之间的距离${\varDelta }_{1}' $和反方位角$ {ba{\textit{z}'}}$,根据TauP函数计算垂向网格点($ h' $,$ {\varDelta }_{1}' $)的P波、S波理论走时,获得P波和S波的理论走时差$ {t}_{\mathrm{S}-\mathrm{P}}' $,并计算网格点的入射角$ {i}_{\mathrm{P}}' $。
3) 单套网格搜索定位。以震相拾取和振幅测量时获取的参数作为搜索标准与网格的走时表参数进行匹配,首先在垂向网格参数信息中搜索$ {t}_{\mathrm{S}-\mathrm{P}}' $与$ {t}_{\mathrm{S}-\mathrm{P}} $的差在±0.03 s内且$ {i}_{\mathrm{P}} $与$ {i}_{\mathrm{P}}' $最接近的网格点,确定地震的震源深度h和临时震中距$ {\varDelta }_{1}' $,将$ {ba{\textit{z}}}' $与测量的反方位角baz的差在±1°范围内的水平网格点作为垂向网格的大致方位,$ {\varDelta }_{1}' $与$ {\varDelta }_{2}' $最接近的水平网格点坐标就是震中的位置(lon,lat)。
4) 交错网格搜索定位。为了提高单台地震定位的精度和搜索稳定性,将水平网格以小格尺寸的1/4 (向东0.005°,向北0.002 5°)为步长沿经度和纬度进行平移,垂向网格以小格尺寸的1/4 (0.25 km)为步长,沿水平方向和垂向进行平移,共移动15次,形成如图3所示的交错网格,将网格的空间采样率提高了16倍。将16个水平网格和16个垂向网格两两组合为一套网格模型,进行上述单套网格搜索定位,可获得256套经度、纬度和深度结果,其均值即为最终的震源位置。最后根据速度模型正演地震射线获得理论P波走时,计算地震事件的发震时刻。
图 3 交错网格的构建(a) 经度-纬度网格的平移过程;(b) 震中距-深度网格的平移过程黑色网格表示起始网格位置,红色、蓝色和绿色网格分别表示每一次移动后的网格位置,三角形表示台站。此示意图按小格尺寸的1/2为步长进行平移Figure 3. The construction of interlaced gird(a) The translation process of the longitude-latitude grid;(b) The translation process of the epicentral distance-depth grid The black grid represents initial grid position,the red,blue and green grids represent the grid position after each move respectively. The triangle represents seismic station. This diagram is translated by 1/2 of the small lattice size per step5) 定位结果可靠性评定。标准差可以反映一组数据的离散程度,通过计算同一套数据的多套不同网格的定位结果的标准差可以有效地评价定位结果的稳定性,本文基于256套网格的定位结果,分别计算了每个地震事件的标准差,标准差越小说明网格模型的差异对结果的影响越小,单台定位结果越稳定。
图2a为CCZ台在国际时(Coodinated Universal Time,缩写为UTC)2017年3月31日10:05:48.688记录的一次近震事件波形,该事件的定位过程如下:首先在图2a中量取P波到时为10:05:48.688,S波到时为10:06:04.712,在图2b中测得P波的三分量振幅uE,uN和uU分别为469 counts,−523 counts和2 684 counts;获得该事件的到时差$ {t}_{\mathrm{S}-\mathrm{P}}$为16.0 s,反方位角$ ba{\textit{z}} $为318.116°,真入射角$ {i}_{\mathrm{P}}= $13.762°。将以上测量参数输入交错网格模型的走时表中搜索后,匹配获得最佳反方位角($ {ba{\textit{z}}}'= $318.048°)、入射角($ {i}_{\mathrm{P}}'= $14.758°)和震源位置(61.782°S,59.784°W,h=132.6 km),波形正演获得的P波走时为18.155 s,发震时刻为10:05:30.533,震中距标准差和深度标准差分别为0.5 km和0.3 km。
1.2.3 震级计算
震级的计算需要使用S波最大振幅、震中距和量规函数。由于不同地区的地震震源深度和地震波衰减特性有较大差异,针对不同地区的震级计算公式往往都不相同,并且使用的仪器特点对量规函数也有影响。因此一般要将宽频带地震记录进行仿真处理再使用相应的量规函数进行校正(刘瑞丰等,2014),但是长城站地区的量规函数目前还未进行测定。
CCZ台使用了频带宽度为50 Hz—60 s的CMG-3ESP地震计,近震S波的主要频段为1—20 Hz,位于该仪器频谱响应的平坦区间内,这说明S波信号的主要频段在该仪器的记录中是真实的,因此可以将灵敏度s应用于本文的震级计算,经计算,该仪器的灵敏度为0.000 79 μm/(s·count)。本研究使用《新地震观测手册》(彼得·鲍曼,2002)3.2.4节的地方震级公式和量规函数进行震级的计算,地方震级计算式如下:
$$ {M}_{\mathrm{L}}=\mathrm{lg}{A}_{\mathrm{m}\mathrm{a}\mathrm{x}}-\mathrm{lg}{A}_{0} \text{,} $$ (3) 式中$ {A}_{\mathrm{m}\mathrm{a}\mathrm{x}}= ( {{A}_{\mathrm{E}} + {A}_{\mathrm{N}}} ) s/{2} $,AE和AN分别为东西向和南北向S波记录的最大振幅,$ -\mathrm{lg}{A}_{0} $是量规函数。
按照该方法,我们挑选了2015—2017年ANSS目录中4次地震重新计算了震级,利用该目录提供的地震位置计算震中距,并从连续波形中量取该地震事件的S波振幅。通过二者震级的对比(表2)可见,相对于ANSS目录中的体波震级,本方法计算的震级(ML)略偏大,最大偏差为0.4级,与ANSS目录中的体波震级偏差不大。
表 2 ANSS目录震级和本方法计算结果的对比Table 2. The magnitude comparison of ANSS catalog and calculation result of this method发震时刻(UTC) 南纬/° 西经/° 深度/km Mb ML (本文震级) 年−月−日 时:分:秒 2 016−05−30 15:35:09.080 61.65 54.70 10.0 4.6 4.9 2 016−07−02 20:02:00.060 61.35 53.93 10.0 4.6 4.9 2 016−07−16 17:30:40.090 60.36 57.56 16.5 4.3 4.7 2 016−10−17 00:13:56.420 63.29 61.91 10.0 4.9 4.9 2. 结果
本研究对CCZ台在2015—2017年记录的三分量连续波形资料进行了人工拾取,分别识别出827,334和291个近震事件,共计1 452次地震,对其中信噪比较高的地震拾取了直达P波和S波的到时,量取其P波三分量振幅和S波水平分量最大振幅,并进行单台定位和结果可靠性评定,获得了112个近震定位结果(表3和图4)。87%的事件震中距和深度标准差均小于网格划分尺寸1.0 km,说明本方法中不同网格模型的差异对结果的影响较小,单台定位结果稳定可靠。
表 3 2015—2017年CCZ台单台定位结果Table 3. Single station location results of CCZ Station during 2015−2017发震时刻(UTC) 南纬/° 西经/° 深度/km ML 震中距标准差/km 深度标准差/km 年−月−日 时:分:秒 2 015−01−09 02:52:57.142 62.32 58.33 47.0 2.1 0.6 0.8 2 015−01−11 18:55:35.171 62.65 59.49 46.6 3.2 0.3 0.5 2 015−01−13 02:26:02.496 62.41 58.36 44.3 2.3 0.4 0.6 2 015−01−17 15:08:20.536 62.69 59.48 44.9 3.0 0.4 0.5 2 015−01−19 18:21:52.508 62.25 59.30 43.6 2.0 1.0 0.3 2 015−01−26 23:24:17.910 62.67 59.56 36.9 3.4 0.3 0.6 2 015−02−01 19:15:40.762 62.69 59.29 2.1 3.5 0.9 1.4 2 015−02−01 21:56:40.560 62.67 59.54 15.7 3.6 0.4 2.2 2 015−02−06 00:16:49.919 62.33 58.98 14.8 2.1 2.4 1.1 2 015−02−22 05:07:40.639 62.33 59.64 45.0 3.6 0.3 0.5 2 015−02−12 01:23:56.584 62.67 59.66 27.3 3.9 0.3 1.4 2 015−02−14 01:52:26.810 62.35 59.29 42.7 1.7 0.7 0.4 2 015−02−16 04:36:05.947 62.61 59.45 54.1 2.7 0.4 0.5 2 015−02−16 09:37:06.332 62.67 59.40 47.1 2.9 0.4 0.5 2 015−03−01 19:38:22.289 62.26 59.44 70.4 3.4 0.8 0.3 2 015−03−08 21:29:07.408 62.64 59.40 48.4 4.5 0.4 0.5 2 015−03−09 17:03:42.891 62.64 59.42 49.1 4.2 0.4 0.5 2 015−03−10 12:17:11.868 62.69 59.53 31.5 3.1 0.3 1.1 2 015−03−27 08:49:32.449 62.67 59.38 41.6 3.7 0.4 0.5 2 015−03−27 11:50:30.010 62.30 59.50 77.1 1.9 0.6 0.3 2 015−04−05 10:38:31.014 62.60 59.61 48.1 3.0 0.4 0.5 2 015−04−07 15:01:46.761 62.69 59.41 43.9 3.6 0.4 0.6 2 015−04−11 07:13:53.799 62.05 59.10 62.9 1.6 0.8 0.3 2 015−04−14 17:24:41.493 62.35 59.41 49.9 2.8 0.6 0.3 2 015−05−09 05:12:09.465 62.45 59.30 48.9 3.0 1.0 0.8 2 015−05−31 11:24:33.141 62.68 59.55 21.4 3.8 0.4 1.5 2 015−07−04 09:14:21.220 62.87 60.16 75.9 3.4 0.3 0.5 2 015−07−05 04:30:40.211 62.83 59.92 91.6 3.6 0.4 0.5 2 015−07−12 08:00:12.972 62.33 60.07 134.1 3.7 0.6 0.3 2 015−07−12 23:11:11.567 62.51 59.98 70.2 3.0 0.4 0.6 2 015−07−20 03:05:51.509 62.64 59.54 43.4 3.8 0.4 0.5 2 015−07−20 05:59:45.273 62.64 59.43 48.4 4.8 0.4 0.5 2 015−07−20 07:04:04.974 62.04 58.50 47.6 2.5 0.9 0.6 2 015−07−20 11:23:14.377 62.68 59.32 48.1 3.5 0.4 0.5 2 015−07−22 11:32:19.363 62.00 59.36 105.8 1.7 0.6 0.3 2 015−07−24 22:09:27.906 62.42 59.51 58.6 2.7 0.6 0.4 2 015−07−27 18:10:25.901 62.28 59.19 26.7 2.2 0.9 0.4 2 015−07−31 18:27:09.665 62.90 60.08 98.9 4.0 0.4 0.5 2 015−08−02 05:54:17.571 61.98 59.22 117.9 1.3 0.7 0.3 2 015−08−02 17:53:21.228 62.12 59.65 76.8 4.5 0.5 0.3 2 015−08−13 17 07:03.185 61.74 58.54 69.3 3.5 0.4 0.5 2 015−08−23 18:01:25.544 62.65 59.51 43.6 3.3 0.4 0.5 2 015−08−24 05:33:26.643 62.66 59.60 34.9 3.1 0.3 0.8 2 015−08−26 07:21:16.635 61.47 58.73 121.4 3.5 0.9 0.6 2 015−08−29 21:02:32.865 62.01 58.98 58.3 1.9 2.0 0.3 2 015−09−04 01:00:05.695 61.72 58.99 100.9 2.7 0.6 0.4 2 015−09−06 03:22:49.571 62.28 59.33 37.1 2.3 0.7 0.3 2 015−09−09 15:43:45.491 61.86 58.37 78.6 3.6 0.5 0.3 2 015−09−12 08:08:14.483 62.28 59.29 59.2 2.7 0.8 0.3 2 015−09−15 21:15:24.730 62.32 58.89 14.0 1.8 1.0 0.7 2 015−09−19 17:37:15.264 62.59 59.69 55.8 3.7 0.3 0.5 2 015−09−26 21:19:39.555 61.99 58.86 59.0 2.8 0.7 0.3 2 015−09−27 00:13:48.033 62.13 59.42 85.1 2.1 0.7 0.3 2 015−10−01 07:15:24.196 62.63 59.59 67.7 4.0 0.4 0.6 2 015−10−06 07:04:50.301 61.75 58.16 83.8 3.1 0.6 0.8 2 015−10−11 00:24:17.187 62.23 59.15 56.2 2.2 1.8 0.3 2 015−10−12 00:09:00.274 61.94 58.25 62.3 2.7 0.7 0.9 2 015−10−30 09:02:58.971 61.93 59.40 117.7 2.6 0.6 0.3 2 015−11−05 10:50:09.692 62.56 59.61 76.2 3.9 0.8 0.5 2 015−12−17 20:10:02.746 62.67 59.86 69.4 3.3 0.3 0.5 2 016−03−05 17:48:06.327 62.28 59.53 55.9 2.5 0.6 0.3 2 016−04−10 21:07:44.857 62.57 59.55 65.9 3.8 0.4 0.5 2 016−05−16 15:46:42.359 62.28 59.41 59.9 3.4 0.7 0.3 2 016−06−07 09:39:09.115 61.77 58.37 51.6 3.6 0.3 0.5 2 016−06−08 16:44:36.534 62.09 59.82 80.8 2.8 0.6 0.4 2 016−06−19 00:34:16.557 62.47 59.57 62.1 2.8 0.4 0.3 2 016−06−21 00:53:16.317 62.65 59.81 71.7 4.3 0.5 0.6 2 016−07−08 06:59:53.855 62.56 59.53 61.6 2.7 0.4 0.4 2 016−07−15 12:52:28.546 62.56 58.33 40.1 4.4 0.3 0.8 2 016−08−11 22:37:06.427 62.27 59.30 57.6 2.2 0.9 0.3 2 016−08−19 13:08:28.771 62.41 59.17 19.6 1.3 0.6 0.6 2 016−08−20 04:48:49.164 62.31 59.61 61.1 2.1 0.6 0.4 2 016−08−22 10:08:51.241 62.00 58.29 61.3 2.4 0.4 0.3 2 016−08−26 16:11:27.861 62.27 59.76 50.9 1.9 0.5 0.5 2 016−09−07 11:30:14.309 61.94 59.75 121.3 2.8 0.6 0.3 2 016−09−16 18:31:24.635 62.51 59.60 43.6 2.9 0.4 0.5 2 016−11−01 23:20:23.770 62.34 59.85 88.4 3.7 0.6 0.4 2 016−11−03 09:46:21.038 62.67 59.84 84.5 4.0 0.3 0.4 2 016−11−11 04:42:38.088 62.36 59.22 33.5 1.8 0.8 0.4 2 016−11−11 20:06:48.364 62.40 59.61 58.8 2.8 0.5 0.3 2 017−01−16 00:38:56.890 62.60 59.48 54.1 4.0 0.4 0.5 2 017−01−16 01:18:34.637 62.62 59.63 44.6 4.0 0.4 0.6 2 017−01−16 01:19:41.900 62.58 59.55 53.9 3.7 0.4 0.5 2 017−01−24 16:03:04.799 62.13 59.37 47.6 1.8 0.8 0.3 2 017−01−25 17:45:15.526 62.25 59.25 32.4 2.2 1.0 0.3 2 017−02−15 04:07:15.207 62.13 59.35 62.1 3.1 0.7 0.3 2 017−02−24 00:32:35.366 61.93 58.65 64.9 3.0 0.6 0.4 2 017−03−01 04:29:05.009 62.54 59.80 47.9 3.5 0.4 0.5 2 017−03−02 08:12:21.446 61.97 59.93 102.3 3.0 0.5 0.3 2 017−03−06 17:54:49.262 62.02 58.41 35.4 2.4 0.4 0.5 2 017−03−07 10:27:33.329 62.33 59.89 74.6 2.7 0.5 0.3 2 017−03−09 16:07:42.708 62.64 59.73 70.6 2.6 0.5 0.6 2 017−03−15 17:59:53.782 62.41 59.67 42.9 2.5 0.4 0.5 2 017−03−31 10:05:30.533 61.78 59.78 132.6 4.1 0.5 0.3 2 017−04−12 08:25:59.180 61.66 58.74 59.7 4.5 0.4 0.5 2 017−05−24 09:54:32.085 62.18 58.81 26.4 1.6 1.5 0.3 2 017−05−30 13:34:00.388 62.27 59.75 76.7 3.1 0.5 0.3 2 017−05−31 20:40:08.639 61.98 58.65 84.6 3.6 0.6 0.3 2 017−07−25 05:34:31.746 62.67 59.82 63.0 3.9 0.3 0.5 2 017−08−09 11:22:15.439 63.00 60.83 77.4 3.8 0.4 0.5 2 017−08−23 03:36:11.790 62.38 59.74 81.8 3.3 0.5 0.3 2 017−08−24 10:41:36.910 62.69 59.72 58.1 3.2 0.3 0.5 2 017−08−29 22:47:54.601 62.31 59.18 35.3 1.7 1.3 0.3 2 017−09−15 16:14:20.470 62.65 59.99 62.9 3.8 0.4 0.5 2 017−10−02 17:14:29.735 62.26 59.16 38.2 1.8 1.0 0.3 2 017−10−03 05:40:27.856 62.17 59.15 38.5 1.9 1.1 0.3 2 017−10−28 00:51:54.113 62.37 59.56 59.3 2.4 0.7 0.4 2 017−11−12 04:58:43.418 62.10 59.14 65.7 2.9 0.9 0.3 2 017−11−13 02:23:55.693 61.96 60.30 103.8 2.5 0.6 0.4 2 017−11−15 13:24:29.718 62.26 59.33 54.3 3.2 0.8 0.3 2 017−11−24 07:37:47.287 62.93 60.30 102.8 4.2 0.4 0.6 2 017−12−28 23:20:23.074 62.26 59.32 53.9 3.3 0.9 0.3 图 4 长城站地区2015—2017年地震分布情况左上角插图给出了2015—2017年ANSS目录中在CCZ台500 km范围内的地震分布,蓝色方框为主图区域Figure 4. Earthquake distribution in the Great Wall Station area from 2015 to 2017The insert at the upper left corner shows the distribution of earthquakes within 500 km from the Great Wall Seismic Station in the ANSS catalog from 2015 to 2017. The blue rectangle is the area of the main figure图4左上角插图给出了ANSS目录同期在CCZ台500 km范围内的地震分布,图4主图显示了本文定位的地震分布。根据到时差与虚波速度(8 km/s)估计,2015—2017年CCZ台记录的1 425次地震的震中距均在500 km范围内,对其中112个信噪比高的近震事件进行定位,得到的最大震级为ML4.8,最小震级为ML1.3,震源深度最深达到134 km。ANSS目录在同期同范围只有18次地震,最大震级为MW5.3,最小震级为Mb4.2,最深的一次地震震源深度为17 km。虽然ANSS目录中的所有事件都能被CCZ台记录到,但是这两个目录中没有重复的定位地震事件,例如发震时刻为2016−05−20 15:35:09.08,震级为Mb4.6的地震,虽然其震级较大,在CCZ台的波形记录信号较强(图5),但是其P波初动信号不清晰,不满足本研究的定位要求。综上所述,CCZ台有效提高了南设得兰群岛地区的地震监测能力,填补了该区域地震目录中的小震,且利用直达P波初动振幅计算入射角也能更好地约束近震的震源深度。
在500 km范围内共监测到1 452次地震,其中定位的112次近震的震中距位于130 km以内。图6和图7分别显示了每天发生地震的数量和可定位地震震级随时间变化的情况。这些结果显示,2015年平均每天发生2—3次地震,在2015年1月11日和2015年7月20日达到14次;曾发生到时差9 s左右的地震丛集。2016—2017年期间,除去没有统计的两个时段,平均每天发生的地震不足两次,地震活动整体比较平稳,仅在2016年11月20日和2017年1月16日分别达到7次和8次,地震目录显示,2016年11月20日发生了到时差约为15 s的小震群,2017年1月16日发生了到时差8 s左右的小震群。统计表3发现,2015—2017年定位的地震震级主要集中在ML1.8—4.5,其中ML1.0—1.9有15次,ML2.0—2.9有36次,ML3.0—3.9有47次,ML4.0—4.9有14次(其中ML4.5以上有4次,包括3次ML4.5地震和1次ML4.8地震)。整体上,2015—2017年长城站地区的地震活动主要特征为在2015年年初和年中最强,2016年下半年略强于2017年,整体强于平静的南极大陆内部。
由图4可以看到,地震呈NE向展布于南设得兰群岛及其附近地区。在南设得兰群岛,地震主要集中在格林尼治岛、罗伯特岛和纳尔逊岛之间的小型海峡以及格林尼治岛南部区域,震源深度介于20—90 km,布兰斯菲尔德海峡中的欺骗岛附近存在震源深度超过60 km的地震,靠近海沟的弧前区域发生过11次震源深度超过100 km的地震,最深的已经达到了134 km。平行于群岛的AA′测线地震震源深度剖面(图8a)显示,在群岛区域内,CCZ台东西20 km范围内的地震震源深度在10—70 km范围内,而更远的区域震源深度基本上都在40 km以上,垂直于群岛的BB′测线地震震源深度剖面(图8b)显示,该地区的震源深度整体偏深,群岛及弧后区域震源深度小于100 km,弧前区域的震源深度较深。
图 8 测线AA′ (a)和BB′ (b)的震源深度分布剖面图图a,b上部为AA′和BB’测线纵截面测线位置的地形高程变化,剖面上投影了AA′测线南北20 km范围内地震和BB′测线所有可定位地震,横轴表示地震到测线的垂线与CCZ台的距离,以台站为原点,向西(北)为负,向东(南)为正Figure 8. Focal depths distribution along profile of survey lines AA′ (a) and BB′ (b)The upper part of the figure is the topographic elevation change along the longitudinal section of the survey lines. The earthquakes within 20 km from the north to the south of AA′ survey line and all locatable earthquakes of BB′ survey line are projected on the profile. The horizontal axis represents the distance between the vertical line of the earthquake to the survey line and the CCZ Station,with the station as the origin,negative to the west (north) and positive to the east (south)3. 讨论
中国南极长城站所在的南设得兰群岛地区位于南极半岛北部,受到南极板块、南美板块、斯科舍板块和菲尼克斯板块相互作用的影响,是冈瓦纳运动在南极地区遗留下来的最后一段俯冲带,具有独特且复杂的构造,存在裂谷、火山和俯冲等作用,是一个活跃的地震带。在该地区由于菲尼克斯板块向南设得兰群岛下方俯冲,俯冲板块上的阿鲁克扩张脊与俯冲边界南设得兰海沟不断碰撞,导致布兰斯菲尔德海峡伴随着火山和岩浆运动逐渐张开,形成了走向相互平行的南设得兰海沟、南设得兰群岛和布兰斯菲尔德海峡这套“沟—弧—盆”构造体系。
本研究在2015—2017年在靠近南设得兰海沟的弧前区域发现了11个震源深度超过100 km的地震事件,最深达到134 km,该地区的岩石圈厚度约为60—70 km (An et al,2015),这些地震发生在岩石圈之下。菲尼克斯板块向南设得兰群岛下的俯冲作用是该区域最主要的构造活动,关于俯冲活动是否仍在继续,目前有两种看法,Barker (1982)通过磁异常研究认为这里的俯冲活动在4 Ma前就已经停止,Robertson Maurice等(2003)和Almendros等(2020)则通过在南设得兰海沟附近发现的较深地震以及逆冲性质的震源机制解确定菲尼克斯板块正在以非常缓慢的速度(<1 cm/a)进行俯冲。Baker和Austin (1998)分析多道地震数据发现,布兰斯菲尔德海峡的伸展可能是由于南设得兰海沟的板块边界从西南方向打开了一个板块窗口从而引发了俯冲的回退。Dziak等(2010)通过研究布兰斯菲尔德海峡区域的地震活动证明了菲尼克斯板块的洋中脊已经消亡,向南设得兰群岛下方的主动俯冲已经停止,但是俯冲板块正在海沟和上地幔中垂直下沉,逆冲机制的地震事件仍可能发生在俯冲板片内部和边界带。1990—2017年全球质心矩张量(CMT)目录显示在靠近海沟的弧前区域存在逆冲走滑性质的事件(图1),且沙克尔顿断裂带所在位置存在一系列左旋走滑性质的震源机制,说明菲尼克斯板块在向南设得兰群岛地区运动。大多数学者认为菲尼克斯板块发生了俯冲回退,本文弧前地震事件的空间分布与上述证据相符,说明靠近海沟的中源深度地震可能是由俯冲活动引起的,这为菲尼克斯板块向南设得兰群岛下的俯冲活动仍在继续提供了地震学证据。但是由于目前此类地震事件较少,暂不能证明菲尼克斯板块的主动俯冲是否停止。
菲尼克斯板块的俯冲活动导致了布兰斯菲尔德海峡的扩张,在此过程中伴随着广泛的火山活动,欺骗岛就是布兰斯菲尔德海峡西南部的一座活火山,该火山在1969年和1970年还发生过两次喷发,后又在1992年、1999年和2015年经历了至少三次的动荡期。Jiménez-Morales等(2022)对欺骗岛2008—2015年期间的地震活动进行了研究,验证了该火山在此期间经历了岩浆在深部侵入并上升到上地壳并引起重复地震活动的过程,Almendros等(2018)对这次动荡期的地震定位结果显示地震事件主要分布于布兰斯菲尔德海峡中,岩浆向东的迁移导致了地震群分布于欺骗岛和利文斯顿岛附近,因此我们认为欺骗岛及其东侧和格林尼治岛南部的地震丛集可能与这次火山活动有关。2020年布兰斯菲尔德海峡东部的Orca海底火山附近也发生了类似的由岩浆活动引起的地震丛集(Olivet et al,2021;Poli et al,2022),由此我们认为发生在布兰斯菲尔德海峡东部的地震可能与火山活动有关。
南设得兰群岛地区处于海沟和海峡中间,受到了俯冲和岩浆活动的影响,以及斯科舍板块的挤压作用(González-Casado et al,2000),该区域的地震较为集中地分布在格林尼治岛、罗伯特岛和纳尔逊岛南北两侧的海岸线附近以及岛屿之间的小型海峡中,乔治王岛北侧也有少许分布。纳尔逊岛附近30 km深度之内的地震可能由与板块挤压有关的壳内应力引起。岛屿之间两个NW-SE向分布的地震条带与地形的吻合表明该区域可能存在一些平行于沙克尔顿断裂带和英雄断裂带的小断裂,它们将南设得兰群岛划分为具有不同岩浆活动和沉积历史的断块(吕文正,吴水根,1989;王公念等,1997),这些近乎平行的断裂可能与深海海底发育的大型断裂带有关(González-Ferrán,1983)。2014—2015年欺骗岛和利文斯顿岛先后经历了由岩浆侵入引发的地震群,该区域的地壳厚度约为30 km (Janik et al,2006;Parera-Portell et al,2021),而格林尼治岛南侧平行于群岛走向的地震条带震源深度普遍大于40 km,最深已经到达90 km,在莫霍面之下,且在空间展布上比Almendros等(2018)的定位结果更偏东北一点,整体比较一致,我们认为这部分地震可能是利文斯顿岛震群的余震,与其一样是由岩浆活动引起的。Robertson Maurice等(2003)认为1997—1999年间南设得兰地区的地震活动显示出了年轻海洋岩石圈的缓慢俯冲活动,罗伯特岛和纳尔逊岛北侧的NW-SE向地震条带可能与其相关。
4. 结论
本研究分析了长城站地震台2015—2017年记录的三分量地震波形,共识别出1 425个近震事件,并对其中112个高信噪比事件进行了定位。通过分析长城站地区的地震活动特征可知,长城站地震台有效提高了南设得兰群岛及其周边地区的地震监测能力,震级下限达到了ML1.3,并且较好地约束了地震事件的震源深度。本研究不仅获得了与岩浆活动相关的地震事件,还在靠近南设得兰海沟的弧前区域获得了11次中源深度地震,直接给出了菲尼克斯板块向南设得兰群岛下的俯冲活动仍然进行着的地震学证据。
项目执行中得到了国家海洋局极地考察办公室、中国极地研究中心等单位和中国南极科学考察队的大力支持,作者在此一并表示感谢。
-
图 1 米林地震震中、台站分布及研究区构造背景图(修改自Xu et al,2012)
Figure 1. Epicentral location of the Milin earthquake,station distribution and tectonic settings of the studied area (revised after Xu et al,2012)
图 2 米林地震序列重定位结果
(a) 震中分布图,图中震源机制解数据来自GCMT (Dziewonski et al,1981;Ekström et al,2012);(b) 走时残差直方图
Figure 2. Relocation result of Milin earthquake sequence
(a) Distribution of Milin earthquake sequence,where the focal mechanism refers to GCMT (Dziewonski et al,1981;Ekström et al,2012) catalog;(b) Histogram of the travel time residuals after relocation
-
白玲,李国辉,宋博文. 2017. 2017年西藏米林6.9级地震震源参数及其构造意义[J]. 地球物理学报,60(12):4956–4963. doi: 10.6038/cjg20171234 Bai L,Li G H,Song B W. 2017. The source parameters of the M6.9 Mainling,Tibet earthquake and its tectonic implications[J]. Chinese Journal of Geophysics,60(12):4956–4963 (in Chinese).
程成,白玲,丁林,李国辉,杨建亚,许强. 2017. 利用接收函数方法研究喜马拉雅东构造结地区地壳结构[J]. 地球物理学报,60(8):2969–2979. Cheng C,Bai L,Ding L,Li G H,Yang J Y,Xu Q. 2017. Crustal structure of Eastern Himalayan Syntaxis revealed by receiver function method[J]. Chinese Journal of Geophysics,60(8):2969–2979 (in Chinese).
方亚如,蔡戴恩,刘晓红,李纪汉,郝晋昇,耿乃光. 1986. 含水岩石破裂前的声发射b值变化[J]. 地震,(2):1–6. Fang Y R,Cai D E,Liu X H,Li J H,Hao J S,Geng N G. 1986. The change of b-value of acoustic emission before ruptures of water-bearing rocks[J]. Earthquake,(2):1–6 (in Chinese).
黄亦磊,周仕勇,庄建仓. 2016. 基于地震目录估计完备震级方法的数值实验[J]. 地球物理学报,59(4):1350–1358. Huang Y L,Zhou S Y,Zhuang J C. 2016. Numerical tests on catalog-based methods to estimate magnitude of completeness[J]. Chinese Journal of Geophysics,59(4):1350–1358 (in Chinese).
李莹甄,殷娜,李小晗. 2014. 不同震级标度转换关系研究概述[J]. 地震工程学报,36(1):80–87. doi: 10.3969/j.issn.1000-0844.2014.01.0080 Li Y Z,Yin N,Li X H. 2014. Review of the conversional relationship for different magnitude scales[J]. China Earthquake Engi-neering Journal,36(1):80–87 (in Chinese).
邵翠茹. 2009. 雅鲁藏布大峡谷地区地震活动性研究[D]. 北京: 中国地震局地球物理研究所: 51–68. Shao C R. 2009. Seismicity of the Yarlung Tsangpo Grand Canyon Region, China[D]. Beijing: Institute of Geophysics, China Earthquake Administration: 51–68 (in Chinese).
王鹏,侯金欣,吴朋. 2017. 2017年九寨沟MS7.0地震序列活动特征分析[J]. 中国地震,33(4):453–462. doi: 10.3969/j.issn.1001-4683.2017.04.001 Wang P,Hou J X,Wu P. 2017. Temporal evolution of the seismicity of the 2017 Jiuzhaigou MS7.0 earthquake sequence[J]. Earthquake Research in China,33(4):453–462 (in Chinese).
王未来,吴建平,房立华,王长在. 2012. 2010年玉树MS7.1地震及其余震的双差定位研究[J]. 中国科学:地球科学,42(7):1037–1046. Wang W L,Wu J P,Fang L H,Wang C Z. 2013. Relocation of the Yushu MS7.1 earthquake and its aftershocks in 2010 from HypoDD[J]. Science China Earth Sciences,56(2):182–191. doi: 10.1007/s11430-012-4450-z
王未来,吴建平,房立华,来贵娟. 2014. 2014年云南鲁甸MS6.5地震序列的双差定位[J]. 地球物理学报,57(9):3042–3051. Wang W L,Wu J P,Fang L H,Lai G J. 2014. Double difference location of the Ludian MS6.5 earthquake sequences in Yunnan Province in 2014[J]. Chinese Journal of Geophysics,57(9):3042–3051 (in Chinese).
张盛峰,吴忠良,房立华. 2014. 双差(DD)定位地震目录能用于地震序列的统计地震学参数计算吗? :云南鲁甸MS6.5地震序列b值的空间分布[J]. 地震地质,36(4):1244–1259. doi: 10.3969/j.issn.0253-4967.2014.04.024 Zhang S F,Wu Z L,Fang L H. 2014. Can the DD-relocated earthquake catalogue be used for the statistical parameters of an earthquake sequence? :A case study of the spatial distribution of b-values for the aftershocks of the 2014 Ludian MS6.5 earthquake[J]. Seismology and Geology,36(4):1244–1259 (in Chinese).
张勇, 许力生, 陈运泰. 2017. 2017年11月18日西藏米林6.9级地震(应急处置科技产品报告)[EB/OL]. [2018-04-20]. http://www.cea-igp.ac.cn/tpxw/276003.html. Zhang Y, Xu L, Chen Y T. 2017. Earthquake summary poster for November 18, 2017, Milin, Tibet, MS6.9 earthquake[EB/OL]. [2018-04-20]. http://www.cea-igp.ac.cn/tpxw/276003.html (in Chinese).
郑钰,杨建思. 2008. 双差算法的剖析及参数对定位的影响[J]. 地震地磁观测与研究,29(3):85–93. doi: 10.3969/j.issn.1003-3246.2008.03.016 Zheng Y,Yang J S. 2008. Analysis of double-difference algorithm and the affect of its parameter in location[J]. Seismological and Geomagnetic Observation and Research,29(3):85–93 (in Chinese).
Aki K. 1965. Maximum likelihood estimate of b in the formula logN=a-bM and its confidence limits[J]. Bull Seismol Soc Am,43:237–239.
Amitrano D. 2003. Brittle-ductile transition and associated seismicity:Experimental and numerical studies and relationship with the b value[J]. J Geophys Res,108(B1):ESE 19-1–ESE 19-15.
Dziewonski A M,Chou T A,Woodhouse J H. 1981. Determination of earthquake source parameters from waveform data for studies of global and regional seismicity[J]. J Geophys Res,86:2825–2852. doi: 10.1029/JB086iB04p02825
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:1–9.
Frohlich C,Davis S D. 1993. Teleseismic b values;Or,much ado about 1.0[J]. J Geophys Res,98(B1):631–644. doi: 10.1029/92JB01891
Gupta H K. 2002. A review of recent studies of triggered earthquakes by artificial water reservoirs with special emphasis on earthquakes in Koyna,India[J]. Earth-Sci Rev,58(3/4):279–310.
Gutenberg B, Richter C F. 1954. Seismicity of the Earth and Associated Phenomena[M]. 2nd ed. Princeton: Princeton University Press: 1–310.
Mogi K. 1962. Study of elastic shocks caused by the fracture of heterogeneous materials and its relations to earthquake phenomena[J]. Bull Earthq Res Inst,Univ Tokyo,40(1):125–173.
Mogi K. 1967. Regional variations in magnitude-frequency relation of earthquakes[J]. Bull Earthq Res Inst,Univ Tokyo,45(2):313–325.
Nakaya S. 2006. Spatio-temporal variation in b value within the subducting slab prior to the 2003 Tokachi-Oki earthquake (M8.0),Japan[J]. J Geophys Res,111(B3):B03311.
Nuannin P,Kulhánek O,Persson L. 2012. Variations of b-values preceding large earthquakes in the Andaman-Sumatra subduction zone[J]. J Asian Earth Sci,61:237–242. doi: 10.1016/j.jseaes.2012.10.013
Papazachos B C. 1974. Dependence of the seismic parameter b on the magnitude range[J]. Pure Appl Geophys,112(6):1059–1065. doi: 10.1007/BF00881508
Rydelek P A,Sacks I S. 1989. Testing the completeness of earthquake catalogues and the hypothesis of self-similarity[J]. Nature,337(6204):251–253. doi: 10.1038/337251a0
Scholz C H. 1968. The frequency-magnitude relation of microfracturing in rock and its relation to earthquakes[J]. Bull Seismol Soc Am,58(9):399–415.
Schorlemmer D,Wiemer S,Wyss M,Jackson D D. 2004. Earthquake statistics at Parkfield:2. Probabilistic forecasting and testing[J]. J Geophys Res,109(B12):B12308. doi: 10.1029/2004JB003235
Schorlemmer D,Wiemer S,Wyss M. 2005. Variations in earthquake-size distribution across different stress regimes[J]. Nature,434(7058):539–542.
Shi Y L,Bolt B A. 1982. The standard error of the magnitude-frequency b value[J]. Bull Seismol Soc Am,72(5):1677–1687.
Smith W D. 1981. The b-value as an earthquake precursor[J]. Nature,289(5794):136–139. doi: 10.1038/289136a0
Waldhauser F,Ellsworth W L. 2000. A double-difference earthquake location algorithm:Method and application to the northern Hayward fault,California[J]. Bull Seismol Soc Am,90(6):1353–1368. doi: 10.1785/0120000006
Waldhauser F,Ellsworth W L. 2002. Fault structure and mechanics of the Hayward fault,California,from double-difference earthquake locations[J]. J Geophys Res,107(B3):ESE 3-1–ESE 3-15.
Wiemer S,Katsumata K. 1999. Spatial variability of seismicity parameters in aftershock zones[J]. J Geophys Res,104(B6):13135–13151. doi: 10.1029/1999JB900032
Wiemer S,Wyss M. 2000. Minimum magnitude of completeness in earthquake catalogs:Examples from Alaska,the western United States,and Japan[J]. Bull Seismol Soc Am,90(4):859–869. doi: 10.1785/0119990114
Wiemer S,Gerstenberge M,Hauksson E. 2002. Properties of the aftershock sequence of the 1999 MW7.1 Hector Mine earthquake:Implications for aftershock hazard[J]. Bull Seismol Soc Am,90(4):1227–1240.
Xu Z Q,Ji S C,Cai Z H,Zeng L S,Geng Q R,Cao H. 2012. Kinematics and dynamics of the Namche Barwa Syntaxis,eastern Himalaya:Constraints from deformation,fabrics and geochronology[J]. Gondwana Res,21(1):19–36. doi: 10.1016/j.gr.2011.06.010
Zeitler P K, Meltzer A S, Brown L, Kidd W S F, Lim C, Enkelmann E. 2014. Tectonics and topographic evolution of Namche Barwa and the easternmost Lhasa block, Tibet[G]//Toward an Improved Understanding of Uplift Mechanisms and the Elevation History of the Tibetan Plateau. Colorado: Geological Society of America, 507: 23–58.