引用本文
张文婷, 季灵运, 朱良玉, 蒋锋云, 徐晓雪. 南天山前陆盆地的一次典型逆冲破裂事件——2020年新疆伽师6.4级地震[J]. 地震地质, 2021,43(2): 394-409.
ZHANG Wen-ting, JI Ling-yun, ZHU Liang-yu, JIANG Feng-yun, XU Xiao-xue. A TYPICAL THRUST RUPTURE EVENT OCCURRING IN THE FORELAND BASIN OF THE SOUTHERN TIANSHAN: THE 2020 XINJIANG JIASHI MS6.4 EARTHQUAKE [J]. SEISMOLOGY AND GEOLOGY, 2021,43(2): 394-409.  
DOI:10.3969/j.issn.0253-4967.2021.02.009
Permissions
Copyright©2021, 《地震地质》编辑部
《地震地质》编辑部 所有
南天山前陆盆地的一次典型逆冲破裂事件——2020年新疆伽师6.4级地震
张文婷1,2), 季灵运1,2),*, 朱良玉1,2), 蒋锋云1,2), 徐晓雪1,2)
1)中国地震局第二监测中心, 西安 710054
2)国家遥感中心地质灾害研究部, 北京 100036
*通讯作者: 季灵运, 男, 1982年生, 博士, 研究员, 主要从事构造大地测量学研究, 电话: 029-85506645, E-mail: dinsar010@163.com。

〔作者简介〕 张文婷, 女, 1993年生, 2018年于长安大学获大地测量专业硕士学位, 工程师, 现主要研究方向为构造大地测量学, 电话: 029-85506645, E-mail: wt_zhang93@163.com。

收稿日期: 2020-09-30 修回日期: 2021-02-19
基金项目: 国家自然科学基金(41904007, 41604015)和国家重点研发计划项目(2017YFC1500102)共同资助
摘要

2020年1月19日, 新疆维吾尔自治区伽师县发生6.4级地震, 震中位于南天山西段的柯坪推覆体附近。 文中基于Sentinel-1 SAR卫星升、降轨数据获取了此次地震的InSAR同震形变场, InSAR结果表明, 地震形变场的EW向跨度约为40km, SN向跨度约为20km, 升轨结果LOS向最大抬升量为5.9cm, 最大沉降量为3.7cm; 降轨结果LOS向最大抬升量约为6.4cm, 最大沉降量约为2cm。 以InSAR观测值为约束, 首先采用均匀滑动模型对发震断层参数进行反演, 为提高计算效率, 在反演前利用四叉树分割方法对原始的高分辨率InSAR观测值进行降采样; 然后在均匀滑动模型反演结果的基础上, 引入基于约束最小二乘的分布式滑动反演方法, 计算出此次地震发震断层面的精细滑动分布。 结果表明, 此次地震的累计地震矩为1.73×1018N·m, 合矩震级 MW6.1, 其发震断层为柯坪断裂。 对区域GPS面应变率、 构造动力背景和区域深浅部构造特征进行了综合分析, 认为此次地震是继1997—2003年伽师6~7级地震震群后发生在南天山强烈逆冲推覆构造环境中的又一次典型的逆冲破裂事件。 这些强震是南天山前陆盆地内不同级别、 不同性质的断裂相继破裂的结果, 共同受控于南天山的逆冲推覆作用。

关键词: 伽师地震; InSAR; 震源机制; 前陆盆地; 南天山
中图分类号:P315.2 文献标志码:A 文章编号:0253-4967(2021)02-0394-16
A TYPICAL THRUST RUPTURE EVENT OCCURRING IN THE FORELAND BASIN OF THE SOUTHERN TIANSHAN: THE 2020 XINJIANG JIASHI MS6.4 EARTHQUAKE
ZHANG Wen-ting1,2), JI Ling-yun1,2), ZHU Liang-yu1,2), JIANG Feng-yun1,2), XU Xiao-xue1,2)
1)The Second Monitoring and Application Center, CEA, Xian 710054, China
2)Department of Geological Hazards, National Remote Sensing Center of China, Beijing 100036, China
Abstract

A MS6.4 earthquake occurred on January 19th, 2020 at Jiashi, Xinjiang, this earthquake is another strong earthquake since the Jiashi MS6~7 earthquake swarm events from 1997 to 2003, and the epicenter was located near the Kalpin nappe in the western part of southern Tianshan. The Kaplin nappe is located in front of southern Tianshan Mountains, which is a thin skinned thrust belt composed of a series of nearly NE-SW thrust nappes under the strong and sustained regenerative orogeny in the Tianshan area. There are some differences in focal positions and fault parameters given by different institutions, therefore in this paper, high resolution InSAR coseismic deformation fields were obtained based on the ascending and descending tracks of Sentinel-1 SAR images to obtain the focal mechanism. The 30m resolution SRTM DEM data is chosen as the external DEM to eliminate the phases caused by topography, the robust Goldstein filtering is applied for phase smoothing, and the Delaunay minimum cost flow method is used for phase unwrapping. The variation range of interference fringes shows that the east-west span of the earthquake deformation field is about 40km, and that of the north-south direction is about 20km, the displacement results show that the maximum uplift displacement is 5.9cm and the maximum subsidence is 3.7cm along the LOS direction of the ascending data, the maximum uplift displacement is 6.4cm and the maximum subsidence is 2cm along the LOS direction of the descending data. And then the InSAR-derived deformation fields are used to obtain the seismogenic mechanism of this earthquake, and to improve the computational efficiency, the quadtree segmentation method is used to desample the original high-resolution InSAR observations before inversion. The coseismic slip distribution of the causative fault was inversed using a uniform sliding inversion method based on a Bayesian approach, and then the fine slip distribution of the fault plane of Jiashi earthquake was inversed using the distributed slip inversion method based on the constrained least squares. It should be noted that the fault plane is set as the shovel shape according to the geometric relationship between the seismogenic fault parameters inverted by uniform sliding and the exposed position of the Kapling Fault on the surface during the distributed slip inversion. According to the difference between the observed and simulated values, it can be seen that the residual error of the inversion model is small, indicating the reliability of the inversion result. The final result shows that the epicenter is located at 39.9°N, 77.28°E and the strike and dip angle of the seismogenic fault is 276° and 10.7°, respectively, the maximum dip slip and strike slip of fault plane is about 0.29m and 0.03m, respectively, which are located at the depth of about 5km underground. The cumulative coseismic moment is 1.73×1018N·m from InSAR inversion, which is equal to the moment magnitude of MW6.1 and the Kalpin Fault is supposed to be the causative fault. Then, regional GPS-derived surface strain rate, tectonic dynamic background, and regional deep and shallow structures were comprehensively analyzed. The results show that the Jiashi MS6.4 earthquake is a typical thrust event that occurred in the thrust nappe of the southern Tianshan. The 2020 Jiashi event and the 1997—2003 Jiashi M6~7 earthquakes swarm are the results of rupture of many faults with different scales and properties. And these events are all controlled by the thrust nappe of southern Tianshan.

Keyword: Jiashi earthquake; InSAR; focal mechanism; foreland basin; southern Tianshan
0 引言

据中国地震台网测定, 2020年1月19日21时27分新疆喀什地区伽师县(39.83° N, 77.21° E)发生了6.4级地震, 震中位于南天山西段的柯坪推覆体附近(图 1), 距离伽师县56km, 震源深度约为16km(孟令媛等, 2020; 聂晓红等, 2020), 此次地震是自1997— 2003年伽师6~7级震群事件以来发生的又一次强烈地震。 多家机构对此次地震的震源机制解进行了解算, 得到的结果均显示此次地震为一次逆冲型地震(表1)。 目前已有学者基于InSAR技术获取了此次地震的同震形变场(温少妍等, 2020), 并反演了断层面滑动分布(李成龙等, 2020)。 本文在此基础上优化了断层面几何构形及滑动分布反演方法, 获取了此次地震更精细的断层面滑动分布, 然后对区域GPS面应变率、 构造动力背景和区域深浅部构造进行了综合分析, 探讨了2020年新疆伽师6.4级地震与1997— 2003年伽师6~7级震群的相互关系。

图 1 伽师地震区域构造
a 黑色细线表示断层(F1柯坪断裂; F2奥兹格尔他乌断裂; F3托克散阿塔能拜勒褶皱-逆断裂带; F4科克布克三山褶皱-逆断裂带; F5奥依布拉克褶皱-逆断裂带; F6喀拉铁克断裂; F7迈丹断裂; F8皮羌断裂; F9巴什托普断裂), 灰色粗线表示二级块体划分, 白色方框为升、 降轨SAR数据覆盖范围, 震源机制解数据来自GCMT(1976— 2020年)①(①https://www.globalcmt.org/。)和USGS(2020年)②(②https://earthquake.usgs.gov/earthquakes/。), 蓝色箭头表示GPS水平速度场; b 黑色实线表示二级块体划分, 红色圆圈表示1900年以来发生的7级以上地震, 红色方框表示图a的位置Fig. 1 Tectonic setting map of the Jiashi earthquake.

表1 不同机构提供的震源机制解 Table 1 Focal mechanism solutions of Jiashi MS6.4 earthquake
1 区域构造背景

2020年新疆伽师6.4级地震的震中位于南天山前陆褶皱冲断带, 其北部为天山褶皱带, 西南为西昆仑-帕米尔高原, 东部和南部为刚性的塔里木盆地(乔学军等, 2007)(图 1)。 受印度板块与欧亚板块持续碰撞的影响, 该区域的构造活动主要表现为天山地块逆冲于塔里木地块之上(田勤俭等, 2006), GPS观测到该区域的地壳缩短率为6~8mm/a(王琪等, 2001)。 在这样的逆冲环境下, 伽师地区中强地震频发, 仅过去的50a内, 距此次震中50km的范围内就发生了13次具有多种震源机制类型的6级以上地震(周仕勇等, 2006; 张竹琪等, 2008)。

从区域相对较小空间尺度的构造背景分析, 2020年新疆伽师6.4级地震的震源区位于南天山山前的柯坪推覆体。 新生代期间, 特别是第四纪以来, 天山地区强烈而持久的再生造山作用不仅形成了分布式逆断裂作用的格局, 使地壳近SN向缩短和垂向隆升, 还在天山南北麓形成逆断裂-褶皱带或背驮式逆断裂推覆体系, 引起天山山脉不断地向外侧扩展, 逆断裂-褶皱带的主体活动期也从山体一侧的中新世早期向前陆盆地逐渐变新(邓起东等, 2000)。 柯坪推覆体的形成及扩展也与滑脱面的伸展有关, 天山主体沿高角度逆断裂向南、 北两侧前陆盆地推挤的过程中, 在盆地下顺寒武系底界面形成了低角度滑脱断裂, 该滑脱断裂的形成和不断向盆地方向的伸展, 导致之上多排褶皱-逆断裂带的逐步形成与推覆体的扩展形成了柯坪推覆体(宋方敏等, 2006)。 推覆体由一系列近NEE-SWW走向、 向S逆冲的薄皮式冲断带组成(Allen et al., 1994), EW长约300km, SN宽60~140km, 北以迈丹断裂与南天山晚古生代造山带为界, 南以柯坪断裂与塔里木盆地为界, 在新生代天山造山带变形向盆地扩展的过程中扮演着极其重要的角色(张耀, 2019)。 跨柯坪断块的人工地震剖线揭示位于基底及其盖层间存在一个总体向天山倾斜的滑脱层, 具有上陡下缓的 “ 犁式” 特点(杨晓平等, 2008), 塔里木盆地的基底沿滑脱面自北向南逆冲推覆, 推覆体南浅北深, 滑脱面上部发育高角度断坡断层, 这些断层或出露地表或隐伏其下, 与滑脱层形成多排近EW向的叠瓦状铲式断层(李杰等, 2015)。 受地壳横向缩短不均匀性的影响, 柯坪推覆体被NNW向的皮羌 “ 类转换断层” (F8)分割为东、 西2段, 此次地震就发生在柯坪推覆体的西段, 该段由多排走向NEE— 近EW向的倒转背斜褶皱、 褶皱前缘断裂及其间背驮式盆地或谷地等组成, 且具有明显的构造分段性特征, 从南到北依次由柯坪塔格逆断裂-背斜(F1)、 奥兹格尔他乌逆断裂-背斜(F2)、 托克散阿塔能拜勒逆断裂-背斜(F3)、 科克布克三山逆断裂-背斜(F4)和奥依布拉克逆断裂-背斜(F5)等组成(徐锡伟等, 2006), 每个逆断裂-背斜以单面山的形式出现, 自新生代晚期向塔里木盆地逆冲推覆, 地壳缩短量约为45km, 地壳缩短率为37%(杨晓平等, 2006), 其中早期形成的褶皱-逆断裂带在北侧, 晚期形成的褶皱-逆断裂带在南侧, 反映推覆体由北向南的扩展(宋方敏等, 2006)。 地貌测量结果显示, 柯坪推覆体西段的逆断裂-褶皱带活动性从北向南逐渐增强, 表现出南天山多排逆断层在前展模式中的活动强度由后排构造向前排构造转移的特点。 从区域地层厚度分析, 认为柯坪逆断裂-褶皱带西段滑脱层的埋深约为9km, 东段滑脱层的埋深约为5km(杨晓平等, 2008), 向N缓倾的滑脱面在天山根部会聚到迈丹断裂上(吴传勇等, 2014)。 2020年伽师6.4级地震的发生为深入理解柯坪推覆体的动力学过程带来了新的启示。

2 2020年新疆伽师6.4级地震的InSAR同震形变场

此次新疆伽师6.4级地震的震中区域人口稀少, 传统大地测量观测站点分布稀疏, 难以获得高分辨率的同震形变场(图 1)。 InSAR技术作为常用的影像大地测量方法, 已被广泛应用于同震形变场的获取及发震断层几何参数及其滑动分布的反演(Massonnet et al., 1998; 乔学军等, 2007; 单新建等, 2009; 张国宏等, 2010; 宋小刚等, 2015; 季灵运等, 2017; 屈春燕等, 2017)。 基于此, 本文利用InSAR技术获取了2020年伽师6.4级地震的同震形变场, 并以此为约束反演其发震断层的几何参数及断层面的滑动分布特征。

本文利用升、 降轨Sentinel-1 SAR卫星数据获取了伽师6.4级地震的同震形变场(表2)。 InSAR数据处理在GAMMA商业软件平台上进行, 首先基于美国NASA发布的30m空间分辨率的SRTM DEM地形数据模拟并消除地形相位, 得到研究区的InSAR同震干涉图(图 2), 然后基于Delaunay三角网的最小费用网络流算法对差分干涉图进行相位解缠, 解缠起始点选取在远离震中的稳定点。 使用远场数据进行二次多项式拟合去除形变场中的残余轨道误差(Fialko et al., 2006), 基于地形相关法估计并去除大气相位延迟误差, 最后根据SAR卫星参数将相位转化为位移, 得到此次地震的升、 降轨LOS向同震形变场(图 3a, d)。

表2 Sentinel-1卫星影像参数 Table 2 Sentinel-1 image parameters

图 2 InSAR同震差分干涉图升(a)、 降轨(b)Fig. 2 Coseismic ascending(a)and descending(b) interferograms from InSAR.

图 3 均匀滑动模型参数直方图Fig. 3 Frequency histograms of uniform slip model parameters.

图 2 中的干涉条纹变化范围显示此次地震形变场EW向跨度约40km, SN向跨度约20km, 图中1个色周对应2.8cm的位移量, 条纹从正变化为负代表地表沿LOS向(雷达视线方向)靠近卫星运动, 干涉图的整体相干性较高且震源区不存在由于形变梯度过大导致的失相干, 表明同震破裂未达到地表。 形变区域内有柯坪和奥兹格尔他乌2条逆冲活动断裂, 干涉条纹的变化显示升、 降轨形变场在柯坪断裂以北(断层上盘)为显著上升, 在奥兹格尔他乌断裂以有北少量下沉, 2条轨道的形变趋势基本一致, 表明此次地震以垂向运动为主导, 符合逆冲型破裂引起的形变特征。 最终获取的同震形变(图 3a, d)结果显示, 升轨结果LOS向最大抬升量为5.9cm, 最大沉降量为3.7cm; 降轨结果LOS向最大抬升量约为6.4cm, 最大沉降量约为2cm, 与前人的结果较为一致(李成龙等, 2020; 温少妍等, 2020)。 结合InSAR形变场分布及已有断层的空间分布, 认为此次地震为典型的逆冲型地震破裂事件, 其发震断层为柯坪断裂。

3 2020年新疆伽师6.4级地震的同震滑动分布
3.1 均匀滑动模型

由于缺少此次地震发震断层的先验信息, 在反演断层面精细化滑动分布前, 需采用均匀滑动模型对发震断层的位置、 震源深度、 走向、 倾角等几何参数进行反演。 本文采用一种基于贝叶斯方法的均匀滑动反演方法, 该反演方法假设地震由单一断层面滑动引起, 发震断层面是一个嵌入在介质均匀且各向同性的弹性半空间中的矩形面(Okada, 1985), 其滑移量和滑动角都均匀分布。 反演时给定较大范围的上、 下界并随机给出相应的初始值, 采用马尔科夫链蒙特卡洛方法对参数空间进行采样, 基于贝叶斯方法构建后验概率函数, 通过迭代不断对参数空间进行优化直至收敛, 进而计算出最佳拟合参数(Bagnardi et al., 2018)。 图 3 所示的蓝色矩形为参数出现的频数, 红色竖线为最优解。

为提高计算效率, 在反演前利用四叉树分割方法对原始的高分辨率InSAR观测值进行降采样, 将采样的形变梯度阈值设为2.8cm, 最终得到1 327个升轨观测值和1 680个降轨观测值, 既保留了形变细节特征又能减少数据冗余。 反演得到的发震断层参数如表3所示, 说明此次地震是以倾滑运动为主的逆冲型地震。 图 3c、 3f为观测值与模拟值间的差值, 可以看出均匀滑动模型反演结果的残差较小, 表明了反演结果的可靠性(图 4)。

表3 均匀滑动反演的断层参数 Table 3 Fault parameters derived by uniform slip inversion

图 4 均匀滑动模型的反演结果
a、 d 升轨和降轨的InSAR同震形变场; b、 e 相应的模拟同震形变场; c、 f 残差图Fig. 4 Simulation results of the uniform-slip model.

3.2 分布式滑动模型

为了获取更为精细的断层面滑动分布, 本文以均匀滑动模型反演得到的发震断层几何参数为基础, 引入基于约束最小二乘的分布式滑动反演方法, 计算伽师地震发震断层面的精细滑动分布。

考虑到震源区的构造背景, 本文在前人研究的基础上(李成龙等, 2020), 对分布式滑动反演的断层面几何构型进行了改进。 首先将均匀滑动模型反演得到的断层长度扩展至40km。 跨柯坪断块的人工地震剖面揭示区域基底及其盖层间存在一个总体向天山倾斜的滑脱层, 具有上陡下缓的 “ 犁式” 特点(杨晓平等, 2008), 滑脱面上部发育高角度断坡断层, 这些断层或出露地表或隐伏其下, 与滑脱层形成多排近EW向的叠瓦状铲式断层(李杰等, 2015)。 因此, 根据均匀滑动反演的发震断层参数及柯坪断裂在地表出露位置间的几何关系, 将断层面设置为铲状, 其倾角从地表到地表之下4.25km深处设为57.8° , 4.25~8km深处设为10.7° , 这样的断层几何参数设置与深部地球物理结果一致。

最后, 将断层面剖分为930个边长为2km的等边三角形单元, 基于约束最小二乘方法计算每个三角形单元的滑动量, 通过最小化模型和观测值之间的残差确定最佳参数。 通过权衡数据拟合程度与断层滑动粗糙度, 对相邻断层片的滑动量施加平滑约束。

与均匀滑动模型结果相比, 分布式滑动模型结果的残差明显减小(图 5)。 图 6 为获取的断层滑动分布图, 结果显示断层面的最大倾滑量约为0.29m, 位于地下约5km深处, 最大走滑量约为0.03m。 计算得到同震的累计地震矩为1.73× 1018N· m, 合矩震级MW6.1。

图 5 分布式滑动模型的拟合结果
a、 d 升轨和降轨的InSAR同震形变场; b、 e 相应的模拟同震形变场; c、 f 残差图Fig. 5 Simulation results of the distributed-slip model.

图 6 伽师地震的同震滑动分布
a 走滑分量; b 倾滑分量Fig. 6 Slip distributions for the Jiashi earthquake.

4 讨论
4.1 基于现今GPS的区域应变率场分析

南天山前陆褶皱冲断带是大陆内部典型的逆冲断裂带, 地质结果和GPS观测数据表明天山南、 北两侧以逆断层-褶皱为代表的薄皮推覆构造自晚第四纪以来活动非常强烈, 这些构造带以水平向的地壳缩短和垂向的构造抬升吸收调节了天山地区至少一半的缩短变形(郑文俊等, 2019)。 此次伽师地震所处的柯坪推覆体作为塔里木盆地周缘惟一一个以寒武系为滑脱层的薄皮冲断带, 在新生代天山造山带变形向盆地内部扩展的过程中扮演着极其重要的角色(张耀, 2019)。 因此, 对柯坪推覆体及周边区域的形变特征进行研究有助于分析此次地震的发震背景。

本文基于GPS水平速度场(Wang et al., 2020)和Shen等(2015)提出的应变计算方法, 处理并获取了柯坪地块及其周边地区的面应变率场(图 7)。 研究区域的面应变率整体为负值, 说明该区域以挤压逆冲变形为主, 与区域所处的逆冲构造背景一致, 区域内发生的中强地震也以逆冲型居多。 值得注意的是, 在塔里木盆地与天山交界带西侧存在一处面应变率显著降低的区域(图 7 中的蓝色区域), 说明此处受到挤压更为强烈, 此次伽师地震的震中位于该区域中应变率梯度变化较大处。 此外, 在该低值区东侧, 为了调节由于柯坪推覆体前展变形速率差异较大引起的横向挤压缩短不均匀性, 派生出NNW向的皮羌 “ 类转换断层” , 该断层由柯坪塔格推覆构造系西段各推覆体前缘逆断裂东端的NNW向弧形拐弯段串联而成, 是一条不连续的、 代表两侧不等量地壳缩短的边界(徐锡伟等, 2006), 该断裂历史上也发生过6.8级地震, 表明其现今可能处于活动状态。

图 7 伽师地震及其邻近地区的面应变率分布图Fig. 7 The plane strain rate field around the Jiashi earthquake and its surrounding areas.

4.2 2020年新疆伽师6.4级地震的发震断层与构造动力背景分析

柯坪推覆体由多条近EW向的逆断裂-背斜带组成, 沿着统一的滑脱层由北向南逆冲推覆(张耀, 2019), 其滑脱面发育在地下5~9km深处, 向N缓倾的滑脱面在天山根部会聚到迈丹断裂上(吴传勇等, 2014), 推覆体南浅北深, 总体向N倾(曲国胜等, 2003)。 根据前文分析, 此次地震的发震断层为铲型结构, 从地表至地下深部由高倾角逐渐过渡为低倾角, 震源深度约为5km, 印证了柯坪地块为薄皮推覆运动的地学认识(陈杰等, 2001)。

图 8 指示了此次地震的发震构造动力背景, 沿过震中的剖线AA'(位置见图 1)绘制出柯坪推覆体的构造剖面图(田勤俭等, 2006; 李杰等, 2015), 将剖面线两侧50km范围内的GPS水平向速率(Wang et al., 2020)投影到垂直于柯坪断裂的方向, 同时提取了剖面线两侧10km范围内的小地震精定位数据及地形数据(图 8)。 可以看出, 垂直于柯坪断裂的GPS运动速率自北向南呈递减趋势, 速率的最大差值约为5mm/a, 说明柯坪推覆体受SN向挤压变形较大。 因此, 此次地震可能是由柯坪地块向塔里木盆地推覆的过程中前方冲断席发生错断所引起的, 且发震断层位于柯坪推覆体前缘小地震较稀疏的区域。 值得注意的是, 西南天山的构造变形不仅分配在推覆体前缘的新生代逆断裂-褶皱带上, 通过山根处的迈丹断裂(F7)全新世活动证据发现推覆体根部断裂现在仍强烈活动, 吸收了一部分缩短变形(吴传勇等, 2014), 从图 7 可以看出其小地震活动也较为稀疏, 该区的地震危险性不容忽视。


图 8 柯坪推覆体的深、 浅部构造剖面图
a 地形与垂直于柯坪断裂走向的GPS速率剖面; b 小地震剖面与深部构造解译。AA'剖面的位置见图 1, 黑色曲线代表F1、 F2、 F6、 F7 4条断层和深部滑脱面, 小地震数据时间范围为1980— 2015年Fig. 8 The shallow and deep tectonic profile of the Kalpin nappe.

4.3 2020年新疆伽师6.4级地震与1997— 2003年伽师6~7级震群相互关系的讨论

1997— 2003年, 南天山前陆环境不同断裂带相继破裂引发了3组6~7级强烈地震(图 1), 震群主要位于塔里木盆地边缘的麦盖提斜坡带与喀什凹陷的过渡边界。 由于该区域位于塔里木盆地西北端、 帕米尔弧型构造东北部以及天山褶皱带柯坪块体的交会部位, 地质构造环境和受力因素复杂(高国英等, 2004), 因此地震的震源机制复杂多变, 既有逆冲型, 也有走滑和正断型(赵翠萍等, 2008)(图 1)。 其中, 1997年1月— 2003年1月伽师震源区地震的震源机制解以走滑型为主, 且有一定数量的地震为正断型, 与区域所处的逆冲环境不一致。 而2003年巴楚-伽师地震的震源机制解则以逆断型为主, 与区域震源机制解类型较为一致(徐锡伟等, 2006)。

地震是新构造运动的表现形式之一, 地震发生的地点、 强度与地壳最新的构造变动样式有着密切的联系。 已有的地震精确定位、 震源机制解、 地壳— 上地幔三维S波速度结构和各种探测成果显示, 1997年伽师强震群的主体发震构造为NNW向下苏洪隐伏断层, 该断层为右旋走滑的NNW向 “ 类转换断层” , 起到分割、 调节西大东小的地壳缩短边界转换作用(刘启元等, 2000; 徐锡伟等, 2006)。 而阿图什1996年6.9级和1998年6.0级地震、 2003年巴楚-伽师地震、 柯坪块体发生的多次6级地震和此次伽师6.4级地震, 无论是震源断错性质还是主压应力方向或震源破裂特征都很一致, 反映出挤压环境条件下的地壳缩短作用, 表现出柯坪块体长期以来的断层运动特性。

图 9 柯坪推覆体及邻近区域的构造背景示意图Fig. 9 Schematic map showing regional tectonic background of the Kalpin nappe and surrounding regions.

从1997— 2003年伽师6~7级地震震群与此次伽师6.4级地震所处的构造环境来看, 均发生在南天山逆冲推覆构造体系的前陆盆地环境中。 对柯坪推覆体的构造动力环境进行简化, 如图 8 所示。 可以看出, 在整体逆冲推覆的构造背景中, 存在主逆断层(柯坪断裂)以及在前陆盆地下方发育的隐伏逆冲断层(F9), 隐伏断裂带的形成可能是由塔里木盆地北缘的非均匀变形及现今构造运动所引起(刘启元等, 2000), 随着逆冲推覆带向盆地发展, 一些断层向前伸展到前陆盆地下方, 被较厚的新生代沉积层覆盖, 形成隐伏断层, 这些隐伏断层也具有孕育中强地震的能力(张竹琪等, 2008)。 在逆冲推覆体上, 由于岩石力学属性存在差异且区域构造应力场在不同部位也存在一定的差异性, 故在大的逆冲推覆背景下可能会伴生发育规模相对较小的正断层, 甚至还伴生发育与主逆断层走向近垂直的横向走滑断层(皮羌断裂、 下苏洪隐伏断层), 这些伴生正断层和走滑断层也都具备发生中强地震的能力。 因此, 根据区域构造动力背景及区域构造分布进行综合分析, 认为1997— 2003年伽师6~7级地震震群与此次伽师6.4级地震的发生机理应为南天山前陆环境中不同级别、 不同性质的断裂相继破裂形成不同类型的破裂事件, 应为同源事件, 共同受控于南天山的逆冲推覆作用。

5 结论

(1)本文利用升、 降轨Sentinel-1 SAR卫星数据获取了2020年新疆伽师6.4级地震的InSAR同震形变场, 并以InSAR形变场为约束, 反演获取了断层面的精细滑动分布。 结果显示, 此次地震的发震断层为柯坪断裂。 断层面最大倾滑量约为0.29m, 位于地下约5km深处, 通过同震释放能量计算得到的震级为MW6.1, 是一次以逆冲为主兼具少量走滑运动的地震破裂事件。

(2)2020年新疆伽师6.4级地震是发生在南天山前陆环境中的一次典型逆冲破裂事件, 与1997— 2003年在该区域发生的多次6~7级地震一起构成系列强震群。 从区域构造动力背景和区域深浅部构造综合分析, 认为1997— 2003年伽师6~7级地震震群与此次伽师6.4级地震的发生机理应为南天山前陆环境中不同级别、 不同性质的断裂相继破裂形成不同类型的破裂事件, 共同受控于南天山的逆冲推覆作用。

致谢 中国地震局地球物理研究所房立华研究员为本研究提供了小地震数据; 欧洲航天局提供了Sentinel-1卫星影像。 在此一并表示感谢!

参考文献
[1] 陈杰, 丁国瑜, Burbank D W, 等. 2001. 中国西南天山山前的晚新生代构造与地震活动[J]. 中国地震, 17(2): 134—155.
CHEN Jie, DING Guo-yu, Burbank D W, et al. 2001. Late Cenozoic tectonics and seismicity in the southwestern Tianshan, China[J]. Earthquake Research in China, 17(2): 34—55(in Chinese). [本文引用:1]
[2] 邓起东, 冯先岳, 张培震, 等. 2000. 天山活动构造 [M]. 北京: 地震出版社.
DENG Qi-dong, FENG Xian-yue, ZHANG Pei-zhen, et al. 2000. Active Tectonics of the Chinese Tianshan Mountains [M]. Seismological Press, Beijing(in Chinese). [本文引用:1]
[3] 高国英, 聂晓红, 夏爱国. 2004. 2003年伽师6. 8级地震序列特征和震源机制的初步研究[J]. 中国地震, 20(2): 179—186.
GAO Guo-ying, NIE Xiao-hong, XIA Ai-guo. 2004. Preliminary study on characteristics of the Jiashi earthquake sequence with MS6. 8 and its focal mechanism solution in 2003[J]. Earthquake Research in China, 20(2): 179—186(in Chinese). [本文引用:1]
[4] 季灵运, 刘传金, 徐晶, 等. 2017. 九寨沟 MS7. 0地震的InSAR观测及发震构造分析[J]. 地球物理学报, 60(10): 4069—4082.
JI Ling-yun, LIU Chuan-jin, XU Jing, #magtechI#, et al. 2017. InSAR observation and inversion of the seismogenic fault for the 2017 Jiuzhaigou MS7. 0 earthquake in China[J]. Chinese Journal of Geophysics, 60(10): 4069—4082(in Chinese). [本文引用:1]
[5] 李成龙, 张国宏, 单新建, 等. 2020. 2020年1月19日新疆伽师县 MS6. 4地震InSAR同震形变场与断层滑动分布反演[J/OL]. 地球物理学进展. kns. cnki. netkcmsdetail/11. 2982. P. 20200608. 1616. 178. html.
LI Cheng-long, ZHANG Guo-hong, SHAN Xin-jian, #magtechI#, et al. 2020. Coseismic deformation and slip distribution of the MS6. 4 Jiashi, Xinjiang earthquake revealed by Sentinel-1A SAR imagery [J/OL]. Progress in Geophysics. kns.cnki.netkcmsdetail/11.2982.P.20200608.1616.178.html. [本文引用:3]
[6] 李杰, 乔学军, 杨少敏, 等. 2015. 西南天山地表三维位移场及断层位错模型[J]. 地球物理学报, 58(10): 3517—3529.
LI Jie, QIAO Xue-jun, YANG Shao-min, #magtechI#, et al. 2015. Detachment fault model characterized for the 3D surface displacement field in the southwestern Tian Shan[J]. Chinese Journal of Geophysics, 58(10): 3517—3529(in Chinese). [本文引用:3]
[7] 刘启元, 陈九辉, 李顺成, 等. 2000. 新疆伽师强震群区三维地壳上地幔S波速度结构及其地震成因的探讨[J]. 地球物理学报, 43(3): 356—365.
LIU Qi-yuan, CHEN Jiu-hui, LI Shun-cheng, #magtechI#, et al. 2000. Passive seismic experiment in the Xingjiang-Jiashi strong earthquake region and discussions on its seismic genesis[J]. Chinese Journal of Geophysics, 43(3): 356—365(in Chinese). [本文引用:2]
[8] 孟令媛, 苑争一, 宋治平, 等. 2020. 2020年1月19日新疆伽师 MS6. 4地震总结[J]. 地震地磁观测与研究, 41(2): 63—89.
MENG Ling-yuan, YUAN Zheng-yi, SONG Zhi-ping, et al. 2020. About Jiashi MS6. 4 earthquake in Xinjiang on Jan. 19, 2020[J]. Seismological and Geomagnetic Observation and Research, 41(2): 63—89(in Chinese). [本文引用:1]
[9] 聂晓红, 刘建明, 向元, 等. 2020. 2020年1月19日新疆伽师 MS6. 4地震前地震活动特征分析[J]. 内陆地震, 34(1): 10—19.
NIE Xiao-hong, LIU Jian-ming, XIANG Yuan, #magtechI#, et al. 2020. Analysis of seismic activity before Jiashi MS6. 4 earthquake in Xinjiang on January 19th, 2020[J]. Inland Earthquake, 34(1): 10—19(in Chinese). [本文引用:1]
[10] 乔学军, 郭利民. 2007. 新疆伽师强震群区的InSAR观测研究[J]. 大地测量与地球动力学, 27(1): 7—13.
QIAO Xue-jun, GUO Li-min. 2007. Study on Jiashi strong earthquake swarm area by InSAR[J]. Journal of Geodesy and Geodynamics, 27(1): 7—13(in Chinese). [本文引用:2]
[11] 屈春燕, 左荣虎, 单新建, 等. 2017. 尼泊尔 MW7. 8地震InSAR同震形变场及断层滑动分布[J]. 地球物理学报, 60(1): 151—162.
QU Chun-yan, ZUO Rong-hu, SHAN Xin-jian, et al. 2017. Coseismic deformation field of the Nepal MS8. 1 earthquake from Sentinel-1A/InSAR data and fault slip inversion[J]. Chinese Journal of Geophysics, 60(1): 151—162(in Chinese). [本文引用:1]
[12] 曲国胜, 李亦纲, 陈杰, 等. 2003. 柯坪塔格推覆构造几何学、 运动学及其构造演化[J]. 地学前缘, 10(S1): 142—152.
QU Guo-sheng, LI Yi-gang, CHEN Jie, #magtechI#, et al. 2003. Geometry, kinematics and tectonic evolution of Kepingtage thrust system[J]. Earth Science Frontiers, 10(S1): 142—152(in Chinese). [本文引用:1]
[13] 单新建, 屈春燕, 宋小刚, 等. 2009. 汶川 MS8. 0地震InSAR同震形变场观测与研究[J]. 地球物理学报, 52(2): 496—504.
SHAN Xin-jian, QU Chun-yan, SONG Xiao-gang, et al. 2009. Coseismic surface deformation caused by the Wenchuan MS8. 0 earthquake from InSAR data analysis[J]. Chinese Journal of Geophysics, 52(2): 496—504(in Chinese). [本文引用:1]
[14] 宋方敏, 闵伟, 韩竹军, 等. 2006. 柯坪塔格推覆体的新生代变形与扩展[J]. 地震地质, 28(2): 224—233.
SONG Fang-min, MIN Wei, HAN Zhu-jun, et al. 2006. Cenozoic deformation and propagation of the Kalpintag fold nappe[J]. Seismology and Geology, 28(2): 224—233(in Chinese). [本文引用:2]
[15] 宋小刚, 申星, 姜宇, 等. 2015. 通过InSAR与GPS数据融合获取汶川地震同震三维形变场[J]. 地震地质, 37(1): 222—231. doi: 103969/j.issn.0253-4967.2015.01.017.
SONG Xiao-gang, SHEN Xing, JIANG Yu, #magtechI#, et al. 2015. Coseismic 3D deformation field acquisition of the Wenchuan earthquake based on InSAR and GPS data[J]. Seismology and Geology, 37(1): 222—231(in Chinese). [本文引用:1]
[16] 田勤俭, 丁国瑜, 郝平. 2006. 南天山及塔里木北缘构造带西段地震构造研究[J]. 地震地质, 28(2): 213—223.
TIAN Qin-jian, DING Guo-yu, HAO Ping. 2006. Seismotectonic study on the west part of the interaction zone between southern Tianshan and northern Tram[J]. Seismology and Geology, 28(2): 213—223(in Chinese). [本文引用:2]
[17] 王琪, 张培震, 牛之俊, 等. 2001. 中国大陆现今地壳运动和构造变形[J]. 中国科学(D辑), 31(7): 529—536.
WANG Qi, ZHANG Pei-zhen, NIU Zhi-jun, #magtechI#, et al. 2001. Present-day crustal movement and tectonic deformation in China continent[J]. Science in China(Ser D), 31(7): 529—536(in Chinese). [本文引用:1]
[18] 温少妍, 李成龙, 李金. 2020. 2020年1月19日新疆伽师 MS6. 4地震InSAR同震形变场特征及发震构造初步探讨[J]. 内陆地震, 34(1): 1—9.
WEN Shao-yan, LI Cheng-long, LI Jin. 2020. Preliminary discussion on characteristics of coseismic deformation field and InSAR seismogenic structure for Xinjiang Jiashi MS6. 4 earthquake on January 19th, 2020[J]. Inland Earthquake, 34(1): 1—9(in Chinese). [本文引用:2]
[19] 吴传勇, 阿里木江, 戴训也, 等. 2014. 西南天山迈丹断裂东段晚第四纪活动的发现及构造意义[J]. 地震地质, 36(4): 976—990.
WU Chuan-yong, Alimujiang, DAI Xun-ye, #magtechI#, et al. 2014. Discovery of the late-Quaternary activity along the eastern segment of Maidan Fault in southwest Tianshan and its tectonic implication[J]. Seismology and Geology, 36(4): 976—990(in Chinese). [本文引用:3]
[20] 徐锡伟, 张先康, 冉勇康, 等. 2006. 南天山地区巴楚-伽师地震( MS6. 8)发震构造初步研究[J]. 地震地质, 28(2): 161—178.
XU Xi-wei, ZHANG Xian-kang, RAN Yong-kang, #magtechI#, et al. 2006. The preliminary study on seismotectonics of the 2003 AD Bachu-jiashi earthquake( MS6. 8), southern Tian Shan[J]. Seismology and Geology, 28(2): 161—178(in Chinese). [本文引用:4]
[21] 杨晓平, 冉勇康, 宋方敏, 等. 2006. 西南天山柯坪逆冲推覆构造带的地壳缩短分析[J]. 地震地质, 28(2): 194—204.
YANG Xiao-ping, RAN Yong-kang, SONG Fang-min, #magtechI#, et al. 2006. The analysis for crust shortening of Kalpin thrust tectonic zone, south-western Tianshan, Xinjiang, China[J]. Seismology and Geology, 28(2): 194—204(in Chinese). [本文引用:1]
[22] 杨晓平, 邓起东, 张培震, 等. 2008. 天山山前主要推覆构造区的地壳缩短[J]. 地震地质, 30(1): 111—131.
YANG Xiao-ping, DENG Qi-dong, ZHANG Pei-zhen, #magtechI#, et al. 2008. Crustal shortening of major nappe structures on the front margins of the Tianshan[J]. Seismology and Geology, 30(1): 111—131(in Chinese). [本文引用:3]
[23] 张国宏, 屈春燕, 宋小刚, 等. 2010. 基于InSAR同震形变场反演汶川 MW7. 9地震断层滑动分布[J]. 地球物理学报, 53(2): 269—279.
ZHANG Guo-hong, QU Chun-yan, SONG Xiao-gang, et al. 2010. Slip distribution and source parameters inverted from co-seismic deformation derived by InSAR technology of Wenchuan MW7. 9 earthquake[J]. Chinese Journal of Geophysics, 53(2): 269—279(in Chinese). [本文引用:1]
[24] 张耀. 2019. 南天山山前柯坪冲断带变形的几何学、 运动学与机制研究[D]. 杭州: 浙江大学.
ZHANG Yao. 2019. Geometry, kinematics and deformation mechanism of the Kepine fold-thrust belt in the South Tian Shan, NW China [D]. Zhejiang University, Hangzhou(in Chinese). [本文引用:3]
[25] 张竹琪, 陈永顺, 林间. 2008. 1997年伽师震群中相邻正断层和走滑断层之间相互应力作用[J]. 中国科学(D辑), 38(3): 334—342.
ZHANG Zhu-qi, CHEN Yong-shun, LIN Jian. 2008. Interaction between adjacent normal faults and strike slip faults in Jiashi earthquake swarm in 1997[J]. Science in China(Ser D), 38(3): 334—342(in Chinese). [本文引用:2]
[26] 赵翠萍, 陈章立, 郑斯华. 2008. 1998—2003年伽师3次不同类型 MS6地震震源破裂过程及短期内余震活动特征[J]. 地球物理学报, 51(4): 1093—1102.
ZHAO Cui-ping, CHEN Zhang-li, ZHENG Si-hua. 2008. Source rupture process of 3 Jiashi MS6 events(1998—2003) and its correlation with the aftershock activity[J]. Chinese Journal of Geophysics, 51(4): 1093—1102(in Chinese). [本文引用:1]
[27] 郑文俊, 张培震, 袁道阳, 等. 2019. 中国大陆活动构造基本特征及其对区域动力过程的控制[J]. 地质力学学报, 25(5): 699—721.
ZHENG Wen-jun, ZHANG Pei-zhen, YUAN Dao-yang, #magtechI#, et al. 2019. Basic characteristics of active tectonics and associated geodynamic processes in continental China[J]. Journal of Geomechanics, 25(5): 699—721(in Chinese). [本文引用:1]
[28] 周仕勇, 姜明明, Robinson R. 2006. 1997年新疆伽师强震群发展过程中发震断层间相互作用的影响[J]. 地球物理学报, 49(4): 1102—1109.
ZHOU Shi-yong, JIANG Ming-ming, Robinson R. 2006. The role of fault interactions in the generation of the 1997 Jiashi strong earthquake swarm, Xinjiang, China[J]. Chinese Journal of Geophysics, 49(4): 1102—1109(in Chinese). [本文引用:1]
[29] Allen M B, Windley B F, Zhang C. 1994. Cenozoic tectonics in the Urumgi-Korla region of the Chinese Tien Shan[J]. International Journal of Earth Sciences, 83(2): 406—416. [本文引用:1]
[30] Bagnardi M, Hooper A. 2018. Inversion of surface deformation data for rapid estimates of source parameters and uncertainties: A Bayesian approach[J]. Geochemistry Geophysics Geosystems, 19(7): 2194—2211. [本文引用:1]
[31] Fialko Y. 2006. Interseismic strain accumulation and the earthquake potential on the southern San Andreas fault system[J]. Nature, 441(7096): 968—971. [本文引用:1]
[32] Massonnet D, Feigl K L. 1998. Radar interferometry and its application to changes in the Earth’s surface[J]. Reviews of Geophysics, 36(4): 441—500. [本文引用:1]
[33] Okada Y. 1985. Surface deformation due to shear and tensile faults in a half-space[J]. Bulletin of the Seismological Society of America, 75(4): 1135—1154. [本文引用:1]
[34] Shen Z, Wang M, Zeng Y, et al. 2015. Optimal interpolation of spatially discretized geodetic data[J]. Bulletin of the Seismological Society of America, 105(4): 2117—2127. [本文引用:1]
[35] Wang M, Shen Z. 2020. Present-day crustal deformation of continental China derived from GPS and its tectonic implications[J]. Journal of Geophysical Research, 125(2): 1—22. [本文引用:2]