Vol.25, NO.2, 2014
Turn off MathJax
Article Contents
Article Navigation >  Journal of Applied Meteorological Science  > 2014  >  25(2): 202-211
Su Yong, Shen Xueshun, Zhang Qian, et al. Application of spline interpolation to physical process feedback accuracy improvement of GRAPES model. J Appl Meteor Sci, 2014, 25(2): 202-211.
Citation: Su Yong, Shen Xueshun, Zhang Qian, et al. Application of spline interpolation to physical process feedback accuracy improvement of GRAPES model. J Appl Meteor Sci, 2014, 25(2): 202-211.
shu

Application of Spline Interpolation to Physical Process Feedback Accuracy Improvement of GRAPES Model

  • 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

  • 5.

    Department of Hydrographical Meteorology, Nanjing University of Information & Technology, Nanjing 210044

  • Received Date: 2013-01-15
  • Rev Recd Date: 2013-12-20
  • Publish Date: 2014-03-31
    Abstract
  • The variable distribution in the vertical direction of GRAPES model's dynamic core adopts Charney-Phillips method. Vertical velocity, potential temperature, water substance are calculated at the whole layer, horizontal velocity and dimensionless pressure are calculated at the half layer, but in physical process, all the variables are placed on the half layer. In order to satisfy the needs of the central difference calculations and better representation of the physical processes in the boundary layer, a nonuniform stratification is adopted, which is dense near the ground, and the higher the more sparse. Therefore, in GRAPES model, linear interpolation is needed to convert variables between whole and half layers before and after the physical process calculation.For the weather prediction model of various international centers, Lorenz layers are used in the physical part and all the variables are on the half layer. Most models also use Lorenz layers in the dynamic core, except for the Unified Model of the UK Meteorological Office, which chooses Charney-Philips layer for dynamic core and uses linear interpolation in dealing with the similar problem of interpolation between whole and half layers.Linear interpolation is relatively simple, but the accuracy is not high, and it will cause deviation especially for lower and higher layers. The cumulative deviation in the temperature and humidity fields will further impact the height and wind fields. In addition, the interpolation process of water substance is also required to ensure monotonic, but the traditional cubic spline interpolation, polynomial interpolation, cannot be guaranteed monotonic, which will bring negative water, instability and other issues.In order to solve the problems above, the traditional cubic spline interpolation method is introduced for potential temperature interpolation in GRAPES model. After some special handling of the boundary value based on the traditional one, a monotonic cubic spline interpolation method is established for water substance, by which the forecast error of potential temperature and humidity fields in the GRAPES model is effectively reduced. The feedback accuracy of physical process is improved, and the model comprehensive performance is also enhanced.
    • FullText(HTML)
    References(21)

  • Fig. 1  Three different algorithms on the sinusoidal wave simulation

    DownLoad: Full-Size Img PowerPoint

    Fig. 2  Three different algorithms on the given distribution simulation

    DownLoad: Full-Size Img PowerPoint

    Fig. 3  The deviation vertical distribution after interpolation at the sample point (30°N, 115°E) over the Yangtze River Plain (a) the potential temperature, (b) the water vapor

    DownLoad: Full-Size Img PowerPoint

    Fig. 4  The regional averaged root mean square error vertical distribution after interpolation of the Northern Hemisphere mid-latitude region (20°—60°N)

    (a) the potential temperature, (b) the water vapor

    DownLoad: Full-Size Img PowerPoint

    Fig. 5  The zonal mean vertical distribution of temperature root mean square error after interpolation (unit:℃)

    (a) linear interpolation, (b) traditional cubic spline interpolation

    DownLoad: Full-Size Img PowerPoint

    Fig. 6  The zonal mean vertical distribution of water vapor root mean square error after interpolation (unit:g·kg-1)

    (a) linear interpolation, (b) monotonic cubilc spline interpolation

    DownLoad: Full-Size Img PowerPoint

    Fig. 7  Average deviation of 72-hour temperature forecast from bulk test against FNL data during 7—20 July 2009(unit:K, shaded for deviation, gray line for forecast field)

    (a) linear interpolation of 700 hPa deviation field, (b) traditional cubic spline interpolation of 700 hPa deviation field, (c) zonal mean deviation vertical profile of linear interpolation, (d) zonal mean deviation vertical profile of traditional cubic spline interpolation

    DownLoad: Full-Size Img PowerPoint

    Fig. 8  Average deviation of 72-hour humidity forecast from bulk test against FNL data during 7—20 July 2009 (unit:g·kg-1, shaded for deviation, gray line for forecast field)

    (a) linear interpolation of 700 hPa deviation field, (b) monotonic cubic spline interpolation of 700 hPa deviation field, (c) zonal mean deviation vertical profile of linear interpolation, (d) zonal mean deviation vertical profile of monotonic cubic spline interpolation

    DownLoad: Full-Size Img PowerPoint

    Fig. 9  The hight anomaly correlation coefficient and root mean square error in the Northern Hemisphere during 7—20 July 2009

    DownLoad: Full-Size Img PowerPoint
  • [1]
    陈德辉, 沈学顺.新一代数值预报系统GRAPES研究进展.应用气象学报, 2006, 17(6): 773-777. doi:  10.11898/1001-7313.20060614
    [2]
    伍湘君, 金之雁, 黄丽萍, 等.GRAPES模式软件框架与实现.应用气象学报, 2005, 16(4): 539-546. doi:  10.11898/1001-7313.20050415
    [3]
    黄丽萍, 伍湘君, 金之雁.GRAPES模式标准初始化方案设计与实现.应用气象学报, 2005, 16(3): 374-384. doi:  10.11898/1001-7313.20050312
    [4]
    李兴良, 陈德辉, 沈学顺.不同垂直坐标系对垂直速度计算的影响.热带气象学报, 2005, 21(3): 265-276. http://www.cnki.com.cn/Article/CJFDTOTAL-RDQX200503005.htm
    [5]
    陈德辉, 杨学胜, 张红亮, 等.多尺度非静力通用模式框架的设计策略.应用气象学报, 2003, 14(4): 452-461. http://qikan.camscma.cn/jams/ch/reader/view_abstract.aspx?file_no=20030456&flag=1
    [6]
    胡江林, 沈学顺, 张红亮, 等.GRAPES模式动力框架的长期积分特征.应用气象学报, 2007, 18(3): 276-284. http://qikan.camscma.cn/jams/ch/reader/view_abstract.aspx?file_no=20070349&flag=1
    [7]
    陈德辉, 薛纪善, 杨学胜, 等.GRAPES新一代全球/区域多尺度统一数值预报模式总体设计研究.科学通报, 2008, 20:2396-2407. doi:  10.3321/j.issn:0023-074X.2008.20.002
    [8]
    徐国强, 陈德辉, 薛纪善, 等.GRAPES物理过程的优化试验及程序结构设计.科学通报, 2008, 53(20): 2428-2434. doi:  10.3321/j.issn:0023-074X.2008.20.006
    [9]
    薛海乐, 沈学顺, 苏勇.地形湍流拖曳力参数化及在GRAPES中的应用.应用气象学报, 2011, 22(2): 169-181. doi:  10.11898/1001-7313.20110206
    [10]
    Davies T, Cullen M J P, Malcolm A J, et al.A new dynamical core for the Met Office's global and regional modeling of the atmosphere.Q J R Meteorol Soc, 2005, 131:1759-1782. doi:  10.1256/qj.04.101
    [11]
    Cote J S, Methot A, Patoine A, et al.The operational CMC-MRB Global Environmental Multiscale (GEM) model.Part Ⅰ:Desgin considerations and formulation.Mon Wea Rev, 1998, 126:1373-1395. doi:  10.1175/1520-0493(1998)126<1373:TOCMGE>2.0.CO;2
    [12]
    Seity Y, Brousseau P, Malardel S, et al.The AROME-France Convective-Scale Operation Model.Mon Mea Rev, 2011, 139:976-991. https://www.researchgate.net/publication/258537490_The_AROME-France_convective-scale_operational_model
    [13]
    宇如聪, 徐幼平.AREM暴雨数值预报模式.北京:中国科学院大气物理研究所, 2004.
    [14]
    Skamarock W C, Klemp J B, Dudhia J, et al.A Descriotion of the Advanced Research WRF Version 3.NCAR Technical Note, 2008. http://www.oalib.com/references/13162859
    [15]
    IFS Documentation-Cy37r2 Operational Implementation 18 May 2011, Part Ⅲ:Dynamics and Numerical Procedures.ECMWF, 2011.
    [16]
    Staniforth A, White A, Wood N, et al.Unified Model Documentation Paper No.15, Joy of U.M.6.3-Model Formulation.UKMO, 2006.
    [17]
    苏勇, 沈学顺.相变修正方案在GRAPES模式标量平流中的应用.气象学报, 2009, 67(7): 1089-1100. http://www.cnki.com.cn/Article/CJFDTOTAL-QXXB200906020.htm
    [18]
    严家明, 刘诗斌, 李辉.线性插值的误差计算方法研究.弹箭与制导学报, 2005, 25(4): 111-112. http://www.cnki.com.cn/Article/CJFDTOTAL-DJZD200504037.htm
    [19]
    阎洪.气候时空数据的样条插值与应用.地理与地理信息科学, 2003, 19(5): 27-31. http://www.cnki.com.cn/Article/CJFDTOTAL-DLGT200305007.htm
    [20]
    狄扬波.应用样条插值方法对湿度场进行客观分析的尝试.热带气象, 1985, 1(1): 74-79. http://www.cnki.com.cn/Article/CJFDTOTAL-RDQX198501008.htm
    [21]
    Kruger C J C.Constrained Cubic Spline Interpolation.2003.http://www.korf.co.uk/spline.pdf.
    • Relative Articles
    • Supplements(0)
    • Cited By
  • 加载中
  • -->

Catalog

    Figures(9)

    Get Citation
    PDF
    XML
    Article views (3322) PDF downloads(1241) Cited by()
    • Received : 2013-01-15
    • Accepted : 2013-12-20
    • Published : 2014-03-31
    Journal Online
    Contact

    © 2020 Journal of Applied Meteorological Science   京ICP备15006967号-1

    Supported by Beijing Renhe Information Technology Co. Ltd

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return