How long can you trust a Starlink TLE? An empirical comparison of SGP4 and high-fidelity propagation against operator-updated truth across a megaconstellation

2026 年 4 月の 24,000 組を超えるスターリンクの TLE 対に関するこの実証的研究は、位置誤差が時間とともにべき乗則に従って増大することを、公的 TLE を用いた高忠実度伝播が運用者エポック残差とモデル整合バイアスのため標準的な SGP4 モデルを上回る性能を一般的に発揮しないことを、そして SGP4 の陳腐化が特定の高度において太陽フラックスと統計的に有意だが較正されていない相関を示すことを明らかにしている。

要約

特定の自動車が、1 時間後、3 時間後、あるいは 1 週間後に高速道路のどこにいるかを予測しようとしていると想像してください。それを助けるための 2 つのツールがあります。

1. 「公式マップ」（SGP4）： これは政府が無料で提供する、シンプルで高速かつ広く利用されている方法です。平均的な交通パターンに基づいて良い推定値を提供する、標準的な GPS アプリのようなものです。
2. 「スーパーコンピュータシミュレーション」（高忠実度）： これは、空気抵抗、自動車の正確な形状、乗客の重量、さらには月の引力に至るまで、あらゆる微小な詳細を考慮する、複雑で物理学的なシミュレーションです。レーシングチームの風洞シミュレーションのようなものです。

この論文が問いかけるのはシンプルな質問です：もし同じ「公式マップ」データから始めるとしたら、「スーパーコンピュータシミュレーション」は実際に自動車がどこにいるかのより良い予測をもたらすのでしょうか？

研究者たちは、低軌道に存在する巨大な車隊のようなスペースXのスターリンク衛星を数千基研究し、その答えを探りました。ここでは、彼らが発見したことをシンプルな比喩を用いて示します。

1. 「鮮度」のルール（データをどのくらい信頼できるか？）

論文によると、「公式マップ」（SGP4）は驚くほど優れていますが、それはごく短い間だけのことです。

• 比喩： 衛星の位置データを天気予報だと考えてみてください。発表から 4 時間後に予報を確認すれば、通常は正確です。しかし、同じ予報を使って 7 日後の天気を予測しようとすれば、それは無用なものになります。
• 発見： スターリンクの場合、「公式マップ」は約4〜6 時間信頼できます。それ以降、誤差は増大し始めます。7 日目には、衛星はマップが示す場所から数十キロメートルも離れた場所にいる可能性があります。研究者たちは、この誤差が（べき乗則のような）予測可能なパターンで増大することを見つけました。つまり、最終更新からどれだけの時間が経過したかによって、データの「古さ」を数学的に推定できるということです。

2. 「スーパーコンピュータ」の驚き（詳細は役立つか？）

「スーパーコンピュータシミュレーション」（高忠実度）は、より多くの物理法則を知っているため、常に勝つだろうと思うかもしれません。しかし、そうではありませんでした。

• 比喩： 10 分後にランナーがどこにいるかを推測しようとしていると想像してください。
  • ツールA（SGP4）： 「時速 10 マイルで走っている」という単純な規則を使います。
  • ツールB（スーパーコンピュータ）： 風、靴の摩擦、筋肉の疲労を考慮する複雑なモデルを使いますが、ぼやけた写真に基づいてランナーの開始速度を推測しなければなりません。
  • 結果： 開始時の写真（公開データ）がぼやけていたため、複雑なモデルは間違った速度からスタートしました。単純な規則（SGP4）の方が実際にはうまく機能しました。なぜなら、それは同じぼやけた写真に「較正」されていたからです。複雑なモデルは、間違った出発点で賢すぎようとしたため、かえって軌道から外れてしまいました。
• 発見： ほとんどの衛星とほとんどの時間枠において、単純な「公式マップ」（SGP4）の方が、複雑なシミュレーションよりも正確でした。複雑なシミュレーションが勝ったのは、ある特定のケースだけでした。それは、最も新しく大型の衛星（v2-mini）について、長い時間（3〜7 日）経過した場合です。その特定のシナリオでは、単純なマップがあまりにもひどく失敗していたため、わずかに欠陥のある複雑なモデルであっても、より良い結果を出せたのです。

3. 「方向」の問題（誤差はどこで発生するか？）

論文は、衛星がどこで誤っていたかを調べました。

• 比喩： 衛星を線路を走る列車だと考えてみてください。誤差が発生するのは、ほとんどが列車が線路から外れる（横方向）か、空へ飛び出す（上下方向）ことによるものではありません。誤差が発生するのは、ほぼ完全に、列車が線路上で早すぎるか遅すぎるという点です。
• 発見： 衛星はほぼ常に正しい「車線」にいましたが、タイミングが数秒から数分ずれていました。これは、最大の誤差源が大気抵抗（空気抵抗）であり、それが衛星を軌道上で減速させたり加速させたりするためです。

4. 「太陽の気象」との関連

研究者たちは、太陽活動（黒点や太陽フレア）が予測を悪化させるかどうかを確認しようとしました。

• 比喩： 大気をスポンジだと考えてみてください。太陽が活発になると、それがスポンジを加熱し、膨張させて「厚く」（密度が高く）します。これにより、衛星はより多くの空気抵抗を感じることになります。
• 発見： 太陽がより活発なとき、予測がわずかに悪化するという手がかりが見つかりましたが、100% の確信を持って証明するにはデータが十分ではありませんでした。それは、雲の中にパターンを見て、嵐が来ることを確認するにはまだ雨が少ないようなものです。

一般の人々への結論

• 短期的には単純なマップを信頼する： 今後数時間以内にスターリンク衛星がどこにいるかを知る必要がある場合、無料の単純なデータ（SGP4）で十分です。
• 過度に複雑にしない： 完璧な出発点（一般市民は持っていません）がない限り、超複雑な物理モデルを使用しても役立ちません。実際には、出発データに含まれる小さな誤差を増幅させるため、むしろ状況を悪化させることが多いです。
• 「新しい」衛星に注意する： 最も新しく大型の衛星は、長期間にわたって単純なマップで追跡するのが困難です。それらの特定の衛星については、複雑なモデルが最終的に優れている可能性がありますが、それは数日間待ってからです。

要約： この論文は、公開データにおいては「少ない方が多い」ことを証明しています。出発情報が完璧でない場合、シンプルで適切に調整されたモデルは、複雑なモデルよりも優れていることが多いのです。最良の戦略は、未来をあまり遠くまで予測しようとするのではなく、データを頻繁に（数時間ごとに）更新することです。

技術概要

技術的概要：Starlink TLE に対する SGP4 と高忠実度軌道伝播の経験的比較

問題定義
低軌道（LEO）メガコンステレーション、特に SpaceX の Starlink の普及により、2 行要素（TLE）データの運用環境は変化している。TLE は学術および商業的な宇宙状況認識（SSA）における事実上の標準であるが、その精度は従来、更新頻度が日次でデータが希薄な古いデータセットに対して評価されてきた。Starlink コンステレーションは 1 日約 6 回 TLE を公開しており、「古さ（staleness）」の問題は週単位ではなく時間単位で生じるものとなっている。公開 TLE がこれらの短期間においてどのように伝播するか、また、軌道決定（OD）の残差を本質的に含む公開 TLE から初期化された場合、標準的な SGP4 解析モデルに対して高忠実度数値伝播器（例：GMAT）が精度向上をもたらすかどうかを理解する上での重要なギャップが存在する。

手法
本研究は、2026 年 4 月の 1 ヶ月間の期間において、501 機の Starlink 衛星にわたる 24,641 組の TLE ページの大規模な経験的スウィープを実施した。データセットは、3 つの高度シェル（540km、550km、560km）と 3 つのプラットフォーム世代（v1.0、v1.5、v2-mini）によって層別化されている。

1. データ構築: 開始 TLE は、衛星・日ごとに 1 回サンプリングされる。「真値（Truth）」は、運用者が後に公開する TLE（「次 TLE」）として定義され、SGP4 を用いてその自身のエポックで評価される。この「真値としての次 TLE」手法は、運用者の内部 OD に対する自己整合性チェックとして機能し、ターゲットエポックにおける OD 残差を含んでいることを認識している。
2. 伝播: 各開始 TLE は、2 つの方法を用いてターゲットエポック（$t_j$）まで前方伝播される。
   • SGP4: 標準的な解析モデル。
   • 高忠実度（Hi-Fi）: NASA の General Mission Analysis Tool（GMAT）を、EGM2008（70×70）重力、NRLMSISE-00 大気抵抗、円錐影太陽放射圧（SRP）、および太陽と月の第 3 体重力を構成して設定したもの。
3. フィルタリング: 半長径のジャンプが 100m を超える区間を含むペアを除外するようにフィルタリングし、比較に対して「無機動」の仮定が成立することを確認する。
4. 分析: 位置誤差（$\|\Delta r\|$）を代理真値に対して計算する。本研究は、誤差の増大を特徴づけるためにべき乗則モデル $\|\Delta r(\Delta t)\| \approx A \Delta t^k$ を適合させる。また、ペアごとの比較を行い、Hi-Fi が SGP4 を上回るケースの割合を決定し、誤差係数 $A$ を太陽フラックス（$F_{10.7}$）に対して回帰分析する。

主要な知見

1. 誤差増大のダイナミクス（仮説 H1）:

   • 位置誤差は、すべてのコホートにおいてべき乗則に従う。
   • SGP4 v1.x: 540km および 560km では亜線形成長（$k \lesssim 1$）を示し、550km でも亜線形である。これは、平均運動バイアスと定数加速度の混合を示唆している。
   • Hi-Fi v1.x: 一般的に超線形成長（$k > 1$）を示し、モデル化されていない定数軌道内加速度（おそらく抵抗に関連するもの）が誤差予算を支配していることを示している。
   • v2-mini コホート: 両方の伝播器に対して超線形成長（$k \approx 1.3\text{--}1.4$）を示すが、Hi-Fi 側は長期（$\Delta t \ge 3$ 日）において SGP4 よりも緩やかに成長する。
   • 規模: 6 時間の地平線において、中央値誤差は約 1km（OD 残差のフロアに制限される）。7 日になると、SGP4 の誤差は約 38km に達するのに対し、Hi-Fi の誤差はプールされた v1.x コホートで約 76km に達する。

2. 伝播器の性能（仮説 H2）:

   • 一般的な結果: 公開 TLE から初期化された高忠実度伝播器は、サンプリングされたすべての地平線（6 時間、1 日、3 日、7 日）において、SGP4 に対する中央値位置誤差を改善しない。SGP4 は約 65–75% のペアで勝利する。
   • メカニズム: この否定的な結果は、3 つの要因に起因する。
     1. OD 残差の支配: 短期間では、誤差は運用者の初期状態の不確実性に支配され、これが両方の伝播器に影響を与える。
     2. カーネルの整合性: 「真値」は SGP4 を用いて構築される。SGP4 の予測は、この真値との構造的な整合性（例：共有された WGS-72 定数、Brouwer 平均要素）から恩恵を受けるのに対し、Hi-Fi 伝播器は動的参照の不一致（TEME から MJ2000Eq への変換および疑似接線軌道初期状態）によって不利益を被る。
     3. 宇宙機特性のバイアス: Hi-Fi 伝播器は、バリスティック係数（$C_D A$）の外部推定値に依存する。これらの公開推定値の誤差は、積分された抵抗誤差として蓄積するのに対し、SGP4 の $B^*$ は TLE を生成するために使用された追跡データに直接適合される。
   • 例外: 長期（$\Delta t \ge 3$ 日）における v2-mini コホートのみが、Hi-Fi が多数のペアで SGP4 を上回る唯一の領域である（7 日で約 55–73% のペアで勝利）。これは、新しい v2-mini 衛星に対する運用者の適合済み $B^*$ が不安定（超線形誤差増大）であるように見えるためであり、Hi-Fi の明示的な抵抗モデリングが SGP4 の平均運動近似を上回る結果となる。

3. 太陽変調（仮説 H3）:

   • SGP4 の古さ係数と $F_{10.7}$ に対する探索的回帰は、560km で従来の有意水準をクリアする正の傾きを示し、熱圏密度勾配と一致する。しかし、中等度の太陽活動の 30 日間の窓は統計的検出力を制限しており、結果は較正された測定というよりは、方向の一貫性として解釈される。

意義と実用的含意
本論文は、公開 Starlink TLE に対する「古さアトラス」を提供し、実務家に高度と世代ごとに層別化された具体的な誤差予算（$A, k$）を提供する。その主な貢献は、直感に反する運用上の発見である。すなわち、ユーザーが改善された初期状態（内部 OD や第三者追跡からのものなど）を保有していない限り、または特に v2-mini 世代の長期的軌道を予測する場合を除き、公開 TLE 入力からの高忠実度伝播への投資は一般的に有益ではないということである。

本研究は、大多数のユースケースにおいて、SGP4 の「真値構築カーネル整合性」と Hi-Fi モデルにおける宇宙機特性の不確実性の蓄積が、高忠実度力モデルの理論的利点を相殺することを確立している。提供されたセルごとのべき乗則パラメータは、将来の高度化伝播器の研究（SGP4-XP や微分可能 SGP4 など）のベンチマークとして、および結合分析における不確実性伝播研究の基準として機能する。

限界
本分析は、単一のコンステレーション（Starlink）、中等度の太陽活動の 30 日間の窓、および「真値としての次 TLE」手法に制限されており、絶対誤差の解釈は運用者の OD 残差によって境界付けられる。v2-mini の抵抗面積は v1.x データからスケーリングされ、適合されたものではないため、モデル化の仮定を導入しているが、感度テストはこのことが主要な結論を変更しないことを示唆している。