蔡林軍 文春暉

摘要：糧食安全是“國之大者”，隨著農(nóng)業(yè)勞動力大量轉(zhuǎn)移，糧食生產(chǎn)韌性是否受到影響及受到的影響有多大，已成為新時期的重要議題?；谥袊Z食主產(chǎn)區(qū)2010—2022年省域面板數(shù)據(jù)，通過構建主產(chǎn)區(qū)糧食生產(chǎn)韌性的綜合指標體系，首先使用熵權與變異系數(shù)綜合賦權法對其進行測度，其次選用固定效應模型實證檢驗農(nóng)業(yè)勞動力轉(zhuǎn)移對糧食生產(chǎn)韌性影響的內(nèi)在邏輯。研究發(fā)現(xiàn)：農(nóng)業(yè)勞動力轉(zhuǎn)移有助于主產(chǎn)區(qū)糧食生產(chǎn)韌性的提升，且對糧食生產(chǎn)韌性的促進效應在較長時間序列上均穩(wěn)健成立。具體而言，農(nóng)業(yè)勞動力轉(zhuǎn)移每增加1%，主產(chǎn)區(qū)糧食生產(chǎn)韌性水平將會平均增加1.712 0%。研究還發(fā)現(xiàn)，北方和南方主產(chǎn)區(qū)農(nóng)業(yè)勞動力轉(zhuǎn)移影響糧食生產(chǎn)韌性水平存在明顯異質(zhì)性，北方主產(chǎn)區(qū)的提升作用更為顯著。因此，進一步引導農(nóng)業(yè)勞動力非農(nóng)轉(zhuǎn)移推動區(qū)域糧食生產(chǎn)韌性水平提升，是新時期保障我國糧食安全的有效途徑。

關鍵詞：糧食主產(chǎn)區(qū)；勞動力轉(zhuǎn)移；糧食生產(chǎn)韌性；糧食安全

中圖分類號：F32? 文獻標識碼：A? 文章編號：2095-5553 （2024） 03-0313-09

Study on the influence of labor transfer on the resilience of food production：

Based on the evidence of major grain producing areas in China

Cai Linjun， Wen Chunhui

（School of Economics， Hunan Agricultural University， Changsha， 410128， China）

Abstract：

Food security is “the biggest of the country”， with the large transfer of agricultural labor force， whether the resilience of food production is affected and how much it is affected has become an important issue in the new era. Based on the provincial panel data of Chinas major grain-producing areas from 2010 to 2022， this paper constructs a comprehensive index system of grain production resilience in major grain-producing areas. Firstly， entropy weight and coefficient of variation are used to measure it. Secondly， a fixed-effect model is used to empirically test the internal logic of the impact of agricultural labor transfer on grain production resilience. The results show that the transfer of agricultural labor contributes to the improvement of the resilience of grain production in the main producing areas， and the promoting effect on the resilience of grain production is stable in a long time series. Specifically， for every 1% increase in agricultural labor transfer， the resilience level of grain production in the main production areas will increase by an average of 1.712 0%. The study also found that there was significant heterogeneity between the northern and southern major producing areas， and the northern major producing areas had a more significant effect on the resilience of grain production. Therefore， it is an effective way to ensure Chinas food security in the new era by further guiding the non-agricultural transfer of agricultural labor and promoting the improvement of the resilience level of regional grain production.

Keywords：major grain producing areas; transfer of labor force; toughness of grain production; food security

0 引言

農(nóng)為國之本，糧足則國安。我國已連續(xù)8年糧食產(chǎn)量穩(wěn)定在6.5×1011 kg以上，為保障糧食安全提振了信心；但不可否認的是，農(nóng)業(yè)勞動力轉(zhuǎn)移正成為我國二元經(jīng)濟結構背景下長期存在的社會現(xiàn)象，大量青壯年勞動力的非農(nóng)轉(zhuǎn)移帶來農(nóng)業(yè)勞動力老齡化、婦女化特征趨勢凸顯，加劇了人們對糧食生產(chǎn)韌性的擔憂。此外，糧食作為特殊商品，具有需求價格彈性小、產(chǎn)業(yè)鏈聯(lián)系數(shù)偏低等特征，使得糧食生產(chǎn)環(huán)節(jié)呈現(xiàn)出脆弱性。在此情景下，如何提高主產(chǎn)區(qū)糧食生產(chǎn)韌性，在生產(chǎn)端筑牢糧食安全防御底線，以保障糧食供給穩(wěn)定增長和糧食質(zhì)量不斷改善，既是一個現(xiàn)實問題，也是一個理論問題。

實現(xiàn)糧食產(chǎn)業(yè)高質(zhì)量發(fā)展與農(nóng)業(yè)農(nóng)村現(xiàn)代化離不開糧食生產(chǎn)韌性的鍛造。何為“韌性”，Holling［1］率先在生態(tài)學領域運用這一概念定義生態(tài)系統(tǒng)所具有的抵抗沖擊、恢復原態(tài)的能力。糧食體系韌性應從個體到國家多層次多角度共同發(fā)力［2］；可見，糧食生產(chǎn)韌性是內(nèi)外部有機結合抵御外部沖擊的能力。學界對于如何鍛造糧食生產(chǎn)韌性的研究較少，其涉及更多的是對提升農(nóng)業(yè)韌性的探討。如羅光強等［3］指出農(nóng)業(yè)的綠色化、農(nóng)業(yè)結構的合理化和高級化是提升農(nóng)業(yè)韌性的關鍵因素，增加糧食供給、擴大地區(qū)市場規(guī)模對農(nóng)業(yè)韌性有促進作用［4］。糧食生產(chǎn)韌性作為農(nóng)業(yè)韌性的一個方面，提升農(nóng)業(yè)韌性的因素為糧食生產(chǎn)韌性的鍛造提供了指引。然而，對于勞動力的非農(nóng)轉(zhuǎn)移，秦立建等［5］認為農(nóng)村留守人員減少了糧食生產(chǎn)中的實際勞動投入對糧食生產(chǎn)會有負面影響；侯孟陽等［6］則強調(diào)農(nóng)業(yè)勞動力轉(zhuǎn)移對農(nóng)業(yè)生產(chǎn)效率具有促進作用。

綜上所述，學者們關于勞動力非農(nóng)轉(zhuǎn)移對糧食生產(chǎn)的影響持有不同的觀點，遺憾的是，他們鮮有關注其對糧食生產(chǎn)綜合因素即韌性的考察，而糧食生產(chǎn)韌性關系到糧食安全能否得到根本保障。鑒于此，本文嘗試做出可能的邊際貢獻在于：從我國的國情出發(fā)，在界定糧食生產(chǎn)韌性內(nèi)涵的基礎上，構建具有本國特色的主產(chǎn)區(qū)糧食生產(chǎn)韌性指標評價體系，在對主產(chǎn)區(qū)13個省域的糧食生產(chǎn)韌性綜合發(fā)展指數(shù)進行測算時，選擇使用熵權與變異系數(shù)綜合賦權法；同時利用我國糧食主產(chǎn)區(qū)2010—2022年的面板數(shù)據(jù)進行實證分析，考察勞動力轉(zhuǎn)移對主產(chǎn)區(qū)糧食生產(chǎn)韌性的影響方向和效果及不同省域勞動力轉(zhuǎn)移影響糧食生產(chǎn)韌性的區(qū)域差異，并基于研究結論提煉相應政策建議，提高糧食生產(chǎn)韌性鍛造政策實施的匹配度，為主產(chǎn)區(qū)糧食生產(chǎn)韌性穩(wěn)步提升提供理論與現(xiàn)實指導。

1 理論分析與研究假說

從本質(zhì)上看，主產(chǎn)區(qū)糧食安全是在各種生產(chǎn)要素的綜合作用下，糧食供給穩(wěn)定增長及糧食質(zhì)量持續(xù)改善［7］。鍛造糧食生產(chǎn)韌性是糧食生產(chǎn)環(huán)節(jié)必備的要求，本文所討論的糧食生產(chǎn)韌性是指糧食生產(chǎn)過程中生產(chǎn)要素受到外界沖擊時，要素分配組合能迅速調(diào)整及向新均衡路徑轉(zhuǎn)變以實現(xiàn)應對沖擊的能力；其主要是抵抗能力、復原能力兩種能力的結合。農(nóng)業(yè)生產(chǎn)中存在三大核心要素，即勞動力、土地、資本［8］，那么農(nóng)業(yè)勞動力轉(zhuǎn)移對主產(chǎn)區(qū)糧食生產(chǎn)韌性的影響方向及效果如何？機制如何？正是本文所考察的核心內(nèi)容。農(nóng)業(yè)勞動力轉(zhuǎn)移可通過以下途徑作用于主產(chǎn)區(qū)糧食生產(chǎn)韌性。

1.1 農(nóng)業(yè)勞動力轉(zhuǎn)移促進農(nóng)戶家庭收入增加

現(xiàn)有研究認為勞動力非農(nóng)轉(zhuǎn)移增加了農(nóng)戶工資性收入，進而對農(nóng)戶總收入產(chǎn)生直接影響和促進作用［9， 10］。農(nóng)戶收入增加，一方面，推動糧食生產(chǎn)過程資本投入量增加，提高了農(nóng)戶對農(nóng)用機械、化肥等生產(chǎn)要素和優(yōu)質(zhì)糧種的購買能力；另一方面，能促進農(nóng)業(yè)生產(chǎn)區(qū)域交通體系的完善及農(nóng)業(yè)生產(chǎn)效率的提升［11， 12］。與此同時，糧食生產(chǎn)作為資源密集型產(chǎn)業(yè)，分工越來越細化迂回生產(chǎn)鏈不斷拉長，糧食生產(chǎn)發(fā)展面臨新挑戰(zhàn)、新業(yè)態(tài)、新形式，對資本將有更大需求，正因如此非農(nóng)轉(zhuǎn)移獲得的工資性收入能緩解其變化下的流動性約束，推進糧食生產(chǎn)過程各層次主體間的上下聯(lián)動，最終實現(xiàn)糧食生產(chǎn)韌性能力的提升。

1.2 農(nóng)業(yè)勞動力轉(zhuǎn)移促進農(nóng)地規(guī)?；?jīng)營

農(nóng)業(yè)勞動力轉(zhuǎn)移為實現(xiàn)耕地規(guī)?；?jīng)營奠定了良好基礎，給糧食生產(chǎn)帶來了規(guī)模經(jīng)濟。農(nóng)地經(jīng)營規(guī)模越大，種植結構“趨糧化”效應越顯著［13］；具體而言，在農(nóng)業(yè)機械化水平上升地推動下，農(nóng)地規(guī)?；?jīng)營不僅能讓土地的利用率得到進一步提高，而且可實現(xiàn)糧食生產(chǎn)由粗放型向集約型方式的高質(zhì)量生產(chǎn)轉(zhuǎn)變，即勞動力轉(zhuǎn)移通過直接作用于主產(chǎn)區(qū)農(nóng)地流轉(zhuǎn)，推進農(nóng)地規(guī)?；?jīng)營為開展機械化作業(yè)創(chuàng)造了條件，從而推動區(qū)域糧食生產(chǎn)韌性鍛造；同時大型機械的“深松翻”和“少免耕”技術能夠改善土壤結構，能起到糧食增產(chǎn)作用［14］。此外，我國轉(zhuǎn)移的農(nóng)業(yè)勞動力中，存有一部分季節(jié)性勞動力轉(zhuǎn)移［15］，因此會呈現(xiàn)出一定的周期性，這部分勞動力在外地賺取非農(nóng)收入后，又會在農(nóng)忙時節(jié)流回農(nóng)村，不會給糧食生產(chǎn)造成顯著不利影響。

1.3 農(nóng)業(yè)勞動力轉(zhuǎn)移加大地區(qū)糧食生產(chǎn)韌性水平差異

我國糧食生產(chǎn)具有土地空間廣、農(nóng)業(yè)企業(yè)多、農(nóng)戶規(guī)模大等特征，屬于典型的大農(nóng)業(yè)生產(chǎn)。主產(chǎn)區(qū)省域間自然稟賦、經(jīng)濟基礎、技術支持及政策安排上存在較大差異，13個糧食主產(chǎn)區(qū)中，北方糧食主產(chǎn)區(qū)有7個，南方糧食主產(chǎn)區(qū)有6個，相較于北方主產(chǎn)區(qū)，南方主產(chǎn)區(qū)多以丘陵山區(qū)為主，耕地“細碎化”不利于農(nóng)地規(guī)?；?jīng)營和機械化作業(yè)，由此造成農(nóng)業(yè)用地撂荒問題凸顯和種植結構“非糧化”現(xiàn)象；而種植結構調(diào)整程度、勞動力機械替代難易度等因素均會一定程度上制約農(nóng)業(yè)勞動力轉(zhuǎn)移對糧食生產(chǎn)韌性的提升作用［16， 17］。

綜上，本文提出以下兩個需要進一步借助實證來檢驗的研究假說。

假說1：農(nóng)業(yè)勞動力轉(zhuǎn)移有助于主產(chǎn)區(qū)糧食生產(chǎn)韌性的提升。

假說2：農(nóng)業(yè)勞動力轉(zhuǎn)移促進主產(chǎn)區(qū)糧食生產(chǎn)韌性提升作用呈現(xiàn)區(qū)域異質(zhì)性。

2 模型設定、變量說明與數(shù)據(jù)來源

2.1 模型設定

為檢驗農(nóng)業(yè)勞動力轉(zhuǎn)移對主產(chǎn)區(qū)糧食生產(chǎn)韌性的影響，根據(jù)機理分析和假設，本文構建以下基準回歸模型

ProduceI，T=α0＋α1TRI，T＋α2ControlsI，T＋μI＋φI＋εI，T（1）

式中：ProduceI，T——熵權與變異系數(shù)綜合賦權法測算出的主產(chǎn)區(qū)糧食生產(chǎn)韌性水平，具體包括生產(chǎn)硬性約束、經(jīng)濟基礎條件、經(jīng)濟增長狀況、生態(tài)治理能力、技術創(chuàng)新進步；TRI，T——主產(chǎn)區(qū)農(nóng)業(yè)勞動力轉(zhuǎn)移情況；ControlsI，T——影響主產(chǎn)區(qū)糧食生產(chǎn)韌性的控制變量組，主要包括主產(chǎn)區(qū)城鎮(zhèn)化水平、經(jīng)濟發(fā)展水平、非農(nóng)經(jīng)濟占比、農(nóng)產(chǎn)品生產(chǎn)價格指數(shù)、財政支農(nóng)水平、農(nóng)民專業(yè)合作社數(shù)量及農(nóng)業(yè)經(jīng)營收入在純收入中所占的比值；α0、α1、α2——待估計的參數(shù)；μI——省份I固定效應；φT——年份T固定效應；εI，T——隨機擾動項。

2.2 變量描述和數(shù)據(jù)說明

被解釋變量是糧食生產(chǎn)韌性水平。本文在借鑒相關韌性研究的基礎上，嘗試對糧食生產(chǎn)韌性構建出多層次的綜合性評價指標體系。參考郝愛民等［18］研究，將主產(chǎn)區(qū)糧食生產(chǎn)韌性劃分為抵抗力和復原力2個一級指標，將抵抗力劃分為生產(chǎn)硬性約束、經(jīng)濟基礎條件，將復原力劃分為經(jīng)濟增長狀況、生態(tài)治理能力、技術創(chuàng)新進步的5個二級指標，建構3個層次16個三級指標的評價體系，測度主產(chǎn)區(qū)糧食生產(chǎn)韌性水平見表1。

層次分析法、熵權法、主成分分析法等均是測度糧食生產(chǎn)韌性最為常用的方法。為了有效規(guī)避主觀因素給模型造成的影響，同時，為了消除單一客觀賦權的誤差與偏頗，提高測度主產(chǎn)區(qū)糧食生產(chǎn)韌性水平的可靠性；參考谷城等［19］測算產(chǎn)業(yè)鏈韌性的做法，本文采用熵權與變異系數(shù)綜合賦權法測度主產(chǎn)區(qū)糧食生產(chǎn)韌性水平。具體測算步驟如下。

第一，構建評價系統(tǒng)初始數(shù)據(jù)矩陣xij，對指標進行標準化處理。

xij=xij-minxijmaxxij-minxij?xij為正向指標maxxij-xijmaxxij-minxij?xij為負向指標（2）

為避免數(shù)值0對取對數(shù)的影響，本文對標準化后的矩陣進行了平移轉(zhuǎn)換，每個元素均加0.001，將其矩陣記為x0ij。

第二，熵權法測度指標權重。

Pij=x0ij/∑mi=1x0ij（3）

ej=-K∑mi=1PijlnPij（4）

Wj=（1-ej）/（∑nj=1（1-ej）（5）

式中：Pij——j指標下i樣本所占比重；ej——j指標的信息熵；K——調(diào)節(jié)系數(shù)，K=1/lnm；m——樣本個數(shù)；Wj——利用熵權法計算得到的j指標權重。第三，變異系數(shù)法測算指標權重。

Dj=1M∑mi=1（xij-x-j）2（6）

Sj=Dj（7）

CVj=Sj/x-j（8）

Zj=CVj/∑nj=1CVj（9）

式中：Dj——方差；Sj——標準差；x-j——均值；CVj——j指標的變異系數(shù)；Zj——利用變異系數(shù)法得到的j指標權重。

第四，根據(jù)熵權法和變異系數(shù)法得到的兩種指標權重，測度主產(chǎn)區(qū)糧食生產(chǎn)韌性水平。計算公式如式（10）～式（12）所示。

f（ξi）=∑nj=1Wjx0ij（10）

f（τi）=∑nj=1Zjx0ij（11）

Ri=f（ξi）×f（τi）（12）

式中：f（ξi）——熵權法計算出的指標權重；f（τi）——變異系數(shù)法計算出的指標權重；Ri——幾何平均值，基于熵權與變異系數(shù)綜合賦權法測算的主產(chǎn)區(qū)糧食生產(chǎn)韌性水平。

根據(jù)《中國統(tǒng)計年鑒》《中國農(nóng)業(yè)機械統(tǒng)計年鑒》2010—2022年數(shù)據(jù)，利用熵權與變異系數(shù)綜合賦權法，嘗試基于時間序列來測算主產(chǎn)區(qū)糧食生產(chǎn)韌性的動態(tài)發(fā)展趨勢，見表2。

根據(jù)表2數(shù)據(jù)顯示，從時間維度看，在2010—2022年這一時間區(qū)間中，主產(chǎn)區(qū)糧食生產(chǎn)韌性發(fā)展平均水平處于一個穩(wěn)定的上升狀態(tài)，在2010—2016年期間其韌性水平發(fā)展指數(shù)均值只有0.337 0，發(fā)展到2017—2022年間便提升至了0.417 8。站在空間這一維度看，其發(fā)展水平則呈現(xiàn)出了顯著的差異性，其中2010—2022年間糧食生產(chǎn)韌性的平均發(fā)展水平達到了0.375 3，一共有6個省份高于平均值，其中江蘇省最高，韌性水平發(fā)展指數(shù)達到了0.502 1。總體表明，在積極轉(zhuǎn)型構建現(xiàn)代糧食產(chǎn)業(yè)體系的新時代背景下，主產(chǎn)區(qū)作為我國大農(nóng)業(yè)地區(qū)，糧食生產(chǎn)韌性水平穩(wěn)步上升，四川省等依托自然資源稟賦優(yōu)勢，在多年的糧食生產(chǎn)實踐中不斷動態(tài)調(diào)整生產(chǎn)結構模式，逐步向高產(chǎn)量、高產(chǎn)值的優(yōu)勢糧食作物方向轉(zhuǎn)變，有著最快的增長速度；江蘇等省份借助于自身地理優(yōu)勢，集中力量開展農(nóng)業(yè)技術革新，促使智慧農(nóng)業(yè)發(fā)展和科技投入水平增量齊頭并進，推動糧食生產(chǎn)韌性水平提升，其發(fā)展水平位居前列；其他地區(qū)受耕地地理屬性和經(jīng)濟底蘊的雙重約束，其糧食生產(chǎn)韌性發(fā)展的區(qū)域表征也正處于穩(wěn)健上升狀態(tài)。

將農(nóng)業(yè)勞動力轉(zhuǎn)移設定為核心解釋變量。本文借鑒張紅麗等［20］選取農(nóng)業(yè)勞動力轉(zhuǎn)移比作為農(nóng)業(yè)勞動力轉(zhuǎn)移的衡量指標。其測算方法是鄉(xiāng)村從業(yè)人員數(shù)量和農(nóng)林牧漁業(yè)從業(yè)人員數(shù)量的差與鄉(xiāng)村從業(yè)人員數(shù)量之比來表示；此外，另有學者在研究過程中，將農(nóng)業(yè)勞動力轉(zhuǎn)移的衡量指標設定為鄉(xiāng)村、農(nóng)林牧漁業(yè)兩類從業(yè)人員之間存在的數(shù)量差值，本文在穩(wěn)健性檢驗中用此方法作為替代。

除農(nóng)業(yè)勞動力轉(zhuǎn)移之外，主產(chǎn)區(qū)糧食生產(chǎn)韌性還受到諸如經(jīng)濟發(fā)展、產(chǎn)業(yè)結構、城鎮(zhèn)化水平、農(nóng)產(chǎn)品生產(chǎn)價格、財政支農(nóng)水平、農(nóng)民專業(yè)合作及農(nóng)業(yè)經(jīng)營性收入等因素影響，基于此，選取人均實際GDP以控制經(jīng)濟發(fā)展水平；選取非農(nóng)經(jīng)濟占比，來控制地區(qū)產(chǎn)業(yè)結構；選取城鎮(zhèn)化率控制城鎮(zhèn)化水平；選取農(nóng)產(chǎn)品生產(chǎn)價格指數(shù)控制農(nóng)產(chǎn)品生產(chǎn)價格；此外還選取了財政支農(nóng)、農(nóng)民專業(yè)合作社數(shù)量及農(nóng)業(yè)經(jīng)營性收入在純收入的比重來控制財政支農(nóng)水平、農(nóng)民專業(yè)合作水平和農(nóng)業(yè)經(jīng)營性收入，見表3。

2.3 數(shù)據(jù)來源與描述性統(tǒng)計

本文以我國13個糧食主產(chǎn)區(qū)2010—2022年的面板數(shù)據(jù)作為研究樣本區(qū)間，以2010年為基期對相關經(jīng)濟指標進行了消脹處理，《中國統(tǒng)計年鑒》及各省份的統(tǒng)計年鑒是本文研究變量數(shù)據(jù)的主要來源，少量的數(shù)據(jù)缺失使用線性插值法將其補齊。通過觀察表4中對糧食主產(chǎn)區(qū)各個變量開展的描述性統(tǒng)計，可以發(fā)現(xiàn)糧食生產(chǎn)韌性發(fā)展水平各個糧食主產(chǎn)區(qū)之間存在明顯差異，因其主產(chǎn)區(qū)糧食生產(chǎn)韌性發(fā)展水平對數(shù)值的均值為-1.002 0，而標準誤為0.301 9。

3 實證分析

3.1 基準回歸結果分析

為有效考察農(nóng)業(yè)勞動力轉(zhuǎn)移對主產(chǎn)區(qū)糧食生產(chǎn)韌性所產(chǎn)生的具體影響方向及效果，本文嘗試借助于2010—2022年期間我國13個糧食主產(chǎn)區(qū)所產(chǎn)生的全部數(shù)據(jù)樣本，將被解釋變量設定為主產(chǎn)區(qū)糧食生產(chǎn)韌性水平，核心解釋變量設定為農(nóng)業(yè)勞動力轉(zhuǎn)移開展面板回歸估計。在參數(shù)估計前，為了檢驗各個變量數(shù)據(jù)序列的平穩(wěn)性，首先選擇使用LLC檢測法開展相關檢驗，最終得到的檢驗結果顯示各變量數(shù)據(jù)平穩(wěn)；隨后為了進一步檢驗變量數(shù)據(jù)是否存在協(xié)整關系，使用pedroni面板協(xié)整檢驗方法，得到的檢驗結果表示變量之間在1%的顯著性水平下通過了統(tǒng)計檢驗表明存在協(xié)整關系，可以開展面板回歸。此外，經(jīng)過Hausman檢驗表明，對參數(shù)開展估計時使用固定效應模型更具合理性。

針對“農(nóng)業(yè)勞動力轉(zhuǎn)移—糧食生產(chǎn)韌性”的基準關系進行了實證檢驗，其中表5第（1）列控制“年份—省份”固定效應。根據(jù)回歸結果可知，TR對應的估計系數(shù)為正值，并且保持在1%的水平下顯著。將控制變量納入之后，表5第（2）列僅控制了省份，第（3）列對“年份—省份”進行雙向固定，結論依舊穩(wěn)健。根據(jù)前文分析，可能的解釋為：一方面，在農(nóng)業(yè)勞動力轉(zhuǎn)移的推動下，農(nóng)村居民家庭的非農(nóng)業(yè)收入會顯著增加，農(nóng)民收入水平進一步提高，繼而促使農(nóng)業(yè)生產(chǎn)性投入增多，包括農(nóng)業(yè)基礎設施建設、農(nóng)用機械化、化肥施用占比等；另一方面，從農(nóng)村中轉(zhuǎn)移出去的勞動力依舊掌握著土地轉(zhuǎn)用權，能夠為耕地規(guī)?；?jīng)營奠定基礎，給糧食生產(chǎn)帶來了規(guī)模經(jīng)濟；因此，農(nóng)業(yè)勞動力轉(zhuǎn)移的這些影響推進糧食生產(chǎn)過程各層次主體間的上下聯(lián)動，最終提升糧食生產(chǎn)過程總體抵御沖擊、適應沖擊的韌性能力。

此外，從各控制變量來看，地區(qū)經(jīng)濟發(fā)展水平GDP、農(nóng)業(yè)經(jīng)營收入在總收入中的占比GS的估計系數(shù)顯著為正，對應的影響系數(shù)為1.738 1、1.184 1，表明農(nóng)戶家庭收入增加確實有利于諸如農(nóng)用設備、固定資產(chǎn)投資和水土流水治理等基礎性工程建設，以防御眾多負外部性造成糧食生產(chǎn)過程中的不良影響實現(xiàn)其韌性的提升。城鎮(zhèn)化水平UB的估計系數(shù)顯著為負，對應的影響系數(shù)為-3.826 0，可能的原因是，糧食主產(chǎn)區(qū)部分省份限于自然地域約束，機械化水平難以施展，城鎮(zhèn)化水平上升，非農(nóng)人口增加，導致農(nóng)業(yè)出現(xiàn)大面積撂荒，致使糧食播種面積、復種指數(shù)等下降，從而使糧食生產(chǎn)韌性水平降低。非農(nóng)經(jīng)濟占比IC、財政支農(nóng)水平FN估計系數(shù)為正符合經(jīng)濟預期，農(nóng)產(chǎn)品生產(chǎn)價格水平FPI、農(nóng)民專業(yè)合作社數(shù)量NH估計系數(shù)為負，但均不顯著。

3.2 穩(wěn)健性檢驗與內(nèi)生性緩解

3.2.1 替換模型

因本文測度的被解釋變量（主產(chǎn)區(qū)糧食生產(chǎn)韌性的對數(shù)）有著取值上限，具體上限為0，符合受限因變量模型的相關條件，在重估模型（1）時可選擇使用受限變量Tobit模型，并且繼續(xù)選擇使用固定效應。Tobit模型對應的估計結果對應為表6第（1）列。將Tobit、基準回歸兩個模型對應的估計結果進行對比分析，可知農(nóng)業(yè)勞動轉(zhuǎn)移對糧食生產(chǎn)韌性的彈性系數(shù)影響方向一致。表明在更換模型后所得到的結論依舊保持穩(wěn)健。

3.2.2 變換被解釋變量測算方法

為避免因指標測算方法不同對回歸結果帶來影響，利用主成分分析法重新測算主產(chǎn)區(qū)糧食生產(chǎn)韌性水平，同樣使用“年份—省份”雙向固定效應模型進行估計，結果如表6第（2）列所示，農(nóng)業(yè)勞動力轉(zhuǎn)移所對應的彈性系數(shù)為正值，并且在1%的水平下保持顯著。

3.2.3 變換核心解釋變量測算方法

參考鄭晶等［21］對農(nóng)業(yè)勞動力轉(zhuǎn)移的測度方法，用“鄉(xiāng)村從業(yè)人員數(shù)量-農(nóng)林牧漁業(yè)從業(yè)人員數(shù)量”測度，結果如表6第（3）列所示，農(nóng)業(yè)勞動力轉(zhuǎn)移的系數(shù)在1%的顯著性水平下為正，驗證了基準回歸結果的穩(wěn)健性。

3.2.4 內(nèi)生性緩解與工具變量

雖然本文最大程度控制了一部分關鍵性變量，并為了有效緩解因變量遺漏而出現(xiàn)的內(nèi)生性偏誤，選擇使用固定效應模型，但依舊無法徹底消除其帶來的影響，特別是當主產(chǎn)區(qū)糧食生產(chǎn)韌性、農(nóng)業(yè)勞動力轉(zhuǎn)移兩者之間存在互為因果的關系時，可能出現(xiàn)下述問題：農(nóng)業(yè)勞動力轉(zhuǎn)移受到糧食生產(chǎn)韌性帶來的反向影響，如前所述，經(jīng)濟增長會提高糧食生產(chǎn)韌性，可能使該區(qū)域機械化水平上升，誘致農(nóng)業(yè)勞動力剩余，進而加快農(nóng)業(yè)勞動力轉(zhuǎn)移進程。為此，本文嘗試利用工具變量法解決該內(nèi)生性問題，參考黃楓等［22］的做法，使用樣本所在省域除自身之外其他樣本的糧食生產(chǎn)韌性均值（Iv1）和主產(chǎn)區(qū)省域公路里程的對數(shù)（Iv2）作為當期農(nóng)業(yè)勞動力轉(zhuǎn)移的工具變量。此外，為了有效消除極端值影響，嘗試縮尾處理變量2%～98%之外的數(shù)據(jù)；具體估計結果為表6第（4）列。其中，過度識別得到的P值為0.887 1，表明工具變量能夠有效滿足外生性檢驗，且根據(jù)RKF檢驗結果可知弱工具變量這一問題并不存在，因此本文選擇使用的工具變量有著較高的合理性。在工具變量法有效的這一情況下，基準回歸得到的結論依舊保持成立。

至此，本文認為假說1得以實證驗證，即農(nóng)業(yè)勞動力轉(zhuǎn)移有助于主產(chǎn)區(qū)糧食生產(chǎn)韌性提升。

3.3 農(nóng)業(yè)勞動力轉(zhuǎn)移對主產(chǎn)區(qū)糧食生產(chǎn)韌性影響的異質(zhì)性分析

3.3.1 南北方主產(chǎn)區(qū)異質(zhì)性

我國糧食生產(chǎn)環(huán)境存在明顯的地域差異。根據(jù)主產(chǎn)區(qū)各省份自然資源稟賦不同，將糧食主產(chǎn)區(qū)分為北方糧食主產(chǎn)區(qū)和南方糧食主產(chǎn)區(qū)，以進一步考察區(qū)分北方主產(chǎn)區(qū)與南方主產(chǎn)區(qū)后農(nóng)業(yè)勞動力轉(zhuǎn)移對糧食生產(chǎn)韌性的影響異質(zhì)性南方糧食主產(chǎn)區(qū)一共包括6個省份，即湖北、安徽、湖南、江蘇、四川、江西；北方糧食主產(chǎn)區(qū)一共包括7個省份，即河北、河南、吉林、黑龍江、遼寧、山東、內(nèi)蒙古。檢驗結果如表7第（1）、第（2）列所示。北方糧食主產(chǎn)區(qū)農(nóng)業(yè)勞動力轉(zhuǎn)移對糧食生產(chǎn)韌性水平發(fā)展的影響系數(shù)為1.917 8，并在1%的統(tǒng)計水平上顯著；南方糧食主產(chǎn)區(qū)農(nóng)業(yè)勞動力轉(zhuǎn)移對糧食生產(chǎn)韌性水平發(fā)展的影響系數(shù)為0.009 0，但不顯著；可能的原因是，北方多以平原地區(qū)為主，農(nóng)業(yè)勞動力轉(zhuǎn)移規(guī)模越大，農(nóng)地規(guī)模化經(jīng)營帶來的糧食生產(chǎn)規(guī)模經(jīng)濟越顯著，農(nóng)戶越會增加勞動要素的機械投入，同時，轉(zhuǎn)移的勞動力也收獲了工資性收入，最終實現(xiàn)糧食生產(chǎn)抵抗力和復原力提升即韌性的鍛造；南方多以丘陵地區(qū)為主，農(nóng)業(yè)勞動力轉(zhuǎn)移，由于土地“細碎化”，機械化作業(yè)施展成本較大，部分地區(qū)出現(xiàn)農(nóng)業(yè)拋荒現(xiàn)象，削弱糧食生產(chǎn)抵抗力和復原力。

至此，本文認為假設2得以實證驗證，即農(nóng)業(yè)勞動力轉(zhuǎn)移促進主產(chǎn)區(qū)糧食生產(chǎn)韌性提升作用呈現(xiàn)區(qū)域異質(zhì)性。

3.3.2 糧食功能區(qū)異質(zhì)性

截至2022年，我國糧食凈調(diào)出省僅剩5個，分別是黑龍江、河南、吉林、內(nèi)蒙古和安徽，為探究農(nóng)業(yè)勞動力轉(zhuǎn)移對主產(chǎn)區(qū)糧食生產(chǎn)韌性是否在糧食功能區(qū)存在差異，本文將13個糧食主產(chǎn)區(qū)按照糧食凈調(diào)出和非凈調(diào)出省進行劃分，其中非糧食凈調(diào)出省是剩余的8個省份。檢驗結果如表7第（3）、第（4）列所示，可見，不論是糧食凈調(diào)出省還是非糧食凈調(diào)出省農(nóng)業(yè)勞動力轉(zhuǎn)移均對糧食生產(chǎn)韌性有正向影響，但糧食凈調(diào)出省的影響更為顯著。值得注意的是，糧食凈調(diào)出省份中有4個來自北方糧食主產(chǎn)區(qū)，這再次驗證了假說2成立和結果的穩(wěn)健性。

4 結論與政策建議

4.1 結論

本文試圖對農(nóng)業(yè)勞動力轉(zhuǎn)移影響主產(chǎn)區(qū)糧食生產(chǎn)韌性的理論機制予以闡述，并借助于現(xiàn)實數(shù)據(jù)開展實證研究。最終的理論機制表明，糧食生產(chǎn)韌性從整體角度展示了，當生產(chǎn)要素遭受到外界的沖擊時，要素分配組合能迅速調(diào)整及向新均衡路徑轉(zhuǎn)變以實現(xiàn)應對沖擊和自我恢復的能力。農(nóng)業(yè)勞動力轉(zhuǎn)移通過影響農(nóng)戶收入和農(nóng)地流轉(zhuǎn)來鍛造糧食生產(chǎn)韌性。實證檢驗中，本文在針對糧食生產(chǎn)韌性構建綜合評價指標體系時嘗試從復原能力、抵抗能力等層面著手構建，使用熵權與變異系數(shù)綜合賦權法對其進行測度，并利用13個主產(chǎn)區(qū)2010—2022年的相關數(shù)據(jù)考察農(nóng)業(yè)勞動力轉(zhuǎn)移對主產(chǎn)區(qū)糧食生產(chǎn)韌性的影響效應及差異。結果表明：（1）主產(chǎn)區(qū)省份通過充分挖掘地域資源稟賦優(yōu)勢，在時序維度上糧食生產(chǎn)韌性水平呈現(xiàn)顯著上升趨勢。具體而言，2010—2022年13個主產(chǎn)區(qū)糧食生產(chǎn)韌性平均水平從2010—2016年間的0.337 0提升至2017—2022年間的0.417 8。（2）農(nóng)業(yè)勞動力轉(zhuǎn)移對主產(chǎn)區(qū)糧食生產(chǎn)韌性呈現(xiàn)出了顯著的正向影響，在綜合考慮內(nèi)生性問題并進行穩(wěn)健性檢驗后，結論依舊保持成立。具體而言，農(nóng)業(yè)勞動力轉(zhuǎn)移每增加1%，主產(chǎn)區(qū)糧食生產(chǎn)韌性水平將會平均增加1.712 0%。（3）將主產(chǎn)區(qū)進行區(qū)域劃分發(fā)現(xiàn)北方主產(chǎn)區(qū)和南方主產(chǎn)區(qū)因自然資源稟賦差異，農(nóng)業(yè)勞動力轉(zhuǎn)移對糧食生產(chǎn)韌性水平存在明顯異質(zhì)性。

4.2 政策建議

1） 政府應積極引導農(nóng)業(yè)剩余勞動力向非農(nóng)部門轉(zhuǎn)移，促進主產(chǎn)區(qū)糧食生產(chǎn)韌性高水平發(fā)展。主產(chǎn)區(qū)農(nóng)業(yè)要素的機械投入富有彈性，農(nóng)業(yè)剩余勞動力仍有較大存量，要加快帶動勞動力、土地、資本等生產(chǎn)要素市場化改革及戶籍制度安排，促進城鄉(xiāng)資源要素循環(huán)互通雙向流動。繼續(xù)促進農(nóng)業(yè)剩余人口轉(zhuǎn)移能緩解非農(nóng)勞動力供給約束、增加農(nóng)戶工資性收入及改善人地矛盾，推動農(nóng)地規(guī)?；?jīng)營和機械化開展。

2） 加大實施農(nóng)地制度改革和鼓勵農(nóng)地流轉(zhuǎn)，引導土地規(guī)?；?jīng)營，提高農(nóng)用機械化水平。農(nóng)地流轉(zhuǎn)是促進土地規(guī)?；?jīng)營的必要條件，而農(nóng)地規(guī)模化經(jīng)營機械化作業(yè)才能得以施展；此外，還需不斷完善農(nóng)地承包權和經(jīng)營權的長效保護機制，既能解決農(nóng)業(yè)勞動力轉(zhuǎn)移的后顧之憂，也能激勵農(nóng)業(yè)大戶長期投資生產(chǎn)的積極性，達到提高糧食生產(chǎn)韌性的效應。

3） 政府需提升適應于丘陵地區(qū)作業(yè)的農(nóng)機研制投入，破解丘陵地區(qū)資源約束，釋放糧食生產(chǎn)要素活力。南方糧食主產(chǎn)區(qū)地形大多以丘陵為主，不利于機械化作業(yè)，由此致使部分地區(qū)出現(xiàn)農(nóng)地撂荒現(xiàn)象，對糧食生產(chǎn)韌性產(chǎn)生顯著削弱作用，政府應加大農(nóng)機研制投入和農(nóng)業(yè)基礎設施建設，優(yōu)化財政支農(nóng)政策，緩解主產(chǎn)區(qū)資源環(huán)境約束，保持糧食生產(chǎn)長期穩(wěn)定增長促進其韌性的提升。

4） 積極宣傳普及節(jié)肥技術，走綠色農(nóng)業(yè)發(fā)展道路。農(nóng)用化肥投入在糧食生產(chǎn)中有重要作用，但是不合理的化肥施用量會導致土壤質(zhì)量退化、糧食不安全和給生態(tài)環(huán)境造成污染，故需大力推廣節(jié)肥技術、優(yōu)化化肥用量配比及新型有機化肥研發(fā)，走綠色發(fā)展道路。

參 考 文 獻

［1］Holling C S. Resilience and stability of ecological systems ［J］. Annual Review of Ecology and Systematics， 1973， 4（1）： 1-23.

［2］Ji Y， Yu X， Zhong F. Machinery investment decision and off-farm employment in rural China ［J］. China Economic Review， 2012， 23（1）： 71-80.

［3］羅光強， 王煥. 數(shù)字普惠金融對中國糧食主產(chǎn)區(qū)農(nóng)業(yè)高質(zhì)量發(fā)展的影響［J］. 經(jīng)濟縱橫， 2022（7）： 107-117.Luo Guangqiang， Wang Huan. The impact of digital inclusive finance on the high quality development of agriculture in Chinas main grain producing areas ［J］. Economic Review Journal， 2022（7）： 107-117.

［4］董亞妮， 祁春節(jié)， 譚明交， 等. 糧食主產(chǎn)區(qū)農(nóng)業(yè)韌性的空間差異及影響因素［J］. 中國農(nóng)業(yè)資源與區(qū)劃， 2023（11）： 1-15.Dong Yani， Qi Chunjie， Tan Mingjiao， et al. Spatial differences and influencing factors of agricultural resilience in major grain producing areas ［J］. Chinese Journal of Agricultural Resources and Regional Planning， 2023（11）： 1-15.

［5］秦立建， 張妮妮， 蔣中一. 土地細碎化、勞動力轉(zhuǎn)移與中國農(nóng)戶糧食生產(chǎn)——基于安徽省的調(diào)查［J］. 農(nóng)業(yè)技術經(jīng)濟， 2011（11）： 16-23.

［6］侯孟陽， 鄧元杰， 姚順波. 農(nóng)村勞動力轉(zhuǎn)移、化肥施用強度與農(nóng)業(yè)生態(tài)效率： 交互影響與空間溢出［J］. 農(nóng)業(yè)技術經(jīng)濟， 2021（10）： 79-94.Hou Mengyang， Deng Yuanjie， Yao Shunbo. Rural labor transfer， fertilizer use intensity and agro-ecological efficiency： Interaction effects and spatial spillover ［J］. Journal of Agrotechnical Economics， 2021（10）： 79-94.

［7］羅光強. 我國糧食主產(chǎn)區(qū)糧食安全責任實現(xiàn)的路徑及對策［J］. 經(jīng)濟縱橫， 2012（1）： 84-87.

［8］姚成勝， 肖雅雯， 楊一單. 農(nóng)業(yè)勞動力轉(zhuǎn)移與農(nóng)業(yè)機械化對中國糧食生產(chǎn)的關聯(lián)影響分析［J］. 農(nóng)業(yè)現(xiàn)代化研究， 2022（2）： 306-317.Yao Chengsheng， Xiao Yawen， Yang Yidan. The integrated impacts of rural labor transfer and agricultural mechanization on Chinas grain production ［J］. Research of Agricultural Modernization， 2022（2）： 306-317.

［9］李谷成， 李燁陽， 周曉時. 農(nóng)業(yè)機械化、勞動力轉(zhuǎn)移與農(nóng)民收入增長——孰因孰果？［J］. 中國農(nóng)村經(jīng)濟， 2018（11）： 112-127.Li Gucheng， Li Yeyang， Zhou Xiaoshi. Agricultural mechanization， labor transfer and the growth of farmers income： A re-examination of causality ［J］. Chinese Rural Economy， 2018（11）： 112-127.

［10］曾龍， 付振奇. 中國城鄉(xiāng)收入差距對農(nóng)業(yè)生產(chǎn)率的影響——基于農(nóng)村勞動力轉(zhuǎn)移與農(nóng)村居民收入的雙重視角［J］. 江漢論壇， 2021（11）： 15-24.

［11］王子成. 農(nóng)村勞動力外出降低了農(nóng)業(yè)效率嗎？［J］. 統(tǒng)計研究， 2015（3）： 54-61.Wang Zicheng. Does migration make rural household less productive？ ［J］. Statistical Research， 2015（3）： 54-61.

［12］馬曉河， 楊祥雪. 城鄉(xiāng)二元結構轉(zhuǎn)換過程中的農(nóng)業(yè)勞動力轉(zhuǎn)移——基于劉易斯第二轉(zhuǎn)折點的驗證［J］. 農(nóng)業(yè)經(jīng)濟問題， 2023（1）： 4-17.Ma Xiaohe， Yang Xiangxue. Agricultural labor migration in the process of urban-rural dual structure transformation： A validation based on Lewis?second turning point ［J］. Issues in Agricultural Econom， 2023（1）： 4-17.

［13］趙思誠， 許慶， 劉進. 勞動力轉(zhuǎn)移、資本深化與農(nóng)地流轉(zhuǎn)［J］. 農(nóng)業(yè)技術經(jīng)濟， 2020（3）： 4-19.Zhao Sicheng， Xu Qing， Liu Jin. Labor migration， capital deepening， and agricultural land circulation in China ［J］. Journal of Agrotechnical Economics， 2020（3）： 4-19.

［14］侯方安. 農(nóng)業(yè)機械化推進機制的影響因素分析及政策啟示——兼論耕地細碎化經(jīng)營方式對農(nóng)業(yè)機械化的影響［J］. 中國農(nóng)村觀察， 2008（5）： 42-48.

［15］程名望， 黃甜甜， 劉雅娟. 農(nóng)村勞動力外流對糧食生產(chǎn)的影響： 來自中國的證據(jù)［J］. 中國農(nóng)村觀察， 2015（6）： 15-21.Cheng Mingwang， Huang Tiantian， Liu Yajuan. The influence of rural labor outflow on food production： From the evidence of China ［J］. China Rural Survey， 2015（6）： 15-21.

［16］鐘甫寧， 陸五一， 徐志剛. 農(nóng)村勞動力外出務工不利于糧食生產(chǎn)嗎？——對農(nóng)戶要素替代與種植結構調(diào)整行為及約束條件的解析［J］. 中國農(nóng)村經(jīng)濟， 2016（7）： 36-47.

［17］鄭旭媛， 徐志剛. 資源稟賦約束、要素替代與誘致性技術變遷——以中國糧食生產(chǎn)的機械化為例［J］. 經(jīng)濟學（季刊）， 2017， 16（1）： 45-66.Zheng Xuyuan， Xu Zhigang. Endowment restriction， factor substitution and induced technological innovation： A case research on the grain producing mechanization in China ［J］. China Economic Quarterly， 2017， 16（1）： 45-66.

［18］郝愛民， 譚家銀. 數(shù)字鄉(xiāng)村建設對我國糧食體系韌性的影響［J］. 華南農(nóng)業(yè)大學學報（社會科學版）， 2022（3）： 10-24.Hao Aimin， Tan Jiayin. Impact of digital rural construction on food system resilience ［J］. Journal of South China Agricultural University （Social Science Edition）， 2022（3）： 10-24.

［19］谷城， 張樹山. 產(chǎn)業(yè)鏈韌性水平測度、區(qū)域差異及收斂性研究［J］. 經(jīng)濟問題探索， 2023（6）： 123-139.Gu Cheng， Zhang Shushan. Level measurement， regional differences and convergence study of the industrial chain resilience ［J］. Inquiry into Economic Issues， 2023（6）： 123-139.

［20］張紅麗， 李潔艷. 農(nóng)業(yè)技術進步、農(nóng)村勞動力轉(zhuǎn)移與城鄉(xiāng)收入差距——基于農(nóng)業(yè)勞動生產(chǎn)率的分組研究［J］. 華東經(jīng)濟管理， 2020（1）： 67-75.Zhang Hongli， Li Jieyan. The impact of agricultural technology progress， rural labor transfer and the income gap between urban and rural areas： Study and analysis on the grouping of agricultural labor productivity ［J］. East China Economic Management， 2020（1）： 67-75.

［21］鄭晶， 高孟菲. 農(nóng)業(yè)機械化、農(nóng)村勞動力轉(zhuǎn)移對農(nóng)業(yè)全要素生產(chǎn)率的影響研究——基于中國大陸31個?。ㄊ小⒆灾螀^(qū)）面板數(shù)據(jù)的實證檢驗［J］. 福建論壇（人文社會科學版）， 2021（8）： 59-71.

［22］黃楓， 孫世龍. 讓市場配置農(nóng)地資源： 勞動力轉(zhuǎn)移與農(nóng)地使用權市場發(fā)育［J］. 管理世界， 2015（7）： 71-81.

基金項目：國家社科基金重點項目（20AJY012）；國家社科基金項目（17CJY040）；湖南省社會科學基金項目（21JD020）

第一作者：蔡林軍，男，1999年生，湖南桂陽人，碩士研究生；研究方向為糧食安全。E-mail： 3606289531@qq.com

通訊作者：文春暉，男，1981年生，湖南安化人，博士，教授，博導；研究方向為農(nóng)業(yè)經(jīng)濟。E-mail： 2414421393@qq.com