Long-wavelength global gravity field models: GRIM4-S4, GRIM4-C4

Summary.  GFZ Potsdam and GRGS Toulouse/Grasse jointly developed a new pair of global models of the Earth's gravity field to satisfy the requirements of the recent and future geodetic and altimeter satellite missions. A precise gravity model is a prerequisite for precise satellite orbit restitution, tracking station positioning and altimeter data reduction. According to different applications envisaged, the new model exists in two parallel versions: the first one being derived exclusively from satellite tracking data acquired on 34 satellites, the second one further incorporating satellite altimeter data over the oceans and terrestrial gravity data. The most recent “satellite-only” gravity model is labelled GRIM4-S4 and the “combined” gravity model GRIM4-C4. The models are solutions in spherical harmonics and have a resolution up to degree and order 60 plus a few resonance terms in the case of GRIM4-S4, and up to degree/order 72 in the case of GRIM4-C4, corresponding to a spatial resolution of 555 km at the Earth's surface. The gravitational coefficients were estimated in a rigorous least squares adjustment simultaneously with ocean tidal terms and tracking station position parameters, so that each gravity model is associated with a consistent ocean tide model and a terrestrial reference frame built up by over 300 optical, laser and Doppler tracking stations. Comprehensive quality tests with external data and models, and test arc computations over a wide range of satellites have demonstrated the state-of-the-art capabilities of both solutions in long-wavelength geoid representation and in precise orbit computation. Received 1 February 1996; Accepted 17 July 1996     

Evaluation of pre-CHAMP gravity models `GRIM5-S1' and `GRIM5-C1' with satellite crossover altimetry

 The new GFZ/GRGS gravity field models GRIM5-S1 and GRIM5-C1, currently used as initial models for the CHAMP mission, have been compared with other recent models (JGM 3, EGM 96) for radial orbit accuracy (by means of latitude lumped coefficients) in computations on altimetry satellite orbits. The bases for accuracy judgements are multi-year averages of crossover sea height differences from Geosat and ERS 1/2 missions. This radially sensitive data is fully independent of the data used to develop these gravity models. There is good agreement between the observed differences in all of the world's oceans and projections of the same errors from the scaled covariance matrix of their harmonic geopotential coefficients. It was found that the tentative scale factor of five for the formal standard deviations of the harmonic coefficients of the new GRIM fields is justified, i.e. the accuracy estimates, provided together with the GRIM geopotential coefficients, are realistic. Received: 20 February 2001 / Accepted: 24 October 2001     

A comparison and analysis of airborne gravimetry results from two strapdown inertial/DGPS systems

In September 1996 the University of Calgary tested a combination of strapdown inertial navigation systems and differential global positioning system (DGPS) receivers for their suitability to determine gravity at aircraft flying altitudes. The purpose of this test was to investigate the long-term accuracy and repeatability of the system, as well as its potential for geoid and vertical gradient of gravity determination. The test took place during a 3-day period in the Canadian Rocky Mountains over a single 100 × 100 km area which was flown with 10-km line spacing. Two flights were done at 4350 m in E–W and N–S profile directions, respectively, and one at 7300 m with E–W profiles. Two strapdown inertial systems, the Honeywell LASEREF III and the Litton-101 Flagship, were flown side by side. Comparison of the system estimates with an upward-continued reference showed root-mean-square (RMS) agreement at the level of 3.5 mGal for 90- and 120-s filter lengths. The LASEREF III, however, performed significantly better than the Litton 101 for shorter filtering periods of 30 and 60 s. A comparison between the two systems results in an RMS agreement of 2.8 and 2.3 mGal for the 90- and 120-s filters. The better agreement between the two systems is mainly due to the fact that the upward-continued reference has not been filtered identically to the system gravity disturbance estimates. Additional low-frequency differences seem to point to an error in the upward-continued reference. Finally, an analysis of crossover points between flight days for the LASEREF III shows a standard deviation of 1.6 mGal, which is near the noise level of the INS and GPS data. Further improvements to the system are possible, and some ideas for future work are briefly presented. Received: 17 March 1998 / Accepted: 1 February 1999     

A new approach for airborne vector gravimetry using GPS/INS

A new method for airborne vector gravimetry using GPS/INS has been developed and the results are presented. The new algorithm uses kinematic accelerations as updates instead of positions or velocities, and all calculations are performed in the inertial frame. Therefore, it is conceptually simpler, easier, more straightforward and computationally less expensive compared to the traditional approach in which the complex navigation equations should be integrated. Moreover, it is a unified approach for determining all three vector components, and no stochastic gravity modeling is required. This approach is based on analyzing the residuals from the Kalman filter of sensor errors, and further processing with wavenumber coefficient filterings is applied in case closely parallel tracks of data are available. An application to actual test-flight data is performed to test the validity of the new algorithm. The results yield an accuracy in the down component of about 3–4 mGal. Also, comparable results are obtained for the horizontal components with accuracies of about 6 mGal. The gravity modeling issue is discussed and alternative methods are presented, none of which improves on the original approach. Received: 18 April 2000 / Accepted: 14 August 2000     

Unification of vertical datums by GPS and gravimetric geoid models with application to Fennoscandia

The second Baltic Sea Level (BSL) GPS campaign was run for one week in June 1993. Data from 35 tide gauge sites and five fiducial stations were analysed, for three fiducial stations (Onsala, Mets?hovi and Wettzell) fixed at the ITRF93 system. On a time-scale of 5 days, precision was several parts in 10⁹ for the horizontal and vertical components. Accuracies were about 1 cm in comparison with the International GPS Geodynamical Service (IGS) coordinates in three directions. To connect the Swedish and the Finnish height systems, our numerical application utilises three approaches: a rigorous approach, a bias fit and a three-parameter fit. The results between the Swedish RH70 and the Finnish N 60 systems are estimated to −19.3 ± 6.5, −17 ± 6 and −15 ± 6 cm, respectively, by the three approaches. The results of the three indirect methods are in an agreement with those of a direct approach from levelling and gravity measurements. Received: 3 April 1996 / Accepted: 4 August 1997     

基于Laplace方程的垂线偏差法反演全球海域重力异常

海洋重力场模型反演的质量主要依赖于采用测高数据的精度、空间分辨率和数据分布密集程度。本文联合Geosat GM/ERM、ERS-1 GM/ERM、TOPEX/Poseidon、Envisat、Cryosat-2、Jason-1 ERM/GM和SARAL/AltiKa等多种测高观测数据集,深入比较了多种波形重跟踪算法的效果,回波数据重跟踪处理后的沿轨海面高标准差。统计表明,Sandwell算法优于MLE-4算法、Davis阈值法、改进阈值法和β参数拟合法;基于不同测高数据波形重采样的结果给出了沿轨海面梯度计算中低通滤波的参数选择方法,并采用Sandwell提出的垂线偏差法,反演了全球海域1′×1′的重力场模型。检核表明,反演结果与DTU13和SIO V23.1模型检核的差值均方根分别为3.4、1.8 mGal,与NGDC船测数据的检核精度为4~8 mGal,且本文模型在部分典型海区内精度更优。     

基于空域最小二乘法求解GOCE卫星重力场的模拟研究

本文论述了最小二乘过程中有色噪声的处理方法,提出使用AR模型对GOCE梯度观测值中的有色噪声进行时域滤波,数值模拟结果验证了该方法的有效性。利用数值模拟验证了直接求逆方法和PCCG法求解大型法方程的有效性,后者的效率远远高于前者。联合加入噪声（有色噪声和白噪声）的卫星重力梯度张量径向分量观测值Vzz和SST观测值,分别使用空域最小二乘法和SA方法恢复了180阶全球重力场模型,前者求解重力场模型的大地水准面和重力异常在180阶次的精度分别为3.01cm和0.75mGal,优于SA方法求解模型的精度。     

On some problems of the downward continuation of the 5′×5′ mean Helmert gravity disturbance

This research deals with some theoretical and numerical problems of the downward continuation of mean Helmert gravity disturbances. We prove that the downward continuation of the disturbing potential is much smoother, as well as two orders of magnitude smaller than that of the gravity anomaly, and we give the expression in spectral form for calculating the disturbing potential term. Numerical results show that for calculating truncation errors the first 180^∘ of a global potential model suffice. We also discuss the theoretical convergence problem of the iterative scheme. We prove that the 5^′×5^′ mean iterative scheme is convergent and the convergence speed depends on the topographic height; for Canada, to achieve an accuracy of 0.01 mGal, at most 80 iterations are needed. The comparison of the “mean” and “point” schemes shows that the mean scheme should give a more reasonable and reliable solution, while the point scheme brings a large error to the solution. Received: 19 August 1996 / Accepted: 4 February 1998     

Precise geoid determination over Sweden using the Stokes–Helmert method and improved topographic corrections

 Four different implementations of Stokes' formula are employed for the estimation of geoid heights over Sweden: the Vincent and Marsh (1974) model with the high-degree reference gravity field but no kernel modifications; modified Wong and Gore (1969) and Molodenskii et al. (1962) models, which use a high-degree reference gravity field and modification of Stokes' kernel; and a least-squares (LS) spectral weighting proposed by Sj?berg (1991). Classical topographic correction formulae are improved to consider long-wavelength contributions. The effect of a Bouguer shell is also included in the formulae, which is neglected in classical formulae due to planar approximation. The gravimetric geoid is compared with global positioning system (GPS)-levelling-derived geoid heights at 23 Swedish Permanent GPS Network SWEPOS stations distributed over Sweden. The LS method is in best agreement, with a 10.1-cm mean and ±5.5-cm standard deviation in the differences between gravimetric and GPS geoid heights. The gravimetric geoid was also fitted to the GPS-levelling-derived geoid using a four-parameter transformation model. The results after fitting also show the best consistency for the LS method, with the standard deviation of differences reduced to ±1.1 cm. For comparison, the NKG96 geoid yields a 17-cm mean and ±8-cm standard deviation of agreement with the same SWEPOS stations. After four-parameter fitting to the GPS stations, the standard deviation reduces to ±6.1 cm for the NKG96 geoid. It is concluded that the new corrections in this study improve the accuracy of the geoid. The final geoid heights range from 17.22 to 43.62 m with a mean value of 29.01 m. The standard errors of the computed geoid heights, through a simple error propagation of standard errors of mean anomalies, are also computed. They range from ±7.02 to ±13.05 cm. The global root-mean-square error of the LS model is the other estimation of the accuracy of the final geoid, and is computed to be ±28.6 cm. Received: 15 September 1999 / Accepted: 6 November 2000     

Prediction of surface horizontal displacements, and gravity and tilt changes caused by filling the Three Gorges Reservoir

 Horizontal displacements, and gravity and tilt changes induced by filling the Three Gorges Reservoir are modeled using elastic loading Green functions. When the water surface reaches its highest level, the effects become maximum on the reservoir banks. The longitudinal and latitudinal components of the horizontal displacements reach −8.2 and 7.7 mm respectively, gravity is increased by up to 3.4 mGal, and the prime vertical and meridian components of the tilt changes are −7.8 and −17.5 arcseconds respectively. Accordingly, the filling of the reservoir will influence values observed from global positioning system (GPS), gravimetry and tilt measurements in the area. The results given can be used to provide important corrections for extracting earthquake-related signals from observed data. Received: 19 January 2001 / Accepted: 3 September 2001     

星载原子干涉技术用于地球重力场测量及其精度评估

重力梯度卫星GOCE通过搭载静电式重力梯度仪,将全球静态重力场恢复至200阶以上。目前GOCE卫星已结束寿命,亟须发展下一代更高分辨率的卫星重力梯度测量来完善200~360阶的全球静态重力场模型。原子干涉型的重力梯度测量在空间微重力环境下可获得较长的干涉时间,因此具有很高的星载测量精度,是下一代卫星重力梯度测量的候选技术之一。本文针对未来更高分辨率全球重力场测量的科学需求,提出了一种适用于空间微重力环境下的原子干涉重力梯度测量方案,其梯度测量噪声可低至0.85mE/Hz1/2。文中对不同类型的卫星重力梯度测量方案进行了重力场反演精度的对比评估,仿真结果表明,相比于现有静电式卫星重力梯度测量,原子干涉型的卫星重力梯度测量有望将重力场的恢复阶数提升至252~290阶,对应的累积大地水准面误差7~8cm,累积重力异常误差3×10-5 m/s2。     

A model for adjustment of differential gravity measurements with simultaneous gravimeter calibration

 A mathematical model is proposed for adjustment of differential or relative gravity measurements, involving simultaneously instrumental readings, coefficients of the calibration function, and gravity values of selected base stations. Tests were performed with LaCoste and Romberg model G gravimeter measurements for a set of base stations located along a north–south line with 1750 mGal gravity range. This line was linked to nine control stations, where absolute gravity values had been determined by the free-fall method, with an accuracy better than 10 μGal. The model shows good consistence and stability. Results show the possibility of improving the calibration functions of gravimeters, as well as a better estimation of the gravity values, due to the flexibility admitted to the values of the calibration coefficients. Received: 15 November 1999 / Accepted: 31 October 2000     

Evaluation of the third- and fourth-generation GOCE Earth gravity field models with Australian terrestrial gravity data in spherical harmonics

In March 2013, the fourth generation of European Space Agency’s (ESA) global gravity field models, DIR4 (Bruinsma et al. in Proceedings of the ESA living planet symposium, 28 June–2 July, Bergen, ESA, Publication SP-686, 2010b) and TIM4 (Migliaccio et al. in Proceedings of the ESA living planet symposium, 28 June–2 July, Bergen, ESA, Publication SP-686, 2010), generated from the Gravity field and steady-state Ocean Circulation Explorer (GOCE) gravity observation satellite was released. We evaluate the models using an independent ground truth data set of gravity anomalies over Australia. Combined with Gravity Recovery and Climate Experiment (GRACE) satellite gravity, a new gravity model is obtained that is used to perform comparisons with GOCE models in spherical harmonics. Over Australia, the new gravity model proves to have significantly higher accuracy in the degrees below 120 as compared to EGM2008 and seems to be at least comparable to the accuracy of this model between degree 150 and degree 260. Comparisons in terms of residual quasi-geoid heights, gravity disturbances, and radial gravity gradients evaluated on the ellipsoid and at approximate GOCE mean satellite altitude ( $h=250$  km) show both fourth generation models to improve significantly w.r.t. their predecessors. Relatively, we find a root-mean-square improvement of 39 % for the DIR4 and 23 % for TIM4 over the respective third release models at a spatial scale of 100 km (degree 200). In terms of absolute errors, TIM4 is found to perform slightly better in the bands from degree 120 up to degree 160 and DIR4 is found to perform slightly better than TIM4 from degree 170 up to degree 250. Our analyses cannot confirm the DIR4 formal error of 1 cm geoid height (0.35 mGal in terms of gravity) at degree 200. The formal errors of TIM4, with 3.2 cm geoid height (0.9 mGal in terms of gravity) at degree 200, seem to be realistic. Due to combination with GRACE and SLR data, the DIR models, at satellite altitude, clearly show lower RMS values compared to TIM models in the long wavelength part of the spectrum (below degree and order 120). Our study shows different spectral sensitivity of different functionals at ground level and at GOCE satellite altitude and establishes the link among these findings and the Meissl scheme (Rummel and van Gelderen in Manusrcipta Geodaetica 20:379–385, 1995).     

Improvements in height datum transfer expected from the GOCE mission

 One of the aims of the Earth Explorer Gravity Field and Steady-State Ocean Circulation (GOCE) mission is to provide global and regional models of the Earth's gravity field and of the geoid with high spatial resolution and accuracy. Using the GOCE error model, simulation studies were performed in order to estimate the accuracy of datum transfer in different areas of the Earth. The results showed that with the GOCE error model, the standard deviation of the height anomaly differences is about one order of magnitude better than the corresponding value with the EGM96 error model. As an example, the accuracy of the vertical datum transfer from the tide gauge of Amsterdam to New York was estimated equal to 57 cm when the EGM96 error model was used, while in the case of GOCE error model this accuracy was increased to 6 cm. The geoid undulation difference between the two places is about 76.5 m. Scaling the GOCE errors to the local gravity variance, the estimated accuracy varied between 3 and 7 cm, depending on the scaling model. Received: 1 March 2000 / Accepted: 21 February 2001     

Precise ERS-2 orbit determination using SLR, PRARE, and RA observations

The second European Remote Sensing satellite ERS-2, launched in 1995, is supported by satellite laser ranging (SLR), the German precise range and range-rate equipment (PRARE) microwave tracking system, and radar altimetry (RA). SLR, PRARE and crossover RA observations (XO/RA) from May, June, and July 1996 have been analyzed using the GEOSAT software developed at Forsvarets forskningsinstitutt (FFI, The Norwegian defence research establishment). Orbits computed with the JGM3 gravity model and the more recent EGM96 gravity model have been compared with orbits from other analysis centers. Based on these comparisons in addition to the post-fit observation residuals, and results from internal orbital overlaps, we conclude that the radial component of the ERS-2 orbit can be determined with a precision of 6 cm. Received: 30 June 1997 / Accepted: 4 February 1998     

A synthetic Earth for use in geodesy

 A synthetic Earth and its gravity field that can be represented at different resolutions for testing and comparing existing and new methods used for global gravity-field determination are created. Both the boundary and boundary values of the gravity potential can be generated. The approach chosen also allows observables to be generated at aircraft flight height or at satellite altitude. The generation of the synthetic Earth shape (SES) and gravity-field quantities is based upon spherical harmonic expansions of the isostatically compensated equivalent rock topography and the EGM96 global geopotential model. Spherical harmonic models are developed for both the synthetic Earth topography (SET) and the synthetic Earth potential (SEP) up to degree and order 2160 corresponding to a 5′×5′ resolution. Various sets of SET, SES and SEP with boundary geometry and boundary values at different resolutions can be generated using low-pass filters applied to the expansions. The representation is achieved in point sets based upon refined triangulation of a octahedral geometry projected onto the chosen reference ellipsoid. The filter cut-offs relate to the sampling pattern in order to avoid aliasing effects. Examples of the SET and its gravity field are shown for a resolution with a Nyquist sampling rate of 8.27 degrees. Received: 6 August 1999 / Accepted: 26 April 2000     

Comparison of undulation difference accuracies using gravity anomalies and gravity disturbances

Errors are considered in the outer zone contribution to oceanic undulation differences as obtained from a set of potential coefficients complete to degree 180. It is assumed that the gravity data of the inner zone (a spherical cap), consisting of either gravity anomalies or gravity disturbances, has negligible error. This implies that error estimates of the total undulation difference are analyzed. If the potential coefficients are derived from a global field of 1°×1° mean anomalies accurate to ε_Δg=10 mgal, then for a cap radius of 10°, the undulation difference error (for separations between 100 km and 2000 km) ranges from 13 cm to 55 cm in the gravity anomaly case and from 6 cm to 36 cm in the gravity disturbance case. If ε_Δg is reduced to 1 mgal, these errors in both cases are less than 10 cm. In the absence of a spherical cap, both cases yield identical error estimates: about 68 cm if ε_Δg=1 mgal (for most separations) and ranging from 93 cm to 160 cm if ε_Δg=10 mgal. Introducing a perfect 30-degree reference field, the latter errors are reduced to about 110 cm for most separations.     

Ocean tidal parameters from starlette data

Starlette was launched in 1975 in order to study temporal variations in the Earth’s gravity field; in particular, tidal and Earth rotation effects. For the period April 1983 to April 1984 over12,700 normal points of laser ranging data to Starlette have been sub-divided into49 near consecutive 5–6 day arcs. Normal equations for each arc as obtained from a least-squares data reduction procedure, were solved for ocean tidal parameters along with other geodetic and geodynamic parameters. The tidal parameters are defined relative to Wahr’s body tides and Wahr’s nutation model and show fair agreement with other satellite derived results and those obtained from spherical harmonic decomposition of global ocean tidal models.     

Validation of GOCE gravity field models by means of orbit residuals and geoid comparisons

Three GOCE-based gravity field solutions have been computed by ESA’s high-level processing facility and were released to the user community. All models are accompanied by variance-covariance information resulting either from the least squares procedure or a Monte-Carlo approach. In order to obtain independent external quality parameters and to assess the current performance of these models, a set of independent tests based on satellite orbit determination and geoid comparisons is applied. Both test methods can be regarded as complementary because they either investigate the performance in the long wavelength spectral domain (orbit determination) or in the spatial domain (geoid comparisons). The test procedure was applied to the three GOCE gravity field solutions and to a number of selected pre-launch models for comparison. Orbit determination results suggest, that a pure GOCE gravity field model does not outperform the multi-year GRACE gravity field solutions. This was expected as GOCE is designed to improve the determination of the medium to high frequencies of the Earth gravity field (in the range of degree and order 50 to 200). Nevertheless, in case of an optimal combination of GOCE and GRACE data, orbit determination results should not deteriorate. So this validation procedure can also be used for testing the optimality of the approach adopted for producing combined GOCE and GRACE models. Results from geoid comparisons indicate that with the 2 months of GOCE data a significant improvement in the determination of the spherical harmonic spectrum of the global gravity field between degree 50 and 200 can be reached. Even though the ultimate mission goal has not yet been reached, especially due to the limited time span of used GOCE data (only 2 months), it was found that existing satellite-only gravity field models, which are based on 7 years of GRACE data, can already be enhanced in terms of spatial resolution. It is expected that with the accumulation of more GOCE data the gravity field model resolution and quality can be further enhanced, and the GOCE mission goal of 1–2 cm geoid accuracy with 100 km spatial resolution can be achieved.     

The ERS-2 altimetric bias and gravity field enhancement using dual crossovers between ERS and TOPEX/Poseidon

 Dual satellite crossovers (DXO) between the two European Remote Sensing satellites ERS-1 and ERS-2 and TOPEX/Poseidon are used to (1) refine the Earth's gravity field and (2) extend the study of the ERS-2 altimetric range stability to cover the first four years of its operation. The enhanced gravity field model, AGM-98, is validated by several methodologies and will be shown to provide, in particular, low geographically correlated orbital error for ERS-2. For the ERS-2 altimetric range study, TOPEX/Poseidon is first calibrated through comparison against in situ tide gauge data. A time series of the ERS-2 altimeter bias has been recovered along with other geophysical correction terms using tables for bias jumps in the range measurements at the single point target response (SPTR) events. On utilising the original version of the SPTR tables the overall bias drift is seen to be 2.6±1.0 mm/yr with an RMS of fit of 12.2 mm but with discontinuities at the centimetre level at the SPTR events. On utilising the recently released revised tables, SPTR2000, the drift is better defined at 2.4±0.6 mm/yr with the RMS of fit reduced to 3.7 mm. Investigations identify the sea-state bias as a source of error with corrections affecting the overall drift by close to 1.2 mm/yr. Received: 25 May 2000 / Accepted: 24 January 2001