羅芳 夏慶利

摘要：根據(jù)湖北省1980～2009年樣本數(shù)據(jù)，將農(nóng)業(yè)總產(chǎn)值指數(shù)對農(nóng)業(yè)機械總動力、化肥施用量、大牲畜存欄頭數(shù)、生豬存欄頭數(shù)和造林面積等進行回歸分析，其中，造林面積用于表征“低碳”，其余解釋變量代表“高碳”，以此甄別湖北省低碳農(nóng)業(yè)發(fā)展現(xiàn)狀。結(jié)果表明，僅農(nóng)業(yè)機械總動力和化肥施用量對農(nóng)業(yè)總產(chǎn)值指數(shù)有顯著影響，說明湖北省農(nóng)業(yè)仍處于“高碳”發(fā)展路徑，低碳農(nóng)業(yè)有著廣闊的發(fā)展空間。計算了兩個主要因素對農(nóng)業(yè)總產(chǎn)值的貢獻率后，提出了湖北省發(fā)展低碳農(nóng)業(yè)的前景方略。

關(guān)鍵詞：湖北?。坏吞嫁r(nóng)業(yè)；現(xiàn)狀；前景

中圖分類號：F830 文獻標識碼：A 文章編號：0439-8114（2013）06-1460-05

自工業(yè)革命以來，全球平均溫度年均增加約0.97 ℃[1]。據(jù)2007年聯(lián)合國政府間氣候變化專門委員會（IPCC）第四次評估報告預(yù)測，到本世紀末，世界平均溫度將在當(dāng)前17 ℃的基礎(chǔ)上增溫1.1～6.4 ℃[1，2]。地球持續(xù)變暖等氣候變化導(dǎo)致了重大自然災(zāi)害，如持續(xù)干旱、超強颶風(fēng)、暴雪和猛烈的冰雹襲擊等頻繁發(fā)生。引起氣候變化的原因是多方面的，其中最主要的原因之一是溫室氣體的過量排放所導(dǎo)致的溫室效應(yīng)。CO2、CH4和N2O等是與農(nóng)業(yè)密切關(guān)聯(lián)的主要溫室氣體。農(nóng)業(yè)是極度依賴自然生態(tài)環(huán)境的產(chǎn)業(yè)，在極端氣候面前顯得尤為脆弱。對于碳排放而言，農(nóng)業(yè)具有雙重效應(yīng)，既是碳源又是碳匯，因此，發(fā)展低碳農(nóng)業(yè)可以通過碳減排和固碳兩種方式實現(xiàn)。低碳農(nóng)業(yè)是低碳經(jīng)濟在農(nóng)業(yè)發(fā)展中的實現(xiàn)形式。與生態(tài)農(nóng)業(yè)和特色農(nóng)業(yè)不同，低碳農(nóng)業(yè)是為維護全球生態(tài)安全、改善全球氣候條件而在農(nóng)業(yè)領(lǐng)域推廣節(jié)能減排技術(shù)、固碳技術(shù)、開發(fā)生物質(zhì)能源和可再生能源的農(nóng)業(yè)，是以低能耗、低排放、低污染為特征，具備農(nóng)業(yè)生產(chǎn)、安全保障、氣候調(diào)節(jié)、生態(tài)涵養(yǎng)和農(nóng)村金融等多元功能的新型農(nóng)業(yè)[3]。

湖北是農(nóng)業(yè)大省，根據(jù)統(tǒng)計數(shù)據(jù)，湖北省2009年的農(nóng)林牧漁總產(chǎn)值為2 985.2億元，居全國第八位[4]。農(nóng)業(yè)以耕作業(yè)為主，糧食生產(chǎn)居首要地位，是中國重要的糧、棉、油和生豬等的生產(chǎn)基地，糧棉油產(chǎn)量和生豬存欄頭數(shù)分別為2 309.1萬t、48.1萬t、314.2萬t和2 546.1萬頭，在全國的排名都在前十以內(nèi)。為了改善自然生態(tài)環(huán)境，也為了實現(xiàn)中國政府提出的“到2020年中國單位GDP的CO2排放比2005年下降40%～45%，并作為約束性指標納入國民經(jīng)濟和社會發(fā)展中長期規(guī)劃”[5]的目標，湖北省發(fā)展低碳農(nóng)業(yè)刻不容緩。

有的學(xué)者采用實證分析的方法來甄別低碳農(nóng)業(yè)發(fā)展現(xiàn)狀的代表性研究，如漆雁斌等[1]利用中國1985～2008年的樣本數(shù)據(jù)進行農(nóng)業(yè)總產(chǎn)值對化肥施用量、機械總動力和農(nóng)作物總播種面積等的回歸分析，發(fā)現(xiàn)化肥施用量對農(nóng)業(yè)產(chǎn)值有顯著影響，是制約低碳農(nóng)業(yè)發(fā)展的重要因素；馬倫姣[6]根據(jù)湖北省1980～2008年數(shù)據(jù)得出化肥施用量和農(nóng)業(yè)機械總動力既是明顯影響農(nóng)業(yè)總產(chǎn)值的因素，亦是低碳農(nóng)業(yè)發(fā)展的制約因素的結(jié)論。

與上述兩個研究的不同之處是，本研究不僅要甄別湖北省農(nóng)業(yè)的“高碳”依存度，即農(nóng)業(yè)總產(chǎn)值指數(shù)對農(nóng)業(yè)機械總動力、化肥施用量、大牲畜存欄頭數(shù)和生豬存欄頭數(shù)等回歸，而且要考量其對“低碳”的依存度，即將造林面積作為解釋變量予以分析；并計算了顯著影響農(nóng)業(yè)總產(chǎn)值的因素的貢獻率。

1 低碳農(nóng)業(yè)發(fā)展的分析框架

“低碳”農(nóng)業(yè)是相對于現(xiàn)代的“高碳”農(nóng)業(yè)而言的[6]?，F(xiàn)代農(nóng)業(yè)高度依賴于石油、煤和天然氣等能源，具有高產(chǎn)、高效、省時、省力、不施糞肥和經(jīng)濟收益高等特點，對提高農(nóng)業(yè)生產(chǎn)效率和農(nóng)產(chǎn)品產(chǎn)量，以及解決世界人口激增而引起糧食需求猛增這一尖銳問題發(fā)揮了重要的作用。但同時也應(yīng)看到，現(xiàn)代農(nóng)業(yè)因大量使用農(nóng)業(yè)機械、化肥、塑料薄膜、柴油和農(nóng)藥等，使其具有高能耗、高污染和高排放的特點。此外，畜禽飼養(yǎng)也是一個制造碳源的過程。據(jù)2006年11月聯(lián)合國糧農(nóng)組織估計[7]，畜禽的呼吸和消化要產(chǎn)生CO2，每年其釋放總量比汽車多18%；部分牲畜在反芻過程中排放CH4，該排放量占人類活動產(chǎn)生CH4的37%，而CH4形成溫室效應(yīng)的能力是CO2的23倍；牲畜產(chǎn)生的NH3占人類活動產(chǎn)生總量的64%，NH3是形成酸雨的重要因素之一。農(nóng)業(yè)除了具有碳源的作用外，它還擁有碳匯功能。植物通過光合作用可以將大氣中的CO2轉(zhuǎn)化為碳水化合物，并以有機碳的形式固定在植物體內(nèi)或土壤中，植樹造林尤其有助于對大量CO2的吸收。

為了探究湖北省農(nóng)業(yè)發(fā)展對“低碳”農(nóng)業(yè)以及“高碳”農(nóng)業(yè)的依存度，選取農(nóng)林牧漁總產(chǎn)值為因變量，農(nóng)業(yè)機械總動力、農(nóng)用物資使用量、牲畜飼養(yǎng)量和造林面積等作為解釋變量，其中，前3項用來刻畫“高碳”農(nóng)業(yè)的屬性，后1項用來表征“低碳”農(nóng)業(yè)的屬性。通過湖北省的樣本數(shù)據(jù)，如果回歸分析的結(jié)果證明“高碳”屬性變量是顯著的，則說明湖北農(nóng)業(yè)仍處于“高碳”農(nóng)業(yè)的發(fā)展路徑上，低碳農(nóng)業(yè)還有巨大的發(fā)展空間；如果“低碳”屬性變量是顯著的，則表明湖北省農(nóng)業(yè)已走上可持續(xù)發(fā)展的良性軌道。

2 數(shù)據(jù)來源與模型選擇

2.1 數(shù)據(jù)來源

1980～2009年的樣本數(shù)據(jù)來自《1981～2010年湖北統(tǒng)計年鑒》（表1）。為了消除價格波動的影響，農(nóng)業(yè)總產(chǎn)值采用農(nóng)林牧漁總產(chǎn)值指數(shù)來代表，并設(shè)初期為100，即1980=100。農(nóng)用物資本應(yīng)包括化肥、塑料薄膜、柴油和農(nóng)藥等，但后面3種物資只有近期的年鑒有統(tǒng)計數(shù)據(jù)，早期的則沒有，而化肥施用量的數(shù)據(jù)很完備，因此，這里僅采用化肥施用折純量的數(shù)據(jù)。牲畜飼養(yǎng)情況用大牲畜存欄頭數(shù)和生豬存欄頭數(shù)來表示。此外，還有農(nóng)業(yè)機械總動力和造林面積數(shù)據(jù)。表1還給出了農(nóng)業(yè)總產(chǎn)值、農(nóng)業(yè)機械總動力和化肥施用量等的增長率，這些數(shù)值是此表中相應(yīng)數(shù)據(jù)通過簡單計算[即（當(dāng)年值-上年值）/上年值]得出，后面計算各因素對農(nóng)業(yè)總產(chǎn)值的貢獻率將用到。

2.2 模型的選擇

因變量為農(nóng)業(yè)即農(nóng)林牧漁業(yè)的總產(chǎn)值指數(shù)，解釋變量為農(nóng)業(yè)機械總動力、化肥施用折純量、造林面積、大牲畜存欄頭數(shù)和生豬存欄頭數(shù)等。柯布-道格拉斯生產(chǎn)函數(shù)能夠較好地表征農(nóng)業(yè)生產(chǎn)情況，根據(jù)需要采用雙對數(shù)模型的形式，記為模型1，即：

lnYt=B1+B2lnMt+B3lnFt+B4lnWt+B5lnLt+B6lnHt+ut

（1）

式中，t表示時間，即年份，Bi（i=1，2，…，6）為截距及偏彈性系數(shù)，Y為農(nóng)業(yè)總產(chǎn)值指數(shù)（1980=100），M為農(nóng)業(yè)機械總動力（萬kw·h），F(xiàn)為化肥施用折純量（萬t），W為造林面積（k·hm2），L為大牲畜存欄頭數(shù)（萬頭），H為生豬存欄頭數(shù)（萬頭），u為誤差項。

預(yù)期各解釋變量對因變量有正的影響，即各偏彈性系數(shù)為正值。

3 回歸結(jié)果分析

3.1 自相關(guān)性的檢驗與消除

回歸利用統(tǒng)計軟件EViews 6.0，用回歸方程（1）擬合表1中的數(shù)據(jù)，采用OLS（普通最小二乘）法，結(jié)果如表2所示。由于樣本是時間序列數(shù)據(jù)，誤差項可能不滿足經(jīng)典線性回歸關(guān)于無自相關(guān)性的基本假定。假如誤差項存在自相關(guān)，那么OLS估計量盡管是線性的和無偏的，但卻不是有效的[8]。

表2中模型1的D-W值等于0.784，接近于0，表明回歸存在嚴重的正的自相關(guān)。迭代法是處理自相關(guān)的標準方法。在EViews中，迭代法就是增加AR（自回歸）項作為解釋變量，即增加AR（1）消除一階自相關(guān)，增加AR（2）消除二階自相關(guān)，增加AR（1）、AR（2）消除一階和二階自相關(guān)，依此類推。因此，首先在模型1中增加解釋變量AR（1），得到模型2。模型2的D-W值等于1.582，說明AR（1）項的加入很大程度地消除了一階自相關(guān)，但還存在微弱的正的自相關(guān)。繼續(xù)增加AR（2）項，得到模型3。模型3的D-W值等于2.028，接近于2，表明一階和二階自相關(guān)被完全消除。

3.2 估計結(jié)果分析

由于模型3的估計量是BLUE，即最優(yōu)線性無偏估計量，下面的分析僅針對該模型而言。模型3的F值為206.336，其概率值為0.000，說明模型的整體回歸效果尚佳。農(nóng)業(yè)機械總動力M、化肥施用折純量F和AR（1）在5%的顯著性水平上顯著，其他解釋變量即使在15%的顯著性水平上都不顯著。農(nóng)業(yè)總產(chǎn)值指數(shù)的機械動力彈性是0.683，即農(nóng)業(yè)機械總動力增加1%，農(nóng)業(yè)總產(chǎn)值指數(shù)平均增加0.683%。農(nóng)業(yè)總產(chǎn)值指數(shù)的化肥彈性是0.472，即化肥施用折純量增加1%，農(nóng)業(yè)總產(chǎn)值指數(shù)平均增加0.472%。此二彈性之和為1.155（≥1），可見湖北省農(nóng)業(yè)生產(chǎn)處于規(guī)模報酬不變或遞增階段。具體是規(guī)模報酬不變或遞增還需進一步的檢驗，不能就1.155看似大于1而推斷規(guī)模報酬遞增，因統(tǒng)計上它可能等于1。因此，這里做更深一層的檢驗如下：

在模型3回歸方程輸出結(jié)果界面上，選擇View（視圖）菜單的子菜單Coefficient Tests（系數(shù)檢驗），再選擇Wald Test，在限制條件的窗口輸入：

C（2）+C（3）=1 （2）

該方程等價于

H0：規(guī)模報酬不變；H1：規(guī)模報酬遞增 （3）

其中，H0為原假設(shè)，H1為備擇假設(shè)。

Wald檢驗的結(jié)果為，F(xiàn)值等于2.404，其概率值等于0.137，χ2值等于2.404，其概率值等于0.121。若選擇α為15%，則拒絕規(guī)模報酬不變的原假設(shè)；若選擇α為10%，則不拒絕原假設(shè)。由此可見，湖北農(nóng)業(yè)正處于規(guī)模報酬遞增和不變的交匯處，有少許規(guī)模擴張的空間。

上述結(jié)果表明，農(nóng)業(yè)機械總動力和化肥施用折純量等高能耗、高污染和高排放的“高碳”農(nóng)業(yè)生產(chǎn)方式支配著農(nóng)業(yè)總產(chǎn)值，而造林面積這一低碳農(nóng)業(yè)生產(chǎn)方式對農(nóng)業(yè)總產(chǎn)值無顯著影響。換言之，湖北農(nóng)業(yè)依循的是一條“高碳”農(nóng)業(yè)，即粗放式經(jīng)營的發(fā)展路徑；這也表明，低碳農(nóng)業(yè)有著廣闊的發(fā)展空間和巨大的潛力，其前景遠大。

4 各因素對農(nóng)業(yè)總產(chǎn)值的貢獻率

5 湖北省低碳農(nóng)業(yè)發(fā)展前景方略

5.1 降低設(shè)施農(nóng)業(yè)的比重

設(shè)施農(nóng)業(yè)利用人工建造的設(shè)施，使傳統(tǒng)農(nóng)業(yè)逐步擺脫自然條件的束縛，走向現(xiàn)代工廠化農(nóng)業(yè)、環(huán)境安全型農(nóng)業(yè)生產(chǎn)、無毒農(nóng)業(yè)的發(fā)展之路，并且打破了傳統(tǒng)農(nóng)業(yè)的季節(jié)性，實現(xiàn)了農(nóng)產(chǎn)品的反季節(jié)上市，進一步滿足了人們多元化、多層次的消費需求。但是，設(shè)施農(nóng)業(yè)過度依賴工業(yè)生產(chǎn)的農(nóng)用薄膜等。中國農(nóng)膜使用量居世界第一位，2007年使用量高達200萬t；2009年湖北農(nóng)用薄膜使用量達到6.12萬t。由于目前使用的絕大部分農(nóng)膜為不可降解塑料，在土壤中自然分解需要200年以上的時間。每年回收的極少，遺留在土壤中的農(nóng)膜將近一半，嚴重影響到土壤的物理性能和結(jié)構(gòu)，可耕性遭到破壞，地力下降；土壤的保水透水能力降低，削弱了農(nóng)田的抗旱能力，還會影響到種子發(fā)芽、長苗及根系生長發(fā)育，最終導(dǎo)致農(nóng)產(chǎn)品產(chǎn)量下降[10]。因此，應(yīng)大幅度地降低設(shè)施農(nóng)業(yè)的比重，或者盡可能地使用可降解農(nóng)膜代替塑料農(nóng)膜。

5.2 根據(jù)農(nóng)產(chǎn)品的碳貢獻率征收碳排放稅

根據(jù)各個農(nóng)業(yè)部門對碳排放的貢獻情況，對農(nóng)產(chǎn)品征稅。例如，反芻動物瘤胃發(fā)酵、畜禽糞便處理過程中產(chǎn)生的CH4，以及糞便還田利用過程中直接或間接的CO2和N2O排放等，是重要的溫室氣體來源[11]。據(jù)測算，牛的糞便釋放出的CO2是人類的80倍[12]；同等分量的由動物反芻所產(chǎn)生的CH4的溫室效應(yīng)是CO2的20余倍[13]。以此為依據(jù)，如果對畜牧養(yǎng)殖征收碳排放稅，那么相應(yīng)的產(chǎn)品（如牛肉與乳制品）的價格將遠遠高于沒有征稅時。通過稅賦杠桿調(diào)節(jié)畜牧業(yè)產(chǎn)量，將顯著地改變農(nóng)業(yè)生產(chǎn)方式以及人們的飲食習(xí)慣，尤其是那些大量排放CH4的農(nóng)產(chǎn)品產(chǎn)量與需求。根據(jù)孫亞男等[14]對規(guī)?；膛鰷厥覛怏w排放量的估算，平均每頭奶牛排放CO2當(dāng)量為4.9 t/年，湖北省2009年牛存欄數(shù)量為334.01萬頭，那么牛排放的CO2當(dāng)量約為1 636.65萬t。如果通過征稅，使得牛的養(yǎng)殖規(guī)模縮小20%，則每年大約可減少排放CO2當(dāng)量327.33萬t。

5.3 通過測土配方施肥技術(shù)減少N2O的排放

通過測土施肥技術(shù)提高氮肥利用率，避免過量施肥造成N2O排放。據(jù)推算，測土配方施肥將減少農(nóng)田N2O排放3%[15]。采用合理的糞便施用方式降低土地利用過程中溫室氣體的排放。研究發(fā)現(xiàn)，不同的施肥方式下N2O的排放通量存在著明顯的差異。采用撒施后翻耕和條施后覆土處理能有效地抑制NH3揮發(fā)和減少N2O等溫室氣體的排放，減輕對環(huán)境的不良影響[16]。2005～2009年，湖北省累計獲中央財政測土配方施肥補貼資金2.16億元，推廣實施面積1 750萬hm2，增產(chǎn)糧食34.75億kg、棉花283.1萬擔(dān)、油料4.45億kg，節(jié)本增效104.6億元[17]。2010年，國家測土配方施肥項目覆蓋了全省107個縣（市、區(qū)、場），肥料利用率提高1%，每公頃化肥用量比之前節(jié)約15 kg，全省免費為700萬農(nóng)戶提供測土配方施肥技術(shù)服務(wù)，推廣面積達到566.67萬hm2。雖然湖北省的測土配方施肥取得了一定成效，但是實際情況還遠不如人意，試驗和示范的較多，真正普及和自覺行動的少，施肥時依然“我行我素”、“跟著感覺走”，化肥還是沒少用[6]，因此，測土配方施肥還有待進一步落實。

5.4 以清潔能源替代化石燃料促使農(nóng)業(yè)機械低碳化

現(xiàn)代農(nóng)業(yè)的一個重要基石是農(nóng)業(yè)機械化。其對于利用資源、抵御自然災(zāi)害、推廣現(xiàn)代農(nóng)業(yè)技術(shù)、促進農(nóng)業(yè)集約經(jīng)營、增加單產(chǎn)與總產(chǎn)、提高勞動生產(chǎn)率、降低農(nóng)產(chǎn)品成本，以及對于減輕農(nóng)民勞動強度和縮小工農(nóng)差別等起到了積極作用。但是，農(nóng)業(yè)機械化因高度依賴化石燃料而造成對環(huán)境的污染與破壞也是無法彌補和挽回的。因此，后現(xiàn)代農(nóng)業(yè)發(fā)展不能沒有農(nóng)業(yè)機械，但又不能步現(xiàn)代農(nóng)業(yè)高污染、高耗能和高排放的后塵，而是應(yīng)該用清潔能源，如生物乙醇、生物柴油和生物石油等生物燃料替代化石燃料。同時，應(yīng)盡快淘汰已有的高耗能農(nóng)業(yè)機械和設(shè)備，以低耗能的機械和設(shè)備取代之。2009年，湖北機耕、機播、機械植保和機械收獲面積分別為409.9、63.7、387.2和255.5萬hm2，農(nóng)業(yè)機械總動力30 572.4萬kw·h。面對如此大的機械化規(guī)模，地方政府應(yīng)出臺鼓勵性政策，以發(fā)放補貼等形式將農(nóng)業(yè)機械與設(shè)備引導(dǎo)到低碳化的發(fā)展軌道上來。

參考文獻：

[1] 漆雁斌，陳衛(wèi)洪.低碳農(nóng)業(yè)發(fā)展影響因素的回歸分析[J].農(nóng)村經(jīng)濟，2010（2）：19-23.

[2] 漆雁斌，毛婷婷，殷凌霄.能源緊張情況下的低碳農(nóng)業(yè)發(fā)展問題分析[J]. 農(nóng)業(yè)技術(shù)經(jīng)濟，2010（3）：106-115.

[3] 李曉燕，王彬彬.低碳農(nóng)業(yè)：應(yīng)對氣候變化下的農(nóng)業(yè)發(fā)展之路[J].農(nóng)村經(jīng)濟，2010（3）：10-12.

[4] 中華人民共和國國家統(tǒng)計局.2010年中國統(tǒng)計年鑒[Z].北京：中國統(tǒng)計出版社，2011.

[5] 吳一平， 劉向華.發(fā)展低碳經(jīng)濟建設(shè)我國現(xiàn)代農(nóng)業(yè)[J]. 毛澤東鄧小平理論研究，2010（2）：58-65.

[6] 馬倫姣. 湖北省農(nóng)業(yè)低碳發(fā)展制約因素的實證探析[J]. 湖北農(nóng)業(yè)科學(xué)，2011，50（10）：2138-2141.

[7] 陳希勇.四川民族地區(qū)發(fā)展低碳農(nóng)業(yè)的制約因素及對策研究[J]. 貴州民族研究，2010，31（3）：127-131.

[8] 袁建文.計量經(jīng)濟學(xué)實驗教程[M].北京： 科學(xué)出版社，2008.

[9] 張成龍，柴沁虎，張阿玲，等.中國玉米生產(chǎn)的生產(chǎn)函數(shù)分析[J]. 清華大學(xué)學(xué)報（自然科學(xué)版），2009，49（12）：2024-2027.

[10] 趙安民. 論農(nóng)田的“白色污染”與防治途徑[J]. 山西農(nóng)業(yè)大學(xué)學(xué)報，2002（2）：179-183.

[11] OLESEN J E， SCHELDE K， WEISKE A. Modelling greenhouse gas emissions from European conventional and organic dairy farms[J]. Agriculture， Ecosystems and Environment，2006，112（2）：207-220.

[12] 高鐸山，孫 蕾.芻議政策性金融對低碳農(nóng)業(yè)經(jīng)濟的信貸支持[J].金融縱橫，2010（6）：49-51.

[13] 張憲英.我國低碳農(nóng)業(yè)解讀及其發(fā)展路徑初探[D].上海：復(fù)旦大學(xué)，2010.

[14] 孫亞男，劉繼軍，馬宗虎.規(guī)模化奶牛場溫室氣體排放量評估[J].農(nóng)業(yè)工程學(xué)報，2010，26（6）：296-301.

[15] 董紅敏，李玉娥，陶秀萍，等.中國農(nóng)業(yè)源溫室氣體排放與減排技術(shù)對策[J].農(nóng)業(yè)工程學(xué)報，2008，24（10）：269-273.

[16] 李 鑫，巨曉棠，張麗娟，等.不同施肥方式對土壤氨揮發(fā)和氧化亞氮排放的影響[J]. 應(yīng)用生態(tài)學(xué)報，2008，19（1）：99-104.

[17] 張愛虎，胡群中.測土配方五年增效百億元[N].湖北日報，2010-08-25（4）.