Geopotential
Geopotential is the gravitational potential energy per unit mass due to Earth's gravity at a given point, equivalent to the work required to lift a unit mass from sea level to that height against the varying gravitational field, and it has units of square meters per second squared (m²/s²).[1] In atmospheric science, it is mathematically defined as \Phi(z) = \int_0^z g(z') \, dz', where g(z') is the acceleration due to gravity at height z', providing a measure that accounts for gravity's decrease with altitude rather than assuming constant gravity.[2] A related concept, geopotential height (Z), expresses this potential in terms of an equivalent height under standard gravity, calculated as Z = \Phi / g_0, where g_0 = 9.80665 m/s² is the global average gravity at mean sea level; this yields units of geopotential meters (gpm), which approximate geometric height but correct for gravitational variations.[1][2] Geopotential height is crucial in meteorology for analyzing pressure surfaces, such as the 500 hPa level, where it helps map atmospheric circulation patterns and forecast weather systems by providing a consistent vertical coordinate that aligns with the hydrostatic equation.[3] In geodesy, geopotential defines equipotential surfaces where the potential is constant, with the geoid representing the specific surface coinciding with mean sea level at rest, excluding tides and currents; this surface, quantified by a reference value like W_0 = 62,636,853.4 m²/s² globally, serves as the zero-height reference for accurate height measurements and vertical datums.[4] Applications include coastal water modeling, flood risk assessment, and integrating geometric heights with gravity data to support navigation, surveying, and global height systems like the North American-Pacific Geopotential Datum of 2022 (NAPGD2022).[4]Fundamental Concepts
Definition
The geopotential at a point in the Earth's gravity field is defined as the work done per unit mass to move a test mass from a reference equipotential surface, typically the geoid, to that point along a plumb line, equivalent to the line integral of the gravity vector over the path.[4] This represents the potential energy associated with the position in the field, distinguishing it from kinetic or other forms of energy.[5] Mathematically, the geopotential \Phi is given by \Phi = -\int_{\text{geoid}}^{\text{point}} \mathbf{g} \cdot d\mathbf{s}, where \mathbf{g} is the local gravity acceleration vector and d\mathbf{s} is the infinitesimal displacement along the integration path.[6] Since the gravity field is conservative (irrotational and source-free outside the Earth), the value of \Phi is independent of the specific path taken between the reference surface and the point, ensuring that all points on an equipotential surface share the same geopotential value.[5] In contrast to geometric height, which simply measures the radial or plumb-line distance from a reference ellipsoid without regard to gravity variations, geopotential accounts for the non-uniform strength and direction of gravity across the Earth's surface, leading to equipotential surfaces that more accurately reflect physical leveling and water flow behavior.[4] The concept of geopotential emerged in the 19th century within physical geodesy, as scientists sought to model the Earth's irregular figure using gravitational potential theory to explain phenomena like the geoid and isostatic equilibrium.[5] Geophysicists such as George Biddell Airy contributed foundational work by applying potential calculations to assess the gravitational effects of mountain masses on latitude observations, laying groundwork for understanding the Earth's oblate shape and density variations.[7] The total geopotential incorporates a rotational component from Earth's spin, which perturbs the purely gravitational field but is addressed separately in detailed analyses.[4]Units and Properties
The geopotential, denoted as Φ, is a scalar quantity representing the potential energy per unit mass in the Earth's gravity field, with standard SI units of square meters per second squared (m²/s²), equivalent to joules per kilogram (J/kg). This unit arises because the geopotential is defined as the work done against gravity to move a unit mass from a reference surface to a given point, integrating the gravitational acceleration along the path. In practice, it is often expressed in geopotential meters (gpm), a derived unit where the numerical value corresponds to the geopotential divided by the standard gravity constant g₀ = 9.80665 m/s², making 1 gpm equivalent to approximately 9.80665 m²/s² for conversion purposes in atmospheric and geodetic applications.[8][9][10] As a scalar field, the geopotential possesses key physical properties that underscore its role in geophysical modeling: it is conservative, meaning the work to move a unit mass between two points is path-independent, and the associated gravity vector field g = -∇Φ is irrotational (∇ × g = 0) and divergenceless (∇ · g = 0) in regions outside mass concentrations, such as in the atmosphere or vacuum. These characteristics stem from the geopotential satisfying Laplace's equation (∇²Φ = 0) in source-free regions, ensuring harmonic behavior and enabling unique solutions for boundary value problems in gravimetry. Additionally, the geopotential forms closed equipotential surfaces where Φ is constant, with the geoid defined as the specific equipotential surface (often referenced to zero for relative heights) that best approximates global mean sea level in a least-squares sense.[8][6][4] For practical computations, the geopotential relates to geometric height h above the reference surface through the approximation Φ ≈ g₀ h, where g₀ = 9.80665 m/s² serves as the conventional mean gravity value, allowing conversion between potential differences and vertical distances with high accuracy for small heights (e.g., in the troposphere). This relation holds because the geopotential height Z = Φ / g₀ provides a dynamically equivalent measure to geometric height, differing by less than 0.3% at 10 km altitude due to gravity variations. Typical absolute values of the geopotential on the Earth's surface, referenced to the geoid equipotential W₀, are around 62.6 million m²/s², with the conventional global value W₀ = 62,636,853.4 m²/s² adopted by the International Association of Geodesy (IAG); relative to the geoid (set to 0), values increase with height but remain on the order of tens of thousands of m²/s² in the lower atmosphere, while the full-scale absolute magnitude establishes the baseline for polar regions where gravitational contributions dominate.[10][4][4]Physical Components
Gravitational Potential
The gravitational potential V arises from the Newtonian gravitational attraction due to Earth's mass distribution and represents the primary component of the geopotential. For a point mass M at a distance r from the observation point, the potential is given byV(r) = -\frac{GM}{r},
where G is the gravitational constant. This formulation assumes a spherically symmetric mass, but Earth's oblate shape and heterogeneous internal structure require an extension to account for deviations from perfect symmetry.[11] To describe the potential outside Earth, the point-mass expression is generalized using a spherical harmonic expansion, which captures the effects of mass irregularities. The full exterior gravitational potential at a point with geocentric radius r, latitude \phi, and longitude \lambda is expressed as
V(r, \phi, \lambda) = -\frac{GM}{r} \sum_{n=0}^{\infty} \sum_{m=0}^{n} \left( \frac{r_e}{r} \right)^n \bar{P}_{n m}(\sin \phi) \left( \bar{C}_{n m} \cos m\lambda + \bar{S}_{n m} \sin m\lambda \right),
where r_e is the reference equatorial radius, \bar{P}_{n m} are the fully normalized associated Legendre functions, and \bar{C}_{n m}, \bar{S}_{n m} are the fully normalized Stokes coefficients describing the mass distribution (with the n=0, m=0 term yielding the monopole -\frac{GM}{r}). Higher-degree terms (n ≥ 2) reflect the oblateness and lateral variations, with the zonal harmonic \bar{C}_{2 0} \approx -4.84 \times 10^{-4} (corresponding to the unnormalized J_2 \approx 1.083 \times 10^{-3}) dominating the equatorial bulge effect. This expansion converges for r > r_e and is derived from Poisson's integral over Earth's mass density.[11][12] Earth's internal mass distribution—comprising the dense iron-nickel core (about 32% of total mass), the silicate mantle (about 67%), and the thin crust (about 1%)—generates these potential anomalies. The core and mantle contribute primarily to low-degree (long-wavelength) terms due to large-scale density contrasts, such as core-mantle boundary undulations, while crustal thickness variations and density heterogeneities produce high-degree (short-wavelength) anomalies observable in regional gravity data. Global gravity models like EGM2008 integrate satellite altimetry, gravimetry, and tracking data to estimate Stokes coefficients up to degree and order 2159, with additional terms to degree 2190 for finer resolution; modern successors, such as the planned EGM2020, aim to extend this to similar or higher degrees using updated datasets from missions like GRACE-FO. As of 2025, although EGM2020 is pending release, models like GOCO06s and monthly GRACE-FO solutions (e.g., AIUB RL02) provide updated Stokes coefficients using post-2018 data.[12][13][14][15] Outside the Earth's mass, the gravitational potential satisfies Laplace's equation
\nabla^2 V = 0,
ensuring harmonic behavior in the exterior domain, with boundary conditions imposed by the surface gravity and mass continuity at r = r_e. This property allows the spherical harmonic series to uniquely represent the field from surface measurements, facilitating model inversion for internal structure.[11] The total geopotential includes this gravitational component plus a centrifugal term due to rotation.[12]
Centrifugal Potential
The centrifugal potential originates from the fictitious centrifugal force experienced in the Earth's rotating reference frame, acting outward perpendicular to the axis of rotation. This potential is derived from the work done by the centrifugal force and is expressed as \phi_c = -\frac{1}{2} \omega^2 \rho^2, where \omega denotes Earth's angular velocity, valued at $7.292115 \times 10^{-5} rad/s, and \rho represents the perpendicular distance from the rotation axis. In terms of spherical coordinates, \rho = r \sin \theta, with r as the geocentric radius and \theta as the colatitude (angle from the north pole). The explicit form thus becomes \phi_c = -\frac{1}{2} \omega^2 (r \sin \theta)^2. This formulation highlights its dependence on latitude: the potential reaches maximum magnitude at the equator (\theta = 90^\circ, \sin \theta = 1), where \rho is largest, and vanishes at the poles (\theta = 0^\circ or $180^\circ). Its overall magnitude is minor relative to the gravitational potential, constituting approximately 0.2% at Earth's surface near the equator.[11][16] Integrated into the total geopotential as an additive term to the gravitational potential, the centrifugal potential distorts the otherwise spherical equipotentials into oblate spheroids, aligning with Earth's observed equatorial bulge and polar flattening. This effect arises because the outward centrifugal contribution reduces effective gravity most strongly at low latitudes, influencing the shape of the geoid.[11][17]Mathematical Formulation
Total Geopotential
The total geopotential W is defined as the sum of the gravitational potential V, arising from the Earth's mass distribution, and the centrifugal potential \phi_c, resulting from the planet's rotation. This formulation is expressed mathematically as W = V + \phi_c, where V satisfies Laplace's equation outside the Earth's masses and \phi_c = -\frac{1}{2} \omega^2 (x^2 + y^2) in a Cartesian coordinate system aligned with the rotation axis, with \omega denoting Earth's angular velocity.[18][19] The effective gravity vector \mathbf{g} is the negative gradient of the total geopotential, \mathbf{g} = -\nabla W = -\nabla V - \nabla \phi_c, representing the vector sum of gravitational attraction and centrifugal acceleration observed at Earth's surface. This yields the magnitude of effective gravity varying latitudinally from approximately 9.78 m/s² at the equator to 9.83 m/s² at the poles, primarily due to the equatorial bulge and rotational effects.[18][19] The derivation follows from the conservative nature of both potentials, allowing the total force per unit mass to be obtained directly from the gradient without additional terms for stationary observers. The Coriolis acceleration, which depends on velocity, is omitted in this potential formulation as it vanishes for static measurements and does not contribute to the position-dependent field.[18] Equipotential surfaces of the total geopotential, where W = constant, define the effective gravity field, with the geoid corresponding to the specific value W = W_0 \approx 62636853.4 m²/s², adopted as the conventional reference for mean sea level. These surfaces are nearly ellipsoidal but exhibit undulations of up to ±100 m due to local mass anomalies in the Earth's crust and mantle, influencing height systems in geodesy.[19][20]Normal and Disturbing Potentials
In physical geodesy, the total geopotential W at any point is decomposed into a reference normal potential U and a disturbing potential T, such that T = W - U. This separation allows for the isolation of deviations from an idealized reference field, facilitating the analysis of gravitational anomalies. The normal potential U represents the geopotential of a rotating, level ellipsoid of revolution, serving as the standard model for the Earth's gravity field in precision applications. The normal potential U is defined as the sum of the normal gravitational potential V_n and the centrifugal potential \phi_c, expressed as U = V_n + \phi_c. It is constructed to be constant on the surface of the reference ellipsoid, thereby satisfying Bruns' formula, which relates orthometric heights H to the difference between the constant normal potential value U_0 and the actual geopotential W via H = (U_0 - W)/\gamma, where \gamma is the normal gravity. A widely adopted realization is the Geodetic Reference System 1980 (GRS80), defined by parameters including the equatorial radius a = 6{,}378{,}137 m, the geocentric gravitational constant GM = 3.986005 \times 10^{14} m³ s⁻², the dynamical form factor J_2 = 1.08263 \times 10^{-3}, and the angular velocity \omega = 7.292115 \times 10^{-5} rad s⁻¹. On the ellipsoid surface, the normal gravity \gamma is computed using the Somigliana-Pizzetti formula: \gamma(\phi) = \gamma_e \frac{1 + k \sin^2 \phi}{\sqrt{1 - e^2 \sin^2 \phi}}, where \gamma_e = 9.7803267715 m s⁻² is the equatorial normal gravity, k = 0.001931851353 is a constant, e^2 = 0.00669438002290 is the squared first eccentricity, and \phi is the geodetic latitude. The disturbing potential T encapsulates perturbations due to the Earth's non-ellipsoidal mass distribution, such as topographic features like mountains and ocean trenches, which cause deviations from the reference field. Outside the Earth, T is a harmonic function satisfying Laplace's equation and is commonly expanded in spherical harmonics for global modeling: T(r, \phi, \lambda) = \frac{GM}{r} \sum_{l=2}^{\infty} \sum_{m=0}^{l} \left( \frac{a}{r} \right)^l \bar{P}_{lm}(\sin \phi) \left[ \bar{C}_{lm} \cos(m\lambda) + \bar{S}_{lm} \sin(m\lambda) \right], where r is the geocentric radius, \phi and \lambda are latitude and longitude, a is the reference radius (typically the equatorial radius), \bar{P}_{lm} are fully normalized associated Legendre functions, and \bar{C}_{lm}, \bar{S}_{lm} are the fully normalized spherical harmonic coefficients. This expansion links T to observable quantities, such as the geoid undulation N, approximated on the geoid surface as N = T / g_0, where g_0 \approx 9.80665 m s⁻² is the standard normal gravity value. Additionally, the free-air gravity anomaly \Delta g is approximately related to the vertical derivative of T by \Delta g \approx -\partial T / \partial h, where h is the ellipsoidal height, providing a key connection for gravity data interpretation.Geopotential Number
The geopotential number, often denoted as C, represents the difference between the gravity potential on a reference equipotential surface, such as the geoid (W_0), and the gravity potential at a specific point (W). It is defined as C = W_0 - W.[21] This has units of m²/s², conventionally expressed in geopotential units (gpu), where 1 gpu = 10 m²/s², providing a measure of potential difference that accounts for gravitational variations.[22] Traditionally, the geopotential number is measured through a combination of gravimetry, which determines local gravity g, and geometric leveling, which provides height differences dh along a path. The value is computed via the integration C = \int g \, dh from the geoid to the point, approximated in practice as C \approx \bar{g} \Delta H, where \Delta H is the orthometric height difference and \bar{g} is the average gravity along the plumb line; more precise forms include corrections for gravity anomalies along the leveling line.[23] Absolute values of C can also be obtained using astrogeodetic methods, which measure deflections of the vertical through astronomical observations to directly estimate potential differences relative to the geoid.[24] The geopotential number relates closely to orthometric height H, the physical height above the geoid along the plumb line, via the approximation H \approx C / \bar{g}, where \bar{g} is the average gravity along that line; this corrects for variations in g that affect traditional leveling. The normalized form C / \gamma (with \gamma normal gravity) yields the dynamic height in meters.[25] In modern geodesy, refinements incorporate satellite data, such as GPS for precise ellipsoidal heights and GRACE mission gravimetry for global gravity field models, enabling more accurate computation of C by combining these with local measurements to reduce uncertainties in the integration.[23]Simplified Models
Nonrotating Symmetric Sphere
The nonrotating symmetric sphere represents the simplest idealized model for the Earth's geopotential, treating the planet as a homogeneous, nonrotating body to facilitate introductory calculations in geophysics and geodesy. This approximation ignores density variations, oblateness, and rotational effects, focusing solely on the gravitational contribution under spherical symmetry.[26] Key assumptions include a uniform density \rho \approx 5514 kg/m³ throughout the volume, corresponding to the Earth's total mass M = 5.972 \times 10^{24} kg, and a mean radius r_e = 6371 km. With no rotation, the centrifugal potential vanishes (\phi_c = 0), so the geopotential \Phi equals the gravitational potential V.[26] For points exterior to the sphere (r > r_e), the gravitational potential is identical to that of a point mass at the center: V(r) = -\frac{GM}{r}, where GM = 3.986 \times 10^{14} m³/s² is Earth's standard gravitational parameter and G is the Newtonian gravitational constant.[27] Inside the sphere (r \leq r_e), the potential takes a quadratic form: V(r) = -\frac{GM}{2 r_e^3} (3 r_e^2 - r^2) = -\frac{3 GM}{2 r_e} \left( 1 - \frac{1}{3} \left( \frac{r}{r_e} \right)^2 \right). This ensures continuity of V and its gradient at the surface r = r_e.[28] The gravitational acceleration \mathbf{g} = -\nabla V points radially inward. Outside the sphere, its magnitude is g(r) = \frac{GM}{r^2}, constant on concentric spheres of fixed r. Inside, g(r) = \frac{GM r}{r_e^3} = \frac{4\pi G \rho r}{3}, which vanishes at the center and increases linearly with r, remaining constant on spheres of constant radius.[28] Consequently, the equipotential surfaces are concentric spheres centered at the origin, as V depends only on radial distance.[28] This model's potential satisfies Poisson's equation \nabla^2 V = 4\pi G \rho inside the sphere (and Laplace's equation \nabla^2 V = 0 outside), solved by radial integration from the center under spherical symmetry and boundary matching.[26] At the surface, V(r_e) = -\frac{GM}{r_e} \approx -62.5 MJ/kg, establishing a baseline for geopotential height computations.[28]Rotating Ellipsoid Approximation
The rotating ellipsoid approximation extends the basic spherical model by accounting for Earth's oblateness due to rotation, treating the planet as a fluid body in hydrostatic equilibrium with uniform density and rotating at angular velocity ω ≈ 7.292115 × 10^{-5} rad/s.[29] This model assumes the ellipsoid is an equipotential surface of the total geopotential, balancing gravitational and centrifugal forces to determine the figure of equilibrium.[30] The geometric flattening f, defined as f = (a - b)/a where a is the equatorial semi-major axis and b the polar semi-minor axis, is approximately 1/298.257 for the reference ellipsoid approximating the real Earth.[29] However, for this uniform density model, Clairaut's theorem provides the foundational relation for oblateness, linking the surface flattening f to the centrifugal parameter m = ω² a³ / (GM), where G is the gravitational constant and M the Earth's mass; the first-order solution yields f ≈ (5/4) m ≈ 1/230, which overestimates the observed value of ≈1/298 due to the real Earth's central concentration of mass.[31] The gravitational potential V is approximated to first order in f asV \approx -\frac{GM}{r} \left[ 1 - f \left( \frac{3}{2} \sin^2 \phi - \frac{1}{2} \right) \right],
where r is the geocentric distance and φ the geodetic latitude; this captures the oblate perturbation via the second-degree zonal harmonic, with the centrifugal potential φ_c = -\frac{1}{2} ω^2 \rho^2 \cos^2 \phi added, where ρ is the distance from the rotation axis.[32] The total geopotential U = V + φ_c renders the ellipsoid surface equipotential, approximating the geoid with undulations typically below 1 m relative to this mean figure.[30] This approximation yields a polar-equatorial radius difference of approximately 21 km, with a ≈ 6378 km and b ≈ 6357 km, reflecting the ~0.3% ellipticity driven by rotation.[29] The normal gravity acceleration g on the ellipsoid surface varies latitudinally and is given by the International Gravity Formula:
g(\phi) = g_e \left( 1 + \beta \sin^2 \phi - \alpha \sin^2 2\phi \right),
where g_e ≈ 9.7803 m/s² is the equatorial value, β ≈ 0.0053024, and α ≈ 0.0000058; the form equivalently expresses the increase toward the poles due to both gravitational concentration and centrifugal reduction at the equator.[32] This model provides an intermediate-fidelity baseline for geodetic computations, improving upon the nonrotating sphere by incorporating rotational effects while neglecting internal density variations.[30]