高明杰，魯洪威，李婷婷，薄沁箐，羅其友

中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所，北京 100081

馬鈴薯是全世界人民喜愛的糧菜兼用作物，為世界人民戰(zhàn)勝饑餓、提供均衡營養(yǎng)作出了重要貢獻(xiàn)[1]。馬鈴薯抗逆性強(qiáng)，比較效益顯著，目前在我國大部分省(自治區(qū))皆有種植，為我國第四大糧食作物[2]。2015年馬鈴薯作為未來主糧被農(nóng)業(yè)農(nóng)村部正式提出，標(biāo)志著馬鈴薯成為繼稻米、小麥、玉米之后保障國家糧食安全的重要農(nóng)產(chǎn)品[3]。根據(jù)聯(lián)合國糧農(nóng)組織(FAO)官方數(shù)據(jù)顯示，我國在2001年后成為世界馬鈴薯第一大生產(chǎn)國，面積和產(chǎn)量均居世界首位[4]。而我國馬鈴薯生產(chǎn)的不穩(wěn)定性較強(qiáng)，產(chǎn)量在年際之間的大幅波動(dòng)時(shí)有發(fā)生，從產(chǎn)業(yè)鏈協(xié)調(diào)發(fā)展的角度看，馬鈴薯生產(chǎn)的波動(dòng)對全產(chǎn)業(yè)鏈的健康發(fā)展構(gòu)成了挑戰(zhàn)。因此，對我國馬鈴薯生產(chǎn)波動(dòng)影響因素進(jìn)行研究，深入分析各種因素對我國馬鈴薯生產(chǎn)的作用機(jī)制，對于保持我國馬鈴薯產(chǎn)業(yè)平穩(wěn)健康發(fā)展具有十分重要的現(xiàn)實(shí)意義。

農(nóng)業(yè)是自然和社會(huì)經(jīng)濟(jì)再生產(chǎn)的過程[5]，通過對過往研究文獻(xiàn)的梳理，發(fā)現(xiàn)影響產(chǎn)量波動(dòng)的因素主要包括自然條件、科技、社會(huì)經(jīng)濟(jì)和政策等四類。在產(chǎn)量波動(dòng)的時(shí)間序列中往往伴隨著受科技進(jìn)步影響的單產(chǎn)波動(dòng)和農(nóng)戶種植決策導(dǎo)致的面積波動(dòng)[6-8]。農(nóng)業(yè)產(chǎn)量的高低與當(dāng)年自然條件的狀況密切相關(guān)，而氣候條件的不可預(yù)知性就決定了農(nóng)民在農(nóng)業(yè)生產(chǎn)的過程中應(yīng)對自然條件變化的困難，諸如旱澇災(zāi)害、病蟲害、生長期日照時(shí)長等常對農(nóng)業(yè)生產(chǎn)造成沖擊[9-10]?；?、農(nóng)業(yè)機(jī)械總動(dòng)力、機(jī)耕面積和有效灌溉面積等體現(xiàn)科技水平的變量對總產(chǎn)量波動(dòng)具有十分顯著的影響[11-12]。而其他學(xué)者[13-15]通過實(shí)證研究，發(fā)現(xiàn)了價(jià)格對產(chǎn)量波動(dòng)的影響關(guān)系，進(jìn)一步厘清了影響農(nóng)產(chǎn)品產(chǎn)量波動(dòng)的社會(huì)經(jīng)濟(jì)因素。政策亦是對生產(chǎn)影響因素分析不可忽視的要素。黃季焜等[16]認(rèn)為政策因素對農(nóng)產(chǎn)品生產(chǎn)的影響是直接且深遠(yuǎn)的，往往會(huì)對此后數(shù)年整體的農(nóng)業(yè)產(chǎn)業(yè)發(fā)展產(chǎn)生深遠(yuǎn)影響。

現(xiàn)有研究主要關(guān)注于糧食產(chǎn)量的影響因素分析，而對某些重要農(nóng)產(chǎn)品的單獨(dú)研究較少，對馬鈴薯產(chǎn)量影響的機(jī)制研究更為匱乏，并且在影響因素的選取中往往側(cè)重于自然、科技、社會(huì)經(jīng)濟(jì)和政策某一或幾個(gè)方面，缺乏對農(nóng)產(chǎn)品生產(chǎn)影響因素的系統(tǒng)性研究。馬鈴薯生產(chǎn)有著自身的波動(dòng)規(guī)律與影響機(jī)制，需要有針對性研究才能得到更有價(jià)值的研究成果。本研究通過構(gòu)建較為全面的影響因素指標(biāo)體系，運(yùn)用向量自回歸模型對1982-2018年我國馬鈴薯產(chǎn)量波動(dòng)影響作用機(jī)制進(jìn)行研究，以期為指導(dǎo)我國馬鈴薯生產(chǎn)平穩(wěn)發(fā)展提供決策參考，推動(dòng)馬鈴薯主糧化戰(zhàn)略順利實(shí)施，提升我國糧食安全保障能力。

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

1.1 VAR模型

向量自回歸模型簡稱VAR模型，1980年由克里斯托弗·西姆斯(Christopher Sims)提出[17]。VAR模型是用模型中所有當(dāng)期變量對所有變量的若干滯后變量進(jìn)行回歸，目前此模型已得到廣泛應(yīng)用[18-20]。

一般表達(dá)式為：

Yt=A1Yt-1+A2Yt-2+…+ApYt-p+

B1Xt-1+…+BrXt-r+εt

(1)

式(1)中，Yt是馬鈴薯產(chǎn)量各影響因素m維內(nèi)生變量向量，Xt是各影響因素d維外生變量向量，內(nèi)生變量與外生滯后期分別為p和r階(A1,A2,……Ap)和(B1,B2,……Br)是待估矩陣；εt是隨機(jī)項(xiàng)。

采用向量自回歸模型(VAR)研究從1982-2018年包含自然、科技、社會(huì)和政策等各變量對我國馬鈴薯產(chǎn)量波動(dòng)的影響，進(jìn)而得出各影響因素之間的動(dòng)態(tài)互動(dòng)關(guān)系。

在運(yùn)用時(shí)間序列數(shù)據(jù)建立VAR模型的過程中，如果數(shù)據(jù)不平穩(wěn)，會(huì)導(dǎo)致出現(xiàn)偽回歸結(jié)果，使得分析結(jié)論失去實(shí)際含義[21]。在建模前我們先采用ADF檢驗(yàn)法對不同變量時(shí)序數(shù)據(jù)進(jìn)行平穩(wěn)分析以檢測各時(shí)序數(shù)據(jù)是否平穩(wěn)。在此之前需剔除數(shù)據(jù)異方差的影響，首先對不包括政策變量的其余所有自變量和因變量分別取對數(shù)。在確定完VAR模型的滯后階數(shù)，對模型協(xié)整檢驗(yàn)通過后，還需要檢驗(yàn)VAR模型的穩(wěn)定性，才能確定建立的模型與實(shí)際相符。AR根圖法可以進(jìn)一步檢驗(yàn)方程是否存在一定的穩(wěn)定性。VAR模型不是理論分析模型，在實(shí)際應(yīng)用中，常常不分析變量間的影響，而是分析模型受到特定沖擊時(shí)對經(jīng)濟(jì)系統(tǒng)產(chǎn)生的影響，我們稱之為脈沖響應(yīng)法[22-25]。運(yùn)用脈沖響應(yīng)分析能更直觀地說明各種變量之間的動(dòng)態(tài)交互關(guān)系,發(fā)現(xiàn)在滯后10期后脈沖響應(yīng)量趨于穩(wěn)定,因此，選擇脈沖響應(yīng)的滯后期數(shù)為10。

1.2 數(shù)據(jù)來源與指標(biāo)構(gòu)建

1)數(shù)據(jù)來源。本研究涉及的1982-2018年我國馬鈴薯產(chǎn)量、成災(zāi)面積、耕地面積、化肥施用量、農(nóng)業(yè)機(jī)械總動(dòng)力、機(jī)耕面積、有效灌溉面積、農(nóng)村人均可支配收入、城鎮(zhèn)化水平、和價(jià)格指數(shù)數(shù)據(jù)來自歷年的《中國農(nóng)村統(tǒng)計(jì)年鑒》《中國統(tǒng)計(jì)年鑒》《中國農(nóng)業(yè)統(tǒng)計(jì)年鑒》《新中國60年統(tǒng)計(jì)資料匯編》《中國機(jī)械工業(yè)年鑒》和國家統(tǒng)計(jì)局“年度分省數(shù)據(jù)庫”。涉及的年均降雨量、年均氣溫和日照時(shí)長來自國家氣象信息中心各地區(qū)地面氣象監(jiān)測站數(shù)據(jù)要素集。

2)指標(biāo)構(gòu)建。參考已有文獻(xiàn)對生產(chǎn)波動(dòng)影響因素的選取，根據(jù)指標(biāo)的代表性和可獲得性，從影響馬鈴薯產(chǎn)量的自然、科技、社會(huì)經(jīng)濟(jì)和政策四個(gè)方面構(gòu)建產(chǎn)量波動(dòng)影響因素指標(biāo)體系。其中自然條件指標(biāo)包括年均降雨量X1(mm)、年均氣溫X2(℃)、日照時(shí)長X3(h)、馬鈴薯成災(zāi)面積X4(khm2)、耕地面積X5(khm2)；科技因素指標(biāo)包括馬鈴薯化肥施用量X6(萬t)、馬鈴薯農(nóng)業(yè)機(jī)械總動(dòng)力X7(萬kW)、馬鈴薯機(jī)耕面積X8(khm2)、馬鈴薯有效灌溉面積X9(khm2)；社會(huì)經(jīng)濟(jì)因素指標(biāo)包括農(nóng)村人均可支配收入X10(元)、城鎮(zhèn)化水平X11(%)、農(nóng)業(yè)生產(chǎn)資料價(jià)格指數(shù)X12(%)、農(nóng)業(yè)生產(chǎn)者價(jià)格指數(shù)X13(%)；政策因素指標(biāo)采用虛擬變量X14(0，1)；馬鈴薯產(chǎn)量指標(biāo)采用歷年馬鈴薯總產(chǎn)量Y(折糧/萬t)。

其中馬鈴薯成災(zāi)面積、馬鈴薯化肥施用量、馬鈴薯農(nóng)業(yè)機(jī)械總動(dòng)力、馬鈴薯機(jī)耕面積和馬鈴薯有效灌溉面積指標(biāo)無法從統(tǒng)計(jì)數(shù)據(jù)中直接獲取，本研究進(jìn)行了估算。估算方法為(以馬鈴薯成災(zāi)面積為例)：馬鈴薯成災(zāi)面積=馬鈴薯種植面積/農(nóng)作物種植面積×成災(zāi)面積，其余數(shù)據(jù)估算方法相同。

2 結(jié)果與分析

2.1 中國馬鈴薯產(chǎn)量時(shí)序變化

我國馬鈴薯生產(chǎn)經(jīng)歷了一個(gè)較為曲折的發(fā)展過程，從有官方統(tǒng)計(jì)數(shù)據(jù)的1982年以來我國馬鈴薯總產(chǎn)量從最初的不足500萬t(折糧，下同)發(fā)展到2018年的近1 800萬t，年均增長3.76%。由圖1可知，我國馬鈴薯總產(chǎn)量是一個(gè)波動(dòng)增長的過程，不同年份產(chǎn)量波動(dòng)較為明顯。根據(jù)馬鈴薯產(chǎn)量增速變化可知，1982年以來共有22 a我國馬鈴薯生產(chǎn)處于擴(kuò)張期，占總年數(shù)的61.11%，其余14 a馬鈴薯產(chǎn)量為負(fù)增長，處于收縮期。我國馬鈴薯產(chǎn)量年際間波動(dòng)較為顯著，而影響我國馬鈴薯產(chǎn)量時(shí)序變化的原因并無量化分析的科學(xué)定論。因此，有必要構(gòu)建計(jì)量模型定量分析馬鈴薯產(chǎn)量波動(dòng)影響機(jī)制。

圖1 我國馬鈴薯產(chǎn)量時(shí)序變化

2.2 中國馬鈴薯產(chǎn)量波動(dòng)影響因素分析

1)平穩(wěn)性檢驗(yàn)。由表1可知，原始序列是一階單整的，符合VAR模型構(gòu)建的要求，可以對各變量構(gòu)建向量自回歸模型。

表1 影響因素指標(biāo)平穩(wěn)性檢驗(yàn) Table 1 Influencing factor index stationarity test

2)滯后階數(shù)確定和VAR模型結(jié)果。確定VAR模型的滯后期是關(guān)系模型有效性的關(guān)鍵環(huán)節(jié)，由表2中AIC和SC信息值可知，模型的滯后階數(shù)取值為1。據(jù)此對各變量進(jìn)行一階滯后向量自回歸，VAR模型回歸方程的可決系數(shù)R2基本都大于90%，說明模型的擬合度較好，VAR模型可有效反映各變量間的動(dòng)態(tài)互動(dòng)關(guān)系。

表2 模型滯后階數(shù) Table 2 Model lag order

3)方程穩(wěn)定性檢驗(yàn)。由圖2可見，檢驗(yàn)后發(fā)現(xiàn)所有特征根都位于單位圓內(nèi)，即AR根模的倒數(shù)皆小于1，證明所構(gòu)造的VAR模型具有穩(wěn)定性，因此后續(xù)可進(jìn)行脈沖響應(yīng)研究和方差分解分析。

圖2 VAR模型穩(wěn)定性檢驗(yàn)

4)脈沖響應(yīng)分析。將各變量對我國馬鈴薯產(chǎn)量的影響進(jìn)行脈沖響應(yīng)分析，結(jié)果如圖3所示。

在自然因素影響方面，由圖3A可知，年均降雨量(X1)對馬鈴薯產(chǎn)量有較為明顯的負(fù)向影響，作為一般需水性作物，在馬鈴薯生長期僅需要400 mm左右的降雨即可滿足其正常生長，在播種期降雨多會(huì)導(dǎo)致種薯塊莖壞死，生長期雨水過多往往會(huì)導(dǎo)致晚疫病的發(fā)生。由圖3B可知，年均氣溫(X2)對馬鈴薯產(chǎn)量表現(xiàn)為負(fù)向影響，但影響并不顯著，各地因氣候條件的不同，氣溫略有差別，但年際間的變化并不明顯。由圖3C可知，從長期來看日照時(shí)長(X3)對馬鈴薯產(chǎn)量有正向影響作用。日照時(shí)長關(guān)系到營養(yǎng)成分的積累，進(jìn)而影響產(chǎn)量。由圖3D可知，馬鈴薯成災(zāi)面積(X4)對馬鈴薯產(chǎn)量表現(xiàn)為正向效應(yīng)，馬鈴薯成災(zāi)面積較大，當(dāng)年價(jià)格行情較好，使得馬鈴薯種植獲利空間增大，種植面積有所增加。由圖3E可知，耕地面積(X5)對馬鈴薯產(chǎn)量影響并不明顯，馬鈴薯種植面積占總耕地面積的比例僅為4%左右，且我國耕地面積的減少主要集中在東部地區(qū)的經(jīng)濟(jì)發(fā)達(dá)地區(qū)，而馬鈴薯大多種植在中西部地區(qū)，因此，全國耕地面積的變化對馬鈴薯總產(chǎn)量的影響并不明顯。

實(shí)線為脈沖響應(yīng)函數(shù)圖像，表示馬鈴薯產(chǎn)量對影響因素變量沖擊產(chǎn)生的響應(yīng)；虛線表示正負(fù)兩倍的標(biāo)準(zhǔn)差偏離區(qū)間。The solid line is the image of impulse response function,which represents the response of potato yield to the impact of influencing factors. The dotted line indicates a deviation of plus or minus two times the standard deviation.

在科技因素影響方面，由圖3F可知，化肥施用量(X6)對馬鈴薯產(chǎn)量沖擊后，馬鈴薯產(chǎn)量在第二期出現(xiàn)了顯著增長，第三期后影響效應(yīng)逐漸減弱，說明馬鈴薯化肥施用量增加對馬鈴薯產(chǎn)量的正向影響是較為明顯的，肥料對馬鈴薯增產(chǎn)響應(yīng)效果顯著。由圖3G可知，馬鈴薯(X7)農(nóng)業(yè)機(jī)械總動(dòng)力的影響效應(yīng)為負(fù)，農(nóng)業(yè)機(jī)械的使用往往伴隨著農(nóng)業(yè)生產(chǎn)規(guī)模的增加，表現(xiàn)為農(nóng)業(yè)生產(chǎn)由精耕細(xì)作向粗放式發(fā)展。由機(jī)械總動(dòng)力增加所體現(xiàn)的農(nóng)業(yè)生產(chǎn)集約化水平的提升，對馬鈴薯總產(chǎn)量具有一定的負(fù)面沖擊。由圖3H、3I可知，馬鈴薯機(jī)耕面積(X8)和有效灌溉面積(X9)對馬鈴薯產(chǎn)量在第二期、第三期均表現(xiàn)為正向效應(yīng)，并在此后效應(yīng)逐漸減弱。馬鈴薯機(jī)械化耕作和在生長期內(nèi)進(jìn)行有效的水分灌溉會(huì)對馬鈴薯增產(chǎn)產(chǎn)生一定的正向效應(yīng)，機(jī)耕面積和有效灌溉面積增加表明我國農(nóng)業(yè)現(xiàn)代化水平的整體提升。

在社會(huì)經(jīng)濟(jì)因素影響方面，由圖3J和3K可知，農(nóng)村人均可支配收入(X10)和城鎮(zhèn)化水平(X11)對馬鈴薯產(chǎn)量沖擊后產(chǎn)量在第二期都表現(xiàn)為增加，并在此后漸趨穩(wěn)定。隨著農(nóng)村人均可支配收入的增加，在市場效益的驅(qū)動(dòng)下，農(nóng)民增加了馬鈴薯這種高投入、高產(chǎn)出的作物種植。并且隨著城鎮(zhèn)化水平的提高，馬鈴薯多樣化食品的市場需求逐漸增加。由圖3L和3M可知，農(nóng)業(yè)生產(chǎn)資料價(jià)格指數(shù)(X12)和農(nóng)業(yè)生產(chǎn)者價(jià)格指數(shù)(X13)對馬鈴薯產(chǎn)量沖擊后，馬鈴薯產(chǎn)量分別表現(xiàn)為負(fù)向效應(yīng)和正向效應(yīng)。農(nóng)業(yè)生產(chǎn)資料價(jià)格指數(shù)增加意味著馬鈴薯生產(chǎn)成本的提升，壓縮了馬鈴薯種植利潤，影響農(nóng)戶種植意愿，馬鈴薯產(chǎn)量最終減少。而農(nóng)業(yè)生產(chǎn)者價(jià)格指數(shù)增加，意味著馬鈴薯市場價(jià)格行情向好，在市場利益的驅(qū)動(dòng)下，農(nóng)戶擴(kuò)大了馬鈴薯的種植進(jìn)而增加了馬鈴薯產(chǎn)量。

在政策因素影響方面，由圖3N可知，政府政策(X14)對馬鈴薯產(chǎn)量施加一個(gè)標(biāo)準(zhǔn)差沖擊后，馬鈴薯產(chǎn)量在第二期、第三期開始增加，并在此后影響趨勢逐漸減弱。政府對農(nóng)業(yè)產(chǎn)業(yè)的支持政策往往會(huì)促進(jìn)某一農(nóng)產(chǎn)品的生產(chǎn)進(jìn)而導(dǎo)致農(nóng)產(chǎn)品產(chǎn)量的增加，但在市場需求沒有擴(kuò)大的情況下，產(chǎn)量增加最終都會(huì)使得價(jià)格下降，并導(dǎo)致收益減少，因此影響產(chǎn)量的仍然是市場效益。

5)方差貢獻(xiàn)率分析。由圖4可知，在不同滯后時(shí)期各因素對馬鈴薯生產(chǎn)影響程度有所不同，且各變量方差貢獻(xiàn)率變化趨勢也不盡相同。在各變量中，年均降雨量(X1)在第二期對馬鈴薯產(chǎn)量方差貢獻(xiàn)率最大(3.3%)。這說明在當(dāng)期降雨量較其他因素對馬鈴薯產(chǎn)量的影響更顯著，且降雨量對馬鈴薯生產(chǎn)的影響在第二期達(dá)到最大，隨后逐年影響效應(yīng)減弱。

此外，在第二期馬鈴薯成災(zāi)面積(X4)對馬鈴薯生產(chǎn)的沖擊僅次于降雨量，方差貢獻(xiàn)率為1.63%，說明馬鈴薯成災(zāi)面積對馬鈴薯生產(chǎn)隨后的沖擊影響同樣很大。由于自然氣候因素造成的馬鈴薯成災(zāi)面積對馬鈴薯的產(chǎn)量有較為直接的影響，負(fù)面沖擊往往在當(dāng)年就能顯現(xiàn)出來，遭受了自然災(zāi)害沖擊的馬鈴薯產(chǎn)量會(huì)出現(xiàn)不同幅度的下降。

年均氣溫(X2)對馬鈴薯產(chǎn)量的影響在第二期、第三期達(dá)到最大，此后影響效應(yīng)逐漸減弱。但氣溫在所有變量的方差貢獻(xiàn)率中占比重較小，說明相比于其他影響因素，氣溫在馬鈴薯產(chǎn)量總體波動(dòng)中的影響效應(yīng)非常有限，這與年際間氣溫的穩(wěn)定性密不可分。

日照時(shí)長(X3)對馬鈴薯產(chǎn)量沖擊的方差貢獻(xiàn)率較小，且不同時(shí)期對馬鈴薯產(chǎn)量的影響較為穩(wěn)定。我國幅員遼闊，自然條件空間差異大，日照時(shí)長對個(gè)別地區(qū)馬鈴薯產(chǎn)量的影響較為明顯，而放眼全國這種影響效應(yīng)則較為穩(wěn)定。

耕地面積(X5)和化肥施用量(X6)對馬鈴薯產(chǎn)量的方差貢獻(xiàn)率較小，并且不同時(shí)期的變化較為穩(wěn)定，說明耕地面積和化肥投入對馬鈴薯產(chǎn)量的影響效應(yīng)較為穩(wěn)定。馬鈴薯農(nóng)業(yè)機(jī)械總動(dòng)力(X7)對總產(chǎn)量的沖擊在最初并不顯著，但隨著時(shí)間的推移，影響效應(yīng)有不斷增大的趨勢。隨著以馬鈴薯機(jī)械總動(dòng)力為表征的馬鈴薯現(xiàn)代化耕種水平不斷提升，馬鈴薯產(chǎn)量也會(huì)受到影響，并且現(xiàn)代化水平較其他因素對馬鈴薯產(chǎn)量的影響效應(yīng)更為長久。馬鈴薯機(jī)耕面積(X8)對馬鈴薯產(chǎn)量的影響效應(yīng)較小且穩(wěn)定，說明機(jī)械耕作對馬鈴薯產(chǎn)量的總體影響并不明顯。

馬鈴薯有效灌溉面積(X9)對總產(chǎn)量沖擊的方差貢獻(xiàn)率在最初占比較小，但隨著時(shí)間的推移，方差貢獻(xiàn)率有逐漸擴(kuò)大的趨勢，說明有效灌溉面積對馬鈴薯產(chǎn)量的影響是一個(gè)逐步加深的過程，馬鈴薯生產(chǎn)基礎(chǔ)設(shè)施的建設(shè)與完善對馬鈴薯產(chǎn)量的影響是深遠(yuǎn)且廣泛的。

由圖4可知，以農(nóng)村人均可支配收入(X10)和城鎮(zhèn)化水平(X11)為表征的社會(huì)發(fā)展水平對馬鈴薯產(chǎn)量的影響效應(yīng)較小，說明社會(huì)發(fā)展水平對馬鈴薯產(chǎn)量的影響并不是決定性因素，伴隨著農(nóng)村人均可支配收入的增加和城鎮(zhèn)化水平的提升，以馬鈴薯為主糧的部分經(jīng)濟(jì)不發(fā)達(dá)地區(qū)居民會(huì)減少對馬鈴薯的消費(fèi)，但由于社會(huì)發(fā)展水平的提升，人們對馬鈴薯精深加工品的需求增加，又刺激了高品質(zhì)馬鈴薯的生產(chǎn)，兩種力量兩相抵消，最終表現(xiàn)為社會(huì)發(fā)展水平對馬鈴薯產(chǎn)量的影響并不顯著。

圖4 各變量不同滯后期的方差貢獻(xiàn)率

市場因素中農(nóng)業(yè)生產(chǎn)資料價(jià)格指數(shù)(X12)對馬鈴薯產(chǎn)量的影響效應(yīng)在第三期達(dá)到最大，并在此后影響效應(yīng)逐步減弱。農(nóng)業(yè)生產(chǎn)資料價(jià)格指數(shù)在第三期對馬鈴薯產(chǎn)量的方差貢獻(xiàn)率僅次于降雨量和成災(zāi)面積。說明由農(nóng)業(yè)生產(chǎn)資料價(jià)格指數(shù)所體現(xiàn)的馬鈴薯生產(chǎn)成本變化對馬鈴薯產(chǎn)量在第三期有較為顯著的影響，且對種植決策的影響存在一定的時(shí)滯。農(nóng)業(yè)生產(chǎn)者價(jià)格指數(shù)(X13)對馬鈴薯產(chǎn)量的影響效應(yīng)較小，且方差貢獻(xiàn)率較為穩(wěn)定。政策因素(X14)對馬鈴薯產(chǎn)量波動(dòng)的影響效應(yīng)隨著時(shí)間的推移有增大趨勢，表明政策對馬鈴薯產(chǎn)量有較為長期的影響，并存在一定的時(shí)滯。

3 研究結(jié)論與對策建議

3.1 研究結(jié)論

通過VAR模型對1982-2018年我國馬鈴薯產(chǎn)量波動(dòng)影響因素進(jìn)行量化研究，得到馬鈴薯產(chǎn)量波動(dòng)與自然、科技、社會(huì)和政策系統(tǒng)要素間存在的動(dòng)態(tài)響應(yīng)規(guī)律，具體如下：

1)自然因素對產(chǎn)量波動(dòng)多為負(fù)向影響，經(jīng)濟(jì)、科技和政策以正向影響為主。降雨量增多、年均氣溫增加對馬鈴薯生產(chǎn)有負(fù)向影響，日照時(shí)長增加對馬鈴薯產(chǎn)量有正向作用，耕地面積對產(chǎn)量影響不顯著，化肥施用量對馬鈴薯產(chǎn)量有顯著正向影響，集約化水平對馬鈴薯總產(chǎn)量有負(fù)面沖擊，馬鈴薯機(jī)耕面積、有效灌溉面積和政府政策都為正向效應(yīng)，農(nóng)村人均可支配收入和城鎮(zhèn)化水平都為正向效應(yīng)，農(nóng)業(yè)生產(chǎn)資料價(jià)格指數(shù)和農(nóng)業(yè)生產(chǎn)者價(jià)格指數(shù)分別表現(xiàn)為負(fù)向效應(yīng)和正向效應(yīng)。

2)自然、科技和市場因素對馬鈴薯產(chǎn)量波動(dòng)影響貢獻(xiàn)較大，科技進(jìn)步與政策影響較為長遠(yuǎn)。在第二期降雨量對馬鈴薯生產(chǎn)的影響最大；馬鈴薯成災(zāi)面積對馬鈴薯生產(chǎn)的方差貢獻(xiàn)率較大；年均氣溫對馬鈴薯產(chǎn)量的影響先增后減；日照時(shí)長對產(chǎn)量影響較為有限；耕地面積和化肥投入影響效應(yīng)較為穩(wěn)定；馬鈴薯農(nóng)業(yè)機(jī)械總動(dòng)力開始影響較小，但隨時(shí)間推移影響效應(yīng)增大；機(jī)耕因素貢獻(xiàn)并不顯著；有效灌溉面積影響逐步加深；社會(huì)發(fā)展水平對馬鈴薯產(chǎn)量影響效應(yīng)較?。晦r(nóng)業(yè)生產(chǎn)資料價(jià)格指數(shù)在第三期有較顯著影響；生產(chǎn)者價(jià)格因素影響程度一般；政策的長期影響存在時(shí)滯。

3.2 對策建議

1)尊重馬鈴薯生產(chǎn)波動(dòng)規(guī)律，制定適當(dāng)生產(chǎn)引導(dǎo)政策。了解馬鈴薯生產(chǎn)波動(dòng)客觀規(guī)律，并在尊重馬鈴薯生產(chǎn)波動(dòng)規(guī)律的基礎(chǔ)之上發(fā)揮政策和措施的引導(dǎo)作用。全國各級(jí)農(nóng)業(yè)主管部門應(yīng)充分認(rèn)識(shí)馬鈴薯產(chǎn)量波動(dòng)規(guī)律，采取適當(dāng)?shù)纳a(chǎn)引導(dǎo)措施，在馬鈴薯產(chǎn)量的“低谷期”適當(dāng)采取相應(yīng)措施刺激生產(chǎn)，如馬鈴薯生產(chǎn)者補(bǔ)貼政策、種薯補(bǔ)貼和良種推廣措施。而在馬鈴薯產(chǎn)量波動(dòng)的“高峰期”，為避免出現(xiàn)“薯賤傷農(nóng)”的局面，最大限度平抑薯價(jià)下跌，政府應(yīng)采取一定措施，減少市場中馬鈴薯流通量，如加大政府政策性采購，對馬鈴薯種植戶建冷藏庫入窖給予一定的補(bǔ)貼資金[25]。

2)提高馬鈴薯生產(chǎn)效率水平，提升馬鈴薯種植現(xiàn)代化水平。相比于其他發(fā)達(dá)國家，我國馬鈴薯產(chǎn)業(yè)大而不強(qiáng)，第一約束即為生產(chǎn)效率水平。對于我們這樣一個(gè)耕地資源有限的國家，依賴土地密集投入發(fā)展的產(chǎn)業(yè)是不符合比較優(yōu)勢的，以此來維持的產(chǎn)量增長和產(chǎn)業(yè)發(fā)展也是不可持續(xù)的[26]。因此，應(yīng)不斷提升我國馬鈴薯生產(chǎn)技術(shù)效率，提高單產(chǎn)水平。通過構(gòu)建靈活有效的脫毒種薯推廣應(yīng)用體系，提升病蟲害防治水平，促進(jìn)馬鈴薯生產(chǎn)的適度規(guī)模化，引進(jìn)先進(jìn)生產(chǎn)力要素來推動(dòng)我國馬鈴薯生產(chǎn)效率水平的整體提升。

3)增強(qiáng)抵御自然因素沖擊能力，加強(qiáng)主產(chǎn)區(qū)自然氣象監(jiān)測預(yù)警。自然條件是馬鈴薯生產(chǎn)波動(dòng)的主要影響因素，作為主要分布在中西部地區(qū)的農(nóng)作物，降雨、光照和自然災(zāi)害對其生產(chǎn)有顯著影響。通過各級(jí)農(nóng)業(yè)主管部門與氣象、應(yīng)急管理局等部門的溝通協(xié)調(diào)，建立聯(lián)動(dòng)工作體系[27]。充分利用氣象監(jiān)測部門的預(yù)警系統(tǒng)，對馬鈴薯主產(chǎn)區(qū)(種植區(qū))提供及時(shí)有效的氣象、災(zāi)害預(yù)警預(yù)報(bào)，增強(qiáng)種植戶抵御自然氣象災(zāi)害沖擊的能力。

4)提升消費(fèi)能力，拓展多種馬鈴薯產(chǎn)品消費(fèi)途徑。受飲食習(xí)慣影響，目前相比于國外，國內(nèi)馬鈴薯消費(fèi)的數(shù)量和方式還需要進(jìn)一步提升[28]。應(yīng)充分發(fā)揮政府和市場的消費(fèi)引導(dǎo)能力，加大對馬鈴薯主食化產(chǎn)品的宣傳，引導(dǎo)人們逐漸接受馬鈴薯產(chǎn)品[29]。食品科研機(jī)構(gòu)和企業(yè)應(yīng)加大馬鈴薯產(chǎn)品開發(fā)創(chuàng)新，不斷豐富馬鈴薯產(chǎn)品種類，生產(chǎn)出更多符合我國居民消費(fèi)習(xí)慣的馬鈴薯產(chǎn)品。