Application of the Correction Algorithm to Mass Conservation in GRAPES_GFS

Su Yong; Shen Xueshun; Zhang Qian

doi:10.11898/1001-7313.20160603

Vol.27, NO.6, 2016

Turn off MathJax

Article Contents

Article Navigation > Journal of Applied Meteorological Science > 2016 > 27(6): 666-675

Su Yong, Shen Xueshun, Zhang Qian. Application of the correction algorithm to mass conservation in GRAPES_GFS. J Appl Meteor Sci, 2016, 27(6): 666-675. DOI: 10.11898/1001-7313.20160603.

Citation:

Su Yong, Shen Xueshun, Zhang Qian. Application of the correction algorithm to mass conservation in GRAPES_GFS. J Appl Meteor Sci, 2016, 27(6): 666-675. DOI: 10.11898/1001-7313.20160603.

Su Yong, Shen Xueshun, Zhang Qian. Application of the correction algorithm to mass conservation in GRAPES_GFS. J Appl Meteor Sci, 2016, 27(6): 666-675. DOI: 10.11898/1001-7313.20160603.

Citation:

Su Yong, Shen Xueshun, Zhang Qian. Application of the correction algorithm to mass conservation in GRAPES_GFS. J Appl Meteor Sci, 2016, 27(6): 666-675. DOI: 10.11898/1001-7313.20160603.

PDF( 5847 KB)

Application of the Correction Algorithm to Mass Conservation in GRAPES_GFS

DOI: 10.11898/1001-7313.20160603

Su Yong^1,2,3,
Shen Xueshun^{3
,
,},
Zhang Qian⁴

1.
Chinese Academy of Meteorological Sciences, Beijing 100081
2.
Department of Atmospheric Sciences, Nanjing University of Information & Technology, Nanjing 210044
3.
Numerical Prediction Center of CMA, Beijing 100081
4.
CMA Center for Communication and Outreach, Beijing 100081

Received Date: 2016-04-21
Rev Recd Date: 2016-07-13
Publish Date: 2016-11-30

Abstract

Abstract

The conservation of mass is very important for the dynamic core of numerical model, especially for climate simulation or medium and long term prediction. For the traditional semi-Lagrangian dynamic core, it is difficult to satisfy the mass conservation theoretically, the finite volume method can be used in the semi-Lagrangian continuous equation to satisfy the mass conservation in theory, but the practical application is more complex, and no operational center adopt the method at this time. On the other hand, the traditional semi-Lagrangian method combined with a simple and easy mass correction algorithm is also a good choice.The GRAPES_GFS (Global-Regional Assimilation and PrEdiction System, Global Forecast System) of Numerical Prediction Center of CMA (China Meteorological Administration) faces a problem that the mass loss is obviously in the process of long term integration. The global mean sea level pressure drops about 1 hPa in 10-day run, due to the choice of the continuity equation and the calculation accuracy of the dynamical core. In the early times of development, there is a very simple mass correction algorithm in GRAPES_GFS, however, in the process of batch testing of the system in 2013, it is found that this algorithm will lead to the increase of the height bias at the top level of the model, so the algorithm is closed.Ideas of C-CAM (Climate-Community Atmosphere Model) model are drawn on, to achieve mass conservation through the correction of surface pressure. A method is developed adjusting mass in each grid box according to different weight coefficient, with the basic idea of modifying more in the grid boxes which have larger changes in mass while less in those with smaller changes. The feasibility of this method is verified by a series of experiments, including ideal test, real case prediction and a batch of cycle prediction, illustrating that the method can effectively reduce the bias of height field under the premise of ensuring the mass conservation, and alleviate the problem of underestimating weather systems. It shows that this method has a certain application value in actual operational forecast, and be further studied in calculating the weight coefficient of each point. It can reduce the bias of the height field at the top level of the model, as well as the high latitude in the Southern Hemisphere.
- GRAPES_GFS;
- mass conservation;
- semi-Lagrangian

FullText(HTML)

References(19)

Figures and Tables

Fig. 1 The change of total mass during the 180-day integration

DownLoad: Full-Size Img PowerPoint

Fig. 2 The change of total mass during 30-day integration

DownLoad: Full-Size Img PowerPoint

图 3 Sen_phyon试验积分1 h时刻单步调整的Exner气压 (填色)

(a) 模式面第15层 (箭头为水平风场, 单位：m·s^-1), (b) 纬向平均的垂直剖面 (线条为垂直速度场, 单位：m·s^-1)

Fig. 3 The magnitude of the adjusted Exner function for the Sen_phyon test after 1-hour integration (the shaded)

(a) the fifteenth model layer (arrow is the horizontal wind field, unit：m·s^-1), (b) vertical profiles of zonal mean (contour is the vertical velocity field, unit:m·s^-1)

DownLoad: Full-Size Img PowerPoint

图 4 积分72 h后试验Ctr_phyon和Sen_phyon差异的纬向平均垂直剖面

(a) 高度场，(b) 温度场

Fig. 4 The vertical profile of zonal mean of the difference between Ctr_phyon and Sen_phyon after 72-hour integration

(a) height field, (b) temperature field

DownLoad: Full-Size Img PowerPoint

图 5 北半球高度场均方根误差平均值随时间演变 (单位: gpm)

(a) Ctr2, (b) Sen2相对于Ctr2的差异

Fig. 5 The evolution of the height field mean square error in the Northern Hemisphere with the time (unit:gpm)

(a) Ctr2, (b) difference of Sen2 with respect to Ctr2

DownLoad: Full-Size Img PowerPoint

图 6 高度场第5天的预报相对于FNL分析场偏差的纬向平均垂直剖面 (单位: gpm)

(a) Ctr2, (b) Sen2

Fig. 6 The zonal mean vertical profile of the height bias for the fifth day with respect to the FNL analysis (unit:gpm)

(a) Ctr2, (b) Sen2

DownLoad: Full-Size Img PowerPoint

Fig. 7 500 hPa height bias for the fifth day with respect to the FNL analysis (a) Ctr2, (b) Sen2

DownLoad: Full-Size Img PowerPoint

图 8 高度场的预报相对于分析场的距平相关系数

(a) 北半球, (b) 南半球

Fig. 8 The hight anomaly correlation coefficient for the cycle test

(a) the North Hemisphere, (b) the South Hemisphere

DownLoad: Full-Size Img PowerPoint

Fig. 9 Averaged 5-day forecast of 500 hPa height field (unit:dagpm)

DownLoad: Full-Size Img PowerPoint

Table 1 Instructions for the three groups of tests

试验分类	试验名称	积分时长/d	试验说明
理想试验	Ctr1	180	Rossby_Haurwitz波
理想试验	Sen1	180	Ctr1及质量订正
实际个例预报试验	Ctr_phyoff	30	实际个例，关闭物理过程
	Ctr_phyon	30	实际个例，打开物理过程
	Sen_phyoff	30	Ctr_phyoff及质量订正
	Sen_phyon	30	Ctr_phyon及质量订正
批量循环预报试验	Ctr2	8	批量循环预报试验
批量循环预报试验	Sen2	8	Ctr2及质量订正

DownLoad: Download CSV

References

[1]	Staniforth A, Cote J.Semi-Lagrangian integration schemes for atmospheric models-A review.Mon Wea Rev, 1991, 119:2206-2223. doi: 10.1175/1520-0493(1991)119<2206:SLISFA>2.0.CO;2
[2]	Thuburn J.Some conservation issues for the dynamical cores of NWP and climate models.J Comput Phys, 2008, 227:3715-3730. doi: 10.1016/j.jcp.2006.08.016
[3]	Diamantakis M.The Semi-lagrangian Technique in Atmospheric Modeling:Current Status and Future Challenges.ECMWF Seminar in Numerical Methods for Atmosphere and Ocean Modeling, 2013.
[4]	Ritchie H, Tanguay M.A comparison of spatially averaged eulerian and semi-lagrangian treatments of mountains.Mon Wea Rev, 1995, 124:167-181.
[5]	陈德辉, 沈学顺.新一代数值预报系统GRAPES研究进展.应用气象学报, 2006, 17 (6):773-777. doi: 10.11898/1001-7313.20060614
[6]	胡江林, 沈学顺, 张红亮, 等. GRAPES模式动力框架的长期积分特征.应用气象学报, 2007, 18 (3):276-284. http://qikan.camscma.cn/jams/ch/reader/view_abstract.aspx?file_no=20070349&flag=1
[7]	陈德辉, 杨学胜, 张红亮, 等.多尺度非静力通用模式框架的设计策略.应用气象学报, 2003, 14 (4):452-461. http://qikan.camscma.cn/jams/ch/reader/view_abstract.aspx?file_no=20030456&flag=1
[8]	伍湘君, 金之雁, 黄丽萍, 等.GRAPES模式软件框架与实现.应用气象学报, 2005, 16 (4):539-546. doi: 10.11898/1001-7313.20050415
[9]	刘艳, 薛纪善, 张林, 等.GRAPES全球三维变分同化系统的检验与诊断.应用气象学报, 2016, 27 (1):1-15. doi: 10.11898/1001-7313.20160101
[10]	苏勇, 沈学顺, 张倩, 等.应用样条插值提高GRAPES模式物理过程反馈精度.应用气象学报, 2014, 25 (2):202-211. doi: 10.11898/1001-7313.20140210
[11]	Jablonowski C, Lauritzen P, Nair R, et al.Idealized Test Cases for the Dynamical Cores of Atmospheric General Circulation Models:A Proposal for the NCAR ASP 2008 Summer Colloquium.2008.
[12]	Zerroukat M, Wood N, Staniforth A, et al.An inherently mass-conserving semi-implicit semi-lagrangian discretization of the shallow water equations on the sphere.Q J R Meteorol Soc, 2009, 135:1104-1116. doi: 10.1002/qj.v135:642
[13]	Zerroukat M, Wood N.SLICE-S:A semi-Lagrangian inherently conserving and efficient scheme for transport problems on the sphere.Q J R Meteorol Soc, 2004, 130:2649-2664. doi: 10.1256/qj.03.200
[14]	Lauritzen H, Nair D, Paul A, et al.A conservative semi-Lagrangian multi tracer transport scheme (CSLAM) on the cubed-sphere grid.J Comput Phys, 2009, 229:1401-1424. http://adsabs.harvard.edu/abs/2010JCoPh.229.1401L
[15]	Diamantakis M, Flemming J.Global Mass Fixer Algorithms for Conservative Tracer Transport in the ECMWF Model.Technical Memorandum, 2013. http://www.researchgate.net/profile/Johannes_Flemming/publication/262765450_Global_mass_fixer_algorithms_for_conservative_tracer_transport_in_the_ECMWF_model/links/0c960538d8908cc9a5000000.pdf?disableCoverPage=true
[16]	苏勇, 沈学顺, 彭新东, 等.PRM标量平流方案在GRAPES全球预报系统中的应用.大气科学, 2013, 37 (6):1309-1325. doi: 10.3878/j.issn.1006-9895.2013.12164
[17]	Priestley A.A quasi-conservative version of the semi-Lagrangian advection scheme.Mon Wea Rev, 1993, 121:621-629. doi: 10.1175/1520-0493(1993)121<0621:AQCVOT>2.0.CO;2
[18]	Bermejo R, Conde J.A conservative quasi-monotone semi-Lagrangian scheme.Mon Wea Rev, 2002, 130:423-430. doi: 10.1175/1520-0493(2002)130<0423:ACQMSL>2.0.CO;2
[19]	McGregor L.C-CAM:Geometric Aspects and Dynamical Formulation.CSIRO Atmospheric Research Technical Paper 70, 2005.