2011年日本M_W9.0地震震后形变机制与震源区总体构造特征

本文利用较为完备的球体位错理论,结合4.5年的震后位移数据,优化了2011年日本M_w9.0地震震源区岩石圈弹性层厚度与地幔黏滞性因子,更新了该强震断层余滑时空演化过程.首先,基于日本列岛215个均匀分布的GPS连续观测站震前2年与震后4.5年的观测数据,提取了2011年日本M_w9.0地震引起的震后位移时空变化;接着,依据断层余滑衰减相对较快的特点,利用黏弹性球体位错理论对震后3~4.5年的GPS观测数据进行反复拟合,确定2011年日本M_w9.0地震震源区地幔黏滞性系数和岩石圈弹性层厚度的最优解分别为6×10^(18)Pa·s和30 km;然后,从震后3年内GPS观测数据中剔除地幔黏滞性松弛效应,获取了断层余滑对应的震后位移场;最后,利用基于球体位错理论的反演算法,反演了2011年日本M_w9.0地震断层余滑的时空演化过程.结果表明,2011年日本M_w9.0地震引起的断层余滑在震后半年内变化显著,震后2年主震区域余滑基本停止,断层的两端存在一定的余滑效应,断层余滑的累计矩震级达到8.59;地震后4年,地幔黏滞性松弛效应对震后位移场的贡献在总体上超过断层余滑的贡献.

利用粘弹性球体位错理论研究2011年日本M_W9.0地震引起的震后位移时空变化

基于2011年3月~2015年9月的GPS观测数据,研究了2011年日本M_W9.0地震引起的震后位移时空变化,分析震后4.5a内震后断层余滑和粘滞性松弛对震后位移场的贡献。基于子断层叠加的思想,对Tanaka的粘弹性球体位错理论配套计算程序(环型解部分)加以改进,克服其近场计算精度不足(甚至错误)的缺陷,提高了近场震后位移的计算精度。利用弹性球体位错理论配套计算程序和改进后的粘弹性球体位错理论配套计算程序,分别计算2011年日本M_W9.0地震产生的同震位移,两组结果的差异仅占信号的1%左右,验证了对粘弹性球体位错理论配套计算程序改造的正确性。最后,结合震后3~4.5aGPS观测数据,利用Tanaka的粘弹性球体位错理论确定了2011年日本M_W9.0地震震中周边区域地幔粘滞性系数,其最佳拟合值为6×10^(18)Pa·s。

2011年日本东北大地震（M_W=9.0）震间与震前变形场特征及其对强震预测的启示

强震震前(preseismic)动力学过程的研究对于地震预测具有十分重要的意义,但由于观测资料的限制,目前对强震前孕震区力学状态及其演化过程的认识还非常有限.2011年日本东北9.0特大地震(Tohoku-Oki)发生在GPS观测台站最为密集的地区,为研究特大地震震间(interseismic)与震前的变形状态提供了难得的机会.文中将利用日本东北大地震之前连续的GPS观测资料,分别计算震间与震前的速度场与变形场.通过对比分析发现,日本东北地区(Tohoku)震前的应变状态与震间的有很大的不同,震间的变形主要受到太平洋板块向日本海沟北西西向的俯冲挤压作用所控制,其主压应变以近东西向压缩为主,日本东北地区的运动方向与太平洋板块的运动方向大体一致.但是,临近地震前(震前)日本东北地区的运动方向发生了很大变化,震前30天的连续GPS观测结果显示,速度场的优势方向经常变换,间歇性地出现与太平洋板块运动方向相反的情况.这意味着震前孕震区的力学状态发生了很大的改变.这种变化可能与震前破裂成核或慢滑移及慢地震等过程有关,这些过程将加速或促进大地震的发生,从而为大地震的发生准备了力学条件.值得特别强调的是,这些现象都是可以通过直接观测能够发现的大地震之前的异常现象.由此可见,加密GPS站点进行连续观测,寻找震前变形异常区以及探索异常的物理机制对于地震预测预报有重要的科学意义.

利用同震和震后位移数据联合反演2011年日本M_(W)9.0地震同震断层滑动

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

2011年日本MW9.0地震引发的海啸对地震背景噪声的影响

基于北美沿岸和内陆地震台站的连续地震波形记录,并结合沿岸台站附近布设的DART系统记录的海底压力数据以及预测潮汐数据,利用时频分析和极化分析方法对2011年3月11日日本东北部海域M_(W)9.0大地震所激发的海啸对地震背景噪声所产生的影响予以深入分析。结果显示:海啸对高频噪声(1.3—1.5 Hz)以及短周期双频微地动噪声(0.18—0.4 Hz)的影响较小,但海啸显著增强了长周期双频微地动(0.1—0.15 Hz)、单频微地动(0.05—0.08 Hz)以及地球背景自由振荡(0.004—0.007 Hz)的振幅,且随着噪声频率的降低,这种振幅增强的影响更明显,影响的持续时间也更长;海啸到达近岸时,对附近台站的各频段噪声均有影响,成为各频段噪声的主控能量来源,且其位置在后续过程中会随时间变化。这表明海啸对噪声特征的影响与海啸传播特性有关,即海啸在传播过程中因受水深、海底地形以及近岸地形的反射、衍射等的影响,能量聚集区域随时间而变化,并非均匀地传播到海岸,从而导致了不同频带噪声主极化方向随时间的变化。

日本M_(W)9.0地震前后郯庐断裂带中南段构造应力场动态变化特征

本文以郯庐断裂带中南段及周缘为研究区,计算和收集了1970—2016年1148次地震震源机制解,采用阻尼应力张量方法反演了构造应力场,通过分析应力场时空动态变化过程和地震活动性特征,探讨了2011年3月11日日本9.0级地震对郯庐断裂带中南段及周缘的影响.结果显示:9.0级地震对郯庐断裂带中南段及周缘造成了显著影响,表现为应力场出现扰动变化和地震活动显著增强.以挤压为主的背景应力场呈现出一定的拉张作用,应力结构及应力方向发生变化,最大主应力σ_(1)轴倾伏角增大,中等主应力σ_(2)轴倾伏角减小,最大主应力σ_(1)轴方向发生小幅逆时针偏转.包括研究区在内的中国陆地东部地区地震活动显著增强,中强地震丛集活动.应力场扰动变化和地震起伏活动存在较好的相关性,并先后呈现出3个不同阶段,扰动导致的构造应力作用方式的阶段性改变是地震活动增强的动力机制.结果还显示,应力场扰动变化和地震活动增强开始于9.0级地震前1~2年,揭示了9.0级地震的影响不仅限于同震和震后,在其孕震后期已经对远场的中国陆地东部产生了影响.本文结果可以为日本俯冲带大震对中国陆地构造活动影响的相关研究提供参考依据.

利用球形地球位错理论研究日本M_(W)9.0地震引起的稳态变形与库仑应力变化

稳态地震变形是指由地震引起的、超长时间尺度内震后黏弹性松弛效应的累加,通过对稳态地震变形的研究,可以研究大地震导致的地球总体震后形变场,有利于分析周边地区对该地震的极限响应.本文首先利用球形地球地震位错理论,计算了日本M_(W)9.0地震在日本列岛和中国东部地区产生的稳态水平位移.结果表明,日本M_(W)9.0地震在日本本岛引起的稳态水平位移最大可达到8 m,大体上是同震水平位移的两倍;在中国东北地区产生的稳态水平位移最大值超过0.8 m,在华北地区产生的稳态水平位移最大值为0.5 m左右,是同震水平位移的30~40倍.接着,本文根据前人研究成果,结合地质资料,建立了中国东北和华北地区15条主要断裂带的几何与滑动构造模型,然后利用球形地球位错理论计算了日本M_(W)9.0地震在上述断裂带上引起的同震与稳态库仑应力变化.结果表明,该强震在中国东北和华北地区主要断裂带上引起的同震库仑应力变化均不超过1 kPa;稳态库仑应力变化显著超过同震结果,但除了依兰—伊通断裂带北段和嫩江断裂带南段的一部分以外,其他断裂带上的稳态库仑应力变化均未超过10 kPa这个地震触发阈值.因此日本M_(W)9.0地震对中国东北和华北地区地震危险性影响有限,不会显著改变区内主要断裂带的应力状态和地震活动性.

基于粘弹性球体地球模型的震后位移与重力变化计算软件

在考虑地球曲率、成层结构、可压缩性与自重的前提下,Tanaka等[1-2]提出一套较为完备的粘弹性球体位错理论,可计算全球任意位置由地震产生的同震与震后形变(含位移、重力变化、大地水准面变化)。Gao等[3]给出了与上述理论相匹配、界面友好的计算软件,能计算30个震后时间点对应的震后形变。本研究针对Gao等的软件进行改进,可计算震后任意时间点对应的震后形变。新软件由3个部分组成:1)与32个震后时间点相对应的32套离散格林函数数值框架;2)格林函数插值计算程序,可针对上述32套格林函数数值框架进行插值运算,输出任意震后时间点对应的格林函数数值结果;3)积分计算程序,调用上述格林函数数值结果,计算任意类型地震在地表任意位置产生的同震与震后形变。一般情况下,使用者只需按要求准备辅助文件,提供发震断层模型和观测站位置信息,以及震中周围地区地幔粘滞性因子,先后运行格林函数插值计算程序和积分计算程序,即可计算出目标地震在地表任意位置产生的同震与震后形变。本文基于粘弹球体位错理论与弹性球体位错理论,分别计算2011年日本MW9.0地震引起的远场同震位移,2套结果的高度一致性证明了新程序的正确性。最后,介绍需要注意的若干事项,便于使用者掌握该软件。

震后GPS坐标时序中对数弛豫时间的估计

全球定位系统(global positioning system,GPS)坐标时序去除同震形变和震前稳态速度场后,采用加权非线性最小二乘估计震后对数弛豫时间,可更准确地提取震后对数弛豫项,从而可以分析震后弛豫项对测站位移的独立物理贡献,并为震后余滑和黏滞性松弛效应等现象的分析提供参考。以日本2009—2019年GPS坐标时序为例,估计2011年Mw 9.0地震震后对数弛豫时间,发现不同站点的对数弛豫时间与其震中距关系显著,且服从高斯分布。据此,构建高斯函数加常数模型,可由震中距概略估计震后对数弛豫时间。高斯分布曲线的峰值、峰值位置、半宽度信息、最低位置分别为3.5 a、0 km、262 km、0.5 a,由此得出震后对数弛豫项影响时间大于0.5 a的站点主要集中在震中距约524 km范围内。震后弛豫效应区域分布的差异性显著,对数弛豫时间越长的区域,弛豫项水平位移表征越大,其中存在两个平均弛豫时间2.5 a的中心区域,与震后余滑的中心区域及时间相吻合。