2001年昆仑山地震同震-震后效应与2021年玛多地震关系探讨.pdf

第6 6 卷第7 期2023年7 月屈春燕，赵德政，单新建等.2 0 2 3.2 0 0 1年昆仑山地震同震一震后效应与2 0 2 1年玛多地震关系探讨地球物理学报，6 6（7）：2741-2756,doi:10.6038/cjg2022Q0273.Qu C Y,Zhao D Z,Shan X J,et al.2023.Coseismic and postseismic deformation of the 2001 Mw7.8 Kunlun Mountainearthquake and its loading effect on the 2021 Mw7.3 Madoi earthquake.Chinese J.Geophys.(in Chinese),66(7):2741-2756,doi:10.6038/cjg2022Q0273.地球物理学报CHINESE JOURNAL OF GEOPHYSICSVol.66,No.7Jul.,20232001年昆仑山地震同震一震后效应与2021年玛多地震关系探讨屈春燕，赵德政*，单新建，岳冲，陈晗，吴东霖，刘恋，张国宏中国地震局地质研究所地震动力学国家重点实验室，北京10 0 0 2 9摘要位于青藏高原中北部的巴颜喀拉地块是我国西部近年来的主体地震活动区，一系列Mw7.0以上强震均发生在该次级块体周边，而其北边界东昆仑断裂带是一条长达2 0 0 0 km、规模最大、活动性最强的深大断裂带.2 0 0 1年在东昆仑断裂带中段发生了Mw7.8昆仑山地震，2 0 2 1年5月在其震中东南部大约450 km处巴颜喀拉块体内部一次级断裂上发生了Mw7.3玛多地震.玛多地震对人们以往认为强震更可能发生在巴颜喀拉块体边界断裂上的认识提出挑战，但是也为研究巴颜喀拉块体边界断裂与块体内部次级断裂活动关系、地震触发关系带来机遇.本文利用前期基于2 0 0 1年昆仑山地震后积累的大量InSAR数据获得的震后大范围形变场时空演化图像和库仑应力变化模型，探讨昆仑山地震与玛多地震的关系.InSAR震后观测结果显示：昆仑山地震后沿东昆仑断裂带出现了长达500km的大范围南北不对称震后形变场，其中南盘形变宽度和量级均明显大于北盘，南盘形变宽度达到2 50 km，断层近场相对平均形变速率达到 2 0 mma-1,而且南盘向南衰减梯度小，整体衰减缓慢，意味着震后形变对巴颜喀拉块体形成持续东向加载作用，并将分摊到块体内部的一系列次级断裂上，应力加载增加次级断层的地震危险性.2 0 15一2 0 2 0 年InSAR震间应变率场则显示次级断裂昆仑山口一江错断裂呈高剪切应变率特征.本文计算了昆仑山地震同震破裂和震后形变引起的玛多震区多条SE向次级断裂的累积库仑应力变化，结果显示昆仑山地震同震和震后形变对玛多地震发震断裂（昆仑山口一江错断裂）形成了一定的应力加载.本文认为昆仑山地震同震和长时间尺度震后形变加速了巴颜喀拉块体的东向运动，而断层本身运动学性质和区域应力扰动共同影响了玛多地震的发生。关键词巴颜喀拉块体；东昆仑断裂；昆仑山地震；玛多地震；震后效应doi:10.6038/cjg2022Q0273中图分类号P541，P2 2 3收稿日期2 0 2 2-0 4-2 1，2 0 2 2-10-18 收修定稿Coseismic and postseismic deformation of the 2001 Mw7.8 Kunlun Mountain earthquakeand its loading effect on the 2021 Mw7.3 Madoi earthquakeQU ChunYan,ZHAO DeZheng,SHAN XinJian,YUE Chong,CHEN Han,WU DongLin,LIU Lian,ZHANG GuoHongState Key Laboratory of Earthquake Dynamics,Institute of Geology,China Earthquake Administration,Beijing100029,ChinaAbstract The Bayan Har block in the north-central Tibetan Plateau has been the main seismic基金项目地震动力学国家重点实验室课题（LED2019A02），中国地震局地质研究所基本科研业务费专项（IGCEA2104，IG CEA 18 0 9）和自然科学基金（418 7 2 2 2 9)联合资助。第一作者简介屈春燕，女，196 6 年生，研究员，主要从事大地测量研究E-mail：d q y q u c h y 16 3.c o m*通讯作者赵德政，男，1992 年生，博士，主要从事大地测量研究E-mail:2742region in west China in recent years,and a series of major earthquakes with Mw7.0 haveoccurred around this block.Its northern boundary,the East Kunlun fault,is an active large-scalefault zone with a total length of 2 0 0 0 k m.T h e 2 0 0 1 M w 7.8 K u n lu n M o u n t a i n(H o h Xi l)earthquake occurred in the middle segment of the East Kunlun fault(or Kunlun fault),and theMw7.3 Madoi earthquake occurred on a secondary fault inside the Bayan Har block,450 k msoutheast of the Kunlun Mountain earthquake.The Madoi earthquake challenges the previouslong-held view that strong earthquakes are more likely to occur on the boundary faults of theBayan Har block.This earthquake also provides an opportunity to study the relationship amongthe boundary faults and secondary faults,and earthquakes triggering within the block.This paperexplores the potential relationship between the Kunlun Mountain earthquake and the Madoiearthquake using postseismic InSAR observations and the Coulomb stress change,which areobtained from a large amount of SAR data acquired after the Kunlun Mountain earthquake.TheInSAR postseismic deformation shows a widespread asymmetric pattern along the rupture zone,and the deformation width and magnitude of the southern part are significantly larger than thenorthern one.The deformation width of the southern part reaches 25o km,and the near-fielddeformation rate reaches 20 mm a-1.The decay of displacement rate on the southern part isinsignificant,implying that the postseismic deformation forms a enduring eastward loading effecton the Bayan Har block,and will be partitioned into a series of secondary faults inside the block.The stress loading increases the seismic potential of the secondary faults.The interseismic strainmap during 20152020 shows a high strain rate zone along the secondary fault,the KunlunMountain Pass-Jiangcuo fault.We calculate the cumulative Coulomb stress change caused bythe coseismic rupture and postseismic processes of the Kunlun Mountain earthquake.The resultsshow that the stress variation on the seismogenic fault of the Madoi earthquake,the KunlunPassJiangcuo Fault,is significant.Based on these results,we suggest that the coseismic andpostseismic deformation of the Kunlun Mountain earthquake may enhance the eastward motion ofthe Bayan Har block and thus contribute to the occurrence of the Madoi earthquake.Keywords Bayan Har block;East Kunlun fault;Kunlun Mountain earthquake;Madoi earthquake;引言0青藏高原是中国大陆乃至全世界的主要活动构造区和强震多发区，而巴颜喀拉地块则是青藏高原目前的主体地震活动区.在印度板块和欧亚板块相互碰撞挤压、印度板块以40 50 mma-1的速率向北俯冲推挤的构造背景下，青藏高原内部发育有一系列规模巨大的活动断裂带，构成了高原内部重要次级活动地块的边界（图1;Molnar and Tapponnier，1978).这些边界断裂带以走滑运动的方式吸收了部分板块汇聚的能量，同时控制和主导着青藏高原系列强震的孕育和发生（Tapponnier etal.，198 2；Zhangetal.，2 0 0 4).位于青藏高原中北部、横贯整个中国西部的长条状巴颜喀拉地块，其南北边界分地球物理学报(Chinese J.Geophys.）Postseismic effects66卷别受到规模巨大、活动性强的甘孜一玉树一鲜水河断裂带和东昆仑断裂带的控制，是青藏高原内部构造变形强烈、历史和现代强震频发的重要次级活动块体，对青藏高原的现代地壳运动、断裂活动和强震孕育起主导和控制作用（邓起东等，2 0 14；张培震等，2003；郑文俊等，2 0 16).近30 年来，巴颜喀拉地块的构造运动和地震活动明显增强.自1997 年以来，在我国青藏高原发生的Ms7.0以上大地震，包括2 次Ms8.0以上特大地震，都发生在巴颜喀拉地块的边界断裂带上.其中包括1997 年Ms7.5玛尼地震、2 0 0 1年Mw7.8昆仑山地震、2 0 0 8 年Ms8.0汶川地震、2 0 10 年Ms7.1玉树地震、2 0 13年Ms7.0芦山地震和2 0 17 年Ms7.0九寨沟地震(图1a).这一系列强震的发生均与巴颜喀拉地块的活动密切相关（徐锡伟等，2 0 0 2；邓起东等，7期屈春燕等：2 0 0 1年昆仑山地震同震一震后效应与2 0 2 1年玛多地震关系探讨274381E40N(a)38N84E87E90E93E96E99E102E105E40N38N36N34N32N30N28N81E90E37N(b）36N35N34NP33N32N90E（a）青藏高原二级块体划分、活动断裂分布、区域构造背景略图，红色沙滩球：19 9 7 年以来发生在巴颜喀拉块体边界的7 级以上地震震中及震源机制解；紫色沙滩球：玛多地震震源机制解；白色粗线：二级块体边界；灰色细线：活动断裂；蓝色箭头为震间GPS速率（Wa n g a n d Sh e n，2 0 2 0）；（b）昆仑山地震和玛多地震区域的局部放大图，红色虚线：昆仑山口江错断裂带；红色粗线：东昆仑断裂带历史地震破裂段，红色圆点：玛多地震余震分布，红色和黑色震源球：玛多地震和历史地震震源机制解。Fig.1 Tectonic background of the Bayan Har block and its surrounding areas(a)Blocks,active faults and regional tectonics in the Tibetan Plateau.Red beach ball:focal mechanisms of earthquakes with M7that occurred on the boundary of the Bayan Har block since 1997;purple beach ball:focal mechanism of the Madoi earthquake;whitethick line:block boundary;gray thin line:active faults;blue arrows indicate the interseismic GPS velocity(Wang and Shen,2020);(b)Enlarged map of the Kunlun Mountain earthquake and the Madoi earthquake.Red dotted line:Kunlun Mountain pass-Jiangcuofault;red thick line:historical earthquake ruptures on the East Kunlun fault;red dot:aftershocks of the Madoi earthquake;red and2014；程佳和徐锡伟，2 0 18；单新建等，2 0 17），是该块体边界断裂带运动变形增强和应力应变快速积累的表现，意味着这一块体的构造运动进人新的活跃期，具备了强震孕育发生的构造条件，可能是当前及未来青藏高原的主体地震活动区.特别是其北边界东昆仑断裂带，规模最大，构造活动性更强，对巴颜喀拉块体的运动变形和地震活动性起调节和控制的主导作用（邓起东等，2 0 14）.构成巴颜喀拉块体北边界的东昆仑断裂带是青藏高原内部一条古老的缝合带，是调节青藏高原内36N颜略拉块体2024M10864292E2001,M.7.8昆仓山口江错断翠带2021M.7.3巴颜喀拉块体92E94E图1巴颜喀拉地块及周边区域地震构造背景略图black beach ball:focal mechanisms of the Madoi earthquake and historical earthquakes.34N32N30N28N84E87E94E90E96E东昆仑断裂带96E98E部挤压变形、地壳物质向东逃逸，并将青藏高原划分为南、北两部分的重要地质构造带（He and Chry，2008;Li et al.,2018;Loveless and Meade,2011).该断裂带西起青海、新疆交界处的鲸鱼湖以西，东至甘肃玛曲以东，全长达2 0 0 0 km，总体走向NWW，主要倾向北，倾角6 0 8 5为一高角度左旋走滑断裂带（徐锡伟等，2 0 0 2；Kirbyetal.，2 0 0 7；邓起东等，2 0 14).根据前人研究结果，东昆仑断裂带可以划分为7 个一级地震破裂段，从西到东分别是：鲸鱼湖段、库赛湖段、东大滩一西大滩段、阿拉克湖段、阿尼93E98E96E100E100E99E102E102E102E104E104E105E37N36N35N34N33N32N2744玛卿山段、玛沁段和山段，更宏观地可分为西段、中段和东段（徐锡伟等，2 0 0 2；VanDer Woerd etal.，2 0 0 2).西段一中段（大致以东经9 9 E一10 0 E为界)第四纪以来的长期平均滑动速率在10 1214-12-断层滑动速率变化8642-0848586 878889990919929394995969798991大地测量学速率：Diao et al.(2019)地质学速率：Van derWoerdet al.(2002)图2 地质学及GPS观测得到的东昆仑断裂带沿断层走向滑动速率（修改自Zhaoetal.，2 0 2 2）Fig.2 Fault slip rates along the strike of the East Kunlun fault,suggested by geological investigations and东昆仑断裂带是一条强震活动带，在其西段和中段曾发生过19 0 2 年Ms7.0秀沟地震，19 37 年Ms7.5花石峡地震，19 6 3年Ms7.0阿拉克湖地震及19 9 7 年Ms7.5玛尼地震等强震.而最近的一次即2 0 0 1年Mw7.8昆仑山地震，是有史以来发生在该断裂上的最大震级地震，也是青藏高原在1951年西藏当雄8.6 级地震以来的最大地震.昆仑山地震造成42 6 km长的地表破裂带和约7 8 m的最大左旋位移（徐锡伟等，2 0 0 2；陈杰等，2 0 0 4；杨国华等，2008).相比于东昆仑断裂带的中段和西段，其东段构造活动性较弱，地震活动较少，间接反映出东昆仑断裂带自西向东活动性减弱、断裂带分段活动性明显的特点（图2）2021年5月2 2 日发生在巴颜喀拉块体东部块体内部次级断层上的Mw7.3玛多地震，与早前及历史地震均发生在该块体的边界断裂上形成明显对比，意味着巴颜喀拉块体内部地壳应力的积累和集中不仅发生在边界断裂上，也发生在块体内部的次级断层上，块体边界断裂和内部断裂的构造运动和地震活动存在相互影响、相互调节的作用（Wang et al.，2 0 2 2).包括本文作者在内的众多学者对2 0 0 1年昆仑山地震的震后形变观测研究结果表明（Garthwaiteet al.,2013;Ryder et al.,2007;Wen et al.,2012;Zhaoetal.，2 0 2 2），昆仑山地震的震后形变主要发生在东昆仑断裂带南盘的巴颜喀拉地块内部，震后地球物理学报(Chinese J.Geophys.）2007)(图 2).玛尼段东昆仑断裂带西段11GPS observations(modified from Zhao et al.,2022)66卷mma-1(Van Der Woerd et al.2000),东段滑动速率快速减小，至10 2 E一10 3E塔藏断裂处减小至2mma-1左右（Diao et al.，2 0 19;K i r b y e t a l.，东昆仑断裂带中段东昆仑断裂带东段2001年昆仑山地震2021年玛多地震破裂范围破裂范围T100101102103104Bellet al.(2011)Garthwaite et al.(2013)Kirby et al.(2007)Lin andGuo(2008)巴颜喀拉地块的东向运动加速，而且这种加速随时间衰减缓慢，对块体东部形成持续的东向加载作用，并分摊到块体内部的多条次级断裂上，促进次级断层的应力积累，增加次级断层的地震危险性（万永革等,2 0 0 3;Shan et al.,2015;Xiong et al.,2010).但具体发生在哪条次级断裂上可能与次级断层的几何结构及发育成熟度有关.玛多地震的发生为我们研究巴颜喀拉块体内部次级断层与边界断裂的运动变形协同关系及东昆仑断裂东段的地震活动性特点等科学问题带来机遇。基于此，本文在利用大范围高分辨率InSAR数据研究2 0 0 1年昆仑山地震同震-震后形变及玛多地震同震形变的基础上，探讨昆仑山地震震后效应对巴颜喀拉块体的东向加载作用及其与玛多地震的关系，深化我们对巴颜喀拉块体及其主次断裂运动变形特征及未来地震危险性的认识和理解.1昆仑山地震震后形变场时空演化特征1.1基于ASAR长条带数据的震后形变场2001年Mw7.8昆仑山地震的震后形变效应及其时空演化特征是一个备受关注的科学问题，对这一问题的解决可为深人理解和认识东昆仑断裂带的活动习性、地震周期、流变结构及两盘介质差异等深Wang.and Shen(2019)Li et al.(2005)Zhenget al.(201%)Harkins et al.(2010)7期层次科学问题提供重要约束和依据.大面积、高分辨率、短周期重复观测的InSAR技术为研究大范围动态变化的震后形变时空特征提供了前所未有的技术途径(Garthwaite et al.,2013;Parizzi et al.,202l;Weiss et al.,2020).我们首先利用2 0 0 3 一2 0 10 年ENVISAT/ASAR长条带数据研究了2 0 0 1年昆仑山地震的震后形变演化特征.ASAR长条带数据扫描幅宽10 0km，利用5个相邻轨道长条带数据覆盖西至太阳湖、东至昆仑山口大约40 0 km长的整个地震破裂带，每个轨道上的可用数据条带有17 2 7 景不等，其中位于震中区的T133条带的数据量最大，有2 7景，可用于进行震后形变动态演化精细研究.采用SBAS InSAR技术进行时序数据处理，通过干涉网络构建、远场轨道相位误差拟合及时空滤波等方法精细去除各项残余误差，提取震后形变.计算了震后不同时段和整个观测时段的平均速率场，以及累计位移变化时间序列，分析了沿东昆仑断裂带不同段落震后形变衰减的时间演化过程及大范围震后形变场的整体分布变化特征，部分结果如图3一6.1.1.1震中区的震后形变演化特征图3为昆仑山地震震中区（T133条带）不同时段平均震后形变速率时空演化图像.T133条带数据观测时段较长（2 0 0 3一2 0 11年）且有较多、高质量干涉图可用，能获取比较完整和连续的震后形变时空演化过程.计算了6 个震后不同时段的形变速率场，分别是2 0 0 32 0 0 6、2 0 0 42 0 0 7、2 0 0 52 0 0 8、2 0 0 6 2009、2 0 0 7 2 0 10、2 0 0 8 2 0 11年.可以看出，除了断裂北侧地势陡峭的昆仑山和南侧部分山区、湖泊失相干比较严重外，其余地区震后形变信号空间连续性较好，基本勾勒出库赛湖段震后形变时空演化特征.震中区是同震位移最大的段落，震后形变效应也最为显著，但在2 0 0 3一2 0 11年的观测时段内，该段落处于相对稳定的变形过程，形变速率随时间的衰减并不显著，只在断层近场显示出较弱的南北盘不同衰减过程.南盘震后形变衰减比较缓慢，北盘的震后形变衰减相对比较显著，南盘近场的震后形变速率基本稳定在48 mma-1的范围内,不同阶段形变宽度也基本一致.北盘近场的震后形变衰减幅度在34mma-1,且随震后形变速率的衰减，北盘形变宽度也逐渐向南缩小，显示出明显的南北盘非对称衰减特征，图4为昆仑山地震震中区（T133条带）震后不同阶段3年时段跨断层累积形变量（2 0 0 3一2 0 10屈春燕等：2 0 0 1年昆仑山地震同震一震后效应与2 0 2 1年玛多地震关系探讨2745年），可以看出南盘巴颜喀拉块体的震后形变无论是形变宽度和形变量级均远大于北盘柴达木盆地，南盘形变宽度可达150 km，3年累积最大形变量可达4050 m m,出现在震中以南约2 0 km处；而北盘柴达木盆地形变宽度仅约7 0 km，3年累积最大形变量不超过30 mm.更为突出的是南盘在观测时段后期的3年累积位移大于早期的3年累积位移，意味着南盘形变衰减很缓慢，甚至没有衰减.1.1.2破裂带东端断层分叉处的震后形变演化特征图5和图6 所示为昆仑山地震破裂带东端（T 90 条带）震后不同时段平均形变速率时空演化的平面图和剖面图（剖面宽度为剖线两侧各5km）.此区域是东昆仑断裂带与南侧分支断裂一一昆仑山口一江错断裂的交叉连接处，昆仑山地震时同震破裂自西向东传播到该岔口后向南偏转沿昆仑山口一江错断裂传播约7 0 km后停止（Lasserre et al.，2005；K lin g e r e t a l.，2 0 0 6），说明了该区域主次断层运动的关联性和复杂性.图5中由于分时段形变速率计算中干涉图数量减少，相干性降低，导致不同时段速率图中缺值漏空区域较多，尤其是断裂带北侧昆仑山脉的大片山区和南侧的一部分地势陡峭区域，但整体上依然反映出断裂带南侧巴颜喀拉块体连续的震后形变信号，南盘的形变宽度和形变幅度均远大于北盘（图6），形变宽度达2 0 0 km，震后不同阶段3年时窗累积最大形变量达约7 10 mma-1，并不比震中区的T133条带小，只是震后形变随时间的衰减较震中区明显，但断裂南侧附近衰减依然缓慢，显示出持续震后形变对昆仑山口一江错断裂的长期加载作用.1.1.3昆仑山地震震后形变场大范围整体分布变化特征我们将5个相邻轨道的长条带InSAR速率场拼接，得到覆盖整个昆仑山地震破裂带的大范围震后形变场空间分布变化图像（图7），此图清晰地勾勒出昆仑山地震震后形变场在2 0 0 3一2 0 10 年观测时段内的整体分布形态、空间作用范围及沿断层/跨断层非均匀变化差异等信息.可以看出，震后形变影响范围很大，在东西向沿断层方向长达50 0 km，南北向跨断层方向宽达30 0 km，呈现出明显的南北盘不对称分布特征，南盘的形变范围和形变量级均明显大于北盘，南盘的形变宽度可达2 0 0 2 50km，是北盘的34倍，而且南盘震后形变场向南部远场区和破裂带东西两端的衰减都比较平稳且缓慢，特别是破裂带东部尾端震后形变依然很显著，在2746地球物理学报(Chinese J.Geophys.）92E93E(a)37N66卷94E92E(b)93E94E92E(c)93E94E36N35N34N-2003-2006(d)37N36N20042007e2005-2008)129-6335N-3-634N-920062009图3昆仑山地震震中区（T133条带）不同时段平均震后形变速率时空演化（a）2 0 0 32 0 0 6 年；（b）2 0 0 42 0 0 7 年；（c）2 0 0 52 0 0 8 年；（d）2 0 0 6 2 0 0 9年；（e）2 0 0 7 2 0 10 年；（f）2 0 0 8 2 0 11年；紫色粗线：同震地表破裂带；黑色虚线：图4剖面位置；沙滩球：2 0 0 1年Mw7.8昆仑山地震的震源机制解。Fig.3 Temporal and spatial evolution of the average postseismic displacement rate in dfferent periods aroundthe epicentral area of the Kunlun Mountain earthquake(T133)(a)20032006；（b)2 0 0 42 0 0 7；（c)2 0 0 52 0 0 8；（d)2 0 0 6 2 0 0 9；（e)2 0 0 7 2 0 10；（f)2 0 0 8 2 0 11;Pu r p le t h ic k lin e:c o s e is m icsurface rupture;Black dotted line:section position in Fig.4;Beach ball:focal mechanism of the Mw7.8 Kunlun Mountain earthquake.断层近场150 km范围以内甚至与震中区的形变量级相当.而断层北盘形变主要集中在震中附近区域，形变宽度窄而且向东西两侧衰减快.这些结果表明昆仑山地震的震后形变主要发生在断层南侧的巴颜喀拉块体，而且作用范围大，持续时间长，而对北侧柴达木盆地的影响是有限的，这意味着东昆仑断裂带南北两侧存在巨大的地壳结构和流变性质差异，且南侧巴颜喀拉块体的深部运动是主导性的.2007-2010震形变场（图7)的分布位置，很显然玛多地震正好位于东昆仑断裂带南侧约130 km、昆仑山地震震后东向加载作用最显著和持久的方向上，而且两个形变场边缘相距不超过30 0 km，玛多地震可能就位于昆仑山地震的震后加载区.另外，昆仑山地震同震破裂沿南侧分支断层一一昆仑山口一江错断裂终止，而这条次级断层的东延正是玛多地震的发震断层，20082011对比分析昆仑山地震震后形变场与玛多地震同-127期Fig.4 Cumulative fault-perpendicular displacement(T133)using a 3-year time window in different stagesThe gray thick line indicates the location of the East Kunlun fault,and the number on the right indicates the acquisition date38N(a)37N36N35N34N-图5昆仑山地震破裂带东端断层分叉处（T90条带)不同时段震后平均形变速率紫色线为昆仑山地震地表破裂带，黑色虚线为图6 形变剖面位置.Fig.5Average post-seismic displacement rate at fault bifurcation(T90 strip)aroundthe eastern end of the coseismic rupture zone in different periodsThe purple line is the surface rupture zone of the 2001 Kunlun Mountain earthquake,and the black dotted line is the position of profiles shown in Fig.6.屈春燕等：2 0 0 1年昆仑山地震同震一震后效应与2 0 2 1年玛多地震关系探讨-50-40-30-20-1001020304050250200150100500图4昆仑山地震震中区（T133条带）震后不同阶段3年时窗跨断层累积形变量（2 0 0 32 0 10 年）图中灰色粗线表示东昆仑断裂带位置，右侧数字表示干涉对的主辅影像拍摄日期.after the Kunlun Mountain earthquake epicentral area(20032010)of the primary and second images of the interferometric pairs.94E95E200620092747Displacement/mm20031030200606012003103020060914200401082006122820041118200706212004111820071004200503032007062120050303200711082005092920071004东昆仑断裂带2005092920080918200602162008022120060216200809182006021620090205200610192009020520061019200904162006122820090312200612282009062520061228-200909032007062120090730200706212010050620071004200910082007100420100715200711082010050694E95E(b)2007201094E(c)2008201195E129630-3-6-9-122748地球物理学报(Chinese J.Geophys.）1220062009年920072010年20082011年666卷东昆仑断裂带5343N-6-9-120图6昆仑山地震破裂带东端断层分叉处（T90条带）不同时段震后平均形变速率剖面灰色背景：剖面处地形；剖面位置见上图黑色虚线。Fig.6 Profiles of average postseismic displacement rate(T90)at the fault bifurcation aroundthe eastern end of the coseismic rupture in different periodsGray background:elevation variations along the section;See the black dotted line in Fig.5 for the position of profile.87E102E38N325010090E93E150Distance/km20096E25099E30038N36N34N36N东昆仑断裂带昆仑山日一江错断裂带2021M34N昆仑山地震震后形变/mm60-40-20020406032N87E图7 昆仑山地震震后累积形变场（2 0 0 32 0 10 年）与2 0 2 1年玛多地震同震形变空间分布关系蓝色箭头为震间GPS速率（WangandShen，2 0 2 0），灰色粗线为东昆仑断裂带，黑色虚线表示昆仑山口一江错断裂带，红色震源球分别为2 0 0 1年昆仑山地震和2 0 2 1年玛多地震震源机制解，黑色震源球为历史地震（M5.5)震源机制解。Fig.7 Spatial distribution relationship between cumulative postseismic displacement field of the Kunlun Mountainearthquake(20032010)and coseismic deformation of the 2021 Madoi earthquakeThe blue arrows indicate the interseismic GPS velocity(Wang and Shen,2020).The gray thick line indicates the East Kunlun fault.The black dottedline indicates the Kunlun Mountain pass-Jiangcuo fault.Two red beach balls are the focal mechanism of the 2001 Kunlun Mountain earthquake and进一步说明玛多地震的孕育发生与昆仑山地震及震后形变密切相关.1.2基于Sentinel-1宽幅数据的震后形变场1.2.1昆仑山地震震后形变场及应变场Sentinel-