草地螟一代幼蟲最佳防治期氣象預報模型

唐紅艷,牛寶亮

(內蒙古興安盟氣象局，內蒙古 烏蘭浩特 137400)

草地螟(Loxostegesticticalis)是我國北方農牧業生產的重要害蟲，具有間歇暴發、集中遷移危害的特點[1-2]。草地螟幼蟲暴發時一夜之間便可吃光一片地，在短期內造成毀滅性災害[3]。因此，準確預測草地螟一代幼蟲危害期，提前采取有效防治措施，對最大限度減輕災害損失具有重要意義。

草地螟防治是保證農作物優質高產的重要環節，適期進行防治不僅防治效果好，還能達到少藥、高效、少污染的效果，若適期防治時間延誤一天都會對農業生產造成巨大損失。草地螟初孵幼蟲(1～2齡)主要取食雜草，食量較小，3齡以后向附近農田轉移，開始危害農作物，4～5齡進入暴食階段(占總食量的95%)，是草地螟對農作物的危害盛期[3-4]。因此，草地螟防治的最佳時期應在3齡幼蟲期之前。為了準確掌握草地螟一代幼蟲最佳防治期，前人通過對草地螟成蟲蛾峰日期的監測和分析雌蛾卵巢發育級別，并參考多年平均氣象條件來估算草地螟3齡幼蟲期，從而確定草地螟最佳防治期[5-11]。上述研究大都以成蟲發育情況為依據，根據多年經驗估算草地螟一代幼蟲最佳防治期，考慮當年環境氣象條件對草地螟產卵孵化以及低齡幼蟲發育進度的影響不足，屬于定性預報方法。溫度是昆蟲生長發育的主導因素，已有研究建立了溫度與昆蟲生長發育關系的多種模型[12]，有效積溫在昆蟲發育期和發生代數預測研究中得到廣泛應用[10,13]。事實上，昆蟲的生存、發生、發展不僅僅受溫度影響，還與其他環境氣象條件關系密切[14]。草地螟越冬代成蟲出現蛾峰后，草地螟成蟲產卵持續期、卵孵化速度以及各齡幼蟲的發育速度都受當年氣象條件影響而有所不同，因此每年的最佳防治期也會不同。本研究在分析氣象條件對草地螟發生規律影響的基礎上[15-16]，以監測到越冬代成蟲蛾峰作為初始日，建立蛾峰初始日至3齡幼蟲期持續天數定量氣象預報模型，通過預報越冬代成蟲蛾峰至3齡幼蟲期持續天數來確定草地螟最佳防治期，以期為農業部門及農民及時、科學開展草地螟防治工作提供參考依據。

1 研究區概況

興安盟(44°14′～47°39′ N，119°28′～123°38′ E)位于內蒙古自治區東北部，地處大興安嶺向松嫩平原過渡地帶，國土面積598.15萬hm2，農作物播種面積77.863萬hm2(2011年)，草原面積303.4萬hm2，屬于典型的農牧交錯區。草地資源分布廣泛，分為草甸草原、干草原、低地草甸和沼澤等。牧草種類主要有貝加爾針茅(Stipabaicalensis)、野古草(Arundinellaanomala)、羊草(Leymuschinensis)、柳蒿(Artemisiaintegrifolia)、冷蒿(A.frigida)、灰菜(Chenopodiumalbum)、豬毛菜(Salsolacollina)等。該區屬于中溫帶大陸性季風氣候，年平均氣溫-2.2～6.8 ℃，≥10 ℃活動積溫1 456.9～3 171.2 ℃·d，年降水量361.9～460.0 mm，蒸發量1 135.9～2 046.5 mm，無霜期65～160 d，光照充足。由于草地螟幼蟲屬于雜食性害蟲，半農半牧區的土地利用特點以及分布廣泛的豆科、藜科、百合科等草地螟喜食植物，為草地螟的生存、繁殖提供了良好的環境及植被條件。

2 資料與方法

草地螟發育歷期觀測資料來源于內蒙古興安盟植保站，包括1979-2010年草地螟越冬代成蟲始見期、蛾峰日期、產卵期、卵孵化期天數和各齡幼蟲期等，氣象資料來源于內蒙古烏蘭浩特市氣象局1979-2010年的逐日氣象資料，包括逐日平均氣溫、最高氣溫、最低氣溫、降水量、相對濕度以及日照時數等。農作物播種面積數據來源于興安盟統計局，草原面積數據來源于興安盟農牧業局。由于草地螟的發生具有間歇性，最長連續資料只有1996-2010年，本研究中均采用1996-2010年資料，氣候資料平均值均采用1981-2010年。對旬氣象因子進行膨化處理，生成不同組合的連續時段氣象因子集，與草地螟越冬代成蟲產卵孵化期、1～3齡幼蟲發育期、蛾峰初始日至3齡幼蟲期持續天數進行相關分析和典型年份對比分析，選取相關系數大、又具有一定物理意義的氣象因子，采用最優子集多元回歸方法建立草地螟越冬代成蟲蛾峰初始日至3齡幼蟲期持續天數預報模型。以草地螟越冬代成蟲出現蛾峰作為初始日，通過預測蛾峰初始日至3齡幼蟲期持續天數來確定草地螟最佳防治期。

3 結果與分析

3.1影響越冬代成蟲產卵孵化期持續天數的關鍵氣象因子 在草地螟越冬代成蟲蛾峰日期一定的前提下，草地螟成蟲產卵孵化期持續天數和幼蟲發育期持續天數決定了草地螟最佳防治期。6月中、下旬是興安盟地區草地螟成蟲產卵和幼蟲孵化盛期，產卵孵化期的持續天數受氣象因子影響年際間差異顯著，導致一代幼蟲發生期和最佳防治期也不同。相關分析結果(表1)表明，影響草地螟產卵孵化期持續天數最關鍵的時期是6月，關鍵氣象因子是溫度、降水量和相對濕度。產卵孵化期持續天數與產卵孵化期間的溫度呈負相關，與降水量和相對濕度呈正相關，即高溫縮短產卵孵化期持續天數，低溫延長產卵孵化期持續天數，而低濕縮短產卵孵化期持續天數，高濕延長產卵孵化期持續天數。其中,與6月中、下旬平均溫度的負相關系數通過0.05信度檢驗，與6月中旬的降水量和6月上、中旬的相對濕度正相關系數通過0.01信度檢驗，影響最顯著的因子是6月上、中旬相對濕度。

1996-2008年典型年份統計結果也表明，6月中旬平均最高溫度超過30 ℃的有5年(1997、2000、2004、2007和2008年)，這5年的溫度距平在4～6 ℃，對應的6月中旬平均溫度都超過23 ℃，6月中旬降水量不足10 mm，相對濕度都不足50%，對應上述年份的產卵孵化期持續天數都在4～6 d；而6月中旬最高溫度低于26 ℃的也有5年(1996、1998、2001、2002和2003年)，對應的6月中旬平均溫度都低于21 ℃，6月中旬的降水量都超過15 mm，相對濕度超過60%，對應上述年份的產卵孵化期持續天數都在7～9 d(表2)。說明草地螟產卵孵化期間遇高溫低濕天氣明顯縮短產卵孵化期持續天數，最佳防治期有可能提前；而低溫高濕天氣則延長產卵孵化期持續天數，最佳防治期有可能推遲。分析結果可作為預報輔助因素或模型訂正因素。

表1 草地螟產卵孵化期持續天數與氣象因子相關分析結果Table 1 Correlation analyses between spawning and hatching time and meteorological factors

3.2影響草地螟1～3齡幼蟲發育速度的關鍵氣象因子 相關分析結果(表3)表明，1～3齡幼蟲發育速度與6月中、下旬平均 (最高) 溫度呈負相關，即高溫加快幼蟲發育，縮短1～3齡幼蟲期，一代幼蟲危害期和最佳防治期提前；與6月上、中旬降水量呈正相關，即降水量多延緩幼蟲發育，1～3齡幼蟲期延長，一代幼蟲危害期和最佳防治期推遲；與6月上、中旬日照時數呈負相關。影響1～3齡幼蟲發育的關鍵氣象因子是6月上、中旬降水量，其次為6月中、下旬平均最高溫度。分析結果同樣可作為預報輔助因素或模型訂正因素。

3.3影響草地螟蛾峰至3齡幼蟲期持續天數的關鍵氣象因子 1996-2008年草地螟歷史資料表明，蛾峰至3齡幼蟲期持續天數最短為13 d，最長為28 d，持續時間相對較短。對蛾峰至3齡幼蟲期持續天數與膨化后的氣象因子進行相關分析，結果(表4)表明，持續天數與4月上、中旬日照時數，5月中旬日照時數，6月上、中旬降水量和相對濕度都呈正相關，即降水量多、相對濕度大，持續天數延長，最佳防治期推遲；與6月上、中旬平均 (最高) 溫度呈負相關，即高溫加快幼蟲發育，縮短持續天數，最佳防治期提前。此期影響最顯著的氣象因子是6月上、中旬相對濕度和降水量，相關系數達到0.754 8和0.693 8，均通過0.01信度檢驗，其次是4月中旬至6月中旬平均最高溫度和6月上、中旬平均溫度，相關系數均通過0.05信度檢驗。日照時數的影響相對弱些。

3.4越冬代成蟲蛾峰至3齡幼蟲期持續天數預報模型的建立 草地螟最佳防治期預報方法可采用分別建立產卵孵化期和1～3齡幼蟲期預報模型，也可采用建立蛾峰日至3齡幼蟲期持續天數預報模型。為便于開展草地螟防治期氣象預報服務并減小分別建立模型的預報誤差，本研究采用蛾峰日至3齡幼蟲期持續天數預報模型。以越冬代成蟲出現蛾峰日為初始日期，考慮因子物理意義和預報時效性，選取相關系數相對大的氣象因子，采用最優子集多元回歸方法建立蛾峰初日至3齡幼蟲期持續天數預報模型：

表2 產卵孵化期氣象指標Table 2 Meteorological indexes of during days of spawning and hatching

表3 1～3齡幼蟲期與氣象因子相關分析結果Table 3 Correlation analyseis between the duration days from 1-age to 3-age larvae and meteorological factors

表4 蛾峰至3齡幼蟲期持續天數與氣象因子相關分析結果Table 4 Correlation analyses between the duration days from the first day of peak catch to 3-age larvaeand meteorological factors

Y=37.5+0.137 5X1-0.255 5X2+0.171 5X3+0.559 0X4(R2=0.861 6)

模型通過0.01信度檢驗。式中，Y為蛾峰初日至3齡幼蟲期持續天數，X1為5月中旬日照時數，日照時數多則延長草地螟成蟲的產卵孵化期，日照時數少則縮短草地螟成蟲產卵孵化期；X2為4月中旬至6月中旬平均最高溫度和，溫度高，草地螟成蟲產卵孵化期縮短，溫度低則產卵孵化期延長；X3為6月上旬降水量，降水量少，幼蟲發育速度快，降水量多則幼蟲發育速度慢；X4為6月上、中旬相對濕度和，相對濕度小，幼蟲發育速度快，相對濕度大則幼蟲發育速度慢。

根據模型中氣象因子所處時段，6月下旬初即可預報草地螟蛾峰初日至3齡幼蟲期持續天數。蛾峰初日日期加上蛾峰至3齡幼蟲期持續天數就是草地螟一代幼蟲最佳防治期的開始日期。

3.5模型檢驗 利用1996-2008年的歷史資料進行模型的模擬檢驗(圖1)，有7年(54%)絕對誤差小于1 d，有4年(31%)絕對誤差在1～2 d，只有2年(15%)絕對誤差在2～4 d。平均絕對誤差只有1.2 d，模擬效果較好。

應用2009-2010年實況氣象資料進行試報檢驗，2009年蛾峰初日至3齡幼蟲期持續天數預報結果為25.6 d，實際為28 d，絕對誤差2.4 d；2010年預報結果為17.7 d，實際為15 d，絕對誤差2.7 d。總體看，預報結果基本準確，可以應用到氣象預報服務中。

圖1 蛾峰初日至3齡幼蟲期持續天數預報模型模擬結果Fig.1 Result of forecast model simulation for the duration days from the first day of peak catch to 3-age larvae

4 討論與結論

草地螟適期防治不僅防治效果好，還能夠最大限度減輕災害損失。高溫、少雨、低濕天氣加快越冬代成蟲產卵孵化及幼蟲發育，從而使蛾峰至3齡幼蟲期持續天數縮短，一代幼蟲危害期提前，最佳防治期也相應提前；而低溫、多雨、高濕天氣延緩越冬代成蟲產卵孵化及幼蟲發育，從而使蛾峰至3齡幼蟲期持續天數延長，一代幼蟲危害期推遲，最佳防治期也相應推遲。

前人對草地螟最佳防治期的預測都以定性預報為主[5-11]，估算誤差相對大，本研究以越冬代成蟲出現蛾峰為初始日期，考慮因子物理意義和預報時效性，選取相關系數相對大的氣象因子，采用最優子集多元回歸方法建立了越冬代成蟲蛾峰日至3齡幼蟲期持續天數定量預報模型。預報模型通過0.01信度檢驗，回代檢驗平均誤差僅有1.2 d，2009-2010年試報檢驗平均誤差2.55 d，模型預報誤差較小。研究結果從當年氣象條件出發，實現了定量預報草地螟一代幼蟲最佳防治期的目標，對于今后開展草地螟最佳防治期預報服務具有重要意義。

[1] 張李香,范錦勝,王貴強.中國國內草地螟研究進展[J].中國農學通報,2010,26(1):215-218.

[2] 彭秀全,劉艷,白龍,等.錫林郭勒退化草地植物功能群與草地螟種群的相關性初探[J].草業科學,2010,27(7):18-23.

[3] 崔萬里.草地螟生物學特性的觀察[J].昆蟲知識,1992(5):289-292.

[4] 徐林波,林愛萍,王慧.草地螟的生物學特性及室內毒力測定研究[J].草業科學,2007,24(9):83-85.

[5] 羅禮智,黃紹哲,江幸福,等.我國2008年草地螟大發生特征及成因分析[J].植物保護,2009,35(1):27-33.

[6] 孔海龍，羅禮智，江幸福，等.幼蟲密度對草地螟生長發育及繁殖的影響[J].昆蟲學報，2011，54(12):1384-1390.

[7] 顧成玉,梁艷春,張廣芝.草地螟種群數量變動及預測預報技術的研究[J].病蟲測報,1987(1):26-31.

[8] 羅禮智,李光博.溫度對草地螟成蟲產卵和壽命的影響[J].昆蟲學報,1993,36(4):459-464.

[9] 屈西峰，夏冰.2003年中國草地螟的發生特點和成功治理經驗[J].中國植保導刊,2004(1):22-24.

[10] 羅禮智,李光博.草地螟的有效積溫及其世代區的劃分[J].昆蟲學報,1993,36(3):332-339.

[11] 李柏年.草地螟中短期預報的初步研究[J].植物保護,2005,31(2):69-71.

[12] 時培建，池本孝哉，戈峰.溫度與昆蟲生長發育關系模型的發展與應用[J].應用昆蟲學報，2011，48(5):1149-1160.

[13] 李偉東，肖瓊，艾洪木，等.瓜實蠅的有效積溫及其在中國的年發生代數預測[J].華東昆蟲學報，2008(3)：188-193.

[14] 周華坤,王曉輝,溫軍,等.果洛州瑪沁縣草原毛蟲蟲害發生與氣候因子的相互關系[J].草業科學,2012,29(1):128-134.

[15] 唐紅艷,牛寶亮.影響草地螟一代幼蟲發生的關鍵氣象因子及指標研究[J].中國農學通報,2010,26(18):270-274.

[16] 唐紅艷.草地螟一代幼蟲發生面積的氣象預報模型[J].中國農業氣象,2011,32(1):144-147.