Atmospheric gravity waves (GWs) play a key role in determining the thermodynamical structure of the Earth’s middle atmosphere. Despite the small spatial and temporal scales of these waves, a few high-top general circulation models (GCMs) that can resolve them explicitly have recently become available. This study compares global GW characteristics simulated in one such GCM, the Japanese Atmospheric GCM for Upper-Atmosphere Research (JAGUAR), with those derived from three-dimensional (3-D) temperatures observed by the Atmospheric Infrared Sounder (AIRS) aboard NASA’s Aqua satellite. The target period is from 15 December 2018 to 8 January 2019, including the onset of a major sudden stratospheric warming (SSW). The 3-D Stockwell transform method is used for GW spectral analysis. The amplitudes and momentum fluxes of GWs in JAGUAR are generally in good quantitative agreement with those in the AIRS observations in both magnitude and distribution. As the SSW event progressed, the GW amplitudes and eastward momentum flux increased at low latitudes in the summer hemisphere in both the model and observation datasets. Case studies demonstrate that the model is able to reproduce comparable wave events to those in the AIRS observations with some differences, especially noticeable at low latitudes in the summer hemisphere. Through a comparison between the model results with and without the AIRS observational filter applied, it is suggested that the amplitudes of GWs near the exits and entrances of eastward jet streaks are underestimated in AIRS observations.

An international joint research project, entitled Interhemispheric Coupling Study by Observations and Modelling (ICSOM), is ongoing. In the late 2000s, an interesting form of interhemispheric coupling (IHC) was discovered: when warming occurs in the winter polar stratosphere, the upper mesosphere in the summer hemisphere also becomes warmer with a time lag of days. This IHC phenomenon is considered to be a coupling through processes in the middle atmosphere (i.e., stratosphere, mesosphere, and lower thermosphere). Several plausible mechanisms have been proposed so far, but they are still controversial. This is mainly because of the difficulty in observing and simulating gravity waves (GWs) at small scales, despite the important role they are known to play in middle atmosphere dynamics. In this project, by networking sparsely but globally distributed radars, mesospheric GWs have been simultaneously observed in seven boreal winters since 2015/16. We have succeeded in capturing five stratospheric sudden warming events and two polar vortex intensification events. This project also includes the development of a new data assimilation system to generate long-term reanalysis data for the whole middle atmosphere, and simulations by a state-of-art GW-permitting general circulation model using reanalysis data as initial values. By analyzing data from these observations, data assimilation, and model simulation, comprehensive studies to investigate the mechanism of IHC are planned. This paper provides an overview of ICSOM, but even initial results suggest that not only gravity waves but also large-scale waves are important for the mechanism of the IHC.

Since 2004, following prolonged stratospheric sudden warming (SSW) events, it has been observed that the stratopause disappeared and reformed at a higher altitude, forming an elevated stratopause (ES). The relative roles of atmospheric waves in the mechanism of ES formation are still not fully understood. We performed a hindcast of the 2018/19 SSW event using a gravity-wave (GW) permitting general circulation model that resolves the mesosphere and lower thermosphere (MLT) and analyzed dynamical phenomena throughout the entire middle atmosphere. An ES formed after the major warming on 1 January 2019. There was a marked temperature maximum in the polar upper mesosphere around 28 December 2018 prior to the disappearance of the descending stratopause associated with the SSW. This temperature structure is referred to as a mesospheric inversion layer (MIL). We show that adiabatic heating from the residual circulation driven by GW forcing (GWF) causes barotropic and/or baroclinic instability before the MIL formation, causing in situ generation of planetary waves (PWs). These PWs propagate into the MLT and exert negative (westward) forcing, which contributes to the MIL formation. Both GWF and PW forcing (PWF) above the recovered eastward jet play crucial roles in ES formation. The altitude of the recovered eastward jet, which regulates GWF and PWF height, is likely affected by the MIL structure. Simple vertical propagation from the lower atmosphere is insufficient to explain the presence of the GWs observed in this event.

Observations with high vertical resolution have revealed that power spectra of horizontal wind and temperature fluctuations versus vertical wavenumber m have a universal shape with a steep slope in a high m range, approximately proportional to m^(-3). Several theoretical models explaining this spectral slope were proposed under an assumption of gravity wave (GW) saturation. However, little evidence has been obtained to show that these universal spectra are fully composed of GWs. To confirm the validity of this assumption, two kinds of m spectra are calculated using outputs from a GW-permitting high-top general circulation model. One is the spectra for GWs designated by fluctuations having total horizontal wavenumbers of 21–639. The other is the spectra of fluctuations unfiltered except extracting a linear trend in the vertical that were often analyzed in the observational studies. Comparison between the two shows that GWs dominate the observed spectra only in a higher m part of the steep slope, whereas disturbances other than GWs significantly contribute to its lower m part. Moreover, geographical distributions of the characteristic wavenumbers, slopes, and spectral densities of GW spectra are examined for several divided height regions of the whole middle atmosphere. It is shown that strong vertical shear below the zonal wind jets as well as the wave saturation is responsible for the formation of the steep slopes of GW spectra.