过去首次约会之后：用于锆石铅丢失分析的集成谐和–不谐和比较方法

《Geochemistry, Geophysics, Geosystems》：Past the First Date: Ensemble Concordant–Discordant Comparison for Zircon Pb-Loss Analysis

　　

锆石铀–铅（U–Pb）年代学是地球科学中解决深时过程定年的基石，然而U–Pb体系并不总是表现为封闭计时器。在流体蚀变或热扰动过程中放射性成因Pb的选择性丢失可导致同位素年龄偏移并产生不谐和分析结果。部分锆石数据集可能保存了多期铅丢失事件，然而许多铅丢失模拟方法仅返回单一“最佳”铅丢失时间。谐和–不谐和比较（CDC）检验通过寻找重建的不谐和分析结果与谐和群体最佳匹配的模拟扰动年龄来估计铅丢失时间。在此，研究人员扩展了CDC框架，在蒙特卡洛实现中保留受支持的局部最优解，而非将每次实现缩减为单一最优解，并报告描述峰位运行间可重复性而非地质年龄不确定性的95%稳定性界限。这一重新设计能够解析多期铅丢失事件的时间，恢复叠加的变质、热液或岩浆事件。将该方法应用于具有规定铅丢失历史的模拟锆石U–Pb数据集，研究人员发现集成CDC在所有基准模拟方法中产生最低的总体中位数绝对误差（19 Myr），并为36个真实事件中的27个赋予了年龄估计。六个单事件基准示例的误差最低，90分位绝对误差为11 Myr，在重叠多事件情形中增加至83 Myr。应用于太古宙伊尔干花岗岩时，集成方法解析出454 Ma和1113 Ma两个铅丢失峰，分别与区域泛非冷却和中元古代瓦拉库尔纳大火成岩省一致。

**论文解读：集成谐和–不谐和比较方法用于锆石铅丢失多期事件解析**

**一、研究背景与科学问题**

锆石铀–铅（U–Pb）年代学是地球科学中解析深时过程定年的基石性工具，广泛应用于地壳生长时限约束、沉积物物源示踪以及变质叠加事件的厘定。该方法依赖于²³⁸U–²⁰⁶Pb和²³⁵U–²⁰⁷Pb两个衰变体系之间的年龄一致性：当两个体系在不确定度范围内给出等效年龄时，分析结果被描述为谐和的；反之则为不谐和的。然而，U–Pb体系并非始终表现为封闭计时器。辐射损伤和流体介导的放射性成因Pb迁移可通过将Pb从晶格中移除而扰动这一一致性。在最佳情况下，铅丢失是近期的且对²⁰⁷Pb/²⁰⁶Pb年龄无影响，但在许多情况下不谐和性使解释复杂化，最坏情况下导致数据被直接剔除。然而，不谐和组分可能保留有关流体蚀变历史的宝贵地质信息。

不谐和性可源自多种通常协同作用的过程。冲击诱发的微结构和晶体塑性变形拓宽了扩散通道，使锆石易受后期铅丢失脉冲的影响；α反冲造成的辐射损伤逐渐使富U区域非晶化，一旦这些蜕晶化域相互连通，低于约300 °C的循环流体即可轻易剥离放射性成因Pb。此外，后期U进入或U增益/Pb丢失耦合反应也可能扰动母子体平衡。由于这些机制在不同时间和尺度上作用，锆石群体常记录多次扰动痕迹，产生复杂的不谐和模式。铅丢失通常表现为谐和图中的一条或多条不谐和线，每条近似为上方截距结晶年龄与下方截距扰动年龄之间的弦。

现有的不谐和数据管理策略包括化学研磨热电离质谱法（CA-TIMS）等物理预处理方法和计算手段。然而，多数计算方法假设或强调单一铅丢失事件，使多事件情形难以解析。先前的谐和–不谐和比较（CDC）检验虽可恢复最古老的真实铅丢失事件，但在多阶段情形中存在局限性：当所有分析经历多期铅丢失时，CDC方法可能在虚假年龄处产生次级极小值，且无法区分真实地质信号与统计伪迹。

**二、研究内容与主要技术方法**

研究人员提出了一种集成CDC工作流程，保留局部最优解并以稳定性界限和支持度量总结铅丢失年龄。研究采用以下关键技术方法：（1）谐和性分类：通过百分比不谐和截止或谐和线相交检验将分析预先分为谐和与不谐和两类；（2）蒙特卡洛不确定性传播：依据报道的不确定度对每次分析进行抖动；（3）下截距网格搜索：在固定年龄网格上评估候选铅丢失年龄；（4）不谐和分析重建：将不谐和点通过候选下截距反投影至谐和线获得重建上截距年龄；（5）Kolmogorov–Smirnov（K–S）拟合优度检验与无效年龄惩罚：比较重建年龄分布与谐和年龄分布并施加无效重建比例惩罚；（6）逐峰检测与集成目录构建：采用scipy标准峰查找例程识别局部极大值，通过突出度、分离度和宽度阈值筛选，并经抛物线插值精化峰位，最终通过逐点中位数聚合和可重复性投票生成受支持峰目录。

**三、研究结果**

**3.1 合成基准测试性能**

在Cases 1–7和Tiers A–C中，CDC具有最低的总体中位数绝对误差（19 Myr），并为36个真实事件中的27个赋予年龄估计（DD为16个，DD–PEAKS为22个）。包含率为CDC 47%、DD 8%、DD–PEAKS 17%。

**3.2 单阶段基准（Cases 1–4）**

在所有三个层级中，CDC在中位数绝对误差和包含率方面均优于DD。Tier B中CDC中位数误差最小（3 Myr，67%包含率），但每个层级匹配估计数量有限。DD返回的匹配估计误差较大且界限不包含真实年龄。

**3.3 两阶段基准（Cases 5–7）**

随着层级从A到C散射增加以及两期事件相互靠近，恢复能力下降。CDC返回11个可匹配估计（共18个真实事件），较老1500 Ma事件在所有9个案例–层级面板中均被匹配，较年轻500 Ma事件仅在6A和6B中被匹配。

**3.4 事件分配**

CDC在每个层级中具有最高匹配率，DD–PEAKS通过报告更多峰提高恢复率但也增加未匹配估计。在匹配事件中，CDC在每个层级具有最小中位数绝对误差。

**3.5 区间宽度与参数敏感性**

CDC中位稳定性界限宽度为50 Myr（所有层级），DD为20–30 Myr。在756个非默认案例–层级–参数组合中，745个（98.5%）精确重现默认目录。在六个单事件面板中CDC绝对误差中位数为2 Myr，90分位为11 Myr；多事件面板中90分位误差升至83 Myr。

**3.6 扇形至零点基准（Case 8）**

当不谐和群体不定义共同下截距时，工作流不报告年龄。三个层级均未产生可报告的铅丢失年龄，表明该方法能识别不支持的几何形态。

**3.7 案例研究**

对来自Burtville地体的太古宙花岗岩样品193435，集成CDC解析出454 Ma（稳定性444–484 Ma）和1113 Ma（稳定性826–1255 Ma）两个峰，分别对应泛非冷却和瓦拉库尔纳LIP事件。单一最优值结果1108 Ma的界限（453–1249 Ma）过宽，跨越两个集成峰。

**四、讨论总结**

在单事件数据集中，集成方法与旧版摘要通常返回相同年龄，但额外提供支持度量和未解析结果。在多事件数据集中，可保留旧版缩减为单一值时丢失的受支持次级峰。CDC和DD的一致性可降低方法选择伪迹的可能性；两者不一致时，推断年龄可能对方法固有假设更敏感。CDC依赖谐和子集代表上截距年龄结构，当不谐和性具有年龄依赖性时可能产生偏差。DD不需要谐和参考但受最重置分析强烈影响，推断下截距可能偏向最重置部分。

**五、研究结论**

研究人员提出了集成CDC工作流，将单一最优K–S网格搜索扩展为多峰集成摘要。所报告的峰附有从蒙特卡洛重采样推导的稳定性界限和支持度量，无清晰内部峰的曲线被报告为未解析或近期边界模式，而非强制赋予离散铅丢失年龄。集成CDC总体中位数绝对误差最低（19 Myr），为36个真实事件中的27个分配估计。应用于太古宙Yilgarn花岗岩样品193435时解析出454 Ma和1113 Ma两个铅丢失峰，分别与区域早古生代冷却和中元古代瓦拉库尔纳LIP一致，表明集成摘要可保留单一最优结果所抑制的地质上有意义的结构。CDC方法的明确局限性在于其依赖于谐和子集对不谐和群体上截距年龄结构的代表程度，无效年龄惩罚和峰选择阈值在剖面平坦或事件强重叠时可能使方法趋于保守。在这些限制内，集成工作流提供了一种从单一不谐和锆石群体中识别、总结和报告多个受支持铅丢失年龄的实用方法。重要的是，不谐和数据不必简单过滤剔除，它们可保留关于流体–岩石相互作用的信息，特别是在辐射损伤使颗粒更易受叠加影响从而成为后期事件更敏感记录器的情况下。