# Earth ellipsoid

An Earth ellipsoid is a mathematical figure approximating the shape of the Earth, used as a reference frame for computations in geodesy, astronomy and the geosciences. Various different ellipsoids have been used as approximations.

It is an ellipsoid of revolution, whose short (polar) axis (connecting the two flattest spots called geographical north and south poles) is approximately aligned with the rotation axis of the Earth. The ellipsoid is defined by the equatorial axis a and the polar axis b; their difference is about 21 km or 0,3 per cent. Additional parameters are the mass function J2, the correspondent gravity formula, and the rotation period (usually 86164 seconds).

Many methods exist for determination of the axes of an Earth ellipsoid, ranging from meridian arcs up to modern satellite geodesy or the analysis and interconnection of continental geodetic networks. Amongst the different set of data used in national surveys are several of special importance: the Bessel ellipsoid of 1841, the international Hayford ellipsoid of 1924, and (for GPS positioning) the WGS84 ellipsoid.

## Historical method of determining the ellipsoid

High precision land surveys can be used determine the distance between two places at nearly the same longitude by measuring a base line and a chain of triangles. (Suitable stations for the end points are rarely at the same longitude). The distance Δ along the meridian from one end point to a point at the same latitude as the second end point is then calculated by trigonometry. The surface distance Δ is reduced to Δ', the corresponding distance at mean sea level. The intermediate distances to points on the meridian at the same latitudes as other stations of the survey may also be calculated.

The geographic latitudes of both end points, φs (standpoint) and φf (forepoint) and possibly at other points are determined by astrogeodesy, observing the zenith distances of sufficient numbers of stars. If latitudes are measured at end points only, the radius of curvature at the midpoint of the meridian arc can be calculated from R = Δ'/(|φsf|). A second meridian arc will allow the derivation of two parameters required to specify a reference ellipsoid. Longer arcs with intermediate latitude determinations can completely determine the ellipsoid. In practice multiple arc measurements are used to determine the ellipsoid parameters by the method of least squares. The parameters determined are usually the semi-major axis, $a$, and either the semi-minor axis, $b$, or the inverse flattening $1/f$, (where the flattening is  $f=(a-b)/a$).

Geodesy no longer uses simple meridian arcs, but complex networks with hundreds of fixed points linked by the methods of satellite geodesy.

## Historical Earth ellipsoids

The reference ellipsoid models listed below have had utility in geodetic work and many are still in use. The older ellipsoids are named for the individual who derived them and the year of development is given. In 1887 the English mathematician Col Alexander Ross Clarke CB FRS RE was awarded the Gold Medal of the Royal Society for his work in determining the figure of the Earth. The international ellipsoid was developed by John Fillmore Hayford in 1910 and adopted by the International Union of Geodesy and Geophysics (IUGG) in 1924, which recommended it for international use.

At the 1967 meeting of the IUGG held in Lucerne, Switzerland, the ellipsoid called GRS-67 (Geodetic Reference System 1967) in the listing was recommended for adoption. The new ellipsoid was not recommended to replace the International Ellipsoid (1924), but was advocated for use where a greater degree of accuracy is required. It became a part of the GRS-67 which was approved and adopted at the 1971 meeting of the IUGG held in Moscow. It is used in Australia for the Australian Geodetic Datum and in South America for the South American Datum 1969.

The GRS-80 (Geodetic Reference System 1980) as approved and adopted by the IUGG at its Canberra, Australia meeting of 1979 is based on the equatorial radius (semi-major axis of Earth ellipsoid) $a$, total mass $GM$, dynamic form factor $J_2$ and angular velocity of rotation $\omega$, making the inverse flattening $1/f$ a derived quantity. The minute difference in $1/f$ seen between GRS-80 and WGS-84 results from an unintentional truncation in the latter's defining constants: while the WGS-84 was designed to adhere closely to the GRS-80, incidentally the WGS-84 derived flattening turned out to be slightly different than the GRS-80 flattening because the normalized second degree zonal harmonic gravitational coefficient, that was derived from the GRS-80 value for J2, was truncated to 8 significant digits in the normalization process.[1]

An ellipsoidal model describes only the ellipsoid's geometry and a normal gravity field formula to go with it. Commonly an ellipsoidal model is part of a more encompassing geodetic datum. For example, the older ED-50 (European Datum 1950) is based on the Hayford or International Ellipsoid. WGS-84 is peculiar in that the same name is used for both the complete geodetic reference system and its component ellipsoidal model. Nevertheless the two concepts—ellipsoidal model and geodetic reference system—remain distinct.

Note that the same ellipsoid may be known by different names. It is best to mention the defining constants for unambiguous identification.

Reference ellipsoid name Equatorial radius (m) Polar radius (m) Inverse flattening Where used
Maupertuis (1738) 6,397,300 6,363,806.283 191 France
Plessis (1817) 6,376,523.0 6,355,862.9333 308.64 France
Everest (1830) 6,377,299.365 6,356,098.359 300.80172554 India
Everest 1830 Modified (1967) 6,377,304.063 6,356,103.0390 300.8017 West Malaysia & Singapore
Everest 1830 (1967 Definition) 6,377,298.556 6,356,097.550 300.8017 Brunei & East Malaysia
Airy (1830) 6,377,563.396 6,356,256.909 299.3249646 Britain
Bessel (1841) 6,377,397.155 6,356,078.963 299.1528128 Europe, Japan
Clarke (1866) 6,378,206.4 6,356,583.8 294.9786982 North America
Clarke (1878) 6,378,190 6,356,456 293.4659980 North America
Clarke (1880) 6,378,249.145 6,356,514.870 293.465 France, Africa
Helmert (1906) 6,378,200 6,356,818.17 298.3
Hayford (1910) 6,378,388 6,356,911.946 297 USA
International (1924) 6,378,388 6,356,911.946 297 Europe
Krassovsky (1940) 6,378,245 6,356,863.019 298.3 USSR, Russia, Romania
WGS66 (1966) 6,378,145 6,356,759.769 298.25 USA/DoD
Australian National (1966) 6,378,160 6,356,774.719 298.25 Australia
New International (1967) 6,378,157.5 6,356,772.2 298.24961539
GRS-67 (1967) 6,378,160 6,356,774.516 298.247167427
South American (1969) 6,378,160 6,356,774.719 298.25 South America
WGS-72 (1972) 6,378,135 6,356,750.52 298.26 USA/DoD
GRS-80 (1979) 6,378,137 6,356,752.3141 298.257222101 Global ITRS[2]
WGS-84 (1984) 6,378,137 6,356,752.3142 298.257223563 Global GPS
IERS (1989) 6,378,136 6,356,751.302 298.257
IERS (2003)[3] 6,378,136.6 6,356,751.9 298.25642 [2]

## Mean Earth ellipsoid and reference ellipsoids

A data set which describes the global average of the Earth's surface curvature is called the mean Earth Ellipsoid. It refers to a theoretical coherence between the geographic latitude and the meridional curvature of the geoid. The latter is close to the mean sea level, and therefore an ideal Earth ellipsoid has the same volume as the geoid.

While the mean Earth ellipsoid is the ideal basis of global geodesy, for regional networks a so-called reference ellipsoid may be the better choice.[4] When geodetic measurements have to be computed on a mathematical reference surface, this surface should have a similar curvature as the regional geoid. Otherwise the reduction of the measurements would get small distortions.

This is the reason for the "long life" of former reference ellipsoids like the Hayford or the Bessel ellipsoid, despite of the fact that their main axes deviate by several hundred meters from the modern values. Another reason is a juridical one: the coordinates of millions of boundary stones should remain fixed for a long period. If their reference surface would change, the coordinates themselves would also change.

However, for international networks, GPS positioning or astronautics, these regional reasons are less relevant. As the knowledge of Earth's figure is increasingly accurate, the International Geoscientific Union IUGG usually adopts the axes of the Earth ellipsoid to the best available data.

## References

1. ^ NIMA Technical Report TR8350.2, "Department of Defense World Geodetic System 1984, Its Definition and Relationships With Local Geodetic Systems", Third Edition, 4 July 1997 [1]
2. ^ a b Note that the current best estimates, given by the IERS Conventions, "should not be mistaken for conventional values, such as those of the Geodetic Reference System GRS80 ... which are, for example, used to express geographic coordinates" (chap. 1); note further that "ITRF solutions are specified by Cartesian equatorial coordinates X, Y and Z. If needed, they can be transformed to geographical coordinates (λ, φ, h) referred to an ellipsoid. In this case the GRS80 ellipsoid is recommended." (chap. 4).
3. ^ IERS Conventions (2003) (Chp. 1, page 12)
4. ^ Alexander, J. C. (1985). "The Numerics of Computing Geodetic Ellipsoids". SIAM Review 27 (2): 241. doi:10.1137/1027056. edit