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年内余震集中分布在泽普断裂库仑应力增加区域, 表明皮山地震主震对余震的发生可能具有一定的应力触发作用。     

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

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

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

InSAR观测与小震分布联合分析2021年云南漾濞MS6.4地震应力卸载对周围区域的影响

为分析2021年5月21日云南省漾濞M_S6.4地震后震区应力变化对周围断层的影响,本文通过InSAR技术获得了漾濞地震的同震形变场,并联合小震分布数据建立断层破裂滑动模型,继而通过计算断层面上的同震库仑应力来评估此次地震对周边断层的影响,以便有效地分析地震破裂的时空解析度.结果显示:(1)在升降轨InSAR数据获得的精细同震形变场中,升轨最大视线向形变量约为5.00 cm,降轨最大视线向形变量约为7.80 cm;(2)余震精定位的主震震中为(99.89°E,25.67°N),震源深度为13.29 km,除主震之外震源深度主要集中在5—15 km;通过小震位置拟合出的发震断层走向为NW-SE(316.69°),断层倾角为88.56°,滑动角为177.97°;(3)基于InSAR同震形变场结果及小震拟合断层参数联合反演得到此次地震的断层滑动以右旋走滑为主,升轨断层最大滑动量为0.80 m,对应的深度为8.85 km,平均滑动量为0.22 m,矩震级为M_W6.41;降轨的断层最大滑动量为0.30m,对应的深度为6.88 km,平均滑动量为0.05 ...     

2021年云南漾濞MS6.4地震序列静态应力触发研究

为深度剖析2021年云南漾濞M_S6.4地震对周围断层的影响以及地震序列间的静态库仑应力影响关系,文章基于主震的破裂模型和Okada给出位移对空间偏导数的解析式,首先计算了主震在周围断层上产生的静态库仑破裂应力,结果表明维西—乔后断裂带中段、澜沧江断裂带北端、红河断裂带北端以及怒江断裂带中段库仑应力均有千帕量级的增加。其次,计算了此次地震对周围地区产生的水平应力场及位移场,发现震中东西两侧物质向外流出,南北两侧向震中汇聚;震中南北两侧沉降,东西两侧隆升;产生的应力场呈EW向挤压,NS向拉张,在一定程度上抵消了该区域背景构造应力场。最后计算了前震-主震-余震序列间的静态库仑应力影响,结果表明前震产生的静态库仑应力促进了主震的发生;在2 km、13 km和18.5 km深度附近,触发的余震(前震和主震产生的库仑应力变化为正)比例很高,但在7 km深度处(同震破裂模型中滑移量最大)大部分余震分布在库仑应力负值区(应力影区),考虑到该深度余震与主震震源机制相差较大,因此通过模拟最易错动的断层面作为余震接受断层面,从而计算出最大静态库仑破裂应力,发现应力影区的余震仍有被触发的...     

基于同震库仑应力变化的汶川地震余震频次研究

在计算同震静态库仑应力变化的基础上,基于速率-状态摩擦定律及与同震静态库仑应力变化相关的余震频次估计方法,研究了汶川地震余震区及附近不同区域内与汶川地震同震静态库仑应力触发相关的"直接"余震频次。结果显示,主震破裂带尤其是主震破裂带南段的"直接"余震频次明显低于实际情况,其原因在于同震库仑应力变化导致的主震破裂面上应力水平的降低。研究结果还显示,汶川地震同震库仑应力变化对余震活动持续时间的影响与震级下限有关。对ML4.0以上余震,持续时间约为震后15~16个月;对ML3.5以上余震,接近60个月。在上述时段内,与同震库仑应力变化相关的"直接"余震占全部余震的比例为44.7%~48.6%。这表明,即使在同震库仑应力变化的"有效"作用时段内,主震破裂面上也大约有50%的余震活动不是缘于同震库仑应力变化的影响,这可能与震后余滑及粘弹松弛等时间相关因素的影响有关。     

基于同震库仑应力变化的汶川地震余震频次研究

在计算同震静态库仑应力变化的基础上,基于速率-状态摩擦定律及与同震静态库仑应力变化相关的余震频次估计方法 ,研究了汶川地震余震区及附近不同区域内与汶川地震同震静态库仑应力触发相关的“直接”余震频次。结果显示,主震破裂带尤其是主震破裂带南段的“直接”余震频次明显低于实际情况,其原因在于同震库仑应力变化导致的主震破裂面上应力水平的降低。研究结果还显示,汶川地震同震库仑应力变化对余震活动持续时间的影响与震级下限有关。对ML4以上余震,持续时间约为震后15~16个月;对ML3.5以上余震,接近60个月。在上述时段内,与同震库仑应力变化相关的“直接”余震占全部余震的比例为44.7%～48.6%。这表明,即使在同震库仑应力变化的“有效”作用时段内,主震破裂面上也大约有50%的余震活动不是缘于同震库仑应力变化的影响,这可能与震后余滑及粘弹松弛等时间相关因素的影响有关。     

完全库仑破裂应力变化与云南龙陵震群序列的应力触发

计算和研究了1976年云南龙陵Ms73、Ms74双主震产生的完全库仑破裂应力变化及其对后续震群序列的动态和静态应力触发作用.结果显示,龙陵第2主震受到了第1主震的动态和静态库仑破裂应力的触发作用.龙陵双主震的13个后续强余震中,发生在龙陵三角形块体内部的强余震,90%都受到了第1或第2主震的动态库仑破裂应力触发或静态库仑破裂应力触发作用,发生在三角形块体外部的强余震,2/3受到了第2主震的动态库仑破裂应力触发作用.从触发作用的强度和范围上看,第2主震的动态和静态库仑破裂应力触发作用都大于第1主震.就近场而言,受到龙陵双主震的动态和静态库仑破裂应力触发作用的余震数目相当,但动态库仑破裂应力的触发作用范围大于静态库仑破裂应力.     

2021年云南漾濞MS6.4地震动态应力对后续余震活动的触发作用

选取IRIS远震台站波形数据,反演了云南漾濞M_S6.4地震震源破裂过程,计算了断层破裂在近场产生的动态库仑破裂应力变化,并讨论了主震对近场余震活动的动态应力触发作用。结果显示:动态库仑应力演化过程与震源破裂特征反演结果一致,其大小分布与地震序列分布的疏密程度也具有较好的相关性。主震产生的静态和动态库仑破裂应力均促进余震的发生,但相比静态应力,余震位于库仑破裂应力正值区域的比例提高了21%,余震与动态库仑应力变化的正负区域有更好的一致性,动态应力能更好地解释震后余震分布的空间特征。垂直于地震序列主干10 km处出现小震丛集,该现象可能是由主震产生的动态库仑破裂应力占主导作用所致。定量分析主震对余震的动态应力触发结果显示,主震后一周内M_S4.0以上的8次余震接收点均受到了动态库仑破裂应力的触发作用。      

基于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,与美国地调局和哈佛大学给出的震源机制解基本一致.同震库仑应力导致了虎牙断裂北段延长线的东北和西南两端应力增强,其中塔藏断裂的罗叉段和马磨段未来强震的危险性值得关注.     

Stress drop assessment of the August 8, 2017, Jiuzhaigou earthquake sequence and its tectonic implications*

By using a broadband Lg attenuation model developed for the Tibetan Plateau, we isolate source terms by removing attenuation and site effects from the observed Lg-wave displacement spectra of the M7.0 earthquake that occurred on August 8, 2017, in Jiuzhaigou, China, and its aftershock sequence. Thus, the source parameters, including the scalar seismic moment, corner frequency and stress drop, of these events can be further estimated. The estimated stress drops vary from 47.1 kPa to 7149.6 kPa, with a median value of 59.4 kPa and most values falling between 50 kPa and 75 kPa. The estimated stress drops show significant spatial variations. Lower stress drops were mainly found close to the mainshock and on the seismogenic fault plane with large coseismic slip. In contrast, the highest stress drop was 7.1 MPa for the mainshock, and relatively large stress drops were also found for aftershocks away from the major seismogenic fault and at depths deeper than the zone with large coseismic slip. By using a statistical method, we found self-similarity among some of the aftershocks with a nearly constant stress drop. In contrast, the stress drop increased with the seismic moment for other aftershocks. The amount of stress released during earthquakes is a fundamental characteristic of the earthquake rupture process. As such, the stress drop represents a key parameter for improving our understanding of earthquake source physics.     

利用InSAR获取2020年新疆于田Mw6.3地震同震形变场与断层滑动分布反演

本文基于InSAR技术,利用欧空局Sentinel-1A/B升降轨SAR数据,提取了2020年6月26日新疆于田Mw6.3地震的同震形变场.利用升、降轨同震形变场约束,分别采用MPSO算法和Bayesian方法反演此次地震发震断层均匀滑动的几何参数,并进行对比.然后采用SDM方法获得发震断层非均匀滑动分布,并分析了同震...     

用形变资料反演1976年唐山地震序列的破裂分布

1976年唐山发生了7.8级地震,相继又发生了两次大余震——滦县7.1级地震和宁河6.9级地震.地震发生在观测条件比较好的地区,水准测量和三角测量测得了地震的同震位移场.本研究采用原始水准测量数据,而不是采用根据水准数据处理的地面沉降图像,和三角测量数据反演了该地震序列的破裂分布.模型构建中考虑了滦县地震和宁河地震的断层形态和大小.结果表明,唐山地震主震断层有明显的右旋走滑性质,最大走向滑动错距>6 m,位于断层南段,北段的走滑分量明显小于南段.主震总地震矩达2.58×10²⁰N·m,与地震波反演得到的地震矩的量级相当;滦县地震断层总体表现为左旋正断层,释放地震矩达4.95×10¹⁹N·m;宁河地震断层总体表现为右旋正断层,释放地震矩达3.94×10¹⁹N·m,比地震波反演的地震矩大一个量级.据此可以推测唐山地震的无震滑移主要发生在宁河地震断层的西部上,滑动性质以正断层为主.该结果对于唐山地震序列后的动力学演变过程及余震发生机理有一定参考.     

2014年美国加州纳帕M_W6.1地震断层参数的Sentinel-1A InSAR反演

2014年8月24日,在美国加州旧金山海湾北部的纳帕地区发生了MW6.1地震.发震断层是西纳帕断裂系统中的一部分,但是该断层之前并未被足够重视.本文利用欧洲空间局最近发射成功并刚刚投入使用的Sentinel-1A卫星获取的第一对同震干涉像对(20140807-20140831),得到了该地震的地表同震形变场,结合震后24h内区域GPS同震形变资料作为约束条件,反演了纳帕地震的断层几何参数以及滑动分布.Sentinel-1A干涉结果表明,此次地震造成了明显的地面形变,视线向最大抬升和最大沉降量均达到了10cm.联合反演结果表明,该发震断层的走向为344°,倾角为80°.主要破裂以右旋走滑为主,平均倾滑角为-146.5°,最大倾滑量达到了1.1m,位于地表下约4km,存在明显的滑动亏损现象.此次地震,累计释放地震矩达1.5×1018 N·m,约合矩震级MW6.1.该结果略小于InSAR单独约束结果,可能与Sentinel-1A像对中包含的快速震后形变分量有关.     

2015尼泊尔M_S8.1地震中等余震震源机制研究

本研究利用西藏台网记录的波形数据,采用gCAP方法反演了2015年4月25日尼泊尔MS8.1大震5次中等余震(5.0≤MS≤6.5)及西藏定日MS5.9地震震源机制解.结果显示,6次地震包含2个正断、2个走滑及2个逆冲型地震.其中2个正断型地震位于主震的东北方向,即发震断层的上盘,表明该区域受到主震同震位移的影响,表现出应力拉张的变化特征;2个走滑型地震在主震破裂的东南方向上,说明随着破裂往东南方向延伸,余震的走滑分量增强;另外2个逆冲型地震位于5月12日MS7.5强余震区域,与MS7.5地震的滑移状态一致,可能与主震同震位移引起该区域处于应力挤压状态密切相关.这些结果表明,尼泊尔MS8.1主震发生后,由于同震位移的影响,不同区域处于不同的应力状态,从而使中等余震表现出不同的震源类型.     

2017年MW6.5九寨沟地震激发的同震库仑应力变化及其对周边断层的影响

2017年8月8日21时19分在我国四川省北部阿坝州九寨沟县发生了M_W6.5左旋走滑型地震.该地震发生在青藏高原巴颜喀拉块体东北缘,东昆仑断裂南东段的塔藏断裂、岷江断裂和虎牙断裂的交汇地带.包括此次地震,近年来在巴颜喀拉块体周缘已发生了九次6级以上强震,表明巴颜喀拉块体周缘主要活动断裂上的应力水平仍处于不断调整之中.本文采用库仑应力模型研究2017年M_W6.5九寨沟地震激发的库仑应力变化、该地震与周边地震的应力触发关系以及强震对周边主要活动断裂的应力扰动.强震序列包括周边区域1536-1975年M≥6历史强震和1976-2017年的M_W≥ 6 gCMT地震目录中的强震,共计32个.研究结果表明：（1）2017年M_W6.5九寨沟地震激发的同震库仑应力变化仅在局部范围内超过0.1×10⁵Pa,且75%的余震（~12.7天）受到该地震明显的同震应力触发作用,而其余25%的余震落在应力影区,采用最优破裂面可以进一步提高同震库仑应力变化与余震分布的空间相关性;（2）2008年M_W7.9汶川地震对2017年M_W6.5九寨沟地震的发生有一定的促进作用,在后者震源处激发的同震库仑应力变化为（0.026~0.263）×10⁵Pa,震后黏弹性库仑应力变化为（0.010~0.032）×10⁵Pa.该库仑应力的变化范围取决于汶川地震源断层参数和九寨沟地震接收断层参数.2013年M_W6.6芦山地震对九寨沟地震的发生几乎没有影响（< 0.001×10⁵Pa）;（3）1654年M8.0甘肃天水南地震对九寨沟地震的发生有明显的促进作用,在九寨沟地震震源处激发的同震库仑应力变化为（0.410~1.266）×10⁵Pa,震后库仑应力变化为（0.147~0.490）×10⁵Pa.1879年M8.0甘肃武都地震可能有比1654年M8.0甘肃天水南地震更强的应力触发作用,但也有可能对九寨沟地震的发生起到抑制作用.在选取的8个九寨沟地震接收断层面上,其中6个接收断层面上该地震所激发的同震库仑应力变化为（0.913~2.364）×10⁵Pa,2个接收断层面上该地震所激发的同震库仑应力变化为（－1.326~－0.454 ）×10⁵Pa;在4个接收断层面上震后库仑应力变化为（0.094~1.072）×10⁵Pa,在另外4个接收断层面上震后库仑应力变化为（－1.593~－0.106）×10⁵Pa.1933年四川叠溪地震对九寨沟地震的发生影响较弱,其所激发的同震库仑应力变化为（0.015~0.080）×10⁵Pa,震后库仑应力变化为（－0.029~0.025）×10⁵Pa;（4）九寨沟地震仅在其附近的岷江断裂北段、塔藏断裂和虎牙断裂南段造成较明显的同震库仑应力变化,其分别为0.09×10⁵Pa、（0.14~2.03）×10⁵Pa和0.25×10⁵Pa.而进一步顾及其余31个强震的库仑应力作用则发现,同震库仑应力增加非常显著的主要活动断裂分段为：岷江断裂北段南侧和岷江断裂南段的库仑应力变化分别升高5.6×10⁵Pa和9.8×10⁵Pa.鲜水河断裂北段南侧库仑应力升高23.0×10⁵Pa,鲜水河断裂南段道孚-康定段的北侧库仑应力升高9.0×10⁵Pa,而最南端库仑应力升高3.0×10⁵Pa;龙门山断裂带中段的北侧库仑应力变化为（6.1~7.4）×10⁵Pa,中段库仑应力增加（2.1~11.5）×10⁵Pa;西秦岭北缘断裂东段库仑应力变化为4.4×10⁵Pa;龙日坝断裂北段最北侧的库仑应力变化为2.0×10⁵Pa;小金河断裂北段库仑应力变化为1.7×10⁵Pa;安宁河断裂北段库仑应力变化为1.6×10⁵Pa;（5）由于下地壳和上地幔的黏弹性松弛作用,所有强震在九寨沟地震震后20年造成的黏弹性库仑应力变化在鲜水河断裂、龙门山断裂中段、塔藏断裂以及秦岭南缘断裂西段比较显著,其分别为：（1.0~3.0）×10⁵Pa、2.8×10⁵Pa、（2.3~2.7）×10⁵Pa和0.9×10⁵Pa.但总体上黏弹性库仑应力变化没有改变各断裂上的同震库仑应力变化空间分布.总的库仑应力变化在鲜水河断裂北段南侧和南段的道孚至康定段北侧、龙门山断裂中段北侧、岷江断裂南段和北段南侧、虎牙断裂、塔藏断裂以及西秦岭北缘东段很显著（均超过4×10⁵Pa）.由于库仑应力明显升高可能预示着地震潜在危险性增强,因此这些断裂分段可能将来需要重点加以关注.