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;同震水平位移场模拟结果显示贡嘎错断裂、康西瓦断裂和普鲁断裂等不同位置主应变特征具有差异性,这种差异特征是否影响断裂带以及周围区域的应力构造特征,值得关注.     

1999年台湾集集地震震源破裂过程

使用GPS同震位移资料和远场P波记录,研究了1999年台湾集集地震震源破裂过程．根据地质构造和余震分布引入了一个由弯曲断层面构造的三段“铲状”断层模型．在使用静态GPS位移资料反演集集地震的断层破裂滑动分布时,由于集集地震断层北部近地表破裂的复杂性,在位错模型中考虑拉张分量对地表同震位移的贡献,可更好地同时拟合GPS观测资料的水平和垂向分量．而纯剪切位错弹性半空间模型和分层地壳模型都无法同时拟合水平和垂向GPS观测资料．在此基础上,同时使用静力学同震位移资料和远场地震波形记录,反演集集地震的震源破裂过程．结果表明,一种垂直于断层面的“挤压性”(负)拉张分量几乎集中分布于地震断层的浅部和北部转折处,而这一带地表破裂远较没有(负)拉张分量出现的南部断层复杂．“冒起构造”的数字模拟表明,这种在集集地震破裂转折处及北部断裂带广为出现的典型破裂造成的地表位移可以用具有负拉张分量(挤压)的逆冲断层更好地模拟．而这种负拉张分量(挤压)的分布正是地震破裂性质和几何复杂性的综合反映,震源破裂过程也显示北部转折处破裂在空间和时间上的复杂性.高滑区域与余震分布表现为负相关.     

基于GPS和强震资料反演汶川8.0级地震的同震滑动模型

通过震中附近GPS同震位移资料,采用SDM反演法,应用均匀介质模型和分层地壳结构模型分别反演汶川8.0级地震的同震滑动,并加入强震资料进行反演对比分析,结果表明:两种模型反演的同震滑动分布与发震断层的科考结果吻合,分层地壳结构模型的反演结果整体上要优于均匀地壳结构模型的反演结果;GPS与强震数据分别反演得到的同震位移方...     

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,与地震波反演结果相当.结合同震滑动分布、同震主应变分布、余震分布和震源机制解等特征,推测此次地震破裂极值区累积的能量得到较充分释放,进一步分析得出此次地震在塔藏断裂、岷江断裂和虎牙断裂处产生了一定的应力变化,值得持续关注.     

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左右地震强震记录的同震位移场,并反演震源滑动模型具有一定的可行性,且其同震位移场和滑动模型可为震后灾害快速评估、救援力量分配、余震趋势判定等快速应急响应工作提供参考依据.     

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

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

均匀弹性分层介质模型中的同震位移、应变以及应力

地震发生时产生的同震位移、应变以及应力变化,特别是同震应力的变化,在地震触发等问题的研究中有着重要的意义.本文进一步发展了基于均匀弹性水平层状介质,利用广义反射透射系数矩阵和离散波数计算同震位移的方法,使之可以计算相应的应变、应力以及同震库仑应力变化.可适用于多种情况,接收点可以位于地表以及地表以下,震源类型可以是剪切位错源以及拉张位错源.通过与半无限介质的解析解相比较,结果一致,验证了方法的可靠性以及计算精度,可以用于计算地震之后库仑应力变化,为判断余震分布提供参考.在计算同震位移时,使用了梯形积分与Filon积分相结合的积分方式,即提高了同震位移计算的速度,又保证了计算精度,有利于反演问题研究.     

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

芦山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级地震同震位移场存在的困难和问题,为今后类似工作提供参考.     

基于形变观测分析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级左右地震复发周期比较一致.在实际孕震能力判定的工作中,由于不同性质的断层段在同震过程中会表现更多的组合形式,断层发震能力判定结果存在更多的不确定性,但利用区域形变观测等资料给出震间断层运动特征的研究工作对于断层强震发震能力的判定具有非常重要的实际意义.     

GPS测定的2020年阿拉斯加MW7.8地震震时地表动态变形及永久同震形变场

以UNAVCO公布的阿拉斯加地区GPS 1 Hz和30 s采样观测数据为基础,采用双差定位方法分别对两种采样率的数据进行处理分析,获得2020年7月22日阿拉斯加M_W7.8地震震时地表动态变形及同震三维形变场。结果显示,震中270 km范围内高频GPS震时波形明显,最大振幅达600 mm。根据各个GPS站的动态形变波形振幅及响应时间认为,其振幅和响应时间受地震的破裂传播方向和场地效应影响较大。静态同震位移矢量指向震中,同震位移大小基本符合随震中距离增大而减小的特征,除站点AC13外,其中距离震中最近站点的最大水平位移达26.7 cm。GPS测定的同震形变表明,2020年阿拉斯加地震是发生在阿拉斯加阿留申俯冲大断裂的一次逆冲型地震。     

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连续观测站的东西分量明显减速,有的甚至减至零,是明显的短期前兆异常。此外,本文还进一步讨论了两种同震位移及其成因,并推断地壳水平运动挤压是此次地震的成因。      

利用GPS和InSAR数据反演2011年日本东北MW9.0地震断层的同震滑动分布

本文首先对Envisat/ASAR数据进行干涉处理, 获取2011年日本东北M_W9.0地震的地表InSAR同震形变场; 然后通过对InSAR同震形变数据重采样方法的深入分析, 选择条纹率法结合干涉图的空间相干性对InSAR同震形变数据进行重采样; 最后基于弹性半空间位错模型, 联合InSAR与GPS形变数据, 采用最小二乘法反演发震断层的滑动分布. 研究结果表明: 考虑相干性的条纹率重采样方法, 更适用于形变场中存在除断层外的有限边界、 且形变场范围较大的InSAR数据重采样处理; 断层滑动主要发生在地表以下50 km范围内, 最大滑动量为49.9 m, 矩张量为4.89×1022 N·m, 所对应的矩震级为M_W9.1, 与地震学反演的结果比较吻合.      

日本MW9.0地震震区及其周缘2002—2015年卫星重力变化时空特征

本文在考虑全球陆面数据同化系统陆地水储量变化影响后,利用2002年4月—2015年4月GRACE卫星RL05月重力场模型获取了2011年日本M_W9.0地震震中及其周边区域的重力场信息;然后给出了日本及其周边2003—2015年的年度累积重力变化和差分重力变化,并且利用经验正态函数方法深入分析了该地震过程中同震重力变化对区域重力场的贡献. 结果显示:日本M_W9.0地震前研究区域出现了幅值高达2×10-8 m/s2的异常重力变化,其同震效应的时间和空间特征均出现在第一模态,且同震重力变化和震后区域重力场变化特征显著,这充分表明该地震对区域重力场的影响显著.      

GPS资料反映的日本东北MW9.0地震的孕震特征

在对海洋板块俯冲型地震孕震模式的分阶段变形特征进行分析的基础上,结合同震位错反演结果分析了M_W9.0级地震的同震变形特征,分析过程中通过精度检验讨论了该反演结果的可靠性.通过对震前和同震GPS结果的对比分析,讨论了二者的差异性.GPS应变(率)结果表明,震前日本岛应变积累主要反映了太平洋板块和菲律宾板块的俯冲作用,同震结果表现为指向震源方向拉张应变的释放,对日本岛的主要影响区域介于35°N和43°N之间.震前和同震GPS剖面结果反映的变形特征具有互补性,但量值相差上百倍.震前的GPS速度、应变率剖面和时间序列结果表明靠近日本岛东海岸一侧可能存在变形趋于极限现象.通过对此次地震可能前兆的分析表明,中长期预测方面、震前GPS时间序列的趋势性偏离、前震活动、震源区b值降低等现象在一定程度上反映了此次地震的孕震特征.     

日本熊本MW7.0地震同震位移场和震源滑动模型反演

利用2016年4月16日日本熊本M_W7.0地震震中周围94个近场强震动台的观测资料和新近改进的强震经验基线校正方法SMBLOC,尝试解算并绘制了一个内陆M7左右走滑型地震的同震位移场全貌,并反演了其震源滑动模型.与日本国土地理院(GSI)公布的该地震57个GPS同震位移结果的比较显示,两种完全不同资料、不同解算方法给出的水平同震位移场的最大幅值均为100 cm左右,均呈右旋走滑为主兼具部分正断分量的震源机制.强震最大水平和垂直永久位移分别为104.5 cm和58.0 cm, 分别出现在震中东北侧的KMMH162台和KMM005台.两种资料单独以及联合反演的震源滑动模型均表明,此次地震为北东侧破裂为主并呈双事件特征,且主要滑动均不在初始破裂点附近, 而是集中于第二次事件周围,即距离初始破裂点东北侧约20 km处的走向长约40 km、倾向宽约20 km的范围内.基于强震和GPS模型所得的最大滑动量分别为5.10 m和5.87 m,量级一致,反演矩震级均为M_W7.1左右;主破裂区近地表滑动量比野外调查结果略微偏大,可能与数值效应有关.此外,还利用不同方法得到的解算结果比较了熊本地震特有的12组台间距在3 km以内的GPS-强震台站对各自的三分量同震位移,其结果表明对于M7左右的地震而言,SMBLOC方法解算同震位移时方向和幅值的可靠性下限约为2 cm.