馮愛萍,吳傳慶,王雪蕾*,王洪亮,周亞明,趙 乾

海河流域氮磷面源污染空間特征遙感解析

馮愛萍1,吳傳慶1,王雪蕾1*,王洪亮2,周亞明1,趙 乾1

(1.生態(tài)環(huán)境部衛(wèi)星環(huán)境應(yīng)用中心,國(guó)家環(huán)境保護(hù)衛(wèi)星遙感重點(diǎn)實(shí)驗(yàn)室,北京 100094；2.北京四維智聯(lián)科技有限公司,北京 100094)

以MODIS遙感數(shù)據(jù)為驅(qū)動(dòng),采用以遙感像元為基本模擬單元的DPeRS(Diffuse pollution estimation with remote semsing)面源污染負(fù)荷估算模型,分析了2016年海河流域氮磷面源污染空間分布特征,并對(duì)“十三五”《重點(diǎn)流域水污染防治規(guī)劃》中劃定的海河流域172個(gè)控制單元進(jìn)行面源污染優(yōu)先控制單元分析.結(jié)果表明,2016年海河流域總氮(TN)排放量為13.62萬t,入河量為2.53萬t;總磷(TP)排放量為8152t,入河量為1597t;空間分布上,海河流域中部和南部地區(qū)氮磷面源污染較重,其中河北省片區(qū)氮磷面源污染物排放量及入河量最大;農(nóng)田徑流型是海河流域氮磷面源污染的主要類型,其次城鎮(zhèn)徑流影響也較大;篩選出海河流域TN和TP面源污染優(yōu)先控制單元分別為127和131個(gè),面積占比分別為84.2%和87.0%.

面源污染；遙感；氮磷；海河流域

海河流域是我國(guó)水資源開發(fā)利用程度最高、水污染最嚴(yán)重的流域,對(duì)其第二次水資源評(píng)價(jià)中首次涉及面源污染源調(diào)查,估算了2000年海河流域城鎮(zhèn)地表徑流、化肥農(nóng)藥使用、農(nóng)村生活污水及固體廢棄物、水土流失、分散式飼養(yǎng)畜禽廢水等5種面源污染源的產(chǎn)生量和入河量,得出海河流域TN和TP主要來自于面源的結(jié)論[1].面源污染作為影響飲用水源地水質(zhì)安全的重要污染,引起管理部門的廣泛關(guān)注.對(duì)于管理部門,面源污染存在污染負(fù)荷的量化和污染源的解析兩個(gè)急需解決的問題[2].模型是解決這兩個(gè)問題的有效手段,我國(guó)學(xué)者在面源污染負(fù)荷模型方面的研究較多,多采用SWAT過程機(jī)理模型[3-5]、農(nóng)業(yè)面源年化模型(AnnAGNPS)[6-7]、輸出系數(shù)模型(ECM)[1,8]和平均濃度法[9-11]等模型方法.SWAT及AnnAGNPS等主流面源模型存在所需數(shù)據(jù)量較大、建模過程較繁雜、模擬所需時(shí)間較長(zhǎng)等特點(diǎn),且多應(yīng)用于海河流域的部分區(qū)域[12-13];海河流域全流域的面源污染研究方法多集中于輸出系數(shù)法,已有對(duì)海河流域種植業(yè)、畜禽養(yǎng)殖業(yè)、農(nóng)村生活、農(nóng)業(yè)等方面的面源污染負(fù)荷量估算和污染特征分析[1,8,14-15],而應(yīng)用遙感技術(shù)對(duì)海河流域面源污染空間特征分析的相關(guān)研究較少.

綜上,本研究從管理需求出發(fā),構(gòu)建以遙感像元為基本水文單元的遙感分布式面源污染估算模型—DPeRS模型[16-17],對(duì)海河流域農(nóng)田生產(chǎn)、農(nóng)村生活、城鎮(zhèn)徑流、畜禽養(yǎng)殖和水土流失產(chǎn)生的TN和TP污染負(fù)荷進(jìn)行估算,分析海河流域氮磷面源污染排放量、入河量和空間分布特征,并進(jìn)行面源污染源和空間管理控制分析,為流域面源污染綜合防治部門提供決策支持.

1 材料和方法

1.1 研究區(qū)概況

海河流域是我國(guó)七大流域之一,位于我國(guó)華北地區(qū),介于東經(jīng)112°~120°、北緯35°~43°之間,東臨渤海,西倚太行,南界黃河,北接蒙古高原,涉及北京市、天津市、河北省大部、山西省東部、山東、河南2省北部,以及遼寧省和內(nèi)蒙古自治區(qū)的小部分,包括37個(gè)地市、337個(gè)縣(市、旗),流域面積為32萬km2,包括灤河、海河和徒駭馬頰河三大水系,研究區(qū)位置見圖1.

圖1 研究區(qū)位置示意

1.2 模型方法

1.2.1 DPeRS面源污染估算模型 DPeRS模型以大尺度面源污染模型[18]為基礎(chǔ),耦合遙感技術(shù),采用基于二元結(jié)構(gòu)原理構(gòu)建的遙感分布模型,實(shí)現(xiàn)了流域面源污染負(fù)荷的月尺度估算.DPeRS模型中既考慮了降水、植被蓋度、地形和地貌等自然要素,同時(shí)也考慮了施肥利用效率、人口、牲畜和家禽等社會(huì)經(jīng)濟(jì)要素.模型算法以遙感數(shù)據(jù)為驅(qū)動(dòng),耦合定量遙感模型和生態(tài)水文過程模型對(duì)流域尺度面源污染負(fù)荷的時(shí)空動(dòng)態(tài)進(jìn)行定量分析.DPeRS模型可以概括為5個(gè)污染類型和2個(gè)元素形態(tài),即:農(nóng)田徑流型、城鎮(zhèn)徑流型、農(nóng)村生活型、畜禽養(yǎng)殖型和水土流失型;溶解態(tài)和吸附態(tài)2種污染物形態(tài).具體污染指標(biāo)包括TN、TP、氨氮(NH4+-N)和化學(xué)需氧量(CODCr);模型包括5大模塊:農(nóng)田氮磷平衡核算模塊[19]、植被覆蓋度定量遙感反演模塊[2,20]、溶解態(tài)污染負(fù)荷估算模塊、吸附態(tài)污染負(fù)荷估算模塊和入河模塊;模型耦合的遙感數(shù)據(jù)包括植被覆蓋度和土地利用數(shù)據(jù),作為模型的輸入數(shù)據(jù),植被覆蓋度數(shù)據(jù)用于估算農(nóng)田徑流型和水土流失型氮磷面源污染負(fù)荷,土地利用數(shù)據(jù)在5個(gè)污染類型面源污染負(fù)荷估算中均有運(yùn)用,模型技術(shù)路線見圖2,模型核心算法詳見文獻(xiàn)[16,21].本研究中僅討論氮磷2個(gè)指標(biāo).

1.2.2 面源污染優(yōu)先控制單元篩選方法 十三五《重點(diǎn)流域水污染防治規(guī)劃》[22]中劃定全國(guó)1784個(gè)控制單元,其中海河流域占172個(gè).對(duì)海河流域控制單元進(jìn)行氮磷面源污染優(yōu)先控制單元的篩選,篩選方法具體為:將對(duì)照年(2005、2010、2015,3a的平均)每個(gè)省的兩項(xiàng)指標(biāo)平均值分別作為該省份面源污染優(yōu)先控制的閾值,本省份內(nèi)每個(gè)控制單元的2項(xiàng)指標(biāo)值分別與本省閾值進(jìn)行對(duì)比,面源污染排放負(fù)荷和入河量均大于本省閾值的單元判定為最優(yōu)先控制單元,命名為源頭、入河過程I類優(yōu)先控制單元;排放負(fù)荷大而入河量小或者排放負(fù)荷小而入河量大的控制單元判定為II類優(yōu)先控制單元,分別命名為源頭II類優(yōu)先控制單元和入河過程II類優(yōu)先控制單元,其他情況判定為一般面源污染控制單元.本研究中,對(duì)照年面源污染情況代表“十三五”前面源污染的本底狀況,海河流域TN和TP指標(biāo)的分省判定閾值見表1.

圖2 技術(shù)路線

1.3 數(shù)據(jù)來源與預(yù)處理

采用DPeRS模型模擬的月尺度面源污染估算海河流域年尺度面源污染量,重點(diǎn)監(jiān)測(cè)由降水引起的地表徑流型面源污染.DPeRS模型運(yùn)行需要的數(shù)據(jù)庫(kù)包括遙感數(shù)據(jù)、氣象數(shù)據(jù)、高程數(shù)據(jù)、土壤數(shù)據(jù)和農(nóng)業(yè)統(tǒng)計(jì)數(shù)據(jù)等,具體數(shù)據(jù)來源和處理方法見表2;具體帶入模型的輸入數(shù)據(jù)包括海河流域土地利用、月植被覆蓋度、月降水量、坡度坡長(zhǎng)、土壤類型等.圖3給出遙感數(shù)據(jù)解譯后得到的海河流域土地利用數(shù)據(jù)、定量遙感反演得到的月均植被覆蓋度數(shù)據(jù)、空間插值后的月降水累積數(shù)據(jù)和應(yīng)用ls_cal.aml程序計(jì)算得到的坡長(zhǎng)數(shù)據(jù).

表2 主要空間數(shù)據(jù)

2 結(jié)果與討論

2.1 海河流域農(nóng)田氮磷平衡空間特征

采用輸入輸出法估算海河流域農(nóng)田氮磷平衡量,結(jié)果表明,縣級(jí)核算結(jié)果氮平衡量為-8.65~ 99.60t/km2,磷平衡量為-2.33~54.03t/km2;結(jié)合空間分析技術(shù),將計(jì)算結(jié)果與土地利用數(shù)據(jù)進(jìn)行空間離散,得到2016年海河流域氮磷平衡的空間分布(圖4),結(jié)果表明,海河流域大部分地區(qū)表現(xiàn)為農(nóng)田氮磷營(yíng)養(yǎng)元素盈余狀態(tài),且磷素盈余較氮素更為明顯,流域中部和南部地區(qū)明顯高于北部地區(qū).

王雪蕾等[23]綜合分析了前人研究成果,將農(nóng)田氮磷平衡量劃分了3個(gè)風(fēng)險(xiǎn)等級(jí),氮素和磷素的低風(fēng)險(xiǎn)閾值分別為10,1.1t/km2,海河流域氮磷營(yíng)養(yǎng)平衡整體情況表明海河流域當(dāng)前的農(nóng)田管理方式不利于流域面源污染的防治.在推動(dòng)農(nóng)業(yè)生產(chǎn)發(fā)展的同時(shí),應(yīng)充分考慮區(qū)域特征,實(shí)施有針對(duì)性的管理措施,如提高化肥施用技術(shù)、合理利用有機(jī)肥、加大灌溉設(shè)施投入等,從源頭控制農(nóng)業(yè)面源污染,推動(dòng)農(nóng)業(yè)綠色發(fā)展.

圖4 海河流域農(nóng)田土壤氮磷平衡空間分布

2.2 海河流域氮磷面源污染物空間排放特征

表3 氮磷面源污染排放量統(tǒng)計(jì)

應(yīng)用DPeRS模型對(duì)海河流域2016年氮磷面源污染排放負(fù)荷進(jìn)行月尺度估算,結(jié)果表明:海河流域面源TN和TP平均污染排放負(fù)荷分別為0.43, 0.03t/km2,排放量分別為13.62萬t和8152.29t.海河流域氮磷面源污染排放負(fù)荷的空間分布如圖5所示,結(jié)果表明, 海河流域氮磷面源污染排放負(fù)荷的空間展布類似,流域中部和南部污染排放負(fù)荷較高,北部污染物負(fù)荷較低,具體表現(xiàn)為北京市南部、天津市、河北省片區(qū)中南部、山西省片區(qū)北部、河南和山東省片區(qū)等區(qū)域面源污染較重.氮磷面源污染空間統(tǒng)計(jì)分析結(jié)果表明,河南和山東省片區(qū)面源TN和TP污染排放負(fù)荷均高于流域平均水平,北京市、天津市和遼寧省片區(qū)面源TN污染排放負(fù)荷高于流域平均水平,河北、山西、河南和山東省片區(qū)面源TN和TP污染排放量相對(duì)較大,模擬結(jié)果的空間統(tǒng)計(jì)信息見表3.

圖5 海河流域氮磷面源污染排放負(fù)荷空間分布

2.3 海河流域氮磷面源污染物空間入河特征

表4 氮磷面源污染入河量統(tǒng)計(jì)

基于海河流域2016年氮磷面源污染排放負(fù)荷結(jié)果,結(jié)合當(dāng)年海河流域空間入河系數(shù),估算了2016年海河流域氮磷面源污染入河負(fù)荷,并統(tǒng)計(jì)分析了入河量信息.結(jié)果表明:海河流域面源TN和TP平均污染入河負(fù)荷分別為0.08,0.005t/km2,入河量分別為2.53萬t和1597t.海河流域氮磷面源污染入河負(fù)荷的空間分布如圖6所示,與排放負(fù)荷相比,海河流域入河負(fù)荷相對(duì)較小,這與該流域水資源量少有密切關(guān)系,海河流域氮、磷面源污染入河負(fù)荷的空間展布類似,流域中部和南部存在較高氮磷面源入河負(fù)荷的零星分布,山東省片區(qū)中部氮磷面源污染入河負(fù)荷較為突出.氮磷面源污染空間統(tǒng)計(jì)分析結(jié)果表明,天津市、河南和山東省片區(qū)面源TN和TP污染入河負(fù)荷均高于流域平均水平,北京市面源TN污染入河負(fù)荷高于流域平均水平,河北、山西、河南和山東省片區(qū)面源TN和TP污染入河量相對(duì)較大,模擬結(jié)果的空間統(tǒng)計(jì)信息詳見表4.

2.4 海河流域氮磷面源污染類型分析

對(duì)2016年海河流域農(nóng)田徑流型、城鎮(zhèn)徑流型、農(nóng)村生活型、畜禽養(yǎng)殖型和水土流失型5種氮磷面源污染類型進(jìn)行統(tǒng)計(jì)分析,結(jié)果表明:農(nóng)田徑流是海河流域最主要的氮磷面源污染源,對(duì)于TN,農(nóng)田生產(chǎn)引起的排放量和入河量均占總污染量的92%以上,其他類型順序如下:城鎮(zhèn)徑流型>農(nóng)村生活型>水土流失型>畜禽養(yǎng)殖型;對(duì)于TP,農(nóng)田生產(chǎn)引起的排放量和入河量均占總污染量的78%以上,其他類型順序如下:畜禽養(yǎng)殖型>城鎮(zhèn)徑流型>水土流失型>農(nóng)村生活型,不同類型面源污染排放量統(tǒng)計(jì)結(jié)果詳見表5.從面源污染類型角度分析,海河流域應(yīng)努力改進(jìn)田間施肥技術(shù),提高氮磷肥料利用率,降低其對(duì)水環(huán)境的危害,此外,城鎮(zhèn)徑流引起的氮磷面源污染也需引起重視.

圖6 海河流域氮磷面源污染入河負(fù)荷空間分布

表5 不同污染類型指標(biāo)核算

2.5 海河流域氮磷面源污染空間管控分析

海河流域氮磷面源污染優(yōu)先控制單元篩選結(jié)果表明:TN面源優(yōu)先控制單元共127個(gè),面積占比達(dá)到84.2%,其中,源頭、入河過程I類優(yōu)先控制單元共66個(gè),面積占比49.8%;II類優(yōu)先控制單元共61個(gè),包括19個(gè)源頭II類優(yōu)先控制單元和42個(gè)入河過程II類優(yōu)先控制單元,II類優(yōu)先控制單元面積占比達(dá)到34.4%.TP面源優(yōu)先控制單元共131個(gè),面積占比達(dá)到87.0%,其中,源頭、入河過程I類優(yōu)先控制單元共55個(gè),面積占比43.2%;II類優(yōu)先控制單元共76個(gè),包括30個(gè)源頭II類優(yōu)先控制單元和46個(gè)入河過程II類優(yōu)先控制單元,II類優(yōu)先控制單元面積占比達(dá)到43.8%.空間分布上,氮磷面源污染優(yōu)先控制單元主要分布在海河流域中部和南部地區(qū),源頭、入河過程I類優(yōu)先控制單元和入河過程II類優(yōu)先控制單元分布面積相對(duì)較大,2016年海河流域氮磷面源污染優(yōu)先控制單元空間分布詳見圖7.

2.6 DPeRS模型的適用性及結(jié)果對(duì)比

對(duì)比國(guó)內(nèi)外其他面源污染評(píng)估模型,DPeRS模型在模型結(jié)構(gòu)、運(yùn)行條件和模擬指標(biāo)等幾個(gè)方面具有較大的管理應(yīng)用優(yōu)勢(shì)[2].DPeRS模型以遙感像元為基本模擬單元,與SWAT等模型提出的水文響應(yīng)單元(HRU)相比,在保證模擬精度的前提下極大地提高了面源污染模擬的空間分辨率[16];同時(shí)DPeRS模型中耦合了定量遙感模型,彌補(bǔ)了無資料地區(qū)模型估算的不足;DPeRS模型的參數(shù)設(shè)置為開放模式,可以根據(jù)參數(shù)豐富度進(jìn)行重新構(gòu)架,根據(jù)管理需求完成不同尺度面源污染監(jiān)測(cè)與評(píng)估.此外,DPeRS模型系統(tǒng)可以實(shí)現(xiàn)遙感像元尺度的污染負(fù)荷空間可視化,直觀提供了面源污染空間分布的“關(guān)鍵地區(qū)”,與傳統(tǒng)總量減排核算方法相比,實(shí)現(xiàn)了從“點(diǎn)”到“面”的突破,為因地制宜制定面源污染防治方案提供了技術(shù)支撐[2].

圖7 海河流域氮磷面源污染優(yōu)先控制單元空間分布

采用地面監(jiān)測(cè)手段進(jìn)行面源污染監(jiān)測(cè),范圍小且類型有限,難度大,模型是較為合適的面源污染研究方法.面源污染評(píng)估模型方法多樣,所需數(shù)據(jù)繁雜,存在不確定性和區(qū)域適用性[24-25],即使在同一區(qū)域,研究結(jié)果也相差較大.朱梅等[1]采用輸出系數(shù)法對(duì)海河流域種植業(yè)面源污染負(fù)荷量進(jìn)行估算,2007年海河流域地表徑流TN流失量為51966t、TP為7976t,DPeRS模型估算的結(jié)果為2016年海河流域農(nóng)田徑流型TN排放量12.66萬t、TP排放量6574t;第一次全國(guó)污染源普查數(shù)據(jù)顯示,2007年北京和天津種植業(yè)TN流失量分別為8262.31t和9097.78t, TP流失量分別為387.99t和498.42t[26],DPeRS模型估算的結(jié)果表明2016年北京和天津農(nóng)田徑流型TN排放量分別為5951.13t和7370.66t,TP排放量分別為193.14t和287.06t.此外,時(shí)間尺度不同,也有一定影響.

3 結(jié)論

3.1 DPeRS模型在流域尺度面源污染模擬結(jié)果表明,海河流域2016年TN污染排放負(fù)荷平均為0.43t/km2,排放量為13.62萬t,入河負(fù)荷平均為0.08t/km2,入河量為2.53萬t;TP污染排放負(fù)荷平均為0.026t/km2,排放量為8152t,入河負(fù)荷平均為0.005t/km2,入河量為1597t.

3.2 氮磷面源污染物的空間分布特征表明,對(duì)海河流域水質(zhì)產(chǎn)生影響的面源污染物TN和TP主要集中在海河流域中部和南部地區(qū),河北省片區(qū)面源污染物排放量及入河量最大.

3.3 DPeRS面源污染源分析表明,農(nóng)田徑流是海河流域TN和TP最主要的面源污染源,分別占總污染量的92%和78%以上,且TN和TP污染負(fù)荷與農(nóng)田氮磷平衡相關(guān),此外,城鎮(zhèn)徑流引起的氮磷面源污染也需引起重視,畜禽養(yǎng)殖對(duì)TP面源污染也有一定影響.

3.4 氮磷面源污染空間管控分析結(jié)果表明,海河流域覆蓋了172個(gè)控制單元,TN面源優(yōu)先控制單元共127個(gè),面積占比達(dá)到84.2%;TP面源優(yōu)先控制單元共131個(gè),面積占比為87.0%.氮磷面源污染優(yōu)先控制單元主要分布在海河流域中部和南部地區(qū).

[1] 朱 梅,吳敬學(xué),張希三.海河流域種植業(yè)非點(diǎn)源污染負(fù)荷量估算[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2010,29(10):1907-1915. Zhu M, Wu J X, Zhang X S. Estimation on non-point source pollution loads of crop farming in Hai Basin [J]. Journal of Agro-Environment Science, 2010,29(10):1907-1915.

[2] 王雪蕾,吳傳慶,馮愛萍,等.利用DPeRS模型估算巢湖流域氨氮和化學(xué)需氧量的面源污染負(fù)荷[J]. 環(huán)境科學(xué)學(xué)報(bào), 2015,35(9):2883- 2891. Wang X L, Wu C Q, Feng A P, et al. Application of DPeRS model on estimation of non-point source pollution load of ammonia nitrogen and chemical oxygen demand in Chao Lake basin [J]. Acta Scientiae Circumstantiae, 2015,35(9):2883-2891.

[3] 歐陽(yáng)威,黃浩波,蔡冠清.巢湖地區(qū)無監(jiān)測(cè)資料小流域面源磷污染輸出負(fù)荷時(shí)空特征[J]. 環(huán)境科學(xué)學(xué)報(bào), 2014,34(4):1024-1031. Ouyang W, Huang H B, Cai G Q. Temporal and spatial characteristics of diffuse phosphorus pollution in the watershed without monitoring data at Chaohu Lake [J]. Acta Scientiae Circumstantiae, 34(4):1024- 1031.

[4] 秦耀民,胥彥玲,李懷恩.基于SWAT模型的黑河流域不同土地利用情景的非點(diǎn)源污染研究[J]. 環(huán)境科學(xué)學(xué)報(bào), 2009,29(2):440-448. Qin Y M, Xu Y L, Li H E.SWAT model of non-point source pollution under different land use scenarios in the Heihe river basin [J]. Acta Scientiae Circumstantiae, 2009,29(2):440-448.

[5] 羅 倩,任 理,彭文啟.遼寧太子河流域非點(diǎn)源氮磷負(fù)荷模擬分析 [J]. 中國(guó)環(huán)境科學(xué), 2014,34(1):178-186. Luo Q, Ren L, Peng W Q. Simulation study and analysis of non-point source nitrogen and phosphorus load in the Taizihe Watershed in Liaoning Province [J]. China Environmental Science, 2014,34(1):178- 186.

[6] 李開明,任秀文,黃國(guó)如,等.基于AnnAGNPS模型泗合水流域非點(diǎn)源污染模擬研究[J]. 中國(guó)環(huán)境科學(xué), 2013,33(S1):54-59. Li K M, Ren X W, Huang G R, et al. Simulation of non-point source pollution in Sihe watershed with Ann AGNPS [J]. China Environmental Science, 2013,33(S1):54-59.

[7] 邊金云,王飛兒,楊 佳,等.基于AnnAGNPS模型四嶺水庫(kù)小流域氮磷流失特征的模擬研究[J]. 環(huán)境科學(xué), 2012,33(8):2659-2666. Bian J Y, Wang F E, Yang J, et al. Simulation of nitrogen and phosphorus loss in Siling Reservoir Watershed with AnnAGNPS [J]. Environmental Science, 2013,33(S1):54-59.

[8] 邱 斌,李萍萍,鐘晨宇,等.海河流域農(nóng)村非點(diǎn)源污染現(xiàn)狀及空間特征分析[J]. 中國(guó)環(huán)境科學(xué), 2012,32(3):564-570. Qiu B, Li P P, Zhong C Y, et al. Characteristics and spatial distribution of the rural non-point source pollution in Haihe River Basin. [J]. China Environmental Science, 2012,32(3):564-570.

[9] 李懷恩.估算非點(diǎn)源污染負(fù)荷的平均濃度法及其應(yīng)用[J]. 環(huán)境科學(xué)學(xué)報(bào), 2000,20(4):397-400. Li H E. Mean centration method for estimation of nonpoint source load and its application [J]. Acta Scientiae circumstantial, 2000,20(4): 397-400.

[10] 歐陽(yáng)威,王 瑋,郝芳華,等.北京城區(qū)不同下墊面降雨徑流產(chǎn)污特征分析 [J]. 中國(guó)環(huán)境科學(xué), 2010,30(9):328-332. Ouyang Wei, Wang W, Hao F H, et al. Pollution characterization of urban stormwater runoff on different underlying surface conditions [J]. China Environmental Science, 2010,30(9):328-332.

[11] 李家科,李懷恩,董 雯,等.渭河關(guān)中段典型支流非點(diǎn)源污染監(jiān)測(cè)與負(fù)荷估算[J]. 環(huán)境科學(xué)學(xué)報(bào), 2011,31(7):1470-1478. Li J K, Li H E, Dong W, et al. Monitoring and load estimation of non-point source pollution on typical tributaries in the Guanzhong reach of the Weihe River [J]. Acta Scientiae Circumstantiae, 2011, 31(7):1470-1478.

[12] 朱新軍,王中根,李建新,等.SWAT模型在漳衛(wèi)河流域應(yīng)用研究[J]. 地理科學(xué)進(jìn)展, 2006,25(5):105-111. Zhu X J, Wang Z G, Li J X, et al. Applications of SWAT model in Zhang Wei River Basin [J]. Geography Science Progress, 2006,25(5): 105-111.

[13] 張永勇,王中根,于 磊,等. SWAT水質(zhì)模塊的擴(kuò)展及其在海河流域典型區(qū)的應(yīng)用[J]. 資源科學(xué), 2009,31(1):94-100. Zhang Y Y, Wang Z G, Yu L, et al. Extended Water Quality Moduel of SWAT Model And Its Application To Hai River Basin [J]. Resource Science, 2009,31(1):94-100.

[14] 朱 梅,吳敬學(xué),張希三.海河流域畜禽養(yǎng)殖污染負(fù)荷研究[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2010,29(8):1558-1565. Zhu M ,Wu J X, Zhang X S. Pollutants Loads of Livestock and Poultry Breeding in Hai Basin, China.[J]. Journal of Agro-Environment Science, 2010,29(8):1558-1565.

[15] 朱 梅.海河流域農(nóng)業(yè)非點(diǎn)源污染負(fù)荷估算與評(píng)價(jià)研究[D]. 北京:中國(guó)農(nóng)業(yè)科學(xué)院, 2011. Zhu M, Study on Agricultural NPSL Loads of Haihe Basin and Assessment on its Environmental Impact [D]. Beijing: Chinese Academy of Agricultural Sciences., 2011.

[16] Wang X L, Wang Q, Wu C Q, et al. A method coupled with remote sensing data to evaluate non-point source pollution in the Xin'anjiang catchment of China [J]. Science of the Total Environment, 2012,430: 132-143.

[17] 王雪蕾,蔡明勇,鐘部卿,等.遼河流域面源污染空間特征遙感解析[J]. 環(huán)境科學(xué), 2013,34(10):3788-3796. Wang X L, CAI M Y, Zhong B Q, et al. Research on spatial characteristic of non-point source pollution in Liaohe River Basin [J]. Environmental Science, 2013,34(10):3788-3796.

[18] 郝芳華,楊勝天,程紅光,等.大尺度區(qū)域面源污染負(fù)荷計(jì)算方法[J]. 環(huán)境科學(xué)學(xué)報(bào), 2006,26(3):375-383. Hao F H, Yang S T, Cheng H G, et al. A method for estimation of non-point source pollution load in the large-scale basins of China [J]. Acta Scientiae Circumstantiae, 2006,26(3):375-383.

[19] Wang X L, Feng A P, Wang Q, et al. Spatial variability of the nutrient balance and related NPSP risk analysis for agro-ecosystems in China in 2010 [J]. Agriculture, Ecosystems and Environment, 2014,193:42- 52.

[20] Jiménez-Mu?oz J C, Sobrino J A, Plaza A, et al. Comparison Between Fractional Vegetation Cover Retrievals from Vegetation Indices and Spectral Mixture Analysis: Case Study of PROBA/CHRIS Data Over an Agricultural Area [J]. Sensors, 2009,9(2):768-793.

[21] 王雪蕾,王新新,朱 利,等.巢湖流域氮磷面源污染與水華空間分布遙感解析[J]. 中國(guó)環(huán)境科學(xué), 2015,35(5):1511-1519. Wang X L, Wang X X, Zhu L, et al. Spatial analysis on diffuse pollution and algal bloom characteristic with remote sensing in Chao Lake Basin [J]. China Environmental Science, 2015,35(5):1511-1519.

[22] 環(huán)境保護(hù)部.重點(diǎn)流域水污染防治規(guī)劃 [Z]. 2017. Ministry of Environmental Protection. Water pollution control planning in key watershed [Z]. 2017.

[23] 王雪蕾,王 橋,吳傳慶,等.國(guó)家尺度面源污染業(yè)務(wù)評(píng)估與應(yīng)用示范[M]. 北京:科學(xué)出版社, 2015:10-115. Wang X L, Wang Q, Wu C Q, et al. Assesment and application of national diffuse pollution management [M]. Beijing: Science Press, 2015:10-115.

[24] 劉志紅,Li L T, McVicar T R,等.專用氣候數(shù)據(jù)空間插值軟件ANUSPLIN及其應(yīng)用[J]. 氣象, 2008,34(2):92-100. Liu Z H, Li L T, McVicar T R, et al. introduction of the professional interpolation software for meteorology data: ANUSPLIN [J]. Meteorological Monthly, 2008,34(2):92-100.

[25] 符素華,劉寶元,周貴云,等.坡長(zhǎng)坡度因子計(jì)算工具[J]. 中國(guó)水土保持科學(xué), 2015,13(5):105-110. Fu S H, Liu B Y, Zhou G Y, et al. Calculation tool of topographic factors [J]. Science of Soil and Water Conservation, 2015,13(5):105- 110.

[26] 第一次全國(guó)污染源普查資料編纂委員會(huì).第一次全國(guó)污染源普查資料文集(之四)-污染源普查技術(shù)報(bào)告(下冊(cè))[R].北京:中國(guó)環(huán)境科學(xué)出版社, 2011:8-11. Compilation Committee of data of the first China pollution source census. Anthology of the first China pollution source census (Part IV)- technical report on pollution source census (Volume 2) [R]. Beijing: China Environmental Press, 2011:8-11.

Spatial character analysis on nitrogen and phosphorus diffuse pollution in Haihe River Basin by remote sensing.

FENG Ai-ping1, WU Chuan-qing1, WANG Xue-lei1*, WANG Hong-liang2, ZHOU Ya-ming1, ZHAO Qian1

(1.Sate Environmental Protection Key Laboratory of Satellite Remote Sensing, Ministry of Ecology and Environment Center for Satellite Application on Ecology and Environment, Beijing 100094, China；2.Beijing AutoAi Technology Co.,Ltd., Beijing 100094, China)., 2019,39(7)：2999~3008

Driven by the MODIS data, DPeRS (Diffuse pollution estimation with remote sensing, DPeRS)model was used to analyses the spatial characteristic of nitrogen and phosphorus diffuse pollution in Haihe River basin on pixel scale in 2016. In Watershed Water Pollution Control Planning of the 13th Five-Year Plan published by MEE, PRC (Ministry Ecology and Environment of People’s Republic China), 172 control units were defined. Here the priority control units (PCU) in Haihe River Basin were analyzed. In 2016, the total discharge of total nitrogen (TN) was 2.53×104ton with the total production being 13.62×104ton, the total phosphorus (TP) was 1597 ton with the total production being 8152 ton. The nitrogen and phosphorus diffuse pollution was relatively serious in the central and southern areas of Haihe River Basin, and the max amount of production and discharge of pollutants were appeared in Hebei province. Agriculture was the most important source for nitrogen and phosphorus diffuse pollution, followed by urban runoff. The numbers of PCU to TN and TP were 127 and 131 respectively, which covered the areas of Haihe River Basin 84.3% and 87.0% respectively.

diffuse pollution; remote sensing; nitrogen and phosphorus; Haihe River Basin

X522

A

1000-6923(2019)07-2999-10

馮愛萍(1988-),女,山西晉中人,工程師,碩士,主要從事流域尺度面源污染防治研究.發(fā)表論文8篇.

2018-12-25

國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2016YFD0800903);國(guó)家自然科學(xué)基金資助項(xiàng)目(41871346)

* 責(zé)任作者, 研究員, wxlbnu@163.com