The January 2022 Hunga Tonga–Hunga Ha’apai eruption was one of the most explosive volcanic events of the modern era1,2, producing a vertical plume that peaked more than 50 km above the Earth3. The initial explosion and subsequent plume triggered atmospheric waves that propagated around the world multiple times4. 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 waves5,6 propagating at phase speeds of 318.2 ± 6 m s−1 at surface level and between 308 ± 5 to 319 ± 4 m s−1 in the stratosphere, and gravity waves7 propagating at 238 ± 3 to 269 ± 3 m s−1 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 source8,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.

Abstract. 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.

Abstract. 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.

[1] Three years of gravity wave observations from the High Resolution Dynamics Limb Sounder instrument on NASA's Aura satellite are examined. We produce estimates of the global distribution of gravity wave momentum flux as a function of individual observed wave packets. The observed distribution at the 25 km altitude level is dominated by the small proportion of wave packets with momentum fluxes greater than 0.5 mPa. Depending on latitude and season, these wave packets only comprise 7-25% of observations, but are shown to be almost entirely responsible for the morphology of the observed global momentum flux distribution. Large-amplitude wave packets are found to be more important over orographic regions than over flat ocean regions, and to be especially high in regions poleward of 40 ı S during austral winter. The momentum flux carried by the largest packets relative to the distribution mean is observed to decrease with height over orographic wave generation regions, but to increase with height at tropical latitudes; the mesospheric intermittency resulting is broadly equivalent in both cases. Consistent with previous studies, waves in the top 10% of the extratropical distribution are observed to carry momentum fluxes more than twice the mean and waves in the top 1% more than 10 the mean, and the Gini coefficient is found to characterize the observed distributions well. These results have significant implications for gravity wave modeling.Citation: Wright, C. J., S. M. Osprey, and J. C. Gille (2013), Global observations of gravity wave intermittency and its impact on the observed momentum flux morphology,

[1] Observations from the High Resolution Dynamics Limb Sounder (HIRDLS) instrument on NASA's Aura satellite are used to quantify gravity wave momentum fluxes in the middle atmosphere. The period around the 2006 Arctic sudden stratospheric warming (SSW) is investigated, during which a substantial elevation of the stratopause occurred. Analysis of the HIRDLS results, together with analysis of European Centre for Medium-Range Weather Forecasting zonal winds, provide direct evidence of wind filtering of the gravity wave spectrum during this period. This confirms previous hypotheses from model studies and further contributes to our understanding of the effects of gravity wave driving on the winter polar stratopause.

Abstract. 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)(1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998) and less variance afterward in the ATOVS period (1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(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, ERAInterim, 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.

No abstract

International audienceUsing the world's largest data set of in situ ocean current measurements, combined with a high-resolution topography roughness data set, we use a model-assisted hierarchical clustering methodology to estimate the global lee wave generation rate at the ocean floor. Our analysis suggests that internal wave generation contributes 0.75±0.19 TW (±2 standard deviation) to the oceanic energy budget but with a strong dependence on the Brunt-Väisäla (buoyancy) frequency climatology used. This estimate is higher than previous calculations and suggests that internal wave generation may be a much more significant contributor to the global oceanic mechanical energy budget than had previously been assumed. Our results imply that lee wave generation and propagation may be a dominant sink of at least half and potentially the overwhelming majority of ocean surface wind work on the geostrophic circulation