叢喜東 李琳 齊永峰 張東來

摘 要：為解決牡丹江森林生態(tài)系統(tǒng)服務(wù)功能對(duì)生態(tài)管理需求的重要性缺乏全面評(píng)估的問題，以黑龍江省牡丹江市森林生態(tài)系統(tǒng)為研究對(duì)象，參考《生態(tài)保護(hù)紅線劃定技術(shù)指南》及當(dāng)?shù)刈匀坏乩硖卣鳎ㄟ^水源涵養(yǎng)、水土保持、防風(fēng)固沙及生物多樣性維護(hù)4個(gè)生態(tài)服務(wù)功能評(píng)估指標(biāo)，構(gòu)建生態(tài)系統(tǒng)服務(wù)功能重要性評(píng)價(jià)模型。結(jié)果表明，牡丹江市森林覆蓋率達(dá)67.40%，2000—2020年森林植被凈初級(jí)生產(chǎn)力（NPP）平均為（730.06±56.64） g/（m·a），且呈現(xiàn)出極顯著的增加趨勢(shì)；單項(xiàng)評(píng)價(jià)表明，當(dāng)?shù)厣稚鷳B(tài)系統(tǒng)的重要性級(jí)別最高的為水土保持功能，綜合評(píng)價(jià)表明，生態(tài)系統(tǒng)服務(wù)功能“極重要”等級(jí)面積最大，占森林總面積的48.74%，集中分布的區(qū)域有寧安市西南和穆棱市南部。結(jié)果可為當(dāng)?shù)貎?yōu)化生態(tài)保護(hù)策略，協(xié)調(diào)區(qū)域發(fā)展供科學(xué)依據(jù)。

關(guān)鍵詞：生態(tài)系統(tǒng)服務(wù)功能；森林；牡丹江；評(píng)價(jià)模型；遙感

中圖分類號(hào)：718.55 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1006-8023（2023）06-0075-06

Evaluation on the Importance of Forest Ecosystem Services in Mudanjiang

CONG Xidong， LI Lin， QI Yongfeng， ZHANG Donglai

（1. Heilongjiang Ecology Institute， Harbin 150081， China; 2.Key Laboratory of Forest Ecology and Forestry Ecological Engineering， Heilongjiang Ecology Institute， Harbin 150081， China; 3.Mudanjiang Forest Ecosystem Research Station， Mudanjiang 157500， China）

Abstract：In order to solve lack of comprehensive assessment of forest ecosystem service functions on ecological management needs in Mudanjiang City， taking the forest ecosystem of Mudanjiang City in Heilongjiang Province as the research object， referring to the 'Technical Guidelines for Delineation of Ecological Protection Red Line' and the local physical and geographical characteristics， a model for evaluating the importance of forest ecosystem service functions was constructed through four ecological service function evaluation indicators： water conservation， soil and water conservation， wind prevention and sand fixation and biodiversity conservation. The results showed that the forest coverage rate in Mudanjiang City reached 67.40% and the average net primary productivity of forest vegetation was （730.06 ± 56.64） g/（m·a） from 2000 to 2020， which showing a significant increasing trend. From the single evaluation index， the importance of soil and water conservation function for the local forest ecosystem was the highest. From the comprehensive evaluation results， the area of extremely important grade of the ecosystem service function was the largest， accounting for 48. 74% of the total forest area， which was concentrated in the southwest of Ning'an City and the south of Muling City. The results could provide scientific basis for optimizing local ecological protection strategies and coordinating regional development.

Keywords：Ecosystem service function; forest; Mudanjiang City; evaluation model; remote sensing

收稿日期：2023-03-22

基金項(xiàng)目：中央財(cái)政林業(yè)科技推廣示范項(xiàng)目（黑[2022]TG24）；2020年度省屬科研院所科研業(yè)務(wù)費(fèi)支持項(xiàng)目（A類CZKYF2020A003）；黑龍江省重點(diǎn)研發(fā)計(jì)劃（GA21C030）；黑龍江省重點(diǎn)研發(fā)計(jì)劃指導(dǎo)類項(xiàng)目 （GZ20220117）。

第一作者簡(jiǎn)介：叢喜東，碩士，副研究員。研究方向?yàn)榱謽I(yè)遙感與生態(tài)信息監(jiān)測(cè)、林業(yè)信息工程、林業(yè)智能裝備和林業(yè)物聯(lián)網(wǎng)技術(shù)等。E-mail： xidongs@live.com

引文格式：叢喜東，李琳，齊永峰，等. 牡丹江森林生態(tài)系統(tǒng)服務(wù)功能重要性評(píng)價(jià)[J].森林工程，2023，39（6）：75-80.

CONG X D， LI L， QI Y F， et al. Evaluation on the importance of forest ecosystem services in Mudanjiang[J]. Forest Engineering， 2023， 39（6）：75-80.

0 引言

生態(tài)系統(tǒng)服務(wù)代表生態(tài)系統(tǒng)為人類提供福祉的能力，綜合評(píng)價(jià)生態(tài)系統(tǒng)服務(wù)功能的重要性等級(jí)，識(shí)別重要區(qū)域的空間分布和異質(zhì)性特征，對(duì)于實(shí)施生態(tài)保護(hù)、促進(jìn)區(qū)域可持續(xù)發(fā)展具有重要意義，森林資源開發(fā)促進(jìn)人類福祉與生態(tài)系統(tǒng)保護(hù)之間的平衡是實(shí)現(xiàn)可持續(xù)發(fā)展的關(guān)鍵。自生態(tài)系統(tǒng)服務(wù)概念提出以來，國(guó)內(nèi)學(xué)者就此開展了廣泛研究，劉曦等、肖驍?shù)仍u(píng)估了東北森林生態(tài)系統(tǒng)服務(wù)功能的價(jià)值量，舒航等利用價(jià)值轉(zhuǎn)移模型評(píng)估了森林類保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)功能，黃龍生等研究了濟(jì)南市森林生態(tài)系統(tǒng)服務(wù)功能的空間格局，殷格蘭等以森林為重點(diǎn)，研究了土地覆蓋變化對(duì)生態(tài)系統(tǒng)服務(wù)功能的影響。

作為陸地生態(tài)系統(tǒng)的主體，森林在水源涵養(yǎng)、水土保持和防風(fēng)固沙等生態(tài)服務(wù)方面具有其他生態(tài)系統(tǒng)不可比擬的優(yōu)勢(shì)，然而氣候變化、森林退化、耕地和城市侵占等諸多因素正在威脅著森林資源及其服務(wù)功能的可持續(xù)發(fā)展。牡丹江市位于長(zhǎng)白山脈北坡，森林覆蓋率超過60%，動(dòng)植物資源豐富，森林生態(tài)系統(tǒng)功能對(duì)當(dāng)?shù)厣鷳B(tài)環(huán)境和經(jīng)濟(jì)社會(huì)發(fā)展意義重大，但目前當(dāng)?shù)氐纳稚鷳B(tài)系統(tǒng)服務(wù)功能對(duì)生態(tài)管理需求的重要性尚未得到全面評(píng)估。本研究正是針對(duì)以上問題，參照《生態(tài)保護(hù)紅線劃定技術(shù)指南》提供的技術(shù)方法，基于遙感、氣象、土壤、地形和地理信息等多源數(shù)據(jù)，通過水源涵養(yǎng)、水土保持、防風(fēng)固沙及生物多樣性維護(hù)4個(gè)生態(tài)服務(wù)功能評(píng)估指標(biāo)，構(gòu)建生態(tài)系統(tǒng)服務(wù)功能重要性評(píng)價(jià)模型，對(duì)森林生態(tài)系統(tǒng)服務(wù)功能重要性進(jìn)行綜合評(píng)價(jià)，研究結(jié)果對(duì)當(dāng)?shù)厣鷳B(tài)區(qū)劃保護(hù)和分類管理具有應(yīng)用價(jià)值。

1 研究區(qū)與研究方法

1.1 研究區(qū)概況

牡丹江市位于黑龍江省東南部，如圖1所示，地處中、俄、朝合圍的“金三角”腹地，山地和丘陵是主要地形特征，海拔平均為230 m。研究區(qū)氣候特征為海洋（半濕潤(rùn)型）中溫帶季風(fēng)氣候，雨熱同季，年平均氣溫4.5 ℃，年降水量563.1 mm。研究區(qū)森林覆蓋率超過60%，動(dòng)植物資源豐富，樹種有包括紅松（Pinus koraiensis）、落葉松（Larix gmelinii）、樟子松（Pinus syivestris var. mongolica）、云杉（Picea asperata）、冷杉（Abies fabri）和水曲柳（Fraxinus mandshurica Rupr.）等百余種，野生經(jīng)濟(jì)植物2 200余種，還有東北虎（Panthera tigris altaica）、梅花鹿（Cervus nippon）和狍子（Capreolus pygargus）等珍稀野生動(dòng)物。

1.2 研究數(shù)據(jù)及方法

1.2.1 氣象數(shù)據(jù)

利用中國(guó)科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心（https：//www.resdc.cn/Default.aspx）獲取氣象數(shù)據(jù)，用于水源涵養(yǎng)功能、防風(fēng)固沙功能和生物多樣性維護(hù)功能的重要性評(píng)價(jià)。研究區(qū)內(nèi)包括7個(gè)觀測(cè)站，觀測(cè)站分布如圖1所示。氣象要素為氣溫、降水量、相對(duì)濕度和平均風(fēng)速，本研究采用1991—2020年逐月平均值參與生態(tài)系統(tǒng)服務(wù)功能重要性評(píng)價(jià)。利用ARCGIS 10.2軟件的IDW工具對(duì)觀測(cè)站的各氣象要素進(jìn)行空間分辨率為500 m的插值處理，得到各氣象要素的空間分布圖。

1.2.2 數(shù)字高程模型

通過中國(guó)科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國(guó)際科學(xué)數(shù)據(jù)鏡像網(wǎng)站（http：//www.gscloud.cn）獲取空間分辨率為90 m數(shù)字高程模型（DEM），數(shù)據(jù)版本為SRTM地形產(chǎn)品V4.1。

1.2.3 土壤性質(zhì)數(shù)據(jù)

通過中國(guó)科學(xué)院地理科學(xué)與資源研究所資源環(huán)境數(shù)據(jù)云平臺(tái)（http：//www.resdc.cn/），獲取土壤質(zhì)地、土壤有機(jī)質(zhì)含量等土壤性質(zhì)數(shù)據(jù)。

1.2.4 遙感數(shù)據(jù)

本研究收集4景2020年Landsat8/OLI數(shù)據(jù)用于研究區(qū)土地覆蓋分類，OLI傳感器的空間分辨率為15～30 m。本研究將第2、4、5波段（30 m分辨率）數(shù)據(jù)和歸一化差異植被指數(shù)（NDVI）相結(jié)合，參考國(guó)家標(biāo)準(zhǔn)《土地利用現(xiàn)狀分類》（GB/T 21010—2017）中的分類體系，采用隨機(jī)森林模型，提取研究區(qū)森林的空間分布。將目視解譯提取的28 546個(gè)驗(yàn)證點(diǎn)用于精度驗(yàn)證，利用誤差矩陣法（Error Matrix）和Kappa的分析方法計(jì)算的總體分類精度為92.58%，Kappa系數(shù)為0.914 7。

同時(shí)，本研究應(yīng)用MODIS數(shù)據(jù)產(chǎn)品——MOD17A3的植被凈初級(jí)生產(chǎn)力（NPP）數(shù)據(jù)參與生態(tài)系統(tǒng)服務(wù)功能重要性評(píng)價(jià)，通過美國(guó)NASA LPDAAC（The Land Processes Distributed Active Archive Center）EOS數(shù)據(jù)中心獲取植被指數(shù)數(shù)據(jù)及其對(duì)應(yīng)的質(zhì)量信息數(shù)據(jù)集，選取植被凈初級(jí)生產(chǎn)力（NPP）和植被初級(jí)生產(chǎn)力（GPP）代表植被指數(shù)。

以上數(shù)據(jù)在ENVI軟件的支持下進(jìn)行數(shù)據(jù)重采樣，設(shè)置空間分辨率為500 m，投影轉(zhuǎn)換為Albers投影。

1.2.5 生態(tài)系統(tǒng)服務(wù)功能重要性評(píng)價(jià)指標(biāo)

參照《生態(tài)保護(hù)紅線劃定技術(shù)指南》提供的技術(shù)方法，基于遙感、氣象、土壤、地形和地理信息等多源數(shù)據(jù)，在Arcgis、ENVI和Origion等軟件的支持下，利用NPP定量指標(biāo)評(píng)價(jià)法，從水源涵養(yǎng)、水土保持、防風(fēng)固沙和生物多樣性保護(hù)4個(gè)維度，構(gòu)建生態(tài)系統(tǒng)服務(wù)功能重要性評(píng)價(jià)模型，對(duì)生態(tài)系統(tǒng)服務(wù)功能重要性進(jìn)行評(píng)價(jià)。

水源涵養(yǎng)功能重要性評(píng)價(jià)，計(jì)算公式如下

W=N×F×F×（1-F）。??? （1）

式中：W為生態(tài)系統(tǒng)水源涵養(yǎng)服務(wù)能力指數(shù)；N為區(qū)域多年生態(tài)系統(tǒng)凈初級(jí)生產(chǎn)力平均值；F為評(píng)價(jià)區(qū)域坡度；F為常年降水量；F為土壤滲流能力因子。

水土保持功能重要性評(píng)價(jià)，計(jì)算公式如下

S=N×（1-K）×（1-F）。??? （2）

式中：S為水土保持服務(wù)能力指數(shù)；K為土壤可蝕性因子；F為評(píng)價(jià)區(qū)域坡度。

防風(fēng)固沙功能重要性評(píng)價(jià)，計(jì)算公式如下

S=（1-N）×（1-K）×F×（1-D）。??? （3）

式中：S為防風(fēng)固沙服務(wù)能力指數(shù)；F為多年平均氣候侵蝕力；D為地表粗糙度因子。

生物多樣性保護(hù)功能重要性評(píng)價(jià)，計(jì)算公式如下

S=N×F×F×（1-F）。??? （4）

式中：S為生物多樣性保護(hù)服務(wù)能力指數(shù)；F為常年降水量；F為常年平均氣溫；F為海拔。

1.2.6 生態(tài)系統(tǒng)服務(wù)功能重要性評(píng)價(jià)模型構(gòu)建

以生態(tài)系統(tǒng)服務(wù)功能（水源涵養(yǎng)、水土保持、防風(fēng)固沙和生物多樣性保護(hù)）重要性評(píng)價(jià)指標(biāo)為基礎(chǔ)，構(gòu)建重要性評(píng)價(jià)模型

E=w·W+w·S+w·S+w·S。??? （5）

式中：E為生態(tài)系統(tǒng)功能重要性指數(shù)；w為權(quán)重，利用層次分析法確定，w=0.344 7，w=0.249 5，w=0.157 3，w=0.248 5。由于各評(píng)價(jià)指標(biāo)的量綱不一致，因此本研究將各指標(biāo)歸一化處理到0～1區(qū)間再計(jì)算E。

1.2.7 生態(tài)系統(tǒng)服務(wù)重要性等級(jí)劃分

采用自然斷點(diǎn)法對(duì)各生態(tài)系統(tǒng)服務(wù)功能重要性指標(biāo)進(jìn)行分類操作，按生態(tài)系統(tǒng)服務(wù)值劃分為一般重要、中等重要、重要和極重要4個(gè)等級(jí)。

2 結(jié)果與分析

2.1 牡丹江森林分布空間特征

全市森林總面積為26 127.22 km，森林覆蓋率達(dá)67.40%，見表1，其中海林市和東寧市森林面積較大，綏芬河和牡丹江市較?。粡纳指采w率角度分析，東寧市森林覆蓋率最高，其次是海林市，林口縣和牡丹江市最低。結(jié)合圖1和圖2分析可見，牡丹江森林主要分布在海拔相對(duì)較高的區(qū)域，如海林市、東寧市和綏芬河市的大部，寧安市的南部和東部，林口縣的西北部，以及穆棱市的南部區(qū)域。

2.2 森林生態(tài)系統(tǒng)NPP變化特征

2000—2020年牡丹江森林NPP平均為（730.06±56.64） g/（m·a），其中2020年最高（903.43 g/（m·a）），2011年最低（572.42 g/（m·a））。年際間呈極顯著增加趨勢(shì)，如圖3所示，平均每年增加7.78 g/（m·a），尤其在經(jīng)歷了2011—2012年低谷期以后，從2013年開始以13.48 g/（m·a）的速度顯著增長(zhǎng)（線性擬合為y=13.48x-26 373.95，R=0.550 1，P=0.035）。

2000年以來研究區(qū)森林NPP平均值為297.43～808.09 g/（m·a），NPP值分布頻數(shù)呈明顯的單峰曲線，如圖4（b）所示，峰值在747 g/（m·a）。結(jié)合圖4（a）和圖4（b）分析可見，NPP在700～750 g/（m·a）

的面積最大，占研究區(qū)森林總面積的48.33%，集中分布的區(qū)域在海林市西部和穆棱市南部；其次是NPP>750 g/（m·a）的區(qū)域，占34.15%，集中分布的區(qū)域在海林市西南部、寧安市大部和東寧市大部。

2.3 森林生態(tài)系統(tǒng)功能重要性

2.3.1 單項(xiàng)服務(wù)功能指標(biāo)的重要性評(píng)價(jià)

根據(jù)生態(tài)系統(tǒng)服務(wù)功能重要性評(píng)價(jià)方法，確定各個(gè)服務(wù)功能的重要性等級(jí)，如圖5所示，總體上看處于最重要等級(jí)的為水土保持功能，水源涵養(yǎng)和生物多樣性保護(hù)2個(gè)功能的重要性級(jí)別略低，防風(fēng)固沙重要性級(jí)別最低。分別各單項(xiàng)指標(biāo)分析可見，水源涵養(yǎng)功能重要性以“重要”和“中等重要”級(jí)別為主，分別占研究區(qū)森林總面積的47.50%和35.20%，穿插分布在研究區(qū)大部分區(qū)域；“極重要”級(jí)別最少，占7.78%，集中分布在西南部山區(qū)；水土保持功能重要性也以“重要”和“中等重要”級(jí)別為主，分別占研究區(qū)森林總面積的40.44%和43.75%，但分布格局與水源涵養(yǎng)重要性級(jí)別不同，“重要”級(jí)別分布更為集中；防風(fēng)固沙功能以“一般重要”級(jí)別為主，占研究區(qū)森林總面積的90.68%；在生物多樣性保護(hù)重要性方面，“中等重要”和“重要”級(jí)別占比分別為31.64%和39.36%，“極重要”級(jí)別占比略小，為21.52%。

2.3.2 生態(tài)系統(tǒng)功能重要性綜合評(píng)價(jià)

對(duì)牡丹江森林生態(tài)系統(tǒng)功能重要性進(jìn)行綜合評(píng)價(jià)，結(jié)果如圖6所示，“極重要”等級(jí)面積最大，占森林總面積的48.74%，集中分布的區(qū)域有寧安市西南和穆棱市南部；“重要”等級(jí)的區(qū)域占森林總面積的36.40%，集中分布的地塊在穆棱市東北部與綏芬河北部交界處；“中等重要”等級(jí)區(qū)域占13.07%，主要分布在海林市西部和穆棱市部分區(qū)域；“一般重要”等級(jí)區(qū)域占比最小，為1.78%，呈零星分布狀態(tài)。

3 結(jié)論

本研究通過構(gòu)建生態(tài)系統(tǒng)服務(wù)功能重要性評(píng)價(jià)模型，對(duì)牡丹江森林生態(tài)系統(tǒng)服務(wù)功能重要性進(jìn)行綜合評(píng)價(jià)，主要結(jié)論如下：全市森林總面積為26 127.22 km，其中海林市和東寧市森林面積較大，森林覆蓋率較高；2000年以來森林NPP呈極顯著的線性增加趨勢(shì)，平均每年增加7.78 g/（m·a）；森林生態(tài)系統(tǒng)服務(wù)功能“極重要”等級(jí)中，以水土保持功能為主，占森林總面積的48.74%，集中分布的區(qū)域有寧安市西南和穆棱市南部。

另外，本研究雖然分析了2000年以來森林NPP的變化趨勢(shì)，但未對(duì)NPP增加的驅(qū)動(dòng)因素進(jìn)行深入分析，因此明確氣候、大氣CO濃度變化、氮沉降和森林?jǐn)_動(dòng)等因子對(duì)NPP動(dòng)態(tài)變化的貢獻(xiàn)是今后研究的重要內(nèi)容之一。

【參 考 文 獻(xiàn)】

[1]HAYHA T， FRANZESE P P. Ecosystem services assessment： a review under an ecological-economic and systems perspective[J]. Ecological Modelling， 2014， 289： 124-132.

[2]HAYHA T， FRANZESE P P. Ecosystem services assessment： a review under an ecological-economic and systems perspective[J]. Ecological Modelling， 2014， 289： 124-132.

[3]LIU Y， SHEN D. Assessment on the value of ecosystem service of important protected area for water supply of city[C]//2011 Second International Conference on Mechanic Automation and Control Engineering. July 15-17， 2011， Inner Mongolia， China. IEEE， 2011： 2962-2966.

[4]WANG J F， LI Y， WANG S， et al. Assessment of multiple ecosystem services and ecological security pattern in Shanxi Province， China[J]. IJERPH， 2023.

[5]RUKUNDO E， LIU S L， DONG Y H， et al. Spatio-temporal dynamics of critical ecosystem services in response to agricultural expansion in Rwanda， East Africa[J]. Ecological Indicators， 2018， 89： 696-705.

[6]劉曦，劉經(jīng)偉.東北國(guó)有林區(qū)森林生態(tài)系統(tǒng)服務(wù)功能價(jià)值量的監(jiān)測(cè)與評(píng)估[J].東北林業(yè)大學(xué)學(xué)報(bào)，2020，48（8）：66-71.

LIU X， LIU J W. Monitoring and valuing ecosystem service in northeast state-owned forests in China[J]. Journal of Northeast Forestry University， 2020， 48（8）： 66-71.

[7]肖驍，穆治霖，趙雪雁，等.基于RS/GIS的東北地區(qū)森林生態(tài)系統(tǒng)服務(wù)功能價(jià)值評(píng)估[J].生態(tài)學(xué)雜志，2017，36（11）：3298-3304.

XIAO X， MU Z L， ZHAO X Y， et al. RS/GIS-based evaluation of forest ecosystem service value in Northeast China[J]. Chinese Journal of Ecology， 2017， 36（11）： 3298-3304.

[8]舒航，莊立文，孫曉杰，等.價(jià)值轉(zhuǎn)移模型在森林類保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)功能評(píng)估中的應(yīng)用[J].東北林業(yè)大學(xué)學(xué)報(bào)，2020，48（12）：52-57.

SHU H， ZHUANG L W， SUN X J， et al. Application of value transfer model in the evaluation of ecosystem service function of forest reserve[J]. Journal of Northeast Forestry University， 2020， 48（12）： 52-57.

[9]黃龍生，王兵，牛香，等.濟(jì)南市森林生態(tài)系統(tǒng)服務(wù)功能空間格局研究[J].生態(tài)學(xué)報(bào)，2019，39（17）：6477-6486.

HUANG L S， WANG B， NIU X， et al. Spatial pattern of the ecosystem service function of forests in Jinan City[J]. Acta Ecologica Sinica， 2019， 39（17）： 6477-6486.

[10]殷格蘭，邵景安，郭躍，等.林地資源變化對(duì)森林生態(tài)系統(tǒng)服務(wù)功能的影響：以南水北調(diào)核心水源地淅川縣為例[J].生態(tài)學(xué)報(bào)，2017，37（20）：6973-6985.

YIN G L， SHAO J A， GUO Y， et al. Effects of forest resource changes on ecosystem function： a case study of Xichuan County， the main source of water for south-to-north water transfer in China[J]. Acta Ecologica Sinica， 2017， 37（20）： 6973-6985.

[11]VERMA M， SINGH R， NEGANDHI D. Forest ecosystem： functions， value and management[M]//SANDHU H. Ecosystem Functions and Management. Cham： Springer， 2017： 101-121.

[12]白志強(qiáng)，王文棟，李緩，等.喀納斯自然保護(hù)區(qū)森林生態(tài)系統(tǒng)服務(wù)功能價(jià)值評(píng)估[J].西北林學(xué)院學(xué)報(bào)，2018，33（6）：292-296.

BAI Z Q， WANG W D， LI H， et al. Evaluation of forest ecosystem service function of Kanas National Nature Reserve[J]. Journal of Northwest Forestry University， 2018， 33（6）： 292-296.

[13]崔亞琴，樊蘭英，劉隨存，等.山西省森林生態(tài)系統(tǒng)服務(wù)功能評(píng)估[J].生態(tài)學(xué)報(bào)，2019，39（13）：4732-4740.

CUI Y Q， FAN L Y， LIU S C， et al. Evaluation of forest ecosystem services value in Shanxi Province[J]. Acta Ecologica Sinica， 2019， 39（13）： 4732-4740.

[14]ZHANG Y Z， LIANG S L. Changes in forest biomass and linkage to climate and forest disturbances over Northeastern China[J]. Global Change Biology， 2014， 20（8）： 2596-2606.

[15]WANG X F， YANG B. Divergent tree radial growth at alpine coniferous forest ecotone and corresponding responses to climate change in northwestern China[J]. Ecological Indicators， 2021， 121： 107052.

[16]PATRA S， KUMAR A， SAIKIA P. Deforestation and forests degradation impacts on livelihood security and climate change： Indian initiatives towards its mitigation[M]//Environmental Degradation： Challenges and Strategies for Mitigation. Cham： Springer， 2022： 371-392.

[17]BETTS M G， YANG Z Q， HADLEY A S， et al. Forest degradation drives widespread avian habitat and population declines[J]. Nature Ecology & Evolution， 2022， 6（6）： 709-719.

[18]ZAKKAK S， KAKALIS E， RADOVIC A， et al. The impact of forest encroachment after agricultural land abandonment on passerine bird communities： the case of Greece[J]. Journal for Nature Conservation， 2014， 22（2）： 157-165.

[19]LIU S D， ZHANG J J， GENG Y H， et al. Plausible response of urban encroachment on ecological land to tourism growth and implications for sustainable management， a case study of Zhangjiajie， China[J]. Ecological Indicators， 2021， 132： 108253.

[20]WAN J Z， WANG C J， QU H， et al. Vulnerability of forest vegetation to anthropogenic climate change in China[J]. Science of the Total Environment， 2018， 621： 1633-1641.

[21]WANG X Y， ZHAO C Y， JIA Q Y. Impacts of climate change on forest ecosystems in northeast China[J]. Advances in Climate Change Research， 2013， 4（4）： 230-241.