How might the North American ice sheet influence the northwestern Eurasian climate?

. It is now widely acknowledged that past Northern Hemisphere ice sheets covering Canada and northern Europe at the Last Glacial Maximum (LGM) exerted a strong inﬂuence on climate by causing changes in atmospheric and oceanic circulations. In turn, these changes may have impacted the development of the ice sheets themselves through a combination of different feedback mechanisms. The present study is designed to investigate the potential impact of the North American ice sheet on the surface mass balance (SMB) of the Eurasian ice sheet driven by simulated changes in the past glacial atmospheric circulation. Using the LMDZ5 atmospheric circulation model, we carried out 12 experiments under constant LGM conditions for insolation, greenhouse gases and ocean. In these experiments, the Eurasian ice sheet is removed. The 12 experiments differ in the North American ice-sheet topography, ranging from a white and ﬂat (present-day topography) ice sheet to a full-size LGM ice sheet. This experimental design allows the albedo and the topographic impacts of the North American ice sheet onto


Abstract.
It is now widely acknowledged that past Northern Hemisphere ice sheets covering Canada and northern Europe at the Last Glacial Maximum (LGM) exerted a strong influence on climate by causing changes in atmospheric and oceanic circulations.In turn, these changes may have impacted the development of the ice sheets themselves through a combination of different feedback mechanisms.The present study is designed to investigate the potential impact of the North American ice sheet on the surface mass balance (SMB) of the Eurasian ice sheet driven by simulated changes in the past glacial atmospheric circulation.Using the LMDZ5 atmospheric circulation model, we carried out 12 experiments under constant LGM conditions for insolation, greenhouse gases and ocean.In these experiments, the Eurasian ice sheet is removed.The 12 experiments differ in the North American ice-sheet topography, ranging from a white and flat (present-day topography) ice sheet to a fullsize LGM ice sheet.This experimental design allows the albedo and the topographic impacts of the North American ice sheet onto the climate to be disentangled.The results are compared to our baseline experiment where both the North American and the Eurasian ice sheets have been removed.In summer, the sole albedo effect of the American ice sheet modifies the pattern of planetary waves with respect to the no-ice-sheet case, resulting in a cooling of the northwestern Eurasian region.By contrast, the atmospheric circulation changes induced by the topography of the North American ice sheet lead to a strong decrease of this cooling.In winter, the Scandinavian and the Barents-Kara regions respond differently to the American ice-sheet albedo effect: in response to atmospheric circulation changes, Scandinavia becomes warmer and total precipitation is more abundant, whereas the Barents-Kara area becomes cooler with a decrease of convective processes, causing a decrease of total precipitation.The gradual increase of the altitude of the American ice sheet leads to less total precipitation and snowfall and to colder temperatures over both the Scandinavian and the Barents and Kara sea sectors.We then compute the resulting annual surface mass balance over the Fennoscandian region from the simulated temperature and precipitation fields used to force an ice-sheet model.It clearly appears that the SMB is dominated by the ablation signal.In response to the summer cooling induced by the American ice-sheet albedo, high positive SMB values are obtained over the Eurasian region, leading thus to the growth of an ice sheet.On the contrary, the gradual increase of the American ice-sheet altitude induces more ablation over the Eurasian sector, hence limiting the growth of Fennoscandia.To test the robustness of our results with respect to the Eurasian ice sheet state, we carried out two additional LMDZ experiments with new boundary conditions involving both the American (flat or full LGM) and high Eurasian ice sheets.The most striking result is that the Eurasian ice sheet is maintained under full-LGM North American ice-sheet conditions, but loses ∼ 10 % of its mass compared to the case in which the North American ice sheet is flat.These new findings qualitatively confirm the conclusions from our first series of experiments and suggest that the development of the Eurasian ice sheet may have been slowed down by the growth of the American ice sheet, offering thereby a new understanding of the evolution of Northern Hemisphere ice sheets throughout glacial-interglacial cycles.

Introduction
The last million years is characterized by alternating glacial and interglacial phases.During glacial periods, large ice sheets covering present-day Canada and northwestern Eurasia (Dyke and Prest, 1987;Clark et al., 1993;Peltier, 2004;Svendsen et al., 2004;Lambeck et al., 2006;Tarasov et al., 2012) exerted a strong influence on climate.Several studies highlighted the importance the climatic changes induced by continental-scale ice sheets (e.g.Clark, 1999).It has been recognized that the ice-sheet topography is likely to be the main factor altering the atmospheric circulation in the Northern Hemisphere (Broccoli and Manabe, 1987;Pausata et al., 2011).As an example, the pioneering study carried out by Manabe and Broccoli (1985) with an atmospheric general circulation model shows that the North American ice sheet caused a split of the westerlies.The authors also highlight a larger amplitude of the planetary waves due to the presence of ice sheets.Recently, Ullman et al. (2014) tested the influence of different North American ice-sheet reconstructions on the climate.They showed that a higher American ice sheet leads to a more zonal Atlantic jet, thereby confirming the key role of ice-sheet topography on atmospheric circulation.Using model outputs from the Paleo Modelling Intercomparison Project 2 (PMIP2) (Braconnot et al., 2007), Laîné et al. (2009) and Rivière et al. (2010) highlight a strengthening and an equatorward displacement of the sub-tropical jet-stream during the Last Glacial Maximum (LGM) with respect to the pre-industrial period.Changes in the position and the strength of the North Atlantic jet stream induce changes in the storm tracks, and therefore in precipitation (Hall et al., 1996;Kageyama and Valdes, 2000;Laîné et al., 2009;Rivière et al., 2010).These changes also have an influence on the energy transport and can therefore modify the temperature.
The climatic changes induced by large-scale ice sheets exert an influence on both temperature and precipitation that drive the ice-sheet surface mass balance (SMB), defined as the sum of snow accumulation and snow and ice ablation.Using a simple ice-sheet model based on an idealized geometry coupled to a stationary-wave model, Roe and Lindzen (2001a, b) underline the importance of accounting for the feedbacks between ice sheets and temperatures induced by changes in stationary waves to properly simulate the evolution of an ice sheet.They show that the self-induced temperature anomaly due to an ice sheet leads to a warming over the ice-sheet western part.This may explain the absence of ice over Alaska at the LGM.They also suggest that the stationary waves excited by the North American ice sheet may have contributed to a warming over Europe, influencing the development of the Eurasian ice sheet.In the same way, with a three-dimensional stationary wave model, Liakka et al. (2011) showed that the position of the southern margin of ice sheets strongly depends on the temperature anomalies due to stationary waves, which are modified by the ice sheet itself.More recently, with the use of the CAM3 atmospheric model run under four different climatic contexts (last interglacial, marine isotope stages 5b and 4 (MIS5b and MIS4), and LGM periods), Löfverström et al. (2014) showed how the atmospheric circulation changes induced by the ice sheets could have influenced the growth of the ice sheets themselves.Similarly to Roe and Lindzen (2001a), they show that the summer atmospheric circulation change due to the presence of ice sheets may cause a sufficient warming over Siberia and Alaska to inhibit ice growth.They also perform two experiments under MIS4 conditions to test the influence of one ice sheet on the other.They conclude that the summer temperature anomaly induced by the presence of the MIS4 American ice sheet is too weak to explain the small size of the European ice sheet at that period, but may have contributed to the westward shift of the ice-sheet mass center.Although several studies have been devoted to the mutual influences of the ice sheets (Roe and Lindzen, 2001a;Beghin et al., 2014;Löfverström et al., 2014), no study has specifically investigated the mechanisms through which the American ice sheet may have influenced the European climate, and therefore the European ice-sheet surface mass balance.
The aim of this paper is to investigate the atmospheric response to changes in the North American ice-sheet boundary conditions and to examine how the related atmospheric circulation changes impact the climate over the Scandinavian and Barents-Kara sea regions.To achieve this goal, we use the LMDZ5 atmospheric general circulation model run with different thicknesses of the North American ice sheet taken as boundary conditions, with no ice at all over Eurasia.We investigate the mechanisms by which the American ice sheet may change the surface mass balance of the Eurasian ice sheet.Besides gradually increasing the altitude of the North American ice sheet, the originality of the present study also relies on ice-sheet model simulations.We use climatic fields simulated by the LMDZ5 model as inputs to a three-dimensional ice-sheet model to compute the surface mass balance of the Eurasian ice sheet.The description of the climate and ice-sheet models is given in Sect. 2 as well as the experimental design.The model results, presented in Sects.3 and 4, are followed by a discussion (Sect.5).Section 6 summarizes the main findings of our study.

The atmospheric model
The LMDZ5 atmospheric model used in this study belongs to the LMDZ model family developed at Laboratoire de Météorologie Dynamique (Sadourny and Laval, 1984;Le Treut et al., 1994, 1998, Hourdin et al., 2006).LMDZ is the atmospheric component of the IPSL-CM5A coupled oceanatmosphere model (Dufresne et al., 2013).The dynamical equations are discretized on a longitude-latitude-staggered Arakawa C-grid (Kasahara, 1977).The model ensures the conservation of both enstrophy (square of wind rotational) for baroclinic flows (Sadourny, 1975a, b) and the axisymmetric component of the angular momentum.The model version used in this study has 39 vertical levels and runs on a 96 × 95 model grid resolution (3.75 • ×1.9 • ).A complete description of the model can be found in (Hourdin et al., 2013).

The ice-sheet model
The GRISLI ice-sheet model is a three-dimensional thermomechanical model which simulates the evolution of icesheet geometry (extension and thickness) and the coupled temperature-velocity fields in response to climate forcing.A comprehensive description of the model can be found in (Ritz et al., 2001) and (Peyaud et al., 2007) tions are solved on a cartesian grid (40 km × 40 km).Over the grounded part of the ice sheet, the ice flow resulting from internal deformation is governed by the shallow-ice approximation (Hutter, 1983;Morland et al., 1984).The model deals with ice flow through ice shelves using the shallow-shelf approximation (MacAyeal, 1989).It also predicts the largescale characteristics of the ice streams using criteria based on the effective pressure and hydraulic load.At each time step, the velocity and vertical profiles of temperature in the ice are computed as well as the new geometry of the ice sheet.The isostatic adjustment of bedrock in response to ice load is governed by the flow of the asthenosphere, with a characteristic time constant of 3000 years, and by the rigidity of the lithosphere.The temperature field is computed both in the ice and in the bedrock by solving a time-dependent heat equation.The surface mass balance is defined as the sum between accumulation and ablation computed by the empirical positive degree-day (PDD) method (Reeh, 1991;Fausto et al., 2009).This method assumes that melt rates of snow and ice are linearly related to the number of PDD through degree-day factors for snow and ice materials (Braithwaite, 1984(Braithwaite, , 1995)).

Experimental set-up
In order to investigate the albedo and topography effects induced by the North American ice sheet on the Eurasian climate at the LGM, we carried out 12 simulations under LGM conditions for greenhouse gases (GHGs), insolation, sea-surface temperatures and sea-ice.In this series of simulations, the altitude of the North American ice-sheet surface ranges from that of the present-day surface to 100 % of that used in the PMIP3 LGM experiments.The PMIP3 LGM ice sheets result from a combination of three reconstructions, namely ICE-6G v2.0 (Toscano et al., 2011), GLAC-1 (Tarasov et al., 2012) and ANU (Lambeck, 2001).The way this new reconstruction has been obtained is explained in detail on the PMIP3 website (https://wiki.lsce.ipsl.fr/pmip3/doku.php/pmip3:design:pi:final:icesheet) and in Abe-Ouchi et al. (2015).In our baseline experiment (noIS), the land-ice mask is modified (with respect to PMIP3) to remove both the European and the American ice sheets and the altitude is that of the present day (Fig. 1).In the other simulations, we only remove the European ice-sheet mask.These simulations are referred to as xxdhL, where dhL represents the surface height difference between the PMIP3 LGM Laurentide ice sheet and the present-day surface, and xx represents the percentage of dhL taken into account.To simplify the writing, the American and Eurasian ice sheets are respectively referred to as "Laurentide ice sheet" (LIS) and the "Fennoscandian ice sheet" (FIS) in the following.
The topography of the 00dhL experiment is therefore the same as today, but the land-ice mask is set to the LGM one (Figs. 1 and 2).Greenland and Antarctic ice sheets are the same as in the PMIP3 experimental set-up.Both insolation and GHG boundary conditions are similar to those defined in the PMIP3 protocol: the orbital forcing is taken at 21 kyr BP from (Berger, 1978), while atmospheric GHG concentrations are those recorded in Antarctic ice cores (CO 2 = 185 ppm, (Lüthi et al., 2008); CH 4 = 350 ppb, (Loulergue et al., 2008); N 2 O = 200 ppb, (Spahni et al., 2005)).The LGM land-sea mask is also taken into account, with closed Bering Strait, and land in Hudson Bay and Barents Sea (Fig. 1).The seasurface temperatures (SSTs) and sea-ice cover come from 50 consecutive years of the IPSL PMIP3 LGM run outputs (Kageyama et al., 2013).These oceanic surface conditions are depicted in Fig. 1.In the North Atlantic area, the temperature meridional gradient is steep in the western part of the basin and weakens in its eastern part, as for the pre-industrial period.The sea-ice edge is located south of Iceland in winter and retreats to the northern limits of the Nordic seas in summer.These rather warm North Atlantic Ocean/Nordic sea conditions are due to a strong Atlantic meridional overturning circulation (Kageyama et al., 2013).Using these boundary conditions, the LMDZ5 simulations are run for 50 years, and we study the last 40 years to be at equilibrium.Two-sided Student t tests with a p value = 0.05 (von Storch and Zwiers, 2003) are performed on key diagnostics.

AGCM results
In this section, we focus on how the LIS influences the climatic fields which can have an impact on the surface mass balance of an ice sheet.The key climatic variables used to compute the amount of ablation with the PDD method are the monthly temperatures and the monthly total (solid plus liquid) precipitation.

Impact of the Laurentide ice sheet
Using the comparison between the 00dhL and the noIS experiments, we first consider the impact of the LIS albedo on the 2 m summer air temperature over the Scandinavia-Barents-Kara (SBK) region (Fig. 3, 00dhL experiment, left column).The most striking feature is a cooling over the whole Eurasian continent with a maximum (more than 10 • C) centered over the Barents and Kara seas and a moderate cooling (2 to 4 • C) over Scandinavia.The area experiencing the largest cooling also corresponds to the area with the largest albedo (orange contours).This is a direct consequence of the LGM land-sea mask used in this study for which the Barents and Kara seas are represented as continental areas.This favors snow accumulation and therefore higher albedo values in this area, causing in turn a positive feedback on temperatures.Figure 3 also displays the 2 m mean summer temperature over the SBK region for different altitudes of the LIS.Here, we only give results of selected experiments to simplify the presentation of our results.It appears that the SBK cooling gets smaller when the LIS gets higher.Overall, the simulated summer surface air temperature evolves together with the LIS altitude, with a progressive reduction of the zonal and meridional extent of the largest cooling zone.When the LIS reaches its full LGM size (100dhL-noIS), the summer cooling over the Barents and Kara seas is reduced to 4-6 and to 2 • C at most over Scandinavia.The topography of the LIS has therefore a warming effect on the SBK region for the summer season.In winter, the response of the temperature to the LIS albedo and its evolution with the LIS altitude are different from the summer response (Fig. 3, right panels).First, the LIS albedo induces warmer temperatures (+3 • C) over southern Scandinavia, the British Isles and the present-day North and Baltic seas, and cooler temperatures over an area in the Barents Sea expanding from Svalbard to No-vaya Zemlya (3 • C).Surprisingly, in the 30dhL experiment, the amplitude of these temperature changes (with respect to noIS) is reduced, showing that the response of winter temperatures to changes in LIS altitude is far from being linear.Beyond 60dhL, the cooling zone progressively expands, and spreads across the entire British Isles and Scandinavia-Barents-Kara area when the LIS has its full LGM size.In contrast to the climatic response simulated for the summer season, the LIS topography leads to colder surface air temperatures in winter.
To understand the origin of these contrasted responses, we investigate the processes through which the changes in the altitude of the Laurentide ice sheet modify the atmospheric circulation.

Atmospheric circulation processes
The comparison of the zonal anomaly (00dhL vs. noIS) of the 500 hPa summer geopotential height shows that the sole albedo effect of the LIS is sufficient to drastically change the atmospheric circulation (Fig. 4, 00dhL and noIS, left column).The ridge over the Rockies, clearly visible in noIS, disappears in 00dhL, and the trough over the Labrador sector is more extensive.Ridges over the North Atlantic, the Greenland ice sheet and northern Europe are more developed in 00dhL (with respect to noIS).By contrast, the troughs over Iceland, the Norwegian and Barents seas appearing in noIS are weaker or even vanish in 00dhL.In noIS, the trough over Svalbard implies southerlies over the Barents Sea (Fig. 4, right column) and therefore, warm temperatures in this area (Fig. 3).The weakening of this trough in 00dhL as well as the stronger ridge over Greenland induces northerlies over Svalbard and the Barents Sea (Fig. 4) and therefore colder temperatures (Fig. 3).
The Labrador trough and the North Atlantic ridge become stronger when the LIS gets higher, until it reaches 60 % of its full size.Beyond 60dhL, the Labrador trough and the Atlantic ridge keep more or less the same amplitude.When the LIS gets higher, the Greenland ridge gets weaker, the European ridge remains almost unchanged (Fig. 4), and there is a return of the Iceland-Svalbard-Barents trough.This new trough, centered between Greenland and Svalbard, along with the weakening of the Greenland ridge, brings southerlies again over the northern Barents Sea.The Barents winds shift from northerlies to southerlies when the LIS is 60 % of its full size.The southerlies bring warmer temperatures, explaining the SBK warming with the increase of the LIS height (Figs. 3 and 4).In winter (Fig. 5), the sole albedo effect of the LIS does not drastically change the geopotential anomaly.However, the meridional wind field shows a slight decrease of the southerlies over the Barents Sea, explaining the slight cooling over this region.There is also a slight decrease of the northerlies simulated over the British Isles, the eastern North Atlantic Ocean and southern Scandinavia, consistent with the moderate warming simulated over this region.In the 30dhL experiment, a strengthening of the southerlies is observed (Fig. 5) inducing a warming (with respect to 00dhL) over the SBK region (Fig. 3) and explaining the non-linear behavior on winter temperatures.When the altitude of the LIS further increases, there is a progressive emergence of a trough over the Scandinavian-Barents sector.This trough comes with northerlies over the Norwegian Sea and weaker southerlies over the Barents sea (Fig. 5, right column), consistently with colder temperatures over this area, as seen in Fig. 3.This analysis explains why opposite sign temperature responses are obtained for the summer and winter seasons.Since ablation is rather sensitive to the summer season, more ice is expected over the Fennoscandian area when the sole albedo effect of the LIS is operating (i.e.00dhL experiment).However, the development of an ice sheet in this area also depends on the amount of snow accumulation.Therefore, we next examine hereafter the LIS impact on the SBK total precipitation and snowfall.

Impact of the Laurentide ice sheet
In this section, we investigate the impact of the LIS on the mean annual precipitation and snowfall over northwestern Eurasia (Fig. 6).In response to the LIS albedo effect (00dhL-noIS panel) a slight excess of precipitation is simulated over the northeastern North Atlantic and the southern and west-ern parts of Scandinavia, but less precipitation is simulated over the Barents-Kara sea area.However, over the whole SBK region, the snowfall anomaly is positive.This increase in snowfall is due to the summer signal (not shown) related to the cooling observed over this region (Fig. 3).This cooling is strong enough to allow for an increase of the fraction of solid precipitation and to maintain the snow cover.
When the LIS reaches its full LGM altitude (100dhL-noIS), the positive anomaly of total precipitation over Scandinavia simulated in the 00dhL experiment shifts southward, reaching even the French and the Iberian Peninsula Atlantic coasts.This results in dry conditions over the western flank of the Scandinavian region, both in terms of total and solid precipitation.Concomitantly, the negative anomaly of total precipitation over the Barents-Kara sea region simulated in 00dhL expands westward leading to a drier climate and to a decrease of snowfall over Scandinavia.These results suggest that the Barents and Kara seas and the Scandinavian regions are sensitive to different characteristics of the LIS; while the precipitation over Scandinavia appears to be sensitive to the LIS height, the precipitation anomaly over the Barents and Kara seas mainly appears to result from the LIS albedo effect and is rather insensitive to height.

Large-scale vs. convective precipitation
To disentangle the mechanisms responsible for precipitation changes in these regions, we split up the total precipitation into its large-scale and convective components (Fig. 7).The large-scale component behaves similarly to the total precipitation shown in Fig. 6 with a positive anomaly in the 00dhL experiment simulated over Scandinavia and the eastern North Atlantic sector and the southward shift of this pattern simulated when the LIS gets higher.Hence the drier signal over the western flank of Scandinavia obtained with the full-LGM LIS (100dhL-noIS) can also be found in the large-scale component of the total precipitation.This strongly suggests that the Scandinavian precipitation, as expected, is driven by large-scale processes.
The negative anomaly of precipitation appearing in the Arctic Ocean north of the Barents-Kara area (Fig. 6) is also clearly visible in the large-scale component (Fig. 7); by contrast, the large negative precipitation anomaly south of the Novaya Zemlya clearly arises from convective processes.Over the Barents-Kara area, the origin of the precipitation signal is less clear but seems to be a combination of both convective and large-scale processes.The decrease of convective precipitation negative anomaly (with respect to noIS) over the entire SBK sector when the LIS is higher can be easily explained by the weakening of the summer cooling discussed in section 3.1.
To further investigate the mechanisms governing precipitation over Scandinavia and the Barents-Kara area, we ex- amine the seasonality of precipitation.We first compute the average precipitation above Scandinavia (55-70 • N : 0-20 • E) (Fig. 8a) and the Barents-Kara region (65-90 • N : 20-100 • E) (Fig. 8b), then we compare summer and winter precipitation anomalies with the annual precipitation anomaly.
Results displayed in Figure 8 show that over Scandinavia winter precipitation changes are clearly correlated with annual precipitation (r = 0.99).This means that the annual signal observed over Scandinavia (Figs. 6 and 7) is mainly due to winter precipitation.Conversely, the Barents-Kara sector is rather governed by summer precipitation (at least for LIS altitudes up to 60 % of the PMIP3 one), as shown by the positive correlation (r = 0.98) obtained for the summer season (Fig. 8b).

Precipitation changes over Scandinavia: implication of the jet stream shift
We therefore look for a mechanism prevailing in winter to explain the Scandinavian precipitation differences.Using the PMIP3 model outputs for the LGM and the PI time slices, Beghin et al. (2015) found a correlation between the southward shift of the North Atlantic jet stream from the preindustrial to the LGM periods and the winter precipitation changes over the Iberian Peninsula.Following the same approach, we examine whether a similar mechanism to the one found in (Beghin et al., 2015) may explain the precipitation changes over Scandinavia when the LIS characteristics (albedo and topography) are modified.Similarly to Chavaillaz et al. (2013), we use the 850 hPa zonal wind to define the jet stream.The position of the jet displayed in Fig. 9 corresponds to the position of the 850 hPa meridional wind maximum.The LIS albedo effect induces a slight northward shift (with respect to noIS) of the North Atlantic jet stream (00dhL-noIS panel).As the LIS becomes higher, the jet moves progressively southward and shifts southward of its noIS position between 30 and 60dhL.From 60dhL and beyond, the Scandinavian winter precipitation anomaly becomes more and more negative.This suggests a relationship between the shift of the North Atlantic jet stream and the precipitation over the Scandinavian region.
To confirm this assumption, we plot the shift of the jet ( ϕ) as a function of the precipitation difference between the noIS and the xxdhL experiments over Scandinavia (Fig. 10).ϕ is defined as the difference of latitudes between noIS and xxdhL experiments where the 850 hPa zonal wind is at its highest value.The latitude of this maximum has been found by computing the zonal mean of the zonal wind over the Atlantic basin (50-10 • W).To minimize the uncertainties due to model discretization, we used a parabolic interpolation around the grid point where the mean zonal wind is largest.The maximum of the corresponding parabola yields a better approximation of the exact latitude of the maximum jet wind.This method allows for a more accurate finding of the latitude of the jet maximum (Chavaillaz et al., 2013;Beghin et al., 2015).We find a good linear correlation between the shift of the North Atlantic jet and the Scandinavian precipitation anomaly during the winter season (Fig. 10).The quantification of the range of the precipitation anomalies and of the jet shift is obtained by bootstrapping: a sample of thousand values is obtained by calculating the average of 1000 randomly picked samples of 40-year duration in the noIS and the xxdhL original samples.Linear regressions of the 5th and 95th percentiles of the bootstrapping sample confirm the close link between the shift of the North Atlantic jet stream and the winter precipitation changes over Scandinavia due to the effect of LIS.
The 500 hPa geopotential height (z500) zonal anomaly over the East Laurentide sector (90-60 • W) does not change significantly when a flat LIS is imposed (Figs. 5,11).As the LIS gets higher, the Labrador trough expands southward: Fig. 11 shows that the zero isoline is located at ∼ 43 • N in the noIS and 00dhL experiments and progressively reaches 37 • N for a full-LGM LIS.Over the Atlantic sector (50-10 • W) the Icelandic low/Azores high pattern remain similar in all experiments up to 60dhL and only show a significant change in the full LIS experiment (Figs. 5,11).Indeed, in 100dhL, the Icelandic low's amplitude weakens but simultaneously expands southward (Fig. 11).Consistently with these changes in the z500 zonal anomaly, the 500 hPa west-erlies shift southward over the whole Atlantic sector (Fig. 11, bottom).Again, the largest southward shift for this variable is between the 60 and 100dhL experiments, showing the nonlinearity of the atmospheric circulation to the LIS height.

Precipitation changes over the Barents-Kara area
As previously mentioned in Sect.3.2.2, the mean annual precipitation changes over the Barents-Kara area are dominated by the summer signal when the LIS remains at a relatively low level (i.e. up to 60dhL).This result is confirmed by the comparison of the patterns of annual (Fig. 6), summer and winter (Fig. 9) precipitation anomalies.The amount of summer precipitation is at its lowest level when the sole albedo effect of the LIS is accounted for.It increases with the LIS   altitude in response to warmer temperatures (Fig. 3).Conversely, the summer snowfall amount decreases from 00 to 100dhL due to the progressive lowering of the fraction of solid vs. liquid precipitation.

Summary
To summarize, the LIS albedo effect induces colder summer temperatures over the entire northwestern Eurasian region (SBK area) with respect to the noIS experiment.Over southern Scandinavia, the winter temperatures (Fig. 3) are warmer due to weaker northerlies (Fig. 5).More precipitation is also simulated in winter in response to the poleward shift of the North Atlantic jet stream (Fig. 6).As the Barents-Kara sector is dominated by the summer signal, less precipitation is simulated in response to colder summer temperatures (with respect to noIS), while a larger snowfall amount is obtained.
When the LIS gets higher, the negative anomaly of summer temperatures over the Barents-Kara region becomes weaker due to changes in the atmospheric circulation.In addition, the higher the LIS, the smaller the amounts of annual precipitation and snowfall over the Scandinavian and the Barents-Kara areas, in response to the equatorward shift of the jet.
4 What is the influence of the LIS on the growth of the FIS?

In case of no initial Fennoscandian ice sheet
What are the impacts of the LIS on the surface mass balance of both Scandinavian and Barents-Kara areas?To investigate this question, we use the LMDZ5 monthly temperatures and precipitation fields of each experiment to force the ice-sheet model GRISLI.LMDZ5 climatic fields are downscaled over the GRISLI grid using a bilinear interpolation.Due to the difference of resolution between the atmospheric and the ice-sheet models, temperature is vertically corrected using a linear vertical gradient of 6 K km −1 .Precipitation is also vertically corrected using an exponential function of the temperature (Charbit et al., 2002(Charbit et al., , 2007)).Snowfall is recalculated using the downscaled precipitation and temperature because the more detailed topography of GRISLI allows snowfall when LMDZ5 provides only liquid precipitation.
Figure 12 displays the snow accumulation and the ablation fields computed by GRISLI for the five selected experiments.As expected from the LMDZ5 results of annual precipitation and snowfall, snow accumulation over Scandinavia is larger when the LIS is flat and decreases as the LIS is higher.A similar observation can be made over the Barents-Kara region.Consistently with the simulated summer temperatures (Fig. 3), when the LIS is flat, the ablation is weaker over the entire SBK region, and increases as the LIS gets higher.Nevertheless, even with the full-LGM LIS, the ablation remains weaker than in the noIS experiment.The resulting surface mass balance is shown in Fig. 13a.The similarities between the ablation (Fig. 12b) and the surface mass balance patterns indicate that the surface mass balance is dominated by the ablation.In the absence of LIS, the surface mass balance is positive over only a small part of the SBK area, namely over Svalbard.Under these conditions, no ice sheet grows (Fig. 13b).When the LIS is flat, the surface mass balance is positive over the Barents and Kara seas and over the northern part of Scandinavia, allowing the growth of ice.The resulting FIS is less extended than those provided by the ICE-6G (Toscano et al., 2011) and the ANU (Lambeck, 2001) reconstructions.This is likely due to the absence of the albedo feedback, since our approach is not based on a full coupling method and does not allow to account for the impact of ice-sheet changes onto the climate.As the LIS gets higher, the limit of the positive surface mass balance shifts westward and northward, excluding the Kara Sea and Scandinavia from the positive surface-mass balance area.When the LIS has its full LGM size, the surface mass balance is positive only over the Svalbard.As a consequence, the simulated FIS is smaller.

What if the Fennoscandian ice sheet is initially high?
The aim of this study is not to provide a fully realistic scenario of the evolution of the ice sheets throughout an entire glacial cycle.Rather, based on idealized cases, it is to inves-tigate the impact of the Laurentide ice sheet on the Eurasian ice sheet surface mass balance.In the first part of this study, we have shown that the growth of the LIS prevents the development of the FIS.However, back to a more realistic context in which both ice sheets co-evolved, the relevance of this conclusion may be questioned and the issue of the FIS inception under full LGM LIS conditions may appear somehow artificial.Thus, to make our previous conclusions more robust, we perform two additional LMDZ experiments with new icesheet boundary conditions.In both simulations, the FIS has the elevation computed by the ice-sheet model for the 00dhL run (i.e. the largest FIS elevation of all the GRISLI experiments), while the LIS is either "flat" (00dhLFIS) or has its full LGM altitude (100dhLFIS).Here, our objective is to test the self-consistency of the overall study, by checking that the growth of the FIS is slowed down when the LIS gets higher.We first examine how a high FIS modifies the northwestern Eurasian climate under 00 and 100dhL conditions (i.e. with flat or full LGM LIS ice sheet), in terms of the temperature and precipitation fields used in our modeling approach to compute the ice-sheet surface mass balance.We then use the climatic outputs to force the ice-sheet model.

Response of the atmosphere
Figure 14 displays the summer surface air temperature, which is the main driver of the ablation, for the 00dhL (flat LIS, no FIS), 00dhLFIS (flat LIS, high FIS), 100dhL (high LIS, noIS) and 100dhLFIS (high LIS, high FIS) experiments.The comparison between Fig. 14a and b and c and d illustrate the impact of the Fennoscandian ice sheet on the 2 m sum- mer temperature.As expected, the presence of Fennoscandia induces a large cooling over the ice sheet itself, due to the temperature-altitude relationship and to the albedo effect.In previous sections, we have shown that the cooling observed over the Barents-Kara region is weakened when the LIS gets higher due to a modification of the large-scale atmospheric circulation (Figs.3-4).A similar mechanism is observed when the Fennoscandian ice sheet is high (100dhLFIS experiment): the cooling over the ice sheet is reduced and the warming along the southern margin is amplified (Fig. 14d vs. 14b).These features result from an increase of warm air advection in the Barents-Kara sector and to a decrease of northerly heat flux in the south.The presence of the FIS also induces a split of the North Atlantic jet stream in two branches resulting in northward deflected heat fluxes along the western margin of the ice sheet and the Barents Sea sector (Figs.15b, d and 16b, d).The heat advection modulates the ice-sheet altitude and albedo effects but is not sufficient to cancel out the cooling observed over the ice sheet.
Figure 17 displays maps of snow accumulation differences (after downscaling on the GRISLI grid) between the four experiments computed by GRISLI after downscaling.The presence of the Fennoscandian ice sheet results in an increase of precipitation and snow accumulation over its western flank, along the Norwegian Sea coast (Fig. 17a, c).This is due to the barrier effect of the ice sheet and to the subsequent increase of the vertical wind.This effect is less pronounced when the LIS is high (Fig. 17c) as a result of the southward displacement of the North Atlantic jet stream (Figs. 6,11).A second striking feature is the increase of snow accumulation in the southern part of the Eurasian ice sheet with respect to the no FIS situation.This is mainly due to colder temperatures over the ice sheet, increasing the fraction of total precipitation which falls as snow.In case of a flat LIS, this accumulation zone is confined to the southern margin and drier conditions (with respect to the noFIS case) in the interior of the Eurasian ice sheet are observed, mainly located over the Kara Sea (Fig. 17a), and likely due to the strong barrier effect.By contrast, when the LIS is high (Fig. 17c), differences between FIS and no FIS experiments are larger and expand over more than the southern half of the FIS.This coincides with a larger temperature difference over the FIS area when the LIS is high (Figs.14c, d) compared to the low LIS situation (Fig. 14a, b).
As mentioned above (Sect.3.2.1 and Fig. 6), in the absence of FIS, the total precipitation is driven by large-scale processes over the Scandinavian region and by convection over the Barents-Kara sector.When the LIS is high, the southward shift of the winter jet stream (Fig. 9) and the increase of summer temperature lead to an overall decrease of snowfall over the Eurasian region.When the FIS is present, the situation is rather different and convective precipitation is negligible (not shown) due to the cold temperatures induced by the presence of the FIS.The differences between the 00 and 100dhLFIS experiments are mainly located over the western and southeastern parts of the ice sheet with respectively larger and smaller snowfall amounts simulated in the 00dhLFIS experiment (with respect to 100dhLFIS).The excess of precipitation (total and solid) in the flat LIS/high FIS case is a direct consequence of the increase of both zonal and meridional humidity transports in winter (not shown) favored by the enhanced strength of the 850 hPa horizontal winds (Fig. 18).

Consequences on the simulated Eurasian ice sheets
The climatic outputs simulated in the 00 and 100dhLFIS experiments have been used to compute the surface mass balance of the Eurasian ice sheet at the GRISLI resolution and to force the ice-sheet model with the downscaling procedure described in Sect.4.1.Starting from the initial state provided by the 00dhL GRISLI run with flat (00dhLFIS) or full LGM LIS (100dhLFIS), GRISLI has been integrated during 100 kyr.At the end of the simulation, the simulated ice volume and ice covered area simulated in the 00dhLFIS experiment are respectively 10 and 11 % greater than those obtained with the 100dhLFIS simulation.The differences in simulated ice thickness between both simulations are displayed in Fig. 19 (00-100dhLFIS).In most parts of the simulated ice covered areas, ice thickness differences are fully consistent with the snow accumulation signal (Fig. 17).In particular, large positive anomalies (> 500 m) are found in the southwestern sector, that are greater than 1000 m near the ice sheet margin.In the northern part, lower ice thicknesses obtained in the 00dhLFIS experiment correspond to the area of lower snow accumulation (with respect to 100dhLFIS), as ablation is very small due to summer temperatures largely below 0 • C. Surprisingly, another large positive difference is found in the eastern part which is characterized by a lower amount of snow accumulation in case of flat LIS.(with respect to full LIS).This pattern results from warmer temperatures simulated in the southern and eastern parts of the 100dhLFIS Eurasian ice sheet in response to the decrease of cold air advection (see Sect. 4.2.1).Warmer temperatures found in the south and east of the FIS when the LIS has its full LGM size favor the westward and the northward retreats of the eastern and southern margins of the ice sheet.Since our approach does not account for the feedback of the simulated ice sheets on the climate, the amplitude of the different responses are likely to be underestimated due to an underestimation of the albedo effect.However, this analysis confirms that the growth of the FIS is slowed down when the LIS is high and demonstrates the self-consistency of the overall study.

Glaciation scenario
Geomorphological reconstructions revealed that the Fennoscandian ice sheet eastern limit reached the Putorana Mountains in Siberia during the Early Weichselian (90-80 kyr BP) (Mangerud et al., 1998;Svendsen et al., 2004).According to Svendsen et al. (2004), this eastern limit moved back to the Kara sea around 60 kyr BP, along with the southward extension of the Scandinavian ice sheet (Andersen and Mangerud, 1989;Lundqvist, 1992;Mangerud et al., 1998;Houmark-Nielsen, 1999).At the LGM, the ice sheet finally reached the British Isles, and the Kara ice shelf did not cover the entire Kara sea anymore (Landvik et al., 1998;Ehlers et al., 2004;Svendsen et al., 2004).More recently, P. Beghin et al.: North American ice sheet influence on the Eurasian climate Kleman et al. (2013) combined geological evidences and numerical modeling to reconstruct the paleotopography of past Northern Hemisphere ice sheets during MIS5b (∼ 86 kyr), and MIS4 (∼ 64 kyr) and 2 (LGM).They found that the LIS and the FIS grow in concert in the early part of the last glacial cycle (i.e. until MIS5b), while after MIS4, their evolution is radically different with a further growth of the LIS and a significant shrinkage of the FIS due to its southwestward retreat towards the main precipitation sources.Löfverström et al. (2014) propose to attribute the westward shift of the FIS and the decrease of the Kara ice sheet to warm temperatures over the southeastern part of the Fennoscandian region, induced by the FIS itself, along with the upslope precipitation effect proposed by Sanberg and Oerlemans (1983).In their study, the LIS exerts only a small influence on the retreat of the Kara ice sheet.While Löfverström et al. (2014) used realistic LIS topographies provided by the Kleman's et al. (2013) reconstructions, we used idealized geometries to study separately the impact of the LIS albedo and topography on climate.The results of our first series of xxdhL experiments clearly show that a large but low LIS favors a positive SMB over the SBK sector (Fig. 13, left panels) and therefore the FIS maintenance (Fig. 13, right panels).By contrast, as the LIS gets higher, ablation increases, leading to a smaller FIS.Back to a more realistic glacial context in which both North American and Eurasian ice sheets exist (xxdhLFIS experiments), we also simulate a net decrease of snow accumulation under full LGM LIS conditions (with respect to flat LIS) over the Kara sector combined with an increase over the western part of FIS (Fig. 17d), consistently with the upslope precipitation effect proposed by Sandberg and Oerlemans (1983).The large positive SMB values over the SBK sector when the LIS is low combined with the westward displacement of the Eurasian ice sheet as the LIS gets higher are fully in line with the ice-sheet paleotopography reconstructions provided by Kleman et al. (2013) for the last glacial cycle.The remote interactions highlighted in the present study between the LIS and the FIS systems also suggest that the ablation increase over the Barents and Kara seas when the LIS has its LGM size may have contributed to the LGM FIS small size.Conversely, the larger size of the FIS (Svendsen et al., 2004) during the late Saalian (around 140 ka) with respect to the LGM may have been related to a smaller LIS.
The remote interactions highlighted in the present study between the LIS and the FIS systems also suggest that the larger size of the FIS (Svendsen et al., 2004) during the late Saalian (around 140 ka) with respect to the LGM may have been related to a smaller LIS.
Nevertheless, we must keep in mind that all experiments presented in this study have used the LGM land-sea mask and therefore the Barents-Kara region is considered to be a continental area.The consequences of this assumption on the simulated large-scale atmospheric circulation are likely to be negligible.However, accounting for a marine Barents-Kara area would decrease the positive albedo feedback in this region, and thus the FIS inception (see Sect. 3.1.1).Our experimental setup uses LGM SST and sea-ice boundary conditions provided by the coupled atmosphere-ocean IPSL model (Kageyama et al., 2013).This prevents us from studying the conditions that could have prevailed with a marine Barents-Kara area, because this sector was covered by the FIS in this run, and future studies will be necessary to further investigate the relationships between the LIS and the FIS within more realistic contexts.

Ocean-atmosphere interactions
Our experimental setup has been designed to investigate how the atmospheric circulation changes in response to the LIS albedo and topography effects may impact the climate over the Eurasian region and therefore the evolution of the surface mass balance of the Fennoscandian ice sheet.We have therefore neglected the potential impact of the ocean and of its interactions with the atmosphere Yet, changes in atmospheric circulation related to the presence of continental ice sheets may induce large modifications of the oceanic circulation.For example, Wunsch (2006) proposed a scenario in which phase shifts of the stationary waves caused by the development of ice sheets modify the wind direction and strength leading to abrupt changes in temperature and precipitation fields (and thus to climatic events having similar features as Dansgaard-Oeschger events).These changes eventually produce large-scale modifications of the oceanic circulation.In turn, any change in the large-scale oceanic circulation may influence the evolution of the Northern Hemisphere ice sheets as well as their mutual interactions through the atmosphere-ocean coupling.
Recently, using the Earth system model COSMOS, Zhang et al. (2014) found that a gradual increase of the Northern Hemisphere ice-sheet height is able of causing abrupt transitions of the Atlantic meridonal overturning circulation through a positive atmosphere-ocean-sea-ice feedback involving the North Atlantic westerlies, the barotropic ocean circulation and sea-ice cover in the North Atlantic ocean.In this study, abrupt AMOC transitions from a weak to a strong mode, occurring for intermediate ice-sheet sizes, are accompanied by a noticeable increase of surface air temperatures in the northeastern North Atlantic.This result suggests that the warming effect simulated in our study over the Eurasian region when the LIS gets higher could be therefore amplified through North Atlantic warmer temperatures if oceanatmosphere interactions were accounted for.However, the mechanism highlighted in Zhang et al. (2014) is primarily driven by the simulated northward shift of the westerlies when the Northern Hemisphere ice sheet altitude gradually increases.Under full LGM conditions, a split of the northern westerlies in two zonal-oriented branches appears around the Laurentide ice sheet.The southward displacement of the jet stream simulated in our study when the LIS altitude grad-ually increases well supports the existence of the southern branch as shown in Zhang et al. (2014) for the LGM, but their northern branch remains absent from our all simulations.
Actually, this comparison illustrates the diversity of model responses to changes in the LIS characteristics.The recent inter-comparison study of the North Atlantic jet stream in LGM vs. pre-industrial conditions in the PMIP3 coupled models (Braconnot et al., 2012) carried out by Beghin et al. (2015) clearly shows that the response of the jet stream to LGM conditions is highly model-dependent.This latter study confirms the range of various results found in the recent literature.Using the outputs of four comprehensive atmosphereocean models, Laîné et al. (2009) found a southward displacement of the North Atlantic jet stream under LGM conditions.A similar conclusion can be drawn from the study of Ullman et al. (2014) carried out with the GISS model and two different LIS topographies as LGM ice-sheet boundary conditions: the difference of the simulated jet speed obtained with these different LIS reconstructions (see their Fig. 2) shows a more southward position of the jet when the LIS is higher.As was found in our study, Löfverström et al. (2014) found an equatorward position of the winter zonal winds and a more zonal jet under LGM conditions (with respect to MIS5b and MIS4).On the contrary, Roberts et al. (2014) found that a more zonal structure is simulated with the FA-MOUS model after the occurrence of a Heinrich event, which in their experimental set-up corresponds to a decrease of the Laurentide ice thickness.This indicates a southward shift of the North Atlantic jet stream when the LIS altitude is lower, and, by contrast, a northward displacement when it is higher.
Besides these contrasted findings, other studies did not find any significant influence of the ocean on the atmospheric circulation.As an example, by using the coupled IPSL-CM4 model, Pausata et al. (2011) showed that the position of the North Atlantic jet stream is not substantially influenced by sea-surface temperatures or sea-ice cover, but rather by icesheet height.In the same way, using the CAM3 atmospheric model coupled to a slab ocean, Löfverström et al. (2014) did not find any significant changes in the simulated temperature patterns over Alaska, Asia and the North Atlantic sector when two different oceanic heat transports corresponding to interglacial and LGM conditions are used.
The variety of the model responses calls for further studies to investigate why the results are so model-dependent.It also underlines the need for a more in-depth investigation of the links between large-scale oceanic and atmospheric circulations with different models run under climatic conditions rather different from those of the present-day periods.

Conclusions
The goal of this study was to investigate the atmospheric processes through which the LIS can influence the northwestern Eurasian climate, and therefore how it can influence the FIS surface mass balance.We performed 12 simulations under LGM conditions, but with different ice-sheet configurations.In our baseline experiment, the LIS and the FIS are removed.In the first set of simulations, we only removed the FIS, and the LIS has different heights.In the final series of experiments, we considered the impact of the LIS height on a high FIS.We therefore identified the effect of the LIS albedo and topography on the northwestern European climate.We focused on the fields having an impact on ice-sheet surface mass balance, i.e. the temperature and the precipitation.We showed the following: 1.The LIS albedo acts to decrease the summer temperatures over the Fennoscandian sector.This decrease is amplified by positive snow albedo feedback.This positive feedback is somewhat related to our experimental setup, which imposes land over the Barents-Kara sector instead of water, allowing snow maintenance.The temperature decrease is weaker when the LIS is higher, due to atmospheric circulation changes.
2. In winter, the LIS albedo impact decreases the Barents-Kara temperature, but increases the Scandinavian temperature.The higher the LIS, the colder the temperature over the whole Fennoscandian sector.
3. The effect of LIS albedo also tends to shift the North Atlantic jet stream poleward, and to bring more precipitation and snowfall over Scandinavia.When the LIS is higher, the Atlantic jet shifts equatorward, bringing less precipitation over northern Europe.As a consequence, precipitation and snowfall decrease over Scandinavia and Barents-Kara sectors when the LIS is higher.
4. The analysis of the GRISLI simulations forced by the LMDZ precipitation and temperature fields shows that the LIS albedo favors the growth of the FIS, essentially because of weak ablation.As the LIS gets higher, the FIS gets smaller due to more ablation, in accordance with the LMDZ5 summer temperature response.When the LIS reaches its full-LGM size, ice remains over Svalbard only.
5. Within the context of glacial-interglacial cycles, the FIS inception under full-LGM LIS conditions appears to be unrealistic.Therefore, we performed additional simulations with a high FIS and flat and full-LGM LIS.The results of these new sensitivity experiments show a similar response to the first series of experiments: in the presence of a full LIS, the FIS undergoes a significant mass loss.This confirms that a large LIS counteracts the growth of the FIS although it is still maintained, such as it does within more realistic contexts, and highlights therefore the self-consistency of the overall study.
This study highlights the mechanisms by which the LIS can influence the FIS surface mass balance.The use of P. Beghin et al.: North American ice sheet influence on the Eurasian climate the ice-sheet model allows for the quantification of the impact of climatic changes due to the LIS presence on the FIS.It shows that the albedo of the LIS favors the FIS growth, whereas the LIS topography acts against the FIS development.This conclusion is reinforced by the results obtained with a high FIS.It is worth noting that without a transient simulation of the last glacial period, it is not possible to accurately elaborate the leadlag relationships of the Northern Hemisphere ice sheets during the last glacial period.However, based on our results, a plausible glaciation scenario can be proposed with an early phase of FIS glacial inception favored by the areal expansion of the LIS, and followed by a limited growth of the FIS with a progressive westward retreat of the FIS itself as the LIS altitude gradually increased.Such proposed relationships between ice-sheet build-up demonstrate that this study also provides new insights in the understanding of the various Northern Hemisphere ice-sheet configurations during the different glacial-interglacial cycles.

Figure 1 .
Figure 1.Ice-sheet and oceanic boundary conditions used for the LMDZ experiments.(a) Ice-sheet height (light blue scale) and land height (brown scale) given as boundary conditions for the 5 selected experiments out of 12. Dark blue-brown limits represent the sea-land mask, and brown-light blue limits represent the land-ice mask.(b) Averaged annual sea-surface temperature (shading) and sea-ice extent for boreal summer (red contours) and winter (blue contours) from the PMIP3 IPSL-CM5A-LR simulations.The averages are computed over 50 years which have been used to create the boundary conditions for the LMDZ experiments.

Figure 3 .
Figure 3. 2 m summer (left) and winter (right) temperature without the LIS (noIS) and anomaly of the 2 m temperature between the xxdhL experiments and the noIS experiment; we can therefore quantify the impact of the LIS presence, depending on its topography.The visible surface albedo difference between xxdhL and noIS experiments are represented by orange contours (albedo difference = 0.4) and yellow contours (albedo difference = 0.2).The statistical significance of the patterns of temperature anomalies has been checked using a two-sided Student t test with a p value = 0.05.The white shading indicates areas where the anomaly is not significant.

Figure 5 .
Figure 5. Same as Fig. 4 for the winter season.

Figure 6 .
Figure 6.Same as Fig. 3 for annual precipitation (left column) and annual snowfall (right column).

Figure 7 .
Figure 7. Large-scale (left panels) and convective (right panels) annual precipitation without LIS (noIS experiment, top panels), and anomalies of large-scale and convective precipitation between the xxdhL and noIS experiments.The statistical significance of the patterns of temperature anomalies has been checked using a two-sided Student-t test with a p value = 0.05.The white shading indicates areas where the anomaly is not significant.

Figure 8 .
Figure 8.(a) Scandinavian and (b) Barents-Kara precipitation anomalies (xxdhL-noIS) during winter (circles) and summer (triangles) vs. annual precipitation anomaly.The lower the Laurentide ice-sheet altitude, the lighter grey the symbols are.

Figure 9 .
Figure 9. (shaded) and position of the North Atlantic jet stream without LIS (dashed) and for the xxdhL experiment (solid line).Right column: summer precipitation (shaded) and snowfall (contour every 0.1 mm day −1 , blue contour for positive values).

PFigure 10 .
Figure10.850 hPa winter jet-stream shift (in latitude) vs. winter precipitation anomaly over the Scandinavian region.The range of results has been obtained by bootstrapping.Up and down triangles represent the 5th precipitation percentile vs. the 95th jet position percentile, and the 95th precipitation percentile vs. the 5th jet position percentile respectively.Circles are the average values and crosses indicate the 5th and 95th percentiles.

Figure 11 .
Figure 11.Meridional cross sections of the 500hPa geopotential height anomaly from zonal means over the East Laurentide (90-60 • W, top panel) and the Atlantic sectors (50-10 • W, middle) in meters, and of the 500 hPa zonal wind (in m s −1 ) for the Atlantic sector.

Figure 12 .
Figure 12.(a) Snow accumulation computed by the GRISLI ice-sheet model for the noIS experiment (top panel), and differences of snow accumulation between the xxdhL and the noIS experiments (in m yr −1 ).(b) same as (a) for the ablation (in m yr −1 ).

Figure 13 .
Figure 13.(a) Surface mass balance at the beginning of the simulation (in m yr −1 ) and (b) surface elevation at the end of the simulation computed by the ice-sheet model for the selected experiments (in meters).

Figure 15 .
Figure 15.Same as Fig. 14 for the vertically integrated meridional (shaded) and the horizontal heat transport (arrows) in summer.Numerical values are given in 10 10 J m −1 s −1 .

Figure 16 .
Figure 16.Same as Fig. 14 for the 850 hPa meridional wind (shaded) and the 850hPa horizontal winds (arrows) for the summer season (in m s −1 ).

Figure 18 .
Figure 18.Same as Fig. 14 for the precipitation (mm day −1 ) and 850 hPa horizontal winds in winter.

Figure 19 .
Figure 19.Difference of simulated ice thickness of the Eurasian ice sheet between 00d and 100dhLFIS after 100 kyr (in meters).