地球内部的流动，怎样维持大尺度磁场并改变其极性？

一直存在的磁场，需要一直做什么？

指南针能在地表许多地方给出大致一致的南北方向，说明地球磁场具有跨越行星尺度的组织。但这并不意味着地底藏着一根永久不变的磁棒。外核是能导电的流体，磁场在其中不断被流动改变，又不断因有限电导而损失能量。要让它长期存在，补充必须持续发生。1

这里至少有三个不同的问题。磁场怎样不衰减掉，是供能和再生的问题；为什么大尺度上常由偶极分量主导，是空间组织的问题；为什么这种主导分量会改换符号，是状态变化的问题。 一个模型能回答第一问，不代表已经回答后两问。

“偶极”可以先理解成一种最简单的大尺度场形状，像有两个主要磁极，但不要求真的存在两个孤立磁荷。真实地磁场还包含倾斜、四极和更多细节。某个测站的方向改变、偶极轴暂时偏离、轴向偶极改变符号，也不是完全相同的事件。古地磁记录需要把局部测量与全球场形状联系起来，才能讨论反转。

本文先讲清导电流动如何持续生成磁场，再用一个能实际计算的例子区分“维持”与“换极性”。随后比较液态金属实验、古地磁数据，以及2025—2026年的三维发电机研究。需要保留的一条边界是：这些证据已能约束机制，但还没有选出适用于全部地球历史、且能预测下一次反转的唯一模型。

导电只是条件，能把磁场再生出来才是发电机

想象外核中已经有一点弱磁场。一团导电流体从它旁边流过，不会仅仅像搬运一块带颜色的水那样把磁场原封不动搬走。运动、感应电场和电流要共同满足电磁规律；当不同位置的速度不同，磁场的方向与强度都可能被改变。

在电导率近似均匀、忽略位移电流等通常的磁流体条件下，磁感应方程写成：

这里 $\mathbf B$ 是磁场，$\mathbf u$ 是流速，$\sigma$ 是电导率，$\eta_m$ 是磁扩散率。第一项描述流动怎样改变磁场，第二项描述有限电导造成的磁扩散。它不是“导电液体一动就自然发电”的保证书。1

对不可压缩流动，还可以把同一方程改写为：

左边是跟着流体观察到的变化。右边第一项强调沿磁场方向的速度差：相邻部分受到不同的运动，才可能拉伸、转向和重新组织磁场。整团流体只做刚性平移，不能靠平移本身持续放大磁场；刚体旋转也没有持续拉伸的应变。把“地球会转”直接接成“所以磁场一直存在”，中间漏掉了最重要的流动结构。

真正的难点是闭合再生过程。一次剪切可以把一种磁场分量变成另一种；如果新的分量随后不能通过别的流动环节补回原来的分量，最初的供给最终仍会耗尽。旋转对流模型中的绕转、沿柱方向运动和回流，能使这种分量转换构成循环。具体由哪些环节占主导会随参数改变，不能把某张理想化流线图当成所有地磁发电机的唯一内部结构。1

常用的磁雷诺数把流动与磁扩散作量级比较：

$Rm$ 较大，表示在这一尺度上，流动有时间显著改造磁场，扩散还来不及完全抹平它。可是有改造能力不等于改造得对：流动也可能把场推向更细的结构，使损耗增加，或生成彼此抵消的分量。自持发电机要检验的是，完整感应过程能否在损耗之上保持净增长，并最终进入有界状态。

这也解释了“种子场”的地位。上面的方程在严格 $\mathbf B=0$ 时不会凭空生成非零解；发电机讨论的是弱磁场如何被放大和维持。种子提供起点，长期能量仍来自流动所获得的机械功，不能把最初的一点磁场当作永久能源。

能量来自哪里，旋转和磁场各做什么？

对目前研究最多的对流发电机，驱动力来自外核中的浮力：热状态和成分差异让不同位置的流体密度不同。内核凝固释放潜热，也把部分轻元素留在周围液体中；这两类过程可以推动外核运动。更长期的冷却、内核成长和核幔边界散热共同限定可用能量。没有内核的早期地球需要另外讨论，不能把今天模型中的底部供给原封不动倒推到全部历史。127

旋转主要组织流动。 科里奥利力使流动倾向于沿旋转轴形成关联结构，从而改变再生磁场的几何条件。它对瞬时速度本身不做功，因为 $\mathbf u\cdot(\boldsymbol\Omega\times\mathbf u)=0$。这句话限定于科里奥利项；不等于一概否认潮汐、进动等机械驱动研究可能调用旋转能量。本文的三维证据以对流驱动为主。

磁场也会改造自己的生成者。 流体中的电流受到 $\mathbf J\times\mathbf B$ 的洛伦兹力，流动因而减速、转向或重新分配。开始时可以把磁场当成很弱的被动量，先看给定流动会不会放大它；但一旦场足够强，就必须让这个反作用回到运动方程里。

因此，一个自洽模型不能只写“流动拉伸磁场”：它还要说明磁场长大以后，为什么不会无限长大，或者把整个再生环节彻底破坏。达到统计稳态时，能量可以持续从浮力驱动流入，并以电阻和黏性耗散流出；宏观磁场可以一直存在，而内部能量一直更新。这里的“稳态”通常是长期统计稳定，不是每个位置每一秒都不动。1

更微妙的是，这套过程通常没有指定唯一的磁场符号。把整个磁场改成 $-\mathbf B$，电流也改成 $-\mathbf J$，而洛伦兹力保持不变：

感应方程也保留相应的正负对称。于是，同一套流动与供能条件可以容纳两种相反极性。方程允许两个方向，并不保证系统一定会在它们之间转换；是否有可达路径，还取决于模式耦合、反馈和扰动。

三个变量，怎样让“有磁场”与“换极性”分开？

Gissinger在2012年研究了一个三模式模型。它把复杂场压缩为偶极幅度 $D$、四极幅度 $Q$，以及一个使两种磁模式相互耦合的流动幅度 $V$：5

这里刻意沿用原文的符号：线性增长项 $\mu Q$ 写在 $Q$ 方程中，不能为了让“偶极主导”的故事更顺，就擅自把它换到 $D$ 上。$\nu$ 在此只是一个模式阻尼系数，不是前面物理流体的黏度。

三条式子可以逐个理解。$Q$ 有自身的净增长，也会通过流动与 $D$ 交换；$D$ 自身会衰减，但可以得到 $VQ$ 的补充；$V$ 受到外部驱动 $\Gamma$，同时受到线性阻尼和磁分量乘积 $QD$ 的反馈。最后一项借用了洛伦兹力对磁场为二次的结构。

这已经把一段反馈链走通了：磁场模式改变流动，流动再改变磁场模式。可是 $\mu$ 等系数本来就把更复杂的感应过程压缩进去了，三个方程不从头证明真实外核为何能发电，也不是2026年三维论文直接缩减出来的同一组方程。 它用于看一种模式耦合如何产生极性变化。

本库采用原文中的两组参数，固定 $\nu=0.1$、$\Gamma=0.9$，初态均为 $(Q,D,V)=(1,-1,0.1)$。积分到模型时间6000，丢弃前1000的过渡段：

图1。本库按原文方程计算。横轴是任意模型时间，幅度也没有地球磁场的单位。原文将这一区域的变化解释为混沌吸引子的连接变化；这里展示实际轨迹，不用一张曲线代替混沌的一般证明。

把数值相对容差从 $10^{-10}$ 收紧到 $10^{-11}$ 后，第一组在同一长观察窗内得到32次过零，第二组仍为0次。短时间0—30内，DOP853与Radau两种求解器的最大幅度差约为 $6.2\times10^{-10}$。长时间的具体轨迹和次数没有逐点收敛，本文也没有证明其长期反转率统计已经收敛。因此可以使用“这两组条件呈现不同的极性行为”这个例子，不能把15次或32次换算成地球的反转频率。

反转本身也不是“磁场整体瞬间关机再开机”。第一组基准计算中，所有观察窗内的 $D$ 过零事件，都有 $|Q|>2.29$。在所选两模式表示中，偶极暂时为零时，另一个磁分量仍然存在：

图2。左图是同一次转换附近的两个磁分量，右图是它们共同描出的轨迹。它没有经过“两种磁分量同时为零”的原点。这是模型中的形状变化，不是对某一次真实地磁反转做了全球成像。

这个例子的价值，是让三个容易混在一起的判断分开：有持续供能，不代表固定极性；波动很复杂，不代表一定反转；偶极过零，也不代表所有磁场都消失。

液态金属实验，给这个想法增加了什么依据？

VKS实验在圆柱容器里用两只反向旋转的叶轮驱动液态钠。叶轮转速相同与不同时，装置的某种空间对称性不同，观察到的磁场状态也不同。在所研究的条件范围内，等速状态产生统计上固定极性的场；破坏对称性后，某些转速区间出现周期振荡或不规则反转。6

2010年的分析并非只拿某个探头的一条变号曲线，声称它就是完整磁场反转。研究者利用关于装置中面对称的两个测点，把信号组合成两类不同对称性的磁分量，再考察它们怎样共同变化。结果支持偶极类与四极类分量在反转过程中发生明显的配合：偶极减弱和换号时，另一类分量不能忽略。

这里的“偶极类”“四极类”还包含具有相同对称性的高阶成分，并不是把整个实验场精确截成球谐的一阶和二阶。实验几何、叶轮驱动和边界也不同于地球。它提供的是导电流体中确实可以出现可调控的自持磁场与模式转换，以及对称性干预具有区分力，不是证明真实外核一定按同一个低维机制反转。

对地球来说，我们无法照样拧动一个旋钮。需要从空间和时间结构中寻找约束，不能只比较“都能反转”这一条特征。

古地磁怎样检验地幔是否在影响外核？

地幔运动远比外核流动慢，但两者通过核幔边界交换热量。如果不同位置散热不同，就可能让外核某些区域更容易对流，另一些区域受到抑制。这是一个很具体的边界影响：地幔不必直接推动所有外核流体，改变热量从哪里流出，就可能改变内部组织。2

Biggin等2026年的研究没有只问“模拟里有没有一根大偶极”。它把古地磁数据和模型输出放到四类统计约束下比较。先理解其中最直观的一种。

一块岩石记录了形成或获得磁化时的局部方向。按偶极几何把这个方向反算回一个等效磁极位置，可以得到“虚地磁极”。同一研究地点、多次形成的岩石给出的这些位置有多分散，用角离散 $S$ 表示。它是磁方向变化的统计量，不是外核流速，也不是磁极真实运动的逐帧录像。

下面重绘作者公开的两个时间段数据：0—5百万年前40个研究地点，以及5—23百万年前27个研究地点。浅绿色区域标出研究重点使用的正负30度磁纬度窗口。3

图3。每个点代表一个研究地点汇总出的离散度，不是一块岩石。误差条直接采用公开源数据中的95%区间长度，本库没有重新处理原始试样或重做作者的一万次重采样。

低磁纬度窗口中分别有19和12个地点。用作者的分位数口径计算，$S$ 的中位数约为11.86度和15.37度，四分位距约为4.32度和6.18度。8 两组点的空间分布、误差和数量不同，不能仅凭后者数值较大，就说“更早的地球一定反转得更频繁”。原研究在计算这些指标前还剔除了与游移、反转有关的异常虚地磁极，这些数字本来就不是反转次数的估计。

四类指标分别约束不同方面：

第三行尤其值得停一下。假设大偶极一直为1，而某个非偶极分量在 $+a$ 和 $-a$ 间波动：先平均分量，会得到零；先平方再平均，会得到 $a^2$。所以“长期平均场很规则”可以与“每一时刻都有明显复杂结构”同时成立。这个简单代数例子没有重建地球数据，只说明统计顺序为什么重要。

原研究计算平均场前，已经按轴向偶极的符号把每个时刻的全部磁场系数统一翻转到同一种极性；因此不能误解成作者把正常、反向场直接相加，靠正负抵消得到了结论。其变化统计保留到球谐10阶，平均场比较保留到4阶，观测能约束的空间细节本来就不同。2

与模型相比，真正增加区分力的是哪一项？

研究比较了31组三维模拟，分成6个系列。均匀散热与不均匀散热模型都可能产生大偶极；仅凭平均偶极占优，区分力有限。较有区分力的是：在维持偶极主导时，均匀边界系列的低纬地点差异往往太小，平均非偶极场又过于轴对称；引入地幔层析成像启发的散热差异后，能产生更明显的经度结构。

以下是补充表1中两对相同驱动参数的例子，均采用 $E=3\times10^{-4}$、$Pr=1$、$Pm=3$。只在每对内比较热流边界，$Ra_T$ 是该论文定义下的热驱动参数：

这里 $q^*=(q_{\max}-q_{\min})/\bar q$。时间平均非偶极场的观测约束约为0.19—0.57，因而表中的不均匀边界例子更接近这个范围。表中是作者模拟的已发表输出，不是本库重新跑了三维磁流体方程；各模型运行时长也不同，不能把四行当成等长独立实验。2

这支持的判断应当是：在已测试的模型家族中，加入空间散热差异有助于同时解释更丰富的地磁统计。 它不是数学定理意义上排除了全部均匀边界发电机，也没有穷尽参数来同时满足所有观测指标。

一个直接的竞争解释是采样不均匀：均匀边界产生的场，也可能因为只测了少数地点而显得各处不同。作者按实际地点分布，对每个相关模拟做了12次随机降采样；模型家族的区分有所模糊，但均匀边界仍未同时满足所用的两项方向离散约束。这使结论比单纯“看图很像”更有依据，不过本库没有独立重跑该检验，也不能把12次降采样当成排除了所有记录偏差。

另外，从地震波速度结构到温度、再到热流边界，还需要物理解释。研究假设低速大区域相对偏热，即使其地震性质也可能包含成分影响。近期古地磁模型中的经度图样提供了关联线索，但0—5和5—23百万年前的数据，直接定位同一种经度图样的检验仍因误差过大而未决。更久远数据支持长期存在不对称性的解释，不能升级为“今天每个地幔大区域的位置已经被磁记录证明数亿年没动”。

稳定分层为什么既能保护大尺度结构，又可能改变反转？

“稳定分层”说的是把一小团流体上下移动时，浮力倾向于让它返回，而不是继续翻滚。它限制径向对流，不等于整个区域完全静止，更不等于磁场极性被固定。

Müller等2026年的研究，在旋转球壳中直接求解不可压缩磁流体方程，人为指定外核顶部一层较稳定的温度梯度。模型采用半径比0.35、无滑移机械边界、导电内核和绝缘地幔；稳定层厚度从0变化到流体壳厚度的0.24。主要参数为 $E=10^{-3}$、$Pr=1$、$Pm=10$。4

首先，稳定导电层会过滤一部分快速变化的磁结构。可以用磁扩散的集肤效应理解：对单一频率的振荡，在均匀导体的简化情形下，穿透深度约为 $\delta_{\rm skin}=\sqrt{2\eta_m/\omega}$。频率越高，穿透越浅。完整球壳还涉及空间结构和流动；这条式子只帮助理解为什么深处的短时变化不一定原样出现在上边界。

在该论文的参数系列里，加入稳定层使边界处的偶极成分更突出，偶极主导向多极主导的转变移到更强驱动处。但下一步加上热流的空间差异时，结果并不只有一种：

这些例子共同使用层厚0.24、$Ra/Ra_c=25$、$N_{\max}/\Omega\approx10$；幅度按该论文 $\delta q=(q_{\max}-q_{\min})/(2\bar q)$ 定义。它的分母比前一篇 $q^*$ 多一个2，不能直接拿“5.3%”与“2.3”比较强弱。$Y_{10}$ 的“没有反转”也是有限运行范围内的结果，不是这种图样永远不能反转的定理。

为了理解原因，作者又做了运动学发电机计算：去掉磁场对流动的洛伦兹力，但仍让流动自己随时间演化，再看磁场怎样响应。它不是固定一张速度快照。在较低驱动、赤道对称、两类模式尚未耦合的对照条件下，稳定层变厚，会使偶极与四极模式的增长率更接近；加入对称性破缺后，两者更容易显著相互作用。这给“既保留较强偶极，又能换极性”提供了一种机制解释，但增长率对照与上表高驱动的自洽反转，不是完全相同的运行条件。

但这不等于量出了真实地球的模式增长率。上述 $E$ 和 $Pm$ 与地球参数仍差得很远，稳定层的存在、强度、厚度也不是本次模拟测出来的。论文中把此机制推广到低惯性地球的说法，还需要进一步检验。

另一个前沿解释：来不及补充的磁通量

不能只因模式图看起来通顺，就把其他过程排除。Aubert等2025年提出了一个不同但可以联系起来的解释：近核表面的上涌与水平环流之间，存在补充和搬运的竞争。7

在其模型中，上涌有助于在近赤道区域形成支持既有偶极方向的磁通量结构，表面环流再把它们向极区输送。如果上涌被顶部稳定性明显削弱，水平搬运却仍然较强，磁通量可能在充分积累以前就被搬走。平均偶极因而变弱；如果波动没有同比减弱，相对于平均值的波动就变大，更容易出现极性转换。

这套解释不用预先把“惯性力增强到足以打乱旋转组织”当作反转的必要条件。作者在底部成分浮力驱动的模型系列中改变顶部稳定性，并分解偶极变化中的感应与扩散贡献，得到支持这种近表面运动学机制的结果。其 $Rm$ 约为950—2100，且采用不同的机械边界和浮力分布，与前一篇强稳定层模型不是同一个参数实验。

它还有一个值得保留的限制：用全球平均散热来调反转率时，能得到中间反转率的范围可能很窄；改用横向散热差异可在所测系列里缓解这个问题，但尚不能据此重建地球全部热历史。作者也指出，内核形成以前的反转仍需要另行研究。

模式增长率接近，与近表面补充/搬运失衡，未必是互相排斥的两个“终极答案”。 前者描述模式怎样相互转换，后者定位某些流动过程怎样改变偶极的收支；它们可能在同一完整系统中发生联系。但目前本文引用的研究没有把两套诊断在同一批地球参数模型上全部做完。把它们直接宣布为同一个机制，或仅凭图样相似宣布哪一个获胜，都超过现有比较。

哪些证据能继续缩小答案？

这一问题的难点，不是模型里完全做不出反转，而是很多不同模型都能做出某种反转。因此比较应该逐渐从“有没有变号”走向同时解释场强、空间结构、短时变化与长时历史，并检查改变条件时的响应。

目前较稳的是：外核导电流动、持续供能与磁场反作用构成自持发电机的物理基础；单靠导电性、整体旋转或一次放大，都不足以说明完整过程。液态金属实验和磁流体计算表明，维持大尺度场与锁定极性可以分开。

有支持、但依赖条件的是：核幔散热差异能改变长期场结构，顶部稳定性能改变过滤、模式耦合或近表面磁通量输运。不同模型出现不同反转响应，要求我们比较具体的浮力分布、流动和参数，不能只看“稳定”“更热”“更不均匀”这些标签。

仍未确定的是：真实地球在某次反转中，各过程的相对贡献到底多大；这些贡献怎样随内核成长和地幔演化改变。本文没有建立下一次反转的预报模型，任何把一次局部减弱或某个低维曲线直接转换成倒计时的做法，都缺少从观测到全球状态演化的关键连接。

计算与资料回查

计算结果与数值检查、两组完整观察窗轨迹、古地磁重绘点与原误差长度随在线版提供。复算代码在仓库 research/calculations/geodynamo_polarity.py。

计算检查包括磁场正负对称、非零定点残差、短时求解器比较和一个二次量的变化恒等式。它们能核对这组方程的实现，不能证明真实外核的参数或机制已经被确定。三维论文的流场、全球磁场反演及原始岩石处理均未由本库独立重跑；来源记录注明了实际读取版本与范围。