諶婕妤 楊思雨 胡愛(ài)玲 許峰 黃運(yùn)新

摘要：農(nóng)田氮磷排放是農(nóng)業(yè)面源污染的重要來(lái)源之一，其排放潛力可能受到諸多環(huán)境因素的影響。為了解農(nóng)田氮磷排放潛力的影響因素，選點(diǎn)湖北省武漢市新洲區(qū)西黃村，利用箱線(xiàn)圖和方差分析方法進(jìn)行研究。結(jié)果表明，研究地點(diǎn)土壤氮含量的平均水平為0.974～1.297 g/kg，其中秋季水田含量相對(duì)較高，而夏季旱地含量相對(duì)較低。土壤磷含量平均水平為0.491～0.676 g/kg，其中夏季旱地近居民區(qū)含量相對(duì)較高，而秋季水田遠(yuǎn)居民區(qū)含量相對(duì)較低。氮含量在季節(jié)、水旱作之間存在顯著差異，而在居民區(qū)遠(yuǎn)近之間不存在顯著差異;磷含量在季節(jié)、水旱作和居民區(qū)遠(yuǎn)近之間均存在顯著差異。

關(guān)鍵詞：農(nóng)田土壤;氮含量;磷含量;時(shí)空差異;面源污染

中圖分類(lèi)號(hào)：S151.9?文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：0439-8114（2020）07-0068-04

DOI：10.14088/j.cnki.issn0439-8114.2020.07.015

農(nóng)業(yè)面源污染（Agricultural non-point sources pollution，ANPSP）是指在農(nóng)業(yè)生產(chǎn)活動(dòng)中，氮、磷等營(yíng)養(yǎng)物質(zhì)以及農(nóng)藥等其他有害物質(zhì)通過(guò)雨水沖刷、地表徑流以及農(nóng)田灌溉而造成的水環(huán)境污染[1]。相對(duì)于點(diǎn)源污染而言，ANPSP來(lái)源廣、遷移復(fù)雜，極難準(zhǔn)確測(cè)定和控制[2，3]，而農(nóng)田氮磷排放是ANPSP的主要來(lái)源[4]。農(nóng)田土壤氮磷流失可導(dǎo)致水體富營(yíng)養(yǎng)化，對(duì)地表水造成污染[4]。過(guò)量的氮還可能滲入地下，導(dǎo)致硝酸鹽濃度過(guò)高，對(duì)地下水造成污染[5]。土壤氮磷流失的監(jiān)測(cè)與控制是ANPSP防控的重要課題之一。

過(guò)度施肥是導(dǎo)致土壤氮磷流失的主因[5]。在中國(guó)，為了獲得高產(chǎn)和豐收，農(nóng)業(yè)生產(chǎn)中廣泛存在著過(guò)量使用化肥的現(xiàn)象[4，6-9]。過(guò)度施肥不但造成浪費(fèi)，而且可能導(dǎo)致土壤養(yǎng)分失調(diào)和板結(jié)，影響土壤的長(zhǎng)期使用。更重要的是，過(guò)度施肥容易造成土壤養(yǎng)分流失，影響水體環(huán)境。

江漢平原位于湖北省中南部，是長(zhǎng)江中下游平原的重要組成部分[10，11]。江漢平原具有氣候適宜、水源豐富等得天獨(dú)厚的條件，這些條件對(duì)發(fā)展農(nóng)業(yè)生產(chǎn)十分有利，因此江漢平原不僅是湖北省而且是全國(guó)重要的糧、棉、油生產(chǎn)基地。但另一方面，江漢平原獨(dú)特的水源和水網(wǎng)條件又為ANPSP創(chuàng)造了便利。因此，江漢平原也是ANPSP的重災(zāi)區(qū)[12]。

近年來(lái)，隨著農(nóng)業(yè)技術(shù)的發(fā)展和土地集約化，農(nóng)業(yè)生產(chǎn)中的化肥使用強(qiáng)度和使用方式都在發(fā)生變化。這些變化不可避免地對(duì)土壤和水體環(huán)境造成影響。通過(guò)土壤監(jiān)測(cè)了解實(shí)際變化情況對(duì)ANPSP防控具有現(xiàn)實(shí)意義。為此，本研究以位于江漢平原東部的武漢市新洲區(qū)為例，通過(guò)抽樣調(diào)查了解土壤氮磷含量水平，為進(jìn)一步估算農(nóng)業(yè)面源污染狀況打下基礎(chǔ)。

1?材料與方法

1.1?土樣采集和指標(biāo)檢測(cè)

1.1.1?采樣地點(diǎn)與方法?采樣地點(diǎn)位于湖北省武漢市新洲區(qū)西黃村（圖1）。采樣在夏、秋兩季各進(jìn)行1次，每次采樣17份，按照近居民區(qū)和遠(yuǎn)居民區(qū)分別采集，在近居民區(qū)和遠(yuǎn)居民區(qū)中再分別選取旱地和水田采樣。每份土樣1 000 g，采樣深度為0～10 cm。采集的土樣用密封袋密封，帶回實(shí)驗(yàn)室進(jìn)行總氮、總磷的測(cè)定。

1.1.2?指標(biāo)檢測(cè)?土壤采回后攤開(kāi)風(fēng)干，剔除石塊和植物殘?bào)w，將其研磨后過(guò)60號(hào)篩（孔徑0.25 mm），取0.5～1.0 g（精確到0.001 g）土樣，放入150 mL開(kāi)氏瓶中用重鉻酸鉀-硫酸消化法測(cè)定總氮，再取1 g（精確到0.001 g）土樣置于50 mL三角瓶中，用硫酸-高氯酸消煮法測(cè)定總磷[13]。

1.2?數(shù)據(jù)分析

1.2.1?箱線(xiàn)圖?采用箱線(xiàn)圖分析了總氮、總磷2個(gè)指標(biāo)觀(guān)測(cè)值的變異情況。其中，箱線(xiàn)圖由箱和線(xiàn)組成，包括1個(gè)箱體和2條分別從箱體頂端向上和底端向下的垂直延長(zhǎng)線(xiàn)（觸須）。箱體下邊界代表第一四分位數(shù)，上邊界代表第三四分位數(shù)，箱體內(nèi)部的線(xiàn)段代表中位數(shù)，箱體長(zhǎng)度代表四分位距;觸須終點(diǎn)為觀(guān)測(cè)值的臨界點(diǎn)，代表最大或最小觀(guān)測(cè)值，臨界點(diǎn)之外為異常值。

1.2.2?方差分析?采用多因素方差分析分別對(duì)總氮和總磷2個(gè)指標(biāo)的空間、土壤類(lèi)型和季節(jié)差異以及三者之間的交互影響進(jìn)行了分析。多因素方差分析采用SPSS軟件進(jìn)行。

2?結(jié)果與分析

2.1?氮含量

2.1.1?空間差異?由表1、圖2可知，近居民區(qū)和遠(yuǎn)居民區(qū)總氮水平存在差異，但兩類(lèi)均沒(méi)有異常值，中位數(shù)分別為1.065、1.052 g/kg，均值分別為1.097、1.180 g/kg，標(biāo)準(zhǔn)差分別為0.082、0.100 g/kg。表2方差分析結(jié)果顯示，近居民區(qū)和遠(yuǎn)居民區(qū)的總氮含量不存在顯著差異（P=0.65）。這表明居民生活污水或畜禽養(yǎng)殖產(chǎn)生的氮源對(duì)周邊農(nóng)田土壤氮含量沒(méi)有顯著影響。

2.1.2?水田旱地差異?由表1、圖2可知，水田和旱地總氮水平存在差異，其中水田總氮沒(méi)有異常值，而旱地總氮存在3個(gè)異常值。旱地總氮存在異常值的原因可能是部分采樣點(diǎn)有過(guò)度施肥現(xiàn)象。水田和旱地總氮的中位數(shù)分別為1.244、0.912 g/kg，均值分別為1.244、0.974 g/kg，標(biāo)準(zhǔn)差分別為0.088、0.076 g/kg。由表2可知，將異常值包含在內(nèi)的方差分析顯示，總氮在水田和旱地之間存在顯著差異（P=0.05），說(shuō)明水田和旱地氮肥的施肥量和土壤殘留具有明顯差異。將異常值去掉后的方差分析顯示，總氮含量在水田和旱地之間差異更為顯著（P=0.00）。

2.1.3?季節(jié)差異?由表1、圖2可知，夏季和秋季總氮水平存在一定差異，兩類(lèi)均沒(méi)有異常值，中位數(shù)分別為0.912、1.339 g/kg，均值分別為0.985、1.297 g/kg，標(biāo)準(zhǔn)差分別為0.067、0.099 g/kg。表2方差分析結(jié)果顯示，總氮含量在夏季和秋季存在顯著差異（P=0.00），說(shuō)明夏季和秋季氮肥的施肥量存在明顯差異。

2.1.4?交互作用?表2方差分析結(jié)果顯示，土壤總氮含量在空間類(lèi)型和土地類(lèi)型之間存在顯著的交互作用（P=0.01），但是在空間類(lèi)型和季節(jié)之間、土地類(lèi)型和季節(jié)之間交互作用不顯著。

2.2?磷含量

2.2.1?空間差異?由表3、圖3可知，近居民區(qū)和遠(yuǎn)居民區(qū)總磷的含量存在差異，但兩類(lèi)均沒(méi)有異常值，中位數(shù)分別為0.619、0.435 g/kg，均值分別為0.656、0.507 g/kg，標(biāo)準(zhǔn)差分別為0.055、0.040 g/kg。表4方差分析結(jié)果顯示，總磷含量在近居民區(qū)和遠(yuǎn)居民區(qū)之間存在顯著差異（P=0.00），表明居民生活污水或畜禽養(yǎng)殖產(chǎn)生的磷源對(duì)周邊農(nóng)田磷含量有顯著影響。

2.2.2?水田旱地差異?由表3、圖3可知，水田和旱地總磷含量存在一定差異，但兩類(lèi)均沒(méi)有異常值，中位數(shù)分別為0.490、0.606 g/kg，均值分別為0.516、0.676 g/kg，標(biāo)準(zhǔn)差分別為0.033、0.068 g/kg。表4方差分析結(jié)果顯示，總磷含量在水田和旱地之間存在顯著差異（P=0.00）。

2.2.3?季節(jié)差異?由表3、圖3可知，夏季和秋季土壤總磷含量存在差異，其中夏季沒(méi)有異常值，而秋季總磷存在1個(gè)異常值。夏季和秋季總磷的中位數(shù)分別為0.644、0.431 g/kg，均值分別為0.663、0.491 g/kg，標(biāo)準(zhǔn)差分別為0.047、0.045 g/kg。由表4可知，將異常值包含在內(nèi)的方差分析顯示，總磷在夏季和秋季之間存在顯著差異（P=0.01）;將異常值去掉后的方差分析顯示，總磷在夏季和秋季之間依然存在顯著差異（P=0.01）。

2.2.4?多因素顯著性結(jié)果?由表4方差分析結(jié)果可知，土壤總磷含量在空間類(lèi)型和土地類(lèi)型之間存在顯著的交互作用（P=0.03），但是在空間類(lèi)型和季節(jié)之間、土地類(lèi)型和季節(jié)之間交互作用都不顯著。

3?小結(jié)與討論

研究結(jié)果表明，研究地點(diǎn)土壤氮含量平均水平為0.974～1.297 g/kg，其中秋季水田含量相對(duì)較高，夏季旱地含量相對(duì)較低，氮含量在季節(jié)和水旱作之間存在顯著差異，而在居民區(qū)遠(yuǎn)近之間不存在顯著差異。研究地點(diǎn)土壤磷含量平均水平為0.491～0.676 g/kg，其中夏季旱地近居民區(qū)含量相對(duì)較高，而秋季水田遠(yuǎn)居民區(qū)含量相對(duì)較低;磷含量在季節(jié)、水旱作和居民區(qū)遠(yuǎn)近之間均存在顯著差異。

由于水田旱地的施肥頻率和強(qiáng)度通常不一樣，作物吸收能力也存在差異，所以氮磷含量在水旱作之間存在差異[14，15]。而影響氮磷含量在季節(jié)之間差異的因素除了施肥強(qiáng)度和作物吸收能力外，還有土壤養(yǎng)分流失的程度，由于不同季節(jié)降雨量不同，土壤氮磷流失在季節(jié)之間存在差異[14，16，17]，必然導(dǎo)致土壤氮磷含量之間的差異。

除了土地類(lèi)型和季節(jié)因素外，離居民區(qū)的遠(yuǎn)近也有可能影響土壤氮磷含量，因?yàn)榭拷用駞^(qū)的農(nóng)田更有可能受到畜禽糞便和生活污水的影響。本研究結(jié)果顯示，離居民區(qū)遠(yuǎn)近對(duì)土壤磷含量影響顯著，對(duì)土壤氮含量影響不顯著。由于對(duì)畜禽糞便和生活污水中的污染物含量和進(jìn)入農(nóng)田的機(jī)制并不清楚，尚無(wú)法解釋這一現(xiàn)象。后續(xù)有必要結(jié)合農(nóng)田灌溉水的氮磷含量分析加以研究，弄清具體原因[15，18]。

土壤氮磷流失是ANPSP的主要來(lái)源之一。有研究顯示農(nóng)田土壤氮磷排放有3/4來(lái)自土壤本底，1/4來(lái)自當(dāng)年施肥，當(dāng)年施肥的氮磷流失率分別為2.2%和0.6%[19]。因此，對(duì)ANPSP的影響更多來(lái)自多年施肥的累積效應(yīng)而不是當(dāng)年施肥[8]。本研究的結(jié)果有助于對(duì)這種累積效應(yīng)進(jìn)行量化，計(jì)算輸出系數(shù)，為進(jìn)一步估算研究區(qū)域的氮磷面源污染負(fù)荷奠定基礎(chǔ)[15，20]。

參考文獻(xiàn)：

[1] 祁俊生. 農(nóng)業(yè)面源污染綜合防治技術(shù)[M]. 成都：西南交通大學(xué)出版社，2009.

[2] ZAKI U A，SALMAN S D，GANG D C. Nonpoint source pollution[J]. Water environmental research，2016，88（10）：1594-1619.

[3] RISSMAN A R，CARPENTER S R. Progress on nonpoint pollution： Barriers and opportunities[J]. Daedalus，2015，144（3）：35-47.

[4] 冉江華，黃?潔. 農(nóng)業(yè)面源污染研究現(xiàn)狀及發(fā)展趨勢(shì)[J]. 山西農(nóng)業(yè)科學(xué)，2009，37（3）：7-10.

[5] 呂?耀. 農(nóng)業(yè)生態(tài)系統(tǒng)中氮素造成的非點(diǎn)源污染[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，1998，17（1）：35-39.

[6] 李曉俐. 我國(guó)農(nóng)業(yè)面源污染現(xiàn)狀及綜合治理[J]. 農(nóng)業(yè)災(zāi)害研究，2015（6）：32-34.

[7] 王艷語(yǔ)，苗俊艷. 世界及我國(guó)化肥施用水平分析[J]. 磷肥與復(fù)肥，2016，31（4）：22-23.

[8] 劉欽普. 中國(guó)化肥施用強(qiáng)度及環(huán)境安全閾值時(shí)空變化[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2017，33（6）：214-221.

[9] GOULDING K W T，POUHON P R，WEBSTER C P，et al. Nitrate leaching from the Broadbalk Wheat Experiment，Rothamsted，UK，as influenced by fertilizer and manure inputs and the weather[J]. Soil use and management，2000，16： 244-250.

[10] 曹雋雋，周?勇，吳宜進(jìn)，等. 江漢平原土地利用演變對(duì)區(qū)域徑流量影響[J]. 長(zhǎng)江流域資源與環(huán)境，2013，22（5）：610-617.

[11] 張冰玢，焦利民. 快速城鎮(zhèn)化背景下江漢平原耕地變化分析[J]. 湖北農(nóng)業(yè)科學(xué)，2015，54（18）：4468-4472.

[12] 范先鵬，楊?利，張富林，等. 湖北省農(nóng)田地表徑流氮磷排放規(guī)律研究與防控技術(shù)應(yīng)用[J]. 中國(guó)科技成果，2016（8）：59-60.

[13] 鮑士旦. 土壤農(nóng)化分析[M]. 北京：中國(guó)農(nóng)業(yè)出版社，1999. 50-59.

[14] 陳效民，吳華山，孫靜紅. 太湖地區(qū)農(nóng)田土壤中銨態(tài)氮和硝態(tài)氮的時(shí)空變異[J]. 環(huán)境科學(xué)，2006，27（6）：1217-1222.

[15] 周躍龍，汪懷建，余?輝，等. 應(yīng)用輸出系數(shù)模型對(duì)太湖流域面源污染負(fù)荷測(cè)算研究[J]. 江西農(nóng)業(yè)大學(xué)學(xué)報(bào)，2014，36（3）：678-683.

[16] 張國(guó)梁，章 申. 農(nóng)田氮素淋失研究進(jìn)展[J]. 土壤，1998，30（6）：291-297.

[17] 李學(xué)平，鄒美玲. 農(nóng)田土壤磷素流失研究進(jìn)展[J]. 中國(guó)農(nóng)學(xué)通報(bào)，2010，26（11）：173-177.

[18] 馬永生，張淑英，鄧蘭萍. 氮、磷在農(nóng)田溝渠濕地系統(tǒng)中的遷移轉(zhuǎn)化機(jī)理及其模型研究進(jìn)展[J]. 甘肅科技，2005，21（2）：106-107.

[19] 劉瓊瓊，張?彥，王妍溪，等. 基于控制單元的獨(dú)流減河流域農(nóng)村非點(diǎn)源污染空間特征解析[J]. 環(huán)境污染與防治，2018，40（9）：18-21，26.

[20] 劉亞瓊，楊玉林，李法虎. 基于輸出系數(shù)模型的北京地區(qū)農(nóng)業(yè)面源污染負(fù)荷估算[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2011，27（7）：7-12.