李玲玲,韓瑞玲,,*,張曉燕

1 河北師范大學(xué)地理科學(xué)學(xué)院, 石家莊 050024 2 河北師范大學(xué)家政學(xué)院, 石家莊 050024

航空碳排放主要源自于航空燃油在燃燒時(shí)與空氣混合后排出的CO2、SO2、CO、氮氧化物(Nitrogen Oxide,NOx)、未燃碳?xì)浠衔?Unburned Hydrocarbons,UHC)、碳煙顆粒(Soot)、細(xì)顆粒物(Fine Particulate Matter,PM)、飛行尾跡以及形成的卷云(Aviation Induced Cloudiness,AIC)等污染物[1],不僅影響空氣質(zhì)量且可造成溫室效應(yīng),其中CO2排放量占總排放量的87%[2]?；趪H民用航空組織(International Civil Aviation Organization,ICAO)標(biāo)準(zhǔn)碳排放模型發(fā)現(xiàn),1990年以來國際航空碳排放呈增長趨勢；基于歐洲環(huán)境署(European Environment Agency,EEA)公布的碳排放量計(jì)算方法,預(yù)計(jì)到2025年,國際航空碳排放量將比2005年水平高出110%以上[3]。對中國而言,隨著空中交通需求擴(kuò)展以及民用航空器數(shù)量增長,航空碳排放還將處于持續(xù)增長態(tài)勢[4]。航空運(yùn)輸碳排放占全球人為排放碳總量的2.0%—2.5%[5],中國航空碳排放量占全球人為碳排放量的0.22%,航空碳排放問題備受關(guān)注[6]。中國民航一直致力于低碳綠色發(fā)展,為2060年前國家實(shí)現(xiàn)碳中和目標(biāo)做努力。

針對航空碳排放進(jìn)行有效減排探索是學(xué)者們一直努力的方向。具體分為3點(diǎn)：①使用替代能源促進(jìn)航空碳減排。使用高效率替代燃料或生物燃料是民航業(yè)普遍承認(rèn)的減排措施,但是尚未發(fā)揮作用且未來發(fā)展趨勢不明朗[7]。Staples[8]指出使用可持續(xù)替代航空燃料可以在2050年之前將航空全生命周期的溫室氣體排放量最多減少68.1%,但造價(jià)昂貴、實(shí)現(xiàn)周期較長[9],并且受到生物生產(chǎn)性土地資源[10]的限制而可行性較低[11]。Marais也表明新技術(shù)和可持續(xù)燃料需要更多努力和時(shí)間來發(fā)展和實(shí)施[12]。②有效利用交通結(jié)構(gòu)因素的影響,在航線規(guī)劃、設(shè)計(jì)、管理等方面促進(jìn)結(jié)構(gòu)優(yōu)化,對航空碳減排將大有裨益。調(diào)整終端區(qū)資源分配[13]和機(jī)場的操作程序、設(shè)計(jì)最佳飛機(jī)航線以及優(yōu)化特定跑道和航線的分配等措施盡管效用較小,但可以在短期內(nèi)進(jìn)行操作更改,是目前實(shí)現(xiàn)航空碳減排的有效措施,一般“最清潔”的航線為短途、直達(dá)航線[14]。③政策措施促進(jìn)航空碳減排。Fukui等研究了美國航空燃油稅增加4.3美分[15],將在短期內(nèi)(增稅后1年)減少約0.14%—0.18%碳排放量,González等也證實(shí)了該方法的有效性[16]。2019年,ICAO開始執(zhí)行國際航空碳抵消和減排計(jì)劃(Carbon Offsetting and Reduction Scheme for International Aviation,CORSIA)[17],要求所有重要的航空公司強(qiáng)制記錄國際航班的碳排放量,并從2021年開始執(zhí)行超過排放基準(zhǔn)的碳增長付款政策。

相較而言,在目前的技術(shù)水平條件下,提高航空碳排放效率可行性高于全面使用替代能源,復(fù)雜程度小于航線結(jié)構(gòu)的調(diào)整和減排政策的制定,是行之有效的促進(jìn)航空碳減排手段[18]。航空碳排放效率是影響航空碳排放的重要因素,其主要通過提高飛行器能源消費(fèi)強(qiáng)度而降低碳排放[19]。測度航空碳排放效率的定量研究方法有成本函數(shù)評(píng)估(Cost Function)[20]、隨機(jī)前沿模型(Stochastic Frontier Model,SFA)[21]和數(shù)據(jù)包絡(luò)分析方法(Data Envelopment Analysis,DEA)。其中,DEA使用數(shù)學(xué)規(guī)劃模型評(píng)價(jià)航空碳排放輸入和輸出單位間的相對有效性,經(jīng)不斷演變,現(xiàn)以基于松弛變量模型(Slack Based Measure,SBM)使用最為廣泛[22]。Schefczvk首先將DEA方法引入航空領(lǐng)域效率研究[23]。Arjomandi等[24]應(yīng)用DEA模型,研究了2007—2010年世界48家主要航空公司的環(huán)境和技術(shù)效率,發(fā)現(xiàn)中國和北亞的航空公司技術(shù)效率最高,歐洲的航空公司環(huán)保效率表現(xiàn)最好；低成本航空公司更注重環(huán)保效率,并在技術(shù)上實(shí)現(xiàn)規(guī)模收益遞增。Cui等[25]提出了一種基于虛擬邊界的動(dòng)態(tài)SBM模型,克服了傳統(tǒng)動(dòng)態(tài)SBM模型在識(shí)別有效決策單元方面的不足,并認(rèn)為人均GDP、飛機(jī)平均使用年限、平均運(yùn)輸距離對航空公司能源效率有顯著影響。Wang等[4]采用考慮非期望輸出的全局松弛變量模型(Global Slack Based Measure,GSBM)計(jì)算航空碳排放的靜態(tài)效率,利用全局馬姆奎斯特-盧恩伯格(Global Malmquist-Luenberger,GML)生產(chǎn)率指數(shù)結(jié)合分析了2009—2013年中國13家航空公司的航空碳排放效率,結(jié)果表明中國航空公司碳排放靜態(tài)效率呈倒U型,且碳排放動(dòng)態(tài)效率各不相同。Liu等[26]采用DEA和GML指數(shù)計(jì)算了中國12家航空公司2007—2013年的民航碳排放績效指數(shù),其碳排放效率通過技術(shù)進(jìn)步提高了11.93%,并認(rèn)為航線分布也影響碳排放效率。

可見,已有的宏觀航空碳排放研究多是針對其減排措施及具體影響因素開展的,碳排放效率研究則多集中于微觀的航空公司視角,且多是將研究單元視為獨(dú)立均質(zhì)的個(gè)體,忽略了研究單元碳排放的異質(zhì)性和關(guān)聯(lián)性。中國幅員遼闊,各區(qū)域民航產(chǎn)業(yè)發(fā)展不均衡,航空碳排放在各省域表現(xiàn)也不盡相同。因此有效分析中國航空碳排放及其效率的時(shí)空分異特征,并針對性的提供碳減排方略是需要積極探索的議題。為此,本文首先利用可拓展隨機(jī)性環(huán)境影響評(píng)估模型(Stochastic Impacts by Regression on Population Affluence and Technology,STIRPAT)計(jì)算中國航空碳排放總量,并對碳排放量的時(shí)空演化特征進(jìn)行了判定,進(jìn)而探明不同尺度區(qū)域航空碳排放的空間差異與空間關(guān)聯(lián)表現(xiàn)。在此基礎(chǔ)上,充分考慮空間異質(zhì)性,利用基于非期望產(chǎn)出的超效率SBM模型和標(biāo)準(zhǔn)差橢圓(Standard Deviational Ellipse,SDE)分析法,對中國省域航空碳排放效率進(jìn)行了時(shí)空演變特征分析,并對各省域航空碳排放效率的空間分布及空間格局演化過程進(jìn)行定量描述,進(jìn)而提出差異化的航空碳減排方略。從全國、區(qū)域、省域3個(gè)尺度進(jìn)行航空碳排放及效率時(shí)空差異研究,為構(gòu)建資源消耗低、環(huán)境污染少的民航綠色產(chǎn)業(yè)結(jié)構(gòu)和生產(chǎn)方式,實(shí)現(xiàn)民航綠色循環(huán)低碳發(fā)展,減少中國航空碳排放做積極努力。

1 研究方法與數(shù)據(jù)來源

1.1 STIRPAT模型

STIRPAT模型是研究碳排放問題的常用方法,是IPAT等式經(jīng)過多次改進(jìn)形成的隨機(jī)形式,其允許各影響因素非單調(diào)、不同比例地變化[27],可實(shí)現(xiàn)對各種驅(qū)動(dòng)因素影響程度的分析[28]。本文構(gòu)建的中國航空碳排放總量的STIRPAT模型如下：

(1)

式中,Ii表示第i年中國國內(nèi)航空碳排放總量(t),Pik表示第i年k省的航空運(yùn)輸總量(t),Ai表示第i年每千噸航空運(yùn)輸量運(yùn)送的距離(km),Ti表示第i年每千噸公里油耗(t/103km),e為航空煤油的碳排放因子(kg/kg)。

1.2 區(qū)域碳排放指數(shù)模型

區(qū)域碳排放指數(shù)模型可在考慮外部條件下對各省域航空碳排放區(qū)域差異進(jìn)行分析。通過計(jì)算由人均碳排放指數(shù)(Wp)和單位面積的碳排放指數(shù)(Wa)2個(gè)概念構(gòu)建的區(qū)域碳排放指數(shù)模型得到區(qū)域碳排放指數(shù)[29],公式如下：

(2)

式中,W為碳排放指數(shù)；Wp為人均碳排放指數(shù),Wp=c/cg；Wa為單位面積的碳排放指數(shù),Wa=Ca/Cg；c和Ca分別為人均和單位面積的碳排放量(t,t/km2),cg和Cg分別為應(yīng)對全球氣候變化目標(biāo)的人均、單位面積碳排放量(t,t/km2),按照溫室氣體濃度限制在500 mL/m3當(dāng)量的碳排放軌跡以及全球人居面積8.58×107km2進(jìn)行測算,分別取值2 t、2.33 t/km2；Wp.max、Wa.max是全球人均溫室氣體排放指數(shù)、單位面積溫室氣體排放指數(shù)最大值,基于世界資源研究所(World Resources Institute,WRI)公布的全球溫室氣體排放數(shù)據(jù),為使結(jié)果區(qū)分度更明顯,將Wp.max、Wa.max定為15和35；區(qū)域碳排放等級(jí)評(píng)估按照W<0.20、0.20≤W<0.50、W≥0.50的劃分標(biāo)準(zhǔn)對應(yīng)低排放、中排放、高排放3個(gè)等級(jí)[29]。

1.3 Moran′s I指數(shù)

Moran′s I指數(shù)分析方法可用于統(tǒng)計(jì)分析航空碳排放空間分布特征和揭示空間分布差異,包括全局Moran′s I指數(shù)和局部Moran′s I指數(shù)。全局Moran′s I指數(shù)能有效反應(yīng)區(qū)域整體集聚情況。局部Moran′s I指數(shù)能識(shí)別局部區(qū)域集聚特性,結(jié)合局部Moran′s I散點(diǎn)圖或集聚圖等形式揭示空間異質(zhì)性規(guī)律[30]。公式如下：

(3)

(4)

(5)

其中,

(6)

Ii為第i個(gè)空間單元的局部Moran′s I指數(shù),wij為空間權(quán)重值,a為空間單元總數(shù)。

1.4 基于非期望產(chǎn)出的超效率SBM模型

在計(jì)算得到中國航空碳排放總量的基礎(chǔ)上,選擇考慮非期望產(chǎn)出的超效率SBM模型對中國航空碳排放效率進(jìn)行測算。SBM模型在目標(biāo)函數(shù)中加入松弛變量,有效解決了將航空碳排放總量作為非期望產(chǎn)出時(shí),在效率評(píng)價(jià)中投入、產(chǎn)出的松弛性問題[31],但是存在著決策單元同時(shí)為1的情況,不利于決策單元的相互比較。為此,超效率SBM模型通過修正松弛變量可解決該問題[32]。具體公式如下：

(7)

式中,ρ*為目標(biāo)效率值；x、yd和yu分別為投入、期望產(chǎn)出和非期望產(chǎn)出；向量x-、y-d和y-u分別為投入松弛量、期望的產(chǎn)出松弛量和非期望產(chǎn)出的松弛量；λ為權(quán)重向量；k表示被評(píng)價(jià)的決策單元(DMU)。ρ*=1表示決策單元完全有效,否則表明存在效率損失,需要對投入、產(chǎn)出量進(jìn)行調(diào)整。效率值具體又分為4個(gè)級(jí)別：極低效率(0<ρ*≤0.3),低效率(0.3<ρ*≤0.6),中效率(0.6<ρ*≤0.9),高效率(ρ*>0.9)[32]。

1.5 標(biāo)準(zhǔn)差橢圓分析法

標(biāo)準(zhǔn)差橢圓分析法是空間統(tǒng)計(jì)分析中研究方向分布的經(jīng)典方法之一。本文使用該方法計(jì)算中國航空碳排放效率的空間中心點(diǎn)、扁率、長軸和短軸、方位角等核心參數(shù),探究其空間格局演化過程,其中中心點(diǎn)的移動(dòng)代表中國航空碳排放效率空間分布中心的變化,扁率、長軸和短軸分別代表中國航空碳排放效率的空間分布強(qiáng)度、分布范圍,方位角代表中國航空碳排放效率的空間分布主要發(fā)展趨勢方向[33],其公式分別為：

(8)

(9)

(10)

式中,SEDx和SEDy分別是(xi,yi)距離分布重心的相對坐標(biāo)；θ是橢圓由正北順時(shí)旋轉(zhuǎn)的角度。

(11)

(12)

式中,σx和σy分別為x軸和y軸的標(biāo)準(zhǔn)差。

1.6 研究區(qū)域與研究數(shù)據(jù)

本文選擇全國31個(gè)省級(jí)行政單元(尚未對中國港澳臺(tái)地區(qū)數(shù)據(jù)進(jìn)行統(tǒng)計(jì))為基本研究范圍。由于統(tǒng)計(jì)資料在2004年前未單獨(dú)統(tǒng)計(jì)航空碳排放效率所需的資本投入數(shù)據(jù),為保持指標(biāo)統(tǒng)一性和數(shù)據(jù)分析連貫性,本文的研究期定為2004—2019年?；赟TIRPAT模型進(jìn)行的中國航空碳排放總量研究所需指標(biāo)主要包括“國內(nèi)航空運(yùn)輸總量、每千噸運(yùn)輸量運(yùn)送的距離、每千噸公里油耗”,其中國內(nèi)航空運(yùn)輸總量指標(biāo)為國內(nèi)旅客運(yùn)輸總量與貨郵運(yùn)輸總量之和,每千噸運(yùn)輸量運(yùn)送的距離和每千噸公里油耗指標(biāo)由國內(nèi)運(yùn)輸總周轉(zhuǎn)量、運(yùn)輸總量和航空煤油消耗量計(jì)算得出。數(shù)據(jù)來源于《從統(tǒng)計(jì)看民航》《民航行業(yè)發(fā)展統(tǒng)計(jì)公報(bào)》《民航機(jī)場生產(chǎn)統(tǒng)計(jì)公報(bào)》,航空煤油的碳排放因子e參照中國民航航空局發(fā)布的《民用航空企業(yè)(航空公司)溫室氣體排放報(bào)告補(bǔ)充數(shù)據(jù)表》,其缺省值定為3.15 kg/kg?；诜瞧谕a(chǎn)出的超效率SBM模型進(jìn)行的航空碳排放效率研究主要包括投入和產(chǎn)出指標(biāo)2個(gè)方面。其中投入指標(biāo)一般選取資本投入和勞動(dòng)力投入指標(biāo)[34]。本文資本投入指標(biāo)選取交通運(yùn)輸固定資產(chǎn)投入數(shù)據(jù)作為各省域航空碳排放效率的投入數(shù)據(jù),來源于《中國第三產(chǎn)業(yè)統(tǒng)計(jì)年鑒》；勞動(dòng)力投入指標(biāo)選取航空運(yùn)輸業(yè)從業(yè)人數(shù)數(shù)據(jù),來源于《中國人口和就業(yè)統(tǒng)計(jì)年鑒》。產(chǎn)出指標(biāo)選取航空運(yùn)輸總量作為期望產(chǎn)出、航空碳排放總量作為非期望產(chǎn)出,數(shù)據(jù)來源于《民航機(jī)場生產(chǎn)統(tǒng)計(jì)公報(bào)》及航空碳排放的計(jì)算結(jié)果。

2 航空碳排放時(shí)空演變特征分析

2.1 航空碳排放總量時(shí)序特征分析

2004—2019年,中國航空碳排放總量持續(xù)增長,由2.48×107t增至11.60×107t,年均排放量為6.09×107t(圖1)。航空碳排放增長主要受到航空運(yùn)輸量增長的影響,全行業(yè)運(yùn)輸周轉(zhuǎn)量由2014年的230.99×109t/km增至2019年的1293.25×109t/km,增長了4.6倍,年均增長11.37%。2004—2019年,中國航空碳排放增長率波動(dòng)下降,由30.4%降至6.38%,年均增長率為12.10%。2004年增長率在研究期內(nèi)最高,主要是受到2003年非典事件的滯后影響。其他增長率較低的年份分別受到不同事件影響,如2008年全球金融危機(jī)及燃油稅調(diào)整、2011年日本地震與敘利亞戰(zhàn)爭等突發(fā)事件以及全球通脹壓力加大導(dǎo)致國內(nèi)航空需求變?nèi)酰?017年主要受薩德事件影響嚴(yán)重,中國客貨運(yùn)多條航線運(yùn)力投放下降20%；2019年受到波音737MAX飛機(jī)停飛事件影響,各航空公司的運(yùn)力安排也相應(yīng)減少,航空運(yùn)輸量增速放緩使得航空碳排放量增速放緩。

2004—2019年,中國航空碳排放強(qiáng)度(航空碳排放量與GDP比值)波動(dòng)下降,由153.5 t/108元降至117.1 t/108元,年均強(qiáng)度為121.03 t/108元(圖2),年均增長率為-0.88%。其中,2004—2013年航空碳排放強(qiáng)度主要呈現(xiàn)明顯下降趨勢；2014—2019年航空碳排放強(qiáng)度主要呈現(xiàn)明顯增長趨勢,僅2017、2019年有小幅下降。但總體上看,中國航空業(yè)節(jié)能減排效果明顯。一方面是因?yàn)檠芯科趦?nèi)航空燃油單位能耗逐年降低(0.51降至0.44),有效抑制了航空碳排放強(qiáng)度的快速增長,這與前人研究結(jié)果一致[7]。另外一方面,是因?yàn)?008年歐洲碳排放交易體系(EU-ETS)[35]將中國33家航空公司納入后,中國為應(yīng)對減排壓力制定了多項(xiàng)措施,如2014年,國家發(fā)改委印發(fā)《國家應(yīng)對氣候變化規(guī)劃(2014—2020年)》指出中國將致力于航空產(chǎn)業(yè)低碳發(fā)展；2017年,中國民航局發(fā)布《中國民用航空發(fā)展第十三個(gè)五年規(guī)劃》,提出到2020年中國計(jì)劃實(shí)現(xiàn)碳強(qiáng)度降低40%—50%,并建成綠色民航標(biāo)準(zhǔn)體系；2018年,中國民航局發(fā)布《關(guān)于深入推進(jìn)民航綠色發(fā)展的實(shí)施意見》,進(jìn)一步促進(jìn)民航航空器節(jié)能減排等,均對抑制碳排放過快增長發(fā)揮了作用。

圖1 2004—2019年中國航空碳排放總量及增長率變化圖 Fig.1 Changes in China′s aviation carbon emissions and growth rate during 2004 to 2019

圖2 2004—2019年中國航空碳排放強(qiáng)度及油耗變化圖 Fig.2 Changes in China′s aviation carbon emission intensity and growth rate during 2004 to 2019

2.2 航空碳排放總量空間演變特征分析

2.2.1空間分異分析

圖3 2004—2019年四大經(jīng)濟(jì)區(qū)占全國總碳排放量的比例變化 Fig.3 The proportion of the four economic zones in the country′s total carbon emissions during 2004 to 2019

2004—2019年,航空碳排放總量在中國四大經(jīng)濟(jì)區(qū)的分布總體表現(xiàn)為“東部>西部>中部>東北”的格局(圖3)。其中各地區(qū)碳排放總量占全國比重呈現(xiàn)差異性變化特征：東部地區(qū)下降(66.38%降至54.31%)、西部地區(qū)上升(21.36%升至28.63%)、中部地區(qū)上升(6.78%升至11.13%)、東北地區(qū)變化不大(5.48%升至5.94%)?？梢?東部仍是中國航空碳排放的主要發(fā)生區(qū)域。

根據(jù)2004—2019年中國航空碳排放總量計(jì)算結(jié)果,使用區(qū)域碳排放指數(shù)模型進(jìn)一步評(píng)價(jià)經(jīng)濟(jì)、人口等條件對中國31個(gè)省域的航空碳排放總量的影響并進(jìn)行等級(jí)分區(qū),借助ArcGIS 10.6軟件,選取了2004、2007、2010、2013、2016和2019年的數(shù)據(jù)繪制中國航空碳排放區(qū)域等級(jí)分布圖(圖4)。結(jié)果顯示：2004—2019年,高排放、中排放等級(jí)區(qū)數(shù)量持續(xù)增加,高排放等級(jí)區(qū)數(shù)量由2個(gè)增至8個(gè),中排放等級(jí)區(qū)數(shù)量由3個(gè)增至6個(gè),在東部地區(qū)、少數(shù)中部地區(qū)聚集,且越到研究后期,分布越廣泛。低排放等級(jí)區(qū)由26個(gè)降至16個(gè),主要集中于西部、東北、中部地區(qū)。據(jù)《中國統(tǒng)計(jì)年鑒》數(shù)據(jù)統(tǒng)計(jì)可得,2004—2019年東部地區(qū)經(jīng)濟(jì)、人口總量占全國的平均比重分別為53.04%、37.59%,2019年廣東、江蘇、山東、浙江、福建、上海的經(jīng)濟(jì)規(guī)模位列全國經(jīng)濟(jì)前十位,北京、天津作為直轄市,經(jīng)濟(jì)總量不容小覷；西部地區(qū)經(jīng)濟(jì)、人口總量占全國的平均比重分別為19.01%、27.41%,陜西省經(jīng)濟(jì)規(guī)模位列全國14位,云南、重慶的經(jīng)濟(jì)發(fā)展速度異軍突起,位居全國前5位；中部地區(qū)經(jīng)濟(jì)、人口總量占全國的平均比重分別為19.99%、26.89%,研究后期河南、湖北、湖南省躋身全國經(jīng)濟(jì)前十位；東北地區(qū)經(jīng)濟(jì)、人口總量占全國的平均比重分別為7.98%、8.11%,遼寧省近年來經(jīng)濟(jì)規(guī)模位列全國15位。經(jīng)濟(jì)與人口規(guī)模大,相應(yīng)的則民航運(yùn)輸規(guī)模也旗鼓相當(dāng),如2010—2019年,東部、西部、中部、東北四大經(jīng)濟(jì)區(qū)機(jī)場旅客吞吐量占全國的比重分別為：55.61%、28.09%、10.16%、6.17%,機(jī)場貨郵吞吐量占全國的比重分別為：76.02%、14.84%、5.63%、3.50%,可見,區(qū)域經(jīng)濟(jì)、人口發(fā)展情況與區(qū)域航空碳排放總量總體呈正向相關(guān),說明其對區(qū)域航空碳排放總量有突出影響。

圖4 2004—2019年中國航空碳排放區(qū)域等級(jí)分布圖Fig.4 Distribution of aviation carbon emission levels in China from 2004 to 2019

2.2.2空間關(guān)聯(lián)分析

使用Moran′s I指數(shù)分析中國航空碳排放總量的省域差異化變化情況,以識(shí)別中國航空碳排放集聚分布及演化特征。具體為,2004—2019年中國航空碳排放量全局Moran′s I指數(shù)航空碳排放總體保持正向空間相關(guān)性,空間響應(yīng)指數(shù)由0.014提高到0.029,增幅達(dá)115.56%,表明中國航空碳排放空間集聚狀態(tài)越來越突出。進(jìn)一步,使用ArcGIS 10.6軟件計(jì)算出2004等6個(gè)年份的航空碳排放空間響應(yīng)指數(shù)的局部Moran′s I值,對各個(gè)省域航空碳排放空間響應(yīng)的空間關(guān)聯(lián)關(guān)系判定(圖5)。結(jié)果顯示：①低低集聚,指自身與周邊航空碳排放空間響應(yīng)程度均較低且空間差異較小的省域,二者呈正相關(guān)關(guān)聯(lián)特征。研究期內(nèi)該類區(qū)域由4個(gè)增至7個(gè),主要呈現(xiàn)向西部→西北→東北逐步擴(kuò)大的集聚趨勢,形成了航空碳排放空間弱響應(yīng)的擴(kuò)展區(qū)域,與圖4的研究結(jié)果一致。②低高集聚,是指自身航空碳排放空間響應(yīng)程度較低,而周邊較高,且空間差異較大的省域,二者呈負(fù)相關(guān)關(guān)聯(lián)特征。研究期內(nèi)該類省域主要在上海周邊分布,因?yàn)樯虾］^其周邊省域航空碳排放差異最高,因此在其周圍形成了低高聚集區(qū)。研究期內(nèi)無高高集聚區(qū)、高低集聚區(qū),說明中國高碳排放區(qū)較分散,還未形成省域之間連片的碳排放高值區(qū)。

圖5 2004—2019年中國局部航空碳排放集聚及演化特征圖Fig.5 Agglomeration and evolution of carbon emissions from local aviation in China from 2004 to 2019

3 航空碳排放效率時(shí)空格局演變分析

3.1 航空碳排放效率時(shí)序演變特征分析

2004—2019年,中國航空碳排放效率總體呈現(xiàn)“M”型波動(dòng)上升趨勢(圖6),經(jīng)歷了由低效率向中效率級(jí)別的轉(zhuǎn)變。2004—2011年,中國航空碳排放效率處于[0.3,0.5]區(qū)間,整體處于低效率級(jí)別。其中“M”的第一個(gè)高值出現(xiàn)在2008年,全球金融危機(jī)導(dǎo)致航空運(yùn)力的下降反倒提升了航空碳排放效率,說明航空碳排放投入的降低減少了碳排放,使得投入和產(chǎn)出之間的比例更協(xié)調(diào),提高了航空碳排放效率。2009—2011年處于低谷期,主要因?yàn)橘Y本投入增加但產(chǎn)出持續(xù)增加導(dǎo)致。第二個(gè)高值出現(xiàn)在2012年,主要因?yàn)橥度虢档突卣{(diào)導(dǎo)致。2012—2019年,中國航空碳排放效率集中于[0.6,0.8]區(qū)間,進(jìn)入中效率等級(jí)。到2019年,效率值又有下降,源于投入和產(chǎn)出的雙下調(diào)。

2004—2019年,四大區(qū)域的航空碳排放效率走勢與全國一致。東部地區(qū)的航空碳排放效率在2004年為低效率,之后除2009—2012年外均處于中效率水平；西部地區(qū)航空碳排放效率經(jīng)歷了低效率(2004—2011年)—高效率(2012—2017年)—低效率(2018—2019年)的轉(zhuǎn)變；中部地區(qū)航空碳排放效率波動(dòng)變化較大,除2017—2019年處于中效率水平,其他年份均處于極低效率、低效率水平；東北地區(qū)航空碳排放效率除2018—2019年處于中效率水平外,其他年份均處于極低效率、低效率水平。研究期內(nèi)雖然部分年份西部地區(qū)航空碳排放效率高于東部地區(qū),整體上呈現(xiàn)東部>西部>中部>東北地區(qū)的態(tài)勢,與上文的航空碳排放的區(qū)域排名結(jié)果一致。主要是因?yàn)闁|部地區(qū)經(jīng)濟(jì)發(fā)達(dá)、人口密集,航空運(yùn)輸總量大但是規(guī)模效應(yīng)突出,且很多城市的作為全國核心的航空樞紐,承載了中國民航運(yùn)輸?shù)闹匾蝿?wù),如北京首都國際機(jī)場是中國最重要、運(yùn)輸生產(chǎn)最繁忙的大型國際航空港,自2008年開始年旅客吞吐量就位居全球第2位；廣州白云國際機(jī)場、上海虹橋國際機(jī)場等重要機(jī)場也均位于東部地區(qū),因此東部地區(qū)民航運(yùn)輸產(chǎn)業(yè)投入產(chǎn)出比高,航空碳排放效率高。西部地區(qū)近年來受到西部大開發(fā)戰(zhàn)略、一帶一路戰(zhàn)略影響,多數(shù)城市經(jīng)濟(jì)發(fā)展速度高于全國平均水平,且部分城市如烏魯木齊、西安、成都的機(jī)場本身就是中國西部地區(qū)重要的航空樞紐,民航規(guī)模位居中國前列。中部地區(qū)、東北地區(qū)因缺少樞紐機(jī)場的中轉(zhuǎn)任務(wù),因此民航運(yùn)輸總規(guī)模相對較弱,航空碳排放效率較低,還需要進(jìn)一步協(xié)調(diào)投入與產(chǎn)出比例。

圖6 2004—2019年中國航空碳排放效率變化趨勢Fig.6 The variation trend of China′s aviation carbon emission efficiency from 2004 to 2019

3.2 航空碳排放效率空間格局演化特征分析

3.2.1空間分布特征分析

根據(jù)2004—2019年中國航空碳排放效率計(jì)算結(jié)果,繪制中國航空碳排放效率空間分布圖(圖7)。具體看：2004 年,僅有上海、海南、云南3地航空碳排放效率為高等級(jí)水平,北京、湖南2地為中等級(jí)水平,遼寧、福建、廣東等11地為低等級(jí)水平,黑龍江、吉林、甘肅等15地為極低等級(jí)水平,東部地區(qū)整體高于其他區(qū)域。2007年,僅有黑龍江、吉林、河北、內(nèi)蒙古、安徽、江西6地為極低等級(jí)水平；西部地區(qū)效率提升較快,如甘肅、貴州、青海、陜西、新疆由極低級(jí)向低級(jí)演化；東部地區(qū)廣東、浙江2地由低級(jí)向中級(jí)轉(zhuǎn)變。2010年,陜西由中級(jí)向高級(jí)演變,其他如山西、海南、北京等地航空碳排放效率等級(jí)有短期下滑。2013年,僅河北保持了極低效率水平,其他各地效率提升較為明顯,中部地區(qū)6省域全部進(jìn)階低級(jí)效率水平,西部地區(qū)7省域進(jìn)入中級(jí)效率水平,平均效率(0.71)超過東部地區(qū)(0.66)成為全國效率最高的區(qū)域。2016年,北京、上海、浙江、甘肅、內(nèi)蒙古、陜西、西藏為航空碳排放高等級(jí)區(qū)域,江蘇、新疆、湖北等14地為中等級(jí)區(qū)域,剩余吉林等10地為低等級(jí)區(qū)域,無極低等級(jí)區(qū)域；西部地區(qū)效率仍然最高,區(qū)域間效率差異減小。2019年,全國碳排放效率有不同程度下降,其中西藏降低為極低等級(jí)區(qū)域,低效率區(qū)域有山東等13地,中等級(jí)區(qū)域北京等12地,高等級(jí)區(qū)域僅有上海等5地；東部地區(qū)效率仍保持較高水平,東北地區(qū)效率水平高于西部和中部地區(qū)。

2004—2019年,各省域航空碳排放效率均存在不同程度的提高：①極低效率區(qū)由15個(gè)降至1個(gè),其集聚區(qū)由中國北部地區(qū)向西部地區(qū)轉(zhuǎn)移。②低效率區(qū)由11個(gè)增至13個(gè),中國絕大多數(shù)地區(qū)屬于該類型區(qū),其集聚區(qū)由西北部向東南區(qū)發(fā)展。③中效率區(qū)數(shù)量大幅增長,由2個(gè)增至12個(gè),其集聚區(qū)由中國東南部地區(qū)向中部和北部拓展。④高效率區(qū)數(shù)量呈現(xiàn)先增加后下降的趨勢,2004—2016年由3個(gè)增長到7個(gè),2016—2019年數(shù)量下降為5個(gè),大部分由內(nèi)陸地區(qū)轉(zhuǎn)移到沿海地區(qū),主要為上海、浙江等地。因此,中國省域航空碳排放效率仍以低效率、中效率為主要類型,未來還有待進(jìn)一步提高中國整體航空碳排放效率。

圖7 2004—2019年中國航空碳排放效率區(qū)域差異分布圖Fig.7 Regional distribution map of China′s aviation carbon emission efficiency from 2004 to 2019

3.2.2空間格局演化特征分析

以2004—2019年31個(gè)省域碳排放效率值為權(quán)重,應(yīng)用ArcGIS 10.6軟件繪制中國航空碳排放效率的標(biāo)準(zhǔn)差橢圓圖,以描述中國航空碳排放效率的空間分布格局(圖8),并基于計(jì)算所得的中心點(diǎn)移動(dòng)方向、橢圓展布性(長軸、短軸與扁率)、橢圓方向性(移動(dòng)方向和距離)等參數(shù)變化情況分析其演化特征(圖9)。

圖8 2004—2019年中國航空碳排放效率分布方向性與演化Fig.8 The direction distribution and evolution of Aviation carbon emission efficiency from 2004 to 2019

(1)空間分布中心變化。從橢圓中心點(diǎn)位置看,2004—2019年中國航空碳排放效率的中心集中分布于中國幾何中心(103.30°E,36.64°N)[36]的東、南方向,說明中國東部、南部地區(qū)的航空碳排放效率較高。其中,2004—2013年中心點(diǎn)自東南向西北方向偏移,此階段西部地區(qū)省域的經(jīng)濟(jì)發(fā)展速度逐步增加,航空運(yùn)輸需求增長,促使航空碳排放效率的提高。2014—2019年中國航空碳排放效率中心不斷向東部、北部遷移,表明中國東部和北部的航空碳排放效率在這一時(shí)期增長較為明顯。但研究期內(nèi)中國整體航空碳排放區(qū)域發(fā)展不平衡性突出,其空間格局表現(xiàn)出由東北—西南的方向分布,從橢圓的扁率看,以2013年為節(jié)點(diǎn),中國航空碳排放效率空間分布的方向性經(jīng)歷了弱化到明顯的過程,表明各省域之間的航空碳排放效率差異性由逐漸減小到逐步增大。

(2)空間分布范圍變化。中國航空碳排放效率標(biāo)準(zhǔn)差橢圓分布范圍在2004—2019年呈先增長后小幅下降增長趨勢,標(biāo)準(zhǔn)差橢圓平均面積為6.02×106km2。其中,2004—2013年標(biāo)準(zhǔn)差橢圓面積持續(xù)增長,由7.11×105km2增至7.48×106km2,長軸與短軸分別由2004年的1043.36 km、1216.85 km變化至2013年的1307.51 km、1067.91 km,表明此階段中國航空碳排放效率的空間分布較為分散,期間航空碳排放效率增高的省域增多；長軸變長、短軸變短,中國航空碳排放效率空間分布格局雖以南—北方向?yàn)橹鲗?dǎo),但東—西方向的發(fā)展也愈加明顯。2013—2019年,標(biāo)準(zhǔn)差橢圓面積由7.48×106km2降至6.88×106km2,長軸與短軸分別變化至2019年的1064.61 km、1124.38 km,說明該階段中國航空碳排放效率的空間分布格局變化不大、較為穩(wěn)定,其效率的空間溢出效應(yīng)尚不明顯,空間分布格局仍以南—北方向?yàn)橹鲗?dǎo),東—西方向的影響較弱。

(3)空間旋轉(zhuǎn)角變化。研究期內(nèi)旋轉(zhuǎn)角θ在43°～80°之間變化,其中2004—2013年橢圓的方位角θ由43°逐漸增加至80°,反映中國航空碳排放效率空間分布格局在東—西方向分布的趨勢逐漸明顯,主要因?yàn)楦魇∮蛱寂欧判什町愔鸩娇s小,區(qū)域分布的均衡性逐步凸顯；2013年之后橢圓的方位角θ由80°降至44°,反映中國航空碳排放效率空間分布格局又逐步回到南—北方向分布為主的類型,且空間分布格局基本保持穩(wěn)定,主要由于研究后期東部、東北地區(qū)航空碳排放效率提高拉大了與東—西方向格局的差距。

圖9 中國航空碳排放效率標(biāo)準(zhǔn)差橢圓核心參數(shù)演化圖Fig.9 Evolution diagram of standard deviation ellipse core parameters of China aviation carbon emission efficiency

4 結(jié)論與討論

4.1 結(jié)論

本文對2004—2019年中國航空碳排放的總量的時(shí)空演化特征進(jìn)行了分析,并基此對2004—2019年航空碳排放效率時(shí)空格局演變過程進(jìn)行了研究。結(jié)果顯示：

(1)中國航空碳排放總量持續(xù)增長,主要受到航空客運(yùn)量增長的影響；同期,中國航空碳排放增長率和航空碳排放強(qiáng)度均呈波動(dòng)下降趨勢,說明中國航空業(yè)節(jié)能減排效果明顯。一方面是因?yàn)檠芯科趦?nèi)航空燃油單位能耗逐年降低,另外一方面,反映了中國為應(yīng)對減排壓力制定的多項(xiàng)措施均對抑制碳排放過快增長發(fā)揮了作用。

(2)航空碳排放總量在中國四大經(jīng)濟(jì)區(qū)的分布總體表現(xiàn)為“東部>西部>中部>東北”的格局,各地區(qū)碳排放總量占全國比重呈現(xiàn)差異性變化特征,東部仍是中國航空碳排放的主要區(qū)域。中國航空碳排放總體保持正向空間相關(guān)性,主要呈現(xiàn)低低集聚、低高集聚狀態(tài),中國高碳排放區(qū)較分散,還未形成省域之間連片的碳排放高值區(qū)。

(3)中國航空碳排放效率總體呈現(xiàn)“M”型波動(dòng)上升趨勢,四大區(qū)域的航空碳排放效率走勢與全國一致,也呈現(xiàn)東部>西部>中部>東北地區(qū)的態(tài)勢,主要是因?yàn)闁|部地區(qū)經(jīng)濟(jì)發(fā)達(dá)、人口密集,航空航空運(yùn)輸總量大且規(guī)模效應(yīng)突出,西部地區(qū)近年來受到西部大開發(fā)戰(zhàn)略、一帶一路戰(zhàn)略影響,多數(shù)城市經(jīng)濟(jì)發(fā)展速度高于全國平均水平,部分樞紐城市民航規(guī)模位居中國前列,帶動(dòng)了航空碳排放效率提高。

(4)各省域航空碳排放效率均存在不同程度的提高,極低效率區(qū)、高效率區(qū)數(shù)量不同程度減少,但低效率區(qū)、中效率區(qū)仍是主要類型。中國東部、南部地區(qū)的航空碳排放效率較高,各省域之間的航空碳排放效率差異性由逐漸減小到逐步增大。中國航空碳排放效率空間分布格局以南—北方向?yàn)橹鲗?dǎo)。

4.2 討論

2004—2019年,中國31個(gè)省域中,經(jīng)濟(jì)發(fā)達(dá)區(qū)域的航空碳排放還處于不斷增長趨勢中,且航空重要樞紐的碳排放呈現(xiàn)點(diǎn)狀快速增長,各省域航空碳排效率以低、中水平為主。因此,提高航空碳排放效率是當(dāng)今促進(jìn)航空碳減排的最主要措施。根據(jù)未來中國航空碳排放空間分布格局演變趨勢,及前文關(guān)于航空碳排放效率成因的探討,建議從以下幾方面提升中國各區(qū)域航空碳排放效率：①對于重要區(qū)域航空樞紐,則需要突出其地域輻射功能,充分帶動(dòng)區(qū)域航空發(fā)展的同時(shí),保障機(jī)場終端區(qū)空域的航空流疏解與暢通,提高航班正常率,并減少機(jī)場擁堵的區(qū)域傳導(dǎo),提高空域利用效率以減少航班繞行與復(fù)飛帶來的碳增加。②各省域機(jī)場在積極提高民航運(yùn)輸服務(wù)能力的同時(shí),應(yīng)統(tǒng)籌協(xié)調(diào),明確各地機(jī)場特色與分工、區(qū)分各地航空目標(biāo)市場,優(yōu)化航線結(jié)構(gòu)、提高上座率,減少區(qū)域同質(zhì)性競爭以提質(zhì)增效。③四大區(qū)域航空減排措施如下：東部地區(qū)應(yīng)繼續(xù)完善綠色低碳民航運(yùn)輸服務(wù)體系；西部地區(qū)應(yīng)繼續(xù)利用國內(nèi)、國際航空樞紐的優(yōu)勢發(fā)揮規(guī)模效應(yīng)；中部地區(qū)應(yīng)利用自身經(jīng)濟(jì)發(fā)展優(yōu)勢積極更新民航交通運(yùn)輸產(chǎn)品以提高市場競爭力降低替代性；東北地區(qū)應(yīng)加快經(jīng)濟(jì)振興,提振民航人氣,通過提高運(yùn)輸規(guī)模達(dá)到碳排放優(yōu)化。④全國總體發(fā)展應(yīng)該統(tǒng)籌一致、共同進(jìn)步。雖然,各區(qū)域航空低碳減排策略各有不同,但全國總體上應(yīng)通過新技術(shù)分享、人才隊(duì)伍建設(shè)等,促進(jìn)民航綠色產(chǎn)業(yè)累積的迭代優(yōu)勢；積極構(gòu)建科技含量高、資源消耗低、環(huán)境污染少的民航綠色產(chǎn)業(yè)結(jié)構(gòu)和生產(chǎn)方式,實(shí)現(xiàn)民航綠色循環(huán)低碳發(fā)展,形成民航綠色經(jīng)濟(jì)發(fā)展新的增長點(diǎn)。

盡管本文綜合考慮了中國航空碳排放及其效率的各項(xiàng)影響因素,如航空運(yùn)輸總量、每千噸運(yùn)輸量運(yùn)送的距離、每千噸公里油耗、航空燃油消耗量等,對中國航空碳排放及效率的區(qū)域差異特征和時(shí)空演化過程進(jìn)行了深入分析。但由于數(shù)據(jù)獲取難度較大,有些資本投入和勞動(dòng)力投入等部分因素等未能全面考慮。未來隨著更全面的航空統(tǒng)計(jì)數(shù)據(jù)庫的完善和出行調(diào)查的實(shí)施,可以對中國航空碳排放做更加精確的模擬和預(yù)測。進(jìn)一步,隨著新技術(shù)實(shí)施和新能源的充分應(yīng)用,以及更加嚴(yán)苛的航空碳減排政策的實(shí)施,中國航空碳排放將呈現(xiàn)更多的新的趨勢,為全球碳減排的有效推進(jìn)和學(xué)術(shù)研究提供了重要方向。