2016年8月24日意大利佩鲁贾MW 6.2级地震同震位移场估计和震源滑动模型反演

利用改进的自动经验基线校正方法SMBLOC,对2016年8月24日意大利佩鲁贾MW 6.2级地震震中周围约60 km内的近场强震记录进行基线校正并尝试给出同震位移场,与GPS观测结果进行对比分析,分别独立和联合两种资料反演震源滑动模型,并根据震源模型进一步给出全空间预测位移场分布.研究结果表明:(1)两种不同的资料给出的水平位移场幅值均为cm级,且均表明断层的错动以正断为主.(2)两种同震位移场分别独立和联合反演所得的震源静态滑动范围基本一致,最大滑动均发生在震中东北侧,强震模型表现出明显的双事件特征,较大滑动分布在震中东北侧和东南侧,GPS模型在震中东南侧的滑动相对较小,其双事件特征不明显.两种模型的最大滑动量分别为0.96 m和0.86 m,较为一致,反演的矩震级均在MW 6.3左右.(3)根据震源滑动模型计算所得的佩鲁贾地震全空间预测的水平同震位移场中最大位移分布区域与震后报告中受灾严重的地区基本一致.表明在一定的条件下,利用SMBLOC方法解算震级较小的MW 6.0左右地震强震记录的同震位移场,并反演震源滑动模型具有一定的可行性,且其同震位移场和滑动模型可为震后灾害快速评估、救援力量分配、余震趋势判定等快速应急响应工作提供参考依据.     

芦山7级地震的同震位移估计和震源滑动模型反演尝试

利用自动经验基线校正方法,分析2013年4月20日芦山M_S7.0级地震13个近场强震动台的观测资料,以估算同震位移场分布,并据此反演了震源滑动模型.经与GPS结果比较,两种不同方法给出的芦山7级地震的水平近场同震位移场幅度都不超过cm级,均显示为典型的逆冲型地震（兼有少量左旋走滑错动）.强震最大水平和垂直永久位移分别为4.9 cm和4.4 cm,分别出现在51YAM台和51QLY 台.两种资料反演的震源滑动模型虽显示多事件特征,但主要滑动均集中在第一次事件,即初始滑动点两侧的走向长约30 km、倾向长约25 km的相对集中的较小范围内,强震和GPS模型的最大滑动量分别为1.14 m和1.09 m,较为一致.其余子事件滑动量小且分布零散,不能排除其数值效应的因素.反演矩震级均在M_w6.7左右,地表破裂应该不明显.文章还讨论了目前在我国利用近场强震动记录估计M_w6~7级地震同震位移场存在的困难和问题,为今后类似工作提供参考.     

芦山7级地震的同震位移估计和震源滑动模型反演尝试

利用自动经验基线校正方法,分析2013年4月20日芦山MS7.0级地震13个近场强震动台的观测资料,以估算同震位移场分布,并据此反演了震源滑动模型.经与GPS结果比较,两种不同方法给出的芦山7级地震的水平近场同震位移场幅度都不超过cm级,均显示为典型的逆冲型地震(兼有少量左旋走滑错动).强震最大水平和垂直永久位移分别为4.9cm和4.4cm,分别出现在51YAM台和51QLY台.两种资料反演的震源滑动模型虽显示多事件特征,但主要滑动均集中在第一次事件,即初始滑动点两侧的走向长约30km、倾向长约25km的相对集中的较小范围内,强震和GPS模型的最大滑动量分别为1.14m和1.09m,较为一致.其余子事件滑动量小且分布零散,不能排除其数值效应的因素.反演矩震级均在Mw6.7左右,地表破裂应该不明显.文章还讨论了目前在我国利用近场强震动记录估计Mw6~7级地震同震位移场存在的困难和问题,为今后类似工作提供参考.     

用近场强震动记录快速估计同震位移并反演震源滑动分布

利用自动经验基线校正方法,分析日本2008年岩手-宫城内陆M_w6.9地震震中周围密集强震动观测台网资料,快速解算出了同震位移场分布,并据此反演了震源滑动模型.经与GPS结果比较,两种不同方法给出的同震位移幅值、方向和总体分布特征较为接近.基于相同断层面参数反演的震源模型空间展布形态、主要滑动范围、平均和最大滑动量、滑动方向以及由模型计算的矩震级等均吻合较好,从而验证了方法的可行性.讨论了自动经验基线校正方法尚存在的问题和不足,为今后利用强震资料快速解算M_w6-7级及以上地震的同震位移场并反演震源滑动分布提供参考.     

2016年日本熊本地震破裂时空过程联合反演

为了深入认识2016年4月15日日本熊本地震破裂的复杂性,利用远场体波资料和同震InSAR资料联合反演了此次地震的震源破裂时空过程. 联合反演结果表明:熊本地震的震源破裂持续时间约为25 s,整个破裂过程释放的总标量矩为6.03×1019 N·m,对应于矩震级M_W7.1;同震滑动主要集中分布于浅部,破裂以右旋走滑为主,但在沿倾向0—5 km范围内,破裂呈较强的正断特征;此次地震破裂的最大同震滑动量约为4.9 m,且最大同震位错区位于背离断层走向上、距离起始破裂点约5—10 km的区域;破裂前期（0—7 s）,在倾向上向浅表发生破裂,在走向上向东北和西南两侧扩展;大约7 s后,破裂背离断层走向主要向东北方向扩展. 根据有限断层联合反演结果推测,此次熊本地震破裂可能出露至地表.      

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

利用同震和震后位移数据联合反演2011年日本MW9.0地震同震断层滑动

基于黏弹性球体位错理论, 联合陆地和海底同震GPS数据以及日本本岛330个陆地GPS站点5~10年的震后数据, 反演了日本M_W9.0地震的断层滑动模型, 提升了断层滑动分布在细节上的合理性。 首先, 基于日本本岛330个陆地GPS站点震前2年和震后10年的连续观测数据, 获取了日本M_W9.0地震震后5~10年的年平均位移, 该时段的位移几乎完全由地幔黏弹性松弛效应引起; 接着, 利用黏弹性球体位错理论对震后5~10年的位移进行反复拟合, 确定了日本M_W9.0地震震源及周边地区的地幔黏滞性系数最优解(9.0×10¹⁸ Pa·s)。 然后, 联合同震和震后位移数据, 引入黏弹性位错格林函数, 反演了2011年日本M_W9.0地震的断层滑动分布。 结果表明, 该地震同震破裂的最大值达到了62.72 m, 同震滑动的总地震矩为4.48×10²² Nm, 相应的矩震级为M_W9.03。 由于黏弹性松弛效应引起的震后位移中包含了同震破裂的信息, 基于黏弹性球体地震位错理论, 联合同震和震后位移数据反演断层同震破裂, 有效提高了日本M_W9.0地震断层滑动分布的可靠性。 最后, 本文提出的反演方法为同震观测结果缺乏的大地震震后科考提供了理论支撑: 在大地震发生之后, 即使在同震期间没有足够的观测数据, 也可以在震后通过对震源区的加密观测积累的震后数据, 使用本文提出的反演方法优化同震断层滑动模型。     

2020年阿拉斯加Mw7.8地震震源特征及邻区地震危险性分析

2020年7月22日,在美国阿拉斯加半岛南部发生了一次Mw7.8地震.本文利用远场地震波形记录与近场GPS台站同震位移资料,反演了本次地震震源机制和震源破裂过程.结果表明:阿拉斯加地震为一震源深度为23 km,地震断层倾角约17°,破裂面上最大的同震位移达到914 cm的逆冲型地震.由此得到该地震的地震矩为6.94X1...     

2018年5月28日吉林松原MS5.7地震发震构造分析

利用黑龙江、吉林、内蒙古三省数字地震台网的三分量宽频带波形资料对2018年5月28日吉林省松原市宁江区发生的M_S5.7地震进行了全波形矩张量反演,获取了此次地震的震源机制解和矩心位置;并使用震源-矩心方法讨论了该地震的发震断层。研究结果显示:松原M_S5.7地震的矩震级为M_W5.2,矩心位置为（45.225°N,124.685°E）,矩心深度为7 km。震源机制解参数显示:该地震为走滑型;节面Ⅰ的走向、倾角和滑动角分别为217°,82°和164°;节面Ⅱ的走向、倾角和滑动角分别为309°,74°和8°;双力偶成分占96.4%,方差减少为93%。震源?矩心图显示震源更接近节面Ⅰ,与北东向的扶余—肇东断裂走向及倾角一致,因此,推测扶余—肇东断裂为发震断层。      

2021年日本福岛县冲地震的震源破裂过程分析基于采用经验格林函数方法的波形反演

利用中小震作为经验格林函数,选取0.2~2.0 Hz频段的强震数据进行波形反演,获得了2021年福岛县冲地震的破裂过程。结果表明:该地震的破裂主要集中在断层面距离震源约25 km的区域内,沿震源向东北延伸约5 km,向西南延伸约20 km;在该区域内,识别出两个滑动量集中的区域,均分布在距离震源约15 km西南侧,主要滑动量集中区域最大滑动量约3.2 m,深度基本与震源一致;次要滑动量集中区域略比主要滑动量集中区域浅约18 km。该震源模型对应矩震级M_w7.3,破裂速度为2.4 km/s。通过选择不同的中小震组合进行波形反演,结果对该震源模型特性无显著影响,揭示了该震源模型的稳健性;基于该震源模型合成反演分析中未使用的强震观测台站强震动,获得的合成波形与观测波形有很好的相关性,充分证明了该震源模型时空特征的可靠性。     

2014年新疆于田MS7.3地震同震位移及位错反演研究

利用于田震中300 km范围内的1个GPS连续站和12个GPS流动站数据,解算得到了2014年新疆于田M_S7.3地震地表同震位移,并反演了发震断层滑动分布,探讨此次地震对周边断裂的影响.地表同震位移结果显示,GPS观测到的同震位移范围在平行发震断裂带的北东-南西向约210 km,垂直发震断裂带的北西-南东方向约为120 km,同震位移量大于10 mm的测站位于震中距约120 km以内;同震位移特征整体表现为北东-南西方向的左旋走滑和北西-南东方向的拉张特征,其中在北东-南西方向,I069测站位移最大,约为32.1 mm,在北西-南东方向,XJYT测站位移最大,约为28.1 mm;位错反演结果表明,最大滑动位于北纬36.05°,东经82.60°,位于深部约16.6 km,最大错动量为2.75 m,反演震级为M_W7.0,同震错动呈椭圆形分布,以左旋走滑为主并具有正倾滑分量,两者最大比值约为2.5：1,同震错动延伸至地表,并向北东方向延伸,总破裂长度约50 km,地表最大错动约1.0 m;同震水平位移场模拟结果显示贡嘎错断裂、康西瓦断裂和普鲁断裂等不同位置主应变特征具有差异性,这种差异特征是否影响断裂带以及周围区域的应力构造特征,值得关注.     

基于形变观测分析2011年日本9.0级地震与断层运动间关系

2011年3月11日日本发生9.0级地震,本文以此次地震的震间、同震和震后形变观测为约束,依据不同时段断层运动空间分布特征分析日本海沟地区强震与断层运动间关系.震间日本海沟地区,断层运动闭锁线深度约为60km,闭锁线以上从深到浅依次为断层运动强闭锁段、无震滑移段和弱闭锁段.由同震位错反演结果,2011年日本9.0级地震同震存在深浅两个滑移极值区,同震较浅的滑移极值区(同震位错量10~50m,深度小于30km)震间为断层弱闭锁段;同震较深的滑移极值区(同震位错量10~20m,深度在40km左右)震间为断层强闭锁段;而在两者之间的过渡带同震位错相对较小,震间断层运动表现为无震滑移.震后初期断层运动主要分布在在闭锁线以上的同震较深滑移极值区,而同震较浅的滑移极值区能量释放比较彻底,断层震后余滑量相对较小.依据本文同震和震间断层运动反演结果,震间强闭锁段积累10m同震位错需要100多年时间,与该区域历史上7级地震活动复发周期相当;震间弱闭锁段积累30~50m同震位错约需要300~600年时间,与相关研究给出的日本海沟9级左右地震复发周期比较一致.在实际孕震能力判定的工作中,由于不同性质的断层段在同震过程中会表现更多的组合形式,断层发震能力判定结果存在更多的不确定性,但利用区域形变观测等资料给出震间断层运动特征的研究工作对于断层强震发震能力的判定具有非常重要的实际意义.     

重庆垫江MS4.4地震:一个可能发生在浅地表的地震

为更好地认识2016年8月11日重庆垫江M_S4.4地震的发震成因,本文首先分析了近场地震台站记录到的该地震产生的短周期面波（Rg）强频散现象,该现象表明本次地震具有明显的浅震特征。然后采用CAP方法反演了该地震的震源机制解和震源深度,结果显示本次地震是走向为近EW向的逆冲型地震,矩震级为M_W4.1,最佳震源深度为1—2 km。在此基础上,结合现场震害、烈度等震线长轴方位和地质剖面图等资料,分析认为在区域应力集中加载作用下,震源区沉积盖层内部的软弱地层产生挤压变形,导致地层内斜交的裂隙发生切层失稳滑动,从而诱发了垫江M_S4.4地震,同时排除了震源区NE走向的黄钦垭口地表断层为发震构造的可能。      

COSEISMIC DISPLACEMENT FIELD OF THE WENCHUAN EARTHQUAKE DERIVED FROM STRONG MOTION RECORDS AND APPLICATION IN SLIP INVERSION

The development of high-rate GNSS seismology and seismic observation methods has provided technical support for acquiring the near-field real-time displacement time series during earthquake. But in practice, the limited number of GNSS continuous stations hardly meets the requirement of near-field quasi-real-time coseismic displacement observation, while the macroseismographs could be an important complement. Compared with high-rate GNSS, macroseismograph has better sensitivity, higher resolution(100~200Hz)and larger dynamic range, and the most importantly, lower cost. However, baseline drift exists in strong-motion data, which limits its widespread use. This paper aims to prove the feasibility and reliability of strong motion data in acquiring seismic displacement sequences, as a supplement to high-rate GNSS. In this study, we have analyzed the strong-motion data of Wenchuan M_S8.0 earthquake in Longmenshan fault zone, based on the automatic scheme for empirical baseline correction proposed by Wang et al., which fits the uncorrected displacement by polynomial to obtain the fitting parameters, and then the baseline correction is completed in the velocity sequence. Through correction processing and quadratic integration, the static coseismic displacement field and displacement time series are obtained. Comparison of the displacement time series from the strong motions with the result of high-rate GPS shows a good coincidence. We have worked out the coseismic displacement field in the large area of Wenchuan earthquake using GPS data and strong motion data. The coseismic displacement fields calculated from GPS and strong motions are consistent with each other in terms of magnitude, direction and distribution patterns. High-precision coseismic deformation can provide better data constraint for fault slip inversion. To verify the influence of strong-motion data on slip distribution in Wenchuan earthquake, we used strong motion, GPS and InSAR data to estimate the stress drop, moment magnitude and coseismic slip model, and our results agreed with those of the previous studies. In addition, the inversion results of different data are different and complementary to some extent. The use of strong-motion data supplements the slip of the fault in the 180km segment and the 270~300km segment, thus making the inversion results of fault slip more comprehensive. From this result, we can draw the following conclusions:1)Based on the robust baseline correction method, the use of strong motion data, as an important complement to high-rate GNSS, can obtain reliable surface displacement after the earthquake. 2)The strong motion data provide an effective method to study the coseismic displacement sequence, the surface rupture process and quick seismogenic parameters acquisition. 3)The combination of multiple data can significantly improve the data coverage and give play to the advantages of different data. Therefore, it is suggested to combine multiple data(GPS, strong motion, InSAR, etc.)for joint inversion to improve the stability of fault slip model.     

用多种数据构建2008年汶川特大地震同震位移场

本文主要以GPS、精密水准观测和卫星SAR遥感图像分析2008年汶川特大地震同震位移特征.GPS数据包括:(1)四川盆地和川西高原地区各类国家等级GPS网点复测;(2)沿破裂带国家天文大地网GPS复测.前者推算的同震位移测定精度优于2 cm,后者6~8 cm.SAR遥感资料包括:(1)ALOS 卫星升轨相位干涉图像,精度优于8 cm;(2)ALOS和ENVISAT卫星影像合成的三维位移图,精度优于0.5 m.同震位移场显示,断层下盘(四川盆地)变形总体呈扇形集中指向震中,断层上盘(龙门山)变形总体上呈逆时针旋转态势,最大的实测水平位移5.5 m.汶川、理县、茂县等地测站位移指向破裂带方向,而平武、青川等地测站逐渐转变为平行,乃至远离破裂带方向,与汶川地震逆冲兼走滑的破裂特征一致.断层上盘大幅隆升,下盘靠近断层的区域以下沉为主,远场表现为幅度很小的隆升,垂直升降区域间,有一条与龙泉山断裂带平行的升降过渡带,调节龙泉断层的应力状态.用实测变形场检验多个地震波破裂模型表明,近场(距离断层50 km) 模型形变准确度可达40~50 cm, 远场精度优于5 cm.     

基于InSAR和GPS观测数据的尼泊尔地震发震断层特征参数联合反演研究

利用日本ALOS-2和欧空局Sentinel-1A卫星获得的尼泊尔地震同震形变场,结合GPS同震位移数据,联合反演了断层滑动分布特征和空间展布.结果表明:尼泊尔地震的同震形变场主要集中在150km×100km的范围内,且分为南北两个相邻的形变中心,南形变中心的视线向抬升量约为1.2m,北形变中心的视线向沉降量约为0.8m,均位于发震断层上盘.位于形变抬升区的KKN4和NAST两个GPS站,抬升量和南向运动量均达到了m级,而远离震区的其他GPS台水平和垂直观测量均在1cm以内.联合反演得到的断层位错分布主要集中在沿走向150km,沿倾向70km的范围内,最大滑动量为5.59m,平均滑动量为0.94m.断层面倾角在浅部约为7°,随着深度增加,倾角逐渐变大,到垂直深度20km时倾角接近12°;5月12日MW7.2级余震位于主震破裂区的"凹"型滑动缺损区域;主震破裂区的上边界与MBT空间位置十分吻合,主震破裂区主要集中的MBT以北50~60km处,垂直深度为8~9km,倾角为9°,继续向北时主震破裂面以10°~12°的倾角向深延伸,在18~20km可能与MHT交汇.因此,初步判定MBT为此次地震的发震断层.     

GNSS观测的九寨沟7.0级地震同震位移初步结果

2017年8月8日四川省九寨沟县发生了7.0级地震,中国大陆构造环境网络与北斗地基增强系统的GNSS连续观测共同监测到了此次地震的同震位移（坐标：东向为正,北向为正）,结果显示：3个站点记录到了明显的同震位移,距离震中43 km的九寨沟台站（SCJZ）在东西向的位移为-9.8±1.5 mm,在南北向的位移为3.3±0.7 mm;距离震中65 km的松潘站（SCSP）在东西向的位移为-1.8±0.7 mm,在南北向的位移为-7.7±0.6 mm;距离震中77 km的舟曲站（GSZQ）在东西向的位移为0.4±1.2 mm,在南北向的位移为3.6±0.8 mm.通过同震位移分布特征,可以推测此次地震为一次左旋走滑型事件,引起水平向同震位移大致不超过150 km范围,地震对东南侧的龙门山断裂带影响非常小,对北侧的塔藏断裂和西侧的岷江断裂处引起的同震位移为厘米级.同震位移的反演结果显示：断层面上滑动量主要集中在7 km深度,最大量值约为0.4 m,平均滑动角为-15°,利用滑动分布计算的相应矩震级为M_W6.4,与地震波反演结果相当.结合同震滑动分布、同震主应变分布、余震分布和震源机制解等特征,推测此次地震破裂极值区累积的能量得到较充分释放,进一步分析得出此次地震在塔藏断裂、岷江断裂和虎牙断裂处产生了一定的应力变化,值得持续关注.     

利用GPS观测资料分析2008年于田Ms7.3地震的同震位移及震后形变

2008年3月21日新疆于田发生Ms7.3级地震.本文通过处理、分析GPS数据,得到破裂断层北侧100 km附近的同震位移及震后形变信息.在观测区域GPS点监测到10 mm左右的同震位移,其中最大为南向14 mm,东向5 mm.同震位移呈现一致性的东南向运动特征,证实于田地震存在显著的左旋走滑分量.震后台站向西南方向运...