趙桂梅 耿涌 孫華平

摘要 供給側(cè)改革和高質(zhì)量發(fā)展的適時提出為中國低碳經(jīng)濟發(fā)展提供了新思路，調(diào)整以往過多的行政配置資源帶來的要素扭曲，研究同一個國家不同省份之間資源要素協(xié)同問題的方法可以為提出有效的節(jié)能減排政策提供科學依據(jù)與決策參考?；诖耍瑢ξ覈寂欧艔姸瓤臻g面板數(shù)據(jù)進行測算，提出探索性時空數(shù)據(jù)分析框架，運用GeoDa軟件刻畫我國省際碳排放強度時空關(guān)聯(lián)及其局部空間差異的波動性，將STIRPAT模型和EKC模型相結(jié)合檢驗我國碳排放強度影響因素的空間溢出特性，揭示我國碳排放強度時空演進的一般規(guī)律。研究結(jié)果表明：我國碳排放強度的局部空間自相關(guān)呈現(xiàn)集聚與分異并存的時空分布特征，碳排放強度水平相近的省區(qū)均呈現(xiàn)出先減弱后增強的空間集聚現(xiàn)象，我國省際碳排放強度存在俱樂部收斂的演進特征。我國省際碳排放強度存在空間溢出性特征，我國碳排放強度空間集聚、分異與躍遷的形成與演變是各地區(qū)資源稟賦以及經(jīng)濟發(fā)展水平等多種資源要素協(xié)同作用的結(jié)果。進一步提出分區(qū)域、有重點的區(qū)域協(xié)同減排措施，提升區(qū)域碳排放協(xié)同治理的效率。

關(guān)鍵詞 高質(zhì)量發(fā)展;碳排放強度;空間效應(yīng);傳導機制;協(xié)同減排

碳排放強度指標反映的是經(jīng)濟發(fā)展過程中的資源利用效率以及碳排放效率，碳排放強度指標的高低在一定程度上可以反映一個國家或者地區(qū)的生產(chǎn)技術(shù)效率水平，基于碳排放強度的節(jié)能減排政策有利于形成推動中國經(jīng)濟轉(zhuǎn)型的倒逼機制和長效機制[1-3]。碳排放強度的準確評估是建立碳減排政策的基礎(chǔ)與前提，因此，國家和地區(qū)碳減排政策的制定需要重點關(guān)注各省區(qū)資源要素協(xié)同，降低碳排放強度是推動綠色低碳轉(zhuǎn)型，實現(xiàn)高質(zhì)量經(jīng)濟發(fā)展的重要抓手。《巴黎協(xié)定》明確了2℃溫升控制目標，發(fā)出了全球綠色低碳轉(zhuǎn)型的明確信號。由此，有針對性地剖析中國省際碳排放強度的空間效應(yīng)及其傳導機制，強化區(qū)域能源和氣候政策的精準性、協(xié)同性和實質(zhì)性作用力，不僅是實現(xiàn)碳排放總量達峰目標的核心任務(wù)，更是保障中國經(jīng)濟高質(zhì)量可持續(xù)發(fā)展的必然要求[4-5]。

1? 文獻回顧

目前，國內(nèi)外學者研究文獻主要集中在對中國碳排放量和碳排放強度影響因素進行分解[6-8]，然而，在建立碳排放影響因素的分解模型時，僅僅關(guān)注到各影響因素對本地區(qū)碳排放變動產(chǎn)生的直接影響，常常忽略了各區(qū)域之間存在的空間關(guān)聯(lián)性，沒有深入研究碳排放影響因素的外溢效應(yīng)以及輻射作用[9-10]。隨著計量經(jīng)濟學的發(fā)展，學者們開始運用空間計量經(jīng)濟模型對碳排放總量、人均碳排放、碳排放強度等變動的影響因素展開分析[11-12]。其中，最為廣泛的是對碳排放總量的研究，通過構(gòu)建碳排放空間計量模型，發(fā)現(xiàn)中國碳排放量的梯度分布對地區(qū)之間碳排放空間集聚產(chǎn)生了強化作用[13-14]，進一步采用空間計量經(jīng)濟學的理論與方法揭示中國省域碳減排時空格局的演變機理。這也為后續(xù)對碳排放空間交互作用機制研究提供了理論支持。

由空間集聚所引致的空間效率與區(qū)域公平之間的權(quán)衡問題是經(jīng)濟地理學的前沿研究主題之一。根據(jù)“地理學第一定律”的原則，空間分布中的地理事物或?qū)傩韵嗷リP(guān)聯(lián)，相鄰區(qū)域之間的聯(lián)系更為密切[15]。探索性空間時空數(shù)據(jù)分析（ESTDA）是以空間權(quán)重矩陣為基礎(chǔ)，檢測和量化不同校準點中取樣變化的空間依賴性的統(tǒng)計方法，ESTDA方法能夠有效地揭示任何地理單元與鄰近單元之間的空間關(guān)系以及關(guān)聯(lián)程度。經(jīng)濟發(fā)展質(zhì)量的提高是經(jīng)濟發(fā)展方式轉(zhuǎn)變的前提，其中政府的環(huán)境治理有助于提升中國經(jīng)濟的高質(zhì)量發(fā)展[16]。同時研究發(fā)現(xiàn)制度質(zhì)量和綠色技術(shù)創(chuàng)新都會對能源效率提升具有重要的影響[17]。根據(jù)新結(jié)構(gòu)經(jīng)濟學的分析，有效的市場和有為的政府共同發(fā)揮作用，才能構(gòu)建起經(jīng)濟持續(xù)發(fā)展的良性機制[18]。供給側(cè)改革和高質(zhì)量發(fā)展的適時提出為中國低碳經(jīng)濟發(fā)展提供了新思路，調(diào)整以往過多的行政配置資源帶來的要素扭曲，研究同一個國家不同省份之間資源要素協(xié)同問題的方法可以為提出真實有效的節(jié)能減排政策提供科學依據(jù)與決策參考。

2? 研究方法及數(shù)據(jù)來源

2.1 探索性時空數(shù)據(jù)分析方法

全局空間自相關(guān)性分析是對區(qū)域之間空間自相關(guān)性的總體度量[19]，通常由Z值檢驗和全局Morans I指數(shù)進行判斷各省區(qū)碳排放強度是否存在空間自相關(guān)性，以及其影響程度和作用方向，計算公式如下：

2.2? 空間計量模型設(shè)定

STIRPAT模型既允許將各系數(shù)作為參數(shù)進行估計，又允許對各影響因子進行適當?shù)姆纸夂透倪M[20]。鑒于STIRPAT模型和EKC假說是環(huán)境污染影響研究的基本理論框架，為了研究資源要素空間外溢效應(yīng)在不同區(qū)域的異質(zhì)性，該研究將兩者相結(jié)合來檢驗各省區(qū)資源要素對碳排放空間效應(yīng)的影響。結(jié)合經(jīng)濟系統(tǒng)與碳排放系統(tǒng)的復(fù)雜性和相關(guān)性[21-22]，在實際回歸模型中，除了包括人口規(guī)模（P）、經(jīng)濟發(fā)展水平（GDPPC）、能源強度（EI）等核心解釋變量外，還加入了能源價格（EP）、對外開放水平（OPEN）、能源消費結(jié)構(gòu)（ECS）、產(chǎn)業(yè)結(jié)構(gòu)（IS）和城鎮(zhèn)化水平（UR）等控制變量，將STRIPAT模型擴展為：

其中，i， t分別表示觀測樣本和時間，β1it～β8it為相應(yīng)回歸元系數(shù)， αit為常數(shù)項。為了更客觀地衡量中國碳排放強度空間相關(guān)性和異質(zhì)性，提高空間計量經(jīng)濟模型的可解釋性，在區(qū)域?qū)用娼⒖臻g面板數(shù)據(jù)模型時，結(jié)合了地理空間權(quán)重矩陣和經(jīng)濟權(quán)重矩陣的點積。利用反比大圓距離法構(gòu)造地理空間權(quán)重矩陣，用經(jīng)度和緯度來計算的大圓距離倒數(shù)來度量空間關(guān)系。在STIRPAT擴展模型的基礎(chǔ)上，建立了以下三類空間面板數(shù)據(jù)模型。

該研究以我國30個省區(qū)作為研究對象（由于數(shù)據(jù)限制，香港、澳門、臺灣和西藏不在研究范圍之列），通過構(gòu)造各省區(qū)平均能源消耗的碳排放系數(shù)測算各省區(qū)碳排放總量和碳排放強度，數(shù)據(jù)來源于歷年的《中國統(tǒng)計年鑒》《中國人口統(tǒng)計年鑒》《中國能源統(tǒng)計年鑒》《中國城市（鎮(zhèn)）生活與價格年鑒》《中國人口與就業(yè)統(tǒng)計年鑒》《中國科技統(tǒng)計年鑒》《中國工業(yè)統(tǒng)計年鑒》及各區(qū)域的統(tǒng)計年鑒。

3? 實證分析

3.1? 碳排放空間效應(yīng)分析

全局Morans I指數(shù)從整體上描述我國各省區(qū)碳排放強度的空間相關(guān)性，時間維度的引入可以體現(xiàn)碳排放強度的空間動態(tài)特征。運用GeoDa軟件生成一次鄰接規(guī)則的空間權(quán)重矩陣（Rook），計算出我國碳排放強度的Morans I值，并采用隨機排列法構(gòu)造正態(tài)分布檢驗其顯著性，進一步將全局Morans I指數(shù)分解到各個省區(qū)單元，根據(jù)全局Morans I指數(shù)分析中各階段起始年份同時考慮時間間隔的均衡性，采用空間計量軟件GeoDa095i，繪制主要年份碳排放強度的局域 Morans I散點圖，反映我國30個省區(qū)碳排放強度與各自相鄰近省區(qū)碳排放強度在空間分布上呈現(xiàn)“集聚”（高-高或低-低） 與“分異”（高-低或低-高） 的時空演進特征。如圖1所示，省際碳排放強度存在較強的空間正相關(guān)關(guān)系，各個省區(qū)碳排放強度對與之相鄰地區(qū)的碳排放強度產(chǎn)生影響，同時，其自身的碳排放強度也會受相鄰地區(qū)碳排放強度的影響。我國各省區(qū)碳排放強度的空間分布是集聚和擴散共存的，研究期間呈現(xiàn)空間正相關(guān)性（高-高和低-低）省市的個數(shù)全部大于19個，表現(xiàn)為聚集特征的省區(qū)已經(jīng)超過63.33%，說明我國碳排放強度的時空集聚特征顯著，通過對6個年份的分析結(jié)果也驗證了全局空間相關(guān)性的分析。進一步分析發(fā)現(xiàn)，大多數(shù)省區(qū)處于第一和第三象限，具有較強的穩(wěn)健性，表明各省及鄰近地區(qū)碳排放強度呈正相關(guān)，這些省區(qū)與相鄰接的省區(qū)體現(xiàn)出相似的特征值，證明了各省域的碳排放強度在空間上存在較為顯著的相互依賴性，我國碳排放強度相近省區(qū)的空間集聚表現(xiàn)出不斷增強的趨勢，碳排放強度水平相近的省區(qū)也將會進一步收斂?？v向比較發(fā)現(xiàn)位于第三象限的“低-低”省區(qū)數(shù)量遠遠多于位于第一象限“高-高”省區(qū)，表明低值聚集數(shù)量多于高值集聚區(qū)，是正向空間自相關(guān)性的主要貢獻部分，表明我國近年來推行的節(jié)能減排和產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化政策取得一定成效。

3.2 碳排放空間效應(yīng)的傳導機制分析

為了更全面地反映空間效應(yīng)，首先對資源要素碳排放空間效應(yīng)影響的空間固定效應(yīng)、時間固定效應(yīng)、時空雙固定效應(yīng)以及無固定效應(yīng)的4類傳統(tǒng)面板數(shù)據(jù)模型進行LM檢驗，結(jié)果如表1所示。四種類型空間效應(yīng)模型均通過0.05顯著性水平的空間滯后LM檢驗，其中，時空雙固定效應(yīng)模型的空間誤差LM檢驗結(jié)果最為顯著，進一步對SLPDM和SEPDM的穩(wěn)健性進行LM檢驗，結(jié)果顯示拒絕非空間性的判斷。從時空固定效應(yīng)LR檢驗結(jié)果也表明拒絕無時間固定效應(yīng)以及無空間固定效應(yīng)的原假設(shè)，因此，LR檢驗結(jié)果顯示時空雙固定效應(yīng)模型能夠更加準確地能刻畫資源要素空間溢出效應(yīng)。

基于LM檢驗結(jié)果顯示拒絕非空間模型時，就應(yīng)該慎重選擇SLPDM或者SEPDM，進一步考慮包括SLPDM和SEPDM的空間杜賓面板數(shù)據(jù)模型（SDPDM），對于SDPDM進行Wald檢驗和LR檢驗，從而判斷空間杜賓面板數(shù)據(jù)模型是否可簡化為SLPDM或者SEPDM。根據(jù)以上LM檢驗的結(jié)果，首先選擇建立時空雙固定效應(yīng)的SDPDM進行檢驗。SDPDM圍繞零假設(shè)H0：θ=0以及H0：θ+δβ=0，進一步展開討論。如果以上兩個零假設(shè)均被拒絕，那么說明SDPDM將能夠最好地描述碳排放強度的數(shù)據(jù)變動。根據(jù)Wald和LR檢驗結(jié)果可以發(fā)現(xiàn)，兩個零假設(shè)均被拒絕，由此判斷SLPDM或者SEPDM兩個模型均被拒絕，因此，選擇SDPDM將更加準確地刻畫各省區(qū)資源要素對碳排放空間效應(yīng)影響。進一步采用Hausman檢驗方法判斷，結(jié)果顯示應(yīng)該選擇反映空間隨機效應(yīng)以及時間固定效應(yīng)的SDPDM模型。具體檢驗結(jié)果如表2所示，我國各省區(qū)資源要素對碳排放影響存在著較強的時空外溢效應(yīng)，表明傳統(tǒng)的面板數(shù)據(jù)模型沒有考慮到空間效應(yīng)問題，其結(jié)果將會存在一定的偏誤。

在以碳排放強度的自然對數(shù)為被解釋變量的杜賓模型揭示出各因素對碳排放強度的作用效果。研究結(jié)果表明，人口規(guī)模和城鎮(zhèn)化水平的估計系數(shù)為負，通過5%的顯著性水平檢驗，說明人口總量的增加以及城鎮(zhèn)化進程的進一步加快，對碳排放強度變動產(chǎn)生抑制作用。能源強度和產(chǎn)業(yè)結(jié)構(gòu)估計系數(shù)為正，通過1%的顯著性水平檢驗，說明在重工業(yè)比重增加，產(chǎn)業(yè)結(jié)構(gòu)不合理的情況下，如果地區(qū)能源效率水平不能得到有效提升，將驅(qū)動碳排放強度的增加。人均收入水平估計系數(shù)為正，通過5%的顯著性水平檢驗，說明人均收入水平是碳排放強度變動的主要因素。能源消費結(jié)構(gòu)估計系數(shù)為正，但是并不顯著，說明能源消費結(jié)構(gòu)對碳排放強度變動產(chǎn)生促進作用。能源價格與對外開放水平估計系數(shù)為負，但是并不顯著，說明能源價格的調(diào)整和對外開放水平擴大有助于推動碳排放強度的下降，雖然影響作用有限，但是意義重大。

進一步分析以碳排放強度的自然對數(shù)為被解釋變量和其他因素為解釋變量的空間滯后項，結(jié)果顯示各影響因素產(chǎn)生的空間溢出性存在著較大差異?？臻g滯后人口規(guī)模、空間滯后城鎮(zhèn)化水平與空間滯后人均收入水平自然對數(shù)估計系數(shù)顯著為正，證明存在人口規(guī)模與城鎮(zhèn)化水平的上升將引致相鄰區(qū)域碳排放強度的上升，產(chǎn)生碳排放空間外溢效應(yīng)?？臻g滯后能源強度和空間滯后產(chǎn)業(yè)結(jié)構(gòu)估計系數(shù)為負，但均不顯著，說明技術(shù)進步與產(chǎn)業(yè)結(jié)構(gòu)轉(zhuǎn)型升級將推動鄰近區(qū)域碳排放強度下降，產(chǎn)生的碳排放空間外溢效應(yīng)，但是效果并不顯著。能源價格估計系數(shù)為正，但均不顯著，說明能源價格的調(diào)整將促使鄰近區(qū)域碳排放強度上升，但是產(chǎn)生的碳排放強度空間外溢效應(yīng)并不顯著?？臻g滯后能源消費結(jié)構(gòu)和空間滯后對外開放水平估計系數(shù)為負，通過5%的顯著性水平檢驗，說明對外開放水平的擴大與能源消費結(jié)構(gòu)的調(diào)整有助于驅(qū)動鄰近區(qū)域碳排放強度的下降，呈現(xiàn)比較顯著的空間外溢效應(yīng)。我國在經(jīng)濟全球化以及區(qū)域經(jīng)濟一體化的背景下，對沿海、沿江地區(qū)加大對外開放的力度，同時加快中西部地區(qū)對外開放的步伐，為中國經(jīng)濟發(fā)展向低碳轉(zhuǎn)型奠定了堅實的基礎(chǔ)。

4? 結(jié)論及政策建議

中國特色社會主義進入新時代，新發(fā)展理念統(tǒng)領(lǐng)下的區(qū)域協(xié)調(diào)發(fā)展更注重質(zhì)量和協(xié)同性。對我國碳排放強度空間面板數(shù)據(jù)進行測算，提出探索性時空數(shù)據(jù)分析框架，運用GeoDa軟件刻畫我國各省區(qū)碳排放強度時空關(guān)聯(lián)及其局部空間差異的波動性，將STIRPAT模型和EKC模型相結(jié)合檢驗我國碳排放強度影響因素的溢出性，揭示資源要素對碳排放空間效應(yīng)的傳導機制，進一步提出促進綠色低碳轉(zhuǎn)型的創(chuàng)新政策。研究結(jié)果表明：①梯度發(fā)展模式下碳排放強度的空間集聚促使我國碳排放強度水平相近的省區(qū)將會進一步收斂。碳排放強度水平相近的省區(qū)均呈現(xiàn)出先減弱后增強的空間集聚現(xiàn)象，而且相比于20年前，這種空間集聚現(xiàn)象有所增強，我國省際碳排放強度存在俱樂部收斂的演進特征。②我國省際碳排放強度存在空間溢出性特征，高碳排放強度集聚區(qū)的輻射作用要遠高于低碳排放強度集聚區(qū)，碳排放強度較高的地區(qū)是節(jié)能減排工作的重點。我國碳排放強度的局部空間自相關(guān)呈現(xiàn)集聚與分異并存的時空分布特征，6個時期碳排放強度表現(xiàn)為聚集特征的省區(qū)已經(jīng)超過63.33%，說明我國碳排放強度具有一定的空間鎖定或路徑依賴特征。③我國各省區(qū)資源要素協(xié)同作用是驅(qū)動碳排放強度下降的重要原因。資源要素對碳排放空間效應(yīng)的傳導機制不僅能夠很好的反映相鄰地區(qū)碳排放強度的時空關(guān)聯(lián)性，而且有利于解釋碳排放強度的非均衡空間分布問題。我國碳排放強度空間集聚、分異與躍遷的形成與演變是各地區(qū)資源稟賦以及經(jīng)濟發(fā)展水平等多種資源要素協(xié)同作用的結(jié)果。④協(xié)同發(fā)展的機制創(chuàng)新是實現(xiàn)綠色低碳轉(zhuǎn)型，構(gòu)建推動我國經(jīng)濟高質(zhì)量發(fā)展體制機制的著力點。為實現(xiàn)我國政府向國際承諾的碳排放強度下降以及碳排放總量達峰的低碳經(jīng)濟轉(zhuǎn)型發(fā)展目標，需要將空間關(guān)聯(lián)效應(yīng)納入減排政策制定和政策效果評價當中，加強對碳排放強度“穩(wěn)定”以及“躍遷”的關(guān)鍵省區(qū)進行有效監(jiān)測與治理，進一步采取分區(qū)域、有重點的區(qū)域協(xié)同減排措施，提升區(qū)域碳減排協(xié)同治理的效率。

碳排放強度的準確評估是建立碳減排政策的基礎(chǔ)，因此，國家和地區(qū)碳減排政策的制定需要重點關(guān)注各省區(qū)資源要素協(xié)同，降低碳排放強度是實現(xiàn)高質(zhì)量經(jīng)濟發(fā)展的重要組成部分。為了實現(xiàn)我國在巴黎協(xié)定國家自主貢獻上承諾的國際承諾減排目標，中央以及地方政府需要明晰資源要素對碳排放空間效應(yīng)的作用機理和空間傳導機制，充分發(fā)揮空間溢出效應(yīng)，因地制宜制定差異化碳減排政策。第一，采取分區(qū)域因地制宜制定碳排放強度控制措施，重點控制碳排放強度水平高，增速又較快的空間集聚省區(qū);適當控制碳排放強度水平趨緩，增速較為穩(wěn)定的空間集聚省區(qū)。從區(qū)域碳排放的溢出效應(yīng)來看，對于區(qū)域間碳排放轉(zhuǎn)移，不能產(chǎn)生正向的經(jīng)濟增長效應(yīng)的產(chǎn)業(yè)轉(zhuǎn)移，應(yīng)嚴格執(zhí)行國家相關(guān)政策法規(guī)，淘汰落后產(chǎn)能。第二，碳排放強度空間集聚顯著的省區(qū)亟需加快產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化與升級，改變資源密集與高碳產(chǎn)品出口貿(mào)易格局，突破高碳排放強度的“空間鎖定”與“路徑依賴”。促進傳統(tǒng)重化工業(yè)向低消耗、輕污染、高技術(shù)產(chǎn)業(yè)拓展和結(jié)構(gòu)優(yōu)化升級，延長產(chǎn)業(yè)鏈，提高出口產(chǎn)品附加值，改變傳統(tǒng)以出口高碳產(chǎn)品為主的貿(mào)易格局，減少貿(mào)易中的碳泄露，有效抑制該地區(qū)碳排放強度對碳排放份額的收斂效應(yīng)。利用碳排放控制監(jiān)測形成的倒逼機制，促進該區(qū)域的產(chǎn)業(yè)技術(shù)水平升級與產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化。第三，碳排放強度空間集聚顯著的省區(qū)應(yīng)做好跨區(qū)域功能布局與協(xié)同治理。根據(jù)不同的發(fā)展階段，對溢出熱點集聚區(qū)的碳排放強度減排目標進行嚴格控制，通過建立更加高效的區(qū)域協(xié)調(diào)發(fā)展新機制，控制集聚區(qū)的碳排放，實現(xiàn)集聚區(qū)域內(nèi)部不同功能區(qū)域的分工合作和資源最優(yōu)化組合。依托高效公共交通體系促進集聚區(qū)域一體化發(fā)展，進一步強化地區(qū)碳排放強度對碳排放份額的發(fā)散效應(yīng)。在供給側(cè)改革中，地方政府要讓要素自由流動從而促進資源優(yōu)化配置，同時要從區(qū)域規(guī)劃層面促進各地區(qū)之間在綠色低碳轉(zhuǎn)型上的協(xié)同創(chuàng)新，實現(xiàn)區(qū)域經(jīng)濟的高質(zhì)量發(fā)展。

（編輯：李 琪）

參考文獻

[1] 何建坤， 劉濱. 作為溫室氣體排放衡量指標的碳排放強度分析[J]. 清華大學學報（自然科學版）， 2004， 44（6）：740-743.

[2] 潘家華， 黃承梁， 李萌. 系統(tǒng)把握新時代生態(tài)文明建設(shè)基本方略——對黨的十九大報告關(guān)于生態(tài)文明建設(shè)精神的解讀[J]. 環(huán)境經(jīng)濟， 2017（20）：66-71.

[3]WANG C， WANG F. China can lead on climate change [J]. Science， 2017， 357（6353）：764.

[4] YE F， XIE X， ZHANG L， et al. An improved grey model and scenario analysis for carbon intensity forecasting in the pearl river delta region of China[J]. Energies， 2018， 11（1）：91-103.

[5]陳詩一. 中國碳排放強度的波動下降模式及經(jīng)濟解釋[J]. 世界經(jīng)濟， 2011（4）：124-143.

[6] 林伯強， 毛東昕. 中國碳排放強度下降的階段性特征研究[J]. 金融研究， 2014（8）：101-117.

[7] DONG L， LIANG H. Spatial analysis on Chinas regional air pollutants and CO2， emissions： emission pattern and regional disparity [J]. Atmospheric environment， 2014， 92：280-291.

[8]DOBES L， JOTZO F， STERN D I. The economics of global climate change： a historical literature review[J]. Review of economics， 2014， 65（3）：281-320.

[9] 張翠菊， 張宗益. 中國省域碳排放強度的集聚效應(yīng)和輻射效應(yīng)研究[J]. 環(huán)境科學學報， 2017， 37（3）：1178-1184.

[10]WANG S， CHU C， CHEN G， et al. Efficiency and reduction cost of carbon emissions in China： a nonradial directional distance function method[J]. Journal of cleaner production， 2016， 113：624-634.

[11] 張友國. 經(jīng)濟發(fā)展方式變化對中國碳排放強度的影響[J]. 經(jīng)濟研究， 2010（4）：120-133.

[12] 趙桂梅， 陳麗珍， 孫華平. 基于異質(zhì)性收斂的中國碳排放強度脫鉤效應(yīng)研究[J]. 華東經(jīng)濟管理， 2017， 31（4）：97-103.

[13] AGUIAR A P D， OMETTO J P， NOBRE C， et al. Modeling the spatial and temporal heterogeneity of deforestation：driven carbon emissions： the INPEEM framework applied to the Brazilian Amazon[J]. Global change biology， 2012， 18（11）：3346-3366.

[14] 馮宗憲， 王凱瑩. 中國省域碳強度集群的空間統(tǒng)計分析[J]. 資源科學， 2014， 36（7）：1462-1468.

[15] CHEN Y H H， PALTSEV S， REILLY J M， et al. Longterm economic modeling for climate change assessment[J]. Economic modelling， 2016， 52：867-883.

[16] GENG Y， SARKIS J ， BLEISCHWITZ R . How to globalize the circular economy[J]. Nature， 2019， 565（7738）：153-155.

[17] SUN H，EDZIAH B K，SUN C，et al.Institutional quality， green innovation and energy efficiency [J].Energy policy，2019.DOI： 10.1016/j.enpol.2019.111002.

[18] 林毅夫. 改革開放40年 中國經(jīng)濟如何創(chuàng)造奇跡[J]. 金融經(jīng)濟， 2018（1）：17-18.

[19] SERGIO J R. Spatial empirics for economic growth and convergence[J]. Geographical analysis，2001， 33（3）：195-214.

[20]李子奈， 李鯤鵬. 關(guān)于計量經(jīng)濟學模型隨機擾動項的討論[J]. 統(tǒng)計研究， 2009， 26（2）：62-67.

[21] 張友國， 鄭玉歆. 碳強度約束的宏觀效應(yīng)和結(jié)構(gòu)效應(yīng)[J]. 中國工業(yè)經(jīng)濟， 2014（6）：57-69.

[22] HAO Y， LIAO H， WEI Y M. Is Chinas carbon reduction target allocation reasonable？ an analysis based on carbon intensity convergence[J]. Applied energy， 2015， 142：229-239.