張曉東, 左亞男, 季 民

（山東科技大學(xué)測(cè)繪與空間信息學(xué)院, 山東 青島 266000）

玉米螟（Pyrausta nubilalis）屬鱗翅目（Lepidoptera）螟蛾科（Pyralidae）, 又名玉米鉆心蟲(chóng), 是玉米上的重要害蟲(chóng), 也是世界性害蟲(chóng), 在山東省主要發(fā)生2～3 代, 玉米螟在玉米的各個(gè)生育時(shí)期都可以為害玉米植株的地上部分, 取食葉片、果穗、雄穗。其幼蟲(chóng)孵化后群居于玉米心葉喇叭口處或嫩葉上取食, 葉片被幼蟲(chóng)咬食后, 會(huì)降低其光合效率；雄穗被蛀, 常易折斷, 影響授粉；苞葉、花絲被蛀食, 會(huì)造成缺粒和秕粒；莖稈、穗柄、穗軸被蛀食后, 形成隧道, 破壞植株內(nèi)水分、養(yǎng)分的輸送, 使莖稈倒折率增加, 子粒產(chǎn)量下降。山東省是糧食大省, 也是玉米高產(chǎn)量省份, 玉米產(chǎn)量約占全國(guó)玉米總產(chǎn)量的10%, 為山東省主要種植作物。玉米產(chǎn)量每年因玉米螟都造成不同程度的損失, 因此, 對(duì)山東省玉米螟發(fā)生風(fēng)險(xiǎn)進(jìn)行有效評(píng)估, 對(duì)玉米螟的防治有重要意義。

災(zāi)害風(fēng)險(xiǎn)區(qū)劃是防災(zāi)的重要工具, 在災(zāi)害風(fēng)險(xiǎn)管理和防災(zāi)減災(zāi)方面具有重要的實(shí)際應(yīng)用價(jià)值。賈改秀等[1］在探究甘肅省玉米病蟲(chóng)害風(fēng)險(xiǎn)評(píng)估指標(biāo)的基礎(chǔ)上, 構(gòu)建了甘肅省玉米病蟲(chóng)害風(fēng)險(xiǎn)區(qū)劃指標(biāo), 并進(jìn)行了風(fēng)險(xiǎn)區(qū)劃和分區(qū)評(píng)述。袁福香等[2］整理了吉林省50 個(gè)縣市的主要農(nóng)作物產(chǎn)量面積資料和相關(guān)病蟲(chóng)害資料, 根據(jù)災(zāi)害發(fā)生的頻數(shù)確定病蟲(chóng)害影響指數(shù), 分析各因子的影響程度, 并將各因子極差化處理, 計(jì)算出綜合風(fēng)險(xiǎn)防御指數(shù), 將吉林省主要農(nóng)作物病蟲(chóng)害分成高、中、低風(fēng)險(xiǎn)防御區(qū)進(jìn)行區(qū)劃評(píng)估。李陽(yáng)等[3］利用彭陽(yáng)縣的氣象、地形、土壤類型等數(shù)據(jù), 基于災(zāi)害風(fēng)險(xiǎn)評(píng)估原理, 運(yùn)用相關(guān)分析法、層次分析法, 結(jié)合GIS 技術(shù)對(duì)彭陽(yáng)縣玉米干旱災(zāi)害進(jìn)行風(fēng)險(xiǎn)評(píng)估與區(qū)劃。Singh 等[4］在進(jìn)行農(nóng)業(yè)區(qū)的風(fēng)險(xiǎn)區(qū)劃時(shí), 從影響農(nóng)業(yè)區(qū)抗災(zāi)能力的4 個(gè)方面考慮選取評(píng)價(jià)指標(biāo), 其多維度、多指標(biāo)的討論過(guò)程為本研究評(píng)價(jià)指標(biāo)的分析與選取提供了思路。

目前國(guó)內(nèi)外針對(duì)風(fēng)險(xiǎn)區(qū)劃[1-15］的研究已相對(duì)成熟, 但仍然缺少生物災(zāi)害方向的區(qū)劃研究, 因此需對(duì)玉米螟災(zāi)害的發(fā)生進(jìn)行風(fēng)險(xiǎn)區(qū)劃, 以期為生物災(zāi)害的防治及預(yù)測(cè)提供科學(xué)支撐。

對(duì)研究區(qū)域內(nèi)玉米螟的發(fā)生資料進(jìn)行數(shù)據(jù)提取, 并將玉米螟發(fā)生程度分級(jí)國(guó)家農(nóng)業(yè)行業(yè)標(biāo)準(zhǔn)（NY/T 1611—2008）與山東省植物保護(hù)總站地方標(biāo)準(zhǔn)相結(jié)合, 提出了新的玉米螟發(fā)生程度分級(jí)準(zhǔn)則, 并基于此準(zhǔn)則進(jìn)行發(fā)生程度定級(jí), 統(tǒng)計(jì)發(fā)生程度頻數(shù)；以玉米螟發(fā)生頻數(shù)、玉米產(chǎn)量和玉米播種面積分別作為致災(zāi)因子、承載體因子和孕災(zāi)環(huán)境因子的影響指標(biāo), 基于參數(shù)估計(jì)法建立風(fēng)險(xiǎn)評(píng)估模型, 構(gòu)建玉米螟災(zāi)害風(fēng)險(xiǎn)評(píng)估指標(biāo)體系, 進(jìn)行了山東省玉米螟災(zāi)害的綜合風(fēng)險(xiǎn)評(píng)估與區(qū)劃, 以期為生物災(zāi)害的防治及預(yù)測(cè)提供科學(xué)支撐。

1 研究區(qū)與數(shù)據(jù)處理方法

1.1 研究區(qū)概況

山東省位于中國(guó)東部沿海、黃河下游, 地處34°25′—38°23′N, 114°36′—112°43′E, 東 西 長(zhǎng) 約700 km, 南北寬約420 km, 陸地總面積約15.67萬(wàn)km2, 是玉米生產(chǎn)和消費(fèi)大省, 玉米面價(jià)格和總產(chǎn)量均居全國(guó)前列, 常年種植面積約2.7×106hm2。山東省屬溫暖帶季風(fēng)氣候, 降水集中, 雨熱同季, 年平均氣溫為11～14 ℃, 年平均降水量為550～950 mm, 且集中于夏季, 光照時(shí)數(shù)年均為2 290～2 890 h, 自然條件非常適合玉米生長(zhǎng)。

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

所用數(shù)據(jù)包括山東省下屬16 個(gè)地市的玉米生產(chǎn)資料及玉米螟病蟲(chóng)害發(fā)生資料。其中, 生產(chǎn)資料來(lái)源于《山東省統(tǒng)計(jì)年鑒》, 包括玉米產(chǎn)量及播種面積資料；病蟲(chóng)害資料來(lái)自山東省農(nóng)業(yè)農(nóng)村廳發(fā)布的植保情況, 補(bǔ)充資料來(lái)源于各期刊文獻(xiàn)、新聞媒體及網(wǎng)站等發(fā)布的病蟲(chóng)情報(bào), 資料年限為2007—2019年。主要數(shù)據(jù)來(lái)源見(jiàn)表1。

表1 主要數(shù)據(jù)來(lái)源

1.3 數(shù)據(jù)預(yù)處理

1.3.1 發(fā)生程度分級(jí)準(zhǔn)則 將收集到的數(shù)據(jù)資料進(jìn)行整合, 由于數(shù)據(jù)資料包含被害株率、蟲(chóng)株率、百株蟲(chóng)量、百株卵量等多項(xiàng)指標(biāo), 且部分月份中已經(jīng)包含個(gè)別市的發(fā)生程度, 經(jīng)過(guò)比對(duì), 發(fā)現(xiàn)其定級(jí)標(biāo)準(zhǔn)不符合中華人民共和國(guó)農(nóng)業(yè)行業(yè)標(biāo)準(zhǔn)——玉米螟測(cè)報(bào)技術(shù)規(guī)范（NY/T 1611—2008）, 因此, 在開(kāi)展定級(jí)工作時(shí), 通過(guò)查閱資料得知, ①山東省植物保護(hù)總站的玉米螟發(fā)生程度分級(jí)地方標(biāo)準(zhǔn)；②在農(nóng)作物蟲(chóng)害方面, 被害株率、蟲(chóng)株率、白穗率及卷葉率等被統(tǒng)稱為被害率, 其計(jì)算公式是統(tǒng)一的, 因此在定級(jí)工作中按照統(tǒng)一標(biāo)準(zhǔn)進(jìn)行定級(jí)。為統(tǒng)一定級(jí)指標(biāo), 根據(jù)NY/T 1611—2008, 并與山東省植物保護(hù)總站農(nóng)作物病蟲(chóng)發(fā)生程度分級(jí)地方標(biāo)準(zhǔn)相結(jié)合, 制定了以被害株率為主要?jiǎng)澐忠蜃蛹鞍僦晗x(chóng)量為輔助因子的玉米螟發(fā)生程度分級(jí)準(zhǔn)則（表2）。

表2 玉米螟發(fā)生程度分級(jí)

1.3.2 發(fā)生程度的定級(jí) 將上述數(shù)據(jù)根據(jù)表2 進(jìn)行發(fā)生程度定級(jí), 選用2009 年9 月的數(shù)據(jù)（資料來(lái)自《山東農(nóng)藥信息》）。首先, 將蟲(chóng)株率按照被害株率的分級(jí)指標(biāo)進(jìn)行發(fā)生程度定級(jí), 且當(dāng)數(shù)據(jù)橫跨兩個(gè)級(jí)別時(shí), 將較高級(jí)別定為其發(fā)生程度, 如2009 年9 月煙臺(tái)市玉米螟被害株率為10%～27%, 則其發(fā)生程度被定為2 級(jí)。其中沒(méi)有數(shù)據(jù)的地市定為未發(fā)生（棗莊市、臨沂市、日照市）, 其他指標(biāo)按照表2 進(jìn)行定級(jí), 定級(jí)工作完成后按表3 所示統(tǒng)計(jì)玉米暝發(fā)生情況。

表3 2009 年9 月山東省各地市玉米螟發(fā)生情況統(tǒng)計(jì)

1.3.3 發(fā)生程度頻數(shù) 定級(jí)工作完成后, 分別統(tǒng)計(jì)發(fā)生程度為1、2、3 級(jí)的頻數(shù)V1、V2、V3, 得到2007—2019 年玉米螟發(fā)生頻數(shù)如表4 所示。

2 風(fēng)險(xiǎn)評(píng)估指標(biāo)的確定及計(jì)算

2.1 風(fēng)險(xiǎn)評(píng)估指標(biāo)的確定

風(fēng)險(xiǎn)評(píng)估指標(biāo)體系的建立是進(jìn)行病蟲(chóng)害風(fēng)險(xiǎn)評(píng)估的關(guān)鍵環(huán)節(jié), 風(fēng)險(xiǎn)指標(biāo)的選取應(yīng)遵循科學(xué)性、實(shí)用性、系統(tǒng)性、人本性原則。依照上述原則, 通過(guò)查閱有關(guān)病蟲(chóng)害研究文獻(xiàn)和歷史資料、危害程度、危害頻率, 結(jié)合病蟲(chóng)害的特點(diǎn), 提出病蟲(chóng)害風(fēng)險(xiǎn)評(píng)估指標(biāo)體系[1, 2, 4-9］, 包括致災(zāi)因子、承載體因子及孕災(zāi)環(huán)境因子。從這3 個(gè)方面闡述其對(duì)災(zāi)害的影響, 選擇具有代表性且可以量化的影響因子作為風(fēng)險(xiǎn)區(qū)劃指標(biāo)。

2.1.1 致災(zāi)因子 致災(zāi)因子是指造成糧食作物損失的各種病蟲(chóng)害, 是重要的風(fēng)險(xiǎn)評(píng)估指標(biāo)。本研究選用玉米螟2007—2019 年不同危害程度（輕發(fā)生、中發(fā)生、重發(fā)生）的發(fā)生頻數(shù)（表4）作為致災(zāi)因子的影響指標(biāo)。

表4 2007—2019 年山東省玉米螟發(fā)生頻數(shù)（V）（單位：次）

2.1.2 承載體因子 承載體因子是指當(dāng)病蟲(chóng)害發(fā)生時(shí)農(nóng)作物受病蟲(chóng)害影響的因素, 是致災(zāi)因子作用的實(shí)體。本研究中, 受玉米螟災(zāi)害影響的承載體是玉米的產(chǎn)量。

2.1.3 孕災(zāi)環(huán)境因子 孕災(zāi)環(huán)境是指影響各種病蟲(chóng)害生長(zhǎng)、繁殖的外界環(huán)境條件。本研究所考慮的影響指標(biāo)包括作物的種植面積, 還有導(dǎo)致病蟲(chóng)害發(fā)生的氣象條件, 氣象條件的適宜與否致使病蟲(chóng)害發(fā)生的程度不同, 此指標(biāo)已包含在致災(zāi)因子的影響中, 不再重復(fù)計(jì)算, 因此本研究只選用玉米播種面積作為病蟲(chóng)害風(fēng)險(xiǎn)評(píng)估指標(biāo)。

綜合以上分析, 所建立的玉米螟災(zāi)害風(fēng)險(xiǎn)評(píng)估指標(biāo)體系如圖1 所示。

圖1 玉米螟風(fēng)險(xiǎn)評(píng)估指標(biāo)體系

2.2 風(fēng)險(xiǎn)評(píng)估指標(biāo)的計(jì)算

2.2.1 玉米螟風(fēng)險(xiǎn)指數(shù) 病蟲(chóng)害影響指數(shù)是衡量病蟲(chóng)害對(duì)農(nóng)作物產(chǎn)量損失風(fēng)險(xiǎn)大小的指標(biāo)。利用山東省16 個(gè)地市2007—2019 年玉米螟危害資料, 根據(jù)玉米螟歷年不同發(fā)生程度（1、2、3 級(jí)）的頻數(shù)（V1、V2、V3）（表3）, 利用加權(quán)平均得出玉米螟對(duì)產(chǎn)量的影響程度, 即病蟲(chóng)害影響指數(shù), 也稱作玉米螟風(fēng)險(xiǎn)指數(shù)影響因子。重發(fā)生的產(chǎn)量損失大, 則權(quán)重系數(shù)大, 輕發(fā)生的產(chǎn)量損失小, 則權(quán)重系數(shù)小。風(fēng)險(xiǎn)指數(shù)（R）的計(jì)算公式如下。

式中,R是玉米螟風(fēng)險(xiǎn)指數(shù)；i是各地市（下同）。

根據(jù)表3 及式（1）計(jì)算出R, 根據(jù)各地市R, 利用反距離權(quán)重法（IDW）、掩膜提取及自然斷點(diǎn)分級(jí)法[11］將玉米螟影響指數(shù)分5 個(gè)等級(jí)進(jìn)行區(qū)劃評(píng)估, 得到山東省玉米螟影響指數(shù)區(qū)劃（圖2）。由圖2 可知, 玉米螟災(zāi)害發(fā)生的主要風(fēng)險(xiǎn)地集中在山東省的西部地區(qū), 該地區(qū)是山東省玉米主要種植區(qū)和高產(chǎn)區(qū), 同時(shí)也是玉米螟發(fā)生的高風(fēng)險(xiǎn)區(qū)。其中, 菏澤市、濱州市風(fēng)險(xiǎn)指數(shù)高, 為玉米螟發(fā)生的高風(fēng)險(xiǎn)地區(qū)；濟(jì)南市中西部為較高風(fēng)險(xiǎn)區(qū)。

圖2 玉米螟影響指數(shù)區(qū)劃

2.2.2 玉米產(chǎn)量影響因子 玉米螟發(fā)生區(qū)域產(chǎn)量越高, 損失越大, 則玉米螟危害風(fēng)險(xiǎn)越大, 因此本研究選擇玉米產(chǎn)量作為承載體因子的影響指標(biāo)。產(chǎn)量影響因子的取值為各地市產(chǎn)量占全省總產(chǎn)量的比重, 產(chǎn)量值為近10 年平均值。玉米產(chǎn)量影響因子（K）的計(jì)算公式如下。

式中,Ky為各地市的玉米總產(chǎn)量；∑Ky為全省的玉米總產(chǎn)量。

根據(jù)式（2）計(jì)算出K, 并繪制玉米產(chǎn)量影響區(qū)劃（圖3）。由圖3可知, 山東省西部地區(qū)及濰坊市中部為玉米高產(chǎn)區(qū), 也是玉米螟發(fā)生高風(fēng)險(xiǎn)區(qū), 其中德州市產(chǎn)量最高；次高地區(qū)依次為菏澤市、聊城市及濰坊市。

圖3 玉米產(chǎn)量影響區(qū)劃

2.2.3 玉米播種面積影響因子 播種面積也是影響玉米螟發(fā)生的主要影響因子, 播種面積越大, 病蟲(chóng)害越容易傳播和流行, 則受害風(fēng)險(xiǎn)性越大, 因此本研究選擇播種面積作為孕災(zāi)環(huán)境因子的影響指標(biāo)。播種面積計(jì)算公式為山東省各地市總播種面積占全省總播種面積的比重, 所用播種面積為近10 年的平均值。玉米播種面積影響因子（S）的計(jì)算公式如下。

式中,Sy為各地市的玉米播種總面積；∑Sy為全省的玉米播種總面積。

根據(jù)式（3）計(jì)算出S, 以S代表各地市的孕災(zāi)環(huán)境因子繪制玉米播種面積影響區(qū)劃（圖4）。由圖4可知, 山東省玉米主要種植區(qū)分布在該省西部地區(qū)及濰坊市中部地區(qū), 也是山東省玉米高產(chǎn)區(qū), 同時(shí)也是玉米螟可能發(fā)生的高風(fēng)險(xiǎn)地區(qū)。其中高風(fēng)險(xiǎn)區(qū)位于德州市西北部及菏澤市中西部地區(qū)；聊城市、濰坊市為較高風(fēng)險(xiǎn)區(qū)。

圖4 玉米播種面積影響區(qū)劃

3 玉米螟綜合風(fēng)險(xiǎn)指數(shù)

為綜合評(píng)價(jià)山東省各地市玉米螟發(fā)生風(fēng)險(xiǎn)狀況, 為山東省生物災(zāi)害的防治提供科學(xué)依據(jù), 對(duì)山東省玉米螟災(zāi)害發(fā)生頻數(shù)進(jìn)行綜合風(fēng)險(xiǎn)區(qū)劃[2, 4, 5, 7-9］。

3.1 數(shù)據(jù)歸一化

由于影響玉米螟災(zāi)害風(fēng)險(xiǎn)的指標(biāo)之間存在量級(jí)和量綱上的差異, 為了使不同的數(shù)據(jù)層能夠計(jì)算和比較, 需要將數(shù)據(jù)標(biāo)準(zhǔn)化處理。本研究采用極差化標(biāo)準(zhǔn)方法將單元數(shù)據(jù)轉(zhuǎn)換至[0, 1］范圍內(nèi)的值。設(shè)研究區(qū)域內(nèi)有m個(gè)評(píng)價(jià)指標(biāo), 有n個(gè)評(píng)價(jià)單元, 則有原始評(píng)價(jià)矩陣如下。

用極差標(biāo)準(zhǔn)化方法進(jìn)行如下?lián)Q算。

式中,Xij是第i個(gè)指標(biāo)第j個(gè)單元對(duì)應(yīng)的原始數(shù)據(jù)；max｛Xij｝和min｛Xij｝分別是第i個(gè)指標(biāo)下所有單元數(shù)據(jù)中的最大值和最小值；X′ij是標(biāo)準(zhǔn)化的單元數(shù)據(jù)。

3.2 綜合風(fēng)險(xiǎn)指數(shù)的計(jì)算

將極差化后的各風(fēng)險(xiǎn)指標(biāo)相加, 得到包含原來(lái)幾個(gè)風(fēng)險(xiǎn)要素綜合影響的新指標(biāo), 稱為綜合風(fēng)險(xiǎn)指數(shù), 進(jìn)而繪制綜合風(fēng)險(xiǎn)區(qū)劃, 進(jìn)行區(qū)劃分析。綜合風(fēng)險(xiǎn)指數(shù)（G）的計(jì)算公式如下。

式中,R′i、K′i、S′i分別為致災(zāi)因子、承載體因子、孕災(zāi)環(huán)境因子極差化后的指標(biāo)。

3.3 結(jié)果與分析

利用反距離權(quán)重法（IDW）、掩膜提取及自然斷點(diǎn)分級(jí)法將綜合風(fēng)險(xiǎn)指數(shù)分5 個(gè)等級(jí)進(jìn)行區(qū)劃評(píng)估, 最后得出山東省玉米螟災(zāi)害綜合風(fēng)險(xiǎn)區(qū)劃（圖5）。由圖5 可知, 綜合風(fēng)險(xiǎn)指數(shù)最大出現(xiàn)在菏澤市, 為0.897, 該地產(chǎn)量高, 種植面積大, 玉米螟發(fā)生頻率高, 是玉米螟為害造成損失的高風(fēng)險(xiǎn)區(qū)；最小出現(xiàn)在日照市, 是玉米螟發(fā)生最低風(fēng)險(xiǎn)區(qū)；其中德州市、濰坊市、濱州市都在0.53 以上, 為較高風(fēng)險(xiǎn)區(qū), 主要集中在山東省中西部主要產(chǎn)量區(qū)；日照市、威海市北部及東營(yíng)市、淄博市、棗莊市等地都在0.20 以下, 為玉米螟發(fā)生低風(fēng)險(xiǎn)區(qū), 不易受到玉米螟為害；濟(jì)南市、濟(jì)寧市、臨沂市、泰安市、青島市都處于0.30～0.53, 為玉米螟發(fā)生中風(fēng)險(xiǎn)區(qū), 條件適宜就容易發(fā)生玉米螟為害, 或有局部暴發(fā)玉米螟的風(fēng)險(xiǎn)；全省其他地區(qū)都處于0.20～0.30, 是玉米螟發(fā)生較低風(fēng)險(xiǎn)區(qū), 比較不容易受到玉米螟為害（表5）。

圖5 山東省玉米螟綜合風(fēng)險(xiǎn)區(qū)劃

4 小結(jié)

1）對(duì)研究區(qū)域內(nèi)玉米螟災(zāi)害的發(fā)生資料進(jìn)行了搜集與提取, 并將中華人民共和國(guó)農(nóng)業(yè)行業(yè)標(biāo)準(zhǔn)——玉米螟測(cè)報(bào)技術(shù)規(guī)范（NY/T 1611—2008）與山東省植物保護(hù)總站的玉米螟發(fā)生程度分級(jí)地方標(biāo)準(zhǔn)相結(jié)合, 提出了新的玉米螟發(fā)生程度分級(jí)標(biāo)準(zhǔn), 并在此標(biāo)準(zhǔn)的基礎(chǔ)上進(jìn)行數(shù)據(jù)預(yù)處理。

2）綜合考慮孕災(zāi)環(huán)境因子、致災(zāi)因子及承載體因子三方面, 建立玉米螟災(zāi)害風(fēng)險(xiǎn)評(píng)估指標(biāo)體系, 對(duì)山東省玉米螟災(zāi)害的發(fā)生風(fēng)險(xiǎn)展開(kāi)分析與評(píng)價(jià), 具體選取了玉米螟病蟲(chóng)害影響指數(shù)、玉米產(chǎn)量、玉米播種面積3 項(xiàng)影響指標(biāo), 并根據(jù)影響指標(biāo)將收集到的數(shù)據(jù)通過(guò)參數(shù)估計(jì)法風(fēng)險(xiǎn)評(píng)估模型和反距離權(quán)重法、掩膜提取法、自然斷點(diǎn)分級(jí)法進(jìn)行處理, 繪制出研究區(qū)域內(nèi)各影響因子的區(qū)劃。

3）為綜合評(píng)價(jià)山東省各地玉米螟災(zāi)害發(fā)生的風(fēng)險(xiǎn)狀況, 對(duì)山東省玉米螟進(jìn)行綜合風(fēng)險(xiǎn)區(qū)劃。為消除不同量綱的影響, 將每一指標(biāo)極差化處理, 然后將極差化后的風(fēng)險(xiǎn)指標(biāo)相加, 得到綜合風(fēng)險(xiǎn)指數(shù), 根據(jù)綜合風(fēng)險(xiǎn)指數(shù)繪制出山東省玉米螟災(zāi)害綜合風(fēng)險(xiǎn)區(qū)劃。

4）綜合風(fēng)險(xiǎn)區(qū)劃結(jié)果表明, 山東省玉米螟災(zāi)害發(fā)生的重點(diǎn)風(fēng)險(xiǎn)區(qū)集中在該省中西部地區(qū)。其中, 菏澤市屬于玉米螟發(fā)生的高風(fēng)險(xiǎn)區(qū)；德州市、濰坊市、濱州市及聊城市屬于較高風(fēng)險(xiǎn)區(qū)；濟(jì)南市、濟(jì)寧市、臨沂市、泰安市及青島市為中風(fēng)險(xiǎn)區(qū)；煙臺(tái)市為玉米螟發(fā)生較低風(fēng)險(xiǎn)區(qū)；日照市、威海市、東營(yíng)市、淄博市及棗莊市屬于低風(fēng)險(xiǎn)區(qū)。區(qū)劃結(jié)果對(duì)山東省防災(zāi)減災(zāi)事業(yè)及玉米螟災(zāi)害的防治具有重要的參考價(jià)值。