[論文レビュー] Non linear regularization for helioseismic inversions. Application for the study of the solar tachocline

T. Corbard, Laure Blanc‐Féraud+2 more|arXiv (Cornell University)|Jan 11, 1999

Solar and Space Plasma Dynamics参考文献 13被引用数 29

ひとこと要約

本論文は、太陽のタコクリン層における急激な勾配をよりよく解像するため、従来の線形チホノフ正則化に代わって非線形正則化を導入する。画像処理分野のエッジ保存型正則化を応用することで、事前に関数形を仮定することなく、薄い遷移層を直接回復可能となり、LOWLデータおよびモンテカルロシミュレーションから得られる結果では、タコクリン層の幅が0.05 R⊙未塔であり、以前の推定値0.05 ± 0.03 R⊙よりもはるかに薄いことが示された。

ABSTRACT

Inversions of rotational splittings have shown that there exists at the base of the solar convection zone a region called the tachocline in which high radial gradients of the rotation rate occur. The usual linear regularization methods tend to smooth out any high gradients in the solution, and may not be appropriate for the study of this zone. In this paper we use, in the helioseismic context of rotation inversions, regularization methods that have been developed for edge-preserving regularization in computed imaging. It is shown from Monte-Carlo simulations that this approach can lead directly to results similar to those reached by linear inversions which however required some assumptions on the shape of the transition in order to be deconvolved. The application of this method to LOWL data leads to a very thin tachocline. From the discussions on the parameters entering the inversion and the Monte-Carlo simulations, our conclusion is that the tachocline width is very likely below 0.05R_sun which lowers our previous estimate of 0.05+/- 0.03R_sun obtained from the same dataset (Corbard et al. 1998).

研究の動機と目的

• 線形正則化が急勾配領域（太陽のタコクリン層など）を解像する際に抱える限界を克服すること。
• 画像処理分野の非線形正則化技術を発展・応用し、ヘルモセイズミック逆問題に適用すること。
• 遷移層の関数形を事前に仮定しない形で、LOWLデータからタコクリン層の幅を特定すること。
• 線形逆問題のデコンボリューションに必要な仮定を回避することで、タコクリン層幅の推定誤差を低減すること。

提案手法

• 画像処理分野のエッジ保存型原理に基づく非線形正則化を採用し、1-ノルムペナルティを用いて区分的滑らかさを促進する。
• 基底関数を半径グリッド点に定義することで、離散的最適化問題として逆問題を定式化し、解を区分的線形に保証する。
• 一次微分作用素を二重対角行列として用い、勾配を計算し、空間的変動を重み付けする対角行列を導入する。
• データ適合性と正則化のバランスをとるために、ARTUR（適応的正則化法：チホノフと不偏リスクを組み合わせた手法）を適用する。
• 一般化交差検証（GCV）およびモンテカルロシミュレーションを用いて正則化パラメータをキャリブレーションし、誤差を評価する。
• 最終解を、上部タコクリン層における線形傾向を含む誤差関数にフィットさせ、プロファイル形状をよりよく捉える。

実験結果

リサーチクエスチョン

• RQ1 非線形正則化手法は、線形チホノフ正則化に比べ、ヘルモセイズミック逆問題における急勾配をよりよく保持できるか？
• RQ2 急勾配特徴を事前に関数形を仮定せずに解像した場合、太陽のタコクリン層の真の幅は何か？
• RQ3 上部タコクリン層に線形傾向を含めることで、回転プロファイルのフィットと解釈にどのような影響を与えるか？
• RQ4 モンテカルロシミュレーションを用いることで、従来の手法に比べ、タコクリン層幅の推定誤差はどの程度低減されるか？
• RQ5 この非線形手法は、急勾配構造を含む他のヘルモセイズミック逆問題に一般化可能か？

主な発見

• タコクリン層の幅は0.05 R⊙未塔と推定され、以前の推定値0.05 ± 0.03 R⊙よりもはるかに狭い。
• 非線形正則化手法により、プロファイルの関数形を事前に仮定せず、回転率の急激な遷移を成功裏に回復した。
• 解釈により、上部タコクリン層における0.7 R⊙から0.8 R⊙の間で、回転率が約70 nHz/R⊙の線形増加が確認され、これは線形逆問題とデコンボリューションでは検出できなかった。
• モンテカルロシミュレーションにより、幅の推定値の妥当性が確認され、従来手法に比べて誤差が低減した。
• この手法により、後処理としてのデコンボリューション（関数形の仮定に依存）を回避し、急勾配特徴を直接回復可能となった。
• 本手法は、タコクリン層領域における音速異常などの急峻な構造的特徴を含む他のヘルモセイズミック問題にも適用可能である。