The EISCAT (European Incoherent SCATer) Scientific Association has provided versatile incoherent scatter (IS) radar facilities on the mainland of northern Scandinavia (the EISCAT UHF and VHF radar systems) and on Svalbard (the electronically scanning radar ESR (EISCAT Svalbard Radar) for studies of the high-latitude ionised upper atmosphere (the ionosphere). The mainland radars were constructed about 30 years ago, based on technological solutions of that time. The science drivers of today, however, require a more flexible instrument, which allows measurements to be made from the troposphere to the topside ionosphere and gives the measured parameters in three dimensions, not just along a single radar beam. The possibility for continuous operation is also an essential feature. To facilitatefuture science work with a world-leading IS radar facility, planning of a new radar system started first with an EU-funded Design Study (2005–2009) and has continued with a follow-up EU FP7 EISCAT_3D Preparatory Phase project (2010–2014). The radar facility will be realised by using phased arrays, and a key aspect is the use of advanced software and data processing techniques. This type of software radar will act as a pathfinder for other facilities worldwide. The new radar facility will enable the EISCAT_3D science community to address new, significant science questions as well as to serve society, which is increasingly dependent on space-based technology and issues related to space weather. The location of the radar within the auroral oval and at the edge of the stratospheric polar vortex is also ideal for studies of the long-term variability in the atmosphere and global change. This paper is a summary of the EISCAT_3D science case, which was prepared as part of the EU-funded Preparatory Phase project for the new facility. Three science working groups, drawn from the EISCAT user community, participated in preparing this document. In addition to these working group members, who are listed as authors, thanks are due to many others in the EISCAT scientific community for useful contributions, discussions, and support.

Long-term and continuous observations of mesospheric–lower thermospheric winds are rare, but they are important to investigate climatological changes at these altitudes on timescales of several years, covering a solar cycle and longer. Such long time series are a natural heritage of the mesosphere–lower thermosphere climate, and they are valuable to compare climate models or long-term runs of general circulation models (GCMs). Here we present a climatological comparison of wind observations from six meteor radars at two conjugate latitudes to validate the corresponding mean winds and atmospheric diurnal and semidiurnal tides from three GCMs, namely the Ground-to-Topside Model of Atmosphere and Ionosphere for Aeronomy (GAIA), the Whole Atmosphere Community Climate Model Extension (Specified Dynamics) (WACCM-X(SD)), and the Upper Atmosphere ICOsahedral Non-hydrostatic (UA-ICON) model. Our results indicate that there are interhemispheric differences in the seasonal characteristics of the diurnal and semidiurnal tide. There are also some differences in the mean wind climatologies of the models and the observations. Our results indicate that GAIA shows reasonable agreement with the meteor radar observations during the winter season, whereas WACCM-X(SD) shows better agreement with the radars for the hemispheric zonal summer wind reversal, which is more consistent with the meteor radar observations. The free-running UA-ICON tends to show similar winds and tides compared to WACCM-X(SD).

Ground-based remote sensing of atmospheric parameters is often limited to single station observations by vertical profiles at a certain geographic location. This is a limiting factor for investigating gravity wave dynamics as the spatial information is often missing, e.g., horizontal wavelength, propagation direction or intrinsic frequency. In this study, we present a new retrieval algorithm for multistatic meteor radar networks to obtain tomographic 3-D wind fields within a pre-defined domain area. The algorithm is part of the Agile Software for Gravity wAve Regional Dynamics (ASGARD) and called 3D-Var, and based on the optimal estimation technique and Bayesian statistics. The performance of the 3D-Var retrieval is demonstrated using two meteor radar networks: the Nordic Meteor Radar Cluster and the Chilean Observation Network De Meteor Radars (CONDOR). The optimal estimation implementation provide statistically sound solutions and diagnostics from the averaging kernels and measurement response. We present initial scientific results such as body forces of breaking gravity waves leading to two counter-rotating vortices and horizontal wavelength spectra indicating a transition between the rotational k⁻³ and divergent $k^{-\mathrm{5}/\mathrm{3}}$ mode at scales of 80–120 km. In addition, we performed a keogram analysis over extended periods to reflect the latitudinal and temporal impact of a minor sudden stratospheric warming in December 2019. Finally, we demonstrate the applicability of the 3D-Var algorithm to perform large-scale retrievals to derive meteorological wind maps covering a latitude region from Svalbard, north of the European Arctic mainland, to central Norway.

Meteor radars have become widely used instruments to study atmospheric dynamics, particularly in the 70 to 110 km altitude region. These systems have been proven to provide reliable and continuous measurements of horizontal winds in the mesosphere and lower thermosphere. Recently, there have been many attempts to utilize specular and/or transverse scatter meteor measurements to estimate vertical winds and vertical wind variability. In this study we investigate potential biases in vertical wind estimation that are intrinsic to the meteor radar observation geometry and scattering mechanism, and we introduce a mathematical debiasing process to mitigate them. This process makes use of a spatiotemporal Laplace filter, which is based on a generalized Tikhonov regularization. Vertical winds obtained from this retrieval algorithm are compared to UA-ICON model data. This comparison reveals good agreement in the statistical moments of the vertical velocity distributions. Furthermore, we present the first observational indications of a forward scatter wind bias. It appears to be caused by the scattering center's apparent motion along the meteor trajectory when the meteoric plasma column is drifted by the wind. The hypothesis is tested by a radiant mapping of two meteor showers. Finally, we introduce a new retrieval algorithm providing a physically and mathematically sound solution to derive vertical winds and wind variability from multistatic meteor radar networks such as the Nordic Meteor Radar Cluster (NORDIC) and the Chilean Observation Network De meteOr Radars (CONDOR). The new retrieval is called 3DVAR+DIV and includes additional diagnostics such as the horizontal divergence and relative vorticity to ensure a physically consistent solution for all 3D winds in spatially resolved domains. Based on this new algorithm we obtained vertical velocities in the range of w = ± 1–2 m s⁻¹ for most of the analyzed data during 2 years of collection, which is consistent with the values reported from general circulation models (GCMs) for this timescale and spatial resolution.

Winds measured by lidar from the Aeolus satellite are compared with winds measured by two ground-based radars – MARA in Antarctica (70.77^∘ S, 11.73^∘ E) and ESRAD (67.88^∘ N, 21.10^∘ E) in Arctic Sweden – for the period 1 July–31 December 2019. Aeolus is a demonstrator mission to test whether winds measured by Doppler lidar from space can have sufficient accuracy to contribute to improved weather forecasting. A comprehensive programme of calibration and validation has been undertaken following the satellite launch in 2018, but, so far, direct comparison with independent measurements from the Arctic or Antarctic regions have not been made. The comparison covers heights from the low troposphere to just above the tropopause. Results for each radar site are presented separately for Rayleigh (clear) winds, Mie (cloudy) winds, sunlit (“summer”) and non-sunlit (“winter”) seasons, and ascending and descending satellite tracks. Horizontally projected line-of-sight (HLOS) winds from Aeolus, reprocessed using baseline 2B10, for passes within 100 km of the radar sites, are compared with HLOS winds calculated from 1 h averaged radar horizontal wind components. The agreement in most data subsets is very good, with no evidence of significant biases (<1 m s⁻¹). Possible biases are identified for two subsets (about −2 m s⁻¹ for the Rayleigh winds for the descending passes at MARA and about 2 m s⁻¹ for the Mie winds for the ascending passes at ESRAD, both in winter), but these are only marginally significant. A robust significant bias of about 7 m s⁻¹ is found for the Mie winds for the ascending tracks at MARA in summer. There is also some evidence for increased random error (by about 1 m s⁻¹) for the Aeolus Mie winds at MARA in summer compared to winter. This might be related to the presence of sunlight scatter over the whole of Antarctica as Aeolus transits across it during summer.

The Hunga Tonga–Hunga Ha′apai volcano eruption was a unique event that caused many atmospheric phenomena around the globe. In this study, we investigate the atmospheric gravity waves in the mesosphere/lower-thermosphere (MLT) launched by the volcanic explosion in the Pacific, leveraging multistatic meteor radar observations from the Chilean Observation Network De Meteor Radars (CONDOR) and the Nordic Meteor Radar Cluster in Fennoscandia. MLT winds are computed using a recently developed 3DVAR+DIV algorithm. We found eastward- and westward-traveling gravity waves in the CONDOR zonal and meridional wind measurements, which arrived 12 and 48 h after the eruption, and we found one in the Nordic Meteor Radar Cluster that arrived 27.5 h after the volcanic detonation. We obtained observed phase speeds for the eastward great circle path at both locations of about 250 m s⁻¹, and they were 170–150 m s⁻¹ for the opposite propagation direction. The intrinsic phase speed was estimated to be 200–212 m s⁻¹. Furthermore, we identified a potential lamb wave signature in the MLT winds using 5 min resolved 3DVAR+DIV retrievals.