自重可压缩层状球体地球模型下地震变形算法的优化研究

本文利用龙格库塔法求解球坐标系下微分方程组,基于自重、可压缩和连续分层地球模型研究地震引起的地球变形,优化相关算法.首先对传统的两次同向积分法进行调整,对微分方程的齐次部分进行两次相向积分,分别由地心向上积分至震源和由地表向下积分至震源.接着利用震源间断条件求解震源上方和震源下方通解的组合系数,得到微分方程组的解矢量....     

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

���ڸ�˹�����Ȩ�������о����������ؿ��α�����

???????????????????????????????????????????????????????????????????????????????????о??????????????1999??2009???GPS????????о???????????????????????????????????????????α?????????????????????????????????150 km????????????????????????????滪???????????????仯???????????????????????????????????????????????-??????????????????????????Σ???????????????????????????б??Σ??????“???????”?????????????????α????????????????????????????????С??????α????????????????????????????????????????б???????????????????????????EW??????????????????????????????????????????к????????е?????????????????????????????????????????     

堪察加地区现今地壳运动与变形特征研究

利用俄罗斯堪察加地区1995～2005年的GPS观测数据,研究了该区现今地壳水平运动速度场特征.在球坐标系中解算了各应变率分量,分析了应变率场的空间分布特征,并与地震学和地质学研究结果进行了综合对比分析.结果表明,堪察加半岛北部的微板块边界并不明显,堪察加南部测站运动速度大于中部和北部地区,愈靠近东部板块汇聚区,测站速度越大.从东海岸到西海岸,测站水平速度存在明显的梯度衰减特征,水平运动方向与太平洋板块向西北的俯冲方向基本一致.各应变率分量具有东部海岸大于中部和西海岸、从东至西呈梯度衰减的特点.堪察加大部分地区处于EW和NS向压缩状态,局部存在拉张.面应变率结果显示绝大部分为压缩区;刚性转动结果表明大部分地区表现为顺时针转动,北部地区和南端顺时针旋转性明显.东部有效应变率明显大于西部地区,东西向梯度衰减关系明显.主压应变率明显大于主张应变率,特别是在东海岸地区.主压应变率方向与中等以上地震的主压应力轴在水平方向的投影方向基本一致.地壳变形场在空间分布上的不一致性主要与太平洋板块在堪察加半岛东南侧的俯冲深度、俯冲方位角、俯冲倾角和俯冲带的耦合强度有关.     

利用地基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连续观测数据研究与地震有关的电离层前兆是可能的。     

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, 与地震波形反演结果一致。     

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)测站记录到了明显同震位移, 其水平运动方向为北东向, 与逆冲为主的震源机制解一致。     

中国地震科学实验场BDS-3定位精度和地壳运动初步分析

我国BDS-3卫星导航系统于2020年7月建成,并为全球用户提供服务,自此开始,中国地震科学实验场(简称“实验场”)的GNSS测站相继接收BDS-3卫星数据,目前已经积累了2年多的观测数据,为BDS-3应用于川滇地区地壳运动监测提供了重要平台和数据保障。为评估实验场BDS-3观测数据的精度及其应用于地壳监测的可行性,本文首先根据各测站所记录的BDS-3的B1I和B3I两个频点的多路径效应和信噪比与卫星高度角之间的关系,对BDS-3的数据质量进行评价。然后,基于高精度数据处理软件GAMIT/GLOBK (10.7版)对各测站同期记录的BDS-3和GPS数据分别进行处理,得到BDS-3和GPS坐标时间序列。基于包括线性项、年周期项和半年周期项等分量的函数模型,利用极大似然法分别对BDS-3和GPS时间序列进行拟合估计,得到各测站的线性速度、年周期和半年周期信号,并对拟合结果对比分析,评定BDS-3和GPS的定位精度。最后,讨论BDS-3定位精度的影响因素和BDS-3水平速度场的区域性特征。研究表明：实验场BDS-3与GPS原始数据质量相当,BDS-3坐标时间序列的均方根残差(RMS)的平均...     

ITRF2000无净旋转对中国大陆地壳形变的影响

通过计算全球14个1级板块的总角动量以及ITRF2000速度场相对NNR-NUVELIA的旋转分量，研究了ITRF2000的无净旋转(也称无整体旋转)符合特性。结果表明．虽然ITRF2000的速度场通过50个核心站归算到NNR-NUVELlA上．然而在研究中国大陆区域地壳形变时，它们之间的差异仍是不可忽视的，其差异性表现为：ITRF2000的速度场相对于NNR-NUVELIA模型的速度场表现为顺时针旋转．前者速度模量大于后者；两个速度场的差异量表现为东部大于西部，E-W向差异量最大值达到5mm／a以上，N-S向差异量接近4mm／a。提出了2种消除这两个速度场之间系统差异的方法。     

数字化地震仪系统的标定

介绍了数字地震仪系统传递函数的特征和测定方法,给出了地震仪标定软件的使用方法和计算步骤。结合实际波形数据处理过程,着重对标定计算中参数的设置,以及需要注意的问题进行了讨论和分析。这对于正确快速地获取地震仪的传递函数具有重要的现实意义,有利于国家数字地震台站和区域地震台站地震观测系统的定期标定工作,同时也将促进数字化地震观测资料的应用研究。