摘要:文章对密云水库上游白河流域土壤重金属元素(As、Pb、Ni、Hg、Cr、Zn、Cu、Cd)的空间分布、污染现状及来源进行分析评价。主成分分析和相关性分析结果表明:Cu、Pb和As元素污染来源于矿区金属的开采、选矿和冶炼过程;Zn和Cd主要来源于早期矿石运输汽车尾气的排放及轮胎磨损;Ni和Cr两种元素在该区域的富集主要来源于含重金属垃圾的随意堆放;采用混汞研磨法金矿提炼,导致白河流域矿区周边土壤中Hg元素的富集,造成局部区域出现污染。综合地累积指数及潜在生态危害指数法评价结果表明,Cd、Cr、Zn和Hg的富集程度较高,潜在生态危害较强,应当引起充分关注。
土壤作为生态系统的重要组成部分,是重要的“汇集体”,接纳了大约90%来自环境中各方面的污染物。随着工业化、城市化的快速发展,土壤环境污染问题日渐突出,重金属污染因其具有潜伏性、长期性和不可抗逆性以及危害性等特点而受到高度重视[1]。过量的重金属元素进入土壤不仅直接影响了土壤质量,更重要的是会伴随食物链进入人体而引发疾病[2]。金属矿山的开采与冶炼是我国土壤重金属污染的主要来源之一。密云水库上游白河流域是受金属矿产开发影响的重要区域,该流域内金属矿区分布较多,主要以铁矿为主,另有少部分金矿和铜矿。国内的一些学者研究发现,白河流域土壤中Hg、Pb和Cd超过国家土壤质量标准[3,4,5]。因此对白河流域金属矿区土壤存在的重金属污染风险进行评价显得尤为重要。从环境地球化学角度出发按照河流流域的空间分布特征进行分区采样,分析流域土壤的重金属的空间分布格局及污染情况,全面了解该地区土壤环境质量。旨在为探索该流域重金属来源和迁移转化规律提供科学依据,为白河流域环境实施治理、优化等措施提供坚实的理论支持。
1、材料与方法
1.1研究区域概况
白河流域地处北京市北部,包含延庆、怀柔、密云3个区,位于东经116°17′~117°30′,北纬40°02′~41°4′之间。白河发源于河北省沽源县,全长280km。经赤城县,于白河堡进入北京市延庆区境,在延庆区境内河道长55km,流域面积828.1km²;在怀柔区境内河道长51.7km,流域面积1280.7km²;在密云区境内,河道长20km,该段河床镶嵌在峡谷中,最窄处仅10多米,比降大,水流湍急,切割作用大,多呈“V”字型。
1.2样品采集与分析
土壤样品采集于2018年5~6月份,选择代表性土壤类型、土地利用类型和地理地貌,地形相对平坦、稳定、植被良好的地块,且采样点远离公路、村庄等污染源(300m以上)的位置。每个取样点采用GPS米级定位,记录空间坐标,共采集土壤样本340个,采样点位置分布图如图1所示。土壤样品经过室内风干磨碎,过100目尼龙网筛,经HNO3-HCl-H2O2消解后,采用原电感耦合等离子体发射光谱仪(ICPE-9820E001)分析消解液中重金属含量。其中,Cu、Zn、Ni和Pb元素采用火焰原子吸收法测定;Cd元素采用石墨炉原子吸收法测定;部分样品在100℃下,用HNO3∶HCl(10mL,1∶1V/V)消解2h,用原子荧光光谱仪(AFS830,中国吉天)测定Cr、As和Hg元素含量[6]。
图1白河流域取样点位置图
1.3评价方法及数据处理
(1)地累积指数法
地累积指数法主要用于定量评价沉积物中的重金属污染的程度[7,8]。地累积指数Igeo=log2[Csi/(K×Cin)],式中:Csi表示元素i在沉积物中的含量,Cni表示沉积物中该元素的地球化学背景值(表1),K为考虑到不同地区沉积特征、岩石地质及其他影响差异可能会引起背景值的变动而设定的修正系数(取值为1.5)。由地累积指数评价标准划分污染等级如表2所示。
表1北京地区土壤重金属的背景值(单位:mg/kg)
(2)潜在生态危害指数法
由瑞典科学家Hakanson提出的潜在生态危害指数法[9,10]是根据重金属元素性质及地球化学环境行为特点,从沉积学角度对沉积物或土壤中重金属污染进行评价的一种方法。潜在生态危害指数法为Cif=CsiCin,Eir=Tir×Cfi,RI=∑Eri,其中,Cfi为第i种单个重金属元素的污染系数;Csi是土壤或沉积物重金属元素i的实测浓度(mg/kg);Cni是重金属元素i评价的参比值(表1);Eri表示某重金属元素i的潜在生态危害系数;Tri表示某重金属元素i的毒性响应系数,本研究以Hakanson选定的重金属毒性系数(As:10,Cd:30,Hg:40,Pb:5,Cu:5,Zn:1,Ni:5,Cr:2单位:mg/kg)作为参评依据;RI是某一点沉积物多种重金属综合潜在生态危害指数,各自分级关系标准如表3所示。
表2地累积指数(Igeo)评价分级标准
表3潜在生态危害系数与污染程度之间的关系
(3)数据的处理
通过SPSS24.0进行土壤重金属基本参数统计以及多元统计分析[11,12],通过ArcGis10.2空间分析软件[13]进行土壤重金属空间分布特征分析,运用反距离插值的方法进行制作土壤重金属空间分布图。
2、结果与分析
2.1土壤重金属基本参数
白河流域土壤重金属基本参数统计结果如表4所示,该结果与陈小敏等2015年研究的密云水库上游金矿区土壤重金属含量均值进行比较,发现样品土壤中Zn元素的平均值较高,而Ni元素与Pb元素的平均值低于前人研究。重金属元素浓度值的基本参数分析表明:Zn元素在土壤中积累趋势较显著;Cu、Zn、Pb与Cd元素变异程度较大,表明其受外源因子干扰影响较大,其他元素均属于中等程度变异,表明其受外源因素影响较小。
2.2土壤地累积指数评价
通过对340份土壤样本中重金属含量计算得到其地累积指数(表5),8种重金属的平均值均小于1,其中Ni、Pb、As、Hg和Cd元素平均值均小于0,处于无污染程度,其余各元素处于轻度污染程度。土壤样本污染比例分析表明,Cr元素的污染占比最大,达92.5%,其中有67.4%的土壤样本处在轻度-中度污染程度,有18.4%的土壤样本处于中等污染程度,5.6%的土壤样本处于中等-强污染程度,并且有0.6%的土壤样本处于强污染程度;Zn元素污染比例达56.5%,52.1%的土壤样本处在轻度-中度污染的程度,3.3%的土壤样本处在中等污染程度,0.6%的土壤样本处于中等-强污染程度;Pb元素的污染比例最低,占比仅为10.6%。各元素统计结果表明:Cd和Hg元素的标准差值较大,表明这两种元素的地累积指数的离散程度较大,即两种元素的变异程度较高,说明该流域内一些地区Hg与Cd元素污染较为严重。
表4土壤重金属基本参数统计描述(单位:mg/kg)
表5土壤重金属地累积指数(Igeo)
2.3土壤重金属潜在生态危害风险评价
土壤重金属潜在生态危害系数(Eri)分析表明(表6),除Cd元素以外,其他7种重金属的潜在生态危害系数平均值均小于40,均属于轻度生态危害。340份土壤样本中,Ni元素的潜在危害系数值均小于40,表明该元素属于轻度生态危害,对区域内的土壤潜在危险贡献率较低;Hg元素的平均值均接近中等生态风险危害;Cd元素的平均值处于强度生态危害。按照污染程度分级,Hg、Cd元素潜在生态危害系数达到中等强、很强的比例分别为19.6%、7%、1.4%,24.6%、5.6%、3.6%。研究表明Hg元素与Cd元素在区域内的土壤生态危害的贡献率较高,该结论与陈小敏等之前的研究具有一定的相似性。
潜在危害指数分布图(图2)显示,重金属污染较为严重的区域主要集中在白河支流琉璃河区域内,该地区有较多的金矿与铁矿区,其中,白河上游千家店地区也出现一定区域的重金属污染。
表6土壤重金属生态危害系数(Eri)
2.4白河流域土壤重金属的空间分布
运用ArcGIS10.2软件反距离插值法对8种重金属含量分布进行空间分析(图3),结果表明:Cd、Zn、Pb和Cu4种重金属元素浓度分布较高的区域比较集中,且污染范围较小,Cd与Zn元素主要集中在白河流域的南部地区,Cu与Pb元素主要富集在工作区的西部;As、Ni、Hg和Cr4种重金属元素的污染区域较为广泛,工作区内大部分土壤都受到不同程度的污染,Ni和Cr元素富集区域主要集中在白河流域的东南部地区,As元素富集区域集中在工作区的西南部。从图3中可以看出,Cu、Pb与As3种元素的浓度较高的区域较为相似,区别在于As元素在工作区内的大部分区域均富集一定含量,这可能与农业用地喷洒农药有关;Zn与Cd两种元素浓度含量较高的区域大致相同;Cr与Ni两种元素分布区域大致相同;Hg元素分布明显区别于其他元素,该元素含量较高的区域主要集中在金矿区附近。
图2白河流域土壤重金属潜在危害指数分布图
2.5土壤重金属含量的相关性分析
由土壤中重金属元素之间的相关性可以推测重金属可能的来源,若重金属间具有显著的相关性,说明其可能存在相同的来源,否则来源不同,国内外多名学者采用该方法进行溯源分析[14,15]。表7为研究区域内土壤重金属之间的相关情况。结果表明Cr与Ni元素相关系数为0.907(P<0.01),说明两元素的主要来源具有一致性;Cu与Pb、As元素在土壤中呈现中等程度以上的相关性(P<0.01),其中Cu、Pb、As3种元素属于亲硫元素,容易与土壤中S形成溶解度极小的硫化物沉淀或者被铁锰氧化物及氢氧化物所吸附;Zn与Cd两种元素具有显著的相关性,达0.9以上,表明Zn与Cd两种元素具有相似的来源。
表7土壤重金属含量之间的相关系数
*在0.05级别(双侧),相关性显著;**在0.01级别(双侧),相关性显著
图3白河流域重金属元素空间分布(单位:mg/kg)
2.6土壤重金属含量主成分分析及来源识别
国内外诸多学者[16]使用主成分分析对重金属来源进行鉴别,利用SPSS24.0软件对数据进行KMO检验,得到的统计量值为0.602,Bartlett球度检验相伴概率为0.000,因此,本研究中的数据可以做因子分析,结果如表8所示。
首先,Cu、Pb和As元素在第一主成分上表现较高的载荷,且Cu、Pb和As元素具有显著的相关性,相关系数均在0.5以上,Cu、Pb、As3种元素的空间分布较为相似,主要集中在白河流域的西部和南部地区,该地区分布较多的铜矿与铁矿,所以推测这3种元素有相同的来源,且3种元素含量较大值均分布在矿区附近,3种元素含量的平均值均略高于土壤背景值,说明受到人为因素的影响。因此,Cu、Pb和As元素污染可能来源于矿区金属的开采和运输。
第二主成分表现在Zn和Cd具有较高的正载荷,且Zn与Cd元素的相关性显著,相关系数在0.9以上,重金属的空间分布结果表明,Zn与Cd元素均分布在白河流域工作区的南部,说明二者具有相同的污染来源。Cd与Zn元素的平均含量均高于土壤背景值,说明其受到人为因素的影响。研究表明,常见的来自于交通源的重金属有Zn、Cd、Pb等,矿山开采后需要车辆运输,所以Zn和Cd可能主要来源于汽车尾气排放及轮胎磨损。
第三主成分表现在Cr和Ni具有较高的正载荷,从相关性分析中可以看出,二者具有显著的相关性,且Cr与Ni两种元素均大多集中在白河流域的东南部和北部地区,推测两种元素具有相同的污染源。研究表明,废旧金属垃圾等均会造成Cr与Ni污染,白河流域东南部和北部地区分布较多的村落,所以Cr、Ni元素的污染很可能来源于含重金属垃圾处理不当。
第四主成分表现在Hg元素含量有较高的正载荷,第四主成分贡献率较小,特征值小于1,若只提取三个主成分,则第一主成分Hg分布较高的正载荷,但Hg与Cu、Pb、As元素的相关性系数不高,说明Hg元素与Cu、Pb、As元素可能都来自于大气沉降与雨水径流,但污染具体来源不同。黄兴星等[3]研究指出,密云水库上游地区曾广泛采用混汞研磨法进行金矿的采选,这造成了选厂周边土壤中Hg元素不同程度的富集。白河流域工作区内分布较多的金矿,如南部地区杨树底下金矿、崎峰茶金矿,中部地区后安岭金矿尾矿库等,重金属Hg元素的空间分布结果证实了这一结论。
表8土壤重金属含量主成分分析结果
3、结论
(1)综合地累积指数法和潜在生态危害指数法评价的结果表明,研究区域内Hg、Cr、Zn和Cd重金属元素富集程度较高,应当引起充分关注。潜在生态危害指数分布图表明白河流域内琉璃河地区重金属富集较为严重,且面积较大。
(2)重金属空间分布、相关性分析与主成分分析结果,表明Cu、Pb和As元素污染可能来源于矿区金属的开采、选矿和冶炼过程,Zn和Cd可能主要来源于汽车尾气排放及轮胎磨损,Cr、Ni元素的污染很可能来源于含重金属垃圾处理不当,土壤中Hg元素的严重富集主要是金矿采选所采用的混汞研磨法导致。
(3)选矿方法、矿山开采及运输、堆砌尾矿库等对矿区土壤重金属造成了较为严重的富集,应加强对密云水库上游白河流域相关矿区土壤质量监测及尾矿库的修复治理,避免对流域内土壤及下游密云水库水质造成污染。
参考文献:
[1]韩平,王纪华,冯晓元,等.北京顺义区土壤重金属污染生态风险评估研究[J].农业环境科学学报,2015,34(1):103-109.
[2]孙慧,毕如田,郭颖,等.广东省土壤重金属溯源及污染源解析[J].环境科学学报,2018,38(2):704-714.
[3]黄兴星,朱先芳,唐磊,等.北京市密云水库上游金铁矿区土壤重金属污染特征及对比研究[J].环境科学学报,2012,32(6):1520-1528.
[4]高彦鑫,冯金国,唐磊,等.密云水库上游金属矿区土壤中重金属形态分布及风险评价[J].环境科学,2012,33(5):1707-1717.
[5]陈小敏,朱保虎,杨文,等.密云水库上游金矿区土壤重金属空间分布、来源及污染评价[J].环境化学,2015,34(12):2248-2256.
[6]董大鹏.氢化物发生原子荧光光谱法测定土壤中的汞[J].广州化工,2009(8):176-178.
[7]何东明,王晓飞,陈丽君,等.基于地积累指数法和潜在生态风险指数法评价广西某蔗田土壤重金属污染[J].农业资源与环境学报,2014(2):126-131.
[8]尚林源,孙然好,汲玉河,等.密云水库入库河流沉积物重金属的风险评价[J].环境科学与技术,2011,34(12H):344-348.
[10]卢聪,李涛,付义临,等.基于生物可利用性与宽浓度范围的Hakanson潜在生态风险指数法的创建——以小秦岭金矿区农田土壤为例[J].地质通报,2015(11):2054-2060.
[11]陈陆望,殷晓曦,刘鑫,等.华北隐伏型煤矿地下水水化学演化多元统计分析[J].煤田地质与勘探,2013(6):43-46.
[13]陈怡先,姜小三,王勇,等.基于GIS矿区土壤重金属生态环境及人体健康评价[J].环境科学学报,2018(4):1642-1652.
[14]孙慧,毕如田,郭颖,等.广东省土壤重金属溯源及污染源解析[J].环境科学学报,2018(2):704-714.
田志君,颜常春,石国峰,梁凯旋,韩娟娟.北京市密云水库上游白河流域土壤重金属含量、来源及污染评价分析[J].安庆师范大学学报(自然科学版),2020,26(02):109-116.
基金:北京市矿山地质环境监测系统(一期)年度监测运行(白河流域示范区)项目(PXM2018_158303_000003).
分享:
半干旱风沙草原区植被覆盖度较低,风蚀沙化造成该区域土地生产力下降、土地退化等生态问题[1]。植被可减缓近地表风速,抑制沙尘释放,增加沙面地表粗糙度,从而达到防风固沙的效果[2]。不同粒度的土壤颗粒对养分的吸附性存在差异[3],细粒含量是影响土壤肥力的重要因素,其含量越多,土壤保肥性越强[4,5]。干涸湖底风蚀物主要是土壤中含有丰富养分的细粒物质,这些细粒物质会随着风沙流运移被地表植被拦截,从而形成土壤养分的富集体[6]。
2024-01-16建筑物的立面线划图是建筑物外表面投影到与其平行的投影面上得到的正投影图,其可清晰地反映建筑物的外部特征,是城市市容规划改造等工作的关键依据。因此,研究如何快速准确地采集建筑物外立面三维信息,绘制建筑物的立面线划图,对实现智慧城市数字化建设有着重要意义。
2023-08-30为适应技术发展、行业变革提出的需求,建设化工园区一体化、协同化、智能化管理模式,探索应用三维地理信息技术建设化工园区数字孪生平台,充分发挥基于服务的三维地理信息技术在地上地下空间一体化、业务协同化等方面优势,提供园区地理空间要素的地图展示、空间分析、可视化呈现等服务。
2023-08-30由于迄今为止在技术层面仍然难以准确地预测地震,为了降低地震造成的损失,评价地震可能诱发的灾害形式并进行合理的抗震设计是当前最为有效的应对方法。地震发生时,震源处释放的能量往往以地震波的形式传至地壳表面,大量震害调查结果表明,剪切波和以瑞利波(Rayleigh wave)为代表的面波更具破坏性,其传播形式如图1所示。
2023-07-10结合项目案例,介绍了国土空间生态修复中的生态破坏区域的地理信息监测技术工作思路,对技术方案的技术指标确定、图斑属性结构、技术实施路线进行了论述,并对监测成果进行了示范。实践证明,文中的经验和做法在生态修复及自然资源监测领域能起到很好的示范作用,可在其他自然资源及各部门管理中进行复制、推广应用。
2021-11-26在地矿工作中,地理信息测绘发挥着重要的作用,无论是矿山矿产考察、矿上环境防治,还是矿山开采与地质勘测等,地理信息的测量与采集都必不可少。对矿山地理信息进行全面、准确测绘是确保矿山安全的重要保障,尤其是随着地矿建设、防治等的数字化发展,地理信息测绘的重要性得到了进一步凸显。基于此,本文对测绘地理信息工作在地矿建设中的作用进行了分析、探讨。
2021-08-07随着我国社会经济的快速发展,国土空间规划工作逐渐受到人们的广泛关注和重视。地理信息大数据在国土空间规划中的应用,不仅提高了土地资源的使用率,而且推动了国土空间规划工作的顺利进行。文章以国土空间规划工作为例,分析了地理信息大数据技术在国土空间规划工作中的具体应用,以此来提高土地资源的使用效率,有效保证我国国土空间规划工作的顺利开展。
2021-06-01精准扶贫是国家扶贫战略下精准到每一个贫困人口脱贫,区别于传统方式的扶贫开发形式,基于贫困的精准数据,为每位贫困户建档立卡,做到扶持对象精准、项目安排精准、资金使用精准、措施到户精准、因村派人精准、脱贫成效精准,切实做到扶真贫、真扶贫,去除帮扶工作中存在的盲点,确保在规定时间内精准脱贫。
2020-12-21农业发展是国家发展、经济建设的基本问题,一直是人们关注的焦点。精准高效的掌握农作物空间种植结构不仅是农作物长势监测和灾情评估的基础,对于国家宏观指导农业生产、制定农业政策、规范农业管理也具有重要的现实意义。传统的农作物种植结构统计方法大多是实地调查,存在主观性强,持续时间长,耗费人力、物力、财力等问题。
2020-12-21地理国情监测统计分析是地理国情监测的重要环节,是以地理国情监测数据成果为基础,对自然、人文等要素进行统计和分析,综合反映各类型统计单元内地理国情的空间分布、空间结构、空间关系、地域差异特征。随着地理国情监测的常态化,挖掘数据价值变得越来越重要,如何提升地理国情监测和统计在自然资源管理、生态保护修复等领域中的应用水平,如何调整统计分析工作的重心是促进地理国情监测业务可持续发展急需解决的问题。
2020-12-21我要评论
期刊名称:干旱区地理
期刊人气:2077
主管单位:中国科学院
主办单位:中国科学院新疆生态与地理研究,中国地理学会
出版地方:新疆
专业分类:地质
国际刊号:1000-6060
国内刊号:65-1103/X
邮发代号:58-45
创刊时间:1978年
发行周期:双月刊
期刊开本:大16开
见刊时间:一年半以上
影响因子:0.713
影响因子:0.000
影响因子:0.000
影响因子:1.236
影响因子:0.000
400-069-1609
您的论文已提交,我们会尽快联系您,请耐心等待!
你的密码已发送到您的邮箱,请查看!