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

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

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

中圖分類號(hào) F124.5? 文獻(xiàn)標(biāo)識(shí)碼 A

文章編號(hào) 0517-6611（2024）06-0068-07

doi：10.3969/j.issn.0517-6611.2024.06.016

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

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

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

1.1 研究區(qū)概況

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

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

1.2.1 數(shù)據(jù)來(lái)源。

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

1.2.2 變量選取與處理。

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

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

1.3 研究方法

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

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

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

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

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

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

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

AICc優(yōu)化設(shè)置，時(shí)空距離參數(shù)比值為1，GEP作為被解釋變量，驅(qū)動(dòng)因子作為解釋變量，回歸結(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，表明時(shí)空地理加權(quán)回歸模型能夠較好地測(cè)度解釋變量對(duì)被解釋變量的影響。

2 研究結(jié)果

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

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

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

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

2.2 供給水平的空間分異

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

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

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

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

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

3 討論

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

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

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

4 結(jié)論

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

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

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

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

參考文獻(xiàn)

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

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

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

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

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

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

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

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

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

［10］ 趙士洞，張永民.千年生態(tài)系統(tǒng)評(píng)估報(bào)告集（一）［M］.賴鵬飛，譯.北京：中國(guó)環(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ū)動(dòng)因素及供需研究進(jìn)展［J］.環(huán)境科技，2020，33（1）：74-78.

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

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

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

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

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

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

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

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

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

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

［24］ 沈楊，汪聰聰，高超，等.基于城市化的浙江省灣區(qū)經(jīng)濟(jì)帶碳排放時(shí)空分布特征及影響因素分析［J］.自然資源學(xué)報(bào)，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］ 白瑪卓嘎，肖燚，歐陽(yáng)志云，等. 甘孜藏族自治州生態(tài)系統(tǒng)生產(chǎn)總值核算研究［J］.生態(tài)學(xué)報(bào)，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)總值評(píng)估［J］.生態(tài)學(xué)報(bào)，2019，39（9）：3062-3074.

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

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

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