基于CL（阴极发光）、U/Pb年龄及地球化学数据综合解读锆石物源的前向图像配准方法

《Geochemistry, Geophysics, Geosystems》：Forward Image Registration for Higher Level Interpretation of Zircon Provenance Based on Combined CL, U/Pb Age, and Geochemical Data

　　

碎屑锆石研究通常产生灰度阴极发光（CL）图像、U/Pb年龄及程度不一的微量元素浓度数据。然而，这些数据所包含的综合信息很少被用于更高层次的锆石物源解读，目前仍主要依赖样品尺度的U/Pb年龄谱对比。本研究提出一种新型的前向图像配准软件解决方案，可实现空间配准数据的整合（光学显微成像与CL纹理/强度/对比度/颜色；Pb/U比/年龄；以及微量元素浓度/比值），从而支持基于数据科学的锆石物源分析。该方案基于开源脚本构建，可对接仪器软件、开源生物信息学图像分析软件以及研究人员自主开发的图像处理算法。配备此工具集后，分析人员可在主控数据表中可视化已配准数据集中的全部锆石数据。锆石颗粒的网格展示可按上述变量排序，便于直观发现变量间的相互依赖关系。研究人员以现代大型河流流域为例，说明可通过交互操作所有可用判据来手动定义、对比和发现新的分组与模式。基于年龄、化学组成和CL属性，研究人员为文献中记载的、具有花岗岩源岩的低频锆石群体提供了候选源区假说。此类探索有望为采用聚类-密度图实现更独特且自动化的沉积物物源分析开辟道路。

锆石（ZrSiO₄）是研究沉积物源最广泛使用的重矿物。碎屑锆石研究中对具统计意义的年龄群体的需求日益增长，这推动了高通量电子显微镜（SEM)及激光剥蚀-电感耦合等离子体质谱（LA-ICP-MS)或二次离子质谱（SIMS)等质谱分析技术的发展。然而，传统方法主要依赖U-Pb年龄谱进行物源对比，未能充分利用CL图像纹理、微量元素地球化学等多模态数据的综合信息。

当前研究面临的核心问题包括：（a）数据量大且缺乏结构化；（b）微分析点位选取和数据分析耗时费力；（c）地质记录越古老，保存的证据越有限。此外，现有公共数据库普遍缺乏原始CL图像，制约了数据的充分解读与再利用。为解决上述问题，研究人员开展了本项研究，旨在构建一种前向图像配准软件解决方案，实现显微成像与微分析实验的空间数据配准，并支持基于聚类方法识别具有共同特征和/或来源的锆石群体。该研究发表于《Geochemistry, Geophysics, Geosystems》。

研究人员采用的技术方法主要包括六个关键步骤：以澳大利亚墨累河入海口处的现代河流沉积物样品CA-MR24-1为研究对象，该样品经重液和磁选富集锆石后，置于25 mm圆形树脂靶中抛光制备。第一步，使用Evident VS200岩相切片扫描仪进行全片成像（包括透射平面偏光、正交偏光和反射光），并通过热场发射扫描电镜（FEG-SEM)获取背散射电子（BSE)和彩色CL图像（RGB探测器）。第二步，利用QuPath和ImageJ-Fiji软件配合BioVoxxel插件进行锆石颗粒的半自动分割，包含阈值分割和不规则分水岭分割两个掩膜步骤。第三步，在Fiji中使用BigWarp插件进行图像的像素级配准，实现BSE、CL及光学显微镜图像的精确对齐。第四步，将分割后的颗粒转入MATLAB脚本生成虚拟展示网格，可按CL强度、长宽比等标准排序。第五步，在QuPath界面上进行离线计算机点位规划，结合多模态影像信息优化激光剥蚀点位。第六步，通过平板扫描仪获取样品靶照片，经两步图像配准（CL图像到照片扫描、照片扫描到LA样品台），将虚拟点位坐标转换至激光剥蚀系统，实现自动化分析序列生成。

研究结果部分：

3.1 微分析位点的前向配准：颗粒分割在BSE镶嵌图上获得最佳效果，共识别534颗锆石，其中304颗具有足够大的面积容纳25 μm圆形激光剥蚀斑点。离线点位规划的优势在于允许样品拥有者与LA-ICP-MS专家在分析前协商确定最合适的斑点直径。斑点位置由U-Pb锆石年代学专家根据每颗锆石的主导CL特征进行编辑优化，避开明显的边部增生区和可能的继承核。配准结果显示虚拟点位与实际激光系统反射光图像之间具有高质量的对齐，大幅缩短了微分析实验的设置时间。

3.2 颗粒网格文档编制：获得分析数据后，可视化分割对象的颗粒网格功能得以充分发挥。研究人员初步过滤剔除了4个错误分析其他矿物的点位（Ti>3,000 ppm，Zr<30%）。网格排序功能的展示按9个独立标准进行：U-Pb年龄递增、全颗粒CL强度递减、CL颜色色调（蓝至黄）、内部CL纹理（ResNet50"循环"提取特征）、颗粒形态（长宽比和凸度）以及分析点位地球化学（Th/U、Gd/Yb、Eu异常）。每个网格贴片显示样品ID、分析点位位置（黄色虚线轮廓）和排序准则的值标签。以U-Pb年龄排序网格为例，可快速评估稀有"年龄"是否实为来自不同CL纹理区的混合分析。

3.3 年龄谱：对300颗锆石的Wetherill谐和图显示从约2,700 Ma至98.6 Ma的正常不谐和或谐和年龄，以及9颗强烈不谐和的锆石（指示早于2700 Ma的Pb丢失）。采用10%相对年龄差筛选"谐和"颗粒（n=66），整体不谐和度众数为+14%。核密度估计（KDE)年龄谱显示大量峰期，反映了墨累-达令盆地长期复杂的岩浆和热构造历史。最丰富的KDE峰期位于250 Ma和410 Ma，480-560 Ma之间存在宽泛的峰群簇。

4 讨论：

4.1 选定碎屑锆石群体的检视：研究聚焦于四个主观选择的时间带：（a）最年轻颗粒²⁰⁶Pb/²³⁸U年龄<142.7 Ma（n=13）；（b）二叠纪-三叠纪群体201.4-301.4 Ma（n=31）；（c）格林威尔期群体980-1,300 Ma（n=18）；（d）最老颗粒²⁰⁷Pb/²⁰⁶Pb年龄>1,828 Ma（n=17）。群体#1（最年轻）虽均为中白垩纪年龄，但内部CL结构和颜色差异显著，包含三颗蓝色CL响应锆石（99、110、115 Ma），可能源自惠特桑德火山省。群体#2（二叠纪-三叠纪）以振荡环带锆石为主，常见扇形叠加，显示主导的火成岩成因；基于相似锆石外观和²³⁸U/²⁰⁶Pb年龄（约3 Ma不确定性内），识别出约238 Ma（红色框）和247 Ma（蓝色框）两个三叠纪年龄亚群。群体#3和#4的CL图像显示明亮边部，复杂内部结构无法被单点分析完全捕获，暗示多旋回沉积再循环、变质作用、混合继承和长期沉积居留。形态凸度分析显示群体#1略低于其他群体，与其次棱角状或破碎颗粒组成一致。Hf和Zr/Hf分析表明群体#1可能指示较少演化的源岩浆。地球化学计算为锆石源区提供了初步近似，如群体#3和#4较低的Th/U暗示变质组分；Dy/Yb略高，暗示源区可能存在石榴石。群体#3具有限定良好的负Eu异常和有限的正Ce异常，与低氧逸度岩浆环境相容。仅群体#1的Ce/√(U⁰×Ti)高于0.3，与高于FMQ缓冲剂的岩浆氧化状态一致。

4.2 潜在碎屑物源点的搜寻：鉴于墨累-达令盆地广泛的地质知识，研究人员对明确界定的年龄、年龄-结构及年龄-结构-化学群体寻找沿排水系统的合理候选源区。群体#2的年龄区间在整个盆地中代表性不足，通过岩石类型年龄检索，识别出新英格兰岩基北部的Ruby Creek花岗岩作为群体#2中锆石的合理候选源岩之一。文献样品的SHRIMP U-Pb分析显示年龄组#1（~238.3±3.3 Ma）与样品2164940（241.7±1.7 Ma）一致，年龄组#2（246.3±3.1 Ma）与200701550和200701548（约246 Ma重叠年龄）一致。U浓度和Th/U比值的KDE分布显示对应年龄组间存在广泛重叠（Th/U峰值0.4-0.6）。详细观察可将群体#2年龄组分为两个纹理组：（a）暗核带振荡边部；（b）振荡环带叠加扇形环带。综合考虑年龄、成分和纹理匹配，仅年龄组#2（246.6 Ma）通过了与潜在源岩的联合检验。

4.3 前向配准的约束、挑战与改进机遇：（a）图像分析方面，网格显示目前为"静态"且排序偏好单一标准；动态沉浸式可视化及计算机视觉量化纹理的自动化浏览匹配是后续发展方向。（b）图像处理方面，图像贴片大小对目标对象大小的适应性不足，无法将大颗粒长晶体尖端完全纳入；且无法实现晶体学定向。（c）数据采集方面，Apreo 2帧大小可进一步扩大，需测试不同样品制备方法；激光系统平台精度对大规模自动化至关重要。

5 结论：

前向配准方法是一种用于相关性显微分析、碎屑物源假说生成和群体识别的新型图像分析工具。通过组装高维数据集、构建矿物展示网格以及从选定锆石群体中识别候选源岩的经验，研究人员证实：（1）序列显微分析和微分析实验数据的早期整合产生协同效应，有利于新样品的地质解读，主要要求是涉及灵活（MATLAB、Python、QuPath脚本）和用户友好（QuPath）软件的一致数据采集和图像分析；（2）网格和分析位点空间参考的定制与保存（QuPath多视图界面、MATLAB主控表过滤/排序）简化了预处理和后处理分析过滤、比较及潜在颗粒群体的解释，比较更为便捷因其以有序颗粒和更小文件体量提示用户软件，网格具有互操作性且可基于多个样品和历史实验构建；（3）人类解释者必须参与分析循环以理解和交互复杂纹理的多模态图像（如锆石生长域及其与U-Pb年龄的联系），该交互需要比本文所呈现的更高水平的自动化和分析技术（即使对于小锆石年龄组），鉴于需要检查以放置位点（可选分类）、比较成分和/或纹理差异及定义碎屑群体的图像数量庞大，这引入了潜在的偏差和局限性。最后，使用颗粒网格而非镶嵌图对分析数据质量和人类操作者性能的直接影响尚待基准和行为研究测定。将本概念验证工作开源化的努力有助于整合地质与技术领域的壁垒，并有望扩展至其他应用、矿物、仪器、数据科学方法及编程语言。