天文現象を予測するには,太陽・月・惑星などの運動について知らなければならない.時々刻々と変わる天体の位置や速度について,これまで本書ではいわゆる DE (Development Ephemeris) を採用してきた.DEは米国ジェット推進研究所 (Jet Propulsion Laboratory, JPL) が惑星探査用に編纂した月・惑星の暦であり,国際的にもっとも広く用いられているものである.具体的には昭和60年版 (1985) からはDE200,平成15年版 (2003) からはDE405を導入してきた.
DEは数値積分によって運動方程式を解くことで暦を推算しており,それに必要な初期位置・速度や運動方程式のパラメータである質量・天文定数などは観測値にもっともあうように推定されている.このため,国際天文学連合 (IAU) の勧告する天文定数系とは必ずしも一致しないが,DEの定数が後にIAUの天文定数系で採用されるなど,より精密な定数系といえるだろう.ところが,IAU 2009/2012 天文定数系では近年の月・惑星探査の成果を採り入れられ,この関係が逆転してしまっていた.
そこで,平成28年版からはDE405の後継にあたるDE430を採用することにした次第である.
DE405とDE430のおもな違いは以下のとおりである.
これまで天文単位 au は 1日 (1d = 86400s) あたり 0.01720209895 rad の角速度で太陽の周りを等速円運動する質点の軌道半径と定義されていた.この0.01720209895はガウス引力定数 k とよばれ,19世紀にガウスが定めたものである.周期 T = 2π/k ≈ 365.256898日 ≈ 1年であるから,これがほぼ地球の軌道半径に近いことは容易に理解できるだろう.
この定義を使うと,日心重力定数 GMS は GMS = k2 [ au3 d-2 ] と一定の値になるから,au の具体的な値がわからなくても天体の運動を調べることができる.光学望遠鏡による観測では天体との距離を測ることはできないから,こうするより手段がなかったわけだ.そして,さまざまなパラメータの1つとして au を推定,GMS = k2・au3 / 864002 により GMS を算出すればよい.しかし,測距を含む観測が主流となった現代では GMS もSI単位で直接求めることが可能であり,天文単位は必ずしも必要ではなくなってしまった.また,この方法では au や GMS の変動が何の変化なのかを判断しづらいし,相対論的な考慮も含まれていない.
そこで,2012年のIAU勧告では天文単位 au の再定義がなされ,時刻系によらず上記のような定義値を採用する一方,ガウス引力定数は定数系から外すことになったわけである.
天体別に具体的な黄経・黄緯・動径の違いをグラフにすると以下のようになる.ただし,惑星は日心座標,月は地心座標,いずれもDE405の数値からDE430の数値を引いたものである.暦象年表の表値は赤経で 0.001s (15mas‡),赤緯で 0.01 ″(10mas),距離で 0.0000001 au (15 km) までとなっているから,太陽・月・内惑星では時折四捨五入により違いが出る程度である一方,外惑星でははっきりと違いが出ることがわかる.理科年表は表示桁数が少ないので影響も小さく,暦要項ではほとんど影響は現れない.
最後に,おもなパラメータの違いをまとめておこう.
項目 | DE405 | DE430 | IAU2009 /2012 |
---|---|---|---|
au [km] | 149 597 870.691 | 149 597 870.700 | 149 597 870.700 |
GMS [×1020 m3s-2] | 1.327 124 400 18 | 1.327 124 400 419 | 1.327 124 400 41 |
質量比 (MS / MP) | |||
水星 | 6023600. | 6023682. | 6023600. |
金星 | 408 523.71 | 408 523.719 | 408 523.719 |
地球+月 | 328 900.5614 | 328 900.5598 | 328 900.5596 |
火星 | 3098 708. | 3098 703.59 | 3098 703.59 |
木星 | 1047.348 6 | 1047.348 625 | 1047.348 644 |
土星 | 3497.898 | 3497.9018 | 3497.9018 |
天王星 | 22902.98 | 22902.98 | 22902.98 |
海王星 | 19412.24 | 19412.26 | 19412.26 |
地球と月の質量比 | 0.012 300 038 | 0.012 300 0369 | 0.012 300 0371 |
‡) masとはmilliarcsecond,すなわち角度の1″の1/1000を表わす単位.→本文に戻る
暦象年表2016より