A single-mode approximation for gravity waves in the thermosphere

Abstract In previous work (Knight et al., 2019) a multilayer Fourier method was introduced for obtaining solutions for a partial differential equation (PDE) associated with a dispersion relation for linear gravity waves in a vertically varying, anelastic, nonhydrostatic, viscous, and thermally diffusive atmosphere over some altitude range, given boundary conditions that allow waves to enter at the lower boundary and exit at the top. The previous work described an exact solution method for the dispersion-relation PDE using all wave modes. The current work describes a type of approximate solution, based on a single upgoing gravity-wave mode for each Fourier component, which is computationally cheaper and simpler to implement. Two different approaches to defining the single-mode solution are presented, one based on the Wentzel–Kramers–Brillouin (WKB) method and the other based on a simplified version of the multilayer approximation. The two approaches are shown to be mathematically equivalent. The accuracy of the single-mode approximation is investigated through comparison with the all-mode solution for a realistic example in which gravity waves generated by a convective storm in the troposphere propagate into the upper thermosphere and for which the single-mode approximation gives reasonably accurate results. It is necessary to include imaginary frequency shifts in some Fourier components in order for the different types of modes (e.g., upgoing vs. downgoing and gravity-wave vs. dissipative) to be well-defined in the presence of molecular viscosity. A systematic approach to determining the smallest possible imaginary frequency shift for any given Fourier component is introduced. The new method efficiently suppresses a specific type of numerical artifact associated with large imaginary frequency shifts.