測距・追跡#189

惑星レーダー天文学 — 地上の巨大アンテナで天体そのものを測る

レーダー高度計が探査機から足元を測る技術だったのに対し、今回は地上のDSN 70m級アンテナから小惑星や金星に電波を撃ち、そのエコーから距離・自転・表面性状を読み取る惑星レーダー天文学。距離の4乗で減衰するレーダー方程式の壁と、遅延-ドップラーマッピングによる表面再構成の幾何学を数式で理解する。

前提知識: レーダー高度計 — 探査機が自分の足元までの距離を測る

惑星レーダーレーダー方程式遅延ドップラーマッピング惑星防衛DSN

この回で学ぶこと

前回扱ったレーダー高度計は、探査機に搭載した送受信機で、自分の真下数百km以下の「足元」を測る技術でした。この回で扱う 惑星レーダー天文学(planetary radar astronomy) は、電波を出す場所と測る対象がちょうど逆転します。すなわち、地球上の巨大アンテナ——DSN(Deep Space Network)のゴールドストーン70mアンテナや、かつてのアレシボ305m電波望遠鏡——から、数百万kmから数千万km彼方の小惑星・金星・月・水星に向けて数百kWの電波を送信し、天体表面で散乱されて返ってくる微弱なエコーを受信するという観測形態です。

探査機を送らずに、地上に居ながらにして天体の

距離(往復光行時間から、探査機測距と同じ原理で)
視線速度と自転(エコーのドップラーシフトとドップラー広がりから)
形状・表面性状(遅延-ドップラー2次元マップと偏波・反射率から)

を測定できる、地上局技術の科学応用の極致とも言える分野です。この回では、まず通信リンクとの決定的な違いである距離の4乗減衰をレーダー方程式から定量的に導き、「なぜ地球近傍小惑星は現実的な観測対象で、遠方の惑星は絶望的に困難なのか」を数値で確認します。次に、この分野の中心的な信号処理技術である 遅延-ドップラーマッピング(delay-Doppler mapping) を取り上げ、時間(遅延)と周波数(ドップラー)という2つの軸だけから、望遠鏡では点にしか見えない天体の表面を2次元画像として再構成する幾何学を定式化します。

なお、送受信のハードウェアも波形設計(PN符号による測距や整合フィルタ・パルス圧縮)も、これまでのレッスンで学んだ深宇宙通信の技術とほぼ地続きです。実際、ゴールドストーンの惑星レーダーはDSNの通信・測距インフラそのものを使っています。違うのは、相手が応答してくれる探査機のトランスポンダではなく、**ただ黙って電波を散乱するだけの自然天体(非協力ターゲット)**だという一点であり、この一点がすべての困難の源になります。

直感的な導入: トランスポンダのいない測距

探査機の測距では、地球局が送った信号を探査機のトランスポンダが増幅して送り返してくれました。だからこそ往路と復路のそれぞれで距離の2乗ずつ、合計しても「往路 $1/R^2$ × 復路 $1/R^2$ 」の減衰を、途中の再増幅で一度リセットできたわけです。

惑星レーダーにはこの中継増幅がありません。送信電波は距離 $R$ を伝搬して天体に到達するまでに $1/R^2$ で薄まり、天体表面はそのごく一部を全天に散乱し、その散乱波がまた $1/R^2$ で薄まりながら地球に戻ってきます。往復で $1/R^4$ 。この「4乗」が惑星レーダーの物理を支配します。距離が10倍になると受信電力は1万分の1、100倍になると1億分の1です。

それでも観測が成立するのは、送り手と受け手が地球上の最強クラスの設備だからです。ゴールドストーンのDSS-14(70m、X帯8560MHz、送信電力450kW)、そして2020年に崩落事故で失われるまで世界最強の惑星レーダーだったアレシボ天文台(305m、S帯2380MHz、送信電力約900kW)。数百kWの連続波を口径数十〜数百mのアンテナで細いビームに絞って撃ち、返ってきた $10^{-21}$ W オーダーのエコーを極低雑音受信機と長時間積分で拾い上げる——リンクバジェットの感覚で言えば、通信では絶対に成立させないような「赤字すれすれ」のリンクを、科学のために力ずくで成立させているのが惑星レーダーです。

レーダー方程式の天体スケール適用

導出

送信電力 $P_t$ 、送信アンテナ利得 $G_t$ とすると、距離 $R$ にある天体の位置での電力束密度(単位面積あたりの電力)は

S_{\text{target}} = \frac{P_t G_t}{4\pi R^2}

天体はこのうち レーダー断面積(radar cross section) $\sigma$ [m²] に相当する電力を捕まえて等方的に再放射する、と定義します( $\sigma$ は幾何学的断面積 $\pi R_a^2$ に反射率などを掛けた実効値で、天体の物理断面積そのものではない点に注意)。地球に戻ってきたエコーの電力束密度は

S_{\text{earth}} = \frac{P_t G_t}{4\pi R^2} \cdot \frac{\sigma}{4\pi R^2}

これを実効開口面積 $A_e = G_r \lambda^2 / (4\pi)$ の受信アンテナで受けると、受信電力は

\boxed{\;P_r = \frac{P_t G_t G_r \lambda^2 \sigma}{(4\pi)^3 R^4}\;}

これが レーダー方程式 です。分母に $R^4$ が現れていること、そして通信リンクの受信電力式(フリスの伝達公式、分母は $R^2$ )と比べて $\sigma/(4\pi R^2)$ という因子が余分に掛かった構造になっていることを確認してください。トランスポンダによる再生中継が、いかに $R^4$ を $R^2$ に「折り返して」くれていたかが分かります。

SNRと積分時間

検出可能性を決めるのは受信電力そのものではなくSNRです。受信系の雑音温度を $T_s$ 、エコーの帯域幅を $B$ 、コヒーレント/インコヒーレント積分をあわせた実効積分時間を $T$ とすると、おおまかに

\text{SNR} \propto \frac{P_t G_t G_r \lambda^2 \sigma}{(4\pi)^3 k T_s R^4} \sqrt{\frac{T}{B}}

という比例関係になります(インコヒーレント積分では $\sqrt{T}$ でしか稼げない点が探査機測距との違いです。エコーは天体の自転による位相の乱雑化のため、無限にコヒーレント積分することができません)。積分時間 $T$ は最大でも天体が地球近傍にいる期間・1晩の可視時間で制限されるので、SNRの主導権はほぼ完全に $R^4$ が握っています。

数値で見る「4乗の壁」

具体的に比較しましょう。天体のレーダー断面積は直径の2乗に比例する( $\sigma \sim \hat{\sigma}\pi R_a^2$ 、 $\hat{\sigma}$ は規格化レーダーアルベドで岩石質天体なら0.1前後)ことに注意して、同じ設備で観測した場合の相対SNRを考えます。

ターゲット	典型距離 $R$	相対 $1/R^4$ (金星最接近=1)
接近中の地球近傍小惑星	0.05 AU ( $7.5\times10^6$ km)	$\sim 10^{3}$
月	0.0026 AU ( $3.8\times10^5$ km)	$\sim 10^{9}$
金星(最接近)	0.28 AU ( $4.2\times10^7$ km)	$1$
火星(最接近)	0.52 AU	$\sim 8\times10^{-2}$
木星系(ガリレオ衛星)	4.2 AU	$\sim 2\times10^{-5}$
土星系(タイタン・環)	8.5 AU	$\sim 1\times10^{-6}$

金星最接近を基準にすると、土星系は6桁落ちです。それでもアレシボは1970年代にガリレオ衛星の、1980年代末〜2000年代に土星の環とタイタンのエコー検出に成功していますが、得られるのは円盤全体を積分したスペクトル程度で、後述する精細なマッピングは不可能でした。逆に、直径わずか数百mの小惑星でも、0.05 AUまで近づいてくれれば金星より3桁有利な $1/R^4$ が断面積の小ささ( $\sigma$ で6〜8桁の不利)をかなり相殺してくれます。惑星レーダーの主戦場が「たまたま近くを通る地球近傍小惑星」になっているのは、この $R^4$ の算術の帰結です。

もうひとつ、近距離ターゲットには運用上の面白い問題があります。往復光行時間は金星最接近で約280秒あるので、1本のアンテナで「送信→切替→受信」が間に合います。ところが月(往復2.6秒)や極端に接近した小惑星では、送信を止めて受信に切り替える前にエコーが返ってきてしまう。このため、ゴールドストーンで送信しグリーンバンク100m望遠鏡等で受信するバイスタティック(2局)構成が常用されます。送受を別サイトに分けるのは、地上局運用で学んだ2-way/3-way測距の構図とよく似ています。

遅延-ドップラーマッピング

問題設定: 点にしか見えない天体を「画像化」する

直径1kmの小惑星が $10^7$ km彼方にあるとき、その視直径は $10^{-4}$ ラジアンの1000分の1、約0.02秒角です。これは巨大アンテナのビーム幅(X帯70mで約0.03度)よりはるかに小さく、レーダービームは天体全体をすっぽり覆ってしまいます。つまり角度方向の分解能では天体表面を分解できません。

そこで発想を転換します。エコーには角度以外に2つの測定可能な軸があります。

遅延(delay): 天体の手前側(地球に近い側)からのエコーは早く、縁(リム)に近い部分からのエコーは遅く返ってくる。
ドップラー(Doppler): 天体は自転しているので、自転によって地球に近づいてくる側の表面からのエコーは正に、遠ざかる側は負にドップラーシフトする。

受信エコーを遅延ビンとドップラービンの2次元ヒストグラムに整理したものが 遅延-ドップラーマップ です。以下でこの2軸の幾何学を定式化します。

等遅延線: 同心円

半径 $R_a$ の球形天体を考え、地球方向をz軸に取ります。天体表面で最も地球に近い点(サブレーダー点)を基準にすると、そこから奥行き $z$ だけ深い位置にある表面素片のエコーの相対遅延は

\Delta\tau = \frac{2z}{c}

です(往復なので係数2。レーダー高度計の $R = c\tau/2$ と同じ関係です)。奥行き $z$ が等しい点の集合は、視線に垂直な平面で球を切った切り口、すなわち**サブレーダー点を中心とする同心円(等遅延リング)**になります。遅延はサブレーダー点の $\Delta\tau = 0$ からリムの $\Delta\tau = 2R_a/c$ まで分布し、天体全体の「遅延の深さ」は直径1kmの小惑星で約6.7μsです。送信波形に帯域幅 $B_w$ のPN符号変調やチャープを使えば、遅延分解能 $\sim 1/B_w$ 、すなわち距離分解能

\Delta z = \frac{c}{2B_w}

が得られます。これは前回導いた式そのものです。ゴールドストーンやアレシボの惑星レーダーは数MHz〜数十MHzの変調帯域を使い、 $10^7$ km彼方の天体表面を数m〜数十m刻みの遅延リングに分解します。

等ドップラー線: 自転軸に平行な直線

次にドップラー軸です。天体が角速度 $\Omega = 2\pi/P$ ( $P$ は自転周期)で自転しているとします。視線に垂直な平面(天球面)に自転軸を投影し、その投影軸からの(視線と自転軸の両方に垂直な方向の)距離を $x$ とすると、表面素片の視線方向速度は

v_{\text{los}} = \Omega\, x \cos\delta

( $\delta$ は自転軸が視線に垂直な状態からのずれ、いわゆるサブレーダー点の自転緯度に相当する傾き)。エコーのドップラーシフトは往復分で

\Delta f = \frac{2 v_{\text{los}}}{\lambda} = \frac{4\pi \cos\delta}{\lambda P}\, x

となり、 $x$ だけの1次関数です。したがって等ドップラー線は「投影された自転軸に平行な直線群」になります。ドップラーシフトはリム( $x = \pm R_a$ )で最大となり、エコー全体のドップラー広がり(帯域)は

B_{\text{echo}} = \frac{8\pi R_a \cos\delta}{\lambda P}

数値例: 直径1km( $R_a = 500$ m)、自転周期 $P = 5$ 時間、 $\cos\delta = 1$ 、X帯 $\lambda = 3.5$ cm なら $B_{\text{echo}} \approx 20$ Hz。エコーは搬送波の周りわずか±10Hzに収まる極めて狭帯域の信号で、周波数分解能 $\Delta f = 1/T_{\text{coh}}$ ( $T_{\text{coh}}$ はコヒーレント積分時間)でFFTすれば、数十秒の積分で0.01Hz級、つまり天体表面を数百本のドップラーストリップに分解できます。逆に、この式はエコーの帯域幅を測るだけで、形状を分解できなくても $R_a/P$ が分かることも意味します。後述する金星の自転周期決定はまさにこの原理によるものでした。

2軸を重ねる: 交点と南北の曖昧性

等遅延線(同心円)と等ドップラー線(平行直線)を天体の見かけの円盤の上に重ねると、格子ができます。ある遅延-ドップラービン $(\Delta\tau, \Delta f)$ に落ちるエコー電力は、対応する円と直線の交点にある表面素片からの散乱です。ここで幾何学上の重要な問題に気づきます。円と直線は一般に2点で交わる——投影自転軸に対して北側と南側の2つの表面素片が、まったく同じ遅延・同じドップラーを持つのです。これが遅延-ドップラーマッピング固有の 南北曖昧性(north-south ambiguity) です。

曖昧性の解決には、(1) アンテナビームをわずかに絞って片半球に感度を寄せる(アレシボの月観測)、(2) 天体の見かけの自転軸の向きが観測期間中に変わることを利用し、複数の観測ジオメトリからの整合性で解く、(3) 形状モデルを仮定してマップ系列全体を最小二乗フィットする、といった方法が使われます。現代の小惑星形状推定は(3)の発展形で、数日間にわたる遅延-ドップラー画像の系列から、自転状態と3次元多面体形状モデルを同時推定する逆問題として解かれています。

一点、注意すべき対比があります。光学画像の1ピクセルが表面の1点に対応するのに対し、遅延-ドップラーマップの1ピクセルは「遅延とドップラーが等しい点の集合」であり、空間座標への写像は自転状態を知って初めて確定する間接的な画像です。それでもこの手法の分解能は距離に依存しません。 $10^7$ km先でも数十mの分解能が出るのは、分解能を決めるのが帯域幅と積分時間という「手元で制御できる量」だからです。角度分解能が距離とともに劣化する光学観測との、最も本質的な違いがここにあります。

実務での使われ方

金星: レーダー天文学の原点

1961年、JPLゴールドストーン(および米英の複数グループ)が金星のレーダーエコー検出に初めて成功し、往復光行時間から天文単位(AU)の値をそれまでの光学観測より3桁近く高い精度(誤差 $\sim10^{-6}$ )で決定しました。これは科学的成果であると同時に、マリナー2号以降のすべての惑星間航法の前提となる実務的成果でした。惑星の位置が1000kmずれていれば探査機は届きません。惑星レーダーは深宇宙航法の測位基盤を整備したのです。

続く1962〜64年、エコーのドップラー広がり $B_{\text{echo}} = 8\pi R_a\cos\delta/(\lambda P)$ の測定から、厚い雲に閉ざされた金星の自転が周期約243日の逆行回転であることが発見されました。表面が一切見えない天体の自転を、スペクトルの幅だけから決めたわけです。1970〜80年代にはアレシボとゴールドストーンの遅延-ドップラーマッピングにより金星表面の地形図が作られ、マクスウェル山地などの主要地形はマゼラン探査機(SARによる全球マッピング、1990年代)に先立って地上からのレーダーで発見されています。

地球近傍小惑星: 形状・軌道決定と惑星防衛

現代の惑星レーダーの最重要ターゲットは地球近傍小惑星(NEA)です。役割は大きく2つあります。

軌道の精密決定。 光学観測が与えるのは天球上の角度2成分だけで、視線方向の情報を直接は持ちません。レーダーはまさにその視線方向の距離(精度数十m〜数百m)と視線速度(ドップラー、精度mm/sオーダー)を与えるため、光学のみの軌道解に1回のレーダー観測を加えるだけで、将来位置の予測誤差が桁で縮みます。実例として、小惑星アポフィス(99942 Apophis)は2021年3月の接近時にゴールドストーンとグリーンバンクのバイスタティック観測で軌道が精密化され、懸念されていた2068年の衝突可能性が完全に否定されました。太陽輻射の非等方再放射による微小な軌道変化(ヤルコフスキー効果)が直接検出できるのも、レーダー測距の精度あってのことです。これらはNASAの惑星防衛(planetary defense)プログラムの根幹をなす観測です。

形状・自転・表面性状の推定。 遅延-ドップラー画像系列からの形状インバージョンにより、トータチス(4179 Toutatis)の非主軸回転、1999 KW4に代表される衛星を持つ連星小惑星、ベンヌ(101955 Bennu)のコマ型形状などが明らかにされてきました。ベンヌの地上レーダー形状モデルは、OSIRIS-RExミッションの計画立案に実際に使われ、探査機到着後のその場観測と良く一致することが確認されています。また偏波比(送信円偏波と同回り/逆回りの受信電力比)からは表面の波長スケールの粗さが、反射率からは表面物質の密度・金属量が推定できます。

設備の現在: アレシボ喪失後

アレシボ305m(S帯2380MHz、900kW)は半世紀にわたり最強の惑星レーダーでしたが、2020年12月の支持ケーブル破断による受信プラットフォーム崩落で失われました。現在、定常運用されている大型惑星レーダーは実質的にゴールドストーンDSS-14(70m、X帯8560MHz、450kW)のソーラーシステムレーダーのみで、グリーンバンク100m望遠鏡などを受信側とするバイスタティック観測で補完されています。DSNの本務である探査機通信・測距とアンテナ時間を分け合っており、「通信インフラと科学観測が同じ設備に同居する」というDSNの二面性を象徴する存在です。

月レーザー測距との対比

「地上から天体までの距離を電磁波の往復で測る」という意味で、惑星レーダーとしばしば対比されるのが月レーザー測距(LLR: Lunar Laser Ranging)です。アポロ11・14・15号とソ連のルノホートが月面に置いたコーナーキューブ再帰反射器にレーザーパルスを撃ち、往復時間から月までの距離をmm精度で測る技術です。両者の違いは本質的です。

LLRのターゲットは光を来た方向に返す協力ターゲットであり、リンクは $1/R^4$ ながら反射器の利得で大きく救済される。惑星レーダーの相手は無指向に散乱する自然表面(非協力ターゲット)。
LLRは反射器という「1点」までの距離を測る測地・重力物理の技術(月の秤動、一般相対論の検証)。惑星レーダーは表面全体のエコーから面としての情報(形状・自転・粗さ)を引き出すイメージング・キャラクタリゼーションの技術。
そして反射器は設置に着陸ミッションを要するのに対し、レーダーは行ったことのない天体を今夜観測できる。

つまり両者は競合ではなく、「協力ターゲット×精密測距」と「非協力ターゲット×リモートセンシング」という直交した役割分担にあります。

演習問題

レーダー方程式 $P_r = P_t G_t G_r \lambda^2 \sigma / ((4\pi)^3 R^4)$ を使って、次の受信電力を概算してください。設備: $P_t = 450$ kW、送受信とも利得 $G = 73$ dBi(X帯70mアンテナ)、 $\lambda = 3.5$ cm。ターゲット: レーダー断面積 $\sigma = 10^5\ \text{m}^2$ (直径約1km・アルベド0.13相当)の小惑星、距離 $R = 0.05$ AU( $7.5\times10^9$ m)。得られた $P_r$ をワットで表し、通信リンクの典型的な受信電力( $10^{-18}$ W オーダー)と比較してください。
問1と同じ小惑星が距離 $0.5$ AU にあるとき、受信電力は何分の1になりますか。また、同じ設備で土星の衛星タイタン( $R \approx 8.5$ AU、 $\sigma \approx 2\times10^{11}\ \text{m}^2$ とする)を観測した場合の受信電力を問1の小惑星(0.05 AU)と比較し、「断面積で6桁勝っても距離4乗で負ける」ことを数値で確認してください。
自転周期 $P = 4.5$ 時間、半径 $R_a = 300$ m の小惑星を、X帯( $\lambda = 3.5$ cm)、 $\cos\delta = 1$ で観測します。(a) エコーの全ドップラー幅 $B_{\text{echo}} = 8\pi R_a\cos\delta/(\lambda P)$ を求めてください。(b) この天体を50本のドップラーストリップに分解するために必要な周波数分解能と、それを実現するコヒーレント積分時間 $T_{\text{coh}} = 1/\Delta f$ を求めてください。
遅延-ドップラーマップの南北曖昧性について、等遅延線(同心円)と等ドップラー線(平行直線)の幾何学に基づいて、(a) なぜ曖昧性が生じるのか、(b) 円盤上で曖昧性が生じない点はどこか(交点が1つになる場所を考える)、を説明してください。

まとめと次回予告

惑星レーダー天文学は、深宇宙通信・測距と同じ地上局インフラと信号処理を使いながら、トランスポンダのない自然天体を相手にする観測技術でした。受信電力はレーダー方程式に従って距離の4乗で減衰し、これが「主戦場は地球近傍小惑星、遠方惑星は円盤積分がやっと」という観測対象の選択を決めていること、そして角度分解能では点にしか見えない天体でも、等遅延線(同心円)と等ドップラー線(自転軸に平行な直線)の格子を使った遅延-ドップラーマッピングによって、距離に依存しない数十m級の分解能で表面を再構成できることを見ました。金星の自転周期の発見からAUの精密決定、そしてアポフィスの衝突可能性の棄却まで、この技術は科学と実務(航法・惑星防衛)の両方を支え続けています。

次回は、電波で天体を「調べる」もうひとつの形態である 掩蔽電波科学(radio occultation) に進みます。今回は地上から電波を撃ちましたが、次回は探査機が惑星の縁に隠れる(掩蔽される)瞬間に、探査機の下りリンク電波が惑星大気をかすめて屈折・減衰するさまを利用して、大気の温度・圧力・電子密度プロファイルを測るという、通信リンクそのものを科学計測器に変えてしまう技術です。

参考文献

S. J. Ostro, “Planetary Radar Astronomy,” Reviews of Modern Physics, vol. 65, no. 4, 1993
S. J. Ostro et al., “Radar Observations of Asteroids,” in Asteroids III, University of Arizona Press
A. E. E. Rogers, T. Hagfors 他の遅延-ドップラー理論初期論文群(Journal of Geophysical Research)
M. I. Skolnik, Introduction to Radar Systems, 3rd ed., McGraw-Hill(レーダー方程式の章)
L. A. M. Benner et al., “Radar Observations of Near-Earth and Main-Belt Asteroids,” in Asteroids IV, University of Arizona Press
I. I. Shapiro, “Planetary Radar Astronomy,” IEEE Spectrum(AU決定・金星自転の歴史的経緯)
DSN Telecommunications Link Design Handbook, DSN No. 810-005(Goldstone Solar System Radar)
J. O. Dickey et al., “Lunar Laser Ranging: A Continuing Legacy of the Apollo Program,” Science