白妙琴 馬文江 彭星琳 阿迪力·艾合買(mǎi)提 馬 瓊,2#
(1.塔里木大學(xué)經(jīng)濟(jì)與管理學(xué)院,新疆 阿拉爾 843300;2.新疆南疆經(jīng)濟(jì)社會(huì)發(fā)展研究院,新疆 阿拉爾 843300)
以CO2為代表的溫室氣體排放導(dǎo)致的全球氣候變化已經(jīng)是全世界面臨的一項(xiàng)嚴(yán)峻挑戰(zhàn),如何有效地降低溫室氣體的排放量是亟待解決的重要課題。根據(jù)聯(lián)合國(guó)糧食及農(nóng)業(yè)組織(FAO)估計(jì),中國(guó)畜牧業(yè)產(chǎn)生的CO2排放約占全球CO2排放總量的18%[1]。因此,精確掌握畜牧業(yè)碳排放的狀況和特征,對(duì)于實(shí)現(xiàn)畜牧業(yè)的綠色、低碳發(fā)展有著十分重要的作用。
碳排放強(qiáng)度即單位國(guó)內(nèi)生產(chǎn)總值(GDP)的CO2排放量,是衡量一個(gè)國(guó)家或地區(qū)碳排放量與經(jīng)濟(jì)發(fā)展之間關(guān)系的指標(biāo)[2]。對(duì)于發(fā)展中國(guó)家來(lái)說(shuō),以碳排放強(qiáng)度為衡量標(biāo)準(zhǔn)來(lái)量化減排目標(biāo),而不是以總的碳排放量和人均碳排放量為指標(biāo),更有現(xiàn)實(shí)意義。近年來(lái),學(xué)者對(duì)碳排放進(jìn)行了深入研究,研究?jī)?nèi)容包括碳排放強(qiáng)度的測(cè)度、時(shí)空演變特征與規(guī)律(包括區(qū)域差異、空間集聚性和空間相關(guān)性)和驅(qū)動(dòng)因子。1)碳排放量和碳排放強(qiáng)度的測(cè)度。以反芻動(dòng)物為代表的畜禽在腸道發(fā)酵和糞便處理過(guò)程中,產(chǎn)生了大量的CH4和N2O等溫室氣體[3]。目前國(guó)際上對(duì)碳排放的測(cè)算方法主要有4種,一是聯(lián)合國(guó)政府間氣候變化專(zhuān)門(mén)委員會(huì)(IPCC)系數(shù)法;二是生命周期法;三是投入產(chǎn)出法;四是實(shí)測(cè)法。其中,張廣勝等[4]的研究采用生命周期評(píng)價(jià)法構(gòu)建碳排放測(cè)算體系,研究表明中國(guó)農(nóng)業(yè)碳排放總量呈上升趨勢(shì),農(nóng)業(yè)碳排放強(qiáng)度呈下降趨勢(shì)。郭嬌等[5]對(duì)中國(guó)畜牧業(yè)的碳排放量進(jìn)行了評(píng)價(jià),得出新疆、西藏是CH4減排的重點(diǎn)地區(qū),也是畜牧業(yè)碳減排的主要區(qū)域。2)碳排放時(shí)空間演變特征與規(guī)律。以往學(xué)者主要采用泰爾系數(shù)、錫爾系數(shù)、變異系數(shù)和空間自相關(guān)系數(shù)(全局Moran’s I指數(shù)、局部G系數(shù))等來(lái)分析碳排放水平及其空間關(guān)聯(lián)。吳強(qiáng)等[6]67對(duì)我國(guó)畜牧業(yè)碳排放的時(shí)間和空間演變進(jìn)行了分析,結(jié)果表明我國(guó)畜牧業(yè)的碳排放在時(shí)間和空間上存在著顯著的相關(guān)性,并表現(xiàn)出了顯著的區(qū)域性集聚特點(diǎn)。孫耀華等[7]采用泰爾系數(shù)對(duì)我國(guó)東部、中部、西部和東北各省份的碳排放強(qiáng)度進(jìn)行了比較,得出不同省份的碳排放強(qiáng)度存在著一定程度上的差別,并且隨著時(shí)間的推移,各省份之間的差異呈遞減趨勢(shì)。胡玥[8]把全國(guó)分成了8個(gè)主要的經(jīng)濟(jì)圈,并利用泰爾指數(shù)、全局Moran’s I指數(shù)和冷熱點(diǎn)分析等方法,對(duì)全國(guó)碳排放強(qiáng)度的空間演變特征進(jìn)行了研究,結(jié)果表明不同地區(qū)之間碳排放強(qiáng)度有很大的差別,且差異逐漸擴(kuò)大。3)碳排放驅(qū)動(dòng)因素研究。以往的研究中關(guān)于碳排放驅(qū)動(dòng)因素的分析大致集中于經(jīng)濟(jì)、人口、制度與技術(shù)等因素的考量。程葉青等[9]運(yùn)用空間計(jì)量方法,分析出能源結(jié)構(gòu)、產(chǎn)業(yè)結(jié)構(gòu)、城市化水平等因素對(duì)我國(guó)城市群的碳排放強(qiáng)度變化具有顯著影響;顏艷梅等[10]采用夏普里分析法,對(duì)我國(guó)主要省份的碳排放強(qiáng)度進(jìn)行了分析,發(fā)現(xiàn)經(jīng)濟(jì)發(fā)展水平對(duì)碳排放強(qiáng)度差異貢獻(xiàn)最大,且這種影響隨時(shí)間而遞增。周杰琦[11]通過(guò)面板回歸模型分析碳排放強(qiáng)度的影響因素,得出人均資本水平、城鎮(zhèn)化水平、市場(chǎng)化程度、行業(yè)結(jié)構(gòu)等是導(dǎo)致我國(guó)不同地區(qū)之間碳排放強(qiáng)度差異的重要因素。
以往學(xué)者在該領(lǐng)域已經(jīng)開(kāi)展了大量的相關(guān)工作,但尚未開(kāi)展基于市域?qū)用娴男竽翗I(yè)碳排放的空間分布規(guī)律及驅(qū)動(dòng)因子的研究。新疆是全國(guó)4個(gè)主要畜牧業(yè)產(chǎn)區(qū)之一,由于畜牧業(yè)的持續(xù)擴(kuò)張,畜牧業(yè)碳排放量也在持續(xù)上升,使得本來(lái)就很脆弱的生態(tài)系統(tǒng)受到了破壞。為此,本研究選取新疆14個(gè)地(州、市)作為研究對(duì)象,依據(jù)IPCC系數(shù)法和空間自相關(guān)分析等方法,估算新疆2001—2020年畜牧業(yè)的碳排放量,進(jìn)而探討新疆各地(州、市)畜牧業(yè)的碳排放空間集聚特征與驅(qū)動(dòng)因素,以期為實(shí)現(xiàn)新疆畜牧業(yè)的低碳可持續(xù)發(fā)展提供參考。
畜禽消化、糞肥處理等活動(dòng)所釋放的CH4,糞肥還田所釋放的N2O等是我國(guó)畜牧業(yè)的重要碳源[6]67。本研究從新疆畜牧業(yè)生產(chǎn)實(shí)際出發(fā),選取反芻動(dòng)物(牛、羊)和非反芻動(dòng)物(豬、馬、驢)為碳源,參照胡向東等[12]提出的畜牧業(yè)生產(chǎn)碳排放核算法,按照《省級(jí)溫室氣體清單編制指南》中的碳排放因子,測(cè)算新疆不同地區(qū)的畜牧業(yè)碳排放量。主要牲畜品種CH4和N2O的排放系數(shù)見(jiàn)表1。為方便分析,將CH4、N2O統(tǒng)一轉(zhuǎn)換為CO2當(dāng)量,也就是1 t CH4折合25 t CO2,1 t N2O折合298 t CO2。
表1 主要牲畜品種的排放系數(shù)Table 1 Carbon emission coefficients of major livestock breeds kg/(頭·a)
畜牧業(yè)碳排放強(qiáng)度是指單位農(nóng)業(yè)經(jīng)濟(jì)效益所排放的碳,由于反芻動(dòng)物(牛、羊)和非反芻動(dòng)物(豬、馬、驢)在新疆牲畜養(yǎng)殖中占95%以上,故畜牧業(yè)碳排放強(qiáng)度用牛、羊、豬、馬、驢的碳排放總量除以畜牧業(yè)總產(chǎn)值來(lái)計(jì)算。
空間自相關(guān)分為全局空間自相關(guān)和局部空間自相關(guān)。采用全局Moran’s I指數(shù)對(duì)新疆各地(州、市)碳排放強(qiáng)度的整體關(guān)聯(lián)度和空間分布模式進(jìn)行測(cè)算,采用局部Moran’s I指數(shù)揭示各地(州、市)的畜牧業(yè)碳排放強(qiáng)度與鄰近地(州、市)的相似性或相關(guān)性,區(qū)分不同地區(qū)之間的空間集聚與空間孤立特點(diǎn),檢驗(yàn)畜牧業(yè)碳排放強(qiáng)度的空間異質(zhì)性。
采用2002—2021年新疆統(tǒng)計(jì)年鑒、中國(guó)農(nóng)村統(tǒng)計(jì)年鑒、中國(guó)畜牧業(yè)年鑒,對(duì)新疆各地(州、市)2001—2020年的畜牧業(yè)碳排放量和碳排放強(qiáng)度進(jìn)行分析。將新疆14個(gè)地(州、市)按地理位置劃分為北疆、東疆和南疆,北疆包括烏魯木齊市、克拉瑪依市、昌吉回族自治州(以下簡(jiǎn)稱(chēng)昌吉州)、伊犁哈薩克自治州(以下簡(jiǎn)稱(chēng)伊犁州)、塔城地區(qū)、阿勒泰地區(qū)和博爾塔拉蒙古自治州(以下簡(jiǎn)稱(chēng)博州);東疆包括吐魯番市和哈密市;南疆包括巴音郭楞蒙古自治州(以下簡(jiǎn)稱(chēng)巴州)、阿克蘇地區(qū)、克孜勒蘇柯?tīng)柨俗巫灾沃?以下簡(jiǎn)稱(chēng)克州)、喀什地區(qū)與和田地區(qū)。
由圖1可見(jiàn),2001—2020年新疆畜牧業(yè)碳排放量大體呈波動(dòng)變化趨勢(shì),2020年新疆畜牧業(yè)碳排放量為2 379.56 萬(wàn)t,相較于2001年增加了139.01萬(wàn)t。2016年畜牧業(yè)碳排放量達(dá)到峰值,為2 915.76 萬(wàn)t。而2001—2020年新疆畜牧業(yè)碳排放強(qiáng)度整體呈下降趨勢(shì),由2001年的16.43 t/萬(wàn)元下降至2020年的3.22 t/萬(wàn)元,減少了13.21 t/萬(wàn)元,這表明該時(shí)期畜牧業(yè)碳減排有一定的有效性,有必要進(jìn)行系統(tǒng)分析。李娜等[13]基于2000—2020年面板數(shù)據(jù),測(cè)算了新疆14個(gè)地(州、市)畜牧業(yè)碳排放量,發(fā)現(xiàn)研究期內(nèi)碳排放量整體呈波動(dòng)變化趨勢(shì),與本研究結(jié)果大體一致,但文獻(xiàn)[13]未對(duì)碳排放強(qiáng)度進(jìn)行計(jì)算,也未對(duì)其進(jìn)行驅(qū)動(dòng)力分析。
圖1 2001—2020年新疆畜牧業(yè)碳排放量和碳排放強(qiáng)度Fig.1 Carbon emission amount and intensity of animal husbandry in Xinjiang from 2001 to 2020
2001—2020年新疆各地(州、市)的畜牧業(yè)碳排放量與碳排放強(qiáng)度見(jiàn)表2,碳排放量高值區(qū)(大于200萬(wàn) t/a)主要分布在伊犁州、塔城地區(qū)、阿勒泰地區(qū)、阿克蘇地區(qū)和喀什地區(qū);中高值區(qū)(150萬(wàn)~200萬(wàn) t/a)分布在昌吉州與和田地區(qū);中值區(qū)(100萬(wàn)~150萬(wàn) t/a)分布在巴州;中低值區(qū)(50萬(wàn)~100 萬(wàn) t/a)分布在克州;低值區(qū)(小于50萬(wàn) t/a)分布在烏魯木齊市、克拉瑪依市、博州、吐魯番市和哈密市。有2個(gè)地(州、市)的畜牧業(yè)碳排放強(qiáng)度屬于高強(qiáng)度型(7~10 t/萬(wàn)元),分別為阿勒泰地區(qū)和克州;7個(gè)地(州、市)屬于中強(qiáng)度型(4~7 t/萬(wàn)元),主要分布于伊犁州、塔城地區(qū)、博州、巴州、阿克蘇地區(qū)、喀什地區(qū)與和田地區(qū)。5個(gè)地(州、市)屬于低強(qiáng)度型(0~4 t/萬(wàn)元),主要分布于烏魯木齊市、克拉瑪依市、昌吉州、吐魯番市和哈密市??傮w而言,2001—2020年新疆各地(州、市)畜牧業(yè)碳排放量高值區(qū)多分布于北疆,低值區(qū)多分布于東疆;南疆碳排放強(qiáng)度明顯高于北疆、東疆,說(shuō)明南疆畜牧業(yè)生產(chǎn)效率遠(yuǎn)低于其他地區(qū),可能的原因?yàn)槟辖竽翗I(yè)科技水平較低。
表2 2001—2020年新疆各地(州、市)畜牧業(yè)碳排放量和碳排放強(qiáng)度Table 2 Carbon emission amount and intensity of animal husbandry in various regions (prefectures,cities) of Xinjiang from 2001 to 2020
采用全局Moran’s I指數(shù)對(duì)新疆各地(州、市)的畜牧業(yè)碳排放強(qiáng)度的動(dòng)態(tài)演變特征進(jìn)行刻畫(huà),并對(duì)其時(shí)空相關(guān)性進(jìn)行檢驗(yàn)。表3顯示了研究期間的全局Moran’s I指數(shù)變化,除了2017—2020年,其他年份都是正值,大部分年份通過(guò)了顯著性檢驗(yàn)(p<0.1),這說(shuō)明新疆各地(州、市)之間的畜牧業(yè)碳排放強(qiáng)度存在著明顯的空間集聚和空間正相關(guān)現(xiàn)象。
表3 2001—2020年全局Moran’s I指數(shù)Table 3 Global Moran’s I index from 2001 to 2020
由表3可知,在全局Moran’s I指數(shù)檢驗(yàn)中未通過(guò)顯著性檢驗(yàn)(p<0.1)的年份均處于2010—2020年,為了更好地揭示2010—2020年間新疆畜牧業(yè)發(fā)展過(guò)程中的碳排放空間集聚特征,需要驗(yàn)證其局部空間自相關(guān)性,為此,選擇2010、2020年新疆畜牧業(yè)碳排放強(qiáng)度局部Moran’s I指數(shù)散點(diǎn)圖進(jìn)行分析,位于第一象限的地(州、市)數(shù)目(均為3個(gè))一致,即“高-高”集聚態(tài)勢(shì)的地(州、市)在數(shù)量上相同,位于第二象限和第四象限的地州市分別為由3個(gè)演變?yōu)?個(gè)和2個(gè)演變?yōu)?個(gè),意味著2010—2020年形成“低-高”集聚和“高-低”集聚態(tài)勢(shì)的地(州、市)數(shù)量增多;位于第三象限的地(州、市)由6個(gè)演變?yōu)?個(gè),意味著2010—2020年形成“低-低”集聚的地(州、市)數(shù)量減少。
借鑒唐志忠等[14]和楊紅娟等[15]的研究,并立足于新疆畜牧業(yè)的現(xiàn)實(shí)研究新疆畜牧業(yè)碳排放的影響因素。將畜牧業(yè)碳排放總量和碳排放強(qiáng)度分別作為被解釋變量,選取經(jīng)濟(jì)發(fā)展水平、城鎮(zhèn)化水平、畜牧業(yè)發(fā)展水平、農(nóng)業(yè)機(jī)械化水平和工業(yè)化發(fā)展水平5個(gè)因素考察它們對(duì)新疆畜牧業(yè)碳排放量和碳排放強(qiáng)度的驅(qū)動(dòng)性。其中:經(jīng)濟(jì)發(fā)展水平采用地區(qū)人均GDP來(lái)衡量;畜牧業(yè)發(fā)展水平采用年度畜牧業(yè)生產(chǎn)總值來(lái)衡量;城鎮(zhèn)化水平采用城鎮(zhèn)人口占總?cè)丝诘谋壤齺?lái)衡量;農(nóng)業(yè)機(jī)械化水平和工業(yè)化發(fā)展水平選擇2001—2020年新疆14個(gè)地(州、市)面板數(shù)據(jù)。構(gòu)建時(shí)空雙向固定的固定效應(yīng)模型(見(jiàn)式(1)),實(shí)證分析結(jié)果見(jiàn)表4。
表4 實(shí)證分析結(jié)果1)Table 4 Empirical analysis results
lnCit=α+β1lnPit+β2lnAit+β3lnUit+
β4lnMit+β5lnDit+μit
(1)
式中:Cit為t時(shí)區(qū)域i的畜牧業(yè)碳排放量或碳排放強(qiáng)度;Pit、Ait、Uit、Mit和Dit分別為t時(shí)區(qū)域i的經(jīng)濟(jì)發(fā)展水平、畜牧業(yè)發(fā)展水平、城鎮(zhèn)化水平、農(nóng)業(yè)機(jī)械化水平和工業(yè)化發(fā)展水平面板數(shù)據(jù);α為常數(shù)項(xiàng);β1、β2、β3、β4、β5為待估系數(shù);μit為時(shí)空固定效應(yīng)系數(shù)。
由表4可知,經(jīng)濟(jì)發(fā)展水平、畜牧業(yè)發(fā)展水平、城鎮(zhèn)化水平、農(nóng)業(yè)機(jī)械化水平和工業(yè)化發(fā)展水平均對(duì)新疆畜牧業(yè)碳排放強(qiáng)度具有顯著影響。經(jīng)濟(jì)發(fā)展水平每增加1%,畜牧業(yè)碳排放強(qiáng)度減少0.534%;畜牧業(yè)發(fā)展水平每增加1%,畜牧業(yè)碳排放強(qiáng)度減少0.470%;城鎮(zhèn)化水平每增長(zhǎng)1%,畜牧業(yè)碳排放強(qiáng)度將增加0.126%;農(nóng)業(yè)機(jī)械化水平每增加1%,畜牧業(yè)碳排放強(qiáng)度減少0.152%;工業(yè)化發(fā)展水平每增加1%,畜牧業(yè)碳排放強(qiáng)度增加0.264%;經(jīng)濟(jì)發(fā)展水平、畜牧業(yè)發(fā)展水平和農(nóng)業(yè)機(jī)械化發(fā)展水平對(duì)畜牧業(yè)碳排放強(qiáng)度呈現(xiàn)出負(fù)效應(yīng),這是由于隨著新疆各地區(qū)經(jīng)濟(jì)水平的發(fā)展,生產(chǎn)手段的進(jìn)步以及機(jī)械設(shè)備的改進(jìn)降低了畜牧業(yè)碳排放強(qiáng)度,從而提升了畜牧業(yè)的生產(chǎn)效率。城鎮(zhèn)化水平和工業(yè)化發(fā)展水平對(duì)畜牧業(yè)碳排放強(qiáng)度呈現(xiàn)出正效應(yīng),很可能是隨著城市化程度的加深,農(nóng)村地區(qū)人才的外流速度加快,從而制約了畜牧業(yè)生產(chǎn)效率的提升,導(dǎo)致畜牧業(yè)碳排放強(qiáng)度提高;各地(州、市)的工業(yè)化發(fā)展水平的“虹吸”作用造成了畜牧業(yè)生產(chǎn)的專(zhuān)業(yè)人才大量流失,工業(yè)對(duì)農(nóng)業(yè)的“反哺”作用不明顯,因此導(dǎo)致了畜牧業(yè)碳排放強(qiáng)度的提升。
1) 從總體上來(lái)看,2001—2020年,新疆畜牧業(yè)碳排放量處于波動(dòng)變化趨勢(shì);而畜牧業(yè)碳排放強(qiáng)度呈下降趨勢(shì)。
2) 從各地(州、市)來(lái)看,2001—2020年,新疆各地(州、市)畜牧業(yè)碳排放量空間分布大致為北疆>南疆>東疆;而畜牧業(yè)碳排放強(qiáng)度也表現(xiàn)出明顯的
空間差異,其中克州的畜牧業(yè)碳排放強(qiáng)度居于首位。
3) 從空間集聚特征來(lái)看,新疆各地(州、市)畜牧業(yè)碳排放強(qiáng)度呈現(xiàn)明顯的空間集聚和空間正相關(guān)現(xiàn)象。
4) 從碳排放強(qiáng)度影響因素分析結(jié)果來(lái)看,城鎮(zhèn)化水平和工業(yè)化水平提升對(duì)新疆各地(州、市)畜牧業(yè)碳排放強(qiáng)度降低具有顯著的阻礙作用,經(jīng)濟(jì)發(fā)展水平、畜牧業(yè)發(fā)展水平和農(nóng)業(yè)機(jī)械化水平提升對(duì)新疆各地(州、市)畜牧業(yè)碳排放強(qiáng)度降低具有顯著的推動(dòng)作用。