尼泊尔MW7.9地震震后断层滑动时空特征分析

利用2015年尼泊尔MW7.9地震震后16个站531 d的连续观测数据，基于弹性位错理论，使用主成分反演方法（PCAIM）分析该地震震后断层滑动的时空变化特性。结果表明：1）震后断层滑动服从对数衰减模式；2）余滑主要分布在主震破裂带北部区域，最大滑动量达到20.6 cm，位于地表以下26.7 km处；3）余滑释放的地震矩为3.36×1020 N·m，对应矩震级约为M_W7.6。     

基于Sentinel-1A数据反演2015年新疆皮山MW6.5地震断层滑动分布

使用Sentinel-1A卫星观测的2015年新疆皮山地震震前震后一对SAR影像,提取该地震的同震形变场。结合Okada位错模型,在考虑渐变入射角和解缠基准偏差的情况下,反演本次地震的断层滑动分布。结果表明,皮山地震是小倾角逆冲型地震,断层破裂开始于地表以下约5 km处,最大滑动量为0.76 m,高滑区集中于深度9～14 km,累计释放地震矩达5.14×10¹⁸,相当于矩阵级M_W6.47。     

InSAR约束下的2021年西藏比如MW5.7地震断层几何参数及滑动分布

利用Sentinal-1A卫星SAR数据获取覆盖2021-03-19西藏比如M_W5.7地震的降轨同震视线向形变场,反演此次地震的发震断层几何参数及同震滑动分布。结果表明,此次地震破裂在一个之前未被探测到的隐伏断层上;破裂以正断为主,伴有少量左旋走滑分量;滑动主要分布在2.4～11.2 km深度处,且在6.4 km深度处滑动量达到最大,约为0.17 m;获取的地震矩约为3.6×10¹⁷ Nm,相当于M_W5.7。这些精确的断层参数及震源模型可为评估区域地震危险性及构造演化提供重要参考。     

2015年西藏定日MW5.7地震震源参数估计和静态应力触发研究

采用InSAR形变监测资料计算2015年西藏定日M_W5.7地震同震形变场,反演发震断层几何参数和滑动分布。在此基础上,研究尼泊尔M_W7.8主震对定日地震的静态库仑应力触发影响。综合分析地震滑动机制和构造特征,认为定日断层为西倾隐伏断层。反演结果表明,地震破裂相对集中,主要深度在6~9 km,破裂以正断滑动为主。发震断层走向约178°,倾角约 48°,破裂区长约5 km,宽约5 km,最大滑动量约0.2 m,释放的地震矩约3.7×10¹⁷ N·m,对应矩震级M_W5.6。尼泊尔主震同震库仑应力在定日地震震源处约为0.2 bar,造成藏南申扎-定结拉张地堑应变加载。     

2022年门源MW6.6地震震源破裂滑动分布

基于Sentinel-1和ALOS-2近场InSAR形变,利用有限断层方法反演2022-01-08门源M_W6.6地震的震源破裂滑动分布。结果显示,门源M_W6.6地震的最大滑动量约3.65 m,释放的地震矩约1.56×10¹⁹ Nm。破裂传播至地表,分东西两段,破裂滑移自震中沿断层向SEE和NWW向延伸,在震中的东南侧断层滑动量较大。有别于其他研究结果,本文认为门源地震震源破裂主要发生在东段的浅部,深度不超过8.0 km,西段的破裂最深可达15.0 km。同震形变整体上符合左旋走滑破裂特征,显示了青藏高原内部块体间的相对运动特征。     

利用InSAR技术研究2016年西藏定结Mw5.3地震发震构造

利用Sentinel-1A卫星数据和D-InSAR技术，获得2016-05-22西藏日喀则市定结县M_W5.3地震LOS方向的同震形变场图像，分析InSAR形变结果的变化特征，并以此为约束反演该地震的断层几何参数和同震滑动分布特征。结果表明，定结M_W5.3地震发震断层走向近SN，断面倾向E，倾角约43°，破裂长度约9 km，同震滑动主要集中在2~5 km深度范围内，以正断倾滑为主，最大滑动量为0.24 m，矩震级为5.3级。2016年定结地震发震构造是定结-申扎伸展断裂系中的一条新生盲断层。     

新疆于田MS7.3地震同震与震后形变机制研究

利用InSAR技术获取2008-03-21新疆于田MS7.3地震的同震和震后形变场。同震分布式滑动反演结果表明,同震断层最大滑动量达5.4 m,主要分布在南部断层的0~5 km深度附近,地震以正断错动为主,兼有左旋走滑分量。震后形变结果表明,发震断层北段两侧存在差异性运动,最大累积差异形变在震后782 d达15 cm。进一步分析表明,震后断层余滑可能是震后形变的主要机制。余滑反演结果表明,震后2 a断层余滑量相对较小,滑移区范围明显减小且均位于浅部区域,北部断层能量释放较彻底,南部仍存有少量能量,整体能量基本释放完。     

利用InSAR数据约束反演2017年Mw6.3精河地震同震破裂模型

利用Sentinel-1降轨卫星数据，研究2017-08-09新疆精河M_W6.3地震的同震形变与断层滑动破裂。考虑轨道和大气误差的影响，基于模型改正后的精细同震干涉图显示，视线向最大地表位移约为7 cm，未发现有地表破裂。此外，基于弹性半空间位错模型，采用两步法策略分别反演地震的均匀断层几何参数和同震滑动分布，结果显示，精河地震发震断层以逆冲运动为主，主破裂区位于8~17 km深度，最大滑移量为0.6 m，表明该地震是一次由逆冲型隐伏断层产生的单次破裂事件。InSAR结果给出的地震释放矩能量为2.91×10¹⁸ Nm，相当于矩震级M_W6.25，与USGS、CENC等机构给出的震级一致。     

基于GPS、InSAR和强震数据联合反演2017年九寨沟MS7.0地震同震滑动分布

以多源观测为约束,利用近场GPS、InSAR同震位移和强震观测数据联合反演2017年九寨沟M_S7.0地震发震断层的几何形状和破裂模型,并初步分析其地震成因。结果表明,九寨沟地震的发震断层走向154°,倾角80°,以左旋走滑为主,兼有少量正断分量;地震破裂主要集中在2~12 km深度范围,最大滑移量为1.6 m,位于地下6 km深处,断层在近地表没有滑移;地震释放的地震矩为5.59×10¹⁸ Nm,矩震级为M_W6.44。结合余震分布认为,该地震的发震断层为隐伏的虎牙断裂北段。     

顾及黏弹性松弛效应的2015年尼泊尔MW7.9地震震后早期余滑模型

处理尼泊尔境内及中国藏南的连续GPS数据，获得尼泊尔地震震后1 a的三维形变场。结果表明，震后形变主要发生在尼泊尔北部及中尼边境区域，水平形变以向南运动为主，垂直形变以隆升为主。采用有限断层模型反演的震后余滑分布范围十分广泛，主要集中在同震破裂下倾部分，向下延伸至30 km，平均余滑深度为20 km。震后1 a余滑释放的地震矩为8.8×1019 Nm，等价矩震级为MW7.3，占同震释放总能量的12%。扣除黏弹性松弛效应之后，反演的余滑分布更加集中于同震破裂下倾区域，断层模型底部无余滑，平均余滑深度减小为16 km。经黏弹性松弛效应改正后的余滑模型释放的地震矩为 5.7 ×1019Nm，仅为同震释放能量的8%左右。该余滑模型不但提高了数据的拟合精度，且其分布特征更接近应力驱动的模型。由于震后余滑发生在深部而非浅部未破裂区域，说明断层浅部仍处在闭锁状态，其未来的地震危险性仍需密切关注。     

Sentinel-1 SAR数据约束的2018年谢通门MW5.8地震同震破裂模型

通过Sentinel-1卫星升降轨数据获取谢通门地震的同震形变场,并基于均匀弹性半无限位错模型反演地震的同震滑动分布模型。InSAR同震形变场表明,升降轨视线向最大形变量分别为0.049 m和0.051 m,形变场长轴大致呈南北方向,位于甲岗-定结断裂西侧。通过对倾角和倾向进行格网搜索发现,西倾节面更可能为该地震的发震节面。反演结果表明,滑动分布主要位于2~10 km深度范围内,平均滑动量为0.02 m,最大滑动量为0.10 m,发震断层倾角为47°,平均滑动角为-81.60°,显示该地震以正倾滑动为主。大地测量数据约束的该地震震中为30.27°N、87.75°E,震源深度为6.58 km,释放地震矩为5.056×10¹⁷ Nm,对应矩震级为M_W5.7,与GCMT、USGS公布的震级基本一致。综合分析震中位置和滑动机制认为,甲岗-定结断裂的分支断层为本次谢通门地震的发震断层。     

基于升降轨Sentinel-1数据分析2017-11-12伊拉克MW7.3地震震源参数

利用Sentinel-1卫星提供的多个轨道SAR数据对2017-11-12伊拉克M_W 7.3地震进行研究，通过处理不同轨道数据并比较，获取震区可靠的同震形变场。利用降轨同震形变场，采用数据分辨率约束的四叉树方法降采样后作为约束条件，基于弹性半空间模型（Okada）反演断层几何参数及滑动分布。结果表明，本次地震位于扎格罗斯褶皱冲断带，产生的破裂带长度为44.5 km、宽度为16.1 km；发震机制为走滑逆冲断层，断层滑动角为140.68°，倾角为8°，最大滑动量为3.83 m，地震造成视线向最大抬升和最大下沉为50 cm和40 cm；此次地震震源深度为15.5 km，累计释放地震矩0.97×1020 Nm。     

基于Sentinel-1A的2021年双湖MW5.7地震同震破裂模型

利用Sentinel-1A SAR数据提取2021年西藏双湖县M_W5.7地震同震形变场及2.5D形变场,反演断层滑动分布模型。计算不同节面解为接收断层产生的库仑应力变化差异确定发震构造,并结合余震分布信息评估未来地震风险性。结果表明,地震震中为34.37°N、87.71°E,震源深度6.51 km,发震断层倾向东、走向33°、倾角50°、平均滑动角-74°,以倾滑为主兼有少量左旋走滑分量,最大滑动量0.26 m。短时间内,震区南部地震风险较小,北部则需要结合更多资料进一步分析。本次地震是在羌塘块体持续向东扩张的背景下,受EW向拉伸作用使得黄水湖正断层发生的一次弥散型变形活动,SN向地堑得到进一步扩张。     

玛多MW7.3地震同震滑动分布的三维有限元模拟分析

首次构建2021年玛多M_W7.3地震三维有限元模型,分析同震滑动分布特征。首先解算高精度GNSS同震形变观测数据;然后建立玛多地震三维有限元模型,并用弹性半空间Okada模型验证其准确性;最后以40个近场和远场GNSS同震形变观测数据为约束,利用最小二乘法反演玛多地震同震滑动分布模型,从而模拟同震位移。结果表明,玛多地震引起的破裂主要分布在野马滩和黄河乡附近,最大滑动值约为3.4 m。该结果与现场考察结果一致,能够很好地解释GNSS同震形变观测数据。     

2021年双湖MS5.8地震InSAR同震形变场及断层滑动分布反演

基于Sentinel-1A数据,利用合成孔径雷达差分干涉(DInSAR)技术提取视线向(LOS)同震形变场,利用基于序列蒙特卡罗采样的贝叶斯方法反演发震断层几何参数与断层滑动分布。结果表明,发震断层走向约47.810°±2.218°,倾角约36.664°±2.499°;震中位置为87.715°E、34.365°N,震源深度为5.673 km。同震滑动分布主要集中在3.9~7.5 km深度,最大滑动量约0.42 m,平均滑动角约-66.598°±3.258°。本次地震以正断为主,兼具少量左旋走滑性质,反演得到的矩震级为M_W5.5,略小于USGS提供的震级。     

2023年塔吉克斯坦MW7.2地震InSAR同震形变场与滑动分布

利用Sentinel-1A卫星SAR影像数据,对2023年塔吉克斯坦M_W7.2地震开展同震形变提取,基于弹性位错模型进行断层反演,并以本文反演得到的右旋节面解为接收面,计算不同深度的静态库仑应力。同震形变结果显示,升轨LOS向最大形变量达15 cm,降轨LOS向最大形变量达16 cm。断层反演结果表明,此次地震最优发震断层走向为131.1°、倾角为85.7°,同震主滑移区分布在深度10～30 km范围内,以右旋走滑为主,最大滑移位置位于地下约20 km深度处,滑移量为3.49 m,未破裂至地表,矩震级为M_W7.16。库仑应力结果显示,该区域库仑应力符合帕米尔高原已有的应力场及地质学研究结果,随着深度增加,其影响范围以发震断层为中心向外扩张,且自5 km深度往下,应力加载区逐渐侵蚀应力卸载区,并开始以加载区为主,在约10 km深度处开始发生余震活动,与本次发震断层相邻的2条断层未来短时间内地震风险性较小。     

2019-07 Ridgecrest地震及其余震的双差重定位分析

选取加州理工学院（Caltech）、美国地质调查局（USGS）及南加州地震数据中心（SCEDC）的震相观测报告和数字波形数据,结合波形互相关技术,利用双差定位方法单独和联合反演Ridgecrest M_W6.4前震、M_W7.1主震及震后1个月内M_W≥2.5余震的震源参数。精定位结果显示,2 098次地震事件形成4个事件簇,得到高达97%重定位率的震源参数;地震序列走时残差为65 ms,互相关走时残差为57 ms。精定位后震中位置线性特征更为显著,震群空间分布呈L型和T型特征,余震展布长约90 km,宽约25 km。定量分析地震活动与断层构造的关系,结果表明,Ridgecrest地震序列始于深部并向浅部区域传播,震源深度主要分布在4~8 km范围,呈分段条带特征,揭示Ridgecrest地震序列对该处规模不一的多断层构造组成的复杂断层系造成破坏,小湖断裂带为地震序列的发震断层。     

联合光学和雷达遥感数据的2021年玛多MW7.3地震同震形变提取

利用Sentinel-1升降轨数据,基于D-InSAR技术获取2021-05-22青海玛多M_W7.3地震LOS向形变场。由于地震破裂达到地表,沿断层剖面处出现了干涉条纹不连续、破碎重叠等现象,为准确提取地震发震断层的几何展布,联合地震前后的Landsat-8光学影像数据,基于频率域互相关算法提取地震水平向形变场及断层的几何展布特征。结果显示,玛多地震的地表破裂轨迹长达155.6 km,在首端及末端存在分支破裂,地表主体破裂带可分为3段,段间走向差异较大,整体呈NWW向展布。发震断层主要以左旋走滑为主,推测本次地震由昆仑山口-江错断裂控制。     

2021年漾濞MW6.1地震同震形变提取与断层源模型反演研究

基于D-InSAR技术,收集Sentinel-1卫星升、降轨及精密轨道数据,结合外部DEM,成功提取2021-05-21云南省大理州漾濞县M_W6.1地震的同震形变场。结果显示,升轨卫星视线向最大形变量级约为6.0 cm,降轨卫星视线向最大形变量级约为7.9 cm。以升、降轨InSAR观测结果为约束,对漾濞地震的断层几何参数及滑动分布进行联合反演发现,此次地震的最优发震断层走向为136.6°、倾角为83.1°,断层破裂主要集中在地下2~12 km深度处,最大滑动量约为0.45 m,位于地下7 km深度处,断层在近地表未出现大面积显著滑移现象,表明此次地震未破裂至地表。     

汶川地震震后余滑和震后粘弹性松弛数值模拟分析

采用2010~2015年汶川地区GNSS震后形变资料,利用有限元法建立三维震后粘弹性松弛模型,通过二维格网搜索获得龙门山断裂带上盘最佳弹性层厚度和中下地壳最佳粘滞系数,并分析汶川地震震后2~7 a粘弹性松弛影响下的震后形变特征;然后采用2008~2009年GNSS震后形变资料,根据最佳参数建立粘弹性松弛与余滑组合模型,并与单一余滑模型进行对比,分析汶川地震震后1 a内的形变特征。研究结果表明,根据模拟值计算得到青藏高原东部弹性层的最佳厚度为25 km,中下地壳的最佳粘滞系数为4.0×10¹⁸ Pa·s;震后1 a内组合模型的拟合效果优于单一余滑模型,其中余滑形变占主要成分。