杜書立， 付 微， 趙 軍， 張興義

(中國科學(xué)院 東北地理與農(nóng)業(yè)生態(tài)研究所院 黑土重點實驗室， 黑龍江 哈爾濱 150081)

東北黑土區(qū)是世界上4大片黑土區(qū)之一，是中國最主要的糧食生產(chǎn)基地之一。東北黑土區(qū)開發(fā)雖只有百余年[1]，但由于受到降雨、徑流、地形、地表植被、人類活動、土地利用類型以及耕作管理等眾多因素的影響，土壤侵蝕嚴(yán)重，主要表現(xiàn)在坡耕地黑土層被剝蝕而變薄，侵蝕溝縱橫造成土地支離破碎[2-5]。據(jù)調(diào)查統(tǒng)計，黑土區(qū)每年流失掉的黑土總量達(dá)1×108～2×108m3，其中流失的氮磷鉀養(yǎng)分就相當(dāng)于數(shù)百萬噸化肥，土壤中有機質(zhì)含量比開墾前下降了近1/2[6-7]。黑土區(qū)侵蝕溝已經(jīng)達(dá)到2.5×105條余，侵蝕耕地5.9×105hm2以上[1]。東北黑土區(qū)的土壤侵蝕問題已引起了廣泛的關(guān)注，近年來，國家和地方逐步加大防治力度[8-10]。GIS與RS已成為當(dāng)前國家和區(qū)域水土流失普查的主要手段，剛完成的“第一次全國水利普查”中的全國水土保持普查就是主要借助GIS與RS技術(shù)[11-12]。利用空間數(shù)據(jù)，對侵蝕區(qū)進(jìn)行整體水土保持規(guī)劃、水土流失治理和效益評價，為直觀認(rèn)識水土流失危害和科學(xué)防治提供了重要手段。GIS具有強大的設(shè)計、規(guī)劃和管理決策功能，通過預(yù)測水土流失，宏觀范圍上直觀地體現(xiàn)出侵蝕現(xiàn)狀，并且可以分析不同類型的數(shù)據(jù)，連接不同的模型形式，從而建立合理的GIS評價模型[13]。同時GIS具有非常人性化的可視界面，非常有利于相關(guān)人員進(jìn)行深入分析以及正確決策。

選取東北黑土區(qū)松嫩平原中心地帶海倫市光榮村作為研究區(qū)，光榮村土壤侵蝕區(qū)域多分布在耕地區(qū)域，在東北黑土區(qū)土壤侵蝕和水土流失表現(xiàn)上具有代表性。結(jié)合ArcGIS，SuperMap和SuperObject等軟件，建立土壤侵蝕預(yù)警和決策支持系統(tǒng)，旨在為黑土區(qū)水土流失綜合治理提供數(shù)字化平臺。

1 研究區(qū)概況

松嫩平原黑土區(qū)光榮村位于黑龍江省海倫市，地理坐標(biāo)為47°20′N，127°48′E。耕地面積450 hm2，占總土地面積的85.1%，農(nóng)田地塊723塊，面積0.3～7 hm2不等，主要種植作物為大豆和玉米[14]。光榮村示范區(qū)地貌類型為平川漫崗地，土壤類型為黑土，土層較厚，平均厚度46.6 cm左右；土壤有機質(zhì)含量平均約為48.4 g/kg；當(dāng)前土壤肥力總體比較肥沃，但侵蝕強度大，面蝕和溝蝕均嚴(yán)重，且仍處于發(fā)展趨勢。

2 土壤侵蝕方程及其各參數(shù)的計算

2.1 土壤侵蝕方程

土壤侵蝕方程是本系統(tǒng)的核心部分，根據(jù)前人的研究以及參考眾多參考文獻(xiàn)，選擇張憲奎等[15-18]在RUSLE模型基礎(chǔ)上建立的黑龍江省土壤流失方程作為土壤侵蝕方程，具體形式如下：

A=R×K×L×S×C×P

(1)

式中：A——單位面積多年平均土壤流失量〔t/(hm2·a)〕；R——降雨侵蝕力因子(J/m2)；K——土壤可蝕性因子〔(t·hm2·h)/(hm2·MJ·mm)〕；L——坡長(m)；S——坡度(%)；C——植被與經(jīng)營管理因子；P——水土保持措施因子。

2.2 參數(shù)計算

根據(jù)張憲奎等[18]的研究，降雨侵蝕力因子R為：

R=EI

(2)

式中：E——一次降雨中某時段降雨量產(chǎn)生的動能(J/m2)；I——某時段降雨強度(cm/h)。利用1980年以來海倫地區(qū)的降水?dāng)?shù)據(jù)，運用張憲奎等的方法，計算得出R=1 225.5(J/m2)。

試驗區(qū)光榮村地處黑土帶的中心區(qū)域，根據(jù)張憲奎等[18]的計算，黑土區(qū)土壤可蝕性因子為K=0.26 〔(t·hm2·h)/(hm2·MJ·mm)〕。

對于地形因子LS，國內(nèi)外學(xué)者提出了很多種不同的計算方法[19-23]。本文使用ArcGIS軟件，根據(jù)光榮實驗區(qū)的DEM數(shù)據(jù)計算坡長因子L和坡度因子S。

植被與經(jīng)營管理因子C是RUSLE中最重要的因子，RUSLE用次因子法計算各種土地利用方式的特定C值是建立在植物冠層、雨滴降落高度、地表覆蓋度和糙度、根量及前茬土地利用方式、土壤前期含水量等基礎(chǔ)之上，從而提供更靈活的能描述種植制度和保土耕作措施或牧草輪作等的C值估計方法[24]。根據(jù)光榮地區(qū)的實際種植作物以及其他植被情況，結(jié)合張憲奎等[18]的研究結(jié)果，我們選取植被與經(jīng)營管理因子C分別為0.262 6(大豆)、0.257 8(玉米)、0.06(林地)和0.15(草地)。

水土保持措施因子P，是指不同的土壤侵蝕治理措施下，農(nóng)地土壤流失量與順坡耕作農(nóng)地土壤流失量之比值。根據(jù)張憲奎等[18]的研究，我們選取P因子的數(shù)值為1(順坡壟)和0.029(橫坡壟)，而在林地與草地地類中，并不存壟向，各種地表植被的生長方式以及位置分布是不符合壟向分布的，因此在分布上對水土保持并沒有很好的效果，所以林地和草地的水土保持措施因子P值均為1。

2.3 建立土壤侵蝕管理地塊圖

以地塊為基本單元逐個調(diào)查近3 a來的降雨強度、頻率以及種植作物，將調(diào)查結(jié)果整理并編碼，建立空間檢索，關(guān)聯(lián)到專家知識庫。同樣以地塊為基本單元，對地塊的相應(yīng)信息進(jìn)行采集，包括土壤類型、土地利用類型、壟向、坡向以及土壤中氮磷鉀等重要養(yǎng)分的含量等。

首先，利用ArcGIS軟件對空間分辨率為2.5 m的 ALOS遙感數(shù)據(jù)進(jìn)行解譯得到土地利用圖，同時進(jìn)行矢量化，并與實地采集土壤養(yǎng)分?jǐn)?shù)據(jù)和其他調(diào)查數(shù)據(jù)進(jìn)行關(guān)聯(lián)[14]。然后，根據(jù)光榮村的DEM，利用ArcGIS軟件計算光榮試驗區(qū)各地塊的實際坡長。根據(jù)張憲奎等[18]建立的黑龍江省土壤流失方程作為土壤侵蝕模型，將地類、坡度、坡長等進(jìn)行分類整理。地類分為耕地、林地、草地以及非農(nóng)業(yè)用地，由于侵蝕溝侵蝕嚴(yán)重，在修復(fù)前無法進(jìn)行利用，因此，將侵蝕溝作為獨立的一個類型，不歸于其他地類之中；坡度根據(jù)>0.5°，0.5°～2°，2°～3°，3°～4°，4°～5°，5°～6°，6°～7°以及>7°分級；坡長根據(jù)<100，100～200，200～500及>500 m分級；對于耕地，作物壟向也是土壤侵蝕的重要因素，壟向分為橫坡壟和順坡壟，作為農(nóng)田土壤侵蝕的重要因子。根據(jù)以上分級，去掉不存在的組合部分，共將光榮村試驗區(qū)所有地塊分為34類，作為侵蝕管理單元模塊，從而制定不同的治理措施。表1中列出了這34類不同組合的具體組合形式，即34類地塊各自的地類、坡度范圍、坡長范圍以及壟向。附圖5根據(jù)表1中的不同組合，將光榮村的34類地塊賦以不同的顏色，從而進(jìn)行分類管理和治理。

表1 光榮村土壤侵蝕管理地塊組合及各因子范圍

3 系統(tǒng)開發(fā)及數(shù)據(jù)庫的建立

3.1 開發(fā)方式

可視化二次開發(fā)是現(xiàn)今GIS應(yīng)用開發(fā)的主要發(fā)展方向。GIS的軟件平臺與可視化二次開發(fā)軟件相結(jié)合，具有很強的應(yīng)用性，不僅快速方便，而且易于各個平臺之間的移植。本文是以超圖公司的SuperMap Object為GIS基礎(chǔ)軟件，利用Visual Basic語言作為開發(fā)程序，二者相集成的開發(fā)方式實現(xiàn)。

3.2 空間數(shù)據(jù)庫建立

空間數(shù)據(jù)和屬性數(shù)據(jù)進(jìn)行分離式儲存，通過關(guān)鍵字段進(jìn)行連接調(diào)用。為了方便屬性數(shù)據(jù)的編輯和更新，將屬性數(shù)據(jù)置于關(guān)系數(shù)據(jù)庫系統(tǒng)之中?？臻g數(shù)據(jù)中包括光榮實驗區(qū)的地形數(shù)據(jù)、土地利用數(shù)據(jù)、土壤類型數(shù)據(jù)等，并且利用ALOS影像數(shù)據(jù)解譯得到光榮村的地塊數(shù)據(jù)；屬性數(shù)據(jù)主要包括一些實地調(diào)查數(shù)據(jù)等。分別對空間數(shù)據(jù)和屬性數(shù)據(jù)進(jìn)行整理：空間數(shù)據(jù)配準(zhǔn)、矢量化，并且分類；屬性數(shù)據(jù)按照相同格式進(jìn)行規(guī)范化、標(biāo)準(zhǔn)化。屬性數(shù)據(jù)可以通過SQL數(shù)據(jù)庫進(jìn)行數(shù)據(jù)訪問進(jìn)而進(jìn)行屬性數(shù)據(jù)查詢等，并且可以通過二次開發(fā)語言與空間圖形數(shù)據(jù)進(jìn)行交互性操作。屬性數(shù)據(jù)庫的調(diào)用完全在后臺執(zhí)行，圖形和屬性數(shù)據(jù)的查詢與維護(hù)在同界面下實現(xiàn)。

4 系統(tǒng)設(shè)計與功能

4.1 系統(tǒng)設(shè)計

4.1.1 系統(tǒng)設(shè)計目標(biāo) 光榮示范區(qū)土壤侵蝕預(yù)警決策支持系統(tǒng)是用于為不同土壤條件和不同侵蝕強度下的土壤侵蝕綜合治理的GIS管理平臺。具體目標(biāo)是：實現(xiàn)直觀的可視化的地塊信息查看，主要是包括侵蝕影響因子和土壤養(yǎng)分信息并且提供相應(yīng)地塊的侵蝕強度；為管理者提供宏觀決策分析平臺，使管理者能夠?qū)ο到y(tǒng)進(jìn)行維護(hù)，數(shù)據(jù)更新、查詢、分析；為農(nóng)民用戶提供根據(jù)不同侵蝕情況而制定的侵蝕綜合治理措施，如生物措施，工程措施，合理施肥等。

4.1.2 系統(tǒng)框架 由于使用者層次和專業(yè)不一，因此系統(tǒng)框架設(shè)計時注重以用戶為主，界面整潔、易懂，操作簡單，人性化強。

對于空間數(shù)據(jù)和屬性數(shù)據(jù)，采用分離的方式進(jìn)行儲存，之后在VB中編寫接口，調(diào)用ODBC與數(shù)據(jù)庫源進(jìn)行連接，之后再利用RDO和SQL進(jìn)行數(shù)據(jù)的訪問查看。主要包括數(shù)據(jù)輸入和輸出，數(shù)據(jù)顯示，空間數(shù)據(jù)的分析，專題圖制作，系統(tǒng)配置與維護(hù)，數(shù)據(jù)庫，侵蝕治理措施等模塊(見圖1)。

圖1 土壤侵蝕預(yù)警和決策支持系統(tǒng)集成過程

4.2 系統(tǒng)功能

在進(jìn)行系統(tǒng)功能和界面設(shè)計時，以方便實用為原則，系統(tǒng)界面簡潔，具有較好的實用性。

(1)地圖基本操作。具有基本的GIS功能，包括放大、縮小、自由縮放、移動和全圖顯示。這些功能直觀，且便于操作，適合不同人群使用。

(2)地形查詢。地塊的坡度范圍是侵蝕強度的重要指標(biāo)，通過地形狀況功能，可以清晰簡易地顯示出各個不同坡度范圍的地塊，同時對于已經(jīng)形成的侵蝕溝也會重點標(biāo)注，這樣非常有助于使用者快速定位不同侵蝕強度的地塊和侵蝕溝位置。

(3)侵蝕預(yù)警與治理功能。侵蝕預(yù)警和治理功能是本系統(tǒng)的重點，主要是利用張憲奎等[18]在RUSLE模型基礎(chǔ)上建立的黑龍江省土壤流失方程作為土壤侵蝕方程為驅(qū)動，進(jìn)行侵蝕模數(shù)計算，之后根據(jù)《黑土區(qū)水土流失綜合防治技術(shù)標(biāo)準(zhǔn)》[25]進(jìn)行侵蝕強度分級，并且根據(jù)不同的侵蝕強度制定綜合治理方案。

4.3 系統(tǒng)應(yīng)用

土壤侵蝕預(yù)警功能模塊中主要包括地塊狀況、地塊信息、侵蝕預(yù)警、綜合治理方案、效益評估等功能。本系統(tǒng)針對2個層次的用戶。第一是管理者，他們有宏觀分析和決策的能力，可以利用該平臺對區(qū)域水土流失情況的變化進(jìn)行分析和管理。第二是農(nóng)戶，他們可以利用該系統(tǒng)平臺，了解自家的地塊水土流失的狀況，在耕作時應(yīng)該采取哪些措施更有利于達(dá)到保水保肥、優(yōu)質(zhì)高產(chǎn)的目的。針對不同的用戶，操作系統(tǒng)側(cè)重不同，但盡量簡單化，采取觸摸點擊的的方式操作。輸出的結(jié)果也盡量簡單易懂。

4.3.1 管理者應(yīng)用平臺 該系統(tǒng)結(jié)合不同地塊的坡長、坡度、土地利用方式、土壤性狀和種植管理等具體信息，進(jìn)行土壤侵蝕分類分級。土壤侵蝕強度分級，系統(tǒng)以年平均侵蝕模數(shù)為判別指標(biāo)，在缺少侵蝕模數(shù)實測資料時，根據(jù)有關(guān)侵蝕方式(面蝕、溝蝕)的判別指標(biāo)進(jìn)行[25]。該區(qū)域共形成了34種不同地形地貌、侵蝕強度、土壤養(yǎng)分等組合方案，當(dāng)管理者用戶選中地塊時，即會在預(yù)警提示窗口中顯示出其侵蝕模數(shù)以及侵蝕強度等具體信息，讓使用者直觀地感受到耕種地塊所處水土流失危害的程度。例如，當(dāng)管理者要對分類號為27的地塊進(jìn)行分析和管理時，點擊地塊分類號為27的地塊，系統(tǒng)即會提示此地塊平均坡度為3.78°，坡度范圍3°～4°，侵蝕模數(shù)為44.3 t/(hm2·a)，介于36～48 t/(hm2·a)，屬于土壤侵蝕極強烈區(qū)域，達(dá)到嚴(yán)重預(yù)警。系統(tǒng)同時會給出相應(yīng)的治理措施如工程措施：修筑坡狀梯田，坡頂種植水保林，林側(cè)挖節(jié)流溝；耕作措施：沿坡狀梯田橫坡起壟、秋起壟、秋雙層深施肥、垅底深松，另增加根茬、有機質(zhì)還田等措施的綜合治理模式；生物措施：沿坡每隔50 m修一條寬1～2 m的生物隔離帶，可選種多年生牧草、多年生經(jīng)濟(jì)植物(如石刁柏、胡枝子等)。

對分類號為3的地塊進(jìn)行分析和管理時，如上例操作，系統(tǒng)會提示該地塊坡度范圍為2～3°，屬中度侵蝕區(qū)，極易成為侵蝕嚴(yán)重區(qū)域，提示具體的治理措施：工程措施需采取農(nóng)田防護(hù)林體系和耕作措施，按國家和省級防護(hù)林體系標(biāo)準(zhǔn)營造農(nóng)田防護(hù)林，減輕風(fēng)蝕、水蝕的侵害；耕作措施中主要采取橫坡起壟、秋起壟、秋雙層深施肥、垅底深松，另增加根茬、有機質(zhì)還田等措施的綜合治理模式；暫時還不需要生物措施。

對于侵蝕溝的治理，用戶需選中第33類地塊，系統(tǒng)會提示該地塊為侵蝕溝區(qū)，并提出預(yù)警。系統(tǒng)針對溝沿、溝底和溝坡3個重點部位提出治理措施，如溝沿治理中包括沿溝造林或者播撒一些根系發(fā)達(dá)，固土能力強的灌木或草本植物，并禁止放牧、取土等人為活動，以防止人為破壞；溝底治理中有修筑谷坊，之后種植根系發(fā)達(dá)、根蘗力強、枝葉繁茂、耐水濕的樹種等措施；溝坡治理中包括造林、修建梯田等措施?？傊?，針對34種不同類別的土壤侵蝕區(qū)，系統(tǒng)都分別組合了以往專家的知識和水土流失治理的標(biāo)準(zhǔn)方案，可以實現(xiàn)針對不同類別進(jìn)行分析和綜合治理。

4.3.2 農(nóng)戶應(yīng)用平臺 農(nóng)戶的應(yīng)用主要是獲得地塊土壤侵蝕的信息，綜合治理的決策方案和指導(dǎo)，達(dá)到保水保肥，提高土壤肥力和作物產(chǎn)量的目的。為了方便農(nóng)民用戶使用，我們設(shè)計了兩種方式進(jìn)入系統(tǒng)：(1)農(nóng)戶直接在圖上找到自家地塊點擊進(jìn)入，得到自家地塊的詳細(xì)信息； (2)由于農(nóng)民的知識和對計算機的應(yīng)用有限，一般的漢語拼音都難以輸入，所以為了方便農(nóng)戶的使用，我們對每個農(nóng)戶承包的地塊都進(jìn)行了編號，這樣農(nóng)戶在使用時輸入自家的編號即可進(jìn)入自家地塊，獲得相應(yīng)的地塊信息。主要包括：面積、土壤養(yǎng)分狀況、肥力等級、地塊坡長、平均坡度、前幾年種植信息、產(chǎn)量等信息。并得到根據(jù)農(nóng)戶自家的自然條件，地形地貌，土壤養(yǎng)分的作物種植指導(dǎo)，如坡度較大的區(qū)域，系統(tǒng)提示改順坡壟為橫坡壟，以及如何起壟，是否還應(yīng)該種植生物隔離帶，施肥的方案選擇，增施有機肥和化肥的用量以及施肥方法等。如點擊農(nóng)戶高大海家，系統(tǒng)會自動彈出地塊信息、養(yǎng)分信息和建議概況3個窗體。養(yǎng)分信息窗口中提示高大海家地塊堿解氮含量偏高，有效磷含量低，速效鉀含量低；地塊信息窗口提示農(nóng)戶高大海家地塊的主要信息，包括地類代碼、農(nóng)戶姓名、坡度等；在建議概況窗口中，針對種植作物大豆，有以下推薦施肥方案：建議施尿素30～31.5 kg/hm2，二胺108 kg/hm2，50%硫酸鉀75 kg/hm2。之后提示施肥方法：氮肥、磷肥和鉀肥作為底肥全部一次性施入。

4.4 系統(tǒng)管理與更新

根據(jù)使用功能的不同，只有管理者和設(shè)計者可以進(jìn)行系統(tǒng)的更新。通過進(jìn)入系統(tǒng)后臺，對系統(tǒng)進(jìn)行數(shù)據(jù)維護(hù)，以及對數(shù)據(jù)進(jìn)行編輯和更新，并且可以進(jìn)行治理方案更新和升級，重新制定，根據(jù)不同要求對專題圖的制作以及進(jìn)行一些空間分析工作的更新等。

5 結(jié) 論

應(yīng)用結(jié)果表明，該系統(tǒng)平臺的應(yīng)用，為管理者提供了到農(nóng)戶地塊的，可以定量分析的水土流失綜合治理分析平臺，促進(jìn)了農(nóng)村信息化，農(nóng)業(yè)現(xiàn)代化管理的進(jìn)程。通過系統(tǒng)平臺的工程指導(dǎo)、生物指導(dǎo)和耕作種植指導(dǎo)，使得水土流失的綜合治理取得了很大的成效，農(nóng)戶在應(yīng)用中方便快捷，在保水保肥提高作物產(chǎn)量上有一定程度的改進(jìn)。

(1)能將土壤侵蝕和土壤養(yǎng)分信息綜合在一起，為農(nóng)戶耕作提供決策方案，提高了土地保水保肥的能力。

(2)利用遙感數(shù)據(jù)和實地調(diào)查數(shù)據(jù)相結(jié)合，在GIS平臺上生成以地塊為基本單元的侵蝕治理單元圖。方便管理者和農(nóng)戶使用。

(3)充分利用當(dāng)?shù)亻L期田間定位試驗研究成果和歷年收集數(shù)據(jù)，結(jié)合張憲奎等[8]提出的黑龍江省土壤侵蝕方程，建立土壤侵蝕模型，非常符合當(dāng)?shù)貙嶋H情況，實用性很高。

(4)采用GIS平臺與可視化編程語言二者相集成的開發(fā)方式，具有較好的可移植性、可擴(kuò)展性；調(diào)用空間數(shù)據(jù)庫和土壤侵蝕模型庫，實現(xiàn)針對不同地塊不同侵蝕強度，提出相應(yīng)的治理措施。用戶界面設(shè)計簡潔、友好，易于操作。

[參考文獻(xiàn)]

[1]孟令欽,李勇.東北黑土區(qū)坡耕地侵蝕溝發(fā)育機理初探[J].水土保持學(xué)報,2009,23(1):7-11.

[2]張曉平,梁愛珍,申艷,等.東北黑土水土流失特點[J].地理科學(xué),2006,26(6):687-692.

[3]范昊明,蔡強國,崔明.東北黑土漫崗區(qū)土壤侵蝕垂直分帶性研究[J].農(nóng)業(yè)工程學(xué)報,2005,21(6):8-11.

[4]范昊明,蔡強國,王紅閃.中國東北黑土區(qū)土壤侵蝕環(huán)境[J].水土保持學(xué)報,2004,18(2):66-701.

[5]崔明,張旭東,蔡強國,等.東北典型黑土區(qū)氣候、地貌演化與黑土發(fā)育關(guān)系[J].地理研究,2008,27(3):527-535.

[6]梁愛珍,張曉平,楊學(xué)明,等.東北黑土有機碳的分布及其損失量研究[J].土壤通報,2008,39(3):533-538.

[7]劉興土,閻百興.東北黑土區(qū)水土流失與糧食安全[J].中國水土保持,2009(1):17-19.

[8]王小兵,吳元元,鄧玲.東北黑土區(qū)黑土退化防治與保護(hù)研究[J].資源與產(chǎn)業(yè),2008,10(3):81-83.

[9]崔明,蔡國強,范昊明.東北黑土區(qū)土壤侵蝕研究進(jìn)展[J].水土保持研究,2007,14(5):29-34.

[10]閆業(yè)超,張樹文,岳書平.近40 a黑土典型區(qū)坡溝侵蝕動態(tài)變化[J].農(nóng)業(yè)工程學(xué)報,2010,26(2):109-115.

[11]魏茂盛.移動GIS 在全國第一次水利普查中的應(yīng)用[J].水利信息化,2011,8(4):48-50.

[12]張巖,劉憲春,李智廣,等.利用侵蝕模型普查黃土高原土壤侵蝕狀況[J].農(nóng)業(yè)工程學(xué)報, 2012,28(10):165-171.

[13]盧敏,張展羽.GIS在農(nóng)業(yè)水土領(lǐng)域的應(yīng)用與研究進(jìn)展[J].中國農(nóng)業(yè)水利水電, 2004(9):8-11.

[14]李勇,趙軍.針對農(nóng)戶地塊的施肥決策支持系統(tǒng)的設(shè)計與實現(xiàn)[J].農(nóng)業(yè)工程學(xué)報,2010,26(1):192-196.

[15]陳云明,劉國彬,鄭粉莉,等.RUSLE侵蝕模型的應(yīng)用及進(jìn)展[J].水土保持研究,2004,11(4):80-83.

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

[17]史志華,蔡崇法,丁樹文.基于GIS和RUSLE的小流域農(nóng)地水土保持規(guī)劃研究[J].農(nóng)業(yè)工程學(xué)報,2002,18(4):172-175.

[18]張憲奎,許靖華.黑龍江省土壤流失方程的研究[J].水土保持通報,1992,12(4):1-9.

[19]杜朝正.基于ArcGIS的坡度分析[J].實驗與技術(shù),2009,25(1):17-18.

[20]劉和平,王秀穎,劉寶元.短坡條件下侵蝕產(chǎn)沙與坡長的關(guān)系[J],水土保持學(xué)報,2011,5(2):1-5.

[21]王秀穎,劉和平,劉寶元.變雨強人工降雨條件下坡長對徑流的影響研究[J].水土保持學(xué)報,2010,24(6):1-5.

[22]章文波,劉寶元.基于GIS的中國土壤侵蝕預(yù)報信息系統(tǒng)[J].水土保持學(xué)報,2003,17(2):89-92.

[23]胡剛,伍永秋,劉寶元,等.GPS和GIS進(jìn)行短期溝蝕研究初探:以東北漫川漫崗黑土區(qū)為例[J].水土保持學(xué)報,2004,18(4):16-19.

[24]陳云明,劉國彬,鄭粉莉,等.RUSLE侵蝕模型的應(yīng)用及進(jìn)展[J].水土保持研究,2004,11(4):80-83.

[25]中華人民共和國水利部.SL446—2009黑土區(qū)水土流失綜合防治技術(shù)標(biāo)準(zhǔn)[S].北京：中國水利水電出版社,2009.