李翠漫, 盧 遠(yuǎn), 劉斌濤, 許貴林

(1.廣西師范學(xué)院 地理科學(xué)與規(guī)劃學(xué)院, 南寧 530001; 2.中國(guó)科學(xué)院 水利部成都山地災(zāi)害與環(huán)境研究所, 成都 610041; 3.廣西師范學(xué)院 北部灣環(huán)境演變與資源利用教育部重點(diǎn)實(shí)驗(yàn)室, 南寧 530001)

土壤侵蝕導(dǎo)致大量土壤資源遭到破壞，使土層變薄，肥沃表土大量流失，地表植被遭到嚴(yán)重破壞，導(dǎo)致自然生態(tài)環(huán)境失調(diào)，嚴(yán)重威脅了人們的生產(chǎn)發(fā)展[1]。1990年珠江流域土壤侵蝕現(xiàn)狀調(diào)查顯示[2]，西江流域土壤侵蝕面積53 327.89 km2,占珠江流域土壤侵蝕面積的93.44%，占西江流域土地總面積的15.10%。2013年全國(guó)第一次全國(guó)水利普查公告，廣西土壤水力侵蝕面積5.05萬(wàn)km2，根據(jù)廣西水利廳頒布的廣西水土保持公報(bào)，截至2015年底，廣西累計(jì)完成水土流失治理面積3.34萬(wàn)km2。水土流失已成為西江流域重要生態(tài)問(wèn)題，研究西江流域土壤侵蝕空間分布，對(duì)廣西乃至整個(gè)珠江流域的可持續(xù)發(fā)展和生態(tài)治理環(huán)境治理提供決策參考，意義深遠(yuǎn)。

長(zhǎng)期以來(lái)土壤侵蝕一直是我國(guó)關(guān)注的熱點(diǎn)問(wèn)題，20世紀(jì)30年代Wischmeier等[3]提出的通用土壤流失方程(USLE)得到廣泛推廣和應(yīng)用，但原始的USLE模型是一個(gè)經(jīng)驗(yàn)?zāi)Ｐ?，主要適用于美國(guó)落基山脈以東的農(nóng)業(yè)區(qū)，在其他地區(qū)使用需要進(jìn)行參數(shù)校正。因此我國(guó)不少學(xué)者在研究區(qū)域土壤侵蝕問(wèn)題上，基于USLE模型基礎(chǔ)上進(jìn)行修正并運(yùn)用修正后的模型對(duì)區(qū)域土壤侵蝕定量評(píng)價(jià)，如劉寶元等[4]改進(jìn)建立了中國(guó)土壤流失方程(CSLE)，楊子生[5]開(kāi)發(fā)了云南金沙江流域土壤流失方程，江忠善等[6]建立了中國(guó)坡面水蝕預(yù)報(bào)模型等。

廣西水土保持監(jiān)測(cè)起步較晚，水土保持監(jiān)測(cè)網(wǎng)點(diǎn)少、基礎(chǔ)研究薄弱，難以直接利用定位監(jiān)測(cè)數(shù)據(jù)來(lái)對(duì)廣西西江流域的水土流失態(tài)勢(shì)進(jìn)行評(píng)價(jià)和分析，為了全面了解廣西水土流失狀況，2015年廣西開(kāi)展了土壤侵蝕調(diào)查與制圖工作，本文正是在此項(xiàng)目背景的支持下，基于USLE模型，對(duì)喀斯特分布廣泛的西江流域，引入石漠化修正因子M，建立西江流域土壤流失方程，定量研究土壤侵蝕的空間分布，并探索其水土流失的特征。

1　研究區(qū)概況

西江作為珠江流域的主流，是華南地區(qū)最長(zhǎng)的河流，為中國(guó)第三大河流，長(zhǎng)度僅次于長(zhǎng)江、黃河。全長(zhǎng)約2 214 km,集水面積約為35.3萬(wàn)km2[7]，西江流域橫跨廣西，在廣西境內(nèi)集水面積約為20.24萬(wàn)km2，占廣西土地總面積的85.7%。自西向東主要水系有南盤(pán)江、紅水河、柳江、黔潯江、郁江、右江、左江。

廣西西江流域喀斯特分布廣泛，碳酸鹽巖及碎屑巖為流域內(nèi)的主要地質(zhì)類型，總體地勢(shì)西北高東南低，地貌類型以山地為主，多為小起伏、中起伏山地，山地面積約占流域總面積的64.68%；流域內(nèi)地處亞熱帶氣候區(qū)，氣候溫暖多雨，雨量主要集中在春夏兩季，土壤種類繁多，以紅壤、赤紅壤及棕色石灰土為主，植被類型主要有針葉林、闊葉林、竹林、灌叢、草叢等。流域范圍內(nèi)人口密度大，人類活動(dòng)強(qiáng)烈，且喀斯特地區(qū)地表坎坷嶙峋、土層淺薄貧瘠，土地利用類型多樣，不合理的開(kāi)發(fā)和利用土地，則出現(xiàn)了一系列生態(tài)問(wèn)題，其中以生態(tài)破壞后所產(chǎn)生的石漠化問(wèn)題最為突出。

2　土壤流失方程建立

2.1　土壤流失方程修正

本研究以通用土壤流失方程(USLE)為藍(lán)本，結(jié)合廣西西江流域特點(diǎn)，引入石漠化修正因子M，建立廣西西江流域土壤侵蝕預(yù)報(bào)模型，進(jìn)行土壤侵蝕定量評(píng)估。方程計(jì)算公式為：

A=R×K×LS×C×P×M

(1)

式中：A為土壤侵蝕模數(shù)[t/(km2·a)]；R為降雨侵蝕力因子[MJ·mm/(hm2·h·a)]；K為土壤可蝕性因子[t·km2·h/(hm2·MJ·mm)]；LS為坡長(zhǎng)坡度因子(無(wú)量綱)；C為地表覆蓋與管理因子(無(wú)量綱)；P為水土保持措施因子(無(wú)量綱)；M為石漠化修正因子(無(wú)量綱)。

2.1.1降雨侵蝕力用EI30作為降雨侵蝕能力指標(biāo)需要連續(xù)記錄的降雨過(guò)程資料，在許多地區(qū)都難以獲取這些資料，因此本文利用觀測(cè)值和經(jīng)驗(yàn)建立降雨侵蝕力的簡(jiǎn)易算法，即通過(guò)氣象站收集到的廣西逐日降雨資料(1980—2015年)，采用章文波等[8-9]基于日降雨資料的降雨侵蝕力計(jì)算模型，計(jì)算得出廣西西江流域平均降雨侵蝕力為：9 585.40 MJ·mm/(hm2·h·a)。

2.1.2土壤可蝕性土壤可蝕性估算模型較多，其中EPIC模型[10]在我國(guó)是應(yīng)用較多的土壤可蝕性計(jì)算模型。張科利等[11]研究指出EPIC模型計(jì)算的K值與中國(guó)各地區(qū)實(shí)測(cè)的K值相差較大，并提出K值修正公式為K=-0.01383+0.51575KEPIC。因此本文土壤可蝕性采用EPIC模型并結(jié)合張科利提出修正K值的公式進(jìn)行估算。

2.1.3地形因子地形是決定土壤侵蝕的重要內(nèi)在因素，在土壤流失方程中地形因素采用地形因子表示，地形因子(LS)即坡長(zhǎng)因子(L)和坡度因子(S)。西南土石山區(qū)的地形特征極為破碎，土壤侵蝕環(huán)境比較復(fù)雜、土壤侵蝕空間異質(zhì)性大。因此本文地形因子計(jì)算采用劉斌濤等[12]提出的西南土石山區(qū)坡度因子的修正算法，計(jì)算法的基本形式如下：

(2)

(3)

2.1.4地表覆蓋與管理因子與水土保持措施因子廣西地區(qū)地形復(fù)雜，地域差異性顯著，以往基于USLE模型及其改進(jìn)模型開(kāi)展區(qū)域土壤侵蝕定量評(píng)價(jià)研究中，往往采用統(tǒng)一的C因子賦值方案，因而忽視了C因子在空間上的差異性。因此本文根據(jù)土壤侵蝕環(huán)境在時(shí)空上的差異性和相似性，將研究區(qū)劃分為5個(gè)土壤侵蝕預(yù)報(bào)區(qū)，通過(guò)空間抽樣調(diào)查技術(shù)“由點(diǎn)及面”計(jì)算C，P因子值，圖1為C，P因子計(jì)算流程圖。

圖1　C，P因子計(jì)算流程

地表覆蓋與管理因子(C)與土地利用類型和植被覆蓋度密切相關(guān)。參考楊子生[5]對(duì)農(nóng)作物經(jīng)營(yíng)C因子研究成果，根據(jù)廣西2 319個(gè)野外調(diào)查單元的成果資料，結(jié)合土地利用數(shù)據(jù)，按照不同土壤侵蝕預(yù)報(bào)分區(qū)先對(duì)耕標(biāo)的種植制度、種植方式進(jìn)行準(zhǔn)化賦值，利用landsat 8影像提取的NDVI值按照土壤侵蝕預(yù)報(bào)區(qū)計(jì)算植被覆蓋度，參考蔡崇法等[13]研究的植被覆蓋與管理因子定量關(guān)系，計(jì)算林地、園地、草地的C因子值，對(duì)于沒(méi)有土壤侵蝕的區(qū)域如水域、建設(shè)用地等統(tǒng)一賦值為0，地表完全裸露最容易受到侵蝕區(qū)域如裸地等，統(tǒng)一賦值為1，利用空間技術(shù)，制作全區(qū)C因子值分布圖(圖2)。

圖2　廣西西江流域C因子值分布

水土保持措施因子(P)反映了各類人工水土保持措施抵抗土壤侵蝕的能力。本文是在參考西南地區(qū)水土保持措施P值研究成果[5,14]，結(jié)合第一次全國(guó)水利普查水力侵蝕強(qiáng)度評(píng)價(jià)結(jié)果、第一次全國(guó)水利普野外調(diào)查單元成果、2014年底廣西土地利用變更調(diào)查成果，綜合考慮廣西水土保持工程措施和耕作措施實(shí)際情況，作統(tǒng)一校正后確定研究區(qū)P因子的賦值，利用空間技術(shù)，制作全區(qū)P因子值分布圖(圖3)。

2.1.5修正因子西江流域巖溶地貌分布廣泛，碳酸鹽巖類成土物質(zhì)少，且成土速度十分緩慢，土層淺薄且土體不連續(xù)性，抗外界干擾能力差，植被遭到破壞后自我修復(fù)能力較差，致使石漠化問(wèn)題較為嚴(yán)重。綜合目前土壤侵蝕研究成果[15-18]，喀斯特地區(qū)的土壤侵蝕，裸巖率與土壤相對(duì)流失量的關(guān)系呈現(xiàn)一個(gè)近似符合指數(shù)衰減關(guān)系(圖4)。因此，根據(jù)研究區(qū)的特點(diǎn)本文引入了喀斯特石漠化修正因子M。

圖3　廣西西江流域P因子值分布

石漠化信息提取的關(guān)鍵是多信息的綜合。本文利用地質(zhì)地貌圖，提取研究區(qū)喀斯特與非喀斯特區(qū)域，以巖石因子、地形因子、土地利用因子、植被因子、土壤因子等因子作為石漠化評(píng)價(jià)指標(biāo)，參考胡順光[19]、岳躍民[20]等對(duì)石漠化信息的提取，同時(shí)借鑒國(guó)家重點(diǎn)基礎(chǔ)研究計(jì)劃(973計(jì)劃)項(xiàng)目“典型山地水土要素時(shí)空耦合特征、效應(yīng)及其調(diào)控”研究成果，評(píng)估獲得了廣西地區(qū)的喀斯特石漠化數(shù)據(jù)。結(jié)合2011年底全國(guó)第二次石漠化監(jiān)測(cè)的結(jié)果，與研究區(qū)提取的喀斯特石漠化數(shù)據(jù)進(jìn)行輔助比對(duì)，由此計(jì)算得到了M因子(圖5)。

圖4　M因子與裸巖率關(guān)系

圖5　M因子分布

2.2　土壤侵蝕強(qiáng)度分級(jí)標(biāo)準(zhǔn)

本文將模型計(jì)算出的土壤侵蝕模數(shù)按照《巖溶地區(qū)水土流失綜合治理技術(shù)標(biāo)準(zhǔn)》(SL461—2009)和《南方紅壤丘陵區(qū)水土流失綜合治理技術(shù)標(biāo)準(zhǔn)》(SL657—2014)規(guī)定的容許土壤流失量，即巖溶地區(qū)容許土壤流失量為50 t/(km2·a)，南方紅壤丘陵區(qū)和非巖溶地區(qū)的西南土石山區(qū)容許土壤流失量為500 t/(km2·a)進(jìn)行分級(jí)(表1)。

表1　廣西西江流域土壤侵蝕強(qiáng)度分級(jí)標(biāo)準(zhǔn)　t/(km2·a)

2.3　土壤侵蝕精度驗(yàn)證

2.3.1基于水文輸沙驗(yàn)證水文驗(yàn)證是基于“輸沙量=土壤侵蝕量×泥沙輸移比”的公式，獲取了研究區(qū)8個(gè)水文站點(diǎn)的徑流量和輸沙量資料(2000—2015年)，在參考大量文獻(xiàn)[21-25]的基礎(chǔ)上，得出廣西的泥沙輸移比可近似取0.41，合理估計(jì)了廣西西江流域內(nèi)8個(gè)水文站點(diǎn)的泥沙輸移比，統(tǒng)計(jì)出了多年平均輸沙量與本文土壤侵蝕模型估算結(jié)果進(jìn)行對(duì)比。由表2可知，相對(duì)誤差最大為融水站16.38%，最小為寧明站3.08%，大于15%只有2個(gè)站點(diǎn)，其他6個(gè)站點(diǎn)均小于15%。由此可知推測(cè)輸沙量與實(shí)測(cè)輸沙量呈現(xiàn)高度一致，本文土壤侵蝕流失模型修正科學(xué)合理。

表2　基于水文資料驗(yàn)證結(jié)果對(duì)比分析

2.3.2與鄰省研究成果對(duì)比陳強(qiáng)等[26]在2005年研究滇東南西疇縣退耕還林的水土保持效應(yīng)中，觀測(cè)到坡耕地的土壤侵蝕模數(shù)為183 t/(km2·a)，本文利用模型計(jì)算出廣西與西疇縣鄰近的德?？h、靖西縣喀斯特旱地土壤侵蝕模數(shù)分別為：197.70，250.79 t/(km2·a)，模型計(jì)算結(jié)果與徑流小區(qū)觀測(cè)結(jié)果較為接近。

鄭悅?cè)A等[27]在使用近3年遙感影像對(duì)廣東青蓮小流域水土流失分析中，得出大于5°的平均土壤侵蝕模數(shù)724.86～1 146.72 t/(km2·a)，與本文研究區(qū)內(nèi)坡耕地平均土壤侵蝕模數(shù)912.04 t/(km2·a)較為接近。

3　土壤侵蝕流失特征分析

3.1　土壤侵蝕空間分布特征

本次土壤侵蝕估算研究中，在不考慮重力侵蝕、工程侵蝕等前提下，廣西西江流域土壤侵蝕主要集中分布在柳江水系西南西北部、左江流域中部及紅水河?xùn)|南部(圖6)。平均土壤侵蝕模數(shù)為135.22 t/(km2·a)，土壤侵蝕總量為2 737.72萬(wàn)t/a；從侵蝕面積上看(圖7)，輕度及以上的土壤侵蝕面積為22 802.75 km2，占流域土地總面積的11.26%；其中輕度侵蝕、中度侵蝕、強(qiáng)烈侵蝕、劇烈侵蝕、極劇烈侵蝕面積分別為15 375.45，5 291.16，1 293.14，620.50，222.50 km2，所占比例分別為7.59%，2.61%，0.64%，0.31%，0.11%。

3.2　不同土地利用土壤侵蝕特征

(1) 耕地土壤侵蝕特征分析。廣西西江流域耕地面積為41 650.36 km2，占土地總面積的20.57%，其中坡耕地面積為7 162.84 km2，占耕地面積的17.20%，梯田地的面積為16 118.01 km2，占耕地面積的38.70%。耕地平均土壤侵蝕模數(shù)為256.88 t/(km2·a)，土壤侵蝕量為1 069.92萬(wàn)t/a，侵蝕面積為8 888.14 km2，占耕地面積的21.34%；坡耕地平均土壤侵蝕模數(shù)為912.04 t/(km2·a)，侵蝕量為653.28萬(wàn)t/a,侵蝕面積為4 109.69 km2,占坡耕地面積的57.38%；梯田地平均土壤侵蝕模數(shù)為79.69 t/(km2·a)，土壤侵蝕量為128.45萬(wàn)t/a,侵蝕面積為2 303.88 km2,占梯田梯面積的14.29%。從侵蝕比例上看，坡耕地的土壤侵蝕比例最大。

圖6　廣西西江流域土壤侵蝕強(qiáng)度分布

圖7　各子流域土壤侵蝕面積

(2) 林地土壤侵蝕特征分析。森林在控制水土流失發(fā)揮著明顯的水土保持作用。廣西西江流域范圍內(nèi)，林地面積115 529.48 km2,占土地總面積的57.06%，林地平均土壤侵蝕模數(shù)為72.32 t/(km2·a)，遠(yuǎn)低于全流域的平均土壤侵蝕模數(shù)，土壤侵蝕量為835.56萬(wàn)t/a，侵蝕面積為5 921.15 km2,占林地面積的5.13%。

(3) 園地土壤侵蝕特征分析。園地面積為8 277.22 km2,僅占所有地類總面積的4.09%，園地平均土壤侵蝕模數(shù)為195.59 t/(km2·a)，土壤侵蝕量為161.90萬(wàn)t/a，土壤侵蝕面積為1 322.99 km2,占園地面積的15.98%，園地面積雖小但其土壤侵蝕比例較高。

(4) 草地土壤侵蝕特征分析。草地總面積為24 085.69 km2,占土地總面積的11.90%。草地平均土壤侵蝕模數(shù)為220.60 t/(km2·a)，土壤侵蝕量為531.32萬(wàn)t/a,土壤侵蝕面積為6 226.25 km2,占草地總面積的25.85%。

圖8　不同土地利用類型土壤侵蝕面積

3.3　不同地形地貌土壤侵蝕特征

廣西西江流域喀斯特地貌分布廣泛，地形復(fù)雜多樣。本文將坡度劃分為<3°，3°～5°，5°～8°，8°～15°，15°～25°，25°～35°，>35°共7個(gè)等級(jí)，地貌上劃分為平原、丘陵、臺(tái)地、小起伏山地、中起伏山地、大起伏山地，分別統(tǒng)計(jì)不同等級(jí)、類別下的土壤侵蝕分布狀況。

從表3可以看出，土壤侵蝕面積主要集中分布在8°～15°，15°～25°坡面上，從侵蝕比例上看，8°～15°坡面侵蝕面積比最大，為25.75%，其次為15°～25°坡面，侵蝕面積比為25.09%；從侵蝕量看，15°～25°坡面侵蝕量最大，為769.08萬(wàn)t/a，其次為8°～15°坡面，土壤侵蝕量為721.14萬(wàn)t/a。

從表4可以看出，大起伏山地平均土壤侵蝕模數(shù)最大，為203.41 t/(km2·a)，平原平均土壤侵蝕模數(shù)最小，為59.12 t/(km2·a)，小起伏山地侵蝕面積及侵蝕量最大，分別為11 251.93 km2，1 030.02萬(wàn)t/a。從侵蝕比例上看，丘陵侵蝕面積比最大，為15.82%，其次為起伏山地，侵蝕面積比為13.99%。由此可知小起伏山地、丘陵是水土流失的主要地貌類型。

表3　不同坡度等級(jí)土壤侵蝕情況

表4　不同地貌類型土壤侵蝕情況

4　結(jié) 論

(1) 在引入M因子修正的USLE模型中，模型計(jì)算得出的結(jié)果與基于水文資料驗(yàn)證的結(jié)果基本一致，與研究區(qū)鄰近區(qū)域相關(guān)研究成果較為接近，表明M因子值的修正科學(xué)合理。

(2) 在不考慮風(fēng)力侵蝕、重力侵蝕、工程侵蝕等前提下，廣西西江流域平均土壤侵蝕模數(shù)135.51 t/(km2·a)，侵蝕面積為22 802.75 km2，占流域總面積的11.26%，土壤侵蝕以輕度、中度侵蝕為主。

(3) 不同地類對(duì)土壤侵蝕程度不同。耕地是水土流失的主要地類，其次為草地、林地及園地，其中坡耕地是耕地土壤侵蝕的主要來(lái)源，同時(shí)也是西江流域水土流失的主要來(lái)源。

(4) 地形、地貌影響土壤侵蝕的空間分布。小起伏山地是西江流域水土流失的主要地貌類型；8°～25°的坡面為土壤侵蝕的易發(fā)區(qū)，水土流失最嚴(yán)重。

參考文獻(xiàn)：

[1]李國(guó)英.對(duì)我國(guó)水土保持工作的幾點(diǎn)思考[J].中國(guó)水土保持,1998(2):20-23.

[2]楊德生.珠江流域土壤侵蝕現(xiàn)狀及其特點(diǎn)[J].人民珠江,1990(4):38-39.

[3]Wischmeier W H, Smith D D. Predicting Rainfall-erosion Losses from Cropland East of the Rocky Mountain Guide for Selection of Practices for Soil and Water Conservation[Z]∥Agricultural Handbook 282. Washington, D.C: USDA-ARS,1965.

[4]劉寶元,畢小剛,符素華.北京土壤流失方程[M].北京:科學(xué)出版社,2010.

[5]楊子生.云南省金沙江流域土壤流失方程研究[J].山地學(xué)報(bào),2002,20(S1):1-9.

[6]江忠善,鄭粉莉,武敏.中國(guó)坡面水蝕預(yù)報(bào)模型[J].泥沙研究,2005(1):1-6.

[7]謝淑琴,劉予偉.西江流域水資源及其特性[J].人民長(zhǎng)江,2003(10):44-45.

[8]章文波,謝云,劉寶元.利用日雨量計(jì)算降雨侵蝕力的方法研究[J].地理科學(xué),2002,22(6):705-711.

[9]章文波,付金生.不同類型雨量資料估算降雨侵蝕力[J].資源科學(xué),2003,25(1):35-41.

[10]Sharply A N, Williams J R. EPIC-erosion/Productivity impact calculator: Model documentation[R]. Washington: United States Department of Agriculture, 1990.

[11]張科利,彭文英,樣紅麗.中國(guó)土壤可蝕性值及其估算[J].土壤學(xué)報(bào),2007,44(1):8-13.

[12]劉斌濤,宋春風(fēng),史展,等.西南土石山區(qū)土壤流失方程坡度因子修正算法研究[J].中國(guó)水土保持,2015(8):49-51,77.

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

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

[15]王小燕,李朝霞,徐勤學(xué),等.礫石覆蓋對(duì)土壤水蝕過(guò)程影響的研究進(jìn)展[J].中國(guó)水土保持科學(xué),2011,9(1):115-120.

[16]Martinez M, Lane L J, Fogel M M. Experimental Investigation of Soil Detachment by Raindrop Impacts[Z]. Agricultural Reviews and Manuals ARMW, 1979:152-155.

[17]符素華,路炳軍,葉芝菡.地表礫石對(duì)降雨徑流及土壤侵蝕的影響[J].水土保持學(xué)報(bào),2010,24(2):15-18.

[18]王濟(jì),蔡雄飛,雷麗,等.不同裸巖率下我國(guó)西南喀斯特山區(qū)土壤侵蝕的室內(nèi)模擬[J].中國(guó)巖溶,2010,29(1):1-5.

[19]胡順光,張?jiān)鱿?夏奎菊.遙感石漠化信息的提取[J].地球信息科學(xué)學(xué)報(bào),2010(6):870-879.

[20]岳躍民,王克林,張兵,等.喀斯特石漠化信息遙感提取的不確定性[J].地球科學(xué)進(jìn)展,2011(3):266-274.

[21]文安邦,張信寶,王玉寬,等.云貴高原區(qū)龍川江上游泥沙輸移比研究[J].水土保持學(xué)報(bào),2003,17(4):140-141.

[22]許烔心.三峽水庫(kù)修建前長(zhǎng)江宜昌—武漢段泥沙輸移比及其影響因子[J].山地學(xué)報(bào),2008,26(1):15-25.

[23]景可,焦菊英,李林育,等.中國(guó)南方紅壤丘陵區(qū)土壤侵蝕量與泥沙輸移比研究[J].水利學(xué)報(bào),2010,41(12):1476-1482.

[24]馬琨,王兆賽,陳欣.紅壤坡面產(chǎn)流產(chǎn)沙與養(yǎng)分流失特征研究[J].寧夏農(nóng)學(xué)院學(xué)報(bào),2003,24(2):3-7.

[25]花利忠,賀秀斌,顏昌宙,等.基于AnnAGNPS模型的大寧河流域泥沙輸移比評(píng)價(jià)[J].水土保持通報(bào),2009,29(3):154-158.

[26]陳強(qiáng),常恩福,畢波,等.滇東南巖溶地區(qū)3種退耕還林模式的水土保持效應(yīng)研究[J].水土保持學(xué)報(bào),2006,20(5):1-4.

[27]鄭悅?cè)A,張曉遠(yuǎn),劉協(xié)亭.基于GIS的粵北青蓮水流域水土流失成因分析[J].廣東水利電,2016(5):24-28.