パイオニア・アノマリー — 精密ドップラー追跡が生んだ謎と、系統誤差解析による決着

この回で学ぶこと

トランスポンダの回で、コヒーレントターンアラウンドによる2-wayドップラー計測が、1分程度の積分でμm/sオーダーという驚異的な速度計測精度を達成することを見ました。またクラメール・ラオ下限の回では、その精度の理論的な床が観測時間とSNRでどう決まるかを学びました。

この回では、その「精密すぎる速度計」が歴史上引き起こした最も有名な事件——パイオニア・アノマリー (Pioneer anomaly) ——を題材にします。1972年・1973年に打ち上げられたパイオニア10号・11号のドップラー追跡データは、太陽系の力学モデルで説明できない、太陽方向を向いた約 $8.7\times10^{-10}\ \mathrm{m/s^2}$ の一定の微小減速を示し続けました。これは地表の重力加速度の100億分の1というごく僅かな量ですが、コヒーレントドップラーの精度はそれを有意に検出してしまったのです。「重力の法則が太陽系スケールで破れているのではないか」という新物理の候補として約30年にわたり議論されたこの謎は、最終的に2012年、**探査機自身の熱放射の非対称性(熱リコイル力)**という徹底した系統誤差解析によって決着しました。

この回のゴールは4つです。

1. ドップラー残差の定式化を通じて、アノマリーが観測量(周波数ドリフト)としてどう現れたかを数式で理解する。
2. CRLBの観点からこの検出を定量的に位置づけ、「統計誤差は問題にならないほど小さく、勝負は系統誤差の切り分けにある」ことを数値で確認する。
3. 新物理を主張する前に既知の力(太陽輻射圧・熱リコイル・地球側モデル誤差)を除去する系統誤差切り分けの方法論——それぞれの誤差源が持つ距離依存性・時間依存性・方向という「指紋」で判別する考え方——を学ぶ。
4. Turyshevらの熱モデリングによる解明の結末と、そこから後続ミッションが得た設計上・運用上の教訓を知る。

これは単なる歴史の紹介ではありません。「計測精度が上がるほど、誤差予算は統計誤差ではなく系統誤差に支配される」という状況は、重力科学・精密軌道決定・基礎物理検証を行うすべての現代ミッションに共通する構図であり、パイオニア・アノマリーはその方法論の最良の教科書だからです。

直感的な導入 — 精密すぎる速度計が拾った「余分な減速」

パイオニア10号・11号は、木星・土星をフライバイした後、太陽系を離脱する双曲線軌道に乗りました。この段階の探査機に働く力は、原理的には非常に単純です。太陽と惑星の重力、そしてわずかな太陽輻射圧。これらをすべてモデル化して軌道を積分すれば、地上局で観測されるはずのドップラー周波数が予測できます。

ところが、JPLの軌道決定プログラム(ODP)で実データをフィットすると、観測されたドップラー周波数が予測よりわずかに、しかし一貫してずれ続けていることが分かりました。ずれの向きは「探査機が予測より少しずつ減速している」方向、大きさは片道換算で約 $6\times10^{-9}$ Hz/s——Sバンドダウンリンク約2.3 GHzに対して、8年間積み重なってようやく約1.5 Hzという、気の遠くなるような小さなドリフトです。

ここで2つの事実が効いてきます。第一に、パイオニアはスピン安定方式(約4.8 rpmで機体全体が回転して姿勢を保つ)の探査機で、3軸制御のボイジャーのように姿勢制御スラスタを頻繁に噴かないため、力学モデルを汚す人為的な加速度がほとんどなく、「重力だけで飛ぶ理想的な試験質量」に近い存在でした。第二に、トランスポンダの回で学んだ通り、2-wayコヒーレントドップラーは地上の水素メーザー基準の安定度をそのまま使えるため、この微小なドリフトを雑音に埋もれさせずに追い続けることができました。

つまりパイオニア・アノマリーは、「探査機が理想的に静かで、速度計が異常に正確だったからこそ見えてしまった謎」なのです。精度が2桁悪ければ、誰も気づかないままだったでしょう。

数式で見るアノマリー

ドップラー残差と未知加速度

軌道決定では、観測されたダウンリンク周波数 $f_{\mathrm{obs}}(t)$ と、力学モデル・伝搬モデルから計算した予測周波数 $f_{\mathrm{model}}(t)$ の差、ドップラー残差 (residual)

$$\delta f(t) \equiv f_{\mathrm{obs}}(t) - f_{\mathrm{model}}(t)$$

を監視します。モデルが完全なら、残差は平均ゼロの雑音だけになるはずです。ところがパイオニアの残差には、雑音に加えて時間にほぼ比例して成長する成分が残りました。

視線方向に未知の一定加速度 $a_P$(太陽方向、すなわち減速の向き)が働いていると仮定すると、視線速度の誤差は $\delta v_r(t) = a_P\, t$ と線形に積み上がります。2-wayドップラーの公式 $\Delta f = (2 v_r/c) f_c$ にこれを代入すれば、残差の周波数ドリフト率は

$$\boxed{\;\dot{(\delta f)} = \frac{2 a_P}{c}\, f_c\;}$$

となります(片道換算ならば係数2が外れます)。Anderson らの解析(1998年、2002年)が報告した値は

$$a_P = (8.74 \pm 1.33)\times 10^{-10}\ \mathrm{m/s^2} \quad (\text{太陽方向})$$

であり、観測された片道換算のドリフト率 $(5.99\pm0.01)\times10^{-9}$ Hz/s とオーダーとして整合します(厳密な換算には軌道フィット上の相関や補正が絡むため、単純な一括変換とは数%〜1割程度ずれます。演習1で確認します)。重要なのは、この加速度がパイオニア10号・11号の両方に、太陽距離約20 AUから70 AUまで、ほぼ一定の大きさで現れたことです。

積算効果の大きさ — なぜ「微小」なのに見えるのか

$a_P \sim 8.7\times10^{-10}\ \mathrm{m/s^2}$ という数字がどれほど小さいか、そしてなぜそれでも見えるのかを数値で掴みましょう。1年($t \approx 3.16\times10^{7}$ s)積算すると、

$$\Delta v = a_P\, t \approx 8.74\times10^{-10} \times 3.16\times10^{7} \approx 2.8\times10^{-2}\ \mathrm{m/s} \approx 2.8\ \mathrm{cm/s}$$ $$\Delta \rho = \frac{1}{2} a_P\, t^2 \approx \frac{1}{2}\times 8.74\times10^{-10}\times (3.16\times10^{7})^2 \approx 4.4\times10^{5}\ \mathrm{m} \approx 440\ \mathrm{km}$$

1年で速度に2.8 cm/s、距離に440 kmもの食い違いが積み上がります。2-wayドップラーの単発精度がμm/s〜数十μm/sオーダーであったことを思い出すと、年間2.8 cm/s(= 28,000 μm/s)のずれは統計誤差の数千倍であり、検出自体は疑いようがありません。微小加速度の検出において観測アーク(追跡期間)の長さが決定的に効く——$\Delta v$ は $t$ に、$\Delta\rho$ は $t^2$ に比例して成長する——ことが、この計算から読み取れます。

統計限界の中での位置づけ — CRLBとの接続

では、この検出の統計的な限界はどこにあるのでしょうか。CRLBの回で学んだ枠組みをそのまま応用します。一定加速度の推定は、速度時系列 $v_i = v_0 + a t_i + n_i$($n_i$ は分散 $\sigma_v^2$ の白色雑音)に対する直線の傾きの推定とみなせます。等間隔な $N$ 点の観測を全長 $T$ にわたって取ったときの最小二乗傾き推定の分散は、よく知られた公式

$$\sigma_a^2 = \frac{\sigma_v^2}{\sum_i (t_i - \bar{t})^2} \approx \frac{12\,\sigma_v^2}{N\,T^2}$$

で与えられます(等間隔標本の時刻の2次モーメント $\sum(t_i-\bar t)^2 \approx N T^2/12$ を使いました。ガウス雑音に対してはこれがそのままCRLBに一致します)。現実的な数字として $\sigma_v = 10\ \mu\mathrm{m/s}$、$N = 10^4$ 回の測定、$T = 20$ 年 $\approx 6.3\times10^{8}$ s を入れると、

$$\sigma_a \approx \sqrt{\frac{12}{10^4}}\times\frac{10^{-5}}{6.3\times10^{8}} \approx 5.5\times10^{-16}\ \mathrm{m/s^2}$$

統計的な検出限界は $10^{-16}\ \mathrm{m/s^2}$ 台——アノマリーの値 $8.74\times10^{-10}$ より6桁も下です。しかし実際に報告された不確かさは $\pm1.33\times10^{-10}$、統計限界より5桁以上大きい。この巨大なギャップこそが、この回の核心です。

長期のドップラー追跡において、誤差予算は白色雑音の統計誤差ではなく、系統誤差に完全に支配されます。 実際のドップラー雑音は長時間では白色でなく(太陽プラズマ・対流圏・基準発振器のゆらぎによる有色雑音でアラン偏差に床ができる)、さらに「モデル化されていない力」はどれも残差に相関した構造を刻みます。CRLBは「これ以上は原理的に無理」という床を与えますが、パイオニア・アノマリーの争点はその床のはるか上、「残差に見えているこの構造は、新物理か、それともモデル化し忘れた既知の物理か」というモデルの同定問題だったのです。

系統誤差の切り分け — 「新物理」を主張する前に

誤差源の「指紋」で判別する

未知の加速度が見つかったとき、軌道決定の専門家が最初にすることは、新理論の論文を書くことではなく、既知の力とモデル誤差のリストを一つずつ潰していくことです。このとき強力な武器になるのが、各誤差源が固有に持つ距離依存性・時間依存性・方向という指紋です。候補となる力を一般形

$$a_{\mathrm{model}}(t, r) = a_{\odot}\!\left(\frac{r_0}{r}\right)^{2} + a_{\mathrm{th}}\, e^{-\lambda t} + a_{\mathrm{const}} + (\text{日周・年周項})$$

のように分解してフィットし、どの項が残差を説明するかを競わせるのです。観測されたアノマリーの指紋は「距離20–70 AUでほぼ一定、方向は太陽(地球)向き、両機に共通」でした。以下、主要な容疑者を順に取り調べます。

容疑者1: 太陽輻射圧 — 距離依存性が合わない

太陽光が機体(主に直径2.74 mの高利得アンテナ皿、面積 $A \approx 5.9\ \mathrm{m^2}$)に当たって押す力による加速度は、太陽光度 $L_\odot = 3.846\times10^{26}$ W、機体質量 $m \approx 241$ kg、反射特性をまとめた係数 $\kappa \sim 1.7$ を使って

$$a_{\mathrm{srp}}(r) = \frac{\kappa\, L_\odot A}{4\pi r^2 c\, m}$$

と書けます。数値を入れると、$r = 20$ AU($3.0\times10^{12}$ m)で $a_{\mathrm{srp}} \approx 5\times10^{-10}\ \mathrm{m/s^2}$ とアノマリーに匹敵しますが、40 AUでは $1/4$ の約 $1.2\times10^{-10}$ まで落ちます。輻射圧は $1/r^2$ で減衰するのに対し、アノマリーは20 AUから70 AUまで一定——指紋が合いません。このため解析では、輻射圧が十分小さくなる20 AU以遠のデータを主に使い、$1/r^2$ 項は別パラメータとして推定して分離しました。太陽輻射圧は「20 AU以遠では主犯たりえない」として除外されます。

容疑者2: 地球側・伝搬経路のモデル誤差 — 周期性が合わない

地球の重力場や自転・歳差モデルの誤差、地上局位置の誤差、地球暦(エフェメリス)の誤差は、観測ジオメトリが地球の自転・公転とともに変わるため、残差に日周・年周の周期的な指紋を残します。実際、パイオニアの残差にも小さな年周項・日周項が見つかりましたが、その振幅から逆算される加速度への寄与は $a_P$ よりはるかに小さく、しかも一定項とは直交する(フィット上分離できる)ことが確認されました。太陽コロナや対流圏の伝搬遅延も、太陽離角や仰角に相関する固有の指紋を持つため、較正モデルとの差として評価され、いずれも $10^{-10}\ \mathrm{m/s^2}$ の一定項を作れないことが示されました。さらに、JPLのODPとは独立に開発されたエアロスペース社のソフトウェアCHASMPでも同じアノマリーが再現され、「解析ソフトのバグ」説も棄却されました。

容疑者3: 探査機自身の熱 — 唯一、指紋が合う

最後に残ったのが、探査機自身が放つ熱です。光子は運動量 $p = E/c$ を運ぶため、熱放射が機体から非対称に放出されれば、その反動(熱リコイル力)が機体を押します。完全にコリメートされた放射パワー $P$ が作る力は

$$F = \frac{P}{c}$$

であり、アノマリーを説明するのに必要な力 $F = m\, a_P \approx 241 \times 8.74\times10^{-10} \approx 2.1\times10^{-7}$ N を作るには、正味

$$P = F c \approx 2.1\times10^{-7} \times 3\times10^{8} \approx 63\ \mathrm{W}$$

の「一方向に偏った」放射があれば足ります(ランバート面なら $F = 2P/(3c)$ なので約95 W)。ここで機体の熱収支を見ると、パイオニアは4基のSNAP-19型RTG(放射性同位体熱電発電機)を搭載し、打ち上げ時の熱出力は合計約2.5 kW、機器区画の電気系廃熱も約100 W ありました。つまり必要な63 Wは、機体が捨てている全熱のわずか数%が一方向に偏れば説明できてしまう規模なのです。

しかも指紋が合います。パイオニアの高利得アンテナは常に地球(ほぼ太陽方向)を向くため、機体に固定された非対称性は常に太陽方向の減速として現れます。距離には依存しません。唯一の弱点は時間依存性で、RTGの熱源プルトニウム238は半減期87.7年で減衰し、電源電力の低下とともに廃熱の配分も変わるため、熱リコイルは数十年スケールでゆっくり減衰するはずです。初期の解析(約11年間のデータ)では一定加速度と減衰加速度を区別できず、ここが最後の争点として残りました。

解明 — データ発掘と熱モデリングによる決着(2012年)

決着をつけたのは、新しい観測ではなく古いデータの発掘と、探査機の精密な熱モデルでした。Turyshev と Toth らは、惑星協会の支援も受けて、廃棄寸前だった磁気テープ等からパイオニア10号の約23年分・11号の約11年分に及ぶドップラー追跡データ(ATDF形式のアーカイブ)と、機体各部の温度・電力を含むテレメトリの全記録を復元しました。

このテレメトリで各時点の熱源(RTG熱出力、機器の消費電力)と温度分布を拘束しながら、機体の3次元有限要素熱モデルを構築し、表面各要素からの放射が運ぶ運動量を積分して熱リコイル加速度を計算しました。結論は明快でした。

• 熱リコイル加速度の計算値はアノマリーの観測値と誤差の範囲で一致する(独立に解析した Rievers & Lämmerzahl (2011) の有限要素モデルも約 $7.4\times10^{-10}\ \mathrm{m/s^2}$ という整合的な値を得ました)。
• 寄与の内訳は、機器区画の電気系廃熱の非対称放出と、RTGの熱がアンテナ皿の背面で反射されて前方(地球方向)へ偏って逃げる成分が主で、いずれも機体を太陽方向へ押し戻す向き。
• 拡張されたデータセットを解析し直すと、アノマリーは厳密には一定ではなく時間とともに緩やかに減衰しており、これは放射性崩壊と電源低下による熱源の減衰と整合する。

こうして2012年、Turyshev らは「パイオニア・アノマリーの熱起源を支持する (Support for the thermal origin of the Pioneer anomaly)」と題した論文で、30年来の謎に事実上の終止符を打ちました。新物理は不要で、犯人は探査機自身の熱だったのです。重力法則の修正(MOND等)でアノマリーを説明しようとした理論群は、この結論と、同種の異常が惑星の軌道(はるかに高精度で追跡されている)に見られないことから支持を失いました。

一連の顛末は「がっかりな結末」ではありません。$10^{-10}\ \mathrm{m/s^2}$ という加速度スケールでは探査機の設計図・運用記録・温度テレメトリまでもが「観測装置の一部」になる、ということを実証した点で、精密航法と宇宙基礎物理の方法論を一段引き上げた成果と評価されています。

実務での使われ方

パイオニア・アノマリーの教訓は、現代の深宇宙ミッションの設計・運用に直接組み込まれています。

• 非重力加速度の標準モデル化。現在のJPLやESAの軌道決定では、熱リコイル力は太陽輻射圧と並ぶ標準的なモデル項です。RTGを搭載するニュー・ホライズンズやカッシーニの航法では、RTG熱と機体形状から熱リコイル加速度($10^{-10}\ \mathrm{m/s^2}$ オーダー)をあらかじめ計算して力学モデルに入れており、パイオニアで起きた「モデル化し忘れ」は再発しない体制になっています。
• 熱設計への反映。精密な重力科学を目的とする探査機では、熱放射をできるだけ対称に逃がす(ラジエータ配置の対称化)、あるいは非対称性を地上で精密に較正できるようにする(表面特性・温度センサ配置の文書化)ことが設計要求に入るようになりました。BepiColomboやJUICEの重力科学実験の誤差予算にも、非重力加速度のモデル化誤差が明示的な項として計上されています。
• 追跡データの長期アーカイブ。パイオニアの解明が可能だったのは、1970〜90年代のドップラーデータとテレメトリが(かろうじて)残っていたからです。この経験を踏まえ、DSNの追跡データ(現在はTRK-2-34等の標準形式)やCCSDS勧告に基づくテレメトリアーカイブを長期保存し、将来の再解析に備えることの科学的価値が広く認識されるようになりました。実際、過去の追跡データの再解析は、地球フライバイ・アノマリーの検証や暦の精密化など、今も現役の研究手法です。
• 「探査機を試験質量として使う」研究の方法論。太陽系スケールで重力理論を検証する後続研究(カッシーニのPPNパラメータ測定など)では、パイオニアで確立された「統計誤差ではなく系統誤差の予算表を作り、指紋で切り分ける」という方法論がそのまま踏襲されています。

演習問題

1. 観測された片道換算の周波数ドリフト率 $\dot{(\delta f)} = 5.99\times10^{-9}$ Hz/s とSバンドダウンリンク周波数 $f_c = 2.295$ GHz から、片道の関係式 $\dot{(\delta f)} = (a_P/c) f_c$ を逆に解いて加速度を求め ($a_P = c\,\dot{(\delta f)}/f_c$)、報告値 $8.74\times10^{-10}\ \mathrm{m/s^2}$ と比較してください。オーダーが一致すること、および1割程度の差が残ること(軌道フィット上の補正・相関によるもの)を確認してください。

2. パイオニア10号の質量を $m = 241$ kg として、アノマリーを説明するのに必要な力 $F = m\, a_P$ を求め、(a) 完全にコリメートされた放射 $P = Fc$、(b) ランバート放射 $P = 3Fc/2$ のそれぞれの場合に必要な正味の指向性放射パワーを計算してください。RTG総熱出力約2.5 kWに対する割合(%)も求め、「全廃熱のわずかな非対称で足りてしまう」ことを確認してください。

3. 本文の傾き推定の式 $\sigma_a \approx \sigma_v\sqrt{12/N}/T$ を使い、$\sigma_v = 10\ \mu\mathrm{m/s}$、$N = 10^4$、観測アーク $T = 5$ 年の場合の統計限界 $\sigma_a$ を計算してください。その値を報告された不確かさ $1.33\times10^{-10}\ \mathrm{m/s^2}$ と比較し、誤差予算がなぜ統計誤差ではなく系統誤差に支配されるのか、本文の議論に沿って説明してください。

4. 太陽輻射圧の式 $a_{\mathrm{srp}}(r) = \kappa L_\odot A/(4\pi r^2 c\, m)$($\kappa = 1.7$, $A = 5.9\ \mathrm{m^2}$, $m = 241$ kg)を使って、$r = 10, 20, 40$ AU での $a_{\mathrm{srp}}$ を計算してください(1 AU $= 1.5\times10^{11}$ m)。その結果をもとに、「アノマリーが20–70 AUでほぼ一定だった」という観測事実が太陽輻射圧犯人説をどう棄却するか、また熱リコイル犯人説とはなぜ矛盾しないかを、各誤差源の「指紋」という観点から論じてください。

まとめと次回予告

この回では、コヒーレント2-wayドップラーの超高精度が生んだ歴史的な謎、パイオニア・アノマリーを追いました。ドップラー残差の一定ドリフトとして現れた太陽方向の微小減速 $a_P \approx 8.7\times10^{-10}\ \mathrm{m/s^2}$ は、統計限界(CRLB)からは6桁も余裕のある「楽な」検出である一方、その解釈は系統誤差の同定問題そのものでした。距離依存性・時間依存性・方向という指紋を使って太陽輻射圧・地球側モデル誤差を消去し、最後に残った熱リコイル力を、発掘された数十年分の追跡データとテレメトリに基づく精密熱モデルで定量したことで、2012年に「探査機自身の熱の非対称放出」として決着——この一連の流れは、精密計測時代のミッションにおける誤差解析の模範例です。あわせて、追跡データを長期アーカイブしておくことの科学的価値と、熱放射の対称化・非重力加速度の標準モデル化という後続機への教訓も確認しました。

さて、今回の主役だった2-wayドップラーは地上の水素メーザー基準に支えられていましたが、探査機が地上の助けなしに自律航法を行う場面では、搭載時計そのものの長期安定性が精度の天井になります。次回は、探査機搭載時計(USO・小型原子時計)の周波数ドリフトを長期にわたってどうモデル化し、どう較正するかという「時計のエイジング」の話に進みます。今回学んだ「ゆっくり成長する系統誤差をモデルで捉える」という視点が、そのまま時計のドリフト推定にも活きてきます。

参考文献

• J. D. Anderson, P. A. Laing, E. L. Lau, A. S. Liu, M. M. Nieto, S. G. Turyshev, “Indication, from Pioneer 10/11, Galileo, and Ulysses Data, of an Apparent Anomalous, Weak, Long-Range Acceleration,” Physical Review Letters 81, 2858 (1998)
• J. D. Anderson et al., “Study of the anomalous acceleration of Pioneer 10 and 11,” Physical Review D 65, 082004 (2002)
• S. G. Turyshev, V. T. Toth, G. Kinsella, S.-C. Lee, S. M. Lok, J. Ellis, “Support for the Thermal Origin of the Pioneer Anomaly,” Physical Review Letters 108, 241101 (2012)
• B. Rievers, C. Lämmerzahl, “High precision thermal modeling of complex systems with application to the flyby and Pioneer anomaly,” Annalen der Physik 523, 439 (2011)
• S. G. Turyshev, V. T. Toth, “The Pioneer Anomaly,” Living Reviews in Relativity 13, 4 (2010)
• DSN Telecommunications Link Design Handbook, DSN No. 810-005(Doppler Tracking に関するモジュール)
• J. H. Yuen (ed.), Deep Space Telecommunications Systems Engineering, JPL Publication 82-76