郭展翅,周葆華*,2,趙 寬,胡睿鑫,劉 勤
(1.安慶師范大學(xué)資源環(huán)境學(xué)院,安徽安慶246133;2.皖江流域水環(huán)境保護(hù)與污染控制安徽省教育廳重點實驗室,安徽安慶246133)
鎘(Cd)是有毒有害的、人體非必須的微量元素,其污染具有隱蔽、非逆轉(zhuǎn)和易累積等特性。Cd排放總量在國內(nèi)外都被列為重點控制指標(biāo)[1]。Cd主要由工業(yè)生產(chǎn)、礦山開發(fā)、金屬冶煉、交通運(yùn)輸和現(xiàn)代農(nóng)業(yè)生產(chǎn)等過程產(chǎn)生,并通過各種途徑進(jìn)入農(nóng)田土壤中[2-5]。農(nóng)田土壤Cd污染不僅破壞農(nóng)田生態(tài)系統(tǒng),而且影響農(nóng)作物產(chǎn)品質(zhì)量安全。Cd污染已經(jīng)引起國內(nèi)外政府和科學(xué)家密切關(guān)注[6-7]。我國農(nóng)田土壤Cd的超標(biāo)率為7%,是我國耕地土壤中超標(biāo)率最高的有害重金屬,農(nóng)田土壤Cd污染形勢嚴(yán)峻[8-10]。長江中下游平原是中國水資源最豐富的區(qū)域,土地利用類型主要以農(nóng)業(yè)用地和林業(yè)用地為主,主要作物有水稻、小麥、油菜等。研究長江中下游農(nóng)田Cd含量的變化特征及生態(tài)風(fēng)險具有重要意義。本課題組選擇長江中下游某地典型農(nóng)田為研究對象,采集了33個表層土壤樣品,使用HNO3-HClO4-HF-HCl消解的方法測定了土壤樣品中Cd的全量,基于Tessier五步連續(xù)提取法提取Cd的5種形態(tài),分析Cd的空間分布特征;采用單項污染指數(shù)法和地累積指數(shù)法揭示區(qū)域內(nèi)農(nóng)田土壤Cd的污染程度,運(yùn)用潛在生態(tài)風(fēng)險指數(shù)法評價研究區(qū)域的生態(tài)風(fēng)險[11-13],為該區(qū)域農(nóng)田土壤Cd污染的監(jiān)測和防控提供理論支持。
運(yùn)用2.5 km×2.5 km網(wǎng)格布點,在研究區(qū)內(nèi)利用GPS現(xiàn)場定位和標(biāo)號,選取33個點位,每個定位的采樣點按田字形分別用采樣器在4個頂點和中間采取5個(0~20)cm表層土樣混均成一個農(nóng)田土壤樣品,將混合的土壤放入具有密封性的聚乙烯自封塑料袋內(nèi)密封,然后編號,并做好現(xiàn)場采樣記錄,帶回實驗室準(zhǔn)備做預(yù)處理[14]。土壤主要分為水田和旱地,33個樣品中有17份水田土壤、16份旱地土壤,分別命名為SH1、SH2、SH3、···、SH33。
土壤樣品前處理參照《土壤環(huán)境檢測技術(shù)規(guī)范》(HJ166-2004),采集的樣品放入實驗室中,在室溫條件下自然風(fēng)干,一般為14 d左右。后用干燥木錘、木棍研磨,并除去樣品中植物根系、石子等雜物,過100目篩裝瓶保存于干燥陰涼處。準(zhǔn)確稱取0.1 g樣品于微波消解罐中,經(jīng)HNO3-HClO4-HF-HCl高溫消解后,使用ICP-AES(Thermo Fisher Scientific ICAP6300)測定Cd全量。運(yùn)用標(biāo)準(zhǔn)土樣GSF-3為參照,驗證Cd回收方式的精準(zhǔn)度,實驗中Cd的回收率為97.90%。Tessier五步連續(xù)提取法:使用氯化鎂提取土壤樣品中Cd的可交換態(tài),乙酸鈉和乙酸提取碳酸鹽結(jié)合態(tài),使用鹽酸羥胺與乙酸提取錳氧化物結(jié)合態(tài),氧化氫和硝酸提取有機(jī)結(jié)合態(tài),運(yùn)用微波消解-ICP測定殘渣態(tài)[15],詳細(xì)參數(shù)見表1。
表1 不同形態(tài)Cd的連續(xù)提取法
單項污染指數(shù)法。單項污染指數(shù)法是國內(nèi)外運(yùn)用最為廣泛的判定污染程度方法之一[16]。根據(jù)農(nóng)田重金屬Cd的判定標(biāo)準(zhǔn)值來判斷采樣點位的污染程度,本文采用《土壤環(huán)境質(zhì)量農(nóng)用地污染風(fēng)險管控標(biāo)準(zhǔn)(試行)》(GB15618-2018)中Cd的篩選值(表2),以便直觀地揭示農(nóng)田土壤Cd的污染狀況。第i個采樣點重金屬Cd的單項污染指數(shù)Pi=ci/S,其中ci為第i個重金屬Cd的測試全量值(mg/kg),S為農(nóng)田重金屬Cd是否超標(biāo)的判定標(biāo)準(zhǔn)值。單項污染指數(shù)法Cd單項污染指數(shù)的分級標(biāo)準(zhǔn)見表3。
表2 農(nóng)田土壤Cd風(fēng)險篩選值
表3 Cd單項污染指數(shù)的分級標(biāo)準(zhǔn)
地累積指數(shù)法。地累積指數(shù)法不僅將人為造成的Cd污染、土壤Cd地球化學(xué)背景值考慮在內(nèi),還將地殼運(yùn)動與遷移過程中引入的Cd等因素綜合考慮在內(nèi)[17],公式為
式(1)中,Igeo為地累積指數(shù);ca為第a個采樣點Cd在土壤中的含量值(mg/kg);B為Cd為在土壤中的地球化學(xué)背景值(mg/kg),本研究1.5為地殼運(yùn)動與遷移過程中引入的Cd背景值的變動系數(shù)。以安徽省江淮流域土壤元素Cd的背景值(0.104 mg/kg)來表示[18],Cd地累積指數(shù)(Igeo)的污染程度分級如表4所示。
潛在生態(tài)風(fēng)險指數(shù)法。潛在生態(tài)風(fēng)險指數(shù)法是國內(nèi)外最常用的生態(tài)評價方法之一,其不僅考慮到重金屬Cd的全量,還引入重金屬Cd毒性系數(shù)(取30),使對研究區(qū)Cd的評價更具有可比性和等價性[18-19]。Cd潛在生態(tài)風(fēng)險指數(shù)法公式為
式(2)中,Er為Cd的潛在風(fēng)險系數(shù);ck為第k個采樣點的重金屬Cd的實測含量值(mg/kg);cp為重金屬Cd的研究區(qū)的土壤背景值,本研究選取安徽省江淮流域土壤元素Cd的背景值[11]。Cd的潛在生態(tài)風(fēng)險指數(shù)的評價標(biāo)準(zhǔn)見表5。
表4 Cd地累積指數(shù)(Igeo)的污染程度分級
表5 Cd潛在生態(tài)風(fēng)險指數(shù)的評價標(biāo)準(zhǔn)
表6為兩種農(nóng)田類型土壤Cd含量基本統(tǒng)計。由表6可知,17個水田土壤樣品重金屬Cd的均值為0.61 mg/kg,超過其最高風(fēng)險管控值,是江淮流域背景值的6倍;16個旱地土壤樣品重金屬Cd平均值為0.66 mg/kg,是其背景值的6.3倍,風(fēng)險管控值的3.3倍;水田土壤pH在5.27~6.8之間,旱地土壤pH在5.06~7.36之間。17個水田采樣點中有10個點位超過Cd的篩選值,16個旱地采樣點中有12個點位超過其篩選值,33個樣品總超標(biāo)率為72.73%。變異系數(shù)可反映樣品的離散程度,本研究區(qū)域內(nèi)農(nóng)田土壤Cd總量的變異系數(shù)為64%,屬于中等變異,說明研究區(qū)內(nèi)農(nóng)田土壤Cd分布不均勻,易受外界因素影響[20]。
利用ArcGis10.6軟件分析了Cd的空間分布特征,如圖1所示。結(jié)果表明,農(nóng)田土壤Cd的濃度從高到低依次為北部、東部、西部,濃度高的區(qū)域主要分布在沿江區(qū),這可能與沿江區(qū)域電鍍、塑料、有色金屬冶煉等企業(yè)將含Cd廢物排入土壤或者長江水體有關(guān)[21-22]。
表6 兩種農(nóng)田類型土壤Cd含量基本統(tǒng)計
圖2為研究區(qū)農(nóng)田土壤重金屬Cd形態(tài)值。由圖2可知,研究區(qū)33個農(nóng)田土壤樣品Cd化學(xué)形態(tài)變化趨勢從高到低依次為可交換態(tài)(36.52%)、鐵錳氧化物結(jié)合態(tài)(35.57%)、殘渣態(tài)(19.41%)、碳酸鹽結(jié)合態(tài)(5.28%)、有機(jī)結(jié)合態(tài)(4.71%)??山粨Q態(tài)是對生物營養(yǎng)或毒害影響最關(guān)鍵的形態(tài),具有較強(qiáng)的遷移性,容易被植物吸收,生物毒性較強(qiáng)[23]。農(nóng)田土壤Cd交換態(tài)高說明在該研究區(qū)內(nèi)農(nóng)田重金屬Cd容易被農(nóng)作物吸收,容易通過食物鏈進(jìn)入人體,從而對人體健康造成危害[24]。
圖1 研究區(qū)域重金屬Cd濃度空間分布
圖2 研究區(qū)農(nóng)田土壤重金屬Cd形態(tài)分析
單項污染指數(shù)法。對區(qū)域內(nèi)33個采樣點單項污染評價統(tǒng)計見表7,17個水田土壤采樣點位整體表現(xiàn)為輕度污染,16個旱地整體表現(xiàn)為中度污染,所采集點位中Cd單項污染指數(shù)為0.46~5.07,平均值為1.68,采樣區(qū)域整體上處于輕度污染。
利用ArcGis10.6對研究區(qū)33個采樣點重金屬Cd單項污染指數(shù)的分級標(biāo)準(zhǔn)進(jìn)行了反距離加權(quán)插值,空間分布特征見圖3。研究區(qū)內(nèi)絕大部分區(qū)域?qū)儆谳p度污染,北部沿江區(qū)域?qū)儆谥卸任廴?,這可能與北部沿江區(qū)煤、石油、金屬冶煉等企業(yè)有關(guān)[25-26]。
表7 研究區(qū)Cd單項污染指數(shù)的評價結(jié)果
地累積指數(shù)法。地累積指數(shù)法分析結(jié)果如表8所示。由表8可知,水田、旱地土壤整體污染程度為偏中度污染,33個采樣點位地累積指數(shù)處于-0.12~3.72,平均值為1.85,研究區(qū)域整體上處于偏中度污染。使用ArcGis10.6軟件對研究區(qū)33個采樣點重金屬Cd地累積指數(shù)的分級標(biāo)準(zhǔn)進(jìn)行反距離加權(quán)插值,地累積指數(shù)法污染空間分布特征見圖4,該區(qū)域采樣點位為偏中度污染,這與表8得出的結(jié)論一致。偏重度污染的點位處于沿江區(qū)域,可能是沿江區(qū)域農(nóng)業(yè)活動和電池業(yè)、煤、石油、金屬冶煉等企業(yè)在生產(chǎn)或排放過程中將含鎘的污染源排入長江有關(guān)[26-27]。
表8 研究區(qū)Cd土壤地累積指數(shù)的評價結(jié)果
圖3 農(nóng)田土壤重金屬Cd單項污染指數(shù)風(fēng)險評價的空間分布
圖4 農(nóng)田土壤重金屬Cd潛在風(fēng)險指數(shù)法風(fēng)險評價的空間分布
潛在生態(tài)風(fēng)險指數(shù)法。潛在生態(tài)風(fēng)險指數(shù)法分析結(jié)果如表9所示。由表9可知,水田、旱地土壤整體上均處于高等生態(tài)危害,采樣區(qū)域內(nèi)33個點位Cd潛在生態(tài)風(fēng)險系數(shù)Er的范圍為41.48~594.2,處于中高等至高等危害水平,平均值為148.12,整體上為高等生態(tài)危害,該研究區(qū)內(nèi)所有采樣點都存在不同程度的生態(tài)危害。為進(jìn)一步對研究區(qū)進(jìn)行潛在生態(tài)風(fēng)險評價,使用ArcGis10.6軟件對研究區(qū)33個采樣點重金屬Cd潛在風(fēng)險系數(shù)的分級標(biāo)準(zhǔn)進(jìn)行反距離加權(quán)插值做出空間分布圖(圖5),研究區(qū)中部出現(xiàn)個別極高生態(tài)危害的點位,這可能與該采樣點位的金屬開采、冶煉等有關(guān)[30]。
表9 研究區(qū)Cd潛在生態(tài)風(fēng)險指數(shù)的評價結(jié)果
圖5 農(nóng)田土壤重金屬Cd地累積指數(shù)風(fēng)險評價的空間分布
綜上所述,研究區(qū)內(nèi)農(nóng)田土壤pH介于弱堿性和酸性之間,水田重金屬Cd的平均值為0.61 mg/kg,旱地重金屬Cd的平均值為0.66 mg/kg,均高于其最高風(fēng)險管控值。水田土壤重金屬Cd的超標(biāo)率為64.7%,旱地土壤重金屬Cd的超標(biāo)率為75%,重金屬Cd的濃度從高到低依次為北部、東部、西部,濃度高的區(qū)域主要呈現(xiàn)在沿江區(qū)域。研究區(qū)內(nèi)農(nóng)田土壤重金屬Cd主要以可交換態(tài)存在。
研究區(qū)內(nèi)不同采樣點位均呈現(xiàn)不同程度的生態(tài)風(fēng)險,單項污染指數(shù)法結(jié)果表明,采集點位Cd處于輕污染;地累積指數(shù)法表明,研究區(qū)Cd為污染程度為偏中度,且沿江區(qū)域為中度污染,其余區(qū)域為無污染至輕度污染;潛在生態(tài)風(fēng)險指數(shù)法結(jié)果表明,研究區(qū)整體表現(xiàn)為中高等至高等生態(tài)危害,部分采樣點位存在極高的生態(tài)危害,應(yīng)當(dāng)引起足夠的重視。造成農(nóng)田土壤污染的原因可能與研究區(qū)域長江中下游有色金屬的冶煉、礦石的煅燒、電池業(yè)等企業(yè)生產(chǎn)排放含鎘廢物至土壤和水體有關(guān)。對于Cd污染比較嚴(yán)重的地區(qū),應(yīng)當(dāng)調(diào)整種植結(jié)構(gòu),改種非食用植物,或者通過施加鈍化劑等技術(shù),改善土壤中Cd的濃度,使農(nóng)產(chǎn)品達(dá)到安全食用的標(biāo)準(zhǔn)。