Filling the gap in the IERS C01 polar motion series in 1858.9-1860.9

本論文は、IERS C01 極運動系列における 1858.9-1860.9 年の 2 年間の欠損データを、パラメトリック天文学モデルとデータ駆動型の特異スペクトル解析（SSA）の 2 つの手法を用いて埋め、特に完全なモデルに基づく SSA 手法の結果が好ましいことを示しています。

要約

この論文は、**「地球の自転の記録に、2 年間も『空白のページ』が空いていたので、それを埋めてつなげました」**というお話です。

専門用語を避け、誰でもわかるような例え話を使って説明しますね。

1. 物語の舞台：「地球の揺れ」と「記録帳」

まず、地球はただ回っているだけではありません。氷河が溶けたり、大気が動いたりして、**「自転軸が少し揺れている（極運動）」**のです。これを「ポールの揺れ（Polar Motion）」と呼びます。

科学者たちは、1846 年（明治維新の少し前）から現在まで、この揺れを毎日記録してきました。この記録帳は**「IERS C01」**という世界標準のデータで、非常に貴重です。

2. 問題発生：「2 年間の欠落」

しかし、この記録帳には1858 年から 1860 年までの 2 年間、あるページが破れてなくなっていました。

• なぜ？
  当時、アメリカのワシントン天文台で観測が止まってしまったからです。イギリス（グリニッジ）とロシア（プルコヴォ）の観測所は残っていましたが、それらだけでは「南北方向の揺れ（Yp）」を正確に測るのに必要な角度が足りなかったのです。
• なぜ困る？
  科学者がこのデータを分析する際、**「データが途切れていると、計算が狂ったり、間違ったリズム（周波数）が見えてしまったりする」**からです。まるで、音楽の CD が 2 曲分飛んでいたら、曲の全体像がわからなくなるのと同じです。

3. 解決策：「2 つの天才的な推測」

そこで、この論文の著者たちは、**「失われた 2 年間のデータを、数学的な魔法で復元しよう」**と試みました。2 つの異なるアプローチ（方法）を使いました。

方法 A：「おなじみのリズム」を頼る（パラメトリックモデル）

• イメージ： 「昔から続くリズム」
  地球の揺れには、決まったリズムがあります。

  1. チャンドラー振動（CW）： 約 14 ヶ月で回る大きな揺れ。
  2. 年振動（AW）： 1 年周期で回る季節的な揺れ。
  3. バイアス： 全体の傾き。

  これらが「少しずつ大きくなったり小さくなったりしている」と仮定して、**「前後のデータから、このリズムがどうなっていたかを数式で計算し、欠けた部分を埋める」**という方法です。

  • 例え話： 散歩している人が「1 歩、2 歩、3 歩…」と歩いているのが見えて、4 歩目が見えない。でも「歩幅が少し大きくなっている傾向がある」とわかれば、「4 歩目はここかな？」と推測して埋める感じです。

方法 B：「データそのものが語る」アプローチ（SSA 法）

• イメージ： 「パズルの完成」
  これは、事前に「リズムはこれだ！」と決めるのではなく、**「残っているデータ全体をパズルのように分析し、隠れたパターンを自動で見つけ出す」**という方法です（特異スペクトル分析：SSA）。
  • 例え話： 欠けたパズルのピースがない状態で、周りのピースの形や色を細かく分析し、「ここにはきっとこの色のピースが来るはずだ」と、データが教えてくれるままにピースを埋めていきます。
  • この論文では、X 軸と Y 軸のデータを組み合わせて分析する「複素 SSA（CSSA）」という高度なテクニックを使いました。

4. 結果：「2 つの方法はよく一致した！」

計算の結果、両方の方法で埋めたデータは、**「元の記録の誤差の範囲内で、ほぼ同じ結果」**になりました。

• どちらが優れている？
  著者たちは、「方法 B（SSA 法）」の方が少し優れていると考えました。
  • 理由： 方法 A は「リズムはこれだ」という前提が必要ですが、方法 B はデータが持っている「複雑な秘密」まで見つけ出せるからです。まるで、単純な計算式で推測するより、AI が過去の全データを学習して推測する方が、より現実に近い答えが出やすいようなものです。

5. 結論：「歴史の欠片が揃った」

この研究によって、1846 年から現在までの「地球の揺れ」の記録が、2 年間の空白なく、つなげられました。

• なぜ重要？
  これにより、科学者たちは「地球の自転が、100 年、200 年という長いスパンでどう変化してきたか」を、より正確に研究できるようになります。気候変動や地球内部の構造の変化を理解する上で、この「つなげられた記録」は非常に重要な手がかりになります。

まとめ

この論文は、**「19 世紀の観測記録に空いていた 2 年間の『空白』を、最新の数学の力で見事に埋め、地球の歴史をより完全な物語に仕上げた」**という、科学の「修復作業」の成功物語です。

失われたページを、推測と計算の力で取り戻し、未来の科学者たちが安心して地球の謎を解き明かせるようにしたのです。

技術概要

以下は、arXiv:2412.07868v3「Filling the gap in the IERS C01 polar motion series in 1858.9–1860.9」に基づく技術的な要約です。

論文の概要

本論文は、国際地球回転・基準系サービス（IERS）が提供する地球回転パラメータ（EOP）の標準時系列である「IERS C01」において、1858.9 年から 1860.9 年までの約 2 年間にわたって欠落している極運動（Polar Motion: PM）の Y 座標（Yp）データを、初めて補完する試みを行ったものである。この欠損は、19 世紀の観測データ（特にワシントン天文台の欠如）に起因しており、連続的な時系列解析を困難にしている。

1. 問題の背景と課題

• IERS C01 の重要性: 1846 年から現在までの地球回転の最長かつ信頼性の高い記録であり、長期的な極運動の変動（チャンドラー振動や年周振動など）を研究する上で不可欠である。
• 欠損の問題: Yp 座標において 1858.9〜1860.9 の間に 21 個のデータ点（2 年間）が欠落している。
• 影響: 欠損データは統計的解析、特にスペクトル解析において偽の周波数を導入するリスクがあり、連続的で均等間隔のデータが必要な手法を適用する際に支障をきたす。また、欠損期間中の地球回転モードの進化に関する情報が失われている。
• 原因: 19 世紀の極運動計算にはグリーンウィッチ、プルコヴォ、ワシントンという 3 観測所の緯度観測データが用いられていた。しかし、この期間、ワシントン観測所での観測が欠如していた。グリーンウィッチ（経度 0 度）とプルコヴォ（経度 30 度）は Xp 方向に近い経度に位置するため、これら 2 地点のデータのみでは Yp 方向の情報を正確に推定できず、結果として Yp データが欠落した。

2. 手法（Methodology）

本研究では、欠損データを補完するために、パラメトリックな天文学モデルとデータ駆動型のモデルという 2 つのアプローチを提案・比較検証した。

A. パラメトリック天文モデル（BCA モデルおよび BCA2 モデル）

• 基本構成: 偏位（Bias）、チャンドラー振動（CW）、年周振動（AW）の 3 つの成分から構成される。
• 特徴: 従来の単純な正弦波モデルではなく、CW と AW の振幅が時間とともに線形に変化するという仮定を導入した（式 3, 4）。
• BCA2 モデル: 周波数スペクトル分析で CW 帯域に 2 つのピーク（1.181 年と 1.231 年）が観測されたため、CW 成分を 2 つの異なる周期を持つ項として追加した 13 パラメータモデル（式 5）。
• 推定: 最小二乗法（LS）を用いて、欠損期間を中心とした 12〜24 年の区間でパラメータをフィッティングし、欠損値を算出した。

B. データ駆動型モデル（特異スペクトル解析：SSA）

• 手法: 特異スペクトル解析（Singular Spectrum Analysis, SSA）を用いた。パラメトリックモデルのように事前のモデル構造を仮定せず、データから信号成分を適応的に抽出する。
• 応用: 極運動の物理的性質（Xp と Yp が位相がπ/2 ずれた調和振動である）を考慮し、複素数時系列（X + iY）として扱う**複素 SSA（CSSA）**を採用した。また、多変量 SSA（MSSA）や 1 次元 SSA（1D-SSA）とも比較した。
• 最適化: ウィンドウ長（L）と信号成分の数（モデル次数 r）を人工的な欠損データを用いた検証により最適化し、欠損補填の誤差を最小化するパラメータ組合せを探索した。

3. 主要な結果

• モデルの精度評価:
  • パラメトリックモデル（BCA/BCA2）と SSA モデルの両方で、欠損期間の前後のデータを用いた再構成誤差（RMSE, MAE）を計算した。
  • 誤差の比較: 両手法とも、IERS C01 シリーズに付与されている Yp の不確かさ（約 0.09 秒）と同程度の誤差（0.08〜0.09 秒）で欠損を補填できた。
  • SSA の優位性: 誤差値は BCA モデルの方がわずかに小さかったが、SSA モデルの方が IERS C01 のデータ構造をより包括的に捉えており、事前の仮定（線形振幅変化など）に依存しないため、より信頼性が高いと判断された。
• 補填値の一致:
  • 両手法で得られた補填値は、特に欠損期間の前半（1859 年〜1860 年）でよく一致していた。後半（1860 年〜1861 年）では若干の乖離が見られたが、観測誤差の範囲内で許容可能なレベルであった。
  • 最終的に、7 つの最適な CSSA 設定（L=60, r=4,5; L=119, r=6-9; L=238, r=7）の平均値を推奨される補填値として提示した（Table 4）。

4. 貢献と意義

• 初の試み: IERS C01 系列の 19 世紀における 2 年間の Yp 欠損を体系的に補完する最初の研究である。
• 連続時系列の提供: 補填値を IERS C01 シリーズに追加することで、1846 年以降の完全な連続かつ均等間隔の極運動時系列が実現された。これにより、欠損を無視したり、不適切な補間を行ったりすることによる解析バイアスを排除できる。
• 長期的変動の理解: 19 世紀のデータは、チャンドラー振動の振幅変動（約 66 年周期や 80 年周期）や位相変動を研究する上で極めて重要である。この補填により、1900 年以降のデータのみでは推定が困難だった長期的な地球回転ダイナミクスの解明が促進される。
• 手法の比較: パラメトリックモデルと非パラメトリック（データ駆動）モデルの両方を適用し、その特性と限界を明確にした。特に、SSA が複雑な時系列構造を捉える能力が高いことを実証した。

5. 結論

本研究は、IERS C01 系列の欠損データを SSA を中心とした手法で高精度に補填し、地球回転の長期的研究に不可欠な連続時系列を完成させた。補填されたデータは、観測誤差の範囲内で信頼性が高く、科学者たちが 1846 年からの地球回転変動をより詳細に解析するための基盤を提供する。今後は、19 世紀の歴史的観測データの掘り起こしと再処理を通じて、さらに過去（1820 年代など）への時系列拡張が期待される。