磷灰石U-Pb年代学与热年代学再审视：从纳米尺度到造山尺度

《Earth-Science Reviews》：Thermochronology and the apatite U-Pb geochronometer revisited: From the nanoscale to the mountain-scale

　　

磷灰石U-Pb年代学已成为约束地壳环境中地质过程时间和持续时间的强大工具。然而，若干案例研究表明，当前的解释框架仍存在问题。这些困难部分源于对实验数据及控制放射性成因Pb丢失机制的理解不完整，部分源于对封闭温度等理论概念或存在固定温度界限的部分保留带（PRZ）的误解。不同类别的热历史如何影响放射性成因Pb丢失仍不确定，尤其是在围岩变形和变质促进热流体运动±矿物本身晶体塑性行为的条件下。因此，研究人员审查了：i）与部分保留部分丢失（PRPL）带相关的理论概念；ii）温度-时间（T-t）对Pb丢失的影响，使用目前采用的扩散参数；iii）现有实验数据以及研究人员为何采用表面屏障加晶格（SBPL）模型重新分析以获得磷灰石中Pb的校正扩散参数；iv）使用校正扩散参数的不同类别T-t路径的影响；v）扩散理论应用的实际方面，包括说明解释陷阱的案例研究。研究人员还为热年代学应用提供了实用指南。在若干方面，可以与⁴⁰Ar/³⁹Ar年代学进行有益类比，后者在基于激光 spot 分析的研究与基于温度控制炉逐步加热实验释放气体的时间年龄剖析的研究之间存在类似的张力。研究还表明，隐蔽交代作用可能同时影响共存磷灰石和钾长石中的同位素分布，进一步加强了磷灰石U-Pb与⁴⁰Ar/³⁹Ar年代学之间的联系。

研究背景方面，磷灰石U-Pb年代学已被广泛用于约束地壳环境中地质过程的时间与持续时间，但现有解释框架存在显著问题。部分案例研究（如Flansburg et al., 2026）将磷灰石U-Pb年代计视为中温范围（约450°C）温度计，假定部分保留带（PRZ）温度范围为350–570°C，并据此否认爱琴海变质核杂岩的存在。研究人员通过反证法指出，若变质岩石学关于温压条件的约束成立，则上述PRZ假设必然有误。问题根源在于对实验数据的误解以及对封闭温度、PRZ等理论概念的误用，特别是目前广泛使用的Pb扩散参数（Cherniak et al., 1991摘要值）可能过于不保留。为此，研究人员开展了系统性再分析，采用表面屏障加晶格（SBPL）模型重新解释实验数据，获得更保留的扩散参数，并据此重新评估PRPL带温度界限、不同T-t路径下Pb丢失行为以及磷灰石U-Pb体系的流体敏感性。研究得出重要结论：校正后的扩散参数表明磷灰石对Pb高度保留，PRPL带下限不可能低至350°C；固定温度界限的PRZ不能用作地质温度计；磷灰石U-Pb体系应归类为流体敏感年代计，年龄簇更可能反映流体脉冲和幕式构造热事件而非缓慢冷却。该研究对热年代学解释框架具有根本性修正意义，论文发表于《Earth-Science Reviews》。技术方法方面，研究人员主要采用：对Cherniak et al.（1991）和Watson et al.（1985）磷灰石Pb扩散实验数据进行时间序列再分析，识别扩散率随 elapsed time 的系统性下降以区分单域扩散（SDD）模型与SBPL模型；使用LINEST函数进行渐近线回归确定晶格扩散参数；基于Dodson（1973, 1986）公式和Crank（1975）扩散方程进行数值模拟；引入准稳态（QSS）理论推导理论最大年龄（t_max）和达到准稳态时间（t_qss）；利用Watson and Cherniak（2013）公式评估热脉冲韧性；采用等效半径近似处理柱状几何；通过Python脚本生成热图。样本数据来源于已发表的实验数据集，无新样本队列。研究结果方面，第一，单域与多域扩散模型再分析表明，Cherniak et al.（1991）数据中扩散率随实验持续时间系统性下降，不符合SDD模型预期，而符合序贯SBPL模型；通过剔除低温序列数据获得明确定义的渐近线（R²=0.996，t-统计量>22），据此提出校正扩散参数：活化能E=311±14 kJ mol^?1，频率因子D₀=1.01×10^?4 m² s^?1。第二，与Watson et al.（1985）数据再分析对比表明，三点回归给出E≈323 kJ mol^?1、D₀≈4.59×10^?5 m² s^?1，两套数据边际1σ区间重叠，均显示磷灰石远比先前认为的更保留。第三，活化体积估算表明，V?_Pb=[1–2]Ω_m对应[0.74–1.50]×10^?5 m³ mol^?1，压力使Dodson温度在1.0 GPa下升高约19–39°C、2.0 GPa下升高约38–78°C，相当于0.5–1.0°C km^?1。第四，PRPL带界限重新定义表明，基于Dodson（1986）公式的快速冷却条件下PRPL带范围为[510–645°C]（粒径10–100 μm），基于理论最大年龄的静态条件下100 μm直径球体PRPL带约为481–569°C，均远高于350°C下限。第五，韧性分析表明，后续热脉冲可穿越PRPL带而不完全重置已有年龄群体；等压逐步加热持续10 Ma时，100 μm直径磷灰石需峰值温度超过500°C才能使体积分数损失超过三分之一。第六，变形加速扩散与流体敏感年代计讨论表明，位错阵列、宽核倾斜边界和动态维持渗透性可促进流体介导的隐蔽交代作用，导致磷灰石U-Pb体系重置，这解释了年龄簇与构造事件吻合的现象。讨论部分总结如下：研究人员强调，固定温度界限的PRPL带不应作为地质温度计使用，PRPL带是速率敏感的动态实体，其界限随粒径、冷却速率和P-T-t路径而变化；建议用PRPL带概念取代PRZ概念；实验确定扩散参数时必须进行同温度不同持续时间序列以区分SDD与SBPL模型；磷灰石U-Pb与钾长石⁴⁰Ar/³⁹Ar高保留核域年代学类似，均属流体敏感年代计，隐蔽交代作用和流体辅助晶界迁移可能是年龄重置的主要机制；年龄簇与构造事件吻合需要更高峰值温度或更长热脉冲，并需常规考虑隐蔽交代作用。研究结论翻译如下：基于应用表面屏障加晶格（SBPL）模型对Cherniak et al.（1991）实验数据进行数值再分析，研究人员提出磷灰石中Pb扩散的新参数：活化能E=311±14 kJ mol^?1，频率因子D₀=1.01×10^?4 m² s^?1，log₁₀D₀的1σ不确定度对应D₀范围为[0.213–4.75]×10^?4 m² s^?1。基于与其他晶质材料的比较，以平均原子位体积Ω_m为单位，估算磷灰石中Pb扩散的活化体积V?在[1–2]Ω_m范围内，对应[0.74–1.50]×10^?5 m³ mol^?1。对特定压力-温度-时间（P-T-t）路径效应的定量估计需要基于QTQt等计算机程序对特定微观结构参数进行精确数值模拟。从半定量角度看，校正扩散参数意味着磷灰石Pb扩散的任何部分保留部分丢失（PRPL）带下限不太可能低于约480–520°C（取决于粒径）。快速冷却条件下PRPL带上限的半定量估计也更高，为589–645°C，而对于稳态环境温度，适用544–596°C的较窄范围。压力效应也需考虑。假设平均晶质岩石密度为2700 kg m^?3，岩石静压力梯度约为26.5 MPa km^?1或0.265 kbar km^?1。因此，10和20 kbar压力分别对应约38和75 km深度。预测的压力效应相当于Dodson温度升高约0.5–1.0°C km^?1。磷灰石PRPL带的中低温界限不能假定为较老磷灰石颗粒生长后T-t路径上所经历最高温度的上限。研究人员记录了允许大幅越界的情形。无论如何，先前提出的磷灰石PRPL带低温下限过低。除边缘带微米级Pb丢失外，即使经过显著时间，在低至350°C的温度下也不会记录到可观测效应（当然不包括流体侵入引起的任何效应）。磷灰石U-Pb年龄梯度可用于热年代学，但模拟表明热事件后快速冷却达到的温度比以前认为的更高。通常基于不同封闭温度矿物年龄分布推断的“缓慢冷却”与年龄剖析数据不一致。相反，可能涉及一系列热脉冲，每个热脉冲与变形和/或变质幕相关，并巧合地促进热流体运动，包括岩浆作用释放的流体。U-Pb年龄群体与已知构造事件的吻合需要更高峰值温度或更长热脉冲，并需常规考虑流体运动引起的隐蔽交代作用。