scholarly journals Exploring the impact of atmospheric forcing and basal boundary conditions on the simulation of the Antarctic ice sheet at the Last Glacial Maximum

2020 ◽  
Author(s):  
Javier Blasco ◽  
Jorge Alvarez-Solas ◽  
Alexander Robinson ◽  
Marisa Montoya

Abstract. Little is known about the distribution of ice in the Antarctic ice sheet (AIS) during the Last Glacial Maximum (LGM). Whereas marine and terrestrial geological data indicate that the grounded ice advanced to a position close to the continental-shelf break, the total ice volume is unclear. Glacial boundary conditions are potentially important sources of uncertainty, in particular basal friction and climatic boundary conditions. Basal friction exerts a strong control on the large-scale dynamics of the ice sheet and thus affects its size, and is not well constrained. Glacial climatic boundary conditions determine the net accumulation and ice temperature, and are also poorly known. Here we explore the effect of the uncertainty in both features on the total simulated ice storage of the AIS at the LGM. For this purpose we use a hybrid ice-sheet-shelf model that is forced with different basal-drag choices and glacial background climatic conditions obtained from the LGM ensemble climate simulations of the third phase of the Paleoclimate Modelling Intercomparison Project (PMIP3). For a wide range of plausible basal friction configurations, the simulated ice dynamics vary widely but all simulations produce fully extended ice sheets towards the continental-shelf break. More dynamically active ice sheets correspond to lower ice volumes, while they remain consistent with the available constraints on ice extent. Thus, this work points to the possibility of an AIS with very active ice streams during the LGM. In addition, we find that the surface boundary temperature field plays a crucial role in determining the ice extent through its effect on viscosity. For ice sheets of a similar extent and comparable dynamics, we find that the precipitation field determines the total AIS volume. However, precipitation is deeply uncertain. Climatic fields simulated by climate models show more precipitation in coastal regions than a spatially uniform anomaly, which can lead to larger ice volumes. We strongly support using these paleoclimatic fields to simulate and study the LGM and potentially other time periods like the Last Interglacial. However, their accuracy must be assessed as well, as differences between climate model forcing lead to a range in the simulated ice volume and extension of about 6 m sea-level equivalent and one million km2.

2021 ◽  
Vol 15 (1) ◽  
pp. 215-231
Author(s):  
Javier Blasco ◽  
Jorge Alvarez-Solas ◽  
Alexander Robinson ◽  
Marisa Montoya

Abstract. Little is known about the distribution of ice in the Antarctic Ice Sheet (AIS) during the Last Glacial Maximum (LGM). Whereas marine and terrestrial geological data indicate that the grounded ice advanced to a position close to the continental-shelf break, the total ice volume is unclear. Glacial boundary conditions are potentially important sources of uncertainty, in particular basal friction and climatic boundary conditions. Basal friction exerts a strong control on the large-scale dynamics of the ice sheet and thus affects its size and is not well constrained. Glacial climatic boundary conditions determine the net accumulation and ice temperature and are also poorly known. Here we explore the effect of the uncertainty in both features on the total simulated ice storage of the AIS at the LGM. For this purpose we use a hybrid ice sheet shelf model that is forced with different basal drag choices and glacial background climatic conditions obtained from the LGM ensemble climate simulations of the third phase of the Paleoclimate Modelling Intercomparison Project (PMIP3). Overall, we find that the spread in the simulated ice volume for the tested basal drag parameterizations is about the same range as for the different general circulation model (GCM) forcings (4 to 6 m sea level equivalent). For a wide range of plausible basal friction configurations, the simulated ice dynamics vary widely but all simulations produce fully extended ice sheets towards the continental-shelf break. More dynamically active ice sheets correspond to lower ice volumes, while they remain consistent with the available constraints on ice extent. Thus, this work points to the possibility of an AIS with very active ice streams during the LGM. In addition, we find that the surface boundary temperature field plays a crucial role in determining the ice extent through its effect on viscosity. For ice sheets of a similar extent and comparable dynamics, we find that the precipitation field determines the total AIS volume. However, precipitation is highly uncertain. Climatic fields simulated by climate models show more precipitation in coastal regions than a spatially uniform anomaly, which can lead to larger ice volumes. Our results strongly support using these paleoclimatic fields to simulate and study the LGM and potentially other time periods like the last interglacial. However, their accuracy must be assessed as well, as differences between climate model forcing lead to a large spread in the simulated ice volume and extension.


1993 ◽  
Vol 341 (1297) ◽  
pp. 253-261 ◽  

A two-dimensional model which links the atmosphere, the mixed layer of the ocean, the sea ice, the continents, the ice sheets and their underlying bedrock has been used to test the Milankovitch theory over the last two glacial-interglacial cycles. A series of sensitivity analyses have allowed us to understand better the internal mechanisms which drive the simulated climate system and in particular the feedbacks related to surface albedo and water vapour. It was found that orbital variations alone can induce, in such a system, feedbacks sufficient to generate the low frequency p art of the climatic variations over the last 122 ka. These simulated variations at the astronomical timescale are broadly in agreement with reconstructions of ice-sheet volume and of sea level independently obtained from geological data. Imperfections in the stimulated climate were the insufficient southward extent of the ice sheets and the too small hemispheric cooling at the last glacial maximum . These deficiencies were partly remedied in a further experiment by using the time-dependent atmospheric CO2 concentration given by the Vostok ice core in addition to the astronomical forcing. In this transient simulation, 70% of the Northern Hemisphere ice volume is related to the astronomical forcing and the related changes in the albedo, the rem aining 30% being due to the CO 2 changes. Analysis of the processes involved shows that variations of ablation are more important for the ice-sheet response than are variations of snow precipitation. A key mechanism in the deglaciation after the last glacial maxim um appears to be the ‘ageing’ of snow which significantly decreases its albedo. The other factors which play an important role are ice-sheet altitude, insolation, taiga cover, ice-albedo feedback, ice-sheet configuration (‘continentality’ and ‘desert’ effect), isostatic rebound, CO 2 changes and tem perature-water vapour feedback. Numerical experiments have also been carried out with a one-dimensional radiative-convective model in order to quantify the influence of the CO 2 changes and of the water vapour feedback on the climate evolution of the Northern Hemisphere over the last 122 ka. Results of these experiments indicate that 67% of the simulated cooling at the last glacial maximum can be attributed to the astronomical forcing and the subsequent surface albedo increase, the remaining 33% being associated with the reduced CO 2 concentration. Moreover, the water vapour feedback explains 40% of the simulated cooling in all the experiments done. The transient response of the clim ate system to both the astronomical and CO 2 forcing was also simulated by the LLN (Louvain-la-Neuve) 2.5-dimensional model over the two last glacial-interglacial cycles. It is particularly significant that spectral analysis of the simulated Northern Hemisphere global ice volume variations reproduces correctly the relative intensity of the peaks at the orbital frequencies. Except for variations with timescales shorter than 5 ka, the simulated long-term variations of total ice volume are comparable to that reconstructed from deep sea cores. For example, the model simulates glacial maxima of similar amplitudes at 134 ka BP and 15 ka BP, followed by abrupt deglaciations. The complete deglaciation of the three main Northern Hemisphere ice sheets, which is simulated around 122 ka BP, is in partial disagreement with reconstructions indicating that the Greenland ice sheet survived during the Eemian interglacial. The continental ice volume variations during the last 122 ka of the 200 ka simulation are, however, not significantly affected by this shortcoming.


2021 ◽  
Author(s):  
Lauren Gregoire ◽  
Niall Gandy ◽  
Lachlan Astfalck ◽  
Robin Smith ◽  
Ruza Ivanovic ◽  
...  

<p>Simulating the co-evolution of climate and ice-sheets during the Quaternary is key to understanding some of the major abrupt changes in climate, ice and sea level. Indeed, events such as the Meltwater pulse 1a rapid sea level rise and Heinrich, Dansgaard–Oeschger and the 8.2 kyr climatic events all involve the interplay between ice sheets, the atmosphere and the ocean. Unfortunately, it is challenging to simulate the coupled Climate-Ice sheet system because small biases, errors or uncertainties in parts of the models are strongly amplified by the powerful interactions between the atmosphere and ice (e.g. ice-albedo and height-mass balance feedbacks). This leads to inaccurate or even unrealistic simulations of ice sheet extent and surface climate. To overcome this issue we need some methods to effectively explore the uncertainty in the complex Climate-Ice sheet system and reduce model biases. Here we present our approach to produce ensemble of coupled Climate-Ice sheet simulations of the Last Glacial maximum that explore the uncertainties in climate and ice sheet processes.</p><p>We use the FAMOUS-ICE earth system model, which comprises a coarse-resolution and fast general circulation model coupled to the Glimmer-CISM ice sheet model. We prescribe sea surface temperature and sea ice concentrations in order to control and reduce biases in polar climate, which strongly affect the surface mass balance and simulated extent of the northern hemisphere ice sheets. We develop and apply a method to reconstruct and sample a range of realistic sea surface temperature and sea-ice concentration spatio-temporal field. These are created by merging information from PMIP3/4 climate simulations and proxy-data for sea surface temperatures at the Last Glacial Maximum with Bayes linear analysis. We then use these to generate ensembles of FAMOUS-ice simulations of the Last Glacial maximum following the PMIP4 protocol, with the Greenland and North American ice sheets interactively simulated. In addition to exploring a range of sea surface conditions, we also vary key parameters that control the surface mass balance and flow of ice sheets. We thus produce ensembles of simulations that will later be used to emulate ice sheet surface mass balance.  </p>


2017 ◽  
Vol 13 (9) ◽  
pp. 1243-1257 ◽  
Author(s):  
Lennert B. Stap ◽  
Roderik S. W. van de Wal ◽  
Bas de Boer ◽  
Richard Bintanja ◽  
Lucas J. Lourens

Abstract. Since the inception of the Antarctic ice sheet at the Eocene–Oligocene transition (∼ 34 Myr ago), land ice has played a crucial role in Earth's climate. Through feedbacks in the climate system, land ice variability modifies atmospheric temperature changes induced by orbital, topographical, and greenhouse gas variations. Quantification of these feedbacks on long timescales has hitherto scarcely been undertaken. In this study, we use a zonally averaged energy balance climate model bidirectionally coupled to a one-dimensional ice sheet model, capturing the ice–albedo and surface–height–temperature feedbacks. Potentially important transient changes in topographic boundary conditions by tectonics and erosion are not taken into account but are briefly discussed. The relative simplicity of the coupled model allows us to perform integrations over the past 38 Myr in a fully transient fashion using a benthic oxygen isotope record as forcing to inversely simulate CO2. Firstly, we find that the results of the simulations over the past 5 Myr are dependent on whether the model run is started at 5 or 38 Myr ago. This is because the relation between CO2 and temperature is subject to hysteresis. When the climate cools from very high CO2 levels, as in the longer transient 38 Myr run, temperatures in the lower CO2 range of the past 5 Myr are higher than when the climate is initialised at low temperatures. Consequently, the modelled CO2 concentrations depend on the initial state. Taking the realistic warm initialisation into account, we come to a best estimate of CO2, temperature, ice-volume-equivalent sea level, and benthic δ18O over the past 38 Myr. Secondly, we study the influence of ice sheets on the evolution of global temperature and polar amplification by comparing runs with ice sheet–climate interaction switched on and off. By passing only albedo or surface height changes to the climate model, we can distinguish the separate effects of the ice–albedo and surface–height–temperature feedbacks. We find that ice volume variability has a strong enhancing effect on atmospheric temperature changes, particularly in the regions where the ice sheets are located. As a result, polar amplification in the Northern Hemisphere decreases towards warmer climates as there is little land ice left to melt. Conversely, decay of the Antarctic ice sheet increases polar amplification in the Southern Hemisphere in the high-CO2 regime. Our results also show that in cooler climates than the pre-industrial, the ice–albedo feedback predominates the surface–height–temperature feedback, while in warmer climates they are more equal in strength.


1997 ◽  
Vol 25 ◽  
pp. 145-152 ◽  
Author(s):  
Gilles Ramstein ◽  
Adeline Fabre ◽  
Sophie Pinot ◽  
Catherine Ritz ◽  
Sylvie Joussaume

In the framework of the Paleoclimate Modelling Intercomparison Project (PMIP), simulations of the Last Glacial Maximum (LGM) have- been performed. More than 10 different atmospheric general circulation models (AGCMs) have been used with the same boundary conditions: sea-surface temperatures prescribed by CLIMAP (1981), ice-sheet reconstruction provided by Peltier (1994), change in insolation, and reduced CO2 content. One of the major questions is to investigate whether the simulations of the LGM are in equilibrium with the prescribed ice-sheet reconstruction. To answer this question, we have used two different approaches. First, we analyze the results of a sel of LGM simulations performed with different versions of the Laboratoire de Meteorolo-gie Dynamique (LMD) AGCM and study the hydrologic and snow- budgets over the Laurcntide and Fennoscandian ice sheets. Second, we use the AGCM outputs to force an ice-sheet model in order to investigate its ability to maintain the ice sheets as reconstructed by CLIMAP (1981) or Peltier (1994).


2014 ◽  
Vol 100 ◽  
pp. 87-110 ◽  
Author(s):  
Colm Ó Cofaigh ◽  
Bethan J. Davies ◽  
Stephen J. Livingstone ◽  
James A. Smith ◽  
Joanne S. Johnson ◽  
...  

2001 ◽  
Vol 56 (3) ◽  
pp. 299-307 ◽  
Author(s):  
Isaac J. Winograd

AbstractThe magnitude of late Wisconsinan (post-35,000 yr B.P.) ice-sheet growth in the Northern Hemisphere is not well known. Ice volume at ∼35,000 yr B.P. may have been as little as 20% or as much as 70% of the volume present at the last glacial maximum (LGM). A conservative evaluation of glacial–geologic, sea level, and benthic δ18O data indicates that ice volume at ∼35,000 yr B.P. was approximately 50% of that extant at the LGM (∼20,000 yr B.P.); that is, it doubled in about 15,000 yr. On the basis of literature for the North Atlantic and a sea-surface temperature (SST) data compilation, it appears that this rapid growth may have been forced by low-to-mid-latitude SST warming in both the Atlantic and Pacific Oceans, with attendant increased moisture transport to high latitudes. The SST ice-sheet growth notion also explains the apparent synchroneity of late Wisconsinan mountain glaciation in both hemispheres.


2017 ◽  
Author(s):  
Sarah L. Bradley ◽  
Thomas J. Reerink ◽  
Roderik S. W. van de Wal ◽  
Michiel M. Helsen

Abstract. Observational evidence, including offshore moraines and sediment cores confirm that at the Last Glacial maximum (LGM) the Greenland ice sheet (GrIS) grew to a significantly larger spatial extent than seen at present, grounding into Baffin Bay and to the continental shelf break. Given this larger spatial extent and it is close proximity to the neighboring Laurentide (LIS) and Innuitian Ice sheet (IIS), it is likely these ice sheets will have had a strong non-local influence on the spatial and temporal behaviour of the GrIS. Most previous paleo ice sheet modelling simulations recreated an ice sheet that either did not extend out onto the continental shelf; or utilized a simplified marine ice parametersiation and therefore did not fully include ice shelf dynamics, and or the sensitivity of the GrIS to this non-local signal from the surrounding ice sheets. In this paper, we investigated the evolution of the GrIS over the two most recent glacial-interglacial cycles (240 kyr BP to present day), using the ice sheet-ice shelf model, IMAU-ICE and investigated the influence of the LIS and IIS via an offline relative sea level (RSL) forcing generated by a GIA model. This RSL forcing controlled via changes in the water depth below the developing ice shelves, the spatial and temporal pattern of sub-ice shelf melting, which was parametrised in relation to changes in water depth. In the suite of simulations, the GrIS at the glacial maximums coalesced with the IIS to the north, expanded to the continental shelf break to the south west but remained too restricted to the north east. In terms of an ice-volume equivalent sea level contribution, at the Last Interglacial (LIG) and LGM the ice sheet added 1.46 m and −2.59 m to the budget respectively. The estimated lowering of the sea level by the Greenland contribution is considerably more (1.26 m) than most previous studies indicated whereas the contribution to the LIG high stand is lower (0.7 m). The spatial and temporal behaviour of the northern margin was highly variable in all simulations, controlled by the sub surface melt (SSM), which was dictated by the RSL forcing and the glacial history of the IIS and LIS. In contrast, the southwestern part of the ice sheet was insensitive to these forcing’s, with a uniform response in an all simulations controlled by the surface air temperature (SAT) forcing, derived from ice cores.


Sign in / Sign up

Export Citation Format

Share Document