The Operational Global Icosahedral-Hexagonal Gridpoint Model GME: Description and High-Resolution Tests

The German Weather Service (Deutscher Wetterdienst) has recently developed a new operational global numerical weather prediction model, named GME, based on an almost uniform icosahedral‐hexagonal grid. The GME gridpoint approach avoids the disadvantages of spectral techniques as well as the pole problem in latitude‐ longitude grids and provides a data structure extremely well suited to high efficiency on distributed memory parallel computers. The formulation of the discrete operators for this grid is described and evaluations that demonstrate their second-order accuracy are provided. These operators are derived for local basis functions that are orthogonal and conform perfectly to the spherical surface. The local basis functions, unique for each grid point, are the latitude and longitude of a spherical coordinate system whose equator and zero meridian intersect at the grid point. The prognostic equations for horizontal velocities, temperature, and surface pressure are solved using a semi-implicit Eulerian approach and for two moisture fields using a semi-Lagrangian scheme to ensure monotonicity and positivity. In the vertical direction, finite differences are applied in a hybrid (sigma pressure) coordinate system to all prognostic variables. The semi-implicit treatment of gravity waves presented here leads to a 3D Helmholtz equation that is diagonalized into a set of 2D Helmholtz equations that are solved by successive relaxation. Most of the same physical parameterizations used in the authors’ previous operational regional model, named EM, are employed in GME. Some results from the verification process for GME are provided and GME performance statistics on a Cray T3E1200 as well as on the ECMWF Fujitsu VPP5000 systems are summarized. For the case of the severe Christmas 1999 storm over France and Germany the pronounced sensitivity of the model with respect to the initial state is discussed. Finally, a test case is shown where it is currently possible, though not yet operationally practical, to run GME at 15-km resolution on the VPP5000.

[1]  G. R. Stuhne,et al.  New icosahedral grid-point discretization of the shallow water equation on the sphere , 1999 .

[2]  G. Mellor,et al.  A Hierarchy of Turbulence Closure Models for Planetary Boundary Layers. , 1974 .

[3]  A. Staniforth,et al.  Semi-Lagrangian integration schemes for atmospheric models - A review , 1991 .

[4]  David L. Williamson,et al.  Integration of the barotropic vorticity equation on a spherical geodesic grid , 1968 .

[5]  Peter Lynch,et al.  The Dolph-Chebyshev Window: A Simple Optimal Filter , 1997 .

[6]  P. Frederickson,et al.  Icosahedral Discretization of the Two-Sphere , 1985 .

[7]  Limin Yang,et al.  Development of a global land cover characteristics database and IGBP DISCover from 1 km AVHRR data , 2000 .

[8]  R. Heikes,et al.  Numerical Integration of the Shallow-Water Equations on a Twisted Icosahedral Grid , 1995 .

[9]  U. Schneider,et al.  Comparison of raingauge analyses, satellite-based precipitation estimates and forecast model results , 1996 .

[10]  F. Giraldo.  Lagrange-Galerkin Methods on Spherical Geodesic Grids , 1997 .

[11]  B. Ritter,et al.  A comprehensive radiation scheme for numerical weather prediction models with potential applications in climate simulations , 1992 .

[12]  G. Doms,et al.  The Nonhydrostatic Limited-Area Model LM (Lokal-Modell) of DWD: Part I: Scientific Documentation (Ve , 1999 .

[13]  André Robert,et al.  A stable numerical integration scheme for the primitive meteorological equations , 1981 .

[14]  David A. Randall,et al.  Numerical Integration of the Shallow-Water Equations on a Twisted Icosahedral Grid. Part II. A Detailed Description of the Grid and an Analysis of Numerical Accuracy , 1995 .

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

[16]  Woo-Sun Yang.  Variable Viscosity Thermal Convection at Infinite Prandtl Number in a Thick Spherical Shell , 1997 .

[17]  D. Burridge,et al.  A split semi‐implict reformulation of the Bushby‐Timpson 10‐level model , 1975 .

[18]  John R. Baumgardner,et al.  Three-dimensional treatment of convective flow in the earth's mantle , 1985 .

[19]  François Lott,et al.  A new subgrid‐scale orographic drag parametrization: Its formulation and testing , 1997 .

[20]  Akio Arakawa,et al.  Integration of the Nondivergent Barotropic Vorticity Equation with AN Icosahedral-Hexagonal Grid for the SPHERE1 , 1968 .

[21]  Todd D. Ringler,et al.  Modeling the Atmospheric General Circulation Using a Spherical Geodesic Grid: A New Class of Dynamical Cores , 2000 .

[22]  M. Tiedtke.  A Comprehensive Mass Flux Scheme for Cumulus Parameterization in Large-Scale Models , 1989 .

[23]  Fedor Mesinger.  Chapter 13 - Numerical Methods: The Arakawa Approach, Horizontal Grid, Global, and Limited-Area Modeling , 2000 .

[24]  J. Thuburn.  A PV-Based Shallow-Water Model on a Hexagonal-Icosahedral Grid , 1997 .

[25]  Yoshinobu Masuda,et al.  An Integration Scheme of the Primitive Equation Model with an Icosahedral-Hexagonal Grid System and its Application to the Shallow Water Equations , 1986 .

[26]  M. J. P. Cullen,et al.  Integrations of the primitive equations on a sphere using the finite element method , 1974 .

[27]  J. Louis.  A parametric model of vertical eddy fluxes in the atmosphere , 1979 .

[28]  G. R. Stuhne,et al.  Vortex Erosion and Amalgamation in a New Model of Large Scale Flow on the Sphere , 1996 .