联合InSAR和GPS研究鲜水河断裂带炉霍—道孚段震间运动特征

鲜水河断裂带是位于青藏高原东侧的大型走滑断裂,地震活动性高,其炉霍—道孚段是少数观测到蠕滑现象的段落之一,但目前尚缺乏高分辨率大地测量数据的约束研究.本文综合利用InSAR和GPS观测数据研究炉霍—道孚段的运动特征,提出一种适用于研究震间形变的InSAR和GPS数据融合方法,将InSAR速度场和GPS插值速度场转换至平行断层方向,利用抽样统计方法去除GPS参考框架影响,在数据融合过程中,注重保留远场构造运动特征和近场高分辨率形变信息.建立包括浅层蠕滑、中间闭锁单元和深部滑动的断层运动模型,以融合的形变速度场为约束,利用马尔可夫链蒙特卡罗方法反演炉霍—道孚段断层滑动速率、浅层蠕滑特征和闭锁状态.InSAR速度场剖面分析结果表明：炉霍—道孚段地表蠕滑特征明显,沿走向方向存在差异,炉霍段蠕滑速率明显低于道孚段,表现出从北西至南东速率逐渐增大、减小、再次增大的变化特征,断层左旋蠕滑速率在1.0~5.2 mm·a^-1之间;断层运动模型反演结果表明,炉霍—道孚段蠕滑深度为~9.8 km,浅层蠕滑速率为3.27~4.18 mm·a^-1,断层深部滑动速率为8.12~9.30 mm·a^-1.炉霍段地震复发周期为370~410年,闭锁层厚度~1 km,地震危险性较低.道孚段地震复发周期为59~65年,存在厚度~7.5 km的闭锁区,可积累足够的应变,浅层蠕滑释放的能量为~1.32×10¹⁶ Nm·a^-1,占构造运动加载积累能量的23%~38%,1981年道孚6.9级地震以来,闭锁区近40年积累的地震矩达到3.62×10¹⁷至7.88×10¹⁷ Nm,相当于一次M_W5.6~M_W5.9地震的能量,其地震危险性值得关注.

     

����GPS���й۲�������ȡ������Ϣ�IJ����۽���

???????????GPS???и?????????????????沨????????????????????GPS???????е????????м????????????????????GPS???й?????????????????????????????????????????????????????????????????????????1.5????????????????????;????????3?,?????????????????????     

����GPS�����о���γ�Ȼ�����д��Ļ��

????GPS????о????γ???????д???????????????????????????????????????????????????????????????????????????????????????????????????????????13.7??0.6??5.8??0.5??7.1??1.0??12.7??0.2??25.4??0.5 mm/a???????????Σ?????????????????????????????????????29.2??1.8??14.1??1.0 mm/a?????????????????????????????????????????18.5??1.0??15.2??0.6 mm/a?????????????????????????????10.7??1.1 mm/a????γ????д????????????????????????л??????     

位错引起的大地水准面变化

从动态的观点介绍了大地水准的物理意义。用数值分析方法计算了常见的几种断层位错导致的大地水准面局部变化的分布图像,并着 断层位错的倾角与大地水疰同变化 之间的关系。     

GPS测站位移速率的解算方法

对目前4种常用的利用GPS复测资料解算测站位移速率的方法进行了分析，通过实例对比了4种方法的结果。坐标法可以直接利用多期复测结果得到测站大地坐标，但要求坐标在同一参考框架内，并采用相同的参考椭球；基线方法可以消除或削弱观测系统的系统误差，与参考框架无直接关系，但在平差计算时至少要引入1个起算点，起算点高程方向的误差会削弱其他点平面坐标的精度，并进而影响位移速率的精度；在选定不同基准解算位移速率时相似变换方法是十分有效的，但如果复测结果不在同一参考框架中，其结果存在框架不一致所带来的误差；QOCA方法能够综合利用多种观测数据，可以把区域观测数据和全球观测数据稳健地结合起来，在削弱观测期间突变因素的影响方面比其他方法具有便利之处。     

利用 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陆态网络观测到的同震形变.     

����TRACK����GPS����Ԫ��λ

????GPS??????????GAMIT?е???????λ??飨TRACK??????1 Hz??3 s??????GPS?????????????ж?λ??????????????????????????????????????????????????????30 s??????1Hz???????????????????С??10 km??10~100 km??????????100 km 3????????д??????λ??????????????????????????????????????????￡??????λ???????????????????5 mm??     

2016年青海门源MS6.4地震震前应变积累及同震变形特征

2016年1月21日青海省门源县发生了M_S6.4地震, 发震断裂为冷龙岭北侧断裂, 震中位置与1986年门源6.4级地震相同。 本文收集了本次地震震中及其周边区域1999—2015年GPS观测资料, 解算了GPS速度场、 跨断裂连续观测站基线时间序列和应变率场。 结果显示, 祁连山断裂带为一条宽度约60 km的连续变形带。 在断裂带南侧地壳运动以顺时针旋转为主, 运动量值没有显著差异; 跨过断裂带到达其北部之后, 地壳运动量值明显减小, 显示出该断裂带的强烈活动特征。 冷龙岭断裂左旋走滑速率为6.17±0.41 mm/a, 挤压速率为1.83±0.38 mm/a, 断层闭锁深度为22.1±3.1 km。 利用GPS连续观测站数据解算的地震同震位移显示, 震中西南侧26.8 km处的青海门源(QHME)测站记录到了明显同震位移, 其水平运动方向为北东向, 与逆冲为主的震源机制解一致。     

青藏高原东北缘地壳应变与应力特征研究

青藏高原东北缘是青藏高原隆升和变形的前缘,地壳变形剧烈,研究其地壳运动和变形对理解该区的构造活动特征和动力学机制具有重要的意义。本研究收集1991—2016年的GNSS速度场数据,采用多尺度球面小波方法计算应变率张量,分析主应变率、面应变率和最大剪应变率的空间分布特征。面应变率结果显示青藏高原东北缘大部分区域具有轻微的压缩,平均压缩率小于15 nstrain/a,压缩率较高的区域集中在青藏高原东北缘的边缘地区,其中祁连山断裂带和海原断裂带区域的面压缩率大于20 nstrain/a。东昆仑断裂带、祁连山断裂带、海原断裂带和六盘山断裂带具有较高的剪应变率,最大约为40 nstrain/a。研究还收集青藏高原东北缘1904—2021年2.0级以上地震的震源机制解,通过区域阻尼应力反演方法得到该区域的应力场,最大主应力方向显示青藏高原东北缘35°N以南的区域表现为近EW向的挤压,而35°N以北的区域表现为NE向的挤压。根据最大、最小主应力的倾角把研究区划分为正断、逆冲和走滑3类区域,发现青藏高原东北缘内部的大部分区域和六盘山断裂的北部区域表现出走滑的运动特征,而祁连山断裂带、海原断裂带和西秦岭...     

利用地基GPS观测数据研究与地震有关的电离层异常

利用中国大陆GPS连续站观测数据, 根据电离层单层模型, 计算得到中国大陆上空电离层电子总含量TEC(Total Electron Content)值的时间序列。 选取2001—2008年发生在中国大陆地区的11个M_S≥6.0地震, 分析了这11个地震前后TEC值时间序列的变化特征。 分析结果显示, 地震活动期间距离震中700 km以外的GPS连续观测站上空TEC值的异常现象不明显, TEC值变化基本在0.5 TECU之内。 11个震例中有6个震例TEC值出现了负异常现象, 主要表现在震前1～5 d, 并且观测站距离震中越近, TEC值变化越明显, 少数震例在地震发生期间TEC值明显增加。 研究认为, 利用GPS连续观测数据研究与地震有关的电离层前兆是可能的。