內(nèi)蒙古河套灌區(qū)農(nóng)田非點源氮磷污染負荷估算

王 艷,焦 燕,楊文柱,閆穎超,鄔 宏,靈 靈,史月超

內(nèi)蒙古河套灌區(qū)農(nóng)田非點源氮磷污染負荷估算

王 艷,焦 燕*,楊文柱,閆穎超,鄔 宏,靈 靈,史月超

(內(nèi)蒙古師范大學(xué)化學(xué)與環(huán)境科學(xué)學(xué)院,內(nèi)蒙古自治區(qū)環(huán)境化學(xué)重點實驗室,內(nèi)蒙古 呼和浩特 010022)

為了解內(nèi)蒙古河套灌區(qū)農(nóng)田非點源污染現(xiàn)狀,引入降雨影響因子()和地形影響因子()對傳統(tǒng)的輸出系數(shù)模型(ECM)進行改進,量化該地區(qū)的農(nóng)田非點源氮、磷污染負荷,結(jié)合污染物入河系數(shù)(),對農(nóng)田非點源氮磷入河(烏梁素海)量進行估算,并識別污染物關(guān)鍵源區(qū).結(jié)果表明,河套灌區(qū)2021年農(nóng)田非點源總氮(TN)、總磷(TP)污染負荷為3520.998和407.125t,分別占不同土地利用類型污染負荷總量的79.674%和89.201%;農(nóng)田非點源TN、TP污染負荷空間分布一致,表現(xiàn)為靠近水體、排干的鄉(xiāng)鎮(zhèn)貢獻量最高,依次為新華鎮(zhèn)、新安鎮(zhèn)和塔爾湖鎮(zhèn),TN、TP污染負荷分別為316.658, 236.725, 199.344t/a和36.615, 27.371, 23.050t/a;農(nóng)田非點源TN、TP入河量為387.310t/a和16.285t/a,其中新華鎮(zhèn)、新安鎮(zhèn)、塔爾湖鎮(zhèn)、隆興昌鎮(zhèn)、團結(jié)鎮(zhèn)和蠻會鎮(zhèn)是河套灌區(qū)的農(nóng)田非點源污染重點防治鄉(xiāng)鎮(zhèn);與實測值進行比較,驗證了改進的輸出系數(shù)模型的合理性,相對誤差分別為-4.743%和19.037%.

河套灌區(qū)；農(nóng)田非點源污染；改進的輸出系數(shù)模型；TN；TP

在過去的幾十年里,我國氮肥和磷肥的施用在提高作物產(chǎn)量的同時,一系列農(nóng)業(yè)環(huán)境污染問題也逐漸顯現(xiàn),其中水環(huán)境污染問題尤其嚴(yán)重,水環(huán)境污染主要由點源污染和非點源污染兩部分構(gòu)成[1].隨著點源污染得到有效治理后,非點源污染對水環(huán)境的影響逐漸突出,其中,農(nóng)業(yè)非點源污染已成為水體富營養(yǎng)化的主要污染源[2].《第二次全國污染源普查公報》[3]結(jié)果顯示,2017年農(nóng)業(yè)源水污染物中的化學(xué)需氧量、氨氮、總氮(TN)以及總磷(TP)排放量分別為1067.13,21.62,141.49和21.20萬t.

然而,農(nóng)業(yè)非點源污染的廣泛性、隨機性、潛伏性[4]等特點,增加了研究和控制農(nóng)業(yè)非點源污染的難度,利用模型對農(nóng)業(yè)非點源污染負荷進行定量估算是主要的研究方法之一.輸出系數(shù)模型[5](ECM)由于不考慮污染物運移的復(fù)雜過程,具有參數(shù)少、建模成本低等特點而得到廣泛應(yīng)用,特別是對觀測數(shù)據(jù)和以往研究不足地區(qū)優(yōu)勢愈加明顯[6-7].一些學(xué)者通過對輸出系數(shù)模型進行不同方面的改進[8-10],探究更適合本流域的輸出系數(shù)模型,使模型模擬值與實測值更加吻合,并結(jié)合GIS等技術(shù)對污染物的空間分布特征進行探究.由于污染物被農(nóng)作物、河道的截留等作用而產(chǎn)生滯留[11],并不會完全進入水體.王萌等[12]在分析流域農(nóng)業(yè)非點源污染負荷時,又考慮了污染物的入河系數(shù),減少了非點源污染估算過程中的不確定性.盡管眾多學(xué)者對模型的改進提高了模型估算精度,但針對南方區(qū)域研究較多,對于北方區(qū)域研究較少.另外,將輸出系數(shù)模型應(yīng)用于農(nóng)業(yè)非點源污染負荷估算的研究多集中于污染物輸出或者單一階段的污染物入河估算[13].由于農(nóng)業(yè)非點源對水環(huán)境的污染要經(jīng)歷一定階段[14],故本研究結(jié)合污染物入河系數(shù),綜合污染物輸出量與污染物入河量估算農(nóng)業(yè)非點源污染負荷.

內(nèi)蒙古河套灌區(qū)是中國三大灌區(qū)之一,也是國家和自治區(qū)重要的商品糧油生產(chǎn)基地[15],該地區(qū)農(nóng)業(yè)活動顯著,化肥施用量大,但利用率低[16],其余部分則是隨著農(nóng)田退水和地表徑流等匯入烏梁素海湖泊,容易造成烏梁素海水環(huán)境的惡化.因此,本研究以內(nèi)蒙古河套灌區(qū)為對象,采用考慮降雨和地形影響因子的輸出系數(shù)模型,量化河套灌區(qū)農(nóng)田非點源TN、TP污染負荷,結(jié)合污染物入河系數(shù),進行農(nóng)田氮、磷污染入河量的估算,并識別污染關(guān)鍵源區(qū),為烏梁素海富營養(yǎng)化的防治以及區(qū)域水環(huán)境的保護提供參考.

1 材料與方法

1.1 研究區(qū)概況

河套灌區(qū)位于黃河內(nèi)蒙古段北岸的“幾”字彎上(40°19′～41°18′N,106°20′～109°19′E),南北跨度約50km,東西跨度約250km[17].該區(qū)域?qū)儆诘湫偷拇箨懶詺夂?降雨少、蒸發(fā)強,多年平均氣溫8.6～12℃,多年平均蒸發(fā)量2000～2400mm,多年平均降雨量111～236mm,且降水量分布極不均勻,夏季降水量(6～8月)占全年降水量的63%～70%,僅靠降雨不能滿足河套灌區(qū)農(nóng)作物的用水需求,春季播種(春灌)與夏季生長期(夏灌)都需要依靠灌溉來保證農(nóng)業(yè)生產(chǎn)[18].河套灌區(qū)年均引黃水量約48億m3[19],90%以上的農(nóng)田退水會通過排干匯入烏梁素海,從而引起烏梁素海水體的污染.研究區(qū)具體位置以及農(nóng)田退水排干分布情況見圖1.

圖1 研究區(qū)概況和土地利用類型

1.2 數(shù)據(jù)來源及處理

本研究以2021年為基準(zhǔn)年,所涉及的數(shù)據(jù)包括數(shù)字高程圖(DEM)、土地利用數(shù)據(jù)、氣象數(shù)據(jù)、水文水質(zhì)數(shù)據(jù)等,各類型空間數(shù)據(jù)的地理坐標(biāo)和投影統(tǒng)一,所有的空間操作基于ArcGIS10.7軟件進行,詳細數(shù)據(jù)來源見表1.其中,土地利用數(shù)據(jù)利用ArcGIS進行重分類,分為耕地、林地、草地、水域、建設(shè)用地和未利用地6類.

表1 數(shù)據(jù)來源

1.3 改進的輸出系數(shù)模型

輸出系數(shù)模型是一種基于統(tǒng)計學(xué)的線性數(shù)學(xué)模型,由Johnes和O’Sulivan于1989年正式提出, 隨后Johnes又對其進行了修正,一般表達式為:

式中:為土地利用類型;為污染物類型;為土地利用類型的種類;L為污染物的輸出量,t/a;E為類土地利用類型中污染物的輸出系數(shù),t/(km2·a);A為類土地利用類型的面積,km2;為降雨直接帶來的污染物負荷,t/a.

傳統(tǒng)的輸出系數(shù)模型忽略了降雨、地形等條件對模擬結(jié)果的影響,而降雨、地形和土地利用方式在非點源污染物的輸移過程中起著重要作用[20],陳學(xué)凱等[21]在輸出系數(shù)模型的基礎(chǔ)上,加入降雨、坡度影響因子對模型進行優(yōu)化,并驗證了改進輸出系數(shù)模型的合理性,李華林等[22]通過模擬分析,確定改進的輸出系數(shù)模型模擬結(jié)果更符合流域非點源污染特征.由于河套灌區(qū)干旱少雨,直接由降雨帶來的污染較小,故而在模擬中忽略不計.因此本文構(gòu)建改進的輸出系數(shù)模型,計算公式為:

式中:為降水影響因子;為地形影響因子.

輸出系數(shù)模型模擬出來的結(jié)果是污染物的輸出量,在降水沖刷和坡度作用下向河道輸移的過程中,會受到植被攔截、自然沉降以及各種生化作用的影響,導(dǎo)致產(chǎn)生的污染物不可能全部到達附近排水溝以及烏梁素海,因此需要引入污染物入河系數(shù)來表征污染物遷移過程中的衰減影響[23-24].公式如下:

L

=

ljLj

(3)

式中:L為污染物輸出量,t/a;為污染物入河量,t/a;λ為污染物的入河系數(shù).

圖2 河套灌區(qū)2021年降雨量和降雨影響因子空間分布

1.3.1 降雨影響因子的確定 降雨量對氮、磷流失量具有較為顯著的影響,降雨影響因子主要受時間不均勻性影響因子t和空間不均勻性因子s的影響[25-26],計算公式為:

式中:為研究區(qū)給定年年降雨總量,mm;R為研究區(qū)多年平均年降雨量,mm;r為子流域給定年年降雨總量,mm;r為子流域多年平均年降雨量,mm.

圖3 河套灌區(qū)地形影響因子

利用河套灌區(qū)2021年7個氣象站點的降雨量數(shù)據(jù),經(jīng)過插值得到研究區(qū)年降雨量的空間分布圖(圖2),研究區(qū)2021年降雨量為57.219～192.798mm,整體呈現(xiàn)由東至西逐漸降低的空間分布特征.

根據(jù)巴彥淖爾市統(tǒng)計年鑒獲得研究區(qū)多年平均降雨量159.5mm,以研究區(qū)鄉(xiāng)鎮(zhèn)為輸出單元進行劃分,帶入公式(4)進行計算,插值得到各個鄉(xiāng)鎮(zhèn)的降雨影響因子,研究區(qū)降雨影響因子為0.420～0.694.根據(jù)有關(guān)研究成果,為0.6104[28],對河套灌區(qū)的DEM數(shù)據(jù)進行提取計算,獲得研究區(qū)坡度為1.801°～8.217°,平均坡度為3.693°,帶入公式(6),借助ArcGIS10.7獲得研究區(qū)地形影響因子的空間分布圖(圖3).研究區(qū)地形影響因子為0.649～1.618,空間分布上與坡度分布一致,研究區(qū)西北部的坡度影響因子較大.

1.3.2 地形影響因子的確定 在匯流過程中,坡度主要通過影響徑流量來影響其攜帶的污染物的流失量,且坡度與地面徑流量呈正相關(guān)關(guān)系,坡度越大地表水流速度越大,地表徑流越大.相關(guān)研究表明,徑流量可以表示為坡度的冪函數(shù)與常量的乘積[27],表達式為:

通過徑流量與地形的關(guān)系,則地形影響因子可以表示為:

式中:為徑流量,m3/s;為坡度,°;,為常量;S為各鄉(xiāng)鎮(zhèn)坡度;為研究區(qū)平均坡度.

1.3.3 輸出系數(shù)的確定 確定合理的污染物輸出系數(shù)是輸出系數(shù)模型的關(guān)鍵,輸出系數(shù)的確定有3種常用方法,即文獻分析法、試驗?zāi)M法和水文統(tǒng)計法[29].鑒于河套灌區(qū)無土地利用輸出系數(shù)實驗數(shù)據(jù),本研究采用文獻分析法[30-31], 盡可能選擇與河套灌區(qū)相符或相近研究區(qū)域的輸出系數(shù),計算其平均值,最終得到本研究的不同土地利用類型的輸出系數(shù),見表2.

表2 河套灌區(qū)不同土地利用類型TN、TP輸出系數(shù)

1.4 模型的驗證

將模型模擬值與實測值進行對比分析,通過相對誤差(e)這一指標(biāo)來驗證改進的輸出系數(shù)模型的適用性.計算方法如下:

式中:e為相對誤差;為模擬值;0為監(jiān)測值.若e為正則說明模擬值偏大,若e為負則說明模擬值偏小,e的絕對值越小,則模型模擬效果較好,反之則模型模擬效果較差.

2 結(jié)果與分析

2.1 河套灌區(qū)不同土地利用類型非點源TN、TP污染負荷分析

在選定的TN、TP輸出系數(shù)以及不考慮流域損失的情況下,根據(jù)模型模擬得到河套灌區(qū)不同土地利用類型所產(chǎn)生的非點源污染負荷(圖4).2021年河套灌區(qū)不同土地利用類型下所產(chǎn)生的非點源TN、TP污染負荷分別為4419.251t和456.414t,其中污染物輸出負荷量最高的土地利用類型為耕地,TN、TP污染負荷分別為3520.998t、407.125t,占河套灌區(qū)不同土地利用類型下非點源TN、TP污染負荷的79.674%和89.201%.河套灌區(qū)其他土地利用類型TN污染負荷的貢獻率大小依次為草地>未利用地>水域>建設(shè)用地>林地,TP污染負荷的貢獻率大小依次為草地>建設(shè)用地>未利用地>水域>林地.

圖4 不同土地利用類型下TN、TP污染負荷占比

2.2 河套灌區(qū)農(nóng)田非點源TN、TP污染負荷分析

由圖5可知,河套灌區(qū)農(nóng)田非點源TN、TP污染負荷空間分布規(guī)律基本一致,但整體上分布不均,呈現(xiàn)出污染負荷局部集中、靠近水體、排干的鄉(xiāng)鎮(zhèn)污染負荷高.很可能是由于該區(qū)域的土地利用方式以耕地為主,從事較多的農(nóng)事活動,且農(nóng)田排水溝是農(nóng)業(yè)非點源氮磷等污染物運輸和遷移至水體的主要途徑[44],因此污染負荷高,而灌區(qū)西部污染負荷比較低的區(qū)域則以草地和未利用地居多,農(nóng)用地少,氮磷等營養(yǎng)物質(zhì)的輸入量也相對較低.其中農(nóng)田非點源TN、TP污染負荷最高的幾個鄉(xiāng)鎮(zhèn)依次為新華鎮(zhèn)、新安鎮(zhèn)和塔爾湖鎮(zhèn),TN、TP輸出量分別為316.658, 236.725, 199.344t/a和36.615, 27.371, 23.050t/a,占農(nóng)田非點源TN、TP輸出總量的21.378%.

圖5 農(nóng)田非點源TN、TP污染負荷空間分布

2.3 河套灌區(qū)農(nóng)田非點源TN、TP污染負荷強度分析

為避免流域面積對污染負荷的影響,本文采用污染負荷強度法分析不同鄉(xiāng)鎮(zhèn)的TN、TP污染負荷強度(圖6).TN污染負荷強度在0.020～5.826kg/hm2之間,平均值為3.042kg/hm2,TP污染負荷強度在0.002～0.674kg/hm2之間,平均值為0.352kg/hm2,且TN、TP污染負荷強度空間分布一致.污染負荷強度較高的地區(qū)依次為團結(jié)鎮(zhèn)、新華鎮(zhèn)和蠻會鎮(zhèn),表明這幾個鄉(xiāng)鎮(zhèn)對整個研究區(qū)來說,更易產(chǎn)生氮、磷污染.最高污染負荷強度與最高污染負荷區(qū)域不一致,污染負荷強度最高的區(qū)域為團結(jié)鎮(zhèn),TN、TP污染負荷強度分別為5.826和0.674kg/hm2.最低污染負荷強度地區(qū)為臨河區(qū)各辦事處.分析其原因,團結(jié)鎮(zhèn)和蠻會鎮(zhèn)土地利用以耕地為主,耕地面積分別占鄉(xiāng)鎮(zhèn)總面積的87.726%和82.936%,且總排干流經(jīng)該鄉(xiāng)鎮(zhèn)導(dǎo)致TN、TP污染強度較高,而各辦事處土地利用類型以居民建設(shè)用地為主,幾乎不存在農(nóng)田非點源輸出.

圖6 農(nóng)田非點源TN、TP污染負荷強度空間分布

2.4 河套灌區(qū)農(nóng)田非點源TN、TP入河量分析

圖7 農(nóng)田非點源TN、TP入河量空間分布

各地農(nóng)田非點源TN、TP的入河系數(shù)各有差異,一般南方地區(qū)為0.1～0.3,而北方地區(qū)由于氣候干旱,降雨量低于南方地區(qū),一般取值0.03～0.1.考慮到河套灌區(qū)農(nóng)田的傳統(tǒng)灌溉方式等措施,結(jié)合文獻資料[23-24,45],選取河套灌區(qū)農(nóng)田非點源TN入河系數(shù)為0.11,TP入河系數(shù)為0.04.

圖7為河套灌區(qū)農(nóng)田非點源TN、TP入河量的空間分布圖,河套灌區(qū)農(nóng)田非點源TN、TP的入河量空間分布特征一致,整體上與污染物輸出負荷具有相似的空間分布,均為靠近水體、排干的鄉(xiāng)鎮(zhèn)污染貢獻量大,農(nóng)田非點源TN、TP入河總量分別為387.310和16.285t/a,其中新華鎮(zhèn)、新安鎮(zhèn)、塔爾湖、西小召鎮(zhèn)和隆興昌鎮(zhèn)對TN、TP入湖量的貢獻較大,占污染物入河總量的30.852%,是今后農(nóng)田非點源污染防治的關(guān)鍵區(qū)域.

2.5 模型驗證分析

為了驗證改進輸出系數(shù)模型的合理性,本研究選取六、七排干及皂沙排干匯水區(qū)所在鄉(xiāng)鎮(zhèn)作為模型的驗證區(qū),根據(jù)《五原縣水體達標(biāo)方案》獲取六、七排干、皂沙排干溝的排水量分別為1065.33萬m3、1354.69萬m3和381.3萬m3,將污染物濃度與排水量進行乘積后得到驗證區(qū)TN、TP實測值.

表3 改進的輸出系數(shù)模型模擬精度

結(jié)合公式(7),結(jié)果如表3所示,驗證區(qū)農(nóng)田非點源TN、TP入河量實測值為73.023t/a和2.457t/a,改進輸出系數(shù)模型模擬值為69.559t/a和2.925t/a.模型模擬結(jié)果與實測值較為相近,相對誤差分別為-4.743%和19.037%,且由降雨影響因子和地形影響因子同時改進的輸出系數(shù)模型提高了模型模擬精度,表明所選用的系數(shù)合理可靠,模型具有較好的模擬能力.

3 討論

氮磷流失的數(shù)據(jù)源是基于種植面積的估算,計算子流域的負荷可以有效降低估算帶來的誤差[46],所以本研究以鄉(xiāng)鎮(zhèn)為輸出單元,基于改進的輸出系數(shù)模型并借助GIS技術(shù),進行河套灌區(qū)2021年農(nóng)田非點源TN、TP污染負荷估算.結(jié)果表明,耕地是不同土地利用類型中非點源污染貢獻量最大的土地利用類型,孫海軍等[47]研究了太湖流域北部山區(qū)小流域非點源污染,說明了種植業(yè)是該流域的主要污染源.河套灌區(qū)農(nóng)田非點源TN、TP污染負荷空間分布一致,靠近水體、排干的鄉(xiāng)鎮(zhèn)農(nóng)田非點源TN、TP貢獻量最高,這與王文章等[48]對射洪縣非點源評估的結(jié)果相符.TN、TP污染負荷最高的幾個鄉(xiāng)鎮(zhèn)依次為新華鎮(zhèn)、新安鎮(zhèn)和塔爾湖鎮(zhèn),這是由于該區(qū)域的土地利用方式以耕地為主,耕地面積較大,從事較多的農(nóng)事活動,這與石慶玲等[49]通過探究農(nóng)業(yè)非點源污染負荷對土地利用的響應(yīng)關(guān)系,證明污染負荷與耕地面積所占比重有很大關(guān)系結(jié)論一致.

目前,將輸出系數(shù)模型應(yīng)用于農(nóng)業(yè)非點源污染負荷估算的研究多集中于污染物輸出或者單一階段的污染物入河估算[50],但農(nóng)業(yè)非點源對水環(huán)境的污染要經(jīng)歷一定階段[51],故本文結(jié)合污染物入河系數(shù),將污染物輸出與污染物入河聯(lián)系起來,結(jié)果表明,2021年河套灌區(qū)農(nóng)田非點源TN、TP污染負荷為3520.998t,407.125t,TN貢獻量約為TP的8.65倍,李娜等[52]在長春市新立城水庫的研究中,若不考慮牛的養(yǎng)殖,TN與TP的排放比達到8.32,本結(jié)果與其研究結(jié)果相似.TN、TP污染負荷強度平均值分別為3.042和0.352kg/hm2,吳用[53]利用SWAT模型對烏梁素海流域的非點源污染進行模擬,結(jié)果表明,平原區(qū)流域中、東部氮和磷元素的輸出平均值分別為5.6,2.3kg/hm2和0.25,0.86kg/hm2,本研究與其研究結(jié)果相符.農(nóng)田非點源TN、TP入河量分別為387.310 和16.285t,管玉玲[54]利用排污系數(shù)法估算烏梁素海流域2001～2014年農(nóng)田非點源污染物入河量,其中TN、TP入河量分別為673.18～918.90t和28.27～ 38.59t.總氮總磷入河量比本研究偏高,可能是因為近年來烏梁素海各項水質(zhì)指標(biāo)入河污染負荷均呈現(xiàn)出不同程度的下降趨勢[55],且其研究區(qū)域為整個烏梁素海流域,范圍較大.考慮到研究區(qū)缺乏TN、TP輸出系數(shù)的監(jiān)測數(shù)據(jù),模型所使用的輸出系數(shù)參考其他地理環(huán)境相似流域的平均值,模擬結(jié)果會具有一定的誤差,以上分析表明本研究結(jié)果具有一定的可信度,今后,在相關(guān)數(shù)據(jù)的支撐下,基于本研究方法所得出的結(jié)果將會更加準(zhǔn)確.

4 結(jié)論

4.1 耕地是河套灌區(qū)不同土地利用類型中非點源TN、TP污染負荷量最大的土地利用類型,其他土地利用類型對TN污染負荷的貢獻率大小依次為草地>未利用地>水域>建設(shè)用地>林地,TP污染負荷的貢獻率大小依次為草地>建設(shè)用地>未利用地>水域>林地.

4.2 河套灌區(qū)2021年農(nóng)田非點源TN、TP污染負荷分別為3520.998和407.125t,污染負荷強度平均值分別為3.042 和0.352kg/hm2,污染貢獻量最高的幾個鄉(xiāng)鎮(zhèn)依次為新華鎮(zhèn)、新安鎮(zhèn)和塔爾湖鎮(zhèn),占農(nóng)田非點源TN、TP輸出總量的21.378%.

4.3 結(jié)合污染物入河系數(shù),估算農(nóng)田非點源TN、TP進入烏梁素海的負荷量,河套灌區(qū)農(nóng)田非點源TN、TP的入河量空間分布特征一致,整體上與污染物輸出負荷具有相似的空間分布,均為靠近水體、排干的鄉(xiāng)鎮(zhèn)污染貢獻量大,農(nóng)田非點源TN、TP入河總量分別為387.310和16.285t/a,其中新華鎮(zhèn)、新安鎮(zhèn)、塔爾湖鎮(zhèn)、西小召鎮(zhèn)、隆興昌鎮(zhèn)、團結(jié)鎮(zhèn)和蠻會鎮(zhèn)是河套灌區(qū)農(nóng)田非點源污染的重點污染防治鄉(xiāng)鎮(zhèn).

4.4 基于降雨和地形影響因子修正改進的輸出系數(shù)模型,對農(nóng)田非點源TN、TP污染負荷模擬的相對誤差分別為-4.743%和19.037%,模型模擬能力比較好.

[1] 王先甲,唐金鵬,李長杰,等.瑪納斯河灌區(qū)分水配水管理信息系統(tǒng) [J]. 干旱區(qū)資源與環(huán)境, 2006,20(4):127-132.

Wang X J, Tang J P, LI C J, et al. Management information system of water assignment in arid irrigation region [J]. Journal of Arid Land Resources and Environment, 2006,20(4):127-132.

[2] Su Y, He S, Wang K, et al. Quantifying the sustainability of three types of agricultural production in China: An emergy analysis with the integration of environmental pollution [J]. Journal of Cleaner Production, 2020,252:119650.

[3] 第二次全國污染源普查公報 [J]. 環(huán)境保護, 2020,48(18):8-10.

The Second National Pollution Source Census Bulletin [J]. Environmental Protection, 2020,48(18):8-10.

[4] 楊林章,馮彥房,施衛(wèi)明,等.我國農(nóng)業(yè)面源污染治理技術(shù)研究進展 [J]. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2013,21(1):96-101.

Yang L Z, Feng Y F, Shi W M, et al. Review of the advances and development trends in agricultural non-point source pollution control in China [J]. Chinese Journal of Eco-Agriculture, 2013,21(1):96-101.

[5] Johnes P J. Evaluation and management of the impact of land use change on the nitrogen and phosphorus load delivered to surface waters: the export coefficient modelling approach [J]. Journal of Hydrology, 1996,183(3/4):323-349.

[6] Chen H Y, Teng Y G, Wang J S. Load estimation and source apportionment of nonpoint source nitrogen and phosphorus based on integrated application of SLURP model, ECM, and RUSLE: a case study in the Jinjiang River, China [J]. Environmental Monitoring and Assessment, 2013,185:2009-2021.

[7] Li S S, Zhang L, Du Y, et al. Evaluating phosphorus loss for watershed management: Integrating a weighting scheme of watershed heterogeneity into export coefficient model [J]. Environmental Modeling & Assessment, 2016,21:657-668.

[8] 楊彥蘭,申麗娟,謝德體,等.基于輸出系數(shù)模型的三峽庫區(qū)(重慶段)農(nóng)業(yè)面源污染負荷估算 [J]. 西南大學(xué)學(xué)報(自然科學(xué)版), 2015,37: 112-119.

Yang Y L, Shen L J, Xie D T, et al. Estimation of pollution loads from agricultural nonpoint sources in Three Gorges Reservoir area (Chongqing) based on the export coefficient modeling approach [J]. Journal of Southwest University (Natural Science Edition), 2015,37: 112-119.

[9] 荊延德,張華美,孫笑笑.基于輸出系數(shù)模型的南四湖流域非點源污染輸出風(fēng)險評估 [J]. 水土保持通報, 2017,37(3):270-274.

JING Y D, ZHANG H M, SUN X X. Risk assessment of non-point source pollution in Nansihu lake basin using output coefficient model [J]. Bulletin of soil and water conservation, 2017,37(3):270-274.

[10] 蔡 明,李懷恩,莊詠濤,等.改進的輸出系數(shù)法在流域非點源污染負荷估算中的應(yīng)用 [J]. 水利學(xué)報, 2004,49(7):40-45.

CAI M, LI H E, ZHUANG Y T, et al. Application of modified export coefficient method in polluting load estimation of non-point source pollution [J]. Journal of Hydraulic Engineering, 2004,49(7):40-45.

[11] Junrui, C, Yuemin, G, Z., Z D, et al. Modeling the sources and retention of phosphorus nutrient in a coastal river system in China using SWAT [J]. Journal of Environmental Management, 2021,278: 111556.

[12] 王 萌,王敬賢,劉 云,等.湖北省三峽庫區(qū)1991～2014年農(nóng)業(yè)非點源氮磷污染負荷分析 [J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2018,37(2):294-301.

WANG M, WANG J X, LIU Y, et al. Analysis of nitrogen and phosphorus pollution loads from agricultural non-point sources in the Three Gorges Reservoir of Hubei Province from 1991 to 2014 [J]. Journal of Agro-Environment Science, 2018,37(2):294-301.

[13] Wang, R, Kalin, L. Combined and synergistic effects of climate change and urbanization on water quality in the Wolf Bay watershed, southern Alabama [J]. Journal of Environmental Sciences:china, 2018, 64:107-121.

[14] Kaixia, J, Zhaofu, L, Chuan, L, et al. The reduction effects of riparian reforestation on runoff and nutrient export based on AnnAGNPS model in a small typical watershed, China. [J]. Environmental science and pollution research international, 2019,26:5934-5943.

[15] 丁春蓮,尹曉敏,田志強,等.河套灌區(qū)總排干溝入烏梁素海面源污染負荷的便捷評估方法探討 [J]. 內(nèi)蒙古水利, 2016,37(6):7-8.

Ding C L, Yin X M, Tian Z Q, et al. Study on convenient assessment method of surface pollution load from main drainage ditch into Wuliangsuhai in Hetao Irrigation Distric [J]. Inner Mongolia Water Resources, 2016,37(6):7-8.

[16] 魏國孝,孫繼成,朱 鋒.內(nèi)蒙古河套灌區(qū)農(nóng)業(yè)面源污染及防治對策 [J]. 中國水土保持, 2009,30(8):27-29.

Wei G X, Sun J C, Zhu F. Agricultural non-point source pollution and its control measures in Hetao Irrigation District of Inner Mongolia [J]. Soil and Water Conservation in China, 2009,30(8):27-29.

[17] Wu L, Gao J E, Ma X Y, et al. Application of modified export coefficient method on the load estimation of non-point source nitrogen and phosphorus pollution of soil and water loss in semiarid regions [J]. Environmental Science and Pollution Research International, 2015,22:10647-10660.

[18] 歐陽威,郭波波,張 璇,等.北方典型灌區(qū)不同灌期農(nóng)田系統(tǒng)中氮素遷移特征分析 [J]. 中國環(huán)境科學(xué), 2013,33(1):123-131.

Ouyang W, Guo B B, Zhang X, et al. Transfer characteristics of soil nitrogen in northern typical irrigation area under different irrigation periods [J]. China Environmental Science, 2013,33(1):123-131.

[19] 黃永江,屈忠義.基于SPA-VFS耦合模型的河套灌區(qū)灌溉用水效率分級 [J]. 農(nóng)業(yè)工程學(xué)報, 2021,37(11):85-92.

Huang Y J, Qu Z Y. Classifying irrigation water use efficiency in Hetao Irrigation District using SPA-VFS coupling model [J]. Transactions of the Chinese Society of Agricultural Engineering, 2021,37(11):85-92.

[20] Yang G X, Best E P H, Whiteaker T, et al. A screening-level modeling approach to estimate nitrogen loading and standard exceedance risk, with application to the Tippecanoe River watershed, Indiana [J]. Journal of Environmental Management, 2014,135:1-10.

[21] 陳學(xué)凱,劉曉波,彭文啟,等.程海流域非點源污染負荷估算及其控制對策 [J]. 環(huán)境科學(xué), 2018,39(1):77-88.

Chen X K, Liu X B, Peng W Q, et al. Estimation of and control strategies for pollution loads from non-point sources in the Chenghai watershed [J]. Environment Science, 2018,39(1):77-88.

[22] 李華林,張守紅,于佩丹,等.基于改進輸出系數(shù)模型的非點源污染評估及關(guān)鍵源區(qū)識別:以北運河上游流域為例 [J]. 環(huán)境科學(xué), 2023, 44(11):6194-6204.

Li H L, Zhang S H, Yu P D, et al. Estimation and critical source areas identification of non-point source pollution based on improved export coefficient models: A case study of the upper Beiyun River Basin [J]. Environmental Science, 2023,44(11):6194-6204.

[23] Chen D J, Lu J, Huang H, et al. Stream nitrogen sources apportionment and pollution control scheme development in an agricultural watershed in eastern China [J]. Environmental Management, 2013,52:450-466.

[24] Cheng X, Chen L D, Sun R H, et al. An improved export coefficient model to estimate non-point source phosphorus pollution risks under complex precipitation and terrain conditions [J]. Environmental Science and Pollution Research, 2018,25:20946-20955.

[25] Ding X W, Shen Z Y, Hong Q, et al. Development and test of the export coefficient model in the upper reach of the Yangtze River [J]. Journal of Hydrology, 2009,383(3/4):233-244.

[26] Chen X K, Liu X B, Peng W Q, et al. Non-point source nitrogen and phosphorus assessment and management plan with an improved method in data-poor regions [J]. Water, 2017,10(1):17.

[27] Liu Q Q, Singh V P. Effect of microtopography, slope length and gradient, and vegetative cover on overland flow through simulation [J]. Journal of Hydrologic Engineering, 2004,9(5):375-382.

[28] Li Y, Wang C, Tang H L. Research advances in nutrient runoff on sloping land in watersheds [J]. Aquatic Ecosystem Health & Management, 2006,9(1):27-32.

[29] 張洪波,李 俊,黎小東,等.缺資料地區(qū)農(nóng)村面源污染評估方法研究 [J]. 四川大學(xué)學(xué)報(工程科學(xué)版), 2013,45(6):58-66.

Zhang H B, Li J, Li X D, et al. Study on rural non-point source pollution assessment method of regions with sparse data [J]. Journal of Sichuan University (Advanced Engineering Sciences), 2013,45(6):58- 66.

[30] 胡富昶,敖天其,胡 正,等.改進的輸出系數(shù)模型在射洪縣的非點源污染應(yīng)用研究 [J]. 中國農(nóng)村水利水電, 2019,61(6):78-82.

Hu F C, Ao T Q, Hu Z, et al. Applied research on the improved export coefficient model in non-point source pollution in Shehong County [J]. China Rural Water and Hydropower, 2019,61(6):78-82.

[31] 楊 雯,敖天其,王文章,等.基于輸出系數(shù)模型的瓊江流域(安居段)農(nóng)村非點源污染負荷評估 [J]. 環(huán)境工程, 2018,36(10):140-144.

Yang W, Ao T Q, Wang W Z, et al. Estimation of pollution loads from rural non-point sources in Qiongjiang river basin (Anju district) based on export coefficient modeling approach [J]. Environmental Engineering, 2018,36(10):140-144.

[32] 李強坤.青銅峽灌區(qū)農(nóng)業(yè)非點源污染負荷及控制措施研究 [D]. 西安:西安理工大學(xué), 2010.

Li Q K. Study on the loads and control measures of agricultural non-point source pollution in Qingtongxia irrigation district [D]. Xi'an: Xi'an University of Technology, 2010.

[33] Duan M Q, Du X, Peng W Q, et al. Quantitative assessment of background pollutants using a modified method in data-poor regions [J]. Environmental Monitoring and Assessment, 2020,192:1-15.

[34] 王有霄,鐘萍麗,于 格,等.膠州灣氮?磷非點源污染負荷估算及時空分析 [J]. 中國海洋大學(xué)學(xué)報(自然科學(xué)版), 2019,49(2):85-97.

Wang Y X, Zhong P L, Yu G, et al. Estimation of N,P non-point source pollution loads and analysis of spatial and temporal characteristies in the Jiaozhou Bay [J]. Periodical of Ocean University of China, 2019, 49(2):85-97.

[35] 張愛平.寧夏黃灌區(qū)稻田退水氮磷污染特征研究 [D]. 北京:中國農(nóng)業(yè)科學(xué)院, 2009.

Zhang A P. Study on the characteristics of nitrogen and phosphorus pollution in the return flow of rice paddy in Ningxia irrigation region [D]. Beijing: Chinese Academy of Agriculture Sciences, 2009.

[36] Du X Z, Su J J, Li X Y, et al. Researchers from Chinese academy of sciences detail findings in environmental pollution (Modeling and evaluating of non-point source pollution in a Semi-Arid Watershed: Implications for watershed management) [J]. Ecology Environment & Conservation, 2016,44(3):247-255.

[37] 程 先,陳利頂,孫然好.考慮降水和地形的京津冀水庫流域非點源污染負荷估算 [J]. 農(nóng)業(yè)工程學(xué)報, 2017,33(4):265-272.

Cheng X, Chen L D, Sun H R. Estimation of non-point source pollution loads in the Beijing-Tianjin-Hebei reservoir basin considering precipitation and topography [J]. Transactions of the Chinese Society of Agricultural Engineering, 2017,33(4):265-272.

[38] 李 根,毛 鋒.我國水土流失型非點源污染負荷及其經(jīng)濟損失評估 [J]. 中國水土保持, 2008,39(2):9-11.

Li G, Mao F. Assessment of soil erosion-type non-point source pollution loads and their economic losses in China [J]. Soil and Water Conservation in China, 2008,39(2):9-11.

[39] 宋 嘉,李懷恩,李家科,等.丹江流域陜西段農(nóng)業(yè)非點源污染負荷估算 [J]. 中國農(nóng)村水利水電, 2020,62(11):67-72.

Song J, Li H E, Li J K, et al. Load estimation of agricultural non-point source pollution in Shaanxi Section of Danjiang River Basin [J]. China Rural Water and Hydropower, 2020,62(11):67-72.

[40] 楊淑靜,張愛平,楊正禮,等.寧夏灌區(qū)農(nóng)業(yè)非點源污染負荷估算方法初探 [J]. 中國農(nóng)業(yè)科學(xué), 2009,42(11):3947-3955.

Yang S J, Zhang A P, Yang Z L, et al. Agricultural non-point source pollution in Ningxia irrigation district and preliminary study of load estimation methods [J]. Scientia Agricultura Sinica, 2009,42(11): 3947-3955.

[41] 丁 洋,趙進勇,董 飛,等.媯水河流域農(nóng)業(yè)非點源污染負荷估算與分析 [J]. 水利水電技術(shù), 2020,51(1):139-146.

Ding Y, Zhao J Y, Dong F, et al. Estimation and analysis of agricultural non-point source pollution load in Guishui River Watershed [J]. Water Resources and Hydropower Engineering, 2020,51(1):139-146.

[42] 龐樹江,王曉燕.流域尺度非點源總氮輸出系數(shù)改進模型的應(yīng)用 [J]. 農(nóng)業(yè)工程學(xué)報, 2017,33(18):213-223.

Pang S J, Wang X Y. Application of modified diffuse total nitrogen export coefficient model at watershed scale [J]. Transactions of the Chinese Society of Agricultural Engineering, 2017,33(18):213-223.

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

Liu Y Q, Yang Y L, Li F H. Estimation of pollution loads from agricultural nonpoint sources in Beijing region based on export coefficient modeling approach [J]. Transactions of the Chinese Society of Agricultural Engineering, 2011,27(7):7-12.

[44] 胡宜剛,吳 攀,趙 洋,等.寧蒙引黃灌區(qū)農(nóng)田排水溝渠水質(zhì)特征 [J]. 生態(tài)學(xué)雜志, 2013,32(7):1730-1738.

Hu Y G, Wu P, Zhao Y, et al. Water quality of cropland drainage ditches in the Yellow River irrigation regions of Ningxia and Inner Mongolia, China [J]. Chinese Journal of Ecology, 2013,32(7):1730- 1738.

[45] 段 揚,蔣洪強,吳文俊,等.基于改進輸出系數(shù)模型的非點源污染負荷估算——以嫩江流域為例 [J]. 環(huán)境保護科學(xué), 2020,46(4):48-55.

Duan Y, Jiang H Q, Wu W J, et al. Estimation of non-point source pollution load based on improved export coefficient model : A case study of Nenjiang Watershed [J]. Environmental Protection Science, 2020,46(4):48-55.

[46] 左琳琳.巢湖流域農(nóng)田面源污染負荷估算及其風(fēng)險評估 [D]. 西安:西安科技大學(xué), 2018.

Zuo L L. Farmland non-point source pollution load estimation and risk assessment in Chaohu Lake Basin [D]. Xi'an: Xi'an University of Science and Technology, 2018.

[47] 孫海軍,吳家森,姜培坤,等.浙北山區(qū)典型小流域農(nóng)村面源污染現(xiàn)狀調(diào)查與治理對策 [J]. 中國農(nóng)學(xué)通報, 2011,27(20):258-264.

Sun H J, Wu J S, Jiang P K, et al. A survey about rural non-point pollution in nountain region of northern Zhejiang Province and the countermeasures [J]. Chinese Agricultural Science Bulletin, 2011, 27(20):258-264.

[48] 王文章,敖天其,史小春,等.基于輸出系數(shù)模型的射洪縣農(nóng)村面源污染負荷估算 [J]. 環(huán)境工程, 2018,36(1):173-177.

Wang W Z, Ao T Q, Shi X C, et al. Estimation of pollution loads from rural non-point sources in Shehong County Based on export coefficient modeling approach [J]. Environmental Engineering, 2018, 36(1):173-177.

[49] 石慶玲.烏梁素海流域非點源污染負荷對LUCC的響應(yīng)研究 [D]. 濟南:濟南大學(xué), 2013.

ShiI Q L. The response of non-point source pollution to LUCC in Wuliangsuhai Watershed [D]. Jinan: University of Jinan, 2013.

[50] 桂平婧,王 豐,李善樸,等.基于階段輸出系數(shù)模型的農(nóng)業(yè)非點源污染負荷估算與評價——以四川省為例 [J]. 浙江農(nóng)業(yè)學(xué)報, 2016, 28(1):110-118.

Gui P J, Wang F, Li S P, et al. Estimation and evaluation of agricultural non-point source pollution load by stage export coefficient model: A case study in Sichuan Province [J]. Acta Agricuhurae Zhejiangensis, 2016,28(1):110-118.

[51] Wang J C, Wu Y Q, Hu A Y, et al. Application and establishment model of non-point source pollution based on statistical data[C]. Xi’an:International Symposium on Water Resource and Environmental Protection, 2011.

[52] 李 娜,韓維崢,沈夢楠,等.基于輸出系數(shù)模型的水庫匯水區(qū)農(nóng)業(yè)面源污染負荷估算 [J]. 農(nóng)業(yè)工程學(xué)報, 2016,32(8):224-230.

Li N, Han W Z, Shen M N, et al. Load evaluation of non-point source pollutants from reservoir based on export coefficient modeling [J]. Transactions of the Chinese Society of Agricultural Engineering, 2016,32(8):224-230.

[53] 吳 用.內(nèi)蒙古河套地區(qū)流域非點源污染模擬及區(qū)域污染物總量控制研究 [D]. 呼和浩特:內(nèi)蒙古農(nóng)業(yè)大學(xué), 2017.

Wu Y. Modeling of non-point sources pollution and total amount control of pollutants in Hetao irrigaton district of the Inner Mongolia [D]. Huhehot: Inner Mongolia Agricultural University, 2017.

[54] 管玉玲.烏梁素海流域面源污染負荷入湖量計算 [J]. 內(nèi)蒙古科技與經(jīng)濟, 2020,39(2):49-52.

Guan Y L. Calculation of surface source pollution loads into the lake in the Wuliangsuhai Lake [J]. Inner Mongolia Science Technologr & Economy, 2020,39(2):49-52.

[55] 楊文煥,尹 強,鐘清濤,等.2013～2019年烏梁素海排干入湖污染負荷與湖區(qū)水質(zhì)的響應(yīng)關(guān)系 [J]. 華北水利水電大學(xué)學(xué)報(自然科學(xué)版), 2021,42(1):7-15.

Yang W H, Yi Q, Zhong Q T, et al. Study on the response relationship between the pollution load of arterial drainage of Lake Wuliangsuhai and the water quality of the lake from 2013 to 2019 [J]. Journal of North China University of Water Resources and Electric Power (Natural Science Edition), 2021,42(1):7-15.

Estimation of nitrogen and phosphorus pollution loads from non-point sources in farmland of Hetao Irrigation District, Inner Mongolia.

WANG Yan, JIAO Yan*, YANG Wen-zhu, YAN Ying-chao, WU Hong, LING Ling, SHI Yue-chao

(Inner Mongolia Key Laboratory of Environmental Chemistry, College of Chemistry and Environmental Science, Inner Mongolia Normal University, Hohhot 010022, China)., 2023,43(12)：6551～6560

To understand the current situation of farmland non-point source pollution in the Hetao irrigation district of Inner Mongolia, the traditional export coefficient model (ECM) was improved, and the farmland non-point source nitrogen and phosphorus pollution load was quantified by precipitation impact factor () and terrain impact factor () in the area. The farmland non-point source nitrogen and phosphorus discharged into the river (Wuliangsuhai) was estimated, and the critical source areas of pollution was identified by the pollution river influx coefficient (). The total export load of non-point source total nitrogen (TN) and total phosphorus (TP) from farmland were 3520.998 and 407.125t in the Hetao irrigation district in 2021, accounting for 79.674% and 89.201% of the total pollution load from different land use types, respectively. The spatial distribution of the export loads for TN and TP of farmland non-point sources was consistent, and the contribution in towns near drainage was the highest. The export loads of TN and TP pollution were 316.658, 236.725, 199.344t/a and 36.615, 27.371, 23.050t/a in Xinhua, Xin'an and Talhu township, respectively. The amount of TN and TP discharged into the river were 387.310 and 16.285t/a from non-point sources of farmland. The towns of Xinhua, Xin'an, Talhu, Longxingchang, Tuanjie and Manhui were key areas for pollution prevention and control in the Hetao irrigation district. The rationality of the improved export coefficient model with relative errors of -4.743% and 19.037% had been verified by comparing with measured values.

Hetao irrigation district；farmland non-point source pollution；improved export coefficient model；TN；TP

X502

A

1000-6923(2023)12-6551-10

王 艷,焦 燕,楊文柱,等.內(nèi)蒙古河套灌區(qū)農(nóng)田非點源氮磷污染負荷估算 [J]. 中國環(huán)境科學(xué), 2023,43(12):6551-6560.

Wang Y, Jiao Y, Yang W Z, et al. Estimation of nitrogen and phosphorus pollution loads from non-point sources in farmland of Hetao Irrigation District, Inner Mongolia [J]. China Environmental Science, 2023,43(12):6551-6560.

2023-04-24

國家自然科學(xué)基金資助項目(42175038);內(nèi)蒙古自治區(qū)杰出青年基金資助項目(2022JQ02);內(nèi)蒙古自治區(qū)高等學(xué)校青年科技英才支持計劃項目(NJYT23041);2022年蒙古自治區(qū)重點研發(fā)和成果轉(zhuǎn)化計劃項目(2022YFHH0035);內(nèi)蒙古師范大學(xué)基本科研業(yè)務(wù)費專項資金資助項目(2022JBTD009);內(nèi)蒙古師范大學(xué)高層次人才科研啟動項目(2020YJRC056)

* 責(zé)任作者, 教授, jiaoyan@imnu.edu

王 艷(1996-),女,河南周口人,內(nèi)蒙古師范大學(xué)碩士研究生,主要從事農(nóng)業(yè)非點源污染研究.發(fā)表文章1篇.2469736853@qq.com.