利用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...     

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

本文在法方程层面融合GOCE卫星的V_(xx)、V_(yy)、V_(zz)和V_(xz)重力梯度分量观测数据和GRACE卫星观测数据,采用直接法解算了220阶次的重力场模型Tongji-GOGR2019S.首先利用IIR带通滤波器在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代模型.     

联合星载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时变重力场能力,为我国卫星重力场反演提供了重要的技术支持.     

融合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卫星重力场建模方法研究

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

GOCE卫星数据确定的全球稳态海面地形和海表地转流

利用GOCE卫星235天的观测数据恢复了200阶次的重力场模型SWJTU-GO01S,结合欧空局提供的最新GOCE重力场模型和CNES-CLS 2011平均海面高模型,计算了全球稳态海面地形和海表地转流,并采用GRACE模型、多源数据同化模型和海洋浮标观测数据对GOCE模型的计算结果进行对比分析.结果表明:由于重力场模型精度和分辨率较高,GOCE计算结果所需的滤波半径小于GRACE结果;GOCE和GRACE模型的计算结果与CNES-CLS09稳态海面地形差异的RMS分别为6 cm和7 cm左右;与海洋浮标实测数据对比发现,GOCE和GRACE的计算结果与实测数据差异明显,但GOCE的计算结果优于GRACE结果,而SWJTU-GO01S与DIR-R4和TIM-R4模型在全球范围内具有较好的一致性.整体而言,GOCE比GRACE数据的计算结果可以反映更小尺度地转流,且计算的精度更高;海洋环流结果和水准数据的对比表明SWJTU-GO01S与DIR-R4和TIM-R4模型的精度符合性较好,三者计算的地转流精度基本相当.     

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

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

海潮对卫星重力场恢复的影响

本文讨论了海潮对卫星重力测量的影响问题. 首先介绍了海潮对卫星重力测量影响的基本理论；采用FES02和TPXO6海潮模型计算了海潮负荷对卫星重力结果前60阶的影响；并用两个模型之间的差异作为海潮模型精度的估计量，据此计算了海潮模型误差对卫星重力结果的影响. 与GRACE恢复的重力场精度的比较说明：海潮对重力场40阶以下的影响都超过了目前重力场恢复精度；尽管由于卫星测高技术的发展，海潮模型的精度有了很大的提高，但目前的全球海潮模型用于GRACE重力场恢复的前12阶的改正还是不够精确. 另外，我们也利用中国东海和南海潮汐资料以及FES02海潮模型讨论了中国近海潮汐效应对GRACE观测的影响. 结果说明该影响与海潮模型的误差相当. 这反映了当前海潮模型的不确定度，因此通过结合全球验潮站资料有望提高海潮对卫星重力测量的改正精度.     

卫星重力探测技术的发展

在地球物理勘探领域中,人造地球卫星的发射为重力测量提供了新的途径。与以往探测重力的手段相比,重力卫星的发射大大改善了人们对地球重力场的认识,随着CHAMP、GRACE和GOCE卫星的发射,将把现有静态中长波部分重力场的精度提高1—2个数量级,并提供长波部分重力场随时间变化的信息。卫星重力学对我国的基础测绘服务和国防建设有着重要的实用价值。     

模拟分析低低跟踪模式重力卫星反演地球重力场的精度

本文利用卫星重力反演与模拟软件ANGELS系统(ANalyst of Gravity Estimation with Low-orbit Satellites)对低低跟踪模式的重力卫星的关键载荷精度指标进行了深入分析.模拟结果表明:(1)对短弧长积分法而言,在低低跟踪模式的关键载荷精度指标中,重力场反演精度对星间距离变率精度最为敏感;(2)通过对目前在轨运行GRACE的载荷指标进行分析,发现轨道数据的误差主要影响重力场的低阶部分(约小于25阶),较高阶次部分(约大于26阶)主要受星间距离变率的误差限制;(3)如果下一代低低跟踪模式的重力卫星的目标之一是把重力异常反演精度较GRACE提高约10倍,则在保持轨道高度和GRACE相同的前提下,轨道、星间距离变率和星载加速度计等关键载荷指标需要达到的最低精度分别约为2cm、10nm·s-1和3.0×10-10 m·s-2;(4)轨道精度和混频误差将是影响下一代低低跟踪模式重力卫星重力场恢复能力进一步提高的主要制约因素,距离变率精度和加速度计精度存在盈余.     

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 cm²·s^-1,平均值为22.3 cm²·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~600m深度范围内所引起的混合可达10-2.79 m2·s-1左右,比大洋的统计结果10-5 m2·s-1高出两个数量级以上.而地中海涡旋所引起的湍流混合率可达10-3.44 m2·s-1左右,与大洋统计结果相比高出1.5个数量级左右,并且地中海涡旋下边界的混合要强于上边界,这一特征与前人的研究一致,另外涡旋上边界之上以及侧边界的外侧也具有非常高的混合率.     

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

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