1. Introduction
Occurring in ∼50% of boreal winters, an enhancement of quasi-stationary planetary waves of zonal wavenumbers 1 and 2 can result in sudden stratospheric warmings (SSWs) (e.g., Butler et al. 2015). Upon dissipation, these waves strongly decelerate the stratospheric flow, inducing an overturning mean meridional circulation that adiabatically warms the polar region. The anomalous polar warming causes the stratopause to descend below its climatological altitude (Matsuno 1971). The coinciding wavenumber-1 and/or wavenumber-2 pattern can project onto a dominant mode of climate variability in the troposphere called the Arctic Oscillation, aka the Northern Annular Mode (e.g., Baldwin and Dunkerton 2001). Consequently, after a 10–30-day period following the stratospheric vortex disruption, the jet stream tends to shift equatorward over the Atlantic and anomalously cold conditions prevail over Europe and northeast America.
As the perturbed polar vortex begins to recover from SSW, a new stratopause can reform at an altitude level at least 10 km above its norm (e.g., Manney et al. 2008; Siskind et al. 2010). These “elevated stratopause” SSW events (or ES-SSWs) reflect the strong coupling between the stratosphere and the mesosphere–lower thermosphere (MLT) region. This coupling is exemplified by an unusually strong polar downwelling during the stratopause reformation that transports long-lived tracers found in the thermosphere well into the stratosphere (e.g., Orsolini et al. 2022, 2017). During climatological wintertime conditions, downwelling over the pole is induced by westward GW drag near the stratopause. However, ES-SSW events give rise to the presence of strong traveling planetary (Rossby) waves around the time of ES-SSW onset (e.g., Iida et al. 2014; Limpasuvan et al. 2016). These waves can impose a significant impact on the circulation of the MLT (e.g., Rhodes et al. 2021; Sassi et al. 2016). Limpasuvan et al. (2016) found that strong WPWs in particular dissipate in the MLT and enhance downwelling over the pole. The presence of these waves further indicates the strong connection between the stratosphere and the overlying atmosphere.
Based on the composite of 13 ES-SSW events, Limpasuvan et al. (2016) reported a robust signature of eastward-propagating planetary waves (EPWs) that intensify approximately a week before ES-SSW onset (see their Fig. 10). With roughly a 10-day period, these waves appear over the polar lower mesosphere region. Focusing on the 2009 ES-SSW event, past studies have shed new light on EPWs. Specifically, Iida et al. (2014) attributed the observed EPWs to local manifestations of barotropic/baroclinic instability in the lower mesosphere. EPW production relieved baroclinic instability and induced an eastward acceleration on the background flow, counteracting westward accelerations from dissipating quasi-stationary planetary waves (QSPWs) (Iwao and Hirooka 2021). Rhodes et al. (2021) suggested that these EPWs originate from the overreflection of upward-propagating EPWs from the lower stratosphere as they approach a region of strong wind shear (e.g., Harnik and Heifetz 2007). The overreflection process then produces EPWs that emanate from the unstable region. Alternatively, Song et al. (2020) indicated that asymmetric gravity wave drag (GWD) in the lower mesosphere can locally generate EPWs that propagate downward into the lower stratosphere.
The new stratopause reformation following ES-SSW onset is accompanied by the slow westward-propagating planetary waves (WPWs). These WPWs are attributed mainly to barotropic/baroclinic instability of the westward polar stratospheric wind, fostered by the anomalous polar warming (Chandran et al. 2013b; Limpasuvan et al. 2012; Tomikawa et al. 2012). With periods between 5 and 12 days, the generated WPWs can propagate into the MLT. Their damping causes strong westward forcing above 80 km that can drive a strong polar downwelling and initiate the intense downward transport of long-lived tracers into the stratosphere (Orsolini et al. 2010). Forcing due to WPWs may help promote the MLT’s recovery from ES-SSW and the conditions for the stratopause reformation at an elevated altitude (Limpasuvan et al. 2016).
The exact nature of EPWs and WPWs and their underlying mechanisms with respect to the ES-SSW phase (i.e., before or after onset) remain uncertain. This uncertainty stems from the small number of observed ES-SSW events that hinders our ability to develop a robust picture of these traveling waves. With the advent of satellite observations since 1980, fewer than 30 ES-SSWs had been identified. The aforementioned studies on these waves were largely based on case studies (particularly, the strong 2009 ES-SSW event) or composites of few events (Limpasuvan et al. 2016). As such, traveling waves are typically overlooked in the context of SSWs. However, having more knowledge of these waves and their roles during ES-SSWs may provide new insights into ES-SSWs, which are recognized as playing a major factor in the surface climate.
This study aims to better understand the sources and impacts of PWs of various phase speeds surrounding ES-SSW events. Our objectives are to 1) characterize the background flow, the associated wave structure, and their coevolution and 2) identify plausible sources that lead to the wave appearance. Our study leverages a unique ensemble numerical experimental setup that yields many ES-SSW events to develop a robust picture of EPWs and WPWs with respect to winters without ES-SSWs. With this framework, we evaluate EPWs, QSPWs, and WPWs during three different time periods: before, 0–10 days after, and beyond 10 days after SSW onset. Our results indicate that EPWs generated by GW dissipation propagate within the region of westward stratospheric wind. EPW growth in this region applies an eastward acceleration, acting to restore the eastward stratospheric winds. Our study explores how QSPWs and WPWs interact with the westward stratospheric wind and reaffirms that WPWs are produced from instabilities at the top of the westward stratospheric wind reversal.
2. Background and methods
The Whole Atmosphere Community Climate Model (WACCM), version 4, developed at the National Center for Atmospheric Research (NCAR) is an atmosphere-only global chemistry–climate model extending up to ∼145 km (Marsh et al. 2013). WACCM was run in the specified dynamics configuration (WACCM-SD) with a horizontal resolution of 0.95° latitude × 1.25° longitude, 88 vertical levels, and key dynamical variables output daily. In this configuration, the simulated temperature and dynamics were constrained up to 50 km with 6-hourly Modern-Era Retrospective Analysis for Research and Application (MERRA), version 2, reanalysis (Gelaro et al. 2017). Temperature and wind fields are nudged by 10% every 30 min through a mass-conserving interpolation of MERRA reanalysis onto the WACCM-SD horizontal grid. See Orbe et al. (2020) for a tangential discussion on nudging experiments. A linear transition is applied between the nudged output below 50 km and the overlying, fully interactive, free-running region above 60 km. Run from 1980 to 2013, the WACCM-SD simulation constitutes the “base run” from which ensembles were generated.
a. SSW identification and classification in the base model run
The definition of ES-SSW events varies significantly in previous studies such that the stratopause altitude has a discontinuity of 10 km (Limpasuvan et al. 2016), 15 km (Chandran et al. 2013a), and 18 km (Karami et al. 2023). With these criteria, ES-SSWs had frequencies of 5, 3, and 2 decade−1, respectively. We identified ES-SSW events using the criteria from Limpasuvan et al. (2016). The criteria provide rigorous constraints on the upper stratosphere, which will be key in understanding the interaction of PWs in this region. Winters when none of these criteria are met persistently (i.e., lasting longer than 5 days) will be referred to as “normal winters.” For example, a winter containing an SSW without an elevated stratopause would be classified as neither a normal winter nor a winter containing an ES-SSW. For the remainder of the study, ES-SSWs will be identified simply as SSWs.
Four identified SSW events were selected from the base run (see Table 1) as reference cases for our ensemble experiments (described below). They have the SSW onset dates of 12 February 1984, 9 January 2006, 22 January 2009, and 5 January 2013. These onset dates are in close agreement with observations, as expected, since the model is nudged with observations. The first two SSW events are classified as “displaced type,” characterized by the perturbed polar vortex shifting off the North Pole sufficiently to produce an SSW (Charlton and Polvani 2007; Kuttippurath and Nikulin 2012). The latter two events are classified as “split type,” where the separation of the polar vortex into two distinct vortices results in an SSW (Coy and Pawson 2015; Kuttippurath and Nikulin 2012; Manney et al. 2009).
Number of normal and SSW ensemble members generated with respect to reference SSW onset dates.
b. Ensemble setup
The four selected SSW events from the base run correspond to observations below 50 km (via nudging) and are free-running above 60 km. In setting up our ensemble experiment, we retain all the specified dynamics configurations noted above. However, we only nudge below the lowermost model level with a linear transition to 0.4 km, leaving the higher levels as free-running. For a selected SSW event, we initialize each ensemble member by randomly perturbing the temperature field of the base run at 40 days prior to the event’s reference SSW onset date shown in Table 1. The perturbation amount is below the model’s rounding error of ∼10−14 K (e.g., Kay et al. 2015). The 40-day lead time well exceeds the reported SSW predictability around 20 days (Domeisen et al. 2020; Karpechko 2018) and allows for randomized outcomes.
Using the aforementioned identification criteria, this setup produced ensemble members with both SSW winters and normal winters (termed “SSW members” and “normal members,” respectively). Some ensemble members are neither normal nor SSW members and were excluded from the results. For each selected SSW event, at least 10 normal members and 10 SSW members were generated. The number of SSW and normal members collected for each event is shown in Table 1. In total, 76 SSW winters and 68 normal winters were collected.
c. Data analyses
The Eliassen–Palm (EP) flux and its divergence were computed using the formulation associated with the transformed Eulerian-mean (TEM) equations given in Andrews et al. (1987). For these calculations, 5-day running averages were applied to dependent field variables (wind, temperature, etc.) to remove perturbations with periods < 5 days (or with cx > 46 m s−1 at 60°N).
However, the value of n2 is difficult to composite as it varies widely, often approaching infinity. Instead, a critical level can be determined to occur when
d. Composites and anomalies
Composites of the SSW members from all selected SSW events are made after aligning them with respect to their SSW onset date. The composite of the normal members was aligned with respect to the reference SSW onset date listed in Table 1. The alignment date is referred to as day 0. The risk of bias for a particular set of ensemble members (described by a row in Table 1) was eliminated by averaging the ensemble members in each set first and then averaging the sets together. Anomalies are calculated by subtracting the composited diagnostic across normal members from the composited diagnostic across SSW members. Assuming a Gaussian distribution of the diagnostic values at each location and time for SSW members and normal members (weighted such that each ensemble set has an equal contribution), the area where these distributions overlap was calculated. This area is the probability that a value found at a location and reference day during an SSW would also be found during a normal winter. Subtracting the area from unity gives the probability that, given an SSW occurrence, the diagnostic value will be different from that found during normal winters. Alternatively, it is the probability of being abnormal (Ab) given the occurrence of an SSW, which is henceforth abbreviated as P(Ab|SSW). Values of P(Ab|SSW) greater than 0.5 suggest that the anomaly is more likely than not associated with an SSW. These anomalies are good indicators of SSW occurrence. Although an anomaly may be large, a low P(Ab|SSW) suggests that the values are still large in scenarios when no SSW occurs. Thus, a large anomaly at this location and reference day may be associated with SSWs but is not a good indicator for the occurrence of an SSW.
e. The overreflection perspective
An elegant perspective on the interaction of PWs with atmospheric boundaries comes from Lindzen et al. (1980) in which the zero isopleths of the relative PW phase velocity,
Wave geometries of vertical PW propagation between 20 and 90 km through the evolution of the SSW. In our composite, these scenarios occur (a) before day 0, (b) on days 0–25, and (c) on days 0–10. The scenario in (c) is similar to that in (b) but focuses exclusively on PWs generated by asymmetric GWD. Angled arrows emphasize overreflection but do not suggest any change in wave phase speed.
Citation: Journal of the Atmospheric Sciences 80, 11; 10.1175/JAS-D-22-0266.1
Figure 1a shows the wave geometry before SSW onset, which is similar to the winter climatology. Upward-propagating incident PWs (thin solid arrows) can overreflect if a perturbation is able to tunnel from the turning level (thin dashed line) to the critical level (thick solid line). This results in an overreflected wave below the turning level and a transmitted wave above the critical level. Overreflection and transmission are represented by thicker solid arrows and thin dashed arrows, respectively. The overreflected wave carries more energy than the incident wave; this wave energy growth can be indicated by EP flux divergence.
The evanescent region lies beyond the turning level where PW propagation is not possible. Here, a tunneling perturbation exponentially decays in amplitude. A thinner evanescent region would increase the likelihood of a PW being able to reach a boundary beyond the turning level. Additionally, the wave geometry a PW experiences depends on the zonal phase velocity of the wave. Blue and green arrows represent PWs with westward and eastward zonal phase velocities, respectively. As shown in Fig. 1a, westward (eastward) waves would experience a thicker (thinner) evanescent region, making it more difficult (easier) to tunnel and stimulate PW growth at the critical level. Hence, phase velocity affects where a PW propagates and dissipates. While additional nonconservative considerations (such as asymmetric GW drag) are necessary, the overreflection perspective connects incident, transmitted, and overreflected waves to one another.
Notably, our study uses composites with respect to an SSW onset date. During individual SSW cases, overreflection could occur at different relative times or not occur at all. Our study seeks to evaluate the general trends in PW growth and dissipation during an SSW event.
f. Asymmetric gravity wave drag
A multiwave parameterization is implemented in WACCM such that a momentum flux versus phase speed function is represented by a set of discretized GWs. Separate parameterizations are implemented for GWs sourced from orography, convection, and frontogenesis. At the source level, the GW is launched and momentum flux is distributed directly aloft depending on the background winds. Detailed discussion on WACCM’s GW parameterization can be found in appendix A of Garcia et al. (2007) and in Richter et al. (2010). Even for parameterized GWs in WACCM, the background winds will influence where and how much momentum flux is deposited.
The troposphere may not be the only source mechanism for PWs in the middle atmosphere. The preferential filtering of GWs by the underlying stratospheric winds can result in zonally asymmetric GW drag in the MLT (Smith 2003). Lieberman et al. (2013) showed that the wintertime characteristics of MLT perturbations were qualitatively consistent with a simple model of dissipating GWs generating a wavenumber-1 PW after being filtered by underlying stratospheric PW perturbations. Ultimately, the GW-induced ageostrophic winds would result in a divergence (convergence) pattern mirroring the stratospheric high pressure (low pressure) regions below (Lieberman et al. 2013). Resultantly, PW perturbations would be imprinted on the MLT by flow divergence induced by GWs that survived the wind perturbations in the stratosphere.
3. Results
The stratosphere maintains an eastward flow throughout a normal winter and experiences an eastward-to-westward wind reversal during winters with SSWs. The zonal-mean zonal wind during normal winters and winters with SSWs are shown in Fig. 2. For normal winters (Fig. 2a), the mesospheric zero-wind line is maintained at 82.6 km with a standard deviation of 12.4 km based on the altitude of the mesospheric zero-wind line. This was calculated with zonal wind values collected 41 days surrounding the reference SSW onset date (Table 1) of each normal ensemble member. During winters with SSWs (Fig. 2b), the zero-wind line rapidly descends leading to wind reversal at 1 hPa on day 0. The stratopause (thick black dotted line) is maintained around 60 km during normal winters (Fig. 2a). While the altitude of the stratopause and the timing of its postonset reformation are not consistent between ensemble members, the stratopause discontinuity and altitude change can still be depicted in the composite. During winters with SSWs (Fig. 2b), the stratopause descends rapidly upon SSW onset to ∼50 km. Roughly around 15–25 days after onset, the stratopause reforms at ∼60 km, which results in a >10-km vertical discontinuity.
Height–time composites of
Citation: Journal of the Atmospheric Sciences 80, 11; 10.1175/JAS-D-22-0266.1
Ultimately, the stark contrast between a steady zero-wind line during normal winters and a dramatically descending zero-wind line during SSW winters depends on the stability of the mesospheric flow. The stability of the flow is diagnosed by
Being a part of a composite study, the exact position of critical levels for EPWs, QSPWs, and WPWs varies between ensemble members and cannot be explicitly shown. However, Fig. 3 shows the composited phase speed distribution of geopotential height perturbation amplitude
Phase speed vs time composites of
Citation: Journal of the Atmospheric Sciences 80, 11; 10.1175/JAS-D-22-0266.1
In Fig. 2, blue (green) contours emphasize the −8 (5) m s−1 isotachs that separate the wind regimes containing critical levels for WPWs, QSPWs, and EPWs. EPWs would find their critical level outside green contours (winds greater than 5 m s−1), QSPWs between green and blue contours, and WPWs beyond blue contours (winds less than −8 m s−1). However, the ability for PWs to reach their critical level depends on the stability of the zonal-mean flow.
In Fig. 4a, positive (red-shaded contours) and negative (gray-shaded contours)
Height–time (a),(b) composites and (c) anomalies of
Citation: Journal of the Atmospheric Sciences 80, 11; 10.1175/JAS-D-22-0266.1
Prior to SSW events, the background flow in the middle atmosphere only becomes baroclinically unstable (indicated by a dark red-shaded region in Fig. 4b) due to the curvature of wind shear as the wind transitions from eastward in the mesosphere to westward in the lower thermosphere. This baroclinic instability could be exacerbated by the westward wind bias in the lower thermosphere of WACCM, hypothesized to be due to the omission of secondary or higher-order GWs. However, during times of a more perturbed vortex like an SSW period, model winds are closer to observations (Harvey et al. 2022).
Figure 4c shows positive (negative)
The zonal-mean zonal wind is closely tied to
Figure 4 reveals that three distinct wind configurations are present surrounding an SSW, which can be evaluated across three notable periods. Sections 3a–c address periods of no reversed stratospheric wind layer (before day 0), a thick reversed stratospheric wind layer (days 0–10), and a thin reversed stratospheric wind layer (after day 10), respectively. Subordinate sections are differentiated by PW phase speed such that the origins and impacts of WPWs, QSPWs, and EPWs are assessed for each scenario.
a. No reversed stratospheric wind layer—Before day 0
1) WPW and QSPW interaction with no reversed stratospheric wind layer—Before day 0
EP flux indicates the amount of energy transported by a PW and its convergence (divergence) depicts the deposition (absorption) of energy by the PW in the form of heat and momentum. For the quasigeostrophic case, the meridional and vertical components of EP flux are related to momentum and heat, respectively (Andrews et al. 1987). Therefore, EP flux shows the direction of wave propagation and flux convergence (divergence) indicates wave deposition (growth). Figures 5a–c depict the composited EP fluxes for WPWs (Fig. 5a), QSPWs (Fig. 5b), and EPWs (Fig. 5c) such that the color-filled contours indicate vertical flux, whereas the blue (red) contours show EP flux convergence (divergence). Brown-filled contours in Figs. 5a–c indicate that QSPWs have an exponentially greater upward flux than WPWs or EPWs prior to SSW. This is expected since QSPWs are the main driver of SSWs. Figures 5d–f show positive (negative) EP flux divergence anomalies in solid (dashed) contours. Orange-shaded contours show P(Ab|SSW). Five or more days prior to SSW onset, this probability is low, suggesting that EP flux convergence/divergence is not a good indicator of SSW occurrence. From day −5 to day −1, QSPW flux convergence anomalies and P(Ab|SSW) become larger as QSPWs induce SSW onset.
Height–time plots averaged between 60° and 70°N comparing fluxes of (a),(d) WPWs, (b),(e) QSPWs, and (c),(f) EPWs. (left) Upward (downward) vertical EP fluxes are shown in tan (blue) shadings and are outlined by black regular-sized and thin contours incremented by +0.2 × 2i and −0.2 × 2i kg s−2, respectively, where i can be [1, 2, 3, 4, 5, 6, 7, 8]. EP flux convergence (blue contours) and divergence (red contours) are incremented by 2 m s−1 day−1. (right) Positive (negative) anomalies of EP flux divergence are indicated by solid (dashed) contours and incremented by 5 m s−1 day−1. The probability that the anomaly is abnormal is given by orange-shaded contours. The zero-wind line is indicated by a thick bold contour.
Citation: Journal of the Atmospheric Sciences 80, 11; 10.1175/JAS-D-22-0266.1
Figure 6 shows EP flux for WPWs, QSPWs, and EPWs organized by row, on different reference days, organized by column. Similar to Fig. 5, red (blue) contours represent EP flux convergence (divergence) and the zero-wind line is represented by a thick black contour. Depending on the column, relevant isotachs that approximate the critical level are also included.
Height–latitude plots of (a)–(c) WPW, (d)–(f) QSPW, and (g)–(i) EPW EP flux divergence and convergence on (left) day −5, (center) day 0, and (right) day 5, shown in red and blue contours, respectively, and incremented by 5 m s−1 day−1. The meridional EP flux vector component was scaled by (100πaρ)−1 cosϕ and the vertical component by (aρ)−1 cosϕ. Regions of
Citation: Journal of the Atmospheric Sciences 80, 11; 10.1175/JAS-D-22-0266.1
Similar to winters devoid of SSWs in Fig. 2a, Fig. 2b shows that prior to day 0, a thick layer of negative
Resultantly, PWs around 65°N are inhibited from propagating to their critical level aloft, which exists beyond the turning level. However, the −8, 0, and 5 m s−1 isotachs (blue, thick black, and green contours) show that EPWs experience the thinnest evanescent region. Therefore, it has the best chance of tunneling to its critical layer.
2) EPW interaction with no reversed stratospheric wind layer—Before day 0
As during a normal winter, discussed in the section 3 introduction, the wave geometry at 70 km before SSW onset favors the overreflection of EPWs. This is evidenced by comparing Figs. 5a–c; divergence near the stratopause prior to SSW is evident only in the eastward component of the flux. Figure 6g shows the EPW flux from a height–latitude perspective. EP flux divergence (red contour) occurs near a critical layer (bold green line) within a
Approaching onset, Fig. 5c shows that a persistent upward flux does not only sustain a region of EP flux divergence but coincides with an increasing amount of EP flux divergence at 70 km approaching SSW onset. This is consistent with the overreflection mechanism; as the evanescent region becomes thinner (as shown in Fig. 4a), incident waves more easily overreflect. In Fig. 5c, the EPW EP flux divergence increases in tandem with
Figure 5f at 70 km shows small average EP flux divergence anomalies with P(Ab|SSW) < 0.1. While EP flux divergence for EPWs is a robust feature prior to SSW growth, the flux divergence at a specific height and reference day (with respect to SSW onset) is not a good indicator for SSW occurrence. Additionally, the low P(Ab|SSW) for all PWs prior to onset (in Figs. 5d–f) suggests that the strength of the tropospheric forcing may not be the only factor in producing an SSW. The persistence of upward PW propagating (and increasing wave overreflection) may be just as important in initiating an SSW event.
b. Thick reversed stratospheric wind layer—From day 0 to day 10
1) WPW interaction with a thick reversed stratospheric wind layer—From day 0 to day 10
Days 0–10 are marked by descended westward winds in the stratosphere and eastward winds in the MLT, resulting in a layer of reversed stratospheric winds that has a zonal-mean westward flow persisting for over 2 weeks. WPWs with more westward phase velocities than the zonal wind speed would propagate past the reversed stratospheric winds unencumbered, whereas WPWs with slower phase velocities would experience a critical level and interact with the reversed stratospheric winds (cf. blue arrows in Fig. 1b). Since
In Fig. 2b, WPWs with phase velocities between −8 m s−1 (blue line) and −20 m s−1 would experience their critical level past a turning level and could overreflect. At the upper boundary of the reversed stratospheric winds, the turning level generally remains within the reversed stratospheric winds. Around day 4, the −10 m s−1 isotach rests above a turning level at the upper boundary of the reversed stratospheric winds. Therefore, WPWs slower than roughly −10 m s−1 would be prone to overreflect at both the bottom and top boundaries, each time extracting energy from the background flow and producing an overreflected and transmitted PW at each boundary (as suggested in Fig. 1b). The transmitted components would compound and result in an amplified WPW aloft. As a composite, these wave geometries are approximated by the average
Compared to EPWs and QSPWs, WPWs have a relatively strong upward EP flux (brown shading) above 80 km (cf. Figs. 5a–c). An EP flux convergence region exceeding 20 m s−1 day−1 in Fig. 5a (blue contours) at ∼95 km shows that WPW dissipation has a significant impact on the background wind. This abnormal EP flux convergence is unique to SSWs with P(Ab|SSW) > 0.6 (shown in Fig. 5d). The enhanced upward EP flux and EP flux convergence for WPWs at high latitudes above 80 km are also evident in Fig. 6c.
Figure 7 shows
Height–time plot averaged between 60° and 70°N. (a)–(c) Geopotential height perturbation amplitudes
Citation: Journal of the Atmospheric Sciences 80, 11; 10.1175/JAS-D-22-0266.1
2) QSPW interaction with a thick reversed stratospheric wind layer—Day 0 to day 10
The zero-wind line is the critical level for stationary PWs which contain much of the upward flux relative to traveling PWs (cf. Figs. 5a–c). On day −5 (Fig. 6d), the zero-wind line is embedded in a region of negative
As expected, Fig. 7b shows a QSPW
3) EPW interaction with a thick reversed stratospheric wind layer—From day 0 to day 10
Interestingly, flux divergence for EPWs at 75 km continues to increase after SSW onset (Fig. 5c). This continued flux divergence is surprising from the overreflection perspective since an upward-propagating EPW should encounter a critical level before the turning level. As illustrated in Fig. 1b, the presence of EPWs (green arrow) in the reversed stratospheric winds should not be possible by upward-propagating EPWs since they would experience a critical level and dissipate in the lower stratosphere. Both Figs. 2b and 6i show the positioning of the green contour at the bottom of the reversed stratospheric winds well outside of the gray-shaded region. Nevertheless, Fig. 5c indicates that flux divergence occurs at both 50 and 70 km from days 0 to 5. These two layers of flux divergence forming at the top and bottom of the reversed stratospheric winds are a common feature, appearing in most of our ensembles at varying NH latitudes and heights. In Fig. 5f, a flux divergence anomaly exceeding 10 m s−1 day−1 with a low P(Ab|SSW) suggests that, while the feature is present, it does not consistently occur at a specific space or time relative to SSW onset. Figure 7c shows that elevated values of EPW
Figure 8 shows the composited GWD at days −10 and 4 averaged between 0.1 and 0.01 hPa. GWD is then smoothed by a 15° longitude and 10° latitude running mean. While westward GWD dominates on day −10 due to filtering from eastward winds below, eastward GWD at 0° longitude indicates that some eastward-propagating GWs are able to propagate through the stratospheric wavenumber-1
Polar plots of zonal GWD and zonal wind averaged from 0.1 to 0.01 hPa (a) 10 days before and (b) 4 days after SSW onset. Before compositing, the average phase from 10 to 5 hPa of the wavenumber-1
Citation: Journal of the Atmospheric Sciences 80, 11; 10.1175/JAS-D-22-0266.1
This asymmetric GWD averaged from 60° to 70°N (bold latitude lines in Fig. 8) is further examined in Fig. 9 with respect to
Altitude vs relative phase shift of zonal GWD averaged between 60° and 70°N. Before compositing, the average phase from 10 to 5 hPa of the wavenumber-1
Citation: Journal of the Atmospheric Sciences 80, 11; 10.1175/JAS-D-22-0266.1
On day −10 near 75 km (Fig. 9a), westward GWD caps the strong eastward winds. Expectedly, this GWD is not significantly different than during normal winters, illustrated by low probability values in Fig. 9c. Four days after SSW between 40 and 80 km, the stratospheric low pressure system with strong eastward wind (thick black contours in Fig. 9a) is replaced by a stratospheric high pressure system with strong westward winds (thin black contours in Fig. 9b) that constitute the reversed stratospheric wind region. The westward wind favors the transmission of eastward-propagating GWs, which impose an eastward GWD above the reversed stratospheric winds. This eastward GWD has a large wavenumber-1 component with the greatest GWD occurring between a 0° and 100° phase shift from the stratospheric wavenumber-1 PW ridge. Figure 9d shows an anomaly exceeding 60 m s−1 day−1 near 75 km with a P(Ab|SSW) > 0.7. As noted by Siskind et al. (2010) and Limpasuvan et al. (2012), the capping of the reversed stratospheric winds by eastward GWD is a reliable consequence of SSWs.
Hovmöller diagrams of GWD are overlaid with geopotential height, Ζg (Fig. 10a), and zonal wind (Fig. 10b) contours averaged from 0.1 to 0.01 hPa and from 60° to 70°N. Their comparison shows the interaction of GWD with regions of low and high pressure (indicated by regions of low and high Ζg, respectively). As in Fig. 9, diagnostics show wavenumber-0 and wavenumber-1 features and are composited with respect to the stratospheric wavenumber-1 geopotential. From approximately day 0 to day 10, a region of increasingly high Ζg shifts and replaces a region of low Ζg. This coincides with the formation of the reversed stratospheric winds layer associated with the movement of a high pressure system over the pole (seen as the development of stratospheric negative
Time vs relative phase shift of zonal GWD averaged from 60° to 70°N and from 0.1 to 0.01 hPa. Before compositing, the average phase from 10 to 5 hPa of the wavenumber-1
Citation: Journal of the Atmospheric Sciences 80, 11; 10.1175/JAS-D-22-0266.1
GWs propagate through the stratospheric wavenumber-1 Ζg ridge, shown by the phase shift, and dissipate on the westward (eastward) side of the mesospheric wavenumber-1 Ζg ridge (trough) shown in reference to Ζg contours. GWD counteracts the underlying westward winds (shown in Figs. 9b and 10b). This GWD is significant with an anomaly exceeding 80 m s−1 day−1 and P(Ab|SSW) > 0.8 (Fig. 10c). The forcing manifests a wavenumber-1 EPW, shown in Fig. 10a by an eastward shift in the high pressure system over time between day 0 and day 10.
In order for asymmetric GWD to directly generate an EPW, an eastward migration of the asymmetry would be expected. However, this does not occur. Therefore, the phase speed of the generated EPW must be a result of the wave geometry where GWD occurs. In other words, asymmetric GWD seeds a region of unstable flow. The wave geometry that the PWs are born into allows them to extract energy from the flow and grow, taking on a phase speed equal to the background wind. Resultantly, EPW growth can be seen in Fig. 5c near 80 km. After day 10, eastward GWD subsides as a low pressure system indicated by low values of geopotential height reforms in this region (Fig. 10a).
Interestingly, this interaction manifests as an EPW that has a destructive interference with the high pressure system, resulting in an overall loss of
In Fig. 10, eastward GWD forcing is seen around a 50° eastward shift from the stratospheric wavenumber-1 ridge a few days before onset. In Fig. 5c, EPW growth also begins a few days before onset. Therefore, the EPWs generated from instability prior to SSWs may be part of one contiguous feature, ultimately enhancing the stability of the mesosphere.
In Fig. 5c, EPW flux divergence likewise occurs at the bottom of the reversed stratospheric winds which suggests another region of EPW wave growth. The wind structure and
c. Thin reversed stratospheric wind layer—After day 10
1) WPW and QSPW interaction with a thin reversed stratospheric wind layer—After day 10
Although weaker than between days 0 and 10, values of WPW
Before day 10, QSPWs were absorbed by the reversed stratospheric winds. After day 10 (see Fig. 7b), a region of enhanced
After day 10, the flux convergence of QSPWs in the MLT increases around 90 km (Fig. 5b). As the composited winds weaken and reverse, the likelihood of QSPWs propagating past the reversed stratospheric winds and impacting the MLT increases. Figure 5e shows the QSPW flux convergence anomaly in the MLT. The associated low P(Ab|SSW) likely results from large variations in the positioning of the stratopause or in tropospheric forcing after day 10. In Fig. 7e, positive
2) EPW interaction with a thin reversed stratospheric wind layer—After day 10
In the scenario shown in Fig. 1b, EPWs are unable to propagate past the reversed stratospheric winds and are trapped in the troposphere. In Fig. 7f, a negative
As discussed in section 3b(3), eastward acceleration near the upper boundary of the reversed stratospheric winds is aided by the overreflection of EPWs seeded by GWs. After day 10, EP flux divergence associated with eastward GWD diminishes as the upper boundary of the reversed stratospheric winds descends below 60 km. The larger density below 60 km would reduce the amount of drag produced from GW momentum deposition and therefore reduce the effectiveness of GWD as a source for EPW growth. Without EPW overreflection providing additional aid to vortex recovery, the reversed stratospheric winds may remain in the lower stratosphere for prolonged periods of time.
The descended region of negative
4. Conclusions
The overreflection perspective was applied to explain PW behaviors with various wave geometries. During normal winters, an unstable mesosphere inhibits QS PW absorption near the zero-wind line by establishing a turning level below. Approaching SSW, a positive feedback loop is created by persistent PW interaction with a thinning evanescent region that vertically orients stratospheric PWs and increases the likelihood of overreflection. This overreflection tends to produce waves with eastward phase velocities, illustrated by a persistent EP flux divergence around 70 km in Fig. 5c. The resulting EPW growth acts to prevent the descent of the mesospheric westward wind into the stratosphere, counteracting the effects of upward-propagating QSPWs (Iwao and Hirooka 2021; Rhodes et al. 2021). After sufficient thinning of the evanescent region, the critical level becomes exposed to upward-propagating PWs and results in a rapid descent of the zero-wind line.
After day 0, GWD acts as a source mechanism on the upper boundary of the reversed stratospheric winds. EPWs can become trapped resulting in two layers of EP flux divergence at the top and bottom boundaries of the reversed stratospheric winds. To our knowledge, this is the first time this feature has ever been discussed. Additionally, WPWs can tunnel into the reversed stratospheric winds region, overreflect, and dissipate in the MLT. The production of WPWs from instability was also described by Kinoshita et al. (2010), Tomikawa et al. (2012), and Limpasuvan et al. (2016). The reversed stratospheric winds are maintained in the stratosphere due to tropospheric QSPW forcing suggested by the close proximity of the turning level and zero-wind line in the lower portion of the reversed stratospheric winds. The persistence of tropospheric forcing is just as important as the strength of tropospheric forcing in inducing SSWs, as found in reanalyses and forecast models (Orsolini et al. 2018).
The present study shows that SSW recovery has an evolving wave geometry that can support downward and upward vertical EP flux, properties of both reflective and absorptive SSWs, respectively, noted in Kodera et al. (2016). Additionally, the quicker recovery of reflective SSWs with downward vertical EP flux shown by Kodera et al. (2016) is indicative of overreflection; a period of strong overreflection would act to restore the polar vortex by inducing an eastward acceleration.
Additionally, a significant presence of WPWs in the MLT after SSW onset agrees with Limpasuvan et al. (2016), which identifies instability as their source mechanism. Sassi et al. (2016) also identify WPWs in the MLT and show that they can significantly impact the mean meridional circulation, enhancing upwelling in the tropics and downwelling at the pole. Enhanced polar downwelling can result in a descent of more nitric oxides produced from energetic particle precipitation at high latitudes, and case studies have shown that WPWs play a key role in modulating this descent in the MLT (Harvey et al. 2021; Orsolini et al. 2017). Harvey et al. (2021) found in a case study of the January 2009 SSW that the longitudinal asymmetry of the polar vortex impacts meridional circulation such that descent rates are 5 times larger within a PW trough. While further research is needed, our study offers a mechanism by which the reversed stratospheric winds can modulate PW phase velocities present in the MLT during SSW recovery.
The overreflection perspective implements critical layer theory to create a framework in which to evaluate PW interaction with various boundaries in the middle atmosphere. The present study shows that PW wave geometries are sensitive to their phase speeds, resulting in drastically different interactions with regions of varying winds like the reversed stratospheric winds. While our composite study shows the general influences of PWs relative to SSW onset, case studies may offer further insight into wave dynamics affecting SSW development that are not particularly correlated to the onset date.
Acknowledgments.
CTR was supported by funding from the National Science Foundation (NSF) awards (RUI 1642232; REU 1560210). VL was supported by NSF Intergovernmental Panel Agreement. Edits were completed while CTR held a National Research Council Research Associateship award at the U.S. Naval Research Laboratory. The authors acknowledge the computing and technical support from the NCAR Computation Information Systems Laboratory as well as the Coastal Carolina University Cyberinfrastructure (CCU CI) project, funded in part by NSF Award MRI 1624068. The authors also recognize the assistance of Dr. Stephen Eckermann in this study’s publication.
Data availability statement.
The relevant daily model output can be accessed through the CCU CI at https://mirror.coastal.edu/sce/CTR_SSWEnsembles.
REFERENCES
Andrews, D. G., C. B. Leovy, and J. R. Holton, 1987: Middle Atmosphere Dynamics. Academic Press, 487 pp.
Baldwin, M. P., and T. J. Dunkerton, 2001: Stratospheric harbingers of anomalous weather regimes. Science, 294, 581–584, https://doi.org/10.1126/science.1063315.
Bancalá, S., K. Krüger, and M. Giorgetta, 2012: The preconditioning of major sudden stratospheric warmings. J. Geophys. Res., 117, D04101, https://doi.org/10.1029/2011JD016769.
Butler, A. H., D. J. Seidel, S. C. Hardiman, N. Butchart, T. Birner, and A. Match, 2015: Defining sudden stratospheric warmings. Bull. Amer. Meteor. Soc., 96, 1913–1928, https://doi.org/10.1175/BAMS-D-13-00173.1.
Chandran, A., R. L. Collins, R. R. Garcia, D. R. Marsh, V. L. Harvey, J. Yue, and L. de la Torre, 2013a: A climatology of elevated stratopause events in the Whole Atmosphere Community Climate Model. J. Geophys. Res. Atmos., 118, 1234–1246, https://doi.org/10.1002/jgrd.50123.
Chandran, A., R. R. Garcia, R. L. Collins, and L. C. Chang, 2013b: Secondary planetary waves in the middle and upper atmosphere following the stratospheric sudden warming event of January 2012. Geophys. Res. Lett., 40, 1861–1867, https://doi.org/10.1002/grl.50373.
Charlton, A. J., and L. M. Polvani, 2007: A new look at stratospheric sudden warmings. Part I: Climatology and modeling benchmarks. J. Climate, 20, 449–469, https://doi.org/https://doi.org/10.1175/JCLI3996.1.
Coy, L., and S. Pawson, 2015: The major stratospheric sudden warming of January 2013: Analyses and forecasts in the GEOS-5 data assimilation system. Mon. Wea. Rev., 143, 491–510, https://doi.org/10.1175/MWR-D-14-00023.1.
Dickinson, R. E., 1973: Baroclinic instability of an unbounded zonal shear flow in a compressible atmosphere. J. Atmos. Sci., 30, 1520–1527, https://doi.org/10.1175/1520-0469(1973)030<1520:BIOAUZ>2.0.CO;2.
Domeisen, D. I. V., and Coauthors, 2020: The role of the stratosphere in subseasonal to seasonal prediction: 1. Predictability of the stratosphere. J. Geophys. Res. Atmos., 125, e2019JD030920, https://doi.org/10.1029/2019JD030920.
Garcia, R. R., D. R. Marsh, D. E. Kinnison, B. A. Boville, and F. Sassi, 2007: Simulation of secular trends in the middle atmosphere, 1950–2003. J. Geophys. Res., 112, D09301, https://doi.org/10.1029/2006JD007485.
Gelaro, R., and Coauthors, 2017: The Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2). J. Climate, 30, 5419–5454, https://doi.org/10.1175/JCLI-D-16-0758.1.
Gille, J. C., and L. V. Lyjak, 1986: Radiative heating and cooling rates in the middle atmosphere. J. Atmos. Sci., 43, 2215–2229, https://doi.org/10.1175/1520-0469(1986)043<2215:RHACRI>2.0.CO;2.
Harnik, N., and E. Heifetz, 2007: Relating overreflection and wave geometry to the counterpropagating Rossby wave perspective: Toward a deeper mechanistic understanding of shear instability. J. Atmos. Sci., 64, 2238–2261, https://doi.org/10.1175/JAS3944.1.
Harvey, V. L., S. Datta-Barua, N. M. Pedatella, N. Wang, C. E. Randall, D. E. Siskind, and W. E. van Caspel, 2021: Transport of nitric oxide via Lagrangian coherent structures into the top of the polar vortex. J. Geophys. Res. Atmos., 126, e2020JD034523, https://doi.org/10.1029/2020JD034523.
Harvey, V. L., N. Pedatella, E. Becker, and C. Randall, 2022: Evaluation of polar winter mesopause wind in WACCMX+DART. J. Geophys. Res. Atmos., 127, e2022JD037063, https://doi.org/10.1029/2022JD037063.
Iida, C., T. Hirooka, and N. Eguchi, 2014: Circulation changes in the stratosphere and mesosphere during the stratospheric sudden warming event in January 2009. J. Geophys. Res. Atmos., 119, 7104–7115, https://doi.org/10.1002/2013JD021252.
Iwao, K., and T. Hirooka, 2021: Opposite contributions of stationary and traveling planetary waves in the Northern Hemisphere winter middle atmosphere. J. Geophys. Res. Atmos., 126, e2020JD034195, https://doi.org/10.1029/2020JD034195.
Karami, K., S. Borchert, R. Eichinger, C. Jacobi, A. Kuchar, S. Mehrdad, P. Pisoft, and P. Sacha, 2023: The climatology of elevated stratopause events in the UA-ICON model and the contribution of gravity waves. J. Geophys. Res. Atmos., 128, e2022JD037907, https://doi.org/10.1029/2022JD037907.
Karpechko, A. Y., 2018: Predictability of sudden stratospheric warmings in the ECMWF extended-range forecast system. Mon. Wea. Rev., 146, 1063–1075, https://doi.org/10.1175/MWR-D-17-0317.1.
Kay, J. E., and Coauthors, 2015: The Community Earth System Model (CESM) large ensemble project: A community resource for studying climate change in the presence of internal climate variability. Bull. Amer. Meteor. Soc., 96, 1333–1349, https://doi.org/10.1175/BAMS-D-13-00255.1.
Kinoshita, T., Y. Tomikawa, and K. Sato, 2010: On the three-dimensional residual mean circulation and wave activity flux of the primitive equations. J. Meteor. Soc. Japan, 88, 373–394, https://doi.org/10.2151/jmsj.2010-307.
Kodera, K., H. Mukougawa, P. Maury, M. Ueda, and C. Claud, 2016: Absorbing and reflecting sudden stratospheric warming events and their relationship with tropospheric circulation. J. Geophys. Res. Atmos., 121, 80–94, https://doi.org/10.1002/2015JD023359.
Kuttippurath, J., and G. Nikulin, 2012: A comparative study of the major sudden stratospheric warmings in the Arctic winters 2003/2004–2009/2010. Atmos. Chem. Phys., 12, 8115–8129, https://doi.org/10.5194/acp-12-8115-2012.
Lieberman, R. S., D. M. Riggin, and D. E. Siskind, 2013: Stationary waves in the wintertime mesosphere: Evidence for gravity wave filtering by stratospheric planetary waves. J. Geophys. Res. Atmos., 118, 3139–3149, https://doi.org/10.1002/jgrd.50319.
Limpasuvan, V., J. H. Richter, Y. J. Orsolini, F. Stordal, and O.-K. Kvissel, 2012: The roles of planetary and gravity waves during a major stratospheric sudden warming as characterized in WACCM. J. Atmos. Sol.-Terr. Phys., 78–79, 84–98, https://doi.org/10.1016/j.jastp.2011.03.004.
Limpasuvan, V., Y. J. Orsolini, A. Chandran, R. R. Garcia, and A. K. Smith, 2016: On the composite response of the MLT to major sudden stratospheric warming events with elevated stratopause. J. Geophys. Res. Atmos., 121, 4518–4537, https://doi.org/10.1002/2015JD024401.
Lindzen, R. S., B. Farrell, and K.-K. Tung, 1980: The concept of wave overreflection and its application to baroclinic instability. J. Atmos. Sci., 37, 44–63, https://doi.org/10.1175/1520-0469(1980)037<0044:TCOWOA>2.0.CO;2.
Manney, G. L., and Coauthors, 2008: The evolution of the stratopause during the 2006 major warming: Satellite data and assimilated meteorological analyses. J. Geophys. Res., 113, D11115, https://doi.org/10.1029/2007JD009097.
Manney, G. L., and Coauthors, 2009: Aura microwave limb sounder observations of dynamics and transport during the record-breaking 2009 Arctic stratospheric major warming. Geophys. Res. Lett., 36, L12815, https://doi.org/10.1029/2009GL038586.
Marsh, D. R., M. J. Mills, D. E. Kinnison, J.-F. Lamarque, N. Calvo, and L. M. Polvani, 2013: Climate change from 1850 to 2005 simulated in CESM1(WACCM). J. Climate, 26, 7372–7391, https://doi.org/10.1175/JCLI-D-12-00558.1.
Matsuno, T., 1971: A dynamical model of stratospheric sudden warming. J. Atmos. Sci., 28, 1479–1494, https://doi.org/10.1175/1520-0469(1971)028<1479:ADMOTS>2.0.CO;2.
O’Neill, A., and C. E. Youngblut, 1982: Stratospheric warmings diagnosed using the transformed Eulerian-mean equations and the effect of the mean state on wave propagation. J. Atmos. Sci., 39, 1370–1386, https://doi.org/10.1175/1520-0469(1982)039<1370:SWDUTT>2.0.CO;2.
Orbe, C., and Coauthors, 2020: Description and evaluation of the specified-dynamics experiment in the Chemistry-Climate Model Initiative. Atmos. Chem. Phys., 20, 3809–3840, https://doi.org/10.5194/acp-20-3809-2020.
Orsolini, Y. J., J. Urban, D. P. Murtagh, S. Lossow, and V. Limpasuvan, 2010: Descent from the polar mesosphere and anomalously high stratopause observed in 8 years of water vapor and temperature satellite observations by the Odin Sub-Millimeter Radiometer. J. Geophys. Res., 115, D12305, https://doi.org/10.1029/2009JD013501.
Orsolini, Y. J., V. Limpasuvan, K. Pérot, P. Espy, R. Hibbins, S. Lossow, K. Raaholt Larsson, and D. Murtagh, 2017: Modelling the descent of nitric oxide during the elevated stratopause event of January 2013. J. Atmos. Sol.-Terr. Phys., 155, 50–61, https://doi.org/10.1016/j.jastp.2017.01.006.
Orsolini, Y. J., K. Nishii, and H. Nakamura, 2018: Duration and decay of Arctic stratospheric vortex events in the ECMWF seasonal forecast model. Quart. J. Roy. Meteor. Soc., 144, 2876–2888, https://doi.org/10.1002/qj.3417.
Orsolini, Y. J., J. Zhang, and V. Limpasuvan, 2022: Abrupt change in the lower thermospheric mean meridional circulation during sudden stratospheric warmings and its impact on trace species. J. Geophys. Res. Atmos., 127, e2022JD037050, https://doi.org/10.1029/2022JD037050.
Rhodes, C. T., V. Limpasuvan, and Y. J. Orsolini, 2021: Eastward-propagating planetary waves prior to the January 2009 sudden stratospheric warming. J. Geophys. Res. Atmos., 126, e2020JD033696, https://doi.org/10.1029/2020JD033696.
Richter, J. H., F. Sassi, and R. R. Garcia, 2010: Towards a physically based gravity wave source parameterization in a general circulation model. J. Atmos. Sci., 67, 136–156, https://doi.org/10.1175/2009JAS3112.1.
Sassi, F., H.-L. Liu, and J. T. Emmert, 2016: Traveling planetary-scale waves in the lower thermosphere: Effects on neutral density and composition during solar minimum conditions. J. Geophys. Res. Space Phys., 121, 1780–1801, https://doi.org/10.1002/2015JA022082.
Siskind, D. E., S. D. Eckermann, J. P. McCormack, L. Coy, K. W. Hoppel, and N. L. Baker, 2010: Case studies of the mesospheric response to recent minor, major, and extended stratospheric warmings. J. Geophys. Res., 115, D00N03, https://doi.org/10.1029/2010JD014114.
Smith, A. K., 2003: The origin of stationary planetary waves in the upper mesosphere. J. Atmos. Sci., 60, 3033–3041, https://doi.org/10.1175/1520-0469(2003)060<3033:TOOSPW>2.0.CO;2.
Song, B.-G., H.-Y. Chun, and I.-S. Song, 2020: Role of gravity waves in a vortex-split sudden stratospheric warming in January 2009. J. Atmos. Sci., 77, 3321–3342, https://doi.org/10.1175/JAS-D-20-0039.1.
Tomikawa, Y., K. Sato, S. Watanabe, Y. Kawatani, K. Miyazaki, and M. Takahashi, 2012: Growth of planetary waves and the formation of an elevated stratopause after a major stratospheric sudden warming in a T213L256 GCM. J. Geophys. Res., 117, D16101, https://doi.org/10.1029/2011JD017243.