蔡耀欣 江怡萱 謝祥財(cái)

(福建農(nóng)林大學(xué) 風(fēng)景園林與藝術(shù)學(xué)院，福建 福州 350002)

生態(tài)系統(tǒng)服務(wù)(Ecosystem Services，ES)是指生態(tài)系統(tǒng)與生態(tài)過(guò)程中所產(chǎn)生的，保證人類(lèi)生存的生命支撐產(chǎn)品和服務(wù)[1-2]。土地利用是人與自然互動(dòng)最為密切的部分，土地利用變化(Land Use Change，LUC)是人類(lèi)作用于自然環(huán)境的重要方式，它能夠影響生態(tài)系統(tǒng)產(chǎn)品與服務(wù)提供能力[3-4]。對(duì)土地利用變化進(jìn)行評(píng)估，是評(píng)價(jià)生態(tài)系統(tǒng)服務(wù)價(jià)值(Ecosystem Services Value，ESV)的基礎(chǔ)，也是科學(xué)規(guī)劃生態(tài)環(huán)境建設(shè)的重要前提[5]。

自生態(tài)系統(tǒng)服務(wù)這一概念被提出以來(lái)，在生態(tài)系統(tǒng)服務(wù)概念及其分類(lèi)、價(jià)值量化評(píng)估和應(yīng)用等方面取得了一定的進(jìn)展。1997年Costanza等構(gòu)建了“全球生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量表”及“全球生態(tài)服務(wù)價(jià)值分析模型”，奠定了生態(tài)系統(tǒng)服務(wù)價(jià)值研究的理論基礎(chǔ)[6]。2003年謝高地等根據(jù)Costanza的方法構(gòu)建了符合中國(guó)國(guó)情的當(dāng)量因子表，并確立了中國(guó)單位面積價(jià)值當(dāng)量及計(jì)算方法[7-8]。國(guó)內(nèi)眾多學(xué)者據(jù)此當(dāng)量表展開(kāi)相關(guān)研究，如張敏等對(duì)各類(lèi)生態(tài)服務(wù)價(jià)值間的權(quán)衡與協(xié)同關(guān)系進(jìn)行研究[9]，劉玉卿及楊延成等從景觀格局演變角度分析對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響[10-11]等。2015年，謝高地等對(duì)當(dāng)量因子表進(jìn)行修訂和細(xì)分[12]，使更多學(xué)者能夠準(zhǔn)確地通過(guò)土地利用變化評(píng)估我國(guó)生態(tài)系統(tǒng)服務(wù)功能的價(jià)值，如鄭舒瓊等以行政區(qū)域作為基本單元，研究土地利用對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空差異的影響[13]。相關(guān)研究均表明土地利用變化不僅影響著生態(tài)系統(tǒng)服務(wù)價(jià)值綜合水平，且對(duì)不同類(lèi)別的生態(tài)系統(tǒng)服務(wù)影響程度也不同。

晉江流域是福建省經(jīng)濟(jì)快速增長(zhǎng)的區(qū)域，是泉州市人口最密集、城鎮(zhèn)化水平最高的地區(qū)[14]。目前對(duì)于晉江流域的研究，多集中于水文調(diào)節(jié)[15]、水資源管理[16]以及水環(huán)境治理[17]等方面，加強(qiáng)對(duì)晉江流域生態(tài)系統(tǒng)服務(wù)價(jià)值的研究對(duì)流域內(nèi)生態(tài)文明建設(shè)及社會(huì)經(jīng)濟(jì)可持續(xù)發(fā)展具有重要理論價(jià)值和現(xiàn)實(shí)意義。本文通過(guò)地理信息技術(shù)對(duì)晉江流域2000—2020年間土地利用狀況進(jìn)行解譯，采用自相關(guān)分析、熱點(diǎn)分析等方法研究土地利用變化對(duì)流域生態(tài)系統(tǒng)服務(wù)價(jià)值的影響及二者間的相互關(guān)系，旨在為該區(qū)域土地資源開(kāi)發(fā)利用和生態(tài)系統(tǒng)可持續(xù)發(fā)展與管理提供科學(xué)依據(jù)。

1 研究區(qū)域概況

晉江流域(見(jiàn)圖1)位于福建省泉州市。晉江是泉州市第一大河流，也是福建省內(nèi)僅次于閩江、九龍江的第三大江，源頭至河口全長(zhǎng)182 km，總體自西北向東南流，流經(jīng)永春縣、安溪縣、南安市、鯉城區(qū)、晉江市。晉江流域總面積5 629 km2，流域內(nèi)年平均氣溫17-21℃，季風(fēng)氣候顯著，降水豐沛[18]。晉江上游有東溪和西溪兩大支流。東溪流域面積1 917 km2，河長(zhǎng)120 km，東溪中游有山美水庫(kù)，是流域內(nèi)重要的水利樞紐工程。西溪流域面積3 101 km2，河長(zhǎng)145 km。兩支流于南安市豐州鎮(zhèn)井兜村雙溪口匯合，并于豐澤區(qū)潯浦入海。

圖1 研究區(qū)示意圖

2 數(shù)據(jù)來(lái)源與處理

本研究采用的晉江流域2000、2010和2020年的土地利用數(shù)據(jù)來(lái)源于地理空間數(shù)據(jù)云平臺(tái)(http://www.gscloud.cn)，空間分辨率為30m×30m，數(shù)據(jù)處理平臺(tái)為ArcGIS 10.6。晉江流域單個(gè)當(dāng)量因子的經(jīng)濟(jì)價(jià)值量計(jì)算主要來(lái)自2010—2020年的《福建省統(tǒng)計(jì)年鑒》《泉州市統(tǒng)計(jì)年鑒》等，糧食價(jià)格數(shù)據(jù)來(lái)源于福建糧油局的公開(kāi)數(shù)據(jù)。

3 研究方法

3.1 土地利用變化計(jì)算

土地利用轉(zhuǎn)移矩陣用于描述研究初期和末期區(qū)域內(nèi)不同土地利用類(lèi)型間的轉(zhuǎn)化情況，有助于發(fā)現(xiàn)土地使用種類(lèi)之間的轉(zhuǎn)化方向[19]，表達(dá)式為:

(1)

式中：S為各地類(lèi)面積，n為土地利用類(lèi)型數(shù)，i、j分別為研究初期和末期的土地利用類(lèi)型。

3.2 生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估方法

結(jié)合實(shí)際情況，對(duì)謝高地等2015年最新改進(jìn)的基礎(chǔ)當(dāng)量表進(jìn)行調(diào)整，得出晉江流域單位面積ESV基礎(chǔ)當(dāng)量表(見(jiàn)表1)。根據(jù)統(tǒng)計(jì)年鑒中2010—2020年福建省糧食作物平均單產(chǎn)為5 756.98kg/hm2、糧食平均價(jià)格為2.88元/kg，依據(jù)相關(guān)文獻(xiàn)，1個(gè)標(biāo)準(zhǔn)當(dāng)量因子的ESV大約等于糧食單位產(chǎn)量?jī)r(jià)值的1/7[20]，計(jì)算得出研究區(qū)的單個(gè)當(dāng)量因子的經(jīng)濟(jì)價(jià)值量為2 371.36元/hm2。最終根據(jù)修正后的基礎(chǔ)當(dāng)量和各類(lèi)型土地利用的面積計(jì)算晉江流域的ESV(見(jiàn)表2)。計(jì)算公式為：

表1 晉江流域單位面積ESV基礎(chǔ)當(dāng)量

表2 晉江流域ESV系數(shù)表 單位：元/hm2·a

(2)

ESVj=∑(Ai×VCij)

(3)

式中：ESV為生態(tài)系統(tǒng)服務(wù)的總價(jià)值，ESVj為 研究區(qū)第j項(xiàng)生態(tài)系統(tǒng)服務(wù)功能價(jià)值；Ai為研究區(qū)第i種類(lèi)型的土地面積，VCi為第i種類(lèi)型土地的單位生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)，VCij為第i種土地利用類(lèi)型的第j項(xiàng)生態(tài)系統(tǒng)服務(wù)功能價(jià)值系數(shù)[12]。

3.3 ESV網(wǎng)格計(jì)算

根據(jù)相關(guān)公式[21]，將研究區(qū)劃分為23 078個(gè)互不重疊的正方形格網(wǎng)，并進(jìn)行疊加分析計(jì)算每個(gè)格網(wǎng)的ESV，得出研究區(qū)2000、2010和2020年ESV值，生成ESV空間分布圖(見(jiàn)圖6)，并采用自然間斷點(diǎn)分級(jí)法將各格網(wǎng)ESV值劃分為5個(gè)等級(jí)，等級(jí)5對(duì)應(yīng)最高值。

3.4 生態(tài)系統(tǒng)服務(wù)價(jià)值變化特征分析

3.4.1 熱點(diǎn)分析

熱點(diǎn)分析可用于研究ESV空間變化是否具有高值聚集和低值聚集的現(xiàn)象，以及確定高值區(qū)和低值區(qū)在空間上發(fā)生聚類(lèi)的位置，公式如下[22]:

(4)

(5)

(6)

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

生態(tài)貢獻(xiàn)率表示某一時(shí)間段內(nèi)不同土地利用變化所產(chǎn)生的生態(tài)服務(wù)價(jià)值變化量對(duì)區(qū)域總服務(wù)價(jià)值變化量的影響大小和方向，可用來(lái)揭示影響區(qū)域生態(tài)服務(wù)價(jià)值變化的主要貢獻(xiàn)因子[23]。

(7)

式中：SKT為K類(lèi)生態(tài)系統(tǒng)在時(shí)間段T的生態(tài)服務(wù)貢獻(xiàn)率，△ESVKT為K類(lèi)生態(tài)系統(tǒng)在時(shí)間段T內(nèi)的生態(tài)服務(wù)價(jià)值變化量。

4 結(jié)果與分析

4.1 土地利用變化分析

4.1.1 數(shù)量變化分析

通過(guò)對(duì)各用地類(lèi)型的面積進(jìn)行統(tǒng)計(jì)，得到晉江流域土地利用面積以及變化率(見(jiàn)表3)。從占比看，流域內(nèi)林地面積最大，占整個(gè)研究區(qū)域的44.75%；其次是耕地，占比29.31%。從變化量看，耕地和建設(shè)用地的變化最為突出，其中建設(shè)用地大幅增加，面積擴(kuò)大21 536.1hm2，增幅為48.64%；耕地面積大幅減少，面積縮減18 500.6hm2，減幅為10.12%；水體面積增幅僅次于建設(shè)用地，為22.26%。

表3 2000—2020年晉江流域土地利用面積和變化率 單位：hm2、%

4.1.2 空間變化分析

由圖2可以看出，晉江流域林地分布范圍最廣，主要聚集在北部的德化縣、永春縣以及南部的安溪縣。2000—2020年，晉江中下游河岸兩側(cè)的建設(shè)用地面積不斷增加，晉江入?？谔幗ㄔO(shè)用地持續(xù)擴(kuò)張，這表明晉江對(duì)于推進(jìn)城市化進(jìn)程以及推動(dòng)經(jīng)濟(jì)發(fā)展有著重要的作用。

圖2 晉江流域土地利用分布

4.1.3 土地利用轉(zhuǎn)移分析

由表4可以看出，晉江流域耕地的轉(zhuǎn)出面積最多，為39 522.60 hm2，主要轉(zhuǎn)出目標(biāo)是建設(shè)用地，轉(zhuǎn)出22 970.34 hm2；其次為林地，轉(zhuǎn)出9 165.06 hm2。林地的轉(zhuǎn)入面積最多，為32 374.71 hm2，其中主要轉(zhuǎn)入源為草地，轉(zhuǎn)入面積為21 781.62 hm2。

表4 2000—2020年晉江流域土地利用和生態(tài)系統(tǒng)服務(wù)價(jià)值轉(zhuǎn)移矩陣 單位：hm2、萬(wàn)元

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

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

由表5以及圖3可知，2000—2020年晉江流域ESV呈現(xiàn)先減后增的趨勢(shì)，前10年共減少了1.95億元，后10年共增加了1.59億元，20年間共減少了0.37億元，增幅為-0.18%。從土地利用類(lèi)型來(lái)看，林地是構(gòu)成晉江流域生態(tài)系統(tǒng)服務(wù)價(jià)值的主體，在2020年達(dá)到137.41億元，占總ESV的69.46%；其次是草地和水體，分別占13.64%和9.11%。從土地利用ESV20年的變化來(lái)看，只有水體的ESV呈上升趨勢(shì)，增加了3.35億元，增長(zhǎng)率為22.86%；耕地的ESV減少最多，減少了1.73億元。

表5 2000—2020年晉江流域土地利用ESV變化情況 單位：億元、%

圖3 2000—2020年晉江流域生態(tài)系統(tǒng)服務(wù)價(jià)值總量變化

由圖4可知，晉江流域生態(tài)貢獻(xiàn)率總和為-5.16%，導(dǎo)致ESV減少的地類(lèi)有耕地、林地、草地、濕地、未利用地，生態(tài)貢獻(xiàn)率之和為-52.58%，耕地是流域主要的負(fù)向貢獻(xiàn)因子，生態(tài)貢獻(xiàn)率為-24.50%；水體是晉江流域ESV唯一的正向貢獻(xiàn)因子，生態(tài)貢獻(xiàn)率為47.42%。

圖4 2000—2020年晉江流域生態(tài)貢獻(xiàn)率

從不同服務(wù)類(lèi)型看(見(jiàn)表6、圖5)，氣候調(diào)節(jié)和水文調(diào)節(jié)占比最高，分別占25.84%和23.57%，原因在于提供氣候調(diào)節(jié)最多的林地面積最大，以及水體提供的水文調(diào)節(jié)基礎(chǔ)當(dāng)量遠(yuǎn)遠(yuǎn)大于其他地類(lèi)。所有服務(wù)類(lèi)型中，僅有水資源供給與水文調(diào)節(jié)的ESV在20年間總體處于增加的趨勢(shì)，其中水資源供給增加了18.60%，水文調(diào)節(jié)增加了1.82%。

表6 2000—2020年晉江流域生態(tài)系統(tǒng)服務(wù)ESV變化情況 單位：億元、%

圖5 2000—2020年晉江流域生態(tài)系統(tǒng)服務(wù)功能價(jià)值

4.2.2 生態(tài)系統(tǒng)服務(wù)價(jià)值空間變化分析

由圖6可以看出，ESV高值區(qū)(等級(jí)4、等級(jí)5)主要分布于晉江中下游水域、山美水庫(kù)區(qū)域以及安溪、永春、德化山區(qū)等林地密集區(qū)；中值區(qū)(等級(jí)3)主要分布于中部及北部草地覆蓋率高的丘陵地帶；低值區(qū)(等級(jí)1、等級(jí)2)主要分布于晉江沿岸城鎮(zhèn)化水平較高的城市和鄉(xiāng)村聚落。由表7可知，晉江流域2000—2020年的ESV全局Moran’s I值均大于0、P值均小于0.001，表明流域內(nèi)ESV的空間分布具有明顯的空間正相關(guān)性，在空間上呈現(xiàn)顯著聚集性分布。由圖7可知，高聚集區(qū)主要分布在晉江中下游和山美水庫(kù)區(qū)域，以及北部與西南部的山區(qū)林地；低聚集區(qū)主要聚集在晉江下游的南安市、豐澤區(qū)、鯉城區(qū)等城鎮(zhèn)化水平較高、人口較為密集的區(qū)域。

表7 2000—2020年晉江流域ESV全局空間自相關(guān)顯著性檢驗(yàn)

圖6 晉江流域ESV空間分布

圖7 晉江流域生態(tài)系統(tǒng)服務(wù)價(jià)值LISA分布

由圖8可知，2000—2020年，晉江流域大部分區(qū)域的ESV處于變化非顯著的狀態(tài)，前10年損失冷點(diǎn)區(qū)多于增值熱點(diǎn)區(qū)，主要分布在晉江中下游的城鎮(zhèn)或聚落周?chē)?；?0年增值熱點(diǎn)區(qū)多于損失冷點(diǎn)區(qū)，且明顯地分布于西溪支流與胡洋溪至山美水庫(kù)的區(qū)域，說(shuō)明了自“十二五”規(guī)劃以來(lái)，泉州市實(shí)施的重點(diǎn)江河流域生態(tài)修復(fù)工程取得了顯著的成效。

圖8 晉江流域ESV變化熱點(diǎn)區(qū)域空間分布

5 討論與結(jié)論

5.1 討論

本文采用的謝高地等改進(jìn)后的當(dāng)量因子法，數(shù)據(jù)需求少，可行性高且評(píng)價(jià)內(nèi)容較為全面。但由于數(shù)據(jù)獲取的方法局限，選用2010—2020年福建省糧食產(chǎn)量與糧食價(jià)格數(shù)據(jù)也在一定程度上影響了價(jià)值評(píng)估的準(zhǔn)確性。其次，本文主要以土地利用變化這一要素作為研究對(duì)象，未考慮其他方面的影響因素，如人口、經(jīng)濟(jì)、生物等，后期研究中應(yīng)更多考慮整合諸多因素。

通過(guò)前文分析，2000—2020年晉江流域生態(tài)系統(tǒng)服務(wù)總價(jià)值呈先下降后上升的趨勢(shì)；在用地類(lèi)型中，水體的ESV呈增長(zhǎng)趨勢(shì)，其余地類(lèi)均呈減少趨勢(shì)；在服務(wù)類(lèi)型中，水資源供給與水文調(diào)節(jié)呈增長(zhǎng)趨勢(shì)，其余服務(wù)類(lèi)型均呈減少趨勢(shì)，由此可見(jiàn)水體及其所提供的生態(tài)系統(tǒng)服務(wù)是晉江流域生態(tài)系統(tǒng)服務(wù)價(jià)值的重要正向驅(qū)動(dòng)因子。大量耕地和草地變成了建設(shè)用地，是導(dǎo)致晉江河岸周?chē)ǔ蓞^(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值降低的主要原因；河岸周?chē)剞D(zhuǎn)變?yōu)樗w和林地，是導(dǎo)致晉江流域庫(kù)區(qū)以及下游周?chē)鷳B(tài)系統(tǒng)服務(wù)價(jià)值升高的主要原因。

從空間上看，晉江流域生態(tài)系統(tǒng)服務(wù)價(jià)值的改變都集中于晉江的中下游部分，經(jīng)分析原因有二：(1)晉江中下游地形多為起伏較小的平原，城鎮(zhèn)化水平高，村落聚集、人口密集，是城市發(fā)展的核心區(qū)域，建設(shè)用地的不斷擴(kuò)張導(dǎo)致了2000—2010年ESV的急劇減少；位于晉江上游的地形多為丘陵，城鎮(zhèn)化水平較低，聚落分散、人口稀少，常年保持著耕作的傳統(tǒng)生產(chǎn)方式，用地性質(zhì)改變程度不高，因此對(duì)于ESV的影響較不顯著。(2)由于政策的影響，城市對(duì)流域水生態(tài)環(huán)境的保護(hù)逐見(jiàn)成效——從“十二五”規(guī)劃時(shí)期開(kāi)始，泉州開(kāi)始實(shí)施重點(diǎn)江河流域生態(tài)修復(fù)工程，積極推進(jìn)河岸生態(tài)緩沖帶修復(fù)，促進(jìn)生態(tài)良性循環(huán)，促使山美水庫(kù)區(qū)域以及晉江西溪流域的ESV在2010—2020年得到了顯著提高。

5.2 結(jié)論

(1)土地利用變化層面，2000—2020年晉江流域耕地面積減少最多，共減少了18 500.6 hm2，建設(shè)用地和水體面積不斷增加，分別增加21 536.1 hm2和1 126.26 hm2。從土地利用轉(zhuǎn)移上看，耕地轉(zhuǎn)換為建設(shè)用地的面積最大，林地與草地之間的轉(zhuǎn)換次之。

(2)生態(tài)系統(tǒng)服務(wù)價(jià)值層面，2000—2020年晉江流域總ESV整體呈下降趨勢(shì)，細(xì)分之后前10年為下降趨勢(shì)，后10年為上升趨勢(shì)。林地是晉江流域ESV最主要的組成部分，且20年來(lái)占比保持在70%左右。在空間分布上，ESV高值區(qū)主要分布在晉江上游森林覆蓋率高的山區(qū)以及晉江水體區(qū)域，低值區(qū)主要分布在晉江中下游的城鎮(zhèn)聚集區(qū)。從時(shí)間變化來(lái)看，晉江干流以及西溪流域的ESV上升明顯，形成了增值熱點(diǎn)區(qū)；而城鎮(zhèn)面積的不斷擴(kuò)張與人口的不斷聚集，造成了晉江下游周?chē)纬闪税邏K狀的損失冷點(diǎn)區(qū)。

(3)土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值影響層面，土地利用的變化影響了生態(tài)要素在晉江流域的分布，進(jìn)而影響生態(tài)系統(tǒng)服務(wù)的分布和構(gòu)成。從數(shù)量變化上看，水體是唯一對(duì)晉江流域ESV提供正向貢獻(xiàn)的用地類(lèi)型，大量耕地轉(zhuǎn)換成為水體促使水文調(diào)節(jié)功能大幅提升，從而造成ESV在2010年后顯著增加。從空間變化上看，研究期間晉江水域附近的耕地轉(zhuǎn)換為水體形成了ESV增值區(qū);建成區(qū)附近的耕地以及水域周?chē)莸剞D(zhuǎn)換為建設(shè)用地形成ESV損失區(qū)。

研究表明，晉江流域水體生態(tài)效益的增加是促進(jìn)流域ESV升高的主要原因。僅占流域總面積1.08%的水體，主導(dǎo)著流域生態(tài)系統(tǒng)服務(wù)價(jià)值的變化，表明水體應(yīng)該成為晉江流域環(huán)境保護(hù)和生態(tài)建設(shè)的重心，因此，提高區(qū)域內(nèi)河流的水生態(tài)環(huán)境將有益于提升流域的總體生態(tài)環(huán)境。今后應(yīng)加強(qiáng)水源涵養(yǎng)林建設(shè)與保護(hù)和河岸生態(tài)治理，選取本地優(yōu)勢(shì)樹(shù)種在河道兩側(cè)建設(shè)植被緩沖帶和隔離帶；加快流域及山美水庫(kù)等水源涵養(yǎng)林和水土保持林修復(fù)；對(duì)河道水利設(shè)施建設(shè)等造成的河道破壞地區(qū)，開(kāi)展植樹(shù)造林，恢復(fù)河道自然生境；推進(jìn)泉州灣河口濕地省級(jí)自然保護(hù)區(qū)“退養(yǎng)還灘”。