Volume 47, Issue 1 ee2019GL086577
Research Letter
Open Access

September 2019 Antarctic Sudden Stratospheric Warming: Quasi‐6‐Day Wave Burst and Ionospheric Effects

Y. Yamazaki

Corresponding Author

GFZ German Research Centre for Geosciences, Potsdam, Germany

Correspondence to: Y. Yamazaki,

yamazaki@gfz-potsdam.de

Search for more papers by this author
V. Matthias

Potsdam Institute for Climate Impact Research, Potsdam, Germany

Search for more papers by this author
Y. Miyoshi

Department of Earth and Planetary Sciences, Kyushu University, Fukuoka, Japan

Search for more papers by this author
C. Stolle

GFZ German Research Centre for Geosciences, Potsdam, Germany

Faculty of Science, University of Potsdam, Potsdam, Germany

Search for more papers by this author
T. Siddiqui

GFZ German Research Centre for Geosciences, Potsdam, Germany

Search for more papers by this author
G. Kervalishvili

GFZ German Research Centre for Geosciences, Potsdam, Germany

Search for more papers by this author
J. Laštovička

Institute of Atmospheric Physics CAS, Prague, Czech Republic

Search for more papers by this author
M. Kozubek

Institute of Atmospheric Physics CAS, Prague, Czech Republic

Search for more papers by this author
W. Ward

Department of Physics, University of New Brunswick, Fredericton, New Brunswick, Canada

Search for more papers by this author
D. R. Themens

Department of Physics, University of New Brunswick, Fredericton, New Brunswick, Canada

Search for more papers by this author
S. Kristoffersen

Department of Physics, University of New Brunswick, Fredericton, New Brunswick, Canada

Search for more papers by this author
P. Alken

Cooperative Institute for Research in Environmental Sciences, University of Colorado Boulder, Boulder, CO, USA

National Centers for Environmental Information, NOAA, Boulder, CO, USA

Search for more papers by this author
First published: 08 January 2020
Citations: 15

Abstract

An exceptionally strong stationary planetary wave with Zonal Wavenumber 1 led to a sudden stratospheric warming (SSW) in the Southern Hemisphere in September 2019. Ionospheric data from European Space Agency's Swarm satellite constellation mission show prominent 6‐day variations in the dayside low‐latitude region at this time, which can be attributed to forcing from the middle atmosphere by the Rossby normal mode “quasi‐6‐day wave” (Q6DW). Geopotential height measurements by the Microwave Limb Sounder aboard National Aeronautics and Space Administration's Aura satellite reveal a burst of global Q6DW activity in the mesosphere and lower thermosphere during the SSW, which is one of the strongest in the record. The Q6DW is apparently generated in the polar stratosphere at 30–40 km, where the atmosphere is unstable due to strong vertical wind shear connected with planetary wave breaking. These results suggest that an Antarctic SSW can lead to ionospheric variability through wave forcing from the middle atmosphere.

Plain Language Summary

A sudden stratospheric warming (SSW) is an extreme wintertime polar meteorological phenomenon occurring mostly over the Arctic region. Studies have shown that Arctic SSW can influence the entire atmosphere. In September 2019, a rare SSW event occurred in the Antarctic region, providing an opportunity to investigate its broader impact on the whole atmosphere. We present observations from the middle atmosphere and ionosphere during this event, noting unusually strong wave activity throughout this region. Our results suggest that an Antarctic SSW can have a significant impact on the whole atmosphere system similar to those due to Arctic events.

1 Introduction

A sudden stratospheric warming (SSW) is a large‐scale meteorological phenomenon in the winter stratosphere, which involves a rapid rise in the polar temperature by a few tens of kelvins in several days (Andrews et al., 1987; Labitzke & Van Loon, 1999). An SSW is triggered by an injection of stationary planetary waves (PWs) from the troposphere, which are driven by topography and land‐sea temperature contrasts. PW breaking in the middle atmosphere leads to an acceleration of the zonal mean flow and changes the mean meridional circulation (Matsuno, 1971). Dynamical effects of PW breaking during SSWs are not limited in the stratosphere but are also well extended into the mesosphere and lower thermosphere (Chandran et al., 2014).

According to the definition by the World Meteorological Organization (McInturff, 1978), a “minor” SSW occurs when a large temperature increase is observed in the winter polar stratosphere, at least by 25 K in a week or less. The event is called “major” if the reversal of the zonal mean flow from eastward to westward occurs poleward of 60° latitude at 10 hPa (32 km) or below, along with the reversal of the meridional temperature gradient. The average number of major SSWs is urn:x-wiley:grl:media:grl60100:grl60100-math-00010.6 per winter in the Northern Hemisphere (NH) (Butler et al., 2015; Charlton & Polvani, 2007). In the Southern Hemisphere (SH), the occurrence of an SSW, whether major or minor, is not as frequent as in the NH because of weaker PW forcing due to smaller topographical differences and land‐sea contrasts. In fact, the September 2002 event (Baldwin et al., 2003; Krüger et al., 2005) is the only major SSW observed in the Antarctic.

In the last decade, the aeronomy community has come to the realization that SSWs can be a significant source of ionospheric variability (Chau et al., 2012; Pedatella et al., 2018). In particular, the January 2009 major Arctic SSW, which took place under extremely quiet solar and geomagnetic activity conditions, enabled many studies to attribute observed ionospheric perturbations to the SSW (e.g., Chau et al., 2010; Goncharenko et al., 2010, 2010; Fejer et al., 2010; Liu et al., 2011; Lin et al., 2019; Nayak & Yiğit, 2019; Oyama et al., 2014; Pedatella & Forbes, 2010; Patra et al., 2014; Pancheva & Mukhtarov, 2011; Rodrigues et al., 2011; Yue et al., 2010; Yadav et al., 2017). Most studies concentrated on the dayside low‐latitude region, where the ionospheric response to the SSW was most pronounced. Modeling studies have suggested that atmospheric tides played an important role in driving ionospheric variability during the January 2009 SSW (Fang et al., 2012; Fuller‐Rowell et al., 2011; Jin et al., 2012; Pedatella et al., 2014; Pedatella & Maute, 2015; Sassi et al., 2013; Wang et al., 2014). Tidal waves at altitudes of the ionospheric urn:x-wiley:grl:media:grl60100:grl60100-math-0002 region (95–150 km) are, in large part, from the middle atmosphere, and their amplitudes and phases can change in response to an SSW (Stening et al., 1997). Among different tidal modes, the semidiurnal lunar tide shows a particularly strong and consistent response to SSWs (Pedatella & Maute, 2015; Zhang & Forbes, 2014). Forbes and Zhang (2012) argued that the large semidiurnal lunar tide observed during the January 2009 SSW can arise from resonant amplification associated with the atmospheric Pekeris mode. Enhanced lunar tidal perturbations in the equatorial ionosphere have been reported for a number of SSW events (Fejer et al., 2011; Liu et al., 2019; Siddiqui et al., 2015, 2018; Park et al., 2012).

As mentioned earlier, SSWs rarely occur in the SH, and the ionospheric response to Antarctic SSWs has been largely unexplored. The only exception is the study by Olson et al. (2013), which examined ionospheric variability during the September 2002 major Antarctic SSW. Although Olson et al. (2013) observed multiday variations in the equatorial ionosphere, their association with the SSW remained somewhat uncertain because of high geomagnetic activity during the event. The main objective of this study is to present observations from the ionosphere and middle atmosphere during the recent Antarctic SSW event in September 2019 and note the presence of unusually strong traveling PW activity throughout the atmosphere and ionosphere at this time.

2 Results and Discussion

2.1 September 2019 SSW

Figure 1 gives an overview of the September 2019 SSW. The polar temperature at 10 hPa, obtained from the MERRA‐2 reanalysis (Gelaro et al., 2017), shows a rapid increase from 207.7 K on 5 September to 258.5 K on 11 September 2019 ( urn:x-wiley:grl:media:grl60100:grl60100-math-0003 = 50.8 K/week) (Figure 1a). This is the largest increase in the Antarctic polar temperature per week in the entire MERRA‐2 data set starting from January 1980. The maximum temperature rise during the September 2002 major SSW was urn:x-wiley:grl:media:grl60100:grl60100-math-0004 = 38.5 K/week. Figure 1b presents the vertical structure of the zonal mean zonal wind at 60°S, as derived from geopotential height measurements by the Aura Microwave Limb Sounder (MLS) (Schwartz et al., 2008; Waters et al., 2006). It can be seen that the eastward zonal mean wind first reversed in the upper mesosphere on 2 September 2019, and in the subsequent days, the region of the wind reversal descended to lower layers, reaching 40 km on 18 September 2019. Since the wind reversal did not reach the 10‐hPa level ( urn:x-wiley:grl:media:grl60100:grl60100-math-000532 km), the event is categorized as a minor warming. Figure 1c shows that there was an enhancement in the amplitude of the stationary PW with Zonal Wavenumber (ZW) 1 during 14–20 August 2019 and during 28 August to 5 September 2019. In both cases, the amplitude attained the largest recorded by Aura/MLS since August 2004. The former event can contribute to the SSW by weakening the zonal mean flow, which is often referred to as preconditioning (e.g., Cámara et al., 2017; Limpasuvan et al., 2004; McIntyre, 1982). Forcing due to PW breaking during the latter event is the likely cause of the zonal wind reversal in the middle atmosphere and hence the SSW. No similar enhancement is found in the amplitude of the stationary PW with ZW2.

grl60100-fig-0001
Overview of middle atmosphere dynamics during the September 2019 sudden stratospheric warming. (a) Stratospheric polar temperature at 10 hPa obtained from the MERRA‐2 reanalysis. The thick black line represents the data for 2019, while the thin purple lines correspond to the data for other years during 1980–2018, among which the data for 2002 are highlighted by green for the occurrence of a major SSW. The red dashed line shows the climatological mean. (b) Zonal mean zonal wind at 60°S derived from the geopotential height (GPH) measurements by the Aura Microwave Limb Sounder (MLS) using the method described by Matthias and Ern (2018). (c) Amplitude of the planetary wave with Zonal Wavenumber (ZW) 1 and ZW2 at 65°S and 48‐km altitude from the Aura/MLS GPH. The red and blue solid lines represent ZW1 and ZW2 waves, respectively. The climatological amplitudes of the ZW1 and ZW2 waves are indicated by the dashed lines with corresponding colors. The gray shaded area show the range between the maximum and minimum values of the amplitude of the ZW1 wave observed by Aura/MLS since August 2004.

As a brief summary, the September 2019 Antarctic SSW was a minor warming but it involved an exceptionally strong stationary PW with ZW1 and a large temperature rise. Furthermore, the event took place during the minimum phase of the solar cycle, similar to the January 2009 SSW, and as will be shown later, overall solar and geomagnetic activities were low, which helps identify SSW effects on the ionosphere. Therefore, the September 2019 event provides an excellent (and rare) opportunity to investigate the ionospheric response to an Antarctic SSW, which is not well understood from previous studies.

2.2 Ionospheric Observations by Swarm

European Space Agency's Earth observation mission Swarm (Friis‐Christensen et al., 2006) involves three identical satellites (A, B, and C), equipped with scientific instruments that are suitable for investigating Earth's magnetic field and its source currents (Friis‐Christensen et al., 2008). The three spacecraft were launched into polar orbits on 22 November 2013, and since 17 April 2014, Swarm A and C fly side by side at an altitude of urn:x-wiley:grl:media:grl60100:grl60100-math-0006460 km, while Swarm B flies at urn:x-wiley:grl:media:grl60100:grl60100-math-0007510 km. Figures 2a–2c show the temporal variability of the equatorial electrojet (EEJ) intensity (e.g., Alken et al., 2015), electron density (e.g., Buchert et al., 2015), and total electron content (TEC) (e.g., Park et al., 2017) as observed by Swarm B during 5 September to 5 October 2019. The data used here were collected from the descending parts of the orbit in 11:00–14:00 magnetic local time (MLT) (see also Figure 2g). Figures 2h and 2i show that overall solar and geomagnetic activity levels were low during this time interval, which is typical for solar minimum conditions. Moderately high geomagnetic activity was observed during 27 September to 1 October 2019, which needs to be taken into account when the ionospheric data are interpreted. Unlike the September 2002 Antarctic SSW, which was examined by Olson et al. (2013), severe geomagnetic activity with urn:x-wiley:grl:media:grl60100:grl60100-math-0008 was not observed. The low urn:x-wiley:grl:media:grl60100:grl60100-math-0009 conditions are preferable for the study of SSW effects on the ionosphere. Modeling studies have shown that the ionospheric response to lower atmospheric forcing would be more pronounced under lower solar flux conditions (Fang et al., 2014; Liu & Richmond, 2013).

grl60100-fig-0002
Overview of ionospheric variations during 5 September to 5 October 2019. (a) Longitude versus time plot of the equatorial electrojet (EEJ) intensity derived from magnetic field measurements in the descending orbits of Swarm B. The data are smoothed using a 3‐day and 50° longitude window. (b) Same as (a) except for the electron density at 20° magnetic latitude. (c) Same as (a) except for the total electron content (TEC) at 20° magnetic latitude at the satellite altitude of urn:x-wiley:grl:media:grl60100:grl60100-math-0010510 km. (d) Percent changes in the Swarm B zonal equatorial electric field (EEF) at urn:x-wiley:grl:media:grl60100:grl60100-math-001190° longitudes with respect to the 29‐day running mean at the corresponding longitudes. (e) Same as (d) except for the Swarm B electron density at 20° magnetic latitude. (f) Same as (d) except for the Swarm B TEC at 20° magnetic latitude. (g) Magnetic local time (MLT) at equatorial crossings for the descending orbits of Swarm B. (h) Geomagnetic activity index urn:x-wiley:grl:media:grl60100:grl60100-math-0012. (i) Solar activity index urn:x-wiley:grl:media:grl60100:grl60100-math-0013.

The EEJ is a narrow band of a zonal electric current that flows along the magnetic equator in the dayside urn:x-wiley:grl:media:grl60100:grl60100-math-0014 region ionosphere at 100‐ to 115‐km altitude (e.g., Yamazaki & Maute, 2017). Under solar minimum conditions with low geomagnetic activity and nearly constant solar flux, day‐to‐day variations of the EEJ intensity are dominated by the changes in neutral winds at urn:x-wiley:grl:media:grl60100:grl60100-math-0015 region heights associated with atmospheric waves from the lower layers (Yamazaki et al., 2014) and thus are a good indicator of lower atmospheric influence on the urn:x-wiley:grl:media:grl60100:grl60100-math-0016 region ionosphere. The methods for deriving the EEJ intensity and equatorial zonal electric field (EEF) from Swarm magnetic field measurements are detailed in Alken et al. (2013). Figure 2a reveals that the EEJ variability was dominated by 6‐day variations during this period. The westward phase propagation of the EEJ intensity perturbations with ZW1 can also be seen. The PW spectrum of the EEJ intensity shows the predominance of the westward propagating ZW1 wave with period 6.0  urn:x-wiley:grl:media:grl60100:grl60100-math-0017 0.2 days (see Figure S1 in the supporting information). Similar spatial and temporal variability was found in the equatorial zonal electric field. Figure 2d shows relative changes in the EEF from the 29‐day running mean, calculated separately at urn:x-wiley:grl:media:grl60100:grl60100-math-001890° longitudes. It can be seen that the EEF underwent 6‐day variations of urn:x-wiley:grl:media:grl60100:grl60100-math-001940% that are out of phase for a 180° longitudinal separation. The amplitude varies in the range of 20–70% depending on the longitude. In a recent study, Yamazaki et al. (2018) reported that the EEJ intensity occasionally shows urn:x-wiley:grl:media:grl60100:grl60100-math-00206‐day variations that have characteristics of a westward propagating wave with ZW1. They attributed the EEJ variations to the quasi‐6‐day wave (Q6DW) that was simultaneously observed in the lower thermosphere. The behavior of the EEJ presented in Figure 2a is similar to those reported by Yamazaki et al. (2018). The amplitude of the 6‐day variations in the EEJ, 10–30 mA/m, is also comparable with those in Yamazaki et al. (2018).

The Q6DW is a westward propagating PW with ZW1, which is occasionally observed in the middle atmosphere (e.g., Forbes & Zhang, 2017; Hirota & Hirooka, 1984; Pancheva et al., 2018; Riggin et al., 2006; Talaat et al., 2001; 2002; Wu et al., 1994). It is often regarded as the (1,1) Rossby normal mode, which is predicted by classical atmospheric wave theory (Forbes, 1995; Kasahara, 1976; Madden, 1979, 2007; Salby, 1984), for its ZW, phase speed, and latitudinal structure. The Q6DW can be excited in the troposphere by heating due to moist convection (Miyoshi & Hirooka, 1999). Additionally, the wave can be excited/amplified in the middle atmosphere due to baroclinic/barotropic instability (Lieberman et al., 2003; Liu et al., 2004; Meyer & Forbes, 1997). Zonal wind perturbations of the Q6DW are largest around the equator and can be up to a few tens of meters per second at urn:x-wiley:grl:media:grl60100:grl60100-math-0021 region heights, which is sufficient to cause detectable changes in dayside ionospheric electric fields and currents (Gan et al., 2016; Miyoshi, 1999; Pedatella et al., 2012). These electric field perturbations in the urn:x-wiley:grl:media:grl60100:grl60100-math-0022 region ionosphere are transmitted to the urn:x-wiley:grl:media:grl60100:grl60100-math-0023 region along equipotential magnetic field lines and affect the distribution of low‐latitude urn:x-wiley:grl:media:grl60100:grl60100-math-0024 region plasmas by modulating their urn:x-wiley:grl:media:grl60100:grl60100-math-0025 plasma drift motions. In this way, the Q6DW can affect the urn:x-wiley:grl:media:grl60100:grl60100-math-0026 region plasma density, as first revealed in the 1990s by ionosonde measurements (e.g., Apostolov et al., 1994; Altadill & Laštovička, 1996; Laštovička, 2006). More recent studies based on global TEC maps have established that the Q6DW effect on the plasma density is largest in the afternoon local time sector near the equatorial ionization anomaly crests ( urn:x-wiley:grl:media:grl60100:grl60100-math-002720° magnetic latitudes) (Gu et al., 2014, 2018; Qin et al., 2019; Yamazaki, 2018).

The 6‐day variations can be seen in both electron density (Figures 2b and 2e) and topside TEC (Figures 2c and 2f) at 20° magnetic latitude. (Figure S2 in the supporting information shows the electron density variations at various latitudes.) The variations are consistent with those in the EEJ/EEF (Figures 2a and 2d), indicating electrodynamic coupling between the urn:x-wiley:grl:media:grl60100:grl60100-math-0028 and urn:x-wiley:grl:media:grl60100:grl60100-math-0029 region ionosphere. The response time of the urn:x-wiley:grl:media:grl60100:grl60100-math-0030 region plasma density to a change in the urn:x-wiley:grl:media:grl60100:grl60100-math-0031 region electric field is 2–4 hr (e.g., Stolle et al., 2008; Venkatesh et al., 2015), which would not be visible in the figures. The relative change in the electron density is in the range of 20–40%, which is appreciably larger than that of TEC, 5–10%. This is not surprising as the amplitude of the Q6DW decreases with altitude in the topside ionosphere, as demonstrated by Gu et al. (2018).

The plasma density and TEC data from the ascending parts of the Swarm B orbit (23:00–02:00 MLT) were also examined, but the 6‐day variations were not as evident as the results derived from the descending orbits. Similarly, the ionospheric data (EEJ, electron density, and TEC) from Swarm A, which was flying around 02:00–05:00 MLT (descending orbits) and 14:00–17:00 MLT (ascending orbits), did not show strong 6‐day variations. The electron density variations from Swarm B (ascending) and Swarm A (ascending and descending) are presented in the supporting information (Figure S3). The different behavior of 6‐day variations in different Swarm data sets reflects the fact that the ionospheric response to the Q6DW depends on MLT and height, as well as on magnetic latitude (Gu et al., 2018). Further studies are required to determine the three‐dimensional structure of the 6‐day ionospheric variations during this event, which could clarify the cause of the differences observed between Swarm A and B.

Previous studies found a significant contribution of the semidiurnal lunar tide to ionospheric variability during NH SSWs (e.g., Park et al., 2012), but it is not known whether the lunar tide plays an equally important role during SH SSWs. The semidiurnal lunar variations in the EEJ intensity derived from the Swarm A and B data during 5 September to 5 October 2019 are presented in the supporting information (Figure S4). It is found that the amplitude of the EEJ semidiurnal lunar variation is 17.7  urn:x-wiley:grl:media:grl60100:grl60100-math-0032 2.1 mA/m for Swarm A (14:00–17:00 MLT) and 16.6  urn:x-wiley:grl:media:grl60100:grl60100-math-0033 2.8 mA/m for Swarm B (11:00–14:00 MLT), which is greater than the climatological value of 9.0  urn:x-wiley:grl:media:grl60100:grl60100-math-0034 0.4 mA/m as reported by Yamazaki et al. (2017) for September daytime (08:00–16:00 local solar time) conditions. The phase, which is defined as the lunar time of maximum, is 10.2  urn:x-wiley:grl:media:grl60100:grl60100-math-0035 0.2 hr for Swarm A and 10.0  urn:x-wiley:grl:media:grl60100:grl60100-math-0036 0.4 hr for Swarm B, which is in good agreement with the climatological value of 10.0  urn:x-wiley:grl:media:grl60100:grl60100-math-0037 0.1 hr. Despite the significant enhancement, the lunar variation accounts for only a small part of the observed EEJ variability (cf. Figures 2a and S4). The relative amplitude of the semidiurnal lunar variation in the topside electron density is 9.9  urn:x-wiley:grl:media:grl60100:grl60100-math-0038 0.7% for Swarm A and 11.1  urn:x-wiley:grl:media:grl60100:grl60100-math-0039 0.1% for Swarm B (also shown in Figure S4). Again, these variations are smaller than the 6‐day variations observed during the same period (Figure 2e).

It is noted that since Swarm slowly precesses in local solar time, it is not possible to resolve short‐term variability of solar tides. Changes in upward‐propagating solar tides can occur during SSWs due to changes in the zonal mean atmosphere (Jin et al., 2012; Pedatella & Liu, 2013), tidal sources (Goncharenko et al., 2012; Siddiqui et al., 2019), and tidal interaction with PWs (Chang et al., 2009; Liu et al., 2010; Maute et al., 2014). Possible changes in solar tides during the September 2019 SSW remain to be investigated.

2.3 Q6DW in the Middle Atmosphere

Traveling PWs in the middle atmosphere are examined using the geopotential height data from Aura/MLS. The analysis method was described in detail in the previous work (Yamazaki & Matthias, 2019) and thus is only briefly summarized here. The amplitude urn:x-wiley:grl:media:grl60100:grl60100-math-0040 and phase urn:x-wiley:grl:media:grl60100:grl60100-math-0041 of waves with period urn:x-wiley:grl:media:grl60100:grl60100-math-0042 were derived by fitting the following formula to the data at a given latitude and height:
urn:x-wiley:grl:media:grl60100:grl60100-math-0043(1)
where urn:x-wiley:grl:media:grl60100:grl60100-math-0044 is the universal time, urn:x-wiley:grl:media:grl60100:grl60100-math-0045 is the longitude, and urn:x-wiley:grl:media:grl60100:grl60100-math-0046 is the ZW. Eastward and westward propagating waves correspond to urn:x-wiley:grl:media:grl60100:grl60100-math-0047 and urn:x-wiley:grl:media:grl60100:grl60100-math-0048, respectively. The data were analyzed for each day using a time window that is 3 times the wave period. The 1‐ urn:x-wiley:grl:media:grl60100:grl60100-math-0049 error in the amplitude is typically below 0.05 km.
grl60100-fig-0003
Overview of quasi‐6‐day wave (Q6DW) activity during the September 2019 SSW as derived from the geopotential height (GPH) measurements by the Aura Microwave Limb Sounder (MLS). (a) Amplitude of the westward propagating Zonal Wavenumber (ZW) 1 waves at 45°S and 97‐km altitude. (b) Same as (a) except for the eastward propagating ZW1 waves. (c) Latitude versus height structures of the westward propagating ZW1 wave with period 6.0 days during 10–30 September 2019. The contour lines indicate the amplitude, while the color represents the phase. (d) Amplitude of the Q6DW, defined here as the maximum amplitude of the westward propagating ZW1 waves at periods 5–7 days, at 45°S and 97‐km altitude. (e) Same as (d) except at 80°S and 32‐km altitude. (f) Vertical structure of the Q6DW at 80°S during 10–30 September 2019. (g) Same as (f) except at 45°S. (h–j) Latitude versus height structures of the zonal mean zonal wind. The areas highlighted by the light‐yellow color indicate the regions where the meridional gradient of the quasi‐geostrophic potential vorticity is negative, which is the necessary condition for barotropic/baroclinic instability.

Figures 3a and 3b show the amplitudes for the westward and eastward propagating waves with ZW1 at 45°S in the lower thermosphere at urn:x-wiley:grl:media:grl60100:grl60100-math-005097 km. Enhanced wave activity can be seen in the westward propagating component (Figure 3a) with period 4–7 days during September 2019, which can be identified as the Q6DW. It is consistent with the appearance of 6‐day variations in the ionosphere (Figures 2a–2c). Such enhanced wave activity is not present in the eastward propagating ZW1 component (Figure 3b) or other components with higher ZWs (not shown here). Although studies have found that the amplitude of the Q6DW in the middle atmosphere is greatest during equinoctial months (Forbes & Zhang, 2017; Qin et al., 2019; Yamazaki, 2018), the wave enhancement in September 2019 was exceptional, with the maximum amplitude larger than 0.4 km in the lower thermosphere, which is much larger than the climatological amplitude (0.15 km) or amplitudes recorded during other individual years during 2004–2018 (Figure 3d). Thus, the large‐amplitude Q6DW observed in September 2019 cannot be explained merely as a seasonal effect.

The latitude and height structures of the 6‐day wave during 10–30 September 2019 are presented in Figure 3c. The amplitude and phase were derived at wave period of exactly 6.0 days, so that the phases calculated at different heights and latitudes can be compared. In the mesosphere and lower thermosphere (above 50 km), the amplitude structure is symmetric about the equator with peaks at approximately urn:x-wiley:grl:media:grl60100:grl60100-math-005145° latitudes, and the phase tends to be horizontally uniform with downward phase progression. These features are in conformity with the theoretically expected Q6DW in the presence of the mean winds and dissipation (e.g., Salby, 1981b, 1981a). Below 50 km, the phase progression is poleward as well as downward, especially in the SH, indicating equatorward and upward energy propagation from the high‐latitude region. Using reanalysis data, Gan et al. (2018) demonstrated how the Q6DW generated in the SH high latitude can propagate into the NH, growing to be a global mode in the mesosphere and lower thermosphere under September equinox conditions.

In Figure 3c, there is a region of locally enhanced amplitudes at 70–80°S and 20‐ to 50‐km altitude, which can be regarded as a source of the large‐amplitude Q6DW observed above. The amplification of the Q6DW from the seasonal background in this region is depicted in Figure 3e. Enhanced wave activity is observed in the same region over a wide range of wavenumbers ( urn:x-wiley:grl:media:grl60100:grl60100-math-0052 from urn:x-wiley:grl:media:grl60100:grl60100-math-00533 to 3) and periods ( urn:x-wiley:grl:media:grl60100:grl60100-math-0054 = 3–20 days) (not shown here). A possible explanation for the wave amplification is baroclinic/barotropic instability (Gan et al., 2018; Lieberman et al., 2003; Liu et al., 2004; Meyer & Forbes, 1997), in which waves can rapidly grow by extracting energy from the unstable mean flow. Figure 3f shows that the wave amplification in the polar middle atmosphere is not uncommon around this time of year, but in 2019, it took place at lower altitudes ( urn:x-wiley:grl:media:grl60100:grl60100-math-005530 km) than in other years ( urn:x-wiley:grl:media:grl60100:grl60100-math-005650 km).

Figures 3h–3j illustrate the development of the atmospheric instability. The areas highlighted by the light‐yellow color indicate the regions where the necessary condition for barotropic/baroclinic instability is met; that is, the meridional gradient of the quasi‐geostrophic potential vorticity is negative (e.g., Liu et al., 2004). It can be seen that unstable regions are formed mainly around the edge of the polar vortex due to the strong vertical and horizontal shear in the zonal wind. As the westward mean flow descends to lower layers, the unstable regions at high latitudes (70–80°S) also move down and hence exciting/amplifying waves at lower altitudes compared to other years. As these waves propagate equatorward and upward, the amplitude at 45°S is greater than other years above urn:x-wiley:grl:media:grl60100:grl60100-math-005740 km (Figure 3g). As numerically demonstrated by Salby (1981b), the vertical growth of amplitude is enhanced where the zonal mean zonal wind is weak and eastward relative to the phase speed of the wave. The westward phase speed of the Q6DW is urn:x-wiley:grl:media:grl60100:grl60100-math-005855 m/s at 45°S and urn:x-wiley:grl:media:grl60100:grl60100-math-005913 m/s at 80°S. Thus, the reduced eastward mean flow and the weak wind reversal during the SSW (Figures 3h–3j) provide favorable conditions for the vertical propagation of the Q6DW. Interactions of the Q6DW with tides and gravity waves could also affect the vertical structure of the Q6DW (e.g., Forbes et al., 2018; Meyer, 1999). A better understanding of the Q6DW propagation in the mesosphere and lower thermopshere during the September 2019 SSW would benefit from a more comprehensive analysis of dynamic fields from an atmospheric reanalysis or general circulation model.

For the NH, possible SSW influences on the vertical propagation of traveling planetary waves in the middle atmosphere have been discussed in a number of studies (e.g., Hirooka & Hirota, 1985; Gu et al., 2018; Matthias et al., 2012; Pancheva et al., 2008; Sassi et al., 2012; Yamazaki & Matthias, 2019). In some cases, a strong Q6DW was observed during an SSW (e.g., Gong et al., 2018; Pancheva et al., 2018) but in general, there is no one‐to‐one correspondence between the occurrence of SSW and Q6DW enhancement in the NH (Yamazaki & Matthias, 2019). Modeling studies also found enhanced Q6DW activity following some SSWs, which has been attributed to barotropic/baroclinic instability in the NH high latitude (Chandran et al., 2013; Tomikawa et al., 2012). For the SH, studies are few because of infrequent occurrence of SSWs. Dowdy et al. (2004) and Espy et al. (2005) observed a westward propagating planetary wave with ZW1 and period around 14 days at 70‐ to 100‐km altitude during the September 2002 Antarctic SSW. The present study finds a strong response of the Q6DW in the mesosphere and lower thermosphere during the September 2019 Antarctic SSW. It is possible that the response of traveling planetary waves to Antarctic SSWs varies from event to event. More studies are needed to clarify this point.

3 Summary and Conclusions

An SSW occurred in the SH in September 2019. Although it was a minor warming, it involved an exceptionally strong Wave 1 planetary wave and a large polar temperature enhancement by 50.8 K/week. The event also took place under solar minimum conditions, which is preferable for studying the ionospheric response. Earlier studies focused on the effect of NH SSWs on the ionosphere, and few studies investigated SH cases.

The analysis of ionospheric data from European Space Agency's Swarm mission during the September 2019 SSW reveals prominent 6‐day variations in the dayside low‐latitude region, including 20–70% variations in the equatorial zonal electric field, 20–40% variations in the topside electron density, and 5–10% variations in the topside TEC. These variations are attributed to the Q6DW simultaneously observed in the middle atmosphere. Evidence is also found for enhanced lunar tidal perturbations in the ionosphere, but their amplitudes are relatively small (e.g., less than 15% in the topside electron density).

The amplitude of the Q6DW in the lower thermosphere is more than 0.4 km in geopotential height, which is found to be the largest observed by Aura/MLS in the SH since August 2004 and thus cannot be explained merely as a seasonal effect. The latitudinal and vertical structures of the Q6DW suggest that the waves are excited/amplified in the polar region at 30‐ to 40‐km altitude, where the atmosphere is unstable due to strong vertical shear in the zonal wind connected with planetary wave breaking. As the Q6DW grows in the vertical, the wave attains large amplitudes in the lower thermosphere, which drives ionospheric variability.

These results suggest that a SH SSW can lead to ionospheric variability by altering middle atmosphere dynamics and propagation characteristics of large‐scale waves from the middle atmosphere to the upper atmosphere.

Acknowledgments

We thank the NASA Goddard Earth Sciences (GES) Data and Information Services Center (DISC) (https://disc.gsfc.nasa.gov/) for making the Aura/MLS geopotential height data (DOI: 10.5067/Aura/MLS/DATA2008) and MERRA‐2 data (DOI: 10.5067/QBZ6MG944HW0) available. We also thank the European Space Agency (ESA) for providing the Swarm data. The Swarm EEJ and EEF data can be downloaded from the website (ftp://swarm-diss.eo.esa.int/Advanced/EEF/). The Swarm electron density data can be downloaded from the website (ftp://swarm-diss.eo.esa.int/Level1b/Latest_baselines/EFIx_LP/). The Swarm TEC data can be downloaded from the website (ftp://swarm-diss.eo.esa.int/Level2daily/Latest_baselines/TEC/TMS/). The geomagnetic activity index urn:x-wiley:grl:media:grl60100:grl60100-math-0060 was provided by the GFZ German Research Centre for Geosciences (https://www.gfz-potsdam.de/en/kp-index/). The solar activity index urn:x-wiley:grl:media:grl60100:grl60100-math-0061 was downloaded from the SPDF OMNIWeb database (https://omniweb.gsfc.nasa.gov). This work was supported in part by ESA through Contract 4000126709/19/NL/IS “VERA” and by the Deutsche Forschungsgemeinschaft (DFG) Grant YA‐574‐3‐1.