This article describes the non‐hydrostatic dynamical core developed for the ICOsahedral Non‐hydrostatic (ICON) modelling framework. ICON is a joint project of the German Weather Service (DWD) and the Max Planck Institute for Meteorology (MPI‐M), targeting a unified modelling system for global numerical weather prediction (NWP) and climate modelling. Compared with the existing models at both institutions, the main achievements of ICON are exact local mass conservation, mass‐consistent tracer transport, a flexible grid nesting capability and the use of non‐hydrostatic equations on global domains. The dynamical core is formulated on an icosahedral‐triangular Arakawa C grid. Achieving mass conservation is facilitated by a flux‐form continuity equation with density as the prognostic variable. Time integration is performed with a two‐time‐level predictor–corrector scheme that is fully explicit, except for the terms describing vertical sound‐wave propagation. To achieve competitive computational efficiency, time splitting is applied between the dynamical core on the one hand and tracer advection, physics parametrizations and horizontal diffusion on the other hand. A sequence of tests with varying complexity indicates that the ICON dynamical core combines high numerical stability over steep mountain slopes with good accuracy and reasonably low diffusivity. Preliminary NWP test suites initialized with interpolated analysis data reveal that the ICON modelling system already achieves better skill scores than its predecessor at DWD, the operational hydrostatic Global Model Europe (GME), and at the same time requires significantly fewer computational resources.

Large-eddy simulations (LES) with the newThis is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited. R. Heinze et al.at building confidence in the model's ability to simulate small-to mesoscale variability in turbulence, clouds and precipitation. The results are encouraging: the high-resolution model matches the observed variability much better at small-to mesoscales than the coarser resolved reference model. In its highest grid resolution, the simulated turbulence profiles are realistic and column water vapour matches the observed temporal variability at short time-scales. Despite being somewhat too large and too frequent, small cumulus clouds are well represented in comparison with satellite data, as is the shape of the cloud size spectrum. Variability of cloud water matches the satellite observations much better in ICON than in the reference model. In this sense, it is concluded that the model is fit for the purpose of using its output for parametrization development, despite the potential to improve further some important aspects of processes that are also parametrized in the high-resolution model.

ICON‐A is the new icosahedral nonhydrostatic (ICON) atmospheric general circulation model in a configuration using the Max Planck Institute physics package, which originates from the ECHAM6 general circulation model, and has been adapted to account for the changed dynamical core framework. The coupling scheme between dynamics and physics employs a sequential updating by dynamics and physics, and a fixed sequence of the physical processes similar to ECHAM6. To allow a meaningful initial comparison between ICON‐A and the established ECHAM6‐LR model, a setup with similar, low resolution in terms of number of grid points and levels is chosen. The ICON‐A model is tuned on the base of the Atmospheric Model Intercomparison Project (AMIP) experiment aiming primarily at a well balanced top‐of atmosphere energy budget to make the model suitable for coupled climate and Earth system modeling. The tuning addresses first the moisture and cloud distribution to achieve the top‐of‐atmosphere energy balance, followed by the tuning of the parameterized dynamic drag aiming at reduced wind errors in the troposphere. The resulting version of ICON‐A has overall biases, which are comparable to those of ECHAM6. Problematic specific biases remain in the vertical distribution of clouds and in the stratospheric circulation, where the winter vortices are too weak. Biases in precipitable water and tropospheric temperature are, however, reduced compared to the ECHAM6. ICON‐A will serve as the basis of further development and as the atmosphere component to the coupled model, ICON‐Earth system model (ESM).

ICON (ICOsahedral Nonhydrostatic) is a unified modeling system for global numerical weather prediction (NWP) and climate studies. Validation of its dynamical core against a test suite for numerical weather forecasting has been recently published by Z€ angl et al. (2014). In the present work, an extension of ICON is presented that enables it to perform as a large eddy simulation (LES) model. The details of the implementation of the LES turbulence scheme in ICON are explained and test cases are performed to validate it against two standard LES models. Despite the limitations that ICON inherits from being a unified modeling system, it performs well in capturing the mean flow characteristics and the turbulent statistics of two simulated flow configurations-one being a dry convective boundary layer and the other a cumulustopped planetary boundary layer.

No abstract

Abstract. As part of a broader effort to develop nextgeneration models for numerical weather prediction and climate applications, a hydrostatic atmospheric dynamical core is developed as an intermediate step to evaluate a finitedifference 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.

We present a technique to synthesise telluric absorption and emission features both for in-situ wavelength calibration and for their removal from astronomical spectra. While the presented technique is applicable for a wide variety of optical and infrared spectra, we concentrate in this paper on selected high-resolution near-infrared spectra obtained with the CRIRES spectrograph to demonstrate its performance and limitation. We find that synthetic spectra reproduce telluric absorption features to about 2%, even close to saturated line cores. Thus, synthetic telluric spectra could be used to replace the observation of telluric standard stars, saving valuable observing time. This technique also provides a precise in-situ wavelength calibration, especially useful for high-resolution near-infrared spectra in the absence of other calibration sources.

Abstract. We present the first stage of a new online-coupled global to regional-scale modeling framework for the simulation of the spatiotemporal evolution of aerosols and trace gases. The underlying meteorological model is the new nonhydrostatic model system ICON (ICOsahedral Nonhydrostatic) which allows a local grid refinement with two-way interactions between the grids. We develop the extension ART (Aerosol and Reactive Trace gases) with the goal of simulating interactions between trace substances and the state of the atmosphere. Within this paper, we present the basic equations and give an overview of the physical parameterizations as well as numerical methods we use.First applications of the new model system for trace gases, monodisperse particles and polydisperse particles are shown. The simulated distribution of two very short-lived substances (VSLS), bromoform (CHBr 3 ) and dibromomethane (CH 2 Br 2 ) reflecting the fast upward transport shows a good agreement with observations and previous model studies. Also, the shape of the simulated tropical profiles is well reproduced. As an example for the treatment of monodisperse particles we present the simulated ash plume of the Eyjafjallajökull eruption in April 2010. Here, a novel approach for the source function is applied. The pattern of the simulated distribution of volcanic ash particles shows a good agreement with previous studies. As an example for the treatment of a polydisperse aerosol, where number densities and mass concentrations are accounted for, we simulated the annual emissions of sea salt. We obtain a total emission flux of 26.0 Pg yr −1 and an emission flux of particles with diameter less than 10 µm of 7.36 Pg yr −1 .