基于GRACE数据的尼泊尔MS8.1地震北向重力梯度变化

利用德克萨斯大学空间研究中心发布的GRACE RL05月重力场模型数据,采用300 km扇形滤波,得到2015年尼泊尔M_S8.1地震北向重力梯度的时空分布,之后利用最小二乘拟合方法分析尼泊尔及其邻近区域北向重力梯度的长期变化趋势及研究区内6个特征点的北向重力梯度月变化时间序列,并结合黏弹性分层位错理论对GRACE检测尼泊尔M_S8.1地震北向重力梯度同震变化的可能性进行分析. 研究结果表明:在尼泊尔M_S8.1地震发生前北向重力梯度表现出比较明显的正负异常变化,而该地震即发生在北向重力梯度正负变化的零值线附近;研究区北向重力梯度年变率在印度板块边界及其垂直方向所形成的四象限呈正负相间分布,6个特征点的北向重力梯度在2012年4月出现比较明显的跳变. 由此推断,北向重力梯度的动态变化过程反映了震前区域物质迁移和震后壳幔物质黏滞性调整等问题.      

融合GOCE和GRACE卫星数据的无约束重力场模型Tongji-GOGR2019S

本文在法方程层面融合GOCE卫星的V_xx、V_yy、V_zz和V_xz重力梯度分量观测数据和GRACE卫星观测数据,采用直接法解算了220阶次的重力场模型Tongji-GOGR2019S.首先利用ⅡR带通滤波器在5~41 mHz的重力梯度带宽范围内对约24个月的GOCE重力梯度观测方程进行无相移滤波处理,并组成解算220阶次重力场模型的法方程,各梯度分量根据相对于参考模型统计精度进行定权;然后与13.5 a GRACE数据建立的180阶次Tongji-Grace02s重力场模型的法方程进行叠加,解算了220阶次的无约束纯卫星重力场模型Tongji-GOGR2019S.利用EIGEN-6C4重力场模型、GNSS/水准数据、DTU15重力异常数据以及欧洲区域似大地水准面模型EGG2015等数据对Tongji-GOGR2019S模型精度进行全面的检核评定,结果表明：引入GOCE卫星梯度数据后,高于72阶的位系数精度优于Tongji-Grace02s模型,Tongji-GOGR2019S模型的整体精度接近同阶次的DIR-R6等GOCE卫星第6代模型.     

2016年台湾美浓MW6.4地震震源参数的InSAR和GPS反演

2016年2月6日台湾西南部高雄市美浓区发生了M_W6.4地震.本文结合ALOS2卫星升降轨、Sentinel-1A升轨SAR数据,采用两轨差分干涉技术获取了该区域的同震形变场,形变结果表明震中西北部以抬升为主,最大视线向形变量约为11.2 cm.基于均匀位错模型和多峰值粒子群（MPSO）算法,利用InSAR和GPS形变数据联合反演了美浓地震的断层几何参数,结果表明震源中心位于22.920°N,120.420°E,深度约12 km,发震断层长度约15 km,走向角307°,倾角16.5°,平均滑动角为51.5°,此次地震是以逆冲倾滑兼左旋走滑的破裂模式.利用格网迭代搜索法得到最优倾角为15.7°,GPS和InSAR最优权比为18：1,最优平滑因子为0.06.基于非均匀位错模型,利用非负最小二乘方法进行线性反演,结果显示最大倾滑和走滑量分别为51.7 cm和55.3 cm,对应矩震级为M_W6.38,略小于GCMT （M_W6.4）的结果.通过与已有文献的比较和对该区域断层构造的分析,发现美浓地震的发震断层为单一断层的解释更为合理,我们推测发震断层是位于左镇、后甲里等断层之间的一条东南-西北走向往东北倾斜的盲断层,并初步推测2010年M_W6.3甲仙地震也同该断层有关.     

2008年新疆乌恰Mw6.7地震震源机制与形变特征的InSAR研究

2008年10月5日新疆乌恰M_w6.7级地震发生在南天山、帕米尔高原及塔里木盆地交汇地带,基于地震波反演的震源机制解确定的震源深度存在较大差异.本文利用日本ALOS卫星的PALSAR图像,获得了本次地震的同震形变场,基于卫星视线向（LOS）和方位向（Azimuth）的形变,采用均匀弹性半无限位错模型和有界最小二乘（BVLS）算法,以网格矩形位错元法对发震断层的几何产状、滑移及分布进行了估算,结果表明本次地震以逆断破裂为主,断层面上最大位错量接近3.4 m,形变中心位于73.8040°E,39.5335°N,深度约5 km,震级估算为M_w6.6;地震发生在走向46°,倾角48°的断层上,发震断层长30 km,宽14 km,闭锁深度9 km,符合该地区浅源地震多发的构造特点,发震断层为乌合沙鲁断裂带.InSAR反演的滑移形变主要集中于地下2~7 km,表明乌恰地震为浅源地震,可能与该断层附近历史地震未完全释放的残余应力积累有关.同时,InSAR反演的断层位错分布呈现双破裂特征,震级分别为M_w6.5和M_w6.1,可能与本次地震的主震和余震相对应,也可能是由主震激发而产生的两组破裂.     

智利M_W8.8地震同震重力梯度变化

利用GFZ Release 05卫星重力GRACE观测数据,计算2010年2月27日智利MW8.8逆冲型地震的同震重力和重力梯度变化,分析其分布特征,可知:由GRACE探测到的同震重力变化在断层俯冲区域可达-9.5μGal,断层隆升区域可达+3.5μGal,结果与利用SNREI地球模型的位错理论计算的同震重力变化较一致,说明利用GFZ Release05 DDK5滤波数据,更能精确的反映同震重力场变化;GRACE检测的智利地震同震径向重力梯度变化Trr最大可达-600μE,位于发震断层东侧俯冲区域;通过对同震重力梯度分布特征分析,初步判断发生同震物质迁移的区域范围在断层俯冲区域为(67°—72°W,33°—38°S),在断层隆升区域为(73°—77°W,35°—39°S)。     

利用走时和波形拟合迭代反演阿拉善地壳速度结构及2015年阿拉善左旗5.8级地震震源机制

据中国地震台网测定,北京时间2015年4月15日15时39分,在内蒙古自治区阿拉善左旗（39.8°N,106.3°E）发生M_S5.8地震,震源深度为10 km.地震发生后多家机构对其开展了研究,本文使用喜马拉雅Ⅱ期布设在南北地震带北段的台站观测数据,通过走时反演和波形拟合反演的迭代,获得了该地区地壳一维速度结构,接着利用直达P波观测与理论走时差对震中位置重定位,然后反演地震的最佳双力偶解以及震源深度,最终得到了区域速度结构、地震的三维坐标、发震时刻以及震源机制解.结果显示,此次地震发生于世界时2015年4月15日7时39分26.718s,震中（39.7663°N,106.4304°E）,震源矩心深度18 km,矩震级M_W5.25,节面Ⅰ走向176°,倾角85°,滑动角-180°,节面Ⅱ走向86°,倾角90°,滑动角-5°.结合该区域断裂带构造运动分析,本文认为此次地震是左旋走滑破裂,略带正断分量,断层面是节面Ⅱ,走向为NEE（近E-W）向,发震构造为震中附近的E-W向隐伏断裂.     

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级余震的主要原因之一.     

2017年伊朗MW7.3地震的震后余滑分布及其对2018年MW6.0地震的触发影响

收集了覆盖监测区域的Sentinel-1卫星雷达影像,利用短基线集干涉测量技术提取了2017年伊朗萨波尔扎哈布M_W7.3地震后283天的地表时序形变,通过二步法反演得到其震后余滑分布,之后采用差分干涉技术获取了2018年发生于同一地区的贾万鲁德M_W6.0地震的同震形变场,并将反演所得的发震断层参数作为应力计算的接收断层参数,来分析2017年M_W7.3强震及其震后活动对2018年贾万鲁德M_W6.0地震的触发影响。结果表明:萨波尔扎哈布地震的震后形变主要由孕震断层面的余滑运动所致,震后283天余滑模型的累积滑移量达到0.7 m;2018年贾万鲁德地震的发震断层走向为355.6°,倾角为89.4°,同震断层破裂以右旋走滑为主,兼具部分正断层运动。本文所得的贾万鲁德地震断层平面上的库仑应力变化表明,2017年M_W7.3主震及其震后余滑对2018年M_W6.0地震的发生具有一定的触发效应,M_W6.0地震的发生可能与区域板块的活动性相关。      

新疆天山地区强震前地震条带不清晰的原因初探

新疆地区在强震前地震条带不明显.本文对比分析了新疆天山与川滇地区区域构造背景和地震震源机制解,得到:①新疆天山地区区域地震构造以逆冲型为主,而川滇地区区域地震构造以走滑型为主.②1970年以来的天山地区中强地震以逆断层和走滑型地震类型为主,走滑型地震多数包含较大的倾滑分量;而川滇地区6级强震以走滑型地震类型为主,地震主要以水平错动为主.认为弱震条带适用于以水平剪切错动为主的地震,而新疆天山地区地震倾滑分量较大,因此中强地震前弱震条带图像不清晰,中短期地震预报意义不显著.     

芦山MS7.0级地震InSAR形变观测及震源参数反演

2013年4月20四川省芦山县发生M_S7.0级地震,目前的研究资料表明地震发生在龙门山断裂南段,但地表未发现明显破裂.本研究利用InSAR技术与Radarsat-2雷达数据,获取了芦山地震同震的部分形变场,结果表明,近场区域的LOS位移发生视线向隆升,量值在7 cm左右.随后利用弹性半空间的位错模型反演了断层面参数,综合反演结果及震源机制解最终确定了发震断层的初始模型,以形变场观测数据为约束,基于梯度下降法反演获得了断层面上的滑动分布,反演得到的矩震级为M_w6.45级,断层走向213°,倾角39°~43°,最大滑动位于地表以下约13 km深度位置,最大滑动量0.91 m,平均滑动角71°,整体上仍以逆冲滑动为主,兼具左旋走滑.推测在双石-大川断裂以东12 km处展布一条隐伏断裂,为本次的发震断裂.     

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地震破裂以正断为主, 且带有一定的左旋走滑。     

利用独立成分分析检测2004年和2012年印度洋地震的重力变化

2004年12月苏门答腊发生M_W9.3地震,造成巨大的质量重新分布.利用GRACE卫星月重力场数据计算了研究区域地面1°×1°网格点上的重力变化时间序列,采用主成分分析和独立成分分析两种方法,提取了重力变化的空间与时间特征,结果显示震中两侧区域的重力变化呈两极分布,其中东侧重力下降,西侧重力增加.相较于传统的主成分分析方法,独立成分分析能更好地从原始信号中提取地震的信息,能分解出具有显著阶跃变化的独立成分.除了2004年的苏门答腊大地震外,独立成分分析还分解得到了2012年Andaman地震的特征,与该次地震的空间特征与时间序列基本一致.相较于适合定量分析的多项式拟合方法,独立成分分析更适合大范围区域的定性分析,建议将两种方法相结合,取长补短,从而为GRACE地震监测提供一种更为客观、有效的方法.     

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。     

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及其巨大余震过程中, 喜马拉雅俯冲带局部区域出现地壳的拉伸、 下沉, 一段时间后, 转变为挤压、 抬升。     

联合GPS和GRACE观测研究日本MW9.0地震震后变形机制

利用GPS和GRACE观测数据研究了日本M_W9.0地震的震后变形特征.GPS观测显示,区域震后位移呈现随指数函数变化特征,变化速率符合大森公式的衰减特性;近五年的震后水平位移累积已达到东向60~165 cm,南向20~65 cm的量值,距震中较远站点已超过同震变化量,且震后变形仍然持续.GRACE观测到显著的震后重力变化,地震破裂两侧的重力变化总体均呈上升趋势,但海洋侧的变化速率较快.联合震后余滑和黏弹性位错理论对震后变形进行了模拟,探索了GPS和GRACE观测的综合应用方法.研究发现,综合考虑震后余滑和黏滞性松弛效应可以对日本地震的震后变形做出较合理的解释,震后初期余滑起主要作用,1至2年以后逐渐减弱,黏滞性松弛作用逐渐增强.在震后变形模拟和区域黏滞性结构反演中形成GPS和GRACE观测结合应用的方法,先基于震后GPS形变估算区域黏滞性结构,而后利用GRACE观测修正深部的黏滞系数,并综合利用这两种观测微调浅层黏滞系数,最终确定区域黏滞性结构.基于该方法反演了日本震源区的地幔黏滞性结构,地震断层破裂两侧的流变参数存在差异,大陆侧的地幔顶层黏滞系数在1.0×10¹⁹ Pa·s量级,而海洋侧的则略小于大陆的,在6.0×10¹⁸ Pa·s量级.     

2016年新疆阿克陶地震Sentinel-1 InSAR同震形变特征

2016年11月25日新疆克孜勒苏州阿克陶县发生M_W6.6地震。 本文利用合成孔径雷达差分干涉测量技术, 对Sentinel-1卫星获取的升、 降轨雷达数据进行了处理, 提取了该次地震的同震形变场, 并结合形变场特征与震源机制解, 采用梯度下降法反演发震断层的滑动分布。 结果表明, 升、 降轨LOS向同震形变场在发震断层两侧具有明显不同的形变特征, 主要形变区域分布在断层两侧, 升轨LOS向形变量可达-8.2 cm与11.2 cm, 降轨LOS向形变量可达-21.4 cm与13.1 cm; 反演的升、 降轨干涉形变场与InSAR测量值之间的残差得到有效控制, 大部分的残差介于±5 cm之间; 断层滑动分布主要集中于沿断层面深约2~18 km处, 最大滑动量位于沿断层面深约7 km处可达0.96 m; 平均滑动角约182.29°, 最大滑动处的滑动角约197.13°, 两个滑动分布中心的滑动角均接近180°, 表明阿克陶地震为一典型的右旋走滑破裂性事件; 当剪切模量取32 Gpa时, 反演的发震断层地震矩M₀可达9.75×10¹⁸, 相当于矩震级M_W6.60, 与地震波形反演结果一致。     

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

2016年11月13日新西兰凯库拉(Kaikoura)MW7.8地震同震位移和应力数值分析

2016年11月13日新西兰南岛北端凯库拉（Kaikoura）发生了M_W7.8大地震,造成了强烈的地表变形并引发大面积滑坡和海啸的发生.基于美国地质调查局（USGS）断层滑动模型,建立全球同震横向不均匀并行椭球型地球模型,计算了此次新西兰凯库拉大地震产生的同震形变和应力及库仑应力变化.初步计算结果表明：新西兰凯库拉M_W7.8地震造成断层上盘东北向抬升,下盘西南俯冲;引起发震区域同震位移较大,从凯库拉到坎贝拉（Campbell）以及首都惠灵顿（Wellington）整体上东北向抬升,最大同震水平位移1.2 m,垂直位移1.1 m.此次大地震释放了发震断层上积累的压应力,但增加了发震断层两端的挤压力;同时,同震剪应力变化增加了NE-SW向断层发生右旋滑动的危险性;采用此次地震发震断层参数计算得出的最大库仑应力变化增加区域集中在发震断层两端,可达到MPa量级.当分别采用新西兰北岛Awatere断裂系和南岛Wellington断裂系参数计算库仑应力变化时,发现新西兰北岛和南岛震中以南区域的库仑应力均增加,可触发部分余震的发生.     

基于InSAR的青海大柴旦地震三维同震形变场获取与震源特征分析

本文利用Envisat ASAR的升、降轨和宽幅数据,通过基于先验知识的最小二乘迭代逼近获取大柴旦2次地震的地表三维同震形变.结果表明,2008年MW6.3地震垂直向形变主要发生在断层南盘,以隆升形变为主,最大隆升量约10cm,北盘沉降量小于等于-1cm.东西向形变在南盘呈向东运动的特征,最大运动量约4cm,北盘向西运动,最大运动量约为-2cm.2009年MW5.8地震垂直向形变显示断层南盘抬升的特征,最大抬升量约27cm,北盘最大沉降量约-3cm.东西向形变表现为南盘向东运动,最大约10cm,北盘向西运动,约为-4cm.可以看出这两次地震均表现为逆冲为主,兼少量左旋走滑的震源特征.视线向结果无法判定同震形变的少量走滑特征,而地表三维分量可以有效地识别出少量左旋还是右旋走滑的震源特性.本文以视线向、垂直向、东西向形变量作为约束条件,利用Okada模型正演了2008年地震同震三维形变场.结果显示,采用逆冲兼少量左旋走滑的发震断层参数,视线向、垂直向、东西向正演结果与观测结果吻合.这也表明采用分解后的地表三维同震形变场可以有效地识别出发震断层的少量左旋走滑特征.