中國畜牧業(yè)碳排放地區(qū)差異、動態(tài)演進與收斂分析

張金鑫 王紅玲

摘要：中國在經(jīng)過了40多年的改革開放之后，已經(jīng)成為全球第二大經(jīng)濟體，同時也是全球第一碳排放大國，造成的碳排放量占到全球碳排放量的28%，明確中國畜牧業(yè)碳排放地區(qū)差異、動態(tài)演進規(guī)律和是否存在收斂性是制定切實有效減排政策的必要前提。本文運用泰爾指數(shù)和核密度實證分析表明：中國東部、中部、西部地區(qū)和全國畜牧業(yè)碳排放量均在2006年達到頂峰，畜牧業(yè)碳排放總體差異呈現(xiàn)出波動性小幅度增長，區(qū)域內(nèi)差異成為總體差異的主要來源;全國及東部、中部和西部地區(qū)畜牧業(yè)碳排放核密度估計都存在明顯的右拖現(xiàn)象，核密度函數(shù)中心呈現(xiàn)先右移后左移的現(xiàn)象，全國和中部具有明顯的梯度差異，而東部和西部則呈現(xiàn)“單極化”的趨勢;全國及東部、中部和西部地區(qū)畜牧業(yè)碳排放不存在δ收斂和β收斂。中國畜牧業(yè)碳排放地區(qū)差異顯著，各級政府要重視這個典型事實，深挖差異的來源，因地制宜制定符合本地區(qū)實際情況的政策來促進畜牧業(yè)碳排放量減少，著力改善畜牧業(yè)生產(chǎn)環(huán)境，提升畜牧業(yè)產(chǎn)出效率。

關(guān)鍵詞：生態(tài)環(huán)境保護;畜牧業(yè);碳排放;地區(qū)差異

基金項目：國家社會科學(xué)基金重大項目“氣候智慧型農(nóng)業(yè)碳減排及碳交易市場機制研究”（項目編號：19ZDA085）;國家自然科學(xué)基金項目“氣候智慧型農(nóng)業(yè)項目評價的實證研究——基于世界銀行在中國的示范項目”（項目編號：71871086）

中圖分類號：F307 ? ?文獻標(biāo)識碼：A ? ?文章編號：1003-854X（2020）09-0041-08

一、引言

全球性氣候變暖是一個對人類生死存亡有重要影響的環(huán)境問題。聯(lián)合國政府間氣候變化專門委員會（IPCC，下同）的研究報告指出造成全球氣候變暖的主要原因是全球溫室氣體排放量的持續(xù)增加。各國政府為了可持續(xù)發(fā)展，采取很多措施來減少溫室氣體排放，控制氣候變暖。中國在經(jīng)過了40多年的改革開放之后，已經(jīng)成為全球第二大經(jīng)濟體，同時也是全球第一碳排放大國，造成的碳排放占到全球碳排放的28%①。大量溫室氣體的排放不僅給環(huán)境造成巨大破壞，也成為經(jīng)濟高質(zhì)量發(fā)展的一大阻礙。因此，節(jié)能減排、降耗增效成為中國經(jīng)濟調(diào)整結(jié)構(gòu)、增加效率、提升質(zhì)量的必由之路。

農(nóng)業(yè)部門是僅次于能源部門的人為碳排放源。根據(jù)聯(lián)合國糧農(nóng)組織2013年的測算，畜牧業(yè)造成的碳排放在全球碳排放總量的比例高達14.5%②。我國作為世界畜牧業(yè)生產(chǎn)第一大國，有效減少畜牧業(yè)碳排放不僅僅是負(fù)責(zé)任大國的有力擔(dān)當(dāng)，也是實現(xiàn)我國畜牧業(yè)高質(zhì)量發(fā)展的必然要求。盡管我國政府為減少碳排放量做出了巨大努力，積極履行《巴黎協(xié)定》規(guī)定的減排義務(wù)，并在2016年出臺了《“十三五”控制溫室氣體排放工作方案》，但在畜牧業(yè)碳減排上仍然面臨著巨大壓力。畜牧業(yè)遍及東部、中部和西部地區(qū)31個省（市），不同省（市）由于生產(chǎn)要素、歷史條件及經(jīng)濟發(fā)展水平的差異可能造成畜牧業(yè)碳排放存在地區(qū)差異，有不同的演進趨勢。因此，對畜牧業(yè)碳排放實際情況進行充分的研究——畜牧業(yè)碳排放的動態(tài)分布、不同地區(qū)間是否存在差異、差異的來源以及是否收斂等是我國制定科學(xué)、有效畜牧業(yè)減排政策的必要前提。

目前，學(xué)者對于我國畜牧業(yè)碳排放的研究更多是關(guān)注微觀方面，例如研究不同飼養(yǎng)方式或不同物種的碳足跡等③。這些微觀的研究有助于人們科學(xué)地認(rèn)識中國畜牧業(yè)碳排放問題，但卻不能為制定科學(xué)的畜牧業(yè)宏觀碳減排政策提供直接的參考。另外一些學(xué)者則對中國畜牧業(yè)碳排放進行宏觀層面的研究。他們多是基于這樣的思路：首先基于政府間氣候變化專門委員會（IPCC）推薦的排放因子法 ④ 或生命周期法 ⑤ 測算畜牧業(yè)碳排放;然后通過描述性統(tǒng)計，比較不同地區(qū)或者同一地區(qū)不同時間的畜牧業(yè)碳排放情況，以此來分析中國畜牧業(yè)碳排放的地區(qū)差異⑥。學(xué)者們的研究結(jié)果多數(shù)都顯示出中國畜牧業(yè)碳排放地區(qū)差異明顯，但在分析結(jié)果上僅僅只能顯示出畜牧業(yè)碳排放量較高的省份與較低的省份之間是存在顯著差異的，至于差異的變化趨勢以及來源則鮮有論及。

學(xué)者們對中國畜牧業(yè)碳排放的研究雖然取得了較為豐富的成果，但還有一定的拓展空間。從研究深度來看，現(xiàn)有的研究缺乏一定的深度。對于畜牧業(yè)碳排放的分析，一般只是在測算畜牧業(yè)碳排放之后，簡單進行時間和空間的比較，缺乏深入、系統(tǒng)的分析，例如對畜牧業(yè)碳排放演進規(guī)律的分析，難以對畜牧業(yè)碳排放分布規(guī)律進行科學(xué)的把握。從研究方法來看，現(xiàn)有的研究主要是簡單的描述性統(tǒng)計，這盡管可以在一定程度上發(fā)現(xiàn)我國畜牧業(yè)地區(qū)差異情況，但是難以進行細(xì)致化、系統(tǒng)化地考察。尤其是這種簡單的描述統(tǒng)計無法對畜牧業(yè)碳排放地區(qū)差異來源與收斂性進行有效地分析。因此，有必要采用更科學(xué)、嚴(yán)謹(jǐn)?shù)姆椒▽ξ覈竽翗I(yè)碳排放的現(xiàn)狀和基本特征進行更為深入的研究。

二、研究方法及數(shù)據(jù)來源

（一）畜牧業(yè)碳排放量測算與模型構(gòu)建

可以看出：我國東部、中部、西部地區(qū)和全國畜牧業(yè)碳排放量均在2006年達到頂峰，從總體來看，東部、中部地區(qū)和全國在1997—2017年之間畜牧業(yè)碳排放量是呈下降趨勢的，但西部地區(qū)則呈現(xiàn)上升趨勢的;從三個地區(qū)畜牧業(yè)碳排放的占比來看，在1997—2017年之間，東部地區(qū)畜牧業(yè)碳排放在全國畜牧業(yè)碳排放所占的比例總體上是下降的，中部地區(qū)總體上是略微增長的，西部地區(qū)總體上是增長的。

（二）畜牧業(yè)碳排放地區(qū)差異

下表4和下圖2描述了我國畜牧業(yè)碳排放地區(qū)差異情況。在1997—2017年畜牧業(yè)碳排放泰爾指數(shù)從0.2297增長到0.2456，年均增長率為0.34%，表明畜牧業(yè)碳排放總體差異呈現(xiàn)出波動性增長，但增長的幅度較小。畜牧業(yè)碳排放區(qū)域內(nèi)差異貢獻值呈現(xiàn)出先增長后下降的趨勢，區(qū)域間差異貢獻值呈現(xiàn)先下降后增長的趨勢。在1997—2017年間，畜牧業(yè)碳排放區(qū)域內(nèi)差異對總體差異的貢獻率均超過80%，表明區(qū)域內(nèi)差異成為畜牧業(yè)碳排放地區(qū)差異的主要來源。

從三大區(qū)域來看，東部地區(qū)畜牧業(yè)碳排放泰爾指數(shù)呈現(xiàn)出“增長—下降—增長”趨勢，具體是：在1997—2000年間，東部地區(qū)泰爾指數(shù)從1997年的0.4189增長到2000年的0.4406，增長了0.0217，年均增長率為1.70%;在2001—2010年，東部地區(qū)泰爾指數(shù)整體上呈現(xiàn)下降趨勢，從2001年的0.4251下降至2010年的0.3572，下降了0.0679，年均下降率為1.92%;在2011—2017年，東部地區(qū)泰爾指數(shù)開始觸底反彈呈現(xiàn)持續(xù)增長趨勢，從2011年的0.3686增長到2017年的0.4219，增長了0.0533，年均增長率為2.28%，但總體上東部地區(qū)畜牧業(yè)碳排放泰爾指數(shù)呈現(xiàn)增長趨勢。中部畜牧業(yè)碳排放泰爾指數(shù)在1997—2017年之間是呈現(xiàn)波動的“增長—下降”趨勢，1997—2009年是波動增長階段，2010—2017年是波動增長階段，在2009年泰爾指數(shù)達到最大值0.1601。總體上，東部地區(qū)畜牧業(yè)碳排放泰爾指數(shù)呈微弱增長趨勢，由1997年的0.0936增長到2017年的0.1141，增長了0.0205，年均增長率為1.00%。西部地區(qū)畜牧業(yè)碳排放泰爾指數(shù)在1997—2017年出“平穩(wěn)—波動—平穩(wěn)”的趨勢，其中2007—2011年是波動階段，在這個階段西部地區(qū)畜牧業(yè)碳排放泰爾指數(shù)呈現(xiàn)先升后降的波動，而其他年份泰爾指數(shù)幾乎是處于平穩(wěn)的狀態(tài)，且總體上呈現(xiàn)下降趨勢。從泰爾指數(shù)的均值來看，東部、中部和西部地區(qū)的均值分別為0.4076、0.1262和0.1545，說明東部地區(qū)畜牧業(yè)碳排放地區(qū)差異最大，其次是西部地區(qū)，差異最小的是中部地區(qū)。

（三）畜牧業(yè)碳排放動態(tài)演進規(guī)律

根據(jù)公式（3）、（4），并利用stata軟件可以對中國及東部、中部和西部地區(qū)畜牧業(yè)碳排放進行核密度估計（參見下圖3—6），以分析畜牧業(yè)碳排放的動態(tài)演進規(guī)律。

從總體上來看，中國畜牧業(yè)碳排放呈現(xiàn)以下的演進規(guī)律：首先，存在明顯的右拖現(xiàn)象，且1997—2007年右拖范圍明顯擴大，而2007—2017年右拖范圍明顯縮小，這說明畜牧業(yè)碳排放量高的省（市）和低的?。ㄊ校┲g的差異明顯，且相互之間的差異經(jīng)歷了擴大到縮小的過程;其次，從波峰來看，畜牧業(yè)碳排放的核密度逐漸在右側(cè)形成一個小波峰，說明畜牧業(yè)碳排放逐漸呈現(xiàn)多極化的趨勢;第三，核密度函數(shù)中心呈現(xiàn)先右移后左移的現(xiàn)象，說明各個?。ㄊ校┬竽翗I(yè)碳排放量先增加后減少;最后，核密度主峰高度在1997—2017年經(jīng)歷了先減少后增加，說明各個?。ㄊ校┬竽翗I(yè)碳排放之間的差異先增長后下降。

從下圖4可以看出東部地區(qū)畜牧業(yè)碳排放存在以下的演進規(guī)律：首先，東部地區(qū)畜牧業(yè)碳排放核密度估計存在明顯的右拖現(xiàn)象，且在1997—2017年右拖范圍先擴大后縮小，說明東部地區(qū)畜牧業(yè)碳排放量高的?。ㄊ校┖偷偷氖。ㄊ校┲g的差異明顯，但高的?。ㄊ校┖偷偷氖。ㄊ校┲g的差異程度是先增加后減少;其次，從波峰來看，東部地區(qū)畜牧業(yè)碳排放經(jīng)歷了“雙峰”到“單峰”的演進過程，即1997年核密度估計曲線在右側(cè)有1個明顯的小波峰，隨后，這個小波峰逐漸消失，說明東部地區(qū)畜牧業(yè)碳排放逐漸從多極化演變?yōu)閱螛O化的趨勢;第三，核密度函數(shù)中心呈現(xiàn)先右移后左移的現(xiàn)象，說明東部地區(qū)各個?。ㄊ校┬竽翗I(yè)碳排放量先增加后減少;最后，核密度主峰高度在1997—2017年先減少后增加，但曲線寬度范圍有明顯擴大的趨勢，這說明東部地區(qū)各個?。ㄊ校┬竽翗I(yè)碳排放之間的差異變化可能較小。

從下圖5可以看出中部地區(qū)畜牧業(yè)碳排放存在以下的演進規(guī)律：首先，中部地區(qū)畜牧業(yè)碳排放核密度估計存在明顯的右拖現(xiàn)象，且右拖范圍先擴大后縮小在微弱擴大，但總體上是呈現(xiàn)縮小的趨勢，說明中部地區(qū)畜牧業(yè)碳排放量高的?。ㄊ校┖偷偷氖。ㄊ校┲g有明顯的差異，但高的?。ㄊ校┖偷偷氖。ㄊ校┲g的差異程度在總體上呈現(xiàn)縮小的趨勢;其次，從波峰來看，中部地區(qū)畜牧業(yè)碳排放表現(xiàn)明顯的“雙峰”和“多峰”的現(xiàn)象，說明中部地區(qū)各個?。ㄊ校┬竽翗I(yè)碳排放在1997—2017年呈現(xiàn)多極化的趨勢;第三，核密度函數(shù)中心呈現(xiàn)先右移后左移的現(xiàn)象，說明中部地區(qū)各個?。ㄊ校┬竽翗I(yè)碳排放量先增加后減少;最后，從總體上來看，核密度主峰高度下降，曲線寬度變化較為弱小，說明中部地區(qū)各個?。ㄊ校┬竽翗I(yè)碳排放之間的差異有一定程度的增加趨勢。

從下圖6可以看出西部地區(qū)畜牧業(yè)碳排放存在以下的演進規(guī)律：首先，西部地區(qū)畜牧業(yè)碳排放核密度估計在1997年存在明顯的右拖現(xiàn)象，在2002年、2007年、2012年、和2017年有明顯的右拖現(xiàn)象，且右拖范圍先擴大后縮小，說明西部地區(qū)各個?。ㄊ校┬竽翗I(yè)碳排放差異呈現(xiàn)擴大趨勢，并且量高的省（市）和低的?。ㄊ校┲g有明顯的差異，但高的?。ㄊ校┖偷偷氖。ㄊ校┲g的差異程度呈現(xiàn)先增加后減小的趨勢;其次，從波峰來看，西部地區(qū)畜牧業(yè)碳排放表現(xiàn)出明顯的“單峰”，說明西部地區(qū)各個?。ㄊ校┬竽翗I(yè)碳排放在1997—2017年表現(xiàn)出單一化的規(guī)律;第三，核密度函數(shù)中心呈現(xiàn)先右移后左移的現(xiàn)象，說明西部地區(qū)各個省（市）畜牧業(yè)碳排放量先增加后減少，并總體上呈現(xiàn)較少趨勢;最后，從總體上來看，核密度主峰高度表現(xiàn)出下降的現(xiàn)象，曲線寬度也呈現(xiàn)出擴大現(xiàn)象，說明西部地區(qū)各個?。ㄊ校┬竽翗I(yè)碳排放量之間的差異是呈現(xiàn)增加趨勢的。

（四）中國畜牧業(yè)碳排放收斂性分析

本文采取變異系數(shù)來分析中國畜牧業(yè)碳排放是否存在δ收斂，結(jié)果如下圖7所示。從圖7可以看出東部、中部和西部地區(qū)以及全國畜牧業(yè)碳排放變異系數(shù)呈現(xiàn)出微弱的增長趨勢，畜牧業(yè)碳排放變異系數(shù)由1997年的0.6642增長到2017年的0.6908，年均增長率為0.20%;東部地區(qū)畜牧業(yè)碳排放變異系數(shù)由1997年的0.9719增長到2017年的0.9760，年均增長率為0.02%;中部地區(qū)畜牧業(yè)碳排放變異系數(shù)由1997年的0.4835增長到2017年的0.5185，年均增長率為0.35%;西部地區(qū)畜牧業(yè)碳排放變異系數(shù)由1997年的0.5421增長到2017年的0.5589，年均增長率為0.15%。變異系數(shù)的分析結(jié)果表明中國的東部、中部和西部地區(qū)以及全國畜牧業(yè)碳排放不存在δ收斂。

根據(jù)公式（9）可以計算出絕對β收斂的估計結(jié)果，具體如表5所示。從表5的數(shù)據(jù)可以看出東部、中部和西部地區(qū)以及全國畜牧業(yè)碳排放β值均不顯著，說明東部、中部和西部地區(qū)以及全國畜牧業(yè)碳排放不存在絕對β收斂。

條件β收斂的估計結(jié)果如下表6所示。從表6的數(shù)據(jù)可以看出，β+1的系數(shù)是顯著的，但是β+1>0，這說明東部、中部和西部地區(qū)以及全國畜牧業(yè)碳排放不存在條件β收斂。

四、簡要研究結(jié)論與啟示

隨著中國經(jīng)濟的不斷發(fā)展，居民對于各類畜禽類的產(chǎn)品的需求日益增加，激發(fā)了中國畜牧業(yè)的高速發(fā)展，也帶了一系列的環(huán)境問題。其中，因為全球氣候變暖的局勢日益嚴(yán)峻，畜牧業(yè)碳排放問題備受關(guān)注。本文利用中國31個?。ㄊ校?997—2017年畜牧業(yè)數(shù)據(jù)測算畜牧業(yè)碳排放水平，運用泰爾指數(shù)對畜牧業(yè)碳排放地區(qū)差異進行分析，厘清了中國畜牧業(yè)碳排放地區(qū)差異的來源，并結(jié)合核密度分析方法對畜牧業(yè)碳排放動態(tài)演進規(guī)律進行分析。本文研究的基本結(jié)論如下：

（1）中國東部、中部、西部地區(qū)和全國畜牧業(yè)碳排放量均在2006年達到頂峰;在1997—2017年之間，東部地區(qū)畜牧業(yè)碳排放在全國畜牧業(yè)碳排放所占的比例在總體上是下降的，中部地區(qū)總體上是略微增長，西部地區(qū)總體上是增長。

（2）畜牧業(yè)碳排放總體差異呈現(xiàn)出波動性增長，但增長的幅度較小。區(qū)域內(nèi)差異是中國畜牧業(yè)碳排放地區(qū)差異的主要來源，畜牧業(yè)碳排放區(qū)域內(nèi)差異貢獻值呈現(xiàn)出先增長后下降的趨勢，區(qū)域間差異貢獻值呈現(xiàn)先下降后增長的趨勢。分區(qū)域來看，東部、西部和中部三個地區(qū)畜牧業(yè)碳排放地區(qū)差異平均值依次減小。

（3）全國及東部、中部和西部地區(qū)畜牧業(yè)碳排放核密度估計都呈現(xiàn)存在明顯的右拖現(xiàn)象，核密度函數(shù)中心先右移后左移的現(xiàn)象。全國和中部地區(qū)畜牧業(yè)碳排放核密度估計呈現(xiàn)“雙峰”和“多峰”趨勢，東部地區(qū)則經(jīng)歷“雙峰”向“單峰”的演變，西部地區(qū)則一直呈現(xiàn)“單峰”形態(tài)。核密度分析結(jié)果顯示全國及東部、中部和西部地區(qū)畜牧業(yè)碳排放在1997—2017年出現(xiàn)先增加后減少的動態(tài)變化，且畜牧業(yè)碳排放水平高的省（市）和水平低的?。ㄊ校┲g差異明顯，全國和中部地區(qū)具有明顯的梯度差異，而東部和西部地區(qū)則呈現(xiàn)“單極化”的趨勢。

（4）全國及東部、中部和西部地區(qū)畜牧業(yè)碳排放變異系數(shù)在1997—2017年間沒有縮小，說明全國及東部、中部和西部地區(qū)畜牧業(yè)碳排放不存在δ收斂。絕對收斂和條件收斂的分析結(jié)果顯示全國及東部、中部和西部地區(qū)畜牧業(yè)碳排放不存在δ收斂和β收斂。

基于上述研究結(jié)論，可以有如下啟示：

第一，要清晰認(rèn)識到中國畜牧業(yè)碳排放的地區(qū)差異，減少畜牧業(yè)碳排放量，提升畜牧業(yè)環(huán)境友好度是中國未來相當(dāng)一段時間內(nèi)畜牧業(yè)發(fā)展的基本要求。從本文的研究結(jié)果來看，中國畜牧業(yè)碳排放地區(qū)差異顯著，各級政府要重視這個典型事實，深挖差異的來源，因地制宜制定符合本地區(qū)實際情況的政策來促進畜牧業(yè)碳排放量減少。

第二，全國及東部、中部和西部地區(qū)畜牧業(yè)碳排放已經(jīng)在2006年達到峰值，但存在明顯的梯度格局?；谶@個事實，中國政府需要制定長效機制來保障畜牧業(yè)碳排放平穩(wěn)有序減少，并針對畜牧業(yè)排放的大省制定個性化政策，提高其畜牧業(yè)生產(chǎn)效率，減少畜牧業(yè)碳排放量。

第三，全國及東部、中部和西部地區(qū)畜牧業(yè)碳排放均不存在δ收斂和β收斂，這反映出畜牧業(yè)經(jīng)濟發(fā)展的不平衡性。中國政府需要采取積極的措施來引導(dǎo)畜牧業(yè)高質(zhì)量發(fā)展，著力改善畜牧業(yè)生產(chǎn)環(huán)境，例如通過推廣氣候智慧型畜牧業(yè)發(fā)展綠色、有機畜牧業(yè)，減少反芻動物腸道發(fā)酵和糞便發(fā)酵的溫室氣體排放，提升畜牧業(yè)產(chǎn)出效率。

注釋：

① W. Li， C. An， C. Lu， The Assessment Framework of Provincial Carbon Emission Driving Factors： An Empirical Analysis of Hebei Province， Science of the Total Environment， 2018， 637， pp.91-103.

②⑦⑧ 胡向東、王濟民：《中國畜禽溫室氣體排放量估算》，《農(nóng)業(yè)工程學(xué)報》2010年第1期。

③ M. Zhuang， W. Li， Greenhouse Gas Emission of Pastoralism is Lower than Combined Extensive/intensive Livestock Husbandry： A Case Study on the Qinghai-Tibet Plateau of China， Journal of Cleaner Production， 2017， 147， pp.514-522.

④ 陳蘇、胡浩：《中國畜禽溫室氣體排放時空變化及影響因素研究》，《中國人口·資源與環(huán)境》2016年第7期。

⑤ 孟祥海、程國強、張俊飚、王宇波、周海川：《中國畜牧業(yè)全生命周期溫室氣體排放時空特征分析》，《中國環(huán)境科學(xué)》2014年第8期。

⑥ 姚成勝、錢雙雙、毛躍華、李政通：《中國畜牧業(yè)碳排放量變化的影響因素分解及空間分異》，《農(nóng)業(yè)工程學(xué)報》2017年第12期。

⑨ 劉培芳、陳振樓、許世遠(yuǎn)、劉杰：《長江三角洲城郊畜禽糞便的污染負(fù)荷及其防治對策》，《長江流域資源與環(huán)境》2002年第5期。

⑩ 根據(jù)國家統(tǒng)計局的劃分標(biāo)準(zhǔn)：東部地區(qū)包括北京、天津、河北、遼寧、上海、江蘇、浙江、福建、山東、廣東、廣西、海南共12個省（市）;中部地區(qū)包括山西、內(nèi)蒙古、吉林、黑龍江、安徽、江西、河南、湖北、湖南共9個?。ㄊ校?西部地區(qū)包括四川、重慶、貴州、云南、西藏、陜西、甘肅、寧夏、青海、新疆10個?。ㄊ校?。

作者簡介：張金鑫，湖北大學(xué)商學(xué)院博士研究生，湖北武漢，430062;王紅玲，湖北大學(xué)商學(xué)院教授、博士生導(dǎo)師，湖北武漢，430062。

（責(zé)任編輯 ?陳孝兵）