A global hexagonal C‐grid non‐hydrostatic dynamical core (ICON‐IAP) designed for energetic consistency

This study describes a new global non-hydrostatic dynamical core (ICON-IAP: Icosahedral Nonhydrostatic model at the Institute for Atmospheric Physics) on a hexagonal C-grid which is designed to conserve mass and energy. Energy conservation is achieved by discretizing the antisymmetric Poisson bracket which mimics correct energy conversions between the different kinds of energy (kinetic, potential, internal). Because of the bracket structure this is even possible in a complicated numerical environment with (i) the occurrence of terrain-following coordinates with all the metric terms in it, (ii) the horizontal C-grid staggering on the Voronoi mesh and the complications induced by the need for an acceptable stationary geostrophic mode, and (iii) the necessity for avoiding Hollingsworth instability. The model is equipped with a Smagorinsky-type nonlinear horizontal diffusion. The associated dissipative heating is accounted for by the application of the discrete product rule for derivatives. The time integration scheme is explicit in the horizontal and implicit in the vertical. In order to ensure energy conservation, the Exner pressure has to be off-centred in the vertical velocity equation and extrapolated in the horizontal velocity equation. Test simulations are performed for small-scale and global-scale flows. A test simulation of linear non-hydrostatic flow over a rough mountain range shows the theoretically expected gravity wave propagation. The baroclinic wave test is extended to 40 days in order to check the Lorenz energy cycle. The model exhibits excellent energy conservation properties even in this strongly nonlinear and dissipative case. The Held–Suarez test confirms the reliability of the model over even longer time-scales. Copyright © 2012 Royal Meteorological Society

[1]  R. Sadourny.  The Dynamics of Finite-Difference Models of the Shallow-Water Equations , 1975 .

[2]  A. Hollingsworth,et al.  An internal symmetric computational instability , 1983 .

[3]  Valerio Lucarini,et al.  ENERGETICS OF CLIMATE MODELS: NET ENERGY BALANCE AND MERIDIONAL ENTHALPY TRANSPORT , 2009, 0911.5689.

[4]  A. Simmons,et al.  An Energy and Angular-Momentum Conserving Vertical Finite-Difference Scheme and Hybrid Vertical Coordinates , 1981 .

[5]  Luis Kornblueh,et al.  Sensitivity of Simulated Climate to Horizontal and Vertical Resolution in the ECHAM5 Atmosphere Model , 2006 .

[6]  R. Salmon,et al.  A general method for conserving quantities related to potential vorticity in numerical models , 2005 .

[7]  D. Williamson,et al.  A baroclinic instability test case for atmospheric model dynamical cores , 2006 .

[8]  William C. Skamarock,et al.  Numerical representation of geostrophic modes on arbitrarily structured C-grids , 2009, J. Comput. Phys..

[9]  W. Skamarock,et al.  The stability of time-split numerical methods for the hydrostatic and the nonhydrostatic elastic equations , 1992 .

[10]  D. Durran.  Numerical methods for wave equations in geophysical fluid dynamics , 1999 .

[11]  Oliver Fuhrer,et al.  Numerical consistency of metric terms in terrain-following coordinates , 2003 .

[12]  On the Hamiltonian formulation of the quasi‐hydrostatic equations , 1995 .

[13]  D. Williamson.  The Evolution of Dynamical Cores for Global Atmospheric Models(125th Anniversary Issue of the Meteorological Society of Japan) , 2007 .

[14]  Almut Gassmann.  Inspection of hexagonal and triangular C-grid discretizations of the shallow water equations , 2011, J. Comput. Phys..

[15]  D. Lüthi,et al.  A new terrain-following vertical coordinate formulation for atmospheric prediction models , 2002 .

[16]  William C. Skamarock,et al.  Conservative Transport Schemes for Spherical Geodesic Grids: High-Order Flux Operators for ODE-Based Time Integration , 2011 .

[17]  Shian‐Jiann Lin.  A “Vertically Lagrangian” Finite-Volume Dynamical Core for Global Models , 2004 .

[18]  Erich Becker,et al.  A Consistent Diffusion Dissipation Parameterization in the ECHAM Climate Model , 2006 .

[19]  Akio Arakawa,et al.  Chapter 1 - A Personal Perspective on the Early Years of General Circulation Modeling at UCLA , 2000 .

[20]  M. Sommer,et al.  Energy–Vorticity Theory of Ideal Fluid Mechanics , 2009 .

[21]  J. Thuburn,et al.  Numerical wave propagation on the hexagonal C-grid , 2008, J. Comput. Phys..

[22]  Hui Wan,et al.  Ensemble Held-Suarez test with a spectral transform model: Variability, sensitivity, and convergence , 2008 .

[23]  P. R. Bannon.  Hamiltonian Description of Idealized Binary Geophysical Fluids , 2003 .

[24]  Christiane Jablonowski,et al.  The pros and cons of diffusion, filters and fixers in Atmospheric General Circulation Models , 2011 .

[25]  Matthias Sommer,et al.  A conservative scheme for the shallow‐water system on a staggered geodesic grid based on a Nambu representation , 2009 .

[26]  A. Arakawa,et al.  A Potential Enstrophy and Energy Conserving Scheme for the Shallow Water Equations , 1981 .

[27]  William C. Skamarock,et al.  A unified approach to energy conservation and potential vorticity dynamics for arbitrarily-structured C-grids , 2010, J. Comput. Phys..

[28]  E. Becker,et al.  Nonlinear Horizontal Diffusion for GCMs , 2007 .

[29]  P. Morrison,et al.  Hamiltonian description of the ideal fluid , 1998 .

[30]  Almut Gassmann,et al.  Towards a consistent numerical compressible non‐hydrostatic model using generalized Hamiltonian tools , 2008 .

[31]  M. Suárez,et al.  A proposal for the intercomparison of the dynamical cores of atmospheric general circulation models , 1994 .

[32]  Edward N. Lorenz,et al.  The nature and theory of the general circulation of the atmosphere , 1967 .

[33]  E. Becker.  Frictional Heating in Global Climate Models , 2003 .

[34]  Hirofumi Tomita,et al.  An optimization of the Icosahedral grid modified by spring dynamics , 2002 .

[35]  A. Staniforth,et al.  A new dynamical core for the Met Office's global and regional modelling of the atmosphere , 2005 .

[36]  Jimy Dudhia,et al.  Conservative Split-Explicit Time Integration Methods for the Compressible Nonhydrostatic Equations , 2007 .

[37]  R. Salmon,et al.  A General Method for Conserving Energy and Potential Enstrophy in Shallow-Water Models , 2007 .

[38]  Lynne D. Talley,et al.  Generalizations of Arakawa's Jacobian , 1989 .