# Stokes-Helmert's scheme for precise geoid determination.

ResumenEl esquema Stokes-Helmert para la determinacion precisa del geoide es un trabajo desarrollado en la Universidad de New Brunswick, Canada, por mas de diez anos de investigacion sin interrupcion, en la que han participado destacados cientificos. Este trabajo, presenta en forma sintetica los pasos para obtener un geoide centimetrico, asi como la formulacion matematica en la que se fundamenta este esquema.

Abstract

The Stokes-Helmert's scheme for the precise geoid determination is a work developed in the University of New Brunswick, Canada, by more than ten years of investigation without interruption, in which outstanding scientists have participated. This work, presents in synthetic form the steps to obtain a centimeter geoide, as well as the mathematical formulation on which east scheme is based.

Stokes-Helmert's geoid software

Stokes-Helmert's geoid software (SHGEO) is a scientific software for precise geoid determination based on the Stokes-Helmert theory of determination of the gravimetric geoid. The software has been developed during more then 10 years period under leadership of professor Petr Vanicek at the Department of Geodesy and Geomatics Engineering, University of New Brunswick. Authors of particular programs are: M. Najafi, P. Novak, J. Huang, J. Janak and R. Tenzer. We also have to mention Z. Martinec, A. Kleusberg, L.E. Sjoberg, W.E. Featherstone, W. Sun whose research presented in their papers was incorporated into the SHGEO software. SHGEO software uses various global models (e.g. TUG87, GRIM4-S4, EGM96). These global models play an important role in the geoid computation scheme. Therefore we acknowledge the contribution of all research teams that have developed these or other global models. Reference manual was compiled by R. Tenzer and J. Janak.

Stokes-Helmert's scheme for precise geoid determination

Introduction

This part of the manual gives a brief theoretical overview of the precise geoid determination process. The details can be found in references.

The Stokes-Helmert scheme for determination of the precise geoid can be summarized to the following steps:

1. Formulation of the boundary value problem on the Earth surface

2. Evaluation of the Helmert gravity anomalies on the Earth surface

3. Downward continuation of the Helmert gravity anomalies onto the geoid

4. Stokes's integration (solution to the Stokes's boundary value problem)

5. Transformation of geoidal heights from the Helmert space to the real space.

Formulation of the boundary value problem

The quantity to be solved is the earth's gravity potential U(r,[OMEGA]) on and outside the geoid and the geoid itself. The geoid is the equipotential surface that approximates the mean sea level most closely. The gravity potential on the geoid is denoted by [W.sub.o](r,[OMEGA]) = const. In order to solve this problema normal gravity potential U(r,[OMEGA]) generated by the mean geocentric ellipsoid of revolution is introduce& The normal gravity potential [U.sub.o] on the mean geocentric ellipsoid is chosen to be equal to the earth's potential on the geoid: [U.sub.o] = [W.sub.o].

The difference of the gravity potential W(r,[OMEGA]) and the normal gravity potential U(r,[OMEGA]) defines the disturbing potential T(r,[OMEGA]),

(1.1) T(r,[OMEGA]) = W(r,[OMEGA]) - U(r,[OMEGA])

When atmospheric attraction is neglected, T(r,[OMEGA]) is harmonic outside the Earth and it satisfies the Laplace equation

(1.2) [[DELTA].sup.2]T(r,[OMEGA]) = 0

Once T(r,[OMEGA]) has been solved, the gravity potential W(r,[OMEGA]) can be obtained at any point by adding U(r,[OMEGA]), which can be computed from existing models. Also when T(r,[OMEGA]) is known on the geoid, the vertical separation between the reference ellipsoid and the geoid can be obtained by the Bruns formula

(1.3) N([OMEGA])= T([r.sub.g]([OMEGA]))/[[gamma].sub.o]([OMEGA])

where T([r.sub.g]([OMEGA])) is the disturbing potential on the geoid, and [[gamma].sub.o]([OMEGA]) is the normal gravity on the mean geocentric ellipsoid. The problem is now reduced to the determination of T(r,[OMEGA]) on and outside the geoid.

However, the disturbing potential T(r,[OMEGA]) does not satisfy the Laplace equation inside of topographical masses where the geoid is often located. Therefore in order to satisfy Laplace's equation, all atmospheric and topographical masses have to be removed or condensed on or beneath the geoid. In Helmert's second condensation method the atmospheric and topographical masses are condensed directly onto the geoid.

When atmospheric and topographical masses are condensed as a single layer that is located on the geoid, the Earth gravity field will slightly change. The space obtained after such a condensation is the Helmert space. The quantities given in the Helmert space are denoted by superscript H. Helmert's gravity potential is defined as follows

(1.4) [W.sup.H](r,[OMEGA]) = W(r,[OMEGA]) - [delta][V.sup.t](r,[OMEGA]) - [delta][V.sup.a] (r,[OMEGA])

The residual topographical potential [delta][V.sup.t](r,[OMEGA]) is defined as a difference of the gravitational potential [V.sup.t](r,[OMEGA]) of topographical masses and the gravitational potential [V.sup.ct](r,[OMEGA]) of condensed topographical masses

(1.5) [delta][V.sup.t](r,[OMEGA]) = [V.sup.t](r,[OMEGA]) - [V.sup.ct](r,[OMEGA])

Similarly, the residual atmospheric potential [delta][V.sup.a](r,[OMEGA]) is defined as a difference of the gravitational potential [V.sup.a](r,[OMEGA]) of atmospheric masses and the gravitational potential [V.sup.ca](r,[OMEGA]) of condensed atmospheric masses

(1.6) [delta][V.sup.a](r,[OMEGA]) = [V.sup.a](r,[OMEGA]) - [V.sup.ca](r,[OMEGA])

By subtracting the normal gravity potential U(r,[OMEGA]) from eqn. (1.4), the disturbing potential [T.sup.H](r,[OMEGA]) in Helmert's space becomes

(1.7) [T.sup.H](r,[OMEGA]) = [W.sup.H](r,[OMEGA]) - U(r,[OMEGA])

Helmert's disturbing potential [T.sup.H](r,[OMEGA]) is harmonic above the geoid, so that it satisfies the Laplace equation

(1.8) [[DELTA].sup.2][T.sup.H](r,[OMEGA]) = 0

To determine [T.sup.H](r,[OMEGA]), the boundary value problem of the third kind outside the geoid has to be solved. Therefore, the boundary values on the a-priori unknown geoid are needed. In this problem, the Helmert gravity anomalies on the geoid serve as the boundary values. To find a relation between the disturbing potential and the Helmert gravity anomalies, let us introduce the radial derivative of the Helmert disturbing potential

(1.9) [delta][T.sup.H](r,[OMEGA])/[delta]r = [delta][W.sup.H](r,[OMEGA])/[delta]r - [delta]U(r,[OMEGA])/[delta]r

The negative radial derivative of the Helmert disturbing potential [T.sup.H](r,[OMEGA]) defines the Helmert gravity disturbance [delta][g.sup.H] (r,[OMEGA]), i.e.,

(1.10) - [delta][T.sup.H](r,[OMEGA])/[delta]r = [delta][g.sup.H](r,[OMEGA]) - [[epsilon].sub.[delta]g](r,[OMEGA])

The second term on the right hand side of eqn. (1.10) is the ellipsoidal correction to the gravity disturbance (Vanicek et al., 1999).

Since the geoidal height above the ellipsoid are not usually available, the gravity disturbance is not considered to be a measurable quantity on the surface of the Earth. Therefore Helmert's gravity disturbance [delta][g.sup.H](r,[OMEGA]) has to be transformed to more commonly available quantity, which is the Helmert gravity anomaly [DELTA][g.sup.H](r,[OMEGA]). This transformation is achieved by adding a term [GAMMA](r,[OMEGA]) to the gravity disturbance. This term accounts for the change in normal gravity due to the difference between the geodetic height h([OMEGA]) and the commonly available orthometric height [H.sup.o]([OMEGA]). This expression can be written with a sufficient accuracy as

(1.11) [GAMMA](r,[OMEGA]] = [T.sup.H](r,[OMEGA])/[gamma](r,[OMEGA]) [delta] [gamma](r,[OMEGA])/[delta]n

In solving the boundary value problem it is convenient to introduce the following spherical approximation

(1.12) 1/[gamma](r,[OMEGA]) [delta] [gamma](r,[OMEGA])/[delta]n [congruent to] - 2/R

The error caused by this approximation is called the ellipsoidal correction [[epsilon].sub.n] (r,[OMEGA]) for the spherical approximation.

Substituting eqns. (1.11) and (1.12) into eqn. (1.10), we can formulate the boundary value problem in the Helmert space by the following equation

(1.13) - [delta][T.sup.H](r,[OMEGA])/[delta]r - 2/R [T.sup.H](r,[OMEGA]) = [DELTA][g.sup.H](r,[OMEGA]) + [[epsilon].sub.n], (r,[OMEGA]) - [[epsilon].sub.[delta]g](r,[OMEGA])

This equation represents the fundamental equation of the physical geodesy valid in the Helmert space. It relates the known boundary values [DELTA][g.sup.H](r,[OMEGA]) to the unknown disturbing potential [T.sup.H](r,[OMEGA]) outside and on the geoid. As the boundary values are the Helmert gravity anomalies on the geoid, the next two sections deal with the derivation of this quantity from the measurements.

Evaluation of the Helmert gravity anomalies on the earth surface

Helmert gravity anomaly [DELTA][g.sup.H]([r,([OMEGA])) referred on the earth surface is given by

(1.14) [DELTA][g.sup.H]([r.sub.t]([OMEGA])) = [g.sup.H]([r.sub.t]([OMEGA])) - [gamma]([H.sup.N]([OMEGA]))

where [gamma]([H.sup.N]([OMEGA])) is the normal gravity on the telluriod in the Helmert space (see Figure 1).

[FIGURE 1 OMITTED]

The Helmert gravity [g.sup.H]([r.sub.t]([OMEGA])) on the earth surface is obtained from the observed gravity g([r.sub.t]([OMEGA])), by adding the direct topographical effect and the direct atmospheric effect:

(1.15) [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

The direct topographical effect (DTE) on gravity, given by a radial derivative of the residual gravitational potential [delta][V.sup.t](r,[OMEGA]) of topographical masses, is the gravitational attraction of topographical masses minus the gravitational attraction of condensed topographical masses. It should be evaluated on the earth surface. Analogous, the direct atmospheric effect (DAE) on gravity is gravitational attraction of atmospheric masses minus the gravitational attraction of condensed atmospherical masses.

Normal gravity [gamma]([H.sup.N]([OMEGA])) on the telluriod in Helmert's space is obtained by upward continuation of normal gravity at the corresponding point on the mean geocentric ellipsoid. The height used for this computation should be the normal height [H.sup.N]([OMEGA]) in Helmert's space, which is the height of Helmert's telluroid above the ellipsoid. In practice, however, the heights on gravity observations on the earth surface are orthometric heights [H.sup.o]([OMEGA]) and the upward continuation of normal gravity is computed using [H.sup.o]([OMEGA]) instead of [H.sup.N]([OMEGA]). Therefore, a correction has to be added to the normal gravity. This correction consists of two parts. The first one, which is caused by the different position of the telluroid in the real and Helmert's space, is called the secondary indirect effect. The second part, which is due to the difference between the normal and orthometric height, is called the correction for the orthometric height. The expression for the normal gravity on Helmert's telluroid is then given by

(1.16) [[gamma].sup.H]([H.sup.N]([OMEGA])) = [gamma]([H.sup.o)([OMEGA])) - 2/R [H.sup.o]([OMEGA])[DELTA][g.sup.SB]([r.sub.t]([OMEGA])) - 2/[r.sub.t]([OMEGA]) [delta][V.sup.t](r.sub.t]([OMEGA])) - 2/[r.sub.t]([OMEGA])) [delta][V.sup.a]([r.sub.t]([OMEGA]))

where [DELTA][g.sup.SB]([r.sub.t]([OMEGA])) is the simple Bouguer gravity anomaly.

The second term on the right-hand side of eqn. (1.16) is the geoid-quasigeoid correction to the boundary value problem. The third and fourth terms stands for the secondary indirect topographical effect on gravity (SITE) and the secondary indirect atmospheric effect on gravity (SIAE).

Inserting eqns. (1.14) and (1.15) back into eqn. (1.14) and considering also the free-air gravity anomaly given by

(1.17) [DELTA][g.sup.FA]([r.sub.t]([OMEGA])) = g([r.sub.t]([OMEGA])) - [gamma]([H.sup.o]([OMEGA]))

the boundary value problem in the Helmert space (eqn. (1.14)) takes the following form

(1.18) [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Free-air gravity anomalies are not very smooth, so that are not suitable for interpolation and averaging. Therefore, in practice the smoother complete Bouguer gravity anomalies [DELTA][g.sup.CB] ([r.sub.t]([OMEGA])) are used to produce the mean gravity data. The complete Bouguer anomaly is defined as

(1.19) [DELTA][g.sup.CB](r,[OMEGA]) = [DELTA][g.sup.SB](r,[OMEGA]) + [delta][g.sup.tc](r,[OMEGA]) = [DELTA][g.sup.FA](r, [OMEGA]) - 2[pi][[rho].sub.o]G[H.sup.o]([OMEGA]) + [delta][g.sup.tc](r,[OMEGA])

where [delta][g.sup.tc](r,[OMEGA]) is the gravimetric terrain correction, i.e. a correction for the attraction of the earth's topography relative to the height at the evaluation point.

Downward continuation of Helmert gravity anomaly

To solve the Stokes boundary value problem in the Helmert space the Helmert's gravity anomalies [DELTA][g.sup.H]([r.sub.t] ([OMEGA])) have to be downward continued onto the geoid.

The Helmert disturbing potential is a harmonic function above the Helmert co-geoid. Poisson's solution to Dirichlet's problem of upward continuation of a harmonic function can be applied in finding a solution to the inverse problem, i.e. downward continuation.

The relation between Helmert's gravity anomalies [DELTA][g.sup.H] ([r.sub.g]([OMEGA])) on the geoid and Helmert's gravity anomalies [DELTA][g.sup.H]([r.sub.t] ([OMEGA])) on the earth surface is given by Poisson's integral

(1.20) [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

where K[[r.sub.t]([OMEGA]),[psi]([OMEGA],[OMEGA]'),R] is the Poisson integral kernel. The solution to the problem of downward continuation is then given as a inverse solution of an integral equation of the first kind, where [DELTA][g.sup.H]([r.sub.t]([OMEGA])) is known and [DELTA][g.sup.H] ([r.sub.g]([OMEGA])) is being determined. This integral equation can be solved iteratively according to Jacobi's iteration approach.

Stokes's boundary value problem

According to eqn. (1.3), the geoidal height N([OMEGA]) can be obtained from the Bruns formula if the disturbing gravity potential T([r.sub.g]([OMEGA])) on the geoid is known. The relation between the gravity anomaly on the geoid and the disturbing gravity potential referred on the geoid is given by Stokes theorem.

Solving Stokes's boundary value problem in the Helmert space, the co-geoid heights [N.sup.H]([OMEGA]) is given by the Stokes integral

(1.21) [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

where S([psi]([OMEGA],[OMEGA]')) is the Stokes integral kernel.

Transformation of geoidal heights from the Helmert space to the real space

To obtain the geoid, the co-geoidal heights [N.sup.H]([OMEGA]) have to be transformed from the Helmert space back to the real space. This transformation is done by evaluating the primary indirect topographical effect (PITE) on the geoidal height according to the following equation

(1.22) [N.sup.t]([OMEGA]) = [delta][V.sup.t]([r.sub.g]([OMEGA]))/[[delta].sub.o]([OMEGA])

where [delta][V.sup.t]([r.sub.g]([OMEGA])) is the residual gravitational potential of topographical masses, which is reckoned at the geoid.

The geoid is finally computed by

(1.23) N([OMEGA]) = [N.sup.H]([OMEGA]) + [N.sup.t]([OMEGA])

The primary indirect atmospheric effect as well as the secondary indirect atmospheric effect can be neglected, as it is shown in (Novak, 2000). Equation (1.23) is then final equation in the precise geoid determination process.

References

Heiskanen, W.A.; Moritz, H. (1967). Physical Goedesy, Freeman, San Francisco.

Huang, J.; Vanicek; P.; Novak, P. (2000). "An alternative algorithm to FFF for the numerical evaluation of Stokes's integral," Studia Geophysica et Geodaetica 44, pp. 374-380.

Huang, J.; Vanicek, E; Pagiatakis, S.; Brink, W. (2001). "Effect of topographical mass density variation on gravity and the geoid in the Canadian Rocky mountains," Journal of Geodesy 74 (11-12), pp. 805-815.

-- (2001). Computational Methods for the Discrete Downward Continuation of the Earth Gravity, IAG General Assembly, Budapest.

Janak, J.; Vanicek, P.; Alberts, B. (2001). Point and mean values of topographical effect, the Digital Earth conference, Fredericton, June 25-28.

Kellogg O.D. (1929). Foundations of potential theory, Springer, Berlin, reprinted by Dover, New York in 1953.

Martinec, Z. (1993). Effect of lateral density variations of topographical masses in improving geoid model accuracy over Canada. Contact report for Geodetic Survey of Canada, Ottawa.

Martinec, Z.; Matyska, C.; Grafarend, E.W.; Vanicek, E (1993). "On Helmert's 2nd condensation method," Manuscripta Geodaetica 18, pp. 417-421.

Martinec, Z.; Vanicek, P. (1994). "The indirect effect of topography in the Stokes-Helmert technique for a spherical approximation of the geoid," Manuscripta Geodaetica 19, pp. 213-219.

-- (1994). "Direct topographical effect of Helmert's condensation for a spherical geoid," Manuscripta Geodaetica 19, pp. 257-268.

-- (1996). Formulation of the boundary-value problem for geoid determination with a higher-degree reference field, Geophys. J. Int. 126, pp. 219-228.

Martinec, Z.; Vanicek, P.; Mainvile, A.; Veronneau, M. (1995). "The effect of lake water on geoidal heights," Manuscripta Geodaetica 20, pp. 193-203.

-- (1996). "Evaluation of topographical effects in precise geoid determination from densely sampled heights," Journal of Geodesy, 70 (11), pp. 746-754.

Martinec, Z. (1996). "Stability investigations of a discrete downward continuation problem for geoid determination in the Canadian Rocky Mountains," Journal of Geodesy, 70/11, 805-828.

-- (1998). Geodetic Boundary Value problem, Springer.

Molodenskij, M.S.; Eremeev, V.F.; Yurkina, M.I. (1960). Methods for Study of the External Gravitational Field and Figure of the Earth, translated from Russian by the Israel Program for Scientific Translations for the Office of Technical Services, US Department of Commerce, Washington, DC., USA, 1962.

Moritz, H. (1980). Advanced Physical Geodesy, H. Wichmann, Karlsruhe.

Novak, P.; Vanicek, P.; Martinec, Z.; Veronneau, M. (2001) "The effect of distant terrain on gravity and the geoid," Journal of Geodesy (in press).

Novak, P.; Vanicek, P.; Veronneau, M.; Featherstone, W.E.; Holmes, S.A. (2001). "On the accuracy of modified Stokes's integration in high-frequency gravimetric geoid determination: A comparison of two numerical techniques," Journal of Geodesy 74, pp. 644-654.

Sjoberg, L.E. (1996). "On the error of analytical continuation in physical geodesy", Journal of Geodesy 70, pp. 724-730.

Stokes, G.G. (1849). On the variation of gravity at the surface of the earth, Trans. Cambridge Philos. Soc. VII, pp. 672-695.

Sun, W.; Vanicek, P. (1996). "On the discrete problem of downward Helmert's gravity continuation," Proceedings of Session G7 (Technikues for local geoid determination), Annual meeting of European Geophysical Society, The Hague, May 6-10, 1996, Reports of the Finnish Geodetic Institute, 96:2, pp. 29-34.

--, (1998). "On some problems of the downward continuation of 5' x 5' mean Helmert's gravity disturbance," Journal of Geodesy 72, 7-8, pp. 411-420.

Vanicek, P.; Krakiwsky, E.J. (1986). Geodesy: the Concepts, 2nd rev. ed., North-Holland, Amsterdam.

Vanicek, P.; Kleusberg, A. (1987). "The Canadian geoid--Stokesian approach," Manuscripta Geodaetica 12, pp. 86-98.

Vanicek, P.; Sjoberg, L.E. (1991). Reformulation of Stokes's theory for higher than second degree reference field and modification of integration kernels, JGR 96 (B4), pp. 6529-6539.

Vanicek, P.; Martinec, Z. (1994). "Stokes-Helmert scheme for the evaluation of a precise geoid," Manuscripta Geodaetica 19, pp. 119-128.

Vanicek, P.; Najafi, M.; Martinec, Z.; Harrie, L.; Sjoberg, L.E. (1996). "Higher-order reference field in the generalized Stokes-Helmert scheme for geoid computation," Journal of Geodesy 70 (3), pp. 176-182.

Vanicek, P.; Sun, W.; Ong, P.; Martinec, Z.; Najafi, M.; Vajda, P.; B ter Horst (1996). "Downward continuation of Helmert's gravity," Journal of Geodesy.

Vanicek, P.; Featherstone, W.E. (1998). "Performance of three types of Stokes's kernel in the combined solution for the geoid," Journal of Geodesy 72, 12, pp. 684-697.

Vanicek, P.; Huang, J.; Novak, P.; Veronneau, M.; Pagiatakis, S.; Martinec, Z.; Featherstone, W.E. (1999). "Determination of boundary values for the Stokes-Helmert problem," Journal of Geodesy 73, pp. 180-192.

Vanicek, P.; Novak, P. (1999). "Comparison between planar and spherical models of topography," CGU Annual Meeting, Banff, May 9-13.

Vanicek, P.; Janak, J. (2000). "The UNB technique for precise geoid determination," CGU Meeting, Banff, May 24-26.

Wichiencharoen, C. (1982). The indirect effect on the computation of geoid undulations, Dept. of Geodetic Science, Rep. 336, Ohio State University.

Robert Tenzer *

Juraj Janak **

* University of New Brunswick. Email: rtenzer@unb.ca

** Department of Theoretical Geodesy, Faculty of Civil Engineering, Slovak University of Technology, Radlinskeho 11, 81368 BRATISLAVA, Slovak Republic. E-mail: janak@aries.svf.stuba.sk

Printer friendly Cite/link Email Feedback | |

Author: | Tenzer, Robert; Janak, Juraj |
---|---|

Publication: | Revista Cartografica |

Date: | Jan 1, 2002 |

Words: | 3428 |

Previous Article: | La obra escrita de los ingenieros geografos mexicanos. |

Next Article: | Glosario GPS y terminos afines. |

Topics: |