环境科学  2018, Vol. 39 Issue (3): 1421-1429   PDF    
基于UNMIX模型的矿区周边农田土壤重金属源解析
卢鑫1,2,3, 胡文友1, 黄标1, 李元2, 祖艳群2, 湛方栋2, 邝荣禧4     
1. 中国科学院土壤环境与污染修复重点实验室(南京土壤研究所), 南京 210008;
2. 云南农业大学资源与环境学院, 昆明 650201;
3. 上海市岩土地质研究院有限公司, 上海 200072;
4. 江苏省地质调查研究院, 南京 210018
摘要: 采集了云南省会泽县铅锌矿区周边42个农田土壤样品,测定了14种元素的含量,应用UNMIX模型进行了土壤重金属源解析的研究,并利用ArcGIS地统计分析模块中的反距离加权插值法分析了污染严重的6种重金属的空间分布,进一步验证源解析结果.结果表明:①研究区土壤中Pb、Zn、Cd污染较为严重,均超过当地土壤背景值的数十倍.②UNMIX模型解析出的3个土壤重金属污染来源分别为工业活动造成的人为污染源(源1),源贡献率为16.32%;燃煤和施肥导致的污染源(源2),源贡献率为68.26%;矿山开采导致的人为污染源及土壤母质造成的自然污染源的综合污染源(源3),源贡献率为15.42%.③研究区的农田土壤重金属污染的空间分布格局与当地的土地利用类型和UNMIX模型解析的结果基本吻合,表明UNMIX模型可以很好地应用于矿区周边农田土壤重金属源解析研究.
关键词: 矿山开采      农田土壤      重金属      UNMIX模型      源解析     
Source Apportionment of Heavy Metals in Farmland Soils Around Mining Area Based on UNMIX Model
LU Xin1,2,3 , HU Wen-you1 , HUANG Biao1 , LI Yuan2 , ZU Yan-qun2 , ZHAN Fang-dong2 , KUANG Rong-xi4     
1. Key Laboratory of Soil Environment and Pollution Remediation, Institute of Soil Science, Chinese Academy of Sciences, Nanjing 210008, China;
2. College of Resources and Environment, Yunnan Agricultural University, Kunming 650201, China;
3. Shanghai Geotechnical Engineering & Geology Institute Co., Ltd., Shanghai 200072, China;
4. Geological Survey of Jiangsu Province, Nanjing 210018, China
Abstract: Forty-two farmland soil samples were collected from the mining area of Huize City, Yunnan Province, to identify the sources of metals in the farmland soils, and 14 selected elements in farmland soils were detected. The UNMIX model was applied to identify the sources of the metals. The inverse distance weighted interpolation method of statistical analysis in ArcGIS was used to develop a spatial distribution map of the six severely polluted metals to validate the source apportionment results of the UNMIX model. The results indicated that:① the metals, Pb, Zn, and Cd, in the studied farmland soils were seriously polluted according to the soil background values of Yunnan Province and the Chinese environmental quality standard; ② there were three potential sources of metals according to UNMIX model:source 1 was the anthropogenic sources caused by industrial activities with a source contribution rate of 16.32%; source 2 was the anthropogenic sources caused by coal combustion and fertilization with a source contribution rate of 68.26%; and source 3 was the combined sources caused by mining related activities and soil parent materials with a source contribution rate of 15.42%; and ③ the spatial distribution of selected metals was in accordance with the local land-use and the results of the UNMIX model. Therefore, the UNMIX model can be effectively applied to the source apportionment of heavy metals in farmland soils around the mining area.
Key words: mining      farmland soil      heavy metal      UNMIX model      source apportionment     

土壤是环境的重要组成部分, 是人类赖以生存的自然环境和农业生产的重要资源.矿山资源的开发, 在对国家和地方经济产生巨大推动作用的同时, 也破坏了矿区周边的生态环境, 造成当地的大气、水体、土壤等污染问题[1, 2].矿山开采活动是造成矿区周边农田土壤中重金属污染的主要原因之一, 已引起了广泛的关注.目前大多数学者对采矿引起的土壤重金属污染问题主要集中在矿区周边土壤重金属含量[3~5]和赋存形态及其生物有效性[4, 6, 7]、重金属空间分布[5, 8~10]、污染评价[11~13].而分析矿区周边土壤的重金属污染来源及其相对贡献, 可以更加科学地对污染状况进行环境风险评价, 对制定更加有效的环境管理措施具有重要意义[14, 15].

污染物源解析的研究最早起源于20世纪60年代的美国, 主要采用以污染源为对象的扩散模型和以污染区域为对象的受体模型进行定性或者定量分析[16].扩散模型是根据已知各个污染源排放量、排放源与研究区域距离、污染物的理化性质以及研究区域中风向、风速和湍流等自然环境要素来研究污染源对该区域的影响程度的一种源解析方式[17]; 而受体模型则是针对污染区域作为研究对象, 直接通过分析排放源和受体样品本身的理化性质来分析污染源对受体的影响, 可以避免由于地形或气象数据的复杂性和不确定造成的源解析结果的偏差, 越来越受到广泛的应用[18].当下的受体模型主要有主成分分析法(principal component analysis, PCA)[19, 20]、正定矩阵因子分析法(positive matrix factorization, PMF)[21, 22]、化学质量平衡法(chemical mass balance, CMB)、因子分析法(factor analysis, FA)和UNMIX[23].目前, 大多数受体模型都是应用于大气污染物的源解析研究中, 如PCA、PMF和UNMIX是对空气环境中悬浮颗粒物、大气或水体环境中多环芳烃(PAHs)和多氯联苯(PCBs)等有机污染物进行源解析研究最常用的3种受体模型[24, 25].对于利用受体模型进行土壤中污染物的源解析研究方面, 我国采用更多的是因子分析法(FA)[26~28].

UNMIX模型由美国环境保护署(U. S. environmental protection agency, USEPA)开发, 主要用于对某一区域内对空气质量产生影响的各种空气污染源的贡献情况进行识别和量化, 从而为制定和实施相应的空气质量标准提供科学依据的几种源解析模型之一, 在国内外得到了广泛的应用. UNMIX方法相比其他受体模型, 具有如下特点:它直接以受体含量计算, 不需要经过数据转换, 模型自带的数据系统可以自动去除不合理的数据; 它事先不需要知道污染源的个数、源成分谱等信息, 直接从受体含量出发通过污染物标识物及指纹元素来识别、推断污染源; 它克服了其他受体模型在解析结果给出的负值源贡献的缺点, 使解析结果更加准确、合乎实际. Pancras等[29]通过对在一定时期内密集采集自佛罗里达州坦帕市的可吸入颗粒物PM2.5中的重金属UNMIX源解析, 得出其可能来源为燃油源、燃煤源、生物质燃烧源、海洋气溶胶源、工业排放源、降尘源等8个来源的结论. Callén等[30]利用UNMIX模型对西班牙Zaragoza地区的PM10进行源解析, 发现对该地区空气中PM10的含量影响较大的有工业或交通源、海洋气溶胶源、生物质燃烧源、地壳源和重型车尾气排放源这5种来源. Lang等[31]应用UNMIX模型对胶州湾湿地土壤中的多环芳烃(PAHs)进行源解析得到其可能的来源有石油源、炼焦炉源和燃煤源, 且三者的源贡献分别为43.2%、39.2%和12.6%.冯小琼等[32]综合应用CMB模型、UNMIX模型对香港地区2005~2010年VOCs污染来源的解析后认为溶剂使用、机动车尾气排放和液化石油气(LPG)的大量使用是造成该地区的VOCs污染的主要原因.

作为广泛应用于空气污染物源解析研究的受体模型, UNMIX应用到土壤中污染物的解析研究鲜见报道.本研究应用UNMIX模型, 参考对大气污染物的源解析的研究方法, 解析云南会泽铅锌矿区周边农田土壤金属的来源, 以期为矿区周边农田土壤的重金属污染防控与治理提供科学依据.

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

会泽县位于云南省东北部, 介于东经103°03′~103°55′, 北纬25°48′~27°04′之间, 平均海拔2 120 m, 气候温和, 雨量充沛, 年均气温12.6℃, 年均降雨量858.4 mm, 属温带山地半潮湿气候, 其铅锌矿是川滇黔铅锌成矿区的大型富铅锌矿床的一个典型代表, 作为我国有名的土法炼锌集散地, 会泽地区的采矿历史最早可以追溯到西汉时期[33].历史上长期的没有环保措施的采矿和冶炼行为在当地遗留下了大量的开矿、冶矿废弃地, 在风吹、雨水淋溶的作用下, 废矿废渣中的重金属元素在周围的环境中扩散, 导致矿区周边土壤中Pb、Zn、Cd和As等重金属高度累积[34].

1.2 样品采集与分析测定

为研究矿区周边农田土壤重金属污染情况, 供试土壤取自云南省曲靖市会泽县矿山镇到者海镇之间(103°36′~103°43′E, 26°33′~26°38′N), 从山谷上游的矿山镇农田到逐步接近下游的者海镇矿区周边农田为采样区域, 依次随机布点采样, 共计42个土壤采样点(图 1).在每个采样点周围随机采集5个深度为0~20 cm的农田耕层土壤样品, 按四分法取1.0 kg为土壤样品, 采样过程中用GPS准确记录采样点位置, 所有土样带回实验室, 经室内自然风干, 剔除植物残体和石块, 研磨后过100 mm筛.土壤样品中的各金属元素在经HCl-HF-HNO3-HClO4四酸消解后, 采用原子吸收光谱法(TAS-990, 北京普析通用仪器有限公司)测定. 为了确保分析结果的可靠性和准确性, 每批样品、每个项目均有3个标准样、3个平行的全程空白同步测定.标样测定结果均在允许误差范围内.

图 1 研究区域及采样点分布示意 Fig. 1 Location of the study area and distribution of sampling sites

1.3 UNMIX模型

UNMIX模型的分析方式是一种特征值分析方式, 是一种基于主成分分析的一种多元受体模型, 它通过自模式曲线分辨技术的转换方式, 确保源分析结果的可靠性[35]. UNMIX模型的运算建立在不同污染源对于受体的源贡献是各个源组分的线性组合; 源中组分对受体的贡献为正值; 样品中有一些源贡献很少甚至没有[36].其计算原理可用以下公式表示[23]:

式中, Cij为第i个样品的第j个物种(j=1, …, n)的含量; Fjk为第j个物种在源k(k=1, …, m)中的质量分数, 代表源的组成; Sik代表源k在第i个样品中的总量, 也就是代表源的贡献率大小; E为分析过程的不确定度或者各个源组成的标准偏差.

UNMIX模型运行操作简单, 可用直接读取格式符合的Excel数据文件.导入数据后, UNMIX模型能够利用自带的数据分析功能, 分析数据的统计特征, 在数据处理窗口(data processing)中可以通过(view/edit points and observations)功能观察含量的边缘分布图或者通过(suggest exclusion)功能对不符合的数据进行剔除, 然后通过选择初始物种(select initial species)以及建议附加物种(suggest additional species)功能, 选择最优于UNMIX模型运算的物种.模型解析出的源成分谱需符合可以用该模型进行解释的最低系统要求(拟合相关系数Min Rsq>0.8, 信噪比Min Sig./Noise>2).

1.4 数据的标准化处理

通常情况下, 土壤中不同元素之间含量差异较大, 为了消除这一差异对结果分析的影响, 在将分析所得的数据导入到UNMIX软件进行分析前, 采用离差标准化的形式对数据进行标准化处理, 经过处理后, 数据均无量纲且观察值的数值范围处于0~1之间.离差标准化公式如下:

式中, Xk(k=1, …, m)为进行离差标准化后的数值, Xi为样品初始分析值, Min Xi为最小分析值, Max Xi为最大分析值.

1.5 数据处理

地统计学空间分布图的制作通过ArcGIS(version 9.3)来完成; 土壤金属来源解析由UNMIX 6.0完成; 其他数据处理由Excel 2010完成.

2 结果与分析 2.1 农田土壤金属的统计特征

铅锌矿区及周边42个农田土壤采样点中的14种元素含量的统计特征见表 1.从中可知:农田土壤中Pb、Zn、Cu、Ni、Mn、Cr、Ti、V和Cd的平均含量全部超过了云南省的土壤背景值, 其中Pb、Zn和Cd分别是当地土壤背景值的27.5、25.3和61.3倍, 并且分别超过国家土壤环境质量标准二级标准(GB 15618-1995)的4.5、11.4和44.6倍, 说明这3种元素在研究区农田土壤中明显富集且污染严重. Pb、Zn和Cd这3种重金属元素的变异系数分别为78.2%、91%和80.6%, 均超过50%, 说明Pb、Zn和Cd在研究区域农田土壤中污染程度差别较大, 污染严重的区域应加强治理.

表 1 土壤中不同元素含量的描述性样品统计特征 Table 1 Statistical characteristics of the selected elements in soils

2.2 农田土壤重金属的UNMIX源解析结果

利用UNMIX 6.0软件, 得出42个土壤采样点14种元素在可解析的多个源中的含量分配情况.其中3个源的Min Rsq=0.83, 代表 83%的物种方差可由该模型解释, 大于系统要求的最小值(Min Rsq>0.8);Min Sig./Noise=3.32, 大于系统要求的最低值(Min Sig./Noise>2), 由此可知这3个源得出的解析结果是可信的.

2.2.1 源成分分析

表 2可知, UNMIX解析出土壤重金属的3个主要来源.在源1的源成分谱中K、Cs和Ba含量较高, 同时也高于其余两个源, 表明源1对K、Cs和Ba元素含量贡献占主导作用.自然条件下, 土壤中全钾和有效钾的含量较低, 而研究区域的钾元素含量的平均值并未超过云南的背景值, 说明钾并未明显富集, 从含量范围来看, 超过背景值的个别点位可能是农民在农田大量使用钾肥的过程中造成的.而结合当地实际情况, 研究区所在的会泽县者海镇是以冶金、建材等产业为主的特色工业城镇, Cs和Ba元素超标则是源于工业污染造成的.因此, 源1代表工业活动引起的污染源.

表 2 土壤重金属三源源成分谱 Table 2 Three source components in heavy metal soils

在源2的源成分谱中, V、Ti、Mn、Ni、Cu和Cr元素的含量高于其他元素, 表明源2主要对这6种元素有贡献.土壤中的V、Ti含量都很低, 主要来源于人类活动的排放, 如化石燃料的燃烧和居民燃煤及焚烧垃圾等, 通过大气沉降造成污染[38].土壤中的Ni和Cr也主要源于上述人类活动.也有报道称化肥中Cr含量比较高, 如磷肥中Cr含量在几十到252.0 mg·kg-1[39].含铜矿的开采和冶炼厂三废的排放、含铜农业化学物质(含铜杀真菌剂和化肥)和有机肥(污泥、猪粪、厩肥和堆肥)的施用可使农田土壤含铜量达到原始土壤的几倍乃至几十倍[40].因此, 源2认为是燃煤和肥料施用所导致的污染源.

在源3的源成分谱中, Pb、Zn、Cd和Ca含量很高, 表明源3对这些金属有着重要的贡献.者海镇有大型的铅锌矿区及金属冶炼厂, 矿产资源的开发及金属冶炼, 会有大量的废水、废气、废渣排放进入环境, 使土壤中的Pb、Zn和Cd严重超标. Pb、Zn和Cd的平均值分别是当地土壤背景值的27.5、25.3和61.3倍, 说明该3种元素在土壤中富集并污染严重.同时, 研究区域从者海镇到矿山镇, 中间有县道Z064连接, 会有大量车辆往来.轮胎摩损会产生大量吸附Zn等重金属的颗粒物进入土壤造成污染[41].地壳中Ca的含量也很高, 大多来源于石灰岩和石膏等.一般认为岩石风化对Ca的贡献为自然作用源.因此, 认为源3是矿山开采活动导致的人为污染源及土壤母质造成的自然污染源的综合污染源.

2.2.2 源贡献分析

图 2可知:源2所占比例最大, 贡献率为68.26%, 表明工业污染、燃煤以及施用肥料等活动是区域土壤重金属污染的主要来源; 源1和源3的贡献率相对源2较小, 说明自然污染源、矿产开采活动对区域土壤重金属元素的贡献较弱.按照常理, 矿区开采及活动应该承担该地区主要的污染源贡献, 但使用该模型解析出来的该部分所占的源贡献率仅为15.42%, 为三源贡献率最小值.可能的原因是会泽县当地对矿产资源开采所使用的方法越来越科学, 管理越来越严格, 对环境造成的污染则随之减少.而工业活动、燃煤及肥料施用是该地区主要的污染源, 也是今后治理的重点.

图 2 土壤重金属不同污染源贡献率 Fig. 2 Contribution rate of different sources in heavy metal soils

各采样点对3个源贡献值情况见图 3.从中可知:源1中2、5~10、12、14~17、19~20、22~23、26、38、41~42号采样点对该源的贡献比较大, 这些采样点主要集中在居民居住地及工业密集区, 说明这些采样点都是工业活动比较频繁的地区, 进一步验证源1是工业源.

图 3 各采样点对不同污染源的贡献 Fig. 3 Contributions of each sampling point to different pollution sources

源2中4~5、10、12、18、20~21、23~29、31~34、36、38~39、41~42号采样点对该源的贡献比较大, 这些采样点也是主要集中在居民居住地和工业密集区, 说明这些采样点都是人类活动比较频繁的地区, 进一步验证源2是燃煤和肥料施用所导致的污染源.

源3中2、14~15、17、25~31、38、40号采样点对该源的贡献比较大, 这些采样点有沿道路分布, 并且大部分集中在居民居住地及矿区周边, 说明这些采样点大部分是人类活动比较频繁的地区, 而矿产资源开采和运输是其主要活动, 进一步说明源3主要是矿山开采活动导致的污染源.

表 3可知, 几个源的相关系数都不为零, 说明源成分之间是相互关联的, 在成因上存在一定的联系, 这也是符合实际结果的; 源1和源2的相关系数较大, 大约在0.5左右, 存在中等程度的相关, 说明源1对源2是有影响的.源1和源3源以及源2和源3之间的相关系数在0.1左右, 说明两两之间没有影响.综上, 研究区农田土壤重金属污染是多种因素综合作用的结果.

表 3 不同土壤重金属污染源之间的相关性 Table 3 Correlations among three different sources of heavy metals in soils

2.3 土壤中6种污染严重的重金属的空间分布

根据记录的采样点坐标位置, 结合每个对应点土壤中重金属的含量, 选取了污染严重的Pb、Zn、Cu、Mn、Ni和Cd这6种重金属, 利用ArcGIS中的地统计分析功能, 进行反距离加权差值分析, 得到矿区周边农田土壤重金属含量的空间分布(图 4).

图 4 研究区土壤6种重金属含量的空间分布 Fig. 4 Spatial distributions of the concentrations of selected heavy metals in the study area

图 4中可以看出, 受矿区渣堆及附近冶炼厂的影响, 采样区西南侧靠近矿渣堆附近的农田土壤Pb、Zn、Cu、Mn、Ni和Cd含量均明显高于其他区域, 并且该区域也是者海镇的主城区, 人口最密集的区域, 说明人类活动是该地区主要的污染源, 与源解析的结果一致; 且研究区域的东北侧附近的农田土壤6种重金属的含量也明显高于中部周边农田, 且研究区域的东北侧也是矿山镇的主城区, 再次验证了人为源对研究区域土壤重金属污染的贡献.从图 4还可以看出, 研究区域的西南角(者海镇)到东北角(矿山镇), 有县道Z064相连, 使得沿着道路附近的农田重金属含量要比远离道路的农田重金属含量要高, 这与源解析结果的源3里面的交通运输导致的人为源一致, 进一步验证了源解析结果的可靠性.

3 讨论

研究区域农田土壤中Pb、Zn和Cd高度积累, 且呈现西南偏高东北偏低的趋势, 与该区域的人类活动方式密切相关.而研究区域的西南部是历史悠久的铅锌矿区, 长时间的矿产开发是导致该地区Pb、Zn和Cd这3种重金属元素高度积累的主要原因.但是根据源贡献分析, 该部分作为源3的主要污染来源只占了全部污染源贡献率的15.42%, 说明通过UNMIX模型解析出来的源3并不是该区域农田土壤重金属的主要污染来源.董騄睿等[42]测定了表层土壤、剖面土壤、肥料及大气降尘中重金属含量, 并通过主成分分析的方法对其进行来源解析, 认为农业施肥对土壤重金属的累积具有非常重要的贡献.与UNMIX模型解析出的源2中肥料施用是该源污染来源的主要贡献之一的结果相吻合.

UNMIX模型作为广泛应用于大气污染物源解析的模型, 目前应用在土壤[31, 36]或者沉积物[43]中的研究还相对较少.本研究的结果表明, UNMIX模型可以很好地应用于土壤重金属污染的源解析.但每种源解析方法都会因为自身的特点而存在一定的局限性, 并非解析出来的结果都会与实际情况完全吻合.因此, 综合应用多种源解析模型和手段, 对不同污染物源解析结果进行综合分析和比较, 这样可以减少单一模型解析结果的局限性, 更好地帮助理解源解析结果并使得源解析的结果更加可靠[44~46].

4 结论

(1) 研究区农田土壤中, Pb、Zn和Cd明显富集, 平均含量分别超过国家土壤环境质量标准的4.5、11.4和44.6倍, 处于严重污染水平.

(2) 利用UNMIX模型解析出3个污染源:源1为工业活动污染源, 源贡献率16.32%;源2为燃煤和肥料施用所导致的污染, 源贡献率68.26%;源3为矿山开采导致的人为污染源及土壤母质造成的自然污染源的综合污染源, 源贡献率15.42%.

(3) 利用ArcGIS中的地统计分析功能, 反距离加权差值分析得到的重金属含量空间分布图与UNMIX模型解析的结果及研究区域的土地利用类型和实际情况基本吻合.表明UNMIX模型是一种有效的农田土壤重金属源解析方法.

参考文献
[1] 李静, 俞天明, 周洁, 等. 铅锌矿区及周边土壤铅、锌、镉、铜的污染健康风险评价[J]. 环境科学, 2008, 29(8): 2327-2330.
Li J, Yu T M, Zhou J, et al. Assessment of health risk for mined soils based on critical thresholds for lead, zinc, cadmium and copper[J]. Chinese Journal of Environmental Science, 2008, 29(8): 2327-2330.
[2] 袁永强, 刘丛强. 广西某地金属冶炼废水外溢对农田土壤的污染特征[J]. 环境科学, 2011, 32(11): 3312-3317.
Yuan Y Q, Liu C Q. Pollution of agricultural soils by a wastewater outflow from a metal smelter in Guangxi Zhuang Autonomous Region[J]. Chinese Journal of Environmental Science, 2011, 32(11): 3312-3317.
[3] 李倩, 秦飞, 季宏兵, 等. 北京市密云水库上游金矿区土壤重金属含量、来源及污染评价[J]. 农业环境科学学报, 2013, 32(12): 2384-2394.
Li Q, Qin F, Ji H B, et al. Contents, sources and contamination assessment of soil heavy metals in gold mine area of upstream part of Miyun Reservoir, Beijing, China[J]. Journal of Agro-Environment Science, 2013, 32(12): 2384-2394. DOI:10.11654/jaes.2013.12.011
[4] 黄兴星, 朱先芳, 唐磊, 等. 密云水库上游某铁矿区土壤重金属含量及形态研究[J]. 中国环境科学, 2012, 32(9): 1632-1639.
Huang X X, Zhu X F, Tang L, et al. Studies on the distribution and chemical speciation of heavy metals in a iron mine soil of the upstream area of Miyun Reservoir, Beijing[J]. China Environmental Science, 2012, 32(9): 1632-1639.
[5] 赵永红, 张静, 周丹, 等. 赣南某钨矿区土壤重金属污染状况研究[J]. 中国环境科学, 2015, 35(8): 2477-2484.
Zhao Y H, Zhang J, Zhou D, et al. Mixed heavy metals contamination of tungsten mine area soil in south of Jiangxi Province[J]. China Environmental Science, 2015, 35(8): 2477-2484.
[6] 王其枫, 王富华, 孙芳芳, 等. 广东韶关主要矿区周边农田土壤铅、镉的形态分布及生物有效性研究[J]. 农业环境科学学报, 2012, 31(6): 1097-1103.
Wang Q F, Wang F H, Sun F F, et al. Fractionation and bioavailability of Pb and Cd in agricultural soils around mining area in Shaoguan Guangdong Province, China[J]. Journal of Agro-Environment Science, 2012, 31(6): 1097-1103.
[7] 叶宏萌, 袁旭音, 赵静. 铜陵矿区河流沉积物重金属的迁移及环境效应[J]. 中国环境科学, 2012, 32(10): 1853-1859.
Ye H M, Yuan X Y, Zhao J. Spatial migration and environmental effects of heavy metals in river sediments from in the Tongling mining area, Anhui Province[J]. China Environmental Science, 2012, 32(10): 1853-1859. DOI:10.3969/j.issn.1000-6923.2012.10.019
[8] 任伯帜, 马宏璞, 郑谐, 等. 锰矿区雨水径流中重金属含量及污染水平的空间结构特征[J]. 环境科学学报, 2014, 34(7): 1730-1737.
Ren B Z, Ma H P, Zheng X, et al. Spatial structure characteristics of heavy metal content and pollution levels of rainfall runoff in manganese ore zone[J]. Acta Scientiae Circumstantiae, 2014, 34(7): 1730-1737.
[9] 鲁荔, 杨金燕, 田丽燕, 等. 大邑铅锌矿区土壤和蔬菜重金属污染现状及评价[J]. 生态与农村环境学报, 2014, 30(3): 374-380.
Lu L, Yang J Y, Tian L Y, et al. Status quo and evaluation of heavy metal pollution of vegetables and soils in Dayi Lead-Zinc Mining Area, China[J]. Journal of Ecology and Rural Environment, 2014, 30(3): 374-380.
[10] 于君宝, 刘景双, 王金达, 等. 矿山复垦土壤典型元素时空变化研究[J]. 中国环境科学, 2001, 21(3): 235-239.
Yu J B, Liu J S, Wang J D, et al. The space-time variation of the typical elements content in covering soil of coal mine reclamation area[J]. China Environmental Science, 2001, 21(3): 235-239.
[11] 刘灵飞, 龙健, 万洪富, 等. 贵州喀斯特山区锑冶炼厂对农业土壤污染特征的影响及风险评价[J]. 土壤, 2013, 45(6): 1036-1047.
Liu L F, Long J, Wan H F, et al. Distribution characteristics and risk assessment of heavy metals in agricultural soils in an abandoned antimony smelter in Guizhou Karst Areas[J]. Soils, 2013, 45(6): 1036-1047.
[12] 王莹, 赵全利, 胡莹, 等. 上虞某铅锌矿区周边土壤植物重金属含量及其污染评价[J]. 环境化学, 2011, 30(7): 1354-1360.
Wang Y, Zhao Q L, Hu Y, et al. Survey and contamination assessment of heavy metals in soil and plants around the Pb/Zn mine in Shangyu, Zhejiang Province[J]. Environmental Chemistry, 2011, 30(7): 1354-1360.
[13] 李良忠, 胡国成, 张丽娟, 等. 矿区家庭灰尘中重金属污染及其潜在生态风险[J]. 中国环境科学, 2015, 35(4): 1230-1238.
Li L Z, Hu G C, Zhang L J, et al. The pollution and potential ecological risk assessment of the heavy metals in household dusts from mineral areas[J]. China Environmental Science, 2015, 35(4): 1230-1238.
[14] Gan J, Lee S J, Liu W P, et al. Distribution and persistence of pyrethroids in runoff sediments[J]. Journal of Environmental Quality, 2005, 34(3): 836-841. DOI:10.2134/jeq2004.0240
[15] Delgado-Moreno L, Wu L, Gan J. Effect of dissolved organic carbon on sorption of pyrethroids to sediments[J]. Environmental Science & Technology, 2010, 44(22): 8473-8478.
[16] Henry R C, Lewis C W, Hopke P K, et al. Review of receptor model fundamentals[J]. Atmospheric Environment (1967), 1984, 18(8): 1507-1515. DOI:10.1016/0004-6981(84)90375-5
[17] 张长波, 骆永明, 吴龙华. 土壤污染物源解析方法及其应用研究进展[J]. 土壤, 2007, 39(2): 190-195.
Zhang C B, Luo Y M, Wu L H. Methods for source apportionment of soil pollutants and their advances in application to soil environmental research[J]. Soils, 2007, 39(2): 190-195.
[18] Gordon G E. Receptor models[J]. Environmental Science & Technology, 1988, 22(10): 1132-1142.
[19] Shi G L, Li X, Feng F C, et al. Combined source apportionment, using positive matrix factorization-chemical mass balance and principal component analysis/multiple linear regression-chemical mass balance models[J]. Atmospheric Environment, 2009, 43(18): 2929-2937. DOI:10.1016/j.atmosenv.2009.02.054
[20] Shi G L, Zeng F, Li X, et al. Estimated contributions and uncertainties of PCA/MLR-CMB results: source apportionment for synthetic and ambient datasets[J]. Atmospheric Environment, 2011, 45(17): 2811-2819. DOI:10.1016/j.atmosenv.2011.03.007
[21] Paatero P. Least squares formulation of robust non-negative factor analysis[J]. Chemometrics and Intelligent Laboratory Systems, 1997, 37(1): 23-35. DOI:10.1016/S0169-7439(96)00044-5
[22] 董騄睿, 胡文友, 黄标, 等. 基于正定矩阵因子分析模型的城郊农田土壤重金属源解析[J]. 中国环境科学, 2015, 35(7): 2103-2111.
Dong L R, Hu W Y, Huang B, et al. Source appointment of heavy metals in suburban farmland soils based on positive matrix factorization[J]. China Environmental Science, 2015, 35(7): 2103-2111.
[23] Hopke P K. Recent developments in receptor modeling[J]. Journal of Chemometrics, 2003, 17(5): 255-265. DOI:10.1002/(ISSN)1099-128X
[24] Wang T Y, Lu Y L, Chen C L, et al. Perfluorinated compounds in a coastal industrial area of Tianjin, China[J]. Environmental Geochemistry and Health, 2012, 34(3): 301-311. DOI:10.1007/s10653-011-9422-2
[25] Shi G L, Tian Y Z, Guo C S, et al. Sediment-pore water partition of PAH source contributions to the Yellow River using two receptor models[J]. Journal of Soils and Sediments, 2012, 12(7): 1154-1163. DOI:10.1007/s11368-012-0540-y
[26] 王雄军, 赖健清, 鲁艳红, 等. 基于因子分析法研究太原市土壤重金属污染的主要来源[J]. 生态环境, 2008, 17(2): 671-676.
Wang X J, Lai J Q, Lu Y H, et al. Main source of soil heavy metal pollution based on factor analysis in Taiyuan[J]. Ecology and Environment, 2008, 17(2): 671-676.
[27] 杨忠平, 卢文喜, 刘新荣, 等. 长春市城区表层土壤重金属污染来源解析[J]. 城市环境与城市生态, 2009, 22(5): 29-33.
Yang Z P, Lu W X, Liu X R, et al. Sources identification of heavy metals in urban soil of Changchun based on principal component analysis[J]. Urban Environment & Urban Ecology, 2009, 22(5): 29-33.
[28] 李瑞平, 郝英华, 李光德, 等. 泰安市农田土壤重金属污染特征及来源解析[J]. 农业环境科学学报, 2011, 30(10): 2012-2017.
Li R P, Hao Y H, Li G D, et al. Characteristics and sources analysis of soil heavy metal pollution in Taian City, Shandong, China[J]. Journal of Agro-Environment Science, 2011, 30(10): 2012-2017.
[29] Pancras J P, Vedantham R, Landis M S, et al. Application of EPA Unmix and nonparametric wind regression on high time resolution trace elements and speciated mercury in Tampa, Florida aerosol[J]. Environmental Science & Technology, 2011, 45(8): 3511-3518.
[30] Callén M S, López J M, Mastral A M. Influence of organic and inorganic markers in the source apportionment of airborne PM10 in Zaragoza (Spain) by two receptor models[J]. Environmental Science and Pollution Research, 2013, 20(5): 3240-3251. DOI:10.1007/s11356-012-1241-1
[31] Lang Y H, Li G L, Wang X M, et al. Combination of Unmix and PMF receptor model to apportion the potential sources and contributions of PAHs in wetland soils from Jiaozhou Bay, China[J]. Marine Pollution Bulletin, 2015, 90(1-2): 129-134. DOI:10.1016/j.marpolbul.2014.11.009
[32] 冯小琼, 彭康, 凌镇浩, 等. 香港地区2005-2010年VOCs污染来源解析及特征研究[J]. 环境科学学报, 2013, 33(1): 173-180.
Feng X Q, Peng K, Ling Z H, et al. Source apportionments and characteristics of VOCs from 2005 to 2010 in Hong Kong[J]. Acta Scientiae Circumstantiae, 2013, 33(1): 173-180.
[33] 王宏镔, 文传浩, 谭晓勇, 等. 云南会泽铅锌矿矿渣废弃地植被重建初探[J]. 云南环境科学, 1998, 17(2): 43-46.
Wang H B, Wen C H, Tan X Y, et al. An investigation on the pioneer plant community of slag spoils of lead and zinc mine, Huize, Yunnan Province and a primary conception of revegetation[J]. Yunnan Environmental Science, 1998, 17(2): 43-46.
[34] 邝荣禧, 胡文友, 何跃, 等. 便携式X射线荧光光谱法(PXRF)在矿区农田土壤重金属快速检测中的应用研究[J]. 土壤, 2015, 47(3): 589-595.
Kuang R X, Hu W Y, He Y, et al. Application of portable x-ray fluorescence (PXRF) for rapid analysis of heavy metals in agricultural soils around mining area[J]. Soils, 2015, 47(3): 589-595.
[35] Henry R C. Current factor analysis receptor models are ill-posed[J]. Atmospheric Environment (1967), 1987, 21(8): 1815-1820. DOI:10.1016/0004-6981(87)90122-3
[36] 艾建超, 王宁, 杨净. 基于UNMIX模型的夹皮沟金矿区土壤重金属源解析[J]. 环境科学, 2014, 35(9): 3530-3536.
Ai J C, Wang N, Yang J. Source apportionment of soil heavy metals in Jiapigou goldmine based on the UNMIX model[J]. Environmental Science, 2014, 35(9): 3530-3536.
[37] 国家环境保护局, 中国环境监测总站. 中国土壤元素背景值[M]. 北京: 中国环境科学出版社, 1990: 87, 330-496.
[38] 杨金燕, 唐亚, 李廷强, 等. 我国钒资源现状及土壤中钒的生物效应[J]. 土壤通报, 2010, 41(6): 1511-1517.
Yang J Y, Tang Y, Li T Q, et al. Soil biogeochemistry and resources situation of vanadium in China[J]. Chinese Journal of Soil Science, 2010, 41(6): 1511-1517.
[39] 闫湘, 王旭, 李秀英, 等. 我国水溶肥料中重金属含量、来源及安全现状[J]. 植物营养与肥料学报, 2016, 22(1): 8-18.
Yan X, Wang X, Li X Y, et al. Contents, source and safety status of major heavy metals in water-soluble fertilizers in China[J]. Journal of Plant Nutrition and Fertilizer, 2016, 22(1): 8-18. DOI:10.11674/zwyf.14443
[40] 吴涛, 伍钧. 铜、铅单一及其复合污染对鱼腥草吸收累积铜和铅的影响[J]. 土壤通报, 2008, 39(3): 670-674.
Wu T, Wu J. Effects of single and combined pollution by Cu and Pb on the accumulation of Cu and Pb in the Houttuynia cordata Thunb[J]. Chinese Journal of Soil Science, 2008, 39(3): 670-674.
[41] 郭平, 谢忠雷, 李军, 等. 长春市土壤重金属污染特征及其潜在生态风险评价[J]. 地理科学, 2005, 25(1): 108-112.
Guo P, Xie Z L, Li J, et al. Specificity of heavy metal pollution and the ecological hazard in urban soils of Changchun City[J]. Scientia Geographica Sinica, 2005, 25(1): 108-112.
[42] 董騄睿, 胡文友, 黄标, 等. 南京沿江典型蔬菜生产系统土壤重金属异常的源解析[J]. 土壤学报, 2014, 51(6): 1251-1261.
Dong L R, Hu W Y, Huang B, et al. Sources of heavy metals in soils of a typical vegetable production system along Yangtze River in Nanjing[J]. Acta Pedologica Sinica, 2014, 51(6): 1251-1261.
[43] Lang Y H, Yang W. Source apportionment of PAHs using Unmix model for Yantai costal surface sediments, China[J]. Bulletin of Environmental Contamination and Toxicology, 2014, 92(1): 30-35. DOI:10.1007/s00128-013-1164-7
[44] Song Y, Xie S D, Zhang Y H, et al. Source apportionment of PM2.5 in Beijing using principal component analysis/absolute principal component scores and UNMIX[J]. Science of the Total Environment, 2006, 372(1): 278-286. DOI:10.1016/j.scitotenv.2006.08.041
[45] Song Y, Dai W, Shao M, et al. Comparison of receptor models for source apportionment of volatile organic compounds in Beijing, China[J]. Environmental Pollution, 2008, 156(1): 174-183. DOI:10.1016/j.envpol.2007.12.014
[46] Yang B, Zhou L L, Xue N D, et al. Source apportionment of polycyclic aromatic hydrocarbons in soils of Huanghuai Plain, China: comparison of three receptor models[J]. Science of the Total Environment, 2013, 443: 31-39. DOI:10.1016/j.scitotenv.2012.10.094