Geodesy is the science of measuring and representing in temporally varying 3D of the geometry, gravity, and spatial orientation of the Earth and other astronomical bodies, such as planets or circumplanetary systems when it is called planetary geodesy. [1]

Geodynamical phenomena, including crustal motion, tides, and polar motion, can be studied by designing global and national control networks, applying space geodesy and terrestrial geodetic techniques, and relying on datums and coordinate systems. The job titles are geodesist and geodetic surveyor. [2]

A modern instrument for geodetic measurements using satellites

History

Geodesy began in pre-scientific antiquity, so the very word geodesy comes from the Ancient Greek word γεωδαισία or geodaisia (literally, "division of Earth").

Early ideas about the figure of the Earth held the Earth to be flat and the heavens a physical dome spanning over it. Two early arguments for a spherical Earth were that lunar eclipses appear to an observer as circular shadows and that Polaris appears lower and lower in the sky to a traveler headed South.

Definition

In the English-speaking world, geodesy refers to the science of measuring and representing geospatial information, while geomatics encompasses practical applications of geodesy on local and regional scales, including surveying.

Geodesy in the German-speaking world consists of higher geodesy (höhere Geodäsie or Erdmessung, literally "geomensuration") — concerned with measuring Earth on the global scale, and practical geodesy or engineering geodesy (Ingenieurgeodäsie) that includes surveying — measuring parts or regions of Earth.

For the longest time, geodesy was the science of measuring and understanding Earth's geometric shape, orientation in space, and gravitational field; however, geodetic science and operations are applied to other astronomical bodies in our Solar System also. [1]

To a large extent, Earth's shape is the result of rotation, which causes its equatorial bulge, and the competition of geological processes such as the collision of plates, as well as of volcanism, resisted by Earth's gravitational field. This applies to the solid surface, the liquid surface ( dynamic sea surface topography), and Earth's atmosphere. For this reason, the study of Earth's gravitational field is called physical geodesy.

Geoid and reference ellipsoid

Geoid, an approximation for the shape of the Earth; shown here with vertical exaggeration (10000 vertical scaling factor).
Ellipsoid - a mathematical representation of the Earth. When mapping in geodetic coordinates, a latitude circle forms a truncated cone.
Equatorial (a), polar (b) and mean Earth radii as defined in the 1984 World Geodetic System

The geoid essentially is the figure of Earth abstracted from its topographical features. It is an idealized equilibrium surface of seawater, the mean sea level surface in the absence of currents and air pressure variations, and continued under the continental masses. Unlike a reference ellipsoid, the geoid is irregular and too complicated to serve as the computational surface for solving geometrical problems like point positioning. The geometrical separation between the geoid and a reference ellipsoid is called geoidal undulation, and it varies globally between ±110 m based on the GRS 80 ellipsoid.

A reference ellipsoid, customarily chosen to be the same size (volume) as the geoid, is described by its semi-major axis (equatorial radius) a and flattening f. The quantity f = ab/a, where b is the semi-minor axis (polar radius), is purely geometrical. The mechanical ellipticity of Earth (dynamical flattening, symbol J2) can be determined to high precision by observation of satellite orbit perturbations. Its relationship with geometrical flattening is indirect and depends on the internal density distribution or, in simplest terms, the degree of central concentration of mass.

The 1980 Geodetic Reference System ( GRS 80), adopted at the XVII General Assembly of the International Union of Geodesy and Geophysics ( IUGG), posited a 6,378,137 m semi-major axis and a 1:298.257 flattening. GRS 80 essentially constitutes the basis for geodetic positioning by the Global Positioning System (GPS) and is thus also in widespread use outside the geodetic community. Numerous systems used for mapping and charting are becoming obsolete as countries increasingly move to global, geocentric reference systems utilizing the GRS 80 reference ellipsoid.

The geoid is a "realizable" surface, meaning it can be consistently located on Earth by suitable simple measurements from physical objects like a tide gauge. The geoid can, therefore, be considered a physical ("real") surface. The reference ellipsoid, however, has many possible instantiations and is not readily realizable, so it is an abstract surface. The third primary surface of geodetic interest — the topographic surface of Earth — is also realizable.

Coordinate systems in space

The locations of points in 3D space most conveniently are described by three cartesian or rectangular coordinates, X, Y, and Z. Since the advent of satellite positioning, such coordinate systems are typically geocentric, with the Z-axis aligned to Earth's (conventional or instantaneous) rotation axis.

Before the era of satellite geodesy, the coordinate systems associated with a geodetic datum attempted to be geocentric, but with the origin differing from the geocenter by hundreds of meters due to regional deviations in the direction of the plumbline (vertical). These regional geodetic datums, such as ED 50 (European Datum 1950) or NAD 27 (North American Datum 1927), have ellipsoids associated with them that are regional "best fits" to the geoids within their areas of validity, minimizing the deflections of the vertical over these areas.

It is only because GPS satellites orbit about the geocenter that this point becomes naturally the origin of a coordinate system defined by satellite geodetic means, as the satellite positions in space themselves get computed within such a system.

Geocentric coordinate systems used in geodesy can be divided naturally into two classes:

1. The inertial reference systems, where the coordinate axes retain their orientation relative to the fixed stars or, equivalently, to the rotation axes of ideal gyroscopes; the X-axis points to the vernal equinox
2. The co-rotating reference systems (also ECEF or "Earth Centred, Earth Fixed"), in which the axes are "attached" to the solid body of Earth. The X-axis lies within the Greenwich observatory's meridian plane.

The coordinate transformation between these two systems to good approximation is described by (apparent) sidereal time, which accounts for variations in Earth's axial rotation ( length-of-day variations). A more accurate description also accounts for polar motion as a phenomenon closely monitored by geodesists.

Coordinate systems in the plane

2D grid for elliptical coordinates
A Munich archive with lithography plates of maps of Bavaria

In geodetic applications like surveying and mapping, two general types of coordinate systems in the plane are in use:

1. Plano-polar, with points in the plane defined by their distance, s, from a specified point along a ray having a direction α from a baseline or axis;
2. Rectangular, with points defined by distances from two mutually perpendicular axes, x and y. Contrary to the mathematical convention, in geodetic practice, the x-axis points North and the y-axis East.

One can intuitively use rectangular coordinates in the plane for one's current location, in which case the x-axis will point to the local north. More formally, such coordinates can be obtained from 3D coordinates using the artifice of a map projection. It is impossible to map the curved surface of Earth onto a flat map surface without deformation. The compromise most often chosen — called a conformal projection — preserves angles and length ratios so that small circles get mapped as small circles and small squares as squares.

An example of such a projection is UTM ( Universal Transverse Mercator). Within the map plane, we have rectangular coordinates x and y. In this case, the north direction used for reference is the map north, not the local north. The difference between the two is called meridian convergence.

It is easy enough to "translate" between polar and rectangular coordinates in the plane: let, as above, direction and distance be α and s respectively, then we have

{\displaystyle {\begin{aligned}x&=s\cos \alpha \\y&=s\sin \alpha \end{aligned}}}

The reverse transformation is given by:

{\displaystyle {\begin{aligned}s&={\sqrt {x^{2}+y^{2}}}\\\alpha &=\arctan {\frac {y}{x}}.\end{aligned}}}

Heights

Height measurement using satellite altimetry

In geodesy, point or terrain heights are " above sea level" as an irregular, physically defined surface. Height systems in use are:

Each system has its advantages and disadvantages. Both orthometric and normal heights are expressed in metres above sea level, whereas geopotential numbers are measures of potential energy (unit: m2 s−2) and not metric. The reference surface is the geoid, an equigeopotential surface approximating the mean sea level as described above. For normal heights, the reference surface is the so-called quasi-geoid, which has a few-metre separation from the geoid due to the density assumption in its continuation under the continental masses. [3]

One can relate these heights through the geoid undulation concept to ellipsoidal heights (also known as geodetic heights), representing the height of a point above the reference ellipsoid. Satellite positioning receivers typically provide ellipsoidal heights unless fitted with special conversion software based on a model of the geoid.

Geodetic datums

Because coordinates and heights of geodetic points always get obtained within a system that itself was constructed based on real-world observations, geodesists introduced the concept of a "geodetic datum" (plural datums): a physical (real-world) realization of a coordinate system used for describing point locations. This realization follows from choosing (therefore conventional) coordinate values for one or more datum points. In the case of height data, it suffices to choose one datum point — the reference benchmark, typically a tide gauge at the shore. Thus we have vertical datums, such as the NAVD 88 (North American Vertical Datum 1988), NAP ( Normaal Amsterdams Peil), the Kronstadt datum, the Trieste datum, and numerous others.

In both mathematics and geodesy, a coordinate system is a "coordinate system" per ISO terminology, whereas the International Earth Rotation and Reference Systems Service (IERS) uses the term "reference system" for the same. When coordinates are realized by choosing datum points and fixing a geodetic datum, ISO speaks of a "coordinate reference system", whereas IERS uses a "reference frame" for the same. The ISO term for a datum transformation again is a "coordinate transformation". [4]

Positioning

GPS Block IIA satellite orbits over the Earth
Initial acquisition of GPS signal in 2D

Geopositioning, or simply positioning, is the determination of the location, as defined by a set of geodetic coordinates, of a point on land, at sea, or in space within a coordinate system (point positioning) or relative to another point (relative positioning). One computes the position of a point in space from measurements linking terrestrial or extraterrestrial points of known location ("known points") with terrestrial ones of unknown location ("unknown points"). The computation may involve transformations between or among astronomical and terrestrial coordinate systems. Known points used in point positioning can be GNSS satellites or triangulation points of a higher-order network.

Traditionally, geodesists built a hierarchy of networks to allow point positioning within a country. The highest in this hierarchy were triangulation networks, densified into the networks of traverses ( polygons) into which local mapping and surveying measurements, usually collected using a measuring tape, a corner prism, and the red-and-white poles, are tied.

Commonly used nowadays is GPS, except for specialized measurements (e.g., in underground or high-precision engineering). The higher-order networks are measured with static GPS, using differential measurement to determine vectors between terrestrial points. These vectors then get adjusted in a traditional network fashion. A global polyhedron of permanently operating GPS stations under the auspices of the IERS is the basis for defining a single global, geocentric reference frame that serves as the "zero-order" (global) reference to which national measurements are attached.

Real-time kinematic positioning (RTK GPS) is employed frequently in survey mapping. In that measurement technique, unknown points can get quickly tied into nearby terrestrial known points.

One purpose of point positioning is the provision of known points for mapping measurements, also known as (horizontal and vertical) control. There can be thousands of those geodetically determined points in a country, usually documented by national mapping agencies. Surveyors involved in real estate and insurance will use these to tie their local measurements.

Geodetic problems

Geodetic control mark

In geometrical geodesy, there are two main problems:

First (direct or forward) geodetic problem
Given the coordinates of a point and the directional ( azimuth) and distance to a second point, determine the coordinates of that second point.
Second (inverse or reverse) geodetic problem
Given the coordinates of two points, determine the azimuth and length of the (straight, curved, or geodesic) line connecting those points.

The solutions to both problems in plane geometry reduce to simple trigonometry and are valid for small areas on Earth's surface; on a sphere, solutions become significantly more complex as, for example, in the inverse problem, the azimuths differ going between the two end points along the arc of the connecting great circle.

The general solution is called the geodesic for the surface considered, and the differential equations for the geodesic are solvable numerically. On the ellipsoid of revolution, geodesics are expressible in terms of elliptic integrals, which are usually evaluated in terms of a series expansion — see, for example, Vincenty's formulae.

Observational concepts

Axial tilt (or Obliquity), rotation axis, plane of orbit, celestial equator and ecliptic. Earth is shown as viewed from the Sun; the orbit direction is counter-clockwise (to the left).
Global gravity anomaly animation over oceans from the NASA's GRACE (Gravity Recovery And Climate Experiment)

As defined in geodesy (and also astronomy), some basic observational concepts like angles and coordinates include (most commonly from the viewpoint of a local observer):

• Plumbline or vertical: (the line along) the direction of local gravity.
• Zenith: the (direction to the) intersection of the upwards-extending gravity vector at a point and the celestial sphere.
• Nadir: the (direction to the) antipodal point where the downward-extending gravity vector intersects the (obscured) celestial sphere.
• Celestial horizon: a plane perpendicular to the gravity vector at a point.
• Azimuth: the direction angle within the plane of the horizon, typically counted clockwise from the north (in geodesy and astronomy) or the south (in France).
• Elevation: the angular height of an object above the horizon; alternatively: zenith distance equal to 90 degrees minus elevation.
• Local topocentric coordinates: azimuth (direction angle within the plane of the horizon), elevation angle (or zenith angle), distance.
• North celestial pole: the extension of Earth's ( precessing and nutating) instantaneous spin axis extended northward to intersect the celestial sphere. (Similarly for the south celestial pole.)
• Celestial equator: the (instantaneous) intersection of Earth's equatorial plane with the celestial sphere.
• Meridian plane: any plane perpendicular to the celestial equator and containing the celestial poles.
• Local meridian: the plane which contains the direction to the zenith and the celestial pole.

Measurements

Variations in the gravity field of the Moon, from NASA
Gravity measurement devices, pendulum (left) and absolute gravimeter (right)
A relative gravimeter

The reference surface (level) used to determine height differences and height reference systems is known as mean sea level. The traditional spirit level directly produces such (for practical purposes most useful) heights above sea level; the more economical use of GPS instruments for height determination requires precise knowledge of the figure of the geoid, as GPS only gives heights above the GRS80 reference ellipsoid. As geoid determination improves, one may expect that the use of GPS in height determination shall increase, too.

The theodolite is an instrument used to measure horizontal and vertical (relative to the local vertical) angles to target points. In addition, the tachymeter determines, electronically or electro-optically, the distance to a target and is highly automated or even robotic in operations. Widely used for the same purpose is the method of free station position.

Commonly for local detail surveys, tachymeters are employed, although the old-fashioned rectangular technique using an angle prism and steel tape is still an inexpensive alternative. As mentioned, also there are quick and relatively accurate real-time kinematic (RTK) GPS techniques. Data collected are tagged and recorded digitally for entry into Geographic Information System (GIS) databases.

Geodetic GNSS (most commonly GPS) receivers directly produce 3D coordinates in a geocentric coordinate frame. One such frame is WGS84, as well as frames by the International Earth Rotation and Reference Systems Service ( IERS). GNSS receivers have almost completely replaced terrestrial instruments for large-scale base network surveys.

To monitor the Earth's rotation irregularities and plate tectonic motions and for planet-wide geodetic surveys, methods of very-long-baseline interferometry (VLBI) measuring distances to quasars, lunar laser ranging (LLR) measuring distances to prisms on the Moon, and satellite laser ranging (SLR) measuring distances to prisms on artificial satellites, are employed.

Gravity is measured using gravimeters, of which there are two kinds. First are absolute gravimeters, based on measuring the acceleration of free fall (e.g., of a reflecting prism in a vacuum tube). They are used to establish vertical geospatial control or in the field. Second, relative gravimeters are spring-based and more common. They are used in gravity surveys over large areas — to establish the figure of the geoid over these areas. The most accurate relative gravimeters are called superconducting" gravimeters, which are sensitive to one-thousandth of one-billionth of Earth-surface gravity. Twenty-some superconducting gravimeters are used worldwide in studying Earth's tides, rotation, interior, oceanic and atmospheric loading, as well as in verifying the Newtonian constant of gravitation.

In the future, gravity and altitude might become measurable using the special-relativistic concept of time dilation as gauged by optical clocks.

Units and measures on the ellipsoid

The definition of latitude (φ) and longitude (λ) on an ellipsoid of revolution (or spheroid). The graticule spacing is 10 degrees. The latitude is defined as the angle between the normal to the ellipsoid and the equatorial plane.

Geographical latitude and longitude are stated in the units degree, minute of arc, and second of arc. They are angles, not metric measures, and describe the direction of the local normal to the reference ellipsoid of revolution. This direction is approximately the same as the direction of the plumbline, i.e., local gravity, which is also the normal to the geoid surface. For this reason, astronomical position determination – measuring the direction of the plumbline by astronomical means – works reasonably well when one also uses an ellipsoidal model of the figure of the Earth.

One geographical mile, defined as one minute of arc on the equator, equals 1,855.32571922 m. One nautical mile is one minute of astronomical latitude. The radius of curvature of the ellipsoid varies with latitude, being the longest at the pole and the shortest at the equator same as with the nautical mile.

A metre was originally defined as the 10-millionth part of the length from the equator to the North Pole along the meridian through Paris (the target was not quite reached in actual implementation, as it is off by 200 ppm in the current definitions). This situation means that one kilometre roughly equals (1/40,000) * 360 * 60 meridional minutes of arc, or 0.54 nautical miles. (Note that this is not exactly so as the two units had been defined on different bases, so the international nautical mile is 1,852 m exactly, which corresponds to the rounding of 1,000/0.54 m to four digits).

Temporal changes

Global plate tectonic movement using GPS

Various techniques are used in geodesy to study temporally changing surfaces, bodies of mass, physical fields, and dynamical systems. Points on Earth's surface change their location due to a variety of mechanisms:

• Continental plate motion, plate tectonics [5]
• The episodic motion of tectonic origin, especially close to fault lines
• Postglacial land uplift due to isostatic adjustment
• Mass variations due to hydrological changes, including the atmosphere, cryosphere, land hydrology, and oceans
• Sub-daily polar motion [7]
• Length-of-day variability [8]
• Earth's center-of-mass (geocenter) variations [9]
• Anthropogenic movements such as reservoir construction or petroleum or water extraction
A NASA project manager talks about his work for the Space Geodesy Project, including an overview of its four fundamental techniques: GPS, VLBI, LLR/ SLR, and DORIS.

Geodynamics is the discipline that studies deformations and motions of Earth's crust and its solidity as a whole. Often the study of Earth's irregular rotation is included in the above definition. Geodynamical studies require terrestrial reference frames [10] realized by the stations belonging to the Global Geodetic Observing System (GGOS [11]).

Techniques for studying geodynamic phenomena on global scales include:

References

1. ^ a b Vaníček, P.; Krakiwsky, E.J. (1986). Geodesy: the Concepts. New York, USA: Elsevier. p. 45. ISBN  0444-87775-4. Until a decade or two ago, geodesy was thought to occupy the space delimited by the following definition ( Helmert, 1880, p.3): "Geodesy is the science of measuring and portraying the earth's surface." Then people involved with geodesy began to realize that this definition no longer fully reflected the role contemporary geodesy played and started searching for a new framework. This search probably culminated in the new definition of geodesy, accepted by the National Research Council of Canada (NRC), that we quote here (Associate Committee on Geodesy and Geophysics, 1973): Geodesy is the discipline that deals with the measurement and representation of the earth, including its gravity field, in a three-dimensional time varying space. At the 1975 Grenoble meeting of the Commission on Education of the International Association of Geodesy (see §4.2), a virtually identical definition ( Rinner, 1979) was adopted, except for the inclusion of other celestial bodies and their respective gravity fields.
2. ^ "Geodetic Surveyors". Occupational Information Network. 2020-11-26. Retrieved 2022-01-28.
3. ^ Foroughi, Ismael; Tenzer, Robert (2017). "Comparison of different methods for estimating the geoid-to-quasi-geoid separation". Geophysical Journal International. 210 (2): 1001–1020. doi: 10.1093/gji/ggx221. hdl:. ISSN  0956-540X.
4. ^ (ISO 19111: Spatial referencing by coordinates).
5. ^ Altamimi, Zuheir; Métivier, Laurent; Rebischung, Paul; Rouby, Hélène; Collilieux, Xavier (June 2017). "ITRF2014 plate motion model". Geophysical Journal International. 209 (3): 1906–1912. doi: 10.1093/gji/ggx136.
6. ^ Sośnica, Krzysztof; Thaller, Daniela; Dach, Rolf; Jäggi, Adrian; Beutler, Gerhard (August 2013). "Impact of loading displacements on SLR-derived parameters and on the consistency between GNSS and SLR results" (PDF). Journal of Geodesy. 87 (8): 751–769. Bibcode: 2013JGeod..87..751S. doi: 10.1007/s00190-013-0644-1. S2CID  56017067. Archived (PDF) from the original on 2019-04-28.
7. ^ Zajdel, Radosław; Sośnica, Krzysztof; Bury, Grzegorz; Dach, Rolf; Prange, Lars; Kazmierski, Kamil (January 2021). "Sub-daily polar motion from GPS, GLONASS, and Galileo". Journal of Geodesy. 95 (1): 3. Bibcode: 2021JGeod..95....3Z. doi:.
8. ^ Zajdel, Radosław; Sośnica, Krzysztof; Bury, Grzegorz; Dach, Rolf; Prange, Lars (July 2020). "System-specific systematic errors in earth rotation parameters derived from GPS, GLONASS, and Galileo". GPS Solutions. 24 (3): 74. doi:.
9. ^ Zajdel, Radosław; Sośnica, Krzysztof; Bury, Grzegorz (January 2021). "Geocenter coordinates derived from multi-GNSS: a look into the role of solar radiation pressure modeling". GPS Solutions. 25 (1): 1. doi:.
10. ^ Zajdel, R.; Sośnica, K.; Drożdżewski, M.; Bury, G.; Strugarek, D. (November 2019). "Impact of network constraining on the terrestrial reference frame realization based on SLR observations to LAGEOS". Journal of Geodesy. 93 (11): 2293–2313. Bibcode: 2019JGeod..93.2293Z. doi:.
11. ^ Sośnica, Krzysztof; Bosy, Jarosław (2019). "Global Geodetic Observing System 2015–2018". Geodesy and Cartography. doi: 10.24425/gac.2019.126090.
12. ^ Pearlman, M.; Arnold, D.; Davis, M.; Barlier, F.; Biancale, R.; Vasiliev, V.; Ciufolini, I.; Paolozzi, A.; Pavlis, E. C.; Sośnica, K.; Bloßfeld, M. (November 2019). "Laser geodetic satellites: a high-accuracy scientific tool". Journal of Geodesy. 93 (11): 2181–2194. Bibcode: 2019JGeod..93.2181P. doi: 10.1007/s00190-019-01228-y. S2CID  127408940.