利用GOCE和GRACE卫星观测数据确定静态重力场模型  

Determination of static gravity field model by using satellite data of GOCE and GRACE

作　　者：赵永奇[1,2] 李建成 徐新禹[1,2,4] 苏勇 踪华[5] 魏辉 ZHAO YongQi;LI JianCheng;XU XinYu;SU Yong;ZONG Hua;WEI Hui(School of Geodesy and Geomatics,Wuhan University,Wuhan 430079,China;Key Laboratory of Geospace Environment and Geodesy,Ministry of Education,Wuhan University,Wuhan 430079,China;School of Civil Engineering and Geomatics,Southwest Petroleum University,Chengdu610500,China;Hubei Luojia Laboratory,Wuhan 430079,China;State Key Laboratory of Science and Technology on Aerospace Intelligent Control,Beijing100854,China)

机构地区：[1]武汉大学测绘学院,武汉430079 [2]武汉大学地球空间环境与大地测量教育部重点实验室,武汉430079 [3]西南石油大学土木工程与测绘学院,成都610500 [4]湖北珞珈实验室,武汉430079 [5]宇航智能控制技术国家级重点实验室,北京100854

出　　处：《地球物理学报》2023年第6期2322-2336,共15页Chinese Journal of Geophysics

基　　金：国家自然科学基金(42192533,41721003,42074019,41774020,41804077);湖北珞珈实验室开放基金资助项目(220100045);四川省自然科学基金项目(2022NSFSC1047)联合资助。

摘　　要：高精度静态卫星重力场模型在全球海洋环流研究、全球/区域数字高程基准面确定等领域有重要应用,本文研究仅利用GOCE卫星和联合GRACE卫星观测数据确定高精度高阶次静态重力场模型.利用GOCE卫星全周期高精度引力梯度分量(V_(xx)、V_(yy)、V_(zz)和V_(xz))观测值基于直接最小二乘法构建300阶次的SGG(Satellite Gravity Gradiometry)法方程,并利用卫星跟踪卫星观测值基于点域加速度法构建130阶SST(Satellite-to-Satellite Tracking)法方程,然后利用方差分量估计联合SGG和SST法方程确定300阶次纯GOCE卫星重力场模型GOSG02S.利用全周期GRACE观测数据由动力学方法解算了180阶次的SWPU-GRACE2021S模型,并将其对应法方程与GOCE卫星法方程联合解算了GRACE和GOCE的联合模型WHU-SWPU-GOGR2022S.分别基于XGM2019模型和GPS水准数据对本文解算的三个模型GOSG02S、SWPU-GRACE2021S和WHU-SWPU-GOGR2022S在频域和空域进行了精度分析,结果表明,GOSG02S和WHU-SWPU-GOGR2022S模型与GO_CONS_GCF_2_DIR_R6、GO_CONS_GCF_2_TIM_R6、GO_CONS_GCF_2_SPW_R5、GOCO06s和Tongji-GMMG2021S等使用了GOCE卫星全周期数据的模型精度相当,精度差异基本都在毫米量级;SWPU-GRACE2021S模型在160阶次之前与国际主流GRACE卫星重力场模型ITSG-Grace2018s和Tongji-Grace02s精度相当.The high-precision static satellite gravity field models have important applications in fields such as global ocean circulation research and global/regional digital elevation datum determination.In this paper,we discussed the determination of high-degree static satellite gravity field models with GOCE observation,GRACE observation and the joint of them.We first constructed a 300-degree Satellite Gravity Gradiometry(SGG)normal equation with the high-precision gravitational gradient components V_(xx),V_(yy),V_(zz)and V_(xz)throughout the entire mission of GOCE by the Direct least squares method,and a 130-degree Satellite-to-Satellite Tracking(SST)normal equation with the SST observation data by the point-wise acceleration approach.The 300-degree GOCE-only satellite gravity field model GOSG02Sis determined by combining SGG and SST normal equations with variance component estimation.The 180-degree model SWPU-GRACE2021Sis then determined based on the dynamic approach with the GRACE data throughout the entire mission cycle of 15years,and the normal equation is combined with GOCE normal equation to determine WHU-SWPU-GOGR2022S,ajoint model of GOCE and GRACE.Finally,the XGM2019model and GPS/leveling data are used for precision analysis of GOSG02S,SWPU-GRACE2021Sand WHU-SWPU-GOGR2022Sin frequency domain and space domain respectively.The results show that the accuracy of GOSG02Sand WHU-SWPU-GOGR2022Sis comparable to GO_CONS_GCF_2_DIR_R6,GO_CONS_GCF_2_TIM_R6,GO_CONS_GCF_2_SPW_R5,GOCO06sand Tongji-GMMG2021Sthat use the entire mission data of GOCE satellite and the accuracy differences are in the order of millimeters.The SWPU-GRACE2021S model has the same accuracy below degree/order 160as the international mainstream of GRACE satellite gravity field models,i.e.,ITSG-Grace2018sand Tongji-Grace02s.

关 键 词：GOCE GRACE 直接最小二乘法 加速度法 ARMA 正则化 

分 类 号：P223[天文地球—大地测量学与测量工程]