2017年8月8日四川九寨沟7.0级地震InSAR同震形变场及断层滑动分布反演

基于InSAR技术,利用欧空局升降轨Sentinel-1A/IW宽幅数据,获取了2017年8月8日四川九寨沟7.0级地震InSAR同震形变场,并以升降轨InSAR观测结果为约束,反演了断层滑动分布,基于三种不同接收断层计算了同震库仑应力变化.结果表明,同震形变场发生在塔藏断裂、岷江断裂和虎牙断裂交汇的三角地带,升降轨干涉位移均显示本次地震的形变场影响范围约为50 km×50 km,形变场长轴方向为NW向,升降轨观测的形变量相反,反映断层运动性质以走滑运动为主,升降轨数据观测得到的最大LOS（Line of Sight,视线向）形变量分别为~22 cm和~14 cm.非对称形变场反映出断层两侧的运动差异.反演结果显示,最大滑动量约为1 m,平均滑动角为-9°,矩震级为M_W6.5,地震破裂主要集中在地下1~15 km深度范围内,但整体而言本次地震破裂较为充分,基本将该区域1973年及1976年4次> M_W6.0地震的破裂空区完全破裂.考虑到塔藏断裂和虎牙断裂的运动性质,可初步判定发震断层为虎牙断裂北侧延伸分支.基于三种不同接收断层模型的同震库仑应力变化计算结果反映出该区域以应力释放为主,进一步触发较大走滑型余震的可能性不大.

     

基于GNSS和InSAR约束的2021年玛多MS7.4地震同震滑动分布及应用

2021年5月22日青海省果洛藏族自治州玛多县发生M_S7.4地震,震中位于青藏高原中部的巴颜喀拉块体,这是近20多年来在巴颜喀拉块体周边发生8次M≥7级强震后,块体内部的一次强震,也是汶川地震以来中国大陆发生的最大一次地震,因此该地震的成因及周边地区未来的地震危险性值得重点关注.本文利用震后及时获取的39个近场流动GNSS观测,联合61个GNSS连续观测、Sentinel-1和ALOS-2 InSAR观测获取了本次地震精细的同震形变场,以此为约束,基于均匀弹性半无限位错模型,反演了发震断层的滑动分布,并计算了同震库仑应力变化.GNSS水平同震形变十分显著,断层南北两侧的GNSS点位,最大水平形变分别达0.7 m和-1.2 m,距震中200 km的测点仍有1 cm左右的同震形变.Sentinel-1和ALOS-2的升降轨InSAR同震形变场显示此次地震造成了约160 km长的地表破裂,最大视线向形变分别达0.9 m和1.2 m.同震滑动分布模型显示,发震断层由主段和次段组成,长度分别为170 km和20 km,主段倾向北,倾角85°,平均滑动角为-4.36°,表明玛多地震是一次典型的走滑型地震.次段倾向南,倾角68°,平均滑动角为-11.84°.地震破裂主要集中在0~15 km深度范围,最大滑动量为4.4 m,对应深度6.97 km.反演给出的矩震量为1.61×10²⁰N·m,对应矩震级M_W7.4.主发震断层上存在4个凹凸体,玛多地震是一次不对称双侧破裂事件.结合余震精定位、野外调查及地质资料,我们认为主发震断裂为昆仑山口—江错断裂,东部的次级破裂与主破裂机制不同.同震库仑应力结果显示,东昆仑断裂玛沁段应力有所增加(>0.01 MPa),处于应力加载状态,未来发生强震的危险性较高.     

九寨沟M_s7.0地震的InSAR观测及发震构造分析

2017年8月8日四川省九寨沟县发生M_s7.0地震.本文基于Sentinel-1 SAR影像,利用InSAR技术获取了此次地震的同震形变场,反演获得同震滑动分布,计算了同震位错对余震分布和周边断层的静态库仑应力变化,并对发震构造进行了分析讨论.结果表明:①InSAR同震形变场显示,九寨沟地震造成地表形变最大量级约为20 cm(雷达视线方向),同震形变存在非对称性分布特征.②同震位错以左旋走滑为主,主要发生在4~16 km深度,最大滑动量约为77 cm,位于9 km深处.反演得到的矩震级为Mw6.46.同震错动未破裂到地表.③大部分余震发生在库仑应力增加区.此次地震增加了震中周边地区一些断裂的库仑应力,如东昆仑断裂带东段、龙日坝断裂、虎牙断裂等.④东昆仑断裂东段的未来地震危险性值得关注.⑤九寨沟地震的发震断层为树正断裂,可能是虎牙断裂的北西延伸隐伏部分,此次地震是巴颜喀拉块体南东向运动受到华南块体的强烈阻挡过程中发生的一次典型构造事件.     

2015年新疆皮山MW6.4地震发震断层和滑动分布反演

本文基于Sentinel-1A卫星影像数据提取了2015年皮山M_W6.4地震的同震形变场, 震中北部以隆升为主, 最大抬升量为12.9 cm; 南部以沉降为主, 最大沉降量为5.5 cm。 采用基于单一断层滑动模型的多峰粒子群优化和蒙特卡罗算法, 以LOS向InSAR形变场为约束, 对发震断层的几何模型进行非线性反演。 在此基础上, 联合InSAR和GPS数据, 利用最速下降法反演断层滑动分布。 综合结果表明: 发震断层是顶部埋深约7.4 km的隐伏断裂, 断层面大小为48 km×35 km, 断层走向、 倾角、 断层滑动角分别为111°、 19°、 91°; 断层最大滑动量0.47 m, 位于深度为10.6 km的区域; 累计地震矩3.89×10¹⁸ N·m, 约合矩震级M_W6.33。 最后, 依据主震断层滑移量计算了主震对周围中小断裂的库仑应力扰动变化, 结果显示距离震中最近的泽普断裂受主震影响的库仑应力明显增加; 震后3年内余震集中分布在泽普断裂库仑应力增加区域, 表明皮山地震主震对余震的发生可能具有一定的应力触发作用。     

2021年青海玛多MW7.4地震GNSS同震形变场及其断层滑动分布

2021年5月22日青海玛多M_W7.4地震作为发生在巴颜喀拉块体内部的一次强震,再次引起了人们对该地区地震活动性的强烈关注.本文基于震后GNSS流动观测和区域连续GNSS站资料,解算了106个站点的同震形变及其中17个站点的高频形变波形.同震形变场显示玛多地震具有典型的左旋走滑特征,GNSS观测到的最大同震位移达到1.2 m.GNSS与InSAR数据相符度较高,GNSS提供了准确的近场形变信息.基于GNSS同震形变场,本文反演了断层滑动分布,并计算了发震断层上产生的库仑应力变化.结果表明,发震断层的滑动破裂存在多个凹凸体,破裂分段特征明显且出露地表,与野外地表破裂考察和余震分布吻合,主体破裂位于断层面0~10 km的浅部区域,最大滑动量达到4.6 m,地震矩1.63×10²⁰N·m,矩震级为M_W7.4;发震断层上静态库仑应力增加区域与余震分布具有一致性,说明余震主要是由静态库仑应力加载而触发的.     

2022年青海门源MW6.6地震发震构造——来自InSAR和高分影像约束

2022年1月8日,在青海省门源县海原断裂带冷龙岭断裂西段和托莱山断裂东端发生了M_W6.6地震,发育了明显地震地表破裂带.本文基于欧空局哨兵1号雷达影像,利用InSAR技术获取了门源地震的同震形变场,升降轨InSAR同震视线向位移表现出相反的形变特征,量级达到～60 cm左右,同震形变存在非对称性分布特征,结合高分7号观测解译了近断层地表形变.利用InSAR观测反演获得了发震断层参数及详细滑动分布,计算了同震库仑应力变化,并对发震构造及震中区域未来地震危险性进行了分析讨论.结果表明：门源地震至少有两条断裂发生了破裂,主断层对应冷龙岭断裂西段,InSAR确定的最优断层模型显示主断层东段存在沿走向变化特征,西段则在地质解译断层基础上向西延伸,次断层对应地质解译的托莱山断裂东端,两个断裂组成一个平躺的Y型分布.断层最大滑动量约为～3.7 m,断层浅部存在滑动,表明该地震破裂到了地表,地表破裂长度约19 km.主断层滑动主要集中在0～9 km深度范围;次断层滑动主要集中在0～4 km深度范围,InSAR确定的矩震级为M_W6.6.库仑应力变化结果显示民...     

2022年泸定6.8级地震GNSS同震形变场及其约束反演的破裂滑动分布

2022年9月5日四川省甘孜州泸定县发生6.8级地震,地震发生在川滇菱形块体东边界鲜水河断裂磨西段西侧附近.本文基于震中350 km范围内中国大陆构造环境网络和中国地震科学实验场118个GNSS连续站数据,观测获得了精细的同震形变场,结果显示：同震形变跨发震断层呈空间对称分布,表明本次地震具有显著的左旋走滑特征;记录到的最大同震形变发生在震中距40 km的SYD5(石棉安顺场)站,东西向和南北向形变量分别达到-22.0±1.2 mm和11.6±0.9 mm,震中距100 km以外测站的水平同震形变均小于5 mm;垂向同震形变不显著.结合走向163°和倾角77°的发震断层模型,本文对发震断层面的同震破裂滑动分布进行了反演,结果显示：同震滑动主要集中在主震东南侧余震空区0~15 km深度范围内,且破裂达到了地表,同震释放的地震矩与矩震级为M_W6.57地震相当.结合理论同震形变场、主应变场和邻近区域主要活动断裂库仑应力变化的空间分布特征,本次可能导致处于强闭锁和地震空区的安宁河断裂石棉—冕宁段未来发震风险性增强.

     

基于Sentinel-1A数据的四川长宁MS6.0地震同震形变场分析及断层滑动分布反演

利用大地测量数据研究2019年6月17日四川长宁M_S6.0地震同震形变场特征和发震断层参数, 基于DInSAR技术处理升降轨Sentinel-1A数据获取干涉相位图, 并考虑大气折射效应和余震形变误差实现同震形变场改正。四叉树采样后的形变数据作为反演数据源, 采用弹性半空间分层模型反演发震断层几何面滑动分布。结果表明本次地震发震机制为兼具逆冲和左旋走滑, 矩震级为M_W5.9, 断层破裂尺度达28 km×20 km, 震源深度约9.4 km。升降轨视线向同震形变场在断层两侧呈现形变特征差异, 最大沉降量分别是8.34 cm(升轨)和4.23 cm(降轨), 最大抬升量分别是5.5 cm(升轨)和7.5 cm(降轨); 发震断层走向为302°, 倾角为43°, 平均滑动角为50°, 断层面最大滑动量达到0.28 m。     

利用InSAR技术研究2016年青海门源MW5.9地震同震形变场及断层滑动分布

2016年1月21日, 青海省门源县冷龙岭断裂带附近发生了M_W5.9地震。 基于Sentinel-1A影像, 采用差分干涉雷达测量技术研究了此次地震产生的同震形变场, 结果表明, 门源地震的形变影响范围约20~30 km, 形变态势在升降轨道形变场均显示为隆升, 基本沿冷龙岭断裂呈近似同心圆展布, 推测可能是冷龙岭断裂与民乐—大马营断裂之间的一条逆断层, 沿雷达视线方向最大形变量级约为6 cm。 均匀滑动反演显示门源发震断层长7.3 km, 宽6.2 km, 走向298.6°, 倾角34.5°, 倾向宽度9.5 km, 沿走向滑动量为170 mm, 沿倾向滑动量为460 mm, 矩震级为M_W5.97; 分布式滑动反演显示门源地震以逆冲为主, 兼具少量右旋走滑分量, 滑动量主要集中在沿断层倾向方向, 距离地表5～15 km处, 最大滑动量约0.3 m, 位于断层倾向深度10 km处, 矩震级为M_W5.93。     

基于升降轨Sentinel-1 SAR数据研究2020年新疆于田MS6.4地震震源机制

2020年6月26日新疆于田县发生M_S6.4地震。 本文先利用Sentinel-1A/B卫星升、 降轨SAR数据, 结合InSAR技术提取此次地震的同震形变场; 并以同震形变场为约束, 利用贝叶斯方法的均匀滑动模型反演发震断层的几何参数; 最后基于梯度下降法(Steepest Decent Method, SDM)来确定发震断层更为精细的滑动分布。 研究结果表明: InSAR形变场的分布主要沿北北东方向, 东西方向跨度约40 km, 南北方向跨度约30 km, 形变特征为西侧沉降与东侧隆升, 西侧最大沉降约200 mm, 东侧最大隆升约70 mm。 于田地震发震断层为一条走向187.68°、 倾角59.78°、 滑动角77.76°的隐伏断层; 于田地震的宏观震源深度为6.29 km, 距震级M_W6.19。 断层滑动分布主要集中在沿走向7~21 km与倾向4~11 km的区域内, 平均滑动量约0.2 m, 在沿倾向7~8 km处的最大滑动量约0.97 m, 同震位错主要表现为正断运动。     

GNSS观测的九寨沟7.0级地震同震位移初步结果

2017年8月8日四川省九寨沟县发生了7.0级地震,中国大陆构造环境网络与北斗地基增强系统的GNSS连续观测共同监测到了此次地震的同震位移（坐标：东向为正,北向为正）,结果显示：3个站点记录到了明显的同震位移,距离震中43 km的九寨沟台站（SCJZ）在东西向的位移为-9.8±1.5 mm,在南北向的位移为3.3±0.7 mm;距离震中65 km的松潘站（SCSP）在东西向的位移为-1.8±0.7 mm,在南北向的位移为-7.7±0.6 mm;距离震中77 km的舟曲站（GSZQ）在东西向的位移为0.4±1.2 mm,在南北向的位移为3.6±0.8 mm.通过同震位移分布特征,可以推测此次地震为一次左旋走滑型事件,引起水平向同震位移大致不超过150 km范围,地震对东南侧的龙门山断裂带影响非常小,对北侧的塔藏断裂和西侧的岷江断裂处引起的同震位移为厘米级.同震位移的反演结果显示：断层面上滑动量主要集中在7 km深度,最大量值约为0.4 m,平均滑动角为-15°,利用滑动分布计算的相应矩震级为M_W6.4,与地震波反演结果相当.结合同震滑动分布、同震主应变分布、余震分布和震源机制解等特征,推测此次地震破裂极值区累积的能量得到较充分释放,进一步分析得出此次地震在塔藏断裂、岷江断裂和虎牙断裂处产生了一定的应力变化,值得持续关注.

     

联合强震记录和InSAR/GPS结果的四川九寨沟7.0级地震震源滑动分布反演及其地震学应用

2017年8月8日我国四川九寨沟发生里氏7.0级地震.本研究利用基线校正方法获得距震中100 km范围内9个强震台站同震位移,基于Sentinel-1卫星干涉SAR影像对获取了InSAR同震形变场.结合GPS形变数据,本研究进行了震源滑动模型联合反演,结果显示此次地震整体以走滑运动为主,释放地震矩约为7.60×10¹⁸N·m（~M_W6.52）.通过对比模拟形变场和观测值显示,联合反演结果优于单独基于InSAR形变场的反演结果.静态应力变化计算结果显示断层平均静态应力降为1.07 MPa.反演滑动模型沿走向和倾角方向拐角波数值分别为0.99×10^-4和1.10×10^-4.同震静态库仑应力变化计算结果显示共有83.6%的余震位于库仑应力增加的区域,被主震所触发的余震占总数的77.9%,主震对后续余震具有显著触发作用.强地面运动模拟结果显示模拟结果在烈度分布范围和等级方面与调查烈度符合度很高,模拟结果能够很好地反映断层破裂的方向性效应等特征.本研究计算结果显示九寨沟地震无论是平均静态应力降还是拐角波数均低于同类型地震的平均水平,这可能是造成本次地震强地震动水平相对不高的原因.

     

近震全波形反演2017年九寨沟M7.0地震序列震源机制解

搜集了四川地震台网的波形资料,采用全波形反演2017年8月8日九寨沟M7.0地震序列震源机制解.反演结果显示,九寨沟主震矩震级为M_W6.36,震源深度为22 km,节面Ⅰ走向为150°,倾角为80°,滑动角为-20°;节面Ⅱ走向为244°,倾角为70°,滑动角为-169°.余震主要分布在14~22 km深度范围内,震源机制以走滑型为主,其中正断型地震2个,逆冲型地震2个,走滑型地震24个,混合型地震8个.断层面优势方向为SSE向,与塔藏断裂和虎牙断裂走向基本一致,但与塔藏断裂最南段存在明显差异.倾角变化集中在60°~80°,滑动角主要分布在0°附近,表明九寨沟地震序列主要受SSE走向、近似直立的左旋走滑断层控制.P轴优势方位为SEE向,仰角主要分布在30°以内,与区域应力场基本一致.震源区的机制类型和应力状态均存在空间分段差异.本文推测此次九寨沟M7.0地震序列可能发生在虎牙断裂向北延伸的隐伏断裂上,但不排除地震引起了塔藏断裂南段和虎牙断裂以北隐伏断裂同时破裂的可能.

     

东昆仑断裂带东端的构造转换与2017年九寨沟MS7.0地震孕震机制

2017年四川九寨沟M_S7.0地震是继2008年汶川M_S8.0地震和2013年芦山M_S7.0地震之后,青藏高原东缘在不到十年的时间内发生的第三个震级M_S7.0以上的强震.这次地震发生在东昆仑断裂带东端,作为青藏高原东北缘的一条大型左旋走滑断裂带,东昆仑断裂带与东端其它构造之间的转换关系仍不清楚,因区内地质构造和地形复杂,东昆仑断裂带东端的主要构造仍缺少深入的研究.本文在总结区域地震构造活动特征、历史地震和现代地震基础上,通过东昆仑断裂带东端已有的和最近开展的活动构造定量研究结果,并结合现今GPS变形场资料和2017年九寨沟M_S7.0地震灾害特征分析,发现东昆仑断裂带最东段塔藏断裂上的左旋走滑除了一小部分继续向东传播转移到文县断裂带上外,大部分转化为其南侧的龙日坝断裂带北段、岷江断裂和虎牙断裂上的近东西向地壳缩短,这可能是岷山隆起的构造机制,而2017年九寨沟M_S7.0地震正是左旋走滑的东昆仑断裂带在东端继续向东扩展的结果.

     

Joint inversion of GNSS and teleseismic data for the rupture process of the 2017 <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript>6.5 Jiuzhaigou,China, earthquake

The spatio-temporal slip distribution of the earthquake that occurred on 8 August 2017 in Jiuzhaigou, China, was estimated from the teleseismic body wave and near-field Global Navigation Satellite System (GNSS) data (coseismic displacements and high-rate GPS data) based on a finite fault model. Compared with the inversion results from the teleseismic body waves, the near-field GNSS data can better restrain the rupture area, the maximum slip, the source time function, and the surface rupture. The results show that the maximum slip of the earthquake approaches 1.4 m, the scalar seismic moment is ~ 8.0 × 10¹⁸ N·m (M_w?≈?6.5), and the centroid depth is ~ 15 km. The slip is mainly driven by the left-lateral strike-slip and it is initially inferred that the seismogenic fault occurs in the south branch of the Tazang fault or an undetectable fault, a NW-trending left-lateral strike-slip fault, and belongs to one of the tail structures at the easternmost end of the eastern Kunlun fault zone. The earthquake rupture is mainly concentrated at depths of 5–15 km, which results in the complete rupture of the seismic gap left by the previous four earthquakes with magnitudes >?6.0 in 1973 and 1976. Therefore, the possibility of a strong aftershock on the Huya fault is low. The source duration is ~ 30 s and there are two major ruptures. The main rupture occurs in the first 10 s, 4 s after the earthquake; the second rupture peak arrives in ~ 17 s. In addition, the Coulomb stress study shows that the epicenter of the earthquake is located in the area where the static Coulomb stress change increased because of the 12 May 2017 M_w7.9 Wenchuan, China, earthquake. Therefore, the Wenchuan earthquake promoted the occurrence of the 8 August 2017 Jiuzhaigou earthquake.