陳蘇+胡浩

摘要畜禽業(yè)是關(guān)系國(guó)計(jì)民生的重要產(chǎn)業(yè)，在滿足居民生活物質(zhì)需求的同時(shí)，其生產(chǎn)過程中所伴隨的溫室氣體排放問題已成為不容忽視的環(huán)境問題。本文基于LMDI模型系統(tǒng)分析了1991-2013年中國(guó)畜禽溫室氣體排放時(shí)空變化及其因素貢獻(xiàn)。結(jié)果表明：①?gòu)臅r(shí)間維度來看，1991-2013年，中國(guó)畜禽溫室氣體排放經(jīng)歷了先快速上升后穩(wěn)定上升再波動(dòng)下降的變化特征，總體呈上升趨勢(shì)。經(jīng)濟(jì)效應(yīng)對(duì)畜禽溫室氣體排放促進(jìn)作用最大，而強(qiáng)度效應(yīng)的抑制作用最大，其次是勞動(dòng)力效應(yīng)和結(jié)構(gòu)效應(yīng)。②期間，經(jīng)濟(jì)效應(yīng)促進(jìn)作用的累計(jì)貢獻(xiàn)呈指數(shù)增長(zhǎng)，而強(qiáng)度效應(yīng)抑制作用的累計(jì)貢獻(xiàn)呈倒“U”，勞動(dòng)力效應(yīng)和結(jié)構(gòu)效應(yīng)抑制作用也不斷加強(qiáng)。③從空間維度來看，中國(guó)畜禽溫室氣體排放的區(qū)域集中度較高，四川、河南、山東、云南和內(nèi)蒙古等省（區(qū)、市）畜禽溫室氣體排放一直位居全國(guó)前列。④省域各效應(yīng)作用方向和程度差異顯著，四川、青海和云南強(qiáng)度效應(yīng)抑制作用較大，遼寧、吉林和黑龍江抑制作用較小；山東、四川和黑龍江結(jié)構(gòu)效應(yīng)抑制作用顯著，新疆和青海促進(jìn)作用明顯；四川、河南、內(nèi)蒙古、山東、云南、湖南和河北經(jīng)濟(jì)效應(yīng)促進(jìn)作用較大，天津、上海、海南和北京促進(jìn)作用較??；四川、湖北、江蘇和山東勞動(dòng)力效應(yīng)抑制作用顯著，新疆、黑龍江和內(nèi)蒙古促進(jìn)作用明顯。最后，基于研究結(jié)論，為未來中國(guó)畜禽溫室氣體減排空間發(fā)展策略提出了幾點(diǎn)建議。

關(guān)鍵詞畜禽；溫室氣體；時(shí)空變化；LMDI模型

中圖分類號(hào)S168文獻(xiàn)標(biāo)識(shí)碼A文章編號(hào)1002-2104（2016）07-0093-08doi：10.3969/j.issn.1002-2104.2016.07.012

20世紀(jì)90年代以來，全球氣候變化成為人類經(jīng)濟(jì)社會(huì)可持續(xù)發(fā)展所面臨的重大挑戰(zhàn)，畜禽溫室氣體排放日益受到社會(huì)各界的關(guān)注。聯(lián)合國(guó)糧農(nóng)組織（FAO）2006年的報(bào)告顯示，每年由牛、羊、馬、駱駝、豬和家禽排放溫室氣體的CO2當(dāng)量占全球排放量的18%[1]。而世界觀察研究所2009年的報(bào)告指出，全球牲畜及其副產(chǎn)品排放溫室氣體的CO2當(dāng)量約占全球總排放量的51%[2]，幾乎是FAO估算量的3倍。可見，畜禽已成為重要的溫室氣體排放源，而畜禽溫室氣體主要源于動(dòng)物腸道CH4排放、動(dòng)物糞便處理過程中產(chǎn)生的CH4和N2O[3]，從動(dòng)物類型來看，反芻動(dòng)物產(chǎn)生的溫室氣體排放最多，其次為豬，最少的是雞[4]。

國(guó)內(nèi)外學(xué)者對(duì)畜禽溫室氣體排放量的測(cè)算及其影響因素進(jìn)行了大量研究。在畜禽溫室氣體排放測(cè)算方面，董紅敏[5]等采用OECD的測(cè)算方法對(duì)中國(guó)三個(gè)時(shí)點(diǎn)（1980年、1985年、1990年）的反芻類動(dòng)物CH4排放量進(jìn)行了估算；FAO[1]利用IPCC的方法和系數(shù)，估算了中國(guó)2004年主要畜禽的溫室氣體排放量；Zhou[6]等測(cè)算了中國(guó)1949-2003年畜禽的溫室氣體排放量；胡向東[7]等測(cè)算了中國(guó)2000-2007年以及各省區(qū)2007年畜禽溫室氣體排放量，結(jié)果表明，2000-2007年中國(guó)畜禽溫室氣體排放量總體呈下降趨勢(shì)，各省區(qū)畜禽溫室氣體排放量呈現(xiàn)區(qū)域集中特點(diǎn)；閔繼勝[8]等測(cè)算了中國(guó)1991-2008年以及各省份畜牧業(yè)溫室氣體排放量，結(jié)果表明，1991年以來，中國(guó)畜牧CH4和N2O排放量均呈先升后降的趨勢(shì)；尚杰[9]等測(cè)算了1993-2011年中國(guó)畜禽溫室氣體排放量，結(jié)果表明，中國(guó)畜禽的CH4排放量整體呈波動(dòng)上升趨勢(shì)，N2O排放量持續(xù)增加。在畜禽溫室氣體排放的影響因素方面，譚秋成[10]研究表明，由于技術(shù)進(jìn)步和技術(shù)效率的提高，單位肉類和牛奶排放的溫室氣體均有大幅度下降；陳瑤[11]等研究表明，經(jīng)濟(jì)因素是影響我國(guó)畜牧業(yè)溫室氣體排放的最大因素，短期內(nèi)效率因素是我國(guó)畜牧業(yè)低碳化發(fā)展的最主要誘因，而從長(zhǎng)期來看勞動(dòng)力因素是我國(guó)畜牧業(yè)低碳化發(fā)展的最主要因素；尚杰[9]等研究表明，動(dòng)物腸道發(fā)酵CH4、N2O排放的影響因素主要取決于動(dòng)物種類、飼料特性、飼養(yǎng)方式和糞便管理方式等。

以上研究取得了有價(jià)值的結(jié)論，為本文深入研究提供了重要的參考數(shù)據(jù)和研究方法。但存在以下可以改進(jìn)之處：一是研究對(duì)象大多側(cè)重于國(guó)家層面畜禽溫室氣體排放量的測(cè)算，全面把握中國(guó)畜禽溫室氣體排放變化規(guī)律，不僅從總體上刻畫其演變特征，更要分析區(qū)域差異；二是關(guān)于畜禽溫室氣體排放成因研究未及深入展開，考慮到畜禽溫室氣體排放的區(qū)域差異性，有必要對(duì)各地區(qū)畜禽溫室氣體排放的影響因素進(jìn)行分析，以便找到進(jìn)一步降低畜禽溫室氣體排放的方向和對(duì)策?；诖耍疚臏y(cè)算分析了1991-2013年中國(guó)畜禽溫室氣體時(shí)空變化規(guī)律，并運(yùn)用LMDI模型從溫室氣體排放強(qiáng)度、農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)、農(nóng)業(yè)經(jīng)濟(jì)水平和農(nóng)業(yè)勞動(dòng)力等方面進(jìn)行因素分解，揭示畜禽溫室氣體排放時(shí)空變化的成因。

陳蘇等：中國(guó)畜禽溫室氣體排放時(shí)空變化及影響因素研究中國(guó)人口·資源與環(huán)境2016年第7期1研究方法及數(shù)據(jù)來源

1.1畜禽溫室氣體排放量的測(cè)算方法

畜禽溫室氣體排放主要包括畜禽胃腸道內(nèi)發(fā)酵的CH4、畜禽糞便處理產(chǎn)生的CH4和N2O和畜禽飼養(yǎng)過程中對(duì)化石能源等消耗產(chǎn)生的CO2[12]。鑒于畜禽生產(chǎn)過程中化石能源消耗相關(guān)數(shù)據(jù)的缺乏，本文選取牛、羊、馬、騾、驢、駱駝、生豬、家禽和兔等動(dòng)物作為研究對(duì)象，測(cè)算中國(guó)及各?。▍^(qū)、市）畜禽溫室氣體排放量，其具體的測(cè)算方法如下：

式中，C、CCH4和CN2O分別為畜禽溫室氣體排放量、CH4和N2O排放量；21和310分別為CH4和N2O轉(zhuǎn)化為CO2當(dāng)量的轉(zhuǎn)化系數(shù)；Ni表示第i種畜禽的平均飼養(yǎng)量；αi和βi表示第i種畜禽的CH4和N2O排放因子。由于畜禽飼養(yǎng)周期不同，需要對(duì)畜禽年平均飼養(yǎng)量進(jìn)行調(diào)整，參考胡向東[7]的計(jì)算方法。當(dāng)出欄率大于或等于1時(shí)，畜禽年平均飼養(yǎng)量用出欄量除以365再乘以其生命周期，主要有生豬、家禽和兔，生命周期分別為200天[7]、55天[13]和105天[7]；當(dāng)出欄率小于1時(shí)，畜禽年平均飼養(yǎng)量用本年末的存欄量表示，為消除單個(gè)時(shí)間點(diǎn)的影響，采取畜禽上年年末存欄量和本年末存欄量的平均數(shù)表示。借鑒已有研究關(guān)于各畜禽的溫室氣體排放系數(shù)，CH4排放系數(shù)來源于2006年IPCC國(guó)家間溫室氣體排放指南[14]，N2O排放系數(shù)來源于胡向東[7]，具體的排放系數(shù)見表1。

1.2畜禽溫室氣體排放影響因素的LMDI分解

因素分解方法作為研究事物變化特征及其作用機(jī)理的一種分析框架，在環(huán)境經(jīng)濟(jì)研究中得到廣泛的應(yīng)用。通行的分解方法主要有兩類，一類是指數(shù)分解方法（Index Decomposition Analysis，IDA），另一類是結(jié)構(gòu)分解方法（Structural Decomposition Analysis，SDA）。SDA方法利用投入產(chǎn)出表，以消費(fèi)系數(shù)矩陣為基礎(chǔ)，對(duì)數(shù)據(jù)要求較高；而IDA方法只需部門加總數(shù)據(jù)，適合分解含有較少因素的、包含時(shí)間序列數(shù)據(jù)的模型。IDA方法包括Laspeyres指數(shù)分解與Divisia指數(shù)分解等，但兩者分解不徹底，存在分解剩余項(xiàng)，Ang[15]等在綜合比較了各種IDA方法基礎(chǔ)上，提出了對(duì)數(shù)平均迪氏指數(shù)法（Logarithmic Mean Divisia Index，LMDI），該方法最大特點(diǎn)在于不會(huì)產(chǎn)生分解剩余項(xiàng)，且允許數(shù)據(jù)中包含零值。因此，本文選用LMDI從溫室氣體排放強(qiáng)度、農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)、農(nóng)業(yè)經(jīng)濟(jì)水平和農(nóng)業(yè)勞動(dòng)力等方面量化分解影響畜禽溫室氣體排放的因素[16]。結(jié)合現(xiàn)有研究成果，將畜禽溫室氣體排放分解為：

C=CLS×LSAGRI×AGRIP×P（2）

式（2）中，C為畜禽溫室氣體排放量，LS為畜牧業(yè)產(chǎn)值，AGRI為農(nóng)林牧漁業(yè)總產(chǎn)值，P為農(nóng)業(yè)勞動(dòng)力的數(shù)量。對(duì)各個(gè)分解因素進(jìn)行定義，定義EI=C/LS為畜禽溫室氣體排放強(qiáng)度，即畜禽溫室氣體排放量與畜牧業(yè)產(chǎn)值之比；定義CI=LS/AGRI為農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)，即畜牧業(yè)產(chǎn)值占農(nóng)林牧漁業(yè)總產(chǎn)值比重；定義SI=AGRI/P為農(nóng)業(yè)經(jīng)濟(jì)水平，即農(nóng)業(yè)勞動(dòng)力的人均農(nóng)林牧漁業(yè)產(chǎn)值。則（2）式可進(jìn)一步表述為：

C=EI×CI×SI×P（3）

由于LMDI的“乘積分解”和“加和分解”最終結(jié)果一致，而后者能較為清晰的分解出影響因素，因此，本文采用

放系數(shù)腸道發(fā)酵1.0068.0051.4018.0010.0046.005.000.254-糞便管理3.5016.001.501.640.901.920.160.080.02N2O

排放系數(shù)糞便管理0.531.001.371.391.391.390.330.020.02注：非奶牛取黃牛和水牛的平均值；羊取山羊和綿羊的平均數(shù)；家禽取雞、鴨、鵝和火雞的平均數(shù)?！凹雍头纸狻钡姆椒ǎㄔ敿?xì)推導(dǎo)過程可參閱Ang[17]etc）：

ΔC=Ct-C0=ΔEI+ΔCI+ΔSI+ΔP（4）

式（4）中，C0為基期畜禽溫室氣體排放總量，Ct為T期溫室氣體排放總量，ΔC為畜禽溫室氣體排放總量變化。這種變化可分解為：ΔEI表示單位畜牧業(yè)產(chǎn)值排放溫室氣體變化，即強(qiáng)度效應(yīng)；ΔCI表示單位農(nóng)林牧漁業(yè)總產(chǎn)值的畜牧業(yè)產(chǎn)值變化，即結(jié)構(gòu)效應(yīng)；ΔSI表示人均農(nóng)林牧漁業(yè)總產(chǎn)值變化，即經(jīng)濟(jì)效應(yīng)；ΔP表示農(nóng)業(yè)勞動(dòng)力變化，即勞動(dòng)力效應(yīng)。由此，畜禽溫室氣體變化直接受制于4種因素的變化。其具體表達(dá)式分別為：

若ΔEI、ΔCI、ΔSI和ΔP的系數(shù)為正值，說明該效應(yīng)對(duì)畜禽溫室氣體排放起到促進(jìn)作用，反之，則起到抑制作用。

1.3數(shù)據(jù)來源及整理

本文以生豬、牛、馬、騾、驢、駱駝、羊、兔和家禽為研究對(duì)象，選取30個(gè)?。▍^(qū)、市）（其中重慶市數(shù)據(jù)合并到四川省數(shù)據(jù)內(nèi)）畜禽的出欄量、存欄量、畜牧業(yè)產(chǎn)值、農(nóng)林牧漁業(yè)總產(chǎn)值以及農(nóng)業(yè)勞動(dòng)力數(shù)量等數(shù)據(jù)，這些數(shù)據(jù)來自于《中國(guó)農(nóng)業(yè)年鑒》、《中國(guó)農(nóng)村統(tǒng)計(jì)年鑒》、《中國(guó)畜牧業(yè)年鑒》?？紤]到產(chǎn)值不具有縱向可比性，因此本文中的畜牧業(yè)產(chǎn)值和農(nóng)林牧漁業(yè)總產(chǎn)值以1990年為基準(zhǔn)年，換算為可比的實(shí)際產(chǎn)值。

2結(jié)果分析

2.1中國(guó)畜禽溫室氣體排放時(shí)序變化

2.1.1畜禽溫室氣體排放的階段變化

依據(jù)畜禽溫室氣體排放測(cè)算公式、各個(gè)畜禽溫室氣體排放系數(shù)和畜禽的出欄、存欄相關(guān)數(shù)據(jù)，量化測(cè)算了中國(guó)1991-2013年的畜禽溫室氣體排放情況，并將其轉(zhuǎn)化為CO2當(dāng)量（圖1）。圖1表明，1991-2013年畜禽溫室氣體排放大致分為3個(gè)階段，在此基礎(chǔ)上，各階段溫室氣體排放總量變化及各效應(yīng)的影響程度見表2。

第一階段（1991-1996年），畜禽溫室氣體排放量快速上升。由1991年的2 746.82萬t上升到1996年的3 746.16萬t，增加了999.34萬t。該時(shí)期經(jīng)濟(jì)效應(yīng)是促進(jìn)溫室氣體排放最主要推動(dòng)力為2 254.88萬t；其他對(duì)溫室氣體排放起到抑制作用，其中強(qiáng)度效應(yīng)抑制作用最大，為-939.47萬t，其次是勞動(dòng)力效應(yīng)和結(jié)構(gòu)效應(yīng)，分別為圖11991-2013年中國(guó)畜禽溫室氣體排放

總量變化趨勢(shì)

第二階段（1997-2006年），畜禽溫室氣體排放量穩(wěn)定上升。受金融危機(jī)、通貨緊縮等因素影響，1997年畜禽平均飼養(yǎng)量較上一年大幅度下降，強(qiáng)度效應(yīng)抑制作用為-451.53萬t，經(jīng)濟(jì)效應(yīng)抑制作用為-202.35萬t，實(shí)現(xiàn)了492.17萬t畜禽溫室氣體的減排，隨后逐年增加，到2006年畜禽溫室氣體排放總量達(dá)到峰值，為4 228.50萬t，增加了482.34萬t（需要說明的是：這里峰值出現(xiàn)的時(shí)間與胡向東等測(cè)算的結(jié)果不同，主要原因是后者2006年畜禽數(shù)據(jù)根據(jù)第二次農(nóng)業(yè)普查結(jié)果進(jìn)行了調(diào)整，而本文畜禽數(shù)據(jù)來源于《中國(guó)農(nóng)業(yè)年鑒》，以保證數(shù)據(jù)來源的統(tǒng)一性）。該時(shí)期經(jīng)濟(jì)效應(yīng)對(duì)溫室氣體排放促進(jìn)作用最大，為801.21萬t，其次是強(qiáng)度效應(yīng)，為171.18萬t。勞動(dòng)力效應(yīng)和結(jié)構(gòu)效應(yīng)對(duì)溫室氣體排放起到不同程度的抑制作用，分別為-329.14萬t和-160.91萬t。

第三階段（2007-2013年），畜禽溫室氣體排放總量呈波動(dòng)下降趨勢(shì)。受飼養(yǎng)周期、飼料成本上漲、畜禽疫?。ㄘi藍(lán)耳?。┘澳戏奖?zāi)害等多種因素影響，2007年和2008年散戶平均飼養(yǎng)量顯著下降，強(qiáng)度效應(yīng)抑制作用顯著，分別為-845.23萬t和-731.03萬t，實(shí)現(xiàn)了830.70萬t畜禽溫室氣體的減排。隨后國(guó)家出臺(tái)了一系列支持畜禽轉(zhuǎn)型發(fā)展的政策，中國(guó)畜禽發(fā)展方式在逐年轉(zhuǎn)變，到2013年畜禽溫室氣體排放總量為3 542.48萬t，減少了686.02萬t。該時(shí)期強(qiáng)度效應(yīng)對(duì)溫室氣體排放抑制作用最大，為-1 933.07萬t，其次是勞動(dòng)力效應(yīng)和結(jié)構(gòu)效應(yīng)，分別為-255.96萬t和-133.83萬t；而經(jīng)濟(jì)效應(yīng)促進(jìn)作用顯著，為1 636.84萬t。

總體來看，1991-2013年，經(jīng)濟(jì)效應(yīng)對(duì)畜禽溫室氣體排放促進(jìn)作用最大，為4 692.93萬t；而強(qiáng)度效應(yīng)抑制作用最大，為-2 701.36萬t，其次是勞動(dòng)力效應(yīng)和結(jié)構(gòu)效應(yīng)，分別為-771.85萬t和-424.06萬t。

度呈顯著的波動(dòng)性（見圖2）。從強(qiáng)度效應(yīng)累計(jì)貢獻(xiàn)值演變趨勢(shì)來看，該效應(yīng)對(duì)抑制畜禽溫室氣體排放的貢獻(xiàn)呈倒“U”，且近幾年其抑制作用呈增強(qiáng)趨勢(shì)。1991-1997年，在國(guó)家宏觀調(diào)控和環(huán)境治理影響下，強(qiáng)度效應(yīng)抑制作用不斷加強(qiáng)，累計(jì)減少了1 391.00萬t溫室氣體；1998-2006年，受國(guó)際環(huán)境、高致病性禽流感以及國(guó)內(nèi)農(nóng)業(yè)政策支持乏力等因素影響，規(guī)模化畜禽養(yǎng)殖進(jìn)程緩慢[18]，強(qiáng)度效應(yīng)抑制作用放緩；2007-2013年，隨著畜禽業(yè)以散養(yǎng)模式為主向現(xiàn)代養(yǎng)殖模式（專業(yè)戶模式和規(guī)模化模式）轉(zhuǎn)變，畜禽規(guī)模化養(yǎng)殖推進(jìn)為溫室氣體排放的實(shí)施提供可能[7]，強(qiáng)度效應(yīng)抑制作用呈增強(qiáng)趨勢(shì)，該時(shí)期累計(jì)實(shí)現(xiàn)1 933.07萬t畜禽溫室氣體的減排，占其總效應(yīng)的281%。

勞動(dòng)力效應(yīng)是僅次于強(qiáng)度效應(yīng)，是抑制畜禽溫室氣體排放的另一重要因素。該效應(yīng)累計(jì)貢獻(xiàn)值呈波動(dòng)下降趨勢(shì)，抑制作用越來越明顯。隨著城鎮(zhèn)化和工業(yè)化的深入推進(jìn)，農(nóng)業(yè)比較效益顯著降低，農(nóng)業(yè)勞動(dòng)力不斷轉(zhuǎn)移到非農(nóng)產(chǎn)業(yè)，農(nóng)業(yè)勞動(dòng)力減少導(dǎo)致散養(yǎng)戶大量退出，為畜禽規(guī)?；B(yǎng)殖提供可能；此外，伴隨著畜禽養(yǎng)殖的規(guī)?；l(fā)展和管理模式的不斷創(chuàng)新，對(duì)從事畜禽勞動(dòng)力的素質(zhì)有更高要求，進(jìn)而導(dǎo)致轉(zhuǎn)移更多的畜禽從業(yè)勞動(dòng)力，單位勞動(dòng)力產(chǎn)出大大增加，促進(jìn)了畜禽溫室氣體的減排。1991-2013年，勞動(dòng)力效應(yīng)實(shí)現(xiàn)了771.85萬t畜禽溫室氣體的減排。

結(jié)構(gòu)效應(yīng)累計(jì)貢獻(xiàn)大致呈現(xiàn)低水平徘徊再高水平徘徊再波動(dòng)下降階段性特征，對(duì)畜禽溫室氣體排放的抑制作用也越來越明顯。1991-1997年，結(jié)構(gòu)效應(yīng)對(duì)畜禽溫室氣體排放累計(jì)貢獻(xiàn)處于低水平，年均累計(jì)貢獻(xiàn)為-54.35萬t；1998-2003年，1998年發(fā)生的長(zhǎng)江全流域特大洪災(zāi)，西南地區(qū)、長(zhǎng)江中下游地區(qū)畜禽養(yǎng)殖遭受巨大破壞，全國(guó)畜牧業(yè)產(chǎn)值占農(nóng)業(yè)總產(chǎn)值較1997年下降了2.28%，結(jié)構(gòu)效應(yīng)累計(jì)凈貢獻(xiàn)為-290萬t，隨后幾年受農(nóng)業(yè)結(jié)構(gòu)調(diào)整的影響，畜禽發(fā)展緩慢，結(jié)構(gòu)效應(yīng)累計(jì)貢獻(xiàn)處于較高水平，年均為-269.24萬t；2004-2013年，結(jié)構(gòu)效應(yīng)的抑制作用越來越明顯，但波動(dòng)性較大。主要是因?yàn)?，一是伴隨著農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整，畜牧業(yè)產(chǎn)值占農(nóng)業(yè)總產(chǎn)值由2004年2471%下降到2013年22.10%，下降了2.61%；二是城鎮(zhèn)居民日益增長(zhǎng)的畜禽產(chǎn)品消費(fèi)，畜牧業(yè)在農(nóng)業(yè)結(jié)構(gòu)中的地位進(jìn)一步提升。在這雙重影響下，該時(shí)期結(jié)構(gòu)效應(yīng)的抑制作用波動(dòng)較大。

經(jīng)濟(jì)效應(yīng)累計(jì)貢獻(xiàn)總體上經(jīng)歷了先快速上升再緩慢下降再逐步上升的變化趨勢(shì)。1991-1996年，市場(chǎng)化改革取得重大進(jìn)步，農(nóng)業(yè)得到了快速發(fā)展，經(jīng)濟(jì)效應(yīng)累計(jì)貢獻(xiàn)快速上升，增加了2 254.88萬t畜禽溫室氣體；1997-2000年，受亞洲金融危機(jī)、通貨緊縮及自然災(zāi)害等因素影響，農(nóng)業(yè)發(fā)展外部環(huán)境不佳，經(jīng)濟(jì)效應(yīng)累計(jì)貢獻(xiàn)緩慢下降，減少了502.53萬t畜禽溫室氣體。2001-2013年，經(jīng)濟(jì)效應(yīng)累計(jì)貢獻(xiàn)逐步上升，基本呈指數(shù)增長(zhǎng)的趨勢(shì)，增加了 2 940.57萬t畜禽溫室氣體。主要是因?yàn)?，隨著經(jīng)濟(jì)增長(zhǎng)和人均收入穩(wěn)定提高，城鄉(xiāng)居民膳食結(jié)構(gòu)發(fā)生變化，對(duì)動(dòng)物性食品的消費(fèi)需求不斷增加，從而帶動(dòng)畜牧業(yè)的發(fā)展，畜禽溫室氣體排放不斷增加。由此可見，未來一段時(shí)間內(nèi)，伴隨經(jīng)濟(jì)繼續(xù)平穩(wěn)發(fā)展和城鄉(xiāng)居民收入倍增計(jì)劃的實(shí)施并得到實(shí)現(xiàn)，經(jīng)濟(jì)效應(yīng)依然是導(dǎo)致畜禽溫室氣體排放的最主要因素。

2.2中國(guó)畜禽溫室氣體排放的空間分異

2.2.1畜禽溫室氣體排放的空間比較

由于中國(guó)各?。▍^(qū)、市）資源稟賦差異及畜牧業(yè)結(jié)構(gòu)不同，畜禽溫室氣體排放呈現(xiàn)不同的空間差異，受篇幅限制，本文只列出部分年份畜禽溫室氣體排放位居前10位的?。▍^(qū)、市）（表3）。

從表3可以看出，1991-2013年，畜禽溫室氣體排放大?。▍^(qū)、市）沒有顯著變化，排名前10位?。▍^(qū)、市）畜禽溫室氣體排放量占全國(guó)排放總量的比重約為57%-60%，說明中國(guó)畜禽溫室氣體排放的區(qū)域集中度較高。其中，四川和河南一直占據(jù)中國(guó)畜禽溫室氣體排放前三名，對(duì)畜禽溫室氣體排放貢獻(xiàn)最大。山東、云南和內(nèi)蒙古等?。▍^(qū)、市）的畜禽溫室氣體排放也一直靠前。

2.2.2畜禽溫室氣體排放各效應(yīng)的空間差異

從1991-2013年中國(guó)省域強(qiáng)度效應(yīng)來看（表4），除天津強(qiáng)度效應(yīng)對(duì)畜禽溫室氣體排放起促進(jìn)作用外，各?。▍^(qū)、市）均起到抑制作用。其中，四川、青海和云南規(guī)模化養(yǎng)殖處于發(fā)展階段[18]，強(qiáng)度效應(yīng)提升空間大，從而表現(xiàn)出對(duì)畜禽溫室氣體排放抑制作用顯著，分別為-279.56萬 t、-221.94萬 t和-212.59萬 t。除北京、上海、海南和寧夏因行政區(qū)劃原因，強(qiáng)度效應(yīng)對(duì)畜禽溫室氣體排放抑制作用較小外，遼寧、吉林和黑龍江規(guī)?；笄蒺B(yǎng)殖程度較高，但缺少對(duì)規(guī)模化養(yǎng)殖的畜禽排泄物處理設(shè)施的改進(jìn)[18]，強(qiáng)度效應(yīng)的抑制作用較小，分別為-17.98萬 t、-25.38萬 t和-27.87萬 t；剩余20個(gè)?。▍^(qū)、市）強(qiáng)度效應(yīng)對(duì)畜禽溫室氣體排放抑制作用介于-200～-30萬 t之間。

從結(jié)構(gòu)效應(yīng)來看，山東、四川和黑龍江屬于糧食主產(chǎn)區(qū)，隨著國(guó)家出臺(tái)了一系列促進(jìn)糧食生產(chǎn)的政策，畜牧業(yè)占農(nóng)業(yè)比重不斷下降，分別下降了43.77%、22.51%和

從經(jīng)濟(jì)效應(yīng)來看，各省（區(qū)、市）經(jīng)濟(jì)效應(yīng)對(duì)畜禽溫室氣體排放均起到促進(jìn)作用，但作用強(qiáng)度有差異。四川、河南、內(nèi)蒙古、山東、云南、湖南和河北畜禽溫室氣體排放位居全國(guó)前10位（見表3），屬于畜牧業(yè)大省，但畜禽養(yǎng)殖方式仍以傳統(tǒng)成分占主導(dǎo)，高投入、高排放發(fā)展模式依舊普遍存在，經(jīng)濟(jì)效應(yīng)促進(jìn)作用較大，分別為612.98萬 t、313.64萬 t、271.28萬 t、269.47萬 t、234.54萬 t、220.69萬 t和220.20萬 t；而天津、上海和北京經(jīng)濟(jì)發(fā)展水平相對(duì)較高，但土地面積小，用于養(yǎng)殖空間有限，畜禽養(yǎng)殖方式向集約化、標(biāo)準(zhǔn)化轉(zhuǎn)變[12] ，經(jīng)濟(jì)效應(yīng)促進(jìn)作用較小，分別為10.18萬 t、11.88萬 t和13.97萬 t；海南促進(jìn)作用也較小，為1289萬 t；剩余19個(gè)省（區(qū)、市）對(duì)畜禽溫室氣體排放促進(jìn)作用介于60-200萬 t之間。

從勞動(dòng)力效應(yīng)來看，新疆、黑龍江和內(nèi)蒙古作為全國(guó)畜禽產(chǎn)品的主要來源地，畜禽產(chǎn)品又是勞動(dòng)密集型產(chǎn)品，為滿足日益增加的畜禽產(chǎn)品需求，勞動(dòng)力投入不斷增加，分別增加了172.84萬人、182.7萬人和49.92萬人，勞動(dòng)力效應(yīng)對(duì)畜禽溫室氣體排放促進(jìn)作用顯著，分別為7291萬 t、3113萬 t和1882萬 t；西藏、云南、海南、遼寧、吉林和山西對(duì)畜禽溫室氣體排放促進(jìn)作用介于0-10萬 t之間。四川、湖北、江蘇和山東經(jīng)濟(jì)發(fā)展水平較高，非農(nóng)就業(yè)機(jī)會(huì)多，畜禽養(yǎng)殖比較效益低，勞動(dòng)力大量流出，造成散養(yǎng)戶空欄或轉(zhuǎn)產(chǎn)，為規(guī)?；笄蒺B(yǎng)殖提供了可能，勞動(dòng)力效應(yīng)抑制作用顯著，分別為-17055萬 t、-5610萬 t、-5294萬 t和-4686萬 t；剩余17個(gè)?。▍^(qū)、市）對(duì)畜禽溫室氣體排放抑制作用介于-40-0萬 t之間。

3結(jié)論與討論

本文基于LMDI模型系統(tǒng)分析了1991-2013年中國(guó)畜禽溫室氣體排放時(shí)空變化及其因素貢獻(xiàn)，揭示了強(qiáng)度效應(yīng)、結(jié)構(gòu)效應(yīng)、經(jīng)濟(jì)效應(yīng)和勞動(dòng)力效應(yīng)對(duì)畜禽溫室氣體總效應(yīng)的貢獻(xiàn)，并識(shí)別了不同時(shí)段以及省域畜禽溫室氣體排放量變化的顯著性貢獻(xiàn)因素。結(jié)果表明：

（1）從時(shí)間維度來看，1991-2013年，中國(guó)畜禽溫室氣體排放經(jīng)歷了先快速上升后穩(wěn)定上升再波動(dòng)下降的變化特征，總體呈上升趨勢(shì)。經(jīng)濟(jì)效應(yīng)對(duì)畜禽溫室氣體排放表41991-2013年中國(guó)省域畜禽溫室氣體排放影響因素分解

效應(yīng)和結(jié)構(gòu)效應(yīng)。期間，經(jīng)濟(jì)效應(yīng)促進(jìn)作用的累計(jì)貢獻(xiàn)呈指數(shù)增長(zhǎng)，而強(qiáng)度效應(yīng)抑制作用的累計(jì)貢獻(xiàn)呈倒“U”，是近幾年畜禽溫室氣體增長(zhǎng)趨勢(shì)有所減緩的主要原因，勞動(dòng)力效應(yīng)和結(jié)構(gòu)效應(yīng)抑制作用不斷加強(qiáng)。

（2）從空間維度來看，中國(guó)畜禽溫室氣體排放的區(qū)域集中度較高，四川、河南、山東、云南和內(nèi)蒙古等?。▍^(qū)、市）畜禽溫室氣體排放一直位居全國(guó)前列。省域各效應(yīng)作用方向和程度差異顯著，四川、青海和云南強(qiáng)度效應(yīng)抑制作用較大，遼寧、吉林和黑龍江抑制作用較?。簧綎|、四川和黑龍江結(jié)構(gòu)效應(yīng)抑制作用顯著，新疆和青海促進(jìn)作用明顯；四川、河南、內(nèi)蒙古、山東、云南、湖南和河北經(jīng)濟(jì)效應(yīng)促進(jìn)作用較大，天津、上海、海南和北京促進(jìn)作用較小；四川、湖北、江蘇和山東勞動(dòng)力效應(yīng)抑制作用顯著，新疆、黑龍江和內(nèi)蒙古促進(jìn)作用明顯。

強(qiáng)度效應(yīng)、結(jié)構(gòu)效應(yīng)、經(jīng)濟(jì)效應(yīng)和勞動(dòng)力效應(yīng)空間上的疊加，形成了畜禽溫室氣體排放總效應(yīng)的空間差異。未來中國(guó)畜禽溫室氣體減排的空間發(fā)展策略有以下幾點(diǎn)：①四川、青海和云南等?。▍^(qū)、市）提高畜禽養(yǎng)殖的規(guī)?；⒓s化和標(biāo)準(zhǔn)化，在減少散戶養(yǎng)殖方式同時(shí)降低單位畜禽溫室氣體排放水平，有效提升畜禽養(yǎng)殖產(chǎn)出效率；遼寧、吉林和黑龍江等?。▍^(qū)、市）應(yīng)制定特定性綜合措施，強(qiáng)化畜禽糞便清潔處理技術(shù)的研發(fā)與應(yīng)用。②新疆、青海、云南、陜西和江西等省（區(qū)、市）應(yīng)充分發(fā)揮資源稟賦優(yōu)勢(shì)，優(yōu)化農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)，實(shí)行農(nóng)牧業(yè)有機(jī)結(jié)合型畜牧業(yè)。③四川、河南、內(nèi)蒙古、山東、云南、湖南和河北等?。▍^(qū)、市）要切實(shí)轉(zhuǎn)變農(nóng)業(yè)生產(chǎn)方式，加快推進(jìn)低碳農(nóng)業(yè)發(fā)展，實(shí)現(xiàn)農(nóng)業(yè)生產(chǎn)中經(jīng)濟(jì)、社會(huì)、生態(tài)效益三者統(tǒng)籌兼顧，促進(jìn)畜牧經(jīng)濟(jì)與氣候資源環(huán)境的全面協(xié)調(diào)可持續(xù)發(fā)展。④新疆、黑龍江和內(nèi)蒙古等?。▍^(qū)、市）草地資源豐富、奶牛業(yè)較為發(fā)達(dá)，因此，積極發(fā)展飼料加工業(yè)和牛奶加工業(yè)，推動(dòng)農(nóng)業(yè)勞動(dòng)力轉(zhuǎn)移。

（編輯：尹建中）

參考文獻(xiàn)（References）

[1]FAO. Livestock long shadow[R]. 2006： 97-110.

[2]GOODLAND R， ANHANG J. Livestock and climate change[J]. World watch， 2009， 22（6）： 10-19.

[3]OLESEN J E， SCHELDE K， WEISKE A. Modelling greenhouse gas emissions from european conventional and organic dairy farms[J]. Agrriculture， ecosystems and environment， 2006， 112： 207-220.

[4]IPCC. Greenhouse gas inventory reference manual： Revised 1996 IPCC guidelines for national greenhouse gas inventories[R]. 1996.

[5]董紅敏， 林而達(dá)， 楊其長(zhǎng). 中國(guó)反芻動(dòng)物甲烷排放量的初步估算及減緩技術(shù)[J]. 農(nóng)村生態(tài)環(huán)境， 1995， 11（3）：4-7.[DONG Hongmin， LIN Erda， YANG Qichang. Methane emitted from ruminants in China and mitigation technologies[J]. Rural ecoenvironment， 1995， 11（3）：4-7.]

[6]ZHOU J B， JIANG M M， CHEN G Q. Estimation of methane and nitrous oxide emission from livestock and poltey in China during 1949-2003[J]. Energy policy， 2007， 35：3759-3767.

[7]胡向東， 王濟(jì)民. 中國(guó)畜禽溫室氣體排放量估算[J]. 農(nóng)業(yè)工程學(xué)報(bào)， 2010， 26（10）：247-252.[HU Xiangdong， WANG Jimin. Estimation of livestock greenhouse gases discharge in China[J]. Transactions of the CSAE， 2010， 26（10）： 247-252.]

[8]閔繼勝， 胡浩. 中國(guó)農(nóng)業(yè)生產(chǎn)溫室氣體排放量的測(cè)算[J]. 中國(guó)人口·資源與環(huán)境， 2012， 22（7）： 21-27.[MIN Jisheng， HU Hao. Caculation of greenhouse gases emission from agricultural production in China[J]. China population， resources and environment， 2012， 22（7）： 21-27.]

[9]尚杰， 楊果， 于法穩(wěn). 中國(guó)農(nóng)業(yè)溫室氣體排放量測(cè)算及影響因素研究[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)， 2015， 23（3）： 354-364.[SHANG Jie， YANG Guo， YU Fawen. Agricultural greenhouse gases emissions and influencing factors in China[J]. Chinese journal of ecoagriculture， 2015， 23（3）： 354-364.]

[10]譚秋成. 中國(guó)農(nóng)業(yè)溫室氣體排放： 現(xiàn)狀及挑戰(zhàn)[J]. 中國(guó)人口·資源與環(huán)境， 2011， 21（10）： 69-75.[TAN Qiucheng. Greenhouse gas emission in Chinas agriculture： situation and challenge[J]. China population， resources and environment， 2011， 21（10）： 69-75.]

[11]陳瑤， 尚杰. 中國(guó)畜牧業(yè)脫鉤分析及影響因素研究[J]. 中國(guó)人口·資源與環(huán)境， 2014， 24（3）： 101-107.[CHEN Yao， SHANG Jie. Disconnect analysis and influence factors of animal husbandry in China[J]. China population， resources and environment， 2014， 24（3）： 101-107.]

[12]劉月仙， 劉娟， 吳文良. 北京地區(qū)畜禽溫室氣體排放的時(shí)空變化特征[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)， 2013， 21（7）： 891-897.[LIU Yuexian， LIU Juan， WU Wenlang. Spatiotemporal dynamics of greenhouse gases emissions from livestock and poultry in Beijing Area during 1978-2009[J]. Chinese journal of ecoagriculture， 2013， 21（7）： 891-897.]

[13]劉培芳， 陳振樓， 許世遠(yuǎn)，等. 長(zhǎng)江三角洲城郊畜禽糞便的污染負(fù)荷及其防治對(duì)策[J]. 長(zhǎng)江流域資源與環(huán)境， 2002， 11（5）：456-460.[LIU Peifang， CHEN Zhenlou， XU Shiyuan， et al. Waste loading and treatment strategies on the excreta of domestic animals in the Yangtze Delta[J]. Resources and environment in the Yangtze Basin， 2002， 11（5）：456-460.]

[14]IPCC. IPCC guidelines for national greenhouse gas inventories Volume 4： agriculture， forestry and other land use[R]. Geneva， Switzerland： IPCC， 2006.

[15]ANG B W， LIU F L. A new energy decomposition method： perfect in decomposition and consistent in aggregation[J]. Energy， 2001， 26： 537-548.

[16]郭朝先. 中國(guó)碳排放因素分解： 基于LMDI分解技術(shù)[J]. 中國(guó)人口·資源與環(huán)境， 2010，20（12）： 4-9.[GUO Chaoxian. Decomposition of Chinas carbon emissions： based on LMDI method[J]. China population， resources and environment， 2010， 20（12）： 4-9.]

[17]ANG B W， LIU F L， CHEW E P. Perfect decomposition techniques in energy and environmental analysis[J]. Energy policy， 2003， 31： 1561-1566.

[18]付強(qiáng)， 諸云強(qiáng)， 楊紅新， 等. 2002-2009年中國(guó)規(guī)?；笄蒺B(yǎng)殖量區(qū)域差異及政策建議[J].農(nóng)業(yè)工程學(xué)報(bào)， 2012， 18（9）： 185-191.[FU Qiang， ZHU Yunqiang， YANG Hongxin， et al. Regional differences of livestock and poultry breeding output in scale and policy suggestions during 2002-2009 in China[J]. Transactions of the Chinese Society of Agricultural Engineering， 2012， 18（9）： 185-191.]