中国农田土壤重金属空间分布特征及污染评价

环境修复网讯:土壤是人类赖以生存和发展的重要资源之一,也是陆地生态系统重要的组成部分。近年来, 随着我国城市化进程加快，矿产资源开发、金属加工冶炼、化工生产、污水灌溉以及不合理的化肥农药施用等因素导致重金属在农田土壤中不断富集[1~4]。重金属作为土壤环境中一种具有潜在危害的污染物, 通常不随水淋滤, 也不能被微生物所降解, 具有易积累、难挥发、毒性大和隐蔽性强等特点[5~7]。当土壤遭受重金属污染时，不仅会影响土壤生物的生长发育, 更重要的是重金属元素还会通过食物链以及皮肤接触等方式在人体中积累，严重威胁人体健康[8~10]。文献[11]显示,全国土壤污染物总超标率为16.1%，其中耕地土壤点位超标率为19.4%,其主要污染物包括镉、铬、铅、铜、汞、砷以及部分有机污染物,由此可见我国农田土壤污染形势较为严峻。

鉴于目前农田土壤环境质量堪忧以及重金属污染的长期危害性，近年来农田土壤重金属污染成为国内外学者研究的热点。当前国内研究多集中在典型耕地[12~14]、蔬菜地[15, 16]、果园[17]和工矿区[18, 19]等小尺度区域土壤重金属的富集特征及空间分布, 而对中国农田土壤重金属的污染现状进行系统地评价分析相对较少， 如Wei等[20]汇总了我国部分城市及农田土壤重金属数据, 并对这些数据展开分析评价；宋伟等[21]根据138个典型耕地土壤重金属数据库计算了我国耕地土壤重金属的污染概况, 结果显示全国耕地土壤重金属的污染率为16.67%；张小敏等[22]采用平均值法对我国农田土壤铅、镉、铜、锌和铬这5种重金属的富集情况及空间分布开展研究, 发现我国土壤重金属的空间分布具有明显区域特征；陈奕云等[23]基于文献计量法对全国34个省级行政区农田土壤重金属的污染状况进行分析评价,结果显示全国Pb和Cd的污染较为严重。总体上全国农田土壤重金属的空间分布研究存在以个别样点值代替区域平均值、样点偏少或未展开有效评价等问题, 难以全面反映中国农田土壤重金属的实际分布特征水平。因此，本文基于2002年以来可收集到的中国各行政区农田土壤重金属实测含量数据, 采用地理信息系统(GIS)及地统计法探讨我国各行政区农田土壤不同重金属：Cr(铬)、Cd(镉)、Pb(铅)、Cu(铜)、Zn(锌)、As(砷)和Hg(汞)的污染情况及空间分布特征, 同时采用地累积指数法对其污染状况展开评价。通过准确掌握中国各行政区农田土壤重金属的实际污染状况, 以期为建立农田土壤空间数据库、全面了解中国农田土壤重金属污染现状、分布特征以及对未来农田污染制定有效的防治策略提供科学的依据。

1 材料与方法

1.1 数据来源与预处理

中国农田土壤重金属含量统计数据均通过“某市(县)重金属”关键词搜索，同一地区以最新年份文献为准。排除城镇工业用地、大气降尘以及海洋沉积物等非农田土壤重金属研究，共选出2002年以来公开发表的文献603篇，包括614个典型农田样点重金属实测数据(由于文献资料有限，统计地区不包含西藏自治区、中国香港、澳门特别行政区、中国台湾以及南海诸岛)，具体的样本点分布见图1。其中样本点坐标为文献中直接标注的采样点或通过反转研究区地理坐标位置获取, 重金属数据通过提取文献中对应的农田土壤重金属含量或根据污染指数反推得到。总体上各文献多采用随机采样法，选择在土壤表层(0~10cm或0~20cm)的深度范围内取样，所有样品经室内风干、压碎、过筛和消解后，用电感耦合等离子质谱法以及原子荧光分光光度法等方法测定各样本重金属含量，并且保证样品分析测试的平均偏差小于10%。

西藏自治区、中国香港及澳门特别行政区、中国台湾和南海诸岛的数据暂缺，下同

图 1 农田土壤重金属样本分布示意

由于部分研究中采样点受矿区和工业区等因素影响较重, 可能导致区域某种重金属含量值存在异常, 而异常值对变异函数的模型拟合具有一定的影响。因此为保证分析结果的可靠性, 在展开各项分析前需对重金属样本数据进行检验处理以剔除异常样本，本文采用3倍标准差法处理各文献数据，以X±3S为标准, X表示某重金属所有样点的平均值，S为标准差，超出该范围的数据则需剔除[24]。结果显示,土壤Cr、Cd、Pb、Zn、Cu、As和Hg的异常值数分别为10、14、7、7、9、5和9个。

1.2 研究方法

1.2.1 普通克里金插值

由于土壤重金属的富集程度在区域内存在空间分异, 且部分地区无统计资料，简单利用样本数据展开分析评价难以准确反映土壤实际状况。克里金法(Kriging)实质上是利用区域化变量的原始数据以及变差函数的结构特点去估计无数据区域的变量值，该方法能最大限度地利用样点信息来确定未知样点的估计值，这不仅考虑到无数据点与已知样点的空间位置, 而且还能利用各邻近样点间的位置关系做出线性无偏的最小方差估计, 使评价结果更精确和更符合实际[25,26]。结合前人在土壤重金属空间分布特征方面的研究[27~29]，本文选取克里金算法中的普通克里金法(Ordinary Kriging)确定中国农田土壤重金属的空间分布特征。

1.2.2 土壤污染评价

土壤污染评价方法采用国内外普遍使用地累积指数法(Igeo),该法是由Muller在20世纪60年代提出的一种能够定量评价土壤重金属污染水平的方法, 地累积指数不仅能充分反映土壤各元素的自然分布特征和累积程度, 而且还能判断土壤环境受自然背景和人类活动的影响[30~32],其计算公式为：

式中,Igeo指土壤重金属i的地累积指数；Ci指土壤重金属i的实际含量(mg·kg-1)；Si指土壤重金属i的评价标准值(mg·kg-1), 评价标准采用各省、直辖市和自治区的土壤背景值；1.5是考虑各行政区土壤岩石差异引起背景值变动的修正系数。地累积指数分级标准见表 1, 土壤重金属i的Igeo值越大, 污染越严重；当Igeo值大于0时, 表明土壤重金属主要来源于人类活动而非成土母质等自然因素。

表 1 土壤地累积指数分级标准

1.3 数据分析

利用SPSS 19.0软件对土壤重金属样本数据进行描述性分析和正态分布检验, 地统计分析在GS+9.0软件中完成,克里金插值及空间分布图的制作在ArcGIS 10.2软件中完成。

2 结果与分析

2.1 土壤重金属含量的描述性统计

剔除异常值后的中国农田土壤重金属描述性统计结果见表 2, 其中土壤环境质量标准采用《土壤环境质量农用地土壤污染风险管控标准(试行)》(GB 15618-2018)中农用地土壤污染风险筛选值的较严格项。结果显示, 农田土壤Cr、Cd、Pb、Cu、Zn、As和Hg的样本平均值分别为59.97、0.240、32.73、28.91、86.52、10.40和0.118 mg·kg-1, 其中全国土壤Cr和As的平均含量均未超出中国土壤背景值,而Cd、Pb、Cu、Zn以及Hg的平均含量分别超出中国土壤背景值2.47、1.26、1.11、1.17和1.82倍, 此外各元素的平均含量均未超过土壤环境质量标准。变异系数可以反映区域重金属元素的分布差异,农田土壤7种重金属元素的变异系数大小顺序为：Hg>Cd>Pb>As>Cu>Zn>Cr,均属中等程度变异(0.10<CV<1.00)[33], 其中Cd和Hg的变异系数均超过0.75, 表明Cd和Hg受外界干扰比较大, 含量分布不均匀, 变异性较强；其余5种重金属元素的变异系数较低, 说明它们受外源影响相对较小，从偏度和峰度值来看, 土壤Cd和Hg的偏度和峰度值较高, 表明土壤Cd和Hg呈现高累积状态。

2.2 正态分布检验及半方差函数模型

普通克里金插值是区域化变量的线性估计, 它假设数据呈正态分布, 认为区域化的变量是未知的, 因此在插值前需对样本数据进行正态分布检验, 而对于非正态分布序列则需通过转换才能得到相应序列，正态分布检验的结果表明(表 3), 除Cr元素外, 其余6种重金属含量均偏离正态分布, 因此需采用对数方法进行标准化处理。经检验, 转换后的数据均符合正态分布, 在此基础上采用半方差函数模型计算中国农田土壤7种重金属元素的半方差函数值, 在各向同性的情况下根据残差平方和(RRS)最小以及决定系数(R2)最大的原则获取最佳空间变异函数理论模型[34, 35]。结果显示，只有Cd元素符合球形模型(Spherical)，其余6种重金属元素均为指数模型(Exponential)。从块基比看, 7种重金属元素的C/(C+C0)值均大于0.5，较好反映了各自空间变异的结构特征,为空间拟合的最佳模型。

2.3 中国农田土壤重金属空间分布特征

利用ArcGIS 10.2软件分别对经过预处理的全国各省、直辖市和自治区的土壤重金属含量数据进行普通克里金插值, 得到中国农田土壤重金属含量的空间插值图(图2)。结果显示,云南、四川、上海以及福建和江西交界有大范围的Cr高值区,甘肃、陕西、山西和湖北等地出现次高值区, 而青海和海南两省土壤Cr含量较低；Cd在甘肃中部、新疆和江苏北部、云南、广西与贵州三省交界、河南与湖北交界、安徽与江西交界以及湖南出现高值, 表明这些区域可能存在明显的Cd污染源, 全国其他地区土壤Cd含量较低且分布较为平均；总体上南方土壤Pb含量明显高于北方, 其中以新疆农田土壤的Pb含量最低, 而辽宁西部、陕西中部、重庆、云南、湖南、安徽以及福建等地则有较大范围的高值区；土壤Cu含量在空间上分布较为均匀, 仅在新疆、内蒙古中部、陕西、河南、天津、云南以及湖南等地含量较高；土壤Zn在全国多地出现连片高值区, 包括四川、广西、贵州、湖南、福建、广东、浙江以及云南北部, 东北三省、内蒙古、宁夏和陕西北部地区土壤Zn含量较低；As的空间分布较为复杂, 云南、四川和贵州三省交界土壤As含量最高, 湖南、广西、广东中部、青海南部、新疆和黑龙江北部以及内蒙古东部存在次高值区, 四川中部以及华北和华东地区的土壤As含量较低；Hg的高值区主要分布在浙江东部、广东中部以及福建、湖南和贵州, 其他省份土壤Hg含量相对较低，从空间分布特征情况来看，土壤Cr、Pb、Cu、Zn和As的空间分布特征明显, 而土壤Cd和Hg在空间分布上十分相似。