柳江流域沉积物重金属生态风险评价及来源分析
Ecological risk assessment and source analysis of heavy metals in sediments of Liujiang River Catchment
通讯作者: 李杰(1983-),男,博士,高级工程师,主要从事生态地球化学和勘查地球化学应用研究工作。Email:187292565@qq.com
责任编辑: 蒋实
收稿日期: 2019-04-11 修回日期: 2019-07-18 网络出版日期: 2020-02-20
基金资助: |
|
Received: 2019-04-11 Revised: 2019-07-18 Online: 2020-02-20
作者简介 About authors
钟晓宇(1983-),男,广西南宁人,硕士,工程师,主要从事地球化学研究工作。Email:walkofmelody@163.com 。
以柳江中下游流域沉积物为研究对象,采集91件水系沉积物样品,计算并分析As、Cd、Cr、Cu、Hg、Ni、Pb、Zn等元素地球化学参数特征,对当前污染现状及潜在风险进行预测评价,并解析重金属元素来源。结果显示:柳江流域重金属元素含量分布表现出较强的地域特征;柳江流域污染等级以轻度及偏中度为主,潜在生态风险等级为中低等级,主要影响因子为Cd;As、Cr、Cu、Hg、Ni来源以自然源为主,Cd、Pb以人为源为主,其中Cd自然源贡献率高值区主要分布于龙江流域,工矿业源贡献率高值区主要分布于金城江周边以及柳江柳州市区段。
关键词:
In this study, 91 sediment samples were collected from the middle and lower reaches of Liujiang River. The geochemical characteristics of As, Cd, Cr, Cu, Hg, Ni, Pb, Zn and other elements were analyzed. The current pollution situation and potential risks were predicted and evaluated, and the sources of heavy metals were analyzed. The results show that: The spatial dispersion of elements in the sediments of Liujiang River Catchment is relatively large, and the distribution of heavy metal elements shows strong regional characteristics.Distribution of heavy metal elements in Liujiang River Catchment shows strong regional characteristics.As, Cr, Cu, Hg and Ni are mainly from natural sources, while Cd and Pb are mainly from human sources. The high-value areas of Cd natural source contribution mainly distribute in Longjiang River Catchment, while the high-value areas of industrial and mining sources contribution mainly distribute around Jincheng River and Liuzhou District of Liujiang River.
Keywords:
本文引用格式
钟晓宇, 吴天生, 李杰, 郑国东, 卓小雄, 关东超, 王磊, 莫斌吉.
ZHONG Xiao-Yu, WU Tian-Sheng, LI Jie, ZHEN Guo-Dong, ZHUO Xiao-Xiong, GUAN Dong-Chao, WANG Lei, MO Bing-Ji.
0 引言
柳江属中国珠江水系西江干流第二大支流,全长773.3 km,流域面积5.7万km2,水力资源丰富,不仅为周边区域人类生活、农业灌溉及工业用水来源,亦是黔桂水上交通要道。该区域人口稠密,工农产业发达,坐落有区内重要的综合性工业基地,亦是广西主要经济及粮食作物产地[5]。为查明柳江流域目前污染状况,笔者通过对柳江流域沉积物元素含量地球化学特征分析研究,评价污染现状,预测潜在风险,并开展重金属来源解析,为今后合理开发利用水资源及制定污染防治措施提供依据。
1 研究区概况
研究区地处柳江中下游河段,(部分)涵盖柳州市、来宾市等4市6县。柳江干流中段称为融江,流经都安、柳城等县,河道相对较窄,滩多流急;下游段俗称柳江,流经柳州市、象州县,水势平缓,河曲较发育,水系多呈树枝状,大部分属喀斯特地貌,沿河阶地、丘陵广布。下文分别用融江及柳江代指柳江水系干流中、下游河段。支流则包括金城江、龙江、洛清江、罗秀河。整个柳江水系地表径流大体顺地势自西北向东南,横贯整个区域,大小支流主要从两侧汇入。
2 样品采集及分析测定
2018年6~10月,使用抓斗式取样器均匀采集柳江流域水底0~20 cm沉积物共91件(图1), 采集样品原始质量约3 000 g,样品在室温下自然风干并剔除其中的水草、贝壳、生物碎屑等杂质,物理敲打破碎过20目尼龙筛后装入自封袋送至实验室。分析测试方法如下:使用X射线荧光光谱法(XRF)测定Cr、Pb、Zn、P、S、Al2O3、CaO、Fe2O3,使用原子荧光法(AFS)测定As及Hg,使用全谱直读等离子体光谱法(ICP-OES)测定Cu、Mn、Ni,Cd、Corg则分别用石墨炉原子吸收光谱法(GFAAS)以及容量法(VOL)测定。结果显示,以上分析方法的检出限、检出合格率均满足或优于《多目标区域地球化学调查规范(1∶250 000)》(DZ/T 0258-2014)的要求。
图1
图1
柳江流域沉积物采样点位置
Fig.1
The location of Liujiang Rriver and the distribution of sampling points
3 研究方法
3.1 地累积指数法
式中:Cn为沉积物中元素含量;Bn为背景浓度,采用研究区表层土壤地球化学背景值作为参比值;k为修正系数,一般为1.5。土壤中重金属含量分级见表1。综合污染程度Igeo为取各元素污染最大值。计算式为:Igeo=Max(Igeo)。
表1 地累计指数与污染程度分级
Table 1
污染程度 | 无 | 轻度 | 偏中度 | 中度 | 偏重度 | 重度 | 严重 |
---|---|---|---|---|---|---|---|
Igeo | ≤0 | 0~1 | 1~2 | 2~3 | 3~4 | 4~5 | >5 |
级别 | 0 | 1 | 2 | 3 | 4 | 5 | 6 |
3.2 潜在污染预测
单因子潜在生态风险系数(the potential ecological risk index)计算公式为:
式中:
表2 广西中部地区土壤基准值
Table 2
元素 | As | Cd | Cr | Cu | Hg | Ni | Pb | Zn |
---|---|---|---|---|---|---|---|---|
背景值/10-6 | 24.17 | 459 | 106 | 30.9 | 203.8 | 45.65 | 33 | 117 |
多因子潜在生态危害指数RI(the potential ecological risk index)计算公式为
重金属潜在风险系数及潜在生态危害指数分级标准见表3。
表3 重金属潜在生态风险指数分级标准
Table 3
单因子潜在风险系数 | 单因子风险等级 | 潜在生态危害指数RI | 多因子风险等级 |
---|---|---|---|
<40 | 低 | <150 | 低 |
40~80 | 中 | 150~300 | 中 |
80~160 | 较高 | 300~600 | 较高 |
160~320 ≥320 | 高 极高 | >600 >600 | 极高 极高 |
3.3 绝对主成分分析—多元线性回归模型(APCA-MLR)
该理论首先由Thurston和Spengler在1985年提出[11],能够定量确定各污染源的污染贡献率。该方法最开始主要用于大气污染物的来源解析[12,13],对水体、沉积物中重金属(有机质)来源解析也应用较广[14,15,16],它不但可以定量确定每个变量对每个源的载荷,还可以定量确定源对其重金属的平均贡献量和在每个采样点的贡献量。计算过程为:对所有重金属元素含量进行标准化,经主成分分析后求出归一后的各样本的因子分数,引入0浓度样本,计算出0浓度样本的因子分数,归一后的各样本因子分数减去0 浓度样本的因子分数得到每个重金属元素的绝对主成分因子,用重金属浓度数据对绝对主成分因子做多元线性回归,得到的回归系数可将绝对主成分因子转化为每个污染源对每个样本的浓度贡献:
式中:b0i为对金属元素i做多元线性回归所得的常数项;bpi是源p对重金属元素i的回归系数;APCSp为调整后的因子p的分数;APCSp·bpi 表示源p对Ci的质量浓度贡献;所有样本的APCSp·bpi平均值即为源平均绝对贡献量。
4 结果与分析
4.1 地球化学参数特征
表4 柳江流域沉积物元素参数统计
Table 4
参数 | As | Cd | Cr | Cu | Hg | Ni | Pb | Zn | Mn | P | S | Al2O3 | Fe2O3 | CaO | Corg |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
最小值 | 2.2 | 0.06 | 23.4 | 6.6 | 0.016 | 8.5 | 6.3 | 19.1 | 75.4 | 154 | 92 | 3.55 | 1.45 | 0.07 | 0.15 |
平均值 | 17.6 | 1.04 | 53.7 | 21.6 | 0.12 | 22.5 | 32.7 | 112 | 688 | 422 | 262 | 9.79 | 3.77 | 0.58 | 0.68 |
最大值 | 45.8 | 6.86 | 89.1 | 52 | 0.474 | 38.9 | 169 | 503 | 1275 | 2173 | 2129 | 17.1 | 6.59 | 3.94 | 2.5 |
变异系数 | 0.49 | 0.91 | 0.26 | 0.39 | 0.61 | 0.34 | 0.8 | 0.69 | 0.42 | 0.5 | 0.83 | 0.34 | 0.27 | 0.97 | 0.63 |
注:Corg及氧化物含量单位为%,其他元素含量单位为10-6
柳江不同河段沉积物重金属元素含量相比,Cr、Cu、Ni含量相差不大,As、Cd、Hg、Pb、Pb则差异明显,且分布具有典型的地域特征。金城江Cd、Hg、Pb、Zn含量明显高于其他水系,均值分别为4.07×10-6、0.296×10-6、112×10-6、316×10-6;As高含量区域分布于融江流域,平均含量为35.5×10-6,为柳江流域平均值的1.9倍;重金属元素的低含量区则集中于罗秀河、洛清江流域,仅为柳江流域平均值的20%~98%。不同水系重金属元素含量详见图2。
图2
图2
不同河段重金属元素含量差异对比
Fig.2
Comparison of heavy metal elements contents in different reaches
4.2 污染现状
柳江流域环境质量良好,综合污染等级以轻度及偏中度为主,占比分别为31.87%、32.97%。主要影响因子为Cd,Igeo范围为-3.5~3.3,平均为 0.6;其余元素Igeo均小于0,其中Cr、Ni无污染。总体上,80%的水系沉积物样品中As、Cu、Hg、Pb、Zn无污染。不同河段污染现状相比,柳江支流水系污染现状呈现两极分化;洛清江及罗秀河底积物最为清洁,无污染等级样点分别占85.71%及90%;金城江沉积物污染程度最大,中度及以上污染等级样点占比63.64%。不同河段沉积物重金属元素地累积指数见表5。
表5 柳江各河段沉积物重金属元素地累积指数
Table 5
河段 | As | Cd | Cr | Cu | Hg | Ni | Pb | Zn |
---|---|---|---|---|---|---|---|---|
融江 | -0.06 | -0.02 | -1.45 | -0.83 | -1.32 | -1.40 | -0.81 | -0.73 |
柳江 | -0.77 | 1.05 | -1.51 | -0.93 | -1.14 | -1.42 | -0.62 | -0.45 |
金城江 | -0.51 | 2.33 | -1.48 | -1.24 | -0.16 | -1.65 | 1.04 | 0.68 |
龙江 | -1.20 | 0.69 | -1.48 | -1.25 | -1.60 | -1.83 | -0.78 | -0.81 |
洛清江 | -2.28 | -1.01 | -1.87 | -1.35 | -2.02 | -1.84 | -1.58 | -1.70 |
罗秀河 | -2.23 | -2.00 | -1.86 | -1.38 | -2.96 | -1.95 | -1.00 | -1.92 |
4.3 风险预测
柳江潜在生态风险等级为中低等级,其中单元素Cd潜在风险系数均值为101,整体属于较高风险等级,较高风险等级占比35.16%,中等风险等级占比27.47%。As、Cr、Cu、Hg、Ni、Pb、Zn潜在风险系数值均小于40,除Hg有少量中—较高风险等级区域分布外,其余元素皆为属低生态风险等级。各重金属的生态风险等级整体排序为Cd>Hg>As>Pb>Cu>Ni>Zn>Cr。
融江、柳江、金城江、龙江、罗秀河、洛清江综合风险指数均值分别为107、155、263、130、37.8、56.6。金城江流域属中等风险等级,明显高于其他水系,生态风险压力较大,主要受Cd、Hg因子的影响,采集的11个水系沉积物样点中,有10个样点综合生态风险等级处于中等以上。不同河段风险预测结果见表6。
表6 柳江流域沉积物重金属元素潜在生态风险评价结果
Table 6
河段 | RI | ||||||||
---|---|---|---|---|---|---|---|---|---|
As | Cd | Cr | Cu | Hg | Ni | Pb | Zn | ||
金城江 | 11.23 | 266.33 | 1.12 | 3.31 | 58.15 | 2.58 | 16.95 | 2.70 | 362.37 |
柳江 | 9.31 | 102.90 | 1.07 | 4.10 | 28.83 | 2.89 | 5.15 | 1.16 | 155.41 |
龙江 | 7.02 | 87.52 | 1.12 | 3.71 | 22.51 | 2.35 | 4.87 | 0.97 | 130.07 |
罗秀河 | 3.63 | 13.88 | 0.89 | 3.17 | 8.83 | 2.13 | 4.87 | 0.45 | 37.84 |
洛清江 | 3.39 | 26.49 | 0.86 | 3.21 | 17.12 | 2.25 | 2.72 | 0.51 | 56.56 |
融江 | 14.69 | 53.26 | 1.11 | 4.35 | 25.38 | 2.91 | 4.36 | 0.93 | 106.99 |
平均 | 8.46 | 100.89 | 1.05 | 3.79 | 28.17 | 2.63 | 6.18 | 1.17 | 152.34 |
4.4 因子分析
分析前,应用Kaiser-Meyer-Olkin(KMO)和Bartlett球形方法对水系沉积物数据进行相关矩阵检验,检验结果显示Bartlett显著性均小于0.01,在99%的置信区间内可靠,数据符合结构效度。
柳江水系沉积物元素可提取3个特征值大于1的主成分,提取方法采用主成分法,旋转方式为凯撒正态化等量最大法,累积贡献率为83.769%,基本包含了16个变量的全部信息,特征值分析结果见表7。
表7 柳江流域沉积物元素含量主成分分析结果
Table 7
公因子 | 方差贡献 | 元素 | 因子 | 元素 | 因子 | ||||||
---|---|---|---|---|---|---|---|---|---|---|---|
特征值 | 贡献率 | 累积贡献率% | PC1 | PC2 | PC3 | PC1 | PC2 | PC3 | |||
1 | 8.942 | 55.885 | 55.885 | As | 0.64 | 0.57 | -0.16 | Zn | 0.30 | 0.92 | 0.14 |
2 | 3.293 | 20.584 | 76.469 | Cd | 0.17 | 0.94 | 0.02 | P | 0.46 | 0.01 | 0.58 |
3 | 1.168 | 7.301 | 83.769 | Cr | 0.84 | 0.31 | 0.31 | S | -0.08 | 0.61 | 0.67 |
4 | 0.703 | 4.394 | 88.163 | Cu | 0.90 | 0.07 | 0.26 | Mn | 0.86 | 0.19 | -0.12 |
5 | 0.574 | 3.585 | 91.748 | Hg | 0.25 | 0.90 | 0.13 | Al2O3 | 0.94 | 0.05 | 0.14 |
6 | 0.337 | 2.107 | 93.855 | F | 0.25 | 0.90 | 0.13 | Fe2O3 | 0.90 | 0.27 | 0.26 |
7 | 0.268 | 1.677 | 95.533 | Ni | 0.92 | 0.21 | 0.21 | CaO | 0.02 | 0.70 | 0.11 |
8 | 0.243 | 1.519 | 97.052 | Pb | 0.14 | 0.86 | 0.32 | Corg | 0.46 | 0.38 | 0.74 |
注:提取方法为主成份分析法;旋转法为具有Kaiser 标准化的正交旋转法;旋转在4次迭代后收敛;下划线表示所占因子较大的元素,即某成份的主要组成元素
第一因子贡献率55.85%,表征因子为As、Cr、Cu、Ni、Mn、Al2O3、Fe2O3,其中Cr、Al2O3属造岩元素;As、Cu、Ni、Mn、Fe2O3为亲铜及亲铁元素,自然背景下,其含量变化主要取决于岩石侵蚀、风化、变质等物理化学过程,特别是Mn、Ni,在表生作用的土壤迁移转化过程中,能够保持良好、稳定的亲合性,元素含量差异性往往主要受土壤母质地质年代变化的影响,如Mn被认为是来自古生代和中生代等最古老的地质岩石中,Ni来自于上新世松散沉积等地质区[18]。与广西桂中地区土壤基准值相比,沉积物中As、Cr、Cu、Ni、Mn、Al2O3、Fe2O含量均较低,富集系数仅为0.49~0.73,属贫化状态。综上所述,第一因子主要为自然背景因子,反映地层时代、岩石矿物、地形地貌以及侵蚀风化作用等要素的影响。
第二因子由Cd、Hg、F、Pb、Zn、Cao组成,特征贡献率20.58%。Cd、Hg、Pb、Zn属重金属元素,元素含量分布与人类活动存在密切关系,高值点皆分布于金城江流域河池市下游段,受人类污水排放影响。该区域矿产资源丰富,周边已探明的矿场有54个,矿点有32个,包括铅、锌、锑、锡等有色金属矿,沉积物元素含量与矿产开发“三废”排放存在一定关系。该河段沉积物重金属元素含量较高,其中Cd含量为1.04×10-6,为区域土壤地球化学基准值的2.3倍。前人研究认为,Zn是工业活动的标识元素,其不但用于压铸合金(汽车、轻工业、机械行业)、电池业、印染业、橡胶业、化学工业等,还与其他合金金属广泛应用于电镀、喷涂等行业[19]。因此,综合分析认为第二因子应与人为因素有关,主要反映工矿产业对元素含量的影响。
4.5 源解析
4.5.1 重金属元素污染来源识别
表8 柳江流域沉积物重金属不同来源贡献率
Table 8
元素来源 | As | Cd | Cr | Cu | Hg | Ni | Pb | Zn | |
---|---|---|---|---|---|---|---|---|---|
自然源 | 0.69 | 0.39 | 0.68 | 0.73 | 0.79 | 0.78 | 0.36 | 0.45 | |
人为源 | 工矿业 | 0.14 | 0.44 | 0.07 | 0.07 | 0.02 | 0.04 | 0.52 | 0.33 |
农业 | -0.14 | -0.09 | 0.09 | 0.03 | 0.03 | 0.01 | 0.24 | 0.07 | |
不确定源 | 0.31 | 0.26 | 0.26 | 0.17 | 0.16 | 0.17 | 0.15 | ||
R2 | 0.75 | 0.91 | 0.89 | 0.87 | 0.84 | 0.94 | 0.88 | 0.95 |
研究区不同元素自然源、人为源占比差异明显:As、Cr、Cu、Hg、Ni以自然源为主,贡献率占比68%~79%,人为源贡献均小于20%,不确定源贡献率15%~30%;Zn人为源占比约为40%,其中工矿业源占比33%,农业源占比7%,不确定源占比约15%;Cd、Pb主要以人为源为主,占比分别为44%及66%,可见柳江中Cd含量较大,人为因素起相当重要的作用。通过APCS-MLR模型定量分析重金属来源贡献率,其结果与主分成分析、聚类分析结果一致。
4.5.2 Cd来源空间分布特征
根据前文污染现状评价及潜在风险评价结果,柳江水系中主要污染影响元素及生态风险元素皆为Cd,其主要来源贡献率空间分布见图3。
图3
Cd自然源贡献率高值区主要分布于龙江流域,融江流域柳城县周边分布面积亦较大,该区域属石炭系及二叠系石灰岩地区,侵蚀风化程度较大,成土母质风化为Cd的主要来源;低值区则分布于洛清江及罗秀河区域,丘陵地貌,成土母质多为泥盆系黏土岩,侵蚀风化程度较低,符合受河流冲积、地形地貌影响的结论。Cd工矿业源贡献率高值区主要分布于金城江周边以及柳江柳州市区段,前者位于河池市下游段,受周边矿山开采冶炼的影响较大,后者处于城市周边,与工业生产排放存在一定的相关性;低值区及低背景面积较广,分布区域人类活动均较弱。
5 结论
1) 研究区水系沉积物中各元素空间离散程度较大。重金属元素含量分布表现出较强的地域特征:高值区分布于金城江流域,低值区集中于罗秀河、洛清江流域。
2) 研究区水系沉积物污染等级以轻度及偏中度为主,Cd为主要污染因子。潜在生态风险等级为中低等级,潜在风险因子主要为Cd,其次为Hg。金城江水系为柳江流域中应予以重点关注的河段。
3) 不同重金属元素其自然源、人为源占比差异明显:As、Cr、Cu、Hg、Ni、Zn均以自然源为主,但Zn人为源占比较大,约为40%;Cd、Pb则主要以人为源为主,占比分别为44%、66%;Cd自然源贡献率高值区主要分布于龙江流域,工矿业源贡献率高值区主要分布于金城江周边以及柳江柳州市区段。
参考文献
上海市小城镇河流沉积物重金属污染特征研究
[J].
Pollution character of heavy metals in river sediments from small towns,Shanghai
[J].
Assessment of spatial distribution and potential ecological risk of the heavy metals in relation to granulometric contents of Veeranam lake sediments,India
[J].
北运河源头区沙河水库沉积物重金属污染特征研究
[J].
Characteristics of heavy metal pollution in the sediments from Shahe reservoir,the upper reach of the north canal river
[J].
海河流域中南部河流沉积物的重金属生态风险评价
[J].
Ecological risk assessment of heavy metals in surface sediments in the southern and central Haihe river
[J].
柳江流域生态系统服务价值的影响研究
[J].
The study of ecosystem service values of Liujiang River Basin
[J].
Index of geoaccumulation in sdiments of the Rhine River
[J].
应用地质累积指数评价沉积物中重金属污染:选择地球化学背景的影响
[J].
Applying geo-accumulation index to assess heavy metal pollution in sediment:influence of different geochemical backgrounds
[J].
Spatial variations of heavy metals contamination in sediments from Odiel River(southwest Spain)
[J].
An ecological risk index for aquatic pollution control:A sedi-mentological approach
[J].
Ecological risk assessment of soil
[J].
A auantitative assessment of source contribution to inhalable particulate matter pollution in metropolitan boston
[J].
Source apportionment of PM2. 5 in Beijing using principal component analysis absolute principal component scores and UNMIX
[J].
Spatial-temporal variations, sources, and transport of airborne inhalable metals (PM10) in urban and rural areas of northern China
[J].
Speciation and potential remobilization of heavy metals in sediments of the Taihu Lake, China
[J].
Source apportionment of PAH in hamilton harbour suspended sedi-ments: Comparison of two factor analysis methods
[J].
Spatio-temporal patterns and source apportionment of coastal water pollution in eastern Hong Kong
[J].
日照市土壤重金属来源解析及环境风险评价
[J].
Sources identification and hazardous risk delineation of heavy metals contamination in Rizhao city
[J].
Distribution of some elements in surface soil over the Kavadarci region,Republic of Macedonia
[J].
泉州城市表层土壤中金属元素来源分析
[J].
Sources of metals in the urban topsoil of Quanzhou city
[J].
菜茶果园红壤微生物量磷与土壤磷以及磷植物有效性之间的关系研究
[J].
Study on relationships among microbial biomass P, soil P and plant-availability of P in red soils
[J].
Response of soil microorganisms to the addition of carbon, nitrogen and phosphorus in a forest Rendzina
[J].
不同施肥处理对水稻土颗粒有机碳与磷的影响
[J].
Effects of different fertilization regimes on organic carbon and phosphorus status in particulate fractions of paddy soil
[J].
基于PCA/APCS和地统计学的广州市土壤重金属来源解析
[J].
Source apportionment of soil heavy metals in Guang-zhou based on the PCA/APCS model and geostatistics
[J].
基于受体模型和地统计学相结合的土壤镉污染源解析
[J].
Source apportionment of soil heavy metal Cd based on the combination of receptor model and geostatistics
.[J].
/
〈 |
|
〉 |
