随着合成孔径雷达干涉测量技术(Interferometric Synthetic Aperture Radar, InSAR)的飞速发展, 海量高质量的干涉图使得大面积地表形变监测成为可能。但星载SAR的标准景幅宽有限, 需要将多轨InSAR数据进行拼接来开展大范围地表形变的监测。聚焦于地震震间形变研究中广域形变场重建的应用需求, 文中首先基于模拟数据集定量分析了入射角对多轨道InSAR形变拼接的影响, 讨论了多轨道单一方向观测时InSAR形变场拼接中的主要误差来源; 提出了基于多项式估计的比值法修正入射角的方法; 面向大型走滑断裂, 讨论并分析了转换到地距向(Azimuth Look Direction, ALD)后进行拼接的应用场景和效果。最后, 以青藏高原东南部为实验区域, 利用GPS水平速度场进行参考基准校正, 分别使用比值法和转换到ALD方向对基于哨兵1号卫星重建的3个轨道InSAR数据进行拼接, 获得了大范围、 高精度的InSAR形变速率场。结果表明, 入射角的差异会造成相邻轨道同一区域InSAR震间形变速率场的差异, 当难以对InSAR形变场进行三维分解时, 文中提出的拼接策略可将多轨道的InSAR形变速率场的参考基准统一到同一空间参考框架下, 能够有效抑制入射角的影响, 实现广域InSAR震间形变速率场产品的拼接。

2021年5月21日云南漾濞县城西侧发生M_S6.4地震。 文中以欧洲航空局升、 降轨Sentinel-1 SAR为数据源, 基于D-InSAR技术获取了此次地震的InSAR同震形变场: 升轨LOS(Line of Sight)向最大形变量约为0.07m, 降轨最大LOS向形变量约为0.08m。 以升、 降轨同震形变场数据和GNSS数据为约束, 反演了断层滑动分布。 反演结果表明, InSAR/GNSS数据对断层倾向的约束能力较弱, NE倾向和SW倾向2种断层模型都能够在误差范围内拟合观测数据, 其中: SW倾向断层模型的最大滑动量约为0.8m, NE倾向断层模型的最大滑动量约为0.6m; 其余参数, 如滑动角为-170°、 倾角为80°、 矩震级为M_W6.07均一致, 破裂也都集中在地下2~10km深度范围内。 结合主震发生后3h内的余震精定位结果分析, 我们认为倾向SW的发震断层模型可能更符合实际情况。 根据断层滑动模型、 同震滑动分布并结合破裂运动学特征, 初步判定发震断层为维西-乔后-巍山断裂以西的一条次级断裂。 基于反演得到的断层滑动模型计算了漾濞地震破裂导致周边断裂的库仑应力变化情况, 结果表明漾濞地震对维西-乔后-巍山断裂、 红河中段断裂和红河北段断裂的南段具有显著的应力卸载作用, 初步认为漾濞地震之后周围断裂的地震危险性有所缓解。

文中基于D-InSAR技术, 利用欧空局C波段升降轨哨兵SAR数据, 获取了2021年5月22日青海玛多M_W7.3地震的InSAR同震形变场, 并对同震形变的空间特征、 量级和断层破裂的分段性进行了分析。 哨兵卫星的高质量观测数据清晰地描绘了玛多地震的地表破裂迹线, 地表破裂长度约210km。 为了进一步认识玛多地震断层深部的同震滑动分布特征和发震断层几何性质, 基于升、 降轨InSAR形变场及精定位余震数据确定的断层几何反演了同震滑动分布, 并基于同震库仑应力变化分析了本次地震对周边区域的应力扰动。 结果表明, 玛多M_W7.3地震发生在巴颜喀拉块体内部的一条次级断层上, 且与东昆仑断裂带的主断裂近平行, 结合野外考察、 地质资料和InSAR地表破裂迹线确定的发震断层为昆仑山口-江错断裂。 玛多地震产生的同震形变场空间影响范围广, 形变场的长轴呈NWW向, 升、 降轨观测的形变量符号相反, 结合哨兵数据的观测几何确定发震断层的运动性质以左旋走滑为主。 基于升、 降轨InSAR数据得到的最大视线向(Line of Sight, LOS)形变量约为0.9m。 同震滑动分布模型显示, 整体上断层的走向为276°, 倾角为80°, 倾向NE, 最大滑移量约为6m, 平均滑动角为4°, 矩震级为M_W7.45, 地震主体破裂发生在0~10km深度范围内, 同震破裂至地表, 与野外考察所观测到的广泛分布的同震地表破裂带一致。 以玛多-甘德断裂为接收断层的同震库仑应力模型显示, 玛多地震在玛多-甘德断裂西段附近产生了明显的应力降, 表明玛多地震释放了玛多-甘德断裂的库仑应力, 导致后者的地震危险性可能大大降低; 而在玛多-甘德断裂和昆仑山口-江错断裂交叉的区域存在应力加载效应, 但考虑到玛多地震的余震会释放多余的能量, 该区发生较大地震的危险性可能降低。 以东昆仑断裂为接收断层的同震库仑应力模型显示, 玛多地震在东昆仑断裂东段附近产生了明显的应力加载效应, 结合该段震间较高的应变率特征, 其地震危险性需要引起关注。