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

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

2015年智利Illapel M_W8.3地震同震效应及其对南美大陆地震危险性影响

基于高性能并行有限元方法,建立含地表地形和Moho面起伏的大规模非均匀椭球地球模型,计算了2015年9月17日智利Illapel M_W8.3地震同震效应,同时讨论2010年智利Maule M_W8.8地震和2014年智利Iquique M_W8.1地震对2015年智利Illapel地震的加载作用,进一步分析这3个地震共同作用对南美大陆及周围断层的地震危险性影响.结果表明:对此次智利Illapel地震CEA(中国地震局地球物理研究所)滑动模型计算的地表最大水平位移约为3 m,USGS(美国地调局)模型的最大水平位移约6 m.特大地震影响范围广,南美大陆接近一半区域同震水平位移量级达到0.5 mm;受3次地震共同作用,智利2010年Maule地震南部俯冲带、2015年Illapel地震和2014年Iquique地震之间俯冲带,以及2014年Iquique地震以北俯冲带,受到10 kPa东西向压应力加载;2010年Maule地震和2014年Iquique地震共同加速了2015年Illapel地震滑动;1960年Valdivia M_W9.5地震破裂区、2014年Iquique地震和1995年Antofagasta地震之间区域以及2010年Maule地震和2015年Illapel地震的未破裂区重合处库仑应力变化均超过10 kPa,有可能加速强震的发生.     

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.静态库仑应力结果显示本次地震增强了震源破裂区周边活动断层的库仑应力,并触发了大量余震,余震主要发生在库仑应力增强区域.结合震间闭锁分布、历史地震及库仑应力变化,我们认为未来需要密切关注与磨西断裂交接的安宁河、大凉山断裂以及康定—磨西段的地震危险性.     

2020年西藏尼玛Mw6.3地震InSAR同震形变特征与破裂滑动分布

2020年7月22日(UTC)西藏尼玛发生Mw6.3地震.为准确确定地震的发震断层、反演破裂滑动分布以及评估地震危险性,本文利用合成孔径雷达差分干涉(DInSAR)技术对Sentinel-1A/B数据进行处理,获取视线向(line of sight,LOS)同震形变场,形变场显示以沉降为主,长轴呈NE向,升、降轨最大沉降量分别为～30 cm和～26 cm.采用基于序列蒙特卡罗采样的贝叶斯方法反演发震断层的位置和几何参数,并在此基础上通过最速下降法反演同震破裂滑动分布.反演结果显示:发震断层走向～30.41°±1.25°,倾角～49.52°±1.06°,属东倾断层.震中位置为(86.89°E,33.18°N),对应深度7.23 km.同震破裂只存在一个滑动中心,主要集中在3.4～11.8 km深度,最大滑动量～0.98 m,平均滑动角-74.54°,表明本次地震是以正断为主兼具少量左旋走滑,释放地震矩约为3.86×1018 N·m,对应矩震级Mw6.36,略大于USGS(United States Geological Survey,美国地质调查局)提供的震级.对于深度5 km和10 km位于破裂区南北端,以及深度15 km和20 km位于破裂区东西侧,其库仑破裂应力变化≥0.01 MPa,未来危险性值得更多的关注.     

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

2015年尼泊尔Gorkha M_W7.9地震与Kodari M_W7.3地震InSAR数据反演及其应力触发分析

根据喜马拉雅断裂系的构造形态,采用缓倾角反铲型断层模型模拟MHT上地震破裂部分的坡坪式发震构造。利用Alos-2及Sentinel-1获取的InSAR数据,反演获得了2015年尼泊尔Gorkha地震及其最大余震Kodari地震的同震滑动分布模型。与单独利用Alos-2或Sentinel-1 InSAR数据的反演结果相比,利用Alos-2和Sentinel-1 InSAR数据联合反演能够提供Gorkha地震破裂的更多细节信息,尤其对深部信息的约束更加明显。联合反演得到的破裂深度最大可达24km,穿过了该区域的闭锁线,到达了闭锁和蠕滑的转换区域。反演的断层模型倾角在3°~10°之间,最大滑动量出现在地下17km处,约4.5m。Gorkha地震和Kodari地震发震性质相似,都是发生在MHT断层上的低角度逆冲型地震,其中Gorkha地震略带右旋分量。反演结果还显示,Gorkha地震与Kodari地震的破裂滑动在空间上存在互补性,Kodari地震就发生在Gorkha地震的破裂空区内。通过计算Gorkha地震对Kodari地震发震断层的库仑破裂应力加载,发现Kodari地震震中恰位于库伦破裂应力正负交界区域,库仑破裂应力加载达0.4MPa,表明Kodari地震可能受到了Gorkha地震的触发。     

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后破裂基本趋于停止.推断鲁甸地震破裂在上地壳浅层集中释放了大多数能量是导致灾害严重的主要原因之一.     

2022年青海门源MW6.6地震InSAR同震形变场与震源特征

本文利用合成孔径雷达差分干涉(Differential Interferometric Synthetic Aperture Radar, D-InSAR)技术处理哨兵卫星升降轨影像,获取了2022年青海门源M_S6.9地震的同震变形场,变形结果显示此次地震导致冷龙岭断裂西段两侧约30 km×20 km范围明显变形,卫星视线向最大同震位移达50 cm.基于有限断层模型,本文构造断层滑动分布模型,反演同震变形场数据.模型显示,冷龙岭断裂的破裂以左旋走滑为主,破裂长度约23 km,主要分布在10 km以上深度,整体形态表现为不规则椭圆盘,其中约3.9 m的最大破裂位于5～6 km深度的椭圆中心.基于反演模型估算的地震平均应力降约5.9 MPa,地震标量矩为1.0×10¹⁹ N·m,对应矩震级为M_W6.6.本文计算了1927年古浪地震、1986年及2016年门源地震对本次地震破裂面的同震静态库仑应力扰动,结果表明冷龙岭断裂西段受到历次地震的库仑应力加载作用,库仑应力增加累积达0.12 MPa,对本次发震断层的破裂有促进作用,其中...     

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 ...     

基于GNSS和InSAR约束的2021年玛多MS7.4地震同震滑动分布及应用

2021年5月22日青海省果洛藏族自治州玛多县发生M_S7.4地震,震中位于青藏高原中部的巴颜喀拉块体,这是近20多年来在巴颜喀拉块体周边发生8次M≥7级强震后,块体内部的一次强震,也是汶川地震以来中国大陆发生的最大一次地震,因此该地震的成因及周边地区未来的地震危险性值得重点关注.本文利用震后及时获取的39个近场流动GNSS观测,联合61个GNSS连续观测、Sentinel-1和ALOS-2 InSAR观测获取了本次地震精细的同震形变场,以此为约束,基于均匀弹性半无限位错模型,反演了发震断层的滑动分布,并计算了同震库仑应力变化.GNSS水平同震形变十分显著,断层南北两侧的GNSS点位,最大水平形变分别达0.7 m和-1.2 m,距震中200 km的测点仍有1 cm左右的同震形变.Sentinel-1和ALOS-2的升降轨InSAR同震形变场显示此次地震造成了约160 km长的地表破裂,最大视线向形变分别达0.9 m和1.2 m.同震滑动分布模型显示,发震断层由主段和次段组成,长度分别为170 km和20 km,主段倾向北,倾角85°,平均滑动角为-4.36°,表明玛多地震是一次典型的走滑型地震.次段倾向南,倾角68°,平均滑动角为-11.84°.地震破裂主要集中在0~15 km深度范围,最大滑动量为4.4 m,对应深度6.97 km.反演给出的矩震量为1.61×10²⁰N·m,对应矩震级M_W7.4.主发震断层上存在4个凹凸体,玛多地震是一次不对称双侧破裂事件.结合余震精定位、野外调查及地质资料,我们认为主发震断裂为昆仑山口—江错断裂,东部的次级破裂与主破裂机制不同.同震库仑应力结果显示,东昆仑断裂玛沁段应力有所增加(>0.01 MPa),处于应力加载状态,未来发生强震的危险性较高.     

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

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°.这些在震源机制、能量释放以及持续时间方面的变化都是当地构造和应力环境复杂性的反映.     

2020年6月23日墨西哥MW7.4地震震源特征

2020年6月23日15时29分04秒（UTC）,在墨西哥南部瓦哈卡州发生了一次震级为M_W7.4的地震,我们利用全球地震台网（GSN）和国际数字地震台网联盟（FDSN）台网的长周期和宽频带P波数据反演分析了这次地震的震源机制、震源时间函数以及时空破裂过程.根据反演结果,这次地震的矩心震中位于15.96°N,95.89°W,矩心深度约为22 km;地震持续15 s左右,释放地震矩1.24×10²⁰ N·m,相当于矩震级M_W7.4;破裂过程比较简单,仅有一个走向和倾向方向尺度相当的凹凸体错动,最大位错达8.1 m,位于21 km深处.凹凸体破裂主要沿断层的滑动方向呈双侧破裂,两个优势破裂方向在地表投影的方位分别位于60°和270°左右.综合构造背景、震源位置、余震分布、震源机制以及时空破裂过程,我们相信这次地震是发生在北美大陆板块和太平洋海底板块相互作用的结果.海底板块朝着大约60°左右的方位运动,以大约22°的倾角插入大陆板块,造成一个凹凸体错动,形成了这次地震.     

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甲仙地震也同该断层有关.     

2018年印尼帕卢MW7.5地震——一次超剪切破裂事件

利用反投影方法,使用日本密集台网Hi-net远场垂直分量568条P波资料对2018年9月28日印尼帕卢M_W7.5地震震源破裂过程进行成像,结果显示此次地震的能量释放比较集中,主要集中在10~20 s之间.破裂有两个集中区,破裂峰值分别位于12 s和19 s,最大能量释放区域位于震中南侧约0~50 km内,另一破裂集中区覆盖了帕卢市及周边区域.破裂主要向南侧延展,破裂总长度至少100 km,平均破裂速度约4.1 km·s^-1,属于一次超剪切破裂事件.     

2008年10月西藏当雄MW6.3地震震后形变提取与余滑反演

2008年10月6日西藏当雄发生M_W6.3地震.本文利用震后2008年10月26日至2010年8月22日的16期ENVISAT ASAR数据,通过小基线集干涉测量、误差校正与MInTS（Multiscale InSAR Time Series）技术提取高精度的震后形变场,利用SDM（Steepest Descent Method）方法反演断层震后余滑演化过程,并分析震后余滑与同震滑动的关系.结果表明：当雄M_W6.3地震的近场震后形变场主要位于断层西侧,在时间演化上具有明显的对数函数衰减规律;震后余滑主要集中于断层中南段深0~15 km区间,最大的余滑量约0.07 m,位于断层深约9.28 km处,滑动角约-103°;震后余滑引起的地震矩能量M₀与矩震级M_W在时间演化上具有指数函数递增规律;当剪切模量μ=32 GPa,震后665天余滑释放的地震矩能量约为1.92×10¹⁷N·m,约占同震滑动释放地震矩的4.8%,相当于矩震级M_W5.46;虽然震后余滑已经延伸到断层浅部0~5 km区间,但由于余滑量相对较小,没有改变同震滑动在断层浅部区域的滑动亏损现象,这可能是2010年11月30日该区域又发生M_W5.3级余震的主要原因之一.     

基于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.通过大柴旦地震发震断层和破裂机制综合分析,初步判断发震断层为西南倾向的大柴旦—宗务隆山断裂.     

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

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