2016年青海门源MS6.4地震前重力变化

利用青藏高原东北缘2011-2015年期间的流动重力观测资料,系统分析了区域重力场变化及其与2016年1月21日青海门源M_S6.4地震发生的关系,结合GNSS、水准观测成果和区域地质构造动力环境,进一步研究了区域重力场变化的时空分布特征及其机理.结果表明：（1）测区内重力场异常变化与祁连山断裂带在空间上关系密切,反映沿祁连山断裂带（段）在2011-2015年期间发生了引起地表重力变化效应的构造活动或变形.（2）门源M_S6.4地震前,测区内先出现了较大空间范围的区域性重力异常,到临近发震前显示出相对闭锁的现象,且围绕震中区周围出现四象限分布特征的局部重力变化,地震发生在重力反向变化过程中,并出现显著的四象限分布特征的重力异常变化,其中,青海门源与甘肃天祝一带重力差异变化达100×10^-8m·s^-2以上.（3）区域重力场动态演化大体反映了青藏高原东北缘物质北东流的动态效应,门源震中附近区域地壳受挤压变形显著、面压缩率和重力剧烈变化的特征最为显著.（4）重力场的空间分布及其随时间变化与地壳垂直与水平运动及地质构造活动等观测结果有一定的对应关系,强震易发生在重力变化四象限分布中心地带或正、负异常区过渡的高梯度带上.     

三江地区不同场源重力场特征及其与构造、地震的关系

利用重力场实测资料及最新的高分辨率超高阶重力位模型IGG97LB（720阶）计算了三江地区不同场源的重力场。在分析研究了不同场源重力场特征之后，结合地震层析成像资料，得出了初步结论：（1）三江地区的康滇菱形块体内部的物质是由青藏高原下面物质向南东方向延伸来的，可见康滇菱形块体与青藏高原的关系密切；（2）康滇菱形块体的低速度层厚度大于周边地区；（3）从三江地区重力场特征及历史上地震发生的情况看，地震容易发生在重力异常等值线拐弯处，特别是复杂的等值线拐弯处，更是地震的多发区。     

重力波包在中层大气温度波导中传播的数值模拟

给出了重力波在中层大气温度波导中的反对称形式导制传播的线性理论模型，并采用二维非线性的数值模型对重力波波包在中层大气温度波导中的传播和演变过程进行了模拟研究.模拟的结果表明，下层大气激发的重力波能量进入波导区域后被俘获形成导制传播.重力波在波导内不停地来回反射，垂直方向的自由传播受到限制，能量在波导内沿着水平方向传输，模拟得到的波参数与理论值相近.重力波包在温度波导中传播时伴随着能量泄漏，且能量泄漏的速率随时间变缓，最终总有部分能量被限制在波导区域.重力波在水平方向上传播几百公里后，依然维持着良好的波结构，同时数值模拟也给出了重力波在波导区域内能量密度的时空分布.     

应用GPS连续观测和流动重力资料分析新疆洛浦MS6.0地震的前兆特征

对2012年3月9日新疆洛浦MS6.0地震前陆态网络GPS连续跟踪站数据与喀什—伽师2010年7月～2012年3月期间重力复测数据进行计算分析。结果表明：（1）断层两端运动速率不同,为应力积累创造条件;（2）地震发生前监测区出现区域性的重力异常变化高值区,区域重力场等值线与发震断裂走向平行;震前区域重力场出现增强—减弱—增强过程。     

中国几处城市高楼重力垂直梯度测定

重力垂直梯度测定是应用地球物理学和物理大地测量学的一个重要的实验性课题。研究重力垂直梯度（以下简称垂直梯度）在地球表面的分布与变化,可以增加对地球重力场和地下地质构造及其活动性的了解。1982-1984年,我们在研究垂直梯度测量原理和     

云南地区及邻域卫星重力变化及其与区域大震的关系

利用GRACE卫星重力场数据,计算并分析云南地区及邻域2003—2014年的卫星重力场变化特征,获得3个特征点上的月重力变化时间序列;结合该地区2003—2014年发生的M_S5以上地震分布,探讨区域卫星重力变化与区域大震的关系。结果表明:云南地区及邻域卫星重力场变化对区域M_S5以上地震有一定的反映,地震多发生在区域重力变化的高梯度带上,地震发生前后,区域出现与地震相对应的特征性卫星重力变化,可能是对壳幔物质调整的反映;区域大震前后,不同震区特征点的重力变化特征存在差异,说明分断层、分区域的重力场变化研究对地震预测有重大意义。     

基于变差函数的插值方法计算华北地区重力场变化

定点重复重力测量是获取区域重力场变化的主要手段之一, 重力场特征与地形起伏、 构造走向等因素相关。 以华北地区为例, 考虑区域内地形和构造的北东向分布规律, 从EGM2008重力模型中拟合各向异性变差函数参数, 利用变差函数网格化插值, 对华北地区2009—2013年期间重力场观测数据进行网格化重建, 获取华北地区重力场时空变化结果。 研究结果表明基于区域地形、 构造特征的各向异性变差函数插值方法, 获得的空间重力场变化在重力异常梯级带上更加明显, 重力变化与活动构造分布具有更好的一致性。 本文研究方法对于恢复区域时空重力场异常具有重要意义, 有助于提高应用重力资料划分潜在地震危险区空间位置的精度, 为华北地区震情研判和构建地震预报定量指标体系提供可靠的地球物理场数据。     

中国大陆西部重力场变化与强震危险性关系

本文基于中国大陆西部流动重力测量获得的区域重力场变化数据, 尝试对其进行多种位场变换, 在一定空间尺度约束下定义了6种区域重力场变化指标量. 在此基础上, 以中国大陆2002年、 2005年和2008年等3期数据为研究对象, 计算并统计了中国大陆西部M_S≥6.0地震震中位置处震前重力场变化及各种导出参数值, 尝试初步构建以地震预报为目的的地震重力学科指标体系参数. 研究结果表明, 震前区域重力场变化与地震发生位置没有显著的统计关系, 但是重力场变化的垂直梯度和解析信号模量具有一定的地震预报意义, 适合作为地震前兆异常使用.      

唐山地震(M=7.8)前后的区域重力变化和扩容过程

本文在考虑各种干扰因素(如温度、气压、地下水和台站高程变化等)的同时,详细分析了唐山地震的后区域重力场变化的几个特征。分析结果表明,唐山地震前发生了明显的区域重力场异常。异常有四个特点:(1)重力变化围绕着震中,最大变化值恰好在震中;(2)重力随时间的变化是非线性的;(3)变化最大值出现在1975年年中,而不是在1976年7月份;(4)1977年区域重力场基本恢复到震前水平。不能用台站高程变化和地下水位变化来解释这些异常。作者用扩容模式计算了区域重力变化理论值,计算结果同观测值相符。此外,作者用相应地区的垂直形变资料和地下水观测资料讨论了该模式的合理性。     

应用GPS连续观测和流动重力资料分析新疆洛浦MS6.0地震的前兆特征

对2012年3月9日新疆洛浦Ms6.0地震前陆态网络GPS连续跟踪站数据与喀什—伽师2010年7月～2012年3月期间重力复测数据进行计算分析.结果表明：（1）断层两端运动速率不同,为应力积累创造了条件;（2）地震发生前监测区出现区域性的重力异常变化高值区,区域重力场等值线与发震断裂走向平行;震前区域重力场出现增强—减弱—增强过程.     

龙羊峡震群前后的地形变特征

以龙羊峡附近的流动重力监测资料和精密水准复测资料为基础，分析了1996年12月－1997年3月龙羊峡震群前后该地区的重力场和地壳垂直变化特征。分析认为：震前重力场变化异常，临震前形成重力变化密集带，地震发生在重力变化密集带上，震后重力变化密集带消失，重力场变化平衡，垂直形变平缓。     

2020年伽师MS6.4地震重力与地壳形变综合分析

2020年1月19日新疆维吾尔自治区伽师县发生M_S6.4地震,此次地震为N倾逆冲兼少量走滑型地震,地震的发震构造是南天山柯坪塔格逆冲推覆体最南缘的柯坪断裂,伽师地震的震中位于喀什-伽师流动重力观测网和中国大陆构造环境监测网络GNSS观测网监测区域内。文中选取2015—2021年的流动重力和GNSS观测资料开展分析。采用绝对重力点约束流动重力数据,并进行全网经典平差计算,进而获得重力场变化图像;基于GNSS数据通过GAMIT/GLOBK软件解算获得研究区的运动速率,借助位移与应变的偏导关系获取水平视应变场分布图。之后,分别研究了孕震区的重力场和应变场的变化特征,并综合探讨了孕震区的重力场、 GNSS形变场特征与伽师M_S6.4地震孕震过程的关系。研究表明：1)一年时间尺度的重力场时变图像较好地反映了伽师M_S6.4地震附近区域重力场的系统演化过程,重力变化零等值线拐弯处及“四象限”特征可为地震预报提供有益参考。2)重力场累积变化反映出伽师M_S6.4地震前的重力场变化受到区域大断裂的控制,研究区一系列中强震...     

重力场动态变化与汶川MS8.0地震孕育过程

基于中国大陆1998～2007年（复测周期2～3年）流动重力观测数据,结合GPS、水准观测成果和区域地质构造动力环境,分析研究了汶川8.0级地震区域重力场动态变化演化特征和孕震机理.结果表明：区域重力场动态演化大体反映了青藏高原物质东流的动态效应和汶川大震孕育的中长期（2～10年）信息;汶川大震孕育的显著重力标志为震中西南持续多年的正重力变化（上升）和出现较大规模的重力变化梯级带,前者有利于地震能量的不断积累,后者有利于地震剪切破裂的发生;与地震孕育相关重力场变化总体呈增大—加速增大—减速增大—发震的过程;8年累积重力变化幅差最大约200×10^-8m·s^-2;2001年昆仑山口8.1级地震孕育发生和震后恢复调整,对区域重力场动态变化和汶川大震的孕育发展具有重要影响;松潘—甘孜块体一般呈现负重力变化,可能反映深部壳幔局部上隆、壳内温度较高而膨胀,有利于逆冲或推覆体运动的形成和大震的发生.     

用最小二乘配置法构建局部重力场模型

基于非均匀分布的陆地重力观测数据,重构局部重力场模型是区域重力资料处理与解释的重要环节。本文对比了多种局部重力场建模方法,并以EGM2008模型提供的自由空气重力异常模型重采样数据进行测试,综合比较了不同噪声条件下不同建模方法的实际效果。结果表明:在不同噪声水平下,优选出适合重力位场问题的协方差函数后,最小二乘配置法的建模效果优于其它方法。      

重力场动态变化与汶川M_S8.0大震

地震重力观测主要监测地表或空间重力场时间变化,即动态重力变化。动态重力场变化是地表位移运动（主要为测点垂直位移）和地球内部物质运移或变化的综合反映,是地球内部运动（主要为地壳运动,其表现为物质密度变化和变形的耦合运动）的重要表征物理量。地震的孕育过程包含于地壳运动过程中,因此,通过对动态重力场变化的深入分析,     

彰武地震前后垂直形变场局部反向原因探索

以探索彰武地震前后垂直形变场局部反向原因为目的，利用动态秩亏平差方法处理辽宁地区水准、重力复测资料。依据重力场变化图按网格法计算重力场变化对垂直形变场的影响，绘制消除重力场变化影响的垂直形变场图。     

日本MW9.0地震震区及其周缘2002—2015年卫星重力变化时空特征

本文在考虑全球陆面数据同化系统陆地水储量变化影响后,利用2002年4月—2015年4月GRACE卫星RL05月重力场模型获取了2011年日本M_W9.0地震震中及其周边区域的重力场信息;然后给出了日本及其周边2003—2015年的年度累积重力变化和差分重力变化,并且利用经验正态函数方法深入分析了该地震过程中同震重力变化对区域重力场的贡献. 结果显示:日本M_W9.0地震前研究区域出现了幅值高达2×10-8 m/s2的异常重力变化,其同震效应的时间和空间特征均出现在第一模态,且同震重力变化和震后区域重力场变化特征显著,这充分表明该地震对区域重力场的影响显著.      

青藏高原东北缘近期重力与地壳形变综合分析与研究

利用青藏高原东北缘2011—2014年期间的流动重力与GPS观测资料和1970—2011年水准观测获得的垂直运动背景场资料,结合区域地质构造动力环境和区内发生的强震事件,系统分析了区域重力场变化、地壳三维形变与区域构造环境和强震活动的关系,进一步研究和探讨了区域重力场变化与区域三维地壳运动的时空分布特征及其机理。结果表明:1)测区内重力场变化与GPS观测反映的水平运动和水准观测反映的垂直运动在空间上关系密切,地壳形变总体表现为沿水平运动的方向重力增加,高原山地挤压隆升,地面重力减少;断陷盆地相对伸展下沉,地面重力增加,反映了新构造活动的继承性。2)地壳形变场变化与活动断裂密切相关,重力变化和垂直形变等值线走向总体上与呈NWW向的祁连-海原断裂走向基本一致,水平形变也在祁连-海原断裂附近产生明显的左旋走滑运动。3)强震易发生在具有显著重力变化的活动断裂带上,也易发生在沿活动性断裂的断块垂直差异运动强烈或兼有强走滑运动的地方。2016年门源震中附近区域地壳受挤压变形显著、面压缩率和重力剧烈变化的特征最为显著。2013年岷县漳县6.6级地震发生在重力变化高梯度带拐弯的地区,也是面压缩率变化过渡带和垂直形变过渡带地区。4)文中第一作者等曾基于该区地壳形变资料在岷县、门源地震前做过一定程度的中期预测,尤其是地点预测。基于上述认识,进而强调在研究区形成的一些地壳形变异常部位可能仍存在中长期强震/大震危险背景。     

GOCE卫星监测的中国区域重力异常

介绍了利用最新一代重力卫星GOCE(Gravity field and steady-state Ocean Circulation Explorer)重力场模型数据计算重力异常的原理和方法,采用最新发布的GOCE重力场模型数据(2009年10月-2010年7月)计算了中国区域(70°～130°E,15°～55°N)的重力异常,为了进一步分析重力异常与区域地质构造及地震活动性的对应关系,将计算结果与中国区域地形、地震活动区域等资料进行了对比分析,结果表明利用GOCE重力场模型数据计算的重力异常能够较好的反映区域地质构造分布特征,而且强震震中通常位于重力异常变化剧烈的高梯度带上.     

中国大陆高精度绝对重力观测数据分析和异常计算

高精度重力异常模型是研究地球内部结构的基础地球物理数据.本文依据中国大陆构造环境监测网络项目的 100个高精度绝对重力基准站数据,首先,分别与EGM2008模型和GGMplus模型进行对比,得到了模型之间的差异性特征统计参数;其次,我们对与绝对重力同点位测量的垂直梯度数据,进行统计分析,结果表明重力垂直梯度变化的标准差高达537.2E.通过提取Crust1.0地壳结构模型参数,计算得到了Moho界面起伏引起的垂直梯度变化,分析认为地表垂直梯度的变化主要与地形和浅部密度异常体相关.最后,基于重力场模型和观测数据的统计结果,针对经典重力异常计算公式中的一些流程提出了优化建议.本文结果对于使用区域重力异常模型,研究地壳结构和提高长剖面、小比例尺重力探查项目中的异常计算精度具有重要意义.