Towards inertial-mode helioseismology: Direct sensing of solar rotation at 75 deg latitude and 0.8 Rsun

技术摘要：迈向惯性模日震学

问题陈述
传统 p 模日震学对太阳内部高纬度及对流区深处的较差自转约束仍然不足。虽然 p 模能有效绘制对流区主体的自转分布，但其灵敏度核在太阳表面附近达到峰值，因此关于旋转轴附近约 $20^\circ$ 范围内及 $0.2R_\odot$ 以下区域的信息有限。因此，对 p 模频率分裂的反演在极区产生了显著差异。作者提出，由科里奥利力恢复且与自转内在关联的全球太阳惯性模，为这些难以探测的区域提供了一种有前景的替代诊断手段。具体而言，本研究聚焦于方位波数 $m=1$ 的高纬度惯性模（即“HL1"模），该模在观测到的惯性模中表现出最大的表面振幅。

方法论
本研究采用正演建模方法，利用 HL1 模的测量频率来约束太阳自转。

1. 参考模型：作者采用基于 p 模日震学（Larson & Schou 2018）导出的参考较差自转剖面（$\Omega_{\text{ref}}$），对应时段为 2010–2024 年。
2. 本征值求解器：利用经过验证的本征值求解器（求解器 1），他们在 $0.65R_\odot < r < 0.985R_\odot$ 的计算域内计算了 HL1 模的复本征频率。该求解器已与独立实现（求解器 2）进行了交叉验证。
3. 灵敏度核计算：为了确定局部自转变化如何影响模态频率，作者数值计算了线性灵敏度核 $K(r, \theta)$。这涉及在特定点网格处扰动自转速率（$\delta\Omega$），同时通过纬度熵梯度调整维持热风平衡，并计算由此产生的频率偏移（$\delta\omega$）。
4. 参数化修正：为了超越线性近似，作者在高纬度自转速率中引入了一种平滑的、依赖纬度但与半径无关的参数化修正（$\Delta\Omega$）。他们重新计算了有限振幅扰动下的本征频率，以映射非线性响应。

关键结果

• 模型与观测的不匹配：使用 p 模参考自转计算出的模型本征频率（$\omega_{\text{ref}}/2\pi = -95.6$ nHz）与卡林顿参考系下的观测平均频率（$\omega_{\text{obs}}/2\pi = -87.9 \pm 1.9$ nHz）存在显著差异。约 $7.7$ nHz 的偏差无法通过改变湍流扩散率或超绝热性等背景参数来解决，这表明自转剖面是误差的主要来源。
• 灵敏度定位：HL1 模的线性灵敏度核高度局域化。其在纬度 $75^\circ$（余纬 $15^\circ$）和半径 $0.8R_\odot$ 处达到峰值，纬度方向的半高全宽为 $7^\circ$，半径方向的半高全宽为 $0.13R_\odot$。该模态对辐射内部的灵敏度可忽略不计。
• 自转速率推断：
  • 局部估计：通过对子午面积分核函数，作者推断在 $(0.8R_\odot, 75^\circ)$ 处的局部自转扰动为 $\delta\Omega/2\pi \approx 8.1 \pm 2.0$ nHz。这意味着局部自转速率为 $365.3 \pm 2.0$ nHz，比 p 模参考估计值（$357.2$ nHz）高出 $8.1$ nHz。
  • 全局参数化估计：使用平滑的参数化修正，最佳拟合模型得出的极区自转增强量为相对于参考模型的 $13.6 \pm 3$ nHz，对应极区自转速率为 $364.0$ nHz。
• 与 p 模的比较：目标位置推断出的自转速率显著高于基于 p 模最优局部平均（OLA）反演得出的估计值，突显了惯性模在高纬度区域提供的独特且可能更准确的约束。

意义与主张
本文声称展示了首个空间分辨的惯性模日震学实例。其主要意义在于证明单个惯性模可以直接、绝对地约束对流区主体、远低于表面且 p 模不敏感的高纬度区域的太阳自转。与依赖顺行模和逆行模之间差分测量的 p 模日震学不同，该方法利用单个逆行惯性模的绝对频率。作者得出结论，经过验证的本征值求解器是解释完整惯性模谱的关键工具，未来的工作将旨在利用多个模态同时求解反演问题，以更全面地绘制太阳较差自转图。