环境科学  2024, Vol. 45 Issue (12): 7166-7176   PDF    
基于源导向的锰矿尾矿库周边土壤重金属风险评估
马杰1,2, 王胜蓝1,2, 秦启荧1, 文川勇1, 李名升2, 封雪2     
1. 重庆市生态环境监测中心, 重庆 401147;
2. 中国环境监测总站, 北京 100012
摘要: 以重庆市锰矿尾矿库周边土壤为研究对象, 运用正定矩阵因子分解(PMF)模型对土壤重金属开展源解析, 结合潜在生态风险指数法和健康风险评估模型, 定量探讨不同污染源影响下的土壤潜在生态风险和健康风险, 并确定优先管控要素.结果表明:研究区土壤除Cr外, Mn、Cd、As、Pb、Cu、Zn和Ni含量均值均高于背景值. PMF模型源解析结果表明, 研究区土壤受矿业源、自然源、农业源和工业源影响, 贡献率分别为23.9%、30.2%、18.8%和27.1%.潜在生态风险评估结果表明, 研究区土壤以“强”生态风险为主, 占比为77.8%.健康风险模型评估结果表明, 研究区土壤重金属对儿童和成人存在人体可接受的致癌风险, 非致癌风险可忽略.污染因子、污染源、健康风险和潜在生态风险综合分析表明, 农业源为研究区土壤首要管控源, Cd为首要管控因子;矿业源和工业源为研究区土壤次要管控源, Mn和As为次要管控因子.
关键词: 土壤      重金属      PMF模型      健康风险评估      潜在生态风险      锰(Mn)     
Risk Assessment of Heavy Metals in Soil Surrounding Manganese Tailings Pond Based on Source-oriented Analysis
MA Jie1,2 , WANG Sheng-lan1,2 , QIN Qi-ying1 , WEN Chuan-yong1 , LI Ming-sheng2 , FENG Xue2     
1. Chongqing Ecological and Environmental Monitoring Center, Chongqing 401147, China;
2. China National Environmental Monitoring Centre, Beijing 100012, China
Abstract: To analyze the source apportionment, potential ecological risk, and health risk of heavy metals in soils surrounding a manganese tailings pond in Chongqing, a positive matrix factorization (PMF) model, potential ecological risk index, and health risk assessment model were used. Further, all three models were combined to explore the risks of heavy metals in soils by different pollution sources to determine the priority control factors. The results showed that except for the Cr concentration, the average values of Mn, Cd, As, Pb, Cu, Zn, and Ni concentration were higher than their corresponding background values. Using the PMF model analysis, mining, natural, agricultural, and industrial sources were identified as the determinants for the accumulation of heavy metals in soils, with contribution rates of 23.9%, 30.2%, 18.8%, and 27.1%, respectively. Using potential ecological risk index analysis, the ecological risk was predominantly categorized as "strong" risk, with the rate of 77.8%. The health risk assessment model analysis revealed that, the carcinogenic risks of adults and children were tolerable, whereas the non-carcinogenic risks were acceptable. According to the analysis results of the relationship among heavy metals, pollution sources, potential ecological risk, and health risks showed that agricultural sources were identified as the priority control sources, with Cd as the primary control factor. Mining sources and industrial sources were identified as secondary control sources, with Mn and As as the secondary control factors.
Key words: soil      heavy metals      PMF model      health risk assessment      potential ecological risk      manganese(Mn)     

锰矿作为我国重要的战略矿产资源, 是钢铁工业的基本原料, 也是电池工业、陶瓷工业、化学工业等领域的重要原料[1, 2].但锰矿在开采、洗选、冶炼和加工过程中会产生大量尾矿或废渣, 由于技术落后、环保意识薄弱和管理不到位等原因, 这些尾矿或废渣在地表径流和大气扩散等作用下, 易对周边环境造成污染[3 ~ 5].有研究表明, 锰渣场或尾矿库周边土壤普遍存在不同程度的重金属污染[6 ~ 8], 不仅会影响农作物生长和农产品减产, 还会威胁人体健康, 诱发神经系统疾病, 肝、肾、脑等器官障碍和损伤[9 ~ 12].因此, 开展锰渣场或尾矿库周边土壤重金属风险评估具有重要意义.

鉴于土壤重金属污染成因较多, 除受锰渣场或尾矿库影响外, 还可能受工业排放、污水灌溉、农药化肥使用和成土母岩风化等因素影响[13 ~ 16].因此, 掌握不同源对土壤污染的贡献程度, 才能为后续源头管控提供支撑.当前, 针对土壤污染源解析的研究较多, 其中正定矩阵因子分解(PMF)模型能对因子分解矩阵进行非负约束, 使解析结果更加准确而被广泛应用[7, 17 ~ 19].但将土壤污染源解析和土壤风险评估关联, 定量探讨它们之间关系的研究较少.

重庆市秀山县、贵州省铜仁市松桃县和湖南省湘西州花垣县因地理位置相邻, 锰矿资源丰富, 被称为中国“锰三角”[20].长期以来受锰矿业发展影响, 周边土壤存在较大的污染风险[21].因此, 本文以重庆市秀山县某锰矿尾矿库为例, 测定尾矿库周边土壤重金属含量, 运用潜在生态风险指数法、健康风险评估模型和PMF模型分别开展生态风险评估、人体健康风险评估和污染源解析, 定量探讨三者之间的关系, 综合确定优先管控因子, 以期为研究区土壤污染防治提供科学支撑.

1 材料与方法 1.1 区域概况

研究区坐落于重庆市秀山县西南部, 海拔高度360~460 m.南部为锰矿尾矿库, 占地面积约12 hm2, 处于山谷之间, 堆存周边企业产生的电解锰渣, 堆存量约200 m3, 堆存近10 a, 尾矿库在2022年扩容, 完善了挡渣墙、防渗和污水处理等设施.尾矿库西侧紧邻锰加工厂, 占地面积约6 hm2, 2010年建厂生产, 现已停产, 内部设备设施已拆除, 研究区外有一对外主干道与尾矿库和锰加工厂相连, 交通便利.尾矿库北侧槽谷地带为成片耕地和居民点, 占地面积分别为20 hm2和15 hm2, 耕地以水田为主, 槽谷中间有一条河流自南向北流淌, 槽谷两侧为荒山和林地, 坡度较大(图 1).

图 1 研究区区位和土壤点位分布示意 Fig. 1 Location of the study area and distribution of soil sampling sites

1.2 样品采集与测定

为使土壤点位布设具有代表性, 根据现场情况, 按照《土壤环境监测技术规范》(HJ/T 166-2004)要求, 充分考虑污染途径、污染范围、污染类型等因素, 在仅靠尾矿库和易受尾矿库地表径流影响的耕作田块上布设土壤监测点位[19], 每个相邻点位相距约100 m, 共布设18个土壤采样点位(图 1).现场用双对角线5点混合法采集0~20 cm的表层土壤, 带回实验室经自然风干后, 按照表 1中的方法分析测试.

表 1 土壤分析方法1) Table 1 Methods of soil analysis

1.3 研究方法 1.3.1 潜在生态风险指数法

潜在生态风险指数法是土壤常用的评价方法之一, 综合考虑了重金属性质、生物毒性和生态效应等因素[25, 26], 公式如(1)~(2):

(1)
(2)

式中, RI为综合潜在生态风险指数;;Tri为第i种土壤重金属毒性响应系数[7, 26]Eri为第i种土壤重金属潜在生态风险指数;Si为第i种重庆市土壤背景值[27].

根据马建华等[28]对潜在生态风险指数的优化分级调整原则, 本研究取8项重金属中的最大毒性系数(Cd=30)为Eri第一级界限值, 其他Eri风险级别的界限值依次乘2, RI第一级界限值则是8项重金属的毒性系数总和取整, 其他RI风险级别的界限值依次乘2, 具体见表 2.

表 2 优化后的潜在生态风险评价标准 Table 2 Optimized grading standards of potential ecological risk

1.3.2 土壤健康风险评估模型

土壤健康风险评估模型采用美国环境保护局推荐的模型, 也是当前最常用的健康风险模型[29], 公式如(3)~(7):

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

式中, ADDi,ingest、ADDi,dermal和ADDi,inhal分别为第i种重金属经口、皮肤和呼吸日均摄入量;Ci为第i种土壤重金属的含量;HQ和HI分别为单一和综合非致癌风险指数;CR和TCR分别为单一和综合致癌风险指数;RfDi为第i种土壤重金属的参考剂量;SFi为第i种土壤重金属的斜率因子.有关参数和取值见表 3表 4.

表 3 土壤健康风险模型参数取值1) Table 3 Parameters and values in health risk assessment model

表 4 参考剂量和斜率因子取值1) Table 4 Reference dose and slope factor values

此外, 为加强模型的准确性, 本研究引入蒙特卡罗不确定性分析, 利用Oracle Crystal Ball 11.1.2.4处理数据, 迭代次数和置信水平分别设置为10 000和95%.

1.3.3 正定矩阵因子分解(PMF)模型

PMF模型是将重金属含量矩阵分解为因子贡献和因子残差矩阵[19], 如公式(8):

(8)

式中, Xij为第i个样本中j种重金属的含量, giki种重金属在源k的含量, fkj为源k对第i个样本的贡献, eij为第i个样本中j种重金属的残差.

通过对数据不确定性矩阵U进行迭代计算, 不断分解矩阵X, 选择最优的矩阵GF, 最终得到最小目标函数Q.公式如(9)和(10):

(9)
(10)

式中, MDL为方法检出限(表 1);uij为第i个样本中j种重金属含量的不确定性值;RSD为重金属含量的相对标准偏差;c为实测含量.

1.3.4 源导向的潜在生态风险和人体健康风险

基于PMF模型获取的不同污染源上各重金属贡献率, 并结合研究区潜在生态风险评估结果, 计算不同源对潜在生态风险的贡献率, 如公式(11)和(12):

(11)
(12)

式中, RIj为第j类源的潜在生态风险;RIi为第i种重金属的潜在生态风险;Fij为第i种重金属在第j类源的贡献率;Dj, RI为第j类源的潜在生态风险贡献率.

基于PMF模型获取的不同污染源上各重金属贡献率, 并结合研究区健康风险评估结果, 计算不同源对人体健康风险的贡献率[36], 如公式(13)~(16):

(13)
(14)
(15)
(16)

式中, HQj为第j类源的非致癌风险;Fij为第i种重金属在第j类源的贡献率;HQi为第i种重金属的非致癌风险;Dj, HQ为第j类源的非致癌风险贡献率;CRj为第j类源的致癌风险;CRi为第i种重金属的致癌风险;Dj, CR为第j类源的致癌风险贡献率.

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

研究区土壤重金属与国内土壤环境标准相比(表 5), 研究区仅Mn、Cd和As存在超标准限值的情况, 均值分别为1 162、0.37和24.3mg·kg-1, 超标占比分别为50.0%、38.9%和11.1%, Pb、Cr、Cu、Zn和Ni无超标准限值的情况, 均值分别为44.0、76.4、33.9、106和37.9mg·kg-1.与重庆土壤背景值相比, 研究区As和Cd含量均值分别是背景值的4.86倍和3.40倍, 远高于背景值, Mn、Pb、Cu、Ni和Zn含量均值是背景值的1.18~1.89倍, Cr含量均值则是背景值的0.96倍, 略低于背景值.

表 5 研究区土壤重金属含量统计 Table 5 Statistical characteristics of heavy metal contents in soil of the study area

变异系数常用于反映土壤受人为活动影响的程度[37].研究区Cd变异系数最大, 为40.7%, 其次是Mn, 与Cd变异系数接近, 为40.3%, 然后是As变异系数33.9%, Pb、Cu、Ni和Zn变异系数在21.5%~31.6%, 最后是Cr, 变异系数最小, 为14.8%.说明研究区Cd、Mn和As受人为活动影响较大, Cr受人为活动影响最小.综上, 说明Mn、Cd和As是研究区土壤的主要污染物.

1)数据来源于《农用地土壤环境风险评价技术规定(试行)》(环办土壤函〔2018〕1479号);2)数据来源于GB 15618-2018中的筛选值

从空间分布看(图 2), Mn含量在尾矿库周边最大, 并沿地表径流方向逐步降低, 总体呈南高北低, 说明尾矿库堆存可能导致周边土壤Mn含量升高.Cd含量在水田区域和尾矿库东侧相对较大, 总体呈东高西低, 说明尾矿库堆存可能导致周边和下游土壤Cd含量升高.As含量在原锰加工厂和尾矿库周边均较大, 且随着距离的增加呈递减趋势, 总体呈南高北低, 说明原锰加工厂生产活动和尾矿库长期堆存可能导致As含量升高.Pb、Cr、Cu、Zn和Ni含量分布相似, 在原锰加工厂附近出现最大值, 而其他区域含量总体较小, 说明原锰加工厂生产活动可能导致Pb、Cr、Cu、Zn和Ni含量升高, 但影响范围有限.综上, Mn、Cd、As、Pb、Cr、Cu、Zn和Ni含量在一定程度上均受到锰加工厂或尾矿库影响, 这与有关研究结果基本一致[3, 6 ~ 8].

图 2 土壤重金属含量空间分布 Fig. 2 Spatial distributions of heavy metal contents in soil

2.2 土壤重金属源解析

相关性分析能在一定程度上反映土壤重金属之间的同源性[38].研究区土壤重金属相关性结果表明(表 6), Mn与其他7项重金属均无显著相关性, 说明Mn的主要来源可能与其他重金属不同;Pb、Cr、Cu、Zn和Ni之间均存在显著正相关关系, 除Cr和Pb相关性系数为0.51(P < 0.05)相对偏低外, 其他元素之间相关性系数在0.65~0.91(P < 0.01)之间, 说明Pb、Cr、Cu、Zn和Ni可能存在同源性.此外, Cu和Ni分别与As均存在显著正相关关系, 相关性系数分别为0.56(P < 0.05)和0.50(P < 0.05), 说明Cu、Ni和As也可能存在一定同源性.因此, 推断Cu、Ni来源相对复杂, 受不同污染源的复合影响可能性较大.

表 6 研究区土壤重金属含量相关性分析1) Table 6 Correlation between the heavy metal contents in soil of the study area

利用EPA PMF 5.0进行源解析, 当因子数设置为4个时, 模型效果最佳.结果如[图 3(a)3(b)]所示, 因子1、因子2、因子3和因子4占比分别为23.9%、30.2%、18.8%和27.1%.其中, 因子1中Mn贡献率最高, 为81.4%, 远高于其他7项重金属.从空间分布看[图 2(a)], 距尾矿库越近, Mn含量越大, 并沿地表径流方向逐步降低, 这与张永江等研究结果基本一致[3].此外, 因子1对Cd、As、Pb、Cr、Cu、Zn和Ni等其他7项重金属也存在一定贡献率, 分别为29.4%、17.7%、16.1%、19.3%、0.1%、17.2%和9.8%, 这可能与矿渣伴生其他重金属有关[7, 39, 40].因此, 推测因子1为矿业源, 主要受尾矿长期堆存影响.

(a)各因子贡献率, (b)各重金属贡献率 图 3 研究区土壤重金属贡献 Fig. 3 Contribution ratios of heavy metals in soil of the study area

因子2中Pb、Cr、Cu、Zn和Ni贡献率最高, 分别为45.3%、44.3%、45.8%、46.5%和41.8%, 结合Pb、Cr、Cu、Zn和Ni具有显著正相关关系, 说明具有一定同源性.从变异系数看, 该5项重金属变异系数低于Mn、Cd和As, 尤其是Cr变异系数最低, 为14.8%, 且Cr含量均值与重庆土壤背景值接近, 是背景值的0.96倍, 说明Cr相对稳定, 受非人为活动影响的可能性较大.有研究表明, 土壤Cr主要受成土母质等自然因素影响[33, 36, 41, 42].因此, 推测因子2为自然源.

因子3中Cd贡献率最高, 为56.6%, 从空间分布看[图 2(b)], Cd含量在水田区域和尾矿库东侧相对较大.大量研究表明, 农田中长期施用化肥和农药等是导致土壤重金属含量升高的重要因素之一[15, 43, 44], 如磷肥原料中含有Cd元素, 长期施用会导致土壤中Cd含量增加[45].Cheraghi等[46]研究发现, 氮磷钾复合肥的施用对土壤中Cd含量的增加具有显著贡献.Shi等[47]研究发现, 氮肥的施用是土壤中Cd含量升高的重要因素.通过对研究区的实地调查, 研究区长期种植水稻, 以施用氮肥和磷肥为主, 推测Cd含量在水田区域含量相对较高与农业生产活动有关.此外, Cd含量在尾矿库东侧也相对较大, 如前所述因子1受尾矿长期堆存影响, Cd在因子1上的贡献率为29.4%, 仅次于Mn.因此, 推测因子3为农业源, 主要受农药化肥施用等影响.

因子4中As贡献率最高, 为55.2%, 从空间分布看[图 2(c)~2(h)], As含量在原锰加工厂和尾矿库周边较大, Pb、Cr、Cu、Zn和Ni含量最大值也在原锰加工厂附近, 在因子4中的贡献率也相对较高, 分别为20.5%、22.4%、44.1%、23.7%和35.9%.有研究表明, 工业废水、废渣和废气排放可能导致重金属在土壤中累积[26, 48, 49], 而电解金属锰行业在生产过程中会排放大量废水、废渣和废气, 鉴于早期企业环境意识薄弱和管理不到位等原因, 均可能造成周边土壤环境污染[3, 4, 50, 51].因此, 推测因子4为工业源, 主要受原锰加工厂生产活动影响.

2.3 土壤重金属风险评估 2.3.1 土壤重金属潜在生态风险评估

土壤重金属综合潜在生态风险评估结果表明[图 4(a)], 研究区土壤以“强”生态风险为主, 占比为77.8%;其次是“很强”和“中等”生态风险, 占比均为11.1%.各重金属潜在生态风险评估结果表明[图 4(b)], Cd的潜在生态风险最高, 以“强”生态风险为主, 占比为61.1%, “很强”和“中等”生态风险占比分别为27.8%和11.1%;其次是As, 以“中等”生态风险为主, 占比为66.6%, “强”和“轻微”生态风险占比均为16.7%;最后是Mn、Pb、Cu、Zn和Ni, 潜在生态风险评估结果均为“轻微”生态风险.

(a)综合潜在风险评估结果, (b)各重金属潜在风险评估结果 图 4 研究区土壤潜在生态风险 Fig. 4 Potential ecological risk in soil of the study area

2.3.2 土壤重金属健康风险评估

各重金属对儿童和成人的非致癌风险(HQ)均值表现为(表 7):1.00E+0 > As > Mn > Cr > Pb > Ni > Cu > Cd > Zn, 说明各重金属对儿童和成人的非致癌风险可忽略.综合非致癌风险指数(HI)表明, 成人HI最大值为1.24E-1(HI≤1.00E+0)说明研究区成人非致癌风险可忽略.儿童HI最大值虽然为1.03E+0(HI > 1.00E+0), 但在95%分位值为9.90E-1(HI≤1.00E+0), 说明整体上研究区儿童非致癌风险也可忽略.

表 7 研究区土壤非致癌风险 Table 7 Non-carcinogenic health risk in soil of the study area

各重金属对儿童和成人的致癌风险(CR)均值表现为(表 8):1.00E-4 > As > 1.00E-6 > Cr > Pb > Cd > Ni, 说明As是引起儿童和成人致癌风险的首要元素.综合致癌风险指数(TCR)表明, 儿童和成人TCR最小值分别为8.58E-6和3.15E-6(1.00E-6 < TCR≤1.00E-4), 最大值分别为3.06E-5和1.07E-5(1.00E-6 < TCR≤1.00E-4), 说明研究区土壤重金属对儿童和成人均存在致癌风险, 但致癌风险为人体可接受水平.

表 8 研究区土壤致癌风险 Table 8 Carcinogenic health risk in soil of the study area

综上, 非致癌和致癌风险均表明As是影响人体健康的主要元素, 其原因在于As对人体的毒性更高, 人体长期摄入As可导致肝脏、肾脏和心肌受损, 出现多发性神经炎的症状等, 因此, 模型中As的参考剂量(RfDi)和斜率因子(SFi)取值较大[18, 52].此外, 儿童非致癌和致癌风险风险高于成人, 这可能和儿童的生理和行为特征等有关, 使其对污染物敏感性更高, 与已有研究结果一致[18, 52].

2.4 基于源导向的土壤重金属风险评估

因上述研究结果表明儿童健康风险高于成人, 故以儿童健康风险评估结果为例, 研究区农用地土壤重金属、污染源与健康风险关系如图 5所示.因子4(工业源)为健康风险的首要污染源, 对致癌和非致癌风险的贡献率最高, 分别为53.8%和36.9%, 污染因子以As为主.其次是因子1(矿业源), 是健康风险的次要污染源, 对致癌和非致癌风险的贡献率分别为17.8%和30.2%, 污染因子以Mn为主.源解析结果表明, 因子4(工业源)和因子1(矿业源)分别与原锰工厂生产和尾矿长期堆存有关, 说明涉锰污染源(因子1和因子4)是引起研究区健康风险的主要因素, 其中致癌风险的贡献率为71.6%, 非致癌风险的贡献率为67.1%, 污染因子以As和Mn为主.

图 5 土壤重金属、污染源和健康风险的关系 Fig. 5 Relationship among heavy metals, pollution sources, and health risks in soil

研究区土壤重金属、污染源与潜在生态风险关系如图 6所示.影响潜在生态风险的污染源表现为:因子3(农业源) > 因子1(矿业源) > 因子4(工业源) > 因子2(自然源).农业源是土壤潜在生态风险的首要污染源, 对潜在生态风险的贡献率为39.0%, 污染因子以Cd为主.有研究表明, 土壤中Cd的活性和迁移能力较强, 易被植物根系吸收, 并转移至植物各个器官(茎、叶和果实等), 积累到一定浓度后对植物产生毒害[53 ~ 55], 且农产品存在污染, 经食用后其污染物更易于被人体吸收, 远超通过手口、皮肤和呼吸等3种暴露途径直接摄入对人体健康的影响[52, 56, 57].因此, 应加强研究区农业污染源的管控, 如采取土壤改良、种植结构调整、使用低毒性有机肥代替化学肥料等措施, 保障农用地安全利用[54, 58].此外, 因子1(矿业源)和因子4(工业源)贡献率相近, 对潜在生态风险的贡献率分别为24.2%和23.0%, 故涉锰污染源应作为土壤潜在生态风险的次要污染源, 污染因子以Mn和As为主.

图 6 土壤重金属、污染源和潜在生态风险的关系 Fig. 6 Relationship among heavy metals, pollution sources, and potential ecological risks in soil

基于源导向的潜在生态风险和健康风险评估结果存在一定差异, 结合现场调查情况看, 原锰加工厂已停产, 生产设施设备已全部拆除, 说明不会对研究区土壤造成新的污染.同时尾矿库在2022年扩容, 完善了挡渣墙、防渗和污水处理等设施, 说明切断了尾矿库污染周边土壤的路径.总体而言, 当前研究区涉锰污染源已得到初步管控.此外, 土壤健康风险评估(表 7表 8)和潜在生态风险评估(图 4)结果表明, 涉锰污染源引起的人体致癌风险为可接受水平, 非致癌风险均可忽略, 而土壤潜在生态风险以“强”为主, 说明研究区土壤对人体的健康风险较小, 应重点关注土壤潜在生态风险的影响, 马杰等[18, 36]在之前研究中仅考虑了污染源与健康风险之间的关系[36], 本研究将潜在生态风险考虑其中, 能更有利于管理者综合决策.综上, 根据污染源分别与人体健康风险和潜在生态风险分析可综合确定, 农业源为研究区首要管控源, Cd为首要管控因子;矿业源和工业源为研究区次要管控源, Mn和As为次要管控因子.

3 结论

(1)除Cr外, 研究区Mn、Cd、As、Pb、Cu、Zn和Ni含量均值高于重庆土壤背景值, 空间分布上看, 8项重金属在一定程度上均受到锰加工厂或尾矿库影响.

(2)源解析结果表明, 研究区矿业源贡献率为23.9%, 受锰尾矿长期堆存影响, 污染因子以Mn为主;自然源贡献率为30.2%, 受成土母岩风化等影响, 污染因子以Pb、Cr、Cu、Zn和Ni为主;农业源贡献率为18.8%, 受农药化肥施用等影响, 污染因子以Cd为主;工业源贡献率为27.1%, 受原锰加工厂生产影响, 污染因子以As为主.

(3)潜在生态风险评估结果表明, 研究区土壤以“强”生态风险为主, 占比为77.8%, 污染因子以Cd和As为主.健康风险评估结果表明, 儿童和成人存在人体可接受的致癌风险, 但非致癌风险可忽略, 污染因子以As为主.

(4)综合污染因子、污染源、健康风险和潜在生态风险分析, 农业源为研究区首要管控源, Cd为首要管控因子;矿业源和工业源为研究区次要管控源, Mn和As为次要管控因子.

参考文献
[1] 任辉, 刘敏, 王自国, 等. 我国锰矿资源及产业链安全保障问题研究[J]. 中国工程科学, 2022, 24(3): 20-28.
Ren H, Liu M, Wang Z G, et al. Security of manganese resources and industrial chain in China[J]. Strategic Study of CAE, 2022, 24(3): 20-28.
[2] 程湘, 胡鹏, 张海坤, 等. 锰矿主要类型、分布特点及开发现状[J]. 中国地质, 2021, 48(1): 102-119.
Cheng X, Hu P, Zhang H K, et al. The main types, distribution and current development of manganese ore deposits[J]. Geology in China, 2021, 48(1): 102-119.
[3] 张永江, 李璐, 马双, 等. 典型锰矿区周边农田土壤重金属污染风险评价及其来源分析[J]. 有色金属(冶炼部分), 2023(10): 138-148.
Zhang Y J, Li L, Ma S, et al. Pollution evaluation and source analysis of heavy metals in farmland soils around the typical manganese mining Area[J]. Nonferrous Metals (Extractive Metallurgy), 2023(10): 138-148.
[4] 陈波, 张浩哲. 建设用地土壤锰风险控制值探讨——以重庆市某锰污染地块为例[J]. 土壤通报, 2020, 51(6): 1482-1488.
Chen B, Zhang H Z. Risk control values of soil manganese in construction lands: a case study of a manganese-contaminated site in Chongqing[J]. Chinese Journal of Soil Science, 2020, 51(6): 1482-1488.
[5] 罗洋, 刘方, 任军, 等. 锰渣堆场优势植物对Mn的富集特点及作用机制[J]. 水土保持通报, 2021, 41(4): 46-53, 60.
Luo Y, Liu F, Ren J, et al. Enrichment characteristics and mechanism of Mn by dominant plants in manganese residue disposal sites[J]. Bulletin of Soil and Water Conservation, 2021, 41(4): 46-53, 60.
[6] Luo X, Ren B Z, Hursthouse A S, et al. Potentially toxic elements (PTEs) in crops, soil, and water near Xiangtan manganese mine, China: potential risk to health in the foodchain[J]. Environmental Geochemistry and Health, 2020, 42(7): 1965-1976. DOI:10.1007/s10653-019-00454-9
[7] 余高, 陈芬, 张晓东, 等. 锰矿区周边农田土壤重金属污染特征、来源解析及风险评价[J]. 环境科学, 2023, 44(8): 4416-4428.
Yu G, Chen F, Zhang X D, et al. Pollution characteristics, source analysis, and risk assessment of heavy metals in the surrounding farmlands of manganese mining area[J]. Environmental Science, 2023, 44(8): 4416-4428.
[8] Li Y, Liu K H, Zhu J, et al. Manganese accumulation and plant physiology behavior of Camellia oleifera in response to different levels of nitrogen fertilization[J]. Ecotoxicology and Environmental Safety, 2019, 184. DOI:10.1016/j.ecoenv.2019.109603
[9] Kováčik J, Štěrbová D, Babula P, et al. Toxicity of naturally-contaminated manganese soil to selected crops[J]. Journal of Agricultural and Food Chemistry, 2014, 62(29): 7287-7296. DOI:10.1021/jf5010176
[10] 陆凤, 陈淼, 陈兰兰. 贵州松桃某电解锰企业锰渣重金属污染特征及对植物生长的毒性效应[J]. 科学技术与工程, 2018, 18(5): 124-129.
Lu F, Chen M, Chen L L. Characteristics of heavy metal pollution in manganese residues of electrolytic manganese in Songtao, Guizhou and its toxic effects on plant growth[J]. Science Technology and Engineering, 2018, 18(5): 124-129.
[11] Miah M R, Ijomone O M, Okoh C O A, et al. The effects of manganese overexposure on brain health[J]. Neurochemistry International, 2020, 135. DOI:10.1016/j.neuint.2020.104688
[12] Boudissa S M, Lambert J, Müller C, et al. Manganese concentrations in the soil and air in the vicinity of a closed manganese alloy production plant[J]. Science of the Total Environment, 2006, 361(1-3): 67-72. DOI:10.1016/j.scitotenv.2005.05.001
[13] Liu C X, Wang C X, Li Y C, et al. Spatiotemporal differentiation and geographic detection mechanism of ecological security in Chongqing, China[J]. Global Ecology and Conservation, 2022, 35. DOI:10.1016/j.gecco.2022.e02072
[14] Ma J, Shen Z J, Wang S L, et al. Source apportionment of heavy metals in soils around a coal Gangue heap with the APCS-MLR and PMF receptor models in Chongqing, southwest China[J]. Journal of Mountain Science, 2023, 20(4): 1061-1073. DOI:10.1007/s11629-022-7819-2
[15] Atafar Z, Mesdaghinia A, Nouri J, et al. Effect of fertilizer application on soil heavy metal concentration[J]. Environmental Monitoring and Assessment, 2010, 160(1-4): 83-89. DOI:10.1007/s10661-008-0659-x
[16] Chaoua S, Boussaa S, El Gharmali A, et al. Impact of irrigation with wastewater on accumulation of heavy metals in soil and crops in the region of Marrakech in Morocco[J]. Journal of the Saudi Society of Agricultural Sciences, 2019, 18(4): 429-436. DOI:10.1016/j.jssas.2018.02.003
[17] Cao J F, Li C F, Zhang L X, et al. Source apportionment of potentially toxic elements in soils using APCS/MLR, PMF and geostatistics in a typical industrial and mining city in Eastern China[J]. PLoS One, 2020, 15(9). DOI:10.1371/journal.pone.0238513
[18] 马杰, 葛淼, 王胜蓝, 等. 基于源导向的农用地土壤重金属健康风险评估及优先控制因子分析[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.
[19] 马杰, 沈智杰, 张萍萍, 等. 基于APCS-MLR和PMF模型的煤矸山周边耕地土壤重金属污染特征及源解析[J]. 环境科学, 2023, 44(4): 2192-2203.
Ma J, Shen Z J, Zhang P P, et al. Pollution characteristics and source apportionment of heavy metals in farmland soils around the gangue heap of coal mine based on APCS-MLR and PMF receptor model[J]. Environmental Science, 2023, 44(4): 2192-2203.
[20] 祝滔, 江长胜, 郝庆菊, 等. 重庆秀山锰矿区土壤和植物锰污染调查与评价[J]. 环境科学与技术, 2012, 35(9): 167-172.
Zhu T, Jiang C S, Hao Q J, et al. Investigation and assessment of contaminated soils and plants by Mn in manganese mining area in Xiushan autonomous county of Chongqing[J]. Environmental Science & Technology, 2012, 35(9): 167-172.
[21] Yang Q Q, Li Z Y, Lu X N, et al. A review of soil heavy metal pollution from industrial and agricultural regions in China: pollution and risk assessment[J]. Science of the Total Environment, 2018, 642: 690-700. DOI:10.1016/j.scitotenv.2018.06.068
[22] GB/T 17141-1997, 土壤质量铅、镉的测定石墨炉原子吸收分光光度法[S].
[23] GB/T 22105.2-2008, 土壤质量总汞、总砷、总铅的测定原子荧光法第2部分: 土壤中总砷的测定[S].
[24] HJ 780-2015, 土壤和沉积物无机元素的测定波长色散X射线荧光光谱法[S].
[25] 吴梅, 刘属灵, 袁余洋, 等. 土壤重金属潜在生态风险指数法优化研究——以重庆市城口县为例[J]. 土壤通报, 2023, 54(2): 473-480.
Wu M, Liu S L, Yuan Y Y, et al. Optimization of potential ecological risk index method for soil heavy metals——a case study of Chengkou County, Chongqing City[J]. Chinese Journal of Soil Science, 2023, 54(2): 473-480.
[26] Yang Y Y, Zhang J X, Xiao X, et al. Speciation and potential ecological risk of heavy metals in soils from overlapped areas of farmland and coal resources in northern Xuzhou, China[J]. Bulletin of Environmental Contamination and Toxicology, 2021, 107(6): 1053-1058. DOI:10.1007/s00128-021-03148-0
[27] 成杭新, 李括, 李敏, 等. 中国城市土壤化学元素的背景值与基准值[J]. 地学前缘, 2014, 21(3): 265-306.
Cheng H X, Li K, Li M, et al. Geochemical background and baseline value of chemical elements in urban soil in China[J]. Earth Science Frontiers, 2014, 21(3): 265-306.
[28] 马建华, 韩昌序, 姜玉玲. 潜在生态风险指数法应用中的一些问题[J]. 地理研究, 2020, 39(6): 1233-1241.
Ma J H, Han C X, Jiang Y L. Some problems in the application of potential ecological risk index[J]. Geographical Research, 2020, 39(6): 1233-1241.
[29] USEPA. Soil Screening Guidance: Technical Background Document[R]. Research Triangle Park: USEPA, 1996.
[30] 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
[31] 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
[32] 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
[33] 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
[34] 环境保护部. 中国人群暴露参数手册(成人卷)[M]. 北京: 中国环境出版社, 2013.
[35] 黄剑波, 姜登登, 温冰, 等. 基于蒙特卡罗模拟的铅锌冶炼厂周边农田土壤重金属健康风险评估[J]. 环境科学, 2023, 44(4): 2204-2214.
Huang J B, Jiang D D, Wen B, et al. Contamination and probabilistic health risk assessment of heavy metals in agricultural soils around a lead-zinc smelter[J]. Environmental Science, 2023, 44(4): 2204-2214.
[36] Ma J, Chu L J, Sun J, et al. Health risk assessment of trace metal(loid)s in agricultural soils based on Monte Carlo simulation coupled with positive matrix factorization model in Chongqing, southwest China[J]. Journal of Mountain Science, 2024, 21(1): 100-112. DOI:10.1007/s11629-023-8304-2
[37] Yang Z P, Li X Y, Wang Y, et al. Trace element contamination in urban topsoil in China during 2000-2009 and 2010-2019: Pollution assessment and spatiotemporal analysis[J]. Science of the Total Environment, 2021, 758. DOI:10.1016/j.scitotenv.2020.143647
[38] Li F Y, Fan Z P, Xiao P F, et al. Contamination, chemical speciation and vertical distribution of heavy metals in soils of an old and large industrial zone in Northeast China[J]. Environmental Geology, 2009, 57(8): 1815-1823. DOI:10.1007/s00254-008-1469-8
[39] 董鑫, 胡浩然, 张晓晴, 等. 基于Meta分析的矿区周边土壤重金属污染特征及风险评价[J]. 地学前缘, 2024, 31(2): 93-102.
Dong J, Hu H R, Zhang X Q, et al. The distribution characteristics and ecological risk of soil heavy metals in mining areas based on Meta-analysis[J]. Earth Science Frontiers, 2024, 31(2): 93-102.
[40] 王蕊, 陈楠, 张二喜, 等. 龙岩市某铁锰矿区土壤重金属地球化学空间分布特征与来源分析[J]. 环境科学, 2021, 42(3): 1114-1122.
Wang R, Chen N, Zhang E X, et al. Geochemical patterns and source analysis of soil heavy metals in an iron and manganese ore area of Longyan city[J]. Environmental Science, 2021, 42(3): 1114-1122.
[41] Fei X F, Lou Z H, Xiao R, et al. Contamination assessment and source apportionment of heavy metals in agricultural soil through the synthesis of PMF and GeogDetector models[J]. Science of the Total Environment, 2020, 747. DOI:10.1016/j.scitotenv.2020.141293
[42] Liu H, Anwar S, Fang L, et al. Source apportionment of agricultural soil heavy metals based on PMF model and multivariate statistical analysis[J]. Environmental Forensics, 2024, 25(1-2): 40-48.
[43] Sun B, Zhang L X, Yang L Z, et al. Agricultural non-point source pollution in China: Causes and mitigation measures[J]. Ambio, 2012, 41(4): 370-379.
[44] Wu H X, Ge Y. Excessive application of fertilizer, agricultural non-point source pollution, and farmers' policy choice[J]. Sustainability, 2019, 11(4). DOI:10.3390/su11041165
[45] Qin G W, Niu Z D, Yu J D, et al. Soil heavy metal pollution and food safety in China: Effects, sources and removing technology[J]. Chemosphere, 2021, 267. DOI:10.1016/j.chemosphere.2020.129205
[46] Cheraghi M, Lorestani B, Merrikhpour H. Investigation of the effects of phosphate fertilizer application on the heavy metal content in agricultural soils with different cultivation patterns[J]. Biological Trace Element Research, 2012, 145(1): 87-92.
[47] Shi T R, Zhang Y Y, Gong Y W, et al. Status of cadmium accumulation in agricultural soils across China (1975-2016): From temporal and spatial variations to risk assessment[J]. Chemosphere, 2019, 230: 136-143.
[48] Varol M, Gündüz K, Sünbül M R. Pollution status, potential sources and health risk assessment of arsenic and trace metals in agricultural soils: a case study in Malatya Province, Turkey[J]. Environmental Research, 2021, 202. DOI:10.1016/j.envres.2021.111806
[49] Luilo G B, Othman O C, Mrutu A. Arsenic: A toxic trace element of public health concern in urban roadside soils in Dar es Salaam City[J]. Journal of Materials and Environmental Science, 2014, 5(6): 1742-1749.
[50] 喻旗, 罗洁, 涂文忠. 电解金属锰生产的污染及其治理[J]. 中国锰业, 2006, 24(3): 42-45.
Yu Q, Luo J, Tu W Z. Pollution and harness of EMM production[J]. China's Manganese Industry, 2006, 24(3): 42-54.
[51] 樊欣, 罗静, 杨梦, 等. 电解金属锰行业环境污染及合理防治建议解析[J]. 中国锰业, 2016, 34(3): 136-137.
Fan X, Luo J, Yang M, et al. Electrolytic manganese industry parsing the prevention and control of environmental pollution and reasonable advice[J]. China's Manganese Industry, 2016, 34(3): 136-137.
[52] 马杰, 佘泽蕾, 王胜蓝, 等. 基于蒙特卡罗模拟的煤矸山周边农用地土壤重金属健康风险评估[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.
[53] 张素荣, 王昌宇, 刘继红, 等. 雄安新区西南部土壤重金属污染特征及生态风险评价[J]. 地学前缘, 2021, 28(4): 238-249.
Zhang S R, Wang C Y, Liu J H, et al. Assessments of heavy metal pollution in soils of the southwestern Xiong'an District and its ecological risk[J]. Earth Science Frontiers, 2021, 28(4): 238-249.
[54] 雷星宇, 胡瑶, 邵颖, 等. 土壤镉污染对卷丹百合生长和镉分布的影响[J]. 中国土壤与肥料, 2023(9): 216-223.
Lei X Y, Hu Y, Shao Y, et al. Effects of soil cadmium pollution on the growth and cadmium distribution of Lilium lancifolium Thunb[J]. Soil and Fertilizer Sciences in China, 2023(9): 216-223.
[55] 冯汉青, 杜变变, 王庆文, 等. 镉胁迫下活性炭对小麦幼根的保护作用[J]. 生态学报, 2016, 36(10): 2962-2968.
Feng H Q, Du B B, Wang Q W, et al. The role of activated carbon in protecting the roots of wheat seedlings under cadmium stress[J]. Acta Ecologica Sinica, 2016, 36(10): 2962-2968.
[56] 马杰, 佘泽蕾, 王胜蓝, 等. 重庆市煤矸山周边农产品镉健康风险评价及土壤环境基准值推导[J]. 环境科学, 2023, 44(9): 5264-5274.
Ma J, She Z L, Wang S L, et al. Health risk assessment and environmental benchmark of cadmium in farmland soils around the gangue heap of coal mine, Chongqing[J]. Environmental Science, 2023, 44(9): 5264-5274.
[57] 王蕊, 陈楠, 张二喜. 龙岩市不同利用类型土壤及农作物Pb、Cd和As污染风险与贡献分析[J]. 环境科学, 2023, 44(4): 2252-2264.
Wang R, Chen N, Zhang E X. Pollution risk and contribution analysis of Pb, Cd, and As in soils and crops under different land use types in Longyan city[J]. Environmental Science, 2023, 44(4): 2252-2264.
[58] Beesley L, Moreno-Jiménez E, Gomez-Eyles J L. Effects of biochar and greenwaste compost amendments on mobility, bioavailability and toxicity of inorganic and organic contaminants in a multi-element polluted soil[J]. Environmental Pollution, 2010, 158(6): 2282-2287.