环渤海夏蝗发生程度气象集成预报方法
An Integrated Bio meteorological Forecasting Method for the Occurrence Level of Locust Around the Bohai Sea of China
-
摘要: 东亚飞蝗发生程度与气象条件关系密切,该文使用1980—2008年的飞蝗资料和气象资料,选择了环渤海4种典型(沿海、水库、洼淀、内涝)蝗区,利用秩相关系数法筛选影响飞蝗发生程度的气象因子,确定了影响不同蝗区夏蝗发生程度的气象因子指标集。利用权重修正气象距离法,建立了气象距离指标预报模型;依据飞蝗自身生物学特性,建立了生物学预报模型。在此基础上,基于蝗虫生物学和气象条件影响的共同作用,建立了夏蝗发生程度的生物气象集成长期预报模型。结果表明:在环渤海飞蝗区域,不同类型蝗区影响气象因子有一定差别,集成预报模型趋势预报准确率高于其他模型。Abstract: Locust is an important loss causing pest around the Bohai Sea of China, and its occurrence degree is close related with weather and climate conditions. In order to study the relationship objectively, four typical locust areas around the Bohai Sea are chosen as experimental areas, which are coast, reservoir, depression shallow lake and water logging areas, respectively. Based on the observed climate and locust data from 1980 to 2008, Spearman order correlation method is used to analyze the meteorological factors affecting the occurrence degree of locust. The results show that meteorological factors have accumulated effects on locust. The temperature in July and August has a significant effect on summer locust of the next year in four locust areas, and higher temperature is more favorable for locust growth and reproduction. October is the key time for the egg life of autumn locust in three areas (reservoir, depression shallow lake and water logging areas), and sufficient rainfall is favorable for the development of locust egg, increasing the occurrence degree of summer locust in the next year. In addition, air temperature in winter and spring and rainfall affect the extent of summer locust in four areas before locust eggs hatch and come out of soil. The historical modeled accuracy are 81%—93% and 78%—89% by the Euclidean distance model modified with the weights of meteorological factors and biological model based on locust bio characteristics for predicting the locust extent, respectively. The extended forecast result of last two years is fairly accurate, i.e., one level difference for one area of the former model and one level difference for two areas of the latter model, and correct in all other areas. Then a comprehensive model is established by integrating the meteorological and biological models to forecast the locust occurrence extent and its historical modeled accuracy are 85%—96%. There is only one level difference in one area for the two years extended forecast, showing that the accuracy of integrated model is better than the single models.
-
引言
环渤海蝗区是我国主要蝗区之一,也是东亚飞蝗发生基地,位于我国的东部沿海地区,包括河北省、天津等蝗区。环渤海东亚飞蝗在20世纪90年代后期至21世纪初期发生程度较重,如1998-2002年环渤海先后多次出现高密度蝗群,密度达1000头/m2以上,给当地农业生产带来严重威胁。蝗虫的发生、发展是自身生物特性与其生存环境相互作用的结果,其中气象条件对蝗虫的影响最直接。大量研究表明,在其他条件 (如适生环境、虫源基数) 具备的情况下,气象条件往往成为影响蝗虫发生程度的决定性因素[1-3]。环渤海蝗区为季风气候区,天气气候环境条件对蝗虫灾害发生作用明显[4-6]。
气象因子对飞蝗发生的影响引起了诸多学者的关注[7-9],其中,既有分析气象条件与蝗虫发生程度的定性关系,如温度对于蝗虫地理分布的影响[10],还有分析二者之间的定量关系,并建立预测模型[11-13],如前期干旱、温度偏高对飞蝗发生程度的定量影响[13-14],预报蝗虫发生的气候定量指标[15]。此外,还有从生态学观点出发,建立基于湿热系数的蝗卵孵化死亡率的估计模型[16]。这些研究从不同时间、空间尺度描述了蝗虫发生及与气象条件的关系,并进行定性分析、定量预报,为了解蝗虫的发生规律及预测、防治起到重要作用。蝗虫发生、发展程度除受气象等环境因素的影响外,其自身生物特性也具有一定作用,如影响飞蝗发生程度的每年虫源基数与诸多因素有关,包括飞蝗的防治力度,飞蝗生存条件的生态环境变化如植被长势好坏、天敌的多少等,这些因素共同对飞蝗发生程度起重要作用。因此,气象条件不能完全表述对飞蝗发生程度的影响规律,以此为基础的飞蝗发生程度的气象分析与预报模型在实际应用中常表现出不稳定性。
基于前人对蝗虫发生、发展程度预报的研究,在不考虑异地虫源的前提下 (因蝗虫多年来预防和治理效果显著,异地迁飞现象很少),本文从环境条件与蝗虫生物学特性出发,综合分析气象和飞蝗自身生物潜能对飞蝗发生程度的共同影响与作用,采用数理统计方法,建立相对稳定的飞蝗发生程度的集成预报模型,开展环渤海飞蝗发生程度预报,为蝗虫防治提供及时的预警信息。
1. 资料来源
1.1 飞蝗资料
选择环渤海典型蝗区黄骅、安新、平山、大城,分别代表沿海、洼淀、水库、内涝4个类型蝗区[17],其中属于沿海蝗区的黄骅位于渤海沿岸,海拔在5 m以下,地势低平,地广人稀,耕作粗放,常年杂草丛生,是东亚飞蝗的发生基地;属于洼淀蝗区的安新位于河北省中部,海拔10 m左右,主要分布在白洋淀周围,其洼淀周边多为苇荒地,形成了飞蝗的天然孳生场所,是东亚飞蝗的常发区;属于水库蝗区的平山位于太行山边沿,海拔120m以上,蝗区主要分布在黄壁庄水库周边无人耕作地,因水库水位不稳定,飞蝗的适生面积变化较大;属于内涝蝗区的大城位于河北省广阔平原的低洼农区,海拔5m左右,飞蝗主要发生在夹荒、渠堤、埝埂、沿河小片洼地等特殊环境。环渤海飞蝗每年发生两代,即夏蝗和秋蝗,本文以夏蝗为研究对象。资料来源于上述4县植保站,时间为1980-2008年,其中1980-2006年资料用于建模,2007-2008年资料用于延伸预报。
1.2 气象资料
选用以上4个蝗区对应的县气象站气象资料,包括月和旬的气温、降水、日照、空气相对湿度、土壤温度等要素,资料年代与夏蝗观测资料一致。
2. 夏蝗发生程度气象预报模型
2.1 因子筛选
Spearman秩相关系数分析常用于以等级尺度观测的资料,用其进行影响夏蝗的气象因子筛选,计算公式为
(1) 式 (1) 中,RS为秩相关系数,N为成对数据 (X,Y) 的个数,D表示相应X,Y值所对应秩的差值。计算时将每组数据按照大小、重要性等方法排序,从1到N确定其秩,用数据的秩代替数据的值来进行相关性分析。
通过对4个类型蝗区气象要素和夏蝗发生程度的秩相关系数计算,筛选关键气象因子 (表 1)。从表 1可见,被选入的因子中,不同蝗区有共同的气象因子,分别是前一年7月、8月气温,因前一年7-8月高温可促进秋蝗生长发育和繁殖,为来年夏蝗的发生提供虫源,这是影响夏蝗发生的关键因子。前一年10月降水偏多是诱发3个蝗区 (安新、平山、大城) 夏蝗大发生的共同因子。10月是秋蝗产卵的关键时期,土壤水分对其影响很大。研究表明,土壤的温、湿度状况对蝗卵孵化与否及死亡率高低有直接影响[18],只要蝗卵从环境中吸收到足够的水分,就可以完成发育。
表 1 影响环渤海不同蝗区的主要气象因子及秩相关系数Table 1. ajor meteorological factors and rank correlation coefficient affecting the locust-infested areas around the Bohai Sea黄骅因地势较低,10月降水偏多易形成长时间积水,部分蝗卵受淹死亡,孵化率降低,因此黄骅蝗区前一年10月降水不是夏蝗发生程度较重的主要因子。其他影响因子因蝗区不同而有一定差别,其中冬季气温或地温也是影响夏蝗发生程度的主要因子,同时春季气温和蝗蝻出土前的降水对夏蝗有显著影响,前期气象因子对夏蝗具有明显滞后和累积效应。
2.2 权重修正的气象距离模型
采用文献[18-19]权重修正的气象因子集向量距离方法,其计算方法为:将采用Spearman秩相关系数筛选不同蝗区的气象因子分别进行有序组合,建立气象因子向量集
(2) 建立两种极端状况的参照向量,即极不利于和极有利于夏蝗发生的气象条件向量
(3) (4) 当Riy>0时,
当Riy<0时,
上式中Riy(秩相关系数) 为第i个因子与夏蝗发生程度的相关系数。式 (2) (4) 中,xi为第i个因子的实际值;σ(xi),avg (xi),max (xi),min (xi) 分别为第i个因子的均方差、平均值、历史最大值和历史最小值,均为1980-2006年气象资料。
计算实际气象向量与两个参照向量Φ01,Φ02之间的欧氏距离
(5) (6) 经过归一化处理,则气象条件有利于夏蝗发生的距离比指标
(7) 式 (5),(6),(7) 中,t′i为第i个因子的标准化值,d1,d2分别表示实际向量与易发生和与不易发生夏蝗的气象向量之间的距离,dr值越大,表示气象条件有利于夏蝗的发生,反之气象条件不利于夏蝗的发生,Pi为第i因子在该距离指标中的贡献权重,计算方法为
(8) (9) (10) 式 (8),(9),(10) 中,Riy为第i个因子与夏蝗发生程度的相关系数 (秩相关系数),Rij为第i,j因子之间的相关系数,Rjy为第j个因子与夏蝗发生程度的相关系数,Ui表示第i个因子对夏蝗发生程度的独立影响程度,U为整个因子集的影响程度,Pi为第i个因子的贡献度,即第i个因子的影响在气象条件总影响中所占的比例,各因子贡献度之和
,当Pi<0时,表示第i个因子对夏蝗的影响已被其他因子的交叉相关所包含,在模式中自动剔除,减少因子之间的线性相关性。利用最小误差统计分析方法[19],经历史反查,建立了气象条件向量的距离比分级指标集 (见表 2),当计算结果dr值大于第I级分级指标dI,并小于dr (I+1)时,预报夏蝗发生程度为第I级。据此建立夏蝗发生程度的气象向量距离预报模型。
表 2 基于气象距离比的环渤海不同蝗区夏蝗发生程度等级指标Table 2. he indexes of meteorological Euclidean distance ratio for predicting the locust occurrence extent表 2中数据为夏蝗发生程度为1,2,3,4,5级对应的dr指标值,夏蝗发生程度同植保站统一指标。
2.3 基于夏蝗自身生物潜能发生程度预报模型
飞蝗本身具有极强的适生和繁殖能力,夏蝗发生与前期的发生状态有密切关系,具有生物延续性,其每年秋蝗的残蝗密度是下一代夏蝗发生的虫源基数,对夏蝗的发生程度起着重要作用[4]。用直线回归方法建立夏蝗发生程度的生物潜能预报模型:y=ax+b式中,y为发生程度预报值;x为上一代秋蝗残蝗密度。所建方程 (见表 3) 均通过0.01的显著性检验。
表 3 环渤海不同蝗区夏蝗自身潜能预报模型Table 3. The self-biotic predicting model of locust around the Bohai Sea2.4 夏蝗发生程度的生物-气象集成预报模型
将气象条件影响与夏蝗自身生物潜能结合[16],并进行权重集成,建立夏蝗生物-气象综合预报模型
(11) 式 (11) 中,y为夏蝗发生程度的综合预报结果;yi为i种模型的预报值;b′i为第i种模型的集成权重,模型的拟合率越高,其权重越大,本文以模型的历史拟合率表示
(12) 式 (12) 中,bi为第i种模型的历史拟合率;n为集成模型个数。
表 4、表 5为利用不同模型进行历史拟合及延伸预报的结果。其中,两年延伸预报中有1站相差1级,经过实地调查分析,其原因可能是2007年秋季降水偏多,一些地块不能及时耕作,成为撂荒地,使夏蝗发生程度偏重。总体上,综合预报准确率明显提高,并且使预报方法更趋向于机理化,可供飞蝗发生区进行业务化应用。
表 4 不同预报模型的夏蝗发生程度历史拟合率比较Table 4. Historical fitting accuaracy of three types of predicting models for the locust extent表 5 不同预报模型的夏蝗发生程度延伸预报与实际情况比较Table 5. Forecasting effects for different models of locust3种预报方法的时效略有不同,气象距离模型提前1个月,自身生物潜能可提前3个月,综合预报模型提前1个月。因此,在业务中可将3种模型综合应用。
3. 结论与讨论
1) 应用秩相关系数进行气象因子普查,发现前期气象因子具有滞后和累积效应,4个不同蝗区的前一年秋蝗发生的关键期7-8月气温对夏蝗发生程度有显著影响,气温越高越有利于秋蝗发育和繁殖,对来年夏蝗的发生有明显促进作用;3个蝗区的前一年10月为秋蝗产卵吸水关键期,决定蝗卵能否成活,降水偏多有利于蝗卵完成发育,对来年夏蝗发生有利。冬季和春季的气温、蝗卵孵化出土前的降水对不同蝗区夏蝗发生程度均有一定影响。
2) 环渤海4个蝗区的气象距离指标模型预报拟合率为81%~93%,两年延伸预报 (2007年和2008年) 结果有1蝗区站相差1级,其他站趋势准确;基于自身生物潜能模型,预报拟合率为78%~89%,两年延伸预报结果有两蝗区站相差1级,其他站趋势准确;基于夏蝗群体动态演变的生物特性与气象条件结合的综合集成预报模型,预报拟合率为85%~96%,两年延伸预报结果2008年有1站与实际相差1级。总体上集成预报模型效果优于单一预报模型。
本文根据历史资料,利用秩相关系数法对影响飞蝗的关键气象因子进行筛选,虽然有一定的生物学意义,但若对关键气象因子进行人工气候箱的试验研究[20]来验证将更有说服力;对夏蝗自身潜能预报利用的是前一年秋蝗残蝗基数做为因变量,但夏蝗发生情况还要经历冬季和早春时段,这一时段外界环境条件对其尚有一定的消长作用,因此,还需根据冬季蝗卵越冬死亡率和春季蝗卵孵化率对预报进行订正。本文对夏蝗发生程度进行了试预报,在飞蝗监测预测过程中还要考虑具体地点的实际情况,今后应尝试发生面积的预报及遥感发生地点的监测,将遥感监测与发生程度和面积预报有机结合,以达到飞蝗有效防治目的。
致谢: 中国气象科学研究院霍治国研究员在本文撰写过程中提出了宝贵意见,在此特表感谢。 -
表 1 影响环渤海不同蝗区的主要气象因子及秩相关系数
Table 1 ajor meteorological factors and rank correlation coefficient affecting the locust-infested areas around the Bohai Sea
表 2 基于气象距离比的环渤海不同蝗区夏蝗发生程度等级指标
Table 2 he indexes of meteorological Euclidean distance ratio for predicting the locust occurrence extent
表 3 环渤海不同蝗区夏蝗自身潜能预报模型
Table 3 The self-biotic predicting model of locust around the Bohai Sea
表 4 不同预报模型的夏蝗发生程度历史拟合率比较
Table 4 Historical fitting accuaracy of three types of predicting models for the locust extent
表 5 不同预报模型的夏蝗发生程度延伸预报与实际情况比较
Table 5 Forecasting effects for different models of locust
-
陈怀亮,张弘,李有.农作物病虫害发生发展气象条件及预报方法研究综述.中国农业气象,2007,28(2):212-216. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGNY200702023.htm 姜燕,霍治国,李世奎,等.东亚飞蝗发生的生态环境辩识及监测预测.生态学杂志,2007,26(5):737-742. http://www.cnki.com.cn/Article/CJFDTOTAL-STXZ200705022.htm 帕尔哈提·努尔.气候条件对蝗害形成中的影响.新疆畜牧业,2008(2):51-52. http://www.cnki.com.cn/Article/CJFDTOTAL-XJCM200802032.htm 任春光,张书敏,唐铁朝,等.白洋淀东亚飞蝗发生日趋严重原因分析.植物保护,2002,28(1):35-37. http://www.cnki.com.cn/Article/CJFDTOTAL-ZWBH200201013.htm 季荣,谢宝瑜,李典谟.河北省南大港农场2002年夏蝗发生特点及原因分析.昆虫知识,2002,39(6):430-432. http://www.cnki.com.cn/Article/CJFDTOTAL-KCZS200206008.htm 任春光,唐铁朝.河北省东亚飞蝗发生动态及未来灾变趋势分析.昆虫知识,2003,40(1):80-82. http://www.cnki.com.cn/Article/CJFDTOTAL-KCZS200301020.htm 马世骏.东亚飞蝗发生与气候条件的关系.昆虫学报,1958,8(1):31-40. 高素华,刘玲,郭安红.高温干旱与蝗虫发生研究.自然灾害学报,2005(3):169-174. 杨洪升,季荣,王婷,等.新疆蝗虫发生的大气环流背景及长期预测.生态学杂志,2008,27(2):218-222. http://www.cnki.com.cn/Article/CJFDTOTAL-STXZ200802015.htm 封传红,王思忠,蒋凡,等.温度对四川省甘孜州西藏飞蝗分布的影响.植物保护,2008,34(1):67-75. http://www.cnki.com.cn/Article/CJFDTOTAL-ZWBH200801020.htm 马世骏,丁岩钦,李典漠.东亚飞蝗中长期数量预测的研究.昆虫学报,1965,14(4):319-338. http://www.cnki.com.cn/Article/CJFDTOTAL-KCXB196504000.htm 姚树然,李春强,王贵生.河北夏蝗大发生与气象条件的关系.昆虫知识,2006,43(4):482-486. http://www.cnki.com.cn/Article/CJFDTOTAL-KCZS200604011.htm 张书敏,任春光,王贵生.东亚飞蝗的长期预测预报技术研究.中国植保导刊,2006,26(7):9-11. 孔海江,陆维松,吕国强,等.干旱前期温度偏高对河南省东亚飞蝗发生的影响.南京气象学院学报,2003,26(4):516-524. http://www.cnki.com.cn/Article/CJFDTOTAL-NJQX200304009.htm 邓自旺,周晓兰,倪绍祥,等.环青海湖地区草地蝗虫发生预报的气候指标研究.植物保护,2005,31(2):29-33. http://www.cnki.com.cn/Article/CJFDTOTAL-ZWBH200502006.htm Muketji M K,Gage s H.A model for estimating hatch and mortality of grasshopper egg population based on soil moisture and heat.Annals of the Entam ological Society of America,1978,71(2):183-190. DOI: 10.1093/aesa/71.2.183
张书敏.河北省东亚飞蝗的发生与治理.北京:中国农业出版社,2006. 项静恬.动态和静态数据处理--时间序列和数据统计分析.北京:气象出版社,1996. 魏淑秋,刘桂莲.中国与世界生物气候相似研究.北京:海洋出版社,1994. 刘静,张宗山,张立荣,等.银川枸杞炭疽病发生的气象指标研究.应用气象学报,2008,19(3):333-340. http://qikan.camscma.cn/jams/ch/reader/view_abstract.aspx?file_no=20080355&flag=1
计量
- 文章访问数: 3446
- HTML全文浏览量: 858
- PDF下载量: 1378