ECMWF 数值预报模式简介ppt课件.ppt

《ECMWF 数值预报模式简介ppt课件.ppt》由会员分享，可在线阅读，更多相关《ECMWF 数值预报模式简介ppt课件.ppt（46页珍藏版）》请在三一办公上搜索。

1、,ECMWF 数值预报模式简介,气象预报的四个方面问题,当前的天气或气候信息完备的综合观测系统完善的观测系统；资料信息识别；资料的综合处理（同化）天气或气候的演变规律从资料得到新认识从资料得到新认识；反映大气运动数学物理规律的微分方程组。外力和强迫的变化地形和边界强迫；太阳常数；引力等从已知预报未来的手段完备的数值模 式；承载数值模式计算、显示和通讯平台 驾驭观测、资料分析和数值模式发展和应 用的人才队伍,大气圈,水圈,冰晶圈/冰冻圈,生物圈,岩石圈,人类圈,纷繁多样数值预报产品-什么会有差别？,分辨率不同 模式的表述的物理过程有差异 计算方案、网格、变量分布等再分析资料分析的初值能代表大气实

2、况吗？,Orography and Resolutions,GSM T213 (60km),RSM(20km),MSM(10km),Orographic effects are better captured by higher resolution models. The surface parameters such as Tsurf might be predicted more realistically by those models.,ECMWF 致力于2016-2025 数值预报发展战略： 利用集合预报方法提前2周预报高影响天气事件 提前4周，无缝隙地预报大尺度系统形势和系统移动

3、 提前1年预测全球范围的异常状况,研究和更丰富的知识基于集合预报的分析和预报技术如何实现目标：Observations High resolution ensemble Earth-system Scalability Funding People（引自 Erland Klln，“Earth system modelling for seamless prediction” ECMWF Annual Seminar 2016）,Future Earth System model and assimilation developmentsErland Klln, ECMWF,Ocean mode

4、l (NEMO),Wave model (ECWAM),Atmosphericmodel (IFS),Air density,Gustiness,Neutral wind,Roughness,Currents,Sea surface temperature,All configurations,Ensemble FC,Solar and non solar fluxes, E-P,Towards a fully coupled system(currently only operational in EPS),ORCA1_Z42,TCo1279/TCo639,14km/28km,9km/18k

5、m,Ensemble systems only:- Medium range forecast - Monthly forecast- Seasonal forecast,EveryIFStime step,Everycouplingtime step(1 or 3 hours),Single executable,Ocean model (NEMO),Wave model (ECWAM),Atmosphericmodel (IFS),Air density,Gustiness,Neutral wind,Roughness,Turbulent energy,Stokes drift,Curre

6、nts,Sea surface temperature,All configurations,Stress,Ensemble FC,Solar and non solar fluxes, E-P,Towards a coupled system,ORCA1_Z42,TCo1279/TCo639,14km/28km,9km/18km,Single executable,Operationalfrom day 0 since 2013,Ocean model (NEMO),Ice concentration,Wave model (ECWAM),Atmosphericmodel (IFS),Air

7、 density,Gustiness,Neutral wind,Roughness,Turbulent energy,Stokes drift,Currents,Ice model (LIM),Sea surface temperature,Ice concentration,All configurations,Stress,Ensemble FC,future operational,Solar and non solar fluxes, E-P,Towards a coupled system,ORCA0.25_Z75,Addingactive sea ice model,Impleme

8、ntation:End 2016.,Single executable,Calling sequence of the single executable,14,Simplified flow chart of the coupled model, here two time steps are shown.In reality, IFS/WAM coupling every IFS time step, butCall to NEMO every hour (or 3 hours) with averaged accumulated fluxes.,Impact of Resolution

9、on tropical cyclone forecast,For instance Typhoon Haiyan: forecasts from 4th, 5th and 6th November 2013, 0 UTCall from operational analysis.,Black: estimated from observationsRed: old operational Ensemble resolution (32 km)Blue: old operational HRES configuration (16 km)Green: experimental: new HRES

10、 configuration (10km),Typhoon Haiyan at peak intensity on November 7, 2013,Impact of Coupling on tropical cyclone forecast,For instance Typhoon Neoguri: forecasts from 6 July 2014, 0 UTC,Black: estimated from observationsGreen: old operational HRES configuration (uncoupled) (16km)Red: experimental:

11、16km coupled to NEMO (ORCA025_Z75)Blue: 16km coupled to NEMO + new physics,minMSLP (hPa),900,950,Neoguri affecting Okinawa on July 8, 2014,Climate reanalyses for the coupled earth model,ECMWF coupled Earth model for medium-range weather forecasting,New coupled assimilation system (CERA) for the coup

12、led Earth model:atmospheric and ocean observations assimilated simultaneouslyocean observations can impact atmospheric estimate and converselyCERA-20C reanalysis in production (1900-2010),Karl and Trenberth 2003,Coupled assimilation system (CERA),EDA variational approach with a 24-hour window that a

13、ssimilates simultaneously atmospheric and ocean observations,coupled model computes observation misfits in each outer iteration,atmospheric and ocean increments are computed in parallel to correct the initial state,SST computed in NEMO and constrained by relaxation,analysis dynamically consistent wi

14、th respect to the coupled model,Conclusions:,ECMWF has a coupled atmosphere-wave-ocean circulation forecasting system, currently operational in the Ensemble Prediction System.Work is ongoing on using a higher resolution ocean components (ORCA025z75) planned for end of 2016 in the Ensemble forecasts

15、and later in the High resolution system.There is a clear benefit in coupling the different models, but it creates new challenges in determining what physical parameters need to be exchanged. Furthermore, model parameterisations might need revisiting.,EC模式动力框架,垂直采用坐标，U=ucos ；V=vcos,对比球坐标系中的基本方程组,曲率项,

16、重力不是指向地心，等位势面也不是球面。,非静力学模式 静力学模式,静力学平衡：运动的垂直尺度远小于水平尺度的情况下成立,模式的分类,基于连续方程的模式分类,非静力学模式根据模式中是否包含声波，分为： 滞（非）弹性（anelastic)模式 假定大气不可压，滤掉声波 弹性(elastic)模式 大气可压、由散度预报气压的变化，但声波 需要特殊的处理,根据连续方程的近似程度，滞弹性和弹性模式又可进行不同的分类：,滞弹性模式,不可压缩模式,滞弹性模式,弹性模式,准可压模式,完全可压模式,数值计算方案,大气运动方程组是一套复杂的非线性方程组，目前，还没有理论解存在，只能借助数值方法求解。数值方法有很多

17、，目前气象上用的主要用的是差分方法、谱方法、有限元法等。差分方法：就是在离散的网格点上，以差商代替微商，求解微分方程的方法。,关于有限差分与谱方法的图例,有限差分,谱方法,| | | n-1 n n+1,从物理上考虑，差分的精度还与采样点密即网格，与要素在采样区的变化分布有关。即采样数据是否能有代表性。,中央差,精确值,采样值,混淆误差,离散网格所能表示的最短波长上表中当L/ x =2时，即两倍格距波，R=0，误差100%，因此，离散网柜完全不能表示两倍格距以下的短波。如格距为100km，则你只能预报出200km以上的天气系统。,1.5倍格距波,3倍格距波,最短的波，2倍格距波,1.5倍格距波

18、被当成3倍格距的波混淆误差会把小于2倍格距的波歪曲为2倍格距以上的波。因此混淆误差主要集中在24倍格距间。预报员会在什么情况下遇到“混淆误差”？,拉格朗日和半拉格朗日方法,拉格朗日和半拉格朗日方法是从另外一个角度来求取平流方程的解。 有解析解：不失一般性，考虑 情况下平流方程的解。,可以证明下述两个式子是成立的：,Think of measuring the temperature:Lagrangian: in a balloon, floating with the windEulerian: on the ground, wind blows past,半拉格朗日方法的精度依赖于DP的精度，还依赖于值路径,球坐标系中运动方程的分量形式,二、热力学方程和连续方程,连续方程,GRAPES方程组,GRAPES方程组,引入参考大气：,线性项分解：（U方程为例）：,