農(nóng)田化肥氮磷地表徑流污染風(fēng)險評估

孫鋮，周華真，陳磊，沈珍瑤

(北京師范大學(xué)環(huán)境學(xué)院,北京 100875）

農(nóng)田化肥氮磷地表徑流污染風(fēng)險評估

孫鋮，周華真，陳磊*，沈珍瑤

(北京師范大學(xué)環(huán)境學(xué)院,北京 100875）

以長江中下游地區(qū)所處的六省一市為研究對象,通過對80個市級行政區(qū)化肥施用情況的調(diào)研,估算了農(nóng)田化肥氮磷地表徑流流失量,在耦合農(nóng)業(yè)化肥流失量、降雨和河網(wǎng)密度三種因素的基礎(chǔ)上,提出了農(nóng)田化肥氮磷污染風(fēng)險分級方法并初步識別了重點區(qū)域。結(jié)果表明,近20年來,其化肥施用量呈現(xiàn)上升趨勢,氮磷地表流失量較高的區(qū)域主要集中于湖南省,氮磷流失強(qiáng)度的平均值分別為2.41、0.61 kg·hm-2·a-1。湖南省為六省一市中農(nóng)田化肥氮磷污染的高風(fēng)險區(qū),上海市為低風(fēng)險區(qū),其他五省則以中、低污染風(fēng)險為主。研究結(jié)果有助于實現(xiàn)對農(nóng)田面源污染的風(fēng)險防范,推動我國農(nóng)業(yè)面源污染防控的優(yōu)化升級。

農(nóng)業(yè)面源污染；化肥；氮磷；風(fēng)險評估

保障糧食供應(yīng)是我國經(jīng)濟(jì)社會可持續(xù)發(fā)展的重要基礎(chǔ),化肥施用是糧食增產(chǎn)的重要環(huán)節(jié)［1］。但隨著我國人口政策的改變,近年來我國人口增長速度有所加快,人均耕地面積隨之減少,如何合理施用化肥以保證農(nóng)作物產(chǎn)量成為了我國亟待解決的問題［2］。根據(jù)聯(lián)合國糧食及農(nóng)業(yè)組織(Food and Agriculture Organization of the United Nations,FAO）和中國國家統(tǒng)計局提供的化肥施用量數(shù)據(jù),2013年全球化肥施用總量為1.67×108t,中國農(nóng)業(yè)化肥施用當(dāng)量為5.91×107t,占全球化肥使用總量的35.48%。同時,單位面積化肥施用量超過世界平均水平三倍多,我國農(nóng)田化肥過量施用面臨的環(huán)境風(fēng)險是其他國家無法相比的［3-4］。因此,在我國進(jìn)行化肥施用調(diào)研及時空尺度的分析有益于環(huán)境管理與污染控制。

近年來,由施肥不當(dāng)或過量施肥帶來的環(huán)境污染問題越來越突出。農(nóng)田氮磷隨地表徑流進(jìn)入受納水體,或隨地下淋溶進(jìn)入土壤,引起水體富營養(yǎng)化和土壤污染等農(nóng)業(yè)面源污染問題［5］。而地下淋溶對于南方地區(qū)影響較小,本研究主要針對地表徑流中的農(nóng)田氮磷流失進(jìn)行估算。隨著工業(yè)污染的治理,農(nóng)業(yè)氮磷流失成為水環(huán)境污染的主要來源［6］。農(nóng)田氮磷流失包括地表徑流和地下淋溶兩個主要方面,其中地表徑流為氮磷流失的主要途徑,是農(nóng)田化肥進(jìn)入周圍水體最直接、最快速的方式［7］,但由于其隨機(jī)性強(qiáng)、周期長等特點,治理起來較為困難,因此控制農(nóng)田氮磷流失是控制農(nóng)業(yè)面源污染的重要途徑。目前,我國氮肥利用率僅為30%～35%,磷肥利用率不足20%,遠(yuǎn)低于發(fā)達(dá)國家水平［8］。為此,2016年中央一號文件提出了“加大農(nóng)業(yè)面源污染防治力度,實施化肥農(nóng)藥零增長行動”的目標(biāo),對農(nóng)田化肥調(diào)研以及對地表徑流氮磷流失的估算與評估有助于實現(xiàn)對農(nóng)田面源污染的風(fēng)險防范。

本文以長江中下游所處的六省一市為研究區(qū),以市為單位,開展農(nóng)業(yè)氮磷流失(地表徑流流失）估算方法和風(fēng)險區(qū)域識別研究。該區(qū)域是我國重要的農(nóng)業(yè)生產(chǎn)基地,但近年來亦產(chǎn)生較多的農(nóng)業(yè)面源問題［9］。由于化肥的不合理利用,導(dǎo)致富營養(yǎng)問題頻發(fā),水環(huán)境惡化,水質(zhì)降低［10］。長江中下游區(qū)域遼闊,因此本文在估算流失量的同時,綜合沖刷、污染和河道三個過程進(jìn)行合理地污染風(fēng)險評估,以對研究區(qū)的地表氮磷污染進(jìn)行重點區(qū)域識別。

1 材料與方法

1.1 試驗區(qū)概況

研究區(qū)位于我國東南部,地跨湖南、湖北、江西、安徽、江蘇、浙江、上海六省一市(圖1）。該區(qū)域大部分處于長江中下游平原,總面積約為9×106km2,海拔約為5～100 m,年均溫度約為14～18℃,年降水量約為1000～1500 mm。研究區(qū)屬于亞熱帶季風(fēng)氣候,以沖積平原為主,具有較為優(yōu)越的農(nóng)業(yè)氣候條件［11］。同時長江中下游地區(qū)是我國重要的糧、油、棉生產(chǎn)基地,亦為中國水資源最為豐富的地區(qū),主要河流有長江、漢江、贛江等［12］。大部分城市也是我國重要的工農(nóng)業(yè)基地,是經(jīng)濟(jì)和科技文化發(fā)達(dá)地區(qū),其中湖南、湖北、江蘇省為我國的農(nóng)業(yè)大省,農(nóng)業(yè)生產(chǎn)活動較多［13］。雙季稻是該區(qū)域的主要種植模式,水田亦是造成化肥流失的主要區(qū)域。

圖1 研究區(qū)地理位置Figure 1 The location of study area

1.2 數(shù)據(jù)資料

本研究的基礎(chǔ)數(shù)據(jù)資料主要包括市級行政單元劃分、農(nóng)作物播種面積、市級化肥施用數(shù)據(jù)、省級化肥施用數(shù)據(jù)、流失系數(shù)以及確定流失系數(shù)所需的相關(guān)數(shù)據(jù)。主要數(shù)據(jù)及來源見表1。

1.3 方法介紹

1.3.1 化肥施用量及施用強(qiáng)度的估算方法

本研究從中國統(tǒng)計年鑒中獲取了2000—2015年六省一市農(nóng)田化肥的施用數(shù)據(jù)(市級）。有些年份或地區(qū)數(shù)據(jù)存在缺失,為了提高精確性和可利用性,選取了數(shù)據(jù)相對完善、年份較近的2013年數(shù)據(jù)進(jìn)行化肥施用的空間分布特征分析。同時,對缺失數(shù)據(jù)進(jìn)行估算,估算方法分為以下兩種.

一是只存在化肥施用總量,無氮肥、磷肥和復(fù)合肥施用數(shù)據(jù)的市級行政區(qū),采取以下公式進(jìn)行估算.

二是部分2013年統(tǒng)計數(shù)據(jù)缺失的市級行政區(qū),取相鄰年份(2012、2014）的化肥施用數(shù)據(jù)的平均值進(jìn)行替代。

表1 應(yīng)用數(shù)據(jù)及來源Table 1 Data and sources

農(nóng)田氮磷肥的總施用量包括單一元素肥料和復(fù)合肥中氮/磷有效成分之和。復(fù)合肥中氮磷鉀有效成分采用通用型化肥比例15∶15∶15進(jìn)行計算［20］,即.

農(nóng)田總施氮/磷量(×105t）=氮/磷肥施用量+復(fù)合肥施用量×15%(2）化肥施用強(qiáng)度表示為單位市級行政區(qū)面積的年化肥施用量(kg·hm-2·a-1）,相對于化肥施用總量更能直觀地反映一個地區(qū)的化肥施用狀況。計算公式為.

1.3.2 地表流失量及流失率的計算

本研究以農(nóng)業(yè)部出版的《全國農(nóng)業(yè)面源污染流失系數(shù)手冊》中的農(nóng)田地表徑流流失系數(shù)為基礎(chǔ),對研究區(qū)80個市級行政區(qū)的農(nóng)田氮磷流失系數(shù)進(jìn)行查算。根據(jù)5個要素對各田塊的流失系數(shù)進(jìn)行劃分,即“農(nóng)區(qū)-地形-梯田/非梯田-旱地/水田-種植類型”(例如.模式47為“南方山地丘陵區(qū)-陡坡地-梯田-旱地-露地蔬菜”）,獲得的氮磷流失系數(shù)分別為1.5%、0.85%。根據(jù)中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心獲取到的數(shù)據(jù),對研究區(qū)各種要素進(jìn)行劃分,從而查詢《全國農(nóng)業(yè)面源污染流失系數(shù)手冊》獲得各個模式的地表氮磷流失系數(shù)。并利用公式(4）將各個模式升尺度,獲得各市級行政區(qū)的氮磷流失系數(shù)μ.

式中：βi為某一田塊下某農(nóng)作物種植面積占市級行政區(qū)農(nóng)作物總播種總面積的比例,%；μi為該田塊的氮磷流失系數(shù)。

其中《全國農(nóng)業(yè)面源污染流失系數(shù)手冊》定義流失系數(shù)公式為.

式中：RL為常規(guī)處理下,即正常施肥狀態(tài)下氮/磷肥流失量,×105t；CK為對照處理氮/磷肥的流失量,×105t；FA為氮/磷肥的施用量,×105t。因此,氮磷肥的流失量可由流失系數(shù)與氮磷肥施用量相乘獲得。

1.3.3 面源污染風(fēng)險評價方法的構(gòu)建

本文按照氮磷指數(shù)法的思路,構(gòu)建農(nóng)田面源污染風(fēng)險方法評價,該方法是綜合考慮了影響農(nóng)業(yè)面源污染的主要因子,評價流域內(nèi)不同地區(qū)發(fā)生氮磷流失危險性高低的一種方法［21］。該方法考慮三種影響氮磷流失過程的因子,即主要考慮沖刷過程(坡面遷移能力）、污染過程(污染物多少）和河道過程(造成水污染的可能性）。其中,沖刷過程主要考慮降雨因子,污染過程主要考慮氮磷流失量,河道過程主要考慮河網(wǎng)密度。降雨數(shù)據(jù)以研究區(qū)的年平均降雨量數(shù)據(jù)為準(zhǔn),氮磷流失量以本文計算為準(zhǔn),而河網(wǎng)密度則以水系數(shù)據(jù)為基礎(chǔ)經(jīng)簡化計算獲得(將河網(wǎng)長度除以行政區(qū)面積）。

面源農(nóng)田污染風(fēng)險等級構(gòu)建如下.首先,市級行政區(qū)的農(nóng)田氮磷流失量、河網(wǎng)密度、降雨等因素的數(shù)值進(jìn)行歸一化處理,見公式(6）和公式(7）；其次,將各因子歸一化的值相乘,計算氮磷潛在危險性指數(shù),為了簡化計算過程,本研究將各影響因子賦予相同的權(quán)重；最后,由高到低進(jìn)行排序,并將各行政區(qū)的風(fēng)險值分為若干等級,在此基礎(chǔ)上識別面源污染風(fēng)險的重點區(qū)域。

式中：xi為各個因子的實際值；xmin為80個市級行政區(qū)中的最小值；xmax為80個市級行政區(qū)中的最大值；x*i為各個因子歸一化處理后所得的值。最終得到相應(yīng)的污染分級方案如表2所示。

2 結(jié)果與討論

2.1 化肥施用的時空分布特征分析

2.1.1 化肥施用的時間變化特征

長江中下游地區(qū)六省一市施肥總量的時間變化趨勢如圖2所示。自1996年至2015年20年內(nèi),施肥總量波動不大,總體呈現(xiàn)上升趨勢,施肥總量范圍在1.233×108～1.513×108t之間,約占全國化肥施用總量的30%。由此說明,研究區(qū)所處的六省一市一直為我國的農(nóng)業(yè)中心。自2011年起,施肥總量達(dá)到1.5×108t,2015年再次下降至1.4×108t,這與2014年我國提出的適時開征化肥稅引導(dǎo)農(nóng)民減少化肥的過度施用有關(guān)［22］。其中,對于各類化肥施用來說,氮、磷、鉀肥施用量基本保持不變,波動較小,氮、磷、鉀、復(fù)合肥所占比例的均值分別為47.97%、16.07%、9.67%、26.29%?？梢钥闯?氮肥施用量最高,鉀肥施用量最低,而復(fù)合肥施用呈現(xiàn)明顯上升趨勢。在農(nóng)田化肥所含有的三種元素中,氮、磷仍是農(nóng)田化肥所含的主要元素,也是造成農(nóng)業(yè)面源污染的潛在關(guān)鍵元素,因此氮磷流失是本研究所要關(guān)注的重點內(nèi)容［23］。同時,糧食產(chǎn)量自2003年至2012年呈現(xiàn)上升趨勢,至2012年糧食產(chǎn)量已達(dá)1.508×109t,而當(dāng)年我國糧食總產(chǎn)量為5.896×109t,實現(xiàn)糧食產(chǎn)量“九連增”,而研究區(qū)糧食總產(chǎn)量約為全國的25.58%［24］。同樣的,2012年至2015年研究區(qū)的糧食產(chǎn)量也呈現(xiàn)上升趨勢,而1996年至2003年糧食產(chǎn)量呈現(xiàn)下降趨勢。由此可以看出,化肥施用量的增加在一定程度上促進(jìn)糧食增產(chǎn),但化肥施用和糧食產(chǎn)量沒有呈現(xiàn)完全的正相關(guān),合理的化肥施用更能保障糧食產(chǎn)量和環(huán)境管理。

2.1.2 化肥施用的空間分布特征

長江中下游地區(qū)幅員遼闊,各個省市的社會經(jīng)濟(jì)和自然環(huán)境狀況存在明顯差異。因此,本研究利用地理信息系統(tǒng)進(jìn)行空間分布特征分析,如圖3所示。為了更直觀地顯示市級行政區(qū)的氮磷施用狀況,采用化肥施用強(qiáng)度作為評估指標(biāo)。對于氮肥施用強(qiáng)度來說,湖南省是施用強(qiáng)度最大的省份,東部地區(qū)都為高氮肥施用強(qiáng)度區(qū),西部地區(qū)處于較高氮肥施用強(qiáng)度地區(qū)；相比之下,江西省和上海市的大部分區(qū)域氮肥施用強(qiáng)度都處于較低水平；磷肥施用強(qiáng)度與氮肥施用強(qiáng)度呈現(xiàn)出明顯差異。湖南省東部和湖北省的中部屬于磷肥流失強(qiáng)度較高的區(qū)域,對這幾個市級行政區(qū)依次排序為湘潭市＞益陽市＞常德市＞宜昌市＞岳陽市＞襄樊市；而浙江省、江西省和上海市的大部分地區(qū)磷肥施用強(qiáng)度較低。綜上所述,湖南省為六省一市中氮肥施用強(qiáng)度最大的省份,而磷肥施用強(qiáng)度較高的為湖南湖北兩省的部分區(qū)域。因此,對于肥料施用控制管理,應(yīng)具體位置具體分析,同一個省份不同市級單位的化肥施用量仍具有一定的差異。同時,上海市的氮磷施用強(qiáng)度一直處于較低水平,這與上海城市化水平較高、農(nóng)田面積比例較小有關(guān)。

圖2 1996年至2015年的化肥施用狀況Figure 2 The variation of fertilizers application from 1996 to 2015

表2 農(nóng)田氮磷污染風(fēng)險分級方案Table 2 The risk classification scheme for farmland nitrogen and phosphrus

2.2 農(nóng)田化肥氮磷流失特征

圖3 2013年農(nóng)田氮磷施用強(qiáng)度空間分布Figure 3 The spatial distribution of nitrogen and phosphorus fertilizer intensity

不同的流失率將直接導(dǎo)致氮磷流失量的不同,農(nóng)田化肥氮磷流失量的空間分布如圖4所示?？傮w來說,湖南省是氮磷流失的主要區(qū)域。對于氮肥流失量,湖南省處于高流失量的市級行政區(qū)排序為.益陽市(7.66×103t）＞岳陽市(7.31×103t）＞衡陽市(7.16×103t）＞永州市(6.05×103t）＞長沙市(5.37×103t）＞邵陽市(5.26×103t）,其次湖南省的郴州市、株洲市、湘潭市以及湖北省的襄樊市和江蘇省的徐州市、鹽城市屬于次高氮流失區(qū)；相比之下,江西省和浙江省的大部分區(qū)域和上海市都屬于氮肥流失較少的區(qū)域?？傮w上來看,研究區(qū)氮肥流失量處于0.521×103～1.214×103t范圍的居多。對于磷肥流失量,湖南省具有較高流失量的市級行政區(qū)排序為.常德市(2.69×103t）＞益陽市(1.97×103t）＞岳陽市(1.81×103t）＞衡陽市(1.73×103t）,湖南省的其他市級行政區(qū)磷流失量也屬于較高水平；同時,其他五省一市的磷流失量基本小于1000 t,尤其是浙江省的全部區(qū)域磷的流失量小于200 t?？傮w上來看,研究區(qū)磷肥流失量處于0.154×103～0.381× 103t范圍的居多。

圖4 2013年農(nóng)田氮磷地表流失量的空間分布Figure 4 The spatial distribution of the amount of nitrogen and phosphorus surface loss

氮磷流失強(qiáng)度分布如圖5所示。對于氮肥施用強(qiáng)度來說,研究區(qū)流失強(qiáng)度的平均值為2.41 kg·hm-2·a-1；而湖南省大部分地區(qū)都屬于氮肥流失強(qiáng)度較高地區(qū),氮肥施用強(qiáng)度遠(yuǎn)高于平均水平,對這幾個市級行政區(qū)依次排序為常德市＞岳陽市＞長沙市＞株洲市＞衡陽市＞永州市＞郴州市＞婁底市＞邵陽市,湖北省的襄樊市和宜昌市、江蘇省的北部地區(qū)氮肥流失強(qiáng)度也處于較高水平；而江西省大部分地區(qū)氮肥流失強(qiáng)度都處于較低水平。對于磷肥施用強(qiáng)度,研究區(qū)流失強(qiáng)度的平均值為0.61 kg·hm-2·a-1。與氮肥流失強(qiáng)度相同,湖南省大部分地區(qū)都屬于磷肥流失強(qiáng)度較高地區(qū),其中湖南省的東北部屬于流失強(qiáng)度最高的地區(qū),對這幾個市級行政區(qū)依次排序為湘潭市＞益陽市＞常德市＞岳陽市＞長沙市＞株洲市；而安徽省、江西省和江蘇省的大部分地區(qū)磷肥流失強(qiáng)度較低。綜上所述,湖南省為氮肥、磷肥流失強(qiáng)度最大的一個省份,同時氮肥、磷肥的流失強(qiáng)度在空間分布上稍有不同。這與不同地區(qū)的種植模式有關(guān),也與氮磷元素在自然界中的形態(tài)有關(guān),氮以溶解態(tài)為主,磷以吸附態(tài)為主,氮肥與磷肥相比更易流失［25］。

圖5 2013年農(nóng)田地表氮磷流失強(qiáng)度的空間分布Figure 5 The spatial distribution of the nitrogen and phosphorus surface loss intensity

2.3 農(nóng)田面源污染風(fēng)險評估

利用1.3.3的方法對長江中下游地區(qū)六省一市開展農(nóng)田面源污染重點區(qū)域識別,風(fēng)險分級結(jié)果如圖6所示?？傮w上來看,大部分高流失強(qiáng)度區(qū)域亦是高污染風(fēng)險區(qū),但有些區(qū)域,例如江蘇省的蘇州市,磷肥流失強(qiáng)度為0.351 kg·hm-2·a-1,卻處于較高流失風(fēng)險區(qū),這與蘇州市臨近長江,河網(wǎng)密度較大有關(guān)。同時,氮磷地表流失污染存在著一定的空間差異,但湖南省中東部地區(qū)同時處于氮磷流失的高風(fēng)險區(qū),上海市和安徽省的北部地區(qū)同時處于氮磷流失的低風(fēng)險區(qū)。對于農(nóng)田氮污染風(fēng)險分區(qū),湖南省的11個市級行政單元、安徽省的銅陵市和蕪湖市、江蘇省的淮安市和宿遷市、湖北省的鄂州市屬于高污染風(fēng)險區(qū)；湖南省的其他地區(qū)均屬于較高污染風(fēng)險區(qū),其他五省也均有分布；相比之下,上海市、安徽省北部和江西省的大部分地區(qū)都處于氮的低污染風(fēng)險區(qū)；同時,高污染風(fēng)險區(qū)約占研究區(qū)面積的20.27%,湖南省高污染風(fēng)險區(qū)約占總面積的75.19%。對于磷地表徑流污染風(fēng)險分區(qū),湖南省的11個市級行政單元、安徽省的安慶市和銅陵市、湖北省的鄂州市、江西省的南昌市屬于高污染風(fēng)險區(qū)；湖南省的其他地區(qū)均屬于較高污染風(fēng)險區(qū),除浙江省外其他四省也均有分布；相比之下,上海市、安徽省北部和南部、江西省東部部分區(qū)域處于氮的低污染風(fēng)險區(qū)；同時,高污染風(fēng)險區(qū)約占研究區(qū)面積的20.07%,湖南省高污染風(fēng)險區(qū)與氮肥高污染風(fēng)險區(qū)相同,約占總面積的75.19%。綜上所述,湖南省為氮磷流失的高風(fēng)險區(qū)域,上海市為氮磷流失的低風(fēng)險區(qū)域,其他五省以中、低污染風(fēng)險為主。因此,長江中下游地區(qū)所處的六省一市污染防治宜“以防為主,以治為輔”,但對于湖南省的高污染風(fēng)險區(qū)域應(yīng)加強(qiáng)控制措施［12］。

圖6 2013年農(nóng)田氮磷地表流失污染風(fēng)險分級Figure 6 The risk level of surface runoff loss of nitrogen and phosphorus fertilizer in 2013

3 結(jié)論

本研究以長江中下游所處的六省一市為研究區(qū),調(diào)研了區(qū)域內(nèi)市級行政區(qū)化肥數(shù)據(jù)及農(nóng)作物播種面積,以農(nóng)田氮磷流失系數(shù)為基礎(chǔ)計算出對應(yīng)的流失量,并利用地理信息系統(tǒng)進(jìn)行了降雨量和河網(wǎng)密度的空間分析,確定了農(nóng)田污染風(fēng)險等級。該方法可直觀地比較各市級行政區(qū)的氮磷流失量,識別農(nóng)田污染高風(fēng)險區(qū),為農(nóng)田面源防控和水環(huán)境治理提供理論基礎(chǔ)。

案例研究結(jié)果表明,近20年來,研究區(qū)化肥施用量呈現(xiàn)上升趨勢。湖南省是氮肥施用強(qiáng)度和流失強(qiáng)度最大的省份。湖南、湖北則為磷肥施用強(qiáng)度較高的區(qū)域,該區(qū)域氮磷流失強(qiáng)度的平均值分別為2.41、0.61 kg·hm-2·a-1。從風(fēng)險分級的角度,農(nóng)業(yè)化肥對周邊水體的高污染風(fēng)險區(qū)為湖南省,上海市為低風(fēng)險區(qū)域,其他五省則以中、低污染風(fēng)險為主。

［1］朱兆良,金繼運.保障我國糧食安全的肥料問題［J］.植物營養(yǎng)與肥料學(xué)報,2013,19(2）.259-273.

ZHU Zhao-liang,JIN Ji-yun.Fertilizer use and food security in China［J］.Plant Nutrition and Fertilizer Science,2013,19(2）.259-273.

［2］張嘉云,李志鑫,喬志遠(yuǎn),等.我國農(nóng)田化肥施用現(xiàn)狀［J］.北方農(nóng)業(yè)學(xué)報,2016,44(3）.118-119.

ZHANG Jia-yun,LI Zhi-xin,QIAO Zhi-yuan,et al.The application status of farmland fertilizer in China［J］.Journal of Northern Agriculture, 2016,44(3）.118-119.

［3］劉欽普.江蘇氮磷鉀化肥使用地域分異及環(huán)境風(fēng)險評價［J］.應(yīng)用生態(tài)學(xué)報,2015,26(5）.1477-1483.

LIU Qin-pu.Regional difference of NPK fertilizers application and environmental risk assessment in Jiangsu Province,China［J］.Chinese Journal of Applied Ecology,2015,26(5）.1477-1483.

［4］向平安,周燕,黃璜,等.氮肥面源污染控制的綠稅激勵措施探討——以洞庭湖區(qū)為例［J］.中國農(nóng)業(yè)科學(xué),2007,40(2）.330-337.

XIANG Ping-an,ZHOU Yan,HUANG Huang,et al.Discussion on the green tax stimulation measure of nitrogen fertilizer non-point source pollution control.Using the Dongting Lake area［J］.Scientia Agriculture Sinica,2007,40(2）.330-337.

［5］Wang A,Tang L H,Yang D W,et al.Spatio-temporal variation of net anthropogenic nitrogen inputs in the upper Yangtze River basin from 1990 to 2012［J］.Science ChinaEarth Sciences,2016,59(11）.1-13.

［6］張俊卿.密云水庫流域森林植被對地表徑流非點源污染的影響［D］.北京.北京林業(yè)大學(xué),2008.

ZHANG Jun-qing.Effect of forest vegetation on the non-point source pollution of surface runoff in Miyun reservoir watershed［D］.Beijing.Beijing Forestry University,2008.

［7］范麗麗,沈珍瑤,劉瑞民.不同降雨-徑流過程中農(nóng)業(yè)非點源污染研究［J］.環(huán)境科學(xué)與技術(shù),2008,31(10）.5-8.

FAN Li-li,SHEN Zhen-yao,LIU Rui-min.Agricultural non-point source pollution in different rainfall runoff［J］.Environmental Science &Technology,2008,31(10）.5-8.

［8］楊憲龍,同延安,路永莉,等.農(nóng)田氮肥利用率計算方法研究進(jìn)展［J］.應(yīng)用生態(tài)學(xué)報,2015,26(7）.2203-2212.

YANG Xian-long,TONG Yan-an,LU Yong-li,et al.Research advances in the calculating method of nitrogen use efficiency(NUE）in cultivated lands［J］.Chinese Journal of Applied Ecology,2015,26(7）. 2203-2212.

［9］劉松,朱建強(qiáng),田皓.長江中下游地區(qū)的主要水問題與對策［J］.長江大學(xué)學(xué)報(自然科學(xué)版）,2012,9(1）.42-46.

LIU Song,ZHU Jian-qiang,TIAN Hao.The main water problem and measuresinthemiddleandlowoftheYangtzeRiver［J］.Journal ofYangtze University(Nat Sci Edit）,2012,9(1）.42-46.

［10］姚瑞華,趙越,王東,等.長江中下游流域水環(huán)境現(xiàn)狀及污染防治對策［J］.人民長江,2014(增刊1）.45-47.

YAO Rui-hua,ZHAO Yue,WANG Dong,et al.Water environment situation and pollution measurements in the middle and low of the Yangtze River basin［J］.Yangtze River,2014(Suppl 1）.45-47.

［11］賴敏,王偉力,郭靈輝.長江中下游城市群農(nóng)業(yè)面源污染氮排放評價及調(diào)控［J］.中國農(nóng)業(yè)資源與區(qū)劃,2016,37(8）.1-11.

LAI Min,WANG Wei-li,GUO Ling-hui.Nitrogen emissions evaluation and control of agricultural non-point source pollution in city group in the middle and low of the Yangtze River［J］.Chinese Journal of Agricultural Resources and Regional Planning,2016,37(8）.1-11.

［12］何書樵,鄭有飛,尹繼福.近50年長江中下游地區(qū)降水特征分析［J］.生態(tài)環(huán)境學(xué)報,2013,22(7）.1187-1192.

HE Shu-qiao,ZHENG You-fei,YIN Ji-fu.An analysis on precipitation characteristics over middle and lower reaches of Yangtze River in the last 50 years［J］.Ecology and Environmental Sciences,2013,22(7）. 1187-1192.

［13］任天志,劉宏斌,范先鵬,等.全國農(nóng)田面源污染流失系數(shù)手冊［M］.北京.中國農(nóng)業(yè)出版社,2015.8-77.

REN Tian-zhi,LIU Hong-bin,FAN Xian-peng,et al.Emission coefficient manual of agricultural non-point source pollution in China［M］. Beijing.China Agriculture Press,2015.8-77.

［14］中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心.地形地貌數(shù)據(jù)庫［DB/OL］.［2017-06-01］.http://www.resdc.cn.

Data Center for Resources and Environmental Sciences,Chinse Academy of Science(RESDC）.Topography database［DB/OL］.［2017-06-01］. http://www.resdc.cn.

［15］中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心.自然地理分區(qū)數(shù)據(jù)庫［DB/ OL］.［2017-06-01］.http://www.resdc.cn.

Data Center for Resources and Environmental Sciences,Chinse Academy of Science(RESDC）.Natural geographical partition database［DB/ OL］.［2017-06-01］.http://www.resdc.cn.

［16］中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心.土地利用數(shù)據(jù)庫［DB/OL］.［2017-06-01］.http://www.resdc.cn.

Data Center for Resources and Environmental Sciences,Chinse Academy of Science(RESDC）.Land use database［DB/OL］.［2017-06-01］. http://www.resdc.cn.

［17］中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心.中國氣象要素空間插值數(shù)據(jù)庫［DB/OL］.［2017-06-01］.http://www.resdc.cn.

Data Center for Resources and Environmental Sciences,Chinse Academy of Science(RESDC）.China meteorological elements spatial interpolation database［DB/OL］.［2017-06-01］.http://www.resdc.cn.

［18］中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心.全球100萬基礎(chǔ)地理數(shù)據(jù)庫［DB/OL］.［2017-06-01］.http://www.resdc.cn.

Data Center for Resources and Environmental Sciences,Chinse Academy of Science(RESDC）.Global 1∶1000000 Fundamental geographic database［DB/OL］.［2017-06-01］.http://www.resdc.cn.

［19］張衛(wèi)峰,李亮科,陳新平,等.我國復(fù)合肥發(fā)展現(xiàn)狀及存在的問題［J］.磷肥與復(fù)肥,2009,24(2）.14-16.

ZHANG Wei-feng,LI Liang-ke,CHEN Xin-ping,et al.The present status and existing problems in China′s compound fertilizer development［J］.Phosphate&Compound Fertilizer,2009,24(2）.14-16.

［20］張平,高陽昕,劉云慧,等.基于氮磷指數(shù)的小流域氮磷污染風(fēng)險評價［J］.生態(tài)環(huán)境學(xué)報,2011,20(增刊1）.1018-1025.

ZHANG Ping,GAO Yang-xin,LIU Yun-hui,et al.Nitrogen and phosphorus pollution risk assessment in small catchment based on the nitrogen and phosphorus index［J］.Ecology and Environmental Sciences, 2011,20(Suppl 1）.1018-1025.

［21］熊冬洋.對我國適時開征化肥稅的思考［J］.稅務(wù)與經(jīng)濟(jì),2014(2）. 97-100.

XIONG Dong-yang.Reflection on timely taxation in China［J］.Taxation and Economy,2014(2）.97-100.

［22］程波,張澤,陳凌,等.太湖水體富營養(yǎng)化與流域農(nóng)業(yè)面源污染的控制［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2005,24(增刊1）.118-124.

CHENG Bo,ZHANG Ze,CHEN Ling,et al.Eutrophication of Taihu Lake and pollution from agricultural non-point sources in Lake Taihu Basin［J］.Journal of Agro-Environment Science,2005,24(Suppl 1）. 118-124.

［23］胡傳景.糧食產(chǎn)量九連增耕地保護(hù)不能松［J］.浙江國土資源,2013 (3）.43-46.

HU Chuan-jing.The increasion of food production in nine years,strict enforcement of the cultivated land protection［J］.Land and Resources of ZhejiangProvince,2013(3）.43-46.

［24］Otto R,de Freitas Júnior J C M,Zavaschi E,et al.Combined application of concentrated vinasse and nitrogen fertilizers in sugarcane. Strategies to reduce ammonia volatilization losses［J］.Sugar Tech,2017, 19(3）.248-257.

The pollution risk assessment of nitrogen and phosphorus loss in surface runoff from farmland fertilizer

SUN Cheng,ZHOU Hua-zhen,CHEN Lei*,SHEN Zhen-yao
(Shool of Environment,Beijing Normal University,Beijing 100875,China）

Nitrogen and phosphorus loss in farmland fertilizer is the main cause of environmental water pollution.Using six provinces and one city in the middle and lower reaches of the Yangtze River as a study area,and by investigating on the data of fertilizer of 80 municipals, we estimated the loss of nitrogen and phosphorus in surface runoff from farmland fertilizer.Based on the combination of three factors,namely fertilizer loss,rainfall,and drainage density,we put forward a scheme to classify the pollution risk from fertilizer nitrogen and phosphorus, and identified the key areas of concern preliminarily.The results show that fertilizer application has increased during the past 20 years,and areas with high nitrogen loss and phosphorus surface runoff are concentrated in Hunan Province.The average nitrogen and phosphorus loss intensities were 2.41 and 0.61 kg·hm-2·a-1,respectively.The zone with a high risk of pollution was Hunan Province,Shanghai City was the low risk zone,and the other five provinces were dominantly middle-and low-risk areas.These results are beneficial to preventing risk from agricultural non-point source pollution and promoting the optimization of control measures in China.

agricultural non-point source pollution;fertilizer;nitrogen and phosphorus;environmental risk assessment

X820.4

A

1672-2043(2017）07-1266-08

10.11654/jaes.2017-0800

孫鋮,周華真,陳磊,等.農(nóng)田化肥氮磷地表徑流污染風(fēng)險評估［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2017,36(7）.1266-1273.

SUN Cheng,ZHOU Hua-zhen,CHEN Lei,et al.The pollution risk assessment of nitrogen and phosphorus loss in surface runoff from farmland fertilizer［J］.Journal of Agro-Environment Science,2017,36(7）.1266-1273.

2017-06-06

孫鋮(1994—）,女,吉林延邊人,碩士研究生,從事農(nóng)業(yè)面源污染研究工作。E-mail:13756892980@163.com

*通信作者：陳磊E-mail:chenlei1982bnu@bnu.edu.cn

國家自然科學(xué)基金項目(51579011,51409003）

Project supported：The National Natural Science Foundation of China(51579011,51409003）