子午線弧
子午線弧（しごせんこ、Meridian arc）とは、測地学において、地球表面または地球楕円体に沿った子午線（経線）の弧を指す。子午線は楕円弧で南北方向に延びる測地線となる。
天文学において、2地点の天文緯度測定と子午線弧の長さとを結合することで地球の円周・半径を決定した。その始まりは、紀元前3世紀のエジプトのエラトステネスで、地球が球体であることを定量的に示した。
緯度差1分に相当する子午線弧長は、海里の定義にも参考にされた。
Contents
エラトステネスの子午線弧
アレクサンドリアの科学者エラトステネスによる測定は、地球の大円周長を計算した最初であった。彼は、夏至の正午において、太陽が古代エジプトの都市シエネ（現在のアスワン）で天頂を通過するということを知っていた。一方で、彼は自身の測定結果から、彼の居住地であるアレクサンドリアで、同時刻の太陽天頂距離が天球大円周長の1/50であるということも日時計が作る角度(7.2°)によって既知としており、天球と地球は同心であることから、アレクサンドリアがシエネの真北にあるならばアレクサンドリア-シエネ間の距離は地球の大円周長の1/50でなければならないと結論づけた。隊商の往来日数のデータを使って、彼はアレクサンドリア-シエネ間の距離を5,000スタディアであると推定した。
この結果は250,000スタディアの地球周長を意味し、単位スタディオンをアッティカスタディオン (185m) と仮定すると、これは46,250kmに相当し、現在の値から約16%大きい。しかし、エラトステネスがエジプトスタディオン (157.5m) を使ったとすれば、彼の測定値は 39,375km（わずか1%程度の誤差）であることが分かる。いずれにしても、幾何設定と古代の状況を斟酌すれば、16%の誤差は称賛に値するものである。
シエネは、正確にアレクサンドリアの真南にはなく、太陽の軌道は想定よりも0.5°傾いていた。また、ナイル川に沿って、または、砂漠を行旅することからの陸路の距離はおよそ10%程度の誤差があったとされる。
エラトステネスによる地球形状の見積もりは、その後何百年もの間受け入れられた。およそ150年後にポセイドニオスが同様の方法によりアレクサンドリア-ロドス島間の緯度差を測定するとともに、子午線弧長を船の速度と航海の期間から仮想的に割り出し、地球周長の算出を試みた。
中世から近世にかけての子午線弧の測量
8世紀に入ると中国でも子午線の計測が行われた。玄宗より新暦編纂の勅命を受けた僧一行は、鉄勒から交州にかけての測量を実施し、緯度1度の子午線弧長を351里80歩（約123.7km）と算出した。この算定と実際との誤差は11パーセントである。9世紀前期には、アッバース朝第7代カリフであるアル＝マアムーンの命により、アル＝フワーリズミーがシンジャール平原において実施した、角度測量によって多少良い結果が算出された。ヨーロッパでは、それまで子午線弧長測量が行われた記録が残っておらず、14世紀にジョン・マンデヴィルが編纂したとされる"The Travels of Sir John Mandeville" （大場正史訳「東方旅行記」, 東洋文庫第19巻, 平凡社, 1964, ISBN 9784582800197）において地球が球形であることが言及されている程度であったが、16世紀になって、もともと医師、生理学者であり、天文学、数学にも関心を持ったジャン・フェルネルが、経度がほぼ等しいパリ-アミアン間の緯度差を1度とみなした上で、荷車の車軸の回転数からその子午線弧長を決定したことを、著書"Ioannis Fernelii Ambianatis Cosmotheoria, libros duos complexa" (1528)に書き記している。
1615年には三角測量によるものとしては最初の子午線弧長測量がヴィレブロルト・スネルにより行われたが、測量結果には数パーセントの誤差があった。その約半世紀後の1669年にジャン・ピカールが本格的な三角測量を行い、緯度差1度に相当する子午線弧長を0.3%程度の精度で測定した。しかしながら、この頃辺りまでは地球の形状はあくまでも真球であるという前提の下に議論が行われていた。
フランス科学アカデミー遠征隊のペルーとラップランドへの派遣
ピカールによる測量以降、測量精度が向上するにつれて、地球の正確な形状についての問題が顕在化し、地球は正確には真球より回転楕円体と考えるべきとの意見が多くなったが、長球なのか扁球なのかについて議論が分かれていた。ジャック・カッシーニは、1713年に自らが行ったダンケルク-ペルピニャン間の測量結果を『地球の大きさと形状』（De la grandeur et de la figure de la terre、1720年）に取りまとめ、この結果とルネ・デカルトの渦動説から、地球が南北に長い長球であることを提唱した。一方では、振り子時計をパリから赤道付近へ持ってゆくと遅くなるというジャン・リシェによる報告からの推測により、アイザック・ニュートンが発表した万有引力の理論から赤道方向に長い扁球であると主張する学者も多数いた。
これを受け、18世紀半ば（1735年～1740年）には、フランス科学アカデミーが、地球楕円体の形状の論争に決着をつけるために赤道近傍と北極近傍の子午線弧長を比較した。この測量事業は、ピエール・ブーゲ、ルイ・ゴダン、シャルル＝マリー・ド・ラ・コンダミーヌ、ピエール・ルイ・モーペルテュイ及びアントニオ・デ・ウジョーアらによってペルー（現在のエクアドル）[1]とラップランド（トルネ谷）で実行された。
測量結果は2地域の同緯度差での子午線弧長に対する有意差を示し、極付近の弧長が赤道付近の弧長よりも大きいというものであった。これは赤道付近のほうが極付近よりも曲率が大きいことを示唆しており、1687年にニュートンが彼の著書『自然哲学の数学的諸原理』の第3巻において提唱したとおり、地球の数学的形状は扁球として解釈できることが確認された。カッシーニが得た測量結果が不正確であったことは、彼の弟子ともいうべきニコラ・ルイ・ド・ラカーユが1739年から2年を費やして再測量を行うことにより確認された。
18世紀後半にかけて、フランス科学アカデミーによってダンケルク-バルセロナ間の子午線弧長の測量が行われ、メートルの定義のために使われた。
伊能忠敬
日本では伊能忠敬が第二次測量（1801年）の結果から緯度1度に相当する子午線弧長を28.2里と導き出している。
子午線弧長の計算
子午線弧長の計算は地図投影法、特に横メルカトル図法（ガウス・クリューゲル図法）において重要な役割を果たす。
赤道から地理緯度 [math]\varphi\,[/math] までの子午線弧長 [math]S(\varphi)\,[/math] は、楕円積分が含まれているため、初等関数では表すことができないが、[math]\varphi\,[/math] の一次単項式と [math]\varphi\,[/math] の偶数倍を位相とする正弦高調波の無限級数の一般式で書き表すことができる。またこれを指定した次数で打ち切れば有限級数の形で近似計算に用いることができる。オイラーは1755年に第三離心率 [math]e^{\prime\prime}\,[/math] の二乗を微小量として用いて無限級数の一般式を得た。
第一離心率を用いた表式
地球楕円体の長半径を [math]a\,[/math]、第一離心率を [math]e\,[/math]、子午線曲率半径[2]を [math]M_\varphi = \frac{a(1-e^2)}{(1-e^2\sin^2\varphi)^{3/2}}\,[/math] とするとき、赤道から地理緯度 [math]\varphi\,[/math] までの子午線弧長 [math]S(\varphi)\,[/math] は以下のように与えられる。
- [math] S(\varphi)=\int_0^\varphi M_\theta\mathrm{d}\theta=a(1-e^2)\Pi(e^2;\varphi,e) [/math]
歴史的に広く用いられてきた [math]S(\varphi)\,[/math] の無限級数式は、ジャン＝バティスト・ジョゼフ・ドランブルが [math]e^2[/math] を微小量として用い1799年に公表したものである[3]。[math]e^8[/math] で打ち切った近似式は下記となる。
- [math] \begin{align}S(\varphi)\approx &\;a(1-e^2)\left\{\left(1+\frac{3}{4}e^2+\frac{45}{64}e^4+\frac{175}{256}e^6+\frac{11025}{16384}e^8\right)\varphi\right. \\ &\ -\frac{1}{2}\left(\frac{3}{4}e^2+\frac{15}{16}e^4+\frac{525}{512}e^6+\frac{2205}{2048}e^8\right)\sin 2\varphi \\ &\ +\frac{1}{4}\left(\frac{15}{64}e^4+\frac{105}{256}e^6+\frac{2205}{4096}e^8\right)\sin 4\varphi \\ &\ -\frac{1}{6}\left(\frac{35}{512}e^6+\frac{315}{2048}e^8\right)\sin 6\varphi \\ &\ +\frac{1}{8}\left.\left(\frac{315}{16384}e^8\right)\sin 8\varphi\right\}\\ \end{align} [/math]
しかしながら、これはヘルメルトの式などに比べると、共通係数として[math](1-e^2)[/math]を括り出していることが原因で[4]、[math]\left(\cdots\right)[/math] 内で [math]e^2[/math] の冪乗の級数の収束性が劣り、精密計算には多くの項数を必要とする。
第三扁平率を用いた表式
フリードリヒ・ヴィルヘルム・ベッセルは1825年に更成緯度[math]\beta=\tan^{-1}\left(\sqrt{1-e^2}\tan\varphi\right)[/math] で表した子午線弧長 [math]S(\beta)[/math] に対して第三扁平率 [math]n=\frac{1-\sqrt{1-e^2}}{1+\sqrt{1-e^2}}[/math] を微小量として用い、二項定理を利用しフーリエ級数展開を行った一般式を得た[5]。その級数係数は [math]n[/math] の偶数もしくは奇数冪乗の冪級数となる。
- [math] \begin{align} S(\beta) &= \frac{a}{1+n} \int_0^\beta \sqrt{1-2n\cos2\beta+n^2} \, d\beta \\ &=\frac{a}{1+n}\left( c_0 \beta + \sum_{l=1}^\infty (-1)^l \frac{c_l}{l} \sin 2 l \beta \right) . \\ c_l &= \sum_{k=0}^\infty a_k \, a_{k+l} .\\ a_k &= \binom{1/2}{k} n^k = (-1)^{k+1} \frac{\left(2k-3\right)!!}{\left(2k\right)!!} n^k . \end{align} [/math]
ここで、[math]j!![/math] は [math]j[/math] の二重階乗を表す。ただしこの式は子午線弧長の計算には広くは用いられなかった。なお一般式ではないがベッセルは、求長緯度 [math]\mu = \frac{\pi}{2}\,\frac{S}{S\!\left(\pi/2\right)}[/math] で [math]\beta[/math] を表す逆関数に当たる級数展開も示している。
楕円積分の関係式及び [math]n[/math] の符号反転を考えると、地理緯度 [math]\varphi[/math] で [math]S[/math] を表した一般式が得られる。これらの級数の収束性は他に知られている計算式よりも優れている。
- [math] \begin{align} S(\varphi) &=\frac{a}{1+n} \int_0^\varphi \frac{\left(1-n^2\right)^2}{\left(1+2n\cos2\varphi+n^2\right)^{\frac{3}{2}}} \, d\varphi \\ &=\frac{a}{1+n}\left( \int_0^\varphi \sqrt{1+2n\cos2\varphi+n^2} \, d\varphi - \frac{2n \sin2\varphi}{\sqrt{1 + 2n \cos2\varphi + n^2}} \right) \\ &=\frac{a}{1+n}\left( c_0 \varphi + \sum_{l=1}^\infty \frac{c_l}{l} \sin 2 l \varphi - \frac{2n \sin2\varphi}{\sqrt{1 + 2n \cos2\varphi + n^2}} \right) \end{align} [/math]
これらの無限級数は、含まれる [math]n[/math] の次数を [math]l_\max[/math] で打ち切れば有限級数となる。すなわち、[math]c_l[/math] を下記のように近似することになる。
- [math] c_{l,\,\textrm{approx}} = \begin{cases} \sum_{k=0}^{\lfloor (l_{\max}-l)/2\rfloor} a_k \, a_{k+l} & (l \le l_{\max}) \\ 0 & (l\gt l_{\max}) \end{cases} [/math]
ただし、[math]\lfloor x\rfloor[/math] は床関数（[math]x\,[/math] を超えない最大の整数）を表すものとする。
ヘルメルト・ベッセルの式
ベッセルはまた1837年に上記の [math]S(\varphi)[/math] の第一式に対しても二項定理の手法で級数展開を得た。括り出された共通係数は[math]a(1-n)^2(1+n)[/math]だった。
さらに、1880年にフリードリヒ・ロベルト・ヘルメルトが、括り出す共通係数を前節と同じである[math]\frac{a}{1+n}[/math]へ変更し、[math]n^4[/math] で打切られた近似式を提示した[6]。
- [math] \begin{align} S(\varphi)\approx &\;\frac{a}{1+n}\left\{\left(1+\frac{n^2}{4}+\frac{n^4}{64}\right)\varphi-\frac{3}{2}\left(n-\frac{n^3}{8}\right)\sin 2\varphi\right. \\ &\ \left.+\frac{15}{16}\left(n^2-\frac{n^4}{4}\right)\sin 4\varphi-\frac{35}{48}n^3\sin 6\varphi+\frac{315}{512}n^4\sin 8\varphi\right\}\\ \end{align} [/math]
これは一般式にするならば下記となる。
- [math] \begin{align} S(\varphi) &=\frac{a}{1+n}\left( c_0 \varphi + \sum_{l=1}^\infty \left(1 - 4 l^2\right) \frac{c_l}{l} \sin 2 l \varphi \right) . \end{align} [/math]
しかし前節の一般式と比べるならば [math]\frac{-2n \sin2\varphi}{\sqrt{1 + 2n \cos2\varphi + n^2}}[/math] の項[7]も級数展開したことは収束性を悪くしており、[math](1 - 4 l^2)[/math] の乗数に現れている。
またヘルメルトによる導出過程は一般論としては不備があり、一般式の導出・証明には至らないものだった。しかしヘルメルトの式は簡潔で精度も良いため近似式としては普及した。
河瀬の式
一般式としてのヘルメルトの式の証明自体については長い間放置されてきていたが、 最終的に2009年に河瀬和重により証明が行われた。
その際に用いられた一般式は、二項定理を経由するものではなく、ゲーゲンバウアー多項式による級数展開を利用し一種類の無限和に集約された形であった[8]。
- [math] \begin{align} S(\varphi)&=\frac{a}{1+n}\sum_{j=0}^\infty\left(\prod_{k=1}^j\varepsilon_k\right)^2\left\{\varphi+\sum_{l=1}^{2j}\left(\frac{1}{l}-4l\right)\sin 2l\varphi\prod_{m=1}^l\varepsilon_{j+(-1)^m\lfloor m/2\rfloor}^{(-1)^m}\right\} \\ \end{align} [/math]
ここで、[math]\varepsilon_i=3n/2i-n\,[/math] である。上式で [math]j=2\,[/math] まで取れば、ヘルメルトの提示した近似式が得られる[9][10]。級数を [math]j=J\,[/math] で打ち切れば、[math]n\,[/math] について [math]2J\,[/math] 次までで打ち切った近似式が得られることになる。
脚注
- ↑ 18世紀においては、エクアドルという国はまだ存在していなかった。当該地域は、当時スペインの管轄下に置かれており、後のキト市となる“キト特別行政区”と呼ばれていた。1830年に独立を果たした際に国の名称として採用された“エクアドル共和国”（「エクアドル」にはスペイン語で『赤道』の意味がある）には、“赤道付近の地域”として選ばれたこの地において実施されることとなった、フランス測地測量事業の名声が影響していると考えられている。
- ↑ 子午線曲率半径は平面曲線（楕円）の幾何学的性質から初等的に求められる。例えば、Rapp, R, (1991): Geometric Geodesy, Part I, §3.5.1, pp. 28–32参照。
- ↑ この式は日本でも広く用いられ、昭和61年版から平成21年版までの理科年表（地学部）にも掲載されていた。
- ↑ 共通係数 [math](1-e^2)[/math] を括り出さずに級数に組み込むか、もしくは[math]\left(1-\frac{1}{4}e^2\right)[/math] を括り出すなどで、収束性は改善される。
- ↑ 二項定理を利用した級数展開は、
- [math] \begin{align} & \left(1+2n\cos\theta+n^2\right)^\alpha \\ &= \left(1+n\, e^{i\theta}\right)^\alpha \left(1+n\, e^{-i\theta}\right)^\alpha \\ &= \left(\sum_{k=0}^\infty \binom{\alpha}{k} n^k e^{ik\theta} \right) \left(\sum_{k=0}^\infty \binom{\alpha}{k} n^k\, e^{-ik\theta} \right) \\ &= \sum_{k=0}^\infty \binom{\alpha}{k}^2 n^{2k} + 2 \sum_{l=1}^\infty \sum_{k=0}^\infty \binom{\alpha}{k} \binom{\alpha}{k+l}\, n^{2k+l} \cos l\theta \\ \end{align} [/math]
- ↑ ヘルメルトの提示では実際には式の形にまとまっていなかったが、1912年にヨハン・ハインリヒ・ルイ・クリューゲルがヘルメルトの結果を式の形に取りまとめている。
- ↑ この項は、不完全楕円積分項の [math]\varphi [/math] に関する二階微分に等しいので、級数展開形では乗数 [math]- 4 l^2[/math] が得られる。
- ↑ ゲーゲンバウアー多項式を利用した級数展開は、二項定理を利用した級数展開の和の取りまとめ方を変えることでも同様の結果が得られるが、
- [math] \begin{align} & \left(1+2n\cos\theta+n^2\right)^\alpha \\ &= \sum_{k=0}^\infty (-n)^k C_k^{(-\alpha)}(\cos\theta) \\ &= \sum_{k=0}^\infty (-n)^k \sum_{l=0}^k \binom{k-l-\alpha-1}{k-l} \binom{l-\alpha-1}{l} \cos (k-2l)\theta \\ &= \sum_{j=0}^\infty\left(\prod_{k=1}^j \nu_k^\alpha\right)^2\left(1+2\sum_{l=1}^{2j}\cos 2l\theta\prod_{m=1}^{l}\left(\nu_{j+(-1)^m\lfloor m/2\rfloor}^\alpha\right)^{(-1)^m}\right) \\ \end{align} [/math]
- ↑ 平成23年版の理科年表から、それまで掲載されていたドランブルの近似式に取って代わり、河瀬の一般式とヘルメルトの近似式が掲載されている。
- ↑ 同じ考え方に立てば、ベッセルが1825年に得た [math]S(\beta)[/math] 及び1837年に得た [math]S(\varphi)[/math] を次のように書き下すこともできる。
- [math] \begin{align} S(\beta)&=\frac{a}{1+n}\sum_{j=0}^\infty\left(\prod_{k=1}^j\bar{\varepsilon}_k\right)^2\left(\beta+\sum_{l=1}^{2j}\frac{\sin 2l\beta}{l}\prod_{m=1}^{l}\bar{\varepsilon}_{j+(-1)^m\lfloor m/2\rfloor}^{(-1)^m}\right) ,\\ S(\varphi)&=\frac{a(1-n^2)^2}{1+n}\sum_{j=0}^\infty\left(\prod_{k=1}^j \delta_k\right)^2\left(\varphi+\sum_{l=1}^{2j}\frac{\sin 2l\varphi}{l}\prod_{m=1}^{l}\delta_{j+(-1)^m\lfloor m/2\rfloor}^{(-1)^m}\right) . \end{align} [/math]
参考文献
- Euler, L. (1755). “Élémens de la trigonométrie sphéroïdique tirés de la méthode des plus grands et plus petits”. Mémoires de l'Académie Royale des Sciences de Berlin 1753 9: 258–293. Figures.
- 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
- Bessel, F. W. (1825): Ueber die Berechnung der geographischen Längen und Breiten aus geodätischen Vermessungen, Astronomische Nachrichten 4, 241–254
- Bessel, F. W. (1837): Bestimmung der Axen des elliptischen Rotationssphäroids, welches den vorhandenen Messungen von Meridianbögen der Erde am meisten entspricht, Astronomische Nachrichten, 14, 333–346
- 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, 44–48
- Krüger, L. (1912): Konforme Abbildung des Erdellipsoids in der Ebene, Veröffentlichung Königlich Preuszischen geodätischen Institutes, Neue Folge, 52, Druck und Verlag von B. G. Teubner, Potsdam, 12
- Florence, Trystram (2001). L'épopée du méridien terrestre (Le procès des étoiles). ISBN 978-2277220138.
- Florence, Trystram (1983/07). 地球を測った男たち. リブロポート. ISBN 978-4845700974.
- 河瀬和重 (2009): 緯度を与えて赤道からの子午線弧長を求める一般的な計算式, 国土地理院時報, 119, 45–55
- 飛田幹男, 河瀬和重, 政春尋志 (2009): 赤道からある緯度までの子午線長を計算する3つの計算式の比較, 測地学会誌, 55(3), 315–324