Meridian arc

From Wikipedia, the free encyclopedia
Jump to: navigation, search

In geodesy, a meridian arc measurement is the distance between two points with the same longitude. Two or more such determinations at different locations then specify the shape of the reference ellipsoid which best approximates the shape of the geoid. This process is called the determination of the Figure of the Earth. The earliest determinations of the size of a spherical Earth required a single arc. The latest determinations use astro-geodetic measurements and the methods of satellite geodesy to determine the reference ellipsoids.

The Earth as a sphere[edit]

Main article: Earth radius

Early estimations of Earth's radius are recorded from Egypt in 240 BC, and from Baghdad caliphs in the 9th century, but it was the Alexandrian scientist Eratosthenes who first calculated the circumference a reasonably good approximate value for the radius. He knew that on the summer solstice at local noon the sun goes through the zenith in the ancient Egyptian city of Syene (Assuan). He also knew from his own measurements that, at the same moment in his hometown of Alexandria, the zenith distance was 1/50 of a full circle (7.2°).

Assuming that Alexandria was due north of Syene, Eratosthenes concluded that the distance between Alexandria and Syene must be 1/50 of Earth's circumference. Using data from caravan travels, he estimated the distance to be 5000 stadia (about 500 nautical miles)—which implies a circumference of 252,000 stadia. Assuming the Attic stadion (185 m) this corresponds to 46,620 km, or 16% too great. However, if Eratosthenes used the Egyptian stadion (157.5 m) his measurement turns out to be 39,690 km, an error of only 1%. Syene is not precisely on the Tropic of Cancer and not directly south of Alexandria. The sun appears as a disk of 0.5°, and an estimate of the overland distance traveling along the Nile or through the desert couldn't be more accurate than about 10%.

Eratosthenes' estimation of Earth’s size was accepted for nearly two thousand years. A similar method was used by Posidonius about 150 years later, and slightly better results were calculated in AD 827 by the Gradmessung[citation needed] of the Caliph al-Ma'mun.

The Earth as an ellipsoid[edit]

Comment: early literature uses the term oblate spheroid to describe a sphere "squashed at the poles". Modern literature uses the term "ellipsoid of revolution" although the qualifying words "of revolution" are usually dropped. An ellipsoid which is not an ellipsoid of revolution is called a tri-axial ellipsoid. Spheroid and ellipsoid are used interchangeably in this article.

The eighteenth century[edit]

In 1687 Newton had published in the Principia a proof that the earth was an oblate spheroid (of flattening equal to 1/230).[1] This was disputed by some, but not all, French scientists. A meridian arc of Picard was extended to a longer arc by Cassini (J.D.) over the period 1684–1718. The arc was measured with at least three latitude determinations, so they were able to deduce mean curvatures for the northern and southern halves of the arc, allowing a determination of the overall shape. The results indicated that the Earth was a prolate spheroid (with an equatorial radius less than the polar radius). (The history of the meridian arc from 1600 to 1880 is fully covered in the first chapter of Geodesy by Alexander Ross Clarke.[2]).

To resolve the issue, the French Academy of Sciences (1735) proposed expeditions to Peru (Bouguer, Louis Godin, de La Condamine, Antonio de Ulloa, Jorge Juan ) and Lappland (Maupertuis, Clairaut, Camus, Le Monnier, Abbe Outhier, Celsius). (The expedition to Peru is described on the page French Geodesic Mission and that to Lappland is described on the page Torne Valley.) The resulting measurements at equatorial and polar latitudes confirmed that the earth was best modelled by an oblate spheroid, supporting Newton.

By the end of the century the French arc had been remeasured and extended from Dunkirk to the Mediterranean. (By Delambre). It was divided into five parts by four intermediate determinations of latitude. By combining the measurements together with those for the arc of Peru, ellipsoid shape parameters were determined and the distance between the equator and pole along the Paris Meridian was calculated as 5130762 toise (as specified by the standard toise bar in Paris). Defining this distance as exactly 10,000,000 m led to the construction of a new standard metre bar as 0.5130762 toise. (See Clarke,[2] pp18–22).

The nineteenth and twentieth centuries[edit]

In the 19th century, many astronomers and geodesists were engaged in detailed studies of the Earth's curvature along different meridian arcs. The analyses resulted in a great many model ellipsoids such as Plessis 1817, Airy 1830, Bessel 1830, Everest 1830, and Clarke 1866. A comprehensive list of ellipsoids is given under Earth ellipsoid.

Meridian distance on the ellipsoid[edit]

The determination of the meridian distance, that is the distance from the equator to a point at a latitude \varphi on the ellipsoid is an important problem in the theory of map projections, particularly the Transverse Mercator projection. Ellipsoids are normally specified in terms of the parameters defined above, a, b, f,  but in theoretical work it is useful to define extra parameters, particularly the eccentricity, e, and the third flattening n. Only two of these parameters are independent and there are many relations between them:


\begin{align}
 f&=\frac{a-b}{a}, \qquad  e^2=f(2-f), \qquad n=\frac{a-b}{a+b}=\frac{f}{2-f},\\
b&=a(1-f)=a(1-e^2)^{1/2},\qquad  e^2=\frac{4n}{(1+n)^2}.
\end{align}


The meridian radius of curvature can be shown[3] to be equal to

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

so that the arc length of an infinitesimal element of the meridian is dm = M(\varphi) \, d\varphi (with \varphi in radians). Therefore the meridian distance from the equator to latitude \varphi is


\begin{align}
m(\varphi) &=\int_0^\varphi M(\varphi) \, d\varphi
= a(1 - e^2)\int_0^\varphi \bigl(1 - e^2 \sin^2 \varphi \bigr)^{-3/2} \, d\varphi.
\end{align}

The distance formula is simpler when written in terms of the parametric latitude,


m(\varphi) = b\int_0^\beta
\sqrt{1 + e'^2\sin^2\beta}\,d\beta,

where \tan\beta = (1-f)\tan\varphi and e'^2 = e^2/(1-e^2). The distance from the equator to the pole, the quarter meridian, is


m_p = m(\pi/2).\,

Even though latitude is normally convined to the range [-\frac12\pi,\frac12\pi], all the formulas given here apply to measuring distance around the complete meridian ellipse (including the anti-meridian). Thus the ranges of \varphi, \beta, and the rectifying latitude \mu, are unrestricted.

Relation to elliptic integrals[edit]

The above integral is related to a special case of an incomplete elliptic integral of the third kind. In the notation of the online NIST handbook[4] (Section 19.2(ii)),


m(\varphi)=a\big(1-e^2\big)\,\Pi(\varphi,e^2,e).

It may also be written in terms of incomplete elliptic integrals of the second kind (See the NIST handbook Section 19.6(iv)),


\begin{align}
m(\varphi) &= a\biggl(E(\varphi,e)-\frac{e^2\sin\varphi\cos\varphi}{\sqrt{1-e^2\sin^2\varphi}}\biggr) \\
&= b E(\beta, ie').
\end{align}

The quarter meridian can be expressed in terms of the complete elliptic integral of the second kind,

m_p=aE(e)=bE(ie').

The calculation (to arbitrary precision) of the elliptic integrals and approximations are also discussed in the NIST handbook. These functions are also implemented in computer algebra programs such as Mathematica[5] and Maxima.[6]

Series expansions[edit]

The above integral may be expressed as an infinite truncated series by expanding the integrand in a Taylor series, performing the resulting integrals term by term, and expressing the result as a trigonometric series. In 1755, Euler[7] derived an expansion in the third eccentricity squared.

Expansions in the eccentricity (e)[edit]

Delambre in 1799[8] derived a widely used expansion on the eccentricity squared,


  m(\varphi)=\frac{b^2}a\bigl(D_0\varphi+D_2\sin 2\varphi+D_4\sin4\varphi
          +D_6\sin6\varphi+D_8\sin8\varphi+\cdots\bigr),

where


\begin{align}
D_0 &= \textstyle  1
        + \frac{3}{4} e^2
        + \frac{45}{64} e^4
        + \frac{175}{256} e^6
        + \frac{11025}{16384} e^8 + \cdots, \\
D_2 &= \textstyle - \frac{3}{8} e^2
        - \frac{15}{32} e^4
        - \frac{525}{1024} e^6
        - \frac{2205}{4096} e^8 - \cdots, \\
D_4 &= \textstyle \frac{15}{256} e^4
        + \frac{105}{1024} e^6
        + \frac{2205}{16384} e^8 + \cdots, \\
D_6 &= \textstyle - \frac{35}{3072} e^6
        - \frac{105}{4096} e^8 - \cdots, \\
D_8 &= \textstyle \frac{315}{131072} e^8 + \cdots.
\end{align}

Rapp[9] gives a detailed derivation of this result.

Expansions in the third flattening (n)[edit]

Series with considerably faster convergence can be obtained by expanding in terms of the third flattening, n instead of the eccentricity. This are related by


e^2 = \frac{4n}{(1+n)^2}.

In 1837, Bessel obtained one such series,[10] which was put into a simpler form by Helmert,[11][12]


m(\varphi)=\frac{a+b}2\bigl(H_0\varphi+H_2\sin 2\varphi+H_4\sin4\varphi+H_6\sin6\varphi+H_8\sin8\varphi+\cdots\bigr),

with


\begin{align}
H_0 &= \textstyle 1
        + \frac{1}{4} n^2
        + \frac{1}{64} n^4 + \cdots, \\
H_2 &= \textstyle - \frac{3}{2} n
        + \frac{3}{16} n^3 + \cdots,&
H_6 &= \textstyle - \frac{35}{48} n^3 + \cdots, \\
H_4 &= \textstyle \frac{15}{16} n^2
        - \frac{15}{64} n^4 - \cdots,\qquad&
H_8 &= \textstyle \frac{315}{512} n^4 - \cdots.
\end{align}

Because n changes sign when a and b are interchanged and because the initial factor, \frac12(a+b), is constant under this interchange, half the terms in the expansions of H_{2k} vanish.

The series can be expressed with either a or b as the initial factor by writing, for example, \frac12(a+b) = a/(1+n) = a(1-n+n^2-n^3+n^4-\cdots), and expanding the result as a series in n. Even though this results in more slowly converging series, such series are used in the specification for the transverse Mercator projection by the National Geospatial Intelligence Agency[13] and the Ordnance Survey of Great Britain.[14]

Series in terms of the parametric latitude[edit]

Bessel [15] derived an expansion of the meridian distance in terms of the parametric latitude \beta in connection with his work on geodesics,


m(\varphi)=\frac{a+b}2\bigl(B_0\beta+B_2\sin 2\beta+B_4\sin4\beta+B_6\sin6\beta+B_8\sin8\beta+\cdots\bigr),

with


\begin{align}
B_0 &= \textstyle  1
        + \frac{1}{4} n^2
        + \frac{1}{64} n^4 + \cdots = H_0,\\
B_2 &= \textstyle - \frac{1}{2} n
        + \frac{1}{16} n^3 + \cdots, &
B_6 &= \textstyle - \frac{1}{48} n^3 + \cdots, \\
B_4 &= \textstyle - \frac{1}{16} n^2
        + \frac{1}{64} n^4 + \cdots, \qquad&
B_8 &= \textstyle - \frac{5}{512} n^4 + \cdots.
\end{align}

Because this series provides an expansion for the elliptic integral of the second kind, it can be used to write the arc length in terms of the geographic latitude as


\begin{align}
m(\varphi)&=\frac{a+b}2\biggl(B_0\varphi-B_2\sin 2\varphi+B_4\sin4\varphi-B_6\sin6\varphi+B_8\sin8\varphi-\cdots \\
&\qquad-\frac{2n \sin2\varphi}{\sqrt{1 + 2n \cos2\varphi + n^2}}
 \biggr) .
\end{align}

Generalized series[edit]

The above series, to eighth order in eccentricity or fourth order in third flattening, provide millimetre accuracy. With the aid of symbolic algebra systems, they can easily be extended to sixth order in the third flattening which provides full double precision accuracy for terrestrial applications.

Delambre[8] and Bessel[15] both wrote their series in a form that allows them to be generalized to arbitrary order. The coefficients in Bessel's series can expressed particularly simply


B_{2k} =
\begin{cases} c_0, & \text{if }k = 0, \\
c_k/k, & \text{if } k > 0,
\end{cases}

where


c_k = \sum_{j=0}^\infty
\frac{(2j-3)!!\, (2j+2k-3)!!}{(2j)!!\, (2j+2k)!!} n^{k+2j}

and k!! is the double factorial, extended to negative values via the recursion relation: (-1)!! = 1 and (-3)!! = -1.

The coefficients in Helmert's series can similarly be expressed generally by


H_{2k} = (-1)^k (1-2k)(1+2k) B_{2k}.

This result was conjected by Helmert[16] and proved by Kawase.[17]

The factor (1-2k)(1+2k) results in poorer convergence of the series in terms \varphi compared to the one in \beta.

The quarter meridian is given by


m_p = \frac{\pi(a+b)}4 c_0 = \frac{\pi(a+b)}4 \sum_{j=0}^\infty
\biggl(\frac{(2j-3)!!}{(2j)!!}\biggr)^2 n^{2j},

a result which was first obtained by Ivory.[18]

Numerical expressions[edit]

The trigonometric series given above can be conveniently evaluated using Clenshaw summation. This method avoids the calculation of most of the trigonometric functions and allows the series to be summed rapidly and accurately. The technique can also be used to evaluate the difference m(\varphi_1) - m(\varphi_2) while maintaining high relative accuracy.

Substituting the values for the semi-major axis and eccentricity of the WGS84 ellipsoid gives


\begin{align}
m(\varphi)&=\bigl(\frac{111\,132.95255}{1^\circ}\varphi
         -16\,038.509\sin 2\varphi
         +16.833\sin4\varphi
         -0.022\sin6\varphi
         +0.00003\sin8\varphi\bigr) \,\mathrm{metres} \\
&= \bigl(\frac{111\,132.95255}{1^\circ}\beta
         -5\,346.170\sin 2\beta
         -1.122\sin4\beta
         -0.001\sin6\beta
         -0.5\times10^{-6}\sin8\beta\bigr) \,\mathrm{metres}
\end{align}

For the WGS84 ellipsoid the quarter meridian is


 m_p= \frac{\pi(a+b)}4 c_0 =10\,001\,965.729 \,\mathrm{metres}.

The perimeter of a meridian ellipse is 4m_p=2\pi(a+b) c_0. Therefore \frac12(a+b)c_0 is the radius of the circle whose circumference is the same as the perimeter of a meridian ellipse. This defines the mean Earth radius as 6\,367\,449.146\,\mathrm{m}.

On the ellipsoid the exact distance between parallels at \varphi_1 and \varphi_2 is m(\varphi_1)-m(\varphi_2). For WGS84 an approximate expression for the distance \Delta m between the two parallels at one half of a degree from the circle at latitude \varphi is given by


  \Delta m=(111\,133 - 560\cos 2\varphi) \,\mathrm{metres}.

The inverse meridian problem for the ellipsoid[edit]

In some problems, we need to be able to solve the inverse problem: given m, determine \varphi. This may be solved by Newton's method, iterating


\varphi_{i+1} = \varphi_i - \frac{m(\varphi_i) - m}{M(\varphi_i)},

until convergence. A suitable starting guess is given by \varphi_0 = \mu where


\mu = \frac{\pi}2 \frac m{m_p}

is the rectifying latitude. Note that it there is no need to differentiate the series for m(\varphi), since the formula for the meridian radius of curvature M(\varphi) can be used instead.

Alternatively, Helmert's series for the meridian distance can be reverted to give[19][20]


\varphi = \mu + H'_2\sin2\mu + H'_4\sin4\mu + H'_6\sin6\mu + H'_8\sin8\mu + \cdots

where


\begin{align}
H'_2 &= \textstyle \frac{3}{2} n
         - \frac{27}{32} n^3 + \cdots,&
H'_6 &= \textstyle \frac{151}{96} n^3 + \cdots, \\
H'_4 &= \textstyle \frac{21}{16} n^2
         - \frac{55}{32} n^4 + \cdots,\qquad&
H'_8 &= \textstyle \frac{1097}{512} n^4 + \cdots.
\end{align}

Similarly, Bessel's series for m in terms of \beta can be reverted to give[21]


\beta = \mu + B'_2\sin2\mu + B'_4\sin4\mu + B'_6\sin6\mu + B'_8\sin8\mu + \cdots,

where


\begin{align}
B'_2 &= \textstyle \frac{1}{2} n
         - \frac{9}{32} n^3 + \cdots,&
B'_6 &= \textstyle \frac{29}{96} n^3 - \cdots, \\
B'_4 &= \textstyle \frac{5}{16} n^2
         - \frac{37}{96} n^4 + \cdots,\qquad&
B'_8 &= \textstyle \frac{539}{1536} n^4 - \cdots.
\end{align}

Legendre[22] showed that the distance along a geodesic on an spheroid is the same as the distance along the perimeter of an ellipse. For this reason, the expression for m in terms of \beta and its inverse given above play a key role in the solution of the geodesic problem with m replaced by s the distance along the geodesic and \beta replaced by \sigma the arc length on the auxiliary sphere.[15][23] The requisite series extended to sixth order are given by Karney,[24] Eqs. (17) & (21), with \epsilon playing the role of n and \tau playing the role of \mu.

See also[edit]

References[edit]

  1. ^ Isaac Newton: Principia, Book III, Proposition XIX, Problem III, translated into English by Andrew Motte. A (searchable) modern translation is available at 17centurymaths. Search the following pdf file for 'spheroid'.
  2. ^ a b Clarke, Alexander Ross, 1880: Geodesy. Clarendon Press. Recently republished at Forgotten Books
  3. ^ Osborne, P (2013)The Mercator Projections (Chapter 5)
  4. ^ F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark, editors, 2010, NIST Handbook of Mathematical Functions (Cambridge University Press).
  5. ^ Mathematica guide: Elliptic Integrals
  6. ^ Maxima, 2009, A computer algebra system, version 5.20.1.
  7. ^ Euler, L. (1755). "Élémens de la trigonométrie sphéroïdique tirés de la méthode des plus grands et plus petits" [Elements of spheroidal trigonometry taken from the method of maxima and minima]. Mémoires de l'Académie Royale des Sciences de Berlin 1753 (in French) 9: 258–293. Figures. 
  8. ^ a b Delambre, J. B. J. (1799): Méthodes Analytiques pour la Détermination d'un Arc du Méridien; précédées d'un mémoire sur le même sujet par A. M. Legendre, De L'Imprimerie de Crapelet, Paris, 72–73
  9. ^ Rapp, R, 1991: Geometric Geodesy, Part I, §3.6, pp. 36–40.
  10. ^ Bessel, F. W. (1837). "Bestimmung der Axen des elliptischen Rotationssphäroids, welches den vorhandenen Messungen von Meridianbögen der Erde am meisten entspricht" [Estimation of the axes of the ellipsoid through measurements of the meridian arc]. Astronomische Nachrichten (in German) 14 (333): 333–346. Bibcode:1837AN.....14..333B. doi:10.1002/asna.18370142301. 
  11. ^ Helmert, F. R. (1880): Die mathematischen und physikalischen Theorieen der höheren Geodäsie, Einleitung und 1 Teil, Druck und Verlag von B. G. Teubner, Leipzig, § 1.7, pp. 44–48. English translation available at http://geographiclib.sourceforge.net/geodesic-papers/helmert80-en.html
  12. ^ Krüger, L. (1912): Konforme Abbildung des Erdellipsoids in der Ebene. Royal Prussian Geodetic Institute, New Series 52, page 12
  13. ^ J. W. Hager, J.F. Behensky, and B.W. Drew, 1989. Defense Mapping Agency Technical Report TM 8358.2. The universal grids: Universal Transverse Mercator (UTM) and Universal Polar Stereographic (UPS)
  14. ^ A guide to coordinate systems in Great Britain, Ordnance Survey of Great Britain.
  15. ^ a b c Bessel, F. W. (2010). "The calculation of longitude and latitude from geodesic measurements (1825)". Astron. Nachr. 331 (8): 852–861. arXiv:0908.1824. doi:10.1002/asna.201011352. English translation of Astron. Nachr. 4, 241–254 (1825), §5. 
  16. ^ Helmert (1880), op. cit., §1.11
  17. ^ Kawase, K. (2011): A General Formula for Calculating Meridian Arc Length and its Application to Coordinate Conversion in the Gauss-Krüger Projection, Bulletin of the Geospatial Information Authority of Japan, 59, 1–13
  18. ^ Ivory, J. (1798). "A New Series for the Rectification of the Ellipsis". Transactions of the Royal Society of Edinburgh 4 (2): 177–190. doi:10.1017/s0080456800030817. 
  19. ^ Helmert (1880), op. cit., §1.10
  20. ^ Adams, Oscar S (1921). Latitude Developments Connected With Geodesy and Cartography, (with tables, including a table for Lambert equal area meridional projection). Special Publication No. 67 of the US Coast and Geodetic Survey. A facsimile of this publication is available from the US National Oceanic and Atmospheric Administratio (NOAA) at http://docs.lib.noaa.gov/rescue/cgs_specpubs/QB275U35no671921.pdf, p. 127
  21. ^ Helmert (1880), op. cit., §5.6
  22. ^ Legendre, A. M. (1811). Exercices de Calcul Intégral sur Divers Ordres de Transcendantes et sur les Quadratures [Exercises in Integral Calculus] (in French). Paris: Courcier. p. 180. OCLC 312469983. 
  23. ^ Helmert (1880), op. cit., Chap. 5
  24. ^ Karney, C. F. F. (2013). "Algorithms for geodesics". Journal of Geodesy 87 (1): 43–42. arXiv:1109.4448. Bibcode:2013JGeod..87...43K. doi:10.1007/s00190-012-0578-z (open access). Addenda. 

External links[edit]

Category:Geodesy Category:Lines of longitude