Impact of geomagnetic storm intensity on thermospheric density and orbital decay of LEO satellites: a G1–G5 per-event analysis

技術要約：地磁気嵐の強度による熱圏密度の変化と低軌道（LEO）衛星の軌道減衰への影響

問題提起
低軌道（LEO）衛星コンステレーションの急速な拡大に伴い、熱圏密度の正確な特性把握に対する運用上のニーズが高まっている。高度600 km以下の衛星にとって、大気抵抗は重力以外の主要な摂動要因であり、地磁気嵐は熱圏密度を著しく増大させ、それによって軌道減脱を加速させる。2022年2月のスターリンク衛星38機の喪失や、2024年10月のスターリンク衛星の予定外の再突入といった過去の事例は、深刻な運用リスクを実証している。既存の統計的手法は多くの場合、重畳エポック解析（SEA）に依存しているが、これは多くの事象の応答を平均化するため、即時の軌道予測や衝突回避に不可欠な鋭いピーク応答を平滑化してしまう可能性がある。さらに、経験的な熱圏モデルは、激しい嵐の発生時における密度の増大を正確に再現することに苦慮する場合がある。

手法
本研究では、「事象ごとのピーク（per-event peak）」解析手法を採用し、太陽周期23、24、25（2003年–2024年）にわたるG1からG5（最大$K_p$値 5から9）までの22件の地磁気嵐から、瞬時的なピーク応答を分離した。データセットには、2024年5月のガノン（Gannon）G5嵐が含まれている。

• データソース： 非摂動の中性質量密度および精密軌道データは、GRACE（2002–2017）およびGRACE-FO（2018–現在）の双子衛星の加速度計測定から導出された。地磁気指数（$K_p$, $Dst$, $SYM-H$, $AU$, $AL$）および太陽風パラメータは、NASA/GSFC OMNIWebから取得した。
• 事象の選定： 以下の3つの基準に基づき、22件の事象を選定した。(1) 最大$K_p \ge 5$（G1以上であること）、(2) 嵐の突発開始（SSC）の48時間前から嵐のピーク後48時間まで連続した加速度計データが存在すること、(3) 前に発生した嵐の主相と重複していないこと。選定には、G1–G4の各カテゴリーに5件ずつ、および利用可能な2件のG5事象（2003年10月および2024年5月）が含まれる。
• 処理プロセス：
  1. 軌道の選定： 各事象について、「嵐の軌道」は嵐のピーク（最小$Dst$）に最も近い軌道と定義し、「静穏参照軌道」は、少なくとも24時間の連続した$K_p \le 3$の後に発生したSSCに最も近い軌道とした。
  2. 減衰率の算出： 生の10秒間高度データを軌道ごとに平均化し、半日周変動を除去するために24時間の移動平均を用いて平滑化した。静穏期と嵐の期間に対して線形適合を行い、高度減衰率（$R_{quiet}$ および $R_{storm}$）を決定した。
  3. 増幅指標： 嵐の時の軌道減衰の増幅は、パーセンテージの増加として計算された：$\frac{R_{storm} - R_{quiet}}{R_{quiet}} \times 100$。
  4. 統計モデリング： 3つの利用可能な予測変数（静穏時の減衰率、衛星運用高度、最大$K_p$）を用いて、嵐の時の減衰率を予測する経験的多重線形回帰（MLR）モデルを開発した。

主な結果

• 強度との単調なスケーリング： 平均軌道減衰率は、嵐の強度とともに単調に増加する。観測された増幅率は、77.2% (G1)、98.5% (G2)、138.1% (G3)、149.6% (G4)、および166.5% (G5)であった。
• 決定的なモジュレーターとしての高度： 衛星運用高度は、独立した重要な要因として作用する。低高度の事象（例：約350 km）では、比較的弱い地磁気強制下においても強い抵抗応答を示し、高高度における平均的な増幅率をしばしば上回った。
• 事象特有の変動性： 強度カテゴリー内でも顕著なばらつきが存在する。特に、2つのG5事象（2003年10月と2024年5月）は、同一の$K_p$値と類似の高度であったにもかかわらず、減衰増幅において約2.6倍の差を示した。著者らは、2024年のガノン嵐における高い増幅を、背景となる太陽活動レベルが低いこと（静穏時の減衰率という分母が小さいこと）、および連続するCME（>30時間）によって駆動された持続時間の長さによるものと考えている。
• MLRモデルの性能： 実験的モデル $D_s = 1.55D_q - 0.29(h_0 - 450) + 7.39Kp_{max} - 35.60$ は、全22事象に対して決定係数（$R^2$）0.91、および平方根平均二乗誤差（RMSE）12.3 m/dayを達成した。このモデルは、衛星固有の物理パラメータ（質量、断面積、抗力係数）を必要とせずに、嵐の強度と高度への依存性を捉えることに成功した。

意義および主張
本論文は、既存のSEAベースの解析に対する補完的なベンチマークとして本研究を位置づけている。事象ごとのピーク応答に焦点を当てることで、平均化の手順によって不明瞭になる鋭い信号を保持しており、経験的な熱圏モデルに対するより直接的な検証リファレンスを提供している。本研究は、データセットとシンプルな、運用上アクセスしやすいMLRツールを提供しており、これにより衛星オペレーターは、標準的な追跡データ（高度、静穏時の減衰）と宇宙天気指数（$K_p$）のみを使用して、嵐の時の減衰率を推定できる。著者らは、これらの結果が、現在の太陽周期25の文脈において、拡大するLEO集団に対する改良された軌道予測および衝突回避戦略の開発をサポートすると主張している。

著者らは、MLRが2つの特定の衛星プラットフォーム（GRACE/GRACE-FO）から導出されており、それらが類似の物理的特性を持っていることを指摘し、モデルの汎用性に関して控えめな範囲を維持している。彼らは、このツールをより多様なLEO衛星集団に適用するためには、衛星の質量、断面積、および速度を組み込むための将来的な作業が必要であることを認めている。