2013年4月20日四川芦山MW6.7 (MS7.0)地震参数的测定

2013年4月20日四川芦山M_W6.7(M_S7.0)地震发生后, 中国地震台网中心(CENC)发布了地震速报参数. 该文利用中国国家地震台网97个台站的资料对地震速报参数进行了修订, 得出: 四川芦山M_W6.7地震的发震时刻为北京时间8时2分47.5秒(世界时间0时2分47.5秒), 震中位置为30.30°N、 102.99°E, 震源深度17 km. 该地震的面波震级为M_S7.0, 短周期体波震级为m_b6.0, 中长周期体波震级为m_B7.0; 利用波形反演的方法计算了震源机制解, 得到的最佳双力偶解的参数分别为节面Ⅰ: 走向17°/倾角48°/滑动角80°; 节面Ⅱ: 走向212°/倾角43°/滑动角101°, 矩震级为M_W6.7. 中国地震台网中心发布本次地震为面波震级M_S7.0, 而美国地质调查局(USGS)国家地震信息中心(NEIC)发布为矩震级M_W6.6. 为了消除这种差别, 建议我国也应将矩震级作为对外发布的首选震级, 使震级的发布与国际接轨.      

2016年日本本州东岸近海MS7.2地震前后的地壳运动

2016年11月22日日本本州东岸近海发生东日本M_W9.0大地震的M_S7.2强余震。本文利用美国内华达大学内华达大地测量实验室网站获得此次大地震周围共30多个GPS连续观测站及其它台站的IGS08全球参考框架坐标时间序列,采用北京附近的GPS连续观测站BJSH作为区域位移参考框架的核心站,获取了此次M_S7.2强余震的同震水平位移和区域参考框架位移时间序列,得到了此次强余震前后的位移时空变化图像。结果显示:尽管本州东岸近海M_S7.2地震作为东日本M_W9.0地震的强余震,受其震后形变的影响强烈,但其地壳水平形变的前兆规律与已观测到的大地震一致;不同的是东日本M_W9.0地震前的垂直位移无积累,而本州东岸近海M_S7.2地震前后的垂直位移保持M_W9.0地震后均匀而缓慢的衰减变化;临震前震中附近的GPS连续观测站的东西分量明显减速,有的甚至减至零,是明显的短期前兆异常。此外,本文还进一步讨论了两种同震位移及其成因,并推断地壳水平运动挤压是此次地震的成因。      

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

利用高频GPS资料研究2016年新西兰凯库拉地震的地表形变及预警震级

2016年11月13日,新西兰凯库拉地区发生了M_W7.8级地震.本文利用1 Hz高频GPS观测数据,基于GAMIT track解算模块,采用主分量分析（PCA）空间滤波方法获取了地震地表形变.同时从trackRTr模拟实时解算的动态位移中,提取P波5秒峰值位移（P_d）和地面峰值位移（PGD）,并根据震级统计回归模型计算预警震级.结果表明：测站的动态形变时间长达2 min,距震源最近的HANM和KAIK站出现二次剧烈形变,震源北部测站的形变幅度大于南部,而高频GPS静态同震形变场表现出先逆冲后走滑的震源机制特征.不同GPS台站的P_d预警震级相差较大,最大震级差为M_W2.5.综合考虑预警震级发布的时效性和可靠性,采用顾及空间分布的四台站PGD联合预警方法,其预警震级在震后23 s达到初始稳定（M_W7.56）,在震后110 s达到最终稳定（M_W7.78）,该震级与USGS矩张量反演震级（M_W7.8）基本一致.     

青海玛多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以上地震需要加强防范.     

2020年阿拉斯加两次MW>7.5地震的应力触发关系研究

地震应力触发理论认为,相邻的大地震之间存在应力触发关系。2020年在阿拉斯加俯冲带发生了两次M_W>7.5地震,分别为7月22日M_W7.8地震和10月19日M_W7.6地震,为研究“双震型”大地震的触发作用提供了机会。利用近场的GPS同震位移和远场的地震波形资料,采用不同的断层参数,反演了M_W7.8地震的4个滑动破裂模型,并结合初始破裂模型,计算了在M_W7.6地震震中位置引起的库仑应力变化。研究结果显示,在5个模型中有4个计算出的应力变化都远远超过了触发阈值0.01 MPa,表明前一个地震对后一个地震有明显的触发作用。根据M_W7.6地震发震断层两侧的闭锁程度变化,以及该区域板块俯冲速率,估算M_W7.8地震使M_W7.6地震提前32年发生。     

2020年西藏尼玛MS6.6和MS4.8地震震源机制测定

北京时间2020年7月23日04时07分,西藏自治区那曲市尼玛县发生M_S6.6地震,震源深度10 km,震中位置为（33.19°N,86.81°E）。主震发生当日18时50分,发生一次M_S4.8强余震,震源深度为10 km。本文基于西藏、青海、新疆区域波形资料,采用ISOLA近震全波形方法对这两次地震进行震源机制反演。结果显示,尼玛M_S6.6主震的最佳断层面解为:节面Ⅰ走向8°/倾角46°/滑动角?93°,节面Ⅱ走向191°/倾角44°/滑动角?87°;矩震级M_W6.4,最佳矩心深度7 km。震源区应力主轴的空间取向为:主压力轴P的方位角220°、倾伏角88°,主张力轴T方位角99°、倾伏角1°。M_S4.8强余震的最佳断层面解为:节面Ⅰ走向12°/倾角47°/滑动角?106°,节面Ⅱ走向214°/倾角45°/滑动角?74°;矩震级M_W5.0,最佳矩心深度6 km。震源区应力主轴的空间取向为:主压力轴P的方位角207°、倾伏角78°,主张力轴T方位角113°、倾伏角1°。震源机制反演结果表明,这两次地震均为以正断型为主的地震事件,与震源区附近先前地震的震源机制有较好的一致性。结合周边地质构造和余震分布,我们认为尼玛M_S6.6地震可能是由位于日干配错断裂和依布茶卡盆地西缘断裂之间的一条正断层活动所引发的。      

GPS观测得到的大地震前兆地壳形变震例

综合介绍2008年汶川大地震以来,GPS观测得到的国内外10多次6—9级,不同构造、不同类型的大地震前兆地壳形变震例:2008年汶川8级大地震、2011年东日本9级巨震、2013年芦山7级,直至2020年6月墨西哥7.4级地震和7月美国阿拉斯加州以南海域7.8级地震等。利用GPS连续观测站区域参考框架水平位移时间序列和水平位移场,特别是水平位移向量时间序列的研究证明,同震水平位移是研究地震前兆形变存在的关键;利用垂直位移和水平位移向量时间序列、同震垂直位移及同震水平位移向量的分解,揭示地震弹性回跳真实方式;提出了符合GPS观测和岩石破裂试验结果的地震压-剪弹性回跳模型;根据已有震例,提出预报不同震级地震的可能性和监测临震前兆形变的GNSS站布设设想。      

基于余震重定位和震源机制解研究青海玛多MS7.4地震序列的发震构造和断裂形态

2021年5月22日02时04分(北京时间),青海果洛州玛多县发生M_S7.4地震,震后余震不断.地质调查和卫星观测对地表断裂痕迹有较好的约束.然而,对于理解区域应力场、地震的产生、传播和终止具有重要意义的地下断层几何结构的约束精度略显不足.利用国家地震台网的连续波形记录,本研究首先基于双差定位法对玛多地震震后25天的余震序列进行重定位,结果显示余震序列大致沿NWW向的江错断裂呈线性分布,位于主震震中两侧,延伸总长~170 km.主震东南侧存在一余震稀疏区,在断裂带东西两端余震分布转向且出现分叉现象,反映出发震断层的复杂几何形态,这与前人研究结果基本一致.进一步采用波形反演方法和P波初动极性反演方法,获得了玛多震源区132个中小余震的震源机制解与震源矩心深度,并基于此对该主余型地震的发震构造与断裂形态进行了初步分析.震源机制解结果表明,玛多M_S7.4主震的发震断裂主要为左旋走滑性质,余震与主震性质整体相同,在断裂带东段存在部分逆冲型余震.震源机制解约束的区域主应力方向约N60°E,与区内整体走滑断裂作用相一致.余震震源深度略微起伏,主要集中在10~12 km,且浅部余震较少,表明浅部应力可能主要通过主震释放,余震深度分布可能限定了主震同震破裂的下边界.玛多主震破裂起始于断裂带走向和倾向发生明显变化的位置,表明断裂带的复杂几何结构可能是此次玛多M_S7.4地震初始破裂空间分布的决定因素.主震破裂结束的两端都有"马尾状"构造(或次级断层),表明这种分叉断层复杂的几何形态可能控制着主震破裂的最终位置.     

2021年5月21日云南漾濞地震及5月22日青海玛多地震的支持向量机震级估算

2021年5月21日及5月22日,云南漾濞县与青海玛多县分别发生破坏性地震,主震震级分别为M_s6.4级与M_s7.4级。本文基于机器学习中的支持向量机方法,以多类型特征参数为输入建立地震预警震级估算模型SVM-M,离线模拟云南漾濞M_s5.6级前震、M_s6.4级主震以及青海玛多M_s7.4级主震的连续震级估算。结果表明:对于云南漾濞M_s5.6级前震,支持向量机方法在首台触发后1s可估算震级为5.6级,且随着首台触发时间的增加,估算震级一直在实际震级附近波动;对于云南漾濞M_s6.4级主震和青海玛多M_s7.4级主震,随着首台触发时间的增加,支持向量机方法对于大震低估问题得到了有效的改善,且震级估算结果逐渐接近实际震级。同时,这3次地震的震级估算离线模拟表明:引入震源距的支持向量机方法（SVM-M1模型）对于震级估算有更好的稳定性,且在地震预警系统的震级估算中有着潜在的应用前景。     

日本东北MW9.0地震的PI模型参数设置与预测效能回溯性检验

以日本局部地区（32.0°N—46.0°N,136.0°E—148.0°E）为研究区域,应用图像信息（PI）方法,获得了不同计算参数模型下包含2011年3月11日日本东北M_W9.0地震的多组预测窗热点图像。以0.5°×0.5°和1.0°×1.0°的网格尺度和5—10年预测窗长为主要参数变量,并以R值和受试者工作特征（ROC）方法检验不同参数模型下PI方法的预测效能。结果表明,多组参数模型下M_W9.0地震所在预测窗内,其震中所在网格及其摩尔邻近网格均曾出现热点图像,表明PI方法可对日本东北M_W9.0地震作出预测。综合R值评分和ROC检验分析可知,网格尺度相对较大、预测窗长相对较长的模型,其预测效果更好。      

2016年青海门源MS6.4地震震前应变积累及同震变形特征

2016年1月21日青海省门源县发生了M_S6.4地震, 发震断裂为冷龙岭北侧断裂, 震中位置与1986年门源6.4级地震相同。 本文收集了本次地震震中及其周边区域1999—2015年GPS观测资料, 解算了GPS速度场、 跨断裂连续观测站基线时间序列和应变率场。 结果显示, 祁连山断裂带为一条宽度约60 km的连续变形带。 在断裂带南侧地壳运动以顺时针旋转为主, 运动量值没有显著差异; 跨过断裂带到达其北部之后, 地壳运动量值明显减小, 显示出该断裂带的强烈活动特征。 冷龙岭断裂左旋走滑速率为6.17±0.41 mm/a, 挤压速率为1.83±0.38 mm/a, 断层闭锁深度为22.1±3.1 km。 利用GPS连续观测站数据解算的地震同震位移显示, 震中西南侧26.8 km处的青海门源(QHME)测站记录到了明显同震位移, 其水平运动方向为北东向, 与逆冲为主的震源机制解一致。     

基于无人机倾斜航空摄影三维点云测量同震倾滑变形研究

2021年5月22日青海玛多发生M_S7.4地震,震源断层错动在地表形成了长达160 km的同震地表破裂。可靠的地震地表破裂带参数是研究震源断层活动机制和评价地震危险性的重要基础。采用无人机倾斜摄影测量技术可以获得高精度的点云数据并产出DOM和DEM数据。通过跨破裂带的地形测量,获取了玛多M_S7.4地震同震地表变形的垂直位移、水平缩短量和水平拉张量等参数。测量结果显示,玛多M_S7.4地震发震断层在不同破裂段具有不同性质和大小的倾滑分量,其中具有压扭性质的野马滩观测点断层垂直位移为0.69～1.01 m,倾向水平缩短量为0.17～0.41 m,倾滑位移为0.71～1.09 m;具有张扭性质的朗玛加合日段断层垂直位移为0.34～0.54 m,倾向水平拉张量为1.99～2.08 m。     

玛曲MS4.4地震震源机制解及其构造意义探讨

2021年5月22日青海省玛多县发生M_S7.4地震,数小时后在距离震中两百多公里的甘肃玛曲县发生了M_S4.4地震。利用甘肃、青海和四川区域测震台网记录的三分向宽频带数字波形资料,反演甘肃玛曲M_S4.4地震的震源机制解,结果显示此次地震活动面走向、倾角和滑动角分别为105.6°、74.1°和-38.7°。参考玛多M_S7.4地震的震源机制解,发现两次地震震源机制解具有较好的一致性,均呈现明显的左旋走滑特征。静态库伦破裂应力改变量分布计算结果表明,玛曲M_S4.4地震震中位置单位面积(m~2)受到来自玛多地震震中方向的拉应力约为0.02 MPa。综合两次地震的震中距、发震时刻和断层分布等情况,初步判断甘肃玛曲M_S4.4地震应为青海玛多M_S7.4强震的一次触发地震。     

基于频谱分析的矩震级测定方法——以山西地震带为例

矩震级M_W是目前量度地震大小最理想的物理量,已成为国际地震学界开展研究的首选震级。对于中强地震事件,矩震级M_W通常可以利用长周期地震波,基于矩张量反演计算获得;对于约3级以下的地震事件,直接进行矩张量反演有时会产生较大误差。本文以山西地震带为例,利用山西地震台网数字波形资料,开展研究区Q值函数模型的分析评估与基于频谱分析的矩震级测定。研究结果表明,基于Q值函数Q(f)=299.4f^0.563计算的频谱M_W更接近矩张量反演的M_W,与Q值模型参数Q₀、α和κ的优化组合结果存在较好吻合;对于3.6≤M_L≤5.1地震事件,基于频谱分析计算得到的矩震级与基于矩张量反演矩震级较接近。在已具有详细的区域非弹性衰减模型基础上,基于频谱分析快速测定中小地震的矩震级具有可行性。     

2021年玛多MS 7.4地震前MS 3、MS 4地震空区嵌套异常分析

利用地震活动图像等传统测震学方法,对玛多M_S7.4地震前的M_S3、M_S4地震空区嵌套异常演化过程作进一步梳理分析,并进行了回顾总结,得出如下认识和结论：(1)该地震空区具有孕育强震的能力且出现视应力高值现象,预示着区内有发震危险性;(2)未来主震震中出现在空区嵌套的交界区附近,这一特征对于认识空区嵌套异常、缩小预测发震区域极其重要;(3)回顾玛多地震前空区演化过程,有2次显著地震事件,使跟踪从中长期转为中短期再转为短期,这可为空区跟踪过程中发震紧迫性的判定不断积累经验。     

云南大理MS6.4地震与意大利拉奎拉MW6.3地震序列的比较研究和危险性分析启示

云南大理M_S6.4地震和意大利拉奎拉（L'Aquila）M_W6.3地震都因未能准确做出危险性预测给人民生命财产造成重大损失。哪些工作可能向政府和公众说明“哪里可能发生地震”、“最大量级多大”、“未来发震趋势如何”等问题都值得深入思考。本文根据天气预报思路,将地震活动按丛集性做分区处理,从各区M-t序列、3级地震活动性、2年来地震能量释放升级趋势、G-R关系等几方面分析了这两次地震主震发生前的地震活动特征,对比分析了它们之间的相似性。认为这两次地震的相似性有以下几点:①震源深度都比较浅。云南大理M_S6.4地震震源深度8 km,意大利拉奎拉M_W6.3地震震源深度8.8 km;②b值相近。云南大理M_S6.4地震b=0.59,意大利拉奎拉M_W6.3地震b=0.61。G-R关系外推震级相近,均为M_GR6左右;③M-t图序列地震能量释放都呈升级趋势,地震活动也相似。平静打破后3级以上地震活跃,都出现震群现象,地震能量释放呈加速状态;④都属于前震-主震-余震型序列。意大利拉奎拉M_W6.3地震主震前最大前震为M_W5.2;云南大理M_S6.4地震主震前最大前震为M_S5.6,且前震震群特征明显。同时,本文讨论了按地震丛集性划分区域的合理性,认为按地震丛集性划分区域更容易把握区域地震活动的特点。这两次地震发生在不同区域,构造差异极大,属不同错动类型,用分区G-R关系和M-t序列分析,都得出较好的结果,再次说明了分区G-R关系和M-t序列分析对判断地震危险性的普适性特征,这为下一步产出“地震危险性云图”提供了依据。      

中强地震能量震级测定

本文根据地震波衰减特性,开展了利用宽频带地震波形数据测定地震波能量E_S和能量震级M_e的方法研究。利用震中距处于20°—98°范围内的宽频带远震P波波形数据,测定了4次国外和4次国内中强震的能量震级M_e,并对其面波震级M_S、矩震级M_W及能量震级M_e进行了分析对比。结果表明:面波震级M_S表示的是地震在某一固定频率所辐射的地震波能量大小;矩震级M_W与地震所产生的断层长度、断层宽度、震源破裂的平均位错量等静态构造效应密切相关;而能量震级M_e反映的是震源动态特征,与地震震源的动力学特性密切相关。由于地震是以地震波形式辐射,能量主要集中在震源谱的拐角频率附近,因此能量震级M_e更适合描述地震的破坏性。由此可见,联合测定面波震级M_S,矩震级M_W和能量震级M_e对于地震定量研究以及地震灾害与风险评估具有重要作用。      

新疆中强地震矩震级的初步研究

整理新疆地震台网2010年1月至2023年2月M_L≥4.0的623个中强地震波形数据，利用CAP方法反演得到相应矩震级，使用回归分析得到矩震级M_W和近震震级M_L之间的关系式为M_W=0.668+0.789M_L,相关系数为0.841,标准差0.272 4。回归的标准化残差近似正态分布，矩震级M_W和近震震级M_L呈正相关关系，表明在新疆地震台网日常工作中使用CAP方法测定中强地震的矩震级是可行的。     

2021年5月22日玛多MS7.4地震震前岩石圈磁场变化特征分析

本文利用2019年和2020年两期的全国流动地磁矢量数据,针对玛多M_S7.4地震震中附近的岩石圈磁场空间分布,分析研究了岩石圈磁场各个分量的震前变化特征。结果表明,玛多M_S7.4地震发生前,震中附近的岩石圈磁场分量均发生了不同程度的变化:震中位于各分量的弱变化区域和零变线附近;震中位于各分量的高梯度带和低梯度带之间。本文证实了地震发生前会引起岩石圈磁场变化这一现象,并总结了震前地磁分量的变化特征,为今后强地震（尤其是M≥7.0地震）的震磁关系研究提供了一个震例参考。