樊玉娜 孫小銀 高程程 王思嬋 廉心怡 解昕昕 劉一飛

摘要：本試驗以山東省白馬河沿岸耕地為研究對象，調(diào)查分析各耕地表層及剖面中重金屬（Cr、Cu、Mn、Ni、Pb、Zn）含量與碳氮含量的空間分布特征，并采用單因子污染指數(shù)（Ei）和潛在生態(tài)風(fēng)險指數(shù)（RI）對其進行風(fēng)險評價。研究結(jié)果顯示，該區(qū)耕地土壤Cr、Cu、Ni、Pb元素含量均低于國家土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)二級標(biāo)準(zhǔn)限制值，有個別剖面Zn元素含量略高，Mn元素低于《農(nóng)用地土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)（征求意見稿）》 標(biāo)準(zhǔn)值；耕地剖面潛在生態(tài)風(fēng)險指數(shù)平均值范圍在12.97～20.59之間，整體基本處于輕微生態(tài)風(fēng)險水平；耕作層有機碳含量較小，碳氮比介于0.90～12.64范圍內(nèi)且多數(shù)在10以下，表明土壤礦化作用強。在治理耕地土壤重金屬污染的同時，應(yīng)注意土壤中各元素含量的均衡，提高耕地土壤生產(chǎn)能力。

關(guān)鍵詞：耕地表層；耕地剖面；重金屬；碳氮比；空間分布；風(fēng)險評價

中圖分類號：S158.4：X53文獻標(biāo)識號：A文章編號：1001-4942（2018）04-0052-07

Abstract Taking the cultivated lands along the Baima River in Shandong Province as the research object， the spatial distribution characteristics of heavy metals（including Cr， Cu， Mn， Ni， Pb， Zn）， carbon and nitrogen in the surface and profiles were surveyed and analyzed. Their risks were evaluated by adopting the single factor pollution index （Ei） and potential ecological risk index （RI）. The results showed that the contents of Cr， Cu， Ni and Pb in cultivated soil were lower than the second grade standard limit values of the National Soil Environmental Quality Standard. The content of Zn in individual profiles were slightly higher and the content of Mn was lower than the level of the Agricultural Land Soil Environment Quality Standards（Exposure Draft）. The potential ecological risk indexes of cultivated land profiles ranged from 12.97 to 20.59， and they were in the slight ecological risk as a whole. The organic carbon content of tilling layer was small. The C/N ratios were in the range of 0.90～12.64 and most were lower than 10， which indicated the strong soil mineralization. In the controll of heavy metal pollution for cultivated land， we should simultaneously pay attention to the balance of elements in soil and improve the productivity of cultivated soil.

Keywords Cultivated land surface； Cultivated land profile； Heavy metal； C/N ratio； Spatial distribution； Risk assessment

土壤是人類賴以生存的自然環(huán)境和農(nóng)業(yè)生產(chǎn)的重要資源，世界面臨的糧食、資源和環(huán)境問題都與土壤密切相關(guān)[1]。隨著我國經(jīng)濟社會的快速發(fā)展，特別是近20年來工業(yè)化、城市化和農(nóng)業(yè)現(xiàn)代化發(fā)展的過程中，土壤質(zhì)量發(fā)生了較大變化，土壤環(huán)境污染問題呈現(xiàn)多樣化、復(fù)雜化和區(qū)域性的發(fā)展態(tài)勢，給我國經(jīng)濟社會發(fā)展帶來了新威脅和挑戰(zhàn)[2]。

目前，國內(nèi)外關(guān)于土壤重金屬污染的研究很多，但大多是土壤重金屬富集遷移研究和土壤重金屬污染現(xiàn)狀及修復(fù)，缺乏對耕地剖面重金屬污染的空間分布特征分析及風(fēng)險評價。分析耕地土壤表層及其剖面的重金屬污染分布與污染風(fēng)險，可以掌握耕地土壤中重金屬污染狀況，并針對研究結(jié)果進行有效的改良，進而確保糧食的安全生產(chǎn)和人體的健康。本研究對山東省白馬河沿岸耕地土壤中Cr、Cu、Mn、Ni、Pb、Zn 六種重金屬元素及其碳氮進行了分析，并采用單因子污染指數(shù)法和潛在生態(tài)風(fēng)險指數(shù)法對剖面進行了風(fēng)險評價，從而為該流域耕地的合理使用提供基礎(chǔ)數(shù)據(jù)和科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

白馬河流域位于山東省濟寧市境內(nèi)，流經(jīng)曲阜、兗州、鄒城、微山4縣（市），流域面積約 936.02 km2。白馬河在地貌上處于魯中南丘陵向魯西南平原的過渡地帶，地形起伏較大，是南四湖主要的入湖河流之一。較大支流有十里溝、石墻河、望云河、大沙河等。流域位于溫帶大陸性季風(fēng)氣候區(qū)，降水多集中在 7—9 月，徑流量年內(nèi)分配不均勻；同時，人類活動作用頻繁，對流域徑流過程產(chǎn)生了不同程度的影響[3]。

1.2 樣品的采集

2016年6月份對山東省白馬河流域進行野外考察和耕地剖面（0～100 cm）的土壤樣品采集，運用土鉆按土壤發(fā)生層每10 cm連續(xù)采集，在10個采樣點共采集97個樣品。土壤剖面樣品的采集點均布設(shè)在白馬河沿岸，采樣點位置如圖1所示。

1.3 樣品的處理方法

將除去石塊等雜物的土壤樣品自然晾干。然后研磨過100目尼龍網(wǎng)篩將其儲存于密封袋中待檢測。

土壤樣品預(yù)處理采用傳統(tǒng)的電熱板消解法。稱取0.5 g樣品放入聚乙烯坩堝中，然后加入氫氟酸5 mL，硝酸6 mL，鹽酸12 mL混合消解。將電熱板調(diào)至230℃放入通風(fēng)廚內(nèi)，待坩堝內(nèi)土樣完全消解后取下，定容至50 mL比色管中并搖勻。采用原子吸收分光光度計分別測定樣品中Cr、Cu、Mn、Ni、Pb、Zn 六種重金屬元素含量。上述試劑均采用優(yōu)級純，溶液配制用去離子水。土壤中C、N含量采用元素分析儀（Vario EL，德國Elementar Analysen systeme GmbH） 測定。

1.4 重金屬污染風(fēng)險評價方法

以國家《土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)》一級標(biāo)準(zhǔn)[4]和《農(nóng)用地土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)（征求意見稿）》標(biāo)準(zhǔn)值[5]為參比值，采用單因子污染指數(shù)（Ei）[6] 和潛在生態(tài)風(fēng)險指數(shù)（RI）評價耕地表層和剖面土壤重金屬污染的潛在生態(tài)風(fēng)險[7]。因Hakanson[6]研究的是8種有機污染物和重金屬的潛在生態(tài)風(fēng)險，而RI的大小與參評污染物的種類和數(shù)量有關(guān)，污染物的數(shù)目越多、毒性越強，RI 值就越大，因此，應(yīng)用 RI 進行生態(tài)風(fēng)險評價時，必須根據(jù)參評污染物的種類和數(shù)量對其進行調(diào)整，本研究根據(jù)文獻[8]中的方法對分級標(biāo)準(zhǔn)做了相應(yīng)調(diào)整[9]，調(diào)整后的重金屬潛在生態(tài)危害指數(shù)評價標(biāo)準(zhǔn)如表1所示。

2 結(jié)果與分析

2.1 白馬河沿岸耕地表層重金屬含量分布

由表2可知，白馬河沿岸耕作層（0～30 cm，下同）土壤樣品6種重金屬含量除S1、S3、S8和S9采樣點Zn略高于國家《土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)》二級標(biāo)準(zhǔn)值以外，其余5種重金屬元素采樣點監(jiān)測值均未超標(biāo)，這說明研究區(qū)土壤對重金屬還有一定的環(huán)境容量[10]。10個采樣點耕作層之間重金屬鎳含量有較大的區(qū)域差異，受人為因素影響的可能性大，其他5種重金屬含量區(qū)域變化小。研究區(qū)各樣點耕作層重金屬平均含量略有差異：S1（白馬河源頭）、S2、S3、S4、S7、S10采樣點重金屬平均含量是Mn>Zn>Cr>Pb>Cu>Ni；S5平均含量為Mn>Cr>Zn>Pb>Cu>Ni；S6平均含量為Mn>Cr>Zn>Cu>Pb>Ni；S8和S9兩地的耕作層重金屬平均含量Mn>Zn>Cr>Cu>Pb>Ni。這些差異可能由于耕作層年年翻作導(dǎo)致不同地點不同重金屬含量無明顯變化規(guī)律。

具體分析表2數(shù)據(jù)發(fā)現(xiàn)，Cr、Cu、Mn、Pb 4種重金屬在不同耕作層的含量最大相差約2倍；Zn的區(qū)域變化大約是3倍；而Ni在所有重金屬中含量最少，但在各耕作層上含量相差9倍還多。造成不同耕作層重金屬含量不等的原因：污水灌溉、肥料和農(nóng)藥的施用、大氣沉降及其他農(nóng)業(yè)活動。

采用單因子指數(shù)法對研究區(qū)耕作層土壤中的6種重金屬污染狀況進行評價[11]。由于耕作層是農(nóng)作物主要生長環(huán)境，為了更好的了解該研究區(qū)耕地是否符合農(nóng)業(yè)生產(chǎn)要求，研究采用《土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)》二級標(biāo)準(zhǔn)（二級標(biāo)準(zhǔn)：為保障農(nóng)業(yè)生產(chǎn)，維護人體健康的土壤限制值）作為參比值[12]，將污染程度劃分為4個等級：清潔（P<1）、輕度污染（1 ≤P<2）、中度污染（2≤P<3）、重度污染（P≥ 3）。研究結(jié)果如表3所示，S8采樣點的Zn污染指數(shù)為2，達到中度污染。S2、S8、S9采樣點的Cu，S1、S2、S4、S7、S8、S10采樣點的Pb以及S1、S2、S3、S7、S9、S10采樣點的Zn污染指數(shù)均大于1，呈現(xiàn)輕度污染水平。其余采樣點耕作層土壤中重金屬污染程度均屬于清潔水平。

2.2 流域耕地表層C、N含量分布

土壤有機碳和氮是土壤養(yǎng)分的重要組成部分，也是植物生長發(fā)育過程中的必需元素[13]。土壤中碳氮元素含量超標(biāo)或不達標(biāo)都會使植物體內(nèi)營養(yǎng)物質(zhì)不平衡。碳氮在土壤中分解和轉(zhuǎn)化，把復(fù)雜的碳素物質(zhì)變成簡單易溶，易于根系和氣孔的吸收[14]。分析研究區(qū)土壤中C、N數(shù)據(jù)發(fā)現(xiàn)，不同樣點耕地表層平均C含量差異明顯，依次為：S9>S7>S2>S1>S6>S3>S10>S5>S8>S4。而不同耕地表層平均N含量差異較小，依次為：S2>S7>S5=S10>S3>S6=S9>S4>S8>S1。

土壤碳氮比的高低對植物生長發(fā)育具有重要影響。如果碳氮比較低，則有利于微生物在有機質(zhì)分解過程中的養(yǎng)分釋放，促進土壤中有效氮的增加；反之，碳氮比較高，則會出現(xiàn)微生物在分解有機質(zhì)的過程中存在氮受限，從而與植物競爭土壤無機氮，不利于植物的生長[15]。由表4可知，山東省白馬河沿岸耕地表層碳氮比范圍在0.90～12.64之間，不同耕地表層碳氮比差異顯著。不同耕地表層碳氮比大小順序依次為：S1>S9>S2>S7>S6>S3>S10>S5>S8>S4。

2.3 沿岸耕地剖面重金屬含量分布特征

不同采樣點土壤剖面中（0～100 cm ，其中S8采樣點土壤剖面深度為0～70 cm）重金屬元素含量垂直分布如圖2所示，可以看出，研究區(qū)各剖面Cr、Cu、Ni、Pb重金屬含量均未超過國家《土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)》二級標(biāo)準(zhǔn)；Mn含量低于《農(nóng)用地土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)（征求意見稿）》標(biāo)準(zhǔn)值；Zn含量除S3（20～30 cm）和S8（10～20 cm）分別為247.55、249.4 μg/g，超出二級標(biāo)準(zhǔn)值外，其余剖面均未超標(biāo)。各剖面重金屬含量基本符合Mn>Zn> Cr， Cu、Ni、Pb三種重金屬含量無明顯差異。

分析各采樣點剖面重金屬垂直變化趨勢（圖2）可以看出，在0～30 cm土壤剖面上，多數(shù)采樣點重金屬含量明顯增加，原因在于表層土壤重金屬隨著淋溶作用下滲并富集于地下20～30 cm處。土壤剖面30 cm以下各段變化沒有明顯規(guī)律，這與各采樣點土壤性質(zhì)差異有關(guān)。下游S10采樣點Zn含量在30 cm以下剖面呈明顯的下降趨勢，其余重金屬含量也同樣隨著土壤剖面深度的增加而呈現(xiàn)出遞減趨勢。

2.4 土壤C、N的剖面分布

總體上來看，耕地表層有機碳質(zhì)量分數(shù)比30 cm以下深度土壤有機碳質(zhì)量分數(shù)相對較高，主要原因可能是耕作表層農(nóng)作物秸稈還田和有機肥的施用。大部分采樣點（S4、S5、S6、S7、S8、S10）土壤有機碳含量隨著剖面深度的增加呈明顯的下降趨勢。不同土壤剖面有機碳含量變化總體呈波動型。土壤剖面中氮質(zhì)量分數(shù)變化總體呈波動型，變化幅度較小。分析圖3發(fā)現(xiàn)，流域內(nèi)各耕地剖面碳氮比依次為：S2>S1>S9>S3>S7>S6>S10>S5>S8>S4。碳氮比在水平分布上差異顯著。而且表土層與底土層土壤碳氮比存在顯著差異。

2.5 耕地重金屬污染風(fēng)險評價

以國家土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)和山東省土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)作參比值，采用Hakanson潛在生態(tài)風(fēng)險指數(shù)法，計算得到白馬河沿岸耕地土壤重金屬在各樣點剖面的單項潛在生態(tài)風(fēng)險指數(shù)（Ei）及綜合潛在生態(tài)風(fēng)險指數(shù)（RI），根據(jù)潛在生態(tài)風(fēng)險分級標(biāo)準(zhǔn)進行了生態(tài)風(fēng)險評價[7]。重金屬元素Cr、Cu、Mn、Ni、Pb和Zn的毒性系數(shù)分別取2、5、1、5、5和1[6]。如表5所示，70%采樣點Cu、Pb單項潛在生態(tài)風(fēng)險指數(shù)介于5～10之間，30%采樣點的指數(shù)大于4，表明研究區(qū)Cu、Pb處于中等生態(tài)危害水平，其余重金屬均屬輕微生態(tài)危害水平。由表6看出，綜合潛在生態(tài)風(fēng)險指數(shù)平均值范圍介于12.97～20.59，S8潛在生態(tài)風(fēng)險指數(shù)平均值大于20，呈現(xiàn)中等生態(tài)危害水平，研究區(qū)其余采樣點基本屬于輕微生態(tài)風(fēng)險水平。

由表6還可看出，S2取樣點0～10 cm土層和S8部分土層RI值大于20，表明這幾個土層生態(tài)風(fēng)險指數(shù)較高，均屬中等生態(tài)風(fēng)險水平。其余研究點雖處于輕微生態(tài)風(fēng)險，但存在個別生態(tài)風(fēng)險指數(shù)較高的剖面，RI值超過19的占到11.34%。沿岸耕地平均潛在生態(tài)風(fēng)險指數(shù)大小為：S8>S9>S1>S2>S7>S5>S6>S4>S3>S10。

3 結(jié)論

白馬河沿岸耕地土壤重金屬Cr、Cu、Ni、Pb含量低于《土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)》二級標(biāo)準(zhǔn)值，有個別剖面Zn含量略高；Mn含量低于《農(nóng)用地土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)（征求意見稿）》標(biāo)準(zhǔn)值。該研究區(qū)整體看符合國家對農(nóng)業(yè)生產(chǎn)的要求。由于耕地表層農(nóng)業(yè)活動頻繁，導(dǎo)致不同耕地各種重金屬含量分布特征不明顯。大部分研究點耕作層重金屬平均含量大小為Mn>Zn>Cr>Pb>Cu>Ni，Mn和Ni在各樣點均是含量最高和最低的元素。耕地潛在生態(tài)風(fēng)險指數(shù)平均值范圍介于12.97～20.59，基本處于輕微生態(tài)風(fēng)險水平，個別剖面處于中等生態(tài)風(fēng)險水平，還有多個剖面處于輕微生態(tài)風(fēng)險最高值。為保持該研究區(qū)優(yōu)良的生態(tài)環(huán)境，需加強對中等生態(tài)風(fēng)險水平的土壤進行治理。

不同土壤耕作層的有機碳質(zhì)量分數(shù)基本都比耕地剖面底層的有機碳質(zhì)量分數(shù)相對較高。不同耕地剖面C、N質(zhì)量分數(shù)變化總體呈波動型。流域內(nèi)不同耕地表層碳氮比依次為：S1>S9>S2>S7>S6>S3>S10>S5>S8>S4。流域內(nèi)耕地剖面碳氮比水平分布差異顯著，碳氮比依次為：S2>S1>S9>S3>S7>S6>S10>S5>S8>S4。

參 考 文 獻：

[1] 李福燕. 海南島農(nóng)用地重金屬污染現(xiàn)狀調(diào)查與評價研究[D]. ?？冢汉Ｄ洗髮W(xué)， 2010.

[2] 劉曉慧.“土十條”將出“新藥”能否根治頑疾？——對我國土壤污染治理發(fā)展趨勢的解讀[N].中國礦業(yè)報，2016-03-01（6）.

[3] 王學(xué)， 張祖陸， 寧吉才. 基于SWAT模型的白馬河流域土地利用變化的徑流響應(yīng)[J]. 生態(tài)學(xué)雜志， 2013， 32（1）： 186-194.

[4] 環(huán)境保護部，國家質(zhì)量監(jiān)督檢驗檢疫總局.GB15618-1995 土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)（修訂）[S]. 1995-07-13.

[5] 環(huán)境保護部，國家質(zhì)量監(jiān)督檢驗檢疫總局. GB15618-201□ 農(nóng)用地土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)（征求意見稿）[S].2015-01-13.

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

[7] 麥麥提吐爾遜·艾則孜， 阿吉古麗·馬木提， 艾尼瓦爾·買買提. 新疆焉耆盆地辣椒地土壤重金屬污染及生態(tài)風(fēng)險預(yù)警[J]. 生態(tài)學(xué)報， 2018，38（3）： 1-12.

[8] 李一蒙， 馬建華， 劉德新， 等. 開封城市土壤重金屬污染及潛在生態(tài)風(fēng)險評價[J]. 環(huán)境科學(xué)， 2015，36（3）： 1037-1044.

[9] 陳鳳， 董澤琴， 王程程， 等. 鋅冶煉區(qū)耕地土壤和農(nóng)作物重金屬污染狀況及風(fēng)險評價[J]. 環(huán)境科學(xué)， 2017，38（10）： 4360-4369.

[10]廖沖， 曾凡萍， 劉澍. 萍鄉(xiāng)市農(nóng)用土壤重金屬含量及其分布特征分析[J]. 中國環(huán)境管理干部學(xué)院學(xué)報， 2015，25（1）： 62-66.

[11]周俊馳， 鐵柏清， 劉孝利， 等. 湖南礦區(qū)縣域耕地重金屬污染空間特征及潛在風(fēng)險評價[J]. 湖南農(nóng)業(yè)科學(xué)， 2017（4）： 75-79.

[12]汪雅各. 農(nóng)業(yè)環(huán)境標(biāo)準(zhǔn)實用手冊[M]. 杭州：浙江大學(xué)出版社， 1991：32.

[13]林麗， 張法偉， 李以康， 等. 高寒矮嵩草草甸退化過程土壤碳氮儲量及C/N化學(xué)計量學(xué)特征[J]. 中國草地學(xué)報， 2012，34（3）： 42-47.

[14]胡濟生， 梁德印. 關(guān)于小麥碳氮營養(yǎng)關(guān)系的商討[J]. 中國農(nóng)業(yè)科學(xué)，1961，2（5）： 26-28.

[15]王建林， 鐘志明， 王忠紅， 等. 青藏高原高寒草原生態(tài)系統(tǒng)土壤碳氮比的分布特征[J].生態(tài)學(xué)報， 2014，34（22）： 6678-6691.