巴基斯坦波特瓦爾高原土壤侵蝕調查報告

陳同德， 朱夢陽， 焦菊英，3， 安韶山，3， 稅軍峰，3， 楊 陽， Fakher Abbas

(1.西北農林科技大學 水土保持研究所 黃土高原土壤侵蝕與旱地農業國家重點實驗室， 陜西 楊凌 712100；2.西北大學 城市與環境學院， 陜西 西安 710127； 3.中國科學院 水利部 水土保持研究所， 陜西 楊凌 712100)

泛第三極面積約2.00×107km2，人口多達3.00×108多，目前面臨氣溫急劇升高[1]、冰川融化[2]和荒漠化加劇[1]等環境問題。泛第三極區域與中國提出的“一帶一路”倡議經濟帶高度重合，了解、研究和解決該區域的資源環境問題是在該倡議實施過程中的重要環節。巴基斯坦是泛第三極南亞國家中重要的一員，調查分析巴基斯坦區域土壤侵蝕規律和治理狀況，將對該國土壤侵蝕防治以及服務泛第三極環境治理具有重要意義。

近年來巴基斯坦人口急劇增加。聯合國人口規劃和預測委員會預測，至2025年，巴基斯坦人口將高達2.60×108人[3]，雖然是南亞地區的農業大國，農業對GDP的貢獻達到25%[4]，但巴基斯坦的耕地僅為2.00×107km2，人口和農耕地均集中分布于印度河及其支流兩岸[5]；而且其76%的土地受到不同程度土壤侵蝕影響，已對該國畜牧業和生態安全造成巨大威脅[6]。人為因素方面，人口增加后毀林開荒、過度放牧和不合理的耕作等活動[6]，加劇土壤侵蝕；在自然因素方面，印度河流域整體降雨量較大，土壤疏松[7]，因此容易發生水蝕。土壤侵蝕不僅引起農耕地土壤表層的營養物質持續不斷地流失，使土地沙化，也會污染下游水體[7]。

目前巴基斯坦的土壤侵蝕研究，是以遙感結合RUSLE等模型分析小流域或者區域尺度的土壤侵蝕時空變化和區域分異規律為主[8]，僅有個別小流域有部分年份的土壤侵蝕觀測數據[9]，并不能完全支撐模型所得結果。另外，巴基斯坦氣候變化劇烈，經濟高速發展，導致區域環境變化迅速，其土壤侵蝕發生的條件、形式和程度發生了改變。因此，本文采用野外調查和3S技術結合的方式，對巴基斯坦的土壤侵蝕現狀、區域分異規律和侵蝕影響因子等方面進行分析，旨在研究該區土壤侵蝕防治對策，為認識泛第三極區域生態環境變化及其治理提供科學依據。

1 調查區域概況

巴基斯坦(23°45′—36°56′N，60°53′—71°03′E)西北部與中國新疆接壤，其東部為印度、西部為阿富汗。巴基斯坦國土面積約8.81×105km2(含巴控克什米爾地區)，南北長約1 600 km，擁有平原、丘陵、高原和山地等多種地貌類型。巴基斯坦大部分地區屬于熱帶干旱和半干旱氣候，處于熱帶季風區西緣，年降雨量為50～1 000 mm，降雨分布空間不均，北部地區年降雨可達2 000 mm[10]。巴基斯坦地形多變，可分為北部高山區、波特瓦爾高原、旁遮普平原、俾路支高原和信德平原等單元，侵蝕地貌也有多種類型(附圖18a—18d)。其中，北部高山區是巴基斯坦最高的區域，也是印度河、杰盧姆河和齊納布河等主要河流的發源地；高山區南部為波特爾高原(32°13′—34°22′N，71°18′—73°57′E)，面積約220 km2，高原中部由平原區和溝壑區組成(附圖19a—19b)，坡度范圍為8%～40%，是巴基斯坦重要的糧食生產地[11]，也是政治文化中心(伊斯蘭堡位于波特爾高原北部)。波特瓦爾高原降雨年際分布不均，約50%～80%的降雨分布于5—10月[11]，因此降雨侵蝕力較大，水土流失嚴重，是巴基斯坦水蝕最嚴重的地區(附圖19c—19d)，總侵蝕量約4.24×108t/a，其中西部(Attock)和東部地區(Jhelum)年平均侵蝕量約8.00×107t/a，而北部地區(Rawalpindi)高達1.60×108t/a[9]；波特瓦爾高原北部區域，即使是覆蓋度較高的林區，部分區域的侵蝕也較嚴重(附圖20a—20d)，另外該區域的塔貝拉水庫，是伊斯蘭堡城市用水以及附近農業用水的主要來源地，印度河每天向該水庫傾注的泥沙量可達50 000 t，嚴重縮短水庫的使用壽命和儲蓄能力[7]。

2 調查內容與方法

2.1 樣點布設及調查路線

波特瓦爾高原伊斯蘭堡附近人口分布密集，地表環境受人為干擾程度較重，土地利用變化較快，因此，基于年降雨量變化(500～1 400 mm)、土地利用的多樣性和地形等。在伊斯蘭堡東北—西南方向布設樣線，以10 km左右的間距選擇了15個面積在0.74～2.14 km2的調查單元(12個小流域，3 個矩形區)進行野外調查(圖1)，調查單元的基本信息如表1所示。

圖1 巴基斯坦波特瓦爾高原土壤侵蝕調查路線和調查單元分布位置表1 波特瓦爾高原土壤侵蝕調查單元基本信息

調查單元編號調查日期年均降雨量/mm海拔高度/m地貌特征主要土地利用120190420520600～700山區,多砂巖林地、耕地220190420500517～526高原區、地勢平坦耕地、建設用地320190420600524～554高原區、調查單元內有寬淺溝壑,面積約占40%耕地、草地420190416640498～525高原區,調查單元東部有溝壑,面積約占15%耕地520190420720508～528高原區耕地、建設用地620190419800517～575高原區、調查單元內有寬淺溝壑,面積約占45%耕地、草地720190419820502～535高原區、調查單元內有寬淺溝壑,面積約占20%耕地、草地820190419800487～524高原區、調查單元內有寬淺溝壑,面積約占5%耕地、草地920190419900543～580高原區、調查單元內多寬淺溝壑,面積約占35%耕地、建設用地1020190418900438～532丘陵區林地、草地11201904171 000499～520河流階地耕地、草地12201904171 100592～816山地區林地、草地13201904181 2001 111～1 385山地區耕地、林地14201904181 2001 161～1 755山地區林地、草地、耕地15201904181 3001 425～1 928山地區林地、草地、耕地

注：降雨量數據來自于Fick, Stephen E和Hijmans, Robert J[12]。

2.2 調查方法和內容

根據第一次中國水利普查水土保持普查中的土地利用的分類標準[13]，利用前期解譯的土地利用和水保措施圖手動記錄和奧維互動地圖APP定位照相，記錄調查單元每個圖斑的土地利用情況，后期基于ArcGIS進行修編，解譯結果如圖1所示。同時，調查并記錄不同地類圖斑的地理位置、土壤侵蝕特征、植被類型與蓋度和農耕地利用現狀等信息，調查共計15個調查單元(圖1)共1 026個圖斑。然后，利用中國土壤流失方程(CSLE)計算每個單元中各圖斑的土壤侵蝕速率t/(km2·a)，以此分析各抽樣單元的土壤侵蝕程度以及空間分布特征。中國土壤侵蝕模型(CSLE)是在美國通用土壤流失方程(USLE)的基礎上，根據中國土壤侵蝕特征、試驗觀測數據和治理技術方法等改進模型。CSLE將USLE的LS算法改為適用于陡坡的算法，將植被覆蓋和管理因子(C)和水土保持措施因子(P)改為水土保持生物措施(B)、工程措施(E)和耕作措施(T)3個因子[14-16]。CSLE已在中國新疆[17](與巴基斯坦接壤)、陜西省[18]和其它不同地形、氣候區[19-21]以及泰國[22]等泛第三極部分區域有成功應用，因此本次調查抽樣單元的土壤侵蝕侵蝕模數計算使用該模型。

2.3 CSLE各因子和流失速率計算方法

本文采用中國第1次水力普查中CSLE各因子的計算方法[13]，計算區域土壤侵蝕力因子(R)和土壤可蝕性因子(R)；利用1弧秒分辨率的SRTM數據和LS_Tool計算工具計算坡度坡長因子(LS)；基于NDVI和PVI數據以及實地調查數據，計算水土保持措施因子(B，E，T)；最終通過CSLE計算得到土壤侵蝕速率(A)。

3 土壤侵蝕特征與土壤侵蝕治理措施

3.1 土地利用類型土壤侵蝕特征

3.1.1 林地土壤侵蝕特征 林地主要分布在巴基斯坦東北部的穆里地區(第13—15調查單元)，喬木林地(主要為馬尾松)可占調查單元總面積的28.87%～69.90%，平均郁閉度在50%以上。平均土壤侵蝕速率在248.17～723.44(t/(km2·a)之間(表2)。自然條件下土壤侵蝕程度較輕，僅部分林間地有片狀侵蝕(圖2a)，部分山坡坡腳發生重力侵蝕(圖2b)；但在生產建設項目區，裸露邊坡侵蝕嚴重(圖2c)，甚至產生大切溝(圖2d，切溝長約20 m，寬1～5 m，深1～3 m)。林間地片狀侵蝕產生原因是馬尾松林地間的草被稀疏，無法對地表實現全覆蓋，降雨匯集后產生漫流，無草被覆蓋的區域，地表的抗侵蝕能力相對較弱，因此發生較大面積剝蝕，甚至在下方產生與坍塌。草被稀疏可能的原因有兩個方面：①此處的年均降雨量在1 200～1 300 mm(表1)，因此馬尾松生長繁茂巨大(個別樹高達40 m多)，馬尾松林與林下草本植物爭奪養分和水分，導致林間草被稀疏；②已有研究證實馬尾松會分泌化感物質，會影響附近的其它植物生長[23]。

圖2 穆里地區林地土壤侵蝕(陳同德攝于20190418)

生產建設項目區域，因植被遭到破壞，土壤失去保護層，在降雨和徑流的直接打擊和沖刷作用下，侵蝕強烈，產生細溝侵蝕和切溝侵蝕，土壤細顆粒物質隨徑流搬運輸移，地下礫石逐漸出露，對自然景觀造成破壞。

坡腳發生重力侵蝕的原因主要包括兩個：①巴基斯坦北部為喜馬拉雅山脈上升地區[24]，因此坡面長、坡度大，應力集中于坡腳，容易發生崩塌；②該區域降雨量大，坡腳處為徑流匯集的區域，水分充足，土層之間的黏結力因水分作用而降低，易在水分作用下發生滑塌。

3.1.2 草地土壤侵蝕特征 在調查范圍內，草地主要分布在古杰爾汗—杰格瓦爾地區(第3—9調查單元)，其面積可占調查單元總面積的16.70%～37.75%，平均土壤侵蝕速率在9.93～306.47 t/(km2·a)之間(表2)。該區域年均降雨量498～580 mm，屬半濕潤氣候。自然條件下草被生長較好，侵蝕程度也較輕，坡面分布片狀侵蝕。是由于在降雨和徑流作用下，無植被覆蓋的區域土壤細顆粒不斷被運移，礫石逐漸出露，出現片狀侵蝕。另外，地表在長期侵蝕作用下，也會產生侵蝕溝道(圖3a)。在人為活動影響下，草地產生較為嚴重的土壤侵蝕，土質道路邊坡，因無植被措施防護，土質松散，降雨容易匯集，產生大量細溝侵蝕，部分區域出現切溝(圖3b)；高草被蓋度的自然山坡因修路被切割，沒有草被防護的區域，失去了植被對土壤的固結作用，不斷發生泄溜、撒落等重力侵蝕，與有草被區域形成鮮明對比(圖3c)；放牧活動可對草被造成嚴重破壞，食草動物啃食植物莖葉，影響植物生長，使地表覆蓋度降低，踩踏使地表土壤緊實，其儲水性能降低，降雨條件下，更容易形成徑流，出現重力侵蝕和切溝侵蝕(圖3d)。

圖3 古杰爾汗-杰格瓦爾地區草地(稀樹)土壤侵蝕(焦菊英攝于20190418-20190419)

3.1.3 建設用地土壤侵蝕特征 巴基斯坦的人口、經濟近年都在飛速增長中[4]，不可避免地會進行開發建設。在本次調查過程中，80%的調查單元均有建設項目，平均土壤侵蝕速率在17.09～1 711.76(t/(km2·a)之間(表2)。在各種項目的建設過程中，會使地表植被遭到破壞，土體遭到嚴重干擾，極易產生水土流失。比如，在第10調查單元(圖1)東側，分布有約5 km2的建設用地，地表原生植被均已被破壞殆盡；排土邊坡(附圖21a)植被稀疏，土質松散，發生了嚴重的細溝、切溝侵蝕，另外，排土邊坡呈多級分布，每一級平臺對下一級而言形成匯水區域，因此越往下的排土邊坡更容易產生股流，形成的細溝、切溝數量越多，深度、寬度也較大；路旁棄土形成的地表植被蓋度雖然較高，但卻形成眾多較大切溝(圖4a，切溝長約30 m，深0.5～2 m，寬0.5～5 m)，主要是由于棄土攤開后形成的地表坡面較長，形成的徑流量較大。加之土質松軟，草被層無法完全吸收水分、分散徑流，形成較大切溝，對當地的土地資源造成嚴重破壞。在第5調查單元(圖1)內，分布有小磚廠群，制磚過程中會大量取土，破壞植被，使得地表裸露，抗侵蝕能力降低，產生切溝侵蝕(附圖21b)；因取土過程中地表不斷被車輛壓實，取土平臺土壤的入滲能力大大降低，取土平臺本身構成較大面積的匯水區域，在降雨條件下，取土平臺邊坡易產生細溝、切溝侵蝕(圖4b)。除以上建設類型，推平地和高壓塔(附圖21c—21d)建設等也有不同程度的土壤侵蝕產生，對當地的水土資源造成了破壞。

圖4 伊斯蘭堡首都區旁遮普地區和門德拉北部地區建設用地侵蝕溝

3.2 調查單元的土壤侵蝕速率

利用CSLE計算得出15個調查單元土壤侵蝕速率的平均值介于51.40～1 338.91 t/(km2·a)之間(表2)。Muhammad等人[7]在恰克瓦爾地區(位于波特瓦爾高原南部)兩個小流域的實測的土壤侵蝕速率變化范圍為479～1 231 t/(km2·a)。S. Ellis等人[25]在穆里地區(位于波特瓦爾高原北部)高海拔砍伐林地和較低海拔過度放牧地實測土壤侵蝕速率約為15 000 t/(km2·a)和500～7 500 t/(km2·a)，這2組數據代表受人為干擾較大的地區的侵蝕速率。目前，波特瓦爾高原僅有的土壤侵蝕實測值約在479～15 000 t/(km2·a)，可見，CSLE模型計算結果基本在這個范圍內。因此，CSLE模型計算結果具有一定可信度，可嘗試在更大范圍內推廣該模型。

3.3 土壤侵蝕治理措施

本次調查過程中發現波特瓦爾高原水土保持措施種類相對較少，主要包括水平梯田、排水溝砌石護坡、小型石谷坊、農田地埂、草被護坡等(圖5)。水平梯田、砌石護坡(圖5a)、和小型石谷坊(圖5b)主要應用于波特瓦爾高原北部穆里山區，此處山勢陡峭，當地人為生產安全，多依山勢開挖水平梯田，水平梯田相鄰坡面一般使用砌石護坡，另有排水溝對水平梯田以及上方匯水進行導排，以防止發生溯源侵蝕損毀水平梯田；小型石谷坊一般修建在集水區出口處，而且每個石谷坊都設置有排水管，達到水走沙留的效果。水平梯田在穆里山區較為常見，大部分的房屋修建在水平梯田上，同時在房屋周圍開墾農田(類似中國的石坎梯田)，因削坡開級，水平梯田穩定性很好，另有砌石或者漿砌石護坡，保證了水平梯田不被徑流沖毀，這類水平梯田兼具生產和生態功能，是山區一類良好的水土保持工程措施。在調查范圍內，發現農田的坡度很小，農田少有侵蝕發生，而且大部分農田修建有質量較好的梯形地梗(圖5c)；這種地埂與中國農田埂有較大區別，中國的農田埂一般呈長條形，而波特瓦爾高原的地埂是農田內側邊長約0.3 m，外側約1.5 m，頂寬約0.5 m的梯形地埂，因梯形地埂外邊長度遠大于內邊長度，因此梯形地埂的穩定性比條形地埂更好。地埂防止了農田土壤侵蝕的發生，對農田內部環境有良好的防護作用。鐵路邊坡防護雖然僅采用了簡單的草被護坡(圖5d)，但效果較好，沒有明顯的土壤侵蝕發生，可能是因為修建鐵路的區域地勢平坦，邊坡較短，匯水面積小，在降雨過程中較難生成股流；另外，火車道一般會鋪約0.2 m的石子(道砟)，降雨進入石子之間的大孔隙，不會形成匯流至鐵道邊坡。

表2 各調查單元的平均土壤侵蝕速率

注：土地利用的土壤侵蝕速率是調查單元內部每一種土地利用的所有斑塊(cell)的土壤侵蝕速率統計計算得到的平均值；平均土壤侵蝕速率是由每個調查單元內部的多個圖斑(cell)的土壤侵蝕速率統計計算得到的平均值。

圖5 波特瓦爾高原水土保持措施類型

4 土壤侵蝕的主要影響因素

波特瓦爾高原土壤侵蝕的影響因素分為自然和人為兩方面。高原本身的自然條件較好，平均降雨量由南到北約500～1 900 mm，即整個高原均為半濕潤和濕潤地區，因此人類影響較小的區域天然植被覆蓋較好，土壤侵蝕程度較輕。但同時，北部山區降雨量大，坡度陡，南部區域降雨集中，降雨侵蝕力大，均為土壤侵蝕的發生提供了潛在的有利條件。當人類活動破壞植被后，將產生劇烈的土壤侵蝕，因此人為活動是影響波特瓦爾高原土壤侵蝕的主導因素。波特瓦爾高原的農業人口比重較大，而且整體較為貧困，對自然資源的依賴程度較高，砍伐森林、過度放牧和開墾農田等引起劇烈侵蝕[6]。更重要的一點是，因為貧困，教育普及率很低[3]，所以無論從政府層面，還是普通民眾中，在生產建設活動過程中沒有足夠的政策或者意識防治土壤侵蝕，保護水土資源。如北部山區目前大量零散無序分布的居民地，因毀林修建房屋、道路，對自然植被造成破壞，在降雨和地形的影響下，建設區域產生較為劇烈的侵蝕；中部和南部區域地勢總體平坦，海拔范圍在500～600 m左右，但人口分布眾多(調查單元3—9)，各類型在建項目產生較嚴重的土壤侵蝕(如圖4和圖5所示，包括松散堆積物、新開挖面和強度擾動地面等)。

5 結論與建議

5.1 結 論

本文對巴基斯坦波特瓦爾高原典型單元的土壤侵蝕特征、植被類型與蓋度和土地利用現狀等信息進行了實地調查，結合CSLE模型，對調查單元的土壤侵蝕進行了定量評價，得到以下結論。

土壤侵蝕主要發生在建設用地，包括開挖邊坡、排土邊坡、土路邊坡和磚廠平臺等人為形成的邊坡和裸露平臺在降雨和徑流的作用下產生溝蝕；農耕地因坡度平緩，田埂規格高，土壤侵蝕特征不明顯；部分自然山坡因過度放牧，發生重力侵蝕和溝蝕。因此，生產建設項目的水土保持規劃、監督和治理工作非常重要。

15個調查單元的平均土壤侵蝕模數介于51.40～1 338.91 t/(km2·a)之間。波特瓦爾高原目前僅有的土壤侵蝕實測值約在479～15 000 t/(km2·a)，而CSLE計算所得結果基本在實測值范圍內，說明計算結果具有一定可信度，可嘗試在更大范圍內推廣該模型。

CSLE模型計算結果準確度評估，應以徑流監測小區所得泥沙數據進行比較，而巴基斯坦土壤侵蝕基礎研究薄弱，現階段并沒有規范的徑流監測小區，后期可在巴基斯坦修建小區，對CSLE模型在該國的適用性進行進一步評估。

5.2 建 議

5.2.1 加強土壤侵蝕基礎研究 土壤侵蝕不僅造成農地草地退化，也使得河流泥沙含量增大(印度河的年輸沙量可達4.35×108t)[7]，影響下游生態安全。巴基斯坦地形和地貌類型多變，降雨量空間分布差異很大，侵蝕空間分異也較大[7]，但目前巴基斯坦沒有全國的土壤侵蝕分區成果。另外，本次調查過程中發現巴基斯坦土壤侵蝕監測數據匱乏，僅有的土壤侵蝕監測小區(位于杰格瓦爾PMAS干旱農業大學研究農場(URF))也不規范。小區沒有完全閉合，無法精確收集降雨后全部的徑流和泥沙，難以準確反映從坡面—流域的土壤侵蝕空間特征。巴基斯坦區域土壤侵蝕研究多利用USLE等其它國家開發的模型，在使用時也少有校準，難以保證結果的準確性，未來應立足本國土壤侵蝕特征開發適用的土壤侵蝕模型。應開展土壤侵蝕長期監測工作，明確各地區土壤侵蝕發生發展規律并分析其主導影響因素，為坡面—流域—區域的土壤侵蝕分區和土壤侵蝕治理工作服務。

5.2.2 重視高原溝壑區溯源侵蝕防治 波特瓦爾高原地勢整體平坦，分布有大量耕地。調查單元土地利用結果顯示，農耕地比例最高可達86.89%。同時，該處由于長期水力侵蝕的結果，高原分布有大量寬淺型的溝壑，目前溝壑面積還在擴張。以杰格瓦爾地區PMAS干旱農業大學研究農場(URF)附近溝壑區為例(圖6)，1號和2號溝頭區面積從2001年(圖6a)的67 100 m2和26 046 m2擴張到2018年(圖6b)的154 838 m2和4 4256 m2(在Googleearth勾繪溝頭區，得到其面積)，平均擴張速率分別可達5 161.06 m2/a和1 071.18 m2/a，溯源侵蝕不斷蠶食農田，危害區域土地資源和糧食安全。這種地貌類型可能與中國黃土高原溝壑區的早期狀態相似，其主要的地貌單元包括塬地和溝壑2部分，其中最大的塬地——董志塬(有“隴東糧倉”之稱)。近年來由于溯源侵蝕，大片農田被損毀，也威脅著塬面人們的生命安全，因此已經引起政府部門和專家學者的關注[26]。目前已開展了大量有關溯源侵蝕的研究和固溝保塬工作。巴基斯坦的溝壑雖然多呈現寬淺型，其深度、寬度遠小于黃土高原的溝壑，但擴張速率依然很大。如果對波特瓦爾高原溝壑區的溝頭不進行治理，地貌類型可能會朝著黃土高原溝壑區的地貌演化，大量的耕地將遭到破壞，最終成為“第二個黃土高原”。因此，亟需治理溯源侵蝕，保護塬面耕地，保證糧食安全。

圖6 杰格瓦爾PMAS干旱農業大學研究農場(URF)附近的溝頭變化(4號調查單元附近)

5.2.3 推行生產建設項目的水土保持規劃、監督和治理工作 巴基斯坦是個經濟和人口正處于高速發展的國家，基礎設施建設還有很大的提升空間[4]。建設過程中會產生開挖地面、松散堆積體、棄土棄渣和擾動地面等多種人為侵蝕類型，從而形成嚴重的土壤侵蝕，影響生態環境。

為避免劇烈人為侵蝕發生，首先需進行生產建設項目水土流失規劃。以巴基斯坦北部為例，該區降雨量充沛，森林繁茂，發育有眾多河流，對下游的工農業生產和經濟發展有著重要影響，因此這一水源涵養區域應規劃為生態保護區，禁止人為破壞。但在調查過程中發現，大量的居民地無序散亂分布于林間(圖7a)，部分開發區域已產生劇烈的土壤侵蝕(圖7b)；因此，大量的開發項目需進行統一的規劃和審批。其次，對建設過程的水土流失進行及時監督和治理，以保護巴基斯坦的水源涵養地區。

如何對各類人為侵蝕進行有效監督和治理？目前巴基斯坦沒有對生產建設項目土壤侵蝕進行監督，防治措施薄弱，但相比較巴基斯坦大量的生產建設項目，現有防治面積和防治措施的種類還遠遠不夠，亟需針對各類建設項目有針對性地進行監管和布設土壤侵蝕防治措施。這方面可以借鑒中國的經驗，在監督管理方面，中國已從以往的文件、報告、現場抽查等方式進入了“3S”監督檢查時代，具有高效、快速、客觀和精確等特點[27]；生產建設項目土壤侵蝕防治措施方面，中國對臨時松散堆積體、公路邊坡和礦區排土場邊坡等不同建設項目產生的易侵蝕地表，根據各地氣候特點，采取稻草簾子[28]、無紡布、遮陽網覆蓋[29]和“沙柳方格+沙打旺”[30]等土壤侵蝕防治措施，并取得了良好的效果。

圖7 穆里地區居民地分布概況

致謝：感謝巴基斯坦PMAS-Arid Agriculture University 對本次合作交流和野外考察的大力支持，特別感謝Shazada Sohail Ijaz教授團隊工作人員在野外調查過程中給予的無私幫助，也感謝巴基斯坦政府對調查人員生命安全的保障工作；感謝西北大學楊勤科教授團隊師生為本研究提供的降雨數據以及對15個調查單元的土地利用圖像解譯以及土壤侵蝕模數計算工作。