長三角地區城鄉居民的食物消費碳排放差異研究

呂夢姍 李 雷，2，3 王成軍，2，3

(1.浙江農林大學經濟管理學院；2.浙江農林大學鄉村振興研究院；3.浙江農林大學碳中和研究院，浙江 杭州 311300)

改革開放以來，中國城鄉居民收入和消費日漸增長，居民的食物消費正從溫飽型消費向發展型消費過渡[1]，食物消費總量顯著提升并呈現多元化趨勢。居民食物消費結構的變化，可能會導致食物消費碳排放的變化，影響資源環境的可持續性。中國作為一個負責任的發展中國家，在2020年提出了“雙碳”目標。作為改善居民健康和協調自然環境系統可持續性的有效途徑，近年來食物系統受到廣泛關注?！吨袊r業產業發展報告2022》指出要加速推進肉類消費向健康營養低碳方向轉型，促進農食系統的碳減排，進而提高中國農食系統的韌性。食物系統作為溫室氣體排放的主要來源，約占全球碳排放的25%～30%，而食物消費又與居民日常生活息息相關[2]，食物作為連接城鄉的載體，起著非常重要的作用，因此向大眾提倡低碳綠色的食物消費方式至關重要[3]。長三角地區作為我國最發達的經濟圈之一，該地區居民快速的生活節奏導致其食物消費模式也在快速變化，這給中國實現“雙碳”目標帶來極大挑戰。長三角地區有責任探索低碳可持續的食物消費路徑，率先推進全面綠色低碳轉型。并且長三角已經開啟探索共同富裕的新模式，該地區城鄉居民收入差距較中國其他地區小，研究該地區城鄉居民的食物消費碳排放變化差異，能為促進城鄉協調發展提供參考，具有較強的現實意義。同時本文的理論意義在于在可持續發展理論的基礎上，豐富了農食系統維度上的可持續發展理論，進一步探討食物消費與碳排放二者的邏輯關系。本文旨在“雙碳”目標背景下探究長三角城鄉居民的食物消費碳排放差異以及不同類型食物的人均消費量對長三角城鄉居民的食物消費碳排放影響，對明確居民食物消費碳排放變化和實現食物消費可持續發展有重要意義。

一、文獻綜述

居民食物消費碳排放，是指居民在日常的食物消費中所產生的二氧化碳量，來源包括食物本身、食物的生產加工運輸、食物的儲存與烹飪等環節[4]。國內外學者就居民食物消費碳排放，運用生命周期法[5]、投入產出法[6]、碳折算系數法[7]等方法進行測算研究。目前關于食物消費碳排放的研究主要集中在以下三個方面：(1)食物消費碳排放的時空演變研究。就中國而言，中國各地食物消費人均碳排放呈現較顯著的民族性、地域空間差異特征[8]。從城鄉分異視角上看，近年來部分學者通過研究發現：城鎮居民在食物消費碳排放總量上高于農村居民。一方面，隨著越來越多的農村人口進入城市，將促進城鎮居民食物消費碳排放的增加，另一方面，城鄉居民食物消費結構多樣化導致食物消費碳排放較之前年份有所下降[9，10]?？傮w上看，中國城鄉居民的食物碳消費同時呈現二元性和趨同性特征[8，11，12]。(2)食物消費碳排放的影響因素研究。有部分學者指出，GDP指數、人口規模、城鎮化、居民收入等因素會對食物消費碳排放帶來不同程度的影響[13-15]。(3)食物消費結構變動與碳排放關系。張翠玲等綜合考慮食物消費的經濟成本最低、水資源消耗最少、耕地資源消耗最少和溫室氣體排放最少4個目標對當前中國居民食物消費結構進行優化[16]；韓金雨等通過VAR模型探尋食物消費結構的升級與農業碳排放之間的關系，發現兩者存在一定的關聯，且具有滯后性，同時指出農業碳排放對植物性食物消費的沖擊更為靈敏[17]。然而目前對全國最發達的經濟區——長三角地區居民的食物消費碳排放的研究比較少見，這不利于準確評估經濟快速發展背景下的中國城鄉居民食物消費碳排放的變動情況，為此，本文從中國長三角地區的1995—2019年居民食物消費數據入手，分析長三角城鄉居民食物消費碳排放整體情況，并基于無條件分位數回歸(UQR)，研究不同類型食物的人均消費量對長三角城鄉居民食物消費碳排放的影響，為實現中國城鄉居民食物消費可持續發展提供參考依據。

二、數據與模型

(一)計量模型的設定

本文設定的計量模型如下，首先是以居民人均食物消費碳排放作為被解釋變量為例：

(1)

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

(2)

(3)

(4)

上式所涉及的具體變量的編碼和含義如表1所示。

表1 變量含義說明

(二)研究方法

無條件分位數回歸(UQR)是對條件分位數回歸(CQR)的補充和拓展。UQR主要應用于勞動經濟學領域，多數學者用此研究不同影響因素對不同群體的收入分布差距的影響[18]。近年來，越來越多的學者將UQR用于其他領域。方超研究信息技術對學生教育的公平性[19]；周春芳和蘇群分析農民工與城鎮職工的就業質量差異[20]；鄭宏遠和李谷成研究數字普惠金融對縣域農業全要素生產率的影響[21]。本文基于Firpo等提供的無條件分位數回歸，利用再中心化影響函數(RIF)進行無條件估計[22]，參考鄭宏遠和李谷成的做法[21]，表示為式(5)：

式1和式5中的LNTPEit均表示人均食物消費碳排放的對數。qτ是人均食物消費碳排放在第τ個分位點的無條件分布，I(LNTPEit≤qτ)是一個區分人均食物消費碳排放比大還是小的二元變量，fLNTPF(qτ)表示在qτ估計的人均食物消費碳排放的概率密度函數。由于本文所利用的是省級面板數據，必須對那些不隨時間變化但可能會影響居民人均食物消費碳排放的因素加以控制?；诖?，本文主要利用BORGEN提出的模型進行估計[23]。

(三)數據來源

本文將食物劃分為7個主要食物類別，選取樣本為1995—2019年長三角地區各省市城鄉居民食物消費的觀察值。其中食物消費量歷年數據取自《中國住戶統計年鑒》《中國統計年鑒》以及各省市統計年鑒。本文的食物消費碳排放系數參照Poore等的食物系統碳排放數據庫的數據[24]，通過整理如表2所示。根據歷年食物消費量和食物消費碳排放系數，通過計算得到某一省份的城鎮或農村居民在某一年的人均食物消費碳排放總量。除此之外，其余的解釋變量則來自《中國人口與就業統計年鑒》與各省市統計年鑒。

表2 食物的碳排放系數

三、實證檢驗與分析

(一)長三角居民食物消費碳排放變化情況

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

圖1 1995-2019年長三角城鄉居民食物消費碳排放總量變化情況

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

圖2 長三角居民1995-2019年人均食物消費碳排放變化趨勢

(二)描述性統計分析

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

(三)無條件分位數回歸

表3 長三角城鄉居民食物消費碳排放的無條件分位數回歸結果

(續表3)

本文考慮城鄉差異，對樣本數據進行回歸，得到以10%、25%、50%、75%、90%為代表的長三角居民食物消費碳排放的無條件分位數固定效應模型回歸結果，如表3所示，可以看出：各類食物的人均消費量大部分能顯著反映長三角城鄉居民食物消費碳排放變化。對于低分位食物消費碳排放群體，糧食類和肉禽類消費均會導致城鄉居民食物消費碳排放增加，蛋類和干鮮瓜果類消費有助于減緩城鄉居民食物消費碳排放。其中低分位食物消費碳排放群體的城鄉差異主要體現在肉禽類消費上，農村居民的肉禽類消費相較于城鎮居民更能導致其食物消費碳排放增加。對于中等分位食物消費碳排放群體，糧食類、肉禽類及水產品消費會導致城鎮居民食物消費碳排放增加，同時糧食類和水產品能減緩農村居民的食物消費碳排放。中等分位食物消費碳排放群體的城鄉差異則主要體現在糧食類消費上。對于高分位食物消費碳排放群體，糧食類和水產品均能導致城鄉居民食物消費碳排放增加，肉類、蛋類和油脂類消費一定程度上減緩城鄉居民的食物消費碳排放。高分位食物消費碳排放群體的城鄉差異主要體現在干鮮瓜果類消費上，并且干鮮瓜果類會導致城鎮居民食物消費碳排放增加。

此外，其他控制變量如人均可支配收入、人口數量、恩格爾系數、老齡化率、人口撫養比等對居民食物消費碳排放均有顯著的異質性影響，考慮到本文重點關注城鄉居民的不同類型食物消費量與其食物消費碳排放的關系，在此不做討論。

本文選取長三角城鄉居民的3類基礎性食物(糧食類、肉禽類、蔬菜類)消費碳排放的無條件分位數回歸結果進行展示，如表4—6所示。

其中表4是針對居民人均糧食類食物消費碳排放進行的實證分析。相較于農村居民，水產品、油脂類、干鮮瓜果類消費增加，會導致城鎮居民的糧食類消費碳排放增加。從無條件分位數回歸上看，對于低分位糧食類消費碳排放群體，蛋類消費會導致長三角城鄉居民糧食類消費碳排放增加，肉禽類、水產品、油脂類消費對城鄉居民糧食類消費碳排放的影響呈相反關系。其中肉禽類和油脂類消費更能引起農村居民糧食類消費碳排放增長，水產品消費則更能引起城鎮居民糧食類消費碳排放增長。從城鄉差異上看，25分位糧食類消費碳排放群體的城鄉差異較大，主要體現在蛋類消費上，農村居民的蛋類消費相較于城鎮居民更能減緩其糧食類消費碳排放。對于中分位糧食類消費碳排放群體，水產品、油脂類、干鮮瓜果類消費對城鄉居民糧食類消費碳排放的影響呈相反關系。其中干鮮瓜果類消費會導致農村居民糧食類消費碳排放增加，水產品和油脂類消費會導致城鎮居民糧食類消費碳排放增加。并且中分位糧食類消費碳排放群體的城鄉差異則主要體現在油脂類消費上。對于高分位糧食類消費碳排放群體，油脂類會導致城鄉居民糧食類消費碳排放均有所增長，肉禽類消費能減緩城鄉居民糧食類消費碳排放。同時，蛋類和干鮮瓜果類對城鄉居民糧食類消費碳排放的影響呈相反關系，蛋類消費對城鎮居民糧食類消費碳排放有減緩作用，干鮮瓜果類消費對農村居民糧食類消費碳排放有減緩作用。并且90分位糧食類消費碳排放群體的城鄉差異也相對較大，主要體現在干鮮瓜果類消費上。綜上，相較于農村居民，城鎮居民通過對不同類型的食物消費，能顯著促進其糧食類消費碳排放。

表4 長三角居民糧食類食物消費碳排放的無條件分位數回歸結果

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

表5 長三角城鄉居民肉禽類食物消費碳排放的無條件分位數回歸結果

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

表6 長三角城鄉居民蔬菜類消費碳排放的無條件分位數回歸結果

(四)穩健性檢驗

為進一步驗證結果的穩健性，本文通過替換和刪除控制變量兩種方法進行穩健性檢驗，具體的做法是將人口撫養比這個變量拆分為老年撫養比和兒童撫養比兩個變量重新估計模型，以及將老齡化率和人口撫養比這兩個控制變量去除后重新估計模型，如表7所示。對控制變量進行處理后，估計結果與系數顯著性與各個模型的總體回歸大體相似，表明本文的估計較為穩健。

表7 穩健性檢驗

(續表7)

四、結論與建議

(一)結論

本文基于1995—2019年長三角城鄉居民7大類食物消費數據，采用無條件分位數回歸探究長三角城鄉居民食物消費對其食物消費碳排放差異的影響。得到以下結論：

第一，長三角居民食物消費碳排放總量逐漸增加并逐漸呈現城鄉二元性，但在人均食物消費碳排放的表現上看二者呈現趨同性。第二，從不同分位的城鄉差異視角上看，低分位食物消費碳排放群體的城鄉差異較大，高分位食物消費碳排放群體的城鄉差異較小。并且低分位食物消費碳排放群體的城鄉差異主要體現在肉禽類消費上，高分位食物消費碳排放群體的城鄉差異主要體現在干鮮瓜果類消費上。第三，對3類基礎性食物而言，相較于農村居民，城鎮居民的食物消費能顯著引發其糧食類消費碳排放增加。與總體情況類似，低分位肉禽類和蔬菜類消費碳排放群體的城鄉差異較大，高分位肉禽類和蔬菜類消費碳排放群體的城鄉差異較小。并且糧食類消費會帶動低分位農村居民和高分位城鎮居民兩個群體的肉禽類和蔬菜類消費碳排放增加。因此，今后需要特別關注低分位農村群體的食物消費碳排放，控制其肉禽類和糧食類消費，以及高分位城鎮群體的食物消費碳排放，控制其干鮮瓜果類和糧食類消費。

(二)建議

為更好地實現“雙碳”目標，基于前文的研究結果，提出以下建議：(1)向大眾提倡健康低碳飲食，在保障居民營養正常攝入的前提下，倡導居民適當控制食物消費量，減少食物浪費，同時改善居民食物消費結構，實現營養健康與環境可持續的雙贏。(2)人口特征變量對居民整體食物消費碳排放的影響值得深入探討。(3)未來需關注中國其他地區或者是貧困地區城鄉居民的食物消費碳排放變化，為全國范圍內的城鄉融合以及共同富裕的推進工作做出貢獻。