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;发震断层上静态库仑应力增加区域与余震分布具有一致性,说明余震主要是由静态库仑应力加载而触发的.     

GNSS观测的2021年5月22日玛多MS7.4地震同震位移及其约束反演的滑动破裂分布

2021年5月22日青海省果洛州玛多县发生M_W7.4地震,此次地震产生的地表破裂在空间上表现出明显的分段特征.本文基于不同来源的GNSS连续观测网数据获取了此次地震的精细三维同震形变场,结果显示:观测到的最大水平位移量达到280 mm,最大垂直形变量仅为25 mm,暗示此次地震的逆冲分量较小;此次地震具有较为明显的左旋走滑特征,同震形变基本对称,在NW-SE向的影响范围更广,该方向上水平同震形变大于3 mm的震中距范围超过500 km.进而,本文以余震精定位结果和GNSS观测到的三维同震形变场为约束,构建了地表破裂线为折线、倾角为85°、倾向西南的断层模型,反演了滑动破裂分布.结果显示:滑动破裂分布在震中两侧不均匀,均破裂到地表,破裂深度达到15 km左右,最大滑移量为4.73 m,计算的矩震级为M_W7.37.该结果与余震精定位结果具有很好的一致性,破裂的极值区正好位于早期余震空区,推测该余震空区未来的发震风险性较低.最后基于反演结果模拟计算了震中区域形变和应变场,结合应变值在断层地表迹线东南侧呈现挤压特征和已有的研究成果,推测此次地震增强了巴颜喀拉块体在东部地区挤压应力的积累特征,导致东部地区发震危险性增强,值得后续跟踪研究.     

2020年伽师MS6.4地震InSAR同震形变与滑动分布特征分析

利用Sentinel-1A升轨和降轨数据,基于D-InSAR技术,获取2020年1月19日伽师M_S6.4地震同震形变场,并结合其他研究机构给出的震源机制解参数和已有研究成果,反演得到伽师地震的发震断层几何特征和滑动分布。研究结果表明,伽师地震同震形变在地表有明显差异;升轨同震形变在卫星视线方向北侧抬升55 mm,南侧下降42 mm;降轨同震形变在卫星视线方面北侧抬升63 mm,南侧下降23 mm。通过反演得到发震断层走向为275°,倾角为20°,地震滑动主要分布在地下5 km处,最大滑动量约为0.32 m,平均滑动角为89.3°,累积地震矩为1.46×1018 N·m,合矩震级M_W6.1,发震构造为具有少量走滑性质的逆冲断裂。从发震构造特征、同震滑动分布推测,伽师地震发震构造是柯坪塔格褶皱带滑脱面以上沉积盖层内的逆冲断裂,支持了柯坪推覆体的薄皮构造模型观点。     

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地震相当.结合理论同震形变场、主应变场和邻近区域主要活动断裂库仑应力变化的空间分布特征,本次可能导致处于强闭锁和地震空区的安宁河断裂石棉—冕宁段未来发震风险性增强.

     

2021年5月21日Ms6.4漾濞地震GNSS同震变形场及其约束反演的破裂滑动分布

2021年5月21日云南大理州漾濞县发生了Ms6.4强震.地震未产生明显的地表破裂,但前震-主震-余震大致呈线状展布于维西—乔后—巍山断裂的西南侧.我们基于覆盖震中区域的高密度连续GNSS观测台阵,观测获得了震区周围50 km范围的高精度、高密度三维同震形变场.结果显示:地震破裂为明显的右旋走滑兼拉张;地表形变场呈现出...     

2021年青海玛多MW7.3地震InSAR的同震形变场、断层滑动分布及其对周边区域的应力扰动

文中基于D-InSAR技术,利用欧空局C波段升降轨哨兵SAR数据,获取了2021年5月22日青海玛多MW7.3地震的InSAR同震形变场,并对同震形变的空间特征、量级和断层破裂的分段性进行了分析.哨兵卫星的高质量观测数据清晰地描绘了玛多地震的地表破裂迹线,地表破裂长度约210km.为了进一步认识玛多地震断层深部的同震滑...     

2022年1月8日青海门源MS6.9地震GNSS同震形变场及其断层滑动分布

2022年1月8日在青海省海北州门源县发生M_S6.9地震,本次地震是继2016年门源M_S6.4地震后冷龙岭断裂周边发生的又一次强震。确定本次地震的破裂分布对分析该地区震害风险具有重要意义。通过收集震中及周边12个GNSS连续站点和震后加密观测的17个流动站点观测资料,获取了震中100 km范围内29个测站的GNSS静态同震形变场,并以此为约束反演了本次地震同震滑动分布。结果显示,近场GNSS观测到的最大形变量可达1.3 m。反演的最优破裂模型显示该地震主破裂区深度位于0～10 km,滑动破裂出露地表,最大滑动量为4.07 m,地震矩释放能量约1.1×10¹⁹ N·m,对应矩震级M_W6.7。门源地震破裂至地表是造成该地区基础设施破坏的直接原因。     

基于Sentinel-1A的新疆阿克陶MS6.7地震同震形变与滑动分布特征

文中利用Sentinel-1A升、降轨数据获取了2016年新疆阿克陶地震的同震形变场。形变场以水平运动为主,形变主要发生于断层南盘;升、降轨最大形变量分别约为12cm和-21cm;基于雷达影像观测右视成像的特点可知阿克陶地震具有右旋走滑的破裂特征,结合形变场形态特征与余震剖面推断认为,阿克陶地震的发震断层为S倾的木吉断裂。基于均匀弹性半空间模型反演双断层面的静态滑动分布,结果显示:2个断层面上各存在1个椭圆状的滑动集中区,破裂基本到达地表,滑动主要发生在沿断层面走向长约50km、沿断层面倾向方向0～20km的范围内,最大滑动量位于约9km深处,量级约为0.7m;西段主要以走滑为主,东段为走滑兼具少量正断性质,反演得到的地震矩约为8.81×10¹⁸N·m,相当于M_W6.57。综合分析认为,本次地震的震源特征为右旋走滑兼有少量正断分量,是发生在公格尔拉张系内拉张环境下的构造地震事件。     

GNSS约束的2022年泸定M6.8地震滑动分布及同震应力变化

2022年9月5日四川泸定县发生M6.8级强震,震中位于鲜水河断裂东南末端,野外地质调查初步结果显示本次地震并未发现明显的地表破裂迹象.本文基于震后科考GNSS流动观测和震中周边连续站观测资料,解算并提取震中90 km范围内31个测站的静态水平向同震位移.结果显示：GNSS同震形变场空间分布呈现明显的左旋走滑特征,观测到的最大水平向同震位移达23 cm,震中距50 km范围内同震位移量普遍大于1 cm.基于GNSS观测资料反演得到的同震滑动分布显示泸定地震地表破裂主要集中在磨西至田湾之间,主破裂深度2～8 km,最大滑动量～1.96 m,地震矩9.25×10¹⁸N·m,对应矩震级M_W6.6.静态库仑应力结果显示本次地震增强了震源破裂区周边活动断层的库仑应力,并触发了大量余震,余震主要发生在库仑应力增强区域.结合震间闭锁分布、历史地震及库仑应力变化,我们认为未来需要密切关注与磨西断裂交接的安宁河、大凉山断裂以及康定—磨西段的地震危险性.     

基于GNSS与InSAR约束的九寨沟MS7.0地震滑动模型及其库仑应力研究

2017年8月8日的九寨沟M_S7.0地震发生在岷江断裂、塔藏断裂及虎牙断裂交汇地区,地处青藏高原东北部的川甘交界地区,位于巴颜喀拉地块的东缘,地质构造复杂,对于九寨沟地震震中位置和发震断层的确定,存在不同意见.本文利用GNSS及升降轨InSAR观测,在获取九寨沟地震同震形变场的基础上,基于均匀弹性半无限位错模型,联合反演了发震断层的滑动分布模型,并计算了同震库仑应力变化.InSAR同震形变场显示,视线向最大沉降量和抬升量分别为0.21 m和0.16 m,形变场长轴为NW向,形变主要集中在断层西侧.距震中40 km和65 km的九寨和松潘两县,水平向的GNSS同震位移分别达14.31 mm和8.22 mm.联合GNSS和InSAR同震形变场反演得到的滑动分布主要集中在沿走向5~33 km,倾向2~20 km的范围内,平均滑动量为0.18 m,最大滑动量为0.91 m.发震断层长40 km,宽30 km,走向155°,倾角81°,滑动角-9.56°.同震位移场及滑移分布模型表明此次地震为一次左旋走滑为主的地震事件,地震破裂并未完全到达地表,与虎牙断裂北段的几何产状和运动学性质更为接近,结合精定位余震的分布,我们确定虎牙断裂北段为此次地震的发震断层,震中位于北纬33.25°,东经103.82°,震源深度10.86 km,矩震量为7.754×10¹⁸ Nm,相应的矩震级为M_W6.5,与美国地调局和哈佛大学给出的震源机制解基本一致.同震库仑应力导致了虎牙断裂北段延长线的东北和西南两端应力增强,其中塔藏断裂的罗叉段和马磨段未来强震的危险性值得关注.

     

2021年青海玛多Mw7.3地震地壳变形及区域孕震环境

北京时间2021年5月22日02时04分,在青海省果洛藏族自治州玛多县发生Mw7.3地震.本文利用震中及邻区2009-2021年GNSS观测资料,研究了此次地震的区域孕震环境、同震和震后初期的变形特征.首先,给出了较高空间分辨率的GNSS速度场,表明震中及其邻区的构造活动以左旋剪切为主,巴颜喀拉块体南北两侧的相对运动是...     

利用空间观测技术研究青海玛多7.4级地震孕育发生变形时空特征

为揭示2021年青海玛多7.4级地震震前变形与同震破裂特征,本文系统收集与处理了青海省CORS站和中国地壳运动观测网络等GNSS资料、以及Sentinel-1和ALOS-2卫星SAR数据,通过多源数据综合处理获取了玛多地震形变场动态结果.震前GNSS形变场结果显示,玛多震源区处于其北侧剪切形变高值区边缘,震前存在左旋剪...     

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年以上.     

汶川地震同震形变场的GPS和InSAR邻轨平滑校正与断层滑移精化反演

为克服InSAR观测汶川地震同震形变场的邻轨不连续问题,提出联合GPS观测值与邻轨平滑约束的同震位移校正方法,采用GPS观测形变去除PALSAR轨道误差引入的残留平地相位,基于形变平滑条件校正邻轨干涉相位的不连续性.ALOS/PALSAR干涉处理结果表明,校正后同震形变场的准确度与平滑性得以显著提高,InSAR高相干点残差达3.6 cm,校正后精度提高约60%,低相干点精度提高约40%,校正后形变场的邻接平滑因子标准差减小约33%,验证了轨道误差校正与邻轨平滑约束方法的准确性与可靠性.进一步基于弹性半空间位错模型的断层滑移反演结果表明,断层滑移主要分布于映秀、北川和青川地区,集中于地壳深度0~16 km范围,最大滑动量（位于北川县城）约为9.0 m,GPS反演模型残差为5.5 cm,InSAR反演模型残差达9.2 cm,InSAR反演精度约有30%的显著提高,由模型反演计算得到的地震矩为8.0469×1020 N·m.     

2021年青海玛多Mw7.5地震的同震变形及断层滑动模型研究

2021年5月22日,中国青海玛多县发生了M_W7.5地震.针对此次地震,国内外多家机构发布的震源机制解和有限断层滑动模型结果大都是基于半无限空间地球模型的位错理论反演得到的,未考虑地球曲率和层状效应的影响.该影响的量级以及其是否可以忽略目前仍是一个未知问题,值得研究.为此,本文利用美国地质调查局(USGS)、中国科学院青藏高原研究所以及中国地震局地球物理研究所提供的三个断层模型,基于弹性半无限空间、均质球、PREM三种地球模型的位错理论,分别计算了地表同震位移场和应变场,并对比分析了这些结果的差异.我们发现基于PREM位错理论的同震位移场与半无限空间模型的对应值差异约为3~28 cm,占PREM位错理论值的10%~30%,应变场差异更大,表明地球曲率和层状效应的影响不可忽略.此外,理论同震形变结果与GNSS和InSAR形变观测数据对比发现,基于PREM模型的理论位移场最接近于观测值.利用InSAR数据和三种地球模型反演得到的有限断层滑动模型存在差异,说明了地球模型的选择对断层模型的反演具有一定影响.本文的结果为今后对此次地震的观测数据物理解释和断层滑动反演提供了理论参考.     

2021年5月21日青海玛多MW7.4地震及余震序列应力降时空分布与孕震环境

2021年5月21日发生在青藏高原巴颜喀拉块体内部的玛多M_W7.4地震是继2008年M_W7.9汶川地震之后在中国陆区发生的最大地震.由于玛多地震发生在板块内部滑移速率较缓慢的次级断裂上,该地震成为研究巴颜喀拉块体内部地震危险性的典型震例.通过计算2021年玛多地震序列中震级大于2.5的地震应力降,我们调查了该地震序列应力释放的时空演化规律.首先,利用已有的青藏高原及其周边地区宽频带Lg波衰减模型,对玛多地震序列产生的地震Lg波观测数据进行了传播路径衰减校正,获得119个地震事件的Lg波震源谱;然后,通过将理论震源函数与观测数据拟合,获得对Lg波标量地震矩、拐角频率和高频下降率的估计;最后,根据标量地震矩和拐角频率计算了2021年玛多地震序列的应力降.结果表明,2021年5月21日玛多地震主震应力降为22 MPa,余震序列的应力降分布范围为0.08~7.5 MPa、中位数为0.39 MPa,平均值为0.88 MPa.应力降与震级具有较强的正相关关系,说明该地区的地震活动可能并不遵从地震自相似理论.换言之,大小地震具有不同的破裂动力学性质.玛多地震序列应力降变化反映了断层面上的应力释放过程具有强烈的非均质性.主震发生之后余震应力降显著降低,但其间仍夹杂着少数具有较高应力降的事件.这些高应力降事件所在区域与主震后破裂面上的闭锁和应力集中区域相对应.

     

青海玛多MS7.4地震对周边活动断裂的库仑应力加载及发震概率增量的计算

2021年5月22日2时4分位于青藏高原巴颜喀拉地块内部青海果洛州玛多县发生了M_S7.4地震(玛多M_S7.4地震),震中位于98.34°E,34.59°N,震源深度17 km(中国地震台网中心测定).震后区域应力调整直接导致了周边断层的库仑应力变化和主震断层外余震的发生,进而影响到区域地震活动性并改变了即发地震概率.本文利用有限断层滑移模型计算了青海玛多地区主要断层的库仑应力变化,结合Dieterich地震发生率模型,给出了库仑应力扰动下超越某一特定震级地震的发震概率公式,得到了受青海玛多M_S7.4地震所产生的库仑应力场的影响周边8条断裂(段)分别发生M_S≥7.0和M_S≥6.0地震的概率变化.结果表明,8条断裂(段)的发震概率均有不同程度的增加.对于发震概率增幅较为明显的甘德南缘断裂、玛多—甘德断裂和西藏大沟—昌马河断裂在震后短时间内(约10年内)发震概率迅速提升,之后趋于平稳,不排除有潜在发生破坏性地震的可能,特别是发生M_S≥6.0地震的可能性.针对达日断裂,由于库仑应力增加造成的影响不明显,短期内发生M_S≥7.0或M_S≥6.0地震的可能性不大,但随着时间的推移,不排除几十年后有潜在发生破坏性地震的可能,尤其要关注发生M_S≥6.0地震的可能性;东昆仑断裂特别是玛沁—玛曲段仍然是未来强震发生的可能区段,需重点关注和防范M_S≥6.0地震乃至M_S≥7.0地震的发生;玉树—甘孜断裂发生M_S≥6.0地震,尤其是发生M_S≥7.0地震的危险性不高,而乌兰乌拉湖—玉树南断裂有潜在发生破坏性地震的风险,对于M_S≥6.0以上地震需要加强防范.