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

1969年阳江6.4级地震发震构造研究

综合研究了阳江地区的地震地质构造、深部构造特征以及形变、断层气测量结果，对1969年7月26日阳江6．4级地震的宏观影响场、震源机制解及余震分布特征进行了分析，认为该次地震的发震构造为NEE向的平冈断裂。     

2016年杂多MS6.2地震的InSAR形变及断层滑动模型

利用欧空局Sentinel-1A卫星的升轨观测,基于InSAR技术获取2016-10-17青海杂多地区MS6.2地震的同震形变场,并分别采用最速降法和三角元法反演地震的位错分布模型。结果表明,本次地震发生在杂多上拉秀断裂带,走向北东东,地震破裂以正断倾滑为主,兼具少许左旋走滑运动分量,最大滑移量0.17 m,震源深度为9 km,平均滑动角约为-76°,最大形变中心位于32.85°N、94.86°E附近,地震位错主要分布在深度8~12 km,与余震的深度分布较为一致,反演给出的矩震级为MW5.9,与USGS和CENC测定的矩震级和面波震级基本一致。     

改则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可以作为监测青藏高原大陆动力学的有效技术。     

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公布的震级基本一致。综合分析震中位置和滑动机制认为,甲岗-定结断裂的分支断层为本次谢通门地震的发震断层。     

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

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

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

利用近震波形研究2016年唐山M4.1地震震源参数及发震断层

基于首都圈地震台网观测数据，利用Hypo2000方法对2016-09-10河北唐山M4.1地震主震震中进行准确测定。采用近震波形反演方法，确定主震的最佳双力偶震源机制解为：节面Ⅰ走向122°、倾角60°、滑动角-42°；节面Ⅱ走向236°、倾角54°、滑动角-142°；矩震级约为M_W4.3。结合近震深度震相sPL确定主震的震源深度约为6 km。使用双差定位法对余震进行重新精确定位，结果显示，余震深度集中分布在4~9 km，整体形态近于铅直，表明发震断层面具有倾向北西、倾角较陡的特点，与节面Ⅱ的性质较为一致，据此推测节面Ⅱ为发震断层面，节面Ⅰ为辅助面。将发震断层面参数与震源区附近断层性质进行对比分析认为，发震断层应为唐山断裂带中唐山-丰南断裂。     

2021年云南漾濞地震InSAR监测及断层滑动分布反演

采用双轨D-InSAR技术获取2021-05-21云南漾濞6.1级地震的升降轨同震形变场,并以此为约束,基于Okada均匀弹性半空间位错模型反演发震断层的几何参数与同震滑动分布,最后通过正演模型预测地震的同震三维形变场,分析地震的形变特征及发震断层的几何形态和构造运动。结果显示,在InSAR形变场中,升轨LOS向最大形变量约10 cm,降轨LOS向最大形变量约13 cm;利用Okada模型反演得出,发震断层长约12.53 km,宽约6.44 km,位于地下3～9 km处,断层走向约136.34°,倾角约83.79°,滑动角约-175.46°;断层的同震滑动主要分布在沿走向9～21 km与沿倾向向下3～9 km范围内,最大滑动位移出现在地下4 km深度处,滑移量达到0.64 m;通过正演模型的地表三维形变预测发现,垂直方向最大抬升约3.22 cm,最大沉降约4.04 cm。总体而言,此次地震的发震断层类型为右旋走滑,且带有正断分量,初步判断地震发生在一条未知的NW-SE向断裂带上,该断裂可能属于维西-乔后-巍山断裂和红河断裂北段的次级或分支断裂。     

九寨沟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同震形变场对提高形变场精度具有一定的作用,反演断层滑动分布的结果较改正前差别不明显,发震断层的属性与虎牙断裂北段的性质基本一致。结合余震重定位结果可以推断,发震断层为虎牙断裂的北向延伸部分,此次地震事件为巴颜喀拉地块南东向扭转与华南地块碰撞的结果,中下地壳粘性流体的差异分布是导致虎牙断裂倾角变化的主要原因。     

台湾美浓MW6.4地震InSAR形变场初步分析

利用欧洲航天局最新发射的宽刈幅、高分辨率Sentinel-1A卫星,第一时间获取2016-02-06台湾美浓MW6.4地震干涉像对,使用D-InSAR技术获取美浓地区的同震形变场。利用震中附近39个同震GPS观测进行对比验证后显示,InSAR获取的同震形变场精度优于1 cm（3σ）。形变发生在宏观震中30 km×30 km范围内,主要表现为沿雷达视线向抬升,最大抬升约12 cm。从形变场空间分布特征可以看出,空间连续性较好,说明宏观震中附近地表未发生明显破裂。宏观震中并不与震中位置重合,而是位于震中西部约15 km处。     

2021年阿拉斯加MW8.2地震震源破裂过程反演

基于GPS近场形变和远场P波,利用有限断层方法反演2021-07-29阿拉斯加M_W8.2地震的震源破裂过程。结果显示,此次地震持续时间约100 s,释放的地震矩约为1.59×10²¹ Nm。断层破裂滑移自震中沿断层向NEE延伸,震中东北侧的断层滑动量较大。基于GPS数据和位错理论模型分析地震的地表形变,同震形变整体上符合海沟特大地震的逆冲断层弹性回跳理论模式。受到断层闭锁的影响,太平洋板块俯冲北美板块,使得上盘北美板块受到长时间持续挤压,地震发生瞬间北美板块转变为拉张松弛状态。     

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,造成藏南申扎-定结拉张地堑应变加载。     

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

2020-07-23西藏尼玛MW6.3地震震源参数InSAR反演与发震构造研究

针对2020-07-23西藏尼玛M_W6.3地震发震构造的复杂性(无明显地表破裂、发震断层倾向存在争议等),利用哨兵-1A升、降轨影像数据获取此次地震的同震形变场,并以此为约束通过设置两组可能性断层滑动模型及结合区域余震分布和地质构造背景,探讨此次地震发震断层的活动方式,在此基础上进一步反演获得断层面滑动分布特征。结果表明,此次地震发震断层为依布茶卡地堑东侧NWW向正断层,断层走向为207°,倾角为33.1°,平均滑动角为-89.27°,以正断倾滑为主;地震破裂主要集中在沿断层面5.4~9.3 km深处(未破裂至地表),最大滑动量约1.7 m,矩震级为6.3级。此次地震的发生是青藏地块现今构造运动变形的结果。     

2016年门源MS6.4强震的发震构造及其对“天祝地震空区”的影响

在综合分析区域活动断裂活动特征、震源机制解、余震精定位等资料的基础上,结合GPS观测数据,推断冷龙岭断裂倾向NE,而非部分早期研究推测的倾向SW。考虑到2016年和1986两次门源M_S6.4强震相似的震源机制与冷龙岭断裂的位置关系以及冷龙岭断裂倾向NE等事实,可推测判定这两次强震的发震构造应为冷龙岭断裂;冷龙岭断裂处于祁连-海原活动断裂系内著名的破裂空段“天祝地震空区”的西段,该地震空区内部的金强河断裂可能为一个次级的地震空区,未来存在6级左右强震危险;冷龙岭断裂目前处于比较活跃的状态,说明该段的应力、应变积累水平较高,未来存在大震可能,且不排除冷龙岭断裂与金强河断裂同时破裂的可能性。