郭軍康，趙 瑾，魏 婷，李艷萍，呂 欣，任 倩，王小紅

（陜西科技大學(xué)環(huán)境科學(xué)與工程學(xué)院，西安 710021）

設(shè)施栽培是我國蔬菜生產(chǎn)的主要方式，為經(jīng)濟(jì)發(fā)展和人民群眾生活水平改善做出了巨大的貢獻(xiàn)，自推進(jìn)“菜籃子”工程建設(shè)以來，設(shè)施菜地面積不斷擴(kuò)大。目前我國蔬菜產(chǎn)業(yè)已形成較大規(guī)模，其產(chǎn)值已成為繼糧食產(chǎn)業(yè)之后農(nóng)業(yè)第二大產(chǎn)業(yè)[1]。2015年，西安市設(shè)施菜地種植面積達(dá)2.8萬hm2，占蔬菜總種植面積40%以上，預(yù)計到2020年，設(shè)施菜地面積將達(dá)到6.67萬hm2，成為蔬菜產(chǎn)業(yè)的主體[2]。但隨著設(shè)施菜地種植年限的延長，不當(dāng)?shù)墓芾泶胧┊a(chǎn)生的負(fù)面影響日益嚴(yán)重，土壤氮磷鉀含量明顯增加，pH降低[3-4]，重金屬累積[5]等問題導(dǎo)致土壤質(zhì)量下降。其中土壤重金屬含量直接關(guān)系到食物安全及人類和動物健康[6]，然而重金屬的生態(tài)危害性不僅與其總量有關(guān)，更大程度上取決于重金屬的形態(tài)特征[7]。重金屬的不同形態(tài)直接影響到重金屬的毒性及其在自然界中的遷移轉(zhuǎn)化過程。

西安市位于關(guān)中平原核心地帶，是我國西部唯一“國家中心城市”。隨著西安城市化進(jìn)程加快，城市郊區(qū)設(shè)施菜地土壤受到污染的風(fēng)險增加[8-10]。研究西部大型城市周邊不同年限設(shè)施菜地的重金屬分布特征、土壤理化性質(zhì)及重金屬形態(tài)變化與污染風(fēng)險評價，能夠?qū)ξ鞑可鷳B(tài)脆弱區(qū)設(shè)施農(nóng)業(yè)土壤質(zhì)量變化和安全利用提供重要的科學(xué)支撐。

為探討西安市郊不同種植年限設(shè)施菜地土壤重金屬的形態(tài)變化及污染程度，本研究以西安市郊不同年限設(shè)施菜地土壤為研究對象，研究設(shè)施菜地理化性質(zhì)，探討了主要重金屬污染元素（Cd、Pb）在設(shè)施菜地中的含量和賦存形態(tài)，采用地積累污染指數(shù)法和潛在生態(tài)風(fēng)險指數(shù)法評價研究區(qū)重金屬污染程度。實驗對研究重金屬環(huán)境生物效應(yīng)及其遷移轉(zhuǎn)化規(guī)律具有重要意義，同時可為設(shè)施菜地的可持續(xù)發(fā)展提供數(shù)據(jù)支持和理論參考。

1 材料與方法

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

西安市地處北緯33°39′～34°45′、東經(jīng)107°40′～109°49′，屬暖溫帶半濕潤易旱季風(fēng)氣候，年平均氣溫為13.6℃，年降水總量在595.9～732.9 mm。本次采集的設(shè)施菜地土壤來自西安市設(shè)施菜地的主要分布地區(qū)（長安區(qū)和高陵區(qū)），設(shè)施類型均為塑料拱棚，土壤為黃土母質(zhì)，土壤類型為塿土。通過實地調(diào)查，研究區(qū)設(shè)施菜地生產(chǎn)類型與管理方式相同，主要生產(chǎn)青菜等，一年種植4～5茬，采用輪作制，肥料主要為有機(jī)肥和復(fù)合肥（磷酸二氫鉀和尿素等），按照180～270 kg·hm-2一茬一施。灌溉水為地下用水。為保證耕作管理措施和土壤本底值的一致性，選取相鄰菜地，并用GPS定位采樣地點和海拔高度（田間調(diào)查問卷數(shù)據(jù)，未發(fā)表）。對選取的設(shè)施菜地進(jìn)行長條形折線法采樣。取0～20 cm層土樣，每棚采集3點，組成混合樣，四分法留樣裝入塑料袋中標(biāo)記密封，總樣本數(shù)為60。樣本帶回實驗室風(fēng)干、研磨、過篩后備用。

1.2 實驗方法

1.2.1 土壤理化性質(zhì)的測定

土壤理化性質(zhì)的測定參考《土壤農(nóng)業(yè)化學(xué)分析方法》[11]：pH值測定采用電位法，水土比為2.5∶1；水解性氮測定采用堿解擴(kuò)散法；有機(jī)質(zhì)測定采用重鉻酸鉀滴定法；堿解氮測定采用堿解擴(kuò)散法；速效磷測定采用0.5 mol·L-1碳酸氫鈉浸提-鉬銻抗比色法；速效鉀測定采用1 mol·L-1醋酸銨浸提-火焰光度比色法。

1.2.2 重金屬總量測定和形態(tài)分析

重金屬總量測定：將土壤樣品研磨，過篩，烘干。稱取0.25 g待測樣品于消化管內(nèi)，加入10 mL硝酸和4 mL氫氟酸，放入電熱消解儀（Digi Block ED54，Lab Tech）進(jìn)行消解，趕酸，直至約1 mL左右，冷卻，用1%HNO3定容至25 mL。按照相同的程序做空白試驗[12]。

重金屬形態(tài)分析：采用歐洲標(biāo)準(zhǔn)測試分析委員會提出的四步分級提取法（簡稱BCR法）[13]。提取形態(tài)依次為酸提取態(tài)（F1）、可還原態(tài)（F2）、可氧化態(tài)（F3）、殘渣態(tài)（F4）。

所有待測的樣品重金屬含量均采用原子吸收光譜儀（AAS；ZEEnit700P/650P，Analytik-Jena，Germa?ny）測定，所有樣品測定3次。為了進(jìn)行質(zhì)量控制，同步分析了國家標(biāo)準(zhǔn)物質(zhì)陜西洛川黃土GBW 07454（GSS-25），結(jié)果顯示各重金屬元素含量相對標(biāo)準(zhǔn)偏差均小于5%。

1.3 土壤環(huán)境質(zhì)量評價方法

參照《溫室蔬菜產(chǎn)地環(huán)境質(zhì)量評價標(biāo)準(zhǔn)》（HJ 333—2006），采用地積累指數(shù)法和潛在生態(tài)風(fēng)險指數(shù)法對西安市郊設(shè)施菜地土壤Cd和Pb污染程度展開評價。

地積累指數(shù)法[14]（Mull指數(shù)）：

式中：Cn是元素n在沉積物中的含量；Bn是沉積物中該元素的地球化學(xué)背景值；k為考慮各地巖石差異可能會引起背景值的變動而取的系數(shù)，k=1.5。

潛在生態(tài)風(fēng)險指數(shù)法（Hakanson）[15]：

本文采用陜西省西安市土壤重金屬元素背景值為參比值[16]。和RI分別為單項重金屬毒性響應(yīng)參數(shù)、潛在生態(tài)風(fēng)險系數(shù)和綜合潛在生態(tài)風(fēng)險指數(shù)。重金屬Pb、Cd的毒性響應(yīng)參數(shù)分別為5和30[14-15]。重金屬污染系數(shù)和生態(tài)危害指數(shù)分級標(biāo)準(zhǔn)見表1。

表1 重金屬生態(tài)風(fēng)險等級劃分標(biāo)準(zhǔn)Table1 Criteria for classification of ecological risk of heavy metals

1.4 數(shù)據(jù)處理

利用Excel 2016和Origin 9.0進(jìn)行數(shù)據(jù)分析和作圖，利用SPSS 21.0進(jìn)行統(tǒng)計分析，一維方差分析法進(jìn)行顯著性檢驗，P＜0.05為差異顯著。

表2 不同種植年限設(shè)施菜地的理化指標(biāo)Table2 Physiochemical indexes of greenhouse vegetable soil with different planting years

2 結(jié)果與討論

2.1 不同種植年限設(shè)施菜地土壤理化性質(zhì)的變化

連年高負(fù)荷耕種使設(shè)施菜地土壤理化性質(zhì)發(fā)生變化（表2）。除干物質(zhì)外，其他理化性質(zhì)呈現(xiàn)顯著性變化，整體上有效鉀、氮、磷含量隨種植年限的增加而增加，在21～25 a內(nèi)達(dá)到最大，比初建棚土壤分別增加了201.1%、131.2%、190.5%。有機(jī)質(zhì)隨著年限的延長呈現(xiàn)先增長后降低的趨勢。這可能與城市化快速發(fā)展和居民對蔬菜消費量日益增加導(dǎo)致設(shè)施菜地土壤耕作頻繁有關(guān)，加上關(guān)中地帶設(shè)施菜地有機(jī)肥的投入量少，偏施氮肥，營養(yǎng)分配不均勻[17]，多茬生長的蔬菜種植過程中土壤有機(jī)質(zhì)得不到有效補(bǔ)充，造成了年限較高的設(shè)施菜地有機(jī)質(zhì)含量下降、水解性氮偏高[18]。這與李玉娣等[19]探討陜西省涇陽縣蔬菜大棚種植年限對土壤養(yǎng)分影響的結(jié)果相一致。而研究區(qū)不同年限設(shè)施菜地土壤pH值在7.58～8.07之間，均低于初建棚時的8.27，表現(xiàn)為隨著種植年限的增加而緩慢下降，呈酸化的趨勢。這與曾希柏等[20]探討山東壽光不同種植年限設(shè)施土壤的酸化與鹽漬化的結(jié)果相似。土壤氮磷鉀以及有機(jī)質(zhì)隨著種植年限的延長而增加，不合理施用設(shè)施菜地農(nóng)藥肥料導(dǎo)致pH降低[21]。

2.2 不同種植年限設(shè)施菜地土壤重金屬含量變化

近年來，設(shè)施菜地土壤重金屬超標(biāo)問題日益受到研究者關(guān)注[19-21]。通過與《溫室蔬菜產(chǎn)地環(huán)境質(zhì)量評價標(biāo)準(zhǔn)》（HJ 333—2006）中重金屬含量對比，西安市不同年限設(shè)施菜地重金屬Cd和Pb總量超標(biāo)率分別為60%和20%（表3）。從不同年限設(shè)施菜地重金屬含量（以5 a為時段）的變化幅度來看，Cd的總量隨著年限延長呈增加趨勢，分別增加了21.4%（6～10 a）、79.4%（11～15 a）、139.3%（16～20 a）、3.4%（21～25 a），在16～20 a內(nèi)增長幅度最高，且均超過其陜西省土壤背景值，11 a后每5年時間段的土壤Cd分別是標(biāo)準(zhǔn)限值的1.53、3.65、3.78倍。同時土壤Pb的總量隨著年限增加而緩慢增加，均超過背景值，除11～15 a，其他年限段分別增加了32.6%、1.9%、61.8%。在11～15 a，Pb總量降低，這可能與土壤環(huán)境有關(guān)。蘆小軍等[22]發(fā)現(xiàn)土壤微生物會促進(jìn)印度芥菜對土壤Pb的吸收，導(dǎo)致土壤Pb含量的下降。

不同年限設(shè)施菜地土壤重金屬元素BCR形態(tài)分析結(jié)果表明，隨著使用年限的增加，土壤Cd、Pb有效態(tài)含量呈遞增趨勢，除6～10 a外，其他年限段Cd有效態(tài)含量分別增加了28.6%、27.8%、117.4%，均超過背景值，在21～25 a內(nèi)增長幅度最高且超過標(biāo)準(zhǔn)限值。而Pb有效態(tài)含量均未超出背景值，在6～10 a和11～15 a內(nèi)分別增加了15.8%和5.2%，使用20 a后有效態(tài)Pb含量降低。這可能與植物吸收有效態(tài)Pb有關(guān)[22]。

與初建設(shè)施菜地土壤相比，土壤重金屬Cd、Pb總量和有效態(tài)隨種植年限延長而有所增加，李一蒙等[23]對河南開封土壤重金屬調(diào)查發(fā)現(xiàn)Cd超標(biāo)問題最為突出。井永蘋等[24]研究山東壽光不同種植年限對設(shè)施菜地土壤重金屬累積量的影響，結(jié)果顯示土壤Cd、Pb含量隨種植年限增加而升高。鄧煥廣等[25]研究表明隨著年限延長酸提取態(tài)As含量呈持續(xù)增加的趨勢，除殘渣態(tài)外，其他形態(tài)As與土壤理化性質(zhì)均存在相關(guān)性。從已有的研究結(jié)果來看，我國各地的設(shè)施土壤重金屬均呈現(xiàn)明顯的累積趨勢，與本研究所得的結(jié)果一致。

表3 不同年限設(shè)施菜地重金屬含量變化（mg·kg-1）Table3 changes of heavy metal contents in greenhouse vegetable fields in different years（mg·kg-1）

2.3 不同種植年限設(shè)施菜地土壤重金屬的形態(tài)分布

土壤中重金屬元素遷移能力的強(qiáng)弱與它們的化學(xué)形態(tài)密切相關(guān)。不同年限的土壤重金屬形態(tài)分布不同（圖1和圖2）；在0～15 a內(nèi)，Cd的形態(tài)特征表現(xiàn)為可還原態(tài)（45.08%）＞酸提取態(tài)（23.76%）＞殘渣態(tài)（18.23%）＞可氧化態(tài)（12.20%）＞水溶態(tài)（0.72%）；而在使用15 a后，Cd的形態(tài)特征表現(xiàn)為殘渣態(tài)（46.40%）＞可還原態(tài)（33.06%）＞酸提取態(tài)（14.12%）＞可氧化態(tài)（6.01%）＞水溶態(tài)（0.41%）。值得注意的是殘渣態(tài)Cd隨著年限增加呈增長趨勢，在16～20 a，增加幅度最高；通過田間調(diào)查發(fā)現(xiàn)，酸提取態(tài)含量與研究區(qū)農(nóng)戶管理模式有關(guān)，即每隔5年添加未被污染的大田土壤從而改善設(shè)施菜地連年耕作造成土壤質(zhì)量下降的現(xiàn)象（田間調(diào)查問卷數(shù)據(jù)，未發(fā)表）?？蛇€原態(tài)和酸提取態(tài)是設(shè)施菜地土壤中Cd的主要存在形態(tài)，酸提取態(tài)與可還原態(tài)之和大于可氧化態(tài)與殘渣態(tài)之和，說明不同年限設(shè)施菜地中Cd的生物有效性大，遷移的風(fēng)險更大。在年限0～15 a內(nèi)，Pb的形態(tài)特征表現(xiàn)為殘渣態(tài)（38.16%）＞酸提取態(tài)（28.42%）＞可還原態(tài)（19.48%）＞可氧化態(tài)（8.80%）＞水溶態(tài)（5.15%）；而使用15 a后，Pb的形態(tài)特征表現(xiàn)為殘渣態(tài)（41.21%）＞可還原態(tài)（29.90%）＞酸提取態(tài)（14.28%）＞可氧化態(tài)（13.11%）＞水溶態(tài)（1.50%）。土壤Pb由酸提取態(tài)向可還原態(tài)轉(zhuǎn)化，說明不同年限設(shè)施菜地土壤理化性質(zhì)的變化會導(dǎo)致Pb生物遷移能力減弱。殘渣態(tài)Pb占總量的百分?jǐn)?shù)比其他形態(tài)高，其范圍為35.42%～42.20%，說明不同年限設(shè)施菜地中Pb的生物遷移能力低，殘渣態(tài)Pb是設(shè)施菜地土壤中Pb的主要存在形態(tài)。隨著年限延長，可還原態(tài)Pb逐漸增加，這可能與土壤中Fe和Mn的氫氧化物特別是Mn的氫氧化物對Pb2+有很強(qiáng)的專性吸附能力有關(guān)[26]。

2.4 不同年限設(shè)施菜地土壤污染評價

通過分析設(shè)施菜地土壤重金屬總量不能全面反映其潛在生態(tài)危害性，因此在地積累污染評價和潛在生態(tài)風(fēng)險評價中采用有效態(tài)含量[27]。根據(jù)1.3中公式計算兩種重金屬的地積累污染指數(shù)和潛在生態(tài)風(fēng)險系數(shù)（表4）。由表4可知，不同年限設(shè)施菜地土壤重金屬Cd的地積累污染指數(shù)依次為21～25 a＞16～20 a＞11～15 a＞1～5 a＞6～10 a，在 11～20 a內(nèi)，0＜Igeo＜1，在21～25 a內(nèi)Igeo＞1，隨著年限增長設(shè)施菜地土壤中Cd的污染程度由無污染到中度污染；而不同種植年限設(shè)施菜地土壤Pb的Igeo均小于0，說明研究區(qū)各年限設(shè)施菜地Pb無污染。根據(jù)重金屬Eri、RI值與風(fēng)險程度的關(guān)系（表1），不同年限設(shè)施菜地Cd的Eri介于33.60～125.07之間，平均值為59.08，其中11～15、16～20、21～25 a年限段有中等生態(tài)風(fēng)險，其他為低生態(tài)風(fēng)險；各年限設(shè)施菜地Pb的Eri均小于40，屬于低生態(tài)風(fēng)險。各年限設(shè)施菜地的RI值介于36.48～128.03之間，平均值為61.95，除21～25 a設(shè)施菜地為中等生態(tài)風(fēng)險，其他都為低生態(tài)風(fēng)險，因此設(shè)施菜地Cd是不同年限設(shè)施菜地的主要風(fēng)險來源，其生態(tài)風(fēng)險系數(shù)對潛在風(fēng)險指數(shù)的貢獻(xiàn)高達(dá)94.36%，重金屬Pb對潛在生態(tài)風(fēng)險指數(shù)的平均貢獻(xiàn)僅占5.64%。即設(shè)施菜地土壤中重金屬的生態(tài)風(fēng)險Cd＞Pb。這與李雪芳[10]研究的西安市市郊菜地重金屬污染趨勢一致。由此可知：不同年限設(shè)施菜地土壤中重金屬Cd的污染較為嚴(yán)重，尤其是高年限設(shè)施菜地，應(yīng)引起有關(guān)部門的高度重視。

圖1 不同年限設(shè)施菜地重金屬Cd形態(tài)分布Figure1 Distribution of Cd speciation in greenhouse vegetable soil with different years

圖2 不同年限設(shè)施菜地重金屬Pb形態(tài)分布Figure2 Distribution of Pb speciation in greenhouse vegetable soil with different years

表4 不同年限重金屬風(fēng)險評價Table4 Risk assessment of heavy metals with different years

2.5 土壤重金屬各形態(tài)與土壤理化性質(zhì)間的相關(guān)性分析

對不同年限設(shè)施菜地土壤中重金屬各形態(tài)之間及其與理化性質(zhì)等指標(biāo)進(jìn)行了相關(guān)性分析（表5和表6），設(shè)施菜地土壤中Cd酸提取態(tài)與可還原態(tài)、殘渣態(tài)、水溶態(tài)（F0）存在顯著性正相關(guān)關(guān)系，其中Cd酸提取態(tài)與可還原態(tài)及水溶態(tài)呈極顯著性相關(guān)，可還原態(tài)與殘渣態(tài)呈極顯著性相關(guān)，表明除可氧化態(tài)Cd外，設(shè)施菜地土壤中Cd的其他形態(tài)之間可相互轉(zhuǎn)化，其中酸提取態(tài)與可還原態(tài)在土壤中轉(zhuǎn)化能力較強(qiáng)，同時有效態(tài)Cd和水溶態(tài)Cd與速效磷、水解性氮呈正相關(guān)，可氧化態(tài)Cd與速效鉀呈極顯著相關(guān)，表明有效態(tài)和水溶態(tài)及可氧化態(tài)Cd極有可能來自設(shè)施菜地種植過程中化肥等肥料的使用[28]。設(shè)施菜地土壤中Pb酸提取態(tài)分別與可還原態(tài)、殘渣態(tài)及水溶態(tài)呈顯著性負(fù)相關(guān)，可還原態(tài)與殘渣態(tài)呈極顯著正相關(guān)，同時可還原態(tài)Pb和殘渣態(tài)Pb與速效磷、水解性氮呈正相關(guān)，表明可還原態(tài)和殘渣態(tài)來源相似，通過實地調(diào)查發(fā)現(xiàn)研究區(qū)設(shè)施菜地周圍無大型工業(yè)污染源，且用地下水灌溉，輸入源主要以磷肥和氮肥為主，因此可能來自設(shè)施菜地施用的肥料。重金屬Cd、Pb不同形態(tài)與土壤pH值的變化表現(xiàn)出顯著的負(fù)相關(guān)性，說明設(shè)施菜地土壤理化性質(zhì)改變，尤其是土壤酸化有利于重金屬Cd與Pb的活化，導(dǎo)致蔬菜暴露污染風(fēng)險增大[29]。

表5 Cd不同形態(tài)與理化性質(zhì)的相關(guān)性分析Table5 Correlation analysis of Cd speciation and physicochemical properties

表6 Pb不同形態(tài)與理化性質(zhì)的相關(guān)性分析Table6 Correlation analysis of Pb speciation and physicochemical properties

3 結(jié)論

（1）設(shè)施菜地土壤理化性質(zhì)（除干物質(zhì)外）隨著種植年限的增加變化顯著。隨著種植年限的延長，表現(xiàn)為設(shè)施菜地土壤中有效磷、有效鉀和水解性氮含量增加，有機(jī)質(zhì)先增加后減少，土壤呈酸化趨勢。

（2）西安市不同種植年限設(shè)施菜地土壤中重金屬Cd和Pb的總量均超過背景值，且隨著年限的增加，在使用設(shè)施菜地11 a后Cd總量超過標(biāo)準(zhǔn)限值，Cd的有效態(tài)含量與總量的趨勢一致，年限21～25 a的土壤超過標(biāo)準(zhǔn)限值。Pb的總量也呈增加趨勢，21～25 a的土壤Pb總量超過標(biāo)準(zhǔn)限值。

（3）隨著年限的增加，Cd酸提取態(tài)與可還原態(tài)之和大于可氧化態(tài)與殘渣態(tài)之和，表明設(shè)施菜地土壤Cd的生物有效性大，遷移的能力較強(qiáng)。而Pb殘渣態(tài)的比例較大，表明殘渣態(tài)Pb是大棚土壤中Pb的主要存在形態(tài)。

（4）隨著年限的增加，設(shè)施菜地土壤中Cd的地積累指數(shù)增加，污染程度由無污染到中度污染；而不同種植年限設(shè)施菜地土壤Pb無污染。設(shè)施菜地土壤中Cd達(dá)中等生態(tài)風(fēng)險；Pb為低生態(tài)風(fēng)險。各年限設(shè)施菜地的Cd生態(tài)風(fēng)險系數(shù)對潛在風(fēng)險指數(shù)的貢獻(xiàn)最高。西安市設(shè)施菜地Cd是不同年限設(shè)施菜地的主要風(fēng)險來源。