秦巴山區退耕還林(草)典型區域土地利用和土壤保持服務變化研究
——以文縣為例

張 強，陸 蔭，楊 青

(甘肅省環境監測中心站，甘肅 蘭州 730020)

1 引言

土壤侵蝕是指土壤或其他地面組成物質在自然營力或人類活動的共同作用下被剝蝕、破壞、分離、搬運和沉積的過程，是威脅生態安全的重要問題之一，已受到國內外學者的廣泛關注[1，2]。健康的生態系統需要一定的侵蝕和沉積，但土壤侵蝕過度會破壞耕地資源，造成一系列自然災害，嚴重影響人類的生產生活活動和經濟社會的可持續發展[3]。土壤保持是生態系統提供的重要服務功能之一，主要包括侵蝕控制和沉積物保持。土壤保持與地形地貌、土壤性質、植被蓋度、降雨強度和農業活動、工程建設等人為因素密切相關[4]。退耕還林(草)等生態工程對土地利用結構變化影響巨大，而土地利用結構變化與土壤保持服務存在強烈耦合關系[5，6]。隴南市文縣處于秦巴山區的最西端，西秦嶺侵蝕構造中高山地區，為甘肅省地質災害最發育，危害最嚴重的縣區之一[7]，自1999年實施退耕還林(草)等生態工程后，文縣的土地利用結構發生明顯變化，土壤保持功能也隨之改變，但目前在該區域研究生態工程的實施對土地利用和土壤保持的影響相對較少。

通用土壤流失方程USLE(Universal Soil Loss Equation)以潛在土壤侵蝕量和實際侵蝕量之差來計算土壤保持量，是目前應用最廣泛的土壤侵蝕模型[8]。但該模型未考慮地塊自身攔截上游沉積物的能力。InVEST模型(The Integrate Valuation of Ecosystem Services and Tradeoffs Tool)是由美國斯坦福大學環境森林研究所、世界自然基金(WWF)和大自然保護協會(TNC)于2007年共同研發的一種生態系統服務與權衡綜合評估模型[9]。該模型的土壤保持模塊考慮了地塊本身攔截上游沉積物的能力，不僅能更加準確地計算地塊的潛在土壤侵蝕量、實際土壤侵蝕量、土壤保持量，而且還能夠結合沉積物清除成本、水庫設計和貼現率來計算沉積物保留價值，為水庫管理者提供決策參考[10]。

鑒于此，本研究以文縣為研究區，采用InVEST模型和ArcGIS平臺，選擇2000年和2015年兩個研究時段，評估區域退耕還林(草)等生態工程實施后土地利用的變化和土壤保持服務的變化，以期為該區域的生態環境治理和土壤侵蝕治理方案制定提供科學依據。

2 研究區概述

文縣隸屬甘肅省隴南市，處于甘肅省與四川省的交界處，地理位置介于E104°16′16″～105°27′29″，N32°35′43″～33°20′36″，總面積約5000 km2[11]。區域地處西秦嶺侵蝕構造中高山地區，地勢西北高東南低，地表起伏大，海拔介于550～4187 m，屬亞熱帶向暖溫帶的過渡地帶，多年平均降水量介于400～1000 mm，地震活動頻繁，地質災害頻發[12]。2018年文縣的三次產業結構比為24.33%∶24.31%∶51.36%，其中，傳統農業發展具有悠久的歷史，但耕地多為陡坡旱地，土壤肥力低，形成粗放落后的農業生產方式；礦產資源豐富，已探明的金屬非金屬礦藏30余種。且隨著人口的增多，新開墾耕地面積增多，水利工程興建，炸山開路，修路建房等一系列人類活動，造成大面積的水土流失，加劇自然災害的發生頻數[7]。

3 研究方法

3.1 土地利用變化分析

利用ArcGIS10.6，開展文縣2000年和2015年兩期土地利用的空間分布及其轉移矩陣研究，分析退耕還林(草)工程實施后的土地利用變化情況，進一步分析了土地利用變化對土壤保持服務的影響。

3.2 土壤保持服務核算

本研究運用InVEST模型，開展文縣2000年和2015年兩期土壤保持服務的變化分析，土壤保持服務包括土壤保持服務的物質量和價值量。InVEST模塊中土壤保持物質量由兩部分組成，一部分是地塊自身攔截的沉積物保留量，用上坡來沙量與泥沙截留率乘積表示。另一部分是因植被覆蓋和水土保持措施而減少的土壤侵蝕量，用潛在土壤侵蝕量與實際土壤侵蝕量的差值表示；土壤保持的價值量涵蓋土壤保持的肥力價值和減少泥沙淤積價值，土壤肥力價值為養分價值和有機物價值之和，減少泥沙淤積價值采用影子工程法測算[13]。

4 數據來源

研究所需的空間數據主要包括土地利用/覆被數據，氣象數據，土壤質地數據，DEM數據，生物物理表、匯水積累量閾值、BorselliK參數、BorselliICo參數和SDR 最大值等InVEST模型所需參數數據。各數據及模型主要參數的獲取路徑如表1所示。

表1 數據來源及其獲取路徑

5 文縣土地利用變化與土壤保持服務變化

5.1 文縣土地利用動態變化分析

根據文縣2000～2015年的土地利用類型面積及百分比、土地利用轉移矩陣(表2、表3)可得，文縣的林地、草地和耕地占全區域土地面積的主要部分，其占比在97%以上。2000～2015年，文縣耕地的轉出面積大于轉入面積，為唯一轉出的土地利用類型，林地、草地、水域、城鄉、工礦、居民用地和未利用土地的轉入面積大于轉出面積，均為轉入土地利用類型。由于退耕還林(草)政策的實施，大部分耕地轉為林地和草地，占耕地轉出面積的70.73%。其中，林地占比最大，多年平均土地面積占土地總面積比例超過46%，2000～2015年林地面積增加11.10 km2，增幅為0.48%；其次占比較大的為草地，其多年平均土地面積占土地總面積比例超過33%，2000～2015年草地面積增加35.33 km2，增幅為2.10%；耕地的多年平均土地面積占土地總面積的比例為17.44%，2000～2015年耕地面積下降65.65 km2，呈明顯的下降趨勢，其轉出比例達7.26%；水域和城鄉、工礦、居民用地、未利用土地的面積占土地總面積的比例較小，不足1%，各土地利用面積均呈明顯的上升趨勢，因其基數較小，其增幅較大，分別為19.71%、39.24%、14.27%。其中，城鄉、工礦、居民用地是所有土地利用類型中增幅最大的，根據文縣統計年鑒，2000～2015年區域總人口約增加2.26萬人，人口的快速增長導致城鄉、工礦、居民用地的需求量增加，城鄉、工礦、居民用地面積的增加可能與人口的增加有關。

表2 2000～2015年文縣土地利用類型百分比 %

表3 2000～2015年間文縣土地利用類型轉移矩陣 km2

從文縣的土地利用空間分布(圖1)來看，2000年和2015年文縣的兩期土地利用空間分布較為一致，各類土地利用呈明顯的空間差異性。其中，林地集中分布在文縣的北部和南部部分區域，草地鑲嵌分布于林地及其邊緣區域；耕地集中分布在區域中部的河谷區域，居民點也分布在河谷區域，且沿河流分布，未利用土地零星分布在區域北部。

圖1 研究區土地利用

林地、草地和耕地是文縣的主要土地利用類型，且為退耕還林還草的主要變化類型，本研究分析了退耕還林還草的空間分布、及其雙向流動的空間分布(圖2)。耕地的轉出面積大于轉入面積，其中，退耕還林還草的斑塊面積較大，在全域均有分布，且轉出耕地多為山地旱地，山地旱地轉為草地的面積占退耕還林還草面積的74.05%，山地旱地轉為林地的面積占退耕還林還草面積的23.26%；林地和草地轉為耕地的區域在北部和中部均有分布，其轉出面積分別為6.66 km2、38.50 km2。林地和草地也存在雙向流動，林地轉為草地的轉出面積大于草地轉為林地的轉出面積，林地轉為草地的區域集中分布在區域西北部，其轉出面積為33.56 km2，草地轉為林地的區域集中分布在區域中東部，其轉出面積為28.68 km2。

圖2 2000～2015年研究區土地資源變化表4 文縣各地類土壤保持量

5.2 文縣土壤保持生態系統服務變化

根據InVEST模型測算可得，2000年、2015年文縣的土壤保持量總量分別為4.87×108t、4.96×108t，區域土壤保持總量增幅為9.38×106t。

如表4所示，各地類中林地、草地和耕地對文縣土壤保持量的貢獻率較高，2000年和2015年兩期的平均貢獻率分別為56.23%、29.86%、11.88%，水域、城鄉、工礦、居民用地和未利用土地的貢獻率較小，其兩期均小于1.05%。其中，林地、草地、城鄉、工礦、居民用地和水域的土壤保持總量呈上升趨勢，林地的土壤保持總量增幅最大，草地次之，2000～2015年的林地、草地的土壤保持總量增幅分別為10.92×106t、2.35×106t，雖然林地面積增幅小于草地，但林地的平均土壤保持量最大，其多年平均土壤保持量為11.85×104t/km2；耕地和未利用土地的土壤保持總量呈下降趨勢，其土壤保持總量降幅分別為4.44×106t、0.17×106t，由于未利用土地面積較小，且耕地的平均土壤保持量較小，其多年平均土壤保持量僅為6.69×104t/km2，導致土壤保持的降幅明顯小于增幅。

一級地類二級地類 平均土壤保持量104 t/km22000年土壤保持總量104 t2015年土壤保持總量104 t耕地水田1.9731.23 30.80 旱地6.786025.99 5582.35 林地有林地12.5120185.70 20103.50 灌木林11.096882.40 6904.05 疏林地11.026.04 1002.51 其他林地9.478.14 162.94 草地高覆蓋度草地10.198319.48 8419.88 中覆蓋度草地7.276203.45 6304.99 低覆蓋度草地2.927.44 60.83 水域河渠11.83342.98 315.58 湖泊9.250.78 0.89 水庫坑塘6.956.24 55.89 灘地7.442.45 10.89 城鄉、城鎮用地1.541.26 2.36 工礦、農村居民點5.77105.87 138.39 居民用地其他建設用地6.55—7.64 未利用土地裸土地14.62—113.44 裸巖石礫地8.8743.66 253.10 其他13.22462.78 123.39

考慮保持土壤肥力價值和減少泥沙淤積價值，采用市場價格法計算可知，文縣2000年和2015年兩期的土壤保持總價值分別為316.26億元、322.36億元，區域土壤保持總價值量增幅為6.09億元。

如表5所示，文縣的各地類中林地、草地和耕地的土壤保持價值量較高，2015年各地類的土壤保持價值總量分別為183.12億元、96.11億元、36.49億元，水域、城鄉、工礦、居民用地和未利用土地的土壤保持價值量較小，2015年各地類的土壤保持價值量分別為2.49億元、0.96億元、3.18億元。與土壤保持總價值量的變化趨勢一致，林地、草地、城鄉、工礦、居民用地和水域的土壤保持價值量呈上升趨勢，耕地和未利用土地的土壤保持價值量呈下降趨勢。

表5 文縣各地類土壤保持價值 萬元

6 結語

為評估退耕還林(草)工程對土壤保持服務的影響，本研究以秦巴山區隴南市文縣為例，采用InVEST模型和ArcGIS平臺，對文縣2000年和2015年兩期的土地利用和土壤保持服務變化進了行分析，結果如下。

(1)文縣的土地利用類型主要由林地、草地和耕地組成，其占比占97%以上。2000～2015年期間，文縣耕地的轉出面積大于轉入面積，為唯一轉出的土地利用類型，林地、草地、水域、城鄉、工礦、居民用地和未利用土地的轉入面積大于轉出面積，均為轉入土地利用類型，且城鄉、工礦、居民用地的增幅最大。

(2)耕地、林地、草地之間存在雙向流轉，其中，退耕還林還草區域多為山地旱地，山地旱地轉為草地的面積占退耕還林還草面積的74.05%，山地旱地轉為林地的面積占退耕還林還草面積的23.26%。

(3)各地類的土壤保持服務能力存在明顯差異，總體呈現出：林地>草地>耕地>未利用地>水域>城鄉、工礦、居民用地的特征。

(4)2000～2015年間的退耕還林(草)工程的實施有效提高了文縣的土壤保持能力，區域的土壤保持物質總量增加了937.53×104t，總價值增加了6.09億元。