田雨露|王佳敏等：硬柱石榴辉岩退变质过程中微量元素迁移规律及其对全岩化学成分的影响

田雨露 , 王佳敏 *, Gültekin Topuz, Aral I. Okay, 刘通 . Trace element mobility and mineral-scale redistribution during retrograde metamorphism of lawsonite-bearing high-pressure rocks from Elekdağ, Turkey[J]. Contributions to Mineralogy and Petrology, 2026, 181: 73. DOI: 10.1007/s00410-026-02358-1 .

俯冲带是连接地球表层与深部物质循环的重要通道。洋壳在俯冲和折返过程中，会经历高压变质、脱水反应及复杂的流体—岩石相互作用，由此引起水和微量元素在矿物、流体及岩石之间不断重新分配。这些过程不仅影响俯冲板片向深部输送物质的效率，也关系到地幔楔交代、岛弧岩浆形成以及全球地球化学循环。因此，一个长期受到关注的科学问题是：俯冲带流体影响下的元素迁移在不同尺度上是如何表现的？

传统上，K、Rb、Ba等大离子亲石元素（LILE）通常被视为流体活动性元素，而Nb、Ta、Zr、Hf等高场强元素（HFSE）则被认为相对稳定。然而，天然高压变质岩记录的元素行为往往比这种简单的“活动—不活动”二分更加复杂。特别是在岩石折返过程中，外来流体可以在岩石尺度上造成明显交代，而峰期矿物分解释放出的元素又可能在矿物反应界面附近被新生矿物迅速重新固定。因此，全岩地球化学组成往往同时包含原岩继承、流体交代和矿物反应等多重信息。

针对这一问题，中国科学院地质与地球物理研究所田雨露博士研究生，在王佳敏副研究员指导下，与刘通副研究员及土耳其伊斯坦布尔技术大学Gültekin Topuz教授和Aral I. Okay教授合作，对土耳其北部Elekdağ变质蛇绿岩中的硬柱石榴辉岩和不同程度退变质的榴辉岩开展了系统研究。综合岩相学、矿物主微量元素分析、相平衡模拟、SIMS金红石U–Pb年代学及微量元素质量平衡计算，重建了这些高压变质岩石从俯冲到折返过程中的P–T–t演化，并进一步追踪不同元素在不同尺度上的迁移方式。

该研究首先揭示了Elekdağ榴辉岩典型的冷俯冲—折返退变质演化过程，并建立了相应的发夹状P–T轨迹，以及变质演化过程中的矿物含量变化。

岩相学观察显示，不同样品保存了程度不同的高压峰期矿物组合及退变质改造特征（图1），保存较好的榴辉岩主要发育石榴子石、绿辉石和硬柱石等高压矿物，其中石榴子石内部广泛发育硬柱石、绿辉石、金红石及多硅白云母等高压矿物包裹体，部分硬柱石仍保留完整的矩形或菱形晶形。随着岩石折返，峰期矿物逐渐发生分解和交代：硬柱石被绿帘石/黝帘石、多硅白云母等矿物组成的假象所替代，金红石边部也可发育榍石反应边。更强烈退变质的样品则形成以蓝闪石、阳起石、绿泥石和绿帘石族矿物为主的退变质矿物组合。这些显微结构完整记录了岩石从硬柱石榴辉岩相峰期、绿帘榴辉岩相、再向蓝片岩相退变的演化过程。相平衡模拟表明，保存较好的榴辉岩峰期变质条件达到约2.3–2.6 GPa和530–590 ℃，对应约75–80 km深度的冷俯冲环境（图2）；蓝片岩相退变质的温压条件为~1.2 GPa和~380 ℃，对应约40 km深度的地壳深度；记录了发夹状的P-T轨迹。SIMS金红石U–Pb定年获得102 ± 14 Ma的年龄，与Elekdağ地区已发表的早白垩世峰期高压变质作用年代学结果一致。

图1 Elekdağ硬柱石榴辉岩的典型显微结构和退变质反应特征。a.斑状变晶结构；b.含矿物包裹体的石榴子石；c.绿辉石及退变质共生体；d.ED68石榴子石包裹体；e.矩形硬柱石包裹体；f.硬柱石假象（Czo+Ph）；g.WT13石榴子石包裹体中的硬柱石及硬柱石假象（Czo+Ph）；h.蓝片岩

图2 Elekdağ硬柱石榴辉岩代表性样品WT13的P–T–t演化轨迹

其次，在建立变质演化框架的基础上，该研究进一步关注流体活动是否改变了岩石原有的地球化学组成。

全岩微量元素结果显示（图3），不同样品之间存在明显差异：部分样品仍保存较典型的MORB型特征，而另一些样品则明显富集K、Rb、Ba、Cs和B等流体活动性元素。矿物原位分析和质量平衡结果进一步表明，多硅白云母是Rb、Ba和Cs等LILE的重要载体，并具有明显的B富集；与此同时，退变质过程中由硬柱石分解形成的绿帘石/黝帘石表现出明显的成分变化，并成为Sr和LREE等元素的重要新生载体，记录了这些元素在矿物反应过程中的重新分配。全岩K2O与Rb、Ba之间的高度相关性也表明，这些元素很可能随同一类富K流体共同进入岩石。上述全岩和矿物尺度证据共同说明，Elekdağ榴辉岩在折返过程中发生了外来流体与岩石之间的物质交换。

图3 榴辉岩全岩及主要矿物的微量元素地球化学特征。a–c.全岩稀土元素和微量元素配分模式；d–f.绿帘石/黝帘石、硬柱石及多硅白云母的微量元素特征

这一认识也带来一个重要问题：现今测得的全岩化学组成是否仍能代表岩石峰期变质时的化学组成？以其中一个高K样品为例，其部分K可能是在折返阶段由外来流体带入。如果直接采用现今全岩成分进行峰期相平衡计算，就可能把后期加入的K“带回”峰期，从而影响对峰期矿物比例和元素预算的判断。

因此研究团队开展了K2O敏感性相平衡模拟，系统改变输入全岩的K2O含量，并计算不同K2O含量条件下的峰期矿物组合和微量元素预算（图4）。结果显示，K2O含量的变化主要控制多硅白云母的稳定性和含量，而对主要高压矿物组合及峰期P–T框架影响相对有限。随着K2O含量升高，计算得到的多硅白云母比例明显增加，Rb、Ba、Cs等LILE的整体预算也随之升高。综合全岩微量元素、多硅白云母原位成分以及K2O敏感性模拟结果，研究认为，Elekdağ榴辉岩在折返过程中经历了外来富K、B水流体的渗入。这些流体向岩石中输入K、Rb、Ba、Cs和B等流体活动性元素，并促进多硅白云母形成，使岩石在宏观尺度上表现出明显的开放体系行为。

图4 不同变质阶段的微量元素组成及K2O含量变化的影响。a–b.不同变质阶段的REE和MORB标准化微量元素组成；c–d.不同K2O含量（0、0.4、0.7和1 wt.%）条件下模拟的峰期REE和微量元素组成

然而，外来流体的进入并不意味着所有微量元素都会发生大尺度迁移。

矿物原位分析和质量平衡结果显示（图5），LREE、Sr以及HFSE表现出了明显不同于K、Rb、Ba和B的行为。硬柱石是俯冲洋壳中重要的含水矿物，也是LREE和Sr的重要储库。在折返过程中，硬柱石逐渐分解并释放这些元素，但LREE和Sr并未大量离开岩石，而是主要被邻近新生的绿帘石/黝帘石重新捕获。特别是WT13样品中的黝帘石表现出系统的成分分带和Sr–Fe变化关系，记录了元素从硬柱石反应域向新生矿物短距离转移的过程。类似地，峰期金红石被榍石交代时，Nb、Ta等HFSE也主要在新旧Ti矿物之间发生重新分配，而未表现出显著的大尺度迁移。

图5 质量平衡计算结果显示不同退变质程度的样品中矿物对微量元素的贡献比例。a-c.基于全岩及矿物主量元素组成计算； d-f.基于相平衡模拟获得的矿物含量计算

这些现象表明，矿物反应可以在微观尺度上形成高度局部化的元素再分配。峰期矿物分解释放出的元素能够借助流体参与反应，但其迁移距离受到反应界面和新生矿物的强烈限制，并通过界面耦合溶解—再沉淀等机制迅速被重新固定。这些结果同时说明，矿物反应不仅控制矿物组合的演化，也影响元素究竟会随流体迁移，还是继续保留在岩石内部。

最后，综合全文研究结果，该研究提出了俯冲带硬柱石榴辉岩折返退变质过程中微量元素迁移的尺度依赖性机制，详细阐述了其对全岩化学成分的影响。在岩石尺度上，外来流体能够输入K、Rb、Ba、B等高度流体活动性元素，使体系表现为开放；而在矿物反应尺度上，LREE、Sr和HFSE等元素则主要在局部反应域内重新分配并被新生矿物捕获，整体仍受到较强的岩石内部缓冲。这种“宏观流体交代—微观矿物缓冲”的耦合过程，是形成复杂全岩地球化学特征的重要原因。

该研究表明，退变质高压岩石最终保存的全岩地球化学组成，是原岩继承、外来流体交代和矿物尺度元素再分配共同作用的结果。因此，仅依据全岩微量元素特征恢复原岩性质或判断构造环境，可能将后期流体和变质反应产生的信号误认为原始岩浆地球化学特征。将全岩地球化学与矿物尺度微量元素、变质P–T–t演化及质量平衡结合，是更准确理解俯冲带流体—岩石相互作用和元素循环的重要途径。