基于非期望產(chǎn)出模型的農(nóng)業(yè)綠色發(fā)展效率與影響因素研究
——以河北省11個地市為例

申思敏，王 亮

（1.首都經(jīng)濟(jì)貿(mào)易大學(xué)，北京 100070；2.西藏大學(xué)財經(jīng)學(xué)院，西藏 拉薩 850000）

1 引言

改革開放以來，特別是實行家庭聯(lián)產(chǎn)承包責(zé)任制以來，我國農(nóng)業(yè)生產(chǎn)規(guī)模、農(nóng)產(chǎn)品生產(chǎn)質(zhì)量都得到了極大程度的提升。我國糧食產(chǎn)量已連續(xù)多年保持在高位水平，穩(wěn)定在0.65萬億公斤以上，糧食安全得到了充分保障。然而，傳統(tǒng)粗放型農(nóng)業(yè)經(jīng)營模式依然主導(dǎo)著現(xiàn)有農(nóng)業(yè)生產(chǎn)運作過程。在多年的農(nóng)業(yè)高投入下，農(nóng)業(yè)生產(chǎn)與自然生態(tài)之間的矛盾日益突出，近年來東北黑土地土層變薄、農(nóng)產(chǎn)品重金屬污染、農(nóng)藥殘留等問題出現(xiàn)在公眾視野。與此同時，農(nóng)業(yè)生產(chǎn)效率低下、農(nóng)業(yè)從業(yè)人員文化素質(zhì)普遍較低、農(nóng)業(yè)勞動者年齡整體偏大等問題也十分突出，并由此導(dǎo)致規(guī)?；⒓s化農(nóng)業(yè)生產(chǎn)普及率不高。過度施肥、過度使用農(nóng)藥、過度依賴地膜等所導(dǎo)致的農(nóng)業(yè)生產(chǎn)過程的面源污染問題亟待解決。同時，在農(nóng)業(yè)生產(chǎn)過程中所引起的水資源浪費和環(huán)境污染嚴(yán)重等綠色發(fā)展問題均已到了非解決不可的關(guān)鍵時刻。

黨的十九大報告強(qiáng)調(diào)，要進(jìn)一步提高國家農(nóng)業(yè)安全保障水平，初步構(gòu)建現(xiàn)代農(nóng)業(yè)體系，推進(jìn)農(nóng)業(yè)綠色全面發(fā)展。同時，國家鄉(xiāng)村振興戰(zhàn)略明確提出，農(nóng)業(yè)農(nóng)村發(fā)展要本著生態(tài)宜居、可持續(xù)發(fā)展的思路展開。習(xí)近平總書記也多次強(qiáng)調(diào)，要保證糧食安全，確保中國人的碗里裝中國糧。糧食安全不僅涉及數(shù)量上的安全可控，更重要的是要確保質(zhì)量安全，只有這樣才能最終實現(xiàn)真正意義上的糧食安全。習(xí)近平總書記在聯(lián)合國大會上明確提出了中國“碳達(dá)峰”和“碳中和”的長遠(yuǎn)發(fā)展目標(biāo)，這對我國優(yōu)化發(fā)展結(jié)構(gòu)提出了更高的要求。從“碳達(dá)峰”到“碳中和”就是要控制溫室氣體排放，最終實現(xiàn)“零排放”。綠色發(fā)展就是要建立在生態(tài)、可持續(xù)基礎(chǔ)之上的以實現(xiàn)“零排放”為目標(biāo)的發(fā)展。農(nóng)業(yè)碳排放主要涉及甲烷、氮氧化物等溫室氣體，據(jù)相關(guān)研究測算，其排放量占總體碳排放量的比例超過1/3①。因此，基于碳足跡視角，將農(nóng)業(yè)生產(chǎn)的碳足跡作為農(nóng)業(yè)綠色發(fā)展中非期望產(chǎn)出的評價指標(biāo)，以此開展農(nóng)業(yè)綠色發(fā)展效率與影響因素研究具有重要的理論意義和現(xiàn)實意義，能夠為下一步農(nóng)業(yè)綠色化、低碳化以及可持續(xù)發(fā)展提供參考。

河北省作為我國重要的糧食主產(chǎn)區(qū)，且作為京津冀一體化的樞紐地區(qū)，其戰(zhàn)略定位的重要性毋庸置疑。然而，受經(jīng)濟(jì)結(jié)構(gòu)的影響，近年來經(jīng)濟(jì)發(fā)展所帶來的環(huán)境問題已經(jīng)成為河北省亟待解決的問題。農(nóng)業(yè)作為河北省經(jīng)濟(jì)發(fā)展過程中的重要支柱性產(chǎn)業(yè)，其綠色可持續(xù)發(fā)展問題不容忽視。對河北省農(nóng)業(yè)綠色發(fā)展效率進(jìn)行現(xiàn)狀判斷以及影響因素分析具有重要的理論意義和現(xiàn)實意義，既有助于河北省乃至京津冀地區(qū)鄉(xiāng)村振興戰(zhàn)略的進(jìn)一步落實，又有助于因地制宜地解決好經(jīng)濟(jì)高質(zhì)量發(fā)展問題，提升河北省經(jīng)濟(jì)發(fā)展質(zhì)量。

2 文獻(xiàn)綜述

“綠色發(fā)展”概念最早由聯(lián)合國開發(fā)計劃署（UNDP）于2002年提出，2011年OECD在界定綠色增長的內(nèi)涵和外延的基礎(chǔ)上提出綠色發(fā)展指南，并呼吁全球各國樹立綠色經(jīng)濟(jì)發(fā)展理念。Barbiroli（2011）、Ling等（2011）認(rèn)為，農(nóng)業(yè)綠色發(fā)展方式應(yīng)同時滿足糧食安全、能源有效使用和技術(shù)創(chuàng)新3個方面［1-2］；而Koohafkan等（2012）則認(rèn)為，建立一套靈活且具有針對性的可持續(xù)發(fā)展體系是實現(xiàn)綠色農(nóng)業(yè)的重要保障［3］。

國內(nèi)對于農(nóng)業(yè)綠色發(fā)展的研究，按照不同的研究視角可以大致分為農(nóng)業(yè)綠色發(fā)展理論內(nèi)涵、農(nóng)業(yè)綠色發(fā)展評價指標(biāo)體系構(gòu)建等方面。其中，于法穩(wěn)（2016）、羅必良（2017）認(rèn)為，綠色農(nóng)業(yè)就是指推動經(jīng)濟(jì)、生態(tài)和社會可持續(xù)發(fā)展的農(nóng)業(yè)發(fā)展模式［4-5］。自鄉(xiāng)村振興戰(zhàn)略規(guī)劃實施以來，越來越多的學(xué)者開始思考農(nóng)業(yè)綠色發(fā)展的具體規(guī)劃與實施戰(zhàn)略問題。張宇（2019）通過分析我國農(nóng)村生態(tài)環(huán)境面臨的三大問題（資源、環(huán)境、資源管理機(jī)制），提出在農(nóng)業(yè)綠色發(fā)展過程中必須推進(jìn)綠色農(nóng)業(yè)產(chǎn)業(yè)化、農(nóng)業(yè)科技創(chuàng)新等關(guān)鍵領(lǐng)域發(fā)展［6］。高尚賓（2019）指出，生態(tài)農(nóng)業(yè)的建設(shè)應(yīng)當(dāng)首先解決農(nóng)業(yè)生產(chǎn)所面臨的資源環(huán)境壓力［7］。對于農(nóng)業(yè)綠色發(fā)展評價體系建設(shè)，國內(nèi)學(xué)者也從各個角度進(jìn)行了相關(guān)研究，其中代表性學(xué)者包括：任運河（2006）、吳丹（2017）、魏琦（2018）、孫煒琳（2019）等，這些學(xué)者大多從資源、環(huán)境、質(zhì)量等多個維度對農(nóng)業(yè)綠色發(fā)展水平進(jìn)行了評價［8-11］；而具體對農(nóng)業(yè)綠色發(fā)展效率進(jìn)行研究的學(xué)者很少。郭付友等（2020）采用DEA-Malmquist指數(shù)法分析了魯南經(jīng)濟(jì)帶35個縣綠色發(fā)展效率，并使用城鄉(xiāng)綠色發(fā)展效率脫鉤狀態(tài)模型、地理探測器等多種計量模型，對2005—2017年魯南經(jīng)濟(jì)帶城鄉(xiāng)綠色發(fā)展效率時空分異特征與驅(qū)動因素進(jìn)行綜合分析［12］。薛蕾等（2020）則通過使用SBM-Un?desirable模型，核算了2004—2015年我國30個?。▍^(qū)市）的農(nóng)業(yè)綠色發(fā)展效率，并基于SAR空間計量模型，驗證了農(nóng)業(yè)產(chǎn)業(yè)集聚對農(nóng)業(yè)綠色發(fā)展效率的空間作用機(jī)制［13］。

然而，以上對于農(nóng)業(yè)綠色發(fā)展的研究，很少基于農(nóng)業(yè)綠色發(fā)展效率角度的省域?qū)嵶C分析，同時對農(nóng)業(yè)綠色發(fā)展效率進(jìn)行分析時忽略了效率值作為一個受限變量的基本特征。因此，本文以河北省為研究對象，創(chuàng)造性地將碳足跡作為農(nóng)業(yè)生產(chǎn)過程中的非期望產(chǎn)出，構(gòu)建了非期望產(chǎn)出的超效率SBM模型，核算了2008—2017年河北省11個地市的農(nóng)業(yè)綠色發(fā)展效率②，并利用面板Tobit模型分析了相應(yīng)的影響因素，以利于河北省農(nóng)業(yè)綠色政策的制定和實施，同時為我國其他地區(qū)農(nóng)業(yè)綠色發(fā)展效率問題研究提供參考范例。

3 研究方法

本文基于碳足跡視角，對農(nóng)業(yè)生產(chǎn)過程中的碳排放和農(nóng)作物生產(chǎn)過程中的碳吸收進(jìn)行核算，將兩者結(jié)合作為碳足跡，并將其視為農(nóng)業(yè)綠色發(fā)展中的非期望產(chǎn)出；同時構(gòu)建了超效率SBM模型，對2008—2017年河北省11個地市的農(nóng)業(yè)綠色發(fā)展效率進(jìn)行了測度。

3.1 碳足跡核算

3.1.1 碳排放量核算。在進(jìn)行碳排放量指標(biāo)核算的過程中，本研究采用IPCC發(fā)布的碳排放系數(shù)法并結(jié)合河北省各個地市的相關(guān)農(nóng)業(yè)生產(chǎn)數(shù)據(jù)，對河北省各個地市碳排放總量進(jìn)行核算，具體的核算公式為：

其中：，Ct表示第t年的碳排放總量，即每年所有碳排放來源（以下簡稱碳源）排放碳量Cit的總和；TYit表示第t年第i種碳源的量；Qit表示碳源的碳排放系數(shù)。本文主要核算了農(nóng)藥、農(nóng)用塑料薄膜、化肥、農(nóng)業(yè)灌溉所耗費的電能、農(nóng)用柴油等五大碳源的排放，各碳源的排放系數(shù)如表1所示。

表1 各碳源碳排放系數(shù)

3.1.2 碳吸收量核算。農(nóng)作物用于光合作用所吸收的碳總量主要由農(nóng)作物進(jìn)行光合作用合成干物質(zhì)所需要的碳量AH（又稱經(jīng)濟(jì)系數(shù)）與農(nóng)作物總生物量BS的乘積決定，具體的核算公式為：

其中：Ew為農(nóng)作物的經(jīng)濟(jì)產(chǎn)量；wi為各種作物的含水量；Hi為農(nóng)作物的碳吸收率。本文根據(jù)河北省農(nóng)業(yè)種植作物的特點，選取種植較為廣泛的以下幾種農(nóng)作物作為碳吸收源進(jìn)行核算（見表2）。

表2 主要農(nóng)作物含水量、經(jīng)濟(jì)系數(shù)、碳吸收率數(shù)值表

3.1.3 農(nóng)業(yè)碳足跡核算。本文參照段華平等［16］的研究方法，具體的核算公式如下：

其中：CEFt表示農(nóng)業(yè)碳足跡；NEPt表示農(nóng)作物固碳潛力；CXt表示農(nóng)作物進(jìn)行光合作用時對碳的吸收量；s表示農(nóng)作物生產(chǎn)使用的耕地面積。

3.2 農(nóng)業(yè)綠色發(fā)展效率核算

傳統(tǒng)的效率測定方法忽視了生產(chǎn)過程中產(chǎn)生非期望結(jié)果的問題，因而對于發(fā)展效率的測度存在一定缺陷。為解決此缺陷造成的發(fā)展效率核算精確度不高問題，Tone提出SBM-DEA模型，將非期望產(chǎn)出作為模型構(gòu)成部分，這提高了效率模型測度的準(zhǔn)確性。但是，此模型中只要存在效率大于1的決策單元，最后的效率輸出值均為1。為克服此問題帶來的效率大于1的決策單元難以彼此比較的問題，Tone又進(jìn)一步改進(jìn)了此模型，提出了超效率SBM模型［17］。本文利用超效率SBM模型，對2008—2017年河北省11個地市的農(nóng)業(yè)綠色發(fā)展效率進(jìn)行核算，具體的核算公式為（4）—（7）：

其中：EFF為農(nóng)業(yè)綠色發(fā)展的效率值；N、M、I分別為綠色農(nóng)業(yè)的投入、產(chǎn)出和非期望產(chǎn)出指標(biāo)；x、y、b分別為3個指標(biāo)的向量；分別為3個指標(biāo)的松弛變量；為k'生產(chǎn)單元在第t'時期的投入、產(chǎn)出、非期望產(chǎn)出值；為決策單元的權(quán)重系數(shù)。

3.3 農(nóng)業(yè)綠色發(fā)展效率影響因素核算

通過超效率SBM模型計算的農(nóng)業(yè)綠色發(fā)展效率取值范圍在［0，1.2］之間，大多數(shù)的數(shù)據(jù)分布于1附近，因而農(nóng)業(yè)綠色發(fā)展效率屬于受限制因變量模型的研究范圍。若使用最小二乘法（OLS）將導(dǎo)致估計結(jié)果有偏和不一致，因而本文基于準(zhǔn)確性和穩(wěn)定性原則，選擇面板Tobit模型進(jìn)行分析。具體的面板Tobit模型為式（8）和（9）：

其中：X表示各個自變量的值；β表示各個自變量的系數(shù)；i=1,2…N，表示研究的個體維度數(shù)量；t=1,2…T表示被觀察個體跨越的時間維度；μit表示個體時期恒量；vit表示殘差項的擾動。

4 研究數(shù)據(jù)來源與變量選取

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

本文選取的碳排放量、農(nóng)業(yè)綠色發(fā)展效率以及農(nóng)業(yè)綠色發(fā)展效率影響因素相關(guān)指標(biāo)數(shù)據(jù)主要來源于《河北農(nóng)村年鑒》《河北經(jīng)濟(jì)年鑒》以及《河北統(tǒng)計公報》?；跀?shù)據(jù)的可獲得性、可比性和完整性等原則，本研究最終選取2008—2017年河北省的相關(guān)數(shù)據(jù)作為研究基礎(chǔ)。

4.2 變量選取

4.2.1 碳足跡核算變量。①碳排放量。本文對于碳排放量的核算，主要以農(nóng)業(yè)生產(chǎn)過程中常用的農(nóng)藥、農(nóng)用塑料薄膜、化肥、農(nóng)業(yè)生產(chǎn)灌溉過程中所耗費的電能以及農(nóng)用機(jī)械在輔助生產(chǎn)過程中使用的柴油等作為碳排放量的主要碳源。

②碳吸收量。農(nóng)作物生產(chǎn)過程中的光合作用會吸收二氧化碳，并將其轉(zhuǎn)化為含碳化合物。在此過程中，部分由于農(nóng)業(yè)生產(chǎn)所產(chǎn)生的碳排放會被農(nóng)作物吸收。為了衡量這種機(jī)制對含碳化合物的吸收作用，本文對小麥、水稻、玉米、谷子等常見的農(nóng)作物碳吸收量均做了核算，并將其加總，求出總的碳吸收量。

4.2.2 農(nóng)業(yè)綠色發(fā)展效率核算變量。對于農(nóng)業(yè)綠色發(fā)展效率值核算，本文選取的農(nóng)業(yè)相關(guān)投入指標(biāo)主要有農(nóng)業(yè)機(jī)械化水平、有效灌溉率以及單位化肥使用量，產(chǎn)出指標(biāo)主要有糧食播種面積單產(chǎn)、油料播種面積單產(chǎn)以及蔬菜播種面積單產(chǎn)；對于農(nóng)業(yè)生產(chǎn)過程的非期望指標(biāo)則利用農(nóng)業(yè)碳足跡進(jìn)行表征。其中，投入指標(biāo)的詳細(xì)核算方法見表3。

表3 農(nóng)業(yè)相關(guān)的主要投入指標(biāo)核算方法

4.2.3 農(nóng)業(yè)綠色發(fā)展效率影響因素。對于農(nóng)業(yè)綠色發(fā)展效率的影響因素，本文主要從社會層面、經(jīng)濟(jì)層面、生態(tài)層面和產(chǎn)業(yè)層面進(jìn)行分析。社會層面主要從勞動力素質(zhì)出發(fā)進(jìn)行考慮，勞動力素質(zhì)的提升有助于農(nóng)業(yè)科技發(fā)展，進(jìn)而提升農(nóng)業(yè)生產(chǎn)效率和農(nóng)業(yè)可持續(xù)能力，本文選取大專及以上學(xué)歷人口占總?cè)丝诘谋戎兀≒er）作為該層面的測度指標(biāo)。經(jīng)濟(jì)層面主要考慮各市的經(jīng)濟(jì)發(fā)展情況，經(jīng)濟(jì)水平的提升能夠帶來更好的農(nóng)業(yè)基礎(chǔ)生產(chǎn)條件和農(nóng)業(yè)發(fā)展環(huán)境，因而選取地區(qū)生產(chǎn)總值（GDP）作為該層面的衡量指標(biāo)。生態(tài)層面主要從各市農(nóng)業(yè)生產(chǎn)中面臨的生態(tài)環(huán)境問題出發(fā)進(jìn)行考慮，農(nóng)業(yè)是對自然環(huán)境依賴性較強(qiáng)的產(chǎn)業(yè)，環(huán)境帶來的挑戰(zhàn)無疑會對其發(fā)展效率產(chǎn)生影響，選取受災(zāi)率（Dia）作為該層面的衡量指標(biāo)。鄉(xiāng)村振興背景下農(nóng)業(yè)產(chǎn)業(yè)化發(fā)展已經(jīng)成為現(xiàn)代化農(nóng)業(yè)產(chǎn)業(yè)鏈建設(shè)的主要途徑，因而對于農(nóng)業(yè)綠色發(fā)展的研究也不能忽略產(chǎn)業(yè)層面的影響，選取產(chǎn)業(yè)化經(jīng)營率（Ind）以及龍頭經(jīng)營組織總數(shù)（Log）作為衡量指標(biāo)，測度農(nóng)業(yè)產(chǎn)業(yè)化經(jīng)營對農(nóng)業(yè)綠色發(fā)展的影響。

5 實證結(jié)果及分析

5.1 河北省農(nóng)業(yè)綠色發(fā)展水平結(jié)果分析

5.1.1 2008—2017年河北省11個地市農(nóng)業(yè)碳足跡變遷。從河北省2008—2017年碳足跡變遷圖來看（如圖1所示），碳足跡整體呈現(xiàn)微弱下降趨勢。這說明隨著經(jīng)濟(jì)水平的提升以及各級政府對于生態(tài)環(huán)境的重視，人們的環(huán)保意識有了一定程度的提高，農(nóng)業(yè)綠色生產(chǎn)條件也有了一定改善。但是，總體而言，這種下降的趨勢并不明顯。這說明在農(nóng)業(yè)投入產(chǎn)出過程中仍存在大量的碳排放，這對農(nóng)業(yè)自然生產(chǎn)條件改善和農(nóng)業(yè)綠色可持續(xù)發(fā)展具有一定的負(fù)向影響。同時，可以看出各地市碳足跡分布并不均衡，且與其經(jīng)濟(jì)發(fā)展?fàn)顩r緊密相關(guān)。經(jīng)濟(jì)發(fā)展水平較高的石家莊市、唐山市碳足跡反而比較高，這說明經(jīng)濟(jì)發(fā)展水平的提升并不能帶來人們對于發(fā)展可持續(xù)性認(rèn)識的提升。這在某種程度上也表明，這些城市的經(jīng)濟(jì)發(fā)展大多依賴于高能耗、高污染的產(chǎn)業(yè)。此外，這10年間，農(nóng)業(yè)綠色生產(chǎn)條件改善程度最大的城市是衡水市，其次是保定市，其他城市未見劇烈的碳足跡波動。

圖1 2008—2017年河北省各地市碳足跡變遷圖

5.1.2 農(nóng)業(yè)綠色發(fā)展效率結(jié)果及空間分布。本研究基于非期望產(chǎn)出的超效率SBM模型，利用Matlab2018b軟件核算了河北省11個地市在2008—2017年的農(nóng)業(yè)綠色發(fā)展效率（如表4所示），并從多個角度對河北省各地市農(nóng)業(yè)綠色發(fā)展效率進(jìn)行比較。

表4 2008—2017年河北省各地市農(nóng)業(yè)綠色發(fā)展效率值

從總體層面看，2008—2017年各個地市的農(nóng)業(yè)綠色發(fā)展效率都有了不同程度的提高（廊坊市除外）。其中，效率提高最明顯的是唐山市，與2008年相比，2017年其農(nóng)業(yè)發(fā)展效率值提升了0.368 8；其次是保定市，農(nóng)業(yè)發(fā)展效率提升了0.362 8。這說明唐山市和保定市具有較好的農(nóng)業(yè)綠色發(fā)展基礎(chǔ)，在政策調(diào)節(jié)下初步取得了一定的農(nóng)業(yè)綠色發(fā)展成效。

從平均層面看，河北省農(nóng)業(yè)綠色發(fā)展效率整體偏低，尤其是石家莊市、邢臺市和滄州市，平均發(fā)展效率在0.550 0～0.560 0之間。石家莊市作為河北省省會城市，其農(nóng)業(yè)綠色發(fā)展方向和發(fā)展結(jié)構(gòu)應(yīng)當(dāng)對其他地市起示范作用，但現(xiàn)實并非如此。這種經(jīng)濟(jì)發(fā)展和綠色發(fā)展相背離的現(xiàn)實應(yīng)當(dāng)引起當(dāng)?shù)卣闹匾暋?/p>

從2017年河北省各個地市的農(nóng)業(yè)綠色發(fā)展效率來看，農(nóng)業(yè)綠色發(fā)展效率最高的城市是秦皇島市，其農(nóng)業(yè)綠色發(fā)展效率值為1.181 7，與其他地區(qū)相比效率較高，表明其農(nóng)業(yè)綠色發(fā)展水平較高。除此之外，承德市和張家口市的農(nóng)業(yè)綠色發(fā)展效率也較高，分別為1.137 9和1.075 7。其他城市的發(fā)展效率均小于等于1，作為省會城市的石家莊市最低，只有0.634 9。這說明河北省整體農(nóng)業(yè)綠色發(fā)展效率仍有待提升，尤其像石家莊市、保定市等經(jīng)濟(jì)發(fā)展水平較高且農(nóng)業(yè)產(chǎn)業(yè)作為經(jīng)濟(jì)重要組成部分的城市，在發(fā)展二三產(chǎn)業(yè)的同時，也應(yīng)注意農(nóng)業(yè)發(fā)展的綠色化和可持續(xù)性問題。

將各個地市10年間的農(nóng)業(yè)綠色發(fā)展效率指數(shù)平均后導(dǎo)入地理信息系統(tǒng)（GIS）中（如圖2所示），結(jié)合空間分類方法，將其簡單分為4類?？梢钥闯觯颖笔∞r(nóng)業(yè)綠色發(fā)展呈現(xiàn)出很強(qiáng)的空間“聚團(tuán)”性，農(nóng)業(yè)綠色發(fā)展程度較高的地市多數(shù)位于河北省的北部，而河北省南部的各個地市在農(nóng)業(yè)綠色發(fā)展方面的表現(xiàn)要明顯差于河北省北部的各個地市。這說明農(nóng)業(yè)綠色發(fā)展存在一定的空間傳播效應(yīng)，通過借鑒鄰接城市實施的相關(guān)農(nóng)業(yè)政策，可以在一定程度上加快本地區(qū)的農(nóng)業(yè)綠色發(fā)展進(jìn)程，這也給政策制定帶來一些啟示。

圖2 河北省各地市農(nóng)業(yè)綠色發(fā)展平均效率分布圖

5.2 河北省農(nóng)業(yè)綠色發(fā)展影響因素研究

針對受限因變量數(shù)據(jù)的面板Tobit模型，可分為面板Tobit固定效應(yīng)模型和面板Tobit隨機(jī)效應(yīng)模型，利用相關(guān)數(shù)據(jù)分別進(jìn)行擬合，兩種模型的擬合結(jié)果如表5所示。根據(jù)擬合結(jié)果可以看出，面板Tobit固定效應(yīng)模型和隨機(jī)效應(yīng)模型各個影響因素系數(shù)除受災(zāi)率之外均比較顯著。為選擇最優(yōu)模型，減少模型擬合帶來的誤差，通過Haus?man檢驗進(jìn)行模型選擇，具體檢驗結(jié)果見表6。由Hausman檢驗的結(jié)果可知，在P=0.000 0的情況下，面板Tobit固定效應(yīng)模型的擬合效果要好于隨機(jī)效應(yīng)模型，因而本文選擇面板Tobit固定效應(yīng)模型的結(jié)果進(jìn)行分析。

表5 兩種面板Tobit模型的回歸結(jié)果

表6 Hausman檢驗結(jié)果

根據(jù)面板Tobit固定效應(yīng)模型的回歸結(jié)果可以看出，所選變量除了受災(zāi)率（Dia）對農(nóng)業(yè)綠色發(fā)展效率無顯著影響之外，其他變量均對農(nóng)業(yè)綠色發(fā)展效率具有一定程度的影響。其中，勞動力的素質(zhì)水平（Per）在1%的顯著水平下對農(nóng)業(yè)綠色發(fā)展效率的提升具有較強(qiáng)的促進(jìn)作用；農(nóng)業(yè)產(chǎn)業(yè)化經(jīng)營率（Ind）在1%的水平下對農(nóng)業(yè)綠色發(fā)展效率有顯著的正向影響。反觀地區(qū)生產(chǎn)總值（GDP）（在1%的水平下）和龍頭經(jīng)營組織總數(shù)（Log）（在1%的水平下）對于農(nóng)業(yè)綠色發(fā)展效率的提升則有負(fù)向影響?？赡艿脑蚴?，河北省經(jīng)濟(jì)發(fā)展仍依賴大量的化石燃料，以能源消耗帶動經(jīng)濟(jì)增長。這種粗放型發(fā)展模式雖然可以在短時間內(nèi)有效帶動經(jīng)濟(jì)高速增長，但對生態(tài)環(huán)境造成的破壞會對經(jīng)濟(jì)長期增長產(chǎn)生負(fù)向作用。受災(zāi)率（Dia）對于農(nóng)業(yè)綠色發(fā)展效率的影響并不顯著，可能的原因在于受災(zāi)率主要受生態(tài)系統(tǒng)自我調(diào)節(jié)的影響，不會受到人類行為過多的干擾。

6 結(jié)論與政策啟示

6.1 結(jié)論

本文根據(jù)河北省2008—2017年的省域面板數(shù)據(jù)，主要討論了2008—2017年河北省各個地市的農(nóng)業(yè)綠色發(fā)展效率變遷以及影響農(nóng)業(yè)綠色發(fā)展效率的主要因素，研究結(jié)論如下：

6.1.1 河北省各地市的碳足跡呈現(xiàn)較微弱的下降趨勢。通過對2008—2017年河北省11個地市碳足跡的核算以及碳足跡均值的核算可以看出，河北省各地市農(nóng)業(yè)碳足跡整體呈現(xiàn)較微弱的下降趨勢，只有少部分地區(qū)（如衡水市）呈現(xiàn)出較明顯的波動趨勢，這說明河北省在農(nóng)業(yè)關(guān)聯(lián)的碳排放方面仍有較大的提升空間。

6.1.2 河北省農(nóng)業(yè)綠色發(fā)展效率有明顯的地域差距。河北省整體農(nóng)業(yè)綠色發(fā)展效率在0.5～1.2之間，整體發(fā)展效率處于良好水平。部分地市農(nóng)業(yè)綠色發(fā)展效率較高，如秦皇島市、張家口市和承德市；但部分地市農(nóng)業(yè)綠色發(fā)展效率較低，尤其是作為省會城市的石家莊市，農(nóng)業(yè)綠色發(fā)展效率明顯低于其他城市，依賴二、三產(chǎn)業(yè)的經(jīng)濟(jì)發(fā)展結(jié)構(gòu)可能是其農(nóng)業(yè)綠色發(fā)展效率低下的原因。

6.1.3 河北省農(nóng)業(yè)綠色發(fā)展效率呈現(xiàn)出地域“聚團(tuán)”現(xiàn)象。由農(nóng)業(yè)綠色發(fā)展效率的空間分布圖可以看出，農(nóng)業(yè)綠色發(fā)展效率呈現(xiàn)出“聚團(tuán)”現(xiàn)象。效率較高的城市如張家口市、承德市，大多分布于河北省北部；效率較低的城市如邢臺市、邯鄲市等，大多分布于河北省南部。這說明河北省農(nóng)業(yè)綠色發(fā)展存在一定的地域特點，針對不同地區(qū)應(yīng)當(dāng)采取不同的農(nóng)業(yè)綠色化發(fā)展模式。

6.1.4 不同經(jīng)濟(jì)指標(biāo)的提升使得農(nóng)業(yè)綠色發(fā)展效率差異化變動。根據(jù)面板Tobit固定效應(yīng)模型擬合結(jié)果可以看出，勞動力素質(zhì)、產(chǎn)業(yè)化發(fā)展對農(nóng)業(yè)綠色發(fā)展效率的提升有較強(qiáng)的正向作用。因此，提升農(nóng)業(yè)從業(yè)者素質(zhì)、加強(qiáng)農(nóng)業(yè)產(chǎn)業(yè)化經(jīng)營力度能有效提高區(qū)域農(nóng)業(yè)綠色發(fā)展效率。而GDP和龍頭經(jīng)營組織總數(shù)這類經(jīng)濟(jì)指標(biāo)的提升則會造成農(nóng)業(yè)綠色發(fā)展效率降低，原因在于其經(jīng)濟(jì)發(fā)展結(jié)構(gòu)中可能有高能源消耗成分的存在。

6.2 政策啟示

結(jié)合以上實證分析結(jié)果，以及農(nóng)業(yè)綠色發(fā)展效率較高地市的發(fā)展經(jīng)驗，本文提出以下政策建議：

6.2.1 因地施策，結(jié)合區(qū)域情況發(fā)展農(nóng)業(yè)綠色化。對于農(nóng)業(yè)綠色發(fā)展效率較高的地區(qū)，應(yīng)當(dāng)保持現(xiàn)有農(nóng)業(yè)綠色發(fā)展勢頭，繼續(xù)深化農(nóng)業(yè)產(chǎn)業(yè)化改革，進(jìn)一步提升農(nóng)業(yè)綠色化水平；對于農(nóng)業(yè)綠色發(fā)展效率中等的地區(qū)，應(yīng)當(dāng)多吸取農(nóng)業(yè)綠色發(fā)展效率較高地區(qū)的發(fā)展經(jīng)驗，穩(wěn)中求進(jìn)，提升其農(nóng)業(yè)綠色化水平；對于農(nóng)業(yè)綠色發(fā)展效率低下且將農(nóng)業(yè)作為其主導(dǎo)產(chǎn)業(yè)的地區(qū)，應(yīng)當(dāng)大力扶持綠色農(nóng)業(yè)相關(guān)產(chǎn)業(yè)鏈建設(shè)，治理和改進(jìn)高耗能、高污染的生產(chǎn)環(huán)節(jié)，同時要加大農(nóng)業(yè)從業(yè)者對于農(nóng)業(yè)綠色化的認(rèn)識，促進(jìn)該地區(qū)農(nóng)業(yè)綠色發(fā)展效率提升。

6.2.2 優(yōu)化現(xiàn)有農(nóng)業(yè)綠色發(fā)展方式，完善綠色農(nóng)業(yè)供應(yīng)鏈體系。政府應(yīng)當(dāng)加強(qiáng)政策引導(dǎo)，如制定生物化肥補貼、有機(jī)農(nóng)藥補貼等政策，有意識地引導(dǎo)農(nóng)業(yè)生產(chǎn)行為，使整個農(nóng)業(yè)生產(chǎn)過程向綠色化、高效化轉(zhuǎn)變；同時，還應(yīng)注重農(nóng)業(yè)生態(tài)技術(shù)研發(fā)與產(chǎn)品推廣，引導(dǎo)農(nóng)村向綠色化、生態(tài)化方向發(fā)展。政府還應(yīng)借助其影響力，幫助對接農(nóng)產(chǎn)品生產(chǎn)方與供應(yīng)方，進(jìn)一步加快農(nóng)業(yè)產(chǎn)業(yè)化、市場化進(jìn)程，構(gòu)建具有現(xiàn)代化特征的綠色農(nóng)業(yè)供應(yīng)鏈體系。

6.2.3 規(guī)范和引導(dǎo)農(nóng)業(yè)從業(yè)者綠色生產(chǎn)行為，提高勞動者素質(zhì)。在農(nóng)業(yè)生產(chǎn)過程中，首先，要提高農(nóng)業(yè)從業(yè)者對農(nóng)業(yè)綠色發(fā)展的認(rèn)識，提高勞動力素質(zhì)水平，可以通過舉辦知識講座、農(nóng)業(yè)綠色生產(chǎn)行為宣傳推廣活動等形式，讓農(nóng)業(yè)從業(yè)者認(rèn)識到農(nóng)業(yè)綠色發(fā)展的重要性；其次，可通過開設(shè)綠色農(nóng)業(yè)技術(shù)課堂，引導(dǎo)農(nóng)業(yè)從業(yè)者改變不良生產(chǎn)行為；最后，還可以通過展示綠色農(nóng)業(yè)發(fā)展成果，激勵更多的微觀主體參與到農(nóng)業(yè)綠色發(fā)展行動中來。

注釋：

①此數(shù)據(jù)來源于聯(lián)合國政府間氣候變化專門委員會（IPCC）和聯(lián)合國糧食及農(nóng)業(yè)組織（FOA）報告。

②其中辛集市核算入石家莊市，定州市核算入保定市。