What is Spherical Trigonometry, Geodesy, and Great-Circle Navigation?
Mathematical Foundation
Laws & Principles
- Why Flat-Map Straight Lines Fail: The Mercator projection preserves angles (conformal) but severely distorts distances and areas at high latitudes. A straight line on a Mercator map represents a rhumb line (constant compass bearing), not a great circle. A flight from Los Angeles to Tokyo that follows a constant westward heading (a rhumb line) is approximately 400 km longer than the great-circle path that curves far north over Alaska. Every fuel tankload saved by flying the correct route is worth tens of thousands of dollars.
- The Polar Routing Effect: Great circles connecting cities in the Northern Hemisphere naturally arc toward the North Pole. JFK→LHR curves over Nova Scotia and Iceland. LAX→NRT curves over Alaska. This is geometrically correct — the pole is the 'short side' of the sphere. In winter, airlines sometimes fly slightly south of the optimal great circle to avoid strong jet stream headwinds, accepting a slightly longer path for better fuel economy.
- Haversine vs. Vincenty Formula: The Haversine formula treats Earth as a perfect sphere (R = 6,371 km). For most aviation and navigation purposes, this is accurate to within 0.3%. For geodetic surveying requiring millimeter precision, the Vincenty formula accounts for Earth's oblateness (equatorial radius = 6,378 km, polar radius = 6,357 km, flattening = 1/298.257). The Haversine formula's main advantage is numerical stability — it avoids the floating-point issues that plague the simpler spherical law of cosines near small distances.
- Nautical Miles and the Arc-Minute Convention: One nautical mile (NM) is defined as exactly one minute of arc of latitude (1/60th of a degree of latitude) along a meridian. Earth's circumference = 360° × 60 NM/° = 21,600 NM = 40,003 km. This makes nautical miles uniquely convenient for great-circle navigation because you can read distance directly from a chart's latitude scale. No conversion factor is needed.
- Antipodal Points: Two points are antipodal if they are exactly opposite each other on the globe (e.g., Madrid and Auckland, or the North Pole and South Pole). The great-circle distance between antipodal points is exactly π × R = 20,015 km (half Earth's circumference). Any path between them is equally a minimum — the concept of 'direction' to an antipodal point is undefined.
Step-by-Step Example Walkthrough
" Calculate the great-circle distance from JFK (40.6413°N, 73.7781°W) to LHR Heathrow (51.4700°N, 0.4543°W). "
- 1. Convert to radians: φ₁ = 40.6413° × π/180 = 0.70919 rad, λ₁ = −73.7781° × π/180 = −1.28741 rad, φ₂ = 51.4700° × π/180 = 0.89844 rad, λ₂ = −0.4543° × π/180 = −0.00793 rad.
- 2. Differences: Δφ = 0.89844 − 0.70919 = 0.18925 rad, Δλ = −0.00793 − (−1.28741) = 1.27948 rad.
- 3. Compute a: sin²(0.18925/2) = sin²(0.09463) = 0.08949² = 0.00894. cos(0.70919) × cos(0.89844) = 0.75988 × 0.62756 = 0.47690. sin²(1.27948/2) = sin²(0.63974) = 0.59620² = 0.35545. a = 0.00894 + 0.47690 × 0.35545 = 0.00894 + 0.16954 = 0.17848.
- 4. c = 2 × atan2(√0.17848, √0.82152) = 2 × atan2(0.42246, 0.90638) = 2 × 0.43614 = 0.87228 rad.
- 5. d = 6371 × 0.87228 = 5,556 km = 3,001 NM = 3,451 statute miles.
- 6. Initial bearing: θ = atan2(sin(1.27948)×cos(0.89844), cos(0.70919)×sin(0.89844) − sin(0.70919)×cos(0.89844)×cos(1.27948)) = atan2(0.95631×0.62756, 0.75988×0.77856 − 0.64938×0.62756×0.29025) ≈ 50.8° — initial heading North-Northeast.