利用区域宽频地震数据反演2021年5月云南漾濞MS6.4地震震源破裂过程

本文基于有限断层模型反演方法,利用区域宽频带数据反演了2021年5月云南漾濞M_S6.4地震的震源破裂过程,结果显示:此次地震的发震断层走向为SE向,主要以右旋走滑为主.破裂主要发生在震源东南侧,最大错动量约为0.55 m,位于深度约9 km处,发生明显破裂的深度约达13 km.此次地震释放的标量地震矩为1.48×10¹⁸N·m,相当于矩震级M_W6.05.地震能量主要在前11 s释放.在深度为6~8 km处破裂速度有明显的变快,可能加剧了地表的震动.     

利用区域宽频地震数据反演2021年5月云南漾濞MS6.4地震震源破裂过程

本文基于有限断层模型反演方法,利用区域宽频带数据反演了2021年5月云南漾濞M_S6.4地震的震源破裂过程,结果显示:此次地震的发震断层走向为SE向,主要以右旋走滑为主.破裂主要发生在震源东南侧,最大错动量约为0.55 m,位于深度约9 km处,发生明显破裂的深度约达13 km.此次地震释放的标量地震矩为1.48×10¹⁸N·m,相当于矩震级M_W6.05.地震能量主要在前11 s释放.在深度为6~8 km处破裂速度有明显的变快,可能加剧了地表的震动.     

2020年Mw7.7加勒比海地震:反投影成像确定的一次超剪切事件

利用反投影方法,使用美国阿拉斯加台网的远场P波垂直分量地震数据,对2020年1月28日M_W7.7加勒比海地震的破裂过程分两个频段(0.2~1.0 Hz,0.5~2.0 Hz)进行成像.结果显示,地震沿西向单侧近线性破裂,破裂规模约为250 km,破裂时长约为98 s,平均破裂速率约为2.55 km·s^-1.破裂的瞬时速率存在明显波动,在大约53~62 s达到约6 km·s^-1,超过了所在深度的S波速度,为超剪切破裂.地震在5~15 s,21~30 s和66~79 s存在三个高频能量释放峰值,相应时段破裂速率较低,破裂能量得以在局部区域集中释放.破裂一直延伸至大洋中脊附近受阻挡而停止,转换断层几乎完全破裂,释放了加勒比板块和北美板块之间由于相对运动而长期积累的应力.发震断层特殊的空间几何特征和长期积累的应变能,为此次超剪切破裂事件的发生提供了断层几何和动力条件.     

区域宽频地震波形反演所揭示的2022年青海门源MS6.9地震震源破裂过程

基于有限断层模型反演方法,利用区域宽频带数据反演了2022年1月青海门源M_S6.9地震的震源破裂过程,并结合地质构造与地震重定位结果判断发震断层走向.综合反演结果表明:此次地震的发震断层走向为WNW向,主要以走滑为主;破裂主要发生在震源两侧,可能存在着双侧破裂,在震后2 s和9 s出现破裂极大值,最大错动量约为1.5 m,位于深度约6km处,发生明显破裂的深度约为16 km,地表破裂长度约20 km;此次地震释放的标量地震矩为1.23×10¹⁹N·m,相当于矩震级M_W6.7,地震能量主要在前15 s释放;发震断层面的倾角为84.6°,接近于垂直,由于破裂范围较大,所以发生明显错动分布的地表投影也长达34 km.     

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左右的断层.印尼地震具有破裂持续时间短、破裂速度慢、高滑动能量带相对集中等显著特点.本研究对进一步增进海洋岩石圈地震的震源特性认识有重要参考意义.     

四川芦山Ms7.0级强烈地震震源运动学特征

使用中国数字地震台网记录的区域宽频带波形,通过频率域和时间域多步反演,研究了2013年四川芦山“4·20”7.0级强烈地震的震源运动学特征.基于点源的震源机制解揭示：地震发震断层面参数分别为走向214°/倾角47°/滑动角96°,表现为一次高倾角的逆冲型事件.矩心在水平方向上位于震中（30.303°N/102.988°E）西南向约4.5 km,矩心深度约17 km.平均总标量地震矩M₀为1.16×10¹⁹ N·m,矩震级M_w约6.6.进一步模拟高达0.5 Hz高频波形,获得了芦山地震破裂过程图像,结果显示：此地震为一次不对称双侧破裂事件.破裂半径约15 km,整个破裂面积为706.7 km²,平均滑动量约0.231 m.破裂在8 s内释放了大多数能量.震后0~3 s内,破裂以孕震点为中心向四周同时扩展,3 s后,破裂表现出明显的方向性,主要向北北东扩展,导致位于震中北东向多数台站视破裂持续时间总体偏小,最小值为4 s.破裂约8 s后基本停止.     

利用区域宽频带数据反演鲁甸MS6.5级地震震源破裂过程

基于有限断层模型反演方法,我们利用区域宽频带数据反演得到了2014年8月3日鲁甸M_S6.5级地震的震源破裂过程.反演结果显示：此次地震的发震断层走向为北北西向,破裂主要以左旋走滑为主,位移主要发生在震源左上方,最大滑动量为0.7 m,模型显示断层破裂可能接近地表,破裂长度约10 km.此次地震释放的标量地震矩为1.97×10¹⁸ N·m,相当于矩震级为M_w 6.1,地震能量主要在前15 s释放.鲁甸地震有四个显著的特点：（1）位移主要集中在浅部,从11 km起破点开始迅速向上传播,大部分位于10 km以上且最大位移位于深度3 km处,从模型来看,破裂可能接近地表,因此地表震动较为强烈;（2）应力降比较大,计算显示释放的同震静态应力降约为2.8 MPa;（3）破裂速度较快,在地表附近超过了2.5 km·s^-1;（4）主震可能发生在一个共轭断层系上.这四个特点可能是导致此次地震造成如此重大人员伤亡和财产损失的最重要的原因.     

2014年云南鲁甸MS6.5地震震源运动学特征

使用中国数字地震台网记录的区域宽频带波形,通过频率域和时间域多步反演,研究了2014年云南鲁甸MS6.5地震震源运动学特征.基于点源的震源机制解揭示:地震发震断层面参数分别为走向342°/倾角83°/滑动角-34°,表现为一次左旋走滑型为主兼有少量正断倾滑分量的事件.质心在水平方向位于震中(103.354°E,27.109°N)东南约5.4km,最佳波形拟合的质心深度约4.4km,平均总标量地震矩M0为2.1×1018N·m,矩震级MW约6.1.破裂过程图像显示:此地震为一次不对称双侧破裂事件,破裂半径约10km,整个破裂面积为227.6km2,平均滑动量约0.16m.破裂在6s内释放了大多数能量,震后0～2s,破裂以孕震点为中心向NW和SE两侧同时扩展,2s后,破裂表现出明显的方向性,主要向SE(沿走向342°相反方向)扩展,故导致SE向多数台站视破裂持续时间总体偏小,最小值为2s.约6s后破裂基本趋于停止.推断鲁甸地震破裂在上地壳浅层集中释放了大多数能量是导致灾害严重的主要原因之一.     

2013年4月20日四川芦山7.0级地震震源破裂特征

2013年4月20日四川芦山发生7.0级地震.本文运用反投影远震P波的方法研究了中心频率为1 Hz的芦山地震震源破裂特征.研究结果显示2013年芦山7.0级地震破裂长度约为20 km,震源破裂时间约为26 s.本文认为在芦山地震的开始阶段(0~4 s)震源的破裂是向震中位置两侧进行的.而芦山地震破裂的第二阶段(5~26 s)是单侧破裂.芦山地震最大能量释放区域位于震中以北.本文对比了运用相同方法研究的发生在同一断裂带上的2008年汶川地震震源破裂特征.发现2013年芦山地震和2008年汶川地震有三点相似之处,即,破裂主要沿北东走向的龙门山断裂带发展;最大能量释放区域没有位于震中;能量都是通过多次子事件来释放的,且第二次能量释放是最大能量释放.对比两次地震破裂区域,可以看出芦山地震的破裂区域是在2008年汶川地震破裂区域的西南端发展的.两次地震的破裂区域占了整个龙门山断裂带的三分之二.     

利用多台阵压缩传感方法反演尼泊尔M_W7.9地震破裂过程

2015年4月25日,尼泊尔地区发生MW7.9地震,震中位于28.1°N,84.7°E.为了详细地研究此次破坏性极强的地震的破裂过程,本文利用多台阵压缩传感方法,使用了阿拉斯加、欧洲和澳大利亚三个台网的共计179个台站的远场P波垂直分量的数据来反演,结果表明本次地震的破裂过程是一个清晰的南东东方向的单侧破裂,破裂尺度约为105km,整体持续时间约为58s.在破裂初始的前15s,能量辐射基本围绕在震源附近,16s后破裂开始向南东东方向以1.9km·s-1的速度破裂.释放能量最大的时间为第38s,位于距震中70km处.该位置从第29秒开始破裂,并持续释放能量长达30s之久.     

A Rupture Model for the 1999 Chi-Chi Earthquake from Inversion of Teleseismic Data Using the Hybrid Homomorphic Deconvolution Method

This study investigates the kinematics of the rupture process of the M _L 7.3 Chi–Chi, Taiwan, earthquake on September 21, 1999. By applying the proposed hybrid homomorphic deconvolution method to deconvolve teleseismic broadband P-wave displacement recordings of the earthquake, this study derives the apparent source time functions (ASTFs) at ten stations located around the epicenter. To further characterize the fault, the kinematic history of the rupture was inverted from ASTFs using a genetic algorithm, coupled with nonlinear iterative technique. The calculated ASFTs reveal that the total rupture event lasted for approximately 27 s. Static slip distribution images indicate that most slip occurred at shallower portions of the fault plane, especially 20–55 km north of the epicenter. The maximum slip reached 20 m at 45 km north of the epicenter, and the average slip throughout the observed rupture area was approximately 2 m. Large asperities on the fault appeared at 25–35 km and 40–50 km north of the hypocenter, and coincided with relatively high rupture velocity. This suggests that the earthquake’s energy may have been released quickly. The rupture velocity decreased upon encountering an asperity, and increased again after passing the asperity. This implies that the rupture required more time to overcome the resistances of the asperities. The maximum rupture velocity was 3.8 km/s, while the average rupture velocity was approximately 2.2 km/s. The rise time distribution suggests that larger slip amplitudes generally correspond to shorter rise times on the subfaults.     

智利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年伊拉佩尔地震的发生,对后者的发生起到了促进作用.     

2017年九寨沟MS7.0地震对周围地区的静态应力影响

基于九寨沟M_S7.0地震的破裂模型及均匀弹性半空间模型,本文计算了该地震在周围主要活动断层上产生的库仑应力变化、在周围地区产生的应力场和位移场和同震库仑应力变化对余震的触发.结果表明：（1）九寨沟地震造成虎牙断裂中段库仑破裂应力有较大增加,已经超过0.01 MPa的阈值,虎牙断裂北段、塔藏断裂中段和岷江断裂北段北部的库仑破裂应力有较大降低,因此尤其要注意虎牙断裂中段的危险性.（2）水平面应力场在该地震震中东西两侧增加（拉张）,张应力起主要作用.在震中南北两侧降低（压缩）,压应力起主要作用.从水平主压和主张应力方向来看,均呈现出条形磁铁的磁场形态.从剖面上的应力场来看,在上盘的面膨胀区域内,大部分点的主张应力方向与地表是垂直的,在其他区域内,主张应力和主压应力均以震中为中心,向外呈辐射状.（3）从地表水平位移场来看,震中东西两侧物质朝震中位置汇聚,南北两侧物质向外流出,在震中处的最大水平位移量达43 mm.从地表垂直位移场来看,震中南北两侧出现明显的隆升,隆升最大值达56.8 mm.震中东西两侧出现明显的沉降,沉降最大值达74.5 mm.从剖面的位移场来看,九寨沟地震为左旋走滑地震,且有一定的正断成分.由分析可以推测该断层破裂在大致22~26 km的深度上就截止了.并推测下盘物质运动的动力来自震源北东东方向（四川块体）深度在6~30 km的下盘下层物质,上盘物质运动的动力来自震源北西西方向（巴颜喀拉块体）深度在0~6 km的上盘上层物质.（4）通过计算不同深度上主震对余震的触发作用可知,主震后的最大余震受到了主震的触发作用,多数其他余震也受到主震的触发作用.主震的发生促使了库仑应力增加地区余震的发生,抑制了一部分库仑应力减少地区余震的发生.     

基于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的库仑应力增加,对此地震有一定的触发作用.     

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%.结合震源位置、震源机制、破裂速度以及辐射效率,我们认为这次地震是由于俯冲板块前缘受到下部地幔物质上浮阻力引起的剪切失稳所致,起初板块内部的脆性破裂表现突出,致使辐射效率较高,后来震源处高温高压下的熔融耗散特征逐渐凸现,致使辐射效率下降.     

Rupture imaging of the 25 April 2015 MW7.9 Nepal earthquake from back-projection of teleseismic P waves

The M_W7.9 Nepal earthquake of 25 April 2015 had over 8, 500 fatalities and was the most destructive earthquake in Nepal since the Bihar-Nepal earthquake in 1934. In this study, we imaged the rupture process of this Nepal event by back-projecting the teleseismic P-wave energy recorded at the three regional networks in Alaska, Australia and Europe. The back-projection images of the three subarrays revealed that the Nepal earthquake propagated along the strike in a southeast direction over a distance of ~ 160–170 km with the duration of ~ 50–55 s. The rupture process was found to be a simple, unilateral event with a near constant velocity of 3.3 km/s. The beam power was mainly distributed in the geographic region just north of Kathmandu and the peak intensity for the source time function curve occurred at about 30 s. The earthquake was destructive due to its occurrence at shallow depth (~ 12–15 km) and the fact that the capital lies in a basin of soft sediment. Additionally, the resonance effect for the longer period waves that occurred in the Kathmandu valley led to destructive aggravation, impacting mainly the taller buildings.