Conservative split-explicit time integration methods for the compressible nonhydrostatic equations

Historically, time-split schemes for numerically integrating the nonhydrostatic compressible equations of motion have not formally conserved mass and other first-order flux quantities. In this paper, split-explicit integration techniques are developed that numerically conserve these properties by integrating prognostic equations for conserved quantities represented in flux form. These procedures are presented for both terrain-following height and hydrostatic pressure (mass) vertical coordinates, two potentially attractive frameworks for which the equation sets and integration techniques differ significantly. For each set of equations, the linear dispersion equation for acoustic/gravity waves is derived and analyzed to determine which terms must be solved in the small (acoustic) time steps and how these terms are represented in the time integration to achieve stability. Efficient techniques for including numerical filters for acoustic and external modes are also presented. Simulations for several idealized test cases in both the height and mass coordinates are presented to demonstrate that these integration techniques appear robust over a wide range of scales, from subcloud to synoptic.