2. 河南省濮阳市植保站
2. Puyang Plant Protection Station
棉铃虫的发生发展与危害程度,主要取决于虫源、食物源、防治技术与措施、气象条件等多项因素。就宏观情况而言,在棉铃虫常发区,虫源、食物源普遍存在;防治技术与措施年际间也不会有大的改变,这两个条件都不会引起年际间棉铃虫发生与危害程度出现大的差别;气象条件年际间差别对棉铃虫的发生、发展与危害程度影响较大。正因为如此,人们常用气象条件制作棉铃虫预报,气象工作者还建立了一些棉铃虫气象指标和预报方法[1-4],但以往的指标与方法,对棉铃虫危害期的气象条件考虑较多,对虫源形成期的气象条件考虑较少。其实,棉铃虫的发生及危害程度具有年际、代际间传承的特点,前期的虫源条件对夏季棉铃虫发生危害程度影响较大。作者根据棉铃虫周年活动规律,将棉铃虫周年繁育期划分为秋季预蛹期、越冬期、春季羽化期、夏季危害期4个阶段,分阶段建立棉铃虫分析评价指标与预报模式,取得较好效果。
1 资料来源与处理气象资料来自当地气象部门,棉铃虫资料来自豫北各市植保部门。对比分析各地虫情资料发现,尽管各地虫情资料时间序列长短不一,并多有缺漏现象,但仍可看出,豫北各市历年棉铃虫发生趋势比较一致,故在建立豫北棉铃虫发生程度时间序列时,采用资料年代较长(1976—2005年),又较完整的新乡市棉铃虫资料代之,其余地市的资料仅作参考使用。棉铃虫发生程度采用植保部门《棉铃虫调查测报规范》5级表示法,1级为轻度发生,2级为中度偏轻发生,3级为中度发生,4级为中度偏重发生,5级为重度发生,各代棉铃虫的发生程度不尽相同,取其中发生程度最重一代的发生程度为年发生程度。
2 棉铃虫繁育与气象条件棉铃虫在豫北地区一年发生4代,9月下旬当气温降到20℃以下时,老熟幼虫开始在5~15cm深处的土层筑土室进入冬前滞育期。秋季降温早,棉铃虫滞育也早,反之,则滞育晚。滞育阶段棉铃虫滞育蛹如遇到一段时期的高温天气过程(例如9或10月份气温偏高),造成部分蛹羽化,而后遇到秋季低温、霜冻,将导致全部羽化个体死亡。忽高忽低的剧烈变温是导致该阶段棉铃虫死亡的主要气候原因。
11月初,当气温降到10℃时,棉铃虫便以蛹态进入越冬期,棉铃虫越冬的关键因素是最冷月气温的高低,冷冬不利于棉铃虫越冬,暖冬有利于棉铃虫越冬。
翌年4月下旬气温升至15℃以上时,越冬蛹开始羽化为成虫,5月上旬为羽化盛期,第一代幼虫主要危害小麦。春季气温回升早,棉铃虫羽化早,反之,则羽化晚。影响春季羽化阶段棉铃虫存活的关键因素是春季天气的不稳定程度,尤其是倒春寒、低温、霜冻、降雨天气将大大降低棉铃虫成虫羽化率,春季霜冻终止日期出现越晚、强度越强,降水越多,棉铃虫死亡率越高。
第一代成虫盛发期为6月上、中旬,产卵盛期在6月中、下旬。6月底为第二代(棉田第一代)幼虫危害盛期,成虫盛发期在7月上、中旬,产卵盛期在7月中旬。7月下旬为第三代幼虫危害盛期,成虫盛发期在8月上、中旬。8月下旬到9月上旬为第四代幼虫危害盛期。在豫北第二、三代幼虫危害棉花最严重,第四代主要发生在贪青的棉田。此期棉铃虫发育的最适温度为25~30℃,相对湿度为70%~95%。棉铃虫适宜于半干旱的气候条件,若降水过多,造成土壤过湿或板结不利于棉铃虫入土化蛹, 同时降水强度大时还可冲刷掉部分虫卵与幼虫,据观察,伴有大风的阵雨常可冲掉60%~70%的卵,对棉铃虫的种群增长及危害起到抑制作用[5]。
3 棉铃虫气象指数的建立 3.1 阶段划分根据棉铃虫的繁育规律,并结合棉铃虫气象条件分析,将棉铃虫周年活动期划分为4个阶段。第一阶段为9月中旬至10月下旬,称作秋季预蛹期,此期的气象条件影响着进入越冬期蛹量的多少;11月上旬至翌年3月下旬,为第二阶段,称作越冬期,此期气象条件直接影响越冬蛹的存活量;4月上旬至5月中旬为第三阶段,称为春季羽化期,此期的气象条件对越冬蛹的成虫羽化率影响较大;第四个阶段为5月下旬至9月上旬,称为夏季危害期,此期的气象条件适与否,对2、3、4代棉铃虫的繁育量影响较大。前3个阶段的气候条件主要影响当年虫源基数,是制作棉铃虫预报的重要基础。夏季的气象条件则直接影响虫害发生与危害程度,是制作棉铃虫预报的主要对象。
3.2 气象因子的选取根据棉铃虫各阶段活动规律及其对气候环境的不同要求,通过棉铃虫发生程度与气象因子的相关分析,分阶段选取气候因子,共选取有统计学意义与昆虫学意义的因子11个,其中有4个因子通过α≥0.01的相关显著性检验,有2个因子通过α≥0.05的相关显著性检验,其余5个因子虽未通过检验,但相关系数较大,仍可作进一步处理应用(见表 1)。
在筛选因子时,虽兼顾了统计学与昆虫学意义两个方面,但若用这些因子直接建立预报模式或建立棉铃虫气象指标,其回代与外延精度仍不能得到满意效果。为提高预报精度,需对所选因子进行重新组合。影响棉铃虫发生程度的气象因素较多,各种因子对棉铃虫的影响程度各不相同,其量纲也不同,往往不能进行直接组合,为解决这一问题,需要对各气象因子作无量纲化处理。对与棉铃虫发生程度呈正相关的因子,用式(1)进行计算[4]:
(1) |
对与棉铃虫发生程度成负相关的因子用式(2)进行计算:
(2) |
式中:Yij为规格化因子,Xij为气象因子,Xjmax和Xjmin分别为Xij的最大值和最小值,i为年代序数,j为因子序数,下同。处理后的归一化因子理论最大值为1,最小值为0,因子消除了量纲差异,按前面划分的4个阶段,分别用式(3)进行加权组合:
(3) |
式中Aj为权重系数,由各因子与棉铃虫发生程度之间的相关系数求得,k为气象指数序数,m为各组合式中组合因子的个数。
本文将Rik定名为棉铃虫气象指数(见表 2)。棉铃虫气象指数系由与其发生程度密切相关的气象因子按照相关系数的大小组合而成,其值大,表示气象条件有利棉铃虫发生,反之则不利棉铃虫的发生。通过不同发病年份的气象指数分析和回归分析,可进一步制定出棉铃虫气象指数指标并建立棉铃虫发生程度预报模式,依此对各阶段棉铃虫气象条件作定量分析评价与预报。
归一化棉铃虫气象指数,系由与棉铃虫发生程度密切相关的气象因子,以相关系数为权重组合而成,犹如提取了各因子的有用成分(事实是各因子之间在一定程度上存在着互补作用),不但使各因子原有的生物学意义得到保留,还使其与棉铃虫发生程度之间的相关性明显增强,Ri1、Ri2、Ri3、Ri4与棉铃虫发生程度之间的相关系数分别达到0.8332、0.4660、0.6061和0.4299,经检验,前三个阶段的相关系数通过了α≥0.01的相关显著性检验,后一阶段则通过α≥0.05的相关显著性检验。将Ri1、Ri2、Ri3、Ri4按时段顺序再次作加权组合,并使后一阶段的组合中包含有前期的有效气象信息。组合后,得到三组气象指数Ri5、Ri6、Ri7(注:Ri5为Ri1、Ri2加权集成值,Ri6为2Ri1、Ri2、Ri3加权集成值,Ri7为Ri1、Ri2、Ri3、Ri4加权集成值),这些组合因子与棉铃虫发生程度之间的相关性更为显著,相关系数分别达到0.7943、0.8378和0.8574,均通过α≥0.01的显著性检验。将Ri1、Ri5、Ri6、Ri7分别命名为棉铃虫秋季预蛹期、越冬期、春季羽化期、夏季危害期归一化气象指数。
1976—2005年,秋季预蛹期气象指数Ri1的平均值为0.58,最大值为0.85,最小值为0.35。棉铃虫发生程度所对应的气象指数范围分别为:1级0.38~0.46;2级0.35~0.46;3级0.49~0.74;4级0.62~0.67;5级0.69~0.85。据此制定出棉铃虫不同发生年份,秋季预蛹期气象指数指标为:1、2级Ri1≤0.50;3、4级0.51≤Ri1≤0.70;5级Ri1≥0.68。用此指标评价历年秋季预蛹期的气象条件,30年中符合年份有26年,不相符的年份有4年,符合率为84%。
越冬期气象指数Ri5的平均值为0.60,最大值为0.79,最小值为0.32。棉铃虫各级发生程度所对应的气象指数范围分别为:1级0.39~0.51;2级0.32~0.52;3级0.53~0.77;4级0.65~0.72;5级0.70~0.78。据此制定出棉铃虫不同发生年份,越冬期气象指数指标为:1、2级Ri5≤0.52;3、4级0.53≤Ri5≤0.69;5级Ri5≥0.70。用此指标评价历年越冬期的气象条件,30年中符合年份有25年,不相符的年份有5年,符合率为83%。
春季羽化期气象指数Ri6的平均值为0.54,最大值为0.78,最小值为0.31。棉铃虫不同发生程度所对应的气象指数范围分别为:1级0.31~0.45;2级0.34~0.46;3级0.44~0.69;4级0.56~0.72;5级0.63~0.78。据此制定出棉铃虫不同发生年份,春季羽化期气象指数指标为:1、2级Ri6≤0.45;3、4级0.46≤Ri6≤0.64;5级Ri6≥0.65。用此指标评价历年越冬期的气象条件,30年中符合年份有25年,不相符的年份有5年,符合率为83%。
夏季危害期气象指数Ri7的平均值为0.604,最大值为0.80,最小值为0.40。棉铃虫不同发生程度所对应的气象指数范围分别为:1级0.40~0.53;2级0.41~0.54;3级0.51~0.69;4级0.64~0.71;5级0.69~0.80。据此制定出棉铃虫不同发生年份,夏季危害期气象指数指标为:1、2级Ri7≤0.54;3、4级0.55≤Ri7≤0.69;5级Ri7≥0.70。用此指标评价历年越冬期的气象条件,30年中符合年份有25年,不相符的年份有5年,符合率为83%。
将以上统计分析结果列于表 3,表中各时段的气象指数指标值,是根据历年棉铃虫发生程度,对照各时段气象指数的实况值而制定的,各时段的指标相互独立,每一时段内虫害等级间的指标相互衔接,无交叉现象。秋、冬、春、夏4个阶段的气象指数统计分析结果表明,气象指数与棉铃虫发生程度密切相关,各阶段棉铃虫气象指数指标的回代准确率均在83%以上,此指标可以用于棉铃虫气象条件的分析评价与预报。
将各阶段入选的气象因子作回归分析,可组建多时效棉铃虫预报模式,模式通式为:
式中z为棉铃虫发生程度的预报值,a、b为回归系数,x为归一化气象因子。
用归一化气象因子和气象指数,逐时段建立棉铃虫发生程度预报方程,共建立4组方程。考虑到前后期气象条件对棉铃虫影响的累积效应,将每一方程的预报、回代结果引入下一阶段的预报方程中作预报因子使用,这样既充分利用了与棉铃虫密切相关的前期气象信息,也使预报、回代精度得到了逐步提高, 其形式见表 4。每时段有两个方程,①为归一化气象因子组成的预报方程,②为气象指数(归一化气象因子的加权集成值)所组成的预报方程,预报结论取两方程的平均值。
经检验,各个方程的复相关系数在0.7780~0.8780之间,F检验值在9.5846~63.5625之间,均通过回归检验。若以回代误差小于1为准,30年中只有一年的绝对误差大于1,回代准确率为96.7%,历年回代情况见图 1。
将2005年度和2006年度各阶段的气象因子代入与之对应的指标和模式中进行试报应用,取得较好效果(见表 5,6,7)。
将表 5中的归一化气象因子和表 6中的归一化气象指数,分别代入表 4中的相关方程,把每组方程的预报值作等权集成,即可作出不同时效的棉铃虫发生程度预报,由表 7可以看出,各时段的试报结果均正确。
5 结论(1) 影响棉铃虫繁育的气象因子较多,将诸多气象因子作无量纲化处理,不但便于因子间相互组合,还可提高模式的预报精度。
(2) 棉铃虫的发生程度具有世代传承的特点,每年的发生危害程度不仅与发生危害期的气象条件有关,还与前期的气候条件及上年的虫原残留量有关。将棉铃虫周年繁育期划分为秋季预蛹期、越冬期、春季羽化期、夏季危害期,既突出了阶段性特点,又便于把棉铃虫周年活动期的气候变化作为一个连续体,分阶段组建预报模式,预报效果较好。
(3) 棉铃虫多时效预报模式的预报时效超前,不但对于棉铃虫的防御有积极意义,而且对农药的生产、调运与购进也有指导作用。
(4) 豫北各市气候背景基本一致,每年棉铃虫的发生趋势也基本一致,在制作棉铃虫发生程度预报时,可用以点代面的方法,通过一站资料的深入分析,作出全区的预报。
(5) 对组合因子的生物学意义只作了初步推断分析,有待通过生物学试验,进一步确定其准确的生物学意义。
[1] |
周志香, 王志伟, 刘文平, 等. 棉铃虫综合因素预报模型研究[J]. 农业气象, 2000, 21(4): 38-42. |
[2] |
于玲, 杨彬云, 相云, 等. 不同气候年型棉铃虫气象指标系统的研究[J]. 农业气象, 2000, 21(3): 27-32. |
[3] |
侯敬和, 刘厚赞. 棉铃虫发生程度长期预报方法[J]. 气象, 1997, 23(2): 41-44. DOI:10.7519/j.issn.1000-0526.1997.02.011 |
[4] |
张建华, 李迎春. 新疆棉铃虫越冬气候指标的分析[J]. 气象, 2000, 23(2): 41-44. DOI:10.7519/j.issn.1000-0526.2000.02.010 |
[5] |
河南农业大学. 农业昆虫学[M]. 郑州: 河南农业大学出版社, 1985: 219-225.
|
[6] |
贾金明. 黄河中下游小麦赤霉病气象指数的建立与应用[J]. 气象, 2002, 28(3): 50-53. DOI:10.7519/j.issn.1000-0526.2002.03.011 |