The January 2022 Hunga Tonga–Hunga Ha’apai eruption was one of the most explosive volcanic events of the modern era^1,2, producing a vertical plume that peaked more than 50 km above the Earth³. The initial explosion and subsequent plume triggered atmospheric waves that propagated around the world multiple times⁴. A global-scale wave response of this magnitude from a single source has not previously been observed. Here we show the details of this response, using a comprehensive set of satellite and ground-based observations to quantify it from surface to ionosphere. A broad spectrum of waves was triggered by the initial explosion, including Lamb waves^5,6 propagating at phase speeds of 318.2 ± 6 m s⁻¹ at surface level and between 308 ± 5 to 319 ± 4 m s⁻¹ in the stratosphere, and gravity waves⁷ propagating at 238 ± 3 to 269 ± 3 m s⁻¹ in the stratosphere. Gravity waves at sub-ionospheric heights have not previously been observed propagating at this speed or over the whole Earth from a single source^8,9. Latent heat release from the plume remained the most significant individual gravity wave source worldwide for more than 12 h, producing circular wavefronts visible across the Pacific basin in satellite observations. A single source dominating such a large region is also unique in the observational record. The Hunga Tonga eruption represents a key natural experiment in how the atmosphere responds to a sudden point-source-driven state change, which will be of use for improving weather and climate models.

Two of the most basic parameters generated from a reanalysis are temperature and winds. Temperatures in the reanalyses are derived from conventional (surface and balloon), aircraft, and satellite observations. Winds are observed by conventional systems, cloud tracked, and derived from height fields, which are in turn derived from the vertical temperature structure. In this paper we evaluate as part of the SPARC Reanalysis Intercomparison Project (S-RIP) the temperature and wind structure of all the recent and past reanalyses. This evaluation is mainly among the reanalyses themselves, but comparisons against independent observations, such as HIRDLS and COSMIC temperatures, are also presented. This evaluation uses monthly mean and 2.5° zonal mean data sets and spans the satellite era from 1979–2014. There is very good agreement in temperature seasonally and latitudinally among the more recent reanalyses (CFSR, MERRA, ERA-Interim, JRA-55, and MERRA-2) between the surface and 10 hPa. At lower pressures there is increased variance among these reanalyses that changes with season and latitude. This variance also changes during the time span of these reanalyses with greater variance during the TOVS period (1979–1998) and less variance afterward in the ATOVS period (1999–2014). There is a distinct change in the temperature structure in the middle and upper stratosphere during this transition from TOVS to ATOVS systems. Zonal winds are in greater agreement than temperatures and this agreement extends to lower pressures than the temperatures. Older reanalyses (NCEP/NCAR, NCEP/DOE, ERA-40, JRA-25) have larger temperature and zonal wind disagreement from the more recent reanalyses. All reanalyses to date have issues analysing the quasi-biennial oscillation (QBO) winds. Comparisons with Singapore QBO winds show disagreement in the amplitude of the westerly and easterly anomalies. The disagreement with Singapore winds improves with the transition from TOVS to ATOVS observations. Temperature bias characteristics determined via comparisons with a reanalysis ensemble mean (MERRA, ERA-Interim, JRA-55) are similarly observed when compared with Aura HIRDLS and Aura MLS observations. There is good agreement among the NOAA TLS, SSU1, and SSU2 Climate Data Records and layer mean temperatures from the more recent reanalyses. Caution is advised for using reanalysis temperatures for trend detection and anomalies from a long climatology period as the quality and character of reanalyses may have changed over time.

Atmospheric gravity waves (GWs) are key drivers of the atmospheric circulation, but their representation in general circulation models (GCMs) is challenging, leading to significant biases in middle atmospheric circulations. Unresolved GW momentum transport in GCMs must be parameterized, but global directional GW observations are needed to constrain this. Here we present an 18‐year climatology of directional stratospheric GW momentum flux (GWMF) from global AIRS/Aqua 3‐D satellite observations during 2002 to 2019. Striking hemispheric asymmetries are found at high latitudes, including dramatic reductions and reversals of GWMF during sudden stratospheric warmings. During Southern Hemisphere winter, a lateral convergence of GWMF toward 60°S is found that has no Northern Hemisphere counterpart. In the tropics, we find that zonal GWMF in AIRS measurements is strongly modulated by the semiannual oscillation (SAO) but not the quasi‐biennial oscillation (QBO). Our results provide guidance for future GW parameterizations needed to resolve long‐standing biases in GCMs.

Gravity waves (GWs) transport momentum and energy in the atmosphere, exerting a profound influence on the global circulation. Accurately measuring them is thus vital both for understanding the atmosphere and for developing the next generation of weather forecasting and climate prediction models. However, it has proven very difficult to measure the full set of GW parameters from satellite measurements, which are the only suitable observations with global coverage. This is particularly critical at latitudes close to 60° S, where climate models significantly under-represent wave momentum fluxes. Here, we present a novel fully 3-D method for detecting and characterising GWs in the stratosphere. This method is based around a 3-D Stockwell transform, and can be applied retrospectively to existing observed data. This is the first scientific use of this spectral analysis technique. We apply our method to high-resolution 3-D atmospheric temperature data from AIRS/Aqua over the altitude range 20–60 km. Our method allows us to determine a wide range of parameters for each wave detected. These include amplitude, propagation direction, horizontal/vertical wavelength, height/direction-resolved momentum fluxes (MFs), and phase and group velocity vectors. The latter three have not previously been measured from an individual satellite instrument. We demonstrate this method over the region around the Southern Andes and Antarctic Peninsula, the largest known sources of GW MFs near the 60° S belt. Our analyses reveal the presence of strongly intermittent highly directionally focused GWs with very high momentum fluxes (∼ 80–100 mPa or more at 30 km altitude). These waves are closely associated with the mountains rather than the open ocean of the Drake Passage. Measured fluxes are directed orthogonal to both mountain ranges, consistent with an orographic source mechanism, and are largest in winter. Further, our measurements of wave group velocity vectors show clear observational evidence that these waves are strongly focused into the polar night wind jet, and thus may contribute significantly to the missing momentum at these latitudes. These results demonstrate the capabilities of our new method, which provides a powerful tool for delivering the observations required for the next generation of weather and climate models.

Nearly all general circulation models significantly fail to reproduce the observed behaviour of the southern wintertime polar vortex. It has been suggested that these biases result from an underestimation of gravity wave drag on the atmosphere at latitudes near 60° S, especially around the "hot spot" of intense gravity wave fluxes above the mountainous Southern Andes and Antarctic peninsula. Here, we use Global Positioning System radio occultation (GPS-RO) data from the COSMIC satellite constellation to determine the properties of gravity waves in the hot spot and beyond. We show considerable southward propagation to latitudes near 60° S of waves apparently generated over the southern Andes. We propose that this propagation may account for much of the wave drag missing from the models. Furthermore, there is a long leeward region of increased gravity wave energy that sweeps eastwards from the mountains over the Southern Ocean. Despite its striking nature, the source of this region has historically proved difficult to determine. Our observations suggest that this region includes both waves generated locally and orographic waves advected downwind from the hot spot. We describe and use a new wavelet-based analysis technique for the quantitative identification of individual waves from COSMIC temperature profiles. This analysis reveals different geographical regimes of wave amplitude and short-timescale variability in the wave field over the Southern Ocean. Finally, we use the increased numbers of closely spaced pairs of profiles from the deployment phase of the COSMIC constellation in 2006 to make estimates of gravity wave horizontal wavelengths. We show that, given sufficient observations, GPS-RO can produce physically reasonable estimates of stratospheric gravity wave momentum flux in the hot spot that are consistent with measurements made by other techniques. We discuss our results in the context of previous satellite and modelling studies and explain how they advance our understanding of the nature and origins of waves in the southern stratosphere.

The quasi‐biennial oscillation (QBO) of tropical stratospheric winds was disrupted during the 2019/20 Northern Hemisphere winter. We show that this latest disruption to the regular QBO cycling was similar in many respects to that seen in 2016, but initiated by horizontal momentum transport from the Southern Hemisphere. The predictable signal associated with the QBO's quasi‐regular phase progression is lost during disruptions and the oscillation reemerges after a few months significantly shifted in phase from what would be expected if it had progressed uninterrupted. We infer from an increased wave‐momentum flux into equatorial latitudes seen in climate model projections that disruptions to the QBO are likely to become more common in future. Consequently, it is possible that in the future, the QBO could be a less reliable source of predictability on lead times extending out to several years than it currently is.