恩施州生態(tài)產(chǎn)品價值核算及驅(qū)動因素分析

譚宇輝 韋正崢 張紹山 羅玉梅 王永飛 張勁

摘要 恩施州是長江中上游重要生態(tài)屏障，具有重大的生態(tài)保護價值和生態(tài)產(chǎn)品價值（生態(tài)系統(tǒng)生產(chǎn)總值，GEP）開發(fā)潛質(zhì)。基于功能價值法核算恩施州2011—2020年的GEP，利用時空地理加權(quán)回歸模型分析其變化特點及主要驅(qū)動因素，結(jié)果表明：恩施州GEP總量呈增長趨勢，年平均增長率2.50%，物質(zhì)供給價值增速＞文化服務(wù)價值增速＞調(diào)節(jié)服務(wù)價值增速；恩施州 GEP與GDP雙增長，經(jīng)濟社會與生態(tài)環(huán)境保護協(xié)調(diào)發(fā)展，各縣市間生態(tài)產(chǎn)品供給水平具有明顯的時空異質(zhì)性；年降水量、林地面積、人均GDP、城鎮(zhèn)化率對恩施州GEP有顯著影響；年降水量、林地面積、人均GDP與GEP正相關(guān)；城鎮(zhèn)化率對GEP的影響由負相關(guān)變?yōu)檎嚓P(guān)，具有臨界閾值效應(yīng)。

關(guān)鍵詞 生態(tài)產(chǎn)品價值；功能價值法；時空異質(zhì)性；恩施州

中圖分類號 F124.5? 文獻標(biāo)識碼 A

文章編號 0517-6611（2024）06-0068-07

doi：10.3969/j.issn.0517-6611.2024.06.016

開放科學(xué)（資源服務(wù)）標(biāo)識碼（OSID）：

Analysis on the Value and Driving Force of Ecosystem Product in Enshi Prefecture

TAN Yu-hui1，WEI Zheng-zheng2，ZHANG Shao-shan3 et al

（1.Enshi Ecology and Environment Bureau， Enshi，Hubei 445000；2.Policy Research Center for Environment and Economy， Ministry of Ecology and Environment of the Peoples Republic of China， Beijing 100029;3.Enshi Environmental Monitoring Station， Enshi， Hubei 445000）

Abstract Enshi Prefecture is an important ecological barrier in the middle and upper reaches of the Yangtze River，having significant ecosystem protection value and potential for the utilization and realization of GEP .This paper uses the functional value method to calculate the GEP in Enshi Prefecture from 2011 to 2020，analyzes its change trend，and explores its driving factors. The results show that，the total amount of GEP in Enshi Prefecture has shown an increasing trend，with an average annual growth rate of 2.5%，the growth rate of material supply value is greater than that of cultural service value，and the growth rate of regulated service value is greater;Enshi Prefectures GEP and GDP have achieved dual growth，with coordinated development of economy，society，and ecological environment protection，the supply level of ecological products among counties and cities exhibits significant spatiotemporal heterogeneity;The annual precipitation，forest area，per capita GDP，and urbanization rate have a significant impact on the GEP of Enshi Prefecture;The annual precipitation，forest area，per capita GDP are positively correlated with GEP;The impact of urbanization rate on GEP changes from negative correlation to positive correlation，with a critical threshold effect.

Key words Gross ecosystem product;Functional value method;Spatiotemporal heterogeneity;Enshi Prefecture

生態(tài)產(chǎn)品價值，又稱生態(tài)系統(tǒng)生產(chǎn)總值（gross ecosystem product，GEP） ，指一定區(qū)域的生態(tài)系統(tǒng)在一定時間內(nèi)向人類所提供的最終產(chǎn)品與服務(wù)等惠益的經(jīng)濟價值總和［1］，包括生態(tài)系統(tǒng)供給服務(wù)價值（生態(tài)系統(tǒng)提供的物質(zhì)產(chǎn)品價值）、調(diào)節(jié)服務(wù)價值和文化服務(wù)價值［2］。研究GEP對開展區(qū)域生態(tài)環(huán)境保護、生態(tài)功能區(qū)劃、環(huán)境經(jīng)濟核算、生態(tài)補償及探索實現(xiàn)生態(tài)產(chǎn)品價值具有重要意義［3］。目前，不少學(xué)者通過能值分析法、當(dāng)量因子法、功能價值法以及模型模擬等方法從不同空間尺度、行政單元開展GEP評估研究。其中，功能價值法因方法成熟，易于理解，符合公眾的認知和判斷，評估結(jié)果可信度高等優(yōu)勢，在中小尺度的生態(tài)系統(tǒng)服務(wù)價值評估中得到廣泛應(yīng)用［4］。

恩施州土壤種類繁多，森林資源豐富，氣候條件優(yōu)越，環(huán)境質(zhì)量優(yōu)良，是我國為數(shù)不多的未大規(guī)模工業(yè)開發(fā)的區(qū)域之一，屬于國家主體功能區(qū)規(guī)劃中的武陵山區(qū)生物多樣性及水土保持生態(tài)功能區(qū)［5］，是長江中上游重要生態(tài)屏障，當(dāng)?shù)厣鷳B(tài)系統(tǒng)在供給優(yōu)質(zhì)農(nóng)林牧漁水資源等生態(tài)物質(zhì)產(chǎn)品的同時，也提供了諸多的生態(tài)系統(tǒng)服務(wù)，例如土壤保持、水源涵養(yǎng)、洪水調(diào)蓄、空氣凈化、水質(zhì)凈化、氣候調(diào)節(jié)、固碳釋氧、景觀游憩、生物保育等，具有重要的生態(tài)、經(jīng)濟和文化價值。根據(jù)李克勤等［6］的研究，2014年恩施州生態(tài)價值（包括生態(tài)資源本身及其所提供的產(chǎn)品和服務(wù)的價值）為10.74萬億元，價值巨大，然而針對恩施州多年的歷史性GEP評估一直缺乏，從而導(dǎo)致對區(qū)域生態(tài)產(chǎn)品理解認識不夠，生態(tài)產(chǎn)品價值開發(fā)實現(xiàn)不足。核算恩施州2011—2020年GEP，利用時空地理加權(quán)回歸（GTWR）模型分析GEP的主要驅(qū)動因子，可對構(gòu)建完善區(qū)域GEP核算評估體系、探索實現(xiàn)生態(tài)產(chǎn)品價值提供一定的科學(xué)參考。

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

1.1 研究區(qū)概況

恩施州位于湖北省西南部武陵山區(qū)，下轄2市6縣8個縣級行政單位，全域24 060 km2，有“八山半水分半田”之稱，如圖1所示。境內(nèi)絕大部分是山地，平均海拔1 000 m，中亞熱帶季風(fēng)氣候區(qū)，年平均氣溫15～17 ℃，年降水量為1 000～1 900 mm；有11個土類236個土種［6］，含硒土壤占全州可利用面積的95.60%，表層土壤硒元素含量背景值0.41 mg/kg，是全球最大天然富硒區(qū)［7］；全州植被覆蓋率達69.23%［8］，有各類自然保護地37處，占全州土地面積12.07%［9］；境內(nèi)河流眾多，5 km以上的河流382條，水能資源理論儲量509萬kW［8］，自然生態(tài)資源豐富。

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

1.2.1 數(shù)據(jù)來源。

采用多源數(shù)據(jù)核算GEP，其中第一、三產(chǎn)業(yè)、農(nóng)林牧漁產(chǎn)品物質(zhì)量、常住人口、旅游收入等相關(guān)數(shù)據(jù)來自《恩施州統(tǒng)計年鑒》（2011—2020年）；恩施州土地資源類型、面積等數(shù)據(jù)來自自然資源部門調(diào)查數(shù)據(jù)；降雨量等氣象數(shù)據(jù)來源于中國氣象數(shù)據(jù)網(wǎng)（http：//data.cma.cn），地表水資源量來自《恩施州水資源公報》（2011—2020年）；自然保護區(qū)、A級景區(qū)、水電裝機等數(shù)據(jù)分別來自林業(yè)、旅游、水利部門的統(tǒng)計數(shù)據(jù)；CPI指數(shù)及農(nóng)林牧漁等價格數(shù)據(jù)來源于中國經(jīng)濟社會大數(shù)據(jù)研究平臺（http：//data.cnki.net/）和國家統(tǒng)計局網(wǎng)站（http：//www.stats.gov.cn/）；所用圖件均來自湖北省地理信息公共服務(wù)平臺（https：//hubei.tianditu.gov.cn/standardMap），審圖號為鄂S（2022）005號，底圖無修改。

1.2.2 變量選取與處理。

人類活動、氣候變化等因素對生態(tài)系統(tǒng)的結(jié)構(gòu)和功能有重大影響，更多情況下，氣候變化與非氣候變化因子同時發(fā)生作用［10-13］?？紤]數(shù)據(jù)可獲得性、指標(biāo)之間不存在共線性且可量化等因素，通過專家打分法，選取對生態(tài)產(chǎn)品價值的影響權(quán)屬較重的地表水資源量、觀光資源量（A級景區(qū)及非A級景區(qū)自然保護地面積）、土壤侵蝕模數(shù)、年降水量、林地面積、耕地面積、水域面積、建設(shè)用地面積、人均GDP、第一產(chǎn)業(yè)總值、旅游總收入、居民消費水平、常住人口數(shù)、城鎮(zhèn)化率14個指標(biāo)作為備選驅(qū)動因子變量。

為避免偽回歸現(xiàn)象，在對所選變量標(biāo)準(zhǔn)化處理后，進行多重共線性檢驗。剔除方差膨脹因子（VIF）大于10的變量，最終選出年降水量、林地面積、人均GDP、城鎮(zhèn)化率4個驅(qū)動因子變量，各變量的VIF分別為：8.276、9.357、1.715、1.746。

1.3 研究方法

1.3.1 生態(tài)產(chǎn)品價值核算方法。

采用功能價值法分步核算恩施州生態(tài)產(chǎn)品價值，在分別核算生態(tài)系統(tǒng)的物質(zhì)供給價值、調(diào)節(jié)服務(wù)價值、文化服務(wù)價值的基礎(chǔ)上，三者求和計算得出生態(tài)產(chǎn)品價值總量。計算公式見表1。

1.3.2 時空地理加權(quán)回歸模型。

時空地理加權(quán)回歸（GTWR）模型在地理加權(quán)回歸（GWR）模型的基礎(chǔ)上引入時間因子，在解決截面數(shù)據(jù)樣本數(shù)量有限的同時，兼顧考慮了時間與空間的非平穩(wěn)性，能有效估計因子參數(shù)［24］，模型如下：

Yi=β0（ui，vi，ti）+pk=1βk（ui，vi，ti）Xik+εi（1）

式中：X與Y分別為解釋變量與被解釋變量，（ui，vi，ti）為第i個樣本單元的時空坐標(biāo)；P為解釋變量個數(shù)；β0（ui，vi，ti）為截距項；βk（ui，vi，ti）為第k個解釋變量的估計系數(shù)；εi為模型殘差。

該研究以ArcGIS10.8軟件為平臺，采用Huang等［25］制作的GTWR插件實現(xiàn)時空地理加權(quán)回歸，其中帶寬采用

AICc優(yōu)化設(shè)置，時空距離參數(shù)比值為1，GEP作為被解釋變量，驅(qū)動因子作為解釋變量，回歸結(jié)果相關(guān)參數(shù)為：

Bindwidth 0.283 2、Sigma 27.931 0、Residual Squqres 62 411.11、AICc 816.25、R2 0.988 3、Adjusted R2 0.987 2、Spatio-temporal Distance Ratio 2.504 6。從擬合優(yōu)度看，R2與矯正后的R2均高于0.98，表明時空地理加權(quán)回歸模型能夠較好地測度解釋變量對被解釋變量的影響。

2 研究結(jié)果

2.1 生態(tài)產(chǎn)品價值年際變化

年際變化能夠較好地反映恩施州生態(tài)產(chǎn)品價值的變化特點，該研究選取2011、2014、2017和2020 4個時間節(jié)點分析GEP的年際變化情況，結(jié)果如表2所示。

2011、2014、2017和2020年恩施州GEP總量分別為 6 055.08億、6 384.79億、7 114.01億和7 417.01億元，逐年增長，2011—2020年的年均增長率為2.50%。2011—2014年，GEP總量年均增長率為1.82%，低于年均增長率；2014—2017年，GEP總量年均增長率為3.81%，高于年均增長率；2017—2020年，GEP總量年均增長率為1.42%，低于年均增長率，且低于2011—2014年的平均增長率。2011—2020年GEP增速呈現(xiàn)先增后減的“A”型趨勢，不同時段增長速度有明顯差異。

2011—2020年，恩施州生態(tài)系統(tǒng)物質(zhì)供給價值、調(diào)節(jié)服務(wù)價值、文化服務(wù)價值的年均增長率分別為12.95%、1.67%、5.97%，物質(zhì)供給價值增速>文化服務(wù)價值增速＞調(diào)節(jié)服務(wù)價值增速。從結(jié)構(gòu)上看，恩施州GEP主要由調(diào)節(jié)服務(wù)提供，占比80%以上，呈現(xiàn)波動下降的態(tài)勢；物質(zhì)供給和文化服務(wù)對恩施州GEP的貢獻較低，占比均在15%以下，呈現(xiàn)波動上升的態(tài)勢。2011—2020年，恩施州GEP內(nèi)部結(jié)構(gòu)呈現(xiàn)出動態(tài)調(diào)整優(yōu)化的特點。從增速上看，2011—2017年文化服務(wù)價值年均增速快于物質(zhì)供給價值和調(diào)節(jié)服務(wù)價值年均增速，物質(zhì)供給價值年均增速呈波動加速狀態(tài)，調(diào)節(jié)服務(wù)價值年均增速相對較慢；2017—2020年文化服務(wù)價值年均增速低于物質(zhì)供給價值，呈現(xiàn)負增長，原因在于新冠疫情導(dǎo)致人員流通減少、旅游經(jīng)濟下滑，恩施州旅游業(yè)及相關(guān)服務(wù)產(chǎn)業(yè)遭受重創(chuàng)，文化服務(wù)價值增長停滯。

2.2 供給水平的空間分異

用單位面積生態(tài)產(chǎn)品價值（即每平方公里生態(tài)產(chǎn)品的價值量，地均生態(tài)產(chǎn)品價值）衡量恩施州各縣市生態(tài)系統(tǒng)的生態(tài)產(chǎn)品產(chǎn)出水平，利用ArcGIS10.8對各縣市不同時間段的地均生態(tài)產(chǎn)品價值（GEP）、地均物質(zhì)供給價值（WZGJ）、地均調(diào)節(jié)服務(wù)價值（TJFW）、地均文化服務(wù)價值（WHFW）進行可視化，并采用自然間斷法分級顯示，如圖2所示。

2011—2020年，各縣市地均生態(tài)產(chǎn)品價值總量、地均物質(zhì)供給價值、地均調(diào)節(jié)服務(wù)價值、地均文化服務(wù)價值逐年上升，生態(tài)產(chǎn)品供給水平呈增長趨勢。各縣市地均生態(tài)產(chǎn)品價值總量在23.46萬～33.85萬元/hm2、地均物質(zhì)供給價值0.86萬～2.52萬元/hm2、地均調(diào)節(jié)服務(wù)價值18.91萬～29.00萬元/hm2、地均文化服務(wù)價值1.60萬～9.60萬元/hm2變化，縣市間生態(tài)產(chǎn)品供給水平存在明顯差異；整體來看，巴東縣、鶴峰縣生態(tài)產(chǎn)品供給水平最高；細分來看，恩施市、宣恩縣、來鳳縣生態(tài)物質(zhì)供給水平相對較高，巴東縣、鶴峰縣、咸豐縣、來鳳縣生態(tài)調(diào)節(jié)服務(wù)供給水平相對較高，恩施市生態(tài)文化服務(wù)供給水平高于其他縣市。

2.3 生態(tài)產(chǎn)品影響因素分析

將GEP作為被解釋變量，年降水量（NJSL）、林地面積（NDMJ）、人均GDP（RJGDP）、城鎮(zhèn)化率（CZHL）作為解釋變量，用GTWR模型分析解釋變量與被解釋變量間的相關(guān)性，對局部樣本可視化，如圖3所示。

總體上看，各縣市不同影響因子的回歸系數(shù)各不相同并具有顯著影響（P＜0.05），且在時間序列上發(fā)生變化；8縣市年降水量、林地面積、人均GDP的回歸系數(shù)均為正數(shù)，對各縣市GEP具有顯著的正向促進效應(yīng)，隨著年降水量、林地面積、人均GDP的增加，各縣市GEP相應(yīng)增大；2011—2020年，城鎮(zhèn)化率因子的回歸系數(shù)由負變正，說明隨著城鎮(zhèn)化率的增加，其對GEP的影響逐漸由負影響變?yōu)檎虼龠M，具有閾值效應(yīng)。對比來看，利川市、咸豐縣、來鳳縣年降水量回歸系數(shù)絕對值相對較大，恩施市、建始縣、巴東縣絕對值相對較?。焕ㄊ?、咸豐縣、來鳳縣林地面積回歸系數(shù)相對較大，恩施市、建始縣、巴東縣回歸系數(shù)相對較??；利川市、恩施市、建始縣人均GDP回歸系數(shù)相對較大，巴東縣、鶴峰縣、來鳳縣的回歸系數(shù)相對較小；2011—2016年8縣市城鎮(zhèn)化率回歸系數(shù)為負數(shù)，2017—2020年8縣市回歸系數(shù)為正數(shù)。

3 討論

2011—2020年，恩施州年平均GEP為6 658.5億元，其中物質(zhì)供給價值占比為5.43%、調(diào)節(jié)服務(wù)價值占比82.31%、文化服務(wù)價值占比12.26%，GEP內(nèi)部結(jié)構(gòu)表現(xiàn)為調(diào)節(jié)服務(wù)價值＞文化服務(wù)價值＞物質(zhì)供給價值，與歐陽志云等［4］在浙江麗水市、趙寅成等［26］在安徽六安市、白瑪卓嘎等［27］在四川甘孜州、牟雪潔等［28］在北京延慶區(qū)、董天等［29］在內(nèi)蒙古鄂爾多斯市的同類研究結(jié)果類似。歐陽志云等［4］提出，根據(jù)GDP與GEP的增長趨勢和關(guān)系，可以評估一個地區(qū)生態(tài)文明建設(shè)進展及人與自然和諧共生的狀態(tài)與趨勢。恩施州GDP由2011年的458.42億元增長到了2020年的1 117.70億元［17］，年均增長率15.98%，GEP與GDP雙增長，GDP增速大于GEP增速，全州社會經(jīng)濟與生態(tài)環(huán)境保護協(xié)調(diào)發(fā)展，生態(tài)文明建設(shè)成效顯著，與2019、2021年恩施州先后獲得國家生態(tài)文明建設(shè)示范州、“綠水青山就是金山銀山”實踐創(chuàng)新基地命名相印證。

該研究得出年降水量、林地面積、人均GDP、城鎮(zhèn)化率因子對恩施州生態(tài)產(chǎn)品價值具有顯著影響，與國內(nèi)其他研究得出的自然氣象、土地利用、經(jīng)濟水平、社會人口等因素對GEP有影響［30-32］的結(jié)論趨于一致。GTWR模型分析結(jié)果顯示，不同時間不同縣市年降水量、林地面積、人均GDP、城鎮(zhèn)化率因子的回歸系數(shù)存在差異，印證了生態(tài)系統(tǒng)服務(wù)驅(qū)動因素及其作用機理的復(fù)雜性，與鄧煌煒等［13］的研究結(jié)論相似。利川市、咸豐縣、來鳳縣年降水量與GEP的相關(guān)性較強，主要原因是上述縣市是恩施州傳統(tǒng)的農(nóng)牧物質(zhì)產(chǎn)品產(chǎn)出大縣，而降水對農(nóng)作物生長和農(nóng)業(yè)生產(chǎn)具有重要作用，間接影響生態(tài)物質(zhì)供給價值。

林地是恩施州最主要的用地類型，提供了大量的林業(yè)產(chǎn)品和調(diào)節(jié)服務(wù)產(chǎn)品，林地面積與GEP總體正相關(guān)，但受林地占比、林相品種、林地等級質(zhì)量的不同，以及天然林保護、退耕還林和石漠化治理等生態(tài)治理差異，回歸系數(shù)在縣市間表現(xiàn)出空間異質(zhì)性。人均GDP作為經(jīng)濟發(fā)展水平的風(fēng)向標(biāo)，回歸系數(shù)均為正，說明經(jīng)濟發(fā)展水平對GEP有正向促進作用，可能有兩個原因，第一是經(jīng)濟發(fā)展水平影響生態(tài)物質(zhì)產(chǎn)品及生態(tài)功能服務(wù)單位價格，第二是經(jīng)濟繁榮與產(chǎn)業(yè)發(fā)展互促，如農(nóng)林牧漁產(chǎn)品的精深加工和旅游文化產(chǎn)業(yè)的繁榮會促進放大生態(tài)產(chǎn)品價值；因經(jīng)濟產(chǎn)業(yè)發(fā)展?fàn)顩r差距，各縣市人均GDP對GEP的影響表現(xiàn)出空間異質(zhì)性。城鎮(zhèn)擴張會減少耕地、林地、水域等能提供重要調(diào)節(jié)服務(wù)功能的生態(tài)用地，從而導(dǎo)致生態(tài)產(chǎn)品價值降低；但城鎮(zhèn)化又會帶來產(chǎn)業(yè)繁榮和經(jīng)濟發(fā)展，對生態(tài)產(chǎn)品價值具有增益放大效應(yīng)；2011—2016年，恩施州城鎮(zhèn)率快速增長，由32.53%增長至41.88%，城鎮(zhèn)化增益效應(yīng)小于城鎮(zhèn)化占用生態(tài)用地負面影響，城鎮(zhèn)化率與GEP負相關(guān)；2017—2020年，恩施州城鎮(zhèn)化率增速變緩，只從41.88%增長至45.56%，城鎮(zhèn)化增益效應(yīng)大于其負面影響，城鎮(zhèn)化率與GEP正相關(guān)；因各縣市發(fā)展基礎(chǔ)、城鎮(zhèn)化路徑和速度不同等原因，城鎮(zhèn)化率回歸系數(shù)表現(xiàn)出時空差異。

4 結(jié)論

利用多源數(shù)據(jù)，結(jié)合統(tǒng)計分析、空間可視化等手段，核算了恩施州2011—2020年的GEP，并運用時空地理加權(quán)回歸模型分析GEP的主要影響因素。具體結(jié)論如下：

（1）2011—2020年，恩施州年平均GEP為6 658.5億元，結(jié)構(gòu)上調(diào)節(jié)服務(wù)價值＞文化服務(wù)價值＞物質(zhì)供給價值；GEP呈增長趨勢，年均增長率2.50%，物質(zhì)供給價值增速＞文化服務(wù)價值增速＞調(diào)節(jié)服務(wù)價值增速；GEP與GDP雙增長，區(qū)域經(jīng)濟社會與生態(tài)環(huán)境保護協(xié)調(diào)發(fā)展。

（2）2011—2020年，恩施州各縣市地均生態(tài)產(chǎn)品價值、地均物質(zhì)供給價值、地均調(diào)節(jié)服務(wù)價值、地均文化服務(wù)價值分別在23.46萬～33.85萬、0.86萬～2.52萬、18.91萬～29.00萬和1.60萬～9.60萬元/hm2間波動變化，且不同縣市間存在明顯差異，具有空間異質(zhì)性。

（3）2011—2020年，恩施州年降水量、林地面積、人均GDP、城鎮(zhèn)化率對GEP具有顯著影響（P＜0.05），不同因子的回歸系數(shù)在縣市間存在空間差異。其中，年降水量、林地面積、人均GDP對GEP正向促進；城鎮(zhèn)化率因子的回歸系數(shù)由負變正，表明對于恩施州這種城鎮(zhèn)化率較低的重點生態(tài)功能區(qū)，隨著城鎮(zhèn)化率的提升，其對GEP的影響逐漸由負面影響變?yōu)檎虼龠M，表現(xiàn)出閾值效應(yīng)。

參考文獻

［1］ 歐陽志云，靳樂山.面向生態(tài)補償?shù)纳鷳B(tài)系統(tǒng)生產(chǎn)總值核算方法［M］.北京：科學(xué)出版社，2017：85.

［2］ 張林波，虞慧怡，李岱青，等.生態(tài)產(chǎn)品內(nèi)涵與其價值實現(xiàn)途徑［J］.農(nóng)業(yè)機械學(xué)報，2019，50（6）：173-183.

［3］ 謝高地，張彩霞，張雷明，等.基于單位面積價值當(dāng)量因子的生態(tài)系統(tǒng)服務(wù)價值化方法改進［J］.自然資源學(xué)報，2015，30（8）：1243-1254.

［4］ 歐陽志云，林亦晴，宋昌素.生態(tài)系統(tǒng)生產(chǎn)總值（GEP）核算研究：以浙江省麗水市為例［J］.環(huán)境與可持續(xù)發(fā)展，2020，45（6）：80-85.

［5］ 湖北省人民政府.湖北省主體功能區(qū)規(guī)劃［Z］. 2012.

［6］ 李克勤，張惠遠.湖北恩施州、十堰市生態(tài)價值研究［M］.北京：中國統(tǒng)計出版社，2017.

［7］ 程水源.《硒與恩施》科普叢書：硒與產(chǎn)業(yè)［M］.長沙：湖南科學(xué)技術(shù)出版社，2016.

［8］ 恩施州統(tǒng)計局，恩施州調(diào)查隊，恩施州調(diào)查監(jiān)測分局.恩施州統(tǒng)計年鑒2020［M］.恩施：恩施州統(tǒng)計局，2020.

［9］ 恩施州林業(yè)局.恩施州自然保護地融合優(yōu)化預(yù)案［Z］. 2020.

［10］ 趙士洞，張永民.千年生態(tài)系統(tǒng)評估報告集（一）［M］.賴鵬飛，譯.北京：中國環(huán)境科學(xué)出版社，2007：B39-B47.

［11］ SUN X，TANG H J，YANG P，et al.Spatiotemporal patterns and drivers of ecosystem service supply and demand across the conterminous United States：A multiscale analysis［J］.Science of the total environment，2020，703：1-17.

［12］ WU X，LIU S L，ZHAO S，et al.Quantification and driving force analysis of ecosystem services supply， demand and balance in China［J］.Science of the total environment，2019，652：1375-1386.

［13］ 鄧煌煒，廖振良.生態(tài)系統(tǒng)服務(wù)驅(qū)動因素及供需研究進展［J］.環(huán)境科技，2020，33（1）：74-78.

［14］ 馬國霞，於方，王金南，等.中國 2015 年陸地生態(tài)系統(tǒng)生產(chǎn)總值核算研究［J］.中國環(huán)境科學(xué)，2017，37（4）：1474-1482.

［15］ 任志遠，李晶.陜南秦巴山區(qū)植被生態(tài)功能的價值測評［J］.地理學(xué)報，2003，58（4）：503-511.

［16］ 景兆鵬，馬友鑫.云南省西雙版納地區(qū)生態(tài)系統(tǒng)服務(wù)價值的動態(tài)評估［J］.中南林業(yè)科技大學(xué)學(xué)報，2012，32（9）：87-93.

［17］ 趙同謙，歐陽志云，王效科，等.中國陸地地表水生態(tài)系統(tǒng)服務(wù)功能及其生態(tài)經(jīng)濟價值評價［J］.自然資源學(xué)報，2003，18（4）：443-452.

［18］ 從日征，王兵，牛香，等.陜西省森林生態(tài)系統(tǒng)凈化大氣環(huán)境功能價值評估［J］.西北林學(xué)院學(xué)報，2017，32（5）：75-82.

［19］ 姚鑫，楊桂山.自然濕地水質(zhì)凈化研究進展［J］.地理科學(xué)進展，2009，28（5）：825-832.

［20］ 王莉雁，肖燚，歐陽志云，等.國家級重點生態(tài)功能區(qū)縣生態(tài)系統(tǒng)生產(chǎn)總值核算研究：以阿爾山市為例［J］.中國人口·資源與環(huán)境，2017，27（3）：146-154.

［21］ 武文婷，夏國元，包志毅.杭州市城市綠地固碳釋氧價值量評估［J］.中國園林，2016，32（3）：117-121.

［22］ 張曉利，馬力，魯小珍，等.游憩價值評價方法探討：以鳳陽山自然保護區(qū)為例［J］.中國人口·資源與環(huán)境，2011，21（S1）：213-216.

［23］ 王兵，鄭秋紅，郭浩.基于Shannon-Wiener指數(shù)的中國森林物種多樣性保育價值評估方法［J］.林業(yè)科學(xué)研究，2008，21（2）：268-274.

［24］ 沈楊，汪聰聰，高超，等.基于城市化的浙江省灣區(qū)經(jīng)濟帶碳排放時空分布特征及影響因素分析［J］.自然資源學(xué)報，2020，35（2）：329-342.

［25］ HUANG B，WU B，BARRY M.Geographically and temporally weighted regression for modeling spatio-temporal variation in house prices［J］.International journal of geographical information science，2010，24（3）：383-401.

［26］ 趙寅成，孫雷，岳正波，等.安徽省六安市生態(tài)系統(tǒng)生產(chǎn)總值核算研究［J］.安徽農(nóng)業(yè)科學(xué)，2021，49（4）：73-77.

［27］ 白瑪卓嘎，肖燚，歐陽志云，等. 甘孜藏族自治州生態(tài)系統(tǒng)生產(chǎn)總值核算研究［J］.生態(tài)學(xué)報，2017，37（19）：6302-6312.

［28］ 牟雪潔，王夏暉，張簫，等. 北京市延慶區(qū)生態(tài)系統(tǒng)生產(chǎn)總值核算及空間化［J］.水土保持研究，2020，27（1）：265-274，282.

［29］ 董天，張路，肖燚，等. 鄂爾多斯市生態(tài)資產(chǎn)和生態(tài)系統(tǒng)生產(chǎn)總值評估［J］.生態(tài)學(xué)報，2019，39（9）：3062-3074.

［30］ 趙曉冏，蘇軍德，王建，等.甘肅省生態(tài)服務(wù)供需關(guān)系及影響因子研究［J］.中國環(huán)境科學(xué)，2021，41（10）：4926-4941.

［31］ 虎陳霞，郭旭東，連綱，等.長三角快速城市化地區(qū)土地利用變化對生態(tài)系統(tǒng)服務(wù)價值的影響：以嘉興市為例［J］.長江流域資源與環(huán)境，2017，26（3）：333-340.

［32］ 姚成勝，朱鶴健，呂晞，等.土地利用變化的社會經(jīng)濟驅(qū)動因子對福建生態(tài)系統(tǒng)服務(wù)價值的影響［J］.自然資源學(xué)報，2009，24（2）：225-233.