结合主成分分析法(PCA)和正定矩阵因子分解法(PMF)的鄱阳湖丰水期表层沉积物重金属源解析*

匡荟芬，胡春华，吴根林，陈 苗

(南昌大学资源环境与化工学院，鄱阳湖环境与资源利用教育部重点实验室，南昌 330031)

重金属是典型的累积型污染物，其具有不可降解性、显著生物毒性和持久性特征，对环境质量构成了潜在的生态风险[1]. 随着我国工业化、城镇化和农业现代化的发展及污染物排放量的增加，沉积物重金属污染已经成为湖泊等水体最为突出的环境问题之一[2-3]. 对于湖泊流域生态环境保护，识别沉积物中重金属的污染源，并定量确定其相对贡献率变得更加重要. 现阶段源解析方法应用最为广泛的是以污染区域为研究对象的受体模型法[4]，模型通过受体和污染源样品的化学组分分析来确定污染源对受体的贡献值，常用的受体模型包括化学质量平衡法(CMB)、主成分分析/因子分析(PCA/FA)、PMF(正定矩阵因子分解法)等[4-7].

其中，PCA是开展污染物源解析研究应用较早且最广泛的方法，它优于CMB的主要原因是不需要通常很难获得的源配置文件，运行较简便，且CMB的应用大多针对于大气污染物[8]；部分学者虽将PCA应用于重金属贡献率研究中，比如林承奇等[9]将PCA应用于九龙江近岸表层沉积物重金属源解析中，得出3个污染源：农业生产活动和采矿活动的复合源、化石燃料燃烧源和岩石风化来源，贡献率分别为27.3%、31.8%和13.4%，但源解析结果总是不可避免地得到负值导致物理意义上的损失. PMF是美国环境保护局推荐的一种定量源解析模型，它的优点是不需要详细的源成分谱信息，可较好处理缺失及不精确的数据，对因子分解矩阵进行非负约束，使得到的源成分谱和源贡献率不会出现负值，结果更加合理[10-13]，已被广泛应用至大气、水体和沉积物等的污染物的源解析[14-16]，在沉积物重金属方面国内研究甚少；但它要求受体样本数据量大，且事先要设定因子数，多次运行程序，会对结果产生一定的误差，因此，确定因子数至关重要[17-18]. 董騄睿等[4]通过PMF法探究南京城郊农田土壤重金属的主要来源，得出4个污染源：农业源、降尘源、工业源和自然源，贡献率分别为30.8%、33.0%、25.4%和10.8%，并验证了该方法的可行性. 而冯小琼等[5]在研究香港地区2005-2010年VOCs污染来源解析时建议采用两种或以上的模型对受体点数据进行污染来源解析，以确保结果的准确性.

鄱阳湖流域内矿产资源丰富，主要有饶河上游的德兴铜矿、信江流域的永平铜矿以及赣江上游的大余钨矿等，长期开采与冶炼产生的“三废”不仅对矿区生态环境产生了影响，而且污染物通过河流汇入鄱阳湖，鄱阳湖生态环境也因此易遭受破坏，其重金属污染问题一直是人们关注的焦点. 以往关于鄱阳湖沉积物重金属元素的研究集中在鄱阳湖某一区域或“五河”入湖口几个采样点4种(Cu、Zn、Cd和Pb)或最多8种(Cr、Ni、Cu、Zn、Cd、Pb、Hg和As)重金属的时空分布和风险评价方面，如伍恒赟等[19]对鄱阳湖沉积物7种重金属(Cd、Hg、As、Cu、Pb、Cr和Zn)含量进行空间分布特征分析并评价其潜在生态风险，且通过重金属之间的相关性分析得出，Hg、Cu、Pb、Zn 4种金属元素具有同源性，而Cr与这几种元素来源不同，但与Cd、As来源相似，Cu、Hg、Pb等重金属主要来自乐安河流域工业排放；杨期勇等[20]对鄱阳湖北部湖区沉积物8种重金属(Cd、Zn、Cu、Cr、Ni、Pb、As和Hg)含量和空间分布特征进行分析并评价其潜在生态风险；而对鄱阳湖表层沉积物重金属污染来源及其相对贡献率尚不清楚. 因此，本文通过采集鄱阳湖丰水期表层沉积物样品来减小不同季节采样的差异从而尽可能降低源解析结果的不确定性，基于重金属元素V、Cr、Co、Ni、Cu、Zn、Sr、Mo、Cd、Sb、W、Pb、Hg和As的含量数据，分析重金属的污染特征，采用PCA定性地推测潜在的污染源，PMF识别重金属的污染源并量化其相对贡献，以期为鄱阳湖沉积物重金属污染防治及水体生态环境保护等提供科学依据.

1 材料与方法

1.1 研究区域

鄱阳湖(28°24′～29°46′N, 115°49′～116°46′E)是中国第一大淡水湖，位于江西省北部、长江中下游南岸，也是中国第二大湖. 湖泊承纳江西省境内赣江、抚河、信江、饶河和修水五大水系来水，经调蓄后，通过北部的狭窄出口进入长江，为长江流域一个重要的季节性吞吐型的浅水湖泊；水域面积在夏季丰水期达到3700 km2以上，而在秋冬季枯水期缩减至1000 km2以下. 鄱阳湖平原地区是中国重要的商品粮基地，江西省是中国建国以来从未间断向国家贡献粮食的两个省份之一，也是中国矿产资源大省(中国主要有色、稀有金属、稀土矿产基地之一，占全国有色金属产业的1/10；江西省重金属国有控股企业中有色冶炼和有色金属矿采选业占总数的63%[21]). 乐安河(饶河南支)中下游的德兴铜矿是亚洲最大的露天铜矿，乐安河流域内德兴市还有其他矿产资源，如朱砂红/富家坞铜钼矿、银山铅锌矿和花桥金矿矿区等，另建有银山铅锌矿、金矿冶炼厂等有色金属加工企业和建材化工等行业；赣江上游流经赣南有色金属采矿区、大余钨矿等矿产资源丰富的地区，中游流经吉泰盆地水稻种植区，下游流经南昌市等城市[22]；信江中游有永平铜矿和贵溪冶炼厂，下游流经余干县后注入鄱阳湖；抚河上游有铀矿，下游流经南昌县等. 鄱阳湖流域主要的铜矿、铅锌矿、钨矿等分布详见图1.

1.2 样品采集与分析

于2018年5-7月从鄱阳湖区使用抓斗式采样器采集了88个表层沉积物样品(0～10 cm)(图1)，将采集的沉积物样品装入干净的可密封聚乙烯袋中，注明采样地点和样品编号等，于0～4℃下保存，并运回实验室进行预处理. 冷冻后的沉积物样品置于干燥、通风、阴凉处自然风干，剔除样品中的残根、贝壳等杂物，用木棒研压，于玛瑙研钵中研磨后过100目(孔径0.149 mm)尼龙筛，筛下样装入干净的可密封聚乙烯袋中，-4℃ 保存于冰箱备用.

图1 鄱阳湖表层沉积物采样点分布

沉积物样品使用硝酸-氢氟酸-高氯酸(HNO3-HF-HClO4)的方法消解[23]，所有样品均在核工业北京地质研究院分析测试中心分析，采用德国Finnigan MAT公司生产的ELEMENT XR型高分辨电感耦合等离子体质谱仪(HR-ICP-MS)测定试液中V、Cr、Co、Ni、Cu、Zn、Sr、Mo、Cd、Sb、W、Zr和Pb含量. 在通过原子荧光法(AFS)测定沉积物中的Hg含量之前，首先用王水(HNO3∶HCl=1∶3)水浴消解[24]；沉积物中的As采用《土壤质量 总汞、总砷、总铅的测定 原子荧光法 第2部分：土壤中总砷的测定》(GB/T 22105.2-2008). 为保证分析结果的准确性，分析过程中平均每10个样品随机抽取1个样品作1个平行样(重复做3次实验)，同时做空白实验并以水系沉积物成分分析标准物质(GBW07309)进行质量控制. 各重金属元素测定结果的相对标准偏差均小于5%，符合国家规定的精密度要求. 所有样品的15种重金属含量均有检出.

1.3 重金属污染评价与源解析方法

1.3.1 富集系数法 重金属易富集于细粒沉积物中，受到“粒度效应”的影响，而富集系数法可以有效校正沉积物粒度和矿物组成变化对重金属含量的影响. 一般用于重金属标准化的元素有：Al、Fe、Mn、Sc、Ti和Zr等，本文选取惰性元素Zr作为参比元素，采用富集系数法区分重金属的自然与人为污染并评判其污染程度. 富集系数EF的计算公式为：

EF=(M/X)样品/(M/X)背景

(1)

式中，(M/X)样品是目标沉积物中重金属元素与参比元素含量的比值，(M/X)背景是背景样品中重金属元素与参比元素含量的比值. 理论上，EF>1，即存在人为污染. 根据Sutherland[25]的划分标准：EF<2，为无污染到弱污染；2≤EF<5，为中等污染；5≤EF<20，为重污染；EF≥20，为极重污染.

1.3.2 主成分分析法和正定矩阵因子分解法 源解析采用主成分分析法(PCA)和正定矩阵因子分解法(PMF)，PCA采用数学降维或特征提取方法，将原来众多的且具有一定相关性的原始变量进行线性变换，提取出数目较少且彼此间互不相关的重要变量，使用较少的有代表性的因子来解释众多变量的主要信息，并推测有关污染源的信息[26]. PMF是由Paatero等[27]提出的一种改进的用于源分配的因子分析受体模型. PMF将矩阵X(浓度矩阵)分解为因子矩阵G(或源贡献矩阵)、矩阵F(源组合物矩阵)和残差矩阵E，它的基本原理基于以下等式：X=GF+E. 其中矩阵F定义为与m个化学同源物有关的p源的化学组成矩阵，G矩阵描述了p源(因子)对n个取样点的贡献[28]. 根据USEPA PMF 5.0用户指南，公式如下所示：

(2)

1.4 数据处理

用Excel 2010和Origin 2018软件对重金属含量数据进行分析处理以及作图，计算了诸如平均值、最小值、最大值，标准偏差和变异系数之类的描述性数据；使用SPSS 17.0对鄱阳湖表层沉积物中重金属含量进行Pearson相关性分析和主成分分析(PCA)；根据记录的采样点坐标位置，结合每个对应点沉积物中重金属的含量，利用ArcGIS 10.2软件中的地统计分析功能，进行反距离加权插值分析，得到鄱阳湖丰水期表层沉积物样点分布图和地统计学空间插值图；沉积物重金属来源解析由USEPA PMF 5.0完成.

2 结果与讨论

2.1 表层沉积物重金属含量及空间分布特征

2.1.1 表层沉积物重金属含量的统计特征 从鄱阳湖丰水期88个表层沉积物样品重金属含量统计结果来看(表1)，Mo、Sb和W这3种重金属含量均呈偏态分布，经转换后基本符合正态分布，因此可用几何均值来表征. 除V和Cr外，Cd、Mo、Hg、Cu、Pb、Zn、W、Sr、As、Ni、Co和Sb的平均含量均超过了江西省土壤的背景值，分别为背景值的5.7、2.2、1.9、1.8、1.5、1.5、1.4、1.3、1.3、1.2、1.0和1.0倍；另外，Cd、Hg、Cu、Mo元素含量的最小值分别是其背景值的2.7、1.2、1.1、1.0倍. V和Cr的平均含量未超过江西省土壤背景值以及Co和Sb的平均含量与背景值相当的原因可能是与土壤和沉积物自身的性质有关，导致江西省土壤各元素中V、Cr、Co和Sb的背景值偏高. Cd、Hg、Cu、Mo、Pb、Sr和Zn超出背景值的比例相对较高，分别为100%、100%、100%、100%、97%、97%和93%. V、Cr、Co、Ni和Sb的EF值均小于1，而Cu、Zn、Sr、Mo、W、Pb、Hg和As的EF值在1～2之间，Cd的EF值最大，为4.93，说明V、Cr、Co、Ni、Cu、Zn、Sr、Mo、Sb、W、Pb、Hg和As为未污染至弱污染水平，而Cd属于中等污染水平，接近于重污染水平. 采用《土壤环境质量 农用地土壤污染风险管控标准》(GB 15618-2018)对鄱阳湖沉积物重金属污染进行评价，从平均值来看，Cr、Ni、Cu、Zn、Pb、Hg和As的含量均低于我国农用地土壤污染风险筛选值，而Cd的含量高于风险筛选值，但低于风险管制值，表明可能存在沉积物污染风险；所有采样点中Cd含量超过农用地土壤污染风险筛选值的百分比为51%. 总体而言，Cd的污染相对较严重.

表1 鄱阳湖丰水湖表层沉积物中重金属含量的描述性统计*

与先前(2011年)对鄱阳湖表层沉积物重金属的研究[30]相比，本研究Cr和Pb的含量较低，Cu、Zn、Cd、Hg和As的含量有所升高，说明鄱阳湖表层沉积物中Cu、Zn、Cd、Hg和As的污染状况有所加重. 将鄱阳湖丰水期表层沉积物重金属含量与第二、三大淡水湖(洞庭湖[30]、太湖[31-32])公布的数据进行比较，显然，本研究鄱阳湖表层沉积物中Cr、Cu、Zn、Cd、Pb、Hg和As的含量均低于洞庭湖(88.29、47.48、185.25、4.65、60.99、0.157 和29.71 mg/kg)；Cu、Cd和Pb的含量低于太湖(39.22、2.88和74.48 mg/kg)，而Cr、Zn和As的含量高于太湖(41.5、111.06和5.94 mg/kg). 这些比较的结果表明，洞庭湖表层沉积物中Cr、Cu、Zn、Cd、Hg和As均高于其他鄱阳湖和太湖，这可能与湖泊的地质背景差异有关.

2.1.2 表层沉积物重金属的空间分布特征 变异系数(CV)从侧面反映了沉积物中重金属元素的区域空间变异性，V、Cr、Co、Ni、Cu、Zn、Sr、Pb和Hg的变异系数较小(<27%)，属于弱变异，空间变异性相对较小；Mo、W、Sb、Cd和As的变异系数较大(>30%)，表现为中等变异，其中，Mo含量表现出最大的空间变异性(CV=46%). Cu、Zn、Sr、Pb、Hg和As的平均含量均超过背景值，但空间变异性相对较小，说明研究区域Cu、Zn、Sr、Pb、Hg和As的污染分布较均匀. 从图2可以清楚地看出重金属的空间分布特征十分明显，Cr、Cu、Zn、Sr、Pb、Hg和As的空间分布十分相似，主要特征表现为鄱阳湖区南部、西南部和东南部重金属含量高，低值区则位于湖区西北部，尤其是Zn和Hg空间分布呈现从南湖区到北湖区依次递减的规律. Cr、Cu、Pb和Hg含量的峰值均出现在饶河入湖口，Zn含量的最高值在赣江入湖口，而As、Sr含量的最高值分别出现在抚河、信江尾闾，研究表明，“五河”输入是入湖污染负荷的主要来源，其占污染负荷总量的80%左右. “五河”的沉积物是铜运输的主要工具，在“五河”流域的工业活动的影响下，特别是乐安河(饶河南支)中、下游的德兴铜矿和信江中游的永平铜矿矿山开采产生的含重金属的废水、废石堆的风化和矿石浮选厂的排放，它们占沉积物总负荷的22.1%，占鄱阳湖铜输入的35.6%以上[33]. 重金属含量的高值出现在饶河、赣江、信江和抚河入湖口也可能与信江流经的贵溪市有大型的有色金属冶炼厂、抚河上游的铀矿、赣南有色金属采矿区、赣江南支贯穿南昌市后携带大量工业废水以及水土流失严重造成土壤中的重金属流失有关，而修水流域林地覆盖面积大，农业面源污染及工业生产企业较少，外源污染输入污染较少[34-35]，表明高含量的重金属与河流入湖污染物有关，源于“五河”来水的污染是鄱阳湖沉积物中重金属含量水平及其区域性分布差异的主要原因.

图2 鄱阳湖丰水期表层沉积物重金属含量的空间分布

图2也揭示了Co、Ni、Mo和Sb含量高度相似的空间变化趋势，高值区基本上均位于鄱阳湖区南部、东北部和修水入湖附近的区域，低值区则位于湖区中部和北部. Co和Ni的平均含量与背景值接近，变异系数也小，空间分布较均匀，说明研究区域基本没有受到Co和Ni的人为污染. Mo的平均含量超出背景值较多，且空间变异性也最大，说明其含量集中分布，可能人为输入污染对其含量有较大影响. Sb的平均含量与江西省土壤背景值相当，但其变异系数较大，是因为Sb在修水入湖区域的含量(3.200 mg/kg)明显高于其他采样点，而修水流域内的矿产资源相对于饶河、信江、抚河及赣江较少，所以重金属污染相对于4条河流较轻[34]，而调查采样时发现可能与此区域的老爷庙水域处的地势特殊，水动力比较复杂，水体是多种水体混合而成有关[36]. V的空间差异不明显，平均含量也低于背景值，分布相对较均匀，其高值区域也分布在修水入湖附近的区域. Cd和W的空间变异性相对较大，主要原因是Cd含量在北部的入长江口附近(1.380 mg/kg)远高于其他地点，很明显，湖区北部Cd含量较高，且整个湖区黄色及以上区域超过农用地土壤污染风险筛选值(其中红色代表高含量，而蓝色代表低含量)；而W在赣江南支入湖口附近的含量(21.10 mg/kg)高于其他采样点，可能是由于钨矿如大余钨矿的开采冶炼在赣江上游[22]，正在重点推进钨资源精深加工产业集群建设的南昌国家高新技术产业开发区也坐落于赣江南支边上导致的.

表2 沉积物重金属元素主成分旋转载荷1)

2.2 表层沉积物重金属的来源解析

2.2.1 表层沉积物重金属的主成分分析 在进行PCA之前，应用Kaiser-Meyer-Olkin(KMO)和Bartlett球形度检验来确定数据对PCA的适用性，一般认为KMO>0.7时适合因子分析. 根据检验结果，KMO值为0.795，Bartlett球形度检验统计量的观测值为614.005，自由度为91，检验的显著性概率为0.000，小于显著性水平0.05，表明各重金属元素之间相关性较强，适合做因子分析. 为了更好地解释主成分，使用varimax归一算法对它们进行旋转，并选择旋转后的数据分析. 鄱阳湖丰水期表层沉积物重金属主成分分析结果见表2，通过varimax旋转以后，根据特征值大于1的原则，筛选出4个因子，解释了69%的变量信息. 表3列出了鄱阳湖丰水期表层沉积物14种重金属元素的Pearson相关系数.

由表2可知，PC1解释了26%的方差变异，Pb、Hg、Zn、Cu和As具有较高的正载荷，Cd和Cr具有较低载荷；且Pb、Hg、Zn、Cu、As和Cd两两之间存在极显著正相关(P<0.01)，Hg和As的相关系数最高，表明它们有相似的空间分布格局，均表现为高含量区域主要发生在湖区南部、西南部和东南部，Pb、Hg、Zn、Cu和As的空间分布亦是如此. 陈翠华等[37]在湖区东南部德兴地区采集水系沉积物样品分析中发现，Hg和As含量高的地区主要集中在德兴铜钼和银山铅锌矿区. 德兴铜矿区是亚洲最大的露天矿铜矿山，铜储量占中国铜储量的17.5%；德兴地区乐安河中下游还分布着朱砂红/富家坞铜钼矿和银山铅锌矿及其冶炼厂，经调查和统计资料显示[19]，德兴矿区已探明铅锌矿有数百万吨. 在铜矿、铅锌矿开采和冶炼过程中，某些硫化物(如闪锌矿、方铅矿、黄铜矿等)会被氧化释放出Pb、Hg、Zn、Cu、As和Cd进入环境[38]. 据不完全统计，每天从德兴铜矿排出超过4×104t矿山酸性废水[39]，废水中含有大量的Cu、Zn、As、Pb和Cd，含量分别为102.1、2.38、0.402、0.424和0.130 mg/L，大大超过《地表水环境质量标准》(GB 3838-2002)V类水中Cu、Zn、As、Pb和Cd的标准限值(分别为1.0、2.0、0.1、0.1 和0.01 mg/L)，废水流经鄱阳湖的主要支流-乐安河，最终都将排入鄱阳湖，表明鄱阳湖丰水期沉积物中Pb、Hg、Zn、Cu、As、Cd和Cr主要来源为铅锌矿、铜矿开采，选矿及其冶炼等生产活动中产生的废水、废渣等. 因此，第一主成分代表了矿业和工业活动等人为来源.

PC2解释了24%的总方差，Ni、Co、Sb和V在该成分的载荷较高(>0.7)，与它们相比，Mo占有较低的载荷，且Ni与Co、Sb之间存在极显著的强相关关系(P<0.01，r>0.7)，Mo与Ni、Co之间存在极显著的相关关系(P<0.01，r>0.5)，表明它们有相似的分布规律. 一般Ni源于金属冶炼废水的排放[40]，钢铁冶炼废水中含有大量的Ni、Co[41]和V[42]元素，矿产废渣中含有丰富的Sb元素[22]；其中，环境中的V也可能来源于成土母质[43]. 但本研究中Ni、Co和Sb的平均含量与背景值接近，V的平均值也未超过江西省土壤背景值，变异系数(12%)最小，总体上受人为活动较小，可认为主要受背景值控制，说明这些重金属主要来源于土壤母质造成的自然源，可认为第二主成分代表自然来源.

表3 不同重金属之间的相关系数

PC3相对贡献率为10%，Sr在该成分的载荷最高(0.85)，其次为Mo(0.53). Sr与Mo之间存在极显著的相关关系(P<0.01)，有研究表明，Sr来源于金属冶炼[44]，Mo与朱砂红/富家坞铜钼矿尾矿有关，重金属矿冶采选后废弃污染场地、尾矿废渣堆存场地众多[40]，其所在的德兴铜矿矿区建有亚洲第一大尾矿坝，已造成207 km2尾矿堆积区，因此，第三主成分代表了矿山开采所致的尾矿和废渣污染源.

PC4占总体方差的9%，W在该成分的载荷最高(0.92)，Cd和Hg也有较低载荷，暗示了Cd和Hg来源的多源性，且Cd与Hg之间存在极显著的相关关系(P<0.01). Dai等[45]认为Cd与鄱阳湖周围农田中肥料和农药的大量施用有关. Yang等[46]测定了肥料中的Cd和Hg含量分别为0.009～2.58和0.005～0.69 mg/kg. 而江西省的传统农业大县-南昌县和余干县均在鄱阳湖周边，农业是农民收入的主要来源，由于将大量农药和肥料等农业化学品投入土壤以获得更高的产量，这些农业用品的广泛和重复使用而导致Cd和Hg残留于土壤中，通过降雨和地表径流将其带入河流和湖泊，表明Cd和Hg来源于农业活动. 江西省辖区内朱溪钨铜矿三氧化钨资源量为2.86×106t，是世界上最大的钨铜矿[21]，此外，钨矿集中分布在赣江上游区，赣州市大余钨矿是我国最重要的钨矿山之一. 因此，可认为第四主成分主要受农业活动的影响.

2.2.2 基于PMF的重金属来源及相对贡献率 将88个表层沉积物样品的14种金属的含量数据和与这些含量相关的不确定度数据输入EPA PMF 5.0软件进行分析. 14种重金属均被归类为“Strong”(S/N>1)，总变量自动变为“Weak”. 设定不同因子数(3～7)，采用Robust模式，将运行次数设置为20次，选择“随机起始种子数”选项，最终确定4个因子能合理解释污染源的类别. 所有沉积物样本的绝大部分残差在-3～3之间. 总变量预测值与观测值之间相关性良好(r2=0.931)，表明表层沉积物重金属被PMF很好地分配. PMF解析出的表层沉积物重金属的4因子(来源)成分谱图如图3所示，4个因子的相对贡献率分别为38%、19%、28%、14%.

图3 沉积物重金属的源成分谱图

因子1具有高负荷的As(60%)、Sr(57%)、Cd(46%)、Hg(45%)、Cu(44%)、Pb(40%)、W(40%)和Zn(38%)，因子载荷与PC1相似，如上所述，金属矿石的采矿和冶炼的工业活动被认为是较重要的金属污染源之一，乐安河沿河流域的德兴铜矿、银山铅锌矿及其冶炼厂等，特别是上游和中游段，Cu、Zn、Pb年排放量就分别达141、69.9、0.5 t[47]；伍恒赟等[19]根据江西省环境监测中心站重点行业企业污染源调查可知，2006-2012年乐安河流域企业每年排放重金属工业废水均超过300万t，含6种主要重金属(As、Cd、Hg、Cu、Pb、Zn)，这些重金属工业废水排入乐安河后流入鄱阳湖，最终使鄱阳湖沉积物成为乐安河重金属的汇，因此，因子1是矿业和工业活动等人为来源.

因子2主要是Cr(50%)和W(44%). 有报道称化肥中Cr含量比较高，一般氮肥、磷肥及过磷酸钙中Cr含量分别为5～3000、30～3000和50～250 mg/kg；据统计，世界上通过肥料输入到土壤中的Cr为0.03～0.38 t/a[48]，过量施用磷肥会造成Cd和Cr的积累. 鄱阳湖地区丰水期正值水稻播种插秧季节，且赣江、饶河流域上游流域为江西省粮食主产区，氮肥和化肥使用量高[34]，而鄱阳湖流域特殊的地貌极易造成水土的流失. 据统计[49]，鄱阳湖流域平均土壤侵蚀总量达1.66×108t，相当于流失土壤5×104hm2，土壤中肥料极易随地表径流进入鄱阳湖，导致面源污染负荷增加，从而使Cr在赣江、抚河、信江和饶河入湖口聚集. W与PC4中的W同来源于钨矿，因此，因子2与PC4结果一致，均受农业活动的影响.

因子3由Zn(48%)大量加权，Hg(36%)、As(34%)和Pb(31%)略微加权. Zn与铅锌矿尾矿有关，且矿产废渣中含有丰富的As元素[50]，经调查和统计资料显示[19]，德兴银山矿区已探明铅锌矿矿产储量分别高达43万t Pb和52万t Zn[51]，乐安河下游还分布银山铅锌冶炼厂，因此，因子3为矿山开采所致的尾矿和废渣污染源. 因子4主要负载在Mo(37%)、Ni(33%)、Sb(32%)、Co(31%)和V(28%)，得到的因子载荷与PCA结果中的PC2一致，代表了自然来源.

2.3 结合PCA和PMF的重金属来源及相对贡献率

PCA和PMF 2种源解析方法的来源分配和相对贡献率见表4，源解析结果中污染来源一致，相互印证. 两种方法均解析出贡献率相当高的矿业和工业活动来源，贡献率分别为26%和38%；尾矿和废渣、农业活动、自然来源是2种方法解析出的其他3个源，PCA和PMF解析得到的相对贡献率分别为10%、9%、24%和28%、19%、14%. 两种方法对重金属污染物来源贡献率的预测略有差别，这主要是由于在方法应用过程中所选择的不确定性参数和选为变量的污染物之间存在差异所致[52]. 通过计算PMF预测的和实验的重金属含量之间的回归进行内部比较，可用于确定模型是否适用[53-54]，结果表明，重金属含量总变量预测值与观测值之间相关性良好(R2=0.931)，说明沉积物重金属被PMF很好地分配，PMF可以很好地应用于沉积物重金属污染的源解析研究并分配污染源对每种重金属的贡献. 先前的研究[23,55]中，使用PMF也发现了污染物的预测浓度和实验浓度之间这种良好的相关性. 与其他学者采用主成分分析法研究鄱阳湖重金属污染源的源解析结果[45]相比，主成分分析法源解析结果中源识别数量的减少和其他信息的丢失可能是由于没有考虑数据的不确定性或没有非负性约束，由表2和图3可知，主成分载荷值为负值，且PCA解析结果虽与源贡献相关，但是并不呈比例，使得累积贡献率少于100%；而PMF源成分谱图中的值均为正值，所以PMF方法可以同时将测量不确定性和非负性约束纳入其计算过程，依赖于更具物理意义的假设，可以提供更好的结果.

本研究综合PCA和PMF的源解析结果得出，鄱阳湖丰水期沉积物重金属有4个主要来源，分别为矿业和工业活动、尾矿和废渣、农业活动、自然来源，源相对贡献率分别为38%、28%、19%、14%. 流域内的矿业和工业活动对鄱阳湖表层沉积物重金属累积的贡献比较明显，表明重金属污染主要受矿业和工业活动等外界人类活动的影响. 矿业和工业活动产生的重金属经南湖区4条主要支流-赣江、抚河、信江和饶河携带后释放到径流水体并迁移汇入鄱阳湖，最终悬浮物中的重金属也都转入沉积物；且Cr、Cu、Zn、Sr、Pb、Hg和As 7种重金属元素在这4条主要支流入湖口处都显示出高度相似的趋势，即高含量聚集在河流入湖口附近，这可能是因为丰水期湖泊水量较大，水位较高，4河来水进入湖泊后，水体流速骤降[45]，因此携带的重金属易于沉积，导致各重金属含量高且分布集中. 为防止沉积物重金属含量继续增加而造成更大的生态风险，鉴于鄱阳湖5大河流矿区工业废水不达标排放，应加强对企业排污的监控和治理力度以及地表径流(“五河”)的重金属入湖通量控制，有针对性地进行环境污染控制.

表4 PCA和PMF的来源和相对贡献率

3 结论

1)除V和Cr外，Co、Ni、Cu、Zn、Sr、Mo、Cd、Sb、W、Pb、Hg和As的平均含量均超过了江西省土壤的背景值；Cd、Hg、Cu、Mo、Pb、Sr和Zn超出背景值的比例相对较高，且所有采样点中Cd含量超过农用地土壤污染风险筛选值的百分比为51%；V、Cr、Co、Ni、Cu、Zn、Sr、Mo、Sb、W、Pb、Hg和As为未污染至弱污染水平，而Cd属于中等污染水平，接近于重污染水平. 总体而言，Cd的污染相对较严重.

2)重金属的含量分布具有显著的区域特征，其中Cr、Cu、Zn、Sr、Pb、Hg和As的空间分布十分相似，表现为在赣江、抚河、信江和饶河入湖口附近区域含量较高，而Co、Ni、Mo和Sb明显在湖区南部、东北部和修水入湖附近聚集，Cd和W的空间变异性相对较大，V的含量分布相对较均匀.

3)PCA和PMF解析结果都表明鄱阳湖丰水期表层沉积物重金属受矿业和工业活动、尾矿和废渣、农业活动、自然来源4种来源的共同影响，源相对贡献率分别为38%、28%、19%、14%；PMF的重金属含量总变量预测值与观测值之间相关性良好，说明PMF可以很好地应用于沉积物重金属污染的源解析研究. 两种方法的源解析结果相互印证，具有很好的一致性.