论文信息：Li, M., Xu, Y., Liu, L., Yang, B., Zhang, Y., Zhu, Y., & Liu, S. (2024). Physical link between effective viscosity and electrical resistivity for dislocation creep in upper mantle and its application in Northwest Xinjiang, China. Geophysical Journal International, 240, 1295–1307. https://doi.org/10.1093/gji/ggae438

上地幔的有效黏度是理解岩石圈动力学与板块运动的核心参数，但在野外观测中直接"测量"它却极其困难。好在近二十年大地电磁（MT）方法给出了一条诱人的路径：既然电阻率剖面可以刻画地下的温度、含水量与部分熔融，那么它能否进一步"换算"出黏度分布？过去的研究大多建立在"电阻率与黏度正相关"的经验假设上，其物理根基仍不够扎实。本文（第一作者为李蔓，浙江大学地球科学学院，我们课题组的合作成果之一）从点缺陷化学与电中性原理出发，为位错蠕变控制的上地幔建立了一条电阻率—有效黏度之间的微观物理桥梁，并成功应用于新疆西北部的塔里木—天山—准噶尔地区。下面梳理这篇 GJI 论文的核心思想与结果。

为什么要把黏度与电阻率联系起来

有效黏度的经验公式（Karato & Wu, 1993；Xu et al., 2018）形如

η = σ / ε̇ = A · dᵐ · Cwʳ · fO₂^q · σⁿ⁻¹ · exp[(E + P·V)/RT],

其中 σ 为应力、ε̇ 为应变率、d 为粒径、Cw 为含水量、E 与 V 分别为激活能与激活体积。而岩石的电阻率 ρ 具有几乎相同的数学形式（也是温度的指数函数，依赖含水量与氧逸度）。两者同形，正是 MT 数据"换算"黏度的诱人之处。

但关键在于为何二者会联系起来。物理根源在于：地幔岩石的蠕变由固体中最慢的粒子（硅）沿位错的扩散速率决定（爬升控制的位错蠕变），而电导则由带电载流子中最快的粒子（小极化子、质子等）的运动主导。一个决定变形、一个决定导电，看似风马牛不相及——然而固体中正负电荷的浓度必须处处满足电中性条件。由于静电相互作用能很大，硅空位浓度一旦改变，载流子浓度必然同步变化。正是这一电中性纽带，把两种性质截然不同的输运过程"绑定"在了一起。

方法：从微观点缺陷到宏观换算公式

作者先论证了上地幔中位错蠕变是主导变形机制（依据地幔捕虏体的组构定向、地震各向异性与热—力学数值模拟）。在爬升控制的位错蠕变下，应变率正比于硅空位浓度：[V_Si]（即 ε̇ ∝ [V_Si]）；而电导率则正比于载流子浓度 N（σe ∝ N）。只要建立 [V_Si] 与 N 之间的关系，就能把黏度与电阻率在宏观层面联系起来——这一步靠的是点缺陷化学与电中性条件。

按照含水量（更准确地说是水活度/氢缺陷）的高低，作者给出了两个理想模型：

  • 无水模型：主导缺陷为三价铁 Fe·M 与 M 位空位 V’‘M，电中性条件 [Fe·M] = 2[V’‘M]。此时变形受硅空位扩散限制，导电由小极化子（Fe·M）扩散主导，得到 [V_Si] ∝ [Fe·M]²，即指数 m = 2。
  • 含水模型：主导缺陷变为 Fe·M 与 H·M，电中性条件 [Fe·M] = [H·M]。变形受 (3H)’_Si 扩散限制，导电由水相关缺陷 H· 主导，得到 [(3H)’_Si] ∝ [H·]^(5/3)，即 m = 5/3。

代入公式并取恒定应变率，最终得到电阻率—有效黏度的换算关系（记为 η_eff/η0 = (ρ/ρ0)^(m/n)）。其中 n 为应力指数（位错蠕变一般 n ≈ 3），m 由表 1 按干/湿条件分别取 2 或 5/3，从而幂指数 m/n 分别约为 2/3 与 5/9。η0 与 ρ0 为参考黏度与参考电阻率（取自实验室实测，用 Karato & Jung (2003) 的蠕变参数和 Gardés et al. (2014) 的电导参数标定）。相比以往需要地球动力学建模或用经验系数标定的做法，这套方案更接近地幔的真实物化状态，物理图像也更清楚。

缺陷浓度随含水量示意
图1　橄榄石中各类点缺陷浓度随含水量的变化示意（引自论文 Figure 1）

随含水量变化，主导带电缺陷由低水时的三价铁 Fe·M 与 M 位空位 V_M，转为高水时的 Fe·M 与 H·M。这一电中性框架是连接蠕变与电导、进而换算黏度的关键。

主要结果：把电阻率断面"译"成黏度分布

为检验该模型的适用性，作者把新疆西北部塔里木—天山—准噶尔地区的最新三维电阻率结构（Liu et al., 2024）中沿剖面提取的电阻率断面换算为有效黏度分布。计算中温度取自东亚岩石圈三维热结构模型（Sun et al., 2022），应变率取自 GPS 观测（10⁻¹⁵ s⁻¹）。

有效黏度分布结果
图2　由电阻率断面换算得到的上地幔有效黏度（引自论文 Figure 3，为论文应用核心结果）

(a)–(c) 分别为无水、含水和组合模型换算得到的有效黏度剖面；(d)(e) 为取自 Deng et al. (2017) 的“软”“硬”岩石圈强度作为对比。图(a)、(d) 中五个三角标记为与图 Fig.5 逐点对比的位置。

主要定量结果包括：

  • 大尺度正相关：无论采用何种电导/蠕变机制组合，大尺度上有效黏度都随 MT 成像的电阻率升高而升高。
  • 横向差异：塔里木与准噶尔顶部岩石圈地幔的有效黏度整体高于天山地区；而天山下方的岩石圈底部黏度又略高于两大盆地。
  • 量级：无水模型（V’’_Si/Fe·M）给出的黏度范围为 10²⁰–10²²·⁵ Pa·s；含水模型（(3H)’_Si/H·）为 10¹⁹·⁷–10²²·⁸ Pa·s。
  • 机制转换：组合模型在深度上大致以 1300 K 等温面为界，浅部取含水、深部取无水机制——这与实验室观测到的橄榄石电导在约 1300 K 附近由质子导电切换为小极化子导电的规律相吻合。

讨论：稳健性、不确定性与物理边界

作者通过与独立数据对比，论证了结果的稳健性。与实验室外推黏度剖面（图4）相比，文中换算结果的量级与随深度递减趋势均接近；与 Deng & Tesauro (2016) 的"软/硬"岩石圈强度模型（图5）比对也能在 MT 剖面的五个点位处良好吻合，进一步确认了电阻率换算黏度的物理有效性。同时，作者系统考察了几类不确定来源：

  • 含水量影响：大陆岩石圈橄榄石含水约 60–100 ppm，而下部岩石圈几乎无水；在相同电阻率下，含水模型（(3H)’_Si/H·）会高估下部岩石圈的黏度，因此应用时应按深度分辨干湿状态。
  • 应变率影响：由 ∂(lg η_eff)/∂(lg ε̇) = (1−n)/n ≈ −2/3（n≈3）可知，应变率提高一个数量级，黏度约下降 0.67 个对数单位（约降至原来的 0.2 倍）。天山褶皱带应变率约 10⁻¹⁵ s⁻¹，而塔里木与准噶尔仅约 3.2×10⁻¹⁶ s⁻¹；考虑二者差异后，对黏度的扰动仍小于 0.5 个对数单位。
  • 温度模型影响：不同温度数据库（如 Deng 模型比 Sun 模型约高 200 K）对无水/含水模型结果的影响都小于 0.5 个对数单位，但会显著改变组合模型中干湿模型的分界深度。

作者也坦诚地界定了方法的适用边界：不能直接用于熔体连通域——当熔体体积分数超过约 0.5%（连通阈值）时，电阻率与黏度都将受熔体相主导，公式失效；此时或许可借助已有的熔体经验模型"换算"，但熔体比例目前难以精确获知，其结果只能作为参考。

结语与展望

这项工作从微观点缺陷化学与电中性原理出发，为位错蠕变控制的上地幔建立了一套有明确物理依据的电阻率—有效黏度换算框架，并成功把新疆西北部的 MT 电阻率成像"译"成了上地幔黏度结构，结果与实验室外推和独立强度模型在误差范围内一致。相比以往依赖经验假设或动力学建模标定的做法，本文方案物理图像更清晰，且能区分干/湿不同物化条件下的不同换算系数。

展望而言，电阻率作为定量的约束，有望更精细地揭示地幔黏度的横向与纵向变化，帮助理解大陆动力学演化。随着矿物物理实验的推进与地球化学、地球物理数据的积累，该模型也将被进一步修正与完善——例如对各电导/变形机制之间 T–P 相界更精确的刻画，以及对地幔含水量空间分布的进一步约束。

延伸阅读与引用

  • Karato, S., & Jung, H. (2003). Effects of pressure on high-temperature dislocation creep in olivine. Philosophical Magazine, 83(3), 401–414.
  • Gardés, E., Gaillard, F., & Tarits, P. (2014). Toward a unified hydrous olivine electrical conductivity law. Geochemistry, Geophysics, Geosystems, 15(12), 4984–5000.
  • Liu, L., & Hasterok, D. (2016). High-resolution lithosphere viscosity and dynamics revealed by magnetotelluric imaging. Science, 353(6307), 1515–1519.
  • Xu, Y., Zhang, A., Yang, B., et al. (2018). Bridging the connection between effective viscosity and electrical conductivity through water content in the upper mantle. Scientific Reports, 8, 1771.
  • Liu, S., et al. (2024). Footprints of continental amalgamation from three-dimensional resistivity structure: insights from the Tianshan and surrounding basins, China. Tectonophysics, 874, 230234.

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