环境科学  2024, Vol. 45 Issue (5): 2983-2994   PDF    
基于源导向和蒙特卡洛模型的广东省某城市土壤重金属健康风险评估
陈莲1, 邹子航1,2, 张培珍2, 王雨菡1,2, 王振江1,3, 林森1, 唐翠明1,3, 罗国庆1,3, 钟建武1, 李智毅1, 王圆1     
1. 广东省农业科学院蚕业与农产品加工研究所, 广州 510610;
2. 广东海洋大学电子信息与工程学院, 湛江 524000;
3. 农业农村部华南都市农业重点实验室, 广州 510610
摘要: 以广东省某市为研究区域, 通过研究其表层土壤重金属含量及空间分布特征, 明确土壤重金属污染情况和优先控制因子, 为该市土壤重金属污染防治提供基础数据. 对该市221个土壤样品中的重金属含量特征进行分析, 并通过潜在健康风险评估(HRA)模型、PMF受体模型和蒙特卡洛模型进行潜在健康风险评估、来源解析和主控因子分析. 结果发现, 该市土壤重金属ω(As)、ω(Hg)、ω(Cd)、ω(Pb)、ω(Cr)、ω(Cu)、ω(Ni)和ω(Zn)的平均值分别为18.16、0.43、1.46、68.57、98.34、64.19、26.53和257.32 mg·kg-1, 变异程度为中高度变异. 除Ni元素以外, 土壤其余重金属元素含量均已一定程度超过广东省土壤背景值, 且Cd和Zn的含量已超过国家二级标准, 重金属污染已较为严重;重金属来源主要为工业源, 自然母质、铅蓄电池制造、交通、人为耕作和农药化肥投入也对土壤重金属累积有不可忽视的影响;土壤重金属对儿童和成人均存在一定程度的可耐受致癌健康风险, 非致癌风险可以忽略不计;儿童的潜在健康风险大于成人, 主要暴露途径均为经口摄入;农药化肥投入源和As应该作为该市土壤重金属影响健康风险的主控因子, 其次为混合源和Cr;重金属空间分布特征及相对污染程度存在差异, 应深化分区监测管控, 加强土壤的污染防治, 减少土壤重金属的人为输入.
关键词: 土壤      重金属      健康风险评估      蒙特卡洛      PMF模型      主控因子     
Health Risk Assessment of Heavy Metals in Soils of a City in Guangdong Province Based on Source Oriented and Monte Carlo Models
CHEN Lian1 , ZOU Zi-hang1,2 , ZHANG Pei-zhen2 , WANG Yu-han1,2 , WANG Zhen-jiang1,3 , LIN Sen1 , TANG Cui-ming1,3 , LUO Guo-qing1,3 , ZHONG Jian-wu1 , LI Zhi-yi1 , WANG Yuan1     
1. Research Institute of Sericulture and Agricultural Products Processing, Guangdong Academy of Agricultural Sciences, Guangzhou 510610, China;
2. School of Electronic Information and Engineering, Guangdong Ocean University, Zhanjiang 524000, China;
3. South China Urban Agriculture Key Laboratory of the Ministry of Agriculture and Rural Affairs, Guangzhou 510610, China
Abstract: Taking a city in Guangdong Province as the research area, the concentration and spatial distribution characteristics of heavy metals in the surface soil were studied to clarify the situation of soil heavy metal pollution and priority control factors, providing basic data for the prevention and control of soil heavy metal pollution in the city. The content characteristics of heavy metals in 221 soil samples in the city were analyzed, and the potential health risk assessment and source analysis were carried out through the Monte Carlo model, the potential health risk assessment (HRA) model, and the PMF receptor model. It was found that heavy metals ω(As), ω(Hg), ω(Cd), ω(Pb), ω(Cr), ω(Cu), ω(Ni), and ω(Zn) in the soil of the city were 18.16, 0.43, 1.46, 68.57, 98.34, 64.19, 26.53, and 257.32 mg·kg-1, respectively, with a moderate to high degree of variation. Except for Ni concentration, the soil concentrations of other heavy metal elements exceeded the background values of soil in Guangdong Province to a certain extent, and the concentrations of Cd and Zn exceeded the national secondary standards, resulting in severe heavy metal pollution; the main sources of heavy metals were industrial sources, and natural parent materials, lead battery manufacturing, transportation, artificial cultivation, and pesticide and fertilizer inputs also had an undeniable impact on the accumulation of heavy metals in the soil. Heavy metals in the soil had a certain degree of tolerable carcinogenic health risk for both children and adults, whereas non-carcinogenic risks could be ignored. The potential health risk of children was greater than that of adults, and the main exposure route was through oral intake. The input sources of pesticides and fertilizers and As should be the main controlling factors for the health risks of heavy metals in the city's soil, followed by mixed sources and Cr. There were differences in the spatial distribution characteristics and relative pollution levels of heavy metals, and it is necessary to deepen zoning monitoring and control, strengthen soil pollution prevention and control, and reduce human input of heavy metals in soil.
Key words: soil      heavy metals      health risk assessment      Monte Carlo      PMF model      master factor     

随着经济的快速发展, 城市产业结构不断调整, 人为活动影响剧烈, 土壤重金属污染成为日趋严重的生态环境问题[1]. 土壤重金属污染来源广泛, 包括采矿、冶炼、制革和染料等工业排放的“三废”及汽车尾气排放、农药和化肥的施用等[2~5]. 重金属在土壤中积累到一定程度不仅导致土壤退化、农作物产量和品质降低, 而且通过径流和淋洗作用污染地表水和地下水, 恶化水文环境, 并可通过直接接触、食物链等途径危及人类的生命和健康[6~8]. 我国土壤重金属污染问题极为严重, 全国污染点位超标率达16.1%以上[9, 10]. 广东省是我国重金属污染防治重点省份, 珠三角40%以上的农田菜地重金属污染超标, 其中10%属严重超标[11]. 因此, 土壤重金属污染问题极为普遍, 已成为全球关注的环境焦点.

城市不同功能区在经济发展、产业结构、人口密度及交通状况的差别导致土壤重金属污染的空间特征及潜在风险具有显著差异[12], 如张军等[13]研究发现宝鸡市城区土壤重金属含量、风险等级在空间维度上随着城区面积的扩展, 出现污染区域由城区西南部向东部迁移的特征;夏子书等[14]采用PMF模型对宁南山区小流域土壤研究发现重金属主要来自:施肥和工业排放混合源、农业活动源、自然母质源、农药使用和矿业开发排放混合源和大气沉降源;郑永立等[15]研究发现广东省揭阳市榕城区重金属主要来自交通、成土母质、大气沉降和工业源, 各重金属元素对人体的非致癌健康风险可忽略不计, 对成年和儿童均存在一定程度的可耐受致癌健康风险;郭志娟等[16]研究发现, 雄安新区整体土壤重金属含量属于无污染状态, 各重金属元素对成年的非致癌健康风险可以忽略不计, 部分样点对儿童存在非致癌健康风险, 对成年和儿童均存在一定程度的可耐受致癌健康风险. 但目前对城市土壤重金属污染研究多集中于特定功能区及单个或多个行政区[17~22], 土壤重金属的污染特征及潜在健康风险评估的研究仍存在空白. 基于此, 本研究整合分析某市土壤重金属砷(As)、汞(Hg)、铬(Cr)、镉(Cd)、铜(Cu)、锌(Zn)、铅(Pb)和镍(Ni)的污染水平, 采用健康风险评估模型结合蒙特卡洛模型评估潜在致癌和非致癌健康风险及风险概率, 采用PMF模型对重金属元素进行源解析, 综合健康风险评估和源解析结果评估主控因子, 以期为该市土壤重金属污染管控与防治提供科学依据和理论支持, 推动土壤环境改善进程.

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

研究区域为广东省的重点城市之一, 是中国极具代表性的特大城市, 世界著名港口城市, 发展历史悠久, 人口密集程度高, 2022年该市常住人口已超过1 800万, 交通网络稠密, 车辆保有量已超过300万辆, 经济发展成熟度高, 2022年人均GDP已超过10万元, 研究区历史悠久, 文化底蕴丰富, 居民生产活动对研究区土壤造成了较大生态污染.

1.2 样品布点与采集

在ArcMap 10.8中通过“创建渔网”功能生成规则的10 km × 10 km网格, 对该市进行样点布设, 同时考虑了该市人口密度和土地利用情况, 并结合现场实际情况, 避开沟渠、林带、田梗和路边旧房基等不适合设点区域, 用手持GPS定位仪记录实际采样点坐标位置, 导入ArcMap中生成实际采样点位图. 在人为活动比较密集的区域, 增加采样的密度, 在森林等受人为活动影响较小的区域, 适当降低采样密度, 如图 1所示, 采样时间为2020年12月29日至2021年1月26日. 采样时, 利用“五点收集法”, 剥除表面覆盖物后, 各采样点采集3个0~20 cm表层土壤样品, 均匀混合制成1个土壤样品, 每个样品1 kg, 装入聚乙烯塑料袋送至实验室. 土壤样品自然风干, 去除植物根系、石块及人造物等, 置于BJ-300多功能粉碎机中研磨30 s后过100目0.15 mm孔径筛, 取100 g土样装入聚乙烯样品袋备用.

图 1 土壤采样点位示意 Fig. 1 Soil sampling sites

1.3 土壤重金属含量测定

对土壤样品依次用HNO3-HCl-HF-HClO4溶液消解, 样品处理完成后送至南京大学现代分析中心进行检测. 测试指标包括As、Hg、Cr、Cd、Cu、Zn、Pb和Ni, 其中As和Hg的检测方法为原子荧光光谱测定法[23, 24], Cd采用美国PerkinElmer公司生产的Elan9000电感耦合等离子体质谱仪测定, Cr、Cu、Zn、Ni和Pb采用美国PerkinElmer公司生产的Optima 5300DV型电感耦合等离子体发射光谱仪进行测定[14]. 同时设置空白对照样以保证分析结果准确性.

1.4 蒙特卡洛健康风险评估

根据美国环保署(USEPA)推荐的风险评估方法, 采用健康风险评估(HRA)模型对土壤环境健康风险进行评估, 土壤重金属往往通过3种途径(经口摄入、皮肤接触和吸入土壤颗粒物)被人体摄入并对人体健康具有致癌风险和非致癌风险[25], 计算公式为:

(1)
(2)
(3)

式中,ADDi,经口摄入、ADDi,皮肤接触和ADDi,呼吸摄入分别为人体通过口直接摄入、皮肤接触和呼吸摄入某重金属的日均暴露剂量;Ci为土壤中某重金属的含量,IngR和InhR分别为每日的土壤摄入量和空气吸入量,CF为转换因子,EF为暴露频率,ED为暴露年限,BW为平均体重,SA为暴露皮肤表面积,SL为皮肤粘着度,ABS为皮肤吸收因子,PEF为地表排放因子,AT为平均暴露时间.

(4)
(5)
(6)
(7)

式中, HQi为某重金属的非致癌健康风险指数, HQ为该样点总的非致癌健康风险指数, HQ /HQi ≤ 1表示重金属的非致癌健康风险可以忽略不计, 反之则表示存在非致癌健康风险[16, 26];RfDi为某重金属的日均摄入量;CRi为某重金属的致癌健康风险指数, CR为该样点总的致癌健康风险指数, 当CRi/CR ≤ 1 × 10-6说明致癌健康风险可以忽略不计, 1 × 10-6 < CRi/CR ≤ 1 × 10-4时, 说明存在可耐受致癌健康风险, CRi/CR > 1 × 10-4时说明存在不可耐受致癌风险. 土壤重金属健康风险评估模型中的暴露参数和不同暴露途径的参考剂量, 均参考自中国生态环境部发布的HJ25.3-2019中风险评估模型参数推荐值和《中国人群暴露参数手册(儿童卷)概要》[26]. 本次研究具体暴露参数和剂量参考值见表 1表 2.

表 1 基于蒙特卡洛的重金属健康风险暴露参数 Table 1 Health risk exposure parameters for heavy metals based on Monte Carlo

表 2 重金属不同暴露途径的参考剂量和致癌斜率因子1) Table 2 Reference dose and carcinogenic slope factor for different exposure pathways of heavy metals

1.5 源解析

采用USEPA推荐的PMF 5.0受体模型对该市土壤重金属的来源进行解析, 其原理是将重金属含量数据原始矩阵X通过因子化分解成因子贡献G和因子残差矩阵F, 然后基于各重金属来源的特征差异确定因子贡献率[32].

(8)
(9)
(10)

式中, Xij为第i个样本中重金属j的含量, fkj为重金属j在来源k中的含量, gik为来源k对第i个样本的贡献率, eij为残差矩阵;Q为通过不断分解矩阵X选择最优的矩阵GF得到的最小目标函数;uij为第i个样本中重金属j含量的不确定性, RSD为重金属含量的相对标准偏差, MDL为方法检出限, Ni、Cu、Cr、Zn、Cd、As、Pb和Hg的检出限分别为1.0、1.0、0.1、2.0、0.01、0.2、0.1和0.005 mg·kg-1.

1.6 基于源导向的健康风险评估

基于PMF模型解析获取各重金属元素对不同污染源的贡献结果和健康风险评估获取的非致癌风险和致癌风险结果, 计算不同重金属元素对非致癌风险和致癌风险的贡献率[28], 计算公式为:

(11)
(12)
(13)
(14)

式中, Bn, HQ为第n类污染源中的非致癌健康指数, Bn, CR为第n类污染源中的致癌健康指数, Cin为第n类污染源中重金属i的含量, Dn, HQ为第n类污染源中非致癌健康指数的贡献率, Dn, CR为第n类污染源中致癌健康指数的贡献率.

1.7 数据分析方法

采用Excel表格进行数据整理, SPSS 21进行数据分析和统计性描述, ArcMap 10.8用于空间特征分析与作图, Origin 2021进行相关图形绘制, PMF 5.0进行源解析.

2 结果与分析 2.1 土壤重金属含量特征及空间分布特征

该市土壤中重金属含量的统计性描述如表 3所示. 整个城市的ω(As)、ω(Hg)、ω(Cd)、ω(Pb)、ω(Cr)、ω(Cu)、ω(Ni)和ω(Zn)的平均值分别为18.16、0.43、1.46、68.57、98.34、64.19、26.53和257.32 mg·kg-1;除Ni元素以外, 其余重金属元素土壤含量的平均值均已超过广东省土壤背景值, 且Cd和Zn的含量已超过国家二级标准, 重金属污染已较为严重. 不同土地利用类型中, 重金属Cd、Cu和Ni的含量平均值由高到低为:工矿用地 > 城镇用地 > 耕地 > 林地, 重金属Cr和Pb的含量平均值由高到低为:城镇用地 > 工矿用地 > 林地 > 耕地, 重金属Zn和As的含量平均值由高到低依次为:城镇用地 > 工矿用地 > 耕地 > 林地, 重金属Hg的含量平均值由高到低分别为:城镇用地 > 林地 > 耕地 > 工矿用地, 城镇用地和工矿用地重金属污染程度较高. 变异系数能反映元素在土壤中分布的均匀性和变异程度[33], 重金属Cr属于中度变异, 其余重金属元素属于高度变异, 说明该市土壤受人为活动干扰严重. 在ArcMap 10.8中利用地统计模块的反距离插值分别对8种重金属进行空间特征分析, 结果如图 2. Ni和Cd污染中心值点基本重合, 主要为该市中部地区;As和Hg的高值区主要集中在西部区域和中部部分区域, Cr和Zn的高值区主要集中在中南部区域, 多个重金属元素在南部个别区域出现高值区, 说明该区域已存在点源污染, 可能会对该区域内生态环境和居民生活带来一定程度的危害, 需提高关注度.

表 3 整体土壤重金属含量特征/mg·kg-1 Table 3 Characteristics of heavy metal content in overall soil/mg·kg-1

图 2 某市土壤不同重金属元素空间分布特征 Fig. 2 Spatial distribution characteristics of different heavy metal elements in soil of a certain city

2.2 源解析 2.2.1 重金属相关性分析

相关性分析可以在一定程度上反映土壤中重金属之间的同源性, 相关系数越大, 表明来自同一污染源的可能性就越大[34]. 该市土壤重金属分析结果如表 4. 除Ni与Cd、Pb、Hg和As, Hg与Cd、Cr、Cu、Pb和Zn, As与Cr、Cu和Zn之间无明显相关性之外, 其余重金属元素之间均有一定程度相关;Ni⁃Cr的相关性系数大于0.5, 说明两种重金属元素存在部分同源性;Zn⁃Cu和Zn⁃Cr的相关性系数接近0.5, 说明Zn和Cu, Zn和Cr之间也部分存在同源性, Cu⁃Cr的相关性系数为0.2左右, 同源性相对较弱, 综上所述, 除Cr⁃Ni外, 多个重金属元素之间的相关性相对较弱, 说明污染来源较为复杂, 土壤受不同来源污染综合影响.

表 4 土壤重金属含量相关性分析1) Table 4 Correlation analysis of heavy metal content in soil

2.2.2 PMF受体模型

采用PMF 5.0受体模型对该市整体表层土壤重金属来源进行解析, 根据信噪比的值将8种重金属都设置为strong, 将8种重金属的总体值设置为weak, 选择3~5因子运行50次后, 选择预测值和实际值相差最小的因子数作为最佳运行因子, 最佳运行因子为5. 运行结果如图 3图 4. 将PMF模型中导出的土壤重金属污染源各因子贡献值的点位数据导入ArcMap中, 进行反距离插值后得到结果如图 5. 因子1的主要特征物为As(87.6%), 杀虫剂和化肥的高强度使用会导致土壤中As的积累[34, 35], 结合图 5中因子1的空间分布高值区基本为种植区和绿地, 以及实际调研, 高值点主要集中在绿地、农田以及村庄垃圾堆肥点, 所以因子1可以判断为农药化肥投入源;因子2的贡献成分比较丰富, 对Cd、Cr和Ni均有较大贡献, Cr和Ni的含量平均值基本和广东省土壤背景值接近, 因此土壤中重金属Cr和Ni主要受其成土母质的影响[36], 而Cd的含量平均值已超过广东省背景值的12.3倍, Cd在土壤中的富集程度较高, 水稻土对Cd的富集程度相对其他土壤要高得多, 因此, 土壤中重金属Cd主要受到人为耕作的影响. 电镀、金属冶炼、燃料和化学工业也会产生镉废水排放, 结合图 5中因子2的空间分布, 中高值区基本覆盖该市大部分区域且包含农田区域, 所以因子2可以判断为一个自然母质、工业和人为耕作的混合源;因子3主要特征物为Cu和Zn, Cu和Zn的平均含量分别为广东省土壤背景值2倍和2.65倍, 金属冶炼加工、机械制造和电镀等工业会导致Cu和Zn的废水废气排放, 造成土壤Cu和Zn的富集[15, 37], 高值点附近存在电线厂、家具厂、金属制品厂等汇集的工业区, 结合图 5中因子3的空间分布贡献值偏高区域基本为工业园区域, 所以因子3可以判断为工业源;因子4的主要特征物为Hg(82.0%), Hg的含量平均值已超过广东省土壤背景值4.51倍, Hg主要来源于交通源, 交通污染会造成土壤中Hg的富集[32, 38], 结合图 5中因子4的贡献值空间分布, 偏高值区域较均匀分布于该市中心城区即车辆流动密集的地方, 所以因子4是交通源;因子5的主要特征物为Pb(80.9%), 铅蓄电池在生产过程中会产生含Pb废气、废水和固体废弃物. 早年该类企业的生产工艺、环保设施落后, 在一定程度上影响了周边土壤环境质量[39], 结合图 5中因子5的高值点分布, 该区域某铅蓄电池制造厂生产过程中对周边土壤中Pb的积累可能造成了显著影响, 所以因子5可以判断为铅蓄电池制造源;各因子贡献空间分布情况与该市土地利用背景基本符合.

图 3 土壤重金属PMF源成分谱 Fig. 3 Source composition spectrum of heavy metal PMF in soil

图 4 PMF模型各因子占比 Fig. 4 Proportion of various factors in PMF model

图 5 PMF模型土壤污染源各因子贡献值空间分布 Fig. 5 Spatial distribution of contribution values of various factors of soil pollution sources in the PMF model

2.3 基于蒙特卡洛的土壤健康风险评价 2.3.1 非致癌健康风险评估

某市土壤重金属非致癌健康风险评估结果如图 6表 5. 8种重金属非致癌健康风险指数(HQi)均值大小表现为:As > Cr > Pb > Cu > Cd > Hg > Ni > Zn;两类人群的非致癌健康风险表现为:儿童的非致癌健康风险指数 > 成人;不同接触方式的非致癌健康风险表现为:经口摄入 > 呼吸摄入 > 皮肤接触. 综合非致癌健康风险指数(HQ)表明:成人和儿童的综合非致癌健康风险指数均值分别为5.31E-02和6.54E-01, 最大值均小于1, 说明该市土壤重金属对人体带来的非致癌健康风险可以忽略不计.

图 6 蒙特卡洛的土壤重金属综合非致癌健康风险指数概率分布 Fig. 6 Monte Carlo probability distribution of comprehensive non⁃carcinogenic health risk index of soil heavy metals

表 5 蒙特卡洛的土壤重金属非致癌健康风险指数评估结果 Table 5 Monte Carlo assessment results of non⁃carcinogenic health risk index of heavy metals in soil

2.3.2 致癌健康风险

该市土壤重金属致癌健康风险评估结果如图 7表 6. 2种致癌重金属的致癌健康风险指数(CRi)平均值大小表现为:As > Cd;两种重金属对儿童的致癌健康风险均大于对成年, 与儿童的生理和行为等特征以及儿童对污染的敏感性更高有关[30];Cd对成年和儿童的CR的最大值范围均在1E-04~1E-06之间, 均值分别为4.70E-06和1.36E-05, 说明Cd对两类人群均存在一定程度的可耐受致癌健康风险;As对成年和儿童的CR的最大值范围均在1E-06~1E-04之间, 均值分别为1.54E-06和4.46E-06, 说明As对两类人群存在一定程度的可耐受致癌健康风险;Cd和As会对人体多个系统功能造成相对更大的危害[30, 40, 41].

图 7 蒙特卡洛土壤重金属综合致癌健康风险指数概率分布 Fig. 7 Probability distribution of comprehensive carcinogenic health risk index of heavy metals in soil in Monte Carlo

表 6 蒙特卡洛的土壤重金属致癌健康风险指数评估结果 Table 6 Monte Carlo assessment results of health risk index for heavy metal carcinogenesis in soil

2.3.3 敏感度分析

敏感度分析用于反映各参数对于风险结果的影响程度, 敏感度值越大表明对风险结果的影响程度越大[40]. 该市土壤重金属致癌健康风险(CR)和非致癌健康风险(HQ)的敏感度分析结果如图 8. 致癌和非致癌健康风险中敏感度值最大的均为As, 致癌健康风险中As对于成年和儿童的敏感度分别为64.7%和74.5%, 非致癌健康风险中As对于成年和儿童的敏感度分别为44.4%和55.5%, 加强管控As在土壤中的累积, 将有助于降低土壤对于人体的潜在健康风险, 在8种重金属中, As应作为主控元素, 其次是Cd、Pb和Cr;BW在致癌和非致癌健康风险中敏感度均为负值, 说明增长体重将有助于降低健康风险, 这与前人研究结果一致[26, 30, 42].

and non⁃carcinogenic health risks 图 8 致癌与非致癌健康风险敏感度分析 Fig. 8 Sensitivity analysis of carcinogenic

2.3.4 基于源导向的土壤重金属优先控制因子分析

基于儿童健康风险指数整体大于成年健康风险指数, 选用儿童的健康风险模型结合PMF模型解析结果对该市土壤主控因子进行分析, 结果如图 9. 因子3工业源是该市土壤重金属最主要的累积来源, 贡献率为47.6%, 主要由Pb、Cu和Zn为主, 王小莉等[43]对开封市的研究表明, Zn、Cu和Pb污染区主要集中在炼锌厂, 化肥厂等企业的周围. 尚婷婷等[44]对白银市城郊周边农田区域土壤的研究表明, 铅锌冶炼厂, 铜业公司等企业产生的“三废”进入周边环境, 可能造成周边土壤的重金属富集;引起人体健康危害最主要的因素是因子1农药化肥投入源, 主要由As贡献(贡献率为96%), 对重金属的贡献率为8.3%, 引起的致癌和非致癌健康风险的贡献率分别为66%和46%, 为最主要致癌和非致癌健康风险来源, 可能是由于不同金属元素对人体的影响存在差异, 模型中参考剂量和斜率因子取值不同. 有研究表明, As的毒性更高, 人体长期摄入As可导致肾、肝和膀胱等内脏器官的癌变[28, 44, 45], 因此须避免过度使用农药和化肥, 加强绿色种植的推广, 以控制和减轻对土壤的持续污染.

百分数表示贡献率, 曲线宽度对应贡献值的大小 图 9 土壤重金属、因子与健康风险的关系 Fig. 9 Relationship between heavy metals, factors, and health risks in soil

3 结论

(1)广东省某市土壤重金属ω(As)、ω(Hg)、ω(Cd)、ω(Pb)、ω(Cr)、ω(Cu)、ω(Ni)和ω(Zn)的平均值分别为18.16、0.43、1.46、68.57、98.34、64.19、26.53和257.32 mg·kg-1, 变异程度为中高度变异. 除Ni元素以外, 其余重金属元素土壤含量的平均值均已超过广东省土壤背景值, 且Cd和Zn的含量已超过国家二级标准, 重金属污染已较为严重.

(2)基于相关性分析和PMF模型结果表明, 该市重金属来源复杂, 主要受工业源影响, 自然母质、铅蓄电池制造、人为耕作、交通和农药化肥投入也对土壤重金属累积有不可忽视的影响, 污染因子主要为As、Cr和Pb.

(3)基于蒙特卡洛模拟的健康风险评估结果表明, 该市土壤重金属对儿童和成年均存在一定程度的可耐受致癌健康风险, 非致癌风险可以忽略不计;儿童的潜在健康风险大于成年, 主要暴露途径均为经口摄入.

(4)基于源导向的主控因子分析结果表明:农药化肥投入源和As应该作为该市土壤重金属影响健康风险的主控因子, 其次为混合源和Cr.

参考文献
[1] Xiang M T, Li Y, Yang J Y, et al. Heavy metal contamination risk assessment and correlation analysis of heavy metal contents in soil and crops[J]. Environmental Pollution, 2021, 278. DOI:10.1016/j.envpol.2021.116911
[2] Pekey H, Doğan G. Application of positive matrix factorisation for the source apportionment of heavy metals in sediments: A comparison with a previous factor analysis study[J]. Microchemical Journal, 2013, 106. DOI:10.1016/j.microc.2012.07.007
[3] Amuno S, Niyogi S, Amuno M, et al. Heavy metal bioaccumulation and histopathological alterations in wild Arctic hares (Lepus arcticus) inhabiting a former lead-zinc mine in the Canadian high Arctic: A preliminary study[J]. Science of the Total Environment, 2016, 556: 252-263. DOI:10.1016/j.scitotenv.2016.03.007
[4] Huang Y, Wang L Y, Wang W J, et al. Current status of agricultural soil pollution by heavy metals in China: A meta-analysis[J]. Science of the Total Environment, 2019, 651: 3034-3042. DOI:10.1016/j.scitotenv.2018.10.185
[5] Liu R M, Guo L J, Men C, et al. Spatial-temporal variation of heavy metals' sources in the surface sediments of the Yangtze River Estuary[J]. Marine Pollution Bulletin, 2019, 138: 526-533. DOI:10.1016/j.marpolbul.2018.12.010
[6] Wang C L, Zou X Q, Feng Z Y, et al. Distribution and transport of heavy metals in estuarine⁃inner shelf regions of the East China Sea[J]. Science of the Total Environment, 2018, 644: 298-305. DOI:10.1016/j.scitotenv.2018.06.383
[7] Lin L, Li C, Yang W J, et al. Spatial variations and periodic changes in heavy metals in surface water and sediments of the Three Gorges Reservoir, China[J]. Chemosphere, 2020, 240. DOI:10.1016/j.chemosphere.2019.124837
[8] Chai M W, Li R L, Tam N F Y, et al. Effects of mangrove plant species on accumulation of heavy metals in sediment in a heavily polluted mangrove swamp in Pearl River Estuary, China[J]. Environmental Geochemistry and Health, 2019, 41(1): 175-189. DOI:10.1007/s10653-018-0107-y
[9] 程俊伟, 蔡深文, 黄明琴. 贵州毛石钼矿区土壤及表层沉积物中重金属形态分析及生态风险评价[J]. 水土保持研究, 2021, 28(1): 353-359.
Cheng J W, Cai S W, Huang M Q. Analysis on speciation and risk assessment of heavy metals in soils and surface sediments of Maoshi molybdenum mining area in Guizhou[J]. Research of Soil and Water Conservation, 2021, 28(1): 353-359.
[10] 谢龙涛, 潘剑君, 白浩然, 等. 基于GIS的农田土壤重金属空间分布及污染评价——以南京市江宁区某乡镇为例[J]. 土壤学报, 2020, 57(2): 316-325.
Xie L T, Pan J J, Bai H R, et al. GIS-Based spatial distribution and risk assessment of heavy metals in farmland soils: a case study of a town of Jiangning, Nanjing[J]. Acta Pedologica Sinica, 2020, 57(2): 316-325.
[11] 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
[12] 张施阳. 基于GIS的上海市不同功能区土壤重金属污染评价及健康风险评估[J]. 环境工程技术学报, 2022, 12(4): 1226-1236.
Zhang S Y. Assessment of soil heavy metal pollution and health risk in different functional areas of Shanghai city based on GIS[J]. Journal of Environmental Engineering Technology, 2022, 12(4): 1226-1236.
[13] 张军, 董洁, 梁青芳, 等. 宝鸡城区土壤重金属空间分布特征及生态风险评价[J]. 干旱区资源与环境, 2018, 32(10): 100-106.
Zhang J, Dong J, Liang Q F, et al. Spatial distribution and risk assessment of heavy metals in urban soils of Baoji city[J]. Journal of Arid Land Resources and Environment, 2018, 32(10): 100-106.
[14] 夏子书, 白一茹, 王幼奇, 等. 基于PMF模型的宁南山区小流域土壤重金属空间分布及来源解析[J]. 环境科学, 2022, 43(1): 432-441.
Xia Z S, Bai Y R, Wang Y Q, et al. Spatial distribution and source analysis of soil heavy metals in a small watershed in the mountainous area of Southern Ningxia based on PMF model[J]. Environmental Science, 2022, 43(1): 432-441.
[15] 郑永立, 温汉辉, 蔡立梅, 等. 基于PMF模型的县域尺度土壤重金属来源分析及风险评价[J]. 环境科学, 2023, 44(9): 5242-5252.
Zheng Y L, Wen H H, Cai L M, et al. Source analysis and risk assessment of heavy metals in soil of county scale based on PMF model[J]. Environmental Science, 2023, 44(9): 5242-5252.
[16] 郭志娟, 周亚龙, 王乔林, 等. 雄安新区土壤重金属污染特征及健康风险[J]. 中国环境科学, 2021, 41(1): 431-441.
Guo Z J, Zhou Y L, Wang Q L, et al. Characteristics of soil heavy metal pollution and health risk in Xiong'an New District[J]. China Environmental Science, 2021, 41(1): 431-441.
[17] 梁敏静, 熊凡, 曾经文, 等. 广州郊区三类工业企业周边农田土壤重金属污染及生态风险评价[J]. 广东农业科学, 2021, 48(7): 103-110.
Liang M J, Xiong F, Zeng J W, et al. Heavy metal pollution and ecological risk assessment of farmland soil around three types of industrial enterprises in Guangzhou suburb[J]. Guangdong Agricultural Sciences, 2021, 48(7): 103-110.
[18] 吴俭, 邓一荣, 林龙勇, 等. 广州市建设用地土壤污染风险管控和修复现状、问题与对策[J]. 环境监测管理与技术, 2021, 33(3): 1-4, 14.
Wu J, Deng Y R, Lin L Y, et al. Risk control and soil remediation of construction land in Guangzhou: current situation, problems and countermeasures[J]. The Administration and Technique of Environmental Monitoring, 2021, 33(3): 1-4, 14.
[19] 王钊, 郭宇, 石晓龙. 广州市某工业区土壤重金属污染评价[J]. 环境科学与技术, 2020, 43(S2): 214-221.
Wang Z, Guo Y, Shi X L. Assessment of heavy metal pollution in soil of an industrial zone in Guangzhou[J]. Environmental Science & Technology, 2020, 43(S2): 214-221.
[20] 雷国建, 陈志良, 刘千钧, 等. 广州郊区土壤重金属污染程度及潜在生态危害评价[J]. 中国环境科学, 2013, 33(S1): 49-53.
Lei G J, Chen Z L, Liu Q J, et al. The assessments of polluted degree and potential ecological hazards of heavy metals in suburban soil of Guangzhou city[J]. China Environmental Science, 2013, 33(S1): 49-53.
[21] 何健, 陈海珍, 程世芬, 等. 广州主城区土壤重金属Pb污染特征及健康风险评价[J]. 现代预防医学, 2014, 41(2): 212-215, 218.
He J, Chen H Z, Cheng S F, et al. Evaluation of the characteristics and health risks of Pb pollution in soils in downtown Guangzhou[J]. Modern Preventive Medicine, 2014, 41(2): 212-215, 218.
[22] 黄银华. 广州市农业表层土壤重金属污染及其与景观格局的多尺度关系[D]. 广州: 广州大学, 2015.
[23] GB/T 22105.1-2008, 土壤质量总汞、总砷、总铅的测定原子荧光法第1部分: 土壤中总汞的测定[S].
GB/T 22105.1-2008, Soil quality—Analysis of total mercury, arsenic and lead contents in soils—Atomic fluorescence spectrometry—Part 1: Analysis of total mercury contents in soils[S].
[24] GB/T 22105.2-2008, 土壤质量总汞、总砷、总铅的测定原子荧光法第2部分: 土壤中总砷的测定[S].
GB/T 22105.2-2008, Soil quality—Analysis of total mercury, arsenic and lead contents in soils—Atomic fluorescence spectrometry—Part 2: Analysis of total arsenic contents in soils[S].
[25] 刘海, 魏伟, 黄健敏, 等. 长江流域(安徽段)土壤-作物系统重金属污染特征及健康风险评价[J]. 环境科学, 2023, 44(3): 1686-1697.
Liu H, Wei W, Huang J M, et al. Heavy metal pollution characteristics and health risk assessment of soil-crops system in Anhui section of the Yangtze River basin[J]. Environmental Science, 2023, 44(3): 1686-1697.
[26] HJ 25.3-2019, 建设用地土壤污染风险评估技术导则[S].
HJ 25.3-2019, Technical guidelines for risk assessment of soil contamination of land for construction[S].
[27] Yang S Y, Zhao J, Chang S X, et al. Status assessment and probabilistic health risk modeling of metals accumulation in agriculture soils across China: a synthesis[J]. Environment International, 2019, 128: 165-174. DOI:10.1016/j.envint.2019.04.044
[28] 马杰, 葛淼, 王胜蓝, 等. 基于源导向的农用地土壤重金属健康风险评估及优先控制因子分析[J]. 环境科学, 2024, 45(1): 396-406.
Ma J, Ge M, Wang S L, et al. Health risk assessment and priority control factors analysis of heavy metals in agricultural soils based on source-oriented[J]. Environmental Science, 2024, 45(1): 396-406.
[29] Chen R H, Chen H Y, Song L T, et al. Characterization and source apportionment of heavy metals in the sediments of Lake Tai (China) and its surrounding soils[J]. Science of the Total Environment, 2019, 694. DOI:10.1016/j.scitotenv.2019.133819
[30] Huang J L, Wu Y Y, Sun J X, et al. Health risk assessment of heavy metal(loid)s in park soils of the largest megacity in China by using Monte Carlo simulation coupled with positive matrix factorization model[J]. Journal of Hazardous Materials, 2021, 415. DOI:10.1016/j.jhazmat.2021.125629
[31] Huang Y N, Dang F, Li M, et al. Environmental and human health risks from metal exposures nearby a Pb-Zn-Ag mine, China[J]. Science of the Total Environment, 2020, 698. DOI:10.1016/j.scitotenv.2019.134326
[32] 沈智杰, 李杰芹, 李彩霞, 等. 基于APCS-MLR和PMF模型的赤泥堆场周边耕地土壤重金属污染源解析[J]. 环境科学, 2024, 45(2): 1058-1068.
Shen Z J, Li J Q, Li C X, et al. Pollution source apportionment of heavy metals in cultivated soil around a red Mud yard based on APCS-MLR and PMF model[J]. Environmental Science, 2024, 45(2): 1058-1068.
[33] 高梦绯, 郑顺安, 刘昌华, 等. 基于多因素融合的耕地土壤重金属污染风险评价[J]. 环境工程, 2023, 41(8): 233-241.
Gao M F, Zheng S A, Liu C H, et al. Risk assessment of heavy metal pollution in farmland soil of county based on multi-factors[J]. Environmental Engineering, 2023, 41(8): 233-241.
[34] 蒋玉莲, 余京, 王锐, 等. 渝东南典型地质高背景区土壤重金属来源解析及污染评价[J]. 环境科学, 2023, 44(7): 4017-4026.
Jiang Y L, Yu J, Wang R, et al. Source analysis and pollution assessment of soil heavy metals in typical geological high background area in Southeastern Chongqing[J]. Environmental Science, 2023, 44(7): 4017-4026.
[35] 徐福银, 胡艳燕, 包兵, 等. 重庆市不同功能区绿地土壤Hg分布特征及其污染评价[J]. 土壤通报, 2015, 46(5): 1253-1258.
Xu F Y, Hu Y Y, Bao B, et al. Distribution characters and pollution evaluation of Hg in urban soil of different functional area in Chongqing[J]. Chinese Journal of Soil Science, 2015, 46(5): 1253-1258.
[36] 陈丹青, 谢志宜, 张雅静, 等. 基于PCA/APCS和地统计学的广州市土壤重金属来源解析[J]. 生态环境学报, 2016, 25(6): 1014-1022.
Chen D Q, Xie Z Y, Zhang Y J, et al. Source Apportionment of soil heavy metals in Guangzhou based on the PCA/APCS model and geostatistics[J]. Ecology and Environmental Sciences, 2016, 25(6): 1014-1022.
[37] 方嘉, 何影, 黄乃涛, 等. 基于PMF模型的农田土壤重金属源暴露风险综合评价: 以浙江省某电子垃圾拆解区为例[J]. 环境科学, 2023, 44(7): 4027-4038.
Fang J, He Y, Huang N T, et al. Integrated analysis on source-exposure risk of heavy metals in farmland soil based on PMF model: A case study in the e-waste dismantling area in Zhejiang province[J]. Environmental Science, 2023, 44(7): 4027-4038.
[38] 任静华, 范健, 何培良, 等. 典型工矿企业周边农田土壤重金属污染评价及来源解析[J]. 环境污染与防治, 2023, 45(6): 855-861, 869.
Ren J H, Fan J, He P L, et al. Pollution characteristics and source apportionment of heavy metals in soils around a typical industrial enterprises[J]. Environmental Pollution & Control, 2023, 45(6): 855-861, 869.
[39] 黄昌辰, 温汉辉, 蔡立梅, 等. 广东揭阳市土壤铜的分布特征及其影响因素[J]. 现代地质, 2021, 35(1): 293-300.
Huang C C, Wen H H, Cai L M, et al. Distribution characteristics and influencing factors of soil Cu in Jieyang city, Guangdong province[J]. Geoscience, 2021, 35(1): 293-300.
[40] 马杰, 佘泽蕾, 王胜蓝, 等. 基于蒙特卡罗模拟的煤矸山周边农用地土壤重金属健康风险评估[J]. 环境科学, 2023, 44(10): 5666-5678.
Ma J, She Z L, Wang S L, et al. Health risk assessment of heavy metals in agricultural soils around the gangue heap of coal mine based on Monte Carlo simulation[J]. Environmental Science, 2023, 44(10): 5666-5678.
[41] 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
[42] 王静, 魏恒, 潘波. 中国农田土壤Cd累积分布特征及概率风险评价[J]. 环境科学, 2023, 44(7): 4006-4016.
Wang J, Wei H, Pan B. Accumulation characteristics and probabilistic risk assessment of Cd in agricultural soils across China[J]. Environmental Science, 2023, 44(7): 4006-4016.
[43] 王小莉, 陈志凡, 魏张东, 等. 开封市城乡交错区农田土壤重金属污染及潜在生态风险评价[J]. 环境化学, 2018, 37(3): 513-522.
Wang X L, Chen Z F, Wei Z D, et al. Heavy metal pollution and potential ecological risk assessment in agricultural soils located in the peri-urban area of Kaifeng City[J]. Environmental Chemistry, 2018, 37(3): 513-522.
[44] 尚婷婷, 张亚群, 周静, 等. 白银市城郊周边农田土壤重金属污染特征及来源解析[J]. 环境生态学, 2023, 5(4): 19-26.
Shang T T, Zhang Y Q, Zhou J, et al. Pollution characteristics and source apportionment of heavy metals in farmland soil, Baiyin suburb[J]. Environmental Ecology, 2023, 5(4): 19-26.
[45] Smith A H, Goycolea M, Haque R, et al. Marked increase in bladder and lung cancer mortality in a region of Northern Chile due to arsenic in drinking water[J]. American Journal of Epidemiology, 1998, 147(7): 660-669. DOI:10.1093/oxfordjournals.aje.a009507