聶 彎，謝彥明

（1.中國(guó)社會(huì)科學(xué)院研究生院，北京 102488；2.西南林業(yè)大學(xué)，云南 昆明 650224）

農(nóng)作物病蟲(chóng)害是我國(guó)農(nóng)業(yè)面臨的一項(xiàng)重要難題。受氣候變化、農(nóng)藥使用強(qiáng)度、生物多樣性變化等因素的影響，近年來(lái)，我國(guó)農(nóng)作物病蟲(chóng)害發(fā)生呈現(xiàn)出兩大特點(diǎn)：一是發(fā)生面積逐年增加，由2003年的30 000萬(wàn)公頃次上升到2013年的60 800萬(wàn)公頃次，并且每年還在以8.85%的速度增長(zhǎng)；二是病蟲(chóng)害的種類逐年增加，有關(guān)研究表明［1］，20世紀(jì)50～70年代，每年約有10種農(nóng)作物病蟲(chóng)害暴發(fā)，80～90年代每年約有15種農(nóng)作物病蟲(chóng)害暴發(fā)，21世紀(jì)以來(lái)農(nóng)作物病蟲(chóng)害暴發(fā)的品種增加到30 種左右。農(nóng)作物病蟲(chóng)害發(fā)生的不斷增加給農(nóng)業(yè)生產(chǎn)造成很大威脅，每年因農(nóng)作物病蟲(chóng)害導(dǎo)致的糧食作物產(chǎn)量損失數(shù)以億計(jì)，有關(guān)資料顯示［2］，2006年，僅水稻、小麥、玉米和大豆4種主要糧食作物的實(shí)際產(chǎn)量損失就高達(dá)127.08億kg，是2000年損失的1.21倍。迄今為止，無(wú)論是水稻、小麥、玉米、大豆等主要糧食作物還是蔬菜、水果等園藝作物的病蟲(chóng)害的發(fā)生都呈現(xiàn)出逐年增加的態(tài)勢(shì)。由防治病蟲(chóng)害導(dǎo)致的生態(tài)環(huán)境問(wèn)題也日益凸顯出來(lái)，典型的就是農(nóng)藥的過(guò)量使用帶來(lái)的一系列問(wèn)題。長(zhǎng)期以來(lái)，通過(guò)農(nóng)藥防治病蟲(chóng)害的發(fā)生是我國(guó)最主要的病蟲(chóng)害防治技術(shù)。農(nóng)藥是一種化學(xué)試劑，具有殺死病蟲(chóng)害的作用，但由于缺乏對(duì)農(nóng)藥防治農(nóng)作物病蟲(chóng)害的準(zhǔn)確預(yù)測(cè)預(yù)報(bào)，導(dǎo)致我國(guó)農(nóng)藥利用率不足10%，并且在不同區(qū)域存在較大的效果差異。一方面，由于農(nóng)藥利用率極低，導(dǎo)致農(nóng)民盲目通過(guò)加大農(nóng)藥使用強(qiáng)度以期控制農(nóng)作物病蟲(chóng)害的發(fā)生，但其對(duì)農(nóng)作物病蟲(chóng)害發(fā)生面積的控制作用到底如何，是否存在邊際遞減效用，有待進(jìn)行深入研究；另一方面，農(nóng)藥殘留通過(guò)食物鏈進(jìn)入人體，對(duì)人類的生存和發(fā)展產(chǎn)生嚴(yán)重的威脅，亟須開(kāi)展相應(yīng)研究。

針對(duì)農(nóng)藥使用強(qiáng)度對(duì)農(nóng)作物病蟲(chóng)害發(fā)生面積的影響研究，國(guó)外學(xué)者研究較早，并且取得了豐碩的研究成果。Sivaprakasam等［3］早在1981年就植物保護(hù)對(duì)高粱病蟲(chóng)害發(fā)生的復(fù)雜影響進(jìn)行研究，結(jié)果表明，植物保護(hù)可以減輕高粱病蟲(chóng)害的發(fā)生，進(jìn)而促進(jìn)農(nóng)作物的增產(chǎn)。Hardwick等［4］在1996年關(guān)于氣候變化對(duì)英國(guó)主要農(nóng)作物病蟲(chóng)害的發(fā)病率、嚴(yán)重程度和分布的影響進(jìn)行了研究，結(jié)果表明，氣候變化使得本土農(nóng)作物病蟲(chóng)害攻擊性的嚴(yán)重程度增加。1997年，Chakraborty等 就氣候變化對(duì)澳大利亞農(nóng)作物病蟲(chóng)害發(fā)生的影響進(jìn)行了全面而系統(tǒng)的研究，結(jié)果表明，燃燒化石燃料、大規(guī)模砍伐森林和其他人類活動(dòng)已經(jīng)導(dǎo)致全球變暖，并且對(duì)澳大利亞植物病蟲(chóng)害的發(fā)生和地理分布產(chǎn)生影響。此后，Knight等［6］、Massaad 等［7］就氣候變化對(duì)農(nóng)作物病蟲(chóng)害發(fā)生的影響進(jìn)行了深入研究，結(jié)果表明，氣候變化，如降雨短缺、炎熱、干旱、白天與黑夜之間的溫度波動(dòng)以及極端天氣對(duì)農(nóng)作物病蟲(chóng)害的種類和地理傳播等具有重要影響。近年來(lái)，Abudulai等［8］和 Danso等［9］從農(nóng)作物管理方式方面，探究它們對(duì)農(nóng)作物病蟲(chóng)害發(fā)生面積的影響。國(guó)內(nèi)關(guān)于農(nóng)作物病蟲(chóng)害發(fā)生面積的影響因素的研究起步較晚?；糁螄?guó)等［1］從2000年開(kāi)始就氣象因素對(duì)農(nóng)作物病蟲(chóng)害的影響進(jìn)行研究，并且就降水減少、暴雨、洪澇、高溫、干旱、等氣象因素對(duì)農(nóng)作物病蟲(chóng)害的影響進(jìn)行了深入研究，結(jié)果表明，一定區(qū)域和時(shí)段的降水、高溫、干旱、臺(tái)風(fēng)等多種氣候變化因素均會(huì)對(duì)農(nóng)作物病蟲(chóng)害的發(fā)生產(chǎn)生影響。趙淼等［10］、尹朝靜等［11］針對(duì)影響農(nóng)作物病蟲(chóng)害發(fā)生的氣候變化因素進(jìn)行了深入研究。目前，國(guó)內(nèi)有關(guān)影響農(nóng)作物病蟲(chóng)害發(fā)生的氣象因素的研究已經(jīng)非常成熟，研究成果也相當(dāng)豐富。就農(nóng)藥使用強(qiáng)度對(duì)農(nóng)作物病蟲(chóng)害發(fā)生面積的影響的相關(guān)研究較多，徐曉鵬［12］研究了農(nóng)戶施藥行為的變遷過(guò)程。郭利京等［13］對(duì)農(nóng)戶生物農(nóng)藥施用意愿進(jìn)行了研究。應(yīng)瑞瑤等［14］等反向分析了農(nóng)作物病蟲(chóng)害防治對(duì)農(nóng)藥施用強(qiáng)度的影響。

綜合以上國(guó)內(nèi)外關(guān)于農(nóng)作物病蟲(chóng)害發(fā)生面積影響因素的研究以及就農(nóng)藥使用強(qiáng)度對(duì)農(nóng)作物病蟲(chóng)害發(fā)生面積的影響的相關(guān)研究可以看出，大部分研究都集中在探究氣候變化對(duì)農(nóng)作物病蟲(chóng)害發(fā)生面積的影響、評(píng)價(jià)農(nóng)藥使用強(qiáng)度以及農(nóng)藥對(duì)農(nóng)作物病蟲(chóng)害的作用機(jī)理和技術(shù)改進(jìn)方面，針對(duì)農(nóng)藥使用強(qiáng)度對(duì)農(nóng)作物病蟲(chóng)害發(fā)生面積的影響的研究尚不多見(jiàn)。在研究方法方面，以往的研究主要采用定性的描述方法，缺乏定量分析，即使運(yùn)用了計(jì)量分析方法，但是限于當(dāng)時(shí)計(jì)量經(jīng)濟(jì)學(xué)發(fā)展尚不成熟，計(jì)量模型也多是簡(jiǎn)單的普通最小二乘估計(jì)。在此背景下，我們基于已有研究成果，以邊際效用遞減理論為指導(dǎo)，采用面板數(shù)據(jù)隨機(jī)效應(yīng)模型就農(nóng)藥使用強(qiáng)度對(duì)農(nóng)作物病蟲(chóng)害發(fā)生面積的影響進(jìn)行實(shí)證研究，并且提出針對(duì)性的對(duì)策建議，以期為我國(guó)防治農(nóng)作物病蟲(chóng)害提供理論依據(jù)。

1 模型設(shè)定與變量說(shuō)明

1.1 模型設(shè)定

靜態(tài)面板數(shù)據(jù)模型主要有混合估計(jì)模型、固定效應(yīng)模型和隨機(jī)效應(yīng)模型3種。首先通過(guò)F檢驗(yàn)在混合估計(jì)模型和固定效應(yīng)模型中選擇了固定效應(yīng)模型，然后，通過(guò)LM檢驗(yàn)在混合估計(jì)模型和隨機(jī)效應(yīng)模型中選擇了隨機(jī)效應(yīng)模型，最后，通過(guò)Hausman檢驗(yàn)在固定效應(yīng)模型和隨機(jī)效應(yīng)模型中選擇了隨機(jī)效應(yīng)模型，具體的模型設(shè)定如下：

式中，i代表30個(gè)?。ㄊ小⒆灾螀^(qū)），t（t =2003，…2013）代表年份，lnnyqd代表農(nóng)藥使用強(qiáng)度的對(duì)數(shù)值，dyx代表農(nóng)作物多樣性指數(shù)，lnsr代表農(nóng)村居民家庭人均純收入的對(duì)數(shù)值，β1、β2、β3分別為農(nóng)藥使用強(qiáng)度、作物多樣性和農(nóng)村居民家庭人均純收入的系數(shù)，εit為隨機(jī)擾動(dòng)項(xiàng)。

1.2 變量說(shuō)明

1.2.1 因變量 選取農(nóng)作物病蟲(chóng)害發(fā)生面積（fsl）作為因變量，并且采用當(dāng)年農(nóng)作物病蟲(chóng)害發(fā)生面積的對(duì)數(shù)作為農(nóng)作物病蟲(chóng)害發(fā)生面積的衡量指標(biāo)。

1.2.2 自變量 選取農(nóng)藥使用強(qiáng)度（nyqd）作為自變量，并且用當(dāng)年農(nóng)作物農(nóng)藥使用量除以當(dāng)年農(nóng)作物播種面積來(lái)衡量，這種方法得到了向濤［15］等多數(shù)學(xué)者的支持。

1.2.3 控制變量 選取農(nóng)作物多樣性（dyx）和農(nóng)村居民家庭人均純收入（sr）作為控制變量。其中，農(nóng)作物多樣性指數(shù)采用Shannon信息指數(shù)計(jì)算得到，計(jì)算公式為H = -∑pilnpi，式中，H為多樣性指數(shù)，pi代表第i種作物的種植面積比例，為確保形式上的一致，規(guī)定ln = 0，這種做法得到了張蕾等［14］多數(shù)學(xué)者的支持。作物多樣性指數(shù)越大，說(shuō)明該?。ㄊ?、自治區(qū)）的農(nóng)作物種植的種類越多，并且各種作物種植面積分布越均勻，反之則說(shuō)明該省（市、自治區(qū)）農(nóng)作物種植的種類越少，各種農(nóng)作物種植面積分布越不均勻。農(nóng)村居民家庭人均純收入采用當(dāng)年農(nóng)村居民家庭人均純收入的對(duì)數(shù)來(lái)衡量。

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

數(shù)據(jù)均來(lái)自歷年《中國(guó)農(nóng)業(yè)年鑒》和《中國(guó)農(nóng)村統(tǒng)計(jì)年鑒》。由于西藏地區(qū)的數(shù)據(jù)缺失嚴(yán)重，故將其剔除，最終得到2003—2013年全國(guó)30個(gè)?。ㄊ?、自治區(qū)）的330個(gè)樣本作為研究樣本。

2 變量特征描述

2.1 變量的描述性統(tǒng)計(jì)

從表1可以看出，近10多年來(lái)，我國(guó)農(nóng)作物病蟲(chóng)害發(fā)生面積的對(duì)數(shù)的均值較大，表明我國(guó)農(nóng)作物病蟲(chóng)害較嚴(yán)重，農(nóng)作物多樣性的對(duì)數(shù)的均值較小，表明我國(guó)農(nóng)作物多樣性較低。

表1 全部樣本變量的描述性統(tǒng)計(jì)

表2 東、中、西部地區(qū)變量的描述性統(tǒng)計(jì)

從表2可以看出，東、中、西部地區(qū)的農(nóng)作物病蟲(chóng)害發(fā)生面積、農(nóng)藥使用強(qiáng)度、作物多樣性水平、農(nóng)村居民家庭人均純收入的描述性統(tǒng)計(jì)結(jié)果支持全部樣本的描述性統(tǒng)計(jì)結(jié)果，但是各地區(qū)之間存在差異。其中，中部地區(qū)農(nóng)作物病蟲(chóng)害發(fā)生面積的均值最大、標(biāo)準(zhǔn)差最小，表明中部地區(qū)農(nóng)作物病蟲(chóng)害發(fā)生最嚴(yán)重；東部地區(qū)和中部地區(qū)農(nóng)藥使用強(qiáng)度均值均大于西部地區(qū)農(nóng)藥使用強(qiáng)度，并且方差均小于西部地區(qū)，說(shuō)明東部地區(qū)和中部地區(qū)農(nóng)藥使用強(qiáng)度普遍高于西部地區(qū)農(nóng)藥使用強(qiáng)度；西部地區(qū)作物多樣性均值最大，但是標(biāo)準(zhǔn)差即不是最小也不是最大，很難從總體上判斷區(qū)域差異性；東部地區(qū)農(nóng)村居民家庭人均純收入均值最大，但方差也最大，說(shuō)明東部地區(qū)農(nóng)村居民家庭人均純收入的狀況并不是十分理想。

2.2 變量的地區(qū)差異分析

2.2.1 各地區(qū)農(nóng)作物病蟲(chóng)害發(fā)生面積的差異 2013年，全國(guó)農(nóng)作物病蟲(chóng)害發(fā)生面積最嚴(yán)重的省份為山東省，高達(dá)5 587萬(wàn)公頃次，是農(nóng)作物病蟲(chóng)害發(fā)生面積最輕的北京市的38倍多。河南省居第2位，高達(dá)5 267.5萬(wàn)公頃次，湖南省以5 052.4萬(wàn)公頃次的高農(nóng)作物病蟲(chóng)害發(fā)生面積高居第3位。青海是西部地區(qū)當(dāng)中農(nóng)作物病蟲(chóng)害發(fā)生面積最小，也是全國(guó)農(nóng)作物病蟲(chóng)害發(fā)生面積僅高于北京市的地區(qū)，其農(nóng)作物病蟲(chóng)害發(fā)生面積為159.5萬(wàn)公頃次。包括天津、上海、海南、寧夏、福建、重慶和貴州地區(qū)的農(nóng)作物病蟲(chóng)害發(fā)生面積均未超過(guò)1 000萬(wàn)公頃次。總體而言，各地區(qū)農(nóng)作物病蟲(chóng)害發(fā)生面積的差異較大，并且具有區(qū)域集聚的特征，表現(xiàn)在中部地區(qū)的農(nóng)作物病蟲(chóng)害發(fā)生面積較東部地區(qū)和西部地區(qū)更甚，支持描述性統(tǒng)計(jì)的結(jié)果（圖1）。

圖1 2013年各地區(qū)農(nóng)作物病蟲(chóng)害發(fā)生面積比較

2.2.2 各地區(qū)農(nóng)作物農(nóng)藥使用強(qiáng)度的差異 2013年，全國(guó)農(nóng)藥使用強(qiáng)度最高的地區(qū)是海南省，高達(dá)51.26 kg/hm2，接近第2名浙江省農(nóng)藥使用強(qiáng)度的兩倍，是農(nóng)藥使用強(qiáng)度最輕的寧夏的24倍多。福建省農(nóng)藥使用強(qiáng)度高居全國(guó)第3位，為25.22 kg/hm2，廣東省農(nóng)藥使用強(qiáng)度為23.43 kg/hm2，居第4位。包括甘肅、江西、北京、江蘇等11個(gè)?。ㄊ校┑霓r(nóng)藥使用強(qiáng)度在10～20 kg/hm2，其他省（市、自治區(qū)）農(nóng)藥使用強(qiáng)度均在10 kg/hm2以下?？傮w而言，各地區(qū)農(nóng)作物農(nóng)藥使用強(qiáng)度的差異較大，并且具有明顯的地域集中分布特征，其中，東部和中部地區(qū)的農(nóng)藥使用強(qiáng)度普遍較高，而西部地區(qū)的農(nóng)藥使用強(qiáng)度普遍較低，與描述性統(tǒng)計(jì)結(jié)果一致（圖2）。

2.2.3 各地區(qū)農(nóng)作物多樣性的差異 2013年，全國(guó)農(nóng)作物多樣性指數(shù)最高的地區(qū)是新疆維吾爾自治區(qū)，其生物多樣性指數(shù)為1.39，湖北省農(nóng)作物多樣性指數(shù)高居第2位，為1.38，廣西壯族自治區(qū)居第3位，為1.36。黑龍江省的農(nóng)作物多樣性指數(shù)最低，僅為0.28，吉林的農(nóng)作物多樣性指數(shù)僅高于黑龍江，為0.49?？傮w而言，各地區(qū)農(nóng)作物多樣性指數(shù)較低，地區(qū)差異較小，但是沒(méi)有明顯的地域集中分布特征，與描述性統(tǒng)計(jì)結(jié)果一致（圖3）。

圖2 2013年各地區(qū)農(nóng)作物農(nóng)藥使用強(qiáng)度比較

圖3 2013年各地區(qū)農(nóng)作物多樣性指數(shù)比較

2.2.4 各地區(qū)農(nóng)村居民家庭人均純收入的差異 2013年，全國(guó)農(nóng)村居民家庭人均純收入最高的地區(qū)為上海，高達(dá)19 595元，北京僅次于上海，農(nóng)村居民家庭人均純收入為18 337.5元。包括浙江、天津、江蘇、廣東、福建、山東和遼寧地區(qū)的農(nóng)村居民家庭人均純收入均超過(guò)了10 000元。甘肅省農(nóng)村居民家庭人均純收入居全國(guó)最后1位，低至5 107.8元，貴州省居倒數(shù)第2位，農(nóng)村居民家庭人均純收入為5 434元。其他?。ㄊ小⒆灾螀^(qū)）的農(nóng)村居民家庭人均純收入均在6 000～10 000元之間?？傮w而言，全國(guó)農(nóng)村居民家庭人均純收入地區(qū)差異明顯，并且具有區(qū)域集中分布的特征，東部地區(qū)遠(yuǎn)高于西部地區(qū)（圖4）。

3 面板模型檢驗(yàn)及隨機(jī)效應(yīng)回歸結(jié)果分析

3.1 面板單位根檢驗(yàn)

為了避免可能出現(xiàn)的虛假回歸現(xiàn)象，需要對(duì)面板數(shù)據(jù)進(jìn)行單位根檢驗(yàn)。面板數(shù)據(jù)的單位根檢驗(yàn)方法主要有LLC、IPS、ADF-F等，為了提高檢驗(yàn)的效力，運(yùn)用這3種檢驗(yàn)對(duì)全部樣本和東部地區(qū)、中部地區(qū)以及西部地區(qū)的樣本均進(jìn)行單位根檢驗(yàn)，并且對(duì)檢驗(yàn)回歸項(xiàng)中有趨勢(shì)項(xiàng)和無(wú)趨勢(shì)項(xiàng)分別進(jìn)行檢驗(yàn)，結(jié)果見(jiàn)表3。

圖4 2013年各地區(qū)農(nóng)村居民家庭人均純收入比較

表3 面板單位根檢驗(yàn)結(jié)果

如表3所示，全部樣本、東部地區(qū)樣本、中部地區(qū)樣本和西部地區(qū)樣本的農(nóng)作物病蟲(chóng)害發(fā)生面積的對(duì)數(shù)、農(nóng)藥使用強(qiáng)度的對(duì)數(shù)、作物多樣性的水平值均通過(guò)了單位根檢驗(yàn)，是平穩(wěn)序列，農(nóng)村居民家庭人均純收入的對(duì)數(shù)的水平值無(wú)法拒絕單位根檢驗(yàn)，為非平穩(wěn)序列，進(jìn)行一階差分后在1%的顯著性水平下拒絕了3種單位根檢驗(yàn)，表明農(nóng)村居民家庭人均純收入的一階差分為平穩(wěn)過(guò)程。

3.2 全部樣本的隨機(jī)效應(yīng)模型回歸結(jié)果分析

從表4模型設(shè)定后的檢驗(yàn)結(jié)果可以看出，固定效應(yīng)模型優(yōu)于混合估計(jì)模型，而Hausman檢驗(yàn)則表明隨機(jī)效應(yīng)模型優(yōu)于固定效應(yīng)模型，因此對(duì)隨機(jī)效應(yīng)回歸模型的回歸結(jié)果進(jìn)行分析和說(shuō)明。

表4 全部樣本3種面板模型回歸結(jié)果

農(nóng)藥使用強(qiáng)度對(duì)農(nóng)作物病蟲(chóng)害發(fā)生面積的影響顯著為正，說(shuō)明農(nóng)藥使用強(qiáng)度越大，農(nóng)作物病蟲(chóng)害的發(fā)生面積越大。農(nóng)藥使用強(qiáng)度每增加1個(gè)單位會(huì)顯著增加農(nóng)作物病蟲(chóng)害發(fā)生面積0.2308個(gè)單位。近幾年，我國(guó)農(nóng)作物病蟲(chóng)害的發(fā)生面積逐年加重，通過(guò)加大農(nóng)藥使用強(qiáng)度來(lái)治理和預(yù)防農(nóng)作物病蟲(chóng)害的需求也日益迫切，進(jìn)而導(dǎo)致我國(guó)農(nóng)作物農(nóng)藥使用嚴(yán)重超標(biāo)，不僅對(duì)農(nóng)業(yè)生態(tài)環(huán)境造成嚴(yán)重的壓力，也對(duì)人類健康產(chǎn)生巨大威脅。但是，與傳統(tǒng)認(rèn)識(shí)相反，農(nóng)藥使用強(qiáng)度的一味增加并不能有效抑制農(nóng)作物病蟲(chóng)害的發(fā)生，反而會(huì)促進(jìn)農(nóng)作物病蟲(chóng)害的發(fā)生。究其原因，是農(nóng)作物病蟲(chóng)害對(duì)農(nóng)藥使用會(huì)產(chǎn)生適應(yīng)性反應(yīng)，并且通過(guò)產(chǎn)生抗體等多種方式來(lái)抵抗農(nóng)藥的殺害。隨著農(nóng)藥使用強(qiáng)度的增加，這種適應(yīng)性反應(yīng)非但不會(huì)減弱，反而會(huì)增強(qiáng)，進(jìn)而導(dǎo)致農(nóng)作物病蟲(chóng)害的發(fā)生面積隨農(nóng)藥使用強(qiáng)度的增加而增加，這與邊際效用遞減規(guī)律相契合。由此可見(jiàn)，通過(guò)增加農(nóng)藥使用強(qiáng)度來(lái)抵制農(nóng)作物病蟲(chóng)害發(fā)生面積的方法行不通，亟待探索新的防治農(nóng)作物病蟲(chóng)害發(fā)生的辦法。

在控制變量方面，作物多樣性的系數(shù)為負(fù)，說(shuō)明作物多樣性對(duì)農(nóng)作物病蟲(chóng)害發(fā)生面積的影響為負(fù)，作物多樣性的增加可以減少農(nóng)作物病蟲(chóng)害的發(fā)生面積。究其原因，主要是作物多樣性可以保證生態(tài)系統(tǒng)的自動(dòng)調(diào)節(jié)能力，進(jìn)而有效抵制病蟲(chóng)害的發(fā)生和發(fā)展，從而減少農(nóng)作物病蟲(chóng)害的發(fā)生。但是，近幾年來(lái)，由于經(jīng)濟(jì)作物種植面積的持續(xù)增加，而非經(jīng)濟(jì)作物種植面積的持續(xù)減少，甚至有的地區(qū)停止種植非經(jīng)濟(jì)作物，使得我國(guó)作物多樣性銳減，這將助長(zhǎng)農(nóng)作物病蟲(chóng)害的發(fā)生，不利于我國(guó)農(nóng)業(yè)的健康可持續(xù)發(fā)展。農(nóng)村居民家庭人均純收入的系數(shù)為正，但是在混合估計(jì)和固定效應(yīng)模型中，農(nóng)村居民家庭人均純收入的系數(shù)有正有負(fù)，出現(xiàn)這種現(xiàn)象的原因可能是樣本量太小。

3.3 東部、中部和西部地區(qū)隨機(jī)效應(yīng)模型回歸結(jié)果的比較分析

東部地區(qū)和中部地區(qū)農(nóng)藥使用強(qiáng)度對(duì)農(nóng)作物病蟲(chóng)害發(fā)生面積的影響與全部樣本隨機(jī)效應(yīng)回歸結(jié)果一致，顯著為正。在東部地區(qū)，農(nóng)藥使用強(qiáng)度每增加1個(gè)單位會(huì)顯著增加農(nóng)作物病蟲(chóng)害的發(fā)生面積0.3818個(gè)單位；在中部地區(qū)，農(nóng)藥使用強(qiáng)度每增加1個(gè)單位會(huì)顯著增加農(nóng)作物病蟲(chóng)害的發(fā)生面積0.3931個(gè)單位?？梢钥闯觯胁康貐^(qū)每增加1個(gè)單位農(nóng)藥使用強(qiáng)度對(duì)農(nóng)作物病蟲(chóng)害發(fā)生面積的促進(jìn)作用比東部地區(qū)更甚。相比之下，西部地區(qū)農(nóng)藥使用強(qiáng)度對(duì)農(nóng)作物病蟲(chóng)害發(fā)生面積的正向影響表現(xiàn)不顯著，可能的原因是西部地區(qū)各?。ㄊ小⒆灾螀^(qū)）的農(nóng)藥使用強(qiáng)度總體小于東部和中部地區(qū)農(nóng)藥使用強(qiáng)度，并且西部地區(qū)各?。ㄊ?、自治區(qū)）的農(nóng)藥使用強(qiáng)度差異較東部地區(qū)和西部地區(qū)更大（表5）。

表5 東、中、西部地區(qū)隨機(jī)效應(yīng)模型回歸結(jié)果

在控制變量方面，東部地區(qū)和西部地區(qū)作物多樣性對(duì)農(nóng)作物病蟲(chóng)害發(fā)生面積的影響均為負(fù)，與全部樣本隨機(jī)效應(yīng)回歸結(jié)果一致，但是，中部地區(qū)作物多樣性的回歸系數(shù)卻為正，可能原因是中部地區(qū)樣本量太少。農(nóng)村居民家庭人均純收入對(duì)農(nóng)作物病蟲(chóng)害發(fā)生面積的影響顯著為正，與全部樣本隨機(jī)效應(yīng)回歸結(jié)果一致。

4 結(jié)論與建議

通過(guò)對(duì)2003—2013年我國(guó)農(nóng)藥使用強(qiáng)度對(duì)農(nóng)作物病蟲(chóng)害發(fā)生面積的影響進(jìn)行隨機(jī)效應(yīng)回歸分析，并且將全國(guó)（除西藏）劃分為東部地區(qū)、中部地區(qū)和西部地區(qū)分別進(jìn)行隨機(jī)效應(yīng)回歸分析，研究結(jié)果表明，全部樣本、東部地區(qū)樣本和中部地區(qū)樣本均顯著支持農(nóng)藥使用強(qiáng)度的增加非但不能減少農(nóng)作物病蟲(chóng)害發(fā)生面積，反而會(huì)促進(jìn)農(nóng)作物病蟲(chóng)害的發(fā)生；全部樣本、東部地區(qū)樣本、中部地區(qū)樣本以及西部地區(qū)樣本均支持作物多樣性對(duì)農(nóng)作物病蟲(chóng)害發(fā)生面積具有抑制作用。

基于上述認(rèn)識(shí)，為了有效控制農(nóng)作物病蟲(chóng)害的發(fā)生面積，可采用如下三方面措施：一是轉(zhuǎn)變農(nóng)作物病蟲(chóng)害防治方式，鼓勵(lì)農(nóng)民科學(xué)用藥，用好藥，少用藥，并且加大力度支持高效農(nóng)藥的研發(fā)，從效率上提高農(nóng)藥對(duì)農(nóng)作物病蟲(chóng)害發(fā)生的防治作用；二是加大力度保護(hù)農(nóng)作物多樣性，通過(guò)加強(qiáng)農(nóng)業(yè)供給側(cè)結(jié)構(gòu)改革，引導(dǎo)農(nóng)民豐富農(nóng)作物種植結(jié)構(gòu)，調(diào)整農(nóng)業(yè)生產(chǎn)結(jié)構(gòu)為生態(tài)友好、結(jié)構(gòu)合理、健康發(fā)展的生產(chǎn)結(jié)構(gòu)，從農(nóng)作物自身源頭開(kāi)始，防止農(nóng)作物病蟲(chóng)害發(fā)生面積的增加；三是東部地區(qū)和中部地區(qū)應(yīng)盡早減少農(nóng)藥使用強(qiáng)度、加強(qiáng)農(nóng)作物多樣性保護(hù)以防治農(nóng)作物病蟲(chóng)害發(fā)生面積的加重。

參考文獻(xiàn)：

［1］霍治國(guó)，李茂松，王麗，等. 氣候變暖對(duì)中國(guó)農(nóng)作物病蟲(chóng)害的影響［J］. 中國(guó)農(nóng)業(yè)科學(xué)，2012，45（10）：1926-1934.

［2］霍治國(guó)，李茂松，李娜，等. 季節(jié)性變暖對(duì)中國(guó)農(nóng)作物病蟲(chóng)害的影響［J］. 中國(guó)農(nóng)業(yè)科學(xué)，2012，45（11）：2168-2179.

［3］Sivaprakasam K，Balasubramanian G. Impact of plant protection on the pest and disease complex of sorghum［J］. Pestology，1981，27（1）：15-35.

［4］Hardwick N V，Armstrong A C，Ellis S A，et al.The impact of climate change on crop pests and diseases［J］. Aspects of Applied Biology，1996，19（1）：11-17.

［5］Chakraborty S，Murray G M，Magarey P A，et al. Potential impact of climate change on plant diseases of economic significance to Australia［J］. Australasian Plant Pathology，1998，27（1）：15-35.

［6］Knight B E A，Wimshurst A A，Alford D V，et al. Impact of climate change on the geographical spread of agricultural pests，diseases and weeds［C］. Plant Protection & Plant Health in Europe：Introduction & Spread of Invasive Species，Held at Humboldt University，Berlin，Germany，June.2005.

［7］Massaad R，Jomaa I，Hobeika C，et al. Climate change：potential impact on plant diseases，pests and physiology in Lebanon［J］. Historical Materialism，2012，12（3）：315-331.

［8］Abudulai M，Kusi F，Seini S S，et al. Effects of planting date，cultivar and insecticide spray application for the management of insect pests of cowpea in northern Ghana［J］. Crop Protection，2017，100：168-176.

［9］Danso J K，Osekre E A，Manu N，et al. Moisture content，insect pests and mycotoxin levels of maize at harvest and post-harvest in the Middle Belt of Ghana［J］. Journal of Stored Products Research，2017，74：46-55.

［10］趙淼，趙闖等. 近20年來(lái)我國(guó)農(nóng)作物病蟲(chóng)害時(shí)空變化特征［J］. 北京大學(xué)學(xué)報(bào)（自然科學(xué)版），2015 （5）：965-975.

［11］尹朝靜，李谷成，范麗霞，等. 氣候變化、科技存量與農(nóng)業(yè)生產(chǎn)率增長(zhǎng)［J］. 中國(guó)農(nóng)村經(jīng)濟(jì)，2016（5）：16-28.

［12］徐曉鵬. 農(nóng)戶農(nóng)藥施用行為變遷的社會(huì)學(xué)考察—— 基于我國(guó)6省6村的實(shí)證研究［J］. 中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào)（社會(huì)科學(xué)版），2017，34（1）：38-45.

［13］郭利京，趙瑾. 認(rèn)知沖突視角下農(nóng)戶生物農(nóng)藥施用意愿研究—— 基于江蘇639戶稻農(nóng)的實(shí)證［J］. 南京農(nóng)業(yè)大學(xué)學(xué)報(bào)（社會(huì)科學(xué)版），2017，17（2）：123-133，154.

［14］應(yīng)瑞瑤，徐斌. 農(nóng)作物病蟲(chóng)害專業(yè)化防治服務(wù)對(duì)農(nóng)藥施用強(qiáng)度的影響［J］. 中國(guó)人口·資源與環(huán)境，2017，27（8）：90-97.

［15］向濤，綦勇. 糧食安全、食品安全與貿(mào)易—— 基于農(nóng)藥使用強(qiáng)度的跨國(guó)面板數(shù)據(jù)分析［J］. 國(guó)際貿(mào)易問(wèn)題，2014（7）：33-41.