Comparison of Major Sudden Stratospheric Warming Impacts on the Mid-Latitude Mesosphere Based on Local Microwave Radiometer CO Observations in 2018 and 2019

: In this paper, a comparison of the impact of major sudden stratospheric warmings (SSWs) in the Arctic in February 2018 (SSW1) and January 2019 (SSW2) on the mid-latitude mesosphere is given. The mesospheric carbon monoxide (CO) and zonal wind in these two major SSW events were observed at altitudes of 70–85 km using a microwave radiometer (MWR) at Kharkiv, Ukraine (50.0 ◦ N, 36.3 ◦ E). Data from ERA-Interim and MERRA-2 reanalyses and Aura Microwave Limb Sounder measurements were also used. It is shown that: (i) The di ﬀ erences between SSW1 and SSW2, in terms of local variability in zonal wind, temperature, and CO in the stratosphere and mesosphere, were clearly deﬁned by the polar vortex (westerly in cyclonic circulation) and mid-latitude anticyclone (easterly) migrating over the MWR station, therefore; (ii) mesospheric intrusions of CO-rich air into the stratosphere over the Kharkiv region occurred only occasionally, (iii) the larger zonal wave 1–3 amplitudes before SSW1 were followed by weaker polar vortex recovery than that after SSW2, (iv) the strong vortex recovery after SSW2 was supported by earlier event timing (midwinter) favoring vortex cooling due to low solar irradiance and enhanced zonal circulation, and (v) vortex strengthening after SSW2 was accompanied by wave 1–3 ampliﬁcation in March 2019, which was absent after SSW1. Finally, the inﬂuence of the large-scale circulation structures formed in individual major SSW events on the locally recorded characteristics of the atmosphere is discussed.


Microwave Radiometer
The Kharkiv MWR measures the CO profile and zonal wind velocity in the upper stratosphere and mesosphere using an emission line of 115.3 GHz. The radiometer can continuously provide a vertical profile of the mesospheric CO layer, excluding times of adverse weather conditions. The radiometer receiving system is a double-sideband, frequency-switched heterodyne receiver system for measurement of atmospheric CO at 115.3 GHz [43]. The MWR data allows for the retrieval of wind speeds from the Doppler shift of the CO line emission and CO volume mixing ratio (VMR) profiles using radiative transfer calculations with the Qpack package version 1.0.93 [44,45] and the forward model of the Atmospheric Radiative Transfer Simulator (ARTS) [46,47]. This method provides average values of the zonal wind speed at altitudes in the range of 70-85 km with a daily integration time scale. The accuracy of the MWR measurement is mainly determined by the radiometric sensitivity of the receiving system, which depends on the temperature of the radiometer, the temperature of the troposphere, and the signal integration time. The VMR CO error of the microwave measurement does not exceed 30% depending on the weather conditions. More details on the MWR measurements have been presented by Wang et al. [26].
The results of microwave measurements of CO and zonal wind in the mid-latitude mesosphere at 70-85 km altitudes, which is still not adequately covered by ground-based observations, are useful for improving our understanding of the impact of SSWs in this region. In this paper, the daily average zonal wind speed at altitudes of 70-85 km and CO profiles in the mesosphere and upper stratosphere during the January-March 2018 (SSW1) [26] and December 2018-January 2019 (SSW2) intervals are discussed and compared. These data cover the main stages of the major SSW1 and SSW2 events.

Data from Other Databases
Daily data sets from ERA-Interim atmospheric reanalysis [48] were used from Reference [40] to plot local zonal wind velocity altitude profiles and zonal wind field in the NH lower mesosphere and stratosphere. The downloaded ERA-Interim data were compared with the MWR observations. Zonal mean zonal wind and zonal wave amplitudes in geopotential height at 10 hPa were analyzed using the MERRA-2 reanalysis data from the National Aeronautics and Space Administration Goddard Space Flight Center, Atmospheric Chemistry and Dynamics Laboratory (NASA GFC ACDL) site [41].
Aura MLS measurements of the air temperature and CO volume mixing ratio were also analyzed. The MLS CO data accuracy is quoted as +20% to +50% at pressure levels from 1 to 0.0022 hPa [42,[49][50][51].
The description of the data analysis methods used has been given by Wang et al. [26] in their Supplementary Materials.

Results: The Local SSW Effects Over the Mid-Latitude Station
The major SSWs were distinguished from minor events by requiring a reversal (from westerly to easterly) of the zonal winds at 10 hPa, 60 • latitude, and an increase in the zonal mean temperature poleward of 60 • at 10 hPa [5]. In February 2018, easterly zonal winds appeared (up to -20 m s −1 ) appeared and then returned to westerly (at +10 m s −1 ) in March (Figure 1a). On January 2, 2019, there was a wind reversal of up to -10 m s −1 .
Following this, the vortex became strong again in February-April 2019, with the greatest westerly wind being over 50 m s −1 , which was higher than that in the preceding November-December (Figure 1b; see also Butler et al. [8]). Both SSWs had a similar persistence of the mean zonal wind reversal in the stratosphere-19 days in 2018 and 21 days in 2019 [8]-but differed in post-SSW westerly wind persistence (compare blue and red rectangles in Figure 1). The final warming occurred in late April in both events.  [41]. Note the difference in the easterly and westerly persistence and velocity change in the SSW and post-SSW time intervals (outlined by blue and red rectangles, respectively).
Following this, the vortex became strong again in February-April 2019, with the greatest westerly wind being over 50 m s −1 , which was higher than that in the preceding November-December (Figure 1b; see also Butler et al. [8]). Both SSWs had a similar persistence of the mean zonal wind reversal in the stratosphere-19 days in 2018 and 21 days in 2019 [8]-but differed in post-SSW westerly wind persistence (compare blue and red rectangles in Figure 1). The final warming occurred in late April in both events.
The local microwave observations at the Kharkiv site, combined with the MLS and reanalysis data, show wide-ranging daily variability in CO, zonal wind, and temperature in the mesosphere and stratosphere during the SSW1 and SSW2 events. The observed local CO variability can be explained mainly by horizontal and vertical air mass redistribution due to planetary wave activity [26].
Horizontal replacement of the CO-rich polar vortex air by the CO-poor air from more equatorward latitudes (outside the vortex) led to a significant mesospheric CO decrease over the station during SSW1, whereas downward motion caused an enhancement of stratospheric CO at about 30 km [26]. The contributions of large-scale processes to the local atmosphere variability over the station during SSW1 and SSW2 are compared below.

Planetary Wave Activity
A comparison of the planetary wave activity in geopotential height during the SSW1 and SSW2 events is presented in Figure 2. The planetary wave amplitudes were calculated from MERRA-2 daily mean geopotential height data at 60°N, 10 hPa, using Fourier decomposition. In the pre-warming period, the wave 1 amplitudes reached about 1700 m [26] and 1500 m (Figure 2), respectively. At the onset of SSW1, the amplitude of wave 1 decreased, and that of wave 2 increased to about 1000 m in February 2018 [26].  [41]. The SSW2 event is bounded by red vertical  [41]. Note the difference in the easterly and westerly persistence and velocity change in the SSW and post-SSW time intervals (outlined by blue and red rectangles, respectively).
The local microwave observations at the Kharkiv site, combined with the MLS and reanalysis data, show wide-ranging daily variability in CO, zonal wind, and temperature in the mesosphere and stratosphere during the SSW1 and SSW2 events. The observed local CO variability can be explained mainly by horizontal and vertical air mass redistribution due to planetary wave activity [26].
Horizontal replacement of the CO-rich polar vortex air by the CO-poor air from more equatorward latitudes (outside the vortex) led to a significant mesospheric CO decrease over the station during SSW1, whereas downward motion caused an enhancement of stratospheric CO at about 30 km [26]. The contributions of large-scale processes to the local atmosphere variability over the station during SSW1 and SSW2 are compared below.

Planetary Wave Activity
A comparison of the planetary wave activity in geopotential height during the SSW1 and SSW2 events is presented in Figure 2. The planetary wave amplitudes were calculated from MERRA-2 daily mean geopotential height data at 60 • N, 10 hPa, using Fourier decomposition. In the pre-warming period, the wave 1 amplitudes reached about 1700 m [26] and 1500 m (Figure 2), respectively. At the onset of SSW1, the amplitude of wave 1 decreased, and that of wave 2 increased to about 1000 m in February 2018 [26].  [41]. Note the difference in the easterly and westerly persistence and velocity change in the SSW and post-SSW time intervals (outlined by blue and red rectangles, respectively).
Following this, the vortex became strong again in February-April 2019, with the greatest westerly wind being over 50 m s −1 , which was higher than that in the preceding November-December (Figure 1b; see also Butler et al. [8]). Both SSWs had a similar persistence of the mean zonal wind reversal in the stratosphere-19 days in 2018 and 21 days in 2019 [8]-but differed in post-SSW westerly wind persistence (compare blue and red rectangles in Figure 1). The final warming occurred in late April in both events.
The local microwave observations at the Kharkiv site, combined with the MLS and reanalysis data, show wide-ranging daily variability in CO, zonal wind, and temperature in the mesosphere and stratosphere during the SSW1 and SSW2 events. The observed local CO variability can be explained mainly by horizontal and vertical air mass redistribution due to planetary wave activity [26].
Horizontal replacement of the CO-rich polar vortex air by the CO-poor air from more equatorward latitudes (outside the vortex) led to a significant mesospheric CO decrease over the station during SSW1, whereas downward motion caused an enhancement of stratospheric CO at about 30 km [26]. The contributions of large-scale processes to the local atmosphere variability over the station during SSW1 and SSW2 are compared below.

Planetary Wave Activity
A comparison of the planetary wave activity in geopotential height during the SSW1 and SSW2 events is presented in Figure 2. The planetary wave amplitudes were calculated from MERRA-2 daily mean geopotential height data at 60°N, 10 hPa, using Fourier decomposition. In the pre-warming period, the wave 1 amplitudes reached about 1700 m [26] and 1500 m (Figure 2), respectively. At the onset of SSW1, the amplitude of wave 1 decreased, and that of wave 2 increased to about 1000 m in February 2018 [26].  [41]. The SSW2 event is bounded by red vertical the station region (50 • N). Note that the SSW evolution occurred with (without) wave 2 amplification in the first (second) event.
In Figure 3, the latitude-time cross-sections of planetary waves 1-3 amplitudes in geopotential height at 10 hPa in NH for the winters of 2017-2018 (December-April, left) and 2018-2019 (November-March, right) are presented. The differences between the winters of 2017-2018 and 2018-2019 are very noticeable. In winter 2017-2018, all three wave components were more active and persistent than in 2018-2019 during the pre-warming time interval (Figure 3, left and right, respectively). At the same time, wave activity similarly weakened during and immediately after SSWs. This was likely caused by the easterly limitation for upward wave propagation during the SSW and the general seasonal decrease of wave activity later in spring.
lines. The dashed red vertical line marks the beginning of zonal wind reversal from the local microwave radiometer (MWR) measurements in December 2018.
In early January 2019 (SSW2), the wave 1 amplitude also decreased to about 1000 m. However, the wave 2 amplitude was very small at approximately 100 m ( Figure 2). Figures 1 and 2 characterize the conditions of the middle stratosphere (10 hPa) at the polar vortex edge (60°N), which could affect the station region (50°N). Note that the SSW evolution occurred with (without) wave 2 amplification in the first (second) event.
In Figure 3, the latitude-time cross-sections of planetary waves 1-3 amplitudes in geopotential height at 10    The wave 1 and wave 2 peaks were nearly anticorrelated in time, indicating an interaction between them and the zonal flow. In the observed 'wave 1-wave 2 anticorrelation, wave 1 oscillations in response to 'wave-mean flow' interaction played a role in causing the decline of wave 2 to coincide with amplification of wave 1 due to 'wave-wave' interaction [52]. The anticorrelation between wave 1 and wave 2 in the stratosphere during major SSW events has been known for a long time due to observations [53] and modeling [52]. The wave 3 activity changed independently of waves 1 and 2, possibly because wave 3, on an interannual timescale, generally represents a Pacific-North America (PNA) pattern and is less associated with the annular mode [54]. The effects of wave 3, which interacted with waves 1-2, have been shown by Shi et al. [55]. The authors analyzed the modulation effects of wave 3 on wave 1-2 activity fluxes in the stratosphere during SSW in 2005. Wave 3 can enhance (weaken) the poleward wave activity, which produces conditions for upward (downward) propagation of wave 1-2 fluxes into the stratosphere. Wave 1-2 fluxes are suppressed, and wave 3 is able to enter upward into the stratosphere, causing warming to stagnate [55]. The behavior of wave 3 was closely related to the evolution of the Atlantic blocking high and the nearby structures of the temperature and pressure fields. Amplifications of wave 3 in SSW2 in Figure 3 (right) avoided those in wave 1 and wave 2, suggesting an interaction between waves 1-3.
Wave 1-3 activity in winter covered both polar and mid-latitude zones. At the MWR station latitude (50 • N; white dashed line in Figure 3), the role of the higher wave number increased due to the equatorial shift of maximum amplitude noted above. Particularly, the relative contribution of wave 3 to the stratosphere disturbance at 50 • N became more significant during the weakening of waves 1 and 2 (compare Figure 3e,f and Figure 3a-d, respectively, around white dashed line).

Zonal Wind Variability in the Mesosphere and Stratosphere
The MWR observations at 70-85 km over Kharkiv showed that, in SSW2, the mesospheric westerlies started to decrease (between 18 m s −1 on 2 January 2019 and −12 m s −1 on 7 January 2019), with reversal from westerly to easterly around 4 January ( Figure 4c). For comparison, the zonal wind variations in SSW1 are reproduced in Figure 4a, based on Figure 5a of Wang et al. [26]. The westerly in SSW2 recovered after 12 January (instead of 23 January by the SSW criterion for 60 • N, 10 hPa, Figure 1b) and increased to 30 m s −1 on 15 January. There was a similar latitudinal shift of maximum amplitude between waves 1, 2, and 3 in both events: 60-80°N, 50-70°N, and 50-60°N, respectively. There was a similar latitudinal shift of maximum amplitude between waves 1, 2, and 3 in both events: 60-80°N (polar vortex region), 50-70°N (around vortex edge at 60°N, red dashed line in Figure 3), and 50-60°N (sub-vortex region), respectively.
The wave 1 and wave 2 peaks were nearly anticorrelated in time, indicating an interaction between them and the zonal flow. In the observed 'wave 1-wave 2′ anticorrelation, wave 1 oscillations in response to 'wave-mean flow' interaction played a role in causing the decline of wave 2 to coincide with amplification of wave 1 due to 'wave-wave' interaction [52]. The anticorrelation between wave 1 and wave 2 in the stratosphere during major SSW events has been known for a long time due to observations [53] and modeling [52]. The wave 3 activity changed independently of waves 1 and 2, possibly because wave 3, on an interannual timescale, generally represents a Pacific-North America (PNA) pattern and is less associated with the annular mode [54]. The effects of wave 3, which interacted with waves 1-2, have been shown by Shi et al. [55]. The authors analyzed the modulation effects of wave 3 on wave 1-2 activity fluxes in the stratosphere during SSW in 2005. Wave 3 can enhance (weaken) the poleward wave activity, which produces conditions for upward (downward) propagation of wave 1-2 fluxes into the stratosphere. Wave 1-2 fluxes are suppressed, and wave 3 is able to enter upward into the stratosphere, causing warming to stagnate [55]. The behavior of wave 3 was closely related to the evolution of the Atlantic blocking high and the nearby structures of the temperature and pressure fields. Amplifications of wave 3 in SSW2 in Figure 3 (right) avoided those in wave 1 and wave 2, suggesting an interaction between waves 1-3.
Wave 1-3 activity in winter covered both polar and mid-latitude zones. At the MWR station latitude (50°N; white dashed line in Figure 3), the role of the higher wave number increased due to the equatorial shift of maximum amplitude noted above. Particularly, the relative contribution of wave 3 to the stratosphere disturbance at 50°N became more significant during the weakening of

Zonal Wind Variability in the Mesosphere and Stratosphere
The MWR observations at 70-85 km over Kharkiv showed that, in SSW2, the mesospheric westerlies started to decrease (between 18 m s -1 on 2 January 2019 and −12 m s −1 on 7 January 2019), with reversal from westerly to easterly around 4 January (Figure 4c). For comparison, the zonal wind variations in SSW1 are reproduced in Figure 4a, based on Figure 5a of Wang et al. [26]. The westerly in SSW2 recovered after 12 January (instead of 23 January by the SSW criterion for 60°N, 10 hPa, Figure 1b) and increased to 30 m s −1 on 15 January.  mid-stratosphere (at 10 hPa, 60°N, Figure 1b), providing evidence that it was a localized feature in the atmospheric circulation. Polar maps of the zonal wind from the ERA-Interim reanalysis [40] ( Figure 5) may clarify how the hemispheric wind patterns affected the local observations. The white circle in Figure 5 shows the location of the MWR site Kharkiv relative to the westerly maximum outlined by the black contour. Zero wind velocity is shown by the white contour, which separates the easterly (light to dark blue) from the westerly (light green-yellow-red). In the mid-stratosphere, a strong local westerly dominated up to mid-January, after which there was strong deceleration toward zero speed (white contour '0′ in mid-January at dashed horizontal line for 32 km in Figure 4d). The maps in Figure 5 show that Kharkiv was below a spiral (latitudinal-varying) structure of strong zonal wind from December 2018-early January 2019, which turned cyclonical toward the pole. This spiral-jet is a known feature and is due to planetary wave breaking, with a spiral direction due to the breaking of anticyclonic waves [56]. Due to the equatorward shift of the westerly maximum, the Kharkiv region was surrounded by near-zero zonal winds on 8 and 20 January 2019 ( Figure 5c).
As noted above, the two local wind reversals appeared at the stratopause level (horizontal line at 50 km in Figure 4d). Figure 5b confirms that, due to displacements of the westerly maximum, Kharkiv was below the easterly area on 25 December 2018 and 8 January 2019. The location of  (Figure 4d). This reversal event did not appear from the mean zonal wind in the mid-stratosphere (at 10 hPa, 60 • N, Figure 1b), providing evidence that it was a localized feature in the atmospheric circulation. Polar maps of the zonal wind from the ERA-Interim reanalysis [40] (Figure 5) may clarify how the hemispheric wind patterns affected the local observations. The white circle in Figure 5 shows the location of the MWR site Kharkiv relative to the westerly maximum outlined by the black contour. Zero wind velocity is shown by the white contour, which separates the easterly (light to dark blue) from the westerly (light green-yellow-red).
In the mid-stratosphere, a strong local westerly dominated up to mid-January, after which there was strong deceleration toward zero speed (white contour '0 in mid-January at dashed horizontal line for 32 km in Figure 4d). The maps in Figure 5 show that Kharkiv was below a spiral (latitudinal-varying) structure of strong zonal wind from December 2018-early January 2019, which turned cyclonical toward the pole. This spiral-jet is a known feature and is due to planetary wave breaking, with a spiral direction due to the breaking of anticyclonic waves [56]. Due to the equatorward shift of the westerly maximum, the Kharkiv region was surrounded by near-zero zonal winds on 8 and 20 January 2019 ( Figure 5c).
As noted above, the two local wind reversals appeared at the stratopause level (horizontal line at 50 km in Figure 4d). Figure 5b confirms that, due to displacements of the westerly maximum, Kharkiv was below the easterly area on 25 December 2018 and 8 January 2019. The location of Kharkiv below the easterly area in the lower mesosphere (at the 64 km level) occurred on the same dates (Figure 5a; 25 December and 8 January).
Thus, the maps for 25 December 2018 and 8 January 2019 in Figure 5a,b explain how zonal wind reversals over Kharkiv in the middle stratosphere-lower mesosphere in December and January ( Figure 4d) related to the changes in horizontal wind patterns. In general, Figure 5 shows that the polar vortex structure became more fragmented below the stratopause (Figure 5c) than above the stratopause (Figure 5a). It can also be seen that the polar vortex at the stratopause and in the lower mesosphere became more much stronger between the onset and end of the SSW2 event (red to dark red between 2 and 20 January in Figure 5a,b).
The two wind reversals in the lower mesosphere ( Figure 4d) occurred a few days later than in the upper mesosphere (Figure 4c), as was similarly observed during the SSW1 (Figures 4b and 4a, respectively). At the Kharkiv latitudes, the wind reversal in the SSW2 did not penetrate as low into the middle stratosphere (Figure 4d) as that in SSW1 (Figure 4b). Therefore, in terms of time delay with respect to the upper mesosphere, it was impossible to compare the two events at these altitudes.

Temperature Profile Changes
The changes in the vertical temperature profile in the Kharkiv region from the MLS measurements [42]  Kharkiv below the easterly area in the lower mesosphere (at the 64 km level) occurred on the same dates (Figure 5a; 25 December and 8 January). Thus, the maps for 25 December 2018 and 8 January 2019 in Figure 5a,b explain how zonal wind reversals over Kharkiv in the middle stratosphere-lower mesosphere in December and January (Figure 4d) related to the changes in horizontal wind patterns. In general, Figure 5 shows that the polar vortex structure became more fragmented below the stratopause (Figure 5c) than above the stratopause (Figure 5a). It can also be seen that the polar vortex at the stratopause and in the lower mesosphere became more much stronger between the onset and end of the SSW2 event (red to dark red between 2 and 20 January in Figure 5a,b).
The two wind reversals in the lower mesosphere ( Figure 4d) occurred a few days later than in the upper mesosphere (Figure 4c), as was similarly observed during the SSW1 (Figures 4b and 4a, respectively). At the Kharkiv latitudes, the wind reversal in the SSW2 did not penetrate as low into the middle stratosphere (Figure 4d) as that in SSW1 (Figure 4b). Therefore, in terms of time delay with respect to the upper mesosphere, it was impossible to compare the two events at these altitudes.

Temperature Profile Changes
The changes in the vertical temperature profile in the Kharkiv region from the MLS measurements [42]   It appears that the change between the elevated and descended stratopause over the station region (Figure 6a) was influenced by the local transition from westerly to easterly due to replacement of the polar cyclone by mid-latitude anticyclone (Figure 5a,b), which, in turn, was associated with the action of zonal wave 1 (Figure 3b). Climatology has demonstrated similar anomalies in the geographical structure of polar winter stratopause temperature and height, with It appears that the change between the elevated and descended stratopause over the station region (Figure 6a) was influenced by the local transition from westerly to easterly due to replacement of the polar cyclone by mid-latitude anticyclone (Figure 5a,b), which, in turn, was associated with the action of zonal wave 1 (Figure 3b). Climatology has demonstrated similar anomalies in the geographical structure of polar winter stratopause temperature and height, with respect to the location of the polar vortices and anticyclones in both observations [15,16,20] and numerical modeling [17,22].
The temperature altitude-time distribution shows that the overall regional temperature in SSW2 was 10-15 • C higher (Figure 6a) than in the SSW1, as shown by Figure 6 in Wang et al. [26]. The regional temperature maximum in the stratopause region was relatively stable around the 50 km height in SSW1 but was broader (see [26]) compared to that of SSW2 ( Figure 6) in both regional and zonal mean data. The downward propagation of the increased temperature anomaly between the upper and lower stratosphere observed in the first event [24] was not observed in the second event ( Figure 6).
The regional stratospheric temperature in the descended stratopause region increased sharply by 15-20 • C during the SSW2 event, accompanied by a sharp decrease in the mesospheric temperature (Figure 6a), in agreement with the significant anticorrelation between mesospheric and stratospheric temperatures during SSWs (see, e.g., Siskind et al. [57]). Similar to SSW1 [26], the zonal mean temperature changes for 47.5-52.5 • N in the SSW2 were not as significant and sharp (Figure 6b) as the regional ones (Figure 6a). Figure 6b confirms that the variability in the zonal means in the mid-latitude atmosphere was significantly smoothed out compared to the local and regional variability (Figures 4 and 6a). The latter strongly depended on the zonal asymmetry between the polar vortex and mid-latitude anticyclones, their displacements, and possible vortex fragmentation, displaying relationships between the zonal wave 1-3 amplitudes.

CO Variability
In this section, we analyze changes in the stratosphere and mesosphere over the Kharkiv site According to the MWR data, the local CO maximum was between 90 and 100 km, with peak values of about 18 ppmv in mid-December 2018 and at the beginning of SSW2, 2-5 January 2019 (Figure 7a). The CO concentration decreased by a third between the pre-warming (~15 ppmv) and post-warming (~10 ppmv) periods. This change at the maximum of the CO layer differed from that in SSW1, in which more stable CO levels (around~15 ppmv) were observed in winter-spring (January-March) 2018, although the CO peak at 18 ppmv was the same (Wang et al. [26], Figure 3a).
A difference also existed at the lower edge of the mesospheric CO layer, as shown by the 6 ppmv level (white curve in Figure 7b). The sharp elevation of the 6 ppmv contour (by about 8 km) observed in SSW1 [26] was completely absent in SSW2 (Figure 7b).
At the same time, a general slow increase in the height of the 6 ppmv boundary, starting at 75 km, was observed in both events. This increase was due to the seasonal trend in CO concentration between winter and summer at this altitude, as shown by Solomon et al. [27] and Lee et al. [58].
As can be seen from the CO maps near the height of the CO layer maximum (at 86 km), Kharkiv was located inside the CO-rich polar area during SSW2 (Figure 8a-c) and outside of it during the post-warming period (Figure 8d). This explains the sharp post-warming CO decrease in Figure 7a. In contrast, the Kharkiv location was in the CO-poor area at the CO layer maximum in SSW1, which resulted in the CO decrease just during the SSW1 [26].
Thus, the change in the position of the edge of the mesospheric polar vortex-which bounds the high-CO area-relative to the ground-based instrument may play a more important role in modulation of the observed CO level than the general evolution of the SSW itself. The 0.5 ppmv level in SSW2 varied over the same altitude range of 50-65 km, but was several kilometres higher in mean altitude (Figure 7c) than that in SSW1 [26]. At the stratopause level around 50 km, a lower (higher) CO amount was observed in the first (second) half of January 2019 (Figure 7d), again determined by the Kharkiv location being outside (inside) the CO-rich area (Figure 8i-l).
Regionally averaged MLS CO data showed a post-warming CO decrease by about 10 ppmv, which is generally consistent with the local MWR data ( Figure 7a) described above. The 6 ppmv and 0.5 ppmv levels also demonstrated very close variability in the local MWR data (white curves in Figure 7b,c) and the regional MLS data (white curves in Figure 7f,g). The changes in stratospheric CO during the beginning of SSW2 were unexpectedly different in the two data types (Figure 7d,h, respectively). The satellite data for the Kharkiv region in the first half of January showed the presence of mesospheric CO levels up to 0.4-0.6 ppmv in the middle-upper stratosphere (30-50 km, Figure 7h). A similar CO increase to 0.1-0.2 ppmv in the stratosphere was seen in the zonal mean MLS data (Figure 7l).
The difference between the MWR and MLS data can be explained by a small fragment of CO-rich air near Kharkiv (Figure 8n). This fragment did not appear in the local MWR record (Figure 7d), but it was partially covered by the regional MLS data (Figure 7h) averaged over a 5-degree latitudinal zone (47.5-52.5 • N) and a 20-degree longitudinal sector (26-46 • E), as well as by zonal mean data (Figure 7l). Generally, local and regional CO variability from MWR and MLS observations (Figures 4c, 7a-h and 8) is consistent with change in the polar vortex structure seen from the ERA-Interim reanalysis ( Figure 5). Thus, the change in the position of the edge of the mesospheric polar vortex-which bounds the high-CO area-relative to the ground-based instrument may play a more important role in modulation of the observed CO level than the general evolution of the SSW itself. The white circle shows the location of the MWR Kharkiv site relative to the high and low CO content areas, which are separated by black contours. Note that Kharkiv fell under the high CO area in the mesosphere during SSW2 (b,c,f,g), which was displaced relative to Kharkiv afterward (d,h). This displacement explains the stepwise decrease in the post-warming CO level from the MWR data in Figure 7a and from regional MLS data in Figure 7e. The 0.5 ppmv level in SSW2 varied over the same altitude range of 50-65 km, but was several kilometres higher in mean altitude (Figure 7c) than that in SSW1 [26]. At the stratopause level around 50 km, a lower (higher) CO amount was observed in the first (second) half of January 2019 (Figure 7d), again determined by the Kharkiv location being outside (inside) the CO-rich area (Figure 8i-l).
Regionally averaged MLS CO data showed a post-warming CO decrease by about 10 ppmv, which is generally consistent with the local MWR data ( Figure 7a) described above. The 6 ppmv and 0.5 ppmv levels also demonstrated very close variability in the local MWR data (white curves in Figure 7b,c) and the regional MLS data (white curves in Figure 7f,g). The changes in stratospheric CO during the beginning of SSW2 were unexpectedly different in the two data types (Figures 7d,h,  respectively). The satellite data for the Kharkiv region in the first half of January showed the presence of mesospheric CO levels up to 0.4-0.6 ppmv in the middle-upper stratosphere (30-50 km, Figure 7h). A similar CO increase to 0.1-0.2 ppmv in the stratosphere was seen in the zonal mean MLS data (Figure 7l).
The difference between the MWR and MLS data can be explained by a small fragment of CO-rich air near Kharkiv (Figure 8n). This fragment did not appear in the local MWR record ( Figure  7d), but it was partially covered by the regional MLS data (Figure 7h) averaged over a 5-degree latitudinal zone (47.5-52.5°N) and a 20-degree longitudinal sector (26-46°E), as well as by zonal mean data (Figure 7l). Generally, local and regional CO variability from MWR and MLS observations (Figures 4c, 7a-h, and 8) is consistent with change in the polar vortex structure seen from the ERA-Interim reanalysis ( Figure 5). Zonal mean MLS data showed an even sharper mesospheric CO maximum (near 97 km, Figure 7i) than local MWR and regional MLS data (Figure 7a,e, respectively). In addition, zonal means revealed the downward motion of increased CO amount in the lower mesosphere and stratosphere (red dashed lines in Figure 7k,l) with velocities of about 220-250 m day −1 (~280 m day −1 from the regional MLS data in Figure 7h). These estimates are in agreement with the results for the pre-warming period of SSW1 ( [26], Figure 3), in which descent tendency in the post-warming period has not been identified. The red dashed lines in Figure 7k,l show that CO-rich mesospheric air took about 2 months to penetrate from the lower mesosphere (~60 km) into the upper stratosphere (~40 km) by crossing the stratopause (~50 km). From ERA-Interim data for the Kharkiv region, the somewhat higher velocity of 500 m day −1 was observed in the post-warming descent of the zonal wind maximum in the same altitude range (white dashed line in Figure 4d).

Discussion
In Section 3, it was shown that local zonal wind, air temperature, and the CO amount observed above the ground-based station during a SSW are dependent on the evolution of the large-scale circulation influenced by planetary wave activity. The SSW and post-SSW effects in the stratosphere and mesosphere are largely determined by the pre-warming vortex state and wave activity [8,9,39,59].

Zonal Waves, Zonal Wind, and Temperature
The stronger pre-warming wave activity in SSW1 (Figure 3) caused more complete destruction of the vortex, without significant following recovery (Figure 1a), in contrast to SSW2 (Figure 1b). In addition, SSW1 occurred in February, when the polar vortex tends to weaken seasonally (black curve in Figure 1). The timing of SSW2 occurred in January-a month and a half earlier-when minimal solar radiation was reaching the Arctic, allowing for enhanced radiative cooling in the polar region and a post-warming transition to a strong vortex [39]. The westerly zonal wind at 10 hPa, 60 • N in March 2019 reached 40-50 m s −1 (Figure 1b; Lee and Butler [39]). These conditions were favorable for tropospheric planetary wave propagation into the stratosphere, as confirmed by the increase of wave 1-3 amplitudes in March 2019 (Figure 3b,d,f). Similar wave amplification did not appear in the post-SSW1 period in March 2018 (Figure 3a,c,e) because of a much weaker westerly zonal wind of 5-10 m s −1 (Figure 1a).
The advantages of our results are the detailed identification and comparison of SSW1 and SSW2 in terms of the locally observed changes in relation to variability in the large-scale circulation structures compared in other works [8,9,39]. As shown in Sections 3.2 and 3.4, local mid-latitude measurements of the mesospheric wind and CO were strongly influenced by the location of the polar vortex and mid-latitude anticyclone, as well as their fragments, relative to the MWR station. A similar dependence was found in the local and regional stratopause temperature and height (Section 3.3).
Strong anticyclonic circulation adverts polar vortex air toward the equator, causing zonal asymmetries in temperature and chemical conditions in polar and mid-latitude zones due to zonal wave 1 dominance [60]. The polar maps in Figure 5 show the variable roles of cyclonic (westerly, red) and anticyclonic (easterly, blue) zonal circulation over the station (white circle) in the locally observed variations of the zonal wind velocity (Figure 4c). Note that small, well-defined fragments of vortex and anticyclone remnants, which usually exist after an SSW [16], appeared in the onset of SSW2 (2 January 2019 in Figure 5), when wind reversal started (Figure 4c). They disappeared with the polar vortex recovery at the stratopause and in the lower mesosphere (20 January in Figure 5a,b).
The stratopause temperature (height) over the Kharkiv region in the pre-warming period was lower (higher) than that during SSW2 (Figure 6a), due to the transition from westerly (in cyclonic circulation) to easterly (in anticyclonic circulation ( Figure 5)). This is consistent with climatology, which has demonstrated similar differences in winter stratopause temperatures and heights with respect to the location of the polar vortices and anticyclones [22], as well as with case study results [16], composite event analysis [15], and model data [23,31].
In general, the differences between SSW1 and SSW2, in terms of local and regional wind and temperature over the MWR station at Kharkiv, were largely determined by causal relationships between the locations of the polar vortex, the mid-latitude anticyclone, and the fragments formed under planetary wave influence during the individual event.

Descent of the Mid-Latitude CO Anomalies
Due to the sharp meridional CO gradient at the polar vortex edge in winter, zonal asymmetry of the vortex accompanied by latitudinal displacements due to wave effects should dramatically affect local CO densities [27]. When comparing close dates, one can see that the CO-rich areas in Figure 8 are located poleward of the vortex edge outlined by black contours for zonal wind maximum in Figure 5 (i.e., inside the polar vortex). As in SSW1 [26], the vortex asymmetry, deformation, and fragmentation in the mesosphere led to appearance of high polar CO levels in the mid-latitudes of 30-50 • N (see contours of 16 ppmv in Figure 8a-d and of 6 ppmv in Figure 8e-h).
The CO descent occurred mainly through the core of the highly displaced vortex, accompanied by an intrusion of CO-rich air from the mesosphere into the mid-stratosphere, which appears to be a defining signature of winters with major SSWs [31]. The mesospheric intrusions resulted in the increase of CO levels to 0.1-0.4 ppmv in the mid-latitude stratosphere around 50 • N (30-50 km; Figure 7h,l for regional and zonal mean data, respectively). It is worth noting that the typical CO level in the mid-latitude stratosphere in winter is about 0.01-0.02 ppmv [27,29,30,61]. It should also be noted that the typical lifetime of CO in the stratosphere is about 10-20 days [28] and may be about 1 week in the presence of sunlight [27]. In the mesosphere, its lifetime is about 2 months, and in regions of total darkness, it is extremely long, and CO is approximately conserved [28]. As the SSWs in 2018 and 2019 evolved near the minimum in the 11-year solar activity change between cycles 24 and 25 (see, e.g., in [62]), the influence of solar irradiance on the CO lifetime [58] did not differ in the two events.
Due to its relatively low lifetime in the stratosphere, the increased CO amount in December and early January (Figure 7h,l) was continuously maintained by downward CO transport from the mesosphere. As noted in Section 1, the descent of CO from the mesosphere to the stratosphere induced by planetary wave activity may occur as low as 25-30 km [16,[26][27][28][29][30][31].
The CO descent velocities in SSW2 (220-280 m day −1 (Figure 7h,k,l)) are consistent with those observed in SSW1 with MWR at Kharkiv by Wang et al. [26] and in previous SSW events from ground-based and satellite measurements [30,59,63,64]. The CO-rich air was concentrated more uniformly inside the polar vortex region in the mesosphere (Figure 8a-h), was very heterogeneously distributed in relatively small areas at the stratopause level (Figure 8i-l) and was present only in small fragments in the stratosphere (Figure 8m-p). Therefore, registration of the CO descent in the stratosphere by local, regional, or zonal mean data depends on whether these fragments fall into the field of view of the ground-based instrument or into a certain longitudinal sector and latitudinal zone when satellite data are used. The results of Figure 8m-p and those presented in Figure 4m-p by Wang et al. [26] provide evidence that intrusions of the CO-rich air from the mesosphere into the stratosphere occur in spatially limited areas in high-and mid-latitude zones.
The 6 ppmv level in the mesosphere, on the contrary, showed a steady rise in altitude (between about 70 km and 85 km (Figure 7b,f,j, Wang et al. [26], Figure 3)). As the CO-rich air at these altitudes was replaced by CO-poor air (seen also in Figure 8a-h), a gradual decrease in the CO level took place during December-February (January-March in Wang et al. [26]). This tendency was consistent with the seasonal decrease in mesospheric CO level between winter and summer [27,58].

Conclusions
MWR measurements over Kharkiv, Ukraine (50.0 • N, 36.3 • E), located geographically in the mid-latitude region of the Eastern Europe, were presented in this paper. The local and regional manifestations of the stratospheric warming event in January 2019 (SSW2) were compared with those occurring in February 2018 (SSW1; [26]). The comparison was based on both original data from local MWR records and regional and zonal mean data from reanalyses and Aura MLS measurements. These results have advantages over earlier studies, presenting a detailed identification of the origins of the locally observed anomalies. This was achieved through analysis of the large-scale structure migration over the station, which occurred due to the interaction of the polar vortex with the mid-latitude anticyclone and could be seen from the zonal wind pattern in the ERA-Interim reanalysis data. The characteristic changes and distinctions in the mesosphere, stratosphere, and at the stratopause level between the pre-and post-warming periods (including the SSW events) were established.
The local variability of CO was explained by zonal and meridional migrations of the CO-rich polar vortex air mass, accompanied by changes in vortex asymmetry, deformation, and fragmentation. Fragments of the vortex or sub-vortex area, which appeared over Kharkiv during the SSW evolution, determined the CO levels recorded by MWR and MLS. Different locations of the vortex edge relative to the Kharkiv location resulted in the main difference between the major SSW1 and SSW2 events in terms of CO change in the pre-warming, warming, and post-warming periods. During the main phase of the SSW2 event, the CO layer maximum within the polar vortex was located over Kharkiv (Figure 8a-c) and moved away in the post-warming period. This explained the sharp post-warming CO decrease observed for this event (Figure 8d), unlike the situation for SSW1, during which Kharkiv was generally located outside the vortex edge at the time of the warming. A general slow increase in the height of the lower edge of the CO layer (at the 6 ppmv level), starting at 75 km, was observed during both events. This increase was consistent with the seasonal trend in CO concentration between winter and summer at this altitude, as shown earlier by Solomon et al. [27] and Lee et al. [58].
Fewer and fewer fragments of CO-rich air appear with decreasing height between the mesosphere and stratosphere. Therefore, changing the position of the polar vortex, which is unevenly filled with CO, affects whether the descent of CO in such fragments into the stratosphere will be observed by a ground-based instrument. The results of this work and those of Wang et al. [26] provide evidence that mesospheric intrusions of CO-rich air into the stratosphere in the major SSW1 and SSW2 events occurred only in spatially limited areas in the high-and mid-latitudes.
The mesospheric CO evolution in the major SSW1 and SSW2 events was broadly consistent between the MWR and MLS data. This was due to the more uniform filling of the CO vortex and the lesser effect of vortex fragmentation on local observation compared to the stratospheric CO.
In addition, we note that, due to the reduction in the 11-year solar activity between cycles 24 and 25, the change in solar irradiance did not cause a significant difference in the CO lifetime between the major SSW events in 2018 and 2019.
It was shown that: (i) The larger zonal wave 1-3 amplitudes before SSW1 were followed by a weaker polar vortex recovery than that after SSW2; (ii) the strong vortex recovery after SSW2 was supported by earlier event timing (midwinter); [39]; (iii) vortex strengthening after SSW2 was accompanied by wave 1-3 amplification in March 2019 (Figure 3b,d,f), which did not occur after SSW1 (Figure 3a,c,e); and (iv) the differences between SSW1 and SSW2, in terms of variability in zonal wind, temperature, and CO in the stratosphere and mesosphere, were clearly defined by the polar vortex (westerly in cyclonic circulation) and mid-latitude anticyclone (easterly) which migrated over the station (Figures 5 and 8).
These results demonstrate that the large-scale circulation structures formed during the individual SSW events under the influence of planetary waves may play a major role in the locally observed variability in zonal wind, air temperature, and trace gases.