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

高精度静态卫星重力场模型在全球海洋环流研究、全球/区域数字高程基准面确定等领域有重要应用,本文研究仅利用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精度相当.

     

新一代GRACE重力卫星反演地球重力场的预期精度

基于低低卫卫跟踪模式,本文主要探讨利用动力学法融合精密轨道数据和星间测距或距离变率数据求解地球重力场的基本原理与方法,该方法既可对两颗低低跟踪卫星的初始状态误差进行有效校正,也可充分利用低轨卫星轨道所包含的低频重力场信息.为探讨适合我国国情的低低跟踪模式下的重力卫星指标,本文以不同星载设备精度指标的组合进行模拟计算,模拟结果显示:(1)把GRACE卫星的星间距离变率指标提高一个量级,其余指标保持与GRACE卫星设计指标一致时,可使地球重力场的精度获得同量级的提高;(2)若星间距离变率为1.0×10^-8 m·s^-1,轨道高度为300 km,加速度计精度为3.0×10^-10 m·s^-2,轨道精度为0.03 m, 星间距离100 km,与利用GRACE的设计指标反演出的重力场精度相比,可提高约121倍,并建议我国未来低低跟踪重力卫星计划参考此指标.     

融合GOCE和GRACE卫星数据的无约束重力场模型Tongji-GOGR2019S

本文在法方程层面融合GOCE卫星的V_xx、V_yy、V_zz和V_xz重力梯度分量观测数据和GRACE卫星观测数据,采用直接法解算了220阶次的重力场模型Tongji-GOGR2019S.首先利用ⅡR带通滤波器在5~41 mHz的重力梯度带宽范围内对约24个月的GOCE重力梯度观测方程进行无相移滤波处理,并组成解算220阶次重力场模型的法方程,各梯度分量根据相对于参考模型统计精度进行定权;然后与13.5 a GRACE数据建立的180阶次Tongji-Grace02s重力场模型的法方程进行叠加,解算了220阶次的无约束纯卫星重力场模型Tongji-GOGR2019S.利用EIGEN-6C4重力场模型、GNSS/水准数据、DTU15重力异常数据以及欧洲区域似大地水准面模型EGG2015等数据对Tongji-GOGR2019S模型精度进行全面的检核评定,结果表明：引入GOCE卫星梯度数据后,高于72阶的位系数精度优于Tongji-Grace02s模型,Tongji-GOGR2019S模型的整体精度接近同阶次的DIR-R6等GOCE卫星第6代模型.

     

融合GOCE和GRACE卫星数据的无约束重力场模型Tongji-GOGR2019S

本文在法方程层面融合GOCE卫星的V_xx、V_yy、V_zz和V_xz重力梯度分量观测数据和GRACE卫星观测数据,采用直接法解算了220阶次的重力场模型Tongji-GOGR2019S.首先利用ⅡR带通滤波器在5~41 mHz的重力梯度带宽范围内对约24个月的GOCE重力梯度观测方程进行无相移滤波处理,并组成解算220阶次重力场模型的法方程,各梯度分量根据相对于参考模型统计精度进行定权;然后与13.5 a GRACE数据建立的180阶次Tongji-Grace02s重力场模型的法方程进行叠加,解算了220阶次的无约束纯卫星重力场模型Tongji-GOGR2019S.利用EIGEN-6C4重力场模型、GNSS/水准数据、DTU15重力异常数据以及欧洲区域似大地水准面模型EGG2015等数据对Tongji-GOGR2019S模型精度进行全面的检核评定,结果表明：引入GOCE卫星梯度数据后,高于72阶的位系数精度优于Tongji-Grace02s模型,Tongji-GOGR2019S模型的整体精度接近同阶次的DIR-R6等GOCE卫星第6代模型.     

卫星重力梯度观测数据L1级构建方法

高精度重力梯度观测数据L1级构建的系统方法是推进我国自主重力卫星任务重要的基础数据处理技术.本文以GOCE卫星L1级数据预处理技术和关键载荷原始数据为参考,面向我国发展的梯度测量卫星的任务需要,系统研究并初步实现了卫星重力梯度观测数据L1级构建方法,主要包括加速度计电压数据转换、多星敏感器联合姿态数据的角速度重建、卫星重力梯度分量构建等技术内容.计算结果表明,加速度计超灵敏轴精度为10^-10~10^-11m·s^-2·Hz^-1/2,达到重力梯度仪设计精度要求;多星敏感器联合解算最佳姿态角速度w_y、w_z在10~100 mHz内精度约提升1个量级,其精度约达到10^-5 rad·s^-1·Hz^-1/2量级,能够有效抑制低精度角速度分量在坐标系转换中导致的噪声传播;基于维纳滤波方法恢复的角速度在5~100 mHz频段内的平方根功率谱密度提升了（5.21~6.56）×10^-11 rad·s^-1·Hz^-1/2,显示了基于高精度角速度解算重力梯度分量的必要性;构建重力梯度各分量计算值与全球重力场和海洋环流探测器（GOCE）官方公布的重力梯度分量精度相当,其梯度张量的迹在20~100 mHz频段范围内约为10 mE·Hz^-1/2,验证了本文构建方法的有效性.研究工作可为下一步我国推进实施民用重力梯度测量卫星任务提供自主的原始数据处理技术支撑与储备.

     

青藏高原东北缘ALOS AW3D30高程模型和EIGEN-6C4自由空气重力异常模型精度分析

利用EGM96全球重力场模型对ALOS AW3D30全球数字高程模型的高程系统进行转换,结合实测GNSS数据,分析ALOS AW3D30模型在青藏高原东北缘地区的精度。利用实测重力数据计算青藏高原东北缘地区的自由空气异常,并与EGIEN-6C4自由空气异常模型进行对比,分析其精度。结果表明：青藏高原东北缘地区ALOS AW3D30模型与264个GNSS测点的差异标准差为3.4 m,该地区ALOS AW3D30精度优于4 m;青藏高原东北缘地区自由空气重力异常整体呈负值,局部为正值,变化范围为-177×10^-5～166×10^-5 m·s^-2,实测重力异常与模型结果空间分布基本一致,存在局部差异;EIGEN-6C4自由空气重力异常变化范围为-163×10^-5～142×10^-5 m·s^-2,实测结果与模型结果差异为-119×10^-5～63×10^-5 m·s^-2,平均值为-20×10^-5<...     

单加速度计模式下的GOCE卫星重力场建模方法研究

GOCE卫星由于加速度计的特殊安装方式,其非保守力主要由普通模式的组合加速度提供,使得单个加速度计的特征更难提取.本文首次采用实测数据,研究了单加速度计模式下的高低跟踪数据处理.利用GOCE任务2009年(2009-11—2009-12)的实测数据,分别以GOCE卫星梯度仪坐标系三个坐标轴正向的加速度计为研究对象,利用1 s间隔的高采样轨道数据,采用动力法同时进行卫星重力场建模和加速度计的精密校准.为了克服两极地区的数据缺失对重力场模型低次系数的影响,即所谓的极空白问题,引入同期GRACE卫星的观测数据,采用方差分量估计方法,建立了GRACE/GOCE卫星跟踪卫星重力场模型WHU-GRGO-SST.该模型完全到100阶次,经6169个美国GPS水准点数据检验,在同阶次上与EGM2008和GGM05S的精度水平相同.分析发现,GOCE卫星的加速度计偏差参数存在显著的漂移,也显示了单加速度计模式处理GOCE高低跟踪数据的优势.本文的研究成果为建立静态高分辨率、高精度的GRACE/GOCE重力场模型提供了更严密的模型与技术方案,同时也为GOCE卫星梯度仪校准,以及梯度数据的深入分析提供了重要的参考信息.     

基于残余星间速度法精确和快速反演下一代GRACE Follow-On地球重力场

基于新型残余星间速度法（RIRM）反演了120阶GRACE Follow-On地球重力场. 第一,由于GPS定轨精度相对较低,通过将激光干涉测距仪的高精度残余星间速度（测量精度10^-7 m·s^-1）引入残余轨道速度差分矢量的视线分量构建了新型RIRM观测方程. 第二,基于2点、4点、6点和8点RIRM公式对比论证了最优的插值点数. 如果相关系数和采样间隔一定,随着插值点数的增加,卫星观测值的信号量被有效加强,而卫星观测值的误差量也同时增加. 因此,6点RIRM公式是提高下一代地球重力场精度的较优选择. 第三,相关系数对地球重力场精度的影响在不同频段表现为不同特性. 随着相关系数的逐渐增大,地球长波重力场精度逐渐降低,而地球中长波重力场精度逐渐升高. 第四,基于6点RIRM公式,通过30天观测数据和采样间隔5 s,分别利用星间速度和残余星间速度观测值,在120阶次处反演下一代GRACE Follow-On累计大地水准面精度为1.638×10^-3 m和1.396×10^-3 m. 研究结果表明：（1）残余星间速度观测量较星间速度对地球重力场反演精度更敏感;（2）GRACE Follow-On地球重力场精度较GRACE至少高10倍.     

利用解析法有效快速估计将来GRACE Follow-On地球重力场的精度

本文首次利用解析法有效快速估计了将来GRACE（Gravity Recovery and Climate Experiment） Follow-On地球重力场的精度. 第一,基于功率谱原理分别建立了新的GRACE Follow-On卫星激光干涉星间测量系统星间速度、GPS接收机轨道位置和轨道速度以及加速度计非保守力误差影响累计大地水准面的单独和联合解析误差模型. 第二,利用提出的GRACE卫星关键载荷匹配精度指标和美国喷气推进实验室（JPL）公布的GRACE Level 1B实测精度指标的一致性,以及估计的GRACE累计大地水准面精度和德国波兹坦地学研究中心（GFZ）公布的EIGEN-GRACE02S地球重力场模型实测精度的符合性,验证了本文建立的解析误差模型是可靠的. 第三,论证了GRACE Follow-On卫星不同关键载荷匹配精度指标和轨道高度对地球重力场精度的影响. 在360阶处,利用轨道高度250 km、星间距离50 km、星间速度误差1×10^-9m/s、轨道位置误差3×10^-5m、轨道速度误差3×10^-8m/s和非保守力误差3×10^-13m/s²,基于联合解析误差模型估计累计大地水准面的精度为1.231×10^-1 m. 本文的研究不仅为当前GRACE和将来GRACE Follow-On地球重力场精度的有效快速确定提供了理论基础和计算保证,同时对国际将来GRAIL（Gravity Recovery and Interior Laboratory）月球卫星重力测量计划的成功实施具有重要的参考意义.     

联合星载GPS和KBRR星间测速数据反演GRACE时变重力场模型

高精度GRACE卫星时变重力场反演一直是卫星重力测量中的难题.为了恢复高精度的时变地球重力场模型,本文联合GRACE卫星的星载GPS和KBR星间测速观测数据,在对GRACE卫星进行精密定轨的同时,解算出60阶月平均地球重力场模型.通过对GRACE卫星的定轨精度、星载GPS相位和KBR星间测速数据的拟合残差以及时变地球重力场模型解算精度等分析,表明:(1)与美国宇航局喷气推进实验室(JPL)发布的约化动力学精密轨道相比,本文确定GRACE卫星轨道三维位置误差小于5 cm.(2)星载GPS相位数据拟合残差为5~8 mm,KBR星间测速数据拟合残差为0.18~0.30μm·s~(-1).(3)解算的月平均重力场模型与美国德克萨斯大学空间研究中心(CSR)、德国地学研究中心(GFZ)和JPL发布的RL05模型精度接近,时变信号在全球范围内具有很好的空间分布一致性.通过计算亚马逊流域和长江流域的水储量变化,本文与上述三个机构的计算结果无明显差异,且相关系数均达0.9以上.可见,本文建立的卫星轨道与重力场同解算法具有反演高精度GRACE时变重力场能力,为我国卫星重力场反演提供了重要的技术支持.     

On High-Frequency Topography-Implied Gravity Signals for a Height System Unification Using GOCE-Based Global Geopotential Models

National height reference systems have conventionally been linked to the local mean sea level, observed at individual tide gauges. Due to variations in the sea surface topography, the reference levels of these systems are inconsistent, causing height datum offsets of up to ±1–2 m. For the unification of height systems, a satellite-based method is presented that utilizes global geopotential models (GGMs) derived from ESA’s satellite mission Gravity field and steady-state Ocean Circulation Explorer (GOCE). In this context, height datum offsets are estimated within a least squares adjustment by comparing the GGM information with measured GNSS/leveling data. While the GNSS/leveling data comprises the full spectral information, GOCE GGMs are restricted to long wavelengths according to the maximum degree of their spherical harmonic representation. To provide accurate height datum offsets, it is indispensable to account for the remaining signal above this maximum degree, known as the omission error of the GGM. Therefore, a combination of the GOCE information with the high-resolution Earth Gravitational Model 2008 (EGM2008) is performed. The main contribution of this paper is to analyze the benefit, when high-frequency topography-implied gravity signals are additionally used to reduce the remaining omission error of EGM2008. In terms of a spectral extension, a new method is proposed that does not rely on an assumed spectral consistency of topographic heights and implied gravity as is the case for the residual terrain modeling (RTM) technique. In the first step of this new approach, gravity forward modeling based on tesseroid mass bodies is performed according to the Rock–Water–Ice (RWI) approach. In a second step, the resulting full spectral RWI-based topographic potential values are reduced by the effect of the topographic gravity field model RWI_TOPO_2015, thus, removing the long to medium wavelengths. By using the latest GOCE GGMs, the impact of topography-implied gravity signals on the estimation of height datum offsets is analyzed in detail for representative GNSS/leveling data sets in Germany, Austria, and Brazil. Besides considerable changes in the estimated offset of up to 3 cm, the conducted analyses show that significant improvements of 30–40% can be achieved in terms of a reduced standard deviation and range of the least squares adjusted residuals.     

Comparison of global geopotential models from the champ and grace missions for regional geoid modelling in Turkey

The continuous efforts on establishment and modernization of the geodetic control in Turkey include a number of regional geoid models that have been determined since 1976. The recently released gravimetric Geoid of Turkey, TG03, is used in geodetic applications where GPS-heights need to be converted to the local vertical datum. To reach a regional geoid model with improved accuracy, the selection of the appropriate global geopotential model is of primary importance. This study assesses the performance of a number of recent satellite-only and combined global geopotential models (GGMs) derived from CHAMP and GRACE missions’ data in comparison to the older EGM96 model, which is the underlying reference model for TG03. In this respect, gravity anomalies and geoid heights from the global geopotential models were compared with terrestrial gravity data and low-pass filtered GPS/levelling data, respectively. Also, five new gravimetric geoid models, computed by the Fast Fourier Transform technique using terrestrial gravity data and the geopotential models, were validated at the GPS/levelling benchmarks. The findings were also compared with the validation results of the TG03 model. The tests showed that as it was expected any of the high-degree combined models (EIGEN-CG03C, EIGEN-GL04C, EGM96) can be employed for determining the gravity anomalies over Turkey. In the west of Turkey, EGM96 and EIGEN-CHAMP03S fit the GPS/levelling surface better. However, all the tested GGMs revealed equal performance when they were employed in gravimetric geoid modelling after de-trending the gravimetric geoid model with corrector surface fitting. The new geoid models have improved accuracy (after fit) compared to TG03.     

基于GRACE Follow-On卫星重力梯度法精确反演地球重力场

由于GRACE Follow-On双星系统等效于基线长为星间距离的一维水平重力梯度仪,因此本文基于GRACE Follow-On卫星重力梯度法开展了精确和快速反演下一代地球重力场的可行性论证研究. 研究结果表明：第一,基于GRACE Follow-On卫星重力梯度法（GFO-SGGM）,利用卫星轨道参数（轨道高度250 km、星间距离50 km、轨道倾角89°、轨道离心率0.001）、关键载荷测量精度（星间距离10-6 m、星间速度10-7 m·s-1、星间加速度10-10 m·s-2、轨道位置10-3 m、轨道速度10-6 m·s-1、非保守力10-11 m·s-2）、观测时间30天和采样间隔10 s反演了120阶地球重力场,在120阶处累计大地水准面精度为9.331×10-4 m. 第二,在120阶内,利用将来GRACE Follow-On双星反演地球重力场精度较现有GRACE双星平均提高61倍,因此GRACE Follow-On卫星重力梯度法是进一步提高地球重力场反演精度的优选方法. 第三,下一代GRACE Follow-On计划较当前GRACE计划的优点如下：轨道高度更低（200~300 km）、载荷精度更高（10-7 ~10-9 m·s-1）和星间距离更短（50~100 km）.     

A comparison of GOCE gravitational models with EGM2008

Four new gravity field models from GOCE, two of them combined with GRACE, are compared here with EGM2008. The objectives are to look into the differences in consecutive ranges of the spherical harmonic expansion globally as well as in selected geographical regions and in the regions of the various data sources used for EGM2008. In general, GOCE is able to contribute to improved global gravity models in the spherical harmonic range between 120 and 200 (and above). The agreement between EGM2008 and the GOCE models is very good in well-surveyed regions such as North America, Europe and Australia, with geoid RMS-differences on the order of 4–6 cm. In other regions, where the surface gravity data available for the development of EGM2008 were poor, such as South America, Africa, South-East Asia or China the RMS-differences are on a level of 30 cm. Here GOCE leads to a significant improvement. These findings are confirmed by the analysis of the areas of the various EGM2008 data sources. In the regions of the so-called “fill-in” data of EGM2008 RMS-geoid height differences are high. In Antarctica GOCE also gives important improvements in terms of spatial resolution and accuracy. In general, the agreement between EGM2008 and the GOCE-models up to degree and order (d/o) 200 is good, with a global (excluding the polar gaps of GOCE orbits, throughout) geoid difference RMS of 11 cm, in the ocean areas 8 cm and 20 cm in the continental areas. GOCE models are better suited for ocean circulation studies because no prior ocean information enters into the data reduction process, as it is the case when deducing gravity anomalies from an altimetric mean sea surface. On the other hand, the good consistency between GOCE-models and EGM2008 in ocean areas very likely indicates that the influence of ocean circulation information on EGM2008 is rather small. The four tested GOCE models behave similarly except at the highest latitudes where GOCE lacks data due to its orbit inclination of 96.5° and some form of regularization which has to be applied.     

基于镭同位素分布的黄海和东海垂直混合速率计算

在黄海和东海采样测定了水体中的镭同位素分布,用平流扩散模型描述镭同位素分布,最小二乘方法计算了垂直涡动扩散系数和上升流或下降流流速.结果给出北黄海中部、南黄海中部、浙江沿岸和台湾北部海域存在上升流,流速分别为0.46×10-3cm·s-1、0.17×10-3~1.39×10-3cm·s-1、2.02×10-3~3.04×10-3 cm·s-1和1.06×10-3~2.51×10-3 cm·s-1.北黄海中部和东海东北部存在下降流.流速分别为-2.30×10-3 cm·s-1和-0.61×10-3~-2.10×10-3 cm·s-1.计算同时给出的垂直涡动扩散系数为5.84~48.2 cm2·s-1,平均值为22.3 cm2·s-1.北黄海和浙江沿岸上升流流速与文献的结果一致;北黄海中部存在下降流与文献的结论一致.本研究结果与文献结果一致是对所建立的方法的肯定,也是对文献研究结果的支持.     

RTM-based omission error corrections for global geopotential models: Case study in Central Europe

The aim of this paper is to evaluate the effects of residual terrain model (RTM) on potential and on gravity and to point out how significant can the omission error of global geopotential models (GGMs) be and how it can influence their testing. The RTM for Central Europe is computed in the spherical approximation. The topography is modelled by spherical tesseroids and the gravitational effect of the topography is obtained as a sum of their partial gravitational effects. A detailed picture of RTM in Slovakia is shown. The testing of GOCE (Gravity Field and Steady-State Ocean Circulation Explorer) global geopotential models in Central Europe published earlier is re-evaluated with the more rigorous omission error estimation. Experimental results show significantly better agreement between the gravity anomalies computed from global geopotential models with the omission-error estimation and gravity anomalies obtained from the direct measurements. On the other hand, for height anomalies such an improvement is not observed. The results are discussed in context of the other previously published studies.     

基于时空域混合法利用Kaula正则化精确和快速解算GOCE地球重力场

为了研究卫星重力梯度技术对中高频地球重力场反演精度的影响,本文基于时空域混合法,利用Kaula正则化反演了250阶GOCE地球重力场.模拟结果表明:第一,时空域混合法是精确和快速求解高阶地球重力场的有效方法;第二,Kaula正则化是降低正规阵病态性的重要方法;第三,基于改进的预处理共轭梯度迭代法可快速求解大型线性方程组...     

地震海洋学方法在海洋混合参数提取中的研究与应用——以南海内波和地中海涡旋为例

混合是海洋中普遍存在的一种海水运动形式,对多个海洋学分支的研究具有重要的影响.随着物理海洋学的研究重心从大尺度向中小尺度现象过渡,近年来混合问题的研究重心也逐渐转向了中小尺度现象.内波与中尺度涡都是非常重要的中小尺度物理海洋学现象,对海洋能量在不同尺度中的级联发挥着重要的作用.本文基于地震海洋学研究了海洋混合参数的提取方法,并以南海内波和地中海涡旋为例进行了计算和分析.结果显示,南海内波在200~600 m深度范围内所引起的混合可达10-2.79 m2·s-1左右,比大洋的统计结果10-5 m2·s-1高出两个数量级以上.而地中海涡旋所引起的湍流混合率可达10-3.44m2·s-1左右,与大洋统计结果相比高出1.5个数量级左右,并且地中海涡旋下边界的混合要强于上边界,这一特征与前人的研究一致,另外涡旋上边界之上以及侧边界的外侧也具有非常高的混合率.     

三门湾沿海声层析潮流观测实验

2009年9月6日至9日在三门湾进行了沿海声层析(Coastal Acoustic Tomography,CAT)潮流观测实验.实验由7台沿海声层析仪组网进行,并分别由渔船定点抛锚于7个站位.实验期间,还进行了定点ADCP(Acoustic Doppler Current Profiler)观测.通过建立逆模式对声传播时间差进行解析,引入权重因子,用L-curve法确定阻尼因子的最佳值,继而根据阻尼最小二乘法得到流速的最佳解.根据逆模式得到的流速分布可知该海区的潮流以半日潮(M2)为主,M2潮流椭圆呈东南-西北走向,潮流基本都是顺着水道,即涨潮为西北流向,退潮为东南流向.西北向与东南向最大流速分别为1.03 m·s-1和1.09 m·s-1.实验期间该区域的余流是从湾外流入湾内,平均流速约为0.05 m·s-1.CAT与定点ADCP流速的东分量和北分量的均方差均小于0.18 m·s-1.这样大面积的潮流和余流水平分布的同步观测,用传统观测手段很难实现.通过以上结果可以得出,沿海声层析技术可以作为一种新的测流方法对强潮海区进行大面积潮流观测,可在我国沿海的海洋环境监测等方面发挥重要作用.