现代河流碎屑磷灰石原位地球化学分析——对长江物源示踪研究的启示

林旭，陈济鑫，吴中海，刘海金，张洋

1.三峡大学土木与建筑学院，湖北宜昌 443002

2.中国地质科学院地质力学研究所，北京 100081

3.东华理工大学地球科学学院，南昌 330032

0 引言

河流是塑造地表形态，串联地球岩石圈、水圈、大气圈物质循环的重要介质[1-4]。因而保存于河流中的碎屑沉积物（砾、砂、泥）对理解构造变形、气候变化和地表侵蚀作用之间的联系具有广泛意义[5-12]。“亚洲水塔”青藏高原与多个亚洲边缘海盆地及亚洲内陆湖泊组成了规模巨大的源—汇系统[13-14]。认识这一源—汇系统的发展过程，包括河流砾、砂、泥的产生、搬运与沉积过程，不仅为探索地球海陆交互作用等关键过程提供基础支撑，也能提供水系演化的信息[15-23]。其中尤以亚洲最大河流，发源于青藏高原东部，流经云贵高原、四川盆地、江汉盆地和长江中下游平原，最终流入太平洋西部陆架边缘海（东海）的长江最受关注[4,11,24]（图1）。而长江的发育和形成综合反映东亚中—新生代以来构造体制转换、气候变化等重要信息，对中国三大地势阶梯的建立具有重要指示意义，也是探索青藏高原构造—气候—沉积系统科学研究的重要组成部分[25-31]。因而，国内外研究者从环境磁学[11]、全岩地球化学[18-20,24]、重矿物组成[7,15,25,30]、黏土矿物组成[31]、单颗粒矿物年龄谱对比[8-10,32]等方面，对其源汇过程开展了广泛研究。然而，全岩地球化学组成更多代表沉积物均一化的物源信息；矿物组合研究，有时因流域内缺乏标型矿物，增加了物源解析的难度；锆石矿物硬度大，可以经历多期次物质沉积循环过程，仅仅依靠单一的U-Pb 年龄谱对比，并不能很好的限定潜在物源区。

近年来，在河流沉积物物源示踪研究中，越来越多的学者开始对某一矿物（例如锆石和磷灰石）同时进行原位物源信息（U-Pb年龄，Sr-Nd同位素等）的获取[33-39]，已被证明在界定物源区，识别构造活动引起的水系变化等方面十分有效。进而避免分析单一矿物引起物源示踪结果的不确定性。尽管磷灰石在重矿物中的含量十分低，在搬运过程中也容易出现磨蚀，在大型河流物源示踪研究中会存在代表性不强的现象。但在大河流域的个别流域，比如流域内存在磷灰石成分特别高的母岩，或者不同支流流域中磷灰石含量差异巨大等情况，可以起到很好的物源示踪效果[34-35]。分析河流下游或者沉积盆地中的碎屑磷灰石的原位地球化学元素组成，将其与上游河流或潜在源区进行对比，从而建立源汇关系。近年来，磷灰石这一物源示踪方法被广泛应用于河流物源示踪[40-42]、盆地古流向重建[36]等方面。但在长江流域这项工作还未开展。因而，我们对长江上游主要干支流开展碎屑磷灰石原位微量元素和Sr同位素分析，将金沙江、大渡河、岷江、嘉陵江支流与丰都段长江干流的结果进行对比，来探索这一方法在长江上游物源示踪研究的可行性。

图1 长江流域位置分布图Fig.1 Location of Yangtze River Basin

1 研究背景

长江发源于青藏高原的唐古拉山脉，在上海市崇明岛以东注入东海，连接了地球上最高的高原和最大的大洋，全长约6 300 km，流域面积约1.8×106km2，是亚洲最大河流。从源头到宜昌是其上游河段。伴随着晚新生代以来青藏高原的隆升，在亚洲东亚和南亚季风气候控制下，长江流域产生并携带巨量的水沙进入太平洋西岸陆架海，对流域生态环境和边缘海的海洋环境产生重大影响，在全球变化中扮演重要角色，是国内外众多重大研究计划的靶区[4,43-44]。

宜昌以上长江流域面积为1.0×106km2，1956—2007 年的年均径流量达到4.3×1011m3，年均输沙量为4.5×108t[45]。长江上游的四大支流金沙江（图2）、岷江（包含大渡河）、嘉陵江和乌江流域面积分别为5×105km2、2.2×105km2、1.6×105km2、0.87×105km2，总面积8.8×105km2，占长江上游流域面积的87.9%。金沙江年均输沙量为2.4×108t，占长江宜昌水文站的54.2%。岷江年输沙量为0.47×108t，占宜昌水文站的10.4%。嘉陵江年输沙量（1.1×108t）占宜昌站的24.1%。乌江流域是长江南岸最大支流，但主要流经云贵高原的喀斯特高山地区，与上述河流相比，河流流域面积和年均输沙量小（0.05×108t），占宜昌站的5.6%。因而乌江并不是宜昌以上长江干流的主要物质输送河流，本次研究未对其开展样品分析。

图2 长江上游干支流流域面积和1956—2007 年的平均输沙量柱状图（数据来源于文献[45]）Fig.2 Column chart of basin area and mean sediment discharge 1956-2007 in the upper reaches of the Yangtze River (data from reference [45])

2 样品来源和研究方法

2.1 样品采集

对宜昌以上长江干支流进行样品采集，每个点大约采集3～5 kg 样品（图1）。为了采样更具有代表性，每一个样品由多个采样点组成，同时避免人为污染源。金沙江样品在宜宾市上游金沙江干流河漫滩采集（28°44′36″N，104°35′50″E）。大渡河样品采自白岩坎镇河漫滩（29°32′46″N，103°44′18″E），岷江样品采自乐山市岷江下游河漫滩（29°37′44″N，103°45′11″E）。嘉陵江干流河漫滩样品采自施家梁镇（29°47′36″N，106°28′25″E）。长江上游干流样品采自丰都河漫滩（29°52′30″N，107°27′33″E）。上述样品主要采集砂级沉积物。

将野外采集回来的碎屑样品经重砂分析、磁性分选等一系列过程，将磷灰石分离出来，并在双目显微镜下进行人工挑选提纯。每个样品随机挑选＞300颗制成环氧树脂靶，并对靶片进行表面抛光处理。然后对所有样品进行背散射图像拍摄，圈定磷灰石颗粒的分析位置，避开包裹体和裂隙部位，提高分析精度。

2.2 磷灰石微量元素分析

磷灰石原位微量元素和Sr同位素分析在南京宏创地质勘查技术服务有限公司微区分析实验室使用激光剥蚀—电感耦合等离子体质谱仪（LA-ICPMS）完成。激光剥蚀平台采用Resolution SE 型193 nm深紫外激光剥蚀进样系统（Applied Spectra，美国），配备S155 型双体积样品池。质谱仪采用Agilent 7900型电感耦合等离子体质谱仪（Agilent，美国），采用束斑直径50 μm、剥蚀频率10 Hz、能量密度3.5 J/cm2、扫描速度3 μm/s 的激光参数剥蚀NIST 612，调节气流以获得高信号强度（238U 约为6×105 cps）、低氧化物产率（ThO/Th＜0.2%）。选用100 μm 束斑线扫NIST 610 对待测元素进行P/A 调谐。磷灰石样品固定在环氧树脂靶上，抛光后在超纯水中超声清洗，分析前用分析级甲醇擦拭样品表面。采用5 个激光脉冲对每个剥蚀区域进行预剥蚀（剥蚀深度～0.3 μm），以去除样品表面可能的污染。在束斑直径50 μm、剥蚀频率5 Hz、能量密度2 J/cm2的激光条件下分析样品。数据处理采用Iolite 程序[46]，以NIST 610 作为外标，43Ca 作为内标计算微量元素含量，内标值假定化学计量式。

采用L1、L2、L4、L4、H1、H2 和H3 法拉第杯组合接收Kr、Rb、Er、Yb和Sr相关质量数的信号。标准溶液NBS 987的87Sr/86Sr测定值为0.710 245±0.000 025，与TIMS 测定值0.710 236 在不确定度范围内一致[47]。测量前，测定并校正法拉第杯的电子噪音背景，以及放大器的接收增益。在束斑直径50 μm、剥蚀频率8 Hz、能量密度6 J/cm2的激光条件下分析样品。通常采集20 s的气体空白，35～40 s的信号区间，信号积分时间0.262 s。剔除88Sr 信号强度低于0.35 V 低信号点以及超出两倍标准偏差的异常点。使用86Sr/88Sr=0.119 4，按指数规律对其他同位素比值的仪器质量分馏进行在线校正。每隔10～15 个样品测试两次Durango磷灰石标样，监控仪器状态。

磷灰石微量元素含量可直接使用非矩阵多维标度MDS（Multidimensional Scaling）统计分析[48]，它是基于K-S统计方法，利用数据之间的相似性对数据进行成分分析，这种相似性分析早已被广泛应用到了碎屑锆石U-Pb年龄、磷灰石微量元素组成、重矿物数据分析中，获得了较好的效果。因此，我们采用该方法对宜昌以上长江干支流碎屑磷灰石的微量元素组成进行MDS 统计分析，进行磷灰石微量元素含量相似/向异性分析，从而判断各个支流的磷灰石组成与丰都段长江干流的磷灰石的相关关系。

3 实验结果

磷灰石颗粒的阴极发光图像如图3 所示。金沙江、大渡河、岷江、嘉陵江和丰都段长江干流的碎屑磷灰石颗粒以次棱角和棱角长柱状为主，说明这些磷灰石没有经历多期次的沉积循环过程，反映了各自流域内最近磷灰石颗粒的外形特征。磷灰石微量元素和Sr同位素分析结果见表1。可以看出，金沙江的磷灰石P 元素含量，尤其是平均值（178 486）与大渡河、岷江、嘉陵江和长江干流明显不同。金沙江的磷灰石Sr元素平均值含量为1 006×10-6，超过其他各个河流，而其稀土元素总量（∑REE）又是本次所有河流中最低的（506.5×10-6）。在87Sr/86Sr 比值结果中，来自龙门山的三条河流（大渡河、岷江和嘉陵江）要比金沙江和丰都长江干流的高。因而，总体来说金沙江下游的磷灰石原位微量元素组成与龙门山的河流和长江干流明显不同。

在P 元素和Sc 元素含量（图4a）、P 元素含量和87Sr/86Sr 比值二维散点图中（图4b），金沙江的分布区域与大渡河、岷江和嘉陵江的明显不同，而丰都段长江干流样品与大渡河、岷江和嘉陵江样品彼此间形成重叠区域。在∑REE和P元素含量二维散点图中（图4c），金沙江下游的磷灰石颗粒分布区域没有和其他流域重叠在一起。在Sr元素含量和87Sr/86Sr比值形成的二维散点图中（图4d），尽管各个流域的分布区域相对集中，但依然能看到金沙江形成特定的分布区域。

图3 磷灰石原位地球化学分析阴极发光图像Fig.3 In situ geochemical analysis of apatite cathodoluminescence images

表1 磷灰石原位地球化学数据统计表Table 1 Descriptive statistics of in situ geochemical data for apatite grains

在MDS分析结果中（图5），金沙江的磷灰石P元素含量与其他流域的磷灰石P 元素含量判别距离较远，而嘉陵江、岷江和大渡河的磷灰石P 元素含量与丰都段长江干流具有相近性。据此判断，长江上游干流的磷灰石主要来自发源于龙门山的河流。

图4 磷灰石地球化学分析结果二维散点图（a）P vs.Sc；（b）P vs. 87Sr/86Sr；（c）∑REE vs.P；（d）Sr vs. 87Sr/86SrFig.4 Two-dimensional scatter plots of (a) P and Sc; (b) P and87Sr/86Sr; (c) ∑REE and P;(d) Sr and87Sr/86Sr contents of fluvial detrital apatite grains

图5 磷灰石MDS 判别图Fig.5 MDS identification diagram of apatite grains

4 讨论

磷灰石富集P 和Sr 等不相容元素，其含量变化由全岩SiO2含量和所在熔体的分配系数控制，因而母岩的分异程度是磷灰石化学组成的主要控制因素[49]。此外，磷灰石中的87Rb经过β-衰变后变成87Sr，在剥蚀、搬运和沉积过程中最初的87Sr/86Sr 比值不发生变化[33]。而这一比值在不同地区、不同时代岩体有明显差异，这是将其运用于大河物源示踪研究的基础。稀土元素（REE）由15个镧系元素（La～Lu）组成，由于它们具有相似的物理和化学性质，在地质体中都倾向于成组出现，在自然界REE 的分配形成存在差异[50-51]，这一差异也被广泛用来开展沉积物物源示踪研究。因而，将磷灰石的P、Sr 和REE 等元素的含量或比值与其他元素做散点图的方法用于物源示踪研究最为普遍[49-51]。对比碎屑磷灰石原位地球化学组成发现，丰都段长江干流的磷灰石颗粒的P、Sr、Sc和∑REE 元素含量与大渡河、岷江和嘉陵江的样品存在重叠，而与金沙江下游的磷灰石颗粒分布区域明显不同（图4）。另外，由于P元素是磷灰石的特征元素，对其进行MDS 相似/向异性分析更具有代表性。金沙江下游磷灰石与丰都段长江干流磷灰石的P元素组成判别距离较远，而丰都段长江干流磷灰石P元素组成和嘉陵江、岷江和大渡河的P元素组成判别距离相近（图5）。这说明丰都段长江物质主要来自发源于龙门山的河流，而金沙江对其碎屑物质的贡献不大。但长江流域的沉积物以轻矿物为主，作为重矿物的磷灰石所占比例低，将其用于长江水系的物源示踪研究时，需要考虑这一影响因素。

图6 （a）研究区已有研究采样点分布图，石英10Be[26,52]，重矿物[27]，锆石U-Pb 年龄[53]；（b）金沙江年径流量和年输沙量变化图（据文献[45]修改）Fig.6 (a) Map of sampling points in the study area, quartz10Be[26,52], heavy minerals[27], zircon U-Pb age[54];(b) Variation of annual runoff and sediment discharge in Jinsha River (modified from reference [45])

从石英宇成核素10Be 物源示踪和地表侵蚀速率结果[26,52]来看（图6a），金沙江流域的侵蚀速率较低，长江上游干流物质主要来自嘉陵江、大渡河和岷江；Vezzoliet al.[27]的重矿物组合和主成分分析结果表明嘉陵江成为长江上游干流碎屑物质的主要输送河流。长江上游黏土矿物物源示踪结果显示，大渡河、嘉陵江、岷江对长江干流的物质贡献大[31]。上述结果与水文站观察到的输沙量结果吻合[45（]图6b）。尽管碎屑锆石U-Pb 年龄[53]和磁铁矿化学组成[54]物源示踪结果表明，长江上游干流的碎屑物质除了主要来自嘉陵江外，金沙江的作用也不可忽视。但考虑到随着金沙江中下游溪洛渡和向家坝水库的建设并投入使用，其输沙量减少至0.02×108t[55]，泥沙主体部分被拦截在这些梯级水库内部，导致其输沙量甚至已经低于乌江（0.05×108t）[45]，对长江干流的泥沙贡献率显著下降。Heet al.[53]分析的锆石U-Pb 年龄封闭温度高，具有强抗风化性，先前保存于河床的碎屑锆石后期经历多期次沉积再循环过程；王中波等[54]的采样时间早于金沙江下游溪洛渡和向家坝水库投入使用的时间（2012 年）。因而，长江干流依然能出现金沙江以前输入的物质信号。此外，尽管嘉陵江[56]、岷江[57]（包含大渡河）的输沙量由于人类活动、气候变化和修建水库显著下降，但下降的幅度低于金沙江，嘉陵江和岷江依然有相对充足的碎屑物质供给。嘉陵江及岷江进入丰都段长江的距离远小于金沙江，发源于龙门山的河流输入的碎屑物质“稀释”了金沙江的碎屑物质信号，导致丰都段长江干流现代河流物质主要受发源于龙门山河流的影响。因而，利用长江流域现代河流碎屑磷灰石原位地球化学组成，可以很好地区分龙门山和青藏高原东南缘的河流物质信号。利用这一方法在江汉盆地的新生代地层中识别金沙江的物质信号，可能具有很好的物源示踪效果，但该方法无法进一步有效区分发源于龙门山的各条河流的物质信号。因此，在今后利用该方法开展长江流域物源示踪研究时，结合其他方法可以达到精准建立长江源汇关系的研究目的。

5 结论

通过对长江上游干支流开展碎屑磷灰石原位微量元素和87Sr/86Sr 比值分析，在二维散点图中发现金沙江的分布区域与大渡河、岷江和嘉陵江的明显不同，而丰都段长江干流样品与大渡河、岷江和嘉陵江样品彼此间形成重叠区域；进一步结合MDS 判断图结果揭示，丰都段长江干流的物质主要来自龙门山。磷灰石原位地球化学分析方法可以有效区分发源于青藏高原东南的金沙江和发源于龙门山的大渡河、岷江和嘉陵江的物质。但其无法进一步有效区分发源于龙门山的各条河流的物质信号。建议将这一方法与其他物源示踪方法联合使用，可以有效建立长江的源汇关系。

致谢 衷心感谢审稿人的建议。