scholarly journals Investigating the evolution of major Northern Hemisphere ice sheets during the last glacial-interglacial cycle

2009 ◽  
Vol 5 (3) ◽  
pp. 329-345 ◽  
Author(s):  
S. Bonelli ◽  
S. Charbit ◽  
M. Kageyama ◽  
M.-N. Woillez ◽  
G. Ramstein ◽  
...  

Abstract. A 2.5-dimensional climate model of intermediate complexity, CLIMBER-2, fully coupled with the GREMLINS 3-D thermo-mechanical ice sheet model is used to simulate the evolution of major Northern Hemisphere ice sheets during the last glacial-interglacial cycle and to investigate the ice sheets responses to both insolation and atmospheric CO2 concentration. This model reproduces the main phases of advance and retreat of Northern Hemisphere ice sheets during the last glacial cycle, although the amplitude of these variations is less pronounced than those based on sea level reconstructions. At the last glacial maximum, the simulated ice volume is 52.5×1015 m3 and the spatial distribution of both the American and Eurasian ice complexes is in reasonable agreement with observations, with the exception of the marine parts of these former ice sheets. A set of sensitivity studies has also been performed to assess the sensitivity of the Northern Hemisphere ice sheets to both insolation and atmospheric CO2. Our results suggest that the decrease of summer insolation is the main factor responsible for the early build up of the North American ice sheet around 120 kyr BP, in agreement with benthic foraminifera δ18O signals. In contrast, low insolation and low atmospheric CO2 concentration are both necessary to trigger a long-lasting glaciation over Eurasia.

2009 ◽  
Vol 5 (2) ◽  
pp. 1013-1053 ◽  
Author(s):  
S. Bonelli ◽  
S. Charbit ◽  
M. Kageyama ◽  
M.-N. Woillez ◽  
G. Ramstein ◽  
...  

Abstract. A 2.5-dimensional climate model of intermediate complexity fully coupled with a 3-dimensional thermo-mechanical ice sheet model is used to simulate the evolution of major Northern Hemisphere ice sheets during the last glacial-interglacial cycle and to investigate the ice sheets responses to both insolation and atmospheric CO2 concentration. This model reproduces the main phases of advance and retreat of Northern Hemisphere ice sheets during the last glacial cycle, although the amplitude of these variations is less pronounced than those based on sea level reconstructions. At the last glacial maximum, the simulated ice volume is 52.5×1015 m3 and the spatial distribution of both the American and Eurasian ice complexes is in reasonable agreement with observations, with the exception of the marine parts of these former ice sheets. A set of sensitivity studies has also been performed to assess the sensitivity of the Northern Hemisphere ice sheets to both insolation and atmospheric CO2. Our results suggest that the decrease of summer insolation is the main factor responsible for the early build up of the North American ice sheet around 120 kyr BP, in agreement with benthic foraminifera δ18O signals. In contrast, low insolation and low atmospheric CO2 concentration are both necessary to trigger a long-lasting glaciation over Eurasia.


2012 ◽  
Vol 6 (6) ◽  
pp. 4897-4938 ◽  
Author(s):  
S. Charbit ◽  
C. Dumas ◽  
M. Kageyama ◽  
D. M. Roche ◽  
C. Ritz

Abstract. Since the original formulation of the positive-degree-day (PDD) method, different PDD calibrations have been proposed in the literature in response to the increasing number of observations. Although these formulations provide a satisfactory description of the present-day Greenland geometry, they have not all been tested for paleo ice sheets. Using the climate-ice sheet model CLIMBER-GRISLI coupled with different PDD models, we evaluate how the parameterization of the ablation may affect the evolution of Northern Hemisphere ice sheets in the transient simulations of the last glacial cycle. Results from fully coupled simulations are compared to time-slice experiments carried out at different key periods of the last glacial period. We find large differences in the simulated ice sheets according to the chosen PDD model. These differences occur as soon as the onset of glaciation, therefore affecting the subsequent evolution of the ice system. To further investigate how the PDD method controls this evolution, special attention is given to the role of each PDD parameter. We show that glacial inception is critically dependent on the representation of the impact of the temperature variability from the daily to the inter-annual time scale, whose effect is modulated by the refreezing scheme. Finally, an additional set of sensitivity experiments has been carried out to assess the relative importance of melt processes with respect to initial ice sheet configuration in the construction and the evolution of past Northern Hemisphere ice sheets. Our analysis reveals that the impacts of the initial ice sheet condition may range from quite negligible to explaining about half of the LGM ice volume depending on the representation of stochastic temperature variations which remain the main driver of the evolution of the ice system.


2019 ◽  
Vol 65 (252) ◽  
pp. 645-661 ◽  
Author(s):  
LU NIU ◽  
GERRIT LOHMANN ◽  
SEBASTIAN HINCK ◽  
EVAN J. GOWAN ◽  
UTA KREBS-KANZOW

ABSTRACTThe evolution of Northern Hemisphere ice sheets through the last glacial cycle is simulated with the glacial index method by using the climate forcing from one General Circulation Model, COSMOS. By comparing the simulated results to geological reconstructions, we first show that the modelled climate is capable of capturing the main features of the ice-sheet evolution. However, large deviations exist, likely due to the absence of nonlinear interactions between ice sheet and other climate components. The model uncertainties of the climate forcing are examined using the output from nine climate models from the Paleoclimate Modelling Intercomparison Project Phase III. The results show a large variability in simulated ice sheets between the different models. We find that the ice-sheet extent pattern resembles summer surface air temperature pattern at the Last Glacial Maximum, confirming the dominant role of surface ablation process for high-latitude Northern Hemisphere ice sheets. This study shows the importance of the upper boundary condition for ice-sheet modelling, and implies that careful constraints on climate output is essential for simulating realistic glacial Northern Hemisphere ice sheets.


1997 ◽  
Vol 25 ◽  
pp. 333-339 ◽  
Author(s):  
Philippe Huybrechts ◽  
Stephen T’siobbel

A quasi-three-dimensional (3-D) climate model (Sellers, 1983) was used to simulate the climate of the Last Glacial Maximum (LGM) in order to provide climatic input for the modelling of the Northern Hemisphere ice sheets. The climate model is basically a coarse-gridded general circulation (GCM) with simplified dynamics, and was subject to appropriate boundary conditions for ice-sheet elevation, atmospheric CO2concentration and orbital parameters. When compared with the present-daysimulation, the simulated climate at the Last Glacial Maximum is characterized by a global annual cooling of 3.5°C and a reduction in global annualprecipitation of 7.5%, which agrees well with results from other, more complex GCMs. Also the patterns of temperature change compare fairly with mostother GCM results, except for a smaller cooling over the North Atlantic and the larger cooling predicted for the summer rather than for the winter over Eurasia.The climate model is able to simulate changes in Northern Hemisphere tropospheric circulation, yielding enhanced westerlies in the vicinity of the Laurentide and Eurasian ice sheets. However, the simulated precipitation patterns are less convincing, and show a distinct mean precipitation increase over the Laurentide ice sheet. Nevertheless, when using the mean-monthly fields of LGM minus present-day anomalies of temperature and precipitation rate to drive a three-dimensional thermomechanical ice-sheet model, it was demonstrated that within realistic bounds of the ice-flow and mass-balance parameters, veryreasonable reconstructions of the Last Glacial Maximum ice sheets could be obtained.


2020 ◽  
Author(s):  
Kenji Izumi ◽  
Paul Paul Valdes ◽  
Ruza Ivanovic ◽  
Lauren Gregoire

<p>The Last Glacial Maximum (LGM; 21,000 yr before present) is a target period of the paleoclimate simulations in the Coupled Model Intercomparison Project Phase 6 – the Paleoclimate Modeling Intercomparison Project Phase 4 (CMIP6-PMIP4) because of abundant paleoenvironmental data in continental, ice, and marine indicators. The LGM was a period of low atmospheric trace gases when large ice sheets covered over North America and Scandinavia. Paleoclimate reconstructions and modeling studies suggest that the Northern Hemisphere climate differed from today.</p><p>In this study, we used the coupled atmosphere and ocean model HadCM3B-M1 in order to investigate the impacts of the main LGM boundary condition changes, in particular, the ICE-6G_C, GLAC-1D, and PMIP3 ice-sheet reconstructions following the PMIP4 protocol, on the mean state of the climate over the Northern Hemisphere. First, we check the surface albedo forcing and feedback with a simplified partial derivative method and assess the surface temperature changes and their composition using a simple surface energy balance equation. Then, we investigate how patterns of stationary waves, westerly jet precipitation over the Northern Hemisphere change in response to the LGM ice-sheet configuration. Finally, we implement a paleo data-model comparison for validation of the large-scale climate changes over the Northern Hemisphere at the LGM. The wintertime stationary waves have the largest amplitude and different responses among the experiments, while stationary waves in summer are weak and similar responses. The LGM simulation with the ICE-6G_C better captures features of the LGM climate, but compared to the reconstructions, the climate model tends to overestimate cooling in summer and underestimate cooling in winter and simulate wetter conditions over the Northern Hemisphere.  </p>


2007 ◽  
Vol 3 (3) ◽  
pp. 423-438 ◽  
Author(s):  
A. Abe-Ouchi ◽  
T. Segawa ◽  
F. Saito

Abstract. The ice sheet-climate interaction as well as the climatic response to orbital parameters and atmospheric CO2 concentration are examined in order to drive an ice sheet model throughout an ice age cycle. Feedback processes between ice sheet and atmosphere are analyzed by numerical experiments using a high resolution General Circulation Model (GCM) under different conditions at the Last Glacial Maximum. Among the proposed processes, the ice albedo feedback, the elevation-mass balance feedback and the desertification effect over the ice sheet were found to be the dominant processes for the ice-sheet mass balance. For the elevation-mass balance feedback, the temperature lapse rate over the ice sheet is proposed to be weaker than assumed in previous studies. Within the plausible range of parameters related to these processes, the ice sheet response to the orbital parameters and atmospheric CO2 concentration for the last glacial/interglacial cycle was simulated in terms of both ice volume and geographical distribution, using a three-dimensional ice-sheet model. Careful treatment of climate-ice sheet feedback is essential for a reliable simulation of the ice sheet changes during ice age cycles.


1997 ◽  
Vol 25 ◽  
pp. 333-339 ◽  
Author(s):  
Philippe Huybrechts ◽  
Stephen T’siobbel

A quasi-three-dimensional (3-D) climate model (Sellers, 1983) was used to simulate the climate of the Last Glacial Maximum (LGM) in order to provide climatic input for the modelling of the Northern Hemisphere ice sheets. The climate model is basically a coarse-gridded general circulation (GCM) with simplified dynamics, and was subject to appropriate boundary conditions for ice-sheet elevation, atmospheric CO2 concentration and orbital parameters. When compared with the present-daysimulation, the simulated climate at the Last Glacial Maximum is characterized by a global annual cooling of 3.5°C and a reduction in global annualprecipitation of 7.5%, which agrees well with results from other, more complex GCMs. Also the patterns of temperature change compare fairly with mostother GCM results, except for a smaller cooling over the North Atlantic and the larger cooling predicted for the summer rather than for the winter over Eurasia.The climate model is able to simulate changes in Northern Hemisphere tropospheric circulation, yielding enhanced westerlies in the vicinity of the Laurentide and Eurasian ice sheets. However, the simulated precipitation patterns are less convincing, and show a distinct mean precipitation increase over the Laurentide ice sheet. Nevertheless, when using the mean-monthly fields of LGM minus present-day anomalies of temperature and precipitation rate to drive a three-dimensional thermomechanical ice-sheet model, it was demonstrated that within realistic bounds of the ice-flow and mass-balance parameters, veryreasonable reconstructions of the Last Glacial Maximum ice sheets could be obtained.


2003 ◽  
Vol 37 ◽  
pp. 173-180 ◽  
Author(s):  
Chris Zweck ◽  
Philippe Huybrechts

AbstractMechanisms that determine time-dependent changes of the marine ice margin in dynamic ice-sheet models are important but poorly understood. Here we derive an empirical formulation for changes in the marine extent when modelling the Northern Hemisphere ice sheets over the last glacial cycle in a three-dimensional thermomechanically coupled ice-sheet model. We assume that the strongest control on changes in marine extent is ice calving, and that the variable most crucial to calving is water depth. The empirical marine-extent relationship is tuned so that the major marine-retreat history of the Laurentide and Eurasian ice sheets is modelled accurately in time and space. We find that this empirical treatment relating marine extent to water depth is sufficient to reproduce the observations, and discuss the implications for the physics of marine margin changes and the dynamics of the Northern Hemisphere ice sheets since the Last Glacial Maximum.


2013 ◽  
Vol 7 (2) ◽  
pp. 681-698 ◽  
Author(s):  
S. Charbit ◽  
C. Dumas ◽  
M. Kageyama ◽  
D. M. Roche ◽  
C. Ritz

Abstract. Since the original formulation of the positive-degree-day (PDD) method, different PDD calibrations have been proposed in the literature in response to the increasing number of observations. Although these formulations generally provide a satisfactory description of the present-day Greenland geometry, they have not all been tested for paleo ice sheets. Using the climate-ice sheet model CLIMBER-GRISLI coupled with different PDD models, we evaluate how the parameterisation of the ablation may affect the evolution of Northern Hemisphere ice sheets in the transient simulations of the last glacial cycle. Results from fully coupled simulations are compared to time-slice experiments carried out at different key periods of the last glacial period. We find large differences in the simulated ice sheets according to the chosen PDD model. These differences occur as soon as the onset of glaciation, therefore affecting the subsequent evolution of the ice system. To further investigate how the PDD method controls this evolution, special attention is given to the role of each PDD parameter. We show that glacial inception is critically dependent on the representation of the impact of the temperature variability from the daily to the inter-annual time scale, whose effect is modulated by the refreezing scheme. Finally, an additional set of sensitivity experiments has been carried out to assess the relative importance of melt processes with respect to initial ice sheet configuration in the construction and the evolution of past Northern Hemisphere ice sheets. Our analysis reveals that the impacts of the initial ice sheet condition may range from quite negligible to explaining about half of the LGM ice volume depending on the representation of stochastic temperature variations which remain the main driver of the evolution of the ice system. The main findings of this paper underline the need for conducting studies with high resolution climate models coupled to detailed snow models to better constrain the temporal and spatial variations of the PDD parameters. The development of such approaches could improve the calibration of the PDD formulation which is still widely used in climate-ice sheet studies.


2014 ◽  
Vol 10 (4) ◽  
pp. 1489-1500 ◽  
Author(s):  
N. Korhonen ◽  
A. Venäläinen ◽  
H. Seppä ◽  
H. Järvinen

Abstract. Earth system models of intermediate complexity (EMICs) have proven to be able to simulate the large-scale features of glacial–interglacial climate evolution. For many climatic applications the spatial resolution of the EMICs' output is, however, too coarse, and downscaling methods are needed. In this study we introduce a way to use generalized additive models (GAMs) for downscaling the large-scale output of an EMIC in very different climatological conditions ranging from glacial periods to current relatively warm climates. GAMs are regression models in which a combination of explanatory variables is related to the response through a sum of spline functions. We calibrated the GAMs using observations of the recent past climate and the results of short time-slice simulations of glacial climate performed by the relatively high-resolution general circulation model CCSM (Community Climate System Model) and the regional climate model RCA3 (Rossby Centre regional Atmospheric climate model). As explanatory variables we used the output of a simulation by the CLIMBER-2 (CLIMate and BiosphERe model 2) EMIC of the last glacial cycle, coupled with the SICOPOLIS (SImulation COde for POLythermal Ice Sheets) ice sheet model, i.e. the large-scale temperature and precipitation data of CLIMBER-2, and the elevation, distance to ice sheet, slope direction and slope angle from SICOPOLIS. The fitted GAMs were able to explain more than 96% of the temperature response with a correlation of >0.98 and more than 59% of the precipitation response with a correlation of >0.72. The first comparison with two pollen-based reconstructions of temperature for Northern Europe showed that CLIMBER-2 data downscaled by GAMs corresponded better with the reconstructions than did the bilinearly interpolated CLIMBER-2 surface temperature.


Sign in / Sign up

Export Citation Format

Share Document