Search

The Hunga Tonga–Hunga Ha′apai volcano erupted on 15 January 2022, launching Lamb waves and gravity waves into the atmosphere. In this study, we present results using 13 globally distributed meteor radars and identify the volcanogenic gravity waves in the mesospheric/lower thermospheric winds. Leveraging the High-Altitude Mechanistic general Circulation Model (HIAMCM), we compare the global propagation of these gravity waves. We observed an eastward-propagating gravity wave packet with an observed phase speed of 240 ± 5.7 m s−1 and a westward-propagating gravity wave with an observed phase speed of 166.5 ± 6.4 m s−1. We identified these waves in HIAMCM and obtained very good agreement of the observed phase speeds of 239.5 ± 4.3 and 162.2 ± 6.1 m s−1 for the eastward the westward waves, respectively. Considering that HIAMCM perturbations in the mesosphere/lower thermosphere were the result of the secondary waves generated by the dissipation of the primary gravity waves from the volcanic eruption, this affirms the importance of higher-order wave generation. Furthermore, based on meteor radar observations of the gravity wave propagation around the globe, we estimate the eruption time to be within 6 min of the nominal value of 15 January 2022 04:15 UTC, and we localized the volcanic eruption to be within 78 km relative to the World Geodetic System 84 coordinates of the volcano, confirming our estimates to be realistic.

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.

Atmospheric gravity waves (GWs) play an important role in atmospheric dynamics but accurately representing them in general circulation models (GCMs) is challenging. This is especially true for orographic GWs generated by wind flow over small mountainous islands in the Southern Ocean. Currently, these islands lie in the “grey zone” of global model resolution, where they are neither fully resolved nor fully parameterised. It is expected that as GCMs approach the spatial resolution of current high-resolution local-area models, small-island GW sources may be resolved without the need for parameterisations. But how realistic are the resolved GWs in these high-resolution simulations compared to observations? Here, we test a high-resolution (1.5 km horizontal grid, 118 vertical levels) local-area configuration of the Met Office Unified Model over the mountainous island of South Georgia (54∘ S, 36∘ W), running without GW parameterisations. The island's orography is well resolved in the model, and real-time boundary conditions are used for two time periods during July 2013 and June–July 2015. We compare simulated GWs in the model to coincident 3-D satellite observations from the Atmospheric Infrared Sounder (AIRS) on board Aqua. By carefully sampling the model using the AIRS resolution and measurement footprints (denoted as model sampled as AIRS hereafter), we present the first like-for-like comparison of simulated and observed 3-D GW amplitudes, wavelengths and directional GW momentum flux (GWMF) over the island using a 3-D S-transform method. We find that the timing, magnitude and direction of simulated GWMF over South Georgia are in good general agreement with observations, once the AIRS sampling and resolution are applied to the model. Area-averaged zonal GWMF during these 2 months is westward at around 5.3 and 5.6 mPa in AIRS and model sampled as AIRS datasets respectively, but values directly over the island can exceed 50 mPa. However, up to 35 % of the total GWMF in AIRS is actually found upwind of the island compared to only 17 % in the model sampled as AIRS, suggesting that non-orographic GWs observed by AIRS may be underestimated in our model configuration. Meridional GWMF results show a small northward bias (∼20 %) in the model sampled as AIRS that may correspond to a southward wind bias compared to coincident radiosonde measurements. Finally, we present one example of large-amplitude (T′≈15–20 K at 45 km altitude) GWs at short horizontal wavelengths (λH≈30–40 km) directly over the island in AIRS measurements that show excellent agreement with the model sampled as AIRS. This suggests that orographic GWs in the full-resolution model with T′≈45 K and λH≈30–40 km can occur in reality. Our study demonstrates that not only can high-resolution local-area models simulate realistic stratospheric GWs over small mountainous islands but the application of satellite sampling and resolution to these models can also be a highly effective method for their validation.

An important source of polar stratospheric clouds (PSCs), which play a crucial role in controlling polar stratospheric ozone depletion, is the temperature fluctuations induced by mountain waves. These enable stratospheric temperatures to fall below the threshold value for PSC formation in regions of negative temperature perturbations or cooling phases induced by the waves even if the synoptic-scale temperatures are too high. However, this formation mechanism is usually missing in global chemistry–climate models because these temperature fluctuations are neither resolved nor parameterised. Here, we investigate in detail the episodic and localised wintertime stratospheric cooling events produced over the Antarctic Peninsula by a parameterisation of mountain-wave-induced temperature fluctuations inserted into a 30-year run of the global chemistry–climate configuration of the UM-UKCA (Unified Model – United Kingdom Chemistry and Aerosol) model. Comparison of the probability distribution of the parameterised cooling phases with those derived from climatologies of satellite-derived AIRS brightness temperature measurements and high-resolution radiosonde temperature soundings from Rothera Research Station on the Antarctic Peninsula shows that they broadly agree with the AIRS observations and agree well with the radiosonde observations, particularly in both cases for the “cold tails” of the distributions. It is further shown that adding the parameterised cooling phase to the resolved and synoptic-scale temperatures in the UM-UKCA model results in a considerable increase in the number of instances when minimum temperatures fall below the formation temperature for PSCs made from ice water during late austral autumn and early austral winter and early austral spring, and without the additional cooling phase the temperature rarely falls below the ice frost point temperature above the Antarctic Peninsula in the model. Similarly, it was found that the formation potential for PSCs made from ice water was many times larger if the additional cooling is included. For PSCs made from nitric acid trihydrate (NAT) particles it was only during October that the additional cooling is required for temperatures to fall below the NAT formation temperature threshold (despite more NAT PSCs occurring during other months). The additional cooling phases also resulted in an increase in the surface area density of NAT particles throughout the winter and early spring, which is important for chlorine activation. The parameterisation scheme was finally shown to make substantial differences to the distribution of total column ozone during October, resulting from a shift in the position of the polar vortex.

Atmospheric gravity waves play a key role in the transfer of energy and momentum between layers of the Earth's atmosphere. However, nearly all general circulation models (GCMs) seriously under-represent the momentum fluxes of gravity waves at latitudes near 60∘ S, which can lead to significant biases. A prominent example of this is the “cold pole problem”, where modelled winter stratospheres are unrealistically cold. There is thus a need for large-scale measurements of gravity wave fluxes near 60∘ S, and indeed globally, to test and constrain GCMs. Such measurements are notoriously difficult, because they require 3-D observations of wave properties if the fluxes are to be estimated without using significant limiting assumptions. Here we use 3-D satellite measurements of stratospheric gravity waves from NASA's Atmospheric Infrared Sounder (AIRS) Aqua instrument. We present the first extended application of a 3-D Stockwell transform (3DST) method to determine localised gravity wave amplitudes, wavelengths and directions of propagation around the entire region of the Southern Ocean near 60∘ S during austral winter 2010. We first validate our method using a synthetic wavefield and two case studies of real gravity waves over the southern Andes and the island of South Georgia. A new technique to overcome wave amplitude attenuation problems in previous methods is also presented. We then characterise large-scale gravity wave occurrence frequencies, directional momentum fluxes and short-timescale intermittency over the entire Southern Ocean. Our results show that highest wave occurrence frequencies, amplitudes and momentum fluxes are observed in the stratosphere over the mountains of the southern Andes and Antarctic Peninsula. However, we find that around 60 %–80 % of total zonal-mean momentum flux is located over the open Southern Ocean during June–August, where a large “belt” of increased wave occurrence frequencies, amplitudes and fluxes is observed. Our results also suggest significant short-timescale variability of fluxes from both orographic and non-orographic sources in the region. A particularly striking result is a widespread convergence of gravity wave momentum fluxes towards latitudes around 60∘ S from the north and south. We propose that this convergence, which is observed at nearly all longitudes during winter, could account for a significant part of the under-represented flux in GCMs at these latitudes.