天津市郊区土壤重金属的污染评价与来源解析

师荣光，张又文，3，许萌萌，郑向群，赵宗山*

（1.农业农村部环境保护科研监测所，天津 300191；2.中国科学院青岛生物能源与过程研究所，生物基材料重点实验室，青岛266101；3.沈阳农业大学土地与环境学院，沈阳 110866）

土壤质量与农业生产、食品安全、人类健康以及生态安全息息相关[1]。作为土壤中主要的污染物质之一，土壤重金属具有生物积累性、持久性、遗传毒性及难降解性等特点[2]，对人类健康和生态系统的可持续发展造成极其不利的影响，因而受到了国内外学者的广泛关注与研究。除成岩源外，土壤重金属主要来源于工业和城市废物排放、采矿活动、化肥和农药的不当使用以及污水灌溉等[3-4]。已有文献报道，截止2014年，美国地区有70%以上的土壤（数十万采样点）遭受到污染物的侵害[5]，而中国同样也有数万公顷的耕地被重金属污染[5-6]。因此研究土壤重金属的含量特征及其污染状况（包括污染水平和生态风险）对土壤重金属的有效防控非常重要。

天津市地处中国的“首都经济圈”，作为直辖市、国家中心城市、环渤海地区经济中心，工业、农业、运输业较发达，人口密度较大[7-8]。自20世纪中叶后，天津市郊农田土壤长期使用污水灌溉，使得天津成为我国五大污灌区之一；另外随着工业的发展和城市的扩张，天津市郊土壤重金属污染日益严重[9]。西青区位于天津市西南部，是我国最大的副食品生产基地之一，蔬菜、渔业是当地的支柱产业。滨海新区位于天津东部沿海地区、环渤海经济圈的中心地带，拥有世界吞吐量第四的天津港，工业企业较为集中[10]，经济发展潜力巨大。快速的工业化、农业集约化及城市人口的迅速增长加速了重金属对土壤和环境的污染[5]。许萌萌等[11]对天津郊区农田土壤重金属进行了含量分析、污染及来源评价和雨水沥浸探讨，但样品数量相对较少，且集中于单一的农田土壤，难以代表研究区的综合土壤状况。

为了全面调查天津市郊土壤的综合污染状况，本文以天津市郊区——西青区和滨海新区（包括原大港区、汉沽区、塘沽区）为研究区域，采集了农田、道路绿化带、民宅及厂区周边等多种用地类型的286个土壤样品，旨在：（1）对研究区土壤重金属的含量与空间分布进行分析；（2）以新颁布的、较为严格的农用地土壤重金属风险筛选值[12]和天津土壤重金属背景值[13]为基准，进一步开展对研究区土壤重金属的污染水平和生态风险的评估；（3）通过多元统计分析方法，根据7种重金属之间的内在相关性，探讨重金属的可能来源，以期为研究区土壤重金属监控管理提供参考。

1 材料与方法

1.1 研究区域概况

天津市属温带季风气候，四季分明，降雨主要集中在每年的6—8月，年平均降雨量为577.8 mm[11]。天津地区土壤类型属于潮土，pH平均值为7.81，有机质平均值为2.25%。采样区域约2000 km2，介于北纬38.56°～39.33°和东经116.91°～118.05°（图1）。为了使样品能够更好地代表研究区域内整体的土壤状况，我们选取了多种用地类型的土壤进行采集，包括农田（耕地、园地、林地）、芦苇地、鱼塘/水塘/水库/河岸边土地、道路绿化带、民宅周围土地、油田/工厂/盐场附近土地等，并根据各位点所处位置的用地类型，将其分为农业用地、居民用地和工业用地，并进行分类探讨。

1.2 样品采集与预处理

土壤样品的采样时间为2017年7月中旬至8月上旬。采样点的设置采用网格布点和随机布点相结合的方式，通过土壤采集器收集0～20 cm深度的表层土壤。为了提高样品的代表性，将每个采样单元采集的5份土壤样品组成一个混合样品，每个混合土样的质量为1 kg。采集的土壤样品总数为286个，其中农业、居民和工业用地的样品分别为154、89个和43个。将样品分别标记并包好后，低温转移至实验室。土样去除植物残体、砂砾等组分后，经室温风干、研磨、过100目尼龙筛后，密封常温保存待用。

1.3 样品分析与质量保证

将定量的土壤样品与混酸（HCl∶HNO3∶HClO4∶HF=1∶2∶1∶1）充分混合后，放置于微波消解仪中消解，消解溶液经赶酸后，于0.45 μm的滤膜过滤，定容至一定体积，利用ICP-MS（Agilent 7900）测定溶液中7种重金属（Cr、Ni、Cu、Zn、As、Pb和Cd）的含量。土壤重金属的测定分析采取平行样和标准物质加入法（国家土壤标准样品GSS-8、GSS-23和GSS-27）进行质量控制。Cr、Ni、Cu、Zn、As、Pb和Cd的回收率分别为72%～87%、81%～94%、84%～89%、80%～95%、90%～131%、85%～109% 和 90%～106%；As的精密度为7.52%，其余6种重金属的精密度均小于5%。

图1 研究区地理位置及采样点示意Figure 1 Location of the study area and the sampling sites

1.4 土壤重金属污染评价

1.4.1 地质累积指数

地质累积指数[14-15]又称地积累指数（Geo-accumulation index），是以自然地质过程造成的背景值为基准，评估人为活动对重金属污染的影响，见公式1。本文采用1990年中国环境监测总站编写的《中国土壤元素背景值》中天津地区7种重金属的土壤背景统计值[7]（平均值，见表1）。

式中：Igeo为重金属的地质累积指数；Ci和Bi分别为重金属i的实测值和背景值，mg·kg-1；1.5为补偿系数，通常表示由成岩效应引起的自然波动[16]和轻微的人为输入[15]。地质累积指数可分为7个级别，Igeo＜0、0≤Igeo＜1、1≤Igeo＜2、2≤Igeo＜3、3≤Igeo＜4、4≤Igeo＜5和Igeo≥5；污染级别为0～6级，分别表示无污染、无污染到中度污染、中度污染、中度污染到强污染、强污染、强污染到极强污染和极强污染。

1.4.2 污染负荷指数

土壤中重金属的污染负荷指数（Pollution load index，PLI）涵盖土壤中7种重金属的污染因子，可以用于评估土壤重金属的整体污染状态，是一种土壤重金属污染简单、有效、综合的表征手段[2,17]。PLI的计算公式见公式（2）和公式（3）。

式中：Cif为重金属i的含量因子；Ci和Bi同公式（1）；在本文中，n=7。污染负荷指数可分为4个级别，PLI≤1、1＜PLI≤2、2＜PLI≤3、PLI＞3，分别表示无污染、中度污染、重度污染和极强污染。

1.4.3 内梅罗综合指数

土壤中7种重金属的风险水平可以通过单因子指数法（Uni-index，Pi）（公式4）和内梅罗综合指数法（Nemerow Index，Pc）[11]（公式5）进行评价。

式中：Pi为重金属i的单因子指数；Ci为土壤重金属含量，mg·kg-1；Si为农用地土壤重金属风险筛选值（GB 15618—2018），mg·kg-1（当Ci≤Si时，则认为该土壤污染物对农产品质量安全、农作物生长或土壤生态环境的风险低）[12]。

式中：Pc为重金属的内梅罗综合指数；Pmean为单因子指数平均值；Pmax为单因子指数最大值。内梅罗综合指数分为 5 个级别，Pc≤0.7、0.7＜Pc≤1.0、1.0＜Pc≤2.0、2.0＜Pc≤3.0和Pc＞3.0，分别表示安全、警戒、轻度污染、中度污染和重度污染。

1.4.4 生态风险指数

土壤重金属的生态风险指数可以用于评估生物群落对有毒金属的敏感性及重金属引起的潜在生态风险[2]。本文采用单一的生态风险因子（公式6）和综合的生态风险指数（Ecological risk index，RI）（公式7）进行评价。

式中：Ei为重金属i的生态风险因子；Ti为重金属i的毒性因子；Cif为重金属i的含量因子（同公式2）。Cr、Ni、Cu、Zn、As、Pb和Cd的毒性因子分别为2、5、5、1、10、5和30[18]。生态风险因子分为5个等级，Ei＜40、40≤Ei＜80、80≤Ei＜160、160≤Ei＜320 和Ei≥320，分别表示低风险、中度风险、较高风险、高风险和超高风险。生态风险指数分为4个级别，RI＜65、65≤RI＜130、130≤RI＜260和RI≥260，分别表示低风险、中度风险、较高风险和超高风险[2，18]。

1.5 数据分析

利用Origin 9.0软件进行三类用地类型的土壤重金属含量图的绘制；利用Surfer 8.0软件进行土壤重金属含量经纬度分布图的绘制；利用IBM SPSS statistics 20软件进行土壤重金属含量的独立样本T检验和多元统计分析（包括相关性分析、主成分分析和聚类分析）。

2 结果与讨论

2.1 土壤重金属含量分析

表1为研究区域土壤重金属的含量统计值。采样点各重金属含量与北京郊区农田土壤相比偏高或相当[19]。与天津污灌区相比，Pb的含量值高一倍，Cr、Ni和Cu的含量相当，Zn和Cd的含量偏低[20]。以上结果表明京津地区土壤重金属含量比较接近，差异不明显。与重金属的农用地土壤污染风险筛选值[12]对比可知（表1），除As外，其他6种重金属基本处于风险筛选值之下。以该值为参考值，5种金属超标率排序情况为As（18.53%）＞Cd（2.10%）＞Zn（1.05%）＞Ni（0.70%）＞Cu（0.35%）。这表明研究区土壤中Cr、Ni、Cu、Zn、Pb和Cd对农用地来讲，风险较低，处于较安全的水平；而As元素的超标情况较为严重。

为了直观地区分不同用地类型中土壤重金属的差异，我们将土壤用地类型分为三类，如表1所示，三类用地类型的重金属含量并无明显差异。通过对重金属进行两两的独立样本检验来鉴别三类重金属是否存在显著差异，该显著性差异检验结果同样显示，三类土壤类型的土壤重金属之间不存在显著性差异（表1）。对超标较为严重的As来讲，三类用地类型的As超标率依次为居住用地（28.09%）＞工业用地（25.58%）＞农业用地（11.04%）。

表1 土壤中7种重金属的含量统计值、对应的土壤背景值和风险筛选值及显著性差异概率Table 1 Content statistics values,background values and risk screening values of seven heavy metals in soils

2.2 土壤重金属的空间分布

图2天津市郊区土壤中7种重金属的空间分布图Figure 2 Spatial distributions of seven heavy metals in soils from the suburb of Tianjin

图2 为研究区土壤中7种重金属含量的空间分布图。以农用地风险筛选值为黄色填充色，低于该值显示为灰色和绿色，高于该值则显示为粉色和红色。由图2可以看出，研究区域内未见明显的点源污染。与其他6种重金属相比，As污染情况较为严重（图2e），且污染区域主要集中在研究区的南部（原大港区）和东北部（原汉沽区），表现为较为明显的面源污染。大港区为工业区，该地聚集了许多工业企业（包括大港油田），汉沽区邻近港口且农田较多，农药使用[21]、污水灌溉[22]以及频繁的工业活动[23]，导致当地土壤中As污染较为普遍。土壤Ni、Cu、Zn和Pb在个别采样点存在超标污染情况，这可能主要是受到人为活动的影响[24]。土壤Cr和Pb在整个研究区内未出现超标情况，表明其对该区域土壤生态系统风险较低。

2.3 土壤重金属污染评价

土壤污染状况通常以对应的地学背景值或环境标准值为参考进行评估。本文采用多种评价手段（包括地质累积指数、污染负荷指数、内梅罗综合指数和生态风险指数）对土壤重金属污染进行综合、全面的评价。

2.3.1 地质累积指数

研究区土壤中7种重金属的地质累积指数（Igeo）平均值排序为 Cd（0.22）＞As（0.10）＞Pb（-0.12）＞Zn（-0.38）＞Ni（-0.73）＞Cu（-0.76）＞Cr（-0.88）（表2）。土壤中As和Cd的Igeo平均值介于0和1之间，而其他5种重金属的Igeo平均值小于0，这表示As和Cd在整个研究区的土壤中处于无污染到中度污染的水平，Cr、Ni、Cu、Zn和Pb在土壤中处于无污染的水平。各重金属的Igeo值主要集中在无污染（Igeo＜0）水平（Cr、Ni、Cu、Zn、As、Pb和Cd分别为98.60%、97.90%、92.31%、80.77%、43.36%、67.83%和35.31%）和无污染到中度污染（0＜Igeo＜1）水平（Cr、Ni、Cu、Zn、As、Pb和Cd分别为1.40%、1.40%、6.99%、16.78%、44.06%、31.82%和52.80%）。综上所述，从整体情况看，半数以上的土壤采样点中As和Cd有新的人为污染源的加入，同时可认为其他5类重金属未对土壤造成明显的人为污染[2]。

2.3.2 污染负荷指数

土壤的污染负荷指数（PLI）涵盖土壤中7种重金属的污染因子，可以更加直观、综合地评估研究区内的土壤污染水平。研究区内土壤的PLI范围为0.34～2.59（均值1.21±0.32），表明该区域土壤整体处于中度污染的水平。根据PLI的分类依据，大部分的土壤样品（74.48%）处于中度污染的水平，24.48%和1.05%的土壤样品则分别处于无污染和重度污染的水平。

农业、居民和工业用地的PLI范围分别为0.55～2.54（均值 1.24±0.30）、0.62～1.86（均值 1.18±0.28）和0.35～2.59（均值1.14±0.43），该数据显示三类用地类型对土壤的污染影响不明显。研究区内土壤PLI的空间分布见图3a，PLI高值点主要位于原塘沽区和西青区：塘沽某园区内（2.59）、小枣树林（2.54）、油田边荒地（2.07）、西青大沽排污河旁某农田菜地（1.97）。高值点污染物主要来源于工业活动、污水渗漏、农药化肥的使用等。从整体来看，研究区域污染水平差异不大，这表明当地没有明显的点源污染，大气沉降和降雨可能是该地区土壤重金属的主要传输方式[25-28]。

2.3.3 内梅罗综合指数

内梅罗综合指数以农用地风险筛选值[6]为参考值，可以用于评估单一重金属对土壤的风险水平，结果如表2所示。研究区土壤重金属内梅罗综合指数排序为 Cd（2.22）＞As（1.56）＞Zn（1.49）＞Ni（1.04）＞Cu（0.87）＞Cr（0.72）＞Pb（0.45）。结果显示研究区内土壤Cd处于中度风险水平；As、Zn和Ni处于轻度风险水平；Cu和Cr处于警戒水平；Pb处于安全水平。对比三类用地类型土壤的内梅罗污染指数结果，农业用地的Cd处于中度风险水平，居民用地和工业用地分别处于安全和警戒水平；研究区内土壤中的其他6种重金属基本处于相同或相近的水平，结果再次证实三类用地类型对土壤重金属污染无明显差异。该结果与上述2.3.1和2.3.2所得结论基本一致。

表2 土壤中7种重金属的地质累积指数、内梅罗综合指数和生态风险因子Table 2 Geo-accumulation indexes,Nemerow indexes and ecological risk factors of seven heavy metals in soils

2.3.4 生态风险指数

经过对土壤中单一重金属对土壤的风险评价后，本文还采用生态风险指数进行了重金属对土壤生态风险的整体评估。土壤生态风险指数（RI）范围为39.63～660.81（均值为100.87±57.76），表明研究区土壤重金属基本处于低风险水平（92.31%）。另外，有6.30%、0.70%和0.35%的位点分别处于中度风险、较高风险和超高风险。土壤RI值的空间分布见图3b，RI高值点主要为耕地、林地，且Cd含量较高，这可能是杀虫剂、杀菌剂等的不当使用造成的[25]。引起土壤潜在生态风险的重金属元素主要为Cd，其Ei平均值为62.74（表2），表明土壤中Cd主要处于中度风险水平。其中有27.62%和55.59%的土壤采样点Cd处于低风险和中度风险水平；另外，13.99%、1.75%和1.05%的采样点中Cd对土壤具有较高风险、高风险和超高风险。全部位点的Cr、Cu、Zn和Pb的Ei值均小于40，因此可认为这4种重金属对土壤具有很低的生态风险。综上所述，7种重金属中，Cd对土壤造成的生态风险最为严重，因此需要重点监控。

2.4 重金属的来源解析

为了更好地监测、控制重金属对土壤的污染，本文采用多元统计分析法（相关性分析、主成分分析和聚类分析）探究重金属之间的关联，并以此推测土壤重金属的来源与途径。

2.4.1 相关性分析

重金属之间的相关性分析通常作为其来源鉴别的依据之一[29-30]。根据土壤中重金属的Spearman相关性分析（表3）可以得出，除As外，Cr、Ni、Cu、Zn、Pb和Cd两两均具有显著的相关性（P＜0.01），表明土壤中这6种重金属具有相同或相似的来源[29]。As元素较为特殊，与Cr呈现显著相关的关系，与Ni、Pb不相关，与Cu、Zn、Cd则呈现负相关的关系。该结果与天津某区域农田重金属的相关性分析结果较为一致[11]，表明两区域重金属空间变化影响因素可能相同，污染来源相似或较为相关。

2.4.2 主成分分析

为了进一步明确研究区土壤中7种重金属的来源，本文进行了重金属的主成分分析。首先，我们对数据进行了KMO和Bartlett检验，KMO值为0.673，Bartlett值为719.186，P值＜0.05，表明数据基本适合主成分分析。主成分分析得到了两个提取因子，成分1以 Zn（0.999）、Pb（0.573）、Cu（0.542）和 Cd（0.512）为主导，其整体贡献率为67.01%；成分2的贡献率占比为21.67%，主要反映的是Ni（0.883）、Cr（0.872）和As（0.200）的作用（图4a）。上述分析可以得出，Pb、Cu和Cd可能具有相同的来源，Zn与这3种金属来源相似；As、Cr和Ni均为成岩元素[31]，且Ni和 Cr的污染源更为相近，As与Ni和Cr的来源差异较大。该结果与重金属的相关性分析结果较为相近。

2.4.3 聚类分析

图3 土壤重金属的污染载荷指数（PLI）和生态风险指数（RI）的空间分布图Figure 3 Spatial distributions of pollution load index and ecological risk index of heavy metals in soil

为了更全面、可靠地分析土壤中7种重金属的关系，本文还进行了重金属的系统聚类分析（图4b）。本研究中，7种重金属被分为了三簇，第一簇为Cr和Ni，第二簇为Zn、Pb、Cu和Cd，As单独为第三簇。距离簇可以用于表示重金属之间的关联程度，距离簇上的值越低，关联就越显著[32]。第一簇中Cr和Ni的距离值（值为1）非常小，表明两者来源非常相似，有着显著的关联；As与前两簇的距离值（值为25）很大，表明As与其他6种重金属的关联小，来源差异较大。该结果与相关性分析和主成分分析具有较好的一致性。综合3种统计分析方法可以得出，Cr和Ni具有相同或极其相似的来源；Zn、Pb、Cu和Cd具有较为相似的来源；As与其余6种重金属的来源差异较大。已有文献中介绍天津农田中As与其余6种重金属来源差异较大，不同的是As主要为自然源[11]。

研究区西青区是天津市主要的蔬菜产区之一，且具有多年的污泥填埋和污水灌溉历史[33]。由于传统污水处理方法无法有效去除污水中的重金属，使得大量重金属通过灌溉方式进入田地。滨海新区为北方制造业、国际航运中心和国际物流中心，拥有世界第四大吞吐量的天津港，聚集着大量的工业企业和丰富的油田、天然气资源。大量的工业生产和运输活动会不可避免地对土壤环境产生一定的污染。Cr和Ni为成岩元素，且在土壤中的含量与当地的土壤环境背景值接近，空间分布结果显示Cr和Ni未出现或仅存在个别点的超标情况，表明土壤中Cr和Ni主要为自然源（成岩过程）；除Cu外，Zn、Pb和Cd在土壤中的含量均高于环境背景值，空间分布结果显示Cu、Zn和Cd存在个别采样点的超标情况，表明除自然源以外，土壤中Zn、Pb、Cu和Cd还有少量人为源的输入（包括工业活动、交通运输[25]、城市来源[26]和农业活动[31，34]）。As的含量明显高于背景值，且采样区域土壤中As的含量超标情况较为严重，表明除成岩过程外，As的主要来源包括农业化学品的使用[35]、污水灌溉、工业生产等人为因素。

3 结论

（1）研究区内土壤重金属未见较为明显的点源污染，且用地类型对土壤重金属的存在水平无明显影响。

（2）4种土壤重金属污染评价方法表明研究区土壤重金属整体处于中度污染和低风险水平；Cd和As处于轻度或中度的污染水平、低等或中等风险的水平；Cr、Ni、Cu、Zn和Pb则整体处于无污染和低风险水平。

表3 土壤中7种重金属含量的Spearman相关性分析Table 3 Spearman correlation analysis of seven heavy metals in soils

图4 土壤中7种重金属的主成分图（a）和聚类分析的系统树图（b）Figure 4 Principal component diagram（a）and dendrogram（b）of seven heavy metals in soils

（3）多元统计分析结果显示，Cr和Ni主要来源为自然源；Zn、Pb、Cu和Cd受自然源和人为源的共同影响；As与以上6种重金属的来源差异较大。

（4）相较于其他5种重金属，Cd和As对研究区土壤的污染贡献最大，含有人为来源的Cd和As的土壤样品数量占总样品的一半以上，因此需要对当地土壤中的Cd和As进行重点监控。