There are few observational techniques for measuring the distribution of kinetic energy within the mesosphere with a wide range of spatial and temporal scales. This study describes a method for estimating the three‐dimensional mesospheric wind field correlation function from specular meteor trail echoes. Each radar echo provides a measurement of a one‐dimensional projection of the wind velocity vector at a randomly sampled point in space and time. The method relies on using pairs of such measurements to estimate the correlation function of the wind with different spatial and temporal lags. The method is demonstrated using a multistatic meteor radar data set that includes ≈10⁵ meteor echoes observed during a 24‐hr time period. The new method is found to be in good agreement with the well‐established technique for estimating horizontal mean winds. High‐resolution correlation functions with temporal, horizontal, and vertical lags are also estimated from the data. The temporal correlation function is used to retrieve the kinetic energy spectrum, which includes the semidiurnal mode and a 3‐hr period wave. The horizontal and vertical correlation functions of the wind are then used to derive second‐order structure functions, which are found to be compatible with the Kolmogorov prediction for spectral distribution of kinetic energy in the turbulent inertial range. The presented method can be used to extend the capabilities of specular meteor radars. It is relatively flexible and has a multitude of applications beyond what has been shown in this study.

We describe in this study the analysis of small and large horizontal-scale gravity waves from datasets composed of images from multiple mesospheric airglow emissions as well as multistatic specular meteor radar (MSMR) winds collected in early November 2018, during the SIMONe–2018 (Spread-spectrum Interferometric Multi-static meteor radar Observing Network) campaign. These ground-based measurements are supported by temperature and neutral density profiles from TIMED/SABER (Thermosphere, Ionosphere, Mesosphere Energetics and Dynamics/Sounding of the Atmosphere using Broadband Emission Radiometry) satellite in orbits near Kühlungsborn, northern Germany (54.1^∘ N, 11.8^∘ E). The scientific goals here include the characterization of gravity waves and their interaction with the mean flow in the mesosphere and lower thermosphere and their relationship to dynamical conditions in the lower and upper atmosphere. We have obtained intrinsic parameters of small- and large-scale gravity waves and characterized their impact in the mesosphere via momentum flux (F_M) and momentum flux divergence (F_D) estimations. We have verified that a small percentage of the detected wave events is responsible for most of F_M measured during the campaign from oscillations seen in the airglow brightness and MSMR winds taken over 45 h during four nights of clear-sky observations. From the analysis of small-scale gravity waves (λ_h < 725 km) seen in airglow images, we have found F_M ranging from 0.04–24.74 m² s⁻² (1.62 ± 2.70 m² s⁻² on average). However, small-scale waves with F_M > 3 m² s⁻² (11 % of the events) transport 50 % of the total measured F_M. Likewise, wave events of F_M > 10 m² s⁻² (2 % of the events) transport 20 % of the total. The examination of large-scale waves (λ_h > 725 km) seen simultaneously in airglow keograms and MSMR winds revealed amplitudes > 35 %, which translates into F_M = 21.2–29.6 m² s⁻². In terms of gravity-wave–mean-flow interactions, these large F_M waves could cause decelerations of F_D = 22–41 m s⁻¹ d⁻¹ (small-scale waves) and F_D = 38–43 m s⁻¹ d⁻¹ (large-scale waves) if breaking or dissipating within short distances in the mesosphere and lower thermosphere region.

Measurements of kinetic energy in vortical and divergent fluctuations in the mesosphere and lower thermosphere can be used to study stratified turbulence (ST) and gravity waves. This can be done using horizontal correlation functions of the fluctuating component of velocity. This study introduces a novel method for estimating these correlation functions using radars that observe Doppler shifts of ionized specular meteor trails. The technique solves the correlation functions directly on a longitudinal‐transverse‐up coordinate system, assuming axial symmetry. This procedure is more efficient and leads to smaller uncertainties than a previous approach. The new technique is applied to a year‐long data set from a multistatic specular meteor radar network in Germany, to study the annual variability of kinetic energy within turbulent fluctuations at 87–93 km of altitude. In monthly averages, the kinetic energy is found to be nearly equipartitioned between vortical and divergent modes. Turbulent fluctuations maximize during the winter months with approximately 25% more energy in these months than at other times. The horizontal correlation functions are in agreement with the inertial subrange of ST, exhibiting a 2/3 power law in the horizontal lag direction, with an outermost scale of ST to be about 380 km. This suggests that horizontal correlation functions could be used to estimate turbulent energy transfer rates.

A virtual meteor radar system based on the upper‐atmosphere extension of the high‐resolution ICOsahedral Non‐hydrostatic general circulation model is constructed to validate multistatic specular meteor radar (SMR) analyses. The virtual radar system examines the validity of mean winds and gradients estimation techniques used in multistatic SMRs. The study is motivated by unexpected mean values and tide‐like features recently observed in the vertical velocities estimated from multistatic SMRs at different latitudes in the mesosphere and lower thermosphere. The proposed analysis confirms multistatic SMR systems' excellent capability to measure the horizontal mean wind components and gradient terms. It is also found that multistatic SMRs can estimate mean vertical winds if they have an amplitude greater than ±2 m/s. Due to the smoothing inherent to the model results, these results should be treated as lower bounds to the error incurred using real data. Hourly variability in vertical velocity estimates up to ±1–2 m/s in the observed vertical winds are due to contamination by small‐scale horizontal structures in the horizontal winds.

Specular meteor radars (SMRs) have significantly contributed to the understanding of wind dynamics in the mesosphere and lower thermosphere (MLT). We present a method to estimate horizontal correlations of vertical vorticity (Q_zz) and horizontal divergence (P) in the MLT, using line‐of‐sight multistatic SMRs velocities, that consists of three steps. First, we estimate 2D, zonal, and meridional correlation functions of wind fluctuations (with periods less than 4 hr and vertical wavelengths smaller than 4 km) using the wind field correlation function inversion (WCFI) technique. Then, the WCFI's statistical estimates are converted into longitudinal and transverse components. The conversion relation is obtained by considering the rotation about the vertical direction of two velocity vectors, from an east‐north‐up system to a meteor‐pair‐dependent cylindrical system. Finally, following a procedure previously applied in the upper troposphere and lower stratosphere to airborne wind measurements, the longitudinal and transverse spatial correlations are fitted, from which Q_zz, P, and their spectra are directly estimated. The method is applied to a special Spread spectrum Interferometric Multistatic meteor radar Observing Network data set, obtained over northern Germany for seven days in November 2018. The results show that in a quasi‐axisymmetric scenario, P was more than five times larger than Q_zz for the horizontal wavelengths range given by ∼50–400 km, indicating a predominance of internal gravity waves over vortical modes of motion as a possible explanation for the MLT mesoscale dynamics during this campaign.

In recent years, multistatic specular meteor radars (SMRs) have been introduced to study the Mesosphere and Lower Thermosphere (MLT) dynamics with increasing spatial and temporal resolution. SMRs, compared to other ground-based observations, have the advantage of continuously measuring the region between 80 and 100 km independent of weather, season, or time of day. In this paper, frequency spectra of MLT horizontal winds are explored through observations from a campaign using the SIMONe (Spread-spectrum Interferometric Multistatic meteor radar Observing Network) approach conducted in northern Germany in 2018 (hereafter SIMONe 2018). The 7-day SIMONe 2018 comprised of fourteen multistatic SMR links and allows us to build a substantial database of specular meteor trail events, collecting more than one hundred thousand detections per day within a geographic area of $$\sim$$∼ 500 km $$\times$$× 500 km. We have implemented two methods to obtain the frequency spectra of the horizontal wind components: (1) Mean Wind Estimation (MWE) and (2) Wind field Correlation Function Inversion (WCFI), which utilizes the mean and the covariances of the line of sight velocities, respectively. Monte Carlo simulations of a gravity wave spectral model were implemented to validate and compare both methods. The simulation analyses suggest that the WCFI helps us to capture the energy of smaller scale wind fluctuations than those capture with MWE. Characterization of the spectral slope of the horizontal wind at different MLT altitudes has been conducted on the SIMONe 2018, and it provides evidence that gravity waves with periods smaller than 7 h and greater than 2 h dominate with horizontal structures significantly larger than 500 km. In the future, these analyses can be extended to understand the significance of small-scale fluctuations in the MLT, which were not possible with conventional MWE methods. Graphical Abstract