We present the results of fully resolved direct numerical simulations of monochromatic gravity waves breaking in the middle atmosphere. The simulations are initialized with optimal perurbations of the gives waves. Given a wavelength of 3 km, the required grid sizes range up to 3.6 billion computational cells, depending on the necessary domain size and the turbulence intensity. Our results provide an insight into the mechanics of gravity wave breaking they will be of great value for the validation of lower order methods for the prediction of wave breaking.

To reduce the computational costs of numerical studies of gravity wave breaking in the atmosphere, the grid resolution has to be reduced as much as possible. Insufficient resolution of small-scale turbulence demands a proper turbulence parameterization in the framework of a large-eddy simulation (LES). The authors validate three different LES methods-the adaptive local deconvolution method (ALDM), the dynamic Smagorinsky method (DSM), and a naïve central discretization without turbulence parameterization (CDS4)-for three different cases of the breaking of well-defined monochromatic gravity waves. For ALDM, a modification of the numerical flux functions is developed that significantly improves the simulation results in the case of a temporarily very smooth velocity field. The test cases include an unstable and a stable inertia-gravity wave as well as an unstable high-frequency gravity wave. All simulations are carried out both in three-dimensional domains and in two-dimensional domains in which the velocity and vorticity fields are three-dimensional (so-called 2.5D simulations). The results obtained with ALDM and DSM are generally in good agreement with the reference direct numerical simulations as long as the resolution in the direction of the wave vector is sufficiently high. The resolution in the other directions has a weaker influence on the results. The simulations without turbulence parameterization are only successful if the resolution is high and the level of turbulence is comparatively low.

The differentially heated rotating annulus is a widely studied tabletop-size laboratory model of the general mid-latitude atmospheric circulation. The two most relevant factors of cyclogenesis, namely rotation and meridional temperature gradient are quite well captured in this simple arrangement. The radial temperature difference in the cylindrical tank and its rotation rate can be set so that the isothermal surfaces in the bulk tilt, leading to the formation of baroclinic waves. The signatures of these waves at the free water surface have been analyzed via infrared thermography in a wide range of rotation rates (keeping the radial temperature difference constant) and under different initial conditions. In parallel to the laboratory experiments, five groups of the MetStröm collaboration have conducted numerical simulations in the same parameter regime using different approaches and solvers, and applying different initial conditions and perturbations. The experimentally and numerically obtained baroclinic wave patterns have been evaluated and compared in terms of their dominant wave modes, spatio-temporal variance properties and drift rates. Thus certain "benchmarks" have been created that can later be used as test cases for atmospheric numerical model validation.

Durran's pseudo-incompressible equations are integrated in a mass and momentum conserving way with a new implicit turbulence model. This system is soundproof, which has two major advantages over fully compressible systems: the Courant-Friedrichs-Lewy (CFL) condition for stable time advancement is no longer dictated by the speed of sound and all waves in the model are clearly gravity waves (GW). Thus, the pseudo-incompressible equations are an ideal laboratory model for studying GW generation, propagation, and breaking. Gravity wave breaking creates turbulence that needs to be parameterized. For the first time the adaptive local deconvolution method (ALDM) for implicit large-eddy simulation (ILES) is applied to non- Boussinesq stratified flows. ALDM provides a turbulence model that is fully merged with the discretization of the flux function. In the context of non-Boussinesq stratified flows this poses some new numerical challenges-the solution of which is presented in this text. In numerical test cases the authors show the agreement of the results with the literature (Robert's hot-cold bubble test case), they present the sensitivity to the model's resolution and discretization, and they demonstrate qualitatively the behavior of the implicit turbulence model for a 2D breaking gravity wave packet.