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

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

1 引言

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

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

2 理論基礎(chǔ)與研究假說

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

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

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

假說2：結(jié)構(gòu)視角耕地多功能對(duì)農(nóng)業(yè)經(jīng)濟(jì)增長(zhǎng)可能起負(fù)向作用。

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

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

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

3.1 計(jì)量模型

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

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

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

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

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

3.2 耕地多功能測(cè)度

3.2.1 評(píng)價(jià)指標(biāo)體系的構(gòu)建

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

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

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

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

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

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

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

3.2.2 耕地多功能測(cè)度方法

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

綜合指數(shù)法：

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

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

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

表1 耕地多功能評(píng)價(jià)指標(biāo)說明Tab.1 Multifunctional farmland indexes description

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

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

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

4 研究結(jié)果

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

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

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

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

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

重點(diǎn)發(fā)展區(qū)、農(nóng)產(chǎn)品主產(chǎn)區(qū)和生態(tài)功能區(qū)總量視角耕地多功能對(duì)農(nóng)業(yè)經(jīng)濟(jì)增長(zhǎng)的影響系數(shù)分別為2.469、2.009、1.418，影響程度表現(xiàn)為重點(diǎ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分別只有重點(diǎn)發(fā)展區(qū)的1/3和1/2[11]，社會(huì)經(jīng)濟(jì)條件相對(duì)較差，并且重點(diǎn)發(fā)展區(qū)區(qū)位條件更優(yōu)越，居民的耕地多功能需求更多樣，單位數(shù)量耕地多功能顯化的農(nóng)業(yè)產(chǎn)出價(jià)值較高，對(duì)農(nóng)業(yè)經(jīng)濟(jì)增長(zhǎng)的影響程度更大?？傮w上，總量視角耕地多功能利用對(duì)農(nóng)業(yè)產(chǎn)出的影響系數(shù)遠(yuǎn)遠(yuǎn)大于耕地面積對(duì)農(nóng)業(yè)產(chǎn)出的影響系數(shù)[2-4]，可能是農(nóng)業(yè)經(jīng)濟(jì)增長(zhǎng)的新途徑。

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

4.2 結(jié)構(gòu)視角耕地多功能對(duì)農(nóng)業(yè)經(jīng)濟(jì)增長(zhǎng)的影響

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

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

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

表3 結(jié)構(gòu)視角耕地多功能模型估計(jì)結(jié)果Tab.3 Estimation results of models with multifunctional farmland in structure

5 結(jié)論與建議

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

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

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

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

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

：

［1］洪銀興. 論經(jīng)濟(jì)增長(zhǎng)方式轉(zhuǎn)變的基本內(nèi)涵［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］吳玉鳴. 中國(guó)區(qū)域農(nóng)業(yè)生產(chǎn)要素的投入產(chǎn)出彈性測(cè)算——基于空間計(jì)量經(jīng)濟(jì)模型的實(shí)證［J］. 中國(guó)農(nóng)村經(jīng)濟(jì)，2010（6）：25 - 37.

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

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

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

［8］宋小青，歐陽(yáng)竹. 耕地多功能內(nèi)涵及其對(duì)耕地保護(hù)的啟示［J］. 地理科學(xué)進(jìn)展，2012，31（7）：859 - 868.

［9］宋小青，吳志峰，歐陽(yáng)竹. 1949年以來中國(guó)耕地功能變化［J］. 地理學(xué)報(bào)，2014，69（4）：435 - 447.

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

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

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

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

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

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

［16］宋小青，歐陽(yáng)竹. 中國(guó)耕地多功能管理的實(shí)踐路徑探討［J］. 自然資源學(xué)報(bào)，2012，27（4）：540 - 551.

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

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

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

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

［21］豐雷，蔣妍，黃曉宇. 土地利用結(jié)構(gòu)與社會(huì)經(jīng)濟(jì)結(jié)構(gòu)——基于中國(guó)不同用地類型區(qū)域的實(shí)證分析［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ì)量經(jīng)濟(jì)學(xué)（第二版）［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］陳美球，王光遠(yuǎn). 農(nóng)民對(duì)耕地非生產(chǎn)性功能的認(rèn)識(shí)及其量化實(shí)證研究——基于1065份問卷調(diào)查［J］. 中國(guó)土地科學(xué)，2013，27（3）：10 - 16.

［29］劉欽普. 中國(guó)化肥投入?yún)^(qū)域差異及環(huán)境風(fēng)險(xiǎn)分析［J］. 中國(guó)農(nóng)業(yè)科學(xué)，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.