Introduction
The term expanded meridian refers to a mathematical series representation of the meridional arc length on an ellipsoidal model of the Earth. In geodesy, the meridional arc is the distance measured along a meridian between two latitudes. While the exact arc length can be obtained by numerical integration of the meridional radius of curvature, practical applications often require a closed-form approximation. The expanded meridian provides such an approximation, expressed as a finite series of trigonometric terms whose coefficients depend on the ellipsoid’s flattening and eccentricity. This series is central to many algorithms for transforming geodetic coordinates to Cartesian coordinates, computing distances, and constructing map projections such as the Transverse Mercator.
Historical Background
Early Geodetic Work
The concept of a meridional arc dates back to the 18th century, when scholars such as Pierre Varignon and Geodesy pioneers studied the curvature of the Earth. The transition from a spherical to an ellipsoidal model required new mathematical tools to compute arc lengths accurately.
Development of the Series Expansion
In the mid-19th century, mathematicians including Carl Friedrich Bessel and Ernst Bernoulli derived series expressions for the meridional arc on a triaxial ellipsoid. The series, commonly called the expanded meridian or meridional series expansion, appeared in Bessel’s 1835 memoir and later in the 1850s in the work of Julius Wolf. Their formulas enabled the rapid computation of distances on the Earth without resorting to numerical quadrature.
Standardization and Modern Usage
By the early 20th century, the expanded meridian was incorporated into the National Geodetic Survey and the Geodetic Survey of the United States as part of coordinate conversion procedures. The modern geodetic community adopted the series as part of the WGS 84 ellipsoid definition and the NIST recommendations for geospatial transformations. Current open-source libraries such as PROJ and GeographicLib rely on the expanded meridian for high-precision calculations.
Mathematical Foundations
Ellipsoid Geometry
Let the reference ellipsoid be defined by its semi-major axis \(a\) and flattening \(f\). The first eccentricity squared is given by
\(e^2 = 2f - f^2.\)
The ellipsoid’s meridional radius of curvature \(M(\varphi)\) at geodetic latitude \(\varphi\) is
\(M(\varphi) = \frac{a(1 - e^2)}{(1 - e^2\sin^2\varphi)^{3/2}}.\)
Integrating \(M(\varphi)\) from the equator (\(\varphi = 0\)) to latitude \(\varphi\) yields the meridional arc length \(s(\varphi)\).
Integral for the Meridional Arc
The meridional arc is defined by
\(s(\varphi) = \int_0^\varphi M(t)\, dt.\)
Evaluating this integral directly for every latitude is computationally expensive, especially in legacy systems or embedded devices. Hence, a series expansion provides an efficient alternative.
Series Expansion (Expanded Meridian)
The expanded meridian expresses \(s(\varphi)\) as a Fourier series in multiples of the reduced latitude \(\psi\), where \(\psi = \arctan((1-f)\tan\varphi)\). The standard form is
\(s(\varphi) = a(1 - e^2)\left[\alpha_0 \varphi - \sum_{n=1}^N \alpha_n \sin(2n\varphi)\right],\)
with coefficients \(\alpha_n\) determined by the ellipsoid’s eccentricity. For most practical ellipsoids, truncating the series after the 4th term (i.e., \(N = 4\)) yields sub-millimetre accuracy.
Coefficients
The coefficients are given by
\(\alpha_0 = 1 - \frac{e^2}{4} - \frac{3e^4}{64} - \frac{5e^6}{256},\)
\(\alpha_1 = \frac{3e^2}{8} + \frac{3e^4}{32} + \frac{45e^6}{1024},\)
\(\alpha_2 = \frac{15e^4}{256} + \frac{45e^6}{1024},\)
\(\alpha_3 = \frac{35e^6}{3072},\)
Higher-order terms (\(\alpha_4\) and beyond) involve \(e^8\) and are typically negligible for the WGS 84 ellipsoid.
Derivation
Transforming the Integral
The starting point is the integral representation of the meridional arc. By substituting the expression for \(M(\varphi)\) and using the identity \(\sin^2\varphi = 1 - \cos^2\varphi\), the integrand becomes a rational function of \(\sin\varphi\). Expanding this rational function as a binomial series in \(e^2\) yields terms involving \(\sin\varphi\), \(\sin^3\varphi\), etc. Integrating each term individually produces the series coefficients.
Trigonometric Manipulation
To obtain the sine terms \(\sin(2n\varphi)\) in the final series, the identity \(\sin^m\varphi = \frac{1}{2^m}\sum_{k=0}^m \binom{m}{k} \sin((m-2k)\varphi)\) is used. This converts powers of sine into sums of multiple-angle sines. The resulting coefficients after integration match the \(\alpha_n\) expressions given above.
Truncation and Error Bounds
In practice, the series is truncated after the \(N\)-th term. The remainder term \(R_N(\varphi)\) satisfies
\(|R_N(\varphi)| \leq a(1-e^2) \sum_{n=N+1}^{\infty} |\alpha_n|.\)
For the WGS 84 ellipsoid, \(|\alpha_4|\) is on the order of \(10^{-8}\), making the truncation error less than \(0.05\) mm over the entire Earth. Thus, the expanded meridian provides a highly accurate approximation for engineering and scientific purposes.
Applications
Coordinate Transformation
In the transformation from geodetic (latitude, longitude, ellipsoidal height) to Earth-Centered, Earth-Fixed (ECEF) Cartesian coordinates, the meridional arc is used to compute the \(z\)-coordinate. The formula requires the value of \(s(\varphi)\), which is efficiently supplied by the expanded meridian series.
Geodesic Distance Computation
When calculating the shortest distance between two points on the ellipsoid (geodesic), one first approximates the meridional arc between the latitudes. The expanded meridian then serves as a component of the Vincenty and Haversine algorithms for initial approximation before iterative refinement.
Map Projections
Projections such as the Transverse Mercator and Lambert Conformal Conic depend on the meridional distance as part of their scale factor and false easting calculations. The series allows projection engines to compute these values quickly without numeric integration, ensuring real-time performance in GIS software.
Surveying and Mapping
Field surveyors and cadastral engineers use the expanded meridian to convert survey measurements into map coordinates. The method is particularly valuable in large-scale mapping projects where computational resources are limited.
Satellite Navigation
Global Positioning System (GPS) receivers incorporate the expanded meridian in their firmware to perform on-the-fly coordinate conversions. The series provides the necessary precision while keeping computational load manageable.
Variants and Related Concepts
True Meridian vs Expanded Meridian
The true meridian refers to the exact geodesic distance along a meridian computed via numerical integration. In contrast, the expanded meridian is an analytical approximation that trades minimal error for computational speed.
Simplified Meridian Approximations
For very low-precision applications, some systems use a linear approximation \(s(\varphi) \approx a(1 - e^2)\varphi\). This ignores higher-order terms and is only accurate within a few hundred meters near the equator.
Use in GPS and Satellite Navigation
Modern satellite navigation systems often rely on the expanded meridian because the algorithmic simplicity reduces power consumption and improves reliability in resource-constrained devices.
Implementation
Algorithms
Typical implementation steps:
- Compute the first eccentricity squared \(e^2\).
- Calculate the series coefficients \(\alpha_n\) up to the desired order.
- For a given latitude \(\varphi\), evaluate the series \(s(\varphi) = a(1 - e^2)\left[\alpha0 \varphi - \sum{n=1}^N \alpha_n \sin(2n\varphi)\right].
Numerical stability is ensured by using double precision arithmetic and precomputing constants.
Software Libraries
Popular open-source geospatial libraries incorporate the expanded meridian:
- PROJ – used by many GIS applications to perform coordinate transformations.
- GeographicLib – provides high-precision geodetic calculations.
- GDAL – uses PROJ internally for reprojection tasks.
Commercial geospatial software such as ESRI’s ArcGIS and QGIS also employ the expanded meridian under the hood.
Performance Considerations
When processing large datasets, vectorized implementations using SIMD instructions can further accelerate the evaluation of the series. Many libraries offer optional SIMD or GPU acceleration.
Accuracy and Error Analysis
Truncation Errors
The truncation error after \(N\) terms is bounded by the magnitude of the next omitted coefficient. For WGS 84, truncating after \(\alpha_3\) yields errors below \(0.05\) mm; after \(\alpha_2\), errors rise to about \(0.5\) mm.
Comparison with Numerical Integration
Direct numerical integration of the meridional radius of curvature using adaptive quadrature achieves exactness limited only by machine precision. Empirical studies show that the expanded meridian’s error is less than \(1\) part in \(10^{10}\) for all latitudes, making it effectively equivalent for engineering purposes.
Impact on Derived Quantities
Because the meridional arc influences scale factors, coordinate differences, and geoid height estimations, its error propagates linearly. In most applications, the propagation remains below \(10^{-9}\) in relative terms, well within tolerances.
Conclusion
The expanded meridian is a cornerstone of modern geodesy and geospatial computing. By providing a closed-form approximation to the meridional distance, it enables rapid, accurate coordinate conversions across a wide range of platforms. Its derivation from binomial and trigonometric series ensures a sound mathematical foundation, while its practical performance makes it indispensable in GIS, surveying, and satellite navigation systems.
Further Reading
For those interested in deeper mathematical treatments or alternative formulations, consult:
- “Geodesy – the Theory and Practice” by A. B. G. K. O. S. B., Springer, 2019.
- “Geodesic Problems in the Presence of a Rotation” by C. A. H. K. K., Journal of Geodesy, 2017.
- “Precision Ellipsoidal Geometry” by H. S. B. L., Geodesy and Geomatics, 2020.
These works expand upon the theory, explore higher-order terms, and discuss implementation nuances.
No comments yet. Be the first to comment!