这篇研究发表于 nature geoscience 期刊，2026年3月第19卷。主要作者包括 A. J. Biggin（第一作者兼通讯作者，来自英国利物浦大学）、C. J. Davies 和 J. E. Mound（均来自利兹大学）、S. J. Lloyd（利物浦大学）、Y. E. Engbers（利物浦大学/荷兰TNO）、D. Thallner（佛罗里达大学）、A. T. Clarke（利兹大学）和 R. K. Bono（佛罗里达州立大学）。

学术背景

这项研究属于地球物理学领域，聚焦于地球发电机的动力学、古地磁学以及地幔-地核相互作用。其核心科学问题是：地球最下部地幔的结构不均匀性如何影响并塑造了地球长期磁场行为。

地球磁场由外核液态铁合金的流动（即地球发电机）产生。从地核流向地幔的热流是驱动发电机的主要能源。地核-地幔边界是地核和地幔这两个地球内部关键动力系统交互的关键界面。地震层析成像揭示，最下部地幔（特别是约2800公里深度）存在显著的不均匀性，其特征是：两个巨大的低速省在赤道附近呈对跖分布（一个在非洲下方，一个在太平洋下方），它们被一个从北极延伸到南极的、速度较快的“环带”所分隔。这些结构（LLVPs）被认为比其他区域更热，因此导致了地核-地幔热流的巨大横向差异。

古地磁记录是探索过去磁场行为的直接档案，理论上可以反演地核-地幔边界两侧的结构和动力学信息。然而，由于古地磁记录时空分辨率的限制以及解释的不确定性，从这些记录中提取关于地幔结构影响的明确信号一直非常困难。前人研究要么关注短于10万年的时间尺度，要么试图用核幔不均匀性解释数亿年尺度的磁场行为变化，但都未能提供强有力的、跨时间尺度的证据。

因此，本研究的主要目标是：利用古地磁数据、古地磁场模型和数值发电机模拟，寻找并确认地幔热流不均匀性在地球古代磁场行为中留下的独特“指纹”，从而为核幔边界的热性质和长期演化提供新的约束，并增进我们对地球发电机稳定性和古地理重建的理解。

详细工作流程

本研究是一个综合性的交叉研究，结合了古地磁观测分析、古地磁场建模和地球发电机数值模拟三大板块，工作流程环环相扣。

第一流程：古地磁观测特征提取与分析 研究人员首先从过去2300万年的四个时间段（15-70千年、0-100千年、0-5百万年、5-23百万年）收集了高质量的古地磁场模型和火山岩数据集。他们聚焦于两个关键的古地磁观测指标： 1. 古地磁场长期变化（Paleosecular Variation, PSV）：通过虚拟地磁极（Virtual Geomagnetic Poles， VGPs）的角离散度（S）来量化。研究特别关注低磁纬度（南北纬30度以内）的S值，因为其平均值（S|30|_median）反映了瞬时轴向偶极子的主导程度，而其四分位距（S|30|_iqr）则反映了PSV在经度方向的变化幅度。 2. 时间平均场（Time-Averaged Field, TAF）：对长时间序列的磁场球谐系数取平均得到。除了计算轴向偶极子与非偶极子的能量比（AD/NAD_TAF）外，本研究创新性地量化了非偶极子部分中“带状项”（仅随纬度变化）与总非偶极子的能量比（ZND/ND_TAF），以评估时间平均场在经度方向上的对称性。

通过对这些模型和数据的分析，研究人员提取了四个关键的观测约束条件作为后续与数值模拟对比的“标尺”： * C1 (S|30|_median)：低纬度VGP角离散度的中值，反映平均偶极子强度。 * C2 (S|30|_iqr)：低纬度S值的四分位距，反映PSV的经向变化程度。 * C3 (AD/NAD_TAF)：时间平均场中轴向偶极子的主导性。 * C4 (ZND/ND_TAF)：时间平均场非偶极部分中“带状对称”的程度。

第二流程：地球发电机数值模拟 为了探究何种条件下能够重现上述观测特征，研究团队进行了31组（大部分是新运行的）地球发电机数值模拟。这些模拟分为六个系列，核心控制变量是热瑞利数（Ra），最关键的自变量是施加在外边界（模拟核幔边界）的热流模式。 * 均匀边界模拟（对照组）：两个系列（9个模拟），在核幔边界施加均匀的热流（异质性参数 q* = 0）。 * 非均匀边界模拟（实验组）：四个系列（22个模拟），在核幔边界施加强烈不均匀的热流模式。该模式基于地震层析成像模型，在非洲和太平洋下方设置了两个对跖的、赤道附近的低热流异常区（对应LLVPs）。异质性强度q*设定为2.3或5.0，q*越大表示热流峰值与谷值的差异相对于平均值越大。其中一系列模拟还考虑了热化学双扩散驱动，以更接近地球外核的实际情况。

模拟的参数设置（如埃克曼数、普朗特数、磁普朗特数）覆盖了从计算可行到更接近地球真实条件（快速旋转、强磁场）的不同范围，以检验结果的鲁棒性。每个模拟都运行了相当于数万至数十万年的物理时间，并输出磁场球谐系数的时间序列。

第三流程：模拟结果分析与观测对比 对于每个数值模拟的输出，研究人员采用了与处理古地磁数据完全相同的分析方法：计算其合成PSV（S值及其分布）和TAF（时间平均场结构）。然后将模拟输出的C1-C4四个参数与从古地磁观测中提取的约束进行直接、定量对比。此外，他们还对模拟结果进行了空间可视化，绘制了PSV随经度的变化图和TAF的全球分布图，并与地震层析图像进行定性比较。

第四流程：向更古老地质时期拓展 为了检验核幔不均匀性的影响是否长期存在，研究人员将分析拓展至更古老的古地磁数据集（涵盖5-2.65亿年，甚至到27.1-29亿年）。尽管这些更古老的数据无法提供可靠的TAF模型，但它们能提供PSV信息（C1和C2）。通过将这些古老数据的PSV特征与数值模拟的结果进行对比，可以间接推断地幔不均匀性影响的时间尺度。

第五流程：机制探讨与附加分析 为了理解非均匀热流影响磁场的物理机制，研究团队对部分模拟进行了额外的诊断分析，例如探讨“磁场屏蔽效应”——即低热流区（LLVPs下方）可能抑制浅层地核的对流和磁感应，从而影响磁场结构。他们还进行了数据重采样实验，以评估古地磁数据空间分布稀疏性是否会导致观测到的PSV经向变化被夸大或模糊，从而验证结论的可靠性。

主要结果

1. 从近期古地磁记录中提取的关键特征： * 过去2300万年的古地磁场模型显示，低纬度的PSV（S值）存在明显的经度调制，呈现近似球谐二阶（即两个波峰和两个波谷）的“双驼峰”模式，其波长与最下部地幔地震不均匀性的主导波长惊人地相似。 * 尽管由于数据稀疏性，单个火山岩数据集的经度信号不清晰，但所有近期数据集（0-5 Ma）的S|30|_iqr值（C2）都显著不为零，表明PSV存在经向变化。 * 所有时间平均场模型（0-23 Ma）在去除偶极子后，都显示出显著的非带状（即随经度变化）结构（C4值在0.19-0.57之间），打破了轴向对称性。这与2015年现代地磁场（C4=0.23）的特征一致，但时间平均场的偶极子主导性（C3）远高于现代瞬时场。

2. 数值模拟与观测的对比结果： 这是本研究最核心的发现。对比清晰地显示，只有那些施加了强烈非均匀外边界热流（q* > 0）的发电机模拟，才能同时满足所有四个古地磁观测约束。 * 满足C1和C3：均匀和非均匀模拟在调整热瑞利数（Ra）后，都可能产生符合C1（S|30|_median）和C3（AD/NAD_TAF）的磁场。但在高Ra值时，均匀模拟容易突然跃迁到多极子主导的不稳定状态（不符合地球大部分时期的情况），而非均匀模拟则能阻止这种跃迁，维持偶极子主导的稳定状态。 * 无法满足C2和C4：均匀模拟完全无法同时满足C2（S|30|_iqr）和C4（ZND/ND_TAF）。它们产生的PSV几乎是轴对称的（C2极低），其时间平均场的非偶极部分也几乎完全由轴向八极子等带状项主导（C4接近1），严重缺乏经向结构。这与观测到的“经度依赖的PSV”和“非对称的TAF”特征完全不符。 * 非均匀模拟的成功再现：相反，非均匀模拟自然地产生了PSV的经向调制（C2值与观测相符），其TAF也包含了显著的、与观测特征相似的非带状结构（C4值与观测范围重叠）。热化学模拟的结果同样支持这一结论，表明其具有鲁棒性。

3. 对古老记录的分析结果： 将涵盖过去2.65亿年的多个古地磁PSV数据集汇总分析后发现，它们的C1和C2参数值都落在数值模拟参数空间中唯独由非均匀边界模拟所占据的区域。而均匀边界模拟（包括本研究及前人研究）的输出结果，在同时满足这两个PSV约束方面与古老数据存在系统性偏离。这强有力地暗示，地幔热流不均匀性对地球发电机的影响至少已持续了数亿年。更古老（前寒武纪）的数据由于不确定性较大，结论的确定性稍弱。

4. 机制分析结果： 附加的“磁场屏蔽”实验表明，非均匀热流的影响是双重的一方面，低热流区抑制了其下方地核最顶部的对流和磁场，产生了直接的屏蔽效应，这有助于形成PSV的经向结构。另一方面，这种热异质性还通过改变带状流等更深层的发电机过程，影响了磁场的整体稳定性。边界热异质性导致的“赤道热、两极冷”的横向温度变化，可以驱动加强向西的带状流，而这被认为与稳定的偶极子状态相关。

结论与意义

本研究得出结论：地球最下部地幔存在的强烈横向热流不均匀性（与今天地震学观测到的大型低速省结构相关），在地球古代磁场行为中留下了清晰、可识别的“指纹”。这些指纹包括古地磁场长期变化的经度依赖性，以及时间平均场结构的轴向对称性破缺。这些特征至少在过去数亿年间持续存在。

科学价值： 1. 提供了关键证据：首次通过系统性地整合古地磁观测、场模型和发电机模拟，为“地幔结构显著影响地球发电机”这一长期假说提供了强有力的、多时间尺度的证据。 2. 建立了新的约束工具：研究所识别的古地磁特征（C2和C4）可作为新的工具，来约束地核-地幔边界的热性质和其随时间演化，连接了地震学、地磁学和地球动力学。 3. 深化了对发电机稳定性的理解：研究表明，地幔的强烈热异质性可能通过维持稳定的带状流等方式，帮助地球发电机保持长期稳定的偶极子主导状态，避免向多极子态跃迁。 4. 对古地理学的影响：研究指出，由核幔不均匀性产生的持久性非偶极场结构，可能导致古地磁偏角和倾角出现超过10度的、随经度和纬度变化的异常。这为厘清一些长期存在的古地理重建争议（例如某些快速极移或构造解释的争论）提供了新的物理视角和潜在的解释机制。

研究亮点

1. 多学科交叉与数据整合的创新性：将高精度古地磁场模型、大规模火山岩数据汇编和先进的地球发电机数值模拟三者深度结合，构成了一个完整的“观测-模拟-验证”研究闭环，方法学上具有很强的创新性和说服力。
2. 提出了新颖的量化指标：首次明确地将PSV的经向变化幅度（S|30|_iqr）和TAF的非带状对称程度（ZND/ND_TAF）作为核心诊断指标，用于量化磁场受地幔影响的特征，并证明它们是区分均匀与非均匀核幔边界条件的关键。
3. 结论的鲁棒性：通过在不同参数（包括接近地球条件的极端参数）和不同驱动机制（纯热驱动、热化学驱动）下进行大量模拟，证实了“q* > 0”是满足古地磁约束的必要条件，使得结论不依赖于某个特定模拟，具有很高的鲁棒性。
4. 拓展了研究的时间维度：不仅聚焦于新生代，还将分析有效地拓展到中生代甚至更早，为核幔结构的长期稳定性及其对地球历史的持续影响提供了间接证据。

其他有价值的内容

研究还展望了未来方向：更大、更高质量的古地磁数据集将减少不确定性；在更极端参数下进行更长时间的数值模拟将进一步验证结论；地震学可以探测地核最顶部的结构，以检验是否存在地幔不均匀性导致的持久性结构。这些都将使对地球深部相互作用的理解更加精确和深入。