基于InSAR分析2015年尼泊尔MW7.8地震形变场特征及断层滑动分布反演

2015年4月25日尼泊尔发生了M_W7.8地震, 本文基于震前、 震后两景Sentinel-1A雷达影像, 采用D-InSAR两轨差分干涉法提取了此次地震的同震形变场。 结果显示, 同震形变场位于喜马拉雅造山带—主边界逆冲断裂(MBT)和主前锋逆冲断裂(MFT)附近, 形变场整体表现为自西北向往东南方向延伸近150 km的纺锤形包络状, 以大面积隆起抬升形变为主, 视线向最大隆升形变达1.18 m, 抬升区北侧存在一小沉陷区, 以InSAR观测值定位同震最大形变中心。 基于均匀介质弹性半空间模型(Okada模型)与InSAR观测数据反演了断层滑动分布。 反演结果表明该地震属于典型逆冲型地震, 发震断层为主喜马拉雅逆冲断裂(MHT), 同震破裂从主喜马拉雅逆冲断裂(MHT)向上沿着主前锋逆冲断裂(MFT)传递。 基于InSAR同震形变场局部形变细节, 结合震区地质背景、 断裂分布及断层运动特征, 获得了同震破裂拟出露地表迹线。     

GPS和InSAR同震形变约束的尼泊尔M_W7.9和M_W7.3地震破裂滑动分布

2015年4月25日尼泊尔爆发MW7.9地震,继而引发5月12日MW7.3级余震,GPS、InSAR监测到震源区及周边大范围同震形变.本文以国内外的GPS和InSAR同震形变为约束,考虑喜马拉雅断裂带岩石圈垂向分层和横向差异的影响,反演主喜马拉雅逆冲断裂在这次主震和余震中破裂面形状和滑动分布.结果显示,主震从USGS确定的震中位置向东偏南延伸100km以上,破裂地面迹线与主前缘逆冲断裂迹线基本一致.破裂面倾角约7°~11°,大部分破裂集中在深度8~20km,同余震分布深度一致.主震最大滑动量约6.0~6.6m,位于14km深处.余震破裂集中在震中附近30km范围内,填补了主震东部破裂空区,最大滑动约3.6~4.6 m,位于13km深.深度20km以下基本没有破裂.地壳介质不均匀性对破裂滑动分布的影响较大,介质不均匀模型的观测值不符值比各向同性弹性半空间模型降低10%以上.本文地震破裂模型特征与地震反射剖面、以及根据震间期大地测量数据反演的喜马拉雅深部蠕滑剖面极其相似.跨喜马拉雅断裂剖面的震间形变量与地震破裂滑移量直接相关.以此推算,尼泊尔中部大震原地复发周期在300年以上.     

基于InSAR同震形变场反演汶川Mw7.9地震断层滑动分布

通过综合分析2008年5月12日汶川地震野外地震地质考察的地表破裂带空间分布及分段资料,结合InSAR干涉形变场资料,构建了五段断层几何结构模型,该模型与野外地震地质考察结果在多数分段上基本一致;基于此五段断层模型,运用敏感性迭代拟合算法反演了汶川地震InSAR同震形变场,获得了断层滑动分布及部分震源参数.结果表明,基于余震精定位获得的地震断层倾角模型模拟的同震形变场与InSAR形变场吻合较好,且残差较小;反演的滑动分布主要集中于地下0～20 km,最大滑动量分别位于北川及青川等地区,最大可达到10 m;沿SW-NE走向,断层面的滑动方向主要以右旋兼逆冲形式为主,在汶川及都江堰地区以强烈的逆冲为主兼有一定右旋走滑分量,在北川及映秀地区以逆冲兼右旋运动为主,在平武及青川等地区则逐渐过渡为以右旋运动为主兼有一定的逆冲分量,其中汶川地区的平均滑动角为97°,北川地区的平均滑动角为119°,青川地区平均滑动角为138°.反演矩张量为7.7×10²⁰ N·m,矩震级达M_w7.9.     

基于InSAR同震形变场反演汶川M_w7.9地震断层滑动分布

通过综合分析2008年5月12日汶川地震野外地震地质考察的地表破裂带空间分布及分段资料,结合InSAR干涉形变场资料,构建了五段断层儿何结构模型,该模型与野外地震地质考察结果在多数分段上基本一致;基于此五段断层模型,运用敏感性迭代拟合算法反演了汶川地震InSAR同震形变场,获得了断层滑动分布及部分震源参数.结果表明,基于余震精定他获得的地震断层倾角模型模拟的同震形变场与InSAR形变场吻合较好,且残差较小;反演的滑动分布主要集中于地下0～20 km,最大滑动量分别位于北川及青川等地区,最大可达到10 m;沿SW-NE走向,断层面的滑动方向主要以右旋兼逆冲形式为主,在汶川及都江堰地区以强烈的逆冲为主兼有一定右旋走滑分量,在北川及映秀地区以逆冲兼右旋运动为主,在平武及青川等地区则逐渐过渡为以右旋运动为主兼有一定的逆冲分量,其中汶川地区的平均滑动角为97°,北川地区的平均滑动角为119°,青川地区平均滑动角为138°.反演矩张量为7.7×10~(20)N·m,矩震级达M_w7.9.     

2015年尼泊尔Mw7.8地震震源机制InSAR反演及强地面运动模拟

2015年4月25日,在尼泊尔中部发生了M_w7.8地震.本文利用ALOS-2和SENTINEL-1A宽幅数据获取了该地震大范围的同震形变场,并反演了该地震断层破裂的几何特征及运动机制,继而以此为约束资料反演地震强地面运动.InSAR结果显示本次地震造成了巨大的地表形变,LOS向最大抬升量达到1.3 m,最大下沉量达到0.7 m.震源机制反演得到的最优的滑动分布模型表明,断层的走向为291°,倾角为7.6°,倾滑主要分布在深度为12~18 km范围,主倾滑分布范围在长度上达到了140 km,该范围内的平均倾滑角为95°.本次地震最大倾滑量达到5.3 m,位于深度15 km处.累计释放地震矩达 6.5×10²⁰N·m,约合矩震级M_w7.8.该地震发生在印度与欧亚板块俯冲逆冲界面之间,发震构造推断为主喜马拉雅逆冲断裂,属于典型的喜马拉雅型——低角度逆断层型强震.以该滑动分布模型参数为基础利用随机振动的有限断层模型进行尼泊尔地震的强地面运动模拟,结果显示最大地震烈度为Ⅸ度,烈度分布的范围及烈度等级与USGS模型结果对比具有很高的符合度.     

芦山MS7.0级地震InSAR形变观测及震源参数反演

2013年4月20四川省芦山县发生M_S7.0级地震,目前的研究资料表明地震发生在龙门山断裂南段,但地表未发现明显破裂.本研究利用InSAR技术与Radarsat-2雷达数据,获取了芦山地震同震的部分形变场,结果表明,近场区域的LOS位移发生视线向隆升,量值在7 cm左右.随后利用弹性半空间的位错模型反演了断层面参数,综合反演结果及震源机制解最终确定了发震断层的初始模型,以形变场观测数据为约束,基于梯度下降法反演获得了断层面上的滑动分布,反演得到的矩震级为M_w6.45级,断层走向213°,倾角39°~43°,最大滑动位于地表以下约13 km深度位置,最大滑动量0.91 m,平均滑动角71°,整体上仍以逆冲滑动为主,兼具左旋走滑.推测在双石-大川断裂以东12 km处展布一条隐伏断裂,为本次的发震断裂.     

2020年1月19日新疆伽师县MS6.4级地震InSAR同震形变场与断层滑动分布反演

利用欧空局Sentinel-1A SAR数据,重建了2020年1月19日新疆伽师县MS6.4地震InSAR同震形变场.以升降轨InSAR形变场联合约束,反演了发震断层参数与同震滑动分布.结果 表明,地震同震形变场(长轴走向近EW向)发生在柯坪塔格褶皱带与奥兹格尔它乌褶皱带之间的区域内,升降轨观测的视线向(LOS)形变量...     

2015年尼泊尔地震三维发震构造及地震危险性研究

2015年4月25日尼泊尔M_W7.8地震发生在喜马拉雅造山带上,人们普遍认为该地震不足以释放该造山带上累积的能量,但对该地区后续地震危险性评价多基于二维或是假三维的形变反演计算结果.本研究从2015年尼泊尔地震主震与其最大余震M_W7.3地震之间的关系出发,着重分析讨论了尼泊尔地震的时间和空间的非均匀性,结合震源机制解、地壳速度结构、精定位后的余震分布及InSAR反演结果,建立了三维发震构造模型,利用非线性摩擦有限元方法,对一个地震周期内断层摩擦行为和块体变形进行了模拟,将计算结果和地表同震形变、形变反演的同震破裂、历史地震时空演化进行对比,在确认该三维模型可靠性的基础上,讨论了该地区后续地震的可能位置,认为在1934年比哈—尼泊尔M_W~8.1地震以东区域,还存在发生大地震的可能,在最大余震M_W7.3地震东南部位,还存在发生中大地震的可能.     

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地震的破裂空区完全破裂.考虑到塔藏断裂和虎牙断裂的运动性质,可初步判定发震断层为虎牙断裂北侧延伸分支.基于三种不同接收断层模型的同震库仑应力变化计算结果反映出该区域以应力释放为主,进一步触发较大走滑型余震的可能性不大.     

基于多期DInSAR数据的2016年定结地震序列发震断层特征研究

2016年5月22日,在西藏定结县发生四次MW4~5地震,研究本序列地震的发震断层几何和运动特征对于认识周边活动断裂性质具有重要意义.由于发震地区偏远,且观测台网分布稀疏,本研究采用星载雷达干涉测量(DInSAR)技术进行了同震形变场重建,但是定结地震震级较小,单干涉像对获取的形变场受相位噪声影响较大.为了解决这一问题,本研究基于时间序列Sentinel-1A干涉数据生成多期同震与非同震干涉图,并利用叠加平均法对本次定结地震同震形变场进行重建,提取了定结2016年5月22日多次地震产生的同震累计整体形变场.基于InSAR同震形变场和区域地质特征,研究进行了滑动分布反演,确定其主要贡献的发震断层几何参数及滑动分布:断层走向为188°,倾角为43°,平均滑动角为78°,发震断层的运动性质以正断为主兼具少量左旋走滑分量,滑动主要集中在断层垂直深度0~9km处,最大滑动量约为25cm,位于断层倾向深度3km处,反演得到的矩震级为MW5.58.本研究结果表明采用星载InSAR叠加平均技术可以较好地压制相位噪声,有效提取此类中小型浅源地震同震微弱形变场.最后,我们认为本次定结地震与藏南拆离断层与申扎-定结断层的活动密切相关.     

2021年青海玛多7 4级地震的同震变形分析

2021年5月22日,青海省玛多县发生了 7.4级地震,该地震发生在巴颜喀拉地块北部边界东昆仑断裂带以南约70 km,属于块体内部断裂带地震.根据中国大陆构造环境监测网络提供资料,距离震源30多公里的玛多台站记录到东西向永久位移约25 cm.同时,InSAR也观测到明显的形变场,升轨和降轨的最大相对形变量分别约1.87 m和2.32 m.为了解释这些大地测量观测数据,本文利用该地震的三个断层滑动模型,基于不同地球模型的地震位错理论,计算同震变形场,并分别与GNSS观测数据和InSAR视线向形变量对比分析,结果显示基于InSAR数据反演的断层滑动模型产生的位移场与球形地球模型的理论计算结果最为吻合.进一步,利用较优断层模型计算2021年青海玛多7.4级地震的理论同震位移、大地水准面、重力和应变等变化,该结果为玛多地震的GNSS和重力观测的解释提供理论参考依据.     

九寨沟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.同震错动未破裂到地表.③大部分余震发生在库仑应力增加区.此次地震增加了震中周边地区一些断裂的库仑应力,如东昆仑断裂带东段、龙日坝断裂、虎牙断裂等.④东昆仑断裂东段的未来地震危险性值得关注.⑤九寨沟地震的发震断层为树正断裂,可能是虎牙断裂的北西延伸隐伏部分,此次地震是巴颜喀拉块体南东向运动受到华南块体的强烈阻挡过程中发生的一次典型构造事件.     

联合GNSS和InSAR观测位移反演2008年汶川大地震断层位错模型参数

利用现代空间大地测量技术,尤其是卫星合成孔径雷达干涉测量,能够获取高精度、高空间分辨率的同震和孕震形变,为地震断层形变和破裂机制研究提供了前所未有的机遇。本文介绍了利用大地测量观测数据反演地震断层位错模型参数的贝叶斯反演方法。联合运用2008汶川大地震前后GNSS和InSAR技术观测获得的同震位移,反演了地震断层的几何参数和滑动位错分布。研究结果表明,汶川地震的断层滑动主要集中在倾角较陡的浅部,同时包含逆冲和右旋走滑,其中最大逆冲6.1m,最大右旋6.5m。根据断层滑动分布正演计算得到的上盘同震位移明显小于下盘,预示该断层两侧孕震形变可能存在较大的不对称性。     

基于InSAR同震形变观测反演2010年新西兰南岛Mw7.1 Darfield地震同震破裂分布

2010年9月4日新西兰南岛Canterbury平原发生了M_w7.1地震,震源深度约为10 km.本次地震发生在一条震前不为人所知的断层上.我们利用覆盖整个震区的合成孔径雷达(SAR)观测资料,通过干涉处理分析获得雷达视线向(LOS)同震形变场;以此资料为约束反演了断层的几何参数以及同震破裂分布.结果显示,该地震造成四条相对独立断层的破裂.大部分的地震矩释放发生在Greendale断层(编号1-4),其错动以右旋走滑为主,最大破裂约为8.5 m.其它三条断层中,经过震源的逆冲断层最大破裂为5.1 m (编号6),位于Greendale断层以西的逆冲断层最大破裂为3.5 m (编号5),位于Greendale断层北面的走滑断层最大破裂为1.9 m(编号7).反演的Greendale断层地表滑动与地质调查得到的地表破裂在形态和数值上均吻合较好.本次地震释放的地震矩为5.0×10¹⁹N·m,矩震级为7.1.板块边界带形变场分析表明,Darfield地震的发生受边界带应变分配在该地区残留构造应力场控制,其复杂性体现了区域构造应力场的特点.地震对其周围地区的应力场影响较大,库仑应力增加区与余震分布有一定对应关系,并在2011年Christchurch 6.3级地震发震断层区域造成约0.1bar的库仑应力增加,对此地震有一定的触发作用.     

综合InSAR和应变张量估计2016年MW7.0熊本地震同震三维形变场

利用合成孔径雷达(Synthetic Aperture Radar,SAR)影像提取与地质活动相关的三维地表形变场,对深入理解地质灾害的形成机制及其潜在灾害风险评估非常重要.目前,利用SAR影像的同震三维形变场提取主要利用单个像素点的多次观测构建观测方程,然后基于加权最小二乘(Weighted Least Squares,WLS)方法分解从而获得同震三维形变场,因此该方法缺乏对相邻像素点空间相关性的约束.考虑相邻同震位移点的应力连续性,研究学者提出了顾及大地测量应变张量和卫星形变观测的SAR同震三维形变场方法(Extended Simultaneous and Integrated Strain Tensor Estimation from geodetic and satellite deformation Measurements,ESISTEM).本文以2016年MW7.0熊本地震为例,收集了覆盖此次地震的ALOS-2卫星升降轨影像,利用传统差分InSAR(DInSAR)方法和子孔径雷达干涉测量(Multiple Aperture InSAR,MAI)方法分别对升降轨SAR影像对进行处理,得到视线向(LOS)形变和方位向形变,最后利用ESISTEM方法获取此次地震的三维同震形变场.此外,利用GPS和野外考察观测对本文的三维形变场结果进行结果精度分析.研究结果表明,与传统WLS方法相比,ESISTEM方法不仅能有效抑制奇异像素点对形变结果的干扰,同时对近断层的失相干信号能进行较好的恢复,更有助于解释地表破裂区的地震形变特征和掌握地震发生机制.本文确定的三维同震形变场结果显示主形变区发生在Futagawa断层中部和Hinagu断层最北端,最大水平位移为2m,抬升为0.55m.断层破裂以NE-SW走向的右旋走滑为主兼有部分正断成分.应变张量分析表明发震断层处受到了明显的收缩力和剪切力的作用.     

GPS观测结果反映的尼泊尔M_w7.8地震孕震特征

针对2015年4月25日尼泊尔Mw7.8地震的孕震特征,本文首先对覆盖尼泊尔及周边地区的5套GPS水平速度场结果进行了融合,得到了近似统一参考框架下的速度场结果;在此基础上通过对此次地震震源区及周边地区的速度场、应变率场、基线时间序列分析,识别了震前变形特征.GPS应变率场结果显示,喜马拉雅主边界断裂存在大范围挤压应变积累,震源区处于近南北向应变积累高值过渡区.跨喜马拉雅构造带的GPS基线时间序列结果表现为持续缩短现象,表明印度板块与欧亚板块之间的持续挤压变形特征,2012年以来的缩短增强现象反映了印度板块对青藏块体的推挤增强作用明显.距离震中较近的西藏南部GPS同震位移结果以南向运动为主且指向震中,反映了青藏高原存在逆冲应变释放现象.综合此次尼泊尔地震前变形和同震应变释放特征,认为此次地震的孕震区域和同震应变释放区域均较大,将会对青藏高原的地壳变形与强震孕育产生深远影响.     

尼泊尔8.1级地震后河北省数字形变观测同震响应研究

北京时间2015年4月25日14时11分26.3秒在尼泊尔境内发生Ms 8.1地震。通过分析河北省各形变台站的震时观测资料,研究各形变测项对尼泊尔地震的同震响应特征,结果表明,形变测项均记录到此次地震,同震响应形态分为振荡型、振荡-阶变型、阶变型三类;不同台站相同类型仪器的同震响应幅度有很大差别;同一台站同一仪器不同方向上记录波形幅度相差较大;地震对应变场的影响以EW向为主;记录到阶变型形变波的台站主要集中在河北省的北部,南部较少;认为尼泊尔地震后河北省北部地区的地震活动可能有增强的趋势,南部相对平静。     

汶川M_S8.0地震孕震机理的形变证据与模型推演

该文展示了汶川Ms8.0地震区域南部多年流动水准测量给出的形变隆起图像、震前川滇地区GPS给出的区域尤旋转基准下的两期水平形变图像、发震断裂带上耿达跨断层短水准震前变化以及GPS给出的震源区同震水平形变和垂直形变图像.经综合推演,给出了一个板内逆冲型强烈地震孕育发生的模型,简称为"压升冲降"模型,可使上述各异常变化的形变图像得到比较统一和一致的解释.     

漾濞MS6.4地震同震形变特征及发震构造探讨

2021年5月21日云南省漾濞县发生MS6.4地震.文中基于升、降轨Sentinel-1 SAR影像,利用InSAR技术获取了此次地震的同震形变场,反演获取了发震断层的精细滑动分布,计算了区域应变分配及同震位错引起的周边各断裂上的库仑应力变化,对发震构造及周边断裂的地震危险性进行了讨论.结果表明:InSAR同震形变场显...     

用连续GPS与远震体波联合反演2015年尼泊尔中部M_S8.1地震破裂过程

由于印度-欧亚板块碰撞,位于板块边界带的喜马拉雅地区大震频繁,但对其活动性的认识仍十分有限.2015年4月25日尼泊尔中东部地区时隔80年再次发生8级地震,为研究板缘地震提供了一次难得机遇.本文用西藏和尼泊尔的GPS连续观测数据和全球分布的远震地震波记录联合反演此次特大地震的破裂过程,结果显示此次地震发生在印度板块与青藏高原接触边界面——喜马拉雅主滑脱断层上.北倾11°、近东西(295°)走向的断层面破裂约100km长(博卡拉到加德满都),130km宽(从加德满都深入我国西藏吉隆县),破裂以逆冲滑动为主,平均幅度达到2.4m,释放的地震矩高达9.4×1020 N·m.反演结果还显示,震源体主要破裂分布深度范围为5~25km,应无地表破裂,属于一次盲地震.基于GPS资料推测的地壳现今运动速率及1833年地震的震源位置,我们推测地震在此次地震破裂区域复发的周期可能为150~200a,而极震区以南的深部滑脱断层仍保持闭锁,未来仍有导致灾害性大震的可能性.