2013年南斯科舍海岭MW7.8地震的多点震源机制反演

2013年11月17日,在南极南奥克尼群岛北、南极板块与斯科舍板块之间发生了一次M_W7.8级地震（2013年南斯科舍海岭M_W7.8地震）,我们利用全球分布的长周期和宽频带地震记录反演确定了这次地震随时间和空间变化的震源机制,验证了提出的一种多点震源机制反演的新方法.首先利用长周期记录的W震相反演了这次地震的矩心矩张量解并利用体波提取了视震源时间函数,同时利用台阵反投影技术从宽频带记录中获得了这次地震的高频源的时空分布,然后基于矩心矩张量解、视震源时间函数以及高频源的时空分布,实现了采用新方法对2013年南斯科舍海岭M_W7.8地震的多点震源机制反演.矩心矩张量解表明,地震矩心在44.50°W/60.18°S,矩心深度19 km,半持续时间49 s,释放标量地震矩4.71×10²⁰ N·m,发震断层走向104°,倾角54°,滑动角8°.视震源时间函数清楚地揭示了地震矩随时间变化的方位依赖性,总体上可以将时间过程分为前60 s和后50 s两个阶段,但前60 s可细分为两次子事件.根据台阵反投影结果,这次地震为沿海沟从西到东的单侧破裂,破裂长度达311 km,可以分为5次子事件,能量释放的峰值点依次为13 s、30 s、51 s、64 s和84 s,平均破裂速度分别为0.6 km·s^-1、2.6 km·s^-1、2.3 km·s^-1、2.8 km·s^-1和3 km·s^-1.多点震源机制反演显示,5次子事件的矩震级分别为M_W7.57,M_W7.48,M_W6.80,M_W7.53和M_W7.08,半持续时间依次为21 s,17 s,6 s,16 s和8 s,走向分别为95°,105°,81°,98°和98°,倾角依次为57°,49°,86°,46°和64°,滑动角-9°,1°,-17°,13°和-4°.这些在震源机制、能量释放以及持续时间方面的变化都是当地构造和应力环境复杂性的反映.     

2018年阿拉斯加湾MW7.9地震震源复杂性

2018年1月23日,在美国阿拉斯加湾海域发生了一次M_W7.9地震.震源机制解表明这次地震以走滑为主,可能发生在近东西向或南北向的陡倾角断层上,早期余震并非线型展布.我们利用视震源时间函数分析确定了此次地震的总体破裂方向,并结合余震的空间展布特征构建了相互交叉的双断层模型,进而通过联合反演远场P波和SH波数据获得了此次地震的时空破裂过程.视震源时间函数分析表明总体破裂方向既非东西也非南北,而且反演结果表明,两个断层上都发生了错动,总体破裂时间~50 s,释放标量地震矩~8.11×10²⁰ Nm.震源时间函数表现出多事件特征,且两个断层破裂的时间过程也不相同.破裂首先在南北向断层的南端开始,很快触发了东西向断层,最后终止于南北向断层的北端.每个断层都具有相当的时空复杂性,位错分布很不均匀.东西向断层具有三个凹凸体,一个位于震源附近,其他两个位于断层两端.南北向断层有两个凹凸体,均位于断层北段,最大滑动量~5.0 m就出现在这里.发生最大位错的南北向断层延伸至阿拉斯加海沟,增加了触发阿拉斯加海沟其他断层发生破裂的可能性.     

基于InSAR和远场地震波联合反演2008年MW6.3大柴旦地震震源破裂过程

利用InSAR同震形变升、降轨数据和远场地震波数据,基于均方根最小与标量地震矩最小双重约束下的模拟退火方法,联合反演2008年11月10日大柴旦M_W6.3地震震源破裂过程.结果表明,2008年大柴旦地震震源破裂过程整体表现为沿倾向方向从深部破裂起始点处向上往地表传播,且破裂未到达地表;在前7 s内,滑动沿西北和东南两个方向传播,7 s后主要沿东南方向传播,破裂过程时间持续约为11 s,同震滑动分布主要集中在地下10~20 km范围内,最大滑动量达-0.71 m;反演结果揭示本次地震为西南倾高角度兼具微量走滑分量的逆冲破裂事件,反演矩张量为3.96×10¹⁸N·m,矩震级约M_W6.37.通过大柴旦地震发震断层和破裂机制综合分析,初步判断发震断层为西南倾向的大柴旦—宗务隆山断裂.     

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

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

日本海沟俯冲带MW9.0地震震源区应力场演化分析

于2011年3月11日发生在日本东北部的M_W9.0级逆冲型板间地震是日本有地震记录以来震级最大的一次地震.本研究基于NIED F-net矩张量解目录中的震源机制解,选取两个长轴相互垂直的矩形区域进行应力场2D反演,获取了日本海沟俯冲带地区应力场的空间及时间分布图像.结果表明：主震前,俯冲带地区应力状态在空间上大体趋于一致,即应力轴（P轴、σ₁轴及S_Hmax轴）系统性地倾向板块汇聚方向,P轴、σ₁轴倾角整体偏缓（<30°）,且远离震源区及日本海沟东侧区域内的应力轴倾角普遍大于主震震源区内应力轴倾角;主震前,受2003年5月26日在宫城县北部发生的M_W7.0地震影响,位于M_W9.0地震震源区西北侧的应力场出现明显扰动,σ₁轴倾向顺时针偏转150°~180°,并于之后大体恢复至震前状态,同期其他地区没有明显变化,这种情况可能和主震断层局部（深部）的前兆性滑动有关;主震后,距离震源区较远处应力场变化不大,主震震源区内应力场发生显著改变,P轴及σ₁轴均以大角度（>60°）倾伏于板块汇聚方向,S_Hmax轴顺时针偏转60°~90°且在日本海沟附近普遍平行于海沟轴.这项研究以时空图像的方式展示了大地震前应力场变化的特点,反映了大地震孕震过程中构造与地震的相互作用,对于理解大地震孕震过程有重要意义.     

2016年台湾美浓MW6.4地震震源参数的InSAR和GPS反演

2016年2月6日台湾西南部高雄市美浓区发生了M_W6.4地震.本文结合ALOS2卫星升降轨、Sentinel-1A升轨SAR数据,采用两轨差分干涉技术获取了该区域的同震形变场,形变结果表明震中西北部以抬升为主,最大视线向形变量约为11.2 cm.基于均匀位错模型和多峰值粒子群（MPSO）算法,利用InSAR和GPS形变数据联合反演了美浓地震的断层几何参数,结果表明震源中心位于22.920°N,120.420°E,深度约12 km,发震断层长度约15 km,走向角307°,倾角16.5°,平均滑动角为51.5°,此次地震是以逆冲倾滑兼左旋走滑的破裂模式.利用格网迭代搜索法得到最优倾角为15.7°,GPS和InSAR最优权比为18：1,最优平滑因子为0.06.基于非均匀位错模型,利用非负最小二乘方法进行线性反演,结果显示最大倾滑和走滑量分别为51.7 cm和55.3 cm,对应矩震级为M_W6.38,略小于GCMT （M_W6.4）的结果.通过与已有文献的比较和对该区域断层构造的分析,发现美浓地震的发震断层为单一断层的解释更为合理,我们推测发震断层是位于左镇、后甲里等断层之间的一条东南-西北走向往东北倾斜的盲断层,并初步推测2010年M_W6.3甲仙地震也同该断层有关.     

内蒙古中西部地区震源参数研究

选取内蒙古测震台网2009—2016年内蒙古中西部M_L≥2.8中小地震波形记录,采用多台联合反演方法,计算该区中小地震震源谱参数,获得地震矩M₀、矩震级M_W、应力降Δσ和震源尺度r,利用线性回归,分析近震震级M_L、地震矩M₀、矩震级M_W、应力降Δσ和震源尺度r之间的关系。结果表明,各参数之间存在一定线性关系,内蒙古中西部地区应力降模型为增加应力降模型（ISD）。     

渭河断陷盆地及邻区环境剪应力场及其与震源参数关系研究

利用陕西数字地震台网的数字地震资料研究了渭河断陷盆地及邻区的震源参数和环境剪应力场, 结果表明该地区地壳内环境剪应力处于较低水平, 平均值为12.7×10⁵ Pa, 并据此分析了该地区的地震活动性, 建立并探讨了环境剪应力与震级M_L、 地震矩M₀、 矩震级M_W的关系, 讨论了环境剪应力与震源深度间的关系。     

2018年斐济MW8.2深震:高频辐射过程和发生原因

2018年8月19日,在斐济东部海域563 km深处发生了一次M_W8.2地震.我们首先挑选位于美国阿拉斯加地区的131个宽频带台站构成台阵,选用垂直分量0.5~2 Hz的高频信号,利用广义台阵反投影技术对这次地震的破裂过程进行了成像,然后基于破裂速度对地震的辐射效率进行了估计.结果表明,这次地震总体上呈单侧破裂,破裂方位在3.0°左右,破裂总长度约51 km,持续时间22 s,平均破裂速度为2.5 km·s^-1.但能量释放有2次高峰,形成两次子事件.第一次为前10 s,峰值在7 s左右,破裂速度为2.9 km·s^-1,辐射效率为45%.第二次为10~22 s,峰值在15 s左右,破裂速度为1.6 km·s^-1,辐射效率为26%.结合震源位置、震源机制、破裂速度以及辐射效率,我们认为这次地震是由于俯冲板块前缘受到下部地幔物质上浮阻力引起的剪切失稳所致,起初板块内部的脆性破裂表现突出,致使辐射效率较高,后来震源处高温高压下的熔融耗散特征逐渐凸现,致使辐射效率下降.     

2022年四川泸定MS6.8地震震源基本特征

2022年9月5日,在四川省甘孜州泸定县发生M_S6.8地震.本研究收集区域台网震相数据、全球地震台网(GSN)、国际数字地震台网联盟(FDSN)与德国地学中心GEOFON台网的宽频带P波数据,利用双差定位、矩心矩张量解反演、有限断层波形反演和视震源时间函数分析等方法,分析了此次地震震源的基本特征.重定位结果显示,余震分布具有丛集性,3个地震丛分别集中在震中附近以及相距30 km左右的东南端和西北端,总体上西北端地震较浅,东南端地震较深.矩心矩张量反演表明,矩心位于29.55°N,102.14°E,深度16 km,释放地震矩1.0068×10¹⁹ N·m,相当于矩震级M_W6.60,双力偶成分占88%,是一次近纯走滑的地震事件.结合余震分布可以断定,地震发生在走向163°、倾角77°(倾向西南),滑动角为-5°的断层面.有限断层反演显示,破裂区主要由两部分构成,破裂起始点及其周围是主要破裂区;另一破裂区位于其东南,总体表现为从西北向东南的单侧破裂,最大滑动量约1.4 m,位于起始破裂点附近.从矩心矩张量反演和有限断层反演得到...     

智利M_W8.3地震与M_W8.8地震的震源过程及其引起的库仑应力特征

2015年9月17日6时54分32秒(北京时间)智利中部伊拉佩尔附近(震中31.57°S,71.67°W)发生了一次M_w8.3大地震,在此次地震震中以南约500 km处的马乌莱地区曾于2010年2月27日14时34分11秒发生过一次M_w8.8强震(震中36.12°S,72.90°W),两次地震余震分布区之间有约75 km的地震空区.本文利用远场体波与面波波形,基于有限断层模型,反演了这两次地震的震源破裂过程.结果显示这两次地震均为逆冲型大地震,2015年伊拉佩尔M_w8.3地震的平均滑动角度为107°,平均滑动量为2.43 m,平均破裂速度为1.82 km·s~(-1),标量地震矩为3.28×10~(21)Nm,95%的标量地震矩在104 s内得到了释放.最大滑动量约8 m,位于沿走向75 km,深度8 km处.2010年马乌莱M_w8.8地震的平均滑动角度为109°,平均滑动量为4.95 m,平均破裂速度1.90 km·s~(-1),标量地震矩为1.86×10~(22)Nm,95%的标量地震矩在121 s内得到了释放.最大滑动量约12.5 m,位于沿走向100 km,深度21 km处.2015年伊拉佩尔M_w8.3地震浅部更大的滑动量应该是其引起了较大海啸的一个原因.基于破裂滑动分布,我们计算了这两次地震引起的周边俯冲带上静态库仑应力变化,结果显示两次地震均显著增加了周边俯冲带上的库仑应力,2010年马乌莱地震使得2015.年伊拉佩尔地震震源区附近的库仑应力增加了(0.01~0.15)×10~5Pa,从应力积累的角度看,2010年马乌莱地震有利于2015年伊拉佩尔地震的发生,对后者的发生起到了促进作用.     

2016年印尼苏门答腊岛海域MW7.8地震震源运动学特征

根据中国和全球地震台网记录的波形记录,采用W震相矩张量反演、反投影分析及有限断层模型反演方法,研究了2016年3月2日印尼7.8级地震破裂过程,分析讨论印尼地震震源运动学特征.结果表明：此地震为一次对称的双侧破裂走滑型事件,北北东─南南西向的断层节面（走向5°/倾角85°）为发震断层面.标量地震矩约6.19×10²⁰ Nm,矩震级为7.79,最大的滑动量约11 m,位于破裂起始点北东,沿着断层走向约30 km处.破裂平均速度2.0~2.2 km·s^-1,破裂持续时间35 s,破裂在5~25 s内释放的能量,约占总能量的97%.最终形成了总长度90 km左右的断层.印尼地震具有破裂持续时间短、破裂速度慢、高滑动能量带相对集中等显著特点.本研究对进一步增进海洋岩石圈地震的震源特性认识有重要参考意义.     

青藏高原东北缘断层系统的大地震迁移概率及断层滑动速度的分段特征

青藏高原东北缘是青藏高原隆升、生长及变形前缘.区域地震活动频繁,且地震在其主要断层带之间时空迁移.为了研究区域大地震在主要断层带之间的迁移规律与概率,以及主要断层带大地震破裂的时空分布特征,本文建立了青藏高原东北缘地区的三维黏弹塑性有限元模型,模拟了区域断层系统的地震循环,得到了人工合成的万年时间尺度的地震目录.根据模拟的地震目录,并结合古地震数据,计算分析了大地震(MW≥7)在研究区各个主要断层带之间的迁移概率,探讨了黏度、高程、统计时间长度等因素对大地震在各主要断层带之间的迁移概率和大地震在各主要断层带上的发生概率的影响,并且初步调查了海原断层带和香山天景山断层带的大地震破裂时空分布特征.研究结果显示:继区域最近两次大地震(1920年海原断层带上的M8.5海原大地震和1927年香山天景山断层带上的M8古浪大地震)之后,下一次大地震(MW≥7)发生在海原断层上的概率最大,约为51%～81%;其次是在香山天景山断层上,概率约为9%～37%.模型结果显示,不同的青藏高原中下地壳上地幔黏度大小,对大地震在各个断层带之间的迁移规律和迁移概率的影响较小;而研究区的高程载荷对地震迁移则有显著的影响:高程载荷易于使得海原断层地震活动减弱及香山天景山断层的地震活动增强.研究结果也显示了青藏高原东北缘地区主要断层带的地震活动与断层滑动速率分布的分段性显著;大地震在断层带上的破裂位置并不固定,呈现不均匀性;并暗示了断层几何形状对地震活动、断层滑动速率分布与大地震破裂位置的控制作用.     

2010智利莫尔M_W8.8地震对其破裂区及附近地区地震活动性的影响

强震的发生与地壳中应力增加有关。库仑破裂应力被用于研究由于某次强震发生所产生的应力变化。2010年智利莫尔M_W8.8地震发生之后,其破裂区以南地区没有发生M≥6.0地震,而在其以北地区于2015年在智利发生了伊亚佩尔M_W8.3级地震。在这两个地区地震活动性悬殊。本文利用库仑破裂应力研究了2010年智利莫尔M_W8.8地震在这两个地区产生的应力变化。结果表明,2010年智利莫尔M_W8.8地震在其破裂区以北地区产生的应力变化为正值,有利于断层错动发震,而在其破裂区以南地区产生的应力变化为负值,具有抑制断层错动发震的作用。     

2011年日本MW9.0地震对沂沭断裂带及其两侧地区地壳运动的同震影响研究

2011年日本M_W9.0地震（简称日本地震）后沂沭断裂带及其两侧地区地震活动显著增强,研究日本地震对该地区地壳运动及地震潜势的影响十分必要.为此,本文通过112个连续GPS观测站获取了研究区高空间分辨率的日本地震同震形变场并得到如下认识：（1）8个定点地球物理观测的同震响应验证了本文同震形变场的可靠性;日本地震的东向拉张使研究区整体上处于张性同震应变状态,但存在局部挤压区域,其中莱州湾至海州湾的挤压条带穿过沂沭断裂带并对断裂带南北两段产生了不同的同震作用,对南段具有拉张作用,对北段产生挤压作用;（2）同震形变场在鲁东隆起和鲁西断块产生了显著的剪应变,地震b值显示上述区域的构造应力在日本地震后增强,因此同震形变场可能改变了这些区域的应力特征;（3）地震矩张量叠加分析显示,同震形变场短期内对鲁西断块、鲁东隆起区和沂沭断裂带南段累积了地震矩,可能有助于上述区域在日本地震以后的地震活动增强;日本地震对沂沭断裂带北段的地震矩具有释放作用,或许是该区域地震活动减弱的原因.     

用连续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,而极震区以南的深部滑脱断层仍保持闭锁,未来仍有导致灾害性大震的可能性.     

基于区域地震波形的2017年新疆精河MS6.6地震破裂方向性及发震构造研究

2017年8月9日的新疆精河MS6.6地震是近年来天山北缘发生的最大地震,震中位于由多条逆冲断层组成的库松木契克断裂带内.由于震源较深、构造形变复杂、区域地震台站相对稀疏,仅根据震源机制解、余震分布和InSAR观测结果等难以直接判定发震构造.本文针对倾滑型地震发展了一种基于区域地震波形的破裂方向性测定方法,利用余震作为参考地震进行路径校正,根据主震和参考地震的波形时移差和Pn-Pg到时差分别确定主震在水平方向和深度方向的破裂尺度,进而推断同震破裂的延展方向和延伸尺度.本文在反演了主震的点源参数后,应用新发展的方法测定了地震的破裂方向性.点源反演结果显示,精河地震是一个发生在中地壳的高角度逆冲地震,矩震级约6.2,质心深度21km,震源持续时间5.5s,两个双力偶节面分别为102°/45°/106°(NP1)和259°/47°/74°(NP2).破裂方向性分析结果显示,地震的破裂面为南倾的NP1节面,地震沿着破裂起始点向西南方向、向下破裂,总破裂长度约11.5km,其中,沿深度的破裂范围约7km,沿水平的破裂范围约9km,平均破裂速度约2.1km·s-1.综合区域地质资料、卫星影像等判定本次地震的发震断层为精河南断层,地震可能只破裂了断层的下段(17~25km),并未破出地表.