非糧化對糧食綠色全要素生產(chǎn)率的門檻效應(yīng)研究*

薛選登，谷秀云

（河南科技大學(xué)經(jīng)濟學(xué)院，洛陽 471000）

0 引言

推進(jìn)農(nóng)業(yè)綠色發(fā)展是農(nóng)業(yè)發(fā)展觀的一場深刻革命。農(nóng)業(yè)農(nóng)村部印發(fā)《農(nóng)業(yè)綠色發(fā)展技術(shù)導(dǎo)則（2018-2030年）》中指出農(nóng)業(yè)生產(chǎn)的“三個轉(zhuǎn)變”，即：從注重數(shù)量為主向數(shù)量質(zhì)量效益并重轉(zhuǎn)變，從注重生產(chǎn)功能為主向生產(chǎn)生態(tài)功能并重轉(zhuǎn)變，從注重單要素生產(chǎn)率提高為主向全要素生產(chǎn)率提高為主轉(zhuǎn)變[1]，這表明我國對于農(nóng)業(yè)生產(chǎn)既注重糧食安全，更兼顧農(nóng)業(yè)發(fā)展質(zhì)量。因此，想要改善化肥、農(nóng)藥、農(nóng)膜、除草劑等傳統(tǒng)投入要素所引致的環(huán)境污染和溫室氣體排放嚴(yán)重的現(xiàn)象，提高農(nóng)業(yè)綠色全要素生產(chǎn)率是時代所需。但與穩(wěn)定實現(xiàn)農(nóng)業(yè)綠色發(fā)展目標(biāo)相悖的是，交織農(nóng)業(yè)生產(chǎn)成本上漲，種糧利潤空間受擠，理性農(nóng)戶出于自身利益考慮，部分地區(qū)的耕地非糧化現(xiàn)象愈發(fā)凸顯，那么耕地非糧化對農(nóng)業(yè)綠色生產(chǎn)是否造成影響？影響“度”如何？不同區(qū)域有無差異？文章以糧食生產(chǎn)為例，實證探究非糧化對糧食綠色全要素生產(chǎn)率的影響。

糧食綠色全要素生產(chǎn)率（Green Total Factor Productivity，GTFP）較傳統(tǒng)全要素生產(chǎn)率（Total Factor Productivity，TFP）來講，可以有效地兼顧資源與環(huán)境約束，反映糧食生產(chǎn)的真實效率。對比以往綠色全要素生產(chǎn)率的相關(guān)文獻(xiàn)可以發(fā)現(xiàn)[2-7]，學(xué)者們多從以下角度開展研究，即界定綠色全要素生產(chǎn)率的相關(guān)概念、對比全要素生產(chǎn)率與綠色全要素生產(chǎn)率的測算方法、拓展研究宏微觀層面數(shù)據(jù)、分析具體的影響因素等，在學(xué)習(xí)吸收已有眾多文獻(xiàn)的基礎(chǔ)上，該文嘗試通過剖析具體影響因素，增加該部分研究的邊際貢獻(xiàn)。

很多學(xué)者在相關(guān)研究中都直接或間接地認(rèn)為“非糧化”是作為農(nóng)業(yè)綠色生產(chǎn)的一個重要的影響因素，立足農(nóng)業(yè)結(jié)構(gòu)、生態(tài)效益和生產(chǎn)質(zhì)量等角度進(jìn)行探討。常偉和馬詩雨[8]提出“非糧化”現(xiàn)象是指在農(nóng)地規(guī)模流轉(zhuǎn)過程中，轉(zhuǎn)入方將轉(zhuǎn)入的農(nóng)地大部分用于非糧食作物的種植或從事非農(nóng)經(jīng)營的現(xiàn)象，而種糧成本高、效益低等是導(dǎo)致“非糧化”的原因。陳衛(wèi)洪和漆雁斌[9]研究發(fā)現(xiàn)稻谷種植面積和大牲畜飼養(yǎng)數(shù)量對減少農(nóng)業(yè)源甲烷的排放具有重要作用，表明不同農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)形式確實對生態(tài)環(huán)境存在影響，因此考慮到農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)轉(zhuǎn)變也可能對農(nóng)業(yè)綠色全要素生產(chǎn)率產(chǎn)生影響。而金芳[10]研究發(fā)現(xiàn)新時期確應(yīng)重視農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)的合理調(diào)整在農(nóng)業(yè)綠色生產(chǎn)率增長中的重要作用，以提升農(nóng)業(yè)整體發(fā)展質(zhì)量。宋戈、白小艷、高佳[11]綜合運用負(fù)外部效益測算方法，發(fā)現(xiàn)研究區(qū)耕地非糧化糧食安全負(fù)外部效益和耕地非糧化生態(tài)環(huán)境負(fù)外部效益的空間分異規(guī)律基本一致。郭珍[12]提出隱性的耕地功能性流失，如優(yōu)質(zhì)耕地偽生態(tài)化、非糧化、過度利用化與耕地邊際化等，對耕地數(shù)量、質(zhì)量、生態(tài)的影響不容忽視。文獻(xiàn)表明，耕地非糧化產(chǎn)生的主要原因是由于生產(chǎn)成本高，收益低，而耕地非糧化所帶來的農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)形式轉(zhuǎn)變，加之耕地非糧化對生態(tài)環(huán)境的潛在外部性，將很大程度上對糧食綠色全要素生產(chǎn)率產(chǎn)生影響。

綜上所述，現(xiàn)有文獻(xiàn)存在一些有待完善之處：一是缺乏非糧化與糧食綠色全要素生產(chǎn)率的直接相關(guān)研究。以往研究對“非糧化”現(xiàn)象與農(nóng)業(yè)綠色生產(chǎn)之間的關(guān)系探討多圍繞各自角度出發(fā)，對兩者之間的作用機制仍待挖掘探討。二是缺乏考慮區(qū)域異質(zhì)性條件影響。由于未綜合考慮不同地區(qū)經(jīng)濟發(fā)展水平、資源稟賦差異等擾動條件，區(qū)域同質(zhì)性假設(shè)所得結(jié)論并不能很好的指導(dǎo)實踐，尤其是沒有重視農(nóng)業(yè)決策主體所受種糧成本收益的各自差異，所帶給區(qū)域非糧化趨勢增減的重要影響。鑒于此，該文利用2008—2018年26個省（市、自治區(qū)）份的面板數(shù)據(jù)，在測算綜合考慮期望產(chǎn)出與非期望產(chǎn)出的糧食綠色全要素生產(chǎn)率的基礎(chǔ)上，實證考察非糧化對糧食綠色全要素生產(chǎn)率的作用機制，是否存在以種糧利潤為門檻的門檻效應(yīng)？若受門檻值影響，驅(qū)動抑或阻礙？

1 非糧化對糧食綠色全要素生產(chǎn)率影響的機理分析

前蘇聯(lián)教科書中曾提到“粗放型增長”和“集約型增長”，兩種增長方式相較而言，前者依靠要素投入驅(qū)動，后者更加注重生產(chǎn)效率與效益驅(qū)動，即全要素生產(chǎn)率的增長。1980年在索洛模型單一強調(diào)勞動和資本要素的基礎(chǔ)上，新增長理論拓展提出將知識和人力資本等作為模型中的獨立要素納入考量，而這一理論創(chuàng)新也更加貼合實際。總而言之，農(nóng)業(yè)生產(chǎn)需要綜合考慮多類投入要素，通過改善投入產(chǎn)出關(guān)系實現(xiàn)農(nóng)業(yè)全要素生產(chǎn)率增長，然而非糧化的發(fā)展會影響要素投入、要素組合以及利用效率。因此，該文選擇從要素投入結(jié)構(gòu)、種植結(jié)構(gòu)兩個視角出發(fā)，分析非糧化對糧食綠色全要素生產(chǎn)率的作用機制。

非糧化改變農(nóng)業(yè)要素投入結(jié)構(gòu)，進(jìn)而影響糧食綠色全要素生產(chǎn)率?，F(xiàn)代市場經(jīng)濟條件下，諸如勞動力、資本等生產(chǎn)要素均會受利益驅(qū)使優(yōu)化資源配置，從低效率或低收益的部門自由流動至高效率或高收益的部門。在農(nóng)業(yè)生產(chǎn)過程中，受要素流動與重新配置的影響，糧食綠色全要素生產(chǎn)率在某種程度上也將發(fā)生變化。近年來，耕地非糧化現(xiàn)象凸顯，大量原本種植糧食的農(nóng)民，由于后受種糧利潤或者土地流轉(zhuǎn)成本的影響，轉(zhuǎn)出土地后選擇從事瓜、蔬、花、果等非糧生產(chǎn)，農(nóng)業(yè)生產(chǎn)的要素投入結(jié)構(gòu)也隨即產(chǎn)生變化。一是蔬菜、水果等非糧化產(chǎn)業(yè)一般為勞動密集型，吸納勞動力較多，相較于高度機械化發(fā)展的糧食產(chǎn)業(yè)而言，該非糧產(chǎn)業(yè)所產(chǎn)生的碳排放將得到一定程度緩解。二是從土地生產(chǎn)解放出來的青年勞動力擇業(yè)進(jìn)入城市務(wù)工，留守種植糧食的勞動力結(jié)構(gòu)將呈現(xiàn)老齡化等特征，且已有相關(guān)研究指出伴隨農(nóng)業(yè)勞動力老齡化程度的加深，糧食綠色全要素生產(chǎn)率在前期會呈現(xiàn)下降趨勢[13]。三是蔬菜、水果等非糧化產(chǎn)業(yè)較之糧食種植而言，其周期性大大縮短，對于化肥、農(nóng)藥、種子、農(nóng)膜等多類農(nóng)資的投入需求將直線上升，生產(chǎn)中產(chǎn)生更多的碳排放，且耕地非糧化的高收益會促使理性小農(nóng)做出自身利益最大化的決策，產(chǎn)生輻射帶動效應(yīng)，加劇周邊地區(qū)“非糧化”，進(jìn)一步影響糧食綠色全要素生產(chǎn)率。

非糧化改變農(nóng)業(yè)種植結(jié)構(gòu)，進(jìn)而影響糧食綠色全要素生產(chǎn)率。耕地非糧化所帶來農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)的調(diào)整，往往伴隨產(chǎn)業(yè)內(nèi)部不同細(xì)分部門和要素利用效率的變化，這均會影響農(nóng)業(yè)碳排放等非期望產(chǎn)出，進(jìn)而引起糧食綠色全要素生產(chǎn)率的變化。一是耕地非糧化后，土地轉(zhuǎn)入經(jīng)營主體將原有種糧耕地改為觀光農(nóng)業(yè)、畜牧養(yǎng)殖業(yè)等，此時田間基礎(chǔ)設(shè)施將被整改，如水利灌溉、河道遷改、道路修設(shè)等，土壤層便會因此遭到破壞。一旦流轉(zhuǎn)期限到期或承包戶因經(jīng)營不善等原因中途離開，后續(xù)土地修復(fù)工作將十分困難，甚至無法再復(fù)耕，進(jìn)而對糧食綠色全要素生產(chǎn)率增長產(chǎn)生長久不利影響[14]。二是耕地資源的稀缺性表明，伴隨耕地非糧化現(xiàn)象凸顯，種糧耕地面積將進(jìn)一步縮減。而耕地面積與種糧勞動力同時減少，將不利于農(nóng)業(yè)發(fā)展適度規(guī)模生產(chǎn)，從另一方面來看，這將倒逼農(nóng)民在有限土地上更高效地投入資金、技術(shù)和管理等，以控制成本，提高產(chǎn)量，獲得生產(chǎn)效益，間接地提高糧食綠色全要素生產(chǎn)率。

由此可見，非糧化對糧食綠色全要素生產(chǎn)率的影響富有雙面性，既有負(fù)向的阻礙作用，又有正向拉動作用，且其作用程度與方向很可能受區(qū)域異質(zhì)性條件如種糧利潤的影響。非糧化對糧食綠色全要素生產(chǎn)率的影響存在以種糧利潤為門檻值的門檻效應(yīng)，這是該研究提出的重要理論假設(shè)。

2 研究方法

2.1 面板門檻模型設(shè)計

為驗證種糧利潤是否擾動非糧化對糧食綠色全要素生產(chǎn)率的作用機制，根據(jù)Hansen[15]，以種糧利潤為門檻變量，以非糧化作為門檻依賴變量，設(shè)計面板門檻模型探究非糧化對糧食綠色全要素生產(chǎn)率的門檻效應(yīng)。

設(shè)定單一面板門檻模型為:

雙重面板門檻模型為

模型（1）（2）包括被解釋變量Yit，個體固定效應(yīng)ui，影響系數(shù)β和θ，門檻依賴變量Xit，指標(biāo)函數(shù)I(·)，門檻變量qit，門檻值γ，控制變量zit，隨機擾動項εit，樣本i，及時間t。

首先需確定門檻值γ的大小，借助Hansen[16]的GridSearch方法將門檻依賴變量qit按最高和最低1%的去除，然后連續(xù)對備選門檻值進(jìn)行格柵化處理，并分別計算回歸模型的殘差平方和S1(γ)，以S1(γ)最小為判斷依據(jù)確定門檻值。最后，對門檻效應(yīng)顯著性和門檻值真實性進(jìn)行檢驗。門檻效應(yīng)顯著性檢驗采用wald檢驗方法。單一門檻模型原假設(shè):H0：β1=β2，如果P＜0.05，則認(rèn)為門檻效應(yīng)顯著。門檻值真實性檢驗參照Hansen[15]以最大似然法進(jìn)行。

2.2 變量選取及測算方法

2.2.1 被解釋變量

（1）糧食綠色全要素生產(chǎn)率（GTFP）測算指標(biāo)。①投入指標(biāo)。包括第一產(chǎn)業(yè)從業(yè)人員、糧食播種面積、農(nóng)業(yè)機械總動力、農(nóng)藥使用量、化肥折純量和有效灌溉面積。參照閔銳和李谷成[17]，通過權(quán)重系數(shù)將糧食生產(chǎn)投入從廣義的農(nóng)業(yè)投入中剝離，權(quán)重系數(shù)A=糧食播種面積/農(nóng)作物播種面積，權(quán)重系數(shù)B=(農(nóng)業(yè)產(chǎn)值/農(nóng)林牧漁總產(chǎn)值)×A。②產(chǎn)出指標(biāo)。糧食產(chǎn)量為期望產(chǎn)出，碳排放量為非期望產(chǎn)出。

（2）糧食綠色全要素生產(chǎn)率（GTFP）測算方法。SBM-GML指數(shù)法是將Global Malmquist-Luenberger生產(chǎn)率指數(shù)法(簡稱GML指數(shù)法)與SBM模型相結(jié)合，通過非徑向非角度的當(dāng)期和全局方向性距離函數(shù)來處理非合意產(chǎn)出的一種方法[18]。

2.2.2 門檻依賴變量

非糧化（flh）該文結(jié)合相關(guān)研究文獻(xiàn)與非糧化內(nèi)涵，選取農(nóng)作物總播種面積與糧食作物播種面積，計算得到非糧作物播種面積，參考閆建偉[19]定義非糧化=非糧作物播種面積/糧食作物播種面積，統(tǒng)籌兼顧糧食種植面積與非糧作物種植規(guī)模的變化。

2.2.3 門檻變量

種糧利潤(lr)該文以種糧收益（稻谷、玉米、小麥）反映歷年各地農(nóng)戶種糧利潤差異。根據(jù)前文假設(shè)可見，在不同種糧利潤的干擾下，非糧化對糧食綠色全要素生產(chǎn)率的影響呈非線性。

2.2.4 控制變量

（1）農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)（jg）。利用種植業(yè)總產(chǎn)值/農(nóng)林牧漁總產(chǎn)值表示，各地區(qū)農(nóng)業(yè)生產(chǎn)資源稟賦條件不同，既是優(yōu)勢也為約束，受經(jīng)濟發(fā)展水平與地形地勢條件的影響，因地制宜選擇農(nóng)業(yè)生產(chǎn)方式，因此，不同的農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)對糧食綠色全要素生產(chǎn)率產(chǎn)生不同的影響。（2）土地質(zhì)量水平（tz）。利用化肥施用量/農(nóng)作物種植面積表示，化肥使用可以提高糧食生產(chǎn)率，而過度使用必然導(dǎo)致邊際效果遞減，對土地質(zhì)量水平產(chǎn)生不利影響。（3）財政支農(nóng)力度（cz）。利用農(nóng)林水事務(wù)支出/財政一般預(yù)算支出表示，可以更好地衡量地區(qū)政府行為對于當(dāng)?shù)剞r(nóng)業(yè)生產(chǎn)環(huán)境的支持力度，不同的農(nóng)業(yè)投入對應(yīng)糧食生產(chǎn)基礎(chǔ)設(shè)施建設(shè)效果不同，進(jìn)而對糧食綠色全要素生產(chǎn)率影響程度不同。（4）自然環(huán)境（sz）。利用受災(zāi)面積/總播種面積表示，以此衡量農(nóng)業(yè)生產(chǎn)的自然條件。農(nóng)業(yè)生產(chǎn)天然存在弱質(zhì)性，使得自然條件的變化對糧食產(chǎn)量有直接影響。

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

鑒于數(shù)據(jù)的科學(xué)性、可獲得性及DEA對異常數(shù)據(jù)的敏感性，該文選取我國2008—2018年26個省（市、自治區(qū)）份的具體投入與產(chǎn)出面板數(shù)據(jù)，所需原始數(shù)據(jù)均為官方統(tǒng)計，摘錄于《中國統(tǒng)計年鑒》《全國農(nóng)產(chǎn)品成本收益資料匯編》等。

3 結(jié)果與分析

3.1 我國糧食綠色全要素生產(chǎn)率整體分析

從時序演變方面來看，我國糧食綠色全要素生產(chǎn)率在樣本期內(nèi)呈現(xiàn)波動上升趨勢，其分解指數(shù)技術(shù)進(jìn)步是糧食綠色全要素生產(chǎn)率增長的主要驅(qū)動力，如圖1所示。分階段來看，糧食綠色全要素生產(chǎn)率于2012—2014年首次出現(xiàn)下降，考慮到主要是因為技術(shù)進(jìn)步和技術(shù)效率進(jìn)入上升的疲軟期，尤其前沿技術(shù)進(jìn)步下降明顯。2015—2017年技術(shù)效率出現(xiàn)樣本期內(nèi)最低值，導(dǎo)致糧食綠色全要素生產(chǎn)率再次出現(xiàn)下降。隨著農(nóng)業(yè)科技在糧食生產(chǎn)中的不斷應(yīng)用，農(nóng)業(yè)生產(chǎn)技術(shù)革新比生產(chǎn)技術(shù)效率改進(jìn)在短期內(nèi)更具有效益，加之長期以來，農(nóng)業(yè)改革主要側(cè)重于農(nóng)業(yè)管理和制度層面，從而對生產(chǎn)技術(shù)效率的革新產(chǎn)生抑制，說明我國糧食依舊處于一種重要素投入而忽視效率改進(jìn)的階段，發(fā)展綠色集約型農(nóng)業(yè)亟需對此改進(jìn)。

圖1 2008—2018年我國糧食綠色全要素生產(chǎn)率及其分解指數(shù)時序演變特征

從各省域空間差異來看，2008—2018年各地糧食綠色全要素生產(chǎn)率均實現(xiàn)破“1”增長。其中，新疆居首位，為1.431，即樣本期內(nèi)其糧食綠色全要素生產(chǎn)率以年均43.09%的速度增長。廣東省最低，僅為1.009，考慮資源環(huán)境約束下的糧食綠色全要素生產(chǎn)率年均增速僅為0.95%。在26個?。ㄊ小⒆灾螀^(qū)）份中，只有新疆、內(nèi)蒙古、浙江省、山東省、湖北省等地的糧食綠色全要素生產(chǎn)率是在技術(shù)進(jìn)步和效率改進(jìn)雙重驅(qū)動下進(jìn)步的，屬于集約型增長的發(fā)展方式。而廣西、四川省、江蘇省、廣東省等10個省區(qū)基本都是依靠技術(shù)進(jìn)步拉動，說明這些地區(qū)在糧食生產(chǎn)過程中仍存在技術(shù)效率不高、資源消耗嚴(yán)重問題，今后必須進(jìn)一步重視效率的提高，兼顧技術(shù)革新和效率改進(jìn)。

從分區(qū)空間差異來看，糧食主銷區(qū)相較于糧食主產(chǎn)區(qū)、產(chǎn)銷平衡區(qū)以及全國平均值而言，其糧食綠色全要素生產(chǎn)率增長更為顯著。經(jīng)濟發(fā)達(dá)的糧食主銷區(qū)主要是由于稟賦差異，如勞動力成本居高，更側(cè)重于技術(shù)效率方面改進(jìn)，其技術(shù)效率帶來的糧食綠色全要素生產(chǎn)率增長驅(qū)動力遠(yuǎn)高于糧食主產(chǎn)區(qū)和產(chǎn)銷平衡區(qū)。相較而言，地形多屬平原的糧食主產(chǎn)區(qū)，資源稟賦較好，在發(fā)展糧食綠色生產(chǎn)中更注重農(nóng)業(yè)綠色生產(chǎn)技術(shù)的推廣應(yīng)用；而產(chǎn)銷平衡區(qū)多位于中西部高原山區(qū)，自然環(huán)境較為惡劣，不利于糧食作物生長，更大程度上需要依靠農(nóng)藥、化肥提升產(chǎn)量，便在一定程度上對于效率改進(jìn)產(chǎn)生抑制作用。

3.2 變量描述性統(tǒng)計

所有變量的基本情況匯總結(jié)果見表1所示。

表1 變量描述性統(tǒng)計

3.3 變量平穩(wěn)性檢驗

所有變量先后通過同根檢驗法LLC和異根檢驗法IPS，除糧食綠色全要素生產(chǎn)率和土地質(zhì)量水平一階差分后平穩(wěn)，其余變量均直接通過單位根檢驗，所有變量平穩(wěn)，結(jié)果見表2所示。

表2 單位根檢驗

3.4 門檻效應(yīng)檢驗及門檻值確定

為排除異方差影響、增強可比性，對所有變量取對數(shù)處理。為檢驗非糧化對糧食綠色全要素生產(chǎn)率是否存在以種糧利潤為門檻的門檻效應(yīng)，先模擬設(shè)定雙重門檻模型，采用Bootstrap循環(huán)檢驗500次，P值未通過顯著性水平檢驗，然后將模型模擬設(shè)定為單一門檻，此時門檻效應(yīng)顯著，見表3所示，故選用單一門檻模型。然后根據(jù)Hansen[15]的LR檢驗確定種糧利潤（lr）單一門檻值為7.638，見表4所示。

表3 門檻效應(yīng)檢驗

表4 門檻值及置信區(qū)間

3.5 模型穩(wěn)健性檢驗及結(jié)果分析

該文通過對比回歸，以檢驗?zāi)Ｐ徒Y(jié)果的穩(wěn)健性，見表5所示。第一個模型為不考慮控制變量的單一門檻回歸模型。第二個模型加入控制變量。第三個模型為固定效應(yīng)模型。從第一個模型回歸結(jié)果來看，當(dāng)種糧利潤水平小于7.638時，非糧化對糧食綠色全要素生產(chǎn)率的影響系數(shù)為-0.206，但并不顯著，當(dāng)種糧利潤水平大于7.638時，影響系數(shù)提高到-0.331，且具有顯著負(fù)向抑制糧食綠色全要素生產(chǎn)率的作用。第二個模型回歸結(jié)果顯示，加入控制變量后非糧化對糧食綠色全要素生產(chǎn)率的抑制程度進(jìn)一步提升，當(dāng)種糧利潤水平大于7.638時，非糧化在5%的顯著性水平上對糧食綠色全要素生產(chǎn)率呈負(fù)向抑制作用。第三個模型經(jīng)Hausman檢驗后設(shè)定為雙固定效應(yīng)模型，回歸結(jié)果再次顯示非糧化對糧食綠色全要素生產(chǎn)率呈顯著抑制作用，且抑制程度有所提升；控制變量的影響和第二個模型的影響程度和方向基本一致。由此可見，第二個模型具有一定的穩(wěn)健性，非糧化與糧食綠色全要素生產(chǎn)率之間為反向關(guān)系，且確實存在以種糧利潤門檻的門檻效應(yīng)，當(dāng)種糧利潤水平大于7.638時，這種反向抑制作用將大幅顯著提升。因此，重點分析加入控制變量的單一面板門檻模型。

表5 穩(wěn)健性檢驗

在第二個模型中，當(dāng)種糧利潤水平小于7.638時，估計系數(shù)為-0.240，但此時回歸結(jié)果并不顯著，說明非糧化對糧食綠色全要素生產(chǎn)率的負(fù)向抑制作用還沒有凸顯，一旦跨過門檻值，非糧化每提高1%，糧食綠色全要素生產(chǎn)率將顯著降低0.348%，由此可見，非糧化與糧食綠色全要素生產(chǎn)率之間受制于種糧利潤。其內(nèi)在原因在于，如果某地區(qū)種糧利潤水平較低，那么該地區(qū)農(nóng)村勞動力的種糧意愿會受到影響，在土地流轉(zhuǎn)過程中會綜合考慮土地流轉(zhuǎn)價格以及種糧成本收益，而選擇非糧種植或者從事非農(nóng)活動，從而抑制糧食全要素生產(chǎn)率增長。反之，如果某地區(qū)種糧利潤較高邁過門檻值后，那么該地區(qū)的農(nóng)村勞動力可以更好地從事糧食種植和規(guī)?；r(nóng)業(yè)生產(chǎn)，此時非糧化每降低1%，糧食綠色全要素生產(chǎn)率較以往將顯著大幅增長。此外，財政支農(nóng)力度的影響為顯著正向拉動，每提高1%，糧食綠色全要素生產(chǎn)率相應(yīng)提高0.385%，是值得關(guān)注的影響因素，表明地區(qū)財政支農(nóng)等相關(guān)補貼政策進(jìn)一步優(yōu)化農(nóng)業(yè)生產(chǎn)環(huán)境，對糧食綠色全要素生產(chǎn)率的提高起到積極作用。自然環(huán)境狀況的影響為顯著負(fù)向抑制，說明考察期內(nèi)各地區(qū)的糧食生產(chǎn)在一定程度上仍受自然災(zāi)害的制約，進(jìn)一步加強農(nóng)業(yè)基礎(chǔ)設(shè)施建設(shè)與農(nóng)業(yè)保險政策等預(yù)警性措施尤為重要。

3.6 耦合度分析

非糧化對糧食綠色全要素生產(chǎn)率的負(fù)向作用受種糧利潤影響，而非糧化與種糧利潤間互相影響，它們之間的協(xié)調(diào)度會是如何？為此該文構(gòu)建非糧化與種糧利潤之間的耦合協(xié)調(diào)度模型[20]。為消除量綱差異，對非糧化和種糧利潤標(biāo)準(zhǔn)化處理得到U1和U2的值。

3.6.1 耦合度評價模型

根據(jù)物理學(xué)上的耦合思想構(gòu)建非糧化與種糧利潤的耦合度評價模型為：

式（3）中，C代表耦合度，U1和U2分別代表標(biāo)準(zhǔn)化處理后的非糧化(flh)和種糧利潤(lr)。C值越接近0，耦合度越小，越接近1，耦合度越高。但耦合的性質(zhì)無法通過耦合度的數(shù)值來判斷，兩個子系統(tǒng)的綜合水平都低也會出現(xiàn)耦合度較高的情況，因此需要評價兩個系統(tǒng)的耦合協(xié)調(diào)度情況。

3.6.2 耦合協(xié)調(diào)度模型

耦合協(xié)調(diào)度可以綜合衡量耦合協(xié)調(diào)互動水平，構(gòu)建耦合協(xié)調(diào)度模型為：

式（4）中，D代表非糧化與種糧利潤的耦合協(xié)調(diào)度程度，C表示耦合度，T指代綜合協(xié)調(diào)度指數(shù)，α和β是兩個系統(tǒng)的權(quán)重系數(shù)，經(jīng)咨詢專家，令α=β=0.5。最終得到的耦合協(xié)調(diào)度指數(shù)，并把2008—2018年26個?。ㄊ?、自治區(qū)）份糧食綠色全要素生產(chǎn)率的平均值制作在同一表格中，見表6所示。

表6 各區(qū)域省份糧食綠色全要素生產(chǎn)率與耦合協(xié)調(diào)度對比變化

為了更清晰地展示分區(qū)域?qū)Ρ茸兓厔荩刂谱骷Z食綠色全要素生產(chǎn)率與耦合協(xié)調(diào)度的分區(qū)域?qū)Ρ茸兓厔輬D，如圖2所示。在減少黑龍江省與甘肅省兩個極值省份的影響后，可以更明顯看出非糧化與種糧利潤的耦合協(xié)調(diào)度與糧食綠色全要素生產(chǎn)率變動具有一致性，隨著各區(qū)域耦合協(xié)調(diào)度值的降低，其對應(yīng)省份糧食綠色全要素生產(chǎn)率也呈下降趨勢，即受種糧利潤干擾導(dǎo)致非糧化對糧食綠色全要素生產(chǎn)率的影響呈非線性。

圖2 糧食綠色全要素生產(chǎn)率與耦合協(xié)調(diào)度的分區(qū)域?qū)Ρ茸兓厔?/p>

4 結(jié)論

4.1 研究結(jié)果

該文首先針對非糧化與糧食綠色全要素生產(chǎn)率之間的作用機理進(jìn)行詳細(xì)論述，其次通過收集2008—2018年我國26個?。ㄊ?、自治區(qū)）份面板數(shù)據(jù)，構(gòu)建面板門檻模型，實證探究非糧化對糧食綠色全要素生產(chǎn)率影響是否存在以種糧利潤為門檻的門檻效應(yīng)，主要研究結(jié)果如下。

（1）2008—2018年全國及其分區(qū)域的糧食綠色全要素生產(chǎn)率皆呈波動上升趨勢，且糧食主銷區(qū)的綠色全要素生產(chǎn)率增長最為顯著。其中，全國及糧食主產(chǎn)區(qū)、產(chǎn)銷平衡區(qū)主要依靠技術(shù)進(jìn)步，而糧食主銷區(qū)則是受技術(shù)效率驅(qū)動，這說明技術(shù)效率在一定程度上彌補了當(dāng)?shù)丶夹g(shù)進(jìn)步方面的動力不足。

（2）種糧利潤作為非糧化影響糧食綠色全要素生產(chǎn)率的干擾因素，門檻值為7.638。當(dāng)種糧利潤不超過門檻值時，非糧化對糧食綠色全要素生產(chǎn)率的抑制作用并不凸顯；一旦跨過門檻值后，非糧化的趨勢增長將對糧食綠色全要素生產(chǎn)率產(chǎn)生顯著的負(fù)向阻礙作用，即非糧化的差別控制可以很大程度上保護(hù)糧食綠色全要素生產(chǎn)率。

（3）分區(qū)域來看，各省份非糧化與種糧利潤之間的耦合協(xié)調(diào)度，雖在數(shù)值上存在一定差異，但其變化趨勢與糧食綠色全要素生產(chǎn)率的變動趨勢基本一致。因此，在通過差別控制非糧化促進(jìn)糧食綠色全要素生產(chǎn)率提升的過程中，須將種糧利潤協(xié)同保護(hù)放在重要位置，效果才能更加顯著。

4.2 建議

根據(jù)該研究的結(jié)果，非糧化作為糧食綠色全要素生產(chǎn)率的重要影響因素，且這種影響存在一個“度”（即門檻）；另鑒于不同區(qū)域各省份非糧化與種糧利潤之間的耦合協(xié)調(diào)度與糧食綠色全要素生產(chǎn)率密切相關(guān)，故不同區(qū)域的農(nóng)業(yè)綠色發(fā)展思路將存在差異。

（1）差別控制非糧化，低門檻與高門檻地區(qū)采取不同應(yīng)對措施。其一，低門檻地區(qū)可建立非糧化預(yù)警機制，進(jìn)行動態(tài)監(jiān)測。典型的低門檻地區(qū)為產(chǎn)銷平衡區(qū)和糧食主銷區(qū)，各自可根據(jù)經(jīng)濟與農(nóng)業(yè)發(fā)展?fàn)顩r、資源稟賦條件如地形地勢等，科學(xué)合理設(shè)定非糧水平彈性范圍，通過事前監(jiān)控、事中管理、事后解決的管控模式，及時監(jiān)測非糧化動態(tài)，確保糧食耕種面積在合理范圍內(nèi)即可。其二，高門檻地區(qū)的非糧化應(yīng)嚴(yán)格控制，切實保護(hù)糧食綠色全要素生產(chǎn)率。典型的高門檻省份包括山東、河南、江蘇、湖北、安徽等，由糧食主產(chǎn)區(qū)幾乎全部覆蓋。但整體而言，糧食主產(chǎn)區(qū)實際上是高門檻與低門檻省份并存的狀態(tài)。因此，在糧食主產(chǎn)區(qū)可根據(jù)實際發(fā)展需要進(jìn)一步細(xì)分糧食功能區(qū)，如糧食種植專用地和糧經(jīng)兩用地。即在嚴(yán)格控制高門檻省份非糧化趨勢、充分確保其糧食種植專用地面積的基礎(chǔ)上，部分低門檻省份可以根據(jù)自身資源稟賦和農(nóng)業(yè)發(fā)展程度同步建立糧經(jīng)兩用地，以期最大程度上提升糧食綠色全要素生產(chǎn)率。

（2）協(xié)同保護(hù)種糧利潤，提升非糧化與種糧利潤之間的耦合協(xié)調(diào)度。分區(qū)域來看，產(chǎn)銷平衡區(qū)耦合協(xié)調(diào)度在0.5以下的4個省份：貴州、山西、甘肅、云南；糧食主產(chǎn)區(qū)耦合協(xié)調(diào)度在0.5以下的2個省份：黑龍江、吉林；糧食主銷區(qū)耦合協(xié)調(diào)度整體表現(xiàn)最好。但整體來看，耦合協(xié)調(diào)度在0.5以下的6個省份存在明顯共性特征，即低門檻省份，且糧食綠色全要素生產(chǎn)率排序靠后。以產(chǎn)銷平衡區(qū)為例，往往農(nóng)業(yè)具規(guī)模，但其糧食生產(chǎn)效率和效益低，二三產(chǎn)業(yè)反哺受到限制，且在種糧成本上行壓力下，農(nóng)戶種糧利潤受到擠壓，需要一定的非糧化來改善收入。簡言之，低門檻省份糧食綠色全要素生產(chǎn)率提升的過程中，更應(yīng)協(xié)同保護(hù)種糧利潤。其一，國家要加大實施綜合性收入補貼、生產(chǎn)補貼等支農(nóng)政策力度，合理引導(dǎo)農(nóng)戶價格預(yù)期，保障小農(nóng)生產(chǎn)效益。其二，依舊不斷完善農(nóng)村農(nóng)業(yè)公共基礎(chǔ)設(shè)施建設(shè)，如水利灌溉、物流交通、信息網(wǎng)絡(luò)、農(nóng)技推廣以及氣象預(yù)警等多方位基礎(chǔ)設(shè)施建設(shè)，降低農(nóng)戶私人種糧成本和種糧風(fēng)險，從源頭上降低耕地非糧化趨勢，才能更好地實現(xiàn)糧食綠色生產(chǎn)發(fā)展。