利用 Yabuki＆Matsu'ura 反演方法计算2011年日本东北地区太平洋海域Mw9.0级地震同震滑动分布

Yabuki & Matsu'ura反演方法是利用ABIC最佳模型参数选取方法和平滑的滑动分布作为约束条件,由形变观测数据计算发震断层滑动分布.本文基于日本列岛同震GPS观测数据和发震断层曲面构造模型,利用Yabuki&Matsu'ura反演方法计算2011年日本东北地区太平洋海域Mw9.0级地震的发震断层同震滑动分布.反演结果表明,断层面上的最大滑动量为35 m,较大滑动分布在浅于30 km的震源中心上部,最大破裂集中在20 km深度的地方,其地震矩约为3.63×1022N·m,对应的矩震级为Mw9.0.模拟结果显示Yabuki&Matsu'ura反演方法更适用于倾角低于40°的断层模型反演.最后,本文基于上述方法获得的发震断层滑动模型,利用地球体位错理论正演计算该地震在中国及其邻区产生的远场形变,正演计算结果基本可以解释由中国GPS陆态网络观测到的同震形变.     

利用Yabuki & Matsu'ura反演方法计算 2011年日本东北地区太平洋海域Mw9.0级地震同震滑动分布

Yabuki & Matsu'ura反演方法是利用ABIC最佳模型参数选取方法和平滑的滑动分布作为约束条件,由形变观测数据计算发震断层滑动分布.本文基于日本列岛同震GPS观测数据和发震断层曲面构造模型,利用Yabuki & Matsu'ura反演方法计算2011年日本东北地区太平洋海域Mw9.0级地震的发震断层同震滑动分布.反演结果表明,断层面上的最大滑动量为35 m,较大滑动分布在浅于30 km的震源中心上部,最大破裂集中在20 km深度的地方.其地震矩约为3.63×1022N·m,对应的矩震级为Mw9.0.模拟结果显示Yabuki & Matsu'ura反演方法更适用于倾角低于40°的断层模型反演.最后,本文基于上述方法获得的发震断层滑动模型,利用地球体位错理论正演计算该地震在中国及其邻区产生的远场形变,正演计算结果基本可以解释由中国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, 与地震学反演的结果比较吻合.      

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

基于黏弹性球体位错理论,联合陆地和海底同震GPS数据以及日本本岛330个陆地GPS站点5～10年的震后数据,反演了日本Mw9.0地震的断层滑动模型,提升了断层滑动分布在细节上的合理性.首先,基于日本本岛330个陆地GPS站点震前2年和震后10年的连续观测数据,获取了日本Mw9.0地震震后5～10年的年平均位移,该时段的...     

利用InSAR获取2020年新疆于田Mw6.3地震同震形变场与断层滑动分布反演

本文基于InSAR技术,利用欧空局Sentinel-1A/B升降轨SAR数据,提取了2020年6月26日新疆于田Mw6.3地震的同震形变场.利用升、降轨同震形变场约束,分别采用MPSO算法和Bayesian方法反演此次地震发震断层均匀滑动的几何参数,并进行对比.然后采用SDM方法获得发震断层非均匀滑动分布,并分析了同震...     

利用InSAR资料反演缅甸Mw6.8地震断层滑动分布

2011年3月24日缅甸东北部发生M_w6.8级地震.本文利用覆盖该地区的升降轨ALOS PALSAR数据,获取了该次地震的同震形变场,并采用灰度配准技术获取了其地表破裂位移.针对影像中因轨道不精确造成的非线性长波长误差,本文采用二次多项式曲面法予以去除,获取了更为精确的同震形变场.最后,基于弹性半空间形变模型反演了该地震的断层滑动分布.结果表明,该地震断层滑动以左旋走滑为主,兼具少许的倾滑运动分量,断层滑动主要集中分布在断层面0~10 km深度范围,最大滑动量达5 m,位于地表以下5 km深处.反演获得的地震标量矩为2.49×10¹⁹N·m,震级约为M_w6.8级.     

超导重力仪检测2011年日本东北Mw9.0级地震前的重力扰动信号

对大地震前的扰动现象的研究有助于认识地震孕育的动力学过程,震前重力扰动已成为关注的热点之一. 对大地震前的扰动现象的研究有助于认识地震孕育的动力学过程,震前重力扰动已成为关注的热点之一. 为检验日本M_w9.0级地震是否存在震前扰动现象,本研究利用全球超导重力仪记录到的地震前后7天内20组秒采样数据进行分析. 经潮汐、大气改正等处理去除仪器的漂移及残余潮汐效应,得到非潮汐重力变化曲线.结果表明大部分振幅大于30×10^-8 m·s^-2的曲线反映了全球M_w≥6级地震引起的高频波动信号,其中11组数据在3月9日M_w7.3级前震之前出现了扰动现象.震前扰动可分解为三个频段,其中,低于0.1 Hz和高于0.18 Hz的分量分别反映了地震波动信号及非构造信息,中间频段(0.118~0.18 Hz)信号能够较大程度地压制地震波动信号、并同时保留异常扰动信息.它的振幅在3月7日10时之前基本保持约1×10^-8 m·s^-2,之后开始逐渐增大,到3月9日7.3级前震前后达到最大,此后振荡衰减,振幅保持约(5~10)×10^-8 m·s^-2,直至主震发生.中间频段信号的变化特征与主震前的应力迁移过程以及实验记录到的地震成核过程有许多相似之处;不过,震前重力异常是否与主震前的应力加速积累有关,仍待进一步研究.     

利用1-HzGPS波形数据反演2011年日本北大地震震源破裂过程

利用1-Hz GPS数据反演2011年3月日本东北大地震的时空破裂过程,通过理论波形图与观测波形图的拟合,以及ABIC准则的判定,反演得到滑移分布的稳定估计.结果表明:破裂集中在震源附近,并在震中以下30 km处达到最大位错量,约72 m,在岩手和福岛靠近海岸线处有两个较小的破裂区,断层的南半段破裂很小,同时也是余震最为密集的地方.整个破裂过程由3个破裂阶段组成,发震后断层从震源向周围缓慢破裂,在50～60 s破裂延伸至地表(海沟);60～90 s在震中下方迅速发生强烈破裂,形成了最大位错区,这一阶段主要是倾向上的双边破裂;90～120 s在断层最北西处和断层南半段发生较小破裂,整个破裂过程持续大约120 s.反演得到的地震矩为3.8×1022 Nm,相应的矩震级为MW 9.0.     

2011年日本Mw9.0地震引起的环渤海地区井水温度变化分析

2011年3月11日的日本M_w9.0地震在中国引起了大范围的井水温度同震变化,而水温震后持续变化的井孔多数分布在井网密度较大、距震中较近的环渤海地区。本文分析了环渤海地区水温震后变化的特征和机理。结果表明,水温震后变化形态为上升,变幅0. 005~0.976°C,空间分布在同震位移较大、张性应变较明显和地震能量密度较大的区域;依据同井水位资料,一些井孔的水温和水位震后变化特征较一致,其水温震后升高可能是地震波增大含水层渗透性的结果;另一些井孔的水位震后变化不显著,其水温震后升高可能与地震波增大井区大地热流有关。     

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

利用超导重力数据检测日本Mw9.0地震的低频自由振荡及谱线分裂

准确估计低频自由振荡及谱线分裂是约束地球内部结构和改进地球模型的重要手段.本文利用四个不同台站的超导重力观测数据系统研究了日本M_w9.0大地震激发的低于1.5 mHz自由振荡及谱线分裂.研究结果表明:(1)选取适当的数据长度,超导重力观测数据可以检测出低于1.5 mHz除₁S₁以外的所有自由振荡;(2)除₀S₂、₀S₃、₀S₀、₂S₁、₃S₁、₁S₂和₀T₂外,重点探测出₃S₂、₀S₄和₁S₄谱线分裂的所有谱峰;(3)与PREM模型理论频率相比,₀S₀观测频率平均向右偏移0.354×10^-3mHz,说明PREM理论模型中地幔底部参数与真实地球可能存在微小偏差;(4)₃S₂的谱线分裂率r为1.485267,比PREM理论谱线分裂宽度约宽50%,表明PREM中地球内核中部介质参数可能存在一定误差,需要进一步改善.另外,quasi-₀T₂的r为1.254206,比PREM理论谱线分裂宽度约宽25%.     

2016年新疆阿克陶MS6.7地震同震形变特征和断层滑动分布反演

本文针对2016年11月25日在新疆阿克陶发生的地震, 用差分干涉测量技术(D-InSAR)对3种不同观测模式的升、 降轨数据进行处理, 提取了多视线向的同震形变场; 根据不同模式的LOS向形变量, 构建形变分解模型, 将其分解为垂直向形变量和沿断层走向形变量; 结合同震形变场特征与震源机制解, 采用单断层模型, 利用梯度下降法(SDM), 以Multi-LOS向形变进行约束, 反演了阿克陶地震的同震滑动分布特征。 研究结果表明, 升、 降轨LOS向同震形变场在发震断层两侧具有不同的形变特征, 发震断层走向近EW向; LOS向形变量分解表明, 此次地震破裂以右旋走滑为主; 滑动分布反演的形变残差介于0~5 cm之间, 发震断层的滑动量主要位于2~16 km深部, 最大滑动量可达1.02 m, 位于断层面深部5.83 km处, 最大滑动量处的滑动角为185.24°; 平均滑动角为181.32°, 平均滑动量为0.12 m; 滑动分布反演也证明该地震为右旋走滑破裂事件, 与LOS向形变分解结果一致; 当剪切模量μ=3.2 GPa时, 反演得到的地震矩震级约M_W6.6。     

2010年1月12日海地 Mw7.0地震InSAR同震形变观测及同震滑动分布反演

2010年1月12日GMT时间21时53分,在海地境内(72.57°W,18.44°N)发生了Mw7.0地震.文中利用干涉合成孔径雷达(InSAR)方法获得了覆盖整个震区的高精度形变观测资料,用以研究该地震的发震机理.采用ALOs PALSAR数据,分析了轨道、大气等误差源对干涉信号的影响,最终获得了雷达视线向(LOS...     

2022年青海门源MW6.9地震同震形变及断层滑动分布反演

2022年1月8日,青海海北州门源县发生M_W 6.9地震,震中位于青藏高原东北缘祁连地震带上冷龙岭断裂和托莱山断裂的交汇处.门源地震活动强烈,造成地表破裂明显,因此研究门源地震的发震机理,对评估周边区域及主要断裂的地震危险性具有重要意义.本文基于D-InSAR技术,利用升降轨Sentinel-1A SAR数据,获取门源地震的同震形变场.结果显示,同震形变主要集中在冷龙岭断裂和托莱山断裂的交汇处,形变长轴整体呈NWW-SEE,同震引发的升轨隆升形变量0.40 m,沉降量0.65 m,降轨隆升形变量0.80 m,沉降量0.70 m,升降轨视线向(LOS)形变呈现符号相反大小相近的特征,断裂运动以左旋走滑为主.以升降轨同震形变信息为约束条件,基于弹性半空间位错理论,采用两步法进行反演,获取门源地震断层的几何参数.结合该区域的地质构造得到断层面上的精细滑动分布,结果表明,断层破裂延伸至地表,破裂迹线长度达22 km.断层滑动主要集中在地下2～12 km,最大滑动量为4.2 m,位于地下7 km处.断层走向约为109°,倾角约为82°,释放的地震矩为2.67×10^...     

基于方差分量估计的2015年尼泊尔MW7.8地震同震滑动分布

基于三角位错的断层自动剖分方法,联合ALOS-2 InSAR和GPS同震位移数据反演了2015年尼泊尔M_W7.8地震的发震断层参数及滑动分布.反演中基于赫尔默特方差分量估计方法确定InSAR、GPS水平向与垂直向的3种观测位移数据的相对权比.反演结果显示,方差分量估计方法定权得到的最大滑动量5.5 m,大于验前方差定权给出的5.1 m的结果,同时前者显示同震破裂沿东南方向传播至接近M_W7.3最大余震区域,然后分别继续向北和东南传播约30 km,并形成一个滑动空区（slip gap）,而后者的结果没有表现出这种显著的滑动空区特征.方差分量估计定权得到的滑动模型能更好地解释观测到的InSAR同震形变场,其拟合残差标准差较验前方差定权的结果减小了2.5 cm,GPS的东和北分量的残差标准差分别减小了0.4 cm和0.9 cm.与两种不同的基于矩形位错的方法相比,利用方差分量估计的三角位错断层自动剖分方法的滑动误差均值和标准差相对最小,其中均值为0.037 m,标准差为0.028 m,表明该方法得到的滑动模型具有相对更好的无偏性和有效性.

     

2011年日本9.0级地震前地震空间相关长度增长现象

使用单键群方法研究了2011年3月11日日本本州东北部海域9.0级地震前地震空间相关长度变化。结果表明,在本次地震发生前约两年,5级以上地震的空间相关长度开始出现明显增长,变化趋势可以用幂次率进行很好的拟合。用没有7级以上强震发生的时段内的地震目录,采用相同方法进行了计算,没有发现明显的幂次率增长现象,说明强震前5级以上地震的空间相关长度的幂次率增长是地震孕育过程中震源区的物理性质特征的表现之一。     

利用GPS和水准测量资料反演2001年昆仑山口西8.1级地震的同震滑动分布

通过分析 2 0 0 1年昆仑山口西MS8.1地震前后的GPS资料和 1 979— 2 0 0 2年的水准测量资料 ,获得了地表同震位移场。利用这些同震位移数据 ,以震后详细野外地质调查破裂数据为约束 ,反演了破裂断层面上的同震滑动分布。结果表明 ,破裂下界深度为 1 4 .2～ 2 1km(70 %置信区间 ) ,最佳破裂深度 1 7km。虽然在太阳湖段和主破裂带西端的中间区域未观测到地表破裂 ,但反演的结果表明此区域存在 2～ 3m左右的左旋水平走滑 ,与InSAR资料分析得到的结果一致。地表以下的破裂西端止于太阳湖段 ,而东端似乎在地表破裂迹线以东 30km范围内仍存在 1 .5～ 2 .0m的左旋滑动。反演的垂直位移表明断层在东经 93°以西部分大体表现为南盘上升 ,而东经 93°以东部分表现为北盘上升。由大地测量和地表破裂调查得到的地震矩释放为 6 .1× 1 0 2 0 N·m ,与地震波资料的反演结果基本一致     

2022年泸定6.8级地震GNSS同震形变场及其约束反演的破裂滑动分布

2022年9月5日四川省甘孜州泸定县发生6.8级地震,地震发生在川滇菱形块体东边界鲜水河断裂磨西段西侧附近.本文基于震中350 km范围内中国大陆构造环境网络和中国地震科学实验场118个GNSS连续站数据,观测获得了精细的同震形变场,结果显示：同震形变跨发震断层呈空间对称分布,表明本次地震具有显著的左旋走滑特征;记录到的最大同震形变发生在震中距40 km的SYD5(石棉安顺场)站,东西向和南北向形变量分别达到-22.0±1.2 mm和11.6±0.9 mm,震中距100 km以外测站的水平同震形变均小于5 mm;垂向同震形变不显著.结合走向163°和倾角77°的发震断层模型,本文对发震断层面的同震破裂滑动分布进行了反演,结果显示：同震滑动主要集中在主震东南侧余震空区0~15 km深度范围内,且破裂达到了地表,同震释放的地震矩与矩震级为M_W6.57地震相当.结合理论同震形变场、主应变场和邻近区域主要活动断裂库仑应力变化的空间分布特征,本次可能导致处于强闭锁和地震空区的安宁河断裂石棉—冕宁段未来发震风险性增强.

     

基于GPS观测分析日本9.0级地震同震位错与近场形变特征

2011年3月11日日本本州宫城县东海岸近海发生MW9.0级地震,本文在对GPS同震位移场分布及误差特征分析的基础上,反演了同震位错分布.误差分析结果表明震源北西向300 km、北北西向550 km、南西向700 km范围内的同震位移量值明显大于误差,可以为位错反演提供有效的地表位移约束.沿震源北西向GPS剖面结果和位错反演位移剖面结果均表明同震近场位移符合指数衰减特征.位错反演结果表明,日本9.0级地震最大同震位错为25.8 m,位于震中附近;位错量大于10 m的同震破裂集中在震中附近400 km范围内;日本海沟南段同震位错量相对较小,此次地震为日本海沟地区典型逆冲型地震.根据此次9.0级地震和该地区以往强震破裂空间分布特征,此次9.0级地震破裂既体现了强震原地复发的特点,又体现了强震破裂的填空性.     

2008年于田MW7.2地震滑动分布反演及三维同震形变场解算

基于ENVISAT ASAR升降轨数据, 利用InSAR获取2008年于田M_W7.2地震同震形变场; 采用SDM反演本次地震断层滑动分布; 使用PSCMP正演获取于田M_W7.2地震南北向模拟形变量, 并结合升降轨同震形变场, 解算三维同震形变场。 同震形变场分析表明, 2008年于田M_W7.2地震以正断为主, 且带有走滑运动特征, 破裂带走向为NNE向。 同震滑动分布反演结果显示, 断层沿走向被分为4段F1、 F2、 F3、 F4, 其滑动分布集中在0~14 km区间, 以F2、 F3段为主, 最大滑动量约5.31 m, 位于F2段深部2.76 km处; 沿破裂带走向, 左旋走滑位移与垂直位移比值有增大的趋势; 反演获得的地震矩M₀=5.58×10¹⁹N·m, 相当于矩震级M_W7.1。 三维同震形变场解算结果显示, 断层上盘整体表现为沉降, 断层下盘整体表现为隆升, 且沉降量明显大于隆升量, 表明地震以正断破裂为主; 除靠近断裂带中上部表现为向东南运动外, 上盘整体上表现为向西南运动; 断层下盘则整体表现为向东北运动, 证明破裂兼有左旋走滑运动。 滑动分布反演、 正演与三维同震形变场解算结果皆表明, 于田M_W7.2地震破裂以正断为主, 且带有一定的左旋走滑。