Middle Atmosphere Research Articles

Limb Temperature Observations in the Stratosphere and Mesosphere Derived from the OMPS Sensor

Molecular scattering (Rayleigh scattering) has been extensively used from the ground with lidars and from space to observe the limb, thereby deriving vertical temperature profiles between 30 and 80 km. In this study, we investigate how temperature can be measured using the new Ozone Mapping and Profiler Suite (OMPS) sensor, aboard the Suomi NPP and NOAA-21 satellites. The OMPS consists of three instruments whose main purpose is to study the composition of the stratosphere. One of these, the Limb Profiler (LP), measures the radiance of the limb of the middle atmosphere (stratosphere and mesosphere, 12 to 90 km altitude) at wavelengths from 290 to 1020 nm. This new data set has been used with a New Simplified Radiative Transfer Model (NSRTM) to derive temperature profiles with a vertical resolution of 1 km. To validate the method, the OMPS-derived temperature profiles were compared with data from four ground-based lidars and the ERA5 and MSIS models. The results show that OMPS and the lidars are in agreement within a range of about 5 K from 30 to 80 km. Comparisons with the models also show similar results, except for ERA5 beyond 50 km. We investigated various sources of bias, such as different attenuation sources, which can produce errors of up to 120 K in the UV range, instrumental errors around 0.8 K and noise problems of up to 150 K in the visible range for OMPS. This study also highlighted the interest in developing a new miniaturised instrument that could provide real-time observation of atmospheric vertical temperature profiles using a constellation of CubeSats with our NSRTM.

Signature of a zonally symmetric semidiurnal tide during major sudden stratospheric warmings and plausible mechanisms: a case study

Sudden Stratospheric Warming (SSW) is a winter phenomenon initiated primarily by the enhanced stationary planetary waves (SPWs), characterized by an increase in polar stratospheric temperature by a few tens of kelvin for a few days. Wave-wave non-linear interaction can produce secondary waves, with sum and difference frequencies of the primary wave frequencies. The sun-synchronous semidiurnal tide is a major component at mid and high latitude middle atmosphere, which non-linearly interacts with the dominant SPW in the stratosphere to produce the zonally symmetric semidiurnal tide component (S0), as observed during two boreal SSWs. The zonally symmetric distribution of ozone has also potential to excite the S0 component by absorption of solar ultraviolet radiation as evident during a rare Austral SSW. Overall, the present study sheds light on the dominant generation mechanisms involved in the S0 enhancement during the SSW.

On the рossibility of multichannel optical backscattering sondes for joint balloon and lidar studies of the aerosol composition of the middle atmosphere

Aerosol backscattering sondes in the practice of aerological sounding, along with lidar observations, are used at night to study and monitor polar stratospheric clouds, tropospheric and stratospheric aerosol, cirrus clouds, pyroconvection, volcanic aerosol, as well as to verify remote methods and means of ground-based and satellite-based aerosol observations. For aerosol sondes, a simple two-wave measurement technique is used, which makes it possible to diagnose changes in aerosol composition by color index. The possibilities of the two-wave technique have limitations, which are discussed in this article. Aerological sounding combined with lidar observations expands the wavelength range for multi-wavelength studies, and direct measurements of atmospheric temperature increase the accuracy of aerosol sensing. The paper considers the application of 3 or more wavelenght techniques. Data from probe measurements using wavelengths of 470, 528, 850 and 940 nm and lidar sensing at wavelengths of 355 and 532 nm are presented.

Studying infrasound propagation in the middle atmosphere with ICON and UA-ICON: comparison with the IFS and ground-based remote sensing

Infrasound signals are used to monitor various anthropogenic and natural sources. In particular, infrasound is one of the four verification technologies used by the International Monitoring System (IMS) of the Comprehensive Test-Ban Treaty organisation (CTBTO). To determine accurate source locations and energies, an accurate model of wind and temperature up to the lower thermosphere is necessary, hence operational NWP products are of great importance for routine infrasound monitoring activities. However, many of these models focus on tropospheric conditions, and the middle atmosphere (MA) is not well represented. UA-ICON is an upper atmosphere version of the ICON model that provides modelled atmospheric parameters up to 150 km. First, to assess ICON and IFS operational analysis products, comparisons to lidar observations are made. The main differences between both products were analysed with respect to winds and temperatures in the MA. Second, UA-ICON simulations outputs were compared to ICON and IFS fields to demonstrate the increased wave activity above ~30 km with UA-ICON. The added value of UA-ICON with respect to ICON and IFS for infrasound propagation simulations is discussed. The comparisons between the instrumental observations and the models will be presented, as well as comparisons between modelled and measured infrasound propagation.

Deep learning methods for modeling infrasound transmission loss in the middle atmosphere

Accurate modeling of infrasound transmission losses (TLs) is essential to assess the performance of the global International Monitoring System (IMS) infrasound network. Among existing propagation modeling tools, parabolic equation method (PE) enables TLs to be finely modeled, but its computational cost does not allow exploration of a large parameter space for operational monitoring applications. To reduce computation times, Brissaud et al. (2022) explored the potential of convolutional neural networks (CNNs) trained on a large set of regionally simulated wavefields (>1000 km from the source) to predict TLs with negligible computation times compared to PE simulations. However, this new method shows difficulties in upwind conditions, especially at low frequencies, and causal issues with winds at large distances from the source affecting ground TLs close to the source. In this study, we have developed an optimized CNN network designed to minimize prediction errors while predicting TLs from globally simulated combined temperature and wind fields spanning over propagation ranges of 4000 km. Our approach enhances the previously proposed one by implementing key optimizations that improve the overall architecture performances. The implemented model predicts TLs with an average error of 20 dB in the whole frequency band (0.1-4 Hz) and explored realistic atmospheric scenarios.

Internal gravity waves in the middle atmosphere and their impact on infrasound propagation across the International Monitoring System

Infrasound technology is used to monitor atmospheric events in order to guarantee compliance with the Comprehensive Nuclear-Test-Ban Treaty (CTBT). Over the 60 infrasound stations initially planned, 53 are operational to date as part of the International Monitoring System (IMS). As opposed to the waveform-related technologies that are used by the IMS to monitor the two other media of the geosphere (the underground and the oceans), infrasound waves are travelling in a continuously varying propagation medium. From minutes to hours and from a few hundreds of meters in the vertical to thousands of kilometers in the horizontal, atmospheric perturbations affect detection capabilities of the IMS at various scales. These perturbations need to be taken into account, for infrasound analysis to be able to provide accurate localization and characterization of sources of interest. In particular, internal gravity waves (atmospheric buoyancy waves) pose a serious challenge to atmospheric modelling because of the unresolved involved processes. They can dramatically affect transmission losses across the globe. Through numerical experiments and using high resolution modelling outputs of the atmosphere, backed up by atmospheric measurements, we demonstrate how gravity waves impact infrasound attenuation and drive signal amplification across the IMS.

Global-scale gravity wave analysis methodology for the ESA Earth Explorer 11 candidate CAIRT

Abstract. In the past, satellite climatologies of gravity waves (GWs) have initiated progress in their representation in global models. However, these could not provide the phase speed and direction distributions needed for a better understanding of the interaction between GWs and the large-scale winds directly. The ESA Earth Explorer 11 candidate CAIRT could provide such observations. CAIRT would use a limb-imaging Michelson interferometer resolving a wide spectral range, allowing temperature and trace gas mixing ratio measurements. With the proposed instrument design, a vertical resolution of 1 km, along-track sampling of 50 km, and across-track sampling of 25 km in a 400 km wide swath will be achieved. In particular, this allows for the observation of three-dimensional (3D), GW-resolving temperature fields throughout the middle atmosphere. In this work, we present the methodology for the GW analysis of CAIRT observations using a limited-volume 3D sinusoidal fit (S3D) wave analysis technique. We assess the capability of CAIRT to provide high-quality GW fields by the generation of synthetic satellite observations from high-resolution model data and comparison of the synthetic observations to the original model fields. For the assessment, wavelength spectra, phase speed spectra, horizontal distributions, and zonal means of GW momentum flux (GWMF) are considered. The atmospheric events we use to exemplify the capabilities of CAIRT are the 2006 sudden stratospheric warming (SSW) event, the quasi-biennial oscillation (QBO) in the tropics, and the mesospheric preconditioning phase of the 2019 SSW event. Our findings indicate that CAIRT would provide highly reliable observations not only of global-scale GW distributions and drag patterns but also of specific wave events and their associated wave parameters. Even under worse-than-expected noise levels of the instrument, the resulting GW measurements are highly consistent with the original model data. Furthermore, we demonstrate that the estimated GW parameters can be used for ray tracing, which physically extends the horizontal coverage of the observations beyond the orbit tracks.

Transient Tropopause Waves

Abstract Flight-level airborne observations have often detected gravity waves with horizontal wavelengths near the tropopause. Here, in situ and remote sensing aircraft data of these short gravity waves trapped along tropopause inversion layer and collected during a mountain-wave event over southern Scandinavia are analyzed to quantify their spectral energy and energy fluxes and to identify nonstationary modes. A series of three-dimensional numerical simulations are performed to explain the origin of these transient wave modes and to investigate the parameters on which they depend. It turns out that mountain-wave breaking in the middle atmosphere and the subsequent modification of the stratospheric flow are the key factors for the occurrence of trapped modes with . In particular, the intermittent and periodic breaking of mountain waves in the lower stratosphere forms a wave duct directly above the tropopause, in which the short gravity waves are trapped. The characteristics of the trapped, downstream-propagating waves are mainly controlled by the sharpness of the tropopause inversion layer. It could be demonstrated that different settings for optimizing the numerical solver have a significantly smaller influence on the solutions.

A review of global long-term changes in the mesosphere, thermosphere and ionosphere: A starting point for inclusion in (semi-) empirical models

The climate of the upper atmosphere, including the mesosphere, thermosphere and ionosphere, is changing. As data records are much more limited than in the lower atmosphere and solar variability becomes increasingly dominant at higher altitudes, accurate trend detection and attribution is not straightforward. Nonetheless, observations reliably indicate that, on average, the mesosphere has been cooling, the density in the thermosphere has been decreasing, and ionospheric layers have been shifting down. These global mean changes can be largely attributed to the increase in CO2 concentration, which causes cooling and thermal contraction in the middle and upper atmosphere. The decline in thermosphere density is particularly relevant from a practical viewpoint, as this reduces atmospheric drag and thereby increases orbital lifetimes and the build-up of space debris. Long-term changes in the ionosphere can have further practical implications and are not only driven by the increase in CO2 concentration, but also by changes in the Earth’s magnetic field. The empirical models that are mostly used to inform applications in industry on the state of the upper atmosphere, as well as being widely used in science, do not yet properly account for long-term trends in the mesosphere, thermosphere and ionosphere. This is problematic when long-term future projections are needed or models rely strongly on older data. This review provides an overview of the main evidence of long-term trends observed in the mesosphere, thermosphere and ionosphere, together with the latest insights on what causes these trends. It is hoped that this may serve as a starting point to include long-term trends in (semi-) empirical models to benefit all users of these models. We also offer some thoughts on how this could be approached.

Analysis of a secondary 16-day planetary wave generation through nonlinear interactions in the atmosphere

Using a nonlinear model of the general circulation of the middle and upper atmosphere (MUAM), spatio-temporal structures of planetary waves (PWs) during boreal winter were studied. Modeling of global atmospheric circulation was performed for January–February. Despite the tropospheric PW sources shaped in the model, the phenomenon of 16-day PW excitation arise out of internal atmospheric sources in the southern lower thermosphere was discovered. To explain the observed phenomenon, a number of numerical experiments were carried out according to different scenarios with a selective turning (on/off) tropospheric sources of PW individual modes (having periods of 4–16 days) in the model. Also, the evolution of perturbed potential enstrophy for a 16-day PW, as well as the contribution of nonlinear interactions between individual PW to it, was studied. This made it possible for the first time to demonstrate explicitly the process of generation a secondary 16-day PW as a result of the nonlinear interconnection of 4-day and 5-day PWs.Graphical

Dual-Wavelength Interferometric Detection Technology for Wind and Temperature Fields in the Martian Middle and Upper Atmosphere Based on LCTF

Dual-Wavelength Interferometric Detection Technology for Wind and Temperature Fields in the Martian Middle and Upper Atmosphere Based on LCTF

Thermal tides in the middle atmosphere at mid-latitudes measured with a ground-based microwave radiometer

Abstract. In recent decades, theoretical studies and numerical models of thermal tides have gained attention. It has been recognized that tides have a significant influence on the dynamics of the middle and upper atmosphere; as they grow in amplitude and propagate upward, they transport energy and momentum from the lower to the upper atmosphere, contributing to the vertical coupling between atmospheric layers. The superposition of tides with other atmospheric waves leads to non-linear wave–wave interactions. However, direct measurements of thermal tides in the middle atmosphere are challenging and are often limited to satellite measurements in the tropics and at low latitudes. Due to orbit geometry, such observations provide only a reduced insight into the short-term variability in atmospheric tides. In this paper, we present tidal analysis from 5 years of continuous observations of middle-atmospheric temperatures. The measurements were performed with the ground-based temperature radiometer TEMPERA (TEMPErature RAdiometer), which was developed at the University of Bern in 2013 and was located in Bern (46.95° N, 7.45° E) and Payerne (46.82° N, 6.94° E). TEMPERA achieves a temporal resolution of 1–3 h and covers the altitude range between 25–50 km. Using an adaptive spectral filter with a vertical regularization (ASF2D) for the tidal analysis, we found maximum amplitudes for the diurnal tide of approximately 2.4 K, accompanied by seasonal variability. The maximum amplitude was reached on average at an altitude of 43 km, which also reflected some seasonal characteristics. We demonstrate that TEMPERA is suitable for providing continuous temperature soundings in the stratosphere and lower mesosphere with a sufficient cadence to infer tidal amplitudes and phases for the dominating tidal modes. Furthermore, our measurements exhibit a dominating diurnal tide and smaller amplitudes for the semidiurnal and terdiurnal tides in the stratosphere.

Effect of microphysics scheme and data assimilation on hydrometeor and radiative flux simulations in the Arctic.

Although clouds are a major factor influencing atmospheric environments in the Arctic, numerical simulations of Arctic clouds are uncertain. In this study, the effects of microphysics scheme and data assimilation (DA) on the simulation of clouds, hydrometeors and radiative fluxes in the Arctic were investigated using the polar weather research and forecasting (WRF) model and three-dimensional variational DA. Compared with the WRF 5-class (WSM5) microphysics scheme, when the Morrison double-moment (Morrison) scheme was used, the simulated amount of cloud ice water decreased by approximately 68%. In contrast, the amount of water vapour, cloud liquid water, snow and rain in the atmosphere increased. With DA, the amount of water vapour increased, leading to increased hydrometeors. The cloud liquid water increased in the middle and low atmospheres when Morrison was used, whereas it increased in the low atmosphere when DA was used. The increase in cloud liquid water by using Morrison resulted in a decrease in the downward short-wave radiative flux at the surface, whereas using DA increased the downward long-wave radiative flux. Changing the microphysics scheme induced redistribution of the region and amounts of hydrometeors, whereas DA induced an increase in hydrometeors in specific regions by adding observation information to the model states.

An investigation into Venusian atmospheric chemistry based on an open-access photochemistry-transport model at 0–112 km

Atmospheric chemistry plays a crucial role in the evolution of climate habitability on Venus. It has been widely explored by chemistry-transport models, but some characteristics are still poorly interpreted. This study is devoted to developing an open-access chemistry-transport model spanning both the middle and lower atmospheres of Venus. It provides a scheme for the structure of the chemistry, especially for the sulfur and oxygen, and investigates the influence of the cloud diffusivity and the SO2 dissolution that are adopted in the clouds. The developed model is based on the VULCAN framework and was updated with the state-of-the-art Venusian atmospheric chemistry. It includes vertical eddy diffusion retrieved recently with the Venus Express observations, and it resolves radiative transfer containing gas absorption and scattering, Mie scattering of the cloud droplets, and absorption of the unknown UV absorber. The obtained abundance profiles of SO, SO2, CO, COS, O, O2, O3, HCl, and NO are in overall agreement with the observations. The results show that the increase in cloud diffusivity has slight effects on the chemical structure. The SO2 mainly dissolves in 50–90 km and evaporates below the clouds. The rapid dissolution-release cycle is responsible for the large upward flux of SO2 at 58 km. At around 70 km, SO has a significant peak that is larger than that of previous studies by an order of magnitude, and S and SO2 also show slight increases. They are attributed to the buffering effects of liquid SO2 in the clouds. O2 is significantly eliminated by SO in this layer. We emphasize the superior regulation of the sulfur cycle on O2 at 70 km and its potential contributions to the long-standing problem of the overestimated O2 abundance.

Effect of atmospheric density stratification on the generation of water waves by wind

The critical level instability mechanism for the generation of water waves by wind is re-examined for the situation when the atmosphere is density stratified. The density stratification is confined to the middle and upper atmosphere and then two cases are investigated. In case (A) no internal gravity waves are generated in the upper atmosphere and the effect of the density stratification is very small. In case (B) vertically propagating internal gravity waves form in the upper atmosphere and travel to infinity causing an energy loss, thus inhibiting the critical level instability in the lower atmosphere. Both cases are examined quantitatively for a logarithmic wind shear profile.

A Puzzling Quasi‐Periodic Variability in the Tropical Middle Atmosphere

AbstractThe Quasi‐Biennial Oscillation and the Semiannual Oscillation have been identified to be the leading modes of variability in the tropical middle atmosphere. With reanalysis data and independent rocket soundings from a low latitude site, we report the existence of yet another variability in the tropical lower mesosphere which is primarily evident as easterly bursts in zonal winds during the months of May‐July. It occurs with a variable interval of 2–5 yrs in the late 20th century and 7–9 yrs in the early 21st century. These Quasi‐Periodic Easterly Bursts are found to have remote influences on the Antarctic polar vortex as well as residual circulation in the lower mesosphere. We identify a potential causative mechanism for the easterly bursts that involve enhanced cross equatorial advection of momentum as well as gravity wave drag. A close association with Quasi Biennial Oscillation winds is observed, however, cause of the observed periodicity remains elusive.

The Climatology of Gravity Waves over the Low-Latitude Region Estimated by Multiple Meteor Radars

Atmospheric gravity waves (GWs) can strongly modulate middle atmospheric circulation and can be a significant factor for the coupling between the lower atmosphere and the middle atmosphere. GWs are difficult to resolve in global atmospheric models due to their small scale; thus, GW observations play an important role in middle atmospheric studies. The climatology of GW variance and momentum in the low-latitude mesosphere and lower thermosphere (MLT) region are revealed using multiple meteor radars, which are located at Kunming (25.6°N, 103.8°E), Sanya (18.4°N, 109.6°E), and Fuke (19.5°N, 109.1°E). The climatology and longitudinal variations in GW momentum fluxes and variance over the low-latitude region are reported. The GWs show strong seasonal variations and can greatly control the mesospheric horizontal winds via modulation of the quasi-geostrophic balance and momentum deposition. The different GW activities between Kunming and Sanya/Fuke are possibly consistent with the unique prevailing surface winds over Kunming and the convective system over the Tibetan Plateau according to the European Centre for Medium-Range Weather Forecasts (ECMWF), Reanalysis v5 (ERA5) data, and outgoing longwave radiation (OLR) data. These findings provide insight for better understanding the coupling between the troposphere and mesosphere.

Mesospheric Ozone Depletion during 2004–2024 as a Function of Solar Proton Events Intensity

Solar proton events (SPEs) affect the Earth’s atmosphere, causing additional ionization in the high-latitude mesosphere and stratosphere. Ionization rates from such solar proton events maximize in the stratosphere, but the formation of ozone-depleting nitrogen and hydrogen oxides begins at mesospheric altitudes. The destruction of mesospheric ozone is associated with protons with energies of about 10 MeV and higher and will strongly depend on the intensity of the flux of these particles. Most studies investigating the impact of SPEs on the characteristics of the middle atmosphere have been based on either simulations or reanalysis datasets, and some studies have used satellite observations to validate model results. We study the impact of SPEs on cold-season ozone loss in both the northern and southern hemispheres using Aura MLS mesospheric ozone measurements over the 2004 to 2024 period. Here, we show how strongly SPEs can deplete polar mesospheric ozone in different hemispheres and attempt to evaluate this dependence on the intensity of solar proton events. We found that moderate SPEs consisting of protons with an energy of more than 10 MeV and a flux intensity of more than 100 pfu destroy mesospheric ozone in the northern hemisphere up to 47% and in the southern hemisphere up to 33%. For both hemispheres, the peak of winter ozone loss was observed at about 76 km. In the northern hemisphere, maximum winter ozone loss was observed on the second day after a solar proton event, but in the southern hemisphere, winter ozone depletion was already detected on the first day. In the southern hemisphere, mesospheric ozone concentrations return to pre-event levels on the ninth day after a solar proton event, but in the northern hemisphere, even on the tenth day after a solar proton event, the mesospheric ozone layer may not be fully recovered. The strong SPEs with a proton flux intensity of more than 1000 pfu lead to a maximum winter ozone loss of up to 85% in the northern hemisphere, and in the southern hemisphere winter, ozone loss reaches 73%.

Movement of decaying quasi-2-day wave in the austral summer-time mesosphere

The quasi-2-day wave (Q2DW) is a large temperature disturbance in the austral summer-time mesosphere. Its decay-phase movement and phase-speed are poorly understood. Q2DW events observed by the microwave limb sounder (MLS) onboard the NASA Aura satellite reveal that during the temperature Q2DW’s decay-phase, the Q2DW temperature disturbance is still substantial, but its phase-speed reduces exponentially, on average, from around −70 to −20 m/s in around 30-days. Observations also reveal significant interannual variability in these phase-speed values. Q2DW events simulated by the extended Canadian middle atmosphere model (eCMAM) reveal that during the Q2DW decay-phase, the temperature Q2DW phase-speed is strongly correlated with the summer mesosphere easterly jet (SMEJ) wind values. eCMAM also shows that the phase-speed interannual variabilities partially reflect the interannual variabilities of the SMEJ. Planetary-wave diagnostics indicate that during the Q2DW’s decay phase, there is still an active transfer of momentum from the SMEJ into Q2DW. The model simulations therefore indicate that the SMEJ plays a substantial role in the decaying temperature Q2DW phase-speed variabilities observed by MLS. This adds to our knowledge of how the SMEJ affects the Q2DW.

Design of a three-frequency Rayleigh lidar for simultaneous temperature and wind measurements.

This study proposes what we believe to be a novel high-spectral-resolution three-frequency Rayleigh lidar for simultaneously measuring middle atmosphere temperature and wind. The temperature and wind could be retrieved without assuming an external reference temperature, as typical for a traditional Rayleigh Doppler lidar. Adopting a similar idea used in sodium temperature/wind lidar, this system alternatively emits laser pulses at three frequencies. It receives the corresponding Rayleigh backscattered signals filtered by an iodine cell as a frequency discriminator. The three frequencies are optimized based on the spectral characteristics resulting from the convolution of the pulse laser lineshape convolved Rayleigh scattering signal with iodine molecular absorption spectrum. A two-dimensional calibration curve for temperature and wind ratio is then generated from the theoretical calculation of the final convoluted spectra and used to retrieve temperature and wind simultaneously. Simulated with the return signals collected by a current broadband Rayleigh lidar (30-inch telescope and 15 W output laser power), the temperature and wind uncertainties with resolutions of 1 km and 1 hr are estimated to be 0.4 K and 0.35 m/s, respectively, at 30 km and increase to 16.3 K and 8.1 m/s at 70 km.