Bifford P. Williams

Energy and momentum exchange among tides, gravity waves (GWs), and turbulence shapes the mesosphere and lower thermosphere, yet it is commonly treated as a one‐way cascade from large‐scale waves to small‐scale dissipation. Using coordinated sodium resonance and Rayleigh LiDARs at Poker Flat Research Range, Chatanika, Alaska, on 12 April 2022, we resolve three Kelvin–Helmholtz billow events near 87–90 km and simultaneously detect short‐period gravity‐wave packets (10–30 min; vertical wavelengths 10–20 km) propagating away from these instabilities. Wavelet analyses show that wave power increases following the peak of the instability, consistent with GWs radiated during billow breakdown. To test causality and quantify the coupling, we perform high‐resolution three‐dimensional compressible simulations constrained by the observed winds and temperatures; the model reproduces the observed billow evolution, the characteristic fishbone pattern in altitude–time, and the lagged growth of wave energy. Momentum‐ and heat‐flux divergences diagnose two linked pathways: tidal shear drives Kelvin–Helmholtz instability and turbulent mixing (downscale transfer), while the instability radiates GWs that export energy and momentum to larger scales (upscale transfer) and measurably modifies the semidiurnal tide. These results provide observationally constrained evidence for a wave–instability feedback loop in the middle atmosphere, motivating parameterizations that allow small‐scale instability and turbulence to generate, as well as dissipate, atmospheric waves.

The Polar Mesospheric Cloud (PMC) Turbulence experiment performed optical imaging and Rayleigh lidar PMC profiling during a 6‐day flight in July 2018. A mosaic of seven imagers provided sensitivity to spatial scales from ∼20 m to 100 km at a ∼2‐s cadence. Lidar backscatter measurements provided PMC brightness profiles and enabled definition of vertical displacements of larger‐scale gravity waves (GWs) and smaller‐scale instabilities of various types. These measurements captured an interval of strong, widespread Kelvin‐Helmholtz instabilities (KHI) occurring over northeastern Canada on July 12, 2018 during a period of significant GW activity. This paper addresses the evolution of the KHI field and the characteristics and roles of secondary instabilities within the KHI. Results include the imaging of secondary KHI in the middle atmosphere and multiple examples of KHI “tube and knot” (T&K) dynamics where two or more KH billows interact. Such dynamics have been identified clearly only once in the atmosphere previously. Results reveal that KHI T&K arise earlier and evolve more quickly than secondary instabilities of uniform KH billows. A companion paper by Fritts et al. (2022), https://doi.org/10.1029/2021JD035834 reveals that they also induce significantly larger energy dissipation rates than secondary instabilities of individual KH billows. The expected widespread occurrence of KHI T&K events may have important implications for enhanced turbulence and mixing influencing atmospheric structure and variability.

NRLMSIS® 2.0 is an empirical atmospheric model that extends from the ground to the exobase and describes the average observed behavior of temperature, eight species densities, and mass density via a parametric analytic formulation. The model inputs are location, day of year, time of day, solar activity, and geomagnetic activity. NRLMSIS 2.0 is a major, reformulated upgrade of the previous version, NRLMSISE‐00. The model now couples thermospheric species densities to the entire column, via an effective mass profile that transitions each species from the fully mixed region below ~70 km altitude to the diffusively separated region above ~200 km. Other changes include the extension of atomic oxygen down to 50 km and the use of geopotential height as the internal vertical coordinate. We assimilated extensive new lower and middle atmosphere temperature, O, and H data, along with global average thermospheric mass density derived from satellite orbits, and we validated the model against independent samples of these data. In the mesosphere and below, residual biases and standard deviations are considerably lower than NRLMSISE‐00. The new model is warmer in the upper troposphere and cooler in the stratosphere and mesosphere. In the thermosphere, N₂ and O densities are lower in NRLMSIS 2.0; otherwise, the NRLMSISE‐00 thermosphere is largely retained. Future advances in thermospheric specification will likely require new in situ mass spectrometer measurements, new techniques for species density measurement between 100 and 200 km, and the reconciliation of systematic biases among thermospheric temperature and composition data sets, including biases attributable to long‐term changes.

Polar mesospheric clouds (PMCs) occur in the summer near 82 ‐85km altitude due to seasonal changes of temperature and humidity. However, water vapor and associated PMCs have also been observed associated with rocket exhaust. The effects of this rocket exhaust on the temperature of the upper mesosphere are not well understood. To investigate these effects, 220 kg of pure water was explosively released at 85 km as part of the Super Soaker sounding rocket experiment on the night of January 25–26, 2018 at Poker Flat Research Range (65°N, 147°W). A cloud formed within 18 s and was measured by a ground‐based Rayleigh lidar. The peak altitude of the cloud appeared to descend from 92 to 78 km over 3 min. Temperatures leading up to the release were between 197 and 232 K, about 50 K above the summertime water frost point when PMCs typically occur. The apparent motion of the cloud is interpreted in terms of the expansion of the explosive release. Analysis using a water vapor radiative cooling code coupled to a microphysical model indicates that the cloud formed due to the combined effects of rapid radiative cooling (∼25 K) by meter‐scale filaments of nearly pure water vapor (∼1 ppv) and an increase in the frost point temperature (from 150 to 200 K) due to the high concentration of water vapor. These results indicate that water exhaust not only acts as a reservoir for mesospheric cloud production but also actively cools the mesosphere to induce cloud formation.

On the night of 18–19 October 2018, sodium resonance lidar measurements show the presence of overturning in the mesospheric sodium layer. Two independent tracers, sodium mixing ratio and potential temperature, derived from resonance and Rayleigh lidar measurements, reveal that vertical spreading of the sodium mixing ratio contours and a layer of convective instability coincide with this overturning. Analysis of lidar measurements also reveals the presence of gravity waves that propagate upward, are saturated, and dissipate at the height of the convective instability. The vertical spreading is analyzed in terms of turbulent diffusive transport using a model based on material continuity of sodium. Estimates of the turbulent eddy diffusion coefficient, K, and energy dissipation rate, ε are derived from the transport model. The energy dissipated by the gravity waves is also calculated and found to be sufficient to generate the turbulence. We consider three other examples of overturning, instability and spreading on the nights of: 17–18 February 2009, 25–26 January 2015, and 8–9 October 2018. For all four events we find that the values of K (∼1,000 m²/s) are larger and the values of ε (∼10–100 mW/kg) are of similar magnitude to those values typically reported by ionization gauge measurements. These examples also reveal that higher levels of turbulent mixing are consistently found in regions of lower stability.

Polar mesospheric cloud (PMC) imaging and lidar profiling performed aboard the 5.9‐day PMC Turbo balloon flight from Sweden to northern Canada in July 2018 revealed a wide variety of gravity wave (GW) and instability events occurring nearly continuously at approximately 82 km. We describe one event exhibiting GW breaking and associated vortex rings driven by apparent convective instability. Using PMC Turbo imaging with spatial and temporal resolution of 20 m and 2 s, respectively, we quantify the GW horizontal wavelength, propagation direction, and apparent phase speed. We identify vortex rings with diameters of 2–5 km and horizontal spacing comparable to their size. Lidar data show GW vertical displacements of ±0.3 km. From the data, we find a GW intrinsic frequency and vertical wavelength of 0.009 ± 0.003 rad s⁻¹ and 9 ± 4 km, respectively. We show that these values are consistent with the predictions of numerical simulations of idealized GW breaking. We estimate the momentum deposition rate per unit mass during this event to be 0.04 ± 0.02 m s⁻² and show that this value is consistent with the observed GW. Comparison to simulation gives a mean energy dissipation rate for this event of 0.05–0.4 W kg⁻¹, which is consistent with other reported in situ measurements at the Arctic summer mesopause.

The Polar Mesospheric Cloud Turbulence (PMC Turbo) instrument consists of a balloon‐borne platform which hosts seven cameras and a Rayleigh lidar. During a 6‐day flight in July 2018, the cameras captured images of Polar Mesospheric Clouds (PMCs) with a sensitivity to spatial scales from ~20 m to 100 km at a ~2‐s cadence and a full field of view (FOV) of hundreds of kilometers. We developed software optimized for imaging of PMCs, controlling multiple independent cameras, compressing and storing images, and for choosing telemetry communication channels. We give an overview of the PMC Turbo design focusing on the flight software and telemetry functions. We describe the performance of the system during its first flight in July 2018.

Though the Kelvin‐Helmholtz instability (KHI) has been extensively observed in the mesosphere, where breaking gravity waves produce the conditions required for instability, little has been done to describe quantitatively this phenomenon in detail in the mesopause and lower thermosphere, which are associated with the long‐lived shears at the base of this statically stable region. Using trimethylaluminum (TMA) released from two sounding rockets launched on 26 January 2018, from Poker Flat Research Range in Alaska, the KHI was observed in great detail above 100 km. Two sets of rocket measurements, made 30 min apart, show strong winds (predominantly meridional and up to 150 ms⁻¹) and large total shears (90 ms⁻¹ km⁻¹). The geomagnetic activity was low in the hours before the launches, confirming that the enhanced shears that triggered the KHI are not a result of the E‐region auroral jets. The four‐dimensional (three‐dimensional plus time) estimation of KHI billow features resulted in a wavelength, eddy diameter, and vertical length scale of 9.8, 5.2, and 3.8 km, respectively, centered at 102‐km altitude. The vertical and horizontal root‐mean‐square velocities measured 29.2 and 42.5 ms⁻¹, respectively. Although the wind structure persisted, the KHI structure changed significantly with time over the interval separating the two launches, being present only in the first launch. The rapid dispersal of the TMA cloud in the instability region was evidence of enhanced turbulent mixing. The analysis of the Reynolds and Froude numbers (Re = 7.2 × 10³ and Fr = 0.29, respectively) illustrates the presence of turbulence and weak stratification of the flow.

Two successive mesospheric bores were observed over northeastern Canada on 13 July 2018 in high‐resolution imaging and Rayleigh lidar profiling of polar mesospheric clouds (PMCs) performed aboard the PMC Turbo long‐duration balloon experiment. Four wide field‐of‐view cameras spanning an area of ~75 × 150 km at PMC altitudes captured the two evolutions occurring over ~2 hr and resolved bore and associated instability features as small as ~100 m. The Rayleigh lidar provided PMC backscatter profiling that revealed vertical displacements, evolving brightness distributions, evidence of instability character and depths, and insights into bore formation, ducting, and dissipation. Both bores exhibited variable structure along their phases, suggesting variable gravity wave (GW) source and bore propagation conditions. Both bores also exhibited small‐scale instability dynamics at their leading and trailing edges. Those at the leading edges comprised apparent Kelvin‐Helmholtz instabilities that were advected downward and rearward beneath the bore descending phases extending into an apparently intensified shear layer. Instabilities at the trailing edges exhibited alignments approximately orthogonal to the bore phases that resembled those seen to accompany GW breaking or intrusions arising in high‐resolution modeling of GW instability dynamics. Collectively, PMC Turbo bore imaging and lidar profiling enabled enhanced definition of bore dynamics relative to what has been possible by previous ground‐based observations, and a potential to guide new, three‐dimensional modeling of bore dynamics. The observed bore evolutions suggest potentially important roles for bores in the deposition of energy and momentum transported into the mesosphere and to higher altitudes by high‐frequency GWs achieving large amplitudes.

Polar mesospheric cloud (PMC) imaging and lidar profiling performed aboard the 5.9 day PMC Turbo balloon flight from Sweden to northern Canada in July 2018 revealed a wide variety of gravity wave (GW) and instability events occurring nearly continuously at approximately 82 km. We describe one event exhibiting GW breaking and associated vortex rings driven by apparent convective instability. Using PMC Turbo imaging with spatial and temporal resolution of 20 m and 2 s, respectively, we quantify the GW horizontal wavelength, propagation direction, and apparent phase speed, and we identify vortex rings with diameters of 3-5 km and horizontal spacing of ~5 km. Lidar data show GW vertical displacements of ±0.3 km. From the data, we find a GW intrinsic frequency and vertical wavelength of 0.009 ± 0.003 rad s-1 and 9 ± 4 km, respectively. We show that these values are consistent with the predictions of numerical simulations of idealized GW breaking. We estimate the momentum deposition rate per unit mass during this event to be 0.04 ± 0.02 m s-2 and show that this value is consistent with the observed GW. Comparison to simulation gives a mean energy dissipation rate for this event of 0.05-0.4 W kg-1, which is consistent with other reported in-situ measurements at the Arctic summer mesopause.