陳澤華 焦思 余愛華 馮滟媚 湯若禹

摘 要：為分析不同土壤重金屬污染評價方法的適用性，以南京市不同功能區(qū)表層土壤為研究對象，通過PXRF法測定重金屬Cr、As、Ni、Pb、Zn、Cu在不同功能區(qū)土壤中的含量。利用單項污染指數(shù)法、地質(zhì)累積指數(shù)法、內(nèi)梅羅綜合指數(shù)法、污染負(fù)荷指數(shù)法和潛在生態(tài)危害指數(shù)法對其污染水平評估。結(jié)果表明：單項污染指數(shù)法和地質(zhì)累積指數(shù)法的評價結(jié)論完全一致，其污染程度為Cr>Pb>Zn>Cu>Ni>As。污染負(fù)荷綜合指數(shù)法與潛在生態(tài)危害指數(shù)法的評價結(jié)果基本一致，由大到小順序為：工業(yè)區(qū)、休閑區(qū)、交通繁忙區(qū)、農(nóng)田區(qū)。內(nèi)梅羅綜合指數(shù)法與潛在生態(tài)危害指數(shù)法對各個功能區(qū)綜合評價結(jié)果也基本一致，分別為輕微污染和存在輕度生態(tài)風(fēng)險。

關(guān)鍵詞：土壤重金屬;污染;評價方法

Abstract：In order to analyze the applicability of different soil heavy metal pollution assessment methods， the surface soils of different functional areas in Nanjing were used as research objects. The contents of heavy metals Cr， As， Ni， Pb， Zn and Cu in different functional areas were determined by PXRF method. The pollution level was evaluated by the single pollution index method， the geo-accumulation index method， the Nemerow comprehensive index method， the pollution load index method and the potential ecological risk index method. The results show that the evaluation results of the single pollution index method and the geo-accumulation index method were completely consistent， and the pollution degree was Cr>Pb>Zn>Cu>Ni>As. The pollution load comprehensive index method and the potential ecological risk index method were basically consistent， and are industrial area > leisure area > heavy traffic area > farmland area. The Nemerow comprehensive index method and the potential ecological risk index method were also basically consistent with the results of a comprehensive evaluation of each functional area， which were slightly polluted and had slight ecological risks.

Keywords：Soil heavy metal; pollute; evaluation methodology

0 引言

土壤污染已逐漸成為環(huán)境污染的主要問題，其中尤以土壤重金屬污染更為突出。由于土壤其自身的結(jié)構(gòu)、特性以及在生態(tài)環(huán)境中扮演的特殊角色，使得土壤污染與大氣污染、水污染相差甚遠(yuǎn)，也更加復(fù)雜[1]。

土壤重金屬污染是指在人類活動過程中將重金屬轉(zhuǎn)移到土壤中，致使土壤中重金屬含量明顯升高，同時造成生態(tài)環(huán)境惡化。土壤重金屬污染是一種由無機物造成的污染，因而具有普遍性、表層積聚性、隱蔽性和不可逆性[2]。根據(jù)2015年中國耕地地球化學(xué)調(diào)查報告顯示，土壤重金屬中度污染、重度污染或超標(biāo)的點位比例占2.5%，覆蓋面積233萬 hm2;輕微污染、輕度污染或超標(biāo)的點位比例占5.7%，覆蓋面積527萬 hm2[3]。城市土壤中重金屬不同程度的累積，對城市居民的健康和城市的發(fā)展都構(gòu)成了潛在的威脅[4]。

目前，評價土壤重金屬污染的方法很多，每種評價方法都有一定的優(yōu)勢與特點，但同時又有一定的局限性，因此為探析每種土壤重金屬污染評價方法的差異性和適用性，本研究以南京市各功能區(qū)表層土壤為研究對象，以單項污染指數(shù)法、地質(zhì)累積指數(shù)法、內(nèi)梅羅綜合指數(shù)法、污染負(fù)荷指數(shù)法和潛在生態(tài)危害指數(shù)法對土壤重金屬污染進(jìn)行評價，以期為土壤重金屬污染評價方法的選擇及對南京市表層土壤重金屬污染狀態(tài)的了解提供一定依據(jù)和參考。

1 材料與方法

1.1 樣品采集與分析

根據(jù)南京市各行政區(qū)相關(guān)資料以及最新行政區(qū)劃圖，確定本研究采樣的大致目標(biāo)范圍，綜合考慮研究區(qū)域內(nèi)地形起伏、土地實際利用情況、土壤類型及城市功能區(qū)的分布（工業(yè)區(qū)、休閑區(qū)、農(nóng)田區(qū)和交通繁忙區(qū)）等因素，選取不同行政區(qū)中的不同功能區(qū)進(jìn)行土壤重金屬污染狀況調(diào)查，其中工業(yè)區(qū)3個，選取了南京煉油廠、南京汽輪機廠和新城科技園;休閑區(qū)5個，選取了棲霞山風(fēng)景區(qū)、紫金山風(fēng)景區(qū)、清涼山公園、武定門公園和河西奧體中心休閑區(qū);交通繁忙區(qū)6個，選取了仙林中心、新街口、新莊立交、古平崗立交、雙橋門立交和集慶門大街;農(nóng)田區(qū)1個，選取八卦洲上的農(nóng)田。所選取的功能區(qū)分別位于南京市主要5個城區(qū)，盡可能做到平均分布且具有一定的代表性。同時根據(jù)《土壤環(huán)境質(zhì)量 建設(shè)用地土壤污染風(fēng)險管控標(biāo)準(zhǔn)（試行）》（GB 36600—2018）和《土壤環(huán)境質(zhì)量 農(nóng)用地土壤污染風(fēng)險管控標(biāo)準(zhǔn)（試行）》（GB 15618—2018）中控制的6種重金屬砷（As）、鉛（Pb）、銅（Cu）、鋅（Zn）、鎳（Ni）和鉻（Cr） 為分析元素。

通過網(wǎng)格布點法確定采樣點和取樣數(shù)，使取樣點分布相對均勻，采樣點盡可能遍布整個功能區(qū)。采樣點通過GPS精確定位，土壤采樣集中在0～20 cm的深度。利用美國Innov公司生產(chǎn)的TRIME-INNOV-X土壤重金屬探測儀原位測定各功能區(qū)表層（0～20 cm）土壤中Cr、Zn、Cu、Pb、Ni、As的含量，在同一個采樣點儀器測量時間為30 s，重復(fù)測量3次，取平均值。

土壤pH采用電位法測定[5]，土壤重金屬的平均值、變異系數(shù)等描述性統(tǒng)計分析采用Excel 2013和SPSS 10.0軟件計算。

1.2 評價方法

1.2.1 單項污染指數(shù)法

單項污染指數(shù)法是針對某一污染物的污染程度進(jìn)行評價的方法，是其他環(huán)境質(zhì)量評價方法和質(zhì)量等級評價的基礎(chǔ)[6]，可以直觀地反映出該區(qū)域中主要的污染因子和其他單個重金屬的污染程度。

單項污染指數(shù)法計算公式為：

1.2.2 內(nèi)梅羅綜合指數(shù)法

單項污染指數(shù)法只能反映不同種類的單個污染物的污染程度，不能完全和充分地反映土壤整體的污染程度。因此就需要將各個重金屬元素的污染指數(shù)結(jié)合起來進(jìn)行評價，內(nèi)梅羅綜合指數(shù)法是目前被廣泛使用的一種綜合指數(shù)法[7]。

內(nèi)梅羅綜合指數(shù)法計算公式為：

1.2.3 地質(zhì)累積指數(shù)法

地質(zhì)累積指數(shù)法也被稱為Muller指數(shù)，是由德國海德堡大學(xué)沉積物研究所的科學(xué)家Muller[8]在20世紀(jì)60年代首次提出，最初它是一種研究水環(huán)境沉積物中重金屬污染程度的定量指標(biāo)。在評估土壤中重金屬污染時除了考慮到人為污染的因素和土壤環(huán)境背景值，還應(yīng)當(dāng)考慮到自然成巖作用是否會引起土壤背景值的變動。地質(zhì)累積指數(shù)法公式為：

在研究某一區(qū)域內(nèi)土壤中重金屬元素的地質(zhì)累積指數(shù)時，應(yīng)該用這個區(qū)域內(nèi)土壤自身背景值作為計算地質(zhì)累積指數(shù)的地球化學(xué)背景值，本研究采用南京市土壤重金屬背景值[9]作為參比，見表3。重金屬地質(zhì)累積指數(shù)（Igeo） 分級與污染程度的關(guān)系，見表4。

1.2.4 污染負(fù)荷指數(shù)法

污染負(fù)荷指數(shù)法是Tomlinson等提出來的一種針對重金屬污染分級的評價方法[10]，該方法被廣泛應(yīng)用于土壤和河流沉積物重金屬污染的評價。

污染負(fù)荷指數(shù)法評價模式如下。

首先計算單個元素的最高污染系數(shù)（CFi）公式為：

1.2.5 潛在生態(tài)危害指數(shù)法

潛在生態(tài)危害指數(shù)法由瑞典科學(xué)家Hakanson提出[11]，本方法根據(jù)沉積學(xué)原理同重金屬相關(guān)性質(zhì)結(jié)合，通過沉積轉(zhuǎn)換規(guī)律建立出土壤或沉積物中重金屬污染評價的模型。潛在生態(tài)危害指數(shù)RI（公式中用RI表示）計算方法公式為：

2 結(jié)果分析

2.1 不同功能區(qū)重金屬含量及數(shù)量統(tǒng)計

由于土壤重金屬污染程度與不同地區(qū)相關(guān)功能存在著密切的聯(lián)系，不同功能區(qū)對于土壤重金屬的貢獻(xiàn)也不盡相同，因此不同功能區(qū)土壤中重金屬（Cr、As、Ni、Pb、Zn、Cu）的污染情況存在著很大差異。對不同功能區(qū)1 600 多個樣點數(shù)據(jù)進(jìn)行分析，得出各區(qū)域6種重金屬的平均含量，見表8。

以《土壤環(huán)境質(zhì)量 建設(shè)用地土壤污染風(fēng)險管控標(biāo)準(zhǔn)（試行）》（GB 36600—2018）為標(biāo)準(zhǔn)，計算土壤中6種重金屬的單項污染指數(shù)。結(jié)果表明，Cr、Ni、Cu、Zn 、Pb 5種重金屬的單項污染指數(shù)均略大于1，存在輕度污染，As 平均含量低于南京市土壤背景值。這也與徐銘焓等[12]在2014年得到的南京市不同功能區(qū)路域表土重金屬污染評價相接近，除As外的其余5種重金屬為輕度污染。

2.2 單項污染指數(shù)法和地質(zhì)累積指數(shù)法

單項污染指數(shù)與地質(zhì)累積指數(shù)都是對單個重金屬元素的污染狀況進(jìn)行評價。但兩種方法使用的參考背景值不相同，單項污染指數(shù)法參考背景值為本區(qū)域的土壤重金屬背景值，本研究參考南京市土壤重金屬背景值（表3）為參考值，其重金屬單項污染指數(shù)見表10。地質(zhì)累積污染指數(shù)評價方法考慮到自然成巖作用，可能會引起土壤背景值發(fā)生變動，考慮成巖作用后的南京土壤重金屬背景值見表11。根據(jù)表8不同功能區(qū)6種重金屬平均含量，計算得各個樣點地質(zhì)累積指數(shù)（Igeo）值見表12。

單項污染指數(shù)法是針對某一污染物的污染程度進(jìn)行評價的方法，也是其他環(huán)境質(zhì)量評價方法和質(zhì)量等級評價的基礎(chǔ)。因此在綜合評價方法中一般會采用較為客觀的單項污染指數(shù)作為基礎(chǔ)。由表10可知，4種不同功能區(qū)的重金屬元素單項污染程度排列為：Cr>Pb>Zn>Cu>Ni>As。由表12可知，不同功能區(qū)6種重金屬元素的地質(zhì)累積指數(shù)排列為：Cr>Pb>Zn>Cu>Ni>As，與單項污染指數(shù)排列順序一致。

2.3 內(nèi)梅羅綜合指數(shù)法、污染負(fù)荷指數(shù)法和潛在生態(tài)危害指數(shù)法

在對單個重金屬元素污染評價結(jié)果的基礎(chǔ)上，為了能夠初步研究并得出不同功能區(qū)整體的土壤重金屬污染程度及分布特征，需要采用綜合的評價方法進(jìn)行評價。本研究分別采用內(nèi)梅羅綜合指數(shù)法、污染負(fù)荷指數(shù)法和潛在生態(tài)危害指數(shù)法進(jìn)行評價，其評價結(jié)果分別見表13—表16。同時對3種方法進(jìn)行比較，見表17。

內(nèi)梅羅綜合指數(shù)均處在輕度污染分界值1.0至重度污染分界值3.0以內(nèi)，污染等級為輕度污染。按照不同功能劃分后土壤重金屬內(nèi)梅羅綜合指數(shù)從大到小排列依次為：休閑區(qū)、工業(yè)區(qū)、農(nóng)田區(qū)、交通繁忙區(qū)。

污染負(fù)荷指數(shù)表明南京市不同功能重金屬的污染狀況整體上較嚴(yán)重，污染等級為Ⅰ，污染負(fù)荷指數(shù)為1.26，屬于中等污染。按照不同功能劃分后土壤重金屬污染負(fù)荷指數(shù)從大到小排列依次為：工業(yè)區(qū)、休閑區(qū)、交通繁忙區(qū)、農(nóng)田區(qū)。

不同重金屬元素的潛在生態(tài)風(fēng)險程度由于考慮到了重金屬毒性響應(yīng)系數(shù)，因此其順序為：As>Pb>Cu>Ni>Cr>Zn;與單項指數(shù)法和地質(zhì)累積指數(shù)法結(jié)果存在差異。南京市不同功能區(qū)所有樣點中重金屬生態(tài)風(fēng)險系數(shù)全部為Ⅰ級（Ei<40），綜合生態(tài)風(fēng)險指數(shù)RI 為31.72，為輕微生態(tài)風(fēng)險水平。其中，As 對研究區(qū)RI 貢獻(xiàn)率最大，應(yīng)引起重視。4種不同功能區(qū)中根據(jù)潛在生態(tài)危害指數(shù)RI數(shù)值全部小于150，對生態(tài)的危害全部屬于輕微危害。按照不同功能劃分后潛在生態(tài)危害指數(shù)從大到小排列依次為：工業(yè)區(qū)、休閑區(qū)、交通繁忙區(qū)、農(nóng)田區(qū)。

綜上所述，單項污染指數(shù)和地質(zhì)累積指數(shù)的評價結(jié)果在不同功能區(qū)的重金屬元素排列順序一致，其中單項污染指數(shù)得到污染情況相較于地質(zhì)累積指數(shù)更為嚴(yán)重，在使用地質(zhì)累積指數(shù)法評價土壤重金屬時沒有出現(xiàn)重污染及以上污染等級，這與其他城市，例如：贛州市[13]、長春市[14]，土壤重金屬的污染水平相近。

由綜合評價方法結(jié)果可知：工業(yè)區(qū)受到的污染程度較大，休閑區(qū)、農(nóng)田和交通繁忙區(qū)土壤污染程度較小，南京市整體上存在一定程度的土壤重金屬污染，與全國大多省會城市土壤環(huán)境相同[15]。對南京市不同功能區(qū)土壤重金屬污染進(jìn)行污染等級評價時，既要考慮多種重金屬協(xié)同作用效應(yīng)，各個不同功能區(qū)地域的影響，又要考慮不同的評價方法，這樣能更加全面準(zhǔn)確地反映實際污染情況。

3 結(jié)論與討論

通過對南京市不同功能區(qū)表層土壤重金屬污染狀況進(jìn)行分析，進(jìn)一步探討5種評價法的適用性及異同。其結(jié)論如下。

（1）南京市主要城區(qū)土壤重金屬Cr、Ni、Cu、Zn、Pb平均質(zhì)量分?jǐn)?shù)均高于南京市土壤背景值，為背景值的1.04～1.91倍。其中，As不存在污染。

（2）單項指數(shù)法和地質(zhì)累積指數(shù)法評價結(jié)果完全一致，其單項重金屬的污染特征表現(xiàn)為：Cr>Pb>Zn>Cu>Ni>As，其中Cr、Cu、Zn在4種功能區(qū)中處于輕微污染等級，As在4種功能區(qū)中處于無污染等級。

（3）內(nèi)梅羅綜合指數(shù)與污染負(fù)荷指數(shù)得到的結(jié)果存在差異，內(nèi)梅羅綜合評價結(jié)果由高到低順序為：休閑區(qū)、工業(yè)區(qū)、農(nóng)田區(qū)、交通繁忙區(qū)。從評價土壤污染程度的過程中發(fā)現(xiàn)，內(nèi)梅羅綜合指數(shù)法選擇單項污染指數(shù)法的評價結(jié)果可以避免人為主觀因素的影響[16]，這樣不僅能夠較為全面地評價不同功能區(qū)土壤重金屬的污染程度，還能突出不同功能區(qū)中相對于南京市土壤背景值的最高濃度重金屬元素對土壤環(huán)境質(zhì)量的影響，也能削弱由于使用各重金屬的平均濃度值帶來對環(huán)境污染評價的作用。污染負(fù)荷指數(shù)評價結(jié)果由高到低順序為：工業(yè)區(qū)、休閑區(qū)、交通繁忙區(qū)、農(nóng)田區(qū)。由于兩種計算方法的不同，累積計算比累加計算更客觀，同時緩解了最大單項重金屬污染指數(shù)，污染負(fù)荷指數(shù)更為準(zhǔn)確。

（4）污染負(fù)荷指數(shù)法與潛在生態(tài)危害指數(shù)法所得出的評價結(jié)論基本一致，其污染等級由高到低分別為：工業(yè)區(qū)、休閑區(qū)、交通繁忙區(qū)、農(nóng)田區(qū)。潛在生態(tài)危害指數(shù)法相較于前兩種綜合評價法，不僅考慮了單一土壤重金屬含量，而且綜合考慮了多種元素之間的相互作用以及土壤環(huán)境對重金屬污染不同特性等因素，更為全面。

（5）3種重金屬污染綜合指數(shù)評價結(jié)果表明，南京市不同功能區(qū)均處于3種方法最低污染等級?？傮w而言，根據(jù)不同方法得出結(jié)果都指出南京市不同功能區(qū)土壤整體重金屬污染水平較低，尚不存在較大的生態(tài)風(fēng)險。

兩種單因子評價方法即單項污染指數(shù)法和地質(zhì)累積指數(shù)法的結(jié)果沒有因不同的評價尺度而得出不同的結(jié)論。單項污染指數(shù)法的評價結(jié)果指數(shù)與重金屬含量具有線性關(guān)系;而地質(zhì)累積指數(shù)法評價土壤污染程度時，其評價結(jié)果與重金屬的分布情況呈現(xiàn)非線性關(guān)系，計算南京土壤重金屬地質(zhì)累積指數(shù)時，既考慮了自然成巖作用引起土壤背景值的變化，同時考慮到了人為因素的影響，因此地質(zhì)累積指數(shù)法彌補了單項指數(shù)法的不足。為了克服兩種方法自身的局限性對評價結(jié)果的影響，將這兩種污染指數(shù)結(jié)合起來，綜合評定土壤樣本對不同功能區(qū)中單項重金屬指數(shù)的污染等級將更為合適。

由于這兩種方法的側(cè)重不同，評價結(jié)果也略有差別。單項指數(shù)法得到的結(jié)果輕度污染為污染的最輕級別，同時地質(zhì)累積指數(shù)法中無—中污染也為最輕級別，兩者綜合后可用輕微污染來描述，因此將兩種評價方法的結(jié)果綜合起來見表18。

由污染負(fù)荷指數(shù)法與內(nèi)梅羅綜合指數(shù)法得到的結(jié)果對比可知，工業(yè)區(qū)污染指數(shù)在兩種計算方法中得到的結(jié)果分別排在第1和第2，也僅為中等污染;因為在內(nèi)梅羅綜合指數(shù)評價農(nóng)田區(qū)時考慮到污染最大值的作用，且由于Pb在農(nóng)田區(qū)指數(shù)較高，所以農(nóng)田區(qū)污染指數(shù)大于交通繁忙區(qū)污染指數(shù)，污染負(fù)荷指數(shù)法中農(nóng)田區(qū)污染等級最低，但兩種方法都是最低污染等級，為清潔土壤;在兩種方法中休閑區(qū)污染指數(shù)均高于交通繁忙區(qū)污染指數(shù)。內(nèi)梅羅綜合指數(shù)法與其他方法結(jié)果存在一定出入，造成這種現(xiàn)象的原因可能是內(nèi)梅羅指數(shù)法中受到某種重金屬元素最大含量的影響，導(dǎo)致其結(jié)果產(chǎn)生偏差，污染負(fù)荷指數(shù)得到的結(jié)果與內(nèi)梅羅綜合指數(shù)進(jìn)行對比得到結(jié)果存在差異，在污染負(fù)荷指數(shù)法中將所有出現(xiàn)的重金屬元素采用同樣的方法累積計算，累積比累加更為客觀，這在一定程度上有效地避免了單一元素對綜合評價的影響，評價方法不僅考慮了單一土壤重金屬含量，而且綜合考慮了多種元素之間的相互作用，更加客觀現(xiàn)實，消除了內(nèi)梅羅指數(shù)法的部分缺陷。

潛在生態(tài)危害指數(shù)法與污染負(fù)荷指數(shù)法兩種方法得到的結(jié)果相同，但是潛在生態(tài)危害指數(shù)在考慮了土壤重金屬的平均含量的基礎(chǔ)上，將重金屬的生態(tài)效應(yīng)和毒理學(xué)知識聯(lián)系起來，綜合考慮了重金屬的毒性在土壤中危害，對于評價區(qū)域內(nèi)重金屬污染水平，消除區(qū)域差異影響，劃分出重金屬潛在危害的程度有一定價值[17-18]。

單項污染指數(shù)法適用于簡單快速地評價土壤重金屬污染情況，可以直接反映出每一種重金屬元素的污染情況。地質(zhì)累積指數(shù)法則是考慮到成巖作用后，相較于單項污染指數(shù)法更為精細(xì)，現(xiàn)在也已經(jīng)被廣泛地接受。內(nèi)梅羅綜合指數(shù)法作為評價環(huán)境整體綜合方法具有一定的局限性和不足，不適用于某一重金屬污染情況較嚴(yán)重的情形。污染負(fù)荷指數(shù)法通過累積的計算方法，適合應(yīng)用于研究土壤重金屬在時間上和空間上變化。而潛在生態(tài)風(fēng)險指數(shù)法則結(jié)合了重金屬毒理學(xué)，更適合于評價重金屬潛在危害以及對人類的影響。

因此，對于類似于南京的人口密集且功能區(qū)分布復(fù)雜的城市進(jìn)行土壤重金屬污染評價，要根據(jù)實際評價地區(qū)的情況、評價結(jié)果的精度要求以及是否需要評估重金屬對人類的潛在危害等多種因素綜合考慮，再選擇使用不同的評價方法，以使評價結(jié)論更加準(zhǔn)確，同時滿足研究的需要。

【參 考 文 獻(xiàn)】

[1]章立佳.上海城市土壤重金屬空間變異結(jié)構(gòu)和分布特征[D].上海：上海師范大學(xué)，2011.

ZHANG L J. Spatial variability and distribution of heavy metals in urban soils of Shanghai[D]. Shanghai： Shanghai Normal University， 2011.

[2]余愛華，盧秀琳，周舒宇，等.城市不同功能區(qū)土壤重金屬特性分析：以南京市玄武區(qū)為例[J].森林工程，2014，30（6）：27-32.

YU A H， LU X L， ZHOU S Y， et al. Characteristics of heavy metals in the soil of different urban areas-a case study of Xuanwu district in Nanjing[J]. Forest Engineering， 2014， 30（6）： 27-32.

[3]中國地質(zhì)調(diào)查局.中國耕地地球化學(xué)調(diào)查報告[R].北京：中國地質(zhì)調(diào)查局，2015.

China Geological Survey. Report on geochemical investigation of cultivated land in China[R]. Beijing： China Geological Survey， 2015.

[4]SPAHIC M P， SAKAN S， CVETKOVIC， et al. Assessment of contamination， environmental risk， and origin of heavy metals in soils surrounding industrial facilities in Vojvodina， Serbia[J]. Environmental Monitoring and Assessment， 2018， 190（4）： 208.

[5]國家生態(tài)環(huán)境部.HJ 962-2018，土壤pH值的測定 電位法[S].北京：中國標(biāo)準(zhǔn)出版社，2018.

Ministry of ecology and environment of the peoples Republic of China. HJ 962-2018， determination of soil pH by potentiometry[S]. Beijing： China Standards Press， 2018.

[6]趙慧，何博，王鐵宇，等.我國南方典型城市土壤重金屬污染特征及源匯關(guān)系分析[J].環(huán)境科學(xué)學(xué)報，2019，39（7）：2231-2239.

ZHAO H， HE B， WANG T Y， et al. Pollution characteristics of heavy metals and source-sink relationship in typical city of the South China[J]. Acta Scientiae Circumstantiae， 2019， 39（7）： 2231-2239.

[7]柳云龍，章立佳，韓曉非，等.上海城市樣帶土壤重金屬空間變異特征及污染評價[J].環(huán)境科學(xué)，2012，33（2）：599-605.

LIU Y L， ZHANG L J， HAN X F， et al. Spatial variability and pollution assessment of heavy metals in soils of Shanghai urban transect[J]. Environmental Science， 2012， 33（2）： 599-605.

[8]MULLER G. Index of geoaccumulation in sediments of the Rhine river[J]. Geojoumal， 1969， 2（3）：108-118.

[9]中國科學(xué)院土壤背景值協(xié)作組.北京、南京地區(qū)土壤中若干元素的自然背景值[J].土壤學(xué)報，1979，16（4）：319-328.

Soil Background Value Collaboration Group of Chinese Academy of Sciences. The natural background values of some trace elements in the important soil types of Beijing and Nanjing areas[J]. Acta Pedologica Sinica， 1979， 16（4）： 319-328.

[10]范拴喜.土壤重金屬污染與控制[M].北京：中國環(huán)境科學(xué)出版社，2011.

FAN S H. Soil heavy metal pollution and control [M]. Beijing： China Environment Science Press， 2011.

[11]HAKANSON L. An ecological risk index for aquatic pollution control， a sedimentological approach[J]. Water Research， 1980， 14（8）： 975-1001.

[12]徐銘焓， 林濤，余愛華，等.南京市不同功能區(qū)路域表土重金屬污染評價及相關(guān)性分析[J].南京林業(yè)大學(xué)學(xué)報（自然科學(xué)版），2014，38（6）：65-71.