劉海偉，宗 浩，王海云，劉文濤，楊舉田*，石 屹，張平德，劉 莉

?

臨沂植煙土壤重金屬空間分布特征與生態(tài)健康風(fēng)險(xiǎn)評(píng)價(jià)

劉海偉1，宗 浩2，王海云1，劉文濤2，楊舉田2*，石 屹1，張平德2，劉 莉2

（1.中國農(nóng)業(yè)科學(xué)院煙草研究所，農(nóng)業(yè)部煙草生物學(xué)與加工重點(diǎn)實(shí)驗(yàn)室，青島 266101；2.山東臨沂煙草有限公司，山東 臨沂 276000）

土壤重金屬空間分布與風(fēng)險(xiǎn)評(píng)價(jià)對(duì)土壤重金屬污染防治具有重要的指導(dǎo)意義。以傳統(tǒng)適宜烤煙種植區(qū)山東省臨沂市為研究對(duì)象，通過布點(diǎn)采樣和地理信息系統(tǒng)制圖，分析了土壤重金屬的空間分布以及地區(qū)、土類間的差異，采用地積累指數(shù)、污染指數(shù)分析了土壤重金屬富集污染水平，并對(duì)土壤重金屬潛在生態(tài)風(fēng)險(xiǎn)和人體健康風(fēng)險(xiǎn)進(jìn)行了評(píng)價(jià)。結(jié)果表明，臨沂植煙耕層土壤重金屬As、Cd、Cr、Cu、Hg、Ni、Pb和Zn平均含量分別為3.85、0.09、47.56、14.45、0.07、18.04、17.82和29.21 mg/kg；空間分布上，北部、中部和西南地區(qū)土壤重金屬含量相對(duì)較高；8種重金屬中僅Cd和Hg在部分土壤中表現(xiàn)出富集和污染；土壤重金屬潛在生態(tài)風(fēng)險(xiǎn)大小順序?yàn)镠g>Cd>As>Ni>Pb>Cu>Cr>Zn，其中部分土壤Cd和Hg對(duì)當(dāng)?shù)丨h(huán)境具有顯著潛在生態(tài)風(fēng)險(xiǎn)；而土壤中所有重金屬均沒有顯著的人體健康風(fēng)險(xiǎn)。因此，臨沂植煙土壤存在一定程度的Cd和Hg污染風(fēng)險(xiǎn)，應(yīng)采取措施進(jìn)行控制。

煙草；煙區(qū)；農(nóng)田土壤；潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)；人體健康風(fēng)險(xiǎn)評(píng)價(jià)

重金屬因其高毒性、持久性、生物富集性已成為阻礙農(nóng)業(yè)可持續(xù)發(fā)展最突出的問題之一。數(shù)十年來，隨著工業(yè)化和城鎮(zhèn)化的快速發(fā)展，土壤污染程度加重，積聚于土壤中的重金屬通過食物鏈對(duì)生態(tài)環(huán)境和人類健康均造成潛在威脅[1]。

近年來，國內(nèi)外均開展了大量的土壤重金屬污染水平和風(fēng)險(xiǎn)評(píng)價(jià)工作[2]?，F(xiàn)有評(píng)價(jià)方法中，基于土壤背景值的方法有地積累指數(shù)法、富集因子法、污染因子法、污染負(fù)荷指數(shù)法等，基于土壤質(zhì)量標(biāo)準(zhǔn)的方法有污染指數(shù)法等，而潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)則更考慮到了不同重金屬毒性的差異[3]。近來土壤重金屬人體健康風(fēng)險(xiǎn)評(píng)價(jià)的研究也卓然興起，特別是礦區(qū)、工業(yè)區(qū)、農(nóng)業(yè)區(qū)的土壤[2-4]。人體健康風(fēng)險(xiǎn)評(píng)價(jià)，包括致癌風(fēng)險(xiǎn)和非致癌風(fēng)險(xiǎn)，是通過對(duì)不同暴露途徑的評(píng)估進(jìn)行的。人體對(duì)土壤重金屬的3種暴露途徑有浮塵吸入、經(jīng)口攝入或皮膚直接接觸土壤，有專家認(rèn)為在評(píng)價(jià)時(shí)經(jīng)口攝入和接觸是主要的暴露途徑[3]，然而實(shí)際操作中，吸入也經(jīng)常作為一種直接的暴露途徑納入到風(fēng)險(xiǎn)評(píng)價(jià)過程中[4]。近期有學(xué)者提出土壤和農(nóng)產(chǎn)品綜合質(zhì)量指數(shù)法[2]，綜合考慮土壤背景值、土壤限量標(biāo)準(zhǔn)與農(nóng)作物限量標(biāo)準(zhǔn)等參比系列，但對(duì)于尚沒有重金屬限量標(biāo)準(zhǔn)頒布的農(nóng)產(chǎn)品（如煙草）及其產(chǎn)地暫無法使用此方法進(jìn)行評(píng)價(jià)。

煙草（L.）是國內(nèi)外廣泛種植的經(jīng)濟(jì)作物，近年其質(zhì)量安全及產(chǎn)地環(huán)境生態(tài)風(fēng)險(xiǎn)狀況也逐漸引起關(guān)注。因此，本研究調(diào)查采集了山東省臨沂市植煙耕層土壤樣品，檢測(cè)了其中砷（As）、鎘（Cd）、鉻（Cr）、銅（Cu）、汞（Hg）、鎳（Ni）、鉛（Pb）和鋅（Zn）的含量，分析了重金屬在土壤中的富集和污染水平，對(duì)比了各區(qū)域的空間分布差異，并評(píng)估了潛在生態(tài)風(fēng)險(xiǎn)和人體健康風(fēng)險(xiǎn)水平，對(duì)當(dāng)?shù)責(zé)熖锿寥乐亟饘倏刂凭哂兄匾闹笇?dǎo)意義。

1 材料與方法

1.1 研究區(qū)域

臨沂市（117°24′~119°11′E，34°22′~36°13′N），位于山東省東南部，氣候?qū)倥瘻貛Ъ撅L(fēng)大陸性氣候，全年平均氣溫13.6 ℃，年平均降水量797 mm，年平均日照為2288 h[5]，在全國煙草種植區(qū)劃上屬于魯中南低山丘陵烤煙適宜區(qū)[6]。

1.2 樣品采集與檢測(cè)

依據(jù)臨沂市各縣煙草種植面積和煙葉產(chǎn)量的空間分布情況，在市境內(nèi)9個(gè)縣共確定了100個(gè)典型采樣點(diǎn)，其中蒼山、費(fèi)縣、莒南、臨沭、蒙陰、平邑、郯城、沂南和沂水樣點(diǎn)數(shù)分別為6、15、10、6、18、7、4、14和20個(gè)（圖1）。同一田塊耕作層（0~20 cm）隨機(jī)多點(diǎn)采樣，組成混合農(nóng)化土樣，經(jīng)分揀雜物、自然風(fēng)干、研磨后過100目篩備用。土壤樣品加入HNO3、H2O2和HF，微波法消解，用電感耦合等離子質(zhì)譜（ICP-MS）測(cè)定重金屬[7-8]。

圖1 研究區(qū)與采樣點(diǎn)位置示意圖

1.3 評(píng)價(jià)方法

1.3.1 地積累指數(shù) 地積累指數(shù)（I）由MüLLER[9]發(fā)明：

其中，C是元素土壤含量（mg/kg），B是元素土壤背景值（mg/kg），源于中國環(huán)境監(jiān)測(cè)中心主持的國家第一次土壤調(diào)查[10]。按I值，分7個(gè)不同土壤污染程度：＜0，無污染；0~1輕微污染；1~2輕度污染；2~3中度污染；3~4重度污染；4~5極重污染；≥5，極度污染。

1.3.2 污染指數(shù) 單因子污染指數(shù)（P）：

P=C/S(2)

式中，S為重金屬的評(píng)價(jià)標(biāo)準(zhǔn)（mg/kg）。本研究采用煙草產(chǎn)地環(huán)境技術(shù)條件標(biāo)準(zhǔn)NY/T 852—2004[11]和土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)GB 15618—1995標(biāo)準(zhǔn)[12]。內(nèi)梅羅綜合污染指數(shù)（）計(jì)算及評(píng)價(jià)標(biāo)準(zhǔn)參照土壤環(huán)境監(jiān)測(cè)技術(shù)規(guī)范HJ/T 166—2004[13]的方法。

1.3.3 潛在生態(tài)風(fēng)險(xiǎn)指數(shù) 潛在生態(tài)風(fēng)險(xiǎn)公式如下[14]：

其中，EI是元素的單一潛在風(fēng)險(xiǎn)因子；T是毒性反應(yīng)因子，As、Cd、Cr、Cu、Hg、Ni、Pb和Zn分別為10、30、2、5、40、5、5和1。是多元素下整體潛在生態(tài)風(fēng)險(xiǎn)指數(shù)，即所有元素EI的和，根據(jù)和值的大小，生態(tài)風(fēng)險(xiǎn)級(jí)：＜40或＜150，無風(fēng)險(xiǎn)；40~80或150~300輕度風(fēng)險(xiǎn)；80~160或300~600中度風(fēng)險(xiǎn)；160~320或≥600重度風(fēng)險(xiǎn)；≥320極重度風(fēng)險(xiǎn)。

1.3.4 健康風(fēng)險(xiǎn)評(píng)價(jià) 人體健康風(fēng)險(xiǎn)評(píng)價(jià)方法參考美國環(huán)保署相關(guān)方法[15-17]，暴露途徑考慮3種，因此日平均暴露量（）由日均經(jīng)口攝入量（ing）、日均皮膚接觸量（derm）和日均呼吸吸入量（inh）組成，各計(jì)算公式如下，其中主要參數(shù)見表1。

非致癌風(fēng)險(xiǎn)用危害商（）表示，公式如下：

其中參數(shù)及數(shù)值見表2。

危害指數(shù)（）用來表征多種化學(xué)物質(zhì)的非致癌潛在風(fēng)險(xiǎn)，是所有考慮評(píng)價(jià)化學(xué)物質(zhì)的和。若小于1，則基本不會(huì)對(duì)人體產(chǎn)生顯著的健康風(fēng)險(xiǎn)，反之則存在風(fēng)險(xiǎn)。

致癌風(fēng)險(xiǎn)機(jī)率（）則計(jì)算如下：

表1 健康風(fēng)險(xiǎn)評(píng)價(jià)中日平均暴露量計(jì)算相關(guān)參數(shù)

表2 不同暴露途徑重金屬參考濃度（RfD）和致癌斜率（SF）

注：a，數(shù)據(jù)來自《污染場(chǎng)地風(fēng)險(xiǎn)評(píng)估技術(shù)導(dǎo)則（HJ 25.3-2014）》[7]。

Note: a, the values were obtained from the Technical Guidelines for Risk Assessment of Contaminated Sites (HJ 25.3-2014)[7].

其中參數(shù)及數(shù)值也見表2。低于1×10?6被視為無致癌風(fēng)險(xiǎn)，而大于1×10?4則被認(rèn)為有顯著致癌風(fēng)險(xiǎn)，介于兩者之間通常被認(rèn)為是可接受風(fēng)險(xiǎn)范圍。

1.4 數(shù)據(jù)分析

應(yīng)用IBM SPSS 19.0進(jìn)行數(shù)據(jù)分析，方差分析采用Ducan法（＜0.05）；空間分布圖使用Arc GIS 10.0繪制，插值采用克里金法（Kriging）；其他圖表則使用Microsoft Excel 2010或Origin Pro 8.0繪制。

2 結(jié) 果

2.1 土壤重金屬含量

臨沂植煙土壤重金屬含量統(tǒng)計(jì)分析結(jié)果及山東省土壤重金屬背景值列于表3。臨沂植煙土壤As、Cd、Cr、Cu、Hg、Ni、Pb和Zn平均含量分別為3.85、0.09、47.56、14.45、0.07、18.04、17.82和29.21 mg/kg。變異系數(shù)Hg和As較高，說明地區(qū)間分布不均衡。檢驗(yàn)表示8種重金屬均正態(tài)分布，峰度和偏度值也證明這一點(diǎn)。重金屬作為土壤礦物成分，成土母質(zhì)和成土過程對(duì)其含量有著基本的影響，土壤元素背景值能較好反映這種自然因素的影響[20]。從表3可見，土壤As、Cr、Cu、Ni、Pb和Zn平均值均小于山東土壤元素背景值，Cd平均值接近，而Hg平均值則相當(dāng)于背景值的近兩倍，說明除自然因素外，人為因素對(duì)土壤中Hg和Cd的變異影響較大。

2.2 土壤重金屬空間分布

2.2.1 重金屬空間分布和區(qū)域差異分析 圖2和圖3是臨沂土壤重金屬空間分布圖和區(qū)域含量方差比較，分別從空間可視化角度和數(shù)據(jù)角度解析地區(qū)間差異，結(jié)果表明，臨沂北部（蒙陰）、中部（費(fèi)縣）和西南（蒼山）一種或多種重金屬含量較高，東南（郯城、臨沭）重金屬含量普遍相對(duì)較低。

2.2.2 土壤類型差異分析 根據(jù)土壤采樣點(diǎn)全球定位系統(tǒng)（GPS）信息，通過世界土壤數(shù)據(jù)庫（Harmonized World Soil Database，HWSD）[21]查找對(duì)應(yīng)土壤類型信息。發(fā)現(xiàn)按土壤診斷系統(tǒng)分類，臨沂植煙土壤主要有石灰性沖積土（FLc）、飽和始成土（Cme）、弱發(fā)育淋溶土（LVh）、石灰性淋溶土（LVk）、飽和粗骨土（RGe）、潛育淋溶土（LVg）和石灰性粗骨土（RGc）7個(gè)類別。樣品數(shù)量少于5個(gè)的土類不計(jì)入。不同土類間元素比較見表4。土壤Cr、Hg、Pb和Zn不同土壤類別之間無顯著差異，As含量CMe顯著高于RGc，Cd含量CMe和LVh顯著高于RGc，Cu含量FLc顯著高于LVk、LVg和RGc，CMe也顯著高于RGc，Ni含量FLc和CMe顯著高于RGc。因此，臨沂植煙土壤從診斷土類看，石灰性沖積土（FLc）和飽和始成土（CMe）重金屬含量較高，而石灰性粗骨土（RGc）含量較低。

表3 土壤重金屬含量描述性統(tǒng)計(jì)分析

注：a，數(shù)據(jù)源于國家第一次土壤污染調(diào)查[10]；b，*表示<0.05，**表示< 0.01。

Note: a, data origins from the first national soil pollution investigation of China in 1980s; b, * for<0.05 and ** for< 0.01.

圖2 臨沂植煙土壤重金屬空間分布圖

注：相同字母表示同一元素在不同地點(diǎn)無顯著性差異（<0.05）。

Note：The same letter indicates that the same element has no significant difference at<0.05 in different locations.

圖3 臨沂各植煙縣市土壤重金屬含量比較

Fig. 3 Contents of heavy metals in soils from different tobacco-growing counties of Linyi

表4 臨沂不同土壤類型重金屬含量差異

注：每行數(shù)據(jù)中相同字母表示無顯著性差異（<0.05）。

Note: The data in the same row marked with the same letter do not differ significantly at<0.05.

2.3 土壤重金屬富集與污染評(píng)價(jià)

2.3.1 地積累指數(shù) 臨沂植煙土壤重金屬地積累

指數(shù)數(shù)據(jù)見圖4。依據(jù)地積累指數(shù)值，超過93%樣品的As、Cr、Cu、Ni、Pb和Zn處于未污染狀態(tài)，而Cd和Hg未污染樣品比例則分別為81%和73%，其余樣品重金屬大多處于輕微和輕度污染，僅7個(gè)樣品的Hg達(dá)到或超過中度污染，從數(shù)據(jù)分布看屬于點(diǎn)源污染。

2.3.2 污染指數(shù) 從污染指數(shù)（圖4）看，As、Cu和Zn均在超過99%的樣品中達(dá)到清潔水平，其余尚清潔水平，而Pb在全部樣品中達(dá)清潔水平。Cd、Cr、Hg和Ni在超過94%的樣品達(dá)到清潔，分別有1、2、2、3個(gè)樣品達(dá)到輕度污染水平，特別是Hg，有一個(gè)樣品達(dá)到重度污染。依據(jù)內(nèi)梅羅綜合污染指數(shù)，臨沂植煙土壤清潔、尚清潔和輕度污染的樣品比例分別為93%、3%和3%，另有一個(gè)樣品因受Hg單項(xiàng)污染影響而達(dá)到重度污染水平。

2.4 土壤重金屬風(fēng)險(xiǎn)評(píng)價(jià)

2.4.1 潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià) 由圖4可知，所有樣品As、Cr、Cu、Ni、Pb和Zn的生態(tài)風(fēng)險(xiǎn)因子均小于40，說明臨沂植煙土壤這6種元素不存在生態(tài)風(fēng)險(xiǎn)。Cd生態(tài)風(fēng)險(xiǎn)無、輕、中水平的樣品比例分別為76%、19%和5%，Hg相應(yīng)3個(gè)水平樣品比例分別為71%、12%和8%，還有5%和4%樣品的Hg分別達(dá)到重度和極重度生態(tài)風(fēng)險(xiǎn)。通過單元素生態(tài)風(fēng)險(xiǎn)因子計(jì)算整體潛在生態(tài)風(fēng)險(xiǎn)指數(shù)（圖4），說明臨沂植煙土壤有84%樣品重金屬整體無顯著生態(tài)風(fēng)險(xiǎn)，10%樣品處于輕度風(fēng)險(xiǎn)，4%樣品處于中度風(fēng)險(xiǎn)，而有2%樣品處于重度生態(tài)風(fēng)險(xiǎn)。代杰瑞等[22]發(fā)現(xiàn)，山東東部表層土壤重金屬潛在生態(tài)風(fēng)險(xiǎn)Hg＞Cd＞As＞Pb＞Cu＞Cr＞Zn，綜合潛在風(fēng)險(xiǎn)特征也是主要受Hg和Cd元素制約。

2.4.2 健康風(fēng)險(xiǎn)評(píng)價(jià) 從表5看出，臨沂土壤8種重金屬的值和值均低于安全水平，說明土壤重金屬基本沒有非致癌風(fēng)險(xiǎn)。經(jīng)口攝入風(fēng)險(xiǎn)商ing最高，成人和兒童各占的比例為79.7%和86.2%，說明經(jīng)口攝入是土壤重金屬人體健康風(fēng)險(xiǎn)的最主要暴露途徑，這與此前其他區(qū)域研究相一致[23]。

8種土壤重金屬對(duì)成人和兒童的危害指數(shù)順序如下：As＞Cr＞Pb＞Ni＞Cu＞Cd＞Hg＞Zn。兒童比成人風(fēng)險(xiǎn)大可能是因?yàn)閮和男袨楦菀讛z入重金屬，比如吸吮手指等[24]。

同值一樣，土壤重金屬對(duì)兒童的致癌風(fēng)險(xiǎn)也高于成人（表5）。4種致癌重金屬的值順序如下：Cr＞As＞Pb＞Cd，其中Cd和Pb對(duì)成人和兒童的值均低于10?6，說明兩種重金屬不存在顯著的人體健康風(fēng)險(xiǎn)，而As和Cr對(duì)成人和兒童的值介于1×10?4和1×10?6之間，處于可接受的致癌風(fēng)險(xiǎn)范圍內(nèi)。所以，臨沂植煙土壤重金屬不存在顯著的人體健康風(fēng)險(xiǎn)。而且本文中用于分析的是土壤重金屬的全量，土壤中生物有效性含量低于全量，因此人體健康風(fēng)險(xiǎn)評(píng)價(jià)的結(jié)果可能略有高估。

注：圖中上下實(shí)心圓分別表示最大和最小值，中間五角星代表平均值，距形上中下橫線分別表示數(shù)據(jù)的3個(gè)四分?jǐn)?shù)，即75%位數(shù)、中位數(shù)和25%位數(shù)（n=100）。

表5 臨沂植煙土壤重金屬健康風(fēng)險(xiǎn)評(píng)價(jià)

3 討 論

3.1 臨沂植煙土壤重金屬的污染和風(fēng)險(xiǎn)評(píng)價(jià)

對(duì)土壤重金屬污染或風(fēng)險(xiǎn)水平進(jìn)行評(píng)價(jià)的各種參數(shù)，最初設(shè)計(jì)目的不盡相同，分析問題的側(cè)重點(diǎn)不同，選擇的參比含量也有所差異，因此各個(gè)評(píng)價(jià)參數(shù)在某個(gè)層面可能具有其獨(dú)到的優(yōu)勢(shì)，在另外的層面可能存在一定的不足。就本文采用的方法而言，地積累指數(shù)重在反映重金屬分布的自然變化特征及人為影響，但其側(cè)重單一金屬。目前應(yīng)用廣泛的污染指數(shù)法，能明確地反映出各種污染物單一或共同對(duì)土壤環(huán)境的影響，但未考慮土壤中各種污染物毒性的差別，只能反映污染的程度而難以反映污染的質(zhì)變特征[25]。潛在生態(tài)危害指數(shù)法綜合考慮了生物有效性及相對(duì)貢獻(xiàn)與地理空間差異等特點(diǎn)，但其加權(quán)帶有一定的主觀性。健康風(fēng)險(xiǎn)評(píng)價(jià)主要估算污染物與負(fù)面健康效應(yīng)之間的關(guān)系，使公眾更易于理解土壤污染狀況對(duì)生物的影響，但對(duì)土壤本身污染程度的反映不及其他評(píng)價(jià)方法，而且根據(jù)不同區(qū)域的不同情況確定評(píng)價(jià)參數(shù)值尤為重要。各個(gè)參數(shù)基于的參比系統(tǒng)是其評(píng)價(jià)側(cè)重點(diǎn)的基礎(chǔ)，典型的參比是土壤背景值或土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)。基于土壤背景值而得到的評(píng)價(jià)結(jié)果在表現(xiàn)環(huán)境和生態(tài)意義上略顯不足[26]，而單一采用土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)則過分強(qiáng)調(diào)其統(tǒng)一性，不能適應(yīng)我國土壤多樣化的特點(diǎn)。所以，在評(píng)價(jià)過程中針對(duì)性地選擇不同的參數(shù)，可以從多個(gè)角度對(duì)土壤重金屬水平進(jìn)行評(píng)判。本文通過地積累指數(shù)可發(fā)現(xiàn)Cd和Hg相較其他元素在臨沂植煙土壤中略有富集，通過污染因子分析發(fā)現(xiàn)土壤樣品重金屬基本處于尚清潔水平，而在圖4中更突顯了Cd和Hg的潛在生態(tài)風(fēng)險(xiǎn)及對(duì)整體風(fēng)險(xiǎn)的影響，不過臨沂植煙土壤重金屬并不存在顯著的人體健康風(fēng)險(xiǎn)。

3.2 臨沂植煙土壤重金屬的空間分布規(guī)律

土壤類型和土壤利用方式是影響土壤重金屬含量的主要因素，究其原因，前者是土壤重金屬的自然來源——成土母巖母質(zhì)和成土過程的反映，后者則是土壤重金屬差異產(chǎn)生的人為因素。不同利用方式下對(duì)土壤環(huán)境的影響存在顯著差異[27]。本文中，臨沂作為山東傳統(tǒng)煙區(qū)，數(shù)十年來，煙區(qū)肥料、農(nóng)藥等農(nóng)用物資由當(dāng)?shù)叵嚓P(guān)主管部門進(jìn)行統(tǒng)一采購，施用方式和用量也根據(jù)生產(chǎn)規(guī)范進(jìn)行，所以臨沂植煙土壤人為影響在區(qū)域內(nèi)造成的變異要遠(yuǎn)小于其他農(nóng)業(yè)利用類型，因此使得自然因素的影響更顯得突出。自然因素中，地形主要是通過引起物質(zhì)、能量的再分配而間接地作用于土壤的形成過程[28]。臨沂地勢(shì)西北高東南低，沿沂河、沭河往北的沂沭斷裂帶在地勢(shì)上將臨沂分成西北、東南兩部分，因此在圖2重金屬空間分布中也能明顯看出一條縱貫?zāi)媳钡姆纸缇€的存在，分界線以東的莒南、臨沭和位于南部沖積平原的郯城土壤重金屬含量顯著較低，特別是Cr、Hg、Ni、Pb，在空間分布上可能受地質(zhì)自然因素的影響相對(duì)較強(qiáng)。沂山、蒙山、尼山3個(gè)山脈在臨沂境內(nèi)西部由北而南排列，中間有河道分割，也對(duì)土壤重金屬北部、中部和西南3個(gè)重點(diǎn)富集區(qū)域的形成起到了決定性作用。

3.3 臨沂植煙土壤重金屬的來源和控制重點(diǎn)

區(qū)別于工業(yè)和城市土壤重金屬來源的復(fù)雜性和礦區(qū)土壤的單一性[3,29-30]，農(nóng)田土壤重金屬主要來源為母質(zhì)、大氣沉降、肥料和灌溉水等。本項(xiàng)目組前期[8]通過主成分分析與聚類分析發(fā)現(xiàn)山東植煙土壤Cr和Ni含量變異主要受成土母質(zhì)控制，As、Cd、Cu、Pb和Zn的變異主要受人為工農(nóng)業(yè)活動(dòng)控制，而Hg的變異主要受大氣沉降控制；土壤重金屬外源貢獻(xiàn)的計(jì)算還發(fā)現(xiàn)，山東植煙土壤重金屬主要外源是大氣沉降，其次是肥料，而控制的重點(diǎn)是Cd和Hg，與本研究通過生態(tài)和健康風(fēng)險(xiǎn)分析得到的結(jié)論是一致的。

另外，本項(xiàng)目組對(duì)臨沂煙葉重金屬檢測(cè)分析表明，煙葉對(duì)Cd的平均富集系數(shù)為8.8，而對(duì)其他重金屬的平均富集系數(shù)均小于1，說明煙草是Cd的單一富集植物且在葉部富集，這與項(xiàng)目組室內(nèi)培養(yǎng)試驗(yàn)的結(jié)論也是一致的[31-32]。正因如此，Cd是臨沂煙草重金屬研究和控制的重中之重。

4 結(jié) 論

對(duì)山東省臨沂市植煙土壤重金屬As、Cd、Cr、Cu、Hg、Ni、Pb和Zn的污染與風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果表明，臨沂植煙土壤93%樣品重金屬整體處于清潔水平，僅Cd和Hg存在輕度富集或污染；84%土壤樣品重金屬整體無顯著生態(tài)風(fēng)險(xiǎn)，10%樣品處于輕度風(fēng)險(xiǎn)；土壤中8種重金屬生態(tài)風(fēng)險(xiǎn)大小順序?yàn)镠g>Cd>As> Ni>Pb>Cu>Cr>Zn，且對(duì)成人和兒童均無顯著的健康風(fēng)險(xiǎn)；Cd和Hg是植煙土壤今后重點(diǎn)控制的重金屬元素，特別是臨沂北部、西南和中部部分地區(qū)。

[1] LUO L, MA Y, ZHANG S, et al. An inventory of trace element inputs to agricultural soils in China[J]. Journal of Environmental Management, 2009, 90(8): 2524-2530.

[2] 王玉軍，吳同亮，周東美，等. 農(nóng)田土壤重金屬污染評(píng)價(jià)研究進(jìn)展[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2017，36（12）：2365-2378.

WANG Y J, WU T L, ZHOU D M, et al. Advances in soil heavy metal pollution evaluation based on bibliometrics analysis[J]. Journal of Agro-Environment Science, 2017, 36(12): 2365-2378.

[3] LI Z, MA Z, VAN DER KUIJP T J, et al. A review of soil heavy metal pollution from mines in China: Pollution and health risk assessment[J]. Science of the Total Environment, 2014: 468-469, 843-853.

[4] HU W, HUANG B, HE Y, et al. Assessment of potential health risk of heavy metals in soils from a rapidly developing region of China[J]. Human & Ecological Risk Assessment An International Journal, 2016, 22(1): 211-225.

[5] 張磊，趙海軍，潘婕，等. 臨沂市旅游氣候資源分析與評(píng)價(jià)[C]//中國氣象學(xué)會(huì). 第32屆中國氣象學(xué)會(huì)年會(huì)，2015：252-256.

ZHANG L, ZHAO H J, PAN J, et al. Analysis and Evaluation of Tourism Climatic Resources in Linyi[C]// China Meteorological Society. The proceedings of 32nd annual meeting of China Meteorological Society, 2015: 252-256.

[6] 王彥亭，謝劍平，李志宏. 中國煙草種植區(qū)劃[M]. 北京：科學(xué)出版社，2010：61.

WANG Y T, XIE J P, LI Z H. Regionalization of tobacco planting in China[M]. Beijing: Science Press, 2010: 61.

[7] 環(huán)境保護(hù)部科技標(biāo)準(zhǔn)司. 污染場(chǎng)地風(fēng)險(xiǎn)評(píng)估技術(shù)導(dǎo)則：HJ 25.3—2014[S]. 北京：中國環(huán)境出版社，2004: 16-17.

Department of Science, Technology and Standards, Ministry of Environmental Protection, PRC. Technical guidelines for risk assessment of contaminated sites: HJ 25.3-2014[S]. Beijing: China Environment Press, 2004: 16-17.

[8] LIU H, ZHANG Y, ZHOU X, et al. Source identification and spatial distribution of heavy metals in tobacco growing soils in Shandong province of China with multivariate and geostatistical analysis[J]. Environmental Science and Pollution Research, 2017, 24: 5964-5975.

[9] MüLLER G. Schwermetalle in den Sedimenten des Rheins-Ver?nderungen seit 1971[J]. Umschau in Wissenschaft und Technik, 1979, 79 (24): 778-783 (in German).

[10] 中國環(huán)境監(jiān)測(cè)總站. 中國土壤元素背景值[M]. 北京：中國環(huán)境科學(xué)出版社，1990：329-380.

China National Environmental Monitoring Centre (CNEMC). The background values of elements in soils of China[M]. Beijing: China Environmental Science Press, 1990: 329-380.

[11] 農(nóng)業(yè)部. 煙草產(chǎn)地環(huán)境技術(shù)條件：NY/T 852—2004 [S]. 北京：中國標(biāo)準(zhǔn)出版社，2004：2.

Ministry of Agriculture, PRC. Environmental requirement for growing area of tobacco: NY/T 852-2004[S]. Beijing: China Standard Press, 2004: 2.

[12] 國家環(huán)境保護(hù)總局. 土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)GB 15618—1995[S]. 北京：中國標(biāo)準(zhǔn)出版社，1995：2.

State Environmental Protection Administration of China. Environmental quality standard for soils: GB 15618—1995[S]. Beijing: China Environment Press, 1995: 2.

[13] 國家環(huán)境保護(hù)總局. 土壤環(huán)境監(jiān)測(cè)技術(shù)規(guī)范HJ/T 166—2004 [S]. 北京：中國環(huán)境出版社，2004：16-17.

State Environmental Protection Administration of China. The technical specification for soil environmental monitoring: HJ/T 166-2004[S]. Beijing: China Environment Press, 2004: 16-17.

[14] HAKANSON L. An ecological risk index for aquatic pollution control: a sedimentological approach[J]. Water Research, 1980, 14(8): 975-1001.

[15] US EPA. Risk assessment guidance for superfund. Volume I. Human health evaluation manual (Part A)[M]. Washington, DC: Office of Solid Waste and Emergency Response, 1989.

[16] US EPA. Exposure factors handbook (Final report, 1997)[M]. Washington, DC: Office of Solid Waste and Emergency Response, 1997.

[17] US EPA. Baseline human health risk assessment[M]. Washington, DC: Office of Solid Waste and Emergency Response, 2001.

[18] US EPA. Supplemental guidance for developing soil screening levels for Superfund sites[M]. Washington, DC: Office of Solid Waste and Emergency Response, 2002.

[19] ZHOU X, CAO Z, MA Y, et al. Concentrations, correlations and chemical species of PM2.5/PM10 based on published data in China: Potential implications for the revised particulate standard[J]. Chemosphere, 2016, 144: 518-526.

[20] 陳懷滿. 土壤—植物系統(tǒng)中的重金屬污染[M]. 北京：科學(xué)出版社，1996.

CHEN H M. Heavy metals pollution in soil-plant system[M]. Beijing: Science Press, 1996.

[21] FAO/IIASA/ISRIC/ISSCAS/JRC. Harmonized World Soil Database (version 1.2)[DB]. FAO, Rome, Italy and IIASA, Laxenburg, Austria, 2012.

[22] 代杰瑞，胡雪平，王增輝. 土壤重金屬污染及潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)——以山東半島藍(lán)色經(jīng)濟(jì)區(qū)為例[J]. 安全與環(huán)境學(xué)報(bào)，2013，13（1）：134-137.

DAI J R, HU X P, WANG Z H. Characteristics of heavy metals pollution in soil and potential ecological risk assessment—a case study of Blue Economic Zone of Shandong Peninsula[J]. Journal of Safety and Environment, 2013, 13(1): 134-137.

[23] XIAO Q, ZONG Y, LU S. Assessment of heavy metal pollution and human health risk in urban soils of steel industrial city (Anshan), Liaoning, Northeast China[J]. Ecotoxicology and Environmental Safety, 2015, 120: 377-385.

[24] ZHAO L, XU Y, HOU H, et al. Source identification and health risk assessment of metals in urban soils around the Tanggu chemical industrial district, Tianjin, China[J]. Science of the Total Environment, 2014: 468-469, 654-662.

[25] 范拴喜，甘卓亭，李美娟，等. 土壤重金屬污染評(píng)價(jià)方法進(jìn)展[J]. 中國農(nóng)學(xué)通報(bào)，2010，26（17）：310-315.

FAN S X, GAN Z T, LI M J, et al. Progress of Assessment Methods of Heavy Metal Pollution in Soil[J]. Chinese Agricultural Science Bulletin, 2010, 26(17): 310-315.

[26] 劉鳳枝，師榮光，徐亞平，等. 耕地土壤重金屬污染評(píng)價(jià)技術(shù)研究——以土壤中鉛和鎘污染為例[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2006，25（2）：422-426.

LIU F Z, SHI R G, XU Y P, et al. The study of assessment technology for farm land soil heavy metal pollutions[J]. Journal of Agro-Environment Science, 2006, 25(2): 422-426.

[27] HU Y, CHENG H. Application of stochastic models in identification and apportionment of heavy metal pollution sources in the surface soils of a large-scale region[J]. Environmental Science & Technology, 2013, 47(8): 3752-3760.

[28] 張鳳榮. 土壤地理學(xué)[M]. 北京：中國農(nóng)業(yè)出版社，2002.

ZHANG F R. Soil geography[M]. Beijing: China Agriculture Press, 2002.

[29] 戴彬，呂建樹，戰(zhàn)金成，等. 山東省典型工業(yè)城市土壤重金屬來源、空間分布及潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)[J]. 環(huán)境科學(xué)，2015，36（2）：507-515.

DAI B, LYU J S, ZHAN J C, et al. Assessment of sources，spatial distribution and ecological risk of heavy metals in soils in a typical industry-based city of Shandong Province，Eastern China[J]. Environmental Science, 2015, 36(2): 507- 515.

[30] 范明毅，楊皓，黃先飛，等. 典型山區(qū)燃煤型電廠周邊土壤重金屬形態(tài)特征及污染評(píng)價(jià)[J]. 中國環(huán)境科學(xué)，2016，36（8）：2425-2436.

FAN M Y, YANG H, HUANG X F, et al. Chemical forms and risk assessment of heavy metals in soils around a typical coal-fired power plant located in the mountainous area[J]. China Environmental Science, 2016, 36(8): 2425-2436.

[31] LIU H, WANG H, MA Y, et al. Role of transpiration and metabolism in translocation and accumulation of cadmium in tobacco plants (L.)[J]. Chemosphere, 2016, 144: 1960-1965.

[32] 劉海偉，石屹，梁洪波. 煙草和卷煙中重金屬遷移分配的研究進(jìn)展[J]. 中國農(nóng)業(yè)科技導(dǎo)報(bào)，2013，15（2）：153-158.

LIU H W, SHI Y, LIANG H B. Research progress on transfer and distribution of heavy metal in tobacco and cigarette[J]. Journal of Agricultural Science and Technology, 2013, 15 (2): 153-158.

Spatial Distribution, Potential Ecological and Health Risk Assessment of Heavy Metals in Surface Soils from Tobacco Cultivating Areas of Linyi, Shandong Province, China

LIU Haiwei1, ZONG Hao2, WANG Haiyun1, LIU Wentao2, YANG Jutian2*, SHI Yi1, ZHANG Pingde2, LIU Li2

(1. Tobacco Research Institute of Chinese Academy of Agricultural Sciences, Key Laboratory of Tobacco Biology and Processing, Ministry of Agriculture, Qingdao 266101, China; 2. Shandong Linyi Tobacco Co., Ltd., Linyi, Shandong 276000, China)

The spatial distribution and risk assessment of heavy metals are the foundation of soil remediation. The heavy metals (As, Cd, Cr, Cu, Hg, Ni, Pb, and Zn) in the surface soils of tobacco (L.) fields in Linyi City of Shandong Province, a traditional tobacco growing area, were collected and determined, and the difference between counties and soil types was analyzed, with spatial distribution maps generated by the geographic information system. The enrichment of heavy metals in soil was assessed by index of geo-accumulation and pollution index, and the risk assessment included potentially ecological risk and health risk assessment. The results showed that, (1) the mean heavy metal concentrations in Linyi tobacco growing soils were 3.85, 0.09, 47.56, 14.45, 0.07, 18.04, 17.82 and 29.21 mg/kg, respectively; (2) the concentrations of heavy metals in soil samples from north, middle and southwest of Linyi were relatively higher than those from other areas; (3) some soil samples were enriched with Cd and/or Hg, which were the only two metals posing a potentially high ecological risk to the local environment, and the order of ecological risk for the heavy metals was Hg>Cd>As>Ni> Pb>Cu>Cr>Zn; (4) The health risk assessment showed no considerable non-carcinogenic or carcinogenic risks for children and adults from these elements.

tobacco; tobacco growing area; agricultural soil; potential ecological risk assessment; health risk assessment

S572.01

1007-5119（2018）04-0041-09

10.13496/j.issn.1007-5119.2018.04.006

公益性行業(yè)（農(nóng)業(yè)）科研專項(xiàng)“煙草增香減害關(guān)鍵技術(shù)研究與示范”（201203091）；中國煙草總公司山東省公司科技計(jì)劃項(xiàng)目“山東濃香低害煙葉研究與開發(fā)”（201101）、“臨沂煙葉脈間失綠成因分析及綜合解決措施研究”（201620）；中國農(nóng)業(yè)科學(xué)院科技創(chuàng)新工程“煙草栽培與調(diào)制創(chuàng)新團(tuán)隊(duì)”（ASTIP-TRIC03）

劉海偉（1981-），男，博士，副研究員，主要研究方向?yàn)檗r(nóng)業(yè)環(huán)境污染治理。E-mail：heaveyliu@163.com。

，E-mail：lyycyjt@163.com

2018-01-24

2018-06-15