AbstractThe mesosphere and lower thermosphere (MLT) is a fluid framework whose multiscale dynamics is determined by a superposition of non‐linear processes and by the interplay of gravity waves and turbulent motions. A thorough comprehension of this atmospheric region requires substantial observational infrastructure, needed to resolve and disentangle its complex dynamics. State‐of‐the‐art observational methods struggle to accurately capture mesoscale dynamics due to the inherent difficulty to perform observations at MLT altitudes. A majority of the observational methods rely on assumptions such as homogeneity, smoothness of the prognostic fields, or zero vertical wind velocities, which may not hold in the upper atmosphere at the mesoscales. In this study, we introduce a novel machine learning‐based approach HYPER (HYdrodynamic Point‐wise Environment Reconstructor), designed to characterize MLT dynamics. HYPER utilizes a physics‐informed neural network to project sparse Doppler meteor detections into four‐dimensional time‐series arrays containing the Cartesian components of the velocity field. This method combines meteor radar observations with the physics prescribed by the Navier‐Stokes equations. The validation of HYPER was conducted through a series of benchmarks on numerical data and the application of our algorithm on actual meteor radar observations, all of which yielded realistic approximations of the reconstructed physical fields. This innovative approach represents a significant step toward an accurate characterization of the MLT dynamics, overcoming the limitations of existing methods, and providing valuable insights into the behavior of this poorly accessible region of the atmosphere.

The NASA Ionospheric Connection Explorer (ICON) was launched in October 2019 and has been observing the upper atmosphere and ionosphere to understand the sources of their strong variability, to understand the energy and momentum transfer, and to determine how the solar wind and magnetospheric effects modify the internally-driven atmosphere-space system. The Far Ultraviolet Instrument (FUV) supports these goals by observing the ultraviolet airglow in day and night, determining the atmospheric and ionospheric composition and density distribution. Based on the combination of ground calibration and flight data, this paper describes how major instrument parameters have been verified or refined since launch, how science data are collected, and how the instrument has performed over the first 3 years of the science mission. It also provides a brief summary of science results obtained so far.

Atmospheric gravity waves play an important role in driving the dynamics of the Mesosphere and Lower Thermosphere and the basic structure of this region is determined by momentum deposition of these waves. Mesospheric bores are a type of non‐linear response that cause the amplification of gravity wave, due to trapping, that is characterized by a propagating step‐like jump followed by undulating waves. They require a stable layer or duct to travel horizontally with little attenuation thereby capable of transporting wave energy and momentum over larger distances. We present a prominent bright undular bore event observed in the mesospheric O( 1 S), O 2 , and OH emission layers on 16 March 2021 over Germany. A striking feature of this observation is the capture of bore's rapid dissipation around the center of the imager's field of view. The vertical temperature profile obtained from the satellite data indicates the presence of temperature inversion layer which acted as a thermal duct for the bore propagation. In addition, we have performed idealized two dimensional direct numerical simulations (DNS) of Navier‐Stokes equations under Boussinesq approximation. The DNS results reproduce many important characteristics of the observed airglow event like the nonlinear wave‐steepening, number of trailing waves, and its dissipation by implementing a thermal duct and a wave‐like perturbation. Furthermore, the DNS results also indicate that the duct width and amplitude of the initial perturbation have a considerable effect on the bore morphology.

The focus of this chapter is the signal processing in atmospheric radar imaging and in particular its application to study polar mesospheric summer echoes (PMSEs). Atmospheric received signals are stochastic and result from the scattering of non-homogeneous atmospheric irregularities. Atmospheric radar systems used in atmospheric radar imaging operate at very high frequencies (VHFs) due to the frequency-dependent scattering properties of the atmospheric irregularities. This results in huge antenna arrays, and it is crucial to make the best use of the still rather limited number of antennas. Applying standard ARI techniques in combination with multiple-input multiple-output (MIMO), we have been able to significantly improved the angular resolution of PMSE observations. This improvement has allowed us to study for the first time a Kelvin–Helmholtz instability event in four dimensions. We have also used MIMO on simple ARI application, i.e., a single scattering target on a given range and frequency, in combination with compressed sensing to develop a novel multistatic radar system called SIMONe to study the mesosphere and lower thermosphere dynamics with unprecedented temporal and spatial resolutions.

Abstract 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

The polar summer mesosphere is the Earth's coldest region, allowing the formation of mesospheric ice clouds. These ice clouds produce strong polar mesospheric summer echoes (PMSE) that are used as tracers of mesospheric dynamics. Here, we report the first observations of extreme vertical drafts (+/- 50 ms-1) in the mesosphere obtained from PMSE, characterized by velocities more than five standard deviations larger than the observed vertical wind variability. Using aperture synthesis radar imaging, the observed PMSE morphology resembles a solitary wave in a varicose mode, narrow along propagation (3-4 km) and elongated (>10 km) transverse to propagation direction, with a relatively large vertical extent (similar to 13 km). These spatial features are similar to previously observed mesospheric bores, but we observe only one crest with much larger vertical extent and higher vertical velocities.

A new incoherent scatter radar called EISCAT 3D is being constructed in northern Scandinavia. It will have the capability to produce volumetric images of ionospheric plasma parameters using aperture synthesis radar imaging. This study uses the current design of EISCAT 3D to explore the theoretical radar imaging performance when imaging electron density in the E region and compares numerical techniques that could be used in practice. Of all imaging algorithms surveyed, the singular value decomposition with regularization gave the best results and was also found to be the most computationally efficient. The estimated imaging performance indicates that the radar will be capable of detecting features down to approximately 90×90 m at a height of 100 km, which corresponds to a ≈0.05∘ angular resolution. The temporal resolution is dependent on the signal-to-noise ratio and range resolution. The signal-to-noise ratio calculations indicate that high-resolution imaging of auroral precipitation is feasible. For example, with a range resolution of 1500 m, a time resolution of 10 s, and an electron density of 2×1011m-3, the correlation function estimates for radar scatter from the E region can be measured with an uncertainty of 5 %. At a time resolution of 10 s and an image resolution of 90×90 m, the relative estimation error standard deviation of the image intensity is 10 %. Dividing the transmitting array into multiple independent transmitters to obtain a multiple-input–multiple-output (MIMO) interferometer system is also studied, and this technique is found to increase imaging performance through improved visibility coverage. Although this reduces the signal-to-noise ratio, MIMO has successfully been applied to image strong radar echoes as meteors and polar mesospheric summer echoes. Use of the MIMO technique for incoherent scatter radars (ISRs) should be investigated further.

Abstract This paper presents for the first time results on winds, tides, gradients of horizontal winds, and momentum fluxes at mesosphere and lower thermosphere altitudes over southern Patagonia, one of the most dynamically active regions in the world. For this purpose, measurements provided by SIMONe Argentina are investigated. SIMONe Argentina is a novel multistatic specular meteor radar system that implements a Spread‐spectrum Interferometric Multistatic meteor radar Observing Network (SIMONe) approach, and that has been operating since the end of September 2019. Average counts of more than 30,000 meteor detections per day result in tidal estimates with statistical uncertainties of less than 1 m/s. Thanks to the multistatic configuration, horizontal and vertical gradients of the horizontal winds are obtained, as well as vertical winds free from horizontal divergence contamination. The vertical gradients of both zonal and meridional winds exhibit strong tidal signatures. Mean momentum fluxes are estimated after removing the effects of mean winds using a 4‐h, 8‐km window in time and altitude, respectively. Reasonable statistical uncertainties of the momentum fluxes are obtained after applying a 28‐day averaging. Therefore, the momentum flux estimates presented in this paper represent monthly mean values of waves with periods of 4 h or less, vertical wavelengths shorter than 8 km, and horizontal scales less than 400 km.

This dataset includes meteor observations and wind estimates taken as part of the SIMONe 2018 campaign in northern Germany on November 5, 2018. The files are in netCDF-4 format and follow CF conventions (https://cfconventions.org/). We recommend loading the data using the xarray Python package. The SIMONe 2018 campaign ran from November 2, 2018 through November 9, 2018 in northern Germany. The radar network consisted of two pulsed transmitters in Juliusruh and Collm and a five-element interferometric MIMO-CW transmitter located in Kühlungsborn. Monostatic receiver stations co-located with the pulsed transmitters and six additional receiver stations located in Mechelsdorf, Breege, Neustrelitz, Guderup, Salzwedel, and Bornim were used to form a total of two monostatic and ten bistatic links. The data from these individual links were then processed to detect specular meteor echoes and estimate their parameters, including Doppler shift. The Doppler shifts, imposed by movement of the meteor trail due to the neutral winds, were then used to estimate the 4-D wind field. More details about the campaign can be found in Vierinen et al. (2019). Details for the wind field estimation can be found in Volz et al. (submitted).

The polar summer mesosphere is the Earth’s coldest region, allowing the formation of mesospheric ice clouds, potentially linked to climate change. These clouds produce strong radar echoes that are used as tracers of mesospheric dynamics. Here we report the first observations of extreme vertical drafts in the mesosphere, characterized by velocities larger than 40 m/s, i.e., more than five standard deviations larger than the observed wind variability. The morphology seems to resemble mesospheric bores, however the scales observed are much larger. Powerful vertical drafts, intermittent in space and time, emerge also in direct numerical simulations of stratified flows, predicting non-Gaussian statistics of vertical velocities. This evidence suggests that mesospheric bores might result from the interplay of gravity waves and turbulent motions. Our extreme event is interpreted as a mesospheric "super-bore", impacting mesospheric mixing and ice-formation, and would potentially impact planning of sub-orbital flights, and the investigation of biological material in the near space.

We present a method for estimating incident photoelectrons' energy spectra as a function of altitude by combining global scale far‐ultraviolet (FUV) and radio‐occultation (RO) measurements. This characterization provides timely insights important for accurate interpretation of ionospheric parameters inferred from the recently launched Ionospheric Connection Explorer (ICON) observations. Quantification of photoelectron impact is enabled by the fact that conjugate photoelectrons (CPEs) directly affect FUV airglow emissions but not RO measurements. We demonstrate a technique for estimation of photoelectron fluxes and their spectra by combining coincident ICON and COSMIC2 measurements and show that a significant fraction of ICON‐FUV measurements is affected by CPEs during the winter solstice. A comparison of estimated photoelectron fluxes with measured photoelectron spectra is used to gain further insights into the estimation method and reveals consistent values within 10–60 eV.

Mesoscale dynamics in the mesosphere and lower thermosphere (MLT) region have been difficult to study from either ground- or satellite-based observations. For understanding of atmospheric coupling processes, important spatial scales at these altitudes range between tens and hundreds of kilometers in the horizontal plane. To date, this scale size is challenging observationally, so structures are usually parameterized in global circulation models. The advent of multistatic specular meteor radar networks allows exploration of MLT mesoscale dynamics on these scales using an increased number of detections and a diversity of viewing angles inherent to multistatic networks. In this work, we introduce a four-dimensional wind field inversion method that makes use of Gaussian process regression (GPR), which is a nonparametric and Bayesian approach. The method takes measured projected wind velocities and prior distributions of the wind velocity as a function of space and time, specified by the user or estimated from the data, and produces posterior distributions for the wind velocity. Computation of the predictive posterior distribution is performed on sampled points of interest and is not necessarily regularly sampled. The main benefits of the GPR method include this non-gridded sampling, the built-in statistical uncertainty estimates, and the ability to horizontally resolve winds on relatively small scales. The performance of the GPR implementation has been evaluated on Monte Carlo simulations with known distributions using the same spatial and temporal sampling as 1 d of real meteor measurements. Based on the simulation results we find that the GPR implementation is robust, providing wind fields that are statistically unbiased with statistical variances that depend on the geometry and are proportional to the prior velocity variances. A conservative and fast approach can be straightforwardly implemented by employing overestimated prior variances and distances, while a more robust but computationally intensive approach can be implemented by employing training and fitting of model hyperparameters. The latter GPR approach has been applied to a 24 h dataset and shown to compare well to previously used homogeneous and gradient methods. Small-scale features have reasonably low statistical uncertainties, implying geophysical wind field horizontal structures as low as 20–50 km. We suggest that this GPR approach forms a suitable method for MLT regional and weather studies.

Atmospheric structures due to gravity waves, turbulence, Kelvin Helmholtz instabilities, etc. in the mesosphere are being studied with a varying of ground-based and satellite-based instruments. At scales less than 100 km, they are mainly studied with airglow imagers, lidars, and radars. Typical radar observations have not been able to resolve spatial and temporal ambiguities due to the strength of radar echoes, the size of the system, and/or the nature of the atmospheric irregularities. In this work we observed spatially and temporally resolved structures of PMSE with unprecedented horizontal resolution, using the improved radar imaging accuracy of the Middle Atmosphere Alomar Radar System (MAARSY) with the aid of a multiple-input multiple output (MIMO) technique. The studies are performed in both the brightness of the mesospheric echoes and their Doppler velocities. The resolutions achieved are less than 1 km in the horizontal direction, less than 300m in altitude, and less than 1 minute in time, in an area of ~15km x 15km around 85km of altitude. We present a couple of wavelike monochromatic events, one drifting with the background neutral wind, and one propagating against the neutral wind. Horizontal wavelengths, periods, and vertical and temporal coverage of the events are described and discussed. A theory of stratified turbulence is employed in the present study. In particular, it is shown that the structure that propagates with the background wind is a large-scale turbulent KHI event. Some important turbulence characteristics, such as a turbulent dissipation rate, buoyancy Reynolds number, and Froude number, support our conclusion.

Abstract. In recent years, multistatic specular meteor radars (SMRs) have been introduced to study the Mesosphere and Lower Thermosphere (MLT) dynamics. In this paper, the statistics of mesoscale MLT power spectra 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 seven-day SIMONe 2018 comprised of fourteen multistatic SMR links and allows to build a substantial database of specular meteor trail events, collecting more than one hundred thousand detections per day within a geographic area of ~ 500 km x 500 km. The two methods we propose to obtain the power spectra in frequency range are (1) Wind field Correlation Function Inversion (WCFI), which utilizes two-point correlations of specular meteor observations, and (2) Mean Wind Estimation (MWE), which determines the MLT winds and gradients from specular meteor observations. Monte Carlo simulations of a gravity wave spectral model were implemented to validate and compare both methods. The simulation analyses suggest that the WCFI is the viable option among them to study the second-order statistics of the MLT winds that helps to capture the energy of small-scale wind fluctuations. Characterization of the spectral slope at different MLT altitudes has been conducted on the SIMONe 2018, and it provides evidence that gravity waves with periods smaller than seven hours and greater than two hours are dominated by waves with horizontal wavelength significantly larger than 500 km, which might be associated to secondary gravity waves. We believe that the presented methods can help us bridge the observational gap between large and small-scale mesospheric wind fluctuations and also improve the capabilities of SMRs.

The mesosphere and lower thermosphere (MLT) region is dominated globally by dynamics at various scales: planetary waves, tides, gravity waves, and stratified turbulence. The latter two can coexist and be significant at horizontal scales less than 500 km, scales that are difficult to measure. This study presents a recently deployed multistatic specular meteor radar system, SIMONe Peru, which can be used to observe these scales. The radars are positioned at and around the Jicamarca Radio Observatory, which is located at the magnetic equator. Besides presenting preliminary results of typically reported large‐scale features, like the dominant diurnal tide at low latitudes, we show results on selected days of spatially and temporally resolved winds obtained with two methods based on: (a) estimation of mean wind and their gradients (gradient method), and (b) an inverse theory with Tikhonov regularization (regularized wind field inversion method). The gradient method allows improved MLT vertical velocities and, for the first time, low‐latitude wind field parameters such as horizontal divergence and relative vorticity. The regularized wind field inversion method allows the estimation of spatial structure within the observed area and has the potential to outperform the gradient method, in particular when more detections are available or when fine adaptive tuning of the regularization factor is done. SIMONe Peru adds important information at low latitudes to currently scarce MLT continuous observing capabilities. Results contribute to studies of the MLT dynamics at different scales inherently connected to lower atmospheric forcing and E‐region dynamo related ionospheric variability.

We present and characterize in time and three spatial dimensions a Kelvin-Helmholtz Instability (KHI) event from polar mesospheric summer echoes (PMSE) observed with the Middle Atmosphere Alomar Radar System. We use a newly developed radar imaging mode, which observed PMSE intensity and line of sight velocity with high temporal and angular resolution. The identified KHI event occurs in a narrow layer of 2.4 km thickness centered at 85 km altitude, is elongated along north-south direction, presents separation between billows of similar to 8 km in the east-west direction, and its billow width is similar to 3 km. The accompanying vertical gradients of the horizontal wind are between 35 and 45 m/s/km and vertical velocities inside the billows are +/- 12 m/s. Based on the estimated Richardson (<0.25), horizontal Froude (similar to 0.8), and buoyancy Reynolds (similar to 2.5 x 10(4)) numbers, the observed event is a KHI that occurs under weak stratification and generates strong turbulence.

Earth and Space Science Open Archive This preprint has been submitted to and is under consideration at Earth and Space Science. ESSOAr is a venue for early communication or feedback before peer review. Data may be preliminary.Learn more about preprints preprintOpen AccessYou are viewing the latest version by default [v3]First studies of mesosphere and lower thermosphere dynamics using a multistatic specular meteor radar network over southern PatagoniaAuthorsJuan FedericoConteiDJorge L.ChauiDJuan MiguelUrcoiDRalphLatteckiDJuhaVierinenJacoboSalvadorSee all authors Juan Federico ConteiDCorresponding Author• Submitting AuthorLeibniz-Institute of Atmospheric Physics at the Rostock UniversityiDhttps://orcid.org/0000-0002-9247-702Xview email addressThe email was not providedcopy email addressJorge L. ChauiDLeibniz-Institute of Atmospheric Physics at the Rostock UniversityiDhttps://orcid.org/0000-0002-2364-8892view email addressThe email was not providedcopy email addressJuan Miguel UrcoiDLeibniz Institute of Atmospheric PhysicsiDhttps://orcid.org/0000-0002-8894-3294view email addressThe email was not providedcopy email addressRalph LatteckiDLeibniz-Institute of Atmospheric Physics (LG)iDhttps://orcid.org/0000-0002-0001-7473view email addressThe email was not providedcopy email addressJuha VierinenUniversity of Tromsoview email addressThe email was not providedcopy email addressJacobo SalvadorCEILAPview email addressThe email was not providedcopy email address

Since the 1950s, specular meteor radars (SMRs) have been used to study the mesosphere and lower thermosphere (MLT) dynamics. Atmospheric parameters derived from SMRs are highly dependent on the number of detected meteors and the accuracy of the meteors' locations. Recently, incoherent and coherent multiple-input-multiple-output (MIMO) radar approaches combined with waveform diversity have been proposed to increase the number of detected meteors and to improve time, altitude, and horizontal resolution of the estimated wind fields. The incoherent MIMO approach refers to the addition of new transmit sites (widely separated), whereas the coherent MIMO refers to the addition of new transmit antennas in the same site (closely separated). In both the cases, a different pseudorandom sequence is transmitted from each antenna element. Unfortunately, the addition of new transmit antennas with different code sequences degrades the performance of conventional signal recovery algorithms. This is a consequence of the cross-interference between the transmitted waveforms, making it worse as the number of transmitters increases. In this article, we propose a signal recovery approach based on compressed sensing, taking advantage of the sparse nature of specular meteor echoes. The approach allows the exact recovery of weak echoes even in interference environments. Besides the advantage of the proposed approach to recover the meteor signal, we discuss the optimal selection of the transmitted waveforms and the minimum code length required for exact recovery. Additionally, we propose a modification of the orthogonal matching pursuit algorithm used in sparse problems to make it applicable in real-time analysis of large data. The success of the proposed approach is corroborated using Monte Carlo simulations and real data from a multi-static spread spectrum meteor radar network installed in northern Germany.