Geodesy
Fundamentals Geodesy Geodynamics Geomatics History
Concepts Geographical distance Geoid Figure of the Earth (radius and circumference) Geodetic coordinates Geodetic datum Geodesic Horizontal position representation Latitude / Longitude Map projection Reference ellipsoid Satellite geodesy Spatial reference system Spatial relations Vertical positions
Technologies Global Nav. Sat. Systems (GNSSs) Global Pos. System (GPS) GLONASS (Russia) BeiDou (BDS) (China) Galileo (Europe) NAVIC (India) Quasi-Zenith Sat. Sys. (QZSS) (Japan) Discrete Global Grid and Geocoding
Standards (history) NGVD 29 Sea Level Datum 1929 OSGB36 Ordnance Survey Great Britain 1936 SK-42 Systema Koordinat 1942 goda ED50 European Datum 1950 SAD69 South American Datum 1969 GRS 80 Geodetic Reference System 1980 ISO 6709 Geographic point coord. 1983 NAD 83 North American Datum 1983 WGS 84 World Geodetic System 1984 NAVD 88 N. American Vertical Datum 1988 ETRS89 European Terrestrial Ref. Sys. 1989 GCJ-02 Chinese obfuscated datum 2002 Geo URI Internet link to a point 2010 International Terrestrial Reference System Spatial Reference System Identifier (SRID) Universal Transverse Mercator (UTM)
v t e

The geoid (/ˈdʒiː.ɔɪd/) is the shape that the ocean surface would take under the influence of the gravity of Earth, including gravitational attraction and Earth's rotation, if other influences such as winds and tides were absent. This surface is extended through the continents (such as with very narrow hypothetical canals). According to Gauss, who first described it, it is the "mathematical figure of the Earth", a smooth but irregular surface whose shape results from the uneven distribution of mass within and on the surface of Earth.^[1] It can be known only through extensive gravitational measurements and calculations. Despite being an important concept for almost 200 years in the history of geodesy and geophysics, it has been defined to high precision only since advances in satellite geodesy in the late 20th century.

All points on a geoid surface have the same geopotential (the sum of gravitational potential energy and centrifugal potential energy). The force of gravity acts everywhere perpendicular to the geoid, meaning that plumb lines point perpendicular and bubble levels are parallel to the geoid. Being an equigeopotential means the geoid corresponds to the free surface of water at rest (if only gravity and rotational acceleration were at work); this is also a sufficient condition for a ball to remain at rest instead of rolling over the geoid. Earth's gravity acceleration (the vertical derivative of geopotential) is thus non-uniform over the geoid.^[2] The geoid undulation or geoidal height is the height of the geoid relative to a given reference ellipsoid. The geoid serves as a coordinate surface for various vertical coordinates, such as orthometric heights, geopotential heights, and dynamic heights (see Geodesy#Heights).

Description

The geoid surface is irregular, unlike the reference ellipsoid (which is a mathematical idealized representation of the physical Earth as an ellipsoid), but is considerably smoother than Earth's physical surface. Although the "ground" of the Earth has excursions on the order of +8,800 m (Mount Everest) and −11,000 m (Marianas Trench), the geoid's deviation from an ellipsoid ranges from +85 m (Iceland) to −106 m (southern India), less than 200 m total.^[3]

If the ocean were of constant density and undisturbed by tides, currents or weather, its surface would resemble the geoid. The permanent deviation between the geoid and mean sea level is called ocean surface topography. If the continental land masses were crisscrossed by a series of tunnels or canals, the sea level in those canals would also very nearly coincide with the geoid. In reality, the geoid does not have a physical meaning under the continents, but geodesists are able to derive the heights of continental points above this imaginary, yet physically defined, surface by spirit leveling.

Being an equipotential surface, the geoid is, by definition, a surface upon which the force of gravity is perpendicular everywhere. That means that when traveling by ship, one does not notice the undulations of the geoid; the local vertical (plumb line) is always perpendicular to the geoid and the local horizon tangential to it. Likewise, spirit levels will always be parallel to the geoid.

Simplified example

Earth's gravitational field is not uniform. An oblate spheroid is typically used as the idealized Earth, but even if the Earth were spherical and did not rotate, the strength of gravity would not be the same everywhere because density varies throughout the planet. This is due to magma distributions, the density and weight of different geological compositions in the Earth's crust, mountain ranges, deep sea trenches, crust compaction due to glaciers, and so on.

If that sphere were then covered in water, the water would not be the same height everywhere. Instead, the water level would be higher or lower with respect to Earth's center, depending on the integral of the strength of gravity from the center of the Earth to that location. The geoid level coincides with where the water would be. Generally the geoid rises where the earth material is locally more dense, which is where the Earth exerts greater gravitational force and pulls more water from the surrounding area.

Shape

The geoid undulation, geoid height, or geoid anomaly is the height of the geoid relative to a given ellipsoid of reference. The undulation is not standardized, as different countries use different mean sea levels as reference, but most commonly refers to the EGM96 geoid.

Relationship to GPS/GNSS

In maps and common use, the height over the mean sea level (such as orthometric height) is used to indicate the height of elevations while the ellipsoidal height results from the GPS system and similar GNSS.

The deviation ${\displaystyle N}$ between the ellipsoidal height ${\displaystyle h}$ and the orthometric height ${\displaystyle H}$ can be calculated by

$${\displaystyle N=h-H}$$

(An analogous relationship exists between normal heights and the quasigeoid.)

So a GPS receiver on a ship may, during the course of a long voyage, indicate height variations, even though the ship will always be at sea level (neglecting the effects of tides). That is because GPS satellites, orbiting about the center of gravity of the Earth, can measure heights only relative to a geocentric reference ellipsoid. To obtain one's orthometric height, a raw GPS reading must be corrected. Conversely, height determined by spirit leveling from a tide gauge, as in traditional land surveying, is closer to orthometric height. Modern GPS receivers have a grid implemented in their software by which they obtain, from the current position, the height of the geoid (e.g., the EGM-96 geoid) over the World Geodetic System (WGS) ellipsoid. They are then able to correct the height above the WGS ellipsoid to the height above the EGM96 geoid. When height is not zero on a ship, the discrepancy is due to other factors such as ocean tides, atmospheric pressure (meteorological effects), local sea surface topography, and measurement uncertainties.

Relationship to mass density

The surface of the geoid is higher than the reference ellipsoid wherever there is a positive gravity anomaly (mass excess) and lower than the reference ellipsoid wherever there is a negative gravity anomaly (mass deficit).^[5]

This relationship can be understood by recalling that gravity potential is defined so that it has negative values and is inversely proportional to distance from the body. So, while a mass excess will strengthen the gravity acceleration, it will decrease the gravity potential. As a consequence, the geoid's defining equipotential surface will be found displaced away from the mass excess. Analogously, a mass deficit will weaken the gravity pull but will increase the geopotential at a given distance, causing the geoid to move towards the mass deficit. The presence of a localized inclusion in the background medium will rotate the gravity acceleration vectors slightly towards or away from a denser or lighter body, respectively, causing a bump or dimple in the equipotential surface.^[6]

The largest absolute deviation can be found in the Indian Ocean Geoid Low, 106 meters below the average sea level.^[7]

Gravity anomalies

Variations in the height of the geoidal surface are related to anomalous density distributions within the Earth. Geoid measures thus help understanding the internal structure of the planet. Synthetic calculations show that the geoidal signature of a thickened crust (for example, in orogenic belts produced by continental collision) is positive, opposite to what should be expected if the thickening affects the entire lithosphere. Mantle convection also changes the shape of the geoid over time.^[8]

Determination

Calculating the undulation is mathematically challenging.^[9][10] This is why many handheld GPS receivers have built-in undulation lookup tables^[11] to determine the height above sea level.

The precise geoid solution by Vaníček and co-workers improved on the Stokesian approach to geoid computation.^[12] Their solution enables millimetre-to-centimetre accuracy in geoid computation, an order-of-magnitude improvement from previous classical solutions.^{[13][14][15][16]}

Geoid undulations display uncertainties which can be estimated by using several methods, e.g., least-squares collocation (LSC), fuzzy logic, artificial neural networks, radial basis functions (RBF), and geostatistical techniques. Geostatistical approach has been defined as the most-improved technique in prediction of geoid undulation.^[17]

Temporal change

Recent satellite missions, such as the Gravity Field and Steady-State Ocean Circulation Explorer (GOCE) and GRACE, have enabled the study of time-variable geoid signals. The first products based on GOCE satellite data became available online in June 2010, through the European Space Agency (ESA)'s Earth observation user services tools.^[18][19] ESA launched the satellite in March 2009 on a mission to map Earth's gravity with unprecedented accuracy and spatial resolution. On 31 March 2011, the new geoid model was unveiled at the Fourth International GOCE User Workshop hosted at the Technical University of Munich, Germany.^[20] Studies using the time-variable geoid computed from GRACE data have provided information on global hydrologic cycles,^[21] mass balances of ice sheets,^[22] and postglacial rebound.^[23] From postglacial rebound measurements, time-variable GRACE data can be used to deduce the viscosity of Earth's mantle.^[24]

Spherical harmonics representation

Spherical harmonics are often used to approximate the shape of the geoid. The current best such set of spherical harmonic coefficients is EGM2020 (Earth Gravity Model 2020), determined in an international collaborative project led by the National Imagery and Mapping Agency (now the National Geospatial-Intelligence Agency, or NGA). The mathematical description of the non-rotating part of the potential function in this model is:^[25]

$${\displaystyle V={\frac {GM}{r))\left(1+{\sum _{n=2}^{n_{\text{max))))\left({\frac {a}{r))\right)^{n}{\sum _{m=0}^{n)){\overline {P))_{nm}(\sin \phi )\left[{\overline {C))_{nm}\cos m\lambda +{\overline {S))_{nm}\sin m\lambda \right]\right),}$$

where ${\displaystyle \phi \ }$ and ${\displaystyle \lambda \ }$ are geocentric (spherical) latitude and longitude respectively, ${\displaystyle {\overline {P))_{nm))$ are the fully normalized associated Legendre polynomials of degree ${\displaystyle n\ }$ and order ${\displaystyle m\ }$, and ${\displaystyle {\overline {C))_{nm))$ and ${\displaystyle {\overline {S))_{nm))$ are the numerical coefficients of the model based on measured data. The above equation describes the Earth's gravitational potential ${\displaystyle V\ }$, not the geoid itself, at location ${\displaystyle \phi ,\;\lambda ,\;r,\ }$ the co-ordinate ${\displaystyle r\ }$ being the geocentric radius, i.e., distance from the Earth's centre. The geoid is a particular equipotential surface,^[25] and is somewhat involved to compute. The gradient of this potential also provides a model of the gravitational acceleration. The most commonly used EGM96 contains a full set of coefficients to degree and order 360 (i.e., ${\displaystyle n_{\text{max))=360}$), describing details in the global geoid as small as 55 km (or 110 km, depending on the definition of resolution). The number of coefficients, ${\displaystyle {\overline {C))_{nm))$ and ${\displaystyle {\overline {S))_{nm))$, can be determined by first observing in the equation for ${\displaystyle V}$ that for a specific value of ${\displaystyle n}$ there are two coefficients for every value of ${\displaystyle m}$ except for ${\displaystyle m=0}$. There is only one coefficient when ${\displaystyle m=0}$ since ${\displaystyle \sin(0\lambda )=0}$. There are thus ${\displaystyle (2n+1)}$ coefficients for every value of ${\displaystyle n}$. Using these facts and the formula, ${\textstyle \sum _{I=1}^{L}I={\frac {1}{2))L(L+1)}$, it follows that the total number of coefficients is given by

$${\displaystyle \sum _{n=2}^{n_{\text{max))}(2n+1)=n_{\text{max))(n_{\text{max))+1)+n_{\text{max))-3=130317}$$

${\displaystyle n_{\text{max))=360}$

For many applications, the complete series is unnecessarily complex and is truncated after a few (perhaps several dozen) terms.

Still, even higher resolution models have been developed. Many of the authors of EGM96 have published EGM2008. It incorporates much of the new satellite gravity data (e.g., the Gravity Recovery and Climate Experiment), and supports up to degree and order 2160 (1/6 of a degree, requiring over 4 million coefficients),^[26] with additional coefficients extending to degree 2190 and order 2159.^[27] EGM2020 is the planned follow-up of 2020 (now overdue), containing the same number of harmonics generated with better data.^[28]