溫莉 周廷剛 劉曉璐 吳浪

摘要 [目的]研究更簡(jiǎn)便、更精確的模型計(jì)算重慶市的水土流失量。[方法]以RS、GIS、eCognition、SPSS等軟件為平臺(tái)，將基于RUSLE模型計(jì)算得到的2002年重慶市各區(qū)、縣的年水土流失量數(shù)據(jù)與官方公布數(shù)據(jù)進(jìn)行回歸分析，建立重慶市土壤侵蝕模型。[結(jié)果]基于RUSLE模型計(jì)算得到的2002年重慶市各區(qū)、縣的年水土流失量數(shù)據(jù)相對(duì)于官方公布數(shù)據(jù)，平均精度達(dá)72.38%，在多緩坡的區(qū)、縣精度較高，在多陡坡的區(qū)、縣精度較低。由建立的重慶市土壤侵蝕模型計(jì)算得到的2002年重慶市10個(gè)區(qū)、縣的年水土流失量數(shù)據(jù)相對(duì)于官方公布數(shù)據(jù)，模擬精度達(dá)85.73%，較好地改善了RUSLE模型在陡坡地區(qū)精度較差的不足。重慶市土壤侵蝕模型相對(duì)于RUSLE模型，在重慶市10個(gè)驗(yàn)證數(shù)據(jù)的水土流失量計(jì)算中，計(jì)算結(jié)果的平均精度提高了10.6%。[結(jié)論] 重慶市土壤侵蝕模型相對(duì)于RUSLE模型在重慶市水土流失量的計(jì)算中具有更好的實(shí)用性，可為重慶市水土流失量的計(jì)算提供參考。

關(guān)鍵詞 模型；RS；GIS；；RUSLE； 水土流失；重慶市

中圖分類號(hào) S157 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 0517-6611（2018）27-0130-04

The Study of Soil Erosion Models in Chongqing Based on RS and GIS

WEN Li1，ZHOU Tinggang1，2，LIU Xiaolu1 et al

（1.School of Geographical Sciences，Southwest University，Chongqing 400715；2.Chongqing Engineering Research Center for Remote Sensing Big Data Application，Chongqing 400715）

Abstract [Objective] To explore a simple and more accurate model to calculate the amount of soil and water loss in Chongqing.[Method] Based on the RUSLE model，this paper took advantage of RS，GIS，eCognition，SPSS and other software.First of all，six factors that affect soil and water loss were calculated and the distribution of soil and water loss of Chongqing in 2002 was obtained by applying the overlay analysis function of GIS.Then，calculating the amount of soil and water loss in various districts and counties of Chongqing in 2002 and making regression analysis with official published data，finally established the model of soil erosion in Chongqing.[Result]With the standard of official data ，the average accuracy of soil and water loss in districts and counties of Chongqing in 2002 which was calculated based on the RUSLE model reached to 72.38%.The accuracy was more accurate in districts which had more gentle slopes，and was poorer in steep districts.With the standard of official data，the simulation accuracy of soil and water loss in districts and counties of Chongqing in 2002 which was calculated based on the established model of soil erosion in Chongqing，reached to 85.73%.The established model of soil erosion in Chongqing improved the deficiencies of RUSLE model in steep lands.Compared with the RUSLE model，the established soil erosion model in Chongqing showed that the average accuracy of calculated results was improved by 10.6% in calculation the amount of soil and water loss in 10 validated data in Chongqing.[Conclusion] Compared with the RUSLE model，the model of soil erosion in Chongqing improves the calculation precision in the areas which have large gradient，and has better practicability in the calculation of soil and water loss in Chongqing.

Key words Model；RS；GIS；RUSLE；Soil and water loss；Chongqing

基金項(xiàng)目 國(guó)務(wù)院三峽辦三峽后續(xù)工作庫(kù)區(qū)生態(tài)與生物多樣性保護(hù)專項(xiàng)“重慶庫(kù)區(qū)重要支流（涪陵—巫山段）水生生境狀況調(diào)查與評(píng)估”（500000201BB5200002）。

作者簡(jiǎn)介 溫莉（1993—），女，四川瀘州人，碩士研究生，研究方向：遙感與地理信息系統(tǒng)應(yīng)用。*通訊作者，教授，博士，從事遙感與地理信息系統(tǒng)應(yīng)用研究。

收稿日期 2018-05-21

水土流失是我國(guó)面臨的三大生態(tài)環(huán)境問題之一[1]。全國(guó)水土流失面積占國(guó)土面積的50.49%[2]。水土流失會(huì)造成土壤肥力下降，大量的徑流泥沙會(huì)導(dǎo)致水庫(kù)淤積、河道阻塞等問題[3]。我國(guó)地域廣闊，自然和經(jīng)濟(jì)社會(huì)條件復(fù)雜，水土流失分布廣、面積大、類型多，區(qū)域水土流失的監(jiān)測(cè)和治理存在較大差異[4]。針對(duì)不同區(qū)域的自然和經(jīng)濟(jì)社會(huì)狀況，提出適應(yīng)于當(dāng)?shù)氐耐寥懒魇ьA(yù)測(cè)模型，是提高土壤流失預(yù)測(cè)預(yù)報(bào)精度的關(guān)鍵，也是防治水土流失的前提。

RUSLE土壤流失模型是一個(gè)通用的經(jīng)驗(yàn)?zāi)Ｐ?，適用于緩坡耕地土壤流失量的計(jì)算。對(duì)于我國(guó)大部分地區(qū)，用RUSLE模型計(jì)算土壤流失量仍有較大誤差。因此，國(guó)內(nèi)學(xué)者以通用土壤方程為基礎(chǔ)，建立了適合于當(dāng)?shù)氐耐寥狼治g預(yù)測(cè)模型。張憲奎等[5]以RUSLE土壤流失模型為基礎(chǔ)，通過對(duì)7年大量試驗(yàn)數(shù)據(jù)的統(tǒng)計(jì)分析，對(duì)模型中各因子的求算方法和數(shù)值進(jìn)行全面的修改，得到了適用于黑龍江省的土壤流失方程。劉斌濤等[6]針對(duì)蘆山地震災(zāi)區(qū)的特殊情況，建立了蘆山地震災(zāi)區(qū)土壤流失模型（LSSLE）。林素蘭等[7]以USLE土壤流失模型為基礎(chǔ)，通過11年的徑流小區(qū)試驗(yàn)，得到了遼北低山丘陵區(qū)坡耕地的土壤侵蝕方程，為該地區(qū)坡耕地治理提供了科學(xué)依據(jù)。筆者以重慶市作為研究區(qū)，以RUSLE土壤流失模型為基礎(chǔ)，計(jì)算2002年重慶市各區(qū)、縣的水土流失量，將基于RUSLE模型計(jì)算得到的2002年重慶市各區(qū)、縣的年水土流失量與官方公布數(shù)據(jù)進(jìn)行回歸分析與驗(yàn)證，建立重慶市土壤侵蝕模型。該模型的建立，可以及時(shí)有效地預(yù)測(cè)、分析重慶市水土流失現(xiàn)狀，為開展水土流失治理工作和政府及有關(guān)部門的決策提供科學(xué)依據(jù)。

1 研究區(qū)概況與研究方法

1.1 研究區(qū)概況 重慶市位于我國(guó)西南部，長(zhǎng)江上游地區(qū)，位于105°11′～110°11′E、28°10′～32°13′N，是青藏高原與長(zhǎng)江中下游平原的過渡地帶，幅員面積8.24萬km2。地貌以丘陵、山地為主，坡地面積較大，有“山城”之稱。重慶年平均氣溫16～18 ℃，年平均降水量較豐富，大部分地區(qū)在1 000～1 350 mm，降水集中在5—9月，占全年總降水量的70%左右。土壤以水稻土和紫色土為主。易發(fā)生水土流失的自然條件，加之不合理的土地利用，使重慶市的水土流失十分嚴(yán)重。

1.2 研究方法

1.2.1 數(shù)據(jù)來源。

重慶市及其周邊鄰省共35個(gè)站點(diǎn)的1992—2002年日降水量數(shù)據(jù)來源于中國(guó)氣象數(shù)據(jù)共享網(wǎng)（http：//data.cma.cn/）；遙感影像來源于https：//glovis.usgs.gov/；DEM數(shù)據(jù)來源于http：//www.gscloud.cn/，選用的是GDEMV2 30 m空間分辨率的產(chǎn)品；土壤數(shù)據(jù)來源于中國(guó)科學(xué)院南京土壤研究所；重慶市水土流失數(shù)據(jù)來源于重慶市水利局。

1.2.2 降雨侵蝕力因子。

降雨侵蝕力（R）為單位降雨作用地表使其產(chǎn)生水土流的能量[8]，是降雨引起土壤侵蝕的潛在能力，它從定量的角度評(píng)價(jià)了降雨對(duì)土壤顆粒分離和搬運(yùn)的作用，是通用土壤流失方程（USLE）及其修正模型（RUSLE）中六大基本因子之一[9]。RUSLE模型中計(jì)算降雨侵蝕力的經(jīng)典算法所需要的數(shù)據(jù)很難獲取且計(jì)算復(fù)雜，因此該研究采用盧喜平等[10]提出的重慶地區(qū)降雨侵蝕力計(jì)算模型求解R值，采用的數(shù)據(jù)是1992—2002年的降水?dāng)?shù)據(jù)。

R=5.249×F1.205

F= 12 i=1 pi p ×pi （1）

式中，R為年降雨侵蝕力，（MJ ·mm）/（hm2·h·a）；p為年平均降雨量，mm；pi為各月平均降雨量，mm。

1.2.3 植被覆蓋度因子。植被覆蓋度因子（C）指在相同的土壤、坡度和降雨條件下，有特定植被覆蓋或田間管理土地上的土壤流失量與實(shí)施清耕、無覆蓋裸露休閑地上的土壤流失量之比，其值在0～1，無量綱[11]。以2002年的TM影像計(jì)算NDVI，計(jì)算公式：

NDVI=（ρNIR+ρRED）/（ρNIR+ρRED） （2）

式中，NDVI為歸一化植被指數(shù)；ρNIR為近紅外波段；ρRED為可見光紅光波段。

在NDVI的基礎(chǔ)上，應(yīng)用公式（3）直接計(jì)算植被覆蓋度因子。在大尺度區(qū)域利用遙感植被指數(shù)（NDVI）計(jì)算 C 值是最有效、最實(shí)用的方法[12]，其計(jì)算公式：

C=exp -2× NDVI 1-NDVI （3）

1.2.4 坡度坡長(zhǎng)因子。

坡度坡長(zhǎng)因子（SL因子）用來衡量地形對(duì)土壤侵蝕的影響?；跂鸥竦臄?shù)字高程模型（DEM）可以實(shí)現(xiàn)基于像元的坡度和坡長(zhǎng)因子的計(jì)算。首先基于DEM生成重慶市坡度圖，然后采用卜兆宏等[13]根據(jù)實(shí)測(cè)數(shù)據(jù)得出了適用于我國(guó)S因子的計(jì)算方法，其計(jì)算公式：

S=0.621 1×203.1576×sinθ （4）

式中，S為坡度因子，無綱量；θ是坡度，單位是弧度。

基于DEM數(shù)據(jù)，根據(jù)Wischmeier等[14]在RUSLE模型中的坡長(zhǎng)因子（L）計(jì)算方法[10]，其計(jì)算公式：

L=（λ/22.13）α （5）

α=β/（β+1） （6）

β=（sinθ/0.0896）/[3×（sinθ）0.8+0.56] （7）

λ=D/cosθ

式中，L為坡長(zhǎng)因子；λ為水平坡長(zhǎng)；α為坡長(zhǎng)指數(shù)；β代表細(xì)溝侵蝕與細(xì)溝間侵蝕的比率；θ是坡度，D為像元尺度坡長(zhǎng)。

1.2.5 土壤可蝕性因子。

土壤可蝕性因子是通用土壤流失方程系列模型（USLE、RUSLE、CSLE）中一個(gè)必要參數(shù)，常用K表示[15]。以重慶市土壤分布圖為底圖，根據(jù)重慶市水土保持局官方公布的重慶市土類土壤可蝕性因子K值對(duì)底圖進(jìn)行賦值。其賦值見表1。

1.2.6 水土保持因子。

水土保持因子（P）指采取一定措施后土壤流失量與順坡種植時(shí)土壤流失量的比值[16]。對(duì)于因子P的計(jì)算目前尚無定量公式，一般選取經(jīng)驗(yàn)值，P取0～1，0代表基本上不發(fā)生侵蝕的地區(qū)；1代表未采取任何控制措施的地區(qū)[17]。

將2002年的Landsat-5 TM影像進(jìn)行預(yù)處理和監(jiān)督分類，并將重慶市土地利用類型分為林地、草地、耕地、工礦和建筑用地、裸地、水體6類，得到土地利用類型圖。參考美國(guó)農(nóng)業(yè)部手冊(cè)703號(hào)和重慶市水土保持局官方公布的重慶水土保持因子P值及其相關(guān)文獻(xiàn)對(duì)重慶市2002年土地利用類型圖的各土地利用類型進(jìn)行賦值。其賦值：林地1.00，草地1.00，耕地 0.03，工礦和建筑用地 0，水體 0，裸地 1.00。

1.2.7 重慶市年土壤侵蝕量（A）的估算。

年土壤侵蝕量A的估算采用RUSLE模型，其計(jì)算公式

A=224.2×R×K×S×L×C×P （8）

式中，A為年平均土壤流失量，t/（km2 ·a）。

1.2.8 重慶市土壤侵蝕模型的建立。

把基于RUSLE模型得到的2002年重慶市38個(gè)區(qū)、縣的年土壤侵蝕總量數(shù)據(jù)分為試驗(yàn)數(shù)據(jù)和驗(yàn)證數(shù)據(jù)。選取28個(gè)數(shù)據(jù)為試驗(yàn)數(shù)據(jù)，以官方數(shù)據(jù)（官方數(shù)據(jù)來自重慶市水利局公布的水土保持公報(bào)）建立擬合模型，選取10個(gè)數(shù)據(jù)為驗(yàn)證數(shù)據(jù)，并通過驗(yàn)證數(shù)據(jù)對(duì)其進(jìn)行精度評(píng)價(jià)。建立的重慶市土壤侵蝕模型：

A=λ×224.2×R×K×S×L×C×P+φ （9）

式中，λ、φ為模型參數(shù)。

安徽農(nóng)業(yè)科學(xué) 2018年

2 結(jié)果與分析

2.1 基于RUSLE模型的土壤侵蝕總量數(shù)據(jù)與官方數(shù)據(jù)的對(duì)比

為了檢驗(yàn)該研究計(jì)算方法的可靠性和探究RUSLE模型的不足，將基于RUSLE模型的土壤侵蝕總量數(shù)據(jù)與官方數(shù)據(jù)進(jìn)行對(duì)比（官方數(shù)據(jù)中各因子的計(jì)算方法與該研究的計(jì)算方法不一致），結(jié)果表明，官方數(shù)據(jù)與基于RUSLE模型的土壤侵蝕總量數(shù)據(jù)的一致性較高，二者一致性變化為45.54%～89.84%，平均精度達(dá)72.38%。表明該研究計(jì)算水土流失各影響因子的方法可靠，計(jì)算結(jié)果精度較高，可以在重慶市使用。土壤侵蝕總量在多緩坡的區(qū)、縣誤差較小，誤差較大的主要出現(xiàn)在奉節(jié)縣、城口縣、云陽(yáng)縣、開州區(qū)、南川區(qū)、江津區(qū)等地勢(shì)崎嶇的區(qū)、縣，主要原因是RUSLE模型適合緩坡耕地水土流失量的計(jì)算，對(duì)于坡度較大的地區(qū)，計(jì)算結(jié)果精度較差。

2.2 重慶市土壤侵蝕模型的建立

重慶市多山地丘陵，地勢(shì)崎嶇，RUSLE模型不完全適合重慶市水土流失量的計(jì)算。為了改善RUSLE模型在陡坡地區(qū)計(jì)算精度的不足，減小數(shù)據(jù)的計(jì)算誤差，把基于RUSLE模型得到的沙坪壩區(qū)、九龍坡區(qū)、酉陽(yáng)縣、奉節(jié)縣、合川區(qū)、江津區(qū)、長(zhǎng)壽區(qū)、云陽(yáng)縣、巫山縣、江北區(qū)的年土壤侵蝕總量數(shù)據(jù)作為驗(yàn)證數(shù)據(jù)，其他28個(gè)區(qū)、縣的年土壤侵蝕總量數(shù)據(jù)作為試驗(yàn)數(shù)據(jù)，建立重慶市土壤侵蝕模型。將試驗(yàn)數(shù)據(jù)與重慶市水利局2002年水土保持公報(bào)的官方數(shù)據(jù)進(jìn)行回歸分析，判定系數(shù)R2為0.752，擬合度較高，不被解釋的變量較少?；貧w方程顯著性檢驗(yàn)的概率 為0，小于顯著性水平0.05，則認(rèn)為系數(shù)不同時(shí)為0時(shí)，被解

釋變量與解釋變量全體的線性關(guān)系顯著，可建立線性方程。線性回歸模型的模型參數(shù)λ=0.826，φ=9.211，擬合線性方程為Y=0.826X+9.211（圖1）。由此可得到重慶市土壤侵蝕模型為A=0.826×224.2×R×K×S×L×C×P+9.211。

2.3 重慶市土壤侵蝕模型的驗(yàn)證 由該研究得到的重慶市土壤侵蝕模型計(jì)算年土壤侵蝕量A的方法相對(duì)簡(jiǎn)單，數(shù)據(jù)資料獲取相對(duì)便捷，但模型存在地域性，必須經(jīng)過驗(yàn)證，才可在

重慶市使用。將重慶市官方數(shù)據(jù)分別與應(yīng)用重慶市土壤侵蝕模型計(jì)算得到的沙坪壩區(qū)、九龍坡區(qū)、酉陽(yáng)縣、奉節(jié)縣、合川區(qū)、江津 區(qū)、長(zhǎng)壽區(qū)、云陽(yáng)縣、巫山縣、江北區(qū)的年土壤侵蝕總量數(shù)據(jù)和基于RUSLE模型計(jì)算得到的以上10個(gè)區(qū)、縣的土壤侵蝕總量數(shù)據(jù)進(jìn)行相關(guān)性分析，相關(guān)度分別為0.842和0.721，官方數(shù)據(jù)與基于重慶市土壤侵蝕模型的土壤侵蝕總量數(shù)據(jù)的相關(guān)性比官方數(shù)據(jù)與基于RUSLE模型的土壤侵蝕總量數(shù)據(jù)的相關(guān)性提高了0.121。對(duì)三者進(jìn)行誤差分析，結(jié)果表明，官方數(shù)據(jù)與基于重慶市土壤侵蝕模型的土壤侵蝕總量數(shù)據(jù)一致性為65.25%～88.99%，平均精度達(dá)85.73%。官方數(shù)據(jù)與基于RUSLE模型的土壤侵蝕總量數(shù)據(jù)的一致性為50.51%～81.63%，平均精度為75.13%，平均精度提高了10.6%（圖2）。將官方數(shù)據(jù)與基于重慶市土壤侵蝕模型的土壤侵蝕總量數(shù)據(jù)進(jìn)行線性擬合，決定性系數(shù)R2為0.870（圖3），官方數(shù)據(jù)與基于RUSLE模型的土壤侵蝕總量數(shù)據(jù)進(jìn)行線性擬合，決定性系數(shù)R2為0.730（圖4），表明官方數(shù)據(jù)與基于重慶市土壤侵蝕模型的土壤侵蝕總量數(shù)據(jù)的擬合度較好，2組數(shù)據(jù)無顯著差異，基于重慶市土壤侵蝕模型的土壤侵蝕總量數(shù)據(jù)比基于RUSLE模型的土壤侵蝕總量數(shù)據(jù)更精確、可靠。

3 結(jié)論

該研究以RUSLE方程為基礎(chǔ)，采用更易獲取的數(shù)據(jù)和更簡(jiǎn)單的方法來計(jì)算影響重慶市水土流失各因子。將基于RUSLE方程計(jì)算得到的2002年重慶市各區(qū)、縣的年水土流失量與官方公布數(shù)據(jù)進(jìn)行回歸分析與驗(yàn)證，建立重慶市土壤侵蝕模型，結(jié)果表明，基于RUSLE模型計(jì)算得到的2002年重慶市各區(qū)、縣的年水土流失量相對(duì)于官方公布數(shù)據(jù)，模擬精度達(dá)72.38%，數(shù)據(jù)一致性較好。但在奉節(jié)縣、城口縣、云陽(yáng)縣、開州區(qū)、南川區(qū)、江津區(qū)等地勢(shì)較崎嶇的區(qū)、縣誤差較大。由建立的重慶市土壤侵蝕模型計(jì)算得到的2002年重慶市10個(gè)區(qū)、縣的年水土流失量數(shù)據(jù)相對(duì)于官方公布數(shù)據(jù)，平均精度達(dá)85.73%，與基于RUSLE模型的土壤侵蝕總量數(shù)據(jù)相比，平均精度提高了10.6%，減輕了地形對(duì)水土流失量的影

響，提高了陡坡地區(qū)水土流失量的計(jì)算精度。

參考文獻(xiàn)

[1]

王剛，張秋平，鄭海金，等.1987-2013年江西省水土流失趨勢(shì)及其社會(huì)經(jīng)濟(jì)驅(qū)動(dòng)力分析[J].生態(tài)科學(xué)，2017，36（3）：115-120.

[2] 錢登峰，莊曉暉，張博.高寒區(qū)凍融侵蝕類型及驅(qū)動(dòng)力分析[J].中國(guó)水土保持，2014（6）：16-17.

[3] 田衛(wèi)堂，胡維銀，李軍，等.我國(guó)水土流失現(xiàn)狀和防治對(duì)策分析[J].水土保持研究，2008，15（4）：204-209.

[4] 王治國(guó)，張超，紀(jì)強(qiáng)，等.全國(guó)水土保持區(qū)劃及其應(yīng)用[J].中國(guó)水土保持科學(xué)，2016，14（6）：101-106.

[5] 張憲奎，徐靖華，盧秀琴，等.黑龍江省土壤流失方程的研究[J].水土保持通報(bào)，1992，12（4）：1-9，18.

[6] 劉斌濤，宋春風(fēng)，史展，等.蘆山地震災(zāi)區(qū)土壤流失方程研究[J].長(zhǎng)江科學(xué)院院報(bào)，2016，33（1）：15-19.

[7] 林素蘭，黃毅，聶振剛，等.遼北低山丘陵區(qū)坡耕地土壤流失方程的建立[J].土壤通報(bào)，1997，28（6）：251-253.

[8] ELTAIF N I，GHARAIBEH M A，ALZAITAWI F，et al.Approximation of rainfall erosivity factors in North Jordan[J].Pedosphere，2010，20（6）：711-717.

[9] 周琳，許武成，尹義星.近61年四川省降雨侵蝕力的時(shí)空變化趨勢(shì)[J].水土保持研究，2018，25（4）：8-14.

[10] 盧喜平. 紫色土丘陵區(qū)降雨侵蝕力模擬研究[D].重慶：西南大學(xué)，2006.

[11] 張巖，劉寶元，史培軍，等.黃土高原土壤侵蝕作物覆蓋因子計(jì)算[J].生態(tài)學(xué)報(bào)，2001，21（7）：1050-1056.

[12] 陳學(xué)兄.基于遙感與GIS的中國(guó)水土流失定量評(píng)價(jià)[D].楊凌：西北農(nóng)林科技大學(xué)，2013：89.

[13]卜兆宏，唐萬龍，潘賢章.土壤流失量遙感監(jiān)測(cè)中GIS像元地形因子算法的研究[J].土壤學(xué)報(bào)，1994，31（3）：322-329.

[14] WISCHMEIER W H，SMITH D D.Predictingrainfall erosion losses：A guide to conservation planning[M]. Washington，DC：Science and Education Administration，United States Department of Agriculture， 1978.

[15] OSTOVARI Y.Modification of the USLE K factor for soil erodibility assessment on calcareous soils in Iran[J].Geomorphology，2016，273：385-395.

[16] PRASANNAKUMAR V，VIJITH H，ABINOD S，et al.Estimation of soil erosion risk within a small mountainous sub-watershed in Kerala，India，using Revised Universal Soil Loss Equation（RUSLE） and geo-information technology[J].Geoscience frontiers，2012，3（2）：209-215.

[17] 韓佳.基于RS與GIS的東川地區(qū)水土流失研究[D].北京：中國(guó)地質(zhì)大學(xué)，2010：50.