[论文解读] 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

一句话总结

本文提出了一种非线性正则化方法用于日震反演，以更好地解析太阳较差自转层中的陡峭梯度，取代传统的线性Tikhonov正则化。该方法借鉴图像处理中的边缘保持正则化技术，直接恢复薄过渡层，而无需预先假设函数形式，从而从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（自适应正则化结合Tikhonov与无偏风险）方法以平衡数据拟合与正则化。
• 通过广义交叉验证（GCV）和蒙特卡洛模拟校准正则化参数，以评估不确定性。
• 将最终解拟合到上较差自转层的误差函数中，并引入线性趋势，以更好地捕捉剖面形状。

实验结果

研究问题

• RQ1 与线性Tikhonov正则化相比，非线性正则化方法是否能更有效地保持日震反演中的陡峭梯度？
• RQ2 在不预先假设函数形式的前提下，解析高梯度特征后，太阳较差自转层的真实宽度是多少？
• RQ3 在上较差自转层中引入线性趋势如何影响旋转速率剖面的拟合与解释？
• RQ4 与以往方法相比，蒙特卡洛模拟在多大程度上降低了较差自转层宽度估计的不确定性？
• RQ5 该非线性方法是否可推广至其他涉及高梯度结构的日震反演问题？

主要发现

• 较差自转层宽度估计为小于0.05 R⊙，显著窄于以往估计的0.05 ± 0.03 R⊙。
• 非线性正则化方法成功恢复了旋转速率的陡峭过渡，而无需预先假设剖面的具体函数形式。
• 该解揭示了在上较差自转层从0.7 R⊙到0.8 R⊙之间，旋转速率呈约70 nHz/R⊙的线性增加，该特征在传统线性反演与去卷积方法中无法检测到。
• 蒙特卡洛模拟证实了宽度估计的稳健性，并相比早期方法显著降低了不确定性。
• 该方法可直接恢复高梯度特征，消除了对依赖于假设函数形式的反演后去卷积的依赖。
• 该方法适用于其他涉及陡峭结构特征的日震问题，例如较差自转层区域的声速异常。