1. Introduction
The planetary-scale stationary wave plays an important role for regulating weather and climate in the Northern Hemisphere extratropics during the winter. Constructive interference with transient waves accentuates the effect of the stationary wave. In particular, it has been shown that stationary–transient wave constructive interference plays an important role in the extratropical response to tropical forcing (Fletcher and Kushner 2011; Smith et al. 2011; Garfinkel et al. 2012; Goss et al. 2016, hereafter GFL). It has been proposed that constructive interference is also likely to play a crucial role in maintaining equable climate conditions (i.e., where the equator-to-pole temperature gradient is relatively weak) even during winter (Lee 2012, 2014). The basis of this equable climate theory is that the stationary wave is a major contributor to the poleward energy flux in mid- to high latitudes during boreal winter (Peixoto and Oort 1992); zonal asymmetry in the climatological tropical latent heating, while independent of the equator-to-pole temperature gradient, is important for forcing the stationary wave (Hoskins and Karoly 1981; Ting and Held 1990; Ting 1996; Held et al. 2002; Chang 2009); and there is a large untapped reservoir of zonal available potential energy (Lorenz 1955), which in principle can be unleashed by disturbances from outside of the extratropical baroclinic zone, such as the tropics (Lee 2014).
GFL hypothesized that constructive interference would be preceded by tropical heating anomalies that reinforce the climatological zonal asymmetry, and that it would be followed by anomalous warming in high latitudes. To test this hypothesis, GFL constructed an index, referred to as the stationary wave index (SWI), which is positive (negative) for constructive (destructive) interference that measures the degree to which the instantaneous daily 300-hPa streamfunction field matches with its winter climatology. It was found that the SWI tends to be positive 7–10 days after an enhanced warm pool convection, which enhances the zonal asymmetry of the climatological heating field in the tropics. Positive SWI days are typically followed by Arctic warming 8–10 days later. Examining the stationary wave responses to individual forcing terms, presented by Held et al. (2002), we find that the responses to diabatic heating fields—tropical and extratropical heating individually—exhibit circulation patterns that promote warmth in the Arctic. Specifically, the wave solutions are composed of a southerly flow over the Bering Strait and Greenland and Norwegian Seas where most of the warm, moist air intrusions occur (Woods et al. 2013; Baggett et al. 2016; Woods and Caballero 2016). The southerly flow is especially prominent in the response to the extratropical heating.
Extratropical heating and tropical heating are unlikely to be independent of each other, however. The results of Baggett et al. (2016) raise the possibility that tropical heating leads to Arctic warming not just by its direct excitation of poleward and upward propagating planetary-scale waves, but also by amplification of these planetary-scale waves through tropically induced extratropical diabatic heating anomaly. Therefore, we hypothesize that during SWI days, there is a latent heat–circulation relay process that takes place as follows: tropical latent heating anomalies excite extratropical circulation anomalies; the extratropical circulation then drives extratropical latent heating anomalies by transporting moisture; the resulting extratropical latent heating, in turn, drives additional extratropical circulation. This hypothesized sequence is somewhat akin to how a hybrid vehicle operates in the sense that the tropical latent heating is analogous to gasoline fuel put into the vehicle and the extratropical latent heating is analogous to electric power generated by running the vehicle. This hypothesis of a latent heat–circulation relay is also supported by the studies of Willison et al. (2013) and Papritz and Spengler (2015) who showed that enhanced synoptic-scale wave activity is often associated with additional latent heat release.
In this study, we test the relay hypothesis in the context of the SWI by performing observational analyses and initial-value model calculations. Specifically, we address the following two questions: 1) During SWI events, what circulation and temperature anomalies are induced by the individual heating anomalies over the tropics and extratropics? 2) During SWI events, is there any evidence that tropical heating anomalies lead the extratropical heating anomalies, and then the stationary–transient wave interference? In section 2, we describe the data, analysis methods, and the setup of the idealized model experiments. The results of our observational analyses are presented in section 3, while the model results are discussed in section 4. Last, section 5 provides a discussion and summary.
2. Data and methods
a. Data
Daily data of zonal wind, meridional wind, temperature, specific humidity, vertically integrated eastward/northward moisture flux, and surface pressure have been acquired from European Centre for Medium-Range Weather Forecasts (ECMWF) interim reanalysis (ERA-Interim; Dee et al. 2011). All variables are on a grid with 2.5° × 2.5° horizontal resolution and with 23 pressure levels, except for the vertically integrated moisture flux and surface pressure. The time period examined here is December–February (DJF) for years 1979–2014. For our study, latent heating generated by convection and large-scale condensation are important variables. Because these individual latent heating variables are unavailable in the ERA-Interim data, we use diabatic heating data provided by the Japanese 55-year Reanalysis (JRA-55; Kobayashi et al. 2015) with a 2.5° × 2.5° horizontal resolution and 37 pressure levels.
There are differences in the climatological diabatic heating among different reanalysis datasets (Ling and Zhang 2013; Wright and Fueglistaler 2013). However, Clark and Feldstein (2019, manuscript submitted to J. Atmos. Sci.) show excellent agreement in nonradiative diabatic heating between ERA-Interim and JRA-55 in their North Atlantic Oscillation (NAO) composites. This comparison is possible because the ERA-Interim provides total diabatic heating and radiative heating; nonradiative diabatic heating was computed by subtracting radiative heating from the total diabatic heating. In the JRA-55 data, the nonradiative diabatic heating was computed by summing the convective heating Qcnv large-scale condensational heating Qlrg, and subgrid-scale vertical diffusion. Also, Zhang et al. (2017) calculated the composite values of nonradiative diabatic heating against outgoing longwave radiation over the western Pacific at 400 hPa from ERA-Interim, JRA-55, and the Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2), and found that the agreement between ERA-Interim and JRA-55 stood out. As a further test to ascertain that our results are insensitive to the choice of the dataset, we computed the stationary wave index (defined in section 2b) using JRA-55 and compared it with the same index computed using the ERA-Interim data. The correlation coefficient between the two time series turned out to be 0.964. In addition, the 300-hPa streamfunction anomaly composites (not shown) from the JRA-55 dataset are essentially indistinguishable from those from the ERA-Interim (Fig. 1). These results lend confidence that JRA-55 Qcnv + Qlrg can be used to investigate the relationship between the latent heating and the circulation.

(a)–(g) Total 300-hPa streamfunction (contours, interval of 1.5 × 107 m2 s−1) and anomalies (shading), and (h)–(n) vertically averaged latent heating anomalies during the SWI+ days. Dotted areas indicate statistical significance at the 10% level. Statistical significance is evaluated by employing a Monte Carlo simulation with 1000 random samples.
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1

(a)–(g) Total 300-hPa streamfunction (contours, interval of 1.5 × 107 m2 s−1) and anomalies (shading), and (h)–(n) vertically averaged latent heating anomalies during the SWI+ days. Dotted areas indicate statistical significance at the 10% level. Statistical significance is evaluated by employing a Monte Carlo simulation with 1000 random samples.
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1
(a)–(g) Total 300-hPa streamfunction (contours, interval of 1.5 × 107 m2 s−1) and anomalies (shading), and (h)–(n) vertically averaged latent heating anomalies during the SWI+ days. Dotted areas indicate statistical significance at the 10% level. Statistical significance is evaluated by employing a Monte Carlo simulation with 1000 random samples.
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1
b. Stationary wave index and diabatic heating index
For the projection domains (λi, θj), we define the tropics as 30°S–30°N, 0°–360°E, and the North Pacific (North Atlantic) domain as latitudes between 30° and 70°N (30° and 80°N) and longitudes between 150° and 260°E (280° and 360°E). The North Pacific (North Atlantic) domain is indicated by the yellow (green) box in Fig. 1j. As was shown by GFL, tropical diabatic heating associated with the SWI tends to peak during lag days −10 to 0 relative to the peak of SWI, where lag day 0 corresponds to the days when SWI exceeds one standard deviation. Therefore, for the tropical domain,

Lag composites of daily heating indices (black for tropical, red for North Pacific, and blue for North Atlantic) during (left) SWI+ and (right) SWI− days. Row are composites during (a),(b) entire SWI days, (c),(f) SWI and tropical heating index >1σ days, (d),(g) SWI and North Pacific heating index >1σ days, and (e),(h) SWI and North Atlantic heating index >1σ days. A 3-day running mean was applied to the heating fields shown in (c)–(h). A Monte Carlo simulation with 1000 random samples is performed for the statistical significance test, and the thick lines indicate statistical significance at the 5% level.
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1

Lag composites of daily heating indices (black for tropical, red for North Pacific, and blue for North Atlantic) during (left) SWI+ and (right) SWI− days. Row are composites during (a),(b) entire SWI days, (c),(f) SWI and tropical heating index >1σ days, (d),(g) SWI and North Pacific heating index >1σ days, and (e),(h) SWI and North Atlantic heating index >1σ days. A 3-day running mean was applied to the heating fields shown in (c)–(h). A Monte Carlo simulation with 1000 random samples is performed for the statistical significance test, and the thick lines indicate statistical significance at the 5% level.
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1
Lag composites of daily heating indices (black for tropical, red for North Pacific, and blue for North Atlantic) during (left) SWI+ and (right) SWI− days. Row are composites during (a),(b) entire SWI days, (c),(f) SWI and tropical heating index >1σ days, (d),(g) SWI and North Pacific heating index >1σ days, and (e),(h) SWI and North Atlantic heating index >1σ days. A 3-day running mean was applied to the heating fields shown in (c)–(h). A Monte Carlo simulation with 1000 random samples is performed for the statistical significance test, and the thick lines indicate statistical significance at the 5% level.
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1
c. Binning procedure–based heating indices
Last, to answer the first question posed in section 1, we divided the SWI days into multiple bins by sorting them based on the magnitude of the three heating indices. Because tropical heating anomalies occur first (Fig. 2), the tropical heating index is used as the first criterion for creating the bins. For example, all SWI+ days are divided into three bins, with the first bin corresponding to those SWI+ days with the top one-third of the
The North Pacific (Atlantic) domain is then used as the second (third) criterion. For these extratropical heating criteria, only two bins are used to retain a sufficient number of SWI days in each of the bins; within the
Summarizing the binning procedure, SWI+ days are divided into multiple bins and, among these bins, we analyze
Acronym definitions after binning procedures and the number of DJF days in each bin out of 3249 all DJF days.


For the purpose of analyzing time sequence of events, we identify individual events based on the following procedures. For each bin, an SWI event is defined as a 15-day time interval that contains at least one SWI day that satisfies the binning criteria and the day of the maximum (or minimum for destructive interference) SWI value. The eighth day within the 15-day interval coincides with the maximum or minimum SWI value, and is defined as the SWI lag 0 days. With this definition, consecutive SWI events are separated from each other by at least 7 days. For the statistical significance of the composite values presented in Figs. 1–5, Monte Carlo simulations are performed. For composites based on N number of events, we generated 1000 composites, with each composite consisting of N randomly chosen events with each event having a time interval extending from lag −10 to lag +10 days (Fig. 1) and from lag −20 to lag +20 days (Figs. 2–5). A null distribution was then constructed from the 1000 random composites for each lag day, and the p value of the observed composite was computed based on the distribution. Again following GFL, to account for the fact that the consecutive days are not independent of each other, for Fig. 1, we divided N number of days by 3.38 and rounded to the nearest integer. In Fig. 1, the stippled region indicates statistical significance at the p < 0.10 level, while in Figs. 2–5, the p < 0.05 level is indicated.
d. Model experiment setup
To test the causal relationships that emerge from our diagnostic analyses, we utilize the spectral dynamic core from the NOAA Geophysical Fluid Dynamics Laboratory (GFDL). The same model setup as in Baggett et al. (2016) is used in this study: a horizontal resolution of triangular 42, a vertical resolution of 28 sigma levels, a damping time scale of 0.1 days at its smallest scale for a fourth-order horizontal diffusion, Newtonian cooling, and Rayleigh friction parameterized as in Held and Suarez (1994). The climatological DJF values of zonal wind, meridional wind, temperature, and surface pressure are used as the initial background state. For passive tracer experiments, we use the zonal-mean climatological specific humidity field. To ensure that the zonally varying DJF climatological state is a solution to the model equations, we add a forcing term that is obtained by integrating the model by one time step starting from the climatological state. This forcing term prevents the model from drifting from the initial state (the DJF climatological state) unless additional forcing (e.g., diabatic heating) is added. Additional details on the model setup and limitations can be found in Franzke et al. (2004). In the control experiment, the model is integrated without any additional forcing, whereas for perturbation experiments, the model is forced by time-dependent (lag days −10 to +7)
3. Observational analysis
a. Relationship between circulation and heating anomalies during SWI
We first examine the 300-hPa streamfunction evolution for SWI+ days (Figs. 1a–g). Because GFL provides a detailed description of the same fields, we present just a few key highlights pertinent to this study. It can be seen that the streamfunction anomalies (shading) in Fig. 1 grow and decay over a period of about two weeks. For the SWI+, by construction, the positive (negative) anomalies coincide with the climatological ridges (troughs). The opposite is the case for the SWI composites (see Fig. S1 in the online supplemental material). Poleward of 50°N, there are two ridges, one centered at the eastern end of Gulf of Alaska and the other centered over the British Isles. As such, associated with these ridges, there are southerly flows over the two ocean corridors to the Arctic Ocean, one over the Bering Sea and the other over the Norwegian Sea.
Figures 1h–n illustrate the composite heating anomalies during SWI+ days. When there is constructive interference, consistent with GFL, there is enhanced warm-pool convection and suppressed convection in central tropical Pacific, reinforcing the zonal asymmetry of the climatological heating field. In the extratropics, which was not examined by GFL, a positive anomaly is seen over the central North Pacific, while a negative anomaly is found over the eastern North Pacific and western North America. These heating anomalies appear at lag day −6 and dissipate during positive lag days (Figs. 1i–m). We also see heating anomalies over the northeastern Atlantic and cooling over subtropical Atlantic shortly after the North Pacific heating is excited. In the case of destructive interference, SWI−, streamfunction and heating anomalies with the opposite signs are found over the tropics and extratropics (see Fig. S1).
Synthesizing the results of the streamfunction and heating field composites, we highlight two findings: First, the SWI+ (SWI−) heating anomaly field reinforces (dampens) the zonally asymmetric component of climatological diabatic heating (e.g., Fig. 8 of Held et al. 2002). This relationship suggests that a strengthening of the zonally asymmetric component of the climatological diabatic heating leads to an amplification of the climatological stationary wave. Because diabatic heating is a major driver of the climatological stationary wave, this transient heating–SWI relationship is consistent with stationary wave theory (e.g., Held et al. 2002; Chang 2009). Second, as was briefly discussed in the previous section, the tropical heating anomalies tend to precede the extratropical heating anomalies by several days, suggesting that the SWI extratropical heating anomalies are, at least in part, driven by the circulation driven by the tropical heating anomalies. As a first step to explore this possibility, we next examine lead–lag relationships between tropical and extratropical heating.
b. Relationship between tropical and extratropical heating anomalies during SWI
Figure 2 shows the heating index composites associated with SWI+ (left column) and SWI− (right column). We first examine composites over all SWI+ and SWI− days, that is,
As a way to evaluate the extent of the association among these three heating anomalies, the aforementioned composite analysis was repeated, first with the condition that
c. Four types of heating anomalies during SWI and associated circulation pattern
Before presenting results from the model calculations, we first examine the structure of the heating and circulation anomalies for the following cases:

Composites of latent heating anomalies averaged from lag days −10 to +3 in (a)
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1

Composites of latent heating anomalies averaged from lag days −10 to +3 in (a)
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1
Composites of latent heating anomalies averaged from lag days −10 to +3 in (a)
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1
The zonal asymmetry in the heating anomalies suggests that the SWI events might be influenced by the Madden–Julian oscillation (MJO; Madden and Julian 1971), which is a 30–60-day tropical variability with eastward-propagating convection. By calculating the composites of the real-time multivariate MJO (RMM; Wheeler and Hendon 2004) indices for each bin, we found that in the
Figure 4 shows the 300-hPa streamfunction anomalies in each bin, averaged from lag days −1 to +1. We overlay vectors of vertically integrated moisture flux poleward of 30°N. In all four bins of SWI+, Ψ′ (shading) reinforces Ψ (black contours). However, there are some differences. Figures 4a–d reveal that intense moisture transport into the Arctic is present in the

Composites of 300-hPa streamfunction anomalies (shading) averaged from lag days −1 to +1 in the same subsets as in Fig. 3, and total 300-hPa streamfunction (contours, interval of 1.5 × 107 m2 s−1). Overlaid vectors represent vertically integrated moisture flux (kg m−1 s−1) north of 30°N. The vectors whose magnitude is less than 10 kg m−1 s−1 are omitted, and reference vector is 70 kg m−1 s−1. A Monte Carlo simulation with 1000 random samples is performed for the statistical significance test, and the dotted areas indicate statistical significance at the 5% level.
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1

Composites of 300-hPa streamfunction anomalies (shading) averaged from lag days −1 to +1 in the same subsets as in Fig. 3, and total 300-hPa streamfunction (contours, interval of 1.5 × 107 m2 s−1). Overlaid vectors represent vertically integrated moisture flux (kg m−1 s−1) north of 30°N. The vectors whose magnitude is less than 10 kg m−1 s−1 are omitted, and reference vector is 70 kg m−1 s−1. A Monte Carlo simulation with 1000 random samples is performed for the statistical significance test, and the dotted areas indicate statistical significance at the 5% level.
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1
Composites of 300-hPa streamfunction anomalies (shading) averaged from lag days −1 to +1 in the same subsets as in Fig. 3, and total 300-hPa streamfunction (contours, interval of 1.5 × 107 m2 s−1). Overlaid vectors represent vertically integrated moisture flux (kg m−1 s−1) north of 30°N. The vectors whose magnitude is less than 10 kg m−1 s−1 are omitted, and reference vector is 70 kg m−1 s−1. A Monte Carlo simulation with 1000 random samples is performed for the statistical significance test, and the dotted areas indicate statistical significance at the 5% level.
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1
Figure 5 presents temperature anomalies averaged vertically between 700 and 1000 hPa. During the 5 days following the SWI events (left column), every bin shows warming at least either over the Pacific or the Atlantic sector of Arctic Ocean. In later lag days, however, Arctic warming is observed only in the

Composites of lower-tropospheric (700–1000 hPa) temperature anomalies during SWI+ days averaged (a)–(d) from lag days 0 to +5 and (e)–(h) from lag days +10 to +20 in the same subsets for SWI+ days as in Figs. 3a–d. Gray shading denotes the region where surface pressure is below 700 hPa. A Monte Carlo simulation with 1000 random samples is performed for the statistical significance test, and the dotted areas indicate statistical significance at the 5% level.
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1

Composites of lower-tropospheric (700–1000 hPa) temperature anomalies during SWI+ days averaged (a)–(d) from lag days 0 to +5 and (e)–(h) from lag days +10 to +20 in the same subsets for SWI+ days as in Figs. 3a–d. Gray shading denotes the region where surface pressure is below 700 hPa. A Monte Carlo simulation with 1000 random samples is performed for the statistical significance test, and the dotted areas indicate statistical significance at the 5% level.
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1
Composites of lower-tropospheric (700–1000 hPa) temperature anomalies during SWI+ days averaged (a)–(d) from lag days 0 to +5 and (e)–(h) from lag days +10 to +20 in the same subsets for SWI+ days as in Figs. 3a–d. Gray shading denotes the region where surface pressure is below 700 hPa. A Monte Carlo simulation with 1000 random samples is performed for the statistical significance test, and the dotted areas indicate statistical significance at the 5% level.
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1
4. Model results
In this section, we use the model (see section 2d for details) to test the causal relationships suggested by the observational analysis presented in the previous section. Specifically, we first ask the extent to which the SWI+ heating anomalies can explain the associated circulation and temperature anomalies. Next, we ask if the tropical heating anomalies can lead to the North Pacific and North Atlantic heating anomalies. Because water vapor is the central ingredient of latent heating, to make headway toward addressing the second question, we use a passive tracer to represent moisture. This approach is supported by the earlier studies that show that large-scale advection can account for much of the water vapor distribution in the free troposphere (Pierrehumbert et al. 2007, and references therein; Baggett et al. 2016; Ming and Held 2018).
a. Circulation and temperature response
For brevity, we only present results for the SWI+ experiments because the model responses to the SWI− heating anomalies are by and large opposite to those of the SWI+ heating. We also confine our analysis to the

The 300-hPa streamfunction anomaly averaged from model days 10 to 12. The model is forced with the heating composite of the
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1

The 300-hPa streamfunction anomaly averaged from model days 10 to 12. The model is forced with the heating composite of the
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1
The 300-hPa streamfunction anomaly averaged from model days 10 to 12. The model is forced with the heating composite of the
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1
Figure 7 shows the vertically averaged (700–1000 hPa) temperature response averaged from model days 11 to 16, which corresponds to the lag days 0 to +5 in observation (Fig. 5a). In Fig. 7a, tropical heating causes warming over eastern Siberia and Greenland and cooling over northern Canada and the Norwegian Sea. The North Pacific heating leads to warming over a broad swath of area poleward of 50°N, ranging from northern North America to Scandinavia, and cooling over eastern Siberia and much of the contiguous United States (Fig. 7b). Figure 7c shows that the North Atlantic heating results in warming (cooling) over the Barents–Kara Seas and northern Europe (eastern Greenland and the Mediterranean Sea). Figure 7d shows the sum of the three model solutions. Compared with the observation (Fig. 5a), we see a reasonable agreement except over northeastern Canada and Greenland where the model solution shows warming while the observations show cooling. We found that horizontal temperature advection is the primary contributor to the temperature anomalies (not shown). Overall, the model temperature solutions indicate that tropical heating is important for warming the Pacific sector of the Arctic, while extratropical heating is more important for warming the Atlantic sector.

The lower-tropospheric (700–1000 hPa) temperature anomaly averaged from model days 11 to 16. The same heating composites as in Fig. 6 are used to force the model, but only over (a) the tropical, (b) the North Pacific, and (c) the North Atlantic domain. (d) Summation of (a)–(c). Gray shading denotes the region where surface pressure is below 700 hPa.
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1

The lower-tropospheric (700–1000 hPa) temperature anomaly averaged from model days 11 to 16. The same heating composites as in Fig. 6 are used to force the model, but only over (a) the tropical, (b) the North Pacific, and (c) the North Atlantic domain. (d) Summation of (a)–(c). Gray shading denotes the region where surface pressure is below 700 hPa.
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1
The lower-tropospheric (700–1000 hPa) temperature anomaly averaged from model days 11 to 16. The same heating composites as in Fig. 6 are used to force the model, but only over (a) the tropical, (b) the North Pacific, and (c) the North Atlantic domain. (d) Summation of (a)–(c). Gray shading denotes the region where surface pressure is below 700 hPa.
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1
b. Passive tracer response
We next address the question of whether tropical heating anomalies can lead to the North Pacific and North Atlantic heating anomalies. Figures 8a and 8b show passive tracer fields simulated by the model. The initial condition of the tracer is the three-dimensional climatological specific humidity during boreal winter. During model days 10–12, it can be seen that the circulation driven by tropical heating transports anomalous tracer over the central North Pacific and western North America, as well as over the east coast of North America and the northern subtropical North Atlantic (Fig. 8a). Figure 8b shows that North Pacific heating also contributes to the positive tracer anomaly over northwestern North America and the northeastern North Atlantic.

(a),(b) The vertically integrated tracer anomaly averaged over model days (a) 10–12 and (b) 11–13. The same heating composites as in Fig. 6 are used to force the model, but only over (left) the tropical and (right) the North Pacific domain. (c),(d) The vertically integrated condensational heating computed from the tracer anomalies in (a) and (b), respectively (see section 4b and appendix for details).
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1

(a),(b) The vertically integrated tracer anomaly averaged over model days (a) 10–12 and (b) 11–13. The same heating composites as in Fig. 6 are used to force the model, but only over (left) the tropical and (right) the North Pacific domain. (c),(d) The vertically integrated condensational heating computed from the tracer anomalies in (a) and (b), respectively (see section 4b and appendix for details).
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1
(a),(b) The vertically integrated tracer anomaly averaged over model days (a) 10–12 and (b) 11–13. The same heating composites as in Fig. 6 are used to force the model, but only over (left) the tropical and (right) the North Pacific domain. (c),(d) The vertically integrated condensational heating computed from the tracer anomalies in (a) and (b), respectively (see section 4b and appendix for details).
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1
Because this passive tracer transport represents moisture transport, the model condensational heating can be estimated (see appendix for the method). The results (Figs. 8c and 8d) show horizontal structures that resemble the composite heating field (Fig. 3a), supporting that the lead–lag relationships in Fig. 2 are indeed causal. The model result also suggests a positive feedback process where the North Atlantic heating drives a circulation that can reinforce the heating that drove the circulation in the first place. This positive feedback is reminiscent of the mechanism of Hoskins and Valdes (1990), but it differs in the sense that the feedback suggested here is through the interplay between latent heating and circulation while the mechanism of Hoskins and Valdes (1990) is self-maintaining through latent heating, circulation, and baroclinity. When tropical and extratropical heating are imposed together, anomalous tracer transport poleward of 60°N is found mostly over the two ocean corridors where anomalous moisture transport into the Arctic Ocean occur (not shown). This result shows that when forced by latent heating, moisture transport, via an increase in downward infrared radiation, is likely to reinforce temperature changes caused by dry dynamics (Yoo et al. 2012).
To quantify the tropics–North Pacific–North Atlantic linkage that involves heating and “moisture” transport, we compute area-weighted averages of the net tracer anomaly over the North Pacific and North Atlantic domains forced by tropical, North Pacific, and North Atlantic heating. Figure 9a shows that, for the Pacific domain, tropical heating is the main driver of tracer transport, whereas the net effect of extratropical heating is nearly zero. For the Atlantic domain illustrated in Fig. 9b, the tracer anomaly forced by North Pacific heating is predominant. To further evaluate the resemblance between the model tracer anomaly and the reanalysis heating anomaly, we compute the projection of the condensational heating estimated from tracer anomalies, shown in Figs. 8c and 8d, onto the composites of heating anomalies

(a),(b) Time series of the area-weighted average of the vertically integrated tracer anomaly. (c),(d) Time series of projection of the model condensational heating anomaly computed from the model tracer onto the observed latent heating anomaly field in the
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1

(a),(b) Time series of the area-weighted average of the vertically integrated tracer anomaly. (c),(d) Time series of projection of the model condensational heating anomaly computed from the model tracer onto the observed latent heating anomaly field in the
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1
(a),(b) Time series of the area-weighted average of the vertically integrated tracer anomaly. (c),(d) Time series of projection of the model condensational heating anomaly computed from the model tracer onto the observed latent heating anomaly field in the
Citation: Journal of the Atmospheric Sciences 76, 9; 10.1175/JAS-D-18-0371.1
5. Summary and conclusions
In this study, we addressed the following two questions: 1) During SWI events, what circulation and temperature anomalies are induced by the individual heating anomalies over the tropics and extratropics? 2) During SWI events, is there any evidence that extratropical heating is excited by waves forced by tropical heating?
To address the first question, we binned the SWI+ days by ranking the magnitude of tropical and extratropical latent heating. Our analyses reveal that there are different flavors of stationary wave interferences. Simultaneous enhancements in zonal asymmetries in both tropical and extratropical heating
To address the second question, we first examined the lead–lag relationships among the tropical and the two extratropical heating indices. The result indicates that while the latent heating anomalies in these three domains can occur by themselves, they tend to occur together within 7–10 days of each other, with the tropical heating anomaly leading the North Pacific heating anomaly, which in turn is followed by the North Atlantic anomaly. This finding suggests that not only the answer to this second question be positive, but also that the circulation driven by the North Pacific heating enhances latent heating over the North Atlantic domain. We tested this heating–circulation relay hypothesis with the model by computing condensational heating from the model’s passive tracer, which represents specific humidity. This initial-value passive tracer experiment indeed supports the heating–circulation relay hypothesis that emerged from the observational analysis: tropical latent heating → circulation anomalies → latent heating in the North Pacific → circulation anomalies → latent heating in the North Atlantic.
The relay picture is a reminder that diabatic heating is not only an important driver of the atmospheric circulation (e.g., Sutcliffe 1951; Hoskins and Valdes 1990), but it also reveals that the extratropical diabatic heating is dependent on tropical heating. The implication is that while the atmospheric response to the individual components of the heating is linear, they are not independent of each other. Therefore, the impact of the tropical heating should not be readily dismissed even for a circulation feature attributable to extratropical heating.
Acknowledgments
We acknowledge discussion with Steven Feldstein and Joseph Clark throughout this study and Steven Feldstein for his comments on this manuscript. We also acknowledge comments by anonymous reviewers. This research is supported by National Science Foundation Grants AGS-1455577, AGS-1822015 and OPP-1723832.
APPENDIX
Estimation of Condensational Heating from Model Passive Tracer
REFERENCES
Baggett, C., S. Lee, and S. B. Feldstein, 2016: An investigation of the presence of atmospheric rivers over the North Pacific during planetary-scale wave life cycles and their role in Arctic warming. J. Atmos. Sci., 73, 4329–4347, https://doi.org/10.1175/JAS-D-16-0033.1.
Bohren, C. F., and B. A. Albrecht, 1998: Atmospheric Thermodynamics. Oxford University Press, 198 pp.
Chang, E. K., 2009: Diabatic and orographic forcing of northern winter stationary waves and storm tracks. J. Climate, 22, 670–688, https://doi.org/10.1175/2008JCLI2403.1.
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.
Feldstein, S. B., 2002: Fundamental mechanisms of the growth and decay of the PNA teleconnection pattern. Quart. J. Roy. Meteor. Soc., 128, 775–796, https://doi.org/10.1256/0035900021643683.
Fletcher, C. G., and P. J. Kushner, 2011: The role of linear interference in the annular mode response to tropical SST forcing. J. Climate, 24, 778–793, https://doi.org/10.1175/2010JCLI3735.1.
Franzke, C., S. Lee, and S. B. Feldstein, 2004: Is the North Atlantic Oscillation a breaking wave? J. Atmos. Sci., 61, 145–160, https://doi.org/10.1175/1520-0469(2004)061<0145:ITNAOA>2.0.CO;2.
Garfinkel, C. I., S. B. Feldstein, D. W. Waugh, C. Yoo, and S. Lee, 2012: Observed connection between stratospheric sudden warmings and the Madden-Julian oscillation. Geophys. Res. Lett., 39, L18807, https://doi.org/10.1029/2012GL053144.
Goss, M., S. B. Feldstein, and S. Lee, 2016: Stationary wave interference and its relation to tropical convection and Arctic warming. J. Climate, 29, 1369–1389, https://doi.org/10.1175/JCLI-D-15-0267.1.
Held, I. M., and M. J. Suarez, 1994: A proposal for the intercomparison of the dynamical cores of atmospheric general circulation models. Bull. Amer. Meteor. Soc., 75, 1825–1830, https://doi.org/10.1175/1520-0477(1994)075<1825:APFTIO>2.0.CO;2.
Held, I. M., M. Ting, and H. Wang, 2002: Northern winter stationary waves: Theory and modeling. J. Climate, 15, 2125–2144, https://doi.org/10.1175/1520-0442(2002)015<2125:NWSWTA>2.0.CO;2.
Hoskins, B. J., and D. J. Karoly, 1981: The steady linear response of a spherical atmosphere to thermal and orographic forcing. J. Atmos. Sci., 38, 1179–1196, https://doi.org/10.1175/1520-0469(1981)038<1179:TSLROA>2.0.CO;2.
Hoskins, B. J., and P. J. Valdes, 1990: On the existence of storm-tracks. J. Atmos. Sci., 47, 1854–1864, https://doi.org/10.1175/1520-0469(1990)047<1854:OTEOST>2.0.CO;2.
Kobayashi, S., and Coauthors, 2015: The JRA-55 Reanalysis: General specifications and basic characteristics. J. Meteor. Soc. Japan, 93, 5–48, https://doi.org/10.2151/jmsj.2015-001.
Lee, S., 2012: Testing of the tropically excited Arctic warming mechanism (TEAM) with traditional El Niño and La Niña. J. Climate, 25, 4015–4022, https://doi.org/10.1175/JCLI-D-12-00055.1.
Lee, S., 2014: A theory for polar amplification from a general circulation perspective. Asia-Pac. J. Atmos. Sci., 50, 31–43, https://doi.org/10.1007/s13143-014-0024-7.
Ling, J., and C. Zhang, 2013: Diabatic heating profiles in recent global reanalyses. J. Climate, 26, 3307–3325, https://doi.org/10.1175/JCLI-D-12-00384.1.
Lorenz, E. N., 1955: Available potential energy and the maintenance of the general circulation. Tellus, 7, 157–167, https://doi.org/10.3402/tellusa.v7i2.8796.
Madden, R. A., and P. R. Julian, 1971: Detection of a 40–50 day oscillation in the zonal wind in the tropical Pacific. J. Atmos. Sci., 28, 702–708, https://doi.org/10.1175/1520-0469(1971)028<0702:DOADOI>2.0.CO;2.
Manabe, S., J. Smagorinsky, and R. F. Strickler, 1965: Simulated climatology of a general circulation model with a hydrological cycle. Mon. Wea. Rev., 93, 769–798, https://doi.org/10.1175/1520-0493(1965)093<0769:SCOAGC>2.3.CO;2.
Ming, Y., and I. M. Held, 2018: Modeling water vapor and clouds as passive tracers in an idealized GCM. J. Climate, 31, 775–786, https://doi.org/10.1175/JCLI-D-16-0812.1.
Overland, J., K. Wood, and M. Wang, 2011: Warm Arctic–cold continents: Climate impacts of the newly open Arctic sea. Polar Res., 30, 15787, https://doi.org/10.3402/polar.v30i0.15787.
Papritz, L., and T. Spengler, 2015: Analysis of the slope of isentropic surfaces and its tendencies over the North Atlantic. Quart. J. Roy. Meteor. Soc., 141, 3226–3238, https://doi.org/10.1002/qj.2605.
Peixoto, J. P., and A. H. Oort, 1992: Physics of Climate. American Institute of Physics, 520 pp.
Pierrehumbert, R. T., H. Brogniez, and R. Roca, 2007: On the relative humidity of the atmosphere. The Global Circulation of the Atmosphere, T. Schneider and A. H. Sobel, Eds., Princeton University Press, 143–185.
Smith, K. L., P. J. Kushner, and J. Cohen, 2011: The role of linear interference in northern annular mode variability associated with Eurasian snow cover extent. J. Climate, 24, 6185–6202, https://doi.org/10.1175/JCLI-D-11-00055.1.
Sutcliffe, R. C., 1951: Mean upper contour patterns of the Northern Hemisphere—The thermal-synoptic view pattern. Quart. J. Roy. Meteor. Soc., 77, 435–440, https://doi.org/10.1002/qj.49707733309.
Ting, M., 1996: Steady linear response to tropical heating in barotropic and baroclinic models. J. Atmos. Sci., 53, 1698–1709, https://doi.org/10.1175/1520-0469(1996)053<1698:SLRTTH>2.0.CO;2.
Ting, M., and I. M. Held, 1990: The stationary wave response to a tropical SST anomaly in an idealized GCM. J. Atmos. Sci., 47, 2546–2566, https://doi.org/10.1175/1520-0469(1990)047<2546:TSWRTA>2.0.CO;2.
Wheeler, M. C., and H. H. Hendon, 2004: An all-season real-time multivariate MJO index: Development of an index for monitoring and prediction. Mon. Wea. Rev., 132, 1917–1932, https://doi.org/10.1175/1520-0493(2004)132<1917:AARMMI>2.0.CO;2.
Willison, J., W. A. Robinson, and G. M. Lackmann, 2013: The importance of resolving mesoscale latent heating in the North Atlantic storm track. J. Atmos. Sci., 70, 2234–2250, https://doi.org/10.1175/JAS-D-12-0226.1.
Woods, C., and R. Caballero, 2016: The role of moist intrusions in winter Arctic warming and sea ice decline. J. Climate, 29, 4473–4485, https://doi.org/10.1175/JCLI-D-15-0773.1.
Woods, C., R. Caballero, and G. Svensson, 2013: Large-scale circulation associated with moisture intrusions into the Arctic during winter. Geophys. Res. Lett., 40, 4717–4721, https://doi.org/10.1002/grl.50912.
Wright, J. S., and S. Fueglistaler, 2013: Large differences in reanalyses of diabatic heating in the tropical upper troposphere and lower stratosphere. Atmos. Chem. Phys., 13, 9565–9576, https://doi.org/10.5194/acp-13-9565-2013.
Yoo, C., S. Lee, and S. B. Feldstein, 2012: Arctic response to an MJO-like tropical heating in an idealized GCM. J. Atmos. Sci., 69, 2379–2393, https://doi.org/10.1175/JAS-D-11-0261.1.
Zhang, K., W. J. Randel, and R. Fu, 2017: Relationships between outgoing longwave radiation and diabatic heating in reanalyses. Climate Dyn., 49, 2911–2929, https://doi.org/10.1007/s00382-016-3501-0.