苍梧M_S5.4地震和呼图壁M_S6.2地震震前重力非潮汐变化特征

分别对苍梧M_S5.4地震震前的梧州台和呼图壁M_S6.2地震震前的库尔勒台连续重力观测数据进行分析研究。结果表明,2次地震前梧州台和库尔勒台分别观测到10个月和7个月左右的重力非潮汐持续性上升变化,月均变化速度分别为5.0μGal/月和5.1μGal/月,累计变化分别约为50μGal和35μGal。经现场核实和理论计算排除环境干扰影响,认为2次地震震前数月的重力非潮汐持续上升变化与孕震过程有关,可为构建连续重力观测地震预测指标体系提供重要参考。     

沈阳连续重力非潮汐变化异常分析

沈阳台定点重力非潮汐变化数据稳定,但自2017-06以来数据曲线连续5个月明显上升,且速率较大,幅度约为5μGal/月,打破了自2015年断电后重归的稳定状态。通过对仪器观测系统、区域地质特征、GPS资料及周边环境等分析发现,GPS垂直分量于同年5月开始下降,GPS垂向变化引起的重力量值-1.714μGal/mm远小于重力变化;抽水的干扰仅产生0.2181μGal的变化,远小于5μGal/月的量值变化。对GPS垂向引起的重力变化进行定量计算分析,并对水井抽水进行正演模型计算,排除抽水和垂向重力变化的影响,重力非潮汐变化量值约为14.7819μGal。研究发现,外界影响变化不足以引起重力的非潮汐异常。综合分析认为,该重力非潮汐变化异常为地震前兆异常,为后续地震预报提供了可靠依据。     

灯塔M5.1地震震前重力变化特征研究

利用辽中南地区2011~2013年流动重力观测资料,深入分析研究2013年灯塔M5.1地震震前的重力空间演化特征及变化机理,结果表明:1)辽中南地区重力场异常变化与灯塔M5.1地震在时空上具有很好的对应性;2)震前1 a重力场四象限特征明显,震前0.5 a出现50μGal的正负反向变化,并伴有重力变化高梯度带,重力点值时序上升和下降趋势明显,重力剖面呈剧烈的上升-下降波动变化;3)结合2011~2012年重力变化背景场和区域地质构造分析认为,灯塔地震可能是由辽中南地区地壳深部物质迁移和地壳构造运动共同作用引起。     

十三陵地震台gPhone重力仪的仪器性能与水文响应分析

对安置于十三陵地震台的两台gPhone重力仪(109、118号)在2013-04~08的连续观测数据进行潮汐分析和提取重力残差处理。结果表明,gPhone重力仪监测到降水导致μGal量级的重力变化;gPhone重力仪获取的二阶重力残差变化与距离十三陵地震台15km的沙河地震台静水位变化之间存在4d左右延迟的正相关对应关系,显示了地下水位观测对重力观测解释的必要性。因此,有必要对我国的连续重力观测网补充地下水位观测,提高连续重力观测的数据处理精度。     

利用GRACE卫星数据提取苏门答腊地震同震和震后重力变化

基于GRACE重力卫星数据,采用P3M6去相关滤波技术、350 km半径的高斯平滑技术与差分方法,提取2004年苏门答腊Mw9.3地震的同震和震后重力场时空演化过程。该地震断层下盘同震重力变化达到4.6 μGal,上盘同震重力变化最大值为-8.3 μGal,该结果与球体位错理论计算结果在量级和分布形态上具有较好的一致性。震后7 a断层下盘重力变化累计达到5.6 μGal,超过同震变化信号|而断层上盘震后重力变化只有0.6 μGal,仅为同震变化信号的11%。该现象表明,断层下盘所在海洋地层的粘滞性较低,地层相对柔软,其地层在震后由于较大的粘滞性松弛效应而持续变形,并伴随较大的重力变化|而断层上盘所在的地层接近陆地地壳特性,地层相对脆硬,震后受地层的粘滞性响应而导致的重力变化相对较小。     

原平地震的地面重力变化特征

利用山西地区的流动重力观测资料,分析2009年和2016年原平2次地震不同时间尺度的重力场演变特征。结果表明：1）0.5 a尺度重力场变化图像较好地反映了这2次地震震前变化较小→反向变化→在重力变化高梯度带零等值线附近发震→震后调整的演化过程;2）不同时间尺度重力场均显示,地震发生在重力变化梯度带零等值线附近,地震前重力变化量约为40～50 μGal。     

青海门源M6.4地震前重力场变化特征及其动力机制分析

利用河西地区2011~2015年相对重力资料,研究2016-01-21青海门源M6.4地震前区域重力场时空演变特征及其动力机制。结果表明,河西重力场持续多年正、负交替变化,地震前震中区域负异常持续增强并形成重力变化高梯度带,重力正负差异变化最大达到110 μGal。地震发生在重力变化零等值线附近,该等值线分布与祁连山北缘断裂、冷龙岭断裂走向基本一致。     

2012年苏门答腊MW8.6级地震震后效应模拟

基于弹性-粘弹性分层半空间中矩形位错理论和地震波反演的断层模型,结合研究区地壳-上地幔平均波速分层结构,利用PSGRN/PSCMP软件模拟计算2012年苏门答腊MW8.6地震产生的同震效应以及考虑2004年MW9.1地震效应的震后地表形变和重力变化,同时给出震后形变和重力的年变化率。结果表明,同震形变和重力变化显示发震断层左旋走滑特征,其变化主要集中于断层附近;海水质量的重新分布对同震垂直位移和重力变化影响显著;震后50 a年均变化率揭示了地表形变和重力变化的过程, 其图像明显受2004年MW9.1地震震后效应影响,呈现复杂的正负6区域分布;震后50 a近场年均形变量达10 mm/a,年均重力变化量达1.2 μGal/a,而远场年均形变量低于1 mm/a,年均重力变化量低于0.1 μGal /a;水平形变在震后400 a内变化明显,变化率逐步减小,800 a后基本稳定;垂直位移和重力变化在震后100 a内变化尤为显著,变化率快速衰减,400 a后基本稳定不变。     

重力卫星观测俯冲带特大地震同震重力场变化

利用去相关的DDK3滤波器对得克萨斯空间中心（CSR）发布的GRACE RL06月重力场模型数据进行处理,采用多项式拟合时间序列分析方法成功提取了2004年苏门答腊M_W9.3、2010年智利M_W8.8和2011年日本东北M_W9.0三个地震事件的同震重力变化、大地水准面变化、垂线偏差变化及重力梯度变化。结果表明,3个震例的同震重力变化范围分别为-15.5~6.5 μGal、-9.1~2.1 μGal和-11.1~4.2 μGal;同震大地水准面变化范围为-5.9~0.8 mm、-3.0~0.8 mm 和 -3.2~0.5 mm;垂线偏差NS向变化范围为 -1.2~2.2 mas、-0.9~1.0 mas 和 -1.1~1.4 mas,EW向变化范围为 -1.8~1.0 mas、-0.8~0.8 mas和 -0.7~1.0 mas;重力梯度各分量数值中,rr分量同震变化最大,其次是rθ分量。3次地震的同震信号空间分布均表现为：同震大地水准面和重力变化信号呈非对称两极分布;垂线偏差呈负-正-负或正-负-正三极分布;重力梯度变化信号呈较复杂的多极分布。观测结果与现有的球体位错理论计算结果较一致,可作为重力卫星数据约束震源机制和地球粘度结构的可靠信号来源。     

利用EOF分析方法研究康定MS6.3地震前后重力时空分布异常特征

对康定MS6.3地震前后的流动重力观测数据采用正交经验函数(EOF)分析方法进行处理,分析主要模式与资料时段的MS5.0以上地震震中位置的关系,探讨主要模式代表的物理含义。结果表明,川滇地区7期累积重力观测资料中大于200 μGal的数据是偶然误差的概率大于95%。第二模式规格向量和累积重力变化均具有相似的分区特征,2013~2014年6次MS6.0以上地震和第二模式规格向量正负转换带有超过83.3%的对应率,具备重力长期变化模式。第一模式地震前后规格值序列明显的同震信号、空间相似的4象限分布以及2014年4次MS6.0以上地震和第二模式规格向量正负转换带100%的对应率表明,该模式代表重力短期变化模式。康定MS6.3地震前第一模式包含4象限特点,震中位置规格值不断增大,且地震前后变化从负转为正。     

����GRACE����������ȡ���Ŵ�������ͬ�����������仯

????GRACE?????????????????P3M6?????????????350 km?????????????????????????2004?????????Mw9.3????????????????????????????????????????????????仯??4.6 μGal??????????????仯?????-8.3 μGal???y????????λ????????????????????????????н????????????7 a????????????仯????5.6 μGal?????????仯???|?????????????????仯???0.6 μGal????????仯????11%????????????????????????????????????????????????????????????????????????Ч??????????Σ???????????????仯|???????????????????????????????????????????????????????????μ??????仯????С??     

基于主成分分析的相对重力观测零漂和固体潮提取

利用主成分分析PCA提取2014-08-15~23CG-5相对重力仪静态相对观测的零点漂移和固体潮。结果显示,PCA提取的零点漂移值与最小二乘线性拟合结果非常一致,其差仅在10-2μGal/d量级;PCA估算的固体潮值与CG-5重力仪内部软件提供的固体潮模型计算值统计结果基本一致,两者之差均小于8μGal,均方根均小于5μGal。     

青藏高原大尺度地表流体的重力效应特征分析

基于2002~2017年三维大气再分析模型、水文模式及GRACE卫星重力数据,采用负荷格林函数和卫星重力反演方法分析青藏高原大尺度地表流体的重力效应特征。结果表明,青藏高原地区大气扣除季节项前后的重力效应峰对峰变化介于1.9~10.7 μGal和1.1~4.6 μGal,其负荷效应呈现出显著的季节性特征;土壤水和积雪扣除季节项前后的重力效应峰对峰变化介于0.9~9.9 μGal和0.7~8.3 μGal,其负荷效应呈现出显著的季节性和年际变化特征;GRACE时变重力场扣除季节项前后的重力变化峰对峰变化介于4.4~24.1 μGal和3.4~18.3 μGal,考虑观测误差、构造运动和信号泄漏误差的影响,其显著的季节性和年际变化反映了多种地表流体的综合影响。     

零漂改正对中国地壳运动观测网络重力数据处理的影响

基于中国地壳运动观测网络2000年流动重力观测数据,分析18台LCR-G型相对重力仪的零漂特性。结果显示,重力仪零漂率在某些时段变化较大,最大达10 μGal/h以上。根据重力仪零漂改正的不同处理方法,利用传统平差模型与改进模型对观测数据进行平差计算,并对结果进行比较。在相同先验中误差的情况下,传统平差模型的后验中误差为39 μGal,重力值平均精度为26.5 μGal,而改进模型的后验中误差为16.5 μGal,重力值平均精度为11.2 μGal,远高于传统平差模型。在数据预处理阶段按闭合观测分段进行零漂改正,可以有效消除重力仪零漂剧烈变化的影响,明显优于将仪器零漂率当作固定参数的计算方法。