2020年1月19日新疆伽师 MS6.4地震序列的活动特征和发震构造
崔仁胜1,2), 赵翠萍2),*, 周连庆2), 陈阳2)
1)中国地震局地球物理研究所, 北京 100081
2)中国地震局地震预测研究所, 北京 100036
*通讯作者: 赵翠萍, 女, 1967年生, 研究员, 主要从事地震学和诱发地震研究, E-mail: zhaocp@cea-ies.ac.cn

〔作者简介〕 崔仁胜, 男, 1985年生, 现为中国地震局地球物理研究所地球物理学专业在读博士研究生, 副研究员, 主要从事地震学研究, E-mail: crs@ief.ac.cn

摘要

文中利用新疆地震台网固定台站和流动台站的观测资料研究了2020年1月19日新疆伽师 MS6.4地震序列的活动特征和发震构造。 重定位结果显示, 地震序列呈NNW和近EW 2个优势方向分布, 其中NNW向余震条带长约20km, 近EW向余震条带长约40km。 NNW向条带的发震断层倾角较陡, 倾向W; 近EW向条带的发震断层西侧倾角较陡, 向E逐渐变缓, 略S倾。 MS6.4主震深度为14.9km, 地震序列优势分布于10~20km深度, 近EW向的余震深度呈现由西向东逐渐变浅的特征。 震源机制解结果表明, MS6.4主震和2次5级以上余震的震源机制以逆冲为主, 2020年1月18日 MS5.7前震以走滑为主。 利用ISOLA方法对 MS6.4主震进行多源模型反演, 结果显示主震由走滑和逆冲性质的2个子事件组成。 推测 MS6.4地震发生在柯坪逆冲推覆构造带滑脱面以下的中下地壳, 主震由NNW向W倾的高角度隐伏断裂和近EW向S倾构造先后破裂组成。

关键词: 伽师 MS6.4地震; 重定位; 震源机制解; 发震构造
中图分类号:P315.2 文献标志码:A 文章编号:0253-4967(2021)02-0329-16
SEISMICITY FEATURE AND SEISMOGENIC FAULT OF THE MS6.4 EARTHQUAKE SEQUENCE ON JANUARY 19, 2020 IN JIASHI, XINJIANG
CUI Ren-sheng1,2), ZHAO Cui-ping2), ZHOU Lian-qing2), CHEN Yang2)
1)Institute of Geophysics, China Earthquake Administration, Beijing 100081, China
2)Institute of Earthquake Forecasting, China Earthquake Administration, Beijing 100036, China
Abstract

The 2020 MS6.4 Jiashi earthquake occurred on January 19, preceded by an MS5.7 foreshock on January 18. These two earthquakes occurred close in space and time raising the question of the relationship between the two events. Using the observation data recorded by fixed stations and temporal stations of Xinjiang seismic network, the seismicity feature and the seismogenic fault of the Jiashi MS6.4 earthquake sequence are studied in this paper. We relocated the Jiashi earthquake sequence from January 18 to August 31, 2020, and obtained the relocations of 1 460 earthquakes by the double-difference algorithm. The high-precision earthquake catalog reveals detailed spatial and temporal evolution of the earthquake sequence. The relocations show that the MS6.4 earthquake is located at 39.835°N, 77.148°E, and the focal depth is 14.9km. The earthquake sequence is distributed in two dominant directions, one is NNW direction, the other is near EW direction. The length of the NNW earthquake belt is about 20km, and the length of the near EW earthquake belt is about 40km. The dip angle of the seismogenic fault of the NNW earthquake belt is steep, dipping to the west. The dip angle of the seismogenic fault of the near EW earthquake belt is steeper in the west, and gradually becomes more gentle from west to east, dipping to the south slightly. The main shock( MS6.4) and the foreshocks including the MS5.7 event occurred along the NNW earthquake belt. A large number of aftershocks occurred along the near EW earthquake belt, and two aftershocks above M5 occurred at the eastern side of the EW earthquake belt. The aftershocks on the south side of the main shock are rare, perhaps affected by the hard blocks of the Tarim Basin. The aftershocks distribution clearly illuminates a near EW-striking structure, likely the extension of the NNW-striking fault activated during the initial sequence. The dominant depth of the earthquake sequence is between 10km and 20km, the focal depth of aftershocks along the near EW direction is gradually shallower from west to east. We determined the focal mechanism solutions of the MS≥5.0 earthquakes by the CAP method. The results of focal mechanism inversion show that the focal mechanism of the main shock and two aftershocks above MS5 are mainly thrusting, and the MS5.7 foreshock is mainly strike-slip. We also determined the moment tensor solution of the main shock using ISOLA method as a single-source. The focal mechanism solutions of the main shock obtained by the two methods are consistent. The moment tensor solution of the main shock has a large non-double couple component, which proves that the rupture process is very complex. By inversion of the main shock using ISOLA method as a multi-source, the main shock, which was reported as a single event, is instead composed of two sub-events, a strike-slip rupture and the second thrust rupture. Within 4s, a strike-slip earthquake triggered a second large rupture on a thrust fault. The first rupture is consistent with the mechanism of the MS5.7 foreshock, and the second rupture is consistent with thrust-faulting mechanisms in the ensuing aftershock sequence. By analyzing the data of spatial distribution and focal mechanism of the earthquake sequence, it is speculated that the Jiashi MS6.4 earthquake occurred in the middle and lower crust below the detachment layer of the Kalpin thrust tectonic zone. The occurrence of the main shock is caused by the joint action of the two faults, the NNW-striking fault with a high dip angle and the near EW-striking fault dipping south. The MS6.4 rupture initiated the adjacent previous NNW-striking rupture of the MS5.7 event, extending the earlier rupture both to the NNW and EW directions. The MS6.4 earthquake is the result of the interaction between the two blocks, the south Tianshan Mountains and the Tarim Block.

Keyword: Jiashi MS6.4 earthquake; earthquake relocation; focal mechanism; seismogenic fault
0 引言

据中国地震台网测定, 北京时间2020年1月18日新疆喀什地区伽师县发生MS5.7地震, 1月19日21时27分再次发生MS6.4地震, MS6.4地震的震中位于(39.8° N, 77.2° E), 震源深度为16km, 这是继2003年MS6.8地震之后伽师地区发生的又一次6级地震。 本次地震的最大烈度为Ⅷ 度, 地震造成1人死亡, 2人轻伤, 部分房屋在地震中开裂和倒塌等。 主震后1h和2月21日分别发生了MS5.2和MS5.1强余震。

本次伽师MS6.4地震发生在塔里木盆地和南天山交会的区域, 该区域受南天山、 塔里木和帕米尔构造块体共同作用。 塔里木块体与南天山相向对冲, 在山前形成一系列N倾、 EW向延伸的逆冲断层和数排背斜褶皱组成呈弧形的柯坪塔格推覆构造, 地质构造复杂, 新构造运动强烈(徐锡伟等, 2006)。

如图 1 所示, 该区域内的主要构造有EW走向的柯坪断裂(F1)、 奥兹格尔他乌断裂(F2)、 八盘水磨断裂(F4)、 NE向的迈丹断裂(F7)以及NNW向的皮羌断裂(F3)。 柯坪断裂和奥兹格尔他乌断裂为2条近EW向的平行断裂, 南侧的柯坪断裂为倾向NW、 倾角约30° 的逆冲断裂, 北侧的奥兹格尔他乌断裂为倾向NNE、 倾角约50° 的逆冲断裂(何文渊等, 2002)。 以右旋走滑的皮羌断裂为界, 在皮羌断裂以东断裂多呈走滑性质, 走向以NE为主; 皮羌断裂以西断裂以逆冲为主, 多为EW走向(沈军等, 2006)。 杨晓平等(2006)的研究表明柯坪塔格推覆褶皱带是一种薄皮构造, 其滑脱面具有东浅西深的特点, 皮羌断裂西侧滑脱面的深度为7~9km, 东侧约为5km。 Li等(2020)的研究表明柯坪塔格断层分段变形从西向东急剧减少, 皮羌断裂以西的缩短率为2.5~2.7mm/a, 以东约为0.3mm/a。 该区域现今的构造格局是新生代时期印度板块与欧亚板块碰撞的产物, 其强震活动多受控于较为年轻的柯坪断裂。 独特的地质构造和应力使得该区域成为中强地震(群)频发区域, 历史上发生过数次中强地震, 仅1997年1— 4月就在伽师地区很小的范围内发生了7次6级强震群活动, 1998年伽师地区发生了2次6级地震和阿图什6.0级地震, 2003年2月24日巴楚-伽师再次发生MS6.8地震, 前人对这几次强震开展了一系列研究(刘启元等, 2000; 周仕勇等, 2001; 郭飙等, 2002; 王卫民等, 2005; 黄媛等, 2006; 赵翠萍等, 2008a, b)。 与前几次强震发生在盆山交界区域的塔里木盆地一侧不同, 此次MS6.4地震及其序列发生在塔里木盆地北侧的柯坪塔格山前, 暂未发现伽师MS6.4地震同震导致的地表破裂, 此次伽师地震的发震断层究竟与哪条断裂有关等科学问题目前尚无定论。

图 1 区域构造与历史地震分布
五角星为伽师MS6.4地震震中, 红色圆圈为1996年至今的伽师及周边6级以上地震震中, 三角形为地震台站, 其中L6505和L6513为主震后增加的流动台站, 黑色圆圈为主要城市。 F1柯坪断裂; F2奥兹格尔他乌断层; F3皮羌断裂; F4八盘水磨断裂; F5衣木干它乌断裂; F6喀拉铁克断裂; F7迈丹断裂; F8塔拉斯-费尔干纳断裂; F9西昆仑北缘断裂
Fig. 1 The regional tectonic map and the distribution of historical earthquakes.

本文使用新疆地震台网固定台站和震后布设的流动台站观测数据, 对伽师MS6.4地震序列进行重新定位, 反演了主震及3次5级以上地震的震源机制解, 研究了伽师地震序列的活动特征和发震构造。

1 数据和方法

由于伽师震源区的监测能力相对较弱, MS6.4地震发生后, 新疆地震局在震中附近架设了2个流动台(L6505和L6513), 结合区域地震观测台站(图 1)分析了2020年1月18日— 8月31日的观测数据, 共得到3 356个地震事件。 伽师地震序列的M-T图(图 2a)显示余震活动衰减较慢。 为获得较为可靠的定位结果, 本文选用观测台站数量> 7的地震事件参加重定位, 经筛选得到1 990个地震事件的观测资料。 选择震中距< 250km的到时数据, 并剔除离散走时数据。 图 2b是经挑选后的P波和S波观测走时数据, 其中包括P波数据17 752条, S波数据16 519条。 参与定位的台站共有15个, 对地震序列有较好的空间包围(图 1)。 采用区域地震台网的宽频带波形资料反演震源机制, 选择震中距< 350km的观测台站且有效记录4次5级以上地震事件(含MS6.4主震)的波形资料。 本文进行重定位和震源机制解反演所使用的一维速度模型如图 2c所示, 该模型综合了前人的上地壳速度结构反演结果(刘启元等, 2000; 赵翠萍等, 2008a)。

图 2 a 伽师地震序列的M-T图(2020年1月1日— 8月31日); b 经挑选后的P波和S波观测走时曲线; c 本文使用的速度模型Fig. 2 The magnitude-time plot of Jiashi earthquake sequence from January 1 to August 31, 2020(a), the selected travel time curves of P and S wave(b)and velocity model used in the study(c).

本文采用双差定位方法(Waldhauser et al., 2000)对伽师地震序列进行重定位, 该方法使用相对走时残差反演地震震源位置, 能够消除震源到台站的传播路径效应, 减少地壳速度模型不准确对定位结果的影响。 双差定位方法在国内外中强地震序列重定位研究中有着广泛的应用(黄媛等, 2008; 张广伟等, 2014; 房立华等, 2018)。 余震序列的高精度定位结果能较好地勾画出地震的发震构造。

本文采用CAP方法(Zhu et al., 1996)反演伽师地震序列中4次5级以上地震的震源机制解。 CAP方法假定震源模型为双力偶模型, 将观测波形分为长周期体波和面波2部分, 分别与理论波形进行拟合从而得到最佳震源机制解(Zhao et al., 1994; Zhu et al., 1996), 其中理论地震图使用频率-波数法计算(Saikia, 1994; Zhu et al., 1996)。 CAP方法已经在汶川、 芦山等地震序列震源机制研究中得到广泛应用(郑勇等, 2009; 吕坚等, 2013; Han et al., 2014; 罗钧等, 2015; 罗艳等, 2015)。

考虑到前震和余震分布所反映的伽师MS6.4地震震源过程和发震构造的复杂性, 本文采用ISOLA方法(Zahradnik et al., 2005, 2014)进一步分析MS6.4地震的震源过程。 该方法将此次地震分解为空间和时间上的多个子事件, 利用格点搜索求取达到最佳波形拟合的子事件发生时间、 位置和震源机制解, 从而探讨伽师震源区复杂的构造环境和断层展布对地震破裂过程的影响。 前人使用ISOLA方法对复杂地震事件进行反演, 并取得了较好的结果(Hicks et al., 2015; Liu et al., 2020), 其中格林函数使用离散波数方法计算得到(Bouchon, 1981)。

2 重定位和地震活动的时空过程

选择观测台站数量> 7的1 990个地震事件开展重定位研究。 图 3a为重定位前伽师地震序列的平面分布图; 利用双差定位方法进行重定位后最终得到1 460个地震的重定位结果, 图 3b为重定位后伽师地震序列的平面分布图, 其EW向、 SN向及深度的平均定位偏差分别为1.14km、 1.43km和1.90km; 重新定位后的地震走时平均均方根残差较重定位前显著降低, 由0.838s下降为0.275s。 本次伽师MS6.4主震的重定位结果为(39.835° N, 77.148° E), 震源深度为14.9km。 由图 3b所示的重定位结果可以看出, 地震呈现2个条带分布, NNW向地震条带相对较短, 约为20km; 近EW向地震条带为优势分布区, 延展长度约为40km。 主震和MS5.7前震的位置相近, 分布在NNW向条带的南侧, NNW向主震南侧仅有少量地震活动, 2次5级以上余震分布在近EW向条带的东侧。 由图 3e— h所示的深度剖面可以看出, 地震主要分布在5~25km深度范围, 优势分布在10~20km之间。 由AA'剖面可以看出, 地震深度由西到东逐渐变浅。 从BB'剖面可以看出南侧主震及余震深度集中在约15km处, 靠近北侧的地震分布在10~20km深度。 穿过主震震中垂直于BB'的CC'剖面显示NNW向的地震活动范围较窄, 发震构造的倾角较大, 倾向W。 由垂直于AA'剖面的DD'、 EE'和GG'深度剖面可以看出, 近EW向的地震活动构造略微S倾。 西侧DD'剖面附近的地震活动分布较窄、 倾角较陡, 由西向东深度变浅, 地震分布也逐渐变宽、 倾角变缓。

图 3 重定位前后伽师地震序列分布及深度剖面
a 重定位前伽师地震序列的平面分布; b 重定位后伽师地震序列的平面分布; c 重定位后伽师地震序列AA'剖面图; d 重定位后伽师地震序列BB'剖面图; e 重定位后伽师地震序列CC'剖面图; f 重定位后伽师地震序列DD'剖面图; g 重定位后伽师地震序列EE'剖面图; h 重定位后伽师地震序列GG'剖面图
Fig. 3 Distribution of initial and relocated earthquake sequences and the vertical cross section of the relocated earthquake sequence.

为分析主震前后地震活动的时空变化特征, 本文按照不同时间段展示地震序列活动的分布情况。 如图 4a、 b可见, 主震发生前48h, 2020年1月18日MS5.7前震序列主要在NNW向发生, 余震很少且沿NNW向形成尺度较小的条带; 直至主震发生, 近EW向没有明显的地震活动。 主震发生后48h, 大量的余震活动形成近EW向条带(图 4c, d), 且震源深度由西向东逐渐变浅, 并于主震后1h在近EW向破裂的东侧发生了MS5.2余震。 主震发生后48~96h(图 4e, f), 余震活动较之前有明显衰减。 图 4e显示这一时期的余震活动主要近EW向展布, NNW向也有少量的地震活动。 根据地震序列活动时空分布特征可以推测MS5.7前震活动主要沿NNW向破裂, 当主震发生后, 破裂向N延展一段距离后沿近EW向展布。

图 4 不同时间段的地震平面及深度剖面图
a 主震前48h地震活动的平面分布图; b 主震前48h地震活动的深度剖面图; c 主震后48h地震活动的平面分布图; d 主震后48h地震活动的深度剖面图; e 主震后48~96h地震活动的平面分布图; f 主震后48~96h地震活动的深度剖面图
Fig. 4 Relocated earthquakes in map view and along the section AA' at different periods of time.

结合重定位结果所展示的地震时空分布特征推测, MS5.7前震的破裂面为NNW向, 断层面倾角较陡, 倾向W; MS6.4主震的破裂面由早期的NNW向转为近EW向, 近EW向断层面西侧的倾角较陡, 向E逐渐变缓, 断层面略微S倾。 推测MS6.4主震的发生受NNW向W倾断裂和近EW向S倾高角度断裂共同控制。

3 震源机制解

基于区域宽频带台站的三分量波形, 本文采用CAP方法反演MS6.4主震的震源机制解。 由图 5 所示的波形拟合图和深度搜索图可以看出观测波形和理论波形拟合较好, 说明主震的震源机制解较可靠, 主震震源机制解的节面I和节面Ⅱ 的走向、 倾角和滑动角分别为84° 、 67° 、 111° 和220° 、 31° 、 50° , 矩震级 MW 5.87, 最佳拟合深度为14.5km, 具有明显的逆冲性质。 本文反演了另外3次5级以上地震的震源机制解, 详细结果见表1和图 6。 除MS5.7前震为走滑事件外, 另外2次余震的震源机制均为逆冲性质, 与主震较为一致。 可以看出, MS5.7前震与主震明显不同, 主震与MS5.2、 MS5.1余震的震源机制解一致, MS5.2、 MS5.1地震是发生在近EW向构造上的强余震。 考虑余震分布的高倾角特征, 推测节面I为MS6.4地震的同震破裂面。

图 5 a 伽师MS6.4主震的震源机制解及波形拟合图; b 不同深度的震源机制解结果
波形拟合图中黑线为观测波形, 红线为理论波形; 左侧大写字母为台站名; 波形下方的数字为理论波形、 观测波形的相对移动时间(s)和拟合相关系数(%)
Fig. 5 Focal mechanism and waveform fitting diagram of the Jiashi MS6.4 earthquake(a) and the focal mechanism solutions at each source depth(b).

表1 MS≥ 5.0地震的震源机制解 Table 1 Focal mechanisms of the MS≥ 5.0 earthquakes

图 6 伽师地震序列MS≥ 5.0地震的震源机制解
圆圈为3级以上地震震中, 颜色代表距主震的时间
Fig. 6 Focal mechanisms of the MS≥ 5.0 earthquakes in Jiashi area.

MS5.7前震和MS6.4主震的位置相近, 均分布在NNW方向的地震条带, 但二者的震源机制解却明显不同。 重定位结果显示, MS6.4主震形成了2个走向分别为NNW和EW向的余震条带, MS6.4主震位于NNW条带, 与EW向展布的余震也有一定距离。 为分析MS6.4主震复杂的发震构造和破裂过程, 本文进一步使用ISOLA方法对MS6.4主震开展矩张量解反演研究。 使用单源模型(Single-point source model)对MS6.4主震进行全波形反演(频带范围为0.02~0.04Hz), 得到主震的矩震级为MW5.86, 震源机制解节面I和节面Ⅱ 的走向、 倾角和滑动角分别为72° 、 68° 、 113° 和203° 、 31° 、 45° 。 该结果与CAP方法反演的主震震源机制解基本一致。 主震矩张量解中的双力偶成分(DC)和补偿线性偶极子(CLVD)分别占50.1%和-49.9%, 显示MS6.4主震具有显著的非双力偶成分, 这表明MS6.4主震的震源破裂过程较为复杂。

以主震位置为中心将震源区划分为9× 11的步长为5km的水平网格, 固定震源深度为15km, 使用多源模型(Multi-point source model)进行全波形反演(频带范围为0.02~0.08Hz)。 通过格点搜索得到2个相关系数较高、 时间相距约4s的子事件(图 7, 8), 表明MS6.4主震可能由2个子事件组成。 图 7 和表2显示子事件1位于起始破裂点北侧25km处, 震级MW5.67, 该事件具有走滑性质, 2个节面均近直立, 其中走向近SN向的节面与NNW向余震的走向一致, 也与MS5.7前震的震源机制一致; 子事件2位于起始破裂点北侧15km、 东侧10km处, 震级MW5.80, 该事件具有明显逆冲性质, 与2次5级强余震的震源机制一致, 2个节面的走向均为近EW向, 节面Ⅰ 基本直立, 节面Ⅱ 的倾角较缓。 2个子事件合成产生的波形拟合情况见图 8, 可见多数台站的拟合情况较好, KSZ和ATS等台的部分波形拟合效果不理想可能与南天山复杂的地下结构有关。

图 7 基于多源模型的MS6.4主震矩张量解反演结果
a 子事件1; b 子事件2。 五角星为主震震中, 等值线为波形拟合相关系数, 最大的沙滩球为最优解
Fig. 7 Moment tensor inversion of the MS6.4 earthquake based on the multi-point-source model.

图 8 最优解的观测波形(黑)和理论波形(红)拟合结果Fig. 8 Observed(black)and synthetic(red)waveforms for the optimum solution.

表2 2个子事件的震源机制解 Table 2 Focal mechanisms of two subevents
4 讨论和结论

本文采用双差定位方法获得了伽师MS6.4地震序列1 460个事件的重定位结果。 结果显示, 伽师MS6.4主震的重定位震中位于(39.835° N, 77.148° E), 震源深度为14.9km。 伽师地震序列活动呈现2个优势方向, 其中NNW向条带延展约20km, 近EW向条带延展约40km。 MS5.7前震至主震发生前的地震活动主要沿NNW向分布, 主震发生后在NNW和EW方向同时产生大量余震活动。 NNW向断层面的倾角较陡, 倾向W; 近EW向断层面西侧的倾角较陡, 向E逐渐变缓, 断层面略微S倾。 震源深度主要集中在5~25km, 优势分布于10~20km, NNW向地震分布较深, 在EW走向由西向东呈现震源深度逐渐变浅的特征。

基于CAP方法和ISOLA方法反演的伽师MS6.4主震的震源机制解一致, 揭示此次地震是以逆冲错动为主的机制。 主震的矩张量解同时显示主震具有较大的非双力偶成分, 表明主震的破裂过程较为复杂。 我们进一步对主震开展基于多源模型的矩张量反演, 结果揭示了主震可能是由时间相差约4s的2个子事件组成, 第1个子事件为走滑性质, 机制与MS5.7前震一致, 其近SN向节面与余震活动形成的NNW向条带吻合, 矩震级MW5.67; 第2个子事件为逆冲性质, 与2次强余震的机制一致, 其节面与近EW向余震的优势走向一致, 矩震级MW5.80。 该破裂过程与主震发生前、 后的地震序列时空分布特征一致。

结合重定位结果和震源机制特征, 我们推测伽师MS6.4地震的发震构造由相交的2条断裂构成, 破裂始于NNW向的断裂, 在向N扩展的过程中与近EW向的断裂交会, 触发了近EW向断裂的活动。 2条断裂共同参与这次地震, 导致余震分布也呈现出2个优势方向。 由于震中附近并没有已知的NNW向断裂, 推测存在1条NNW向高倾角走滑性质的隐伏断裂, 可能是切穿柯坪逆冲推覆体— 中下地壳的基底断裂, 具有与皮羌断裂类似的 “ 类转换断层” 性质(徐锡伟等, 2006), 调节柯坪逆冲推覆体褶皱带的横向不均匀性。 震中附近近EW向的柯坪断裂为倾向NW、 倾角约30° 的逆冲断裂(何文渊等, 2002)。 由本文的伽师序列重定位结果可以看出, 近EW向逆冲事件的发震构造为S倾的高角度断裂, 与柯坪断裂的倾向相反; 柯坪断裂为滑脱面深度7~9km的逆冲推覆构造, 而主震重定位结果及矩心深度均为14~15km, 显示其孕震深度已到达柯坪逆冲推覆体滑脱面下的中下地壳, 可能与深部地壳中的基底断裂有关, 推测存在1条近EW向的S倾构造, 伽师MS6.4地震的初始破裂是在NNW向主压应力作用下发生的。

根据哈佛大学提供的1996— 2003年6级以上地震震源机制解(图 1), 塔里木盆地内的1997年伽师震群和1998年2次6级地震的震源机制解多为走滑或正断型, 而同样发生在盆地内的2003年的巴楚-伽师MS6.8地震为逆冲型, 表明该区域存在挤压与拉张并存的构造背景(胥颐等, 2000)。 周仕勇等(2001)根据震源精确定位结果推测1997年的伽师强震群可能存在NNW走向的发震构造; 赵翠萍(2006)根据1997年震群活动特征推测该强震群序列总体上呈现NE和NNW向的断层活动。 相对于1997年的强震群, 1998年2次强震向NE迁移, 余震分布呈2条平行的NE向条带展布。 2003年的巴楚-伽师MS6.8地震发生在柯坪塔格推覆构造南缘一条近EW走向尚未出露的逆断层上(徐锡伟等, 2006); 伽师北部1996年和1998年的2次阿图什6级地震靠近柯坪断裂, 为纯逆冲型, 断层面走向近EW, 与柯坪断裂的走向一致。 可见, 2020年1月19日伽师MS6.4地震与1997— 1998年伽师震群以及2003年巴楚-伽师MS6.8地震的发震断裂明显不同, 此次地震是由NNW向高角度W倾的隐伏断裂和近EW向S倾的次级断裂先后破裂所构成的。

深地震测深资料表明, 塔里木盆地的结晶基底界面埋深为10~12km, 而康拉德界面的埋深为22~25km, 这2个界面由塔里木盆地向南天山急剧变浅(张先康等, 2002), 中下地壳界面起伏不稳定。 本文的重定位结果表明, 伽师MS6.4地震的震源正是位于地壳结构不稳定的过渡带, 伽师地震序列的深度主要集中在10~20km, 地震活动深度上受到结晶基底界面和康拉德界面限制, 主要分布在两者之间, 表明地壳结构不稳定可能是该区域地震多发的影响因素。 塔里木块体向N插入南天山, 两者交会区域的地壳出现强烈的不均匀变形, 岩石物性和地震波速都存在着明显差异(胥颐等, 2000), 与主要断裂构成了复杂的耦合关系, 形成该区域中强地震频繁发生的构造环境。 本次伽师MS6.4地震的发生是南天山与塔里木两大块体作用的产物, 近EW向分布的余震序列震源深度逐渐变浅的原因可能是中下地壳界面急剧变浅导致的, NNW向主震南侧余震活动较少, 推测可能是受到坚硬的塔里木盆地物质阻挡所致。

致谢 本研究使用了新疆地震台网的观测数据、 Waldhauser的双差定位程序、 朱露培教授的CAP程序和Zahradní k的ISOLA程序; 审稿专家和主编对本文提出了宝贵的修改意见。 在此一并表示感谢!

参考文献
[1] 房立华, 吴建平, 苏金蓉, . 2018. 四川九寨沟 MS7. 0地震主震及其余震序列精定位[J]. 科学通报, 63(7): 649662.
FANG Li-hua, WU Jian-ping, SU Jin-rong, et al. 2018. Relocation of mainshock and aftershock sequence of the MS7. 0 Sichuan Jiuzhaigou earthquake[J]. Chinese Science Bulletin, 63(7): 649662(in Chinese). [本文引用:1]
[2] 郭飙, 刘启元, 陈九辉, . 2002. 新疆伽师强震区余震序列的地震台阵定位[J]. 地震地质, 24(2): 199207.
GUO Biao, LIU Qi-yuan, CHEN Jiu-hui, et al. 2002. Relocation of earthquake sequences using joint hypocenter determination method: Portable seismic array study in Jiashi region, Xinjiang[J]. Seismology and Geology, 24(2): 199207(in Chinese). [本文引用:1]
[3] 何文渊, 李江海, 钱祥麟, . 2002. 塔里木盆地柯坪断隆断裂构造分析[J]. 中国地质, 29(1): 3743.
HE Wen-yuan, LI Jiang-hai, QIAN Xiang-lin, et al. 2002. Analysis of fault structures in the Kalpin fault uplift, Tarim Basin[J]. Geology in China, 29(1): 3743(in Chinese). [本文引用:2]
[4] 黄媛, 吴建平, 张天中, . 2008. 汶川8. 0级大地震及其余震序列重定位研究[J]. 中国科学(D辑), 38(10): 12421249.
HUANG Yuan, WU Jian-ping, ZHANG Tian-zhong, et al. 2008. Relocation of the M8. 0 Wenchuan earthquake and its aftershock sequence[J]. Science in China(Ser D), 38(10): 12421249(in Chinese). [本文引用:2]
[5] 黄嫒, 杨建思, 张天中. 2006. 2003年新疆巴楚-伽师地震序列的双差法重新定位研究[J]. 地球物理学报, 49(1): 162169.
HUANG Yuan, YANG Jian-si, ZHANG Tian-zhong. 2006. Relocation of the Bachu-Jiashi, Xinjiang earthquake sequence in 2003 using the double difference location algorithm[J]. Journal of Geophysics, 49(1): 162169(in Chinese). [本文引用:1]
[6] 刘启元, 陈九辉, 李顺成, . 2000. 新疆伽师强震群区三维地壳上地幔S波速度结构及其地震成因的探讨[J]. 地球物理学报, 43(3): 356365.
LIU Qi-yuan, CHEN Jiu-hui, LI Shun-cheng, et al. 2000. Passive seismic experiment in Xinjiang Jiashi strong earthquake region and discussion on its seismic genesis[J]. Chinese Journal of Geophysics, 43(3): 356365(in Chinese). [本文引用:2]
[7] 罗钧, 赵翠萍, 周连庆. 2015. 2013年8月香格里拉德钦-得荣 MS5. 9地震序列震源机制与应力场特征[J]. 地球物理学报, 58(2): 424435.
LUO Jun, ZHAO Cui-ping, ZHOU Lian-qing. 2015. Focal mechanisms and stress field of the Shangri-La Deqen, Yunnan Province-Derong, Sichuan Province MS5. 9 earthquake sequence in August, 2013[J]. Chinese Journal of Geophysics, 58(2): 424435(in Chinese). [本文引用:1]
[8] 罗艳, 赵里, 曾祥方, . 2015. 芦山地震序列震源机制及其构造应力场空间变化[J]. 中国科学(D辑), 45(4): 538550.
LUO Yan, ZHAO Li, ZENG Xiang-fang, et al. 2015. Focal mechanisms of the Lushan earthquake sequence and spatial variation of the stress field[J]. Science in China(Ser D), 45(4): 538550(in Chinese). [本文引用:1]
[9] 吕坚, 王晓山, 苏金蓉, . 2013. 芦山7. 0级地震序列的震源位置与震源机制解特征[J]. 地球物理学报, 56(5): 17531763.
Jian, WANG Xiao-shan, SU Jin-rong, et al. 2013. Hypocentral location and source mechanism of the MS7. 0 Lushan earthquake sequence[J]. Chinese Journal of Geophysics, 56(5): 17531763(in Chinese). [本文引用:1]
[10] 沈军, 陈建波, 王翠, . 2006. 2003年2月24日新疆巴楚-伽师6. 8级地震发震构造[J]. 地震地质, 28(2): 205212.
SHEN Jun, CHEN Jian-bo, WANG Cui, et al. 2006. The seismogenic tectonics of the MS6. 8 Bachu-Jiashi, Xinjiang earthquake on Feb. 24, 2003[J]. Seismology and Geology, 28(2): 205212(in Chinese). [本文引用:1]
[11] 王卫民, 李丽, 赵连锋, . 2005. 2003年2月24日新疆伽师 MS6. 5地震震源破裂过程研究[J]. 地球物理学报, 48(2): 343351.
WANG Wei-min, LI Li, ZHAO Lian-feng, et al. 2005. Rupture process of Jiashi, Xinjiang earthquake( MS6. 5)of Feb. 24, 2003[J]. Chinese Journal of Geophysics, 48(2): 343351(in Chinese). [本文引用:1]
[12] 徐锡伟, 张先康, 冉勇康, . 2006. 南天山地区巴楚-伽师地震( MS6. 8)发震构造初步研究[J]. 地震地质, 28(2): 161178.
XU Xi-wei, ZHANG Xian-kang, RAN Yong-kang, et al. 2006. The preliminary study on seismotectonics of the 2003 AD Bachu-Jiashi earthquake( MS6. 8), southern Tianshan[J]. Seismology and Geology, 28(2): 161178(in Chinese). [本文引用:3]
[13] 胥颐, 刘福田, 刘建华, . 2000. 天山地震带的地壳结构与强震构造环境[J]. 地球物理学报, 43(2): 184193.
XU Yi, LIU Fu-tian, LIU Jian-hua, et al. 2000. Crustal structure and tectonic environment of strong earthquakes in the Tianshan earthquake belt[J]. Chinese Journal of Geophysics, 43(2): 184193(in Chinese). [本文引用:2]
[14] 杨晓平, 冉勇康, 宋方敏, . 2006. 西南天山柯坪逆冲推覆构造带的地壳缩短分析[J]. 地震地质, 28(2): 194204. [本文引用:1]
[15] YANG Xiao-ping, RAN Yong-kang, SONG Fang-min, et al. 2006. The analysis for crust shortening of Kalpin thrust tectonic zone, south-western Tianshan, Xinjiang, China[J]. Seismology and Geology, 28(2): 194204(in Chinese). [本文引用:1]
[16] 张广伟, 雷建设, 梁姗姗, . 2014. 2014年8月3日云南鲁甸 MS6. 5地震序列重定位与震源机制研究[J]. 地球物理学报, 57(9): 30183027.
ZHANG Guang-wei, LEI Jian-she, LIANG Shan-shan, et al. 2014. Relocations and focal mechanism solutions of the 3 August 2014 Ludian, Yunnan MS6. 5 earthquake sequence[J]. Chinese Journal of Geophysics, 57(9): 30183027(in Chinese). [本文引用:1]
[17] 张先康, 赵金仁, 张成科, . 2002. 帕米尔东北侧地壳结构研究[J]. 地球物理学报, 45(5): 665671.
ZHANG Xian-kang, ZHAO Jin-ren, ZHANG Cheng-ke, et al. 2002. Crustal structure at the northeast of the Pamirs[J]. Chinese Journal of Geophysics, 45(5): 665671(in Chinese). [本文引用:1]
[18] 赵翠萍. 2006. 1997—2003年新疆伽师震源区特征的地震学方法研究[D]. 北京: 中国地震局地球物理研究所.
ZHAO Cui-ping. 2006. Seismological studies on the characteristics of Jiashi source region from 1997 to 2003 [D]. Institute of Geophysics, China Earthquake Administration, Beijing(in Chinese). [本文引用:1]
[19] 赵翠萍, 陈章立, 郑斯华, . 2008 a. 伽师震源区中等强度地震矩张量反演及其应力场特征[J]. 地球物理学报, 51(3): 782792.
ZHAO Cui-ping, CHEN Zhang-li, ZHENG Si-hua, et al. 2008a. Moment inversion of moderate earthquakes and the locally perturbed stress field in the Jiashi source region[J]. Chinese Journal of Geophysics, 51(3): 782792(in Chinese). [本文引用:2]
[20] 赵翠萍, 陈章立, 郑斯华. 2008 b. 1998—2003年伽师3次不同类型 MS6地震震源破裂过程及短期内余震活动特征[J]. 地球物理学报, 51(4): 10931102.
ZHAO Cui-ping, CHEN Zhang-li, ZHENG Si-hua. 2008b. Source rupture process of 3 Jiashi MS6 events(1998—2003) and its correlation with the aftershock activity[J]. Chinese Journal of Geophysics, 51(4): 10931102(in Chinese). [本文引用:1]
[21] 郑勇, 马宏生, 吕坚, . 2009. 汶川地震强余震( MS≥5. 6)的震源机制解及其与发震构造的关系[J]. 中国科学(D辑), 39(4): 413426.
ZHENG Yong, MA Hong-sheng, Jian, et al. 2009. Source mechanism of strong aftershock of( MS≥5. 6)of Wenchuan earthquake and the implication for seismotectonic[J]. Science in China(Ser D), 39(4): 413426(in Chinese). [本文引用:1]
[22] 周仕勇, 许忠淮, 陈晓非. 2001. 伽强震群震源特征及震源机制力学成因分析[J]. 地球物理学报, 44(5): 654662.
ZHOU Shi-yong, XU Zhong-huai, CHEN Xiao-fei. 2001. Analysis on the source characteristics of the 1997 Jiashi swarm, western China[J]. Chinese Journal of Geophysics, 44(5): 654662(in Chinese). [本文引用:1]
[23] Bouchon M. 1981. A simple method to calculate Green’s functions for elastic layered media[J]. Bulletin of the Seismological Society of America, 71(4): 959971. [本文引用:1]
[24] Han L, Zeng X, Jiang C, et al. 2014. Focal mechanisms of the 2013 MW6. 6 Lushan, China earthquake and high-resolution aftershock relocations[J]. Seismological Research Letters, 85(1): 814. [本文引用:1]
[25] Hicks S P, Rietbrock A. 2015. Seismic slip on an upper-plate normal fault during a large subduction megathrust rupture[J]. Nature Geoscience, 8(12): 955960. [本文引用:1]
[26] Li A, Ran Y K, Francisco G, et al. 2020. Segmentation of the Kepingtage thrust fault based on paleoearthquake ruptures, southwestern Tianshan, China[J]. Natural Hazards: Journal of the International Society for the Prevention and Mitigation of Natural Hazards, 103(1): 13851406. [本文引用:1]
[27] Liu J, Zahradnik J. 2020. The 2019 MW5. 7 Changning earthquake, Sichuan Basin, China: A shallow doublet with different faulting styles[J]. Geophysical Research Letter, 47(4): 19. doi: 10.1029/2019GL085408. [本文引用:1]
[28] Saikia C K. 1994. Modified frequency-wavenumber algorithm for regional seismograms using Filon’s quadrature: Modelling of Lg waves in eastern North America[J]. Geophysical Journal International, 118(1): 142158. [本文引用:1]
[29] Waldhauser F, Ellsworth W L. 2000. A double-difference earthquake location algorithm: Method and application to the Northern Hayward Fault, California[J]. Bulletin of the Seismological Society of America, 90(6): 13531368. [本文引用:1]
[30] Zahradnik J, Serpetsidaki A, Sokos E, et al. 2005. Iterative deconvolution of regional waveforms and a double-event interpretation of the 2003 Lefkada earthquake, Greece[J]. Bulletin of the Seismological Society of America, 95(1): 159172. [本文引用:1]
[31] Zahradnik J, Sokos E. 2014. The MW7. 1 Van Eastern Turkey, earthquake 2011: Two-point source modelling by iterative deconvolution and non-negative least squares[J]. Geophysical Journal International, 196(1): 522538. [本文引用:1]
[32] Zhao L S, Helmberger D V. 1994. Source estimation from broadband regional seismograms[J]. Bulletin of the Seismological Society of America, 84(1): 91104. [本文引用:1]
[33] Zhu L P, Helmberger D V. 1996. Advancement in source estimation techniques using broadband regional seismograms[J]. Bulletin of the Seismological Society of America, 86(5): 16341641. [本文引用:3]