Application of High-precision UAV Aerial Survey in the Detailed Study of Surface Rupture of Maduo MS7.4 Earthquake in 2021
-
摘要: 北京时间2021年5月22日,青海省果洛州玛多县发生MS7.4地震,震中位于巴颜喀拉地块内部,根据震源机制解和野外地表破裂调查确定发震构造为以左旋走滑运动为主的江错断裂。本研究利用大疆Phantom 4 RTK无人机在震后采集大量地表破裂照片,采用集成SfM(Structure from Motion)算法的PhotoScan软件处理获得高分辨率DEM和正射影像,同时结合野外实地考察对研究区地表破裂的分布特征及断错地貌类型进行详细解译。利用基于MATLAB语言开发的位移测量软件LaDiCaoz,限定玛多地震在研究区产生的左旋走滑位移约为0.4 m。地表破裂精细化解译显示,在左旋右阶阶区发育小规模的挤压鼓包和里德尔共轭剪切破裂,在左旋左阶阶区发育走向为N40°~50°E,宽度达数十厘米的张裂缝带,指示发震构造的左旋走滑性质。本研究为震后基于无人机摄影测量技术快速提取地表破裂的定量参数和进行地表破裂精细化研究提供了可行、高效和科学的技术方法。Abstract: On May 22, 2021 (Beijing time), an Ms 7.4 earthquake occurred in Maduo county, Guoluo prefecture, Qinghai province. The epicenter was located in the Bayan Har Block, and the seismogenic structure was determined to be Jiangcuo fault with left-lateral strike-slip movement based on the focal mechanism solution and field investigation of surface rupture. In this study, the DJI Phantom 4 RTK UAV was used to collect a large number of photos of the surface rupture after the earthquake, and PhotoScan software integrated with the Structure from Motion (SfM) algorithm was used to obtain DEM and ortho images with high resolution. At the same time, combined with field investigation, the distribution characteristics of surface rupture and the types of displaced landform in the study area were interpreted in detail. Using the graphical dislocation measurement software LaDiCaoz developed based on Matlab, the left-lateral strike-slip displacement of Maduo earthquake in the study area is limited to 0.4 m. The results of the detailed interpretation of surface fractures show that there are small-scale compressional swelling and Riedel conjugate shear cracks in the left-lateral right-stepping zone and tensional fracture zones with a strike of N40°~50°E and a width of tens of centimeters in the sinistral left-steeping zone, indicating the sinistral strike-slip nature of the seismogenic fault. This study provides a feasible, efficient, and scientific-technical method for the rapid extraction of quantitative parameters of surface rupture and the detailed analysis of surface rupture based on UAV photogrammetry technology after large earthquakes.
-
Key words:
- UAV aerial survey /
- Surface rupture /
- Coseismic displacement
-
引言
跨断层形变监测是一种观测物理意义明确、测量精度高、便于组织实施以及数据产品稳定可靠的地震前兆监测方法之一(薄万举等, 1994, 1998;李杰等,2010;张晶等,2011;张希等,2014;马伶俐等,2015;李腊月等,2016;孙启凯等,2017)。通过多期次、短距离的跨断层水准和基线测量,能够准确记录断层活动的位移量,监测活动断裂的运动特征,更好地反映断裂带及其所处区域的形变和运移特征(牛安福等,2013)。针对川滇地块的地壳形变、运动学特征及其形成机制开展的研究工作产出了大量成果,但是,目前关于这些问题也存在不同认识(乔学军等,2004;武艳强等,2011;洪敏等,2014;吴中海等,2015;白玉柱等,2016)。本文利用川滇地块南部附近的跨断层观测资料,对区域断裂带的现今形变变化进行分析,研究川滇地块南部的整体运动学特征,综合提取断层活动的异常信息,为区域地震危险性分析提供一定的参考。
1. 区域地质概况
川滇地块位于青藏高原的东缘,南北地震带的南段,大地构造上处于扬子板块、松潘-甘孜地块、兰坪-思茅地块以及冈底斯-腾冲地块等构造单元的拼合汇聚地带,在地质历史上经历了多次不同时期不同规模的构造运动与岩浆活动,构造属性复杂。该区域地震活动强烈而频繁,地震与构造运动之间存在复杂的因果关系,是研究地震活动与断层运动关系的热点地区(图 1)。
自新生代以来,印度板块与亚欧板块的相互碰撞拼贴,造成地壳强烈缩短增厚,青藏高原得以快速隆升(Tapponnier等,2001;许志琴等,2011)。伴随着两大板块的俯冲拼合,高原也发生着同构造挤压伸展,物质和能量向北西、北东、东、南东等多个方向转移释放。这种持续的俯冲碰撞使得两大板块拼合区域的喜马拉雅造山带、雅鲁藏布江缝合带以及高原内部地块与周缘区域构造活动强烈,并形成一系列大断裂带,如北西走向的鲜水河断裂、红河断裂及近南北走向的小江断裂、安宁河-则木河断裂带等。很多学者对此区域开展了研究,并取得了大量成果(Wang等,1998;Burchfiel等,2003;陈长云等,2016)。其中,处于川滇地块南部的活动断裂,是本文关注的重点。
2. 跨断层监测概况
川滇地块南部断裂所在的云南地区的跨断层监测场地从20世纪80年代开始观测,至今30余年,主要为流动观测场地,共有11处,分别是楚雄、峨山、建水、剑川、丽江、石屏、通海、下关、永胜、宜良和羊街,主要分布于金沙江-红河断裂带、小江断裂带、南华-楚雄-建水断裂带及其次级断裂(剑川断裂、洱海断裂、通海-峨山断裂、建水-石屏断裂、丽江断裂)上(表 1),多为综合观测场地。每个场地包括2条短水准和2条短基线测线,分别各有1条直交和斜交测线,通过这种组合方式来共同监测场地所跨断层的三维运动变化。从场地观测历史来看,该区域的跨断层观测基本是每月1期,个别年份每2个月观测1期,数据连续且完整。
表 1 研究区跨断层监测场地概况Table 1. Basic information of the cross-fault monitoring sites in the area场地 经度/°E 纬度/°N 所跨断层 断层产状 断层性质 楚雄 101.52 25.03 楚雄断裂 NE∠58° 右旋挤压 峨山 102.50 24.12 通海-峨山断裂 S∠60° 右旋挤压 建水 102.81 23.63 石屏-建水断裂 S∠70° 左旋挤压 剑川 99.93 26.63 剑川断裂 W∠75° 左旋拉张 丽江 100.25 26.88 丽江断裂 SW∠70° 左旋挤压 石屏 102.45 23.77 石屏-建水断裂 NE∠70° 右旋拉张 通海 102.72 24.00 通海-峨山断裂 NE∠68° 右旋挤压 下关 100.30 25.62 红河断裂 NE∠60° 右旋拉张 永胜 100.72 26.70 程海断裂 W∠35° 左旋挤压 宜良 103.03 25.02 马街-南羊街断裂 SE∠68° 左旋挤压 羊街 103.10 25.55 小江断裂西支 SE∠73° 左旋挤压 3. 数据计算及分析
将收集到的研究区30余年的跨断层观测数据进行预整理,对数据格式做统一标准化处理。首先,消除跨断层观测资料中明确因干扰引起的突跳变化,对观测资料进行连续性检验,剔除重复的数据。然后,对于日常观测中因各种原因造成缺测的数据采用线性插值法修补,得到周期统一的、可靠的形变数据。对数据的处理分析采用以下方法:
云南区域内的跨断层场地同时布设有两条同桩的基线和水准测线,可采用下面3个公式来计算断层活动的三维分量:
$$\Delta S=\frac{\Delta {{L}_{1}}\text{sin}{{\alpha }_{2}}-\Delta {{L}_{2}}\sin {{\alpha }_{1}}}{\sin ({{\alpha }_{2}}-{{\alpha }_{1}})}$$ (1) $$\Delta R=\frac{\Delta {{L}_{1}}-\Delta S\cos {{\alpha }_{1}}}{\text{sin}{{\alpha }_{1}}}$$ (2) $$\Delta H=\frac{1}{2}(\Delta {{H}_{1}}+\Delta {{H}_{2}})$$ (3) 其中,$\Delta H$为水准垂直变化量(以上盘相对于下盘的下降为正),$\Delta L$为基线伸缩量(以伸长为正),$\Delta S$为断层活动水平走滑量(两盘做顺扭运动为正),$\Delta R$为断层活动水平张压量,${{\alpha }_{1}}$和${{\alpha }_{2}}$分别代表由断层线方向逆时针转动至与两条测线方向重合或者平行时的角度。当$\Delta S$为正值时,表示右旋活动,$\Delta R$为正值时,表示张性活动,若均为负值则表示断层活动性质相反(薄万举等,1998)。
利用MATLAB对计算得出的数据绘制曲线图,单幅曲线图上方均标注有其代表的场地名及运动方式(图 2)。通过曲线图可以发现跨断层场地所观测的不同断裂表现出了跟地质背景相同或者不同的运动学特征,其中丽江、石屏和下关场地所跨断层的现今活动特征跟地质背景一致;楚雄、剑川和永胜场地所跨断层在走滑和拉张特征上都表现出跟背景特征不一致的现象;峨山、建水和通海场地所观测的断层在走滑性质上和地质背景相符,而在张压性质上相反。根据现今运动特征,可将场地分为3种类型:
(1)运动特征与地质背景一致。丽江场地所跨的丽江次级断裂,在地质背景上是一处表现为左旋走滑特征的活动断裂,兼具挤压特点,现今计算的三维活动量,依然表现出相同特征,其中断层的水平张压量在2012年后出现大幅异常变化。石屏场地所跨的石屏-建水断裂,在地质背景上是一处右旋拉张特征的活动断裂,现今右旋运动速率在2013年后虽然有所放缓,但是仍然与地质背景一致。下关场地所处的是一条表现为右旋走滑兼有拉张特征的深断裂,2003年水平走滑量曾出现右旋剪切转折,现今曲线重回向上趋势,继续保持右旋走滑特征。羊街和宜良场地所处的小江断裂带西支断裂,现今数据仍显示出以继承性的左旋走滑兼具挤压特征为主。
(2)运动特征与地质背景不一致。楚雄场地所跨断层为具有右旋挤压背景特征的楚雄断裂,在2008年后现今运动特征以强烈左旋拉张为主,但自2014年左旋和拉张的运动速率都出现大幅放缓的异常现象。剑川场地所处断裂为红河断裂带北段的剑川断裂,是一处左旋走滑兼具拉张特征的活动断裂,现今三维运动特征表现为右旋走滑兼具挤压特征,其中的垂直位移量自2015年至今出现了大幅抬升的异常变化;水平走滑量自2011年后也出现大幅抬升,右旋运动加速。永胜场地所监测的断裂是具有左旋挤压特征的程海深断裂,从计算所得的三维活动量曲线可以看出,该断裂现在运动学特征为右旋拉张,其垂直位移量曲线从2014年开始大幅异常抬升,水平走滑量自1994年开始稳步抬升。
(3)与地质背景相比,走滑特征一致,张压特征相反。峨山场地和通海场地共同监测曲江断裂的不同部位,该断裂在地质背景上具有右旋挤压特征,两者现今走滑特征与背景一致,但是张压量曲线表现出相反的拉张运动的异常特点。其中,峨山场地走滑分量在1998年后大幅抬升,之后趋于稳定的右旋运动状态;通海场地的垂直位移量在1992年后出现转折,断层表现为张性活动,走滑分量在2003年大幅下降,断层剪切活动速率减小,之后趋于稳定,依然显示右旋运动。建水场地所跨的石屏-建水断裂的建水段,断层性质为左旋挤压,从计算得出的水平走滑量参数来看,现今的走滑运动特征与地质背景一致,垂直位移量在1999年后发生趋势改变,断层活动由压性转变为相反的拉张特征。
4. 讨论与结论
在印度洋板块和亚欧板块的互相作用下,巨大的青藏高原造就了地球的“第三极”,这种作用力至今仍在持续。板块的碰撞拼合产生了强烈的构造活动,许多大规模断裂带应运而生(Xu等,2016),同时,也孕育了一系列强震。对此一些学者做过许多研究,并取得大量有价值的科研成果(Tapponnier等,2001;许志琴等,2011;陈长云等,2016)。尽管目前对青藏高原的陆内侧向运动形式、时间和逃逸量存在着争议,但是对位于其东缘的川滇地块发生的晚新生代南东向的逃逸和围绕东喜马拉雅构造结顺时针的旋转有着一致的认识(乔学军等,2004;武艳强等,2011;洪敏等,2014;吴中海等,2015)。
通过对位于川滇地块南部主要断裂上的跨断层场地三维形变量的分析,可以看出部分断裂表现出继承性特征,如现今下关和石屏场地的三维活动量表现出跟以往相同的持续平稳的活动速率。而另一部分断裂表现出了异常于背景活动特征的现象,如楚雄、永胜和剑川场地在走滑分量上都表现出了反向活动特点。研究区具有以水平运动为主、垂直运动为辅的背景特征,现今看来,依然如此。通过计算分析得到的宜良、羊街、楚雄以及丽江场地的形变结果来看,羊街、宜良场地所控制的小江断裂带分支为川滇地块东部边界断裂,现呈现左旋走滑兼具压性的特征;楚雄场地监控的楚雄-南华断裂位于川滇地块的南边部,表现出左旋拉张的运动特征;丽江场地所处的丽江-小金河断裂作为川滇地块南北部的分界断裂,现今以左旋挤压运动为主。可以看出,川滇地块南部区域整体上表现出以水平顺时针旋转为主的运动方式。根据GPS观测资料得到的该区域现今水平运动速度场也显示了相同的运动特点(乔学军等,2004;武艳强等,2011;洪敏等,2014;王伶俐等,2015)。造成这种运动方式的深层次动力学机制,应该源于印度洋板块俯冲至亚欧板块下,碰撞拼合后巨大高耸的青藏高原崛起,在超级重力作用下,物质开始向周缘,尤其是向东缘流动,在遇到横亘的北东向龙门山造山带时受阻,流动速率下降,被消化、吸收,再次转移(张家声等,2003;崔效锋等,2006)。其中,相当一部分能量和物质向南东运移,规模相对较小的川滇活动地块随之发生南东向旋转,同时起着调节周围各大板块和地块之间运动平衡的“交通环岛”的作用。另外,印度洋板块至今持续强烈地北向运动,跟扬子板块形成了大尺度上的右旋运动,也推动着川滇地块南东向旋转运移。
根据观测场地所跨断层的三维运动特征及计算所得曲线可以看出,在左旋的小江断裂带、右旋的红河断裂带和右旋的南华-楚雄-建水断裂带共同汇聚的峨山、通海和建水场地区域,原先的地壳形变表现出压性特征,现今由于断裂性质的异常变化,该区域场地集体表现出拉张的运动特点。滇西北区域的永胜、丽江和剑川场地所跨断层也表现出了一定的异常运动特征,剑川断裂和程海断裂呈水平反向运动,除此之外,包括丽江场地所跨的丽江断裂在内,3条断裂的垂向运动特征也存在大幅异常变化。云南地区中强震发生前的一定时段或震时,区域内跨断层场地观测到的水准和基线数据多数会出现较明显的变化,且部分场地形变量会发生速率的急变或转折,如文山5.3级、孟连西7.3级、姚安6.5级、丽江7.0级地震震前和震时,部分场地就表现出了加速继承性运动和转折性运动的特征(张兴华等,1996;付虹等,1997;李忠华等,1998;施顺英等,2005;李瑞莎等,2009;刘强等,2010;王永安等,2011)。同时,云南地区2015年至今存在长时间缺少4级以上中强地震的背景,不符合区域以往的地震活动特征,现今应力积累较高。自2016年云龙5.0级地震后至今,区域中强震活动已平静300天左右;滇西北地区自2001年永胜6.0级地震至今,已14.7年未发生6级以上地震;滇中南的石屏、建水等区域从2001年江川5.1级地震起算,5级地震平静期已长达15.5年。基于此,该区域在中短期尺度可能存在一定中强震的危险。
综上所述,得到以下几点结论:
(1)川滇地块南部跨断层场地的形变特征表现出继承和异常于地质背景特征的现象。① 丽江、石屏、下关、宜良和羊街场地所跨断层的现今活动特征与地质背景一致。丽江断裂和小江断裂表现出持续地左旋挤压运动特征;石屏-建水断裂石屏段表现出右旋拉张特征;红河断裂表现出跟地质背景相同的右旋拉张特征。② 楚雄、剑川和永胜场地所跨断层在走滑和张压特征上都出现与背景特征相背的现象。楚雄断裂、剑川断裂和程海断裂均表现出相反的运动特征。③ 峨山、建水和通海场地所监测的断层现今活动状态与地质背景相比,在走滑性质上相符,张压性质上相反。通海-峨山断裂现今运动特征,表现出与地质背景相同的右旋走滑和与地质背景相反的拉张特征。石屏-建水断裂建水段现今依然呈左旋走滑运动,但是张压性质发生反转,表现出张性运动的活动特点。
(2)川滇地块南部整体上表现出以水平顺时针旋转为主的运动方式。
(3)根据跨断层观测资料所得出的断层三维运动特征及曲线可知,滇西北及滇中的峨山、通海、石屏和建水等区域地壳形变特征存在异常,应关注该区域中短期尺度内可能存在的中强震危险。
致谢: 对编辑和审稿老师的辛苦审阅和宝贵建议,表示衷心感谢。 -
-
[1] 陈立春, 王虎, 冉勇康等, 2010. 玉树MS7.1级地震地表破裂与历史大地震. 科学通报, 55(13): 1200—1205.Chen L. C., Wang H., Ran Y. K., et al., 2010. The MS7.1 Yushu earthquake surface rupture and large historical earthquakes on the Garzê-Yushu Fault. Chinese Science Bulletin, 55(31): 3504—3509. (in Chinese) [2] 程佳, 徐锡伟, 甘卫军等, 2012. 青藏高原东南缘地震活动与地壳运动所反映的块体特征及其动力来源. 地球物理学报, 55(4): 1198—1212 doi: 10.6038/j.issn.0001-5733.2012.04.016Cheng J. , Xu X. W. , Gan W. J. , et al. , 2012. Block model and dynamic implication from the earthquake activities and crustal motion in the southeastern margin of Tibetan Plateau. Chinese Journal of Geophysics, 55(4): 1198—1212. (in Chinese) doi: 10.6038/j.issn.0001-5733.2012.04.016 [3] 邓起东, 2007. 中国活动构造图. 北京: 地震出版社.Deng Q. D., 2007. Map of active tectonics in China. Beijing: Seismological Press. (in Chinese) [4] 邓起东, 闻学泽, 2008. 活动构造研究——历史、进展与建议. 地震地质, 30(1): 1—30 doi: 10.3969/j.issn.0253-4967.2008.01.002Deng Q. D. , Wen X. Z. , 2008. A review on the research of active tectonics——history, progress and suggestions. Seismology and Geology, 30(1): 1—30. (in Chinese) doi: 10.3969/j.issn.0253-4967.2008.01.002 [5] 邓起东, 高翔, 陈桂华等, 2010. 青藏高原昆仑—汶川地震系列与巴颜喀喇断块的最新活动. 地学前缘, 17(5): 163—178Deng Q. D. , Gao X. , Chen G. H. , et al. , 2010. Recent tectonic activity of Bayankala fault-block and the Kunlun-Wenchuan earthquake series of the Tibetan Plateau. Earth Science Frontiers, 17(5): 163—178. (in Chinese) [6] 邓起东, 程绍平, 马冀等, 2014. 青藏高原地震活动特征及当前地震活动形势. 地球物理学报, 57(7): 2025—2042 doi: 10.6038/cjg20140701Deng Q. D. , Cheng S. P. , Ma J. , et al. , 2014. Seismic activities and earthquake potential in the Tibetan Plateau. Chinese Journal of Geophysics, 57(7): 2025—2042. (in Chinese) doi: 10.6038/cjg20140701 [7] 高帅坡, 2017. 基于无人机摄影测量技术的活动构造定量参数提取研究. 北京: 中国地震局地质研究所.Gao S. P. , 2017. A quantitative parameters extraction study of active tectonics based on UAV photogrammetry technology. Beijing: Institute of Geology, China Earthquake Administrator. (in Chinese) [8] 荆帅军, 帅向华, 甄盟, 2019. 基于无人机倾斜影像的三维建筑物震害精细信息提取. 地震学报, 41(3): 366—376Jing S. J. , Shuai X. H. , Zhen M. , 2019. Fine information extraction of 3D building seismic damage based on unmanned aerial vehicle oblique images. Acta Seismologica Sinica, 41(3): 366—376. (in Chinese) [9] 康文君, 徐锡伟, 于贵华等, 2020.2种基于Matlab平台的断层位移测量软件对比分析——以阿尔金断裂东段为例. 地震地质, 42(3): 732—747 doi: 10.3969/j.issn.0253-4967.2020.03.013Kang W. J. , Xu X. W. , Yu G. H. , et al. , 2020. Comparison study of two kinds of codes to measure fault-offsets based on Matlab: a case study on eastern Altyn Tagh Fault. Seismology and Geology, 42(3): 732—747. (in Chinese) doi: 10.3969/j.issn.0253-4967.2020.03.013 [10] 李金香, 常想德, 姚远等, 2019. 无人机技术在新疆塔县地震的应用及实现. 华南地震, 39(3): 57—64Li J. X. , Chang X. D. , Yao Y. , et al. , 2019. Application and realization of UAV technology in the earthquake of Tajik Autonomous County of Taxkorgan, Xinjiang. South China Journal of Seismology, 39(3): 57—64. (in Chinese) [11] 李云, 徐伟, 吴玮, 2011. 灾害监测无人机技术应用与研究. 灾害学, 26(1): 138—143 doi: 10.3969/j.issn.1000-811X.2011.01.026Li Y. , Xu W. , Wu W. , 2011. Application research on aviation remote sensing UAV for disaster monitoring. Journal of Catastrophology, 26(1): 138—143. (in Chinese) doi: 10.3969/j.issn.1000-811X.2011.01.026 [12] 李智敏, 李文巧, 李涛等, 2021.2021年5月22日青海玛多MS7.4地震的发震构造和地表破裂初步调查. 地震地质, 43(3): 722—737Li Z. M., Li W. Q., Li T., et al., 2021. Seismogenic fault and coseismic surface deformation of the Maduo MS7.4 earthquake in Qinghai, China: a quick report. Seismology and Geology, 43(3): 722—737. (in Chinese) [13] 刘超, 雷启云, 余思汗等, 2021. 基于无人机摄影测量技术的地震地表破裂带定量参数提取——以1709年中卫南M7½地震为例. 地震学报, 43(1): 113—123 doi: 10.11939/jass.20200039Liu C. , Lei Q. Y. , Yu S. H. , et al. , 2021. Using UAV photogrammetry technology to extract the quantitative parameters of earthquake surface rupture zone: a case study of the southern Zhongwei M7½ earthquake in 1709. Acta Seismologica Sinica, 43(1): 113—123. (in Chinese) doi: 10.11939/jass.20200039 [14] 刘静, 陈涛, 张培震等, 2013. 机载激光雷达扫描揭示海原断裂带微地貌的精细结构. 科学通报, 58(1): 41—45 doi: 10.1360/972012-1526Liu J. , Chen T. , Zhang P. Z. , et al. , 2013. Illuminating the active Haiyuan fault, China by airborne light detection and ranging. Chinese Science Bulletin, 58(1): 41—45. (in Chinese) doi: 10.1360/972012-1526 [15] 陆博迪, 孟迪文, 陆鸣等, 2011. 无人机在重大自然灾害中的应用与探讨. 灾害学, 26(4): 122—126 doi: 10.3969/j.issn.1000-811X.2011.04.024Lu B. D. , Meng D. W. , Lu M. , et al. , 2011. Application and exploration of unmanned aerial vehicle in major natural disasters. Journal of Catastrophology, 26(4): 122—126. (in Chinese) doi: 10.3969/j.issn.1000-811X.2011.04.024 [16] 潘家伟, 白明坤, 李超等, 2021.2021年5月22日青海玛多MS7.4地震地表破裂带及发震构造. 地质学报, 95(6): 1655—1670Pan J. W. , Bai M. K. , Li C. , et al. , 2021. Coseismic surface rupture and seismogenic structure of the 2021-05-22 Maduo (Qinghai)MS7.4 earthquake. Acta Geologica Sinica, 95(6): 1655—1670. (in Chinese) [17] 孙稳, 何宏林, 魏占玉等, 2019. 基于无人机航测获取高分辨率DEM数据的断层几何结构精细解译与分析——以海原断裂唐家坡为例. 地震地质, 41(6): 1350—1365 doi: 10.3969/j.issn.0253-4967.2019.06.003Sun W. , He H. L. , Wei Z. Y. , et al. , 2019. Interpretation and analysis of the fine fault geometry based on high-resolution DEM data derived from UAV photogrammetric technique: a case study of Tangjiapo site on the Haiyuan Fault. Seismology and Geology, 41(6): 1350—1365. (in Chinese) doi: 10.3969/j.issn.0253-4967.2019.06.003 [18] 铁瑞, 王俊, 贾连军等, 2016. 强震地震数据统计及其地表破裂特性研究. 世界地震工程, 32(1): 112—116Tie R. , Wang J. , Jia L. J. , et al. , 2016. Data statistics of strong-moderate earthquake and characteristics research of ground rupture. World Earthquake Engineering, 32(1): 112—116. (in Chinese) [19] 闻学泽, 2018. 巴颜喀拉块体东边界千年破裂历史与2008年汶川、2013年芦山和2017年九寨沟地震. 地震学报, 40(3): 255—267 doi: 10.11939/jass.20170211Wen X. Z. , 2018. The 2008 Wenchuan, 2013 Lushan and 2017 Jiuzhaigou earthquakes, Sichuan, in the last more than one thousand years of rupture history of the eastern margin of the Bayan Har block. Acta Seismologica Sinica, 40(3): 255—267. (in Chinese) doi: 10.11939/jass.20170211 [20] 徐锡伟, 陈文彬, 于贵华等, 2002.2001年11月14日昆仑山库赛湖地震(MS8.1)地表破裂带的基本特征. 地震地质, 24(1): 1—13Xu X. W., Chen W. B., Yu G. H., et al., 2002. Characteristic features of the surface ruptures of the Hoh Sai Hu (Kunlunshan) earthquake (MS8.1), northern Tibetan Plateau, China. Seismology and Geology, 24(1): 1—13. (in Chinese) [21] 徐锡伟, 陈桂华, 2018. 活动断层避让问题探讨与建议. 城市与减灾, (1): 8—13. doi: 10.3969/j.issn.1671-0495.2018.01.004 [22] 姚鑫, 周振凯, 李凌婧等, 2017.2017年四川九寨沟Ms7.0地震InSAR同震形变场及发震构造探讨. 地质力学学报, 23(4): 507—514Yao X. , Zhou Z. K. , Li L. J. , et al. , 2017. InSAR co-seismic deformation of 2017 Ms7.0 Jiuzhaigou earthquake and discussions on seismogenic tectonics. Journal of Geomechanics, 23(4): 507—514. (in Chinese) [23] 张培震, 邓起东, 张国民等, 2003. 中国大陆的强震活动与活动地块. 中国科学D辑, 33(S1): 12—20.Zhang P. Z., Deng Q. D., Zhang G. M., et al., 2003. Active tectonic blocks and strong earthquakes in the continent of China. Science in China Series D: Earth Sciences, 46(S2): 13—24. [24] 张培震, 徐锡伟, 闻学泽等, 2008.2008年汶川8.0级地震发震断裂的滑动速率、复发周期和构造成因. 地球物理学报, 51(4): 1066—1073Zhang P. Z. , Xu X. W. , Wen X. Z. , et al. , 2008. Slip rates and recurrence intervals of the Longmen Shan active fault zone and tectonic implications for the mechanism of the May 12 Wenchuan earthquake, 2008, Sichuan, China. Chinese Journal of Geophysics, 51(4): 1066—1073. (in Chinese) [25] 中国地震局地球物理研究所, 2021. 2021年5月22日青海果洛州玛多县7.4级地震应急科技支撑简报. (2021-05-22)[2021-11-01]. http://www.cea-igp.ac.cn/kydt/278249.html. [26] Deng Q. D. , Wu D. M. , Zhang P. Z. , et al. , 1986. Structure and deformational character of strike-slip fault zones. Pure and Applied Geophysics, 124(1—2): 203—223. [27] Fossen H. , 2010. Structural geology. Cambridge: Cambridge University Press. [28] Johnson K. , Nissen E. , Saripalli S. , et al. , 2014. Rapid mapping of ultrafine fault zone topography with structure from motion. Geosphere, 10(5): 969—986. doi: 10.1130/GES01017.1 [29] Klinger Y. , Etchebes M. , Tapponnier P. , et al. , 2011. Characteristic slip for five great earthquakes along the Fuyun fault in China. Nature Geoscience, 4(6): 389—392. doi: 10.1038/ngeo1158 [30] Lin Z. , Kaneda H. , Mukoyama S. , et al. , 2013. Detection of subtle tectonic–geomorphic features in densely forested mountains by very high-resolution airborne LiDAR survey. Geomorphology, 182: 104—115. doi: 10.1016/j.geomorph.2012.11.001 [31] Oskin M. E. , Arrowsmith J. R. , Corona A. H. , et al. , 2012. Near-field deformation from the El Mayor-Cucapah earthquake revealed by differential LIDAR. Science, 335(6069): 702—705. doi: 10.1126/science.1213778 [32] USGS (United States Geological Survey), 2021. M7.3-southern Qinghai, China-Moment tensor. (2021-05-21)[2021-11-01]. https://earthquake.usgs.gov/earthquakes/eventpage/us7000e54r/executive. [33] Wallace R. E. , 1998. The geology of earthquakes. Eos, Transactions American Geophysical Union, 79(9): 115. [34] Wang S. F., Fan C., Wang G., et al., 2008. Late Cenozoic deformation along the northwestern continuation of the Xianshuihe fault system, Eastern Tibetan Plateau. GSA Bulletin, 120(3—4): 312—327. [35] Zhou Y. , Parsons B. , Elliott J. R. , et al. , 2015. Assessing the ability of Pleiades stereo imagery to determine height changes in earthquakes: a case study for the El Mayor‐Cucapah epicentral area. Journal of Geophysical Research, 120(12): 8793—8808. [36] Zielke O. , Arrowsmith J. R. , Ludwig L. G. , et al. , 2010. Slip in the 1857 and earlier large earthquakes along the Carrizo Plain, San Andreas fault. Science, 327(5969): 1119—1122. doi: 10.1126/science.1182781 [37] Zielke O. , Arrowsmith J. R. , 2012. LaDiCaoz and LiDARimager—MATLAB GUIs for LiDAR data handling and lateral displacement measurement. Geosphere, 8(1): 206—221. doi: 10.1130/GES00686.1 -