旧暦の計算《時刻系編》
本ページは旧暦の計算のサブページです。
本ページでは、時刻の様々な表現方法（時刻系）の差異について記載する。
（あくまで旧暦の計算に必要な範囲の記載に留め、厳密な説明は省略する）
目次
地方真太陽時
時刻の定義として最も素朴なものが地方真太陽時 LAT である。
これは、その地点で太陽が南中する時刻を「正午」とすることによって定まる時刻である。
太陽が南中する時刻の間隔（すなわち地方真太陽時における1日の長さ）は概ね 86,400 秒であるが、季節により変動する（下図）。
これは、太陽の進む速さが一定でないことと、太陽の軌道面（黄道面）に対して地球の地軸が約 23°26′ 傾いていることに因る。
地方真太陽時は、寛政暦・天保暦において作暦の基準として使われていた時刻系である。
地方平均太陽時
地方真太陽時は太陽の動きに即してはいるが、1日の長さが一定ではなく不便である。
そこで、地方真太陽時の1日の長さを均した時刻系である地方平均太陽時 LMT が広く用いられている。
太陽の南中時刻（地方真太陽時における正午 12:00）を地方平均太陽時で示したものが下図である。
この時刻と地方平均太陽時の正午 12:00 との差を均時差 ET(t) と呼ぶ。
同じ時刻 t を地方真太陽時で表した値 tLAT と地方平均太陽時で表した値 tLMT との関係は次式の通りとなる。
均時差 ET(t) は、太陽の平均黄経 λ̄☉(t)、太陽の真黄経 λ☉(t)、地球の赤道傾斜角 ε (= 23°26′) を用いて次式で表される。
平均黄経と真黄経については朔の計算編を参照。
但し、上式において arctan(tan λ☉(t)/cos ε) は λ☉(t) と同じ象限の角度（すなわち atan2(sin λ☉(t), cos ε cos λ☉(t))）を採るものとする。
なお、本稿において黄経の基準（黄経 0°）は冬至点の方向である。
旧暦の計算に使われる地方平均太陽時
一口に「地方平均太陽時」と言っても、観測地点の経度に依って時刻の値は異なる（時差）。
旧暦の計算に関係する地方平均太陽時には以下のものがある。
- 京都平均太陽時
- 明治改暦まで作暦の基準が置かれていた京都における地方平均太陽時。
本稿では、東経 135°45′45″（京都御苑内測候所）における時刻を採用し、GMT+9:03:03 に等しいものとして扱う。 - 東京平均太陽時
- 明治改暦以降、国際子午線会議（1884年）まで日本の子午線の基準が置かれていた東京における地方平均太陽時。
本稿では、東経 139°45′15″（江戸城天守台）における時刻を採用し、GMT+9:19:01 に等しいものとして扱う。 - 中央標準時 JST（明石平均太陽時）
- 東経 135° における地方平均太陽時。1888年より日本の標準時と定められている。
GMT+9:00:00 に等しい。 - 日本夏時刻 JDT
- 1948–1951年夏季に施行された夏時刻（サマータイム）。中央標準時を1時間進めたもの。
GMT+10:00:00 に等しい。 - グリニッジ標準時 GMT
- 経度 0°（本初子午線）における地方平均太陽時。世界時 UT1 及び UTC の元になった時刻。
本稿では、1971年以前は UT1 に、1972年以降は UTC に等しいものとして扱う。[注1]
この他にも日本で使用されていた時刻系はあるが（中央標準時制定以前の各地の地方平均太陽時や、1896–1937年に宮古列島以西で適用されていた西部標準時等）、旧暦の計算に関係するものではないので本稿では扱わない。
UT1 と UTC
平均太陽時は真太陽時を均した時刻系ではあるが、それでも1日の長さは一定とはならない。
これは、地球の自転周期が不規則に増減することに因る。
とは言え、時計の針の進みを地球の自転周期の不規則な変動に同期させるのは現実的ではない。
そこで、実用性を確保しつつ平均太陽時と歩調を合わせた時刻系として協定世界時 UTC が定められている。
これに対し、観測に基づく真のグリニッジ標準時を世界時 UT1 と呼ぶ。
今日我々が常用している時刻系は、この協定世界時を基準に定義されているものである。
（例えば、現在の日本時間 JST は UTC を9時間進めた時刻 UTC+9:00:00 に正確に等しい）
2022年以前の協定世界時
2022年まで、協定世界時 UTC は以下のように定められていた。
- UTC における1日の長さは基本的には丁度 86,400 秒（国際単位系秒）である。
- UT1 に対して UTC が 0.9 秒以上進む場合、余分な秒を1秒挿入して UTC の進みを遅らせる。（正の閏秒）
正の閏秒を適用（挿入）した日に限り、UTC における1日の長さは 86,401 秒となる。 - UT1 に対して UTC が 0.9 秒以上遅れる場合、時刻を1秒切り詰めて UTC の進みを早める。（負の閏秒）
負の閏秒を適用（削除）した日に限り、UTC における1日の長さは 86,399 秒となる。
協定世界時は、「秒」の刻みを一定（国際単位系秒）としつつ、時刻の値が世界時 UT1 に対して ±0.9 秒の範囲に収まるように定めた時刻系であった。
2022年以降の協定世界時
2022年、国際度量衡総会にて、閏秒による協定世界時 UTC の調整は当面行わないことが決まった。
すなわち、機械的に丁度 86,400 秒を以て1日とするということである。
今後どのように協定世界時 UTC と平均太陽時との同期を取っていくのか（或いは取っていかないのか）については、2026年の国際度量衡総会にて提案がある模様。
地球時
協定世界時は閏秒の挿入/削除により1日の長さが変動するので、地球の自転とは無関係な天体の運行の計算等には不向きである。
そのような計算には、以下で定まる地球時 TT という時刻系が用いられる。
- TT における1日の長さは常に丁度 86,400 秒（国際単位系秒）である。
- 1977年1月1日 00:00:00 (UTC) = 1977年1月1日 00:00:48.184 (TT) である。
ΔT = UT1 − TT
地球時 TT は地球の自転と同期していないため、世界時 UT1 との乖離は不規則に増減する。
この乖離を ΔT と書く。（下図）
同じ時刻 t を地球時 TT で表した値 tTT と世界時 UT1 で表した値 tUT1 との関係は次式の通りとなる。
ΔT' = UTC − TT
協定世界時 UTC と地球時 TT の時刻の乖離（1977年時点で 48.184 秒）は、閏秒を挿入する度に広がり、削除する度に狭まる。
本稿ではこの乖離を ΔT′ と呼ぶことにする。（下図）
同じ時刻 t を地球時 TT で表した値 tTT と世界時 UTC で表した値 tUTC との関係は次式の通りとなる。
1972年に現行の UTC が導入されてから2022年までは |ΔT − ΔT′| < 0.9 秒 が保たれていた。
2022年以降は当面 ΔT′ = 69.184 秒で固定されている。
ユリウス年数
天体の運行を計算する際は、次式で定められる J2000.0 と呼ばれる基準時刻からの経過時間の形で時刻を表すことが多い。
| J2000.0 | = | 2000年1月1日 12:00:00.000 (TT) |
| = | 2000年1月1日 11:58:55.816 (UTC) | |
| = | 2000年1月1日 20:58:55.816 (JST) |
この J2000.0 からの経過時間を 31,557,600 秒（365.25 日）を単位として表したものがユリウス年数 JY(t) である。
時刻 t が J2000.0 より前の時刻である場合は JY(t) は負の値となる。
| JY(t) | := | J2000.0 から t までの経過秒数/31,557,600 秒 |
| = | (tTT − 2000年1月1日 12:00)/31,557,600 秒 |
暦の上の日時は地球時 TT ではなく地方真太陽時 LAT や地方平均太陽時 LMT（現在は JST）であるから、旧暦を計算する際は実際には以下の形の式を用いることになる。
（以下の式において、時刻の差は閏秒等を考慮せず単純に日数と時刻の差により算出される時間を表す）
| JY(t) = | ((tLAT − 2000年1月1日 12:00) + ΔT − TZ − ET(t))/31,557,600 秒 | （寛政暦・天保暦） | |
| ((tLMT − 2000年1月1日 12:00) + ΔT − TZ)/31,557,600 秒 | （明治改暦から現行 UTC 導入まで） | ||
| ((tJST − 2000年1月1日 12:00) + ΔT′ − TZ)/31,557,600 秒 | （現行 UTC 導入以降） |
ここで、TZ は UTC（或いは GMT）に対する LAT 乃至 LMT の時差である。
| TZ = | +9:03:03（京都） | （1872年以前） | |
| +9:19:01（東京） | （1873–1887年） | ||
| +9:00:00（明石） | （1888年以降、但し1948–1951年夏季を除く） | ||
| +10:00:00（夏時刻） | （1948–1951年夏季） |
計算例1: 1989年7月6日3時（グレゴリオ暦・中央標準時）
例として、1989年7月6日 03:00（グレゴリオ暦・中央標準時）のユリウス年数を計算する。
- 2000年1月1日 12:00 から 1989年7月6日 03:00 までの時間は −3831 日 9:00:00
- 1989年7月当時 ΔT′ = 56.184 秒
- 協定世界時 UTC と中央標準時 JST との時差は TZ = +9:00:00
であるから、1989年7月6日 03:00（グレゴリオ暦・中央標準時）のユリウス年数は次式の通りとなる。
= ((1989年7月6日 03:00 − 2000年1月1日 12:00) + ΔT′ − TZ)/31,557,600 秒
= (−3831 日 9:00:00 + 56.184 秒 − 9:00:00)/31,557,600 秒
= (−3831 日 17:59:03.816)/31,557,600 秒
= (−331,063,143.816 秒)/31,557,600 秒
= −10.49075797
計算例2: 1854年2月13日12時（グレゴリオ暦・京都真太陽時）
別の例として、1854年2月13日 12:00（グレゴリオ暦・京都真太陽時）のユリウス年数を計算する。
- 2000年1月1日 12:00 から 1854年2月13日 12:00 までの時間は −53282 日
- NASA の略算式に拠ると、1854年2月当時 ΔT = 7.286 秒
- 協定世界時 UTC と京都時間との時差は TZ = +9:03:03
- 均時差は ET(1854年2月13日 12:00 (LAT)) = (52.8273° − arctan(tan 54.1398°/cos 23°26′))·1 日/360° = −868.836 秒
であるから、1854年2月13日 12:00（グレゴリオ暦・京都真太陽時）のユリウス年数は次式の通りとなる。
= ((1854年2月13日 12:00 − 2000年1月1日 12:00) + ΔT − TZ − ET(1854年2月13日 12:00 (LAT)))/31,557,600 秒
= (−53282 日 + 7.286 秒 − 9:03:03 − (−868.836 秒))/31,557,600 秒
= (−53282 日 8:48:26.878)/31,557,600 秒
= (−4,603,596,506.878 秒)/31,557,600 秒
= −145.87917037
参考資料
- 暦計算室（国立天文台）
- NASA GODDARD SPACE FLIGHT CENTER ECLIPSE WEB SITE（米国立航空宇宙局）
- The Astronomical Data Portal @ UK Hydrographic Office（英国王立航海暦局）
脚注
- ^ 1.
- 世界時 UT1 が導入されたのは1956年であるが、それ以前についても遡及的に GMT = UT1 であったものとして扱う。
- ^ 2.
- HMNAO の予測値は、LOD の予測値として掲載されている式 lod = +1.72t − 3.5 sin(2π(t + 0.75)/14) (t は1825年からの経過年数) を積分して算出した。
- ^ 3.
- 略算式を用いて太陽の平均黄経 λ̄☉(t) 及び真黄経 λ☉(t) を求めるにはユリウス年数 JY(t) が必要となるので、素直に考えると無限循環に陥ってしまう。とは言え、地方真太陽時と地方平均太陽時の差は然程大きくないことから、一旦地方真太陽時の時刻 tLAT をその侭地方平均太陽時の時刻と見做して（すなわち ET(t) を無視して）太陽の平均黄経及び真黄経を求めてしまえばそれでよい。こうして求まる ET(t) の値は真の ET(t) の値に対して最大0.2 秒程度の誤差を含むが、これは寛政暦・天保暦の精度に対して小さいので問題無い。