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

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年玛多MW7.3地震同震形变提取

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

基于升降轨Sentinel-1 SAR影像研究精河MS6.6地震震源机制

采用Sentinel-1A/B卫星升降轨SAR数据获取了2017-08-09新疆精河MS6.6地震的同震形变场，并以InSAR形变场为约束，利用均匀滑动模型和分布式滑动模型反演了发震断层的滑动分布。结果表明，升降轨InSAR形变场均显示为隆升，形变影响范围约30 km × 40 km；发震断层走向约为83.8°，倾角约为40.6°，滑动角约为89.1°，震源深度约为20.7 km，矩震级为MW6.35，同震位错以逆冲运动为主兼有极少量左旋走滑分量。精河地震的发震断层为库松木楔克山前断裂，该次地震是新疆北天山地区逆冲断裂带深部发生错动的结果。     

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资料反演宁洱Ms6.4地震震源参数及库仑应力变化

首先利用ALOS PALSAR数据，通过D-InSAR技术获取2007-06-03云南宁洱MS6.4地震的同震形变场，然后基于Okada弹性半空间位错模型反演该地震的断层几何以及精细滑动分布，最后计算宁洱地震后周边断层的静态库仑应力变化。结果表明，形变主要集中在西盘，最大视线向形变量为51.6 cm；反演得到的震源位置为23.05°N、101.02°E，深度3 km，断层走向145°，倾向49.5°，平均滑动角153°，发震断层为NNW向普洱断裂，断层活动以右旋走滑为主，兼具逆冲分量；断层面最大滑动量为1.2 m，反演得到的震级为MW619。基于库仑应力场发现，磨黑断裂处于库仑应力增加区域，而2014年景谷地震位于负值区域。结合实地考察资料和反演结果表明，宁洱地震为浅源地震，但断层并未出露地表。     

基于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向地堑得到进一步扩张。     

基于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。结合余震分布认为,该地震的发震断层为隐伏的虎牙断裂北段。     

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

利用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年定结地震发震构造是定结-申扎伸展断裂系中的一条新生盲断层。     

�ٰ�2004��Ms6.7��2005��Ms6.5�����InSAR�۲⼰�α����

???????Envisat??????????????2004??7??12???2005??4??8?????????Ms6.7??Ms6.5?????InSAR????α??????????????ε????????????20 km??20 km??Χ???EW?????????Σ??????η???19.0 cm??30.5 cm???????????????????????????ε???????????????????????,????λ??83.71??E,30.70??N,??Mw6.1,83.72??,E30.52??N????Mw6.2??????????????? 1 m??1.4 m?????ε??????????NW?????3?????????????NS?????????-????????????λ     

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地震序列对该处规模不一的多断层构造组成的复杂断层系造成破坏,小湖断裂带为地震序列的发震断层。     

����Ms6.6�����InSAR�۲⼰�ϲ�λ����

????3??Envisat?????????,?????2008??10??6?????????Ms6.6??????????????α?,????????????4????????????δ?????????????α?,?????????????????20 km??20 km??Χ??????0.3 m????????Ρ??α??????????????????????????????????????t???????????λ???????????????????????????λ?????????????????????????????????????????????λ???3 m, ????λ??90.374??E, 29.745??N????Mw6.35??????????????189?????60????????????????′???????ε?μ?????????????12 km,??11 km,???9.5 km???????????????????????????????????????????????й??     

1923年仁达地震与1973年炉霍地震的破裂模型及触发关系研究

采用1923年仁达MS7.3地震和1973年炉霍MS7.4地震地质考察的同震破裂数据,基于弹性半空间位错模型和最小二乘方法反演历史地震破裂模型。以仁达地震的同震滑动分布为扰动源,采用PSGRN/PSCMP程序计算其在断层面上产生的同震和震后库仑应力变化。结果表明,仁达地震最大滑动量约为3.3 m,矩震级为MW7.0;炉霍地震最大滑动量约为3.7 m,矩震级为MW7.3。2次地震的最大滑动量均位于浅部0~5 km,且震源机制相似,以左旋走滑为主。仁达地震同震及震后效应引起炉霍地震震中区域的库仑应力增强约29.79 kPa,达到地震触发阈值10 kPa。     

台湾美浓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处。     

俯冲带地震循环的数值模拟——以日本Tohoku MW9.0地震为例

基于滑移弱化摩擦准则,以日本2011年Tohoku MW9.0地震为例,建立一个以物理规律控制地震循环过程的二维有限元数值模型。结果显示,参考模型在1 000 a间的6次大地震表现出特征地震的规律,地震重复周期约161 a,单位破裂长度地震矩为1.13×1020 Nm/km,在两次大地震中间会发生一次5.62×1018 Nm/km的小地震。参考模型的数值结果与地表同震GPS位移、震间GPS速度分布具有较好的一致性。模型弹性参数的不确定性对同震以及震间形变的影响有限,模型粘性参数主要影响震间形变场。数值计算也显示,假设震间形变仅由断层运动规律所决定,那么在一个地震周期内,模型空间重力异常基本上随时间均匀变化,在大陆一侧距海沟100 km处可达-370 μGal;速度场的变化主要发生在震后约5 a的时间内,此后基本保持稳定增加。     

�����ݶ��½�����������Mw7.1����ͬ�𻬶��ֲ�

????InSAR??????ENVISAT ASAR?????????????2008??3??21??????7.1??????????α?????????????????????????????????????????????????????????????????????????????????????????????????????????????λ??????0~14 km???????????λ???ж?F2??2.05 km?????5.47 m??????F1????????????????????????????ж?F2?????F3-N/F3??S??????????????????????????????????????????????ж?????Σ????λ???????С?????????????Mo=5.16??10 ¹⁹ N??m?????????7.07????