資源稟賦、勞動力轉(zhuǎn)移對耕地利用功能轉(zhuǎn)型的影響研究
2023-10-24 15:17:43陳莎楊潤佳葉艷妹薛繼兵
中國土地科學(xué) 2023年2期
陳莎 楊潤佳 葉艷妹 薛繼兵
摘要：研究目的：探討耕地資源稟賦、農(nóng)業(yè)勞動力轉(zhuǎn)移對耕地利用功能轉(zhuǎn)型的影響機制，推動耕地利用轉(zhuǎn)型和多功能管理。研究方法：面板門檻模型分析。研究結(jié)果：（1）資源稟賦、農(nóng)業(yè)勞動力轉(zhuǎn)移對耕地利用功能綜合水平和耕地利用功能形態(tài)均存在門檻效應(yīng)，長期來看勞動力轉(zhuǎn)移對耕地利用功能綜合水平的提升趨于有利，對耕地利用功能形態(tài)的影響呈現(xiàn)從“正向到負向再到負向略有緩解”的變化特征；（2）當(dāng)納入耕地資源稟賦考慮時，勞動力轉(zhuǎn)移對耕地利用功能綜合水平和耕地利用功能形態(tài)的影響從“負”轉(zhuǎn)變?yōu)椤罢保译S著耕地資源稟賦的上升正向影響不斷凸顯。研究結(jié)論：要因勢利導(dǎo)，把握本地區(qū)勞動力轉(zhuǎn)移的發(fā)展趨勢，通過農(nóng)業(yè)技術(shù)創(chuàng)新和農(nóng)業(yè)生產(chǎn)政策的完善合理調(diào)控勞動力結(jié)構(gòu)及農(nóng)業(yè)要素投入。同時，要善于發(fā)揮本地區(qū)耕地資源稟賦的優(yōu)勢并認識其與社會經(jīng)濟發(fā)展的內(nèi)在耦合性和階段差異性，在勞動力轉(zhuǎn)移過程中加強全域土地綜合整治和生態(tài)保護修復(fù)，推動適度規(guī)模經(jīng)營并挖掘發(fā)展特色農(nóng)業(yè)產(chǎn)業(yè)，以實現(xiàn)耕地利用轉(zhuǎn)型和功能優(yōu)化升級。
關(guān)鍵詞：耕地功能轉(zhuǎn)型；面板門檻回歸；耕地可持續(xù)利用；功能形態(tài)
中圖分類號：F301.2 文獻標志碼：A 文章編號：1001-8158（2023）02-0071-11
基金項目：浙江省哲學(xué)社會科學(xué)規(guī)劃課題項目（20NDJC140YB）。
耕地利用關(guān)乎全球糧食安全和可持續(xù)發(fā)展目標的實現(xiàn)，一直是國內(nèi)外學(xué)者和社會各界關(guān)注的焦點[1]。當(dāng)前，我國人口眾多、人均耕地資源稀缺的基本地情沒有改變，全球氣候變化帶來的風(fēng)險干擾致使區(qū)域農(nóng)業(yè)自然資源脆弱性不斷增強，快速城鎮(zhèn)化、工業(yè)化出現(xiàn)的耕地非農(nóng)化、鄉(xiāng)村主體老弱化、農(nóng)村空廢化、環(huán)境污損化等“鄉(xiāng)村病”問題制約著鄉(xiāng)村振興戰(zhàn)略和耕地可持續(xù)利用[2]。促進耕地利用功能擴展、升級使之與經(jīng)濟社會發(fā)展相適應(yīng)，對耕地資源的高效、生態(tài)和可持續(xù)利用意義重大。
土地利用轉(zhuǎn)型被視為土地利用/覆蓋變化（LUCC）研究的新視角，自21世紀初被引入中國以來，在概念內(nèi)涵、理論假說與框架以及區(qū)域?qū)嵶C等方面都取得了長足的進展[3-5]。耕地利用轉(zhuǎn)型作為土地利用轉(zhuǎn)型研究的重要部分，在概念與內(nèi)涵上延續(xù)了土地利用轉(zhuǎn)型的研究成果，其實質(zhì)是耕地利用形態(tài)在長時間序列變化下的趨勢性轉(zhuǎn)折[6]，目前該領(lǐng)域的研究主要集中在耕地利用轉(zhuǎn)型研究框架[7]、轉(zhuǎn)型路徑[8]、驅(qū)動機制[9-10]及其對農(nóng)業(yè)生產(chǎn)[11]、農(nóng)業(yè)經(jīng)濟增長[12]的影響等方面，涉及生態(tài)環(huán)境響應(yīng)、城鄉(xiāng)關(guān)系互動、鄉(xiāng)村轉(zhuǎn)型發(fā)展等多個重要議題和領(lǐng)域[13-15]。一般認為，耕地利用轉(zhuǎn)型可以從空間轉(zhuǎn)型（包括數(shù)量、景觀格局、經(jīng)營格局方面的變化）和功能轉(zhuǎn)型（耕地提供的產(chǎn)品和服務(wù)功能方面的變化）兩方面進行診斷[4]，進一步地，耕地利用功能轉(zhuǎn)型又可以從耕地功能綜合水平和功能形態(tài)兩方面進行探討[16]。已有研究多關(guān)注耕地的經(jīng)濟貢獻、糧食生產(chǎn)、社會保障、生態(tài)服務(wù)等功能，通過指標替代法[17]、物質(zhì)量法[18]、價值量法[19]等對耕地的功能量進行核算，運用相關(guān)分析[20]、均方根誤差[21]、耦合協(xié)調(diào)度[17]等方法探討耕地多功能協(xié)同或者沖突的結(jié)構(gòu)性特征，但從耕地利用功能綜合水平和功能形態(tài)兩個方面探討耕地利用功能轉(zhuǎn)型及其內(nèi)在機制的研究還鮮有報道。
耕地系統(tǒng)可被視為由自然生態(tài)子系統(tǒng)與社會經(jīng)濟子系統(tǒng)耦合互動而成的復(fù)雜系統(tǒng)，受到人類活動、社會經(jīng)濟發(fā)展、管理制度等多因素的影響和制約[22]，耕地利用功能演變受自然稟賦[23]、社會經(jīng)濟發(fā)展[24]、政策干預(yù)等因素的影響[25]，在當(dāng)前緊張的人地關(guān)系下，農(nóng)業(yè)人口轉(zhuǎn)移和耕地資源稟賦是影響耕地利用功能轉(zhuǎn)型的關(guān)鍵因素和重要驅(qū)動力[11，26]。此外，已有研究多采用基于線性相關(guān)假設(shè)的多元回歸方法探討各項驅(qū)動因素與耕地功能水平的關(guān)系，忽視了耕地利用功能轉(zhuǎn)型可能會依賴于耕地資源稟賦條件和勞動力非農(nóng)轉(zhuǎn)移而呈現(xiàn)非線性關(guān)系，而以“門檻回歸”為代表的非線性計量方法為這一問題提供了恰當(dāng)?shù)姆治龉ぞ撸谙嚓P(guān)領(lǐng)域研究中體現(xiàn)出一定的優(yōu)勢[27]。
浙江省屬于經(jīng)濟較發(fā)達地區(qū)，但耕地資源少，分布散，且近年來表現(xiàn)出空間分離度增加趨勢[28]。另一方面，浙江省憑借市場經(jīng)濟的先發(fā)優(yōu)勢，有序推進農(nóng)業(yè)轉(zhuǎn)移人口市民化，近年來農(nóng)業(yè)轉(zhuǎn)移人口使浙江城鎮(zhèn)化率提高了30多個百分點①[29]。一些學(xué)者對浙江省耕地多功能時空演變特征[30]和權(quán)衡關(guān)系[31]等進行了分析，但尚未從耕地功能綜合水平和功能結(jié)構(gòu)形態(tài)兩個方面探討耕地利用功能轉(zhuǎn)型的機制，特別是耕地資源稟賦和農(nóng)業(yè)人口轉(zhuǎn)移對其產(chǎn)生的影響。鑒于此，針對當(dāng)前研究不足和耕地利用功能轉(zhuǎn)型的重要性，本文在分析耕地利用功能轉(zhuǎn)型受耕地稟賦、勞動力轉(zhuǎn)移影響機制的基礎(chǔ)上，運用面板門檻模型，分別以勞動力轉(zhuǎn)移和耕地資源稟賦為門檻變量，考察浙江省69個縣區(qū)耕地利用功能轉(zhuǎn)型的非線性演變特征。研究結(jié)論可為完善耕地可持續(xù)利用和糧食安全政策提供決策參考。
1 耕地利用功能轉(zhuǎn)型的理論解釋
1.1 耕地利用功能轉(zhuǎn)型的內(nèi)涵
從耕地利用功能的“整體性”和“結(jié)構(gòu)性”內(nèi)涵出發(fā)，本文認為耕地利用功能轉(zhuǎn)型有兩層含義：第一，凝結(jié)在耕地資源上的生態(tài)系統(tǒng)綜合服務(wù)出現(xiàn)趨勢性變化轉(zhuǎn)折，即在宏觀政策和消費需求的誘導(dǎo)下，達到生產(chǎn)有效、環(huán)境友好、人地和諧的目標，提升耕地功能的綜合價值；第二，側(cè)重耕地功能結(jié)構(gòu)形態(tài)的轉(zhuǎn)變，關(guān)注耕地利用模式、空間結(jié)構(gòu)、耕地功能的此強彼弱方面的根本性變化。依據(jù)“誘致性生產(chǎn)替代理論”[32]，隨著社會經(jīng)濟不同發(fā)展階段的演進，耕地利用投入要素的相對價格和產(chǎn)品（服務(wù)）的相對價格會出現(xiàn)變化，進而導(dǎo)致耕地各項功能的相對強弱不同，呈現(xiàn)差異性的耕地利用功能形態(tài)。
鑒于此，本文將耕地利用功能轉(zhuǎn)型劃分為功能綜合水平和功能形態(tài)兩個方面。耕地功能綜合水平，即探究耕地利用產(chǎn)生的“共生功能束”的綜合價值，一般考慮耕地功能的生產(chǎn)、生態(tài)、社會和文化等功能之和[31]；耕地功能形態(tài)，隨著社會經(jīng)濟發(fā)展從初級向中高級階段演變，耕地利用功能形態(tài)表現(xiàn)出從“社會生態(tài)型”到“經(jīng)濟社會型”再到“生態(tài)經(jīng)濟型”的轉(zhuǎn)型特征[32]，大量研究證據(jù)表明[20，33]，隨著社會經(jīng)濟水平的提升，耕地的作物生產(chǎn)和社會保障功能的貢獻度有所降低，文化休閑功能逐步增長，生態(tài)調(diào)節(jié)功能出現(xiàn)“恢復(fù)性”提升，據(jù)此本文認為生態(tài)調(diào)節(jié)、環(huán)境保育和文化休閑功能是較高階段的主導(dǎo)功能，物質(zhì)生產(chǎn)功能、社會保障功能是較低階段的主導(dǎo)功能。
1.2 耕地稟賦、勞動力轉(zhuǎn)移對耕地利用功能轉(zhuǎn)型的影響機制分析
1.2.1 勞動力轉(zhuǎn)移對耕地利用功能轉(zhuǎn)型的影響機制分析
目前關(guān)于農(nóng)業(yè)勞動力析出對耕地利用的影響研究主要關(guān)注耕地利用集約度，根據(jù)“粗放化假說”，當(dāng)大量農(nóng)業(yè)勞動力轉(zhuǎn)移到競爭性行業(yè)，勞動力成本上升，會造成農(nóng)業(yè)收縮，土地利用或表現(xiàn)出集約度下降即粗放經(jīng)營，或面積收縮即退出生產(chǎn)[34]。粗放利用、拋荒棄耕會導(dǎo)致耕地的物質(zhì)生產(chǎn)功能下降，由于減少了化肥施用減少了碳排放，同時耕地作為農(nóng)村生活保障的地位重要性明顯削弱。而根據(jù)“集約化假說”，留守勞動力仍承擔(dān)較多的農(nóng)業(yè)生產(chǎn)活動，依賴匯款等非農(nóng)收入通過追加化肥等投資來替代外出的勞動損失[35]，能夠繼續(xù)維持或提升農(nóng)業(yè)生產(chǎn)，也有學(xué)者認為耕地集約將經(jīng)歷勞動集約、資金集約、技術(shù)集約的持續(xù)性、階段性轉(zhuǎn)變[36]，故認為勞動力轉(zhuǎn)移對耕地物質(zhì)產(chǎn)出功能沒有明顯影響甚至有所提升，但由于化肥等農(nóng)業(yè)物質(zhì)的大量施用會降低耕地的凈固碳能力，可能產(chǎn)生農(nóng)村生態(tài)環(huán)境的負外部性。另一方面，勞動力轉(zhuǎn)移有助于耕地流轉(zhuǎn)從而促進適度規(guī)模經(jīng)營[37]，而集中連片經(jīng)營的耕地更有可能提供良好的景觀美學(xué)體驗進而提升社會福祉[38]。
1.2.2 耕地稟賦對耕地利用功能轉(zhuǎn)型的影響機制分析
本文在勞動力轉(zhuǎn)移對耕地利用功能轉(zhuǎn)型影響中納入對耕地稟賦條件的考慮，影響機制見圖1。在勞動力轉(zhuǎn)移的過程中，對耕地資源稟賦條件較好的地區(qū)而言，更容易促成生產(chǎn)要素轉(zhuǎn)型升級和有效運行，通過流轉(zhuǎn)經(jīng)營、政策支持、產(chǎn)業(yè)升級等方式影響耕地利用功能綜合水平和耕地利用功能形態(tài)。具體來看，一方面，高耕地稟賦地區(qū)通過規(guī)模經(jīng)營、集中整理等能有效實現(xiàn)規(guī)模經(jīng)營和集約利用，減少由于耕地細碎化、邊際化所致的拋荒、閑置現(xiàn)象；另一方面，耕地資源稟賦優(yōu)越的地區(qū)往往具有更高的農(nóng)業(yè)產(chǎn)出，在國家嚴格耕地保護的主題背景下，耕地資源稟賦突出的地區(qū)擁有顯著的農(nóng)業(yè)發(fā)展戰(zhàn)略優(yōu)勢，將影響當(dāng)?shù)匾幌盗械母乇Ｗo和農(nóng)業(yè)支持政策的制定，有效推動產(chǎn)業(yè)結(jié)構(gòu)升級、生產(chǎn)規(guī)?；娃r(nóng)業(yè)效益提升[39-40]；此外，隨著社會經(jīng)濟水平的提升，公眾對于綠色生態(tài)產(chǎn)品的意識和需求提升，擁有良好耕地稟賦的地區(qū)具備發(fā)展景觀農(nóng)業(yè)、休閑農(nóng)業(yè)、生態(tài)農(nóng)業(yè)等新型業(yè)態(tài)的天然優(yōu)勢，有利于推動耕地功能形態(tài)轉(zhuǎn)型。
2 研究方法
2.1 變量選取與處理
2.1.1 被解釋變量
2.1.2 解釋變量
本文的解釋變量是基于具體研究目的和相關(guān)研究進行篩選確定[11，30，48]。以各區(qū)縣人均耕地面積（pCLA）來衡量耕地資源稟賦。以勞動力非農(nóng)就業(yè)水平（ALT）反映勞動力轉(zhuǎn)移程度，以上二者為核心解釋變量。同時結(jié)合可拓展隨機性環(huán)境影響評估（STIRPAT）模型，以人均生產(chǎn)總值（GDPP）衡量富裕程度，以農(nóng)業(yè)機械投入強度（MII）衡量技術(shù)水平，并引入復(fù)種指數(shù)（MCI）、非糧作物播種面積占比（NGP）和財政支農(nóng)水平（FSA）3項因素對標準STIRPAT模型進行擴展。表1 給出了所有解釋變量的測算方式和描述性統(tǒng)計，偏度和峰度計算結(jié)果顯示大部分變量原始數(shù)據(jù)不符合正態(tài)分布。
2.2 耕地利用功能轉(zhuǎn)型的門檻回歸模型
門檻效應(yīng)是指當(dāng)某一個經(jīng)濟參數(shù)達到特定的臨界值后，引起的另一個經(jīng)濟參數(shù)發(fā)生方向上或者數(shù)量上變化的結(jié)構(gòu)突變現(xiàn)象，該經(jīng)濟參數(shù)的臨界值即為門檻值[31]。門檻回歸是確定導(dǎo)致經(jīng)濟系統(tǒng)結(jié)構(gòu)發(fā)生突變的臨界值的有效工具，并一定程度上避免了主觀判斷臨界值導(dǎo)致的統(tǒng)計誤差和估計偏誤[48]。為探索農(nóng)村勞動力轉(zhuǎn)移和耕地資源稟賦對耕地利用功能轉(zhuǎn)型影響的門檻效應(yīng)，設(shè)定面板門檻回歸模型如下：
以一門檻值模型為基礎(chǔ)，還可以擴展得到兩個及以上門檻的多門檻值模型。兩門檻模型計算過程與一門檻類似，是在第一個門檻值確定的情形下，估計第二個門檻值。實證分析中常采用多門檻值模型進行門檻效應(yīng)檢驗并估計具體門檻值。
2.3 數(shù)據(jù)來源和預(yù)處理
本文以2000年、2005年、2010年、2015年浙江省69個區(qū)縣的面板數(shù)據(jù)為研究樣本，耕地面積數(shù)據(jù)來源于浙江省歷年土地利用變更調(diào)查數(shù)據(jù)庫，耕地利用功能測度過程和影響因素所涉及的農(nóng)業(yè)投入產(chǎn)出數(shù)據(jù)、人口數(shù)據(jù)、社會經(jīng)濟數(shù)據(jù)主要來源于歷年的《浙江統(tǒng)計年鑒》，部分來源于各地級市統(tǒng)計年鑒和各區(qū)縣統(tǒng)計局，氣象數(shù)據(jù)來源于中國氣象數(shù)據(jù)服務(wù)共享網(wǎng)（http：//data.cma.cn），土壤深度與孔隙率數(shù)據(jù)來源于1∶100萬中國土壤數(shù)據(jù)庫；縣域中心距最近地級市轄區(qū)中心直線距離數(shù)據(jù)通過百度地圖的測距功能獲取。耕地斑塊聚集度、耕地景觀多樣性數(shù)據(jù)通過Fragstats 4.2軟件計算得到。
數(shù)據(jù)預(yù)處理包括：第一，針對耕地面積數(shù)據(jù)，為了消除第二次全國土地調(diào)查前后數(shù)據(jù)突變，將第二次全國土地調(diào)查公布的 2009年耕地面積減去當(dāng)年土地利用變更調(diào)查中的耕地面積，將差值加到2000年、2005年耕地面積上，完成了對數(shù)據(jù)的校正。第二，針對部分區(qū)縣個別年份個別數(shù)據(jù)的缺失，以歷年數(shù)據(jù)為基礎(chǔ)通過移動平均、趨勢外推等方法進行數(shù)據(jù)重建。第三，為消除價格指數(shù)對農(nóng)業(yè)產(chǎn)值的影響，將各年的農(nóng)業(yè)產(chǎn)值轉(zhuǎn)化為2000年可比價。第四，為了使統(tǒng)計數(shù)據(jù)范圍與現(xiàn)有行政區(qū)范圍一致，本文對部分行政區(qū)作市轄區(qū)合并處理，由于舟山市嵊泗縣耕地面積較小，暫不納入該地區(qū)，最終以69 個區(qū)縣為研究樣本。
3 結(jié)果與分析
3.1 面板回歸結(jié)果
（1）變量多重共線性檢驗。運用方差膨脹因子VIF分析對文中的各解釋變量進行多重共線性因子檢驗，方差膨脹因子VIF檢驗結(jié)果見表2。結(jié)果顯示，所有變量中最大的VIF值為3.40，遠小于10。因此，結(jié)合相關(guān)系數(shù)和方差膨脹因子VIF的分析結(jié)果，可認為變量間不存在多重共線性。
（2）面板數(shù)據(jù)回歸結(jié)果。研究運用多種回歸方法分析變量間的數(shù)量關(guān)系，包括混合回歸（Pooled）、固定效應(yīng)模型（FE）和隨機效應(yīng)模型（RE）。為了得到更準確的估計結(jié)果，采用LM 檢驗和Hausman 檢驗從上述三個模型中選擇最優(yōu)模型。本文分別以耕地利用功能綜合水平（MFCI）和耕地利用功能形態(tài)（FSI）為被解釋變量，回歸結(jié)果見表3。根據(jù)F檢驗、LM檢驗、Hausman檢驗的結(jié)果，本文應(yīng)采用“固定效應(yīng)模型（FE）”進行分析。
（3）耕地利用功能綜合水平回歸結(jié)果。當(dāng)被解釋變量為耕地利用功能綜合水平lnMFCI時，lnpCLA、lnALT的回歸系數(shù)均為負值，說明資源稟賦、農(nóng)業(yè)勞動力轉(zhuǎn)移對耕地利用功能綜合水平產(chǎn)生了負向影響。lnMCI、lnNGP、lnFSA的回歸系數(shù)均為正，說明當(dāng)復(fù)種指數(shù)越高、非糧作物種植比例越高以及財政支農(nóng)水平越高時，會促進耕地利用功能綜合水平提升；人均生產(chǎn)總值lnGDPP、農(nóng)業(yè)機械投入強度lnMII對耕地利用功能綜合水平的影響均未通過顯著性檢驗。
（4）耕地利用功能形態(tài)回歸結(jié)果。當(dāng)被解釋變量為耕地利用功能形態(tài)lnFSI時，與lnMFCI類似的是，lnpCLA、lnALT的回歸系數(shù)均為負，說明資源稟賦、農(nóng)業(yè)勞動力轉(zhuǎn)移對耕地利用功能形態(tài)轉(zhuǎn)型產(chǎn)生了負向影響。且lnpCLA、lnALT回歸系數(shù)的絕對值遠大于被解釋變量為lnMFCI的回歸系數(shù)。非糧作物種植面積占比lnNGP對耕地功能形態(tài)高級化產(chǎn)生了正向影響，lnGDPP、lnMII、lnMCI和lnFSA對耕地功能形態(tài)轉(zhuǎn)型的影響均未通過顯著性檢驗?？赡艿慕忉屓缦拢阂皇悄壳耙匝a充耕地面積為主要目標的耕地保護投入（如高標準農(nóng)田建設(shè)、耕地整治等）不足以改變以生產(chǎn)功能為主的耕地功能結(jié)構(gòu)形態(tài)，二是“數(shù)量、質(zhì)量、生態(tài)”三位一體的耕地保護成效可能具有滯后性，有待后續(xù)研究進一步檢驗。
3.2 耕地利用功能綜合水平面板門檻回歸結(jié)果
采用面板門檻回歸模型對變量間是否存在非線性關(guān)系進行檢驗。當(dāng)耕地利用功能綜合水平lnMFCI為被解釋變量時，分別以lnALT和lnpCLA為門檻變量，對不存在門檻值、存在單一門檻值以及兩個門檻值分別進行檢驗，直到結(jié)果不顯著為止。運用Stata15.1統(tǒng)計軟件，通過 “自助法”（Bootstrap）反復(fù)抽樣300次從而得出檢驗統(tǒng)計量對應(yīng)的P值，判斷是否存在門檻效應(yīng)，檢驗結(jié)果見表4。
表4結(jié)果顯示，當(dāng)lnALT為門檻變量時，F(xiàn)統(tǒng)計量在一門檻模型中，至少在5%的水平下顯著，表明此模型中存在一個門檻值，門檻值的估計結(jié)果為-0.910 3，95%置信區(qū)間為（ -1.010 8， -0.856 8）。圖（2a）為該門檻值的似然比函數(shù)圖。當(dāng)lnpCLA為門檻變量時，F(xiàn)統(tǒng)計量在一門檻和二門檻的模型中均在10%水平下顯著，因此模型中存在兩個門檻值，門檻值的估計結(jié)果分別為5.673 8和6.636 3。圖2（b）和圖2（c）分別為2個門檻值5.673 8、6.636 3的似然比函數(shù)圖。似然比函數(shù)圖可以更為清晰地理解門檻值的估計和置信區(qū)間的構(gòu)造過程，圖中虛線為LR值在5% 顯著水平下的臨界值，虛線以下的區(qū)域構(gòu)成門檻值的95% 置信區(qū)間。圖2（b）和圖2（c）顯示LR值在95% 漸進有效置信區(qū)間[5.670 8，5.674 1]和[6.632 5，6.640 5]內(nèi)接近于零，檢驗結(jié)果無法拒絕門檻估計值為其真實值一致估計量的原假設(shè)，由此可斷定模型估計存在雙重門檻效應(yīng)。
面板門檻模型的參數(shù)估計結(jié)果見表5。當(dāng)農(nóng)業(yè)勞動力轉(zhuǎn)移處于初級階段時（lnALT≤-0.910 3），其對耕地利用功能綜合水平的影響因素為-0.101；當(dāng)農(nóng)業(yè)勞動力轉(zhuǎn)移處于較高水平時（lnALT>-0.910 3），回歸系數(shù)變?yōu)?.222，且都至少在5%水平下顯著。由此可知，農(nóng)業(yè)勞動力轉(zhuǎn)移對耕地利用功能綜合水平的不利影響會隨著農(nóng)業(yè)勞動力轉(zhuǎn)移的推進而逐漸下降，當(dāng)越過門檻值-0.910 3時，會促進耕地利用功能綜合水平提升。其余控制變量由于與面板數(shù)據(jù)回歸結(jié)果（表3）差別不大，此處不再贅述，下同。
以lnpCLA為門檻變量，當(dāng)lnpCLA≤5.673 8時，勞動力轉(zhuǎn)移水平lnALT對耕地利用功能綜合水平提升的影響系數(shù)為-0.163，并在1%水平下通過顯著性檢驗；當(dāng)5.673 8
3.3 功能結(jié)構(gòu)形態(tài)面板門檻回歸結(jié)果
以耕地利用功能形態(tài)lnFSI為被解釋變量，分別對lnAPT和lnpCLA是否存在門檻效應(yīng)進行檢驗。檢驗結(jié)果見表6。
由表6可知，當(dāng)分別以lnALT和lnpCLA為門檻變量時，均可得到以下結(jié)論，F(xiàn)統(tǒng)計量在一門檻和二門檻模型檢驗中，均至少在10%水平下顯著，因此模型中存在兩個門檻值。表7報告了門檻值估計結(jié)果。與表7相對應(yīng)，圖3和圖4分別為lnALT兩個門檻值-0.917 6、-0.199 6和lnpCLA兩個門檻值6.028 5、6.621 8在95%置信區(qū)間下的似然比函數(shù)圖，2個門檻值均位于臨界值7.35虛線以下的區(qū)域構(gòu)成的95% 置信區(qū)間內(nèi)，由此可以認為上述門檻值是真實有效的。
表8報告了面板門檻回歸的參數(shù)估計結(jié)果。當(dāng)農(nóng)業(yè)勞動力轉(zhuǎn)移處于初級階段時（lnALT≤-0.917 6），其對耕地利用功能形態(tài)的影響因素為0.267；當(dāng)農(nóng)業(yè)勞動力轉(zhuǎn)移處于中高水平時（-0.917 6
以lnpCLA為門檻變量時，當(dāng)lnpCLA≤6.028 5時，勞動力轉(zhuǎn)移水平lnALT對耕地利用功能形態(tài)提升的影響系數(shù)為-0.526，并在1%水平下通過顯著性檢驗；當(dāng)6.028 5
4 結(jié)論與討論
4.1 結(jié)論
本文將耕地利用功能轉(zhuǎn)型劃分為功能綜合水平和功能形態(tài)兩個方面，通過面板門檻回歸方法分析了勞動力轉(zhuǎn)移、耕地資源稟賦對耕地利用功能轉(zhuǎn)型的門檻特征以及各變量的影響程度。主要結(jié)論如下：第一，無論是耕地利用功能綜合水平還是功能形態(tài)，勞動力轉(zhuǎn)移水平、耕地資源稟賦均對其表現(xiàn)出顯著的門檻效應(yīng)；第二，以勞動力轉(zhuǎn)移為門檻變量時，勞動力轉(zhuǎn)移水平的增長對耕地利用功能綜合水平的提升趨于有利，對耕地利用功能形態(tài)轉(zhuǎn)型的影響呈現(xiàn)出從“正向到負向再到負向略有緩解”的轉(zhuǎn)變過程，可認為勞動力轉(zhuǎn)移是推進區(qū)域耕地利用程度、產(chǎn)出效益和利用方式轉(zhuǎn)變的基礎(chǔ)，但由于耕地各項功能之間的變化程度不同，即便同為增長其增長程度也并不同步，導(dǎo)致耕地功能形態(tài)轉(zhuǎn)型受勞動力轉(zhuǎn)移水平的影響表現(xiàn)出更多樣的區(qū)間差異性；第三，當(dāng)納入耕地資源稟賦考慮時，勞動力轉(zhuǎn)移水平對耕地利用功能綜合水平和耕地功能形態(tài)的影響均趨于有利，隨著耕地資源稟賦的上升，會不斷促進耕地利用功能綜合水平提升和功能形態(tài)高級化。
4.2 討論
研究結(jié)論對于完善相關(guān)政策具有重要啟示：首先，當(dāng)前仍有大量的農(nóng)村勞動力向外轉(zhuǎn)移，地方政府應(yīng)通過農(nóng)業(yè)技術(shù)創(chuàng)新、土地利用及農(nóng)業(yè)支持政策的創(chuàng)新，合理調(diào)控勞動力結(jié)構(gòu)及農(nóng)業(yè)要素投入，促進農(nóng)村勞動力收入需求、農(nóng)業(yè)生產(chǎn)與生態(tài)安全等目標的相互協(xié)調(diào)；其次，要善于發(fā)揮本地區(qū)耕地資源稟賦的優(yōu)勢并認識其與社會經(jīng)濟發(fā)展的內(nèi)在耦合性和階段差異性，一方面，城市化背景下仍需有效控制耕地流失，通過耕地“占補平衡”“進出平衡”守牢耕地穩(wěn)定利用紅線；另一方面，應(yīng)認識到以補充耕地面積為主要目標的耕地保護不足以改變以生產(chǎn)功能為主導(dǎo)的耕地功能形態(tài)，需要加強山水林田湖草系統(tǒng)治理邏輯下的耕地資源優(yōu)化配置，通過全域土地綜合整治和生態(tài)保護修復(fù)在耕地功能總量增長的同時，優(yōu)化功能形態(tài)并促進功能的多樣化協(xié)同。此外，對于勞動力轉(zhuǎn)移過程中頻發(fā)的耕地撂荒、耕地利用邊際化問題，應(yīng)積極推動農(nóng)戶復(fù)耕、規(guī)模流轉(zhuǎn)、委托代管以及集體助耕等多途徑治理。通過培育新型農(nóng)業(yè)經(jīng)營主體，促進適度規(guī)模經(jīng)營并挖掘發(fā)展特色農(nóng)業(yè)產(chǎn)業(yè)，推動農(nóng)業(yè)要素資源重新配置和耕地資源稟賦價值轉(zhuǎn)化，充分發(fā)揮農(nóng)村勞動力轉(zhuǎn)移、耕地資源稟賦對提升耕地利用功能綜合水平和功能形態(tài)的合力作用。
本文還存在一些不足。首先，本文所采用的實際衡量指標存在一定的偏差，農(nóng)業(yè)勞動力轉(zhuǎn)移原則上應(yīng)為農(nóng)村的農(nóng)業(yè)人口轉(zhuǎn)移到非農(nóng)產(chǎn)業(yè)的人數(shù)，本文出于數(shù)據(jù)處理的便捷考慮，采用勞動力非農(nóng)就業(yè)水平來簡化替代。其次，人均耕地面積只是耕地資源稟賦的一個方面，其他因素如耕地質(zhì)量、連片度、地形、光溫條件等本文并未納入，以上因素與人文和社會經(jīng)濟要素的互動關(guān)系并進而引起耕地利用功能轉(zhuǎn)型的內(nèi)在機制是下一步研究的重點。最后，本文重點關(guān)注了勞動力轉(zhuǎn)移對耕地利用功能轉(zhuǎn)型的影響，但耕地功能的結(jié)構(gòu)性改變亦可能對勞動力轉(zhuǎn)移產(chǎn)生“反作用”，例如耕地經(jīng)濟功能的強化可能對勞動力產(chǎn)生擠出作用，但就目前浙江省耕地利用功能轉(zhuǎn)型的特點和勞動力轉(zhuǎn)移的水平來看，耕地利用功能形態(tài)的變化總體上較為緩慢，勞動力轉(zhuǎn)移依然主要受非農(nóng)就業(yè)的比較效益驅(qū)動，或許尚未進入到“反作用”產(chǎn)生的階段，因此在本文中暫不考慮，但這一問題確實有待進一步探討。
參考文獻（References）：
[1] VIANA C M， FREIRE D， ABRANTES P， et al. Agricultural land systems importance for supporting food security and sustainable development goals： a systematic review[J] . Science of the Total Environment， 2022， 806： 150718. doi： 10.1016/j.scitotenv.2021.150718.
[2] ZHOU Y， LI X H， LIU Y S. Cultivated land protection and rational use in China[J] . Land Use Policy， 2021， 106： 105454. doi： 10.1016/j.landusepol.2021.105454.
[3] 龍花樓. 土地利用轉(zhuǎn)型的解釋[J] . 中國土地科學(xué)，2022，36（4）：1 - 7.
[4] 宋小青. 論土地利用轉(zhuǎn)型的研究框架[J] . 地理學(xué)報，2017，72（3）：471 - 487.
[5] 田俊峰， 王彬燕， 程利莎， 等. 政策主導(dǎo)下的區(qū)域土地利用轉(zhuǎn)型過程與機制——以中國東北地區(qū)為例[J] . 地理研究，2020，39（4）：805 - 821.
[6] 蔣夢凡， 李智國， 李杰， 等. 耕地利用轉(zhuǎn)型突變點檢測及其時空特征——以昆明市東川區(qū)為例[J] . 中國土地科學(xué)，2022，36（3）：86 - 95.
[7] 曲藝， 龍花樓. 中國耕地利用隱性形態(tài)轉(zhuǎn)型的多學(xué)科綜合研究框架[J] . 地理學(xué)報，2018，73（7）：1226 - 1241.
[8] 吳宇哲， 沈欣言. 中國耕地保護治理轉(zhuǎn)型：供給、管制與賦能[J] . 中國土地科學(xué)，2021，35（8）：32 - 38.
[9] 史洋洋， 呂曉， 郭貫成， 等. 基于GIS和空間計量的耕地利用轉(zhuǎn)型時空格局及其驅(qū)動機制研究[J] . 中國土地科學(xué)，2019，33（11）：51 - 60.
[10] 柯善淦， 崔?，摚?盧新海， 等. 耕地利用綠色轉(zhuǎn)型的時空格局及其驅(qū)動機制研究——以湖北省為例[J] . 中國土地科學(xué)，2021，35（12）：64 - 74.
[11] 劉光盛， 王紅梅， 盧德偉， 等. 農(nóng)村勞動力析出及稟賦依賴對耕地集約利用的影響[J] . 中國農(nóng)業(yè)資源與區(qū)劃，2021，42（7）：58 - 65.
[12] 盧新海， 唐一峰， 易家林， 等. 基于空間計量模型的耕地利用轉(zhuǎn)型對農(nóng)業(yè)經(jīng)濟增長影響研究[J] . 中國土地科學(xué)，2019，33（6）：53 - 61.
[13] 張茜茜， 廖和平， 巫芯宇， 等. 鄉(xiāng)村振興背景下的“人、地、業(yè)”轉(zhuǎn)型空間差異及影響因素分析——以重慶市渝北區(qū)為例[J] . 西南大學(xué)學(xué)報（自然科學(xué)版），2019，41（4）：1 - 9.
[14] 趙寧寧， 郭炎， 李志剛， 等. 中部地區(qū)鄉(xiāng)村生產(chǎn)要素協(xié)同轉(zhuǎn)型的時空格局及影響因素——以湖南岳陽為例[J] . 自然資源學(xué)報，2021，36（12）：3170 - 3185.
[15] SONG B J， ROBINSON G M， ZHOU Z X. Agricultural transformation and ecosystem services： a case study from Shaanxi Province， China[J] . Habitat International， 2017， 69： 114 - 125.
[16] 胡偉艷， 朱慶瑩， 張安錄， 等. 總量與結(jié)構(gòu)視角耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響——以湖北省為例[J] . 中國土地科學(xué)，2018，32（5）：62 - 70.
[17] JIANG G H， WANG M Z， QU Y B， et al. Towards cultivated land multifunction assessment in China： applying the“influencing factors-functions-products-demands”integrated framework[J] . Land Use Policy， 2020， 99： 104982. doi： 1016/j. landusepol. 2020.104982.
[18] TURNER K G， ODGAARD M V， B CHER P K， et al. Bundling ecosystem services in Denmark： trade-offs and synergies in a cultural landscape[J] . Landscape and Urban Planning， 2014， 125： 89 - 104.
[19] 高進云， 楊微， 喬榮鋒. 天津市郊區(qū)休閑農(nóng)莊農(nóng)地游憩價值評估——以西青區(qū)楊柳青莊園和水高莊園為例[J] .資源科學(xué)，2014，36（9）：1898 - 1906.
[20] 朱慶瑩， 胡偉艷， 趙志尚. 耕地多功能權(quán)衡與協(xié)同時空格局的動態(tài)分析——以湖北省為例[J] . 經(jīng)濟地理，2018，38（7）：143 - 153.
[21] 方瑩， 王靜， 孔雪松， 等. 耕地利用多功能權(quán)衡關(guān)系測度與分區(qū)優(yōu)化——以河南省為例[J] . 中國土地科學(xué)，2018，32（11）：57 - 64.
[22] 宋爽， 王帥， 傅伯杰， 等. 社會—生態(tài)系統(tǒng)適應(yīng)性治理研究進展與展望[J] . 地理學(xué)報，2019，74（11）：2401 - 2410.
[23] 史小祺， 李陽兵. 貴州省近40年耕地功能轉(zhuǎn)型評價及演變差異分析[J] . 中國巖溶，2018，37（5）：722 - 732.
[24] 史洋洋， 呂曉， 黃賢金， 等. 江蘇沿海地區(qū)耕地利用轉(zhuǎn)型及其生態(tài)系統(tǒng)服務(wù)價值變化響應(yīng)[J] . 自然資源學(xué)報，2017，32（6）：961 - 976.
[25] 包智明， 曾文強. 生計轉(zhuǎn)型與生態(tài)環(huán)境變遷——基于云南省Y村的個案研究[J] . 云南社會科學(xué)，2021（2）：158 - 164.
[26] 郭陽， 徐志剛. 耕地流轉(zhuǎn)市場發(fā)育、資源稟賦與農(nóng)地規(guī)模經(jīng)營發(fā)展[J] . 中國農(nóng)村經(jīng)濟，2021（6）：60 - 75.
[27] HANSEN B E. Threshold effects in non-dynamic panels： estimation， testing and inference[J] . Journal of Econometrics， 1999， 93（2）： 345 - 368.
[28] 苗苗， 迮寒露， 吳永華. 浙江省耕地利用轉(zhuǎn)型時空演變特征分析[J] . 浙江農(nóng)業(yè)學(xué)報，2021，33（4）：753 - 760.
[29] 黃祖輝， 朋文歡， 米松華. 有序推進農(nóng)業(yè)轉(zhuǎn)移人口市民化：浙江實踐與思路對策[J] . 浙江經(jīng)濟，2015（14）：40 - 43.
[30] 陳莎， 楊潤佳， 李冠. 浙江省耕地多功能時空格局變化、影響因素及空間溢出效應(yīng)[J] . 農(nóng)業(yè)工程學(xué)報，2022，38（16）：21 - 32.
[31] 朱從謀， 李武艷， 杜瑩瑩， 等. 浙江省耕地多功能價值時空變化與權(quán)衡—協(xié)同關(guān)系[J] . 農(nóng)業(yè)工程學(xué)報，2020，36（14）：263 - 272.
[32] 宋小青， 李心怡. 區(qū)域耕地利用功能轉(zhuǎn)型的理論解釋與實證[J] . 地理學(xué)報，2019，74（5）：992 - 1010.
[33] 楊雪， 談明洪. 近年來北京市耕地多功能演變及其關(guān)聯(lián)性[J] . 自然資源學(xué)報，2014，29（5）：733 - 743.
[34] 王雨濛， 楊志海. 勞動力轉(zhuǎn)移、資本投入與耕地集約利用[J] . 江西財經(jīng)大學(xué)學(xué)報，2013（5）：86 - 93.
[35] GRAY C L， BILSBORROW R E. Consequences of outmigration for land use in rural Ecuador[J] . Land Use Policy， 2014， 36： 182 - 191.
[36] 杜國明， 薛濡壕， 于鳳榮. 耕地集約利用轉(zhuǎn)型的理論解析[J] . 資源科學(xué)，2022，44（3）：425 - 435.
[37] 李曉方. 發(fā)達國家勞動力轉(zhuǎn)移與土地流轉(zhuǎn)模式對我國的啟示[J] . 農(nóng)業(yè)經(jīng)濟，2022（6）：107 - 109.
[38] FAGERHOLM N， OTEROS-ROZAS E， RAYMOND C M， et al. Assessing linkages between ecosystem services， landuse and well-being in an agroforestry landscape using public participation GIS[J] . Applied Geography， 2016， 74： 30 - 46.
[39] 張一達， 劉學(xué)錄， 任君， 等. 基于耕地多功能權(quán)衡與協(xié)同分析的耕地利用轉(zhuǎn)型研究——以北京市為例[J] . 中國農(nóng)業(yè)資源與區(qū)劃，2020，41（6）：25 - 33.
[40] 向敬偉， 李江風(fēng)， 曾杰. 鄂西貧困縣耕地利用轉(zhuǎn)型空間分異及其影響因素[J] . 農(nóng)業(yè)工程學(xué)報，2016，32（1）：272 - 279.
[41] 谷家川， 查良松. 皖江城市帶農(nóng)作物碳儲量動態(tài)變化研究[J] . 長江流域資源與環(huán)境，2012，21（12）：1507 - 1513.
[42] 陳麗， 郝晉珉， 王峰， 等. 基于碳循環(huán)的黃淮海平原耕地固碳功能研究[J] . 資源科學(xué)，2016，38（6）：1039 - 1053.
[43] 謝光輝， 韓東倩， 王曉玉， 等. 中國禾谷類大田作物收獲指數(shù)和秸稈系數(shù)[J] . 中國農(nóng)業(yè)大學(xué)學(xué)報，2011，16（1）：1 - 8.
[44] 田云， 張俊飚， 李波. 中國農(nóng)業(yè)碳排放研究：測算、時空比較及脫鉤效應(yīng)[J] . 資源科學(xué)，2012，34（11）：2097 -2105.
[45] 賀亞亞， 田云， 張俊飚. 湖北省農(nóng)業(yè)碳排放時空比較及驅(qū)動因素分析[J] . 華中農(nóng)業(yè)大學(xué)學(xué)報（社會科學(xué)版）， 2013（5）：79 - 85.
[46] 閔繼勝， 胡浩. 中國農(nóng)業(yè)生產(chǎn)溫室氣體排放量的測算[J] . 中國人口·資源與環(huán)境，2012，22（7）：21 - 27.
[47] 杜勇， 稅偉， 孫曉瑞， 等. 海灣型城市生態(tài)系統(tǒng)服務(wù)權(quán)衡的情景模擬——以福建省泉州市為例[J] . 應(yīng)用生態(tài)學(xué)報，2019，30（12）：4293 - 4302.
[48] 侯孟陽， 姚順波. 中國農(nóng)村勞動力轉(zhuǎn)移對農(nóng)業(yè)生態(tài)效率影響的空間溢出效應(yīng)與門檻特征[J] . 資源科學(xué)，2018，40（12）：2475 - 2486.
Impact of Resource Endowment and Rural Labor Transfer on Farmland Use Multifunction Transition： An Analysis Based on Threshold Regression
CHEN Sha1， YANG Runjia2， YE Yanmei2， XUE Jibing1
（1. School of Public Administration， Zhejiang University of Finance & Economics， Hangzhou 310018， China； 2. School of Public Affairs， Zhejiang University， Hangzhou 310058， China）
Abstract： The purpose of this study is to explore the influencing mechanism of farmland endowment and rural labor transfer on farmland use multifunction transition in Zhejiang Province. The panel threshold model was used to demonstrate the threshold effects of farmland endowment and rural labor transfer on farmland use multifunction transition and other quantitative relationships among selected variables. The results show that： 1）farmland endowment and rural labor transfer pose threshold effects on both farmland multifunction comprehensive level and functional pattern. In the long run， rural labor transfer tends to be more favorable for the promotion of multifunction comprehensive level， and it poses periodic influence of “positive to negative and slightly negative” on functional pattern. 2）Taking farmland endowment into account， the impact of labor transfer on comprehensive multifunction level and functional pattern changes from “negative” to “positive”， and the positive impact becomes increasingly prominent with the increase of farmland endowment. In conclusion， it is necessary to grasp the developing trend of labor transfer. The improvement in agricultural technology and a series of supporting policy is needed for adjusting labor force structure and agricultural factor input so as to promote the spatial optimization and functional upgrading of agriculture. Meanwhile， the advantages of farmland resource endowment should be exploited， and the inherent coupling and stage difference between farmland resource endowment and socioeconomic development should be investigated. In the process of labor transfer， the government needs to strengthen the comprehensive land consolidation and ecological protection， and to develop new forms of agricultural industry to realize the optimization of farmland use pattern and multifunction transition.
Key words： multifunction transition； panel threshold regression； sustainable farmland use； functional pattern
（本文責(zé)編：張冰松）
