Fact-checked by Grok 2 weeks ago

Geodesy

Geodesy is the scientific discipline that measures and represents the Earth's geometric shape, orientation in space, gravity field, and their variations over time, providing a foundational framework for understanding planetary geometry and dynamics.^[1]^[2] The field originated in ancient Greece, where scholars like Eratosthenes estimated the Earth's circumference around 240 BCE using geometric principles and observations of the sun's angle, laying early groundwork for quantifying the planet's size and form.^[3] By the 18th century, expeditions such as those led by Pierre Louis Maupertuis in Lapland and Charles Marie de La Condamine at the equator confirmed the Earth's oblate spheroid shape through arc measurements, resolving debates on whether it was perfectly spherical.^[3] The 20th century marked a shift with the advent of satellite technology, enabling global-scale observations that refined models like the World Geodetic System (WGS 84), which defines an ellipsoid with an equatorial radius of approximately 6,378 kilometers and a flattening ratio of 1/298.257.^[3]^[4] Geodesy encompasses several interconnected branches, including geometric geodesy, which focuses on determining the Earth's surface coordinates and shape using techniques like triangulation and satellite ranging; physical geodesy, which studies the gravity field and the geoid—an irregular equipotential surface approximating mean sea level—to account for mass distribution and deflections; and dynamic geodesy, which examines the Earth's rotation, polar motion, and temporal changes such as crustal deformations.^[3]^[2] These branches rely on precise datums, such as the North American Datum of 1983 (NAD 83), to establish reference frameworks for consistent positioning worldwide.^[3] In modern practice, geodesy employs space-based tools like the Global Positioning System (GPS), which uses a constellation of satellites to achieve centimeter-level accuracy in positioning, alongside methods such as very long baseline interferometry (VLBI) and satellite laser ranging for monitoring geophysical phenomena.^[2]^[3] Applications span navigation and mapping, where it supports the U.S. National Spatial Reference System for boundary delineation and shoreline charting; hazard assessment, including earthquake prediction and volcanic monitoring through millimeter-scale detection of ground motion; and climate studies, tracking sea-level rise, ice mass loss, and post-glacial rebound at rates of several millimeters per year.^[1]^[2]^[4]

Fundamentals

Definition and Scope

Geodesy is the scientific discipline that deals with the measurement and representation of the Earth's geometry, physics, and temporal variations, including its shape, gravity field, and rotation, to determine accurate spatial positions on and near the surface.^[5] This encompasses monitoring changes in these properties over time and space, extending to other celestial bodies as well.^[5] The field integrates observations from ground-based and space-based techniques to model the Earth as a dynamic system.^[1] The scope of geodesy is divided into several key branches that address distinct aspects of Earth's properties and measurements. Physical geodesy focuses on the gravity field and the geoid, which represents the equipotential surface of Earth's gravity.^[6] Geometric geodesy deals with the precise determination of Earth's shape, size, and spatial orientation through coordinate systems.^[6] Dynamic geodesy examines Earth's rotation, polar motion, and tidal effects, capturing temporal variations.^[6] Satellite geodesy, a cross-cutting branch, employs space-based observations to enhance global-scale measurements across these areas.^[7] Geodesy maintains strong interdisciplinary connections with fields such as surveying for precise land mapping, geophysics for studying tectonic movements, oceanography for sea-level monitoring, and navigation for positioning systems.^[8] These links enable applications in creating accurate topographic maps, supporting global navigation satellite systems, and tracking environmental changes like ice mass loss and sea-level rise.^[8] A core concept in modern geodesy is achieving positioning accuracy to sub-millimeter levels, particularly through Global Navigation Satellite Systems (GNSS) in static configurations, which underpin high-precision geodetic networks worldwide.^[9]

Historical Development

The origins of geodesy trace back to ancient civilizations, where early scholars sought to determine the Earth's size and form through geometric measurements. In the 3rd century BCE, Eratosthenes, chief librarian at Alexandria, calculated the Earth's circumference by comparing the angle of the sun's rays at noon in Syene (modern Aswan) and Alexandria, assuming a distance of about 500 miles between the cities; his estimate of roughly 25,000 miles was remarkably close to the modern value of 24,901 miles for the equatorial circumference.^[3] Greek astronomers like Hipparchus further advanced the field in the 2nd century BCE by developing the concepts of latitude and longitude, using stellar observations to establish a grid system for mapping positions on the Earth's surface.^[3] Roman contributions, such as those by Ptolemy in the 2nd century CE, refined these ideas into a geocentric model that influenced cartography for centuries, though it underestimated the Earth's circumference at around 18,000 miles.^[3] During the 17th and 18th centuries, theoretical and empirical advances solidified the understanding of Earth's shape as an oblate spheroid. Isaac Newton proposed in his 1687 Principia Mathematica that gravitational forces would flatten the Earth at the poles due to its rotation, hypothesizing an oblateness that sparked international debate.^[3] To test this, the French Academy of Sciences organized expeditions in the 1730s and 1740s: Pierre Louis Moreau de Maupertuis led a team to Lapland in 1736–1737 to measure a northern arc of the meridian, while Charles Marie de La Condamine's group surveyed a southern arc near the equator in Peru from 1735 to 1744; their findings confirmed a degree of latitude was shorter at the equator than at the poles, validating the oblate spheroid model with a flattening of about 1/300.^[3] The 19th century marked the institutionalization of geodesy through mathematical refinements and national mapping efforts. Carl Friedrich Gauss revolutionized geodetic computations in the early 1800s by developing the least squares adjustment method, which minimized errors in large-scale triangulation networks and became a cornerstone for precise surveying.^[3] In the United States, the Coast Survey was established in 1807 under Ferdinand Hassler, the nation's first scientific agency dedicated to coastal mapping and geodetic control; it introduced rigorous triangulation techniques across North America, laying the groundwork for continental-scale datums.^[10] Alexander Ross Clarke's 1866 ellipsoid, with a semi-major axis of 6,378,206 meters and flattening of 1/295, provided a reference model widely adopted for North American surveys, influencing the North American Datum of 1927.^[3] In the 20th century, technological innovations expanded geodesy's precision and scope. Following World War II, inertial navigation systems (INS), initially developed for military applications like the V-2 rocket, were adapted for geodetic surveying in the 1950s and 1960s; devices such as the Litton INS enabled rapid, autonomous position fixes in remote areas, achieving accuracies of tens of meters over short baselines through gyroscopic and accelerometric measurements.^[11] The launch of Sputnik 1 in 1957 by the Soviet Union initiated the satellite era, as its orbital perturbations provided the first global data on Earth's gravitational field, refining the oblateness estimate to a flattening reciprocal of 298.^[12] The U.S. Navy's Transit system, operational from the early 1960s, used Doppler shifts from polar-orbiting satellites to determine positions with 20-meter accuracy, supporting the World Geodetic System 1960 (WGS 60) datum.^[12] The 1970s and 1980s saw the rise of the Global Positioning System (GPS), transforming geodesy into a global, real-time discipline. Conceived in the 1970s by the U.S. Department of Defense, GPS built on Transit technology with a constellation of medium-Earth-orbit satellites broadcasting precise timing signals; initial development included launches of prototype Block I satellites starting in 1978, achieving full operational capability by 1995 but enabling geodetic applications from the mid-1980s with sub-meter static positioning via carrier-phase measurements.^[13] This shift from ground-based to space-based methods allowed for centimeter-level accuracy in tectonic studies and datum realizations, culminating in the WGS 84 reference frame.^[13] Into the late 20th and early 21st centuries, satellite missions advanced gravity field modeling. The Gravity Recovery and Climate Experiment (GRACE), launched in 2002 by NASA and the German Aerospace Center, used twin satellites in a low-Earth orbit to measure inter-satellite distance variations caused by gravity anomalies, providing monthly global maps of Earth's mass distribution with resolutions down to 300–400 km and sensitivities to changes as small as 1 cm of water equivalent.^[14] Its successor, GRACE Follow-On (GRACE-FO), launched in 2018, continued these observations with enhanced laser interferometry for inter-satellite ranging, maintaining data continuity through 2025. By 2025, recent milestones include the integration of artificial intelligence (AI) in geodetic data processing, where machine learning algorithms enhance the quality of observations from GNSS and gravimetry by automating outlier detection and noise reduction, as pursued by the International Association of Geodesy (IAG)'s AI4G initiative to improve product accuracy for Earth system monitoring.^[15] Advancements in quantum gravimetry, leveraging cold-atom interferometry, have enabled portable sensors achieving 10^{-8} m/s² precision for field measurements, supporting next-generation satellite missions like proposed quantum-enhanced gradiometers for finer resolution of geoid undulations.^[16]^[17]

Earth's Geometric and Gravitational Models

Geoid

The geoid is defined as the equipotential surface of the Earth's gravity field that best approximates mean sea level in the absence of currents, tides, and atmospheric effects, extending continuously under landmasses as a hypothetical surface of constant gravitational potential.^[18] This surface undulates irregularly due to uneven mass distribution within the Earth, such as variations in crustal density and mantle convection, resulting in deviations from a smooth sphere or ellipsoid.^[19] Physically, the geoid arises from the total gravity potential W, which combines the gravitational potential V generated by the Earth's mass (external to the surface) and the centrifugal potential \Phi due to Earth's rotation. On the geoid, W equals a constant value W_0, typically corresponding to global mean sea level.^[18] The relationship between the geoid and heights is described by Bruns' formula, which expresses the orthometric height H above the geoid at a point P as

H = \frac{W_0 - W_P}{\gamma},

where W_P is the gravity potential at P and \gamma is the normal gravity (approximating actual gravity on the geoid).^[20] This formula underscores the geoid's role as the zero-level reference for orthometric heights, which measure elevation relative to sea level along plumb lines.^[21] Geoid determination involves several methods, each leveraging different observations of the gravity field. Gravimetric methods compute geoid undulations using gravity anomalies through Stokes' integral, which integrates the anomalous gravity over the Earth's surface after removing long-wavelength effects from global models and restoring topographic contributions via the remove-compute-restore technique.^[22] Astrogeodetic methods derive local geoid heights from astronomic observations of vertical deflections, combining them with gravity data for high-resolution regional models, particularly in areas with sparse satellite coverage like mountainous terrain.^[23] Satellite-based approaches, such as those from the Gravity Recovery and Climate Experiment (GRACE), its successor GRACE Follow-On (GRACE-FO, ongoing as of 2025), and Gravity Field and Steady-State Ocean Circulation Explorer (GOCE) missions, provide global gravity field data that inform high-degree spherical harmonic models like the Earth Gravitational Model 2008 (EGM2008), achieving resolutions down to about 10 km.^[24]^[25]^[26]^[27] Geoid undulations, or the vertical separations from a reference ellipsoid, typically range from -100 to +100 meters globally, reflecting imbalances in the gravity field.^[28] These undulations can be decomposed into components due to anomalous mass distributions (gravity effects) and topographic loading, with the latter often corrected using digital elevation models to isolate the true geopotential signal.^[19] The geoid serves as the fundamental reference for orthometric heights, enabling the conversion of satellite-derived ellipsoidal heights to physically meaningful elevations above sea level, which is essential for hydrology, engineering, and environmental monitoring.^[29] It also plays a critical role in unifying global geodetic datums by providing a consistent vertical framework that bridges regional height systems and facilitates accurate realization of international reference frames.^[30]

Reference Ellipsoid

The reference ellipsoid is a mathematical model of an oblate spheroid that approximates the Earth's overall shape for geodetic purposes, defined by its semi-major axis a (the equatorial radius) and flattening f = (a - b)/a, where b is the semi-minor axis (the polar radius).^[31] This model simplifies computations by representing the Earth as a smooth, rotationally symmetric surface, with typical parameters including a \approx 6378 km and b \approx 6357 km to reflect the planet's equatorial bulge due to rotation.^[32] The dynamical form factor, such as the second-degree gravitational harmonic J_2, further refines the model by incorporating rotational effects, with J_2 = 1.08263 \times 10^{-3} in the Geodetic Reference System 1980 (GRS80).^[33] The surface of the reference ellipsoid is mathematically described by the equation

\frac{x^2 + y^2}{a^2} + \frac{z^2}{b^2} = 1,

where the z-axis aligns with the Earth's rotation axis, providing a geometric foundation for positioning calculations.^[34] Historical reference ellipsoids were developed based on arc measurements and gravitational data to fit regional or global observations. The Bessel ellipsoid of 1841, with a = 6377397.155 m and $1/f = 299.1528128, was widely used for European and Asian networks due to its fit to continental measurements.^[35] The Clarke ellipsoids, proposed by Alexander Ross Clarke, include the 1866 version (a = 6378206.4 m, $1/f = 294.97869821) optimized for North American surveys, and the 1880 version (a = 6378249.2 m, $1/f = 293) for global applications.^[36]^[37] A modern standard, the WGS84 ellipsoid (adopted in 1984), uses a = 6378137 m and f = 1/298.257223563, derived from satellite and gravitational data for worldwide consistency.^[38] As the foundational geometric surface in geodesy, the reference ellipsoid serves as the basis for map projections, which transform its curved coordinates into flat representations, and for coordinate transformations between different systems.^[39] It defines ellipsoidal heights measured perpendicular to this surface, which differ from orthometric heights tied to the geoid by the ellipsoidal-geoid separation.^[40] Contemporary refinements to reference ellipsoids, such as those in the International Terrestrial Reference Frame (ITRF), build on the GRS80 parameters and integrate Global Positioning System (GPS) observations to achieve sub-centimeter accuracy in realizations like ITRF2020, enabling precise global positioning.^[41]

Reference Systems and Datums

Coordinate Systems

In geodesy, coordinate systems provide the mathematical framework for specifying positions on or near the Earth's surface in three dimensions. The two primary types are Cartesian and geodetic coordinates, both defined relative to a reference ellipsoid that approximates the Earth's shape. Cartesian coordinates, often in the Earth-Centered, Earth-Fixed (ECEF) system, use orthogonal axes (X, Y, Z) with the origin at the Earth's center of mass; the Z-axis aligns with the Earth's rotation axis (positive toward the North Pole), the X-axis points toward the prime meridian at the equator, and the Y-axis completes the right-handed system.^[42] Geodetic coordinates, in contrast, consist of latitude φ (angle from the equator to the point along the ellipsoid normal), longitude λ (angle from the prime meridian in the equatorial plane), and ellipsoidal height h (distance along the normal from the ellipsoid surface).^[43] Transformations between geodetic and Cartesian coordinates are essential for computations in positioning and navigation. The conversion from geodetic (φ, λ, h) to ECEF Cartesian (X, Y, Z) coordinates is given by:

\begin{align*} X &= (N + h) \cos \phi \cos \lambda, \\ Y &= (N + h) \cos \phi \sin \lambda, \\ Z &= \left( N (1 - e^2) + h \right) \sin \phi, \end{align*}

where N = \frac{a}{\sqrt{1 - e^2 \sin^2 \phi}} is the prime vertical radius of curvature, a is the semi-major axis of the ellipsoid, and e is the eccentricity.^[42] These equations account for the ellipsoidal geometry, ensuring accurate representation of positions up to the precision required for global applications. The inverse transformation from Cartesian to geodetic coordinates is more complex and typically solved iteratively due to the nonlinearity introduced by the ellipsoid.^[44] For applications involving dynamics, such as satellite orbits, celestial coordinate systems like the Earth-Centered Inertial (ECI) frame are used. The ECI system has its origin at the Earth's center but is fixed relative to distant stars, with axes aligned to the mean equator and equinox at a reference epoch (e.g., J2000.0); unlike ECEF, it does not rotate with the Earth.^[45] Transformations between ECI and ECEF incorporate Earth's rotation, including sidereal time, as well as long-term effects like precession (slow axial wobble over 26,000 years) and nutation (smaller oscillations with periods up to 18.6 years due to lunar and solar torques).^[45] These adjustments ensure inertial stability for space-based geodesy. In regional or local contexts, where curvature effects are minimal, planar approximations simplify computations by projecting the ellipsoidal surface onto a flat plane. The Universal Transverse Mercator (UTM) system divides the Earth into 60 zones, each 6° wide in longitude, using the Transverse Mercator projection to minimize distortion (scale error less than 0.1% within each zone, spanning approximately 670 km at the equator). Coordinates in UTM are expressed as easting (X) and northing (Y) in meters, facilitating map-based work; similar systems, like state plane coordinates in the U.S., adapt projections for even smaller areas to further reduce errors. The World Geodetic System 1984 (WGS84) serves as the global standard for coordinate systems, particularly in GPS applications, defining an ECEF frame aligned to the Earth's center of mass with parameters such as semi-major axis a = 6,378,137 m and flattening f = 1/298.257223563.^[38] A key distinction within WGS84 is between geocentric latitude (angle from the equator to the position vector from the Earth's center) and geodetic latitude (angle along the ellipsoid normal), which differ by up to 11.5 minutes of arc at mid-latitudes due to the Earth's oblateness, converging to equality at the equator and poles.^[46] Coordinate transformations between different systems or datums often employ similarity models to account for discrepancies in origin, orientation, and scale. The Helmert (7-parameter) transformation, widely used for datum shifts, applies three translations (ΔX, ΔY, ΔZ), three rotations (small angles Rx, Ry, Rz), and a scale factor (1 + s) to align frames, modeled as:

\begin{pmatrix} X' \\ Y' \\ Z' \end{pmatrix} = (1 + s) R \begin{pmatrix} X \\ Y \\ Z \end{pmatrix} + \begin{pmatrix} \Delta X \\ \Delta Y \\ \Delta Z \end{pmatrix},

where R is the rotation matrix.^[44] The Bursa-Wolf model, a variant of this approach, applies the scale and rotations to position vectors before translation, offering equivalent results for small parameters but preferred in some standards for its formulation in geocentric coordinates. These methods achieve sub-meter accuracy in global transformations when parameters are derived from networks of control points.^[44]

Geodetic Datums

A geodetic datum serves as a reference framework that defines the relationship between a coordinate system and the physical Earth, primarily through a set of parameters specifying the position, orientation, and scale of a reference ellipsoid relative to the Earth's surface. These parameters typically include the origin (often defined by the coordinates of an initial point), scale (related to the ellipsoid's dimensions), and orientation (determined by azimuths or rotation angles), enabling the accurate mapping of positions on or near the Earth. For horizontal datums, which focus on latitude and longitude, the framework is anchored to the ellipsoid's surface, while three-dimensional datums extend this with six geocentric parameters: three translations for the origin, three rotations for orientation, and one scale factor.^[47]^[32]^[48] Horizontal geodetic datums provide the foundation for measuring positions across the Earth's surface, with notable examples including the North American Datum of 1927 (NAD27) and the North American Datum of 1983 (NAD83). NAD27, developed using the Clarke 1866 ellipsoid, relied on a network of about 26,000 control points primarily from triangulation and astronomic observations, with its origin at Meades Ranch, Kansas, and an orientation tied to that location's azimuth; however, its ellipsoid center is offset by approximately 236 meters from modern global models due to inconsistencies in the North American network. In contrast, NAD83 employs the Geodetic Reference System 1980 (GRS80) ellipsoid, which better approximates the global Earth figure, and is realized through over 250,000 control points incorporating satellite data for improved accuracy and plate-fixed alignment, while the International Terrestrial Reference Frame (ITRF) extends this globally with a geocentric origin at the Earth's center of mass. These shifts from regional to global realizations reflect advancements in precision, reducing distortions in continental-scale mapping.^[47]^[49]^[50] As of 2025, the U.S. National Spatial Reference System (NSRS), which includes NAD83, is undergoing modernization by the National Geodetic Survey (NGS). Preliminary beta products for the new North American Terrestrial Reference Frame of 2022 (NATRF2022), replacing NAD83, were released on June 17, 2025, incorporating time-dependent models for tectonic plate motion and GNSS data for enhanced accuracy. The full operational release is expected in 2026.^[51] The realization of geodetic datums occurs through dense networks of control points whose positions are precisely determined and tied to the defined parameters, ensuring consistency in coordinate assignments. Historically, these networks were established via ground-based methods like triangulation, but modern realizations leverage Global Navigation Satellite Systems (GNSS) such as GPS for global coverage and sub-centimeter accuracy, allowing datums to be maintained with high fidelity across continents. This GNSS-based approach facilitates the integration of millions of points into a unified framework, minimizing local distortions and supporting applications from surveying to navigation.^[47]^[52]^[53] Datum shifts arise from discrepancies between different reference frameworks, often driven by geophysical processes like plate tectonics, which cause continental drift at rates of several centimeters per year, necessitating transformations that account for translations, rotations, and scaling. For instance, converting coordinates from NAD27 to NAD83 involves systematic offsets up to tens of meters due to the ellipsoids' differing centers and the inclusion of new observations. Transformation grids, such as the NADCON (North American Datum Conversion) tool developed by the National Geodetic Survey, model these shifts using interpolated grids derived from over 150,000 control points, enabling accurate conversions for U.S. applications without requiring full reobservation.^[54]^[55] International standards for geodetic datums are embodied in the International Terrestrial Reference System (ITRS), a conceptual framework maintained by the International Earth Rotation and Reference Systems Service (IERS), with realizations like the ITRF updated approximately every five to six years to incorporate new observations and refine parameters. The ITRS employs a no-net-rotation condition, minimizing cumulative rotation relative to the Earth's crustal plates by constraining periodic motions in a core network of stations, thus preserving the datum's stability amid tectonic movements. These updates ensure global consistency, with ITRF2020, for example, integrating GNSS, satellite laser ranging, and very long baseline interferometry data for millimeter-level precision.^[56]^[57]^[58] Horizontal geodetic datums integrate with vertical datums, which are gravity-based references like mean sea level, to form complete three-dimensional positioning systems, though the primary focus remains on anchoring horizontal coordinates to the ellipsoid for planar measurements.^[47]

Vertical Positioning

Heights

In geodesy, heights represent the vertical position of points relative to specific reference surfaces, with three primary types: ellipsoidal, orthometric, and normal heights. The ellipsoidal height h is the geometric distance measured along the normal to a reference ellipsoid from its surface to the point of interest, typically determined using satellite-based methods like GNSS.^[59]^[29] Orthometric height H is the physical distance along the plumb line (direction of gravity) from the geoid—an equipotential surface approximating mean sea level—to the point, providing a reference meaningful for gravitational potential.^[59]^[29] Normal height H^* measures the distance from the quasi-geoid, a surface that closely follows the geoid but is defined using normal gravity to avoid irregularities in actual gravity measurements.^[59] The relationships between these heights account for deviations between the reference surfaces. The fundamental relation is h = H + N, where N is the geoid undulation, the separation between the ellipsoid and geoid along the plumb line.^[59]^[29] For normal heights, h \approx H^* + \zeta, with \zeta as the height anomaly representing the quasi-geoid to ellipsoid separation; differences between orthometric and normal heights are typically on the order of millimeters to a few centimeters, even in mountainous regions.^[59] Orthometric heights incorporate gravity corrections using actual observed gravity g along the plumb line, unlike normal heights which use a theoretical normal gravity \gamma.^[29] Heights are measured through geometric and gravimetric principles. Geometric methods, such as trigonometric leveling, involve observing angles and distances to compute height differences relative to a reference surface, providing direct ellipsoidal or approximate orthometric values.^[59] Gravimetric methods determine heights via differences in gravitational potential, integrating gravity observations to model separations like N or \zeta.^[29] The orthometric height difference \Delta H between two points is defined as the integral of gravity along the plumb line path, approximated as

\Delta H = \int \frac{g \, ds}{\gamma},

where g is the actual gravity, ds is the element of distance along the path, and \gamma is the mean normal gravity; this accounts for the geopotential number and ensures consistency with equipotential references.^[29] Challenges in height determination arise from environmental variations and gravitational irregularities. Tidal effects cause fluctuations in sea level and the geoid, introducing biases up to decimeters in local systems.^[29] Deflections of the vertical—the angular difference between the plumb line and ellipsoid normal—can reach up to about 0.03° (100 arcseconds) in extreme cases, such as in the Himalayas, due to mass anomalies, complicating alignments between geometric and physical heights.^[60] These heights play critical roles in hydrology for modeling water flow along equipotential surfaces and in aviation for accurate altimetry relative to mean sea level.^[61] Modern techniques achieve sub-centimeter accuracy by combining GNSS for ellipsoidal heights with global gravity models like EGM2008, enabling precise orthometric determinations without extensive leveling networks.^[62]^[23]

Vertical Datums

A vertical datum serves as a standardized reference surface for measuring heights, typically defined relative to mean sea level (MSL) or an equipotential surface of the Earth's gravity field, enabling consistent determination of elevations above or below this reference.^[63] In geodesy, it is realized through a network of benchmarks with assigned heights, often tied to the geoid, which approximates MSL as an equipotential surface where gravity potential is constant.^[64] This reference is essential for applications like surveying, mapping, and engineering, as it accounts for the irregular gravity field influencing vertical positions.^[65] Globally, the International Association of Geodesy (IAG) promotes a unified vertical reference system based on a conventional zero potential surface with a fixed gravity potential value, W₀ ≈ 62636853.4 m²/s², defining the ideal geoid for worldwide height consistency.^[66] The European Vertical Reference System (EVRS), realized through the Unified European Levelling Network (UELN), exemplifies this approach as a zero-tidal system that integrates leveling data, gravity observations, and GNSS-derived heights across Europe, with EVRF2007 providing normal heights relative to this potential surface.^[67] This framework supports kinematic height references aligned with the International Terrestrial Reference System (ITRS), facilitating continent-wide vertical control.^[68] Regionally, vertical datums are tailored to specific areas, often building on historical leveling networks. In the United States, the North American Vertical Datum of 1988 (NAVD88) was established through a continent-wide adjustment of leveling surveys originating from the National Geodetic Vertical Datum of 1929 (NGVD29), with heights held fixed at one benchmark while incorporating data from 26 tide gauges to approximate MSL.^[69] Similarly, the Australian Height Datum (AHD), adopted in 1971, defines zero height as the MSL determined from 30 tide gauges around the mainland, serving as the national reference for vertical control despite known distortions from long-wavelength geoid undulations.^[70] Realization of vertical datums traditionally relies on tide gauges to establish local MSL, with networks of benchmarks connected via differential leveling to propagate heights inland.^[50] Modern hybrid systems increasingly incorporate Global Navigation Satellite Systems (GNSS) for ellipsoidal heights and gravimetry to model the geoid, enabling orthometric heights via the relation H = h - N, where H is orthometric height, h is ellipsoidal height, and N is geoid undulation; this approach, as in the U.S. GRAV-D project, aims for centimeter-level accuracy without extensive re-leveling.^[71] Challenges in maintaining vertical datums include the impacts of global sea level rise, currently averaging 3.4 mm per year since 1993, which progressively misaligns tide gauge-based references and necessitates periodic updates to preserve accuracy.^[72] For instance, NAVD88 exhibits biases up to 1 meter relative to modern geoid models due to crustal motions and sea level changes, prompting the National Geodetic Survey to develop the North American-Pacific Geopotential Datum of 2022 (NAPGD2022) as a gravity-based replacement incorporating dynamic components for improved temporal stability. As of November 2025, NAPGD2022 remains in beta testing phase, with full implementation expected soon.^[73] Integration of vertical datums with horizontal ones occurs through three-dimensional networks, where GNSS provides unified positioning, and gravity anomalies ensure consistency between the geoid and reference ellipsoid by correcting for mass distribution effects. This holistic approach, as in NAPGD2022 paired with NATRF2022, aligns vertical and horizontal components in a single reference frame, enhancing overall geodetic reliability for applications like coastal management.^[69]

Positioning and Measurement Techniques

Global Positioning Methods

Global Navigation Satellite Systems (GNSS) form the cornerstone of modern global positioning in geodesy, enabling precise determination of positions anywhere on Earth through satellite signals. The primary GNSS constellations include the Global Positioning System (GPS) operated by the United States, GLONASS by Russia, Galileo by the European Union, and BeiDou by China, each comprising dozens of satellites in medium Earth orbit that broadcast ranging signals for user receivers.^[74]^[75] These systems operate on the principle of trilateration, where a receiver measures distances to at least four satellites to solve for its three-dimensional position and receiver clock offset, with carrier-phase measurements providing the highest precision by tracking the phase of the satellite's radio signal wavelength.^[76]^[77] The fundamental measurement in GNSS is the pseudorange, which approximates the distance between satellite and receiver but includes biases. The pseudorange \rho is given by \rho = c \Delta t = \| \mathbf{r}^s - \mathbf{r}^u \| + c(\delta t^u - \delta t^s) + d_{ion} + d_{trop} + \epsilon, where c is the speed of light, \Delta t is the signal travel time, \mathbf{r}^s and \mathbf{r}^u are satellite and user positions, \delta t^u and \delta t^s are clock offsets, d_{ion} and d_{trop} represent ionospheric and tropospheric delays, and \epsilon includes other errors like multipath.^[78] Positioning accuracy varies by technique: Real-Time Kinematic (RTK) achieves centimeter-level precision using nearby reference stations to resolve carrier-phase ambiguities, while Precise Point Positioning (PPP) delivers decimeter-level accuracy globally by applying precise satellite orbits and clock corrections without local references.^[79]^[80] Key error sources include ionospheric delays, which can reach meters during solar maxima and are mitigated via dual-frequency observations, and multipath reflections from surfaces that distort signals by up to tens of centimeters.^[80]^[81] Beyond GNSS, space-oriented techniques like Very Long Baseline Interferometry (VLBI) and Satellite Laser Ranging (SLR) contribute to global positioning by providing independent validations and ties to reference frames. VLBI uses a network of radio telescopes to observe distant quasars, measuring baseline vectors between stations with millimeter accuracy to determine Earth orientation parameters such as polar motion and UT1, essential for linking terrestrial and celestial frames.^[82]^[83] SLR employs ground-based lasers to range to retroreflector-equipped satellites like LAGEOS, yielding station coordinates with sub-centimeter precision and supporting orbit determinations that calibrate GNSS.^[84]^[85] These methods underpin critical geodetic applications, including maintenance of the International Terrestrial Reference Frame (ITRF), where GNSS and SLR data combine to realize a stable global coordinate system with velocities tied to plate motions.^[86] In plate tectonics monitoring, continuous GNSS networks detect crustal deformations at rates of 2-10 cm per year, as seen in subduction zones and rift areas, enabling forecasts of seismic hazards.^[87]^[88] Recent advances through 2025 emphasize multi-GNSS integration, combining signals from GPS, GLONASS, Galileo, and BeiDou to improve satellite visibility and reduce dilution of precision, achieving up to 20% better convergence times in PPP compared to single-system use.^[89] Additionally, proposed Low Earth Orbit (LEO) constellations, such as Starlink, with thousands of satellites at altitudes below 600 km, could enhance global coverage for GNSS augmentation by providing denser signals for real-time kinematic positioning, potentially enabling ubiquitous centimeter accuracy even in challenging environments.^[90]

Local Measurement Methods

Local measurement methods in geodesy focus on terrestrial techniques for establishing high-precision control networks over regional areas, typically spanning tens to hundreds of kilometers, where ground-based observations provide the necessary accuracy for mapping and infrastructure projects. These methods emphasize direct measurements of angles, distances, and heights using specialized instruments, forming the foundation for local coordinate systems before integration with broader reference frameworks. Traditional approaches include triangulation, trilateration, and leveling, which have evolved with technological advancements to enhance efficiency and precision.^[91] Triangulation involves measuring angles between points to form a network of triangles, starting from a known baseline distance. Theodolites are employed to record horizontal and vertical angles with accuracies up to 1 arcsecond, enabling position determination via the law of sines and enabling network extension across landscapes. This method was historically dominant for large-scale control surveys due to its reliance on angular precision rather than extensive distance measurements.^[91]^[92] Trilateration complements triangulation by focusing solely on distance measurements, using electronic distance measurement (EDM) instruments to compute positions from intersections of spheres centered at known points. EDM devices achieve accuracies of 1:20,000 over distances up to 5 km, making trilateration simpler and more cost-effective for dense networks where visibility for angles is limited.^[91]^[93] Leveling determines height differences through differential observations with spirit levels and rods, establishing vertical control relative to a datum like mean sea level. Backsight and foresight readings calculate height of instrument and elevation changes, with circuits closed to verify consistency. This technique is essential for orthometric heights in local surveys.^[94] Modern instrumentation has streamlined these processes. Total stations integrate theodolites with EDM and electronic data recording, measuring angles and distances simultaneously up to several kilometers with millimeter-level precision, depending on the instrument and conditions, facilitating automated computations for traversing and detailed surveys.^[95] For deformation analysis, Interferometric Synthetic Aperture Radar (InSAR) uses radar phase differences from satellite imagery to map surface displacements at millimeter resolution over local areas, particularly effective for monitoring subsidence in urban environments. Gravimeters, such as absolute models like the FG5 (accuracy ±2 μGal), conduct anomaly surveys to detect gravity variations tied to subsurface mass changes, supporting geoid refinements in regional networks.^[96]^[97] Error sources in these methods include atmospheric refraction, which bends light rays in leveling sights, introducing elevation discrepancies up to several millimeters over 60 m, and Earth's curvature, which causes systematic drops in line-of-sight observations. Refraction is mitigated by balancing sight lengths or applying corrections based on temperature models, while curvature effects are isolated and adjusted in unbalanced setups. Least squares adjustment optimizes network solutions by minimizing residuals in overdetermined systems, accounting for errors in both observations and design matrices via weighted total least squares for enhanced reliability in coordinate transformations.^[98]^[99] A key correction in leveling accounts for curvature:

\delta = \frac{d^2}{2R}

where \delta is the correction (in meters), d is the sight distance (in meters), and R = 6371 km is Earth's mean radius; for d = 60 m, \delta \approx 0.00028 m.^[98] These methods apply to cadastral surveying, where precise boundary delineations support land registration, and subsidence monitoring, detecting urban rates of 1-10 mm/yr, such as 8.8 mm/yr in the Anambra Metropolis linked to groundwater extraction.^[100] Integration with Global Navigation Satellite Systems (GNSS) forms hybrid networks, combining local terrestrial data with satellite positions for sub-centimeter control, while post-2020 LiDAR advancements via terrestrial scanners produce dense 3D point clouds for deformation modeling, refined by B-splines to handle large datasets in high-resolution mapping.^[101]

Geodetic Computations

Geodetic Problems

Geodetic problems encompass the fundamental computational tasks in geodesy for determining positions, distances, and orientations on the Earth's ellipsoidal surface. These problems arise when transforming measurements between points using coordinate systems such as latitude and longitude, focusing on geodesics—the shortest paths on the ellipsoid. The two primary problems are the direct and inverse geodesic calculations, which form the basis for more complex network adjustments and positioning. The direct geodetic problem involves computing the coordinates of an endpoint given the starting position (latitude φ₁ and longitude λ₁), initial azimuth α₁, and geodesic distance s. This requires solving for the endpoint (φ₂, λ₂) and final azimuth α₂, accounting for the ellipsoid's curvature through series expansions or iterative methods. Vincenty (1975) provided a compact iterative solution using nested equations for efficient programming, applicable to distances from centimeters to over 20,000 km with a maximum error of 0.01 mm in distance and 0.00005″ in azimuth.^[102] The inverse geodetic problem, conversely, determines the geodesic distance s and azimuths α₁ and α₂ given two points (φ₁, λ₁) and (φ₂, λ₂). Vincenty's method employs an iterative approach on an auxiliary sphere to refine the longitude difference λ until convergence on the angular distance σ, starting with an initial guess λ ≈ L (where L = λ₂ - λ₁), using reduced latitudes u₁ = \atan((1-f) \tan \phi_1) and u₂. The central angular distance is computed as:

\sin^2 \sigma = (\cos u_1 \sin \lambda)^2 + (\cos u_1 \sin u_2 - \sin u_1 \cos u_2 \cos \lambda)^2

\cos \sigma = \sin u_1 \sin u_2 + \cos u_1 \cos u_2 \cos \lambda

with iteration on λ via ellipsoidal corrections involving eccentricity until change is negligible. Once converged, the auxiliary angular distance Δσ is found as:

\Delta \sigma = B \sin \sigma \left\{ \cos 2\sigma_m + B \left[ \cos \sigma (-1 + 2 \cos^2 2\sigma_m) - \frac{1}{6} B \cos 2\sigma_m (-3 + 4 \sin^2 \sigma) (-3 + 4 \cos^2 2\sigma_m) \right] \right\}

where σ_m = (σ₁ + σ)/2, B = 1/1024 [256 + U² (-128 + U² (74 - 47 U²))], and U² = \cos^2 \alpha (a² - b²)/b² with α the equatorial azimuth. The geodesic distance is then s = b (σ - Δσ), with azimuths from:

\tan \alpha_1 = \frac{\cos u_2 \sin \lambda}{\cos u_1 \sin u_2 - \sin u_1 \cos u_2 \cos \lambda}

\tan \alpha_2 = \frac{\cos u_1 \sin \lambda}{\sin u_1 \cos u_2 + \cos u_1 \sin u_2 \cos \lambda}

This yields sub-millimeter accuracy for most global distances.^[102] Other geodetic problems include intersection and resection. The intersection problem computes the position where two geodesics, defined by known points and azimuths, meet on the ellipsoid; it is solved iteratively using Newton-Raphson methods combined with numerical integration for the eccentricity correction, without arc length restrictions. Resection determines an unknown point's position from measured distances or angles to known points, often requiring similar iterative adjustments on the ellipsoid to handle curvature effects. Sjöberg (2002) extended these to constant-azimuth arcs, normal sections, and full geodesics.^[103] Map projection distortions represent another class of geodetic problems, involving the quantification of errors in transforming ellipsoidal coordinates to a plane, affecting scale, shape, area, and direction. These distortions are analyzed using metrics like Tissot's indicatrix, which maps infinitesimal circles to ellipses on the projection; conformal projections preserve angles but distort areas, while equal-area projections maintain sizes at the expense of shapes. Snyder (1987) details computational methods for over 20 projections, including formulas for scale factors k = ds'/ds (where ds' is projected and ds is geodesic distance), essential for minimizing errors in large-scale mapping.^[104] Modern numerical methods enhance solution accuracy and robustness. Karney's (2013) algorithm refines the direct and inverse problems using sixth-order series in flattening f, Newton's method for inversion, and auxiliary sphere mappings, achieving round-off errors under 15 nanometers on the WGS84 ellipsoid—far surpassing Vincenty's 0.5 mm threshold for practical applications. It handles all point pairs, including near-antipodal cases where Vincenty's iteration may fail.^[105] Challenges in geodetic computations include polar singularities, where azimuths become undefined near the poles, requiring special handling like auxiliary coordinate rotations to avoid division by zero in longitude-based formulas. Karney (2013) mitigates this through robust convergence criteria independent of longitude. In large-scale networks, collinearity—when points align nearly linearly—leads to ill-conditioned matrices in least-squares adjustments, amplifying errors; this is addressed by regularization techniques or additional constraints to ensure network stability.^[105]

Units and Measures on the Ellipsoid

In geodesy, measurements on the ellipsoidal surface of the Earth require specific units and conventions that account for the oblate spheroid geometry, distinguishing them from spherical or planar approximations. These include angular measures for latitude and longitude, linear distances along curved paths, radii of curvature for local approximations, and area computations that preserve ellipsoidal properties. The reference ellipsoid, such as the WGS 84 model, serves as the baseline for these calculations, with the semi-major axis defined precisely in meters to ensure global consistency.^[106] Angular units in ellipsoidal geodesy primarily use degrees for geodetic latitude (φ) and longitude (λ), though radians are employed in computational formulas for precision. Latitude φ is the angle between the equatorial plane and the normal to the ellipsoid at a point, ranging from -90° to 90°. For certain computations, such as geodesic projections or auxiliary sphere mappings, the reduced latitude (β), also known as parametric latitude, is used; it is defined by tan β = \sqrt{1 - e^2} \tan \phi, where e is the eccentricity of the ellipsoid, simplifying trigonometric relations by treating the ellipsoid as a sphere of stretched meridians.^[31]^[107] Linear distances on the ellipsoid are measured along geodesics, which are the shortest paths analogous to great circles on a sphere but following the ellipsoidal surface. A geodesic distance between two points is the arc length along the great ellipse connecting them, computed using iterative algorithms that solve the direct or inverse geodesic problem. This differs from the Euclidean chord length, which is the straight-line distance through the ellipsoid; for example, on the WGS 84 ellipsoid, the chord length is shorter by a factor approaching cos(Δ/2) for small angular separations Δ, but geodesics are essential for accurate surveying over long distances.^[105]^[108] Local curvature on the ellipsoid is quantified by two principal radii: the meridian radius of curvature M(φ) and the prime vertical radius N(φ). The meridian radius M(φ), which describes curvature in the north-south direction along a meridian, is given by

M(\phi) = \frac{a (1 - e^2)}{(1 - e^2 \sin^2 \phi)^{3/2}},

where a is the semi-major axis and e² = 1 - (b/a)² with b the semi-minor axis. The prime vertical radius N(φ), perpendicular to the meridian in the east-west direction, is

N(\phi) = \frac{a}{\sqrt{1 - e^2 \sin^2 \phi}}.

For the WGS 84 ellipsoid (a = 6,378,137 m, 1/f = 298.257223563 where f is flattening), M(φ) varies from about 6,335,439 m at the equator to 6,399,594 m at the poles, while N(φ) ranges from 6,378,137 m at the equator to larger values toward the poles. These radii enable approximations of small areas as developable surfaces for mapping.^[109]^[107] Area measures on the ellipsoid often employ the authalic sphere, an auxiliary sphere with the same surface area as the ellipsoid, to facilitate equal-area projections. The authalic radius c is c = a √(1 - (1/3)e² - (3/5)e⁴ + ... ) ≈ 6,371,007 m for WGS 84, preserving total surface area through authalic latitude β where sin β = √((1 - e²)/(1 - e² sin² φ)) sin φ. The total surface area of the Earth's reference ellipsoid is approximately 510 million km², closely matching observed values and enabling projections like the Albers equal-area conic.^[110] The fundamental unit of length is the meter, redefined in 1983 by the International Bureau of Weights and Measures as the distance light travels in vacuum in 1/299,792,458 of a second, fixing the speed of light at exactly 299,792,458 m/s for reproducible precision. This definition underpins the WGS 84 ellipsoid parameters, ensuring GPS-derived positions use consistent metric scaling without ambiguity in long-range computations.^[111]^[106] For points above the ellipsoid, such as ground stations, distance calculations incorporate ellipsoidal height h via corrections to geodesic formulas; the effective distance scales by factors like (R + h)/R, where R is the local radius (e.g., N(φ) + h for east-west), reducing measured slopes or baselines by up to 0.1% at h = 100 m to account for the expanded geometry. These adjustments are critical in high-precision GPS applications to align ellipsoidal and topographic distances.^[112]^[113]

Dynamic Aspects

Observational Concepts

Geodetic observations form the foundation of spatial and gravitational measurements used to model Earth's shape, orientation, and gravity field. These observations are categorized into geometric, which determine positions and distances through direct spatial measurements; gravimetric, which quantify gravitational anomalies to infer subsurface mass distributions; and astronomic, which measure deflections of the vertical to link local gravity directions with celestial references.^[114]^[115]^[116] Error propagation in geodetic analysis quantifies how uncertainties in observations affect estimated parameters, typically represented by variance-covariance matrices that capture both variances and covariances among variables. The Gauss-Markov model underpins least-squares adjustments in geodesy, assuming observations as linear combinations of unknown parameters plus random errors with zero mean and known covariance structure, yielding the best linear unbiased estimator.^[117]^[118]^[119] In this framework, the adjusted values \hat{x} are computed via the normal equations solution:

\hat{x} = (A^T P A)^{-1} A^T P l

where A is the design matrix relating observations to parameters, P is the weight matrix inversely proportional to observation variances, and l denotes the vector of observations. This formulation minimizes the weighted sum of squared residuals, ensuring optimal parameter estimates under the model's assumptions.^[120]^[121] Stochastic models describe the random components of geodetic observations, distinguishing between white noise, random walk processes, and systematic effects. For Global Navigation Satellite Systems (GNSS), receiver and satellite clocks are often modeled as random walks to account for their integrated phase errors over time, while tropospheric delays—arising from atmospheric refraction—are treated as systematic biases mitigated through random walk parameterization with process noise tuned to weather model predictions.^[122]^[123]^[124] Quality measures assess the reliability of geodetic adjustments, with confidence ellipsoids visualizing the uncertainty volumes around estimated positions derived from the variance-covariance matrix at a specified probability level, such as 95%. The chi-squared test evaluates residuals by comparing the sum of weighted squared residuals to its expected value, detecting inconsistencies or model misspecifications if the test statistic exceeds critical thresholds.^[125]^[126]^[127] For dynamic surveys involving time-varying phenomena, Kalman filtering provides an advanced recursive framework for real-time processing, sequentially updating state estimates with new observations while propagating uncertainties through prediction and correction steps, thus enabling efficient handling of kinematic trajectories in GNSS-based monitoring.^[128]^[129]

Temporal Changes

Temporal changes in Earth's geodetic parameters arise from dynamic processes that alter the planet's shape, gravity field, and rotation over various timescales, ranging from hours to centuries. These variations, monitored primarily through space-based geodesy such as satellite gravimetry and Global Navigation Satellite Systems (GNSS), reflect interactions between the solid Earth, oceans, atmosphere, and cryosphere. Understanding these changes is essential for refining reference frames and interpreting long-term trends in global positioning.^[130] Solid Earth tides represent periodic deformations of the Earth's crust induced by the gravitational attractions of the Moon and Sun, causing vertical displacements of up to 30–40 cm and horizontal shifts of several centimeters. These tides occur with diurnal and semidiurnal periods, superimposed on longer-term signals, and are modeled using Love numbers—nondimensional parameters that quantify the Earth's elastic response to tidal forcing, with the vertical Love number h_2 typically around 0.6 for the Earth. Observations from superconducting gravimeters and GNSS confirm that these deformations must be corrected in geodetic measurements to achieve millimeter-level accuracy.^[131]^[132] Post-glacial rebound, or glacial isostatic adjustment, involves the ongoing uplift of Earth's crust in regions previously burdened by ice sheets during the last glacial maximum, approximately 20,000 years ago. In Scandinavia, uplift rates range from 3–4 mm/yr along the coastline to nearly 10 mm/yr in the Bothnian Bay, driven by viscoelastic relaxation of the mantle following ice melt. Similarly, around Hudson Bay in North America, rates reach up to 10 mm/yr, reflecting the former extent of the Laurentide Ice Sheet. These motions, detected via GNSS and satellite altimetry, contribute to regional sea level variations and are modeled using viscoelastic Earth structures to isolate contemporary climate signals.^[133]^[134] Tectonic motions encompass gradual plate movements and abrupt coseismic displacements, continuously tracked by networks like the International GNSS Service (IGS). Interseismic strain accumulation occurs at rates of several millimeters to centimeters per year along plate boundaries, while earthquakes produce sudden slips; for instance, the 2011 Tohoku-oki event (Mw 9.0) generated coseismic displacements exceeding 5 m horizontally along Japan's coast due to thrust faulting on the Pacific-North American plate boundary. IGS stations provide dense time series for modeling these deformations, enabling the estimation of fault slip distributions and postseismic relaxation over years.^[135]^[130]^[136] Climate-driven changes in sea level and mass redistribution are quantified through missions like the Gravity Recovery and Climate Experiment (GRACE) and its follow-on, distinguishing between steric effects (thermal expansion of seawater) and barystatic changes (additions or losses of water mass). GRACE data reveal accelerating mass loss from the Greenland Ice Sheet, averaging about 280 Gt/yr from 2006–2018, equivalent to roughly 0.8 mm/yr of global sea level rise, primarily from surface melting and iceberg calving. Recent GRACE-FO data indicate Greenland mass loss averaged ~270 Gt/yr from 2002–2023. Barystatic contributions from land ice dominate recent sea level acceleration, comprising around 1.7 mm/yr since 1993, while steric effects account for 1.2–1.8 mm/yr over the same period; overall, global mean sea level has risen at 3.7 mm/yr (3.2–4.2 mm/yr) from 2006–2018, with the rate accelerating to approximately 4.5 mm/yr as of 2023–2024. These signals, including groundwater depletion and ocean mass gain, are isolated via forward modeling of gravity anomalies.^[137]^[138]^[139]^[140]^[141] Variations in Earth's rotation include polar motion and changes in length of day (LOD), influenced by angular momentum exchanges among the solid Earth, atmosphere, oceans, and core. The Chandler wobble, a free nutation of the rotation axis relative to the crust, has a period of approximately 435 days and amplitude up to 0.2 arcseconds, excited by atmospheric and oceanic loading. Polar motion also features annual and semiannual components of about 100 milliarcseconds, driven by seasonal mass redistributions, alongside secular drift of 3.5 milliarcseconds per year. LOD fluctuations, on intraseasonal to decadal scales, arise from axial angular momentum variations, with seasonal amplitudes of several milliseconds linked to atmospheric winds and El Niño events; decadal changes of up to 7 milliseconds reflect core-mantle interactions. These are monitored via space geodesy and modeled using excitation functions from global fluid models.^[142] Modeling temporal geodetic changes relies on time series analysis of GNSS and gravimetry data to separate trends, periodic signals, and noise. Techniques such as spline interpolation capture non-linear deformations, like postseismic transients or seasonal hydrological loading, while accounting for stochastic processes in velocity estimation. Future projections to 2100, integrated into IPCC assessments, forecast amplified changes under emissions scenarios; for example, global mean sea level medians are 0.44 m (likely range 0.32–0.62 m) under SSP1-2.6 or 0.77 m (0.63–0.98 m) under SSP5-8.5 relative to 1995–2014, with Greenland contributing 0.06–0.13 m and post-glacial rebound modulating regional rates by up to 10 mm/yr. These models incorporate viscoelastic adjustments and climate forcings for scenario-based predictions.^[143]^[144]^[137]^[145]

Applications and Modern Advances

Geophysical and Environmental Applications

Geodesy contributes significantly to geophysics through the analysis of gravity anomalies, which reveal subsurface mass distributions useful for mineral exploration. Isostatic residual gravity anomalies isolate local geological features by removing the effects of crustal thickness variations and long-wavelength trends, enabling the identification of potential ore deposits and structural traps. For instance, in New Mexico, these anomalies highlight basins and uplifts associated with mineral resources, aiding targeted exploration efforts.^[146] Additionally, geodetic gravity measurements support models of crustal thickness via isostatic compensation, where the relationship between topography and gravity anomalies indicates variations in lithospheric structure. Techniques such as the analysis of geoid-to-topography ratios and spectral admittance functions have been applied to estimate crustal thicknesses on Earth and other planets, providing insights into tectonic processes.^[147] In environmental science, geodesy facilitates monitoring of sea level rise using satellite altimetry, which measures ocean surface heights to quantify global changes driven by climate warming. Missions like GRACE and GRACE-FO further enhance this by detecting gravity variations linked to mass redistribution, including contributions from ice melt to sea level. These satellites have tracked groundwater depletion in regions like California's Central Valley, revealing long-term declines in terrestrial water storage due to drought and agricultural overuse, with mass loss rates exceeding 20 gigatons per year in peak drought periods.^[148]^[149] Geodetic techniques are vital for hazard assessment, particularly in volcanic monitoring where interferometric synthetic aperture radar (InSAR) detects surface deformations signaling unrest. InSAR has revolutionized volcano geodesy since the 1990s, capturing millimeter- to centimeter-scale displacements over broad areas, often identifying 1-10 cm of pre-eruptive inflation or deflation that informs eruption forecasting. For tsunamis, high-resolution bathymetric data derived from geodetic surveys integrates with topographic models to simulate wave propagation and inundation, improving hazard mapping and evacuation planning in coastal zones.^[150]^[151] Climate applications of geodesy include assessing ice sheet mass balance through gravity and altimetry observations, which reveal accelerating losses from Greenland and Antarctica. GRACE data indicate Antarctic ice sheet mass loss of approximately 150 gigatons per year over recent decades, contributing substantially to global sea level rise. Land subsidence monitoring via geodetic methods, such as InSAR, also links to the carbon cycle by tracking deformation over storage sites, where CO2 injection can cause measurable ground movement affecting reservoir integrity.^[152] Beyond natural systems, geodesy supports navigation by integrating inertial navigation systems (INS) with global navigation satellite systems (GNSS), enhancing accuracy in dynamic environments like urban or airborne settings through error correction via geodetic reference frames. In urban planning, 3D city models built from geodetic data enable simulation of development scenarios, optimizing land use and infrastructure placement with precise elevation and positional information.^[153]^[154] Recent advances to 2025 incorporate artificial intelligence for anomaly detection in geodetic datasets, improving the identification of subtle gravity or deformation signals in noisy environments. Quantum sensors for field gravimetry offer higher sensitivity for rapid deployment in disaster zones, enabling real-time monitoring of crustal responses during events like earthquakes.^[155]

Notable Geodesists and Milestones

One of the earliest notable figures in geodesy was Eratosthenes (c. 276–194 BC), a Greek mathematician and geographer who accurately estimated the Earth's circumference around 240 BC by measuring the angle of the sun's rays at noon in Alexandria and Syene (modern Aswan), assuming the sun was directly overhead in Syene and using the known distance between the cities. His calculation yielded approximately 252,000 stadia, equivalent to about 39,000–46,000 km depending on the stade length, which was within 2–15% of the modern value. Hipparchus (c. 190–120 BC), another Greek astronomer, advanced geodetic concepts by developing a system of latitude and longitude grids, dividing the Earth into a spherical coordinate framework that enabled more precise mapping and positional astronomy.^[156] In the 18th and 19th centuries, Jean-Baptiste Delambre and Pierre Méchain conducted the French meridian arc measurement from 1792 to 1798, a triangulation survey along the Paris meridian to determine the Earth's shape and size, which confirmed its oblate spheroid form and provided key data for the meter definition. Their work yielded an arc length that helped establish the Earth's flattening ratio as approximately 1/334. Hermann von Helmholtz (1821–1894) contributed to physical geodesy through his work on potential theory in the mid-19th century, applying mathematical principles to describe the gravitational potential of the Earth as a rotating fluid body, influencing later models of geoid determination. In the 20th century, Martin Hotine (1898–1968), a British geodesist, pioneered three-dimensional geodesy during World War II and postwar mapping efforts, developing the Hotine transformation for coordinate conversions between ellipsoids and introducing concepts for mapping the third dimension in large-scale surveys. Weikko A. Heiskanen (1895–1971) and Helmut Moritz co-authored the seminal textbook Physical Geodesy in 1967, which systematized the field by integrating gravitational potential theory, isostasy, and geoid computation methods, becoming a standard reference for generations of geodesists. Richard H. Rapp (b. 1937) advanced satellite geodesy in the late 20th century through his work on gravity field modeling, using data from missions like Seasat and Geosat to develop high-resolution Earth gravity models such as EGM96, which improved global height systems. Among modern geodesists, Michael G. Sideris (b. 1955) has made significant contributions to numerical methods for gravity field determination, developing least-squares collocation techniques and software for geoid modeling from satellite and terrestrial data, as detailed in his influential papers on Stokes's integral approximations. Wolfgang Torge (1936–2020), a prominent German geodesist, led efforts in establishing global gravity networks, including the coordination of the International Gravity Standardization Net 1971 (IGSN71), enhancing the accuracy of absolute gravity measurements worldwide. Notable women in geodesy include Irene K. Fischer (1907–2009), who developed the World Geodetic System 1966 and advanced geoid modeling, and Gladys West (b. 1930), whose mathematical modeling of Earth's shape contributed to the development of the Global Positioning System (GPS).^[157] Key milestones include the founding of the International Association of Geodesy (IAG) in 1862 at the initiative of Johann Jacob Baeyer, which promoted international cooperation in geodetic measurements and remains the oldest scientific association in Earth sciences. The International Federation of Surveyors (FIG) was established in 1878 in Paris to foster global standards in surveying and mapping, evolving into a key body for geomatics professionals. The Global Positioning System (GPS) achieved full operational capability in 1995, revolutionizing geodesy by enabling precise global positioning at centimeter-level accuracy for civilian use. The Gravity Field and Steady-State Ocean Circulation Explorer (GOCE) satellite, launched by the European Space Agency in 2009, provided the first global high-resolution gravity gradient map, advancing physical geodesy with data on Earth's mean sea surface. In contemporary geodesy up to 2025, Richard B. Langley has been a leading figure in Global Navigation Satellite Systems (GNSS), contributing to ambiguity resolution techniques and multi-GNSS integration, as evidenced by his development of the RTKLIB software enhancements and numerous IAG publications. Prestigious awards like the William Bowie Medal, established by the American Geophysical Union in 1939, recognize outstanding contributions to geophysics, including in geodesy.

References

[1]
What is geodesy? - NOAA's National Ocean Service
Jun 16, 2024 · Geodesy is the science of accurately measuring and understanding the Earth's geometric shape, orientation in space, and gravity field.
[2]
What is Geodesy? - SERC (Carleton)
Dec 14, 2015 · Geodesy is the science of accurately measuring the Earth's size, shape, orientation, mass distribution and how these vary with time.
[3]
[PDF] Geodesy for the Layman - National Geodetic Survey
Mar 16, 1984 · INTRODUCTION. CHAPTERS: I History of Geodesy. II Figure of the Earth. III Geodetic Surveying Techniques. IV Geodetic Systems. V Physical Geodesy.
[4]
About | School of Earth Sciences - The Ohio State University
Geodesy is the science of determining the size and shape of the Earth (including its temporal variation) using measurements primarily (today) of distance, time, ...
[5]
IAG Statutes - IAG - Geodesy
(a) Geodesy is the discipline that deals with the measurement and representation (geometry, physics, and temporal variations) of the Earth and other celestial ...
[6]
IAG Commissions - IAG - Geodesy
IAG Commission 1 focuses on defining, realizing, and improving geodetic reference systems to support studies on Earth's rotation, gravity field, sea level ...
[7]
International Association of Geodesy (IAG) - UNOOSA
The classical definition of Geodesy is the science concerned with the shape, size, and the gravity field of the Earth. Geodesy today is much more than that.
[8]
About Geodesy | Measuring Earth and Its Impact on Daily Life
The task of geodesy is to determine the size, shape, gravity field and rotation of the Earth as a function of space and time.
[9]
[PDF] User Guidelines for Single Base Real Time GNSS Positioning
Apr 1, 2014 · The base station should use a ground plane, choke ring, or a current high quality, geodetic, multipath rejecting antenna while the rover ...
[10]
The 200th Anniversary of the Survey of the Coast | National Archives
Feb 10, 2023 · Ferdinand Hassler founded the Survey of the Coast in 1807 and introduced precise geodetic surveying methods never before used in North America.
[11]
[PDF] geodetic surveying 1940 - 1990
The most important development in the history of surveying occurred early in the 1980's when the NAVSTAR. Global Positioning System (GPS) became operational.
[12]
[PDF] Satellite Geodesy 1958-1964 - NASA Technical Reports Server
Much progress has been made in satellite geodesy during the years. 1958 to 1964, although it is likely that sub- stantial improvements will follow the launch of ...Missing: GPS | Show results with:GPS
[13]
[PDF] GPS and Space-Based Geodetic Methods
In the run-up to the space age in 1955, scientists at the Naval Research Laboratory first pro- posed the application of satellite observations to geodesy. By.
[14]
GRACE | Mission
The GRACE twin satellites, launched 17 March 2002, are making detailed measurements of Earth's gravity field changes and revolutionizing investigations.Missing: geodesy | Show results with:geodesy
[15]
Artificial Intelligence for Geodesy | AI4G Focus Area – GGOS
The Focus Area aims to explore the potential of AI and machine learning methods in improving the quality and accuracy of geodetic observations and products.Artificial Intelligence For... · Introduction To Ai4g · Objectives
[16]
Benefits of Future Quantum Accelerometers for Satellite Gravimetry
Sep 1, 2024 · We investigate the benefits of future quantum accelerometers based on cold atom interferometry (CAI) on current and upcoming satellite gravity mission concepts.2.1 Accelerometers · 3 Sst Simulations And... · 3.2 Sst Simulation Result
[17]
What is the geoid? - National Geodetic Survey - NOAA
Jan 30, 2001 · Here is the one currently adopted at NGS: geoid: The equipotential surface of the Earth's gravity field which best fits, in a least squares ...
[18]
Geoids - an overview | ScienceDirect Topics
The geoid is defined as an equipotential surface of the Earth's gravity field that corresponds to the undisturbed global sea level, serving as a height ...
[19]
[PDF] Geoid Surfaces and Theory - National Geodetic Survey
A geoid height is the ellipsoidal height from an ellipsoidal datum to a geoid ... Bruns Formula: Height Relationships: NHh. +. = R. T r. T gg. Q. P. 2.
[20]
The Bruns formula in three dimensions | Journal of Geodesy
Jan 4, 2014 · The Bruns formula is generalized to three dimensions with the derivation of equations expressing the height anomaly vector or the geoid ...
[21]
Precise geoid computation using Stokes-Helmert's scheme and strict ...
A geoid determination using Stokes' formula requires that all topographic masses be removed to a layer on or below the geoid [1]. To satisfy this condition ...
[22]
Precise Geoid Determination in the Eastern Swiss Alps Using ... - NIH
Nov 2, 2024 · Astrogeodetic observations are the most suitable method for geoid determination in mountainous regions such as the Swiss Alps where satellite- ...<|separator|>
[23]
GRACE | Gravity Recovery and Climate Experiment Mission: GFZ
The primary science objective of the GRACE mission was to measure the Earth's gravity field and it's time variability with unprecedented accuracy.
[24]
GOCE (Gravity field and steady-state Ocean Circulation Explorer)
GOCE was a geodynamics and geodetics mission with the objective of determining the stationary gravitational field. The mission ended in October 2013 after ...
[25]
NGA - Office of Geomatics
Online calculator to determine EGM96 and EGM08 geoid height based on latitude and longitude inputs.EGM96 Geoid Calculator · WGS 84 · Precise Imagery · Coordinate Systems
[26]
[PDF] converting gps height into navd88 elevation with the geoid96 geoid ...
Globally, geoid heights can range from about +75 to -100 meters. In the conterminous United States, geoid heights vary from -8 to -53 meters. In addition to ...
[27]
[PDF] Datums, Heights and Geodesy - National Geodetic Survey
General focus here is on the development of geoid height models to relate datums – not the use of these models in determining GPS-derived orthometric heights.
[28]
[PDF] Why Do We Need a Proper Geoid?
May 8, 2009 · For the geoid as a height datum, the second definition is usually chosen. The reason is that it is easier to realise than the first definition ...
[29]
[PDF] Geometric Reference Systems in Geodesy - Knowledge Bank
traditional horizontal geodetic datum is based on eight parameters: three to define its origin, three to define its orientation, and two to define the ellipsoid ...
[30]
Chapter IV GEODETIC SYSTEMS
Selection of the reference ellipsoid provides the radius and flattening factors. The simplest means of obtaining the other factors to establish the geodetic ...
[31]
[PDF] Geodetic Reference System 1980
The reference ellipsoid is defined to enclose the whole mass of the earth, including the atmosphere; as a visualization, one might, for instance, imagine the ...
[32]
Reference Frames in GNSS - Navipedia - GSSC
Feb 23, 2021 · The main parameters of the reference ellipsoid are nearly the same as those defined by ITRS with the semi major axis (a), second degree ...
[33]
[PDF] STATE PLANE COORDINATES Washington State - WA DNR
Mar 13, 2007 · CLARKE 1866 a = 6,378,206.4 m 1/f = 294.97869821. GEODETIC REFERENCE ... BESSEL 1841 a = 6,377,397.155 m 1/f = 299.1528128. Page 9. THE ...
[34]
[PDF] NOAA Technical Memorandum NOS NGS 82
Feb 26, 2020 · For reference, the Clarke 1866 ellipsoid has a = 6,378,206.4 m (exact) and e2 = 0.00676865799729... (not exact). Page 9. 8 which might come from ...
[35]
[PDF] Geomatics Guidance Note number 7, part 2 – September 2019 - IOGP
Sep 2, 2019 · Ellipsoid: Clarke 1866 a = 6378206.400 metres = 20925832.16 US survey feet ... Ellipsoid: Clarke 1880 (IGN) a = 6378249.2 metres. 1/f ...
[36]
World Geodetic System 1984 (WGS 84) - NGA - Office of Geomatics
Defining Parameters: WGS 84 identifies four defining parameters. These are the semi-major axis of the WGS 84 ellipsoid, the flattening factor of the Earth, the ...
[37]
Coordinate systems, map projections, and transformations—ArcGIS ...
Ellipsoidal coordinate systems reference a mathematically derived spheroidal or ellipsoidal surface. Since they are calculated on a mathematical model ...
[38]
14. Ellipsoids | The Nature of Geographic Information - Dutton Institute
Ellipsoids are commonly used as surrogates for geoids so as to simplify the mathematics involved in relating a coordinate system grid with a model of the ...
[39]
[PDF] NOAA Technical Report NOS NGS 62
Apr 20, 2021 · According to the IAG Global Geodetic Observing System (GGOS) requirements (Plag and Pearlman, 2009), the ITRF is intended to be accurate to 1.0 ...
[40]
Ellipsoidal and Cartesian Coordinates Conversion - Navipedia - GSSC
The ( x , y , z ) ECEF cartesian coordinates can be expressed in the ellipsoidal coordinates ( φ , λ , h ) , where φ and λ are, respectively, the latitude ...<|separator|>
[41]
[PDF] Coordinate Systems - UNOOSA
ECEF to Geodetic and vice versa. X = N + ℎ cos𝜑 cos𝜆. Y = N + ℎ cos𝜑 sin𝜆. Z = N 1 − e. 2. + ℎ sin𝜑. 𝜑 = Latitude. 𝜆 = Longitude. H = Height above Ellipsoid.
[42]
[PDF] A MANUAL FOR GEODETIC COORDINATE TRANSFORMATIONS ...
Jul 3, 1977 · The equations requiredto solve certain coordinate and associated error transformation problems are stated, the notation used is explained, and a ...
[43]
[PDF] Reference Frames Coordinate Systems
– The center of any inertial frame is ALWAYS the solar system barycenter. ... • Except for attitude derived from the 1976 IAU Earth precession model and.
[44]
Universal Transverse Mercator Coordinates - Tools
The Universal Transverse Mercator Coordinate (UTM) system provides coordinates on a world wide flat grid for easy computation.
[45]
[PDF] Coordinates - Naval Postgraduate School
For the ellipsoidal model of the earth, the ellipsoid can be used as the reference surface. This is done and the resulting height is called ellipsoidal height ...<|separator|>
[46]
A tutorial on datums - VDatum - NOAA
Unlike a horizontal datum, a three dimensional datum provides the foundation for accurate determination of ellipsoidal heights.<|control11|><|separator|>
[47]
[PDF] Coordinate systems used in geodesy - National Geodetic Survey
Origin: Is at the center of the reference ellipsoid used for defining the datum in question. 2. w-axis: Coincides with the semiminor axis b of the reference ...
[48]
The Problem with NAD27 | Pima County, AZ
The reference ellipsoid used for NAD 27 was the Clarke ellipsoid of 1866, that of NAD 83 the the Geodetic Reference System of 1980 (GRS 80 ) ellipsoid. The two ...
[49]
[PDF] 4 Survey Datums - Caltrans
The mathematical reference surface used to represent NAD 83 is an ellipsoid called the Geodetic Reference System of 1980 (GRS 80). GRS 80 is a world-wide model.<|separator|>
[50]
[PDF] User Guidelines for Single Base Real Time GNSS Positioning
Height Modernization practices (see http://www.ngs. noaa.gov/heightmod/ ) can produce 2 cm local orthometric height accuracy from static GPS procedures. It is ...Missing: refinements | Show results with:refinements
[51]
[PDF] NOAA Technical Report NOS NGS 60
“The network accuracy of a control point is a value that represents the uncertainty in the coordinates of the control point with respect to the geodetic datum.
[52]
Geodetic datum - Wikipedia
This phenomenon is called datum shift or, more generally, datum transformation, as it may involve rotation and scaling, in addition to displacement. Because ...
[53]
NADCON - North American Datum Conversion Utility - Tools
NADCON was developed in order to facilitate conversion between the datums,. The grids used by the program are based on more than 150,000 horizontal control ...Missing: plate tectonics
[54]
[PDF] NOAA Technical Memorandum NOS NGS-50
A complete horizontal datum consists of: (1) all of the parameters necessary to define a particular coordinate system, and (2) a set of control points whose ...
[55]
ITRF2014: A new release of the International Terrestrial Reference ...
Jul 18, 2016 · As regard to the rotation parameters, we chose the minimum constraint approach, i.e., imposing no net periodic rotation conditions on a set of ...Missing: biennial | Show results with:biennial
[56]
Itrf2020 - IGN
The orientation of the ITRF2020 seasonal signals, expressed in both the CM and CF frames, is defined in such a way that there is no net seasonal rotation of ...Missing: biennial | Show results with:biennial
[57]
[PDF] Comparison and Assessment of Three ITRS Realizations
Jun 12, 2021 · Abstract: A terrestrial reference frame (TRF) is derived based on historical geodetic data and is normally updated every 5–6 years.Missing: biennial | Show results with:biennial
[58]
Height Systems | Orthometric, Ellipsoidal & Normal Heights
### Summary of Height Systems in Geodesy
[59]
Deflection of the Vertical Survey - National Geodetic Survey - NOAA
Local deflection angles can approach many 10s of arc-seconds, and can be determined with an accuracy of 0.1 arc-seconds.
[60]
[PDF] DISPLAYING GEOMETRIC ALTITUDE RELATIVE TO MEAN SEA ...
Jul 9, 2009 · Likewise, we refer to “orthometric height” (height of a point relative to the geoid) here as. “geometric altitude relative to MSL.” 4 ...
[61]
An Improved Approach to the Generation of Hybrid Geoids
May 16, 2017 · The accuracy of the satellite models has been reported to be better than 1 cm to degree and order 70 (Tapley et al., 2005). Based on those ...
[62]
The Vertical Datum: Global Positioning Tutorial
Aug 12, 2024 · The vertical datum is a collection of points with known heights above or below mean sea level, used to measure land sinking.
[63]
Datums and Reference Frames - National Geodetic Survey - NOAA
Dec 7, 2023 · A geodetic datum or reference frame is an abstract coordinate system with a reference surface (such as sea level) that serves to provide known ...Missing: realization | Show results with:realization
[64]
[DM-05-051] Vertical (Geopotential) Datums
Vertical datum – Any level surface (as, for example, mean se level) taken as a surface of reference from which to reckon elevations. WGS84 – World Geodetic ...
[65]
[PDF] A global vertical datum defined by the conventional ... - DiVA portal
Abstract. The geoid, according to the classical Gauss–Listing definition, is, among infinite equipotential surfaces of the Earth's gravity.
[66]
[PDF] EVRF2007 as Realization of the European Vertical Reference System
Since 1994 the IAG Sub-commission for Europe. (EUREF) has enhanced the Unified European Leveling. Network (UELN) and defined a European Vertical. Reference ...
[67]
[PDF] European Vertical Reference System (EVRS) 2007 - EUREF
Normal heights are equivalent to geopotential numbers. c). The EVRS is a zero tidal system1, conforming to the. IAG Resolutions No 16 adopted in Hamburg in 1983.
[68]
[PDF] NGS Onward to 2022: Replacing NAD83 & NAVD88
North American Vertical Datum 1988 (NAVD88) Shortcomings. Approximate Geoid ... Extent of NAPGD2022 Gravimetric Geoid Model (GEOID2022). American Samoa.
[69]
Australian Height Datum | Geoscience Australia
Dec 14, 2017 · The datum surface is that which passes through mean sea level at the 30 tide gauges and through points at zero AHD height vertically below the ...
[70]
[PDF] Implementing a Vertical Datum Referenced to a Gravimetric Geoid ...
May 10, 2012 · In order to address the deficiencies of the NAVD 88, NGS has been engaged in a National Height. Modernization Program, enabling better access to ...
[71]
Sea Level Rise - Florida Climate Center
Nearly half of this rise in sea levels has occurred since 1993, and the rate of global mean level rise since 1993 has been approximately 3.4 mm per year. It ...
[72]
New Datums - National Geodetic Survey - NOAA
New Datums: Replacing NAVD 88 and NAD 83 · Easy access to horizontal and vertical reference systems · Improved horizontal accuracy · Improved vertical accuracy ...FAQs · Naming Convention · Get Prepared · Track Our ProgressMissing: realization | Show results with:realization
[73]
[PDF] Introduction to GPS and other Global Navigation Satellite Systems
Jun 8, 2017 · BeiDou will provide three carrier frequencies foreseen to be interoperable with other systems ... conjunction with GPS, GLONASS or Galileo, will ...
[74]
Galileo General Introduction - Navipedia - GSSC
Jul 10, 2025 · Galileo is designed to be fully interoperable with GPS, GLONASS, and BeiDou, providing enhanced reliability and precision.
[75]
Chapter 5: Land Surveying and GPS - Penn State
The Global Positioning System (GPS) employs trilateration to calculate the coordinates of positions at or near the Earth's surface.
[76]
[PDF] Chapter 9 Stand-Alone Positioning
Code and Carrier-Phase Measurement Model Each GNSS signal has two components that are tracked by the receiver: the code and the carrier. The code is modulated ...
[77]
The Pseudorange Equation | GEOG 862 - Dutton Institute
On the left-hand side, the pseudorange measurement equals rho, the true range, plus all of these other factors.
[78]
PPP Fundamentals - Navipedia - GSSC
Dec 18, 2021 · In order to achieve the expected centimeter-level accuracy, PPP relies on very accurate error models that are not generally required when ...Missing: sources | Show results with:sources
[79]
Lesson 8: Real-Time Global Positioning System Surveying
Real-time positioning is built on the idea that, with the important exceptions of multipath and receiver noise, these GPS/GNSS error sources are correlated.
[80]
Real-Time PPP-RTK Performance Analysis Using Ionospheric ... - NIH
In this study, we present and analyze the real-time PPP-RTK user performance using ionospheric corrections from multi-scale regional networks.Missing: sources | Show results with:sources
[81]
VLBI | Earth - NASA
VLBI is unique in its ability to define an inertial reference frame and to measure the Earth's orientation in this frame.
[82]
IVS Home Page
The International VLBI Service for Geodesy and Astrometry (IVS) provides a service to support geodetic, geophysical and astrometric research and operational ...
[83]
ILRS Home Page
Scientific products derived using SLR and LLR data include precise geocentric positions and motions of ground stations, satellite orbits, components of Earth's ...List of Stations · About ILRS · Meetings · Satellite MissionsMissing: coordinates | Show results with:coordinates<|separator|>
[84]
Network | Site Information - ILRS
May 22, 2025 · SLR Station Coordinates (SLRF2020); SLR System Eccentricities (Up, North, East and X, Y, Z), maintained by the ILRS CB and by the CDDIS; SLR ...
[85]
[PDF] Realization of Terrestrial Reference Frame at the observation level ...
Combines techniques (GPS + SLR + VLBI) at the observation level. • Dispenses with traditional ground ties for SLR/GPS. • Connection between SLR and GPS is ...Missing: applications | Show results with:applications
[86]
Earth Sciences - Navipedia - GSSC
Sep 7, 2018 · The global network of IGS stations can provide continuous information about motions of the continents or smaller tectonic plates. In addition at ...
[87]
GPS And Tectonics | How GPS Works - GPS Spotlight - UNAVCO.org
GPS aids us in better understanding our planet by allowing us to measure how the surface of the Earth moves. Much of this motion causes earthquakes, builds ...Missing: applications | Show results with:applications
[88]
[PDF] 2024 IGS Technical Report - Index of /pub
We have made progress in the inclusion of more and more systems into our operational product lines. With the regular reprocessing efforts the. IGS strongly ...
[89]
RTK from space: a novel and low-cost GNSS augmentation ...
Mar 19, 2025 · In this paper, we proposed a novel and low-cost Global Navigation Satellite System (GNSS) augmentation technique named “Real Time Kinematic (RTK) from space”
[90]
5.5 Land Surveying and Conventional Techniques for Measuring ...
We will now introduce two techniques that surveyors use to create and extend control networks (triangulation and trilateration) and two other techniques used ...Missing: terrestrial | Show results with:terrestrial
[91]
[PDF] CHAPTER 2 POSITIONING - International Hydrographic Organization
The first two methods accomplish the basic control networks, primary or of inferior order, via triangulation, trilateration and traverse operations. Afterwards, ...
[92]
5: Land Surveying and GPS - Geosciences LibreTexts
Aug 11, 2020 · Trilateration is an alternative to triangulation that relies upon distance measurements only. Electronic distance measurement technologies make ...
[93]
[PDF] Chapter 1 Surveying - USDA
Geodetic surveying measures all elevations from a level surface, and because all instruments use a straight reference plane, this requires all instrument ...
[94]
Introduction to EDM, Theodolite and Total station
May 24, 2021 · There are three methods for angular measurements: Normal method; Repetition method; Direction method ( Reiteration method). Components of Total ...
[95]
[PDF] Modern Ground Networks for Space Geodesy Applications - cddis
Aug 13, 2010 · Provides unique information on the deformation due to natural hazards (volcanoes, landslides, earthquakes, etc.) • At right is an InSAR map of ...
[96]
[PDF] NOAA Technical Report NOS NGS 69 - National Geodetic Survey
Here, we present a foundational introduction to the various types of geophysical phenomena that are causing both size and shape change to the geoid, geodetic ...
[97]
[PDF] Refraction Tests by the - National Geodetic Survey
... least-squares adjustments to determine elevations for each survey. (See appendixes. A and B.) ELEVATION DIFFERENCES. CONTAINING. REFRACTION ERRORS. Each set of ...
[98]
Use of Total Least Squares Adjustment in Geodetic Applications
Mar 16, 2024 · The method takes into account the errors on the side of the unknowns or in the matrix of adjustment model A , also with the possibility of ...
[99]
(PDF) Detection of Crustal Deformation Using a Gravimetric and ...
Aug 6, 2025 · The InSAR displacement maps reveal subsidence and uplift rates up to 8.8 mm/year, and 7.7 mm/year, respectively. This observed deformation in ...
[100]
Recent Advances in Modelling Geodetic Time Series and ... - MDPI
Recent advances include modeling geodetic time series using instruments like GNSS, GRACE, InSAR, and TLS to monitor geodynamics, climate change, and natural ...
[101]
[PDF] Direct and inverse solutions of geodesics - National Geodetic Survey
This paper presents compact, iterative direct and inverse solutions for geodesics on an ellipsoid, using nested equations for elliptic terms and only three ...
[102]
Intersections on the sphere and ellipsoid
The problems of intersection and resection in two dimensions are well known in geodetic surveying, and their solutions can be found in most textbooks on ...
[103]
[PDF] Map projections--a working manual - USGS Publications Warehouse
Nov 12, 1987 · This manual describes map projections, useful for those interested in their philosophy, history, or mathematics, and includes non-mathematical ...
[104]
[PDF] Algorithms for geodesics - arXiv
The goal of this paper is to adapt the geodesic methods of Helmert (1880) and his predecessors to modern comput- ers. The current work goes beyond Vincenty in ...
[105]
NGA Geomatics - WGS 84
### Summary of WGS 84 Use of Meter and Ellipsoid Parameters Relevant to Units
[106]
[PDF] GEODETIC POSITION COMPUTATIONS
Two classical geometric geodetic computation problems are treated; they are called the direct and inverse geodetic problems. There are various approaches ...
[107]
Geodetic Uses of Artificial Satellites - Astrophysics Data System
... chord distance), and not the distance along the terrestrial ellipsoid (length of the geodesic). Next let us consider conditions imposed on the directions. A ...
[108]
[PDF] Radius of the Earth - Radii Used in Geodesy
The radius of curvature for the latitude is found by extending the line perpendicular to the ellipsoid (the geodetic vertical) down until it hits the polar axis ...
[109]
[PDF] General Theory of Equivalent Projections - National Geodetic Survey
Jan 1, 2006 · In order to take account of the spheroid, we decided to project the authalic sphere as a further study in the theory of map projection. Let ...
[110]
Meter | NIST - National Institute of Standards and Technology
May 9, 2019 · This definition also locked the speed of light at 299,792,458 meters per second in a vacuum. Length was now no longer an independent standard ...Missing: WGS84 | Show results with:WGS84
[111]
Elevation Factor | GEOG 862 - Dutton Institute
It is possible to refine the calculation of the elevation factor by using an average of the actual radial distances from the center of the ellipsoid to the end ...
[112]
Height Above Ellipsoid (HAE) Calculations - Blue Marble Geographics
Inverse distance calculation using the geodesic method for points that contain ellipsoidal height values are calculated on an ellipsoid that has been expanded ...
[113]
[PDF] Determination of an Orthometric Height Profile in the Okavango ...
Four types of observations can be used for the determination of the geoid: - gravimetric,. - astronomical,. - satellite or. - geodetic observations. In general ...
[114]
[PDF] CHAPTER 2 POSITIONING
These local Geoids are gravimetric and independent from Geoid undulation values obtained from combining GPS and geometric levelling observations; they are ...
[115]
[PDF] Determination of Astro-Geodetic Vertical Deflections using Digital ...
May 10, 2012 · In the scope of this study we discuss the astro-geodetic observations performed on a test station using a CCD camera, Schmidt-Cassegrain type ...<|separator|>
[116]
[PDF] Notes on Adjustment Computations - School of Earth Sciences
Dec 26, 2018 · 3.1 The Least-Squares Solution Within the Gauss-Markov Model . ... 3.2 Estimated Variance Component Within the Gauss-Markov Model .
[117]
[PDF] Statistical uncertainty and error propagation
Mar 27, 2014 · This lecture is about some of the basics of uncertainty: the use of variances and covari- ances to express statistical uncertainty (“error ...
[118]
[PDF] NOAA Technical Memorandum NOS NGS 74 On Least Squares ...
May 1, 2018 · In any case, we do not recommend the least-squares solution within the. Gauss-Markov Model with Stochastic Constraints described in Section 2.
[119]
[PDF] Adjustment Computations - School of Earth Sciences
Feb 4, 2021 · Equations. 4.1 Model Definition. In the least-squares adjustment within the model of condition equations, the un- known parameters ξ are not ...
[120]
Least Squares | myGeodesy
Lecture notes on the topic of least squares with plane surveying applications (222 pages). These notes contain 10 chapters plus an appendix on matrix algebra.
[121]
Optimum stochastic modeling for GNSS tropospheric delay ...
Dec 23, 2016 · We propose to take benefit from numerical weather prediction models to define optimum random walk process noise.
[122]
An Improved Stochastic Model for the Geodetic GNSS Receivers ...
Feb 18, 2024 · The receiver tracking error stochastic (RTES) model can improve GNSS precise point positioning (PPP) performance under ionospheric ...
[123]
(PDF) Optimum stochastic modeling for GNSS tropospheric delay ...
Dec 23, 2016 · In GNSS data processing, the station height, receiver clock and tropospheric delay (ZTD) are highly correlated to each other.
[124]
[PDF] Adjustment Computations - School of Earth Sciences
Jan 19, 2021 · function of the random errors was required for the derivation of their least-squares prediction. Page 126. 116. CHAPTER 5. THE GAUSS-MARKOV ...
[125]
A new relationship between the quality criteria for geodetic networks
Aug 4, 2018 · The goal of this paper is to present a new relationship between the quality criteria for geodetic networks.Missing: collinearity | Show results with:collinearity
[126]
[PDF] 5 Classifications of Accuracy and Standards - Caltrans
If a network has a standard error of unit weight close to 1.0, no high observation residuals, and still does not pass the Chi-square test, the network should be.
[127]
Real-time deformation monitoring with GPS and Kalman Filter
Jun 24, 2014 · The main purpose of this research is to develop a real-time GPS monitoring system with the aid of a Kalman Filter for use in an active tectonic region.
[128]
Evaluation of linear Kalman filter processing geodetic kinematic ...
The purpose of this paper is to assess the efficiency of the linear Kalman filter as a method for the estimation of kinematic process observed with ...
[129]
Supporting the ITRF - International GNSS Service (IGS)
The IGS publicly available geodetic products, together with the reference frame stations allow further densifying and providing a universal access to the ITRF.
[130]
Tidal Love and Shida numbers estimated by geodetic VLBI - PMC
Tidal forces induced by the gravitational attraction of the Moon and Sun cause a periodic deformation of the Earth, called Earth tides. The total displacement ...
[131]
A rotational ellipsoid model for solid Earth tide with high precision
Nov 18, 2024 · Solid Earth tide represents the response of solid Earth to the lunar (solar) gravitational force ... Love numbers for a rotating spheroidal Earth ...
[132]
Glacial isostatic adjustment and the radial viscosity profile from ...
Nov 8, 2002 · We observe an elliptical region of positive uplift velocities, from 3 to 4 mm yr−1 around the Scandinavian coastline and increasing to almost 10 ...
[133]
[PDF] Observation of glacial isostatic adjustment in ''stable'' North America ...
Jan 26, 2007 · Vertical velocities show present-day uplift. ( 10 mm/yr) near Hudson Bay, the site of thickest ice at the last glacial maximum. The uplift ...
[134]
Slip-History database :: 2011 Tohoku-oki Earthquake
The coast of Japan sprang east, into the ocean, in some places by as much as 5 m (16 feet), the motions are indicated by the blue arrows. Ocean bottom GPS ...
[135]
Tectonic motion monitoring of Antarctica using GPS measurements
Aug 6, 2025 · In this study, data obtained from IGS stations between 2004 and 2012 ... motion including those estimated from space geodetic measurements.
[136]
[PDF] Chapter 9: Ocean, Cryosphere and Sea Level Change
Chapter 9 covers Ocean, Cryosphere and Sea Level Change, including oceans, sea ice, and sea level change. It includes sections on ocean surface, heat, salinity ...
[137]
Greenland | Vital Signs - NASA Sea Level Change Portal
Since 2002, Greenland has lost ice, causing sea levels to rise by about 2.5 inches (63mm). The loss is measured in gigatonnes of ice mass.Missing: 400 barystatic
[138]
Barystatic sea level change observed by satellite gravimetry - PNAS
Jun 30, 2025 · We derive three decades of barystatic sea level (BSL) change estimates since 1993, providing a long record of satellite gravimetry-based BSL estimates.
[139]
[PDF] 3.09 Earth Rotation Variations – Long Period
Atmospheric angular momentum fluctuations, length-of-day changes and polar motion. ... Earth's polar motion by atmo- spheric angular momentum variations ...
[140]
Detection and interpretation of the time-varying seasonal signals in ...
The time-varying periodic variations in Global Navigation Satellite System (GNSS) stations affect the reliable time series analysis and appropriate ...
[141]
Velocity estimation of GNSS coordinate time-series accounting for ...
SUMMARY. We revisit the challenge of estimating velocities from global navigation satellite system (GNSS) coordinate time-series by incorporating stochasti.2 Materials And Methods · 2.1 Trajectory Model And... · 3 Test On Synthetic...
[142]
[PDF] ISOSTATIC RESIDUAL GRAVITY ANOMALIES OF NEW MEXICO
The Bouguer gravity anomaly map indicates a strong regional gradient due to changing crustal thickness, which obscures the anomalies of geologic interest.Missing: isostasy | Show results with:isostasy
[143]
[PDF] 10.05 Gravity and Topography of the Terrestrial Planets
These methods include crustal thickness modeling, the analysis of geoid/ topography ratios (GTRs), modeling of the spectral admittance and correlation functions ...
[144]
GRACE (Gravity Recovery And Climate Experiment) - eoPortal
GRACE is an international cooperative US-German dual-minisatellite SST (Satellite-to-Satellite Tracking) geodetic mission with the overall objective to obtain ...<|separator|>
[145]
GRACE and GRACE-FO track California's land water changes
Aug 7, 2022 · GRACE and GRACE-FO map changes in Earth's gravity from water movement, showing California's long-term water decline due to drought and ...Missing: altimetry | Show results with:altimetry
[146]
Volcano geodesy using InSAR in 2020: The past and next decades
Feb 12, 2022 · The study of volcano deformation has grown significantly through they year 2020 since the development of interferometric synthetic aperture ...Missing: tsunami bathymetry
[147]
Integrating Topographic and Bathymetric Data for High-Resolution ...
Jan 19, 2023 · This article focuses on the creation of three-dimensional digital elevation models (DEMs) of Earth's surface to support tsunami modeling, inundation mapping, ...Missing: geodesy volcano deformation InSAR<|separator|>
[148]
Ice Sheet Mass Changes over Antarctica Based on GRACE Data
Oct 11, 2024 · This study primarily examines changes in the Antarctic ice sheet's mass balance as independently derived from GRACE and GRACE Follow-On ...
[149]
Street Smart: 3D City Mapping and Modeling for Positioning with ...
Jul 1, 2015 · A particle-filter-based positioning method using a 3D map to rectify the errors created by multipath and non-line-of-sight signals on the positioning result.Missing: planning | Show results with:planning<|separator|>
[150]
3D Viewpoint Management and Navigation in Urban Planning - MDPI
3D geovisualization is essential in urban planning as it assists the analysis of geospatial data and decision making in the design and development of land ...Missing: geodesy | Show results with:geodesy
[151]
The Potential of Hydrogeodesy to Address Water‐Related and ...
Nov 2, 2024 · Satellite-based geodetic observations can help track freshwater availability by measuring the temporal variation of geometry and gravity over ...