Exceptionally persistent Eurasian cold events and their stratospheric link

Persistent boreal winter cold spells (PCEs) can heavily strain the economy and significantly impact everyday life. While sudden stratospheric warmings are considered a precursor for Eurasian (EUR) cold events, these temperature extremes may occur during the full range of stratospheric variability. We investigate PCEs relative to the prevailing stratospheric polar vortex regime before their onset, with a particular focus on extremely weak (SSW) and strong (SPV) stratospheric winds by performing (lagged) composite analysis based on ERA5 reanalysis. On average, SPV PCEs that are concentrated over central EUR, are colder, shorter and set in more abruptly compared to SSW PCEs. A quasi-stationary, mid-tropospheric anticyclone over the Arctic Ocean that blocks warm air advection toward EUR is connected to the canonical downward progression of the negative North Atlantic Oscillation for SSW PCEs. In contrast, during SPV PCEs, the anticyclone is part of a Rossby wave having an origin co-located with negative wave activity flux anomalies over and being influenced by stratospheric wave reflection toward the North Atlantic. Its slow east-ward propagation is likely related to Arctic surface warming and unusually weak zonal winds over EUR.

From a remote perspective, Eurasian weather has been shown to be strongly influenced by the circulation over the North Atlantic (e.g. Barnett et al. 1984) whose extendedrange predictability can be inferred from the dynamic state of the stratosphere (e.g. Domeisen et al. 2020). In winter, increased radiative cooling over the pole facilitates strong temperature contrasts between the mid-latitude and polar stratosphere, resulting in strong, quasi-zonal, circumpolar westerlies (Waugh and Polvani 2010;Waugh et al. 2017). Vertically propagating planetary waves from the troposphere may deposit easterly momentum onto the stratospheric polar vortex (McIntyre and Palmer 1983) via wave breaking (Charney and Drazin 1961), thereby disturbing the geostrophically balanced flow and fostering an abrupt increase of polar cap temperatures (e.g. Matsuno 1971;O'Neill et al. 2015). Positive polar cap temperature as well as geopotential height anomalies, corresponding to the negative phase of the Northern Annular Mode, may migrate downward into the troposphere (e.g. Dunkerton 1999, 2001; Thompson et al. 2002;Hitchcock and Simpson 2014;Kodera et al. 2016) through wavemean-flow interaction, manifesting at the surface as the negative phase of the Arctic Oscillation and the North Atlantic Oscillation (NAO) (e.g. Baldwin and Dunkerton 2001;Hitchcock and Simpson 2014). Cold Northern European and Eurasian surface temperature anomalies may develop in response (Thompson and Wallace 2001;Kretschmer et al. 2018aKretschmer et al. , 2018bKing et al. 2019;Finke and Hannachi 2022), that have an increased predictability through this stratosphere-troposphere coupling (Statnaia et al. 2022). Further, a reduction in meridional temperature, and thus, geopotential height gradient due to an unusually warm Arctic may decrease large-amplitude Rossby wave progression over the North Atlantic, potentially increasing the persistence of e.g. downstream cold spells (Francis and Vavrus 2012). Chen et al. (2021) found Ural blocking to be more persistent when Barents-Kara sea warming is weak, while a strong warming yields a weakened stratospheric polar vortex with a negative NAO response and subsequent downstream progression of Ural blocking.
A non-zonally symmetric configuration (Kodera et al. 2008(Kodera et al. , 2013(Kodera et al. , 2016Kretschmer et al. 2018b) or an unusually strong polar vortex in the high-latitude stratosphere (e.g. Perlwitz and Graf 2001;Harnik 2003, 2004), resulting from undisturbed radiative cooling, instead, may constrain the upward wave propagation by allocating a layer that reflects planetary waves back down to the troposphere where they may influence the structure of planetary waves (Perlwitz and Harnik 2004). The consequent channeling of planetary waves in the troposphere fosters the development of the positive phase of the Northern Annular Mode and NAO (e.g. Baldwin and Dunkerton 2001;Scaife et al. 2005) corresponding to warm surface temperature anomalies in northern Europe and eastern Asia (Thompson et al. 2002).
Stratospheric polar vortex variability in the sense of changes in mid-stratospheric zonal mean zonal wind, can thus yield surface temperature changes in Eurasia, i.e. cold (warm) temperature anomalies have been linked to a weak (strong) stratospheric polar vortex. However, Eurasian cold spells may also occur during episodes with unusually strong stratospheric vortex conditions. We explore common and distinctive statistical and dynamical characteristics of persistent Eurasian cold spells after strong and weak polar vortex events. The top-to-bottom approach has been used extensively to study the stratospheretroposphere connection. Complementing, our study is based on a bottom-up approach, i.e. analysis based on the time evolution of the tropospheric event, thus enlarging the comprehensive understanding of mechanisms inducing Eurasian cold spells, ultimately aiming for an improvement of their extended-range prediction.

Data
In the present study, daily averages of 6-hourly ERA5 reanalysis data (Hersbach et al. 2020) are considered. The data consist of sea level pressure (SLP), geopotential, 3D wind and 2m air temperature (T2m), for the 40-year time period from 1979 to 2019 for the November to March (NDJFM) season on a 1 • × 1 • regular latitude-longitude grid and 37 pressure levels from 1000hPa to 1hPa. The 29th February is excluded from the analysis.

Study Region
The objective of the study is to identify and classify cold extremes and investigate their potential linkage to stratospheric variability. An empirical orthogonal function (EOF) analysis has been applied to the 5-day running mean T2m anomalies with respect to their climatology, i.e. 40year daily mean, spatially bounded between the 20 • N and 90 • N latitude, to identify a geographical region of interest in Eurasia. Composites of T2m anomalies for the first 4 modes of variability, explaining 11.65%, 11.05%, 7.67% and 6.37%, respectively, for the corresponding standardized principle components (PC) with values higher than 1, are shown in Fig. 1. The black box, located over a vast region of Eurasia (EUR) [50 • N -65 • N; 15 • E -130 • E], identifies the study region that was chosen based on the main regions of variability corresponding to EOF1, EOF3 and EOF4. Accordingly, T2m anomaly composites based on days with PC1 > 1 (a), PC3 > 1 (c) and PC4 > 1 (d) reflect cold temperatures over parts of Eurasia. The region over the eastern US, emerging in PC2 > 1 (b), PC3 > 1 and PC4 > 1, is left for future research.

Construction of Indices and Event Classification
The daily, 5-day running mean T2m is latitude-weighted by multiplication of each grid point value by the squareroot of the cosine of latitude before averaging over EUR to construct an index that is shown in Fig. 1e. We find a distribution with significant negative skewness (-0.23), indicating a longer tail on the cold temperature side and reflecting the intermittent occurrence of extreme cold events. A persistent cold event (PCE) is defined if the T2m index is below one standard deviation (−2.53K; dark blue values in Fig. 1e) for at least 6 consecutive days of which the first day of threshold exceedence is taken as the first event day. If the last day of one event and the first day of the following event are less than 5 days apart and, additionally, the T2m anomaly does not increase to -0.25 standard deviations, the events are counted as one, yielding 58 events in total. No major changes of the overall results occur when varying these two conditions slightly.
The strength of the stratospheric polar vortex is evaluated using the 10hPa zonal mean zonal wind at 60 • N as an index shown in Fig. 1f, whose distribution is slightly flat (kurtosis: 2.78) and negatively skewed (skewness: -0.18). The onset day of an extreme stratospheric vortex is taken when the zonal mean wind turns easterly, commonly associated with a sudden stratospheric warming (SSW), or exceeds 48m/s, i.e. a strong stratospheric polar vortex (SPV).
The cold events at the surface are divided into different groups depending on the occurrence of a stratospheric vortex event in the 40 days prior to their onset. Thus, we classify 39 neutral, 9 SSW PCEs, and 11 SPV PCEs that are listed in Table 1 including the date of the stratospheric event (first column), time to PCE (second column), date of PCE (third column) and length of PCE (fourth column). Note, that one of the events (06.03.1980) is counted as both SPV PCE and SSW PCE and two of the SSW PCEs share the same detected stratospheric signal prior to their onset. The type of SSW events (vortex split (S) or vortex displacement (D)) evaluated by Hall et al. (2021) have been matched with the detected SSW onset days in this study (Table 1b; fifth column). Among the detected SSW PCEs are 5 vortex split and 3 vortex displacement events. Mitchell et al. (2013) and Kretschmer et al. (2018b), for example, find SSW associated with the latter type to induce unusually cold temperatures in North America while an exceptional cooling over Eurasia is rather associated with the former event type. However, Maycock and Hitchcock (2015) shows that the difference across responses is weak. Therefore, and additionally keeping the sample size of our study in mind, we have not further distinguished these two groups in the analysis. The gap (second column) between SPV/SSW date (first column) and PCE onset date (third column) as well as the length of the latter (fourth column) are given in days. The type of SSW is given in the fifth column in (b)

Statistical and Dynamical Analysis
The quantile regression (QR) model (Koenker and Hallock 2001) has been used to analyze trends in the distribution of extreme cold T2m anomalies corresponding to the lowermost quantiles 0.01 and 0.05 as well as the median for comparison. As described in Elsner et al. (2008), gridded fields of T2m anomalies are used as the response variable and time as the predictor variable. Composites of various variables for SSW PCE and SPV PCE onset days as well as for negative and positive lags are analysed to expose common and distinctive characteristics and patterns. Significance at the 10% level is evaluated based on 10000 bootstrap samples (Efron and Tibshirani 1993) unless stated differently.
The Plumb-flux (Plumb 1985), a diagnostic tool allowing for the derivation of the three dimensional flux of waveactivity (WAF), is almost parallel to the group velocity and thus energy propagation of the flow in the slowly varying limit. The three dimensional vector is given by Eq. 1.
Here, p (p/1000hP a) is the pressure, φ and λ are latitude and longitude, a and Ω are the radius and rotation rate of the earth, N the Brunt-Vaisälä frequency, Ψ the perturbation streamfunction and z = −H ln p the height with constant scale height H . Positive values of the vertical component correspond to upward propagating waves.
The vertical configuration of zonal stratospheric winds may provide for a reflective surface that redirects upward travelling planetary waves back down to the troposphere. The reflective state of the stratosphere can be evaluated by following the approach in Nath et al. (2014), who use the meridionally averaged (58 • N-74 • N) zonal wind to calculate the time evolution of the reflective index across longitudes with Eq. 2.
The stratosphere is reflective in case the zonal wind in 2hPa (U 2) is less than in 10hPa (U 10), i.e. negative U r .

Persistent Eurasian Cold Spells
Using QR to evaluate trends in gridded T2m anomaly fields with a focus on extremely cold EUR temperatures, we find that for the 0.01 (Fig. 2a) and 0.05 ( Fig. 2b) quantiles, the most eastern and most western parts of the study region show a positive trend that is substantially higher than what the median QR ( Fig. 2c)   left-hand y-axis) including 90% confidence interval and T2m absolute value index composite (grey; right-hand y-axis). f) Histogram of PCE's durations exceeding -8K and, interestingly, partly overlap with the region covered by the negative QR coefficients for the 0.01 and 0.05 quantile. The composite time evolution of the EUR T2m anomaly (absolute value) index for 60 days before to 60 days after the onset is shown as the blue (grey) line in Fig. 2e, including the 90% confidence interval as shading, to analyse the extent and change of intensity of PCEs with time. The dashed line indicates the threshold of one negative standard deviation corresponding to the T2m anomaly index. On average the T2m drops to negative anomalies about 10 days before, reaching peak absolute values of about -13 • C with anomalies remaining negative until 50 days after the PCE onset. Note that the index is not significantly different from zero when the confidence interval overlaps the zero line, which is the case between -3-days to +10-days lag. This can partially be reasoned by varying duration among the PCEs, shown in a Histogram in Fig. 2f. We find that 28 PCEs last between 6 and 11 days, another 17 between 12 and 17 days and 7 events between 18 and 23 days. Further, three PCEs are longer than three weeks, i.e. between 24 and 29 and another three even longer than one month with one being exceptionally long-lasting starting on the 8th December 2009 until the 12th February 2010, i.e. 67 days.

Characteristics Depending on the Stratosphere
Following the bottom-up perspective, we group the PCEs depending on whether they are preceded by an SPV (11 events), SSW (9 events) or a neutral state in the stratosphere (39 events) to evaluate differences across groups. Figure 3a shows the distribution of PCE's by month. In total, most PCEs develop in December (14) and January (13), followed by November (12), March (10) and lastly February (9) with the least amount of detected PCEs. More SPV PCEs occur during December and January compared to February and March. In contrast, during February and March SSW PCEs are more common than during December and January. Interestingly, in February, PCEs are more common to emerge after an SSW event than a neutral state of the PCEs can thus occur during the full range of stratospheric variability. In the following we show common and distinct characteristics of each group. The duration of neutral PCEs from the first day of detection (x-axis) is given in Fig. 3b. Of 39 events in total, 18 are shorter than or equal to 10 days, 18 between 10 and 20 days and 3 events longer than 20 days. Consistent with the overall positive trend in outer quantiles over Eurasia (see e.g. Fig. 2b), fewer events are detected during the second half of the time series, i.e. between 1999-2019. Similarly, for SSW PCEs (SPV PCEs) we find 6 (7) events in the first and 3 (3) in the second half of the time series. In the following, we focus on SPV PCEs and SSW PCEs. An analysis of neutral PCEs can be found in the supplementary information. Figure 3c shows the time between SPV and the onset of the PCE (transparent) and the duration of the corresponding PCE (solid color). Here, the date on the x-axis marks the day of the initial stratospheric SPV signal. Three PCEs start within 10 days after a SPV, which then last around 10 days. Six PCEs start between 20 and 30 days after a SPV with durations ranging between 6 and 33 days. Lastly, two PCEs begin after 37 and 40 days after a SPV, lasting for 12 and 16 days, respectively. Figure 3d shows the same as Fig. 3c but for PCEs following SSWs. Within 11 days after an SSW event, 5 PCEs have been detected that have durations between 11 and 26 days. The remaining 4 PCEs developed 4 to around 5.5 weeks after an SSW of which two last less than 2 weeks and two longer than three weeks.
Noticeably, 8 out of 11 SPV PCE events start more than 20 days after the stratospheric signal while in case of SSW PCEs 5 out of 9 start no later than 11 days after the stratospheric signal.

Evolution
In the next step, we investigate differences in the evolution of extent and intensity of SSW PCEs and SPV PCEs with time. Figure 4a-b show a lagged composite of the EUR T2m anomaly (absolute value) index from 60 days before to 60 days after the onset day as the coloured (grey) line including its 90% confidence interval as shading. In case of SPV PCEs (Fig. 4a), the T2m index is mainly above average until about 5 days before the PCE onset. Thereafter, the T2m index is below one standard deviation for about 10 days and below the mean for almost one month. At higher lags, the T2m index is mostly above the mean. Further, the composite temperature anomalies drop to more negative values for SPV PCEs than SSW PCEs (Fig. 4b) corresponding to around -15.5 • C for the former and around -13.5 • C for the latter case. In case of SSW PCEs, the T2m index is lower than its mean the entire period before the drop below one standard deviation. While overlapping with the zero line, the confidence interval, noticeably decreases between around Fig. 4 Lagged T2m anomaly index composite (coloured; left-hand yaxis) including 90% confidence interval and T2m absolute value index composite (grey; right-hand y-axis) for SPV PCEs (a) and SSW PCEs (b). Lagged evolution of T2m anomalies meridionally averaged over the study area for SPV PCEs (c) and SSW PCEs (d) 40 and 33 days before the PCE. On average, the T2m stays lower than the PCE threshold for around 15 days, lower than its mean for 40 days, and mostly afterwards.
Further investigating the longitudinal evolution from -60days to +60-days lag, composite T2m anomalies, latitudeweighted and meridionally averaged over the study region, are shown in Fig. 4c-d. Vertical dashed lines mark the zonal borders of the study region. Only significant values are shown. In case of SPV PCEs (Fig. 4c), we find positive anomalies to prevail across EUR for negative lags. Further, from 0-days to +10-days lag, negative T2m anomalies are most dominant in the central part. The eastward tilt of the anomaly pattern may hint at an eastward propagation of the anomalies. Thus, the farther east in the study region, the longer the negative anomalies persist, ranging between less than 10 days on the western side to almost 20 days on the eastern side. In contrast, for SSW PCEs (Fig. 4d), we find negative T2m anomalies in eastern EUR around 60 days and in the central part around 30 to 40 days before the PCEs start. Further, at about 20 days before the PCE onset, significant negative anomalies develop in western EUR. During the PCE, while anomalies in eastern EUR are colder, the ones on the western side are more persistent (up until +20-days to +25-days lag) with a slight eastward propagation. Figure 5 shows SPV PCE composite anomalies for -10days, -5-days, 0-days, +5-days and +10-days lag (left to right) for T2m (first row), SLP (second row), z500 (third row) and 250hPa zonal wind (u250) (fourth row). Only significant anomalies are shown. Positive T2m anomalies in the central eastern part accompanied by negative SLP anomalies in the center of the study area emerge 10 days before the PCE onset (first column). At the same time, a z500 wave train stretches from North Africa to the Fig. 5 Lagged composite T2m (first row), SLP (second row), z500 (third row) u250 (fourth row) anomalies for -10-days, -5-days, 0-days, +5-days and 10-days lag (left to right) relative to the onset of SPV PCEs. Only significant anomalies are shown eastern part of the study region. Further, we find negative (positive) z500 and SLP anomalies over the southern tip of Greenland (western North Atlantic), accompanied by accelerated central North Atlantic u250. Five days before the event (second column), positive T2m anomalies in the study region fade while negative ones develop over its western half. The North Atlantic height anomalies merge with the Rossby wave train found at -10-days lag with positive (negative) z500 centers over the eastern and western (central) part of EUR. At 0-days lag (third column) we see a strong zonal acceleration of the u250 over the North Atlantic and the southern border of the study area. Further, the positive heights over the Arctic Ocean expand with positive SLP coinciding with negative T2m over central northern EUR. Strong positive T2m anomalies develop over the Arctic Ocean in the Barents sea area that fade at +10days lag (fifth column). From the 0-days to +5-days lag (fourth column) we find a slowly eastward propagating wave. Together with positive SLP over EUR, the negative T2m anomalies in the study region strengthen and expand at +5-days lag. A zonal shift of T2m at +10-days lag, covering particularly the central and eastern part of EUR, reflect the eastward tilt identified in Fig. 4c.
For SSW PCEs, the pattern evolution is shown in Fig. 6. At -10-days lag (first column), we do not find significant T2m anomalies in most parts of EUR except for the far west. The latter negative anomalies are accompanied by positive (negative) SLP and z500 anomalies over Iceland (central Europe). Five days later (second column), negative T2m anomalies develop north of the northern border as well as the western part of the study region. Positive SLP and z500 anomalies strengthen and extend further toward EUR. Additionally we find a positive z500 anomaly over the Ural region as well as two negative ones on the most western and most eastern sides of EUR. Together with negative SLP anomalies over the central North Atlantic, both negative SLP and z500 anomalies over central Europe and eastern EUR as well as positive z500 over the Arctic Ocean strengthen at 0-days lag and are accompanied by a southward shift of the jet (third column) that coincides with negative z500 anomalies over EUR. Concurrent T2m anomalies stretch across northern EUR before moving further south, covering almost the entire study region 5 days later (fourth column) with particularly strong values over its central part. Concomitantly, a negative NAO signature in SLP and z500 together with a southward shift of the jet stream is visible over the North Atlantic at +5-days and +10-days lag (fourth and fifth column). The latter can be additionally observed over eastern EUR. With increasing lag, the positive SLP anomalies together with the center of strongest negative T2m anomalies migrate further toward the eastern part of Eurasia.
Concluding, in both cases, the T2m anomalies affect the whole study region but set in more abruptly in the SPV PCE case. Further, a strong positive quasi-stationary z500 anomaly with a barotropic structure develops over the Arctic, yielding a southward shift of the jet over Eurasia and thus blocking warm air advection from the South, thus facilitating the PCE. In the SSW PCE case, Ural blocking is present before the PCE onset and the Arctic anticyclonic anomaly is connected to positive z500 anomalies over Greenland related to the negative NAO phase. In contrast, co-located with warm T2m anomalies, the positive z500 anomaly over the Arctic is included in a zonally propagating Rossby wave train originating in the North Atlantic in the SPV PCE case.

The Stratospheric Link
In the next step we want to see how and if the stratosphere can be linked to the PCEs. On that account, composites of latitude-weighted geopotential heights and temperature anomalies averaged over all grid points between 60 • N and 75 • N in the former as well as north of 70 • N in the latter case are shown in Fig. 7a-d.
For SPV PCEs (Fig. 7a) we find significant negative height anomalies in the upper stratosphere from around -25-days lag. A downward migration of the signal is observed but the tropopause is not crossed before about 30 days after the PCE. However, the signal reaches the approximate tropopause level (horizontal line) close to 0days lag. Similarly, we see downward progressing cold polar cap temperature anomalies (Fig. 7c) from around -35-days to 0-days lag that pass the tropopause only at around -20days lag. Interestingly, we find positive anomalies in the troposphere and surface around the onset of the SPV PCEs that likely reflect the near surface warming over the Barents sea ( Fig. 5; first row).
In case of SSW PCEs (Fig. 7b), significant positive geopotential height anomalies are observed from about -20days lag. A downward migration of this signal is clearly visible, reaching the surface at about -20-days lag as well as between -5-days and +10-days lag. This is consistent with the development of the negative phase of the NAO, visible in Fig. 6. Additionally, a warming over the polar cap, identified at around -20-days lag, progresses downward into the troposphere but without significant connection to the surface (Fig. 7d).
In the right constellation, the lower stratosphere can be in a reflective state, i.e. instead of being absorbed, upward traveling planetary waves are reflected back to the troposphere where they accumulate and influence the circulation. Following Nath et al. (2014), we calculated the reflective index of the stratosphere as a function of longitude and time (Fig. 7e-f) with negative values indicating that reflection is possible. Only significant values are shown. In the SPV PCE case (Fig. 7e), the stratosphere is significantly non-reflective from across most longitudes from -20-days to +20-days lag except between around 100 • W to 140 • W. Across the latter longitudes, the stratosphere is dominantly reflective up until +5-days lag. In contrast, in case of SSW PCEs (Fig. 7f), we only find significant positive values across most longitudes and lags indicating prominent wave absorption. However, before the PCE onset, values are dominantly non-significant.
To further investigate upward and downward propagation of planetary waves, Fig. 8 shows the latitude-weighted and then averaged between 58 • N and 74 • N vertical WAF component (shading) as well as WAF in the x-z domain (arrows) across longitudes for -10-days, -5-days, 0-days, +5-days and +10-days lag (from top to bottom) for SPV PCEs (left column) and SSW PCEs (right column). The vectors are scaled with 15 and 0.5 in the x and z direction, respectively. Consistent with the reflective index, we indeed find downward reflected waves between about 120 • W to around 55 • W for -10-days to +5-days lag for the SPV PCE case that converge with upward directed waves in the mid-to upper troposphere, thereby piling up energy and thus potentially influencing the North Atlantic circulation. A minor downward WAF at similar location but spatially limited is found for SSW PCEs at -5-days lag. In both cases, upward wave propagation into the stratosphere is prominent over eastern EUR and the Pacific.
The horizontal distribution of absolute values and anomalies of the 3D WAF with vertical component as shading and horizontal component (x-y) as arrows at 500hPa (first and second column) and 150hPa (third and fourth column) composited for SPV PCEs at -10-days, -5-days, 0days, +5-days and +10-days lag (top to bottom) is shown in Fig. 9. for -60-days to +60-days lag. Horizontal black line marks the approximated tropopause level. Reflective index U r for SPV PCEs (e) and SSW PCEs (f) for -20-days to +30-days lag  At -10-days to -5-days lag, stationary waves dominantly propagate up-and, between 30 • N and 60 • N southeastward over East Asia, the North Pacific and western North America in both heights, as well as over the North Atlantic in 500hPa, which is particularly more (less) than usual over the Gulf of Alaska (Central North Atlantic and eastern Siberia/EUR) at -5-days lag. Western EUR experiences unusually strong upward wave propagation in 500hPa from -5-days to +5-days lag. Additionally, we find an eastward flux over EUR, the Aleutians and several zonal stretches  Fig. 9 but for SSW PCEs between 60 • N and 75 • N for both heights and throughout the analysed lags. Interestingly, the North Atlantic eastward flux is unusually weak from -5-days to +5-days lag which may hint at slower than usual eastward propagating Waves. At -5-days, we see unusually strong eastward fluxes, locally confined over the Barents sea, towards the region covered by the anticyclonic anomaly. In 150hPa, negative WAF anomalies, i.e. downward propagating waves, are found over northern Canada and Greenland between -10-days lag to +5-days lag. At 0-days lag, unusually strong down-and eastward directed North American WAF appears in 500hPa, seemingly part of a wave train with alternating positive and negative vertical WAF anomalies over the North Atlantic. Concurrently, upward directed anomalous WAF over central (150hPa) and western (500hPa) EUR strengthens. Figure 10 shows the same as Fig. 9 but for SSW PCE events. At -10-days to -5-days lag, we find unusual upward propagating stationary wave over the central east Pacific, and, different to the SPV PCE case, over the North Atlantic in both heights. At the same time, southeastward directed WAF are found over the North Atlantic and central east Pacific while eastward wave propagation is found between about 55 • N to 75 • N and 90 • E to 90 • W. At -5-days and 0-days lag, we find negative (positive) vertical WAF anomalies at 150hPa (500hPa) over western EUR. From 0-days to +5-days lag, upward directed WAF over eastern EUR is unusually weak in both heights. Simultaneously, upward directed WAF over the central (both heights) and western (500hPa) North Pacific is unusually weak. Positive (negative) WAF anomalies are found at 500hPa (150hPa) over the eastern North Atlantic at +5-days lag (+5-days and +10-days lag). Simultaneously, eastward wave propagation over Northern Europe is unusually weak. At +10-days lag positive vertical WAF anomalies over Canada and negative anomalies over central EUR emerge at 150hPa while at 500hPa weak (unusually strong) upward wave propagation over Southern Greenland (central Europe and Central to eastern North America) appears.
Concluding, in the SPV PCE case, we find a reflective stratosphere over the eastern US consistent with downward high-latitude wave propagation converging in the midtroposphere as well as downward propagation over the central North Atlantic of which the latter is co-located with positive height anomalies included in the wave train slowly propagating toward EUR and unusually strong eastward wave fluxes over the Barents sea before the PCE onset. For SSW PCEs we find the upper stratosphere to be mainly nonreflective but with downward wave propagation from the lower stratosphere that is at similar location but weaker as in the SPV PCE case.

Discussion and Conclusion
The negative impacts of PCEs make their early prediction, preferably on S2S-timescales, a necessity researchers are striving for. The study region EUR is chosen as it overlaps with the area with highest explained variability of T2m across the Northern Hemisphere based on an EOF analysis of daily winter-time data. Central EUR T2m extremes show a decreasing trend, i.e. more than -0.1 • /yr, based on QR using low quantiles that cannot be identified using the median. This trend, together with a concomitant positive trend over the Barents/Kara sea is consistent with previous studies linking sea ice loss and related reduction in meridional temperature gradient (e.g. Outten and Esau 2012) to blocking over and cold air advection toward EUR (e.g. Mori et al. 2014;Yao et al. 2017).
The T2m index, as a latitude-weighted average across EUR, is negatively skewed, indicating the recurrent presence of extreme cold events. Based on the stratospheric state in the 40 days before their onset, we find PCEs to occur during the full range of stratospheric variability, but SPV PCEs dominate in December/January and SSW PCEs in February/March. Conspicuously, 8 out of 11 SPV PCEs start more than 20 days after the initial stratospheric signal in 10hPa while 5 out of 9 SSW PCEs start no later than 11 days after. In both cases, the T2m anomalies affect the whole study region, but SPV PCEs are colder and set in more rapidly after a period of above average T2m. In contrast, EUR is exposed to unusual cold T2m before SSW PCEs set in. Consistent with the AO response given in (e.g. Kim and Son 2016), coldest anomalies are located along the northern border of EUR during the onset of SSW PCEs, which then migrate eastward. Negative anomalies, however, are most persistent in western EUR. Contrariwise, SPV PCEs endure longest in eastern EUR during their evolution.
In both cases, a strong positive, quasi-stationary z500 anomaly with a barotropic structure over the Barents/Kara sea and northern Siberia and a downstream negative anomaly in mid-latitudes yields a southward shifted jet stream over EUR, thereby, inhibiting (promoting) warm (cold) air advection from the South (North). The remote link to the North Atlantic, however, differs across SSW and SPV PCEs. In the former case, we find this anticyclonic anomaly to be connected to positive heights over Greenland as part of the negative NAO phase progressing downward from the stratosphere, reaching the surface 5 days before the PCE onset. Simultaneous unusually weak eastward wave fluxes over Northern Europe and Ural blocking are likely to be involved in generating the cold temperatures. In contrast to absorbing upward traveling planetary waves in the SSW PCE case, the stratosphere is reflective over the North American sector before SPV PCE events. Stratospheric reflection is consistent with geopotential height and polar cap temperature anomalies not progressing into the troposphere (Perlwitz and Harnik 2004). Upward propagating planetary waves over the Aleutians are reflected toward the western North Atlantic, where they converge with upward propagating waves in the mid-troposphere. Subsequent accumulation of wave energy potentially influences the tropospheric circulation (e.g. Perlwitz and Harnik 2004). We further find locally confined, unusually strong eastward wave fluxes into the Barents sea anticyclone. The latter is connected to a slowly propagating zonal Rossby wave train whose North Atlantic origin is approximately co-located with negative vertical WAF anomalies. While local, tropospheric processes can generate the negative height anomalies over the Barents/Kara sea (Inoue et al. 2012), the above results suggest a remote influence in the sense of stratospheric wave reflection toward the North Atlantic and a subsequent modulation of z500 anomalies upstream of the Arctic anticyclone that act to reinforce the latter. Previous studies congruently linked the North Atlantic circulation with blocking over the Eurasian sector (e.g. Luo et al. 2016;Chen et al. 2021;Murto et al. 2022). While differing in exact location but consistent with our study, Luo et al. (2016) and Chen et al. (2021) link North Atlantic negative height anomalies and intensified wind speeds, as found in the present study, to the slowly eastward propagation of a Rossby wave train including an anticyclonic anomaly over EUR. Further, we find concomitant warming over the Barents Sea region and weakening of the mean westerly flow over northern EUR to be coherent with findings by Yao et al. (2017), who connect a reduced meridional temperature gradient to decreased Eurasian wind speeds and thus increased stationarity and persistence of the anticyclonic anomaly that is linked to EUR cold extremes.
Concluding, SSW PCEs and SPV PCEs have distinct and common local characteristics and remote influences. In both cases, the ascertained link to winter-time variability of the stratospheric polar vortex provides potential to enhance their extended-range predictability.