生态环境学报 2018, 27(8): 1547-1555 www.jeesci Ecology and Environmental Sciences E-mail: editor@jeesci 基金项目：国家自然科学基金项目（41606053）；广东省自然科学基金项目（2016A030310364）；广东省教育厅青年创新人才项目（2016KQNCX098）；岭南师范学院人才引进专项项目（ZL1607）
作者简介：罗松英（1985年生），女，讲师，博士，主要从事土壤重金属污染、地球化学等研究。E-mail:luosongying@163 收稿日期：2018-04-06 及生态风险评价
罗松英，王嘉琦，周敏，叶嘉雯，柯思茵，陈东平，赖晓清
大丈夫的私房钱岭南师范学院地理系，广东 湛江 524048
摘要：近年来随着沿海地区工农业的迅速发展，大量污染物汇集于河口与海湾区，使得红树林湿地重金属污染也日趋严重。为探究东海岛红树林湿地表层（0~15 cm ）土壤重金属污染现状，于2017年6月沿东海岛海岸线典型红树林分布区布设了7个站位（由东南向西北方向依次为YC 、SYC 、JJC 、XWC 、
TMHHK 、DCMT 和DTSD ）进行采样，测定样品中As 、Cd 、Cr 、Cu 、Hg 、Ni 、Pb 、Zn 8种重金属元素质量分数并分析其空间分布特征，采用潜在生态危害指数进行污染及潜在生态风险评价，结合相关性分析和主成分分析方法探讨其重金属来源。结果表明，（1）东海岛红树林湿地表层土壤8种重金属的质量分数平均值均超过雷州半岛土壤环境背景值，但未超过国家土壤环境质量一级标准限值。（2）单因子污染指数显示，东海岛红树林湿地表层土壤各重金属元素总体上到达中等污染程度。（3）潜在生态风险评价结果显示，东海岛红树林湿地表层土壤重金属污染总体上属于中等潜在生态风险；其中Hg 潜在生态风险参数最高，平均值为105.524，属于强潜在生态风险，其余7种重金属元素均属于轻微潜在生态风险。（4）从潜在生态风险空间分布特征看，东海岛红树林表层土壤重金属空间差异较明显，其中XWC 站位的潜在生态风险指数为351.690，为强潜在生态风险（主要来自Hg 、Cd 和As 污染），YC 、TMHHK 、DCMT 和DTSD 4个站位为中等潜在生态风险，SYC 和JJC 表现为轻微潜在生态风险。（5）统计分析结果表明，东海岛红树林湿地表层土壤中As 、Cr 、Cu 、Ni 、Pb 和Zn 这6种重金属元素之间具有强相关性，结合实地调查结果推测其来源主要与水产养殖排污及农业面源污染有关；而Cd 和Hg 之间的相关性不明显，推测其来源分别与船舶污染和生活排污有关。 关键词：空间分布；生态风险评价；重金属；主成分分析；东海岛 DOI: 10.16258/jki.1674-5906.2018.08.022
中图分类号：X820.4 文献标志码：A 文章编号：1674-5906（2018）08-1547-09
引用格式：罗松英, 王嘉琦, 周敏, 叶嘉雯, 柯思茵, 陈东平, 赖晓清. 2018. 湛江东海岛红树林湿地表层土壤重金属空间分布特征及生态风险评价[J]. 生态环境学报, 27(8): 1547-1555.
LUO Songying, WANG Jiaqi, ZHOU Min, YE Jiawen, KE Siyin, CHEN Dongping, LAI Xiaoqing. 2018. Spatial distribution and ecological risk assessment of heavy metals in the surface soils of mangrove wetland in Donghai Island, Zhanjiang [J]. Ecology and Environmental Sciences, 27(8): 1547-1555.
红树林是分布在海陆交替过渡地带的湿地生态系统，受陆源包括人类活动以及海域环境信息的叠加影响（李瑞利等，2012；Liu et al.，2014）。由于红树林生境潮流速度较低，土壤沉积物具有粘粒含量高、有机质含量高、硫化物含量高、还原性强等特征，能够大量接受来自潮汐和地表径流所携带的重金属污染物，往往成为潜在的重金属富集区（程皓等，2009；Bosire et al.，2005）。重金属污染具有隐蔽性、毒性高、不易分解等特点，可通过食物链的放大作用最终危害人体健康，因此重金属污染备受关注（Nath et al.，2013）。对于整个红树林生态系统而言，土壤沉积物是重金属污染物的主要
贮存库，绝大部分的重金属沉淀固定在土壤沉积物
中（戴纪翠等，2009）。红树林土壤重金属含量的高低与其生境理化性质、土壤重金属背景值高低和受附近城市排放重金属污染的程度等密切相关（刘景春，2006；李柳强等，2008）。国内外对不同海
岸红树林重金属污染的研究结果表明，不同红树林区重金属污染程度不同，受陆域环境影响较大的红树林湿地重金属潜在生态风险有逐渐加强的趋势（郭菊兰等，2013；Sandilyan et al.，2014；Li et al.，2015；Miola et al.，2016）。
东海岛位于广东省湛江市东南部，地势平坦，
岛上红树林湿地面积约1475 hm 2，
属于湛江红树林
1548 生态环境学报第27卷第8期（2018年8月）
国家级自然保护区范围（王中丙，2013；朱小丹，2016）。自1992年东海岛开发试验区设立以来，岛上填海项目逐渐增多，潮间带区域逐渐被吞食，红树林湿地面积大幅减少（李晓敏，2008）。随着钢铁、化工、石化等项目相继落户东海岛，海岛周边生态环境压力不容忽视。同时由于东海岛及周边地区城市化和工业化进程的急剧加速，工业废水、城市生活污水、养殖排污、船舶污染和农业面源污染对生态环境造成持续性压力，生长于潮间带环境的红树林湿地重金属污染问题日益凸显（刘加飞，2013）。目前，针对东海岛红树林湿地土壤重金属的研究也少见报道，仅见刘加飞（2013）对东海岛潮间带表层沉积物重金属污染进行过研究，其研究对象并未专门针对红树林湿地生态系统，且仅对Cu、Zn、Cd、Pb、Cr这5种重金属元素进行污染评价，研究不够系统全面。
本研究通过分析东海岛红树林湿地土壤重金属含量及其空间分布特征，采用潜在生态危害指数法和多元统计法对红树林湿地土壤重金属污染现状进行评价及来源解析，以期为东海岛红树林环境污染监测与生态保护提供科学参考。
1 材料与方法
1.1 研究区概况
湛江东海岛（20°55'~21°05'N，110°11'~ 110°33'E）位于湛江湾与雷州湾之间，全岛面积约为289.49 km2，海岸线为126 km，是中国第五大岛屿（王中丙，2013）。属于北热带季风海洋性气候区，年平均气温23 ℃，年均降水雨量1800~2000
mm，长夏无寒冬，热带气旋是影响该地区的主要灾害性天气。潮汐类型为不规则半日潮，波浪以风浪为主，多年平均波高为0.9 m（陈则实，1999）。红树林主要分布在东海岛南岸、西岸以及北岸靠西段，主要红树植物为白骨壤（Avicennia marina），少量的秋茄（Kandelia candel）、无瓣海桑（Sonneratia apetala）和红海榄（Rhizophora stylosa）。西岸红树林常见分布于河口，周围有大量的水产养殖池，养殖废水直接排入红树林；南岸红树林主要分布于内湾或背风区，面向雷州湾，有少量水产养殖池或盐场。北岸有少量红树林分布，位于东参码头和对岸小岛上。自东海岛成立开发实验区以来，该地的交通、通讯、电力、海水养殖、钢铁、石化、农副产品加工、饲料、城郊型农业等临
港工业和资源型工业的规模和效益不断扩大和提高，已经成为开发区经济的规模产业（刘加飞，2013）。据调查资料显示，湛江湾和雷州湾的主要污染来源有城市生活污水、工业废水、养殖排污、船舶污染和农业面源污染（刘芳文等，2015）。
1.2 样品采集与研究方法
1.2.1 样品采集
影子价格综合考虑东海岛红树林分布面积、潮位变化及人类活动影响等因素，于2017年6月沿着海岸线在红树林典型分布区布设了7个采样站位，从东南向西北方向分别是YC、SYC、JJC、XWC、TMHHK、DCMT和DTSD。采样站位分布如图1所示，每个站位采集土壤样品2~3个，共15个样品。采样时先用塑料小铲刮去表层枯枝残叶，然后使用直径为
图1 湛江东海岛红树林湿地采样站位分布
Fig. 1 Distribution of sampling stations in mangrove wetlands of Donghai Island in Zhanjiang
罗松英等：湛江东海岛红树林湿地表层土壤重金属空间分布特征及生态风险评价 1549
负载均衡
6 cm 的PVC 管压入淤泥中采集表层0~15 cm 土壤样品，每个样点采用梅花采样法取5个土样混合，采集约1 kg 土样，所采集的样品装入密封的乙烯封口袋中，带回实验室。 1.2.2 样品处理及测试方法
野外采集的土壤样品经过自然风干、研磨、过筛（100目）处理后，测定砷（As ）、镉（Cd ）、铬（Cr ）、铜（Cu ）汞（Hg ）、镍（Ni ）、铅（Pb ）、锌（Zn ）8种重金属元素的质量分数。测试地点为广州澳实矿物实验室，其测试过程为：称取两份样品，经过王水消解和四酸（高氯酸、硝酸、和盐酸）消解处理后，分别采用电感耦合等离子体发射光谱（ICP-AES ）和电感耦合等离子体质谱仪（ICP-MS ）进行测试，综合取值为最后的检测结果（陈碧珊等，2017）。对每3个样品进行平行测定，同时测定空白样。其中ICP-AES 型号是Varian VISTA ，ICP-MS 型号是Agilent 7700x 。微量分析的标样准确度RSD 小于10%，回收率高于90%。
洛克王国!圣龙骑士1.2.3 研究方法
本研究采用潜在生态危害指数法对东海岛红树林土壤重金属进行污染及潜在生态风险评价，该方法相对快速简便，综合考虑了重金属的毒性、迁移转化规律、区域敏感性以及重金属区域背景值的差异，可以综合反映沉积物中重金属对生态环境的影响潜力，并给出危害程度定量划分方法，是适合异源沉积物之间评价的有效手段（邓利等，2014）。潜在生态危害指数法计算公式如下：
（1）单因子污染指数
阳光聊天室
f n
=/i i i
C C C （1） 式中，f i
C 为重金属i 的污染系数；C i 为重金属
i 的实测浓度含量；n i
C 为重金属的评价参比值。东海岛典型地带性土壤为砖红壤性红壤，故本研究采用广东省土壤的环境背景值中砖红壤中的重金属含量作为背景值。
（2）某一区域重金属i 的潜在生态危害系数
r r f =i i i E T C × （2）
式中，r i
E 为重金属i 的潜在生态危害系数；r
i
T 为重金属毒性响应系数（toxic response coefficient ）；r i T 反映了重金属的毒性水平及生物对重金属污染的敏感程度，其毒性响应系数见表1（徐争启等，2008）。
（3）多种重金属潜在生态风险指数
r r f RI ==m m
i i i i i E T C ×∑∑ （3）
式中，RI （the potential ecological risk index ）
为多种重金属潜在生态风险指数；r i
E 为金属i 的潜
在生态危害系数；r i
T 为重金属毒性响应系数，反映重金属的毒性水平及生物对金属污染的敏感程度。重金属污染潜在生态风险指标和分级关系如表2所示（黄建国，2007）。 1.2.4 数据处理
多元统计分析通过SPSS 22.0软件完成，采用相关性和主成分分析法进行重金属来源分析。相关性分
析是指对两个或多个具备相关性的变量元素进行分析，从而衡量两个变量因素的相关程度。通过重金属间的相关系数可更加准确地表明其来源途径的相似程度（雷凌明等，2014）。主成分分析是研究变量之间相互关系的多变量分析技术，通过对原变量的线性变换，形成了一系列新的相互独立变量，称之为主成分（孙钦帮等，2017）。在进行主成分分析之前，采用KMO 和Bartlett 法对原始数据集进行因子分析适宜性检验（李玉等，2006）。KMO 值为0.700，大于0.5；Bartlett 的球形度检验相伴概率为0.000，小于显著性水平0.05，表明原始数据集适合进行因子分析（姚晓飞等，2011）。
表1 重金属毒性系数
Table 1 The toxicity factor of heavy metals
Heavy metal
As
Cd
Cr
Cu
Hg
Ni
Pb
Zn
r i
T 10 30 2 5 40 5 5 1
表2 重金属污染潜在生态风险指标和分级分析
Table 2 Indices and grades of potential ecological risk index of heavy metal pollution
Enrichment coefficient of contamination (f二次开发
i C )
Enrichment pollution
degree
The potential ecological risk factor (r i
E )
Ecological risk pollution degree The Potential ecological risk index (RI) Total potential ecological
risk degree
<1 Slight <40 Slight RI<150 Slight 1≤f i C <3 Medium 40≤r i
E <80 Medium 150≤RI<300 Medium 3≤f i
C <6 Strong 80≤r i
E <160 Strong 300≤RI<600 Strong
f
i C
≥6 Very strong 160≤r i
E <320
r i
E ≥320
Very strong Fortissimo
600≤RI<1200 RI ≥1200
Very strong Fortissimo
1550 生态环境学报 第27卷第8期（2018年8月）
2 结果与讨论
2.1 土壤重金属质量分数特征及空间分布
由表3可知，As 、Cd 、Cr 、Cu 、Hg 、Ni 、Pb 、Zn 等8种元素重金属质量分数平均值超过雷州半
岛土壤环境背景值（许炼烽等，1996）；8种元素重金属质量分数平均值在《土壤环境质量标准》（GB15618—1995）一级标准范围内，但YC 、XWC 和TMHHK 站位采样区的个别元素，如As 、Hg 、Pb 质量分数超过《土壤环境质量标准》一级标准值。从重金属质量分数范围看（表3），8种重金属元素质量分数范围差异较大，As 为 6.800~26.600 mg·kg -1，Cd 为0.023~0.068 mg·kg -1，Cr 为21.000~ 74.000 mg·kg -1，Cu 为5.700~21.000 mg·kg -1，Hg 为0.037~0.200 mg·kg -1，Ni 为8.700~ 32.000 mg·kg -1，Pb 为16.000~45.800 mg·kg -1，Zn 为32.000~112.000
mg ·kg -1。
变异系数可以反映各重金属元素的离散程度，且变异系数越大说明土壤重金属污染受人为活动干扰越强烈（赵秀峰等，2010）。从变异系数上看（表3），东海岛红树林湿地表层土壤中重金属元素Pb 的变异系数最小，为33.91%，离散程度较低，反映该元素分布相对较均匀，推测目前受人为活动干扰不明显；As 、Cd 、Cr 、Cu 、Ni 、Zn 的变异系数均在40%以上，分布较不均匀，说明受到一定程度的人为活动干扰；Hg 的变异系数最大，为62.16%，分布最为不均，推测其受人为活动干扰也
较强烈。
从重金属质量分数空间分布上看，7个采样站位重金属质量分数差异较大，选取土壤中超标状况比较严重的4种重金属元素As 、Hg 、Pb 和Zn 进行分析，结果如图2所示。其中As 质量分数平均值在7个
采样站位中均高于雷州半岛土壤环境背景值，以西部靠近水产养殖地的XWC 和TMHHK 2个采样站位为最高，超出国家土壤环境质量标准值范围；Hg 含量平均值在7个采样站位中也均高于雷州半岛土壤环境背景值，以西南岸XWC 采样站位为最高，超出国家标准值范围；Pb 除了在SYC 、JJC 和DCMT 这3个采样站位外，在其他4个采样站位的质量分数平均值均接近或高于雷州半岛土壤环境背景值，其中YC 和TMHHK 采样站位的Pb 质量分数平均值均超出国家标准范围；Zn 质量分数平均值在7个采样站位中均高于雷州半岛土壤环境背景值，以东南岸YC 采样站位为最高，但未超出国家标准值范围。
2.2 表层土壤重金属潜在生态风险分析及空间分布 2.2.1 表层土壤重金属污染分析及空间分布
东海岛红树林湿地表层土壤不同重金属元素单因子污染评价结果如表4所示，As 、Cd 、Cr 、Cu 、Hg 、Ni 、Pb 、Zn 8种金属元素均达到中等污染。8种重金属元素污染程度排序如下：Zn>Hg> As>Ni>Cu>Cr>Pb>Cd 。其中，西南岸XWC 站位8
表4 东海岛红树林湿地土壤重金属单因子污染指数
Table 4 The single factor pollution index of heavy metals in the surface soil of mangrove wetland in Donghai Island
Station Enrichment coefficient of contamination (f i
C )
As Cd Cr Cu Hg Ni Pb Zn YC 2.556 1.559 1.642 2.300 2.133 2.359 1.509 4.333 SYC 1.481 0.794 0.775 0.900 1.567 0.920 0.774 1.746 JJC 1.389 0.912 0.719 0.963 1.367 0.820 0.782 1.833 XWC 3.148 1.324 1.411 1.750 6.233 1.958 1.338 3.119 TMHHK 3.537 0.706 1.840 2.063 2.367 2.332 1.534 4.024 DCMT 1.889 1.324 1.024 1.375 2.600 1.239 0.991 2.381 DTSD 2.537 1.382 1.051 1.538 2.200 1.293 1.150 2.667 Average value
2.362 1.143 1.209 1.556 2.638 1.560 1.154 2.872
表3 东海岛红树林湿地表层土壤重金属元素质量分数
Table 3 The mass fraction of heavy metals in the surface soil of mangrove wetland in Donghai Island
Heavy metal
Minimum value/ (mg·kg -1)
Maximum value/
(mg·kg -1)
Average value/ (mg·kg -1)
Standard deviation/
(mg·kg -1)
Variable coefficient/%
National standard/ (mg·kg -1)
Background value/
(mg·kg -1)
As 6.800 26.600 12.500 5.478 43.82 15.000 5.400 Cd 0.023 0.068 0.039 0.016 41.03 0.200 0.034 Cr 21.000 74.000 43.466 17.458 40.16 90.000 36.140 Cu 5.700 21.000 12.500 5.011 40.09 35.000 8.000 Hg 0.037 0.200 0.074 0.046 62.16 0.150 0.030 Ni 8.700 32.000 17.200 7.984 46.42 40.000 10.980 Pb 16.000 45.800 27.000 9.155 33.91 35.000 23.400 Zn
32.000 112.000 60.800 24.763 40.73 100.000 21.000
罗松英等：湛江东海岛红树林湿地表层土壤重金属空间分布特征及生态风险评价 1551
种重金属元素的单因子污染指数最高，Hg 为很强污染程度，As 和Zn 为强污染程度；Cd 、Cr 、Cu 、Ni 和Pb 5种元素为中等污染程度。东南岸JJC 和SYC 站位8种元素的单因子污染指数普遍较低，这两个区域污染程度相近，除As 、Hg 、Zn 3种元素达到中等污染程度级别外，其余5种元素均表现为轻微污染级别。
从超标状况比较严重的As 、Hg 、Pb 、Zn 4种重金属元素污染程度空间分布看（图3），其污染程度与元素质量分数的高低表现一致。其中，As 在东海岛西部采样站位XWC 和TMHHK 污染最严重，为强污染，其他5个采样站位为中等污染；Hg 以西南岸XWC 站位最严重，为很强污染，其他6个采样站位为中等污染；Pb 除在SYC 、JJC 和DCMT 这3个采样站位为轻微污染外，其他4个采样站位为中等污染；Zn 在所有红树林采样站位中的污染程度普遍最高，其中YC 、XWC 和TMHHK 为强污染，其他4个采样站位为中等污染。
通过野外实地调查及对比土壤环境背景值发现，东海岛红树林湿地土壤重金属污染程度与土壤重金属背景值高低及周围人为活动排放重金属程度密切相关。采样站位XWC 和TMHHK 中这4种重金属元素
污染均较为明显，主要与临近居民点和养殖鱼塘有关；XWC 站位野外现场可见有多处生活垃圾倾倒点及养殖鱼塘排污口，而TMHHK 站位则位于入海河口，周围分布着大量养殖鱼塘，且东海岛西部以农业和水产养殖业为主，入海河口成为各类污染物的汇集处（李柳强等，2008）。采样站位YC 的Zn 为强污染，其余3种元素为中等污染，主要与该采样区农业种植和水产养殖业有关。由此推测，东海岛红树林土壤重金属污染主要来自水产养殖排污、农业面源污染及生活排污，这与刘芳文等（2015）研究结果一致。东海岛由养殖排污、农业面源污染及生活排污引起的重金属污染问题应当引起重视，应注意避免人为直接排污及合理施用农药化肥。
2.2.2 表层土壤重金属潜在生态风险特征及空间分布
东海岛红树林湿地表层土壤重金属生态风险危害指数计算结果如表5所示。从单个重金属元素的潜在生态风险指数看，8种重金属元素中Hg 潜在生态风险指数最高，平均值在80~160之间，属
Background value
National standard
图2 东海岛7个采样站位主要超标重金属元素质量分数与国家土壤环境质量标准及环境背景值的对比
Fig. 2 Comparison of the mass fraction of the excessive heavy metals in seven sampling stations in
Donghai Island with national soil environment quality
level standard and environmental background values
510152025Y C
S Y C
J J C
X W C
T M H H K
D C M T
D T S D
ω(A s )/(m g ·k g -1)
Sampling stations
0.020.040.060.080.10.120.140.160.180.2
Y C
S Y C
J J C
X W C
T M H H K
D C M T
D T S D
ω(H g )/(m g ·k g -1)
Sampling stations
0.0
5.010.015.020.025.030.0
35.040.0Y C
S Y C
J J C
X W C
T M H H K
D C M T
D T S D
ω(P b )/(m g ·k g -1)
Sampling stations 0
102030405060708090100Y C
S Y C
J J C
X W C
T M H H K
D C M T
D T S D
ω(Z n )/(m g ·k g -1)
Sampling stations
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议