Mesoscale Gravity Waves Research Articles

Radar Signatures and Surface Observations of Elevated Convection Associated with Damaging Surface Winds

Abstract Identifying radar signatures indicative of damaging surface winds produced by convection remains a challenge for operational meteorologists, especially within environments characterized by strong low-level static stability and convection for which inflow is presumably entirely above the planetary boundary layer. Numerical model simulations suggest the most prevalent method through which elevated convection generates damaging surface winds is via “up–down” trajectories, where a near-surface stable layer is dynamically lifted and then dropped with little to no connection to momentum associated with the elevated convection itself. Recently, a number of unique convective episodes during which damaging surface winds were produced by apparently elevated convection coincident with mesoscale gravity waves were identified and cataloged for study. A novel radar signature indicative of damaging surface winds produced by elevated convection is introduced through six representative cases. One case is then explored further via a high-resolution model simulation and related to the conceptual model of up–down trajectories. Understanding the processes responsible for, and radar signature indicative of, damaging surface winds produced by gravity wave coincident convection will help operational forecasters identify and ultimately warn for a previously underappreciated phenomenon that poses a threat to lives and property. Significance Statement We identified unique radar and observational signatures of thunderstorms that produce damaging surface winds through a recently discovered mechanism. The radar and observational signatures can be used to issue warnings to protect lives and property in situations where damaging winds were previously unexpected. Key observational signatures include associated increases in surface pressure, sustained wind, and wind gust magnitudes, as well as little to no change or an increase in surface temperature. In addition, base radar data exhibit a divergence signature, including in regions of little or no detectable precipitation. Additional study is needed to answer why some atmospheric environments are supportive of the unique damaging-wind-producing mechanism while others are not.

Effects of mesoscale gravity waves on sporadic E simulated by a one-dimensional dynamic model

Effects of mesoscale gravity waves on sporadic E simulated by a one-dimensional dynamic model

Daily and Subdaily Wind and Divergence Variations Observed by a Shipborne Doppler Radar off the Southwestern Coast of Sumatra

Abstract This study investigated the daily cycle of the wind and divergence fields observed off the southwestern coast of Sumatra during a field campaign of the Years of the Maritime Continent pilot study. An algorithm was developed to retrieve kinematic variables from the single-Doppler data collected aboard the Research Vessel Mirai from 24 November to 13 December 2015. The observed daily cycles of the wind and divergence fields consisted of diurnal, semidiurnal, and short-term variations. Diurnal wind variation was characterized by deep and three-dimensional circulation. There was an approximate phase locking of the semidiurnal variation to the diurnal variation, both in the wind and divergence fields. The short-term wind variation occurred at a time scale of ∼1–3 h, and this pattern was associated with density currents or mesoscale gravity waves. Up to 73% of the daily vertical motion variance can be attributed to the diurnal and semidiurnal vertical motion variations with comparable strengths. Concurrently, precipitation propagated offshore in phase with density currents and mesoscale gravity waves. Our results suggest that diurnal and semidiurnal wind variations dominate the daily evolution of precipitation, whereas density currents and mesoscale gravity waves control offshore propagation. Additionally, it appears that the daily precipitation cycle is modulated by multiple-time-scale wind variabilities of less than a day, which is also responsible for the development of strong nocturnal convection off the southwestern coast of Sumatra. Significance Statement To improve our understanding of the daily wind and divergence cycle off the southwestern coast of Sumatra, we examined wind data collected by a shipborne Doppler radar. The observed daily cycles of the wind and divergence fields consisted of diurnal and semidiurnal variations, as well as a 1–3-h variation associated with a density current or mesoscale gravity wave. Our results suggest that diurnal and semidiurnal wind variations dominate the daily evolution of precipitation, whereas density currents and mesoscale gravity waves control offshore propagation. Thus, we highlight the role of multiple-time-scale wind variabilities of less than a day in modulating the daily precipitation cycle off the southwestern coast of Sumatra.

Ozone–gravity wave interaction in the upper stratosphere/lower mesosphere

Abstract. The increase in amplitudes of upward propagating gravity waves (GWs) with height due to decreasing density is usually described by exponential growth. Recent measurements show some evidence that the upper stratospheric/lower mesospheric gravity wave potential energy density (GWPED) increases more strongly during the daytime than during the nighttime. This paper suggests that ozone–gravity wave interaction can principally produce such a phenomenon. The coupling between ozone-photochemistry and temperature is particularly strong in the upper stratosphere where the time–mean ozone mixing ratio decreases with height. Therefore, an initial ascent (or descent) of an air parcel must lead to an increase (or decrease) in ozone and in the heating rate compared to the environment, and, hence, to an amplification of the initial wave perturbation. Standard solutions of upward propagating GWs with linear ozone–temperature coupling are formulated, suggesting amplitude amplifications at a specific level during daytime of 5 % to 15 % for low-frequency GWs (periods ≥4 h), as a function of the intrinsic frequency which decreases if ozone–temperature coupling is included. Subsequently, the cumulative amplification during the upward level-by-level propagation leads to much stronger GW amplitudes at upper mesospheric altitudes, i.e., for single low-frequency GWs, up to a factor of 1.5 to 3 in the temperature perturbations and 3 to 9 in the GWPED increasing from summer low to polar latitudes. Consequently, the mean GWPED of a representative range of mesoscale GWs (horizontal wavelengths between 200 and 1100 km, vertical wavelengths between 3 and 9 km) is stronger by a factor of 1.7 to 3.4 (2 to 50 J kg−1, or 2 % to 50 % in relation to the observed order of 100 J kg−1, assuming initial GW perturbations of 1 to 2 K in the middle stratosphere). Conclusively, the identified process might be an important component in the middle atmospheric circulation, which has not been considered up to now.

Radar observations of winds, waves and tides in the mesosphere and lower thermosphere over South Georgia island (54° S, 36° W) and comparison with WACCM simulations

Abstract. The mesosphere and lower thermosphere (MLT) is a dynamic layer of the earth's atmosphere. This region marks the interface at which neutral atmosphere dynamics begin to influence the upper atmosphere and ionosphere. However, our understanding of this region and our ability to accurately simulate it in global circulation models (GCMs) is limited by a lack of observations, especially in remote locations. To this end, a meteor radar was deployed from 2016 to 2020 on the remote mountainous island of South Georgia (54∘ S, 36∘ W) in the Southern Ocean. In this study we use these new measurements to characterise the fundamental dynamics of the MLT above South Georgia including large-scale winds, solar tides, planetary waves (PWs), and mesoscale gravity waves (GWs). We first present an improved method for time–height localisation of radar wind measurements and characterise the large-scale MLT winds. We then determine the amplitudes and phases of the diurnal (24 h), semidiurnal (12 h), terdiurnal (8 h), and quardiurnal (6 h) solar tides at this latitude. We find very large amplitudes up to 30 m s−1 for the quasi 2 d PW in summer and, combining our measurements with the meteor SAAMER radar in Argentina, show that the dominant modes of the quasi 5, 10, and 16 d PWs are westward 1 and 2. We investigate and compare wind variance due to both large-scale “resolved” GWs and small-scale “sub-volume” GWs in the MLT and characterise their seasonal cycles. Last, we use our radar observations and satellite temperature observations from the Microwave Limb Sounder to test a climatological simulation of the Whole Atmosphere Community Climate Model (WACCM). We find that WACCM exhibits a summertime mesopause near 80 km altitude that is around 10 K warmer and 10 km lower in altitude than observed. Above 95 km altitude, summertime meridional winds in WACCM reverse to poleward, but this not observed in radar observations in this altitude range. More significantly, we find that wintertime zonal winds between 85 to 105 km altitude are eastward up to 40 m s−1 in radar observations, but in WACCM they are westward up to 20 m s−1. We propose that this large discrepancy may be linked to the impacts of secondary GWs (2GWs) on the residual circulation, which are not included in most global models, including WACCM. These radar measurements can therefore provide vital constraints that can guide the development of GCMs as they extend upwards into this important region of the atmosphere.

The Mesoscale Gravity Wave Response to the 2022 Tonga Volcanic Eruption: AIRS and MLS Satellite Observations and Source Backtracing

AbstractOn 15 January 2022, the Hunga Tonga–Hunga Ha'apai volcano erupted violently. This exceptional event excited a manifold of atmospheric waves. Here, we focus on the mesoscale part of the wave spectrum. About 8.5 hr after the eruption a strong atmospheric gravity wave (GW) was observed in the stratosphere by the satellite instruments Atmospheric Infrared Sounder (AIRS) and Microwave Limb Sounder (MLS) in the vicinity of Tonga. By ray‐tracing, we confirm the eruption as the source of this GW event. We determine the wave characteristics of the GW in terms of horizontal and vertical wavelengths and GW momentum flux. The strength of the GW is compared to the usual Southern Hemisphere flux values during this week. The event is comparable to the strongest convective events considering MLS, and exceptionally strong considering AIRS, which observes faster waves only.

Studies on the Competition Between Homogeneous and Heterogeneous Ice Nucleation in Cirrus Formation

AbstractCirrus ice crystals are produced heterogeneously on ice‐nucleating particles (INPs) and homogeneously in supercooled liquid solution droplets. They grow by uptake of water molecules from the ice‐supersaturated vapor. The precursor particles, characterized by disparate ice nucleation abilities and number concentrations, compete for available vapor during ice formation events. We investigate cirrus formation events systematically in different temperature and updraft regimes, and for different INP number concentrations and time‐independent nucleation efficiencies. We consider vertical air motion variability due to mesoscale gravity waves and effects of supersaturation‐dependent deposition coefficients for water molecules on ice surfaces. We analyze ice crystal properties to better understand the dynamics of competing nucleation processes. We study the reduction of ice crystal numbers produced by homogeneous freezing due to INPs in both, individual simulations assuming constant updraft speeds and in ensemble simulations based on a stochastic representation of vertical wind speed fluctuations. We simulate and interpret probability distributions of total nucleated ice crystal number concentrations, showing signatures of homogeneous and heterogeneous nucleation. At typically observed, mean updraft speeds (≈15 cm s−1) competing nucleation should occur frequently, even at rather low INP number concentrations (<10 L−1). INPs increase cirrus occurrence and may alter cirrus microphysical properties without entirely suppressing homogeneous freezing events. We suggest to improve ice growth models, especially for low cirrus temperatures (<220 K) and low ice supersaturation (<0.3).

Mesoscale Gravity Waves and Midlatitude Weather: A Tribute to Fuqing Zhang

Abstract Over the course of his career, Fuqing Zhang drew vital new insights into the dynamics of meteorologically significant mesoscale gravity waves (MGWs), including their generation by unbalanced jet streaks, their interaction with fronts and organized precipitation, and their importance in midlatitude weather and predictability. Zhang was the first to deeply examine “spontaneous balance adjustment”—the process by which MGWs are continuously emitted as baroclinic growth drives the upper-level flow out of balance. Through his pioneering numerical model investigation of the large-amplitude MGW event of 4 January 1994, he additionally demonstrated the critical role of MGW–moist convection interaction in wave amplification. Zhang’s curiosity-turned-passion in atmospheric science covered a vast range of topics and led to the birth of new branches of research in mesoscale meteorology and numerical weather prediction. Yet, it was his earliest studies into midlatitude MGWs and their significant impacts on hazardous weather that first inspired him. Such MGWs serve as the focus of this review, wherein we seek to pay tribute to his groundbreaking contributions, review our current understanding, and highlight critical open science issues. Chief among such issues is the nature of MGW amplification through feedback with moist convection, which continues to elude a complete understanding. The pressing nature of this subject is underscored by the continued failure of operational numerical forecast models to adequately predict most large-amplitude MGW events. Further research into such issues therefore presents a valuable opportunity to improve the understanding and forecasting of this high-impact weather phenomenon, and in turn, to preserve the spirit of Zhang’s dedication to this subject.

Poleward transport of African dust to the Iberian Peninsula organized by a barrier jet and hydraulic jumps: Observations and high-resolution simulation analyses

Most air quality stations in Spain exceeded the European Union's daily PM10 limit due to the February 2016 Saharan dust outbreak, which resulted from two successive dust storms in Northwest Africa. This study identifies the meso-β/meso-γ-scale dynamical processes responsible for developing these dust storms using observations and high-resolution Weather Research and Forecasting model coupled with Chemistry simulations. The results revealed that the first dust storm was associated with a strong barrier jet (BJ). The BJ formed on the southeastern foothills of the Saharan Atlas Mountains (SAM) when an easterly-northeasterly low-level Mediterranean flow was blocked by a stably stratified layer close to the SAM. The BJ intensified just after sunrise on 20 February and the associated near-surface peak winds organized the first dust storm. The second dust storm was linked to a mesoscale gravity wave (MGW) and hydraulic jumps. A long-lived westward propagating MGW was triggered by a downslope flow interacting with the stable layer near the northeastern edge of the Tinrhert Plateau in eastern Algeria. When this MGW crossed the Tademaït Plateau, hydraulic jumps formed on its lee side. The strong winds accompanying these hydraulic jumps formed the second dust storm on 21 February. The lifted dust extended over a depth of 2–3 km in the growing daytime boundary layer and was advected poleward by the southerly/southeasterly mid-tropospheric winds. Our results underline the importance of resolving terrain-induced mesoscale processes to understand dust storm dynamics, which are difficult to represent in coarse-resolution numerical models.

Diurnal Characteristics of Gravity Waves over the Tibetan Plateau in 2015 Summer Using 10-km Downscaled Simulations from WRF-EnKF Regional Reanalysis

Diurnal variations of gravity waves over the Tibetan Plateau (TP) in summer 2015 were investigated based on high-resolution downscaled simulations from WRF-EnKF (Weather Research and Forecasting model and an ensemble Kalman filter) regional reanalysis data with particular emphasis on wave source, wave momentum fluxes and wave energies. Strong diurnal precipitations, which mainly happen along the south slope of the TP, tend to excite upward-propagating gravity waves. The spatial and temporal distributions of the momentum fluxes of small-scale (10–200 km) and meso-scale (200–500 km) gravity waves agree well with the diurnal precipitation distributions. The power spectra of momentum fluxes also show that the small- and meso-scale atmospheric processes become important during the period of the strongest rainfall. Eastward momentum fluxes and northward momentum fluxes are dominant. Wave energies are described in terms of kinetic energy (KE), potential energy (PE) and vertical fluctuation energy (VE). The diurnal variation and spatial distribution of VE in the lower stratosphere correspond to the diurnal rainfall in the troposphere.

Variability of the Stratospheric Quasi-Biennial Oscillation and Its Wave Forcing Simulated in the Beijing Climate Center Atmospheric General Circulation Model

AbstractIt is well known that the stratospheric quasi-biennial oscillation (QBO) is forced by equatorial waves with different horizontal/vertical scales, including Kelvin waves, mixed Rossby–gravity (MRG) waves, inertial gravity waves (GWs), and mesoscale GWs, but the relative contribution of each wave is currently not very clear. Proper representation of these waves is critical to the simulation of the QBO in general circulation models (GCMs). In this study, the vertical resolution in the Beijing Climate Center Atmospheric General Circulation Model (BCC-AGCM) is increased to better represent large-scale waves, and a mesoscale GW parameterization scheme, which is coupled to the convective sources, is implemented to provide unresolved wave forcing of the QBO. Results show that BCC-AGCM can spontaneously generate the QBO with realistic periods, amplitudes, and asymmetric features between westerly and easterly phases. There are significant spatiotemporal variations of parameterized convective GWs, largely contributing to a great degree of variability in the simulated QBO. In the eastward wind shear of the QBO at 20 hPa, forcing provided by resolved waves is 0.1–0.2 m s−1 day−1 and forcing provided by parameterized GWs is ~0.15 m s−1 day−1. On the other hand, westward forcings by resolved waves and parameterized GWs are ~0.1 and 0.4–0.5 m s−1 day−1, respectively. It is inferred that the eastward forcing of the QBO is provided by both Kelvin waves and mesoscale convective GWs, whereas the westward forcing is largely provided by mesoscale GWs. MRG waves barely contribute to the formation of the QBO in the model.

Modulation of precipitation by a mesoscale gravity wave in northern Italy: Detection and analysis

AbstractA 6 h event of increased and gravity‐wave modulated stratiform rain with a period of approximately 1 h is described and two wave generation mechanisms (geostrophic adjustment and shear instability) are investigated. The predictions of a linear stability analysis of the flow (an unstable mode of about 100 km wavelength originating from the shear instability of the tropospheric jet and propagating in a duct encompassing almost the entire troposphere) are shown to agree very well with the measurements from a microbarometer network and wind‐profiler soundings.

A Stochastic Representation of Temperature Fluctuations Induced by Mesoscale Gravity Waves

AbstractUbiquitous mesoscale gravity waves generate high cooling rates important for cirrus formation. We make use of long‐duration balloon observations to devise a probabilistic model describing mesoscale temperature uctuations (MTF) away from strong wave sources. We define background conditions based on observed probability distributions of temperature and underlying vertical wind speed fluctuations. We show theoretically that MTF are subject to damping at a rate near the Coriolis frequency when the vertical wind speed fluctuations are autocorrelated over a fraction of a Brunt‐Väisälä period. We find that for background wave activity, a decrease in temperature of 1K translates into cooling rate standard deviations and mean updraft speeds of 4–8Kh−1 and ≈ 10–20 cms−1, respectively, depending on latitude and stratification. We introduce an effective Coriolis frequency to generate cooling rates in equatorial regions consistent with balloon data. Above ice saturation, MTF are large enough to affect ice crystal nucleation. Our results help constrain uncertainty in aerosol‐cirrus interactions, provide insights to better meet challenges in comparing measurement data with model simulations, and support the development of cutting‐edge ice cloud schemes in global models.

The Impact of Mesoscale Gravity Waves on Homogeneous Ice Nucleation in Cirrus Clouds

AbstractEffects of a spectrum of mesoscale gravity waves on homogeneous aerosol freezing in midlatitude cirrus are studied by means of parcel model simulations that are driven by random vertical wind speeds constrained by balloon measurements. Stochastic wave forcing with mean updraft speeds of 5–20 cm/s leads to substantial nucleated ice crystal number concentrations (ICNC) of 0.1–1 cm−3 in situations with slow large‐scale cooling, which by itself would generate fewer ice crystals. The stochastic nature of wave‐driven air parcel temperatures enhances ICNC even further, but the times required to reach freezing conditions unsupported by large‐scale cooling may vary widely. In the presence of wave forcing, ice crystals with low ICNC (<1–10 L−1) are also generated by homogeneous freezing, albeit only rarely. Comparisons with aircraft measurements suggest significant influences of heterogeneous ice‐nucleating particles and ice crystal sedimentation on ICNC, but quantifying their individual contributions remains elusive.

Damping of gravity waves by kinetic processes in Jupiter's thermosphere

We employ a linearized full-wave model to show that the rovibrational damping (hereafter, RV damping) of gravity waves is an important source of atmospheric heating, which needs to be taken into account when modeling the wave processes in the jovian thermospheres. We find that RV damping of internal gravity waves becomes strong when the wave frequencies approach the thermal collision frequency. RV damping, compared to the damping by molecular viscosity and thermal diffusivity, can effectively lower the damping altitudes for waves with large frequencies, and requires larger upward wave energy fluxes to produce prescribed wave amplitudes for these waves. Particularly, the energy deposited by RV damping, viscous and thermal dissipation of the two waves identified in measurements by the Atmospheric Structure Instrument (ASI) on the Galileo Probe can produce a maximum temperature about 40 K higher than those in previous studies. We further demonstrate RV damping effect on the mesoscale gravity wave observed during the New Horizons flyby on Jupiter. Under the influence of RV damping, this wave has the vertical wavelength and damping altitude comparable to those of wave 2 retrieved from the ASI observations. The energy dissipation of these three observed waves, if coexist, should provide sufficient heating to produce the observed mean temperature between 357 km and 700 km above 1 bar.

Banded Convective Activity Associated With Mesoscale Gravity Waves Over Southern China

AbstractBanded convective activity that occurred near the south coast of China on 30 January 2018 was investigated through convection‐allowing simulations using a nonhydrostatic mesoscale model. The simulations capture reasonably well the observed characteristics of this event. The convective bands are found to be closely related to an episode of mesoscale gravity waves propagating northeastward with a wave speed of around 12 m/s and a primary wavelength of about ~40–50 km. Further analyses and sensitivity experiments reveal that the environment provides a wave duct for these gravity waves, with a thick low‐level stable layer below 850 hPa capped by a low‐stability reflecting layer with a critical level. The strength and depth of the low‐level stable layer determine the intrinsic phase speed and wavelength of the ducted gravity waves. In the sensitivity tests that the stable layer depth is reduced, the wave characteristics change according to what are predicted with the wave duct theory. The convective bands collocate and propagate in phase with the peak updraft regions of the gravity waves, suggesting strong interactions of convection and gravity waves, in which the ducted gravity waves can trigger and modulate convection, while latent heating from convection enhances the waves. In essence, both wave ducting and wave‐convection interaction are jointly responsible for the banded convective activity.

On the Momentum Budget of the Quasi-Biennial Oscillation in the Whole Atmosphere Community Climate Model

Abstract This study documents the contribution of equatorial waves and mesoscale gravity waves to the momentum budget of the quasi-biennial oscillation (QBO) in a 110-level version of the Whole Atmosphere Community Climate Model. The model has high vertical resolution, 500 m, above the boundary layer and through the lower and middle stratosphere, decreasing gradually to about 1.5 km near the stratopause. Parameterized mesoscale gravity waves and resolved equatorial waves contribute comparable easterly and westerly accelerations near the equator. Westerly acceleration by resolved waves is due mainly to Kelvin waves of zonal wavenumber in the range k = 1–15 and is broadly distributed about the equator. Easterly acceleration near the equator is due mainly to Rossby–gravity (RG) waves with zonal wavenumbers in the range k = 4–12. These RG waves appear to be generated in situ during both the easterly and westerly phases of the QBO, wherever the meridional curvature of the equatorial westerly jet is large enough to produce reversals of the zonal-mean barotropic vorticity gradient, suggesting that they are excited by the instability of the jet. The RG waves produce a characteristic pattern of Eliassen–Palm flux divergence that includes strong easterly acceleration close to the equator and westerly acceleration farther from the equator, suggesting that the role of the RG waves is to redistribute zonal-mean vorticity such as to neutralize the instability of the westerly jet. Insofar as unstable RG waves might be present in the real atmosphere, mixing due to these waves could have important implications for transport in the tropical stratosphere.

Limited angle tomography of mesoscale gravity waves by the infrared limb-sounder GLORIA

Abstract. Three-dimensional measurements of gravity waves are required in order to quantify their direction-resolved momentum fluxes and obtain a better understanding of their propagation characteristics. Such 3-D measurements of gravity waves in the lowermost stratosphere have been provided by the airborne Gimballed Limb Observer for Radiance Imaging of the Atmosphere (GLORIA) using full angle tomography. Closed flight patterns of sufficient size are needed to acquire the full set of angular measurements for full angle tomography. These take about 2 h and are not feasible everywhere due to scientific reasons or air traffic control restrictions. Hence, this paper investigates the usability of limited angle tomography for gravity wave research based on synthetic observations. Limited angle tomography uses only a limited set of angles for tomographic reconstruction and can be applied to linear flight patterns. A synthetic end-to-end simulation has been performed to investigate the sensitivity of limited angle tomography to gravity waves with different wavelengths and orientations with respect to the flight path. For waves with wavefronts roughly perpendicular to the flight path, limited angle tomography and full angle tomography can derive wave parameters like wavelength, amplitude, and wave orientation with similar accuracy. For waves with a horizontal wavelength above 200 km and vertical wavelength above 3 km, the wavelengths can be retrieved with less than 10 % error, the amplitude with less than 20 % error, and the horizontal wave direction with an error below 10∘. This is confirmed by a comparison of results obtained from full angle tomography and limited angle tomography for real measurements taken on 25 January 2016 over Iceland. The reproduction quality of gravity wave parameters with limited angle tomography, however, depends strongly on the orientation of the waves with respect to the flight path. Thus, full angle tomography might be preferable in cases in which the orientation of the wave cannot be predicted or waves with different orientations exist in the same volume and thus the flight path cannot be adjusted accordingly. Also, for low-amplitude waves and short-scale waves full angle tomography has advantages due to its slightly higher resolution and accuracy.

Interactions between Mesoscale and Submesoscale Gravity Waves and Their Efficient Representation in Mesoscale-Resolving Models

AbstractAs present weather forecast codes and increasingly many atmospheric climate models resolve at least part of the mesoscale flow, and hence also internal gravity waves (GWs), it is natural to ask whether even in such configurations subgrid-scale GWs might impact the resolved flow and how their effect could be taken into account. This motivates a theoretical and numerical investigation of the interactions between unresolved submesoscale and resolved mesoscale GWs, using Boussinesq dynamics for simplicity. By scaling arguments, first a subset of submesoscale GWs that can indeed influence the dynamics of mesoscale GWs is identified. Therein, hydrostatic GWs with wavelengths corresponding to the largest unresolved scales of present-day limited-area weather forecast models are an interesting example. A large-amplitude WKB theory, allowing for a mesoscale unbalanced flow, is then formulated, based on multiscale asymptotic analysis utilizing a proper scale-separation parameter. Purely vertical propagation of submesoscale GWs is found to be most important, implying inter alia that the resolved flow is only affected by the vertical flux convergence of submesoscale horizontal momentum at leading order. In turn, submesoscale GWs are refracted by mesoscale vertical wind shear while conserving their wave-action density. An efficient numerical implementation of the theory uses a phase-space ray tracer, thus handling the frequent appearance of caustics. The WKB approach and its numerical implementation are validated successfully against submesoscale-resolving simulations of the resonant radiation of mesoscale inertia GWs by a horizontally as well as vertically confined submesoscale GW packet.

Observation of Kelvin–Helmholtz instabilities and gravity waves in the summer mesopause above Andenes in Northern Norway

Abstract. We present observations obtained with the Middle Atmosphere Alomar Radar System (MAARSY) to investigate short-period wave-like features using polar mesospheric summer echoes (PMSEs) as a tracer for the neutral dynamics. We conducted a multibeam experiment including 67 different beam directions during a 9-day campaign in June 2013. We identified two Kelvin–Helmholtz instability (KHI) events from the signal morphology of PMSE. The MAARSY observations are complemented by collocated meteor radar wind data to determine the mesoscale gravity wave activity and the vertical structure of the wind field above the PMSE. The KHIs occurred in a strong shear flow with Richardson numbers Ri < 0.25. In addition, we observed 15 wave-like events in our MAARSY multibeam observations applying a sophisticated decomposition of the radial velocity measurements using volume velocity processing. We retrieved the horizontal wavelength, intrinsic frequency, propagation direction, and phase speed from the horizontally resolved wind variability for 15 events. These events showed horizontal wavelengths between 20 and 40 km, vertical wavelengths between 5 and 10 km, and rather high intrinsic phase speeds between 45 and 85 m s−1 with intrinsic periods of 5–10 min.