《数值天气预报》课件：第一讲数值预报介绍.pptx

1、李江南李江南2016年年2月月 15989069278QQ:894897370非线性多尺度多圈层相变（跃变）多种强迫和反馈动力、物理和化学综合大气演变过程的复杂性 要掌握在地球系统多圈层相互作用影响下大气的动力、物理和化学的演变规律，依据已知的信息，预知大气未来的演变过程和状态。 希望能准确告诉人们未来的大气状态。希望能准确告诉人们未来的大气状态。但气象问题的复杂性使得人们无法通过自身大脑的直接思维去准确掌握大气的演变趋势，也没有一个科学家能以理论方式直接求解数学物理方程得到大气未来的状态，也难于直接通过物理实验来客观地认识大气的演变规律。l当前的天气或气候信息完备的综合观测系统l天气或气候的

2、演变规律从资料得到新认识l外力和强迫的变化地形和边界强迫；太阳常数；引力l从已知预报未来的手段完备的数值模式；承载数值模式计算、显示和通讯平台驾驭观测、资料分析和数值模式发展和应用的人才队伍 一是天气预报，二是气候预报一是天气预报，二是气候预报 经常不断变化着的大气状态，既是一定时间和空间内的大气状态，也是大气状态在一定时间间隔内的连续变化。 指发生在大气中发生的各种自然现象，即某瞬时内大气中各种气象要素（如风、云、雾、雨、雪、霜、雷、雹等）空间分布的综合表现。 一定地区的天气现象随时间的变化过程。 长时间内气象要素和天气现象的平均或统计状态，时间尺度为月、季、年、数年到数百年以上。气候以冷、

3、暖、干、湿这些特征来衡量，通常由某一时期的平均值和离差值表征。 天气、气候现象都是地球大气运动的结果。天气现象更侧重于描述短期内大气具体的演变过程而气候现象主要是指长期、代表大气各种指标平均值的变化过程。 不同对象的数值预报所使用的技术方案与预报产品都有很大区别，例如一般的数值天气预报关注并预报具体的天气过程与气象要素的演变，而短期气候模式预测则关注、预报月与季节尺度的冷暖、旱涝趋势而不是具体的天气过程。 天气天气和和气候气候模式研究的区别模式研究的区别研究对象研究对象：单一大气单一大气 耦合系统耦合系统关注目标关注目标：瞬时结果瞬时结果 平均状态平均状态关键变量关键变量：降水降水 温度温度关

4、注内容关注内容：演变过程演变过程 能量收支能量收支 精度精度 稳定和守恒稳定和守恒水汽过程水汽过程：凝结降水凝结降水 云辐射强迫云辐射强迫 云辐射反馈云辐射反馈决定因素决定因素：初值初值 强迫和反馈强迫和反馈同化同化 物理和化学规律物理和化学规律 大气并没有向我们展现出类似潮汐变化那种类型的周期性，因此我们无法采用像潮汐预报那样的方法来预报天气。我们没有发现大气在某一特定时刻的状态与另一时刻的状态之间存在一系列简单的因果关系。运用数值方法的天气预报运用数值方法的天气预报， 1951, Charney（1917-1981）数值天气预报的巨匠）数值天气预报的巨匠之一之一151820年，世界上诞生了

5、第一张天气图。1855年，法国人莱伐尔用天气图追索克里米亚战争时出现的风暴。1860年，英国人弗兹罗用电话收集天气报告，绘成天气图。第一次世界大战时,地面天气图已成为天气预报的工具。二战以后，随着高空探测的增多和世界高空探测站的逐渐加密，天气分析从二维扩展为三维。地面观测初始时期地面观测初始时期极锋学说应用时期极锋学说应用时期探空发展时期探空发展时期斜压理论发展时期斜压理论发展时期卫星观测初始时期卫星观测初始时期气象雷达推广应用时期气象雷达推广应用时期动力数值天气预报应用时期动力数值天气预报应用时期动力统计预报应用时期动力统计预报应用时期人工智能系统应用时期人工智能系统应用时期 现代气象学的研

6、究（探测基础上）已把大量的物理、化学定律引入到对天气系统各种演化过程的描述上来。因此，现在我们可以用一套微分方程组来表示大气的运动和变化（数学模型）（质量守恒、能量守恒、动量守恒和水汽守恒）。 20 世纪,在世界(国际) 气象组织的组织推动下,世界各国开展广泛合作,逐步建立了遍布全球的气象观测站网,人们根据观测资料分析绘制过去时刻天气变化的天气图,并以此为基础发展了天气学(定性) 预报方法以及结合数理统计学的统计天气预报方法。但是,由于这些方法并未反映大气内部的运动变化规律由于这些方法并未反映大气内部的运动变化规律,所以其预报水平的提高受到很大的限制。因此,要进一步提高天气预报准确率,必须从大

7、气内部的物理规律(大气动力学、热力学、质量守恒、水相变化规律等) 出发,建立数学物理模型,用数学物理的方法,并借助巨型计算机技术,预测未来天气的变化,这就是我们常说的“数值天气预报”。21理解理解“数值数值”二字二字将支配大气运动的偏微分方程，以观测值为初始条件，利用计算机的数值解法(时间积分)求得的理论预报称为数值预报。因为控制方程是非线性的，得不到解析解，只能用数值求解的方法，将实况观测值作为初始条件(初值)，利用电子计算机求解(时间积分)，所以加上“数值”两字。Meteorologists need relationships (equations) between variables

8、you want to know (e.g., v, T, p, z, q) and the forcing mechanisms that cause changes in these variables.Example:What Is a Numerical Weather Prediction Model? Or Example of a Prognostic Equation:数数字化，定字化，定量，客量，客观观 预报是建立在一组对大气运动规律了解了的数学物理预报是建立在一组对大气运动规律了解了的数学物理方程组上。方程组上。代表了数字化、定量客观化的预报方向。代表了数字化、定量客观化的

9、预报方向。可以预报天气的极值变化。可以预报天气的极值变化。初值初值：从站点值到格点值的插值：从站点值到格点值的插值四维资料同化是尽可能利用不同时刻的多种观测资料得到一个协调的适合模式特点的能较准确代表大气或气候系统初始状态的模式初试场。模式动力框架模式动力框架 从世界上发达国家/ 地区及中国在用业务模式动力框架方面的主要技术特征来看,已由多尺度多模式动力框架的阶段进入了多尺度一体化模式动力框架的阶段。即过去一个业务数值预报中心根据不同的预报对象(如全球中期天气预报、区域短期天气预报、中尺度天气预报、台风/ 飓风预报) ,分别建立了基于不同模式动力框架的数值天气预报业务模式。而今,各国业务数值预

10、报中心已逐步建立了可用于不同尺度预报对象的、基于多尺度一体化模式动力框架的数值天气预报业务模式。(1) 20 世纪90 年代中期以来,各国的全球和区域/ 中尺度模式水平和垂直分辨率都有明显提高,全球模式水平分辨率已由100120 km 普遍提高到目前的60 km 左右,模式垂直层数一般都多于30层;传统的区域模式已被中尺度模式所取代,中尺度业务模式水平分辨率目前普遍在10 km 左右。(2) 对于空间离散,有限区域模式多采用欧拉差分格式,全球模式则多采用拉格朗日差分格式;对于时间离散,半隐式和将快慢波分离(分裂- 显式)的时间分裂是较常用的时间积分方案。将半隐式和半拉格朗日方法结合一起的方案已

11、在业务数值预报模式中得到了最广泛的使用。(3) 模式动力框架的改进已不再是单一进行,而是与模式物理过程的改进同步进行,必须考虑模式动力过程与模式物理过程之间的协调性,以及模式的整体性。( 4) 主要发达国家和中国都正在致力于研发各自的新业务模式非静力多尺度一体化模式或非静力中尺度模式。模式物理过程参数化模式物理过程参数化 模式物理过程参数化方案的优化和改进对数值模式预报准确率的提高起着关键的作用。近年来, 人们对大气物理机制的认识,有助于模式物理过程参数化方案的改进,使数值预报模式更进一步完善。(1) 在全球中期数值预报模式中,发展和改进格点尺度云(凝结物) 预报方案已同改进积云对流参数化方案

12、一样受到重视。 过去,因全球数值预报模式分辨率比较低,模式中的大尺度(即格点尺度) 凝结过程的处理比较简单, 人们更多的注意力放在对积云对流参数化方案的发展和改进上。由于在此过程中未考虑云的生成过程,使得模式大气中缺少了水物质(云水、雨水、冰、雪等) 的拖曳作用,当模式分辨率逐渐提高时垂直速度会出现虚假增长,从而导致模式降水量的虚假增加。这一问题在模式的水平格距减小到100 km 以下后逐渐变得较为突出。(2) 区域/ 中尺度模式中云物理过程显式预报方案的细化和逼真描述成为改进模式降水预报的重点。 随着业务区域/ 中尺度模式分辨率的提高,更多的注意力被集中在细化和逼真地描述云物理过程的显式预报

13、方案上,期望改善降水预报效果。(3) 陆面过程的描述及陆气耦合的改进在气候数值模拟研究中的重要作用不容置疑。近些年来的研究又显示,陆面过程方案及其细化描述在中短期高分辨率数值天气预报模式中的作用也是非常重要的。模式程序软件模式程序软件 随着模式分辨率的提高,计算量呈几何级数式地增长,对计算机资源的需求也随之急速增加。大规模并行计算机是能满足这一需求。气象数值预报模式运行要达到最高的计算效率,模式程序的并行化是不可避免的。 另一方面,数值预报模式性能的完善化,使得模式程序规模进一步扩大,加之全球/ 有限区、天气气候一体化模式的提出,模式程序更进一步复杂化,使得模式程序的研制、运行、维护、发展变得

14、更加困难。因此，模式大型程序的设计应按现代软件的规范来进行：(1) 标准化。模式程序每一单元(主程序、子程序、函数、循环语句、变量名称、注释行等) 的编写必须按给定的一体化编程标准严格进行。现在随着模式程序规模越来越大,参与模式编程的人越来越多,模式程序编写的标准化越来越重要。(2) 模块化。按数值模式的功能和算法(如一个积云对流参数化方案、一个辐射方案、一个时间差分方案、一个插值计算方案等) ,将标准化的模式程序单元组合成一个个模块,以便人们可以“插- 拔”式选择不同的模块装配成一套完整数值预报模式,用于不同的预报目的。(3) 并行化。随着模式性能的完善、分辨率的提高,模式的计算量也呈几何级

15、数增加,需要更高性能的巨型计算机才能实现数值模式的大规模科学计算。高性能计算机的有效使用必须采用并行化的模式程序。初初始始状状态态未未来来状状态态数值模式对演变规律的对演变规律的综合理论认识综合理论认识外力和强迫外力和强迫模式结构模式结构数值预报流程图数值预报流程图气象观测资料资料同化数值预报模式模式后处理 数值预报释用产品分发 中央计算机群系统通信及计算机传输网络通信及计算机传输网络资料同化资料同化 逐时刻逐时刻“推算推算”未来天气演变未来天气演变数值天气预报原理数值天气预报原理空空 间间 格格 点点 逐一格点求解逐一格点求解一组数学方程一组数学方程数数 值值 模模 式式 天天 气气 演演

16、变变 图图 客观分析和资料同化系统主要解决从观测资料到模式可用资料的字转换。初始化过程去除资料中的干扰。数值模式一组数理方程组数理方程组(干模式）物理过程参数化描述大气中辐射、行星边界层、积云对流、海气相互作用、微量气体等（湿模式）。数值产品的后处理包括分析诊断再加工各类产品、图形化 显示。数值预报分类方法数值预报分类方法根据预报对象的时间尺度 分： 短期数值天气预报（13天） 中期数值天气预报（410天） 气候模式预测等（10天以上）根据预报的空间范围与尺度 全球数值预报（T213L31） 区域数值预报（HLAFS） 中尺度数值预报（WRF）新一代多尺度通用数值预报体系（GRAPES）数值预

17、报产品的特点数值预报产品的特点种类繁多 1.基本气象要素 2.基于基本气象要素通过计算间接获得的反映大气热力和动力瞬时状态的物理量。 3.模式直接输出的近地面层天气要素 4.专项预报产品（台风路径、降水等）空间、时间分辨率高时空分布的连续性好 数值预报的分析产品在空间分布及时间演变的连续性上比根据实际观测资料人工分析要好，特别是在实况资料稀少的海洋和高原地区。预报误差特征极其复杂产生误差的原因产生误差的原因分析误差 目前，观测系统并不完全按照天气预报的要求建立的，而且观测资料包含各种不同类型、不同分布密度、不同观测频率和观测精度。基于这种不完善的观测系统基础，所得到的资料同化分析场与真实大气之

18、间必然存在差异，这种来自分析场上的误差导致了模式计算上的误差。模式误差 模式的水平和垂直分辨率不够精细，物理过程参数化不够完善，难免有这种或那种的假定或简化，很难完全描写真实大气特征而造成误差。 1904年,挪威学者V.Bjerknes在世界上首次对数值天气预报理论作了非常明确的表述,认为大气的未来状态原则上完全由大气的初始状态、已知的边界条件和大气运动方程、质量守恒方程、状态方程、热力学方程所共同决定。换句话说,在给定大气初始状态和边界条件下,通过求解描述大气运动变化规律的数学物理方程组,可以把未来的天气“较精确地计算出来”。当时困难： 1）方程组非线性，无法解出解析解。 2）没有好的观测资

19、料，无法确定大气的“真实”初值。 英国科学家Richardson首次提出直接用“数值积分”这些方程的方法求解这一问题:取未经简化的完全原始方程,水平格距200 km ,垂直4 层(层顶200 hPa) ,中心位于德国,1910 年5 月20 日07 世界时的观测作初值。他借助一把10 英寸的滑动式计算尺,制作出了世界上第一张6 h 地面气压数值预报图(时间积分为1910 年5 月20 日0410 世界时) 。可是,这张地面气压预报图“预报”的6 h 气压变化为146 hPa ,实际观测气压几乎没有多大变化,从精度上看该预报毫无“参考”价值,而且其计算时间花了将近一个月,从时效上也已毫无“预报”

20、意义了(他估计用这样的计算工具,要从时效上做出有“预报”意义的天气预报来,必须64 000 人同时进行模式计算才行) 。可以说,Richardson 的首次数值天气预报是失败的。L. F. Richardson差分法。问题：差分法。问题：1 1）观测资料缺乏，大气初始状态无法确定。）观测资料缺乏，大气初始状态无法确定。2 2）对于原始方程，差分的稳定性理论没有。）对于原始方程，差分的稳定性理论没有。3 3）计算速度无法解决。）计算速度无法解决。RichardsonRichardson实验失败的主要原因实验失败的主要原因 方案过于普遍化，对大气波动和数值计算中的一些基本理论问题认识不够。大气变化

21、实际上是多种波动的叠加，其中代表天气意义的变化的波动是那些波长比较长、变化比较缓慢的波动，他们可能是其他一些更短的波动的平均结果。短波意味着时空尺度小、移动变化快。因此，如果我们对预报的初值不加处理或处理不好，都会造成初始状态的不平衡，由此产生快速移动的重力波，随着时间积分的进行，重力波的发展掩盖了预报中气象信号的初始变率。而且Richardson的计算方案也没有考虑应满足计算稳定条件，因此计算结果只能是导致“计算崩溃”。重力波的振幅虚假地增大，完全掩盖了天气波动的发展。 观测网的极端稀缺和计算量的巨大也是阻碍Richardson成功的两大因素。当时一战前，人们对大气的了解主要还是依靠放置于地

22、面上的观测仪器所获得的。无线电通信、无线电探测、航空飞行都还不发达。对高层大气的了解基本是空白。只是依靠有限的高山站的数据推测得出。美国人曾经使用风筝携带仪器探测大气。受资料和通信的限制，实时预报是不可能的，即使只是根据事后收集的数据作实验计算，恶劣的初值和边界条件就足以导致计算失败。Richardson工作的意义 虽然，Richardson的实验结果使人沮丧，以至其后20多年里数值预报无人问津，但现代数值预报的发展确实是按照Richardson当年的思路在发展 。他不仅确立了现代数值预报的基本观念，也提出了解决这一问题的具体方法。 1950年Charney(美）等人数值预报的成功正是在充分认

23、识Richardson失败的教训基础上取得的。数值预报开始于上个世纪数值预报开始于上个世纪2020年代，以年代，以19221922年英国数学家年英国数学家RichardsonRichardson的的Weather Weather Prediction by Numerical ProcessPrediction by Numerical Process一书为标志。人们把一书为标志。人们把RichardsonRichardson的工作认的工作认为是现代数值预报的开始，也称为为是现代数值预报的开始，也称为数值预报发展的第一个里程碑。数值预报发展的第一个里程碑。 Courant、Friedrichs

24、和和Lewy， 1928提出差分方法的稳定性理论：即对于线性平流方程提出差分方法的稳定性理论：即对于线性平流方程在用中央差分法时，时间步长和空间步长必须满足：在用中央差分法时，时间步长和空间步长必须满足：, 1xtc这就是这就是CFL准则。数值预报仍无法实现：方程未简化！准则。数值预报仍无法实现：方程未简化！ 20 世纪40 年代前后大气科学取得了重大突破,人们揭示出了大气中存在着3 大波动:声波、重力波、天气慢波。Rossby 波是重要的天气学波动,是根据著名的Rossby大气长波理论提出来的。这些为数值天气预报滤波模式发展奠定了大气科学理论基础。到第二次世界大战后,地面和高空观测密度、范围

25、大大增加,并出现大容量、高速电子计算机,为数值天气预报模式发展又提供了可靠的初值条件和有力的计算手段与工具。Rossby 1939：长波理论：控制大气大尺度运动的方程是正压无辐散涡度方程：长波理论：控制大气大尺度运动的方程是正压无辐散涡度方程：,0fVt 长波理论解决了数值预报的理论问题，为数值预报的实现奠定了理论基础。 Charney 等在1950 年, 借助美国世界首台电子计算机ENIAC(Elect ronic Numerical Integrator and Computer) ,用滤掉(或不包括) 重力波和声波的准地转平衡(quasi- geost rophic) 滤波一层模式,成功

26、地制作出了500hPa 高度场形势24 h 预报,从而开创了数值天气预报滤波模式时代。图1世界上第一张成功的数值天气预报图：500mb高度场和涡度场(引自气象简编) 1950年，Charney等人计算出了历史上第一张数值预报天气图。这一结果的公布被人认是数值预报发展的第二个里程碑。 继Charney 等的成功之后,Rossby 返回欧洲瑞典领导一个研究小组,也成功地利用瑞典制造的、当时世界上强大的“BESK”计算机,再现了Charney 等的数值预报试验。4 年后的1954 年瑞典在世界上率先开始了业务(实时) 数值天气预报,较之美国开始业务数值天气预报早了6个月 。从这一年开始,数值天气预报

27、从纯研究探索走向了业务应用,同时也意味着地球科学首先由大气科学开始从定性研究向定量研究迈出了坚实的第一步。 准地转滤波模式对于研究认识副热带大尺度大气动力过程是很有用的,但是它太简化,它的精度不足于可以使数值天气预报研究应用不断发展,因而,用原始方程模式取而代之就成了最可能的选择。 要从滤波模式走到原始方程模式必须逾越两道障碍:一、是如何获取足够精度的初始水平散度场的问题,而水平散度又不是气象观测变量;二、是如何选择满足计算稳定条件的时间步长,这意味着若时间步长过短、对计算机能力要求过高而影响其可行性。 Charney 通过小时间步长、初始水平散度取为零的正压原始方程模式的成功试验证明了原始方

28、程模式用于数值天气预报中是可行的。 另外，Charney对非绝热和摩擦项、水汽凝结过程、辐射过程、湍流过程等物理过程的重要性和作用进行了一一讨论。随后的研究也越来越重视次网格物理过程的参数化影响问题。从20 世纪60 年代中期起,次网格物理过程参数化的重要性得到了确定,也逐步走向成熟。 与此同时,包含有简单物理过程参数化方案的、较完善的原始方程数值天气预报全球模式也在逐渐形成。1965 年,Smagorinsky 等提出了当时较高分辨率的9 层大气环流模式, 数值试验结果表明该模式的设计构造是成功的,这是数值天气预报模式业务应用10 年后,在数值天气预报模式设计上取得的重大突破,为现代数值天气

29、预报模式的研究与应用奠定了重要基础。 1975年欧洲科学技术合作计划开始，其中为加强数值天气预报很快成立了欧洲中期天气预报中心（European Centre for Medium-Range Weather Forecasts，简称ECMWF。第一个实时中期预报于1979年6月完成，1979年8月1日开始业务中期天气预报，标志着数值天气预报走向成熟。数值模式的先驱们经过了一个世纪赋有挑战性的数值模式的先驱们经过了一个世纪赋有挑战性的创新工作，奠定了现代数值模式发展的基础。创新工作，奠定了现代数值模式发展的基础。1904 -1950s：数值模式的酝酿期1950s-1960s：数值模式架构的完善

30、期1960s-1980s：大气环流模式的发展1980s-1990s：完善物理过程和资料同化1990s之后 ：耦合模式的发展基于物理规律的数值预报理论的发展，使人类可以利用计算机重现或预测发生在自然界的天气变化过程。这是地球科学由“定性”走向“定量”的重大进步。数值预报理论和技术水平越来越高，应用领域越来越广泛，数值预报技术已被认为是未来解决天气预报、气候预测等问题的根本科学途径。数值天气预报理论和应用是过去一个多世纪以来地球科学的最重大进步和成就之一。数值模式是天气预报不可缺少的和气候预测最具潜力的工具。数值模拟是加强对天气和气候系统演变规律认识所不可替代的工具，是提高理论认识的重要途径。对天

31、气气候演变规律的准确认识都会在数值模式中得到体现。数值模式的整体水平是一个气象强国气象综合水平的集中体现。数值模式的重要地位数值模式的重要地位促进数值预报迅速发展的因素促进数值预报迅速发展的因素：探空技术及先进的探测技术的发展为预报提供了必要食粮。 通讯技术的发展为业务预报资料的收集提供的必要的手段。 动力气象和天气学的发展为各种模式的发展提供了理论。 计算机和计算技术的发展是NWP实现的必备工具。 1.预报范围扩大：时间上：短期、中期、长期、气候、古气候恢复。空间上：局地、全球、空间2.大气中的各种物理过程的描述更加细致。地形、辐射、行星边界层、积云对流、海气相互作用、微量气体。3.计算理论

32、不断改进、资料处理更加完善。4.预报准确率不断提高、产品更加丰富。因此，数值预报已成为现代天气预报的基础。气候模式发展近几年的发展趋势近几年的发展趋势: :精细的短期预报。用风暴尺度模式预报灾害性天气（提供有预报技巧的强对流天气预报）较完善的资料同化系统。可以最大限度地从观测中提取有用信息，尤其是可以把卫星雷达等非直接观测的遥感信息加入到资料同化系统中。改进中期预报，尤其是通过使用集合预报方法改进中期预报的水平。全球模式的分辨率将可提高至1015 km ,可预报天数将比目前提高12 d ;有限区域短期天气预报模式的分辨率将可提高至210 km ,可提前1248 h 提供逐小时的滚动天气预报。对

33、由集合预报指示的预报误差快速增长区，增加观测。发展更为完善的海-地-气耦合模式。发展各种更高效的计算方法。为政府决策及公众服务提供更多、更细致、更广泛的预报参考。诸如空气污染、公路状况、河流洪水预报、污染物扩散、健康保险等。这里提到的许多方面都是属于非传统的气象服务领域。而数值预报是我们可以依靠的最主要手段。kmms1hr1d1m1y10yseasonalclimatesynoptic临近预报短期中期长期季、年气候typhoon 现代数值天气预报正向两方面发展：一方面试图解决大尺度、全球、季、气候预报；另一方面，局地的、小尺度问题。 中国是气象数值预报起步较早的国家之一,早在20 世纪5060

34、 年代初,中国第一代数值预报专家就开始了数值预报模式及相关计算方法的研究,并建立了试验预报系统。但从总体上说,数值预报的发展缓慢。 改革开放以后,中国的气象数值预报才进入了真正大发展的时期。1980 年国家气象中心利用中国自行研制的亚欧区域短期预报模式(简称为A 模式) ,开始发布日常48 h 形势预报,标志着中国气象数值预报进入业务实用阶段。“六五”期间由国家气象中心、中国科学院大气物理研究所、北京大学组成的联合数值预报中心建立了北半球和亚洲区域模式系统(简称为B 模式) ,将中国数值预报业务向前推进了一步。中期数值预报模式中期数值预报模式（李泽椿、李泽椿、 纪立人纪立人等）等）2层全球格点

35、模式层全球格点模式（曾庆存、梁信忠曾庆存、梁信忠等）等）4层海洋环流模式层海洋环流模式（张学洪张学洪等）等）HLAFS 模式模式（张玉玲、郭肖容张玉玲、郭肖容等）等）p- 混合坐标模式混合坐标模式（钱永甫钱永甫等）等）YH模式模式（颜宏颜宏等）等）TL系列模式系列模式（薛纪善薛纪善等）等）REM -坐标坐标（宇如聪、曾庆存宇如聪、曾庆存等）等）海气耦合模式海气耦合模式（曾庆存、张学洪曾庆存、张学洪等）等）短期气候预测动力模式系统（短期气候预测动力模式系统（BCC，IAP等）等）耦合气候系统模式（耦合气候系统模式（IAP等）等）国内国内80-90年代的数值模式发展年代的数值模式发展国内自主研制的

36、非静力平衡中尺度模式不多,且多限于研究模式1991 1992 1995 1996 1997 1998120h48h168h48h240h48h320km/9L200km /15L200km /16L110km /15L120km /19L55km /20LOIOIOI OI OI OICYBER992CYBER992CRAY-C92CRAY-C92CRAY-C92CRAY-C92200220032004240h48h48h60km 60km /31L/31L27km 27km /20L/20L20km-20km-1km1kmOIOIOIOIOIIBM-SPIBM-SPIBM-SPIBM-SPI

37、BM-SPIBM-SP时间时间模式模式时效时效分辨率分辨率分析同化分析同化计算机计算机80年代年代(AB模式模式)90年代年代(全球、区域全球、区域)我国目前主要运行的业务模式T213L31(全球中期预报谱模式)T639L60HLAFSMM5(有限区域格点模式)WRF(MM5升级）Grapes（有中国自主知识产权的格点模式）T代表采用三角形谱系数截断，沿纬圈213的波，垂直分31层研究成果：开发建立了我国新一代全球研究成果：开发建立了我国新一代全球/ /区区域一体化多尺度通用数值预报系统域一体化多尺度通用数值预报系统GRAPESGRAPES集合预报集合预报 GRAPES的应用的应用大气科学研究

38、大气科学研究沙尘暴预报沙尘暴预报实时业务预报实时业务预报观测试验模拟观测试验模拟 热带气旋预报热带气旋预报闪电模拟闪电模拟 教学培训教学培训 气候模拟气候模拟 GRAPESGRAPES“家族家族”GRAPES_Meso GRAPES_Meso 中尺度预报模式中尺度预报模式GRAPES_GlobalGRAPES_Global全球预报模式全球预报模式GRAPES_TCMGRAPES_TCM台风预报模式台风预报模式GRAPES_TMMGRAPES_TMM热带气象预报模式热带气象预报模式GRAPES_DAMGRAPES_DAM沙尘气溶胶预报模式沙尘气溶胶预报模式GRAPES_SWIFTGRAPES_S

39、WIFT临近预报模式临近预报模式GRAPES_MEPSGRAPES_MEPS中尺度集合预报模式中尺度集合预报模式GRAPES_LMGRAPES_LM雷电研究模式雷电研究模式主要数值天气预报业务体系主要数值天气预报业务体系各系统的侧重功能各系统的侧重功能模式预报的主要侧重点T639L60T639L60预报预报1-10天天内全球内全球大尺度环流大尺度环流的演变、的演变、大范围降水大范围降水的的发生发展等，如高空槽脊、副热带高压，地面高、低发生发展等，如高空槽脊、副热带高压，地面高、低气压和冷、暖锋面，大范围雨带等气压和冷、暖锋面，大范围雨带等GRAPES-MesoGRAPES-Meso对对中国区域

40、的降水中国区域的降水进行短期（进行短期（60小时以内小时以内）预报，）预报，如降水发生区域、降水强度、降水出现时段等如降水发生区域、降水强度、降水出现时段等 台风模式台风模式专门针对发生在西太平洋专门针对发生在西太平洋-我国沿海的台风我国沿海的台风，进行其，进行其中心位置、移动路径的中短期（中心位置、移动路径的中短期（96小时内小时内）预报）预报 中期中期T213T213集合预报集合预报对未来对未来1-10天，特别是第天，特别是第5-10天期间、全球大尺度环天期间、全球大尺度环流和大范围降水发展演变的可能性（即流和大范围降水发展演变的可能性（即概率概率）进行预）进行预报报 中尺度中尺度( (基

41、基于于MM5,27-MM5,27-9km)9km)预报预报京津京津及其它地区及其它地区气象要素气象要素短时间（短时间（48小时小时）的）的连续变化，如气温、风、相对湿度、降水量等间隔连续变化，如气温、风、相对湿度、降水量等间隔3小时的预报小时的预报 到目前为止到目前为止我国的业务数值预报模式系统仍然是以引进国外为主！我国天气预报所用的数值预报产品仍然以国外的模式结果为主！我国自主发展的数值模式中针对我国特殊气象问题的物理过程还很少！我国数值模式的资料同化能力很低！中国中国-D3中国中国-D5中国中国-D7、 分别T213L312003年7月份3、5、7天预报500hPa距平相关系数欧洲中心业务

42、模式欧洲中心业务模式500hPa的距平相关系数的距平相关系数5天天的形势预报准确率达的形势预报准确率达85%以上！可以上！可用预报达用预报达一周一周！ 集合预报是国际上二十世纪九十年代初投入业务运行的数值天气预报新技术。集合预报1969 年首先由Epstein 从大气运动的随机性角度提出的，其理论基础是蒙特卡罗统计试验法。 传统的数值预报是将观测值作为初值，通过求解大气动力学方程和热力学方程组，就可以客观、定量地做出天气预报。 但大气探测的误差是非常大的。地球表面71% 是海洋，在广阔的洋面上几乎没有观测站；在2 9 % 的陆地上，高原和沙漠以及南极洲上也几乎没有观测站，再加上观测仪器本身的误

43、差以及资料分析、同化处理中引入的误差等。数值预报用的初值非常不准确。以不准确的初始场做数值预报，是造成预报不准确的一个重要原因。 集合预报方法是把数个由少许不同的初值计算得到的数值预报加以平均，作为预报。其基本思想是通过平均，减少由上述的误差源产生的各自决定论的预报误差部分，而保留有意义的情报（信息），从而使预报的均方根误差与各自预报的均方根误差相等或更小。 集合预报突破了传统数值预报观念，认为大气的初值是不确定的，在初始场上加少则几种、多则几百种各种小扰动，使得初始场成为某种概率密度函数。所以集合预报被认为开创了数值预报的新纪元。一个好的数值预报系统应是具备高质量观测资料、好的同化系统、好的

44、预报模式、强大的计算机、好的集合预报系统。 70-80年代的研究探索和数值试验后，至90年代初，NCEP和ECMWF于1992年12月先后建立了各自的集合数值预报业务系统,这标志着集合数值预报已进入了实际业务应用的成熟阶段。现在，在美国、欧洲等发达国家的数值预报业务体系中，集合预报已占据了非常重要的位置（如ECMWF的集合预报产品已占据所有数值预报产品的60%以上），而且集合预报产品已在气象台站得到了广泛地应用。我国集合数值天气预报系统建立了多套集合数值天气预报系统：建立了多套集合数值天气预报系统：u1996年5月建立了基于T63L16模式的全球中期全球中期集合预报系统（采用时间滞后平均法（L

45、AF）生成扰动初值）； u1999年建立了32个成员的基于T106L19模式的全球中期集合预报系统（采用奇异向量法（SVs）生成扰动初值），并于2000年3月投入业务应用；u2006年12月15个成员的新一代全球中期全球中期T213L31集合数值天气预报系统（采用增长模繁殖法（BGM）生成扰动初值）实现了准业务化运行 。u已研发基于WRF的区域短期区域短期集合数值天气预报系统，提供奥运气象服务。 全球全球面条图面条图全球区域一体化的集合预报全球区域一体化的集合预报产品加工和图形显示系统产品加工和图形显示系统谱格变化、高斯格点转化谱格变化、高斯格点转化GRIB格式格式侧重强对流指数产品侧重强对流

46、指数产品全球模式输出量：全球模式输出量：谱系数、规约高斯格点场区域区域区域模式输出量：区域模式输出量：标准经纬格点场平均值平均值+发散度图发散度图概率图概率图烟羽图烟羽图集合预报产品大致分三类：1 1、集合平均或集合中值预报，是集合预报最初级的应用。、集合平均或集合中值预报，是集合预报最初级的应用。一般可能比单个预报准确，但对大气不稳定而可能出现一般可能比单个预报准确，但对大气不稳定而可能出现的分叉且多平衡态无能为力。的分叉且多平衡态无能为力。2 2、大气预报的可信度预报，通常用集合预报成员间的发、大气预报的可信度预报，通常用集合预报成员间的发散程度来度量。如面条图。散程度来度量。如面条图。3

47、 3、概率预报，概率分布包含了该集合预报系统所能提供、概率预报，概率分布包含了该集合预报系统所能提供的所有信息，最大程度地包含了实际大气可能发生的各的所有信息，最大程度地包含了实际大气可能发生的各种情况。种情况。天气预报的挑战大气是个大气是个浑沌浑沌 ( (Chaotic)Chaotic) 系统系统线性线性是指量与量之间的正比关系，用直角坐标形象线性是指量与量之间的正比关系，用直角坐标形象地画出来，是一根直线。在线性系统中，部分之和地画出来，是一根直线。在线性系统中，部分之和等于整体，描述线性系统的方程遵从叠加原理，即等于整体，描述线性系统的方程遵从叠加原理，即方程的不同解加起来仍然是解。方程

48、的不同解加起来仍然是解。 非线性非线性则指整体不等于部分之和，叠加原理失效，非线性则指整体不等于部分之和，叠加原理失效，非线性方程的两个解之和不再是方程的解。非线性方程的两个解之和不再是方程的解。非线性系统非线性系统的主要特征的主要特征非线性系统非线性系统的主要特征是当系统中的参量发生极微的主要特征是当系统中的参量发生极微小变化时，在一些关节上，可以引起系统运动形式小变化时，在一些关节上，可以引起系统运动形式的定性改变，在对外界激励的响应上，则表现为出的定性改变，在对外界激励的响应上，则表现为出现与外界激励有本质区别的行为现与外界激励有本质区别的行为 混沌理论的本质，简单的个体遵循简单的规律，

49、混沌理论的本质，简单的个体遵循简单的规律，相互作用可以建立复杂和不可预测的行为。相互作用可以建立复杂和不可预测的行为。浑沌系统的特性浑沌系统的特性不是无规律的变化，也不是完美的规律性变化初始很小的不确定性，随时间演进将成为后续变化很大的不确定性由一定的非線性作用導致，在確定性系統中 出現極其複雜、貌似無規的運動。93LorenzLorenz奇异吸引子奇异吸引子(strange attractor)(strange attractor)轨迹图轨迹图Lorenz動力方程式 dx / dt = -(x y) dy / dt = -xz + x y dz / dt = xy bz x, y, z :

50、速度、溫度、溫度梯度 ,b : 確定的控制參數運動軌跡 確定性 繞A、B兩點 隨機性 圈數、大小94从三个不同角度从三个不同角度(x-y-z)(x-y-z)看到的看到的LorenzLorenz奇奇异吸引子异吸引子進入混沌倍週期分岔進入混沌 以生物種群繁衍為例 Xn+1 = aXn ( 1 Xn ) 0 Xn 0 Xn：某一年的數量 ； a：控制參數，生殖率進入混沌混沌中的普適性(Feigenbaum） - 結構普適性 單峰映射函數 進入混沌的過程具有相同的分岔 結構。 - 測度普適性 寬度比極限值2.5029 間距比極限值4.6692 碎形倍週期分岔中的對初始條件的敏感依賴性差之毫釐，失之千里