磷灰石裂变径迹技术对中国云南省西北部杨拉斑岩型铜矿成矿后低温历史的约束

《Journal of Geochemical Exploration》：Apatite fission-track constraints on the post-mineralization low-temperature history of the Yangla skarn Cu Deposit, NW Yunnan, China

　　• 短的磷灰石裂变径迹（9.1 μm）显示出Jiangbian时期的存在记录了APAZ区域的居住情况，而非冷却过程。• 重模型将Jiangbian时期的温度定为1.2亿至1.3亿摄氏度，并在1000万年前降至60摄氏度以下。• Lunong的径迹更长（11.9 μm），定义了一个部分退化的区域。• Equal-Dpar对之间的差异为3.8 μm和4700万年：动力学解释无法说明这种对比。• F4断层两侧的岩块在退化程度上存在1.5至2.5公里的差异；Lunong是更深的目标区域。

1. 引言
矽卡岩矿床是全球铜资源的重要组成部分，矿床形成与随后保存之间的耦合是经济地质学中的一个长期研究主题（Yang等人，2026a；Yang等人，2026b）。一个矿床是否能够作为经济目标持续存在，取决于其成矿后的热历史和退化过程，这些过程控制着矿体的垂直保存、叠加氧化和富集的程度，以及深层原生矿的可勘探性（Yuan等人，2020；Yang等人，2025a）。尽管人们越来越认识到成矿后热构造过程的重要性（Reiners和Brandon，2006；Yuan，2016；Xu等人，2021a），但对三江特提斯造山带内大型矽卡岩铜矿床的系统低温热年代学研究仍然很少。尽管已经积累了大量的关于成矿同期岩浆-热液过程的地质年代学和地球化学数据（Zhu等人，2015；Deng等人，2024），但对成矿后冷却路径、退化程度和保存状态的定量约束仍然不足（He等人，2023；Ou等人，2021；Yang等人，2019），这限制了保存状况的评估和深部勘探目标的预测（Li等人，2021；Du等人，2017）。
Yangla矽卡岩Cu矿床位于中国西南部滇西北部的金沙江缝合带中部，含有约1.3亿至1.5亿吨矿石，平均铜品位为1.03%（Zhu等人，2015）。该矿床位于花岗闪长岩和碳酸盐岩的接触带内，其构造受到缓倾逆断层和陡倾正断层共同控制（Du等人，2017）。主要的矿化阶段被严格限定在晚三叠世，钼矿ReOs年龄为2.318±1.3至2.333±3.0百万年，花岗闪长岩锆石UPb年龄为2.331±1.4至2.356±1.2百万年（Yang等人，2014；Zhu等人，2015），矿化流体的均一化温度为300–600摄氏度（Zhu等人，2015；Du等人，2019）。尽管已经系统地描述了成矿过程，包括矿化温度和压力条件以及矿化流体的组成和演化，但成矿后的冷却和退化历史仍然知之甚少。在Yangla，这一问题由于多阶段的热历史而更加复杂：除了约2.33亿年的Cu矿化之外，该矿床还记录了一个约3000万年后的晚期W事件（Li等人，2024）。这两个热事件之间大约2亿年的间隔表明，在中新生代期间，矿区经历了多次构造-热扰动。然而，每个热事件在多大程度上覆盖了原始矿体，以及这些事件是否在低温热年代学信号中有记录，目前尚未确定。此外，当前的采矿目标深度为200–800米，但成矿后的累积退化程度这一评估深层矿体保存状况的关键参数，仅通过地质推断得出，缺乏独立的热年代学约束。另外，该矿床位于青藏高原东南缘的构造和地貌活跃区域。区域性的中新世退化脉冲（Zhao等人，2009；Tian等人，2015；Cao等人，2019；He等人，2023）是否在矿床尺度上表现出来，以及NE走向的F4断层两侧是否存在差异性退化，这些问题都尚未解决（Liu-Zeng等人，2018；He等人，2023）。
磷灰石裂变径迹（AFT）热年代学对60至120摄氏度的温度范围非常敏感，由于其校准良好的部分退火动力学，是约束矿物沉积物低温冷却历史的理想工具（Ketcham等人，2007；Reiners和Brandon，2006；Tao等人，2023）。我们整合了从Jiangbian和Lunong矿区采集的七个花岗岩样品的AFT年龄和磷灰石成分数据，以回答三个问题：（1）Yangla矿化区是否可以识别出新近纪的低温记录，以及它与已发表的区域性冷却记录有何关联？（2）这两个矿区的对比AFT年龄模式具有何种地质意义？（3）这些观察结果如何有助于评估构造因素对矿体保存的影响？通过将区域热年代学背景与矿区内的观察结果联系起来，本研究探讨了空间年龄差异如何在差异性退化和矿化岩石体积保存的评估中发挥作用。

2. 地质背景
2.1. 区域地质
Yangla矽卡岩Cu矿床位于三江古特提斯构造域的中部，位于金沙江缝合带（JSZ）内，西侧为Zhongza块，东侧为Qamdo-Simao块（Wang等人，2024a；Kou等人，2026）（图1a）。它是中国西南部滇西北部多金属成矿带的重要组成部分（Hou等人，2009；Hou等人，2023）。JSZ保留了古特提斯洋向西俯冲并最终闭合的完整地质记录，并经历了从晚古生代到新生代的多阶段构造演化（Hou等人，2024）。这种演化包括了系列构造-岩浆事件，包括洋壳俯冲、弧-大陆碰撞、碰撞后的伸展以及随后的走滑再活化，这些事件共同控制了该地区中等到长英质侵入体的沉积和热液流体通道的构造（Zhao等人，2006；Deng等人，2025；Hou等人，2024；Wang等人，2024b；Wang等人，2025a）。

该地区的岩浆活动呈现出明确的两阶段特征，每个阶段对应一个不同的成矿系统。晚二叠世至晚三叠世的岩浆活动产生了一系列花岗闪长岩侵入体，包括Baimaxueshan（2.55–2.48亿年）、Ludian（2.34–2.14亿年）和Jiaren（约2.30亿年）侵入体，它们分布在大致N-S走向的Yangla断层东侧（He等人，2018；Wang等人，2024a）（图1b）。相比之下，始新世的岩浆侵入体（约4300–3300万年前）主要分布在金沙江断层和Ailaoshan-Red河走滑剪切带，并与区域性的大规模斑岩-矽卡岩铜-金-钼矿床有遗传关联，与主Yangla矿化阶段无直接关系（Xie等人，2022）。这种两阶段的岩浆模式反映了不同时间期的地质动力学环境和相应的成矿系统。该地区的多个大型矿床，包括Yulong、Yangla、Xuejiping、Pulang、Lanping和Beiya，都与连续的俯冲和碰撞阶段及其相关的岩浆-热液过程紧密相关。总体而言，这些矿床示例了在整个古特提斯洋闭合周期内发展的多阶段、多类型的成矿特征（Deng等人，2017；Liu等人，2025）。

2.2. Yangla矽卡岩Cu矿床
Yangla矽卡岩Cu矿床是一个在古特提斯金沙江分支俯冲及随后的闭合过程中形成的碰撞-造山接触变质系统（Li等人，2021）。该矿床由从北到南分布的七个矿块组成：Beiwu、Nilv、Jiangbian、Linong、Lunong、Tongjige和Jiaren，其中Linong和Lunong含有主要的经济矿体（图2）。Yangla矽卡岩Cu矿床形成于晚三叠世（约2.33–2.30亿年）的印支期花岗闪长岩岩浆-热液活动，矿化流体温度为300–600摄氏度，估计的沉积深度约为2–5公里（Zhu等人，2015；Du等人，2019）。这些矿体（KT1–KT6）在花岗闪长岩-碳酸盐岩接触带内形成了一系列层控到接触平行的透镜体（图3）。随后，该系统受到了晚三叠世片麻岩岩脉（约2.14–2.05亿年）和渐新世（约30–29百万年）WSb热液事件（流体温度约为150–280摄氏度）的叠加影响，这些事件受到NE走向断层的控制，并主要集中在Lunong块（Liu等人，2021；Li等人，2024；Deng等人，2025；Sun等人，2025）。这些跨越约2亿年的多阶段热事件可能被AFT系统记录下来，是本研究的重点；详细的矿床生成、蚀变分带和矿物组合描述见其他文献（Jia等人，2020；Li等人，2021；Deng等人，2024）。

2.3. 新生代构造、侵蚀和气候背景
Yangla矿床的保存和退化受到青藏高原东南缘新生代演化的控制。随着印度-亚洲碰撞的开始（约5500万年前；Tapponnier等人，2001），该地区首先经历了始新世–渐新世的左旋剪切作用（约4000–3000万年前）和晚始新世–渐新世的地壳缩短和逆冲作用，以Yulong逆冲带为代表（约2800–2000万年前；Cao等人，2019；Wu等人，2024；Wang等人，2025b；Zhao等人，2025a），随后是中中新世加速退化阶段，伴随着Ailaoshan-Red河剪切带的运动反转变为右旋走滑，金沙江断层系统的重新活化，以及Three Rivers峡谷的深切作用（Tian等人，2015；He等人，2023；Shen等人，2024；Zhao等人，2025b）。这种构造活动与中中新世气候最适宜期（约1700–1400万年前）的亚洲季风降水加剧相耦合，增强了金沙江、澜沧江和怒江地区的河流侵蚀和侵蚀效率（Nie等人，2018）。Yangla矿区位于这个构造和气候活跃的走廊内，两侧是N-S走向的金沙江和Yangla断层，以及NE走向的次级断层，这为下面的矿床规模冷却和差异性退化提供了侵蚀-构造背景。

3. 样本和分析方法
3.1. 采样策略和样本制备
为了系统地约束Yangla矽卡岩Cu矿床的成矿后热演化和退化历史，对来自两个主要矿区Jiangbian和Lunong的岩石样品进行了磷灰石裂变径迹（AFT）分析。从与Cu矿化空间相关的印支期侵入体中选择了七个新鲜的、粗粒度的含磷灰石样品，包括细粒花岗岩、花岗闪长岩和花岗斑岩。三个样品（JB14-1、JB14-2、JB14-3）来自Jiangbian矿区，海拔高度在2997–3287米之间，涵盖了所有已记录的侵入岩性。四个样品（LU14-1、LU14-2、LU14-3、LU14-7）来自Lunong矿区，海拔高度在3052–3359米之间，同样涵盖了完整的岩性谱。这两个样品系列位于NE走向的F4断层的两侧，海拔高度范围相当（分别为290米和307米），因此块之间的任何系统差异不能仅归因于海拔。野外采样优先选择新鲜、完整的基岩露头，排除了高度风化的区域、与断层相关的角砾岩体以及显示明显热液蚀变的间隔，以便所测量的年龄和径迹长度分布能够反映侵入体的热历史，而不是局部流体-岩石相互作用。
矿物分离在南京顺克检测技术服务有限公司按照标准协议进行。样品通过颚式破碎机进行破碎，并筛选至60目（大约0.25毫米）的粗粒度。通过连续的重液分离和磁选实现了矿物的富集。随后，在双筒显微镜下手工挑选出高质量的磷灰石颗粒，选择那些具有发育良好的晶面、高透明度且没有可见夹杂物和裂纹的样品，以确保裂变径迹的均匀保存。

3.2. 磷灰石裂变径迹分析
裂变径迹（AFT）测年是在自然资源部的活动构造与地质安全关键实验室、中国地质科学院岩力学研究所进行的，采用了LA-ICP-MS/FT分析技术（Yang等人，2023年）。选定的磷灰石颗粒被安装在玻璃载玻片上的环氧树脂中，磨平以暴露内部晶面，并抛光至约1微米的表面光洁度。抛光后的样品在5.5摩尔/升的HNO?溶液中于20摄氏度下浸泡20秒以显示自发性裂变径迹。蚀刻后，表面通过真空溅射镀膜法覆盖5-7纳米的金膜，以提高表面反射率并减少成像过程中的内部反射伪影（Gleadow等人，1986年）。

使用配备高分辨率数字成像的Zeiss Axio Imager M2m自动裂变径迹分析系统测量了自发性径迹密度和Dpar（平行于c轴的裂变径迹蚀刻坑直径）。如果抛光表面与晶体学c轴完全平行、径迹分布均匀且无裂纹，则选择这些颗粒进行分析。样品和颗粒坐标被系统记录下来，并拍摄高分辨率光学图像，以便在相同颗粒位置进行后续的激光烧蚀分析。

单个颗粒的238U浓度通过激光烧蚀-电感耦合等离子体质谱（LA-ICP-MS）技术测定，使用的是Coherent GeoLas HD 193纳米激光系统和Agilent 7900四极杆ICP-MS设备。分析参数优化如下：激光光斑直径32微米，重复频率5赫兹，能量密度5.0焦耳/平方厘米，烧蚀时间50秒。烧蚀物质通过高纯度氦气与少量氮气的混合气体输送到质谱仪中。监测的同位素有43Ca、238U、232Th和Pb。NIST-610和NIST-612标准参考玻璃作为外部校准标准，43Ca作为内部标准以校正仪器漂移和元素分馏效应。分析过程包括每八个样品分析使用一个Durango磷灰石参考材料和两个标准玻璃测量，以确保持续的分析质量控制（Pang等人，2019年）。

单个颗粒的AFT年龄是通过绝对（无zeta）LA-ICP-MS/FT校准（Vermeesch，2009年；Cogné等人，2020年）计算得出的，该方法将原位238U（NIST-610/612外部标准；43Ca内部标准）与自发性径迹密度相结合，而不依赖于外部探测器（Yang等人，2025b）。每八个未知样品分析一次的Durango磷灰石标准用于推导和验证特定会话的校准因子：重复分析得到的LA-ICP-MS/FT年龄为29.7-32.3百万年，与其参考年龄（31.4±0.5百万年）一致，并确认了会话间的准确性（Yang等人，2023年）。

在用于年龄确定的相同样品上，使用Yang等人（2023年）描述的程序测量了受限裂变径迹长度和Dpar。仅测量了Zeiss Axio Imager M2m系统下1000倍放大倍数下可见的双端都清晰的水平径迹内的受限径迹（TINT）；拒绝了径迹内的裂隙（TINCLE）和与表面相交的径迹。记录了每个径迹与晶体学c轴的角度，并按照Ketcham（2007年）的方法将长度投影到c轴上；所有逆向建模都使用了c轴投影的值，表1中报告了两个测量出的平均径迹长度。每个样品测量了23到89条受限径迹。

表1. 研究样品的磷灰石裂变径迹（AFT）分析结果。

样品ID 经度 纬度 海拔（米） 岩性 颗粒数（个/克） 每颗粒平均计数面积（平方厘米） ρ（×10^5厘米^-2） U（微克/克） P（χ^2%） 中心年龄±1σ（百万年） 长度（微米）（个） DPar
Jiangbian矿块 JB14–199°07′44″E 28°56′38″N 31 98 细粒花岗岩 32 245 3.47E-05 2.37 29.40 75% 16.2±2.1 8.1±1.8 (89) 1.83
JB14–299°07′45″E 28°56′05″N 29 97 辉长岩 30 270 5.46E-05 1.69 26.95 54% 12.9±1.6 10.4±1.5 (62) 1.99
JB14–399°07′47″E 28°56′26″N 32 87 花岗闪长岩 22 218 4.92E-05 1.96 23.55 21% 17.4±2.6 9.7±1.4 (23) 1.64
Lunong矿块 LU14–199°06′46″E 28°54′40″N 31 27 细粒花岗岩 17 32 2.25E-05 0.88 7.62 93% 26.8±9.4 10.5±1.4 (47) 2.93
LU14–299°06′59″E 28°54′39″N 32 11 花岗闪长岩 26 113 1.89E-05 2.58 12.70 11% 45±10 13.0±1.0 (52) 2.20
LU14–399°06′50″E 28°54′37″N 33 59 辉长岩 19 481.66E-05 1.40 5.37 22% 63±24 11.9±1.1 (56) 1.86

单个颗粒的AFT年龄是通过激光烧蚀-电感耦合等离子体质谱（LA-ICP-MS）技术结合Agilent 7900四极杆ICP-MS确定的。分析参数优化如下：激光光斑直径32微米，重复频率5赫兹，能量密度5.0焦耳/平方厘米，烧蚀时间50秒。烧蚀物质通过高纯度氦气与少量氮气的混合气体输送到质谱仪中。监测的同位素有43Ca、238U、232Th和Pb。NIST-610和NIST-612标准参考玻璃作为外部校准标准，43Ca作为内部标准以校正仪器漂移和元素分馏效应。分析序列包括每八个样品分析使用一个Durango磷灰石参考材料和两个标准玻璃测量，以确保持续的分析质量控制（Pang等人，2019年）。

单个颗粒的AFT年龄是通过绝对（无zeta）LA-ICP-MS/FT校准（Vermeesch，2009年；Cogné等人，2020年）计算的，该方法将原位238U（NIST-610/612外部标准；43Ca内部标准）与自发性径迹密度相结合，而不依赖于外部探测器（Yang等人，2025b）。每八个未知样品分析一次的Durango磷灰石标准用于推导和验证特定会话的校准因子：重复分析得到的LA-ICP-MS/FT年龄为29.7–32.3百万年，与其参考年龄（31.4±0.5百万年）一致，并确认了会话间的准确性（Yang等人，2023年）。

使用Yang等人（2023年）描述的程序，在用于年龄确定的相同样品上测量了受限裂变径迹长度和Dpar。仅测量了Zeiss Axio Imager M2m系统下1000倍放大倍数下可见的双端都清晰的水平径迹内的受限径迹（TINT）；拒绝了径迹内的裂隙（TINCLE）和与表面相交的径迹。记录了每个径迹与晶体学c轴的角度，并按照Ketcham（2007年）的方法将长度投影到c轴上；所有逆向建模都使用了c轴投影的值，表1中报告了两个测量出的平均径迹长度。每个样品测量了23到89条受限径迹。

表1. 研究样品的磷灰石裂变径迹（AFT）分析结果。

样品ID 经度 纬度 海拔（米） 岩性 颗粒数（个/克） 每颗粒平均计数面积（平方厘米） ρ（×10^5厘米^-2） U（微克/克） P（χ^2%） 中心年龄±1σ（百万年） 长度（微米）（个） DPar
Jiangbian矿块 JB14–199°07′44″E 28°56′38″N 31 98 细粒花岗岩 32 245 3.47E-05 2.37 29.40 75% 16.2±2.1 8.1±1.8 (89) 1.83
JB14–299°07′45″E 28°56′05″N 29 97 辉长岩 30 270 5.46E-05 1.69 26.95 54% 12.9±1.6 10.4±1.5 (62) 1.99
JB14–399°07′47″E 28°56′26″N 32 87 花岗闪长岩 22 218 4.92E-05 1.96 23.55 21% 17.4±2.6 9.7±1.4 (23) 1.64

LU14–199°06′46″E 28°54′40″N 31 27 细粒花岗岩 17 32 2.25E-05 0.88 7.62 93% 26.8±9.4 10.5±1.4 (47) 2.93
LU14–299°06′59″E 28°54′39″N 32 11 花岗闪长岩 26 113 1.89E-05 2.58 12.70 11% 45±10 13.0±1.0 (52) 2.20
LU14–399°06′50″E 28°54′37″N 33 59 辉长岩 19 481.66E-05 1.40 5.37 22% 63±24 11.9±1.1 (56) 1.86
LU14–799°06′35″E 28°54′30″N 30 52 花岗闪长岩 29 140 5.23E-05 0.98 14.49 30% 13.7±2.6 12.0±1.1 (31) 1.79

注：中心年龄和不确定性以±1σ报告，使用RadialPlotter（Vermeesch，2009年）根据LA-ICP-MS/FT协议（Yang等人，2023年）计算得出。ρ = 自发性裂变径迹表面密度（径迹/厘米^2）。Ns = 所有测定颗粒中计数的自发性径迹总数。长度（微米）（个）：未校正c轴投影的平均受限径迹长度，括号内为测量的径迹数量；“—”表示可测量的受限径迹少于三个，无法进行统计MTL确定。Dpar = 平行于晶体学c轴的平均蚀刻坑直径（微米），用作裂变径迹退火行为的动力学代理；在与自发性径迹计数相同的抛光颗粒表面上测量。

3.3. 逆向热历史建模
时间-温度（t-T）历史使用HeFTy v. 1.9.3（Ketcham，2018年）进行反演，输入为单个颗粒年龄分布、c轴投影的受限长度分布和样本平均Dpar，以及Ketcham（2018年）的多动力学退火模型，初始径迹长度为16.3微米。对于每个样品，在从140摄氏度开始的搜索空间内测试了100,000条随机生成的t-T路径，开始时间设定为超过该样品的最老单个颗粒年龄，结束时间为当前表面温度15±5摄氏度。没有强制要求在任何特定时间冷却，因此冷却的时间和速率是由数据驱动的。采用Kuiper统计量评估拟合优度（GOF）；根据Ketcham（2007年），GOF>0.05的路径被认为是“可接受的”，GOF>0.5的路径被认为是“良好的”。由于AFT反演不是唯一的，因此解释基于整个可接受拟合包络的共同特征，而不是任何单个最佳拟合路径。

4. 结果
4.1. AFT年龄和受限径迹长度分布
三个Jiangbian样品的中心年龄紧密聚集。样品JB14-1的年龄为16.2±2.1百万年（n=32颗粒；MSWD=0.82；P(χ^2)=0.75）（图4a），JB14-2的年龄为12.9±1.6百万年（n=30；MSWD=0.95；P(χ^2)=0.54）（图4c），JB14-3的年龄为17.4±2.6百万年（n=22；MSWD=1.2；P(χ^2)=0.21）（图4e）。所有三个样本的χ^2检验都在5%的水平上通过，因此每个样本内的单个颗粒年龄与单一群体一致（Galbraith等人，1990年；Vermeesch，2009年），尽管JB14-3中测定颗粒的数量较少，导致分辨率较低。有限混合模型解析出每个样本中的单一年龄组分（比例=100%），峰值年龄分别为16.2±2.1百万年、12.9±1.6百万年和17.9±2.4百万年，这些年龄在相应中心年龄的不确定性范围内；较小的差异反映了中心年龄背后的随机效应模型，该模型允许过度分散。这三个年龄在±1σ范围内重叠，跨度约为4.5百万年。单个颗粒的238U浓度在每个样本内的颗粒间有所不同，JB14-1中为19.5至140 ppm，JB14-2中为18.3至35.6 ppm，JB14-3中为16.2至30.4 ppm；这种颗粒间的异质性是岩浆磷灰石的典型特征，需要通过单个颗粒的U测定来进行无偏年龄计算。

下载：下载高分辨率图像（916KB）
下载：下载全尺寸图像

图4. Jiangbian矿块的磷灰石裂变径迹数据。（a, c, e）JB14–1、JB14–2和JB14–3的单个颗粒AFT年龄的径向图。径向尺度表示单个颗粒年龄（百万年），x轴表示精度（t/σ），y轴表示标准化估计值；点按单个颗粒的238U浓度着色。每个样本的中心年龄、测定颗粒数量（n）、MSWD和P(χ^2)以及通过有限混合模型解析出的峰值年龄分别给出。（b, d, f）相同样品的受限径迹长度分布。

相比之下，这三个样品的受限径迹长度分布显著缩短。平均径迹长度（MTL）分别为8.1±1.8微米（n=89；JB14-1；图4b），10.4±1.5微米（n=62；JB14-2；图4d）和9.7±1.4微米（n=23；JB14-3；图4f）。所有三个分布都是单峰且宽度较大，标准差为1.4至1.8微米，并且呈负偏态，JB14-1和JB14-3的表观年龄-海拔关系（AER）适用于Jiangbian和Lunongore区块的AFT中心年龄。水平条表示±1σ的中心年龄不确定性，虚线是对每个区块进行的普通最小二乘拟合（Jiangbian，n = 3；Lunong，n = 4）。有三个限制因素影响了这些梯度的解释：首先，Jiangbian的采样起伏仅为290米，Lunong为307米，这比通常认为能够产生有意义剥露率的垂直范围小一个数量级（Yuan等人，2007年）。其次，每个阵列仅包含三个或四个样本，因此决定系数的高度较低。第三，样本在每个区块内横向分布，而不是形成一个单一的垂直剖面，因此不能排除区块内结构的相对垂直偏移（van der Beek等人，2010年）。更根本的是，有限的轨迹长度数据表明，这些梯度不能直接解读为剥露率。年龄-海拔斜率只有在样本通过闭合窗口时迅速冷却的情况下才接近剥露率，即跨越起伏的时间相对于年龄范围较短。这两个区块都不满足这一条件：Jiangbian的平均轨迹长度为8.1至10.4微米，Lunong为10.5至13.0微米，记录了在磷灰石部分退火区内的长时间停留。在这种条件下，Jiangbian的年龄阵列被压缩，而Lunong的年龄阵列被大幅扩展，Lunong的0.006公里/百万年（相当于6米/百万年，持续49百万年）对于一个活跃的剥露造山边缘来说速度过慢，而更像是剥露部分退火区的特征。因此，这里的表观梯度被视为对这两个阵列的定性描述，定量冷却和剥露率是通过逆热历史模型得出的。

4.3. 磷灰石微量元素和稀土元素（REE）特征
对磷灰石进行了原位微量元素分析，特别关注了以球粒陨石为标准的稀土元素（REE）分布模式。Jiangbian矿区块的三个样本（JB14-1、JB14-2和JB14-3）的REE分布曲线分别显示在图9(a-c)中，Lunong矿区块的四个样本（LU14-1、LU14-2、LU14-3和LU14-7）的REE分布曲线显示在图10(a-d)中；每条曲线代表一个单点（或单颗粒）的分析结果。Jiangbian矿区块的磷灰石通常显示出一致的REE分布模式，不同分析点的曲线大致平行，表明晶体化学成分的整体均匀性。JB14-1、JB14-2和JB14-3样本在LaNd段显示出相对较高的浓度，在Eu位置有明显的负异常，并在整个GdLu段逐渐减少，点间分散主要表现为标准化值的大小变化，在重REE（GdLu）段尤为明显，特别是在JB14-2和JB14-3中。这些模式，特别是强烈的Eu负异常，是磷灰石从相对还原的斜长石分异花岗岩熔体中结晶的特征（Belousova等人，2002年；Chu等人，2009年），并不表明随后发生了成分变化。

图9. 以球粒陨石为标准的Jiangbian矿区块磷灰石稀土元素（REE）分布模式，标准化为Sun和McDonough（1989年）的CI球粒陨石值。(a–c) Jiangbian矿区块样本：(a) JB14–1, (b) JB14–2, (c) JB14–3。所有三个样本都显示出LREE富集、明显的Eu负异常以及系统性的HREE贫化，大致平行的曲线表明了统一的岩浆结晶环境。

图10. 以球粒陨石为标准的Lunong矿区块磷灰石稀土元素（REE）分布模式，标准化为Sun和McDonough（1989年）的CI球粒陨石值。(a) LU14–1; (b) LU14–2; (c) LU14–3; (d) LU14–7。样本LU14–1和LU14–3显示出相对规则的模式，具有负Eu异常，而LU14–2和LU14–7显示出强烈的不规则、锯齿状模式，具有LREE贫化和HREE富化，以及明显的元素间大小变化，这被解释为多阶段热液叠加的结果。

Lunong矿区块的磷灰石显示出更大的变异性，无论是在样本之间还是样本内部（Bau和Dulski, 1995; Irber, 1999）。LU14-1（图10a）显示出随着原子序数增加而单调减少的模式，并有明确的负Eu异常。样本LU14-3（图10c）显示出逐渐增加的轻稀土元素（LREE）段（LaSm）、明显的负Eu异常、在Gd处达到相对最大值，然后在重稀土元素（HREE）段向Lu方向系统性地减少。相比之下，LU14-2（图10b）和LU14-7（图10d）显示出La和LREE的贫化、HREE的富集，以及Ce和Lu之间不规则的段尺度曲率，相邻镧系元素之间存在显著差异，绝对REE浓度也有相当大的点间变化。这种向上凸起的 segmentation 是镧系元素四重效应的形态，在磷灰石中是由含氟和氯的水流相互作用产生的，而不是由岩浆分异引起的（Bau和Dulski, 1995; Irber, 1999）。因此，REE数据沿着与裂变径迹数据相同的边界划分了样本组。Jiangbian的磷灰石在成分上是均匀的，并保留了岩浆REE模式，与其狭窄的Dpar范围一致，表明它们来自单一的动力学群体。四个Lunong样本中有两个记录了流体介导的成分变化，这与记录中的Oligocene WSb热液事件仅限于Lunong矿区块是一致的（Li等人，2024年；Deng等人，2025年）。

4.4. 逆热历史模型
对七个样本中的五个样本进行了逆模型计算（图11）；JB14-3和LU14-7由于封闭轨迹太少（n = 23和31）而无法进行反演，因此被排除。在这五个样本中，JB14-2、LU14-2和LU14-3既再现了测量的年龄也再现了长度分布，而JB14-1和LU14-1虽然再现了长度，但不满足年龄标准。在每个可接受的范围内有两个共同特征：没有样本需要通过闭合窗口进行快速单调冷却，每个样本都需要在磷灰石部分退火区（APAZ）内或紧邻其上有一个较长的时间间隔；Jiangbian样本在任何时间都比Lunong样本热40到60°C。

图11. 使用Ketcham等人（2007年）的多动力学退火模型，通过HeFTy（Ketcham等人，2007年）得到的五个Yangla磷灰石样本的逆热历史模型。左侧面板（a, c, e, g, i）将测量的封闭轨迹长度分布（直方图，带有核密度估计）与模型预测进行比较；右侧面板（b, d, f, h, j）显示了相应的时间-温度路径。JB14-1的 length分布得到很好的再现（模型MTL = 9.1 ± 1.4 μm，长度GOF = 0.43），但没有满足年龄标准的路径（模型年龄19.3 Ma，年龄GOF = 0.02；图11a, b），模型高估了测量的中心年龄16.2 ± 2.1 Ma。JB14-2提供了最清晰的解决方案（模型年龄13.2 ± 0.8 Ma，年龄GOF = 0.60；模型MTL = 10.6 ± 1.8 μm，长度GOF = 0.85；图11c, d），其良好拟合的路径将其置于最近的Oligocene到最早的Miocene期间，之后需要连续冷却通过APAZ，并在约10 Ma时降至60°C以下；从大约13 Ma开始从高于120°C的快速冷却路径超出了良好拟合范围，因此其中心年龄12.9 ± 1.6 Ma无法记录通过闭合等温线的过程。LU14-1同样拟合了其长度（模型MTL = 10.5 ± 1.4 μm，长度GOF = 0.87），但未能拟合其年龄（21.5 ± 3.1对比26.8 ± 9.4 Ma，年龄GOF = 0.03；图11e, f），在这里低估了测量值。其余两个Lunong样本描述了系统上更冷和更浅的历史：LU14-2（模型年龄43.9 ± 2.6 Ma，年龄GOF = 0.63；模型MTL = 11.9 ± 1.2 μm，长度GOF = 0.83；图11g, h）从未超过大约80°C，并在大约55至80°C之间遵循近乎等温的路径，直到大约10 Ma后的最终冷却；而LU14-3（模型年龄57.9 ± 2.2 Ma，年龄GOF = 0.37；模型MTL = 11.9 ± 1.1 μm，长度GOF = 0.72；图11i, j）在大约70至57 Ma之间从大约125°C冷却到大约60°C，然后在Cenozoic的其余时间内保持在60至70°C，处于APAZ的上限。因此，区块之间的模型温度差异从Paleogene时期就存在，并不是由最后的冷却事件单独产生的。

JB14-1和LU14-1的失败是系统性的而非随机的，两个模型都没有产生良好的拟合范围。这两个模型都不用于地质解释，Jiangbian的历史仅基于JB14-2。因为AFT反演不是唯一的，并受到施加的搜索空间和约束条件的限制（最佳拟合路径包含的信息少于包围它们的范围）。JB14-1、JB14-2和LU14-1的最佳拟合路径在最终冷却前的大约5到3 Ma时，温度有短暂的上升（图11b, d, f）。在25到28°C/km的古老地热梯度下，这将需要大约2.5到3.5公里的上新世埋藏和重新揭露，而Yangla地区没有保存这样的地质或结构证据，在相同范围内存在质量相当的单调冷却路径；因此，这种上升被视为模型伪影。因此，以下解释依赖于整个可接受拟合范围共有的特征，即每个样本在给定时间的温度区间以及两个矿区块的相对顺序。

5. 讨论
5.1. 通过磷灰石裂变径迹热年代学记录的矿化后冷却
Yangla矽卡岩Cu矿床的主要矿化阶段被限制在晚三叠世（大约233 Ma；Yang等人，2014年；Zhu等人，2015年），矿形成流体的均化温度为300–600°C（Du等人，2019年），远高于AFT系统的有效闭合温度（Gallagher等人，1998年；Ketcham等人，1999年）。因此，AFT年龄不能记录初始的岩浆后冷却，而是记录了每个样本通过AFT闭合窗口的较晚时间（Tao等人，2017年，Tao等人，2019年）。在这个共同的起点上，两个矿区块记录了截然不同的信号：Jiangbian样本结合了紧密聚集的早到中Miocene年龄（17.4–12.9 Ma）和显著缩短的封闭轨迹，而Lunong样本结合了广泛分散的年龄（63–13.7 Ma）和系统较长的轨迹。因为表观年龄和轨迹长度都因退火而减少，这种对比不仅限制了每个区块冷却的时间，还限制了它从何温度冷却，逆模型量化了这一差异。综合地质年代学和AFT数据揭示了Yangla矽卡岩Cu矿床的多阶段矿化后热历史，从晚三叠世岩浆-热液矿化到Oligocene热液叠加，再到不同的Miocene冷却和剥露（图12）。

图12. 从已发表的地质年代学和本研究的AFT数据编译的Yangla矽卡岩Cu矿床矿化后热历史的温度-时间框架。y轴表示温度（°C）；x轴表示地质年龄（Ma）。水平条表示每个热事件或热年代学系统活动窗口的推断时间和持续时间；数据来源：Chen等人，1999年；Yang等人，2012年；Meng等人，2016年；Li等人，2024年；Deng等人，2025年。

晚三叠世到古近纪的演变（阶段I和II）构成了两个区块不同Cenozoic历史的共同前奏。阶段I（大约233–230 Ma）包括同期碰撞的岩浆活动和Cu矿化。阶段II（大约227–30 Ma）是一个从晚三叠世到古近纪的漫长后碰撞到后造山时期，在此期间，矿床经历了多次构造-热扰动，包括大约227 Ma的贫瘠石英斑岩、大约214–205 Ma的细粒岩脉、Yanshanian岩浆再活化，以及初期印度-亚洲碰撞的远场效应（大约55 Ma；Tapponnier等人，2001年）。尽管SE边缘在晚Eocene–Oligocene期间经历了缩短和推覆相关的剥露（Cao等人，2019年；Wu等人，2024年），但Jiangbian AFT系统在这个区间内没有记录到可分辨的快速冷却特征，表明这个区块在Miocene之前没有通过闭合；JB14-2的逆模型将样本在这个区间末期置于120–130°C，仍在闭合窗口以下，并位于磷灰石部分退火区（APAZ）内，因此短轨迹长度记录了在那里的长时间停留而不是后来的再加热事件。相反，古近纪的信号仅保存在较少重置的Lunong颗粒中，其模型将LU14-3在大约57 Ma时置于60°C以下，LU14-2在整个Cenozoic期间保持在80°C以下。

5.1.1. Jiangbian矿区块：延长的APAZ停留时间和持续的中Miocene冷却
三个Jiangbian样本产生了相对聚集的中心年龄（17.4–12.9 Ma），与共享的中Miocene低温冷却记录相符合。然而，它们的年龄一致性本身并不能证明快速冷却或排除更早的退火。具有相似轨迹保留行为的磷灰石群体，在不同的热停留、重热和随后冷却组合后，可能会产生相似的表观年龄。有限的轨迹长度数据消除了这种不确定性。平均轨迹长度分别为8.1 ± 1.8微米、10.4 ± 1.5微米和9.7 ± 1.4微米（样本数量分别为89、62和23；合并后为9.1微米），这些长度与未退火的磷灰石中诱导轨迹的初始长度（约16.3微米；Donelick等人，2005年）相比显著减小，且分布范围较广且呈负偏态。快速冷却通过闭合窗口后，轨迹长度变为13–15微米，分布呈现狭窄而对称的状态，因为这些轨迹是在较低温度下形成的，并且在随后的过程中很少缩短（Gleadow等人，1986年；Ketcham等人，2007年）。因此，观察到的分布与快速中新世冷却情况不符，而是表明大部分轨迹是在样品位于磷灰石部分退火区（APAZ）期间积累的，并在此过程中逐渐缩短的。JB14-2的逆模型量化了这一过程。最佳拟合路径表明，样品在最近的上新世到最早的中新世期间位于120–130°C，需要以约5°C/百万年的速度持续冷却，到约1000万年前温度降至60°C以下；而从约1300万年前开始的从高于120°C的快速冷却则不符合最佳拟合范围。因此，江边地区的中心年龄被下调，这是由于在较长的时间间隔内轨迹积累和退火共同作用的结果，而非单一冷却事件的结果。据此，第四阶段的定义被重新调整为约2200万至1000万年的时间区间，对应于样品通过APAZ的长期冷却过程，而不是1700万至1300万年间的短暂脉冲冷却（图13d）。这种紧密的集群结果与这一解释一致：由于江边地区的磷灰石形成了一个统一的动力学群体（Dpar值为1.64–1.99微米），并且所有颗粒的U和REE元素分布均匀，因此无论APAZ停留时间多长，均能得到连贯的单颗粒年龄。因此，年龄的连贯性限制了动力学的均匀性，而非冷却速率。大约3000万年前的W–Sb事件早于这些年龄，如下文中对吕农地区的讨论所示，该事件仅在该地区引发了局部时间重置，而没有影响江边地区的年龄记录，这与模型中JB14-2位于120–130°C的要求相符，低于该事件发生的150–280°C的流体温度。江边地区相对一致的REE元素模式并不排除更早的热扰动，正如变化较大的吕农地区模式并未显示出AFT时间重置一样。

下载：下载高分辨率图片（572KB）
下载：下载全尺寸图片

图13. 藏东南部高原边缘和三江地带已发表的低温热年代学年龄（AFT、ZFT、AHe、ZHe）的地域汇编。包含坐标、年龄和来源的完整数据集见补充表S1。缩写说明：缝合/地块边界：BNS，邦公-怒江缝合带；KLS，昆仑缝合带；CSSZ，昌宁-孟连缝合带；ARSZ，阿老山-红河剪切带；断层/剪切带：GYF，甘孜-玉树断层；XSHF，仙水河断层；LTF，理塘断层；ZDF，中甸断层；BTF，八塘断层；XJF，小江断层；DLSF，德钦-丽江-香格里拉断层；YJSZ，英江剪切带；GLSZ，高黎贡剪切带；NTHF，南汀河断层；逆冲带：LMSTB，龙門山逆冲带；MJF，岷江断层。局部热年代学剖面/地块：YYTB，沿岸源剖面/地块；JQTB，建川剖面/地块。河流：金沙江、澜沧江（湄公河）、怒江、雅龙江、大渡河。地块/区域：松潘-甘孜、强塘、拉萨、喜马拉雅、长江。

关于图例中颜色的解释，请参考本文的网络版本。

5.1.2. 吕农矿块：AFT年龄分散与多阶段热叠加
四个吕农样本的中心年龄分别为2680万±94万、4500万±1000万、6300万±240万和1370万±260万，覆盖的时间范围（约4900万年）远大于江边样本的约450万年。尽管如此，所有四个样本的MSWD（平均轨迹长度分散度）均在0.54–1.4的范围内，因此这种分散反映了样本间的真实差异，而非分析或统计误差；卡方检验排除了这四个样本具有共同年龄的可能性。单个颗粒年龄与Dpar值无关，这进一步排除了Dpar异质性作为分散原因的可能性，因为在吕农地区，Dpar值范围较宽（1.34–2.93微米），但单个颗粒的年龄并未显示与此相关的依赖性。另外两个观察结果表明，这种分散反映了不同的热路径，而非轨迹保留能力的差异。首先，吕农地区的受限轨迹长度系统性地长于江边地区的轨迹（10.5–13.0微米，合并后为11.9微米），表明尽管年龄较长，但累积退火作用较弱。其次，中心年龄随样本海拔高度的增加而单调增加，这是按顺序通过APAZ的样本所预期的关系，而不是随机叠加的群体所表现出的。因此，适用于LU14-2和LU14-3的逆模型都表明，这两个样本在新生代大部分时间内温度低于80°C；LU14-1的模型未通过年龄标准，未被采用；而LU14-7的年龄（1370万±260万）与江边样本的年龄在统计上无法区分。

长期以来，这种分散被认为与第三阶段（约3000万至2900万年前）的Oligocene WSb热液重新活动有关。流体包裹体微测温表明，形成辉石的流体温度约为150–280°C（Li等人，2024年），这一范围与AFT闭合窗口重叠且在局部超出该窗口。从表面上看，这样的温度可能足以重置AFT系统，但中新世的江边地区年龄数据表明该事件仅局限于特定区域，并非整个区域同时发生。有两个观察结果支持这一点：首先，WSb重新活动受到NE方向断层的控制，并集中在吕农地区（Li等人，2024年；Deng等人，2025年）；沿断层通道的加热效应会随距离减弱，因此高温并未扩散到江边地区。其次，这两个地区的情况正好体现了这一点：近端的吕农地区保留了分散的、部分重置的年龄和异常的REE元素分布，而远端的江边地区则保持了均匀的U、REE元素和Dpar值，以及一个集中的磷灰石群体。热历史模型对这一事件的范围进行了更严格的限制。LU14-3的中心年龄为6300万±240万，平均轨迹长度为11.9微米，其采用的模型要求整个Oligocene期间的温度在55–70°C之间（图11i, j）。如果温度持续升高到150–280°C，其轨迹会完全退火，年龄也会被重置；因此，保留的前Oligocene年龄表明WSb流体仅局限于吕农地区的特定结构通道中，并未提升整个矿块的温度。因此，Oligocene事件是针对个别样本的局部叠加效应，而非整个地区年龄分散的原因。

在颗粒尺度上，LU14-2和LU14-7中异常的REE元素分布反映了流体与岩石的相互作用，可能导致卤素成分和REE元素替代的局部变异，从而影响有效的轨迹保留行为（Carlson等人，1999年）。然而，这种机制处于次要地位，因为它通过相同的成分控制机制起作用，且在两个样本中均未观察到年龄与Dpar值的相关性。因此，吕农地区的年龄分散主要是由于样本间真实的热路径差异造成的，表现为从晚白垩世到中新世期间依次通过APAZ的过程，而磷灰石的成分异质性和局部流体叠加起到了次要调节作用。整个地块因此保持了部分退火区域的特征，而非单一冷却事件的结果。

5.2. 西藏高原东南缘构造演化控制下的差异性剥露
已发表的研究表明，西藏东南部的新生代剥露过程涉及构造变形和地表过程的空间变化，不能用单一的同步冷却轨迹来描述（Clark和Royden，2000年；Tapponnier等人，2001年；Clark等人，2004年；Lei等人，2022年；Wu等人，2024年）。Cao等人（2019年）将Oligocene至早期中新世的地形形成与逆冲作用联系起来，而Nie等人（2018年）提出增强的季风降水是中新世湄公河切割的驱动因素；对沿湄公河两侧山体的建模强调除了切割作用之外的不同岩体抬升历史（Ou等人，2021年），Shen等人（2024年）记录了1700万至1100万年前上游金沙江地区的快速剥露。在这种 variation 的区域背景下，Yangla地区的数据提出了两个问题：沉积物尺度的年龄记录是否反映了区域性的中新世加速现象，以及相距约300米高度差的两个地块是如何达到不同温度水平的。

5.2.1. 剥露强度与年龄-海拔梯度的限制
两个地块的表观年龄-海拔梯度（分别为0.064和0.006千米/百万年）不能直接视为剥露速率，因为这两个地块都没有快速通过共同的闭合等温线冷却，且都保留了表明长期停留在APAZ的轨迹长度。因此，剥露强度是根据逆模型得出的。最佳拟合路径表明，JB14-2在约2200万年前位于120–130°C，相对于当时的地表低3.7–4.6千米（地热梯度为25–28°C/千米；He等人，2023年），到约1000万年前温度降至60°C以下（低约1.8千米）。因此，江边地区自早中新世以来已剥露了约4千米，平均剥露速率为0.2千米/百万年，这与西藏东部边缘的新生代剥露速度（Clark等人，2006年；Liu-Zeng等人，2018年；Shen等人，2021年）一致，是表观年龄-海拔梯度的三倍。在同一时间区间内，公认的吕农模型认为LU14-2的温度低于80°C，LU14-3的温度为55–70°C，相当于1.5–2.6千米。因此，这两个地块的剥露深度相差约1.5–2.5千米，这一差异仅凭年龄数据无法确定。

已发表的AFT、ZFT、AHe和ZHe年龄数据汇编（图13；补充表S1）将这种对比置于区域背景中：最年轻的年龄（<1500万）集中在切割明显的三江峡谷和主要断层沿线，而古近纪至中生代的年龄则保存在剥露较少的内陆地区。江边地区的年龄属于主要的15–800万年龄群体，而分散的吕农地区年龄属于次要的古近纪群体，因此这两个地块涵盖了整个区域范围，它们的对比反映了整个地带的剥露梯度。

5.2.2. 逆冲作用与河流切割的相互作用
岩石抬升、地表抬升和剥露是相互独立的因素（England和Molnar，1990年），在这里区分它们至关重要，因为Yangla地区的两个主要成因——古近纪缩短作用和新近纪切割作用——作用于不同的平衡机制。问题不在于哪个过程在构造上更活跃，而在于哪个过程 responsible for 去除了江边样本冷却过程中经过的约4千米的覆盖层。藏东南缘地区的古高程数据使得这种区分成为可能，因为这些数据显示地表抬升和中新世剥露在时间上是解耦的。来自Yangla矿区周围新生代盆地的稳定同位素和古植物学证据表明，该地区在晚Eocene至早Oligocene期间已经达到了接近现代的高度：建川地区约2.6千米（约4000万年前），Liming地区约2.7千米（约3700万年前）（Li等人，2015年；Hoke等人，2014年），以及Markam地区在晚Eocene至早Oligocene期间从约2.9千米上升到约3.9千米（Su等人，2019年）。尽管个别古高程数据的不确定性较大，且一些重建认为Eocene时期的抬升较低，主要抬升发生在Neogene时期（Gourbet等人，2017年；Wang等人，2018年；Botsyun等人，2019年），但整体来看，这些数据一致支持MIocene之前的地形升高。地壳增厚和逆冲作用相关的岩石抬升（以玉龙逆冲带为例，约28–20百万年前；Cao等人，2019年；Wu等人，2024年）显著早于江边地区的冷却期（约2200万至1000万年前），并在样本开始离开APAZ之前就已经形成了当前的地形格局。

这种变化是由河流作用引起的。由于快速河流切割的开始本身就是基底下降和侵蚀加速的标志（Clark等人，2005年；Ouimet等人，2010年），三江峡谷的低温度热年代学记录提供了相关证据。金沙江及其支流记录了约3000万、2200万、1800万和1500万年前的多次切割事件，其中最近的一次在上游表现为1700万至1100万年的快速剥露（Shen等人，2024年），而澜沧江、大渡江和雅龙江系统在这一时期的切割几乎是同步的（Zhao等人，2026年）。中中新世气候最适宜期（约1700万至1400万年）的季风加强进一步加剧了这些峡谷的侵蚀作用（Nie等人，2018年）。江边的年龄范围与大约1500万年的金沙脉冲一致，且无法区分。在我们的区域汇编中，最年轻的冷却年龄追踪的是切割的峡谷，而不是推覆带。然而，仅靠切割本身无法解释大约5000万年的年龄差异、平均轨迹长度2.7微米的差异以及在大约300米的地形起伏和单一的东北向构造中1.5-2.5公里的剥露深度差异。因此，F4上的阶梯状地形记录了叠加在区域侵蚀信号上的差异性垂直滑移。我们认为，古近纪的缩短作用是地形形成的前提条件，中中新世金沙河的切割是矿体剥露的直接原因，而F4上的逆向滑移则是将这种剥露作用分配到江边和龙农两大块之间的机制。

5.3. 矿床保存方式及其对深部勘探的启示
流体包裹体微热测年法将杨拉地区原生矿化流体的捕获压力限制在1.32-2.10千巴之间，均质化温度为373-409摄氏度（Yang等人，2014年；Zhu等人，2015年）。假设矿化流体的压力介于岩石静压（约26.5兆帕·千米^-1）和静水压力（约9.8兆帕·千米^-1）之间的梯度，矿化的古深度估计为2-5千米，这符合中等到浅层的矽卡岩型铜矿床的特征。这一深度被用作剥露预算的上限（Meinert等人，2005年）。

5.3.1. 矿化后的冷却和剥露模式及其区域类比
因此，杨拉矿床经历了两阶段的保存过程：从晚三叠世到古近纪，矿床在APAZ（矿化作用带）内长期停留，净剥露有限；随后在新生代出现了差异性剥露。与三江地块和西藏东部的类似岩浆-热液型铜矿床相比，这一过程具有独特性。根据本文提供的AFT（后热液事件）地质年代数据和区域地质证据，杨拉矿床的矿化后热演化可以被描述为一个两阶段的保存模型：首先是在矿化深度附近长期停留，净剥露有限；随后在新生代出现差异性剥露。在主要铜矿化事件（约2.33-2.30百万年前）之后，矿体进入了后碰撞期到后造山期，这一时期持续了大约2亿年，并伴随着印支期花岗闪长岩的侵入。在这段长时间内，矿床经历了晚三叠世的片岩脉形成（2.14-2.05百万年前；Meng等人，2016年）和渐新世的局部WSb热液再活化（30-29百万年前；Li等人，2024年），但江边地块并未通过AFT闭合窗口，其温度一直保持在APAZ范围（60-120摄氏度）内（Reiners和Brandon，2006年）。与三江地区和西藏东部的典型岩浆-热液型铜矿床相比，杨拉矿床的保存路径显得非常独特。例如，玉龙斑岩型CuMo矿床（约4100万年前）的锆石和磷灰石（UTh/He）年龄分别为37.5百万年和33.4百万年，表明该矿床在矿化后的5-8百万年内从岩浆冷却状态迅速转变为接近地表的条件，这与大约33-30百万年前西藏高原东缘的快速抬升事件相吻合（Li等人，2012年）。其保存方式属于“近同期剥露型”，即矿体在相对较短的时间内暴露于地表，并经历了显著的氧化和表生富集作用（Zhao等人，2025年）。岗底斯地区的曲龙斑岩型Cu矿床则表现出较慢的剥露响应（约1360-1240万年前；Hou等人，2009年；Dai等人，2013年），矿化后的剥露距离仅为1.3-2.8千米，属于典型的“晚期有限剥露型”，其中深部的原生矿体保存良好。熊村斑岩型Cu矿床（侏罗纪矿化）在20-13百万年内经历了4-6千米的快速剥露，与杨拉地区的江边矿体剥露幅度相当（Zhou等人，2019年）。相比之下，杨拉矽卡岩型Cu矿床的矿化年龄最早（约2.33百万年），但其AFT冷却年龄主要集中在中中新世（约1560万年前）；这个约2.15亿年的时间间隔是该地区同类矿床中记录最长的，反映了矿体在矿化后长时间停留于矿化深度附近，而在新生代之前的净剥露很小。

5.3.2. 剥露预算和原生矿体的保存深度
两个独立的深度约束提供了江边地块矿化后剥露预算的依据：一个是约233百万年时的流体包裹体矿化古深度为2-5千米；另一个是约22百万年时的模型计算深度为3.7-4.6千米。由于两者重叠，说明该地块在模型计算的冷却期开始时已经回到了或接近其矿化深度，其整个剥露过程都发生在早中新世之后。我们强调，这仅限定了约22百万年时的深度及之后的剥露过程；并不要求整个2.33-30百万年间的覆盖层厚度保持恒定，因为AFT系统对其闭合等温线以上的埋藏-剥露路径不敏感。因此，不能排除其间有更深层次的埋藏随后恢复到矿化深度的情况，这会增加累积剥露出地表的距离，而不仅仅是净剥露距离。

龙农地块的剥露过程较为温和。公认的模型表明，LU14-2和LU14-3自古近纪以来位于最上层APAZ内或之上，距离为1.5-2.6千米，因此该地块已经接近其当前的构造高度，而江边地块此时仍深约4千米。约1.5-2.5千米的剥露深度差异决定了两大块原生矿体的相对保存状况。这些数据与区域记录一致：Shen等人（2024年）记录了金沙河断层带熊松段在17-11百万年间的剥露距离为4千米；Liu-Zeng等人（2018年）报告了紧邻杨拉的德钦地区从中新世至今的剥露距离为3-5千米。所有深度转换都假设了恒定的地热梯度和几乎不变的地形变化，但在这种条件下这种状况不太可能持续2000万年。

这种演化过程在图14a-d的概念模型中得到了总结。矿床形成于晚三叠世（约2.33-2.30百万年），处于同碰撞构造环境中，随后从侏罗纪到古近纪一直位于APAZ内，在此期间裂变径迹经历了持续的部分退火和缩短；这一时期的覆盖层厚度并未被AFT系统确定，可能存在变化。剥露过程在始新世-渐新世期间开始，在印度-亚洲碰撞的远场作用下进行，最终在中新世达到了构造抬升和金沙河快速切割的顶峰，此时江边样品在约22-10百万年间通过了APAZ，记录了现代西藏高原东南缘高山峡谷地貌的热特征。

图14. 杨拉矽卡岩型Cu矿床矿化后历史的概念框架。(a) 晚三叠世与花岗岩侵入和矽卡岩形成的矿化。(b) 从晚三叠世到古近世，在磷灰石部分退火区（APAZ，60-120摄氏度）内长期停留。(c) 渐新世期间，龙农地块内沿东北向通道的WSb热液再活化。(d) 从早中新世到现在，F4断层两侧的差异性剥露，江边地块自约22百万年以来已剥露了约4千米，而龙农地块仅剥露了1.5-2.6千米。

5.3.3. 对深部勘探的启示
上述差异性剥露为深部勘探提供了定量空间框架：东北向的F4断层将两个从不同深度剥露出来的地块分隔开来，这两块地块的剥露深度相差约1.5-2.5千米。类似的断层控制分异现象也在艾劳山-红河剪切带（Wang等人，2025a）和丽江-大竹断层系统（Cao等人，2019）中有记录。由此可以得出以下三点启示：
(1) 江边矿体：经过约4千米的剥露后，原本位于2-5千米深度的矿体现在暴露于接近或低于矿化系统的原始基底的位置，因此上部的氧化和表生富集部分大部分已被移除（Cooke等人，2005年）。勘探应优先考虑矿化接触带的侧向和向下延伸部分以及隐蔽的矽卡岩段，而不是额外的垂直延伸部分，因为后者可能有限。
(2) 龙农矿体：仅经历了约1.5-2.6千米的剥露，因此暴露出了相同矿化系统的较浅层次，因此在当前地表之下仍保留了大量的原生矽卡岩矿化。主要成分为黄铜矿-磁铁矿-黄铁矿的组合基本未发生变化，该地块在当前采矿深度（200-800米）以下的1000-2000米深度内仍具有良好的原生矽卡岩矿化潜力。
(3) F4断层系统：自印支期矿化以来，该断层一直作为构造-热液通道发挥作用，并在中新世的WSb热液再活化过程中起到了关键结构控制作用（Li等人，2024年）。作为划分两大块矿体不同剥露程度的边界，它标志着氧化矿体和保存良好矿体之间的转变，其向下和沿走向的延伸部分值得优先勘探。通过关联其独特的气层单元、侵入边界和矿化层位，可以测试其位移是否重复、偏移或截断了勘探目标。

6. 结论
(1) 三个江边样品的年龄紧密聚集（17.4-12.9百万年），且退火轨迹较长（平均MTL 9.1微米，n=174），这与快速冷却过程不符。对JB14-2的逆向建模（年龄GOF=0.60）表明，该地块在约22百万年时温度为120-130摄氏度，并在约10百万年后降至60摄氏度以下，因此这些是记录了在磷灰石部分退火区内长时间停留的混合年龄。
(2) 四个龙农样品的年龄分散（63-13.7百万年），但退火轨迹较长（平均MTL 11.9微米，n=186），年龄随海拔升高而单调递增。接受的模型要求矿体在晚白垩世到古近纪期间通过了部分退火区，并随后保持在80摄氏度以下，因此该地块保留了一个剥露后的部分退火区轮廓；渐新世的WSb流体必须被限制在离散的通道内，因为LU14-3保留了前渐新世的年龄。
(3) 退火动力学无法解释这种年龄差异：单个颗粒的年龄与两块样品的Dpar无关，JB14-1和LU14-3的Dpar虽然相似，但平均轨迹长度相差3.8微米，年龄相差47百万年。因此，明显的年龄-海拔梯度（0.064和0.006千米百万年^-1）描述的是两组数据的特征，而非剥露速率。
(4) 模型显示，江边地块在约22百万年时的深度为3.7-4.6千米，龙农地块为1.5-2.6千米，这两种地块在300米的地形起伏和F4断层上的距离相差约1.5-2.5千米。基本上，所有江边地块的剥露都发生在早中新世之后，与17-11百万年金沙河的上部切割同时发生。矿化与最终冷却之间的约2.15亿年间隔是三江地区记录中最长的，而剥露较少的龙农地块在当前采矿深度（200-800米）之下保留了更大的原生矽卡岩矿体范围。

CRedI作者贡献声明
Li Yang：撰写——原始草稿、方法学、调查、正式分析、数据管理、概念化。
Zhidan Zhao：可视化、监督、资金获取、正式分析。
Jing Yang：监督、软件应用、调查、资金获取、正式分析、数据管理、概念化。
Jing Meng：调查、资金获取、正式分析、数据管理。
Wanming Yuan：可视化、验证、监督、数据管理。