On the ellipsoidal correction to the spherical Stokes solution of the gravimetric geoid

The solutions of four ellipsoidal approximations for the gravimetric geoid are reviewed: those of Molodenskii et al., Moritz, Martinec and Grafarend, and Fei and Sideris. The numerical results from synthetic tests indicate that Martinec and Grafarends solution is the most accurate, while the other three solutions contain an approximation error which is characterized by the first-degree surface spherical harmonic. Furthermore, the first 20 degrees of the geopotential harmonic series contribute approximately 90% of the ellipsoidal correction. The determination of a geoid model from the generalized Stokes scheme can accurately account for the ellipsoidal effect to overcome the first-degree surface spherical harmonic error regardless of the solution used.     

The nullfield method for the ellipsoidal Stokes problem

The ellipsoidal Stokes problem is one of the basic boundary-value problems for the Laplace equation which arises in physical geodesy. Up to now, geodecists have treated this and related problems with high-order series expansions of spherical and spheroidal (ellipsoidal) harmonics. In view of increasing computational power and modern numerical techniques, boundary element methods have become more and more popular in the last decade. This article demonstrates and investigates the nullfield method for a class of Robin boundary-value problems. The ellipsoidal Stokes problem belongs to this class. An integral equation formulation is achieved, and existence and uniqueness conditions are attained in view of the Fredholm alternative. Explicit expressions for the eigenvalues and eigenfunctions for the boundary integral operator are provided. Received: 22 October 1996 / Accepted: 4 August 1997     

On the combination of gravity anomalies and gravity disturbances for geoid determination in Western Australia

The geoid gradient over the Darling Fault in Western Australia is extremely high, rising by as much as 38 cm over only 2 km. This poses problems for gravimetric-only geoid models of the area, whose frequency content is limited by the spatial distribution of the gravity data. The gravimetric-only version of AUSGeoid98, for instance, is only able to resolve 46% of the gradient across the fault. Hence, the ability of GPS surveys to obtain accurate orthometric heights is reduced. It is described how further gravity data were collected over the Darling Fault, augmenting the existing gravity observations at key locations so as to obtain a more representative geoid gradient. As many of the gravity observations were collected at stations with a well-known GRS80 ellipsoidal height, the opportunity arose to compute a geoid model via both the Stokes and the Hotine approaches. A scheme was devised to convert free-air anomaly data to gravity disturbances using existing geoid models, followed by a Hotine integration to geoid heights. Interestingly, these results depended very weakly upon the choice of input geoid model. The extra gravity data did indeed improve the fit of the computed geoid to local GPS/Australian Height Datum (AHD) observations by 58% over the gravimetric-only AUSGeoid98. While the conventional Stokesian approach to geoid determination proved to be slightly better than the Hotine method, the latter still improved upon the gravimetric-only AUSGeoid98 solution, supporting the viability of conducting gravity surveys with GPS control for the purposes of geoid determination. AcknowledgementsThe author would like to thank Will Featherstone, Ron Gower, Ron Hackney, Linda Morgan, Geoscience Australia, Scripps Oceanographic Institute and the three anonymous reviewers of this paper. This research was funded by the Australian Research Council.     

On the ellipsoidal corrections to gravity anomalies computed using the inverse Stokes integral

The formulas of the ellipsoidal corrections to the gravity anomalies computed using the inverse Stokes integral are derived. The corrections are given in the integral formulas and expanded in the spherical harmonics series. If a coefficient model such as the OSU91A is given, the corrections can be easily computed. Received: 19 August 1996 / Accepted: 28 September 1998     

Ellipsoidal geoid computation

Modern geoid computation uses a global gravity model, such as EGM96, as a third component in a remove–restore process. The classical approach uses only two: the reference ellipsoid and a geometrical model representing the topography. The rationale for all three components is reviewed, drawing attention to the much smaller precision now needed when transforming residual gravity anomalies. It is shown that all ellipsoidal effects needed for geoid computation with millimetric accuracy are automatically included provided that the free air anomaly and geoid are calculated correctly from the global model. Both must be consistent with an ellipsoidal Earth and with the treatment of observed gravity data. Further ellipsoidal corrections are then negligible. Precise formulae are developed for the geoid height and the free air anomaly using a global gravity model, given as spherical harmonic coefficients. Although only linear in the anomalous potential, these formulae are otherwise exact for an ellipsoidal reference Earth—they involve closed analytical functions of the eccentricity (and the Earths spin rate), rather than a truncated power series in e². They are evaluated using EGM96 and give ellipsoidal corrections to the conventional free air anomaly ranging from –0.84 to +1.14 mGal, both extremes occurring in Tibet. The geoid error corresponding to these differences is dominated by longer wavelengths, so extrema occur elsewhere, rising to +766 mm south of India and falling to –594 mm over New Guinea. At short wavelengths, the difference between ellipsoidal corrections based only on EGM96 and those derived from detailed local gravity data for the North Sea geoid GEONZ97 has a standard deviation of only 3.3 mm. However, the long-wavelength components missed by the local computation reach 300 mm and have a significant slope. In Australia, for example, such a slope would amount to a 600-mm rise from Perth to Cairns.     

The ellipsoidal corrections to the topographic geoid effects

The topographic effects by Stokes formula are typically considered for a spherical approximation of sea level. For more precise determination of the geoid, sea level is better approximated by an ellipsoid, which justifies the consideration of the ellipsoidal corrections of topographic effects for improved geoid solutions. The aim of this study is to estimate the ellipsoidal effects of the combined topographic correction (direct plus indirect topographic effects) and the downward continuation effect. It is concluded that the ellipsoidal correction to the combined topographic effect on the geoid height is far less than 1 mm. On the contrary, the ellipsoidal correction to the effect of downward continuation of gravity anomaly to sea level may be significant at the 1-cm level in mountainous regions. Nevertheless, if Stokes formula is modified and the integration of gravity anomalies is limited to a cap of a few degrees radius around the computation point, nor this effect is likely to be significant.AcknowledgementsThe author is grateful for constructive remarks by J Ågren and the three reviewers.     

A new strategy for the atmospheric gravity effect in gravimetric geoid determination

Prior to Stokes integration, the gravitational effect of atmospheric masses must be removed from the gravity anomaly g. One theory for the atmospheric gravity effect on the geoid is the well-known International Association of Geodesy approach in connection with Stokes integral formula. Another strategy is the use of a spherical harmonic representation of the topography, i.e. the use of a global topography computed from a set of spherical harmonics. The latter strategy is improved to account for local information. A new formula is derived by combining the local contribution of the atmospheric effect computed from a detailed digital terrain model and the global contribution computed from a spherical harmonic model of the topography. The new formula is tested over Iran and the results are compared with corresponding results from the old formula which only uses the global information. The results show significant differences. The differences between the two formulas reach 17 cm in a test area in Iran.     

Transforming ellipsoidal heights and geoid undulations between different geodetic reference frames

Transforming height information that refers to an ellipsoidal Earth reference model, such as the geometric heights determined from GPS measurements or the geoid undulations obtained by a gravimetric geoid solution, from one geodetic reference frame (GRF) to another is an important task whose proper implementation is crucial for many geodetic, surveying and mapping applications. This paper presents the required methodology to deal with the above problem when we are given the Helmert transformation parameters that link the underlying Cartesian coordinate systems to which an Earth reference ellipsoid is attached. The main emphasis is on the effect of GRF spatial scale differences in coordinate transformations involving reference ellipsoids, for the particular case of heights. Since every three-dimensional Cartesian coordinate system ‘gauges’ an attached ellipsoid according to its own accessible scale, there will exist a supplementary contribution from the scale variation between the involved GRFs on the relative size of their attached reference ellipsoids. Neglecting such a scale-induced indirect effect corrupts the values for the curvilinear geodetic coordinates obtained from a similarity transformation model, and meter-level apparent offsets can be introduced in the transformed heights. The paper explains the above issues in detail and presents the necessary mathematical framework for their treatment. An erratum to this article can be found at     

Terrain corrections to power 3 in gravimetric geoid determination

In precise geoid determination by Stokes formula, direct and primary and secondary indirect terrain effects are applied for removing and restoring the terrain masses. We use Helmert's second condensation method to derive the sum of these effects, together called the total terrain effect for geoid. We develop the total terrain effect to third power of elevation H in the original Stokes formula, Earth gravity model and modified Stokes formula. It is shown that the original Stokes formula, Earth gravity model and modified Stokes formula all theoretically experience different total terrain effects. Numerical results indicate that the total terrain effect is very significant for moderate topographies and mountainous regions. Absolute global mean values of 5–10 cm can be reached for harmonic expansions of the terrain to degree and order 360. In another experiment, we conclude that the most important part of the total terrain effect is the contribution from the second power of H, while the contribution from the third power term is within 9 cm. Received: 2 September 1996 / Accepted: 4 August 1997     

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     

A new method for computing the ellipsoidal correction for Stokes's formula

 This paper generalizes the Stokes formula from the spherical boundary surface to the ellipsoidal boundary surface. The resulting solution (ellipsoidal geoidal height), consisting of two parts, i.e. the spherical geoidal height N ₀ evaluated from Stokes's formula and the ellipsoidal correction N ₁, makes the relative geoidal height error decrease from O(e ²) to O(e ⁴), which can be neglected for most practical purposes. The ellipsoidal correction N ₁ is expressed as a sum of an integral about the spherical geoidal height N ₀ and a simple analytical function of N ₀ and the first three geopotential coefficients. The kernel function in the integral has the same degree of singularity at the origin as the original Stokes function. A brief comparison among this and other solutions shows that this solution is more effective than the solutions of Molodensky et al. and Moritz and, when the evaluation of the ellipsoidal correction N ₁ is done in an area where the spherical geoidal height N ₀ has already been evaluated, it is also more effective than the solution of Martinec and Grafarend. Received: 27 January 1999 / Accepted: 4 October 1999     

Meissel-Stokes核函数应用于区域大地水准面分析

为提高区域大地水准面计算精度,基于EGM2008地球重力场位系数模型分析Meissel-Stokes核函数、截断误差系数以及截断误差。选取实验区,采用移去-恢复法评价Meissel-Stokes核函数计算大地水准面的精度。结果表明:Meissel-Stokes核函数及其截断误差系数收敛速度快;截断误差小且稳定。在积分半径不易扩展的情况下,应用Meissel-Stokes核函数计算区域大地水准面,比标准Stokes计算大地水准面精度略高。     

Performance of three types of Stokes's kernel in the combined solution for the geoid

When regional gravity data are used to compute a gravimetric geoid in conjunction with a geopotential model, it is sometimes implied that the terrestrial gravity data correct any erroneous wavelengths present in the geopotential model. This assertion is investigated. The propagation of errors from the low-frequency terrestrial gravity field into the geoid is derived for the spherical Stokes integral, the spheroidal Stokes integral and the Molodensky-modified spheroidal Stokes integral. It is shown that error-free terrestrial gravity data, if used in a spherical cap of limited extent, cannot completely correct the geopotential model. Using a standard norm, it is shown that the spheroidal and Molodensky-modified integration kernels offer a preferable approach. This is because they can filter out a large amount of the low-frequency errors expected to exist in terrestrial gravity anomalies and thus rely more on the low-frequency geopotential model, which currently offers the best source of this information. Received: 11 August 1997 / Accepted: 18 August 1998     

Construction of Green's function for the Stokes boundary-value problem with ellipsoidal corrections in the boundary condition

Green's function for the boundary-value problem of Stokes's type with ellipsoidal corrections in the boundary condition for anomalous gravity is constructed in a closed form. The `spherical-ellipsoidal' Stokes function describing the effect of two ellipsoidal correcting terms occurring in the boundary condition for anomalous gravity is expressed in O(e <sup>2</sup> <sub>0</sub>)-approximation as a finite sum of elementary functions analytically representing the behaviour of the integration kernel at the singular point ψ=0. We show that the `spherical-ellipsoidal' Stokes function has only a logarithmic singularity in the vicinity of its singular point. The constructed Green function enables us to avoid applying an iterative approach to solve Stokes's boundary-value problem with ellipsoidal correction terms involved in the boundary condition for anomalous gravity. A new Green-function approach is more convenient from the numerical point of view since the solution of the boundary-value problem is determined in one step by computing a Stokes-type integral. The question of the convergence of an iterative scheme recommended so far to solve this boundary-value problem is thus irrelevant. Received: 5 June 1997 / Accepted: 20 February 1998     

A computational scheme to model the geoid by the modified Stokes formula without gravity reductions

In a modern application of Stokes formula for geoid determination, regional terrestrial gravity is combined with long-wavelength gravity information supplied by an Earth gravity model. Usually, several corrections must be added to gravity to be consistent with Stokes formula. In contrast, here all such corrections are applied directly to the approximate geoid height determined from the surface gravity anomalies. In this way, a more efficient workload is obtained. As an example, in applications of the direct and first and second indirect topographic effects significant long-wavelength contributions must be considered, all of which are time consuming to compute. By adding all three effects to produce a combined geoid effect, these long-wavelength features largely cancel. The computational scheme, including two least squares modifications of Stokes formula, is outlined, and the specific advantages of this technique, compared to traditional gravity reduction prior to Stokes integration, are summarised in the conclusions and final remarks. AcknowledgementsThis paper was written whilst the author was a visiting scientist at Curtin University of Technology, Perth, Australia. The hospitality and fruitful discussions with Professor W. Featherstone and his colleagues are gratefully acknowledged.     

On the solvability of the Stokes pseudo-boundary-value problem for geoid determination

A new form of boundary condition of the Stokes problem for geoid determination is derived. It has an unusual form, because it contains the unknown disturbing potential referred to both the Earth's surface and the geoid coupled by the topographical height. This is a consequence of the fact that the boundary condition utilizes the surface gravity data that has not been continued from the Earth's surface to the geoid. To emphasize the `two-boundary' character, this boundary-value problem is called the Stokes pseudo-boundary-value problem. The numerical analysis of this problem has revealed that the solution cannot be guaranteed for all wavelengths. We demonstrate that geoidal wavelengths shorter than some critical finite value must be excluded from the solution in order to ensure its existence and stability. This critical wavelength is, for instance, about 1 arcmin for the highest regions of the Earth's surface. Furthermore, we discuss various approaches frequently used in geodesy to convert the `two-boundary' condition to a `one-boundary' condition only, relating to the Earth's surface or the geoid. We show that, whereas the solution of the Stokes pseudo-boundary-value problem need not exist for geoidal wavelengths shorter than a critical wavelength of finite length, the solutions of approximately transformed boundary-value problems exist over a larger range of geoidal wavelengths. Hence, such regularizations change the nature of the original problem; namely, they define geoidal heights even for the wavelengths for which the original Stokes pseudo-boundary-value problem need not be solvable. Received 11 September 1995; Accepted 2 September 1996     

The IAG approach to the atmospheric geoid correction in Stokes' formula and a new strategy

The well-known International Association of Geodesy (IAG) approach to the atmospheric geoid correction in connection with Stokes' integral formula leads to a very significant bias, of the order of 3.2 m, if Stokes' integral is truncated to a limited region around the computation point. The derived truncation error can be used to correct old results. For future applications a new strategy is recommended, where the total atmospheric geoid correction is estimated as the sum of the direct and indirect effects. This strategy implies computational gains as it avoids the correction of direct effect for each gravity observation, and it does not suffer from the truncation bias mentioned above. It can also easily be used to add the atmospheric correction to old geoid estimates, where this correction was omitted. In contrast to the terrain correction, it is shown that the atmospheric geoid correction is mainly of order H of terrain elevation, while the term of order H ² is within a few millimetres. Received: 20 May 1998 / Accepted: 19 April 1999     

Review and numerical assessment of the direct topographical reduction in geoid determination

Recent papers in the geodetic literature promote the reduction of gravity for geoid determination according to the Helmert condensation technique where the entire reduction is made in place before downward continuation. The alternative approach, primarily developed by Moritz, uses two evaluation points, one at the Earths surface, the other on the (co-)geoid, for the direct topographic effect. Both approaches are theoretically legitimate and the derivations in each case make use of the planar approximation and a Lipschitz condition on height. Each method is re-formulated from first principles, yielding equations for the direct effect that contain only the spherical approximation. It is shown that neither method relies on a linear relationship between gravity anomalies and height (as claimed by some). Numerical tests, however, show that the practical implementations of these two approaches yield significant differences. Computational tests were performed in three areas of the USA, using 1×1 grids of gravity data and 30×30 grids of height data to compute the gravimetric geoid undulation, and GPS/leveled heights to compute the geometric geoid undulation. Using the latter as a control, analyses of the gravimetric undulations indicate that while in areas with smooth terrain no substantial differences occur between the gravity reduction methods, the Moritz–Pellinen (MP) approach is clearly superior to the Vanicek–Martinec (VM) approach in areas of rugged terrain. In theory, downward continuation is a significant aspect of either approach. Numerically, however, based on the test data, neither approach benefited by including this effect in the areas having smooth terrain. On the other hand, in the rugged, mountainous area, the gravimetric geoid based on the VM approach was improved slightly, but with the MP approach it suffered significantly. The latter is attributed to an inability to model the downward continuation of the Bouguer anomaly accurately in rugged terrain. Applying the higher-order, more accurate gravity reduction formulas, instead of their corresponding planar and linear approximations, yielded no improvement in the accuracy of the gravimetric geoid undulation based on the available data.     

Improved convergence rates for the truncation error in gravimetric geoid determination

 When Stokes's integral is used over a spherical cap to compute a gravimetric estimate of the geoid, a truncation error results due to the neglect of gravity data over the remainder of the Earth. Associated with the truncation error is an error kernel defined over these two complementary regions. An important observation is that the rate of decay of the coefficients of the series expansion for the truncation error in terms of Legendre polynomials is determined by the smoothness properties of the error kernel. Previously published deterministic modifications of Stokes's integration kernel involve either a discontinuity in the error kernel or its first derivative at the spherical cap radius. These kernels are generalised and extended by constructing error kernels whose derivatives at the spherical cap radius are continuous up to an arbitrary order. This construction is achieved by smoothly continuing the error kernel function into the spherical cap using a suitable degree polynomial. Accordingly, an improved rate of convergence of the spectral series representation of the truncation error is obtained. Received: 21 April 1998 / Accepted: 4 October 1999