總量與結構視角耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響
——以湖北省為例

胡偉艷，朱慶瑩，張安錄，George Ayea Karnwie，李景旺

（華中農(nóng)業(yè)大學公共管理學院，湖北 武漢 430070）

1 引言

生產(chǎn)要素對農(nóng)業(yè)經(jīng)濟增長的影響歷來是經(jīng)濟理論研究的核心問題[1]，國內外大量文獻對中國的研究表明，耕地作為物質性資產(chǎn)對農(nóng)業(yè)經(jīng)濟增長的作用逐漸變小甚至為負[2-4]。中國工業(yè)化和城鎮(zhèn)化進程中城鎮(zhèn)空間變化引發(fā)了大量征地問題，優(yōu)質耕地大量減少[5]，農(nóng)業(yè)經(jīng)濟增長陷入雙重窘境，迫切需要進行耕地利用方式轉型。中國政府于 2007 年以來強力推進“耕地紅線”、“生態(tài)紅線”劃定和“生態(tài)文明建設”，2017年中央一號文件提出將推行農(nóng)業(yè)生產(chǎn)環(huán)境治理和綠色生產(chǎn)方式、增強農(nóng)業(yè)可持續(xù)發(fā)展能力作為供給側改革的重要內容，這些政策要求農(nóng)業(yè)生產(chǎn)向多功能利用方向轉變。作為農(nóng)業(yè)經(jīng)營最主要的元素和載體，耕地的多功能利用可達成農(nóng)業(yè)的多功能[6]，可能是培育農(nóng)業(yè)經(jīng)濟增長的新動能新途徑。

黨的十八大和十八屆三中、四中全會指出：當前，中國經(jīng)濟增長進入新常態(tài)，要素的總量投入對經(jīng)濟增長固然重要，但合理的要素利用結構與方式對經(jīng)濟增長更為關鍵，農(nóng)業(yè)領域更是如此。目前，許多學者在耕地多功能的內涵[7-8]、指標體系的建構[9]、耕地多功能評價[10]以及時空演變特征[11-12]等領域側重于耕地多功能總量上的研究，已有文獻開始探討耕地多功能權衡與協(xié)同或者沖突的結構性研究[10]，但從總量和結構兩個方面探討耕地多功能利用對農(nóng)業(yè)經(jīng)濟增長影響的研究還鮮有報道。從總量與結構視角研究耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響對深化農(nóng)業(yè)農(nóng)村可持續(xù)發(fā)展理論[13]、提升農(nóng)業(yè)經(jīng)濟增長質量和效率[14]、破解中國“三農(nóng)”等實踐問題[13]具有重要意義。湖北省作為中部地區(qū)的核心省份和全國重要的糧棉油商品生產(chǎn)基地，于2011年印發(fā)《湖北省主體功能區(qū)劃》，設定重點開發(fā)區(qū)、農(nóng)產(chǎn)品主產(chǎn)區(qū)和生態(tài)功能區(qū)。對其他地區(qū)的研究表明不同主體功能區(qū)耕地多功能的顯化及主導功能存在較大差異[12]，由此引起的耕地多功能總量與結構也存在明顯差異。本文以湖北省為例，在對耕地多功能利用對農(nóng)業(yè)經(jīng)濟增長進行理論分析的基礎上，擴展C-D道格拉斯生產(chǎn)函數(shù)模型，利用70個縣2000—2014年統(tǒng)計數(shù)據(jù)，從總量和結構視角探討不同主體功能區(qū)耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響，為差異化實施耕地多功能利用與管理，促進農(nóng)業(yè)經(jīng)濟持續(xù)發(fā)展和增長提供參考。

2 理論基礎與研究假說

耕地多功能是指耕地系統(tǒng)提供滿足人類生存與發(fā)展所需產(chǎn)品與服務的各種能力[15]，包括提供糧食、蔬菜、油料等生產(chǎn)功能；調節(jié)大氣、水文、氣候等生態(tài)功能；提供農(nóng)民生活保障和就業(yè)、維護社會糧食安全等社會功能和提供耕作風景、開敞空間景觀等文化功能[16]，包括總量上的耕地多功能和結構上的耕地多功能[8，10]。其中，總量耕地多功能是指一定區(qū)域耕地各類產(chǎn)品或服務顯化數(shù)量的加權和[9-10]，對農(nóng)業(yè)經(jīng)濟增長的影響機制為：生產(chǎn)功能提供農(nóng)產(chǎn)品，對農(nóng)業(yè)經(jīng)濟效益產(chǎn)生直接影響，生產(chǎn)功能越強，農(nóng)業(yè)經(jīng)濟增長越明顯；社會功能提供農(nóng)民就業(yè)機會、降低外出務工的風險，提高農(nóng)民福利水平，間接影響農(nóng)業(yè)生產(chǎn)的效益與效率[17]；生態(tài)功能增強能夠改善耕地本底質量，鞏固農(nóng)業(yè)再生產(chǎn)基礎，減少自然災害的發(fā)生和要素投入，降低農(nóng)業(yè)產(chǎn)出損失和生產(chǎn)成本[18]；景觀文化功能吸引游客觀光休閑，增加游客對農(nóng)產(chǎn)品等的價值認同[19]。劉衛(wèi)東指出綜合功能的耕地利用具有明顯的優(yōu)勢[20]，耕地多功能之間融合協(xié)同，總量視角耕地多功能越強，能促進農(nóng)業(yè)經(jīng)濟增長。據(jù)此提出假說1。

假說1：總量視角耕地多功能對農(nóng)業(yè)經(jīng)濟增長可能起正向作用。

結構視角耕地多功能是指一定區(qū)域耕地各類功能之間的比例關系或組合狀況，體現(xiàn)耕地功能類型的多樣性與數(shù)量的均勻度[10，21]。理論上，耕地功能之間的競爭或權衡，會引起耕地功能總量的減少，進而可能對農(nóng)業(yè)經(jīng)濟增長產(chǎn)生負向影響，甚至降低農(nóng)民福祉；相反，如果耕地功能之間存在協(xié)同，則會增加耕地功能總量，可能會對農(nóng)業(yè)經(jīng)濟增長產(chǎn)生正向影響，從而增加農(nóng)民福利[22]。實證研究表明耕地生產(chǎn)功能與生態(tài)環(huán)境功能之間存在競爭關系[23]；生產(chǎn)功能與社會功能、社會功能與生態(tài)環(huán)境功能之間呈負向關聯(lián)[10]。結構視角耕地多功能越強，各類耕地功能數(shù)量越均勻，功能之間的競爭或權衡更劇烈，耕地多功能總量減少。結合假說1提出假說2。

假說2：結構視角耕地多功能對農(nóng)業(yè)經(jīng)濟增長可能起負向作用。

湖北省不同主體功能區(qū)耕地資源稟賦、區(qū)位條件、經(jīng)濟發(fā)展水平、公眾對耕地多功能的需求等自然條件、社會經(jīng)濟條件存在很大差異[12]。其中，重點發(fā)展區(qū)社會經(jīng)濟條件優(yōu)越，區(qū)位條件好，可達性較強，耕地生態(tài)及景觀文化功能價值顯化較大；農(nóng)產(chǎn)品主產(chǎn)區(qū)土壤質量、氣候、地貌、耕地總量等稟賦條件好，農(nóng)戶對耕地的依賴程度高，耕地的生產(chǎn)功能和社會功能的價值顯化較大；生態(tài)功能區(qū)耕地量較少且距離市區(qū)較遠，耕地生態(tài)功能較好但其價值并未得到充分顯化[12]。基于以上分析及經(jīng)驗事實，不同主體功能區(qū)數(shù)量或結構方面耕地多功能顯化對農(nóng)業(yè)產(chǎn)出影響不同，其中重點發(fā)展區(qū)最大，其次是農(nóng)產(chǎn)品主產(chǎn)區(qū)，最后是生態(tài)功能區(qū)。據(jù)此形成假說3。

假說3：無論是數(shù)量還是結構視角，不同主體功能區(qū)耕地多功能對農(nóng)業(yè)經(jīng)濟增長影響程度存在差異，可能表現(xiàn)為：重點發(fā)展區(qū)＞農(nóng)產(chǎn)品主產(chǎn)區(qū)＞生態(tài)功能區(qū)。

3 計量模型、變量測度與數(shù)據(jù)來源

3.1 計量模型

C-D生產(chǎn)函數(shù)模型是大多數(shù)學者研究經(jīng)濟增長的邏輯起點，新古典增長模型認為，經(jīng)濟增長主要由資本與勞動兩要素內生決定，技術進步外生影響。本文考慮耕地多功能，將新古典經(jīng)濟增長模型進行擴展。參考相關研究[4]，假定某區(qū)域的農(nóng)業(yè)產(chǎn)出水平由勞動力、化肥、機械動力、耕地多功能4類因素決定，模型設定的具體形式如下：

式（1）中：Cpvit為農(nóng)業(yè)產(chǎn)出水平，用種植業(yè)產(chǎn)值（108元）衡量；Laborit為農(nóng)業(yè)勞動力投入量，用農(nóng)業(yè)從業(yè)人員（104人）衡量；Fertit為化肥資本投入，用化肥投入（t）衡量；Mechit為機械資本投入，用機械總動力（kW）衡量；M fi clit為耕地多功能，從總量和結構兩個視角進行測度，具體測度方法見下文；A為效率系數(shù)，代表廣義技術水平；euit為隨機干擾項 ，衡量未納入模型的因素對農(nóng)業(yè)經(jīng)濟增長的影響；α，β1，β2 為待估參數(shù)，分別代表勞動、化肥、機械投入對農(nóng)業(yè)經(jīng)濟增長的產(chǎn)出彈性；γ為耕地多功能對農(nóng)業(yè)經(jīng)濟增長的產(chǎn)出彈性，用來衡量耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響；t代表年份，t = 2000，2001，…，2014；i代表縣域，i = 1，2，…，70。

為避免異方差，對式（1）兩邊同時取自然對數(shù)后估計參數(shù)，估計模型如下：

式（2）為耕地多功能對農(nóng)業(yè)經(jīng)濟增長影響的實證計量模型。

本文采用Eviews 8.0軟件對序列進行相同單位根檢驗LLC和不同根單位根檢驗Fisher-ADF[24]檢驗，驗證各序列的平穩(wěn)性；采用Pedroni[25]提出的基于殘差的7種面板數(shù)據(jù)協(xié)整檢驗方法對序列進行協(xié)整檢驗；采用husman檢驗確定個體影響類型，設定正確的參數(shù)估計模型。

3.2 耕地多功能測度

3.2.1 評價指標體系的構建

耕地多功能測度的關鍵在于耕地多功能的分類及其評價指標體系的構建。本文在借鑒耕地多功能內涵與分類研究文獻的基礎上[7-9，12，26-27]，考慮實際情況，將湖北省耕地功能分為生產(chǎn)功能、社會功能、生態(tài)功能三類，構建多功能評價指標體系，并用熵權法確定權重，結果如表1所示。

生產(chǎn)功能：耕地生產(chǎn)功能指耕地生產(chǎn)活動的產(chǎn)出[7]。選取糧食單產(chǎn)水平、復種指數(shù)、耕地地均產(chǎn)值、土地墾殖率4項指標來衡量。其中，糧食單產(chǎn)直接反映耕地系統(tǒng)生產(chǎn)力大小，復種指數(shù)反映耕地利用水平，耕地地均產(chǎn)值直接反映耕地經(jīng)濟產(chǎn)出水平和反映生產(chǎn)力的高低，土地墾殖率反映縣域土地開發(fā)利用程度和耕地在所有土地利用類型結構中的比重，體現(xiàn)耕地系統(tǒng)結構和生產(chǎn)功能的穩(wěn)定性。

社會功能：耕地社會功能主要體現(xiàn)為糧食安全保障功能和就業(yè)保障功能，糧食安全保障功能是指滿足區(qū)域人口糧食需求的保障能力，社會就業(yè)保障功能指承載農(nóng)村剩余勞動力的能力[28]。選取反映糧食數(shù)量安全的人均糧食供給量、糧食作物面積占比和反映質量安全的糧食安全品質保障系數(shù)三項指標來衡量糧食安全功能；選取地均勞動力承載指數(shù)、農(nóng)業(yè)從業(yè)人員占農(nóng)業(yè)人口比重和反映農(nóng)業(yè)機械對勞動替代能力的農(nóng)業(yè)機械化水平三項指標來衡量就業(yè)保障功能。

生態(tài)功能：耕地生態(tài)功能指維護農(nóng)田生態(tài)系統(tǒng)平衡和農(nóng)村土壤健康的能力[28]。它不僅受到自身系統(tǒng)多樣性的影響，而且受到人類施用化肥等土地利用行為的外部影響。選取農(nóng)田系統(tǒng)多樣性、化肥環(huán)境污染風險指數(shù)（包括氮肥環(huán)境污染風險指數(shù)、磷肥環(huán)境污染風險指數(shù)、鉀肥環(huán)境污染風險指數(shù)）指標來進行衡量。其中，農(nóng)田系統(tǒng)多樣性參考宋小青等的研究[9]，計算公式為：

式（3）中：Hi代表農(nóng)田系統(tǒng)多樣性指數(shù)；Pi為各類農(nóng)作物播種面積占總播種面積比例；n為研究區(qū)的農(nóng)作物的種類數(shù)，結合湖北省實際情況，選取油菜、大豆、棉花 、水稻、小麥、玉米6類典型農(nóng)作物進行計算。

化肥環(huán)境污染風險指數(shù)參考劉欽普的研究[29]，計算公式如下：

式（4）中：Ri代表第i種化肥的環(huán)境污染風險指數(shù)；Fi代表單位面積耕地第i種化肥投入量；Ti為第i種化肥的安全閾值；Wi為第i種化肥污染環(huán)境的風險權重。參考劉欽普[29]的研究成果，氮肥、磷肥、鉀肥的Ti值分別設置為125 kg/hm2、62.5 kg/hm2、62.5 kg/hm2，Wi值分別設置為0.648、0.230、0.122。

3.2.2 耕地多功能測度方法

目前，總量視角耕地多功能測度方法主要有綜合指數(shù)法[27]、標準化和法[30]，綜合指數(shù)法將耕地多功能作為整體，考慮各個指標在整體中的相對重要性；標準化和法將各子功能作為整體，將各子功能標準化后的值進行加總，不將各指標納入多功能整體進行考慮。結構視角耕地多功能測度方法主要借鑒多樣性指數(shù)測度方法，包括香農(nóng)多樣性指數(shù)[31]和辛普森倒數(shù)指數(shù)[32-33]，其中香農(nóng)多樣性指數(shù)能夠較好的衡量類型的多樣性，但在反映數(shù)量均勻度等結構特征方面存在較大缺陷，而辛普森倒數(shù)指數(shù)能夠較好的刻畫耕地多功能數(shù)量的均勻度特征[33]。根據(jù)前文對總量和結構視角耕地多功能概念的界定，本文采用綜合指數(shù)法和辛普森倒數(shù)指數(shù)測度耕地多功能，計算公式如下所示：

綜合指數(shù)法：

式（5）中：Fθi為θ地區(qū)第i年的耕地總功能的得分值；n 為指標個數(shù)；Wj為第j項指標權重，采用熵權法確定；X'θij表示θ地區(qū)第i年第j個指標的標準化值。（其中，負向指標首先采用倒數(shù)法Xθij= 1/xθij進行歸一化處理，然后與正向指標統(tǒng)一按照公式X'θij= Xθij/Xmax進行標準化，x代表原始指標值，X代表歸一化后的指標值。）

辛普森倒數(shù)指數(shù)法：

式（6）中：SRIθi為θ地區(qū)第i年的辛普森倒數(shù)指數(shù)，反映耕地功能類型的多樣性與數(shù)量的均勻度，SRIθi值越大，耕地功能類型越多樣，數(shù)量越均勻；Fθij為θ地區(qū)第i年第j項耕地功能的得分值；N為研究區(qū)耕地功能的種類數(shù)。

表1 耕地多功能評價指標說明Tab.1 Multifunctional farmland indexes description

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

本文數(shù)據(jù)為湖北省70縣域2000—2014年的統(tǒng)計數(shù)據(jù)。年末耕地面積、糧食作物的播種面積、農(nóng)業(yè)從業(yè)人員數(shù)量、農(nóng)村從業(yè)人員數(shù)量、氮肥、磷肥和鉀肥施用量、農(nóng)業(yè)機械總動力、糧食作物總產(chǎn)量、農(nóng)業(yè)總產(chǎn)值等數(shù)據(jù)來源于2001—2015年《湖北省農(nóng)村統(tǒng)計年鑒》；地區(qū)總面積、種植業(yè)產(chǎn)值、各縣人口總量、、農(nóng)業(yè)人口數(shù)量等來自2001—2015年《各縣市統(tǒng)計年鑒》及《湖北省統(tǒng)計年鑒》，化肥施用安全閾值標準以及氮肥、磷肥和鉀肥環(huán)境污染權重相關數(shù)據(jù)來源于文獻參考[9，29]。個別縣市個別年份的數(shù)據(jù)缺失，以相鄰年份數(shù)據(jù)為基礎，采用移動平均、趨勢外推的方法進行數(shù)據(jù)重建；為扣除價格因素影響，所有農(nóng)業(yè)總產(chǎn)值以2000年為基期進行平減。

圖1 研究區(qū)農(nóng)業(yè)產(chǎn)出與綜合功能指數(shù)變化圖Fig.1 Changes of agricultural output and comprehensive functional index in study area

4 研究結果

4.1 總量視角耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響

根據(jù)總量視角耕地多功能的測度結果，繪制研究區(qū)和不同主體功能區(qū)耕地綜合功能指數(shù)和農(nóng)業(yè)產(chǎn)出隨時間的變化圖（圖1—圖3）。從圖1看，研究區(qū)農(nóng)業(yè)產(chǎn)出始終處于增長狀態(tài)；耕地綜合功能指數(shù)總體呈現(xiàn)“U”型的變化：2000—2007年，呈現(xiàn)下降—上升—再下降—再上升的曲折變化，但總體上處于下降狀態(tài)；2007—2014年期間也有一定波動，但總體呈現(xiàn)上升狀態(tài)；這與宋小青關于全國耕地多功能總量利用轉型的結論相似[9]，湖北省耕地功能數(shù)量利用也于2007年左右開始轉型。從圖2看，2000—2014年，不同主體功能區(qū)農(nóng)業(yè)產(chǎn)出始終處于增長狀態(tài)，產(chǎn)值及其增長速率整體表現(xiàn)為重點發(fā)展區(qū)大于農(nóng)產(chǎn)品主產(chǎn)區(qū)大于生態(tài)功能區(qū)。從圖3看，不同主體功能區(qū)綜合功能指數(shù)及其變化差異顯著：數(shù)量上，表現(xiàn)為重點發(fā)展區(qū)和農(nóng)產(chǎn)品主產(chǎn)區(qū)較大，生態(tài)功能區(qū)較小，結合前文分析可知，這是由于耕地多功能在不同主體功能區(qū)耕地多功能顯化程度的差異所致；變化趨勢上，大多數(shù)重點發(fā)展區(qū)表現(xiàn)為增加趨勢，農(nóng)產(chǎn)品主產(chǎn)區(qū)變化平穩(wěn)，生態(tài)功能區(qū)表現(xiàn)為下降趨勢。鑒于此，后文將分為2000—2007年、2007—2014年、重點發(fā)展區(qū)、農(nóng)產(chǎn)品主產(chǎn)區(qū)、生態(tài)功能區(qū)5類情況作為5個估算模型進行回歸分析，探討研究區(qū)轉型前后和不同主體功能區(qū)總量視角耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響。

圖2 湖北省不同主體功能區(qū)農(nóng)業(yè)產(chǎn)出變化圖Fig.2 Changes of agricultural output in different main functional areas in Hubei Province

圖3 湖北省不同主體功能區(qū)綜合功能指數(shù)變化圖Fig.2 Changes of comprehensive function index in different main functional areas in Hubei Province

表2 總量視角耕地多功能模型估計結果Tab.2 Estimation results of models with multifunctional farmland in quantity

5個模型的穩(wěn)定性檢驗表明所有序列在水平值檢驗時接受存在單位根的原假設，對一階差分后的各序列進行單位根檢驗發(fā)現(xiàn)，各序列在5%水平下拒絕原假設，所有序列滿足非平穩(wěn)的一階單整要求。協(xié)整性檢驗發(fā)現(xiàn)所有統(tǒng)計量均在10%的顯著水平下拒絕“不存在協(xié)整關系”原假設，表明非平穩(wěn)序列l(wèi)nY、lnL、lnK1、lnK2、lnMficl之間存在長期協(xié)整關系。husman檢驗顯示，所有相伴概率均在水平上顯著，拒絕原假設，選用固定效應模型。鑒于本文估計研究區(qū)各時段及各主體功能區(qū)耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響，不考慮截面和時間變量對參數(shù)的影響，采用固定效應變截距模型進行估計。

數(shù)量視角耕地多功能模型參數(shù)估計及檢驗結果如表2所示。從變量的回歸系數(shù)來看，所有變量值都在以下，表明解釋變量顯著；各個模型的值均通過了顯著水平為1%的檢驗，表明模型整體擬合較好；各個模型調整的R2均較大，說明建立的固定效應變截距模型與數(shù)據(jù)的相關擬合度較好。回歸結果能夠較好分析總量視角耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響。

從表2看，總量視角耕地多功能對農(nóng)業(yè)經(jīng)濟增長影響為正，與假說1相吻合。結合理論基礎分析可知，這是由于耕地多功能總量的增加會改善耕地及其他農(nóng)業(yè)生產(chǎn)要素或農(nóng)產(chǎn)品的供給數(shù)量、質量與效率，進而提高農(nóng)業(yè)生產(chǎn)效率和效益,促進農(nóng)業(yè)經(jīng)濟增長。2000—2007年和2007—2014年耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響系數(shù)分別為0.947、2.374，2007—2014年影響系數(shù)較2000—2007年要大，這是由于2007—2014年間對應的社會經(jīng)濟發(fā)展水平較2007—2010年要高，單位數(shù)量耕地多功能顯化的農(nóng)業(yè)產(chǎn)出價值更大，農(nóng)業(yè)生產(chǎn)效率與效益提升更明顯，正向影響效果更大。

重點發(fā)展區(qū)、農(nóng)產(chǎn)品主產(chǎn)區(qū)和生態(tài)功能區(qū)總量視角耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響系數(shù)分別為2.469、2.009、1.418，影響程度表現(xiàn)為重點發(fā)展區(qū)大于農(nóng)產(chǎn)品主產(chǎn)區(qū)大于生態(tài)功能區(qū)，與假說3相吻合。這是由于湖北省生態(tài)功能區(qū)和農(nóng)產(chǎn)品主產(chǎn)區(qū)的人均GDP分別只有重點發(fā)展區(qū)的1/3和1/2[11]，社會經(jīng)濟條件相對較差，并且重點發(fā)展區(qū)區(qū)位條件更優(yōu)越，居民的耕地多功能需求更多樣，單位數(shù)量耕地多功能顯化的農(nóng)業(yè)產(chǎn)出價值較高，對農(nóng)業(yè)經(jīng)濟增長的影響程度更大。總體上，總量視角耕地多功能利用對農(nóng)業(yè)產(chǎn)出的影響系數(shù)遠遠大于耕地面積對農(nóng)業(yè)產(chǎn)出的影響系數(shù)[2-4]，可能是農(nóng)業(yè)經(jīng)濟增長的新途徑。

勞動力對農(nóng)業(yè)經(jīng)濟增長的產(chǎn)出彈性為負，表明勞動力資源的再投入已經(jīng)不能成為湖北省農(nóng)業(yè)經(jīng)濟增長的源泉，湖北省農(nóng)業(yè)生產(chǎn)已經(jīng)處于勞動力“過密化”或是“內卷化”階段。化肥和機械投入對農(nóng)業(yè)經(jīng)濟增長的產(chǎn)出彈性為正，說明通過增加化肥和機械動力投入可以促進農(nóng)業(yè)經(jīng)濟的增長。

4.2 結構視角耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響

根據(jù)結構視角耕地多功能的測度結果，繪制辛普森倒數(shù)指數(shù)與農(nóng)業(yè)產(chǎn)出隨時間的變化圖（圖4—圖5）。從圖4看，研究區(qū)辛普森倒數(shù)指數(shù)與綜合功能指數(shù)的變化相反，呈倒“U”型狀態(tài)：2000—2007年間，辛普森倒數(shù)指數(shù)波動頻繁，但整體呈增加狀態(tài)；2007—2014年間，辛普森倒數(shù)指數(shù)一直呈下降狀態(tài)。結合圖1可知，結構視角耕地多功能的增加會引起耕地多功能總量減少，與理論分析吻合；類似于總量視角耕地多功能，結構視角耕地多功能利用也于2007年發(fā)生轉型。從圖5看，辛普森倒數(shù)指數(shù)及其變化在不同主體功能區(qū)差異顯著：數(shù)量上，總體表現(xiàn)為生態(tài)功能區(qū)大于農(nóng)產(chǎn)品主產(chǎn)區(qū)大于重點發(fā)展區(qū)；變化趨勢上，大多數(shù)重點發(fā)展區(qū)和農(nóng)產(chǎn)品主產(chǎn)區(qū)呈下降趨勢，生態(tài)功能區(qū)呈上升趨勢。基于以上描述與分析，后文將分研究區(qū)2000—2007年、2007—2014年、重點發(fā)展區(qū)、農(nóng)產(chǎn)品主產(chǎn)區(qū)、生態(tài)功能區(qū)5類情況作為5個估算模型進行回歸分析，探討研究區(qū)轉型前后和不同主體功能區(qū)結構視角耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響。

圖4 研究區(qū)農(nóng)業(yè)產(chǎn)出和辛普森倒數(shù)指數(shù)變化圖Fig.4 Changes of agricultural output and Simpson’s reciprocal index in the study area

圖5 湖北省不同主體功能區(qū)辛普森倒數(shù)指數(shù)變化圖Fig.5 Changes of Simpson’s reciprocal index in different main functional areas in Hubei Province

結構視角5個模型各序列也通過了平穩(wěn)性及協(xié)整性檢驗，經(jīng)huasman檢驗后同樣采用固定效應變截距模型進行估計，估計及其檢驗結果如表3所示。回歸結果同樣能夠較好分析耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響。從表3看，結構視角耕地多功能對農(nóng)業(yè)經(jīng)濟增長影響為負，與假說2相吻合。結合理論基礎分析可知，可能是由于耕地各類功能的數(shù)量越均勻，功能之間的競爭與權衡程度更高，致使耕地多功能總量減少，進而影響耕地及其他農(nóng)業(yè)生產(chǎn)要素或農(nóng)產(chǎn)品等的供給數(shù)量、質量與效率，降低農(nóng)業(yè)經(jīng)濟增長的效率和效益，對農(nóng)業(yè)經(jīng)濟增長產(chǎn)生負向影響。2000—2007年和2007—2014年耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響系數(shù)分別為-0.421、-1.058。2007—2014年間負向影響程度較2000-2007年要大，這是由于2007—2014年對應的社會經(jīng)濟發(fā)展水平較2000—2007年更高，單位結構耕地多功能引起農(nóng)業(yè)生產(chǎn)效率與效益下降更明顯，農(nóng)業(yè)產(chǎn)出損失價值更大，負向影響效果更大。

重點發(fā)展區(qū)、農(nóng)產(chǎn)品主產(chǎn)區(qū)和生態(tài)功能區(qū)結構視角耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響系數(shù)分別為-1.011、-1.135、-0.618，負向影響程度排序為農(nóng)產(chǎn)品主產(chǎn)區(qū)大于重點發(fā)展區(qū)大于生態(tài)功能區(qū)，基本符合假說3。這是由于社會經(jīng)濟發(fā)展水平越高，耕地資源稟賦越優(yōu)，單位結構耕地多功能對應的農(nóng)業(yè)產(chǎn)出價值更大；其中，農(nóng)產(chǎn)品主產(chǎn)區(qū)的負向影響程度大于重點發(fā)展區(qū)，可能的原因是耕地資源稟賦相比于社會經(jīng)濟發(fā)展水平和區(qū)位條件對農(nóng)業(yè)經(jīng)濟增長的影響效應更大。勞動、機械和資本要素對農(nóng)業(yè)經(jīng)濟增長的影響與上述數(shù)量視角類似。

表3 結構視角耕地多功能模型估計結果Tab.3 Estimation results of models with multifunctional farmland in structure

5 結論與建議

本文以湖北省為例，從總量和結構兩個視角測度耕地多功能，運用擴展的C-D生產(chǎn)函數(shù)模型分析了總量與結構視角不同時段和不同主體功能區(qū)耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響，得到以下結論。

（1）總量上的耕地多功能在2000—2014年間總體呈“U”型變化，而結構上的耕地多功能總體呈倒“U”型變化，2007年左右均開始轉型。不同主體功能區(qū)總量和結構上的耕地多功能存在顯著差異：重點發(fā)展區(qū)總量上耕地多功能較強且處于增強狀態(tài)，結構上最弱且一直減弱；農(nóng)產(chǎn)品主產(chǎn)區(qū)總量上耕地多功能較弱但變化平穩(wěn)，結構上較弱且持續(xù)減弱。生態(tài)功能區(qū)總量上的耕地多功能最弱且不斷減弱，結構上最強且不斷增強。

（2）總量視角耕地多功能對農(nóng)業(yè)經(jīng)濟增長的正向影響顯著，2007—2014年的正效應大于2000—2007年，耕地多功能增強將促進農(nóng)業(yè)經(jīng)濟的增長，但在不同的主體功能區(qū)存在顯著的區(qū)域差異，重點發(fā)展區(qū)的效應大于農(nóng)產(chǎn)品主產(chǎn)區(qū)，而生態(tài)功能區(qū)最小。

（3）結構視角耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響統(tǒng)計上顯著為負，2007—2014年的負效應大于2000—2007年，耕地多功能增強將阻礙農(nóng)業(yè)經(jīng)濟的增長，并且在不同的主體功能區(qū)亦存在顯著差異，結構視角的負效應表現(xiàn)為：農(nóng)產(chǎn)品主產(chǎn)區(qū)大于重點發(fā)展區(qū)大于生態(tài)功能區(qū)，多功能之間的多樣化協(xié)同效應還未得到發(fā)揮。

基于上述研究結論，提出以下建議：（1）注重耕地多功能總量和結構的辯證關系，當前階段以加強總量的耕地多功能為主，應在總量增長的同時，提高功能的多樣化協(xié)同，促進農(nóng)業(yè)經(jīng)濟的持續(xù)增長。（2）實施區(qū)域差異化的耕地多功能利用政策，進一步顯化各主體功能區(qū)的主導功能，發(fā)展二三產(chǎn)業(yè)，推動地區(qū)產(chǎn)業(yè)融合，提升社會經(jīng)濟發(fā)展水平。重點發(fā)展區(qū)靠近城市距消費市場近，城市居民對耕地的生態(tài)景觀文化功能需求大，地方政府可有序建設一批市民農(nóng)園、休閑農(nóng)場等進一步顯化耕地的生態(tài)景觀文化功能；農(nóng)產(chǎn)品主產(chǎn)區(qū)土壤質量好，在推動農(nóng)地規(guī)模經(jīng)營增強耕地生產(chǎn)功能的同時，推行秸稈還田、測土施肥等模式和技術，合理控制農(nóng)業(yè)勞動力數(shù)量轉移，增強耕地的生態(tài)功能和社會功能。生態(tài)功能區(qū)結構視角的耕地多功能性較強且處于增強階段，此類型區(qū)的耕地承擔了更多的非商品性功能與服務，但這種非商品性功能與服務的價值并未得到較好顯化，應加強農(nóng)地生態(tài)補償，積極培育農(nóng)業(yè)市場、耕地生態(tài)系統(tǒng)服務市場，顯化農(nóng)產(chǎn)品的“綠色”價值。（3）農(nóng)戶是耕地利用的最直接主體，應加強對農(nóng)戶的科學施肥、科學種植等技術培訓，增強其管理農(nóng)地的知識和能力以增強農(nóng)地的多功能性并促進其價值的顯化，從而推動農(nóng)業(yè)經(jīng)濟持續(xù)增長。

：

［1］洪銀興. 論經(jīng)濟增長方式轉變的基本內涵［J］. 管理世界，1999（4）：15 - 22.

［2］LIO M, LIU M C. Governance and agricultural productivity: a cross-national analysis［J］. Food Policy，2008，33（6）：504 -512.

［3］HU B, MCALEER M. Estimation of Chinese agricultural production efficiency with panel data［J］. Mathematics and Computers in Simulation，2014，68（5）：474 - 483.

［4］吳玉鳴. 中國區(qū)域農(nóng)業(yè)生產(chǎn)要素的投入產(chǎn)出彈性測算——基于空間計量經(jīng)濟模型的實證［J］. 中國農(nóng)村經(jīng)濟，2010（6）：25 - 37.

［5］郁建興，高翔. 從行政推動到內源發(fā)展：中國農(nóng)業(yè)農(nóng)村的再出發(fā)［M］. 北京：北京師范大學出版社，2013.

［6］WILSON G A. The spatiality of multifunctional agriculture: a human geography perspective［J］. Geoforum，2009，40（2）：269 - 280.

［7］姜廣輝，張鳳榮，孔祥斌，等. 耕地多功能的層次性及其多功能保護［J］. 中國土地科學，2011，25（8）：42 - 47.

［8］宋小青，歐陽竹. 耕地多功能內涵及其對耕地保護的啟示［J］. 地理科學進展，2012，31（7）：859 - 868.

［9］宋小青，吳志峰，歐陽竹. 1949年以來中國耕地功能變化［J］. 地理學報，2014，69（4）：435 - 447.

［10］楊雪，談明洪. 近年來北京市耕地多功能演變及其關聯(lián)性［J］. 自然資源學報，2014，29（5）：733 - 743.

［11］羅成，蔡銀鶯. 湖北省農(nóng)產(chǎn)品主產(chǎn)區(qū)耕地資源功能的時空演變［J］. 經(jīng)濟地理，2016，36（3）：153 - 161.

［12］施園園，趙華甫，鄖文聚，等. 北京市耕地多功能空間分異及其社會經(jīng)濟協(xié)調模式解釋［J］. 資源科學，2015，37（2）：247 - 257.

［13］胡偉艷，魏安奇，趙志尚，等. 農(nóng)地多功能供需錯位與協(xié)同作用研究進展及趨勢［J］. 中國土地科學，2017，31（3）：89 - 97.

［14］向敬偉，李江風. 貧困山區(qū)耕地利用轉型對農(nóng)業(yè)經(jīng)濟增長質量的影響［J］. 中國人口·資源與環(huán)境，2018，28（1）：71 - 81.

［15］OECD P E, SHOBAYASHI M. Multifunctionality: towards an analytical framework［M］. Paris：OECD， 2001.

［16］宋小青，歐陽竹. 中國耕地多功能管理的實踐路徑探討［J］. 自然資源學報，2012，27（4）：540 - 551.

［17］霍雅勤，蔡運龍，王瑛. 耕地對農(nóng)民的效用考察及耕地功能分析［J］. 中國人口·資源與環(huán)境，2004（3）：107 - 110.

［18］王利文. 中國北方農(nóng)牧交錯帶生態(tài)用地變化對農(nóng)業(yè)經(jīng)濟的影響分析［J］. 中國農(nóng)村經(jīng)濟，2009（4）：80 - 85.

［19］柯新利，李紅艷，劉榮霞. 武漢市耕地景觀游憩功能與可達性的空間匹配格局［J］. 長江流域資源與環(huán)境，2016，25（5）：751 - 760.

［20］劉衛(wèi)東. 耕地多功能保護問題研究［J］. 國土資源科技管理，2008（1）：1 - 5.

［21］豐雷，蔣妍，黃曉宇. 土地利用結構與社會經(jīng)濟結構——基于中國不同用地類型區(qū)域的實證分析［J］. 地域研究與開發(fā)，2011，30（1）：110 - 115.

［22］HASSAN R, SCHOLES R, ASH N. Ecosystems and human well-being: current state and trends［J］. Journal of Bacteriology，2005，1（5）：1387 - 1404.

［23］SU C H. Study on spatiotemporal change of ecosystem services and its anthropogenic driving mechanisms a case study in Yanhe Watershed［D］. Beijing: Graduate University of Chinese Academy of Sciences，2011．

［24］李子奈. 計量經(jīng)濟學（第二版）［M］. 北京：高等教育出版社，2005：360.

［25］PEDRONI P. Panel cointegration: asymptotic and finite sample properties of, pooled time series tests with an application to the PPP hypothesis［J］. Department of Economics Working Papers, 2004, 20（3）： 597 - 625.

［26］楊雪，談明洪. 北京市耕地功能空間差異及其演變［J］.地理研究，2014，33（6）：1106 - 1118.

［27］LUUK F, FILOMENA D, IRMGARD E. A conceptual framework for the assessment of multiple functions of aggro-ecosystems: a case study of Tra’s-os-Montes olive groves［J］.Journal of Rural Studies，2009，25（1）：141 - 155.

［28］陳美球，王光遠. 農(nóng)民對耕地非生產(chǎn)性功能的認識及其量化實證研究——基于1065份問卷調查［J］. 中國土地科學，2013，27（3）：10 - 16.

［29］劉欽普. 中國化肥投入?yún)^(qū)域差異及環(huán)境風險分析［J］. 中國農(nóng)業(yè)科學，2014，47（18）：3596 - 3605.

［30］RODRGUEZ-LOINAZ G, ALDAY J G, ONAINDIA M. Multiple ecosystem services landscape index: a tool for multifunctional landscapes conservation［J］. Journal of Environmental Management，2015，147：152 - 163.

［31］PLIENINGER T, BILLING C, OHNESORGE B, et al.Exploring futures of ecosystem services in cultural landscapes through participatory scenario development in the Swabian Alb, Germany［J］. Ecology and Society，2013，18（3）：261 - 272.

［32］RAUDSEPP-HEARNE C, PETERSON G D, BENNETT E M.Ecosystem service bundles for analyzing trade-offs in diverse landscapes［J］. Proceedings of the National Academy of Sciences，2010，107（11）：5242 - 5247.

［33］STU¨RCK J, VERBURG P H. Multifunctionality at what scale? a landscape multifunctionality assessment for the European Union under conditions of land use change［J］.Landscape Ecology，2017，32（3）：481 - 500.