On the geoid–quasigeoid separation in mountain areas

The separation between the reference surfaces for orthometric heights and normal heights—the geoid and the quasigeoid—is typically in the order of a few decimeters but can reach nearly 3 m in extreme cases. The knowledge of the geoid–quasigeoid separation with centimeter accuracy or better, is essential for the realization of national and international height reference frames, and for precision height determination in geodetic engineering. The largest contribution to the geoid–quasigeoid separation is due to the distribution of topographic masses. We develop a compact formulation for the rigorous treatment of topographic masses and apply it to determine the geoid–quasigeoid separation for two test areas in the Alps with very rough topography, using a very fine grid resolution of 100 m. The magnitude of the geoid–quasigeoid separation and its accuracy, its slopes, roughness, and correlation with height are analyzed. Results show that rigorous treatment of topographic masses leads to a rather small geoid–quasigeoid separation—only 30 cm at the highest summit—while results based on approximations are often larger by several decimeters. The accuracy of the topographic contribution to the geoid–quasigeoid separation is estimated to be 2–3 cm for areas with extreme topography. Analysis of roughness of the geoid–quasigeoid separation shows that a resolution of the modeling grid of 200 m or less is required to achieve these accuracies. Gravity and the vertical gravity gradient inside of topographic masses and the mean gravity along the plumbline are modeled which are important intermediate quantities for the determination of the geoid–quasigeoid separation. We conclude that a consistent determination of the geoid and quasigeoid height reference surfaces within an accuracy of few centimeters is feasible even for areas with extreme topography, and that the concepts of orthometric height and normal height can be consistently realized and used within this level of accuracy.     

On the indirect effect in the Stokes–Helmert method of geoid determination

The classical integral formula for determining the indirect effect in connection with the Stokes–Helmert method is related to a planar approximation of the sea level. A strict integral formula, as well as some approximations to it, are derived. It is concluded that the cap- size truncated integral formulas will suffer from the omission of some long-wavelength contributions, of the order of 50 cm in high mountains for the classical formula. This long-wavelength information can be represented by a set of spherical harmonic coefficients of the topography to, say, degree and order 360. Hence, for practical use, a combination of the classical formula and a set of spherical harmonics is recommended. Received: 10 March 1998 / Accepted: 16 November 1998     

On the equivalence of spherical splines with least-squares collocation and Stokes’s formula for regional geoid computation

This work is an investigation of three methods for regional geoid computation: Stokes’s formula, least-squares collocation (LSC), and spherical radial base functions (RBFs) using the spline kernel (SK). It is a first attempt to compare the three methods theoretically and numerically in a unified framework. While Stokes integration and LSC may be regarded as classic methods for regional geoid computation, RBFs may still be regarded as a modern approach. All methods are theoretically equal when applied globally, and we therefore expect them to give comparable results in regional applications. However, it has been shown by de Min (Bull Géod 69:223–232, 1995. doi: 10.1007/BF00806734) that the equivalence of Stokes’s formula and LSC does not hold in regional applications without modifying the cross-covariance function. In order to make all methods comparable in regional applications, the corresponding modification has been introduced also in the SK. Ultimately, we present numerical examples comparing Stokes’s formula, LSC, and SKs in a closed-loop environment using synthetic noise-free data, to verify their equivalence. All agree on the millimeter level.     

High-resolution regional geoid computation without applying Stokes’s formula: a case study of the Iranian geoid

A new theory for high-resolution regional geoid computation without applying Stokess formula is presented. Operationally, it uses various types of gravity functionals, namely data of type gravity potential (gravimetric leveling), vertical derivatives of the gravity potential (modulus of gravity intensity from gravimetric surveys), horizontal derivatives of the gravity potential (vertical deflections from astrogeodetic observations) or higher-order derivatives such as gravity gradients. Its algorithmic version can be described as follows: (1) Remove the effect of a very high degree/order potential reference field at the point of measurement (POM), in particular GPS positioned, either on the Earths surface or in its external space. (2) Remove the centrifugal potential and its higher-order derivatives at the POM. (3) Remove the gravitational field of topographic masses (terrain effect) in a zone of influence of radius r. A proper choice of such a radius of influence is 2r=4×10⁴ km/n, where n is the highest degree of the harmonic expansion. (cf. Nyquist frequency). This third remove step aims at generating a harmonic gravitational field outside a reference ellipsoid, which is an equipotential surface of a reference potential field. (4) The residual gravitational functionals are downward continued to the reference ellipsoid by means of the inverse solution of the ellipsoidal Dirichlet boundary-value problem based upon the ellipsoidal Abel–Poisson kernel. As a discretized integral equation of the first kind, downward continuation is Phillips–Tikhonov regularized by an optimal choice of the regularization factor. (5) Restore the effect of a very high degree/order potential reference field at the corresponding point to the POM on the reference ellipsoid. (6) Restore the centrifugal potential and its higher-order derivatives at the ellipsoidal corresponding point to the POM. (7) Restore the gravitational field of topographic masses ( terrain effect) at the ellipsoidal corresponding point to the POM. (8) Convert the gravitational potential on the reference ellipsoid to geoidal undulations by means of the ellipsoidal Bruns formula. A large-scale application of the new concept of geoid computation is made for the Iran geoid. According to the numerical investigations based on the applied methodology, a new geoid solution for Iran with an accuracy of a few centimeters is achieved.Acknowledgments. The project of high-resolution geoid computation of Iran has been support by National Cartographic Center (NCC) of Iran. The University of Tehran, via grant number 621/3/602, supported the computation of a global geoid solution for Iran. Their support is gratefully acknowledged. A. Ardalan would like to thank Mr. Y. Hatam, and Mr. K. Ghazavi from NCC and Mr. M. Sharifi, Mr. A. Safari, and Mr. M. Motagh from the University of Tehran for their support in data gathering and computations. The authors would like to thank the comments and corrections made by the four reviewers and the editor of the paper, Professor Will Featherstone. Their comments helped us to correct the mistakes and improve the paper.     

Topographic effects by the Stokes–Helmert method of geoid and quasi-geoid determinations

The topographic potential and the direct topographic effect on the geoid are presented as surface integrals, and the direct gravity effect is derived as a rigorous surface integral on the unit sphere. By Taylor-expanding the integrals at sea level with respect to topographic elevation (H) the power series of the effects is derived to arbitrary orders. This study is primarily limited to terms of order H ². The limitations of the various effects in the frequently used planar approximations are demonstrated. In contrast, it is shown that the spherical approximation to power H ² leads to a combined topographic effect on the geoid (direct plus indirect effect) proportional to H˜² (where terms of degrees 0 and 1 are missing) of the order of several metres, while the combined topographic effect on the height anomaly vanishes, implying that current frequent efforts to determine the direct effect to this order are not needed. The last result is in total agreement with Bjerhammar's method in physical geodesy. It is shown that the most frequently applied remove–restore technique of topographic masses in the application of Stokes' formula suffers from significant errors both in the terrain correction C (representing the sum of the direct topographic effect on gravity anomaly and the effect of continuing the anomaly to sea level) and in the term t (mainly representing the indirect effect on the geoidal or quasi-geoidal height). Received: 18 August 1998 / Accepted: 4 October 1999     

Regional geoid computation by least squares modified Hotine’s formula with additive corrections

Geoid and quasigeoid modelling from gravity anomalies by the method of least squares modification of Stokes’s formula with additive corrections is adapted for the usage with gravity disturbances and Hotine’s formula. The biased, unbiased and optimum versions of least squares modification are considered. Equations are presented for the four additive corrections that account for the combined (direct plus indirect) effect of downward continuation (DWC), topographic, atmospheric and ellipsoidal corrections in geoid or quasigeoid modelling. The geoid or quasigeoid modelling scheme by the least squares modified Hotine formula is numerically verified, analysed and compared to the Stokes counterpart in a heterogeneous study area. The resulting geoid models and the additive corrections computed both for use with Stokes’s or Hotine’s formula differ most in high topography areas. Over the study area (reaching almost 2 km in altitude), the approximate geoid models (before the additive corrections) differ by 7 mm on average with a 3 mm standard deviation (SD) and a maximum of 1.3 cm. The additive corrections, out of which only the DWC correction has a numerically significant difference, improve the agreement between respective geoid or quasigeoid models to an average difference of 5 mm with a 1 mm SD and a maximum of 8 mm.     

The direct topographical correction in gravimetric geoid determination by the Stokes–Helmert method

 The direct topographical correction is composed of both local effects and long-wavelength contributions. This implies that the classical integral formula for determining the direct effect may have some numerical problems in representing these different signals. On the other hand, a representation by a set of harmonic coefficients of the topography to, say, degree and order 360 will omit significant short-wavelength signals. A new formula is derived by combining the classical formula and a set of spherical harmonics. Finally, the results of this solution are compared with the Moritz topographical correction in a test area. Received: 27 July 1998 / Accepted: 29 March 2000     

On the Non-linear Motion of IGS Station

With the daily SINEX files of the IGS, the time series of IGS stations are obtained using an independently developed software under generalized network adjustment models with coordinate patterns. From ...     

The spheroidal fixed–free two-boundary-value problem for geoid determination (the spheroidal Bruns' transform)

The target of the spheroidal Gauss–Listing geoid determination is presented as a solution of the spheroidal fixed–free two-boundary value problem based on a spheroidal Bruns' transformation (“spheroidal Bruns' formula”). The nonlinear spheroidal Bruns' transform (nonlinear spheroidal Bruns' formula), the spheroidal fixed part and the spheroidal free part of the two-boundary value problem are derived. Four different spheroidal gravity models are treated, in particular to determine whether they pass the test to fit to the postulate of a level ellipsoidal gravity field, namely of Somigliana–Pizzetti type. Received: 4 May 1999 / Accepted: 21 May 1999     

On non-linear low–low SST observation equations for the determination of the geopotential based on an analytical solution

In satellite data analysis, one big advantage of analytical orbit integration, which cannot be overestimated, is missed in the numerical integration approach: spectral analysis or the lumped coefficient concept may be used not only to design efficient algorithms but overall for much better insight into the force-field determination problem. The lumped coefficient concept, considered from a practical point of view, consists of the separation of the observation equation matrix A=BT into the product of two matrices. The matrix T is a very sparse matrix separating into small block-diagonal matrices connecting the harmonic coefficients with the lumped coefficients. The lumped coefficients are nothing other than the amplitudes of trigonometric functions depending on three angular orbital variables; therefore, the matrix N=B ^T B will become for a sufficient length of a data set a diagonal dominant matrix, in the case of an unlimited data string length a strictly diagonal one. Using an analytical solution of high order, the non-linear observation equations for low–low SST range data can be transformed into a form to allow the application of the lumped concept. They are presented here for a second-order solution together with an outline of how to proceed with data analysis in the spectral domain in such a case. The dynamic model presented here provides not only a practical algorithm for the parameter determination but also a simple method for an investigation of some fundamental questions, such as the determination of the range of the subset of geopotential coefficients which can be properly determined by means of SST techniques or the definition of an optimal orbital configuration for particular SST missions. Numerical results have already been obtained and will be published elsewhere. Received: 15 January 1999 / Accepted: 30 November 1999     

Generalization of the Lambert–Lagrange projection

The Lagrange projection represents conformally the terrestrial globe within a circle. This is achieved by compressing the latitude and longitude and by applying the new coordinates into the equatorial stereographic projection. The same concept can be generalized to any conformal projection, although the application of this technique to other analytical functions is less known. In this work, the general Lambert–Lagrange projection formula is proposed and the application of the modified coordinates is discussed on projections: stereographic, conformal conic and Gauss–Schreiber. In general, the results are merely a curiosity, except for the case of Gauss–Schreiber, where the use of coordinates with altered scale can be applied in the optimization of conformal projections.     

Unification of New Zealand’s local vertical datums: iterative gravimetric quasigeoid computations

New Zealand uses 13 separate local vertical datums (LVDs) based on geodetic levelling from 12 different tide-gauges. We describe their unification using a regional gravimetric quasigeoid model and GPS-levelling data on each LVD. A novel application of iterative quasigeoid computation is used, where the LVD offsets computed from earlier models are used to apply additional gravity reductions from each LVD to that model. The solution converges after only three iterations yielding LVD offsets ranging from 0.24 to 0.58 m with an average standard deviation of ±0.08 m. The so-computed LVD offsets agree, within expected data errors, with geodetically levelled height differences at common benchmarks between adjacent LVDs. This shows that iterated quasigeoid models have a role in vertical datum unification.     

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     

Determination of the boundary values for the Stokes–Helmert problem

 The definition of the mean Helmert anomaly is reviewed and the theoretically correct procedure for computing this quantity on the Earth's surface and on the Helmert co-geoid is suggested. This includes a discussion of the role of the direct topographical and atmospherical effects, primary and secondary indirect topographical and atmospherical effects, ellipsoidal corrections to the gravity anomaly, its downward continuation and other effects. For the rigorous derivations it was found necessary to treat the gravity anomaly systematically as a point function, defined by means of the fundamental gravimetric equation. It is this treatment that allows one to formulate the corrections necessary for computing the `one-centimetre geoid'. Compared to the standard treatment, it is shown that a `correction for the quasigeoid-to-geoid separation', amounting to about 3 cm for our area of interest, has to be considered. It is also shown that the `secondary indirect effect' has to be evaluated at the topography rather than at the geoid level. This results in another difference of the order of several centimetres in the area of interest. An approach is then proposed for determining the mean Helmert anomalies from gravity data observed on the Earth's surface. This approach is based on the widely-held belief that complete Bouguer anomalies are generally fairly smooth and thus particularly useful for interpolation, approximation and averaging. Numerical results from the Canadian Rocky Mountains for all the corrections as well as the downward continuation are shown. Received: 9 March 1998 / Accepted: 16 November 1998     

On the influence of the ground track on the gravity field recovery from high–low satellite-to-satellite tracking missions: CHAMP monthly gravity field recovery using the energy balance approach revisited

In this paper, the influence of the ground track coverage on the quality of a monthly gravity field solution is investigated for the scenario of a high–low satellite- to-satellite tracking mission. Data from the CHAllenging Minisatellite Payload (champ) mission collected in the period April 2002 to February 2004 has been used to recover the gravity field to degree and order 70 on a monthly basis. The quality is primarily restricted by the accuracy of the instruments. Besides, champ passed through a 31/2 repeat mode three times during the period of interest resulting in an insufficient spatial sampling and a degraded solution. Contrary to the rule of thumb by Colombo (The global mapping of gravity with two satellites, Publications on Geodesy, vol 7(3), Netherlands Geodetic Commission, The Netherlands, 263 pp, 1984), see also Wagner (J Geod 80(2): 94–103, 2006), we found that the monthly solutions themselves could be recovered to about degree 30, not 15. In order to improve the monthly gravity solutions, two strategies have been developed: the restriction to a low degree, and the densification of the sampling by the introduction of additional sensitive measurements from contemporaneous satellite missions. The latter method is tested by combining the champ measurements with data from the Gravity Recovery And Climate Experiment (grace). Note that the two grace satellites are considered independent here, i.e. no use is made of the K-band ranging data. This way, we are able to almost entirely remove the influence of the ground track leaving the accuracy of the instruments as the primary restriction on the quality of a monthly solution. These findings are especially interesting for the upcoming swarm-mission since it will consist of a similar configuration as the combined champ and (grace) missions.     

On Definition of Geograghic Variable Scaling

IntroductionThe map is a basic form of geographic informationvisualization[1]. To provide space attributes or geo-graphic orders is the basic function of a map. Incartography, according to the different measure ofphenomenal quantitative attribute, four fo…