Investigation and evaluation of heavy metal pollution in cultivated land in Huishui County, Guizhou Province
-
摘要:
为调查贵州省惠水县耕地土壤重金属污染现状, 采集了5013个表层土壤样品, 测定Cd、Hg、As、Pb、Cr、Cu、Ni、Zn含量和pH值, 选用单项污染指数法和地累积指数法对农作物-土壤重金属污染现状进行了评价, 并进行相关性分析和主成分分析。研究结果表明: ①土壤重金属元素As、Cd、Cr、Cu、Hg、Ni、Pb、Zn含量超标率分别为8.74%、34.51%、3.83%、0.56%、1.56%、0.24%、0.04%、1.60%, 其中Cd(k1=1.46)、Hg(k1=1.60)均值超过贵州省土壤A层背景值; ②单项污染物指数均值从小到大顺序为Cu(0.20) < Pb(0.22) < Ni(0.26) < Hg(0.27) < Zn(0.31) < Cr(0.44) < As(0.50) < Cd(1.70), 其中Cd、As的Pi均值较大; 重金属元素地累积指数均值从小到大顺序为Cu(-1.24) < Ni(-1.18) < As < Pb(-0.99) < Zn(-0.93) < Cr(-0.71) < Hg(-0.05) < Cd(0.04), Ig污染程度较高的为Cd和Hg; ③相关性分析表明As、Hg与其他6种重金属元素之间相关性较弱, 暗示具有不同来源特征; Cd与Cr、Ni、Pb、Zn相关性较强, 表明可能具有同源性; 主成分分析表明Cd、Pb、Zn、Cr、Ni来源于自然源与交通源, 主要受到上石炭统大铺组+黄龙组、下石炭统、上泥盆统革老河组+高坡场组地层与近SN向、EW向断裂构造的控制; Cu、Ni来源于农业因子与自然源; As和Hg来源于自然源与燃煤源。
Abstract:This paper has collected 5013 surface soil samples to determine the contents of Cd, Hg, As, Pb, Cr, Cu, Ni, Zn, and pH values, and has chosen single factor pollution index method and ground accumulation index method to evaluate the status of heavy metal pollution in crop-soil and has carried out some correlation analysis and principal component analysis to study the current situation of heavy metal pollution in cultivated soil of Huishui County, Guizhou Province. The results show that: ①The excess rates of soil heavy metal elements As, Cd, Cr, Cu, Hg, Ni, Pb, and Zn were 8.74%, 34.51%, 3.83%, 0.56%, 1.56%, 0.24%, 0.04%, and 1.60% respectively. And the mean values of Cd(k1=1.46) and Hg(k1=1.60) exceeded the background values of soil layer A of Guizhou Province. ② The order of average individual Pi from small to large is presented as: Cu(0.20) < Pb(0.22) < Ni(0.26) < Hg(0.27) < Zn(0.31) < Cr(0.44) < As(0.50) < Cd(1.70), and the mean Pi of Cd and As was larger. The mean value of Ig of heavy metal elements from small to large is Cu(-1.24) < Ni(-1.18) < As, Pb(-0.99) < Zn(-0.93) < Cr(-0.71) < Hg(-0.05) < Cd(0.04), and Cd and Hg have higher Ig pollution degree. ③ The correlation analysis has showed that the correlations between As, Hg and the other six heavy metals are weak, which indicate that they have different sources while Cd is strongly correlated with Cr, Ni, Pb, and Zn, indicating that it may have homology. The principal component analysis shows that Cd, Pb, Zn, Cr, and Ni originated from natural and transportation sources, and were mainly controlled by upper Carboniferous Dapu Formation-Huanglong Formation(C1-2d+C2h), Lower Carboniferous(C1j+C1s+C1x+C1dw) and Upper Devonian Gelaohe-Gaopochang Formation(D3g+D3gp) and SN-trending and EW-trending faults. The Cu and Ni come from agricultural factors and natural sources, and the As and Hg come from natural sources and coal-burning sources.
-
Keywords:
- soil /
- heavy metal /
- investigation /
- evaluation /
- Huishui County /
- Guizhou Province
-
为有效履行中国自然资源科学管理需要,推进国土空间现代化治理,自然资源部于2020年印发了《自然资源调查监测体系构建总体方案》,提出地表基质的概念,地表基质是孕育和支撑森林、草原、水、湿地等各类自然资源的基础物质(自然资源部,2020a)。按照地表基质层、地表覆盖层、管理层三大数据层,自然资源部明确了自然资源统一调查监测的思路,并全面启动了全国地表基质调查工作。作为一项新型探索性工作,一些学者从地表基质的内涵、分类方案、调查思路、技术方法等方面进行了深入思考和研究(葛良胜等,2020;殷志强等,2020;贾磊等,2022;李响等,2023),为地表基质相关研究提供了借鉴和参考。地表基质作为地球关键带的重要组成部分,在自然资源多要素物质循环和能量转移转化中起着重要的关键纽带作用,也是植被根系所能触达的基础物质,对植被生长具有复杂的作用过程(Richter et al., 2009;Banwart et al., 2013;丁永建等,2014;张甘霖等,2021)。受成土母质、气候、地形等因素影响,地表基质结构特征、元素组成、养分含量等存在一定的空间分布差异,表现出明显的空间异质性特征,进而影响植被群落的组成、构建和生物量(李新荣,2005;邵方丽,2012;赵倩等,2020)。由于地表基质空间异质性的复杂性,目前对于地表基质与植被生长作用机理的认识依然不够(姚晓峰等,2022;杨顺华等,2024),造成山丘区土地的不合理开发、林地衰退、水土流失等一系列资源环境问题,严重制约着自然资源的可持续利用和发展(刘硕等,2019;赵辛金等,2020)。因此,充分研究地表基质空间异质性特征,了解地表基质养分组成和分布规律,阐明地表基质对植被生态的影响机理,对于进一步认识地球关键带结构及地表作用规律,支撑国土空间生态环境修复具有重要意义(马腾等,2020;殷志强等,2022)。
本文结合近年地表基质层调查工作取得的认识,以宁波山丘区地质基质为研究对象,采取剖面研究和样品测试的方法,从地表基质层空间结构异质性和元素特征异质性角度,分析了地表基质异质性产生的原因和影响因素,阐述了地表基质空间异质性对植被群落构建和空间分布的影响机理,探讨了3种典型地表基质对植被生长的适宜性条件,可为植被生态系统修复提供建议。
1. 研究区概况及研究方法
宁波市地处长江三角洲南翼,浙江省东部沿海,地理位置东经120°55'~122°16'、北纬28°51'~30°33',属于亚热带季风气候,年均降水量1480 mm,年均气温16.4℃。植被以人工针叶林、毛竹林、经济林为主,少量次生阔叶林,地带性的亚热带常绿阔叶林植被残存于边远山区(蓝雪春,2017),其中针叶林是主要植被类型,广布于中低山、丘陵山地;阔叶林是典型的地带性植被类型,主要分布在宁海、北仑、颍州等地;灌丛以壳斗科、胡桃科、杜鹃花科、金缕梅科为主。
研究区位于太平洋陆缘火山岩带的西南部浙闽粤中生代火山岩带北段,丽水−余姚深断裂以东(舒良树,2012)。自晚中生代白垩纪以来,区内火山活动十分强烈,先后经历了燕山期和喜马拉雅期2次构造岩浆活动,直至更新世仍见其活动的踪迹。岩石建造以中生代中酸性火山岩为主,其次为新生代基性火山岩(余明刚等,2021),常见岩石类型有花岗岩、玄武岩、流纹岩、砂岩等。土壤以红壤为主,主要由花岗岩、流纹岩和玄武岩风化形成。依据《地表基质分类方案(试行)》的原则(自然资源部,2020b),并结合相关研究成果(贾磊等,2022),研究区地表基质划分为4大类12个小类,其中花岗岩、流纹岩和玄武岩类残坡积地表基质分布最广(图1)。
本文采取剖面研究和采样测试分析的方式,分析地表基质的空间结构特征和元素变化特征。空间结构主要调查地表基质的剖面结构、质地、岩石裂隙发育程度,质地在野外采用过筛法,微粒组分送相关实验室,采用激光粒度仪进行检测,依据《中国土壤系统分类土族和土系划分标准》进行土质分类(张甘霖,2013)。样品主要包括土质化学样品、质地样品和岩石样品,采样位置见图1,每个点位分别采集6~10个样品,共采集土质化学样品30件,质地样品14件,岩石样品3件,并对其N、P、K、Ca、Na等24种主量和微量元素的含量与变化特征进行测试分析(骆珊等,2021),具体测定方法见表1。
表 1 各类指标测定方法Table 1. Measurement methods for indicators检测项目 分析方法 检测仪器 Al2O3、CaO、K2O、MgO、MnO、Na2O、SiO2、TFe2O3、P2O5、TiO2 X射线荧光光谱法 X射线荧光
光谱仪FeO 基准重铬酸钾溶液滴定 — As、Cd、Co、Cu、
Mo、Pb、V 、Zn电感耦合等离子体质谱法 等离子质谱仪 Cl 粉末压片-X射线
荧光光谱法X射线荧光
光谱仪F 离子选择电极法 离子选择性电极 Hg 蒸汽发生-冷原子荧光
光谱法原子荧光光谱仪 N 凯氏法 — 2. 结果分析
2.1 空间结构异质性
地表基质的空间结构演化与岩石风化成土密切相关,同时母质岩石或者沉积物的性质在很大程度上决定着地表基质的理化性质(Li et al., 2016;裴小龙等,2020)。本文主要对研究区花岗岩、流纹岩和玄武岩3类岩石建造的地表基质剖面结构特征和质地特征进行了分析(图2)。
从剖面结构看,花岗岩类地表基质主要为残坡积砂土(PM01),垂向结构为雏形土−母质层,厚度30~50 cm,向下在坡度较缓的区域增厚,整体风化层较深,其浅部非构造节理为锯齿状非贯通型节理,节理密度中等,平均5~10 个/m2,张开度多小于5 mm,以紧闭型为主,风化的残积物及土壤等充填物较少。流纹岩类地表基质上坡区主要为残坡积砂土(PM02),30~50 cm,垂向结构为雏形土−母质层;在缓坡区为残坡积壤质砂土、壤土,土质层较厚,垂向结构为淋溶层−淀积层−母质层,土壤类型为红壤,风化层深度整体较浅,并发育平直型、锯齿状的贯通性节理,节理密度平均为2~3个/m2,张开度集中在1.0~5.0 mm和大于5 mm,属于中等张开程度,风化残积物多属于断续充填。玄武岩类地表基质(PM03)土壤发育程度较高,陡坡区垂向结构为残坡积壤质砂土,厚度30~100 cm,风化深度较厚;在缓坡为残坡积壤土,土质层较厚,垂向结构为淋溶层−淀积层−母质层,土壤类型为红壤,风化层深度整体较浅。玄武岩浅部发育层状节理,节理密度平均为10~15个/m2,深部发育柱状节理,多为贯通性节理,密度平均为3~5 个/m2。
从质地分析结果(小于2 mm粒度)看,花岗岩粒度分布较宽,在2~20 μm有一个主峰值,在0.2~1 μm和50~200 μm分别有较小的峰值,在野外的筛分法中,花岗岩地表基质表层粒度明显较粗,成分主体由难风化的石英颗粒构成。流纹岩粒度分布较稳定,整体粒度峰值为2~20 μm,表层与下部层位差异不明显。玄武岩表层与母质层主峰在20 μm附近,淀积层整体在0.2~1 μm附近,粒度偏细(图3)。
综上,3类典型母质岩石形成的地表基质空间结构特征差异较明显,表现出较强的空间异质性。其中,花岗岩母质形成的地表基质土质粗,易形成砂土,裂隙、孔隙发育,深度大,透水性及含水性强;流纹岩母质形成的地表基质土质较粗,易形成壤质砂土、壤土,风化深度浅,裂隙、孔隙相对发育,浅表透气好,水分充足;玄武岩母质形成的地表基质土质细,易形成壤土、粘质壤土,节理发育,孔隙不发育,通气性差。
2.2 元素特征异质性
山丘区地表基质的空间异质性主要受基岩结构与矿物成分影响(Clair et al., 2015;Xavier et al., 2015),同时岩石风化成土与元素的迁移过程直接影响地表基质的养分状况(Katemaher et al., 2019)。因此,探讨地表基质的元素异质性,需要从母质岩石的背景元素、岩石风化强度及元素迁移的特点进行分析,总结元素异质性特征及产生的原因。
2.2.1 地质背景元素
通过对比不同地质背景的元素变化,分析地表基质元素异质性特征。从主量元素特征看,花岗岩、流纹岩富Si、K,贫Fe、Mg,而基性岩恰相反,富Fe、Mg,贫Si、K,土质元素A、B层也表现为同样的富集特点,即岩石元素含量偏高时,土质中的含量相对偏高,Na、Ca受淋滤作用影响,A层含量普遍偏低(图4−a)。从微量元素特征看,玄武岩中微量元素含量普遍偏高,仅Mo含量相对较低;花岗岩、流纹岩类Se、Mo相对较高,其他元素含量相对较低,A、B、C层中微量元素具有相似的富集关系,非金属元素多数表现为表层含量偏高,这与生物作用密切相关(图4−b)。综合主量和微量元素特征看,地表基质自下而上具有明显的继承性。
本文依据《土地质量地球化学评价规范》(DZ/T 0295—2016),对剖面不同表层土壤的8种养分元素(Co、Ca、Cu、Zn、Mg、Mn、Mo、Fe)进行对比分析(图5)。从不同地质背景土壤养分元素看,花岗岩土壤中多数元素较缺乏或缺乏,仅Mo丰富;流纹岩土壤中Mo、Zn丰富,Fe、Mn中等—较缺乏,其他元素均缺乏;玄武岩土壤Fe、Mg、Co、Cu、Zn等元素均为丰富,Ca较缺乏,Mo缺乏。总体来说,玄武岩风化形成的土壤养分元素较丰富,而花岗岩、流纹岩形成的土壤养分元素整体较缺乏。
2.2.2 风化强度
研究表明,岩石的地球化学特征直接影响地表基质各层位的元素含量,目前用来描述衡量岩石或风化壳化学风化强度的指标,最常用的是化学蚀变指数CIA(包志伟等,1992;Priscia et al., 1999;刘成禹等,2011)。
CIA=[Al2O3/(Al2O3+CaO∗+Na2O+K2O)]×100 (1) 式中,氧化物按摩尔百分比计算,CaO*指存在于硅酸盐中的CaO。
3种典型岩石类型的化学风化强度具有一定差异(表2),玄武岩化学风化强度值明显高于花岗岩、流纹岩,意味着玄武岩更容易释放岩石中的营养元素,易于形成地表基质的养分。
表 2 地表基质风化强度CIA指数Table 2. CIA index of weathering strength of surface substrates点位 层位 CIA 点位 层位 CIA D106 淋溶层 73.09 D206 淋溶层 77.98 母质层 79.68 淀积层 83.36 D204 淋溶层 74.95 母质层 92.29 母质层 77.03 D302 淋溶层 85.86 D303 淋溶层 97.62 淀积层 87.74 母质层 95.96 母质层 98.26 2.2.3 元素迁移系数
本文以母岩元素含量为基础,分析元素的迁移特征,进而分析不同元素的迁移富集规律。迁移系数K由公式(2)计算获得(宁晓波等,2009;刘孜等,2020)。
K=Cn/CD (2) 式中,K为迁移系数,Cn为土壤各发生层元素含量,CD为母岩元素含量。其中,当K>1时表示富集,K<1时表示贫瘠。KA、KC、KR分别代表表层、母质层、基岩层的K值。
通过统计不同点位的元素迁移特征发现,地表基质中向上富集的元素主要有N、Cl、K及较多的微量元素,Na、Ca、Mg等元素流失量最大(图6)。As、Cd、Pb等重金属元素明显富集,N的K值最高,在10~30之间;Na和P元素表层微弱富集,母质层明显贫瘠。在玄武岩缓坡区个别点位,由于存在汇水区,表现为Na、Ca富集。
2.3 地表基质空间异质性对植被生态的影响
地表基质的空间结构和元素分异是岩石风化、水文和生物化学长期共同作用的结果(Sharma et al., 2000;Li et al., 2016;Russell et al., 2022)。岩石风化时,裂隙不断崩解形成微粒,地表基质的空间结构在横向上表现为粗骨土—壤土的变化,垂向上表现为粒度向上不断变细的特征。元素在不同地质作用下会发生不同的迁移富集行为,造成元素分异的特征。风化作用下,易风化矿物元素优先释放,如Fe、Mg、Na、Ca等;水文作用下,Na、K、Ca、Mg等活泼金属的离子态受流水作用而流失;生物化学作用下,N、Cl、F、K等元素逐渐向表层富集。总体看,生物所需元素和重金属元素会向表层富集(图7),活泼金属元素向下游迁移。
岩石母质的元素背景和岩石风化程度一定程度上控制着地表基质的养分状况,进而影响植被的群落构建和空间分布。地表基质的空间结构直接影响表层土的地气和水分交换、微生物群落类型和有机质组成,进而影响植被根系的生长和分布(黄金廷等,2008;赵建飞等,2022)。基岩的抗风化矿物的浓度调节其蓄水能力,进而直接影响植被的空间分布(Russell et al., 2022)。基岩中抗风化矿物浓度越高,颗粒越大,越影响地表基质中水分的分布;当基岩中抗风化矿物含量少时,易风化矿物将释放更多的养分,地表基质层中的养分元素含量也随之升高(汪振立等,2002;顾尚义等,2003;Roering et al., 2010)。从水分与养分的供给看,花岗岩风化形成的土壤,土质以砂土居多,水分蒸发快,水肥易流失,浅表缺水,深部含水,植被必须耐旱或者根系可以达到潜水/裂隙水才能生长,适宜松、杉、灌木等植被。流纹岩以粉土居多,下坡区表层砂壤土具有良好的通气性,下层壤土、粘壤土具有一定保水性,且风化层较浅,但部分养分元素相对缺乏,适宜槭、栎、竹、茶等多类植被。玄武岩整体保水性较好,矿质营养元素丰富,但土壤中粘粒含量高,不利于植被向下层扎根,且部分土壤通气性差,对植被根系具有要求,适宜杉、竹、茶等多类植被(表3)。
表 3 三种典型地表基质的植被适宜条件Table 3. Vegetation suitability conditions for three typical ground substrate地表基质
类型养分特征 水分制约 优势根系 适宜植被 花岗岩类 缺乏,富K、Mo 深层水分
较丰富垂直根系 松、杉、
灌木等流纹岩类 缺乏—中等,富K、Zn、P 上坡相对缺水,下
坡水分供给丰富垂直根系、辐射根系 槭、栎、竹、茶等多类
植被玄武岩类 多数元素丰富 水分充沛,
中层隔水辐射根系、串珠根系 杉、竹、茶等多类植被 3. 讨 论
在自然环境中,因受地形、降水、水文等因素的影响,地表基质结构特征和元素特征的空间分布往往不是均质的,表现出空间分布上的复杂性和不均匀性,具有明显的空间异质性特征。地表基质与其他资源要素的空间异质性具有一定的相似性,由时间和空间维度上的综合作用形成,具有空间梯度性和空间缀块性特征(Li et al., 1995;张志永,2017)。空间异质性可以理解成尺度的函数,具有尺度的表现性,随着尺度的改变,其空间异质性也会发生改变(Titus et al., 2002;王志述等,2015),意味着离开尺度去讨论其空间异质性将毫无意义。同时,地表基质空间异质性为促进各资源要素寻求和占据其适宜生态条件,形成多元空间格局提供了基础条件,这对优化景观空间配置、保护生物多样性和维持生态系统稳定具有重要意义(徐武美等,2015;刘绿怡,2018)。
地表基质空间异质性影响着植被生态,如植被群落的构建、组成、分布等。同时,植被生态对地表基质空间异质性具有一定的反馈作用,如植被的物种组成、营养积累和贮留能力会影响地表基质元素的空间分布格局,此外,植被的退化将加剧土壤的风蚀和营养元素的迁移,进而影响地表基质元素空间异质性(胡忠良等,2009;刘晨,2015)。因此,地表基质空间异质性对植被的影响是一种复杂的相互作用过程,本文仅从地质学角度定性探讨了地表基质对植被的影响机理,具有一定的局限性,今后还需加大研究力度,定量分析其影响规律和反馈机制,对于进一步掌握现有植被生态状态并预测其动态变化具有重要意义。
4. 结 论
本文以典型剖面阐述了3类地表基质的空间异质性,结合岩石的结构特征和元素组成,分析了地表基质空间异质性产生的原因,进而探讨了不同植被与地表基质的适宜性。
(1)宏观尺度下,地表基质的异质性受岩石风化、生物地球化学、水文等作用共同影响。岩石风化作用和水文作用影响基质的结构、厚度和粒度变化。受化学风化、水文作用及生物化学作用的影响,不同元素表现出不同的地球化学行为,相似的迁移行为促使性质相近的元素具有类似的分布规律。
(2)景观尺度下,岩石地质建造控制着景观尺度下的地表基质空间异质性。基岩中抗风化矿物浓度越高,颗粒越大,其形成的土质越粗;当基岩中某些元素含量越高,其形成的土质元素含量越高。
(3)地表基质主要影响水分和养分的分布,进而影响着植被类型。因地表基质的空间结构制约着水分空间分布,而元素的异质性表现为不同的养分状况,水分和养分的空间分布和丰度差异,制约着山丘区的植被生长。
致谢:感谢长三角宁波地区自然资源地表基质层调查项目组提供的基础资料。感谢审稿专家的宝贵意见和建议,极大地提升了本文的质量,感谢团队成员的辛勤工作和努力付出。
-
图 1 惠水土壤取样位置图(贵州省有色金属和核工业地质勘查局等,2019)
Figure 1. Soil sampling map of Huishui County
图 2 惠水地质图(贵州省有色金属和核工业地质勘查局等, 2019)
Figure 2. Geological sketch of Huishui County
图 3 惠水县岗度-龙塘山泥盆系—下石炭统对比图(贵州省地质调查院,2017)
Figure 3. The correlation map of Devonian to Lower Carboniferous of Gangdu-Longtangshan in Huishui
表 1 单项污染指数分级标准
Table 1 Pollution degree classification of single pollution index
指数 程度 影响水平 ≤0.7 安全 土壤、农作物清洁 0.7~1 警戒 土壤、农作物尚清洁 1~2 轻度污染 土壤、农作物受到轻度污染 2~3 中度污染 土壤、农作物受到中度污染 >3 重度污染 土壤、农作物受到重度污染 表 2 土壤环境质量农用地土壤污染风险管控标准(试行)(GB15618—2018)
Table 2 Soil pollution risk control standards for agricultural land of soil eavironmental quality(trial)(GB15618—2018)
mg/kg 项目 Cd Hg As Pb Cr Cu Ni Zn 水田 其他 水田 其他 水田 其他 水田 其他 水田 其他 果园 其他 风险筛选值 pH≤5.5 0.3 0.3 0.5 1.3 30 40 80 70 250 150 150 50 60 200 5.5 < pH≤6.5 0.4 0.3 0.5 1.8 30 40 100 90 250 150 150 50 70 200 6.5 < pH≤7.5 0.6 0.3 0.6 2.4 25 30 140 120 300 200 200 100 100 250 pH>7.5 0.8 0.6 1 3.4 20 25 240 170 350 250 200 100 190 300 表 3 惠水县耕地土壤元素及指标地球化学参数统计
Table 3 Statistical results of cultivated soil elements and index geochemical parameters in Huishui County
项目 As Cd Cr Cu Hg Ni Pb Zn pH 最大值 659.00 15.00 496.00 773.00 9.120 523.00 205.00 624.00 8.78 最小值 0.18 0.01 8.76 0.24 0.004 1.23 3.76 6.46 3.92 平均值 15.13 1.02 87.95 20.29 0.160 25.88 26.50 78.49 5.84 标准差 17.31 1.70 48.13 21.85 0.230 17.03 11.98 49.59 0.91 变异系数 1.14 1.67 0.55 1.08 1.440 0.66 0.45 0.63 0.16 中位数 12.15 0.82 104.80 29.90 0.120 31.35 27.95 86.10 6.66 贵州省土壤A层背景值 20.00 0.70 95.90 32.00 0.100 39.10 35.20 99.50 6.20 中国土壤A层背景值 11.70 0.10 61.00 22.60 0.100 26.90 26.00 74.20 6.70 k1 0.76 1.46 0.92 0.63 1.600 0.66 0.75 0.79 0.94 k2 1.29 10.20 1.44 0.90 1.600 0.96 1.02 1.06 0.87 注:pH为无量纲,重金属元素单位为mg/kg;k1=平均值/贵州省土壤A层背景值,k2=平均值/贵州省土壤A层背景值 表 4 不同成土母岩类型、不同成因类型、不同土壤类型中重金属含量平均值
Table 4 Average value of soil heavy metal content in different parent rock types, different genetic types and different soil types
mg/kg 元素 As Cd Cr Cu Hg Ni Pb Zn 样品数 碳酸盐岩 16.29 1.18 90.32 17.98 0.17 25.48 27.07 80.61 3760 石英砂岩 10.91 0.31 67.25 13.59 0.10 21.80 21.90 62.17 327 粉砂岩 12.57 0.63 85.86 28.51 0.10 28.86 26.19 75.40 178 泥岩 11.22 0.54 85.28 34.82 0.12 29.44 25.53 75.20 648 元素 As Cd Cr Cu Hg Ni Pb Zn 样品数 黄壤 12.29 0.37 71.02 17.84 0.13 21.65 22.07 61.43 1711 石灰土 20.20 2.06 115.13 20.41 0.19 31.17 32.28 103.74 1793 水稻土 12.27 0.49 73.55 22.07 0.16 23.75 24.62 66.51 1255 紫色土 11.09 0.43 77.22 31.62 0.09 28.80 23.52 71.85 154 元素 As Cd Cr Cu Hg Ni Pb Zn 样品数 水田 11.35 0.61 75.13 17.62 0.14 22.13 23.93 65.52 3000 旱地 21.35 1.75 109.88 24.18 0.19 32.20 31.04 100.63 1787 果园 16.53 0.48 83.75 30.62 0.12 26.82 22.99 73.66 106 茶园 19.60 0.66 73.57 18.15 0.14 18.34 23.50 69.87 20 表 5 土壤重金属单因子指数及分级
Table 5 Single pollution index and classification of soil
项目 重金属 As Cd Cr Cu Hg Ni Pb Zn Pi 最大值 21.97 25.00 2.48 7.73 15.20 5.23 1.71 2.50 Pi 最小值 0.01 0.02 0.04 0.00 0.01 0.01 0.03 0.03 Pi 平均值 0.50 1.70 0.44 0.20 0.27 0.26 0.22 0.31 单因子指数统计/% As Cd Cr Cu Hg Ni Pb Zn 分级 安全 79.83 52.70 87.69 97.87 96.91 98.32 99.82 95.47 警戒 11.39 12.79 8.38 1.56 1.56 1.44 0.14 2.95 轻度 7.16 14.62 3.91 0.54 0.88 0.22 0.04 1.56 中度污染 0.96 5.53 0.02 0.00 0.38 0.00 0.00 0.02 重度污染 0.66 14.36 0.00 0.04 0.28 0.02 0.00 0.00 表 6 土壤重金属地累积指数(Ig)统计结果
Table 6 Geo-accumulation index(Ig)of heavy metals in soil
项目 As Cd Cr Cu Hg Ni Pb Zn 最大值 4.46 3.92 1.78 4.01 5.79 3.16 1.97 2.06 最小值 -7.38 -6.63 -4.04 -7.64 -5.37 -5.57 -3.80 -4.54 平均值 -0.99 0.04 -0.71 -1.24 -0.05 -1.18 -0.99 -0.93 地累积指数占比/% As Cd Cr Cu Hg Ni Pb Zn 严重污染 0.00 0.00 0.00 0.00 0.02 0.00 0.00 0.00 重度污染 0.02 0.00 0.00 0.02 0.12 0.00 0.00 0.00 偏重度污染 0.00 1.40 0.00 0.02 0.40 0.02 0.00 0.00 中度污染 0.26 5.21 0.00 0.00 0.84 0.00 0.00 0.02 偏中度污染 1.34 7.12 0.54 0.70 3.77 0.06 0.04 0.60 轻度污染 0.70 21.15 10.89 5.25 25.11 3.83 3.81 6.54 无污染 98.88 77.50 88.57 94.02 69.74 96.09 96.15 92.84 表 7 土壤重金属相关系数
Table 7 Correlation coefficients between heavy metals in the soil
项目 As Cd Cr Cu Hg Ni Pb Zn pH As 1 Cd 0.312 1 Cr 0.388 0.630* 1 Cu 0.158 0.172 0.232 1 Hg 0.433* 0.242 0.188 0.117 1 Ni 0.338 0.545* 0.582* 0.670* 0.220 1 Pb 0.402* 0.669* 0.645* 0.298 0.296 0.605* 1 Zn 0.384 0.775* 0.644* 0.445* 0.315 0.781* 0.754* 1 pH 0.143 0.348 0.284 0.215 0.108 0.401* 0.373 0.436* 1 注:*代表P < 0.01 表 8 土壤重金属主成分分析结果
Table 8 Results of principal component analysis of soil heavy metals
项目 F1 F2 F3 As 0.53 -0.57 0.13 Cd 0.79 -0.09 -0.06 Cr 0.77 -0.11 -0.15 Cu 0.52 0.43 -0.20 Hg 0.41 -0.53 0.47 Ni 0.84 0.25 -0.18 Pb 0.84 -0.06 0.02 Zn 0.92 0.06 -0.06 pH 0.52 0.34 0.20 方差百分比% 44.77 12.14 10.25 累积贡献率% 44.77 56.92 67.16 注:F1—第一主成分;F2—第二主成分;F3—第三主成分 -
Acosta J A, Faz A, Martínez-marttínez S, et al. Enrichment of metals in soils subjected to different land uses in typical Mediterranean environment(Murcia City, southeast Spain)[J]. Applied Geochemistry, 2011, 26(3): 405-414. doi: 10.1016/j.apgeochem.2011.01.023
Bhuiyan M A H, Karmaker S C, Doza B M, et al. Enrichment, sources and ecological risk mapping of heavy metals in agricultural soils of dhaka district employing SOM, PMF and GIS methods[J]. Chemosphere, 2021, 263(128339): 1-14.
Cong J Y, Long H Y, Zhang Y, et al. Ecological environment response of benthic foraminifera to heavy metals and human engineering: A case study from Jiaozhou Bay, China[J]. China Geology, 2022, 5(1): 12-25.
Karim Z, Qureshi B A, Mumtaz M, et al. Heavy metal content in urban soils as an indicator of anthropogenic and natural influences on landscape of Karachi—A multivariate spatio-temporal analysis[J]. Ecological Indicators, 2014, 42: 20-31. doi: 10.1016/j.ecolind.2013.07.020
Li G, Sun G X, Ren Y, et al. Urban soil and human health, a review[J]. European Journal of Soil Science, 2018, 69: 196-215. doi: 10.1111/ejss.12518
Nanos N, Rodr ĺ guez Martin J A. Multiscale analysis of heavy metal contents in soils: Spatial variability in the Duero river basin(Spain)[J]. Geoderma, 2012, 189/190: 554-562. doi: 10.1016/j.geoderma.2012.06.006
Nemerow N L. Scientific stream pollution analysis[M]. New York: McGraw-Hill Companies, 1974.
Zou J Y, Song Z F, Cai K. Source apportionment of topsoil heavy metals and associated health and ecological risk assessments in a typical Hazy City of the North China Plain[J]. Sustainability, 2021, 13(10046): 1-14.
鲍丽然, 邓海, 贾中民, 等. 重庆秀山西北部农田土壤重金属生态健康风险评价[J]. 中国地质, 2020, 47(6): 1625-1636. https://www.cnki.com.cn/Article/CJFDTOTAL-DIZI202006003.htm 邓良基, 凌静, 张世熔, 等. 四川旱耕地生产、生态问题及水土流失综合治理研究[J]. 水土保持学报, 2002, 16(2): 8-11. https://www.cnki.com.cn/Article/CJFDTOTAL-TRQS200202002.htm 董秀茹, 刘浩洋, 刘洪彬. 基于耕地资源质量分类的辽宁省耕地土壤条件及空间分布特征分析[J]. 土壤通报, 2021, 52(5): 1020-1027. https://www.cnki.com.cn/Article/CJFDTOTAL-TRTB202105002.htm 杜梅, 张强英, 任培, 等. 西藏年楚河流域农用地土壤重金属分布与生态风险评价[J]. 环境工程技术学报, 2022, 12(5): 1618-1625. https://www.cnki.com.cn/Article/CJFDTOTAL-HKWZ202205032.htm 范云飞, 杨茗钛, 金来福, 等. 贵州省惠水县涟江大坝耕地质量现状及开发利用建议[J]. 贵州地质, 2020, 37(4): 490-496. https://www.cnki.com.cn/Article/CJFDTOTAL-GZDZ202004010.htm 冯志刚, 刘威, 张兰英, 等. 贫Cd碳酸盐岩发育的土壤Cd的富集与超常富集现象-以贵州岩溶区为例[J]. 地质通报, 2022, 41(4): 533-544. http://dzhtb.cgs.cn/gbc/ch/reader/view_abstract.aspx?file_no=20220402&flag=1 顾思博, 周金龙, 曾妍妍, 等. 新疆民丰县农田土壤重金属污染特征与生态风险评价[J]. 干旱区资源与环境, 2019, 33(12): 90-95. https://www.cnki.com.cn/Article/CJFDTOTAL-GHZH201912014.htm 贵州省地质调查院. 中国区域地质志: 贵州志[M]. 北京: 地质出版社, 2017: 1-592. 贵州省有色金属和核工业地质勘查局, 核资源地质调查院. 贵州省惠水县耕地质量地球化学调查评价报告[R]. 2019. 国家环境保护局. 中国土壤环境背景值[M]. 北京: 中国环境科学出版社, 1990: 330-369. 胡明. 大荔县农田土壤重金属分布特征与污染评价[J]. 干旱区资源与环境, 2014, 28(1): 79-84. https://www.cnki.com.cn/Article/CJFDTOTAL-GHZH201401014.htm 黄勇, 段续川, 袁国礼, 等. 北京市延庆区土壤重金属元素地球化学特征及其来源分析[J]. 现代地质, 2022, 36(2): 634-644. https://www.cnki.com.cn/Article/CJFDTOTAL-XDDZ202202024.htm 李鹏, 张惠娟, 徐莉, 等. 麦玉轮作区农田土壤重金属调查及评价[J]. 农业环境科学学报, 2022, 41(1): 46-54. https://www.cnki.com.cn/Article/CJFDTOTAL-NHBH202201006.htm 李伟, 高海涛, 张娜, 等. 拉萨市城区土壤重金属分布特征及生态风险评价[J]. 环境工程技术学报, 2022, 12(3): 869-877. https://www.cnki.com.cn/Article/CJFDTOTAL-HKWZ202203030.htm 梁立成, 余树全, 张超, 等. 浙江省永康市城区土壤重金属空间分布及潜在生态风险评价[J]. 浙江农林大学学报, 2017, 34(6): 972-982. https://www.cnki.com.cn/Article/CJFDTOTAL-ZJLX201706002.htm 廖启林, 华明, 金洋, 等. 江苏省土壤重金属分布特征与污染源初步研究[J]. 中国地质, 2009, 36(5): 1163-1174. https://www.cnki.com.cn/Article/CJFDTOTAL-DIZI200905022.htm 刘娟, 张乃明, 于鸿, 等. 沘江两岸耕地土壤重金属径流迁移模拟研究[J]. 安全与环境学报, 2021, 21(4): 1823-1831. https://www.cnki.com.cn/Article/CJFDTOTAL-AQHJ202104055.htm 骆珊, 张德明, 卢定彪, 等. 乌蒙山区毕节市耕地土壤微量元素丰缺评价及其影响因素[J]. 地质通报, 2021, 40(9): 1570-1583. http://dzhtb.cgs.cn/gbc/ch/reader/view_abstract.aspx?file_no=20210916&flag=1 吕柏楠, 王超, 师华定, 等. 基于受体模型与地统计的耕地土壤重金属污染源解析[J]. 环境科学研究. 2021, 34(12): 2962-2969. https://www.cnki.com.cn/Article/CJFDTOTAL-HJKX202112031.htm 生态环境部, 中华人民共和国国家市场监督管理总局. 土壤环境质量农用地土壤污染风险管控标准(试行)(GB 15618—2018)[S]. 2018. 田威, 李娜, 倪才英, 等. 江西省稻渔系统中土壤和稻谷重金属污染特征及健康风险评价[J]. 生态毒理学报, 2021, 16(3): 331-339. https://www.cnki.com.cn/Article/CJFDTOTAL-STDL202103030.htm 邬光海, 王晨昇, 陈鸿汉. 内蒙古废弃钨钼矿区周围土壤重金属污染生态环境评价及成因分析[J]. 中国地质, 2020, 47(6): 1838-1852. https://www.cnki.com.cn/Article/CJFDTOTAL-DIZI202006020.htm 息朝庄, 张鹏飞, 吴林锋, 等. 土壤重金属污染现状调查与评价: 以贵州省惠水县涟江高效农业园区为例[J]. 湖南城市学院学报(自然科学版), 2022, 31(4): 51-56. https://www.cnki.com.cn/Article/CJFDTOTAL-HNCG202204009.htm 息朝庄, 吴林锋, 张鹏飞, 等. 贵州省惠水土壤-灌溉水-雨水-大气降尘中Cd、As等微量元素特征及来源讨论[J]. 中国地质, 2023, 50(1): 192-205. https://www.cnki.com.cn/Article/CJFDTOTAL-DIZI202301015.htm 谢国雄, 应金耀, 章明奎. 大气沉降与施肥方式对梨园重金属平衡的影响[J]. 中国农学通报, 2019, 35(16): 88-94. https://www.cnki.com.cn/Article/CJFDTOTAL-ZNTB201916014.htm 叶嘉敏, 余厚平, 简敏菲, 等. 鄱阳湖流域农田重金属污染的生态风险评估[J]. 江西师范大学学报: 自然科学版, 2016, 40(4): 429-436. https://www.cnki.com.cn/Article/CJFDTOTAL-CAPE201604020.htm 叶霖, 刘铁庚. 贵州都匀牛角塘富镉锌矿床中镉的分布及福存状态探讨[J]. 矿物学报, 2001, 21(1): 115-118. https://www.cnki.com.cn/Article/CJFDTOTAL-KWXB200101018.htm 张连科, 李海鹏, 黄学敏, 等. 包头某铝厂周边土壤重金属的空间分布及来源解析[J]. 环境科学, 2016, 37(3): 1139-1146. https://www.cnki.com.cn/Article/CJFDTOTAL-HJKZ201603049.htm 中华人民共和国国土资源部. 多目标区域地球化学调查规范(1: 250 000)DZ/T 0258-2014[S]. 2014. 周墨, 唐志敏, 张明, 等. 赣州地区土壤-水稻系统重金属含量特征及健康风险评价[J]. 地质通报, 2021, 40(12): 2149-2158. http://dzhtb.cgs.cn/gbc/ch/reader/view_abstract.aspx?file_no=20211217&flag=1 周殷竹, 王彪, 刘义, 等. 青海囊谦县城周边农耕区土壤质量地球化学评价及富硒土地利用分区[J]. 干旱区资源与环境, 2020, 34(10): 93-101. https://www.cnki.com.cn/Article/CJFDTOTAL-GHZH202010014.htm 周永超, 孙慧兰, 陈学刚, 等. 绿洲城市伊宁市表层土壤重金属污染特征及其生态风险评价[J]. 干旱区资源与环境, 2019, 33(2): 127-133. https://www.cnki.com.cn/Article/CJFDTOTAL-GHZH201902019.htm