黃和平，李亞麗，楊斯玲,2＊

（1.江西財(cái)經(jīng)大學(xué)經(jīng)濟(jì)學(xué)院，江西南昌 330013；2.華東交通大學(xué)土木建筑學(xué)院，江西南昌 330013）

引言

人類活動(dòng)已經(jīng)極大地改變了地球環(huán)境[1]，城鎮(zhèn)作為人口流與物質(zhì)流的匯集地[2]，其CO2排放量已占到全球總量的80%以上[3]。隨著城鎮(zhèn)化進(jìn)程的加快，大量的農(nóng)村人口向城鎮(zhèn)轉(zhuǎn)移，居民食物消費(fèi)已經(jīng)對(duì)城鎮(zhèn)生態(tài)系統(tǒng)產(chǎn)生重要影響[4,5]。食物消費(fèi)涉及生產(chǎn)、加工、運(yùn)輸、流通、處理等環(huán)節(jié)，這些過(guò)程均伴隨著能源的投入和溫室氣體的排放[6,7]。有研究顯示，與食物消費(fèi)相關(guān)的溫室氣體的排放量已達(dá)到人類活動(dòng)產(chǎn)生的溫室氣體的25%左右[8]。隨著經(jīng)濟(jì)發(fā)展和居民生活水平的提高，居民的食物消費(fèi)行為和方式發(fā)生了重大變化[9]，相應(yīng)地引起了食物碳排放結(jié)構(gòu)的改變。因此，清晰地了解城鎮(zhèn)居民食物消費(fèi)碳排放量、結(jié)構(gòu)構(gòu)成和變化趨勢(shì)，有利于提高居民節(jié)約食物保護(hù)環(huán)境意識(shí)，也有助于針對(duì)性改善食物消費(fèi)對(duì)生態(tài)環(huán)境的影響。

1 文獻(xiàn)綜述

國(guó)內(nèi)外諸多學(xué)者對(duì)食物消費(fèi)碳排放的研究主要集中在碳排放測(cè)度上，從生產(chǎn)和消費(fèi)兩個(gè)視角進(jìn)行計(jì)算，基于生產(chǎn)角度的較多[10-12]，而基于消費(fèi)角度的較少[13-15]?；谙M(fèi)角度測(cè)算的食物消費(fèi)碳排放可分為直接碳排放和間接碳排放。食物消費(fèi)直接碳排放指消耗食物本身產(chǎn)生的碳排放，主要采用碳折算系數(shù)法計(jì)算，Weber 等[16]運(yùn)用碳折算系數(shù)法測(cè)算了美國(guó)不同食物類型的碳排放，認(rèn)為紅肉消費(fèi)產(chǎn)生的碳排放是蔬菜類消費(fèi)碳排放的將近4 倍，Tilman 等[17]得到每單位肉類消費(fèi)產(chǎn)生的碳排放是谷物類消費(fèi)碳排放的10倍，曹志宏等[18]運(yùn)用碳折算系數(shù)法測(cè)算了不同食物類型的碳排放，認(rèn)為中國(guó)居民食物消費(fèi)以糧食、蔬菜和豬肉為主，三者共占2016 年全國(guó)食物碳排放總量的65.71%。食物間接碳排放指食物生產(chǎn)、加工、運(yùn)輸、流通、儲(chǔ)藏等環(huán)節(jié)消耗能源產(chǎn)生的碳排放，主要基于生命周期角度采用碳折算系數(shù)法和投入產(chǎn)出模型計(jì)算，Pelletier 等[19]基于生命周期角度，測(cè)算了食物生產(chǎn)、加工、運(yùn)輸、流通等消費(fèi)過(guò)程產(chǎn)生的間接碳排放，國(guó)內(nèi)學(xué)者也較多基于生命周期過(guò)程測(cè)算食物消費(fèi)間接碳排放[20,21]；部分學(xué)者運(yùn)用投入產(chǎn)出模型測(cè)算食物消費(fèi)間接碳排放[22-24]?；谏芷诮嵌鹊奶寂欧畔禂?shù)法能夠測(cè)算生產(chǎn)、加工、運(yùn)輸?shù)雀麟A段消耗能源產(chǎn)生的碳排放，計(jì)算簡(jiǎn)便，較為實(shí)用[25]；但在中國(guó)，省級(jí)食物消費(fèi)能源數(shù)據(jù)難以獲取，且碳排放系數(shù)法假設(shè)能源碳排放因子不發(fā)生改變，難以反映技術(shù)進(jìn)步，而投入產(chǎn)出模型不僅可以反映技術(shù)進(jìn)步，還可以測(cè)度部門間碳排放的轉(zhuǎn)移，計(jì)算較為準(zhǔn)確[26]。食物消費(fèi)碳排放在研究尺度上，涉及宏觀、中觀和微觀層次，宏觀層次主要基于全國(guó)數(shù)據(jù)測(cè)算碳排放，如高晶等[27]對(duì)中國(guó)城鄉(xiāng)居民食物消費(fèi)產(chǎn)生的碳排放進(jìn)行對(duì)比測(cè)算，認(rèn)為城鎮(zhèn)居民食物消費(fèi)碳排放總量高于農(nóng)村居民；中觀層次上基于省級(jí)、市級(jí)數(shù)據(jù)對(duì)食物消費(fèi)碳排放進(jìn)行測(cè)算，余廣霞等[28]、謝晟宇等[29]分別對(duì)安徽省、昆明市的食物消費(fèi)碳排放進(jìn)行測(cè)算；微觀層次上主要基于微觀調(diào)研數(shù)據(jù)測(cè)算行業(yè)碳排放，朱強(qiáng)等[30]基于微觀調(diào)研數(shù)據(jù)研究了高校食堂浪費(fèi)及外賣的碳足跡，認(rèn)為高校餐飲食物消費(fèi)碳足跡較高。

綜上所述，國(guó)內(nèi)外學(xué)者運(yùn)用不同的方法從不同尺度對(duì)食物消費(fèi)碳排放進(jìn)行測(cè)算，取得了諸多成果。但仍可以從以下幾點(diǎn)深入：一是目前宏觀和中觀尺度的研究主要基于截面數(shù)據(jù)分析中國(guó)城鄉(xiāng)或某個(gè)?。▍^(qū)、市）的食物消費(fèi)碳排放，而?。▍^(qū)、市）之間橫向和縱向?qū)Ρ确治鲚^為缺乏，因此，難以反映食物消費(fèi)產(chǎn)生的碳排放在時(shí)空上的差異；二是在研究方法上，較多基于生命周期角度運(yùn)用碳折算系數(shù)測(cè)算食物消費(fèi)的間接碳排放，忽略了技術(shù)進(jìn)步，而投入產(chǎn)出模型不僅可以反映經(jīng)濟(jì)發(fā)展和技術(shù)進(jìn)步，且計(jì)算結(jié)果較為準(zhǔn)確，涵蓋了行業(yè)間轉(zhuǎn)移的碳排放；三是在研究?jī)?nèi)容上，較多分析時(shí)間序列上城鄉(xiāng)飲食行為產(chǎn)生的差異，較少考慮食物消費(fèi)碳排放的結(jié)構(gòu)和空間分布特征，而明確碳排放的結(jié)構(gòu)特征，針對(duì)性減少碳排放，對(duì)低碳城市建設(shè)具有重要參考價(jià)值。因此，本文在采用碳折算系數(shù)和投入產(chǎn)出模型計(jì)算中國(guó)城鎮(zhèn)居民1990—2018年食物消費(fèi)碳排放的基礎(chǔ)上，運(yùn)用探索性空間數(shù)據(jù)分析（ESDA）方法，借助Geoda1.4.6 軟件，對(duì)中國(guó)除西藏、香港、澳門、臺(tái)灣外的30 個(gè)?。▍^(qū)、市）城鎮(zhèn)居民食物消費(fèi)碳排放的時(shí)空演變特征進(jìn)行分析。

2 研究方法與數(shù)據(jù)來(lái)源

2.1 中國(guó)城鎮(zhèn)居民食物消費(fèi)碳排放的測(cè)算方法

本文研究對(duì)象為食物消費(fèi)碳排放，由食物直接碳排放、家庭間接碳排放和產(chǎn)業(yè)間接碳排放三部分組成。食物直接碳排放指食物本身所含有的碳排放；家庭間接碳排放指烹飪和儲(chǔ)藏食物消耗能源產(chǎn)生的碳排放；產(chǎn)業(yè)間接碳排放指食物消費(fèi)過(guò)程中生產(chǎn)、加工、流通等不同環(huán)節(jié)和不同行業(yè)間轉(zhuǎn)移產(chǎn)生的碳排放。

2.1.1 人均食物碳排放

人均食物直接碳排放測(cè)算采用綜合碳折算系數(shù)法，計(jì)算見(jiàn)下面公式：

式中，Wit為i地區(qū)第t年的人均食物消費(fèi)直接碳排放；Wcj為第j種食物的直接碳排放量；Rj為第j種食物的綜合碳折算系數(shù)；Kj為第j種食物的料肉轉(zhuǎn)化比；r1為糧食的直接碳排放系數(shù)；rj為第j種食物的直接碳折算系數(shù)；Cfj、Cpj、Ccj分別為第j種食物中所含有的脂肪、蛋白質(zhì)以及碳水化合物量；fj、pj、cj分別為第j種食物脂肪、蛋白質(zhì)和碳水化合物所對(duì)應(yīng)的含碳量，其含量數(shù)據(jù)來(lái)源于食物成分表。根據(jù)智靜等[31]計(jì)算的不同食物直接碳折算系數(shù)和動(dòng)物產(chǎn)品的料肉轉(zhuǎn)化比，得到各種食品的綜合碳折算系數(shù)（表1）。

表1 不同種類食物的綜合碳折算系數(shù)（單位：kgC/kg）

2.1.2 家庭間接碳排放測(cè)算

人均家庭的間接碳排放主要用于食物冷藏和炊事所消耗能源產(chǎn)生的碳排放。由于研究區(qū)大部分地區(qū)缺乏炊事消耗的能源數(shù)據(jù)，根據(jù)各?。▍^(qū)、市）統(tǒng)計(jì)年鑒數(shù)據(jù)可知，城鎮(zhèn)地區(qū)家庭常住人口多為三口之家，因此，炊事耗能的碳排放可借鑒曹志宏等[20]的碳折算系數(shù)法進(jìn)行計(jì)算（表2），記為Eit，表示i地區(qū)第t年食物烹飪階段碳排放。

居民家庭用于冷藏食物所消耗的電能產(chǎn)生的碳排放計(jì)算見(jiàn)下式：

表2 加工單位質(zhì)量食物的能耗折算系數(shù)（單位：tC/t）

式中，Sit為i地區(qū)第t年儲(chǔ)藏食物消耗電力產(chǎn)生的碳排放；根據(jù)冰箱電耗限定值及能量效率等級(jí)中的相關(guān)規(guī)定，每臺(tái)電冰箱的年耗電量為507.37 kW·h，電力的碳消耗折算系數(shù)為0.23 kgC/kg[20]，折合116.70 kgC/kg；Nit為i地區(qū)每百戶家庭擁有的電冰箱數(shù)量；nit為i地區(qū)第t年的平均戶規(guī)模。

2.1.3 產(chǎn)業(yè)間接碳排放測(cè)算

投入產(chǎn)出法根據(jù)投入產(chǎn)出表計(jì)算產(chǎn)品的直接消耗系數(shù)和完全消耗系數(shù)，從而測(cè)算行業(yè)的直接和間接碳排放，主要適用于研究碳排放在區(qū)域間、行業(yè)間的轉(zhuǎn)移，通過(guò)編制投入產(chǎn)出表和建立相應(yīng)的線性代數(shù)方程體系，構(gòu)建一個(gè)模擬現(xiàn)實(shí)的國(guó)民經(jīng)濟(jì)結(jié)構(gòu)和社會(huì)再生產(chǎn)過(guò)程的經(jīng)濟(jì)數(shù)學(xué)模型。測(cè)算人均食物消費(fèi)產(chǎn)業(yè)間接碳排放的投入產(chǎn)出模型如下[23]：

式中，Cit為i地區(qū)第t年產(chǎn)業(yè)間接碳排放；CIit為投入產(chǎn)出表中四大產(chǎn)業(yè)（農(nóng)林牧漁業(yè)、食品加工業(yè)、交通運(yùn)輸業(yè)和郵政倉(cāng)儲(chǔ)業(yè)、住宿餐飲業(yè)和批發(fā)零售業(yè)）對(duì)應(yīng)的能源碳排放強(qiáng)度，即產(chǎn)業(yè)碳排放量αEit（α為標(biāo)準(zhǔn)煤碳排放系數(shù)，Eit為四大產(chǎn)業(yè)消耗的標(biāo)準(zhǔn)煤矩陣）和該產(chǎn)業(yè)總產(chǎn)出Xit的比值；Ait為投入產(chǎn)出表中的直接消耗系數(shù)矩陣；I為對(duì)應(yīng)維度的單位矩陣；(I-Ait)-1是列昂惕夫逆矩陣，表示各個(gè)產(chǎn)業(yè)單位產(chǎn)出所需的所有產(chǎn)品的完全投入；Yit為城鎮(zhèn)居民與食品相關(guān)的最終消費(fèi)支出。

2.1.4 人均食物消費(fèi)碳排放測(cè)算

人均食物消費(fèi)碳排放是由人均食物直接碳排放、人均家庭間接碳排放和人均產(chǎn)業(yè)間接碳排放3 部分組成，計(jì)算見(jiàn)下式：

式中，F(xiàn)it為i地區(qū)第t年人均食物消費(fèi)碳排放。

各省（區(qū)、市）總碳排放計(jì)算見(jiàn)下式：

式中，Tit為i地區(qū)第t年食物消費(fèi)總碳排放（萬(wàn)t），Pit為i地區(qū)第t年城鎮(zhèn)常住人口（萬(wàn)人）。

2.2 空間探索性分析

我國(guó)不同?。▍^(qū)、市）之間的碳排放存在較大的空間流動(dòng)性，導(dǎo)致區(qū)域間存在潛在的空間相關(guān)性，本文運(yùn)用ESDA 對(duì)我國(guó)各省（區(qū)、市）城鎮(zhèn)居民食物消費(fèi)總碳排放的空間相關(guān)性進(jìn)行研究?？臻g自相關(guān)研究某一屬性值在相關(guān)區(qū)域空間上的關(guān)聯(lián)集聚程度，分為描述指定屬性在整個(gè)空間聚集性的全局空間自相關(guān)，和進(jìn)一步分析具體局部單元屬性與周邊相鄰空間集聚性的局域空間自相關(guān)[32]。

2.2.1 全局空間自相關(guān)

全局相關(guān)性檢驗(yàn)應(yīng)用于研究整個(gè)區(qū)域空間單元之間的關(guān)聯(lián)性，常用莫蘭指數(shù)（Moran’sI）表示。全局Moran’sI指數(shù)見(jiàn)下式[33]：

式中，n為研究對(duì)象的個(gè)數(shù)，本文中n=30；Xi與Xj分別表示i、j區(qū)域的觀測(cè)值；Wij為空間權(quán)重矩陣；為觀測(cè)值的平均值。在顯著性水平下，若Moran’sI值為正，表示城鎮(zhèn)居民食物消費(fèi)總碳排放整體呈顯著空間集聚；若Moran’sI值為負(fù)，表示碳排放整體呈顯著空間分異。對(duì)于空間是否存在自相關(guān)性，常采用統(tǒng)計(jì)檢驗(yàn)的方式，即Z檢驗(yàn)，Z值為正且顯著時(shí)，表明存在正的空間自相關(guān)，Z值為負(fù)且顯著時(shí)，表明存在空間分異，Z值為零時(shí)，觀測(cè)值呈獨(dú)立隨機(jī)分布。檢驗(yàn)見(jiàn)下式[32]：

式中，Z（I）為標(biāo)準(zhǔn)化統(tǒng)計(jì)量；E（I）為期望值；VAR（I）為方差。

2.2.2 局部空間自相關(guān)

全局Moran’sI指數(shù)僅從整體分析區(qū)域城鎮(zhèn)居民食物消費(fèi)總碳排放的空間差異，不能有效闡述局部的空間集聚強(qiáng)度，還需要進(jìn)一步進(jìn)行局部空間自相關(guān)分析，常用統(tǒng)計(jì)量Local Moran’sI指數(shù)表示。Local Moran’sI見(jiàn)下式：

式中，Ii為L(zhǎng)ocal Moran’sI指數(shù)，其余指標(biāo)含義同式（10）。

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

本文中城鎮(zhèn)居民食物消費(fèi)的種類主要包括：糧食、蔬菜、植物油、瓜果、豬肉、牛羊肉、禽肉、鮮蛋、水產(chǎn)品、酒、奶及奶制品共11 類。具體指標(biāo)和數(shù)據(jù)來(lái)源見(jiàn)表3。

表3 測(cè)算數(shù)據(jù)指標(biāo)和數(shù)據(jù)來(lái)源

食品加工業(yè)能源消耗根據(jù)食品加工業(yè)產(chǎn)值占工業(yè)產(chǎn)值比例，從工業(yè)能源消耗量中剝離出來(lái)；直接消耗系數(shù)根據(jù)《中國(guó)地區(qū)投入產(chǎn)出表》中的計(jì)算方法（直接消耗系數(shù)=直接投入/總投入）測(cè)算，直接投入和總投入數(shù)據(jù)來(lái)源于《中國(guó)地區(qū)投入產(chǎn)出表》；1990—2002 年、2003—2007 年、2008—2018 年分別運(yùn)用2002 年、2007 年、2012 年《中國(guó)地區(qū)投入產(chǎn)出表》計(jì)算的直接消耗系數(shù)，其中，2002 年、2007 年《中國(guó)地區(qū)投入產(chǎn)出表》中交通運(yùn)輸業(yè)及倉(cāng)儲(chǔ)業(yè)和郵政業(yè)部門分離，批發(fā)和零售貿(mào)易業(yè)、住宿和餐飲業(yè)分離，在數(shù)據(jù)處理中將兩者加總合并為交通運(yùn)輸和郵政倉(cāng)儲(chǔ)業(yè)、住宿餐飲和批發(fā)零售業(yè)，2012 年將批發(fā)和零售、住宿和餐飲部門數(shù)據(jù)合并為住宿餐飲和批發(fā)零售業(yè)。個(gè)別缺省數(shù)據(jù)采用數(shù)據(jù)平滑處理的方法進(jìn)行補(bǔ)充。與食物消費(fèi)相關(guān)的支出以1990 年為不變價(jià)，運(yùn)用居民消費(fèi)價(jià)格平減指數(shù)法消除價(jià)格影響。

3 實(shí)證分析

3.1 城鎮(zhèn)居民食物消費(fèi)碳排放的時(shí)間演變分析

受于篇幅限制，本文選取1990 年、2000 年、2010 年和2018 年對(duì)30 個(gè)?。▍^(qū)、市）的人均食物、人均家庭、人均產(chǎn)業(yè)、人均碳排放和總碳排放的時(shí)空演變趨勢(shì)進(jìn)行分析（表4 和表5）。

從表4 可以看出，各?。▍^(qū)、市）的人均間接碳排放在研究期內(nèi)呈增長(zhǎng)趨勢(shì)，其中人均家庭間接碳排放增長(zhǎng)幅度居中，人均產(chǎn)業(yè)間接碳排放增長(zhǎng)幅度最大且始終處于主導(dǎo)地位；大部分地區(qū)的人均直接碳排放在研究期內(nèi)呈增長(zhǎng)趨勢(shì)，但增長(zhǎng)較為緩慢，北京、江蘇、廣東增長(zhǎng)幅度較大，而上海、浙江、河南、陜西、甘肅、青海、寧夏地區(qū)呈下降趨勢(shì)。北京、江蘇、廣東地區(qū)人均直接碳排放增長(zhǎng)幅度較大的原因可從食物消費(fèi)量的原始數(shù)據(jù)探知，北京、江蘇、廣東地區(qū)糧食消費(fèi)數(shù)量有所下降，但蔬菜、鮮瓜果和動(dòng)物性食物消費(fèi)數(shù)量顯著增長(zhǎng)，總體上人均直接碳排放顯著增多，這些地區(qū)經(jīng)濟(jì)發(fā)達(dá)，居民對(duì)肉類、瓜果、蔬菜類的需求增多，反映了該地區(qū)飲食結(jié)構(gòu)的改變；上海、浙江、河南、陜西、甘肅、青海、寧夏地區(qū)糧食、蔬菜消費(fèi)量顯著下降，動(dòng)物性食物的消費(fèi)增長(zhǎng)不顯著，只有瓜果、鮮奶類消費(fèi)增長(zhǎng)幅度較大，總體上人均直接碳排放呈下降趨勢(shì)，但產(chǎn)生機(jī)理不同，上海、浙江地區(qū)經(jīng)濟(jì)發(fā)達(dá)，外來(lái)務(wù)工人員較多，居民收入偏高，外出就餐較多；而河南、陜西、甘肅、青海、寧夏處于內(nèi)陸地區(qū)，水資源相對(duì)缺乏，水產(chǎn)品消耗偏少，城鎮(zhèn)居民體力消耗量有所下降，因此，糧食、蔬菜消費(fèi)量下降，且由于地區(qū)經(jīng)濟(jì)發(fā)展較為落后，即使隨著經(jīng)濟(jì)發(fā)展，居民收入增多，但隨著其他消費(fèi)的增長(zhǎng)和通貨膨脹，居民在動(dòng)物性食物的消費(fèi)上仍然較為節(jié)儉。

人均家庭間接碳排放是由烹飪和儲(chǔ)藏食物消耗能源導(dǎo)致的，而烹飪消耗的能源與食物消耗數(shù)量有關(guān)，依據(jù)上文分析，人均食物消費(fèi)數(shù)量增長(zhǎng)幅度較小，因此，炊事能源產(chǎn)生的碳排放增長(zhǎng)也相對(duì)不顯著；人均家庭間接碳排放的增長(zhǎng)主要是由儲(chǔ)藏食物消耗電能導(dǎo)致的，經(jīng)濟(jì)發(fā)展使得居民購(gòu)買冰箱（冰柜）的數(shù)量顯著增加，而隨著居民生活觀念的轉(zhuǎn)變，父母與子女分居現(xiàn)象增多，居民家庭平均戶規(guī)模呈下降趨勢(shì)，根據(jù)式（5）計(jì)算的人均家庭間接碳排放呈增長(zhǎng)趨勢(shì)。

經(jīng)過(guò)計(jì)算，30 個(gè)?。▍^(qū)、市）的能源利用強(qiáng)度（能源消耗量/產(chǎn)業(yè)產(chǎn)值）皆呈下降趨勢(shì)，反映出經(jīng)濟(jì)發(fā)展帶動(dòng)了技術(shù)進(jìn)步。但是，大部分地區(qū)的人均產(chǎn)業(yè)碳排放呈幾何倍數(shù)增長(zhǎng)，北京、上海、寧夏地區(qū)的人均產(chǎn)業(yè)碳排放在大多年份始終處于前列。其中北京、上海地區(qū)經(jīng)濟(jì)發(fā)達(dá)，第三產(chǎn)業(yè)占比最高，與食物消費(fèi)相關(guān)的基礎(chǔ)設(shè)施投入較為完善，四大產(chǎn)業(yè)產(chǎn)值始終處于前列，且居民在食物消費(fèi)上的支出也處于前列；而寧夏地區(qū)人均產(chǎn)業(yè)碳排放偏高的可能原因是當(dāng)?shù)丶夹g(shù)水平相對(duì)落后，能源利用不充分，第三產(chǎn)業(yè)發(fā)展較為薄弱，與食物消費(fèi)相關(guān)的基礎(chǔ)設(shè)施投入比例相對(duì)偏低，產(chǎn)出相對(duì)偏低。大部分地區(qū)人均產(chǎn)業(yè)碳排放的劇烈增長(zhǎng)期是在2010 年以后，主要是因?yàn)橹袊?guó)經(jīng)濟(jì)處于中高速發(fā)展期，對(duì)基礎(chǔ)設(shè)施建設(shè)的重視、交通運(yùn)輸?shù)谋憷?、互?lián)網(wǎng)經(jīng)濟(jì)的快速發(fā)展等，使得居民與食物相關(guān)的支出增多，人均產(chǎn)業(yè)間接碳排放增多。海南地區(qū)的人均產(chǎn)業(yè)間接碳排放始終處于全國(guó)較低水平，海南地區(qū)雖然四大產(chǎn)業(yè)產(chǎn)值相對(duì)偏低，但能源消耗量相對(duì)偏低，且第三產(chǎn)業(yè)較為發(fā)達(dá)，產(chǎn)出效率高。安徽、福建、江西、山東、河南、湖北、湖南、廣西、云南、陜西等地區(qū)的人均產(chǎn)業(yè)間接碳排放大多數(shù)年份處于較低水平，這與當(dāng)?shù)亟?jīng)濟(jì)發(fā)展水平偏低，居民食物相關(guān)支出也相對(duì)偏低有關(guān)。

從表5 可以看出，中國(guó)30 個(gè)省（區(qū)、市）人均碳排放和總碳排放皆呈顯著增長(zhǎng)趨勢(shì)，人均碳排放的增長(zhǎng)幅度寧夏＞新疆＞上海＞浙江＞青海＞北京，增長(zhǎng)幅度最小的是甘肅地區(qū)。寧夏、新疆增長(zhǎng)幅度大主要是人均產(chǎn)業(yè)碳排放顯著增加導(dǎo)致的，兩地的能源消耗強(qiáng)度雖呈下降趨勢(shì)，但由于其原本的能源消耗水平就都處于全國(guó)較高水平，且兩地的第三產(chǎn)業(yè)發(fā)展較為落后，與食物消費(fèi)相關(guān)的基礎(chǔ)設(shè)施投入較低，產(chǎn)出較低，技術(shù)水平較為落后。因此，寧夏、新疆、青海等經(jīng)濟(jì)發(fā)展較為落后、地廣人稀地區(qū)的首要任務(wù)是加快人口城鎮(zhèn)化建設(shè)，在集中建設(shè)、節(jié)約資源的同時(shí)，促進(jìn)當(dāng)?shù)亟?jīng)濟(jì)發(fā)展。上海、浙江、北京地區(qū)經(jīng)濟(jì)發(fā)達(dá)，技術(shù)較為先進(jìn)，能源利用強(qiáng)度較低，人均碳排放增長(zhǎng)幅度卻較高可能與這些地區(qū)居民在與食物消費(fèi)相關(guān)支出過(guò)高有關(guān)。甘肅地區(qū)人均碳排放增長(zhǎng)幅度最小，與人均產(chǎn)業(yè)碳排放增長(zhǎng)最小密切相關(guān)，甘肅地區(qū)在1990年能源消耗強(qiáng)度偏高，隨著經(jīng)濟(jì)發(fā)展和技術(shù)進(jìn)步，能源利用強(qiáng)度呈下降趨勢(shì)，與食品消費(fèi)相關(guān)的支出雖然呈增長(zhǎng)趨勢(shì)，但與其他地區(qū)相比，消費(fèi)支出偏低，直接能源消耗系數(shù)處于較高水平且變化不顯著，總體來(lái)說(shuō)，人均產(chǎn)業(yè)碳排放變化不顯著?？偺寂欧旁鲩L(zhǎng)幅度廣東＞浙江＞山東＞江蘇＞河北＞上海，增長(zhǎng)幅度最小的是吉林地區(qū)，總碳排放與人均碳排放變化趨勢(shì)不一致的主要原因是城鎮(zhèn)常住人口的變動(dòng)，廣東、浙江、山東、江蘇、上海等經(jīng)濟(jì)發(fā)達(dá)地區(qū)更能吸引人口的流入，人口集聚能進(jìn)一步促進(jìn)經(jīng)濟(jì)增長(zhǎng)，因此，經(jīng)濟(jì)發(fā)達(dá)地區(qū)城鎮(zhèn)常住人口增長(zhǎng)較快，總碳排放變化較多。吉林、甘肅、青海地區(qū)總碳排放變化較小，除與人均碳排放變化較小有關(guān)之外，還與當(dāng)?shù)亟?jīng)濟(jì)發(fā)展較為落后，無(wú)法吸引人流的集聚，且人口流失較為嚴(yán)重密切相關(guān)。因此，這些地區(qū)首要任務(wù)是提高當(dāng)?shù)亟?jīng)濟(jì)發(fā)展水平，政府應(yīng)出臺(tái)政策吸引外來(lái)投資，同時(shí)鼓勵(lì)校企合作，解決人口就業(yè)、社會(huì)保障等問(wèn)題，留存人才。

表4 1990—2018年中國(guó)城鎮(zhèn)居民食物消費(fèi)的人均直接、家庭間接、產(chǎn)業(yè)間接碳排放（單位：kg/人）

表5 1990—2018年中國(guó)城鎮(zhèn)居民食物消費(fèi)的人均碳排放和總碳排放

3.2 城鎮(zhèn)居民食物消費(fèi)碳排放的空間格局演變分析

3.2.1 總體空間格局演變特征

本文基于4 階鄰接矩陣，利用Geoda1.4.6 分析軟件，計(jì)算1990—2018 年中國(guó)省際城鎮(zhèn)居民食物消費(fèi)總碳排放的全局Moran’sI指數(shù)，分析總體格局演化特征（表6）。從表6 可以看出，除2003 年、2004年、2005 年、2006 年、2007 年和2009 年外，其余年份總碳排放Moran’sI指數(shù)均在10%或5%的水平上通過(guò)顯著性檢驗(yàn)，說(shuō)明城鎮(zhèn)居民食物消費(fèi)總碳排放大部分年份呈現(xiàn)正的空間相關(guān)性、高（低）相鄰省際單元相對(duì)集聚模式。從整體上看，食物消費(fèi)總碳排放Moran’sI指數(shù)值呈“M”形趨勢(shì)，并根據(jù)其數(shù)值的波動(dòng)情況，空間集聚趨勢(shì)可大致分為三個(gè)階段。1990—2000 年Moran’sI指數(shù)大體呈上升趨勢(shì)，表明食物消費(fèi)總碳排放的空間集聚效應(yīng)呈增強(qiáng)趨勢(shì)；2001—2009年Moran’sI指數(shù)呈下降趨勢(shì)，且大部分年份未通過(guò)顯著性水平檢驗(yàn)，不存在顯著的空間集聚效應(yīng)，省域食物消費(fèi)總碳排放的Moran’sI指數(shù)既受鄰近省份的影響，也與當(dāng)?shù)鼐用耧嬍诚M(fèi)習(xí)慣、當(dāng)?shù)鼗A(chǔ)設(shè)施便利程度相關(guān)，兩者綜合效應(yīng)導(dǎo)致Moran’sI指數(shù)不顯著。2010—2018 年Moran’sI指數(shù)呈下降趨勢(shì)，但皆通過(guò)了顯著性水平檢驗(yàn)，表明存在正的空間集聚效應(yīng)，但空間關(guān)聯(lián)性減弱。

3.2.2 局部空間格局演化特征

本文進(jìn)一步對(duì)其局部空間集聚格局演化特征進(jìn)行探究。選取1990 年、2000 年、2010 年和2018 年4 個(gè)典型年份對(duì)中國(guó)城鎮(zhèn)居民食物消費(fèi)總碳排放的局部Moran’sI指數(shù)進(jìn)行計(jì)算（表7）。從表7 可以看出，1990 年、2000 年、2010 年和2018 年4 個(gè)年份僅存在高高型和低高型兩種類型。1990 年高高集聚區(qū)主要集中在東北和東部偏北沿海地區(qū)，這些地區(qū)當(dāng)年為資源大?。▍^(qū)、市），能源利用強(qiáng)度偏高，經(jīng)濟(jì)較為發(fā)達(dá)，與食物消費(fèi)相關(guān)的工業(yè)制造、交通運(yùn)輸、餐飲住宿等基礎(chǔ)設(shè)施投入較多，人均食物消費(fèi)碳排放偏高，同時(shí)經(jīng)濟(jì)發(fā)展存在集聚效應(yīng)，經(jīng)濟(jì)發(fā)達(dá)地區(qū)吸引人口集聚，進(jìn)一步提高了碳排放。上海處于低高型區(qū)，碳排放低值被周圍高值包圍，上海食物消費(fèi)碳排放在1990年處于低值的可能原因是能源利用強(qiáng)度偏低，與食物消費(fèi)相關(guān)的基礎(chǔ)設(shè)施投入比例偏低，且當(dāng)年上海人口相對(duì)較少，雖然居民在食物消費(fèi)相關(guān)上的支出偏高，綜合作用導(dǎo)致總體碳排放偏低。2000 年、2010 年和2018 年，隨著經(jīng)濟(jì)的發(fā)展，經(jīng)濟(jì)發(fā)達(dá)地區(qū)在發(fā)生變化，總碳排放高高型區(qū)域一直集中在上海、江蘇、浙江地區(qū)，這些地區(qū)處于東部沿海，經(jīng)濟(jì)發(fā)達(dá)，吸引眾多外來(lái)人口的集聚，人口集聚進(jìn)一步促進(jìn)經(jīng)濟(jì)增長(zhǎng)，雖然隨著技術(shù)進(jìn)步，能源利用強(qiáng)度降低，但居民在食物消費(fèi)上的需求增多，與食物消費(fèi)相關(guān)的基礎(chǔ)設(shè)施投入比例增多，與食物消費(fèi)相關(guān)的支出顯著增長(zhǎng)，因此，這些地區(qū)自2000 年以后始終處于高高區(qū)；而天津在2000 年以后始終處于低高型區(qū)，天津的食物消費(fèi)總碳排放在2000 年以后處于較低水平，除了與技術(shù)進(jìn)步，還與北京的發(fā)展吸引人口的流入有關(guān)。安徽、福建與周圍地區(qū)相比，皆處于經(jīng)濟(jì)相對(duì)不發(fā)達(dá)區(qū)域，周圍地區(qū)人口的集聚造成當(dāng)?shù)厝丝诘牧魇А?/p>

表6 1990—2018年中國(guó)城鎮(zhèn)居民食物消費(fèi)總碳排放全局自相關(guān)Moran’s I 指數(shù)

表7 中國(guó)城鎮(zhèn)居民食物消費(fèi)總碳排放的局部空間相關(guān)性聚類分布

4 研究結(jié)論與政策建議

4.1 研究結(jié)論

本文運(yùn)用碳折算系數(shù)法和投入產(chǎn)出模型分別測(cè)度了1990—2018 中國(guó)城鎮(zhèn)居民食物消費(fèi)的直接碳排放、家庭間接碳排放和產(chǎn)業(yè)間接碳排放，并利用空間探索性分析方法對(duì)總碳排放的空間特征進(jìn)行分析，主要結(jié)論如下：

（1）中國(guó)30 個(gè)?。▍^(qū)、市）的人均間接碳排放在研究期內(nèi)呈增長(zhǎng)趨勢(shì)，人均家庭間接碳排放增長(zhǎng)幅度居中，人均產(chǎn)業(yè)間接碳排放增長(zhǎng)幅度最大且始終處于主導(dǎo)地位；大部分地區(qū)的人均直接碳排放增長(zhǎng)較為緩慢。

（2）各省（區(qū)、市）人均碳排放和總碳排放皆呈顯著增長(zhǎng)趨勢(shì)，人均碳排放的增長(zhǎng)幅度寧夏＞新疆＞上海＞浙江＞青海＞北京，甘肅最??；總碳排放增長(zhǎng)幅度廣東＞浙江＞山東＞江蘇＞河北＞上海，吉林最小。

（3）中國(guó)城鎮(zhèn)居民食物消費(fèi)的總碳排放在大部分年份呈正的空間相關(guān)性，1990—2000 年顯著性呈上升趨勢(shì)，2001—2009 年顯著性下降，且該區(qū)間大部分年份不顯著，2010—2018 年顯著性呈下降趨勢(shì)，整體上呈“M”形波動(dòng)。局部空間相關(guān)性結(jié)果表明，1990年、2000 年、2010 年和2018 年4 個(gè)典型年份的食物消費(fèi)總碳排放存在高高型和低高型兩種，且2000 年及以后高高型穩(wěn)定在上海、江蘇、浙江地區(qū)。

4.2 政策建議

根據(jù)以上研究結(jié)果，本文對(duì)減少城鎮(zhèn)居民食物消費(fèi)碳排放的政策建議如下：

（1）鼓勵(lì)大家庭集聚生活，合理搭配膳食結(jié)構(gòu)。從食物消費(fèi)的直接碳排放結(jié)果可以看出，居民的飲食結(jié)構(gòu)在發(fā)生變化，由以“植物性消費(fèi)”為主向“植物性+動(dòng)物性消費(fèi)”方式轉(zhuǎn)變，在動(dòng)物性食物消費(fèi)上的數(shù)量顯著增加，不僅增加了食物消費(fèi)的碳排放，也會(huì)帶來(lái)肥胖、高血脂等問(wèn)題，損害人體健康。因此，應(yīng)該在生活中葷素結(jié)合，合理搭配膳食結(jié)構(gòu)。大家庭集聚生活不僅可以形成家庭規(guī)模優(yōu)勢(shì)，減少炊事能源消耗和冰箱個(gè)數(shù)，降低家庭間接碳排放，且大家庭集聚生活可以培養(yǎng)感情，相互照顧，在一定程度上可以緩解我國(guó)目前面臨的“少無(wú)所養(yǎng)，老無(wú)所依”的問(wèn)題。

（2）促進(jìn)產(chǎn)業(yè)結(jié)構(gòu)升級(jí)，加快技術(shù)進(jìn)步。從城鎮(zhèn)居民食物消費(fèi)的產(chǎn)業(yè)碳排放分析可以發(fā)現(xiàn)，能源利用強(qiáng)度高的地區(qū)食物消費(fèi)碳排放偏高，因此應(yīng)該提高當(dāng)?shù)丶夹g(shù)水平。政府在加大科研投入力度的同時(shí)，可以鼓勵(lì)外來(lái)企業(yè)在當(dāng)?shù)亟◤S，鼓勵(lì)廠商引入先進(jìn)設(shè)備，提倡校企合作，提高人才待遇等手段促進(jìn)當(dāng)?shù)丶夹g(shù)水平的提高。在結(jié)果分析中可以發(fā)現(xiàn)，寧夏、新疆等經(jīng)濟(jì)發(fā)展較為落后的地區(qū)，居民食物消費(fèi)的人均產(chǎn)業(yè)碳排放也相對(duì)偏高，除因當(dāng)?shù)丶夹g(shù)水平有限外，也與當(dāng)?shù)爻擎?zhèn)規(guī)劃布局相關(guān)建設(shè)投入有關(guān)，更深層次原因是與當(dāng)?shù)禺a(chǎn)業(yè)結(jié)構(gòu)相關(guān)，第二產(chǎn)業(yè)基礎(chǔ)設(shè)施投入比重較大，第三產(chǎn)業(yè)發(fā)展落后。

（3）縮小地區(qū)經(jīng)濟(jì)發(fā)展水平差距，加強(qiáng)區(qū)域經(jīng)濟(jì)合作，促進(jìn)資源流轉(zhuǎn)。根據(jù)食物消費(fèi)總碳排放和空間相關(guān)性分析結(jié)果可以發(fā)現(xiàn)，總碳排放高的地區(qū)一般為經(jīng)濟(jì)較為發(fā)達(dá)的地區(qū)，且高高集聚區(qū)也為經(jīng)濟(jì)發(fā)達(dá)區(qū)域，經(jīng)濟(jì)發(fā)展優(yōu)勢(shì)可以吸引人口流入，帶動(dòng)當(dāng)?shù)亟?jīng)濟(jì)發(fā)展，同時(shí)也具有一定的空間溢出效應(yīng)，對(duì)周圍地區(qū)的經(jīng)濟(jì)發(fā)展產(chǎn)生影響。經(jīng)濟(jì)發(fā)展可以提高技術(shù)水平，目前主要是經(jīng)濟(jì)發(fā)展決定了人口的流動(dòng)方向。因此，經(jīng)濟(jì)落后地區(qū)應(yīng)該結(jié)合自身優(yōu)勢(shì)，加強(qiáng)與周圍優(yōu)勢(shì)地區(qū)的經(jīng)濟(jì)合作，促進(jìn)資源要素如資金、技術(shù)、資源等的空間流轉(zhuǎn)，緩解經(jīng)濟(jì)發(fā)達(dá)地區(qū)的人口壓力。