scholarly journals Climate and ice sheet evolutions from the last glacial maximum to the pre-industrial period with an ice sheet – climate coupled model

2021 ◽  
Author(s):  
Aurélien Quiquet ◽  
Didier M. Roche ◽  
Christophe Dumas ◽  
Nathaëlle Bouttes ◽  
Fanny Lhardy

Abstract. The last deglaciation offers an unique opportunity to understand the climate – ice sheet interactions in a global warming context. In this paper, to tackle this question, we use an Earth system model of intermediate complexity coupled to an ice sheet model covering the Northern Hemisphere to simulate the last deglaciation and the Holocene (26–0 ka BP). We use a synchronous coupling every year between the ice sheet and the rest of the climate system and we ensure a closed water cycle considering the release of freshwater flux to the ocean due to ice sheet melting. Our reference experiment displays a gradual warming in response to the forcings, with no abrupt changes. In this case, while the amplitude of the freshwater flux to the ocean induced by ice sheet retreat is realistic, it is sufficient to shut down the Atlantic meridional overturning from which the model does not recover within the time period simulated. However, with reduced freshwater flux we are nonetheless able to obtain different oceanic circulation evolutions, including some abrupt transitions between shut-down and active circulation states in the course of the deglaciation. The fast oceanic circulation recoveries lead to abrupt warming phases in Greenland. Our simulated ice sheet geometry evolution is in overall good agreement with available global reconstructions, even though the abrupt sea level rise at 14.6 kaBP is underestimated, possibly because the climate model underestimates the millenial- scale temperature variability. In the course of the deglaciation, large-scale grounding line instabilities are simulated both for the Eurasian and North American ice sheets. The first instability occurs in the Barents-Kara seas for the Eurasian ice sheet at 14.5 kaBP. A second grounding line instability occurs circa 12 kaBP in the proglacial lake that formed at the southern margin of the North American ice sheet. With additional asynchronously coupled experiments, we assess the sensitivity of our results to different ice sheet model choices related to surface and sub-shelf mass balance, ice deformation and grounding line representation. While the ice sheet evolutions differ within this ensemble, the global climate trajectory is only weakly affected by these choices. In our experiments, only the abrupt shifts in the oceanic circulation due to freshwater fluxes are able to produce some millenial-scale variability since no self-generating abrupt transitions are simulated without these fluxes.

2021 ◽  
Vol 17 (5) ◽  
pp. 2179-2199
Author(s):  
Aurélien Quiquet ◽  
Didier M. Roche ◽  
Christophe Dumas ◽  
Nathaëlle Bouttes ◽  
Fanny Lhardy

Abstract. The last deglaciation offers an unique opportunity to understand the climate–ice-sheet interactions in a global warming context. In this paper, to tackle this question, we use an Earth system model of intermediate complexity coupled to an ice sheet model covering the Northern Hemisphere to simulate the last deglaciation and the Holocene (26–0 ka). We use a synchronous coupling every year between the ice sheet and the rest of the climate system and we ensure a closed water cycle considering the release of freshwater flux to the ocean due to ice sheet melting. Our reference experiment displays a gradual warming in response to the forcings, with no abrupt changes. In this case, while the amplitude of the freshwater flux to the ocean induced by ice sheet retreat is realistic, it is sufficient to shut down the Atlantic meridional overturning circulation from which the model does not recover within the time period simulated. However, with reduced freshwater flux we are nonetheless able to obtain different oceanic circulation evolutions, including some abrupt transitions between shut-down and active circulation states in the course of the deglaciation. The inclusion of a parameterisation for the sinking of brines around Antarctica also produces an abrupt recovery of the Atlantic meridional overturning circulation, absent in the reference experiment. The fast oceanic circulation recoveries lead to abrupt warming phases in Greenland. Our simulated ice sheet geometry evolution is in overall good agreement with available global reconstructions, even though the abrupt sea level rise at 14.6 ka is underestimated, possibly because the climate model underestimates the millennial-scale temperature variability. In the course of the deglaciation, large-scale grounding line instabilities are simulated both for the Eurasian and North American ice sheets. The first instability occurs in the Barents–Kara seas for the Eurasian ice sheet at 14.5 ka. A second grounding line instability occurs ca. 12 ka in the proglacial lake that formed at the southern margin of the North American ice sheet. With additional asynchronously coupled experiments, we assess the sensitivity of our results to different ice sheet model choices related to surface and sub-shelf mass balance, ice deformation and grounding line representation. While the ice sheet evolutions differ within this ensemble, the global climate trajectory is only weakly affected by these choices. In our experiments, only the abrupt shifts in the oceanic circulation due to freshwater fluxes are able to produce some millennial-scale variability since no self-generating abrupt transitions are simulated without these fluxes.


2020 ◽  
Vol 234 ◽  
pp. 106223 ◽  
Author(s):  
April S. Dalton ◽  
Martin Margold ◽  
Chris R. Stokes ◽  
Lev Tarasov ◽  
Arthur S. Dyke ◽  
...  

2019 ◽  
Vol 5 (8) ◽  
pp. eaav8754 ◽  
Author(s):  
Daniel P. Lowry ◽  
Nicholas R. Golledge ◽  
Nancy A. N. Bertler ◽  
R. Selwyn Jones ◽  
Robert McKay

Modern observations appear to link warming oceanic conditions with Antarctic ice sheet grounding-line retreat. Yet, interpretations of past ice sheet retreat over the last deglaciation in the Ross Embayment, Antarctica’s largest catchment, differ considerably and imply either extremely high or very low sensitivity to environmental forcing. To investigate this, we perform regional ice sheet simulations using a wide range of atmosphere and ocean forcings. Constrained by marine and terrestrial geological data, these models predict earliest retreat in the central embayment and rapid terrestrial ice sheet thinning during the Early Holocene. We find that atmospheric conditions early in the deglacial period can enhance or diminish ice sheet sensitivity to rising ocean temperatures, thereby controlling the initial timing and spatial pattern of grounding-line retreat. Through the Holocene, however, grounding-line position is much more sensitive to subshelf melt rates, implicating ocean thermal forcing as the key driver of past ice sheet retreat.


2012 ◽  
Vol 5 (5) ◽  
pp. 1273-1295 ◽  
Author(s):  
D. Pollard ◽  
R. M. DeConto

Abstract. The formulation of a 3-D ice sheet-shelf model is described. The model is designed for long-term continental-scale applications, and has been used mostly in paleoclimatic studies. It uses a hybrid combination of the scaled shallow ice and shallow shelf approximations for ice flow. Floating ice shelves and grounding-line migration are included, with parameterized ice fluxes at grounding lines that allows relatively coarse resolutions to be used. All significant components and parameterizations of the model are described in some detail. Basic results for modern Antarctica are compared with observations, and simulations over the last 5 million years are compared with previously published results. The sensitivity of ice volumes during the last deglaciation to basal sliding coefficients is discussed.


2021 ◽  
Vol 15 (7) ◽  
pp. 3329-3354
Author(s):  
Trevor R. Hillebrand ◽  
John O. Stone ◽  
Michelle Koutnik ◽  
Courtney King ◽  
Howard Conway ◽  
...  

Abstract. Chronologies of glacier deposits in the Transantarctic Mountains provide important constraints on grounding-line retreat during the last deglaciation in the Ross Sea. However, between Beardmore Glacier and Ross Island – a distance of some 600 km – the existing chronologies are generally sparse and far from the modern grounding line, leaving the past dynamics of this vast region largely unconstrained. We present exposure ages of glacial deposits at three locations alongside the Darwin–Hatherton Glacier System – including within 10 km of the modern grounding line – that record several hundred meters of Late Pleistocene to Early Holocene thickening relative to present. As the ice sheet grounding line in the Ross Sea retreated, Hatherton Glacier thinned steadily from about 9 until about 3 ka. Our data are equivocal about the maximum thickness and Mid-Holocene to Early Holocene history at the mouth of Darwin Glacier, allowing for two conflicting deglaciation scenarios: (1) ∼500 m of thinning from 9 to 3 ka, similar to Hatherton Glacier, or (2) ∼950 m of thinning, with a rapid pulse of ∼600 m thinning at around 5 ka. We test these two scenarios using a 1.5-dimensional flowband model, forced by ice thickness changes at the mouth of Darwin Glacier and evaluated by fit to the chronology of deposits at Hatherton Glacier. The constraints from Hatherton Glacier are consistent with the interpretation that the mouth of Darwin Glacier thinned steadily by ∼500 m from 9 to 3 ka. Rapid pulses of thinning at the mouth of Darwin Glacier are ruled out by the data at Hatherton Glacier. This contrasts with some of the available records from the mouths of other outlet glaciers in the Transantarctic Mountains, many of which thinned by hundreds of meters over roughly a 1000-year period in the Early Holocene. The deglaciation histories of Darwin and Hatherton glaciers are best matched by a steady decrease in catchment area through the Holocene, suggesting that Byrd and/or Mulock glaciers may have captured roughly half of the catchment area of Darwin and Hatherton glaciers during the last deglaciation. An ensemble of three-dimensional ice sheet model simulations suggest that Darwin and Hatherton glaciers are strongly buttressed by convergent flow with ice from neighboring Byrd and Mulock glaciers, and by lateral drag past Minna Bluff, which could have led to a pattern of retreat distinct from other glaciers throughout the Transantarctic Mountains.


2021 ◽  
Author(s):  
Natalya Gomez ◽  
Michael Weber ◽  
Peter Clark ◽  
Jerry Mitrovica ◽  
Holly Han

<p>A longstanding hypothesis for near-synchronous evolution of global ice sheets over ice-age cycles invokes an interhemispheric sea-level forcing whereby sea-level rise due to ice loss in the Northern Hemisphere in response to insolation and greenhouse gas forcing causes grounding-line retreat of marine-based sectors of the Antarctic Ice Sheet (AIS). Recent studies have shown that the AIS experienced substantial millennial-scale variability during and after the last deglaciation, with several times of recorded increased iceberg flux and grounding line retreat coinciding, within uncertainty, with well documented global sea-level rise events, providing further evidence of this sea-level forcing. However, the sea level changes associated with ice sheet mass loss are strongly nonuniform due to gravitational, deformational and Earth rotational effects, suggesting that the response of the AIS to Northern Hemisphere sea-level forcing is more complicated than previously modelled.</p><p>We adopt an ice-sheet model coupled to a global sea-level model to show that a large or rapid Northern Hemisphere sea-level forcing enhances grounding-line advance and associated mass gain of the AIS during glaciation, and grounding-line retreat and AIS mass loss during deglaciation. Relative to models without these interactions, including the Northern Hemisphere sea-level forcing leads to a larger AIS volume during the Last Glacial Maximum (about 26,000 to 20,000 years ago), subsequent earlier grounding-line retreat and millennial-scale AIS variability throughout the last deglaciation. These findings are consistent with geologic reconstructions of the extent of the AIS during the Last Glacial Maximum and subsequent ice-sheet retreat, and with relative sea-level change in Antarctica. </p>


2012 ◽  
Vol 5 (2) ◽  
pp. 1077-1134 ◽  
Author(s):  
D. Pollard ◽  
R. M. DeConto

Abstract. The formulation of a 3-D ice sheet-shelf model is described. The model is designed for long-term continental-scale applications, and has been used mostly in paleoclimatic studies. It uses a hybrid combination of the scaled Shallow Ice and Shallow Shelf Approximations for ice flow. Floating ice shelves and grounding-line migration are included, with parameterized ice fluxes at grounding lines that allows relatively coarse resolutions to be used. All significant components and parameterizations of the model are described in some detail. Basic results for modern Antarctica are compared with observations, and simulations over the last 5 million yr are shown to be similar to previously published results using an earlier model version. The sensitivity of ice retreat during the last deglaciation to basal sliding coefficients is discussed.


2019 ◽  
Vol 5 (12) ◽  
pp. eaaw2610 ◽  
Author(s):  
Margaret S. Jackson ◽  
Meredith A. Kelly ◽  
James M. Russell ◽  
Alice M. Doughty ◽  
Jennifer A. Howley ◽  
...  

Atmospheric greenhouse gas concentrations are thought to have synchronized global temperatures during Pleistocene glacial–interglacial cycles, yet their impact relative to changes in high-latitude insolation and ice-sheet extent remains poorly constrained. Here, we use tropical glacial fluctuations to assess the timing of low-latitude temperature changes relative to global climate forcings. We report 10Be ages of moraines in tropical East Africa and South America and show that glaciers reached their maxima at ~29 to 20 ka, during the global Last Glacial Maximum. Tropical glacial recession was underway by 20 ka, before the rapid CO2 rise at ~18.2 ka. This “early” tropical warming was influenced by rising high-latitude insolation and coincident ice-sheet recession in both polar regions, which lowered the meridional thermal gradient and reduced tropical heat export to the high latitudes.


Geosciences ◽  
2021 ◽  
Vol 11 (7) ◽  
pp. 266
Author(s):  
Wei Liu ◽  
Zhengyu Liu ◽  
Shouwei Li

We explore the change in Southern Ocean upwelling during the last deglaciation, based on proxy records and a transient climate model simulation. Our analyses suggest that, beyond a conventional mechanism of the Southern Hemisphere westerlies shift, Southern Ocean upwelling is strongly influenced by surface buoyancy forcing and the local topography. Over the Antarctic Circumpolar Current region, the zonal mean and local upwelled flows exhibited distinct evolution patterns during the last deglaciation, since they are driven by different mechanisms. The zonal mean upwelling is primarily driven by surface wind stress via zonal mean Ekman pumping, whereas local upwelling is driven by both wind and buoyancy forcing, and is tightly coupled to local topography. During the early stage of the last deglaciation, the vertical extension of the upwelled flows increased downstream of submarine ridges but decreased upstream, which led to enhanced and diminished local upwelling, downstream and upstream of the submarine ridges, respectively.


2017 ◽  
Vol 10 (3) ◽  
pp. 1383-1402 ◽  
Author(s):  
Paolo Davini ◽  
Jost von Hardenberg ◽  
Susanna Corti ◽  
Hannah M. Christensen ◽  
Stephan Juricke ◽  
...  

Abstract. The Climate SPHINX (Stochastic Physics HIgh resolutioN eXperiments) project is a comprehensive set of ensemble simulations aimed at evaluating the sensitivity of present and future climate to model resolution and stochastic parameterisation. The EC-Earth Earth system model is used to explore the impact of stochastic physics in a large ensemble of 30-year climate integrations at five different atmospheric horizontal resolutions (from 125 up to 16 km). The project includes more than 120 simulations in both a historical scenario (1979–2008) and a climate change projection (2039–2068), together with coupled transient runs (1850–2100). A total of 20.4 million core hours have been used, made available from a single year grant from PRACE (the Partnership for Advanced Computing in Europe), and close to 1.5 PB of output data have been produced on SuperMUC IBM Petascale System at the Leibniz Supercomputing Centre (LRZ) in Garching, Germany. About 140 TB of post-processed data are stored on the CINECA supercomputing centre archives and are freely accessible to the community thanks to an EUDAT data pilot project. This paper presents the technical and scientific set-up of the experiments, including the details on the forcing used for the simulations performed, defining the SPHINX v1.0 protocol. In addition, an overview of preliminary results is given. An improvement in the simulation of Euro-Atlantic atmospheric blocking following resolution increase is observed. It is also shown that including stochastic parameterisation in the low-resolution runs helps to improve some aspects of the tropical climate – specifically the Madden–Julian Oscillation and the tropical rainfall variability. These findings show the importance of representing the impact of small-scale processes on the large-scale climate variability either explicitly (with high-resolution simulations) or stochastically (in low-resolution simulations).


Sign in / Sign up

Export Citation Format

Share Document