干旱半干旱區坡面覆被格局的水土流失效應研究進展

高光耀，傅伯杰，呂一河，劉 宇，3，王 帥，周 繼

(1.中國科學院生態環境研究中心城市與區域生態國家重點實驗室，北京 100085;2.黃土高原土壤侵蝕與旱地農業國家重點實驗室，中國科學院水利部水土保持研究所，陜西楊凌 712100;3.中國科學院遙感應用研究所農業與生態遙感研究室，北京 100101)

土地利用/土地覆被結構與生態過程研究是對綜合自然地理學的深化和發展，是地表過程研究的前沿［1－2］。水土流失是世界性的生態環境問題，我國是世界上水土流失最嚴重的國家之一，而黃土高原則是我國土壤侵蝕最嚴重的地區。水土流失使土地退化，生產力下降，生態環境惡化。由此產生的大量泥沙淤塞江河湖泊，加劇洪水災害，影響和制約區域經濟和可持續發展。除自然因素外，不合理的土地利用是造成土壤侵蝕的重要驅動因子。土地利用與水土流失的關系一直是國內外普遍關注的重要科學問題［3－5］。

植被恢復與建設是減緩或防治水土流失的重要措施。例如在我國黃土高原實施的小流域綜合治理和退耕還林(草)工程，旨在通過黃土高原地區的植被恢復與生態環境建設達到控制水土流失的目標。要想有效控制水土流失，一方面要選擇合理的植物物種及其搭配，另一方面也要合理地設計植被空間分布格局。特別是在干旱半干旱地區，水分制約加上人類活動等眾多因素影響，植被恢復不可能達到理想的完全覆蓋狀況，而會在空間上呈現斑點(塊、簇)狀或條帶狀分布等離散特征，在坡面上會形成裸地與植被斑塊的鑲嵌、不同植被類型斑塊的組合與裸地的鑲嵌以及植被的條帶分布等典型覆被格局［5－7］。例如，在黃土高原溝坡等坡度較大、生態系統退化較嚴重的地段，植被以隨機性斑塊化格局為主;而在坡度較小、侵蝕和生態系統退化程度較低的坡面，植被格局以人為干擾和自然因素綜合作用形成的帶狀特征為主［4］。在地中海干旱半干旱區，人為干預和土地棄耕等交替活動的驅動使植被群落經過不同階段的退化、更新，最終形成鑲嵌式的空間結構［7－8］。顯然，具有不同覆被格局的坡面，其水土流失特征會有明顯差異。坡面是各種地表過程發生發展的重要地理單元，也是流域/區域的基本構成單元。坡面尺度水土流失過程的研究對探討大尺度上土地利用格局演變的水土流失效應具有重要作用，而且坡面的土壤侵蝕防治是流域和區域尺度水土保持工作的關鍵，具有重要的理論和實踐意義。因此，干旱半干旱區坡面尺度上土地覆被格局的水土流失效應研究一直是國內外的研究熱點。

基于文獻調研，本文從植被斑塊對水土流失的影響、坡面覆被格局對水土流失的影響和耦合覆被格局與水土流失的手段與方法3個方面對干旱半干旱區坡面覆被格局與水土流失關系的相關研究進行綜述與探討，歸納需要進一步研究的問題，旨在從過程和機制上更加深刻認識覆被格局對水土流失過程的影響機理，為今后開展土地覆被結構與水土流失相關研究提供一定參考依據，并為坡面植被恢復格局設計、有效控制水土流失和改善干旱半干旱區的水資源狀況等提供理論依據。

1 植被斑塊對水土流失的影響

斑塊尺度植被對水土流失有著直接和間接的影響。首先，植被以及枯枝落葉層能減輕雨滴的直接打擊以及對土壤的剝離，有效的攔截地表徑流，減緩其流速，很大程度上減弱或消除降雨和徑流侵蝕能量，減少侵蝕發生的機會［5，9］。植被對徑流泥沙的影響關鍵在于改變了降雨侵蝕動能。植被生長引起植物高度增加、根系伸長、葉面積增加、覆蓋度增大、植被結構發生變化，植被對降雨的阻截再分配作用因此發生變化，改變了降雨對地表的侵蝕擊濺力［10］。其次，植被及枯落物通過改善表土結構、物理化學和水文性質以及微地形，進而影響水土流失，如降低土壤可蝕性，增加入滲能力等［5，10］。植被影響水土流失的主要因素有植被類型、植被覆蓋度、植被枯枝落葉層、植被高度和植被根系等［5，10－11］。

不同的植被類型，有不同的葉片形態和層次結構，各個層次的形態也有顯著差異，進而對水土流失的影響會有不同。大葉片植被更容易將雨滴匯聚成更大的水滴，增強雨滴濺蝕。具有多層結構的植被群落比單層植被更能保護土壤，減輕侵蝕強度［12］，而林草復合的植被覆蓋結構比純林或純草覆蓋具有更強的水土保持功能［13］。植被生長形態也具有一定的作用，而植被演替能引起群落結構發生改變，從而改變降雨對地表的侵蝕能力。國內外學者對不同植被類型、不同植被的層次及形態結構對土壤侵蝕的控制作用進行了大量研究［5，9－11］。

植被蓋度是影響水土流失的關鍵因素［14］。植被覆蓋度的增加會加強對降雨的攔截，降低雨滴動能，進而減弱降雨侵蝕力。實驗結果表明，徑流量和產沙量隨著覆蓋度增加而呈線性［9，15］或指數［9，16－17］減小。很多研究也探討了有效植被蓋度的問題，認為只有達到一定蓋度之后植被才能起到有效減輕土壤侵蝕的作用［15，18－19］。葉面積指數也是評價植被水土保持效益的重要指標。與植被覆蓋度相比，植被葉面積指數反映了冠層葉片空間結構，在表征土壤侵蝕、評價植被水土保持效果上更為穩定和可靠［13］。

植被高度的變化也會顯著影響土壤侵蝕。當植被覆蓋度不變而高度，尤其是冠層距地面高度變化時，由于冠層下雨滴直徑和動能變化，濺蝕量會發生變化。研究表明只有在一定高度下，植被覆蓋才能有效地降低雨滴動能，植被過高其冠層匯集的雨滴能量可能更大，對地表的濺蝕更強［5，10］。Morgan［20］的研究表明，當草本高度大于0.3 m時，植被保護土壤不被雨滴濺蝕的能力降低，當大于1 m時，雨滴的濺蝕會超過自然雨滴。蔡強國和陳浩［21］的實驗資料表明當植被高度為3—4 m時植被已無消除雨滴動能的功能。

植被枯枝落葉層是控制土壤侵蝕的重要因素?？萋湮镉行p弱雨滴對土壤顆粒的分散動力，具有較大的水分截持能力，影響土壤入滲和蒸發，同時會改變地表糙率等地形結構，顯著降低水流挾沙能力，減少徑流匯集及泥沙搬運?？萋湮锏慕M成結構、持水能力、厚度、蓄積量等因素是影響土壤侵蝕的重要指標［22］。雖然具有相似的覆蓋度，不同植被由于枯落物特性不同，土壤侵蝕程度會不同。例如，以叢生禾草為主的混合草地的枯落物厚度和粒度要大于草坪草地，使得前者的土壤侵蝕程度要低［23］。

除了植被的地上部分，植物根系對土壤侵蝕也有重要影響。Gyssels等［24］與Gyssels和Poesen［25］在比利時黃土帶對谷類作物和草類植被的土壤侵蝕效應研究表明，根系密度和秧苗密度的增加均能導致侵蝕速率呈指數形式降低，特別是在植物生長初期，根系對土壤侵蝕具有重要的控制作用。李勇［26］在黃土高原的研究表明，植物根系能提高土壤抗沖性和入滲能力?？偟膩砜?，植被從地上部分的冠層到地下部分的根系，即植被的垂直結構，都對水土流失有著直接或間接的作用。

2 覆被格局對水土流失的影響

植被是土壤侵蝕的重要影響因素，而覆被格局即各種覆被類型包括植被、枯落物和裸地等的空間分布和數量結構對坡面水土流失具有直接的控制作用。覆被空間分布格局影響地表徑流的匯集和攜沙能力，改變了水土流失過程的連續性［5，27］。在以超滲產流為主的干旱半干旱地區，裸地往往由于土壤結皮的存在而降低入滲率，增加地表徑流，而有植被覆蓋的土壤則因土壤特性的改善，入滲率增大，成為降水吸收區，從而導致徑流泥沙源匯區的產生［10－11］。植被空間分布格局改變了徑流和泥沙運移路徑的連通性，對水土流失作用顯著［28－31］。

Boer和Puigdefábregas［32］的模擬研究發現裸露區域和植被的空間組織對暴雨徑流和侵蝕的預測有顯著影響。Boix－Fayos等［33］在西班牙東南部半干旱區的長期監測研究表明，相同規格的監測樣地因覆被格局不同，總流失量可相差9倍之多。Bautista等［30］的研究表明，在半干旱環境，植被空間結構屬性和植物功能多樣性是比表層土壤屬性更好的徑流和侵蝕解釋變量。盡管植被覆蓋度和生物量是水文模擬中常用的植被屬性，而這些變量不足以預測半干旱景觀的徑流和侵蝕量［32］。鑒于覆被格局對產流產沙的重要影響，在水土流失研究中，格局對過程的作用日益得到重視。針對覆被格局與水土流失過程的關系研究逐漸增多和不斷深入，景觀生態學的格局分析方法逐漸應用到其中。

國際上的很多相關研究都是在定位觀測基礎上，分析植被斑塊以及裸地－植被斑塊鑲嵌格局的水土流失效應［28，33］。Reid等［34］在美國新墨西哥州北部研究了林下植被斑塊、林間植被斑塊和林間裸地斑塊等3種覆被斑塊類型的水土流失效應。結果表明，裸露斑塊的產流和產沙量要大于林間植被斑塊，產流產沙最小的是林下植被斑塊，而且裸露斑塊是水土流失的源，其他兩種斑塊類型起到了匯的作用［34］。以色列荒漠化灌叢區的實驗研究也表明，植被斑塊越少，越易造成水土流失，增加發揮土壤和水分“匯”功能的植被或其他斑塊并對它們合理布局，有助于生態系統的保護和恢復［35］。Puigdefábregas［28］分析了干旱區斑塊和立地尺度上植被的非均勻分布特征對徑流和泥沙的影響。Boer和Puigdefábregas［32］的模擬研究發現坡面空間結構化的植被格局(植被斑塊和裸露斑塊鑲嵌分布)比均一化植被格局下的產流產沙量大。Ludwig等［36］在澳大利亞東北部開展的微型小區和坡面尺度觀測研究表明，對于具有相近草本蓋度、坡度和土壤類型的兩個對比坡面，植被空間格局粗?；矣休^大裸地斑塊的坡面，其侵蝕速率是植被格局細粒化坡面的40多倍，后者可以直接由微型小區的觀測結果線性尺度上推，而前者的實測值是線性遞推結果的2.5倍。這主要是由于植被格局粗?；旅娴母脖话邏K結構未能有效降低徑流速度和留滯泥沙，說明坡面植被斑塊的空間離散分布使得其生態水文過程存在跨尺度的非線性特征。Bautista等［30］在西班牙地中海半干旱景觀坡面尺度上對植被空間格局和功能多樣性與坡面水文功能的關系進行監測分析，發現徑流量和斑塊密度呈負相關，而且產流產沙量均隨植被空間格局的粗粒化而增加。Bartley等［17，29］在澳大利亞北部半干旱區牧場小流域內的坡面上長達6a的監測試驗表明，當坡面上分布有大量裸地小斑塊時，其徑流量比具有較低植被覆蓋度但裸露斑塊較少的坡面產流量高6—9倍，而土壤流失量差別可達60倍，特別是當低覆蓋度植被斑塊靠近坡底和溝道時，坡面產流產沙量顯著增加。該研究充分表明裸地斑塊及其分布位置在坡面徑流侵蝕中的重要作用，而在坡底位置布設中、高覆蓋度植被斑塊可以有效阻截泥沙，減少對河網的泥沙輸入。

坡面上植被的條帶分布特征也引起了國際上很多學者的注意。植被條帶類型、寬度、覆蓋度和位置是影響其控制水土流失的重要因素［37］。Ludwig等［38］在澳大利亞東部半干旱區坡面上的模擬結果表明，去除灌叢斑塊后，其徑流截持率比條帶狀格局降低約25%，而條紋狀或線狀格局的灌木地比點狀鑲嵌格局的截流能力提高約8%。Valentin等［39］通過總結相關研究結果表明，坡面植被－裸地帶狀格局形成產流－聚流系統，有利于截獲有限的水資源，構成限制土壤侵蝕的天然屏障，而且這種帶狀格局下的生物量往往高于點狀或單一點狀或散亂分布格局。Rey［40］研究了坡底草本和灌叢等低矮植被對上坡徑流泥沙的攔截作用，指出植被屏障蓋度達到20%就可以有效攔截上坡的侵蝕泥沙。Martníez Raya等［41］在西班牙東南部坡地的山杏林間種植不同類型的植被條帶，4a的徑流小區定位觀測結果表明植被條帶能有效減少徑流和產沙，而且不同植被條帶類型控制水土流失的效果差異明顯。

與國際上的研究相比，國內關于坡面覆被格局與水土流失過程關系的實驗研究顯得比較薄弱。蘇敏等［42］在黃土丘陵溝壑山坡地的小區內進行草灌帶狀間作試驗和草糧帶狀輪作試驗，評價了不同植被配置方式的水土保持效應，結果表明草灌帶間作配置方式的水保性能要強于草糧帶間作。游珍等［43］在寧夏固原的自然荒草坡面小區內，從坡頂或坡底開始剪除一定長度的地表植被，形成不同的坡面植被分布，研究其對降雨侵蝕的影響，結果認為:在相同面積條件下，位于坡底的植被比位于坡頂的植被保水作用高2.4倍，保土作用高2.8倍。李勉等［44］采用室內放水沖刷試驗方法，研究了坡面不同草被覆蓋度及空間配置下坡溝系統的侵蝕產沙變化過程及特征，結果表明坡面草被空間配置方式的產沙量大小順序為坡上部＞坡中部＞坡下部。沈中原［45］在寧夏固原的荒草坡地上，通過人工處理的方式，研究了草地聚集結構(坡頂、坡中、坡底)、帶狀結構、隨機結構和網絡結構等多種格局對水土流失效應的影響，結果表明五種植被格局坡面產沙量大小順序為坡頂聚焦＞坡底聚焦＞坡中聚焦＞帶狀格局＞隨機格局，產沙量大小順序為帶狀格局＞坡中聚焦＞坡頂聚焦＞隨機格局＞坡底聚焦。徐海燕等［46］基于安塞試驗站的坡面小區實驗，研究坡耕地與草地不同配置方式的侵蝕產沙特征，所考慮的土地利用組合方式為上坡是谷子，下坡是不同類型的草本或撂荒地，結果表明谷子－撂荒地配置方式的減沙效應最好。朱冰冰等［18］以植被覆蓋完整、均一分布的自然荒草坡為原狀坡面，從坡頂向下剪除植被，形成不同覆蓋度的小區，通過人工模擬降雨試驗，研究草本植被覆蓋度對坡面降雨徑流侵蝕的影響，并確定臨界植被覆蓋度為60%—80%。在整個坡面上，傅伯杰等在黃土丘陵區研究了不同土地利用結構對土壤水分分布的影響［4］，并通過137Cs示蹤量化了不同土地利用組合的土壤侵蝕強度差異，結果表明坡中部位林地和草地結構能夠抑制土壤侵蝕，坡面土壤侵蝕量受土地利用類型分布和坡位共同影響［47］?？梢钥闯?，國內關于覆被格局對水土流失作用的實測研究大都集中在不同農作物的種植方式以及不同植被的人工規則組合對土壤侵蝕的影響，且多以定性或半定量分析為主，缺少對自然坡面覆被格局－水土流失過程關系定量化的機理研究。

分析以上國內外相關進展可以看出，覆被格局影響水土流失的關鍵是坡面徑流泥沙源匯區的連通性和空間分布［7，28，33］。植被斑塊、枯落物、碎石和地表坑洼作為地表物質流的截留、阻礙功能單元，阻滯產流產沙，是徑流和泥沙的主要匯區域，而植被斑塊間的裸露地表則是徑流和泥沙的源區域［31，48］。源匯區的連通性與降雨－徑流－侵蝕過程密切相關，該連通性以地表土壤顆粒的剝離和沿地表的移動為目標過程，以地表水文連接性為基礎，這與以生物過程為對象的景觀連接性存在明顯區別［31］。高強度的降雨、陡峭的地形以及粗糙度低的裸露土壤等有利于增強水文連接性，而較強的徑流、泥沙源區的連通性和坡面與溝道系統的有效連接有助于侵蝕產沙和泥沙輸送［31，48］。反之，匯單元在空間的相互連通及其控制區域與坡底溝道的連通則制約著坡面泥沙向溝道的輸送，可以有效減少坡面侵蝕的發生和出口輸沙量［31，48］。

對于源匯區的空間配置對水土流失的影響，可得到如下一般性規律:1)抑制侵蝕功能強的“匯”單元越靠近坡底，越有利于阻止泥沙的輸出。當以阻止泥沙輸出為目的時，應盡可能使抑制侵蝕功能強的單元靠近坡底分布［17，29，37，40］;2)“源”越破碎，越不利于侵蝕產沙，“源”連通性越強，侵蝕產沙越容易。當以防止坡面侵蝕產沙為目的時，應盡量采取分散徑流、泥沙“源”以降低其連通性的配置格局［30，36］;3)泥沙“源”距離出口越近，泥沙越容易輸出，與出口直接連通的“源”單元越多，坡面泥沙輸出量越大［31，48］;4)某單元接受上坡徑流供給越多，且在下坡路徑上遇到的阻礙越少，則該點泥沙越容易到達出口，對水土流失的貢獻越大［31，48］。

3 耦合覆被格局與水土流失的手段與方法

為定量表征覆被格局對水土流失的影響，建立兩者的耦合關系，目前的主要研究思路是:在對覆被格局進行準確描述基礎上，建立耦合水土流失過程的覆被格局指數和采用耦合格局信息的徑流與侵蝕模型。

3.1 覆被制圖

準確提取地表覆被信息，對覆被格局進行精確描述，獲取高分辨率的覆被分布圖是研究覆被格局對水土流失影響的重要前提。目前的坡面尺度覆被制圖方法主要有人工調查法［30］、照相法［27，49－50］和遙感影像法［17］。

人工調查方法工作量大，得到的覆被信息比較粗略，難以得到準確的覆被分布圖。照相法采用支架或氣球懸吊的方法將高分辨率相機固定在離地面一定高度(一般小于10 m)對坡面進行拍照，每張照片覆蓋的范圍可達到幾十平方米，大大節省了工作量，而且分辨率一般為幾厘米，要遠高于傳統的航拍方法。照片在GIS平臺上進行處理后得到覆被分布圖，可以精確獲得坡面各柵格單元的覆被類型和覆蓋度等重要信息，從而對覆被空間分布結構進行準確描述。高分辨率的遙感影像是開展大面積覆被制圖工作的有效途徑。Bartly等［17］利用Quickbird衛星影像資料(分辨率為2.4 m2)對澳大利亞北部半干旱區面積為11930，2031 m2和2861 m2的3個坡面進行覆被制圖和格局分析。

3.2 耦合水土流失過程的格局指數

景觀格局指數是描述覆被空間分布格局的有效途徑，是表達景觀格局－生態過程關系的重要工具。隨著景觀生態學理論和空間信息獲取與分析方法的發展，發展出眾多的景觀指數研究覆被格局與土壤侵蝕過程的關系，并基于景觀指數與土壤侵蝕過程變量關系的統計分析，探討土壤侵蝕對格局的響應。常用的格局指數包括空間和非空間指數，主要從圖形幾何或空間拓撲的角度描述覆被類型比例配置及其鑲嵌的空間結構、分布以及統計特征。劉宇等［51］闡述了FRAGSTATS 3.3軟件計算的邊界/斑塊密度、形狀、連接性、多樣性等4個方面12個常用的景觀指數在景觀格局－土壤侵蝕關系研究中的意義，分析了指數應用的局限性和原因，指出景觀數據屬性、景觀指數本身性質和土壤侵蝕過程的復雜性使得常規景觀指數與土壤侵蝕表征變量間存在不確定性。景觀指數變化的土壤侵蝕意義不明確，多種景觀指數在解釋土壤侵蝕狀況時會相互矛盾。以上局限性使得難以通過常規景觀指數來指示土壤侵蝕特征，而缺乏對土壤侵蝕過程機理的考慮是常規景觀指數在土壤侵蝕研究中存在不足的主要原因。因此，構建面向土壤侵蝕過程的格局指數將覆被空間組織結構與水土流失過程聯系起來并建立定量關系顯得尤為必要。

土壤侵蝕過程包括原位土壤顆粒的剝離(產沙過程)和泥沙通過輸送和沉積過程實現的再分配過程。裸地等徑流、泥沙產生區在水流路徑和泥沙輸送通道方向上的連通性對土壤侵蝕過程發揮著重要的控制作用，而植被等覆被類型則對徑流、泥沙則起著重要的阻滯作用。鑒于源匯區連通性在覆被格局影響土壤侵蝕過程中的關鍵作用，許多學者通過對水沙從源向匯的輸移通道空間連通性的描述，嘗試建立面向水土流失過程的格局指數來研究覆被格局－水土流失過程關系，將覆被格局特征與水土流失過程定量聯系起來，增強指數對土壤侵蝕過程的表征能力［48，51－52］。Jaeger［53］采用描述徑流、泥沙“源”的裸露斑塊分離度指數量化植被分布格局對產流、侵蝕能力的影響。植被斑塊發揮著攔截上坡來水來沙的功能，為了將植被分布格局與水土流失過程有效耦合，Imeson和Prinsen［52］認為植被－裸地鑲嵌格局可以指示水土流失源匯區域的幅度、空間分布和連通性，確立了空隙度、裸地斑塊破碎度、植被斑塊上坡坡長和裸地斑塊連通度4種格局指數定量化表述植被的空間分布狀況。該序列指標量化了裸地或植被斑塊即徑流泥沙源匯區的分布范圍和連通性，用來表征水土流失格局和過程。

Ludwig等［48］提出了基于水土流失過程描述植被斑塊空間分布的方向性滲透指數(DLI)，該指數通過徑流和泥沙的源、匯之間的歐氏距離來描述小區/坡面整體的物質留滯能力。DLI指數反映了植被斑塊間裸露區域的粒度和連接性，對植被斑塊的空間排列方式和形狀變化反應敏感，能較好地區分具有不同水土流失狀況的覆被格局［48］。Bautista等［30］在西班牙東南部(地中海半干旱區)的坡面徑流小區試驗研究表明，除極端暴雨情況外，徑流量與DLI指數具有較好的線性相關性，而產沙量與DLI呈指數函數關系，DLI指數可以較好地表征小區的水土流失狀況。但是，DLI指數忽略了植被類型之間的差別，將植被斑塊完全視為物質的匯，不能正確反映多種植被類型鑲嵌格局下的物質流失，而且DLI指數對柵格單元數量敏感，沒有考慮柵格到出口的距離對坡地徑流和泥沙的貢獻。

Mayor等［31］在坡面小區觀測試驗的基礎上，發展了綜合考慮植被空間分布和地形影響的匯流路徑長度指數(Flowlength)。該指數通過植被斑塊之間、地形洼地之間以及植被斑塊與地形洼地之間水流路徑長度的平均值來定量反映徑流路徑之間的連通性［31］。小區監測數據證實，Flowlength指數能夠區分具有相似覆蓋度的不同植被格局的連通性差異程度，指數隨裸地斑塊的粒度增加而增加，與產流和產沙量間均呈較好的線性相關性［31］。Flowlength指數體現了徑流、泥沙沿地表的運移過程和植被、地形洼地的阻滯功能，但將植被對徑流、泥沙的阻滯視為完全攔截。事實上，植被斑塊本身的特征(類型、蓋度等)決定了其攔截效率，更多時候徑流、泥沙會以一定的比例通過植被斑塊。此外，Flowlength指數僅僅反映裸露區域之間的連通程度，沒有考慮裸露區域的面積大小。

針對以上DLI指數和Flowlength指數的不足，劉宇［54］在DLI指數中引入反應柵格內斑塊的土壤保持功能及其與出口距離的權重系數，在Flowlength指數中為每個柵格單元引入反映覆被類型水土保持能力、坡度和植被覆蓋度的權重系數，更真實地反映了不同植被類型對徑流、泥沙的阻滯作用，建立改進的DLI和Flowlength指數，利用在黃土高原羊圈溝小流域開展的小區觀測試驗對指數的適用性進行了驗證，并模擬分析了不同坡面覆被空間配置格局的土壤侵蝕效應。

3.3 耦合格局信息的模型模擬

模型模擬是定量研究水土流失過程的重要手段。以往的坡面土壤侵蝕模型如 RUSLE/USLE［55］和WEPP［56］等大部分都是基于單一的土地利用類型，缺少對不同土地覆被格局水土流失效應的模擬研究。采取適當的方法，將坡面覆被格局信息融入到土壤侵蝕模型中，充分反映由于生物及非生物因素引起的植被以及土壤屬性空間分布結構對徑流泥沙等水文過程的影響，加強模型對土地覆被的表征，將會有效提高模型的預測精度，捕捉格局變化導致的水土流失效應變化，揭示水土流失對覆被格局變化的內在響應機理。這是耦合覆被格局與水土流失的重要手段，也是土壤侵蝕模擬的前沿方向，國際上已經開始了此方面的相關研究。

Muller等［57］以運動波坡面流模型和Smith入滲模型為基礎，根據實測結果采用統計學、地統計學和隨機方法對模型中植被覆蓋度、土壤入滲率和導水率等參數的非均勻分布和空間連接性特征進行定量表征，通過模型參數的空間分布反映覆被格局對水土流失的影響。模型采用的參數分布結果反映了坡面覆被格局的內在連接性特征，能對坡面徑流量和坡面產流的空間分布狀況進行較好的模擬，但模型結構比較復雜，需要較多點上的資料來構建模型參數的分布模型。

Arnau－Rosalén等［27］將西班牙地中海半干旱景觀的兩個南北向對比坡面的地表覆被圖劃分成規則的柵格單元并明確各單元的覆被類型，利用微型小區的人工模擬降雨試驗得到每種覆被類型的產流規律。在此基礎上，采用Horton模型計算每個柵格單元的入滲和產流量，并在整個坡面上進行連續的分布式模擬，得到地表徑流的空間分布，結果表明:由于坡面中部存在高覆蓋度的草本和裸露巖石碎片，形成徑流吸收區，阻礙了徑流路徑在整個坡面的連續性，使得徑流在上坡和下坡位匯流集中;南坡裸地、嵌入地中的巖石碎片、結皮等覆被類型具有較好的水文連通性，使得南坡產流區比例和產流量要高于北坡。該研究表明，將地表覆被制圖與徑流小區試驗有機結合，并利用經驗性模型對坡面產流進行分布式模擬，是研究坡面覆被格局對徑流、侵蝕影響的簡單有效手段。

Boer和Puigdefábregas［32］在100 m×100 m的假定坡面上，設定植被的平均蓋度、標準差和自相關長度等指標，根據非條件高斯模擬得到坡面上植被蓋度的空間分布，并將具有相同平均蓋度但均勻分布的坡面作為對照，采用基于過程的空間分布式模型LISEM對次降雨事件下坡面的徑流和產沙過程進行模擬，結果表明坡面植被覆蓋和土壤屬性的空間分布結構化特征對水土流失具有重要影響。在模擬過程中，為建立植被和土壤參數的協同關系，假定一些重要參數如曼寧系數、土壤孔隙度、飽和導水率等與植被覆蓋度呈簡單的線性關系，而且僅僅考慮植被蓋度的空間分布，沒有考慮不同植被類型之間的差異性。

Frot和van Wesemael［58］利用基于專家的分布式模型STREAM對西班牙半干旱區3個不同覆被格局草地坡面的產流空間分布和出口徑流量進行模擬。該模型根據植被和裸巖兩者覆蓋度的不同組合劃分不同等級的生態水文響應單元，對坡面進行柵格單元劃分并確定其水文響應單元類型，對每個柵格單元建立徑流和入滲之間的平衡方程，再根據DEM進行匯流計算得到出口徑流量。每一類水文響應單元的穩滲率和初損值根據研究區的相關實驗結果，采用專家分析和模型率定的方法確定。STREAM模型反映了不同水文響應單元的植被覆蓋度和土壤入滲特性的空間差異，對坡面水文非均質性能較好地描述，模擬結果與實測值吻合較好。但模型沒有考慮不同植被類型的水文性質差異，而且模型地域性很強，即使屬于同一氣候帶且土壤質地相似的區域，模型也需要重新校正和檢驗后才能應用［59］。

分析以上相關進展可以看出，目前耦合坡面覆被格局與水土流失的模型均采用分布式模擬的方式，對產流和產沙過程的描述既有經驗模型，也有物理過程模型，而對覆被格局影響的表征主要是通過在模型中描述植被以及土壤屬性參數的空間分布特征和不同水文響應單元的連通性來實現，水文響應單元也主要是根據覆被特征劃分。以上這種隱式表達覆被格局影響的模擬方法要得到比較準確的模擬結果，對模型參數和支撐數據的要求較高。如果能將反映覆被格局的參數直接嵌入到模型中，可以將覆被格局信息與模型更加簡單地耦合起來，對于研究覆被格局與水土流失過程關系具有重要意義。

4 問題與展望

針對干旱半干旱區坡面覆被格局與水土流失關系，國內外學者以斑塊尺度上植被的產流產沙特征觀測為基礎，重點分析了植被類型、植被的垂直結構和形態特征的影響;基于實驗觀測研究了坡面上不同覆被類型、裸地－植被斑塊鑲嵌格局和植被條帶格局的水土流失規律，探討了坡面植被空間分布格局和徑流泥沙源匯區的連通性對水土流失的影響;為表征覆被格局的影響，以坡面的精確覆被制圖為基礎，通過基于水土流失過程的覆被格局指數和耦合覆被格局信息的徑流與侵蝕模型來定量研究覆被格局的水土流失效應。本文對需要進一步研究的熱點問題和未來的發展趨勢總結如下:

(1)覆被格局的動態變化

植被生長的季節動態變化和年際間變化會伴隨著植被形態特征及其分布格局的改變，對水文響應有重要影響。另一方面，在干旱半干旱區，降雨和地形等因素的綜合作用會使坡面植被斑塊等水土流失匯功能單元的強度發生變化，甚至轉化成源，具有動態性［60］。因此，無論是從植被格局的分布特征還是從它對水土流失的功能來看，坡面覆被格局都存在一定的動態變化規律。以往的研究大多從靜態的角度去分析覆被格局對水土流失過程的控制作用，較少在長期監測基礎上以動態系統的觀點研究格局的變化特征及其對水土流失造成的影響。

(2)覆被格局與水土流失過程的相互作用

覆被格局影響水土流失過程，而徑流泥沙產生可以改變土壤特性及微地理環境，進而反作用于植被生長、更新、演替及植被分布，引起覆被格局變化。覆被格局與水土流失過程構成一個相互作用和相互影響的反饋調節系統。以往的研究主要集中在植被對水土流失的控制上，而關于水土流失對植被及其分布格局影響的研究較少，兩者之間的相互作用關系需要加強研究。

(3)尺度問題

景觀格局對水土流失的作用具有尺度依賴性。在不同尺度上，景觀格局與水土流失過程相互作用的主要影響因素不同，在斑塊尺度上主要是植被特性和土壤屬性，在坡面尺度上主要是覆被分布格局和地形因子，在流域尺度上更為復雜，主要有植被、土壤分布和土地利用方式以及地形、氣象和人類活動干擾等。本文所綜述的僅是坡面尺度上的相關研究工作。從小區、坡面尺度到流域或區域尺度，水土流失有著既相互關聯又有明顯區別的生物物理控制過程，使得各尺度間具有連續性，但由于徑流的非連續性以及地表的異質性等復雜因素影響，小區、坡面的觀測結果很難直接上推至流域尺度。因此，需要系統理解植被與徑流泥沙之間的聯系與反饋機制及其隨尺度變化的規律性，明確不同尺度上水土流失過程的關鍵影響因子和驅動因素，發展尺度上推方法。

(4)覆被格局與水土流失的耦合模型

基于水土流失過程的覆被格局指數與耦合格局信息的水土流失模型是定量分析覆被格局與水土流失關系的有效手段，這兩方面的研究剛剛起步，是未來的重要發展趨勢。構建過程意義明確且簡單實用的格局指數，將覆被格局與地理－水文過程聯系起來，這是景觀格局－土壤侵蝕相互作用研究在方法上的需求。對模型而言，應以分布式水文模型為基礎，將覆被格局分布及各覆被類型參數耦合到模型中，將覆被格局的動態信息與產流產沙過程相結合，加強對格局與徑流泥沙反饋系統的耦合，建立真正意義上的覆被格局－水土流失過程耦合模型。

［1］ Leng S Y，Song C Q.Review of land surface geographical process study and prospects in China.Advances in Earth Science，2005，20(6):600－606.

［2］ Fu B J，Zhao W W，Chen L D.Progress and perspective of geographical－ecological processes.Acta Geographica Sinica，2006，61(11):1123－1131.

［3］ Pimentel D，Harvey C，Resosudarmo P，Sinclair K，Kurz D，McNair M，Crist S，Shpritz L，Fitton L，Saffouri R，Blair R.Environmental and economic costs of soil erosion and conservation benefits.Science，1995，267(5201):1117－1123.

［4］ Fu B J，Chen L D，Qiu Y，Wang J，Meng Q H.Land Use Structure and Ecological Process in the Hilly Area of Loess Plateau，China.Beijing:The Commercial Press，2002.

［5］ Xu X L，Ma K M，Fu B J，Liu X C，Huang Y，Qi J.Research review of the relationship between vegetation and soil loss.Acta Ecologica Sinica，2006，26(9):3137－3143.

［6］ Tongway D J，Valentin C，Seghieri J.Banded Vegetation Patterning in Arid and Semiarid Environments:Ecological Processes and Consequences for Management.New York:Springer，2001.

［7］ Boix－Fayos C，Martínez－Mena M，Arnau－Rosalén E，Calvo－Cases A，Castillo V，Albaladejo J.Measuring soil erosion by field plots:understanding the sources of variation.Earth－Science Reviews，2006，78(3/4):267－285.

［8］ Gallego Fernández J B，Rosario García Mora M，Novo F G.Vegetation dynamics of Mediterranean shrublands in former cultural landscape at Grazalema Mountains，South Spain.Plant Ecology，2004，172(1):83－94.

［9］ Garcia－Estringana P，Alonso－Blázquez N，Marques M J，Bienes R，Alegre J.Direct and indirect effects of Mediterranean vegetation on runoff and soil loss.European Journal of Soil Science，2010，61(2):174－185.

［10］ Zhang Z Q，Wang S P，Sun G，Xie B Y.Runoff and sediment yield response to vegetation change at multiple scales:a review.Acta Ecologica Sinica，2006，26(7):2356－2364.

［11］ Ludwig J A，Wilcox B P，Breshears D D，Tongway D J，Imeson A C.Vegetation patches and runoff－erosion as interacting ecohydrological processes in semiarid landscapes.Ecology，2005，86(2):288－297.

［12］ Casermeiro M A，Molina J A，de la Cruz Caravaca M T，Hernando Costa J，Hernando Massanet M I，Moreno P S.Influence of scrubs on runoff and sediment loss in soils of Mediterranean climate.Catena，2004，57(1):91－107.

［13］ Sun J J，Yu D S，Shi X Z，Gu Z J，Zhang W T，Yang H.Comparison of between LAI and VFC in relationship with soil erosion in the red soil hilly region of south China.Acta Pedologica Sinica，2010，47(6):1060－1066.

［14］ Zhang G H，Liang Y M.A summary of impact of vegetation coverage on soil and water conservation benefit.Research of Soil and Water Conservation，1996，3(2):104－110.

［15］ Vásquez－Méndez R，Ventura－Ramos E，Oleschko K，Hernández－Sandoval L，Parrot J F，Nearing M A.Soil erosion and runoff in different vegetation patches from semiarid Central Mexico.Catena，2010，80(3):162－169.

［16］ Moreno－de las Heras M，Merino－Martín L，Nicolau J M.Effect of vegetation cover on the hydrology of reclaimed mining soils under Mediterranean－Continental climate.Catena，2009，77(1):39－47.

［17］ Bartley R，Corfield J P，Abbott B N，Hawdon A A，Wilkinson S N，Nelson B.Impacts of improved grazing land management on sediment yields，Part 1:hillslope processes.Journal of Hydrology，2010，389(3/4):237－248.

［18］ Zhu B B，Li Z B，Li P，You Z.Effect of grass coverage on sediment yield of rain on slope.Acta Pedologica Sinica，2010，47(3):401－407.

［19］ Jiao J Y，Wang W Z，Li J.Effective cover rate of woodland and grassland for soil and water conservation.Acta Phytoecologica Sinica，2000，24(5):608－612.

［20］ Morgan R P C.Vegetative－based technologies for erosion control//Skokes A，Spanos I，Norris J E，Cammeraat E，eds.Eco－and Ground Bio－Engineering:The use of Vegetation to Improve Slope Stability.Springer，2007:265－272.

［21］ Cai Q G，Chen H.The effect of vegetation cover on reduction of rainfall erosion//Chen Y Z，ed.The Source of the Coarse Sediment in the Middle Yellow River and the Mechanism of Soil and Sediment Production.Beijing:Meteorology Press，1989.

［22］ He S X，Li X Y，Mo F，Zhou B，Gao G L.The water conservation study of typical forest ecosystems in the forest transect of eastern China.Acta Ecologica Sinica，2011，31(12):3285－3295.

［23］ Thurow T L，Blackburn W H，Taylor C A Jr.Hydrologic characteristics of vegetation types as affected by livestock grazing systems，Edwards Plateau，Texas.Journal of Range Management，1986，39(6):505－509.

［24］ Gyssels G，Posen J，Nachtergaele J，Govers G.The impact of sowing density of small grains on rill and ephemeral gully erosion in concentrated flow zones.Soil and Tillage Research，2002，64(3/4):189－201.

［25］ Gyssels G，Poesen J.The importance of plant root characteristics in controlling concentrated flow erosion rates.Earth Surface Processes and Landforms，2003，28(4):371－384.

［26］ Li Y.Vegetation Root and Soil Anti－Scourability in Loess Plateau.Beijing:Science Press，1995.

［27］ Arnau－Rosalén E，Calvo－Cases A，Boix－Fayos C，Lavee H，Sarah P.Analysis of soil surface component patterns affecting runoff generation:an example of methods applied to Mediterranean hillslopes in Alicante(Spain).Geomorphology，2008，101(4):595－606.

［28］ Puigdefábregas J.The role of vegetation patterns in structuring runoff and sediment fluxes in drylands.Earth Surface Processes and Landforms，2005，30(2):133－147.

［29］ Bartley R，Roth C H，Ludwig J，McJannet D，Liedloff A，Corfield J，Hawdon A，Abbott B.Runoff and erosion from Australia's tropical semi－arid rangelands:influence of ground cover for differing space and time scales.Hydrological Processes，2006，20(15):3317－3333.

［30］ Bautista S，Mayor á G，Bourakhouadar J，Bellot J.Plant spatial pattern predicts hillslope runoff and erosion in a semiarid Mediterranean landscape.Ecosystems，2007，10(6):987－998.

［31］ Mayor á G，Bautista S，Small E E，Dixon M，Bellot J.Measurement of the connectivity of runoff source areas as determined by vegetation pattern and topography:a tool for assessing potential water and soil losses in drylands.Water Resources Research，2008，44:W10423，doi:10.1029/2007WR006367.

［32］ Boer M，Puigdefábregas J.Effects of spatially structured vegetation patterns on hillslope erosion in a semiarid Mediterranean environment:a simulation study.Earth Surface Processes and Landforms，2005，30(2):149－167.

［33］ Boix－Fayos C，Martínez－Mena M，Calvo－Cases A，Arnau－Rosalén E，Albaladejo J，Castillo V.Causes and underlying processes of measurement variability in field erosion plots in Mediterranean conditions.Earth Surface Processes and Landforms，2007，32(1):85－101.

［34］ Reid K D，Wilcox B P，Breshears D D，MacDonald L.Runoff and erosion in a Pi?on－Juniper woodland:influence of vegetation patches.Soil Science Society of America Journal，1999，63(6):1869－1879.

［35］ Shachak M，Sachs M，Moshe I.Ecosystem management of desertified shrublands in Israel.Ecosystems，1998，1(5):475－483.

［36］ Ludwig J A，Bartley R，Hawdon A A，Abott B N，McJannet D.Patch configuration non－linearly affects sediment loss across scales in a grazed catchment in north－east Australia.Ecosystems，2007，10(5):839－845.

［37］ Gumiere S J，Bissonnais Y L，Raclot D，Cheiron B.Vegetated filter effects on sedimentological connectivity of agricultural catchments in erosion modelling:a review.Earth Surface Processes and Landforms，2011，36(1):3－19.

［38］ Ludwig J A，Tongway D J，Marsdem S G.Stripes，strands or stipples:modelling the influence of three landscape banding patterns on resource capture and productivity in semi－arid woodlands，Australia.Catena，1999，37(1/2):257－273.

［39］ Valentin C，d'Herbès J M，Poesen J.Soil and water components of banded vegetation patterns.Catena，1999，37(1/2):1－24.

［40］ Rey F.Effectiveness of vegetation barriers for marly sediment trapping.Earth Surface Processes and Landforms，2004，29(9):1161－1169.

［41］ Martníez Raya A，Durán Zuazo V H，Francia Martínez J R.Soil erosion and runoff response to plant－cover strips on semiarid slopes(SE Spain).Land Degradation and Development，2006，17(1):1－11.

［42］ Su M，Lu Z F，Zhang X C，Li G X.The appraisal and analysis on benefits of water and soil conservation under the different planting pattern in loess hilly area.Bulletin of Soil and Water Conservation，1990，10(4):46－52.

［43］ You Z，Li Z B，Jiang Q F.Study on the effect of vegetation patterns on the slope on the rainfall erosion.Journal of Sediment Research，2005，(6):40－43.

［44］ Li M，Yao W Y，Chen J N，Ding W F，Yang J F，Li L，Yang C X.Impact of different grass coverages on the sediment yield process in the slopegully system.Acta Geographica Sinica，2005，60(5):725－732.

［45］ Shen Z Y.Study on the Effect of Vegetation Slope Pattern on Soil and Water Loss［D］.Shaanxi Xi'an:Xi'an University of Technology，2006.

［46］ Xu H Y，Zhao W W，Zhu H F，Zhu J.Characteristics of soil loss under different land use combination of sloping farmland and grassland in the hilly－gully region of Loess Plateau.Science of Soil and Water Conservation，2009，7(3):35－41.

［47］ Fu B J，Wang Y F，Lu Y H，He C S，Chen L D，Song C J.The effects of land－use combinations on soil erosion:a case study in the Loess Plateau of China.Progress in Physical Geography，2009，33(6):793－804.

［48］ Ludwig J A，Eager R W，Bastin G N，Chewings V H，Liedloff A C.A leakiness index for assessing landscape function using remote sensing.Landscape Ecology，2002，17(2):157－171.

［49］ Lesschen J P，Cammeraat L H，Kooijman A M，van Wesemael B.Development of spatial heterogeneity in vegetation and soil properties after land abandonment in a semi－arid ecosystem.Journal of Arid Environments，2008，72(11):2082－2092.

［50］ Ruiz Sinoga J D，Martinez Murillo J F.Hydrological response of abandoned agricultural soils along a climatological gradient on metamorphic parent material in southern Spain.Earth Surface Processes and Landforms，2009，34(15):2047－2056.

［51］ Liu Y，Lü Y H，Fu B J.Implication and limitation of landscape metrics in delineating relationship between landscape pattern and soil erosion.Acta Ecologica Sinica，2011，31(1):267－275.

［52］ Imeson A C，Prinsen H A M.Vegetation patterns as biological indicators for identifying runoff and sediment source and sink areas for semi－arid landscapes in Spain.Agriculture Ecosystems and Environment，2004，104(2):333－342.

［53］ Jaeger J A G.Landscape division，splitting index，and effective mesh size:new measures of landscape fragmentation.Landscape Ecology，2000，15(2):115－130.

［54］ Liu Y.Effect of Landscape Pattern Dynamics on Soil Erosion［D］.Beijing:The Graduate School of the Chinese Academy of Sciences，2011.

［55］ Renard G R，Foster G R，Weesies G A.Prediction rainfall erosion by water:A guide to conservation planning with the revised universal soil loss equation(RUSLE).USDA Agricultural Handbook，No.703，1997.

［56］ Flanagan D C，Gilley J E，Franti T G.Water Erosion Prediction Project(WEPP):development history，model capabilities，and future enhancements.Transactions of the American Society of Agricultural and Biological Engineers，2007，50(5):1603－1612.

［57］ Mueller E N，Wainwright J，Parsons A J.Impact of connectivity on the modeling of overland flow within semiarid shrubland environments.Water Resources Research，2007，43:W09412，doi:10.1029/2006WR005006.

［58］ Frot E，van Wesemael B.Predicting runoff from semi－arid hillslopes as source areas for water harvesting in the Sierra de Gador，southeast Spain.Catena，2009，79(1):83－92.

［59］ Evrard O，Cerdan O，van Wesemael B，Cauvet M，Le Bissionnais Y，Raclot D，Vandaele K，Andrieux P，Bielders C.Reliability of an expertbased runoff and erosion model:application of STREAM to different environments.Catena，2009，78(2):129－141.

［60］ Wainwright J，Parsons A J，Schlesinger W H，Abrahams A D.Hydrology－vegetation interactions in areas of discontinuous flow on a semi－arid bajada，Southern New Mexico.Journal of Arid Environments，2002，51(3):319－338.

參考文獻:

［1］ 冷疏影，宋長青.陸地表層系統地理過程研究回顧與展望.地球科學進展，2005，20(6):600－606.

［2］ 傅伯杰，趙文武，陳利頂.地理－生態過程研究的進展與展望.地理學報，2006，61(11):1123－1131.

［4］ 傅伯杰，陳利頂，邱揚，王軍，孟慶華.黃土丘陵溝壑區土地利用結構與生態過程.北京:商務印書館，2002.

［5］ 徐憲立，馬克明，傅伯杰，劉憲春，黃勇，祁建.植被與水土流失關系研究進展.生態學報，2006，26(9):3137－3143.

［10］ 張志強，王盛萍，孫閣，謝寶元.流域徑流泥沙對多尺度植被變化響應研究進展.生態學報，2006，26(7):2356－2364.

［13］ 孫佳佳，于東升，史學正，顧祝軍，張文太，楊浩.植被葉面積指數與覆蓋度定量表征紅壤區土壤侵蝕關系的對比研究.土壤學報，2010，47(6):1060－1066.

［14］ 張光輝，梁一民.植被蓋度對水土保持功效影響的研究綜述.水土保持研究，1996，3(2):104－110.

［18］ 朱冰冰，李占斌，李鵬，游珍.草本植被覆蓋對坡面降雨徑流侵蝕影響的試驗研究.土壤學報，2010，47(3):401－407.

［19］ 焦菊英，王萬中，李靖.黃土高原林草水土保持有效蓋度分析.植物生態學報，2000，24(5):608－612.

［21］ 蔡強國，陳浩.植被覆蓋對降雨濺蝕的影響 //陳永宗.黃河粗泥沙來源及侵蝕產沙機理研究文集.北京:氣象出版社，1989:41－47.

［22］ 賀淑霞，李敘勇，莫菲，周彬，高廣磊.中國東部森林樣帶典型森林水源涵養功能.生態學報，2011，31(12):3285－3295.

［26］ 李勇.黃土高原植物根系與土壤抗沖性.北京:科學出版社，1995.

［42］ 蘇敏，盧宗凡，張興昌，李夠霞.黃土丘陵區不同種植方式水保效益的分析評價.水土保持通報，1990，10(4):46－52.

［43］ 游珍，李占斌，蔣慶豐.坡面植被分布對降雨侵蝕的影響研究.泥沙研究，2005，(6):40－43.

［44］ 李勉，姚文藝，陳江南，丁文峰，楊劍鋒，李莉，楊春霞.坡面草被覆蓋對坡溝侵蝕產沙過程的影響.地理學報，2005，60(5):725－732.

［45］ 沈中原.坡面植被格局對水土流失影響的實驗研究［D］.西安:西安理工大學，2006.

［46］ 徐海燕，趙文武，朱恒峰，朱婧.黃土丘陵溝壑區坡耕地與草地不同配置方式的侵蝕產沙特征.中國水土保持科學，2009，7(3):35－41.

［51］ 劉宇，呂一河，傅伯杰.景觀格局－土壤侵蝕研究中景觀指數的意義解釋及局限性.生態學報，2011，31(1):267－275.

［54］ 劉宇.景觀格局變化的土壤侵蝕效應研究［D］.北京:中國科學院研究生院，2011.