子午線弧

子午線弧（しごせんこ、Meridian arc）とは、測地学において地球表面または地球楕円体に沿った子午線（経線）の弧を指す。子午線は楕円弧で南北方向に延びる測地線となる。
天文学において、2地点の天文緯度測定と子午線弧の長さとを結合することで地球の円周・半径を決定した。その始まりは、紀元前3世紀のエジプトのエラトステネスで、地球が球体であることを定量的に示した。
緯度差1分に相当する子午線弧長は、海里の定義にも参考にされた。
エラトステネスによる子午線弧長の推定.
アレクサンドリアの科学者エラトステネスによる測定は、地球の大円周長を計算した最初であった。彼は、夏至の正午において、太陽が古代エジプトの都市シエネ（現在のアスワン）で天頂を通過するということを知っていた。一方で、彼は自身の測定結果から、彼の居住地であるアレクサンドリアで、同時刻の太陽天頂距離が天球大円周長の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年）には、フランス科学アカデミーが、地球楕円体の形状の論争に決着をつけるために赤道近傍と北極近傍の子午線弧長を比較した。この測量事業は、ピエール・ブーゲ、ルイ・ゴダン、シャルル＝マリー・ド・ラ・コンダミーヌ、ピエール・ルイ・モーペルテュイ及びアントニオ・デ・ウジョーアらによってペルー（現在のエクアドル）とラップランド（トルネ谷）で実行された。
測量結果は2地域の同緯度差での子午線弧長に対する有意差を示し、極付近の弧長が赤道付近の弧長よりも大きいというものであった。これは赤道付近のほうが極付近よりも曲率が大きいことを示唆しており、1687年にニュートンが彼の著書『自然哲学の数学的諸原理』の第3巻において提唱したとおり、地球の数学的形状は扁球として解釈できることが確認された。カッシーニが得た測量結果が不正確であったことは、彼の弟子ともいうべきニコラ・ルイ・ド・ラカーユが1739年から2年を費やして再測量を行うことにより確認された。
18世紀後半にかけて、フランス科学アカデミーによってダンケルク-バルセロナ間の子午線弧長の測量が行われ、メートルの定義のために使われた。
伊能忠敬による子午線弧の測量.
日本では伊能忠敬が第二次測量（1801年）の結果から緯度1度に相当する子午線弧長を28.2里と導き出している。
子午線弧長の計算.
地球楕円体に基づく子午線弧長の計算は地図投影法、特に横メルカトル図法（ガウス・クリューゲル図法）において重要な役割を果たす。
赤道から地理緯度 formula_1 までの子午線弧長 formula_2 は、楕円積分が含まれているため、初等関数では表すことができないが、formula_1 の一次単項式と formula_1 の偶数倍を位相とする正弦高調波の無限級数の一般式で書き表すことができる。またこれを指定した次数で打ち切れば有限級数の形で近似計算に用いることができる。
第三離心率を用いた一般式.
オイラーは1755年に第三離心率 formula_5 の二乗を微小量として用いて無限級数の一般式を得た。
第一離心率を用いた表式.
地球楕円体の長半径を formula_6、第一離心率を formula_7として、子午線曲率半径は formula_8 となる。赤道から地理緯度 formula_1 までの子午線弧長 formula_2 は以下のようにformula_11の部分積分で与えられる。
歴史的に広く用いられてきた formula_2 の無限級数一般式は、ジャン＝バティスト・ジョゼフ・ドランブルが1799年に公表し、共通係数として率直に formula_14を括り出し、formula_15 を微小量として級数展開したものである。
しかしながら、これはヘルメルトの式などに比べると、係数 formula_17 の formula_18 内に formula_19の項が現れ、多くの項数を必要とする。また共通係数としてformula_20を括り出していることが原因で、formula_21 内で formula_15 の冪乗の級数の収束性が劣る。
第三扁平率を用いた表式.
更成緯度で表した表式.
フリードリヒ・ヴィルヘルム・ベッセルは1825年に更成緯度formula_23 で表した子午線弧長 formula_24 に対して、第三扁平率 formula_25を用い、共通係数としてformula_26 を括り出し微小量としてformula_27を用いて二項定理を利用しフーリエ級数展開を行った一般式を得た。その級数係数は formula_27 の偶数もしくは奇数冪乗の冪級数となる。
ここで、formula_30 は formula_31 の二重階乗を表す。ただしこの式は子午線弧長の計算には広くは用いられなかった。なお一般式ではないがベッセルは、求長緯度 formula_32 で formula_33 を表す逆関数に当たる級数展開も示している。
地理緯度で表した表式.
ここで楕円積分の関係式及び formula_27 の符号反転を考えると、地理緯度 formula_35 で formula_36 を表した一般式が得られる。これらの級数の収束性は他に知られている計算式よりも優れている（級数展開にformula_27の奇数次項が現れないなど）。
これらの無限級数は、含まれる formula_27 の次数を formula_40 で打ち切れば有限級数となる。すなわち、formula_41 を下記のように近似することになる。
ただし、formula_43 は床関数（formula_44 を超えない最大の整数）を表すものとする。
ヘルメルト・ベッセルの式.
ベッセルはまた1837年に上記の formula_45 に対しても同じく二項定理の手法で級数展開一般式を得た。括り出された共通係数はformula_46だった。
さらに、1880年にフリードリヒ・ロベルト・ヘルメルトが、括り出す共通係数を前節と同じformula_26へ変更し、formula_48 で打切られた近似式を提示した。
これは一般式にするならば下記となる。
しかしながら前節の一般式と比べるならば formula_51 の項も級数展開したことは収束性を悪くしており、formula_52 の乗数に現れている。
加えて、ヘルメルトによる導出過程は一般論としては不備があり、一般式の導出・証明には至らないものだった。しかしヘルメルトの式は簡潔で精度も良いため近似式としては普及した。
河瀬の式.
一般式としてのヘルメルトの式の証明自体については長年放置されていたが、
最終的に2009年に河瀬和重により証明が行われた。
その際に用いられた一般式は、二項定理を経由するものではなく、ゲーゲンバウアー多項式による級数展開を利用し一種類の無限和に集約された形であった。
ここで、formula_54 である。上式で formula_55 まで取れば、ヘルメルトの提示した近似式が得られる。級数を formula_56 で打ち切れば、formula_57 について formula_58 次までで打ち切った近似式が得られることになる。