The decline in Barents Sea ice has been implicated in forcing the “warm-Arctic cold-Siberian” (WACS) anomaly pattern via enhanced turbulent heat flux (THF). This study investigates interannual variability in winter [December–February (DJF)] Barents Sea THF and its relationship to Barents Sea ice and the large-scale atmospheric flow. ERA-Interim and observational data from 1979/80 to 2011/12 are used. The leading pattern (EOF1: 33%) of winter Barents Sea THF variability is relatively weakly correlated (r = 0.30) with Barents Sea ice and appears to be driven primarily by atmospheric variability. The sea ice–related THF variability manifests itself as EOF2 (20%, r = 0.60). THF EOF2 is robust over the entire winter season, but its link to the WACS pattern is not. However, the WACS pattern emerges consistently as the second EOF (20%) of Eurasian surface air temperature (SAT) variability in all winter months. When Eurasia is cold, there are indeed weak reductions in Barents Sea ice, but the associated THF anomalies are on average negative, which is inconsistent with the proposed direct atmospheric response to sea ice variability. Lead–lag correlation analyses on shorter time scales support this conclusion and indicate that atmospheric variability plays an important role in driving observed variability in Barents Sea THF and ice cover, as well as the WACS pattern.
Arctic sea ice has decreased in all seasons during recent decades. The largest declines in areal extent have occurred during summer and early autumn (up to 10% decade−1; Serreze et al. 2007), but the thicker multiyear ice cover is shrinking rapidly in winter as well (Comiso 2012). Wintertime Arctic sea ice area declines are concentrated in the Barents Sea (Serreze et al. 2009; Screen and Simmonds 2010a; Parkinson and Cavalieri 2012). The diminishing Arctic sea ice cover has been accompanied by near-surface warming in the high latitudes, particularly during the winter season (e.g., Screen and Simmonds 2010b), but the midlatitudes have recently experienced some anomalous winters with severe cold spells and above-normal snow cover in parts of Europe, Russia, and North America (e.g., Cohen et al. 2007; Cattiaux et al. 2010; Ghatak et al. 2010; Guirguis et al. 2011; Orsolini et al. 2012; Coumou and Rahmstorf 2012). The cold Eurasian surface temperatures are part of what has been termed the “warm-Arctic cold-Siberian” (WACS) anomaly pattern (or slight variations thereof; see Inoue et al. 2012; Mori et al. 2014). The relationship between this WACS pattern and Barents Sea ice conditions is the focus of this study, which probes reanalysis data for clues about how the two are linked.
Various mechanisms have been proposed to explain how the WACS pattern is generated. Early on, Honda et al. (2009) pointed to the fact that, where sea ice retreats, we expect enhanced ocean-to-atmosphere energy loss in the form of turbulent heat fluxes (THFs). They suggested that such a THF perturbation triggers a stationary Rossby wave train that amplifies the climatological-mean wintertime high pressure over Siberia and enhances northerly cold-air advection over eastern Eurasia. Other studies focus on the role of the THF perturbation in modifying the near-surface meridional temperature gradient, leading to altered cyclone pathways and a cold anticyclonic flow anomaly north of the Eurasian continent (Inoue et al. 2012), or weakened zonal winds and reduced heat transport from the North Atlantic to the interior of the Eurasian continent (Petoukhov and Semenov 2010; Outten and Esau 2012). A more global (less regional) alternative suggests that the modified (weaker) meridional temperature gradient causes the atmospheric flow to be wavier, thus leading to more extreme midlatitude weather in general (Francis and Vavrus 2012).
A different class of studies examines the phenomenon of cold Eurasian winters independently of Arctic sea ice. These studies attribute cold winters to persistent atmospheric blocking but not necessarily in connection with declining ice cover. Wintertime blocking patterns in recent cold years (e.g., 2005/06, 2009/10) are shown to be linked to sudden stratospheric warmings and/or sea surface temperature anomalies (Scaife and Knight 2008; Croci-Maspoli and Davies 2009; Cattiaux et al. 2010), with Rossby wave trains playing a possible amplifying role over Siberia (Takaya and Nakamura 2005; Park et al. 2011; Cheung et al. 2013; Sato et al. 2014). While such atmospheric circulation changes could be forced by sea ice changes (Mori et al. 2014), they could also be an expression of natural internal variability or reflect a direct atmospheric response to external radiative forcing.
The range of proposed mechanisms, including some that do not explicitly require a role for sea ice, has led to some debate over details of the relationship between the WACS pattern and sea ice cover. Some studies suggest that the WACS pattern is a delayed wintertime atmospheric response to reduced sea ice cover in the preceding autumn (Honda et al. 2009; Francis et al. 2009; Overland and Wang 2010; Overland et al. 2011) possibly via an indirect response involving the stratosphere (Jaiser et al. 2012; Orsolini et al. 2012; Cohen et al. 2014; Kim et al. 2014; García-Serrano et al. 2015; King et al. 2015). Others suggest that it is a simultaneous response to winter Barents Sea ice anomalies (Petoukhov and Semenov 2010; Hori et al. 2011; Inoue et al. 2012; Outten and Esau 2012). The link to winter sea ice anomalies has been shown to be more robust in at least one comparison (Tang et al. 2013). Furthermore, while some studies focus on interannual covariability in sea ice and Eurasian winters (Honda et al. 2009; Petoukhov and Semenov 2010; Hori et al. 2011; Inoue et al. 2012), others examine potential linkages to strong negative trends in sea ice cover (Outten and Esau 2012; Mori et al. 2014). Most recently, a series of observational and modeling studies has questioned whether there is adequate evidence that sea ice influences atmospheric circulation, blocking, and Eurasian winter temperatures at all (Hopsch et al. 2012; Barnes 2013; Screen and Simmonds 2013; Barnes et al. 2014; Wallace et al. 2014; Woollings et al. 2014).
Despite unresolved issues about the mechanism, seasonality, and time scale, many of the existing hypotheses attempting to explain the WACS pattern and its variability share a common starting point: an enhancement of ocean-to-atmosphere THFs in late autumn to winter as a direct consequence of reduced sea ice. These enhanced THFs are hypothesized to create a large-scale atmospheric response leading to the WACS pattern. This study will explore the robustness of this causal chain in reanalysis data by examining the THF covariability associated with both Barents Sea ice conditions and cold Eurasian winters.
Studies using AGCMs provide qualified support for the causal chain summarized in the previous paragraph. Movements of the sea ice edge are found to affect the spatial distribution of the THF field, producing a distinctive dipole signature with enhanced ocean-to-atmosphere THF in regions of ice loss next to the new ice edge and reduced THF next to the old ice edge (e.g., Alexander et al. 2004; Deser et al. 2004; Magnusdottir et al. 2004; Singarayer et al. 2006; Seierstad and Bader 2009; Budikova 2009; Deser et al. 2010). However, the large-scale atmospheric response to these THF anomalies is weak compared to natural variability (Screen et al. 2014); varies considerably among AGCM studies (Bader et al. 2011); exhibits nonlinear behavior (e.g., Lau and Holopainen 1984; Peng et al. 1997; Deser et al. 2007, 2010; Liptak and Strong 2014); and depends strongly on the magnitude, location, and polarity of the sea ice changes (e.g., Honda et al. 1999; Magnusdottir et al. 2004; Deser et al. 2004; Kvamstø et al. 2004), as well as the season (e.g., Alexander et al. 2004; Bhatt et al. 2008; Seierstad and Bader 2009; Deser et al. 2010).
Furthermore, AGCM experiments cannot capture the two-way interactions that link the ocean and the atmosphere in nature. The uncoupled AGCM setup allows us to investigate the atmospheric response to variability in the surface ocean (sea surface temperature and sea ice). It does not account for the other direction of coupling, known to be important on weekly to interannual time scales (Bjerknes 1964; Cayan 1992; Deser et al. 2000; Gulev et al. 2013), namely, the atmospheric driving of surface ocean variability via changes in near-surface winds, temperature gradients, and humidity gradients across the ocean–atmosphere interface (Fang and Wallace 1994; Prinsenberg et al. 1997; Deser et al. 2000; Rigor et al. 2002; Sorteberg and Kvingedal 2006; Wu and Zhang 2010; Germe et al. 2011; Wettstein and Deser 2014). As a result, uncoupled experiments often exhibit unrealistically strong and persistent ocean-to-atmosphere fluxes and reduced internal atmospheric variability (Barsugli and Battisti 1998). The THF anomalies in AGCM studies are by definition caused by the prescribed surface ocean conditions, and thus they cannot be used as evidence for the proposed causal chain.
In the current study, we explore the coupled interactions between Barents Sea ice cover, THF, and the large-scale atmospheric circulation during winter to further examine the causal chain linking Barents Sea ice conditions to the WACS pattern. We focus on THF in the Barents Sea, the region with the most pronounced winter sea ice variability, as the coupling agent between the ocean and atmosphere. Specifically, we explore whether Barents Sea THF variability is driven by sea ice or atmospheric variability in the observational record. Whereas previous observational studies investigated the WACS pattern in association with a few extremely low Barents Sea ice years (Honda et al. 2009; Hori et al. 2011; Inoue et al. 2012) or in association with trends in Barents Sea ice (Outten and Esau 2012; Outten et al. 2013), we focus on interannual variability during all years of the satellite record, excluding the linear trend. By removing the trend, we maximize the statistical degrees of freedom and isolate mechanisms that can produce the WACS pattern independent of lower-frequency climate change. This allows us to evaluate a candidate null hypothesis that the WACS pattern can be attributed at least in part to internal atmospheric variability. The corollary of our detrending decision is that the results of this study are not directly comparable to an examination of coupled interactions in the context of long-term, externally forced climate change.
2. Data and methods
We use data from the European Centre for Medium-Range Weather Forecasts interim reanalysis (ERA-Interim; Dee et al. 2011) from 1979/80 to 2011/12. We focus analyses on the 3-month winter season from December to February (DJF), but the results and conclusions are also valid over an extended winter season (October–February), a point we will return to in the discussion section. Monthly and daily THF, sea surface temperature (SST), surface air temperature (SAT), 10-m near-surface wind in the zonal (U) and meridional (V) directions, sea level pressure (SLP), and geopotential height at 500 hPa (Z500) are analyzed in the 1.5° × 1.5°-resolution ERA-Interim. The THF is defined as the sum of the surface latent and sensible heat fluxes. Positive THF values indicate heat transfer from the ocean to the atmosphere. Sea ice concentration (SIC) is estimated from passive microwave satellite data on a 25 km × 25 km grid (NSIDC; Cavalieri et al. 1996). This satellite-based SIC is independent from, but consistent with, SIC obtained from ERA-Interim for the Barents Sea as a whole (r = 0.98 for a raw correlation using monthly data, r = 0.97 when linearly detrended) and also in terms of the geographic distribution of the sea ice edge (Fig. 1). All results and conclusions are similar when SIC from ERA-Interim is used instead of the NSIDC data. We define two regions (both indicated in Fig. 6) for the analyses: the Barents Sea (67°–81°N, 12°– 73°E) and the Eurasian region (45°–81°N, 12°–120°E). The Eurasian region includes mid- and high-latitude continental Eurasia and the adjacent Arctic seas. Note that the Barents Sea region is wholly subsumed within the northwest quadrant of the larger Eurasian region. The results in this study are not very sensitive to the exact definition of the regions, several of which were tested and all of which produced qualitatively similar results. The Barents Sea ice index (ICEBar) is defined as the detrended area-averaged SIC over the Barents Sea region, with positive ICEBar corresponding to less Barents Sea ice.
Unless otherwise stated, results in this study are based on monthly anomalies calculated by removing the linear trend from each month individually and then removing the climatological-mean seasonal cycle from the detrended data. Standard empirical orthogonal function (EOF)/principal component (PC) analysis is performed on monthly anomalies of 1) THF over the Barents Sea and 2) SAT over the Eurasian region. Gridpoint values of the anomaly fields are weighted by the square root of the cosine of latitude (equivalent to area weighting the variance) before performing the EOF analysis. EOFs presented in this study are well separated from each other and lower-order EOFs according to the criterion described by North et al. (1982). The resulting spatial patterns are referred to as EOFs and the associated time series are referred to as PCs. Broader climate relationships are explored by regressing monthly anomalies of relevant variables onto standardized PCs and ICEBar, such that the resulting regression maps are displayed in physical units per standard deviation of the index.
Later in this study we create “synthetic” time series of the WACS pattern at daily (WACSdsyn) and seasonal (WACSssyn) time scales. In both cases, the synthetic time series is created by projecting Eurasian SAT anomaly patterns at these two time scales onto the second EOF of Eurasian DJF SAT variability (i.e., what is defined in the results section as the WACS pattern). For WACSdsyn, the linear trend and climatological daily mean are removed to create the anomaly patterns, and these are used in a variety of lead–lag correlation and regression analyses to provide additional context to our evaluation of the causal chain described in the introduction section. For WACSssyn, only the climatological seasonal mean is removed to create the anomaly patterns, and we assess the variance associated with the trend and in the linearly detrended WACSssyn time series.
The conclusions of the current study have been found to be robust when other satellite-based [OAFlux (Yu et al. 2008) and HOAPS (Andersson et al. 2011)] and reanalysis-based [NCEP1 (Kalnay et al. 1996)] THF products are used instead of ERA-Interim data. The OAFlux and HOAPS products are thought to give improved heat flux estimates compared to reanalysis products, but a disadvantage is that they mask out data in areas with sea ice concentration above 50%, which makes it impossible to examine comprehensively the associations between THF and SIC variability.
Because ocean-to-atmosphere THF has been widely invoked as the intermediary for communicating sea ice changes to the atmosphere (cf. references in the introduction section), we first examine the leading patterns of winter (DJF) THF variability and their association with sea ice cover in the Barents Sea. Figure 2 shows the leading Barents Sea THF PCs along with the ICEBar index defined in section 2. The first PC (PC1) is associated with 33% of the total DJF variability in Barents Sea THF, but it is relatively weakly correlated with the ICEBar time series (r = 0.30). The second PC (PC2) is associated with less of the variability (20%), but it is more strongly correlated with the ICEBar (r = 0.60). The relationships between Barents Sea THF PCs and ICEBar are robust not only for DJF but also throughout the extended [October–February (ONDJF)] winter season (Table 1).
The spatial patterns associated with ICEBar and the THF PCs support the idea that THF PC2 is more connected to sea ice cover than THF PC1. The THF signature of ICEBar is familiar (Fig. 3a)—strongly enhanced ocean-to-atmosphere THF over areas of reduced SIC accompanied by moderate reductions in THF directly to the south, indicating a northward shift in intense turbulent heat loss from the ocean associated with the sea ice edge (cf. Deser et al. 2004; Magnusdottir et al. 2004; Alexander et al. 2004; Deser et al. 2010). Recall that increases in ICEBar correspond to reductions in Barents-averaged SIC. Similar THF and SIC signatures appear in association with the second pattern of Barents Sea THF variability (THF EOF2; Fig. 3c), suggesting that THF EOF2 captures at least some aspects of the atmospheric response to sea ice reductions. The leading pattern (THF EOF1), however, tells a different story (Fig. 3b). THF EOF1 describes a broad, Barents-wide reduction in THF (reduced ocean-to-atmosphere heat loss in the yellow box) that is associated with relatively weak sea ice reductions. The THF–sea ice relationship (reduced THF with reduced sea ice) in Fig. 3b is incompatible with the expected relationship associated with an atmospheric response to sea ice reduction (enhanced THF with reduced sea ice).
The strength of winter THF depends on near-surface wind speeds and on the temperature difference between the (warm) ocean and (cold) atmosphere (Maykut 1982; Schröder et al. 2003; Alexander et al. 2004). If THF is enhanced as a result of a reduction in sea ice cover (as represented by ICEBar and THF PC2), then the ocean provides more heat to the atmosphere and we expect warm SAT anomalies. This relationship is shown in Figs. 3a and 4d and 3c and 4f for ICEBar and THF PC2, respectively, meaning that the variability described by the two indices is thus far consistent with an atmospheric response to Barents Sea ice perturbations. In contrast, both THF and sea ice cover are reduced with THF PC1 increases (Fig. 3b), while SAT anomalies are also warm (Fig. 4e). This suggests that the THF reductions associated with THF PC1 are likely driven by atmospheric changes, in particular a warming of the overlying atmosphere that reduces the ocean–atmosphere temperature difference. Oceanic (SST) correlations with THF PC1 are generally much weaker than the corresponding atmospheric (SAT) correlations (cf. Figs. 4b and 4e), which is to be expected if the ocean is responding to atmospheric changes. In comparison, THF PC2 exhibits spatially coherent oceanic and atmospheric patterns and slightly stronger correlations with SST than with SAT (cf. Figs. 4c and 4f), in concert with an atmospheric response to ocean changes, such as Barents Sea ice variability. Together, the analyses in Figs. 2–4 suggest an active role for the atmosphere in driving a substantial portion of THF variability over the Barents Sea.
Before exploring the atmosphere’s role further, we check the robustness of the sea ice–related THF variability (i.e., THF PC2) in individual winter months (Fig. 5). The localized dipole-like structure in THF appears in all months, and is consistently associated with enhanced ocean-to-atmosphere THF over areas of reduced SIC (Figs. 5a–c) and warmer SATs over the Barents Sea (Figs. 5d–f). On a hemispheric scale, the SAT regressions in all months show a generally cooler Eurasian continent similar to the WACS pattern described by Inoue et al. (2012), but the cooling is only statistically significant in December (Figs. 5d–f). In contrast, there is substantial variation in the large-scale atmospheric circulation patterns associated with THF PC2 (Fig. 5g–i). There are few robust features in the SLP regressions across all months, and we do not observe the expected linear atmospheric response to enhanced local heating—a negative SLP anomaly directly downstream of a heat source (e.g., Hoskins and Karoly 1981). We see some evidence of northerly cold-air advection over Siberia, in agreement with the circulation pattern described by Honda et al. (2009) and other studies listed in the introduction section. Similar results are obtained if the ICEBar index is used instead of THF PC2, particularly during midwinter.
The regression analyses in Fig. 5 show that sea ice–related THF variability is present in all winter months and is consistently associated with warming over the Barents Sea, but the large-scale circulation patterns vary substantially. As established in the introduction section, the variation from month to month could be due to nonlinearities in the dynamical response of the atmosphere or to internal variability, which is known to play a role even in the presence of strong forcing (Deser et al. 2004, 2012; Wettstein and Deser 2014). Internal variability is especially important on larger spatial scales; so even if the regional signature of THF PC2 can be understood as a simple atmospheric response to sea ice variability, hemispheric signatures like the WACS pattern seem to be more complicated.
Following this idea, we consider variability associated with the large-scale WACS pattern (Inoue et al. 2012; Mori et al. 2014) identified independently of Barents Sea ice variability. The WACS pattern emerges consistently as the second EOF (20%) of monthly Eurasian SAT variability in all winter months (Figs. 6d–f).1 The SAT and SLP regressions onto the WACS index are substantially stronger and generally more robust across all winter months than the comparable regressions onto THF PC2 (Figs. 5d–i) or ICEBar (not shown). In its positive phase, the WACS pattern exhibits robust out-of-phase SAT anomalies between warmer temperatures over the Barents–Kara Seas and colder temperatures over Eurasia. It is consistently associated with a prominent SLP high over the Eurasian continent and southerly near-surface wind anomalies over the Barents Sea from December to January (Figs. 6g–i). When Eurasia is cold, Barents Sea ice concentration is somewhat lower than average, but the associated Barents Sea THF anomalies are, on average, negative (Figs. 6a–c) for each winter month. This is the opposite of what would be expected if the WACS pattern were a direct atmospheric response to reduced Barents Sea ice. Thus, atmospheric variability appears not only to shape the WACS pattern itself, but also to contribute substantially to observed THF variability over the Barents Sea and large-scale SAT variability.
Causality is better explored by examining lead–lag relationships between the WACS pattern and atmospheric or oceanic indicators on shorter time scales. Figure 7 shows lagged correlations between a variety of indicators and a daily index of the WACS pattern (WACSdsyn as defined in the data and methods section). On daily time scales, the Barents-averaged THF (THFdBar; blue curve) leads the WACS pattern, with reduced ocean-to-atmosphere heat fluxes peaking 2 days before the WACS pattern. This echoes the message from the monthly analysis (Fig. 6), particularly that the WACS pattern does not appear to be a direct atmospheric response to reduced Barents Sea ice (and increased THF) on short (daily) or long (monthly) time scales. Note that because THF PC1 represents a basinwide pattern of THF variability in the Barents Sea, it is very similar to THFdBar except for a change in sign (see Fig. 3b). The lagged correlation curves associated with the daily version of these two indices are also very similar except for a change in sign (r = −0.96; see also Fig. 3b). On the other hand, the daily version of the THF PC2 index, which represents sea ice–related THF variability, exhibits low correlations with WACSdsyn at all lags (red curve). The lagged correlation between ICEdBar and WACSdsyn (black curve) peaks at zero lag, but it is stronger with ICEdBar lagging rather than leading, casting further doubt on the idea that the WACS pattern is a direct consequence of reduced Barents Sea ice.
Finally, southerly wind anomalies over the Barents region (VdBar; green curve) also peak 2 days before WACSdsyn. Supposing that these wind anomalies are associated with an anomalous atmospheric circulation pattern, this could suggest 1) advection of warm air from the south into the Barents Sea region, 2) a basinwide reduction of THF as a result of a correspondingly weakened temperature contrast between the warmer ocean and colder atmosphere, and 3) a more slowly evolving decrease in Barents Sea ice cover driven mechanically by southerly winds pushing the sea ice edge northward. This chain of events is consistent with the lead–lag relationships in Fig. 7 and suggests, along with previous results, that the WACS pattern does not originate directly from enhanced THF caused by Barents Sea ice reduction.
The corresponding spatial patterns at various lags (Fig. 8) provide some final clues about how the WACS pattern, THF, and atmospheric circulation over the Barents Sea are linked. The patterns are consistent with the results shown in Fig. 7: the WACS pattern is associated with negative Barents Sea THF anomalies (reduced ocean-to-atmosphere THF) and a relatively weak reduction in the Barents Sea ice cover (Figs. 8a–c), which together are inconsistent with the notion that the WACS pattern is a direct atmospheric response to reduced Barents Sea ice. The spatial distribution and diminishing intensity of THF anomalies over time (Figs. 8a–c) are instead consistent with atmospheric changes driving ocean-to-atmosphere THF variability (e.g., Bjerknes 1964; Cayan 1992; Gulev et al. 2013). The associated near-surface signature exhibits southerly and southwesterly winds (Figs. 8d–f), consistent with warm temperature advection over the Barents Sea. The large-scale surface (SLP) and midtropospheric (Z500) anomalies (Figs. 8g–i) resemble patterns associated with wave propagation and winter atmospheric blocking (e.g., Park et al. 2011; Cheung et al. 2013) and, to some degree, with the east Atlantic pattern (Wallace and Gutzler 1981; Smoliak 2009; Wettstein and Wallace 2010). It is also noteworthy that the strongest Z500 anomalies occur prior to the strongest (lag 0) correlation between the WACS pattern and Barents Sea ice reductions (cf. Figs. 7 and 8g–i), further supporting the interpretation that the large-scale atmospheric wave exists prior to any atmospheric changes initiated by Barents Sea ice (and the associated THF) variability. Our interpretation of the collective results is that the WACS pattern reflects a large-scale pattern of intrinsic coupled climate variability that is related to both the atmospheric forcing of and the atmospheric response to Barents Sea ice variability.
In the previous section, we investigated interactions between the observed Barents Sea ice cover and large-scale atmospheric circulation with a focus on ocean-to-atmosphere THF as the coupling agent between the underlying ocean and the overlying atmosphere. The goal was to evaluate a candidate null hypothesis that the observed interannual variability in the winter WACS pattern might primarily be an expression of atmospheric variability rather than an atmospheric response to sea ice variability. Our results suggest that we cannot reject this null hypothesis, and that atmospheric circulation variability does in fact play a critical role in determining the SLP, SAT, THF, and even the sea ice distribution both over the Barents Sea and farther afield.
A number of studies mentioned in the introduction section suggest that the Barents Sea ice variability in October and November (ON) causes a delayed response in the winter atmospheric circulation (Honda et al. 2009; Francis et al. 2009; Overland and Wang 2010; Overland et al. 2011; Jaiser et al. 2012; Orsolini et al. 2012; Cohen et al. 2014; Kim et al. 2014; García-Serrano et al. 2015; King et al. 2015). Because the previous section focuses on synchronous analyses during the DJF season, it does not address the possibility of this autumn-to-winter linkage. To explore this issue, Fig. 9 shows lagged regressions between the winter (DJF)-mean WACS index with THF, SIC, SST, and SLP anomalies in the months leading up to winter. The DJF WACS pattern is indeed associated with negative SIC anomalies (in the Barents, Kara, Laptev, and east Siberian Seas) and positive THF anomalies in the vicinity of the SIC anomalies during October and November, but many of these anomalies are insignificant.
The October–November regression results (Figs. 9a,b) are consistent with the suggestion that negative SIC/positive THF anomalies during autumn can influence the wintertime circulation but still do not provide adequate grounds for a direct causal interpretation. The geographical distribution of the fields (SIC, THF, SLP, and 10-m winds) is transitory over the late fall and early winter and the THF anomalies are displaced from the strongest SIC anomalies, suggesting an important role for nonlinear ice-circulation effects, internal atmospheric variability, or both. Many studies also point to more complicated, indirect pathways linking SIC and circulation via autumn stratospheric anomalies that propagate downward to create the tropospheric signals later in winter (Jaiser et al. 2012; Orsolini et al. 2012; Cohen et al. 2014; Kim et al. 2014; García-Serrano et al. 2015; King et al. 2015). While the downward propagation of stratospheric anomalies is clear in observations, the question of whether they are caused by sea ice anomalies is not yet settled (see discussion in King et al. 2015). Given the apparent complexity of the mechanisms and the relatively few degrees of freedom, a clear diagnosis of the causal linkage between autumn ice anomalies and midwinter WACS anomalies is perhaps unobtainable using reanalysis.
Sea ice loss is prominent in the observed record and is also consistently projected for the future within coupled GCM experiments forced by anthropogenic increases in greenhouse gases. The associated atmospheric response is reflected in various aspects of climate, both locally and remotely (e.g., Serreze et al. 2009; Simmonds and Keay 2009; Honda et al. 2009; Francis et al. 2009; Screen and Simmonds 2010a,b; Overland and Wang 2010; Overland et al. 2011; Blüthgen et al. 2012; Jaiser et al. 2012; Orsolini et al. 2012; Screen et al. 2013; Mori et al. 2014). To more completely evaluate the relationship between the WACS pattern, Barents Sea ice, and Barents Sea THF, we performed a full set of parallel analyses on raw (not detrended) variables.
The main conclusions of the current study remain when trends are retained (not shown). Consistent with Mori et al. (2014), we find trends toward a positive DJF WACS pattern and increased Barents Sea ice loss, both of which are particularly evident from 2004 onward. Attributing causality to this relationship is difficult, however, because of few degrees of freedom and because the correlation between the raw seasonal time series, r = 0.57 (0.47 for the monthly time series), represents substantially more shared variance than in the corresponding detrended analysis, r = 0.39 (0.33 for the monthly time series). More importantly, both synchronous and lagged regressions onto the WACS pattern are largely similar between raw and detrended analyses (not shown). In addition, an intraseasonal analysis like in Fig. 9 but performed on raw (not detrended) data returns very similar results to those shown in Fig. 9, suggesting the trend is perhaps not critically important in understanding the DJF WACS pattern.
An analysis of the total seasonally averaged variance associated with the linear trend and in linearly detrended data suggests that at least a third of the observed 1979/80 to 2011/12 variance in the winter Barents Sea ice is associated with the linear trend (Table 2). The trend accounts for a small fraction of the variance in the DJF WACS index (8%) and hardly any of the variance in winter Barents Sea THF (1%). The mismatch in the relative importance of the trend further suggests it is unlikely that the WACS pattern results from a direct linear atmospheric response to sea ice loss or variability. In short, our decision to explore the interannual variability allows for a focus on the mechanistic relationships between the atmospheric circulation anomalies, sea ice variability, and the connection of both to the upper ocean using the maximum number of degrees of freedom.
The goal of identifying the atmospheric response to Barents Sea ice variability in the fully coupled system motivated our study, and it is critical for understanding our results in the context of previous work. Studies based on observations find compelling correlations between sea ice and atmospheric circulation variability (e.g., Honda et al. 2009; Hori et al. 2011; Inoue et al. 2012; Outten and Esau 2012; Mori et al. 2014), but causality is challenging to determine in nature because of the many complicated forcing/response interactions and the relatively few years of data. Studies based on uncoupled modeling experiments do isolate an atmospheric response by design—they prescribe an ice anomaly under an atmospheric model and examine the resulting changes to the atmospheric circulation. While there is undoubtedly a WACS-like response to prescribed Barents Sea ice reductions, the particular result is somewhat model dependent (Honda et al. 2009; Petoukhov and Semenov 2010; Orsolini et al. 2012; Mori et al. 2014). It is also unclear how much of the atmospheric response is due to the prescribed sea ice loss as compared with other surface (SST, snow cover) changes, and how much is directly attributable to changes in radiative forcing. Further complicating the task is the presence of strong natural variability and rich two-way ocean–atmosphere interactions in nature (Deser et al. 2004, 2012; Wettstein and Deser 2014; Screen et al. 2014). Our focus on the surface turbulent heat fluxes allows us to use the reanalysis-based data, which represent the fully coupled system, without having to assume that observed atmospheric circulation changes are purely a response to sea ice forcing. Our results are generally more consistent with studies attributing cold Eurasian winters to atmospheric circulation variability (e.g., Hurrell 1995; Jeong and Ho 2005; Zhang et al. 2008; Cattiaux et al. 2010; Park et al. 2011; Sato et al. 2014) than with those attributing them to Arctic sea ice variability. Relationships between ocean circulation, SST, Barents Sea ice, and Barents Sea THF (Årthun et al. 2012; Sato et al. 2014) are certainly intriguing, but cleanly assessing their relative importance in a coupled WACS response is inherently difficult with low-frequency variability and so few years of data.
The implications of this observational study must be viewed strictly in the context of the fully coupled climate system on interannual time scales. In this context, the ocean, ice, and atmosphere are continuously interacting and none of the three components is obviously dominant. Atmospheric circulation variability is driven by many factors, sea ice being just one of these (see also Deser et al. 2012; Wettstein and Deser 2014). The results collected here demonstrate that the proposed direct link between variability in Barents Sea ice and the WACS pattern is at best an incomplete characterization of nature. Our study suggests that a broader and more rigorous diagnosis of the causes of Eurasian continental temperature extremes is warranted.
Our study examines relationships between Barents Sea ice, turbulent heat flux (THF), and the “warm-Arctic cold-Siberian” (WACS) pattern. Our main findings are as follows:
The leading pattern (EOF1) of winter Barents Sea THF variability appears to be driven primarily by atmospheric rather than sea ice variability. The sea ice–related THF variability appears as EOF2.
The sea ice–related THF variability (EOF2) is robust over the entire winter season, but its link to cold anomalies over Siberia is not.
The WACS pattern is a robust pattern of temperature variability (EOF2 of winter Eurasian SAT) and is associated with weak reductions in Barents Sea ice cover but also reduced ocean-to-atmosphere THF. The sense of these relationships is inconsistent with the proposed pathways of the WACS pattern being a direct atmospheric response to Barents Sea ice reductions.
The role of Barents Sea ice variability in forcing the WACS pattern appears to be minor in the detrended analyses performed here. One consistent interpretation of the results is that large-scale atmospheric circulation variability produces southerly wind anomalies that push the Barents Sea ice edge northward, advect warm air into the region, and reduce Barents Sea THF, and that are ultimately associated with cold Siberian temperatures downstream. The analyses leave open the possibility that Barents Sea ice variability drives an atmospheric response that feeds back on preexisting atmospheric circulation anomalies.
Although CL and JJW are listed alphabetically, they contributed equally to the manuscript. We wish to thank J. Kristiansen, P. Kushner, three anonymous reviewers, and the editor for their helpful comments. We acknowledge the European Centre for Medium-Range Weather Forecasts for providing the ERA-Interim data and the National Snow and Ice Data Center for providing the sea ice data. This work was supported by the NordForsk TRI program under project GREENICE 61841 (CL, JJW), the Research Council of Norway under projects CLIMARC 225032 and Bjerknes Compensation 227137 (SAS), and the Centre for Climate Dynamics (SKD) at the Bjerknes Centre (CL). JJW acknowledges support from the G. Unger Vetlesen Foundation.
The first EOF (40%) of Eurasian SAT is associated with the temperature signature of the North Atlantic/Arctic Oscillation (cf. Thompson and Wallace 1998).