Analyzing links between simulated Laptev Sea sea ice and atmospheric conditions over adjoining landmasses using causal-effect networks
We investigate how sea ice interacts with the atmosphere over adjacent landmasses in the Laptev Sea region as a step towards a better understanding of the connection between sea ice and permafrost. We identify physical mechanisms as well as local and large-scale drivers of sea-ice cover with a focus on one region with highly variable sea-ice cover and high sea-ice productivity: the Laptev Sea region. We analyze the output of a coupled ocean–sea-ice–atmosphere–hydrological-discharge model with two statistical methods. With the recently developed causal-effect networks we identify temporal links between different variables, while we use composites of high- and low-sea-ice-cover years to reveal spatial patterns and mean changes in variables.
We find that in the model local sea-ice cover is a driven rather than a driving variable. Springtime melt of sea ice in the Laptev Sea is mainly controlled by atmospheric large-scale circulation, mediated through meridional wind speed and ice export. During refreeze in fall thermodynamic variables and feedback mechanisms are important – sea-ice cover is interconnected with air temperature, thermal radiation and specific humidity. Though low sea-ice cover leads to an enhanced southward transport of heat and moisture throughout summer, links from sea-ice cover to the atmosphere over land are weak, and both sea ice in the Laptev Sea and the atmospheric conditions over the adjacent landmasses are mainly controlled by common external drivers.
To better understand both the mechanisms behind and the strength of the interaction between sea ice and land, we explore links between sea ice and the atmosphere over land and identify local and large-scale drivers of sea-ice cover in the Laptev Sea. Sea ice interacts with the atmosphere on different scales. However, while links from sea ice to large-scale atmospheric processes have been shown (e.g., Simmonds, 2015; Luo et al., 2017; Screen et al., 2018; Samarasinghe et al., 2019), the strongest coupling to the atmosphere is local (Screen and Simmonds, 2010; Screen et al., 2013). Sea ice influences near-surface temperatures by changing the local energy budget and regulating the moisture and energy which enter the lower atmosphere (Screen and Simmonds, 2010; Screen et al., 2013). This effect is more predominant in fall than in spring (Serreze et al., 2009; Serreze and Barry, 2011; Screen et al., 2012). Additionally, downward radiation plays a role in changing the net surface fluxes and thereby the surface temperature. Downward radiation has been associated with moisture fluxes from midlatitudes into the Arctic. The moisture fluxes show a positive trend in recent decades (Serreze and Barry, 2011; Lee et al., 2017). Little attention has been focused on the physical mechanisms through which variability in sea ice influences the atmosphere over land. Nevertheless, from prior research we know that sea ice can exert such an influence on land (Lawrence et al., 2008; Ogi et al., 2016). Changes in the atmosphere over land which are attributed to declining sea ice lead to various responses in the permafrost landscapes, ranging from increased methane emissions (Parmentier et al., 2013, 2015) to changes in vegetation productivity (Bhatt et al., 2008; Macias-Fauria et al., 2017) and in vegetation composition (Post et al., 2013). Thus, a better understanding of the connection between sea ice and land is valuable, especially since sea ice and the permafrost covering adjacent landmasses are both highly vulnerable to climate change. In this paper, we aim for a better understanding of the physical mechanisms behind the connection of sea ice to the atmosphere over land.
As mentioned, sea ice has a strong impact on the energy balance of the ocean surface, giving rise to several feedback mechanisms, such as lapse-rate feedback (Pithan et al., 2013; Pithan and Mauritsen, 2014), water–vapor feedback (Francis and Hunter, 2007) or ice–albedo feedback, which has been identified as a major control on Arctic temperatures (Deser et al., 2000; Serreze et al., 2009; Serreze and Barry, 2011; Graversen et al., 2014). These interconnections all contribute to a strong seasonality and interannual variability of sea ice, especially in the marginal seas of the Arctic Ocean (Deser et al., 2000). Because different processes might overlay each other when looking at the Arctic as a whole, we focus on one region where we expect a comparably strong connection of sea ice and the adjacent land: the Laptev Sea is one of the key contributors to net sea-ice production in the Arctic (Bareiss and Görgen, 2005; Bauer et al., 2013) and shows large year-to-year variability (Haas and Eicken, 2001) (Fig. 2). To identify the main processes between sea ice and the atmosphere over permafrost, we include as many variables as possible in our analysis (for an overview of all included variables, see Fig. 1 and Table A1). Because observations are sparse in space and time and not available for all relevant variables, we use model output. The big advantage of such a model over observational studies is that we can analyze a very large range of variables in a physically consistent system and in high resolution, both spatially and temporally. In contrast to reanalysis, we can run the model with the same forcing several times and can thus produce more data of stable climatic conditions. We are additionally able to compare different timescales and analyze the interactions on a monthly and daily scale. Previous studies have shown that unusually strong storm activities can change the state of Arctic sea ice in the long run (Screen et al., 2011; Simmonds and Rudeva, 2012, 2014). Such features on short timescales, for example the appearance of cyclones, can not be seen in an analysis based on monthly means only. However, for all model and reanalysis studies, it is important to keep in mind that the knowledge we can gain from looking at a larger scale is only as good as our understanding of the underlying process we depict in the model.
We use forcing from the time period of 1950 to 1989. In this period and in the Laptev Sea, we do not yet observe a general downward trend of sea ice. We run the model repetitively to improve statistical power. On the thus obtained 160 years of model output we employ two statistical methods. The first is called causal-effect networks, a recently developed method (Runge et al., 2012, 2014, 2015), which has been successfully applied by Kretschmer et al. (2016) to analyze Arctic drivers of midlatitude winter circulation. This method allows us to (a) identify the important links in an unbiased way and (b) differentiate whether two variables are subject to a common external forcing or which variable is forcing the other. Building on the results from the causal-effect networks, we group model years with exceptionally high and low sea-ice cover. These composites reveal spatial patterns and mean changes in variables, allowing us to gain a deeper physical understanding.
In the following, we start with introducing the causal-effect networks and the composite analysis. Then, we analyze the impact of sea-ice cover on the atmosphere over land as well as the drivers of sea-ice cover during the onset of the melting season and during refreeze. Finally, we put our results in the wider context.
In order to understand links between Laptev Sea sea ice and regional climatic conditions, we analyze the output of a regionally coupled ocean–sea-ice–atmosphere–hydrological-discharge model with two complementary methods. The model consists of the global ocean–sea-ice model MPIOM (Jungclaus et al., 2013) coupled to the regional atmosphere model REMO (Jacob and Podzun, 1997; Jacob, 2001), which in our setup covers most of the Northern Hemisphere. The atmosphere model has a horizontal resolution of approximately 55 km and 27 vertical levels. In the global ocean model, the grid poles are located over North America and Russia, leading to a horizontal resolution of up to 5 km within the Arctic Ocean. The model has 40 vertical levels with varying depth. The regional atmosphere model and the ocean model in the uncoupled domain are forced with model output from the global model MPIOM and ECHAM5 for the time period 1950–1989 and the model was run repetitively. Details on the model setup can be found in Sein et al. (2015) and Niederdrenk et al. (2013) and on the experimental design in Niederdrenk et al. (2016).
The model simulations have been validated against observations by Niederdrenk et al. (2013) and show a realistic mean Arctic climate for this time period. Also, the variability in sea-ice extent and thickness is captured well, being on the lower edge of observations. Due to its high resolution the regional model simulates more realistically than a global model the sea-ice transport within the Arctic (Niederdrenk et al., 2016). For our analysis, we run 40 years of forcing repetitively, so we can use 160 years of model output in total. Because the model melts ice directly from the ice edge, as it is not able to simulate landfast ice and polynyas well, it shows less ice in the Laptev Sea than observations (see Fig. 4). Nevertheless, the variability of the model lies within the observational records. For the chosen time period the output does not show a drift in sea-ice cover in the Laptev Sea (see Fig. 4) as sea-ice decline accelerated only in the 1990s. Thus, we can use the model output to examine underlying links between the climate and Laptev Sea sea ice assuming no or only little disturbance through model drift or anthropogenic forcing. Sea-ice cover is the quantity we use to represent sea ice in our analysis since we expect that sea-ice cover has a larger impact on atmosphere–ocean exchange processes than, e.g., sea-ice volume. Sea ice covers the Laptev Sea completely during winter, and most variability can be observed during the summer months (see Fig. 2). We therefore focus on the summer season to understand what influences sea-ice cover and what is influenced by sea-ice cover.
As a first method we use causal-effect networks, implemented in the Python package TiGraMITe (Runge et al., 2012, 2014, 2015). Causal-effect networks are an algorithm for causal discovery: the algorithm finds causal links in a dataset without a priori knowledge on physical mechanisms. Causal-effect networks determine how a perturbation moves through a set of one-dimensional time series. Each time series represents one variable, for example precipitation, and the temporal path of the perturbation is considered a causal link. The procedure to find links and gauge their strength is divided into two steps. In the first, relevant causal and contemporary links for each time series are identified. In the second step, the strength of these links is quantified. To identify the links of a target time series, the correlation between this target time series and, one after the other, all other potentially driving time series are evaluated. For each variable, first the direct correlation is computed and then, in an iterative manner, the partial correlation by including all possible other time series. For each target time series, all potentially driving time series are also shifted back in time, as the signal in the driven time series will lag behind the signal in the driving time series. This shift is increased from one time step to a maximum time lag τm. Only if a link remains significant no matter which subset of time series was included in the correlation analysis, we add a time series to the list of preliminary drivers of the target time series. This procedure is repeated for all time series in the dataset, so that we know all preliminary drivers for all time series. In a second step, these preliminary drivers are used to reevaluate the link strength between each pair of time series by applying a multiple linear regression. We compute the multiple linear regression between a time series; the preliminary drivers of this time series; and, iteratively, one other time series at, one after the other, all possible time lags. The results of all multiple linear regressions are then summarized in a causal regression matrix of dimension (N, N, τm), where N is the number of all time series we use and τm is the maximum time lag. Values in this matrix range from zero to one and are similar to a partial correlation coefficient. The threshold when a link is considered significant was set to 0.2, such that in our networks links that are roughly 2 standard deviations stronger than the mean link between all time series are considered significant. If we find a lagged link between two time series, we call it “causal”. We also consider “contemporary” links, where the time lag is zero. For these links, it is not possible to determine the direction of information flow. Note that the concept of causation in causal-effect networks is related to Granger causality (Granger, 1969), which tests whether it is possible to predict the future development of one time series from the past development of another time series (Runge et al., 2014). Additionally, the first step of this algorithm is an adapted PC algorithm (Spirtes et al., 2000) based on the idea that to identify causation we need to exclude common drivers (Pearl, 2000; Runge et al., 2012). Compared to many other frameworks making use of Granger causality, the above algorithm is very computationally efficient even on high dimensions (large N⋅τm) because using results of the first step drastically reduces the complexity of the second step (Runge et al., 2019).
The method becomes more reliable when the set of variables is large. To get a complete picture, we select all variables given in Fig. 1 and Table A1. We spatially reduce the data by either computing averages or sums over the area of interest (red box in Fig. 2). If the variable has a “per-area” unit, we weight the sum using the grid-cell area. Most of the variables are nondirectional scalar fields, but we split vector fields, such as wind, into their zonal (positive sign: west to east) and meridional (positive sign: south to north) component and reduce the components separately. There are not only horizontal but also vertical fluxes. These we either integrate vertically over the atmosphere or only consider surface fluxes from the atmosphere to land, ocean or ice. Lastly, for ice transport, we compute the gross import and export out of and into the Laptev Sea region. The boxes for the atmosphere and ocean model are given in Fig. 3.
To resolve effects on different scales, we analyze both monthly and daily means. When analyzing monthly-mean time series, we consider the past year (τm=12 months) as maximal time lag, while we allow a lag of 1 week (τm=7 d) in the daily setup. We look at the connection between land and sea ice especially during June–September when vegetation is photosynthesizing and sea-ice cover is low and variable. This variability accentuates the differences between high and low sea-ice-cover years, which is important for the composite analysis. We additionally investigate what drives sea-ice evolution. To do so, we investigate the phases of strong changes in early ice melt and ice refreeze in spring and fall. While the Laptev Sea is still completely ice covered in March in most years, sea-ice cover starts declining in April in most years, so we consider the months April, May and June for ice melt. Starting in September, ice cover grows from its minimum extent to a nearly complete cover in November. Consequently, for the ice refreeze we focus on the months of September, October and November. Since our goal is to connect land and sea-ice cover, we only consider the atmosphere over landmasses adjacent to the Laptev Sea (see brown area in Fig. 3) for the summer period (June–September). In contrast, for early ice melt (April–June) and ice refreeze (September–November), we also want to understand what drives sea ice. So for these two cases, we look at the atmosphere over land and ocean together (entire shaded area in Fig. 3).
With causal-effect networks, we can disentangle temporal relationships. To also understand spatial patterns and quantify dependencies, we extend our analysis by using composites of years with exceptionally high and low summer sea-ice cover in the Laptev Sea between June and October. These months coincide with the months we chose for the summer causal-effect network but are extended by October to include more of the refreezing phase, which only starts sometime during September. We select all years which deviate from the mean sea-ice state by more than 1.3 standard deviations. We select 1.3 standard deviations as a threshold (a) to only use years with significant (larger than 1 standard deviation) anomalies in sea-ice cover and (b) to even out the number of years contributing to the two composites as well as possible. With 1.3 standard deviations we get 17 years with high sea-ice cover and 7 with low (see Fig. 4).
As our main focus lies on the interaction between sea ice and land, we start our analysis by looking at the summer months when sea-ice variability is largest. Afterwards we will cover spring and fall to find drivers of the melt onset and the refreezing of ice.
To understand how sea-ice cover influences the adjacent land, we first look at causal-effect networks during the summer months. Here, both on a monthly and on a daily scale, we only consider the atmosphere over land and not over the ocean. In this way, direct interactions between the sea surface and the atmosphere are excluded. In Fig. 5a, we see the results of the causal-effect network analysis using monthly means throughout the summer months. Links with a strength below a threshold of 0.2 are neglected as insignificant in our analysis. Out of the 16 variables used in the setup, 4 are connected to sea-ice cover. The strength of the links to sea ice is weak (average strength: 0.25): the (solely contemporary) connections the variables have between each other are much stronger (mean strength: 0.57). For clarity, these connection are not shown in the figure. There is only one causal link, from meridional wind speed to sea-ice cover, while all other links are contemporary. All these links have a negative sign: higher sea-ice cover is associated with a decrease in the strongly interconnected variables specific humidity and thermal downward radiation at the surface, as well as with 2 m air temperature and 10 m meridional wind speed. Sea-ice cover has fewer links on the daily scale (see Fig. 5b) than on the monthly scale. It only links with meridional wind speed, and the only causal link is pointed towards, not away from, sea-ice cover, meaning that sea-ice cover is the driven variable.
If we investigate wind with the composites by looking at mean wind speeds over the Laptev Sea region, we see that wind direction is variable during the summer months, and there is a slight excess of southward wind. While in the mean about 56 % of the days during June to September show winds pointing southward, about 53 % of the days in the low-sea-ice-cover years and 64 % of the days in the high-sea-ice-cover years have a mean southward wind direction: high sea-ice cover coincides with a higher fraction of southward wind and vice versa for low sea-ice cover. More precisely, for high-sea-ice-cover years southward transport of moisture and heat is enhanced (see Fig. 6 and Table 1). The median of the net flux is lower compared to the mean, and this southward shift even leads to a change in sign for moisture transport from net-northerly to net-southerly transport compared to the mean. For low-sea-ice-cover years, on the other hand, both northward and southward transport of moisture and heat increases. Again the median is lower than the median of all years for net fluxes. This indicates that the bulk of transport also shifts towards more southward transport, but, for low sea-ice cover, there is a higher spread in the distribution visible in both the net and the gross fluxes. Using a two-sample Kolmogorov–Smirnov test we find that for both moisture and heat transport the high-sea-ice-cover and low-sea-ice-cover distributions are significantly different (p values < 0.001).
3.2 Early melt season of sea ice
Looking at monthly drivers of the early melt of sea ice in spring and including both the atmosphere over the ocean and over land, we find in total 10 variables that have links to sea-ice cover (Fig. 7a). The gross of connections is contemporary (nine variables). Two causal links are found, both of them directed towards sea-ice cover and with a lag time of 1 month. The lagged (causal) links are on average weaker than the contemporary ones. Returning to the contemporary links, the strongest links connect sea-ice cover to meridional wind speed, sea-ice export and meridional moisture transport, respectively.
On a daily scale (Fig. 7b), zonal wind speed and sea-ice export are still two of the three variables with the strongest connection to sea-ice cover. The third strongest connection now is with sea-ice import. These three directional variables also have causal links influencing sea-ice cover with a time lag of 1 and 2 d. Two other variables with causal links pointing towards sea-ice cover are the zonal vector components of wind speed and of moisture transport. Both have a 1 d lag to sea-ice cover, and there is no contemporary link between sea-ice cover and either of the two. On a monthly scale, they have neither a contemporary nor a causal link. Links between sea-ice cover and sea-level pressure, surface longwave downward radiation, and specific humidity on the other hand are significant on a monthly but not on a daily scale.
Since the strongest connections that sea-ice cover exhibits during early ice melt are connected to meridional transport and wind, we inspect spatial patterns using the composites for sea-level pressure in May (Fig. 8): we generally have a high over Greenland, but it is stronger for high sea-ice cover than for other years. For high sea-ice cover, we also observe a more distinct polar high, as opposed to low sea-ice cover years, which are characterized by a strong high-pressure zone centered over the Chukchi Sea as well as a low over the Kara Sea. Locally over the Laptev Sea, these differences on the large scale lead to contrasting pressure fields. As air flows around pressure systems, we also have vastly different wind patterns in our area of interest. For high sea-ice cover, wind in the Laptev Sea predominantly blows from the northeast, carrying air from similar or higher latitudes. For low sea-ice cover on the other hand, we have strong transpolar wind caused by the bipolar pressure field. Consequently, we also have a pronounced transpolar drift of sea ice originating roughly in the Laptev Sea and pointing towards Greenland. Over the Laptev Sea, the average wind direction is pointing northward for low sea-ice cover, and the air masses arrive from the southeast, thus carrying comparably warmer air. In the monthly ice-melt causal-effect network (Fig. 7a), this warm air is reflected in a contemporary link between air temperature and sea-ice cover.
3.3 Refreezing of the sea in fall
The link between air temperature and sea-ice cover is also present during refreeze, when this is the strongest link of sea-ice cover to any variable (Fig. 7c) in the causal-effect network. As for spring, it was set up to include the atmosphere over both ocean and land. During ice refreeze, the links from sea-ice cover to sea-ice import and export as well as meridional wind speed are weaker than during the melt onset. Heat and moisture transport do not have significant links to sea-ice cover. The autocorrelation of sea-ice cover, on the other hand, increased from 0.39 to 0.49. The auto-links of air temperature (strength: 0.40) and snow height (strength: 0.53) are also exceptionally strong compared to the other displayed variables. The total number of variables that are significantly connected to sea-ice cover decreases. The number drops from 10 during ice melt to 8 during refreeze. Apart from meridional heat and moisture transport, the sensible surface heat flux and sea-level pressure are replaced by snow height and the latent surface heat flux. Both of these variables have only causal and no contemporary links to sea-ice cover. The latent surface heat flux on the other hand is a dependent variable, with a time lag of 3 months. The lower sea-ice cover is, the higher the absolute fluxes from ocean to atmosphere are. The dependency of sea-ice export on sea-ice cover also changes: the contemporary and the causal link from sea-ice export to sea-ice cover during ice melt are replaced by a causal link from sea-ice cover to sea-ice export – more sea-ice cover leads to more export in the next month. In general, the number of causal links increases in the ice refreeze monthly-mean causal-effect network compared to melt. In addition, during refreeze sea-ice cover is also a preceding variable – during melt sea-ice cover is always the driven variable.
4.1 Drivers of sea-ice cover
Tracing influences of sea-ice cover on the atmosphere over land with causal-effect networks is not straightforward: sea-ice cover couples more strongly to atmospheric variables integrated over both land and ocean (Fig. 7a, c) than to the atmosphere over land (Fig. 5). For summer months, where we look solely at the atmosphere over land, we detect no causal link from sea-ice cover to any variable. The contemporary links that sea-ice cover has are weaker and smaller in number than in the causal-effect networks of melt and refreeze, where we included the atmosphere over the ocean as well: sea ice has a stronger impact on the atmosphere directly over the ocean. This is in line with prior findings that the impact of sea ice on the atmosphere is mostly local (Screen et al., 2012, 2013).
A causal link we still observe is from meridional wind speed to sea-ice cover. Moreover, meridional wind speed is the only variable connected to sea-ice cover in the daily-mean summer causal-effect network, when only the atmosphere over land was considered. During early melt, when using both the atmosphere over land and over the ocean, we observe the same dependence of sea-ice cover on meridional wind speed.
In spring, we find a counter-intuitive link from air temperature to sea-ice cover (higher temperatures lead to more sea ice, Fig. 7a), which is most likely an artifact of the causal-effect network analysis because of our choice of time series: we include downward longwave radiation in our analysis, which is strongly connected to air temperature. But higher air temperatures coincide also with increased upward longwave radiation, which leads to an energy loss at the surface. This upward longwave radiation was not included in the causal-effect network. We think that within the analysis the positive link from temperature to sea ice is caused by an overcompensation: the additional energy available at the surface through downward radiation is subtracted from the influence of air temperature on ice, but the energy radiated upward is not. This only becomes important for the causal link rather than the contemporary connection because for causal links we can subtract the influence of downward radiation at the present and prior time steps reinforcing the overcompensation.
Besides the link from air temperature to sea-ice cover, variables connected to the pressure field have the strongest links to sea-ice cover on both a monthly and daily scale. Those include local sea-level pressure and particularly the north/south components, like meridional wind speed, but also meridional moisture transport and sea-ice export. These variables are all strongly interconnected with each other on the monthly scale (Fig. 7a). While the variables have only contemporary links to sea-ice cover on monthly scale, they are driving sea-ice cover with a lag of 1 to 2 d in the daily-mean causal-effect network (Fig. 7b). This is in line with a prior observation-based study by Krumpen et al. (2013), in which sea-ice export from February to May has been singled out as one driver of sea-ice cover. The main local variable determining sea-ice cover in the Laptev Sea during early melt is thus meridional wind, which is also one main driver of sea-ice export. The stronger the northward flow, the less sea-ice cover. We can support this with our findings in the May composites (Fig. 8): for high-sea-ice-cover years, wind blows from the northeast, carrying colder, less humid air, and, for low-sea-ice-cover years, winds are blowing from the southeast, enhancing transpolar wind and sea-ice export and bringing warmer, more humid air. (Meridional) wind depends on air pressure. The large-scale patterns of the May pressure field resemble the patterns of previous, mostly observational or reanalysis-based, studies (e.g., Wang et al., 2009; Jaiser et al., 2012; Overland et al., 2012). For high sea-ice cover, similar to previous results from Luo et al. (2019b), sea-level pressure patterns resemble the negative phase of the Arctic Oscillation index (AO) – a high over the central Arctic Ocean (Wang et al., 2009). We observe a pronounced high over Greenland. This might hint to Greenland blocking, an event which has been linked to a negative North Atlantic Oscillation index (NAO) (thus also a negative AO). Greenland blocking has been found to cause strong melt of the Greenland ice sheet (Overland et al., 2012) as well as an enhanced influx of moisture in the North Atlantic region (Yang and Magnusdottir, 2017). This influx then enhances the water–vapor feedback. Instead of a positive AO pattern, low-sea-ice-cover years match an Arctic Dipole (AD) pattern (Fig. 8, similar to Wang et al., 2009; Overland et al., 2012): we observe a high over the Chukchi and Beaufort seas with even significantly higher pressure south of the Chukchi Sea and a low over the Barents Sea. The occurrence of the AD pattern fits well with the causal-effect-network-based result: that meridional wind is strongly connected to sea-ice cover in the Laptev Sea. We have pronounced sea-ice export because the AD pattern leads to anomalies in the meridional wind and a stronger transpolar drift of sea ice, while the AO is more connected to zonal wind promoting either convergence or divergence of sea ice in the Arctic (Wang et al., 2009). Advection of Siberian air masses has been linked to extremely low sea-ice cover in the Laptev Sea before, for example by Haas and Eicken (2001). Apart from enhanced ice drift, the authors stress the importance of atmospheric heat content carried into the Laptev Sea. Additionally, temperature fluctuations caused by advection have been identified as a dominant driver of sea-ice variability by Olonscheck et al. (2019). This agrees well with our findings since, during melt, meridional wind, air temperature and meridional heat transport all have strong links to sea-ice cover. In contrast, we find that the moisture content of the air transported into the Laptev Sea region plays a larger role than the heat content of the air. This is in agreement with a satellite- and reanalysis-based study by Yang and Magnusdottir (2017), who found that influxes of humid air into high latitudes from the North Atlantic lead to extremely low sea-ice cover in the Greenland, Barents and Kara seas.
We conclude that, during the melt phase, large-scale circulation systems exert a strong influence on the evolution of sea-ice cover. Since melting continues through summer (Kim et al., 2016) until September, the prevailing link from meridional wind speed to sea-ice cover portrays this connection.
During refreeze in fall, meridional wind speed is still important (Fig. 7b) but loses its dominance. While the connection between meridional wind speed and sea-ice cover in the monthly-mean causal-effect network for summer is on the same magnitude as the other connections of sea-ice cover, meridional wind speed has the weakest connection to sea ice among the significant variables during refreeze. Instead, state variables connected to feedback mechanisms, like ice–albedo feedback or water–vapor feedback, are increasingly important and also strongly interconnected. We can already identify these feedback mechanisms during early melt, even though they are less dominant then. In the fall monthly-mean causal-effect network (Fig. 7c), we observe a cluster of specific humidity, downward longwave surface radiation and air temperature. Specific humidity and downward longwave radiation are physically closely related as moist air reflects more longwave radiation than dry air, giving raise to the water–vapor feedback. Additionally, the autocorrelation of sea-ice cover and air temperature rises during refreeze. Since the strongest connection between air temperature and sea-ice cover is contemporary and the autocorrelation reaches back one time step, this can be interpreted as part of the strong interdependence of these variables. Sea-ice cover (increase) is also driving both latent surface heat flux and sea-ice export during refreeze in contrast to ice melt, when sensible heat fluxes and sea-ice export are driving sea ice (loss). The transition from mainly circulation to mainly feedback driven sea-ice cover is in agreement with Barton and Veron (2012), who found moisture-related feedback mechanisms to be a strong driver of Laptev Sea sea ice in September and October and in agreement with Arctic-wide findings that sea ice has a stronger impact on its surrounding in fall than in spring (Rigor et al., 2002; Serreze et al., 2009; Screen et al., 2013; Parmentier et al., 2015)
4.2 Influence of sea-ice cover on the atmosphere over land
What does this mean for the connection from sea-ice cover to the atmosphere over land? In summer, we observe a contemporary connection between sea-ice cover and specific humidity over land (Fig. 5). This link is also present in the ice-melt and refreeze monthly-mean causal-effect networks, where we consider the atmosphere over land and ocean, in contrast to the summer, where we focus on the atmosphere over land. To better understand this connection, we inspect the composites. For high-sea-ice-cover years, we observe more days with on average southward wind – this is in compliance with our earlier findings: ice export is minimized, and wind tends to carry colder air into the Laptev Sea region. More southward wind also means more southward transport of heat and moisture (Fig. 6 and Table 1). When open water areas are reduced compared to the mean state, both humidity and air temperature over the ocean are also reduced. So, the increase in gross southward transport of both moisture and heat can be attributed to the increase in air masses flowing south. At the same time, gross northward transport decreases. For low-sea-ice-cover years, there are fewer days with southward wind than usual. Still, we have a higher gross southward transport of heat and moisture because the air is warmer and moister over open water. More heat and moisture gets transported per air volume. In contrast to high-sea-ice-cover years, we also have an increase in transport northward, corresponding to our earlier findings: low-sea-ice-cover years are caused by above-average northward wind and an influx of warm and moist air from the south. If we compute net moisture and heat fluxes for the low-sea-ice-cover years, the median is only marginally shifted towards an increase in southward transport for both variables. It is remarkable that the increase in southward transport more than balances the increase in northward transport we already identified above as the primary cause of above-average sea-ice cover decrease. We conclude that the contemporary links in the monthly-mean summer causal-effect network incorporate causality in both directions: a decrease in sea-ice cover, as in the ice-melt period, is most likely initiated by increased meridional wind pointing north. But the increases in temperature and moisture over land, which then accompany low sea ice, are both the cause, when transported northward, and result of the diminished ice cover due to consequently enhanced southward transport of heat and moisture in air. This is in line with a model study by Lawrence et al. (2008), who found that rapid sea-ice loss leads to warming on land. In contrast to their analysis, we break down this connection from sea ice to land and identify the responsible atmospheric variables. On the spatial (Laptev Sea) and temporal (days to months) scales analyzed here, we do not find a connection via large-scale circulation starting with changes in sea ice. However, we do not rule out that changes in whole-Arctic sea-ice cover or in other regions, such as the Kara and Barents Sea, can have an impact on land through large-scale circulation; such connections between sea ice and atmospheric circulation patterns have been shown before (e.g., Ogi et al., 2016; Yao et al., 2017; Luo et al., 2019a; Li et al., 2020). Concerning local and direct links from sea ice to land in the Laptev Sea region, we find that over land the impact of sea ice is smaller than the impact of large-scale drivers, which simultaneously lead to changes in local sea-ice cover.
Concerning our research question on the drivers of the sea-ice variability, we conclude the following from our above model analysis.
Drivers of sea-ice-cover variability differ during the early melt (April–June) and refreeze (September–November) period. During early melt, meridional wind and connected directional variables, like ice export, are the prime drivers, all of them determined by the large-scale circulation. We find that high sea-ice cover often coincides with a negative Arctic Oscillation pattern, while low sea-ice cover coincides with an Arctic Dipole-like pressure pattern. During refreeze, in line with prior research (Rigor et al., 2002; Serreze et al., 2009; Screen et al., 2013), we conclude that thermodynamic feedback mechanisms related to temperature and moisture determine the speed of ice formation.
Regarding the connection between sea ice and the state of the atmosphere over land, we find the following.
Sea-ice cover is primarily connected to the atmosphere directly above it, and signals wash out quickly as we move to the land. Nevertheless, during low-sea-ice-cover years southward heat and especially moisture transport is enhanced. Heat and moisture are both variables which have a strong influence on the carbon budget as they affect both photosynthesis and respiration.
We show that the impact of sea ice on moisture is at least as strong as on temperature. Studies which only focus on temperature overlook an important pathway of information, especially since the coupling is weak.
A possible explanation for the lack of causal links between the atmosphere over land and sea-ice cover lies within temporal scales: the daily variations in the signal emitted by sea-ice cover have a significantly smaller impact than other forcings on the atmosphere over land. Otherwise, we would detect links on a daily scale. A potential slow but persistent forcing from sea-ice cover on the other hand might be diluted on a monthly scale due to the fast interactions within the atmosphere especially because sea-ice cover exerts a remote forcing. We might still detect a link for the slow forcing but more likely several atmospheric variables and a contemporary instead of a causal link.
The above analysis was done using the climate of the last century before the onset of strong changes in the Arctic. The observed loss of sea ice in the Arctic as a whole has been projected to lead to changes in the general circulation pattern, namely via a shift to more negative Arctic Oscillation and North Atlantic Oscillation indices, which favor high sea-ice cover and thus exert a negative feedback (Jaiser et al., 2012). No clear connection to anthropogenic forcing has yet been drawn for the Arctic Dipole, but from reanalysis we know that there is a higher persistence of this pattern in the 21st century compared to the last (Overland and Wang, 2010). Depending on which of these two shifts persists over the other we can expect a faster or slower decrease in sea ice in the Laptev Sea.
Since sea ice is generally decreasing, we expect that southward heat and moisture transport from the open ocean to the land will increase. Higher temperatures lead to deeper thawing of the seasonally frozen soil layers during summer due to enhanced southward transport of heat and moisture in air. This increases soil-organic-matter decomposition in the then-unfrozen soils, enhancing greenhouse gas emissions, but temperature is also one of the parameters steering the rate of photosynthesis, to name two of many processes which will be altered. Higher air moisture improves photosynthetic rates (high moisture leads to more efficient carbon fixation); changes the amount of (both short- and longwave) incoming solar radiation due to stronger reflection; or, if higher humidity in the air leads to higher precipitation and water tables, changes the pathway of soil-organic-matter decomposition (Shaver et al., 2000). A general warming and an enhanced hydrological cycle are key features of global climate change (Huntington, 2006; Stocker et al., 2013). In our model study we find that lower-than-usual sea ice in the Laptev Sea causes warming and an increase in air moisture over land, which might add to the abovementioned trends. Nevertheless, we found the link from sea ice to land to be weak under stable conditions, and, if this relation holds under different conditions, we expect climate change over land to be driven primarily by large-scale circulation. Though a general decrease in sea ice in the whole Arctic has an effect on large-scale circulation (Jaiser et al., 2012; Li and Wang, 2013), which then affects the Arctic landmasses, the local and direct link from sea-ice cover to adjacent landmasses appears marginal.
a Sea-ice export and import are computed by summing the gross positive and negative values of transects at the outer borders of the areas indicated by the masks in Fig. 3. b To estimate the influence of the variable on the atmosphere over land, we calculated the flow through the southern edge of the area of interest (as defined in Fig. 3), instead of using the mean over the whole area, as we do for most other variables. c Vertically integrated over all atmospheric model layers.
Primary data and scripts used in this study are archived by the Max Planck Institute for Meteorology and can be obtained by contacting email@example.com. The model output is available through Niederdrenk and Mikolajewicz (2014) at https://cera-www.dkrz.de/WDCC/ui/cerasearch/entry?acronym=DKRZ_LTA_899_ds00003 (last access: 18 November 2020).
All authors designed the study together, ZR conducted the data analysis and prepared the manuscript, and ALN provided the model simulations and edited the manuscript.
The authors declare that they have no conflict of interest.
We thank Thomas Kleinen and two anonymous reviewers for their helpful comments on the manuscript. We thank Marlene Kretschmer for her introduction to causal-effect networks.
This work was partially funded by the DFG Cluster of Excellence CliSAP (grant no. EXC177).
The article processing charges for this open-access publication were covered by the Max Planck Society.
This paper was edited by David Schroeder and reviewed by two anonymous referees.
Bareiss, J. and Görgen, K.: Spatial and temporal variability of sea ice in the Laptev Sea: Analyses and review of satellite passive-microwave data and model results, 1979 to 2002, Global Planet. Change, 48, 28–54, https://doi.org/10.1016/j.gloplacha.2004.12.004, 2005. a
Barton, N. P. and Veron, D. E.: Response of clouds and surface energy fluxes to changes in sea-ice cover over the Laptev Sea (Arctic Ocean), Clim. Res., 54, 69–84, https://doi.org/10.3354/cr01101, 2012. a
Bauer, M., Schröder, D., Heinemann, G., Willmes, S., and Ebner, L.: Quantifying polynya ice production in the Laptev Sea with the COSMO model, Polar Res., 32, 20922, https://doi.org/10.3402/polar.v32i0.20922, 2013. a
Bhatt, U. S., Alexander, M. A., Deser, C., Walsh, J. E., Miller, J. S., Timlin, M. S., Scott, J., and Tomas, R. A.: The Atmospheric Response to Realistic Reduced Summer Arctic Sea Ice Anomalies, Arctic Sea Ice Decline: Observations, Projections, Mechanisms, and Implications, 180, 91–110, https://doi.org/10.1029/180gm08, 2008. a
Deser, C., Walsh, J. E., and Timlin, M. S.: Arctic sea ice variability in the context of recent atmospheric circulation trends, J. Climate, 13, 617–633, https://doi.org/10.1175/1520-0442(2000)013<0617:Asivit>2.0.Co;2, 2000. a, b
Granger, C. W. J.: Investigating causal relations by econometric models and cross-spectral methods, Econometrica: journal of the Econometric Society, 37, 424–438, https://doi.org/10.2307/1912791, 1969. a
Graversen, R. G., Langen, P. L., and Mauritsen, T.: Polar Amplification in CCSM4: Contributions from the Lapse Rate and Surface Albedo Feedbacks, J. Climate, 27, 4433–4450, https://doi.org/10.1175/Jcli-D-13-00551.1, 2014. a
Haas, C. and Eicken, H.: Interannual variability of summer sea ice thickness in the Siberian and central Arctic under different atmospheric circulation regimes, J. Geophys. Res.-Oceans, 106, 4449–4462, https://doi.org/10.1029/1999jc000088, 2001. a, b
Jacob, D.: A note to the simulation of the annual and inter-annual variability of the water budget over the Baltic Sea drainage basin, Meteorol. Atmos. Phys., 77, 61–73, https://doi.org/10.1007/s007030170017, 2001. a
Jaiser, R., Dethloff, K., Handorf, D., Rinke, A., and Cohen, J.: Impact of sea ice cover changes on the Northern Hemisphere atmospheric winter circulation, Tellus Series A, 64, 11595, https://doi.org/10.3402/tellusa.v64i0.11595, 2012. a, b, c
Jungclaus, J. H., Fischer, N., Haak, H., Lohmann, K., Marotzke, J., Matei, D., Mikolajewicz, U., Notz, D., and von Storch, J. S.: Characteristics of the ocean simulations in the Max Planck Institute Ocean Model (MPIOM) the ocean component of the MPI-Earth system model, J. Adv. Model. Earth Syst., 5, 422–446, https://doi.org/10.1002/jame.20023, 2013. a
Kim, K.-Y., Hamlington, B. D., Na, H., and Kim, J.: Mechanism of seasonal Arctic sea ice evolution and Arctic amplification, The Cryosphere, 10, 2191–2202, https://doi.org/10.5194/tc-10-2191-2016, 2016. a
Kretschmer, M., Coumou, D., Donges, J. F., and Runge, J.: Using Causal Effect Networks to Analyze Different Arctic Drivers of Midlatitude Winter Circulation, J. Climate, 29, 4069–4081, https://doi.org/10.1175/Jcli-D-15-0654.1, 2016. a
Krumpen, T., Janout, M., Hodges, K. I., Gerdes, R., Girard-Ardhuin, F., Hölemann, J. A., and Willmes, S.: Variability and trends in Laptev Sea ice outflow between 1992–2011, The Cryosphere, 7, 349–363, https://doi.org/10.5194/tc-7-349-2013, 2013. a
Lawrence, D. M., Slater, A. G., Tomas, R. A., Holland, M. M., and Deser, C.: Accelerated Arctic land warming and permafrost degradation during rapid sea ice loss, Geophys. Res. Lett., 35, L11506, https://doi.org/10.1029/2008gl033985, 2008. a, b
Lee, S., Gong, T., Feldstein, S. B., Screen, J. A., and Simmonds, I.: Revisiting the Cause of the 1989–2009 Arctic Surface Warming Using the Surface Energy Budget: Downward Infrared Radiation Dominates the Surface Fluxes, Geophys. Res. Lett., 44, 10654–10661, https://doi.org/10.1002/2017GL075375, 2017. a
Li, F. and Wang, H. J.: Autumn Sea Ice Cover, Winter Northern Hemisphere Annular Mode, and Winter Precipitation in Eurasia, J. Climate, 26, 3968–3981, https://doi.org/10.1175/Jcli-D-12-00380.1, 2013. a
Li, M., Luo, D., Simmonds, I., Dai, A., Zhong, L., and Yao, Y.: Anchoring of atmospheric teleconnection patterns by Arctic sea ice loss and its link to winter cold anomalies in East Asia, Int. J. Climatol., 1–12, https://doi.org/10.1002/joc.6637, 2020. a
Luo, B., Luo, D., Wu, L., Zhong, L., and Simmonds, I.: Atmospheric circulation patterns which promote winter Arctic sea ice decline, Environ. Res. Lett., 12, 054017, https://doi.org/10.1088/1748-9326/aa69d0, 2017. a
Luo, B., Wu, L., Luo, D., Dai, A., and Simmonds, I.: The winter midlatitude-Arctic interaction: effects of North Atlantic SST and high-latitude blocking on Arctic sea ice and Eurasian cooling, Clim. Dynam., 52, 2981–3004, https://doi.org/10.1007/s00382-018-4301-5, 2019a. a
Luo, D., Chen, X., Overland, J., Simmonds, I., Wu, Y., and Zhang, P.: Weakened Potential Vorticity Barrier Linked to Recent Winter Arctic Sea Ice Loss and Midlatitude Cold Extremes, J. Climate, 32, 4235–4261, https://doi.org/10.1175/JCLI-D-18-0449.1, 2019b. a
Macias-Fauria, M., Karlsen, S. R., and Forbes, B. C.: Disentangling the coupling between sea ice and tundra productivity in Svalbard, Sci. Rep., 7, 1–10, https://doi.org/10.1038/s41598-017-06218-8, 2017. a
Niederdrenk, A. L., Sein, D. V., and Mikolajewicz, U.: Interannual variability of the Arctic freshwater cycle in the second half of the twentieth century in a regionally coupled climate model, Clim. Dynam., 47, 3883–3900, https://doi.org/10.1007/s00382-016-3047-1, 2016. a, b
Niederdrenk, L. and Mikolajewicz, U.: Wechselwirkungen zwischen verschiedenen Komponenten des arktischen Klimasystems (bm0899), available at: http://cera-www.dkrz.de/WDCC/ui/Compact.jsp?acronym=DKRZ_lta_899 (last access: 18 November 2020), 2014. a
Ogi, M., Rysgaard, S., and Barber, D. G.: The influence of winter and summer atmospheric circulation on the variability of temperature and sea ice around Greenland, Tellus Series A, 68, 14, https://doi.org/10.3402/tellusa.v68.31971, 2016. a, b
Olonscheck, D., Mauritsen, T., and Notz, D.: Arctic sea-ice variability is primarily driven by atmospheric temperature fluctuations, Nat. Geosci., 12, 430–434, https://doi.org/10.1038/s41561-019-0363-1, 2019. a
Overland, J. E. and Wang, M. Y.: Large-scale atmospheric circulation changes are associated with the recent loss of Arctic sea ice, Tellus Series A, 62, 1–9, https://doi.org/10.1111/j.1600-0870.2009.00421.x, 2010. a
Overland, J. E., Francis, J. A., Hanna, E., and Wang, M. Y.: The recent shift in early summer Arctic atmospheric circulation, Geophys. Res. Lett., 39, L19804, https://doi.org/10.1029/2012gl053268, 2012. a, b, c
Parmentier, F. J. W., Christensen, T. R., Sorensen, L. L., Rysgaard, S., McGuire, A. D., Miller, P. A., and Walker, D. A.: The impact of lower sea-ice extent on Arctic greenhouse-gas exchange, Nat. Clim. Change, 3, 195–202, https://doi.org/10.1038/Nclimate1784, 2013. a
Parmentier, F. J. W., Zhang, W. X., Mi, Y. J., Zhu, X. D., van Huissteden, J., Hayes, D. J., Zhuang, Q. L., Christensen, T. R., and McGuire, A. D.: Rising methane emissions from northern wetlands associated with sea ice decline, Geophys. Res. Lett., 42, 7214–7222, https://doi.org/10.1002/2015GL065013, 2015. a, b
Pearl, J.: Causality: Models, reasoning and inference, Cambridge University Press, Cambridge, United Kingdom, 2000. a
Pithan, F., Medeiros, B., and Mauritsen, T.: Mixed-phase clouds cause climate model biases in Arctic wintertime temperature inversions, Clim. Dynam., 43, 289–303, https://doi.org/10.1007/s00382-013-1964-9, 2013. a
Post, E., Bhatt, U. S., Bitz, C. M., Brodie, J. F., Fulton, T. L., Hebblewhite, M., Kerby, J., Kutz, S. J., Stirling, I., and Walker, D. A.: Ecological consequences of sea-ice decline, Science, 341, 519–524, https://doi.org/10.1126/science.1235225, 2013. a
Rigor, I. G., Wallace, J. M., and Colony, R. L.: Response of sea ice to the Arctic oscillation, J. Climate, 15, 2648–2663, https://doi.org/10.1175/1520-0442(2002)015<2648:Rositt>2.0.Co;2, 2002. a, b
Runge, J., Heitzig, J., Marwan, N., and Kurths, J.: Quantifying causal coupling strength: A lag-specific measure for multivariate time series related to transfer entropy, Phys. Rev. E, 86, 15, https://doi.org/10.1103/PhysRevE.86.061121, 2012. a, b, c
Runge, J., Petoukhov, V., and Kurths, J.: Quantifying the Strength and Delay of Climatic Interactions: The Ambiguities of Cross Correlation and a Novel Measure Based on Graphical Models, J. Climate, 27, 720–739, https://doi.org/10.1175/Jcli-D-13-00159.1, 2014. a, b, c
Runge, J., Petoukhov, V., Donges, J. F., Hlinka, J., Jajcay, N., Vejmelka, M., Hartman, D., Marwan, N., Palus, M., and Kurths, J.: Identifying causal gateways and mediators in complex spatio-temporal systems, Nat. Commun., 6, 8502, https://doi.org/10.1038/ncomms9502, 2015. a, b
Runge, J., Nowack, P., Kretschmer, M., Flaxman, S., and Sejdinovic, D.: Detecting and quantifying causal associations in large nonlinear time series datasets, Sci. Adv., 5, eaau4996, https://doi.org/10.1126/sciadv.aau4996, 2019. a
Samarasinghe, S. M., McGraw, M. C., Barnes, E. A., and Ebert-Uphoff, I.: A study of links between the Arctic and the midlatitude jet stream using Granger and Pearl causality, Environmetrics, 30, e2540, https://doi.org/10.1002/env.2540, 2019. a
Screen, J. A., Simmonds, I., and Keay, K.: Dramatic interannual changes of perennial Arctic sea ice linked to abnormal summer storm activity, J. Geophys. Res., 116, D15105, https://doi.org/10.1029/2011JD015847, 2011. a
Screen, J. A., Simmonds, I., Deser, C., and Tomas, R.: The atmospheric response to three decades of observed Arctic sea ice loss, J. Climate, 26, 1230–1248, https://doi.org/10.1175/Jcli-D-12-00063.1, 2013. a, b, c, d, e
Screen, J. A., Bracegirdle, T. J., and Simmonds, I.: Polar climate change as manifest in atmospheric circulation, Curr. Clim. Change Rep., 4, 383–395, https://doi.org/10.1007/s40641-018-0111-4, 2018. a
Sein, D. V., Mikolajewicz, U., Groger, M., Fast, I., Cabos, W., Pinto, J. G., Hagemann, S., Semmler, T., Izquierdo, A., and Jacob, D.: Regionally coupled atmosphere-ocean-sea ice-marine biogeochemistry model ROM: 1. Description and validation, J. Adv. Model. Earth Syst., 7, 268–304, https://doi.org/10.1002/2014ms000357, 2015. a
Serreze, M. C., Barrett, A. P., Stroeve, J. C., Kindig, D. N., and Holland, M. M.: The emergence of surface-based Arctic amplification, The Cryosphere, 3, 11–19, https://doi.org/10.5194/tc-3-11-2009, 2009. a, b, c, d
Shaver, G. R., Canadell, J., Chapin, F. S., Gurevitch, J., Harte, J., Henry, G., Ineson, P., Jonasson, S., Melillo, J., Pitelka, L., and Rustad, L.: Global warming and terrestrial ecosystems: A conceptual framework for analysis, Bioscience, 50, 871–882, https://doi.org/10.1641/0006-3568(2000)050[0871:Gwatea]2.0.Co;2, 2000. a
Spirtes, P., Glymour, C. N., Scheines, R., and Heckerman, D.: Causation, prediction, and search, MIT press, Cambridge, Massachusetts, United States, 2000. a
Stocker, T. F., Qin, D., Plattner, G. K., Alexander, L. V., Allen, S. K., Bindoff, N. L., Bréon, F. M., Church, J. A., Cubasch, U., Emori, S., Forster, P., Friedlingstein, P., Gillett, N., Gregory, J. M., Hartmann, D. L., Jansen, E., Kirtman, B., Knutti, R., Kumar, K. K., Lemke, P., Marotzke, J., Masson-Delmotte, V., Meehl, G. A., Mokhov, I. I., Piao, S., Ramaswamy, V., Randall, D., Rhein, M., Rojas, M., Sabine, C., Shindell, D., Talley, L. D., Vaughan, D. G., and Xie, S. P.: Technical Summary, pp. 33–115, Cambridge University Press, Cambridge, United Kingdom, 2013. a
Titchner, H. A. and Rayner, N. A.: The Met Office Hadley Centre sea ice and sea surface temperature data set, version 2:1. Sea ice concentrations, J. Geophys. Res.-Atmos., 119, 2864–2889, https://doi.org/10.1002/2013JD020316, 2014. a
Wang, J., Zhang, J. L., Watanabe, E., Ikeda, M., Mizobata, K., Walsh, J. E., Bai, X. Z., and Wu, B. Y.: Is the Dipole Anomaly a major driver to record lows in Arctic summer sea ice extent?, Geophys. Res. Lett., 36, 5, https://doi.org/10.1029/2008GL036706, 2009. a, b, c, d
Yang, W. C. and Magnusdottir, G.: Springtime extreme moisture transport into the Arctic and its impact on sea ice concentration, J. Geophys. Res.-Atmos., 122, 5316–5329, https://doi.org/10.1002/2016JD026324, 2017. a, b
Yao, Y., Luo, D., Dai, A., and Simmonds, I.: Increased quasi stationarity and persistence of winter Ural blocking and Eurasian extreme cold events in response to Arctic warming. Part I: Insights from observational analyses, J. Climate, 30, 3549–3568, https://doi.org/10.1175/Jcli-D-16-0261.1, 2017. a