Identifying drivers controlling the synchronicity of Heinrich-type ice sheet surges from the European and North American ice sheets

Author(s):  
Clemens Schannwell ◽  
Marie-Luise Kapsch ◽  
Uwe Mikolajewicz ◽  
Florian Ziemen

<p>Heinrich-type ice sheet surge events are among the most prominent signals in the paleoclimate data records. Even though these events have previously been intensely studied, it still remains an open question whether the cyclic ice sheet surges are triggered by internal ice dynamics, climate forcing, or a combination of the two. In simulations of the last deglaciation using the fully-coupled Max Planck Institute Earth System Model, surges from the European and North American ice sheets often occur in synchronicity. This model behaviour is in agreement with observations from sediment cores that find a similar pattern in the isotopic fingerprint of the deposited ice-rafted detritus. The synchronicity indicates that climate forcing is playing an important role in initiating ice sheet surges. In this study, we use the coupled ice-sheet-solid earth model PISM-VILMA in a northern hemispheric setup to investigate the modelled synchronicity of the surge events. More specifically, we perform an ensemble of simulations to study if the modelled synchronicity is a direct result of one of the surge locations causing other surge locations to be a activated as well. Moreover, we aim to investigate whether previously suggested trigger mechanisms such as regional changes in sea level or ocean temperatures are indeed key processes in controlling the synchronicity of these surge events.</p>

2018 ◽  
Author(s):  
Niall Gandy ◽  
Lauren J. Gregoire ◽  
Jeremy C. Ely ◽  
Christopher D. Clark ◽  
David M. Hodgson ◽  
...  

Abstract. Uncertainties in future sea level projections are dominated by our limited understanding of the dynamical processes that control instabilities of marine ice sheets. A valuable case to examine these processes is the last deglaciation of the British-Irish Ice Sheet. The Minch Ice Stream, which drained a large proportion of ice from the northwest sector of the British-Irish Ice Sheet during the last deglaciation, is well constrained, with abundant empirical data which could be used to inform, validate and analyse numerical ice sheet simulations. We use BISICLES, a higher-order ice sheet model, to examine the dynamical processes that controlled the retreat of the Minch Ice Stream. We simulate retreat from the shelf edge under constant "warm" surface mass balance and subshelf melt, to isolate the role of internal ice dynamics from external forcings. The model simulates a slowdown of retreat as the ice stream becomes laterally confined at a "pinning-point" between mainland Scotland and the Isle of Lewis. At this stage, the presence of ice shelves became a major control on deglaciation, providing buttressing to upstream ice. Subsequently, the presence of a reverse slope inside the Minch Strait produces an acceleration in retreat, leading to a "collapsed" state, even when the climate returns to the initial "cold" conditions. Our simulations demonstrate the importance of the Marine Ice Sheet Instability and ice shelf buttressing during the deglaciation of parts of the British-Irish Ice Sheet. Thus, geological data could be used to constrain these processes in ice sheet models used for projecting the future of our contemporary ice sheets.


2002 ◽  
Vol 48 (160) ◽  
pp. 70-80 ◽  
Author(s):  
Gerard H. Roe

AbstractThe interaction between ice sheets and the rest of the climate system at long time-scales is not well understood, and studies of the ice ages typically employ simplified parameterizations of the climate forcing on an ice sheet. It is important therefore to understand how an ice sheet responds to climate forcing, and whether the reduced approaches used in modeling studies are capable of providing robust and realistic answers. This work focuses on the accumulation distribution, and in particular considers what features of the accumulation pattern are necessary to model the steady-state response of an ice sheet. We examine the response of a model of the Greenland ice sheet to a variety of accumulation distributions, both observational datasets and simplified parameterizations. The predicted shape of the ice sheet is found to be quite insensitive to changes in the accumulation. The model only differs significantly from the observed ice sheet for a spatially uniform accumulation rate, and the most important factor for the successful simulation of the ice sheet’s shape is that the accumulation decreases with height according to the ability of the atmosphere to hold moisture. However, the internal ice dynamics strongly reflects the influence of the atmospheric circulation on the accumulation distribution.


2021 ◽  
Author(s):  
Uwe Mikolajewicz ◽  
Olga Erokhina ◽  
Marie-Luise Kapsch ◽  
Clemens Schannwell ◽  
Florian Ziemen

<p>It is challenging to model the last deglaciation, as it is characterized by abrupt millennial scale climate events, such as ice-sheet surges, that are superimposed on long-term climate changes, such as a global warming and the decay of a substantial part of the glacial ice sheets. Within PMIP, several groups have simulated the last deglaciation with CMIP-type models prescribing ice sheets from reconstructions. Whereas this type of simulations accounts for the effects of ice-sheet changes including meltwater release on climate, the prescribed ice sheet evolution is typically not consistent with the simulated climate evolution. Here we present a set of deglacial simulations that include fully interactive ice sheets that respond to changes in the climate. The setup also allows for feedbacks between ice sheets and climate and , hence, allows for a more realistic representation of the mechanisms of the last deglaciation, as the simulated climate and ice sheet changes are fully consistent..</p><p>The model consists of the coarse resolution set-up of MPI-ESM coupled to the ice sheet model mPISM (Northern Hemisphere and Antarctica) and the solid earth model VILMA. The model includes interactive icebergs and an automated calculation of the land-sea mask and river routing directions. A set of synchronously coupled simulations were started from an asynchronously coupled spin-up at 26ky and integrated throughout the deglaciation into the Holocene. The only prescribed external forcing are atmospheric concentrations of greenhouse gases and earth orbital parameters. One goal of this ensemble was to find the optimal combination of model parameters for the simulation of the deglaciation.</p><p>The model simulates the decay of the ice sheets, the rise of sea level, the flooding of shelf seas and the opening of passages. A large fraction of the ice sheet retreat is due to dynamical events (e.g. the final decay of the ice sheets on Barents Shelf or the Hudson Bay). Superimposed on the relatively slow glacial/interglacial transition are abrupt climate changes, triggered for example by recurrent ice sheet surges. These surges correspond to Heinrich Events tand result in a weakening of the AMOC. Three source regions for ice sheet surges occur during these simulations: from the Laurentide ice sheet through Hudson Strait, from the Laurentide ice sheet northward directly to the Arctic ocean, and from the Fennoscandian ice sheet into the Norwegian Sea. The characteristic climate response shows a large dependence on the surge location.</p><p>The simulated changes in strength of the AMOC are except for millennial-scale reduction events only moderate. However, during glacial periods, brine release is the central process for deep water formation in both hemispheres, in contrast to the Holocene. dDuring the deglaciation the ventilation of the deep ocean is strongly reduced, leading to a strong increase of the simulated deep water ages. This effect lasts longest in the deep North Pacific and extends in some simulations into the Holocene.</p>


2018 ◽  
Vol 12 (11) ◽  
pp. 3635-3651 ◽  
Author(s):  
Niall Gandy ◽  
Lauren J. Gregoire ◽  
Jeremy C. Ely ◽  
Christopher D. Clark ◽  
David M. Hodgson ◽  
...  

Abstract. Uncertainties in future sea level projections are dominated by our limited understanding of the dynamical processes that control instabilities of marine ice sheets. The last deglaciation of the British–Irish Ice Sheet offers a valuable example to examine these processes. The Minch Ice Stream, which drained a large proportion of ice from the northwest sector of the British–Irish Ice Sheet during the last deglaciation, is constrained with abundant empirical data which can be used to inform, validate, and analyse numerical ice sheet simulations. We use BISICLES, a higher-order ice sheet model, to examine the dynamical processes that controlled the retreat of the Minch Ice Stream. We perform simplified experiments of the retreat of this ice stream under an idealised climate forcing to isolate the effect of marine ice sheet processes, simulating retreat from the continental shelf under constant “warm” surface mass balance and sub-ice-shelf melt. The model simulates a slowdown of retreat as the ice stream becomes laterally confined at the mouth of the Minch strait between mainland Scotland and the Isle of Lewis, resulting in a marine setting similar to many large tidewater glaciers in Greenland and Antarctica. At this stage of the simulation, the presence of an ice shelf becomes a more important control on grounded ice volume, providing buttressing to upstream ice. Subsequently, the presence of a reverse slope inside the Minch strait produces an acceleration in retreat, leading to a “collapsed” state, even when the climate returns to the initial “cold” conditions. Our simulations demonstrate the importance of the marine ice sheet instability and ice shelf buttressing during the deglaciation of parts of the British–Irish Ice Sheet. We conclude that geological data could be applied to further constrain these processes in ice sheet models used for projecting the future of contemporary ice sheets.


1992 ◽  
Vol 338 (1285) ◽  
pp. 235-242 ◽  

The prediction of short-term (100 year) changes in the mass balance of ice sheets and longer-term (1000 years) variations in their ice volumes is important for a range of climatic and environmental models. The Antarctic ice sheet contains between 24 M km 3 and 29 M km 3 of ice, equivalent to a eustatic sea level change of between 60m and 72m. The annual surface accumulation is estimated to be of the order of 2200 Gtonnes, equivalent to a sea level change of 6 mm a -1 . Analysis of the present-day accumulation regime of Antarctica indicates that about 25% ( ca. 500 Gt a -1 ) of snowfall occurs in the Antarctic Peninsula region with an area of only 6.8% of the continent. To date most models have focused upon solving predictive algorithms for the climate-sensitivity of the ice sheet, and assume: (i) surface mass balance is equivalent to accumulation (i.e. no melting, evaporation or deflation); (ii) percentage change in accumulation is proportional to change in saturation mixing ratio above the surface inversion layer; and (iii) there is a linear relation between mean annual surface air tem perature and saturation mixing ratio. For the A ntarctic Peninsula with mountainous terrain containing ice caps, outlet glaciers, valley glaciers and ice shelves, where there can be significant ablation at low levels and distinct climatic regimes, models of the climate response must be more complex. In addition, owing to the high accumulation and flow rates, even short- to medium -term predictions must take account of ice dynamics. Relationships are derived for the mass balance sensitivity and, using a model developed by Hindmarsh, the transient effects of ice dynamics are estimated. It is suggested that for a 2°C rise in mean annual surface tem perature over 40 years, ablation in the A ntarctic Peninsula region would contribute at least 1.0 mm to sea level rise, offsetting the fall of 0.5 mm contributed by increased accum ulation.


2015 ◽  
Vol 11 (10) ◽  
pp. 1467-1490 ◽  
Author(s):  
P. Beghin ◽  
S. Charbit ◽  
C. Dumas ◽  
M. Kageyama ◽  
C. Ritz

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 full-size 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.


The Holocene ◽  
2018 ◽  
Vol 28 (7) ◽  
pp. 1160-1172 ◽  
Author(s):  
Samuel E Kelley ◽  
Jason P Briner ◽  
Sandy L O’Hara

The observational record of ice margin position reveals asynchrony in both the timing and magnitude of Greenland Ice Sheet (GrIS) margin fluctuations and illustrates the complex reactions of ice sheets to climatic perturbations. In this study, we reconstruct the timing and pattern of middle- and late-Holocene GrIS margin fluctuations at two locations, ~190 km apart, in central West Greenland using radiocarbon-dated sediment cores from proglacial-threshold lakes. Our results demonstrate that deglaciation occurs at both sites during the early Holocene, with the ice sheet remaining in a smaller-than-present ice margin configuration until ~500 years ago when it readvanced into lake catchments at both sites. At our northern site, Sermeq Kujatdleq, the late-Holocene advance of the GrIS approached maximum position during the past 280 years, with the culmination of the advance occurring at AD 1992–1994, and modern retreat was underway by AD 1998–2001. In contrast, field and observational evidence suggest that the GrIS at our southern site, Nordenskiöld Gletscher, has been advancing or stable throughout the 20th century. These results, in conjunction with previous work in the region, highlight the asynchronous nature of late-Holocene advances and subsequent modern retreat, implying that local variability, such as ice velocity or ice dynamics, is responsible for modulating ice margin response to changes in climate on these decadal to centennial timescales. Additional high-resolution records of past ice sheet fluctuations are needed to inform and more accurately constrain our predictions of future cryosphere response to changes in climate.


2021 ◽  
Author(s):  
Marlene Klockmann ◽  
Marie-Luise Kapsch ◽  
Uwe Mikolajewicz

<p><span>Coupled climate models have produced very different states of the Atlantic Meridional Overturning Circulation (AMOC) in simulations of the Last Glacial Maximum (LGM). In particular, many of them failed to capture the shoaling of the upper AMOC cell, which was indicated by reconstructions. In sensitivity simulations with the Max-Planck-Institute Earth System Model (MPI-ESM) we found that the glacial AMOC response is the sum of two large opposing effects: a strengthening and deepening of the upper cell in response to the glacial ice sheets and a weakening and shoaling of the upper cell in response to the low glacial greenhouse gas concentrations. The magnitude of the respective effects likely depends on the background climate, the ice sheet reconstruction used, and model specifics such as the representation of brine release in the Southern Ocean. </span></p><p><span>Transient simulations of the deglaciation with two differently tuned versions of MPI-ESM and two different ice-sheet reconstructions differ strongly in their respective AMOC states during the LGM. These simulations, together with selected PMIP3 and PMIP4 LGM simulations, provide a good opportunity to compare the effect of different ice sheet reconstructions on the glacial AMOC. We compare key variables such as water mass properties, salt transport and Southern Ocean sea-ice formation across this ensemble of opportunity with the aim of increasing our understanding of the role of ice sheets in the glacial AMOC response.</span></p>


1995 ◽  
Vol 21 ◽  
pp. 103-110 ◽  
Author(s):  
G. S. Boulton ◽  
N. Hulton ◽  
M. Vautravers

A numerical model is used to simulate ice-sheet behaviour in Europe through the last glacial cycle. It is used in two modes: a forward mode, in which the model is driven by a proxy palaeoclimate record and the output compared with a geological reconstruction of ice-sheet fluctuation; and an inverse mode, in which we determine the climate function that would be required to simulate geologically reconstructed ice-sheet fluctuations. From these simulations it is concluded that extra-glacial climates may be poor predictors of ice-sheet surface climates, and that climatic transitions during the glacial period may have been much more rapid and the intensity of warming during the early Holocene much greater than hitherto supposed. Stronger climate forcing is required to drive ice-sheet expansion when sliding occurs at the bed compared with a non-sliding bed. Sliding ice sheets grow more slowly and decay more rapidly than non-sliding ice sheets with the same climate forcing.


Sign in / Sign up

Export Citation Format

Share Document