基于WRF-Flexpart的一次褐飞虱回迁过程模拟研究

褐飞虱是影响亚洲地区水稻生产的迁飞性害虫,其远距离迁飞与种群分布会受到大气环流和局地天气条件的影响。采用WRF-Flexpart耦合模式、GIS空间分析和Vincenty方位角/距离估计等方法,选取2009年9月30日—10月7日发生在中国15个植保站的褐飞虱迁入过程,并根据模式输出、GIS分析和轨迹计算结果探讨了多种大气物理胁迫共存下大气动力场、降水、温度场和湿度场对褐飞虱秋季回迁过程及其降落虫量分布的影响。研究结果表明:(1)此次重大迁飞过程中,大气动力场、降水、温度场和湿度场在不同地区、不同时间段对褐飞虱种群的迁入和降落起了不同的作用,因而造成了迁入虫量分布的时空差异。(2)从各虫源地迁出的褐飞虱种群是向着温暖而湿润的地方迁飞的,迁飞方向主要受所经区域盛行风向的影响,迁飞距离由途径区的风速所决定。(3)当迁飞种群途径局地降水较多的区域时,初期迁入量并不多,经过短暂的一段时间后降虫量会突然增加。(4)在一定的温度范围内,大气湿度条件会影响褐飞虱的迁入虫量:褐飞虱迁入虫量在近地面温度高于26℃、相对湿度(925 hPa)大于70%且下沉气流较强的地区分布最多,而在近地面温度低于24℃、相对湿度(925 hPa)小于50%、无明显下沉气流的地区未发现有褐飞虱的迁入。     

上海地区2007年褐飞虱的后期迁入和虫源地的个例分析

褐飞虱后期迁入是指在8月下旬到9月上旬水稻生育后期出现的异常迁入。分析了2007年上海地区褐飞虱后期迁入种群的虫源性质、迁飞降落过程、大气背景场以及虫源地情况,结果表明:(1)2007年8月下旬到9月初上海地区出现的褐飞虱灯诱高峰为外地迁入峰,迁飞峰次分别为8月2223日,8月2731日。(2)8月2223日降虫主要受台风"圣帕"外围倒槽气流的影响;8月2731日褐飞虱的迁入主要受降雨过程的影响。(3)回推轨迹分析表明:2007年8月下旬到9月初上海地区后期迁入的虫源地主要分布在浙江东南沿海稻区(8月2223日),浙江中、北部(8月2731日),江南稻区北部包括皖南、赣北(8月29日)。(4)以怀宁为代表的皖南虫源地的顺推轨迹分析表明:8月下旬该区早中稻进入生育后期,食料条件恶化,褐飞虱长翅型成虫大量外迁,皖南稻区存在为上海提供虫源的对应关系。(5)虫源地的种植制度、虫源基数和合适的气象条件是导致后期迁入的重要因素。     

大气低温胁迫对中国褐飞虱年内初始迁入的影响

为了揭示大气低温胁迫对中国褐飞虱年内初始迁入的影响,更好地预警来自境外的褐飞虱早期迁入,通过统计2000-2017年中国华南、西南两个稻区褐飞虱年内的始见期和首次迁入峰日,逆推其迁飞轨迹和虫源区,分析大气温度场对迁飞过程产生的作用,比较了褐飞虱在我国不同稻区、不同年内初始迁入期受大气低温胁迫产生的作用差异。结果表明：（1）近年来褐飞虱初始迁入中国的时间提前,初始迁入华南稻区的时间比西南稻区早,华南稻区始见期提早可能与褐飞虱种群越冬北界北移有关。（2）西南稻区褐飞虱年内初始迁入的境外虫源主要来自缅甸,华南稻区年内初始迁入的境外虫源主要来自越南和老挝的中北部。（3）对盛行迁飞层的温场分布研究表明,在褐飞虱年内初始迁入过程中低温屏障发生的概率约为54.4%,迁入当晚降虫地的平均低温强度为13.45℃,平均降温幅度为1.88℃。其中低温胁迫在华南稻区表现更为显著（发生概率为58.3%,平均强度为13.18℃）,在始见期表现得更明显（发生概率为70.6%,平均强度为12.53℃）。     

2009年云南省白背飞虱早期迁入种群的虫源地范围与降落机制

对2009年云南省江城、金平、西畴、师宗、彝良5个站点4-5月出现的白背飞虱灯下高峰日进行轨迹分析,并结合部分田间调查数据、峰期对应的大气环流背景、峰日风温场等气象因子,明确2009年云南白背飞虱早期迁入种群的虫源地范围和降落机制。研究结果表明:(1) 2009年云南4-5月份前期迁入飞虱虫源主要来自缅甸,部分虫源来自老挝与越南,极少数虫源来自印度东部因帕尔地区。(2) 低温屏障、风场切变以及垂直气流扰动是造成白背飞虱集中降落的最主要气象因子。(3) 缅甸虫源区偏西气流盛行以及冬春稻种植面积扩大而导致的虫源基数增加是造成2009年云南白背飞虱大量迁入的根本原因。     

2007年广西早稻田褐飞虱发生动态及虫源分析

褐飞虱Nilaparvata lugens(Stal)是我国水稻上的重要害虫,近年来在全国连续大暴发,给我国水稻生产造成严重的损失。广西是境外虫源迁入我国的第一站,分析其早稻田褐飞虱的种群动态和迁飞规律,对全国褐飞虱的预测预报和防治工作意义重大。采用田间系统调查与雌虫卵巢系统解剖的方法研究了2007年南宁地区早稻田的褐飞虱发生动态和世代种群特征,分析全区灯诱数据确定迁入高峰期,并运用美国NOAA网站的HYSPLIT平台对2007年广西褐飞虱早期的6段迁入高峰期进行了数值模拟和虫源地分析。结果表明:(1)2007年广西南宁地区褐飞虱在晚熟品种田的发生虫量明显高于早熟品种田,发生、危害期明显延长,为我国北方其他稻区的褐飞虱的迁入积累了充足的虫源基数。(2)确定了南宁地区的早稻田褐飞虱各发生世代的虫源性质。第2代属于大部迁入型;第3代前期属于本地繁殖,部分迁入型,后期属部分迁出型;第4代属于大部分迁出型。(3)4月下旬-5月上旬,虫源地主要分布在越南、老挝中部的地区(18°-20°N);5月中下旬,虫源地主要分布在越南、老挝较北的地区(20°-22°N),部分来自我国海南省和雷州半岛;而6月上旬桂北地区大规模迁入的虫源地主要分布在越南东北部、广西西南部;低空急流为广西褐飞虱大规模迁入提供了运载气流,而当时的降雨是迫使褐飞虱集中降落的主要原因。     

桂东北稻区第七代褐飞虱迁飞规律及虫源分析

褐飞虱是我国水稻上的一种重要迁飞性害虫,在桂东北稻区不能越冬。近年来,作为次害代的第7代褐飞虱在桂东北稻区频繁暴发,桂东北稻区作为褐飞虱南北往返迁飞的重要通道,研究阐明该稻区第7代褐飞虱的迁飞规律及虫源分布,对第7代褐飞虱的早期预警和有效防控意义重大。利用灯光诱捕的方法研究了2007—2010年桂东北稻区褐飞虱的迁飞动态,并运用大气质点轨迹分析平台HYSPLIT(Hybrid Single-Particle Lagrangian Integrated Trajectory)和气象图形处理软件GrADS(Grid Analysis and Display System)对褐飞虱迁入、迁出峰次进行轨迹分析和气象背景分析。结果表明:(1)桂东北稻区第7代褐飞虱每年都会出现2—3次迁飞高峰,虫量高峰期主要集中在9月底至10月上旬。(2)迁入虫源主要来自江西中西部和湖南中东部稻区,从本地起飞的部分褐飞虱飞行24 h或36 h后能到达广西南部沿海、越南北部和海南等越冬区。(3)925 hPa的温度高则褐飞虱虫量高;盛行东北气流较多的年份,褐飞虱发生严重;微弱的下沉气流和无雨天气有利于褐飞虱的迁飞。     

复杂地形下褐飞虱迁飞的数值模拟:个例研究

褐飞虱是影响我国水稻生产的重要迁飞性害虫之一,它的爆发严重影响水稻的生长,并诱发水稻病害,导致水稻减产。研究复杂地形条件下褐飞虱的迁出虫源地、空中迁飞轨迹、降落区,探明复杂地形对褐飞虱迁飞的影响机制,对害虫测报与防治、农业防灾减灾和保障我国粮食安全具有重要的理论意义和实践价值。为了分析覆盖复杂地形的大气动力场、温度场、湿度场对褐飞虱迁入及降落分布的影响,采用WRF-Flexpart耦合轨迹计算模式和GIS空间分析等方法,对2008年9月30日—10月3日发生在广东曲江、肇庆、梅县三站和湖北宜昌、安徽东至、江西吉安三站的褐飞虱迁入和降落过程进行了数值模拟,以揭示复杂地形条件下褐飞虱种群降落和密度分布的时空变化结构。(1)广东曲江、肇庆、梅县三站的后向轨迹模拟结果显示,迁入三站的虫源均来自该站点的西北方向。地形较高且复杂多变时,褐飞虱难以穿越且迁飞距离较短、方向多变。(2)湖北宜昌、安徽东至、江西吉安三站的前向轨迹模拟结果显示,当站点附近的山脉较低且有山谷通道,褐飞虱沿山脉顺风方向迁飞,且迁飞距离较远。当山脉地势较高且没有明显的山谷,则褐飞虱遇到山脉阻挡而转向造成种群滞留。(3)褐飞虱迁飞种群的密度沿山脉走向呈带状分布,山坡较为陡峭、断崖显著时,向远离山体的方向迁飞。若山脉由多个山岭构成,则褐飞虱可从其峡谷穿越,密度分布较为分散。(4)强水平气流有利于褐飞虱的远距离迁飞,下沉气流对褐飞虱降落起着重要的作用,当有强下沉气流且气温较高时,有利于褐飞虱种群的降落,并聚集形成高密度迁入区。(5)垂直方向上,在一定的温度范围内,褐飞虱趋向于在暖层中迁飞。褐飞虱密度沿河谷地带呈带状分布且密度高值区多分布在较温暖的地区。秋季,褐飞虱降落区多分布在相对湿度50%左右的区域。模拟的褐飞虱迁飞轨迹、迁飞方位角和迁飞距离等与实际发生的褐飞虱迁飞路径和降虫区之间偏差较小,该模拟方法较大程度地提高了我国迁飞性害虫的业务预报水平。未来拟提高测报褐飞虱虫情数据的时空分辨率,以为得到较高准度和精度的模拟结果。     

云南勐海早稻白背飞虱种群动态及虫源地分析

【目的】白背飞虱Sogatella furcifera(Horváth)是勐海稻区重要的迁飞性害虫,对当地水稻的生产造成严重威胁。本试验旨在研究勐海早稻白背飞虱种群发生动态与虫源地位置,明确迁飞种群发生规律,为实施害虫的异地预测预报提供科学依据。【方法】利用HYSPLIT轨迹分析模型和Gr ADS软件,分析2000-2016年灯诱峰期迁入种群的虫源地分布和气象背景情况;通过田间系统调查和雌成虫卵巢系统解剖,研究2013年勐海稻区白背飞虱的种群动态以及各发生世代的虫源性质。【结果】勐海稻区4月的虫源地主要集中在缅甸中部,部分来自于缅甸南部和金三角地区,5月的虫源地集中在缅甸东部,泰国、老挝、越南的北部也能提供部分虫源;灯诱峰期的天气背景主要为风切变、下沉气流、低温屏障和降水;第2代白背飞虱为危害勐海早稻的主害代,时间为4月中旬至5月中旬。【结论】明确了勐海早稻迁入种群的虫源地分布及主害代发生时期,为白背飞虱在当地的防治和异地预测预报提供了技术支撑。     

褐飞虱的降落机制

应用MM5中尺度数值模拟模式和GRADS气象图形软件,对1999年6月22日～7月2日一次褐飞虱迁飞降落过程及其大气背景场进行了个例研究,并侧重对灵川、东至、徽州3个站点的降虫机制做了分析。结果表明：（1）此次迁飞过程有两个明显的降虫区：天柱．灵川．贺县一线范围内（降虫区Ⅰ）和三都-天柱-黔阳-修水-东至-徽州-奉贤-线（降虫区Ⅱ）;（2）低空急流为褐飞虱北迁提供了运载气流,降虫区位于降水区边缘;（3）地形降水胁迫害虫降落是降虫区Ⅰ褐飞虱降落的主要原因,梅雨锋降水导致了褐飞虱集中降落在降虫区Ⅱ;（4）虫源的空间分布差异和不同的天气条件导致了东至、徽州地区不同的降虫分布。     

滇东南地区白背飞虱早期迁入种群的虫源地与降落机制

【目的】滇东南地区位置特殊,地形复杂,是白背飞虱Sogatella furcifera(Horváth)迁入云南的重要通道之一,早期迁入种群虫源及其降落机制的阐述将为白背飞虱的异地预测提供理论依据。【方法】(1)用HYSPLIT对滇东南地区富宁、麻栗坡、丘北和广南4个站点2008—2012年4—5月白背飞虱迁入高峰日做轨迹分析以阐明其虫源区分布。(2)用Gr ADS分析NCEP的再分析格点数据以形成各迁入峰期的风温场,阐释迁入峰得以形成的气象背景和降落机制。【结果】(1)滇东南白背飞虱早期迁入种群虫源地主要分布在越南和老挝北部,部分来自于缅甸中北部及掸邦高原、泰国北部、滇西和滇西南。其中,4月份以缅甸和金三角为主,5月份则主要来自越南和老挝北部。(2)影响白背飞虱降落的因子中,风切变、低温屏障、降雨和下沉气流分别占75.9%、56.5%、31.8%和25.2%。【结论】滇东南白背飞虱早期迁入虫源主要来自中南半岛,风切变、低温屏障、降雨和下沉气流四者共同作用形成了白背飞虱在滇东南的频繁迁入。     

Simulation of the Long Range Migration of Brown Planthopper,Nilaparvata lugens (Stal), by Using Boundary Layer Atmospheric Model and the Geographic Information System

Brown planthopper (Nilaparvata Lugens) is a migrant pest which can migrate from tropical area into subtropical and temperate area every year and caused a lot of damage in rice production. To understand the migration process of this pest, the Boundary Layer Atmospheric (BLAYER) model and Geographic Information System (GIS) were used to analyze the migration waves from June to July in 1997, 1998 and 1999 in South Korea. The simulation results showed: 1) Each migration wave had different mass distribution and immigration area at different time; 2) The vertical air current value distributed in BPH taking-off and landing area was about several centimeters per second, which is lower than BPH flight ability. 3) The trajectory route showed that BPH had different migration routes emigrated from different source areas according to the weather system. The main source of BPH immigrated in South Korea was from the East part of Guangdong Province and South East part of Fujian Province, People's Republic of China. 4) Comparing with the different migration heights (733 m, 1,348 m and 1,963 m above ground level), BPH mass distributed more northwestward in low height than in high height. 5) BPH mass movement also gave another evidence to identify the light trap data in late July.     

中南半岛前期异常气候条件对中国南方稻区褐飞虱灾变性迁入的影响及其预测模型

近30多年来,气候变化对中国褐飞虱的灾变性迁入带来了明显的影响,为了进一步了解虫源地的异常气候变化对我国褐飞虱迁入量的影响,利用1980—2016年中国各植保站提供的褐飞虱虫情资料及同期美国国家环境预测中心(NCEP)提供的全球气象再分析资料,分析了中国褐飞虱境外主要虫源地中南半岛前期异常气候条件对中国南方稻区褐飞虱发生程度的影响,并对褐飞虱发生等级与影响其迁飞的气象因子进行了相关性分析,筛选出关键预报因子,分别应用支持向量机(SVM)、BP神经网络和多元回归分析3种方法对代表站点褐飞虱年发生等级进行了预测,并比较了3种预测模型的优劣。结果表明:(1)中南半岛气候异常区主要分布在北部,异常气候的发生次数在中南半岛呈现出北高南低的特征,并从北向南呈环状递减。(2)中南半岛前期温度偏高(暖冬、暖春)、相对湿度偏大(湿冬、湿春),易引起褐飞虱在中国南方稻区的偏重以上发生;若中南半岛前期气候偏冷(冷冬、冷春)、偏干(干冬、干春),则常导致褐飞虱在中国南方稻区的偏轻以下发生。(3)通过比较3种模型的历史回代率和预测准确率,发现3种模型对褐飞虱的发生程度均有一定的预测能力,其中SVM模型的预测效果最好,BP神经网络次之,多元线性回归模型最差,表明SVM更加适用于生产实际中的褐飞虱发生程度预测。     

七角星蜡蚧在云南省的分布与危害

[目的]七角星蜡蚧(半翅目:蚧次目:蚧科:蜡蚧属)是一种外来入侵害虫,自2013年在云南省西双版纳州勐仑镇发现以来,中国大陆再无系统性的调查。近年来作者在云南省收集蚧虫期间,在多地发现七角星蜡蚧为害,现对其在云南省的发生情况进行报道,以期为其监测和防治提供参考。[方法]2016年10月—2021年9月,对云南省15州(市)的64个县(区市)进行了系统调查,详细研究了七角星蜡蚧的分布、寄主及一些生物学习性。[结果]七角星蜡蚧在云南省西双版纳州的景洪市、勐腊县(勐仑镇和磨憨镇)、勐海县(打洛镇),普洱市的宁洱县,德宏州的盈江县(那邦镇),红河州的河口县、金平县(勐拉乡和者米乡)4州(市)7县(市)9地有分布,新增8个分布点。除宁洱县外,其余发生地均位于边境线上,表明其可能为境外多点入侵。在云南该虫危害寄主植物共4科4属5种,包括杧果和莲雾2种经济树种。七角星蜡蚧多寄生在叶片背面,沿叶脉分布,1年发生多代。在河口县的龙船花上,该虫每雌产卵38~50粒,虫口密度达185~250头·叶^-1。[结论]七角星蜡蚧在云南省零星分布,尚未造成严重损害,但有向内扩散的趋势,应加以警惕。     

中国囊萼花属一新记录种——茎花囊萼花

首次报道了囊萼花属(Cyrtandromoea Zoll.)的1个中国新记录种——茎花囊萼花(Cyrtandromoea grandis Ridl.),该种原记载产于泰国、马来半岛和苏门答腊岛,2019年在中国云南省勐海县布朗山地区发现有分布。茎花囊萼花与囊萼花(Cyrtandromoea grandiflora C.B.Clarke)最为相似,但前者聚伞花序生于老茎上,每一花序3～18朵花,二者容易区别。     

基于两种轨迹模型的褐飞虱迁飞轨迹比较研究

异地预测是迁飞性害虫发生预测的重要内容,迁飞轨迹模拟和预测是能较好地反映害虫迁飞时空动态的一种异地预测方法。褐飞虱作为我国水稻生产上的一种重要迁飞性害虫,其迁飞轨迹的准确预报,可为其灾变预警和有效防控提供科学依据。为了比较选择一些准确性好、分辨率高、易于推广应用的害虫迁飞轨迹模型,选取2006年7月初发生在湖南省洪江市的一次褐飞虱重大北迁过程作为典型个例,运用中尺度天气研究和预报模式WRF,结合NCEP气象再分析数据,利用HYSPLIT和FLEXPART两种轨迹计算模式对褐飞虱迁飞轨迹进行模拟,并验证模型模拟和计算的准确度和精确度。研究结果表明:(1)WRF-HYSPLIT和WRF-FLEXPART两种轨迹计算模式在虫源地、迁飞路径(迁飞方位角和走向)、迁飞高度、迁飞速率和迁飞距离计算上总体趋势一致,但存在一定的差异,后者的起伏变化大于前者。(2)尽管两种耦合模式在调用WRF模式输出的预报场物理变量方面有很多相同之处,但WRF-FLEXPART耦合模式在运行计算过程中比WRF-HYSPLIT耦合模式多考虑了对流参数、地表胁迫和各种地形参数,因而能更全面地反映中尺度天气过程(特别是对流性天气过程)对昆虫起飞、空中飞行和降落的动力作用,也能更真实地反映地表物理过程、大气湍流结构和地形起伏对褐飞虱种群迁飞的影响。(3)从褐飞虱种群对生境和取食条件选择上看,两种模式模拟的各高度迁入种群的虫源区、迁飞路径和降落地都是合理的、准确的。但从褐飞虱迁出、空中飞行和降落所处的三维流场来看,WRF-FLEXPART模式轨迹走向与盛行气流方向的吻合度要明显高于WRF-HYSPLIT模式。(4)两种模式均可作为业务工具在迁飞性害虫测报中推广应用。     

江苏沿江和江淮区褐飞虱前期迁入量与太平洋海温场的遥相关及其可能机制

为筛选出有效的长期预测因子,对江苏通州和江苏淮安褐飞虱前期迁入量与太平洋海温场遥相关时空分布分别进行了相关分析及其稳定性检验。结果表明： 与褐飞虱前期迁入量连续稳定显著相关的海温区,在空间上主要分布在中太平洋和南太平洋（72.4%）,在时间上主要分布在前两年和前一年（82.7%）。分别分析褐飞虱前期迁入量与当年4-8月上旬各旬降水量、平均温度、当年1-8月500 hPa西太平洋副高月平均指数之间的两两相关关系,并比较这三者与前期太平洋海温场遥相关时空分布的异同。根据对两地的个例分析表明：前一冬春季赤道中东太平洋海温的变化引起当年6月副高面积指数的变化,影响江苏通州当年6月下旬降水,从而影响江苏通州褐飞虱的前期迁入量;前两年春夏季赤道中东太平洋海温的变化引起当年5月副高北界位置的变化,影响江苏淮安当年7月上旬降水,从而影响江苏淮安褐飞虱的前期迁入量。在此基础上,对褐飞虱前期迁入量与太平洋海温场遥相关的可能机制提出“海温→大气环流（副高）→气候→褐飞虱前期迁入"的一般假想模式。     

褐飞虱持续暴发成因——以江西泰和县为例

为了明确泰和县近年来晚稻田褐飞虱Nilaparvata lugens(Stl)连年暴发的规律。通过分析泰和县2008—2013年褐飞虱的田间调查及灯诱数据,并结合气象条件,研究了泰和县晚稻田褐飞虱暴发的主要影响因素。结果表明:(1)危害早稻田的褐飞虱以2、3代为主,其田间种群数量的多少主要受春季迁入虫量的限制,本地早稻收割后的滞留虫源对晚稻种群数量的影响很小。(2)危害晚稻田的褐飞虱以5、6代为主,后期迁入虫量、秋季温度以及晚稻田初期外来迁入虫量是影响晚稻田后期暴发的关键因子。(3)不同年份间晚稻田褐飞虱暴发类型不同。2008和2009年属后期迁入虫量与"暖秋"共同影响类型;2010年是"暖秋"影响类型;2011和2012年是后期迁入虫量影响类型;2013年是晚稻田早期迁入虫量积累影响类型。后期迁入虫量、秋季温度以及晚稻田早期迁入虫量是导致泰和县晚稻田褐飞虱连年暴发的主要原因。这些结果为褐飞虱暴发规律的阐明及可持续治理提供了理论依据和参考价值。     

气候变暖对中国褐飞虱越冬北界的影响

基于GIS技术,利用1951—2010年中国南方地区289个气象站点观测资料和褐飞虱历史越冬考察资料,通过采用典型年褐飞虱越冬实况调查资料进行检验,对前人研究的褐飞虱越冬北界指标进行优选,最终选择1月10℃等温线作为中国褐飞虱越冬北界位置确定的指标,分析了最近60年来气候变暖背景下中国褐飞虱越冬北界年代际变化以及前后30年的变化。结果表明:在气候变暖背景下,由于暖冬年份频次的增加,使得每隔10年、30年褐飞虱越冬北界平均界限均不同程度北移,与冷冬年相比,暖冬年褐飞虱越冬北界位置可向北推移2～3个纬度;年代际变化以20世纪80年代最为明显,80年代之后的30年较之前的30年北移显著;空间北移较明显的地区为云南中部和福建南部,其中云南境内最大位移处60年约向北移动了近90km。     

Teleconnection between the early immigration of brown planthopper (Nilaparvata lugens Stål) and ENSO indices: implication for its medium- and long-term forecast

Teleconnection between the early immigration of brown planthopper (BPH) and El Nino/Southern Oscillation (ENSO) indices from January of two years previously to the current June was investigated to make long-term forecast. The teleconnection results were as follows: ENSO indices which were significantly correlated with the early immigration of BPH were primarily sea surface temperature anomalies (SSTA) in N3, N4 and N3.4 regions, accounting for 71.8% of the total. Significant ENSO indices from two years and one year before the immigration events had a proportion of about 84%, while those in the current year only accounted for 16.7%. There was significantly negative correlation between the early immigration of BPH and SSTA in each Nino region from two years before to the previous spring, whereas there was significantly positive correlation between these two factors during the period from the previous winter to the current spring. The significant correlation between the early immigration of BPH and SSTA in each Nino region in the last summer and autumn did not show any obvious tendencies. The relationship between the southern oscillation index (SOI) and the early immigration of BPH was opposite to that between the immigration and SSTA in each Nino region. The above mentioned significant ENSO indices were used as key factors to build forecasting models for the early immigration of BPH by step-wise multiple linear regression analysis. Finally, 12 integrated forecasting models were obtained, which could make predictions 3–27 months ahead and had a predictive accuracy of 88.9%.