廣西石漠化地區(qū)耕地“非糧化”特征刻畫(huà)及歸因識(shí)別

韋燕飛，魏 鋒，童新華

（1.南寧師范大學(xué) 自然資源與測(cè)繪學(xué)院，廣西 南寧 530100；2.南寧師范大學(xué) 地理科學(xué)與規(guī)劃學(xué)院，廣西 南寧 530100）

0 引言

民為國(guó)基，谷為民命。糧食事關(guān)國(guó)運(yùn)民生，糧食安全是國(guó)家安全的重要基礎(chǔ)。由于糧食生產(chǎn)空間與城鄉(xiāng)發(fā)展空間矛盾日益加深，保障國(guó)家糧食安全的耕地空間面臨被壓縮的現(xiàn)實(shí)困境。數(shù)據(jù)表明，中國(guó)耕地流失占世界耕地流失總量的1?4，糧食生產(chǎn)空間生產(chǎn)潛力總量自改革開(kāi)放以來(lái)約下降了1 ×107t［1］。然而，除了耕地面積的顯性流失之外，工商資本下鄉(xiāng)逐利、農(nóng)業(yè)生產(chǎn)專(zhuān)業(yè)化格局的形成［2］等因素推動(dòng)耕地非糧化帶來(lái)的糧食生產(chǎn)力隱形損失則更為突出?！胺羌Z化”體現(xiàn)為由糧食作物生產(chǎn)調(diào)整為經(jīng)濟(jì)作物種植，是基于農(nóng)業(yè)內(nèi)部比較收益而主動(dòng)進(jìn)行的生產(chǎn)結(jié)構(gòu)調(diào)整。在理想情況下，耕地“非糧化”是追求經(jīng)濟(jì)發(fā)展與耕地保護(hù)相統(tǒng)一的理性選擇，而過(guò)度“非糧化”則是資本逐利性?xún)?nèi)生推動(dòng)的結(jié)果。為此，科學(xué)認(rèn)知“非糧化”的正負(fù)外部性，深入剖析“非糧化”的形成機(jī)理對(duì)于新時(shí)期保證國(guó)家糧食安全具有重要的現(xiàn)實(shí)意義。

目前，國(guó)內(nèi)學(xué)者對(duì)耕地“非糧化”的研究主要集中在以下幾方面：“非糧化”形成機(jī)制研究，從農(nóng)地流轉(zhuǎn)［3，4］、種糧效率低下［5］、成本收益［6－8］、工商資本下鄉(xiāng)［9］、農(nóng)戶(hù)分化［10，11］、政府管理與制度［12］、農(nóng)業(yè)經(jīng)營(yíng)主體的資源稟賦和競(jìng)爭(zhēng)能力［13］等不同維度探析推動(dòng)“非糧化”的形成及過(guò)度“非糧化”的演進(jìn)；耕地“非糧化”負(fù)外部效益的研究，主要體現(xiàn)為壓縮種糧面積［14］、糧食產(chǎn)量減少［15］、用途變化導(dǎo)致的土壤層破壞［16，17］、耕地 利 用 呈 現(xiàn) 出“邊 際 化”等 不 良 趨 勢(shì)［18］；研究尺度的選擇上，涵蓋全國(guó)［19］、省級(jí)單元［20］、鄉(xiāng)鎮(zhèn)級(jí)單元［21］、糧食生產(chǎn)典型區(qū)［22－24］、家庭農(nóng)場(chǎng)［25］。

盡管現(xiàn)有研究已涉及“非糧化”的多個(gè)維度，但在研究尺度的選擇上，多集中于全國(guó)、省域尺度，而對(duì)省域內(nèi)部縣級(jí)之間的空間差異研究則相對(duì)較少，對(duì)微觀(guān)層面的分析更有待深化。因此，本文以廣西巖溶地區(qū)石漠化綜合治理區(qū)域45 個(gè)縣作為研究單元，基于統(tǒng)計(jì)數(shù)據(jù)對(duì)該區(qū)域“非糧化”水平及其類(lèi)型的時(shí)空演化特征進(jìn)行分析，在此基礎(chǔ)上，采用空間計(jì)量模型探究影響耕地“非糧化”的影響因子，從微觀(guān)層面為國(guó)家糧食安全和耕地保護(hù)政策調(diào)整提供科學(xué)參考。

1 研究區(qū)概況、研究方法和數(shù)據(jù)來(lái)源

1.1 研究區(qū)概況

根據(jù)《南方丘陵山地帶生態(tài)保護(hù)和修復(fù)重大工程建設(shè)規(guī)劃（2021—2035 年）》劃定的巖溶地區(qū)石漠化綜合治理區(qū)域，目前區(qū)域內(nèi)現(xiàn)有水土流失面積約2.74× 106hm2，石漠化土地面積1.8 × 106hm2，重度和極重度石漠化土地占45%，涉及湖南、廣西、貴州省共計(jì)69 個(gè)縣，其中廣西段包括南寧、柳州市等9個(gè)地級(jí)市45 個(gè)縣。廣西段屬于中亞熱帶季風(fēng)濕潤(rùn)氣候區(qū)，雨量充沛、氣候溫和，但由于石漠化嚴(yán)重，導(dǎo)致生態(tài)系統(tǒng)極其脆弱，水土流失敏感度高。截至2020年，廣西段石漠化綜合治理區(qū)域的主要糧食作物為稻谷和玉米，主要經(jīng)濟(jì)作物為油料、甘蔗、木薯和蔬菜，農(nóng)作物總播種面積2.463 × 106hm2，糧食種植面積為1.173×106hm2。

圖1 廣西石漠化地區(qū)地理位置Figure 1 Geographical location of rocky desertification area in Guangxi

1.2 研究方法

耕地“非糧化”的概念界定及其測(cè)算?？茖W(xué)界定耕地“非糧化”概念是認(rèn)知并管控耕地“非糧化”現(xiàn)象的重要前提。目前國(guó)內(nèi)外學(xué)界普遍較為認(rèn)同的耕地“非糧化”概念具有以下兩種內(nèi)涵：從狹義上說(shuō)，耕地的“非糧化”是指在耕地上種植糧食作物以外經(jīng)濟(jì)作物的行為，如蔬菜、水果、花卉等；從廣義上說(shuō)，耕地的“非糧化”則是指在耕地從事一切與糧食作物無(wú)關(guān)的、但仍屬于農(nóng)業(yè)生產(chǎn)范疇的行為［11］。在非糧化面積和非糧化率的測(cè)算方面，有學(xué)者認(rèn)為可通過(guò)“糧作比”［26］、“非糧食種植流轉(zhuǎn)耕地面積占流轉(zhuǎn)農(nóng)地面積比例”［27］等指標(biāo)來(lái)描述地區(qū)“非糧化”水平。本文依據(jù)《關(guān)于防止耕地非糧化穩(wěn)定糧食生產(chǎn)的意見(jiàn)》，同時(shí)為了保證數(shù)據(jù)的可得性和廣西傳統(tǒng)糧食種植習(xí)慣，避免蔬菜生產(chǎn)有耕地和設(shè)施農(nóng)用地、水果生產(chǎn)有耕地和園地等多個(gè)交叉來(lái)源此類(lèi)問(wèn)題，僅將稻谷、小麥、玉米列入糧食作物計(jì)算范圍，將除這三類(lèi)的耕地種植行為統(tǒng)一視為非糧化的具體表現(xiàn)。計(jì)算公式為：

式中：Ang為非糧化面積；Ac為耕地面積；Ag為糧食播種面積；I 為復(fù)種指數(shù)；Rng為非糧化率。

空間探索性研究方法（ESDA）。ESDA 是探索研究區(qū)域各項(xiàng)指標(biāo)的空間分布格局和相互作用機(jī)理的一種空間分析方法，主要用于檢驗(yàn)具有空間位置的某要素的觀(guān)測(cè)值是否顯著地與其相鄰空間點(diǎn)上的觀(guān)測(cè)值相關(guān)聯(lián)［29］。包括全局空間自相關(guān)和局部空間自相關(guān)，通常采用Moran′s I 指數(shù)表示空間要素的自相關(guān)程度。Moran′s I 指數(shù)的值域?yàn)椋郏?，1］，小于0說(shuō)明負(fù)相關(guān)，等于0 說(shuō)明不相關(guān)，大于0 說(shuō)明正相關(guān)。通常用Z 值檢驗(yàn)自相關(guān)顯著性，當(dāng)Z ＞1.96或Z ＜－1.96 時(shí)，變量的空間自相關(guān)顯著，否則不顯著。

全局Moran′s I 指數(shù)的表達(dá)式為：

空間滯后模型（Spatial lag model，SLM）。該模型主要探究的是各變量在某一地區(qū)發(fā)揮作用時(shí)是否存在空間溢出效應(yīng)，也被稱(chēng)為空間自回歸模型。空間滯后模型考慮周邊區(qū)域的被解釋變量對(duì)研究區(qū)的被解釋變量的影響，以空間滯后項(xiàng)的形式加入計(jì)量方程，形成的空間計(jì)量模型。SLM 模型表達(dá)式如下［30］：

式中：ε～N（0，σ2In）；y 為被解釋變量矩陣；X為解釋變量矩陣；W 為鄰接權(quán)重矩陣；ρ為空間效應(yīng)系數(shù)；參數(shù)β反映自變量對(duì)因變量的影響程度；Wy為空間滯后因變量，屬于內(nèi)生變量，反映空間距離對(duì)區(qū)域行為的作用。

空間誤差模型（Spatial error model，SEM）。該模型主要解決的是隨機(jī)誤差之間的空間自相關(guān)問(wèn)題，盡可能減小測(cè)量周邊區(qū)域時(shí)產(chǎn)生的不確定性因素對(duì)本地相關(guān)活動(dòng)產(chǎn)生的影響，也被稱(chēng)為空間自相關(guān)模型?？臻g誤差模型考慮周邊區(qū)域的解釋變量對(duì)研究區(qū)的解釋變量的影響，以空間誤差項(xiàng)的形式加入計(jì)量方程，形成的空間計(jì)量模型。SEM 模型表達(dá)式如下［30］：

式中：ε為隨機(jī)干擾向量；λ為空間誤差相關(guān)系數(shù)，度量了樣本中觀(guān)測(cè)值相互的空間依賴(lài)作用；Wε為空間滯后誤差項(xiàng)；μ為服從正態(tài)分布的隨機(jī)誤差項(xiàng)向量；其他同公式（6）。

1.3 數(shù)據(jù)來(lái)源

以廣西巖溶地區(qū)石漠化綜合治理區(qū)域45 個(gè)縣為研究單元，土地利用數(shù)據(jù)選用中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心Landsat 8?TM 數(shù)據(jù)（空間分辨率為15m×15m），2020 年農(nóng)作物總播種面積、糧食作物播種面積、地區(qū)生產(chǎn)總值等社會(huì)經(jīng)濟(jì)數(shù)據(jù)來(lái)自《廣西統(tǒng)計(jì)年鑒》《中國(guó)縣域統(tǒng)計(jì)年鑒》和各地級(jí)市政府門(mén)戶(hù)網(wǎng)站資料。通過(guò)查閱年鑒直接獲取或間接計(jì)算得到指標(biāo)的所需數(shù)值，部分缺失數(shù)據(jù)通過(guò)SPSS 26 軟件中多重插補(bǔ)辦法予以補(bǔ)充。

2 結(jié)果及分析

2.1 廣西耕地“非糧化”的總體特征

本文基于Arcgis 10.2 軟件中的自然斷點(diǎn)分級(jí)法，將廣西石漠化地區(qū)2020 年的非糧化面積和非糧化率的空間分布情況進(jìn)行可視化分析。

從耕地非糧化面積視角來(lái)看，廣西石漠化地區(qū)2020 年“非糧化”面積為1.324×106hm2，總體上形成以南寧、來(lái)賓和柳州市隸屬縣份為核心的高值地區(qū)并以百色市隸屬縣份為代表的低值地帶（圖2）。整體上該區(qū)域的非糧化面積均值為5.76hm2，興賓區(qū)和武鳴縣的非糧化面積分別達(dá)到1.335× 105hm2和1.053 × 105hm2，為該區(qū)域的最高值，而合山市和凌云縣的非糧化面積僅為0.65hm2和0.69hm2，為該區(qū)域的最低值。其中，興賓區(qū)和武鳴縣的糧食產(chǎn)量下降趨勢(shì)顯著，非糧化形勢(shì)最為嚴(yán)峻。興賓區(qū)為典型的巖溶溶蝕地貌，土壤薄而貧瘠，地表水流失嚴(yán)重，屬干旱地區(qū)，糧作條件局限性大。武鳴縣丘陵占土地總面積63.5%，地塊規(guī)模小而分散，不利于大規(guī)模機(jī)械化耕作，同時(shí)武鳴區(qū)隸屬?gòu)V西省會(huì)南寧市，受土地流轉(zhuǎn)成本相對(duì)較高、地方政府財(cái)政壓力大等外部因素的驅(qū)動(dòng)，傾向于種植經(jīng)濟(jì)效益更大的非糧作物。合山市種植業(yè)以水稻等糧食作物為主，同時(shí)由于該區(qū)域行政面積較小，因此合山市非糧化現(xiàn)象并不突出。凌云縣地處云貴高原的延伸部分，山地面積占全縣總面積的93. 32%，耕地僅占全區(qū)的5.47%，這導(dǎo)致凌云縣耕作面積有限，間接降低了該區(qū)域的耕地非糧化面積。

圖2 廣西石漠化地區(qū)非糧化面積空間分布Figure 2 Spatial distribution of non－ grainized areas in rocky desertification areas in Guangxi

從耕地非糧化率視角出發(fā)，廣西石漠化地區(qū)2020年平均非糧化率為37.25%，總體上形成以中段東西走向?yàn)楹诵牡母咧档貐^(qū)，并以左段和中間南北走向?yàn)楹诵牡牡椭档貛В▓D3）。其中，蒙山縣和巴馬瑤族自治縣的非糧化率分別達(dá)到62.33%和58.09%，為廣西石漠化區(qū)域的最高值，而上林縣和融水苗族自治縣分別為9.16%和10.91%，為該區(qū)域的最低值。蒙山縣在經(jīng)歷早期地殼運(yùn)動(dòng)的背景下，地質(zhì)構(gòu)造相對(duì)復(fù)雜，開(kāi)展大規(guī)模的糧食作物生產(chǎn)難度大，因此農(nóng)業(yè)以水果等經(jīng)濟(jì)作物為主要代表，農(nóng)民人均種植水果收入持續(xù)增長(zhǎng)，進(jìn)一步鼓勵(lì)農(nóng)民“棄糧種果”的經(jīng)濟(jì)行為。巴馬瑤族自治縣素有“八山一水一分田”之稱(chēng)，石山地占30%，山多地少的地質(zhì)狀況使得土地尤為珍貴，因此農(nóng)作物選擇更傾向收益水平更高的經(jīng)濟(jì)作物。上林縣于2010 年列入滇桂黔石漠化片區(qū)縣，后經(jīng)石漠化綜合治理工程開(kāi)展，2020年經(jīng)濟(jì)作物播種面積1.83 × 104hm2，同比下降7.07%，非糧化趨勢(shì)遏制成效顯著。融水苗族自治縣土地以山地為主，故有“九山半水半分田”之說(shuō)，山地占土地面積的85.48%。全縣耕地范圍較小，農(nóng)業(yè)以糧食種植為主，故非糧化率也相對(duì)較低。

圖3 廣西石漠化地區(qū)非糧化率空間分布Figure 3 Spatial distribution of non－ grainization rates in rocky desertification areas in Guangxi

2.2 廣西耕地非糧化的空間格局

將廣西石漠化地區(qū)45 個(gè)縣區(qū)非糧化面積和非糧化率作為觀(guān)測(cè)變量，采用Global Moran′s I 指數(shù)判定不同尺度非糧化的空間特征。從圖4 可見(jiàn)：縣域?qū)用娴母胤羌Z化面積和非糧化率的Moran′s I ＞0，Moran′s I 分別為0.1 978 61 和0.1 159 59，且全部通過(guò)了P ＜0.001 的顯著性檢驗(yàn)，表明在空間上存在一定的相關(guān)性，但相關(guān)性不強(qiáng)。

圖4 廣西石漠化地區(qū)耕地非糧化面積與非糧化率Global Moran′s I 散點(diǎn)圖Figure 4 Global Moran′s I scatter plot of cultivated land non－ grain area and non－ grain rate in rocky desertification area of Guangxi

在Global Moran′s I 的基礎(chǔ)上，進(jìn)一步引入Local Moran′s I 分析各縣區(qū)的集聚模式和相互間的空間狀態(tài)（圖5）。

圖5 廣西石漠化地區(qū)縣級(jí)單元耕地非糧化率與非糧化面積空間集聚Figure 5 Spatial agglomeration of non－ grain conversion rate and non－ grain area of cultivated land in 111 county－ level units in Guangxi

由圖5a可知，非糧化面積在空間上形成了以百色—河池市為核心的低—低集聚區(qū)，涉及7 個(gè)縣級(jí)單元，這個(gè)分區(qū)內(nèi)的縣級(jí)單元具有明顯的空間相關(guān)性，即相鄰區(qū)域的非糧化面積均相對(duì)較小，這類(lèi)區(qū)域多為糧食產(chǎn)區(qū)，非糧化面積相對(duì)較少；以南寧—來(lái)賓—柳州市為核心的高—高集聚區(qū)，涉及5 個(gè)縣級(jí)單元，相鄰區(qū)縣間非糧化面積相互關(guān)系強(qiáng)，主要?dú)w結(jié)于地區(qū)石漠化程度、經(jīng)濟(jì)條件等因素的推動(dòng)；以上林—合山市為代表低—高集聚區(qū)，即自身非糧化面積相對(duì)較小而相鄰縣域非糧化面積相對(duì)較大，主要是由于自身耕地少、區(qū)位特殊，導(dǎo)致糧作面積碎片化程度高；其余縣域耕地非糧化面積的空間相關(guān)性均不顯著。

由圖5b可知，非糧化率的局部自相關(guān)特征與耕地非糧化差別較大，主要是以高—高集聚區(qū)為主，橫跨來(lái)賓、桂林和梧州市，共計(jì)5 個(gè)縣級(jí)單元，這些單元的非糧化率較高且與臨近單元影響密切，主要是因?yàn)槭潭容^高，導(dǎo)致耕地碎片化程度高，無(wú)法開(kāi)展規(guī)模化生產(chǎn)，糧作收益難以覆蓋機(jī)械化經(jīng)營(yíng)成本而促使非糧化水平的提升；以樂(lè)業(yè)、凌云和鳳山縣為代表的低—低集聚區(qū)，這與非糧化面積的低—低集聚區(qū)空間分布大體一致，主要是因?yàn)楫a(chǎn)糧區(qū)域的非糧化管控更為嚴(yán)格，更加側(cè)重耕作的社會(huì)屬性而非經(jīng)濟(jì)屬性；以巴馬瑤族自治縣和大新縣為代表的高—低集聚區(qū)，體現(xiàn)為非糧化率顯著高于周邊縣域單元；以象州縣為代表的低—高聚集區(qū)，體現(xiàn)為自身非糧化率低且與鄰近單元的特征差異顯著；其余縣域耕地非糧化率的空間相關(guān)性均不顯著。

2.3 驅(qū)動(dòng)因子識(shí)別及地域分異

驅(qū)動(dòng)因子的初步選?。夯谇叭说南嚓P(guān)研究和數(shù)據(jù)的可獲取性，從自然資源稟賦、經(jīng)濟(jì)發(fā)展水平、要素投入水平和居民生活水平4 個(gè)方面選取17 個(gè)具有代表性的指標(biāo)作為解釋變量，將耕地非糧化面積和耕地非糧化率作為被解釋變量，構(gòu)建符合廣西石漠化地區(qū)現(xiàn)實(shí)情況的耕地非糧化影響因素指標(biāo)體系（表1）。

表1 廣西耕地非糧化影響因素指標(biāo)體系Table 1 Index system of influencing factors of cultivated land non－ grain in Guangxi

驅(qū)動(dòng)因子的共線(xiàn)性檢驗(yàn)：對(duì)預(yù)選因子進(jìn)行共線(xiàn)性檢驗(yàn)，基于SPSS軟件通過(guò)Z－score標(biāo)準(zhǔn)化方法對(duì)影響因子進(jìn)行標(biāo)準(zhǔn)化處理，接著采用Pearson 相關(guān)性分析進(jìn)行多重共線(xiàn)性診斷，保留VIF 值小于10，容差大于0.1 的影響因子，結(jié)果表明農(nóng)作物總播種面積、糧食作物播種面積、耕地面積和設(shè)施農(nóng)業(yè)占地面積的VIF 值大于10，容差小于0.1，存在共線(xiàn)性，同時(shí)GDP、農(nóng)村居民恩格爾系數(shù)、單位面積第一產(chǎn)業(yè)產(chǎn)值和城鄉(xiāng)居民收入比未通過(guò)顯著性檢驗(yàn)（P ＞0.05），因此將這8 個(gè)影響因子進(jìn)行剔除。廣西石漠化地區(qū)耕地非糧化與影響因素的共線(xiàn)性診斷結(jié)果，如表2 所示。

表2 影響因素共線(xiàn)性診斷分析結(jié)果Table 2 Influencing factors collinearity diagnostic analysis results

模型選擇及估計(jì)結(jié)果：為了選擇有統(tǒng)計(jì)學(xué)意義顯著性的空間計(jì)量模型，采用經(jīng)典回歸OLS 模型對(duì)計(jì)算得到的數(shù)據(jù)進(jìn)行檢驗(yàn)，用OLS 回歸后的殘差進(jìn)行LMLAG LMERR檢驗(yàn)。根據(jù)拉格朗日乘子檢驗(yàn)原理，如果LM－Error 和LM－Lag 這兩個(gè)統(tǒng)計(jì)量 均不顯著，則選擇OLS模型作為最終模型；若只有1 個(gè)統(tǒng)計(jì)量顯著，那么LM－Error 統(tǒng)計(jì)量顯著則指向空間誤差模型（SEM），而LM－Lag 統(tǒng)計(jì)量顯著則指向空間滯后模型（SLM）。結(jié)果表明，LM－ Error 的顯著性?xún)?yōu)于LM－ Lag，故選擇空間誤差模型作為最佳模型，并運(yùn)用GeoDA1.10 對(duì)截面數(shù)據(jù)進(jìn)行空間回歸檢驗(yàn)。

從統(tǒng)計(jì)學(xué)顯著性檢驗(yàn)結(jié)果看，9 個(gè)指標(biāo)均通過(guò)P ＜0.01 置信水平的顯著性檢驗(yàn)，這表明農(nóng)村居民最低生活保障人數(shù)、第一產(chǎn)業(yè)增加值、城鎮(zhèn)化水平、耕地年退化率、水土流失率、單位耕地農(nóng)業(yè)機(jī)械總動(dòng)力、人均耕地面積、第一產(chǎn)業(yè)對(duì)經(jīng)濟(jì)貢獻(xiàn)值和復(fù)種指數(shù)對(duì)耕地非糧化面積和耕地非糧化率均有顯著的驅(qū)動(dòng)作用，同時(shí)二者的R2接近1，說(shuō)明模型的擬合程度較高，這驗(yàn)證了廣西石漠化地區(qū)耕地非糧化影響因素框架的正確性（表3）。

表3 空間滯后模型SEM選擇檢驗(yàn)及估計(jì)結(jié)果Table 3 SEM selection test and estimation results

2.4 耕地非糧化面積驅(qū)動(dòng)力分析

從非糧化面積的回歸系數(shù)來(lái)看，第一產(chǎn)業(yè)增加值（35.5332）、城鎮(zhèn)化水平（24.0685）、農(nóng)村居民最低生活保障人數(shù)（18.0384）和人均耕地面積（15.9446）顯著高于其他指標(biāo)，表明自然資源稟賦、經(jīng)濟(jì)發(fā)展水平和居民生活水平對(duì)耕地“非糧化”面積的正向影響最大，而單位耕地機(jī)械總動(dòng)力（－3.8331）、水土流失率（－3.4592）和耕地年退化率（－ 2.3078）指標(biāo)數(shù)值為負(fù)，表明對(duì)耕地非糧化面積的負(fù)面影響主要集中在自然資源稟賦和要素投入水平。

從自然資源稟賦看，人均耕地面積和復(fù)種指數(shù)對(duì)耕地“非糧化”具有驅(qū)動(dòng)作用，其中人均耕地面積的驅(qū)動(dòng)屬性最為顯著，表明人均耕地面積每增加1個(gè)單位，廣西石漠化地區(qū)的耕地非糧化面積就會(huì)上升15.9446 個(gè)單位，而耕地年退化率和水土流失率均為負(fù)值，表明該地區(qū)耕地質(zhì)量的下降會(huì)對(duì)耕地“非糧化”的趨勢(shì)產(chǎn)生抑制作用，但耕地年退化率和水土流失率同為負(fù)向指標(biāo)，因此不能簡(jiǎn)單地把能約束“非糧化”現(xiàn)象的指標(biāo)定義為科學(xué)收斂“非糧化”趨勢(shì)的有效途徑。

從經(jīng)濟(jì)發(fā)展水平看，第一產(chǎn)業(yè)增加值和城鎮(zhèn)化水平顯著性最強(qiáng)，二者每增加1 個(gè)單位，“非糧化”面積就會(huì)增加35.5332、24.0685 個(gè)單位。從空間視角出發(fā)，廣西石漠化地區(qū)中南部區(qū)域非糧化面積大，涉及區(qū)縣隸屬于南寧、柳州市，城市定位側(cè)重經(jīng)濟(jì)、政治層面而非主要糧食生產(chǎn)承載地，而城鎮(zhèn)化水平作為城市經(jīng)濟(jì)水平的表征，城市經(jīng)濟(jì)發(fā)展程度越高，則農(nóng)業(yè)比較效益越低，農(nóng)民種糧意愿越低，“非糧化”面積擴(kuò)大趨勢(shì)更加顯著。第一產(chǎn)業(yè)對(duì)經(jīng)濟(jì)的貢獻(xiàn)值同樣對(duì)耕地非糧化面積的擴(kuò)大具有正向作用，這表明經(jīng)濟(jì)利益的驅(qū)動(dòng)同樣會(huì)導(dǎo)致區(qū)域耕地非糧化面積的擴(kuò)大。

從要素投入水平看，單位耕地機(jī)械總動(dòng)力的數(shù)值為負(fù)，對(duì)耕地“非糧化”現(xiàn)象同樣具備一定程度的抑制作用。由于耕作機(jī)械化作為現(xiàn)代農(nóng)業(yè)科技要素投入的主要表現(xiàn)形式，有利于糧食作物的規(guī)?；a(chǎn)，然而在石漠化地區(qū)當(dāng)農(nóng)業(yè)生產(chǎn)過(guò)程中不具備大規(guī)模使用農(nóng)業(yè)機(jī)械的條件時(shí)，糧食作物播種面積就會(huì)呈現(xiàn)碎片化，耕地糧作功能被削弱。從務(wù)農(nóng)人員的角度出發(fā)，選擇農(nóng)業(yè)機(jī)械需求較低的經(jīng)濟(jì)作物是更為理性的行為，并且經(jīng)濟(jì)收益的驅(qū)動(dòng)會(huì)放大這種農(nóng)作行為。

從居民生活水平看，農(nóng)村居民最低生活保障人數(shù)的增多對(duì)耕地非糧化面積的擴(kuò)大具有顯著的驅(qū)動(dòng)效應(yīng)，表現(xiàn)為農(nóng)村居民最低生活保障人數(shù)每增加1個(gè)單位，耕地非糧化面積就會(huì)增長(zhǎng)18.0384 個(gè)單位。從本質(zhì)上講，這是經(jīng)濟(jì)性的耕作行為，石漠化地區(qū)是貧困高發(fā)區(qū)，地區(qū)的居民生活水平越低，在一定程度上會(huì)更傾向于選擇經(jīng)濟(jì)作物來(lái)提升收益以改善生活水平。

2.5 耕地非糧化率驅(qū)動(dòng)力分析

從非糧化率的回歸系數(shù)來(lái)看，第一產(chǎn)業(yè)增加值（18.1466）、農(nóng)村居民最低生活保障人數(shù)（10.1484）、復(fù)種指數(shù)（6. 0254）和單位面積第一產(chǎn)業(yè)產(chǎn)值（8.0254）對(duì)非糧化率增長(zhǎng)的驅(qū)動(dòng)效應(yīng)最顯著，表明經(jīng)濟(jì)發(fā)展水平、自然資源稟賦和居民生活水平對(duì)耕地非糧化率的正向影響最大，同時(shí)耕地年退化率（6.0047）和水土流失率（5.1170）對(duì)非糧化率的提高也具備驅(qū)動(dòng)效應(yīng)，而單位耕地機(jī)械總動(dòng)力（－4.0344）則有利于非糧化率的減低，表明要素投入水平的提高在一定程度上會(huì)導(dǎo)致耕地非糧化率的下降。其余指標(biāo)雖顯著性不強(qiáng)，但仍具備一定的參考性，也能一定程度反映與耕地“非糧化”水平的相關(guān)性。

從自然資源稟賦看，水土流失率、耕地年退化率、人均耕地面積和復(fù)種指數(shù)對(duì)耕地非糧化率的提高均具有促進(jìn)作用，主要是由于水土流失和耕地退化均會(huì)導(dǎo)致人均可利用耕地大幅減少，在耕地要素不足的情況下，少量的糧食作物產(chǎn)量帶來(lái)的收入難以維系務(wù)農(nóng)家庭的開(kāi)支，因此驅(qū)使務(wù)農(nóng)人員選擇收益水平更高的經(jīng)濟(jì)作物，同時(shí)復(fù)種程度越高，則非糧化率會(huì)也會(huì)同步增大。

從經(jīng)濟(jì)發(fā)展水平看，第一產(chǎn)業(yè)增加值和單位面積第一產(chǎn)業(yè)產(chǎn)值對(duì)耕地非糧化率的提高具有顯著的驅(qū)動(dòng)作用。由于耕地非糧化作用的直接結(jié)果是經(jīng)濟(jì)收益的提升，因此從系統(tǒng)的角度來(lái)看，這兩者與耕地非糧化率在一定程度上會(huì)產(chǎn)生相互作用并且同向發(fā)展，即非糧化率越高，第一產(chǎn)業(yè)產(chǎn)值提升的越快，反之亦然。從保障國(guó)家糧食安全的角度出發(fā)，單純追求經(jīng)濟(jì)利益的耕作行為與我國(guó)的糧食安全政策相駁，與國(guó)家總體安全的要求不符。

從要素投入水平看，單位耕地農(nóng)業(yè)機(jī)械總動(dòng)力的增長(zhǎng)對(duì)降低非糧化率具有一定的成效。由于大部分糧食作物對(duì)農(nóng)業(yè)機(jī)械的依賴(lài)性較高，在農(nóng)業(yè)機(jī)械投入不足的背景下，務(wù)農(nóng)人員更青睞對(duì)機(jī)械要素依賴(lài)程度較低的經(jīng)濟(jì)作物，因此適度加大農(nóng)業(yè)機(jī)械投入量不僅會(huì)提升務(wù)農(nóng)人員耕種經(jīng)濟(jì)作物的積極性，降低務(wù)農(nóng)的人工成本，同時(shí)也有利于提高務(wù)農(nóng)效率，科學(xué)管理農(nóng)作物生長(zhǎng)環(huán)境。

從居民生活水平看，農(nóng)村居民最低生活保障人數(shù)基數(shù)的擴(kuò)大也會(huì)一定程度上導(dǎo)致非糧化率的提高，這與此指標(biāo)導(dǎo)致非糧化面積增大的邏輯大體一致，即為改善生活條件而選擇在耕地上種植收益較高的經(jīng)濟(jì)作物。從這個(gè)角度出發(fā)，國(guó)家應(yīng)該更多的考慮改善基層務(wù)農(nóng)人民的收入水平和結(jié)構(gòu)，在收入得到保障的背景下，也能在一定程度緩解“非糧化”現(xiàn)象的發(fā)展趨勢(shì)。

3 結(jié)論、討論與建議

3.1 結(jié)論

本文從廣西石漠化地區(qū)非糧化面積和非糧化率的特征進(jìn)行描述分析，并通過(guò)空間誤差模型識(shí)別耕地“非糧化”的關(guān)鍵驅(qū)動(dòng)因子。主要結(jié)論如下：①?gòu)V西石漠化地區(qū)非糧化總體特征的刻畫(huà)分為兩方面：非糧化面積方面，廣西石漠化地區(qū)2020 年平均非糧化面積為29418.27hm2，總體上形成以南寧、來(lái)賓和柳州市隸屬縣份為核心的高值地區(qū)，并以百色市隸屬縣份為代表的低值地帶；非糧化率方面，廣西石漠化地區(qū)2020 年平均非糧化率為37.25%，總體上形成以中段（東西走向）為核心的高值地區(qū)并以左段和中間（南北走向）為核心的低值地帶。②以廣西石漠化地區(qū)非糧化面積和非糧化率作為觀(guān)測(cè)變量，通過(guò)全局莫蘭指數(shù)表明縣域?qū)用娴母胤羌Z化面積和非糧化率在空間上存在一定的相關(guān)性，在此基礎(chǔ)上，運(yùn)用局部莫蘭指數(shù)深入分析縣域的集聚模式，發(fā)現(xiàn)非糧化面積與非糧化率的空間自相關(guān)特征差別較大，同時(shí)非糧化面積的集聚效應(yīng)比非糧化率更顯著，非糧化率的空間狀態(tài)相對(duì)離散。③廣西石漠化區(qū)域耕地“非糧化”驅(qū)動(dòng)力方面，從非糧化面積的回歸系數(shù)來(lái)看，第 一 產(chǎn) 業(yè) 增 加 值（35. 5332）、城 鎮(zhèn) 化 水 平（24.0685）、農(nóng)村居民最低生活保障人數(shù)（18.0384）和人均耕地面積（15.9446）對(duì)非糧化面積增長(zhǎng)的驅(qū)動(dòng)效應(yīng)最顯著；從非糧化率的回歸系數(shù)來(lái)看，第一產(chǎn)業(yè)增加值（18.1466）、農(nóng)村居民最低生活保障人數(shù)（10.1484）、復(fù)種指數(shù)（6.0254）和單位面積第一產(chǎn)業(yè)產(chǎn)值（8.0254）對(duì)非糧化率增長(zhǎng)具有顯著的驅(qū)動(dòng)效應(yīng)?？傮w而言，第一產(chǎn)業(yè)增加值、農(nóng)村居民最低生活保障人數(shù)均對(duì)非糧化面積和非糧化率的提高具有顯著影響，國(guó)家應(yīng)持續(xù)關(guān)注石漠化地區(qū)自然資源限制所誘發(fā)的貧困問(wèn)題，平衡經(jīng)濟(jì)發(fā)展—糧食安全的關(guān)系，正視經(jīng)濟(jì)收入欠保障的務(wù)農(nóng)人員對(duì)于耕地非糧化趨勢(shì)的推動(dòng)作用。

3.2 討論

本文是以廣西石漠化地區(qū)作為研究對(duì)象，從縣域尺度展開(kāi)耕地非糧化研究，評(píng)價(jià)單元較小，雖然評(píng)價(jià)結(jié)果精確度更高、差異性更顯著，但由于部分?jǐn)?shù)據(jù)并未統(tǒng)計(jì)到縣級(jí)單元，存在指標(biāo)數(shù)據(jù)缺失的約束性，因此會(huì)影響評(píng)價(jià)指標(biāo)體系的科學(xué)性和完整性。同時(shí)，本文在進(jìn)行驅(qū)動(dòng)力分析時(shí)，把廣西石漠化地區(qū)作為整體，尚未觸及縣級(jí)單元耕地非糧化歸因的地域差異性，接下來(lái)的研究中將進(jìn)一步深入分析各縣級(jí)單元耕地非糧化的具體成因，以為地域性政策的提出提供理論參考。

從實(shí)踐來(lái)看，種植糧食作物比較利益低下、土地流轉(zhuǎn)成本攀升及部分地區(qū)不合理的鄉(xiāng)村振興措施對(duì)耕地非糧化局面的形成具有助推作用。從形式來(lái)看，主要包括耕地撂荒、農(nóng)村發(fā)展經(jīng)濟(jì)的要求和經(jīng)濟(jì)利益驅(qū)使。因此，“非糧化”整治并非簡(jiǎn)單地清理與“執(zhí)法”，應(yīng)注重以下三方面工作的開(kāi)展：第一，協(xié)同好耕地“非糧化”與“藏糧于地”的關(guān)系；第二，既要提高地方政府激勵(lì)糧食生產(chǎn)的內(nèi)生動(dòng)力，也要提高經(jīng)營(yíng)者的糧食生產(chǎn)積極性；第三，客觀(guān)認(rèn)識(shí)不同耕地“非糧化”現(xiàn)象對(duì)糧食安全的影響。

3.3 建議

基于上述結(jié)論，提出以下建議：①提高糧食種植效益。一是通過(guò)打造地方糧食品牌，提升產(chǎn)前研發(fā)和產(chǎn)后精深加工環(huán)節(jié)延伸投入力度以?xún)?yōu)化糧食產(chǎn)業(yè)鏈和供應(yīng)鏈；二是加入對(duì)農(nóng)業(yè)基礎(chǔ)設(shè)施的公共投入力度；三是加強(qiáng)農(nóng)業(yè)科技投入力度。②規(guī)范農(nóng)村土地流轉(zhuǎn)市場(chǎng)。一是引導(dǎo)耕地流轉(zhuǎn)費(fèi)處于合理水平，制定全國(guó)性耕地流轉(zhuǎn)價(jià)格指導(dǎo)標(biāo)準(zhǔn)，避免由于過(guò)高或過(guò)低的耕地流轉(zhuǎn)費(fèi)對(duì)市場(chǎng)產(chǎn)生負(fù)面影響；二是強(qiáng)化耕地流轉(zhuǎn)后用途管制，首先是制度層面的約束并輔以相應(yīng)的監(jiān)管配套措施，其次可建立農(nóng)村土地流轉(zhuǎn)管理信息系統(tǒng)，明確記錄土地流轉(zhuǎn)前后的信息變更。