我國(guó)馬鈴薯區(qū)域格局演變及其影響因素分析

張 爍 羅其友* 馬力陽(yáng)

(1.中國(guó)農(nóng)業(yè)科學(xué)院 農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所，北京 100081；2.西安科技大學(xué) 馬克思主義學(xué)院，西安 710054)

馬鈴薯是我國(guó)第四大糧食作物，對(duì)保障我國(guó)糧食安全和鄉(xiāng)村脫貧攻堅(jiān)具有重要作用[1-2]。近年來(lái)，馬鈴薯產(chǎn)業(yè)快速發(fā)展[3-4]，種植面積不斷擴(kuò)大。由于一些地區(qū)盲目擴(kuò)種、跟風(fēng)種植，導(dǎo)致馬鈴薯市場(chǎng)競(jìng)爭(zhēng)加劇，價(jià)格波動(dòng)加大，部分地區(qū)市場(chǎng)優(yōu)勢(shì)不明顯、種植規(guī)模縮減，馬鈴薯區(qū)域格局正在發(fā)生變化。因此，開(kāi)展馬鈴薯區(qū)域格局演變研究及其影響因素分析，對(duì)優(yōu)化馬鈴薯生產(chǎn)區(qū)域布局、促進(jìn)馬鈴薯產(chǎn)業(yè)持續(xù)發(fā)展具有重要意義。

已有馬鈴薯區(qū)域格局演變的相關(guān)文獻(xiàn)從不同角度研究了馬鈴薯區(qū)域分布特征及變化趨勢(shì)，主要研究?jī)?nèi)容包括區(qū)域比較優(yōu)勢(shì)分析[5-7]、種植區(qū)域集中程度揭示[8-9]、生產(chǎn)重心演變軌跡和種植區(qū)域的擴(kuò)展趨勢(shì)研究等[10-11]。在此基礎(chǔ)上，進(jìn)一步分析區(qū)域分布變化的成因，結(jié)果顯示影響馬鈴薯空間分布的自然因素包括氣候條件和自然災(zāi)害，社會(huì)經(jīng)濟(jì)因素包括市場(chǎng)需求、價(jià)格、生產(chǎn)技術(shù)水平、比較收益和政策等[8,10,12]。主要研究方法包括定性分析[6，9]和定量分析[13-14]，并且對(duì)影響因素進(jìn)行了區(qū)域?qū)Ρ确治鯷12]。然而，已有研究大都從省域尺度分析區(qū)域分布特征及變化趨勢(shì)，缺少縣域尺度的研究。在區(qū)域格局演變的成因方面，影響因素包括了自然、社會(huì)和經(jīng)濟(jì)因素，且研究方法不斷完善，但實(shí)證研究較少，尤其缺少比較不同因素在不同地區(qū)和不同時(shí)間段的重要程度。因此，本研究擬基于全國(guó)馬鈴薯縣級(jí)數(shù)據(jù)，研究馬鈴薯種植格局演變趨勢(shì)，并進(jìn)一步研究其影響因素，以期為馬鈴薯種植合理布局提供參考，為馬鈴薯區(qū)劃研究提供理論依據(jù)。

1 數(shù)據(jù)來(lái)源與研究方法

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

本研究以全國(guó)作為研究的區(qū)域范圍(統(tǒng)計(jì)數(shù)據(jù)未包含港澳臺(tái)地區(qū)。下同。)，使用1990—2016年我國(guó)馬鈴薯縣級(jí)種植面積、產(chǎn)量和單產(chǎn)數(shù)據(jù)研究區(qū)域格局演變趨勢(shì)。同時(shí)，采用2011—2018年作為時(shí)間樣本，選取黑龍江、遼寧、陜西、甘肅、青海、寧夏、新疆、內(nèi)蒙古、河北、山西、重慶、四川、貴州、云南14 個(gè)省(市、自治區(qū))作為地區(qū)樣本，利用馬鈴薯省級(jí)數(shù)據(jù)研究區(qū)域格局演變的影響因素。馬鈴薯縣級(jí)種植面積和脫毒種薯應(yīng)用數(shù)據(jù)源于國(guó)家馬鈴薯產(chǎn)業(yè)技術(shù)體系，省級(jí)單產(chǎn)數(shù)據(jù)源于國(guó)家統(tǒng)計(jì)局，用工數(shù)量、化肥費(fèi)用和種子費(fèi)用數(shù)據(jù)源于《全國(guó)農(nóng)產(chǎn)品成本收益資料匯編》(2012—2019)[15]、全國(guó)氣象監(jiān)測(cè)站降水和氣溫月值數(shù)據(jù)源于國(guó)家氣象局。

1.2 研究方法

1.2.1核密度分析

依據(jù)Arc GIS提供的核密度分析工具[16]，運(yùn)用核密度工具探討馬鈴薯種植區(qū)域的空間分布特征，核密度分析法主要是借助移動(dòng)的柵格單元對(duì)空間范圍內(nèi)點(diǎn)要素的密集程度進(jìn)行估計(jì)[17]，計(jì)算公式如下：

(1)

式中：D(xi，yi)為空間位置(xi，yi)處的核密度值；r為距離衰減閾值；u為與位置(xi，yi)的距離≤r的要素點(diǎn)數(shù)；k函數(shù)則表示空間權(quán)重函數(shù)；d表示當(dāng)前要素點(diǎn)與(xi，yi)2點(diǎn)之間的歐氏距離。

1.2.2全局自相關(guān)

全局空間自相關(guān)是對(duì)地理要素屬性值在整個(gè)區(qū)域的空間特征描述，選取全局莫蘭指數(shù)來(lái)測(cè)度馬鈴薯種植的全局空間自相關(guān)程度，計(jì)算公式如下[18-19]：

(2)

(3)

式中：E(I)、VAR(I)分別表示莫蘭指數(shù)的期望值和方差。全局莫蘭指數(shù)的取值介于[-1，1]，在給定的顯著性水平下，當(dāng)I>0，表明存在空間正相關(guān)，即馬鈴薯種植面積較高(或較低)的地區(qū)在空間上顯著聚集；當(dāng)I<0，表示空間負(fù)相關(guān)，即該區(qū)域與周邊地區(qū)馬鈴薯種植面積存在空間差異；當(dāng)I=0，表示空間不相關(guān)，即馬鈴薯在空間上隨機(jī)分布。

1.2.3固定效應(yīng)模型

為了研究我國(guó)馬鈴薯區(qū)域格局演變的影響因素，本研究基于2011—2018年14 個(gè)省的面板數(shù)據(jù)，構(gòu)建區(qū)域格局演變影響因素的實(shí)證模型。通過(guò)對(duì)比混合OLS、隨機(jī)效應(yīng)模型和固定效應(yīng)模型的結(jié)果，從而判斷出使用固定效應(yīng)模型更合適。使用混合OLS方法的基本假設(shè)是核心控制變量不存在個(gè)體效應(yīng)，而檢驗(yàn)結(jié)果顯示10 個(gè)地區(qū)的P值在10%的水平下顯著，基本能夠證明存在個(gè)體效應(yīng)，混合OLS模型不適用。隨后對(duì)固定效應(yīng)和隨機(jī)效應(yīng)進(jìn)行比較，檢驗(yàn)結(jié)果顯示P值小于0.01，拒絕隨機(jī)效應(yīng)模型，基本可以判斷固定效應(yīng)模型更加適用。因此本研究使用固定效應(yīng)模型，并且在證明個(gè)體效應(yīng)和時(shí)間效應(yīng)都較為顯著的基礎(chǔ)上，采用雙向固定效應(yīng)模型進(jìn)行研究。實(shí)證計(jì)量模型如下[21-22]：

Y=β0+β1x1it+β2x2it+β3x3it+
β4x4it+β5x5it+μit+φt+εit

(4)

式中：i表示省份；t表示年份；Y表示總產(chǎn)量；β0為常數(shù)項(xiàng)；x1表示氣候指數(shù)；x2表示用工數(shù)量；x3表示化肥費(fèi)用；x4表示種子費(fèi)用；x5表示脫毒種薯應(yīng)用率；μi用于控制地區(qū)固定效應(yīng)；φt用于控制時(shí)間固定效應(yīng)；εit表示誤差項(xiàng)。

1.3 變量選取

本研究在解釋區(qū)域格局演變成因的變量中，將馬鈴薯產(chǎn)量作為被解釋變量，分別從氣候因素、生產(chǎn)投入要素和技術(shù)水平方面選取5 個(gè)解釋變量，包括氣候指數(shù)、用工數(shù)量、脫毒種薯應(yīng)用率、種子費(fèi)用和化肥費(fèi)用。河北、山西、內(nèi)蒙古代表華北，遼寧、黑龍江代表東北，重慶、四川、貴州、云南代表西南，陜西、甘肅、青海、寧夏、新疆代表西北。

1)氣候指數(shù)(x1)：氣候?qū)︸R鈴薯種植具有決定性作用，已有研究成果表明，降水量和氣溫對(duì)馬鈴薯生育期及產(chǎn)量具有重要影響[23-24]。本研究利用氣候的指數(shù)公式，測(cè)算氣候因素對(duì)馬鈴薯產(chǎn)量的影響。其計(jì)算公式為：

式中：Weaj表示j年氣候指數(shù)；Pj表示j年馬鈴薯生育期平均降水量；Tj表示j年馬鈴薯生育期平均氣溫[25]。

2)用工數(shù)量(x2)：近年來(lái)人工成本及其所占比例不斷上升，2018年四川、云南人工成本所占比例超過(guò)60%，已經(jīng)成為總成本的重要組成部分。為解決多重共線(xiàn)性問(wèn)題，本研究將人工成本替換為用工數(shù)量。用工數(shù)量是指馬鈴薯生產(chǎn)過(guò)程中生產(chǎn)者(包括其家庭成員)和雇傭工人直接勞動(dòng)的天數(shù)。

3)化肥費(fèi)用(x3)：化肥是馬鈴薯生產(chǎn)重要投入要素，化肥費(fèi)用是指實(shí)際施用的各種化肥費(fèi)用之和。

4)種子費(fèi)用(x4)：種子是馬鈴薯重要生產(chǎn)資料，種子費(fèi)用是指馬鈴薯種植實(shí)際使用的種薯費(fèi)用。

5)脫毒種薯應(yīng)用率(x5)：本研究選取脫毒種薯應(yīng)用率代表技術(shù)水平。脫毒種薯應(yīng)用率是指脫毒種薯應(yīng)用面積在總種植面積中的比例。

2 結(jié)果與分析

2.1 馬鈴薯種植面積變化分析

統(tǒng)計(jì)局?jǐn)?shù)據(jù)顯示，2016年，我國(guó)馬鈴薯種植面積和產(chǎn)量分別達(dá)到480.24 萬(wàn)hm2和1 698.55 萬(wàn)t，相比1990年，分別增加68%和162%。其中華北、東北、西北和西南地區(qū)種植面積所占比例達(dá)到89%左右，是我國(guó)馬鈴薯主產(chǎn)區(qū)。西南和西北地區(qū)種植面積分別達(dá)到217.01、107.66 萬(wàn)hm2，增幅分別為120%和74%，面積所占比例分別為45%、22%。華北地區(qū)種植面積為77.65 萬(wàn)hm2，增加了31%，東北地區(qū)面積為26.82 萬(wàn)hm2，下降了16%，面積所占比例分別降至16%、6%。

2.2 馬鈴薯格局演變

1990年我國(guó)馬鈴薯種植格局整體呈現(xiàn)“單核心+外圍小中心”的分布特征(圖1(a))。核心種植區(qū)域分布在城口、巫溪、巫山、云陽(yáng)、開(kāi)州、恩施等；“外圍小中心”區(qū)域分布在綏化、齊齊哈爾、大同、烏蘭察布、榆林、張家口、天水、定西、平?jīng)龅取?/p>

2000年種植格局仍然呈現(xiàn)“單核心+外圍小中心”的分布特征(圖1(b))，但隨著烏蘭察布、大同、張家口等地區(qū)核密度的增加，核心區(qū)域轉(zhuǎn)移到了蒙中、晉北和冀西北；重慶、恩施、達(dá)州等地區(qū)核密度下降，由“核心區(qū)”變?yōu)椤巴鈬≈行膮^(qū)”。

2010年種植格局呈“單核心+外圍小中心”分布(圖1(c))，由于固原、定西和平?jīng)龅鹊貐^(qū)核密度逐漸增加，甘東南和寧南成為“核心區(qū)”；隨著烏蘭察布和張家口等地區(qū)核密度下降，蒙中、晉北和冀西北逐漸退出“核心區(qū)”；齊齊哈爾和綏化的核密度下降，黑西南逐漸退出“外圍中心區(qū)”。

2016年，馬鈴薯種植表現(xiàn)為“雙核心+外圍小中心”的分布特征(圖1(d))。近年來(lái)，畢節(jié)、昭通、曲靖、六盤(pán)水、宜賓、遵義等地區(qū)核密度不斷增加，核心區(qū)域轉(zhuǎn)移到云東北和貴西北。定西、平?jīng)觥足y、蘭州、固原和中衛(wèi)等地區(qū)核密度增加，成為第2 個(gè)強(qiáng)核密度區(qū)域。綏化、張家口、大同、朔州、榆林等地區(qū)核密度逐漸下降。

綜上，1990—2016年，我國(guó)馬鈴薯種植格局整體呈現(xiàn)出由“單核心+外圍小中心”向“雙核心+外圍小中心”演變的態(tài)勢(shì)。1990—2016年甘東南、寧南、云東北和貴西北馬鈴薯種植區(qū)域的核密度逐漸增加，蒙中、晉北、冀西北和黑西南逐漸下降。可能的原因是西南一些地區(qū)在氣候、經(jīng)濟(jì)收益、生產(chǎn)技術(shù)和政策等方面具有相對(duì)優(yōu)勢(shì)。西南地區(qū)氣候溫和、水量豐富，馬鈴薯多分布于高原、丘陵，自然條件適合馬鈴薯生長(zhǎng)[26]，一些地區(qū)大規(guī)模發(fā)展冬季馬鈴薯，冬季馬鈴薯商品性和賣(mài)相普遍較好，具有明顯的市場(chǎng)優(yōu)勢(shì)[27-28]。西南地區(qū)馬鈴薯凈利潤(rùn)高于其他糧食作物和油菜等，農(nóng)戶(hù)生產(chǎn)積極性高，大面積冬閑田為馬鈴薯規(guī)模種植提供了土地資源。當(dāng)?shù)赝茝V的免耕覆蓋稻草馬鈴薯栽培技術(shù)，簡(jiǎn)單易行且能夠節(jié)約成本，農(nóng)戶(hù)認(rèn)可度高。云南、貴州等地區(qū)為了脫貧致富相繼出臺(tái)馬鈴薯產(chǎn)業(yè)扶持政策[29-31]。東北地區(qū)馬鈴薯的利潤(rùn)遠(yuǎn)低于水稻，農(nóng)戶(hù)生產(chǎn)積極性不高[31-32]。西北地區(qū)氣候冷涼，馬鈴薯種植氣候適宜性較好。近幾年華北地區(qū)馬鈴薯連作，土傳病害加重發(fā)生，加之冬季馬鈴薯一般2 月份上市，調(diào)節(jié)鮮薯上市時(shí)間、增加市場(chǎng)供給同時(shí)，對(duì)北方一些省份馬鈴薯銷(xiāo)售造成一定沖擊。

目前，我國(guó)馬鈴薯種植集中分布在4 大熱點(diǎn)區(qū)域，分別是以甘肅南部為核心的西北地區(qū)，以云貴川渝為核心的西南地區(qū)，以?xún)?nèi)蒙古中部為核心的華北地區(qū)和以黑龍江為核心的東北地區(qū)，整體由“單核心+外圍小中心”向“雙核心+外圍小中心”演變。西北和西南地區(qū)核密度增加，逐漸發(fā)展為“核心區(qū)”。華北和東北地區(qū)核密度逐漸下降。

2.3 馬鈴薯空間自相關(guān)分析

2.3.1全局空間自相關(guān)分析

以馬鈴薯縣域種植面積數(shù)據(jù)為基礎(chǔ)，運(yùn)用Arc GIS進(jìn)行全局自相關(guān)分析，結(jié)果顯示(表1)：馬鈴薯種植面積的全局莫蘭指數(shù)均大于0，且P值為0，結(jié)果通過(guò)了顯著性檢驗(yàn)，說(shuō)明馬鈴薯種植表現(xiàn)為集聚特征，具有空間正相關(guān)關(guān)系，即種植規(guī)模較大的縣相互臨近，種植規(guī)模較小的縣相互臨近。同時(shí)，莫蘭指數(shù)從1990年的0.677 5下降到2010年的0.494 2，2016年上升到0.598 1，說(shuō)明馬鈴薯種植的空間相關(guān)性先減弱后增強(qiáng)，空間集聚程度先降后升。

表1 馬鈴薯種植面積的全局莫蘭指數(shù)Table 1 Global Moran index of potato planting area

3.3.2局部空間自相關(guān)分析

進(jìn)一步進(jìn)行局部空間自相關(guān)分析，將馬鈴薯種植的集聚類(lèi)型分為高—高、低—低、高—低和低—高等4 類(lèi)。從LISA集聚圖(圖2)來(lái)看，1990—2016年我國(guó)馬鈴薯種植面積變化表現(xiàn)出明顯的局部分異，大部分縣域落在了高—高、低—低和低—高類(lèi)型區(qū)。

高—高類(lèi)型區(qū)為馬鈴薯種植熱點(diǎn)地區(qū)，該類(lèi)型區(qū)內(nèi)縣域內(nèi)和其相鄰縣域馬鈴薯種植面積差異較小，集聚程度都較高。1990年該類(lèi)型區(qū)分布在337 個(gè)縣，包括齊齊哈爾、海拉爾、集寧、包頭、呼和浩特、朔州、天水、定西、隴南、重慶、恩施、西昌、昭通、遵義、畢節(jié)等地區(qū)。2000年減少至314 個(gè)縣，哈爾濱、延安、中衛(wèi)、甘南等地區(qū)分布范圍減小，四平、廣安、寶雞、承德等地區(qū)范圍擴(kuò)大。2010年增至324 個(gè)，綏化、興安盟、松原、大慶、承德、太原、大同等地區(qū)分布范圍減少，慶陽(yáng)、吳忠、自貢、遂寧、內(nèi)江、樂(lè)山、南充、眉山、宜賓、黔南等地區(qū)范圍擴(kuò)大。2016年增至349 個(gè)，黑河、錫林郭勒、張家口、忻州、隴南等地區(qū)空間分布范圍縮小，寶雞、延安、固原、襄陽(yáng)、黔西南等地區(qū)分布范圍擴(kuò)大。由此可知，該類(lèi)型分布范圍縮減的地區(qū)主要分布在我國(guó)東北和華北地區(qū)，范圍擴(kuò)大的主要分布于西南地區(qū)。隨著時(shí)間的變化，高—高類(lèi)型區(qū)域逐漸由東北和華北向西北和西南地區(qū)集聚。

低—低類(lèi)型區(qū)為馬鈴薯種植冷點(diǎn)地區(qū)，縣域自身和其相鄰縣域馬鈴薯種植集聚程度都較小。1990年該類(lèi)型區(qū)分布在542 個(gè)縣，主要包括山東、江蘇、安徽、江西、商丘、周口、杭州、衢州、河源、韶關(guān)等地區(qū)。2000年減少到413 個(gè)，東營(yíng)、煙臺(tái)、威海、寧波、紹興、金華、衢州、武漢、黃石等地區(qū)分布范圍減小，許昌、漯河、信陽(yáng)等地區(qū)分布范圍擴(kuò)大。2010年增至478 個(gè)，杭州、嘉興等地區(qū)范圍縮小，邯鄲、濟(jì)南、東營(yíng)、德州、聊城、濱州、安陽(yáng)、鶴壁、新鄉(xiāng)、濮陽(yáng)等地區(qū)范圍擴(kuò)大。2016年增至690 個(gè)，天津、石家莊、唐山、秦皇島、邢臺(tái)、保定、滄州、廊坊、衡水、上海、臺(tái)州、麗水等地區(qū)分布范圍擴(kuò)大。該類(lèi)型區(qū)基本分布在華東地區(qū)的部分縣，呈現(xiàn)繼續(xù)向東南沿海和華北地區(qū)延伸的趨勢(shì)。

圖2 1990(a)、2000(b)、2010(c)和2016(d)年我國(guó)縣域馬鈴薯種植面積變化LISA圖Fig.2 LISA maps of the change of potato planting area of China in 1990 (a)， 2000 (b)， 2010 (c) and 2016 (d)

低—高類(lèi)型區(qū)縣域自身馬鈴薯種植集聚程度較低，相鄰縣域集聚程度較高，形成局域異質(zhì)“冷點(diǎn)”。1990—2016年該類(lèi)型區(qū)數(shù)量由204減少到187個(gè)縣，西南地區(qū)和內(nèi)蒙古中部的分布范圍縮減，西北地區(qū)的范圍擴(kuò)大。隨著時(shí)間的變化，低—高型區(qū)域逐漸由分散分布向西北地區(qū)集聚，并且呈現(xiàn)向四川中部延伸的趨勢(shì)。高—低類(lèi)型區(qū)縣域自身馬鈴薯種植集聚程度較高，縣域之間種植面積變化空間差異較大，集聚程度較低，形成局域異質(zhì)“熱點(diǎn)”。1990—2016年該類(lèi)型區(qū)的數(shù)量從19增加到23 個(gè)，整體分布比較分散且范圍小，基本分布在部分省的邊界，在東部省份的分布范圍有緩慢增加的趨勢(shì)。

綜上，高—高、低—低集聚區(qū)分布范圍擴(kuò)大，低—高集聚區(qū)范圍縮小。高—高集聚區(qū)由相對(duì)分散分布逐漸向西北和西南地區(qū)集中，低—低集聚區(qū)基本分布在華東地區(qū)部分縣，呈現(xiàn)繼續(xù)向東南沿海和華北地區(qū)延伸的趨勢(shì)，低—高集聚區(qū)逐漸由分散分布向西北地區(qū)集中并且呈現(xiàn)繼續(xù)向四川中部延伸的趨勢(shì)。以上分析證明東北和華北地區(qū)馬鈴薯種植的集聚程度逐漸下降，西北和西南地區(qū)集聚程度逐漸提高。

3 馬鈴薯區(qū)域格局演變的影響因素分析

3.1 分樣本回歸結(jié)果分析

3.1.1空間分樣本回歸

通過(guò)空間分樣本回歸，對(duì)不同地區(qū)馬鈴薯生產(chǎn)的影響因素進(jìn)行對(duì)比分析。結(jié)果顯示(表2)，2011—2018年，氣候指數(shù)、用工數(shù)量和脫毒種薯應(yīng)用率對(duì)4 大地區(qū)馬鈴薯生產(chǎn)均有顯著影響。氣候指數(shù)和用工數(shù)量對(duì)華北影響最大，影響系數(shù)分別為0.213 4和1.023 9；對(duì)西南影響最小，影響系數(shù)分別為0.164 3和0.792 7。脫毒種薯應(yīng)用率對(duì)西南影響最大，對(duì)華北影響最小，系數(shù)分別為0.084 1和0.093 1。

地區(qū)樣本回歸結(jié)果說(shuō)明，近幾年技術(shù)水平、氣候和勞動(dòng)力投入對(duì)主產(chǎn)區(qū)馬鈴薯種植區(qū)域格局具有顯著影響。氣候和勞動(dòng)力投入對(duì)華北地區(qū)的影響大于其他地區(qū)，技術(shù)水平對(duì)西南地區(qū)的影響大于其他地區(qū)，但影響程度的區(qū)域差異較小。

表2 馬鈴薯區(qū)域格局演變影響因素地區(qū)分樣本回歸結(jié)果Table 2 Factors affecting the evolution of potato regional pattern

3.1.2時(shí)間分樣本回歸

通過(guò)時(shí)間分樣本回歸，分析不同時(shí)間段馬鈴薯生產(chǎn)影響因素的變化情況。結(jié)果顯示(表3)，2011—2014年和2015—2018年氣候指數(shù)、用工數(shù)量、脫毒種薯應(yīng)用率對(duì)馬鈴薯生產(chǎn)具有顯著影響。氣候的影響系數(shù)從0.158 1降至0.156 8，脫毒種薯應(yīng)用率的系數(shù)從0.086 2降至0.071 3，用工數(shù)量的系數(shù)從1.055 0增加到1.063 9，高于其他因素。

時(shí)間樣本回歸結(jié)果說(shuō)明，2011—2018年技術(shù)水平、氣候和勞動(dòng)力投入對(duì)馬鈴薯種植區(qū)域格局具有顯著影響。氣候和技術(shù)水平影響程度有所下降，勞動(dòng)力投入影響程度增加，且高于其他因素，說(shuō)明勞動(dòng)力的投入是這段時(shí)期馬鈴薯種植區(qū)域格局的重要影響因素。

表3 馬鈴薯區(qū)域格局演變影響因素時(shí)間分樣本回歸結(jié)果Table 3 Factors influencing the evolution of potato regional pattern

3.1.3空間與時(shí)間分樣本回歸

為了得出綜合性結(jié)論，再次對(duì)華北、東北、西南和西北地區(qū)進(jìn)行時(shí)間分樣本回歸，結(jié)果顯示(表4)，2011—2014年和2015—2018年氣候指數(shù)、用工數(shù)量、脫毒種薯應(yīng)用率對(duì)各地區(qū)馬鈴薯生產(chǎn)具有顯著影響。各地區(qū)氣候指數(shù)的影響系數(shù)下降，用工數(shù)量的影響系數(shù)上升，除華北地區(qū)，脫毒種薯應(yīng)用率影響系數(shù)下降。從影響因素區(qū)域?qū)Ρ惹闆r來(lái)看，氣候指數(shù)、用工數(shù)量和脫毒種薯應(yīng)用率對(duì)各地區(qū)馬鈴薯生產(chǎn)具有顯著影響，華北地區(qū)氣候指數(shù)和用工數(shù)量的系數(shù)大于其他地區(qū)，2011—2014年系數(shù)分別為0.201 8和1.286 5，2015—2018年分別為0.198 4和1.294 9；脫毒種薯應(yīng)用率對(duì)西南的影響大于同時(shí)期其他地區(qū)，兩段時(shí)期影響系數(shù)分別為0.092 7和0.078 1。

以上說(shuō)明，近幾年氣候、勞動(dòng)力投入和技術(shù)水平對(duì)我國(guó)主產(chǎn)區(qū)馬鈴薯種植區(qū)域格局具有顯著影響，勞動(dòng)力投入的影響程度大于其他因素。在時(shí)間上，技術(shù)水平和氣候的影響程度下降，勞動(dòng)力投入的影響程度上升。在空間上，氣候和勞動(dòng)力投入對(duì)華北的影響最大，技術(shù)水平對(duì)西南的影響最大，但影響程度區(qū)域差異較小。

4 結(jié)論和政策建議

本研究采用空間統(tǒng)計(jì)分析方法和縣級(jí)統(tǒng)計(jì)數(shù)據(jù)，定量分析了我國(guó)1990—2016年馬鈴薯區(qū)域格局演變趨勢(shì)。基于2011—2018年14 個(gè)省的面板數(shù)據(jù)，運(yùn)用雙向固定效應(yīng)模型分析區(qū)域格局演變的影響因素。得出以下結(jié)論：1)我國(guó)馬鈴薯種植集中分布在4 大熱點(diǎn)區(qū)域，分別是以甘肅南部為核心的西北地區(qū)，以云貴川渝為核心的西南地區(qū)，以?xún)?nèi)蒙古中部為核心的華北地區(qū)和以黑龍江為核心的東北地區(qū)，整體由“單核心+外圍小中心”向“雙核心+外圍小中心”演變。西北和西南地區(qū)核密度增加，逐漸發(fā)展為“核心區(qū)”。華北和東北地區(qū)核密度逐漸下降；2)我國(guó)馬鈴薯種植的空間分布表現(xiàn)為顯著集聚，集聚程度先減后增，主要表現(xiàn)出高—高、低—高和低—低集聚類(lèi)型，東北和華北地區(qū)馬鈴薯種植的集聚程度逐漸下降，西北和西南地區(qū)集聚程度逐漸提高。高—高集聚區(qū)由相對(duì)分散分布逐漸向西北和西南地區(qū)集中，低—低集聚區(qū)基本分布在華東地區(qū)，呈現(xiàn)繼續(xù)向東南沿海和華北地區(qū)延伸的趨勢(shì)，低—高集聚區(qū)逐漸由分散分布向西北地區(qū)集中并且呈現(xiàn)繼續(xù)向四川中部延伸的趨勢(shì)；3)我國(guó)馬鈴薯區(qū)域格局演變主要受氣候、勞動(dòng)力投入和技術(shù)水平等要素的影響，其中勞動(dòng)力投入的影響程度最大。短期內(nèi)各影響因素的影響程度會(huì)因時(shí)間和空間的差異而有所不同，但差異較小。

基于以上研究結(jié)論，本研究提出以下政策建議：

1)制定全國(guó)馬鈴薯綜合區(qū)劃方案，為優(yōu)化空間布局提供參考；2)大范圍推廣脫毒種薯、提升馬鈴薯種植技術(shù)水平有助于提高馬鈴薯產(chǎn)量和質(zhì)量；3)建議實(shí)施良種和機(jī)械等補(bǔ)貼。目前馬鈴薯產(chǎn)業(yè)已經(jīng)成為部分地區(qū)脫貧增收的支柱產(chǎn)業(yè)，但相比三大主糧，政策支持力度小，政策的支持有利于調(diào)動(dòng)生產(chǎn)積極性，促進(jìn)馬鈴薯產(chǎn)業(yè)發(fā)展；4)加快建設(shè)馬鈴薯大數(shù)據(jù)平臺(tái)，實(shí)時(shí)傳遞價(jià)格和市場(chǎng)供求信息，以合理調(diào)控馬鈴薯生產(chǎn)規(guī)模、規(guī)避市場(chǎng)風(fēng)險(xiǎn)和保證薯農(nóng)收益。