呂夢姍 李 雷,2,3 王成軍,2,3
(1.浙江農(nóng)林大學經(jīng)濟管理學院;2.浙江農(nóng)林大學鄉(xiāng)村振興研究院;3.浙江農(nóng)林大學碳中和研究院,浙江 杭州 311300)
改革開放以來,中國城鄉(xiāng)居民收入和消費日漸增長,居民的食物消費正從溫飽型消費向發(fā)展型消費過渡[1],食物消費總量顯著提升并呈現(xiàn)多元化趨勢。居民食物消費結(jié)構(gòu)的變化,可能會導(dǎo)致食物消費碳排放的變化,影響資源環(huán)境的可持續(xù)性。中國作為一個負責任的發(fā)展中國家,在2020年提出了“雙碳”目標。作為改善居民健康和協(xié)調(diào)自然環(huán)境系統(tǒng)可持續(xù)性的有效途徑,近年來食物系統(tǒng)受到廣泛關(guān)注。《中國農(nóng)業(yè)產(chǎn)業(yè)發(fā)展報告2022》指出要加速推進肉類消費向健康營養(yǎng)低碳方向轉(zhuǎn)型,促進農(nóng)食系統(tǒng)的碳減排,進而提高中國農(nóng)食系統(tǒng)的韌性。食物系統(tǒng)作為溫室氣體排放的主要來源,約占全球碳排放的25%~30%,而食物消費又與居民日常生活息息相關(guān)[2],食物作為連接城鄉(xiāng)的載體,起著非常重要的作用,因此向大眾提倡低碳綠色的食物消費方式至關(guān)重要[3]。長三角地區(qū)作為我國最發(fā)達的經(jīng)濟圈之一,該地區(qū)居民快速的生活節(jié)奏導(dǎo)致其食物消費模式也在快速變化,這給中國實現(xiàn)“雙碳”目標帶來極大挑戰(zhàn)。長三角地區(qū)有責任探索低碳可持續(xù)的食物消費路徑,率先推進全面綠色低碳轉(zhuǎn)型。并且長三角已經(jīng)開啟探索共同富裕的新模式,該地區(qū)城鄉(xiāng)居民收入差距較中國其他地區(qū)小,研究該地區(qū)城鄉(xiāng)居民的食物消費碳排放變化差異,能為促進城鄉(xiāng)協(xié)調(diào)發(fā)展提供參考,具有較強的現(xiàn)實意義。同時本文的理論意義在于在可持續(xù)發(fā)展理論的基礎(chǔ)上,豐富了農(nóng)食系統(tǒng)維度上的可持續(xù)發(fā)展理論,進一步探討食物消費與碳排放二者的邏輯關(guān)系。本文旨在“雙碳”目標背景下探究長三角城鄉(xiāng)居民的食物消費碳排放差異以及不同類型食物的人均消費量對長三角城鄉(xiāng)居民的食物消費碳排放影響,對明確居民食物消費碳排放變化和實現(xiàn)食物消費可持續(xù)發(fā)展有重要意義。
居民食物消費碳排放,是指居民在日常的食物消費中所產(chǎn)生的二氧化碳量,來源包括食物本身、食物的生產(chǎn)加工運輸、食物的儲存與烹飪等環(huán)節(jié)[4]。國內(nèi)外學者就居民食物消費碳排放,運用生命周期法[5]、投入產(chǎn)出法[6]、碳折算系數(shù)法[7]等方法進行測算研究。目前關(guān)于食物消費碳排放的研究主要集中在以下三個方面:(1)食物消費碳排放的時空演變研究。就中國而言,中國各地食物消費人均碳排放呈現(xiàn)較顯著的民族性、地域空間差異特征[8]。從城鄉(xiāng)分異視角上看,近年來部分學者通過研究發(fā)現(xiàn):城鎮(zhèn)居民在食物消費碳排放總量上高于農(nóng)村居民。一方面,隨著越來越多的農(nóng)村人口進入城市,將促進城鎮(zhèn)居民食物消費碳排放的增加,另一方面,城鄉(xiāng)居民食物消費結(jié)構(gòu)多樣化導(dǎo)致食物消費碳排放較之前年份有所下降[9,10]。總體上看,中國城鄉(xiāng)居民的食物碳消費同時呈現(xiàn)二元性和趨同性特征[8,11,12]。(2)食物消費碳排放的影響因素研究。有部分學者指出,GDP指數(shù)、人口規(guī)模、城鎮(zhèn)化、居民收入等因素會對食物消費碳排放帶來不同程度的影響[13-15]。(3)食物消費結(jié)構(gòu)變動與碳排放關(guān)系。張翠玲等綜合考慮食物消費的經(jīng)濟成本最低、水資源消耗最少、耕地資源消耗最少和溫室氣體排放最少4個目標對當前中國居民食物消費結(jié)構(gòu)進行優(yōu)化[16];韓金雨等通過VAR模型探尋食物消費結(jié)構(gòu)的升級與農(nóng)業(yè)碳排放之間的關(guān)系,發(fā)現(xiàn)兩者存在一定的關(guān)聯(lián),且具有滯后性,同時指出農(nóng)業(yè)碳排放對植物性食物消費的沖擊更為靈敏[17]。然而目前對全國最發(fā)達的經(jīng)濟區(qū)——長三角地區(qū)居民的食物消費碳排放的研究比較少見,這不利于準確評估經(jīng)濟快速發(fā)展背景下的中國城鄉(xiāng)居民食物消費碳排放的變動情況,為此,本文從中國長三角地區(qū)的1995—2019年居民食物消費數(shù)據(jù)入手,分析長三角城鄉(xiāng)居民食物消費碳排放整體情況,并基于無條件分位數(shù)回歸(UQR),研究不同類型食物的人均消費量對長三角城鄉(xiāng)居民食物消費碳排放的影響,為實現(xiàn)中國城鄉(xiāng)居民食物消費可持續(xù)發(fā)展提供參考依據(jù)。
本文設(shè)定的計量模型如下,首先是以居民人均食物消費碳排放作為被解釋變量為例:

(1)
其中i為省份(直轄市),t為年份,α為截距項,j表示不同的食物種類;β1-β12分別表示解釋變量的估計系數(shù),度量各影響因素對長三角城鄉(xiāng)居民食物消費碳排放的影響;μi、θt、εit分別為空間效應(yīng)、時間效應(yīng)和隨機擾動項。式1的被解釋變量是居民人均食物消費碳排放,以其對數(shù)形式表示。主要解釋變量是各類食物的人均消費量,控制變量包括人均可支配收入、人口數(shù)量和恩格爾系數(shù),除恩格爾系數(shù)以外,其他幾項均用對數(shù)形式表示。同時選取老齡化率和人口撫養(yǎng)比(老年撫養(yǎng)比與兒童撫養(yǎng)比之和)作為人口特征變量。最后設(shè)置時間與地區(qū)虛擬變量,同時考慮個體與時間固定效應(yīng)。其次本文進一步分析各類食物消費的碳排放所受影響,由于文章篇幅有限,選取的被解釋變量以3類基礎(chǔ)性食物(糧食類、肉禽類、蔬菜類)的人均消費碳排放為例。變量的設(shè)置與式1相似,構(gòu)建的計量模型如式(2)—(4)所示:
(2)
(3)
(4)
上式所涉及的具體變量的編碼和含義如表1所示。

表1 變量含義說明
無條件分位數(shù)回歸(UQR)是對條件分位數(shù)回歸(CQR)的補充和拓展。UQR主要應(yīng)用于勞動經(jīng)濟學領(lǐng)域,多數(shù)學者用此研究不同影響因素對不同群體的收入分布差距的影響[18]。近年來,越來越多的學者將UQR用于其他領(lǐng)域。方超研究信息技術(shù)對學生教育的公平性[19];周春芳和蘇群分析農(nóng)民工與城鎮(zhèn)職工的就業(yè)質(zhì)量差異[20];鄭宏遠和李谷成研究數(shù)字普惠金融對縣域農(nóng)業(yè)全要素生產(chǎn)率的影響[21]。本文基于Firpo等提供的無條件分位數(shù)回歸,利用再中心化影響函數(shù)(RIF)進行無條件估計[22],參考鄭宏遠和李谷成的做法[21],表示為式(5):
式1和式5中的LNTPEit均表示人均食物消費碳排放的對數(shù)。qτ是人均食物消費碳排放在第τ個分位點的無條件分布,I(LNTPEit≤qτ)是一個區(qū)分人均食物消費碳排放比大還是小的二元變量,fLNTPF(qτ)表示在qτ估計的人均食物消費碳排放的概率密度函數(shù)。由于本文所利用的是省級面板數(shù)據(jù),必須對那些不隨時間變化但可能會影響居民人均食物消費碳排放的因素加以控制。基于此,本文主要利用BORGEN提出的模型進行估計[23]。
本文將食物劃分為7個主要食物類別,選取樣本為1995—2019年長三角地區(qū)各省市城鄉(xiāng)居民食物消費的觀察值。其中食物消費量歷年數(shù)據(jù)取自《中國住戶統(tǒng)計年鑒》《中國統(tǒng)計年鑒》以及各省市統(tǒng)計年鑒。本文的食物消費碳排放系數(shù)參照Poore等的食物系統(tǒng)碳排放數(shù)據(jù)庫的數(shù)據(jù)[24],通過整理如表2所示。根據(jù)歷年食物消費量和食物消費碳排放系數(shù),通過計算得到某一省份的城鎮(zhèn)或農(nóng)村居民在某一年的人均食物消費碳排放總量。除此之外,其余的解釋變量則來自《中國人口與就業(yè)統(tǒng)計年鑒》與各省市統(tǒng)計年鑒。

表2 食物的碳排放系數(shù)
長三角居民的食物消費碳排放總量呈逐步上升的趨勢,從1995年的2137.91億千克二氧化碳當量增至2019年的3419.88億千克二氧化碳當量。圖1為長三角城鄉(xiāng)居民整體食物消費碳排放總量近年來的變化情況。其中,城鎮(zhèn)居民的食物消費碳排放呈穩(wěn)步上升趨勢,農(nóng)村居民的食物消費呈弱N型趨勢,并且城鎮(zhèn)居民的食物消費碳排放在2005年開始超過農(nóng)村居民。表明城鄉(xiāng)二元結(jié)構(gòu)的人口流動極大地促進了城鎮(zhèn)居民食物消費的增加。但由于居民整體生活水平的改善,農(nóng)村居民的食物消費碳排放在2010年前后也緩慢增加。

圖1 1995-2019年長三角城鄉(xiāng)居民食物消費碳排放總量變化情況
長三角居民人均食物消費碳排放計算結(jié)果如圖2所示。整體上看,居民人均食物消費碳排放近25年增長了30.94%。從城鄉(xiāng)差異上看,2010年以前農(nóng)村居民的人均食物消費碳排放大于城鎮(zhèn)居民,2010年以后城鎮(zhèn)居民的人均食物消費碳排放超過了農(nóng)村居民,到2019年農(nóng)村居民食物消費碳排放又反超城鎮(zhèn)居民。但從整體上看,城鎮(zhèn)居民的人均食物消費碳排放與農(nóng)村居民相當,并且二者的差距也在逐漸減小。結(jié)合上文所提到的長三角城鄉(xiāng)居民食物消費碳排放總量變化情況,反映出長三角城鄉(xiāng)居民食物消費碳排放同時具有二元性和趨同性。

圖2 長三角居民1995-2019年人均食物消費碳排放變化趨勢
從長三角地區(qū)居民近25年的主要食物消費的整體上看,居民人均食物消費碳排放為1256.23kg CO2e,總體呈正偏態(tài)分布。從頻數(shù)分布上看,36%的食物消費人均碳排放小于1200 kg CO2e,54%的食物消費人均碳排放在1200—1500 kg CO2e之間,僅有10%的食物消費人均碳排放超過1500 kg CO2e。長三角地區(qū)農(nóng)村居民的人均食物消費碳排放略高于城鎮(zhèn)居民。從不同食物類型上看,除了人均糧食類消費碳排放農(nóng)村居民高于城鎮(zhèn)居民外,在其他食物類型上城鎮(zhèn)居民的食物消費碳排放均高于農(nóng)村居民。并且長三角城鄉(xiāng)居民消費在肉禽類食物消費上產(chǎn)生的碳排放都是最多,城鎮(zhèn)居民消費油脂類食物產(chǎn)生的碳排放最少,農(nóng)村居民消費干鮮瓜果類產(chǎn)生的碳排放最少。與上述情況相同,除了糧食類食物,城鎮(zhèn)居民在其他種類的食物消費量都要大于農(nóng)村居民。城鎮(zhèn)居民的收入水平與人口數(shù)量高于農(nóng)村居民,而在恩格爾系數(shù)、老齡化率和人口撫養(yǎng)比這三方面,農(nóng)村居民要高于城鎮(zhèn)居民。

表3 長三角城鄉(xiāng)居民食物消費碳排放的無條件分位數(shù)回歸結(jié)果

(續(xù)表3)
本文考慮城鄉(xiāng)差異,對樣本數(shù)據(jù)進行回歸,得到以10%、25%、50%、75%、90%為代表的長三角居民食物消費碳排放的無條件分位數(shù)固定效應(yīng)模型回歸結(jié)果,如表3所示,可以看出:各類食物的人均消費量大部分能顯著反映長三角城鄉(xiāng)居民食物消費碳排放變化。對于低分位食物消費碳排放群體,糧食類和肉禽類消費均會導(dǎo)致城鄉(xiāng)居民食物消費碳排放增加,蛋類和干鮮瓜果類消費有助于減緩城鄉(xiāng)居民食物消費碳排放。其中低分位食物消費碳排放群體的城鄉(xiāng)差異主要體現(xiàn)在肉禽類消費上,農(nóng)村居民的肉禽類消費相較于城鎮(zhèn)居民更能導(dǎo)致其食物消費碳排放增加。對于中等分位食物消費碳排放群體,糧食類、肉禽類及水產(chǎn)品消費會導(dǎo)致城鎮(zhèn)居民食物消費碳排放增加,同時糧食類和水產(chǎn)品能減緩農(nóng)村居民的食物消費碳排放。中等分位食物消費碳排放群體的城鄉(xiāng)差異則主要體現(xiàn)在糧食類消費上。對于高分位食物消費碳排放群體,糧食類和水產(chǎn)品均能導(dǎo)致城鄉(xiāng)居民食物消費碳排放增加,肉類、蛋類和油脂類消費一定程度上減緩城鄉(xiāng)居民的食物消費碳排放。高分位食物消費碳排放群體的城鄉(xiāng)差異主要體現(xiàn)在干鮮瓜果類消費上,并且干鮮瓜果類會導(dǎo)致城鎮(zhèn)居民食物消費碳排放增加。
此外,其他控制變量如人均可支配收入、人口數(shù)量、恩格爾系數(shù)、老齡化率、人口撫養(yǎng)比等對居民食物消費碳排放均有顯著的異質(zhì)性影響,考慮到本文重點關(guān)注城鄉(xiāng)居民的不同類型食物消費量與其食物消費碳排放的關(guān)系,在此不做討論。
本文選取長三角城鄉(xiāng)居民的3類基礎(chǔ)性食物(糧食類、肉禽類、蔬菜類)消費碳排放的無條件分位數(shù)回歸結(jié)果進行展示,如表4—6所示。
其中表4是針對居民人均糧食類食物消費碳排放進行的實證分析。相較于農(nóng)村居民,水產(chǎn)品、油脂類、干鮮瓜果類消費增加,會導(dǎo)致城鎮(zhèn)居民的糧食類消費碳排放增加。從無條件分位數(shù)回歸上看,對于低分位糧食類消費碳排放群體,蛋類消費會導(dǎo)致長三角城鄉(xiāng)居民糧食類消費碳排放增加,肉禽類、水產(chǎn)品、油脂類消費對城鄉(xiāng)居民糧食類消費碳排放的影響呈相反關(guān)系。其中肉禽類和油脂類消費更能引起農(nóng)村居民糧食類消費碳排放增長,水產(chǎn)品消費則更能引起城鎮(zhèn)居民糧食類消費碳排放增長。從城鄉(xiāng)差異上看,25分位糧食類消費碳排放群體的城鄉(xiāng)差異較大,主要體現(xiàn)在蛋類消費上,農(nóng)村居民的蛋類消費相較于城鎮(zhèn)居民更能減緩其糧食類消費碳排放。對于中分位糧食類消費碳排放群體,水產(chǎn)品、油脂類、干鮮瓜果類消費對城鄉(xiāng)居民糧食類消費碳排放的影響呈相反關(guān)系。其中干鮮瓜果類消費會導(dǎo)致農(nóng)村居民糧食類消費碳排放增加,水產(chǎn)品和油脂類消費會導(dǎo)致城鎮(zhèn)居民糧食類消費碳排放增加。并且中分位糧食類消費碳排放群體的城鄉(xiāng)差異則主要體現(xiàn)在油脂類消費上。對于高分位糧食類消費碳排放群體,油脂類會導(dǎo)致城鄉(xiāng)居民糧食類消費碳排放均有所增長,肉禽類消費能減緩城鄉(xiāng)居民糧食類消費碳排放。同時,蛋類和干鮮瓜果類對城鄉(xiāng)居民糧食類消費碳排放的影響呈相反關(guān)系,蛋類消費對城鎮(zhèn)居民糧食類消費碳排放有減緩作用,干鮮瓜果類消費對農(nóng)村居民糧食類消費碳排放有減緩作用。并且90分位糧食類消費碳排放群體的城鄉(xiāng)差異也相對較大,主要體現(xiàn)在干鮮瓜果類消費上。綜上,相較于農(nóng)村居民,城鎮(zhèn)居民通過對不同類型的食物消費,能顯著促進其糧食類消費碳排放。

表4 長三角居民糧食類食物消費碳排放的無條件分位數(shù)回歸結(jié)果
表5是針對居民人均肉禽類食物消費碳排放進行的實證分析,結(jié)果表明糧食類消費能顯著減少長三角城鄉(xiāng)居民的肉禽類消費碳排放,水產(chǎn)品消費能顯著引發(fā)長三角城鄉(xiāng)居民的肉禽類消費碳排放增長。從無條件分位數(shù)回歸上看,對于低分位肉禽類消費碳排放群體,糧食類、蛋類、水產(chǎn)品、油脂類消費對城鄉(xiāng)居民肉禽類消費碳排放的影響呈相反關(guān)系。其中糧食類、蛋類、水產(chǎn)品消費會導(dǎo)致農(nóng)村居民肉禽類消費碳排放增加,油脂類消費引起城鎮(zhèn)居民肉禽類消費碳排放增長。并且低分位肉禽類消費碳排放群體的城鄉(xiāng)差異較大,主要體現(xiàn)在糧食類和油脂類消費上。對于中高分位肉禽類消費碳排放群體,糧食類、蛋類、干鮮瓜果類消費對城鄉(xiāng)居民肉禽類消費碳排放的影響呈相反關(guān)系,這三類食物消費均會引起城鎮(zhèn)居民肉禽類消費碳排放增長,而會減緩農(nóng)村居民肉禽類消費碳排放。中高分位的肉禽類消費碳排放群體的城鄉(xiāng)差異較小,主要受到油脂類、蔬菜類、糧食類消費的影響。

表5 長三角城鄉(xiāng)居民肉禽類食物消費碳排放的無條件分位數(shù)回歸結(jié)果
表6是居民人均蔬菜類食物消費碳排放實證分析,結(jié)果表明除了肉禽類和干鮮瓜果類對農(nóng)村居民的蔬菜類消費碳排放有減緩作用外,其他食物消費會導(dǎo)致長三角城鄉(xiāng)居民的蔬菜類消費碳排放增長。從無條件分位數(shù)回歸上看,對于低分位蔬菜類消費碳排放群體,糧食類對城鄉(xiāng)居民蔬菜類消費碳排放的影響呈相反關(guān)系,糧食類消費對農(nóng)村居民蔬菜類消費碳排放有增加作用,對城鎮(zhèn)居民蔬菜類消費碳排放有減緩作用。對于中分位蔬菜類消費碳排放群體,糧食類、蛋類、油脂類消費均會導(dǎo)致城鄉(xiāng)居民蔬菜類消費碳排放增長,肉禽類消費會減緩城鄉(xiāng)居民蔬菜類消費碳排放。從城鄉(xiāng)差異上看,中低分位蔬菜類消費碳排放群體的城鄉(xiāng)差異較大,主要體現(xiàn)在糧食類消費上,農(nóng)村居民的糧食類消費相較于城鎮(zhèn)居民更能導(dǎo)致其蔬菜類消費碳排放增長。對于高分位蔬菜類消費碳排放群體,蛋類、水產(chǎn)品、油脂類消費能引發(fā)長三角城鄉(xiāng)居民蔬菜類消費碳排放增長,肉禽類消費對城鄉(xiāng)居民蔬菜類消費碳排放有減緩作用。高分位蔬菜類消費碳排放群體的城鄉(xiāng)差異較小,主要體現(xiàn)在糧食類和肉禽類消費上,農(nóng)村居民在這兩類食物上的消費相較于城鎮(zhèn)居民更能減緩其蔬菜類消費碳排放。綜上,糧食類消費會導(dǎo)致低分位農(nóng)村居民群體和高分位城鎮(zhèn)居民群體蔬菜類消費碳排放的增長,并且蛋類和油脂類消費會導(dǎo)致長三角城鄉(xiāng)居民蔬菜類消費碳排放增長。

表6 長三角城鄉(xiāng)居民蔬菜類消費碳排放的無條件分位數(shù)回歸結(jié)果
為進一步驗證結(jié)果的穩(wěn)健性,本文通過替換和刪除控制變量兩種方法進行穩(wěn)健性檢驗,具體的做法是將人口撫養(yǎng)比這個變量拆分為老年撫養(yǎng)比和兒童撫養(yǎng)比兩個變量重新估計模型,以及將老齡化率和人口撫養(yǎng)比這兩個控制變量去除后重新估計模型,如表7所示。對控制變量進行處理后,估計結(jié)果與系數(shù)顯著性與各個模型的總體回歸大體相似,表明本文的估計較為穩(wěn)健。

表7 穩(wěn)健性檢驗

(續(xù)表7)
本文基于1995—2019年長三角城鄉(xiāng)居民7大類食物消費數(shù)據(jù),采用無條件分位數(shù)回歸探究長三角城鄉(xiāng)居民食物消費對其食物消費碳排放差異的影響。得到以下結(jié)論:
第一,長三角居民食物消費碳排放總量逐漸增加并逐漸呈現(xiàn)城鄉(xiāng)二元性,但在人均食物消費碳排放的表現(xiàn)上看二者呈現(xiàn)趨同性。第二,從不同分位的城鄉(xiāng)差異視角上看,低分位食物消費碳排放群體的城鄉(xiāng)差異較大,高分位食物消費碳排放群體的城鄉(xiāng)差異較小。并且低分位食物消費碳排放群體的城鄉(xiāng)差異主要體現(xiàn)在肉禽類消費上,高分位食物消費碳排放群體的城鄉(xiāng)差異主要體現(xiàn)在干鮮瓜果類消費上。第三,對3類基礎(chǔ)性食物而言,相較于農(nóng)村居民,城鎮(zhèn)居民的食物消費能顯著引發(fā)其糧食類消費碳排放增加。與總體情況類似,低分位肉禽類和蔬菜類消費碳排放群體的城鄉(xiāng)差異較大,高分位肉禽類和蔬菜類消費碳排放群體的城鄉(xiāng)差異較小。并且糧食類消費會帶動低分位農(nóng)村居民和高分位城鎮(zhèn)居民兩個群體的肉禽類和蔬菜類消費碳排放增加。因此,今后需要特別關(guān)注低分位農(nóng)村群體的食物消費碳排放,控制其肉禽類和糧食類消費,以及高分位城鎮(zhèn)群體的食物消費碳排放,控制其干鮮瓜果類和糧食類消費。
為更好地實現(xiàn)“雙碳”目標,基于前文的研究結(jié)果,提出以下建議:(1)向大眾提倡健康低碳飲食,在保障居民營養(yǎng)正常攝入的前提下,倡導(dǎo)居民適當控制食物消費量,減少食物浪費,同時改善居民食物消費結(jié)構(gòu),實現(xiàn)營養(yǎng)健康與環(huán)境可持續(xù)的雙贏。(2)人口特征變量對居民整體食物消費碳排放的影響值得深入探討。(3)未來需關(guān)注中國其他地區(qū)或者是貧困地區(qū)城鄉(xiāng)居民的食物消費碳排放變化,為全國范圍內(nèi)的城鄉(xiāng)融合以及共同富裕的推進工作做出貢獻。