1. 引言
土壤作为人类赖以生存的物质基础 [1],对人类社会生存发展具有重大的意义。土壤质量的好坏直接关系到人类的生存状况,土壤质量问题正在不断得到世界范围内的关注,土壤质量的保持和提高已日益成为当今国际土壤学界、环境科学界、地学界共同关注的课题。
近年来,在高密度人口活动下土壤重金属元素相关领域研究成果丰硕,这些成果主要涉及土壤重金属元素的污染与防治、区域土壤重金属分布、大学校园土壤重金属污染以及土壤重金属元素来源等方面。如Ahemad等 [2] 对植物体和土壤中重金属特征进行了研究,分析了重金属元素的污染程度和来源进行了分析和评价;宋智萍等 [3] 对山西师范大学校园内土壤Pb、Cr、Zn、Cu和Ni五种重金属元素的污染状况进行了研究分析,发现这五种重金属元素都有不同程度的富集但没有达到污染的程度;刘红侠等 [4] 对徐州市北郊农业土壤重金属污染情况研究分析,发现其农业土壤5种重金属平均含量均远高于该市农业土壤背景值,Zn为污染警戒等级;朱伟等 [5] 对镇江城市道路沉积物中重金属污染的来源进行分析,发现重金属主要来源与交通污染和生活污染源有关,通过CER值评价RDS中重金属的人为污染程度,Pb、Zn、Cr、Ni和Cu均表现出一定程度的人为影响;周广柱等 [6] 对兖矿地区的土壤环境质量进行分析,研究结果表明,一部分地区表土中镍、镉等重金属含量较高,与矿山距离较近的地区均受到重金属不同程度的污染,并有可能对兖矿地区农业生产构成威胁。
宿州学院教育园区位于宿州市东部,地处宿州火车站与宿州高铁站之间,交通便利。基于土壤重金属超标的危害性和近年来土壤重金属污染不断恶化的现实情况以及国家建设生态文明的重大战略部署,本文旨在研究与探讨高密度人口活动下土壤重金属元素的特征及其与人为活动的联系,通过对各个功能区中土壤各种元素成分分布情况进行分析,从而有效的推测出土壤中元素与人类高密度活动的联系,提出建议用以控制以及降低土壤中有害元素,提高城市土壤质量并为土壤环境的保护、污染防治和生态建设提供参考或帮助。
2. 材料与方法
2.1. 研究区简介
研究区为宿州学院教育园区,坐落在安徽省宿州市区东约10 km处,校园北邻天地阳光学习学校,西侧和南侧均为农田和村庄,东侧为学府大道,学府大道南接宿灵公路北连宿州大道,交通便利。该校区在校师生及职工约2万人,研究区面积约51.43万平方米,人口密度约3.89万人/平方千米,属于人口密集区。如图1,校园整体呈南北走向的长方形,师生公寓、食堂及浴室等生活建筑共20栋,总体分布于校园西北侧;教学楼、实验楼及体育场馆等教学场所20余处,总体分布于校园东南侧;二冬大道和闵贤路呈南北向贯穿校园。校园内地表水系有珍珠湖和柳溪,珍珠湖面积约7760 m2,深2~3 m,北与柳溪以涵洞相连通;柳溪源于宿州市郊沱河支流自西向东注入小黄河,东西横穿校园中部约5000 m。
2.2. 样品采集与测试
选取高密度人口活动的宿州学院教育园区为研究区,依据研究分析需求在教学楼、主干道、河流和湖泊、学生宿舍、食堂等附近布置采样点23个。采样点分布见图2。采样时使用干净的塑料铲将表层杂物除去,垂直地将土壤切开,采取0~20 cm深度的新鲜地表土壤0.5~1 kg装入密封样品袋,分别编号1~23,用铅笔注明采样地点、日期和采样名称,且采样时避免接触铁质物品和使用铁质器械。
待样品采集完毕,将土样摊开于干净的塑料布上,在通风、干燥和避免阳光照射的室内自然风干。样品风干时,需用一张标签纸放在样品中防止混匀。风干过程中要经常翻动土样加速干燥以及用不锈钢镊子剔除土壤以外的动植物残体,而后对风干后的土样按照研究分析要求进行破碎、研磨,研磨为粉状后过200目的尼龙筛,充分混匀后,装入密封袋备用。制备好的样品要妥善贮存,避免日晒、潮湿和高温的污染。注意每次处理样品前必须对使用过的尼龙筛、相关研磨工具进行清洁。
将预处理好的样品进行压片,先将模具安装好,土样取两勺放于模具中,轻压使其平整;再用硼酸粉末充填模具周围,轻轻垂直拿起模具避免土样被带出;然后放入垫片和T型圆柱。将模具整体慢慢的放置于压片机上,旋紧螺旋,关紧气阀,加压至20 MPa,停留2分钟左右;慢慢打开放气阀,使压力缓慢下降到0,拧开螺旋取出模具,打开底座,压出模块,小心取出压好的样片。注意每次处理土样前均要使用酒精清洁模具整体。处理好的样品用XRF荧光光谱仪(型号:EDX5500H)进行检测,测试精度为0.01 ppm。测试过程中,每测10个土壤样品需测试一次标样校正。
2.3. 统计与分析
样品检测结果统计如表1所示。由统计可知,所测土壤中3种元素的含量Mn含量最高,含量范围385.19~698.07 mg∙kg−1,平均含量为492.43 mg∙kg−1;Zn含量范围46.1~116.22 mg∙kg−1,平均含量为63.12 mg∙kg−1;Hg含量最低,含量范围0.01~0.55 mg∙kg−1,平均含量为0.12 mg∙kg−1。
Table 1. Test results of soil heavy metal content on the campus of Suzhou University (mg/kg)
表1. 宿州学院校园土壤重金属含量测试结果(mg/kg)
如表2所示,3种元素含量的变异系数Mn最大,达到0.41;其次为Zn元素的变异系数为0.23;Hg的变异系数最小,为0.01。Mn变异系数较大,则说明该元素在表层土壤中的空间分布差异较大,受到高密度人口活动的影响比较大;Hg的变异系数较小,则说明该元素在表层土壤中的空间分布差异较小,受到高密度人口活动的影响极小;Zn元素则介于前两者中间。但是,不能片面地只通过变异系数的大小来推测各种元素受高密度人口活动影响的程度,应该考虑到其他的影响因素以及从多方面去分析考量。
Table 2. Statistics of test results
表2. 检测结果统计表
注:中国土壤元素背景值源于文献 [7],土壤环境质量标准源于文献 [8]。
综合表1和表2,Mn元素的含量均值略高于其背景值,Hg元素的含量均值是其背景值的三倍,Zn元素的含量均值略低于其背景值。一般来说,某个元素的正异常越大、异常点越多,则其受高密度人口活动影响的程度越大。
根据《土壤环境质量标准》规定的土壤环境质量二级标准即为保障农业生产,维护人体健康的土壤限制值,Mn、Hg和Zn的元素含量均值都在土壤环境质量二级标准限制值以下。
2.4. 污染评价方法
2.4.1. 地质积累指数法
地质累积指数(Igeo)通常称为Muller指数,不仅考虑了自然地质过程造成的背景值的影响,而且也充分注意了人为活动对重金属污染的影响,因此,该指数不仅反映了重金属分布的自然变化特征,而且可以判别人为活动对环境的影响,是区分人为活动影响的重要参数。但该方法不能评价某一种土壤重金属元素含量分布的综合污染指数或多种土壤重金属元素含量分布的综合污染指数,只适用于分析某研究对象中单独的土壤重金属元素污染程度。公式如下 [9]:
(1)
式中:Cn为样品中元素n的浓度;Bn为背景浓度;1.5为修正指数,通常用来表征沉积特征、岩石地质及其它影响。如表3所示,地质累积指数可分为7个级别,级别越高污染程度越大。
Table 3. Geological accumulation index (Igeo) grading criteria
表3. 地质累积指数(Igeo)分级标准
Mn、Hg和Zn的地质累积指数(Igeo)计算结果见表4,统计结果见表5。Mn元素的所有测试值污染级别全为0级,皆为无污染;Hg元素无污染占60.87%,无污染–中度污染占8.70%,中度污染占13.04%,中度污染–强污染占13.04%,强污染占4.35%,4号样品测值为强污染;Zn的所有测值主要以无污染为主,无污染–中度污染次之。22号样品采样位置为珍珠湖东北向的河床,该点的Mn的地积累指数达到最大值,为−0.05。4号样品采样位置为三食堂,该点的Hg的地积累指数达到最大值,为3.20;地积累指数级别仅次于4号的是19号样品采样位置为柳溪懿行桥正下方的河床,该点的Hg的地积累指数将近3。Zn的地积累指数在21个测点测值上为负,只有18号和19号测点样品的地积累指数超过零。有14个测点的Mn、Hg和Zn的地积累指数级别皆为0级。
Table 4. Results of geo accumulation index calculations
表4. 地积累指数计算结果
Table 5. Statistical table of geo accumulation index results
表5. 地积累指数结果统计表
通过使用地积累指数对数据的分析,可以发现在三食堂采样的4号的Hg的地积累指数级别最高,人为活动大的18号、19号等采样点的地积累指数级别较高,位于校园偏僻地方、非人车流量大的15号、16号等采样点的地积累指数级别较低。
2.4.2. 内梅罗综合指数法
只反映不同种类的单个污染物的污染程度是不行的,要完全和充分地反映土壤整体的污染程度。因此就需要将各个重金属元素的污染指数结合起来进行评价,综合污染指数兼顾了单因子污染指数平均值和最高值,可以突出污染较重的重金属污染物的作用。但同时可能会夸大或缩小一些重金属污染物的影响作用,导致对土壤综合污染程度分级的不精确,对土壤元素含量分布特征的分析也会有偏差。内梅罗综合指数法是目前被广泛使用的一种综合指数法 [10]。内梅罗综合指数法计算公式为 [11]:
(2)
式中:P是采样点的综合污染指数,Pimax为i采样点地积累指数中的最大值,
为地积累指数加权平均值。公式如下 [12]:
(3)
式中:
为权重,Pi为单项测量值。根据Dalway J. Swaine [13] 对重金属的分类,Mn、Hg和Zn的权重都为1。土壤综合污染程度分级标准见表6。
Table 6. Classification standard of soil comprehensive pollution degree
表6. 土壤综合污染程度分级标准
如表7所示,23个样品综合污染程度1级安全的有16个,该16个样品Mn和Zn的地积累指数级别皆为0级,且Hg的地积累指数级别大部分为0级;2级污染警戒的只有一个,为22号样品;3级轻度污染的有4个;4号和19号样品为4级中度污染,其Hg元素的地积累指数级别都较高。
Table 7. Calculation results and statistics of Nemero composite index
表7. 内梅罗综合指数计算结果与统计
内梅罗综合指数的分析表明,可以发现在三食堂采样的4号样品的Hg元素的土壤综合污染程度最强,校园河流柳溪的18号、19号等采样点的土壤综合污染程度较强,位于校园偏僻地方、非人车流量大的8号、15号和16号等采样点的土壤综合污染程度较弱。
由此可见,虽然地积累指数法与内梅罗综合指数法的侧重方向和计算方法有所不同,但是这两种方法的分析结果可以达到相互验证的效果。
3. 结果与分析
3.1. 描述性统计
首先对23个土壤样点的3种重金属元素的检测结果(表1)和检测结果统计表(表2)进行分析,了解样品数据的分布状况。变异系数 [13] 反映了采样总体中各样点之间平均变异程度。其中,样点之间重金属元素含量差异最大的是Hg (变异系数为1.33),最小的是Mn (变异系数为0.14),3种重金属元素含量的总体平均变异程度由小到大依次为Mn、Zn和Hg。而且3种重金属元素含量的极差均比较大,表明土壤重金属元素含量存在着明显的空间分布差异 [14]。
3.2. 空间分布特征
克里金(Kriging)插值法是利用区域化变量的原始数据和变异函数的结构特点,对未采样点的区域化变量的取值进行线性无偏最优估计的一种方法,也是一种在许多领域都很有用的地质统计格网化方法,采用普通克里金插值,对不符合正态分布的采取反距离加权插值进行空间可视化表达 [14] [15]。Surfer8.0基于windows操作系统的三维绘图软件,可以轻松制作基面图、数据点位图、分类数据图、等值线图、线框图、地形地貌图、趋势图、矢量图以及三维表面图等 [16]。Surfer8.0软件自带有克里金插值功能,则在确定克里金插值模型后,利用Surfer8.0软件绘制出土壤重金属元素Mn、Hg和Zn的空间区域分布点位图、等值线图,清晰地呈现出各重金属元素含量的空间分布情况。土壤重金属含量的空间分布特征如图3~5所示。
由图3可知Mn相对富集于20号、22号和23号采样点所围成范围内和以及2号、10号和17号采样点附近;在2号、4号、5号、9号、12号、13号和16号采样点所在的区域出现相对贫化。由图4可知相对富集于8号和19号采样点所围成范围内以及4号和18号采样点附近;在研究区的西部、西南、
Figure 3. Spatial distribution of Mn elements (mg∙kg−1)
图3. Mn元素空间分布(mg∙kg−1)
Figure 4. Spatial distribution of Hg elements (mg∙kg−1)
图4. Hg元素空间分布(mg∙kg−1)
Figure 5. Spatial distribution of Zn elements (mg∙kg−1)
图5. Zn元素空间分布(mg∙kg−1)
西北区域以及中部的3号、9号、10号、17号和21号采样点出现相对贫化。由图5可知Zn的空间分布相对于Mn和Hg总体上较集中,在研究区中部区域含量较高,相对富集于17号和19号采样点所围成范围内。由图3~5可知,17号采样点的Hg元素和19号采样点的Mn元素与另二者相比都呈现出明显的富集与贫化的两极差异。
3.3. 人为活动与元素富集的关系
根据国家土壤元素背景值和土壤环境质量标准进行多方面分析,全部样点重金属元素含量低于国家土壤环境质量二级标准,总体上受外源物质影响和富集程度较低。Hg元素的环境质量指数异常较大,且变异系数较大,说明某些样点受高密度人口活动影响较大,富集较为明显。而且Mn、Hg和Zn元素含量较高的样点主要集中在人口密度较大,交通活动比较繁忙的功能区附近。
通过土壤重金属的描述性统计分析、地积累指数法、内梅罗综合指数法和元素空间分布特征对Mn、Hg和Zn在研究区内不同区域的污染程度及相对富集程度进行了分析,土壤环境质量总体上较好,有轻微的污染程度。Mn、Hg和Zn在地表土壤中的元素含量均受到高密度人口活动下一定程度的影响,高密度人口活动对Mn、Hg和Zn的影响都较大。
4号采样点位于三食堂,样品土壤潮湿,土质疏松,Hg元素含量最高,所以Hg出现较强的负异常;18号和19号采样点位于校园河流柳溪,样品实际为河流底泥,所以Hg出现较强的正异常。Mn、Hg和Zn均表现为正异常的9、10、11、12号采样点等处,9、11、12号采样点日常人流量和车流量都很大,且都位于停车位附近;10号采样点位于操场,空间相对较小,人流量较大。几乎全部采样点的Mn和Zn元素的污染程度较小,均表现为正异常。
4. 结论
(1) 研究区地表土壤重金属元素含量的统计结果和空间分布特征表明,重金属Hg元素含量的变异系数和极差均较大,空间分布差异显著。
(2) 从地累积指数法内和梅罗综合指数法可以看出3种重金属元素的污染程度顺序为Hg > Zn > Mn,宿州学院教育园区的地表土壤已受到不同程度的重金属污染,处于轻度污染水平,校园重金属污染问题亟需解决;
(3) 对宿州学院教育园区校区地表土壤中Mn、Hg和Zn含量的检测结果进行了污染评价分析,发现地表土壤中重金属元素含量与分布特征会受到高密度人口活动的影响,不同元素受高密度人口活动影响的程度也存在差异,Mn、Hg和Zn三种元素受人为活动的影响都较大,一般表现为重金属元素富集。
基金项目
安徽省大学生创新创业训练项目(201910379162)、安徽省高校学科(专业)拔尖人才学术资助(gxbjZD2020091)共同资助。
参考文献
NOTES
*第一作者。
#通讯作者。