论文阅读：Comparison of Line-of-Sight Magnetograms Taken by the SDO/HMI and SOHO/MDI

一、简介

HMI数据包含了所有的MDI观测值，但具有更高的空间和时间分辨率以及更好的数据质量。为了将MDI的研究和空间天气模拟推广到HMI，有必要对两种仪器的观测数据进行比较。本文主要比较了SDO/HMI和SOHO/MDI的视向磁图。

• 对于磁场信号，MDI比HMI的大1.4倍左右，随着中心到边缘的距离变化而变化
• 对于随机噪声的上界，HMI如表所示。

1"resolution HMI	720s	6.3Mx/cm^(-2)
	45s	10.2Mx/cm^(-2)

没有表现出p模泄漏。在12h和14h时都表现出周期性。
MDI采用重新标定（零点偏移和磁通密度的修正）后的数据。

4"resolution MDI	1min	26.4Mx/cm^(-2)
	5min	16.2Mx/cm^(-2)

对于存在的噪声，可能是由于左右的圆分析器的光学畸变引起的粒化和p模泄漏。

• 对于line-of-sight magnetic-field全景图（synoptic charts），表现的比较均匀。HMI和MDI分别为2.3和5.0。HMI没有表现出明显的周期性。

本文的组织如下：

• 第2部分对HMI视向磁图进行了详细的描述，包括推导磁图所用的算法、噪声水平以及问题。
• 第3节是对MDI磁图进行的修正，例如对零点偏移的修正和对空间依赖的磁通密度低估的修正。
• MDI和HMI磁图的比较见第4节。

二、HMI Line-of-Sight Magnetograms

2.1 相机介绍

HMI是全日面4096 × 4096像素。空间分辨率为1′′，像元尺寸为0.5′′。谱线为光球中形成的FeI 6173，滤波轮廓的宽度为76 mA。该仪器有两个CCD摄像机：“正面摄像机"和"侧面摄像机”。
前置相机在左旋圆偏振( LCP )和右旋圆偏振( RCP )两种偏振状态下采集跨越谱线的6个波长的滤波图，图像之间间隔3.75秒。获取一组完整的12个滤波图需要45。导出多普勒图和视线磁图。
侧面相机专门用于测量矢量磁场。当前帧列表需要135秒来获得6个波长位置的6个偏振态的滤波图。由这些测量值计算出斯托克斯参量[ I、Q、U、V]，并进一步得到翻转、去模糊后的矢量磁场的反演。

2.2 视向磁图

HMI视向磁图是由LCP和RCP的滤波片差值导出的多普勒速度。每个滤波图都通过去除宇宙射线碰撞、填充数据间隙、校正光学畸变等方法进行了校正。
对于45秒磁图，LCP和RCP的滤波片考虑了太阳自转，并插值到目标时刻，然后将它们组合产生多普勒图。时间插值通常使用6个点：只要有可能，目标时间之前3次，目标时间之后3次。因此时间间隔为6 × 45秒。对每个谱片施加一个权重，根据谱片在时间上与目标时间的接近程度而变化。速度算法分别应用于LCP和RCP。
对于720秒的磁图，在计算多普勒速度之前，对配准后的单色像图进行时间平均。这种平均需要在720秒时间窗口前后增加额外的单色像图，因为在平均之前需要将滤波图插值到一个规则的、更精细的网格上。平均法采用余弦切趾的矩形框，其 FWHM为720秒。
最后应用速度算法。在实际应用中，所有的过程实际上是结合在一起的。

2.3 磁噪声

按照MDI磁图( Liu , Zhao , & Hoeksema , 2004)的方法，粗略地估计了HMI视向磁图的噪声水平。假设磁图中的噪声具有高斯分布，并使用高斯函数来拟合磁图低场像素的分布。

取最佳高斯拟合的σ (标准差)为噪声水平，45秒磁图的噪声水平为10.2 Mx/cm^(-2)， 720秒磁图的噪声水平为6.3Mx/cm^(-2)。采用同样的方法，MDI 1分钟和5分钟磁图的噪声分别为26.4Mx/cm^(-2) 和16.0Mx/cm^(-2)。当然，噪声被高估了，因为这个分布包含了真实的太阳弱场区域的观测。
对全日面上的噪声也进行了估计。

噪声相对于磁盘中心看起来对称，向边缘逐渐增加。图中的水平条纹表明，用这种方法产生的噪声也包括磁场信号。
还估计了噪声对中心到边缘距离的依赖性。

该方法计算得到的这种中心到边缘的噪声变化可能部分来自于磁场的真实信号：当距离边缘越近时，200高斯内包含更真实的视线场，用于拟合高斯分布。因此拟合高斯噪声的σ大于真实噪声。

2.4 周期

由于轨道速度的24小时周期性和仪器校准中的不确定度，在强HMI磁场中观察到日变化和半日变化。相应的视向磁图不准确可能造成12小时和24小时周期，而更好的校准曲线可以最大限度地减少这种影响。

由周期性引起的变化平均约为信号的2.7 %。
当平均场的功率谱减去一段时间内的平均值，结果没有明显的周期性。这意味着单独将太阳自转和轨道速度结合起来并没有引入显著的周期性。

三、MDI Magnetograms

3.1 相机介绍

MDI的主要观测量，即多普勒速度，由4个单色像图得到。视线磁图是利用LCP和RCP的两个多普勒图的差来确定的。MDI可以在两种空间分辨率下进行观测。高分辨率路径的视场为11弧分× 11弧分，以0.62角秒的像素尺寸略微向磁盘中心的北部偏移。通常情况下，MDI以96分钟的速度产生全日面磁图。有时，当高速率遥测可用时，MDI运行运动程序，以提供全日面或高分辨率磁图，频率为1分钟。获得一幅磁图大约需要30秒钟。MDI在其96分钟的频率时间序列中产生了两种类型的全日面磁图。一分钟磁图由在30秒时间间隔内收集的一组滤波图产生。五分钟磁图是在五分钟时间间隔内获得的五个一分钟磁图的平均值。

3.2 全日面磁图标准化

3.2.1偏移量

相机快门曝光时间的微小不确定性会引起测量结果的误差。这种快门噪声为每个成像可观测量增加了一个近乎均匀的背景偏移量。
对于MDI磁图，使用这种算法( Liu , Zhao , & Hoeksema , 2004)来消除这种偏移。这种方法对于改进大尺度磁场的测量和建模是有效的。

3.2.2场通量密度的低估

MDI数据的通量密度值在活动区中比ASP数据低0.64倍，在半影区和半影区低0.69倍。Tran等人( 2005 )通过比较MDI和Wilson山天文台( Mount Wilson Observatory，MWO )的全日面磁图，在较弱的磁场区域解决了这个问题。主要利用同时测量磁敏感谱线和MDI谱线中MWO处的视线磁场进行综合分析，交叉关联。利用因子[ δ-1 = 4.5-2.5 × ( 1 - cos2 ρ)]对MWO Fe I 5250磁图进行线轮廓饱和校正。这个因子最近被修正为δ - 1 = 4.15-2.82 × ( 1 - cos2 ρ) - ( Ulrich et al , 2009)。

实际用于校正MDI全日面磁图的比值是由图6所示的比值图经新的饱和因子校正后得到的平滑掩膜。

3.2.3非均匀噪声

由于滤波图是组合在一起的，所以MDI数据中的失真不能确定。CCD的右下方是噪声较大的地方。
图中显示了太阳磁盘上的噪声，数据来自于2007年3月拍摄的约32 000幅1分钟磁图。对于每个像素，使用高斯函数来拟合32 000次测量的低场像素的分布。

图中显示了高斯函数的标准差，从磁盘中心到边缘逐渐增加。在右下方区域最高。中位数为26.4 Mx cm - 2

3.2.4饱和现象

太阳黑子本影强度低和15位数值精度的限制，导致MDI磁图中某些黑子本影出现饱和现象。由于MDI强度数据也可以在较低的频率下获得，因此可以利用已建立的连续体强度与磁场强度( I - B关系) 之间的关系来校正磁图中的饱和现象。

3.2.5耀斑期间的磁场反转

在大耀斑期间，MDI磁象图中有时会出现磁场反向的现象。这种反转是暂时的，只在耀斑期间观察到。通常反转发展出与耀斑亮带相同的模式。它不太可能是磁场的真实变化。到目前为止，还没有方法对这种影响进行修正。

四、Comparison of HMI and MDI Magnetic Field Observations

4.1 比较

在逐像素的基础上，比较了HMI和校准的MDI磁图。方法如下所述。

• 首先通过卷积二维高斯函数来降低HMI的空间分辨率，使其与MDI的空间分辨率相当。使用的高斯函数的FWHM为4.7 HMI像素，并被15个HMI像素截断。
• 然后，将每个两角秒的MDI像素映射到太阳表面。从相应的模糊HMI磁图中平均出包含在该区域内的HMI像素。仅部分包含在区域内的HMI像素通过该分数进行调整。通过这种方式，生成了一个MDI代理像素。将这些代理像素组合在一起，我们得到了一个MDI代理磁图。
• 对于每个MDI磁图，校正了畸变和位置角[ P角]的偏移。取0.22 ° (顺时针方向)作为MDI磁图的P角偏移量。以P角偏移校正的MDI磁像图为参考，通过计算MDI磁像图与HMI磁像图之间的相关系数来确定可能需要应用于HMI磁像图的旋转角度。当HMI 45 s磁图顺时针旋转0.05 时，相关系数达到最大；当HMI 720 s磁图顺时针旋转0.10 时，相关系数达到最大。

4.2 检查

HMI 8小时45秒磁图(左)和MDI高分辨率磁场数据(右)的功率谱。
功率谱是由靠近磁盘中心的一个安静区域制成的。MDI谱图中的脊线表明磁图中存在显著的p模泄漏。这可能是由于左、右圆分析器的光学畸变不同造成的。HMI磁场数据的功率谱中没有发现p模脊。

对太阳边缘的数据质量也进行了检验。MDI 1分钟磁图(顶部)、HMI 45秒磁图(中部)和HMI 720秒磁图(底部)从北60 °到北极的太阳区域。许多在MDI磁图中不可见的小尺度磁性元素在HMI磁图中清晰可见。

新学的单词

filtergram 单色像
longitudinal field 纵场
transverse field横场
synoptic chart 全景图