青藏高原GRACE卫星重力长期变化

本文采用最新的GRACE(Gravity Recovery and Climate Experiment)(RL05)数据,通过水文模型(Global Land Data Assimilation System, GLDAS与Climate Prediction Center, CPC)扣除土壤水及雪水的影响,利用Paulson提供的冰川模型结果扣除GIA(Glacial Isostatic Adjustment)的影响,采用尺度因子的方法减少数据处理过程中误差的影响,最终基于最小二乘计算方法得到2003—2013中国及周边地区长期性重力异常变化情况.结果发现青藏高原有较为明显的重力上升信号,我们认为该信号可能由印度板块俯冲欧亚板块导致青藏高原地壳增厚所引起.接着依据GPS观测结果和艾黎均衡假说构建了地壳形变模型并通过直立长方体模型予以正演模拟分析.以班公湖—怒江断裂带为界将青藏高原划分为南北两大区块,结果显示青藏高原重力异常大致以0.2 μGal·a-1的速率在递增,小于GRACE得到的0.3±0.08 μGal·a-1的增长速率(对应于地壳增厚速率约3 mm·a-1),剩余未解释部分可能与湖水、冰川因素、冻土因素等有关.该结果对于认识青藏高原隆升动力学有一定参考意义.     

GRACE时变重力位系数误差处理滤波方法研究

由于GRACE时变重力位系数中高阶次项存在较大的误差,奇(偶)项阶数之间存在着相关性误差,直接利用GRACE时变重力场模型数据反演地表质量变化时,会使结果出现严重的条带噪声,必须采用最佳滤波函数进行空间平滑.本文从滤波理论出发,依据信噪比最大准则和模拟误差比较分析了不同滤波方法的优缺点,并以华北平原为例进行了数值分析....     

基于GRACE RL05数据研究智利地震长期震后重力变化特征

利用GRACE卫星重力可对地震引起的大范围重力变化进行观测,并从重力数据中发现主要的变化特征.发生于2010年的M_W8.8智利地震震级较高,可观测到震中附近广泛的同震和震后长期重力变化.本文基于GRACE RL05 Level-2时变重力场数据,对2010年智利地震的同震和震后长期变化进行了计算.对同震变化的计算发现,智利地震引起的同震变化极值达-5 μGal,而本文为减小水文信号的干扰而采用的3年平均的方法可以获得良好的效果.在对震后重力变化的计算中发现,智利地震震后在2011—2016年间的重力变化存在先增大后逐渐衰减的过程.对震后变化的拟合表明,智利地震震中附近有约1μGal的震后重力变化,震后变化的特征时间约1.1年.同时,在智利地震中未出现较明显的两个震后变化阶段（短期、长期）.

     

基于卫星重力的青藏高原东缘多尺度时变重力场研究

利用经过去相关滤波处理的GRACE时变重力场模型获得了青藏高原东缘2003—2012年的卫星年重力变化图像,并针对该区域近年发生的三次特大地震,结合震前及震后月重力场变化图像,分析与强震有关的卫星重力场变化特征。从区域年重力变化图像可以看出,三次大震均发生在年重力变化较低的时段内,震前小幅值变化可能是地震发生的中短期前兆;从汶川地震和玉树地震发生前后的月重力场变化图像可以发现,发震前后断层附近的重力变化模式发生变化,这可能印证了震后位场变化恢复理论;从汶川地震前后的龙门山断层附近点上的周重力变化趋势可以明显发现,汶川地震发生（第20周）后近9周的时间,断层东西侧呈现了相反的重力变化特征,这可能是对震后壳幔物质调整过程的反映。     

航空重力测量数据的小波滤波处理

构造三类连续小波函数对航空重力测量数据进行小波滤波处理. 三类连续小波函数分别用于对测量数据在某一空间尺度（或时间尺度）上的低通,一阶求导和二阶求导滤波. 着重介绍三类连续小波函数的构造原理与过程,并说明其相对于传统的数字滤波器的优势. 对系统的技术参数（滤波器窗口宽度参数δ 和尺度参数s）进行了调试实验. 实算结果显示了方法的可行性和有效性.     

GRACE和地面重力测量监测到的中国大陆长期重力变化

自2002年以来,GRACE卫星探测计划可提供高精度的时变地球重力场,用以探测地球系统的物质分布.自1998年中国大陆重力监测网建立以来,利用FG5绝对重力仪和LCR-G型相对重力仪每2年对该网进行重复测量获取重力场时变信息.基于此,本文利用GRACE和地面重力测量获得了中国大陆重力场的长期年变率,利用位错理论根据USGS发布的断层模型计算了2008年汶川Ms8.0级地震的同震重力变化并进行了300 km高斯滤波.GRACE卫星重力和地面重力结果均表明华北地区地下水流失严重,在绝对重力基准站上,GRACE卫星重力与绝对重力变化率较为一致,汶川区域的地面重力变化结果可视为大地震前兆信息.     

基于GRACE卫星和地表流动重力的汶川MS 8.0地震前重力时空变化

利用GRACE卫星重力场模型和地表流动重力观测资料,计算2008年汶川M_S 8.0地震发生前6年的重力变化,对卫星和流动重力段差结果与卫星重力反映的重力场动态变化特征进行研究,结果表明：①GRACE卫星重力段差受滤波半径影响显著,与地表流动重力观测结果相比,在重力变化数值上差异较大,在变化率上较为一致;②在汶川地震孕育阶段,川滇地区重力等值线呈“增大—减速增大—减小”的特征,震前2年形成近似垂直于龙门山断裂带的重力变化梯度带。     

重力卫星GRACE检测出2010年智利Mw8.8地震的同震重力变化

本研究通过对重力卫星GRACE观测数据的处理,采用去相关加300 km半径的高斯平滑,成功地提取了2010智利M_w8.8地震所产生的重力变化信号,最大变化幅度达到7 μGal,并且与位错理论计算结果具有较好的一致性.这是继GRACE检测出2004苏门答腊M9.3地震重力变化后的又一个卫星观测地震的例证,说明GRACE具有检测出M<9.0量级地震的能力,为利用GRACE研究地震以及其更广泛的应用提供了可靠的依据.     

基于GRACE卫星RL05数据的南极冰盖质量变化分析

CSR(Centre for Space Research)最近发布了RL05数据,其空间分辨率、精度和周期变化特性等都优于RL04数据.本文采用300 km的扇形滤波及P5M11去相关滤波削弱南北条带等重力场模型误差,并采用Paulson2007模型进行冰川均衡模型改正,利用CSR RL05与RL04数据计算分析了南极2002年到2012年的质量变化序列及其变化趋势的空间分布特性,并选取8个特征点进一步分析了其质量变化序列.同时,对CSR、JPL(Jet Propulsion Laboratory)、GFZ (GeoForschungsZentrum)三个机构发布的RL05数据采用相同的滤波方法进行计算,得到整个南极的质量变化分别为-195.7±20.5 Gt/a、-203.8±23.1 Gt/a、-133.2±29.9 Gt/a,对全球海平面变化的影响分别为0.54±0.06 mm/a、0.56±0.06 mm/a、0.37±0.09 mm/a.     

利用GRACE观测数据研究苏门答腊区域的黏滞性结构

本文尝试利用卫星重力观测资料和震后黏弹性松弛理论研究苏门答腊地区的区域流变结构,为更好地认识区域地球动力学环境提供依据．利用GRACE卫星重力资料,计算了2004年苏门答腊Mw9．3地震的同震及震后的重力变化．计算中使用平滑半径为500km的高斯滤波器．结果显示苏门答腊Mw9．3地震破裂东侧以陆地为主的上盘同震下沉量很大,造成约9x10＾-8mS＾－2的重力下降阶变,而西侧处于海水下的下盘重力同震上升较小约2x10＾－8mS＾－2,但其震后上升较快．流变结构对岩石的变形有很大的影响,是地球动力学数值模拟取得可靠科学结果的基础．本文尝试了利用卫星重力变化时间序列来反演苏门答腊地区的黏滞性结构;即基于GRACE时变重力场,利用自重力、黏弹性、平面分层模型模拟了该地震的同震和震后变形,并将获得的空间固定点的重力变化与GRACE重力场及点位时间序列进行比较,估计该地区的黏滞性系数在1．0x10＾18PaS的量级,且断层两侧的流变参数存在差异．最后结合苏门答腊区域的构造特点讨论了黏滞性系数的影响因素．     

起伏地形条件下流动重力数据处理方法研究

以地震监测为目的流动重力重复测量获取的重力年变信号通常在十几到几十微伽量级.但由于重力测点所处的地形高程不同,将空间测点直接插值得到的重力异常空间图像会叠加由于地形起伏引起的异常畸变.本文通过6面体模型定量计算了不同规模和埋深场源的地形效应;对青藏高原东缘地区测网的重力观测数据进行了分析处理,并通过位场曲面延拓校正技术消除了地形影响.结果表明,位场曲面延拓校正技术可以有效地改善重力数据质量,有利于更准确地圈定和认识异常的空间形态特征.     

航空重力数据样条平滑滤波试验研究

采用B样条平滑滤波技术,对GT-1A航空重力勘查系统的测量数据进行了滤波试验研究,结果表明,通过选择合适的样条、平滑等参数,B样条平滑滤波技术在航空重力数据的低通滤波数据处理中可以发挥重要作用,获得与GT-1A系统滤波结果几乎同样满意的滤波效果,二者滤波结果差值的标准偏差可达±0.4 mGal.     

航空重力测量数据的滤波与处理

我国首套航空重力测量系统(CHAGS)于2002年研制成功.本文设计了CHAGS观测数据的滤波处理流程,着重讨论了比力及各项改正的滤波计算方法,提出了'两步滤波法'.基于某航空重力测量实测数据,对测线重力异常的计算精度作了比较和分析,重复测线计算结果表明,CHAGS确定测线重力异常的内符合精度优于±3mGal,相应的波长分辨率为8km;与地面实测数据向上延拓数值比较表明,对于同样的波长分辨率,外部符合精度优于±4mGal.     

重复绝对重力测量观测的滇西地区和拉萨点的重力变化及其意义

根据近年来在滇西地区和西藏拉萨绝对重力点上的一些新的绝对重力重复观测结果 ,对1990年以来在这些点上绝对重力重复观测结果进行了分析. 滇西地区的5个绝对重力观测 点的近10年的重复观测结果表明,大部分观测时段没有出现与地球动力学事件有关的重力变 化,只有丽江和洱源2个绝对重力观测点的观测结果显示丽江地震前后有变化. 为了研究该 重力变化的原因,本文正演模拟了出现重力变化期间丽江地震同震位错引起的重力变化,模 拟结果与实际观测结果有较好的一致性. 西藏拉萨近10年以来的重复观测结果给出了拉萨点 的绝对重力值以-1.82±0.9×10-8m·s-2·a-1速率下降,这从重力学 的角度反映出青藏高原的隆升. 对拉萨点的重力变化机制进行了探讨,根据印度板块与欧亚 大陆俯冲模型计算的拉萨点重力变化速率与观测到的重力变化速率较一致,表明现今拉萨地 区重力变化是由于印度板块与欧亚大陆俯冲所引起. 根据印度板块向欧亚大陆俯冲的位错模 型计算的拉萨点的重力变化与地表垂直位移的关系,将重力变化转换为拉萨点的隆升速率为 8.7mm/a.      

利用GRACE卫星数据研究汶川地震前后重力场的变化

介绍了利用GRACE卫星数据计算和分析重力场变化的方法.以EIGEN-GRACE02S重力场模型为背景,对高斯平滑处理前后中国大陆月重力场的变化作了比较,认为取平均半径为666 km作高斯平滑,能够得到合理的重力变化结果.以3个月时间尺度计算了2008年汶川地震前后中国大陆的重力场变化,将结果与2006年和2007年相同时间尺度的变化结果进行比较,发现三者相符合,并且同震重力的负变化可用地壳膨胀模型理论来解释.     

起伏地形条件下流动重力数据处理方法研究

当前在我国大陆基于固定测点的流动重力重复测量是获取区域性重力场变化信号的主要观测手段.通过定期的重复测量,可以了解重力异常变化的空间特征,这些变化信息与区域构造变形和多种动力学过程密切相关.     

利用GRACE数据反演日本Mw9.0地震区域黏滞性

利用CSR(Center for Space Research)发布的GRACE RL05月重力场模型数据,通过水文模型GLDAS(Global Land Data Assimilation System)和CPC(Climate Prediction Center)扣除土壤水和雪水的影响,根据冰川模型扣除GIA(Global Isostatic Adjustment)的影响,采用P3M6去相关滤波和300 km扇形滤波,基于最小二乘拟合的方法得到日本M_W 9.0地震的同震及断层上下盘两个特征点重力变化时间序列,利用PSGRN/PSCMP模型对日本M_W 9.0地震区域黏滞性进行了反演,并计算了同震及震后5年研究区域重力变化的空间分布.结果表明,扣除土壤水和冰川均衡调整因素的影响,同震重力变化为-5.2×10~(-8)～2.9×10~(-8) ms~(-2);两个特征点在震后5年重力均增加,下盘重力增加较大;日本M_W 9.0地震区域黏滞性横向差异较明显,断层上下盘的地幔黏滞系数分别为2.5×10~(18) Pa·s、5.0×10~(17) Pa·s时,与GRACE观测值较接近,综合考虑断层上下盘的震后重力变化,区域黏滞系数大约为1.5×10~(18) Pa·s.     

gPhone重力仪与GRACE观测的陆地水季节性重力变化对比

利用中国大陆构造环境监测网络中张家口、淮北、溧阳重力台的gPhone连续重力观测数据,经过Tsoft软件数据预处理、固体潮、海潮负荷、极移等主要潮汐信号的扣除、局部气压改正,以及分段拟合仪器漂移后,得到观测数据的非潮汐信号,该信号时间序列呈现明显的季节性,主要反映陆地水储量变化;在此基础上,将地面重力非潮汐信号的月均值与重力卫星GRACE的观测结果进行比较。结果发现,仪器观测较稳定的张家口站和淮北站,地面重力和卫星重力观测的季节性信号具有一定的一致性;而溧阳站两种观测的相关性较差,可能与仪器观测的稳定性及台站周边区域环境因素有关。     

重力卫星GRACE Mascon产品的应用研究进展与展望

2002年重力卫星GRACE的成功发射极大地促进了地球科学多个领域,包括全球海平面变化、极地冰盖与高山冰川消融、水文以及固体地球等多个领域的发展.然而,GRACE观测数据主要是以球谐系数的形式给出,需要应用者进行一系列预处理才可以得到对应的物理量.为了克服此困难,也为了提高GRACE恢复重力场地空间分辨率,相关机构在近些年推出了新一代GRACE观测数据产品,即Mascon产品.该产品的初衷是便于非大地测量和地球物理专业的人使用,比如水文学家、海洋学家,它无需进行任何后处理过程,使用上更加方便.然而,尽管Mascon产品以较高的空间分辨率(如1°)给出,但是,该产品的应用范围以及其实际的分辨率等都是科学家们非常关注的问题.目前科学家们已经对该产品在不同流域尺度以及不同应用领域上的适用性问题进行了系统性地评估.本文综合介绍了Mascon产品的基本原理和方法、三家Mascon产品的差异,并梳理了该产品和球谐系数产品之间在一些具体物理问题的应用中的适用性以及应该注意的问题,为广大科研工作者提供科学依据和使用参考.     

GRACE卫星观测到的与汶川Ms8.0地震有关的重力变化

利用GRACE卫星重力资料,计算了中国大陆及周边的卫星重力时变场和地表密度变化分布,获取了具有代表性的点位区域的每月重力变化时间序列.同时获得了WUSH、LHAS、KUNM、LUZH站相对于区域参考框架的GPS位移时间序列.卫星重力观测结果显示喜马拉雅弧形带的重力在2004年苏门答腊M_w9.3地震后快速下降, 2006~2008年尤为明显,西域地块西北边界带上震后重力下降也较为显著;而沿青藏高原北至东边界2007年出现明显的重力上升沿构造边界的弧形分布,且2008年南北地震带中南段重力上升变化显著.这些苏门答腊地震后的重力变化趋势到汶川地震发生后才开始改变.GPS位移结果显示四个台站均记录到苏门答腊大地震的同震信号,震后WUSH、LHAS、KUNM站水平位移向量出现明显的运动趋势改变,且一直持续到2008年汶川M_s8.0地震的发生.GRACE卫星揭示的青藏高原及周边地表质量的变化为解释汶川地震的动力机制提供了新的观测途径和资料.本文结合区域构造运动的特点和GPS位移,对GRACE观测的时变重力场特征及汶川地震的动力机制进行了初步解释和讨论.