Limb Shift of the Fe I 6569 Ã line on the Sun

技術要約：太陽における Fe I 6569 Å 線のアームシフト（Limb Shift）

問題提起
太陽の光球には、よく知られた「アームシフト（limb shift）」、すなわち、スペクトル線が太陽円盤の縁（リム）付近で中心部に対して見かけ上の赤方偏移を示す現象が存在する。この現象は主に対流運動に起因している。すなわち、明るい上昇流であるグラニュール（青方偏移）が円盤中心部ではスペクトルを支配するが、リム付近では投影効果によってその見かけの視線速度が減少するため、赤方偏移した粒間領域（intergranular lanes）が顕在化する。しかし、この効果の正確な測定は、空間データの平均化手法に依存する。二つの異なるアプローチが存在する：**スペクトル平均化法（手法1）**は、速度を算出する前にスペクトル線プロファイルを空間的に平均化する手法であり、**速度平均化法（手法2）**は、各ピクセルに対してまずドップラー速度を算出し（ドップラーグラムを作成）、その後、それらの速度を平均化する手法である。これらの手法は、低解像度の観測機器によるピクセル畳み込み効果により、歴史的には同様の結果をもたらしてきたが、高空間解像度がこれら二つの手法間の乖離に与える影響については、これまで十分に探求されてこなかった。さらに、観測パイプラインでは本質的な対流青方偏移が無視されることが多く、これが系統的なバイアスを生じさせている。

手法
本研究では、観測データと数値シミュレーションを組み合わせた二角的なアプローチを用いて、Fe I 6569.2 Å 線のリムシフトを調査する：

1. 観測データ： データは、シーイングの影響を受けない分光観測を提供する中国の Hα 太陽探査機（CHASE）ミッションから取得された。本研究では、高いスペクトルおよび空間解能（1.2 アーセカンド）を維持するために、ビンニングを行わずに Fe I バンドに焦点を当てた 2022 年 8 月 9 日の特定のスキャンデータを利用した。太陽円盤は、ヘリオセントリック角（$\mu = \cos \theta$）に基づく同心円状のリングに分割された。
2. 数値シミュレーション： 3D 放射磁気流体力学（MHD）モデルである Bifrost コードを用いて合成データを生成した。シミュレーションボックス（24 $\times$ 24 $\times$ 17 Mm$^3$）には 400 個以上のグラニュールが含まれている。スペクトル線プロファイルは、局所熱平衡（LTE）および完全再配分（CRD）の仮定の下、RH1.5D コードを用いて合成された。
3. 解析手法：
   • 手法 1（スペクトル平均化）： 各ヘリオセントリック角においてスペクトル線プロファイルを空間的に平均化し、その後、平均化されたプロファイルからドップラー速度を導出する。
   • 手法 2（速度平均化）： すべてのピクセルに対してドップラー速度を計算してドップラーグラムを生成し、その後、ヘリオセントリック角にわたってこれらの速度を平均化する。
   • 較正： 両方のデータセットに対し、全域的な速度場（太陽の自転、子午線循環、衛星の運動）および重力赤方偏移を除去する処理を行った。シミュレーションと観測の比較により、系統的な対流青方偏移のオフセットを決定した。
4. 解像度解析： 合成データの空間解像度を、ネイティブのシミュレーション解像度（0.066″）から様々な観測能力に合わせてガウス型点広がり関数（PSF）を用いて系統的に劣化させた。

主要な結果

• 対流青方偏移の定量化： 合成データと観測データの比較により、円盤中心において約 0.36 km/s の系統的な対流青方偏移が特定された。これは、標準的な CHASE 波長較正パイプラインでは考慮されていない。
• 解像度依存性：
  • CHASE の空間解像度（1.2″）において、手法 1 と手法 2 はほぼ同一のリムシフト曲線を生じさせ、観測された傾向（中間角 $\mu \approx 0.65$ で青方偏移が増大した後、リム付近で急速に赤方偏移へ遷移する）と一致する。
  • **高空間解像度（< 1″）**において、二つの手法は著しく乖離する。手法 2 は古典的なリムシフトプロファイルから逸脱し、円盤中心からリムに向かって推測される青方偏移が増大してから減少するという異常な傾向を示す。対照的に、手法 なる 1 は古典的な挙動と一致し続ける。
  • 空間解像度が低下するにつれ（PSF 幅 > 2.0″）、手法 2 の結果は手法 1 に収束し、空間的な平滑化が信号を平均化するため、二つの手法が等価になることを示している。
• スペクトル解像度： リムシフト曲線は、スペクトル解像度が 0.016 Å（$R > 400,000$）より良ければ、スペクトル解像度に対して比較的鈍感である。
• 物理的起源： 本研究では、リムシフトを二つの成分に分解した：
  1. コントラスト寄与： 速度と線深（明るさ）の相関によって引き起こされる。明るいグラニュールは手法 1 においてより大きな重みを持つ。この寄与はリムに向かって減少し、負になることさえある。
  2. ドップラーグラム寄与： 手法 2（等重み付け）においても存在する。これは密度不均一性（約 100 m/s の赤方偏移を導入する）および波打ち効果（corrugation effects）（$\tau=1$ 面が平坦ではないこと）に起因する。水平速度が視線方向に投影されることで生じる波打ち効果は、リムに向かって支配的な青方偏移を生み出す。

意義および結論
本論文は、平均化手法の選択は極めて重要であり、計測機器の空間解像度に依存すると結論付けている。

• **低空間解像度（> 2.0″）**の観測機器の場合、手法 1 と手法 2 の区別は無視でき、どちらの手法もリムシフトの測定やドップラーグラムの補正に適用可能である。
• 高解像度観測機器（例：DKIST, SST/CRISP）の場合、両手法は根本的に異なる結果をもたらす。手法 1（スペクトル平均化）は、伝統的なリムシフトの定義および古典的な観測結果と一致する。一方、手法 2（速度平均化）は、ドップラーグラムの寄与と波打ち効果に支配された異なる物理的実態を明らかにするものであり、それは歴史的に定義されてきた「リムシフト」とは異なる可能性がある。
• 本研究は、先行文献における不一致が、使用される線の違い、バイセクター強度の違い、あるいは低解像度機器による暗黙的な平均化効果に起因している可能性を示唆している。
• これらの知見は、高解像度の太陽観測を解釈するための枠組みを提供する。手法 1 はリムシフト曲線自体の測定に適しており、一方で手法 2 は、結果として得られる曲線が古典的なプロファイルとは異なることを理解した上で、他の現象（差動自転など）を分離するためにリムシフト効果を除去することを目的とする研究に適している可能性がある。

著者らは、これらの結論が Fe I 656 9 Å 線および Bifrost シミュレーションモデルに特有のものであることを注記しており、これらの知見を一般化するためには、追加の光球線を用いた系統的な調査が必要であるとしている。