刘善科等-JP&CG：花岗质岩石中锂同位素行为的界面动力学机制

锂（Li）是稀有金属成矿和花岗质岩浆演化过程中的关键不相容元素，锂同位素（δ 7 Li）被广泛用于示踪岩浆结晶、熔体—流体分离和水岩相互作用。然而，许多花岗岩—伟晶岩体系中常出现“Li含量越高、δ 7 Li值越低”的反相关现象（图1），这一特征难以用传统的平衡结晶或分离结晶模型解释（图2）。与此同时，花岗质矿物中Li与K同属碱金属元素，但二者是否受相同晶体化学机制控制仍缺乏清晰认识，是制约锂同位素示踪花岗质岩浆演化的关键问题。

图1 全球花岗岩—伟晶岩体系中锂含量与锂同位素组成的关系。全球花岗岩和伟晶岩样品显示出显著的Li含量和δ 7 Li变化范围（a, b），部分代表性体系表现出“Li含量升高而δ 7 Li降低”的反相关趋势（c, d）

图2 平衡结晶和分离结晶模型对花岗岩—伟晶岩Li–δ 7 Li演化趋势的约束。基于岩石学模拟和矿物—熔体Li分配/同位素分馏参数建立的平衡结晶（EQ）与分离结晶（FC）模型表明，传统结晶过程能够导致残余熔体中Li富集，但难以单独再现天然样品中广泛观察到的δ 7 Li–Li反相关关系。这说明经典岩浆分异模型不足以解释花岗岩—伟晶岩体系中的锂同位素行为

围绕这一科学问题，中国科学院地质与地球物理研究所刘善科等从全岩—矿物两个尺度开展系统研究。一方面，整合全球花岗岩和伟晶岩Li含量与δ 7 Li数据，结合AlphaMELTS软件模拟平衡/分离结晶过程，检验传统结晶模型是否能够产生观测到的δ 7 Li–Li反相关。另一方面，以中国七个花岗质地区的钾长石为对象，测定其Li同位素组成，并与已发表的K同位素、Al/Si摩尔比数据相结合，比较Li、K同位素在同一矿物体系中的响应机制。研究表明：

1. 平衡结晶和分离结晶均难以单独再现花岗岩—伟晶岩体系中广泛出现的δ 7 Li–Li反相关（图2）。相比之下，当晶体快速生长而熔体扩散补给不足时，晶体—熔体界面可形成富Li边界层；由于 6 Li 扩散快于 7 Li，界面附近可产生瞬时Li同位素分馏，从而使低Li矿物或早期结晶相记录相对较重的δ 7 Li，而高Li残余熔体或晚期体系表现为较低δ 7 Li（图3）。

2. 钾长石的分析结果表明，δ 7 Li与Al/Si摩尔比解耦，而K同位素则与Al/Si控制的晶体化学环境密切相关，说明Li同位素主要反映不相容元素的传输受限和界面边界层过程，K同位素则更多记录矿物结构和晶体化学分馏（图4）。

该系列研究提出了花岗质岩石中“Li同位素边界层动力学控制—K同位素晶体化学控制”的新认识，为解释花岗岩—伟晶岩体系中Li同位素异常分馏提供了统一框架，也为利用Li–K同位素联合识别岩浆分异、快速结晶、熔体—流体演化过程提供了新的理论依据。

图3 晶体—熔体界面扩散驱动的锂同位素分馏机制。在快速晶体生长和扩散补给受限条件下，不相容元素Li可在晶体—熔体界面形成富集边界层。由于轻同位素 6 Li 相对于 7 Li扩散更快，界面附近可产生瞬时同位素分馏，使晶体或近界面熔体记录相对较重的δ 7 Li，而远离界面的残余熔体则可能演化出较低δ 7 Li。该机制为解释花岗岩—伟晶岩体系中Li富集与δ 7 Li降低的反相关现象提供了动力学框架

图4 花岗质钾长石中Li–K同位素解耦及其不同分馏机制。钾同位素主要受Al/Si比值调控的晶体化学环境控制，反映K在长石结构位中的近似平衡分馏；相比之下，Li作为高度不相容微量元素，其同位素组成与Al/Si比值解耦，更容易受到晶体—熔体界面边界层、扩散动力学和局部生长条件的影响。Li–K同位素解耦揭示了花岗质岩石中主量结构元素与微量不相容元素在同位素分馏机制上的根本差异

研究成果发表于JP&CG：

（1）刘善科，苏本勋，李文君. Kinetic isotope fractionation of lithium in granite–pegmatite systems: The role of diffusion at crystal–melt interfaces[J]. Journal of Petrology, 2026. DOI: 10.1093/petrology/egag039 .

（2）刘善科，苏本勋，王世莹，李文君. Al/Si ratio controls on δ 41 K in granitic K-feldspar: Decoupling Li-K isotopes and contrasting fractionation controls[J]. Chemical Geology, 2026, 713: 123417. DOI: 10.1016/j.chemgeo.2026.123417 .