攀西某芒果產(chǎn)區(qū)土壤重金屬污染源解析

王哲，蔣虎，羅瑩，賈文靜，曾秋平，張振龍，張家千，馮喜楊

(1.西南科技大學(xué)環(huán)境與資源學(xué)院，綿陽 621010；2.中國科學(xué)技術(shù)大學(xué)地球與空間科學(xué)學(xué)院，合肥 230000；3.中國科學(xué)院、水利部成都山地災(zāi)害與環(huán)境研究所，成都 610041；4.中國科學(xué)院大學(xué)，北京 100049)

攀枝花地區(qū)由于其氣候溫暖、光照充足、降雨適宜等優(yōu)越自然條件，使得該地區(qū)所產(chǎn)芒果肉質(zhì)細膩，果味醇香，產(chǎn)量高，是中國的主要芒果產(chǎn)地[1]。然而，攀枝花以礦產(chǎn)資源為主導(dǎo)的經(jīng)濟發(fā)展方式和長期的農(nóng)業(yè)生產(chǎn)活動可能會引起土壤中重金屬的富集，影響芒果的質(zhì)量與經(jīng)濟的可持續(xù)發(fā)展。攀枝花地區(qū)礦產(chǎn)資源豐厚，在礦產(chǎn)資源開采、運輸和選冶等生產(chǎn)活動中不可避免地會產(chǎn)生大量的尾礦渣、廢水、廢氣及粉塵，致使重金屬不同程度地富集于土壤[2]，潛在威脅該地區(qū)芒果產(chǎn)業(yè)鏈的持續(xù)性發(fā)展。同時，芒果種植過程中農(nóng)藥、化肥的長期施用也會導(dǎo)致殘留的重金屬富集。隨著重金屬在土壤中富集到一定程度，不僅會使土壤肥力流失，而且會影響農(nóng)產(chǎn)品的質(zhì)量與品質(zhì)，經(jīng)食物鏈被人體攝入，影響身體發(fā)育與健康[3]。

為了厘清芒果產(chǎn)區(qū)的土壤及芒果中重金屬的富集情況，中外學(xué)者開展了大量研究。汪婷等人以攀枝花市鹽邊縣某芒果種植區(qū)為研究區(qū)分析了土壤中重金屬元素的空間分布特征，并評估了研究區(qū)土壤重金屬污染程度，發(fā)現(xiàn)研究區(qū)土壤中Cd 潛在生態(tài)危害較強；楊定清等[4]測定了攀枝花芒果基地土壤和芒果中的鎳含量，發(fā)現(xiàn)西區(qū)和仁和芒果基地土壤鎳含量顯著高于中國其他地區(qū)土壤鎳平均含量，并指出芒果中平均鎳含量與各基地土壤鎳平均含量之間呈現(xiàn)出正相關(guān)趨勢。Chata等[5]調(diào)查了尼日利亞某市場的4個芒果品種的重金屬污染狀況，并測定了芒果中Pb、Cu、Fe、Mn等重金屬含量，發(fā)現(xiàn)僅有Pb超過了世衛(wèi)組織和聯(lián)合國糧食及農(nóng)業(yè)組織的重金屬限值；Rahimzadeh等[6]測定了伊朗不同基因型不同成熟度的芒果栽培樣品中的Pb、Cr、Cd、As等重金屬含量，發(fā)現(xiàn)未成熟芒果中重金屬含量顯著高于熟芒果，而不同基因型的芒果樣品中重金屬殘留量與品種之間的關(guān)系不顯著，并指出所有芒果樣品中重金屬殘留量低于可接受安全性的最大可接受水平。Silva等[7]測定了巴西某芒果種植園土壤和葉片中Cu、Cr、Fe、Zn、Mn、Ni和Pb等重金屬的濃度，發(fā)現(xiàn)芒果栽培易導(dǎo)致土壤表層Cu和Zn積累，表層和次表層Ni、Pb、Mn和Fe的含量降低。Sam-Uket等[8]測定了尼日利亞某金礦區(qū)附近的芒果植物根、莖、葉、果實等部位的重金屬(Pb、Cd、Au和Fe)，發(fā)現(xiàn)芒果植株中的重金屬含量普遍高于允許的限度，認為該礦區(qū)附近的芒果果實不適合人類食用。以上研究以評估芒果產(chǎn)區(qū)重金屬污染現(xiàn)狀為主，表明芒果種植區(qū)周圍地理環(huán)境、人類生產(chǎn)活動、種植方式的不同會影響區(qū)內(nèi)土壤及芒果中重金屬的積累狀況，缺乏芒果產(chǎn)區(qū)的重金屬污染源解析相關(guān)的研究，亟待厘清芒果種植區(qū)土壤重金屬的富集與工農(nóng)業(yè)生產(chǎn)活動之間關(guān)系。因此，現(xiàn)結(jié)合地理信息系統(tǒng)(geographic information system，GIS)空間分析、相關(guān)性分析以及正定矩陣因子分析(PMF源解析法)等方法，以米易縣某芒果種植區(qū)土壤為研究對象聚焦于芒果種植，分析區(qū)內(nèi)土壤重金屬污染現(xiàn)狀及其來源，以期為當(dāng)?shù)卣贫ㄍ寥乐亟饘傩迯?fù)策略提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況與樣品采集

研究區(qū)位于攀枝花市東北部的米易縣芒果種植區(qū)，面積約0.37 km2。在氣候上，屬于干熱河谷立體氣候，平均年降雨1 112.6 mm，最低氣溫-2.4 ℃，最高氣溫4.03 ℃，平均氣溫約為19.7 ℃。在距離研究區(qū)3 km范圍內(nèi)，西方有一個以鐵礦石為主要原料的鐵精粉廠。為了便于芒果運輸與種植，居民建筑區(qū)選址靠近種植區(qū)，區(qū)內(nèi)道路網(wǎng)密集，交通方便，如圖1所示。

為查明芒果種植區(qū)的土壤重金屬污染現(xiàn)狀及來源，采集了27個表層土樣并記錄采樣點位置?？紤]到取樣的地形限制，如陡坎、建筑區(qū)等，依據(jù)土壤監(jiān)測規(guī)范[9]布設(shè)采樣點位置，如圖1所示，采樣點間距100 m左右，每個土壤樣品由采樣點10 m范圍內(nèi)的5處土樣(各400 g)均勻混合，四分后留取對角土樣1 kg后形成并裝入聚乙烯塑料袋密封。

圖1 研究區(qū)位置

1.2 實驗測量

在測定重金屬含量前，土樣需要經(jīng)過3個處理階段，包括除濕、篩分以及消解。土壤樣品在25 ℃的自然條件下風(fēng)干(除濕階段)后進入篩分階段。在篩分階段，將風(fēng)干土壤中的異物去除，使用木棍碾碎并均勻混合，均勻四等分后選取對角位土壤過孔徑為2 mm的尼龍篩并留取200 g土樣，然后進一步研磨并篩分過孔徑為0.15 mm的尼龍篩并為下一階段留取30 g土樣。在消解階段，精確稱取0.150 0 g備用土樣置于消解罐中并加入5 mL硝酸(ρHNO3=1.42 g/mL)，浸潤0.5 h去除有機物。然后，再加入2 mL過氧化氫(ρH2O2=1.49 g/mL)和2 mL的氫氟酸(ρHF=1.49 g/mL)，蓋緊罐蓋，安裝好高溫高壓消解罐后，置于電熱恒溫鼓風(fēng)干燥箱中以180 ℃密閉消解4 h。最后，在冷卻到室溫后，移到比色管中，填充2%的硝酸溶液，直到容量達到25 mL，移送至西南科技大學(xué)分析檢測中心以電感耦合等離子體發(fā)射光譜法測量樣品中重金屬含量，并整理匯總用于數(shù)據(jù)統(tǒng)計與分析。剩余土樣按固液比1∶2.5在燒杯中均勻混合后的用于pH的測定，測定方法為玻璃電極法[10]。

1.3 源解析法

PMF(positive matrix factorization)分析法也稱正定矩陣因子分析法，能夠?qū)⒏呔S變量簡化并轉(zhuǎn)變?yōu)閹讉€綜合因子，無需測定復(fù)雜的原譜，在明確解析源貢獻率方面具有優(yōu)勢[11-12]。通過PMF分析法，將測量的土壤重金屬濃度數(shù)據(jù)對應(yīng)的原始矩陣分解為兩個因子矩陣(源貢獻率矩陣F和源成分譜矩陣G)和一個殘差矩陣D，如式(1)～式(3)所示。

(1)

式(1)中：Xij為第i個樣品第j種元素的濃度，mg/kg；Fik為源k對第i個樣品的貢獻；Gkj為源k中第j個元素的濃度，mg/kg；Dij為殘差矩陣；i=1,2,…,N;j=1,2,…,M。

為了獲得矩陣G和F的最優(yōu)解，PMF通過逐步分解原始矩陣的迭代計算使Q值最小化，計算公式為

(2)

式(2)中：uij為不確定度矩陣中的元素。

在利用PMF分析污染源前，需要輸入濃度數(shù)據(jù)和不確定度數(shù)據(jù)，不確定度由式(3)和式(4)計算得出。

若c≤MDL,有

(3)

若c>MDL,有

(4)

式中：error為相對標準偏差；c為元素濃度，mg/kg；MDL為方法檢出限。

2 結(jié)果與分析

2.1 土壤酸堿性及重金屬含量特征分析

研究表明，芒果適宜種植在土壤pH為5.5～7.5的中性偏酸土壤中，否則芒果的產(chǎn)量與質(zhì)量將會降低[13-14]。并且，不同的土壤pH也可能影響重金屬在土壤中的有效性和毒性。通過室內(nèi)測量，研究區(qū)采集的27個土壤樣品的pH在4.89～8.09，均值為5.95，92.59%的采樣點pH<7，土壤呈現(xiàn)中性偏酸性，適宜芒果的種植。

基于重金屬含量和pH描述性統(tǒng)計結(jié)果(表1)，依據(jù)四川省土壤元素背景值[15]與《土壤環(huán)境質(zhì)量標準》(GB 15618—2018)[16](表2)計算了各重金屬的超標情況，如表3所示。超標率計算結(jié)果顯示：與四川省土壤元素背景值對比，超標率排序為Cd > Cu > Ni > Pb > Cr > Zn；與《土壤環(huán)境質(zhì)量標準》(GB 15618—2018)的污染風(fēng)險篩選值對比，超標率排序為Cd > Cu > Cr > Ni > Pb > Zn。Cd的超標率高達100%，大幅高出其他重金屬的超標率，同時Cd濃度均值與四川省背景值的比值遠遠大于其他重金屬。因此，Cd對研究區(qū)內(nèi)的土壤重金屬污染起控制性作用。As、Hg含量測定顯示，研究區(qū)土壤中As、Hg含量極少，部分樣品的含量甚至因低于儀器檢測限無法測定，能夠檢測出的重金屬含量均低于農(nóng)用地土壤污染風(fēng)險篩選值，故未納入污染源分析。此外，通過對研究區(qū)內(nèi)的芒果果肉處理、分析后發(fā)現(xiàn)果肉內(nèi)中這6種金屬含量遠遠低于《食品安全國家標準食品中污染物限量》(GB2762—2017)[17]關(guān)于新鮮水果或蔬菜中重金屬含量的要求，表明重金屬元素未在芒果果肉中富集，芒果農(nóng)產(chǎn)品是安全的，但其他農(nóng)產(chǎn)品的超標情況待采樣分析。

表2 土壤重金屬背景值與篩選值

表3 重金屬含量點位超標情況

2.2 相關(guān)性分析

相關(guān)性分析是源解析常用并易于實現(xiàn)的方法，相關(guān)性系數(shù)越高暗示重金屬來源一致的可能性越高[18]。相關(guān)性分析顯示：除Pb外，其余重金屬之間均有較好的相關(guān)性；其中，Cr與Ni的皮爾遜相關(guān)系數(shù)在0.95以上，具有很強的同源性。詳細相關(guān)性系數(shù)如表4所示。

表4 相關(guān)性分析結(jié)果

2.3 空間分布特征

基于ArcGIS Desktop平臺[19]，多次調(diào)試Geostatistical Analyst模塊中的數(shù)據(jù)預(yù)處理參數(shù)，在研究區(qū)范圍內(nèi)對各重金屬含量執(zhí)行反距離加權(quán)插值，選取平均標準差最接近0且均方根標準化預(yù)測誤差最接近1的結(jié)果，據(jù)此分析土樣中的重金屬元素含量的空間分布規(guī)律并繪制結(jié)果圖，如圖2所示。

圖2 研究區(qū)土壤重金屬含量空間分布

通過對研究區(qū)土壤重金屬含量空間分布分析結(jié)果可以看出：Cd、Zn、Ni均呈現(xiàn)北、西北和南方濃度較高的特征；Pb呈現(xiàn)出濃度高值區(qū)向車輛流量大的公路附近聚集的現(xiàn)象，并且靠近居民區(qū)處濃度達到最高；Cr、Cu呈現(xiàn)出北、西北方濃度較高，南方較低的特征。

2.4 重金屬源解析

2.4.1 篩除離散值

離散值的存在會影響PMF源解析法計算貢獻率的準確性[20]。因此，有必要在相關(guān)性分析的基礎(chǔ)上對相關(guān)性呈顯著正相關(guān)的重金屬元素組繪制散點圖，篩選離群值，如圖3所示。分析結(jié)果表明M9和M21為需要篩除的離群值點。

圖3 散點圖

2.4.2 重金屬源解析

基于PMF分析模型，通過多次調(diào)節(jié)因子數(shù)使目標函數(shù)Q取得最小值，最終確定5個因子，其對各重金屬的源貢獻率和源成分譜如表5和表6所示。據(jù)此，結(jié)合土壤重金屬空間分布、元素間相關(guān)性、PMF分析結(jié)果等信息分析了研究區(qū)土壤重金屬的來源，具體分析如下。

表5 源貢獻率表

表6 源成分譜表

隨著城市化進程的加快，生活資料越加豐富多樣，影響了土壤中鋅的積累[21]。以往研究顯示，金屬制品或鍍金屬制品等生活垃圾中金屬離子的溶出，以及在使用與廢棄過程中重金屬離子的釋放，致使生活垃圾成為重金屬污染直接來源之一[22]，并且含鋅的生活垃圾是土壤中Zn的重要貢獻者。龍於洋等[23]指出新鮮垃圾填埋層的垃圾體是土壤中Zn的重要輸入來源。安曉雯等[24]研究發(fā)現(xiàn)大連市城市生活垃圾衛(wèi)生填埋場鋅含量嚴重超標，分析認為電池、電器元件等電子生活垃圾和部分工業(yè)垃圾是生活垃圾填埋場鋅污染嚴重的原因之一。楊淑英等[25]通過分析生活垃圾中重金屬元素的污染特性，認為垃圾主要重金屬來源可能是弱酸性食品殘渣對金屬制品的腐蝕和溶解以及清掃的塵土，鋅污染來源于電池和塵土?，F(xiàn)場考察發(fā)現(xiàn)，當(dāng)?shù)赜幸粋€堆積生活垃圾以及廢棄農(nóng)產(chǎn)品的垃圾堆積處。結(jié)合PMF分析結(jié)果，源1對Zn元素的貢獻率最高，高達41.47%，對Zn的貢獻率起主導(dǎo)作用，故推測源1為生活源。

Gmochowska等[26]認為Pb的來源可分自然源和人為來源，分別源于成巖母質(zhì)和采礦、冶煉、化石燃料燃燒等人類生產(chǎn)經(jīng)濟活動。趙多勇等[27]運用二元模型分析大氣降塵和成巖母質(zhì)對耕種區(qū)鉛的貢獻率，成巖母質(zhì)母質(zhì)的貢獻率大于50%。張宇豪等[28]在以往研究的基礎(chǔ)上統(tǒng)計遭受鉛污染土壤的污染源的出現(xiàn)頻率，發(fā)現(xiàn)出現(xiàn)頻率最高的污染源是汽車尾氣和燃煤源。其余研究也表明與機動車尾氣排放和汽車磨損排放相關(guān)的交通源是鉛(Pb)的來源之一[29-30]。研究區(qū)位于農(nóng)村，為了便于芒果運輸，區(qū)內(nèi)交通網(wǎng)絡(luò)密集，呈現(xiàn)出交通網(wǎng)密集程度由居民區(qū)附近向四周逐漸降低的現(xiàn)象，如圖1所示。通過GIS的空間插值模塊對項目區(qū)域內(nèi)土壤中Pb的濃度插值可視化，結(jié)果顯示：Pb濃度的分布情況與交通網(wǎng)密集程度的分布高度一致，如圖2(e)所示。除此之外，PMF分析表明源2對Pb的貢獻率高達32.09%，而對Cd、Cu、Ni和Zn的貢獻率均小于6.5%(Cr為19.24%)，表明源2是對研究區(qū)的主要貢獻重金屬元素為Pb。綜上所述，推測源2為交通源。

研究區(qū)地處攀西地區(qū)，該地區(qū)礦產(chǎn)資源豐富，礦石熔煉量巨大[31]。在金屬冶煉期間如果廢渣處理不當(dāng)，鋅礦、銅礦及其他礦石冶煉等生產(chǎn)活動容易導(dǎo)致鎘、鉻、銅、鎳、鋅等重金屬進入環(huán)境并積累[32-33]。研究區(qū)2公里范圍內(nèi)有1家鐵精粉廠，其生產(chǎn)原料正是富含鐵礦石，位于研究區(qū)西方，如圖1所示。PMF分析表明源3對Cd、Cr、Cu、Ni和Zn都有一定的貢獻率，對Cu的貢獻率高達54.33%。同時，距鐵粉廠越遠，土壤重金屬(Cd、Cr、Cu、Ni、Zn)含量呈現(xiàn)降低趨勢，如圖4所示。故推斷源3為工業(yè)源。

圖4 A點到C點沿線重金屬濃度曲線

PMF分析表明，源4對元素Cd、Pb、Zn的貢獻率均大于25%。根據(jù)表1的描述性統(tǒng)計結(jié)果，Cd、Pb、Zn的變異系數(shù)相對較小，這表明其空間變異程度較小，沒有顯著受到人類活動的影響[34]。因此，推斷源4源于土壤母質(zhì)，屬于自然源。

由于攀西地區(qū)屬南亞熱帶氣候，是中國著名的熱作區(qū)之一，是中國四川省農(nóng)業(yè)資源開發(fā)潛力最大，農(nóng)產(chǎn)品特色最為突出的區(qū)域。尤其在“八五”期間，米易縣是國家重點建設(shè)的農(nóng)業(yè)基地[35]。在農(nóng)產(chǎn)品生產(chǎn)活動中，農(nóng)民追求高產(chǎn)量，為農(nóng)作物施加大量的無機、有機肥料，為土壤中重金屬的富集提供了外來輸入途徑。王美等[36]通過查閱大量文獻，發(fā)現(xiàn)在農(nóng)產(chǎn)品生產(chǎn)過程中，土壤重金屬的富集與磷肥有密切關(guān)系(Cd表現(xiàn)尤其突出)，與氮、鉀肥關(guān)系不大，并且施用畜禽肥等有機肥對土壤重金屬含量的影響大于化肥。吳建繁等[37]指出以磷礦為原料生產(chǎn)的磷含有Cd、Cr等重金屬元素。黃青青等[38]收集并測定中國5個磷礦富集省的磷肥樣品，結(jié)果表明Cd有明顯的超標現(xiàn)象。余垚等[39]指出長期大量施用重金屬含量過高的磷肥有可能造成重金屬累積。Jiao等[40]指出在正常種植活動下施用磷肥對受施土壤中微量元素積累的影響是有限的和局部的。邢素麗等[41]實驗分析發(fā)現(xiàn)施入畜禽糞便有機肥普遍增加了土壤中重金屬Cd、Cr、Cu的含量。Uprety等[42]分析了長期施用有機肥和無機肥對耕地的重金屬微元素的影響，發(fā)現(xiàn)只使用磷肥的農(nóng)田土壤中As、Cd和Cr含量最高，僅使用家禽糞便的農(nóng)田土壤中Cu、Mn和Ni含量最高，僅使用糞水的農(nóng)田土壤中Pb和Zn含量最高，這意味著家禽肥為土壤中Cu、Ni的富集提供了外來輸入來源。根據(jù)現(xiàn)場考察結(jié)果，居民地在研究區(qū)海拔較低處，項目區(qū)需要施肥的面積廣闊，為了提高種植產(chǎn)品的質(zhì)量，農(nóng)民傾向于使用磷肥以及含水率較低易于搬運的家禽糞等有機肥料，導(dǎo)致Cd、Cr、Cu以及Ni進入土壤并積累，與源5對Cd、Cr、Cu以及Ni有更大的貢獻率相呼應(yīng)。因此，推測源5為農(nóng)業(yè)源。

3 結(jié)論

(1)統(tǒng)計結(jié)果顯示，土樣中Cd、Cr、Cu、Ni、Pb和Zn的平均值含量分別為3.64、132.06、67.35、57.35、26.86、63.726 mg/kg，是四川省背景值的21.51、1.44、2.33、1.64、1.03、0.76倍。根據(jù)農(nóng)用地土壤污染風(fēng)險篩選值，Cd、Cr、Cu、Ni、Pb和Zn的點位超標率分別為100.00%、33.33%、37.04%、29.63%、0%、0%，但在芒果的果肉中未檢測到重金屬超標，表明土壤中的重金屬污染未富集到芒果中。

(2)基于PMF分析方法，結(jié)合數(shù)值統(tǒng)計、相關(guān)性分析以及空間插值結(jié)果，識別出研究區(qū)土壤重金屬的5個污染源，包括生活源、交通源、工業(yè)源、自然源以及農(nóng)業(yè)源，為該地區(qū)土壤重金屬污染的修復(fù)提供了參考依據(jù)。

(3)Cd對研究區(qū)內(nèi)的土壤重金屬污染起控制性作用，主要來源于土壤母質(zhì)、農(nóng)業(yè)生產(chǎn)活動，而與人類生產(chǎn)活動相關(guān)的污染源對Cd的總貢獻率是土壤母質(zhì)來源的2.3倍，芒果種植區(qū)土壤重金屬的富集與經(jīng)濟作物的生產(chǎn)活動之間有一定的關(guān)系。

(4)研究區(qū)土壤中Cu主要來源于工農(nóng)業(yè)生產(chǎn)活動，其貢獻率共計97.19%。

(5)研究區(qū)土壤中Cr、Ni主要來源于芒果種植過程中施肥等農(nóng)業(yè)生產(chǎn)活動。