Local geoid determination combining gravity disturbances and GPS/levelling: a case study in the Lake Nasser area, Aswan, Egypt

 The use of GPS for height control in an area with existing levelling data requires the determination of a local geoid and the bias between the local levelling datum and the one implicitly defined when computing the local geoid. If only scarse gravity data are available, the heights of new data may be collected rapidly by determining the ellipsoidal height by GPS and not using orthometric heights. Hence the geoid determination has to be based on gravity disturbances contingently combined with gravity anomalies. Furthermore, existing GPS/levelling data may also be used in the geoid determination if a suitable general gravity field modelling method (such as least-squares collocation, LSC) is applied. A comparison has been made in the Aswan Dam area between geoids determined using fast Fourier transform (FFT) with gravity disturbances exclusively and LSC using only the gravity disturbances and the disturbances combined with GPS/levelling data. The EGM96 spherical harmonic model was in all cases used in a remove–restore mode. A total of 198 gravity disturbances spaced approximately 3 km apart were used, as well as 35 GPS/levelling points in the vicinity and on the Aswan Dam. No data on the Nasser Lake were available. This gave difficulties when using FFT, which requires the use of gridded data. When using exclusively the gravity disturbances, the agreement between the GPS/levelling data were 0.71 ± 0.17 m for FFT and 0.63 ± 0.15 for LSC. When combining gravity disturbances and GPS/levelling, the LSC error estimate was ±0.10 m. In the latter case two bias parameters had to be introduced to account for a possible levelling datum difference between the levelling on the dam and that on the adjacent roads. Received: 14 August 2000 / Accepted: 28 February 2001     

On the combined adjustment of a gravimetrically determined geoid and GPS levelling stations

Summary A new basic geodetic network using the GPS technique is now being set up in France. There will be altogether 1000 benchmarks connected to the French levelling network. Obviously, the GPS levelling points are not dense enough to produce a national levelling reference surface. A gravimetrically determined geoid has therefore been proposed to be used for the interpolation between the GPS levelling points. However, because of long-wavelength errors, we consider that a gravimetric geoid does not have sufficient accuracy. A regression by fitting the gravimetrically determined geoid to the GPS levelling points is generally proposed. Unfortunately, this country-wide geoid fitting work cannot eliminate local deformations in the geoid, which happen in areas where there are errors or shortages of gravity or DTM data. This paper proposes and discusses a combined adjustment method. The principle is to divide up the geoid into small pieces and then to adjust them to the GPS levelling points locally with constraint conditions for the common points of the adjacent pieces. In order to benefit from the advantages of the high resolution and high relative accuracy of the gravimetric geoid, as well as the high absolute accuracy of the GPS levelling points, we establish respectively a relative observation equation for the difference of the gravimetric geoid undulation and an absolute observation equation for the GPS levelling points. Finally, we adjust the observation equations as a whole. Several global and local systematic errors are also taken into account and some special cases, such as adjustment in groups and blunder detection, are also discussed.     

A new hybrid geoid model for Japan, GSIGEO2000

 The latest gravimetric geoid model for Japan, JGEOID2000, was successfully combined with the nationwide net of GPS at benchmarks, yielding a new hybrid geoid model for Japan, GSIGEO2000. The least-squares collocation (LSC) method was applied as an interpolation for fitting JGEOID2000 to the GPS/leveling geoid undulations. The GPS/leveling geoid undulation data were reanalyzed in advance, in terms of three-dimensional positions from GPS and orthometric heights from leveling. The new hybrid geoid model is, therefore, compatible with the new Japanese geodetic reference frame. GSIGEO2000 was evaluated internally and independently and the precision was estimated at 4 cm throughout nearly the whole region. Received: 15 October 2001 / Accepted: 27 March 2002 Acknowledgments. Messrs. Toshio Kunimi and Tadashi Saito at the Third Geodetic Division of the Geographical Survey Institute (GSI) mainly carried out the computations of most of the updated leveled heights. With regard to the reanalysis of GPS data, the discussions with Messrs. Yuki Hatanaka and Shoichi Matsumura of GSI were of great help in building the analysis strategy. Messrs. Kazuyuki Tanaka and Hiromi Shigematsu collaborated in the preparatory stages of GPS data computation. The authors' thanks are extended to these colleagues. Some plots were made by GMT software (Wessel and Smith 1991). Correspondence to: Y. Kuroishi     

配置法在局部大地水准面精化中的实例研究

从最小二乘配置方法的基本原理出发,以我国某地区范围内1km分辨率的大地水准面高模型数据为例,根据实用公式计算了试验区大地水准面高的协方差值后,采用多项式函数模型和高斯函数模型分别拟合了该地区大地水准面高的局部协方差函数,并对试验区内18个检核点做了推估计算。根据推估值(Nfit)与实测值(NGPSL)的比较分析表明,虽然多项式协方差函数模型略优于高斯协方差函数模型,但它们都能以厘米级的精度拟合局部大地水准面,这表明了配置法用于精化厘米级大地水准面的有效性。     

How to handle topography in practical geoid determination: three examples

 Three different methods of handling topography in geoid determination were investigated. The first two methods employ the residual terrain model (RTM) remove–restore technique, yielding the quasigeoid, whereas the third method uses the classical Helmert condensation method, yielding the geoid. All three methods were used with the geopotential model Earth Gravity Model (1996) (EGM96) as a reference, and the results were compared to precise global positioning system (GPS) levelling networks in Scandinavia. An investigation of the Helmert method, focusing on the different types of indirect effects and their effects on the geoid, was also carried out. The three different methods used produce almost identical results at the 5-cm level, when compared to the GPS levelling networks. However, small systematic differences existed. Received: 18 March 1999 / Accepted: 21 March 2000     

On the adjustment of combined GPS/levelling/geoid networks

A detailed treatment of adjustment problems in combined global positioning system (GPS)/levelling/geoid networks is given. The two main types of `unknowns' in this kind of multi-data 1D networks are usually the gravimetric geoid accuracy and a 2D spatial field that describes all the datum/systematic distortions among the available height data sets. An accurate knowledge of the latter becomes especially important when we consider employing GPS techniques for levelling purposes with respect to a local vertical datum. Two modelling alternatives for the correction field are presented, namely a pure deterministic parametric model, and a hybrid deterministic and stochastic model. The concept of variance component estimation is also proposed as an important statistical tool for assessing the actual gravimetric geoid noise level and/or testing a priori determined geoid error models. Finally, conclusions are drawn and recommendations for further study are suggested. Received: 9 September 1998 / Accepted: 8 June 1999     

Minimization and estimation of geoid undulation errors

The objective of this paper is to minimize the geoid undulation errors by focusing on the contribution of the global geopotential model and regional gravity anomalies, and to estimate the accuracy of the predicted gravimetric geoid.The geopotential model's contribution is improved by (a) tailoring it using the regional gravity anomalies and (b) introducing a weighting function to the geopotential coefficients. The tailoring and the weighting function reduced the difference (1) between the geopotential model and the GPS/levelling-derived geoid undulations in British Columbia by about 55% and more than 10%, respectively.Geoid undulations computed in an area of 40° by 120° by Stokes' integral with different kernel functions are analyzed. The use of the approximated kernels results in about 25 cm () and 190 cm (maximum) geoid errors. As compared with the geoid derived by GPS/levelling, the gravimetric geoid gives relative differences of about 0.3 to 1.4 ppm in flat areas, and 1 to 2.5 ppm in mountainous areas for distances of 30 to 200 km, while the absolute difference (1) is about 5 cm and 20 cm, respectively.A optimal Wiener filter is introduced for filtering of the gravity anomaly noise, and the performance is investigated by numerical examples. The internal accuracy of the gravimetric geoid is studied by propagating the errors of the gravity anomalies and the geopotential coefficients into the geoid undulations. Numerical computations indicate that the propagated geoid errors can reasonably reflect the differences between the gravimetric and GPS/levelling-derived geoid undulations in flat areas, such as Alberta, and is over optimistic in the Rocky Mountains of British Columbia.Paper presented at the IAG General Meeting, Beijing, China, August 8–13, 1993.     

区域似大地水准面精化的方法探讨与精度分析

阐述与GPS/水准高程拟合相关的基本理论及区域似大地水准面精化的原理与方法,并结合区域的具体实例,重点比较二次多项式曲面拟合法、三角剖分双线性内插法、加权平均推估法的精度情况,从中得到结论。实验表明,在没有足够重力数据的支持下,三角剖分双线性内插法获得比较理想的精度,利用该方法精化似大地水准面获得的高程可以代替四等水准。     

青藏高原中部似大地水准面精化的应用研究

介绍了精化似大地水准面的常用方法。采用青海和西藏三级GPS大地控制网测量数据,利用移动曲面法、改进的移动曲面法、二次曲面拟合法和改进的二次曲面拟合法对青藏高原中部似大地水准面进行了拟合计算。结果表明,利用改进的移动曲面法和改进的二次曲面法能有效地提高本地区似大地水准面的精化效果,精度优于0.5 m。     

The strategy for the determination of a new local geoid in China

The strategy for the determination of a new local geoid in China is: (1) This task should be fulfilled in a short time as possible; (2) The new local geoid can be with different resolution and accuracy, which depend on the needs of the social and economical development in different regions; (3) The newly established height anomaly control network, i.e. a gridded GPS levelling network should combine as possible with the existing astro-gravimetric levelling network; (4) Global geopotential model, DTM and surface gravity data are used to densify the resolution of the local geoid with remove-restore technique. The accuracy and resolution possibly obtained nowadays within the capability of China are discussed. The required accuracies for GPS positioning, levelling and coordinate transformation are suggested.     

Undulation and anomaly estimation using Geos-3 altimeter data without precise satellite orbits

The paper describes results obtained from the processing of 53 Geos-3 arcs of altimeter data obtained during the first weeks after the launch of the satellite in April, 1975. The measurement from the satellite to the ocean surface was used to obtain an approximate geoid undulation which was contaminated by long wavelength errors caused primarily by altimeter bias and orbit error. This long wavelength error was reduced by fitting with a low degree polynomial the raw undulation data to the undulations implied by the GEM 7 potential coefficients, in an adjustment process that included conditions on tracks that cross. The root mean square crossover discrepancy before this adjustment was ±12.4 meters while after the adjustment it was ±0.9 m. These adjusted undulations were used to construct a geoid map in the Geos-3 calibration area using a least squares filter to remove remaining noise in the undulations. Comparing these undulations to ones computed from potential coefficients and terrestrial gravity data indicates a mean difference of 0.25 m and a root mean square difference of ±1.92 m. The adjusted undulations were also used to estimate several 5^o, 2^o, and 1^o anomalies using the method of least squares collocation. The resulting predictions agreed well with known values although the 1^o x 1^o anomalies could not be considered as reliably determined.     

A strategy for determining the regional geoid by combining limited ground data with satellite-based global geopotential and topographical models: a case study of Iran

The computation of regional gravimetric geoid models with reasonable accuracy, in developing countries, with sparse data is a difficult task that needs great care. Here we investigate the procedure for gathering, evaluating and combining different data for the determination of a gravimetric geoid model for Iran, where limited ground gravity data are available. Heterogeneous data, including gravity anomalies, the high-resolution Shuttle Radar Topography Mission global digital terrain model and different global geopotential models including recently published Gravity Recovery and Climate Experiment models, are combined through least-squares modification of the Stokes formula. The new gravimetric geoid model, IRG04, agrees considerably better with GPS/levelling than any of the other recent local geoid model in the area. Its RMS fit with GPS/levelling is 0.27 m and 3.8 ppm in the absolute and relative view, respectively. The relative accuracy of IRG04 is four times better than the most recently published global and regional geoid models available in this area. This progress shows the practical potential of the method of least-squares modification of Stokes’s formula in combination with heterogeneous data for regional geoid determination     

提高GPS高程拟合模型精度的有效方法

在调兵山地区,选择76大地高和正常高都已知的点作为原始数据,采用多项式拟合加权平均法建立该地区大地水准面,包括泰勒二次多项式加权平均法和双二次多项式加权平均法。通过实例发现,泰勒二多项式曲面加权平均法和双二次多项式加权平均法对高程异常的短波项有很好的改进效果,拟合残差值比较大的点采用泰勒二次多项式曲面加权平均法和双二次多项式加权平均法残差值变小,大大提高了拟合模型的精度,改善了拟合模型。     

Models for Fitting Gravimetric Geoids and GPS Results

The aim of this investigation is to study how to use a gravimetric(quasi) geoid for levelling by GPS data in an optimal way.The advent of precise geodetic GPS has made the use of a technique possible,which might be called GPS- gravimetric geoid determination.In this approach,GPS heights above the reference ellipsoid are determined for points whose levelled (orthometric) height H is above sea level people have already surveyed;for these points,we thus have the values of the geoid undulation N.These values are then used to constrain the geoid undulations N‘ obtained from the gravimetric solution.     

Fitting gravimetric geoid models to vertical deflections

Regional gravimetric geoid and quasigeoid models are now commonly fitted to GPS-levelling data, which simultaneously absorbs levelling, GPS and quasi/geoid errors due to their inseparability. We propose that independent vertical deflections are used instead, which are not affected by this inseparability problem. The formulation is set out for geoid slopes and changes in slopes. Application to 1,080 astrogeodetic deflections over Australia for the AUSGeoid98 model shows that it is feasible, but the poor quality of the historical astrogeodetic deflections led to some unrealistic values.     

Effect of the SRTM global DEM on the determination of a high-resolution geoid model: a case study in Iran

Any errors in digital elevation models (DEMs) will introduce errors directly in gravity anomalies and geoid models when used in interpolating Bouguer gravity anomalies. Errors are also propagated into the geoid model by the topographic and downward continuation (DWC) corrections in the application of Stokes’s formula. The effects of these errors are assessed by the evaluation of the absolute accuracy of nine independent DEMs for the Iran region. It is shown that the improvement in using the high-resolution Shuttle Radar Topography Mission (SRTM) data versus previously available DEMs in gridding of gravity anomalies, terrain corrections and DWC effects for the geoid model are significant. Based on the Iranian GPS/levelling network data, we estimate the absolute vertical accuracy of the SRTM in Iran to be 6.5 m, which is much better than the estimated global accuracy of the SRTM (say 16 m). Hence, this DEM has a comparable accuracy to a current photogrammetric high-resolution DEM of Iran under development. We also found very large differences between the GLOBE and SRTM models on the range of −750 to 550 m. This difference causes an error in the range of −160 to 140 mGal in interpolating surface gravity anomalies and −60 to 60 mGal in simple Bouguer anomaly correction terms. In the view of geoid heights, we found large differences between the use of GLOBE and SRTM DEMs, in the range of −1.1 to 1 m for the study area. The terrain correction of the geoid model at selected GPS/levelling points only differs by 3 cm for these two DEMs.     

基于Matlab的GPS高程异常拟合

由于似大地水准面是一个不规则的曲面,它无法用一个精准的曲面来模拟,这就使得GPS只能提供给我们高精度的大地高,而不是工程中需要的正常高,严重影响了GPS3维定位的应用发展,使其提供维坐标的优越性未能得到充分发挥。在借鉴国内外研究成果的基础上,针对我国国土面积广,地形尤为复杂,似大地水准面起伏比较大的特点,结合工程实例,对大面积区域内GPS高程异常拟合作了相应的分析和研究。     

Robust estimation and optimal selection of polynomial parameters for the interpolation of GPS geoid heights

The polynomial interpolation of least squares and interpolation moving least squares based on control stations with known GPS (global positioning system) ellipsoidal heights and levelling orthometric heights are the most often used methods for the interpolation of the geoid heights. But in their applications there occur two problems: one lies in selecting the suitable polynomial parameters; the other in reducing the influences of some possibly abnormal data points. To solve both of the problems, without emphasizing a sound theoretical basis, a heuristic solution with the help of robust estimation technique and optimization criteria for the regression equation is presented. Through two actual numerical examples it is shown that the new solution concept is efficient and can be realized easily on computers. Received: 23 May 1996 / Accepted: 27 March 1997     

Geoid undulation modelling and interpretation at Ladak, NW Himalaya using GPS and levelling data

Fast and accurate relative positioning for baselines less than 20 km in length is possible using dual-frequency Global Positioning System (GPS) receivers. By measuring orthometric heights of a few GPS stations by differential levelling techniques, the geoid undulation can be modelled, which enables GPS to be used for orthometric height determination in a much faster and more economical way than terrestrial methods. The geoid undulation anomaly can be very useful for studying tectonic structure. GPS, levelling and gravity measurements were carried out along a 200-km-long highly undulating profile, at an average elevation of 4000 m, in the Ladak region of NW Himalaya, India. The geoid undulation and gravity anomaly were measured at 28 common GPS-levelling and 67 GPS-gravity stations. A regional geoid low of nearly −4 m coincident with a steep negative gravity gradient is compatible with very recent findings from other geophysical studies of a low-velocity layer 20–30 km thick to the north of the India–Tibet plate boundary, within the Tibetan plate. Topographic, gravity and geoid data possibly indicate that the actual plate boundary is situated further north of what is geologically known as the Indus Tsangpo Suture Zone, the traditionally supposed location of the plate boundary. Comparison of the measured geoid with that computed from OSU91 and EGM96 gravity models indicates that GPS alone can be used for orthometric height determination over the Higher Himalaya with 1–2 m accuracy. Received: 10 April 1997 / Accepted: 9 October 1998     

The northern European geoid: a case study on long-wavelength geoid errors

 The long-wavelength geoid errors on large-scale geoid solutions, and the use of modified kernels to mitigate these effects, are studied. The geoid around the Nordic area, from Greenland to the Ural mountains, is considered. The effect of including additional gravity data around the Nordic/Baltic land area, originating from both marine, satellite and ground-based measurements, is studied. It is found that additional data appear to increase the noise level in computations, indicating the presence of systematic errors. Therefore, the Wong–Gore modification to the Stokes kernel is applied. This method of removing lower-order terms in the Stokes kernel appears to improve the geoid. The best fit to the global positioning system (GPS) leveling points is obtained with a degree of modification of approximately 30. In addition to the study of modification errors, the results of different methods of combining satellite altimetry gravity and other gravimetry are presented. They all gave comparable results, at the 6-cm level, when evaluated for the Nordic GPS networks. One dimensional (1-D) and 2-D fast Fourier transform (FFT) methods are also compared. It is shown that even though methods differ by up to 6 cm, the fit to the GPS is essentially the same. A surprising conclusion is that the addition of more data does not always produce a better geoid, illustrating the danger of systematic errors in data. Received: 4 July 2001 / Accepted: 21 February 2002