電離圏シンチレーション — 太陽活動が地球近傍のあらゆる衛星リンクを揺さぶるとき

この回で学ぶこと

太陽合の回では、探査機が太陽の背後に隠れる時期に電波が太陽コロナという乱流的なプラズマ媒質を貫くことで、群遅延の変動と位相シンチレーションが生じることを見ました。そこで導入した電子数積分量TEC(Total Electron Content)と、プラズマ分散関係から導かれる $\Delta\phi_{\text{corona}}\propto \text{TEC}/f$ という位相への効き方は、実は太陽コロナに限った話ではありません。電波が電離したプラズマを通過する限り、同じ物理が働きます。

そして私たちは、太陽コロナよりずっと身近な——というより、常にその中で衛星通信をしている——プラズマ媒質をもう1つ知っています。地球自身の電離圏(ionosphere)、高度およそ60kmから1000km程度に広がる、太陽紫外線・X線によって大気が部分的に電離した層です。太陽合が起きるのは探査機の軌道配置がたまたま太陽の裏に回り込む、数年に一度・数週間だけの特殊な幾何学的条件でしたが、電離圏は地球を回るすべての衛星、そして地上のすべてのGNSS受信機にとって、電波が地球に出入りするたびに必ず通過する、避けようのない媒質です。

この回では、太陽コロナで組み立てたTECとプラズマ分散の枠組みを、電離圏という空間スケール(太陽半径ではなく地球半径、光年ではなくキロメートル)も時間スケール(合の数週間ではなく、太陽フレア・CMEに応答する数分〜数日)もまったく異なる現象に適用し直します。特に、電離圏の電子密度が滑らかではなく不均一な構造(プラズマバブルなど)を持つときに生じる電離圏シンチレーション——振幅と位相の両方が乱される現象——を扱い、その強度を表す実務指標である振幅シンチレーション指数 $S_4$ を導入します。そして、これがGNSS測位のサイクルスリップや静止通信衛星のリンク設計にどう効いてくるかを見ていきます。

直感的導入: 「揺らぐガラス」は太陽の外側にもある

太陽合の回で、太陽コロナを「揺らぐガラス」に例えました。星がまたたくのは光そのものが揺らぐのではなく、途中の媒質(大気やプラズマ)の屈折率が不均一に揺らいでいるからだ、という比喩です。

この比喩を最後まで正直に辿ると、実は身近な星のまたたきの原因は太陽コロナではなく、地球自身の大気です。可視光であれば主に対流圏の温度・密度のゆらぎが原因ですが、電波(特にVHF〜L帯、GNSSが使う周波数帯を含む)にとって支配的な「揺らぐガラス」は、地球大気の中でもプラズマ状態にある層、すなわち電離圏です。つまり私たちは今回、太陽コロナという特殊な状況限定の擾乱源から、地球を回るすべての衛星・すべての地上局が常に付き合っている、より身近な擾乱源に話を移すことになります。

ここで区別しておくべき点が2つあります。1つは、マルチパスフェージングの回で扱った現象との違いです。マルチパスは、地面や建造物といった幾何学的な反射体によって生じる、複数経路の電波が干渉し合う現象でした。電離圏シンチレーションは反射ではなく、プラズマそのものの電子密度が空間的に不均一であることによる、屈折・回折の効果です。経路は本質的に1本のままですが、その1本の経路が通過する媒質の性質が場所ごとに揺らいでいる、という点で、太陽コロナの位相シンチレーションと同じ系統の物理現象です。

もう1つは、太陽コロナとの違いです。太陽コロナのTECは太陽からの距離(SEP角)にほぼ滑らかに依存する、比較的予測しやすい量でした。これに対し電離圏は、地球の自転・磁場・大気潮汐・太陽風との相互作用が絡み合う複雑な系であり、TECの空間分布そのものが**不均一な構造(irregularity)**を持ちます。この不均一構造こそが、電離圏を通過する電波に振幅・位相の両方の激しいランダム変動——シンチレーション——を引き起こす主因です。以下ではこれを定式化します。

電離圏TEC: 空間・時間スケールの違いを押さえる

太陽合の回で導入したTECの定義とプラズマ分散の関係式をそのまま流用します。視線に沿った電子密度の積分

$$\text{TEC} \equiv \int_{\text{視線}} n_e(l)\, dl \qquad [\text{個}/\text{m}^2]$$

と、それによる群遅延・位相前進

$$\Delta R_{\text{group}} \approx 40.3\,\frac{\text{TEC}}{f^2}\ \ [\text{m}], \qquad \Delta\phi \approx -\frac{2\pi\times 40.3}{c}\cdot\frac{\text{TEC}}{f}\ \ [\text{rad}]$$

という関係式は、プラズマ分散関係 $\omega^2=\omega_p^2+c^2k^2$ という媒質によらない普遍的な物理から出てくるものなので、太陽コロナだろうと地球電離圏だろうと形は変わりません。変わるのは、$n_e$ の値そのものと、それが分布する空間・時間スケールです。

太陽コロナのTECは、太陽表面近くを直撃する深い合で $10^{18}\text{–}10^{19}\ \text{個/m}^2$ のオーダーに達することもありましたが、これは太陽半径 $R_\odot\approx7\times10^8$ mという巨大なスケールにわたる積分の結果でした。一方、地球電離圏の鉛直TECは、通常時でおよそ $10^{16}\text{–}10^{18}\ \text{個/m}^2$ のオーダーで、電離圏物理では慣習的に $1\ \text{TECU} = 10^{16}\ \text{個/m}^2$ という単位を使います。積分区間はわずか高度60〜1000km程度、水平規模で言えば地球半径 $R_E\approx6{,}378$ km程度のスケールに収まる、太陽コロナよりずっとコンパクトな媒質です。

しかし積分区間が短くても電子密度そのものが桁違いに濃い(電離圏最大電子密度は $10^{11}\text{–}10^{12}\ \text{m}^{-3}$ 程度で、太陽風領域のコロナ外縁部よりずっと高密度)ため、TECの絶対値としては太陽合の穏やかな時期と同程度かそれ以上になり得ます。GPSのL1周波数($f=1.575$ GHz)で典型的な鉛直TEC(たとえば静穏時の日中で $20\ \text{TECU}=2\times10^{17}\ \text{個/m}^2$)による群遅延を見積もると、

$$\Delta R_{\text{group}} \approx 40.3\times\frac{2\times10^{17}}{(1.575\times10^9)^2} \approx 3.25\ \text{m}$$

となり、GNSS単独測位の誤差要因として無視できない大きさであることが分かります(これがGNSS受信機が2周波観測でこの遅延を較正する理由であり、太陽合の回で見た多周波数較正と同じ原理です)。ただし、この滑らかなTECによる測位バイアスは較正可能な系統誤差であり、これから扱う不均一構造によるシンチレーションとは性質が異なる問題です。

電子密度不規則構造とシンチレーションの物理

電離圏の電子密度は、大局的には滑らかに高度・緯度・時刻に依存しますが、その上に**空間スケール数十cm〜数kmの不均一な構造(irregularity)が重なっています。代表例が、赤道域の夜間に発生するプラズマバブル(equatorial plasma bubble)**で、下層の重い中性大気と上層の軽いプラズマの境界にレイリー・テイラー不安定性が働き、低密度のプラズマの塊が電離圏を上向きに突き抜けるように成長する現象です。極域では発生機構がまったく異なり、太陽風からのエネルギー流入(粒子の降り込み)とオーロラ電流系によって電子密度が細かく構造化されます。

これらの不均一構造を電波が通過すると、視線ごとに屈折率(したがって位相速度)がわずかに異なるため、波面が場所によって不規則に歪みます。この歪んだ波面が、電離圏から地上受信機までの距離 $z$ を伝搬する間に回折によって振幅の揺らぎに変換される——これが電離圏シンチレーションの物理的なメカニズムです。

この伝搬過程は、電離圏を「薄い位相スクリーン」とみなす弱散乱理論(Rytov近似)で扱うのが標準的です。位相スクリーンを通過した直後の位相擾乱 $\phi(\mathbf{x})$ の空間パワースペクトル密度を $\Phi_\phi(\kappa)$($\kappa$: 波面に垂直な面内の空間波数)とすると、地上までの距離 $z$ を伝搬したあとの対数振幅(強度の揺らぎに相当)のパワースペクトルは、

$$\Phi_\chi(\kappa) = \Phi_\phi(\kappa)\,\sin^2\!\left(\frac{\kappa^2 z}{2k}\right), \qquad k=\frac{2\pi f}{c}$$

という形になることが知られています。この $\sin^2(\cdot)$ の因子はフレネルフィルタと呼ばれ、波数 $\kappa$ がフレネルスケール $\kappa_F\sim\sqrt{k/z}$(等価的にフレネル半径 $r_F=\sqrt{\lambda z}$ の逆数程度)より十分小さい大規模構造($\kappa \ll \kappa_F$)からの寄与をほぼ遮断し、フレネルスケール前後の中小規模構造からの寄与を選択的に強調します。直感的には、大規模な不均一構造は波面全体をなだらかに傾ける(位相を変える)だけで、地上に届くまでに焦点を結ぶような回折は起こしにくいのに対し、フレネルスケール程度の細かい構造は、伝搬中に建設的・破壊的な干渉を起こして振幅そのものを激しく変動させる、という違いです。

この構造から、位相シンチレーションと振幅シンチレーションでは周波数依存性が異なることが導かれます。位相シンチレーションの分散 $\sigma_\phi^2$ は、$\Delta\phi\propto\text{TEC}/f$ の関係がそのまま揺らぎ成分にも当てはまるため、太陽コロナの場合と同じく

$$\sigma_\phi \;\propto\; \frac{1}{f}$$

というおおむね緩やかな周波数依存性を持ちます。一方振幅シンチレーションは、フレネルフィルタが波数 $k\propto f$ に依存する形で効いてくるため、位相よりも急な周波数依存性を持ちます。電離圏の不均一構造が典型的に持つべき乗型の空間スペクトル(観測的にはスペクトル指数 $p\sim2$–$3$ 程度がよく報告されます)のもとでフレネルフィルタを積分すると、振幅シンチレーションの標準偏差(あるいは後述の $S_4$ 指数)はおおむね

$$S_4 \;\propto\; f^{-p'} \qquad (p' \approx 1.5,\ \text{弱散乱領域の典型値})$$

という、位相より急な減衰則に従うことが弱散乱理論と観測の両面から知られています(この指数 $p'$ は不均一構造のスペクトル形状や幾何に依存する近似値であり、太陽コロナの節で見た $\text{TEC}\propto\psi^{-(\nu-1)}$ の指数 $\nu$ と同じく、媒質のモデルに依存するパラメータです)。位相は $1/f$、振幅は概ね $f^{-1.5}$ で減衰するという、この2つの異なる周波数依存性の存在が、後述する高周波数帯(Ka帯以上)でのシンチレーション対策の物理的根拠になります。

シンチレーション指数 $S_4$: 実務での標準指標

振幅シンチレーションの強さを定量化する実務標準の指標が、$S_4$ 指数です。受信強度(信号電力に比例する量)を $I(t)$ とし、その平均で正規化した揺らぎの標準偏差として定義されます。

$$\boxed{\ S_4 \equiv \sqrt{\dfrac{\langle I^2\rangle - \langle I\rangle^2}{\langle I\rangle^2}}\ }$$

$\langle\cdot\rangle$ は通常60秒程度の時間窓での平均を表します。$S_4=0$ はまったく揺らぎのない理想的な受信状態、$S_4$ が大きいほど振幅の落ち込み(フェード)が激しいことを意味します。弱いシンチレーション($S_4\lesssim0.3$)では受信強度の分布はおおむね対数正規分布や上で導いたRytov近似に従いますが、不均一構造が強く発達し多重散乱が卓越する強散乱(saturation)領域では、マルチパスフェージングの回で見たレイリー分布に近づき、$S_4$ は理論上1に漸近します(強い集束効果によって一時的に1をわずかに超える観測例も報告されています)。実務上は次のような目安で強度を分類することが多いです。

$S_4$ の範囲	目安
$S_4 < 0.3$	弱いシンチレーション
$0.3 \le S_4 < 0.6$	中程度
$S_4 \ge 0.6$	強いシンチレーション

位相側にも対応する指標があり、受信位相 $\phi(t)$ の(一定期間内での、緩やかなトレンドを除去した後の)標準偏差 $\sigma_\phi$(単位ラジアン、あるいは度)がよく使われます。GNSS受信機の分野では、しばしば1分間の位相標準偏差(15秒ごとにサンプリングした値の標準偏差)を $\sigma_\phi$ として報告する慣習があります。

地理的・時間的な発生傾向

電離圏シンチレーションは地球上どこでも一様に起きるわけではなく、はっきりとした地理的な偏りがあります。

赤道異常帯(equatorial anomaly region): 地磁気赤道の南北およそ $\pm15^\circ$–$20^\circ$ の帯状の領域で、特に日没後2〜3時間にプラズマバブルが発生しやすく、南米・アフリカ・東南アジアの一部地域では強いシンチレーション($S_4>0.6$)が頻繁に観測されます。これはレイリー・テイラー不安定性が日没直後の急激な電離圏電場の変化によって成長しやすいためです。

極域・オーロラ帯: 地磁気緯度でおよそ $60^\circ$–$75^\circ$ 付近のオーロラ帯や極冠域では、太陽風からのエネルギー流入や粒子降り込みによって電子密度が細かく構造化され、赤道域とは異なる発生機構でシンチレーションが起きます。極域のシンチレーションは磁気嵐の活動度と強く相関し、地磁気活動が高い時期に頻度・強度ともに増大します。

中緯度: 通常は電離圏が比較的滑らかで、シンチレーションの発生頻度は低いですが、大規模な磁気嵐が発生した際には中緯度でも一時的に強いシンチレーションが観測されることがあります。

これら地理的な傾向に加えて、太陽活動の11年周期との強い相関があります。太陽活動が活発な極大期には、太陽からのEUV(極端紫外線)放射が強く電離圏の電子密度そのものが底上げされるだけでなく、太陽フレアやCME(コロナ質量放出)によって引き起こされる磁気嵐の頻度も増えるため、赤道・極域いずれの発生機構も活性化し、シンチレーションの発生頻度・強度が顕著に増加します。逆に太陽活動極小期にはシンチレーション自体は大きく減少しますが、ゼロにはなりません。この11年周期は、太陽合の回で扱った太陽コロナの擾乱がSEP角という幾何学的条件に強く支配されていたのとは対照的に、太陽そのものの活動度に直接支配される時間依存性です。

実務での使われ方

GNSS測位精度とサイクルスリップ。 GNSS搭載航法の回で扱ったように、GNSS受信機は擬似距離とキャリア位相の連続的な追尾によって測位・軌道決定を行います。強い位相シンチレーションが発生すると、受信機の追尾ループが位相の急激な変動に追従しきれず、キャリア位相の整数値バイアスの連続性が失われるサイクルスリップが発生しやすくなります。これは太陽合の回でPLLのループ帯域幅とコロナ位相雑音のスペクトルの兼ね合いとして論じたのと同じ構造の問題で、電離圏の不均一構造の移動速度(プラズマバブルの東西ドリフト速度、典型的に $100$–$150\ \text{m/s}$ 程度)がシンチレーションの時間変動スケールを決め、それが受信機のトラッキングループ帯域幅と競合します。強いサイクルスリップが連続すると、精密単独測位(PPP)やRTK(リアルタイムキネマティック)測位のような高精度測位手法は数分〜数十分にわたって精度を大きく落とすか、測位そのものが一時的に不能になることがあります。米国連邦航空局(FAA)のWAAS(Wide Area Augmentation System)をはじめとする衛星航法補強システム(SBAS)は、赤道域のシンチレーションを電離圏の脅威モデル(ionospheric threat model)に明示的に組み込み、補強信号の完全性(インテグリティ)保証の設計に反映させています。

静止通信衛星でのKa帯採用の物理的根拠。 前節で見た通り、位相シンチレーションは $1/f$、振幅シンチレーション($S_4$)は概ね $f^{-1.5}$ という、より急な周波数依存性を持ちます。したがって、L帯(GNSS、$1$–$2$ GHz)やC帯・Ku帯を使う通信サービスと比べ、Ka帯(約$27$–$40$ GHz)以上の高周波数帯を使う静止通信衛星は、電離圏の不均一構造による振幅の落ち込みをはるかに小さく抑えられます。周波数比をたとえばL帯($1.5$ GHz)からKa帯($30$ GHz)まで20倍取ると、$S_4$ は理論上 $20^{-1.5}\approx1/89$ まで低減される計算になり、赤道域の強いシンチレーション環境下でもKa帯であればリンクマージンへの影響が限定的になります。ただしKa帯は降雨減衰の回で見た通り降雨減衰には弱く、電離圏シンチレーションへの耐性と降雨減衰への脆弱性というトレードオフの中で周波数帯が選ばれていることを押さえておく必要があります。

宇宙天気予報サービスの運用上の役割。 電離圏シンチレーションは太陽フレア・CMEという突発的な太陽活動に応答して発生・増減するため、事前の予測と即時の監視が運用上重要です。米国海洋大気庁の宇宙天気予報センター(NOAA SWPC)や欧州宇宙機関(ESA)の宇宙天気サービスネットワークは、太陽X線フラックス・太陽風パラメータ・地磁気指数(Kp指数など)をもとに電離圏擾乱・シンチレーションの発生確率を予報として発信しています。また米空軍研究所(AFRL)が運用してきたSCINDA(Scintillation Network Decision Aid)のような専用の地上観測網は、赤道域を中心にGNSS信号の $S_4$・$\sigma_\phi$ をリアルタイムで監視し、通信・測位運用者に警報を提供する役割を担っています。深宇宙探査機の運用で太陽合の時期にあらかじめ運用を切り替えるのと同様に、地球近傍のGNSS・通信衛星の運用者も、これらの宇宙天気予報を参照して、シンチレーションが激しくなる時間帯・地域でのミッションクリティカルな測位・通信オペレーションを避けるといった運用上の判断を行っています。

演習問題

1. 静穏時のGPS L1周波数($f=1.575$ GHz)での鉛直TECが $20$ TECU($1\ \text{TECU}=10^{16}$ 個/m$^2$)であるとき、式 $\Delta R_{\text{group}}\approx40.3\,\text{TEC}/f^2$ を用いて群遅延(距離換算)を計算してください。同じTECがKa帯($f=30$ GHz)の信号を通過した場合の群遅延も求め、両者を比較してください。

2. ある観測点で赤道域のプラズマバブル通過時にL帯で $S_4=0.7$(強いシンチレーション)が観測されたとします。振幅シンチレーションが弱散乱領域で近似的に $S_4\propto f^{-1.5}$ に従うと仮定し、同じ電離圏条件下でKa帯($L$帯の周波数の20倍)を使った場合の $S_4$ のおおよその値を見積もってください。またこの結果が、L帯を使うGNSSとKa帯を使う静止通信衛星とで、シンチレーションへの耐性がどう異なるかを説明してください。

3. 電離圏シンチレーションの発生に、(a) 赤道異常帯での日没後のプラズマバブル、(b) 極域・オーロラ帯での粒子降り込み、という2つの異なる物理機構があることを踏まえ、なぜ中緯度地域では一般にシンチレーションが起きにくいのか、この回で述べた発生機構の地理的偏りをもとに考察してください。

4. 太陽合の回でのコロナ位相雑音と、この回での電離圏位相シンチレーションは、どちらもPLL/GNSS受信機のトラッキングループにサイクルスリップのリスクをもたらす点で共通しています。両者が引き起こすトラッキング障害の「発生条件」がそれぞれ何に支配されているか(太陽合はSEP角という幾何学的条件、電離圏シンチレーションは何か)を比較し、運用上どちらが予測・回避しやすいと考えられるか、自分の言葉で論じてください。

まとめと次回予告

太陽合の回で組み立てたTECとプラズマ分散の枠組みは、太陽コロナだけでなく、地球自身の電離圏という別の空間・時間スケールを持つ媒質にもそのまま適用できることを見ました。電離圏の電子密度が滑らかではなく不均一な構造(プラズマバブルなど)を持つとき、電波の波面はフレネル回折を通じて振幅・位相の両方が乱される電離圏シンチレーションが生じ、その強さは振幅シンチレーション指数 $S_4$ で定量化されます。位相シンチレーションが太陽コロナと同じ $1/f$ の周波数依存性を持つのに対し、振幅シンチレーションはフレネルフィルタの効果でより急な $f^{-1.5}$ 程度の依存性を持ち、これが静止通信衛星でKa帯が有利になる物理的根拠になっていることを確認しました。またシンチレーションが赤道異常帯・極域で特に強く、太陽活動の11年周期と強く相関するという地理的・時間的な傾向、そしてGNSS測位のサイクルスリップという実務上の課題を、GNSS搭載航法の回の議論と接続する形で見ました。

次回は視点をさらに運用の現場に移し、複数の衛星・デブリが接近し合う際に必要となる衝突回避運用における通信要求に軽く触れます。今回見た宇宙天気による通信品質の変動が、まさにクリティカルな軌道変更コマンドを送る必要があるタイミングと重なったときにどう備えるか、という運用上の難しさを覗いていきます。

参考文献

• Yeh, K. C., and Liu, C. H., “Radio wave scintillations in the ionosphere,” Proceedings of the IEEE, 70(4), 324–360 (1982)
• Rino, C. L., “A power law phase screen model for ionospheric scintillation, 1. Weak scatter,” Radio Science, 14(6), 1135–1145 (1979)
• Kintner, P. M., Ledvina, B. M., and de Paula, E. R., “GPS and ionospheric scintillations,” Space Weather, 5, S09003 (2007)
• Basu, S., Basu, Su., et al., “Scintillations, plasma drifts, and neutral winds in the equatorial ionosphere after sunset,” Journal of Geophysical Research, 101(A12), 26795–26809 (1996)
• ITU-R Recommendation P.531, Ionospheric propagation data and prediction methods required for the design of satellite services and systems
• NOAA Space Weather Prediction Center, Ionospheric Scintillation (公開資料)
• Davies, K., Ionospheric Radio, Peter Peregrinus Ltd.