土地利用變化對太行山區水源涵養功能的影響研究

摘要：太行山是京津冀城鎮群和華北平原的生態屏障和重要的水源涵養地，為探究土地利用變化對水源涵養功能的影響，基于土壤、土地利用、DEM、氣象等數據，借助InVEST模型定量評估太行山區的水源涵養功能，并分析土地利用變化對水源涵養功能的影響。結果表明：（1）2000—2020年太行山區多年平均水源涵養總量為53.67×108" m3，呈現出W型的變化趨勢；水源涵養深度均呈現出由南向北、由東向西逐漸降低的空間分布特征。（2）不同地類的水源涵養功能強弱差異明顯，依次為草地＞林地＞未利用地＞耕地＞建設用地＞水域。（3）林地和耕地與水源涵養功能呈現顯著正相關，而草地、建設用地和未利用地與水源涵養功能呈現顯著負相關。

關鍵詞：InVEST模型；水源涵養；土地利用；時空變化；太行山

中圖分類號：X171.1；K930文獻標志碼：A文章編號：16735072（2025）02016208

Study on the Impact of Land Use Changeson Water Conservation Function in Taihang Mountains

LI Qinglan，WANG Chengwu，WANG Zhoufeng，XIE Liang

（School of Geoscience and Technology，Southwest Petroleum University，Chengdu Sichuan 610500）

Abstract：The Taihang Mountains are an ecological barrier and an important water conservation area for the BeijingTianjinHebei urban agglomeration and the North China Plain.To explore the impact of land use changes on the water conservation function，based on soil，land use，DEM，and meteorological data，we quantitatively assessed the water conservation function of the Taihang Mountains with the help of the InVEST model，and explored the impact of land use changes on the water conservation function.The results show that：（1） From 2000 to 2020，the total average annual water conservation in the Taihang Mountains was 5367×108 m3，showing a Wshaped trend.The depth of the water conservation all presents the spatial distribution of gradually decreasing from the south to the north and from the east to the west.（2） The water conservation functions of different land use types differ significantly，with the following order of magnitude：grassland gt; forest land gt; unused land gt; cropland gt; construction land gt; water.（3） Forest land and cropland show significant positive correlation with the water conservation function，while grassland，construction land and unused land show significant negative correlation with the water conservation function.

Keywords：InVEST model; water conservation; land use; spatial and temporal variation; Taihang Mountains

水源涵養功能是重要的生態服務功能，在生態系統與區域水循環中發揮著水源供給、調節徑流、水土保持等關鍵作用。土地利用變化是自然環境和人類活動綜合作用的結果，直接影響著生態系統的演化過程，進而驅使生態系統的結構和功能發生改變。土地利用變化主要通過改變下墊面條件間接改變區域蒸發、入滲、產流等水文過程，從而使區域水源涵養功能發生變化［1］。水源涵養功能與土地利用變化聯系密切，相互影響，并共同影響著一個地區的社會經濟和生態系統。因此，研究水源涵養功能對土地利用變化的響應對區域土地資源和水資源的管理具有重要意義。

近年來，空間分析技術在生態水文過程研究中得到了廣泛應用。國內外學者運用VIC（Variable Infiltration Capacity）［2］、SCSCN（Soil Conservation Service Curve Number）［3］、SWAT（Soil and Water Assessment Tool）［4］和InVEST（Integrated Valuation of Ecosystem Services and Tradeoffs）等模型來探討區域水源的涵養功能。其中，InVEST模型具有廣泛的適用性和強大的動態性，能夠基于GIS平臺實現生態服務功能的可視化與空間分析。在水源涵養功能的模擬和評估分析中，InVEST模型具有一定的優勢，被國內外學者廣泛應用，如利用InVEST模型評估加納普拉河［5］、伊朗半干旱森林［6］和埃塞俄比亞Winike［7］等流域的供水生態系統服務功能，驗證了模型在流域尺度的準確性和適用性；基于InVEST模型模擬流域［811］、山區［12］、市縣［1314］等不同尺度的水源涵養功能及其時空分布，并分析和探討了水源涵養對土地利用與氣候變化的響應。以上研究表明，InVEST模型在不同地理環境的研究區中均取得了較好的應用效果，是量化評估水源涵養的有效工具。隨著人類活動對生態環境的干擾日益加劇，對土地利用變化導致的生態系統價值響應的研究變得尤為重要。而目前的研究多聚焦于水源涵養功能的時空變化及其影響因素［67，910，13］，對土地利用與水源涵養之間耦合關系的研究相對較少［1415］，需要進一步探討。值得一提的是，InVEST模型可以根據山區生態系統的特點和實際情況進行參數本地化調整，將生態系統服務功能轉化為具體的數值，從而定量地評估水源涵養的效益。這有助于更深入了解不同土地利用類型的水源涵養功能，并更好地進行土地利用與水源涵養之間耦合關系的系統研究。

中國的山地面積約占陸地國土空間的64.89%［16］，是自然生態系統的重要組成部分，是國家可持續發展的重要生態基底。山地生態系統在水源涵養、土壤保持、氣候調節、生物多樣性維持、物質生產資料供給等方面具有突出的作用［17］，山地生態系統提供的功能不僅服務于山區本身的發展，更在國家生態安全格局中起著骨架作用。太行山區是華北平原和京津冀地區重要的生態屏障和水源涵養功能區，但受自然環境變化以及人類活動的影響，區域水土流失嚴重。目前，有學者對太行山區土地利用景觀格局時空變化［18］、氣候與植物覆蓋度的耦合關系［19］、土壤侵蝕時空分異［20］、水資源供需關系［21］等進行了研究，但缺少土地利用變化與水源涵養功能之間耦合關系的研究。為此，本文運用InVEST模型，從水源涵養功能變化的角度分析評估土地利用變化對其的影響，以期為太行山區土地利用規劃管理、水資源合理配置提供參考。

1研究區概況

太行山區（110°19′—116°40′E，34°35′—40°46′N）是中國東部地區的重要山脈和地理分界線，位于山西省與華北平原之間（圖1），縱跨北京、河北、山西、河南4個省市［22］，面積約為13.48萬km2，平均海拔約為861 m。該區屬于暖溫帶半濕潤大陸性季風氣候，氣溫從南到北逐漸降低，年平均氣溫約為10 ℃；南部地區年均降水量為550～650 mm，北部地區為400～550 mm［21］，降水主要集中在7—9月，這段時間暴雨較多從而引起嚴重的水土流失；主要土地利用類型為耕地、林地和草地，三者占總面積的90%左右；主要植被類型為溫帶落葉闊葉林；主要土壤類型是山地褐土；太行山區是黃河流域和海河流域的重要水源區，孕育了眾多河流，包括漳河、沁河、丹河、子房河、峪河、淇河等大型河流。

2研究方法

2.1數據來源本文所采用的數據主要為土地利用數據、氣象數據、土壤數據、生物物理系數表、Zhang系數、流速系數和地形指數。所有柵格數據空間分辨率均統一重采樣為1 km，并投影為WGS_1984_UTM_Zone_50坐標系。

土地利用數據采用中國科學院資源環境科學數據中心（https：//www.resdc.cn/）提供的2000、2005、2010、2015、2020年5期土地利用遙感監測數據，空間分辨率為1 km。根據研究需求，使用ArcGIS軟件將原始土地利用遙感監測數據的25個二級地類重分類為耕地、林地、草地、水域、建設用地、未利用地6個一級地類。

氣象數據包括降水和潛在蒸散量，均來源于國家地球系統科學數據中心（http：//www.geodata.cn）提供的逐月數據集。使用ArcGIS軟件將nc數據轉換為柵格數據，再利用空間分析工具將逐月數據進行年度累加合成得到年數據。

土壤數據包括土壤深度、植物可利用含水量和土壤飽和導水率，來源于世界土壤數據庫（https：//www.fao.org/home/en/）里提供的HWSD1.2土壤數據集。其中，土壤深度可以直接從該數據集中獲取，植物可利用含水量則利用該數據集提供的土壤中砂粒、粉粒、粘粒以及有機碳含量計算得出［9］，土壤飽和導水率是基于該數據集中的土壤砂粒、粘粒含量等計算獲得［23］。

生物物理系數表（表1）包括土地利用類型（LULC_desc）、代碼（lucode）、植被根系深度（root_depth）、蒸散系數（Kc）及實際蒸散發AET計算公式（LULC_veg），其中蒸散系數參考前人研究結果［2426］和InVEST用戶指南得出。

Zhang系數表征區域降水格局和水文地質特征的氣候季節因子，其值介入1～10。根據太行山區相關參考文獻發現當Zhang系數取5.86時，InVEST模型模擬產水量的效果較好，具有顯著的線性關系［21］。

流速系數以美國農業部自然資源保護局提供的國家工程手冊中的流速坡度景觀表格里的數據為基準值，再乘以1000得到。

地形指數是基于DEM數據，使用ArcGIS軟件中坡度工具計算百分比坡度柵格數據，最后利用公式計算得到［23］。

2.2土地利用轉移矩陣

土地利用轉移矩陣能夠直觀展現研究區域在不同時期不同地類之間的面積變化和轉移方向［27］。采用土地利用類型轉移矩陣可以定量分析太行山區2000—2020年各地類的面積變化及轉移去向，為太行山區土地資源規劃管理提供相應的數據支撐，幫助決策者更好地理解和規劃土地利用。

2.3水源涵養量計算

InVEST模型的產水量模塊是基于水量平衡原理，再結合降水、潛在蒸散發、土壤深度、植被可利用含水量和土地利用等參數計算獲得區域產水量［28］。本文采用InVEST模型3.10.2版本的產水量模塊來計算太行山區的產水量。在得到產水量的基礎上，綜合考慮研究區不同地類的流速系數、地形指數和土壤飽和導水率等因素，計算出太行山區的水源涵養量［911］。

2.4情景模擬法

為了更準確地探究土地利用變化對水源涵養功能的影響，需要排除氣候因素對研究結果產生的干擾。

2.5Pearson相關分析

使用SPSS軟件計算Pearson相關系數來分析不同土地利用類型與水源涵養之間的相關性。Pearson相關系數取值在-1～1，其絕對值大小代表不同土地利用類型對水源涵養的影響程度，越大代表相關性越強［15］。

3結果與分析

3.1土地利用變化分析2000—2020年太行山區土地利用空間分布（圖2）和面積變化（表3）結果表明：（1）從土地利用類型構成看，以耕地、林地、草地為主，多年平均面積占比分別為36.13%、28.97%和26.75%；其次依次為建設用地（5.19%）和水域（1.36%）；未利用地面積最小，僅占0.13%。（2）從土地利用空間分布上看，太行山區耕地主要分布在四周及南部地區，林地主要分布在東北部和西南部地區，草地主要分布在東南部、中部及中北部地區，建設用地則分布在經濟發達的城市地區，水域和未利用地零散分布在研究區內。（3）從土地利用變化趨勢看，20年來各地類的面積變化特征表現為“三增三減”。具體來說，林地、水域和建設用地面積增加，耕地、草地和未利用地面積減少，其中，耕地面積減少最多，由50 729.64 km2減至46 680.6 km2，建設用地面積的增幅最大，為87.23%，由5652.29 km2增至10 582.61 km2。（4）從土地利用類型轉換方式看，研究區內各地類之間轉換頻繁。耕地是太行山區近20年來轉出面積最多的地類，主要轉為草地、建設用地和林地，建設用地轉入面積遠大于轉出面積。

3.2太行山區水源涵養空間格局

2000—2020年太行山區產水量的模擬結果顯示（表4）：太行山區多年平均產水深度和產水總量分別為121.76 mm和165.09×108 m3。在2000—2020年，產水深度和產水總量均呈現先減少后增加、再減少后增加的W型趨勢。到了2020年，產水深度和產水總量達到最大值，分別為16648 mm和225.47×108 m3。而最低值出現在2005年，分別為96.56 mm和131.05×108 m3。

2000—2020年太行山區的水源涵養量結果顯示（表4）：多年平均水源涵養深度和水源涵養總量分別為39.59 mm和53.67×108 m3，2005年水源涵養深度和水源涵養總量達到最低值，分別為29.47 mm和40.00×108 m3，最高值出現在2020年，分別為57.78 mm和78.25×108 m3。在這20年里，水源涵養總量以1.50×108 m3·a-1的速率增長?？梢钥闯觯春B量的變化趨勢與產水量變化一致，呈現出W型變化趨勢。這種W型變化趨勢主要受到氣候變化、人類活動和自然環境等多種因素的綜合影響。在太行山區這一特定地理環境中，降水量的W型變化可能是導致水源涵養量和產水量呈現W型變化趨勢的主要原因。

2000—2020年太行山區水源涵養深度空間格局具有一致性，均呈現出由南向北、由東向西逐漸降低的空間分布特征（圖3）。水源涵養深度作為水源涵養功能的量化指標，能夠很好地表征單元面積上水源涵養量的大小，進而揭示不同地區水源涵養功能的強弱［29］。太行山區的水源涵養深度高值區主要集中在晉城市東部、長治市東南部和忻州市五臺山等地區，這些地區的降水量較為豐富［30］，同時林地和草地的覆蓋率較高，導致蒸散發量相對較少。水源涵養深度的低值區面積較大且主要分布在太行山北部地區，主要受到降水量的顯著影響?？傮w來看，太行山地區降水量和潛在蒸散發量的差值越高的地區，水源涵養功能也相應越高；反之，差值越低的地區，水源涵養功能則相對較低。

3.3不同土地利用類型的水源涵養量

不同地類水源涵養深度反映不同土地利用類型的水源涵養功能強弱。2000—2020年不同土地利用類型的水源涵養深度與水源涵養總量的統計結果顯示（表5）：太行山區土地利用類型水源涵養深度為草地＞林地＞未利用地＞耕地＞建設用地＞水域，表明草地的水源涵養功能最強，其水源涵養深度達46.84 mm，水域的水源涵養功能最弱。林地、草地和耕地是水源涵養總量的主要貢獻地類，它們分別占研究區水源涵養總量的32.78%、31.52%和30.83%。相比之下，建設用地和未利用地提供的水源涵養量僅占總量的3.37%，而水域的水源涵養量始終為0。盡管未利用地的水源涵養功能強于耕地和建設用地，但由于耕地和建設用地的面積遠大于未利用地，因此未利用地的水源涵養總量遠遠低于耕地和建設用地。這表明各地類的水源涵養總量不僅取決于水源涵養功能，還受到分布面積的影響。

3.4土地利用變化對水源涵養功能的影響

將水源涵養情景模擬法的模擬結果與2000年實際情況進行對比（表6），不同情景下的產水深度和產水總量表現出以下趨勢：2000年＞情景1＞情景2＞情景3＞情景4。而對于水源涵養深度和水源涵養總量，則呈現出以下趨勢：情景1＞2000年＞情景2＞情景3＞情景4。

太行山區土地利用類型面積與水源涵養的相關系數表明（表7）：除水域之外，其他地類面積與水源涵養深度和水源涵養總量之間均存在顯著相關性。

土地利用的水源涵養深度表征不同用地類型水源涵養功能的強弱，水源涵養總量表征不同用地類型水源涵養對區域水源涵養貢獻的大小。因此，在水源涵養功能方面，耕地和林地與水源涵養功能呈顯著正相關關系，未利用地、草地和建設用地與水源涵養功能呈顯著負相關關系；相關性大小依次為：未利用地＞耕地＞草地＞建設用地＞林地。在區域水源涵養貢獻方面，耕地、林地和草地與水源涵養貢獻呈顯著正相關關系，建設用地和未利用地與水源涵養貢獻呈顯著負相關關系；相關性大小依次為：未利用地＞耕地＞建設用地＞林地＞草地。

4結論與討論

本文基于InVEST模型的年產水量模塊，結合太行山區的各類數據進行參數本地化調整，進而對該區的水源涵養功能進行評估分析。結果表明，太行山區不同土地類型的水源涵養功能存在顯著差異。土地利用變化主要通過改變地表下墊面，進而影響土壤下滲量、蒸散發量和理化性質等因素，從而影響水源涵養深度［14，31］。植被根系和土壤結構能夠有效地攔截和吸收降水，減緩水流的流速，從而減少水土流失并增強水源涵養功能［12］。草地和林地由于其常年植被覆蓋，具有強大的植被截留、下滲和蓄水能力，因此具有較高的水源涵養功能［3234］。在城鎮化過程中，應重視草地和林地的開發和保護，充分發揮其水源涵養作用。耕地由于作物種植季節的影響，有時會出現植被覆蓋度低、土地裸露的情況，因此水源涵養功能相對較低。建設用地會導致地表下墊面硬化，不透水面面積增加，進而改變下墊面的滲透性、滯水性等特性，對區域水源涵養功能產生影響［3132］。

土地利用變化對區域水源涵養具有重要的影響作用。水源涵養受到氣候和人類活動等多方面的影響［9，1213］。鑒于氣候因素存在不可控性，土地利用變化作為人類可控因素，研究其對水源涵養的影響更具有實踐價值［15，32］。土地利用變化不僅改變了區域地表結構，還對生態系統結構和過程產生了深遠影響，從而影響生態系統服務功能［35］。多年來，太行山區的主要土地利用類型為耕地、林地和草地，合計占據研究區總面積的近九成。山區耕地與建設用地呈現出負相關關系，建設用地面積不斷增加。在2000—2015年，建設用地面積的增長速度相對緩慢，在2015—2020年，增長速度則急劇上升。耕地面積變化趨勢恰好與之相反。這種現象的出現主要是由于城鎮化進程的快速推進以及社會經濟的高速發展，導致周邊大量的耕地和草地被迫轉化為滿足城鎮新增人口生活需求的建設用地［32］。

深入探究土地利用變化對水源涵養功能的影響，不僅有助于協調區域水土資源矛盾，改善生態安全狀況，而且能夠為生態系統管理提供科學依據，從而促進可持續發展。太行山區作為華北平原重要的生態屏障，對京津冀城市群的社會經濟發展起著不可或缺的支撐作用。因此，維護太行山區的水源涵養功能，實施科學的土地資源利用政策具有重要的意義。由于InVEST模型在計算柵格單元產水量時未能考慮地表水和地下水對產水量的影響，可能對研究結果的準確性產生一定的影響。因此，在未來的研究中要進一步加強InVEST模型的參數驗證，以提高模型在不同區域的適應性和可靠性。

參考文獻：

［1］馬尚鈺.鄭州市土地利用/覆被變化預測及其對水源涵養功能的影響［D］.鄭州：華北水利水電大學，2022.

［2］楊思敏，權全，徐家隆，等.延河流域土地利用變化過程中水源涵養功能研究［J］.中國水土保持，2022，485（8）：3336.

［3］杜世勛，郭新亞，榮月靜.基于Budyko假設和SCSCN模型的河源區水源涵養功能研究［J］.水土保持研究，2018，25（1）：147152.

［4］劉景紅，鄭曉，樊俊美，等.基于SWAT模型的渾河中上游水源涵養服務價值評估［J］.應用生態學報，2021，32（11）：39053912.

［5］BESSAH E，RAJI A O，TAIWO O J，et al.Hydrological responses to climate and land use changes：the paradox of regional and local climate effect in the Pra River Basin of Ghana［J］.Journal of Hydrology：Regional Studies，2020，27：100654.

［6］DANESHI A，BROUWER R，NAJAFINEJAD A，et al.Modelling the impacts of climate and land use change on water security in a semiarid forested watershed using InVEST［J］.Journal of Hydrology，2021，593：125621.

［7］ANESEYEE A B，SOROMESSA T，ELIAS E，et al.Evaluation of water provision ecosystem services associated with land use/cover and climate variability in the Winike Watershed，Omo Gibe Basin of Ethiopia［J］.Environmental Management，2022，69（2）：367383.

［8］劉美娟，仲俊濤，王蓓，等.基于InVEST模型的青海湖流域產水功能時空變化及驅動因素分析［J］.地理科學，2023，43（3）：411422.

［9］周雪彤，孫文義，穆興民，等.1990—2020年三江源水源涵養能力時空變化及影響因素［J］.生態學報，2023，43（23）：98449855.

［10］柳冬青，曹二佳，張金茜，等.甘肅白龍江流域水源涵養服務時空格局及其影響因素［J］.自然資源學報，2020，35（7）：17281743.

［11］許麗婷，劉海紅，黃麗潔，等.2000—2020年汾河流域生態環境與水源涵養時空變化［J］.干旱區研究，2023，40（2）：313325.

［12］寧亞洲，張福平，馮起，等.秦嶺水源涵養功能時空變化及其影響因素［J］.生態學雜志，2020，39（9）：30803091.

［13］王冶，薛忠財，王瑞豐，等.1990—2017年承德市水源涵養功能時空變化格局與影響因素研究［J］.生態與農村環境學報，2022，38（9）：11381146.

［14］王陽，魏建兵，程雨露，等.渾河上游新賓縣土地利用演變對水源涵養能力影響評價［J］.水土保持研究，2023，30（4）：169176.

［15］李威，呂思思，趙祖倫，等.土地利用對流域水源涵養及水質凈化的影響：以烏江流域為例［J］.生態學報，2023，43（20）：83758389.

［16］文安邦，湯青，歐陽朝軍，等.中國山地保護與山區發展：回顧與展望［J］.中國科學院院刊，2023，38（3）：376384.

［17］高會，付同剛，梁紅柱，等.太行山區生態系統服務冷熱點區域識別及其權衡/協同關系分析［J］.中國生態農業學報（中英文），2022，30（7）：10451053.

［18］魏靜，劉麗麗，王紅云，等.1990—2020年太行山區土地利用景觀格局時空變化［J］.中國生態農業學報（中英文），2022，30（7）：11231133.

［19］范晨雨，景海濤，王莉，等.太行山區氣候時空變化及其對植被覆蓋度的影響［J］.水土保持研究，2020，27（3）：146152.

［20］王豐，劉金銅，付同剛，等.基于RUSLE模型的太行山區土壤侵蝕時空分異特征及影響因子研究［J］.中國生態農業學報（中英文），2022，30（7）：10641076.

［21］朱建佳，劉金銅，梁紅柱，等.太行山區水資源供需關系的垂直梯度特征［J］.應用生態學報，2019，30（2）：472480.

［22］杜仲瑞.太行山漫記［M］.北京：學苑出版社，2017.

［23］歐陽欣.沱江流域水源涵養功能時空變化格局及未來情景模擬［D］.成都：成都理工大學，2023.

［24］任甜甜.大連市水源供給與水土保持生態系統服務功能評估［D］.大連：遼寧師范大學，2021.

［25］劉曄.永定河流域生態系統服務功能評估及預測［D］.邯鄲：河北工程大學，2023.

［26］王曉晨.基于InVEST模型的陜西秦嶺（渭河流域）生態系統服務評估［D］.西安：西北大學，2023.

［27］曹宇鵬，方江平.基于Fragstats的日喀則市土地利用景觀格局分析［J］.測繪與空間地理信息，2022，45（9）：6972.

［28］陳駿宇，劉鋼，白楊.基于InVEST模型的太湖流域水源涵養服務價值評估［J］.水利經濟，2016，34（2）：2529.

［29］馬靖宣，金曉媚，張緒財，等.基于InVEST模型的張承地區水源涵養功能時空變化特征［J］.水文地質工程地質，2023，50（3）：5464.

［30］楊倩，劉登峰，孟憲萌，等.1960—2017年山西省降水和氣溫的時空變化特征分析［J］.人民珠江，2019，40（6）：2733.

［31］喬亞軍，張慧，韓曉盈，等.基于地理空間角度的黑龍江省水源涵養功能變化的驅動因素分析［J］.生態學報，2023，43（7）：27112721.

［32］劉樹鋒，陳記臣，關帥，等.基于InVEST模型的楊溪河流域土地利用變化對水源涵養功能的影響［J］.科學技術與工程，2022，22（12）：47464751.

［33］呂樂婷，李青，楊勇.遼寧省2001—2020年水源涵養量變化及其影響因素［J］.水土保持通報，2022，42（1）：290296.

［34］薛健，李宗省，馮起，等.1980—2017年祁連山水源涵養量時空變化特征［J］.冰川凍土，2022，44（1）：113.

［35］孫孝平，李雙，余建平，等.基于土地利用變化情景的生態系統服務價值評估：以錢江源國家公園體制試點區為例［J］.生物多樣性，2019，27（1）：5163.