环境科学  2026, Vol. 47 Issue (9): 6433-6442   PDF    
广西某水源地水土重金属空间分布与生态风险评价
赵武强1,2, 胡俊良1,2, 李剑锋1,2,3, 李堃1,2, 程顺波1,2, 汤朝阳1,2     
1. 中国地质调查局武汉地质调查中心(中南地质科技创新中心),武汉 430205;
2. 中国地质调查局花岗岩成岩成矿地质研究中心,武汉 430205;
3. 辽宁师范大学地理科学学院,大连 116029
摘要: 为探讨水源地保护区水土重金属污染、生态风险及溯源问题,采集广西大同江源头保护区1~3级水系表层土壤和水体样品254件和66件,联用单项污染指数、内梅罗综合污染指数、潜在生态风险指数及主成分分析法对其开展系统地评价工作. 结果表明:研究区土壤重金属平均含量由高到低依次为:Cr>Zn>Pb>As>Cu>Ni>Cd>Hg,16件样品单元素超过《土壤环境质量农用地土壤污染风险管控标准》(GB 15618-2018)筛选值(pH<5.5);河水中不同重金属元素均值远低于地表Ⅰ类水标准. 大同江干流土壤中Cu、Cd和Hg具有分段富集的规律,受到母岩(黄洞口组1~4段)自然风化过程的控制;少数样品存在Cr、As、Pb及Cu、Hg明显偏高的现象,可以与村镇及农田对应;六典金矿区出现多处As污染,系为矿山开发导致. 土壤样品内梅罗综合污染指数(PN)为0.11~10.54,处于警戒值、轻度污染的样品占比为33.86%和49.61%;河水单项污染指数整体远低于1,PN为0.144~1.733,位于安全区的样本占比为96.96%;对比显示,土壤背景值的选取问题是造成上述差异的根本原因. 溯源分析揭示区内Cu、Pb、Zn、Ni和Cr主要为自然源,Cd和Hg主要为生活-农业源,而As主要为矿业源. 总之,广西大同江源头保护区内重金属富集主要受控于基岩风化成土过程,整体处于清洁、无污染的状态.
关键词: 水土重金属      空间分布      生态风险      水源地      大同江      广西     
Spatial Distribution and Ecological Risk Assessment of Heavy Metals in Soil and Water in a Water Source Area of Guangxi
ZHAO Wu-qiang1,2 , HU Jun-liang1,2 , LI Jian-feng1,2,3 , LI Kun1,2 , CHENG Shun-bo1,2 , TANG Chao-yang1,2     
1. Wuhan Center of Geological Survey, China Geological Survey (Central South Innovation Center of Geological Science and Technology), Wuhan 430205, China;
2. Research Center of Granite Diagenesis and Metallogenic Geology, China Geological Survey, Wuhan 430205, China;
3. College of Geographical Sciences, Liaoning Normal University, Dalian 116029, China
Abstract: To investigate the issues of heavy metal pollution in soil and water, ecological risks, and pollution source identification in the water source protection area, 254 surface soil samples and 66 water samples were collected from the 1st to 3rd grade water systems in the source protection area of the Datong River, Guangxi Zhuang Autonomous Region. A systematic assessment was conducted on these samples by integrating the single pollution index, Nemerow comprehensive pollution index, potential ecological risk index, and principal component analysis. The results showed that the average contents of heavy metals in the soil of the study area followed the order: Cr>Zn>Pb>As>Cu>Ni>Cd>Hg. Among them, 16 samples exceeded the screening values (pH<5.5) specified in the Soil Environmental Quality Risk Control Standard for Soil Contamination of Agricultural Land (GB 15618-2018). The average contents of various heavy metals in river water were much lower than the standards for Class Ⅰ surface water. In the soil of the main stream of Datong River, Cu, Cd, and Hg exhibited a pattern of segmented enrichment, which was controlled by the natural weathering process of parent rocks (Members 1-4 of the Huangdongkou Formation). A small number of samples showed significantly high concentrations of Cr, As, Pb, as well as Cu and Hg, which could be correlated with the distribution of villages, towns, and farmlands. Multiple occurrences of As pollution were detected in the Liudian gold mining area, which were attributed to mining activities. For the soil samples, the Nemerow comprehensive pollution index (PN) ranged from 0.11 to 10.54, with samples at the warning level and mild pollution level accounting for 33.86% and 49.61%, respectively. For the river water samples, the single pollution indices were generally far below 1, and the PN values ranged from 0.144 to 1.733, with 96.96% of the samples falling into the safe category. Comparative analysis indicated that the selection of soil background values was the fundamental cause of the aforementioned differences. Source identification analysis revealed that Cu, Pb, Zn, Ni, and Cr in the study area were mainly from natural sources; Cd and Hg were primarily from domestic-agricultural sources; and As was mainly from mining sources. In conclusion, the enrichment of heavy metals in the source protection area of the Datong River, Guangxi Zhuang Autonomous Region, is mainly controlled by the weathering of bedrock and soil-forming processes, and the overall quality of the study area is in a clean and unpolluted state.
Key words: heavy metals in soil and water      spatial distribution      ecological risk      water source area      Datong River      Guangxi     

土壤和水体中重金属含量受自然过程和人为活动的综合影响[1~4],表现出明显的时空差异[5,6]. 高背景值地区,自然成土过程即可使重金属富集到对环境和生态系统造成有害影响的水平;而人类长期不规范的生活、生产活动可能导致重金属活化、迁移和富集[7~9];两者均可能导致土壤质量、生产力下降及水质恶化[10]. 生理学研究表明,重金属具有一定的生物毒性;可经由生物摄入途径(消化道、呼吸系统及表皮渗透)侵入人体,对后者身体健康构成潜在威胁[11~15]. 考虑到重金属污染具有易累积、潜伏期长、难降解的特点[16~18],对其动态监测和治理逐渐成为环境科学和地球科学研究的热点[19,20]. 然而,当前对不同研究区/功能区的重金属调查工作多数针对土壤展开,且对其空间分布、污染及风险评价多以数值方法插值直接圈定等值线的形式予以表征;尤其以地形起伏较大的华南丘陵区为大比例尺研究对象时,上述方式在一定程度上忽略了地貌形态的阻隔效应,对其污染范围、程度存在潜在的夸大. 此外,重金属在土壤中存在不同的赋存状态,水文条件对其分布亦有重要的影响[21~24];因此,配套分析水体中重金属含量具有极强的科学性、互补性和有效性[20].

大同江源头保护区是广西壮族自治区重要的水源地和生态保护屏障. 考虑到河水与土壤同为此类功能区的组成要素,对其开展系统的重金属调查与评价工作具有重要的现实与指导意义. 研究区地形起伏较大,属于丘陵地貌;地层及岩性组成简单,以陆源碎屑岩为主;土地类型以林地为主,道路、居民区和农业土地沿水系分布,在其南、北两段分别有停采的小型金矿床. 因此,从水系视角系统探讨区内水土重金属相关问题具有相对便捷和科学的优势,且益于溯源工作;以期为保障水源地的生态环境安全提供依据.

1 材料与方法 1.1 研究区概况

大同江发源于金秀县观音山北麓,流域面积1 139 km2,干流总长106 km,是广西贵港市平南县第二大河流;流经马练乡、同和镇,与大黎河汇合后注入蒙江[图 1(a)和1(b)]. 大同江源头保护区主体位于马练瑶族乡,其主流为研究区内唯一的三级水系;上游及其支流源自大瑶山国家级自然保护区,居民点和少量农业土地均沿大同江及其支流分布,整体以林地为主[图 1(c)]. 内区森林资源丰富,平均气温19.5~20.2 ℃,年降雨量1 562~1 993 mm,属亚热带季风气候. 研究区地层自北向南依次为寒武系黄洞口组1~4段,由灰绿色长石石英砂岩、不等粒砂岩与粉砂质页岩呈不等厚互层构成;上述岩石组合在区内构成向上由粗变细的4个沉积旋回. 区内断裂构造以NNE向为主,控制了区内二级水系(A、B)和多数一级水系的展布. 此外,在大同江东侧支流的一级水系b、c、i和m处具有两处小型金矿区,因环保及水源地要求已关闭.

图 1 大同江源头保护区水土采样位置示意 Fig. 1 Sampling locations of soil and water in the source protection area of the Datong River

1.2 样品采集与分析

结合区内水文条件,笔者将一级水系依次进行编号(a~r),二级水系主要为A、B两处,大同江干流为三级水系. 本次工作采集土壤样品254件,采样深度为0~20 cm;采样位置见图 1(c). 土壤样品采用多点采样法获取,去除植物根系和碎石等杂质,留取1kg装入样品袋. 干燥后研磨过100目尼龙筛,取筛下粒级样品150 g备用. 样品测试在中国冶金地质总局山东局测试中心完成,Cu、Pb、Zn、Cd、Ni和Cr含量采用电感耦合等离子体质谱法(ICP-MS)测定,As和Hg含量使用原子荧光光谱法(AFS)测定;Cu、Pb、Zn、Cr、Ni和As元素检出限为0.01 mg·kg-1,Cd和Hg元素检出限分别为0.001 mg·kg-1. 所有样品分析测试一级标准物质合格率为100%,重复样合格率均>90%,样品RSD<1%.

本文66件地表水体样品取样位置见图 1(c). 采样点坐标采用便携式GPS记录,每件样品均在样点周围3次取样混合而成. 样品测试在中国冶金地质总局山东局测试中心完成,先将水样置入Millipore硝酸纤维素过滤器进行过滤;然后使用ICP-MS 7700x分析微量元素质量浓度,本文分析的微量元素包括Cu、Pb、Zn、Cr、Ni、Cd、As和Hg,其实际检测限分别为0.1、0.1、0.2、0.2、0.02、0.001、0.2和0.01 μg·L-1. 为了保证测得数据的准确性,每测定10个样品检测一次标样,样品RSD<10%.

1.3 评价方法 1.3.1 单项污染指数法

计算公式为:Pi=Ci/Cn

式中,Pi为单项污染指数,Ci为元素i的实测含量,Cn为元素i的背景值. 本文分别选取广西壮族自治区土壤背景值[25]和《土壤环境质量农用地土壤污染风险管控标准》(GB 15618-2018)筛选值[26]进行计算与对比;地表水体样品背景值选取地表水环境质量1类水标准值[27]. 根据单项污染指数可以将污染划分为:未污染(Pi<1)、轻度污染(1≤Pi<2)、中度污染(2≤Pi<3)和重度污染(Pi≥3)这4个等级.

1.3.2 内梅罗综合污染指数法

内梅罗综合污染指数法是在单项污染指数基础上将不同元素所占的权重值纳入计算,考虑了不同元素对水土环境造成不同影响的现实问题,克服了平均法各种污染元素分担的缺陷,能反映污染最严重的元素给土壤环境造成的危害[28,29]. 其计算公式为:

PN=Pave2+Pmax22

式中,PN为内梅罗综合污染指数,Pave为单项污染指数平均值,Pmax为单项污染指数最大值. 内梅罗综合污染指数将污染分为:安全(PN<0.7)、警戒值(0.7≤PN<1)、轻度污染(1≤PN<2)、中度污染(2≤PN<3)及重度污染(PN≥3)这5个等级.

1.3.3 潜在生态风险评价

潜在生态危害指数法是将生态效应、毒理效应和环境效应联系在一起,可系统、全面的评估土壤中元素的污染状况. 计算公式为:

RI=∑i=1nEri=∑i=1nTri·Cfi=∑i=1nTri·CsiCni

式中,Cfi为元素i的污染系数,Csi为元素i的实测含量,Cni为元素i的土壤背景值,Eri为元素i单项潜在生态风险指数,Tri为元素i的毒性响应系数,RI表示综合潜在生态风险指数. 土壤样品选择《土壤环境质量农用地土壤污染风险管控标准》(GB 15618-2018)筛选值[26]为土壤元素背景值,水体样品选取地表水环境质量1类水标准[27]作为背景值. 重金属毒性响应系数及潜在生态危害指数分级见文献[20].

1.3.4 主成分分析法

PCA可实现对高维变量进行降维效果,以反映多指标变量的主要信息[18]. 主要步骤如下:①建立原始变量矩阵,进行标准化处理;②计算标准化数据的相关系数矩阵及其特征根,确定主成分个数;③确定主成分,结合资料进行溯源分析.

1.3.5 PMF模型

正定矩阵因子分解模型(PMF)是一种多元因子分析模型;能够将样本的数据矩阵X分解为因子贡献矩阵G和因子谱矩阵F,定量计算样本因子贡献. 此方法针对未知源成分谱分析,采用非负约束和定量化结合,广泛应用于解析重金属溯源及其贡献[30~33].

2 结果与讨论 2.1 重金属空间分布特征

土壤样品重金属含量分析结果见表 1和图 2(a). 其中,Pb、As和Hg平均含量略超出广西壮族自治区土壤环境背景值;Cu、Pb、Zn、Cr、Ni、Cd、As和Hg样本超标率分别为:1.18%(3件)、54.72%(139件)、3.94%(10件)、28.35%(72件)、5.12%(13件)、0.78%(2件)、35.98%(66件)和42.91%(109件). 参考《土壤环境质量农用地土壤污染风险管控标准》(GB 15618-2018)筛选值(pH<5.5),仅有部分样品存在Pb、Cr、Ni和As含量超过筛选值的情况,分别占比0.39%(1件)、0.79%(2件)、0.39%(1件)和6.30%(16件),以As超标率为最高;余下元素均未超标. 土壤重金属变异系数由大到小依次为:As>Hg>Cd>Ni>Pb>Cu>Zn>Cr,仅As的变异系数为1.37>1,明显大于其它元素.

表 1 研究区土壤重金属含量与Pi统计1) Table 1 Statistics of soil heavy metal content and Pi in the study area

图 2 研究区水土重金属含量与Pi箱线图 Fig. 2 Boxplots of heavy metal contents in soil and water and Pi in the study area

地表水体样品重金属质量浓度分析结果见表 2和图 2(b). 其中,部分样品Cr6+给出40 μg·L-1的分析值,多数样品质量浓度低于检测下限;全部样品Ni质量浓度低于检测下限,在后续分析中不予讨论. 参考地表水环境Ⅰ类水质标准,Cu、Pb、Zn、Cd、As和Hg平均含量均明显优于Ⅰ类水要求;Pb、Zn、As和Hg样本超标率分别为1.52%(1件)、3.03%(2件)、1.52%(1件)和1.52%(1件),其中,1件样品As元素达到重度污染水平,余下4件超标样品为轻度污染. 水体样品重金属变异系数由大到小依次为:As>Cu>Zn>Cd>Pb>Hg,其中,As元素变异系数为7.015,同样明显大于其它元素.

表 2 研究区水体重金属质量浓度与Pi统计1) Table 2 Statistics of heavy metal content and Pi in the water body of the study area

众所周知,同一水系或连通水系上下游的重金属分布可能伴随溶解-迁移-富集等成因联系[19,20],从水系视角亦可以将水土重金属分布进行耦合讨论,这是本研究的优势;其难点在于表征重金属空间分布问题,而选取水系剖面是比较方便、直观的表现形式. 为探讨大同江源头保护区代表性地段的重金属元素的空间分布规律,笔者选取三级水系(A-A′)和两个一、二级水系(B-B′、C-C′)剖面(图 1)统计其土壤重金属含量曲线(图 3). 结果表明:① A-A′剖面中Zn、Cr、Ni,Pb、As,Cu、Cd和Hg 3组元素的变化趋势几乎一致,Cu、Cd和Hg组合的分段富集规律与黄洞口组地层1~4段的界线完全一致,反映不同岩石(母岩)对土壤中重金属丰度的控制,换而言之,这是地质背景的自然源信号;其它元素组合变化整体平稳,波动幅度略小,同样折射其自然源属性[29,34];余下5种元素存在单个样品偏高甚至超过筛选值的情况,可以与居民点和农田对应,暗示局部存在人为污染的情况[35~38]. ② B-B′剖面系由二级水系A、一级水系b组成,后者位于六典金矿区;以发育多处As污染为特征,与A-A′剖面相比显示As元素明显富集,系为矿山开发导致[39~44];Cu和Hg元素的变化趋势与其它两个剖面相比明显缩窄,是因该剖面方向与地层延伸方向一致,即:母岩未发生明显的周期变化. ③ C-C′剖面由二级水系B、一级水系d组成,以多处村庄和农田为特色;在地球化学视角,以最显著的单点异常为特征;Cu、Cd和Hg元素的变化趋势可以与A-A′剖面对比. ④水体重金属分析显示,一级水系m的上游处于清洁状态,流经三联金矿区后重金属含量增加,而后快速衰减,汇流至大同江即为无污染;这以分布规律与土壤剖面完全一致;暗示六典金矿区对大同江源头保护区的重金属富集影响有限[12]. ⑤在三联金矿区,土壤中的重金属含量(一级水系m)整体低于上述3个剖面;但在河水样品中,存在一个重金属含量最高的样品,支持重金属污染具有隐蔽性和复杂性的特点[18,29].

图 3 研究区土壤重金属含量空间分布 Fig. 3 Spatial distribution of heavy metal contents in soil in the study area

2.2 重金属污染评价

单项污染指数均值从大到小依次为:Pb>Hg>As>Cr>Ni>Zn>Cu>Cd,前3项属轻度污染,后5项处于未污染区间. 统计显示,仅Pb整体以轻度污染为主,Hg等7个元素以未污染为主[表 1和图 2(c)];具体而言,Cu、Zn和Cd未污染样品占比为98.82%、96.06%和99.21%;Cr和Ni未污染样品占比为71.65%和94.88%,轻度污染占比26.38%和4.72%;Pb、As和Hg含有4个污染等级的样品,未污染样品占比为45.28%、74.02%和57.09,轻度污染占比为53.54%、17.72%和37.01%,中度污染占比为0.79%、2.36%和4.72%,余下为重度污染占比. 总之,区内土壤重金属污染程度很低,Pb、Hg、As是相对严重的污染因子. 需要补充的是,如以农用地土壤污染风险管控标准计算单项污染指数,区内仅存在个别污染;因此,采用广西壮族自治区背景值计算单项污染指数更利于反映本区土壤重金属富集和污染程度[图 2(c)]. 区内地表水体样品单项污染指数均值依次为:Hg>Zn>Cu>As>Cd>Pb,均值最大的Hg元素Pi=0.243,远低于1. 统计显示,66件水环境样品中8种元素均以未污染为主,除低于检测下限的Cr和Ni外,Pb和Cd处于未污染状态;Cu、Zn和Hg分别有1、2和1件样品为轻度污染,As有1件样品为中度污染[表 2和图 2(d)]. 因此,区内土壤、水体整体处于清洁、安全的状态.

大同江源头保护区土壤中重金属内梅罗综合污染指数为0.11~10.54,平均值为1.28;处于安全区间样本20件,占比为7.87%;处于警戒值样本86件,占比为33.86%;处于轻度污染区间样本126件,占比为49.61%;处于中度污染样本11件,占比为4.33%;处于重度污染样本11件,占比4.33%[图 4(a)]. 区内水体中重金属内梅罗综合污染指数介于0.144~1.733,平均0.233;处于安全区间样本64件,占比为96.96%;处于警戒值样本1件,占比为1.52%;处于轻度污染区间样本1件,占比为1.52%[图 4(b)]. 在空间分布视角,土壤中重金属富集的样品大多处于村镇;有研究表明,人类生活、生产排放的废水、固体废物未经科学处理等可能导致土壤遭受重金属污染[45~49];部分富集样品采自村镇附近的水系,证实农业活动中的有机肥料、化肥及农药的不合理使用同样可能引起重金属在土壤和水体中过度富集[10,50];一级水系a、b、d源头处即出现重度污染的样品,说明岩石风化的自然过程也可引起重金属富集[20,29]. 此外,六典金矿区土壤重金属的富集程度略高于三联金矿区,河水反映前者未遭受重金属污染,但研究区内仅有的1个轻度污染的水样品却位于三联矿区的三级水系m;具体机制有待后续研究.

图 4 研究区水土重金属PN和RI空间分布特征 Fig. 4 Spatial distribution characteristics of PN and RI for heavy metals in soil and water in the study area

基于大同江源头保护区不同水系土壤和地表水体中重金属含量、空间分布及污染评价结果,结合区域地质、地形和水文条件可知:个别土壤样品存在Pb、Ni、Cr和As超过农业土地筛选值标准,可能受到居民生活、农业及矿业活动的综合影响所致. 考虑到不同水系代表一个小的封闭系统,水体重金属显示区内几乎未发生重金属污染;因此,上述土壤孤点式超标对周边影响有限. 换而言之,区内居民生活与农业生产虽导致重金属出现局部富集现象,多数情况并未达到有害的水平[12]. 图 4(a)显示研究区土壤和水体中重金属分别呈轻度污染和安全的状态,表现明显的解耦;主要原因是二者分别基于广西壮族自治区土壤重金属背景值和地表水环境质量Ⅰ类水标准值计算作图. 考虑到重金属在土壤中存在多种赋存状态,进入人体的重金属也无法被人体全部吸收[19,29];因此,大同江源头保护区土壤与水体重金属整体处于安全、可控的状态.

2.3 潜在生态风险评价

土壤单元素潜在生态风险均值由大到小依次为:Cd(13.02)>Hg(5.04)>As(4.00)>Pb(1.89)>Ni(1.52)>Cu(1.35)>Cr(0.99)>Zn(0.21),除了As元素有1件样品为中度生态风险以外;余下全部样品均属轻微生态风险. 计算表明,254件土壤样品综合潜在生态风险指数(RI)为11.07~83.09,平均值为28.03. RI属轻度生态风险样本232件,占比为91.34%;RI属中度生态风险样本21件,占比为8.27%;RI属重度生态风险样本1件,占比为0.39% [图 4(b)]. 其中,1件重度生态风险样品和1件中度生态风险样品处于六典金矿区,与矿业开发有关;余下中度生态风险样品位于村镇、农田或水系起点区域,分别受到人类活动和岩石风化过程控制. 河水单元素潜在生态风险平均值由大到小依次为:Hg(9.73)>Cd(1.09)>Cu(1.3)>As(0.42)>Pb(0.13)>Zn(0.12),除1件样品Hg元素属中度生态风险外,余下样品均属轻微生态风险. 计算表明,全部水环境样品综合潜在生态风险指数(RI)为8.62~54.81,平均值为12.07;均属轻度生态风险[图 4(b)]. 整体而言,综合潜在生态风险指数(RI)采用农用地土壤污染风险管控标准进行计算明显降低了生态风险程度,与按地表水Ⅰ类水标准计算的水体内梅罗综合污染等级和潜在生态风险指数趋于一致.

2.4 重金属污染溯源

为探究区内重金属的来源归属,本文对研究区土壤、河水中重金属元素进行主成分分析,Kaiser-Meyer-Olkin检验结果为0.789和0.766,Bartlett球形检验结果为0.001和<0.000;证实本文所测数据符合主成分分析法的要求,分析结果见图 5. 土壤样品主要共提取3个特征值:① PC1方差贡献率累积32.0%,具有较高的贡献率;主要变量为Cu、Pb、Zn、Cr和Ni;② PC2方差贡献率累积24.6%,主要变量为Zn、Cd和Hg;③ PC3方差贡献率累积14.6%,主要变量为Pb和As [图 5(a)]. 河水样品主要共提取3个特征值:① PC1方差贡献率累积38.5%,具有较高的贡献率;主要变量为Cu、Zn、Cr和Hg;② PC2方差贡献率累积17.7%,主要变量为Pb、Cd、As和Hg;③ PC3方差贡献率累积16.4%,主要变量为As[图 5(b)]. 整体而言,大同江源头保护区土壤、水体重金属具有类似的主控因素.

图 5 研究区水土重金属主成分分析 Fig. 5 Analysis of the main components of heavy metals in the soil and water of the study area

正定矩阵因子分解(PMF)模型可以通过合理设置关键参数(如因子个数、迭代次数和误差阈值等),确保来源解析结果的科学性与可靠性. 本文分析同样判定区内土壤重金属主要有3个来源(图 6),分别贡献51.13%、33.67%和15.19%. 其中,因子1载荷元素主要为Cd和Hg,贡献率为97.29%和79.32%,余下6项重金属元素中最高贡献率的元素为Zn,与土壤、水体重金属主成分PC2组合一致;鉴于其富集区在空间上与村镇及周边农业活动有关,解释为生活-农业源. 因子2载荷元素主要为Cu、Pb、Zn、Ni和Cr,贡献率分别为51.38%、43.83%、32.21%、54.49%和55.28%,与主成分PC1组合基本一致;结合重金属空间分布与污染评价结果,该组元素因黄洞口组碎屑岩自然风化-蚀变成土作用而富集. 因子3主要载荷元素为As,占比为86.77%,与主成分PC3一致;解释为矿业源. 在地质学视角,岩石矿化和蚀变强度达到工业品位即为矿床,尚未达到工业品位的岩石则为矿化体;通常,矿区的矿化、蚀变范围大于工业矿体,这些岩石风化成土之后,在局部地形、水文条件作用下可能重金属富集区域进一步扩大[20,29];此外,矿山开采和尾矿处理不当等过程会加快蚀变岩石与矿石矿物的风化速率,会加速重金属的释放过程[10~12,49];六典金矿区及附近的土壤重金属的分布与污染可能是上述过程综合作用的结果,最好的指示元素为As.

图 6 基于PMF模型的研究区周边土壤重金属溯源解析 Fig. 6 Source apportionment of soil heavy metals around the study area based on PMF model

3 结论

(1)研究区254件土壤样品重金属平均含量由高到低依次为:Cr>Zn>Pb>As>Cu>Ni>Cd>Hg,139件、72件样品Pb、Cr超出广西壮族自治区土壤环境背景值,16件样品的Pb、Cr、Ni、As含量超过《土壤环境质量农用地土壤污染风险管控标准》(GB 15618-2018)筛选值(pH<5.5);66件河水样品中全部元素均值远低于地表Ⅰ类水标准.

(2)土壤重金属空间分布研究显示,大同江干流Zn、Cr、Ni,Pb、As,Cu、Cd和Hg 3组元素的变化趋势几乎一致,Cu、Cd、Hg元素分段富集规律与黄洞口组地层1~4段的界线完全一致,反映自然源的影响;余下元素存在单个样品明显偏高的特征,可与村镇及农田对应;三级水系b位于六典金矿区,以发育多处As污染为特征,系为矿山开发导致.

(3)单项污染指数揭示区内土壤中Pb以轻度污染为主,余下元素以未污染为主;Pb、Hg和As是相对严重的污染因子;内梅罗综合污染指数(PN)为0.11~10.54,处于警戒值和轻度污染的样品占比为33.86%和49.61%;水体中单项污染指数整体远低于1,PN为0.144~1.733,位于安全区间的样本占比96.96%;造成水、土壤污染结果出现解耦的原因为土壤背景值的选取问题.

(4)主成分和PMF模型揭示区内Zn、Ni和Cr为主要自然源,As和Pb主要为矿业源,而Cu、Cd和Hg则为生活-农业源. 大同江源头保护区土壤与地表水体中重金属富集主要受控于自然过程,整体处于安全状态.

参考文献
[1] Saavedra-Mella F, Liu Y J, Southam G, et al. Phosphate treatment alleviated acute phytotoxicity of heavy metals in sulfidic Pb-Zn mine tailings[J]. Environmental Pollution, 2019, 250: 676-685. DOI:10.1016/j.envpol.2019.04.100
[2] Sun X D, Liu Y, Guo S, et al. Interregional supply chains of Chinese mineral resource requirements[J]. Journal of Cleaner Production, 2021, 279. DOI:10.1016/j.jclepro.2020.123514
[3] 周子寒, 张云霞, 许仁智, 等. 广西岩溶地区铅锌矿区周边优势植物重金属富集特性及应用潜力[J]. 环境科学, 2025, 46(10): 6621-6629.
Zhou Z H, Zhang Y X, Xu R Z, et al. Accumulation and application potentials of dominant plants on Pb-Zn mines in karst region, Guangxi[J]. Environmental Science, 2025, 46(10): 6621-6629. DOI:10.13227/j.hjkx.202409126
[4] Hou Y R, Zhao Y Y, Lu J L, et al. Environmental contamination and health risk assessment of potentially toxic trace metal elements in soils near gold mines-a global meta-analysis[J]. Environmental Pollution, 2023, 330. DOI:10.1016/j.envpol.2023.121803
[5] Liu Z, Du Q Q, Guan Q Y, et al. A Monte Carlo simulation-based health risk assessment of heavy metals in soils of an oasis agricultural region in northwest China[J]. Science of the Total Environment, 2023, 857. DOI:10.1016/j.scitotenv.2022.159543
[6] 刘展. 河西地区重点工农业区农田土壤重金属污染及风险研究[D]. 兰州: 兰州大学, 2023.
Liu Z. Study on heavy metal contamination and risk of agricultural soils in key industrial and agricultural areas in the Hexi region[D]. Lanzhou: Lanzhou University, 2023.
[7] 甘杰, 邹霖, 陈贝贝, 等. 湘南地区土壤重金属污染分布特征及特定源-生态风险评估[J]. 环境科学, 2025, 46(1): 430-441.
Gan J, Zou L, Chen B B, et al. Distribution and source-specific ecological risk analysis of soil heavy metals in South Hunan Province[J]. Environmental Science, 2025, 46(1): 430-441. DOI:10.13227/j.hjkx.202401210
[8] Yang S Y, Sun L J, Sun Y F, et al. Towards an integrated health risk assessment framework of soil heavy metals pollution: theoretical basis, conceptual model, and perspectives[J]. Environmental Pollution, 2023, 316. DOI:10.1016/j.envpol.2022.120596
[9] Sheng D R, Meng X H, Wen X H, et al. Contamination characteristics, source identification, and source-specific health risks of heavy metal(loid)s in groundwater of an arid oasis region in Northwest China[J]. Science of the Total Environment, 2022, 841. DOI:10.1016/j.scitotenv.2022.156733
[10] 马常莲. 新疆若羌县绿洲带地下水与农用地土壤环境质量及重金属健康风险评价[D]. 乌鲁木齐: 新疆农业大学, 2023.
Ma C L. Environmental quality and health risk assessment of heavy metals of groundwater and agricultural land in oasis belt of Ruoqiang County, Xinjiang[D]. Urumqi: Xinjiang Agricultural University, 2023.
[11] Al Osman M, Yang F, Massey I Y. Exposure routes and health effects of heavy metals on children[J]. BioMetals, 2019, 32(4): 563-573. DOI:10.1007/s10534-019-00193-5
[12] Bernard A. Cadmium & its adverse effects on human health[J]. The Indian Journal of Medical Research, 2008, 128(4): 557-564.
[13] Duruibe J O, Ogwuegbu M O C, Egwurugwu J N. Heavy metal pollution and human biotoxic effects[J]. International Journal of Physical Sciences, 2007, 2(5): 112-118.
[14] Zahra N, Kalim I, Mahmood M, et al. Perilous effects of heavy metals contamination on human health[J]. Pakistan Journal of Analytical & Environmental Chemistry, 2017, 18(1): 1-17.
[15] Lu J, Wu J, Zhang C, et al. Possible effect of submarine groundwater discharge on the pollution of coastal water: occurrence, source, and risks of endocrine disrupting chemicals in coastal groundwater and adjacent seawater influenced by reclaimed water irrigation[J]. Chemosphere, 2020, 250. DOI:10.1016/j.chemosphere.2020.126323
[16] 李剑锋, 冯李霄. 湖南某锡矿区土壤重金属污染及健康风险评价[J]. 中国地质, 2023, 50(3): 897-910.
Li J F, Feng L X. Health risk assessment of heavy metal pollution in soil of a tin mining area in Hunan Province[J]. Geology in China, 2023, 50(3): 897-910.
[17] 李剑锋, 冯李霄, 陈希清, 等. 大义山东南部土壤重金属分布特征及其风险评价[J]. 环境工程技术学报, 2023, 13(1): 287-294.
Li J F, Feng L X, Chen X Q, et al. Heavy metal distribution characteristics of soils in southeastern Dayi Mountain and its risk evaluation[J]. Journal of Environmental Engineering Technology, 2023, 13(1): 287-294.
[18] 冯李霄. 湖南某矿区土壤重金属与稀土元素污染及健康风险评价[D]. 大连: 辽宁师范大学, 2023.
Feng L X. Evaluation of heavy metal and rare earth element contamination and health risk of soil in a mining area in Hunan Province[D]. Dalian: Liaoning Normal University, 2023.
[19] 陈希清, 秦拯纬, 李剑锋, 等. 湘南骑田岭南坡环境地球化学特征研究[J]. 华南地质, 2021, 37(3): 298-312.
Chen X Q, Qin Z W, Li J F, et al. The environmental geochemical characteristics of the south slope of Qitianling Mountain, southern Hunan Province[J]. South China Geology, 2021, 37(3): 298-312.
[20] Wang L, Li J F, Liu X W, et al. Heavy metal pollution in water and soil and associated health risks in a tin mining region of Hunan Province, Central China[J]. China Geology, 2025, 8(2): 314-324.
[21] Wu Y F, Li X, Yu L, et al. Review of soil heavy metal pollution in China: spatial distribution, primary sources, and remediation alternatives[J]. Resources, Conservation and Recycling, 2022, 181. DOI:10.1016/j.resconrec.2022.106261
[22] 付蓉洁, 辛存林, 于奭, 等. 石期河西南子流域地下水重金属来源解析及健康风险评价[J]. 环境科学, 2023, 44(2): 796-806.
Fu R J, Xin C L, Yu S, et al. Analysis of heavy metal sources in groundwater and assessment of health risks: an example from the southwest sub-basin of the Shiqi River[J]. Environmental Science, 2023, 44(2): 796-806. DOI:10.13227/j.hjkx.202203144
[23] Wen X H, Lu J, Wu J, et al. Influence of coastal groundwater salinization on the distribution and risks of heavy metals[J]. Science of the Total Environment, 2019, 652: 267-277. DOI:10.1016/j.scitotenv.2018.10.250
[24] 庞彩艳, 文琦, 丁金梅, 等. 黄河上游流域生态系统服务变化及其权衡协同关系[J]. 生态学报, 2024, 44(12): 5003-5013.
Pang C Y, Wen Q, Ding J M, et al. Ecosystem services and their trade-offs and synergies in the upper reaches of the Yellow River basin[J]. Acta Ecologica Sinica, 2024, 44(12): 5003-5013.
[25] 中国环境监测总站. 中国土壤元素背景值[M]. 北京: 中国环境科学出版社, 1990.
[26] GB 15618-2018, 土壤环境质量农用地土壤污染风险管控标准(试行)[S].
[27] GB 3838-2002, 地表水环境质量标准[S].
[28] Xia P H, Ma L, Sun R G, et al. Evaluation of potential ecological risk, possible sources and controlling factors of heavy metals in surface sediment of Caohai Wetland, China[J]. Science of the Total Environment, 2020, 740. DOI:10.1016/j.scitotenv.2020.140231
[29] 赵武强, 李剑锋, 李堃, 等. 广西某锡矿区周边土壤重金属污染及风险评价[J/OL]. 中国地质, 1-12. https://link.cnki.net/urlid/11.1167.P.20250902.1625.010, 2025-09-03.
Zhao W Q, Li J F, Li K, et al. Heavy metal pollution and risk assessment of soils around a tin mining area in Guangxi[J/OL]. Geology in China, 1-12. https://link.cnki.net/urlid/11.1167.P.20250902.1625.010, 2025-09-03.
[30] 朱舒, 李舒怡, 鱼洋, 等. 关中平原农田土壤剖面重金属空间分布、影响因素、风险评价及源解析[J]. 环境科学, 2025, 46(4): 2501-2511.
Zhu S, Li S Y, Yu Y, et al. Spatial distribution, influencing factors, risk assessment, and source analysis of heavy metals in soil Profile of farmland in the Guanzhong plain[J]. Environmental Science, 2025, 46(4): 2501-2511. DOI:10.13227/j.hjkx.202403016
[31] 李雨茹, 段伟杰, 董嘉臻, 等. 基于PMF模型和地理探测器的金陵河下游土壤重金属源风险评价[J]. 环境科学, 2025, 46(4): 2512-2524.
Li Y R, Duan W J, Dong J Z, et al. Risk assessment of soil heavy metal sources in the lower reaches of the Jinling River based on PMF model and Geodetector[J]. Environmental Science, 2025, 46(4): 2512-2524. DOI:10.13227/j.hjkx.202404001
[32] Li Z Y, Ma Z W, van der Kuijp T J, et al. A review of soil heavy metal pollution from mines in China: pollution and health risk assessment[J]. Science of the Total Environment, 2014, 468-469: 843-853. DOI:10.1016/j.scitotenv.2013.08.090
[33] 柴磊, 王新, 马良, 等. 基于PMF模型的兰州耕地土壤重金属来源解析[J]. 中国环境科学, 2020, 40(9): 3919-3929.
Chai L, Wang X, Ma L, et al. Sources appointment of heavy metals in cultivated soils of Lanzhou based on PMF models[J]. China Environmental Science, 2020, 40(9): 3919-3929.
[34] 王菲, 吴泉源, 吕建树, 等. 山东省典型金矿区土壤重金属空间特征分析与环境风险评估[J]. 环境科学, 2016, 37(8): 3144-3150.
Wang F, Wu Q Y, Lü J S, et al. Spatial characteristics and environmental risk of heavy metals in typical gold mining area of Shandong Province[J]. Environmental Science, 2016, 37(8): 3144-3150. DOI:10.13277/j.hjkx.2016.08.041
[35] Chen L, Zhou M X, Wang J Z, et al. A global meta-analysis of heavy metal(loid)s pollution in soils near copper mines: evaluation of pollution level and probabilistic health risks[J]. Science of the Total Environment, 2022, 835. DOI:10.1016/j.scitotenv.2022.155441
[36] 黄颖. 不同尺度农田土壤重金属污染源解析研究[D]. 杭州: 浙江大学, 2018.
Huang Y. The exploring of heavy metal pollution source apportionment in various scale of agricultural soils[D]. Hangzhou: Zhejiang University, 2018.
[37] 杨雳, 白宗旭, 薄文浩, 等. 中国农田土壤重金属污染分析与评价[J]. 环境科学, 2024, 45(5): 2913-2925.
Yang L, Bai Z X, Bo W H, et al. Analysis and evaluation of heavy metal pollution in farmland soil in China: a meta-analysis[J]. Environmental Science, 2024, 45(5): 2913-2925. DOI:10.13227/j.hjkx.202305253
[38] Liu F T, Wang X Q, Dai S, et al. Spatial variations, health risk assessment, and source apportionment of soil heavy metals in the middle Yellow River Basin of northern China[J]. Journal of Geochemical Exploration, 2023, 252. DOI:10.1016/j.gexplo.2023.107275
[39] Liu B L, Ai S W, Zhang W Y, et al. Assessment of the bioavailability, bioaccessibility and transfer of heavy metals in the soil-grain-human systems near a mining and smelting area in NW China[J]. Science of the Total Environment, 2017, 609: 822-829. DOI:10.1016/j.scitotenv.2017.07.215
[40] Sun Z H, Xie X D, Wang P, et al. Heavy metal pollution caused by small-scale metal ore mining activities: a case study from a polymetallic mine in South China[J]. Science of the Total Environment, 2018, 639: 217-227. DOI:10.1016/j.scitotenv.2018.05.176
[41] 孙建伟, 贾煦, 刘向东, 等. 豫西金矿集区矿业活动对周边农田土壤重金属影响研究[J]. 岩矿测试, 2023, 42(1): 192-202.
Sun J W, Jia X, Liu X D, et al. Influence of mining activities in the gold ore concentration area in western Henan on the heavy metals in surrounding farmland soil[J]. Rock and Mineral Analysis, 2023, 42(1): 192-202.
[42] 冯娟, 艾昊, 陈清敏, 等. 秦岭山区某金属矿区土壤重金属污染评价及迁移路径解析[J]. 岩矿测试, 2023, 42(6): 1189-1202.
Feng J, Ai H, Chen Q M, et al. Evaluation and migration path analysis of soil heavy metal pollution in a metal mining area of Qinling Mountain[J]. Rock and Mineral Analysis, 2023, 42(6): 1189-1202.
[43] Meng C Y, Wang P, Hao Z L, et al. Ecological and health risk assessment of heavy metals in soil and Chinese herbal medicines[J]. Environmental Geochemistry and Health, 2022, 44(3): 817-828. DOI:10.1007/s10653-021-00978-z
[44] Han Q, Liu Y, Feng X X, et al. Pollution effect assessment of industrial activities on potentially toxic metal distribution in windowsill dust and surface soil in central China[J]. Science of the Total Environment, 2021, 759. DOI:10.1016/j.scitotenv.2020.144023
[45] Karimyan K, Alimohammadi M, Maleki A, et al. Human health and ecological risk assessment of heavy metal(loid)s in agricultural soils of rural areas: a case study in Kurdistan Province, Iran[J]. Journal of Environmental Health Science and Engineering, 2020, 18(2): 469-481. DOI:10.1007/s40201-020-00475-y
[46] Wu G F, Wang L L, Yang R, et al. Pollution characteristics and risk assessment of heavy metals in the soil of a construction waste landfill site[J]. Ecological Informatics, 2022, 70. DOI:10.1016/j.ecoinf.2022.101700
[47] Zhou H, Chen Y, Yue X M, et al. Identification and hazard analysis of heavy metal sources in agricultural soils in ancient mining areas: a quantitative method based on the receptor model and risk assessment[J]. Journal of Hazardous Materials, 2023, 445. DOI:10.1016/j.jhazmat.2022.130528
[48] Kumari M, Kumar A, Bhattacharya T. Assessment of heavy metal contamination in street dust: concentrations, bioaccessibility, and human health risks in coal mine and thermal power plant complex[J]. Environmental Geochemistry and Health, 2023, 45(10): 7339-7362. DOI:10.1007/s10653-023-01695-5
[49] Macklin M G, Thomas C J, Mudbhatkal A, et al. Impacts of metal mining on river systems: a global assessment[J]. Science, 2023, 381(6664): 1345-1350. DOI:10.1126/science.adg6704
[50] 李军, 高占栋, 马利邦, 等. 黄河上游白银段煤矿区周边土壤重金属累积特征、污染风险与溯源解析[J]. 环境科学, 2025, 46(4): 2475-2485.
Li J, Gao Z D, Ma L B, et al. Accumulation characteristics, pollution risk, and source analysis of soil heavy metals around a coal mining area in the Baiyin section of the upper Yellow River[J]. Environmental Science, 2025, 46(4): 2475-2485. DOI:10.13227/j.hjkx.202404060