流域土地利用變化及其對生態(tài)系統(tǒng)服務(wù)價(jià)值的影響-以唐河流域山西段為例

段蘭蘭，賀麗，薛占金，劉菊，王琰，嚴(yán)俊霞

（1.山西大學(xué) 黃土高原研究所，山西 太原 030006；2.山西大學(xué) 資源與環(huán)境工程研究所，山西 太原 030006；3.山西省林業(yè)和草原科學(xué)研究院，山西 太原 030006）

0 引言

生態(tài)系統(tǒng)服務(wù)是指為維持生命生存和發(fā)展而從生態(tài)環(huán)境中獲取的食物和原材料的產(chǎn)品和服務(wù)［1］，其永續(xù)發(fā)展對人類發(fā)展和進(jìn)步具有重要影響［2-3］。近年來，對于生態(tài)系統(tǒng)服務(wù)價(jià)值（Ecosystem Service Value，ESV）的研究已有許多［4-6］。早在1997 年，Costanza 提出ESV 分類方法，并建立了量化評估體系［7-8］。在千年生態(tài)系統(tǒng)評估（MA）的基礎(chǔ)上，國外學(xué)者評估了不同流域的水資源、土壤保持等生態(tài)系統(tǒng)服務(wù)價(jià)值［9-10］。與國外相比，我國生態(tài)系統(tǒng)服務(wù)的研究起步較晚，但結(jié)合不同區(qū)域的生態(tài)特征，融合地理學(xué)、生態(tài)學(xué)等方法，開展了不同區(qū)域［11-13］和不同空間［14-17］相關(guān)研究，其中單位面積服務(wù)價(jià)值當(dāng)量因子研究［18］最具代表性。

土地是生存之本、發(fā)展之基、生態(tài)之要。“人類世”（指地球的最近代歷史）以來，人類活動對土地生態(tài)系統(tǒng)干擾日益加劇，土地生態(tài)系統(tǒng)服務(wù)衰減［19］。近年來，土地利用對生態(tài)系統(tǒng)服務(wù)的影響受到國內(nèi)外學(xué)者的關(guān)注，采用主成分分析、空間自相關(guān)理論、引入“三生空間”概念（生產(chǎn)空間、生活空間和生態(tài)空間的簡稱）等方法分析了土地利用類型及生態(tài)系統(tǒng)服務(wù)時(shí)空演變［5-6，9-10］、影響維度［13-16］。已有成果主要集中在土地類型變化劇烈的城市地區(qū)［14-16］，對處于重要生態(tài)功能區(qū)、水源涵養(yǎng)區(qū)的流域上游生態(tài)系統(tǒng)服務(wù)關(guān)注較少。流域生態(tài)系統(tǒng)是由自然、經(jīng)濟(jì)和社會組成的一種復(fù)合生態(tài)系統(tǒng)。通過衡量流域ESV，可以為流域管理和生態(tài)保護(hù)提供理論依據(jù)［20-21］。土地利用變化在一定程度上不僅對生態(tài)系統(tǒng)的類型和分布范圍產(chǎn)生影響，而且影響不同生態(tài)系統(tǒng)服務(wù)類型的強(qiáng)度，從而使總ESV 發(fā)生改變。唐河是大清河的源頭之一，其水土環(huán)境狀況直接關(guān)系到大清河下游的雄安新區(qū)和京津冀地區(qū)可持續(xù)發(fā)展。因此本研究以唐河流域山西段為例，利用2000、2010、2020 年3 期土地利用數(shù)據(jù)，分析2000-2020 年唐河流域土地利用類型的變化特征及其對ESV 影響，確定其主要驅(qū)動因素，為唐河流域提高生態(tài)環(huán)境質(zhì)量、恢復(fù)生態(tài)系統(tǒng)服務(wù)提供科學(xué)依據(jù)。

1 研究區(qū)概況與數(shù)據(jù)來源

1.1 研究區(qū)概況

唐河發(fā)源于山西省渾源縣千佛嶺鄉(xiāng)，在靈丘縣下北泉村進(jìn)入河北省淶源縣。在山西省內(nèi)，唐河流域地理位置為：39°09′N—39°38′N，113°42′E—114°30′E，河長為96 km，流域面積為2193 km2（圖1）。地勢呈西北高、東南低的特點(diǎn)。氣候?yàn)闇貛О敫珊荡箨懶约撅L(fēng)氣候，1 月平均氣溫為-12 ℃，7 月平均氣溫為22 ℃，≥10 ℃的積溫平均為2600 ～3100 ℃，多年平均溫為6.5 ℃～7.0 ℃，年均降水量500 mm～550 mm，年均蒸發(fā)量1000 mm～1050 mm［22］，多年平均無霜期為73 d～217 d。區(qū)內(nèi)優(yōu)勢植物有華北落葉松（Larixprincipis-rupprechtii）、油松（PinustabuliformisCarr）、側(cè)柏（Platycladusorientalis）、楊樹（Populusspp.）、刺槐（Robiniapseudoacacia）、繡線菊（Spiraeasalicifolia）等。土壤類型以潮土和褐土為主。

圖1 唐河流域山西段地理位置Fig.1 Geographical position of Tanghe river basin in Shanxi province

1.2 數(shù)據(jù)來源與處理

本文所用的土地利用類型數(shù)據(jù)主要來源于國家基礎(chǔ)地理信息中心（http：//www.ngcc.cn/ngcc/），分辨率為30 m×30 m。首先利用ArcGIS 10.2 對2000、2010 和2020 年3 期土地利用數(shù)據(jù)進(jìn)行處理，具體步驟為：①利用ArcToolbox 設(shè)為空函數(shù)工具，將柵格數(shù)據(jù)去除黑邊；②利用Mosaic to new raster 進(jìn)行合并；③利用Extract by the mask將合并后的數(shù)據(jù)進(jìn)行裁剪，得到研究區(qū)范圍和土地利用數(shù)據(jù)；④利用Reclassify tool 將土地類型分為6 類：林地、耕地、草地、水域、建設(shè)用地和未利用地，根據(jù)屬性表和空間分辨率計(jì)算各地類面積；⑤利用Raster calculator 得到土地轉(zhuǎn)移矩陣。

年均氣溫?cái)?shù)據(jù)來自worldclim.org；年均降水量數(shù)據(jù)來源于國家氣象科學(xué)數(shù)據(jù)中心（http：//data.cma.cn/）；單位糧食產(chǎn)量數(shù)據(jù)來源于《山西省統(tǒng)計(jì)年鑒》，糧食價(jià)格來源于《中國農(nóng)產(chǎn)品價(jià)格調(diào)查年鑒》；GDP、城鎮(zhèn)人口和農(nóng)村人口來源于山西省和大同市統(tǒng)計(jì)年鑒。

2 研究方法

2.1 土地利用變化分析

（1）凈變化量表示特定區(qū)域特定時(shí)間段內(nèi)某一地類的量的變化，計(jì)算公式為：

式中：Ci，Uia和Uib分別為第i種地類在研究期的凈變化量（km2）、研究初期面積（km2）、研究末期面積（km2）。

（2）土地利用動態(tài)度表示特定區(qū)域特定時(shí)間段內(nèi)不同地類的變化速率，其絕對值的大小代表了面積變化速率的幅度，數(shù)值越大，速率越快。分為單一和綜合土地利用動態(tài)度2 類，計(jì)算公式［23］分別為：

式中：K為某一地類的土地利用動態(tài)度（%）；Ua為該地類在研究初期面積（km2）；Ub為該地類在研究末期面積（km2）；T為研究時(shí)段（a）。

式中：LUi為第i種地類在研究初期面積（km2）；ΔLUi-j為研究期內(nèi)第i種地類轉(zhuǎn)變?yōu)槠渌仡惷娣e（km2，絕對值）；T為研究時(shí)段（a）。

（3）土地利用轉(zhuǎn)移矩陣表示不同地類之間轉(zhuǎn)移方向和數(shù)量的變化［24］，計(jì)算公式為：

式中：Aij為第i種地類在研究時(shí)段內(nèi)演替為第j種地類面積（km2）；i（i=1，2，…，n），j（j=1，2，…，n）分別為研究初期和末期的地類，本研究中n=6。

2.2 生態(tài)系統(tǒng)服務(wù)價(jià)值

2.2.1 生態(tài)系統(tǒng)服務(wù)價(jià)值評估

本文在Costanza 等［9］研究基礎(chǔ)上，參考謝高地等制定的價(jià)值當(dāng)量因子表［25-27］，將其作為修正系數(shù)，分別對6 種地類的ESV 進(jìn)行修正。通過查閱統(tǒng)計(jì)年鑒，計(jì)算得出唐河流域山西段的平均糧食單位產(chǎn)量為5 734 kg/hm2及山西省的糧食收購均價(jià)為2.78 元/kg，全國平均糧食單產(chǎn)市場價(jià)格的1/7 等于1 個當(dāng)量因子的價(jià)值［26］，所以該流域1 個當(dāng)量因子的經(jīng)濟(jì)價(jià)值量為2 277.22 元/hm2。根據(jù)研究區(qū)1 個當(dāng)量因子的價(jià)值可以計(jì)算出研究區(qū)單位面積ESV（表1）［28］。

表1 唐河流域山西段不同土地利用類型生態(tài)系統(tǒng)當(dāng)量因子（yuan/hm2）Table 1 Ecosystematic equivalent factors of different land use in Tanghe river basin in Shanxi province (yuan/hm2)

唐河流域山西段各時(shí)期的ESV 計(jì)算公式為：

式中：ESV為生態(tài)系統(tǒng)服務(wù)的總價(jià)值（元）；Ai為第i種地類面積（hm2）；VCi為第i種地類的單位面積生態(tài)服務(wù)價(jià)值（元/hm2·a-1）。

2.2.2 生態(tài)貢獻(xiàn)率

生態(tài)貢獻(xiàn)率是指在特定時(shí)間段內(nèi)不同地類ESV 變化量與總ESV 變化量的比值。通過計(jì)算生態(tài)貢獻(xiàn)率可以確定影響研究區(qū)ESV 變化的主要因子。計(jì)算公式為：

式中：Si，ΔESVi分別為研究時(shí)段第i種地類的貢獻(xiàn)率（%）和生態(tài)服務(wù)價(jià)值變化量（元）；i為地類。

2.2.3 人為影響綜合指數(shù)

人類活動在景觀格局的變化中發(fā)揮著重要的作用。為了分析人類活動對ESV 的影響，本文利用人為影響綜合指數(shù)定量分析一定區(qū)域內(nèi)景觀總體受人類干擾的強(qiáng)度［29］。計(jì)算公式為：

式中：HAI表示人為影響綜合指數(shù)；N為地類的種類；Ai為第i種地類面積（hm2）；Pi為第i種地類所反映的人為影響強(qiáng)度系數(shù)；TA為研究區(qū)總面積（hm2）。

為了使計(jì)算結(jié)果更準(zhǔn)確，本文通過Lohani清單法、Leopold 矩陣法和Delphi 法分別計(jì)算得出不同地類的人為影響強(qiáng)度系數(shù)，再通過求三者的平均值，最終得出研究區(qū)不同地類的人為影響強(qiáng)度系數(shù)Pi（表2）。

表2 唐河流域山西段不同地類人為影響強(qiáng)度系數(shù)Table 2 Intensity coefficients of anthropogenic influence of different land use in Tanghe river basin in Shanxi province

HAI值取值范圍為0 ～ 1，值越大意味著人為活動對景觀組分干擾越大，HAI＞0.80、0.60＜HAI＜ 0.80、0.40＜HAI≤ 0.60、0.20 ＜HAI≤ 0.40、0 ＜HAI≤ 0.20 分別表示景觀組分受到人類干擾程度為高、較高、中、較低和低。

3 結(jié)果與分析

3.1 土地利用變化

3.1.1 土地利用類型的時(shí)空變化

由表3 可知，20 年來研究區(qū)內(nèi)每種土地利用類型的面積均發(fā)生了變化。可以看出，林地和耕地是研究區(qū)的主要土地利用類型，二者約占總面積的50%以上（圖2）。林地多集中在研究區(qū)的外圍山區(qū)，其面積從2000 年的722.02 km2（占34.90%），增加到2020 年的739.91 km2（占35.80%）。耕地多分布于研究區(qū)中部，20年間面積從736.96 km2減少到690.87 km2。2000-2020 年草地面積變化不大，減少了8.51 km2；建設(shè)用地面積增加了34.20 km2。

表3 2000、2010和2020年唐河流域山西段不同土地類型面積變化Table 3 Area change in different land use in Tanghe river basin in Shanxi province in 2000, 2010 and 2020

圖2 2000、2010和2020年唐河流域山西段土地利用類型Fig.2 Land use type of Tanghe river basin in Shanxi province in 2000, 2010 and 2020

3.1.2 土地利用動態(tài)度

從單一土地利用動態(tài)度來看，2000-2010、2010-2020、2000-2020 年研究區(qū)內(nèi)動態(tài)度變化較大的是水域和建設(shè)用地，分別為-5.00%和3.47%、86.00%和8.72%、19.00%和7.61%；林地、草地、耕地、未利用地的動態(tài)度較小，均低于0.40%（表4）。2010-2020 年不同土地利用動態(tài)度差異較大，為0～86.00%，主要原因是水域和建設(shè)用地的土地利用動態(tài)度較高，分別為86.00%和8.72%。整體來看，20 年間唐河流域山西段綜合土地利用動態(tài)度為0.13%，其中2000-2010、2010-2020 年的綜合土地利用動態(tài)度分別為0.11%和0.25%，后者變化更劇烈。

表4 2000-2020年唐河流域山西段土地利用動態(tài)度（%）Table 4 Dynamic of different land use in Tanghe river basin in Shanxi province during 2000-2020 （%）

3.1.3 土地利用面積轉(zhuǎn)移

從2000-2010、2010-2020、2000-2020 年研究區(qū)土地利用轉(zhuǎn)移矩陣（表5-7）可以看出，2000-2010 年草地轉(zhuǎn)入的土地面積最多（58.34 km2），其次是林地（35.55 km2）；轉(zhuǎn)出的面積中，草地仍為最多（43.09 km2），其次為耕地（41.48 km2）。2010-2020 年林地轉(zhuǎn)入的面積最多，為85.08 km2，其次是草地（80.94 km2）；轉(zhuǎn)出的面積中，草地最多（103.72 km2），其次是林地（62.38 km2）。在兩個時(shí)間段轉(zhuǎn)入和轉(zhuǎn)出面積中，水域和未利用地都較少。整體來看，20 年間草地轉(zhuǎn)入和轉(zhuǎn)出的面積均最多，分別為119.17 km2和110.66 km2。

表5 2000-2010年唐河流域山西段土地利用轉(zhuǎn)移矩陣（km2）Table 5 Land use transfer matrix of Tanghe river basin in Shanxi province from 2000 to 2010 (km2)

表6 2010-2020年唐河流域山西段土地利用轉(zhuǎn)移矩陣（km2）Table 6 Land use transfer matrix of Tanghe river basin in Shanxi province from 2010 to 2020 (km2)

表7 2000-2020年唐河流域山西段土地利用轉(zhuǎn)移矩陣（km2）Table 7 Land use transfer matrix of Tanghe river basin in Shanxi province from 2000 to 2020 (km2)

在2000-2010、2010-2020、2000-2020 年，草地、林地和耕地三者之間發(fā)生了明顯的土地利用轉(zhuǎn)移。其中草地向林地轉(zhuǎn)移的面積最多，分別為28.60、79.75、88.37 km2；其次是林地轉(zhuǎn)為草地，分別為28.60、54.76、63.94 km2。

3.2 生態(tài)系統(tǒng)服務(wù)價(jià)值變化

3.2.1 各地類生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)間變化

2000、2010、2020 年研究區(qū)總ESV 分別為1 800 807.88、1 803 853.41 和1 822 512.73 萬元（表8）。在不同土地利用類型中，林地、草地和耕地具有較高的ESV 貢獻(xiàn)率，分別為70.19%、26.18%、3.61% 和69.66%、26.81%、3.51% 和71.08%、25.49%、3.34%，建設(shè)用地和未利用地的ESV 較低。2000-2020 年研究區(qū)的總ESV呈逐漸上升趨勢，增加了21 704.85 萬元，增幅為0.01%，其中林地、水域和建設(shè)用地的ESV分別增加了31 320.50、1 076.23 和155.77 萬元；耕地和草地提供的ESV 均有所降低，分別減少了4 061.88 和6 837.77 萬元。2010 年與2000 年相比，除草地提供的ESV 增加了12 255.11 萬元外，耕地、林地和水域提供的ESV 均有所降低，分別減少了1 653.81、7 407.32 和148.45 萬元，研究區(qū)總ESV 增加了3 405.53 萬元。2020 年與2010 年相比，林地、水域、建設(shè)用地和未利用地提供的ESV 分別增加了38 727.83、1224.68、155.77 和52.00 萬元，耕地和草地提供的ESV均有所減低，分別減少了2 408.06 和19 092.88萬元，研究區(qū)總ESV 增加了18 659.32 萬元（表8）?？傮w來看，20 年間研究區(qū)ESV 增加的主要原因是林地、水域和建設(shè)用地面積的增加。

表8 2000-2020年唐河流域山西段各地類生態(tài)系統(tǒng)服務(wù)價(jià)值Table 8 ESV of different land use in Tanghe river basin in Shanxi province from 2000 to 2020

3.2.2 單項(xiàng)生態(tài)系統(tǒng)服務(wù)類型價(jià)值變化

從一級服務(wù)功能來看，2000、2010、2020 年占總ESV 的比例從高到低依次為調(diào)節(jié)、支持、供給和文化服務(wù)。其中占比最高的是調(diào)節(jié)服務(wù)，均在60% 左右，其次為支持服務(wù)，占比在25%左右，供給和文化服務(wù)的占比較低，約為0.05%。2000-2020 年，調(diào)節(jié)、支持和文化服務(wù)呈持續(xù)增加趨勢，2010 年比2000 年分別增加了2 356.70、714.76 和229.40 萬元，2020 年比2010年分別增加了14 040.45、3 740.26 和903.22 萬元，2020 年比2000 年分別共增加了16 397.15、4 455.02 和1 132.62 萬元。供給服務(wù)呈持續(xù)減少趨勢，2010 年比2000 年減少255.33 萬元，2020 年比2010 年減少24.61 萬元，20 年共減少了279.94 萬元（圖3）。

圖3 2000—2020年唐河流域山西段生態(tài)系統(tǒng)一級服務(wù)功能價(jià)值及其變化Fig.3 ESV and its change in Tanghe river basin in Shanxi province from 2000 to 2020

從二級服務(wù)功能來看，2000-2010 年除食物生產(chǎn)和原材料生產(chǎn)價(jià)值降低了218.37 萬元和36.96 萬元以外，其他二級服務(wù)功能的ESV 增加了195.52 萬元～967.45 萬元；2010-2020 年僅食物生產(chǎn)價(jià)值降低了390.68 萬元，其他二級服務(wù)功能的ESV 增加了366.06 萬元～4 666.28 萬元。過去20 年間原材料生產(chǎn)價(jià)值呈先減少后增加趨勢；食物生產(chǎn)價(jià)值呈持續(xù)降低趨勢，年均減少30.45 萬元；氣體調(diào)節(jié)、氣候調(diào)節(jié)、水源涵養(yǎng)、廢物處理、保持土壤、維持生物多樣性和提供美學(xué)景觀呈持續(xù)增加趨勢，每年分別增加86.63、354.87、270.85、107.51、93.94、128.81 和56.63 萬元（表9）。

表9 2000-2020年唐河流域山西段土地生態(tài)系統(tǒng)單項(xiàng)服務(wù)價(jià)值Table 9 Individual ESV of different land use in Tanghe river basin in Shanxi province from 2000 to 2020

3.2.3 各地類對生態(tài)系統(tǒng)服務(wù)價(jià)值變化的生態(tài)貢獻(xiàn)率

2000-2010、2010-2020、2000-2020 年研究區(qū)ESV 變化貢獻(xiàn)率較高的是耕地、林地和草地，分別為6.65%、29.78%、49.28%和3.00%、48.22%、23.77% 和6.23%、48.03%、10.49%（表10），3 個地類生態(tài)貢獻(xiàn)率合計(jì)在98% 以上；建設(shè)用地和未利用地的生態(tài)貢獻(xiàn)率較低。表明影響研究區(qū)總ESV 變化最大的土地利用類型是耕地、林地和草地，三者是研究區(qū)ESV變化的主要貢獻(xiàn)因子和影響因子。

表10 2000-2020年唐河流域山西段各地類對ESV變化的生態(tài)貢獻(xiàn)率（%）Table 10 Ecological contribution to ESV change in different land use in Shanxi province in 2000-2020 (%)

3.3 驅(qū)動力分析

從以上分析可以看出，20 年間唐河流域的總ESV 增加了21 704.85 萬元，這與環(huán)境變化和人類活動有很大的關(guān)系。已有研究表明，自然環(huán)境變化和人類活動在一定程度上影響著物種的正常生長和群落演替［30］。

3.3.1 人為活動

為了更清楚了解人為活動對ESV 的影響。本文首先將研究區(qū)劃分為5 km×5 km 的格網(wǎng)，并分別計(jì)算每個格網(wǎng)的人為影響強(qiáng)度，將其作為每個格網(wǎng)中心點(diǎn)的值；然后通過趨勢分析和正態(tài)檢驗(yàn)，利用ArcGIS 10.2 里Kriging 空間插值方法得到唐河流域山西段2000、2010 和2020年3 期人為影響強(qiáng)度的空間分布圖（圖4）。

圖4 唐河流域山西段三期人為干擾綜合強(qiáng)度分布Fig.4 Distribution of comprehensive intensity of human disturbance over three periods in Tanghe river basin in Shanxi province

從圖4 可以看出，唐河流域山西段人為干擾綜合強(qiáng)度的空間分布主要表現(xiàn)為：中高影響和高影響強(qiáng)度占主要地位，主要分布在研究區(qū)中部；低影響和中影響強(qiáng)度呈小塊狀、零星分布；中低影響強(qiáng)度主要分布在研究區(qū)東南部。過去20 年間低影響和中低影響強(qiáng)度基本保持不變，中影響強(qiáng)度范圍有所擴(kuò)張，中高影響和高影響強(qiáng)度區(qū)域范圍有所減少。進(jìn)一步分析表明，干擾強(qiáng)度高的區(qū)域主要土地利用類型為耕地；干擾強(qiáng)度低、中低的區(qū)域主要土地利用類型是林地和草地；干擾強(qiáng)度大的區(qū)域由大量耕地、草地轉(zhuǎn)為林地，促使ESV 增加，與前述唐河流域山西段ESV 的變化規(guī)律相一致。

3.3.2 自然因素

溫度和降水不僅是影響植物生長最直接的環(huán)境因子，還會影響生態(tài)系統(tǒng)的正常演替。為進(jìn)一步分析引起唐河流域山西段ESV 變化的自然驅(qū)動力，本文利用研究區(qū)的溫度和降水與年份之間的相關(guān)性進(jìn)行回歸分析（圖5）。

圖5 2000-2020年唐河流域山西段氣候變化特征Fig.5 Characteristics of climatic change from 2000 to 2020 in Tanghe river basin in Shanxi province

從圖5 可以看出，過去20 年來唐河流域山西段的年降水量呈上升趨勢，上升率為10.816 mm/10 a；年均溫也呈上升趨勢。主要表現(xiàn)為：全球氣候變暖引起溫度帶北移，對唐河流域山西段的影響主要體現(xiàn)為降水量增加；年降水量增加使土壤含水量上升，植被存活率提高，ESV 提高。

4 討論

過去20 年來，研究區(qū)土地利用類型的結(jié)構(gòu)和分布發(fā)生了明顯變化，其中耕地面積減小，林地、建設(shè)用地和水域面積增加。其原因有：①農(nóng)村人口減少導(dǎo)致耕地面積減少，調(diào)查資料顯示，2000-2020 年研究區(qū)城鎮(zhèn)人口增加了80 901 人，農(nóng)村人口減少了58 134 人。城鎮(zhèn)化過程中，大批的農(nóng)村人口向城鎮(zhèn)遷移，農(nóng)村許多耕地撂荒，從而導(dǎo)致耕地面積減少、建設(shè)用地面積增加；②三北防護(hù)林、退耕還林還草、京津風(fēng)沙源治理等生態(tài)政策和工程的實(shí)施，使大量耕地轉(zhuǎn)變?yōu)榱植莸?，林草地面積增加［31-32］，生態(tài)環(huán)境趨于好轉(zhuǎn)，20 年間研究區(qū)共完成退耕還林還草4 226.23 hm2；③研究區(qū)在原有生態(tài)保護(hù)措施基礎(chǔ)上，通過劃定唐河水源保護(hù)區(qū)、關(guān)停煤礦和鐵礦、低產(chǎn)地（如：坡耕地、沙化地等）改造為草地和林地等手段，實(shí)現(xiàn)了凈水保土增綠的環(huán)境效益，極大地改善了流域生態(tài)環(huán)境功能，提升了流域生態(tài)環(huán)境質(zhì)量。在本研究中，2000-2020 年研究區(qū)總ESV 增加了21 704.85萬元，增幅為0.01%，這主要得益于草地向林地的轉(zhuǎn)移，這與陳明葉等［33］在大清河阜平流域的研究結(jié)果一致。

自然和人為因素均會引起ESV 的變化。在自然因素中，溫度和降水直接影響著生態(tài)系統(tǒng)中的各項(xiàng)服務(wù)功能。Shen 等［34］通過研究京津冀地區(qū)的生態(tài)系統(tǒng)服務(wù)變化，表明植被指數(shù)是其主要影響因子。但自然因素（如：氣溫、降水等）在短期內(nèi)變化幅度不大，對ESV 影響不明顯［35］。趙志剛等［36］研究發(fā)現(xiàn)鄱陽湖生態(tài)經(jīng)濟(jì)區(qū)影響ESV 變化的主要驅(qū)動因子是城鎮(zhèn)化率，其次是城鎮(zhèn)人口、人口密度和GDP。茆長寶等［37］研究發(fā)現(xiàn)影響土地利用變化的主要因素是人口數(shù)量和質(zhì)量。在本研究中，人類活動是影響唐河流域土地利用變化的主導(dǎo)因素，這與祝冰潔等［38］對大清河山區(qū)的研究結(jié)果一致。

與大清河流域已有研究成果相比，本文的計(jì)算結(jié)果與其他研究結(jié)果有一定的差別，主要原因是研究方法不同，本文采用“價(jià)值當(dāng)量因子法”，結(jié)合唐河流域的實(shí)際情況對當(dāng)量因子進(jìn)行修正，使得研究結(jié)果與實(shí)際較為吻合。潘瑩等［39］采用InVEST 模型對大清河流域生態(tài)系統(tǒng)服務(wù)進(jìn)行評估，使用該模型應(yīng)評價(jià)其適用性，且如何選擇模型參數(shù)并進(jìn)行本地化校正是關(guān)鍵問題。

需要指出的是，本研究重點(diǎn)分析了唐河流域山西段土地利用變化及其對ESV 變化的影響，并探究了影響唐河流域山西段ESV 變化的驅(qū)動因素。限于資料的獲取性，本研究重點(diǎn)分析了溫度、降水和人為干擾情況對ESV 變化的影響。未來，應(yīng)考慮地形、坡向、城鎮(zhèn)和農(nóng)村人口數(shù)量、質(zhì)量以及GDP 等因素對ESV 變化的影響。此外，應(yīng)從不同時(shí)空維度深入分析唐河流域山西段生態(tài)系統(tǒng)服務(wù)的權(quán)衡與協(xié)同關(guān)系，探究唐河流域山西段的生態(tài)補(bǔ)償機(jī)制，為制定流域生態(tài)保護(hù)政策提供參考依據(jù)。

5 結(jié)論

（1）2000-2020 年，唐河流域山西段的主要地類是耕地和林地。20 年來，耕地從736.96 km2減少到690.87 km2；草地減少了8.51 km2；林地的面積從722.02 km2增加到739.91 km2；建設(shè)用地增加了34.20 km2；水域和未利用地的面積變化不大。20 年間土地利用轉(zhuǎn)移主要發(fā)生在林地、草地、耕地三者之間，其中草地約有88.37 km2轉(zhuǎn)為耕地。

（2）2000-2020 年唐河流域山西段土地利用變化最大的是水域，動態(tài)度為19.00%，面積增加了0.38 km2；其次是建設(shè)用地，動態(tài)度為7.61%，面積增加了34.20 km2。20 年間，唐河流域山西段綜合土地利用動態(tài)度為0.13%，其中2000-2010、2010-2020 年的綜合土地利用動態(tài)度分別為0.11% 和0.25%，2010-2020 年變化更劇烈。

（3）2000-2020 年唐河流域山西段總ESV呈逐漸增加趨勢，20 年增加了21 704.85 萬元，年均凈增加值達(dá)5 246.98 元/km2。在一級服務(wù)價(jià)值中，除供給服務(wù)呈先減后增趨勢，調(diào)節(jié)、支持和文化服務(wù)所提供的價(jià)值均呈持續(xù)增加趨勢。二級生態(tài)服務(wù)價(jià)值以氣候調(diào)節(jié)和水源涵養(yǎng)服務(wù)為主，分別增加了7 097.30 萬元和5 416.96 萬元。

（4）2000-2020 年研究區(qū)ESV 變化貢獻(xiàn)率較高的是林地、耕地和草地，分別是48.03%、6.23%，10.49%；且三者生態(tài)貢獻(xiàn)率合計(jì)在98%以上，說明林地、耕地和草地是影響研究區(qū)ESV 變化的主要因子。

（5）生態(tài)系統(tǒng)服務(wù)的發(fā)展變化是自然因素和人為活動共同作用的結(jié)果，通過分析唐河流域山西段ESV 變化及其驅(qū)動因素，人為活動是唐河流域山西段ESV 變化的主要驅(qū)動因素。