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条断层未来短时间内地震风险性较小。     

基于GPS的2021年玛多MW7.4地震同震静、动态形变场

利用青海省CORS网GPS连续观测数据获取2021年玛多M_W7.4地震的同震形变场以及震时地表运动状态。结果表明,断层南北两侧的近场测站分别展现出南东向和北西向运动,符合左旋走滑机制。同震形变集中在震中距300 km范围内,震中距150 km以内的站点均能反映出cm级的位移,最大水平向位移为28.3 cm(JDUO站)。高频GPS动态形变提取出的永久位移与静态解算结果相当,其动态波形最大峰值为49.9 cm(KANQ站)。依据震中距和波形初动时间估计得到地震波速为2.8~4.9 km/s,断层东端的站点估算速度大于其他站点,可能与此次破裂的方向性或者破裂传播速度有关。依据震级经验公式,动态波形估算震级在M6.8~7.6区间,拟合平均震级为M7.35±0.15。若在实时条件下,震后70 s内可得到稳定的震级。     

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

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

改则Ms6.9地震InSAR同震形变场模拟

以InSAR同震形变场为约束,采用各向同性弹性半空间位错模型,利用试错法对2008年西藏改则县扎西错\%M\%s6.9地震断层几何形状进行格网搜索寻优模拟。观测值离散化得到1 010个数据总集,并重采样80次得到80个数据子集,每个子集之间存在微小的噪声差异,采用Bootstrap方法进行80次搜索寻优。首先用单断裂矩形断层模型进行模拟,最大残差为30 cm。然后用双断裂矩形断层模型进行模拟,最大残差为15 cm。双断裂模型反映了观测形变场的基本特征,并与地震资料符合。根据80次搜索寻优得到主断裂倾滑量95%置信区间为1.6~2.4 m。模拟表明,地表变形是主、余震两次破裂的结果。实际观测和数值模拟结果表明,InSAR可以作为监测青藏高原大陆动力学的有效技术。     

2021年云南漾濞Mw6.4级地震的InSAR监测与反演

2021年5月21日,云南漾濞地区发生M_w 6.4级地震,然而发震区的活动构造前期研究薄弱。针对此次地震开展监测与震源参数反演研究,对深入掌握漾濞地区及其邻近区域的孕震及发震机理和构造活动特征具有十分重要的科学意义。利用欧洲航天局两个轨道的Sentinel-1A SAR卫星影像获取了升、降轨的同震形变场,并利用Okada矩形弹性位错模型反演了震源参数以及同震的滑动分布。监测结果显示,同震破裂引发的最大隆升和最大下沉形变均超过10 cm。同震位错模型反演结果表明：漾濞M_w 6.4级地震是由长度12.2 km、宽度5.0 km、走向132°、倾角83°的断层引发的;此次地震是一次显著的右旋走滑运动的破裂事件,其平均滑动量为0.19 m,平均滑动角为-153.6°,矩震级为M_w 6.1,与美国地质调查局(USGS)、全球矩心矩张量(GCMT)计划等公布的震源机制解结果一致。同时,同震库伦应力变化计算结果显示,同震引发的库伦应力变化为4.21×10²⁵ dyne·cm^-1,相当于矩震级...     

GPS约束下的2022年台湾MW6.9地震破裂滑动分布

利用9个GPS站水平和垂直向形变数据反演2022年台湾M_W6.9地震破裂滑动分布。结果表明,此次地震破裂至地表,滑移以走滑为主兼逆冲分量,主要沿NNE向延伸,有2个破裂集中区,共释放地震矩约5.73×10¹⁹ Nm。基于GPS观测和位错理论模型对此次地震的地表形变进行分析,认为同震形变整体上符合台湾岛东海岸区域构造运动特征,菲律宾海板块俯冲欧亚板块是中央山脉断裂活动的主要动力来源。     

利用GPS位移数据校正2011日本Tohoku地震的InSAR形变场

轨道误差是InSAR数据处理中的一个重要因素,对从最初的图像配准到最后的高程值或形变值图像的生成都有着重要影响。含有误差的轨道参数造成基线误差以残差条纹的形式存在于干涉图中。本文利用4个条带的61景PALSAR升轨数据,通过InSAR干涉测量得到覆盖2011日本Tohoku地震的雷达视线向的地表同震位移场。通过联合GPS同震位移数据,采用二次多项式曲面拟合的方法去除InSAR同震形变数据中的卫星轨道误差,并分析雷达入射角变化对轨道误差去除的影响。对干涉图中拟合残差进行分析,改正后的InSAR同震形变场精度得到较大提高。     

基于光学遥感影像的2019年加利福尼亚MW7.1和MW6.4双震同震形变测量

利用地震前后的Landsat-8和Sentinel-2光学遥感影像数据,基于频率域互相关算法提取加利福尼亚州MW7.1及MW6.4地震的同震形变场,同时针对形变数据易受轨道误差、条带误差及时间失相干的影响,分别采用最小二乘多项式曲线拟合、改进均值相减等方法去除系统误差项。研究表明,MW7.1主震以右旋走滑为主,地表形变的断层迹线呈NNW走向,长度达55 km,最大滑移量约为2.82 m;MW6.4地震的发震断层迹线呈NE走向,长度达15 km,最大滑移量约为1.05 m,推测2次地震的发震断层分别为NNW向右旋走滑断裂及NE向左旋走滑断裂,二者形成典型的共轭关系。     

2020年土耳其MW5.7地震InSAR同震形变特征与破裂滑动分布

为准确认识2020-06-14土耳其MW5.7地震的发震位置、构造特点以及地震危险性,利用D-InSAR技术对Sentinel-1A数据进行处理,基于GACOS进行大气校正,获得视线向（line of sight,LOS）同震形变场。降轨LOS向同震形变场显示,断层北侧抬升,最大抬升形变量约8.87 cm;南侧沉降,最大沉降量约-7.75 cm。以LOS向同震形变为约束,先采用贝叶斯自举优化法反演发震断层几何参数,然后使用有限断层方法反演地震破裂滑动分布。结果显示,断层走向约257.48°±0.65°、倾角约79.69°±0.98°、滑动角约154.2°±3.8°,震中位置为（40.754°E,39.389°N）,破裂区域长度约8 km、宽度约6 km,破裂的最浅埋深约0.8 km、最大埋深约8.9 km,最大滑动量约0.57 m,对应深度约4.278 km。地震释放地震矩约4.54×1017 Nm,对应矩震级MW5.7,与土耳其灾害和紧急情况管理局公布的结果一致。此次土耳其地震受近东西向的潜伏断层控制,以右旋走滑为主兼具少量的逆冲特性。地震造成部分地区的库伦破裂应力增量超过0.1 bar,这些区域未来的地震危险性值得关注。     

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深度处,断层在近地表未出现大面积显著滑移现象,表明此次地震未破裂至地表。     

利用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等机构给出的震级一致。     

荣县地震InSAR同震与时间序列形变初步分析

借助Sentinel-1A卫星升降轨数据，利用合成孔径雷达差分干涉测量(DInSAR)和Okada均匀滑动模型反演获得四川荣县地区三维同震形变场及荣县地震震源参数。结果表明，此次地震造成的断裂现象大体沿NS向，两侧沉降区域中间伴随抬升现象，最大垂直形变量约18 mm；荣县地震发震断层符合逆断层特性，断层长约9.5 km，宽约2.1 km。采用时间序列分析技术获取该地区2018-11～2019-03各时间段的累计形变，结果发现，该区域在震前形变波动较小，地震是造成累计形变较大的主要原因，震后地表变化趋于稳定。     

基于GPS观测值研究日本M_W9.0地震震后形变机制

基于GPS观测的震后水平位移对2011年日本MW9.0地震的震后形变特征进行研究。震后近5a(截至2015-12),震后水平位移累积达到东向60~165cm,南向20~65cm,距离震中较远处的G104、G105及J192站点观测到的震后位移累积变化已超过同震位移,且震后形变还在持续。联合震后余滑和粘弹性位错理论模拟震后形变,利用这2种作用机制对震后GPS水平形变进行解释。研究表明,震后余滑在震后形变最初阶段起主要作用,但随着时间增长而逐渐衰减,粘滞性松弛作用的贡献随着时间增长而变大,GPS观测到的震后地表形变可由这2种机制结合得到较合理的解释。利用震后GPS形变模拟估算地震区域的地幔粘滞系数在1.5×1019 Pa·s量级。     

利用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断层余滑量相对较小,滑移区范围明显减小且均位于浅部区域,北部断层能量释放较彻底,南部仍存有少量能量,整体能量基本释放完。     

九寨沟Ms7.0地震的InSAR观测及断层破裂研究

为提高九寨沟2017-08-08 MS7.0地震的InSAR同震形变场精度,厘清发震断层的构造形态及形成机制,对形变区相干系数进行统计分析并确定相干性阈值,利用GACOS对形变场进行大气校正,再根据余震分布和地质背景确定发震断层的基本形状,最后基于Okada弹性半空间位错模型反演发震断层的滑动分布。改正后的InSAR同震形变场显示,视线向最大下沉为25 cm,最大抬升为10 cm,分别位于震中西北和东南,形变长轴为北西向,主要形变区位于发震断层西部。改正后的InSAR形变场残差均方根较改正前小,最大滑动量为0.9 m,平均滑动角为-0.5°,破裂主要集中在地下1~20 km范围内,矩震级为6.5,与USGS和GCMT等机构的结果一致。研究表明,利用GACOS改正九寨沟InSAR同震形变场对提高形变场精度具有一定的作用,反演断层滑动分布的结果较改正前差别不明显,发震断层的属性与虎牙断裂北段的性质基本一致。结合余震重定位结果可以推断,发震断层为虎牙断裂的北向延伸部分,此次地震事件为巴颜喀拉地块南东向扭转与华南地块碰撞的结果,中下地壳粘性流体的差异分布是导致虎牙断裂倾角变化的主要原因。