臧　正，鄒欣慶,2,*

1 南京大學(xué)地理與海洋科學(xué)學(xué)院, 南京　210023 2 中國南海研究協(xié)同創(chuàng)新中心， 南京　210093

?

中國大陸省際生態(tài)-經(jīng)濟(jì)效率的時(shí)空格局及其驅(qū)動(dòng)因素

臧正1，鄒欣慶1,2,*

1 南京大學(xué)地理與海洋科學(xué)學(xué)院, 南京210023 2 中國南海研究協(xié)同創(chuàng)新中心， 南京210093

摘要:在參考和借鑒有關(guān)生態(tài)系統(tǒng)服務(wù)價(jià)值理論研究成果的基礎(chǔ)上，提出生態(tài)-經(jīng)濟(jì)產(chǎn)出效率概念及其內(nèi)涵、表征方法，結(jié)合有關(guān)統(tǒng)計(jì)數(shù)據(jù)計(jì)算了2001—2013年中國大陸31省區(qū)生態(tài)-經(jīng)濟(jì)產(chǎn)出效率并就其時(shí)空分布及重心移動(dòng)方向、移動(dòng)距離等進(jìn)行分析，在此基礎(chǔ)上應(yīng)用探索性空間數(shù)據(jù)分析方法及偏最小二乘回歸模型對其空間集聚特征和驅(qū)動(dòng)因素進(jìn)行實(shí)證分析。研究結(jié)果表明：考察時(shí)段內(nèi)中國大陸省際生態(tài)-經(jīng)濟(jì)產(chǎn)出效率自西向東依次提高的三級梯度分布格局趨于形成；省區(qū)之間生態(tài)-經(jīng)濟(jì)產(chǎn)出效率具有全局空間依賴性和空間正相關(guān)特征，顯著高-高集聚和顯著低-高集聚表現(xiàn)相對穩(wěn)定、高-低集聚區(qū)的生態(tài)效率溢出效應(yīng)促使顯著低-低集聚區(qū)逐步消失；人均國民生產(chǎn)總值、單位面積內(nèi)的道路交通里程數(shù)、科研資金投入占GDP的比重等三因素顯著驅(qū)動(dòng)了中國大陸省際生態(tài)-經(jīng)濟(jì)產(chǎn)出效率溢出，但東、中、西部地區(qū)的制約因素不盡相同。

關(guān)鍵詞：生態(tài)經(jīng)濟(jì)學(xué)；生態(tài)效率；實(shí)證分析；時(shí)空格局；驅(qū)動(dòng)力

生態(tài)效率是指人類社會(huì)發(fā)展過程中獲得的經(jīng)濟(jì)價(jià)值增量和由此導(dǎo)致的環(huán)境影響增量的比值[1]，因其能夠?qū)⒖沙掷m(xù)發(fā)展的宏觀目標(biāo)與中、微觀尺度的區(qū)域發(fā)展規(guī)劃或企業(yè)管理有效結(jié)合在一起，使得相關(guān)研究成為生態(tài)學(xué)、經(jīng)濟(jì)學(xué)等學(xué)科領(lǐng)域日益關(guān)注的焦點(diǎn)問題[2- 3]。隨著世界可持續(xù)發(fā)展工商業(yè)聯(lián)合會(huì)(WBCSD)、經(jīng)濟(jì)發(fā)展與合作組織(OECD)、歐洲環(huán)境署(EEA)等國際組織的推廣[4- 6]，各國學(xué)者在不同時(shí)空尺度內(nèi)開展了一系列生態(tài)效率理論與實(shí)踐應(yīng)用研究：在理論方法上，由于研究目的、研究對象或研究視角不同，導(dǎo)致對經(jīng)濟(jì)價(jià)值及環(huán)境影響的度量存在差異，基于投入-產(chǎn)出分析的單一比值法、多指標(biāo)綜合評價(jià)法、數(shù)據(jù)包絡(luò)分析模型等在不同空間尺度或不同行業(yè)得到了實(shí)證[7- 12]，一些研究還針對傳統(tǒng)方法存在的弊端或缺陷進(jìn)行了相應(yīng)改進(jìn)；在實(shí)踐應(yīng)用中，由于廣大中小企業(yè)的資金或技術(shù)限制而難以實(shí)施持續(xù)的生態(tài)效率策略[13]，因此國外大型企業(yè)得以率先開展生態(tài)效率評價(jià)并將其與企業(yè)發(fā)展戰(zhàn)略、決策聯(lián)系起來[14- 15]，研究視角也不斷拓寬并逐漸與溫室效應(yīng)、糧食安全、生物多樣性保護(hù)等有關(guān)全球變化問題聯(lián)系到一起，促使相關(guān)研究成為全球可持續(xù)發(fā)展研究領(lǐng)域的熱點(diǎn)。

生態(tài)效率概念及其評價(jià)方法引入我國之后，在有關(guān)其理論方法、影響因素識(shí)別及區(qū)域?qū)嵶C等方面取得了一系列成果：在理論方法研究方面，劉晶茹等從復(fù)合生態(tài)系統(tǒng)視角提出了復(fù)合生態(tài)效率概念，據(jù)此構(gòu)建了生態(tài)產(chǎn)業(yè)園的生態(tài)效率評價(jià)指標(biāo)體系[16]；季丹、芮俊偉等分別應(yīng)用生態(tài)足跡法針對中國大陸省區(qū)以及小尺度區(qū)域生態(tài)工業(yè)園的生態(tài)效率開展了實(shí)證研究[17- 18]；潘丹等、李健等分別基于非期望產(chǎn)出的SBM模型與非參數(shù)距離函數(shù)模型評價(jià)了中國大陸農(nóng)業(yè)生態(tài)效率和區(qū)域生態(tài)效率[19- 20]。在區(qū)域生態(tài)效率評價(jià)研究方面：張淑英等通過環(huán)境壓力指標(biāo)、應(yīng)用Tobit模型分析了我國區(qū)域工業(yè)生態(tài)效率的影響因素，發(fā)現(xiàn)人均國民生產(chǎn)總值和地理位置具有驅(qū)動(dòng)作用(正面因素)、而外向型經(jīng)濟(jì)結(jié)構(gòu)卻在一定程度上制約(負(fù)面因素)區(qū)域生態(tài)效率提高[21]；陳作成等從系統(tǒng)協(xié)同的視角分析了我國西部地區(qū)產(chǎn)業(yè)系統(tǒng)生態(tài)效率[22]，研究結(jié)果表明資本投入及技術(shù)進(jìn)步能夠促進(jìn)清潔生產(chǎn)效率提高、市場化水平提高起負(fù)面作用，而外商直接投資對生態(tài)環(huán)境的作用則表現(xiàn)為“節(jié)能”而非“減排”；陳傲利用因子分析模型探討了中國區(qū)域生態(tài)效率的差異，結(jié)果表明環(huán)保資金投入及產(chǎn)業(yè)結(jié)構(gòu)調(diào)整對改善區(qū)域生態(tài)效率有積極影響、而以排污費(fèi)為代表的環(huán)境經(jīng)濟(jì)政策的顯著作用尚未發(fā)現(xiàn)[23]；在生態(tài)效率的影響因素研究方面，汪克亮等、崔瑋等分別基于節(jié)能視角和減排視角展開了實(shí)證研究，發(fā)現(xiàn)資源配置不合理、環(huán)境管理水平低下等兩方面因素是制約中國大陸地區(qū)生態(tài)效率提升的主導(dǎo)因素[24]，認(rèn)為通過提高區(qū)域技術(shù)進(jìn)步與技術(shù)效率指數(shù)有助于縮小地區(qū)生態(tài)效率的顯著差異[25]；初善冰等、羅能生等、潘興俠等分別探討了外商直接投資、城鎮(zhèn)化水平及產(chǎn)業(yè)結(jié)構(gòu)對區(qū)域生態(tài)效率差異的影響，發(fā)現(xiàn)外商直接投資對東部地區(qū)提高生態(tài)效率的正向影響較中、西部地區(qū)顯著[26]，城鎮(zhèn)化水平的提高對東、西部地區(qū)生態(tài)效率的影響方向截然相反[27]，而產(chǎn)業(yè)結(jié)構(gòu)不合理是導(dǎo)致我國中部地區(qū)生態(tài)效率低的首要因素[28]；潘興俠從區(qū)域經(jīng)濟(jì)-資源-環(huán)境復(fù)合系統(tǒng)的視角開展了生態(tài)效率綜合評價(jià)研究工作，深入分析了我國區(qū)域生態(tài)效率差異的影響因素及其收斂特征，進(jìn)而從經(jīng)濟(jì)發(fā)展、產(chǎn)業(yè)結(jié)構(gòu)、人力資本、外資利用、環(huán)境政策以及科技實(shí)力等6個(gè)方面提出了提升策略[29]。

縱觀上述文獻(xiàn)及其他代表性研究成果可以看出：以往我國學(xué)者大多沿用了國外提出的生態(tài)效率概念，在核心思想基本一致的基礎(chǔ)上開展了不同時(shí)空尺度的實(shí)證研究[30- 32]，并在中國大陸東、中、西部地區(qū)生態(tài)效率評價(jià)中取得比較一致的認(rèn)識(shí)——概而言之，多數(shù)研究成果普遍支持中國生態(tài)效率具有“東高西低，南優(yōu)北劣”的觀點(diǎn)，進(jìn)而為不同地區(qū)提高生態(tài)效率、突破生態(tài)效率的梯度限制等提供了科學(xué)建議與對策；但是，由于研究視角不同[22, 24- 25]、所用評價(jià)方法不同[17,20]、選取的指標(biāo)不同[16,21]，導(dǎo)致有關(guān)研究結(jié)果存在不確定性、成果之間缺乏可比性等問題[26- 28]；而從國內(nèi)、外對比來看，國外研究已經(jīng)日益重視生態(tài)效率與全球變暖、糧食安全、生物多樣性等重大問題的密切聯(lián)系，為此從學(xué)科交叉視角對傳統(tǒng)生態(tài)效率評估模型及方法進(jìn)行改良和修正，并逐步由簡單的生態(tài)效率評價(jià)轉(zhuǎn)向其驅(qū)動(dòng)機(jī)制研究[3]，而國內(nèi)相關(guān)研究在這方面尚有待深入。鑒于此，本文擬在學(xué)習(xí)和參考國內(nèi)外有關(guān)生態(tài)效率研究成果的基礎(chǔ)上，嘗試提出并界定基于生態(tài)系統(tǒng)服務(wù)理論的生態(tài)-經(jīng)濟(jì)產(chǎn)出效率概念及其內(nèi)涵、表征方法，進(jìn)而應(yīng)用探索性空間數(shù)據(jù)分析方法和偏最小二乘回歸法對有關(guān)年度中國大陸省區(qū)生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的時(shí)空關(guān)聯(lián)格局及其驅(qū)動(dòng)因素展開實(shí)證研究，旨在為豐富和發(fā)展有關(guān)生態(tài)效率評價(jià)方法提供一些借鑒、為有關(guān)地區(qū)和部門制定區(qū)域生態(tài)效率提升策略等提供一些參考。

1研究方法

1.1生態(tài)-經(jīng)濟(jì)效率的綜合表征方法

生態(tài)效率概念的主旨思想可概括為“環(huán)境影響最小化-經(jīng)濟(jì)價(jià)值最大化”，即以最小的環(huán)境代價(jià)獲得最大的經(jīng)濟(jì)價(jià)值[2- 3]。評價(jià)生態(tài)效率的核心和難點(diǎn)是經(jīng)濟(jì)與環(huán)境維度的表征，為此本文在借鑒有關(guān)研究成果的基礎(chǔ)上提出基于生態(tài)系統(tǒng)服務(wù)價(jià)值理論的生態(tài)-經(jīng)濟(jì)產(chǎn)出效率評價(jià)模型，以期從生態(tài)系統(tǒng)服務(wù)與經(jīng)濟(jì)發(fā)展成就的聯(lián)系視角豐富有關(guān)生態(tài)效率研究方法。

1.1.1基于生態(tài)系統(tǒng)服務(wù)價(jià)值理論的生態(tài)-經(jīng)濟(jì)效率評價(jià)

生態(tài)系統(tǒng)服務(wù)是生態(tài)系統(tǒng)通過生態(tài)過程向人類提供的、具有一定數(shù)量和質(zhì)量的生態(tài)產(chǎn)品和服務(wù)，通?？梢杂闷湄泿呕问健⒓瓷鷳B(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行衡量。人類生產(chǎn)和生活過程中每年消耗的生態(tài)系統(tǒng)服務(wù)價(jià)值大部分轉(zhuǎn)移、固化到當(dāng)年的國民經(jīng)濟(jì)與社會(huì)發(fā)展成就之中，因此生態(tài)系統(tǒng)服務(wù)的持續(xù)供給是區(qū)域可持續(xù)發(fā)展的基礎(chǔ)，人類通過生態(tài)系統(tǒng)服務(wù)消費(fèi)結(jié)合勞動(dòng)力、資本投入實(shí)現(xiàn)了從生態(tài)系統(tǒng)服務(wù)到經(jīng)濟(jì)社會(huì)收益的轉(zhuǎn)化。為進(jìn)一步衡量作為自然資本投入到經(jīng)濟(jì)生產(chǎn)過程中的生態(tài)系統(tǒng)服務(wù)與相應(yīng)經(jīng)濟(jì)產(chǎn)出的定量關(guān)系，本文借鑒Johnston、M?ler等提出的服務(wù)-收益函數(shù)及生產(chǎn)函數(shù)概念及定義[33- 34]，提出基于生態(tài)系統(tǒng)服務(wù)價(jià)值理論的區(qū)域生態(tài)-經(jīng)濟(jì)產(chǎn)出效率評價(jià)模型如下：

EEE=GDP/ESV

(1)

式中,GDP為國民生產(chǎn)總值(萬元)，代表區(qū)域國民經(jīng)濟(jì)與社會(huì)發(fā)展成就；ESV代表投入到國民經(jīng)濟(jì)與社會(huì)發(fā)展過程中的全部生態(tài)系統(tǒng)服務(wù)；EEE稱為生態(tài)-經(jīng)濟(jì)產(chǎn)出效率，表征區(qū)域國民經(jīng)濟(jì)與社會(huì)發(fā)展過程中生態(tài)系統(tǒng)服務(wù)到經(jīng)濟(jì)與社會(huì)收益的轉(zhuǎn)化能力(無量綱；本文中也將其簡稱為生態(tài)-經(jīng)濟(jì)效率或生態(tài)效率)，其值越大表明消費(fèi)等量生態(tài)系統(tǒng)服務(wù)產(chǎn)出的經(jīng)濟(jì)與社會(huì)效益越大。

生態(tài)系統(tǒng)服務(wù)價(jià)值既包括物質(zhì)產(chǎn)品、文化娛樂服務(wù)等直接利用價(jià)值，也包括維持生態(tài)平衡、保護(hù)環(huán)境等非直接利用價(jià)值，服務(wù)范圍幾乎涵蓋國民生產(chǎn)各個(gè)部門。本文參考謝高地等提出的中國陸地生態(tài)系統(tǒng)服務(wù)價(jià)值評價(jià)方法[35]，應(yīng)用如下模型核算生態(tài)系統(tǒng)服務(wù)價(jià)值：

(2)

式中,P代表原糧市場糧食平均價(jià)格(萬元/t)，Y代表耕地的平均糧食產(chǎn)量(t/hm2)，兩者可在一定程度上代表特定地區(qū)的物價(jià)水平及農(nóng)林牧漁等第一產(chǎn)業(yè)生產(chǎn)力水平；Eij代表基于耕地的食物生產(chǎn)價(jià)值得出的不同種類生態(tài)系統(tǒng)服務(wù)價(jià)值的當(dāng)量因子(無量綱，i、j分別代表第i類生態(tài)系統(tǒng)中的第j種生態(tài)系統(tǒng)服務(wù)價(jià)值的當(dāng)量因子)，表征各類生態(tài)系統(tǒng)為人類提供產(chǎn)品、服務(wù)的能力；Ai代表第i類生態(tài)系統(tǒng)對應(yīng)占有的土地面積(hm2)，可表征區(qū)域資源稟賦(自然條件及生態(tài)系統(tǒng)狀況等)；ESV代表區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值，由食物生產(chǎn)等9種生態(tài)系統(tǒng)服務(wù)價(jià)值表征(其中非直接利用價(jià)值作為間接投入，替代國民經(jīng)濟(jì)與社會(huì)發(fā)展過程中的生態(tài)環(huán)境治理成本)、并可由6類生態(tài)系統(tǒng)供應(yīng)(萬元，年值)。

1.1.2相對生態(tài)效率的表征方法

以往研究成果表明，特定時(shí)期、特定區(qū)域的生態(tài)效率可能同時(shí)受到經(jīng)濟(jì)發(fā)展水平、產(chǎn)業(yè)規(guī)模與結(jié)構(gòu)、生產(chǎn)力發(fā)展水平、科技進(jìn)步及區(qū)域自然條件等的綜合影響[24- 29]，而由前文提出的式(1)可知本文界定的生態(tài)-經(jīng)濟(jì)產(chǎn)出效率EEE是具有絕對內(nèi)涵的生態(tài)效率表征，為便于比較，進(jìn)一步提出以其標(biāo)準(zhǔn)化值表征不同評價(jià)單元的生態(tài)-經(jīng)濟(jì)產(chǎn)出效率：假設(shè)共有n個(gè)待評價(jià)單元，其中第k個(gè)待評價(jià)單元基于其區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值投入(ESVk)獲得的總經(jīng)濟(jì)產(chǎn)出為GDPk，則其相對生態(tài)-經(jīng)濟(jì)產(chǎn)出效率REEEk(下文簡稱相對生態(tài)效率)可表示為：

REEEk=EEEk/EEEn

(3)

式中,EEEn代表所有評價(jià)單元中的同期生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的最大值與最小值之差，EEEk代表評價(jià)單元k的生態(tài)-經(jīng)濟(jì)產(chǎn)出效率與最小值之差；以此為參考標(biāo)準(zhǔn)得出的相對生態(tài)-經(jīng)濟(jì)產(chǎn)出效率REEEk(無量綱)表征其在所有考察樣本中的相對位置：其值越大，表明其生態(tài)效率相對越優(yōu)，REEEk>1表示其優(yōu)于參考標(biāo)準(zhǔn)；反之若REEEk<1，表明其劣于參考標(biāo)準(zhǔn)。

1.2生態(tài)效率的空間關(guān)聯(lián)格局及其驅(qū)動(dòng)力分析方法1.2.1探索性空間數(shù)據(jù)分析方法

空間計(jì)量經(jīng)濟(jì)學(xué)理論認(rèn)為，特定空間單元上的某種經(jīng)濟(jì)地理現(xiàn)象或其屬性值與鄰近單元上的同一現(xiàn)象或?qū)傩灾凳窍嗷リP(guān)聯(lián)的，彼此具有空間依賴性[36- 38]。生態(tài)系統(tǒng)的地域差異及生態(tài)系統(tǒng)服務(wù)的動(dòng)態(tài)性使其在為人類利用過程中具有空間自相關(guān)性和溢出效應(yīng)：某一地區(qū)的生態(tài)效率提高或降低，將對本地區(qū)及其他地區(qū)產(chǎn)生特定影響。本文據(jù)此引入探索性空間數(shù)據(jù)分析方法對不同空間單元、不同時(shí)期的生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的時(shí)空關(guān)聯(lián)格局進(jìn)行初步分析，以便為進(jìn)一步揭示驅(qū)動(dòng)或制約中國大陸地區(qū)生態(tài)效率提升的影響因素奠定基礎(chǔ)。

探索性空間數(shù)據(jù)分析方法(ESDA)是一系列空間數(shù)據(jù)分析技術(shù)與方法的集合，能夠借助可視化的表達(dá)方法描述數(shù)據(jù)的空間分布特征、識(shí)別空間數(shù)據(jù)的結(jié)構(gòu)及其異常值，并可據(jù)此檢測某些地理現(xiàn)象的空間集聚效應(yīng)、揭示相互之間的空間作用機(jī)制等，因此常用于空間統(tǒng)計(jì)研究[39]。空間權(quán)重矩陣是空間經(jīng)濟(jì)計(jì)量模型的應(yīng)用前提，一般可通過空間鄰接規(guī)則(公共邊界或公共點(diǎn))、多邊形原則、距離最小原則構(gòu)建；考慮到省際大尺度空間限制，本文擬在對多種空間賦權(quán)方案進(jìn)行比較的基礎(chǔ)上，嘗試依據(jù)距離衰減原理、通過各省區(qū)重心之間的距離確定空間權(quán)重矩陣，以期進(jìn)一步提高評價(jià)結(jié)果的可靠性。分析過程應(yīng)用MapInfo10.0(繪制空間分布圖)和GeoDa095i(計(jì)算空間自相關(guān)指數(shù)，繪制空間集聚圖、散點(diǎn)圖)完成，因文章篇幅所限、相關(guān)理論基礎(chǔ)及具體應(yīng)用步驟恕不贅述，請參閱相關(guān)文獻(xiàn)[39- 40]。

1.2.2偏最小二乘回歸分析方法

區(qū)域自然地理?xiàng)l件無疑會(huì)導(dǎo)致生態(tài)系統(tǒng)服務(wù)供給具有空間異質(zhì)性，而經(jīng)濟(jì)與社會(huì)發(fā)展的本底條件、甚至人類個(gè)體或群體的心理與文化因素等也在不同程度上影響區(qū)域生態(tài)系統(tǒng)服務(wù)消費(fèi)總量及消耗強(qiáng)度，因此區(qū)域生態(tài)效率可能同時(shí)受到多重自然-人文因素影響：自然因素導(dǎo)致生態(tài)系統(tǒng)服務(wù)與經(jīng)濟(jì)產(chǎn)出關(guān)系具有一定的非線性特征，而人文因素如文化背景、科技進(jìn)步、市場機(jī)制因素等可能導(dǎo)致生態(tài)效率與區(qū)域經(jīng)濟(jì)和社會(huì)發(fā)展存在更加復(fù)雜的耦合關(guān)系。因此，在初步分析區(qū)域生態(tài)效率的時(shí)空關(guān)聯(lián)及其溢出或收斂特征的基礎(chǔ)上，進(jìn)一步探討不同時(shí)空尺度生態(tài)效率的驅(qū)動(dòng)或制約因素就顯得極具理論與實(shí)踐價(jià)值。

偏最小二乘回歸法(PLS)是基于傳統(tǒng)最小二乘法(OLS)改進(jìn)的新型多元統(tǒng)計(jì)分析方法，可對系統(tǒng)信息進(jìn)行有效分解和篩選，因此適用于多自變量對因變量的回歸建模，特別是變量之間存在多重共線性時(shí)其可靠性更高[41]。本文基于區(qū)域生態(tài)效率與經(jīng)濟(jì)社會(huì)發(fā)展水平及產(chǎn)業(yè)結(jié)構(gòu)、產(chǎn)業(yè)規(guī)模等具有高度相關(guān)性的考慮，引入PLS模型對驅(qū)動(dòng)或制約中國大陸不同地區(qū)生態(tài)-經(jīng)濟(jì)產(chǎn)出效率空間關(guān)聯(lián)格局的影響因素進(jìn)行分析：分別以東、中、西部省區(qū)的相對生態(tài)-經(jīng)濟(jì)產(chǎn)出效率REEE為自變量，同時(shí)參考相關(guān)生態(tài)效率研究文獻(xiàn)所采用的部分評價(jià)指標(biāo)的基礎(chǔ)上[24- 29]，最終選擇經(jīng)過標(biāo)準(zhǔn)化處理的人均國民生產(chǎn)總值(X01)、外商直接投資額(X02)、進(jìn)出口總額(X03)、城鎮(zhèn)人口比例(X04)、技術(shù)市場成交額占GDP比重(X05)、單位國土面積內(nèi)的道路交通里程數(shù)(X06)、環(huán)境保護(hù)與治理資金投入額(X07)、大中專學(xué)生人數(shù)占地區(qū)人口比重(X08)、環(huán)境技術(shù)效率全局參比(X09)及科技研發(fā)投入占GDP比重(X10)等10個(gè)指標(biāo)為條件變量，從地區(qū)經(jīng)濟(jì)發(fā)展水平(X01)、外資投入(X02)及產(chǎn)業(yè)規(guī)模(X03)、技術(shù)進(jìn)步(X05、X10)及環(huán)境效率因素(X07、X09)、區(qū)域本底條件(X04、X06)及文化背景(X08)等幾個(gè)方面進(jìn)行回歸分析；分析過程通過SIMCA-P_11.5軟件實(shí)現(xiàn)，具體操作步驟請參閱相關(guān)文獻(xiàn)[41]。

2結(jié)果分析

2.1數(shù)據(jù)來源及有關(guān)說明

考慮到數(shù)據(jù)的可得性以及文章的篇幅所限，本文暫以中國大陸31個(gè)省區(qū)為研究對象就其2001—2013年生態(tài)效率的時(shí)空演變格局及其驅(qū)動(dòng)因素進(jìn)行實(shí)證分析，有關(guān)數(shù)據(jù)來源及說明如表1所示：各省區(qū)歷年人均國民生產(chǎn)總值、外商直接投資額等影響因素指標(biāo)X01—X10及糧食單產(chǎn)依據(jù)國家統(tǒng)計(jì)局網(wǎng)站提供的有關(guān)數(shù)據(jù)計(jì)算、整理得出。原糧市場價(jià)格分別以國家發(fā)展和改革委員會(huì)網(wǎng)站公布的2012年相關(guān)省區(qū)糧食最低收購價(jià)格及原糧市場谷物平均價(jià)為基礎(chǔ)，結(jié)合有關(guān)年度物價(jià)指數(shù)折算得出當(dāng)年價(jià)；土地利用面積以中國統(tǒng)計(jì)年鑒公布的2005—2008年數(shù)據(jù)為基礎(chǔ)，結(jié)合中國國家資源環(huán)境遙感宏觀調(diào)查與動(dòng)態(tài)研究所、中國國土資源統(tǒng)計(jì)年鑒(光盤版)等提供的增減變動(dòng)數(shù)據(jù)計(jì)算整理得出；由于統(tǒng)計(jì)數(shù)據(jù)缺失，研究內(nèi)容暫未包括香港和澳門特別行政區(qū)以及臺(tái)灣省、南海諸島等。

表1　數(shù)據(jù)來源及說明

東部地區(qū)(省區(qū))包括北京、天津、河北、遼寧、山東、江蘇、上海、浙江、福建、廣東、海南等11省區(qū)；中部地區(qū)(省區(qū))包括吉林、黑龍江、山西、河南、湖北、湖南、安徽、江西等8省區(qū)；西部地區(qū)(省區(qū))包括內(nèi)蒙古、陜西、甘肅、寧夏、新疆、青海、西藏、重慶、四川、貴州、云南、廣西等12省區(qū)。內(nèi)蒙古、黑龍江、吉林、遼寧等4省區(qū)為玉米主產(chǎn)區(qū)，P值以玉米最低收購價(jià)格為準(zhǔn)；河北、江蘇、安徽、山東、河南、湖北等6省區(qū)為小麥主產(chǎn)區(qū)，P值以小麥最低收購價(jià)格為準(zhǔn)；江西、湖南、廣西、四川4省(區(qū))為中晚稻主產(chǎn)區(qū)，P值以中晚稻最低收購價(jià)格為準(zhǔn)；其他17省區(qū)P值以原糧市場的玉米、小麥、稻米價(jià)格的均值為準(zhǔn)

2.2省際生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的時(shí)空演變趨勢分析

依據(jù)式(1)—(2)及有關(guān)數(shù)據(jù)計(jì)算2001—2013年中國大陸31個(gè)省區(qū)的生態(tài)-經(jīng)濟(jì)產(chǎn)出效率及其重心位置、重心移動(dòng)方向與距離等，結(jié)果如圖1(按照3個(gè)梯度聚類顯示；鑒于研究時(shí)段內(nèi)空間分布格局的演變趨勢具有較高連續(xù)性，故僅列出了4個(gè)代表年度)、表2所示。

2.2.1省際生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的時(shí)空分布格局分析

從圖1可看出：2001年大陸東部的北京、天津、河北、遼寧、山東、江蘇、上海、浙江、福建、廣東及海南等11省區(qū)(以下稱東部省區(qū))中，除了北京、天津、上海3個(gè)直轄市的生態(tài)-經(jīng)濟(jì)產(chǎn)出效率EEE高于10(梯度Ⅰ)之外，東部其他8省區(qū)及黑龍江、吉林、山西、河南、安徽、湖北、湖南、江西等8個(gè)省區(qū)(以下稱中部省區(qū))的EEE值整體上(除黑龍江之外)處于(1, 10 ]區(qū)間(梯度Ⅱ)；內(nèi)蒙古等其余12省區(qū)(以下稱西部省區(qū))的EEE值則基本(重慶除外)處于(0, 1 ]區(qū)間(梯度Ⅲ)。整體來看，這一時(shí)期中國大陸省區(qū)的生態(tài)-經(jīng)濟(jì)產(chǎn)出效率以“胡煥庸線(黑河-騰沖)”為界、呈現(xiàn)比較顯著的“東部和中部地區(qū)較高(不含東北、西南兩翼)，西部地區(qū)較低”的空間分布格局。

圖1　中國大陸省區(qū)生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的時(shí)空分布格局(2001—2013年)Fig.1　Temporal and spatial distribution patterns of eco-economic output efficiency in Chinese mainland

圖1顯示2005年東部非直轄市省區(qū)中江蘇率先進(jìn)入生態(tài)效率第一集團(tuán)(梯度Ⅰ)，西部的陜西、四川、貴州及廣西4地則升至梯度Ⅱ集團(tuán)。此后(2009年，圖1)，山東、浙江、廣東以及中部的河南等4省區(qū)也提升至梯度Ⅰ；由于黑龍江的加入，這一時(shí)期中部地區(qū)仍有7地EEE值處在(1, 10 ]區(qū)間；而隨著云南、寧夏兩地加入梯度Ⅱ集團(tuán)，處于梯度Ⅲ的西部省區(qū)在這一時(shí)期已減少至5個(gè)。綜合上述結(jié)果來看，這一時(shí)段內(nèi)中國大陸的生態(tài)-經(jīng)濟(jì)產(chǎn)出效率較2001年有所提升，在整體保持“東中部地區(qū)高—西部地區(qū)低”的特征的同時(shí)，東北及西南兩翼的省際生態(tài)效率持續(xù)提升，由西向東漸次提高的三級梯度分布格局初具雛形。

從圖1d可以看出：與2005、2009年相比，2013年中國大陸東、中部省區(qū)的EEE值繼續(xù)提升，河北與遼寧兩地進(jìn)入梯度Ⅰ集團(tuán)、內(nèi)蒙古也向上提升一級，促使梯度Ⅱ的范圍整體越過“胡煥庸線”，自西向東逐級提升的三級梯度格局基本形成。值得關(guān)注的是，這一時(shí)期西部的甘肅、青海、新疆、西藏等4省區(qū)(以下簡稱“甘肅等西部4省區(qū)”)生態(tài)效率仍未得到有效提升，導(dǎo)致其在整個(gè)研究時(shí)段內(nèi)均處于梯度Ⅲ集團(tuán)。

2.2.2省際生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的重心轉(zhuǎn)移路徑分析

從表2可以看出：研究時(shí)段內(nèi)中國大陸省區(qū)生態(tài)-經(jīng)濟(jì)產(chǎn)出效率EEE的重心始終處在東經(jīng)116.42—116.71°、北緯34.35—34.54°(E116°25′—E116°43′，E34°21′—E34°32′)范圍之內(nèi)(圖1，位于安徽省境內(nèi))。這一結(jié)果表明研究時(shí)段內(nèi)中國大陸省區(qū)的生態(tài)-經(jīng)濟(jì)產(chǎn)出效率重心位置比較穩(wěn)定、整體移動(dòng)速度相對緩慢(年均約3km)，與圖1所隱含的3個(gè)梯度格局演變趨勢相吻合——研究時(shí)段內(nèi)中國大陸生態(tài)效率整體提升，結(jié)果趨于均衡。進(jìn)一步對比表2與圖1可以看出：2001—2006年中國大陸生態(tài)效率重心向北或東北方向移動(dòng)、整體移動(dòng)距離僅為12.05km，2007—2013年向西或西北方向移動(dòng)、移動(dòng)距離為24.92km；而圖1顯示研究時(shí)段內(nèi)中國大陸省區(qū)生態(tài)-經(jīng)濟(jì)產(chǎn)出效率格局的演變特征存在局部差異：中部省區(qū)及甘肅等西部4省區(qū)提升速度緩慢，東部省區(qū)及西部其它8省區(qū)提升速度相對較快。上述結(jié)果表明在東部省區(qū)生態(tài)效率提升空間逐漸縮小的條件下，未來進(jìn)一步提高中國大陸生態(tài)效率的潛力主要在中、西部尤其是甘肅等西部4省區(qū)之間，因此進(jìn)一步分析不同省區(qū)生態(tài)效率的局部空間關(guān)聯(lián)、進(jìn)而識(shí)別不同梯度內(nèi)區(qū)域生態(tài)效率提升的驅(qū)動(dòng)或制約因素，就顯得尤為必要。

2.3省際生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的空間集聚效應(yīng)分析2.3.1省際生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的全局空間自相關(guān)分析

通過全局空間自相關(guān)指數(shù)可以識(shí)別省區(qū)之間生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的相似程度，有助于從整體上支撐前述中國大陸生態(tài)效率時(shí)序演變趨勢的判斷。應(yīng)用GeoDa095i軟件得到不同時(shí)期EEE值的全局空間自相關(guān)指數(shù)Moran′ I(莫蘭指數(shù))及其顯著性水平(P值，采用蒙特卡羅模擬檢驗(yàn))如表3所示(為進(jìn)一步提高評價(jià)結(jié)果的穩(wěn)健性，通過4種空間權(quán)重方案計(jì)算)。

表2　中國大陸省區(qū)生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的重心位置及其移動(dòng)方向、移動(dòng)距離(2001—2013年)

表3　中國大陸省區(qū)生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的全局空間自相關(guān)指數(shù)及其顯著性檢驗(yàn)值

方案Ⅰ—Ⅳ(ProgramⅠ—Ⅳ)分別代表依據(jù)公共邊原則(Rook Contiguity)、公共邊或公共點(diǎn)原則(Queen Contiguity)、多邊形原則(K-Nearest Distance)、距離最小原則(Threshold Distance)確定的空間權(quán)重方案

從表3可以看出，依據(jù)4種空間賦權(quán)方案得出的中國大陸省區(qū)生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的全局莫蘭指數(shù)Moran′ I在整個(gè)研究時(shí)段內(nèi)均為正值并且顯著水平均在5%以上，表明大陸省區(qū)的生態(tài)-經(jīng)濟(jì)產(chǎn)出效率在考察期內(nèi)整體具備空間正相關(guān)特征，即生態(tài)效率高的省區(qū)或生態(tài)效率低的省區(qū)各自呈現(xiàn)空間集聚特征(高-高集聚或低-低集聚)。此外，依據(jù)4種賦權(quán)方案得出的Moran′ I均值顯示：研究時(shí)段內(nèi)中國大陸省區(qū)生態(tài)-經(jīng)濟(jì)產(chǎn)出效率空間正相關(guān)水平的高低次序?yàn)榉桨涪?方案Ⅱ>方案Ⅲ>方案Ⅳ；全局莫蘭指數(shù)整體介于0.1—0.5之間并保持相對穩(wěn)定，表明這種高-高集聚或低-低集聚的趨勢年際變化相對不大。據(jù)此，可進(jìn)一步分析個(gè)別省區(qū)或局部地區(qū)生態(tài)-經(jīng)濟(jì)產(chǎn)出效率提高或降低對其他省區(qū)乃至中國大陸整體產(chǎn)生了何種影響，從而為檢驗(yàn)省區(qū)之間生態(tài)-經(jīng)濟(jì)產(chǎn)出效率是否存在溢出效應(yīng)并為驅(qū)動(dòng)因素識(shí)別提供依據(jù)。為了提高評價(jià)結(jié)果的可靠性，下文依據(jù)方案Ⅳ進(jìn)行局部空間自相關(guān)(LISA)分析。

2.3.2省際生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的局部空間自相關(guān)分析

以經(jīng)過標(biāo)準(zhǔn)化處理的各個(gè)省區(qū)有關(guān)年度的生態(tài)-經(jīng)濟(jì)產(chǎn)出效率(本文中，即為REEE值)為橫軸(X2001、X2005、X2009、X2013)、以其空間滯后值(即由權(quán)重方案Ⅳ決定的與該省區(qū)相鄰的所有省區(qū)REEE平均值，W_X2001、W_X2005、W_X2009、W_X2013)為縱軸構(gòu)建莫蘭散點(diǎn)圖，可以更加清晰地描繪每個(gè)省區(qū)與周邊省區(qū)生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的空間集聚或離散趨勢(第一、第二、第三、第四象限分別對應(yīng)高-高集聚、低-高集聚、低-低集聚、高-低集聚)。因全局莫蘭指數(shù)(表3)僅能體現(xiàn)省際生態(tài)效率的整體水平差異及其分布特征的穩(wěn)定性，所以還需進(jìn)一步通過局部Moran指數(shù)分析顯著相似省區(qū)的空間自相關(guān)程度。為此，應(yīng)用GeoDa095i得出各年度中國大陸31個(gè)省區(qū)在2001、2005、2009、2013年四個(gè)時(shí)期生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的散點(diǎn)圖及其高于5%顯著性水平的LISA集聚圖如圖2所示。

圖2　中國大陸省區(qū)生態(tài)效率的散點(diǎn)圖及LISA集聚圖Fig.2　The scattered point map and LISA concentration map of ecological efficiency of the provinces in Chinese mainland

圖2顯示的是在5%和1%兩個(gè)顯著水平上中國大陸省區(qū)生態(tài)效率的LISA集聚圖和散點(diǎn)圖。從圖2右邊的散點(diǎn)圖可以看出在2001、2005、2009、2013年4個(gè)時(shí)期，中國大陸生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的空間負(fù)相關(guān)省區(qū)個(gè)數(shù)(第二、四象限)始終處于11—14個(gè)之間，這些省區(qū)之間生態(tài)-經(jīng)濟(jì)產(chǎn)出效率具有比較顯著的異質(zhì)性特征并且數(shù)量相對穩(wěn)定；空間正相關(guān)的省區(qū)個(gè)數(shù)則介于17—20之間、高-高集聚的“熱點(diǎn)”省區(qū)數(shù)量(第一象限)少于同期低低集聚的“盲點(diǎn)”省區(qū)數(shù)量(第四象限)，并且各時(shí)期分屬于這兩種集聚特征的省區(qū)數(shù)量也比較穩(wěn)定。圖2在2001與2005年兩個(gè)時(shí)期的集聚圖顯示北京、天津、河北、河南、山東、江蘇、上海等7個(gè)省區(qū)呈現(xiàn)高-高集聚特征，2009年河北從該區(qū)域消失；2001年青海及四川兩省區(qū)構(gòu)成顯著的低-低集聚區(qū)、2005年重慶和廣西兩地加入其中、2009年四川退出，這3個(gè)時(shí)期低-低集聚省區(qū)尚比較顯著(5%)，2013年該集團(tuán)則在整個(gè)中國大陸消失；吉林、山西、湖北、安徽等4地構(gòu)成了2001年中國大陸生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的顯著低-高集聚區(qū)，2005年遼寧加入，2009年河北由高-高集聚區(qū)轉(zhuǎn)入，2013年內(nèi)蒙古和寧夏兩地的加入使得該集團(tuán)省區(qū)個(gè)數(shù)增至8個(gè)。此外，依據(jù)權(quán)重方案Ⅳ，新疆因地域面積大而致“鄰居較少(在既定距離閾值上屬于“無鄰居”省區(qū))”，因此未予評價(jià)；由于整個(gè)研究時(shí)段內(nèi)高-低集聚省區(qū)顯著水平低(Not Significant)，在此也不予分析和討論。

2.4省際生態(tài)-經(jīng)濟(jì)產(chǎn)出效率溢出效應(yīng)的驅(qū)動(dòng)因素分析

對比圖2、圖1所得結(jié)果可知：在甘肅等西部4省區(qū)生態(tài)-經(jīng)濟(jì)產(chǎn)出效率低且提升速度緩慢、中部8省區(qū)效率較高但提升空間整體縮小的情況下，東部省區(qū)及西部其他8省區(qū)生態(tài)效率得到了有效提升，尤其是西部其他8省區(qū)生態(tài)效率的持續(xù)提升，促使研究時(shí)段內(nèi)呈顯著低-低集聚特征的省區(qū)從中國大陸整體消失，表明這些省區(qū)生態(tài)效率具有溢出效應(yīng)。

為了進(jìn)一步識(shí)別促使上述溢出效應(yīng)產(chǎn)生的驅(qū)動(dòng)因素，分別以中國大陸東部省區(qū)、中部省區(qū)、西部省區(qū)的相對生態(tài)效率REEE為因變量、以經(jīng)過標(biāo)準(zhǔn)化處理和共線性診斷的相應(yīng)年度生態(tài)效率影響指標(biāo)X01—X10為自變量構(gòu)建偏最小二乘回歸模型。應(yīng)用SIMCA-P 11.5軟件分析得到東、中、西部3地區(qū)PLS模型的主成分個(gè)數(shù)分別為2、1、3個(gè)，Qh2值分別為0.7128、0.8406、0.9774，R2Y值分別為0.8636、0.8836、0.9082——表明經(jīng)過首次主成分提取的3個(gè)模型交叉有效性(Qh2值)及其對因變量的解釋能力(R2Y值)較高[41]，各地區(qū)PLS模型的擬合效果較好(誤差小于0.5并且不存在特異點(diǎn))。據(jù)此得出各模型的標(biāo)準(zhǔn)回歸系數(shù)及自變量投影重要性(VIP值)，結(jié)果分別如圖3所示。

圖3　各地區(qū)PLS模型的標(biāo)準(zhǔn)回歸系數(shù)及自變量投影重要性Fig.3　The standard regression coefficient and the independent variable of the projection importance of PLS model

從圖3所示的中國大陸東、中、西部省區(qū)相對生態(tài)效率的PLS模型回歸系數(shù)(經(jīng)歸一劃處理的標(biāo)準(zhǔn)回歸系數(shù))可以看出：對各地區(qū)相對生態(tài)-經(jīng)濟(jì)產(chǎn)出效率起正向作用的自變量個(gè)數(shù)介于7—8個(gè)之間，起負(fù)向作用的自變量個(gè)數(shù)為2—3個(gè)；各地區(qū)的指標(biāo)X01、X06、X10的VIP值均大于1，表明人均國民生產(chǎn)總值、單位面積內(nèi)的道路交通里程以及科研投入占GDP的比重等3個(gè)指標(biāo)對提高各地區(qū)生態(tài)-經(jīng)濟(jì)產(chǎn)出效率均有重要促進(jìn)作用，提升這些指標(biāo)可以顯著驅(qū)動(dòng)區(qū)域生態(tài)效率向周邊地區(qū)溢出(回歸系數(shù)為正值)，從而帶動(dòng)相鄰地區(qū)生態(tài)效率的提高，因此研究時(shí)段內(nèi)生態(tài)效率的低-低集聚區(qū)范圍逐漸縮小以至消失(圖2)。

除上述共同特征外，圖3顯示部分指標(biāo)對東、中、西部地區(qū)的作用差異也比較明顯：東部地區(qū)技術(shù)市場成交額占GDP的比重(X05)顯著驅(qū)動(dòng)了區(qū)域相對生態(tài)效率提高，體現(xiàn)了技術(shù)效率因素對該地區(qū)的重要作用；而外商直接投資額占GDP的比重(X02)及環(huán)保資金投入(X07)卻對該地區(qū)生態(tài)效率提升起負(fù)作用(回歸系數(shù)為負(fù)值)，表明外商直接投資(FDI)可能已經(jīng)使“污染物避難所”假說在中國東部地區(qū)成立[48]，環(huán)保資金的大量投入使得經(jīng)濟(jì)發(fā)展對此類投入產(chǎn)生依賴，反而不利于區(qū)域生態(tài)效率提升。除了指標(biāo)X02之外，進(jìn)出口總額(X03)也在一定程度上制約了中部地區(qū)相對生態(tài)效率的提高，表明以資源消耗為主體的外貿(mào)產(chǎn)業(yè)結(jié)構(gòu)及FDI對該地區(qū)生態(tài)效率的負(fù)面貢獻(xiàn)均相對較高，未來應(yīng)當(dāng)大力發(fā)展先進(jìn)制造業(yè)以逐漸消除結(jié)構(gòu)性污染[28]。與東部地區(qū)相對的是，F(xiàn)DI和環(huán)保資金投入驅(qū)動(dòng)西部地區(qū)相對生態(tài)效率提高的作用顯著，其中后者尤為顯著，體現(xiàn)了西部經(jīng)濟(jì)發(fā)展、生態(tài)環(huán)境改善對資金的高度依賴，究其原因可能與東、西部地區(qū)所處的經(jīng)濟(jì)發(fā)展階段差異有關(guān)；而除了指標(biāo)X03之外，城鎮(zhèn)人口比重(X04)對西部地區(qū)相對生態(tài)效率的提高亦起負(fù)向作用，表明該地區(qū)尚處于環(huán)境庫茲涅茨曲線的上升期、即城市化水平的提高使得區(qū)域生態(tài)效率相對下降[27]。

回顧圖1—圖3所得結(jié)果可以看出：各地區(qū)經(jīng)濟(jì)社會(huì)背景不同、發(fā)展階段各異，因此驅(qū)動(dòng)和制約區(qū)域生態(tài)效率提高的因素不同，導(dǎo)致研究時(shí)段內(nèi)中國大陸省區(qū)生態(tài)-經(jīng)濟(jì)產(chǎn)出效率呈現(xiàn)由東向西依次降低的三級梯度特征；近年來東部地區(qū)生態(tài)效率提升緩慢，但其溢出效應(yīng)已經(jīng)帶動(dòng)了中西部省區(qū)尤其是重慶等西部8省區(qū)的生態(tài)效率持續(xù)提升；未來，進(jìn)一步提升中國大陸地區(qū)生態(tài)效率的潛力在“胡煥庸線”西北部的省區(qū)，但有關(guān)地區(qū)不宜盲目照搬東、中部的發(fā)展經(jīng)驗(yàn)，應(yīng)當(dāng)依據(jù)區(qū)域自然地理?xiàng)l件及其所處的特定發(fā)展階段，制定因地制宜、因時(shí)制宜的生態(tài)效率提升策略，以期循序漸進(jìn)地解決區(qū)域人口、資源、環(huán)境與發(fā)展的矛盾問題。

3結(jié)論與討論

生態(tài)效率是學(xué)界聚焦的熱點(diǎn)問題，探討不同時(shí)空尺度的生態(tài)效率格局、演變趨勢及驅(qū)動(dòng)力具有理論意義和現(xiàn)實(shí)意義。本文在參考和借鑒有關(guān)生態(tài)系統(tǒng)服務(wù)價(jià)值理論研究成果的基礎(chǔ)上提出了生態(tài)-經(jīng)濟(jì)產(chǎn)出效率概念內(nèi)涵及其表征方法，評價(jià)過程中兼顧了經(jīng)濟(jì)發(fā)展與資源消耗、環(huán)境損失代價(jià)，因此相對于以往側(cè)重污染物分析的生態(tài)效率研究更全面；為了提高評價(jià)結(jié)果的可靠性，本文引入重心模型分析了2001—2013年中國大陸31個(gè)省區(qū)生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的重心移動(dòng)方向及其移動(dòng)距離、通過4種空間賦權(quán)方法分析了評價(jià)單元的全局空間依賴性，據(jù)此分析其時(shí)序演變趨勢、空間分布格局及其驅(qū)動(dòng)因素，初步得到幾點(diǎn)結(jié)論如下：

(1)研究時(shí)段內(nèi)中國大陸省際生態(tài)-經(jīng)濟(jì)產(chǎn)出效率整體提升，自西向東依次提高的三級梯度分布格局基本形成，但中部8省區(qū)及甘肅、新疆、青海、西藏等西部4省區(qū)生態(tài)效率提升速度相對緩慢，東部省區(qū)及西部其它8省區(qū)提升速度較快。

(2)研究時(shí)段內(nèi)中國大陸省區(qū)生態(tài)-經(jīng)濟(jì)產(chǎn)出效率具有全局空間依賴性和空間正相關(guān)特征，而局部分異特征同時(shí)存在——顯著高-高集聚和顯著低-高集聚表現(xiàn)相對穩(wěn)定，高-低集聚區(qū)的生態(tài)效率溢出效應(yīng)促使顯著低-低集聚區(qū)范圍逐漸縮小以至消失。

(3)人均國民生產(chǎn)總值、單位國土面積內(nèi)的道路交通里程數(shù)、科研資金投入占GDP的比重等3個(gè)因素顯著驅(qū)動(dòng)了中國大陸省區(qū)生態(tài)-經(jīng)濟(jì)產(chǎn)出效率整體提升；制約各地區(qū)生態(tài)效率溢出或提升的主要因素不盡相同，東、中、西部地區(qū)依次為外商直接投資額和環(huán)境保護(hù)與治理資金投入額、外商直接投資額和進(jìn)出口總額、進(jìn)出口總額和城鎮(zhèn)人口比例。

本文以謝高地等提出的中國陸地生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子表為基礎(chǔ)計(jì)算得出大陸31個(gè)省區(qū)歷年生態(tài)系統(tǒng)服務(wù)價(jià)值，與以往的影子工程法、機(jī)會(huì)成本法等相比，計(jì)算過程中分別依據(jù)有關(guān)省區(qū)歷年糧食單產(chǎn)、原糧市場價(jià)格等參數(shù)進(jìn)行客觀核算，可在一定程度上避免不確定性，所得結(jié)果與前人同類研究成果比較接近[49- 50]，進(jìn)一步表明該方法具有一定的推廣應(yīng)用價(jià)值。需要指出的是：本文將各類土地利用類型統(tǒng)一劃分為耕地等6類生態(tài)系統(tǒng)，實(shí)際上并不能涵蓋所有生態(tài)系統(tǒng)類型；而將不同省區(qū)同一生態(tài)系統(tǒng)類型的當(dāng)量因子取等值的處理方法也可能忽略局地生態(tài)系統(tǒng)生物量的空間異質(zhì)性[51- 52]，可能由此導(dǎo)致生態(tài)-經(jīng)濟(jì)產(chǎn)出效率的評價(jià)結(jié)果與客觀事實(shí)略有出入，后續(xù)研究將持續(xù)關(guān)注生態(tài)系統(tǒng)服務(wù)價(jià)值評價(jià)方法的改進(jìn)與修正問題。此外，因本文所界定的生態(tài)效率內(nèi)涵、應(yīng)用的研究方法與以往研究存在差異，加之研究視角、選取的評價(jià)指標(biāo)不同，可能導(dǎo)致得出的生態(tài)效率影響因素與以往研究不盡相同——但是這一結(jié)果愈加體現(xiàn)了區(qū)域生態(tài)效率影響因素眾多、影響機(jī)制復(fù)雜的特征，進(jìn)一步表明應(yīng)用不同方法、基于不同時(shí)空尺度開展多角度生態(tài)效率評價(jià)研究的必要性。

參考文獻(xiàn)(References):

[1]Schaltegger S, Burritt R. Contemporary Environmental Accounting: Issues, Concept and Practice. Sheffield: Greenleaf Publishing, 2000.

[2]呂彬, 楊建新. 生態(tài)效率方法研究進(jìn)展與應(yīng)用. 生態(tài)學(xué)報(bào), 2006, 26(11): 3898- 3906.

[3]尹科, 王如松, 周傳斌, 梁菁. 國內(nèi)外生態(tài)效率核算方法及其應(yīng)用研究述評. 生態(tài)學(xué)報(bào), 2012, 32(11): 3595- 3605.

[4]WBCSD. Eco-efficiency Leadership for Improved Economic and Environmental Performance. Geneva: WBCSD, 1996.

[5]Organization for Economic Cooperation and Development. Eco-Efficiency. Paris: OECD, 1998.

[6]MOLL S, GEE D. Making sustainability accountable: Eco-efficiency, resource productivity and innovation. Proceedings of a Workshop on the Fifth Anniversary of the European Environment Agency (EEA), Copenhagen, 1998: 1- 37.

[7]Nieminen E, Linke M, Tobler M, Beke B V. EU COST Action 628: life cycle assessment (LCA) of textile products, eco-efficiency and definition of best available technology (BAT) of textile processing. Journal of Cleaner Production, 2007, 15(13- 14): 1259- 1270.

[8]Figge F, Hahn T. Sustainable value added—measuring corporate contributions to sustainability beyond eco-efficiency. Ecological Economics, 2004, 48(2): 173- 187.

[9]Kristina D, Paul E. Eco-efficiency trends in the UK steel and aluminum industries: differences between resource efficiency and resource productivity. Journal of Industrial Ecology, 2005, 9(4): 171- 188.

[10]Van Caneghem J, Block C, Cramm P, Mortier R, Vandecasteele C. Improving eco-efficiency in the steel industry: the Arcelormittal Gent case. Journal of Cleaner Production, 2010, 18(8): 807- 814.

[11]Korhonen P J, Luptacik M. Eco-efficiency analysis of power plants: an extension of data envelopment analysis. European Journal of Operational Research, 2004, 154(2): 437- 446.

[12]Dyckhoff H, Allen K. Measuring ecological efficiency with data envelopment analysis. European Journal of Operational Research, 2001, 132(2): 312- 325.

[13]Van Middelaar C E, Berentsen P B M, Dolman M A, de Boer I J M. Eco-efficiency in the production chain of Dutch semi-hard cheese. Livestock Science, 2011, 139(1/2): 91- 99.

[14]Erkko S, Melanen M, Mickwitz P. Eco-efficiency in the Finnish EMAS reports—a buzz word? Journal of Cleaner Production, 2005, 13(8): 799- 813.

[15]Aoe T. Eco-efficiency and ecodesign in electrical and electronic products. Journal of Cleaner Production, 2007, 15(15): 1406- 1414.

[16]劉晶茹, 呂彬, 張娜, 石垚. 生態(tài)產(chǎn)業(yè)園的復(fù)合生態(tài)效率及評價(jià)指標(biāo)體系. 生態(tài)學(xué)報(bào), 2014, 34(1): 136- 141.

[17]季丹. 中國區(qū)域生態(tài)效率評價(jià)——基于生態(tài)足跡方法. 當(dāng)代經(jīng)濟(jì)管理, 2013, 35(2): 57- 62.

[18]芮俊偉, 周貝貝, 錢誼, 孫靖. 生態(tài)工業(yè)園區(qū)生態(tài)效率評估方法研究及應(yīng)用. 生態(tài)與農(nóng)村環(huán)境學(xué)報(bào), 2013, 29(4): 466- 470.

[19]潘丹, 應(yīng)瑞瑤. 中國農(nóng)業(yè)生態(tài)效率評價(jià)方法與實(shí)證——基于非期望產(chǎn)出的SBM模型分析. 生態(tài)學(xué)報(bào), 2013, 33(12): 3837- 3845.

[20]李健, 鄧傳霞, 張松濤. 基于非參數(shù)距離函數(shù)法的區(qū)域生態(tài)效率評價(jià)及動(dòng)態(tài)分析. 干旱區(qū)資源與環(huán)境, 2015, 29(4): 19- 23.

[21]張淑英, 李德山, 劉媛媛. 區(qū)域工業(yè)生態(tài)效率評價(jià)及其影響因素研究. 統(tǒng)計(jì)與決策, 2013, (3): 61- 64.

[22]陳作成, 龔新蜀. 西部地區(qū)產(chǎn)業(yè)系統(tǒng)生態(tài)效率測評與影響因素分析. 中國科技論壇, 2013, (10): 49- 55.

[23]陳傲. 中國區(qū)域生態(tài)效率評價(jià)及影響因素實(shí)證分析——以2000- 2006年省際數(shù)據(jù)為例. 中國管理科學(xué), 2008, 16(S1): 566- 570.

[24]汪克亮, 孟祥瑞, 程云鶴. 技術(shù)的異質(zhì)性、節(jié)能減排與地區(qū)生態(tài)效率——基于2004—2012年中國省際面板數(shù)據(jù)的實(shí)證分析. 山西財(cái)經(jīng)大學(xué)學(xué)報(bào), 2015, 37(2): 69- 80.

[25]崔瑋, 苗建軍, 楊晶. 考慮碳排放的城市非農(nóng)用地生態(tài)效率的實(shí)證研究. 系統(tǒng)工程, 2012, 30(12): 10- 18.

[26]初善冰, 黃安平. 外商直接投資對區(qū)域生態(tài)效率的影響——基于中國省際面板數(shù)據(jù)的檢驗(yàn). 國際貿(mào)易問題, 2012, (11): 128- 144.

[27]羅能生, 李佳佳, 羅富政. 中國城鎮(zhèn)化進(jìn)程與區(qū)域生態(tài)效率關(guān)系的實(shí)證研究. 中國人口·資源與環(huán)境, 2013, 23(11): 53- 60.

[28]潘興俠, 何宜慶. 中部六省生態(tài)效率評價(jià)及其與產(chǎn)業(yè)結(jié)構(gòu)的時(shí)空關(guān)聯(lián)分析. 統(tǒng)計(jì)與決策, 2015, (3): 127- 130.

[29]潘興俠. 我國區(qū)域生態(tài)效率評價(jià)、影響因素及收斂性研究[D]. 南昌: 南昌大學(xué), 2014.

[30]程翠云, 任景明, 王如松. 我國農(nóng)業(yè)生態(tài)效率的時(shí)空差異. 生態(tài)學(xué)報(bào), 2014, 34(1): 142- 148.

[31]楊凱, 王要武, 薛維銳. 區(qū)域梯度發(fā)展模式下我國工業(yè)生態(tài)效率區(qū)域差異與對策. 系統(tǒng)工程理論與實(shí)踐, 2013, 33(12): 3095- 3101.

[32]王恩旭, 武春友. 基于超效率DEA模型的中國省際生態(tài)效率時(shí)空差異研究. 管理學(xué)報(bào), 2011, 8(3): 443- 450.

[33]Johnston R J, Russell M. An operational structure for clarity in ecosystem service values. Ecological Economics, 2011, 70(12): 2243- 2249.

[34]M?ler K G, Aniyar S, Jansson ?. Accounting for ecosystem services as a way to understand the requirements for sustainable development. Proceedings of the National Academy of Sciences of the United States of America, 2008, 105(28): 9501- 9506.

[35]謝高地, 魯春霞, 冷允法, 鄭度, 李雙成. 青藏高原生態(tài)資產(chǎn)的價(jià)值評估. 自然資源學(xué)報(bào), 2003, 18(2): 189- 196.

[36]Tobler W R. A computer movie simulating urban growth in the Detroit region. Economic Geography, 1970, 46: 234- 240.

[37]Goodchild M Fs. The validity and usefulness of laws in geographic information science and geography. Annals of the Association of American Geographers, 2004, 94(2): 300- 303.

[38]楊華, 李小文, 施宏, 趙開廣, 韓麗娟. SARS沿交通線的“飛點(diǎn)”傳播模型. 遙感學(xué)報(bào), 2003, 7(4): 251- 255.

[39]李國平, 王春楊. 我國省域創(chuàng)新產(chǎn)出的空間特征和時(shí)空演化——基于探索性空間數(shù)據(jù)分析的實(shí)證. 地理研究, 2012, 31(1): 95- 106.

[40]萬魯河, 王紹巍, 陳曉紅. 基于GeoDA的哈大齊工業(yè)走廊GDP空間關(guān)聯(lián)性. 地理研究, 2011, 30(6): 977- 984.

[41]于松青, 林盛. 基于偏最小二乘回歸的山東省電力需求預(yù)測分析. 干旱區(qū)資源與環(huán)境, 2015, 29(2): 14- 20.

[42]中華人民共和國中央人民政府. 糧食局等部門部署2012年玉米和大豆臨時(shí)收儲(chǔ)工作. [2015-07-01]. http://www.gov.cn/gzdt/2012- 11/21/content_2271750.htm.

[43]中華人民共和國國家發(fā)展與改革委員會(huì). 關(guān)于印發(fā)2012年小麥最低收購價(jià)執(zhí)行預(yù)案的通知. [2015-07-01]. http://www.gov.cn/zwgk/2012-05/24/content_2143976.htm

[44]中華人民共和國國家發(fā)展與改革委員會(huì). 關(guān)于印發(fā)2012年中晚稻最低收購價(jià)執(zhí)行預(yù)案的通知. [2015-07-01]. http://www.gov.cn/gzdt/2012-09/04/content_2216235.htm

[45]中華人民共和國國家發(fā)展與改革委員會(huì). 7月國內(nèi)市場糧食價(jià)格. 國家發(fā)改委網(wǎng)站, [2015-07-01]. http://www.gov.cn/zfjg/content_726858.htm

[46]中國國家統(tǒng)計(jì)局. 2002—2011年中國統(tǒng)計(jì)年鑒. 國家統(tǒng)計(jì)局網(wǎng)站, [2015-07-01]. http://www.stats.gov.cn/tjsj/ndsj.

[47]中華人民共和國國土資源部. 中國國土資源統(tǒng)計(jì)年鑒. 北京: 地質(zhì)出版社, 2012.

[48]蘇梽芳, 廖迎, 李穎. 是什么導(dǎo)致了“污染天堂”: 貿(mào)易還是FDI?——來自中國省級面板數(shù)據(jù)的證據(jù). 經(jīng)濟(jì)評論, 2011, (3): 97- 104, 116- 116.

[49]陳仲新, 張新時(shí). 中國生態(tài)系統(tǒng)效益的價(jià)值. 科學(xué)通報(bào), 2000, 45(1): 17- 22.

[50]鄭德鳳, 臧正, 趙良仕, 孫才志. 中國省際資源環(huán)境成本及生態(tài)負(fù)荷強(qiáng)度的時(shí)空演變分析. 地理科學(xué), 2014, 34(6): 672- 680.

[51]謝高地, 魯春霞, 成升魁. 全球生態(tài)系統(tǒng)服務(wù)價(jià)值評估研究進(jìn)展. 資源科學(xué), 2001, 23(6): 5- 9.

[52]謝高地, 甄霖, 魯春霞, 肖玉, 陳操. 一個(gè)基于專家知識(shí)的生態(tài)系統(tǒng)服務(wù)價(jià)值化方法. 自然資源學(xué)報(bào), 2008, 23(5): 911- 919.

基金項(xiàng)目:國家基礎(chǔ)科學(xué)人才培養(yǎng)基金項(xiàng)目(J1103408)

收稿日期:2015- 06- 02;

修訂日期:2015- 11- 30

*通訊作者

Corresponding author.E-mail: zouxq@nju.edu.cn

DOI:10.5846/stxb201506021108

Analysis on spatial and temporal pattern and the driving forces of ecological-economic efficiency for provinces on Chinese mainland

ZANG Zheng1, ZOU Xinqing1,2,*

1SchoolofGeographicandOceanographicSciences,NanjingUniversity,Nanjing210023,China2CollaborativeInnovationCenterofSouthChinaSeaStudies,Nanjing210093,China

Abstract：Referring to the research findings of ecosystem service value, a new concept “ecological-economic output” was proposed and defined in this paper the connotation and characterization method of this concept was also given. With the statistical data, ecological- economic output efficiency of 31 provinces of Chinese mainland during 2001—2013 was calculated. The spatial-temporal distribution, gravity center transferring direction and distance were also analyzed. On this base, the spatial agglomeration trait and the driving forces were discussed through spatial data analysis method and partial least squares regression model. The results showed that in Chinese mainland, during the research period, ecological-economic output efficiency at province level showed three-level graded distribution, characterized by increasing from the west to the east. Ecological- economic output efficiency had the trait of overall spatial dependence and spatial positive correlation. Remarkable high-high agglomeration and low-high agglomeration in some regions were stable, and ecological overflow effect of high-low agglomeration made low-low agglomeration vanished eventually. GDP per capita, road lengths per unit area and proportion of scientific research funding from GDP drove the improvement of ecological- economic output efficiency of Chinese mainland. However, the restricting elements were varied from place to place.

Key Words:ecological economics; ecological efficiency; empirical analysis; spatial and temporal pattern; driving forces

臧正，鄒欣慶.中國大陸省際生態(tài)-經(jīng)濟(jì)效率的時(shí)空格局及其驅(qū)動(dòng)因素.生態(tài)學(xué)報(bào),2016,36(11):3300- 3311.

Zang Z, Zou X Q.Analysis on spatial and temporal pattern and the driving forces of ecological-economic efficiency for provinces on Chinese mainland.Acta Ecologica Sinica,2016,36(11):3300- 3311.