Effects of Geomagnetic Cutoff Rigidity Variations during Forbush Decreases

以下は、提示された論文「Effects of Geomagnetic Cutoff Rigidity Variations during Forbush Decreases（フォバッシュ減少中の地磁気カットオフ硬さ変動の影響）」の技術的サマリーです。

1. 研究の背景と問題提起

• フォバッシュ減少（FD）の観測課題: 太陽風擾乱（特に ICME）に伴う銀河宇宙線（GCR）の短期的な減少現象であるフォバッシュ減少（FD）を解明するには、時間分解能の高いプロトン束の測定が不可欠です。地上の中性子モニター（NM）は長期間にわたり高時間分解能（1 時間単位など）で観測を提供していますが、NM の計数率は「地磁気カットオフ硬さ」以上の一次宇宙線に対する積分応答であり、硬さごとのプロトン束を直接得ることはできません。
• 既存手法の限界: 近年、NM 計数率を深層学習モデルを用いて AMS（Alpha Magnetic Spectrometer）のデータで訓練し、時間分解能の高いプロトン束を再構築する手法が提案されています。しかし、強い ICME 事象に伴う地磁気嵐（Geomagnetic Storm）が発生すると、地上の地磁気カットオフ硬さが変動します。
• 核心的な問題: 地磁気カットオフ硬さの変動は、NM 計数率に局所的な異常（過剰または不足）をもたらします。特に、これまでの研究では 10 GV 付近での影響が強調されていましたが、本研究ではこの影響がより低い硬さ（1 GV まで）に及ぶ可能性を示唆しています。これらの地磁気由来のノイズが除去されていない場合、FD の硬さ依存性（スペクトル形状）の推定にバイアスが生じ、物理的な解釈が歪められるリスクがあります。

2. 手法とアプローチ

本研究は、2011 年 5 月から 2019 年 10 月までの期間（AMS の公開された日次プロトン束データが利用可能な期間）を対象に、以下の手法を採用しました。

• データの比較と参照:
  • 対象データ: NM 観測から再構築された「修正前の時間分解能（1 時間）プロトン束」。
  • 参照データ: 宇宙空間で直接硬さを測定する AMS の「日次プロトン束」。AMS は地上の地磁気カットオフの影響を受けないため、安定したスペクトル参照基準として機能します。
• 異常の検出:
  • ICME 事象の期間において、NM 再構築データの 1 日平均値と AMS の日次データを比較します。
  • AMS の不確かさの 3σ 以上を逸脱する硬さビン（rigidity bins）と日付を特定し、これらを「異常区間」として選定します。
  • 選定された連続する区間は、前後 1 日ずつ拡張して処理範囲とします。
• 補正アルゴリズム:
  • AMS の日次スペクトルを、以下のパラメータ化モデルで記述します（式 1）：
    $$\gamma(R) = \gamma_\infty + \frac{A}{R^k}$$
    ここで、$\gamma_\infty$ は高硬さ限界、$A$ は低硬さの偏差、$k$ は硬さに対する減衰率です。
  • 選定された区間において、AMS から得られたスペクトルパラメータ（$\gamma_\infty, A, k$）を日次から時間分解能へ線形補間し、時間的に連続したスペクトル制約を生成します。
  • 補正対象の硬さビンにおけるプロトン束を、残りの正常なビンの NM 再構築データから得られる時間ごとの正規化係数（$C_t$）と、AMS 由来のスペクトル形状を組み合わせて再計算（式 3）します。
  • 重要点: AMS の日次データを直接時間系列に挿入するのではなく、AMS の「スペクトル形状の制約」のみを用いて、NM 由来の時間変動性を保持しつつ局所的な歪みを修正します。

3. 主要な結果

• 低硬さ域への影響の拡大:
  • 従来の研究で 10 GV 付近に限定されていたカットオフ硬さ変動の影響は、実際にはより低い硬さ域に及んでいることが判明しました。具体的には、2015 年 6 月 25 日の ICME 事象において、1 GV 付近（2.15-2.4 GV ビン）から 10 GV 付近（9.26-10.1 GV ビン）まで、NM 再構築データに AMS 基準からの局所的な逸脱（過剰または不足）が観測されました。
• 補正の効果:
  • 2015 年 6 月 25 日の事象を例に、補正前後を比較しました。
    • 低硬さ域（2.15-2.4 GV）: 補正前は AMS 値より低く推定されていましたが、補正により AMS 値に近づきました。
    • 高硬さ域（9.26-10.1 GV）: 補正前は AMS 値より高く推定されていましたが、補正により AMS 値に近づきました。
  • 地磁気嵐のピーク（Dst 指数が最も負の値をとる期間）に、NM 再構築データと AMS 基準との乖離が最大になることが確認されました。
• スペクトルパラメータの安定性:
  • 付録の分析により、AMS 日次データにおいて、FD 発生日と静穏日のスペクトルパラメータ（$\gamma_\infty, A, k$）は、全体的な傾向として同じパラメータ - 束の関係を追従することが示されました。これは、FD 時の異常が「スペクトル形状そのものの変化」ではなく、「地磁気遮蔽の変動に起因する局所的な歪み」であることを裏付けています。

4. 結論と意義

• 結論:
  • 地上 NM 観測に基づく時間分解能の高いプロトン束再構築データには、地磁気嵐に伴うカットオフ硬さ変動に起因する局所的な歪み（1 GV から 10 GV まで）が含まれている。
  • この歪みを除去するために、AMS の日次スペクトルを外部参照として用いた補正手法を提案し、特定の ICME 事象に対して有効であることを実証した。
• 科学的意義:
  • データ品質の向上: 補正されたデータは、短時間スケールでの宇宙線変調の研究において、より一貫性のある基礎データを提供します。
  • 物理的解釈の明確化: 以前は「太陽圏内の変調」として解釈されていた FD のスペクトル変化の一部が、実際には「地磁気遮蔽の変化」に起因する可能性を示し、両者の分離を可能にします。
  • 手法の革新: 現在、時間分解能と硬さカバレッジの両面で AMS と同等の独立した地上データが存在しないため、AMS を「スペクトル形状の制約」として利用し、NM の時間変動性を維持するハイブリッドな補正アプローチは、将来の宇宙線モニタリングにおいて重要な指針となります。

本研究は、地上観測データの限界を認識しつつ、宇宙空間観測データを巧みに活用することで、宇宙線物理学における短時間スケールの現象解明の精度を高めることを目指しています。