2015年尼泊尔地震同震滑动及震后余滑的三角位错模型反演

2015年尼泊尔Mw 7.9级地震发生在印度板块向欧亚板块低角度俯冲的喜马拉雅断裂带上。对该地震的滑动模型进行精化，对于理解喜马拉雅断裂带的孕震模式具有重要意义。采用三角形位错元构建主喜马拉雅断裂“双断坡”几何模型，联合全球定位系统（Global Positioning System，GPS）和合成孔径雷达干涉测量（interferometric synthetic aperture radar，InSAR）资料反演2015年尼泊尔地震同震滑移及震后余滑。结果表明，尼泊尔地震最大同震滑移达到7.8 m，深度为15 km，位于中地壳断坡和浅层断坪的接触部位。不考虑中地壳断坡结构会使反演的最大滑移量偏低。震后余滑主要分布在同震破裂区北侧，释放的地震矩为1.02×1020 N·m，相当于一次Mw 7.3级地震，约占主震释放地震矩的12%。同震库伦应力变化和震间断层闭锁分布均表明，尼泊尔地震破裂区南部宽约60 km的区域仍具有较高的地震危险性。     

门源Mw5.9级地震形变InSAR观测及区域断裂带深部几何形态

利用Sentinel-1A卫星升轨、降轨合成孔径雷达影像数据，提取了2016年门源Mw5.9级地震的高精度合成孔径雷达干涉同震形变场，利用单纯形法和非负最小二乘法反演确定了地震断层几何和滑动分布，并构建了区域断裂带的深部几何形态模型。结果表明，门源Mw5.9级地震同震形变以地表抬升为主，沿升轨、降轨视线向的最大值分别为5.3 cm、7.1 cm；地震断层走向、倾角分别为133°、43°；地震滑动以逆冲为主，主要发生在地下6.14~12.28 km处，最大滑动量约0.5 m，平均滑动角为66.85°，地震矩为1.0×1018 N·m（Mw5.94）；形变观测拟合残差均方根为0.36 cm；区域断裂带的深部几何形态以花状构造为特征，整体倾向南西，门源地震发震断裂为花状构造中未出露地表的盲断层。相关成果能够为研究区域地壳运动与变形、活动断裂与地震孕育发生等提供参考。     

日本2011 Tohoku-Oki Mw 9.0级地震的同震形变及其滑动分布反演:GPS和InSAR约束

利用覆盖整个日本东北部的3个条带的Envisat/ASAR降轨数据和6个条带的ALOS/PALSAR升轨数据,通过二通差分干涉处理,首先得到了2011年日本Tohoku-Oki Mw 9.0级地震的初步同震形变干涉结果,然后,利用GPS同震观测值对InSAR结果进行校正和基准统一。在此基础上,采用弹性半空间矩形位错模型,联合GPS和InSAR观测结果对同震滑动分布进行反演,获取了发震断层的断层滑动分布。结果显示,此次地震的滑动分布主要发生在40～50km深度范围内,最大滑移量为50.3m,释放的能量为3.20×1022 N·m(相当于Mw 8.94级)。     

GPS解算的日本Mw9.0级地震的远场同震地表位移

利用日本Mw9.0地震前后各3d的77个GPS连续站观测资料和JPL精密星历,计算了该地震对周缘各国及相邻板块运动的同震影响。结果显示,日本地震对日本、中国东北部、朝鲜半岛均有较大影响,其同震水平位移影响范围达2 000km,地震造成中国东北东移1～3cm,华北东移3～8mm,朝鲜半岛东移2cm。中国境内的GPS测站的垂直同震位移均呈上升趋势,但是上升幅度非常小。结合日本岛测定的GPS同震位移和约20m高的海啸可以推断,本次日本地震的同震垂直位移具有波状起伏特征,符合逆冲挤压型地震的粘弹性理论模型。鉴于近年地震活动趋势,应该密切关注横贯中国东部的深大断裂带（郯庐断裂带）的断层活动,借助中国陆态网络密集覆盖提供的有利条件,跟踪研究郯庐断裂带附近300km范围的区域地壳现今变形、应力应变场的时变过程等,分析郯庐断裂带的断层活动,捕捉异常信息。     

2016年意大利阿马特里切Mw 6.2地震震源机制InSAR反演

2016年8月24日，意大利中部阿马特里切（Amatrice）地区发生Mw 6.2地震。采用ALOS-2条带模式和SENTINEL-1A宽幅模式的合成孔径雷达（synthetic aperture radar，SAR）数据分别进行SAR差分干涉测量处理，获取了该地震的同震形变场。结果显示，本次地震造成意大利中部地区发生明显的地壳形变，在雷达视线向最大沉降量达19.6 cm。基于合成孔径雷达干涉测量（interferometry synthetic aperture radar，InSAR）和GPS同震形变场数据对此次地震的发震断层进行联合反演，通过改进倾角和平滑系数获取方法，得到了最优滑动分布模型。通过使用单断层模型和双断层模型进行反演可知，双断层模型反演结果优于单断层反演结果，两种模型下反演模型相关系数分别为0.85和0.89，发震断层走向分别为160°和158°，倾角分别为44°和46°，倾滑分布主要位于地下5~7 km，平均倾滑角为-80°，最大倾滑量0.9 m位于地壳深度5 km处，该发震断层是亚平宁冲断带的一部分，为NW-SE向延伸的正断层，断层长约20 km。综合使用地震同震形变场和GPS数据对震源机制进行反演、模拟和分析，获取了高精度的震源参数，可以为分析地震危险性和断层破裂参数等提供数据支持。     

利用InSAR研究西藏拉孜地震同震形变

采用欧空局ERS—1卫星雷达数据获取了1993年Mw 6.2级西藏拉孜地震同震形变场,并采用弹性位错模型反演了该地震的断层参数。结果显示该地震产生的最大地表位移约为12 cm,断层倾角为42°,不支持该断层为低角度正断层的论断。     

升降轨InSAR数据约束下的2007年阿里地震反演分析

2007年5月5日,西藏阿里地区发生了Mw6.1级地震。本文采用Envisat卫星的升、降轨SAR数据获取地震的精确同震地表位移,然后基于弹性半空间位错模型,分别采用MPSO非线性和最小二乘线性反演算法确定了断层几何参数和滑动分布模型。结果表明,分布式滑动分布模型能较好地解释观测到的InSAR地表形变场。地震发震断层为走向158°、倾角43°的西南倾断层,主要的滑动量集中在7~12km的深度,最大滑动量约0.3m,位于9km的深度。反演给出的地震矩为1.24×1018 Nm,与地震学结果相一致。     

InSAR数据约束下2016年和2022年青海门源地震震源参数及其滑动分布

中国青海省门源县于2016年和2022年分别发生了Mw 5.9和Mw 6.7地震,相距不足40 km。利用欧洲空间局Sentinel-1A升降轨雷达影像,采用合成孔径雷达干涉测量(interferometric synthetic aperture radar, InSAR)技术分别获取两次地震的同震地表形变场,进而利用弹性半空间的位错模型确定上述事件的震源参数,基于分布式滑动模型反演确定两次地震断层面上的滑动分布,并探讨2016年门源地震对2022年门源地震的发震影响及触发机制。结果表明,2016年门源地震为逆冲型地震,并未破裂到地表,升、降轨同震形变场沿视线向的最大形变量分别为6.7 cm和7.0 cm,断层的最大滑动量为0.53 m,主要集中在地下4~12 km区域滑动。2022年门源地震同震形变场沿NWW-SEE向破裂,降轨影像最大视线向地表形变量为78 cm,断层的最大滑动值达到3.5 m,处于地下4 km左右,断层滑动分布模型揭示此次地震为左旋走滑型地震;结合冷龙岭断裂的运动性质和几何特征,可初步判定发震断层主要为冷龙岭断裂的西段、且极有可能破裂到了其西北端西侧的托莱山断裂。静态库仑应力触发关系显示,2016年门源地震对2022年门源地震的发生有一定的促进作用。     

2022年泸定Mw 6.6地震InSAR同震形变与滑动分布

2022-09-05,青藏高原东缘的鲜水河断裂上发生了泸定Mw 6.6地震,该地震是鲜水河断裂上40年来发生的最大地震,研究该地震的运动学和同震破裂模式对理解青藏高原东缘构造形变机制和评估鲜水河断裂以及安宁河断裂的地震危险性具有重要意义。利用Sentinel-1和ALOS-2卫星雷达影像,采用合成孔径雷达干涉技术获取了泸定地震的同震形变场,进而基于弹性半空间的位错模型,确定了本次地震发震断层的几何参数和滑动分布。结果表明,泸定地震是一次典型的左旋走滑事件,发震断层西倾,倾角约为72°,走向沿NNW-SSE方向,约为167°;断层破裂主要集中在0~10 km深度,最大滑动发生在约5.8 km深度,约为2.23 m;同震释放的地震矩约为8.74×1018 N·m,相当于矩震级Mw 6.59。通过对震后光学影像解译,发现此次地震诱发的滑坡多集中分布在发震断层西侧,该现象与余震主要集中在断层西侧的结果相一致,可认为是地震上盘效应的体现。     

InSAR数据约束下的2023年赫拉特地震序列发震断层探讨及其建筑物损毁评估

2023年10月7日阿富汗西部赫拉特省在不足1 h内接连发生4次Mw 5.5+的地震,称之为“2023年赫拉特地震序列”,此次地震序列是阿富汗境内过去20多年来遭遇伤亡最严重的地震事件,研究该地震序列的发震断层几何和快速分析损坏建筑物的分布状况对理解Herat断裂系统的构造机制、保障高效救援以及灾后规划重建等工作具有重要科学意义。基于欧洲空间局Sentinel-1升降轨合成孔径雷达影像,利用InSAR技术（interferometric synthetic aperture radar）获取2023年赫拉特地震序列的同震形变场,以升降轨InSAR观测为约束,反演确定发震断层几何和断层滑动分布,并对发震构造进行分析;基于多时相InSAR相干性变化探测方法分析并提取本次震后建筑物损毁代理图（building damage proxy map, BDPM）。结果表明,升降轨同震形变场均位于Herat断裂带和Siakhubulak断裂带之间,升降轨数据观测得到的最大视线向形变量分别约为32.3 cm和58.9 cm。反演结果显示,同震位错以逆冲运动为主兼具少量右旋走滑运动,发震断层北倾,倾角约...     

喜马拉雅主逆冲断层闭锁程度与滑动亏损特征研究

基于融合的GPS速度场结果,使用DEFNODE负位错反演程序估算了喜马拉雅主逆冲断层（the main Himalayan thrust,MHT）的闭锁程度和滑动亏损空间分布,并结合剖面结果分析了断层远、近场的运动特征。结果表明,MHT的闭锁深度基本达到18~24 km,断层面闭锁宽度达到102~136 km,两次历史大地震破裂区域之间的未破裂段落和未发生大地震的段落闭锁深度更深,闭锁断层面更宽,2015年尼泊尔Mw7.8大地震就发生在两次大地震破裂区域之间的段落;MHT总滑动亏损速率和垂直断层挤压滑动亏损速率自东向西逐渐减小,平行断层右旋滑动亏损速率则基本上自东向西逐渐增加;MHT 3条剖面拟合结果也反映出其存在很强的闭锁。根据估算的此次Mw7.8地震的复发周期230年和最近500多年发生的大地震分布,认为MHT整条段落尤其是尼泊尔西部与印度接壤处和可能还没有破裂的不丹地区依然有发生8级大地震的危险。     

汶川地震GPS形变约束的破裂分段特征及滑移

为了解2008年汶川地震破裂分段特征及滑移,采用弹性位错模型和模拟退火算法,数值模拟汶川震区密集的GPS同震形变。结果表明,GPS同震形变场至少需要用铲状的映秀-北川五断裂加上灌县－江油断裂来模拟,该模型对GPS数据的符合程度与汶川地震滑动分布模型相当,基本反映了汶川地震破裂特征。映秀－北川断裂总长255~294km。南段以逆冲为主,分上下两层,最深达30km。中段右旋走滑和逆冲都很大,而北段以右旋走滑为主。映秀段地震周期最短为3000年。破裂上层深度14km可能是龙门山中央断裂的闭锁深度。汶川地震可能源于相邻块体的相对运动和挤压、深部滑脱或浅部闭锁,在薄弱构造处首先爆发地震。     

唐山7.8级地震三维地表形变和断层运动特征

通过系统收集整理1976年唐山7.8级地震前后区域水准和边角网观测资料,在研究发展动态平差程序的基础上,解算得到了唐山地震震前和同震三维形变定量结果与误差分布,并反演得到了断层面同震滑动分布。结果显示,震前阶段丰润-唐山-唐海-昌黎-迁安一带整体隆升,隆升高值区位于唐山附近,震中区的水平形变量不大,唐山菱形地块相对于周围地块逆时针扭转;同震阶段唐山断裂以右旋错动为主,并呈现西升东降的垂向变形特征;反演结果表明多段断层参与了唐山7.8级地震的应变释放,主震发震断层的南段破裂主要发生在深度约6~18 km、长轴约50 km范围,北段破裂主要发生在深度约7~17 km、长轴约30 km范围,滦县地震断层也发生了少量滑动;主震断层最大水平位错为7.82 m、最大垂直位错为2.04 m,等效震级为Mw 7.58。挖掘得到的唐山震前和同震三维形变场定量结果为深入研究该地震的孕育发生机制及地震影响提供了珍贵资料。     

东日本Mw9.0大地震前、同震及震后地壳水平运动

利用“中国地壳运动观测网络(二期)”245个GPS连续观测站及其他站的大量观测结果,采用离震中很远的7个GPS站作为位错参考框架,得到2011年3月11日东日本Mw9.0大地震的同震水平位移和应变。此种参考框架的位移场解,符合弹性位错模型离震中很远且水平位移接近于0的要求,不仅有利于位错模型反演,也有利于合理分析此次大地震同震水平位移的影响范围。地震前后区域参考框架下远场GPS连续观测位移时间序列和同震位移弹性位错模型反演结果表明,东日本大地震同震水平位移既是地震断层破裂的结果,也是远场震前数月或数年水平位移的弹性回跳,据此可以进一步确认这些站震前水平位移前兆异常。由于这些震前异常位移中亦包含了同震位移量(方向相反),因此,弹性位错模型为定量研究大地震前远场地壳运动异常提供了一种理论模型。位错参考框架中的同震位移及其位错模型和区域参考框架的位移时间系列的综合研究有助于对此次大地震地壳运动异常的认识。这种持续时间为数月至数年的弹性形变异常,可以为大地震的中短期预测提供依据。     

基于多源SAR数据的2022年门源Ms 6.9地震同震破裂模型反演研究

2022年1月8日青海省海北州门源县发生Ms 6.9地震,震中位于青藏高原东北缘祁连-海原断裂中段,属历史地震空区,基于多源合成孔径雷达（synthetic aperture radar,SAR）遥感数据研究该地震的破裂模式对理解青藏高原东北缘构造变形机制、应变释放过程以及地震危险性评估具有重要意义。首先利用Sentinel-1数据和合成孔径雷达差分干涉测量（differential interferometry synthetic aperture radar,D-InSAR）技术获取了门源地震的同震形变场,视线（line of sight,LOS）向形变场显示此次地震造成了约20 km长的地表破裂,最大形变约0.75 m;然后基于Sentinel-2卫星数据,利用光学影像配准和相关技术获取了本次地震的东西向同震形变场,最大同震位移达2.5 m;最后基于均匀弹性半无限位错模型,以LOS向形变场为约束反演了断层的滑动分布模型。结果显示,门源地震是一次典型的左旋走滑型地震,地震破裂主要集中在0~10 km深度范围,最大滑动量3.25 m,滑动角10.44°,对应深度4.89 km;反演给出的矩震量为1.07×1019 N·m,对应矩震级Mw 6.6。结合野外考察和地质资料,初步判定发震断裂为冷龙岭断裂,并引起托莱山断裂发生同震滑动。同震库仑应力结果显示,冷龙岭断裂东段和托莱山断裂西段应力状态为加载,未来具有发生强震的风险。