Combining a Geomagnetic Precursor with Dynamic Time Warping for Solar-Cycle Peak and Trajectory Forecasting

太阳并非一盏稳定不变的灯，它是一颗以大约十一年的节奏进行呼吸的躁动恒星。这种被称为太阳周期的“心跳”，决定了出现在太阳表面的黑子数量，并驱动着围绕地球的空间天气强度。当太阳活跃时，它会向地球发送强大的辐射和磁场脉冲，可能干扰卫星、无线电通信和电网。对于必须保护技术的科学家和工程师来说，预测该周期何时达到顶峰以及将如何演变，是一个至关重要的挑战。难点在于时机：在一段新周期的起始阶段，可用的太阳黑子数据非常少，因此很难猜测即将到来的风暴会有多强。相反，随着周期的成熟，黑子的模式变得清晰，但关于风暴总强度的早期线索往往会被遗失或被遗忘。

两位研究人员，Varun Soman 和 Amaranathan Sabarinath，开发了一种结合两种不同类型历史线索的新方法，以弥补这一差距。一种线索来自地球磁场，即使在新的黑子周期完全发展起来之前，它也能作为探测太阳活动的一个敏感探测器。另一种线索则来自回顾过去：寻找与当前正在展开的周期非常相似的前期太阳周期，并利用它们的未来路径作为指南。通过将这两类信息编织在一起，团队创建了一种比仅使用其中一种线索更准确的预测方法。他们的工作在最后四个已完成的太阳周期中进行了测试，结果表明，这种混合方法可以比以往的方法以显著更低的误差预测太阳周期的峰值强度，同时还能提供一个更真实的周期起伏过程图景。

他们方法的核心依赖于一种称为“动态时间规整”（dynamic time warping）的技术，这本质上是一种灵活的形状比较方式。想象一下尝试匹配两个不同形状的山丘；其中一个可能高而窄，而另一个可能矮而宽。一种僵化的比较可能会因为两者的峰值发生的时间不同或持续时间不同，而无法识别出它们在形状上的相似性。动态时间规整允许计算机拉伸或压缩其中一个山丘的时间轴以适应另一个，从而在不强制两者长度相等的情况下，找到曲线之间最佳的匹配。在本研究中，研究人员使用这一工具将正在发展的太阳周期的早期部分形状，与一个完整的历史周期库进行对比。他们并没有强制要求历史周期缩短以匹配当前的数据，而是保持了历史周期的完整原始长度，并让匹配算法找到最佳对齐方式。这一特定的选择被证明是提高预测准确性的最重要因素。

在此方法出现之前，研究人员通常会将历史周期缩短，以匹配预测时刻所拥有的数据长度。这项研究发现，这种重采样（resampling）的做法实际上扭曲了比较过程，导致预测准确度降低。通过保持历史周期的完整性并使用灵活的匹配工具，研究人员显著降低了预测未来太阳黑子路径的平均误差。当他们在太阳周期 21 至 24 上进行测试时，新技术的预测轨迹平均误差约为 20 个平滑黑子，而使用旧的重采样方法时的误差约为 26 个。这种改进并非源于更复杂的数学计算，而仅仅是因为该方法保留了历史样本的真实长度和形状，从而使算法能够选择更好的匹配项。

第二个重大创新在于团队如何结合两种不同类型的预测。源自地球磁场的地磁前兆在给出周期强度早期估计方面表现出色，但它无法提供太多关于曲线形状的信息。而模拟法（analog method）通过寻找相似的过去周期来运作，随着更多黑子数据的积累，它在预测形状方面会变得越来越好，但在初期对于总强度的判断可能存在不确定性。研究人员创建了一个系统，根据每种预测的历史不确定性来分配权重。当周期刚刚开始时，地磁线索具有更大的权重，因为它是关于峰值强度的唯一信息来源。随着周期的发展和更多黑子数据的累积，权重会向模拟法转移，从而提供更详细的周期演化图景。这种动态平衡的行为使得结合后的峰值预测平均误差约为 16 个平滑黑子，相比于仅使用地磁线索产生的 25 个黑子的误差，有了显著提升。

团队解决的另一个实际问题涉及太阳周期的“下降阶段”。当预测结果根据最新的观测数据进行调整时，简单的数学偏移有时会将预测的黑子数量推至零以下，尤其是在周期的衰减年份。由于黑子数量不可能为负数，这会导致产生一个不符合现实的、人为的零值平线。研究人员引入了一条规则来防止这种情况，确保预测轨迹在物理上是可能的，即仅在必要时才通过调整数学计算来保持数值为正。这一改变并不一定会使历史预测变得更精确，但它确保了预测结果保持现实感，且没有出现不可能出现的伪影。

该研究还考察了该方法处理不确定性的能力。通过基于历史误差生成数千条可能的未来路径，团队创建了显示可能结果范围的“包络线”（envelopes）。他们发现，虽然这些新的包络线通常比旧方法的包络线更窄、更精确，但它们并非完美。对于某些周期（如太阳周期 22），该方法难以捕捉快速下降的过程，导致预测与现实之间存在较大的差距。这表明，尽管该方法是一个稳健的工具，但它并不是一颗预知未来的水晶球；太阳的行为仍可能令我们感到意外，特别是当一个周期的行为与现有的历史案例迥异时。研究人员强调，应将这些不确定性带视为经验性的离散度度量，而非保证的极限，因为用于测试的历史周期数量较少，很难确定每一个可能的结果。

展望未来，研究人员为应用此方法于即将到来的第 26 号太阳周期制定了严格的计划。他们决定在周期开始前固定模型的所有规则和参数，以确保在结果已知后不会做出任何更改。这种方法将是对其在实时环境下预测新周期能力的真正考验。目前，该方法已被应用于正在进行的第 25 号太阳周期。对该周期的分析展示了综合估算值是如何演变的：从一个地磁猜想开始，随着更多黑子数据的加入，预测结果逐渐得到精炼。研究结论指出，通过将早期的磁场线索与灵活的历史“孪生兄弟”搜索相结合，科学家现在可以提供更可靠的太阳行为预测，从而更好地应对影响现代世界的空间天气。