TLE(2行軌道要素)の精度と限界 — 数万物体カタログを支える簡易軌道表現

この回で学ぶこと

前回までの航法データメッセージの回では、CCSDSのODMファミリー(OPM・OEM・OMM)という、機関間で軌道情報を厳密に交換するための現代的な標準を学びました。その中でOMMを説明するとき、「いわゆるTLEが運んできたのと同種の情報をCCSDSの標準構造で表現しなおしたもの」という言い方をしました。この回では、その**TLE(Two-Line Element set、2行軌道要素)**そのものを正面から扱います。

TLEは、1970年代からNORAD(北米航空宇宙防衛司令部)系の組織が地球周回物体のカタログ管理に使い続けてきた、わずか2行・各69文字の固定長テキストフォーマットです。2026年現在、追跡・公開されているカタログ物体は稼働衛星・デブリを合わせて数万個にのぼり、その大部分の軌道情報は今なおTLE(またはその内容を包んだOMM)として配信されています。メガコンステレーションの回で見たように、低軌道の物体数が急増する時代において、TLEは衝突回避の第一次スクリーニングを支える事実上の共通言語です。

一方でTLEには、深宇宙探査機の精密軌道決定に慣れた目から見ると驚くほど大きな限界があります。精度はkmオーダーであり、時間が経つほど劣化します。しかもTLEの「軌道要素」は、ニュートン力学の教科書に出てくる瞬時(接触)軌道要素ではなく、SGP4という特定の簡易伝搬モデルとセットでのみ意味を持つ疑似平均要素です。この回では、フォーマットの構造、SGP4との関係、精度の定量的な限界という3つの柱を通して、「TLEをどう正しく使い、どこから先は使ってはいけないか」を理解します。

直感的導入: 数万物体を「1物体2行」で回すという設計思想

まず、なぜこんなに切り詰めたフォーマットが生まれ、生き残っているのかを考えましょう。

radiometric navigationの回や航法データメッセージの回で見てきた深宇宙探査機の軌道決定は、いわば少数精鋭の世界です。対象はたかだか数機、1機ごとに測距・ドップラー・DDORの観測を割り当て、精緻な力学モデル(高次重力場、太陽輻射圧、スラスタ力積…)で重み付き最小二乗推定を回し、状態ベクトルと共分散をOPM/OEMで丁寧に受け渡す。1物体あたりにかける計算・運用コストは莫大ですが、その見返りとしてm〜cmオーダーの精度が得られます。

地球周回物体のカタログ管理は正反対の大量生産の世界です。追跡対象は数万個。その大半は応答してくれない非協力物体(デブリ、運用終了衛星)で、観測はレーダーや光学望遠鏡による受動的な追跡が頼りです。1物体ごとに精密な数値積分モデルを走らせ、巨大なエフェメリスファイルを配布していたら、計算資源も通信帯域も(そしてフォーマットが生まれた1970年代なら、パンチカードの枚数も)まったく足りません。そこで採られたのが、

1. 軌道をごく少数の平均的なパラメータ(軌道要素+大気抵抗係数)に圧縮し、
2. 誰もが同じ軽量な解析的伝搬モデル(SGP4/SDP4)を使って任意時刻の位置速度を復元する

という分業です。TLEはこの前半、つまり「圧縮された軌道パラメータの入れ物」であり、1物体あたりたった138文字(2行×69文字)で1つの軌道を表現します。数万物体のカタログ全体でも数MBに収まる——この徹底した軽さこそが、半世紀にわたってTLEが使われ続けてきた理由です。ただしこの分業には代償があり、それが本回の後半のテーマである「精度の限界」と「SGP4以外で使ってはいけない」という制約です。

定式化その1: TLEフォーマットの構造とチェックサム

実物を見ましょう。以下は広く例題として使われる、国際宇宙ステーション(ISS)の2008年のTLEです。

1 25544U 98067A   08264.51782528 -.00002182  00000-0 -11606-4 0  2927
2 25544  51.6416 247.4627 0006703 130.5360 325.0288 15.72125391563537

各行はきっかり69文字で、フィールドの位置(桁)そのものが意味を持つ固定長ASCII形式です。区切り文字はなく、「何桁目から何桁目までが何のフィールドか」が仕様で定められています。主なフィールドを整理すると次のようになります。

1行目(識別・時刻・力学補助情報)

フィールド	例の値	意味
カタログ番号	25544	NORADカタログ番号(ISS)
分類	U	Unclassified(非機密)
国際識別符号	98067A	1998年67番目の打上げのA号機
エポック	08264.51782528	2008年の通算264.51782528日(小数日でUTC時刻を表す)
$\dot{n}/2$	-.00002182	平均運動の1次時間微分の半分 [rev/day²]
$B^*$	-11606-4	大気抵抗パラメータ(後述の指数表記)
チェックサム	7	行末1桁

2行目(軌道要素本体)

フィールド	例の値	意味
軌道傾斜角 $i$	51.6416	[deg]
昇交点赤経 $\Omega$	247.4627	[deg]
離心率 $e$	0006703	小数点省略: 0.0006703 と読む
近地点引数 $\omega$	130.5360	[deg]
平均近点角 $M$	325.0288	[deg]
平均運動 $n$	15.72125391	[rev/day]
通算周回数	56353	打上げからの周回数
チェックサム	7	行末1桁

バイト数を節約するための独特な省略記法に注意してください。離心率は先頭の「0.」が省略され、$B^*$ の -11606-4 は $-0.11606\times10^{-4}$ を意味します(仮数部の小数点と指数の「e」まで省略した表記)。また、軌道長半径 $a$ はフィールドとして存在せず、平均運動 $n$ から計算します。$n$ を [rev/day] から [rad/s] に直し、ケプラーの第三法則を逆に解いて

$$a = \left(\frac{\mu}{n^2}\right)^{1/3}, \qquad n\,[\mathrm{rad/s}] = \frac{2\pi \, n\,[\mathrm{rev/day}]}{86400}$$

例のISSの値 $n = 15.72125391$ rev/day を代入すると $n \approx 1.1433\times10^{-3}$ rad/s、地球重力定数 $\mu = 398600.4418\ \mathrm{km^3/s^2}$ を使って $a \approx 6731$ km、高度にして約353 kmが得られます。

チェックサム

各行の69文字目はチェックサムで、その行の1〜68文字目について

$$c = \left(\sum_{\text{数字}} d_k + (\text{マイナス記号の個数})\right) \bmod 10$$

つまり「数字はその値を加算、マイナス記号は1として加算、英字・空白・小数点・プラス記号は無視」して10で割った余りを取ります。これはパンチカードやテレタイプ回線での1文字化けを検出するための、時代を感じさせる素朴な誤り検出です。誤り訂正符号のレッスン群で学んだような強力な符号と比べれば検出能力はごく限られますが(2箇所の誤りが相殺すれば検出できません)、フォーマット破損の一次チェックとして今もパーサ実装で必ず検証されます。

定式化その2: SGP4/SDP4と「疑似平均要素」— TLE最大の落とし穴

TLEを理解するうえで最も重要なのは、TLEはそれ単体では軌道を定義しないという点です。TLEは、SGP4/SDP4という特定の解析的伝搬モデルへの入力として設計されており、両者は不可分のペアです。

SGP4の考え方

SGP4(Simplified General Perturbations 4)は、周期が225分未満の(おおむね高度約6000 km以下の)物体向けの解析的軌道伝搬理論で、より高い軌道や深宇宙寄りの物体には月・太陽の重力と共鳴項を追加したSDP4(Simplified Deep-space Perturbations 4)が使われます。数値積分を一切行わず、摂動論に基づく閉じた式の評価だけで任意時刻の位置・速度を返すのが特徴で、1物体あたりの計算コストは数値積分より桁違いに小さく、数万物体への一斉伝搬が現実的になります。

SGP4が解析的に取り込んでいる主な摂動は、地球の扁平による $J_2$ 項(とその高次の一部)、および指数型大気モデルによる抵抗です。$J_2$ 摂動の永年変化(secular variation)は、軌道要素の長期的なドリフトとして次のように書けます。

$$\dot{\Omega} = -\frac{3}{2}\, n J_2 \left(\frac{R_E}{p}\right)^{2} \cos i, \qquad \dot{\omega} = \frac{3}{4}\, n J_2 \left(\frac{R_E}{p}\right)^{2} \left(5\cos^2 i - 1\right)$$

ここで $J_2 \approx 1.08263\times10^{-3}$、$R_E = 6378.137$ km は地球赤道半径、$p = a(1-e^2)$ は半直弦です。先ほどのISSの例($a \approx 6731$ km、$i = 51.64°$)で計算すると $\dot{\Omega} \approx -5.1$ °/day となり、軌道面が1日に約5度も西向きに回転することが分かります。SGP4はこうした永年項に加えて、1周回の中で振動する短周期項・長周期項も解析式で評価します。大気抵抗はTLEの $B^*$ パラメータ

$$B^* = \frac{C_D A}{2m}\,\rho_0 R_E$$

($C_D$: 抵抗係数、$A/m$: 面積質量比、$\rho_0$: 基準大気密度)を通じてモデル化され、平均運動の増加(=高度低下)として現れます。

「SGP4前提の疑似平均要素」という本質

ここからが最大の注意点です。摂動論では、短周期振動をならした**平均要素(mean elements)と、その瞬間の位置速度に厳密に対応する接触要素(osculating elements、瞬時要素)を区別します。TLEに書かれている $i, \Omega, e, \omega, M, n$ は接触要素ではありません。それどころか、一般的な(たとえばBrouwer理論の)平均要素とも厳密には一致しません。TLEの要素は、「SGP4に入力したときに観測と最もよく合うように推定された、SGP4固有の内部表現としての疑似平均要素」**です。カタログ側の軌道決定も「SGP4で伝搬した軌道が観測に合うように要素を推定する」形で行われており、モデルの癖ごと要素の数値に焼き込まれています。

このことの実務上の帰結は明快です。

• TLEの要素をSGP4以外の伝搬モデル(たとえば自作の高精度数値積分器)に入力してはいけません。 TLEの $n$ から求めた $a$ を接触要素として数値積分の初期値に使うと、短周期項の分だけ初期状態が最初からずれます。低軌道ではこのずれだけで位置にして約10 kmオーダーの誤差になり、せっかくの高精度積分器がTLE本来の精度より悪い答えを出すという本末転倒が起きます。
• 逆に、SGP4で伝搬する限り、この内部表現のずれは自動的につじつまが合います。TLEの精度とは「TLE+SGP4」というペアの精度であり、そのペアで使うことが唯一の正しい使い方です。

なお座標系にも注意が必要で、SGP4の出力はTEME(True Equator, Mean Equinox)と呼ばれるやや特殊な座標系で定義されます。前回のOMMの例で REF_FRAME = TEME、MEAN_ELEMENT_THEORY = SGP4 というメタデータが明示されていたのは、まさにこの「モデルと座標系を書かなければ数値が解釈できない」というTLEの性質を、CCSDS標準の枠組みで補ったものです。

定式化その3: 精度の限界と誤差成長

では「TLE+SGP4」ペアの精度はどの程度でしょうか。物体・軌道・観測頻度に依存しますが、低軌道の典型値として広く引用されるのは次のようなオーダーです。

• エポック直後: 位置誤差 1 km 前後(良好な場合で数百m、悪い場合は数km)
• 誤差成長: 1日あたり 1〜3 km 程度ずつ劣化(大気活動が激しい時期の低高度物体ではさらに急速)
• 誤差の方向: 大部分が**進行方向(along-track)**成分

誤差が進行方向に集中する理由は、軌道力学の基本的な性質から定量的に説明できます。軌道長半径の推定誤差 $\Delta a$ は、ケプラーの第三法則を通じて平均運動の誤差

$$\Delta n = -\frac{3}{2}\frac{n}{a}\,\Delta a$$

を生み、これが時間に比例して位相(進行方向の位置)のずれとして蓄積します。

$$\Delta s(t) \approx a\,|\Delta n|\,t = \frac{3}{2}\,n\,\Delta a \, t$$

たとえば低軌道($n \approx 1.14\times10^{-3}$ rad/s)で $\Delta a = 20$ m というわずかな軌道長半径誤差があるだけで、1日後の進行方向誤差は $\Delta s \approx \frac{3}{2}\times1.14\times10^{-3}\times0.02\times86400 \approx 3$ km に達します。さらに大気密度の予測誤差(太陽活動によって密度は数十%単位で変動します)は $\dot{a}$ の誤差、すなわち $\Delta a$ 自体の時間成長を生むため、進行方向誤差は時間の2乗で効いてくる成分も持ちます。TLEが「エポックから離れるほど急速に信用できなくなる」のは、この機構によるものです。

深宇宙精密軌道決定との対比

この数字を、これまで学んできた深宇宙航法と並べてみましょう。radiometric navigationの回で見たように、深宇宙探査機の精密軌道決定では、測距がmオーダー、DDORが数nradの角度精度(1 AUで百km台の横方向精度、地球近傍なら遥かに良い)を達成し、再構成後の軌道はm〜cmオーダーで決まることも珍しくありません。GNSS衛星の精密暦に至っては数cmの世界です。

手法	典型精度	対象物体数	1物体あたりのコスト
TLE + SGP4	km オーダー(劣化あり)	数万	極小(解析式の評価のみ)
運用者の精密軌道暦(OEM)	数十m〜数百m	自社機のみ	中(数値積分+自社観測)
深宇宙精密OD / GNSS精密暦	cm〜m	数機〜数十機	莫大

つまりTLEと精密ODの間には3〜6桁の精度差があります。これは技術の優劣ではなく、設計目標の違いです。TLEは「数万物体の全体を、そこそこの精度で、誰でも軽量に扱えるようにする」ことに最適化されており、その目的においては今なお他に代わるものがない、という位置づけで理解するべきです。

実務での使われ方

Space-TrackとCelesTrakによる公開配信

現在、TLEカタログの一次配信元は米宇宙軍系のSpace-Track.org(第18宇宙防衛飛行隊系の組織が維持するカタログの公開窓口)であり、研究者T. S. Kelso氏が運営するCelesTrakが、長年にわたり整理・加工されたカタログデータの事実上の標準的な入手先として広く使われています。近年は両者とも、前回で学んだCCSDSのOMM形式(KVN/XML/JSON)での配信を拡充しており、「TLEの中身をOMMの器で配る」という移行が進んでいます。固定長69文字という器はカタログ番号の5桁制限(いわゆるAlpha-5対応で暫定延命)など限界が見えており、長期的にはOMMが後継の器になりますが、中身が「SGP4前提の疑似平均要素」であるという本質は器が変わっても変わらない点に注意してください。

コンジャンクション評価における役割分担

メガコンステレーションの回で見たとおり、低軌道の物体密度の増大により、衝突回避(コンジャンクション評価)は衛星運用の日常業務になりました。ここでTLEは第一次スクリーニングの道具として使われます。数万物体×数日先という全数総当たりに近い計算では、軽量なSGP4伝搬が唯一現実的な選択肢だからです。ただしTLEのkmオーダー誤差では、接近距離数百mの判断はできません。そのため実務は次の2段構えになっています。

1. 一次スクリーニング: カタログ全体をTLE(+SGP4)で伝搬し、しきい値(たとえば数km〜数十km)以内に接近する候補ペアを粗く抽出する。
2. 精密評価: 候補となった物体について、運用者が提供する**精密軌道暦(OEM、共分散付き)**と高精度な特別摂動法(数値積分)による軌道を突き合わせ、CDM(Conjunction Data Message)で衝突確率とともに通知し、回避マニューバの要否を判断する。

つまり「TLEでふるいにかけ、OEMで判断する」という役割分担です。Space Data Association(SDA)加盟事業者や大手コンステレーション運用者が自社機の精密エフェメリスを積極的に共有するのは、TLEだけでは最終判断に足る精度がないことの裏返しでもあります。

標準実装という互換性の要

「TLE+SGP4はペア」という性質は、SGP4の実装そのものがカタログの一部であることを意味します。実装ごとに式の打ち切りや定数が異なれば、同じTLEから異なる軌道が出てしまうからです。この問題に対しては、Valladoらが2006年に公開した論文 “Revisiting Spacetrack Report #3” と、それに付随する検証済み公開実装(いわゆるSGP4公式実装。C++/Fortran/Python等に移植され、python-sgp4 などのライブラリの基礎になっています)が事実上の参照実装として機能しており、世界中の運用ソフトウェアがこの実装との一致を検証しています。「フォーマット標準」だけでなく「アルゴリズムの参照実装」まで含めて初めて相互運用が成立するという、航法データメッセージの回で学んだ標準化の議論の、もう一歩踏み込んだ実例です。

演習問題

1. 本文のISSのTLE 2行目 2 25544 51.6416 247.4627 0006703 130.5360 325.0288 15.72125391563537 について、チェックサムの規則(数字は加算、マイナス記号は1、その他は無視、mod 10)を最後の1桁を除いた68文字に適用し、行末のチェックサム 7 が正しいことを確かめてください。また、この方式では検出できない誤りのパターンを1つ挙げてください。

2. 平均運動 $n = 15.72125391$ rev/day から軌道長半径 $a$ を計算し、高度(円軌道近似で $a - 6378$ km)を求めてください($\mu = 398600.4418\ \mathrm{km^3/s^2}$)。さらに、この $a$ は疑似平均要素由来の値であり、実際の瞬時の地心距離とは短周期項の分だけ異なることを踏まえ、この値を「その瞬間のISSの正確な高度」として使ってよいかどうか論じてください。

3. 軌道長半径の推定誤差 $\Delta a = 50$ m を持つ低軌道物体($n \approx 1.14\times10^{-3}$ rad/s)について、進行方向誤差の成長式 $\Delta s(t) \approx \frac{3}{2} n \Delta a\, t$ を用いて、3日後の進行方向誤差を km 単位で求めてください。この誤差は、接近距離しきい値 1 km のコンジャンクション精密判断に耐えるでしょうか。また、一次スクリーニングのしきい値を数十kmと大きく取る理由をこの結果から説明してください。

4. ある学生が「TLEから $a, e, i, \Omega, \omega, M$ を読み取り、それを初期接触要素として自作の高精度数値積分器($J_2$〜高次重力場・大気抵抗込み)に入力すれば、SGP4より正確な予測ができるはずだ」と主張しています。この主張のどこが誤っているかを「疑似平均要素」「短周期項」「TLE生成時の軌道決定がSGP4を前提としていること」の3つの観点から説明し、精度を上げたい場合の正しいアプローチ(どのようなデータを入手すべきか)を前回の内容も踏まえて述べてください。

まとめと次回予告

この回では、地球周回物体カタログの共通言語であるTLEを、フォーマット・伝搬モデル・精度の3つの面から整理しました。TLEは1物体138文字という徹底した軽さで数万物体のカタログ配信を可能にする一方、その軌道要素はSGP4前提の疑似平均要素であって瞬時要素ではなく、必ずSGP4/SDP4とペアで使わなければならないこと。精度はkmオーダーで、大気抵抗と軌道長半径誤差に起因する進行方向誤差が1日数kmのペースで成長すること。そして実務では、TLEによる軽量な一次スクリーニングと、運用者精密軌道暦(OEM)による精密評価という2段構えの役割分担でコンジャンクション評価が成立していることを見ました。深宇宙のm〜cmオーダーの精密軌道決定とは3〜6桁も違う世界ですが、それは「数万物体を全員で回す」ためにあえて選ばれた設計であり、器がOMMに移行しても本質は変わりません。

次回は視点を再び深宇宙に戻し、探査機の航法精度そのものを底上げする鍵として注目される深宇宙搭載可能な小型原子時計技術に触れます。地上局の水素メーザーに頼る往復(2-way)計測から、探査機自身が高安定な時刻を持つことで可能になる片道(1-way)計測への転換が、軌道決定と自律航法をどう変えうるかを見ていきます。

参考文献

• F. R. Hoots, R. L. Roehrich, Models for Propagation of NORAD Element Sets, Spacetrack Report No. 3, 1980
• D. A. Vallado, P. Crawford, R. Hujsak, T. S. Kelso, “Revisiting Spacetrack Report #3”, AIAA 2006-6753
• D. A. Vallado, Fundamentals of Astrodynamics and Applications, Microcosm Press
• CCSDS 502.0-B-3, Orbit Data Messages(OMM形式の規定)
• CelesTrak (T. S. Kelso), https://celestrak.org — TLE/OMMカタログとSGP4関連資料
• Space-Track.org — 米宇宙軍系カタログの公開配信サイト
• O. Montenbruck, E. Gill, Satellite Orbits: Models, Methods, and Applications, Springer