Journal cover Journal topic
Geoscientific Model Development An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

Journal metrics

  • IF value: 4.252 IF 4.252
  • IF 5-year value: 4.890 IF 5-year 4.890
  • CiteScore value: 4.49 CiteScore 4.49
  • SNIP value: 1.539 SNIP 1.539
  • SJR value: 2.404 SJR 2.404
  • IPP value: 4.28 IPP 4.28
  • h5-index value: 40 h5-index 40
  • Scimago H index value: 51 Scimago H index 51
Volume 6, issue 3
Geosci. Model Dev., 6, 735-763, 2013
https://doi.org/10.5194/gmd-6-735-2013
© Author(s) 2013. This work is distributed under
the Creative Commons Attribution 3.0 License.
Geosci. Model Dev., 6, 735-763, 2013
https://doi.org/10.5194/gmd-6-735-2013
© Author(s) 2013. This work is distributed under
the Creative Commons Attribution 3.0 License.

Model description paper 05 Jun 2013

Model description paper | 05 Jun 2013

The ICON-1.2 hydrostatic atmospheric dynamical core on triangular grids – Part 1: Formulation and performance of the baseline version

H. Wan1,4,*, M. A. Giorgetta1, G. Zängl2, M. Restelli3, D. Majewski2, L. Bonaventura3, K. Fröhlich1,2, D. Reinert2, P. Rípodas2, L. Kornblueh1, and J. Förstner2 H. Wan et al.
  • 1Max Planck Institute for Meteorology, Hamburg, Germany
  • 2German Weather Service, Offenbach am Main, Germany
  • 3MOX – Department of Mathematics F. Brioschi, Politecnico di Milano, Milan, Italy
  • 4International Max Planck Research School on Earth System Modelling, Hamburg, Germany
  • *now at: Pacific Northwest National Laboratory, Richland, WA, USA

Abstract. As part of a broader effort to develop next-generation models for numerical weather prediction and climate applications, a hydrostatic atmospheric dynamical core is developed as an intermediate step to evaluate a finite-difference discretization of the primitive equations on spherical icosahedral grids. Based on the need for mass-conserving discretizations for multi-resolution modelling as well as scalability and efficiency on massively parallel computing architectures, the dynamical core is built on triangular C-grids using relatively small discretization stencils.

This paper presents the formulation and performance of the baseline version of the new dynamical core, focusing on properties of the numerical solutions in the setting of globally uniform resolution. Theoretical analysis reveals that the discrete divergence operator defined on a single triangular cell using the Gauss theorem is only first-order accurate, and introduces grid-scale noise to the discrete model. The noise can be suppressed by fourth-order hyper-diffusion of the horizontal wind field using a time-step and grid-size-dependent diffusion coefficient, at the expense of stronger damping than in the reference spectral model.

A series of idealized tests of different complexity are performed. In the deterministic baroclinic wave test, solutions from the new dynamical core show the expected sensitivity to horizontal resolution, and converge to the reference solution at R2B6 (35 km grid spacing). In a dry climate test, the dynamical core correctly reproduces key features of the meridional heat and momentum transport by baroclinic eddies. In the aqua-planet simulations at 140 km resolution, the new model is able to reproduce the same equatorial wave propagation characteristics as in the reference spectral model, including the sensitivity of such characteristics to the meridional sea surface temperature profile.

These results suggest that the triangular-C discretization provides a reasonable basis for further development. The main issues that need to be addressed are the grid-scale noise from the divergence operator which requires strong damping, and a phase error of the baroclinic wave at medium and low resolutions.

Publications Copernicus
Download
Citation
Share