Optimization Algorithm for Determining the Source Surface Radius Based on Parker Solar Probe in situ Measurements from Encounters 1 to 19

技術要約：パーカー・ソーラー・プローブのin situ観測に基づくソース表面半径決定のための最適化アルゴリズム

問題提起
ポテンシャル磁場ソース表面（PFSS）外挿法は、光球磁力計から大規模なコロナ磁場を推定するための標準的な手法である。このモデルにおける重要なパラメータは、磁場が放射状になると仮定される球状境界であるソース表面半径（$R_{ss}$）である。$R_{ss} = 2.5 R_s$ が慣習的な基準値であるが、先行研究では、単一の値ではすべての観測的制約（コロナホール、開いた磁束、磁極性など）を普遍的に満たすことはできないことが示されている。最適な $R_{ss}$ は太陽活動に伴って変化することが知られているが、厳密かつ観測に基づいた値を決定することは依然として困難である。本論文では、エンカウンター1–19におけるパーカー・ソーラー・プローブ（PSP）の in situ 磁場測定値とPFSS外挿との間の不一致を最小化することにより、最適な $R_{ss}$ を決定するという逆問題を扱う。

手法
著者らは、$R_{ss}$ の決定を自由境界フレームワーク内における形状最適化問題として捉えている。

1. 数学的定式化:

   • 前向きのPFSS問題は、混合ディリクレおよびノイマン条件を持つ二階楕円境界値問題として定式化される。著者らは、ソボレフ空間 $H^1(\Omega)$ における前向き解のウェルポーズド性（存在性と一意性）を証明している。
   • 逆問題は、許容可能な領域集合上での目的関数の最小化として構成される。著者らは、許容集合のコンパクト性と目的関数の連続性を証明し、それによって逆問題における大域的最小解の存在を確立した。これらの証明は、球状および非球状のソース表面幾何学の両方に適用される。

2. データと前処理:

   • 入力: グローバル・オシレーション・ネットワーク・グループ（GONG）による日向光球磁力計が内部境界条件として機能する。
   • 観測: エンカウンター1–19におけるPSPの径方向磁場（$B_r$）および径方向太陽風速度（$v_r$）（FIELDSおよびSWEAP機器）。
   • 処理: PSPデータに対して、外れ値除去（四分位範囲基準）および時間平均化（20分間隔）を行う。パーカー・スパイラルによるバックマッピングを連続的に行うため、太陽風速度の欠損値は線形補間によって補完される。
   • 検証: 1 AU付近のAdvanced Composition Explorer（ACE）からの独立した in situ 観測を用いてクロスバリデーションを行う。

3. 最適化アルゴリズム:

   • 目的関数: アルゴリズムは、観測されたPSPの $B_r$ とPFSS予測の $B_r$ の間の平均二乗誤差（MSE）を最小化する。PFSS予測は、宇宙機の位置をパーカー・スパイラルに沿ってソース表面まで遡り、PFSS磁場をサンプリングし、宇宙機の距離まで径方向にスケーリングすること（$r^2 B_r = \text{const}$）で導出される。
   • 正規化: エンカウンターに依存する振幅の違いを取り除くため、MSEの計算前に、観測およびモデルの両方の磁場に対してZスコア正規化を行う。
   • ソルバー: 勾配ベースの反復スキーム（アルゴリズム1）が採用されている。勾配は中心差分を用いて近似される。アルゴリズムは、収束するか最大反復回数に達するまで反復を行う。
   • 検証: アルゴリズムは、解析的なポテンシャル磁場解およびガウスノイズを加えた合成データに対して検証され、規定された $R_{ss}$ を回収できる能力を示した。

4. 多目的解析:

   • 相互トレードオフを理解するために、著者らは2つの追加指標、すなわち磁場スケーリング係数（$\alpha_k$、開いた磁束に関連）と極性予測精度（$P^k_a$）を導入している。
   • パレート分析（アルゴリズム2）を実行し、MSE目的関数と極性予測誤差の間のパレート境界を特定する。これにより、最適化が主に振幅の一致によって駆動されているのか、あるいは極性の正確さによって駆動されているのかを定量化する。

主な結果

• 最適 $R_{ss}$ の傾向: PSPエンカウンター1–19における最適化された $R_{ss}$ は、太陽極小期から第25太陽周期の昇時期にかけて一般に増加する。この傾向は、独立したACEのデータによって裏付けられている。
• 基準との比較: 固定された基準値 $R_{ss} = 2.5 R_s$ と比較して、最適化された解は、磁場の振幅（開いた磁束の不一致）との一致が概ね改善されており、磁極予測の精度を維持または向上させている。
• 太陽極小期 vs. 昇時期:
  • 太陽極小期（初期のエンカウンター）において、最適化は主に振幅の一致（開いた磁束）によって駆動される。観測に一致するように開いた磁束を増やすため、より低い $R_{ss}$ 値が選択される。
  • 昇時期には、パレート境界がより密集するようになり、極性予測の精度が支配的な指標になることを示している。複雑なセクター構造や極性反転をより良く捉えるために、最適な $R_{ss}$ は増加する。
• 局所的な偏差: 特定のエンカウンター（例：9および13）は、一般的な傾向から局所的に逸脱している。エンカウンター13については、複雑な惑星間コロナ質量放出（ICME）を通過したため、複雑な極性反転を捉えるために非常に低い $R_{ss}$（$1.27 R_s$）が選択された。これは、乱れた環境におけるパーカー・スパイラルの仮定の限界を浮き彫りにしている。

意義と主張
本論文は、PFSS逆問題に対する数学的に厳密な基礎を提供し、ソース表面半径の推定における一意の最小化子の存在を証明すると主張している。解析解および in situ データに対して検証された特定の最適化アルゴリズムを構築することにより、本研究は $R_{ss}$ が静的な定数ではなく、太陽活動とともに進化する動的なパラメータであることを示している。

著者らは、単一の「全域的に最適な」 $R_{ss}$ は存在せず、むしろ、特定の観測的制約や優先される評価指標（例：開いた磁束 vs. 極性）に依存することを強調している。また、振幅の一致（開いた磁束）から極性の一致へと、支配的な指標が変化することを明らかにしている。さらに、本論文は「開いた磁束問題」に言及しており、最適化によって in situ 磁束との一致は改善されるものの、コロナホールに基づく磁束推定と in situ 測定との間の根本的な不一致は解消されないことを指摘している。これは、将来的な改善には、より優れた極磁場観測や、より高度なコロナモデルが必要であることを示唆している。本研究は、Solar Orbiterのようなミッションからのデータが入手可能になるにつれ、非球状のソース表面の形状最適化へと発展する枠組みを提示している。