基于地理探測器的土壤重金屬影響因子分析及其污染風(fēng)險(xiǎn)評價(jià)

周偉 ，李麗麗，周旭，石佩琪，黃冬青

1.西南大學(xué)地理科學(xué)學(xué)院，重慶 400715；2.北京大學(xué)城市與環(huán)境學(xué)院，北京 100871；3.中國科學(xué)院地理科學(xué)與資源研究所/資源與環(huán)境信息系統(tǒng)國家重點(diǎn)實(shí)驗(yàn)室，北京 100101；4.重慶交通大學(xué)建筑與城市規(guī)劃學(xué)院，重慶 400074

土壤作為自然地理環(huán)境的重要組成部分，是人類賴以生存的自然環(huán)境和農(nóng)業(yè)生產(chǎn)的重要資源（趙其國等，2015），同時(shí)也為系統(tǒng)內(nèi)部能量流動(dòng)和物質(zhì)循環(huán)提供了基礎(chǔ)（孫晶，2011）。伴隨中國城市化進(jìn)程和工業(yè)化的迅速發(fā)展，導(dǎo)致重金屬在土壤中富集，從而造成了土壤重金屬污染。土壤重金屬在土壤中富集，不僅影響土壤性質(zhì)和肥力，而且通過食物鏈累加等間接威脅人類健康。因此，土壤重金屬污染問題已然成為環(huán)境地球科學(xué)研究熱點(diǎn)（葉琛等，2011），更成為國家政府高度關(guān)注的民生問題。近年來國家發(fā)布一系列政策報(bào)告中均提到土壤污染防治的嚴(yán)峻形勢，尤其2017年發(fā)布的《長江經(jīng)濟(jì)帶生態(tài)環(huán)境保護(hù)規(guī)劃》中提到關(guān)于推進(jìn)重點(diǎn)區(qū)域土壤污染防治，以上海、重慶等為重點(diǎn)區(qū)域，完成土壤污染排，建立土壤污染綜合防治先行區(qū)。2018年重慶市發(fā)布的《重慶市環(huán)境保護(hù)條例》加強(qiáng)土壤污染防治禁止將污染物向土壤環(huán)境轉(zhuǎn)移。

土壤重金屬污染風(fēng)險(xiǎn)評價(jià)和污染來源分析，是土壤重金屬污染防治的重要基礎(chǔ)，國內(nèi)外學(xué)者也開展了廣泛探索。Hossain et al.（2015）利用多源統(tǒng)計(jì)分析方法，對重金屬空間分布和潛在來源進(jìn)行了識(shí)別。有學(xué)者利用地統(tǒng)計(jì)學(xué)、地理信息系統(tǒng)等多種方法進(jìn)行土壤重金屬污染空間分布繪制、風(fēng)險(xiǎn)評價(jià)和污染來源分析（Mamut et al.，2017；Song et al.，2018；Zhang et al.，2018）。也有學(xué)者利用主成分分析和污染負(fù)荷指數(shù)法進(jìn)行土壤重金屬污染來源解析和污染風(fēng)險(xiǎn)評價(jià)（易文利等，2018）。近年來，在重金屬污染來源解析研究中地理探測器開始得到廣泛應(yīng)用，因其具有客觀反應(yīng)地理環(huán)境要素對解釋因子的影響力和空間異質(zhì)性的作用（王勁峰等，2017）。有學(xué)者利用地理探測器揭示變量因子與重金屬污染分布的空間相關(guān)關(guān)系和交互作用（李雨等，2017；肖武等，2018）；齊杏杏等（2019）基于地理探測器定量揭示了影響土壤重金屬含量的主要因素。

重慶市作為長江上游經(jīng)濟(jì)帶的重要核心城市，其土壤污染狀況和土壤環(huán)境質(zhì)量對于保障糧食安全和助推長江上游生態(tài)屏障功能的發(fā)揮具有重要意義。學(xué)者們對重慶市土壤重金屬污染也開展了廣泛研究，如城市綠地（楊杉等，2018）、城市道路（郭靜等，2016），典型農(nóng)業(yè)區(qū)（王金霞等，2018）、茶園（張慧等，2017），礦區(qū)（張延玲等，2018），三峽庫區(qū)消落帶（葉琛等，2011）等，研究集中在土壤重金屬污染狀況和風(fēng)險(xiǎn)評價(jià)。然而針對坡耕地及其土壤重金屬污染影響因素的探索較少。因此，本文以重慶主城區(qū)坡耕地的土壤為研究對象，利用單因子污染指數(shù)法進(jìn)行重金屬污染評價(jià)，基于地理探測器進(jìn)行土壤重金屬污染影響因素和各因子交互分析，為重慶市土壤質(zhì)量評價(jià)和土壤污染修復(fù)提供理論依據(jù)。

1 數(shù)據(jù)與方法

1.1 數(shù)據(jù)來源

1.1.1 野外采樣和土壤重金屬含量測試

采樣區(qū)域?yàn)橹貞c市的江津區(qū)、九龍坡區(qū)、大渡口區(qū)和巴南區(qū)，采樣時(shí)間為2018年10月，設(shè)置了21個(gè)采樣地，每個(gè)樣地設(shè)置3個(gè)采樣點(diǎn)，樣點(diǎn)間間隔1—2 km，最后共獲取63個(gè)土壤樣品。每個(gè)樣點(diǎn)采用五點(diǎn)混合采樣法，采集表層0—20 cm的土樣，每份樣品采集1 kg左右，每個(gè)樣品的取土深度及采集量基本保持一致，采集的樣品用土壤密實(shí)袋密封保存，并做好標(biāo)記帶回。利用相關(guān)GPS儀器對所有采樣點(diǎn)進(jìn)行精準(zhǔn)定位，詳細(xì)記錄區(qū)分采樣點(diǎn)的土壤質(zhì)地、土地利用類型、植被類型等相關(guān)信息，利用相機(jī)對采樣區(qū)域進(jìn)行拍照。研究區(qū)和野外采樣點(diǎn)如圖1。

圖1 研究區(qū)域和采樣點(diǎn)圖Fig.1 Study area and the field sampling sites

將土壤樣品帶回進(jìn)行室內(nèi)自然風(fēng)干，用木棒壓磨去除土壤中的砂石或植物殘?bào)w等雜質(zhì)，研磨并過100目（0.25 mm）的土壤篩，將過篩的每個(gè)土壤樣品分為2份，其中1份用于土壤重金屬元素的測試，另一份密封保存?zhèn)洳?，以防交叉污染。土壤重金屬Cu、Zn、Pb、Cr、Ni元素含量采用電感耦合等離子體質(zhì)譜儀（ICP-MS），Mn采用石墨爐原子吸收分光光度法（GF-AAS），測試實(shí)驗(yàn)是在西南大學(xué)巖溶環(huán)境重慶市重點(diǎn)實(shí)驗(yàn)室進(jìn)行。分析方法的準(zhǔn)確度和精密度采用國家一級土壤標(biāo)準(zhǔn)物質(zhì)（GBW）進(jìn)行檢驗(yàn)。

1.1.2 氣象數(shù)據(jù)

利用中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)（http://cdc.cma.gov.cn）查詢獲得2018年10月重慶市標(biāo)準(zhǔn)氣象站點(diǎn)的平均溫度和相對濕度。利用ArcGIS 10.2空間插值獲得氣象柵格數(shù)據(jù)。

1.1.3 其他影響因子數(shù)據(jù)

土壤質(zhì)地與類型數(shù)據(jù)、土地利用數(shù)據(jù)、地形數(shù)據(jù)、河流與道路數(shù)據(jù)通過查詢中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心（http://www.resdc.cn）獲取相應(yīng)數(shù)據(jù)。

1.2 研究方法

1.2.1 地理探測器

地理探測器是一種基于變量間的空間分異性來探測其內(nèi)在聯(lián)系的空間分析模型（王勁峰等，2017），本文中主要利用地理探測器分析5種重金屬Cu、Zn、Cr、Ni和Pb與11種影響因子的相關(guān)性以及多種影響因子的交互作用。本研究采用地理探測器中的因子探測器、交互作用探測器、風(fēng)險(xiǎn)探測器和生態(tài)探測器全面客觀地分析土壤重金屬污染因素。

因子探測器是探究某個(gè)因子X解釋了目標(biāo)Y的空間分異的能力強(qiáng)弱，用因子的解釋力PD,H進(jìn)行判斷。交互作用探測器是識(shí)別不同影響因子之間的交互作用，即影響因子D1和D2共同作用時(shí)是否會(huì)增加或減弱對因變量Y的解釋力，或這些因子對Y的影響是相互獨(dú)立的。

風(fēng)險(xiǎn)探測器主要用于探測影響因子對土壤重金屬分布是否有顯著的差別，具有風(fēng)險(xiǎn)性，用t統(tǒng)計(jì)量來檢驗(yàn)。生態(tài)探測器是用F統(tǒng)計(jì)量來衡量，即用于比較兩個(gè)影響因子D1和D2對屬性Y的空間分布的影響是否有顯著的差異（王勁峰等，2017）。

1.2.2 單因子指數(shù)法

單因子指數(shù)法可以用來確定重金屬污染物的主次要性及其污染程度。一般以污染指數(shù)Pi來表示，以重金屬含量實(shí)測值和評價(jià)標(biāo)準(zhǔn)相比除去量綱來計(jì)算污染指數(shù)：

式中，Pi為樣品中重金屬元素i的單因子污染指數(shù)；Ci為樣品中重金屬含量實(shí)測值；Si為重金屬元素i的評價(jià)標(biāo)準(zhǔn)。

土壤重金屬污染等級劃分標(biāo)準(zhǔn)參照的是由夏家淇（1996）提出的單因子指數(shù)法分級標(biāo)準(zhǔn)。土壤重金屬污染等級分為4級：污染指數(shù)Pi＜1，污染等級為清潔（1級）；污染指數(shù)1≤Pi＜2，為輕度污染（2級）；污染指數(shù)2≤Pi＜3，為中度污染（3級）；污染指數(shù)Pi≥3，為重度污染（4級）。

2 結(jié)果與分析

2.1 土壤重金屬統(tǒng)計(jì)特征分析

從表1可以看出，研究區(qū)內(nèi)土壤重金屬中Cu、Zn、Cr和Pb的含量均值均高于重慶市土壤背景值（趙麗娟，2008），是背景值的1.3—1.4倍，這說明研究區(qū)土壤已經(jīng)受到了Cu、Zn、Cr和Pb的污染。Ni的含量均值低于重慶市土壤背景值，暫未出現(xiàn)污染現(xiàn)象。

表1 研究區(qū)實(shí)測土壤重金屬含量的統(tǒng)計(jì)特征值Table 1 Statistical characteristic of soil heavy metals in the study area

從土壤中5種重金屬的變異系數(shù)來看，值大小為 Cu＞Pb＞Zn＞Ni＞Cr，其中 Cu 和 Pb 的變異系數(shù)分別為0.57和0.40，屬于高度變異，來源復(fù)雜；Zn、Ni和Cr的變異系數(shù)介于0.2—0.4之間，屬于中等變異。說明這5種重金屬受人為因素影響較明顯，眾多研究也表明，變異系數(shù)與人類活動(dòng)干擾程度呈正比（宋金茜等，2017；張永江等，2015；楊杉等，2018）。

偏度為零表示重金屬元素?cái)?shù)值呈正態(tài)分布，Mcgrath et al.（2004）研究認(rèn)為，自然界表層土壤中微量元素的存在都遵從正態(tài)分布，在母質(zhì)類型相同的土壤中金屬元素分布形態(tài)相同或相近，但人類活動(dòng)會(huì)改變其分布狀態(tài)。其中Cu、Zn、Cr和Pb這4種重金屬元素的偏度出現(xiàn)正值產(chǎn)生重金屬分布形態(tài)的正偏現(xiàn)象，其分布峰為右偏，可能是由于人為活動(dòng)等造成重金屬元素不斷富集的現(xiàn)象，這也與重金屬的變異系數(shù)的特征分析結(jié)果相一致。

峰度是測度數(shù)據(jù)在均值附近的集中程度，表 1中5種重金屬的峰度系數(shù)均大于0，說明其分布的集中程度高于正態(tài)分布，且5種重金屬的集中程度表現(xiàn)為 Cu＞Pb＞Zn＞Cr＞Ni，Cu 的數(shù)據(jù)分布較為集中，Ni的數(shù)據(jù)分布最為分散。

2.2 相關(guān)性分析

從表2中可以看出，研究區(qū)Cu、Zn、Cr、Ni和Pb這5種重金屬之間的相關(guān)關(guān)系較強(qiáng)（采用了63個(gè)土壤樣品數(shù)據(jù)）。除了Ni和Pb在P＜0.05水平上相關(guān)性顯著，Cu、Zn、Cr、Ni、Pb在P＜0.01水平上兩兩之間相關(guān)性極為顯著，容易有復(fù)合污染或同源性。其中Cr與Ni的相關(guān)性系數(shù)最大，為0.86（P＜0.01）；Pb與 Ni的相關(guān)性系數(shù)最小，為 0.30（P＜0.05）；不同重金屬之間的相關(guān)性差異明顯。

表2 不同重金屬元素之間的相關(guān)系數(shù)矩陣Table 2 Correlation coefficient matrix between different heavy metal elements

2.3 重金屬污染評價(jià)

通過單因子污染指數(shù)計(jì)算發(fā)現(xiàn)，研究區(qū)內(nèi)土壤Cu、Zn、Cr和Pb的污染指數(shù)均值大于1，表明土壤中4種重金屬均出現(xiàn)輕度污染；只有Ni平均污染指數(shù)小于1，未出現(xiàn)污染。

從表3中我們可以發(fā)現(xiàn)：Cu和Pb含量的單因子污染指數(shù)在各個(gè)等級中均有分布，雖然大部分分布在清潔與輕度污染之間，仍有近10%的樣點(diǎn)在中度污染與重度污染之間；Zn和Cr含量的單因子污染指數(shù)處于清潔至中度污染之間，其中大部分位于輕度污染階段；Ni在5種重金屬中相較而言最為清潔，Ni含量的樣點(diǎn)處于清潔至輕度污染之間。

表3 土壤重金屬污染單因子指數(shù)評價(jià)結(jié)果Table 3 Evaluation results of single factor index of soil heavy metal pollution

由此可見，Cu和Pb對研究區(qū)土壤環(huán)境質(zhì)量造成的危害最為嚴(yán)重，其原因可能是日常生活用品廢棄物、工業(yè)廢棄物等的堆放掩埋，工業(yè)廢水的亂排放，汽車尾氣等廢氣的不清潔處理等（徐福銀等，2014）。Zn和Cr對研究區(qū)土壤環(huán)境質(zhì)量造成的危害也較為嚴(yán)重，出現(xiàn)Cr污染的部分樣點(diǎn)，分布在巴南區(qū)西北方向和大渡口區(qū)東南方向一帶，原因可能與其周圍有多家化工企業(yè)和礦石資源有關(guān)（肖武等，2018）。Ni對研究區(qū)土壤環(huán)境質(zhì)量造成的影響較小?？偟膩砜矗寥乐亟饘倏蓙碓从谕寥莱赏聊纲|(zhì)、工業(yè)廢水灌溉農(nóng)田、工業(yè)污染物的堆積與排放、含重金屬劣質(zhì)化肥農(nóng)藥的使用和大氣降塵等（吳塹虹等，2007）。

2.4 土壤重金屬污染影響因素的地理探測器分析

2.4.1 因子探測器

利用因子探測器探測11個(gè)環(huán)境影響因子對5種重金屬元素的解釋力PD,H值（表4）。就Pb而言，相對濕度、坡度、土壤類型、土地利用類型、粘粒、距河流距離、距道路距離對土壤中重金屬Pb解釋力較大（PD,H＞0.5）；其中，坡度和距河流距離兩種影響因子對 Pb的解釋力高達(dá) 0.93以上，而且研究區(qū)內(nèi)坡度范圍主要在0°—15°之間，說明坡度對重金屬的影響是極為顯著的，重金屬污染的形成與遷移與地勢有著較大的關(guān)系。這也進(jìn)一步驗(yàn)證了影響因子對重金屬的解釋力PD,H既有一致性，又有差異性（李雨等，2017）。說明坡度和距離河流距離是影響重金屬污染和富集的重要驅(qū)動(dòng)因素，同時(shí)表明河流中重金屬元素富集較多。

表4 不同影響因子對5種土壤重金屬的解釋力PD, H值統(tǒng)計(jì)結(jié)果Table 4 The explanatory power value PD, H of different influence factors on five heavy metals

其余4種重金屬Cu、Zn、Cr、Ni，土壤類型、坡度對其含量的解釋力PD,H值均高于0.652；高程對Zn、Cr、Ni解釋力高達(dá)0.75，說明土壤類型、坡度和高程對其空間分布差異具有重要影響，這與前人研究結(jié)果相一致（陳運(yùn)帷等，2019，齊杏杏等2019），其研究發(fā)現(xiàn)海拔、土壤類型、土壤PH是影響土壤重金屬分布的主要因子。

2.4.2 交互作用探測

環(huán)境中土壤重金屬的空間分布是由多種影響因子共同作用的結(jié)果（李忠義等，2009）。研究中我們利用交互作用探測器探測 11種環(huán)境影響因子對5種重金屬元素空間分布變化的作用。本部分選取污染較嚴(yán)重的Pb為例繪制圖表進(jìn)行3種子探測器分析。

表5顯示了兩種影響因子疊加之后的PD,H值，值越大說明交互作用越強(qiáng)，且非線性增強(qiáng)表示兩種不同影響因子的交互作用大于兩種影響因子獨(dú)立作用的和。模型結(jié)果表明高程與坡度、土地利用類型、土壤類型、距河流距離的交互作用值最大，其中高程與平均氣溫、相對濕度、砂粒、粘粒、粉砂為非線性增強(qiáng)；這說明土壤重金屬 Pb污染源的生成與遷移受地勢和氣候影響較大。這與齊杏杏等（2019）研究結(jié)果類似，其研究發(fā)現(xiàn)對Pb累積量分布差異影響的主要因素是地勢和氣候條件。

表5 不同影響因子對土壤重金屬元素Pb影響的交互作用Table 5 Interaction of different influencing factors on Pb

就Cu而言，高程與其余10個(gè)因子交互作用值均達(dá)到0.81以上。就Zn而言，高程與5個(gè)因子（坡度、砂粒、粘粒、粉砂、距道路距離）的交互作用最強(qiáng)。就Ni而言，高程與6個(gè)因子（坡度、土壤類型、砂粒、粘粒、粉砂、距道路距離），平均氣溫、相對濕度分別與高程、土壤類型、距河流距離、距道路距離交互作用較大。就Cr而言，土壤類型與砂粒、粘粒、粉砂、距河流距離的交互作用，以及平均氣溫、相對濕度分別與高程、土壤類型、距河流距離、距道路距離的交互作用占主導(dǎo)地位。以上結(jié)果表明，高程與其他因子交互作用對于重金屬空間分異具有重要影響，同時(shí)氣候條件和土壤類型也是重金屬累積分布的重要影響因子。

2.4.3 風(fēng)險(xiǎn)探測

風(fēng)險(xiǎn)探測器解釋了各影響因子內(nèi)部不同子區(qū)域不同類別間是否存在顯著性差異（圖 2）。風(fēng)險(xiǎn)探測器是采用顯著性水平為0.05的t檢驗(yàn)，“Y”表示存在顯著性差異，“N”表示不存在顯著性差異（表6）。結(jié)果發(fā)現(xiàn)，土壤重金屬Pb元素與平均氣溫、相對濕度、高程、土地利用類型存在顯著性差異，這說明不同氣候條件、海拔高度和土地利用類型下，Pb元素含量分布差異大。李雨等（2017）研究也發(fā)現(xiàn)平均氣溫、相對濕度對重金屬分布影響高于其他因素。而與坡度、土壤類型、土壤質(zhì)地（砂粒、粘粒、粉砂）、距河流距離、距道路距離無顯著性差異，表明以上因素對Pb的分布差異影響小。

圖2 影響因子空間分布Fig.2 Distribution of influence factors

土壤重金屬Cu與平均氣溫、相對濕度、高程、距河流距離、距道路距離均存在顯著差異性，而與其他 6個(gè)因子間無顯著性差異。Cr、Ni與平均氣溫、相對濕度、高程、距河流距離存在顯著性差異。Zn與高程、土地利用類型存在顯著性差異。表明不同重金屬元素的空間分布對 11個(gè)影響因子的響應(yīng)不同。

2.4.4 生態(tài)探測

生態(tài)探測側(cè)重于比較一種影響因子和另一種影響因子對土壤重金屬的空間分布的影響是否有顯著的差異。采用顯著性水平為0.05的t檢驗(yàn)，若顯著，則記為“Y”；若不顯著，則記為“N”。

生態(tài)探測結(jié)果（表6）表明，高程對Pb的影響與平均氣溫、相對濕度、坡度、土地利用類型、土壤類型、土壤質(zhì)地（砂粒、粘粒、粉砂）、距離河流和道路的距離等影響因子間存在顯著差異，而其他影響因子之間的差異不顯著。說明相比其他8個(gè)因素，地勢、氣候條件和母質(zhì)等對 Pb空間分布異質(zhì)性具有重要影響。

表6 重金屬Pb的生態(tài)探測結(jié)果Table 6 Ecological detection results of heavy metal Pb

對Cu、Cr和Ni，高程、坡度、土地利用類型、土壤類型、土壤質(zhì)地（砂粒、粘粒、粉砂）、距河流距離和距道路距離對其空間分布的影響與其他影響因子間存在顯著差異，而其他影響因子之間的解釋力差異不顯著。對Zn，高程、土壤質(zhì)地（砂粒、粘粒、粉砂）、距河流距離和距道路距離對Zn空間分布的影響與其他影響因子間存在顯著差異。表明地勢、母質(zhì)和人為干擾對以上4種重金屬空間分布有顯著影響，與齊杏杏等（2019）研究結(jié)論類似。

3 結(jié)論

（1）研究區(qū)內(nèi)土壤受到Cu、Zn、Cr和Pb重金屬元素的污染，其平均含量均高于重慶市土壤背景值，是背景值的1.3—1.4倍；從變異系數(shù)和相關(guān)性分析可以看出，5種元素均呈現(xiàn)中等和高度變異且相關(guān)性較強(qiáng)。

（2）單因子污染指數(shù)均值由大到小依次為：Pb、Cu、Zn、Cr、Ni；Cu、Pb為重度污染水平，其余為中度或輕度污染水平。因此Cu、Pb是導(dǎo)致整個(gè)研究區(qū)域土壤處于輕度污染生態(tài)風(fēng)險(xiǎn)的主要因素。

（3）5種重金屬與11種環(huán)境因子的因子探測發(fā)現(xiàn)，高程、坡度和土壤類型對5種土壤重金屬含量的解釋力最顯著。交互作用探測發(fā)現(xiàn)，高程與其他因子交互作用是重金屬空間分布差異的主導(dǎo)因素，氣候條件和土壤類型也是重要影響因子。土壤重金屬空間分布是多種因素共同作用的結(jié)果，主導(dǎo)因子高程、坡度和土壤類型可作為后續(xù)土壤重金屬含量空間預(yù)測模型的輔助變量。