苹果气象服务系统研究——以运城市为例硕士论.doc

成都信息工程学院硕士学位论文苹果气象服务系统研究以运城市为例摘要运城市是山西省的苹果产业基地，是全国苹果优生区。目前苹果面积已达328万亩，年产量达300多万吨，约占全国总产量的18。气象条件是影响苹果产量和品质的重要气候生态因子，研究苹果产业与气候生态环境的关系，对于趋利避害挖掘气候资源潜力，减轻气象灾害危害，加快发展苹果产业，增强运城市苹果市场竞争力，获得更大的经济效益都具有重要的现实意义和深远的战略意义。为了使运城市苹果产业更好的发展，使得苹果的生长发育适合气候变化的趋势，需要深入研究气象影响因子（包含气象条件、气象灾害等）对苹果生长、发育、产量、质量影响及相应对策。运城市的气候条件有适合苹果生长、发育并且高产的优质气候资源，也存在着病虫害和气象灾害等对苹果生长、发育影响的不利因素，例如冬春干旱、花期冻害、冰雹、秋季连阴雨等。本文针对山西省运城市苹果产业实际发展的需求，对影响山西省运城市苹果的上述方面进行了比较全面系统的研究，以期能为提高苹果产业气象服务水平提供更多的科学依据。第一章 绪论1.1 研究背景。1.2 国内外研究现状 1.2.1. 国外相关研究国外对苹果产业化发展非常重视，已进行了大量的苹果气候生态关系研究，美、日、俄等国家在这方面取得了较多研究成果。日本苹果栽培专家根据日本气候生态环境条件的特点，在日本青森县南津轻郡藤崎町进行杂交培育举世闻名的“富士”优等果品2。20世纪60年代开始大力推广富士苹果，该品种被认为是“划时代品种”，是苹果冲向世界的市场竞争的一张王牌。关于气象条件的关系方面，Heinicke对苹果树内的光照分布模式做过研究3；Reinhard Paulesich对形状不同的树冠光照分布状况做过研究4；Michael Glantz研究过苹果不同树型树冠内太阳辐射分布状况5；PaFaHHer研究得出不同地区不同苹果品种需要的正积温9002850，个别品种达到3000左右；HJIpar研究表明在积温30003200左右时，苹果树干最大，产量最高6。A.n.JIocOB研究表明，苹果果实体积与蓄水量之间相关系数达0.96，J.E.Goode研究表明干旱使苹果体积减少50，优果率下降65。在气象灾害方面对苹果的专门研究比较少，但有气象灾害对其他作物的影响，StevenM.Quiring, Timothy N. Papakryiakou分析了加拿大小麦干旱的影响，利用曲线回归模型建立了小麦干旱指数 7 ；Ashok K. Mishra Vijay P. Singh基于概率模型分析干旱时空变化，对干旱进行预测 8 。国外科学家在应用雷达识别冰雹云方面做了大量的工作，Mather把45dBz 高 度0层高度（H0）+3.3km作为地面降雹的判据9。瑞士、法国、意大利等国，将45dBz高度H0+1.4km作为冰雹识别指标10。在果树病虫害与气象条件的关系研究方面，日本走在前列。日本学者森下争产对近100年日本褐飞虱大发生与ElNino的关系进行了研究，得出E1Nino次年为日本褐飞虱大发生年。韩国利用气象卫星收视系统建立了果树病虫害信息可视网络咨询服务系统采用国际互联网，开发果树病虫害监测、预报及管理系统。对于产量预测预报，国外对此研究较早，60年代中期，美国的WGDuncun和RSLoomis与荷兰的Dewit同时提出作物群体模式研究11，现在苏联、美国、荷兰、日本等国家对棉花、小麦、糖用甜菜等作物在农业气象预报方面的研究已取得了明显成果。国外早在20世纪70年代就开始将计算机技术应用到农业生产中，以解决病虫防控和农业生产问题。如：1996年，瑞士开发出谷物预测预报系统EPIPRE，可准确地预报杀虫剂的施用日期，为化学药剂的合理使用提供决策，降低环境污染12。1.2.2. 国内相关研究我国苹果生产从上世纪中期开始起步，80年代迅速发展并从国外引进不少新品种，苹果产量和品质在不断提高。如何合理利用当地气候生态环境优势，充分发挥当地气候资源潜力，建立适宜的栽培管理技术规范，趋利避害，促进苹果产量，尤其是品质提高，已成为各果区关注的热点。深化果业气象服务内容，规范果业气象服务技术体系越来越引起国内外果业管理部门、经销商和种植户的重视。国内对果树与气象条件进行了大量的研究。蒲金涌等13人研究发现，甘肃黄土高原自1971-2005年，年平均气温增加了0.91.8,受此影响，苹果叶芽放、展叶、开花等早期物候期提前了14d左右，叶变色、落叶等后期物候期后延了25d；郭海英等14研究发现随着气候变暖，果树在春季发育的物候期普遍提前，而果树的春季各发育期间隔并未出现缩短的趋势；毛明策等15研究了陕西省长武温度与苹果初花期之间的关系；李美荣16分析了陕西省7个果区气候变化及对苹果开花期的影响,并建立了苹果始花期物候模型；杨小利等17研究陇东黄土高原典型站苹果生长对气候变化的响应。关于苹果花期冻害引起了多位专家学者的关注，已有多位学者对其进行了灾害特征及影响分析、气候适宜性及灾害风险等的研究。关于受冻指标，李健18等以2 、0 、 -2 为不同冻害等级临界气温，采用统计结合试验等方法确定陕西苹果花期冻害主要危害地区为延安和渭北西部果区，主要危害时段为4月中旬；李美荣等以0 、 -2 、 -4 为不同冻害等级临界气温，统计分析得出陕西苹果花期冻害发生频率呈西北向东南减少的特征。董能利19对黄河河道地区果树冻害发生规律及防治对策进行过研究；杨梅宁、王秋萍20分别对苹果花蕾期冻害发生规律及防御对策进行了研究；毛明策等21分别对苹果初花期与温度的关系及气候变暖对苹果物候期的影响进行了研究 ,为防御苹果花期冻害提供气象信息；夏丽花等22利用逐步回归分析建立了短期最低气温预报方程,对果树低温冻害进行短期预报。旱涝气候预测应是气候预测的重点，对旱涝气候规律和预测方法的研究已是当务之急。然而，旱涝灾害特别是降水量的预测是大气科学中最为复杂、最为困难的课题之一。20世纪90年代以来，随着气候系统研究的深入，海气相互作用、陆气相互作用、冰雪大气相互作用等研究领域都取得了令人瞩目的成果。为气候预测提供了新理论和新方法。季劲钧23对陆面物理和生物学过程进行了模拟研究。吴国雄等研制出全球海洋大气陆面系统模式(GOALS)，广泛应用于气候研究。汤懋苍从地气相互作用出发，根据冬季深层地温场与夏季降水场有明显的相关关系，研制出“地气图法”，多年来用于制作汛期旱涝预报，取得了一定的效果。利用场时间序列做区域旱涝气候预测是最诱人的一个应用方面。于是，我们以陆气相互作用为物理基础，结合运城的自然环境和气候特点，建立预测指标，进行预测试验。多年来，国内外许多专家致力于冰雹预报方法的研究。如雷雨大风等。冯晋勤24对福建地区的冰雹天气进行了系统研究，并详尽介绍了冰雹的形成机理和预报方法，为以后开展冰雹研究奠定了基础。河南省把40dBz作为有无冰雹的判据之一，陕北地区用711雷达的45dBz强回波顶高H0+2.3km作为判据来识别冰雹云。这些研究为预测冰雹天气提供了参考。我国学者在20世纪80年代以前,主要是以寻找、观察农作物病虫害发生发展的周边气象条件为主,并以指标法、物候、统计等方法进行了短期预报。进入80年代以后,大尺度的气候背景(如厄尔尼诺事件、海温、副热带高压和大气环流特征量等)与农作物病虫害发生流行的相关关系及其亲合机制研究开始出现,气象条件与害虫生命表、农作物产量损失、防治指标和预测模型的研究得到了加强；80年代末始,农作物病虫害的中、短期预测进入研究应用阶段,长期预测研究进入探索阶段；20世纪90年代以来,由于气候变化的影响,气候变暖对农作物病虫害发生流行影响研究逐步开展。在作物产量预测方面，国内其它省研究状况不尽相同，采取的方法也不同。如：用遥感技术方法预报产量，利用ARMA模型及岭回归即以趋势产量和气象产量分别建模的方法，利用灰色马尔柯夫链动态预测模型模拟趋势产量，应用自然正交方法将趋势产量和气象产量展开，利用交叉谱分析方法，从农业气象角度建立不同作物群体的光合生产模式分析作物生长季(59月)气温的变化对粮食产量高低的影响25-29等等。研究的作物也不尽相同，其中对水稻和小麦的研究较多，包括主要对水稻及冬小麦产量的主要影响因子和影响时期进行研究。利用随机天气模型研究气候变化对我国小麦发育及产量的可能影响，分析小麦和水稻的气象产量特征，根据作物产量形成与农艺性状的关系，对Baier-Robertson的模型作了改进，建立了小麦预报模式，利用化肥施用量、一季稻种植面积和l3月平均温度建立预测模型29-32等等。我国在水稻生长动态模拟研究和遥感技术用于冬小麦估产方面取得了较大进展，但在玉米、大豆等早作物发育期、产量跟踪预报技术方面都还没有重要的突破，更没有形成完整的体系。由于近年农产品的产量预报工作主要集中在小麦和玉米等大田作物，鲜有针对果品等经济作物进行产量预报的具体方法或模式。在我国，农业气象服务系统起步较晚。当前，已建立了一些用于苹果病虫害管理平台，例如，莱阳农学院秦淑莲等33研制的苹果园病虫害测报计算机信息系统；中国农业大学赵朋等34开发的苹果病虫害管理信息系统等。这些系统各适合不同的区域范围， 但目前在山西省尚没有专门的，全面的苹果预警系统在业务服务及生产中使用。1.3 本文主要研究内容根据山西省运城市的气候特点，结合苹果生长发育期的气象观测资料，分析运城市历史苹果生长发育期、产量与气象影响因子的关系，结合国内外优秀的苹果生产经验，综合得出，（1）山西省运城市苹果生长发育期适宜的和不适宜的气候条件指标；（2）对影响运城市苹果产业的花期冻害、干旱、冰雹等的临界指标和对运城市苹果生产的影响进行研究；（3）对运城苹果的病虫害的为害状况、发生规律、盛发期、发生的有利和不利气象条件进行全面分析，归纳主要病虫害发生趋势的预报指标；根据气象观测资料为基础，结合当地实际情况，用非线性回归方法建立预测模式，预测运城市苹果单产，为果农和政府指导提供依据；（4）建立运城市苹果主产县苹果面积产量数据库，运城苹果气象灾害、物候观测数据库，苹果生育期气象指标数据库，苹果气象灾害指标数据库，建立苹果产量预报流程；为运城市苹果生产提供有利的决策支持和服务，对运城苹果产业充分利用气候资源、趋利避害，发展苹果生产起到重要作用。第二章 运城市苹果生长发育期气象指标研究苹果在漫长的形成过程中，已经适应了固定的气候和生态环境，所以对气象条件有特别的要求。气象条件的突变会严重影响苹果的产量和品质。研究苹果的各个生长发育周期的气象条件，对于在生产苹果中合理地利用气候资源、集中优势资源，提高苹果的产量和品质有着非常重要的意义。结合诸多苹果生产中心总结的生产经验，综合列出苹果生长发育期有利和有害的气象条件指标，为苹果生产提供决策服务。2.1 影响苹果生长的基本气象指标2.1.1. 气温气温是影响苹果生育的重要气象指标，它是限制苹果分布的决定性因子，对于苹果的产量的影响至关重要。在苹果年生长周期中，温度与其生长发育有密切关系。当春季气温5时芽开始萌动，15以上时开花，开花期的最适温度为 1718。生长期间(410 月)的平均气温 1218、夏季(68 月)1824最适宜生长。生长季温度 1824时叶片光合能力最强，而呼吸消耗的碳水化合物相对较少；如果高于 24，则消耗增加；高至35时，同化产物则全被消耗。夏季 68 月平均气温 26以上的地区，苹果花芽不易分化，产量很低36。苹果需要一定的低温通过休眠，才能正常生长发育。 苹果通过休眠需要7.2以下低温 9001000 小时，山西省运城市苹果主产区能正常休眠。苹果的花蕾能耐-2.8-3.9的低温，花期为-1.7-2.2，谢花后的幼果不能低于-1.1，成熟的果实能耐短暂的-4-6的低温。果实发育的有效基础温度为10，花后气温低于 5或高于20，均能阻碍细胞分裂。果实发育后期的温度状况影响果实的发育进程和质量。夏秋温度对果实着色、果品质量有明显影响。气温日较差大，夜温较低，可提高果实含糖量，有利于着色。运城市日平均气温大于等于5C的持续天数超过252天，一年内日平均气温大于10C的有效积温一般为38934694C，利于苹果的生育。据分析，运城市苹果生产基地的亩产量和48月份日平均气温大于10C的有效积温有着显著的相关性，相关系数在0.514以上（通过了a=0.05的显著性水平检测）。冬季气温偏低、变化幅度大、寒流经常使树木受冻、抽条或者休眠过早解除，开花不齐，容易受到冻害等；夏季气温偏高或者过低容易使树木营养失调，不易开花结果，病虫害严重，影响产量和质量。通过分析年平均气温、生长期（410月）平均气温和苹果生长温量指数，与我国其他苹果产区对比，我国北部的主要苹果产区年平均气温9.211.9C，生长期（410月）平均气温16.218.7C，温量指数78.495.9C。山西省运城市苹果主要产区年平均气温11.914.2C，生长期（410月）平均气温19.321.7C，大于等于10C积温在3893C以上。完全能够满足苹果生长对温度的要求。2.1.2. 降水量苹果正常的生长发育和结果与环境水分含量有着紧密的关系。水分供给过剩或者太少，都会让果树营养失调。在一般土壤环境下，年降水量在500800mm的地区适宜栽种苹果。据分析，山西省运城市苹果主产区苹果亩产量和48月间降水量有显著的相关关系，相关系数在0.578以上（通过了a=0.05的显著性水平检验）。如果年降水量大于800mm的地方，特别是在高温高湿的环境下，苹果果树生长过于旺盛，苹果质量下降并且容易发生病虫害。降水量如果小于500mm，而且分布不均匀经常出现早春干旱和夏季伏旱，会使树体遭受旱灾。山西省运城市苹果区年降水量大多在500680mm，适宜于苹果的生长发育。2.1.3. 光照苹果是喜光、长日照的树种，需要充足的光照条件。光照对于苹果树的生长结果和果实质量有着显著的影响。如果光照充足、紫外线多，苹果树生长缓和，枝条粗壮，果实含糖量高，着色浓，良品率高。当果实获得全日光照的70%以上时，果实着色全红，获得全日光照的70%40%时，果实部分着色，获得全日光照的40%以下时，果实基本不着色，良品率显著下降。因此，直射到果实上的光照状况是影响果实着色好坏的一个重要气象因素。苹果树适宜的光照条件为年日照时数在1500小时以上。据分析，山西省运城市苹果主产区苹果亩产量与48月日照时数有着显著的相关关系，相关系数大于0.540（通过了a=0.05的显著性水平检验）。山西省运城市苹果区，年日照时数一般1500小时左右，年日照百分率为46%52%，对提高苹果产量和质量非常有利。综上分析，得出苹果影响苹果生长的基本气象指标见表3-1所示。表3-1 影响苹果生长的基本气象指标项目气象指标热量条件1、生长季平均气温为1218；2、日平均气温大于等于5的持续天数大于170天；3、日平均气温大于10的有效积温为25003000。水分条件4、年降水量在500800mm。光照条件5、年日照时数为22002800小时。第三章 运城市苹果气象灾害及病虫害研究近年来，苹果产业的迅速发展受到省政府、运城市政府广泛关注，运城市更是把苹果列为经济发展的支柱产业之一。黄土高原苹果种植区生态气候优势明显，多项气象指标符合苹果优生区指标需求。但是，该区大陆性季风气候特征明显，加之地形复杂，导致气候脆弱，气象灾害比较严重，其中花期冻害、伏旱、连阴雨、冰雹、是影响苹果产量和品质的主要气象灾害。同时由于病虫害对果树也有很深的影响。3.1 花期冻害苹果花期冻害是影响苹果产量、品质的主要气象灾害之一。尤其20世纪80年代以来，随着气候变暖加剧，导致苹果花期提前及遭遇低温冻害的几率和强度增加，使得苹果生产风险进一步加大。选取运城万荣、临猗、平陆、盐湖、永济、夏县、芮城7个气象台站19732013年苹果开花期日极端最低气温和日平均气温，作为设计冻害指数的基本气象资料。同时选取上述7个气象台站20032013年的苹果开花期监测结果，作为计算花期冻害指数的物候资料。苹果花期冻害是指在果树开花期，遇有强降温天气过程，气温达到一定值时，使果树花器遭受冻害的气象灾害。开花期气温降至5或以下，花粉遭受冻害，直接影响授粉、受精和坐果；降至-2.0，花器冻死，无法坐果，造成大量花蕾受冻，遭受冻害的花蕾开花后质量差，坐果率低。因此，为了寻找能够反映花期冻害强度并且物理意义明确的指标来表征花期冻害特点，定义花期冻害指数 其中，Ha为花期冻害指数，为开花期日极端最低气温5的日数（d），该日数越多，花期冻害可能越严重；为开花期平均气温（），该值越小冻害风险可能越大。Ha的物理意义为开花期平均气温为1时，可能发生冻害的低温出现的日数。显然，Ha的大小即反映花期冻害的强弱。Ha越大，冻害越重；反之，则越轻。这里定义的Ha是仅从气象角度反映的花器冻害或花粉冻害的轻重特征。苹果开花期（初花终花）天数是设计花期冻害指数的基本参数。起初将开花期定为10d，得到的冻害指数序列自相关系数较低，无法通过显著性检验。可能是因为开花期天数定得太少，使一些开花时间未能包括进去。后来又以最早的初花期到最晚的终花期时段来确定开花期天数，这样一来，多数果区开花期天数长达2527d。经试验冻害指数序列自相关系数也很低。经反复试验比较，最后以各果区苹果盛花期为中心，时间向前、向后各推10d，将开花期天数定为20d，发现冻害指数序列自相关系数明显提高，一般均能通过显著性检验。各果区苹果开花期天数具体确定结果见表4-1所示。从表1看出，虽然各基地县苹果开花期天数均定为20d，但它们的气候意义是不同的，代表了基地县早春（34月）不同气候时段。不同气候时段统计出的Te5和T显然是不同的。重要参数。一般讲苹果花受冻害的气温指标定为-1.7，但实际上当日极端最低气温5时花粉即遭受冻害。花粉遭受冻害后失去活力不能发芽或冻死，直接造成无法授粉、受精，必然使坐果率下降，给产量和品质造成损失。苹果花器遭受冻害，症状明显，如果把花器遭受冻害视为“显性”冻害，则花粉遭受冻害可视为“隐性”冻害，它受冻症状不明显，常不被人们所察觉，从而得不到重视。所以常规的花期冻害的概念应当予以扩展，花粉遭受冻害也应属花期冻害的一部分。起初将开花期极端最低气温Te指标定为0，后又定为2、4，再用低温出现日数、计算Ha，结果发现Ha序列相当一部分、甚至1/3以上指数值为0或数值很小，使指数分级无法进行。后来将Te指标调为5，统计发现Ha序列数值大为改观，指数分级也能正常进行。另外，多数苹果基地县气象站在平川，而苹果园多在山坡，当气象站Te5时，果园气温往往2，况且树冠气温常低于百叶箱气温，极易出现辐射降温型晚霜冻，使苹果遭受花期冻害。根据运城市果区气候特点和近年来花期冻害发生分布情况，确定花期冻害指数的频率分布，并将频率适当向偏重等级方向倾斜，使Ha整体频率分布呈略偏态分布。这种做法对于减少冻害漏报有利。将7个基地县的花期冻害指数分成重、偏重、中等、偏轻和轻5个等级，详见表4-2所示。利用气象要素自身年变化具有准周期性和自记忆性的自相关关系 ,建立自回归预测模式表4-3所示37-39。为了检验自回归预测模式的预报效果，分别用7个基地县19732010年、19732011年和19732012年的花期冻害指数序列，算得|RK|max为0.3650.618，均能通过=0.05的显著性检验。利用上述预测模式，对2010、2011和2012年的花期冻害指数作独立样本预测试验，共得到21个预测结果。经与实际花期冻害指数比较，其中预测的花期冻害等级与实际的花期冻害等级一致的有8个，相差一级的有6个，相差两级的有4个，相差三级的有3个。若将预测的结果按三级评定：一致的和相差一级的视为预测准确和基本准确，达66.7%，相差两级的视为预测不准确，为19%，相差三级的视为预测错误，为14.3%，该预测模式能较客观地反映出运城苹果生产基地县的花期冻害强度和趋势，具有一定的预测能力。3.2 伏旱伏旱对果树的影响是广泛而深刻的，表现在生长发育各个阶段，如萌发、营养生长，直达具体的生理代谢过程，如光合作用、呼吸作用，水和营养元素的吸收、运输，一些酶的活力和果树体内某些有机物的消长。选取运城万荣、临猗、平陆、盐湖、永济、夏县、芮城7个站19732013年每年7月中旬8月上旬的降水量和无降水日数，作为设计伏旱指数的基本气象资料。伏旱是一种伏期长期干燥少雨稳定的气候现象,表现为缺乏足够的降水,反映的是气候水热平衡特性。常规多用降水量距平来讨论伏旱，但是由于降水距平百分率对于平均值的依赖较大,对于降水空间分布极不均匀的运城市果业基地不宜使用统一的降水距平百分率标准 ,应当确定一个能够反映伏期降水特征、干旱状况和伏旱强度，并且物理意义明确的指标来表征伏旱特点。我们设计了用如下表达式来定义伏旱指数:其中Fh为伏旱指数，NR=0为伏期无降水日数，该日数越大，伏旱强度可能越强；R7中-8上为伏期降水量，该值越小伏旱可能越严重。该表达式的物理意义为伏期每毫米降水需要维持的无降水（伏旱）日数。Fh越大，伏旱越严重；反之则越轻。根据运城市干旱气候背景特点，将7个基地县的伏旱指数分成强、偏强、中等、偏弱和弱5个等级，详见表4-4、表4-5表所示。表4-4 运城苹果基地县伏旱指数分级伏旱等级伏旱强度万荣临猗平陆盐湖永济1弱Fh0.10Fh0.16Fh0.08Fh0.12Fh0.122偏弱0.10Fh0.160.16Fh0.250.08Fh0.110.12Fh0.150.12Fh0.163中等0.16Fh0.300.25Fh0.550.11Fh0.250.15Fh0.300.16Fh0.364偏强0.30Fh0.630.55Fh1.500.25Fh0.600.30Fh0.520.36Fh0.755强Fh0.63Fh1.50Fh0.60Fh0.52Fh0.75伏旱等级伏旱强度夏县芮城1弱Fh0.12Fh0.122偏弱0.12Fh0.180.12Fh0.183中等0.18Fh0.300.18Fh0.304偏强0.30Fh0.600.30Fh0.905强Fh0.60Fh0.90利用气象要素自身年变化具有准周期性和自记忆性的自相关关系，建立自回归预测模式表4-5所示。为了检验自回归预测模式的预报效果，分别用7个基地县19732010年、19732011年和19732012年的伏旱指数序列，算得21个|RK|max为0.3060.620，均能通过=0.05的显著性检验。利用上述的预测模式，对2010、2011和2012年的伏旱指数作独立样本预测试验，共得到21个预测结果，经与实际伏旱指数比较，其中预测的伏旱等级与实际的伏旱等级一致的有6个，相差一级的有9个，相差两级的有5个，相差三级的有1个。若将预测的结果按三级评定：一致的和相差一级的视为预测准确和基本准确，达71.4%，相差两级的视为预测不准确，为23.8%，相差三级的视为预测错误，为4.8%，说明该预测模式能较客观地反映出运城市苹果生产基地县的伏旱强度和趋势，具有一定的预测能力。3.3 连阴雨连阴雨过程常常与低温过程相伴，造成降水雨量大，持续时间长、气温低、光照不足，致使部分果园积水、内涝和湿度过大，极易诱发裂果、烂果、落果和果树病害流行，对即将成熟的果品产量和品质产生影响。选取运城万荣、临猗、平陆、盐湖、永济、夏县、芮城7个站19732013年每年9月中旬10月上旬的3d降水日数和无降水日数，作为设计连阴雨指数的基本气象资料。运城市对连阴雨天气的定义为测站连续4天及以上，日降水量大于或等于0.1毫米，且测站过程降水大于20毫米的降水天气。而在近几年的苹果气象服务中，笔者发现在苹果摘袋、着色和采收时期，尤其是着色和采收期，持续3天的连阴雨，以及连阴雨或连阴天所带来的低温天气即可对苹果着色、果面光滑度产生影响，且影响程度随连阴雨天气长度的增加而增加，但与过程降水量的多寡关系不大。所以为了能够反映连阴雨天气长度对苹果摘袋、着色和采收的影响强度，我们设计如下表达式来定义连阴雨指数:其中为连阴雨指数；为9月中旬10月上旬降水（R0.1）连续3天及以上的日数，该日数越多，连阴雨危害越重；为9月中旬10月上旬无降水日数，该日数越多，连阴雨危害越轻。将7个基地县的连阴雨指数分成重、偏重、中等、偏轻和轻5个等级，详见表4-6所示。表4-6 运城市苹果基地县连阴雨指数分级连阴雨等级连阴雨强度万荣临猗平陆盐湖1弱Lu0.10Lu0.15Lu0.10Lu0.152偏弱0.10Lu0.310.15Lu0.350.10Lu0.350.15Lu0.303中等0.31Lu0.50.35Lu0.600.35Lu0.600.30Lu0.704偏强0.5Lu1.00.60Lu1.200.60Lu1.00.70Lu1.405强Lu1.0Lu1.20Lu1.0Lu1.40连阴雨等级连阴雨强度永济夏县芮城1弱Lu0.18Lu0.10Lu0.122偏弱0.18Lu0.500.10Lu0.300.12Lu0.283中等0.50Lu0.800.30Lu0.580.28Lu0.654偏强0.8Lu1.300.58Lu1.200.65Lu1.155强Lu1.30Lu1.20Lu1.15根据7个基地县19732013年连阴雨指数序列，用上述统计方法建立各基地县连阴雨指数K阶自回归模式见表4-7所示37-39。表4-7 各基地县连阴雨指数K阶自回归模式（=0.05）基地县NKK阶自回归预测模式|RK| maxf|RK|2万荣40300.514170.389临猗40270.345210.352平陆40300.538180.368盐湖40260.394220.354闻喜40300.467180.375新绛40300.459170.388芮城40290.514190.367为了检验自回归预测模式的预报效果，对2010、2011和2012年的连阴雨指数作独立样本预测试验。若将预测的结果按三级评定：一致的和相差一级的视为预测准确和基本准确，达81.5%，相差两级的视为预测不准确，为16.5%，相差三级的视为预测错误，为2%。该预测模式能较客观地反映出运城市苹果生产基地县的连阴雨强度和趋势，具有一定的预测能力。3.4 冰雹灾害 冰雹这种灾害性天气通常也造成果树枝叶、花、果实机械损伤，从而引起各种生理障碍、诱发病虫害等，使果实商品率下降甚至绝收。3.4.1. 资料收集与整理收集近些年来运城冰雹天气过程。同时，选取近些年运城市测站有风暴强中心大于50dBz、风暴垂直积分液态含水量（VIL）大于20 kg/m2、回波顶高（ET）大于8km以上的强风暴，运城市的冰雹过程用于预报方程的建立，和方程回代检验。将冰雹过程的雷达资料进行整理回放，结合冰雹发展的特点，利用历史冰雹过程的天气雷达资料建立天气雷达判别指标；同时，利用常规探空资料统计冰雹过程的0层、-20层、垂直风切变、-20与0层的高度差等参数。利用二次Fisher判别分析方法和模糊逻辑方法开展运城市冰雹自动识别方法的对比分析。利用天气雷达风暴的单体识别和跟踪算法（SCIT）对风暴自动识别的结果，对出现冰雹、强风暴的雷达产品从风暴生成到消亡进行逐时次统计，统计参数分别为：风暴所对应的时间、方位、距离、风暴垂直积分液态含水量（VIL）、风暴最强回波强度（DBZM）、最强回波所对应高度（HT）、风暴顶高（TOP）、风暴移向、移速（FCST MVMT）及回波顶高（ET）。在建立判别方程时，冰雹的雷达参数以出现冰雹时刻为准，无具体时间的则以报道时间前一小时内最强回波出现的时间为准。同时，利用常规探空资料统计每次强天气过程的对流有位能（CAPE）、0层、-20层、垂直风切变、-20与0层的高度差等参数40。对风暴的单体识别和跟踪算法（SCIT）进行检验，我们利用软件对运城出现冰雹天气对应的风暴发展过程中风暴的位置预报进行12、30、42、60分钟的误差检验。检验中根据风暴提供的未来移向、移速计算风暴未来12、30、42、60分钟的方位和距离，并与风暴的实况做比较，计算方位和距离误差。3.4.2. 冰雹判别指标为对冰雹天气进行定性判别，需选取与产生冰雹相关性好的参数，首先选取风暴最强回波强度作为反映风暴强弱的重要指标；风暴垂直积分液态含水量（VIL）表示将反射率因子数据转换成等价的液态水值，是用来判断对流风暴强度的一个十分有用的参量，但它的阈值随季节变化，Amburn研究表明，如果VIL密度（VIL/ET）超过4g/m3，则风暴肯定会产生直径超过2cm的大冰雹，因此，选取VIL密度指标用来判断冰雹天气；由于冰雹增长都在0层高度以上，0层上部的VIL反映的是0以上降水粒子分布的情况，所以，零度层高度以上VIL也能反映冰雹生长过程，因此我们选取VIL密度*（回波顶高-0层高度）来表示冰雹增长的情况；同时，判断大冰雹最有效的方法是检查较高的反射度因子能否扩展到0，特别是-20等温线高度以上，因此，我们选取最强回波高度-0层高度来表示冰雹发展的强度，风暴顶高（-20层高度）表示冰雹云发展的高度。为验证以上参数选取是否适用，首先，先简单分析产生冰雹天气的风暴以及强风暴的风暴垂直积分液态含水量（VIL）、最强回波强度（DBZM）及回波顶高（ET）特征，经统计（见表4-8），产生冰雹天气的风暴VIL值以及风暴的最强回波强度值较大，而冰雹对应风暴的回波顶高差别并不大，但由于VIL用来判断对流风暴强度的的阈值随季节变化，因此，仅以VIL、最强回波强度及回波顶高进行冰雹天气判断是不够的。表4-8 强天气雷达产品参数统计特征值 雷达参数天气类型风暴垂直积分液态含水量（VIL）风暴最强回波强度（DBZM）回波顶高（ET）冰雹4558641112强风暴243756371027再对选取的5个判别指标进行统计分析。以下分别用DBZM、VILM、H0、H20、VILE表示风暴最强回波强度、VIL密度、最强回波高度-0层高度、风暴顶高（-20层高度）和VIL密度*（回波顶高-0层高度）5个判别指标，由表4-9可以看出，冰雹对应风暴及强风暴在这5个判别指标上差别较大，最强回波高度在0层以上，风暴顶高与-20层高度的厚度差(H20)较大，说明对应冰雹天气的风暴发展强度强、高度高。为进一步验证这5个判别指标的可靠性，分别计算这5个指标与冰雹天气的相关系数分别为：0.69、0.77、0.50、0.57和0.62，当相关系数0.4即为显著性相关，0.7时则为高度相关，所以，这5个指标与冰雹天气的相关性都较高，以VIL密度(VILM)0.77相关系数为最大。表4-9冰雹天气雷达产品判别指标统计特征值 雷达参数天气类型DBZMVILMH0H20VILE冰雹644.250.842.9726.72强风暴56.373.22-0.57-0.1913.28通过对强风暴，冰雹过程中对应风暴计算5个识别参数的分布概率密度图，计算过程中强风暴合并为强雷暴，分布情况见图4-1。图4-1 冰雹5个识别参数的分布概率密度3.4.3. 预报方程的建立（1）判别分析法判别分析法是定性分类判别的一种方法，它是根据判别对象不同类别，选择适当的判别因子，利用在不同类别的样本内，寻找判别因子与判别对象的关系，建立判别因子与判别对象不同类别的预报方程，选择适当的判别规则，判别某次多个判别因子观测样本所属的天气类型。常用的判别分析法主要有：距离判别、Fisher判别、Bayes判别、逐步判别法。气象上主要采用Fisher判别、Bayes判别、逐步判别法。一般来说，Bayes判别要求数据分布是多元正态分布，Fisher判别对数据分布没有特殊的要求，逐步判别是在对判别因子进一步筛选的基础上建立判别方程。由于判别分析方法的一个特点即计算量巨大，在统计过程中借用最常用的SPSS统计分析软件进行预报方程的建立。为了便于比较，在预报方程建立前，需对选取的判别指标用Xzi=（XiX）/S （X为平均值、S为标准差）公式进行无量纲化处理，将无量纲化处理后的冰雹天气对应5个判别指标的样本输入SPSS软件，软件将自动输出两组预报方程。在预报中先对风暴进行筛选，若风暴强中心大于50dBz、VIL大于20 kg/m2、回波顶高大于8km，则将相应的预报因子代入冰雹预报方程，比较Y1、Y2大小，若Y1Y2，则预报为冰雹天气；否则预报为无冰雹天气。冰雹预报方程如下：Y2 = VILE * 2.62 + H20 * 1.131+H0 *0.568+DBZM *1.128+VILM *(-1.327)-3.283（2）模糊逻辑法模糊逻辑法是L.A.Zadeh教授于1965年提出的，主要包括三个步骤：模糊化过程、模糊逻辑推理以及精确化计算。所谓的模糊化即将识别参数以隶属函数的方式转化为模糊基，我们选取梯形函数作为隶属函数的基本形式进行模糊化，函数表达式如下：在构造冰雹的隶属函数时，使用参考门限值来自于图3中识别参数的概率密度，对于DBZM， =60，=65；对于VILM，=2，=3.5；对于H0,=0.5，=2；对于H20， =0，=2；对于VILE，=20，=30。下图4-2给出了基于梯形函数构建的5个识别参数隶属函数对应的冰雹风暴的模糊基。模糊逻辑推理相当于求冰雹或强雷暴的条件概率，即：其中，代表风暴类型，即冰雹，代表识别参数，n代表识别参数的个数，为第j个识别参数对第i类型风暴的贡献，是第j个识别参数的权重，我们这里根据识别参数与冰雹的相关系数，取=0.25，=0.3，=0.05，=0.15，=0.25。当计算出条件概率大于0.5时，即判断为冰雹，小于0.5判断为强雷暴。图4-2 基于T函数的冰雹识别参数模糊基设置同理，在构造强风暴的隶属函数时，使用参考门限值来自于图4-1中识别参数的概率密度，对于VIL下降值，=0，=10；对于强回波中心高度下降值，=0，=2；对于H20，=0，=1。图4-3给出了基于梯形函数构建的3个识别参数隶属函数对应的雷雨大风的模糊基。取=0.5，=0.25，=0.25。当计算出条件概率小于0.5判断为强风暴。图4-3 基于T函数的雷雨大风识别参数模糊基设置关键技术：利用Fisher二次判别方法、模糊逻辑法分别建立冰雹天气预报方程。表4-10强风暴雷达产品判别指标统计特征值 雷达参数天气类型强回波中心下降高度VIL下降值H20强风暴-0.04-1.03-0.19通过对强风暴对应风暴计算3个识别参数的分布概率密度图，分布情况见图4-4。对比图4-3和图4-4，冰雹识别参数的分布概率比较集中，说明了冰雹识别参数具有代表性。图4-4 强风暴别参数的分布概率密度关键技术：根据冰雹天气特点及雷达回波特征，建立基于天气雷达产品的冰雹5个判别指标：风暴最强回波强度、VIL密度、最强回波高度-0层高度、风暴顶高-（-20层高度）和VIL密度*（回波顶高-0层高度）5个判别指标。建立基于天气雷达产品的雷雨大风3个判别指标：强回波中心高度、VIL下降值和风暴顶高-（-20层高度）。3.5 苹果病虫害运城苹果主要病虫害有早期落叶病、斑点落叶病、叶螨等。对病虫害的发生规律，盛发期、有利和不利发生的气象条件做了全面的分析，归纳出用来做病虫害的发生规律和气象条件指标。(一) 早期落叶病早期落叶病是能够引起苹果早期落叶的一类病害的总称，包括生理性病害和病理性病害。这里主要介绍病理性引起的落叶病。引起早期落叶病的病害主要有以下几种：生理性病害（由不良的环境条件等引起）、传染性病害（由各种病原菌侵染引起）： 褐斑病 斑点落叶病 圆斑病 灰斑病 轮班病。在运城市发生和危害比较重的主要是褐斑病和斑点落叶病。病虫害发生的气象条件及指标如下表4-11所示：表4-11 苹果早期落叶病病虫害名称:苹果早期落叶病危害部位主要危害叶片,造成大量早期落叶,削弱树势，影响产量和花芽分化发生规律续表4-111褐斑病：5月下旬开始发病，78月盛发，迅雨来的早，雨量多，发病早而重。老树，弱树发病重。2斑点落叶病：5月中旬开始发病，6月中旬进入急增期，6月下旬7月上中旬进入盛发期。春季雨水多，空气湿度在70%以上发病重。危害盛期褐斑病：78月；斑点落叶病：6月下旬7月上中旬；8月。有利气象条件不利气象条件1早春（3月）气温普遍偏高，有利病菌越冬；2早春（3月）降水较多，有利病菌越冬；3迅雨较常年偏早，对发病不利；4 67月降水量偏多，有利该病发生；1早春（3月）气温偏低不利病菌越冬；2早春（3月）降水偏多，对病菌越冬不利；3迅雨较常年偏迟，对发病不利；4 迅雨（67月）降水少，不利该病发生；偏重发生气象条件指标偏轻发生气象条件指标（1） 冬前调查病叶率19.1%（2） 3月平均气温距平0.2；（3） 3月降水量距平-10；（4） 夏季汛雨早，67月降水量距平10%（1） 冬前调查病叶率30%（2） 122平均气温距平0.2（3） 3月上旬5月中旬平均气温距平0.2（4） 3月上旬5月中旬降水距平10%（1） 越冬基数病株率3%（2） 122月平均气温-0.2（3） 3月上旬5月中旬平均气温-0.2（4） 3月上旬5月中旬降水量距平-10%3.6 本章小结(1) 对运城市苹果主要的气象灾害如花期冻害、伏旱、连阴雨、冰雹等建立了其相应的判别指标，根据判别指标进行了灾害等级划分，并建立了预测模型。(2) 对运城苹果主要病虫害为害状况、发生规律，盛发期、有利和不利发生的气象条件做了全面的分析，归纳出用来做病虫害趋势预测技术和气象条件指标，结果编成简表。第四章 苹果产量预报目前，中国主要苹果产区还没有针对苹果产量预测方法的系统性研究成果，仅在2000 年初个别研究者对苹果产量定性预测思路做过少量研究，在2000 年后鲜有研究者在苹果产量预测技术方法上有所突破，无法满足苹果产业发展对气象服务的需求。运城市是苹果生产大市，果业气象服务需求量大，且质量要求，尤其是对制定营销方案而需的苹果产量预测产品的需求更加突出。因此，有必要围绕苹果产量预测的决策需要，并着眼于对苹果产业长期发展的新要求，开展苹果产量定量预测方法研究，逐步建立苹果产量预测业务体系，弥补此业务的空白，以期为运城市合理制定每年苹果营销策略、规划苹果产业发展，提供科学有力的决策依据，为其他省份开展苹果产量预测及其他经济林果产量预测业务及科研提供参考依据。4.1 预测数据来源运城市目前共有7个苹果生产基地县，本研究拟从其中选择万荣县、临猗县、平陆县3个在地理位置、生态环境、气候特点和苹果生产水平具有代表性的基地县，作为苹果产量预测方法研究对象。根据苹果树栽培的气候生态条件，影响苹果产量与品质的关键气象条件，以及运城市苹果主要物候期及成熟采摘时间，确定苹果主要生长期及果实膨大期气候资源为预测因子，包括48 月10积温、48月日照时数、48月降水量和68月降水量。受运城市果业局苹果产量数据限制，各县单产数据为20002014 年，气候资源数据亦配合产量数据，数据序列为20002014年。4.