[論文レビュー] Iterative inversion of synthetic travel times successful at recovering sub-surface profiles of supergranular flows

Jishnu Bhattacharya, Shravan Hanasoge+2 more|arXiv (Cornell University)|Aug 11, 2017

Solar and Space Plasma Dynamics参考文献 34被引用数 3

ひとこと要約

本研究では、太陽の表面下における2次元スーパーグランールーブの流れの垂直断面を回復するための反復的逆問題手法を提示する。流れのストリーム関数をBスプライン基底関数で展開し、反復的に旅行時間の不一致を最小化することで、f-p₄モードを用いて表面下8–10 Mmまでのピーク深さおよびリターンフロー深さを、ノイズなしの条件下でも正確に回復する。これは、点別逆問題に比べて優れた性能を示している。

ABSTRACT

We develop a helioseismic inversion algorithm that can be used to recover sub-surface vertical profiles of 2-dimensional supergranular flows from surface measurements of synthetic wave travel times. We carry out seismic wave-propagation simulations through a 2-dimensional section of a flow profile that resembles an averaged supergranule, and a starting model that has flows only at the surface. We assume that the wave measurements are entirely without realization noise for the purpose of our test. We expand the vertical profile of the supergranule stream function on a basis of B-splines. We iteratively update the B-spline coefficients of the supergranule model to reduce the travel-times differences observed between the two simulations. We carry out the exercise for four different vertical profiles peaking at different depths below the solar surface. We are able to accurately recover depth profiles of four supergranule models at depths up to $8-10\, ext{Mm}$ below the solar surface using $f-p_4$ modes, under the assumption that there is no realization noise. We are able to obtain the peak depth and the depth of the return flow for each model. A basis-resolved inversion performs significantly better than one where the flow field is inverted for at each point in the radial grid. This is an encouraging result and might act as a guide in developing more realistic inversion strategies that can be applied to supergranular flows in the Sun.

研究の動機と目的

• 表面での旅行時間測定値から、太陽の表面下における2次元スーパーグランールーブの流れの垂直断面を回復できるヘルモセイズミック逆問題アルゴリズムの開発。
• 基底関数によるパrameter化を用いた反復的逆問題が、標準的な点別逆問題に比べてスーパーグランールーブ構造の解像度を向上させることの検証。
• 理想的なノイズなし条件下で、主なスーパーグランールーブのパラメータ（ピーク深さおよびリターンフロー深さ）を回復できる能力の評価。
• 現実的なスーパーグランールーブに類似した流れ断面の合成波動伝播シミュレーションを用いて、逆問題手法の妥当性の検証。

提案手法

• 合成スーパーグランールーブの流れ断面を通じた2次元波動伝播シミュレーションを実施し、初期モデルでは表面流れのみを仮定する。
• 流れのストリーム関数をBスプライン基底関数でパラメータ化することで、滑らかで柔軟な再構成が可能となる。
• 合成旅行時間と観測旅行時間の不一致を最小化するように、反復的最適化手順でBスプライン係数を更新する。
• 不一致は、正規化されたL₂ノルムおよびストリーム関数、径方向および水平速度成分の深さ依存誤差指標を用いて定量化する。
• ノイズのない状態でf-p₄モードを用い、ピークが異なる深さ（最大14 Mm）にある4つの合成モデルに対して逆問題を適用する。
• 正規化された不一致の低減度を指標としてモデルの改善を評価し、逆問題で得られた断面と真のモデルとの深さ方向の比較を実施する。

実験結果

リサーチクエスチョン

• RQ1 基底関数によるパrameter化を用いた反復的逆問題は、点別逆問題に比べてスーパーグランールーブの流れの深さ断面をより正確に回復できるか？
• RQ2 f-p₄モードを用いて、どの深さまで垂直速度のピーク深さを信頼性を持って回復できるか？
• RQ3 水平速度が方向を逆転するリターンフローの深さは、正確に再構成できるか？
• RQ4 真のピークが14 Mmに位置する場合、この手法は断面をどの程度正確に回復できるか？（これは、従来の地震的検出限界を超える。）
• RQ5 ノイズのない条件下で、基底関数表現が収束性および精度に与える影響は何か？

主な発見

• 反復的逆問題は、最終反復で正規化L₂不一致が約0.01％にまで低下し、真の垂直断面を正確に回復した。
• すべてのテストモデルにおいて、垂直速度のピーク深さが正確に回復された。特に、表面下8–10 Mmにピークを持つモデルに対しても同様の精度を示した。
• 水平速度が逆転するリターンフローの深さは、良好な精度で回復されたが、正確な断面形状までは完全に再構成できなかった。
• f-p₄モードを用いた場合、ピークが14 Mmに位置するスーパーグランールーブの断面は回復できなかった。これは、このモードセットでは深部構造の解像度に限界があることを示唆している。
• 基底関数によるパラメータ化を用いた逆問題は、点別逆問題を上回る性能を示した。これは、流れを低次元の基底空間にパラメータ化することで、安定性と精度が向上することを示している。
• 今後の実用的逆問題戦略の開発にとって前向きな結果であり、特にノイズが少ないまたはノイズなしの条件下での応用が期待される。