俯冲带是地球深部水循环的关键通道。海洋板片在俯冲前于海沟外侧发生挠曲，形成一系列外隆起断层。海水可沿这些断层进入洋壳和岩石圈地幔，引发橄榄岩蛇纹石化反应，将大量结构水固定在板片中，并随俯冲进入地球深部。该过程不仅影响全球水循环，也与板片强度、俯冲动力学、弧岩浆作用和地震活动密切相关。

长期以来，大洋岩石圈地幔在进入海沟时究竟携带多少水进入俯冲带仍存在较大争议。传统估算往往假定地幔在一定深度范围内发生较均一的蛇纹石化，但近年来地震各向异性、断层渗透率和数值模拟研究表明，蛇纹石化更可能集中在外隆起断层及其邻近破裂带中（图1a）。两类认识导致全球板片地幔水通量估计相差显著。造成这一分歧的关键在于：外隆起断层如何控制流体进入深度、蛇纹石化空间分布和最终水通量，目前仍缺乏动力学约束。

图1 (a) 俯冲带外隆起板块地幔蛇纹石化示意图；(b) Yang et al. (2023)汇总的外隆起地震深度分布 (Mw>5)

针对这一问题，中国科学院地质与地球物理研究所岩石圈演化与环境演变全国重点实验室博士研究生张睿与其导师杨建锋、赵亮，以及合作者张畅、刘通、邓林培和王新，开展了高分辨率地球动力学数值模拟研究。研究团队将板片挠曲断层发育、温度控制的蛇纹石化反应动力学以及蛇纹石化导致的流变弱化过程耦合起来，系统评估外隆起地幔蛇纹石化及其水通量的控制机制。

模型采用黏弹塑性流变，显式刻画板片进入海沟前的弯曲破裂和塑性应变集中过程，并引入蛇纹石化反应动力学，用以描述温度对反应进程的控制（~270°C反应速率达到峰值）。在模型中，流体沿断层控制的通道进入Moho以下地幔，其最大有效进入深度用参数表示。研究团队进一步计算不同、板片年龄和反应速率条件下的蛇纹石化范围和水通量。

模型结果表明，外隆起地幔蛇纹石化并非简单受板片年龄控制，而主要受三个因素调节：流体最大进入深度、断层网络几何形态和蛇纹石化反应动力学。随着增大，蛇纹石化可向更深地幔发展，水通量先升高；当达到约17 km时，蛇纹石化效率达到峰值。继续增加流体进入深度后，蛇纹石化导致的弱化会使应力和应变过度集中到少数主断层中，抑制次级断层发育，使断层间距增大、横向水化范围减小，最终导致总体蛇纹石化体积和水通量下降（图2）。

图2 流体最大穿透深度（）对地幔蛇纹石化水通量的影响

这一结果揭示了蛇纹石化水通量的非线性控制机制：更深的流体进入并不必然带来更大的总水通量。只有当蛇纹石化深度和断层侧向分布之间达到适当平衡时，板片地幔才能实现最高效的水化。该认识解释了为什么局域化断层水化也能产生接近传统均一水化模型的总体水通量。

研究还发现，蛇纹石化与断层发育之间存在正反馈。板片挠曲形成初始断层和破裂带，为流体进入地幔提供通道；流体进入后触发蛇纹石化反应，导致地幔有效黏度降低，进一步促进应变集中和断层发育（图2e-图2h）。该反馈机制有助于在海沟附近形成高程度蛇纹石化的断层带，并控制水化区域的空间分布。

通过将模型结果与全球俯冲带总长度相结合，以=12km作为保守估算的中心值（误差区间取=5–17km），计算得到全球俯冲带水输入通量约为Tg/Myr (相当于将厚度约1 km的地幔完全蛇纹石化，或使厚度约12 km的地幔达到8%的蛇纹石化程度)。该数值表明，板片地幔携带的深部水并不需要依赖大范围均一蛇纹石化即可达到与既有全球水循环估算相当的量级；由密集分布、高度蛇纹石化的断层网络控制的区域水化，可以有效解释不同研究之间的水通量差异（图3）。

图3 (a)地幔蛇纹石化最大深度与地幔蛇纹石化程度的关系；(b)板块年龄与地幔蛇纹石化程度的关系；(c)文献和本研究中蛇纹石化地幔的水通量

该研究提出，外隆起地幔蛇纹石化应理解为一个由断层控制、反应动力学调制、并受水化—变形反馈影响的非线性过程。其意义在于，一方面为调和“均一水化”与“断层局域化水化”两类水通量估算提供了动力学解释；另一方面强调了断层间距、流体进入深度和反应速率在深部水循环中的关键作用。这种“断层局部化”的水化模式意味着，俯冲板片在进入地幔后释放流体也可能是脉冲式、非均匀的，这对于解释弧前岩浆活动的时空分布以及深源地震的成因具有重要的启示意义。未来，采用更真实的二相流模拟，耦合蛇纹石化微观动力学（蛇纹石化相变及体积变化）过程，将能更准确地描述俯冲板片蛇纹石化及水通量。

张睿（博士生）