杞麓湖径流区不同湿地沉积物重金属污染特征及潜在生态风险评价

浦江，张翠萍，刘淑娟，3，杨小燕，赵斌，李淑英，陆轶峰，王媛媛，周元清*

（1.云南大学生态与环境学院，昆明 650091；2.玉溪师范学院污染控制与生态修复研究中心，云南 玉溪 653100；3.云南农业大学动物科学技术学院，昆明 650201；4.云南农业大学资源与环境学院，昆明 650201；5.玉溪市抚仙湖保护开发投资有限责任公司，云南 玉溪 653100）

湖泊富营养化和重金属污染一直是生态学和环境科学研究的重要课题[1]。2016 年以来，随着湖泊污染控制与生态修复工程推进，农业面源污染治理效果显著，湖泊生态环境明显改善，低污染水治理已成为我国湖泊水环境保护的重要组成部分。山水林田湖草生命共同体的理念将湖泊生态环境扩展为集水区-湿地-汇水河流等有机结合的开放水系统。享有“地球之肾”美誉的湿地，具有重要生态系统服务价值和多种生态功能[2]。湿地是各类污染物重要的“汇”[3-4]，也是研究污染物迁移转化的敏感过渡带，作为流域低污染水入湖生态屏障，对拦截、过滤、削减入湖污染物和改善水环境质量至关重要。农业径流区农药化肥大量施用所导致的重金属污染不容忽视。低污染水中重金属经絮凝沉降、氧化还原等过程积累于沉积物中，当环境因素变化时，部分重金属会释放进入水体，并通过食物链富集，危害人类健康[3，5-6]。

重金属作为低污染水中典型污染物，其空间分布和生态风险是湿地生态环境中的关键问题[7]。湿地受人类活动和土地利用方式等影响[8]，沉积物中重金属具有一定空间分异特征，摸清湖滨带低污染水不同类型湿地沉积物重金属空间分布特征和生态风险，对湖泊水环境保护具有重要指导意义[9]。国内外关于湖滨带低污染水的研究，主要集中在水处理技术和氮磷去除[10-11]，JIANG 等[12]基于PMF-RI 和PMF-HHR 方法对湿地土壤进行环境风险评价，但流域低污染水对湖滨带湿地重金属污染的贡献不容忽视[13]。

杞麓湖为云南省九大高原湖泊之一，通海县是云南最大蔬菜种植和集散地[14]，流域内蔬菜花卉复种指数较高，农业面源污染负荷较重[15]，杞麓湖农田径流区是研究低污染水所致重金属污染的典型区。近年来，随着杞麓湖“以渔净水”项目和国家湿地公园恢复工程的实施，湿地成为“山-城-林-田-湖”空间格局的重要节点，湿地面临的污染形势更严峻，保护湿地生态环境尤为迫切。因此，本研究以农业径流区作为探究低污染水对不同类型湿地重金属污染的切入点，摸清杞麓湖农业径流区湿地沉积物重金属空间分布和生态风险，旨在为杞麓湖流域环保预警和低污染水治理提供基础数据。

1 材料与方法

1.1 研究区概况

杞麓湖（102°30′25″～102°52′53″E，23°55′11″～24°14′49″N）属珠江流域南盘江水系，滇中高原封闭型断陷湖泊，径流面积354.94 km2，总人口24万，流域四周高、中部低，湖周为平坝区[15]。湖泊长轴呈东西向新月形，长10.4 km，南北较窄，平均宽3.5 km，湖岸线全长32 km，最大水深6.84 m，平均4 m。湖泊海拔高程多在1 979～2 100 m 之间，属亚热带湿润高原季风气候，流域年均气温15.5 ℃，水源主要来自降雨径流补给，年降雨量800～1 100 mm，蒸发量小，气候湿润[16]。杞麓湖属典型富营养老龄化湖泊，是关系通海县社会经济发展主体的“母亲湖”，湖盆内分布着通海县90%人口，流域内包含4 条主要入湖河流，即红旗河、者湾河、中河和大新河，覆盖了全县7 个乡镇。水质长期处于劣Ⅴ类，水质污染以有机污染和氮、磷为主[17]。

1.2 采样点布设与指标分析

结合杞麓湖流域现状及地形地貌，于2017—2019年8 月（丰水期）在杞麓湖农业径流区设32 个采样点（其中河流湿地n=13、人工湿地n=10、湖泊湿地n=9，图1）。在不同类型湿地样地内，用抓斗式采泥器采集0～30 cm 表层沉积物，用手持式GPS（GARMIN Oregon 739，上海佳明）记录采样点经纬度。每个湿地按梅花形采集5 个以上沉积物样品，混匀为2 kg，去除杂质后装入封口袋，带回实验室制样，待测。土壤重金属Cu、Zn、Pb、Cd、Cr、Ni 用王水回流消解原子吸收法（NY/T 1613—2008）测定；Hg 用原子荧光法（GB/T 220105.1—2008）测定；全氮（TN）用半微量凯氏蒸馏定氮法（NY/T 53—1987）测定；全磷（TP）用钼锑抗比色法（NY/T 88—1988）测定；pH 用PHS-3C 型酸度计测定；有机质（OM）用硫酸与重铬酸钾氧化-外加热-容量法（NY/T 1121.6—2006）测定。实验室加入空白样品和标准物质进行质量控制，各测试项目相对标准偏差（RSD）小于10%，每个样品三组平行实验，平行样间的误差小于10%，按平均值统计。

图1 研究区位置及采样点分布Figure 1 Location of the study area and the sampling point distribution

1.3 评价方法

1.3.1 单因子污染评价

以云南省土壤重金属背景值[18]（表1）作为参比值，用单因子污染指数法评价湿地沉积物中重金属对环境的污染程度[19]。

1.3.2 潜在生态风险评价

采用Hakanson 潜在生态风险指数法[18]评价杞麓湖农业径流区湿地重金属污染，重金属毒性响应系数（）见表1，表示重金属毒性水平和生物对重金属污染的敏感程度，潜在生态风险指数（IR）算式如下：

式中：为重金属i相对于背景值的污染指数；为重金属i的实测值，mg·kg-1；为重金属i的背景参考值，mg·kg-1；为重金属元素i的单因子潜在生态风险指数，和IR值越大，生态危害程度越大，其分类等级见表2。样点污染程度比例=重金属元素不同污染等级样品数/样品总数×100%，样点潜在生态风险比例=潜在生态风险评价等级样品数/样品总数×100%。

表2 污染程度和潜在生态风险程度分类等级Table 2 Degree of pollution and potential ecological risk

1.4 数据统计分析

用Excel 2013进行数据统计，用SPSS 20.0软件进行Pearson 相关性分析和主成分分析，不同处理结果用平均值±标准差（Mean±SD）表示，用Origin 8.0 软件作图。用ArcGIS 15.0 软件绘制采样点分布图，用ArcGIS反距离权重插值法（IDW）分析研究区数据，并绘制重金属空间分布图及潜在生态风险分布图。

2 结果与分析

2.1 农业径流区湿地重金属的污染特征

由表3 可知，杞麓湖农业径流区不同类型湿地沉积物重金属平均含量总体表现为河流湿地＞人工湿地＞湖泊湿地；同类湿地重金属平均含量排序为Zn＞Cr＞Cu＞Ni＞Pb＞Cd＞Hg；河流湿地沉积物Zn 和Ni 标准差较大，表明采样点变异较大，主要是由于湖泊西北部纳家营轧钢厂附近采样点数值偏高。除Pb 外，河流湿地其余6 种重金属平均含量均超过云南省土壤背景值，Cu、Zn、Cd、Cr、Ni、Hg 分别为土壤背景值的1.84、2.10、4.49、1.51、2.23、5.17 倍，人工湿地和湖泊湿地Zn、Cd、Cr和Hg 含量均超过背景值，表明杞麓湖农业径流区低污染水中存在不同程度重金属污染，重金属Cd 污染严重，在河流湿地和人工湿地超背景值的样点比例高达100%。不同类型湿地的变异系数差异明显，河流湿地Ni 和Hg、湖泊湿地Pb 的标准差高于平均值，变异系数较大，表明3 种元素含量空间分布不均匀，离散性相对较大。

表3 研究区域土壤重金属污染特征分析Table 3 Descriptive statistics of heavy metals in soils of the study area

2.2 农业径流区湿地重金属的空间分布特征

杞麓湖农业径流区湿地7 种重金属空间分异特征明显，总体呈西北部高于东南部，自北向南、自西向东递减趋势（图2）。其中，Cu、Zn、Cr 和Cd 空间含量变化大体一致，具有一定相似性，空间分异较大。Cd存在多个点源，高值区较多；Ni和Hg除高值区外整体分布均匀，无明显空间分异，Hg高值区出现在西南部杞麓湖国家湿地公园，Ni 空间分布最均匀，且平均含量略高于云南省土壤背景值；Pb 部分区域存在一定空间分异，除部分区域含量偏高外，整体分布均匀，西部和东部部分样点含量稍高。重金属污染高值主要在湖泊西北部纳家营轧钢厂附近及西南部杞麓湖国家湿地公园。

图2 杞麓湖湿地沉积物重金属的空间分布Figure 2 Spatial distribution of heavy metals in sediment of Qilu Lake wetlands

2.3 农业径流区湿地重金属污染评价

2.3.1 单因子污染指数法

以云南省土壤背景值为参比值，杞麓湖农业径流区湿地沉积物均受到不同程度的重金属污染。除Pb和Ni 外，三种不同类型湿地沉积物重金属Cu、Zn、Cd、Cr 和Hg 的单因子污染指数均大于1，呈现不同程度的富集。其中，河流湿地Cd 和Hg 的污染指数高达8.39和24.02，重度污染占比76.93%和38.46%（图3）。说明Cd 和Hg 富集程度高，呈重度污染，存在较高生态风险。

图3 基于云南省土壤背景值的单因子污染指数Figure 3 Single factor pollution index based on soil background value in Yunnan Province

2.3.2 潜在生态风险指数法

通过反距离权重插值法结合潜在生态风险指数对各风险等级占比进行统计，结果如图4a 所示，农业径流区各采样点潜在生态风险指数平均值为264，属中度生态风险水平。人工湿地和湖泊湿地强生态风险样点约占11%，大部分样点处于轻微和中度生态风险水平，河流湿地强生态风险样点占比高达30.7%，局部出现生态风险高值区，极强生态风险样点占15.39%，潜在生态风险指数最大值达1 241。潜在生态风险空间分布呈西高东低的趋势（图4b），西部地区处于强生态风险水平，东部地区处于中度风险水平，少数采样点周围处于轻微风险水平。其中，极强潜在生态风险区主要分布在西北部杞麓湖纳家营轧钢厂附近及西南部的杞麓湖国家湿地公园。

图4 研究区潜在生态风险指数（IR）等级占比（a）和空间分布图（b）Figure 4 The map of risk level proportion（a）and spatial distribution（b）of potential ecological risk index（IR）in the study area

2.4 农业径流区湿地沉积物重金属空间分异特征及污染源分析

杞麓湖农业径流区湿地表层沉积物各重金属与土壤理化因子相关性分析见表4，结果表明：Cu与Zn、Cr、Ni，Zn与Pb、Ni、Hg，Ni与Cr之间呈极显著相关（P＜0.01），推测Cu、Zn、Cr、Ni、Hg 5种重金属污染具有同源性，7种重金属元素与土壤pH、OM、TN、TP间存在一定相关性，其中TN 与Cd 极显著正相关（P＜0.01），TP 与Pb、Cr 显著负相关（P＜0.05），Zn、Pb、Cr、Ni、Hg 与OM呈极显著正相关（P＜0.01）。

表4 土壤理化因子与重金属含量的相关性Table 4 Pearson correlation of soil physical and chemical factors and heavy metal concentrations

运用SPSS 软件主要提取到PC1 和PC2 两个主成分（表5），累计方差贡献分别为50.77%和68.99%，能较好解析污染源信息。第一主成分（PC1）中Zn、Cr、Ni 载荷较高，分别为0.913、0.727 和0.834，由相关性分析结果可知，这3 种元素具有较好相关性且空间分布特征相似，说明来源可能一致。

表5 土壤重金属元素因子载荷Table 5 Factor matrix of heavy metal elements of soil

第二主成分（PC2）中Cd 和Pb 载荷值分别为0.623和0.390，其中，Cd与其他重金属相关性不显著，说明Cd具有单独来源，且与TN呈极显著正相关。

3 讨论

3.1 杞麓湖径流区不同湿地沉积物重金属分布特征

杞麓湖农业径流区不同类型湿地沉积物重金属平均含量表现为河流湿地＞人工湿地＞湖泊湿地，原因是湿地对重金属具有削减作用，由农业径流区上段城镇生活低污染水、交通运输以及农业生产中化肥农药施用等产生的重金属经地表径流或大气沉降进入河流湿地，一部分重金属在河流湿地沉积物中积累，另一部分重金属先进入人工湿地，经湿地削减后进入湖泊湿地，最后进入湖泊水体。除Pb外，河流湿地其余6 种重金属平均含量均超过云南省土壤背景值，河流湿地和人工湿地Cd 超标严重，超标率达100%。变异系数可以反映各样本数据变量空间上的离散程度，有研究表明，变异系数较小的元素自然来源占主导地位，而变异系数较大的元素则受人为影响较大[20]。Ni、Hg和Pb这3种重金属呈强变异，表明河流湿地Ni和Hg、湖泊湿地Pb 受人为活动干扰严重，主要受径流区农业活动影响，还与流域内钢铁产业及其加工制造业等有关。

杞麓湖农业径流区湿地沉积物Cu、Zn、Cr 和Cd具有相似的空间分布特征，整体呈西北高东南低，自北向南、自西向东递减的趋势。Cd 空间分异较大，存在多个点源，整体为高值分布，Hg高值区出现在西南部杞麓湖国家湿地公园，该区域人类活动密集，环湖公路车辆较多，其产生的尾气中重金属Hg 经大气沉降随雨水进入湿地，导致该区Hg浓度偏高，这主要与通海县坝区密集人类活动和交通运输有关。Ni 空间分布最均匀，且平均含量略高于云南省土壤背景值，其高值贡献区主要来自湖泊西北部纳家营轧钢工业，说明Ni元素主要受成土母质和点源影响。Pb存在一定空间分异，除部分区域含量偏高外，整体分布均匀，西部和东部局部地区含量稍高。

3.2 农业径流区沉积物重金属潜在生态危害解析

湿地沉积物重金属的风险评价一直是湿地生态环境评价的核心方法[21]，基于云南省土壤背景值进行评价，从单因子污染指数看，三种类型湿地Cd 和Hg呈重度污染，富集严重，存在较大生态安全隐患，应优先控制。由于湿地沉积物对Cd 和Hg 具有较强的吸附作用[22]。强潜在生态风险主要贡献区是西北部杞麓湖纳家营轧钢厂及西南部的杞麓湖国家湿地公园，高值区内有红旗河、者湾河、中河和大新河4 条主要入湖河流（图1），土地利用类型主要以农业用地为主，人类活动频繁，地表植被覆盖率低，以农作物为主，因此整个西北部地区潜在生态风险偏高。其中，杞麓湖西南部和西北部坡面是流域内重要水源涵养区，也是生态敏感性和土壤可蚀性较高地区[23]，应得到管理部门重视。

3.3 农业径流区沉积物营养盐与重金属的相关性

依据重金属相关性可间接推测沉积物重金属污染来源和途径，若各元素间相关性显著或极显著，表明在该地区重金属元素具有同源性或呈复合污染[24]。其次，湿地沉积物理化性质也直接影响重金属随低污染水进入湿地的环境行为[25]。本研究中，TN 与Cd 呈极显著正相关（P＜0.01），TP 与Pb、Cr 显著负相关（P＜0.05），Zn、Pb、Cr、Ni、Hg 与OM 呈极显著正相关（P＜0.01），农业径流区低污染水OM 含量丰富，对重金属络合能力强，这可能是造成该径流区表层沉积物中Zn、Pb、Cr、Ni、Hg 含量较高的主要原因。农业径流区湿地沉积物pH 平均值为7.73，呈弱碱性，7 种重金属均与pH呈负相关，其中Cd与pH呈极显著负相关（P＜0.01），这是因为pH 越高，土壤表面的负电荷数量越多，与有机酸形成配合物的机率也越大，导致重金属吸附能力减弱，重金属在土壤中的移动性增强[26]。

主成分分析法（PCA）是研究多个变量相关性的多元统计方法，能较好地根据相关性密切程度进行元素分类，通过各因子中的特征元素判断元素来源[27-28]。本研究第一主成分为混合来源，主要与人类生活和农业生产有关。第二主成分为Cd 和Pb，载荷值分别为0.623 和0.390，其中，Cd 与其他重金属的相关性不显著，说明Cd 具有单独来源，且与TN 呈极显著正相关，表明Cd与农业活动和农药、化肥施用等有关[29]，此外，Cd 污染与工业生产有关，可解释杞麓湖西北部纳家营钢铁产业样点Cd污染严重。结合风险评价和空间分布结果，Cd 应作为该区域污染防控与生态修复重点。Pb 主要在车辆的燃料燃烧和催化剂、铅蓄电池的使用中大量产生[30]，Pb 的高值区主要在杞麓湖国家湿地公园环湖公路附近，说明其主要来自人类活动与交通运输。

4 结论

（1）杞麓湖农业径流区不同湿地类型的重金属含量存在一定差异，表现为河流湿地＞人工湿地＞湖泊湿地，同类湿地中重金属平均含量大小为Zn＞Cr＞Cu＞Ni＞Pb＞Cd＞Hg，Pb、Ni 和Hg 呈强变异性，其中Cd 和Hg超背景值样点比例较高，应重点防控。

（2）ArcGIS 反距离权重插值法分析结果表明，Cd、Cu、Cr、Zn含量具有相似空间分异特征，呈自北向南、自西向东递减趋势，两个含量高值区分别是西北部的纳家营钢铁厂附近湿地和西南部的杞麓湖国家湿地公园，应予优先关注和重点管控。

（3）单因子污染指数结果表明Cd 和Hg 呈重度污染，存在危害风险，其余5 种重金属污染程度较轻。农业径流区潜在生态风险指数介于112～1 241，平均264，属中度风险水平，且生态风险空间分布呈西高东低趋势，与含量空间分布特征一致。

（4）相关性分析和主成分解析结果表明，径流区湿地沉积物各重金属主要受人类农业活动干扰。Cu、Ni、Cr 和Zn 的污染源主要是农业生产活动；Pb 来源主要与道路交通有关，还受工业活动影响；Cd 主要受农业活动和工业活动影响。

研究区采样点的位置和数量对分析结果具有一定程度的影响，本研究缺乏沉积物分层采样和重金属形态分析，还需在后续研究中提高采样密度并进行重金属形态的深入研究。