2.1 Earth Geometry: Geoid, Reference Ellipsoids, and Heights

2.1 Earth Geometry: Geoid, Reference Ellipsoids, and Heights

Core Principle: In geospatial science, height is never a singular, unambiguous measurement. Satellite navigation systems measure geometric height above an idealized mathematical ellipsoid ($h$), whereas water flows and civil infrastructure respond strictly to potential energy differences along the Earth's irregular gravity field—measured as orthometric height ($H$) above the geoid. The mathematical bridge between these systems is the geoid undulation ($N$), expressed through the fundamental formula $h = H + N$.

---

1. Foundations of Geodesy and Earth Geometry

Geodesy is the scientific discipline dedicated to measuring and representing the Earth's geometric shape, orientation in space, and gravitational field, as well as monitoring how these properties change over time. Far from being an inert, static sphere, the Earth is an active, continuously deforming body subject to:

• Tectonic Plate Motion: Lithospheric plates gliding over the asthenosphere at rates of 1 to 10 centimeters per year.
• Post-Glacial Isostatic Rebound: Ongoing vertical uplift of crustal masses (up to 1 centimeter per year in regions such as Hudson Bay and Fennoscandia) relieved from the weight of Pleistocene ice sheets.
• Tidal Deformation: Diurnal solid Earth tides caused by gravitational interactions with the Moon and Sun, displacing the crust vertically by up to 30 centimeters.
• Localized Subsidence and Uplift: Anthropogenic fluid extraction (groundwater, hydrocarbons) and seismic events inducing sudden or gradual multi-meter displacements.

Because geospatial professionals integrate satellite observations, aerial imagery, and ground surveys across wide spatial extents, they must model Earth's shape through mathematically tractable yet physically rigorous approximations.

---

2. Three Approximations of Earth Shape

Geodesists and GIS professionals employ three fundamental models to represent the Earth's shape, each serving distinct mathematical, computational, and physical purposes.

   Topographic Surface (Physical Mountains, Valleys, Oceans)
   ---------------------------------------------------------
        \                  /               /
         \   Plumb Line   /  Deflection   /  Ellipsoidal Normal
          \     (g)      /   Angle (θ)   /        (n)
           v            v               v
   ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Geoid (W = W0, Equipotential)
   =============================================== Reference Ellipsoid (GRS80/WGS84)

The Sphere

The sphere is the simplest mathematical model of the Earth, defined by a single mean radius ($R \approx 6{,}371\text{ km}$). It assumes that the distance from the Earth's center to every point on the surface is identical. While mathematically convenient for spherical trigonometry and global projections at small scales (e.g., continental or world maps smaller than $1:5{,}000{,}000$), the sphere fails to account for planetary rotation. Modeling regional or local geospatial data on a sphere introduces positional errors of up to 21 kilometers (approximately 0.33%) between the equator and the poles, making it entirely unsuitable for engineering design, cadastral surveying, or high-accuracy GIS.

The Oblate Spheroid (Reference Ellipsoid)

Because the Earth rotates around its polar axis approximately once every 24 hours, rotational centrifugal force opposes gravitational attraction at the equator while dropping to zero at the rotational poles. Over geological epochs, this centrifugal acceleration deformed the Earth into an oblate spheroid (an ellipsoid of revolution)—flattened at the poles and bulging at the equator. An ellipsoid provides a smooth, mathematically regular, two-parameter geometric surface upon which horizontal coordinates (geodetic latitude $\phi$ and geodetic longitude $\lambda$) can be calculated with sub-millimeter precision.

The Geoid

While the reference ellipsoid provides a clean geometric baseline, it ignores internal mass variations. The Earth's crust, mantle, and core contain heterogeneous density variations (such as dense oceanic basalt, buoyant granitic continental roots, and deep mantle plumes). These localized mass anomalies warp the Earth's gravity field, causing the true physical surface of equal gravitational potential—the geoid—to undulate irregularly above and below the smooth mathematical ellipsoid. The geoid corresponds physically to the shape that global oceans would take if governed solely by gravity and rotation, undisturbed by winds, tides, currents, and salinity gradients, and extended continuously through continental landmasses.

---

3. Geometric Parameters of Reference Ellipsoids

A reference ellipsoid of revolution is generated by rotating a two-dimensional ellipse about its minor (polar) axis. Defining an ellipsoid requires specifying at least two geometric parameters, typically the semi-major axis and the flattening factor.

Mathematical Formulation

1. Semi-Major Axis ($a$): The equatorial radius, measuring the distance from the ellipsoid center to the equator.
2. Semi-Minor Axis ($b$): The polar semi-axis, measuring the distance from the center to either geographic pole.
3. Flattening Factor ($f$): A dimensionless measure of polar compression relative to the equatorial radius:

$f = \frac{a - b}{a}$

Because Earth's polar flattening is slight, $f$ is customarily expressed as an inverse fraction ($1/f \approx 298.257$).

1. First Eccentricity Squared ($e^2$): A fundamental parameter used in coordinate conversion and projection formulas:

$e^2 = \frac{a^2 - b^2}{a^2} = 2f - f^2$

1. Second Eccentricity Squared ($e'^2$):

$e'^2 = \frac{a^2 - b^2}{b^2} = \frac{e^2}{1 - e^2}$

Historical and Modern Reference Ellipsoids

Over the past two centuries, geodesists have derived various ellipsoids to minimize local discrepancies between the ellipsoid surface and the regional geoid. Prior to the satellite era, ellipsoids were calculated from regional terrestrial triangulation networks, resulting in non-geocentric (locally oriented) ellipsoids such as Clarke 1866. Modern satellite tracking and orbital analysis enabled the derivation of geocentric ellipsoids whose origin coincides with Earth's center of mass.

Technical Observation: The GRS80 (Geodetic Reference System 1980) and WGS84 (World Geodetic System 1984) ellipsoids share the identical semi-major axis ($a = 6{,}378{,}137.0\text{ m}$). Their flattening values differ by only $0.00000016$ due to slight differences in defining physical constants, causing their polar semi-minor axes ($b$) to differ by less than $0.1\text{ millimeter}$—a negligible geometric difference in GIS workflows.

---

4. The Geoid: Physics of Equipotential Surfaces

The geoid is defined as an equipotential surface of the Earth's gravity field that best fits global Mean Sea Level in a least-squares sense. An equipotential surface is a surface across which the gravitational potential energy ($W$) is constant everywhere:

$W(x, y, z) = V(x, y, z) + \Phi(x, y, z) = W_0$

Where:

• $V(x, y, z)$ is the gravitational potential resulting from mass attraction of Earth's materials.
• $\Phi(x, y, z)$ is the centrifugal potential generated by planetary rotation.
• $W_0$ is the specific geopotential constant chosen to define the geoid (adopted as $62{,}636{,}856.0\text{ m}^2\text{s}^{-2}$ in modern standards).

The Plumb Line and Deflection of the Vertical

At any point on Earth, the gravity vector $\mathbf{g}$ acts perpendicular (orthogonal) to the equipotential surface passing through that point. A freely hanging surveyor's plumb line follows this gravity vector. Because mass distribution within the crust and mantle is heterogeneous, the direction of gravity curves slightly with depth and diverges from the mathematical normal to the reference ellipsoid.

• Ellipsoidal Normal: A straight line perpendicular to the geometric reference ellipsoid at a given latitude and longitude.
• Plumb Line: A physical, curved trajectory tangent to the gravity vector $\mathbf{g}$ at every point along its path, perpendicular to the geoid.
• Deflection of the Vertical ($\theta$): The angular divergence between the ellipsoidal normal and the true gravity plumb line. Decomposed into a north-south meridian component ($\xi$) and an east-west prime vertical component ($\eta$):

$\theta = \sqrt{\xi^2 + \eta^2}$

In mountainous areas with massive terrain variations (e.g., the Andes or Colorado Rockies), the deflection of the vertical can reach 30 to 60 arcseconds. This angular discrepancy causes optical leveling instruments (which orient to local gravity via spirit vials) to diverge measurably from pure geometric satellite vectors.

Geoid Undulation ($N$)

The separation distance between the geoid and the reference ellipsoid, measured along the ellipsoidal normal, is termed the geoid undulation or geoid height ($N$):

• If the geoid is above the ellipsoid at a given coordinate, $N$ is positive ($N > 0$).
• If the geoid is below the ellipsoid at a given coordinate, $N$ is negative ($N < 0$).

Across the conterminous United States (CONUS), relative to the GRS80 and WGS84 ellipsoids, the geoid lies entirely beneath the ellipsoid. Values range from approximately $-8\text{ meters}$ in southern Florida to $-53\text{ meters}$ in the northern Rocky Mountains of Idaho and Montana.

---

5. Height Systems and the Fundamental Relationship

Understanding the three height components and their algebraic relationship is one of the most frequently tested concepts on the GISP exam.

                                    Topographic Ground Surface
                                              * [Point P]
                                             /|
                                            / |
                                           /  |
                                          /   |
                      Orthometric Height /    | Ellipsoidal Height
                             (H)        /     |       (h)
                                       /      |
                                      /       |
             Geoid (Mean Sea Level)  ~        |
            ~~~~~~~~~~~~~~~~~~~~~~~~~*        |
                                     |        |
                Geoid Undulation (N) |        |
                 (Negative in CONUS) |        |
                                     v        v
            ==================================*====== Reference Ellipsoid

Defining the Three Heights

1. Ellipsoidal Height ($h$):

   • The geometric distance measured along the ellipsoidal normal from the surface of the reference ellipsoid to the point of interest on the physical surface of the Earth.
   • Generated directly by Global Navigation Satellite System (GNSS) receivers (GPS, GLONASS, Galileo).
   • Purely mathematical; has no inherent physical relationship to gravity, water flow, or hydrostatic pressure.

2. Orthometric Height ($H$):

   • The physical distance measured along the curved plumb line from the geoid to the point of interest on the Earth's surface.
   • Represents true elevation above Mean Sea Level (MSL).
   • Governs fluid mechanics, gravity drainage, runoff flow directions, and civil engineering design.

3. Geoid Undulation / Geoid Height ($N$):

   • The separation distance from the reference ellipsoid to the geoid, measured along the ellipsoidal normal.
   • Modeled via numerical geoid models (e.g., GEOID12B, GEOID18).

The Fundamental Vertical Formula

$h = H + N$

Rearranging to solve for orthometric height (elevation):

$H = h - N$

Rearranging to solve for geoid undulation:

$N = h - H$

Step-by-Step Worked Numerical Examples

Example 1: Calculating Orthometric Elevation in CONUS (Negative $N$)

A civil engineering GIS analyst establishes a control point in Denver, Colorado using dual-frequency GNSS. The receiver outputs an ellipsoidal height $h = 1{,}609.34\text{ m}$. Querying NGS GEOID18 at the station yields a geoid undulation $N = -20.50\text{ m}$. What is the orthometric height ($H$) of the station?

$H = h - N = 1{,}609.34\text{ m} - (-20.50\text{ m}) = 1{,}609.34\text{ m} + 20.50\text{ m} = 1{,}629.84\text{ meters}$

Exam Trap Alert: Notice the double negative! Because $N$ is negative across North America, subtracting a negative number results in an increase. The orthometric elevation ($1{,}629.84\text{ m}$) is greater than the ellipsoidal height ($1{,}609.34\text{ m}$). A common error is calculating $1{,}609.34 - 20.50 = 1{,}588.84\text{ m}$, introducing a 41-meter error.

Example 2: Station with Positive Geoid Undulation ($N > 0$)

A geodetic station located on an oceanic island in the North Atlantic records an ellipsoidal height $h = 45.20\text{ m}$. The regional gravimetric geoid model reports $N = +12.30\text{ m}$. Calculate the orthometric height ($H$):

$H = h - N = 45.20\text{ m} - (+12.30\text{ m}) = 32.90\text{ meters}$

Here, because the geoid swells above the reference ellipsoid ($N > 0$), the orthometric height above Mean Sea Level ($32.90\text{ m}$) is lower than the geometric ellipsoidal height ($45.20\text{ m}$).

---

6. Benchmark Leveling vs. Modern GNSS Elevation Modeling

How elevation data is captured dictates the accuracy and vertical system of the resulting GIS layers.

Differential Spirit Leveling

Classical differential leveling involves an optical or digital leveling instrument positioned midway between two graduated rods. The instrument's optical sight-line is leveled perpendicular to the local gravity vector via a compensator or spirit bubble. Reading backsights and foresights determines the physical elevation difference ($\Delta H$) between adjacent benchmarks along the curved plumb line.

• Strengths: Directly measures orthometric height differences; achieves first-order millimeter precision over kilometer-length baselines; inherently honors gravity and hydraulic flow.
• Limitations: Extremely labor-intensive and slow; requires contiguous lines of sight along road networks; subject to systematic accumulation of refraction and monument settling errors; vulnerable to physical monument destruction.

GNSS + Geoid Modeling

Modern positioning uses multi-constellation GNSS receivers to determine three-dimensional Cartesian vectors ($X, Y, Z$) relative to the geocenter, which convert mathematically to geodetic latitude, longitude, and ellipsoidal height ($h$). To convert geometric $h$ into usable orthometric elevation ($H$), geodesists apply a hybrid geoid model:

$H_{\text{derived}} = h_{\text{GNSS}} - N_{\text{model}}$

Gravimetric vs. Hybrid Geoid Models

• Pure Gravimetric Geoid Models (e.g., USGG2012): Created solely from terrestrial, airborne, and satellite gravity measurements (such as GRACE and GOCE missions). Represents the true physical geopotential surface, but may deviate by tens of centimeters from published leveling benchmarks due to residual long-wavelength gravity errors.
• Hybrid Geoid Models (e.g., GEOID12B, GEOID18): Produced by warping a gravimetric geoid model to fit thousands of high-precision "GPS-on-benchmarks"—physical monuments that have both first-order leveled orthometric heights in NAVD88 and high-accuracy GNSS-derived ellipsoidal heights in NAD83. Hybrid models are specifically engineered for surveying and GIS workflows to allow GNSS receivers to directly reproduce published NAVD88 orthometric heights within 1 to 2 centimeters across CONUS.

---

7. Practical Geospatial Applications and Scenarios

Scenario 1: Stormwater Drainage and Culvert Inversion Failures

A municipal GIS team is modeling a gravity-fed storm sewer network designed with a minimum slope of $0.20%$ ($2\text{ mm}$ drop per meter) across a 1.5-kilometer drainage corridor. A field crew captures the invert elevations of upstream and downstream catch basins using standalone consumer-grade GNSS without applying a geoid model, recording raw ellipsoidal heights ($h$).

Across this 1.5-kilometer stretch, the local geoid undulation slopes by $18\text{ centimeters}$. By relying on ellipsoidal heights ($h$) rather than orthometric heights ($H$), the engineering team introduces an artificial vertical bias greater than the total designed hydraulic drop of the pipeline. The resulting culvert is constructed with an adverse (backward) grade, causing stormwater to pool and flood surrounding neighborhoods during severe precipitation.

Scenario 2: FEMA Floodplain Delineation and Base Flood Elevation (BFE)

In accordance with National Flood Insurance Program (NFIP) guidelines, FEMA Base Flood Elevations (BFEs) are defined exclusively as orthometric elevations referenced to the North American Vertical Datum of 1988 (NAVD88). If an airborne LiDAR contractor delivers digital elevation models (DEMs) referenced to raw ellipsoidal heights ($h$) without performing geoid conversion ($N$), the entire terrain surface will sit 20 to 35 meters below or above its true physical relationship to riverine water levels across North America, corrupting regulatory Special Flood Hazard Area (SFHA) boundaries.

---

8. GISP Exam Traps & Pitfalls

• The Sign Trap: Always use the signed geoid undulation in $H = h - N$. Much of the conterminous United States has negative $N$, but the sign and magnitude depend on location and geoid model; never replace $N$ with its absolute value.
• Water Flows by Gravity, Not Ellipsoid Geometry: Water can physically flow toward a point that has a higher ellipsoidal height if the geoid undulation at that location drops even faster. Fluid flow is strictly governed by potential energy ($H$), never by geometric distance from the ellipsoid ($h$).
• The Geoid is Not a Smooth Mathematical Surface: You cannot represent the geoid with an algebraic polynomial or simple geometric formula. It is an equipotential surface determined by complex, uneven mass distributions within the Earth's crust and mantle, requiring gridded interpolation lookup files (e.g., GEOID18 .bin grids).
• Equating Mean Sea Level Directly with the Geoid: In casual discussion, the geoid is often called "Mean Sea Level." On the GISP exam, remember that true local Mean Sea Level deviates from the geoid by up to 1 to 2 meters globally due to Dynamic Ocean Topography (DOT) caused by persistent oceanic currents (such as the Gulf Stream), salinity gradients, prevailing winds, and water temperature variations.