昆崳山國(guó)家級(jí)自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值演變

張文馨, 范小莉, 王 強(qiáng), 房 用, 時(shí) 良, 梁 玉

(1.山東省林業(yè)科學(xué)研究院, 濟(jì)南 250014; 2.山東省林業(yè)外資與工程項(xiàng)目管理站,濟(jì)南 250014; 3.煙臺(tái)市昆崳山林場(chǎng), 山東 煙臺(tái)264112)

生態(tài)系統(tǒng)服務(wù)是指生態(tài)系統(tǒng)和生態(tài)過程所形成并維持人類生存所必需的環(huán)境條件及獲得的福祉[1-3]。生態(tài)系統(tǒng)服務(wù)價(jià)值(ecosystem service value，ESV)核算是表征一個(gè)地區(qū)生態(tài)環(huán)境變化的重要指標(biāo)，是實(shí)現(xiàn)生態(tài)安全，評(píng)估生態(tài)保護(hù)成效的重要手段[3-5]，基于此，ESV在相關(guān)研究領(lǐng)域成為國(guó)內(nèi)外研究的重點(diǎn)和熱點(diǎn)[6]。1997年，Costanza首次對(duì)全球生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行了核算[7]，推進(jìn)了ESV研究的快速發(fā)展。作為生態(tài)系統(tǒng)服務(wù)價(jià)值主要估算方法之一，基于土地利用/覆被數(shù)據(jù)的面積當(dāng)量法在ESV研究中發(fā)揮了重要作用[8]，尤其是謝高地對(duì)中國(guó)多種陸地生態(tài)系統(tǒng)單位面積ESV當(dāng)量表修訂后[9]，該方法在核算森林、濕地等ESV及分析土地利用對(duì)ESV影響等多個(gè)領(lǐng)域都取得了顯著進(jìn)展，黃木易等對(duì)1970—2015年大別山區(qū)生態(tài)服務(wù)價(jià)值進(jìn)行測(cè)算[10]，鄧楚雄等分析了洞庭湖區(qū)土地生態(tài)系統(tǒng)的服務(wù)價(jià)值時(shí)空變化[11]，王航等分析了土地利用及景觀格局演變對(duì)淮河上游生態(tài)系統(tǒng)服務(wù)價(jià)值的影響[12]，薛冰等評(píng)估了貴州省森林生態(tài)系統(tǒng)服務(wù)功能[13]，婁佩卿等分析了1998—2018年京津冀土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響[14]，這些研究成果為我國(guó)國(guó)土空間管制[15]、重點(diǎn)生態(tài)功能區(qū)管理和國(guó)家生態(tài)安全格局的構(gòu)建[16-17]及自然保護(hù)地生態(tài)補(bǔ)償?shù)萚18]方面提供了有益參考。

昆崳山國(guó)家級(jí)自然保護(hù)區(qū)成立于2008年，是山東省第一個(gè)森林生態(tài)類型的國(guó)家級(jí)自然保護(hù)區(qū)[19]，也是山東省重要生態(tài)功能區(qū)，在保護(hù)生物多樣性、涵養(yǎng)水源、水土保持等方面發(fā)揮了重要作用。昆崳山國(guó)家級(jí)自然保護(hù)區(qū)地處我國(guó)東部沿海的經(jīng)濟(jì)發(fā)達(dá)地區(qū)，經(jīng)濟(jì)的快速發(fā)展和頻繁的人類活動(dòng)對(duì)該區(qū)域的生態(tài)系統(tǒng)造成了較大威脅，也使得土地利用類型發(fā)生了變化。土地作為經(jīng)濟(jì)發(fā)展和人類活動(dòng)的基礎(chǔ)，其類型的變化對(duì)生態(tài)系統(tǒng)服務(wù)功能及價(jià)值產(chǎn)生明顯的影響[11,20-21]。目前針對(duì)昆崳山區(qū)域的生態(tài)系統(tǒng)價(jià)值評(píng)估，主要采用影子工程法、機(jī)會(huì)成本法、市場(chǎng)價(jià)值法及利用森林資源調(diào)查及長(zhǎng)期監(jiān)測(cè)數(shù)據(jù)進(jìn)行核算[22-23]，在區(qū)域尺度上基于土地利用/覆被數(shù)據(jù)的面積當(dāng)量法計(jì)算昆崳山區(qū)域的ESV尚未見報(bào)道。本文研究2008—2018年不同時(shí)期土地利用類型變化背景下保護(hù)區(qū)ESV的變化，掌握自然保護(hù)區(qū)經(jīng)過十年的發(fā)展后其生態(tài)系統(tǒng)服務(wù)功能價(jià)值的動(dòng)態(tài)變化，為自然保護(hù)地整合、生態(tài)效益補(bǔ)償?shù)忍峁┲匾獏⒖?，從而進(jìn)一步明確保護(hù)區(qū)的生態(tài)功能定位，促進(jìn)昆崳山國(guó)家級(jí)自然保護(hù)區(qū)的可持續(xù)發(fā)展。

1 材料與方法

1.1 研究區(qū)概況

昆崳山國(guó)家級(jí)自然保護(hù)區(qū)位于山東半島東部，跨牟平、文登兩區(qū)(市)，地理坐標(biāo)為東經(jīng)121°37′0″—121°51′0″，北緯37°12′20″—37°18′50″。保護(hù)區(qū)氣候?qū)倥瘻貛Ъ撅L(fēng)型大陸性氣候，受太平洋暖濕氣流和西伯利亞干冷氣流控制，具有四季分明，季風(fēng)顯著，雨熱同期，空氣濕潤(rùn)，溫差較小，光照充足等特點(diǎn)。年均氣溫11.9℃，年均降水984.4 mm，年均蒸發(fā)量1 923.4 mm，年均相對(duì)濕度71%，無(wú)霜期200 d左右。昆崳山國(guó)家級(jí)自然保護(hù)區(qū)森林覆蓋率77.07%，屬于典型的森林生態(tài)系統(tǒng)，是我國(guó)天然赤松林面積最大、保護(hù)最為完整的區(qū)域，在世界上也具有重要代表性[22-23]。

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

研究中2008年、2013年和 2018年不同時(shí)期昆崳山的遙感影像分別是由Landsat TM數(shù)據(jù)和 Landsat OLI數(shù)據(jù)組成，分辨率30 m。其中，2008年遙感影像數(shù)據(jù)采用的是陸地衛(wèi)星 TM5數(shù)據(jù)，2013年遙感影像數(shù)據(jù)以及2018年的遙感影像數(shù)據(jù)則采用陸地衛(wèi)星 Landsat OLI數(shù)據(jù)。研究所利用數(shù)據(jù)均從地理空間數(shù)據(jù)云平臺(tái)和中國(guó)科學(xué)院遙感和數(shù)字研究所數(shù)據(jù)共享平臺(tái)獲取。

根據(jù)昆崳山國(guó)家級(jí)自然保護(hù)區(qū)土地利用現(xiàn)狀和研究目標(biāo)，參照土地利用現(xiàn)狀分類標(biāo)準(zhǔn)(GBT21010—2007)，將研究區(qū)的土地利用分為 5類，分別為林地、耕地、建設(shè)用地、水域、裸地及其他用地。

1.3 土地利用動(dòng)態(tài)度分析

土地利用動(dòng)態(tài)度是評(píng)價(jià)一定時(shí)期內(nèi)不同土地利用類型變化速度的指標(biāo)，包括單一土地利用類型動(dòng)態(tài)度和綜合土地利用動(dòng)態(tài)度[24-26]。

單一土地利用類型動(dòng)態(tài)度：

(1)

式中:K為研究期某種土地利用類型動(dòng)態(tài)度(%);Ua和Ub分別為研究期末和期初某種土地利用類型面積(hm2);T為研究時(shí)長(zhǎng)(a)。

綜合土地利用動(dòng)態(tài)度C(%)：

(2)

式中:Ui為研究期初第i類土地利用類型面積(hm2);Ui,j為研究期第i類土地利用類型轉(zhuǎn)為非i類土地利用類型面積的絕對(duì)值(hm2);T為研究時(shí)長(zhǎng)(a);n為研究區(qū)的土地利用類型總數(shù)。

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

本研究的生態(tài)系統(tǒng)服務(wù)價(jià)值核算運(yùn)用當(dāng)量因子法，參照謝高地等人研究結(jié)果，將昆崳山生態(tài)系統(tǒng)服務(wù)類型分為供給服務(wù)、調(diào)節(jié)服務(wù)、支持服務(wù)和文化服務(wù)4個(gè)一級(jí)類型，進(jìn)一步細(xì)化為食物生產(chǎn)、原料生產(chǎn)、水資源供給、氣體調(diào)節(jié)、氣候調(diào)節(jié)、水文調(diào)節(jié)、凈化環(huán)境、土壤保持、維持養(yǎng)分循環(huán)、生物多樣性和美學(xué)景觀11個(gè)二級(jí)類型[9,27]。

根據(jù)謝高地等人制定的價(jià)值當(dāng)量表以及對(duì)價(jià)值當(dāng)量的計(jì)算，結(jié)合我國(guó)不同省份農(nóng)田生態(tài)系統(tǒng)生物量因子[9,28]，根據(jù)昆崳山實(shí)際，對(duì)自然保護(hù)區(qū)各生態(tài)系統(tǒng)的服務(wù)價(jià)值系數(shù)進(jìn)行修正：林地對(duì)應(yīng)針闊混交林，耕地對(duì)應(yīng)旱地，水域?yàn)闈竦睾退档钠骄?，裸地及其他?duì)應(yīng)裸地，建設(shè)用地不考慮生態(tài)系統(tǒng)服務(wù)價(jià)值，從而得出昆崳山單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量(表1)和生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)(表2)。

表1 昆崳山國(guó)家自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子

表2 昆崳山國(guó)家級(jí)自然保護(hù)區(qū)單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù) 元/(hm2·a)

生態(tài)系統(tǒng)服務(wù)價(jià)值計(jì)算引用Costanza等[7]計(jì)算公式為：

(3)

(4)

式中:ESV為昆崳山生態(tài)系統(tǒng)服務(wù)價(jià)值總量(元);Ai為第i類土地利用類型面積(hm2);VCi為第i類土地利用類型對(duì)應(yīng)的價(jià)值系數(shù)(元/hm2·a);n為研究區(qū)的土地利用類型總數(shù)。ESVf為昆崳山生態(tài)系統(tǒng)各類服務(wù)價(jià)值(元);VCfi為各類服務(wù)功能價(jià)值系數(shù)[元/(hm2·a)]。

2 結(jié)果與分析

2.1 土地利用變化分析

2.1.1 不同土地利用類型變化分析 基于遙感解譯的2008—2018年昆崳山土地利用類型及其面積比例見圖1和表3。昆崳山土地利用類型為林地、耕地、建設(shè)用地、水域、裸地及其他。其中林地為最主要的土地利用類型、2008年、2013年、2018年面積分別為10 551 hm2，10 715 hm2，11 554 hm2，面積占比分別為67.75%，68.81%和74.19%，分布范圍較廣；其次為耕地，主要分布于昆崳山的西北方向，2008年、2013年、2018年面積分別為3 458 hm2，3 369 hm2，2 896 hm2，面積占比分別為22.21%，21.63%和18.59%；裸地及其他分布較為分散，2008年、2013年、2018年面積分別為1 238 hm2，1 059 hm2，666 hm2，面積占比分別為7.95%，6.8%和4.27%。

表3 昆崳山國(guó)家級(jí)自然保護(hù)區(qū)2008-2018年土地利用變化

圖1 昆崳山國(guó)家級(jí)自然保護(hù)區(qū)土地利用分類

2008—2018年期間，林地面積增加了1 003 hm2，所占面積比例由67.75%增至74.19%；建設(shè)用地和水域也略有增長(zhǎng)，分別增長(zhǎng)72 hm2，60 hm2，建設(shè)用地所占面積比例由0.87%增至1.33%，水域所占面積比例由1.22%增至1.61%；耕地面積減少562 hm2，所占面積比例由22.21%減少至18.59%；裸地及其他面積減少572 hm2，所占面積比例由7.95%減少至4.27%。

2.1.2 土地利用變化動(dòng)態(tài)度 2008—2018年，昆崳山不同土地利用類型土地利用變化動(dòng)態(tài)度見表4。從單一土地利用動(dòng)態(tài)度(K)來(lái)看，10 a間，林地、建設(shè)用地和水域呈正向變化趨勢(shì)，動(dòng)態(tài)變化分別為0.95%，5.33%和3.14%；耕地、裸地及其他呈負(fù)向變化趨勢(shì)，動(dòng)態(tài)變化度分別為-1.63%和-4.62%，二者是林地、建設(shè)用地和水域面積增長(zhǎng)的主要來(lái)源。2008—2013年林地、耕地、裸地及其他土地利用類型的單一土地利用動(dòng)態(tài)度變化程度低于2013—2018年，但水域和建設(shè)用地在2008—2013年的單一土地利用動(dòng)態(tài)度變化程度高于2013—2018年。

表4 昆崳山國(guó)家級(jí)自然保護(hù)區(qū)2008—2018年期間土地利用變化動(dòng)態(tài)度 %

整體來(lái)看，2008—2018年綜合土地利用動(dòng)態(tài)度(LC)為0.68%，2008—2013年綜合土地利用動(dòng)態(tài)度為1.35%，2013—2018年綜合土地利用動(dòng)態(tài)度為1.02%。

2.1.3 土地利用類型轉(zhuǎn)移特征 根據(jù)昆崳山土地利用變化轉(zhuǎn)移矩陣可知(表5)，2008—2018年，土地利用類型變化主要是裸地及其他、耕地和林地之間的轉(zhuǎn)移，轉(zhuǎn)入面積最多的土地利用類型是林地，10 a間轉(zhuǎn)入面積1 129.93 hm2；轉(zhuǎn)出面積較大的土地利用類型是耕地和裸地及其他，分別轉(zhuǎn)出699.98 hm2，744.37 hm2。2008—2013年及2013—2018年，土地利用類型變化依然是林地和耕地、裸地及其他之間的轉(zhuǎn)移，整體變化趨勢(shì)與十年間是一致的。

表5 不同時(shí)期昆崳山國(guó)家級(jí)自然保護(hù)區(qū)土地利用變化轉(zhuǎn)移矩陣 hm2

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

由表6可知，2008—2018年，昆崳山生態(tài)系統(tǒng)服務(wù)價(jià)值(ESV)增加了3 201.82萬(wàn)元，增長(zhǎng)率達(dá)9.53%。其中，2008—2013年，昆崳山ESV增長(zhǎng)較少，僅增長(zhǎng)1 121.53萬(wàn)元；2013—2018年，昆崳山ESV增長(zhǎng)較為明顯，增長(zhǎng)2 080.29萬(wàn)元。生態(tài)服務(wù)功能一級(jí)分類中，對(duì)昆崳山ESV貢獻(xiàn)率最大的是調(diào)節(jié)服務(wù)價(jià)值，2008年、2013年、2018年分別占到了ESV總值的63.92%，64.29%和64.46%；其次為支持服務(wù)價(jià)值，2008年、2013年、2018年分別占到了ESV總值的24.01%，23.72%和23.74%；文化服務(wù)的比例最低，2008年、2013年、2018年分別占到了ESV總值的4.68%，4.67%和4.71%。

生態(tài)服務(wù)功能二級(jí)分類中，氣候調(diào)節(jié)、水文調(diào)節(jié)、土壤保持價(jià)值的貢獻(xiàn)率較大，其中貢獻(xiàn)率最大的為氣候調(diào)節(jié)(表6)。2008年、2013年、2018年，氣候調(diào)節(jié)的ESV分別為9 278.62萬(wàn)元、9 438.55萬(wàn)元和10 136.58萬(wàn)元，分別占當(dāng)年ESV的27.63%，27.20%和27.56%。各項(xiàng)生態(tài)服務(wù)功能ESV比例各年度變化不大，表明研究區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值格局比較穩(wěn)定。

表6 2008-2018年昆崳山國(guó)家級(jí)自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值

昆崳山2008—2018年期間各土地利用類型生態(tài)系統(tǒng)服務(wù)價(jià)值見表7。不同土地利用類型中，林地和水域的生態(tài)系統(tǒng)服務(wù)價(jià)值呈增長(zhǎng)趨勢(shì)，耕地、裸地及其他呈減少趨勢(shì)。2008年、2013年、2018年，林地ESV貢獻(xiàn)率分別達(dá)到88.58%，87.05%和88.56%，表明林地對(duì)昆崳山生態(tài)系統(tǒng)服務(wù)功能發(fā)揮有著至關(guān)重要的作用。5種不同土地利用類型中，ESV由大到小依次為：林地>水域>耕地>裸地及其他>建設(shè)用地。

表7 2008-2018年期間昆崳山國(guó)家級(jí)自然保護(hù)區(qū)各土地利用類型生態(tài)系統(tǒng)服務(wù)價(jià)值

3 討 論

土地利用格局變化是表征一個(gè)趨于生態(tài)環(huán)境變化的重要指標(biāo)，對(duì)不同土地利用類型的生態(tài)系統(tǒng)服務(wù)價(jià)值有著顯著影響[29-31]。昆崳山自然保護(hù)區(qū)成立以后，為進(jìn)一步加強(qiáng)對(duì)森林資源的保護(hù)，2008—2018年，保護(hù)區(qū)采取了一系列荒山造林、退耕還林、生態(tài)修復(fù)等措施，實(shí)施了“綠水青山生態(tài)修復(fù)”、“大地環(huán)境全面整治”、“原生動(dòng)物種群保育”等工程，同時(shí)加大了基礎(chǔ)設(shè)施建設(shè)的力度，在土地利用變化上表現(xiàn)為昆崳山的林地、建設(shè)用地(主要為防火道路建設(shè))和水域面積呈增長(zhǎng)趨勢(shì)，耕地、裸地及其他面積呈減少趨勢(shì)。林地是昆崳山主要的土地利用類型，2018年林地面積占到了昆崳山總面積的74.19%。2008—2018年昆崳山自然保護(hù)區(qū)土地利用類型變化主要是裸地及其他、耕地和林地之間的轉(zhuǎn)移，面積增加最多的土地利用類型為林地，林地增加1 003 hm2，面積比例增加6.44%，10 a間林地的ESV增長(zhǎng)2 828.13萬(wàn)元，雖然占總ESV的比例略有降低，但仍然是維持昆崳山生態(tài)系統(tǒng)服務(wù)的主體，這與昆崳山作為森林生態(tài)類型的自然保護(hù)區(qū)的生態(tài)定位相一致；面積減少較大的土地利用類型是耕地、裸地及其他，其ESV也呈現(xiàn)出下降的趨勢(shì)；雖然建設(shè)用地在面積上有一定程度上的增長(zhǎng)，但由于整體所占比例較低，土地利用開發(fā)強(qiáng)度不高，其產(chǎn)生的效應(yīng)對(duì)生態(tài)系統(tǒng)造成的負(fù)面效應(yīng)微弱[32]。各土地利用類型變化趨勢(shì)及其ESV的變化，表明昆崳山自然保護(hù)采取的一系列生態(tài)修復(fù)措施實(shí)施效果顯著，昆崳山國(guó)家級(jí)自然保護(hù)區(qū)土地利用結(jié)構(gòu)朝著增強(qiáng)其生態(tài)功能的方向發(fā)展，也印證了適度人為干預(yù)對(duì)森林生態(tài)系統(tǒng)具有一定正向促進(jìn)作用[33]。

當(dāng)量因子法主要是基于各種服務(wù)的價(jià)值當(dāng)量，結(jié)合生態(tài)系統(tǒng)的分布面積對(duì)生態(tài)系統(tǒng)服務(wù)進(jìn)行評(píng)估，小尺度區(qū)域上運(yùn)用當(dāng)量因子法核算生態(tài)系統(tǒng)服務(wù)價(jià)值具有更高的可靠性和實(shí)踐性[34]。運(yùn)用當(dāng)量因子法對(duì)昆崳山2008—2018年的生態(tài)系統(tǒng)服務(wù)價(jià)值及其變化進(jìn)行測(cè)算，對(duì)昆崳山自然保護(hù)區(qū)生態(tài)系統(tǒng)的保護(hù)與可持續(xù)發(fā)展具有指導(dǎo)意義。研究發(fā)現(xiàn)昆崳山ESV為3.5億元/年，2008—2018年間，昆崳山ESV增加3 201.82萬(wàn)元，增長(zhǎng)率達(dá)9.53%。生態(tài)系統(tǒng)服務(wù)功能11個(gè)二級(jí)指標(biāo)中，食物生產(chǎn)功能價(jià)值呈減少趨勢(shì)，其余功能價(jià)值均呈不同程度的增長(zhǎng)趨勢(shì)，其中氣候調(diào)節(jié)和水文調(diào)節(jié)功能價(jià)值增長(zhǎng)明顯。生態(tài)功能是昆崳山的主要服務(wù)功能，其中調(diào)節(jié)服務(wù)是昆崳山的主要生態(tài)功能，其ESV貢獻(xiàn)率在2008年、2013年、2018年分別達(dá)到63.92%，64.29%和64.46%。對(duì)同一個(gè)研究區(qū)域，指標(biāo)體系、評(píng)估方法的不同會(huì)導(dǎo)致評(píng)估結(jié)果的差異[34-35]，王玉濤采用影子工程法、機(jī)會(huì)成本法、市場(chǎng)價(jià)值法等研究方法，得出昆崳山自然保護(hù)區(qū)的生態(tài)系統(tǒng)服務(wù)功能價(jià)值達(dá)4.99億元/年，主要的生態(tài)系統(tǒng)服務(wù)功能為涵養(yǎng)水源和凈化空氣[22]，這與本研究結(jié)果存在一定差異。但按照直接價(jià)值(供給服務(wù))和間接價(jià)值(調(diào)解服務(wù)、支持服務(wù)和文化服務(wù))分類，研究結(jié)果都表明昆崳山生態(tài)系統(tǒng)服務(wù)功能主體是間接服務(wù)，與歐陽(yáng)志云等人的研究結(jié)果一致[36]，這也進(jìn)一步表明了昆崳山主體功能是生態(tài)功能，而不是原材料供給?；谏鷳B(tài)系統(tǒng)服務(wù)價(jià)值的昆崳山生態(tài)服功能分類結(jié)果對(duì)保護(hù)區(qū)生態(tài)保護(hù)及相關(guān)措施的實(shí)施具有一定指導(dǎo)意義，保護(hù)區(qū)應(yīng)合理規(guī)劃昆崳山土地利用格局，嚴(yán)格保護(hù)森林生態(tài)系統(tǒng)，維持其提供生態(tài)功能的能力。

4 結(jié) 論

本研究探究了2008—2018年昆崳山自然保護(hù)區(qū)不同土地利用類型動(dòng)態(tài)變化，運(yùn)用當(dāng)量因子法估算了昆崳山國(guó)家級(jí)自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值。10 a間，不同土地利用類型中，面積增加最多的土地利用類型是林地，林地面積增加1 003 hm2，面積比例增加6.44%；面積減少較大的土地利用類型是耕地、裸地及其他，面積分別減少562 hm2和 572 hm2，面積比例分別下降3.62%和3.68%。10 a間林地、建設(shè)用地和水域呈增長(zhǎng)趨勢(shì)，其他類型土地呈下降趨勢(shì)。從單一土地利用動(dòng)態(tài)度來(lái)看，10 a間，林地、建設(shè)用地和水域動(dòng)態(tài)變化分別為0.95%，5.33%和3.14%；耕地、裸地及其他呈動(dòng)態(tài)變化度分別為-1.63%和-4.62%，土地利用類型變化主要是裸地及其他、耕地與林地之間的轉(zhuǎn)移。

不同土地利用類型中，林地是昆崳山主要的土地利用類型，2018年林地面積占昆崳山總面積的74.19%，林地的生態(tài)系統(tǒng)服務(wù)價(jià)值呈增長(zhǎng)趨勢(shì)，2008—2018年，林地的ESV增長(zhǎng)2 828.13萬(wàn)元，2018年林地ESV貢獻(xiàn)率達(dá)到88.56%。林地對(duì)昆崳山生態(tài)系統(tǒng)服務(wù)功能發(fā)揮有著至關(guān)重要的作用，是維持昆崳山生態(tài)系統(tǒng)服務(wù)的主體。

昆崳山自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值10 a間呈增長(zhǎng)趨勢(shì)，ESV增加了3 201.82萬(wàn)元，增長(zhǎng)率達(dá)9.53%。各類生態(tài)系統(tǒng)服務(wù)功能中，除食物生產(chǎn)功能價(jià)值有所減少，其余功能價(jià)值均呈不同程度的增長(zhǎng)。調(diào)節(jié)服務(wù)價(jià)值10 a間增長(zhǎng)2 242.32萬(wàn)元，占ESV總值的比例提高0.54%。二級(jí)類型中，氣候調(diào)節(jié)和水文調(diào)節(jié)功能價(jià)值增長(zhǎng)最為明顯，氣候調(diào)節(jié)的ESV 10年間增長(zhǎng)857.96萬(wàn)元，水文調(diào)節(jié)的ESV 10年間增長(zhǎng)872.64萬(wàn)元，調(diào)節(jié)服務(wù)是昆崳山自然保護(hù)區(qū)的主要生態(tài)功能。

昆崳山的土地利用和ESV的變化都向著有利于促進(jìn)森林生態(tài)系統(tǒng)生態(tài)功能發(fā)揮的方向發(fā)展，但作為膠東半島重要的生態(tài)屏障，昆崳山生態(tài)功能的發(fā)揮程度對(duì)于維持區(qū)域生態(tài)安全具有重要意義，雖然前期保護(hù)區(qū)通過荒山造林、退耕還林、加大基礎(chǔ)設(shè)施建設(shè)力度、生態(tài)修復(fù)等措施的實(shí)施，有效維護(hù)了昆崳山生態(tài)系統(tǒng)的安全，不斷提升了生態(tài)系統(tǒng)服務(wù)功能，但生態(tài)系統(tǒng)保護(hù)和修復(fù)是一項(xiàng)長(zhǎng)期工程，今后應(yīng)在加大保護(hù)力度基礎(chǔ)上，結(jié)合自然保護(hù)地整合優(yōu)化，進(jìn)一步優(yōu)化自然保護(hù)地布局，合理利用并規(guī)劃土地利用，促進(jìn)昆崳山自然保護(hù)區(qū)生態(tài)系統(tǒng)健康穩(wěn)定發(fā)展。

本研究?jī)H測(cè)算了昆崳山國(guó)家級(jí)自然保護(hù)區(qū)建立后10 a間保護(hù)區(qū)土地利用格局及生態(tài)系統(tǒng)服務(wù)功能的變化，沒有涉及更長(zhǎng)時(shí)間的變化，在下步研究中，將進(jìn)一步擴(kuò)大研究的時(shí)間序列，對(duì)昆崳山國(guó)家級(jí)自然保護(hù)區(qū)的土地利用格局及生態(tài)系統(tǒng)服務(wù)功能的變化開展更為完善的研究，為區(qū)域生態(tài)政策的制定和生態(tài)保護(hù)與建設(shè)工程的實(shí)施提供更多的科學(xué)依據(jù)。