[论文解读] Polar Field Correction for HMI Line-of-Sight Synoptic Data

Xudong Sun|arXiv (Cornell University)|Jan 12, 2018

Distributed and Parallel Computing Systems被引用 8

一句话总结

本文提出一种时空插值方案，用于校正HMI视线合成图与同步太阳磁像数据中缺失或低质量的极区磁场。通过使用每年最大极向观测几何条件下的基准观测数据，并应用五阶切比雪夫多项式拟合与线性时间插值，该方法提升了极区磁场估计的准确性，将75°纬度以上区域的均方根误差降低至0.73–0.99 G，并通过62°–75°纬度范围的1.21倍缩放因子，实现了对未来极区磁场的可靠外推。

ABSTRACT

This document provides some technical notes on the polar field correction scheme for the HMI synoptic maps and daily updated synchronic frames. It is intended as a reference for the new data products and for some minor updates on our previous scheme for MDI (Sun et al. 2011).

研究动机与目标

• 解决由于太阳极区附近视线观测受限导致的HMI极区磁场测量系统性低估与噪声问题。
• 为HMI合成图与每日同步帧开发一种一致且具有物理解释性的校正方案，同时考虑太阳自转与观测几何条件。
• 利用历史数据实现对未来极区磁场强度的可靠外推，提升空间天气与太阳周期预测模型的性能。
• 通过在最大极向观测条件下对齐测量值，校正MDI与HMI数据中潜在的偏差。
• 提供一个随新观测数据不断更新的、带时间戳的确定性校正数据库，确保数据随时间保持一致性。

提出的方法

• 每年基于最大极向观测几何条件下的最佳观测数据构建极区磁场数据库——北极为9月6日02:00 UTC，南极为3月8日02:00 UTC。
• 从完整卡林顿周期合成图中提取纬度高于60°的极区区域，并将卡林顿经度180°对齐至第1800列。
• 将数据重映射至立体投影极坐标系，采用五阶切比雪夫多项式进行平滑与插值，再将结果重映射回原始圆柱等面积（CEA）坐标系。
• 对于每张合成图或同步帧，基于对应极区最近的两个数据库记录（T₀ ≤ T ≤ T₁）进行线性时间插值以应用校正。
• 仅对纬度高于75°的像素应用校正，用插值结果替换原始值，以降低噪声并填补数据空缺。
• 前向外推模型使用62°–75°纬度范围内年份N的平均磁场值，乘以1.21倍缩放因子，以预测年份N+1中75°纬度以上区域的平均磁场值。

实验结果

研究问题

• RQ1 如何系统性地校正因视线观测受限而导致的HMI合成图与同步图中缺失或噪声严重的极区磁场数据？
• RQ2 在最大观测几何条件下，使用基准观测数据，最优的时空插值方法是什么？
• RQ3 该时间插值方案在重建真实极区磁场演变过程中的准确性如何，特别是在高倾角与数据丢失期间？
• RQ4 能否推导出一个一致的缩放因子，从当前观测外推未来极区磁场强度，从而提升太阳周期预测能力？
• RQ5 MDI与HMI极区磁场测量中存在何种偏差？如何通过仪器间对比加以校正？

主要发现

• 在合成图中，75°纬度以上区域，校正后与原始HMI数据之间的均方根（rms）差异为0.73 G（北半球）与0.98 G（南半球）。
• 在每日同步帧中，rms误差为0.70 G（北半球）与0.99 G（南半球），当极区远离地球（b ≤ -5°）时，北半球误差显著上升至1.34 G。
• 校正后的极区磁场值与直接磁像测量值高度一致，验证了该方法的准确性。
• 使用62°–75°纬度范围1.21倍缩放因子的前向外推模型，对未来极区磁场的预测均方根误差为1.00 G，2000年、2004年与2016年存在异常值，归因于磁场的快速变化。
• 校正后的HMI数据中，南半球极区磁场略强于原始值，可能由于2017年3月基准记录存在高估。
• 截至2017年12月，北半球校正数据已确定至2017年9月，南半球至2017年3月，后续数据标注为暂定，将随新数据更新。