J2摂動による軌道変動シミュレーター — 太陽同期軌道設計

地球は完全な球ではなく、赤道が約21km膨らんだ扁平回転体。この J2項（係数 1.08×10⁻³）が人工衛星の軌道面を歳差させ、近地点を回転させます。軌道長半径・離心率・傾斜角を変えると、昇交点赤経 Ω と近地点引数 ω の永年変動率がリアルタイムで分かり、太陽同期軌道（SSO）や凍結軌道の条件を探せます。

パラメータ設定

軌道長半径 a

地球中心からの平均距離。R_E + 高度

離心率 e

0=円軌道、0.5=長楕円軌道

軌道傾斜角 i

赤道面に対する軌道面の角度。90°超は逆行軌道

近地点引数 ω

昇交点から近地点までの角度

昇交点赤経 Ω₀

春分点を基準とした軌道面の方位

伝播日数

日

この期間で累積する Ω・ω の総変動を計算

計算結果

—

軌道周期 (min)

—

dΩ/dt (deg/day)

—

dω/dt (deg/day)

—

Ω総変動 (deg)

—

ω総変動 (deg)

—

太陽同期との差 (deg/day)

—

軌道と摂動の可視化

地球（中央・扁平を強調）の周りを衛星軌道（緑）が周回。軌道面はΩ方向に歳差（青矢印）、軌道楕円の長軸はω方向に回転（橙矢印）します。

dΩ/dt vs 軌道傾斜角 — 太陽同期条件

dω/dt vs 軌道傾斜角 — 凍結軌道条件

理論・主要公式

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

J2摂動による昇交点赤経の歳差率 dΩ/dt と近地点引数の回転率 dω/dt。n：平均運動 √(μ/a³)、p = a(1-e²)：半通径、R_E：地球赤道半径、i：軌道傾斜角。

$$n = \sqrt{\frac{\mu}{a^3}},\qquad T = 2\pi\sqrt{\frac{a^3}{\mu}}$$

平均運動 n と軌道周期 T。μ = 398600.4418 km³/s²（地球の重力定数 GM）。

$$\dot\Omega_{\text{SSO}} = +0.9856\ \text{deg/day} = \frac{360°}{365.25\ \text{day}}$$

太陽同期条件：dΩ/dt が地球の太陽公転と一致する値。i≈98°（LEO）で実現。

J2摂動と軌道変動

🙋

「J2摂動」って、よく宇宙の本で見るけど、何のことなんですか？人工衛星って、ケプラーの法則で楕円軌道をずっと回り続けるんじゃないんですか？

🎓

いい質問だね。ケプラーの法則は「地球が完全な球で、質量が中心に一点集中している」という理想的な条件でしか成り立たない。実際の地球は自転による遠心力で、赤道半径が極半径より約21km膨らんだ扁平な回転体なんだ。この「赤道のお腹のふくらみ」を球面調和関数で展開したときの2次項の係数が J2 = 1.08×10⁻³。月や太陽の引力よりも2〜3桁大きい、衛星軌道の最大の摂動源だよ。

🙋

そんなわずかな膨らみで、軌道がどう変わるんですか？

🎓

主に2つの効果が出る。一つ目は「軌道面そのものが歳差する」こと。地球を北極から見たとき、衛星軌道が赤道を横切る点（昇交点）が、コマの軸のように回転していくんだ。式で書くと dΩ/dt = -3/2·n·J2·(R_E/p)²·cos(i)。傾斜角 i の cos が効いているから、i<90°（順行軌道）では西へ歳差、i>90°（逆行軌道）では東へ歳差する。右のグラフで i=98° を選んでみて。dΩ/dt がちょうど +0.986 deg/day になっているはずだ。

🙋

+0.986 deg/day って…なんか中途半端な数字ですね。何か意味があるんですか？

🎓

これがまさに「太陽同期軌道」の鍵なんだ。地球は太陽の周りを365.25日で1周するから、360°÷365.25日 = 0.9856 deg/day。つまり、軌道面が1年でちょうど1周だけ東へ歳差すれば、衛星から見て太陽と軌道面の角度関係が常に一定に保たれる。これが SSO（Sun-Synchronous Orbit）だよ。LANDSAT や WorldView のような地球観測衛星はこの軌道を使って、世界中どこを撮っても「現地の朝10時半」みたいに同じ太陽高度で撮影できるんだ。影の長さが揃うから、画像比較・変化検出が桁違いにやりやすい。

🙋

なるほど！ 2つ目の効果は何ですか？dω/dt のグラフも左下のスライダーに連動して動いてますね。

🎓

2つ目は「近地点引数 ω の回転」、つまり楕円軌道の長軸そのものが軌道面内で回転していく現象だね。dω/dt = factor·(2 - 2.5sin²i) の式で、sin²i = 0.8 のとき、すなわち i = 63.4° または 116.6° で dω/dt = 0 になる。これを critical inclination（凍結軌道の傾斜角）と呼ぶよ。この傾斜角では楕円が回らないから、近地点をいつも南半球に固定したり、遠地点をいつも北極上空に置いたりできる。旧ソ連の Molniya 軌道（i=63.4°、e=0.74、周期12時間）が有名で、遠地点付近で衛星が「ゆっくり通過」する性質を使って、北極圏の通信中継に使われたんだ。

🙋

じゃあ凍結軌道は GPS や ISS とかでも使われているんですか？

🎓

GPS は半同期軌道（周期12時間、i=55°）で、ISS は i=51.6° の低傾斜角だから、どちらも凍結軌道ではないよ。ISS の場合は「ロシアのバイコヌール宇宙基地からロケット打ち上げできる最大傾斜角」が制約になっている。むしろ気象衛星や火星探査機の周回軌道で凍結軌道は重宝されていて、火星リコネサンス・オービター（MRO）もほぼ凍結軌道に設計されている。軌道力学は「自然の摂動を打ち消す」より「摂動を逆手に取ってミッションに使う」のが醍醐味なんだ。

よくある質問

地球は完全な球体ではなく、自転による遠心力で赤道半径が極半径より約21km膨らんだ扁平回転体（オブレートエリプソイド）です。この扁平を地球重力ポテンシャルの球面調和展開で表したときの2次帯状調和係数が J2 = 1.082635×10⁻³ で、月や太陽の摂動より2〜3桁大きい主要な摂動源です。J2項は人工衛星の軌道面を西または東に歳差させ（昇交点赤経の変動）、軌道楕円の長軸方向も回転させます（近地点引数の回転）。LEO衛星では1日あたり数度の変化が生じるため、軌道設計では必ず考慮します。

太陽同期軌道とは、軌道面の歳差率 dΩ/dt が地球の太陽公転と同じ +0.9856 deg/day になるよう設計された軌道です。これにより衛星は常に同じ地方時で各地点上空を通過でき、地球観測衛星にとって理想的な照明条件が得られます。dΩ/dt = -3/2·n·J2·(R_E/p)²·cos(i) の式から、軌道高度（半径）と傾斜角 i の組合せで +0.9856 deg/day が成立する条件を逆算します。高度500-800km の LEO では傾斜角 97-99° の逆行軌道になります（LANDSAT、WorldView、ALOS等）。本ツールでは「太陽同期との差」が0.05 deg/day 未満なら太陽同期と判定します。

凍結軌道は近地点引数 ω の変動率 dω/dt がゼロになる軌道で、軌道楕円の形状と向きが長期間維持される設計です。dω/dt = factor·(2 - 2.5sin²i) = 0 となる条件から、sin²i = 0.8、すなわち傾斜角 i = 63.4° または i = 116.6°（critical inclination）が解となります。この傾斜角ではJ2摂動による近地点回転が打ち消され、楕円軌道の長軸が同じ方向を向き続けます。代表例は旧ソ連の Molniya 軌道（i=63.4°、e≈0.74、12時間周期）で、近地点を南半球に固定して遠地点を北極上空に長時間滞在させ、高緯度通信に利用しました。

J2摂動式は半通径 p = a(1-e²) を使うため、円軌道（e=0）でも正常に計算できます。ただし e がきわめて小さいときは近地点引数 ω そのものが幾何学的に定義しづらく、ω の歳差を直接観測することは困難です。実務では円軌道に近い衛星でも数値積分で平均要素を追跡し、ω や昇交点経度などの「平均近地点離角」を扱います。また本シミュレーターは secular（永年）項のみを扱い、long-period（長周期）や short-period（短周期）の振動は含みません。実運用では SGP4 や数値積分で全項を考慮するのが標準です。

実世界での応用

地球観測衛星と太陽同期軌道：LANDSAT-8/9、Sentinel-2、WorldView、ALOS-2、Himawari（除く）、Planet社のSkySatなど、ほぼすべての光学地球観測衛星が太陽同期軌道を採用しています。高度約700km、傾斜角98.2°前後で、現地時刻10:00〜10:30に降交点を通過する設計が標準です。同じ太陽高度・同じ影の長さで撮影することで、地表変化検出や植生指数計算が安定します。J2摂動を逆手にとって自然に太陽同期を実現する、軌道力学の最高傑作の一つです。

Molniya軌道とロシア・北極圏通信：傾斜角63.4°、離心率0.74、周期12時間のMolniya軌道は、旧ソ連が高緯度（緯度60°以上）地域への通信中継に開発した代表的な凍結軌道です。遠地点が常に北半球上空にあり、ケプラー第2法則により衛星が遠地点付近をゆっくり通過するため、1日のうち8〜10時間も連続して通信可能です。GEO（赤道上空）からは見えない北極圏で、Molniya 1-Kからロシアの政府通信・テレビ放送中継に長年使用されました。

GPS・GNSS と軌道補正：GPS（半同期軌道、a≈26600km、i=55°）や Galileo、北斗、QZSS も J2摂動の影響を受け、長期的に Ω と ω が変動します。航法精度を維持するため、地上局から各衛星の軌道要素を毎日更新し、軌道暦（ephemeris）として配信しています。さらにスラスタによる軌道保持マニューバ（station keeping）も定期的に実施し、設計軌道面からのずれを補正しています。J2の効果を予測することは、GNSS時刻同期と位置決定の精度の根幹です。

火星・月探査ミッションの周回軌道設計：火星リコネサンス・オービター（MRO）や火星オデッセイは火星のJ2（火星も扁平）の効果を利用して凍結軌道に設計されており、近地点高度の自然な変動を最小化しています。月探査衛星（NASA Lunar Reconnaissance Orbiter）では月の重力場の不均一さ（mascon）がJ2より複雑な摂動を生み、頻繁な軌道修正が必要です。中央天体の重力場特性を理解することが、長期間の科学観測ミッションの成立条件です。

よくある誤解と注意点

まず最大の誤解が、「太陽同期軌道は常に同じ場所を撮影できる」というもの。SSO で揃うのは「地方時」と「太陽高度」だけで、各地点を直下点として撮影できる頻度は別問題です。例えば LANDSAT-8（高度705km、SSO）の再訪周期は16日。同じ地点を真上から撮るのは16日に1回しか巡ってきません。災害観測のように短時間で繰り返し撮影したい場合は、複数の SSO 衛星をコンステレーション化（Planet Labs、Capella Space など）するか、傾斜軌道の SAR 衛星を使います。SSO の利点は「比較しやすい画像」が撮れることであって、「いつでも撮れる」ではありません。

次に、「J2摂動は永久に同じ率で続く」という思い込み。本ツールが計算するのは secular（永年）項のみで、これは時間の1次関数で増える「平均的な変動」です。しかし実際には long-period（長周期）項として ω の周期がおよそ「2π / |dω/dt|」程度で振動し、short-period（短周期）項として軌道周期に同期した小さな振動も重畳します。例えば iSO 軌道で30日後の Ω 総変動は約29.6°ですが、ある時点のスナップショットでは±0.5°程度の short-period 振動を含む値になります。ミッション解析では SGP4/SDP4 や Cowell法による数値積分を使い、J2以外（J3, J4, ...）や大気抵抗・太陽輻射圧も含めて計算するのが標準です。

最後に、「critical inclination 63.4° なら近地点が完全に凍結する」という誤解。J2項に対しては確かに dω/dt = 0 になりますが、J3（地球重力場の3次帯状調和、洋ナシ型を表す項）が小さく無視できない値で残り、e と ω が長周期で結合振動します。真の凍結軌道を作るには i, e, ω の組合せを慎重に選び、J3 摂動による振動を中立点に置く設計が必要で、これは「Brouwer-Lyddane 凍結軌道」と呼ばれます。Molniya軌道も厳密には完全凍結ではなく、定期的なステーションキーピングで離心率と近地点引数を維持しています。摂動理論は近似であり、「completely frozen」を求めるなら高次項の評価が必須です。

J2摂動による軌道変動シミュレーター — 太陽同期軌道設計

J2摂動と軌道変動

よくある質問

実世界での応用

よくある誤解と注意点

使い方ガイド

具体的な計算例

実務での注意点