论文信息：Zhang, Y., Sun, S., Mooney, W. D., & Xu, Y. (2026). Lithospheric magnetization derived from petrological and satellite constraints. Journal of Geophysical Research: Solid Earth, 131, e2025JB032111. https://doi.org/10.1029/2025JB032111

岩石圈磁化强度是理解地球磁场的岩石圈分量、反演地下物质组成与热状态的关键参数。然而，由于观测手段的固有局限，我们长期只能"看到"它的一小部分。本文以一作张壹带领的研究为代表，提出了一种岩石学先验 + 卫星磁测数据联合反演的新框架，重建了覆盖全波段的全球岩石圈感应磁化强度（以垂直积分磁化率 VIS 表征）分布。下面围绕这篇 JGR: Solid Earth 论文，梳理研究动机、方法要点与主要结论。

为什么长波长磁化结构一直难以获得

卫星磁测（CHAMP、Swarm 等）为我们提供了全球尺度的岩石圈磁场观测，如 LCS-1、CHAOS、MF7、CM6 等球谐（SH）模型。但这些模型的功率谱（Lowes, 1974）清楚地表明：岩石圈磁场只有在球谐阶数 16 阶及以上才主导观测量；在更低的阶数（长波长部分），信号被强大得多的地核磁场完全淹没（图1）。这意味着，仅凭卫星数据，我们无法可靠地分离出长波长的岩石圈磁场。

各磁场模型功率谱
图1　各磁场模型的功率谱（SH 1–80 阶，400 km 高度，引自论文 Figure 1）

岩石圈磁场的信号在 SH 16 阶及以上才占主导，更低阶的信号被地核磁场掩盖。这决定了单纯卫星反演无法恢复长波长磁化结构。

此外，还存在一类被称为磁湮灭子（magnetic annihilator）的问题——某些磁化强度分布（例如以固定磁化率磁化的球壳）不产生任何可观测磁场。它使得岩石圈磁化的全球平均值及其长波长特征本质上无法仅靠磁场观测确定。因此，非磁学手段（即地质—岩石学约束）对于确定全球平均磁化强度和长波长结构是不可或缺的。

方法：如何用"岩石学 + 卫星"补齐波长缺口

针对上述问题，本文采用球谐域反演策略，把长波长与短波长分别交给两类不同的约束：

  • 长波长（SH 0–16 阶）：由岩石学先验模型 SM3-SI（Hemant & Maus, 2005）提供约束。该模型基于世界地质构造图、主要岩石类型的实验室磁化率测量以及地壳地震厚度建立，可以给出可靠的全球平均值与长波长磁化变化，从而规避磁湮灭子问题。
  • 短波长（SH 16–80 阶）：由卫星磁数据（采用 CHAOS-8 模型）约束，卫星数据在这些阶数上分辨率高、可信度强。

反演在单个球面等效源层上进行（球面三角剖分生成约 81,920 个三角棱柱单元，垂向积分磁化率 VIS = 磁化率 × 层厚）。目标函数由数据拟合项与模型约束项之和构成，采用 Lévy-Gradient Descent（L-GD） 随机搜索算法迭代求解——该算法不仅能收敛，还能顺带给出反演模型的不确定度估计。合成数据实验表明，恢复的磁场在各波段误差均小于 1 nT，验证了反演框架的有效性。

反演验证图
图2　反演结果的定量验证（引自论文 Figure 4）

(a)(b) 参考与反演长波长 VIS 模型的 SH 系数差与功率谱（0–16 阶）；(c)(d) 观测与反演 Br 异常的 SH 系数差与功率谱（16–80 阶）。差异主要集中于 15–16 阶过渡带，其余阶数拟合良好。

最终获得的 VIS 模型给出了非负、高分辨率的全球岩石圈感应磁化分布——负值仅在约 2.5% 的极少数区域（主要位于海洋）出现，反演 200 次迭代后达到 RMS = 2.87，表明岩石圈感应磁化与观测磁场在 16–80 阶高度自洽。

主要结果：磁化强度如何在构造单元间变化

新模型的最直观产出是图3所示的全球 VIS 分布。它清晰地刻画了磁化强度与大地构造及主要地质省的对应关系：克拉通区域普遍呈现高 VIS，而活动造山带区域显著偏低。

全球VIS模型
图3　反演得到的全球岩石圈感应磁化强度（VIS）分布（引自论文 Figure 6，为论文核心结果）

黑色为大陆与板块边界，红色为板块界线。前寒武纪克拉通（北美、东南欧、西伯利亚、澳洲等）呈现高磁化，造山带与伸展地壳偏低，大洋高原与部分俯冲带亦表现出较强磁化。

对反演模型按地壳类型与基底年龄统计（图4），可以得到定量结论：

  • 地壳类型：盾、平台、盆地（前寒武纪区域）的平均 VIS 高达约 0.8 km·SI；造山带居中（约 0.48 km·SI），伸展地壳与大陆边缘最低（约 0.46 km·SI）；洋壳与洋中脊最低（约 0.4 km·SI），大洋高原（LIPs）可达约 0.61 km·SI。
  • 基底年龄：平均 VIS 从太古宙（0.773 km·SI）向更年轻单元单调递减，至中生代—新生代降至约 0.477 km·SI。
  • 值得注意的反差：岛弧、大陆弧与俯冲带的磁化（约 0.49–0.57 km·SI）显著高于纯岩石学模型 SM3-SI 的预测（此前这些区域的强磁化在岩石学模型中基本缺失），暗示俯冲系统内存在超越地壳的磁源。
不同构造与年龄平均VIS对比
图4　不同地壳类型（左）与基底年龄（右）的面积加权平均 VIS 对比（引自论文 Figure 8）

反演模型（inverted）与岩石学模型 SM3-SI（含其长波长分量）对比。两者在大陆地壳上总体一致，但在岛弧、俯冲带等区域存在显著差异。

利用反演的 VIS 模型可以正演重建完整的岩石圈磁场，包括此前卫星无法直接观测的长波长分量（图5）。重建的长波长磁异常与前寒武纪基底、地壳类型及近期构造活动呈现更密切的相关性，修正了强磁化区域的畸变异常。

重建磁场图
图5　重建的岩石圈磁场（引自论文 Figure 7）

(a) 包含剩磁的完整岩石圈 Br 分量；(b) 感应 Br 分量；(c) 长波长 Br 异常（SH 1–15 阶），400 km 高度计算。前寒武纪基底区域呈现清晰的大尺度正背景，澳洲则显示显著负异常。

对磁源成因的解释与讨论

结合全球热流模型（Lucazeau, 2019）与结晶地壳厚度 ECM1 模型（Mooney et al., 2023），作者对磁化强度的成因给出了物理解释：

  1. 前寒武纪省具有低温热流（<60 mW/m²）与厚地壳（>38 km）——居里等温面深，磁化更强；随着基底年龄变新，热流升高（54–64 mW/m²）、磁化率下降。
  2. 造山带与伸展地壳地壳减薄、磁性矿物减少、居里面抬升，故磁化偏低。
  3. 海洋区的高磁化几乎只与**大洋大水成岩省（LIPs / 洋底高原）**相关——其热流较低、地壳较厚。
  4. 俯冲带与岛弧的强磁化最引人注目：冷板片使居里深度加深（可超过 50 km），板片脱水在地幔楔内引发蛇纹石化，显著提高磁化率，从而形成超越地壳的深部磁源。

需指出的是，第 4 点推断基于磁异常、居里深度估计与蛇纹石化模型之间的间接关联，作者也强调需要借助地震学、大地电磁等其他地球物理手段进一步验证。

结语与展望

这项工作通过岩石学先验约束长波长 + 卫星数据约束短波长的联合反演策略，突破了卫星数据无法直接恢复长波长岩石圈磁场的瓶颈，获得了一张更高分辨率、更准确的全球岩石圈感应磁化分布，并厘清了不同地壳类型与基底年龄的磁化特征及其岩浆—热状态成因。

文章同时指出若干改进方向：当前模型将岩石圈视为单一磁层，未来可考虑上、下地壳分离为多个磁层；海洋剩磁模型的精确性也有待提升（例如考虑倾斜极性边界的剩磁结构）。作为开放获取论文，其数据与源码已托管于 Zenodo（Zhang, Sun, et al., 2025），便于后续研究与复现。

延伸阅读与引用

  • Hemant, K., & Maus, S. (2005). Geological modeling of the new CHAMP magnetic anomaly maps. JGR, 110(B12).
  • Olsen, N., Ravat, D., & Purucker, M. E. (2023). On the determination and interpretation of the lithospheric induced magnetisation. EPSL, 606.
  • Williams, S. E., & Gubbins, D. (2019). Origin of long‐wavelength magnetic anomalies at subduction zones. JGR: Solid Earth, 124(9).
  • Zhang, Y., Sun, S., Mooney, W., & Xu, Y. (2025). Lithospheric magnetization derived from petrological and satellite constraints [Dataset]. Zenodo. https://doi.org/10.5281/zenodo.15129305

本文为对该论文的中文解读，图片均引自论文原文。