Distribution Characteristics and Source Analysis of Soil Heavy Metals in the Liaoyang—Dandong Region
-
摘要:
辽阳—丹东地区是辽宁省重要粮食与经济作物生产区,以往土壤重金属研究主要集中在200cm以浅,掌握土壤0~500cm重金属分布特征及来源,对监测黑土地质量及研究土壤重金属迁移转化意义重大。本文采集了1381件土壤表层(0~20cm)、160件土壤中层(150~250cm)和75件土壤深层(250~500cm)样品,对其中8种重金属元素和Sc等元素含量进行了测定,探讨研究区土壤重金属在不同层分布规律并解析来源。调查结果显示:土壤表层、中层和深层As、Cd、Cr、Hg、Co、Pb、Sc、Zn和Cu平均含量均未超过污染风险筛选值;通过多元统计,Cd和Hg表层富集明显,Pb和Zn表现一定富集趋势,Cd、Hg和Pb在表层分布不均匀;通过空间分析方法和地累积指数,研究区土壤表层中Cd、Hg、Pb和Zn整体表现为轻微污染,中度污染点零星分布在人员集中和采矿活跃区域;土壤中层Pb元素轻微污染,深层元素均未污染;利用Pearson相关性分析、主成分分析(PCA)和人为贡献率(ACR)揭示Pb、As、Co、Cu、Cr、Sc和Zn主要受自然背景的影响,Pb和Zn局部轻微受人为因素的影响,Cd和Hg则受人为活动的影响显著。研究区需要加强对Hg、Cd、Pb和Zn元素在不同地块富集趋势的监测,中层深度需要关注Pb的富集,同时开展区内不同流域土壤的Cd和Hg在不同层位的迁移转化研究。
要点(1)变异系数、富集因子等多元统计方法和地理空间分析,相互印证揭示土壤重金属元素分布特点。
(2)地累积指数评价土壤重金属污染程度,表层土壤Cd、Hg、Pb和Zn表现为轻微污染,中层Pb轻微污染,深层元素均未污染。
(3) Pearson相关性分析、主成分分析和人为贡献率解析土壤重金属来源,Cd和Hg受人为活动的影响显著。
HIGHLIGHTS(1) Multivariate statistical methods such as variation coefficient and enrichment factor and geospatial analysis confirm each other and reveal the distribution characteristics of heavy metal elements in soil.
(2) Soil heavy metal pollution was evaluated by index of geoaccumulation. The surface soil was slightly polluted by Cd, Hg, Pb and Zn, the middle layer of soil was slightly polluted by Pb, and the deep soil was not polluted.
(3) The sources of heavy metals in soil were analyzed by Pearson correlation analysis, principal component analysis and anthropogenic contribution rate, and Cd and Hg were significantly affected by anthropogenic activities.
Abstract:The Liaoyang—Dandong region is an important grain and cash crop production area in Liaoning Province. It is of great significance for monitoring the quality of black soil and studying the migration and transformation of heavy metals in soil to understand the distribution characteristics and sources of heavy metals in soil from 0 to 500cm. Soil samples were collected from different layers, and the contents of 8 heavy metal elements and Sc elements were determined. The average contents did not exceed the pollution risk screening value, Cd and Hg were obviously enriched in the surface soil, and Pb and Zn showed a certain enrichment trend. Cd, Hg, Pb and Zn in the surface soil were slightly polluted, the middle layer of soil was slightly polluted by Pb, and the deep layer of soil was not polluted. Pb, As, Co, Cu, Cr, Sc and Zn were affected mainly by natural background, Pb and Zn were affected slightly by local human factors, and Cd and Hg were affected significantly by human activities. In the future, it will be necessary to strengthen the monitoring of the enrichment trend of Hg, Cd, Pb and Zn elements in different plots. The BRIEF REPORT is available for this paper at http://www.ykcs.ac.cn/en/article/doi/10.15898/j.ykcs.202404070080.
BRIEF REPORTSignificance: The Liaoyang—Dandong area is an important grain producing area in Liaoning Province. Due to the intensive use of black land and the industrial expansion of the city, a heavy metal enrichment trend began to appear in this area, attracting professional attention. The degradation of geochemical properties of black soil often shows heavy metal toxic elements (As, Cd, Cr, Hg, Pb, Ni, etc.) exceeding the normal content in terms of trace elements. The purpose of this study was to explore the distribution characteristics of heavy metal elements in soil in the Liaoyang—Dandong area from a shallow perspective of 500cm, and analyze the sources of these heavy metal elements, to provide scientific basis for the ecological environment protection of black soil in this area.
Methods: 1381 surface (0−20cm), 160 middle (150−250cm) and 75 deep soil samples were collected from the Liaoyang—Dandong area (Fig.1). The contents of heavy metals As, Cd, Cr, Hg, Co, Pb, Zn, Cu and Sc in soil were determined by XRF or ICP-MS/OES or AFS method (Table 1). Multivariate statistical and spatial distribution analysis methods were used to reveal the distribution characteristics of heavy metals in soil, and the geoaccumulation index, which isolates anthropogenic pollution from the total concentration determined within a sample[23], was used to assess the pollution degree of different depths of soil. Pearson correlation analysis, principal component analysis (PCA) and anthropogenic contribution rate (ACR) were used to further explore the potential sources and effects of heavy metals.
Data and results: The variation coefficient of Cd and Hg in the surface soil was greater than or equal to 0.3, and the enrichment factor was greater than 1.9, showing the enrichment rule of the surface soil (Table 2 and Fig.2). The content of heavy metals in soil at different depths varied greatly. As, Cr, Co, Sc and Cu elements were almost identical in the peak area of the topsoil with those in the lower soil, which indicated that these elements had good inheritance characteristics. However, for Cd and Hg, their spatial distribution patterns were significantly different from those of other elements (Fig.3). Pearson correlation coefficient showed that Hg correlation was weak or not correlated in all regions, and Cd was not correlated in sedimentary rocks, moraines and quaternary alluvial regions (Table 3). Hg and Cd in surface soil were greatly affected by human activities. The evaluation results of geoaccumulation index showed that the surface soil was slightly polluted by Cd, Hg, Pb and Zn as a whole, the moderate pollution points were scattered in the areas where people were concentrated and mining was active, and the middle layer of soil was slightly polluted by Pb, which was due to the increase of ion exchange Pb content caused by topsoil acidification. According to the principal component analysis of heavy metal content in the surface and middle layers of soil and the calculation of anthropogenic contribution, As, Co, Cu, Cr, Pb and Zn originated mainly from the natural background, Pb and Zn were slightly affected by local human factors, while Cd and Hg were significantly affected by anthropogenic activities (Table 4, Fig.5). In addition, As was controlled mainly by sedimentary rocks, Pb was controlled mainly by granite, and Hg was partly attributed to metamorphic rocks.
-
土壤是人类活动赖以生存的基础,是农作物生产的主要载体,随工业现代化进程的加快,也影响着区域内物质和能量的重新分配[1]。土壤同时是生态系统的重要组成部分,是环境系统内物质、能量等转化的重要介质,成土母质决定了土壤的性质[2],此外土壤重金属元素还会受到人为活动不同程度的影响,因此研究土壤重金属元素分布特征,厘清重金属元素与成土母质之间的关系,在此基础上可以确定人为活动对重金属元素分布的影响程度[3]。表层土壤重金属包含自然背景的继承和人类活动的输入,一般情况下采用土壤深层重金属含量评价重金属污染背景值[4],成土母质和土壤采样一般限在200cm以浅[5],对于堆积覆盖较厚的地区,深层深度500cm以浅更能反映土壤母质的真实情况。
地学统计和多元统计分析方法是分析重金属元素空间分布特征的有效工具[6],主成分分析法可以区分人为因素还是自然因素对土壤重金属含量的影响[7-9]。研究者利用地累积指数(Igeo)与生态风险指数(IER)探讨了研究区的Ni、As、Cr、Cu、Pb、Cd及Zn的环境污染状况,同时还对其生态风险作出警示[10]。有学者运用土壤重金属静态环境容量的概念(QI)以衡量温州农田土壤的风险程度,并借助地理统计学的技术手段,明确了污染区域,同时也开展不确定性的分析,最后用正定矩阵因子分解模型(PMF)对可能存在的潜在源进行了数量化的分配[11],另有学者运用主成分分析法(PCA)和正定矩阵因子分解法(PMF)对沉积物重金属进行源解析,PCA 和 PMF 解析结果都表明鄱阳湖丰水期表层沉积物重金属受4 种来源的共同影响,其中矿业和工业活动的影响最大,相对贡献率为38%,其次是尾矿和废渣,相对贡献率为 28%,再是农业活动,相对贡献率为19%,最后是自然来源的相对贡献率为14%[12]。
辽阳—丹东地区是辽宁省重要的粮食产地,由于黑土地高强度利用与城市的工业扩张,该区域重金属开始出现富集趋势,部分区域重金属元素富集明显,引起了行业关注[13-14]。有学者选择辽河平原作为一个研究区,研究结果发现成土过程对Cd、Hg、As和Pb的平均贡献率分别为82.7%、85.2%、88.6%和91.7%,人为活动的平均贡献率分别为17.3%、14.8%、11.4%和8.3%,贡献率的空间分布表现为土壤环境背景与污染源的叠加[15]。本研究拟从500cm以浅的角度探究辽阳—丹东地区土壤重金属含量分布特征,并利用皮尔逊(Pearson)相关系数法、主成分分析(PCA)和人为贡献率(ACR)来解析这些重金属的来源,从而为该地区的黑土生态环境保护提供科学依据。
本次研究对象为土壤表层、中层和深层,深度范围为500cm以浅,表层为(0~20cm),中层为150~250cm,深层为250~500cm(含风化壳)。区域尺度上不同岩性地区土壤表层、中层和深层中重金属元素含量和空间分布的对比分析,也可以为区分重金属来源提供有价值的信息,依据1∶20万地质空间数据和土壤成因,将研究区划分为变质岩区、沉积岩区、第四系冰碛物区、第四系冲积物区、花岗岩区和闪长岩区等。钪(Sc)不是重金属,为了表述方便,将Sc合并在重金属范围内一并讨论。
1. 研究区地质概况
研究区位于辽宁省辽阳市和丹东市之间,地处中国东北地区,属于辽东半岛地质构造单元。该地区地质构造复杂,受华北板块、太平洋板块和松花江裂谷构造运动的影响,地形起伏较大,东部本溪—丹东一带以丘陵和低山地貌为主,海拔在0~1500,西部辽阳—鞍山一带地势相对平坦,区内有多条主要河流贯穿其中,部分形成了广阔的冲积平原,为当地农业生产提供了大量水源和肥沃土壤[16],人口密度适中,主要以农民居多,具有悠久的历史和丰富的民族文化[17]。地质构造和矿产资源丰富多样,东边岩浆岩是以花岗岩和闪长岩为主,沉积岩主要分布在研究区本溪和灯塔境内,变质岩主要分布在宽甸、本溪、海城和辽阳县境内,西部以第四系冲积物为主,土壤类型以棕壤、暗棕壤、草甸土、水稻土和褐土为主,为研究重金属元素的分布规律和源解析提供了基础[18]。
2. 实验部分
2.1 样品采集与分析测试
根据研究区土地利用、地质构造和土壤类型,采用1∶25万多目标区域地球化学调查规范采集样品[19-20](图1)。按照4km×4km网格采取土壤表层(0~20cm)样品1381件,按照12km×12km网格采集土壤中层(150~250cm)样品160件,在160个中层采样点中采集土壤深层(250~500cm)样品75件。三层样品每件质量约1.0kg,每8件表层样品中间采取1件中层样品。采集的土壤样品经过晾晒风干,碾碎并通过孔径为20目的尼龙筛进行处理,用聚乙烯袋进行储存,对采样点进行微调,尽量选择便于采集又符合布设的样点,同时标记好采样点、日期、土壤名称、编号等信息[21]。取一份150g左右样品,采用玛瑙无污染样品机器加工,将样品磨碎至200目,每件样品分成两份:一份纸袋装约25g供X射线荧光光谱法(XRF)测试;一份纸袋装约75g供ICP-MS、ICP-OES和AFS法测试[22]。
分析方法如下:①Cd采用电感耦合等离子体质谱法 (ICP-MS,X-Series Ⅱ型,美国ThermoFisher公司);②Sc、Co和Zn采用电感耦合等离子体发射光谱法 (ICP-OES,Icap7400型,美国ThermoFisher 公司);③As、Hg采用原子荧光光谱法(AFS,XGY-1011A型,中国地质科学院地球物理地球化学勘查研究所);④Cr、Cu和Pb采用X射线荧光光谱法(XRF,PW4400/40型,荷兰Panalytical公司)。
2.2 测试数据质量控制
土壤样品分析工作由自然资源部沈阳矿产资源监督检测中心完成,分析过程中严格执行多目标区域地球化学调查(1∶25 万)(DD2005-01)、《生态地球化学评价样品分析技术要求(试行)(DD2005-03)和《地质矿产实验室测试质量管理规范》(DZ 0130—2006)等技术标准。为保证分析测试质量,根据含量从国家一级标准物质GBW07401至GBW07412中每 50 件样品随机插入标准物质及重复样,各元素分析检出限、准确度、精密度(RSD)及报出率见表1。分析数据的准确度合格率达到98%以上,随机抽取5%样品的重复分析表明,相对偏差合格率达到98%以上,分析数据质量达到质量管理要求。
表 1 分析方法质量监控Table 1. Quality control of analysis methods元素 分析方法 检出限(mg/kg) 准确度(△logC) RSD(%) 报出率(%) As AFS 0.20 0.004~0.024 3.06~6.51 100 Cd ICP-MS 0.02 0.004~0.039 3.34~8.88 100 Cr XRF 1.80 0.002~0.011 0.56~3.71 100 Cu XRF 0.90 0.002~0.036 0.59~6.00 100 Hg AFS 0.0003 0.000~0.031 2.69~8.59 100 Pb XRF 1.00 0.000~0.032 0.44~4.90 100 Zn ICP-OES 0.30 0.000~0.016 0.30~8.58 100 Sc ICP-OES 0.30 0.000~0.035 2.50~6.26 100 Co ICP-OES 0.60 0.000~0.035 1.79~6.21 100 2.3 土壤重金属评价方法
2.3.1 地累积指数法
根据土壤重金属元素的含量和分布特征,可以采用土壤重金属地累积指数(Igeo)来评估土壤重金属污染程度。Igeo的计算公式如下所示[23]:
$$ I_{\mathrm{geo}}=\mathrm{log}_2\frac{C_{\mathrm{n}}}{k\times B_{\mathrm{n}}} $$ (1) 式中:Cn是重金属n的实际测量值(mg/kg);Bn则表示该种重金属元素在未受到人为污染时的背景数值。而k则是根据自然成岩过程中背景数值变化所选择的修正指标,k=1.5。
根据地累积指数,对污染严重程度的分类准则包括:当地累积指数Igeo≤0时,表示没有污染;若Igeo在0~1之间,被认为是轻度污染;若Igeo 在1~2之间,则视为中等污染;若Igeo在2~3之间,则是中等强度污染;若Igeo在3~4之间,表明有较强污染;若Igeo≥4,则可认定为强污染;若Igeo≥5,就意味着存在极强污染情况[24]。
2.3.2 人为贡献率
人为贡献率(ACR)是衡量土壤中重金属元素来源的重要参数之一,是Bing于2016提出的地球化学计算模型[25],通过对表层土壤样品进行分析,得出了人为贡献率的具体数值。研究结果显示,人为贡献率反映了土壤中重金属元素的人为活动来源占其总含量的比例,为进一步了解土壤表层中重金属元素的来源提供了重要依据。由于自然界中Sc元素分布广泛,且其分布特征极少受到外界因素的影响,人为贡献率(ACR)是以Sc元素作为参考元素,用来量化土壤中各类重金属元素的人为贡献程度。
计算公式[1]为:
$$ \mathrm{ACR}=100\times\frac{\mathrm{M\mathrm{e}}_{\mathrm{T}}-\mathrm{M}\mathrm{e}\mathrm{_N}}{\mathrm{M}\mathrm{e}_{\mathrm{T}}} $$ (2) $$ \mathrm{\mathrm{Me}_{\mathrm{N}}=\mathrm{Sc}\mathrm{_s}\times\left(\frac{\mathrm{M}\mathrm{e}}{\mathrm{S}\mathrm{c}}\right)_{\mathrm{B}}} $$ (3) 式中:ACR为人为贡献率(%),一般认为ACR值小于40%,则其人为贡献率较低;MeT为土壤重金属元素的实测含量(mg/kg);MeN为土壤重金属元素背景值的修正含量(mg/kg),代表自然成因;Scs为实际测量的Sc元素的含量(mg/kg);$\left( \dfrac{\mathrm{M}\mathrm{e}}{\mathrm{S}\mathrm{c}}\right)_{\mathrm{B}} $为重金属元素和Sc元素背景水平之间的比值。
2.3.3 富集因子
Sc元素因其强大的抗风化性能,将Sc作为比较元素对土壤表层、中层、深层元素含量进行标准化处理[26]。富集因子(EF)计算公式如下所示[27]:
$$ EF=\frac{\left(\dfrac{\mathrm{M}\mathrm{e}_{\mathrm{T}}}{\mathrm{S}\mathrm{c}_{\mathrm{s}}}\right)_{\mathrm{表}\mathrm{层}\mathrm{土}\mathrm{壤}}}{\left(\dfrac{\mathrm{M}\mathrm{e}_{\mathrm{T}}}{\mathrm{S}\mathrm{c}_{\mathrm{s}}}\right)_{\mathrm{深}\mathrm{层}\mathrm{土}\mathrm{壤}}} $$ (4) 式中:MeT为土壤重金属元素的实测含量(mg/kg);Scs为实际测量的Sc元素的含量(mg/kg)。
根据研究区重金属的分布特征,将富集因子分为三个等级:EF<2(无或轻微潜在污染);2<EF<5(中度潜在污染);EF>5(高度潜在污染)。
2.3.4 相关性分析
相关性分析是指对两个或多个具备相关性的变量元素进行分析,从而衡量两个变量因素的相关密切程度,相关性的元素之间需要存在一定的联系或者概率才可以进行相关性分析。
2.3.5 主成分分析
主成分分析(PCA)是一种多元数据的统计学研究方法,它通过对数据矩阵进行线性转换,将初始变量转化为一系列无关联的主要组成部分,以揭示数据之间的内在结构和规律。
2.4 数据处理
数据处理过程中,本课题组首先对获得的原始数据进行筛选,采用Grubbs法去除了可能存在的异常值和错误数据,以确保数据的可靠性。然后,采用地学统计和多元统计方法对筛选后的数据进行计算,包括计算每种重金属的平均值、标准差和变异系数等统计指标,揭示其含量分布特征[28]。
同时,使用GIS和SPSS软件分析了不同层位土壤重金属元素之间的相关性,通过空间分布分析,进一步揭示重金属元素在土壤中的空间分布特征和相互作用规律[29]。为了深入了解重金属元素的来源,本文运用主成分分析(PCA)、人为贡献率(ACR)对可能的污染源进行辨识和解析。
3. 结果与讨论
3.1 土壤重金属元素含量特征
通过对比土壤表层、中层和深层重金属元素的平均含量(表2,图2),发现As、Cd、Hg、Co、Pb、Cu和Zn在表层中的含量高于中层和深层,而Cr则相反;表层Cd、Hg和Pb含量的中值小于平均值,中层和深层As、Cd、Cr、Hg、Co、Pb和Cu含量的中值小于平均值,反映出研究区以上重金属含量明显存在极大值;与辽河流域的土壤地球化学背景值相比,研究区表层土壤重金属含量普遍超过了辽河流域的背景值,但相较农用地土壤污染风险管控标准,均未超过污染风险筛选值。
表 2 土壤重金属元素含量统计Table 2. Statistics of heavy metal element contents in soil统计参数 样品数量 采样深度 As
(mg/kg)Cd
(mg/kg)Cr
(mg/kg)Hg
(mg/kg)Co
(mg/kg)Pb
(mg/kg)Sc
(mg/kg)Zn
(mg/kg)Cu
(mg/kg)平均值 1381 表层 9.643 0.243 70.311 0.063 14.991 39.856 11.695 85.196 24.827 160 中层 8.828 0.084 70.848 0.028 14.186 27.627 11.736 75.194 23.509 75 深层 9.442 0.107 70.709 0.033 14.081 30.484 11.752 76.426 24.246 中值 1381 表层 9.650 0.183 71.678 0.060 15.009 36.867 11.750 87.660 24.917 160 中层 8.435 0.078 70.700 0.025 14.000 27.300 11.800 75.950 22.100 75 深层 8.700 0.088 70.30 0.028 13.80 28.70 11.80 75.70 23.00 最小值 1381 表层 7.294 0.001 43.868 0.020 7.570 18.804 8.483 41.681 15.422 160 中层 1.410 0.021 19.700 0.006 3.250 2.260 4.350 21.700 6.640 75 深层 1.150 0.021 10.70 0.0047 2.260 2.260 3.090 12.00 3.54 最大值 1381 表层 11.79 1.093 142.353 0.181 20.673 99.575 15.357 140.552 32.513 160 中层 32.70 0.300 277 0.061 37.800 58.500 21.0 0 270 50.0 75 深层 49.50 1.480 516 0.250 51.90 231 31.40 402 115 标准偏差 1381 表层 0.984 0.150 13.144 0.019 2.505 14.862 1.483 14.606 3.432 160 中层 4.514 0.039 31.644 0.014 5.234 8.173 3.336 26.258 8.277 75 深层 5.470 0.092 33.293 0.024 5.279 14.082 3.381 27.646 10.446 变异系数 1381 表层 0.102 0.618 0.187 0.300 0.167 0.373 0.127 0.171 0.138 160 中层 0.511 0.464 0.447 0.517 0.369 0.296 0.284 0.349 0.352 75 深层 0.579 0.864 0.471 0.731 0.375 0.462 0.288 0.362 0.431 富集因子 − 表层/中层 1.096 2.903 0.996 2.249 1.060 1.448 1.000 1.137 1.060 − 表层/深层 1.026 2.282 0.999 1.903 1.070 1.314 1.000 1.120 1.029 背景值 − 表层 6.720 0.130 60.000 0.030 11.00 24.000 9.300 54.00 18.700 表层与中层、表层与深层比值富集因子(表2,图2)变化规律基本一致,而Cd、Hg、Pb和Zn两种富集因子均超过了1,说明这4种元素在研究区存在富集趋势,Cd和Hg两种元素的富集因子超过了1.9,表现明显富集,其余重金属元素的两种富集因子均在1左右。表层Cd和Pb变异系数大于中层,Cd变异系数为0.618,属于强度变异,表明局部地区污染源较为明显,Pb和Hg变异系数大于0.25,属于中等强度变异,人为因素对Pb和Hg的富集有一定影响,深层变异系数均大于表层和中层,可能是深层重金属元素含量受风化层影响。
3.2 土壤重金属元素空间分布特征
土壤元素地球化学图能够为识别元素潜在来源、圈定潜在污染区域提供有效的信息[30]。本文利用ArcGis 软件,通过空间克里金插值法(Kriging)绘制了研究区土壤表层和中层重金属元素地球化学图(图3),鉴于表层与中层、表层与深层的富集系数规律一致,本次空间数据只制作表层和中层重金属化学元素图。
图3展示了不同深度土壤层重金属含量空间差异性很大,As、Cr、Co、Sc和Cu元素在表层土壤的高峰区域几乎与其下部土壤一致,这表明这些元素具有良好的传承特性。然而,对于Cd和Hg而言,其空间分布模式显著异于其他元素。这两种元素在表层土壤的高峰地区也明显区别于中层土壤,由此可以推断出表层土壤内的Cd和Hg含量分布受到人类活动的潜在影响。从图3可以看出,表层Cd和Hg的高值区位于辽宁的海城市和辽阳地区,该地区一直以来是农业和采矿业比较活跃的地区,容易引起重金属在表层富集。
通过对变质岩、沉积岩、第四系冰碛物、第四系冲积物和花岗岩中层样品和其相对应的表层样品的数据进行皮尔逊相关性的研究,可以评估出这些样品间的关联度高低。研究区土壤表层和中层重金属元素之间的相关系数列于表3,显示Cr、Sc、Co元素之间存在着明显的正相关关系(相关系数介于0.521~0.891之间,p<0.01),冰碛物的As呈现弱相关,其他都是中等强度相关,Pb呈现中等程度相关,Zn除了变质岩区域,Cu除了冰碛物区域不相关,其他区域表现相关,说明表层土壤中Cr、Sc和Co主要来源于中层土壤,受母质控制,外源输入较小,而在所有类型区域,表层与中层土壤中的Hg呈现不相关或相关性较弱,在沉积岩、冰碛物、第四系冲积物区域,表层与中层土壤中的Cd表现不相关,表层土壤Hg和Cd受人为活动影响较大。
表 3 土壤表层和中层重金属元素之间的相关系数Table 3. Correlation coefficients between heavy metal elements in surface and deep soil元素 Pearson系数 变质岩 沉积岩 第四系冰碛物 第四系冲积物 花岗岩 As 0.491** 0.408* 0.285* 0.530** 0.491** Cd 0.658** 0.112 0.045 0.208 0.613** Cr 0.713** 0.641** 0.742 ** 0.571** 0.596** Hg 0.283 0.027 0.073 0.215 0.236 Co 0.723** 0.774** 0.723** 0.654** 0.600** Pb 0.584** 0.477** 0.481** 0.539** 0.507** Sc 0.620** 0.891** 0.558** 0.521** 0.541** Zn 0.078 0.623** 0.658** 0.391** 0.326* Cu 0.393* 0.899** −0.049 0.506** 0.290* 注:“**”表示在 0.01 级别(双尾)相关性显著;“*”表示在 0.05 级别(双尾)相关性显著。 3.3 土壤不同层位重金属污染情况(地累积指数)
经过计算,辽阳—丹东地区土壤重金属元素的地累积指数(Igeo)统计如图4所示。表层土壤中的As、Cr、Co、Sc、Cu的地累积指数均值都低于0,此外重金属在土壤表层、中层和深层地累积指数按顺序递减,土壤表层为:Hg(0.41)>Pb(0.06)>Zn(0.05)>Cd(0.02)>As(−0.07)>Co(−0.16)>Cu(−0.19)>Sc(−0.27)>Cr(−0.38);土壤中层为:Pb(0.71)>Zn(−0.19)>Sc(−0.31)>Co(−0.32)>Cu(−0.34)>As(−0.38)>Cr(−0.46)>Hg(−0.92)>Cd(−1.37);土壤深层为:Cd(−0.31>As(−0.35)>Zn(−0.36)>Cu(−0.39)> Sc(−0.45)>Pb(0.48)>Co(−0.51)>Hg(−1.18)>Cr(−1.36)。
土壤表层地累积指数除了Pb,其余元素都大于中层,土壤中层地累积指数除了Pb,其余元素均小于0,表层、中层Pb地累积指数分别为0.06和0.71,均为轻度污染。臧传子等[31]研究了不同pH值、有机质含量和土壤黏粒对Pb地球化学行为的影响,认为Pb离子交换态含量随pH值降低而升高,随有机质含量降低而升高,随黏粒含量降低而升高。由于长期的不合理耕作使得东北黑土退化严重,土壤有机质不断减少、土壤持续酸化等问题凸显[32]。研究区表层土壤酸化使Pb离子交换态比例上升,直接导致土壤中离子交换态Pb含量增加,在淋溶作用下活性离子交换态Pb在土壤中层富集;深层土壤重金属的地累积指数都小于0,表明深层没有受到重金属污染,在土壤中层富集的Pb元素在土壤深层没有富集;Igeo分布图显示,土壤表层Cd和Hg有超过60%的采样点地累积指数大于0,说明它们受到了一定程度的污染。对于森林、草甸、灌丛和草原表层土壤,大气零价汞沉降输入的贡献分别为:62.2%±12.5%、45.1%±11.7%、50.9%±9.7%、50.4%±12.7%,大气二价汞沉降输入的贡献分别为:10.2%±4.2%、18.3%±6.5%、14.0%±5.7%、16.5%±7.6%,地质来源汞的贡献分别为:27.7%±8.7%、36.6%±6.8%、35.2%±4.9%、33.1%±10.0%[33]。根据土壤表层、中层和深层Hg地累积指数,从表层到深层依次降低,结合Hg表层空间分布特点,高值区主要集中在鞍山一辽阳一本溪一带,结合地质构造,这一带是辽宁铁矿分布区,也是人类活动密集和森林覆盖区,研究区土壤表层的Hg富集,大部分来至大量Hg废渣经大气沉降输入;Cd重金属分布为面污染形态,研究区有多条河流经过,且属于农业高度利用区域,河流灌溉是土壤表层Cd富集的主要因素。
3.4 土壤重金属来源解析
3.4.1 主成分分析
本次研究中,主成分分析(PCA)用于揭示土壤中重金属元素的来源和污染情况。经过PCA分析发现表层第一主要成分与As、Zn和Cu存在正相关关系,而与Cd呈现负相关关系,第二主要成分与Cr、Pb、Co、Sc和Hg呈负相关关系。
在稳健多元成分数据分析图中,表层土壤前两个主轴解释了61%总累积方差(图5)。此两种主成分代表了土壤重金属的主要影响因素,不同区域取样的元素成分差异比较大。
第一主成分(PC1)对总体变异性的贡献率达到了44.3%,这表明其是决定土壤表层中重金属含量的关键要素。第二个主成分(PC2)对总体变异性的贡献率达到16.7%。
其中,As成分变异性最高,沉积岩区土壤表层的As含量最高,同时该岩区的土壤表层样品也集中在As元素周围,可以看出土壤表层As分布受沉积岩控制。
Co、Cu、Cr、Pb和Zn元素组合在第一主成分上具有显著的影响力,花岗岩区的取样点主要位于该元素组合周围,故这些元素组合在土壤表层中的含量分布特征主要是成土母质所决定的[34],也就是受花岗岩控制。对于Pb而言,由于其在土壤中层的含量高于表层,花岗岩对其控制更为显著,此外,随着研究区土壤表层酸化,土壤表层受人为影响富集的Pb在淋溶作用下发生了迁移,在土壤中层发生了富集。Cd和Hg元素在第二主成分上具有较高的权重。根据这些重金属的空间分布特征,可以看出这些重金属在人类活动的频繁地区及工业矿区的含量相对较高,因此有理由推断出,人为干扰可能是导致该地区土壤表层的Cd和Hg含量显著升高的关键因素。汪春鹏等[13]研究了辽阳地区的土壤重金属分布特点,表明Hg和Cd是最重要的生态环境威胁因素之一,Hg被列为较强的生态风险,Cd属于中等生态风险,这与本次研究的结果一致。此外,Hg元素在第一主成分上也具有较高的权重,变质岩的样品点主要分布Hg元素,可见Hg元素受到变质岩的控制,土壤表层Hg有来自土壤中层的迁移转化。
3.4.2 人为贡献率
人为活动对表层土壤中的重金属含量的影响程度差异较大,由于Sc是参比元素,本次不参与讨论。然而,值得注意的是,如表4所示,在所有采集的样品中,无论是Co还是Cu,其人为贡献率(ACR)皆小于40%,只有极少数点位的As、Cr及Zn的ACR值超出40%,并且这些数值都在60%以下,这意味着表层土壤中As、Co、Cu、Cr和Zn五种元素主要为自然风化成因来源,另外还有19%的Pb和1.3%的Zn的样品点ACR值超过40%,表明Pb和Zn轻微受到人为因素的影响;与此相反,Cd、Hg则受到了较为强烈的人为干扰,大约50%样品中Hg的ACR值高于40%,近30%样品中Cd的ACR值超过40%,近15%样品中Cd的ACR值超过60%。
表 4 表层土壤人为贡献率统计Table 4. Statistics of anthropogenic contribution rate of surface soil人为贡献率
数值范围样品数量(件) As Cd Cr Hg Co Pb Sc Zn Cu ACR>40% 2 398 5 686 0 262 − 18 0 ACR>60% 0 207 0 56 0 43 − 0 0 3.5 本文研究采集土壤的主要特点
本文研究的土壤为500cm以浅深度范围,与以往土壤的调查深度稍微有点不同。本文研究的250~500cm深度范围是以往土壤研究没有涉及的深度,这个深度对应到地表基质层次模型的沉积层[35],主要由风化层组成,相比于土壤中层,本次深层土壤重金属分布特征表现为:深层的8种重金属含量中值均小于表层中值,反映出研究区深层8种重金属元素含量明显存在极大值;表层与深层富集因子相对于表层与中层富集因子,变化规律基本一致,Cd和Hg富集因子略有降低,这与深层采集样品有关,深层包含风化层,风化层重金属的含量对深层土壤具有一定的影响;相比于中层(Pb的地累积指数大于0),深层土壤重金属地累积指数全部小于0,深层土壤均未受到人为因素的影响,因此本文认为表层与深层的比值更适合进行富集因子评价。
4. 结论
研究区土壤表层和中层各重金属含量变异系数和富集因子揭示的重金属分布特征为:土壤表层Cd和Hg变异系数大于0.25,富集因子大于1.9,呈现出土壤表层富集规律;地理空间分析显示As、Cr、Co、Sc和Cu元素在土壤表层的高峰区域几乎与其下部土壤一致,Cd和Hg的空间分布模式显著异于其他元素;地累积指数法评价研究区重金属污染情况,结果显示表层Cd、Hg、Pb和Zn整体处于轻微污染状态,中度污染点零星分布在人员集中和采矿活跃区域,土壤中层Pb表现轻微污染,原因是由于土壤表层酸化引起离子交换态Pb含量增加,在淋溶作用下迁移富集,土壤深层重金属均未受到污染;Pearson相关系数显示,所有区域Hg相关性较弱或不相关,在沉积岩、冰碛物、第四系冲积物区域,Cd表现不相关。
土壤表层Hg和Cd受人为活动影响较大;通过土壤表层重金属含量主成分分析和人为贡献率计算,As、Co、Cu、Cr、Pb和Zn主要来自自然背景的贡献,As主要受沉积岩控制,Pb主要受花岗岩控制,Pb和Zn局部轻微受到人为因素的影响,Cd和Hg受人为活动的影响显著,Hg部分来自变质岩的贡献。土壤中层和深层采样点分层深度不完全一致,样品采集数量不足是本文研究的主要弊端。
致谢:在论文写作过程中,中国地质调查局地球物理调查中心耿国帅高级工程师、于小健工程师、张沥元工程师、王乔助理工程师、杨硕助理工程师在制图和数据方面给予了大量帮助,在此一并表示感谢。
-
表 1 分析方法质量监控
Table 1 Quality control of analysis methods
元素 分析方法 检出限(mg/kg) 准确度(△logC) RSD(%) 报出率(%) As AFS 0.20 0.004~0.024 3.06~6.51 100 Cd ICP-MS 0.02 0.004~0.039 3.34~8.88 100 Cr XRF 1.80 0.002~0.011 0.56~3.71 100 Cu XRF 0.90 0.002~0.036 0.59~6.00 100 Hg AFS 0.0003 0.000~0.031 2.69~8.59 100 Pb XRF 1.00 0.000~0.032 0.44~4.90 100 Zn ICP-OES 0.30 0.000~0.016 0.30~8.58 100 Sc ICP-OES 0.30 0.000~0.035 2.50~6.26 100 Co ICP-OES 0.60 0.000~0.035 1.79~6.21 100 表 2 土壤重金属元素含量统计
Table 2 Statistics of heavy metal element contents in soil
统计参数 样品数量 采样深度 As
(mg/kg)Cd
(mg/kg)Cr
(mg/kg)Hg
(mg/kg)Co
(mg/kg)Pb
(mg/kg)Sc
(mg/kg)Zn
(mg/kg)Cu
(mg/kg)平均值 1381 表层 9.643 0.243 70.311 0.063 14.991 39.856 11.695 85.196 24.827 160 中层 8.828 0.084 70.848 0.028 14.186 27.627 11.736 75.194 23.509 75 深层 9.442 0.107 70.709 0.033 14.081 30.484 11.752 76.426 24.246 中值 1381 表层 9.650 0.183 71.678 0.060 15.009 36.867 11.750 87.660 24.917 160 中层 8.435 0.078 70.700 0.025 14.000 27.300 11.800 75.950 22.100 75 深层 8.700 0.088 70.30 0.028 13.80 28.70 11.80 75.70 23.00 最小值 1381 表层 7.294 0.001 43.868 0.020 7.570 18.804 8.483 41.681 15.422 160 中层 1.410 0.021 19.700 0.006 3.250 2.260 4.350 21.700 6.640 75 深层 1.150 0.021 10.70 0.0047 2.260 2.260 3.090 12.00 3.54 最大值 1381 表层 11.79 1.093 142.353 0.181 20.673 99.575 15.357 140.552 32.513 160 中层 32.70 0.300 277 0.061 37.800 58.500 21.0 0 270 50.0 75 深层 49.50 1.480 516 0.250 51.90 231 31.40 402 115 标准偏差 1381 表层 0.984 0.150 13.144 0.019 2.505 14.862 1.483 14.606 3.432 160 中层 4.514 0.039 31.644 0.014 5.234 8.173 3.336 26.258 8.277 75 深层 5.470 0.092 33.293 0.024 5.279 14.082 3.381 27.646 10.446 变异系数 1381 表层 0.102 0.618 0.187 0.300 0.167 0.373 0.127 0.171 0.138 160 中层 0.511 0.464 0.447 0.517 0.369 0.296 0.284 0.349 0.352 75 深层 0.579 0.864 0.471 0.731 0.375 0.462 0.288 0.362 0.431 富集因子 − 表层/中层 1.096 2.903 0.996 2.249 1.060 1.448 1.000 1.137 1.060 − 表层/深层 1.026 2.282 0.999 1.903 1.070 1.314 1.000 1.120 1.029 背景值 − 表层 6.720 0.130 60.000 0.030 11.00 24.000 9.300 54.00 18.700 表 3 土壤表层和中层重金属元素之间的相关系数
Table 3 Correlation coefficients between heavy metal elements in surface and deep soil
元素 Pearson系数 变质岩 沉积岩 第四系冰碛物 第四系冲积物 花岗岩 As 0.491** 0.408* 0.285* 0.530** 0.491** Cd 0.658** 0.112 0.045 0.208 0.613** Cr 0.713** 0.641** 0.742 ** 0.571** 0.596** Hg 0.283 0.027 0.073 0.215 0.236 Co 0.723** 0.774** 0.723** 0.654** 0.600** Pb 0.584** 0.477** 0.481** 0.539** 0.507** Sc 0.620** 0.891** 0.558** 0.521** 0.541** Zn 0.078 0.623** 0.658** 0.391** 0.326* Cu 0.393* 0.899** −0.049 0.506** 0.290* 注:“**”表示在 0.01 级别(双尾)相关性显著;“*”表示在 0.05 级别(双尾)相关性显著。 表 4 表层土壤人为贡献率统计
Table 4 Statistics of anthropogenic contribution rate of surface soil
人为贡献率
数值范围样品数量(件) As Cd Cr Hg Co Pb Sc Zn Cu ACR>40% 2 398 5 686 0 262 − 18 0 ACR>60% 0 207 0 56 0 43 − 0 0 -
[1] 黄勇, 段续川, 袁国礼. 北京市延庆区土壤重金属元素地球化学特征及其来源分析[J]. 现代地质, 2022, 36(2): 634−644. doi: 10.19657/j.geoscience.1000-8527.2022.02.24 Huang Y, Duan X C, Yuan G L. Geochemical characteristics and sources of heavy metals in soils of Yanqing District, Beijing[J]. Geoscience, 2022, 36(2): 634−644. doi: 10.19657/j.geoscience.1000-8527.2022.02.24
[2] 孙凯, 孙彬彬, 周国华, 等. 福建龙海土壤重金属含量特征及影响因素研究[J]. 现代地质, 2018, 32(6): 197−205. doi: 10.19657/j.geoscience.1000-8527.2018.06.18 Sun K, Sun B B, Zhou G H, et al. Characteristics and influencing factors of heavy metals in soils of Longhai, Fujian[J]. Geoscience, 2018, 32(6): 197−205. doi: 10.19657/j.geoscience.1000-8527.2018.06.18
[3] 汤金来, 赵宽, 胡睿鑫, 等. 滁州市表层土壤重金属含量特征、源解析及污染评价[J]. 环境科学, 2023, 44(6): 3562−3572. Tang J L, Zhao K, Hu R X, et al. Characteristics, source analysis and pollution assessment of heavy metals in surface soil of Chuzhou[J]. Environmental Science, 2023, 44(6): 3562−3572.
[4] 宋运红, 杨凤超, 刘凯, 等. 三江平原耕地土壤重金属元素分布特征及影响因素的多元统计分析[J]. 物探与化探, 2022, 46(5): 1064−1075. doi: 10.11720/wtyht.2022.0048 Song Y H, Yang F C, Liu K, et al. Multivariate statistical analysis of distribution characteristics and influencing factors of heavy metal elements in cultivated soil in Sanjiang Plain[J]. Geophysical and Geochemical Exploration, 2022, 46(5): 1064−1075. doi: 10.11720/wtyht.2022.0048
[5] 姚晓峰, 杨建锋, 左力艳, 等. 地表基质的内涵辨析与调查思路[J]. 地质通报, 2022, 41(12): 2097−2105. doi: 10.12097/j.issn.1671-2552.2022.12.002 Yao X F, Yang J F, Zuo L Y, et al. Connotation analysis and investigation of surface matrix[J]. Geological Bulletin of China, 2022, 41(12): 2097−2105. doi: 10.12097/j.issn.1671-2552.2022.12.002
[6] Yuan G L, Sun T H, Han P, et al. Environmental geochemical mapping and multivariate geostatistical analysis of heavy metals in topsoils of a closed steel smelter: Capital iron & steel factory, Beijing, China[J]. Journal of Geochemical Exploration, 2013, 130(1): 15−21. doi: 10.1016/j.gexplo.2013.02.010
[7] Wang A T, Wang Q, Li J, et al. Geo-statistical and multivariate analyses of potentially toxic elements’ distribution in the soil of Hainan Island (China): A comparison between the topsoil and subsoil at a regional scale[J]. Journal of Geochemical Exploration, 2019, 197: 48−59. doi: 10.1016/j.gexplo.2018.11.008
[8] 王诚煜, 李玉超, 于成广, 等. 葫芦岛东北部土壤重金属分布特征及来源解析[J]. 中国环境科学, 2021, 41(11): 5227−5236. doi: 10.19674/j.cnki.issn1000-6923.20210608.007 Wang C Y, Li Y C, Yu C G, et al. Distribution characteristics and sources of heavy metals in soil of Northeast Huludao[J]. China Environmental Science, 2021, 41(11): 5227−5236. doi: 10.19674/j.cnki.issn1000-6923.20210608.007
[9] 王建明, 施泽明, 郑培佳, 等. 四川铅锌冶炼工业区周边土壤重金属地球化学特征及源解析[J]. 地球与环境, 2023, 51(3): 287−298. doi: 10.14050/j.cnki.1672-9250.2022.050.083 Wang J M, Shi Z M, Zheng P J, et al. Geochemical characteristics and source apportionment of heavy metals in soil around Sichuan lead-zinc smelting industrial zone[J]. Earth and Environment, 2023, 51(3): 287−298. doi: 10.14050/j.cnki.1672-9250.2022.050.083
[10] Zhu Y, An Y F, Li X Y, et al. Geochemical characteristics and health risks of heavy metals in agricultural soils and crops from a coal mining area in Anhui Province, China[J]. Environmental Research, 2023, 241(1): 117670−117680. doi: 10.1016/j.envres,2023.117670
[11] Xia F, Zhao Z F, Niu X, et al. Integrated pollution analysis, pollution area identification and source apportionment of heavy metal contamination in agricultural soil[J]. Journal of Hazardous Materials, 2024, 465(3): 133215.1−133215.10. doi: 10.1016/j.jhazmat.2023.133215
[12] 匡荟芬, 胡春华, 吴根林, 等. 结合主成分分析法(PCA)和正定矩阵因子分解法(PMF)的鄱阳湖丰水期表层沉积物重金属源解析[J]. 湖泊科学, 2019, 32(4): 964−976. doi: 10.18307/2020.0406 Kuang H F, Hu C H, Wu G L, et al. Analysis of heavy metal sources in surface sediments of Poyang Lake in wet period by combining principal component analysis (PCA) and positive definite matrix factorization (PMF)[J]. Lake Science, 2019, 32(4): 964−976. doi: 10.18307/2020.0406
[13] 汪春鹏, 尤建功, 孙浩, 等. 辽阳市土壤重金属含量特征及潜在风险评价[J]. 地质通报, 2021, 40(10): 1680−1687. doi: 10.12097/j.issn.1671-2552.2021.10.010 Wang C P, You J G, Sun H, et al. Characteristics and potential risk assessment of soil heavy metal content in Liaoyang City[J]. Geological Bulletin, 2021, 40(10): 1680−1687. doi: 10.12097/j.issn.1671-2552.2021.10.010
[14] 丁宇雪, 初禹, 金晶泽, 等. 东北地区自然资源监测与黑土退化研究[M]. 武汉: 中国地质大学出版社, 2021: 199−205. Ding Y X, Chu Y, Jin J Z, et al. Study on natural resources monitoring and degradation of black soil in Northeast China[M]. Wuhan: China University of Geosciences Press, 2021: 199−205.
[15] Wu B, Li L L, Guo S H, et al. Source apportionment of heavy metals in the soil at the regional scale based on soil-forming processes[J]. Journal of Hazardous Materials, 2023, 448(1): 130910−130915. doi: 10.2139/ssrn.4091461
[16] 赵秀芳, 张永帅, 冯爱平, 等. 山东安丘地区农业土壤重金属元素地球化学特征及环境评价[J]. 物探与化探, 2020, 44(6): 1446−1454. Zhao X F, Zhang Y S, Feng A P, et al. characteristics and environmental evaluation of heavy metal elements in agricultural soils in Anqiu, Shandong Province[J]. Geophysical and Chemical Exploration, 2020, 44(6): 1446−1454.
[17] Hou S N, Na Z, Lin T. Effect of soil pH and organic matter content on heavy metals availability in maize (Zea mays L. ) rhizospheric soil of non-ferrous metals smelting area[J]. Environmental Monitoring and Assessment, 2019, 10(191): 1−10. doi: 10.1007/s10661-019-7793-5
[18] Ali I, Khan I M, Khan M J, et al. Exploring geochemical assessment and spatial distribution of heavy metals in soils of Southern KP, Pakistan: Employing multivariate analysis[J]. International Journal of Environmental Analytical Chemistry, 2020, 17(102): 1−15. doi: 10.1080/03067319.2020.1804894
[19] 赵岩, 郭常来, 崔健, 等. 辽宁省锦州市北镇农业区土壤重金属分布特征、生态风险评价及源解析[J/OL]. 中国地质(2023-03-14 ).https://kns.cnki.net/kcms/detail/11.1167.P.20230313.1251.002.html. Zhao Y, Guo C L, Cui J, et al. Distribution characteristics, ecological risk assessment and source analysis of soil heavy metals in Beizhen agricultural area, Jinzhou City, Liaoning Province[J/OL]. Geology in China (2023-03-14).https://kns.cnki.net/kcms/detail/11.1167.P.20230313.1251.002.html.
[20] 刘玖芬, 赵晓峰, 侯红星, 等. 地表基质调查分层及分层测试指标体系设计与构建[J]. 岩矿测试, 2024, 43(1): 16−29. doi: 10.15898/j.ykcs.202310080157 Liu J F, Zhao X F, Hou H X, et al. The surface of the substrate layered and layered testing index system design and construction[J]. Rock and Mineral Analysis, 2024, 43(1): 16−29. doi: 10.15898/j.ykcs.202310080157
[21] 赵君, 饶竹, 王鹏, 等. 黑龙江讷河市富锗土壤地球化学特征及影响因素浅析[J]. 岩矿测试, 2022, 9(4): 642−651. doi: 10.3969/j.issn.0254-5357.2022.4.ykcs202204013 Zhao J, Rao Z, Wang P, et al. Heilongjiang nehe rich germanium soil geochemical characteristics and influencing factors of analyses[J]. Rock and Mineral Analysis, 2022, 9(4): 642−651. doi: 10.3969/j.issn.0254-5357.2022.4.ykcs202204013
[22] 赵恒谦, 常仁强, 金倩, 等. 河北西石门铁矿区土壤重金属污染空间分析及风险评价[J]. 岩矿测试, 2023, 42(2): 371−382. doi: 10.15898/j.carolcarrollnki.11-2131/td.,202203290066 Zhao H Q, Chang R Q, Jin Q, et al. Hebei westone door iron mining area of soil heavy metal pollution of spatial analysis and risk evaluation[J]. Rock and Mineral Analysis, 2023, 42(2): 371−382. doi: 10.15898/j.carolcarrollnki.11-2131/td.,202203290066
[23] Williams A J, Antoine J. Evaluation of the elemental pollution status of Jamaican surface sediments using enrichment factor, geoaccumulation index, ecological risk and potential ecological risk index[J]. Marine Pollution Bulletin, 2020, 157: 111288.
[24] 陈泽华, 焦思, 余爱华, 等. 土壤重金属污染评价方法探析——以南京市为例[J]. 森林工程, 2019, 36(3): 28−36. doi: 10.16270/j.cnki.slgc.2020.03.005 Chen Z H, Jiao S, Yu A H, et al. Evaluation methods of soil heavy metal pollution: A case study of Nanjing[J]. Forest Engineering, 2019, 36(3): 28−36. doi: 10.16270/j.cnki.slgc.2020.03.005
[25] Bing H J, Wu Y H, Zhou J, et al. Historical trends of anthropogenic metals in Eastern Tibetan Plateauas reconstructed from alpine lake sediments over the last century[J]. Chemosphere, 2016, 148: 211−219. doi: 10.1016/j.chemosphere.2016.01.042
[26] 张宪依, 庞成宝, 王安婷, 等. 海南岛表层及深层土壤重金属分布特征及源解析[J]. 现代地质, 2020, 34(5): 970−978. Zhang X Y, Pang C B, Wang A T, et al. Distribution and source analysis of heavy metals in surface and deep soil of Hainan Island[J]. Geoscience, 2020, 34(5): 970−978.
[27] Gao L, Wang Z W, Shan J J, et al. Aquatic environ mental changes and anthropogenic activities reflected by the sedi-mentary records of the Shima River, Southern China[J]. Environmental Pollution, 2017, 224(5): 70−81. doi: 10.1016/j.envpol.2016.12.056
[28] 刘兴旺, 苗万里. 基于多元统计和地统计分析法的县域土壤重金属源解析[C]//中国土壤学会土壤环境专业委员会第十九次会议暨农田土壤污染与修复研讨会, 2017. Liu X W, Miao W L. Analysis of heavy metal sources in county soil based on multivariate statistics and geostatistical analysis[C]//The 19th Meeting of the Soil Environment Professional Committee of the Chinese Soil Society and the Symposium on Agricultural Soil Pollution and Remediation, 2017.
[29] Duan X C, Yu H R, Ye T R, et al. Geostatistical mapping and quantitative source apportionment of potentially toxic elements in top- and sub-soils: A case of suburban area in Beijing, China[J]. Ecological Indicators, 2020, 112(5): 106085.1−106085.11. doi: 10.1016/j.ecolind.2020.106085
[30] Xia R, Zhang S Q, Li J, et al. Spatial distribution and quantitative identification of contributions for nutrient and beneficial elements in top- and sub-soil of Huairou District of Beijing, China[J]. Ecological Indicators, 2023, 154: 110853. doi: 10.1016/j.ecolind.2023.110853
[31] 臧传子, 温汉辉, 蔡立梅, 等. 广东省揭阳市土壤铅的空间分布特征及影响因素[J]. 现代地质, 2019, 35(5): 1425−1432. doi: 10.19657/j.geoscience.1000-8527.2021.24 Zang C Z, Wen H H, Cai L M, et al. Spatial distribution characteristics and influencing factors of soil lead in Jieyang City, Guangdong Province[J]. Geoscience, 2019, 35(5): 1425−1432. doi: 10.19657/j.geoscience.1000-8527.2021.24
[32] 刘忆莹, 裴久渤, 汪景宽. 东北典型黑土区耕地有机质与pH的空间分布规律及其相互关系[J]. 农业资源与环境学报, 2019, 36(6): 738−743. Liu Y Y, Pei J B, Wang J K. Spatial distribution and relationship between organic matter and pH in the typical black soil region of Northeast China[J]. Journal of Agricultural Resources and Environment, 2019, 36(6): 738−743.
[33] Liu N T, Cai X Y, Jia L Y, et al. Quantifying mercury distribution and source contribution in surface soil of Qinghai—Tibetan Plateau using mercury isotopes, environmental science & technology[J]. 2023, 57(14): 5903–5912.
[34] Ma Z W, Chen K, Li Z Y, et al. Heavy metals in soils and road dusts in the mining areas of Western Suzhou, China: A preliminary identification of contaminated sites[J]. Journal of Soils and Sediments, 2016, 16: 204−214. doi: 10.1007/s11368-015-1208-1
[35] 裴小龙, 祝晓松, 冯欣, 等. 基于自然资源统一管理的地表基质模型、分类及调查研究[J/OL].地质通报(2024-06-06). https://link.cnki.net/urlid/11.4648.P.20240605.1410.004. Pei X L, Zhu X S, Feng X, et al. Based on unity of natural resources management of surface matrix model, classification and investigation[J/OL]. Geological Bulletin (2024-06-06). https://link.cnki.net/urlid/11.4648.P.20240605.1410.004.