呂 川，張洪鈺，齊 琳

(吉林省環(huán)境科學(xué)研究院，吉林 長(zhǎng)春 130012)

研究表明，目前點(diǎn)源逐漸得到有效控制，非點(diǎn)源污染成為繼點(diǎn)源污染后水環(huán)境的首要污染源［1］。其中農(nóng)業(yè)引起的非點(diǎn)源污染是目前水體污染最重大的問(wèn)題之一［2］。農(nóng)業(yè)生產(chǎn)活動(dòng)中通過(guò)農(nóng)田地表徑流、農(nóng)田排水和地下滲漏，以農(nóng)田中的氮、磷、農(nóng)藥、化肥及其他有機(jī)或無(wú)機(jī)污染物質(zhì)(通過(guò)降水或灌溉等形式)進(jìn)入水體所造成的污染，統(tǒng)稱為農(nóng)業(yè)非點(diǎn)源污染(Nonpoint－source agricultural pollution)［3］。非點(diǎn)源污染的發(fā)生具有隨機(jī)性、廣泛性、模糊性和滯后性等特點(diǎn)［4］，其影響因素也比較廣泛，如土壤、地形、降水、水文條件、土地利用方式等，污染負(fù)荷分布差異性非常顯著，體現(xiàn)在不同景觀單元上差異也十分明顯［5－6］。決定收納水體質(zhì)量的少數(shù)景觀單元輸出的污染物占整個(gè)流域污染負(fù)荷較大比重［7］，因而成為非點(diǎn)源污染物的關(guān)鍵源區(qū)(CSAs:Critical Source Areas)。確定流域內(nèi)農(nóng)業(yè)非點(diǎn)源污染最有效的方法就是針對(duì)不同的農(nóng)田生態(tài)系統(tǒng)進(jìn)行非點(diǎn)源污染風(fēng)險(xiǎn)評(píng)價(jià)［8］。國(guó)外有不少學(xué)者認(rèn)識(shí)到了這一點(diǎn)，進(jìn)行了大量相關(guān)的研究［9－10］。Petersen等［9］對(duì)賓夕法尼亞州進(jìn)行研究采用GIS技術(shù)和農(nóng)業(yè)非點(diǎn)源污染潛力指數(shù)系統(tǒng)(APPI:Agricultural Pollution Potentia Index)來(lái)評(píng)價(jià)非點(diǎn)源污染發(fā)生的潛力，識(shí)別不同流域污染流失的重點(diǎn)發(fā)生區(qū)。Lemunyon等［11］對(duì)農(nóng)業(yè)地區(qū)磷元素流失的危險(xiǎn)性進(jìn)行了評(píng)價(jià)，在多因子相互作用的影響下使用了磷指數(shù)評(píng)價(jià)法。Hughes等［12］對(duì)愛(ài)爾蘭的農(nóng)業(yè)地區(qū)進(jìn)行磷分級(jí)方案存在的研究，提出了地塊尺度磷流失危險(xiǎn)分級(jí)方案和流域尺度的磷流失危險(xiǎn)分級(jí)方案。李琪等［13］對(duì)官?gòu)d水庫(kù)上游妨水河農(nóng)耕區(qū)的磷元素風(fēng)險(xiǎn)性進(jìn)行了評(píng)價(jià)，以流域磷分級(jí)方案為前提，提出了修正的流域尺度磷分級(jí)方案。

本文綜合考慮遼河源頭區(qū)流域影響非點(diǎn)源污染物流失的主要因素，綜合多種數(shù)據(jù)源，采用RS和GIS技術(shù)，篩選非點(diǎn)源污染風(fēng)險(xiǎn)的主要影響因子，從源因子和遷移因子兩方面出發(fā)建立遼河源頭區(qū)流域風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)體系，得出污染流失的風(fēng)險(xiǎn)指數(shù)。根據(jù)風(fēng)險(xiǎn)指數(shù)對(duì)農(nóng)業(yè)非點(diǎn)源污染風(fēng)險(xiǎn)劃分等級(jí)，識(shí)別非點(diǎn)源污染的關(guān)鍵源區(qū)，針對(duì)關(guān)鍵源區(qū)發(fā)生風(fēng)險(xiǎn)的順序提出管理措施和建議，為流域非點(diǎn)源污染控制工作提供依據(jù)，對(duì)流域治理目標(biāo)的實(shí)現(xiàn)起到積極作用。

1 研究區(qū)概況

遼河源頭區(qū)流域位于吉林省西南部，地處東經(jīng) 123°42'－125°31'，北緯 42°34'－44°08'，流域面積15746 km2，占全省總土地面積的8.4%，主要河流有東遼河、西遼河、招蘇臺(tái)河。流域地處松遼平原中部，低山、丘陵和平原兼?zhèn)?，位于遼河流域上游，地勢(shì)由東南向西北緩降，海拔高程611～120 m。東南部為低山丘陵地帶，山高谷寬，山間夾雜小的溝川平地，土質(zhì)肥沃，山體表層土質(zhì)較好。中西部為平原區(qū)，為起伏的臺(tái)地和平緩的平原，東遼河、招蘇臺(tái)河穿行其間，沿河兩岸地勢(shì)低平。氣候?qū)贉貛Т箨懶约撅L(fēng)氣候區(qū)，四季分明，多年平均降水量為545mm，大部分降水量集中在6～9月，約占全年降水量的80%，多年平均蒸發(fā)量為1020mm。地表水資源總量為15.02億m3，人均水資源量為619 m3，是全省人均占有量的45%，水資源相對(duì)短缺。水資源年內(nèi)變化較大，大部分降水量集中在6～9月，約占全年降水量的80%。徑流年內(nèi)分配不均衡，6～9月份為豐水期，3、4、5、10月份為平水期，11月份進(jìn)入枯水期，流量顯著減少，有的河段甚至出現(xiàn)斷流現(xiàn)象。

遼河源頭區(qū)流域地跨四平市、遼源市所轄的7個(gè)市、縣，共計(jì)81個(gè)鄉(xiāng)鎮(zhèn)。分別是:四平市市區(qū)、公主嶺市、梨樹(shù)縣、雙遼市、伊通縣(大孤山鎮(zhèn)、小孤山鎮(zhèn)、靠山鎮(zhèn))，遼源市市區(qū)、東遼縣。2008年，流域所轄7個(gè)市縣總?cè)丝?52.5萬(wàn)人，其中農(nóng)業(yè)人口196.32萬(wàn)人［14］。該流域是吉林省重要的能源基地，同時(shí)也是東北地區(qū)重要的商品糧基地。流域內(nèi)農(nóng)業(yè)產(chǎn)業(yè)化進(jìn)程較快，是吉林省農(nóng)業(yè)種植和牧業(yè)養(yǎng)殖的重點(diǎn)區(qū)域。2008年，流域內(nèi)耕地約占流域面積的42.19%。豐富的糧食資源為發(fā)展畜禽養(yǎng)殖提供了有利條件，2008年，流域內(nèi)共有大牲畜近119.2萬(wàn)頭，小牲畜近470.5萬(wàn)頭。但近年來(lái)，由于流域內(nèi)產(chǎn)業(yè)結(jié)構(gòu)不盡合理，工業(yè)結(jié)構(gòu)性污染和面源比較突出，對(duì)遼河流域水質(zhì)造成較嚴(yán)重污染，使生態(tài)環(huán)境與社會(huì)經(jīng)濟(jì)發(fā)展的矛盾日益尖銳。

2 指標(biāo)體系的建立

2.1 指標(biāo)的選擇

非點(diǎn)源污染的發(fā)生受多種因素的影響:不同的土地利用類(lèi)型影響著水土保持的程度及氮磷元素的來(lái)源量，施肥量及牲畜排泄量直接影響氮磷元素的排放程度，土地距河流的遠(yuǎn)近程度影響非點(diǎn)源污染的遷移擴(kuò)散過(guò)程;非點(diǎn)源污染的發(fā)生還受降雨影響，降雨具有間歇性，其強(qiáng)度大小又受到發(fā)生地的土壤類(lèi)型、土地利用類(lèi)型和地形條件的約束，從而具有顯著的區(qū)域特征，此外，人口的變動(dòng)也直接作用于非點(diǎn)源污染。因此，結(jié)合現(xiàn)有資料，最終選擇遼河源頭區(qū)流域的人口密度因子、年降雨量因子、土地利用類(lèi)型因子、至河流的距離因子、氮磷肥施用量因子、畜禽養(yǎng)殖污染物因子和年徑流深因子7個(gè)影響因子構(gòu)建了農(nóng)業(yè)非點(diǎn)源污染風(fēng)險(xiǎn)評(píng)估指標(biāo)體系。

2.2 指標(biāo)權(quán)重的確定

各影響因子對(duì)農(nóng)業(yè)非點(diǎn)源污染的危害程度不同，因此需要確定各影響因子和權(quán)重以獲得更準(zhǔn)確的污染風(fēng)險(xiǎn)等級(jí)［15］。此次采用層次分析法中的冪法確定影響因子權(quán)重。由于農(nóng)業(yè)非點(diǎn)源污染受多因子共同作用，且具有隨機(jī)性、廣泛性、模糊性和滯后性等特點(diǎn)，因此適合采用層次分析法［16］。分析得到各評(píng)價(jià)因子的權(quán)重(表1)。

表1 遼河源頭區(qū)流域農(nóng)業(yè)非點(diǎn)源污染風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)及權(quán)重Tab.1 Agricultural non-point source pollution risk evaluation and weight in the source areas of Liaohe river watershed

3 非點(diǎn)源風(fēng)險(xiǎn)評(píng)價(jià)模型的建立

3.2 各影響因子的獲取

3.2.1 人口密度因子 人口與非點(diǎn)源污染的關(guān)系十分密切，人口數(shù)量的增長(zhǎng)直接引起了污染物排放量的增多，空氣污染以及水質(zhì)污染等，是影響非點(diǎn)源污染的重要因素之一。本文通過(guò)四平、遼源2009年年鑒得到遼河源頭區(qū)各市縣人口密度，采用反距離空間差值法得到遼河源頭區(qū)人口密度圖(圖1)。

3.2.2 年降雨量因子 降雨是影響地表土壤侵蝕和非點(diǎn)源擴(kuò)散的重要因素之一，而遼河源頭區(qū)流域內(nèi)降雨年際變化大，大部分降雨集中在夏季，對(duì)河流流量變化影響顯著。本文采用研究區(qū)遼源氣象站和梨樹(shù)氣象站的年降雨量數(shù)據(jù)，使用ArcGIS空間分析模塊的插值功能將點(diǎn)數(shù)據(jù)形成趨勢(shì)面，得到遼河源頭區(qū)流域的年降雨量圖層(圖2)。

3.2.3 至河流的距離因子 不同土地距河流的遠(yuǎn)近程度對(duì)非點(diǎn)源污染風(fēng)險(xiǎn)有著相當(dāng)?shù)挠绊?，?dāng)評(píng)價(jià)尺度不斷擴(kuò)大到流域尺度時(shí)，地表徑流和土壤侵蝕為主產(chǎn)生的非點(diǎn)源污染會(huì)隨著匯入水體而造成更大的風(fēng)險(xiǎn)，風(fēng)險(xiǎn)隨著距河流的距離增加而減少，因此距河流距離成為制約非點(diǎn)源污染擴(kuò)散過(guò)程的一個(gè)重要因子［17－18］?；贒EM生成的水系圖，利用ArcGIS軟件的距離制圖功能，得到距離因子圖(圖3)。

圖1 遼河源頭區(qū)流域人口密度因子Fig.1 Population density factor in the source areas of Liaohe river watershed

圖2 遼河源頭區(qū)流域年降雨量因子Fig.2 Annual rainfall factor in the source areas of Liaohe river watershed

圖3 遼河源頭區(qū)流域至河流距離因子Fig.3 To the river distance factor in the source areas of Liaohe river watershed

圖4 遼河源頭區(qū)流域禽畜污染物因子Fig.4 Livestock contaminants factor in the source areas of Liaohe river watershed

3.2.4 畜禽養(yǎng)殖污染物因子 近年來(lái)，隨著四平和遼源地區(qū)經(jīng)濟(jì)的發(fā)展，遼河源頭區(qū)的規(guī)?；B(yǎng)殖場(chǎng)和散戶養(yǎng)殖數(shù)量與污染治理措施不相匹配，畜禽的糞便和污水排放量劇增，使得畜禽養(yǎng)殖業(yè)非點(diǎn)源污染問(wèn)題日益凸顯。根據(jù)2009年四平及遼源統(tǒng)計(jì)年鑒，得到遼河源頭區(qū)流域各地區(qū)畜禽養(yǎng)殖業(yè)情況，見(jiàn)表2。

表2 遼河源頭區(qū)流域各地區(qū)畜禽養(yǎng)殖量Tab.2 Various regions of livestock and poultry in the source areas of Liaohe river watershed

結(jié)合《畜禽養(yǎng)殖業(yè)污染物排放標(biāo)準(zhǔn)》(GB 18596－2001)中的排泄系數(shù)和畜禽糞便中污染物的含量，確定2009年遼河源頭區(qū)各縣市氮磷產(chǎn)生量，并按照《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB－3838－2002)中3類(lèi)標(biāo)準(zhǔn)求出等標(biāo)污染負(fù)荷，折算得到遼河源頭區(qū)的禽畜污染物密度見(jiàn)表3。利用ArcGIS軟件得到該因子圖層(圖4)。

表3 遼河源頭區(qū)流域各縣市禽畜污染物密度Tab.3 Livestock pollutant density of cities and counties in the source areas of Liaohe river watershed

3.2.5 氮磷肥施用量因子 施肥量是影響農(nóng)業(yè)非點(diǎn)源污染的重要因子，適量的施用化肥可以促進(jìn)農(nóng)作物生長(zhǎng)，但超量施用化肥土壤會(huì)養(yǎng)分過(guò)量富集，同時(shí)化肥隨水土流失及地面徑流匯入水體，造成水體富營(yíng)養(yǎng)化。對(duì)遼河源頭區(qū)各縣市采用施肥量平均計(jì)算的方法，將獲取的各縣市氮磷肥使用量中氮磷有效成分除以該地區(qū)面積，作為各縣市的單位面積氮磷含量，得到氮磷肥施用量因子(圖5)。

3.2.6 土地利用類(lèi)型因子 土地利用類(lèi)型影響土壤保持水土的能力和氮磷元素的來(lái)源量［19］，土地利用方式的不同間接決定非點(diǎn)源污染的程度。土地利用類(lèi)型因子采用2009年美國(guó)陸地資源衛(wèi)星Landsat TM影像進(jìn)行遙感解譯，劃分為林地、草地、耕地、交通、住宅、水域以及未利用地。河源頭區(qū)的土地利用類(lèi)型因子見(jiàn)圖6。

3.2.7 年徑流深因子 地表徑流是污染物輸入水體的主要?jiǎng)恿C(jī)制，也是衡量地表徑流量大小的重要依據(jù)。徑流深是指將年徑流量均勻地鋪在整個(gè)流域面積上所相當(dāng)?shù)乃畬由疃龋?0］。本文采用遼河源頭區(qū)流域8個(gè)水文站1999—2009年的日流量數(shù)據(jù)為基礎(chǔ)，計(jì)算得到每個(gè)水文站的11年的年徑流量數(shù)據(jù)，在ArcGIS下插值計(jì)算生成流域年徑流量柵格數(shù)據(jù)。多年平均徑流量與流域面積的比值即為多年平均年徑流深，通過(guò)GIS柵格計(jì)算獲得各個(gè)柵格的年徑流深分布圖(圖7)。

圖5 遼河源頭區(qū)流域氮磷肥施用量因子Fig.5 The amount of nitrogen and phosphorus application factor in the source areas of Liaohe river watershed

圖6 遼河源頭區(qū)流域土地利用類(lèi)型因子Fig.6 Type of land use factor in the source areas of Liaohe river watershed

3.3 非點(diǎn)源污染風(fēng)險(xiǎn)綜合指數(shù)的計(jì)算

對(duì)多種數(shù)據(jù)源進(jìn)行預(yù)處理后，利用Arc/Info的Grid模塊，以30 m×30 m柵格為研究單元，對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，將7個(gè)因子的柵格屬性數(shù)據(jù)進(jìn)行空間疊加得到綜合指數(shù)柵格數(shù)據(jù)文件和屬性表，采用多因子加權(quán)綜合評(píng)分法，計(jì)算每個(gè)柵格的農(nóng)業(yè)非點(diǎn)源風(fēng)險(xiǎn)綜合指數(shù)。

農(nóng)業(yè)非點(diǎn)源風(fēng)險(xiǎn)綜合指數(shù)表示為:

式中:wi為各因子對(duì)農(nóng)業(yè)面源污染風(fēng)險(xiǎn)程度的影響權(quán)重，Pi為各因子標(biāo)準(zhǔn)化后的等級(jí)分值。

分級(jí)標(biāo)準(zhǔn)的確定借鑒了前人的研究成果［21－22］，為反映農(nóng)業(yè)非點(diǎn)源污染風(fēng)險(xiǎn)的貢獻(xiàn)程度，根據(jù)研究區(qū)實(shí)際情況對(duì)各影響因子的分級(jí)標(biāo)準(zhǔn)進(jìn)行校正，最終將7個(gè)因子劃分為4個(gè)級(jí)別(表4)，分布賦予1－4分，不同級(jí)別對(duì)農(nóng)業(yè)非點(diǎn)源污染風(fēng)險(xiǎn)程度貢獻(xiàn)不同，級(jí)別越高貢獻(xiàn)值越大。

表4 農(nóng)業(yè)非點(diǎn)源污染風(fēng)險(xiǎn)影響因子分級(jí)Tab.4 Impact factor classification of agricultural non-point source pollution risk

圖8 遼河源頭區(qū)流域非點(diǎn)源污染風(fēng)險(xiǎn)程度等級(jí)分布圖Fig.8 The degree of risk of non－point source pollution level map in the source areas of Liaohe river watershed

4 結(jié)果與分析

4.1 非點(diǎn)源污染風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果

根據(jù)農(nóng)業(yè)非點(diǎn)源風(fēng)險(xiǎn)綜合指數(shù)公式，通過(guò)GIS的空間分析運(yùn)算功能將各因子的評(píng)價(jià)值圖件處理疊加，得到研究區(qū)內(nèi)農(nóng)業(yè)非點(diǎn)源風(fēng)險(xiǎn)指數(shù)圖，再根據(jù)分級(jí)標(biāo)準(zhǔn)將圖件重分類(lèi)得到流域尺度的農(nóng)業(yè)非點(diǎn)源風(fēng)險(xiǎn)評(píng)價(jià)圖(圖8)。

遼河源頭區(qū)非點(diǎn)源污染發(fā)生的四級(jí)區(qū)即關(guān)鍵源區(qū)占流域面積的21.77%，主要集中在梨樹(shù)南部、公主嶺東部和南部地區(qū)。四級(jí)區(qū)主要分布在位于遼河源頭區(qū)東北部、中部以及東南部丘陵區(qū)，人口密度、畜禽養(yǎng)殖污染物以及施肥因子高風(fēng)險(xiǎn)等級(jí)區(qū)內(nèi)，且是坡耕地的主要分布區(qū)域。三級(jí)區(qū)占流域總面積的37.45%，分布在除雙遼縣外遼河源頭區(qū)大部分地區(qū)。二級(jí)區(qū)占流域面積的30.03%，主要分布在梨樹(shù)縣、四平市的北部和公主嶺市的西南部。二級(jí)和三級(jí)區(qū)位于中部和南部平原丘陵區(qū)，土地利用以旱地為主;地形階梯狀結(jié)構(gòu)明顯，由大部分的平原和小部分的丘陵組成，因此地形因子影響較小，畜禽污染因子和施肥因子影響較大，農(nóng)田養(yǎng)分流失是本區(qū)的主要問(wèn)題。一級(jí)區(qū)僅占流域面積的10.75%，主要分布在雙遼市，地處遼河源頭區(qū)西部地，地勢(shì)平緩，土地利用主要為耕地、住宅用地和未利用地，植被覆蓋率較低。

4.2 管理控制措施

根據(jù)已經(jīng)劃定的遼河源頭區(qū)流域非點(diǎn)源風(fēng)險(xiǎn)等級(jí)區(qū)，本章根據(jù)不同等級(jí)區(qū)的特點(diǎn)和實(shí)際情況，分別制定不同面源污染控制區(qū)域面源污染控制管理控制措施。①畜禽養(yǎng)殖類(lèi)非點(diǎn)源污染控制措施。遼河源頭區(qū)非規(guī)?；酿B(yǎng)殖是控制禽畜養(yǎng)殖污染的關(guān)鍵，其中非規(guī)模化養(yǎng)殖占養(yǎng)殖總量約75%，污染物流失量很大?？蓪⑿笄菁S污制成有機(jī)肥，有效回收氮、磷、鉀等元素，減少畜禽排泄帶來(lái)的非點(diǎn)源污染。②土壤侵蝕類(lèi)非點(diǎn)源污染控制措施。針對(duì)流域土壤侵蝕嚴(yán)重的現(xiàn)狀，對(duì)于坡耕地采用水土保持及生物措施，由于地梗植物帶減緩水土流失和土壤中N、P流失效果較好，可采用3種措施改變土地利用類(lèi)型，分別為:25°以上坡耕地退耕還林;15°～25°坡耕地變?yōu)楣麍@;25°以上荒草坡還林。③岸邊植被過(guò)濾帶。選取不同植被種植在坡耕地的坡面形成植被過(guò)濾袋，使坡面漫流時(shí)的農(nóng)田非點(diǎn)源污染物質(zhì)得到截留和過(guò)濾。同時(shí)，植被過(guò)濾帶作為截留、緩沖、吸收坡面地表徑流的控制措施，可整體改善河流水質(zhì)。

5 結(jié)論

本文根據(jù)遼河流域非點(diǎn)源污染特征，結(jié)合流域氣候特征、土地利用結(jié)構(gòu)、人口分布、化肥農(nóng)藥施用量、水土流失狀況等，以GIS和地統(tǒng)計(jì)學(xué)技術(shù)為手段，采用多因子綜合分析法，得到遼河源頭區(qū)農(nóng)業(yè)非點(diǎn)源污染風(fēng)險(xiǎn)程度評(píng)價(jià)分值，根據(jù)不同風(fēng)險(xiǎn)程度劃分流域非點(diǎn)源污染控制區(qū)并得到了污染關(guān)鍵源區(qū)。結(jié)果表明，遼河源頭區(qū)非點(diǎn)源污染發(fā)生的關(guān)鍵源區(qū)占流域面積的21.77%，主要集中在梨樹(shù)南部、公主嶺東部和南部地區(qū)。關(guān)鍵源區(qū)主要分布在人口密度、畜禽養(yǎng)殖污染物以及施肥因子高風(fēng)險(xiǎn)等級(jí)區(qū)內(nèi)，且是坡耕地的主要分布區(qū)域。針對(duì)關(guān)鍵源區(qū)分析非點(diǎn)源風(fēng)險(xiǎn)高發(fā)原因提出一系列控制措施，根據(jù)研究區(qū)內(nèi)坡耕地水土流失情況嚴(yán)重，提出采取泥沙截留等方面的工程措施和適當(dāng)改變區(qū)域土地利用方式、調(diào)整坡耕地種植結(jié)構(gòu)和種植植被緩沖帶進(jìn)行污染物截留和緩沖;根據(jù)部分鄉(xiāng)鎮(zhèn)畜禽散養(yǎng)量大，對(duì)畜禽糞便處理和消納方式較為單一，提出應(yīng)合理布局畜禽養(yǎng)殖，推廣新型畜禽糞便處理技術(shù)。

本研究探討的農(nóng)業(yè)非點(diǎn)源風(fēng)險(xiǎn)多因子綜合評(píng)價(jià)和關(guān)鍵源區(qū)識(shí)別方法，有助于在流域尺度上直觀表達(dá)農(nóng)業(yè)非點(diǎn)源風(fēng)險(xiǎn)的空間差異，并可確定需要重點(diǎn)加以治理的關(guān)鍵區(qū)域，提出切實(shí)有效的管理控制措施，對(duì)遼河源頭區(qū)流域非點(diǎn)源污染的控制和管理起到積極作用，為源頭區(qū)河流生非點(diǎn)源污染態(tài)修復(fù)和遼河源頭區(qū)水源保護(hù)與水質(zhì)安全提供技術(shù)支撐。

［1］Potter K M，Cubbage F W.Awatershed－scale model for predicting non－point pollution risk in north carolina［J］.Environmental Management，2004，34(1):62－74.

［2］許其功，席北斗，沈珍瑤，等.三峽庫(kù)區(qū)典型小流域氮磷流失特征［J］.環(huán)境科學(xué)學(xué)報(bào)，2007，23(3):41－45.

［3］呂川，陳明輝，馬繼力，等.基于模糊聚類(lèi)分析的吉林省農(nóng)業(yè)非點(diǎn)源污染負(fù)荷現(xiàn)狀評(píng)價(jià)［J］.現(xiàn)代農(nóng)業(yè)科技，2011，8(5):263－268.

［4］Weng Q.Modeling urban growth effects on surface runoff with the integration of remote sensing and GIS［J］.Environmental Management，2001，28(6):737－748.

［5］李振煒，于興修，姚孝友，等.農(nóng)業(yè)非點(diǎn)源污染關(guān)鍵源區(qū)識(shí)別方法研究進(jìn)展［J］.生態(tài)學(xué)雜志，2011，30(12):2907－2914.

［6］李琪，陳利頂，齊鑫，等.流域尺度農(nóng)業(yè)磷流失危險(xiǎn)性評(píng)價(jià)與關(guān)鍵源區(qū)識(shí)別方法［J］.應(yīng)用生態(tài)學(xué)報(bào)，2007，18(9):1982－1986.

［7］呂家瓏.農(nóng)田土壤磷素淋溶及其預(yù)測(cè)［J］.生態(tài)學(xué)報(bào)，2003，23(12):2689－2701.

［8］Sivertun A，Prange L.Non－point source critical area analysis in the Gissel watershed using GIS［J］.Environmental Modelling Software，2003，18:887－898.

［9］劉曉，蘇維詞，王錚，等.基于RRM模型的三峽庫(kù)區(qū)重慶開(kāi)縣消落區(qū)土地利用生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)［J］.環(huán)境科學(xué)學(xué)報(bào)，2012，32(1):248－256.

［10］張維理，冀宏杰，Kolbe H，等.中國(guó)農(nóng)業(yè)面源污染形式估計(jì)及控制對(duì)策Ⅱ.歐美國(guó)家農(nóng)業(yè)面源污染狀況及控制［J］.中國(guó)農(nóng)業(yè)科學(xué)，2004，37(7):1018－1025.

［11］Lemunyon J L，Gilbert R G.Concept and need for a phosphorus assessment tool［J］.Journal of Production Agriculture，1993，6:484－485.

［12］余紅，沈珍瑤.大寧河流域非點(diǎn)源污染不確定性分析［J］.北京師范大學(xué)學(xué)報(bào):自然科學(xué)版，2008，44(1):86－91.

［13］李琪，陳利頂，齊鑫，等.妨水河流域農(nóng)耕區(qū)非點(diǎn)源磷污染危險(xiǎn)性評(píng)價(jià)與關(guān)鍵源區(qū)識(shí)別［J］.環(huán)境科學(xué)，2008，29(1):32－37.

［14］吉林省統(tǒng)計(jì)局.吉林省統(tǒng)計(jì)年鑒［Z］.北京:中國(guó)統(tǒng)計(jì)出版社，2008.

［15］麻素挺，湯潔，林年豐.基于GIS和RS多源空間信息的吉林省西部生態(tài)環(huán)境綜合評(píng)價(jià)［J］.資源科學(xué)，2004，26(4):140－145.

［16］楊育武，湯潔.脆弱生態(tài)環(huán)境指標(biāo)庫(kù)的建立及其定量評(píng)價(jià)［J］.環(huán)境科學(xué)研究，2002，15(4):46－49.

［15］Gburek W J，Sharpley A N，Heathwaite L，et al.Phosphorus management at the watershed scale:A modification of the phosphorus index［J］.Journal of Environmental Quality，2000，29:135－140.

［16］Heathwaite L，Sharpley A N，Gburek W J.A conceptual approach for integrating phosphorus and nitrogen management at watershed scales［J］.Journal of Environmental Quality，2000，29:158－166.

［17］Li Q，Cheng L，Qi X，et al.Assessing field vulnerability to phosphorus loss in Beijing agricultural area using revised field phosphorus ranking scheme［J］.Environment Science，2007，19:977－985.

［18］李瀟瀚，王瑞，卞建民，等.基于GIS的伊通河流域農(nóng)業(yè)非點(diǎn)源污染風(fēng)險(xiǎn)評(píng)估［J］.中國(guó)水土保持，2012，9(1):33－35.

［19］張淑榮，陳利頂，傅伯杰，等.農(nóng)業(yè)區(qū)非點(diǎn)源污染潛在危險(xiǎn)性評(píng)價(jià)——以于橋水庫(kù)流域磷流失為例［J］.第四紀(jì)研究，2003，23(3):262－269.

［20］Gburek W J，Sharpley A N.Hydrologic controls on phosphorus loss from upland agricultural watersheds［J］.Environmental Quality，1998，27:267－277.

［21］李琪，陳利頂，齊鑫，等.流域尺度農(nóng)業(yè)磷流失危險(xiǎn)性評(píng)價(jià)與關(guān)鍵源區(qū)識(shí)別方法［J］.應(yīng)用生態(tài)學(xué)報(bào)，2007，18(9):1982－1986.

［22］Hughes K J，Magette W L，Kurz I.Identifying critical source areas for phosphorus loss in Ireland using field and catchment scale ranking schemes［J］.Journal of Hydrology，2005，304:430－445.