黄河下游典型滩区土壤重金属污染特征及来源解析.pdf

黄河下游典型滩区土壤重金属污染特征及来源解析徐祖奔1,2，伍艳1,2，赵越1,2，乔奥克3，冯帅滔3，刘振延3（1.黄河水利委员会黄河水利科学研究院，郑州450003；2.水利部堤防安全与病害防治工程技术研究中心，郑州450003；3.河南黄科工程技术检测有限公司，郑州450003）摘要：针对黄河滩区土壤重金属污染、来源复杂的问题，该研究以原阳滩区复合污染土壤为例，结合土壤重金属含量空间分布和正定矩阵因子分解（positivematrixfactorization,PMF）模型，探讨滩区土壤重金属富集特征及来源。源解析结果表明，原阳滩区土壤重金属污染受工业源、交通源、自然来源、燃煤污染源和农业源影响，农业源占主导，相对贡献率达 23.5%，其次为工业源、自然来源、交通源、和燃煤污染源。该研究能够准确解析多金属复合土壤复杂成因，可为黄河滩区多金属复合污染土壤的污染源识别提供参考。关键词：重金属；风险评价；正定矩阵因子分解法；源解析；黄河滩区doi：10.11975/j.issn.1002-6819.202302055中图分类号：X53；S159文献标志码：A文章编号：1002-6819(2023)-15-0200-08徐祖奔，伍艳，赵越，等.黄河下游典型滩区土壤重金属污染特征及来源解析J.农业工程学报，2023，39(15)：200-207.doi:10.11975/j.issn.1002-6819.202302055http:/www.tcsae.orgXUZuben,WUYan,ZHAOYue,etal.SourceapportionmentandpollutionofsoilheavymetalsintypicalfloodplainintheLowerReachesoftheYellowRiverJ.TransactionsoftheChineseSocietyofAgriculturalEngineering(TransactionsoftheCSAE),2023,39(15):200-207.(inChinesewithEnglishabstract)doi:10.11975/j.issn.1002-6819.202302055http:/www.tcsae.org0引言土壤污染已经成为全球性的环境问题，其中重金属污染尤其突出。土壤中的重金属污染分为自然源和人为源，自然源是土壤形成过程的内部来源，包括矿物风化和背景源的元素淋滤。而人为源是土壤重金属的外部来源，人类活动，如工业燃煤、采矿、冶炼，农业化肥、农药使用，生活中汽车尾气，电子产品丢弃等，都会导致不同程度的土壤重金属污染。然而，重金属污染具有隐秘性强、持久性和生物蓄积性等特点1，很难从土壤中去除。过量的重金属不仅导致土壤中植物发育受阻，造成作物减产和农产品品质下降，甚至直接或者间接通过多种途径进入人体，导致多种疾病2-3。例如，过量Pb 摄入不仅会导致高血压，损害骨骼、免疫系统、内分泌，降低儿童和成年人的智力潜能，还会影响肾脏和心脏的功能4。过量摄入 Cu、Hg、Cr 也会对健康产生有害影响5。此外，长期接触 Cd 的人更容易患肺癌和骨折。因此，土壤重金属污染与人类生产生活息息相关，防控重金属污染已迫在眉睫6。近年来，人们对土壤重金属污染的关注与研究日益增多，从土壤污染特征、重金属风险评价和污染源解析方法等方面开展了广泛的研究7-8。土壤污染特征主要包括空间分布特征和分层异质性等，常用的方法有地理信息系统9、地统计学方法等10；土壤重金属污染评价用于评估多种重金属含量和形态对动物、植物和人类可能存在的风险，常用的重金属污染评价方法有单因子指数法11、内梅罗综合污染法12、富集因子法13、潜在生态风险指数14和地累积指数15等。污染物源解析主要用于定性确定污染物的来源，并且可以量化污染源的贡献。目前，最常用的重金属源解析方法有化学物质平衡（chemicalmassbalance，CMB）模型16、正定矩阵因子分 解（positive matrix factorization，PMF）模 型17、UNMIX 模型18、多元统计分析模型19等。在这些方法和模型中，PMF 模型不仅可以定量分析重金属，确保非负源贡献值和因子分布，而且具有不需要测量源成分谱、分解矩阵中非负元素、利用数据的标准差进行优化的优点20，因此 PMF 模型已成为美国环保署推荐的源解析方法，并成功地应用于识别和量化土壤中的重金属污染。黄河滩区是黄河两侧重要的水陆交替区域，不仅发挥着泄洪、沉积黄河泥沙的重要功能，还是百万群众赖以生存的场地。随着人类活动的加剧，滩区土壤重金属污染问题突出。因此，了解滩区重金属污染状况，解析滩区重金属污染来源，对于解决滩区重金属污染问题、保护滩区群众健康安全具有重要意义。当前，针对土壤重金属污染来源的研究主要是利用源解析模型进行解析，但是黄河滩区经济活动剧烈，社会活动多元化，导致滩区土壤多重金属污染复合，单一的源解析模型面临巨大的挑战。而基于地理空间插值分析和 PMF 模型结合的方式能够在空间尺度对源解析结果进行优化，进一步提高收稿日期：2023-02-11修订日期：2023-07-11基金项目：中央级公益性科研院所基本科研业务费专项项目（HKY-JBYW-2021-03）作者简介：徐祖奔，研究方向为土壤重金属风险评价及源解析。Email：通信作者：伍艳，博士，正高级工程师，研究方向为土壤环境质量评价及修复。Email：wuyan3000_第39卷第15期农 业 工 程 学 报Vol.39No.152002023年8月TransactionsoftheChineseSocietyofAgriculturalEngineeringAug.2023源解析结果的准确性。同时，现有的研究多聚焦于黄河流域水体和沉积物，对于滩区区域性，较大样本溯源评价研究较少。因此，对黄河滩区土壤重金属的污染特征、风险评估和来源解析开展深入研究，更能有效地保障黄河流域生态保护和高质量发展。原阳滩区是黄河东西走向流经原阳南部形成的宽广滩地，位于新乡市西南隅，南滨黄河，东接封丘，隔河与开封相望，土质肥沃，是良好的作物生产基地，具有黄河流域、农业生产和人类生活三者交汇的典型特征。因此，本文以原阳滩区滩地土壤为研究对象，测定样点土壤中 8 种重金属元素的含量，利用反距离加权（inversedistanceweighted，IDW）空间插值分析土壤重金属含量的空间分布特征，采用地累积指数法和潜在生态风险指数对原阳滩区土壤重金属进行综合评价，并结合 PMF 模型确定研究区污染来源，分析污染主要来源贡献。以期通过对原阳滩区土壤环境质量综合分析，为黄河流域滩区生态保护和土壤污染修复提供科学依据。1材料与方法1.1研究区概况研究区位于河南省原阳县境西南部黄河滩地，地理位置为 1133711415E，34523501N，南濒黄河，北依黄河大堤。原阳滩区横跨辖内桥北乡、韩董庄乡、蒋庄乡、官厂乡、陡门乡 5 个滩区乡，总面积约 200km2。土地利用类型以耕地为主，林地和住宅用地次之，工业用地零星分布。原阳为河南省主要的粮食生产基地，境内农作物有小麦、玉米、花生、大豆、水稻、棉花等。原阳滩区产业结构主要由工业、轻工业、农业结构组成，其境内分布有机械修理厂、塑料厂、装饰厂、玻璃加工厂、废品收购厂、建材厂等。该区域工业废弃物和电子废弃物造成的污染并未得到合理的处理和处置，长期的农业灌溉及农业生产也造成大量的土地污染。而且，被污染的土壤通过降雨径流等方式，成为黄河生态环境的污染来源。1.2土壤样品的采集采样单元以 1km 1km 为网格布点，GPS 记录采样点实际坐标，采集 020cm 表层土壤，样品采集参照GB/T36197-2018土壤质量土壤采样技术指南，采样布点图见图 1，共 81 个样点。采样时，去除根茬、落叶、砾石等杂物，用“对角线四分法”取上下两份土样，每份不少于 1kg 并充分混合均匀。采集后土样装入样品袋中，贴好标签，记录采样点 GPS 坐标，在常温保存备用21。1.3重金属含量的测定样品风干后，过 0.149mm 筛，称取土样 0.10.5g（精度为 0.0001g）置于已经洗净并干燥的 50mL 聚四氟乙烯坩埚之中，采用“硝酸-氢氟酸-高氯酸”体系（GB/T17141-1997）在电加热板上消解，其含量用电感耦 合 等 离 子 质 谱 仪（XSeries-2 型 ICP-MS，美 国ThermoFisher）测定22。在重金属测定过程中，重复测量 3 次取平均值，设定标准空白并其中 20%样品进行了平行试验，同时每隔10 次测量进行一次标样质控，各标准物质回收率均在80%120%之间，标准偏差在 5%以内。113500E113500E113400E113400E3550N3550N3500N3500N34550N34550N0 2 4km原阳县河南省采样点Sampling sites黄河The Yellow River研究区Study areaN图 1研究区概况及采样点位图Fig.1Overviewofstudyareaandlocationofsamplingsites1.4重金属污染评价方法本研究采用德国学者 MULLER 提出23的地累积指数法（geoaccumulationindex，Igeo）对研究区表层土壤重金属污染程度进行评价；采用 HAKANSON24建立的潜在生态风险指数法评价研究区土壤生态风险水平，毒性响应系数Hg（40）Cd（30）Pb=Cu（5）Cr（2）Zn（1），分级标准见表 1。表1地累积指数 Igeo与潜在生态风险指数分级Table1Classificationofgeoaccumulationindex(Igeo)andpotentialecologicalriskindexIgeo污染程度PollutiondegreeEiRI生态危害程度Ecologicalharm0无(0,40)(0,150)轻微01轻度-中等40,80)150,300)中等12中等80,160)300,600)强23中等-强160,320)600,1200)很强34强320,+)1200,+)极强危害45强-极强5极强注：Ei和 RI 分别为单项和综合潜在生态风险指数。Note:EiandRIaresingle,comprehensivepotentialecologicalriskindex,respectively.1.5PMF 源解析QQPMF模型是根据源的组成或指纹量化源样本的贡献，其每个受体原始数据被分解为因子得分和因子载荷的乘积加上残差，其方程见式（1），通过加权最小二乘法进行限定和迭代计算，使目标函数最小化，目标函数的计算方程见式（2）8。xij=pk=1gikfkj+eij,(i=1,2,n;j=1,2,m)（1）Q=ni=1mj=1xijpk=1gikfkjuij2=ni=1mj=1(eijuij)2（2）xijijgikkifkjkjeijeiji式中为第 个样品中第 个元素的测量浓度（mg/kg）；为源 对第 个样品的相对贡献：为源 中第 个元素的浓度（mg/kg）；为残差。为第 个样品中第15期徐祖奔等：黄河下游典型滩区土壤重金属污染特征及来源解析201juijij第 个元素浓度的不确定度。为第 个样品中第 个元素含量的不确定性大小。1.6数据处理与分析土壤重金属含量、重金属风险评价采用 Excel 进行统计和计算；土壤重金属相关性和主成分通过 SPSS 进行分析；采用 ArcGIS10.5 制作土壤采样点概况图和反距离加权空间插值图。2结果与分析2.1土壤重金属含量特征研究区土壤重金属含量统计分析结果如表 2 所示。重金属平均浓度由大到小依次为 Zn（88.25mg/kg）、Cr（73.55mg/kg）、As（29.00mg/kg）、Pb（27.77mg/kg）、Cu（27.05mg/kg）、Ni（25.38mg/kg）、Cd（0.31mg/kg）和 Hg（0.24mg/kg）。与河南省土壤背景值相比（见表 2），除 Ni 外，原阳滩区土壤重金属的平均浓度均高于背景限值，说明表层土壤重金属含量受到外部和人为因素的影响。特别注意的是，Hg 的平均含量约为土壤背景值的 7 倍，说明 Hg 在表层土壤中富集。对比中国土壤环境质量农用地土壤污染风险管控标准（GB15618-2018）（pH 值7.5），只有 As 元素含量均值高于土壤风险筛选值，表明原阳滩区 As 污染严重，存在重金属污染生态风险。由表 2 可知，Cr、Zn、As、Hg 的变异系数（variationcoefficient，CV）较大，说明这些金属的分布不均匀19。表2原阳滩区土壤重金属含量Table2HeavymetalconcentrationsinYuanyangfloodplainsoil指标Index含量 Concentration/(mgkg1)变异系数Variationcoefficient/%BV/(mgkg1)RS/(mgkg1)平均值Mean最小值Min最大值Max标准差SDCr73.5523.26307.8246.9863.963.8250Ni25.387.5145.839.7738.526.7190Cu27.055.0954.8210.1837.619.7100Zn88.2521.95390.6857.6565.360.1300As29.002.0470.1021.9475.711.425Cd0.310.140.570.1033.40.0740.6Pb27.7711.3756.348.6531.119.6170Hg0.2402.470.50206.60.0343.4注：BV 为河南省土壤背景值；RS 为土壤风险筛选值（pH 值7.5）。Note:BVrepresentssoilbackgroundvalueinHenanprovince;RSrepresentsriskscreeningvaluesforsoilcontamination（pHvalue7.5）.2.2土壤重金属含量空间分布特征为了进一步表征土壤重金属的人为来源，采用反距离权重插值（IDW）分析了土壤中重金属的空间分布。8 种重金属在土壤中的浓度分布如图 2 所示，研究区土壤中重金属的空间分布具有明显的特征。整体上看，原阳滩区重金属呈岛状分布，Cu、Ni、Pb、Cd 污染较为广泛，且高值集中在研究区的中部蒋庄乡。Zn、Cr、As和 Hg 高值分布存在空间重叠，说明该区域表层土壤受人类活动的影响显著且很可能存在点源污染。通过资料收集和现场调查，原阳县城附近乡镇拥有电子厂、化肥厂、木材加工厂、机械修理厂、蓄电池修理厂、塑料厂等，在研究区中部为蒋庄乡西北方向，公路旁有一处较大机械修理厂，推测该区域土壤污染原因主要为工业活动产生的重金属经淋溶和大气沉降25。As/(mgkg1)2.0611.7811.7921.5021.5131.2131.2240.9340.9450.6530.3770.0950.6660.36Pb/(mgkg1)11.3817.8017.8124.2224.2330.6430.6537.0637.0743.4949.9256.3443.5049.91Cd/(mgkg1)0.130.190.200.250.260.320.330.380.390.440.510.570.450.50Hg/(mgkg1)00.350.360.700.711.051.061.411.421.762.122.481.772.11Ni/(mgkg1)7.5413.0013.0118.4718.4823.9323.9429.3929.4034.8544.3345.7834.8640.32Cr/(mgkg1)23.2863.9163.92104.54104.55145.17145.18185.79185.80226.42267.06307.68226.43267.05Zn/(mgkg1)22.0374.6874.69127.34127.35180.00180.01232.65232.66285.31337.97390.62285.32337.96Cu/(mgkg1)47.7154.8040.6147.7033.5240.6026.4233.5119.3226.4112.2319.325.1212.22N01234km图 2原阳滩区土壤重金属含量空间分布Fig.2SpatialdistributionofheavymetalconcentrationsinYuanyangfloodplain2.3生态风险评价2.3.1地累积指数法重金属的地累计指数结果如图 3 所示。Cr Ni Cu Zn As Cd Pb Hg20246第1090百分位数范围10th90th percentile1.5倍四分位距内的范围1.5 times Interquartile Range中位线Median line均值Mean value地累计指数Geoaccumulation index 重金属Heavy metals图 3土壤重金属地累计指数统计结果箱式图Fig.3Boxplotsofgeoaccumulationindexforheavymetals8 种重金属的地累积指数均值为 Pb（0.36）Cd（0.40）Cr（0.45）Cu（0.48）Ni（0.49）Zn（0.61）As（1.27）Cd(127.26)As(25.44)Pb(7.08)Cu(6.86)Ni(5.70)Cr(2.31)Zn(1.47)，其中重金属 Hg 的平均潜在生态风险最高，属于很强生态危害；Cd 其次，为强生态危害，其他6种重金属为轻微生态危害。以上结果表明，Hg 和 Cd 是导致该区域生态危害的主要元素，可能会危及当地生态环境和人体健康。研究区各样点重金属的综合潜在生态风险指数（RI）为 90.23057.11，其中轻微、中等、强、很强和极强生态危害所占的比例分别为 29.63%、37.04%、8.64%、16.05%和 8.64%。研究区样点以轻微和中等强度危害为主，其综合潜在生态风险指数平均值为 459.31，属于强生态危害水平。Hg 和 Cd 对综合潜在生态风险指数的贡献率分别为 61.66%和 27.71%，合计 89.37%，表明研究区表层土壤重金属的潜在生态风险主要是由 Hg 和 Cd 引起。对重金属综合潜在生态风险指数进行空间插值分析，见图 4，其中西边桥北乡和中段蒋庄乡重金属污染最为严重，属于很强生态风险水平，说明该区域土壤 Hg 和Cd 污染较重，对环境生态危害最大。轻微危害Slight中等危害Intermediate强危害Strong很强危害Very strong极强危害FortissimoN0 1 2 3 4km图 4土壤重金属潜在生态风险空间分布Fig.4Spatialdistributionofpotentialecologicalriskofsoilheavymetals2.4土壤重金属源解析2.4.1聚类分析对研究区 8 种重金属进行聚类，聚类分析结果如图 5所示，8 种元素可以被划分为四大类，类为 Cd 和 Hg，类为 Ni、Pb、Cu 和 As，类为 Cr，类为 Zn。各金属间 Person相关系数见表 3，结果显示 Ni、Cu、Pb两两之间存在极显著正相关性（P0.98），说明 PMF 模型总体拟合效果较好，因子能够准确反映原始数据信息。因子 1 的主要载荷是 Pb，贡献率为 47.1%。研究表明，铅是交通排放的主要指标，可以通过车辆燃油燃烧、汽车发动机和轮胎摩擦来释放17,27-28。在本研究中，农田道路分布机动车、农业机械设备排放的废气中含铅，通过大气沉降对土壤造成污染。此外，Cd、Zn 和 Ni 浓第15期徐祖奔等：黄河下游典型滩区土壤重金属污染特征及来源解析203度的贡献率分别为 40.8%、36.1%和 33.2%，Cd、Zn 和Ni 在汽车尾气中含量丰富，这些污染物通过大气沉降和空气尘埃吸附在耕作土壤中积累29。此前，GUAN 等30均认为 Pb 是由于燃料燃烧、发动机和催化剂的使用而产生的交通排放的主要标志。因此，因子 1 可以被认为是交通运输的来源。表4重金属的源成分谱及贡献率Table4Sourcesspectraandsourcecontributionofheavymetals元素Elements源成分谱 Sourcesspectra/(mgkg1)源贡献率 Sourcecontribution/%因子 1Factor1因子 2Factor2因子 3Factor3因子 4Factor4因子 5Factor5因子 1Factor1因子 2Factor2因子 3Factor3因子 4Factor4因子 5Factor5Cr021.0931.3918.192.74028.742.824.83.7Ni8.3909.815.981.1233.2038.823.74.4Cu3.982.572.1517.191.1414.79.57.963.64.2Zn26.293.6121.8621.08036.15.030.028.90As0.5302.274.3421.771.807.915.075.3Cd0.1220.046000640.06740.815.4021.422.3Pd13.001.484.922.895.3347.15.417.810.519.3Hg0.0140.4150.046003.087.39.700因子 2 主要载荷元素为 Hg，贡献率高达 87.3%，同时 Cr 的贡献也有 28.7%。研究表明，Hg 的来源复杂，受土壤母质、工业排放（如燃煤）、交通、农业（如高Hg 农药）等多因素影响31。研究区 Hg 的变异系数高达206.6%，为强变异，说明 Hg 污染与人类活动密切相关。研究区的 Hg 平均含量约为土壤背景的 7 倍，说明 Hg 污染为高污染来源，同时由图 2 土壤中 Hg 的空间分布可知，高 Hg 污染呈斑状分布，说明 Hg 的来源大概率来至工业排放。同时 CHEN 等32-33研究表明，土壤表层 Hg的累积受工业排放和燃煤影响较大，结合实地调查，该区域之前分布有玻璃加工厂、钢材厂和烧砖厂等高燃煤工厂，说明 Hg 污染很可能与燃煤或钢材加工引起的大气沉降相关34。因此，因子 2 可以解释为燃煤污染源。因子 3 的主要载荷元素为 Cr 和 Ni，贡献率分别为46.8%和 42.0%，明显高于其他元素。土壤中 Cr 和 Ni 的含量通常与其成土母质中的含量有关，这些元素很大可能参与成岩8,35。由表 2 可知，Ni 的平均含量低于河南省背景值，且变异系数（38.5%）较小，说明受人为活动干扰较小，说明其为自然母质的背景含量。这与张旺等19得出 Cr 和 Ni 主要受成土母质影响的结论一致。因此，因子 3 可以解释为地质背景引起的自然来源。因子 4 对 Cu 和 Zn 的贡献率较高，分别为 63.6%和28.9%，其次是 Ni、Cr 和 Cd。研究表明，Cu、Zn 和 Cd与农业活动有关，肥料、农药和灌溉水对这几种元素都有一定的贡献17,36。LUO 等37-38也证实铜和锌是肥料和农用化学品进入农业土壤的一个重要来源。因此，因子4 可解释为农业源。因子 5 中 As 的贡献率最高为 75.3%，远高于其他因子。研究表明，土壤中 As 的来源主要受土壤母质、大气沉降、农业活动及交通等影响19，研究区 As 的污染浓度平均值为 29.00mg/kg，是土壤背景值的两倍多，且变异系数高达 75.7%，说明 As 污染主要受人类活动影响。结合研究区 As 污染空间分布（图 2），As 污染存在明显的高值分布，说明污染源分布不均呈零散分布，同时因子 1 中 Cd 和 Pb 也分别有 22.3%和 19.3%的贡献，说明 As 与 Cd、Pb 具有同源性。而 As、Cd 和 Pb 的来源通常与工业活动存在很大关系，化石燃料燃烧，金属冶炼等都会向大气中排放 As、Cd 和 Pb，最终以大气沉降的方式进入土壤39-40。实地勘察发现，研究区周边之前存在大量违规乡镇企业（机械厂、建材厂、蓄电池、电子厂、修理站等），虽然现已拆迁，但是由于工业活动长期积累，土壤中 As、Cd 和 Pb 的含量有所增加。因此，因子 5 可以被认为是工业活动的一个来源。综上结果可知，原阳滩区重金属污染来源受多种因素影响，包括工业源、交通源、自然来源、燃煤污染源和农业源，其相对贡献率分别为 16.2%、22.1%、19.4%、18.9%和 23.5%。3结论1）原阳滩区土壤重金属 Hg、Cd、As、Pb、Cu、Zn、Cr 平均含量均超过了河南省背景值，除 Cd 外，7 种元素的重金属含量都低于 GB15618-2018（pH 值为 6.57.5）中的筛选值。Hg 的平均含量最高，约为背景值的7 倍，且表现为高变异性，人类活动干扰显著。2）地累计指数评价结果表明，Hg 污染最为严重，其次为 As 和 Cd。潜在生态风险评价由强到弱依次为：HgCdAsPbCuNiCrZn，重金属Hg具有很强生态危害，Cd具有强生态危害，其他重金属属于轻微生态危害。从空间分布上看，位于原阳滩区西边桥北乡和中段蒋庄乡污染最为严重，且以 Hg 和 Cd 为主要污染。3）聚类分析和 PMF 分析结果表明，原阳滩区土壤重金属污染来源为工业源、交通源、自然来源、燃煤污染源和农业源，其相对贡献率分别为 16.2%、22.1%、19.4%、18.9%和 23.5%。参考文献李鸿博，钟怡，张昊楠，等.生物炭修复重金属污染农田土壤的机制及应用研究进展J.农业工程学报，2020，36(13)：173-185.LIHongbo,ZHONGYi,ZHANGHaonan,etal.MechanismfortheapplicationofbiocharinremediationofheavymetalcontaminatedfarmlandanditsresearchadvancesJ.TransactionsoftheChineseSocietyofAgricultural1204农业工程学报（http:/www.tcsae.org）2023年Engineering(TransactionsoftheCSAE),2020,36(13):173-185.(inChinesewithEnglishabstract)ZENGJ,LUOX,CHENGY,etal.Spatialdistributionoftoxicmetal(loid)s at an abandoned zinc smelting site,SouthernChinaJ.JournalofHazardousMaterials,2022,425:1-9.2QING,NIUZ,YUJ,etal.Soilheavymetalpollutionandfood safety in China:Effects,sources and removingtechnologyJ.Chemosphere,2021,267:1-12.3OBENG-GYASI E,ROOSTAEI J,GIBSON J M.LeaddistributioninurbansoilinaMedium-SizedCity:Household-scaleanalysisJ.EnvironmentalScience&Technology,2021,55(6):3696-3705.4AHMADW,ALHARTHYRD,ZUBAIRM,etal.Toxicandheavy metals contamination assessment in soil and water toevaluatehumanhealthriskJ.ScientificReports,2021,11(1):2045-2322.5YANGQ,LIZ,LUX,etal.Areviewofsoilheavymetalpollution from industrial and agricultural regions in China:Pollution and risk assessmentJ.Science of the TotalEnvironment,2018,642:690-700.6ZHANG Y,LI T,GUO Z,et al.Spatial heterogeneity andsource apportionment of soil metal(loid)s in an abandonedlead/zincsmelterJ.JournalofEnvironmentalSciences,2023,127:519-529.7FEIX,LOUZ,XIAOR,etal.Contaminationassessmentandsource apportionment of heavy metals in agricultural soilthrough the synthesis of PMF and GeogDetector modelsJ.ScienceoftheTotalEnvironment,2020,747:141293.8JINH,ZHIHONGP,JIAQINGZ,etal.Sourceapportionmentand quantitative risk assessment of heavy metals at anabandoned zinc smelting site based on GIS and PMFmodelsJ.JournalofEnvironmentalManagement,2023,336:1-11.9QUM,CHENJ,HUANGB,etal.Enhancingapportionmentof the point and diffuse sources of soil heavy metals usingrobust geostatistics and robust spatial receptor model withcategorical soil-type dataJ.Environmental Pollution,2020,265:114964.10WANGYB.Ecologicalriskidentificationandassessmentofland remediation project based on GIS technologyJ.Environmental Science and Pollution Research,2023,30:70493-70505.11WEIJ,ZHENGXG,LIUJT.ModelinganalysisofheavymetalevaluationincomplexgeologicalsoilbasedonNemerowindexmethodJ.Metals,2023,13(2):439.12LI Y Y,ZHOU H D,GAO B,et al.Improved enrichmentfactor model for correcting and predicting the evaluation ofheavymetalsinsedimentsJ.ScienceoftheTotalEnvironment,132021,755:142437.LISC,XIAOW,ZHAOYL,etal.Incorporatingecologicalriskindexinthemulti-processMCREmodeltooptimizetheecologicalsecuritypatterninasemi-aridareawithintensivecoal mining:A case study in northern ChinaJ.Journal ofCleanerProduction,2020,247:119143.14FEI X,XIAO R,CHRISTAKOS G,et al.Comprehensiveassessment and source apportionment of heavy metals inShanghai agricultural soils with different fertility levelsJ.EcologicalIndicators,2019,106:1-9.15ZHANGQF,SHIB,SUGJ,etal.ApplicationofahybridGEM-CMB model for source apportionment of PAHs in soil ofcomplex industrial zoneJ.Journal of Hazardous Materials,2023,445:130565.16ANAMANR,PENGC,JIANGZ,etal.IdentifyingsourcesandtransportroutesofheavymetalsinsoilwithdifferentlandusesaroundasmeltingsitebyGISbasedPCAandPMFJ.ScienceoftheTotalEnvironment,2022,823:153759.17周蓓蓓，郭江，陈晓鹏，等.基于 unmix 模型的安徽大矾山废弃矿区土壤重金属源解析J.农业工程学报，2021，37(24)：240-248.ZHOU Beibei,GUO Jiang,CHEN Xiaopeng,et al.SourceapportionmentofsoilheavymetalsinabandonedminingareasinDafanMountainofAnhuiProvincebasedontheUNMIXmodel.J.TransactionsoftheChineseSocietyofAgriculturalEngineering(TransactionsoftheCSAE),2021,37(24):240-248.(inChinesewithEnglishabstract)18张旺，高珍冉，邰粤鹰，等.基于 apcs-Mlr 受体模型的贵州喀斯特矿区水田土壤重金属源解析J.农业工程学报，2022，38(3)：212-219.ZHANG Wang,GAO Zhenran,TAI Yueying,et al.Sourceanalysis of the heavy metals in paddy field soils in KarstminingareasofGuizhouusingAPCS-MLRreceptormodelJ.TransactionsoftheChineseSocietyofAgriculturalEngineering(TransactionsoftheCSAE),2022,38(3):212-219.(inChinesewithEnglishabstract)19ZHOU L,ZHAO X,MENG Y,et al.Identification prioritysourceofsoilheavymetalspollutionbasedonsource-specificecologicalandhumanhealthriskanalysisinatypicalsmeltingand mining region of South ChinaJ.Ecotoxicology andEnvironmentalSafety,2022,242:113864-113864.20王海洋，韩玲，谢丹妮，等.矿区周边农田土壤重金属分布特征及污染评价J.环境科学，2022，43(4)：2104-2114.WANGHaiyang,HANLing,XIEDanni,etal.Distributioncharacteristicso