王纯枝 霍治国 张 蕾 郭安红 黄 冲 陆明红
1)(国家气象中心, 北京 100081)2)(中国气象科学研究院, 北京 100081)3)(南京信息工程大学气象灾害预报预警与评估协同创新中心, 南京 210044)4)(全国农业技术推广服务中心, 北京 100125)
小麦是我国重要的粮食作物之一,北方地区是我国冬小麦的主产区,其产量对我国粮食安全占有举足轻重的地位,而小麦蚜虫是危害小麦产量和品质的主要虫害[1-3],研究我国北方主要小麦种植区域麦蚜气象适宜度预报预警技术,对做好气象保障、降低小麦产量损失、确保粮食安全具有重要意义。小麦蚜虫俗称油虫、腻虫、蜜虫,属同翅目(Hemiptera),蚜科(Aphidoidea),可对小麦进行刺吸危害,麦长管蚜和禾谷缢管蚜是影响我国小麦生产的最主要害虫[4-5],麦长管蚜则是危害的优势种[3,6]。
我国北方小麦主产区以麦长管蚜(多为害上部叶片,抽穗灌浆期集中穗部为害)和麦二叉蚜(喜在苗期为害)发生数量最多,为害最重,其直接为害以成蚜、若蚜吸食小麦叶片、茎秆和嫩穗的汁液,使被害部位形成条斑、枯萎,植株生长停滞甚至整株枯死,影响小麦光合作用及营养吸收、传导,穗期易造成麦粒不饱满或形成秕粒,使千粒重降低,造成减产;间接为害是麦蚜传播小麦病毒病[3],如为害较大的小麦黄矮病。小麦蚜虫在南方无越冬期,在北方麦区以无翅胎生雌蚜在麦株基部叶丛或土缝内越冬,或者以卵在麦苗枯叶上、杂草上、茬管中、土缝内越冬。麦蚜的越冬、发育和繁殖与气象条件密切相关[7-10],如麦长管蚜通常喜中温不耐高温,适宜温度范围为13~25℃,适宜空气相对湿度范围为40%~80%[11];麦二叉蚜则喜干怕湿,5 d平均气温16~25℃和空气相对湿度35%~67%为其适宜的繁殖条件[12-13]。
气候条件是小麦蚜虫暴发成灾的主要影响因子[14]。暖冬可使小麦蚜虫越冬存活率增加,发生期提前[15],暖春有利于虫害扩展速度加快[16],气候变暖可使小麦蚜虫虫害发生趋势加重[11]。几乎所有大范围流行性农作物重大病虫害的发生发展均与气象条件密切相关[7,10],温度、空气湿度等直接与害虫的生殖、行为、种群演变、分布和发生程度有关[9-10,17]。目前,基于传统病虫害预测预报方法包括经验预测法、实验预测法和统计预测法[18-19]。有学者开始运用现代非线性理论,将人工神经网络、相空间重构预测法[20]、小波分析[21]、马尔可夫链[22]、支持向量机和局部支持向量回归[1]、模糊认知图(FCM)算法[8]等方法运用于蚜虫等害虫预测预报中,但由于算法复杂、训练过程长、维数较高、参数较多或不能满足全局算法一致等局限性,与传统的统计预测法相比,在业务中应用有限,统计预测法因组建模型相对简单、参数获取简便、应用推广便捷等优点,仍广泛应用[3,23-25]。因此,本研究根据我国北方主要小麦种植区域1958年以来的小麦蚜虫发生面积、发生程度和气象资料,结合小麦发育期区域分布规律和麦蚜自身生物学特性,采用相关分析法、逐步回归法等,研究小麦蚜虫年发生程度与上年冬季至当年6月上旬气象因子间的关系,筛选与小麦蚜虫发生程度相关显著的关键气象因子及时段,构建北方小麦主产区不同区域的小麦蚜虫气象适宜度预报模型,用于评估危害程度,为政府和生产部门及时采取防范措施、有效防控提供科技支撑和决策依据。
研究区为我国北方冬小麦主产区8个省(市),包括河北、北京、天津、山西、河南、山东、江苏和安徽。1949—2018年各省小麦蚜虫年发生面积、发生程度等级、实际损失及2012—2018年始发期、盛发期资料来源于全国农业技术推广服务中心,其中发生程度等级序列连续性偏差,山西1949—1972年资料整体欠缺。1981—2018年冬小麦常年发育期和不同年份发育期资料来源于国家气象中心。气象资料来源于国家气象中心,由于气象资料从1958年开始有较完整记录,故采用长度为1958—2018年,要素包括上年12月至当年6月上旬的逐日平均气温、最高气温、最低气温、降水量、日照时数、空气相对湿度、大雨(暴雨)日数等。另外,引进了温雨系数指标C,即
C=P/T。
(1)
式(1)中,P为月或旬累积降水量(单位:mm),T为月或旬平均气温(单位:℃)。
考虑到各省虫害实际发生程度等级资料序列不完整,而发生面积序列相对完整,为使建立的小麦蚜虫气象等级指标具有业务实用性,参照于彩霞等[26]对稻飞虱发生面积率的划分标准,计算1958—2018 年各省小麦蚜虫发生面积平均值,以平均值上下波动50% 为等级间隔,对小麦蚜虫发生面积进行轻、偏轻、偏重和重4个等级划分,1级表示轻发生(对应气象等级为不适宜),2级表示偏轻发生(对应气象等级为较适宜),3级表示偏重发生(对应气象等级为适宜),4级表示重发生(对应气象等级为非常适宜)。小麦蚜虫测报国家标准[27]将蚜虫发生程度划分为5级(轻发生、偏轻发生、中等发生、偏重发生、大发生),为了保持一致,验证时将中等发生和偏重发生合并为偏重等级,即合并后的偏重发生等级对应发生面积分级为3级,大发生对应发生面积分级为4级,其余等级一一对应不变。表1给出了各省小麦蚜虫发生面积分级及对应发生程度等级划分。经对各省发生面积分级与实际发生程度等级进行一致性分析,发现华北、黄淮等级匹配一致性准确率分别为90.3%,95.8%(华北样本量n=93,黄淮样本量n=72),因此,采用发生面积分级代替实际发生程度等级可行。
表1 小麦蚜虫发生等级划分Table 1 Classification for occurrence area of wheat aphids
气象站点的选择依据中国气象地理区划[28],并利用中国土地利用图[29]和ArcGIS分析工具中的叠加Intersect功能剔除高山、城市开发区等非农田站点,操作在ArcMap10.2环境下完成,最终选取8个省(市)601个气象站点(图1),其中在华北选取246站,包括北京10站、河北133站、天津11站、山西92站;在黄淮选取216站,包括山东106站、河南110站;在苏皖两省选取139站,包括安徽70站、江苏69站。时间均为1958—2018年。对地面气象资料,先将气象站点上年12月至当年6月上旬的逐日资料处理成旬平均和月平均资料,再计算各省(市)气象要素和温雨系数的区域平均,以及上年冬季平均气温、最高气温、最低气温的逐年区域平均值,以便与各省虫害发生面积相对应。
图1 研究区内气象站点分布Fig.1 Distribution of meteorological stations in the study area
1.3.1 相关分析
小麦蚜虫发生程度以年发生面积为依据,因子普查采用Pearson相关分析法,分析年发生程度与上年12月至当年6月上旬不同旬、月时段组合地面气象要素和上年冬季平均气温、最高气温、最低气温间的相关关系,筛选与麦蚜发生程度显著相关的华北、黄淮关键气象因子及影响时段,相关系数的检验采用双尾t检验。
1.3.2 独立性检验
选用主成分识别法进行因子的独立性检验[30-32],剔除共线性因子,建立因子间的相关矩阵R。
(2)
剔除共线性自变量数量的确定方法:rij为因子i和j之间的相关系数,求相关系数矩阵R的特征值λi,若因子间存在共线性,则有特征根接近0。因此,如果有k个特征量值近似等于0,则证明因子之间有k重共线性,需要剔除k个自变量。主成分分析法的具体步骤详见文献[33-34]。
1.3.3 因子归一化
在筛选出关键因子后,由于各因子影响程度不同、量纲不同,不便直接应用,因此,对各关键气象因子采用极差化方法进行标准化处理(归一化法,使数据均分布在0~1之间)。
1.3.4 气象适宜度综合指数和模型构建
利用华北、黄淮分区域归一化处理后的关键气象因子与对应省份小麦蚜虫发生面积划分等级建立华北、黄淮小麦蚜虫气象适宜度预报模型,小麦蚜虫发生发展气象适宜度综合指数计算公式为
(3)
式(3)中,Z为小麦蚜虫发生发展气象适宜度综合指数,C0为常数,Yi为归一化后的第i项关键气象因子,Ki为第i项归一化关键因子的回归系数,n为因子个数。
将1958—2015年华北、黄淮地区资料用于建模和模型拟合回代检验,2016—2018年华北、黄淮及苏皖两省资料用于模型外推预报检验。
2.1.1 关键气象因子筛选
根据小麦蚜虫生物学特性[4-5]和小麦常年发育期规律[35],华北、黄淮小麦蚜虫为害特点见表2。小麦蚜虫主要为害期为冬小麦拔节至乳熟前期,盛发期为冬小麦抽穗开花期;华北麦蚜主要为害期通常为4月上旬至5月中旬,盛发期为4月下旬至5月上旬;黄淮麦蚜主要为害期为3月下旬至5月上旬,盛发期为4月中下旬(表2)。小麦蚜虫为害盛期,蚜虫种类主要是小麦穗期优势种麦长管蚜[5]。分析春季至夏初冬小麦主产区8个省(市)逐旬平均气温变化发现,4月上旬安徽、河南两省平均气温已达到适宜麦长管蚜发生发展的基点气温(13℃)[11],4月中旬至6月上旬8个省(市)逐旬平均气温基本都在麦长管蚜发生发展的适宜生理气象指标范围(13~25℃),仅4月中旬山西省平均气温偏低,说明小麦蚜虫为害盛期气温条件并非影响其发生程度的主导限制因子。因此,探讨冬春季的气象条件对麦长管蚜的影响及其相互关系尤为必要。
表2 北方小麦主产区小麦蚜虫为害特点[4-5,27,35]Table 2 Damage characteristics of wheat aphids in the main wheat growing areas in the northern China(from Reference [4-5,27,35])
首先采用单因子相关分析方法,计算华北、黄淮各省小麦蚜虫发生面积与上年12月至当年6月上旬各旬(月)气象因子包括平均气温、最高气温、最低气温、降水量、日照时数、空气相对湿度、大雨(暴雨)日数、无雨日数、温雨系数以及上年冬季平均气温、最高气温、最低气温等的相关系数,筛选出相关系数达到0.05显著性水平的因子共40个,并进行独立性检验,结合北方小麦主产区蚜虫为害特点(表2),最终筛选出影响华北地区小麦蚜虫年发生程度的关键气象因子8个,分别为上年冬季平均气温、当年3月温雨系数、3月最高气温大于等于25℃的日数、3月下旬日照时数、4月上旬平均气温、4月下旬最高气温大于等于28℃的日数、4月大雨日数和5月上旬空气相对湿度介于40% ~ 80%范围的日数;筛选出影响黄淮地区小麦蚜虫年发生程度的关键气象因子6个,分别为上年冬季平均气温、当年1月下旬降水量、3月上旬空气相对湿度大于80%的日数、3月平均气温、4月温雨系数和4月下旬无雨日数。具体见表3和表4。
由表3、表4可见,小麦蚜虫发生与冬季、春季尤其春季的气象因子密切相关,冬春不同时段平均气温与华北、黄淮小麦蚜虫发生面积均呈显著正相关,冬季气温偏高,利于小麦蚜虫越冬基数提高,春季气温偏高,利于小麦蚜虫发生发展,这与文献[3,8,15,24]的研究结果一致,热量对小麦蚜虫发生发展具有正相关协同作用。3月、4月温雨系数与华北、黄淮小麦蚜虫发生程度呈显著负相关,春季降水偏少、气温偏高,温雨系数偏小,利于苗蚜和穗蚜的发生发展。华北小麦蚜虫发生面积与3月最高气温大于等于25℃的日数和5月上旬空气相对湿度为40%~80%的日数呈显著正相关,与3月温雨系数、3月下旬日照时数、4月下旬最高气温大于等于28℃的日数和4月大雨日数呈显著负相关。黄淮小麦蚜虫发生面积与1月下旬降水量、3月上旬空气相对湿度大于80%的日数和4月温雨系数呈显著负相关,与4月下旬无雨日数呈显著正相关。这是由于小麦蚜虫在不同虫态、不同小麦发育期对空气湿度、日照时数、水分条件的要求存在差异所致。据研究[24],热量和水分条件是影响小麦蚜虫种群消长的关键气象因素,水分因子(降水和湿度)具有反向抑制作用,其中大雨有冲刷作用,大雨日数越多,越利于抑制小麦蚜虫发生发展;相反,无雨日数越多,越有利于小麦蚜虫的暴发流行。
表3 华北地区小麦蚜虫发生程度与关键气象因子相关系数Table 3 Correlation coefficients between the occurrence area of wheat aphids and key meteorological factors in North China
注:①T,Tmx,C,D,S,H,P分别表示平均气温、最高气温、温雨系数、日数、日照时数、空气相对湿度、降水量,各指标后下标数字表示╳月╳旬,如T41,C3分别表示4月上旬平均气温、3月温雨系数,Tw表示冬季平均气温,以此类推。下同。DTmx3,DTmx43,DP4,DH51分别表示3月最高气温大于等于25℃的日数、4月下旬最高气温大于等于28℃的日数、4月大雨日数、5月上旬空气相对湿度介于40%~80%的日数。
②*,**和***分别表示达到0.05,0.01和0.001显著性水平(样本量为217)。
表4 黄淮地区小麦蚜虫发生程度与关键气象因子相关系数Table 4 Correlation coefficients between the occurrence area of wheat aphids and key meteorological factors in the Huanghuai Area
注:P表示降水量。DH31,DP43分别表示3月上旬空气相对湿度大于80%的日数、4月下旬无雨日数,其余各指标后下标紧跟数字同表3类推。*,**和***分别表示达到0.05,0.01和0.001显著性水平(样本量为116)。
2.1.2 小麦蚜虫气象适宜度模型的建立
用华北区域对应省份小麦蚜虫发生面积等级序列与华北归一化的8个关键气象因子进行回归分析,建立华北小麦蚜虫气象适宜度预报模型如下:
Znc=2.304+1.423Y1-0.88Y2-
0.267Y3-1.399Y4+1.343Y5-0.474Y6-
0.091Y7+0.247Y8。
(4)
其中,Znc为在具有一定虫源条件下仅考虑气象条件影响的当年华北小麦蚜虫发生发展气象适宜度指数,Y1,Y2,Y3,Y4,Y5,Y6,Y7,Y8分别为上年冬季平均气温(Tw)、当年3月温雨系数(C3)、3月最高气温大于等于25℃的日数(DTmx3)、3月下旬日照时数(S33)、4月上旬平均气温(T41)、4月下旬最高气温大于等于28℃的日数(DTmx43)、4月大雨日数(DP4)和5月上旬空气相对湿度为40%~80%的日数(DH51)。方程复相关系数为0.682,达到0.001显著性水平,样本量n=217。
用黄淮区域对应省份小麦蚜虫发生面积等级序列与黄淮归一化的6个关键气象因子进行回归分析,构建黄淮小麦蚜虫气象适宜度预报模型如下:
Zhh=0.605+2.085Y1-0.725Y2-1.41Y3+
0.796Y4+0.443Y5+1.667Y6。
(5)
其中,Zhh为当年黄淮小麦蚜虫发生发展气象适宜度指数,Y1,Y2,Y3,Y4,Y5,Y6分别为上年冬季平均气温(Tw)、当年1月下旬降水量(P13)、3月上旬空气相对湿度大于80%的日数(DH31)、3月平均气温(T3)、4月温雨系数(C4)和4月下旬无雨日数(DP43)。方程复相关系数为0.688,达到0.001显著性水平,样本量n=116。
Z值越大表示气象条件对小麦蚜虫发生发展越有利,反之越不利。为方便实际应用,将北方地区小麦蚜虫资料中华北217个样本和黄淮116个样本的小麦蚜虫发生程度从小到大排序,分析虫害发生程度分别为轻(1级)、偏轻(2级)、偏重(3级)和重(4级)(对应小麦蚜虫测报标准[27]中等和偏重合并后)4个不同级别的样本分布,将气象适宜度指数Z值划分为4个级别(表5):1级为气象条件不适宜虫害发生发展,Z<1.5;2 级为气象条件较适宜虫害发生发展,1.5≤Z<2.5;3 级为气象条件适宜虫害发生发展,2.5≤Z<3.5;4 级为气象条件非常适宜虫害发生发展,Z≥3.5。
表5 华北、黄淮小麦蚜虫发生发展气象适宜度指数分级表Table 5 Classification of meteorological suitability index for the occurrence and development of wheat aphids in North China and the Huanghuai Area
利用1958—2015年资料进行模型回代拟合检验,华北和黄淮小麦蚜虫气象适宜度等级与实际发生等级相比,区域模型拟合级别一致平均准确率均在90%以上(表6)。由表6两区域分级拟合结果可知,4个级别拟合准确率均不小于80%,2级(气象条件较适宜)和3级(气象条件适宜)预报准确率均达到90%以上,其中华北2级、3级基本正确的准确率达100%,气象条件较适宜和适宜等级的预报效果比气象条件不适宜(1级)明显好,即随着预报适宜程度增加,预报准确率明显提升;4级(气象条件非常适宜)预报准确率较2级和3级虽有所下降,但华北4级预报准确率在85%以上,黄淮在90%以上,预报效果仍较好。
表6 1958—2015年华北和黄淮小麦蚜虫气象适宜度模型回代拟合准确率Table 6 Hindcast accuracy of forecast models of meteorological suitability for wheat aphids in North China and the Huanghuai Area from 1958 to 2015
注:预报与实际一致为正确,相差1个等级为基本正确,相差2个或2个以上等级为错误。
利用分区域虫害气象适宜度模型对2016—2018年华北和黄淮小麦蚜虫气象适宜度等级进行预报,利用黄淮区域气象适宜度模型对江苏、安徽两省2016—2018年小麦蚜虫气象适宜度等级进行外推预报检验,结果见表7和表8。由表7可以看到,与实际发生面积分级相比,模型对2016—2018年华北和黄淮小麦蚜虫气象适宜度等级预报平均准确率分别为75%和100%;黄淮模型外推至苏皖地区,两省预报3年平均准确率为100%(表8),效果理想。研究还表明:建立的小麦蚜虫气象等级预报模型对蚜虫发生偏重年份预报准确率较偏轻年份更高,说明模型对偏重年份的灾害等级的反映和响应效果更好。
根据2016—2018年农业农村部对虫情实际发生等级监测结果(表7和表8),2016—2018年河北省小麦蚜虫实际发生程度分别为4级、4级和3级,山西、山东、河南和安徽4省分别均为3级,天津市和江苏省分别均为2级、3级、3级,北京均为2级。经检验,气象适宜度等级与小麦蚜虫实际发生等级监测结果相比,模型预测华北、黄淮及安徽、江苏两省基本正确(相差不超过1个等级)的平均准确率分别为91.7%,100%,100%和100%,模型预报效果较好。由表8可以看到,利用黄淮模型外推预报苏皖两省2016—2018年小麦蚜虫气象适宜度等级,与实际发生等级相比,两省预报3年平均基本正确准确率为100%;单独看,两省3年中预报等级每年均正确或基本正确,基本正确准确率均为100%,模型外推预报效果较为理想。因此,建立的分区域模型适用于华北、黄淮和江淮地区小麦蚜虫发生发展气象适宜度等级预报。
表7 2016—2018年华北和黄淮小麦蚜虫气象适宜度预报模型预报准确率Table 7 Extrapolated accuracy of forecast models of meteorological suitability for wheat aphids in North China and the Huanghuai Area from 2016 to 2018
注:预报与实际一致为正确,相差1个等级为基本正确,相差2个或2个以上等级为错误。
表8 2016—2018年黄淮小麦蚜虫气象适宜度预报模型外推预报准确率Table 8 Extension forecasting accuracy of meteorological suitability models for wheat aphids in the Huanghuai Area from 2016 to 2018
注:预报与实际一致为正确,相差1个等级为基本正确,相差2个或2个以上等级为错误。
本文基于1958—2015年北方小麦主产区8个省(市)小麦蚜虫逐年发生面积及对应的601个气象站逐日气象资料,采用相关分析、主成分分析和逐步回归等方法,从上年12月至当年6月上旬各气象要素不同时段的旬值、月值和上年冬季平均气温、最高气温、最低气温中,筛选出影响华北和黄淮小麦蚜虫年发生程度的关键气象因子及其影响时段,构建分区域的小麦蚜虫气象适宜度指数和预报模型,可支持业务应用,得到以下主要结论:
1) 影响华北小麦蚜虫发生发展的关键气象因子有8个,分别是上年冬季平均气温、当年3月温雨系数、3月最高气温大于等于25℃的日数、3月下旬日照时数、4月上旬平均气温、4月下旬最高气温大于等于28℃的日数、4月大雨日数和5月上旬空气相对湿度为40%~80%的日数;影响黄淮麦蚜发生发展的关键气象因子有6个,分别是上年冬季平均气温、当年1月下旬降水量、3月上旬空气相对湿度大于80%的日数、3月平均气温、4月温雨系数和4月下旬无雨日数。
2) 利用构建的分区域小麦蚜虫气象适宜度模型对华北和黄淮麦蚜发生发展气象等级进行预报,2016—2018年3年平均试报准确率均在75%以上,预报效果较好。利用黄淮小麦蚜虫气象适宜度模型预报苏皖两省2016—2018年小麦蚜虫发生等级,预报误差不超过1级,模型异地预报3年苏皖分省等级均基本正确,效果较好,这说明黄淮小麦蚜虫气象适宜度模型对江淮地区麦蚜发生等级有一定的指示性,可用于江淮麦蚜发生等级反演和预报。
3) 模型考察开始时间为上年冬季,结束时间最迟为5月上旬。结合不同区域小麦蚜虫盛发期的时间分布,对华北、黄淮和江淮可从上年冬季起,利用模型监测计算某省或某站小麦蚜虫发生发展气象适宜度指数,判别相应气象等级,进行预测预报,如4月中旬可发布灾害预警,当预报未来一旬或一个月气象条件利于虫害发生时,可将4月下旬至5月上旬气象关键因子在常年值基础上,根据预报波动幅度微调,用于开展麦蚜发生发展预报预警服务。
关键时段关键气象因子对小麦蚜虫发生发展具有决定性作用。以小麦蚜虫发生典型年2014年为例,该年黄淮海大部地区小麦穗期蚜虫偏重发生[36],主要原因是春季气温偏高,3月山东、河南平均气温分别为10.5℃和12.0℃(均高于常年值6.7℃ 和8.6℃),分别位居历史次高和首位;4月上旬河北、北京、天津、山西区域平均气温分别为14.7℃,15.3℃,15.2℃和12.6℃(均高于常年值11.6℃,11.4℃,12.0℃和9.7℃),京津冀气温均位居历史首位,气温显著偏高加快了生物发育进程,小麦发育期和蚜虫发生期均提前,山东小麦蚜虫始发期为当年2月20日,偏早约1个月。2014年关键时段关键气象因子中温度和湿度均在穗期优势种麦长管蚜发生发展的适宜气象条件范围内,高温少雨导致5月初河北中南部、山东中西部及半岛西部和河南北部部分地区出现旱情[37],温高光足,适宜的光、温、水条件和干旱叠加效应,促进了小麦穗期蚜虫的发生为害,造成2014年小麦产量损失达9.186×105t[36],这与前人研究结论吻合,水分缺乏可以提高干旱地区小麦蚜虫潜在的适应能力,小麦蚜虫将耗费更多的时间取食[38]。
关键气象因子的选取是从小麦蚜虫发生发展的适宜生理气象指标和区域气候特征考虑,对进行小麦蚜虫监测预报有较好指示意义。由于不同地区气候条件及小麦蚜虫发育进程和作物发育期本身的差异,小麦蚜虫的气象适宜度在时间和空间上存在一定差异,春季强降雨是影响麦长管蚜种群消减的关键因子之一[38-39],强降雨多不利于小麦蚜虫爆发,而降水偏少、无雨日数偏多及导致的干旱,有利于小麦蚜虫发生发展。全球气候变暖导致的温室效应和极端天气等衍生的干旱等灾害的频频发生[40-43],导致小麦蚜虫对逆境环境适应能力的显著提高[11,44-45]。这也是本研究针对不同区域建立不同的小麦蚜虫气象适宜度综合指标的意义所在。
本文仅考虑在一定虫源基数条件下地面气象条件对小麦蚜虫发生发展的适宜性,但小麦蚜虫是迁飞性害虫,其发生发展除受地面气象条件的重要影响外,还受大尺度环流背景、虫源基数、自身生物学特性、寄主作物及品种、天敌情况、耕作栽培方式、施肥和灌溉水平、田间管理措施及人为防治等综合因素的影响[3,8,24],北方小麦主产区小麦蚜虫发生发展气象适宜度等级仅对地面气象条件的满足程度进行了分析。实际应用中,可结合小麦蚜虫气象适宜度等级监测预报与其他因素综合考虑小麦蚜虫的发生情况,开展相关农业气象业务和服务。