A strong mountain wave, observed over Central Europe on 12 January 2016, is simulated in 2D under two fixed background wind conditions representing opposite tidal phases. The aim of the simulation is to investigate the breaking of the mountain wave and subsequent generation of nonprimary waves in the upper atmosphere. The model results show that the mountain wave first breaks as it approaches a mesospheric critical level creating turbulence on horizontal scales of 8–30 km. These turbulence scales couple directly to horizontal secondary waves scales, but those scales are prevented from reaching the thermosphere by the tidal winds, which act like a filter. Initial secondary waves that can reach the thermosphere range from 60 to 120 km in horizontal scale and are influenced by the scales of the horizontal and vertical forcing associated with wave breaking at mountain wave zonal phase width, and horizontal wavelength scales. Large‐scale nonprimary waves dominate over the whole duration of the simulation with horizontal scales of 107–300 km and periods of 11–22 minutes. The thermosphere winds heavily influence the time‐averaged spatial distribution of wave forcing in the thermosphere, which peaks at 150 km altitude and occurs both westward and eastward of the source in the 2 UT background simulation and primarily eastward of the source in the 7 UT background simulation. The forcing amplitude is ∼2 × that of the primary mountain wave breaking and dissipation. This suggests that nonprimary waves play a significant role in gravity waves dynamics and improved understanding of the thermospheric winds is crucial to understanding their forcing distribution.

We use a nonlinear, fully compressible, two-dimensional numerical model to study the effects of dissipation on gravity wave packet spectra in the thermosphere. Numerical simulations are performed to excite gravity wave packets using either a time-dependent vertical body forcing at the bottom boundary or a specified initial wave perturbation. Three simulation case studies are performed to excite (1) a steady state monochromatic wave, (2) a spectrally broad wave packet, and (3) a quasi-monochromatic wave packet. In addition, we analyze (4) an initial condition simulation with an isothermal background. We find that, in cases where the wave is not continually forced, the dominant vertical wavelength decreases in time, predominantly due to a combination of refraction from the thermosphere and dissipation of the packets' high frequency components as they quickly reach high altitude. In the continually forced steady state case, the dominant vertical wavelength remains constant once initial transients have passed. The vertical wavelength in all simulations increases with altitude above the dissipation altitude (the point at which dissipation effects are greater than the wave amplitude growth caused by decreasing background density) at any fixed time. However, a shift to smaller vertical wavelength values in time is clearly exhibited as high-frequency, long vertical wavelength components reach high altitudes and dissipate first, to be replaced by slower waves of shorter vertical wavelength. Results suggest that the dispersion of a packet significantly determines its spectral evolution in the dissipative thermosphere.

Multiple events during the Deep Propagating Gravity Wave Experiment measurement program revealed mountain wave (MW) breaking at multiple altitudes over the Southern Island of New Zealand. These events were measured during several research flights from the National Science Foundation/National Center for Atmospheric Research Gulfstream V aircraft, utilizing a Rayleigh lidar, an Na lidar, and an Advanced Mesospheric Temperature Mapper simultaneously. A flight on 29 June 2014 observed MWs with horizontal wavelengths of ~80–120 km breaking in the stratosphere from ~10 to 50 km altitude. A flight on 13 July 2014 observed a horizontal wavelength of ~200–240 km MW extending from 20 to 90 km in altitude before breaking. Data from these flights show evidence for secondary gravity wave (SGW) generation near the breaking regions. The horizontal wavelengths of these SGWs are smaller than those of the breaking MWs, indicating a nonlinear generation mechanism. These observations reveal some of the complexities associated with MW breaking and the implications this can have on momentum fluxes accompanying SGWs over MW breaking regions.

A 2‐D nonlinear compressible model is used to simulate a large‐amplitude, multiscale mountain wave event over Mount Cook, NZ, observed as part of the Deep Propagating Gravity Wave Experiment (DEEPWAVE) campaign and to investigate its observable signatures in the hydroxyl (OH) layer. The campaign observed the presence of a λx=200 km mountain wave as part of the 22nd research flight with amplitudes of >20 K in the upper stratosphere that decayed rapidly at airglow heights. Advanced Mesospheric Temperature Mapper (AMTM) showed the presence of small‐scale (25–28 km) waves within the warm phase of the large mountain wave. The simulation results show rapid breaking above 70 km altitude, with the preferential formation of almost‐stationary vortical instabilities within the warm phase front of the mountain wave. An OH airglow model is used to identify the presence of small‐scale wave‐like structures generated in situ by the breaking of the mountain wave that are consistent with those seen in the observations. While it is easy to interpret these feature as waves in OH airglow data, a considerable fraction of the features are in fact instabilities and vortex structures. Simulations suggest that a combination of a large westward perturbation velocity and shear, in combination with strong perturbation temperature gradients, causes both dynamic and convective instability conditions to be met particularly where the wave wind is maximized and the temperature gradient is simultaneously minimized. This leads to the inevitable breaking and subsequent generation of smaller‐scale waves and instabilities which appear most prominent within the warm phase front of the mountain wave.

A combination of ray theory and 2-D time-dependent simulations is used to investigate the linear effects of a time-dependent, vertically, and horizontally inhomogeneous background horizontal wind field on the propagation, refraction, and reflection of small-scale gravity wave packets. Interactions between propagating waves of different scales are likely to be numerous and important. We find that a static medium-scale wave wind field of sufficient amplitude can channel and/or critical-level filter a small-scale wave or cause significant reflection, depending upon both waves' parameters. However, the inclusion of a time-dependent phase progression of the medium-scale wave can reduce energy loss through critical-level filtering by up to ∼70% and can also reduce reflection by up to ∼60% for the cases simulated. We also find that the relative direction of propagation between the small-scale and medium-scale wave can significantly affect small-scale wave filtering. When the phases are progressing in the same horizontal direction, the small-scale wave is far more likely to become trapped and ultimately critical-level filtered than if the phases are propagating in opposite horizontal directions unless reflection occurs first. Considerations of time-dependent winds associated with medium-scale-propagating waves and their directionality are important for assessing the propagation and dispersion of small-scale waves over large horizontal distances.

[1] Abundant short-period, small-scale gravity waves have been identified in the mesosphere and lower thermosphere over Halley, Antarctica, via ground-based airglow image data. Although many are observed as freely propagating at the heights of the airglow layers, new results under modeled conditions reveal that a significant fraction of these waves may be subject to reflections at altitudes above and below. The waves may at times be trapped within broad thermal ducts, spanning from the tropopause or stratopause to the base of the thermosphere ( 140 km), which may facilitate long-range propagation ( 1000s of km) under favorable wind conditions. Citation: Snively, J. B., K. Nielsen, M. P.Hickey, C. J. Heale, M. J. Taylor, and T. Moffat-Griffin (2013), Numerical and statistical evidence for long-range ducted gravity wave propagation over Halley, Antarctica, Geophys. Res. Lett., 40,[4813][4814][4815][4816][4817]

Gravity wave packets excited by a source of finite duration and size possess a broad frequency and wave number spectrum and thus span a range of temporal and spatial scales. Observing at a single location relatively close to the source, the wave components with higher frequency and larger vertical wavelength dominate at earlier times and at higher altitudes, while the lower frequency components, with shorter vertical wavelength, dominate during the latter part of the propagation. Utilizing observations from the Na lidar at Utah State University and the nearby Mesospheric Temperature Mapper at Bear Lake Observatory (41.9°N, 111.4°W), we investigate a unique case of vertical dispersion for a spectrally broad gravity wave packet in the mesopause region over Logan, Utah (41.7°N, 111.8°W), that occurred on 2 September 2011, to study the waves' evolution as it propagates upward. The lidar-observed temperature perturbation was dominated by close to a 1 h modulation at 100 km during the early hours but gradually evolved into a 1.5 h modulation during the second half of the night. The vertical wavelength also decreased simultaneously, while the vertical group and phase velocities of the packet apparently slowed, as it was approaching a critical level during the second half of the night. A two-dimensional numerical model is used to simulate the observed gravity wave processes, finding that the location of the lidar relative to the source can strongly influence which portion of the spectrum can be observed at a particular location relative to a source.

The investigation of atmospheric tsunamigenic acoustic and gravity wave (TAGW) dynamics, from the ocean surface to the thermosphere, is performed through the numerical computations of the 3D compressible nonlinear Navier-Stokes equations. Tsunami propagation is first simulated using a nonlinear shallow water model, which incorporates instantaneous or temporal evolutions of initial tsunami distributions (ITD). Ocean surface dynamics are then imposed as a boundary condition to excite TAGWs into the atmosphere from the ground level. We perform a case study of a large tsunami associated with the 2011 M9.1 Tohuku-Oki earthquake and parametric studies with simplified and demonstrative bathymetry and ITD. Our results demonstrate that TAGW propagation, controlled by the atmospheric state, can evolve nonlinearly and lead to wave self-acceleration effects and instabilities, followed by the excitation of secondary acoustic and gravity waves (SAGWs), spanning a broad frequency range. The variations of the ocean depth result in a change of tsunami characteristics and subsequent tilt of the TAGW packet, as the wave's intrinsic frequency spectrum is varied. In addition, focusing of tsunamis and their interactions with seamounts and islands may result in localized enhancements of TAGWs, which further indicates the crucial role of bathymetry variations. Along with SAGWs, leading long-period phases of the TAGW packet propagate ahead of the tsunami wavefront and thus can be observed prior to the tsunami arrival. Our modeling results suggest that TAGWs from large tsunamis can drive detectable and quantifiable perturbations in the upper atmosphere under a wide range of scenarios and uncover new challenges and opportunities for their observations.