湄公河典型區域旱雨季土地利用及土壤侵蝕差異研究

趙娜 張道熙 韓東楓 陳小娟 邱娟 金瑤 龔昱田

摘要：以湄公河干流老撾沙灣拿吉至柬埔寨桔井段重要支流南蒙河、色丹河、色貢河、色桑河、斯雷博河及其左右岸緩沖區50 km范圍為研究區域，基于Landsat8遙感和河網數據，對比植被覆蓋度以及空間分布特征，分析區域內旱季（1-2月）、雨季（6-9月）土地利用及土壤侵蝕變化。結果表明，季節性降雨導致土地利用變化明顯，旱、雨季土地利用變化主要表現為河道沿岸水體動態變化；旱季水體面積較雨季凈減少43 259 hm2，占總變化面積的93.8%；雨季至旱季耕地轉化占總變化面積38.00% （17 514.44 hm2），其次為灌木地（22.98%，10 590.54 hm2）和裸地（21.77%，10 033.32 hm2）。土地利用類型轉化導致土壤侵蝕變化，旱、雨季土壤極強烈侵蝕和劇烈侵蝕區域變化不大；旱季以輕度侵蝕、中度侵蝕增加為主，凈增加面積分別為4 107 hm2和5 864 hm2，分別占總變化面積的36.9%和52.7%；增加面積主要由無侵蝕類型發展而來，體現在旱季沿岸水體面積的大幅降低。建議在水土保持過程中，改傳統稻田種植為耕作、輪作和培肥一體化的水保措施，增強基礎設施建設，加強河岸帶生態護坡建設與管理，降低自然洪水或人為工程對自然林地和植被的破壞。

關鍵詞：土地利用；土壤侵蝕；旱雨季；湄公河

中圖分類號：X144? ? ? ? 文獻標志碼：A? ? ? ? 文章編號：1674-3075（2024）02-0047-07

瀾滄江-湄公河是世界第六、亞洲第三大河流，發源于中國，在云南出境后，流經緬甸、老撾、泰國、柬埔寨和越南，最終匯入南海（長江水利委員會國際合作與科技局，2016）。近年來，在全球氣候變化背景下，瀾湄流域（特別是下游區域）的洪澇、干旱等自然災害呈突發、多發、并發態勢，已成為制約區域經濟社會可持續發展的重要因素之一。

土壤侵蝕是指地表土壤及母質在自然或人為因素影響下被破壞、剝蝕、搬運和沉積的過程（王占禮，2000；肖繼兵等，2017；史東梅等，2019），目前已成為全球主要的生態環境問題；其危害包括土壤資源退化、土壤肥力下降、淤積抬高河床、加劇洪澇災害、淤塞水庫湖泊、影響開發利用等（劉爭光，2021；Borrelli et al， 2021）。上游瀾滄江區域已有較為詳細的研究（甘淑，2001；Barton et al，2004；姜昀等，2006；姚華榮等，2006；姚華榮和崔保山，2006；李立科，2011）；下游湄公河區域為熱帶季風氣候，雨熱同期，水網密布，土壤肥沃，年累計降水量最低為440 mm，最高為3 500 mm。70%的柬埔寨人口和53%的老撾人口生活在湄公河兩側15 km范圍的河道走廊內，高度依賴于水資源謀生（MRC， 2011）。

隨著湄公河流域內與水資源相關的航運、灌溉、水電和漁業等工農業快速發展，流域水生態系統面臨著諸多挑戰（MRC， 2015）；同時，建設用地、耕地的大幅度增加，也多來源于林地和草地等土地類型的轉化或退化，生態系統服務功能也出現不同程度的衰退（歐陽志云等，1999；畢曉麗和葛劍平，2004；趙軍和楊凱，2007；董哲仁，2008；伍星等，2009；陸穎等，2013）。土壤侵蝕已經是湄公河流域主要環境問題，但相關的研究比較滯后（甘淑等，2001；吳芳等，2019）。因此，探究湄公河重點區域的旱季、雨季土地利用及土壤侵蝕變化具有重要意義。

1? ?材料與方法

1.1? ?區域概況

研究區域位于中南半島（105°10′54″～107°43′43″E，11°48′3″～15°45′18″N），總面積7.07萬km2， 以湄公河干流老撾沙灣拿吉至柬埔寨桔井段重要支流南蒙河（Nam Mun River）、色丹河（Se Done River）、色貢河（Se Kong River）、色桑河（Se San River）、斯雷博河（Sre Pok River）及其左右岸緩沖區50 km的范圍為邊界，主要涉及老撾南部的占巴色省（Chambasse）和阿速坡省（Attapeu）、柬埔寨北部的上丁省（Stung Treng）和桔井省（Kratie）以及泰國東北部的烏汶府（Ubon Ratchatani）（圖1）。區域內有著名的3S小流域、四千美島、西潘墩濕地、湄公河洪泛平原、伊諾瓦底江豚保護區等，生態功能非常重要；常年炎熱，具有很強的季節性降雨特點，年際間年總降雨量的變化幅度一般在±15%，干支流年際間流量的季節性變化較大。

1.2? ?研究方法

基于拼接的30 m分辨率Landsat 8遙感數據和河網數據，勾繪出湄公河干流老撾沙灣拿吉至柬埔寨桔井段主河道及重要支流南蒙河、色丹河、色貢河、色桑河、斯雷博河，并建立距離為50 km的緩沖區。通過分析研究區域及子流域的植被覆蓋度、土壤侵蝕程度，比較2010年旱季（1-2月）和雨季（6-9月）的水域分布、土地利用狀況變化及沿岸帶侵蝕等要素的差異。研究區域土地利用數據來自于湄公河委員會（MRC） http：//portal.mrcmekong.org。

1.3? ?數據分析

1.3.1? ?子流域劃分? ?利用ArcGIS軟件對研究區域數據進行處理，經過填洼、生成流向、水流累積量等分析，生成1185個子流域。使用Raster to Polygon將子流域由柵格轉為矢量數據，如圖2所示。

1.3.2? ?坡度數據? ?使用Slope計算獲得坡度數據，通過Zonal Statistics as Table獲得各子流域范圍內坡度平均值（圖3）。

1.3.3? ?植被覆蓋度（VFC）? ?基于ENVI5.3軟件平臺，對Landsat 8衛星數據進行輻射定標、幾何校正和大氣校正，計算歸一化植被指數（NDVI）。利用研究區域矢量數據生成的ROI建立一個掩膜文件，統計NDVI值的累積概率分布，分別取累積概率為5%和95%的NDI值作為NDVImin和NDVImax。旱季和雨季的歸一化植被指數（NDVI）計算結果如下：NDVImaxdry = 0.560253，NDVImaxrainy = 0.522455，NDVImindry = 0.012260，NDVIminrainy = 0.147679。

將整個地區分為3部分。NDVI<0.012260，VFC取值為0；NDVI>0.560253，VFC取值為1；兩者之間的像元使用公式計算，得到單波段的植被覆蓋度圖像文件。使用Band Math工具計算，公式如下： [FVC =]

式中：FVC表示植被覆蓋度，NDVIsoil為NDVI值累計分布頻率5%的值，NDVIveg為NDVI值累計分布頻率95%的值。

1.3.4? ?土地利用? ?研究區域內具有代表性的土地利用類型為農業用地（耕地）、林地、城鎮用地、濕地等；其中，城鎮用地由于面積有限、分布范圍較小，在子流域尺度上存在大量零值，而濕地部分歸于河網密度指標。因此，研究中土地利用指標為農業用地和林地。使用ArcGIS 10.7中的Tabulate Area工具計算旱季、雨季各土地利用類型在子流域的面積，然后計算出各土地利用類型在子流域中的占比（圖4）。

1.3.5? ?土壤侵蝕? ?采用國內土壤侵蝕分類分級標準（中華人民共和國水利行業標準SL 190-2007），共分為6個侵蝕等級；土壤侵蝕基于植被覆蓋、坡度和土地利用（旱、雨季）數據計算得到。

2? ?結果與分析

2.1? ?旱季和雨季土地利用變化

研究區域旱、雨季土地利用變化如圖5所示。旱季水體面積凈減少約43 259 hm2，占總變化面積的93.8%；耕地凈增加約16 863 hm2，占36.6%，其次為裸地（21.7%）和灌木地（21.2%）（表1）。耕地主要位于研究區域南部湄公河主河道沿岸，灌木主要位于湄公河主河道中北部及支流色桑河、斯雷博河下游沿岸，裸地主要位于湄公河主河道中南部。這些區域涉及著名的3S小流域、老撾四千美島和西潘墩濕地以及柬埔寨的湄公河洪泛平原，是湄公河流域典型的洪水高發區；湄公河流域明顯的季節性降雨特點，導致這些區域土地利用類型的變化也非常明顯。

2.2? ?旱季和雨季土壤侵蝕變化

研究區域旱、雨季土壤侵蝕情況如圖6所示。雨季凈增加約10 543 hm2的無侵蝕區域，占總變化面積的94.8%；旱季輕度、中度、強烈、極強烈和劇烈侵蝕面積均有所增加，以輕度和中度侵蝕增加為主，凈增加面積分別為4 107 hm2和5 864 hm2，分別占總變化面積的36.9%和52.7%（表2）。研究表明，土壤侵蝕變化明顯的區域位于南部湄公河主河道和3S小流域，這些區域是著名的水稻主產區，也是人口稠密的地方。由于氣候變化、基礎設施發展和人口增長的影響，湄公河流域土壤侵蝕有加劇趨勢（MRC， 2015）。

3? ?討論

3.1? ?季節性降雨導致土地利用變化明顯

湄公河流域為熱帶季風氣候，常年炎熱，季節性潮濕。冬季和夏季氣流逆轉，熱帶季風氣候具有很強的季節性降雨特點。湄公河干流桔井監測站和色貢河、色桑河和斯雷博河主要支流監測站數據顯示，8月和9月平均流量都在30 000 m3/s以上，干流巴色站、支流色丹河站和南蒙河在20 000 m3/s以上。上述監測站在旱季（2-4月）月均流量都在5 000 m3/s以下（長江水利委員會國際合作與科技局，2016）；同時，柬埔寨70%的農業人口和老撾53%的農業人口生活在河岸兩側15 km寬的河道走廊內，且主要以種植和捕魚為生（MRC， 2017）；其中水稻（包括雨季稻、浮稻、洪水漫灌稻、旱季稻等）是流域內大部分居民的主要種植作物（MRC， 2019）。在明顯的季節性降雨條件下，雨季和旱季的水稻平均產量存在明顯差異（FAO-RDES 糧農組織-亞洲數據交換系統）。研究區域旱、雨季土地利用變化主要表現為河道沿岸水體動態，旱季水體轉為耕地的占比最大，且耕地主要位于研究區域南部湄公河主河道沿岸。

3.2? ?土地利用類型轉化導致土壤侵蝕變化

土壤侵蝕是湄公河流域主要的環境問題（吳芳等，2019）。總體上看，旱、雨季土壤極強烈侵蝕和劇烈侵蝕的區域變化不大，其他類型變化顯著。旱季中度和輕度侵蝕面積相對于雨季有較大幅度的增加，雨季無侵蝕面積相對于旱季增加顯著。

湄公河流域土壤侵蝕的主導因素是土地利用類型和坡度（吳芳等，2019）；本研究結果也證實了這一點。區域內旱、雨季土地利用類型在一定范圍內發生顯著變化，繼而導致了土壤侵蝕方式和強度的變化。例如，耕地向林地、草地、水域的轉化，將導致侵蝕強度顯著降低；林地向耕地、草地轉化，或者水域向其他土地利用類型的轉化，均使侵蝕強度顯著增加。

吳芳等（2019）基于土壤侵蝕空間特征，對79.5萬km2 湄公河流域的33個子流域進行了4個優先治理等級分析。對照其結果，沙灣拿吉-桔井片區屬于微度和輕度侵蝕，治理等級為第4級和第3級。本研究中，研究區域旱季和雨季無侵蝕-輕度侵蝕面積占比分別為91.20%和94.15%，強烈-極強-劇烈侵蝕面積占比為0.45%和0.53%。相較于湄公河流域其他區域，沙灣拿吉-桔井片區土壤侵蝕強度不高，與該區域裸地占比較少、植被覆蓋度高有關。

3.3? ?加強湄公河流域水保措施和基礎設施建設

沙灣拿吉-桔井片區是湄公河流域典型區域，有著名的3S小流域、四千美島、西潘墩濕地和湄公河洪泛平原等，其生態功能非常重要。區域內人口密集，種植業和漁業發達，同時還具有貿易頻繁、干支流水電開發較快、經濟發展迅猛等特點。因此，在其水土保持過程中，改傳統稻田種植為耕作、輪作和培肥一體化的水保措施，同時增強基礎設施建設，加強河岸帶生態護坡建設與管理，降低自然洪水或人為工程對自然林地和植被的破壞。

參考文獻

畢曉麗，葛劍平，2004. 基于IGBP土地覆蓋類型的中國陸地生態系統服務功能價值評估[J]. 山地學報， 22（1）：48-53.

董哲仁，2008. 河流生態系統結構功能模型研究[J]. 水生態學雜志， 29（1）：1-7.

甘淑，何大明，袁希平，2001. NDVI在瀾滄江流域山區林地覆蓋遙感監測中的應用研究[J]. 林業科學， 37（4）：134-136.

姜昀，高吉喜，歐曉昆，2006. 瀾滄江流域云南段土地利用格局變化及環境影響分析[J]. 環境科學研究， 19（3）：46-51.

李立科，2011. 云南瀾滄江流域植被時空動態及對區域氣候變化的響應[D]. 上海：華東師范大學.

劉爭光，2021. 土壤侵蝕研究進展[J]. 農業與技術雜志， 4（1）：92-96.

陸穎，潘鋒，何大明，2013. 基于分布式水文模型的藤條江河流土壤侵蝕研究[J]. 水生態學雜志， 34（1）：1-6.

歐陽志云，王如松，趙景柱，1999. 生態系統服務功能及其生態經濟價值評價[J]. 應用生態學報， 10（5）：635-640.

史東梅，金慧芳，蔣光毅，2019. 土壤侵蝕對坡耕地耕層質量退化作用及其評價趨勢展望[J]. 農業工程學報， 35（18）：118-126.

王占禮，2000. 中國土壤侵蝕影響因素及其危害分析[J]. 農業工程學報， 16（4）：32-36.

吳芳，朱源，許丁雪，等，2019. 湄公河流域土壤侵蝕空間特征及其優先治理區確定[J]. 生態學報， 39 （13）：4761-4772.

伍星，沈珍瑤，劉瑞民，等，2009. 土地利用變化對長江上游生態系統服務價值的影響[J]. 農業工程學報， 25（8）：236-241.

肖繼兵，孫占祥，劉志，等， 2017. 降雨侵蝕因子和植被類型及覆蓋度對坡耕地土壤侵蝕的影響[J]. 農業工程學報， 33（22）：159-166.

姚華榮，崔保山，2006. 瀾滄江流域云南段土地利用及其變化對土壤侵蝕的影響[J]. 環境科學學報， 26（8）：1362-1371.

姚華榮，楊志峰，崔保山，2006. GIS支持下的瀾滄江流域云南段土壤侵蝕空間分析[J]. 地理研究， 25（3）：421-429.

長江水利委員會國際合作與科技局，2016. 湄公河開發利用與管理[M]. 武漢：長江出版社.

趙軍，楊凱，2007 .生態系統服務價值評估研究進展[J]. 生態學報， 27（1）：346-356.

中華人民共和國水利部， 2008. 土壤侵蝕分類分級標準：SL190-2007[S/OL]. 北京：中國水利水電出版社.[2008-01-04]https：//wenku.so.com/d/0ab7f214901fb46b98fac34ec00a28bb

Barton A P， Fullen M A， Mitchell D J， et al， 2004. Effects of soil conservation measures on erosion rates and crop? productivity on subtropical ultisols in Yunnan Province[J]. China Agriculture Ecosystems and Environment， 104：343-357.

Borrelli P， Alewell C， Alvarez P， et al， 2021. Soil erosion modelling： A global review and statistical analysis[J]. The Science of the Total Environment， DOI：10.31223/X5GS3T

MRC， 2011. Annual Mekong Flood Report 2010[R]. Mekong River Commission， Vientiane， Lao PDR.

MRC， 2015. Annual Mekong Flood Report 2014[R]. Mekong River Commission， Vientiane， Lao PDR.

MRC， 2017. Mekong basin-wide fisheries management and development strategy 2018－2022[R]. The Mekong River Commission，Vientiane， Lao PDR.

MRC， 2019. Final Report-Agricultural land use Monitoring[R]. The Mekong River Commission， Vientiane， Lao PDR.

（責任編輯? ?萬月華）

收稿日期：2021-09-03? ? ? 修回日期：2022-06-20

基金項目：中國-東盟國家防洪抗旱應急管理合作項目。

作者簡介：趙娜，1978年生，女，博士，研究員，主要從事水生生物學研究。E-mail：chinesezhaona@163.com