基于輸出系數(shù)模型的南四湖流域非點(diǎn)源污染輸出風(fēng)險(xiǎn)評(píng)估

荊延德, 張華美, 孫笑笑

(1.山東省高校南四湖濕地生態(tài)與環(huán)境保護(hù)重點(diǎn)實(shí)驗(yàn)室, 山東 濟(jì)寧 273165; 2.曲阜師范大學(xué) 地理與旅游學(xué)院, 山東 日照 276826)

基于輸出系數(shù)模型的南四湖流域非點(diǎn)源污染輸出風(fēng)險(xiǎn)評(píng)估

荊延德1,2, 張華美1,2, 孫笑笑1,2

(1.山東省高校南四湖濕地生態(tài)與環(huán)境保護(hù)重點(diǎn)實(shí)驗(yàn)室, 山東 濟(jì)寧 273165; 2.曲阜師范大學(xué) 地理與旅游學(xué)院, 山東 日照 276826)

[目的] 探究非點(diǎn)源污染物N、P輸出風(fēng)險(xiǎn)在不同子流域、坡度等級(jí)和縣市區(qū)的分布特征,為南四湖流域不同風(fēng)險(xiǎn)區(qū)制定管理方案和控制對(duì)策提供科學(xué)依據(jù)。[方法] 以數(shù)字高程模型(digital elevation model, DEM)為基礎(chǔ)數(shù)據(jù),運(yùn)用ArcGIS軟件,提取南四湖流域的DEM,制作坡度等級(jí)圖,再運(yùn)用水文分析功能,提取水系圖、河流圖、并劃定子流域,結(jié)合2013年南四湖流域土地利用圖,通過運(yùn)用輸出風(fēng)險(xiǎn)模型分析不同的土地利用類型和坡度下氮(N),磷(P)污染的空間變化。[結(jié)果] 氮素的平均風(fēng)險(xiǎn)概率達(dá)到51.67%,磷素的平均風(fēng)險(xiǎn)概率達(dá)到9.14%,南四湖流域非點(diǎn)源污染輸出風(fēng)險(xiǎn)湖東小于湖西;隨坡度增大,高風(fēng)險(xiǎn)區(qū)面積減小,低風(fēng)險(xiǎn)區(qū)面積增加;就不同縣市區(qū)而言,濟(jì)寧市中區(qū)的N風(fēng)險(xiǎn)較小之外,其余縣市區(qū)N風(fēng)險(xiǎn)均較大,P風(fēng)險(xiǎn)較小。[結(jié)論] N是流域最主要的非點(diǎn)源污染物,非點(diǎn)源污染輸出風(fēng)險(xiǎn)大小與土地覆蓋和坡度有密切關(guān)系。

南四湖流域; 輸出系數(shù); 非點(diǎn)源污染; 風(fēng)險(xiǎn)評(píng)估

文獻(xiàn)參數(shù)： 荊延德, 張華美, 孫笑笑.基于輸出系數(shù)模型的南四湖流域非點(diǎn)源污染輸出風(fēng)險(xiǎn)評(píng)估[J].水土保持通報(bào)，2017,37(3)：270-274.DOI:10.13961/j.cnki.stbctb.20170401.001； Jing Yande, Zhang Huamei, Sun Xiaoxiao. 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.DOI:10.13961/j.cnki.stbctb.20170401.001

隨著南四湖流域工業(yè)化、城市化進(jìn)程加快，大量工業(yè)廢水、生活污水、農(nóng)業(yè)回水等未經(jīng)處理注入河流和湖泊中，使得湖泊污染嚴(yán)重[1]。其中農(nóng)業(yè)引起的非點(diǎn)源污染是目前水體污染最重大的問題之一[2]。與點(diǎn)源污染相比，非點(diǎn)源污染的發(fā)生具有隨機(jī)性、廣泛性、模糊性、滯后性等特點(diǎn)[3]，這些特點(diǎn)決定了非點(diǎn)源污染難以預(yù)測(cè)，難以完全真實(shí)計(jì)算[4-6]。輸出系數(shù)模型避開了非點(diǎn)源污染發(fā)生的復(fù)雜過程，不考慮中間過程和內(nèi)在機(jī)制，操作簡(jiǎn)單，所需參數(shù)少，易獲取資料[7]，又有一定的核算精度，被廣泛應(yīng)用于大尺度流域范圍上的非點(diǎn)源污染研究[8]。

輸出系數(shù)模型在中國(guó)不同區(qū)域非點(diǎn)源污染研究上得到廣泛關(guān)注和實(shí)踐應(yīng)用，已有關(guān)于三峽庫區(qū)大寧河流域[9]、拉薩河流域[10]、呼蘭河流域[8]等流域非點(diǎn)源污染輸出風(fēng)險(xiǎn)分析的研究報(bào)道，研究重點(diǎn)大多集中在輸出系數(shù)的確定、污染輸出風(fēng)險(xiǎn)的計(jì)算、空間變化特征、控制對(duì)策制定等方面。但是以上研究大多針對(duì)河流流域的非點(diǎn)源污染輸出風(fēng)險(xiǎn)評(píng)價(jià)研究，對(duì)以湖為中心的流域非點(diǎn)源污染輸出風(fēng)險(xiǎn)研究甚少。據(jù)調(diào)查，針對(duì)山東省南四湖流域的非點(diǎn)源污染方面的研究還比較少，運(yùn)用輸出系數(shù)模型并結(jié)合RS和GIS技術(shù)對(duì)南四湖流域進(jìn)行輸出風(fēng)險(xiǎn)評(píng)估更是未見報(bào)道。本研究擬以南四湖流域?yàn)檠芯繉?duì)象，基于輸出系數(shù)法，主要考慮耕地、林地、草地、水域、建設(shè)用地和未利用土地6種土地利用類型非點(diǎn)源污染，估算南四湖各子流域非點(diǎn)源N，P輸出風(fēng)險(xiǎn)概率，分析非點(diǎn)源污染輸出風(fēng)險(xiǎn)隨坡度和行政區(qū)的變化特征，研究結(jié)果將為南四湖流域社會(huì)經(jīng)濟(jì)發(fā)展與生態(tài)環(huán)境保護(hù)相協(xié)調(diào)提供決策依據(jù)和理論參考。

1 研究區(qū)概況

南四湖是微山湖、昭陽湖、獨(dú)山湖、南陽湖的總稱，其中微山湖面積最大，故又稱微山湖，因位于山東省西南部的濟(jì)寧以南而得名。南四湖為淮河流域第2大淡水湖，是魯南、蘇北地區(qū)的重要水源地，也是國(guó)家南水北調(diào)東線工程重要的輸水通道和調(diào)蓄湖泊。南四湖流域位于北緯34°41′—35°98′，東經(jīng)114°82′—117°83′，入湖河流53條，總流域面積3.17×104km2，地處暖溫帶、半濕潤(rùn)地區(qū)，屬暖溫帶大陸性冬夏季風(fēng)氣候[1],流域內(nèi)多年平均降水量700 mm左右，年平均氣溫14.2 ℃。流域內(nèi)耕地面積占84%，有較大比重，隨著流域內(nèi)工農(nóng)業(yè)的迅速發(fā)展，人為氮、磷輸入量隨之增加，使得南四湖水質(zhì)狀況受到不同程度的影響，非點(diǎn)源污染是導(dǎo)致南四湖水質(zhì)惡化和水體富營(yíng)養(yǎng)化的主要因素。

2 研究方法和數(shù)據(jù)處理

2.1 研究方法

2.1.1 輸出系數(shù)模型 20世紀(jì)70年代初期，北美地區(qū)首先提出了輸出系數(shù)模型，主要用來估算土地利用和湖泊富營(yíng)養(yǎng)之間的關(guān)系[11]，假設(shè)所有土地類型的輸出系數(shù)都一樣，導(dǎo)致結(jié)果和真實(shí)情況有很大誤差，因而限制了模型的應(yīng)用[12]。Johnes等[13-14]不斷對(duì)模型進(jìn)行改進(jìn)，使模型得到了進(jìn)一步完善和廣泛應(yīng)用。輸出系數(shù)模型公式的一般表達(dá)為：

式中： j——污染物類別，本研究中主要為N，P兩種污染物； i——流域內(nèi)土地利用類型，本研究中主要考慮耕地、林地、草地、水域、建設(shè)用地、未利用土地； Eij——污染物j在域內(nèi)第i種土地利用類型上的輸出系數(shù)。公式中Ai表示的含義不同，Lj表示的意義也隨之變化。當(dāng)Ai表示流域內(nèi)第i種土地利用類型面積時(shí)，Lj表示的是污染物j在第i種土地利用上的污染負(fù)荷量(kg/a)；當(dāng)Ai表示流域內(nèi)第i種土地利用類型的面積百分比時(shí)，Lj表示的是流域內(nèi)污染物j帶權(quán)重的平均輸出系數(shù)〔kg/(km2·a)〕，此用法就是輸出風(fēng)險(xiǎn)模型，用于計(jì)算非點(diǎn)源污染輸出風(fēng)險(xiǎn)概率。P——降水中攜帶的污染物的負(fù)荷量(kg/a)，由于本研究區(qū)范圍較大，且地表土地利用覆蓋變化小，以N為主的營(yíng)養(yǎng)物質(zhì)在雨水中含量較少，降雨發(fā)生后研究區(qū)各處降水量強(qiáng)度差別不大，綜上所述，不考慮降水對(duì)研究區(qū)非點(diǎn)源污染帶來的風(fēng)險(xiǎn)。

由于林地相對(duì)于其他土地利用類型而言受到人為影響的作用較小，對(duì)水土流失又有滯留效應(yīng)，一般認(rèn)為林地是氮，磷污染的匯集區(qū)，自然狀況下不會(huì)產(chǎn)生非點(diǎn)源污染輸出風(fēng)險(xiǎn)，以林地的最大輸出系數(shù)作為評(píng)價(jià)閥值，超過閥值的百分比即為風(fēng)險(xiǎn)概率。

2.1.2 輸出系數(shù)的確定 運(yùn)用輸出風(fēng)險(xiǎn)模型的關(guān)鍵在于確定輸出系數(shù)，本研究采用查閱文獻(xiàn)法，得到耕地、林地、草地、水域、建設(shè)用地、未利用土地的輸出系數(shù)。中國(guó)對(duì)此研究[15]集中在長(zhǎng)江上游、松花江、黑河、太湖流域等熱點(diǎn)地區(qū)，但是針對(duì)山東省內(nèi)各流域的研究甚少，因此在確定輸出系數(shù)時(shí)，要考慮到與研究區(qū)自然狀況相近的輸出系數(shù)。參考孟曉云等[16]在云蒙湖流域的研究成果和全國(guó)范圍輸出系數(shù)的平均值，進(jìn)一步確定南四湖流域不同土地利用類型的輸出系數(shù)(表1)。

表1 南四湖流域不同土地利用類型全氮和全磷輸出系數(shù)

2.2 數(shù)據(jù)來源與處理

本研究所需數(shù)據(jù)有南四湖流域的DEM數(shù)據(jù)、

2013年Landsat TM影像和行政區(qū)劃圖(如表2)。主要運(yùn)用的數(shù)據(jù)處理軟件為ENVI 5.1和ArcGIS 10.1軟件。運(yùn)用ENVI 5.1軟件對(duì)遙感圖像進(jìn)行校正、鑲嵌、裁剪及監(jiān)督分類等，從而提取土地利用數(shù)據(jù)，結(jié)合全國(guó)土地資源分類系統(tǒng)和南四湖流域土地利用現(xiàn)狀，根據(jù)自己的研究需要將土地利用類型重新劃分為耕地(83.9%)、草地(0.2%)、林地(5.9%)、水域(5.2%)、建設(shè)用地(4.7%)與未利用土地(0.1%)。運(yùn)用ArcGIS軟件對(duì)下載的DEM數(shù)據(jù)進(jìn)行拼接、裁剪處理得到研究范圍內(nèi)的DEM數(shù)據(jù)，并依據(jù)南四湖流域?qū)嶋H地形特征、土地利用精度和DEM數(shù)據(jù)精度，利用ArcGIS水文分析功能，將研究區(qū)DEM數(shù)據(jù)經(jīng)過洼地填充、流向計(jì)算、匯流累積量計(jì)算及河網(wǎng)提取，根據(jù)自己的研究目的劃分出183個(gè)子流域。

表2 數(shù)據(jù)來源及說明

3 結(jié)果與分析

3.1 非點(diǎn)源污染輸出風(fēng)險(xiǎn)空間分布特征

本研究結(jié)合南四湖流域土地利用類型圖和子流域劃分結(jié)果，計(jì)算出耕地、林地、草地、水域、建設(shè)用地和未利用土地6種土地利用類型在每個(gè)子流域的面積百分比，代入公式求出每個(gè)子流域的輸出風(fēng)險(xiǎn)值，得出2013年南四湖流域總氮，總磷輸出風(fēng)險(xiǎn)空間分布狀況，結(jié)果見圖1。

由圖1可知，南四湖流域總氮，總磷輸出風(fēng)險(xiǎn)在空間分布上有較大差異?？偟骄敵鲲L(fēng)險(xiǎn)概率達(dá)到0.52，總磷平均輸出風(fēng)險(xiǎn)概率僅0.09，氮是流域最主要的非點(diǎn)源污染物。

流域內(nèi)氮風(fēng)險(xiǎn)概率超過0.50的區(qū)域占總流域的77.38%，大部分地區(qū)是中等風(fēng)險(xiǎn)區(qū)，磷風(fēng)險(xiǎn)概率為0的地區(qū)占流域面積的10.95%，研究區(qū)內(nèi)所有子流域的磷輸出風(fēng)險(xiǎn)概率均小于0.2，對(duì)南四湖流域非點(diǎn)源污染貢獻(xiàn)較小。

圖1 南四湖流域非點(diǎn)源污染輸出風(fēng)險(xiǎn)空間分布

從圖1中還可以看到，南四湖流域非點(diǎn)源氮污染輸出風(fēng)險(xiǎn)湖東小于湖西，一方面，氮在環(huán)境中多以溶解態(tài)形式存在，且據(jù)調(diào)查，南四湖湖東入湖徑流總量(2.75×109m3)大于湖西徑流(2.15×108m3)，即使大量污染物從東面入湖，也不會(huì)大幅度增加湖東面污染輸出風(fēng)險(xiǎn)；另一方面，湖東面大部分流域是市區(qū)，耕地占總流域面積小于建設(shè)用地，市區(qū)氮，磷污染的來源主要是生產(chǎn)生活排放，城市污水排放進(jìn)入環(huán)境前經(jīng)過處理，使得污染物含量降低。

由于南四湖流域研究區(qū)氮和磷輸出風(fēng)險(xiǎn)范圍有較大差別，目前沒有明確統(tǒng)一的等級(jí)劃分標(biāo)準(zhǔn)，根據(jù)圖1中氮、磷風(fēng)險(xiǎn)劃分間隔將風(fēng)險(xiǎn)劃分5個(gè)等級(jí)，等級(jí)越高，風(fēng)險(xiǎn)值越大。各風(fēng)險(xiǎn)等級(jí)累計(jì)統(tǒng)計(jì)分析結(jié)果見圖2。結(jié)合圖1和圖2可以看出，氮風(fēng)險(xiǎn)主要集中在0.52～0.56，占到總研究區(qū)面積的65%；磷風(fēng)險(xiǎn)主要集中在0.10～0.20，占總研究區(qū)面積的62%。氮風(fēng)險(xiǎn)遠(yuǎn)高于磷風(fēng)險(xiǎn)，氮風(fēng)險(xiǎn)高于0.50的面積占到77%，磷風(fēng)險(xiǎn)沒有超過0.50?？傮w來說，南四湖流域非點(diǎn)源氮污染較嚴(yán)重。

圖2 南四湖流域非點(diǎn)源污染輸出風(fēng)險(xiǎn)累計(jì)面積比

3.2 基于坡度等級(jí)的流域非點(diǎn)源污染輸出風(fēng)險(xiǎn)特征

運(yùn)用ArcGIS的空間分析工具，從研究區(qū)DEM中提取坡度數(shù)據(jù)，并進(jìn)行坡度分級(jí)，再將坡度分級(jí)數(shù)據(jù)和子流域圖結(jié)合得出每個(gè)子流域的平均坡度，根據(jù)自己研究需要將研究區(qū)劃分為4個(gè)坡度等級(jí)分別為平地[0°～3°)、緩平地[3°～5°)、緩地[5°～8°)、斜坡[8°～15°)，分析南四湖流域非點(diǎn)源污染輸出風(fēng)險(xiǎn)隨坡度的變化特征，結(jié)果如表3，坡度對(duì)于氮，磷風(fēng)險(xiǎn)的空間分布有較大影響。隨著坡度增大，氮輸出風(fēng)險(xiǎn)在等級(jí)1，2，3的面積逐漸增大，風(fēng)險(xiǎn)等級(jí)4，5面積逐漸減??；磷輸出風(fēng)險(xiǎn)各等級(jí)面積隨坡度的變化與氮有類似規(guī)律。且由于在研究區(qū)范圍內(nèi)地形主要以平原為主，土地利用類型以耕地最多，外界氮，磷輸入較多。

表3 南四湖流域非點(diǎn)源污染輸出風(fēng)險(xiǎn)在不同坡級(jí)上的分布

3.3 基于行政區(qū)劃的流域非點(diǎn)源污染輸出風(fēng)險(xiǎn)特征

運(yùn)用ArcGIS軟件，結(jié)合研究區(qū)的行政區(qū)劃圖、土地利用圖，求出每個(gè)縣(區(qū))各土地利用類型面積百分比，代入輸出系數(shù)模型公式，求得每個(gè)縣(區(qū))的輸出風(fēng)險(xiǎn)，分析結(jié)果如圖3所示。由圖3可以看出，在不同縣(區(qū))內(nèi)氮，磷輸出風(fēng)險(xiǎn)有較大差別。每個(gè)縣區(qū)內(nèi)氮輸出風(fēng)險(xiǎn)均大于磷輸出風(fēng)險(xiǎn)。在研究區(qū)范圍內(nèi)的30各縣區(qū)中，濟(jì)寧市中區(qū)污染輸出風(fēng)險(xiǎn)概率最小，氮風(fēng)險(xiǎn)概率不足1%，無磷污染風(fēng)險(xiǎn)；除鄒城區(qū)、棗莊市中區(qū)、山亭區(qū)氮風(fēng)險(xiǎn)概率小于0.5之外，余下26個(gè)縣區(qū)的氮風(fēng)險(xiǎn)概率均0.5～0.6，面積占到整個(gè)研究區(qū)的90%；棗莊市中區(qū)、山亭區(qū)是無磷風(fēng)險(xiǎn)區(qū)，其余縣區(qū)磷風(fēng)險(xiǎn)概率在0.2以內(nèi)。但風(fēng)險(xiǎn)值的跨度與張立坤等[8]研究的呼蘭河流域各縣市區(qū)風(fēng)險(xiǎn)分布跨度有所不同，由于南四湖流域土地利用類型在各縣市區(qū)之間差異較小，使得輸出風(fēng)險(xiǎn)概率相對(duì)集中在一定范圍內(nèi)。

圖3 南四湖流域非點(diǎn)源污染輸出風(fēng)險(xiǎn)在不同縣區(qū)的分布

在本研究中研究的縣市區(qū)數(shù)量比國(guó)內(nèi)其他類似研究包含的縣市區(qū)多，即減少了偶然誤差并增加了地區(qū)之間對(duì)比度，又能從更廣的視角體現(xiàn)流域內(nèi)氮，磷污染的空間現(xiàn)狀。綜合來看，南四湖流域氮污染范圍比較廣，大部分縣市區(qū)已即將超過中等污染程度，不僅對(duì)南四湖的水質(zhì)有破壞作用，還對(duì)以南四湖流域?yàn)橹饕Y源的生產(chǎn)生活造成嚴(yán)重影響。

3.4 南四湖流域非點(diǎn)源污染分區(qū)管理與控制對(duì)策

為了對(duì)不同污染區(qū)域提出有效的管理措施，將整個(gè)南四湖流域劃分為3個(gè)不同等級(jí)的風(fēng)險(xiǎn)區(qū)：輕微風(fēng)險(xiǎn)區(qū)(0.00～0.20)，低風(fēng)險(xiǎn)區(qū)(0.20～0.40)，中風(fēng)險(xiǎn)區(qū)(0.40～0.60)。從以上研究了解到，磷輸出風(fēng)險(xiǎn)值主要集中在0.00～0.20，整個(gè)研究區(qū)都屬于磷輕微風(fēng)險(xiǎn)區(qū)；氮輸出風(fēng)險(xiǎn)值主要集中0.20～0.60，其中低風(fēng)險(xiǎn)區(qū)占研究區(qū)面積的2%，中風(fēng)險(xiǎn)區(qū)占98%，屬于中風(fēng)險(xiǎn)程度的區(qū)域面積較大。由于氮元素是南四湖流域最主要的非點(diǎn)源污染物，因此對(duì)氮污染的治理工作最為緊迫。以氮的輸出風(fēng)險(xiǎn)值為依據(jù)，劃分管理區(qū)，結(jié)果見圖4。

圖4 南四湖流域氮污染治理管理分區(qū)

研究區(qū)有2%的面積在低風(fēng)險(xiǎn)區(qū)，主要分布在研究區(qū)西南的邊緣地區(qū)。土地利用類型以耕地和林地為主，且比例相近，耕地會(huì)增加污染風(fēng)險(xiǎn)，由于離湖區(qū)較遠(yuǎn)，對(duì)湖水影響不大。除了對(duì)該區(qū)實(shí)施保護(hù)策略外，增加對(duì)氮污染源控制、減小污染范圍有著重要作用。中風(fēng)險(xiǎn)區(qū)在研究區(qū)分布范圍大，土地利用類型以耕地和城鎮(zhèn)為主，污染來源較多，農(nóng)藥、化肥過度使用及生活污水排放都會(huì)造成非點(diǎn)源污染。南四湖流域主要以小麥和玉米為糧食作物，在成熟過程中對(duì)農(nóng)藥化肥需求量大，僅靠實(shí)施保護(hù)措施已達(dá)不到控制污染的目的，應(yīng)限制農(nóng)藥使用計(jì)量、提高農(nóng)藥效果、控制化肥使用量，才能減少非點(diǎn)源污染排放。在實(shí)際控制管理中，需要將有限的治理資源投入到對(duì)水體危害可能性最大而范圍相對(duì)較小的區(qū)域，從而降低治理難度并提高治理成效[17]。

4 討論與結(jié)論

(1) 河流流量和土地利用類型對(duì)非點(diǎn)源污染有明顯影響。流量越大，非點(diǎn)源污染輸出風(fēng)險(xiǎn)越低，非點(diǎn)源污染越輕。耕地和草地所占比重越大的地區(qū)，非點(diǎn)源污染輸出風(fēng)險(xiǎn)越高，污染越嚴(yán)重。從而使得南四湖流域非點(diǎn)源污染湖東小于湖西，這與李爽等[18]學(xué)者的觀點(diǎn)并不沖突。本研究氮輸出風(fēng)險(xiǎn)值相對(duì)集中在0.38～0.56，不如其他學(xué)者研究的河流流域氮風(fēng)險(xiǎn)值分散，原因是研究區(qū)范圍較大，地區(qū)之間差距不明顯，使得流域內(nèi)部未出現(xiàn)極端風(fēng)險(xiǎn)區(qū)。

(2) 坡度對(duì)南四湖流域非點(diǎn)源污染有顯著影響。隨著坡度增大，高風(fēng)險(xiǎn)區(qū)面積減小，低風(fēng)險(xiǎn)區(qū)面積增大。研究區(qū)內(nèi)坡度在0°～8°為最容易產(chǎn)生非點(diǎn)源污染的地區(qū)，研究結(jié)論與方廣玲等[10]研究的拉薩河流域低風(fēng)險(xiǎn)等級(jí)主要分布在15°以下自然用地的結(jié)論有所不同，主要是由于不同研究區(qū)其自然環(huán)境和土地利用狀況不同造成的。

(3) 南四湖流域非點(diǎn)源氮污染比磷污染嚴(yán)重，氮是最主要的非點(diǎn)源污染物，這與張立坤等[8]的研究結(jié)果一致。以氮的輸出風(fēng)險(xiǎn)值為依據(jù)，將南四湖流域劃分為不同的等級(jí)風(fēng)險(xiǎn)區(qū)，并針對(duì)不同風(fēng)險(xiǎn)區(qū)制定了相應(yīng)的控制對(duì)策。輕微風(fēng)險(xiǎn)區(qū)以保護(hù)為主，控制輸入到環(huán)境中氮，磷的量；低風(fēng)險(xiǎn)區(qū)主要污染物是氮，在研究區(qū)占有面積較小，因此可以通過種植具有固氮作用的植物，減小氮污染范圍；中風(fēng)險(xiǎn)區(qū)面積較廣，控制污染排放源達(dá)不到治理效果，改變土地利用方式、增加自然綠地面積能有效緩解污染狀態(tài)。

[1] 沈吉,張祖陸,楊麗原,等.南四湖:環(huán)境與資源研究[M].北京:地震出版社,2008.

[2] 許其功,劉鴻亮,沈珍瑤,等.三峽庫區(qū)典型小流域氮磷流失特征[J].環(huán)境科學(xué)學(xué)報(bào),2007,27(2):326-331.

[3] Weng Qihao. Modeling urban growth effects on surface runoff with the intergration of remote sensing and GIS[J]. Environmental Management, 2001,28(6):737-748.

[4] 丁曉雯,劉瑞民,沈真瑤.基于水文水質(zhì)資料的非點(diǎn)源輸出系數(shù)模型參數(shù)確定方法及其應(yīng)用[J].北京師范大學(xué)學(xué)報(bào):自然科學(xué)版,2006,42(5):534-538.

[5] 楊林章,王德建,夏立忠.太湖地區(qū)農(nóng)業(yè)面源污染特征及控制途徑[J].中國(guó)水利,2004(20):29-30.

[6] 蔡明,李懷恩,莊詠濤,等.改進(jìn)的輸出系數(shù)法在流域非點(diǎn)源污染負(fù)荷估算中的應(yīng)用[J].水利學(xué)報(bào),2004(7):1-8.

[7] Ding Xiaowen, Sheng Zhenyao, Hong Qian, et al. Development and test of the export coefficient model in the upper reach of the Yangtze River[J]. Journal of Hydrology, 2010,383(3/4):233.

[8] 張立坤,香寶,胡鈺,等.基于輸出系數(shù)模型的呼蘭河流域非點(diǎn)源污染輸出風(fēng)險(xiǎn)分析[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2014,33(1):148-154.

[9] 田甜,劉瑞民,王秀娟,等.三峽庫區(qū)大寧河流域非點(diǎn)源污染輸出風(fēng)險(xiǎn)分析[J].環(huán)境科學(xué)與技術(shù),2011,34(6):185-190.

[10] 方廣玲,香寶,杜加強(qiáng),等.拉薩河流域非點(diǎn)源污染輸出風(fēng)險(xiǎn)評(píng)估[J].農(nóng)業(yè)工程學(xué)報(bào),2015,31(1):247-254.

[11] Reckhow K H, Simpson J J. A procedure using modeling and error analysis for prediction of lake phosphorus concentration from land use information[J]. Canadian Journal of Fisheries and Aquatic Sciences, 1980,37(9):1439-1448.

[12] 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 modeling approach[J]. Journal of hydrology, 1996, 183(3/4):323-349.

[13] Johnes P J, O′Sullivan P E. The natural history of Slapton Ley nature research ⅩⅤⅢ. Nitrogen and phosphorus losses from the catchment: An export coefficient approach[J]. Field Studies, 1989,7(2):285-309.

[14] Johnes P J. An investigation of the effects of land use upon water quality in the windrush catchment[D]. Oxford: University of Oxford, 1990.

[15] 應(yīng)蘭蘭,侯西勇,路曉,等.我國(guó)非點(diǎn)源污染研究中輸出系數(shù)問題[J].水資源與水工程學(xué)報(bào),2010,21(6):90-95.

[16] 孟曉云,于興修,泮雪芹.云蒙湖流域土地利用變化對(duì)非點(diǎn)源氮污染負(fù)荷的影響[J].環(huán)境科學(xué),2012,33(6):1789-1794.

[17] Strauss P, Leone A, Ripa M N, et al. Using critical source areas for targeting cost-effective best management practices to mitigate phosphorous and sediment transfer at the watershed scale[J]. Soil Use and Management, 2007,23(S1):144-153.

[18] 李爽,張祖陸,孫媛媛.基于SWAT模型的南四湖流域非點(diǎn)源氮磷污染模擬[J].湖泊科學(xué),2013,25(2):236-242.

Risk Assessment of Non-point Source Pollution in Nansihu Lake Basin Using Output Coefficient Model

JING Yande1,2, ZHANG Huamei1,2, SUN Xiaoxiao1,2

(1.KeyLaboratoryofNansihuLakeWetlandEcologicalandEnvironmentalProtectioninUniversitiesofShandong,Jining,Shandong273165,China; 2.CollegeofGeographyandTourism,QufuNormalUniversity,Rizhao,Shandong276826,China)

[Objective] The objective of this study is to explore the distribution characteristics of the output risk of non-point source pollutants nitrogen(N) and phosphorus(P) under different conditions of sub basin, slope grade and geographical demarcation(counties), and to provide scientific basics for administer plan and control countermeasure of different risk areas of Nansihu Lake Basin. [Methods] Based on digital elevation model(DEM) data, the Nansihu Lake Basin DEM was extracted and slope grade was extracted using ArcGIS software. Then applying hydrologic analysis function, we extracted the drainage and river figures and delimited the drainage basin. Combined with the land use map of Nansihu Lake Basin in 2013, the space changes of N, P pollution under different land cover and slope were analyzed using the output risk model. [Results] The average risk probability of N reached 51.67%. The average risk probability of P reached 9.14%. The output risk of non-point source pollution in the east of lake was less than that in the west. With the increase of slope, the high risk area reduced, but the low risk area increased. In terms of geographical demarcation, Jining City was determined as less risk area of N, other counties had higher risk. The risk of P was less for all counties. [Conclusion] Nitrogen is the main basin non-point source pollutants. The output risk value of non-point source pollution has close relationship with land cover and slope.

Nansihu Lake Basin; output coefficient; non-point source pollution; risk assessment

2016-10-09

2016-10-20

教育部人文社會(huì)科學(xué)研究面上項(xiàng)目“基于輸出系數(shù)模型的南四湖流域非點(diǎn)源污染負(fù)荷估算及其對(duì)土地利用變化的響應(yīng)研究”(15YJAZH027)

荊延德(1970—),男(漢族),山東省淄博市人,博士,教授,碩士生導(dǎo)師,主要從事環(huán)境修復(fù)與資源再生方面的研究。E-mail:jingyande@163.com。

張華美(1989—),女(漢族),河南省項(xiàng)城市人,碩士研究生,研究方向?yàn)榄h(huán)境保護(hù)與修復(fù)。E-mail:zhanghuamei12@163.com。

A

1000-288X(2017)03-0270-05

X820.4