生態(tài)文明視閥下中國區(qū)域GDP與ESV的協(xié)同關(guān)系研究

[關(guān)鍵詞] 生態(tài)系統(tǒng)服務(wù)價值；似不相關(guān)回歸模型；時空演進(jìn)；協(xié)同關(guān)系；驅(qū)動因素

［中圖分類號］F127［文獻(xiàn)標(biāo)志碼］A[文章編號］1008-0694（2025）03－0057－18

[基金項目]重慶市哲學(xué)社會科學(xué)創(chuàng)新工程重點(diǎn)項目“生態(tài)產(chǎn)品價值實現(xiàn)機(jī)制研究”（2025CXZD18）；重慶市統(tǒng)計局2025年度統(tǒng)計科研重點(diǎn)項目“重慶市經(jīng)濟(jì)韌性評價研究”（2025ky03 ）；重慶市北碚區(qū)統(tǒng)計局第五次全國經(jīng)濟(jì)普查課題研究重點(diǎn)項目“北碚區(qū)現(xiàn)代化產(chǎn)業(yè)體系的構(gòu)建與測度研究”。

黨的二十大報告指出，要推進(jìn)城鄉(xiāng)融合和區(qū)域協(xié)調(diào)發(fā)展，建立生態(tài)產(chǎn)品價值實現(xiàn)機(jī)制。目前，我國提出了一系列生態(tài)空間類型并制定了不同的管控措施，多數(shù)地區(qū)開始出臺地方核算標(biāo)準(zhǔn)以及核算技術(shù)指南，不僅可以為生態(tài)系統(tǒng)功能提供保障，也為經(jīng)濟(jì)高質(zhì)量發(fā)展提供綠色支撐。但是，隨著我國城鎮(zhèn)化和工業(yè)化的快速發(fā)展，人口集聚、國土開發(fā)建設(shè)等活動加劇，國內(nèi)生產(chǎn)總值（GDP）增長導(dǎo)向型的發(fā)展模式最終造成一定程度的自然環(huán)境污染、生態(tài)系統(tǒng)退化和資源過度消耗等問題，阻礙了生態(tài)和經(jīng)濟(jì)的協(xié)調(diào)可持續(xù)發(fā)展。因此，如何在追求經(jīng)濟(jì)發(fā)展與保護(hù)生態(tài)環(huán)境之間取得最佳平衡已成為一個十分關(guān)鍵且亟待解決的問題，如何推動經(jīng)濟(jì)社會發(fā)展全面綠色轉(zhuǎn)型以及生態(tài)環(huán)境質(zhì)量改善由量變向質(zhì)變也成為各級政府的工作重點(diǎn)與社會關(guān)注的焦點(diǎn)。

一、文獻(xiàn)綜述

現(xiàn)有文獻(xiàn)關(guān)于生態(tài)產(chǎn)品價值的核算方法主要有三種：一是當(dāng)量因子法，1-2]用于核算生態(tài)系統(tǒng)服務(wù)價值（EcosystemServiceValue，ESV），表示通過生態(tài)系統(tǒng)功能直接或間接獲得產(chǎn)品與服務(wù)的價值。20世紀(jì)末，國外學(xué)者Costanza等（1997）、3]Daily（1997）就開始了關(guān)于ESV的研究。4此后，生態(tài)系統(tǒng)服務(wù)價值評估成為國內(nèi)外相關(guān)學(xué)科領(lǐng)域研究的熱點(diǎn)問題。許多學(xué)者也在全球、國家、區(qū)域和流域等不同尺度上計算了ESV。[5-8]二是功能價值法，9]用于核算生態(tài)系統(tǒng)生產(chǎn)總值（Gross EcosystemProduct，GEP），表示生態(tài)系統(tǒng)為人類福祉和經(jīng)濟(jì)社會可持續(xù)發(fā)展提供的最終產(chǎn)品與服務(wù)價值的總和。國內(nèi)學(xué)者開展GEP核算已涉及全球、全國、省域、市域等不同區(qū)域尺度。10-13]三是能值轉(zhuǎn)換法，[14-15]用于將生態(tài)系統(tǒng)中不同類別的能量和物質(zhì)轉(zhuǎn)換為統(tǒng)一標(biāo)準(zhǔn)的能值以進(jìn)行定量分析。但該方法在使用中需標(biāo)明能值基線，在生態(tài)產(chǎn)品價值核算上的應(yīng)用較少。相對功能價值法而言，當(dāng)量因子法的核算體系較為規(guī)范，評價方法和參數(shù)標(biāo)準(zhǔn)統(tǒng)一，運(yùn)算過程較為簡單，數(shù)據(jù)要求也較低。此外，謝高地等（2015）在Costanza等的研究基礎(chǔ)上對反映我國不同生態(tài)系統(tǒng)單位面積ESV的當(dāng)量因子表進(jìn)行了更新和完善，16特別適用于區(qū)域尺度ESV的評估。[17]

國內(nèi)外學(xué)者對生態(tài)系統(tǒng)與經(jīng)濟(jì)系統(tǒng)之間的研究不斷豐富完善。一是理論基礎(chǔ)方面，美國經(jīng)濟(jì)學(xué)家Grossman等（1995）首次提出環(huán)境庫茲涅茨曲線（Environmen-talKuznetsCurve，EKC），指出了環(huán)境污染與經(jīng)濟(jì)發(fā)展存在先增后減的“倒U型”曲線關(guān)系。18]二是研究方法方面，大多數(shù)學(xué)者從地理角度或基于統(tǒng)計數(shù)據(jù)構(gòu)建生態(tài)系統(tǒng)指標(biāo)評價體系并計算生態(tài)產(chǎn)品價值以研究生態(tài)環(huán)境與經(jīng)濟(jì)發(fā)展相互作用的演化，[19-20]較少學(xué)者采用當(dāng)量因子法[21]和功能價值法進(jìn)行研究。[22-23]三是研究內(nèi)容方面，主要包括生態(tài)和經(jīng)濟(jì)的時空變化特征與影響因素、[24-25]耦合協(xié)調(diào)關(guān)系與發(fā)展模式、協(xié)調(diào)發(fā)展路徑等，26-27]并從全國、省際、市域、地區(qū)、流域等多種空間尺度進(jìn)行分析。28-32]目前，關(guān)于生態(tài)環(huán)境與經(jīng)濟(jì)發(fā)展的時空演變、耦合協(xié)調(diào)關(guān)系的研究多集中于農(nóng)業(yè)經(jīng)濟(jì)、旅游經(jīng)濟(jì)和土地利用等方面。33]

綜上所述，學(xué)界對生態(tài)系統(tǒng)與經(jīng)濟(jì)系統(tǒng)的相關(guān)研究比較成熟，但仍存在以下不足：一是現(xiàn)有文獻(xiàn)主要集中在流域以及城市群等大尺度研究區(qū)域，對市、縣以及重要生態(tài)功能區(qū)等小尺度區(qū)域的生態(tài)產(chǎn)品價值核算與研究不夠豐富。二是研究內(nèi)容大多集中于時空演變及耦合協(xié)調(diào)等單方面研究，對區(qū)域空間協(xié)同關(guān)系及其時空分異特征的研究不足，也缺乏對影響GDP與ESV的驅(qū)動因素識別，且解釋變量多集中在社會經(jīng)濟(jì)因素。三是從實證方法看，運(yùn)用空間計量模型進(jìn)行驅(qū)動因素分析的文章尚不多見。本文的貢獻(xiàn)和可能存在的創(chuàng)新點(diǎn)主要有以下三點(diǎn)：一是將時空演進(jìn)、協(xié)同關(guān)系和驅(qū)動因素三者結(jié)合在一起進(jìn)行研究，更全面地揭示經(jīng)濟(jì)發(fā)展和生態(tài)保護(hù)之間的內(nèi)在關(guān)聯(lián)。二是驅(qū)動因素指標(biāo)的選取全面考慮了科學(xué)技術(shù)和生態(tài)環(huán)境等自然地理環(huán)境和社會經(jīng)濟(jì)發(fā)展因素，估計結(jié)果更加準(zhǔn)確。三是采用空間計量模型即似不相關(guān)回歸模型分析驅(qū)動因素，并進(jìn)行時空異質(zhì)性分析，結(jié)果更具針對性。

二、研究區(qū)域與研究方法

1．研究區(qū)域概況

重慶市位于我國西南部，下轄38個區(qū)縣。作為長江流域重要生態(tài)屏障，重慶在水土保持、水源涵養(yǎng)等方面發(fā)揮著關(guān)鍵作用。但是，重慶市又作為生態(tài)環(huán)境高度敏感區(qū)，生態(tài)資源比較脆弱，其經(jīng)濟(jì)的快速發(fā)展帶來了一系列生態(tài)問題，如水土流失、生態(tài)承載力降低、石漠化等，34-36城市建設(shè)用地的擴(kuò)張也給生態(tài)帶來了負(fù)面影響。7]因此，深入研究重慶市生態(tài)系統(tǒng)服務(wù)價值與國內(nèi)生產(chǎn)總值的時空關(guān)系，尤其是研究重慶“一區(qū)兩群”的GDP與ESV變化特征，①對于推動重慶區(qū)域協(xié)調(diào)可持續(xù)發(fā)展具有重要的現(xiàn)實意義。

2．?dāng)?shù)據(jù)來源

本研究的樣本為重慶市38個區(qū)縣，研究時期為2000、2010、2020年，所用數(shù)據(jù)主要包括土地利用數(shù)據(jù)、GDP數(shù)據(jù)和驅(qū)動因素數(shù)據(jù)：（1）土地利用數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心。（2）GDP數(shù)據(jù)來源于相應(yīng)年份的《重慶統(tǒng)計年鑒》。（3）驅(qū)動因素涉及的社會經(jīng)濟(jì)、人口結(jié)構(gòu)、生態(tài)環(huán)境等數(shù)據(jù)來源于相應(yīng)年份的《重慶統(tǒng)計年鑒》、相關(guān)部門的官方網(wǎng)站以及官方媒體公開披露的統(tǒng)計年報，缺失的數(shù)據(jù)通過區(qū)縣統(tǒng)計年鑒、國民經(jīng)濟(jì)和社會發(fā)展統(tǒng)計公報以及EPS數(shù)據(jù)庫進(jìn)行補(bǔ)充；年均氣溫、年降水量數(shù)據(jù)來自中國氣象數(shù)據(jù)網(wǎng)；地勢起伏度所用的高程數(shù)據(jù)來自地理空間數(shù)據(jù)云；歸一化植被指數(shù)（NormalizedDifferenceVegetation In-dex，NDVI）數(shù)據(jù)通過地理空間數(shù)據(jù)云下載。為了增加各年份間的可比性，本文以2020 年為基期通過國內(nèi)生產(chǎn)總值指數(shù)對GDP進(jìn)行平減處理。

3．研究方法

（1）生態(tài)系統(tǒng)服務(wù)價值核算。本研究采用當(dāng)量因子法核算生態(tài)系統(tǒng)服務(wù)價值，根據(jù)重慶市實際情況，參考周李磊等（2023）的研究成果得到重慶市不同生態(tài)系統(tǒng)服務(wù)價值修正系數(shù)。38]由于建設(shè)用地對生態(tài)環(huán)境的影響具有積累效應(yīng)，難以量化其ESV，因此，只將耕地、林地、草地、水域和未利用地五種土地利用類型納入研究中。

式（1）和式（2）中： ESV_i 、ESV分別為第i種生態(tài)系統(tǒng)服務(wù)的價值和總的生態(tài)系統(tǒng)服務(wù)價值； A_i 為第 i 種土地利用類型的面積， VC_i 為第i種生態(tài)系統(tǒng)服務(wù)的價值系數(shù)。

（2）耦合協(xié)調(diào)分析。Tapio解耦模型適用于分析兩種指標(biāo)之間的相關(guān)性，已有學(xué)者將類似的模型用于研究經(jīng)濟(jì)增長與能源消耗、39]區(qū)域經(jīng)濟(jì)與生態(tài)產(chǎn)品價值之間的關(guān)系。40]基于該模型，本文利用 ESV 和GDP 的變化引人耦合指數(shù)以描述重慶市生態(tài)保護(hù)與經(jīng)濟(jì)發(fā)展之間的協(xié)同關(guān)系，具體公式如下：

式（3）、式（4）和式（5）中：e代表ESV與GDP之間的耦合指數(shù)， GDP_pr 代表GDP的變化率， ESV_pr 代表ESV的變化率，pi和pj分別為研究期的開始和結(jié)束年份，pr為研究期間，即pi－pj。為了區(qū)分不同的耦合協(xié)調(diào)類型，參照Guan 等（2022）的研究，41將計算結(jié)果分為八種類型，如表1所示。但由于重慶市GDP 從2000年以來處于增長狀態(tài)，①故不予考慮 ΔGDPlt;0 的情況。

表1耦合協(xié)調(diào)類型

（3）空間自相關(guān)分析。現(xiàn)有研究多采用單變量Moran'sI指數(shù)來描述空間自相關(guān)關(guān)系，雙變量自相關(guān)分析是對傳統(tǒng)方法的進(jìn)一步改進(jìn)，其可以揭示兩個變量之間的空間相關(guān)性。本研究采用雙變量Moran'sI分析GDP與ESV之間的空間集聚和離散模式，即采用全局空間自相關(guān)法（Moran'sI指數(shù)）檢測重慶市GDP與ESV的空間自相關(guān)，采用局部二元空間自相關(guān)法（LISA聚類圖）檢測區(qū)縣間的空間自相關(guān)。全局 ，小于0表示GDP與ESV呈負(fù)相關(guān)，大于0則表示正相關(guān)，如果接近0，則表明存在離散分布。具體計算公式如下：

式（6）、式（7）和式（8）中： W_ij 為空間權(quán)重值， I₁ 為雙變量的全局Moran'sI指數(shù)， I₂ 為雙變量的局部Moran'sI指數(shù)， n 為空間單元個數(shù)， Z_xi 與 Z_yj 分別代表兩個變量在空間單元i、j的標(biāo)準(zhǔn)化值。

（4）驅(qū)動因素分析。除了自然本身的初始資源稟賦，人類活動也會對經(jīng)濟(jì)系統(tǒng)和生態(tài)系統(tǒng)造成影響。結(jié)合重慶市實際狀況，主要從社會經(jīng)濟(jì)、技術(shù)水平、人口結(jié)構(gòu)、自然資源、生態(tài)環(huán)境和農(nóng)業(yè)狀況6個方面選取驅(qū)動因素以考察自然地理環(huán)境因素和社會經(jīng)濟(jì)發(fā)展因素對GDP與ESV協(xié)調(diào)發(fā)展的影響（詳見表2），其中，對經(jīng)濟(jì)發(fā)展水平、投資力度、城鎮(zhèn)居民收入、農(nóng)村居民收入和地方財政支出5個因素作對數(shù)處理。

考慮到解釋變量之間的同期相關(guān)性，即存在某未知因素可能同時影響重慶市GDP與ESV，并且各方程的隨機(jī)誤差項之間也可能存在同期相關(guān)性。采用似不相關(guān)（SUR）聯(lián)立方程進(jìn)行樣本估計，該模型得到的參數(shù)估計值方差更小，可以顯著提升整體估計效率，是分析驅(qū)動因素的重要方法。42-43]具體模型設(shè)定如下：

式（9）和式（10）中，i代表區(qū)縣； 、 lnESV_i，t 為被解釋變量，分別表示第i個區(qū)縣的GDP與ESV水平， X_{k，i，t} 表示第 k 個解釋變量， α₀ 和 β₀ 為截距項，α_k 和 β_k 為待估參數(shù)， 1 u_i，t 表示方程的隨機(jī)誤差項。

表2GDP與ESV的驅(qū)動因素指標(biāo)表

三、GDP與ESV的時空演進(jìn)特征與協(xié)同關(guān)系分析

1．GDP與ESV時空演進(jìn)特征分析

（1）重慶市GDP時空演變特征分析。如表3所示，①從GDP整體空間分布格局來看，剔除物價因素，研究期間重慶市經(jīng)濟(jì)發(fā)展空間不均衡，重慶的西部區(qū)域GDP 明顯高于東部區(qū)域，主要原因在于西部區(qū)域區(qū)位條件優(yōu)越，經(jīng)濟(jì)基礎(chǔ)雄厚，吸引并聚集了較多的高素質(zhì)人才，而東部區(qū)域交通運(yùn)輸不發(fā)達(dá)，人口外流嚴(yán)重，經(jīng)濟(jì)基礎(chǔ)較為薄弱。從地均GDP來看，低值區(qū)主要分布在渝東南和渝東北地區(qū)且數(shù)值均低于全市平均水平，中心城區(qū)的土地面積最少且經(jīng)濟(jì)密度遠(yuǎn)高于全市平均水平。從人均GDP來看，全市呈現(xiàn)出西高東低的特征，僅中心城區(qū)的人均GDP遠(yuǎn)高于全市平均水平。整體而言，重慶市存在著明顯的空間分異特征，GDP數(shù)值總體呈“西高東低，一區(qū)兩群”的分布格局，且高低格局基本穩(wěn)定，初步形成了以渝中和萬州為核心的兩個高值區(qū)。

表3重慶市“一區(qū)兩群”GDP和ESV

（2）重慶市ESV時空演變特征分析。土地利用變化分析。如表4所示，用地類型面積反映耕地和林地是重慶市主要的生態(tài)系統(tǒng)。近20年來，重慶市耕地和草地等生態(tài)系統(tǒng)區(qū)域有所減少，林地、水域以及建設(shè)用地等生態(tài)系統(tǒng)面積有所增加。用途的改變引起了用地類型的變化，各類生態(tài)系統(tǒng)面積的增減表明人類活動對生態(tài)的影響已經(jīng)加強(qiáng)，并且占用了部分農(nóng)田生態(tài)系統(tǒng)，導(dǎo)致耕地等生態(tài)系統(tǒng)面積減少以及服務(wù)價值受損。而ESV的增長得益于林地和水域等生態(tài)系統(tǒng)面積的增多，這與重慶市積極實施天然林資源保護(hù)、退耕還林還草等林業(yè)重點(diǎn)工程密切相關(guān)，張艷軍等（2017）也指出三峽以及水利工程的大規(guī)模建設(shè)提高了水域 ESV［63]用地類型ESV顯示，2020年生態(tài)系統(tǒng)產(chǎn)生的ESV大小為：林地 gt; 耕地 gt; 草地 gt; 水域 gt; 未利用地。耕地和草地出現(xiàn)不同程度的下降，主要原因可能是部分居民生態(tài)意識薄弱，不適當(dāng)?shù)目撤ヒ约伴_發(fā)建設(shè)等活動致使土地退化或資源被占用。

表4重慶市生態(tài)系統(tǒng)用地類型ESV

生態(tài)系統(tǒng)服務(wù)類型變化分析。如表5所示，2000—2020年重慶市ESV逐年提高，實際增長率為 4.48% 。生態(tài)系統(tǒng)調(diào)節(jié)服務(wù)是ESV的核心價值，2020年占比達(dá)到 67.94% 。并且近20年間生態(tài)系統(tǒng)調(diào)節(jié)服務(wù)價值增加了227.01億元，表明提供改善人類生存與生活條件的調(diào)節(jié)服務(wù)水平不斷提升。但重慶市ESV及各個服務(wù)類型ESV的增長速度較低且均呈減弱態(tài)勢，其中，調(diào)節(jié)服務(wù)價值增速下降幅度最高，表明人類活動給重慶市生態(tài)系統(tǒng)帶來負(fù)面影響。在供給服務(wù)價值上，雖然生態(tài)產(chǎn)品供給能力逐漸增強(qiáng)，保障了人們?nèi)找嬖鲩L的物質(zhì)需要，但從整體上看，生態(tài)系統(tǒng)對生態(tài)環(huán)境的正向作用力越來越弱。

表5重慶市生態(tài)系統(tǒng)服務(wù)類型ESV①

區(qū)縣生態(tài)系統(tǒng)服務(wù)價值分析。如表3所示，②從空間分異特征來看，重慶市ESV總體發(fā)展態(tài)勢良好，但區(qū)域差異性顯著，總體呈現(xiàn)“東北、東南高，西部低”的空間格局，具有明顯的“一區(qū)兩群”分布特征，局部呈現(xiàn)高一高集聚和低—低集聚的分布格局。具體而言，ESV高值區(qū)主要分布在渝東北和渝東南地區(qū)，且渝東北ESV高于渝東南，該類地區(qū)擁有較多的山脈，自然資源豐富，具有經(jīng)濟(jì)水平落后和城鎮(zhèn)化水平較低的特點(diǎn)，人類活動對生態(tài)系統(tǒng)造成的負(fù)面影響較小，并且也得益于近年來的石漠化治理、水土保持對生態(tài)環(huán)境建設(shè)和保護(hù)的作用，其生態(tài)系統(tǒng)穩(wěn)定。低值區(qū)主要分布在中心城區(qū)和西部的主城新區(qū)，中心城區(qū)的ESV遠(yuǎn)遠(yuǎn)低于其他片區(qū)，該類地區(qū)處于低山丘陵地帶，社會經(jīng)濟(jì)的快速發(fā)展擾動了區(qū)域生態(tài)系統(tǒng)，并且居民建設(shè)用地向外擴(kuò)張、城市建設(shè)規(guī)模集中連片等也導(dǎo)致ESV偏低。從時間來看，只有渝東南和渝東北地區(qū)空間格局依次發(fā)生變化，ESV進(jìn)一步提升。

從地均ESV來看，主要呈現(xiàn)“西低東高，一區(qū)兩群”的分布格局，且空間集聚程度不斷加深。另外，運(yùn)用趨勢面法分析地均ESV的分異趨勢發(fā)現(xiàn)（詳見圖1），擬合曲線變化差異不大，從南北方向看，各地區(qū)地均ESV從北至南呈現(xiàn)“U型”趨勢，即渝東南、渝東北地區(qū)地均ESV較高；從東西方向看，各地區(qū)地均ESV從西到東逐漸上升，并且在2010 年表現(xiàn)出明顯的“西部下降、東部上升”特點(diǎn)。就人均ESV值而言，人口密度是造成渝東南地區(qū)和中心城區(qū)人均ESV之間顯著差異的主要原因，說明高強(qiáng)度的人類社會經(jīng)濟(jì)活動導(dǎo)致生態(tài)系統(tǒng)受到劇烈干擾，在一定程度上降低了該區(qū)域的ESV水平。

圖1重慶市2000年、2010年和2020年地均ESV趨勢面演化情況

2.GDP與ESV的協(xié)同關(guān)系分析

如表6所示，在整個研究期間，重慶市耦合指數(shù)為0.51，生態(tài)經(jīng)濟(jì)處于低度失調(diào)發(fā)展?fàn)顟B(tài)。但e值由2000—2010年的0.56增長至2010—2020年的0.75，兩者耦合協(xié)調(diào)的階段性特征顯著，即重慶市全面實施林業(yè)等生態(tài)重點(diǎn)工程，推動生態(tài)恢復(fù)，生態(tài)經(jīng)濟(jì)協(xié)調(diào)關(guān)系得到明顯改善。但這兩段時期重慶生態(tài)經(jīng)濟(jì)均處于低度失調(diào)狀態(tài)，表明研究期間生態(tài)保護(hù)與經(jīng)濟(jì)發(fā)展之間的矛盾仍然存在。

表6 重慶市GDP與ESV的耦合指數(shù)

GDP與ESV的協(xié)同關(guān)系具有明顯的時空差異。如表7所示，重慶市生態(tài)經(jīng)濟(jì)協(xié)同發(fā)展存在3種類型，高度協(xié)調(diào)、低度協(xié)調(diào)和低度失調(diào)。從生態(tài)經(jīng)濟(jì)協(xié)調(diào)類型時間演變過程來看，2000—2010年協(xié)調(diào)類型以低度失調(diào)為主，占全部區(qū)縣數(shù)的81.58% ；2010—2020年重慶市協(xié)調(diào)類型變化較為顯著，低度失調(diào)區(qū)縣數(shù)占全部區(qū)縣數(shù)的 63.16% ，數(shù)量顯著下降可能與重慶大力實施生態(tài)環(huán)境保護(hù)政策息息相關(guān)。其中，分布于主城新區(qū)和渝東北地區(qū)的豐都縣、合川區(qū)、江津區(qū)等區(qū)縣由低度協(xié)調(diào)轉(zhuǎn)為低度失調(diào)狀態(tài)，表明這些區(qū)縣在開始快速發(fā)展經(jīng)濟(jì)的同時對生態(tài)環(huán)境的平衡造成了破壞。而有14個區(qū)縣的生態(tài)與經(jīng)濟(jì)正逐步形成良性互動關(guān)系，由失調(diào)轉(zhuǎn)變?yōu)閰f(xié)調(diào)狀態(tài)，其中，城口縣、大渡口區(qū)、南岸區(qū)和萬州區(qū)實現(xiàn)跨階段轉(zhuǎn)變，由低度失調(diào)轉(zhuǎn)變?yōu)楦叨葏f(xié)調(diào)狀態(tài)，表明這4個區(qū)縣在發(fā)展經(jīng)濟(jì)的同時也開始高度重視生態(tài)保護(hù)。整個研究期間的協(xié)調(diào)程度不高，絕大部分區(qū)縣處于低度失調(diào)狀態(tài)。

表7重慶市協(xié)調(diào)類型區(qū)縣數(shù)量情況

結(jié)合空間演化過程來看，本文以耦合指數(shù)1為分界線，將協(xié)調(diào)發(fā)展特征劃分為生態(tài)相對滯后與經(jīng)濟(jì)相對滯后，耦合指數(shù)范圍在0.8—1.2的為同步發(fā)展，使用ArcGIS10.8繪制的空間分布圖詳見圖2。2000—2010年，重慶市絕大多數(shù)區(qū)縣耦合指數(shù)小于0.89，整體協(xié)調(diào)發(fā)展特征表現(xiàn)為生態(tài)相對滯后。生態(tài)相對發(fā)展較好的地區(qū)集中在主城新區(qū)以及渝東北地區(qū)，渝北區(qū)和萬州區(qū)生態(tài)相對滯后特征最為明顯，這兩個區(qū)也是經(jīng)濟(jì)發(fā)展的重點(diǎn)城市。2010—2020年，中心城區(qū)的發(fā)展比較符合EKC的變化趨勢，由生態(tài)相對滯后型向同步發(fā)展轉(zhuǎn)變，西部區(qū)域的主城新區(qū)進(jìn)入EKC上行區(qū)間，由同步發(fā)展向生態(tài)相對滯后型轉(zhuǎn)變。此外，生態(tài)相對發(fā)展較好的地區(qū)向“兩群”方向移動，并表現(xiàn)為經(jīng)濟(jì)相對滯后，該現(xiàn)象與程靜等（2022）提出的耦合指數(shù)具有空間溢出效應(yīng)結(jié)論相符。4]但在2000—2020 年間，重慶市GDP與 ESV的協(xié)同關(guān)系并不強(qiáng)，為絕對性的低度失調(diào)且偏向于生態(tài)相對滯后型發(fā)展。

圖2重慶市耦合指數(shù)空間分布情況

四、GDP與ESV協(xié)同關(guān)系的驅(qū)動因素分析

1．GDP與ESV空間自相關(guān)分析

采用GeoDa1.20軟件進(jìn)行空間自相關(guān)分析，以人均GDP、地均ESV分別表征不同區(qū)域間的經(jīng)濟(jì)發(fā)展水平、生態(tài)環(huán)境狀況。如圖3所示，通過建立空間權(quán)重矩陣，計算出2000年、2010 年和2020 年重慶市全局空間自相關(guān)指數(shù)Moran'sI值分別為-0.465、-0.455和-0.521（ Plt;0.001 ），表明研究期內(nèi)重慶市人均GDP與地均ESV呈顯著空間負(fù)相關(guān)關(guān)系且總體關(guān)系逐漸增強(qiáng)。

圖3 重慶市人均GDP與地均ESVMoran'sI指數(shù)散點(diǎn)圖

本文進(jìn)一步利用LISA聚類圖分析局部空間自相關(guān)，如圖4所示，2000—2020年重慶市人均GDP與地均ESV的局部空間格局具有顯著相似性。在 95% 的置信度水平下可以清晰地將重慶市各區(qū)縣分為3種局部自相關(guān)類型。具體表現(xiàn)為：（1）高一低集聚區(qū)，代表經(jīng)濟(jì)發(fā)展水平較高的區(qū)域會抑制鄰域的生態(tài)環(huán)境狀況，代表區(qū)縣有北碚區(qū)、江北區(qū)等，且該類集聚隨時間推移在重慶市主城都市區(qū)進(jìn)一步擴(kuò)散。（2）低一低集聚區(qū)，代表經(jīng)濟(jì)發(fā)展水平和鄰域生態(tài)的發(fā)展均處于劣勢，但雙低的數(shù)量在研究期間顯著減少，代表區(qū)縣有潼南區(qū)。（3）低一高集聚區(qū)，代表經(jīng)濟(jì)發(fā)展水平相對落后但鄰域擁有較好的生態(tài)環(huán)境狀況，典型的如巫溪、酉陽等渝東北、渝東南的區(qū)縣，且近20年來未發(fā)生較大變化。綜合來看，重慶市西部區(qū)域主要以高一低組合為主，東部區(qū)域以低一高組合為主，總體發(fā)展的均衡性協(xié)調(diào)性有待提高，并且局部空間自相關(guān)結(jié)果表明區(qū)縣經(jīng)濟(jì)與生態(tài)發(fā)展的深層次互動會對鄰域產(chǎn)生一定的輻射作用，也再次證明耦合協(xié)調(diào)度存在明顯的空間溢出效應(yīng)。

圖4重慶市人均GDP與地均ESV雙變量局部LISA聚類圖

2．驅(qū)動因素分析

（1）全樣本估計。如表8所示，B-P檢驗的P值小于0.01，說明選擇SUR方法可以提高估計效率。參考連玉君等（2017）的做法，5]先除去面板數(shù)據(jù)中的個體效應(yīng)，將變換后的數(shù)據(jù)視為截面數(shù)據(jù)，再進(jìn)行SUR估計和組間系數(shù)檢驗。由于本文樣本數(shù)量不大，因此使用了小樣本條件的方法，并進(jìn)行了聚類調(diào)整得到穩(wěn)健標(biāo)準(zhǔn)誤。

限于篇幅，表中未展示參數(shù)的聚類穩(wěn)健標(biāo)準(zhǔn)誤，僅顯示參數(shù)的顯著性P值。如表8所示，在 5% 的顯著性水平下，經(jīng)濟(jì)發(fā)展水平和年降水量為正向驅(qū)動效應(yīng)，即兩者數(shù)值增加會進(jìn)一步提升生態(tài)環(huán)境與經(jīng)濟(jì)發(fā)展的協(xié)調(diào)性，且系數(shù)均在 1% 的水平上存在組間顯著差異。年降水量不僅會影響生態(tài)系統(tǒng)的水資源供給、水文調(diào)節(jié)和生物多樣性等功能，也有利于植被生長和生態(tài)系統(tǒng)修復(fù)，其對ESV的驅(qū)動作用顯著高于GDP。產(chǎn)業(yè)結(jié)構(gòu)為負(fù)向驅(qū)動效應(yīng)，且組間系數(shù)沒有顯著差異，表明第二產(chǎn)業(yè)產(chǎn)值增加在一定程度上意味著工業(yè)企業(yè)的大幅度擴(kuò)張，粗放式的工業(yè)化發(fā)展模式不僅需要消耗大量資源，也會排放廢棄物對生態(tài)造成破壞，從而抑制環(huán)境與經(jīng)濟(jì)的協(xié)調(diào)關(guān)系。此外，張俊等（2023）指出工業(yè)企業(yè)向郊區(qū)和鄉(xiāng)村轉(zhuǎn)移也會對生態(tài)系統(tǒng)服務(wù)產(chǎn)生擠出效應(yīng)。城鎮(zhèn)化率對GDP和ESV的作用呈“U型”的非線性關(guān)系，且組間系數(shù)差異不顯著，可能是由于早期部分區(qū)縣過度依賴土地出讓收人和土地抵押融資推進(jìn)城鎮(zhèn)化建設(shè)，浪費(fèi)了大量生態(tài)用地資源，城鄉(xiāng)建設(shè)用地的擴(kuò)張也侵占和破壞了部分生態(tài)用地，導(dǎo)致GDP和ESV降低，但隨著以人為核心的新型城鎮(zhèn)化與生態(tài)文明建設(shè)戰(zhàn)略的實施，推動了GDP和ESV協(xié)調(diào)發(fā)展。

表8全樣本回歸結(jié)果

注： ①^*** 、**、*分別表示參數(shù)估計在 1% 、 5% 及 10% 的臨界水平上統(tǒng)計顯著，下同； ② Breusch -Pagan 服從卡方分布； ③ 組間系數(shù)差異檢驗的P值顯著則代表組間系數(shù)存在顯著差異。

在 10% 的顯著性水平下，地方財政支出對GDP和ESV有正向驅(qū)動作用，意味著包含有節(jié)能環(huán)保支出的財政投人對生態(tài)經(jīng)濟(jì)的協(xié)調(diào)關(guān)系產(chǎn)生積極作用。人口密度會提高GDP但降低ESV水平，即人口密度的提升雖然會給經(jīng)濟(jì)發(fā)展帶來人力資源和規(guī)模，但也會減弱生態(tài)系統(tǒng)承載力，并且該指標(biāo)越高往往意味著人類活動強(qiáng)度越大，會擠占部分生態(tài)系統(tǒng)空間，并通過非農(nóng)活動消耗更多資源，阻礙生態(tài)與經(jīng)濟(jì)的協(xié)調(diào)發(fā)展。加入人口密度的平方項發(fā)現(xiàn)，①人口密度對GDP和ESV的驅(qū)動作用呈近似的“倒U型”關(guān)系，即人口密度在超過一定水平后會對GDP和ESV造成負(fù)向驅(qū)動，系數(shù)也存在組間顯著差異。結(jié)合前文分析發(fā)現(xiàn)，人口密度對ESV的制約作用更強(qiáng)，生態(tài)環(huán)境更脆弱。

（2）異質(zhì)性分析。年份異質(zhì)性分析。本文考察各研究時期驅(qū)動因素對GDP 和ESV的異質(zhì)性影響。結(jié)果顯示，地方財政支出會對GDP和ESV產(chǎn)生正向驅(qū)動作用且系數(shù)在總體上逐漸增大。人口密度對ESV的影響始終為負(fù)，但隨著時間推移，人口密度對ESV的負(fù)向驅(qū)動效應(yīng)逐漸減弱。年均氣溫對GDP產(chǎn)生正向影響而對ESV產(chǎn)生負(fù)向影響，且隨時間變化呈現(xiàn)“正負(fù)向效果減弱”的態(tài)勢。具體而言，氣溫升高會通過改變生物多樣性和生態(tài)系統(tǒng)功能而對生態(tài)系統(tǒng)產(chǎn)生較大的破壞效應(yīng)。]城鎮(zhèn)化水平僅在2020 年對GDP 和ESV 產(chǎn)生顯著影響，且對 ESV 的負(fù)向驅(qū)動作用增強(qiáng)。

區(qū)域異質(zhì)性分析。驅(qū)動因素對區(qū)域的GDP和ESV發(fā)展有著較為明顯的空間異質(zhì)性。經(jīng)濟(jì)發(fā)展水平在主城都市區(qū)和渝東南地區(qū)對GDP和ESV的協(xié)調(diào)發(fā)展具有正向驅(qū)動作用，對渝東南地區(qū)的影響效果更大。人口密度對渝東北地區(qū)GDP、渝東南地區(qū)ESV的正向驅(qū)動作用均高于主城都市區(qū)。地方財政支出僅對渝東北地區(qū)和渝東南地區(qū)的ESV產(chǎn)生正向驅(qū)動作用。在主城都市區(qū)和渝東北地區(qū)，城鎮(zhèn)化水平對GDP和ESV的驅(qū)動作用呈近似“U型”的非線性關(guān)系，渝東北地區(qū)的城鎮(zhèn)化水平系數(shù)更大，且僅對渝東北地區(qū)ESV的影響處于上升階段。在渝東南地區(qū)，城鎮(zhèn)化水平對GDP具有正向驅(qū)動作用，對ESV的影響為“U型”關(guān)系。

結(jié)合數(shù)據(jù)的可得性，在社會經(jīng)濟(jì)和自然資源等因素的基礎(chǔ)上，本研究對2020年重慶市增加技術(shù)水平、生態(tài)環(huán)境和農(nóng)業(yè)狀況3個方面的指標(biāo)數(shù)據(jù)展開進(jìn)一步分析，如表9所示。首先，從社會經(jīng)濟(jì)方面看，城鎮(zhèn)居民收人對ESV具有負(fù)向驅(qū)動作用，而農(nóng)村居民收入對GDP與ESV具有正向驅(qū)動作用，因此，有必要提升農(nóng)村居民收入水平，縮小城鄉(xiāng)居民收入差距。其次，從技術(shù)水平方面看，Ramp;D經(jīng)費(fèi)投入強(qiáng)度正向驅(qū)動GDP與ESV的協(xié)調(diào)發(fā)展。再次，從自然資源方面看，歸一化植被指數(shù)對ESV有著積極影響，且組間系數(shù)差異顯著。最后，從農(nóng)業(yè)狀況方面看，農(nóng)作物播種面積對ESV的正向驅(qū)動作用更明顯。總體而言，社會經(jīng)濟(jì)等因素對GDP和ESV的影響比自然環(huán)境等因素更大。

表92020年重慶市驅(qū)動因素分析結(jié)果

五、結(jié)論與啟示

本文測算了重慶市2000—2020年三個階段的生態(tài)系統(tǒng)服務(wù)價值，分析了GDP和ESV的演化規(guī)律，綜合運(yùn)用Tapio解耦模型、空間自相關(guān)模型和似不相關(guān)回歸模型等方法，深人剖析了GDP和ESV的協(xié)同關(guān)系以及影響二者關(guān)系的驅(qū)動因素，其主要研究結(jié)論如下：

第一，重慶市ESV不斷提高，但ESV增長速度較低且逐漸減弱，生態(tài)系統(tǒng)調(diào)節(jié)服務(wù)是其核心價值。耕地和草地生態(tài)系統(tǒng)的面積和ESV值呈不同程度的下降。因此，要因地制宜加快區(qū)域?qū)嵤┥鷳B(tài)保護(hù)和修復(fù)項目，強(qiáng)化國土空間規(guī)劃和用途管制，嚴(yán)守耕地保護(hù)紅線，增強(qiáng)生態(tài)系統(tǒng)循環(huán)能力和穩(wěn)定性。

第二，重慶市GDP呈“西高東低”的分布格局，ESV呈“東北、東南高，西部低”的空間格局，均符合“一區(qū)兩群”分布特征。GDP與ESV的負(fù)相關(guān)性逐漸增強(qiáng)，空間集聚特征明顯，重慶市的西部、東部區(qū)域分別以高一低集聚、低一高集聚為主。因此，應(yīng)基于ESV建立生態(tài)財政轉(zhuǎn)移支付政策，從橫向、縱向兩個維度健全市場化、多元化的生態(tài)補(bǔ)償機(jī)制，加大綠色金融支持力度，以綠色金融改革創(chuàng)新試驗區(qū)開展金融支持生態(tài)產(chǎn)品價值實現(xiàn)機(jī)制試點(diǎn)。

第三，GDP與ESV的協(xié)同關(guān)系有著明顯的時空差異。重慶市協(xié)調(diào)類型以低度失調(diào)為主，發(fā)展特征整體表現(xiàn)為生態(tài)相對滯后，但耦合指數(shù)上升，生態(tài)經(jīng)濟(jì)協(xié)同關(guān)系趨于良性。因此，區(qū)縣間要打破行政壁壘，搭建開放共享的綠色發(fā)展成果平臺與生態(tài)信息云平臺，增強(qiáng)跨區(qū)域生態(tài)保護(hù)的聯(lián)合調(diào)控能力。同時，基于協(xié)調(diào)度的溢出效應(yīng)，深化區(qū)域之間的合作，以成渝地區(qū)雙城經(jīng)濟(jì)圈建設(shè)為引領(lǐng)，加快產(chǎn)業(yè)轉(zhuǎn)型升級，推動人才、技術(shù)等綠色發(fā)展要素雙向靈活流動。

第四，經(jīng)濟(jì)發(fā)展水平、年降水量和地方財政支出有利于GDP與ESV的協(xié)調(diào)發(fā)展，產(chǎn)業(yè)結(jié)構(gòu)具有抑制作用。從異質(zhì)性分析來看，地方財政支出的驅(qū)動作用進(jìn)一步加強(qiáng)，經(jīng)濟(jì)發(fā)展水平對渝東南地區(qū)的正向驅(qū)動效應(yīng)最大，城鎮(zhèn)化水平對渝東北地區(qū)ESV的影響正處于上升階段。進(jìn)一步分析表明，縮小城鄉(xiāng)居民收入差距和增加Ramp;D經(jīng)費(fèi)投人正向驅(qū)動GDP與ESV的協(xié)調(diào)發(fā)展。因此，政府要加大節(jié)能減排工作力度，加大對區(qū)域的環(huán)保財政投人力度，建立和落實地方污染治理指導(dǎo)等相關(guān)配套政策和財政激勵機(jī)制，并且要繼續(xù)完善重慶市碳交易平臺政策，動態(tài)優(yōu)化配額分配機(jī)制。同時，要縮小城鄉(xiāng)居民收人差距，提高產(chǎn)業(yè)鏈增值收益，大力發(fā)展鄉(xiāng)村特色產(chǎn)業(yè)。此外，要依據(jù)生態(tài)系統(tǒng)承載力合理控制區(qū)域人口密度，通過城鎮(zhèn)化將人口轉(zhuǎn)移出生態(tài)環(huán)境脆弱地區(qū)以減輕對生態(tài)的負(fù)面影響。

參考文獻(xiàn)：

[1][3] COSTANZA R，D'ARGE R，DE GROOT R， et al. The value of the world's ecosys-tem services and natural capital [J]. Nature，1997，（387）.

[2］謝高地，魯春霞，冷允法，等．青藏高原生態(tài)資產(chǎn)的價值評估［J]．自然資源學(xué)報，2003，（02）.

[4]DAILY G C. Nature's services： Societal dependence on natural ecosystems [M]. Washing-ton DC：Island Press，1997：220-221.

[5] COSTANZA R，DE GROOT R， SUTTON P，et al. Changes in the global value of ecosys-tem services[J].Global Environmental Change，2O14，（26）.

[6] ABDULLAH S， ADNAN M S G，BARUA D， et al. Urban green and blue space changes：A spatiotemporal evaluation of impacts on ecosystem service value in Bangladesh [J].Ecological Infor-matics，2022，（70）.

[7]［45］［50］［51］朱玉鑫，姚順波．基于生態(tài)系統(tǒng)服務(wù)價值變化的環(huán)境與經(jīng)濟(jì)協(xié)調(diào)發(fā)展研究—以陜西省為例［J]．生態(tài)學(xué)報，2021，（09）.

[8] ZHANG S，WANG Y，WANG Y，et al. Spatiotemporal evolution and influencing mecha-nisms of ecosystem service value in the Tarim River Basin，northwest China [J]. Remote Sensing，2023，（03）.

[9]歐陽志云，朱春全，楊廣斌，等．生態(tài)系統(tǒng)生產(chǎn)總值核算：概念、核算方法與案例研究［J]．生態(tài)學(xué)報，2013，（21）.

[10]JIANG HQ，WU W J，WANG JN，et al. Mapping global value of terrestrial ecosystem services by countries [J]. Ecosystem Services，2O21，（11）.

[11]OUYANG Z Y， ZHENG H，XIAO Y，et al. Improvements in ecosystem services from in-vestments in natural capital [J]． Science，2O16，（6292）.

[12］白楊，李暉，王曉媛，等．云南省生態(tài)資產(chǎn)與生態(tài)系統(tǒng)生產(chǎn)總值核算體系研究［J].自然資源學(xué)報，2017，（07）.

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

[14] ODUM H T. Environmental accounting： Emergy and environmental decision making [M].New York： John Wiley amp; Sons，1996：1-380.

[15］喻鋒，李曉波，王宏，等．基于能值分析和生態(tài)用地分類的中國生態(tài)系統(tǒng)生產(chǎn)總值核算研究［J]．生態(tài)學(xué)報，2016，（06）.

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

[17］［48］［49］［62］鄧偉，周渝，張勇，等．重慶市生態(tài)保護(hù)紅線區(qū)生態(tài)系統(tǒng)服務(wù)價值時空演變特征及其驅(qū)動［J]．長江流域資源與環(huán)境，2020，（01）.

[18]GROSSMAN G M，KRUEGER A B. Economic growth and the environment [J]. The Quarterly Journal Of Economics，1995，（02）.

[19］吳玉鳴，張燕．中國區(qū)域經(jīng)濟(jì)增長與環(huán)境的耦合協(xié)調(diào)發(fā)展研究［J］．資源科學(xué)，2008，（01）.

[20］［31］［33］王芳．基于耦合協(xié)調(diào)度模型的生態(tài)系統(tǒng)與經(jīng)濟(jì)系統(tǒng)協(xié)同發(fā)展研究——以京津冀地區(qū)為例［J]．湖北社會科學(xué)，2021，（06）.

[21] CHEN W X， ZENG J， ZHONG M X， et al. Coupling analysis of ecosystem services value and economic development in the Yangtze River Economic Belt： A case study in Hunan Province，Chi-na [J].Remote Sensing，2021，（08）.

[22］[30］韓增林，趙玉青，閆曉露，等．生態(tài)系統(tǒng)生產(chǎn)總值與區(qū)域經(jīng)濟(jì)耦合協(xié)調(diào)機(jī)制及協(xié)同發(fā)展—以大連市為例［J]．經(jīng)濟(jì)地理，2020，（10）.

[23][29][41] GUAN S，LIAO Q，WU W J，et al. Revealing the coupling relationship be-tween the gross ecosystem product and economic growth： A case study of Hubei Province [J]. Sustain-ability，2022，（13）.

[24］曹躍群，趙世寬，楊玉玲，等．重慶市生態(tài)系統(tǒng)服務(wù)價值與區(qū)域經(jīng)濟(jì)增長的時空動態(tài)關(guān)系研究［J]．長江流域資源與環(huán)境，2020，（11）.

[25][47] FU J， ZHANG Q，WANG P，et al. Spatio - temporal changes in ecosystem service value and its coordinated development with economy： A case study in Hainan Province，China [J].Remote Sensing，2022，（04）.

[26］張雅杰，馬國創(chuàng)．荊州市域縣市區(qū)生態(tài)經(jīng)濟(jì)耦合協(xié)調(diào)關(guān)系與發(fā)展模式研究［J]．生態(tài)經(jīng)濟(jì)，2016，（12）.

[27］隋建利，陳豪．生態(tài)足跡視域下環(huán)境與經(jīng)濟(jì)增長協(xié)調(diào)發(fā)展路徑研究［J]．財貿(mào)經(jīng)濟(jì)，2021，（06）.

[28][6O] XU M Y，CHEN C T，DENG X Y. Systematic analysis of the coordination degree ofChina's economy - ecological environment system and its influencing factor [J]. Environmental Sci-ence and Pollution Research，2019，（29）.

[32]張中浩，聶甜甜，高陽，等．長江經(jīng)濟(jì)帶生態(tài)系統(tǒng)服務(wù)與經(jīng)濟(jì)社會發(fā)展耦合協(xié)調(diào)關(guān)聯(lián)時空特征研究［J]．長江流域資源與環(huán)境，2022，（05）.

[34］張可，朱康文，雷波，等．2005—2015年三峽庫區(qū)（重慶段）水土流失敏感性時空分異研究［J]．昆明理工大學(xué)學(xué)報（自然科學(xué)版），2020，（03）.

[35]周健，官冬杰，周李磊．基于生態(tài)足跡的三峽庫區(qū)重慶段后續(xù)發(fā)展生態(tài)補(bǔ)償標(biāo)準(zhǔn)量化研究［J]．環(huán)境科學(xué)學(xué)報，2018，（11）.

[36］魏興萍，楊華．重慶巖溶地區(qū)石漠化分布與地理環(huán)境因素的關(guān)系［J］．重慶師范大學(xué)學(xué)報（自然科學(xué)版），2014，（05）.

[37］張騫，高明，楊樂，等．1988—2013年重慶市主城九區(qū)生態(tài)用地空間結(jié)構(gòu)及其生態(tài)系統(tǒng)服務(wù)價值變化［J]．生態(tài)學(xué)報，2017，（02）.

[38］周李磊，蘇湘媛，向洪莉，等．重慶市生態(tài)系統(tǒng)服務(wù)與人類福祉耦合關(guān)系模擬［J].中國環(huán)境科學(xué)，2023，（05）.

[39]DONG B， ZHANG M，MU HL，et al. Study on decoupling analysis between energy con-sumption and economic growth in Liaoning Province [J]. Energy Policy，2O16，（97）.

[40］魏偉，石培基，魏曉旭，等．中國陸地經(jīng)濟(jì)與生態(tài)環(huán)境協(xié)調(diào)發(fā)展的空間演變［J].生態(tài)學(xué)報，2018，（08）.

[42]楊倫，劉某承，閔慶文，等．哈尼梯田地區(qū)農(nóng)戶糧食作物種植結(jié)構(gòu)及驅(qū)動力分析[J]．自然資源學(xué)報，2017，（01）.

[43］劉自敏，馬靚靚，張婭．中國清潔與非清潔能源的替代彈性估計及結(jié)構(gòu)優(yōu)化研究[J]．環(huán)境經(jīng)濟(jì)研究，2023，（01）.

[44][46]［59]LI HM，JIANG Z M，DONG G H，et al. Spatiotemporal coupling coordina-tion analysis of social economy and resource environment of central cities in the Yelow River Basin[J].Discrete Dynamics in Nature and Society，2021，（01）.

[52]［61］呂德勝，王玨，程振．黃河流域數(shù)字經(jīng)濟(jì)、生態(tài)保護(hù)與高質(zhì)量發(fā)展時空耦合及其驅(qū)動因素［J]．經(jīng)濟(jì)問題探索，2022，（08）.

[53]［54］［55］［56］［57］［58］胡雪．山西省生態(tài)經(jīng)濟(jì)系統(tǒng)協(xié)調(diào)度演變格局及其影響因素［J]．水土保持通報，2023，（02）.

[63］張艷軍，官冬杰，翟俊，等．重慶市生態(tài)系統(tǒng)服務(wù)功能價值時空變化研究［J]．環(huán)境科學(xué)學(xué)報，2017，（03）.

[64］程靜，陳紅翔，韓永貴．三江源地區(qū)生態(tài)經(jīng)濟(jì)系統(tǒng)協(xié)調(diào)度評價及其時空演變特征[J]．水土保持研究，2022，（06）.

[65］連玉君，廖俊平．如何檢驗分組回歸后的組間系數(shù)差異？［J]．鄭州航空工業(yè)管理學(xué)院學(xué)報，2017，（06）.

[66]張俊，汪輝，馮越珺．農(nóng)業(yè)生態(tài)系統(tǒng)服務(wù)價值評價及其驅(qū)動因素：基于動態(tài)調(diào)節(jié)當(dāng)量的實證考察[J]．財貿(mào)研究，2023，（04）.

[67]GARCIA F C，BESTION E，WARFIELD R，et al. Changes in temperature alter the relationship between biodiversity and ecosystem functioning [J]. Proceedings of the National Academy of Sciences，2018，（43）.

（責(zé)任編輯 張筠）

Abstract： Under the background of ecological civilization construction，it is of great significance to explore the relationship between GDP and ecosystem services value（ESV）for regional green and sustainable development. Based on the ecosystem services value data of 38 districts and counties in Chongqing from 2O0O to 2O2O，the spatial-temporal evolution characteristics of GDP and ESV were analyzed，and the synergistic relationship between GDP and ESV and its driving factors were discussed using Tapio decoupling model and seemingly unrelated regresson model. The study found that ： （1）The ESV of Chongqing constantly improved，but the growth rate of ESV was low and gradually declined，and the ecosystem regulation service was its core value.（2） The spatial differentiation characteristics of GDP and ESV were obvious，and the negative correlation between regional GDP and ESV gradually strengthened，with distinct spatial agglomeration features.（3）The synergistic relationship between GDP and ESV exhibited obvious spatial -temporal diffrences，and the overal coordination degree was not high，but the coupling index increased，and the synergistic relationship betwen ecology and economy tended to be benign. （4） The impacts of different driving factors on the coordinated development of GDP and ESV showed significant diferences and spatial-temporal heterogeneity.Economic development level，annual precipitation，and local fiscal expenditure were conducive to the coordinated development of GDP and ESV，while industrial structure had an inhibitory effct.Based on these findings，policy recommendations were proposed from the aspects of providing ecological compensation， building a service -oriented government，and promoting the construction of urban-rural ecological civilization.

Keywords： Ecosystem Services Value； Seemingly Unrelated Regression Model; Spatial - Tempo-ral Evolution； Synergistic Relationship；Driving Factors