Gravity waves, scale asymptotics and the pseudo-incompressible equations

Multiple-scale asymptotics is used to analyse the Euler equations for the dynamical situation of a gravity wave (GW) near breaking level. A simple saturation argument in combination with linear theory is used to obtain the relevant dynamical scales. As a small expansion parameter, the ratio of the inverse of the vertical wavenumber and potential temperature and pressure scale heights is used, which we allow to be of the same order of magnitude here. It is shown that the resulting equation hierarchy is consistent with that obtained from the pseudo-incompressible equations, both for non-hydrostatic and hydrostatic GWs, while this is not the case for the anelastic equations unless the additional assumption of sufficiently weak stratification is adopted. To describe vertical propagation of wavepackets over several atmospheric-scale heights, Wentzel-Kramers-Brillouin (WKB) theory is used to show that the pseudo-incompressible flow divergence generates the same amplitude equation that also obtains from the full Euler equations. This gives a mathematical justification for the use of the pseudo-incompressible equations in the study of GW breaking in the atmosphere for arbitrary background stratification. The WKB theory interestingly even holds at wave amplitudes close to static instability. In the mean-flow equations, we obtain in addition to the classic wave-induced momentum-flux divergences a wave-induced correction of hydrostatic balance in the vertical momentum equation, which cannot be obtained from Boussinesq or anelastic dynamics.