Search

This study explores the abrupt split of the polar vortex in the upper stratosphere prior to a recent sudden stratospheric warming event on 5 January 2021 (SSW21) and the mechanisms of vortex preconditioning by using the Modern-Era Retrospective Analysis for Research and Applications version 2 (MERRA2) global reanalysis data. SSW21 is preceded by the highly distorted polar vortex that was initially displaced from the pole but eventually split at the onset date. Vortex splitting is most significant in the upper stratosphere (1 hPa altitude) accompanied by the anomalous growth of westward-propagating planetary waves (PWs) of zonal wavenumber (ZWN) 2 (WPW2). While previous studies have suggested the East Asian trough as a potential source for the abnormal WPW2 growth, the prominent westward-propagating nature cannot be explained satisfactorily by the upward propagation of the quasi-stationary ZWN2 fluxes in the troposphere. More importantly, WPW2 exhibits an obvious in situ excitation signature within the barotropically and baroclinically destabilized stratosphere, dominated by the easterlies descending from the stratopause containing the WPW2 critical levels. This suggests that the vortex split is attributed to the WPW2 generated in situ within the stratosphere via instability. Vortex destabilization is achieved as the double-jet structure consisting of a subtropical mesospheric core and a polar stratospheric core develops into SSW21 by encouraging the anomalous dissipation of the upward-propagating tropospheric ZWN1 PWs. This double-jet configuration is likely a favorable precursor for SSW onset, not only for the SSW21 but generally for most SSWs, through promoting the anomalous growth of unstable PWs as well as the enhancement of the tropospheric PW dissipation.

In January 2020, unexpected easterly winds developed in the downward-propagating westerly quasi-biennial oscillation (QBO) phase. This event corresponds to the second QBO disruption in history, and it occurred 4 years after the first disruption of 2015/16. According to several previous studies, strong midlatitude Rossby waves propagating from the Southern Hemisphere (SH) during the SH winter likely initiated the disruption; nevertheless, the wave forcing that finally led to the disruption has not been investigated. In this study, we examine the role of equatorial waves and small-scale convective gravity waves (CGWs) in the 2019/20 QBO disruption using Modern-Era Retrospective Analysis for Research and Applications version 2 (MERRA-2) global reanalysis data. In June–September 2019, unusually strong Rossby wave forcing originating from the SH decelerated the westerly QBO at 0–5∘ N at ∼50 hPa. In October–November 2019, vertically (horizontally) propagating Rossby waves and mixed Rossby–gravity (MRG) waves began to increase (decrease). From December 2019, the contribution of the MRG wave forcing to the zonal wind deceleration was the largest, followed by the Rossby wave forcing originating from the Northern Hemisphere and the equatorial troposphere. In January 2020, CGWs provided 11 % of the total negative wave forcing at ∼43 hPa. Inertia–gravity (IG) waves exhibited a moderate contribution to the negative forcing throughout. Although the zonal mean precipitation was not significantly larger than the climatology, convectively coupled equatorial wave activities were increased during the 2019/20 disruption. As in the 2015/16 QBO disruption, the increased barotropic instability at the QBO edges generated more MRG waves at 70–90 hPa, and westerly anomalies in the upper troposphere allowed more westward IG waves and CGWs to propagate to the stratosphere. Combining the 2015/16 and 2019/20 disruption cases, Rossby waves and MRG waves can be considered the key factors inducing QBO disruption.

In February 2016, the descent of the westerly phase of the quasi-biennial oscillation (QBO) was unprecedentedly disrupted by the development of easterly winds. Previous studies have shown that extratropical Rossby waves propagating into the deep tropics were the major cause of the 2015/16 QBO disruption. However, a large portion of the negative momentum forcing associated with the disruption still stems from equatorial planetary and small-scale gravity waves, which calls for detailed analyses by separating each wave mode compared with climatological QBO cases. Here, the contributions of resolved equatorial planetary waves (Kelvin, Rossby, mixed Rossby–gravity (MRG), and inertia–gravity (IG) waves) and small-scale convective gravity waves (CGWs) obtained from an offline CGW parameterization to the 2015/16 QBO disruption are investigated using MERRA-2 global reanalysis data from October 2015 to February 2016. In October and November 2015, anomalously strong negative forcing by MRG and IG waves weakened the QBO jet at 0–5∘ S near 40 hPa, leading to Rossby wave breaking at the QBO jet core in the Southern Hemisphere. From December 2015 to January 2016, exceptionally strong Rossby waves propagating horizontally (vertically) continuously decelerated the southern (northern) flank of the jet. In February 2016, when the westward CGW momentum flux at the source level was much stronger than its climatology, CGWs began to exert considerable negative forcing at 40–50 hPa near the Equator, in addition to the Rossby waves. The enhancement of the negative wave forcing in the tropics stems mostly from strong wave activity in the troposphere associated with increased convective activity and the strong westerlies (or weaker easterlies) in the troposphere, except that the MRG wave forcing is more likely associated with increased barotropic instability in the lower stratosphere.

Effects of realistic propagation of gravity waves (GWs) on distribution of GW pseudomomentum fluxes are explored using a global ray-tracing model for the 2009 sudden stratospheric warming (SSW) event. Four-dimensional (4D; x–z and t) and two-dimensional (2D; z and t) results are compared for various parameterized pseudomomentum fluxes. In ray-tracing equations, refraction due to horizontal wind shear and curvature effects are found important and comparable to one another in magnitude. In the 4D, westward pseudomomentum fluxes are enhanced in the upper troposphere and northern stratosphere due to refraction and curvature effects around fluctuating jet flows. In the northern polar upper mesosphere and lower thermosphere, eastward pseudomomentum fluxes are increased in the 4D. GWs are found to propagate more to the upper atmosphere in the 4D, since horizontal propagation and change in wave numbers due to refraction and curvature effects can make it more possible that GWs elude critical level filtering and saturation in the lower atmosphere. GW focusing effects occur around jet cores, and ray-tube effects appear where the polar stratospheric jets vary substantially in space and time. Enhancement of the structure of zonal wave number 2 in pseudomomentum fluxes in the middle stratosphere begins from the early stage of the SSW evolution. An increase in pseudomomentum fluxes in the upper atmosphere is present even after the onset in the 4D. Significantly enhanced pseudomomentum fluxes, when the polar vortex is disturbed, are related to GWs with small intrinsic group velocity (wave capture), and they would change nonlocally nearby large-scale vortex structures without substantially changing local mean flows.

Convection as one dominant source of atmospheric gravity waves (GWs) has been the focus of investigation over recent years. However, its spatial and temporal forcing scales are not well known. In this work we address this open issue by a systematic verification of free parameters of the Yonsei convective GW source scheme based on observations from the High Resolution Dynamics Limb Sounder (HIRDLS). The instrument can only see a limited portion of the gravity wave spectrum due to visibility effects and observation geometry. To allow for a meaningful comparison of simulated GWs to observations, a comprehensive filter, which mimics the instrument limitations, is applied to the simulated waves. By this approach, only long horizontal-scale convective GWs are addressed. Results show that spectrum, distribution of momentum flux, and zonal mean forcing of long horizontal-scale convective GWs can be successfully simulated by the superposition of three or four combinations of parameter sets reproducing the observed GW spectrum. These selected parameter sets are different for northern and southern summer. Although long horizontal-scale waves are only part of the full spectrum of convective GWs, the momentum flux of these waves is found to be significant and relevant for the driving of the QBO (quasi-biennial oscillation). The zonal momentum balance is considered in vertical cross sections of GW momentum flux (GWMF) and GW drag (GWD). Global maps of the horizontal distribution of GWMF are considered and consistency between simulated results and HIRDLS observations is found. The latitude dependence of the zonal phase speed spectrum of GWMF and its change with altitude is discussed.

An idealized baroclinic instability case is simulated using a ∼ 10 km resolution global model to investigate the characteristics of gravity waves generated in the baroclinic life cycle. Three groups of gravity waves appear around the high-latitude surface trough at the mature stage of the baroclinic wave. They have horizontal and vertical wavelengths of 40–400 and 2.9–9.8 km, respectively, in the upper troposphere. The two-dimensional phase-velocity spectrum of the waves is arc shaped with a peak at 17 m s−1 eastward. These waves have difficulty in propagating upward through the tropospheric westerly jet. At the breaking stage of the baroclinic wave, a midlatitude surface low is isolated from the higher-latitude trough, and two groups of quasi-stationary gravity waves appear near the surface low. These waves have horizontal and vertical wavelengths of 60–400 and 4.9–14 km, respectively, and are able to propagate vertically for long distances. The simulated gravity waves seem to be generated by surface fronts, given that the structures and speeds of wave phases are coherent with those of the fronts.

The momentum forcing of the QBO (quasi-biennial oscillation) by equatorial waves is estimated using recent reanalyses. Based on the estimation using the conventional pressure-level data sets, the forcing by the Kelvin waves (3–9 m s−1 month−1) dominates the net forcing by all equatorial wave modes (3–11 m s−1 month−1) in the easterly-to-westerly transition phase at 30 hPa. In the opposite phase, the net forcing by equatorial wave modes is small (1–5 m s−1 month−1). By comparing the results with those from the native model-level data set of the ERA-Interim reanalysis, it is suggested that the use of conventional-level data causes the Kelvin wave forcing to be underestimated by 2–4 m s−1 month−1. The momentum forcing by mesoscale gravity waves, which are unresolved in the reanalyses, is deduced from the residual of the zonal wind tendency equation. In the easterly-to-westerly transition phase at 30 hPa, the mesoscale gravity wave forcing is found to be smaller than the resolved wave forcing, whereas the gravity wave forcing dominates over the resolved wave forcing in the opposite phase. Finally, we discuss the uncertainties in the wave forcing estimates using the reanalyses.

The characteristics of horizontal divergence induced by typhoon-generated gravity waves (HDTGWs) and the influence of HDTGW on typhoon evolution are investigated based on the simulation results of Typhoon Saomai (2006) using the Weather Research and Forecasting (WRF) model. The power spectral density of HDTGW shows dominant powers at horizontal wavelengths of 20–30 km and at periods of less than 1 h. This is associated with gravity waves generated by vigorous convective clouds in an inner core region of the typhoon. However, the domain-averaged HDTGW in the upper troposphere and lower stratosphere had a spectral peak at 24 h, which is well correlated with the minimum sea-level pressure of the typhoon, especially during a rapidly developing period. The 24 h period of the averaged HDTGW stems from the inertia–gravity waves generated by the convective clouds in the spiral rainbands, and showed no clear association with the thermal tides or the diurnal variation of precipitation.

For a case study of Typhoon Ewiniar performed with a mesoscale model, we compare stratospheric gravity wave (GW) momentum flux determined from temperature variances by applying GW polarization relations and by assuming upward propagating waves, with GW momentum flux calculated from model winds which is considered as a reference. The temperature-based momentum-flux profile exhibits positive biases relative to the reference, which fluctuate significantly with altitude. The vertically-averaged magnitude of the positive biases is about 14% of the reference momentum flux. We found that this deviation from the reference stems from the interference between upward and downward propagating waves. The downward propagating GWs are due mainly to partial reflections of upward propagating waves at altitudes where the background wind and stability change with height. When the upward and downward propagating waves are decomposed and their momentum fluxes are calculated separately from temperature perturbations, the fraction of the momentum flux arising from the downward propagating waves is about 4.5–8.2% of that from the upward propagating waves. The net momentum flux of upward and downward propagating GWs agrees well with the reference from the model wind perturbations. The implications of this study for the GW momentum-flux observations from satellites are discussed.