The Level Ellipsoid


Table of Contents

 

Introduction

The ellipsoid provides a mathematical surface which closely resembles the earth's shape, and is used in surveying with plane coordinate systems and in geodesy. For a point, it also provides two dimensions, latitude and longitude, of the point's three-dimensional position. The third dimension, gravity, can be approximated by assigning an approximate gravitational field to the ellipsoid. This field should closely approximate the actual gravity field of the earth. In so defining this equipotential surface, the direction of a plumb line becomes perpendicular to the ellipsoid and thus is known as the normal gravity field. This surface is sometimes known as the level ellipsoid.

 

The Normal Gravity Field

The normal figure of the earth should guarantee a


The rotational ellipsoid (spheroid) meets these needs. In addition to the standard ellipsoidal parameters of a and f, the terms of total mass, M, and rotational angular velocity are needed to define this surface. The gravity field is developed from the angular rotation and gravitation. The level surface is the surface of the ellipsoid.

Theorem of Stokes: If a body of total mass M rotates with constant angular velocity, w (omega), about a fixed axis, and if S is a level surface of its gravity field enclosing the entire mass, then the gravity potential in the exterior space of S is uniquely determined by M, w, S, a, f.

Figure 3.11 The figure to the right shows a set of confocal ellipsoids with constant linear eccentricity e. The level surfaces coincide with surfaces of equal density and equal pressure. The hydrostatic equilibrium of the normal figure is created by redistribution of the actual masses of the Earth. According to theory, the closed representation of the normal gravity potential may be achieved in a system of ellipsoid coordinates. The point P is then specified by ellipsoidal coordinates u (the semi-minor axis), b (the reduced latitude), and l (the geographic longitude). From these terms, we can derive that the semi-major axis is

The transformation into Cartesian coordinates is

U is composed of gravitational potential and the potential of centrifugal acceleration and is computed as

U = V + F

The gravitational ellipsoid satisfies DV = 0 in the space exterior to the ellipsoid containing mass M. Imposing rotational symmetry on the normal gravity field, the nonzonal terms in the spherical harmonic expansion disappear. In exterior space, closed expression for the potential of normal gravity is

where q is

and q0 = qu=b.


The surfaces of constant potential, that is U = constant, are termed spheropotential surfaces (sherops) are not ellipsoids (excepting the surface U0). The equation for the normal gravity on the ellipsoid is

where a and b are the semi-major and semi-minor axes of the ellipsoid, respectively, are the gravity at the equator and poles, respectively. The equation for the normal gravity field above the ellipsoid is given by Heiskanen and Moritz (1967) as

where h denotes the height above the ellipsoid in meters.



Thereom of Pizzetti

Thereom of Clairaut

 

where b is the gravity flattening, , and e' is the second eccentricity. From the above two theorems, it can be seen that there are only four independent quantities.



Series Expansions in the Normal Gravity Field

The potential of normal gravity is

For n=2 with the substitution of , we get

Solving for r and setting U = U0 gives the radius vector to the level ellipsoid where r = a, or

and the normal gravity with respect to r is

Substituting in a latitude of 0 and 90 degrees, the above equations yield a semi-major axis of a and equatorial gravity of ga. Also from these equations we see that the geometric flattening f is

and the gravity flattening is

where

is the ratio of centrifugal acceleration to the normal gravity at the equator. Using the above equations, an approximation for the theorem of Pizzetti is

and to Clairaut's theorem is

From these formulas we obtain Newton's gravity formula

Thus if two gravity values are known on the ellipsoid and at different geodetic latitudes, ga and b can be computed. Similarly with known values for a and angular velocity, the quantity m can be computed. The following series expansions are also useful

Notice that in the above equations, the definition for the normal gravity potential is defined by a function of a (the earth's semi-major axis), J2 (dynamic form factor), GM (the product of Newton's gravitational constant and the earth's mass), and the earth's angular velocity w. For GRS80 these constants are

  1. a = 6,378,137 m
  2. GM = 3,986,005 ´ 108 m3/sec2
  3. J2 = 108,263 ´ 10-8
  4. w = 7,292,115 ´ 10-11 rad/sec

The derived constants for GRS80 are

  1. b = 6,356,752.3141 m
  2. 1/f = 298.257222101
  3. m = 0.00344978600308
  4. ga = 9.7803267715 m/sec2
  5. gb = 9.8321863685 m/sec2

The equation for the normal gravity field above the ellipsoid is given by Heiskanen and Moritz (1967) as

where h denotes the height above the ellipsoid.


Updated by Charles D. Ghilani, Ph.D.
Last Update: October 04, 2006
Penn State Surveying Program, Copyright © 2000 - 2007