1. Introduction
The Arctic is warming faster than the rest of the globe (Schneider and Held 2001; Screen and Simmonds 2010; Serreze and Barry 2011). Arctic sea ice cover has undergone unprecedented declines (Comiso et al. 2008; Stroeve et al. 2012b), and the subsurface heat stored in the Arctic Ocean has increased in the past decade (Timmermans et al. 2018). The remote influence of these changes is a pressing topic, in large part because of a possible link to recent episodes of extreme weather in the midlatitudes (e.g., Liu et al. 2012; Tang et al. 2013; Petoukhov et al. 2013). One suggestion for how the Arctic affects the midlatitudes is through the stratosphere: sea ice loss is thought to lead to conditions that favor the upward propagation of planetary waves into the stratosphere, which produce a weakening of the stratospheric polar vortex (Garfinkel et al. 2010; Kolstad et al. 2010; Cohen et al. 2014). The polar vortex signal can then be communicated back down to the troposphere to create a negative North Atlantic Oscillation (NAO) response (Baldwin and Dunkerton 2001; Polvani and Waugh 2004), along with its associated surface climate effects (Kolstad et al. 2010; Scaife et al. 2005; Ambaum and Hoskins 2002; Thompson and Wallace 1998).
A number of studies have investigated the existence of this stratospheric pathway linking Arctic sea ice loss to a midlatitude circulation response, but the results are not clear-cut. There is evidence that, in the present climate, fall sea ice loss is associated with a weaker polar vortex in early winter and a negative northern annular mode (NAM) response in late winter to early spring (Jaiser et al. 2013; Feldstein and Lee 2014; Cohen et al. 2014; García-Serrano et al. 2015; King et al. 2016; Yang et al. 2016; Peings and Magnusdottir 2014; Nakamura et al. 2015). Some modeling studies support these findings (Kim et al. 2014; Nakamura et al. 2016), while others report a strengthening of the polar vortex in response to sea ice loss either earlier or later during the cold season (Cai et al. 2012; Scinocca et al. 2009; Screen and Simmonds 2013; Blackport and Screen 2019). These discrepancies could partly be explained by differences in details of how the sea ice cover is reduced in the model experiments, as it is only ice loss in the Atlantic–Barents Sea sector that seems able to weaken the polar vortex (Screen et al. 2018; Sun et al. 2015; Screen 2017; Zhang et al. 2017). Biases in the background state of models could also play a role. For example, both the strength of the stratospheric response and the sign of the NAO response to sea ice loss can depend on the climatological latitude of the jet stream (Smith et al. 2017).
Furthermore, there are open questions about the exact processes connecting sea ice loss to the polar vortex. It is well known that planetary waves originating in the troposphere and propagating upward into the stratosphere can weaken the polar vortex (Matsuno 1971). In fact, breakdowns of the polar vortex are often preceded by short-term poleward heat flux anomalies, or upward wave activity flux anomalies (Polvani and Waugh 2004; Hitchcock et al. 2013a; Cohen and Jones 2011; Garfinkel et al. 2010). Sea ice loss in the Atlantic–Barents Sea sector is thought to enhance high-latitude upward wave activity fluxes by producing tropospheric circulation anomalies that facilitate upward propagation of waves (Jaiser et al. 2013; Kim et al. 2014; Feldstein and Lee 2014; Sun et al. 2015). Some models reproduce aspects of this stratospheric pathway in response to sea ice loss (Nakamura et al. 2016; Hoshi et al. 2017; Zhang et al. 2017). However, other models find that sea ice loss leads to weaker upward fluxes because the weakened equator-to-pole temperature gradient weakens the baroclinic eddy source (Seierstad and Bader 2009), which is a consistent feature across coupled models (Screen et al. 2018; Smith et al. 2017).
Given these seemingly conflicting lines of evidence, it is legitimate to ask whether fall sea ice loss, or associated surface heat fluxes, can in general modulate the large-scale atmospheric circulation to activate the stratospheric pathway in a consistent and detectable way. And if such a stratospheric bridging mechanism is in play, what is its timing and signature? Why is it that models do not agree on the sign of the response? In reality and in realistic experimental setups, it is difficult to isolate the influence of sea ice on the stratosphere from the influence of other factors (Overland et al. 2016), such as forcing from the tropical Pacific (Feldstein and Lee 2014); substantial noise from internal variability (Sun et al. 2015; Blackport and Kushner 2017); long stratospheric persistence (Shepherd 2000); changes in snow cover and sea surface temperature (Cohen et al. 2014); and interference from winter sea ice loss, which has been shown to dominate the midwinter circulation response in some models (Sun et al. 2015). It is also difficult to separate causes and effects, as, for example, transient sea ice loss may be caused by low-frequency atmospheric circulation anomalies, which then are difficult to disentangle from the atmospheric response to sea ice loss (Blackport et al. 2019). These complications greatly reduce the signal-to-noise ratio of the sea ice influence on the atmosphere and hamper analyses and interpretations (Overland et al. 2016).
Idealized process studies offer a way to clarify the atmospheric response to sea ice loss. The stratosphere responds to Arctic forcing in simplified model setups (Wu and Smith 2016; Zhang et al. 2017), but there are outstanding questions regarding the time scale and dynamics of the adjustment in the context of large internal variability and natural persistence in the stratosphere. Here, we use idealized general circulation model (GCM) simulations to study the atmospheric response to short-term Arctic perturbations designed to represent the transient forcing due to fall sea ice loss. A large ensemble of simulations allows us to estimate mean time scales and amplitudes of the response, and to identify robust mechanisms.
2. Method
We use the idealized GCM described in O’Gorman and Schneider (2008) and Frierson et al. (2006). The model solves the primitive equations at a spectral resolution of T85 (i.e., with a transformed grid spacing of 1.4° at the equator), with 30 unevenly spaced σ-coordinate levels and the highest level at σ = 4.2 × 10−5. The GCM uses a two-stream gray radiative transfer scheme, in which the radiative flux is only a function of temperature (i.e., there are no cloud radiative or water vapor feedbacks). The radiative forcing at the top of the atmosphere is an approximation of annual-mean insolation, so that variations of the zonal-mean jet position only arise by internally generated variability. The longwave optical thickness τ of the atmosphere has components that depend linearly and quadratically on pressure, and it is proportional to sin(ϕ)2, where ϕ is latitude, with an optical thickness of 7.2 at the equator and 3 at the pole (Frierson et al. 2006; O’Gorman and Schneider 2008). This produces a tropospheric jet in each hemisphere similar to that seen in Northern Hemisphere winter (Fig. 1), despite the annual-mean insolation. The lack of diabatic heating from shortwave absorption by ozone leads to a lack of separation of the tropospheric jet and the stratospheric polar vortex (Fig. 1), but the tropospheric mean state is in good agreement with the observed winter zonal mean. The annular mode is the dominant mode of variability in the stratosphere and troposphere of the model, and its variability is similar to that seen in more complex models (Figs. S1–S3 in the online supplemental material). A better representation of the stratospheric mean state could have been obtained by adjusting the stratospheric radiative forcing (Kushner and Polvani 2004), which may change the dominant adjustment mechanisms to episodic surface forcing.

(a) Zonal-mean zonal wind in idealized GCM. The vertical axis is a logarithmic σ coordinate axis (equal increments are approximately equal altitude increments). The gray line indicates the tropopause. (b) Standard deviation of zonal-mean zonal wind between 50° and 75°N.
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1

(a) Zonal-mean zonal wind in idealized GCM. The vertical axis is a logarithmic σ coordinate axis (equal increments are approximately equal altitude increments). The gray line indicates the tropopause. (b) Standard deviation of zonal-mean zonal wind between 50° and 75°N.
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
(a) Zonal-mean zonal wind in idealized GCM. The vertical axis is a logarithmic σ coordinate axis (equal increments are approximately equal altitude increments). The gray line indicates the tropopause. (b) Standard deviation of zonal-mean zonal wind between 50° and 75°N.
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
Earth’s surface is represented by a 1-m-thick slab ocean. This shallow representation of a mixed layer results in a shorter decay time scale of surface energy budget anomalies compared to observations. Turbulent surface fluxes at the lower atmospheric boundary are represented by the standard bulk aerodynamic formulas and Monin–Obukhov similarity theory. A tropical ocean energy flux divergence (Q flux) with an amplitude of 60 W m−2 at the equator is added to mimic the equatorial ocean circulation, which helps to obtain a more realistic Hadley cell and tropospheric jet position (Bordoni and Schneider 2008; Levine and Schneider 2011).
a. Representation of Arctic sea ice anomalies and short-term ocean heat loss
Typically, model studies either prescribe sea ice anomalies [see the overview in Gao et al. (2015)], which artificially fixes part of the surface energy balance, or alter the radiative fluxes to remove sea ice (Screen et al. 2018), which still tampers with the surface energy balance but at least allows for an energetically consistent adjustment of the climate system to the perturbation. A more direct approach is to perturb the surface energy balance itself by imposing Q fluxes that mimic the effect of sea ice variability on the surface heat flux (Fig. 2). There are still some complications here because a substantial portion of surface heat flux variability is not due to sea ice changes (Walsh and Johnson 1979; Fang and Wallace 1994; Deser et al. 2010; Woods and Caballero 2016; Sorokina et al. 2016), but rather drives sea ice changes (Smedsrud et al. 2013; Lee et al. 2017; Timmermans et al. 2018; Blackport et al. 2019). But the approach does allow for more control in the experimental setup, since it circumvents discussions about local ice–atmosphere feedbacks.

Composite of daily net ocean cooling in the Barents–Kara Sea (1979–2012). Gray lines indicate fluctuations around the seasonal cycle in fall season (October–January) for a given year. Each line is repositioned to the day (day 0) when upward fluxes first exceed two climatological standard deviations (orange line; ~50 W m−2), which typically occurs in October or November. Thin blue lines are the 3 years of strongest ocean cooling, and the thick black line is the mean of these years. The data are derived from ERA-Interim (Dee et al. 2011) by first subtracting the seasonal cycle and then using a 5-day running mean.
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1

Composite of daily net ocean cooling in the Barents–Kara Sea (1979–2012). Gray lines indicate fluctuations around the seasonal cycle in fall season (October–January) for a given year. Each line is repositioned to the day (day 0) when upward fluxes first exceed two climatological standard deviations (orange line; ~50 W m−2), which typically occurs in October or November. Thin blue lines are the 3 years of strongest ocean cooling, and the thick black line is the mean of these years. The data are derived from ERA-Interim (Dee et al. 2011) by first subtracting the seasonal cycle and then using a 5-day running mean.
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
Composite of daily net ocean cooling in the Barents–Kara Sea (1979–2012). Gray lines indicate fluctuations around the seasonal cycle in fall season (October–January) for a given year. Each line is repositioned to the day (day 0) when upward fluxes first exceed two climatological standard deviations (orange line; ~50 W m−2), which typically occurs in October or November. Thin blue lines are the 3 years of strongest ocean cooling, and the thick black line is the mean of these years. The data are derived from ERA-Interim (Dee et al. 2011) by first subtracting the seasonal cycle and then using a 5-day running mean.
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
Here we prescribe an additional Q flux to mimic temporary sea ice anomalies in the Arctic and the associated perturbation of the surface energy budget. We use a Q-flux convergence/divergence imposed for 20 days (between day −24 and day −5) in an annulus between 70° and 85°N for each ensemble member
The polar boundary layer (BL) stratification affects the upward propagation of perturbations at the surface. A weakening of the BL stratification is observed as part of the long-term trend in the Arctic (Screen and Simmonds 2010). This is expected to continue in the future (Deser et al. 2010). In our GCM, the near-surface static stability in the polar regions is larger than in midlatitudes (Schneider and O’Gorman 2008), but it is still weaker than typical Arctic winter conditions (Dee et al. 2011).
We choose the amplitude of the imposed polar Q fluxes based on an order-of-magnitude analysis of observed surface energy flux anomalies in the Barents–Kara Sea in fall (Fig. 2). The precise value turns out to be relatively unimportant because, as we will discuss, the atmospheric response to surface energy budget perturbations is approximately linear in the perturbation amplitude.
b. Smoothing and output frequency
The model output is sampled as 5-day means and is smoothed between adjacent 5-day blocks. Data at a given time are the centered mean over 10 days. The origin (0) of the time axis is 5 days after the imposed Q flux ends, that is, when no contribution of the Q flux is contained in the centered 10-day mean. The 10-day means centered on day −5 and day −25 contain contributions from the imposed Q fluxes at the beginning and end of the forcing period.
c. Ensemble runs
To isolate the transient response to the polar surface perturbation, a set of N = 50 initial conditions was created from an unperturbed climatological run by saving restart files every 60 days after 1440 days of spinup. From each of the N initial conditions, a pair of ensemble members was started, one (

Estimate of sufficient ensemble member size Nmin for characteristic atmospheric variables by inverting the Student’s t test for differences of means as in Screen et al. (2013). (a) Surface pressure, (b) midlatitude tropopause geopotential height, and (c) midlatitude zonal wind. Results are for members with initial conditions spaced 30 days apart.
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1

Estimate of sufficient ensemble member size Nmin for characteristic atmospheric variables by inverting the Student’s t test for differences of means as in Screen et al. (2013). (a) Surface pressure, (b) midlatitude tropopause geopotential height, and (c) midlatitude zonal wind. Results are for members with initial conditions spaced 30 days apart.
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
Estimate of sufficient ensemble member size Nmin for characteristic atmospheric variables by inverting the Student’s t test for differences of means as in Screen et al. (2013). (a) Surface pressure, (b) midlatitude tropopause geopotential height, and (c) midlatitude zonal wind. Results are for members with initial conditions spaced 30 days apart.
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
To test whether the responses
The climatology and standard deviation at each grid point are calculated from a 3600-day-long unperturbed control integration, started from the initial condition of the first member (after 1440 days of spinup). Figure 1 shows the climatological-mean zonal wind
3. Linearity and symmetry
We first analyzed to what extent the response is linear in the forcing by comparing simulations with positive forcing (Q-flux convergence) and negative forcing (Q-flux divergence) of different amplitudes. Figure 4 shows the zonal-mean responses

Northern annular mode at σ = 0.1 for different ensembles. (a) Red lines represent the ensemble-mean response
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1

Northern annular mode at σ = 0.1 for different ensembles. (a) Red lines represent the ensemble-mean response
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
Northern annular mode at σ = 0.1 for different ensembles. (a) Red lines represent the ensemble-mean response
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
The signal-to-noise ratio of the responses
We tested to what degree zonal asymmetries in forcing shape the response by comparing simulations with the same zonal-mean forcing (66.6 W m−2), but concentrating it to one-third of a latitude circle, that is, using a Q-flux convergence of 200 W m−2 in a 120°-wide longitude sector. The dashed red and blue lines in Fig. 4a show the NAM response in these experiments with positive and negative forcing, and the dashed black line shows the symmetrized response. Within statistical uncertainties, these responses are indistinguishable from the response to the zonally symmetric forcing (Figs. 4a,b). We found similar results for other variables, such as lower-stratospheric zonal wind. Stationary waves, which have been posited to be important on the basis of other modeling studies (Smith et al. 2010, 2011; Hitchcock et al. 2013b) are absent in our model setup. Thus, our results suggest that they are not essential for the response to high-latitude warming in our simulations. However, they may still play a role in the real atmosphere, where stationary waves are prominent in the Northern Hemisphere.
4. Thermal and eddy-driven response
a. Direct thermal response to zonally symmetric heating

Zonal-mean response
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1

Zonal-mean response
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
Zonal-mean response
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1

As in Fig. 5b, but normalized by the standard deviation at each height estimated from the control simulation (Fig. 1b).
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1

As in Fig. 5b, but normalized by the standard deviation at each height estimated from the control simulation (Fig. 1b).
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
As in Fig. 5b, but normalized by the standard deviation at each height estimated from the control simulation (Fig. 1b).
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
b. Eddy-driven stratospheric response
While the initial stratospheric response is a result of thermal-wind adjustments to changes in surface temperature gradients, the persistence of the delayed response needs to be understood in the context of eddy–mean flow interactions.
Figure 7 shows the time and height dependencies of the zonal-mean (Eulerian) meridional velocity

Time slices of the symmetric response
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1

Time slices of the symmetric response
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
Time slices of the symmetric response
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
While the residual circulation captures eddy mass transport and hence gives a better description of tracer transport than the Eulerian mean circulation, here we focus on the Eulerian mean because it dominates the response of the circulation to the forcing (cf. the residual circulation response in appendix B). The Eulerian mean stratospheric circulation responds to the surface forcing primarily by weakening, especially after cessation of the surface forcing (e.g., days 31–60; see left two columns of Fig. 7). That is, there is anomalous southward flow in the subtropics and anomalous northward flow in midlatitudes for about 2 months after the surface forcing has ended. By continuity, this is compensated by anomalous upward flow between around 30° and 50°N and downward flow poleward of 50°N.

(a) Response of vertical wave activity flux just above the midlatitude tropopause between 45° and 75°N at σ = 0.1 (red). (b) Response of the meridional wave activity flux at 45°N between σ = 5 × 10−3 and σ = 0.1 (blue). (c) Response of the meridional wave activity flux divergence (eddy momentum flux convergence) in the lower midlatitude stratosphere between 45° and 75°N and between σ = 5 × 10−3 and σ = 0.1. The gray dashed line shows the climatological value and the colored shading the 95% confidence intervals of the respective variables. The red bars between days −25 and 0 indicate the forcing time. We show the symmetric response
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1

(a) Response of vertical wave activity flux just above the midlatitude tropopause between 45° and 75°N at σ = 0.1 (red). (b) Response of the meridional wave activity flux at 45°N between σ = 5 × 10−3 and σ = 0.1 (blue). (c) Response of the meridional wave activity flux divergence (eddy momentum flux convergence) in the lower midlatitude stratosphere between 45° and 75°N and between σ = 5 × 10−3 and σ = 0.1. The gray dashed line shows the climatological value and the colored shading the 95% confidence intervals of the respective variables. The red bars between days −25 and 0 indicate the forcing time. We show the symmetric response
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
(a) Response of vertical wave activity flux just above the midlatitude tropopause between 45° and 75°N at σ = 0.1 (red). (b) Response of the meridional wave activity flux at 45°N between σ = 5 × 10−3 and σ = 0.1 (blue). (c) Response of the meridional wave activity flux divergence (eddy momentum flux convergence) in the lower midlatitude stratosphere between 45° and 75°N and between σ = 5 × 10−3 and σ = 0.1. The gray dashed line shows the climatological value and the colored shading the 95% confidence intervals of the respective variables. The red bars between days −25 and 0 indicate the forcing time. We show the symmetric response
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
Figure 9 confirms that the polar stratospheric dynamic warming is primarily a result of the dynamic warming associated with anomalous subsidence; dynamic warming associated with eddy fluxes and meridional advection plays a lesser role (see appendix B). The static stability

(a) Thermodynamic budget in the polar stratosphere (75°–90°N, σ = 5 × 10−3 to σ = 0.1), with the total temperature tendency (red), the subsidence warming [right-hand side of (5); blue], and total dynamic warming including meridional advection and eddy fluxes (dashed blue; appendix B). (b) Vertical wind in the polar stratosphere [within the same box as in (a)]. Shading shows 95% confidence intervals. We used the symmetric response
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1

(a) Thermodynamic budget in the polar stratosphere (75°–90°N, σ = 5 × 10−3 to σ = 0.1), with the total temperature tendency (red), the subsidence warming [right-hand side of (5); blue], and total dynamic warming including meridional advection and eddy fluxes (dashed blue; appendix B). (b) Vertical wind in the polar stratosphere [within the same box as in (a)]. Shading shows 95% confidence intervals. We used the symmetric response
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
(a) Thermodynamic budget in the polar stratosphere (75°–90°N, σ = 5 × 10−3 to σ = 0.1), with the total temperature tendency (red), the subsidence warming [right-hand side of (5); blue], and total dynamic warming including meridional advection and eddy fluxes (dashed blue; appendix B). (b) Vertical wind in the polar stratosphere [within the same box as in (a)]. Shading shows 95% confidence intervals. We used the symmetric response
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
In summary, the delayed response of the zonal wind in the stratosphere arises from eddy-driven changes in the stratospheric circulations. They, in turn, flow from the changes in the tropospheric circulation associated with reduced meridional temperature gradients and reduced vertical wave activity propagation. The crucial reason why the stratospheric response is delayed relative to that of the troposphere is that thermal damping time scales in the stratosphere are longer than in the troposphere.
5. Decay and dependence on the mean flow
a. Decay of stratospheric anomalies and effect on troposphere
Stratospheric zonal wind anomalies respond on two time scales. There is a fast adjustment to surface temperature anomalies through thermal-wind balance, (2), and a slower eddy-mediated response involving an overturning circulation and dynamic heating and cooling, (5). In the ensemble mean, the fast response is evident in the surface warming signal, which decays within a month (Fig. 5c and red contours in Figs. 7c and 7f). The eddy-mediated lower-stratospheric warming lasts at least 50 days, along with the weakening of the stratospheric zonal winds, and these signals decay gradually over radiative time scales (Fig. 7i). The amplitude of the lower-stratospheric warming depends mainly on the integrated amount of heat input from the surface and less on the heating rate (see appendix A).
Stratosphere-to-troposphere linkages beyond 60 days, as seen in Baldwin and Dunkerton (2001), are not clearly recognizable in the ensemble mean in our simulations (Figs. 5b,c). This is because the phases of the downward-propagating signals from the stratosphere have a random component and are weak compared to sudden stratospheric warmings, making it difficult to identify coherent tropospheric signals against the background of internal variability.
b. Dependence on mean flow
To understand the ensemble spread of the circulation response to heating, we investigate how the response depends on the mean flow of each ensemble member.
Figure 10 shows lagged regressions of the potential temperature response in the polar cap (75°–90°N) against the fast stratospheric zonal wind response (days 1–30), both normalized by their climatological standard deviations at each level. Positive values of the regression coefficient indicate a polar cap warming with a weakening of the stratospheric winds, that is, positive potential temperature anomalies associated with negative zonal wind anomalies. We pool ensemble members

Regression of polar cap (75°–90°N) potential temperature response S on the earlier stratospheric zonal wind response in the green box (σ = 10−3–10−2, 50°–75°N, days 1–30). Each symmetric response Si is the anomaly from the control member Ci, so a possible dependence of downward propagation on the control simulation variability is not captured in this regression. Positive values indicate warmer potential temperature with a weakened polar vortex. Stippling shows points that are significantly different from zero (5% level).
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1

Regression of polar cap (75°–90°N) potential temperature response S on the earlier stratospheric zonal wind response in the green box (σ = 10−3–10−2, 50°–75°N, days 1–30). Each symmetric response Si is the anomaly from the control member Ci, so a possible dependence of downward propagation on the control simulation variability is not captured in this regression. Positive values indicate warmer potential temperature with a weakened polar vortex. Stippling shows points that are significantly different from zero (5% level).
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
Regression of polar cap (75°–90°N) potential temperature response S on the earlier stratospheric zonal wind response in the green box (σ = 10−3–10−2, 50°–75°N, days 1–30). Each symmetric response Si is the anomaly from the control member Ci, so a possible dependence of downward propagation on the control simulation variability is not captured in this regression. Positive values indicate warmer potential temperature with a weakened polar vortex. Stippling shows points that are significantly different from zero (5% level).
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
The lagged regression analysis shows that more negative anomalies of stratospheric zonal wind lead longer-lasting and more downward-propagating positive polar temperature anomalies (i.e., positive regression coefficients). The regression coefficients reach maximum amplitude in the troposphere and near the surface about 30 days after the early lower-stratospheric zonal wind response (Fig. 10). The early wind response itself depends on the combined effect of the forcing and its (nonlinear) interaction with the background stratospheric zonal wind in any given ensemble member. Members with a weaker stratospheric vortex thus have an increased likelihood of producing significant tropospheric responses (Baldwin and Dunkerton 2001; Hitchcock et al. 2013b). However, the nonlinear interactions also affect whether or not signals propagate downward from the stratosphere, and these are subject to variability generated within individual members.
To examine the sensitivity of the delayed zonal-wind response on the mean-flow sampled by internal variability, we regress the mean zonal wind in the control simulation Ci against the late (negative) zonal wind response Si in the ensemble pair branched off from that control simulation. As a measure of the late response, we use the zonal-wind response in the lower stratosphere averaged over the space–time box spanning 50°–75°N and days 50–75, again normalized by the climatological standard deviation at each level. To establish whether the resulting regression coefficients are significant, we compute a bootstrap PDF of a corresponding regression within the control runs alone, regressing the mean zonal wind in the control simulations Ci against the control zonal wind in the lower-stratospheric space–time box. Significance is then established as in section 5a.
Figure 11 shows the resulting regression coefficients and their significance. The strength of the delayed wind response (green box) is related to the strength of the midlatitude flow during the forcing (days −25 to −5): members that respond to polar surface warming with more weakening of the stratospheric wind between days 50 and 75 tend to be initialized from control states with stronger tropospheric winds. The reason appears to be regression toward the mean, resulting from the model’s natural relaxation toward the mean state. The tropospheric midlatitude wind is related to the meridional temperature gradient in the lower troposphere. Anomalously strong meridional winds are associated with anomalously strong temperature gradients, or, assuming fixed subtropical temperatures, an anomalously cold polar troposphere. The tropospheric temperature at the pole and the meridional gradient naturally relax to the mean state. The natural relaxation toward the mean adds to the response to the forcing.

Correlation of midlatitude zonal wind (50°–75°N) in the control ensemble with negative zonal wind response in a space–time box (50°–75°N, σ = 0.1 to σ = 5 × 10−3, days 50–75) of the forced ensemble S. Positive values indicate positive zonal wind anomalies in C with a negative zonal wind anomaly in the box in S. Stippling indicates that the correlation is significantly larger than the correlation in a control ensemble at the 5% level.
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1

Correlation of midlatitude zonal wind (50°–75°N) in the control ensemble with negative zonal wind response in a space–time box (50°–75°N, σ = 0.1 to σ = 5 × 10−3, days 50–75) of the forced ensemble S. Positive values indicate positive zonal wind anomalies in C with a negative zonal wind anomaly in the box in S. Stippling indicates that the correlation is significantly larger than the correlation in a control ensemble at the 5% level.
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
Correlation of midlatitude zonal wind (50°–75°N) in the control ensemble with negative zonal wind response in a space–time box (50°–75°N, σ = 0.1 to σ = 5 × 10−3, days 50–75) of the forced ensemble S. Positive values indicate positive zonal wind anomalies in C with a negative zonal wind anomaly in the box in S. Stippling indicates that the correlation is significantly larger than the correlation in a control ensemble at the 5% level.
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
6. Discussion and summary
We have used a large ensemble of idealized aquaplanet GCM simulations to determine the circulation response to short-term polar surface heating or cooling. We extracted robust and statistically significant aspects of the circulation response from the large background of internal variability using a large ensemble of simulations. The idealized modeling approach circumvents the complexity of the ocean–atmosphere boundary layer, seasonal changes, and stationary waves (Cohen et al. 2014; Overland et al. 2016). Yet many aspects of the circulation response we find in the idealized setting resemble circulation responses found in more complex models (Sun et al. 2015; Kim et al. 2014; Nakamura et al. 2016; Wu and Smith 2016; Zhang et al. 2017).
We chose an idealized GCM setup that produces a tropospheric mean state resembling that observed, while the stratospheric mean state lacks a separation of the zonal-mean jet. The key dynamical mechanisms we identified in this setup can be broken down into a fast thermal-wind response and a slower eddy-mediated response (Fig. 12). Polar surface heating leads to a warmer surface (red stripes on the surface in Fig. 12) and reduces meridional temperature gradients. By thermal-wind balance, the surface heating leads to weakened zonal winds aloft (Fig. 12, top panel). The tropospheric temperature and thermal-wind responses both decay on a radiative time scale (about 30 days; see appendix A) after the forcing ends. But while the tropospheric response persists, the weakened meridional temperature gradients lead to weakened upward wave activity (EP) fluxes in midlatitudes, which eventually lead to a weakened pattern of horizontal wave activity fluxes (momentum fluxes) in the midlatitude stratosphere (Fig. 12, bottom panel). The weakened pattern of eddy momentum flux convergence and divergence is balanced by a weakened Eulerian mean circulation in the stratosphere, with weakened rising motion (anomalous subsidence) in the mid- and high-latitude stratosphere (Fig. 12 lower panel, green arrows). This implies anomalous dynamic warming of the polar stratosphere, which, in turn, modulates the stratospheric zonal flow by thermal-wind balance. The long relaxation times in the stratosphere lead to relatively long (up to 90 days in our GCM) time scales of the stratospheric response, which can in some cases reconnect to the troposphere, when the background conditions set by internal variability are favorable (Fig. 10).

Schematic of (top) thermal-wind and (bottom) eddy-driven wind adjustments in response to short-term surface warming. Latitude increases from right (0°N) to left (90°N). Time increases from the front (day −30) to the back right (day 100). The tropopause is indicated by the gray layer and corner lines. Gray arrows represent climatological zonal-mean wind strength and position (at day 0 and day 40 for the top panel), with the spatial structure sketched at day 100 in red. Blue arrows illustrate the zonal wind weakening due to warming, either at the surface (red stripe) or in the lower stratosphere (only in the bottom panel). The surface heat flux forcing occurred within the green box at the surface, and the red shading illustrates the associated surface warming. The bottom panel is the same as the top panel, but shows the eddy-driven response. The wave activity flux into the stratosphere weakens (downward-pointing orange arrows), leading to weaker horizontal wave activity fluxes in the stratosphere (black arrows). Colored tubes illustrate the horizontal wave activity diverges (red) in the subtropics and converge (blue) in midlatitudes. Green arrows indicate the anomalous stratospheric Eulerian zonal-mean flow, which warms the polar stratosphere dynamically through anomalous subsidence (red area in the stratosphere) and results in weaker zonal wind in the stratosphere (blue arrows).
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1

Schematic of (top) thermal-wind and (bottom) eddy-driven wind adjustments in response to short-term surface warming. Latitude increases from right (0°N) to left (90°N). Time increases from the front (day −30) to the back right (day 100). The tropopause is indicated by the gray layer and corner lines. Gray arrows represent climatological zonal-mean wind strength and position (at day 0 and day 40 for the top panel), with the spatial structure sketched at day 100 in red. Blue arrows illustrate the zonal wind weakening due to warming, either at the surface (red stripe) or in the lower stratosphere (only in the bottom panel). The surface heat flux forcing occurred within the green box at the surface, and the red shading illustrates the associated surface warming. The bottom panel is the same as the top panel, but shows the eddy-driven response. The wave activity flux into the stratosphere weakens (downward-pointing orange arrows), leading to weaker horizontal wave activity fluxes in the stratosphere (black arrows). Colored tubes illustrate the horizontal wave activity diverges (red) in the subtropics and converge (blue) in midlatitudes. Green arrows indicate the anomalous stratospheric Eulerian zonal-mean flow, which warms the polar stratosphere dynamically through anomalous subsidence (red area in the stratosphere) and results in weaker zonal wind in the stratosphere (blue arrows).
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
Schematic of (top) thermal-wind and (bottom) eddy-driven wind adjustments in response to short-term surface warming. Latitude increases from right (0°N) to left (90°N). Time increases from the front (day −30) to the back right (day 100). The tropopause is indicated by the gray layer and corner lines. Gray arrows represent climatological zonal-mean wind strength and position (at day 0 and day 40 for the top panel), with the spatial structure sketched at day 100 in red. Blue arrows illustrate the zonal wind weakening due to warming, either at the surface (red stripe) or in the lower stratosphere (only in the bottom panel). The surface heat flux forcing occurred within the green box at the surface, and the red shading illustrates the associated surface warming. The bottom panel is the same as the top panel, but shows the eddy-driven response. The wave activity flux into the stratosphere weakens (downward-pointing orange arrows), leading to weaker horizontal wave activity fluxes in the stratosphere (black arrows). Colored tubes illustrate the horizontal wave activity diverges (red) in the subtropics and converge (blue) in midlatitudes. Green arrows indicate the anomalous stratospheric Eulerian zonal-mean flow, which warms the polar stratosphere dynamically through anomalous subsidence (red area in the stratosphere) and results in weaker zonal wind in the stratosphere (blue arrows).
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
It should be noted that this mechanism differs from that of sudden stratospheric warmings (SSW; or vortex breakdown; Matsuno 1971; Polvani and Waugh 2004; Garfinkel et al. 2010; Hitchcock et al. 2013a). It consequently differs from the usual explanations of Arctic to midlatitude leakage (Cohen et al. 2014; Overland et al. 2016). Our experiment shows a relatively slow (20–30 days) and weak (about 5%–10%) stratospheric warming resulting from a weakening of poleward heat fluxes in midlatitudes and subsequent adjustment of the momentum budget (Hirota and Sato 1969; Kushner and Polvani 2004; Limpasuvan et al. 2005; de la Cámara et al. 2018). In contrast, SSWs are rapid, intense events arising from positive heat flux anomalies by planetary waves that directly perturb the polar stratosphere. Some studies show evidence that polar warming can lead to SSWs (Wu and Smith 2016; Kim et al. 2014; Jaiser et al. 2013), but this pathway is weak compared to the midlatitude pathway in our simulations (Fig. B2b in appendix B; see also Wu and Smith 2016; Seierstad and Bader 2009; Peings and Magnusdottir 2014; Feldstein and Lee 2014; Sun et al. 2015; Nakamura et al. 2015). A better understanding of the conditions under which one pathway is more active than the other could help resolve why models differ in the sign and strength of their stratospheric responses to Arctic sea ice loss (Screen et al. 2018). The outlined mechanism can be also detected in more comprehensive models by analyzing momentum-flux divergences, or residual velocities, rather solely heat fluxes reaching the stratosphere.
A few findings from our simulations are especially noteworthy:
Responses with long time scales can be generated through the long memory of the stratosphere, without memory in surface or near-surface conditions. The cumulative effect of eddies with relatively short time scales on the stratospheric circulation is crucial for generating the late response.
Forcings that are zonally asymmetric and energetically equivalent (in the zonal mean) produce similar responses as forcings that are zonally symmetric in our GCM simulations. That is, stationary waves are not essential for the long-term response. However, stationary waves may still play important roles in the response to high-latitude forcings in Earth’s atmosphere and more complex models (e.g., Smith et al. 2010, 2011; Screen 2017; Smith and Scott 2016).
A robust feature of recent coupled model experiments is a weakening of the poleward flank or shift of the midlatitude jet (Screen et al. 2018; Smith et al. 2017; McCusker et al. 2017; Blackport and Kushner 2016; Deser et al. 2015). We suggest this response is consistent with the thermal adjustment shown in section 4 (Fig. 12, top). Coupled model experiments that show this zonal mean thermal-wind adjustment might also show a common stratospheric response to the reduced surface baroclinicity, as demonstrated here and in idealized models with a more realistic topography (Wu and Smith 2016).
The idealized model responses scale with the amount of heat release and peak around when the surface perturbation terminates, when both the thermal and eddy-driven response are simultaneously significant (section 3 and appendix A). In addition to seasonal changes in the jet position and strength, this may help explain sensitivities of the circulation response to the month and duration of sea ice reduction (King et al. 2016; Blackport et al. 2019).
The stratospheric overturning leads to downward motion in high latitudes, which has to be compensated by southward flow in lower levels, as seen in this model (Fig. 7d). This might be related to cold-air outbreaks in the real atmosphere (Kolstad et al. 2010; Overland et al. 2016).
Acknowledgments
The idealized GCM simulations were performed on ETH Zurich’s Euler computing cluster. We thank Noel Keenlyside for facilitating this work in an early stage, and Tobias Bischoff, Farid Ait-Chaalal, Robert Wills, and Ori Adam for support with running and analyzing the GCM as well as discussions of the paper. This work was partially supported by the NFR KLIMAFORSK programme (Project 255027).
APPENDIX A
Change of Perturbation Length
Figure A1 shows the response in four additional ensemble sets in addition to those in Fig. 4. The symmetric response

(a) Polar surface temperature and (b) lower-stratosphere NAM (at σ = 0.1) for four experiments with circumpolar forcing. Black lines show the symmetric response
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1

(a) Polar surface temperature and (b) lower-stratosphere NAM (at σ = 0.1) for four experiments with circumpolar forcing. Black lines show the symmetric response
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
(a) Polar surface temperature and (b) lower-stratosphere NAM (at σ = 0.1) for four experiments with circumpolar forcing. Black lines show the symmetric response
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
APPENDIX B
TEM Equations in σ Coordinates on a Sphere
a. Residual velocities

Time slices of the response of (left) residual meridional velocity, (center) residual vertical velocity, and (right) mean zonal wind for (a)–(c) days −30 to 0, (d)–(f) days 1 to 30, and (g)–(i) and days 31 to 60. (left) Shading shows the response of the residual meridional velocity
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1

Time slices of the response of (left) residual meridional velocity, (center) residual vertical velocity, and (right) mean zonal wind for (a)–(c) days −30 to 0, (d)–(f) days 1 to 30, and (g)–(i) and days 31 to 60. (left) Shading shows the response of the residual meridional velocity
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
Time slices of the response of (left) residual meridional velocity, (center) residual vertical velocity, and (right) mean zonal wind for (a)–(c) days −30 to 0, (d)–(f) days 1 to 30, and (g)–(i) and days 31 to 60. (left) Shading shows the response of the residual meridional velocity
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
b. Wave activity fluxes

(a) Vertical wave activity flux divergence averaged over days 1–30. Gray contours show the climatology, with solid as positive and dashed as negative. The thicker gray line is the zero line, and the black line indicates the tropopause. (b) As in (a), but for horizontal wave activity flux divergence. (c) Horizontal and vertical wave activity flux divergence ∇ ⋅ F [see (B9); shading] and wave activity flux vectors [see (B8)], scaled according to Edmon et al. (1980). Stippling indicates significance of the response at the 5% level.
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1

(a) Vertical wave activity flux divergence averaged over days 1–30. Gray contours show the climatology, with solid as positive and dashed as negative. The thicker gray line is the zero line, and the black line indicates the tropopause. (b) As in (a), but for horizontal wave activity flux divergence. (c) Horizontal and vertical wave activity flux divergence ∇ ⋅ F [see (B9); shading] and wave activity flux vectors [see (B8)], scaled according to Edmon et al. (1980). Stippling indicates significance of the response at the 5% level.
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
(a) Vertical wave activity flux divergence averaged over days 1–30. Gray contours show the climatology, with solid as positive and dashed as negative. The thicker gray line is the zero line, and the black line indicates the tropopause. (b) As in (a), but for horizontal wave activity flux divergence. (c) Horizontal and vertical wave activity flux divergence ∇ ⋅ F [see (B9); shading] and wave activity flux vectors [see (B8)], scaled according to Edmon et al. (1980). Stippling indicates significance of the response at the 5% level.
Citation: Journal of the Atmospheric Sciences 77, 2; 10.1175/JAS-D-19-0133.1
c. Thermodynamic equation
REFERENCES
Ambaum, M. H. P., and B. J. Hoskins, 2002: The NAO troposphere–stratosphere connection. J. Climate, 15, 1969–1978, https://doi.org/10.1175/1520-0442(2002)015<1969:TNTSC>2.0.CO;2.
Ambaum, M. H. P., B. J. Hoskins, and D. B. Stephenson, 2001: Arctic Oscillation or North Atlantic Oscillation? J. Climate, 14, 3495–3507, https://doi.org/10.1175/1520-0442(2001)014<3495:AOONAO>2.0.CO;2.
Andrews, D. G., and M. E. McIntyre, 1976: Planetary waves in horizontal and vertical shear: The generalized Eliassen-Palm relation and the mean zonal acceleration. J. Atmos. Sci., 33, 2031–2048, https://doi.org/10.1175/1520-0469(1976)033<2031:PWIHAV>2.0.CO;2.
Andrews, D. G., and M. E. McIntyre, 1978: Generalized Eliassen-Palm and Charney-Drazin theorems for waves in axisymmetric mean flows in compressible atmospheres. J. Atmos. Sci., 35, 175–185, https://doi.org/10.1175/1520-0469(1978)035<0175:GEPACD>2.0.CO;2.
Baldwin, M. P., and T. J. Dunkerton, 2001: Stratospheric harbingers of anomalous weather regimes. Science, 294, 581–584, https://doi.org/10.1126/science.1063315.
Blackport, R., and P. J. Kushner, 2016: The transient and equilibrium climate response to rapid summertime sea ice loss in CCSM4. J. Climate, 29, 401–417, https://doi.org/10.1175/JCLI-D-15-0284.1.
Blackport, R., and P. J. Kushner, 2017: Isolating the atmospheric circulation response to Arctic sea ice loss in the coupled climate system. J. Climate, 30, 2163–2185, https://doi.org/10.1175/JCLI-D-16-0257.1.
Blackport, R., and J. A. Screen, 2019: Influence of Arctic sea ice loss in autumn compared to that in winter on the atmospheric circulation. Geophys. Res. Lett., 46, 2213–2221, https://doi.org/10.1029/2018GL081469.
Blackport, R., J. A. Screen, K. Wiel, and R. Bintanja, 2019: Minimal influence of reduced Arctic sea ice on coincident cold winters in mid-latitudes. Nat. Climate Change, 9, 697–704, https://doi.org/10.1038/s41558-019-0551-4.
Bordoni, S., and T. Schneider, 2008: Monsoons as eddy-mediated regime transitions of the tropical overturning circulation. Nat. Geosci., 1, 515–519, https://doi.org/10.1038/ngeo248.
Borodina, A., E. M. Fischer, and R. Knutti, 2017: Emergent constraints in climate projections: A case study of changes in high-latitude temperature variability. J. Climate, 30, 3655–3670, https://doi.org/10.1175/JCLI-D-16-0662.1.
Cai, D., M. Dameris, H. Garny, and T. Runde, 2012: Implications of all season Arctic sea-ice anomalies on the stratosphere. Atmos. Chem. Phys., 12, 11 819–11 831, https://doi.org/10.5194/acp-12-11819-2012.
Cohen, J., and J. Jones, 2011: A new index for more accurate winter predictions. Geophys. Res. Lett., 38, L21701, https://doi.org/10.1029/2011GL049626.
Cohen, J., and Coauthors, 2014: Recent Arctic amplification and extreme mid-latitude weather. Nat. Geosci., 7, 627–637, https://doi.org/10.1038/ngeo2234.
Comiso, J. C., C. L. Parkinson, R. Gersten, and L. Stock, 2008: Accelerated decline in the Arctic sea ice cover. Geophys. Res. Lett., 35, L01703, https://doi.org/10.1029/2007GL031972.
Dee, D. P., and Coauthors, 2011: The ERA-Interim reanalysis: Configuration and performance of the data assimilation system. Quart. J. Roy. Meteor. Soc., 137, 553–597, https://doi.org/10.1002/qj.828.
de la Cámara, A., M. Abalos, and P. Hitchcock, 2018: Changes in stratospheric transport and mixing during sudden stratospheric warmings. J. Geophys. Res. Atmos., 123, 3356–3373, https://doi.org/10.1002/2017JD028007.
Deser, C., R. Tomas, M. Alexander, and D. Lawrence, 2010: The seasonal atmospheric response to projected Arctic sea ice loss in the late twenty-first century. J. Climate, 23, 333–351, https://doi.org/10.1175/2009JCLI3053.1.
Deser, C., R. Tomas, and L. Sun, 2015: The role of ocean–atmosphere coupling in the zonal-mean atmospheric response to Arctic sea ice loss. J. Climate, 28, 2168–2186, https://doi.org/10.1175/JCLI-D-14-00325.1.
Edmon, H. J., B. J. Hoskins, and M. E. McIntyre, 1980: Eliassen-Palm cross sections for the troposphere. J. Atmos. Sci., 37, 2600–2616, https://doi.org/10.1175/1520-0469(1980)037<2600:EPCSFT>2.0.CO;2.
Fang, Z., and J. M. Wallace, 1994: Arctic sea ice variability on a timescale of weeks and its relation to atmospheric forcing. J. Climate, 7, 1897–1914, https://doi.org/10.1175/1520-0442(1994)007<1897:ASIVOA>2.0.CO;2.
Feldstein, S. B., and S. Lee, 2014: Intraseasonal and interdecadal jet shifts in the Northern Hemisphere: The role of warm pool tropical convection and sea ice. J. Climate, 27, 6497–6518, https://doi.org/10.1175/JCLI-D-14-00057.1.
Frierson, D. M. W., I. M. Held, and P. Zurita-Gotor, 2006: A gray-radiation aquaplanet moist GCM. Part I: Static stability and eddy scale. J. Atmos. Sci., 63, 2548–2566, https://doi.org/10.1175/JAS3753.1.
Gao, Y., and Coauthors, 2015: Arctic sea ice and Eurasian climate: A review. Adv. Atmos. Sci., 32, 92–114, https://doi.org/10.1007/s00376-014-0009-6.
García-Serrano, J., C. Frankignoul, G. Gastineau, and A. de la Cámara, 2015: On the predictability of the winter Euro-Atlantic climate: Lagged influence of autumn Arctic sea ice. J. Climate, 28, 5195–5216, https://doi.org/10.1175/JCLI-D-14-00472.1.
Garfinkel, C. I., D. L. Hartmann, and F. Sassi, 2010: Tropospheric precursors of anomalous Northern Hemisphere stratospheric polar vortices. J. Climate, 23, 3282–3299, https://doi.org/10.1175/2010JCLI3010.1.
Held, I. M., and T. Schneider, 1999: The surface branch of the zonally averaged mass transport circulation in the troposphere. J. Atmos. Sci., 56, 1688–1697, https://doi.org/10.1175/1520-0469(1999)056<1688:TSBOTZ>2.0.CO;2.
Hirota, I., and Y. Sato, 1969: Periodic variation of the winter stratospheric circulation and intermittent vertical propagation. J. Meteor. Soc. Japan, 47, 390–402, https://doi.org/10.2151/jmsj1965.47.5_390.
Hitchcock, P., T. G. Shepherd, and G. L. Manney, 2013a: Statistical characterization of Arctic polar-night jet oscillation events. J. Climate, 26, 2096–2116, https://doi.org/10.1175/JCLI-D-12-00202.1.
Hitchcock, P., T. G. Shepherd, M. Taguchi, S. Yoden, and S. Noguchi, 2013b: Lower-stratospheric radiative damping and polar-night jet oscillation events. J. Atmos. Sci., 70, 1391–1408, https://doi.org/10.1175/JAS-D-12-0193.1.
Holmes, C. R., T. Woollings, E. Hawkins, and H. de Vries, 2016: Robust future changes in temperature variability under greenhouse gas forcing and the relationship with thermal advection. J. Climate, 29, 2221–2236, https://doi.org/10.1175/JCLI-D-14-00735.1.
Hoshi, K., and Coauthors, 2017: Poleward eddy heat flux anomalies associated with recent Arctic sea ice loss. Geophys. Res. Lett., 44, 446–454, https://doi.org/10.1002/2016GL071893.
Jaiser, R., K. Dethloff, and D. Handorf, 2013: Stratospheric response to Arctic sea ice retreat and associated planetary wave propagation changes. Tellus, 65, 19375, https://doi.org/10.3402/TELLUSA.V65I0.19375.
Kim, B.-M., S.-W. Son, S.-K. Min, J.-H. Jeong, S.-J. Kim, X. Zhang, T. Shim, and J.-H. Yoon, 2014: Weakening of the stratospheric polar vortex by Arctic sea-ice loss. Nat. Commun., 5, 4646, https://doi.org/10.1038/NCOMMS5646.
King, M. P., M. Hell, and N. Keenlyside, 2016: Investigation of the atmospheric mechanisms related to the autumn sea ice and winter circulation link in the Northern Hemisphere. Climate Dyn., 46, 1185–1195, https://doi.org/10.1007/s00382-015-2639-5.
Kolstad, E. W., T. Breiteig, and A. A. Scaife, 2010: The association between stratospheric weak polar vortex events and cold air outbreaks in the Northern Hemisphere. Quart. J. Roy. Meteor. Soc., 136, 886–893, https://doi.org/10.1002/QJ.620.
Kushner, P. J., and L. M. Polvani, 2004: Stratosphere–troposphere coupling in a relatively simple AGCM: The role of eddies. J. Climate, 17, 629–639, https://doi.org/10.1175/1520-0442(2004)017<0629:SCIARS>2.0.CO;2.
Lee, S., T. Gong, S. B. Feldstein, J. A. Screen, and I. Simmonds, 2017: Revisiting the cause of the 1989–2009 Arctic surface warming using the surface energy budget: Downward infrared radiation dominates the surface fluxes: Attribution of Arctic surface warming. Geophys. Res. Lett., 44, 10 654–10 661, https://doi.org/10.1002/2017GL075375.
Levine, X. J., and T. Schneider, 2011: Response of the Hadley circulation to climate change in an aquaplanet GCM coupled to a simple representation of ocean heat transport. J. Atmos. Sci., 68, 769–783, https://doi.org/10.1175/2010JAS3553.1.
Limpasuvan, V., D. L. Hartmann, D. W. J. Thompson, K. Jeev, and Y. L. Yung, 2005: Stratosphere-troposphere evolution during polar vortex intensification. J. Geophys. Res., 110, D24101, https://doi.org/10.1029/2005JD006302.
Liu, J., J. A. Curry, H. Wang, M. Song, and R. M. Horton, 2012: Impact of declining Arctic sea ice on winter snowfall. Proc. Natl. Acad. Sci. USA, 109, 4074–4079, https://doi.org/10.1073/pnas.1114910109.
Matsuno, T., 1971: A dynamical model of the stratospheric sudden warming. J. Atmos. Sci., 28, 1479–1494, https://doi.org/10.1175/1520-0469(1971)028<1479:ADMOTS>2.0.CO;2.
McCusker, K. E., P. J. Kushner, J. C. Fyfe, M. Sigmond, V. V. Kharin, and C. M. Bitz, 2017: Remarkable separability of circulation response to Arctic sea ice loss and greenhouse gas forcing. Geophys. Res. Lett., 44, 7955–7964, https://doi.org/10.1002/2017GL074327.
Nakamura, T., K. Yamazaki, K. Iwamoto, M. Honda, Y. Miyoshi, Y. Ogawa, and J. Ukita, 2015: A negative phase shift of the winter AO/NAO due to the recent Arctic sea-ice reduction in late autumn. J. Geophys. Res. Atmos., 120, 3209–3227, https://doi.org/10.1002/2014JD022848.
Nakamura, T., K. Yamazaki, K. Iwamoto, M. Honda, Y. Miyoshi, Y. Ogawa, Y. Tomikawa, and J. Ukita, 2016: The stratospheric pathway for Arctic impacts on midlatitude climate. Geophys. Res. Lett., 43, 3494–3501, https://doi.org/10.1002/2016GL068330.
O’Gorman, P. A., and T. Schneider, 2008: The hydrological cycle over a wide range of climates simulated with an idealized GCM. J. Climate, 21, 3815–3832, https://doi.org/10.1175/2007JCLI2065.1.
Overland, J. E., and Coauthors, 2016: Nonlinear response of mid-latitude weather to the changing Arctic. Nat. Climate Change, 6, 992–999, https://doi.org/10.1038/nclimate3121.
Peings, Y., and G. Magnusdottir, 2014: Forcing of the wintertime atmospheric circulation by the multidecadal fluctuations of the North Atlantic Ocean. Environ. Res. Lett., 9, 034018, https://doi.org/10.1088/1748-9326/9/3/034018.
Petoukhov, V., S. Rahmstorf, S. Petri, and H. J. Schellnhuber, 2013: Quasiresonant amplification of planetary waves and recent Northern Hemisphere weather extremes. Proc. Natl. Acad. Sci. USA, 110, 5336–5341, https://doi.org/10.1073/pnas.1222000110.
Polvani, L. M., and D. W. Waugh, 2004: Upward wave activity flux as a precursor to extreme stratospheric events and subsequent anomalous surface weather regimes. J. Climate, 17, 3548–3554, https://doi.org/10.1175/1520-0442(2004)017<3548:UWAFAA>2.0.CO;2.
Scaife, A. A., J. R. Knight, G. K. Vallis, and C. K. Folland, 2005: A stratospheric influence on the winter NAO and North Atlantic surface climate. Geophys. Res. Lett., 32, L18715, https://doi.org/10.1029/2005GL023226.
Schneider, T., and I. M. Held, 2001: Discriminants of twentieth-century changes in Earth surface temperatures. J. Climate, 14, 249–254, https://doi.org/10.1175/1520-0442(2001)014<0249:LDOTCC>2.0.CO;2.
Schneider, T., and P. A. O’Gorman, 2008: Moist convection and the thermal stratification of the extratropical troposphere. J. Atmos. Sci., 65, 3571–3583, https://doi.org/10.1175/2008JAS2652.1.
Scinocca, J. F., M. C. Reader, D. A. Plummer, M. Sigmond, P. J. Kushner, T. G. Shepherd, and A. R. Ravishankara, 2009: Impact of sudden Arctic sea-ice loss on stratospheric polar ozone recovery. Geophys. Res. Lett., 36, L24701, https://doi.org/10.1029/2009GL041239.
Screen, J. A., 2017: Simulated atmospheric response to regional and pan-Arctic sea-ice loss. J. Climate, 30, 3945–3962, https://doi.org/10.1175/JCLI-D-16-0197.1.
Screen, J. A., and I. Simmonds, 2010: The central role of diminishing sea ice in recent Arctic temperature amplification. Nature, 464, 1334–1337, https://doi.org/10.1038/nature09051.
Screen, J. A., and I. Simmonds, 2013: Exploring links between Arctic amplification and mid-latitude weather. Geophys. Res. Lett., 40, 959–964, https://doi.org/10.1002/grl.50174.
Screen, J. A., C. Deser, I. Simmonds, and R. Tomas, 2013: Atmospheric impacts of Arctic sea-ice loss, 1979–2009: Separating forced change from atmospheric internal variability. Climate Dyn., 43, 333–344, https://doi.org/10.1007/s00382-013-1830-9.
Screen, J. A., and Coauthors, 2018: Consistency and discrepancy in the atmospheric response to Arctic sea-ice loss across climate models. Nat. Geosci., 11, 155–163, https://doi.org/10.1038/s41561-018-0059-y.
Seierstad, I., and J. Bader, 2009: Impact of a projected future Arctic sea ice reduction on extratropical storminess and the NAO. Climate Dyn., 33, 937–943, https://doi.org/10.1007/s00382-008-0463-x.