GRACE RL06与RL05时变重力场模型数据初步比较分析

简要介绍CSR、GFZ和JPL机构的GRACE RL06时变重力场模型数据,并对比分析RL06和RL05数据的解算模型。从全球陆地水储量变化反演结果、时变重力场模型阶方差和C20项时间序列3个方面,对2004-01~2014-11期间RL06和RL05时变重力场模型数据进行对比分析。结果表明,GRACE RL06时变重力场模型数据质量的和精度较RL05有明显提高,其全球陆地水储量变化反演结果去条带噪声效果更好、信噪比更高;在高阶项部分,RL06模型数据的阶方差小于RL05;RL06模型数据的C20项时间序列幅值变化小于RL05,与SLR所得C20项数据也更接近。相同条件下,采用CSR RL06模型阶方差最小,利用RL06模型所得全球陆地水储量变化反演结果信噪比值最大。     

OpenMP并行算法在卫星重力场模型反演中的应用

利用卫星重力数据反演地球重力场需要解决重力场模型的高效计算问题。分析了最小二乘直接法求解重力场模型涉及的密集型计算任务,基于OpenMP实现了卫星重力场模型直接求解的并行算法。利用30天、5秒采样间隔的沿轨扰动位T和径向扰动重力梯度Tr数据,分别反演了60阶次的卫星重力场模型,计算结果表明,OpenMP并行算法能够有效提高直接法求解卫星重力场模型的计算效率,并具有很好的稳定性。     

基于能量守恒方法的重力场反演快速异构并行算法

在利用卫星跟踪卫星资料解算重力场模型位系数时,其海量观测数据处理以及大型方程组解算过程存在计算效率低下、对平台硬件要求高等问题。针对以上问题,提出一种基于能量守恒方法的重力场反演快速异构并行算法,基于CUDA 在GPU端实现并行计算设计矩阵,结合MKL库与分区平差法、预处理共轭梯度法在CPU端完成低内存消耗下的法方程快速构建与求解,实现重力场模型反演的异构并行计算。运用该算法处理GRACE-FO卫星2020-01-01～06-30期间观测数据,反演获得120阶重力场模型GM-GraceFO2020h;与现有模型以及算法对比分析表明,该算法所得模型与现有GRACE重力场模型精度相当,且相较于传统的串行算法,反演耗时减少98.479%,内存消耗减小1个数量级。     

区域地壳结构差异对GRACE-FO反演中国近10个月陆地水负荷垂直形变的影响

利用CSR、JPL、GFZ三家机构最新发布的新一代重力卫星GRACE-FO（GRACE-Follow On）时变重力场模型揭示2018-06~2019-05我国陆地水储量随时间演变的时空特征,研究利用Hard与PREM模型分析区域地球结构差异对GRACE-FO重力数据反演我国近10个月陆地水负荷垂直形变的影响。结果表明,60阶径向负荷勒夫数的相对差异为-4.267%,前60阶垂直负荷形变最大差异位于云南一带的澜沧江流域地区,幅值约为0.7 mm/a。因此在使用GRACE-FO重力卫星数据估算中国大陆地表垂直负荷形变时,区域地球结构的影响不容忽视。     

SWARM卫星时变重力场反演及其精度

基于SWARM卫星的精密轨道数据,利用短弧积分法解算2015-01～2021-12共84个月的40阶次TVG-SWARM月时变重力场模型,并与ASU、COST-G、IGG和ITSG等机构的月时变重力场模型进行比较。结果表明：1)从大地水准面阶误差与模型位系数误差谱看,不同SWARM模型的低阶位系数精度相当,特别是前10阶均与ITSG-GRACE/GRACE-FO接近;2)不同SWARM模型与ITSG-GRACE/GRACE-FO模型全球陆地水储量变化趋势空间分布具有较好的一致性,在亚马孙流域、格陵兰岛、密西西比河流域和西西伯利亚等区域,TVG-SWARM与ITSG-GRACE/GRACE-FO模型的趋势差值分别为0.23 cm/a、0.27 cm/a、0.57 cm/a和0.47 cm/a,相关系数均达到0.85以上,并与IGG-SWARM模型结果最为接近。本文研究结果证明了TVG-SWARM模型精度可靠,可以用于监测大尺度陆地水储量变化。     

中国大陆及邻GRACE卫星重力变化研究

探讨如何运用GRACE数据进行重力变化的计算和分析.对GRACE重力场位系数的最大阶次选取进行了分析,并利用月卫星重力场模型计算了中国大陆及邻区的重力变化,通过平均半径为800km的高斯平滑,得出了合理的重力变化结果.     

GRACE月重力场模型的地震应用

介绍GRACE卫星的发展状况;根据强地震中长期预报需求，对卫星重力场模型的阶次截取并利用截取后的重力场模型进行重力变化试算，得出了一些有意义的结果。     

星载GPS实时定轨中重力场模型阶次的最优确定方法

分析重力场模型的截断误差及其影响因素,并用多颗不同轨道高度的低轨卫星GPS实测数据进行实时定轨仿真分析,提出实时定轨精度能够达到1.0 m(伪距法实时定轨)和0.5 m(相位法实时定轨)的重力场模型阶次的最优确定方法。使用CHAMP和TerraSAR-X 2颗卫星实测数据,验证了重力场模型阶次的最优确定方法的可行性和有效性。该方法不仅避免了传统经验法确定重力场模型阶次不准确的问题,而且不需要借助精密轨道进行精度评估,便于实时定轨系统的工程应用。     

利用新版GRACE时变模型反演珠江流域陆地水储量变化

本文以CSR05、GFZ05、JPL05、ITSG-Grace2016、CSR06、GFZ06、JPL06及ITSG-Grace2018等8种GRACE时变重力场模型作为研究对象,通过300 km高斯滤波和Swenson方法等数据处理策略去除时变模型的噪声,反演2005-01~2012-12珠江流域的陆地水储量变化,并利用ITSG-Grace2018模型、GLDAS水文模型、降水数据、实测地下水数据等资料对珠江流域的陆地水储量变化进行综合分析。结果表明：1）3大官方机构发布的2006版GRACE时变模型的解算精度较2005版有显著改善,ITSG-Grace2018模型的精度也较ITSG-Grace2016有一定的改善,且新版在不同时变模型间的差异更小;2）ITSG-Grace2018模型反演的珠江流域陆地水储量整体呈上升趋势,其季节性变化特征与GLDAS水文模型、降水数据及地下水测井监测数据具有较好的一致性。     

利用GOCE卫星观测数据反演地球重力场模型

利用GOCE卫星约6个月的重力梯度数据和约1 a的几何轨道数据,联合解算250阶次的地球重力场模型TJGOCE01。GOCE重力梯度数据的低频误差采用ⅡR数字滤波器处理,粗差采用阀值法和移动窗口阀值法组合探测与剔除。直接在梯度仪坐标系中建立GOCE卫星的重力梯度观测方程,采用改进的短弧边值法建立几何轨道观测方程。两类观测值的权根据其先验精度确定,采用Kaula规则约束的正则化方法解算法方程。解算的TJGOCE01模型相对于EIGEN6C2模型在250阶次的大地水准面误差和大地水准面累积误差分别为19.4 mm和177.9 mm。北美地区GPS水准观测数据的检验结果表明,TJGOCE01模型的中误差为0.544 m,略优于欧空局公布的同阶次的第二代时域法和空域法解算的GOCE重力场模型。     

�۲���������һ���������������Ӱ�����

????????????????????????й???????????????????????????????о?????????????????????仯????????????????????????????????????????????????????????????????????????и???????????????2??3 cm?????￡??????????????????????????????????侫?????????????????????????????????????????nm??????????????????????????????????????????????????????????????????????????????????????????У??????????????????????     

EGMTools：重力场模型评估及重力量计算软件

高精度、高分辨率的地球重力场可用于统一全球高程基准及确定大地水准面等,因此全面评估重力场模型精度及其快速导出的各重力量（如高程异常等）具有重要意义。本文利用VS2015+Intel XE 2016平台及OpenMP等并行技术研发一套用于重力场模型评估及相关重力量计算的软件--EGMTools,该软件的主要功能有：1）可实现基于实测点或格网重力异常、高程异常及垂线偏差等数据对重力场模型截断至任意阶次的精度进行评估;2）可实现利用重力场模型计算任意点或格网（SRTM地形表面）重力异常、高程异常及垂线偏差等重力量;3）在全球范围内,可按照一定的格网间隔（如1°×1°）对不同重力场模型所表示的重力量（如高程异常等）按不同的阶次进行比较;4）重力场模型格式转换、阶方差分析及RTM重力量计算的数据预处理等。     

���������CHAMP�������ָ��е�Ӧ��

??????????CHAMP?????????????????????????????н?????1???????????С?????LS??????????CHAMP?????????????????????????????????????????????????????????????????????????????????????CHAMP??????????У???????????????????????????1?GFZ????????????1???? CHAMP????RSO????????????????????????20???????(5??100) m ²s ^-2??С??????????????LS????????????60??60???????????????????EIGEN??CG03C??EIGEN2???????б?????????????????????Ч?????????????????CHAMP????????????????????LS???????????????ε?EGM96??EIGEN1S??EIGEN2????     

基于EGM2008和剩余地形模型的区域似大地水准面精化方法

利用EGM2008计算重力场模型高程异常,结合SRTM和DTM2006.0计算剩余地形模型,进而计算RTM高程异常,以弥补EGM2008重力场模型短波信号缺失的不足。从实测高程异常中减去重力场模型高程异常以及RTM高程异常,得到残余高程异常。利用曲面拟合方法拟合残余高程异常模型,解决高程基准系统偏差问题。通过实际GPS/水准数据计算证明,该方法能显著提高GPS点的高程异常计算精度,达到精化区域似大地水准面的目的。     

低低跟踪重力卫星星间距离变率对区域质量异常敏感度分析

推导了区域质量变化和球谐系数变化量之间的函数关系,根据该关系式可将不同区域质量异常对重力场变化的影响分离出来。设计仿真实验,分析时变重力场模型下星间距离变率相对稳态背景重力场模型下星间距离变率值的残差变化。结果表明,区域质量异常对星间距离变率的影响在空间分布上是有限的,但星间距离变率对区域质量异常具有很高的敏感性。     

��������̽�⼼��������������о�

????????????????￡??????????????????????????4??????????????????????????????y???????????????????????????????????????о???????????????????????????     

EGM96��EGM2008����������ģ�ͼ��㵯�����Ŷ������ıȽ�

????????????????????????????????????????????????????????????36???360???2 160????????????????????????????????????36???360???EGM96??EGM2008??????????????????????????????????????????????????????????????????????????????????????????????????????????EGM96??EGM2008???????????????????????????????????????????EGM2008????????????EGM96????     

����GRACEʱ������������й������ܱߵ�����ˮ�����±仯

????GRACE???????????????2005 ??1~11?μ??й?????????????????±仯???о???????????GRACE??????????????????????к???????????????????????cm???Ч???仯????????÷???????????????????仯?о??????????