岩石圈磁化强度的全球新图像：来自岩石学与卫星磁测数据的联合约束
论文信息：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，验证了反演框架的有效性。 ...