漓江流域土地利用時空變化對生態(tài)系統(tǒng)服務(wù)價(jià)值的影響

李 寧,王浩宇,王金葉,黃良珂

(1.桂林理工大學(xué) a.環(huán)境科學(xué)與工程學(xué)院;b.廣西空間信息與測繪重點(diǎn)實(shí)驗(yàn)室,廣西 桂林 541006;2.桂林理工大學(xué)博文管理學(xué)院,廣西 桂林 541006)

0 引 言

Daily[1]提出了生態(tài)系統(tǒng)服務(wù)價(jià)值相關(guān)理論和研究方法,Costanza等[2]計(jì)算生態(tài)系統(tǒng)服務(wù)價(jià)值的方法引起國內(nèi)外學(xué)者的廣泛關(guān)注,生態(tài)系統(tǒng)服務(wù)價(jià)值的量化對于生態(tài)學(xué)、經(jīng)濟(jì)學(xué)及社會學(xué)的連接具有里程碑的意義。歐陽志云等[3]、謝高地等[4]對生態(tài)系統(tǒng)服務(wù)內(nèi)涵的發(fā)展作出重要的貢獻(xiàn)。朱春全[5]提出把自然生態(tài)系統(tǒng)的生產(chǎn)總值即生態(tài)系統(tǒng)提供的最終產(chǎn)品和服務(wù)價(jià)值的總和納入中國可持續(xù)發(fā)展的評估核算體系。Costanza等[6]將生態(tài)學(xué)、生態(tài)經(jīng)濟(jì)學(xué)和社會學(xué)結(jié)合為一個統(tǒng)一的生態(tài)系統(tǒng)服務(wù)框架,其核心目標(biāo)是造福人類社會。

生態(tài)系統(tǒng)作為生物圈最重要的組成部分,在平衡生命支持系統(tǒng)和環(huán)境方面發(fā)揮著不可替代的作用[7-8]。生態(tài)系統(tǒng)服務(wù)通常被定義為其對人類提供或部分提供有價(jià)值的商品和服務(wù)[9],是指來自生態(tài)系統(tǒng)功能的生物所帶來的益處,這些功能對人類福祉、健康、生計(jì)和生存都是必要的[10-12]。作為理解人類活動與環(huán)境相互作用的最重要指標(biāo)之一[13],土地利用與土地覆蓋變化在全球生態(tài)與環(huán)境變化中扮演著重要的角色[14-16]。生態(tài)系統(tǒng)服務(wù)價(jià)值(ESV)是指用經(jīng)濟(jì)法則對生態(tài)系統(tǒng)的服務(wù)和自然資本進(jìn)行估算。

漓江流域是世界重要的旅游目的地, 形成了以旅游為主線的社會經(jīng)濟(jì)系統(tǒng), 特別是生態(tài)旅游又是以生態(tài)資源環(huán)境為前提進(jìn)行的一系列活動,因此,該區(qū)域的生態(tài)與環(huán)境關(guān)系著流域的可持續(xù)發(fā)展。近年來,隨著城鎮(zhèn)化速度的加快及旅游活動的迅速擴(kuò)張,漓江流域的土地利用景觀格局發(fā)生了較大變化,流域面臨水土環(huán)境污染、生物多樣性降低、土地退化等眾多生態(tài)與環(huán)境問題[17-18]。已有學(xué)者對漓江流域的環(huán)境相關(guān)情況進(jìn)行初步研究[19-21],但是從土地利用與土地覆蓋變化角度研究生態(tài)系統(tǒng)服務(wù)價(jià)值的文獻(xiàn)較少。本文以遙感數(shù)據(jù)為基礎(chǔ),基于生態(tài)與環(huán)境及生態(tài)經(jīng)濟(jì)理論,應(yīng)用3S技術(shù)對漓江流域的土地利用與土地覆蓋的時空變化及生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行研究,為漓江流域的植被保護(hù)、水資源保護(hù)、景觀保護(hù)及土地開發(fā)利用提供科學(xué)依據(jù)。

1 研究區(qū)域概況

漓江流域地處北緯24°38′10″—25°53′59″,東經(jīng)110°07′39″—110°42′57″,流域面積583 150 hm2,涉及廣西桂林市疊彩區(qū)、秀峰區(qū)、七星區(qū)、象山區(qū)、雁山區(qū)全境以及資源縣、興安縣、靈川縣、臨桂縣、陽朔縣、平樂縣的部分區(qū)域(圖1)。漓江是珠江流域的上游,是桂林的母親河,漓江流域蘊(yùn)藏有豐富的歷史遺存、優(yōu)美的自然景觀和良好的生態(tài)與環(huán)境,特別是全球分布面積最大、景觀美學(xué)價(jià)值最高、最具典型性和不可替代性的喀斯特地貌景觀，是桂林山水的靈魂和精華。漓江上游的貓兒山、青獅潭、海洋山森林覆蓋率較高,高海拔地帶的森林植被保存較好,這些植被特別是常綠闊葉林成為漓江源頭的“綠色水庫”。漓江上游地區(qū)山體高大,受低緯度和海拔變化的影響,產(chǎn)生水熱資源的再分配,造就了以中亞熱帶為基帶的暖溫帶、溫帶的垂直立體氣候特征,為上游地區(qū)不同植被的發(fā)育提供了優(yōu)越的條件,是廣西三大林區(qū)之一。耕地和林地是流域內(nèi)兩種主要的土地覆蓋類型。

2 數(shù)據(jù)來源、處理和研究方法

2.1 數(shù)據(jù)來源

從美國地質(zhì)勘探局(USGS)下載條帶號(path)為124和125、 行編號(row)為42和43，云量少于10%的2000年Landsat 5 TM,2010年Landsat 7 ETM+和2016年Landsat 8 OLI 3期的數(shù)據(jù), 每期4景圖像, 共12景圖像。 為最大程度減小誤差, 所有的遙感數(shù)據(jù)大部分取自流域雨量較多的4—8月份, 具體數(shù)據(jù)見表1。

2.2 數(shù)據(jù)處理

利用ENVI 5.3和ArcGIS 10.0進(jìn)行影像處理, 數(shù)據(jù)均使用WGS84坐標(biāo)系, 經(jīng)過波段合成、 拼接、裁剪等預(yù)處理過程。 由于Landsat 7 ETM+機(jī)載掃描行校正器(scan lines corrector, SLC) 于2003年5月31日突然發(fā)生故障,掃描的影像數(shù)據(jù)出現(xiàn)了重疊和部分缺失。因此,本文使用縫隙填充的方式將丟失的數(shù)據(jù)填補(bǔ)出來。修復(fù)算法是自適應(yīng)局部回歸法,原理是把不同獲取時間、與原圖地理坐標(biāo)相同的影像數(shù)據(jù)提取出來。使用這種符合條件的數(shù)據(jù)來填充有壞道的影像,使有壞道的影像可以用于實(shí)驗(yàn)研究。借助Google Earth,使用支持向量機(jī)(SVM)分類方法對影像進(jìn)行監(jiān)督分類,并對分類結(jié)果進(jìn)行小斑塊合并及斑塊平滑等處理。為確定分類的精度和可靠性,借助手持GPS進(jìn)行野外實(shí)地調(diào)查,結(jié)合Google Earth及原始影像獲取真實(shí)的感興趣區(qū)驗(yàn)證樣本,使用混淆矩陣的方法,對分類結(jié)果進(jìn)行精度驗(yàn)證。其中，2000、2010和2016年監(jiān)督分類影像精度驗(yàn)證結(jié)果的分類精度分別為93%、94%、93%,Kappa系數(shù)分別為0.89、0.91和0.92,最后劃分出耕地、林地、水域和建設(shè)用地4種土地生態(tài)類型，并對各土地利用類型面積進(jìn)行計(jì)算。

圖1 漓江流域地理位置Fig.1 Map of the Lijiang River basin

表1 數(shù)據(jù)信息

2.3 計(jì)算方法

謝高地等[22-23]在Costanza等[2]的研究基礎(chǔ)上,根據(jù)中國的實(shí)際情況制定出了中國生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量表和單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值表。由于生態(tài)系統(tǒng)的復(fù)雜性和地域的區(qū)別性,針對不同的地域特點(diǎn),使用同一個標(biāo)準(zhǔn)計(jì)算生態(tài)系統(tǒng)服務(wù)價(jià)值會帶來一定的誤差。因此,本文把每種土地類型與最接近的生態(tài)系統(tǒng)服務(wù)類型連接起來,林地、水域、耕地分別對應(yīng)森林、水體、農(nóng)田。其中建設(shè)用地生態(tài)系統(tǒng)服務(wù)價(jià)值取值為0,基于生態(tài)服務(wù)功能強(qiáng)度與生物量呈線性關(guān)系的假定[23],修訂耕地單位面積食物生產(chǎn)生態(tài)服務(wù)價(jià)值(VC):

VC=(b/B)p,

(1)

其中,VC為修正后的單位面積耕地食物生產(chǎn)生態(tài)服務(wù)價(jià)值；b為耕地生態(tài)系統(tǒng)的生物量因子；B為全國耕地生態(tài)系統(tǒng)生物量因子,其值為1.00；p為耕地生態(tài)系統(tǒng)食物生產(chǎn)的基準(zhǔn)單價(jià)。 漓江流域耕地生物量因子取廣西的均值0.98[24], 即VC=(0.98/1.00)p=0.98p。

生態(tài)服務(wù)劃分為氣體調(diào)節(jié)、氣候調(diào)節(jié)、水源涵養(yǎng)、土壤形成與保護(hù)、廢物處理、生物多樣性維持、食物生產(chǎn)、原材料生產(chǎn)、休閑娛樂共9類。根據(jù)中國陸地生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值當(dāng)量表以及中國不同陸地生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值表[23],修正后的漓江流域單位面積耕地食物生產(chǎn)生態(tài)系統(tǒng)服務(wù)價(jià)值(VC=867.20元/hm2)乘以中國陸地生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值當(dāng)量即可得到流域不同生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值(表2)。

表2 漓江流域生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值

Table 2 Ecological service value of per unit area of the ecosystem in Lijiang River basin 元/hm2

生態(tài)系統(tǒng)服務(wù)功能土地利用類型林地耕地水域建設(shè)用地氣體調(diào)節(jié)3035.2433.600.00氣候調(diào)節(jié)2341.44771.81398.910.00水源涵養(yǎng)2775.04520.3217673.540.00土壤形成與保護(hù)3382.081266.118.6720.00廢物處理1136.031422.2115765.70.00生物多樣性保護(hù)2827.07615.7122159.330.00食物生產(chǎn)86.72867.286.720.00原材料生產(chǎn)2254.7286.728.670.00娛樂文化1110.028.6723763.650.00總 計(jì)18948.325992.3539865.190.00

漓江流域的生態(tài)系統(tǒng)服務(wù)價(jià)值(ESV)按下式計(jì)算:

(2)

(3)

(4)

其中,ESVi代表j種生態(tài)系統(tǒng)服務(wù)功能價(jià)值總量；Si代表i種土地類型面積；VCij代表i類單位面積土地的j種生態(tài)系統(tǒng)服務(wù)價(jià)值；ESVj代表i種土地類型生態(tài)系統(tǒng)服務(wù)價(jià)值總量；ESV代表生態(tài)系統(tǒng)服務(wù)價(jià)值總量。 結(jié)合漓江流域的生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值(表2)及流域的各類土地覆蓋面積(表3)得到漓江流域的各地類生態(tài)系統(tǒng)服務(wù)價(jià)值(表4)及對應(yīng)的9種生態(tài)系統(tǒng)服務(wù)功能價(jià)值(表5)。

3 結(jié)果與討論分析

3.1 土地利用的時空變化

基于遙感解譯數(shù)據(jù),獲取漓江流域2000、2010、2016三年的土地覆蓋分類圖(圖2)。通過對分類圖進(jìn)行統(tǒng)計(jì)分析,獲得2000、2010、2016年土地利用面積(表3)。3個時期各類型土地面積所占比例如圖3所示。2000—2016年從土地利用類型覆蓋及構(gòu)成比例上看,漓江流域土地利用類型主要為林地,林地面積占流域面積的比例為73.53%～70.02%,比例較小的為建設(shè)用地和水域,分別為1.74%～6.36%和1.27%～1.16%。2000—2016年,漓江流域的土地利用類型發(fā)生了不同程度的變化,林地面積先增加后減少,由2000年的428 819 hm2增加到2010年的429 436 hm2,而后又下降到2016年的408 281 hm2,2000—2010年林地面積增加主要是由于林果種植、退耕還林、石漠化治理等措施的實(shí)施[25]。2000—2016年,耕地及水域面積占流域總面積的比例降低,耕地面積比例從23.46%降到22.46%, 水域面積比例從1.27%降到1.16%, 而建設(shè)用地面積占流域總面積的比例增加,特別是臨桂區(qū),增加明顯。這與城市建設(shè)、旅游用地增加、城鎮(zhèn)化發(fā)展以及2007年桂林城市發(fā)展堅(jiān)持“保護(hù)漓江,發(fā)展臨桂,再造一個新桂林”的戰(zhàn)略部署相對應(yīng)。同時表明，隨著經(jīng)濟(jì)的發(fā)展、政策因素及城市的擴(kuò)張,城市人口不斷增加,城市住房需求不斷增加,建設(shè)用地呈現(xiàn)擴(kuò)張趨勢,而林地及水域面積的減少表明，漓江流域的生態(tài)與環(huán)境在不斷退化。

表3 漓江流域土地覆蓋面積

注: 比例為各地類面積占流域總面積占比。

從變化幅度和轉(zhuǎn)移情況看,2000—2016年建設(shè)用地面積增加, 林地、 耕地和水域的面積出現(xiàn)不同程度的減少, 建設(shè)用地增幅最大, 增加了265.74%, 林地減少4.79%,耕地和水域面積分別減少了4.27%和7.95%。漓江流域的土地利用變化趨勢為林地、耕地和水域向建設(shè)用地轉(zhuǎn)移, 這種轉(zhuǎn)移主要是因?yàn)檎叩淖兓?、城?zhèn)化發(fā)展及其他經(jīng)濟(jì)活動的驅(qū)動。

圖2 2000—2016年漓江流域土地覆蓋分類圖Fig.2 Land cover classification map in Lijiang River basin from 2000 to 2016

圖3 2000—2016年地類面積比例Fig.3 Proportion of land type area from 2000 to 2016

從各區(qū)域的土地利用類型變化來看, 建設(shè)用地增加最為明顯的是興安縣、 臨桂區(qū)、 雁山區(qū)以及桂林市中心城區(qū)(疊彩區(qū)、 秀峰區(qū)、 七星區(qū)、 象山區(qū)), 這與興安縣近幾年的人口快速增長、 桂林市“再造一個新桂林”的臨桂新區(qū)建設(shè)、 雁山區(qū)大學(xué)城建設(shè)以及桂林市“強(qiáng)力疏解改造提升老城, 全面實(shí)施城市交通暢通優(yōu)化工程, 加大老城路橋基礎(chǔ)設(shè)施建設(shè)力度”的戰(zhàn)略部署等有很大的關(guān)系。 城鎮(zhèn)快速發(fā)展、 農(nóng)業(yè)種植結(jié)構(gòu)調(diào)整[26]、 工業(yè)化、 旅游業(yè)發(fā)展[27]以及林業(yè)生態(tài)建設(shè)是土地利用變化的主要驅(qū)動因素。

3.2 漓江流域生態(tài)系統(tǒng)服務(wù)價(jià)值

從各土地利用類型生態(tài)系統(tǒng)服務(wù)價(jià)值及總生態(tài)系統(tǒng)服務(wù)價(jià)值來看(表4), 2000—2016年間, 林地生態(tài)系統(tǒng)服務(wù)價(jià)值呈先小幅增加后下降趨勢, 耕地和水域生態(tài)系統(tǒng)服務(wù)價(jià)值一直呈下降趨勢。 從整個研究期各土地利用類型生態(tài)服務(wù)價(jià)值的構(gòu)成來看, 各地類生態(tài)服務(wù)價(jià)值所占比例有所變化, 不同地類的生態(tài)服務(wù)價(jià)值排序?yàn)? 建設(shè)用地<水域<耕地<林地。 17年間, 林地生態(tài)系統(tǒng)服務(wù)價(jià)值所占比例均超過85%, 遠(yuǎn)遠(yuǎn)大于其他土地類型的生態(tài)系統(tǒng)服務(wù)價(jià)值, 說明了林地是流域內(nèi)整個生態(tài)系統(tǒng)服務(wù)價(jià)值的主體部分, 耕地雖然分布廣, 但生態(tài)系統(tǒng)服務(wù)價(jià)值貢獻(xiàn)低。 2000、 2010、 2016年漓江流域生態(tài)系統(tǒng)服務(wù)總價(jià)值分別為92.48、 92.50和 88億元, 2010—2016年生態(tài)系統(tǒng)服務(wù)總價(jià)值呈下降趨勢。 這表明漓江流域的生態(tài)與環(huán)境在城鎮(zhèn)化及其他產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展的過程當(dāng)中逐漸退化。 土地利用變化及生態(tài)系統(tǒng)服務(wù)價(jià)值的降低主要與人類經(jīng)濟(jì)活動有關(guān)[28]。

從各單項(xiàng)生態(tài)服務(wù)功能價(jià)值變化來看(表5),2000—2016年研究區(qū)各項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值基本都出現(xiàn)先增加后減少的趨勢,總體服務(wù)價(jià)值排序?yàn)?土壤形成與保護(hù)>水源涵養(yǎng)>氣體調(diào)節(jié)>生物多樣性保護(hù)>氣候調(diào)節(jié)>原材料生產(chǎn)>廢物處理>娛樂文化>食物生產(chǎn)。土壤形成與保護(hù)、水源涵養(yǎng)、氣體調(diào)節(jié)貢獻(xiàn)率處于前3位,2000、2010、2016年上述3項(xiàng)服務(wù)功能的合計(jì)貢獻(xiàn)率分別為47.31%、47.33%、47.31%,均接近于總體服務(wù)價(jià)值的一半,主要是由于林地作為漓江流域主要土地類型, 具有水源涵養(yǎng)、 氣體調(diào)節(jié)和較高的土壤形成與保護(hù)功能。生物多樣性保護(hù)、原材料生產(chǎn)和娛樂文化的價(jià)值減少,主要是由于這些生態(tài)服務(wù)價(jià)值貢獻(xiàn)較高的林地和水域面積減少造成的。廢物處理和食物生產(chǎn)價(jià)值減少,主要是由于生產(chǎn)價(jià)值系數(shù)較高的水域及耕地面積減少造成的。

表4 2000—2016年漓江流域各土地利用類型生態(tài)系統(tǒng)服務(wù)價(jià)值

Table 4 Ecosystem service value of each land use type in Lijiang River basin from 2000 to 2016 億元

年份統(tǒng)計(jì)類型林地耕地水域建設(shè)用地總服務(wù)價(jià)值2000價(jià)值81.258.203.03092.48比例/%87.858.873.2801002010價(jià)值81.378.113.02092.50比例/%87.978.773.2601002016價(jià)值77.367.852.79088.00比例/%87.918.923.170100

表5 2000—2016年漓江流域各項(xiàng)生態(tài)系統(tǒng)服務(wù)功能價(jià)值

Table 5 Value of various ecosystem service function in Lijiang River basin from 2000 to 2016 億元

生態(tài)系統(tǒng)服務(wù)功能2000年價(jià)值貢獻(xiàn)率/%2010年價(jià)值貢獻(xiàn)率/%2016年價(jià)值貢獻(xiàn)率/%氣體調(diào)節(jié)13.6114.7213.6214.7212.9614.73氣候調(diào)節(jié)11.1212.0211.1312.0310.6012.05水源涵養(yǎng)13.9115.0413.9215.0513.2115.01土壤形成與保護(hù)16.2317.5516.2417.5615.4617.57廢物處理7.988.637.968.617.578.60生物多樣性保護(hù)13.1314.2013.1314.1912.5014.20食物生產(chǎn)1.581.711.571.701.511.72原材料生產(chǎn)9.7910.589.8010.599.3210.59娛樂文化5.135.555.135.554.875.53合 計(jì)92.4810092.5010088.00100

3.3 漓江流域生態(tài)系統(tǒng)服務(wù)價(jià)值的空間分布

在ArcGIS 10.0軟件中,基于1 km×1 km網(wǎng)格分別計(jì)算2000、2010、2016年漓江流域生態(tài)系統(tǒng)服務(wù)價(jià)值的空間差異,將生態(tài)系統(tǒng)服務(wù)價(jià)值分為極低(<45萬元/hm2)、 低(45～90萬元/hm2)、 中(90～135萬元/hm2)、 高(135～180萬元/hm2)和極高(>180萬元/hm2)5個級別(圖4)[29]。流域生態(tài)服務(wù)價(jià)值極低的區(qū)域主要分布在興安縣的東南方向、臨桂區(qū)、雁山區(qū)以及桂林市中心城區(qū)(疊彩區(qū)、 秀峰區(qū)、 七星區(qū)、 象山區(qū)),且極低區(qū)域有擴(kuò)大的趨勢,特別是2010—2016年,擴(kuò)大范圍廣。生態(tài)系統(tǒng)服務(wù)價(jià)值的極低區(qū)域與建設(shè)用地及耕地空間分布相吻合,高和極高區(qū)域與林地及水域空間分布相吻合。從時空分布分析,漓江流域的北部生態(tài)系統(tǒng)服務(wù)價(jià)值逐年擴(kuò)大,中西部的生態(tài)系統(tǒng)服務(wù)價(jià)值逐年減少,低生態(tài)系統(tǒng)服務(wù)價(jià)值區(qū)域逐年向極低生態(tài)系統(tǒng)服務(wù)價(jià)值轉(zhuǎn)移；東南部的部分區(qū)域(陽朔及靈川的部分區(qū)域)及高生態(tài)系統(tǒng)服務(wù)價(jià)值區(qū)域有擴(kuò)大趨勢；中南部及主要水域區(qū)域出現(xiàn)極高生態(tài)面積減少的趨勢。生態(tài)系統(tǒng)服務(wù)價(jià)值的空間分布變化表明,漓江流域的源頭林地保護(hù)初見成效,隨著退耕還林政策的逐步實(shí)施,有部分區(qū)域的生態(tài)系統(tǒng)逐步恢復(fù),但是城鎮(zhèn)化的發(fā)展及人口的增多,造成了建設(shè)用地的增多及部分水域面積的減少,從而導(dǎo)致生態(tài)系統(tǒng)服務(wù)價(jià)值的降低,流域內(nèi)生態(tài)系統(tǒng)的破壞超出了對生態(tài)系統(tǒng)的保護(hù)。

4 結(jié) 論

(1)基于3S技術(shù),研究漓江流域土地利用時空變化對生態(tài)系統(tǒng)服務(wù)價(jià)值的影響,為漓江流域土地利用發(fā)展提供了參考依據(jù)。通過對研究區(qū)域2000—2016年間土地利用圖進(jìn)行變化分析發(fā)現(xiàn),流域土地利用發(fā)生了較大變化,自2007年桂林市“再造一個新桂林”的戰(zhàn)略部署后, 建設(shè)用地增加顯著。 建設(shè)用地增加主要發(fā)生在臨桂區(qū)、桂林市中心城區(qū)及興安縣等。

圖4 2000—2016年漓江流域生態(tài)系統(tǒng)服務(wù)價(jià)值時空分布Fig.4 Temporal and spatial distribution of ecosystem service value in Lijiang River basin from 2000 to 2016

(2)根據(jù)不同的地域特點(diǎn),對流域的單位面積耕地食物生產(chǎn)生態(tài)價(jià)值進(jìn)行修訂,結(jié)合流域土地類型面積計(jì)算生態(tài)系統(tǒng)服務(wù)價(jià)值,結(jié)果發(fā)現(xiàn)生態(tài)系統(tǒng)服務(wù)總價(jià)呈先微弱增加后下降的趨勢,土壤形成與保護(hù)、水源涵養(yǎng)、氣體調(diào)節(jié)這3項(xiàng)生態(tài)系統(tǒng)服務(wù)功能占主要比例。2000—2016年間生態(tài)系統(tǒng)服務(wù)價(jià)值(ESV)的變化表明，研究區(qū)內(nèi)生態(tài)與環(huán)境正在逐步退化,流域土地利用變化,特別是城市擴(kuò)張對生態(tài)系統(tǒng)服務(wù)產(chǎn)生重大影響。

(3)生態(tài)系統(tǒng)服務(wù)價(jià)值極低區(qū)域出現(xiàn)上升趨勢,并主要集中在興安縣東南方向及桂林市中心城區(qū)和臨桂新區(qū)。中值及低值區(qū)域逐漸被極低區(qū)域所替代,部分極高值區(qū)域出現(xiàn)減少趨勢。研究期間,只有北部及東南部的部分區(qū)域出現(xiàn)生態(tài)系統(tǒng)服務(wù)增加的趨勢。隨著發(fā)展,流域的土地利用類型結(jié)構(gòu)發(fā)生變化,建議相關(guān)部門在土地利用的過程中考慮人為因素對生態(tài)系統(tǒng)的影響,在經(jīng)濟(jì)發(fā)展的同時,要做好生態(tài)保護(hù),保護(hù)水域、林地、耕地等生態(tài)價(jià)值系數(shù)高的土地利用類型,以提升流域生態(tài)系統(tǒng)服務(wù)價(jià)值,做到經(jīng)濟(jì)、社會及生態(tài)三者的協(xié)調(diào)可持續(xù)發(fā)展。