掩蔽電波科学 — 電波の屈折で他天体の大気を輪切りにする

この回で学ぶこと

対流圏屈折補正の回では、地球の中性大気が電波を曲げ、遅らせる現象を「測距・ドップラー観測に混入する系統誤差」として扱い、Saastamoinenモデルやマッピング関数を使ってそれを取り除く方法を学びました。また重力科学とドップラー計測の回では、航法の道具だったコヒーレントドップラーを、視点を転換して相手天体の重力場を読み取る科学計測器として使う発想を見ました。

この回は、その2つの発想を掛け合わせます。すなわち、「大気による屈折」という誤差要因を、今度は相手天体の大気を測るための信号そのものとして使うのです。探査機が地上局から見て惑星の背後に隠れる(掩蔽、occultation)直前・直後のわずか数分間、探査機からの電波は必然的にその惑星の大気・電離圏を斜めに貫いて地球に届きます。このとき電波は、地球大気で起きるのとまったく同じ物理——屈折率の高度勾配による経路の湾曲と伝搬の遅れ——によって曲げられ、その曲がりが受信周波数の微小なズレ(ドップラー残差)として観測されます。この残差の時系列から屈折角を求め、Abel変換という積分変換で大気の屈折率プロファイル $n(r)$ を復元し、さらに温度・圧力・電子密度の高度分布へと変換する。これが**掩蔽電波科学(radio occultation, 電波掩蔽)**です。

地球の対流圏補正では屈折率モデルは「与えるもの」でしたが、掩蔽観測では屈折率こそが「求めるもの」です。同じ数式が、符号を変えるように役割を反転させる様子に注目しながら読み進めてください。

直感的導入: 惑星の縁をかすめる電波は「大気の断層写真」を撮っている

探査機が惑星の周回軌道やフライバイ軌道を飛んでいて、その軌道が地球から見て惑星の裏側を通るとします。探査機が惑星の縁(リム)に近づき、隠れる直前の数分間を想像してください。探査機から地球へ向かう電波の視線は、惑星の縁すれすれをかすめ、時間とともにその最接近点(tangent point)がどんどん惑星表面へ向かって降りていきます。つまり視線は、大気の最上層から下層へと、大気を高度方向に輪切りに走査(スキャン)していくのです。これをイングレス(ingress、没入)掩蔽と呼びます。探査機が裏側から再び現れるときには逆に下層から上層へ走査され、これをエグレス(egress、出現)掩蔽と呼びます。

視線が大気を通過している間、電波には3つのことが起こります。

1. 屈折による湾曲: 大気の屈折率は下層ほど大きいので、電波の経路は惑星側に凸に曲がります。曲がった分だけ伝搬経路の幾何が変わり、受信周波数に真空伝搬の予測値からのズレ(ドップラー残差)が生じます。
2. 減衰(吸収・散乱): 大気分子による吸収(たとえば木星大気のアンモニア)や、土星の環の粒子による遮蔽で、受信信号の振幅が下がります。
3. 分散性の周波数依存: 電離圏(プラズマ)を通過する場合、対流圏(非分散性)とは異なり、屈折の大きさが周波数の2乗に反比例します。複数周波数で同時観測すれば、中性大気と電離圏の寄与を分離できます。

鍵になるのは1番目です。ドップラー残差の時系列 $\Delta f(t)$ は、視線の最接近点が大気のどの高度を通っているかに応じて刻々と変わります。高度ごとの屈折の強さが時系列に順番に「記録」されていくので、これを逆向きに解けば屈折率の高度分布が復元できるはずです。ちょうどCTスキャンが、体を貫いたX線の減衰量の集合から体内の断面像を再構成するのと同じ発想の逆問題であり、実際、これから見るAbel変換はCTで使われるRadon変換の球対称版に相当します。

数式定式化 1: 幾何光学近似での屈折角とドップラー残差

ブーゲの定理と屈折角の積分表示

対流圏屈折補正の回で導入したブーゲの定理をそのまま使います。天体中心からの距離 $r$、屈折率 $n(r)$ の球対称大気の中を進む光線に沿って、

$$n(r)\, r \sin\phi(r) = a = \text{一定}$$

が保存されます($\phi$ は動径方向と光線の進行方向のなす角)。この保存量 $a$ をインパクトパラメータと呼びます。光線が大気に入る前の直線経路を延長したときの、天体中心からの距離(を屈折率1で評価したもの)です。

光線が大気を貫通する間に進行方向が曲がる総角度、すなわち屈折角(bending angle) $\alpha$ は、光線の経路に沿った屈折率勾配の積算として、次の積分で与えられます(幾何光学近似、Fjeldbo, Kliore & Eshleman 1971)。

$$\alpha(a) = -2a \int_{r_t}^{\infty} \frac{1}{\sqrt{n^2 r^2 - a^2}}\, \frac{d \ln n}{dr}\, dr$$

ここで $r_t$ は光線の最接近点の半径で、ブーゲの定理より $n(r_t)\, r_t = a$ を満たします。積分の中身を眺めると、(1) 屈折率勾配 $d\ln n/dr$ が大きい高度ほど寄与が大きい、(2) 分母は最接近点 $r=r_t$ でゼロに近づくため、光線は自分の最接近点付近の大気の情報を最も強く拾う、という2つの性質が読み取れます。この第2の性質が「視線の最接近点の降下=高度方向の走査」という直感を数学的に裏付けています。屈折率が高度とともに単調減少する通常の中性大気では $d\ln n/dr < 0$ なので $\alpha > 0$、つまり光線は天体側に曲がります。

観測量はドップラー残差

屈折角 $\alpha$ 自体は直接測れません。地上局が測るのは受信周波数です。探査機と地上局の位置・速度は軌道決定によって高精度に分かっているので、「大気がなかった場合(真空伝搬)」の受信周波数の予測値 $f_{vac}(t)$ が計算できます。実際の観測値との差

$$\Delta f(t) \equiv f_{obs}(t) - f_{vac}(t)$$

が掩蔽観測の生の観測量です。重力科学の回のドップラー残差と定義はまったく同じですが、今回は残差を生む犯人が重力場ではなく大気の屈折です。

電波の経路が屈折角 $\alpha$ だけ曲がると、送信点・受信点での電波の射出・到来方向が真空中の直線から傾き、視線速度への射影が変わることで周波数がずれます。掩蔽面内の幾何で整理すると、小さな $\alpha$ に対して残差は近似的に

$$\Delta f \approx \frac{f}{c}\, v_\perp\, \alpha$$

と書けます。ここで $v_\perp$ は探査機・地上局の相対速度のうち、掩蔽面内で視線に垂直な成分(最接近点を大気中で降下・上昇させる速度に対応)です。実務では、この近似式ではなく、既知の位置・速度ベクトルとブーゲの定理(送信側・受信側の両端で $a$ が共通)を連立させた厳密な幾何方程式を各時刻について解き、観測値 $\Delta f(t_i)$ の一つ一つを組 $(\alpha_i, a_i)$ に変換します。こうして観測の時系列から屈折角プロファイル $\alpha(a)$ が得られます。

Abel変換による逆算

順問題は「$n(r)$ が与えられたら $\alpha(a)$ が積分で求まる」でした。掩蔽観測で解きたいのはその逆、「$\alpha(a)$ の観測から $n(r)$ を求める」逆問題です。幸い、球対称を仮定したこの積分方程式はAbel型と呼ばれる形をしており、**厳密な逆変換(Abel変換)**が存在します。

$$\ln n(a_1) = \frac{1}{\pi} \int_{a_1}^{\infty} \frac{\alpha(a)}{\sqrt{a^2 - a_1^2}}\, da$$

観測された屈折角プロファイル $\alpha(a)$ をこの式に代入して数値積分すれば、各インパクトパラメータ $a_1$ に対する屈折率 $n(a_1)$ が求まり、対応する実際の高度は $r_t = a_1 / n(a_1)$ で復元できます。積分が $a_1$ から無限大まで(すなわち、その光線より上空を通ったすべての光線について)走ることに注意してください。ある高度の屈折率を知るには、それより上層を通過したすべての光線の観測が必要——だからこそ掩蔽観測は、大気の外から表面(または電波が届かなくなる深さ)まで、走査を途切れなく続けることに意味があるのです。

モデルフィッティング(パラメータの最小二乗推定)に頼らず、観測から直接プロファイルが「解ける」という点が、掩蔽観測の際立った強みです。仮定は球対称性(および幾何光学近似が成り立つこと)だけであり、大気の組成や温度構造についての事前モデルを必要としません。

数式定式化 2: 屈折率から大気パラメータへ

中性大気: 数密度・圧力・温度プロファイル

復元された屈折率プロファイルを物理量に変換します。中性大気の屈折率指数 $N = (n-1)\times 10^6$ は分子の分極による寄与であり、分子数密度 $\nu(r)$ に比例します。

$$N(r) = \kappa\, \nu(r)$$

比例係数 $\kappa$ は大気の組成で決まる定数で(分子種ごとの屈折体積の重み付き和)、木星・土星なら水素・ヘリウムの混合比、火星・金星ならCO$_2$主体、タイタンならN$_2$主体として評価します。ここが掩蔽観測に必要な数少ない事前情報です。組成が分かれば $N(r)$ から数密度 $\nu(r)$ が直ちに得られます。

数密度が分かれば、静水圧平衡と理想気体の状態方程式から温度・圧力が芋づる式に求まります。静水圧平衡 $dP/dr = -\nu(r)\, \bar m\, g(r)$($\bar m$: 平均分子質量、$g$: 重力加速度)を上端 $r_{top}$ から積分すると、

$$P(r) = P(r_{top}) + \int_{r}^{r_{top}} \nu(r')\, \bar m\, g(r')\, dr'$$

理想気体 $P = \nu k_B T$ より、

$$T(r) = \frac{P(r)}{k_B\, \nu(r)}$$

となり、温度・圧力の高度プロファイルが得られます。上端の境界条件 $P(r_{top})$ には不定性が残りますが、十分高い高度から積分を始めれば、下層に向かうにつれて積分項が支配的になり境界条件の影響は指数関数的に薄れていきます。ドップラー残差の測定という「周波数の計測」だけから、他天体の温度計・気圧計を高度分解能つきで実現してしまうのが、この手法の鮮やかさです。

電離圏: 電子密度プロファイル

天体の電離圏をプラズマとして扱うと、周波数 $f$ の電波に対する屈折率は

$$n \approx 1 - \frac{f_p^2}{2f^2}, \qquad f_p^2 = \frac{n_e\, e^2}{4\pi^2 \varepsilon_0 m_e}$$

で与えられます($f_p$: プラズマ周波数、$n_e$: 電子密度)。数値を入れて整理すると

$$n - 1 = -\frac{40.3\, n_e}{f^2}$$

($n_e$ は m$^{-3}$、$f$ は Hz)という便利な形になります。中性大気と比べたときの重要な違いは2つです。

1. 符号が負: 電離圏では屈折率が1より小さく、高度勾配も中性大気と逆向きになり得るため、屈折角の符号も反転します。掩蔽の時系列で、まず電離圏由来の小さな負の(または正負に振れる)残差が現れ、続いて中性大気由来の大きな正の残差が現れる、という特徴的なシグネチャが観測されます。
2. 分散性($f^{-2}$ 依存): 2周波数以上で同時観測して残差を比較すれば、$f^{-2}$ でスケールする成分(電離圏)としない成分(中性大気)を分離できます。これは対流圏補正の回で見た「電離層遅延の2周波分離」とまったく同じ論法を、今度は相手天体のプラズマに適用したものです。

Abel変換の枠組みはそのまま使えます。電離圏起源の屈折角プロファイルを逆変換して $n(r)$ を求め、上の式を $n_e$ について解けば、電子密度の高度プロファイル $n_e(r)$ が得られます。火星や金星の電離圏の主要な観測データは、歴史的にこの手法によって蓄積されてきました。

幾何学的制約と軌道設計: 数分間しかない観測窓

掩蔽観測には、他の観測手法にはない厳しい幾何学的制約があります。

観測できるのは掩蔽の起きる瞬間だけです。探査機・天体・地球(受信局)が一直線に並ぶ配置は軌道上の特定の時刻にしか実現せず、しかも視線の最接近点が大気(たとえば厚さ数百km)を通過するのは、イングレス・エグレスそれぞれで数十秒〜数分間にすぎません。最接近点がリムを横切る速度は探査機の軌道速度と幾何で決まり、たとえば毎秒数kmでリムに降りていくなら、200〜300kmの大気を貫く時間は1〜2分です。この短い窓の中に、電離圏から下層大気までの全プロファイルの情報が圧縮されて届きます。

この制約は運用にも波及します。イングレスの途中で探査機は天体の陰に入り、アップリンクが先に途絶します。したがって掩蔽の科学観測は多くの場合、探査機搭載のUSO(超高安定発振器)を周波数基準とする1-wayダウンリンクで行わざるを得ません。トランスポンダによるコヒーレント2-wayが使えないため、観測精度はUSOの周波数安定度(カッシーニ搭載USOでアラン偏差 $\sim 2\times10^{-13}$)が直接律速します。ドップラー残差に要求される分解能はmHz〜数十mHz(Xバンドで $\Delta f/f \sim 10^{-12}$)のオーダーなので、USOの安定度はぎりぎりの綱渡りであり、掩蔽電波科学はUSO技術の最も要求の厳しいユーザーの一つです。

さらに、掩蔽が起きるかどうか・どこで起きるかは軌道設計そのものです。フライバイなら、B平面上の照準点を選ぶことで地球掩蔽の有無・掩蔽される緯度が決まります。周回機なら、軌道面と地球方向の幾何が掩蔽の季節(掩蔽シーズン)を決め、軌道周回ごとにイングレス・エグレスの緯度が少しずつ移動していきます。カッシーニは土星本体・環・タイタンの掩蔽を意図的に起こすよう周回軌道の設計を繰り返し変更し、17年のミッションを通じて広い緯度範囲のプロファイルを蓄積しました。掩蔽電波科学は「観測機器」をほとんど追加せずに済む(通信系そのものが計測器になる)代わりに、軌道設計と運用計画の中に観測を書き込んでおく必要がある手法だと言えます。

実務での使われ方

• **マリナー4号(NASA/JPL、1965年)**による火星フライバイが、惑星探査における電波掩蔽の最初の実施例です。Sバンド電波の掩蔽から火星の地表気圧がわずか4〜7hPa(地球の1%以下)しかないことが判明し、それまでの推定(数十hPa)を覆して、その後の火星着陸機設計(パラシュート設計など)の前提を書き換えました。掩蔽電波科学が「通信系がそのまま科学機器になる」ことを示した歴史的な観測です。
• **ボイジャー1号・2号(NASA/JPL、1979–1989年)**は、木星・土星・天王星・海王星の4巨大惑星すべてとタイタンで掩蔽観測を実施しました。S/Xバンド2周波の残差から中性大気と電離圏を分離し、各惑星の温度—圧力プロファイル(約1hPa〜1bar領域)と電離圏の電子密度分布を初めて系統的に決定しています。ボイジャー1号のタイタン掩蔽は、地表気圧約1.5bar・地表温度約94Kという値を確定させ、タイタンが厚い窒素大気を持つ天体であることを裏付けました。
• カッシーニ(NASA/ESA/ASI、2004–2017年)の電波科学サブシステム(RSS)は、S(2.3GHz)/X(8.4GHz)/Ka(32GHz)の3周波同時ダウンリンクを搭載USO基準で送信できる、掩蔽観測として最も充実した構成でした。土星本体の掩蔽に加えて、土星の環による掩蔽では、電波が環を貫通する際の減衰(光学的深さ)と前方散乱を3波長(13cm・3.6cm・0.94cm)で比較することで、環の粒子サイズ分布(おおむねcm〜m級の粒子が支配的であること)や環の微細構造(密度波・小衛星の作る隙間)を半径方向分解能100m級で明らかにしました。波長より十分大きい粒子はどの波長も同じように遮るのに対し、波長と同程度の粒子は波長ごとに減衰が変わる——回折の基本性質を逆手に取った粒径測定です。
• 火星周回機による電離圏・大気観測: Mars Global Surveyor、Mars Expressなどは多数回の地球掩蔽を積み重ね、火星電離圏の電子密度プロファイル(主ピークは高度約130km、太陽天頂角に依存して $n_e \sim 10^{11}$ m$^{-3}$ 級)と下層大気の温度構造の気候学的データセットを構築しました。火星電離圏の知識の骨格は掩蔽観測によって作られたと言ってよく、これは将来の火星通信・航法インフラ(電離圏遅延補正)の設計にも直結する工学データでもあります。
• GPS掩蔽(GNSS-RO)との比較: 同じ原理を地球大気に適用したのがGNSS電波掩蔽で、低軌道衛星がGPS衛星の電波を大気越しに受信します。台湾・米国の**COSMIC(FORMOSAT-3、2006年)は6機編隊で1日約2,000回、後継のCOSMIC-2(2019年)**は1日約4,000回以上の掩蔽プロファイルを取得し、数値気象予報(NWP)へリアルタイムで同化されています。惑星掩蔽が「1回の掩蔽イベントを数ヶ月がかりで設計する」世界であるのに対し、GNSS-ROは「送信機(GNSS衛星)が常に多数飛んでいる」ため掩蔽が毎日大量に起きる——同じ物理・同じAbel変換でありながら、送信機インフラの有無が観測頻度を6桁近く変えるという、対照的で示唆的な実例です。

演習問題

1. 火星電離圏の主ピークの電子密度を $n_e = 2\times10^{11}\ \text{m}^{-3}$ とする。$n - 1 = -40.3\, n_e / f^2$ を用いて、Sバンド($f=2.3$ GHz)とXバンド($f=8.4$ GHz)それぞれに対する屈折率のずれ $n-1$ と屈折率指数 $N=(n-1)\times10^6$ を計算せよ。また両者の比が $(8.4/2.3)^2$ に一致することを確かめ、この周波数依存性が中性大気と電離圏の寄与の分離にどう使えるか説明せよ。

2. 近似式 $\Delta f \approx (f/c)\, v_\perp\, \alpha$ を用いる。GNSS掩蔽(GPS L1、$f=1.575$ GHz)で、掩蔽面内の垂直相対速度を $v_\perp = 2.7$ km/s、地球下層大気での屈折角を $\alpha = 0.02$ rad(約1.1°)とするとき、ドップラー残差 $\Delta f$ を求めよ。次に、同じ式を火星掩蔽(Xバンド $f=8.4$ GHz、$v_\perp = 3$ km/s、火星の希薄な大気での $\alpha = 2\times10^{-4}$ rad)に適用して $\Delta f$ を求め、両者を比較せよ。

3. ある土星掩蔽で、視線の最接近点がリムに垂直な方向に 2.0 km/s の速さで降下していくとする。電離圏上端(高度3,000km相当)から、電波が深部大気の屈折・吸収で追跡不能になる高度(気圧約1barの高度)までの走査幅を3,200kmとするとき、この掩蔽の観測窓の長さを求めよ。また、この間ドップラー残差を0.1秒ごとにサンプリングすると、高度方向の生のサンプリング間隔は何mになるか。

4. イングレス掩蔽の進行中、探査機へのアップリンクは途中で途絶するため、掩蔽電波科学は搭載USOを基準とする1-wayダウンリンクで行われるのが通例である。(a) 2-wayコヒーレント運用が使えれば地上局の水素メーザー(安定度 $\sim10^{-15}$)を基準にできるのに対し、1-way運用ではなぜUSOの安定度($\sim2\times10^{-13}$)が観測精度を直接律速するのか。(b) Xバンド(8.4GHz)でUSO安定度 $2\times10^{-13}$ が生む周波数の揺らぎは何mHzか。これを本文で述べたドップラー残差の要求分解能と比較し、掩蔽電波科学がUSOに要求する性能水準について論ぜよ。

まとめと次回予告

この回では、対流圏屈折補正で「除去すべき誤差」だった大気屈折を反転させ、相手天体の大気を測る信号として使う掩蔽電波科学を学びました。ブーゲの定理から屈折角の積分表示を導き、観測されるドップラー残差を $(\alpha, a)$ の組に変換し、球対称仮定の下でAbel変換により屈折率プロファイル $n(r)$ をモデルフィッティングなしで直接復元できること、そして中性大気では静水圧平衡経由で温度・圧力プロファイルに、電離圏では $f^{-2}$ の分散性を利用して電子密度プロファイルに変換できることを見ました。観測窓がイングレス・エグレスの数分間しかないという幾何学的制約が、1-way/USO運用や軌道設計への組み込みという運用上の特徴を生むことも確認しました。マリナー4号の火星気圧測定からカッシーニの環の粒径分布、そして毎日数千プロファイルを気象予報に流し込むGNSS-ROまで、「通信系がそのまま科学機器になる」というこの手法の射程は、測距・追跡カテゴリの締めくくりにふさわしい広がりを持っています。

測距・追跡カテゴリはこの回で一区切りです。次回はシステム・運用カテゴリに移り、地上局や探査機の内部で高周波信号を導く導波管のモード理論を扱います。金属の管の中を電磁波がどのような固有モード(TE/TMモード)で伝わるのか、遮断周波数はどこから来るのか——自由空間の伝搬とはまた違った、境界条件が支配する電磁気学の世界に入っていきます。

参考文献

• G. Fjeldbo, A. J. Kliore, V. R. Eshleman, “The Neutral Atmosphere of Venus as Studied with the Mariner V Radio Occultation Experiments,” The Astronomical Journal, 76(2), 123–140, 1971
• A. J. Kliore et al., “Occultation Experiment: Results of the First Direct Measurement of Mars’s Atmosphere and Ionosphere,” Science, 149(3689), 1243–1248, 1965
• G. F. Lindal et al., “The Atmosphere of Titan: An Analysis of the Voyager 1 Radio Occultation Measurements,” Icarus, 53(2), 348–363, 1983
• E. A. Marouf, G. L. Tyler, P. A. Rosen, “Profiling Saturn’s Rings by Radio Occultation,” Icarus, 68(1), 120–166, 1986
• E. R. Kursinski et al., “Observing Earth’s Atmosphere with Radio Occultation Measurements Using the Global Positioning System,” Journal of Geophysical Research, 102(D19), 23429–23465, 1997
• A. J. Kliore et al., “Cassini Radio Science,” Space Science Reviews, 115, 1–70, 2004
• DSN Telecommunications Link Design Handbook, DSN No. 810-005