廖元群，廖富強(qiáng)，萬(wàn)小星，齊述華

(1.江西師范大學(xué) 鄱陽(yáng)湖濕地與流域研究教育部重點(diǎn)實(shí)驗(yàn)室/地理與環(huán)境學(xué)院，江西 南昌 330022；2.江西省水政監(jiān)察總隊(duì)，江西 南昌 330009)

水土流失是在外營(yíng)力(水力、風(fēng)力、重力、凍融張縮力等)作用下，土壤物質(zhì)剝離、遷移或沉積的過(guò)程。水土流失的發(fā)生主要是地形、土壤、降雨、植被及土地利用方式等多因子共同作用的結(jié)果[1-2]。目前，常用于區(qū)域水土流失評(píng)價(jià)的經(jīng)驗(yàn)?zāi)Ｐ陀蠻SLE、RUSLE和CSLE。2018年水利部印發(fā)了《區(qū)域水土流失動(dòng)態(tài)監(jiān)測(cè)技術(shù)規(guī)定(試行)》，要求采用CSLE模型開(kāi)展全國(guó)全覆蓋的水土流失動(dòng)態(tài)監(jiān)測(cè)。地形因子是CSLE模型中的重要參數(shù)，包括坡長(zhǎng)因子(L)和坡度因子(S)，其中：L是指其他條件一致的情況下，某種坡長(zhǎng)坡面土壤侵蝕量與22.13 m坡長(zhǎng)坡面土壤侵蝕量的比值[3]；S是指降雨、土壤、坡長(zhǎng)、地表狀況等條件一致的情況下，某坡度下坡面土壤侵蝕量與坡度為5.14°時(shí)坡面土壤侵蝕量的比值[3]。

研究人員對(duì)地形因子的計(jì)算方法及其尺度依賴(lài)性進(jìn)行了研究。梁曉珍等[4]分別利用分段坡長(zhǎng)法和匯流面積法計(jì)算L、S因子值，并分析了不同方法計(jì)算的L因子值及其對(duì)水土流失評(píng)價(jià)的影響；矯欣航[5]利用同一地形資料空間內(nèi)插重采樣成不同空間分辨率DEM，分析不同分辨率的DEM計(jì)算的L、S因子值對(duì)水土流失評(píng)價(jià)的影響；羊秀娟等[6]利用GDEMDEM和SRTM DEM重采樣為不同分辨率，分析地形因子的分辨率對(duì)水土流失評(píng)價(jià)的影響；李新艷[7]分析了不同地貌區(qū)的宏觀(guān)和微觀(guān)地形因子對(duì)區(qū)域水土流失評(píng)價(jià)的影響；LU Shaojuan et al.[8]比較了5種地形數(shù)據(jù)計(jì)算的L、S因子值。上述成果主要集中于研究不同地形因子提取方法、不同分辨率和不同地貌區(qū)的地形因子與土壤侵蝕的關(guān)系，而對(duì)山地地貌區(qū)不同尺度地形因子及其對(duì)區(qū)域水土流失評(píng)價(jià)結(jié)果的影響研究較少?；谏鲜鲅芯勘尘埃狙芯窟x擇山地特征明顯的崇義縣，分別利用1∶1萬(wàn)、1∶5萬(wàn)地形圖和ASTER GDEM數(shù)據(jù)，比較山地條件下L、S因子的尺度依賴(lài)性及其對(duì)區(qū)域水土流失評(píng)價(jià)結(jié)果的影響，希望能為山區(qū)水土流失動(dòng)態(tài)監(jiān)測(cè)提供參考。

1 研究區(qū)概況

崇義縣地處江西省西南邊陲，隸屬贛州市，位于東經(jīng)113°55′～114°38′、北緯25°24′～25°55′之間，東西長(zhǎng)約73 km，南北寬約59 km，國(guó)土總面積2 197 km2，東與南康區(qū)接壤，南與大余縣和廣東省仁化縣相交，西與湖南省汝城縣、桂東縣毗鄰，北與上猶縣交界。全縣地勢(shì)由西南向東北方向傾斜，按地貌特征分為中山、低山、高丘陵、河谷階地4種類(lèi)型，見(jiàn)圖1。

圖1 研究區(qū)概況

2 研究數(shù)據(jù)與方法

2.1 土地利用數(shù)據(jù)

利用高分一號(hào)衛(wèi)星2 m分辨率光學(xué)遙感影像，依據(jù)《土地利用現(xiàn)狀分類(lèi)》(GB/T 21010—2017)制定的分類(lèi)體系，以人工目視解譯為主，生成崇義縣土地利用空間分布數(shù)據(jù)。

2.2 地形數(shù)據(jù)及L、S因子計(jì)算

在利用崇義縣1∶1萬(wàn)、1∶5萬(wàn)等高線(xiàn)高程數(shù)據(jù)和ASTER GDEM30 m空間分辨率的柵格高程數(shù)據(jù)，通過(guò)反距離加權(quán)算法重采樣，插值生成10 m空間分辨率DEM數(shù)據(jù)的基礎(chǔ)上，利用符素華等[3]研發(fā)的地形因子軟件計(jì)算L、S因子值。該算法充分考慮了坡耕地的特點(diǎn)，結(jié)合MCCOOL提出的公式[9]和LIU Baoyuan提出的公式[10]計(jì)算S因子，計(jì)算公式為

(1)

式中：θ為坡度，(°)；S為坡度因子，取值范圍為[0, 9.99]。

采用FOSTER et al.[11]提出的分段坡長(zhǎng)法計(jì)算柵格中每個(gè)像元的坡長(zhǎng)因子值，計(jì)算公式[12]為

(2)

其中

(3)

上兩式中：Li為第i個(gè)柵格的坡長(zhǎng)因子；λi、λi-1分別為第i個(gè)和第i-1個(gè)坡段的坡長(zhǎng)，m；m為坡長(zhǎng)指數(shù)。

2.3 CSLE模型

采用劉寶元等[12]提出的中國(guó)土壤流失方程CSLE模型，計(jì)算公式為

A=R·K·L·S·B·E·T

(4)

式中：A為土壤侵蝕模數(shù)，t/(hm2·a)；R為降雨侵蝕力因子，MJ·mm/(hm2·h·a)；K為土壤可蝕性因子，t·hm2·h/(hm2·MJ·mm)；B為植被覆蓋與生物措施因子，無(wú)量綱；E為工程措施因子，無(wú)量綱；T為耕作措施因子，無(wú)量綱。

其中：R、K因子數(shù)據(jù)來(lái)自水利部信息中心；采用Landsat-8衛(wèi)星獲取的光學(xué)影像計(jì)算的NDVI和MODIS 16天合成23時(shí)相NDVI數(shù)據(jù)產(chǎn)品，通過(guò)融合法計(jì)算24個(gè)半月尺度的NDVI，再按照像元二分模型計(jì)算植被覆蓋度，分別計(jì)算園地、林地和草地的B因子值；依據(jù)高分衛(wèi)星遙感影像提取的水土保持工程措施確定E因子值，研究區(qū)提取的工程措施主要有土坎水平梯田、水平階(反坡梯田)和地埂，分別賦值為0.084、0.151和0.347；依據(jù)《中國(guó)耕作制度70年》，崇義縣屬于南嶺丘陵山地水田旱地二熟三熟區(qū)，耕地的耕作措施因子賦值0.338。

3 結(jié)果與分析

3.1 不同地形數(shù)據(jù)的地形特征差異

坡度是地形特征的要素，坡度的大小直接影響地表的物質(zhì)流動(dòng)與能量轉(zhuǎn)換的規(guī)模與強(qiáng)度[6]。按照《土壤侵蝕分類(lèi)分級(jí)標(biāo)準(zhǔn)》(SL 190—2007)中的坡度分級(jí)標(biāo)準(zhǔn)，將坡度分為6級(jí)(0°～5°、5°～8°、8°～15°、15°～25°、25°～35°、≥35°)，以1∶1萬(wàn)、1∶5萬(wàn)和ASTER地形數(shù)據(jù)進(jìn)行坡度分析，結(jié)果見(jiàn)圖2。從圖2可以看出：作為山區(qū)縣，崇義縣8°～35°坡度帶分布面積較大；在0°～15°緩坡帶，1∶1萬(wàn)地形數(shù)據(jù)土地面積大于A(yíng)STER和1∶5萬(wàn)；在≥25°陡坡度帶，1∶5萬(wàn)地形數(shù)據(jù)土地面積明顯大于A(yíng)STER和1∶1萬(wàn)；1∶1萬(wàn)和1∶5萬(wàn)高程數(shù)據(jù)中，比例尺越大，計(jì)算的平均坡度越小，陡坡區(qū)分布面積占比越小。

圖2 不同地形數(shù)據(jù)坡度分級(jí)面積分布

依據(jù)崇義縣地形地貌特點(diǎn)，將高程分為4個(gè)等級(jí)(表1)：低丘及河谷階地(<300 m)、高丘(300～500 m)、低山(500～1 000 m)和中山(1 000～3 500 m)。分別統(tǒng)計(jì)各級(jí)高程面積，結(jié)果表明3種地形數(shù)據(jù)在各高程帶的分布面積相近。

表1 不同地形數(shù)據(jù)高程面積分級(jí)

3.2 不同地形數(shù)據(jù)計(jì)算的坡長(zhǎng)因子差異

根據(jù)崇義縣坡長(zhǎng)因子空間分布(圖3)，統(tǒng)計(jì)分析不同地形數(shù)據(jù)計(jì)算的L因子頻率分布(圖4)。坡長(zhǎng)因子分級(jí)面積占比計(jì)算結(jié)果見(jiàn)表2。結(jié)果表明：基于1∶1萬(wàn)、1∶5萬(wàn)和ASTER地形數(shù)據(jù)計(jì)算的L因子平均值分別為1.78、1.69、1.63，離差系數(shù)分別為0.48、0.54、0.62；在L因子1.0～1.5和≥2.5分級(jí)范圍內(nèi)的面積占比均表現(xiàn)為1∶5萬(wàn)>ASTER>1∶1萬(wàn)，且在0.5～1.0的分級(jí)范圍內(nèi)1∶5萬(wàn)、ASTER地形數(shù)據(jù)面積占比均為0；在L因子1.5～2.0分級(jí)范圍內(nèi)的面積占比1∶1萬(wàn)>1∶5萬(wàn)>ASTER。

表2 坡長(zhǎng)因子分級(jí)面積占比

圖3 崇義縣不同地形數(shù)據(jù)計(jì)算的L因子空間分布

圖4 崇義縣不同地形數(shù)據(jù)計(jì)算的L因子頻率分布

3.3 不同地形數(shù)據(jù)計(jì)算的坡度因子差異

根據(jù)崇義縣坡度因子空間分布(圖5)，統(tǒng)計(jì)分析不同地形數(shù)據(jù)計(jì)算的S因子頻率分布(圖6)。坡度因子分級(jí)面積占比計(jì)算結(jié)果見(jiàn)表3。結(jié)果表明：基于1∶1萬(wàn)、1∶5萬(wàn)和ASTER地形數(shù)據(jù)計(jì)算的S因子平均值分別為2.66、3.93、3.04，離差系數(shù)分別為0.56、0.41、0.52；1∶1萬(wàn)地形數(shù)據(jù)S低值區(qū)的面積占比高，1∶5萬(wàn)地形數(shù)據(jù)S高值區(qū)的面積占比高；1∶1萬(wàn)地形數(shù)據(jù)中85.66%的區(qū)域S<7.5，1∶5萬(wàn)地形數(shù)據(jù)中77.50%的區(qū)域S>4.5，ASTER地形數(shù)據(jù)中S值在各分級(jí)范圍內(nèi)分布較為均勻。

圖5 崇義縣不同地形數(shù)據(jù)計(jì)算的S因子空間分布

圖6 崇義縣不同地形數(shù)據(jù)計(jì)算的S因子頻率分布

表3 坡度因子分級(jí)面積比例

3.4 基于不同地形數(shù)據(jù)計(jì)算的水土流失面積差異

根據(jù)《區(qū)域水土流失動(dòng)態(tài)監(jiān)測(cè)技術(shù)規(guī)定(試行)》，各土壤侵蝕因子重采樣為10 m空間分辨率，分別利用不同來(lái)源DEM計(jì)算的地形因子，以CSLE模型計(jì)算土壤侵蝕模數(shù)。根據(jù)《土壤侵蝕分類(lèi)分級(jí)標(biāo)準(zhǔn)》(SL 190—2007)，將土壤侵蝕強(qiáng)度分為微度、輕度、中度、強(qiáng)烈、極強(qiáng)烈和劇烈侵蝕等6個(gè)等級(jí)，其中輕度及以上侵蝕區(qū)為水土流失區(qū)。根據(jù)崇義縣土壤侵蝕強(qiáng)度空間分布(圖7)，統(tǒng)計(jì)不同侵蝕強(qiáng)度面積(表4)，結(jié)果表明：3種地形數(shù)據(jù)評(píng)價(jià)的崇義縣水土流失面積均以輕度侵蝕為主，但地形數(shù)據(jù)的差異對(duì)水土流失強(qiáng)度評(píng)價(jià)的結(jié)果影響顯著，水土流失面積1∶5萬(wàn)>ASTER >1∶1萬(wàn)。

表4 土壤侵蝕面積 km2

圖7 崇義縣不同地形數(shù)據(jù)計(jì)算的土壤侵蝕強(qiáng)度空間分布

利用坡度與土壤侵蝕強(qiáng)度進(jìn)行疊加分析，結(jié)果見(jiàn)圖8。從圖8可以看出，崇義縣以輕度侵蝕為主，主要分布在15°～35°坡度帶；1∶1萬(wàn)和ASTER地形數(shù)據(jù)計(jì)算的土壤侵蝕強(qiáng)度在各坡度帶分布的面積比例較相近；從各坡度帶的土壤侵蝕分布情況看，1∶5萬(wàn)地形數(shù)據(jù)計(jì)算的中度以上侵蝕比1∶1萬(wàn)和ASTER地形數(shù)據(jù)均高一個(gè)坡度帶等級(jí)。結(jié)果說(shuō)明，在崇義縣坡度是影響水土流失強(qiáng)度的重要因素。

圖8 不同坡度帶水土流失面積占比

利用高程與土壤侵蝕強(qiáng)度進(jìn)行疊加分析，結(jié)果見(jiàn)圖9。從圖9可以看出，崇義縣土壤侵蝕主要分布在300～1 000 m高程帶；1∶1萬(wàn)和ASTER地形數(shù)據(jù)計(jì)算的極強(qiáng)烈以上侵蝕在各高程帶分布的面積比例較相近；1∶5萬(wàn)和ASTER地形數(shù)據(jù)計(jì)算的中度和強(qiáng)烈侵蝕在各高程帶分布的面積比例較相近。

圖9 不同高程帶水土流失面積占比

4 結(jié) 論

以崇義縣為例，分別采用1∶1萬(wàn)、1∶5萬(wàn)和ASTER地形數(shù)據(jù)計(jì)算L、S因子值和土壤侵蝕模數(shù)，探討山區(qū)條件下地形因子的尺度依賴(lài)性及其對(duì)水土流失評(píng)價(jià)的影響，得到以下結(jié)論：

(1)1∶1萬(wàn)地形數(shù)據(jù)計(jì)算的L因子主要分布在2.0～2.5分級(jí)范圍內(nèi)，1∶5萬(wàn)地形數(shù)據(jù)計(jì)算的L因子主要分布在1.0～1.5、2.0～2.5分級(jí)范圍內(nèi)，ASTER地形數(shù)據(jù)計(jì)算的L因子主要分布在<0.5、2.0～2.5分級(jí)范圍內(nèi)。

(2)1∶1萬(wàn)地形數(shù)據(jù)計(jì)算的S因子主要分布在中低值區(qū)，1∶5萬(wàn)地形數(shù)據(jù)計(jì)算的S因子主要分布在中高值區(qū)，ASTER地形數(shù)據(jù)在各分級(jí)分布較為均勻。

(3)不同地形因子的尺度效應(yīng)對(duì)水土流失評(píng)價(jià)結(jié)果影響較大，不同地形因子計(jì)算的水土流失面積1∶5萬(wàn)>ASTER>1∶1萬(wàn)。從坡度、高程與土壤侵蝕的關(guān)系看，在不同坡度帶和高程帶中，不同地形因子數(shù)據(jù)計(jì)算的土壤侵蝕空間分布具有一致性，各強(qiáng)度等級(jí)間的土壤侵蝕空間分布存在差異。