A model of the three-dimensional rotating compressible Euler equations on the cubed sphere is presented. The model uses a mixed mimetic spectral element discretization which allows for the exact exchanges of kinetic, internal and potential energy via the compatibility properties of the chosen function spaces. A Strang carryover dimensional splitting procedure is used, with the horizontal dynamics solved explicitly and the vertical dynamics solved implicitly so as to avoid the CFL restriction of the vertical sound waves. The function spaces used to represent the horizontal dynamics are discontinuous across vertical element boundaries, such that each horizontal layer is solved independently so as to avoid the need to invert a global 3D mass matrix, while the function spaces used to represent the vertical dynamics are similarly discontinuous across horizontal element boundaries, allowing for the serial solution of the vertical dynamics independently for each horizontal element. The model is validated against standard test cases for baroclinic instability within an otherwise hydrostatically and geostrophically balanced atmosphere, and a non-hydrostatic gravity wave as driven by a temperature perturbation.

A mixed mimetic spectral element method is applied to solve the rotating shallow water equations. The mixed method uses the recently developed spectral element histopolation functions, which exactly satisfy the fundamental theorem of calculus with respect to the standard Lagrange basis functions in one dimension. These are used to construct tensor product solution spaces which satisfy the generalized Stokes theorem, as well as the annihilation of the gradient operator by the curl and the curl by the divergence. This allows for the exact conservation of first order moments (mass, vorticity), as well as higher moments (energy, potential enstrophy), subject to the truncation error of the time stepping scheme. The continuity equation is solved in the strong form, such that mass conservation holds point wise, while the momentum equation is solved in the weak form such that vorticity is globally conserved. While mass, vorticity and energy conservation hold for any quadrature rule, potential enstrophy conservation is dependent on exact spatial integration. The method possesses a weak form statement of geostrophic balance due to the compatible nature of the solution spaces and arbitrarily high order spatial error convergence.

In a previous article [J. Comp. Phys. 357 (2018) 282–304] [4], the mixed mimetic spectral element method was used to solve the rotating shallow water equations in an idealized geometry. Here the method is extended to a smoothly varying, non-affine, cubed sphere geometry. The differential operators are encoded topologically via incidence matrices due to the use of spectral element edge functions to construct tensor product solution spaces in H(rot), H(div) and L₂. These incidence matrices commute with respect to the metric terms in order to ensure that the mimetic properties are preserved independent of the geometry. This ensures conservation of mass, vorticity and energy for the rotating shallow water equations using inexact quadrature on the cubed sphere. The spectral convergence of errors are similarly preserved on the cubed sphere, with the generalized Piola transformation used to construct the metric terms for the physical field quantities.