沈 虹,張萬順,彭 虹

（1.北京大學(xué)城市與環(huán)境學(xué)院,地表過程分析與模擬教育部重點實驗室,北京100871;2.武漢大學(xué)資源與環(huán)境科學(xué)學(xué)院,武漢430079;3.武漢大學(xué)水利水電學(xué)院,武漢430072;4.北京大學(xué)深圳研究生院,城市人居環(huán)境科學(xué)與技術(shù)重點實驗室,廣東深圳 518055）

非點源污染研究最早始于20世紀(jì)60年代美國和日本所開展的土地利用與湖泊富營養(yǎng)化關(guān)系調(diào)查。而后,隨著計算機、3S技術(shù)的發(fā)展,以及野外實驗、室內(nèi)分析、水質(zhì)監(jiān)測等方法的改進(jìn)完善,描述地表水文過程、土壤侵蝕過程、污染物遷移轉(zhuǎn)化過程等的數(shù)學(xué)模型隨之出現(xiàn)[1],各種不同類型的非點源模型相繼問世并得到廣泛應(yīng)用[2-3]。非點源污染研究也從統(tǒng)計分析向半機理及機理過程轉(zhuǎn)變[4],研究尺度從小流域擴展到中大型流域甚至國家尺度。相對于國外,國內(nèi)研究起步雖晚但發(fā)展迅速,已有的研究案例覆蓋了不同尺度、不同類型的非點源,主要的研究關(guān)注點包括河網(wǎng)水系的分維影響[5]、負(fù)荷的模擬估算及空間分布特征[6]、污染物的遷移轉(zhuǎn)化[7]等。因此,建立定量化模型仍是當(dāng)前非點源污染研究的重點。

根據(jù)負(fù)荷產(chǎn)生和遷移過程,非點源負(fù)荷可以分為溶解態(tài)和顆粒態(tài)兩種類型[8]。溶解態(tài)負(fù)荷主要是指水溶性污染物,顆粒態(tài)負(fù)荷則是指在侵蝕性降雨作用下,被沖刷剝蝕吸附在土壤及泥沙中的污染物,它伴隨著泥沙遷移進(jìn)入河道,最后在河流中解析形成的水體污染[9]。水土流失、土地利用及水肥管理方式是顆粒態(tài)非點源負(fù)荷形成的重要原因,它不僅攜帶了大量地表肥沃的土壤,同時還帶走了大量累積在地表的磷污染物質(zhì)[10-11]。因此,有必要對土壤流失與顆粒態(tài)負(fù)荷的定量關(guān)系進(jìn)行深入研究。目前常用的研究方法是GIS與USLE結(jié)合,已廣泛應(yīng)用于黑河[12]、九龍江[13-14]、巢湖[15]等流域。

漢江作為長江的重要支流,其中下游水質(zhì)狀況對于保障沿岸飲水及工農(nóng)業(yè)生產(chǎn)安全具有重要的意義。然而,近年來隨著漢江中下游水質(zhì)的惡化,非點源污染已成為流域水體污染的主要影響因素之一。目前關(guān)于漢江中下游水土流失以及顆粒態(tài)非點源磷負(fù)荷的研究較少[16-17],研究范圍也主要集中在流域沿岸,多以行政意義界定流域邊界,忽視了污染負(fù)荷在遷移轉(zhuǎn)化過程中所依賴的自然地理要素,且方法以綜合層次判別分析和經(jīng)驗方法為主,缺少定量的模型化研究。因此,本研究基于GIS平臺,采用USLE和負(fù)荷輸出經(jīng)驗方程對漢江中下游土壤侵蝕和顆粒態(tài)非點源負(fù)荷進(jìn)行了定量估算。按照自然地理要素和時空過程一致性原則,以柵格作為基本計算單元,以子流域作為基本的匯總和模擬單元,分別研究了影響土壤侵蝕的五大因素,通過空間疊加方法得到了流域土壤侵蝕狀況,并在此基礎(chǔ)上估算了非點源負(fù)荷及其空間分布特征,以期為流域水土保持和顆粒態(tài)非點源防控提供理論依據(jù)。

1 流域概況及環(huán)境數(shù)據(jù)

1.1 流域概況

本研究采用的流域范圍以自然邊界為準(zhǔn),利用DEM基于集水區(qū)生成,地跨湖北、河南兩省(30°-34°N,100°-114°E),上起丹江口,下至漢江入長江河口,全長616 km,總流域面積6.078 4萬km2,涉及湖北省的26個縣、市、區(qū)和河南省的13個縣、市。流域內(nèi)擁有河南南陽盆地和湖北江漢平原兩大農(nóng)業(yè)區(qū),孕育了唐白河、南河、北河、小清河、蠻河、竹皮河、漢北河等河流水系。近年來漢江中下游先后5次發(fā)生富營養(yǎng)化現(xiàn)象,水質(zhì)狀況日益堪憂,非點源磷污染成為其重要原因之一。

1.2 基礎(chǔ)數(shù)據(jù)及預(yù)處理

本研究中所有的基礎(chǔ)數(shù)據(jù)(表1)均被納入到統(tǒng)一的坐標(biāo)系統(tǒng)下:Albers等面積圓錐投影,并轉(zhuǎn)化為100 m×100 mGrid格式。為了體現(xiàn)流域內(nèi)地理要素的空間異質(zhì)性,在ArcSWAT中對凹陷填充后的DEM進(jìn)行了再分割。經(jīng)過反復(fù)對比分析發(fā)現(xiàn),當(dāng)集水區(qū)面積閾值等于30 000 hm2時,生成的水系與實際情況具有最大相似性,此時共劃分31個子流域,得到流域總面積為6.0784萬km2(圖1);同時,在DEM 的基礎(chǔ)上進(jìn)行坡度、坡向分析,將坡度分為五組0～7.5°,7.5°～12.5°,12.5°～ 17.5°,17.5°～ 22.5°,22.5°～ 90°,合并破碎多邊形得到坡度圖。將坡向按東北、東、東南、南、西南、西、西北、北8個方向進(jìn)行分組得到坡向圖。解譯流域2007年Landsat TM影像,經(jīng)過幾何校正、影像拼接、直方圖匹配、顏色糾正等,得到疏林地、有林地、水田、旱地、建設(shè)用地、水域和未利用地7類土地利用類型;土壤數(shù)據(jù)來源于1∶400萬中國土壤分布圖、中國土壤數(shù)據(jù)庫、湖北省和河南省土種志,包括土壤的空間分布、質(zhì)地組成、有機質(zhì)含量、TP背景值等,制作土壤磷元素含量和其他氣象要素等相關(guān)專題圖。

表1 基礎(chǔ)數(shù)據(jù)及說明

圖1 漢江中下游子流域分布圖

2 研究方法

以柵格為非點源的基本計算單元,以子流域為非點源負(fù)荷匯總的基本單元,基于GIS、RS技術(shù),采用ULSE方程和負(fù)荷輸出經(jīng)驗方程對顆粒態(tài)非點源負(fù)荷量進(jìn)行估算。

2.1 土壤侵蝕模數(shù)估算

USLE(universe soil loss equation)是一種經(jīng)驗統(tǒng)計模型,它由Wischmeier于20世紀(jì)60年代提出。模型的基本思想是侵蝕模數(shù)大小取決于引起侵蝕的力與土壤抵抗降雨沖刷的能力的組合,計算結(jié)果表述為降雨引起的單位面積土壤流失量[18]。目前,將USLE與GIS/RS結(jié)合已經(jīng)成為土壤侵蝕研究的重要方法[19-20]。

式中:A——年均土壤侵蝕量[t/(hm2?a)];α——單位轉(zhuǎn)換系數(shù),為224.2;R——降雨侵蝕力因子[(MJ?mm)/(hm2?h?a)];K代表土壤可蝕性因子[(t?hm2?h)/(hm2?MJ?mm)];LS——坡度坡長因子,無量綱;C——經(jīng)營管理因子,無量綱;P——水土保持措施因子,無量綱。

2.1.1 降雨侵蝕力因子R 降雨侵蝕力旨在評價降雨引起土壤侵蝕的潛在能力,它與降雨量、降雨歷時、降雨強度、雨滴大小等因素有關(guān)[21]。利用研究區(qū)1961-2006年月平均降雨量和年平均降雨量資料,根據(jù)Wischmeier經(jīng)驗公式估算R[22],表達(dá)式為

式中:Pi——月均降雨量(mm);P——年均降雨量(mm)。

2.1.2 土壤可蝕性因子K 土壤可蝕性因子K是土壤侵蝕性的量度,它反映了不同類型土壤抵抗侵蝕能力的高低,與土壤質(zhì)地和有機質(zhì)含量密切相關(guān)。目前K估算主要基于小區(qū)實驗展開,常用的經(jīng)驗估算辦法是Williams等人提出的土壤質(zhì)地、有機質(zhì)含量測定方法。本研究在采集漢江中下游土壤資料、圖件、土壤類型空間分布等資料的基礎(chǔ)上,通過查表和參考眾多文獻(xiàn)研究[23-25],確定了研究區(qū)土壤可蝕性因子K值。

表2 漢江中下游土壤可蝕性因子K值

2.1.3 坡長坡度因子LS 坡度坡長因子代表地形因素變化對土壤侵蝕的影響。坡度、坡長與侵蝕作用力呈現(xiàn)正向關(guān)系,當(dāng)坡度和坡長因子越陡或越長時,徑流速度將越大,相應(yīng)的,在坡面匯集的徑流所引起的侵蝕作用力也將增強。本研究引用王寧等[26]關(guān)于地形因子LS的估算方法,得到坡長坡度因子LS的計算表達(dá)式如下:

式中:LS——坡度坡長因子,無量綱;l——坡長(m);θ——坡度(°);m——坡長指數(shù),其值與坡度相關(guān)。

2.1.4 經(jīng)營管理因子C 經(jīng)營管理因子C是指在相同土壤、坡度和降雨條件下,某種特定土地覆被的土壤流失量與耕作后連續(xù)休閑地的土壤流失量比值,可以反映地表覆蓋、土地利用以及一般管理措施對土壤侵蝕速度的影響。由于經(jīng)營管理因子的影響因素非常多,通常采用簡化處理,僅考慮最主要影響因素植被覆蓋度與經(jīng)營管理因子C的作用關(guān)系[27-28]。資料研究證實,NDVI與植被覆蓋具有良好的線性關(guān)系[29]。因此本研究將遙感解譯的NDVI作為橋梁,利用NDVI推導(dǎo)出植被蓋度(公式4),進(jìn)而建立植被蓋度與經(jīng)營管理因子的經(jīng)驗公式(公式5)。其值介于0～1,當(dāng)植被蓋度為0時,土地裸露,最易造成土壤侵蝕,因此經(jīng)營管理因子 C最大;當(dāng)C≥78.3%時,植被蓋度較高,可以阻擋降雨動能造成的侵蝕,近似認(rèn)為此時經(jīng)營管理最完善,經(jīng)營管理因子C賦值為0。

式中:NDVI(Normalized difference vegetation index)——植被歸一化指數(shù),無量綱;lc——植被蓋度,無量綱;C——經(jīng)營管理因子,無量綱。

2.1.5 水土保持因子P 水土保持因子P指采取水土保持措施后的土壤流失量與順坡種植時的土地流失量比值,反映了水土保持措施對于土壤侵蝕的控制作用,數(shù)值介于0～1。目前,水土保持因子P的求取方法有基于坡度分級賦值和基于不同土地利用賦值等方式。

考慮到研究區(qū)的特點,本研究采用土地利用賦值方法。在查閱流域水保措施資料及參考大量文獻(xiàn)的基礎(chǔ)上[30-31],確定研究區(qū)水土保持因子P值。水土保持因子P值越大表明水保措施越差,土壤潛在侵蝕的可能性越大,例如未利用地通常為裸地,可近似認(rèn)為未采取任何水保措施,因此水土保持因子P值為1;而對于水域,由于不會產(chǎn)生水土流失,因此可視為擁有完善的水保措施,賦值為0。

表3 漢江中下游水土保持因子P值

2.2 顆粒態(tài)氮磷負(fù)荷估算

參照相關(guān)研究確定顆粒態(tài)非點源負(fù)荷的經(jīng)驗輸出模型為[7,13-15,17,21,32]:

式中:Wxp——顆粒態(tài)非點源磷負(fù)荷[kg/(hm2?a)];α——單位換算常數(shù),無量綱;Cp——顆粒態(tài)磷負(fù)荷濃度,來源于中國科學(xué)院南京土壤研究所提供的中國土壤數(shù)據(jù)庫、湖北省和河南省土種志;A——土壤流失量[t/(km2?a)];η——富集系數(shù),無量綱 ;Sd——流域土壤侵蝕遷移到河道的運移比,無量綱。

(1)土壤磷濃度計算。土壤中磷含量是土壤流失的直接補給,它反映了土壤流失的潛在能力。本研究中采用的土壤磷濃度來源于中國土壤數(shù)據(jù)庫、湖北省和河南省土種志(表4)。

表4 土壤中顆粒態(tài)非點源磷負(fù)荷輸出濃度 %

(2)富集比。侵蝕形成的泥沙具有富集養(yǎng)分元素的特性,同雨前表土養(yǎng)分含量相比,侵蝕形成的泥沙養(yǎng)分含量較高,這種現(xiàn)象被稱為泥沙富集作用[33],一般用某養(yǎng)分在侵蝕泥沙中的含量與其在被侵蝕土壤中的含量之比表示。本研究在參閱同類研究的基礎(chǔ)上[12,17,33],確定研究區(qū)富集比為2。

(3)遷移比Sd。顆粒態(tài)非點源磷負(fù)荷伴隨土壤侵蝕、遷移的過程中會出現(xiàn)損耗,只有真正進(jìn)入河道中的污染負(fù)荷才會對流域水質(zhì)造成影響,并且距離河道的遠(yuǎn)近直接影響著運移損耗。根據(jù)長江水利委員會對長江流域長期定位的研究結(jié)果:長江流域的Sd值大約為0.1～0.4[17],據(jù)此確定研究區(qū)內(nèi)平均泥沙輸移比為0.28,最大值0.4,最小值0.1。利用Arc-GIS空間分析工具計算每個柵格單元到最近河道的距離,按照距離倒數(shù)賦值輸移比,得到 Sd空間分布圖。

3 結(jié)果與分析

3.1 土壤侵蝕估算及驗證

2007年研究區(qū)土壤的平均侵蝕模數(shù)為867.55 t/(km2?a),最高侵蝕模數(shù)為子流域 21,高達(dá)5 158.75 t/(km2?a),分布在神農(nóng)架林區(qū)。最低侵蝕模數(shù)為子流域28,低至111.42 t/(km2?a),分布在竹皮河流域附近。按照1997年水利部頒布的《土壤侵蝕分類分級標(biāo)準(zhǔn)》將研究區(qū)土壤侵蝕劃分為八級,得到侵蝕程度分類圖(附圖1)。總體上,漢江中下游流域?qū)儆谳p度侵蝕區(qū),但依然有5.69%的區(qū)域侵蝕模數(shù)大于2 500 t/(km2?a),其中,中度侵蝕1.8%,強度侵蝕1.21%,極強度侵蝕1.32%,劇烈侵蝕1.37%,土壤侵蝕的主要分布區(qū)以中低山和丘陵地帶為主,而土壤侵蝕的最嚴(yán)重區(qū)域則以中高山農(nóng)耕地為主,如神農(nóng)架林區(qū)、房縣、?？怠⒐瘸遣糠謪^(qū)域及嵩縣、內(nèi)鄉(xiāng)等地帶,這與文獻(xiàn)[16]的研究結(jié)論較為一致。

3.2 非點源負(fù)荷空間分布

研究區(qū)顆粒態(tài)非點源磷負(fù)荷為8 496.79 t/a,對比2002年史志華等人[17]在該流域的計算結(jié)果8 893.4 t/a,顆粒態(tài)非點源磷負(fù)荷的總量有所減少,但空間分布具有一致性,支流唐白河、南河、蠻河等非點源磷負(fù)荷的入?yún)R量較高。從子流域單位負(fù)荷的角度分析,子流域21,16,15,12,22,2,24顆粒態(tài)磷單位負(fù)荷較高(附圖2),其中最高為子流域21,單位負(fù)荷高達(dá)8.42 kg/(hm2?a)。具體來說,這些子流域主要分布在三片區(qū)域,一是神農(nóng)架林區(qū)、房縣、?？?、谷城、南漳西北部一帶,平均顆粒態(tài)磷單位負(fù)荷為4.1 kg/(hm2?a);二是嵩縣、內(nèi)鄉(xiāng)一帶,平均顆粒態(tài)磷單位負(fù)荷量為3.0 kg/(hm2?a);三是鐘祥東北部,平均顆粒態(tài)磷單位負(fù)荷為1.45 kg/(hm2?a)。在神農(nóng)架、房縣、?？狄粠щm然有林地面積大,但海拔高,坡度大,平均坡度在22.5°以上,且在經(jīng)濟發(fā)展的帶動下,將坡地改成耕地種植高山玉米、蔬菜等現(xiàn)象較普遍,因此土壤流失最為嚴(yán)重,隨水土流失攜帶的顆粒態(tài)非點源負(fù)荷亦最高。在嵩縣、內(nèi)鄉(xiāng)一帶,其北部平均坡度達(dá)22.5°以上,坡耕地亦是造成土壤流失和非點源負(fù)荷的最重要原因;而在南部,作為南陽盆地的一部分,磷肥的大量施用是顆粒態(tài)負(fù)荷增高的重要原因之一。在鐘祥東北部,山地丘陵是最主要的地貌類型,土地利用方式多為旱地、疏林地、未利用地等,植被覆蓋度較低,導(dǎo)致坡長坡度因子LS、經(jīng)營管理因子C、水土保持措施因子P較高,因此土壤流失產(chǎn)生的非點源負(fù)荷量高,此外,分布在該地區(qū)的石灰土,土壤T P含量較高也是部分原因。自鐘祥以下的漢江下游河段,顆粒態(tài)非點源磷負(fù)荷均較低,平均單位磷負(fù)荷為0.3 kg/(hm2?a),見圖2。

表5 漢江中下游土壤侵蝕模數(shù)分級表

圖2 子流域顆粒態(tài)非點源單位負(fù)荷量和總量關(guān)系

表6 各支流非點源磷負(fù)荷的匯入情況

在實際負(fù)荷遷移轉(zhuǎn)化過程中,顆粒態(tài)負(fù)荷不是直接通過子流域匯入干流,而是由支流攜帶進(jìn)入主干河道,因此將子流域按照主要支流進(jìn)行匯總分析有利于指導(dǎo)支流的非點源防控工作。在本研究中,將31個子流域按構(gòu)成分組為9個主要支流和直接入江10種情形,分別為北河、南河、小清河、唐白河、蠻河、雙河、利河、竹皮河、漢北河和直接入江。從支流入?yún)R的非點源負(fù)荷分析(表6),研究區(qū)內(nèi)顆粒態(tài)非點源磷負(fù)荷主要來源于支流入?yún)R,直接入江的顆粒態(tài)非點源磷負(fù)荷僅占16.87%。唐白河由于流域面積大,同時又流經(jīng)“中州糧倉”南陽盆地,農(nóng)業(yè)磷肥廣泛施用使得流失的土壤中含有大量的顆粒態(tài)非點源磷負(fù)荷,因而成為最大的顆粒態(tài)非點源磷負(fù)荷貢獻(xiàn)者,總量達(dá)到3 079.26 t,貢獻(xiàn)率達(dá)36.24%;南河因流經(jīng)神農(nóng)架、?？怠⒎靠h等土壤流失最嚴(yán)重區(qū),顆粒態(tài)非點源磷負(fù)荷較高,達(dá)到2 454.45 t,貢獻(xiàn)率為28.89%,居第二位。這里山地綿延、地勢陡峻,且坡地耕作,植被漸被破環(huán),侵蝕流失的負(fù)荷量大;蠻河流經(jīng)南漳區(qū)域,磷負(fù)荷量也較高,位居第三。北河、漢北河、利河的顆粒態(tài)非點源磷負(fù)荷貢獻(xiàn)率居中,為2.75%～3.72%。小清河、雙河、竹皮河因面積小,單位負(fù)荷量低,因而對顆粒態(tài)氮磷負(fù)荷總量的貢獻(xiàn)率低,磷的貢獻(xiàn)率均小于等于1.23%。

4 結(jié)論

(1)漢江中下游2007年土壤的平均侵蝕模數(shù)為867.55 t/(km2?a),總體上屬于輕度侵蝕,但依然有5.69%的區(qū)域侵蝕模數(shù)大于2 500 t/(km2?a)。流域內(nèi)由于土壤流失引起的顆粒態(tài)非點源磷負(fù)荷流失量達(dá)8 496.79 t/a,其中南河、唐白河、蠻河等支流的非點源單位負(fù)荷量和總負(fù)荷量均較高,最高單位負(fù)荷量可達(dá)8.42 kg/(hm2?a),與水土流失區(qū)嚴(yán)重區(qū)基本一致。該區(qū)域主要分布在神農(nóng)架林區(qū)、房縣、?？档壬降厍鹆暌约澳详柵璧氐绒r(nóng)業(yè)發(fā)達(dá)地區(qū),與坡耕地耕作、農(nóng)田化肥施用緊密相關(guān),是今后水土流失和非點源防控的關(guān)鍵區(qū),因此在防控非點源工作中,要從顆粒態(tài)負(fù)荷產(chǎn)生的源頭做起,一是要加大森林覆蓋,二是要嚴(yán)禁坡地耕作,防控和減少水土流失。

(2)由于部分河口泥沙等相關(guān)數(shù)據(jù)較難獲取,因此本研究中無法利用實測數(shù)據(jù)對估算的顆粒態(tài)負(fù)荷量進(jìn)行驗證。同時,在本研究中以負(fù)荷的經(jīng)驗輸出公式代替了負(fù)荷復(fù)雜的遷移轉(zhuǎn)化機理,在一定程度上會影響估算的精度,因此有待在今后的工作中進(jìn)一步改進(jìn)和完善。

[1]程紅光,岳勇,楊勝天,等.黃河流域非點源污染負(fù)荷估算與分析[J].環(huán)境科學(xué)學(xué)報,2006,26(3):384-391.

[2]Michele Munafòa,Giuliano Cecchib,Fabio Baioccoa,et al.River pollution from non-point sources:a new simplified method of assessment[J].Journal of Environmental Management,2005,77(2):93-98.

[3]Polyakova V,Faresa A,Kubob D,et al.Evaluation of a non-point source pollution model,AnnAGNPS,in a tropical watershed[J].Environment Modelling and Software,2007,22(11):1517-1627.

[4]張秋玲,陳英旭,俞巧鋼,等.非點源污染模型研究進(jìn)展[J].應(yīng)用生態(tài)學(xué)報,2007,18(8):1886-1890.

[5]榮琨,陳興偉,王林,等.流域非點源污染與水系分維關(guān)系研究[J].水土保持研究,2009,16(2):128-131.

[6]王曉燕,王曉峰,汪清平,等.北京密云水庫小流域非點源污染負(fù)荷估算[J].地理科學(xué),2004,24(2):227-231.

[7]李懷恩,沈晉.流域非點源污染模型的建立與應(yīng)用實例[J].環(huán)境科學(xué)學(xué)報,1997,17(2):141-147.

[8]李繼承.嘉陵江流域非點源污染負(fù)荷模擬研究[D].重慶:重慶大學(xué),2007.

[9]Bardgett R D,Anderson J M,Behan-Pelletier V,et al.The influence of soil biodiversity on hydrological pathways and the transfer of materials between terrestrial and aquatic ecosystem[J].Ecosystems,2001,4:421-429.

[10]程文娟,史靜,夏運生,等.滇池流域農(nóng)田土壤氮磷流失分析研究[J].水土保持學(xué)報,2008,22(5):52-55.

[11]張玉斌,鄭粉莉,武敏.土壤侵蝕引起的農(nóng)業(yè)非點源污染研究進(jìn)展[J].水科學(xué)進(jìn)展,2007,18(1):123-132.

[12]胥彥玲,李懷恩,倪永明,等.基于 USLE的黑河流域非點源污染定量研究[J].西北農(nóng)林科技大學(xué)學(xué)報:自然科學(xué)版,2006,34(3):138-142.

[13]黃金良,洪華生,張珞平,等.基于GIS和 USLE的九龍江流域土壤侵蝕量預(yù)測研究[J].水土保持學(xué)報,2004,18(5):75-79.

[14]黃金良,洪華生,張珞平,等.基于GIS的九龍江流域農(nóng)業(yè)非點源氮磷負(fù)荷估算研究[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2004,23(5):866-871.

[15]汪曉輝,巢湖流域非點源N、P污染排放負(fù)荷估算及控制研究[D].合肥:合肥工業(yè)大學(xué),2006.

[16]史志華,蔡崇法,丁樹文,等.基于GIS的漢江中下游土壤侵蝕強度初步分析[J].水土保持學(xué)報,2001,15(6):14-17.

[17]史志華,蔡崇法,丁樹文,等.基于GIS的漢江中下游農(nóng)業(yè)面源氮磷負(fù)荷研究[J].環(huán)境科學(xué)學(xué)報,2002,22(4):473-477.

[18]曾凌云.基于RUSLE模型的喀斯特地區(qū)土壤侵蝕研究:以貴州紅楓湖流域為例[D].北京:北京大學(xué),2008.

[19]徐清艷,周躍.基于GIS和 USLE的突然侵蝕預(yù)測研究:以云南新平大紅山鐵礦為例[J].水土保持研究,2009,16(4):72-75.

[20]張有全,宮輝力,趙文吉,等.基于GIS和 USLE的密云縣土壤侵蝕評價及空間特征研究[J].水土保持研究,2007,14(3):358-364.

[21]楊勝天,程紅光,步青松,等.全國土壤侵蝕量估算及其在吸附態(tài)氮磷流失匡算中的應(yīng)用[J].環(huán)境科學(xué)學(xué)報,2006,26(3):366-374

[22]Wischmeier W H,Smith D D.Predicting Rainfall Erosion losses:a Guide to Conservation[M].Washington DC:US Department of Agriculture,Science and Education Administration,Agricultural Handbook,1978.

[23]汪東川.基于ArcGIS的坡面產(chǎn)流產(chǎn)沙過程模擬:以重慶市璧山縣為例[D].重慶:西南大學(xué),2005.

[24]鄧良基,侯大斌,王昌全,等.四川自然土壤和旱耕地土壤可蝕性特征研究[J].中國水土保持,2003(7):23-25.

[25]張科利,彭文英,楊紅麗.中國土壤可蝕性值及其估算[J].土壤學(xué)報,2007,44(1):7-13.

[26]王寧,朱顏明,徐崇剛.GIS用于流域徑流污染物的量化研究[J].東北師范大學(xué)學(xué)報:自然科學(xué)版,2002,34(2):92-98.

[27]趙琰鑫,張萬順,王艷,等.基于3S技術(shù)和USLE的深圳市茜坑水庫流域土壤侵蝕強度預(yù)測研究[J].亞熱帶資源與環(huán)境學(xué)報,2007,2(3):23-28.

[28]蔡崇法,丁樹文,史志華,等.應(yīng)用 USLE模型與地理信息系統(tǒng)IDRISI預(yù)測小流域土壤侵蝕量的研究[J].水土保持學(xué)報,2000,14(2):19-24.

[29]高龍華.遙感和GIS支持下流域非點源污染模型研究[D].南京:河海大學(xué),2006.

[30]卜兆宏,孫金莊.水土流失定量遙感方法及應(yīng)用研究[J].土壤學(xué)報,1997,34(3):235-245.

[31]許月卿,邵曉梅.基于GIS和RUSLE的土壤侵蝕量計算:以貴州省貓?zhí)恿饔驗槔齕J].北京林業(yè)大學(xué)學(xué)報,2006,28(4):67-71.

[32]施為光.四川省清平水庫流域非點源污染負(fù)荷計算[J].重慶環(huán)境科學(xué),2002,22(2):33-36.

[33]薛素玲.基于GIS的黑河流域非點源氮磷模擬[D].西安:西安理工大學(xué),2006.