2011年日本MW9.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¹⁸ Pa·s和30 km;然后,从震后3年内GPS观测数据中剔除地幔黏滞性松弛效应,获取了断层余滑对应的震后位移场;最后,利用基于球体位错理论的反演算法,反演了2011年日本M_W9.0地震断层余滑的时空演化过程.结果表明,2011年日本M_W9.0地震引起的断层余滑在震后半年内变化显著,震后2年主震区域余滑基本停止,断层的两端存在一定的余滑效应,断层余滑的累计矩震级达到8.59;地震后4年,地幔黏滞性松弛效应对震后位移场的贡献在总体上超过断层余滑的贡献.     

利用GRACE观测数据研究苏门答腊区域的黏滞性结构

本文尝试利用卫星重力观测资料和震后黏弹性松弛理论研究苏门答腊地区的区域流变结构,为更好地认识区域地球动力学环境提供依据．利用GRACE卫星重力资料,计算了2004年苏门答腊Mw9．3地震的同震及震后的重力变化．计算中使用平滑半径为500km的高斯滤波器．结果显示苏门答腊Mw9．3地震破裂东侧以陆地为主的上盘同震下沉量很大,造成约9x10＾-8mS＾－2的重力下降阶变,而西侧处于海水下的下盘重力同震上升较小约2x10＾－8mS＾－2,但其震后上升较快．流变结构对岩石的变形有很大的影响,是地球动力学数值模拟取得可靠科学结果的基础．本文尝试了利用卫星重力变化时间序列来反演苏门答腊地区的黏滞性结构;即基于GRACE时变重力场,利用自重力、黏弹性、平面分层模型模拟了该地震的同震和震后变形,并将获得的空间固定点的重力变化与GRACE重力场及点位时间序列进行比较,估计该地区的黏滞性系数在1．0x10＾18PaS的量级,且断层两侧的流变参数存在差异．最后结合苏门答腊区域的构造特点讨论了黏滞性系数的影响因素．     

震后GPS观测数据揭示的日本M_W9.0地震周边地区地幔黏滞性结构垂向变化

本文利用大范围的震后GPS数据和黏弹性球形地球位错理论,定量研究了日本M_W9.0地震周边地区地幔黏滞性结构的垂向变化.首先结合陆地和海底的GPS观测数据,以及基于球形地球位错理论格林函数和贝叶斯反演方法,反演了该地震的同震滑动分布,发现其最大错动量高达59m.然后在均一地幔黏滞性结构的假设前提下,确定了震源周边地区地幔黏滞因子的最优解,发现依据该地幔黏滞因子获得的理论远场震后位移和GPS观测结果之间的均方根误差高达0.81cm,不能解释远场观测结果.为解决上述问题,本文对震中周边地区地幔黏滞性结构沿垂向方向进行分层,建立了一个随深度变化的地幔黏滞性构造模型,然后综合利用远近场的GPS数据对该地区地幔黏滞因子进行反演研究,结果表明,震源周边地区岩石圈弹性层厚度最优解为40km,40～220km深度的地幔黏滞因子最优解为6×10~(18)Pa·s,220～670km深度之间的地幔黏滞因子最优解为1.5×10~(19)Pa·s.上述地幔黏滞性构造使远场的均方根误差降为0.12cm,仅为利用均一地幔黏滞性构造所得均方根误差值的15%,大大提高了远场模拟结果的准确性.最后,观测值和模拟值之间的均方根误差分析表明,近场震后形变数据主要约束浅层的地幔黏滞性结构,而远场震后形变数据主要约束深部的地幔黏滞性结构.     

基于2001年MW7.8可可西里地震震后形变模拟研究藏北地区岩石圈流变学结构

青藏高原岩石圈的流变学结构和形变机制是地学界长期争论的重大科学问题.2001年发生在东昆仑断裂带的M_W7.8可可西里地震造成青藏高原北部地区岩石圈构造应力场的很大改变,引起下地壳与上地幔的快速弛豫形变,从而为研究这一问题提供了难得的机会.本研究采用该区域的GPS震后观测,反演这一地区岩石圈的流变学参数并探讨其形变机制.反演所采用的数据来自45个GPS观测点,其中包括一个中国地壳运动观测网络的基准站,数据最长时间跨度达6.4年.大地震震后形变场主要来源于地壳、上地幔的黏弹性松弛与断层面上的震后余滑,因此本研究同时反演介质的黏滞系数和断层的震后余滑.考虑到东昆仑断层南侧的巴颜喀拉-羌塘地区与北侧的柴达木盆地地区具有明显不同的地壳结构,断层南北两侧采用不同的Burgers体流变学结构,其下地壳-上地幔的短期和长期黏滞系数采用网格搜索法获得;断层震后余滑反演则同时施加近似正比于库仑应力的约束.最终结果显示：东昆仑断层北侧柴达木盆地地区下地壳-上地幔短期和长期黏滞系数分别为5×10¹⁸ Pa·s和1.5×10²⁰ Pa·s;东昆仑断层南侧巴颜喀拉-羌塘地区下地壳-上地幔短期和长期黏滞系数分别为1.5×10¹⁸ Pa·s和1.5×10¹⁹ Pa·s.这一结果表明：巴颜喀拉-羌塘地区下地壳-上地幔黏滞系数显著低于柴达木盆地,意味着巴颜喀拉-羌塘地区下地壳可能存在部分熔融,其地壳形变模式更趋近于连续形变,而柴达木盆地形变模式更趋近于块体运动.研究区下地壳长期黏滞系数比下地壳流模型所主张的黏滞系数高2~3个数量级,表明下地壳流在本地区可能不存在.     

利用GRACE数据反演日本Mw9.0地震区域黏滞性

利用CSR(Center for Space Research)发布的GRACE RL05月重力场模型数据,通过水文模型GLDAS(Global Land Data Assimilation System)和CPC(Climate Prediction Center)扣除土壤水和雪水的影响,根据冰川模型扣除GIA(Global Isostatic Adjustment)的影响,采用P3M6去相关滤波和300 km扇形滤波,基于最小二乘拟合的方法得到日本M_W 9.0地震的同震及断层上下盘两个特征点重力变化时间序列,利用PSGRN/PSCMP模型对日本M_W 9.0地震区域黏滞性进行了反演,并计算了同震及震后5年研究区域重力变化的空间分布.结果表明,扣除土壤水和冰川均衡调整因素的影响,同震重力变化为-5.2×10~(-8)～2.9×10~(-8) ms~(-2);两个特征点在震后5年重力均增加,下盘重力增加较大;日本M_W 9.0地震区域黏滞性横向差异较明显,断层上下盘的地幔黏滞系数分别为2.5×10~(18) Pa·s、5.0×10~(17) Pa·s时,与GRACE观测值较接近,综合考虑断层上下盘的震后重力变化,区域黏滞系数大约为1.5×10~(18) Pa·s.     

利用震后黏弹性位错理论研究苏门答腊地震(M_w9.3)的震后重力变化

本文利用2003—2011年的GRACE RL05数据提取了苏门答腊地震(Mw9.3)引起的震后重力变化,发现断层两侧震后重力变化速率存在明显差异,断层下盘总体变化率为0.55μGal/yr,断层上盘为0.16μGal/yr.基于子断层叠加的编程思想,本文将Tanaka的黏弹球体位错理论配套计算程序(简称黏弹位错程序)加以改造,克服了其近场计算精度不足(甚至错误)的缺陷,可用来研究大地震引起的近场震后位移与重力变化.本文利用改造后的黏弹位错程序计算了2004年苏门答腊地震(Mw9.3)产生的同震重力变化,计算结果在空间分布和量级上均与利用弹性位错程序计算获得的结果一致,验证了我们对黏弹位错程序进行改造的正确性.最后,结合GRACE卫星观测数据,本文利用Tanaka的黏弹位错理论研究了苏门答腊地区的地幔黏性因子.结果表明,该地区地幔黏滞性具有显著的横向差异,当发震断层上下两盘的地幔黏滞性系数分别取8×1018 Pa·s和1×1018 Pa·s时,模拟的震后重力变化在总体空间分布和变化趋势上与GRACE卫星观测结果更接近.     

利用同震和震后位移数据联合反演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地震断层滑动分布的可靠性。 最后, 本文提出的反演方法为同震观测结果缺乏的大地震震后科考提供了理论支撑: 在大地震发生之后, 即使在同震期间没有足够的观测数据, 也可以在震后通过对震源区的加密观测积累的震后数据, 使用本文提出的反演方法优化同震断层滑动模型。     

强震后地表变形的动力学机制研究—以1999年台湾集集地震为例

强震后地表变形的动力学机制是地球动力学研究的重要方面,现在普遍认为震后变形主要由断层的震后余滑或由介质的黏弹性松弛所至。1999年台湾集集地震GPS观测系统记录到了空前的资料,为研究震后变形的动力学机制提供了难得的机会。本研究认为集集地震后地表变形由震后断层余滑、下地壳／上地幔的黏弹性松弛、震源区介质的破裂、孔隙弹性回跳、地下流体的运移、介质孔隙度及孔隙压的变化等多种因素共同影响决定。为抓住重点,研究中将介质的破裂、地下流体的运移和孔隙弹性回跳等因素等效为震源区介质的物性变化。文中运用黏弹性有限单元模型（麦克斯威尔体）、利用GPS观测的时间序列资料对震后余滑、地壳／地幔黏度以及等效的震源区介质物性变化进行了反演。反演模型给出了震后余滑的分布及变化特征,反演结果初步显示台湾地区的下地壳／上地幔的黏度分别为2.7×10^18,4.2×10^20Pa·s。此外,反演结果还给出每种影响因素对地表变形的贡献大小,在集集地震后的450d时间里,断层的震后余滑引起的地表变形占总变形的44.6％,下地壳／上地幔的黏性松弛占34.7％,等效的震源区介质的物性变化占20.7％。     

盲逆断层震后变形耦合机制研究——以2017年伊朗MW7.3地震为例

震后变形是地震周期中一个重要阶段,深入研究其产生机制有助于提高对地震演化过程、断层性质及地震危险性评估的认识.2017年伊朗M_W7.3地震InSAR观测资料丰富,为研究扎格罗斯造山带山前褶皱带盲逆冲断层震后变形机制提供了切入点.合理详细的震源模型及高精度的形变观测是开展震后余滑或黏弹松弛研究的首要前提,本文采用多视角InSAR资料联合远场波形数据反演该地震的破裂滑动分布,并解算得到伪三维地表位移,利用InSAR时序分析提取震源区主震后一年半内变形特征,显示震源区震后形变特征与同震位移场类似,LOS向形变速率最大值约8 cm·a^-1,震后变形预示震源区可能存在明显的震后应力调整现象.利用震后半年内形变资料约束的纯运动学余滑模型表明同时期释放能量为矩震级M_W6.7,进一步探讨应力驱动震后余滑及下地壳黏弹性松弛对震后变形的贡献,基于分层黏弹模型的模拟计算表明震后余滑和黏弹松弛效应的耦合模型可以更好地解释地表震后变形特征,其中震后余滑主要分布在同震破裂区的上方浅部,对主震后一年半内震后地表变形起主导作用,震源区壳下黏滞系数量级下...     

基于GPS观测研究2010年青海玉树MS7.1地震震后地壳形变特征及其机制

2010年4月14日青海玉树M_S7.1地震发生在青藏高原东南部甘孜-玉树地震带,在震后7~10天内,我们快速建立了由15个GPS测站组成的跨地震破裂带观测剖面,包括1个连续站,3个半连续站和11个流动站,对所有站进行了240多天的观测,获取了该次地震的震后形变时空特征.采用欧拉矢量和位错模型解算了背景速度场,并从GPS观测的形变场中扣除该分量.采用分层黏弹性位错模型计算余震引起的地表形变,结果表明余震对部分测站的位移造成不可忽视的影响.采用对数模型拟合位移时间序列,表明特征衰减时间为6.7±1.2天.利用最速下降法反演震后余滑时空分布,反演结果表明震后断层活动以左旋滑动为主,断层南盘具有少量的抬升.在空间分布上,余滑主要位于同震破裂区的两侧,西北侧的余滑几乎达到地表,而东南区的余滑基本在同震破裂区的下方,余滑最大的区域位于结古镇东南下方10~20 km的深度范围.随着震后离逝时间的增加,2个余滑区在空间上保持不变,余滑区的面积逐渐扩大.余滑的矩释放为（1.5~5.1）×10¹⁸Nm,相当于1个M_W6.1~6.4地震释放的能量.分层岩石圈黏弹性模型计算的地壳孔隙弹性反弹形变与地表观测值相差较大,不能解释观测到的震后变形.采用麦克斯维尔流变体模型计算下地壳和上地幔松弛引起的地表形变,显示出其对地表形变的贡献较小.GPS观测得到的震后形变所具有的快速衰减特征,以及余滑模型能够较好地拟合GPS地表形变,表明2010年玉树M_S7.1地震后早期阶段的地壳形变主要是由余滑机制决定的.     

GRACE观测的区域重力场变化与大地震关系探讨

利用GRACE重力卫星观测资料, 系统地分析了喜马拉雅俯冲带及周边几次7级以上大地震前后的区域重力场变化特征。 研究发现, 2005年巴基斯坦M_S7.8、 2008年汶川M_S8.0、 2011年缅甸M_S7.2、 2013年芦山M_S7.0和2015年尼泊尔M_S8.1等大地震都发生在重力显著变化区域及周边, 地震前几年开始区域重力在原变化趋势下呈现快速减小特征, 而后, 减小后的重力恢复增加, 大地震均发生在这种重力变化过程中。 大地震前的这种重力变化过程与苏门答腊M_W9.3和东日本M_W9.0地震的同震和震后区域重力变化过程类似, 只是幅度要远远小于这两个地震的同震和震后变化, 而且重力减小恢复较快。 综合分析认为喜马拉雅俯冲带及周边区域几次7级以上大地震前后的区域重力变化, 都与巨大地震过程中, 板块活动调整的弹性变形、 区域地下深部地幔物质运移和热引起的区域气候改变有关。 据此, 为区域重力变化的共同特征提出一种物理过程模式解释。 苏门答腊M_W9.3及其巨大余震过程中, 喜马拉雅俯冲带局部区域出现地壳的拉伸、 下沉, 一段时间后, 转变为挤压、 抬升。     

2008年10月西藏当雄MW6.3地震震后形变提取与余滑反演

2008年10月6日西藏当雄发生M_W6.3地震.本文利用震后2008年10月26日至2010年8月22日的16期ENVISAT ASAR数据,通过小基线集干涉测量、误差校正与MInTS（Multiscale InSAR Time Series）技术提取高精度的震后形变场,利用SDM（Steepest Descent Method）方法反演断层震后余滑演化过程,并分析震后余滑与同震滑动的关系.结果表明：当雄M_W6.3地震的近场震后形变场主要位于断层西侧,在时间演化上具有明显的对数函数衰减规律;震后余滑主要集中于断层中南段深0~15 km区间,最大的余滑量约0.07 m,位于断层深约9.28 km处,滑动角约-103°;震后余滑引起的地震矩能量M₀与矩震级M_W在时间演化上具有指数函数递增规律;当剪切模量μ=32 GPa,震后665天余滑释放的地震矩能量约为1.92×10¹⁷N·m,约占同震滑动释放地震矩的4.8%,相当于矩震级M_W5.46;虽然震后余滑已经延伸到断层浅部0~5 km区间,但由于余滑量相对较小,没有改变同震滑动在断层浅部区域的滑动亏损现象,这可能是2010年11月30日该区域又发生M_W5.3级余震的主要原因之一.     

2011年日本MW9.0地震对沂沭断裂带及其两侧地区地壳运动的同震影响研究

2011年日本M_W9.0地震（简称日本地震）后沂沭断裂带及其两侧地区地震活动显著增强,研究日本地震对该地区地壳运动及地震潜势的影响十分必要.为此,本文通过112个连续GPS观测站获取了研究区高空间分辨率的日本地震同震形变场并得到如下认识：（1）8个定点地球物理观测的同震响应验证了本文同震形变场的可靠性;日本地震的东向拉张使研究区整体上处于张性同震应变状态,但存在局部挤压区域,其中莱州湾至海州湾的挤压条带穿过沂沭断裂带并对断裂带南北两段产生了不同的同震作用,对南段具有拉张作用,对北段产生挤压作用;（2）同震形变场在鲁东隆起和鲁西断块产生了显著的剪应变,地震b值显示上述区域的构造应力在日本地震后增强,因此同震形变场可能改变了这些区域的应力特征;（3）地震矩张量叠加分析显示,同震形变场短期内对鲁西断块、鲁东隆起区和沂沭断裂带南段累积了地震矩,可能有助于上述区域在日本地震以后的地震活动增强;日本地震对沂沭断裂带北段的地震矩具有释放作用,或许是该区域地震活动减弱的原因.     

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

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

Seismologic applications of GRACE time-variable gravity measurements

The Gravity Recovery and Climate Experiment (GRACE) has been measuring temporal and spatial variations of mass redistribution within the Earth system since 2002. As large earthquakes cause significant mass changes on and under the Earth’s surface, GRACE provides a new means from space to observe mass redistribution due to earthquake deformations. GRACE serves as a good complement to other earthquake measurements because of its extensive spatial coverage and being free from terrestrial restriction. During its over 10 years mission, GRACE has successfully detected seismic gravitational changes of several giant earthquakes, which include the 2004 Sumatra–Andaman earthquake, 2010 Maule (Chile) earthquake, and 2011 Tohoku-Oki (Japan) earthquake. In this review, we describe by examples how to process GRACE time-variable gravity data to retrieve seismic signals, and summarize the results of recent studies that apply GRACE observations to detect co- and post-seismic signals and constrain fault slip models and viscous lithospheric structures. We also discuss major problems and give an outlook in this field of GRACE application.