Simulation of the last glacial cycle by a coupled, sectorially averaged climate—ice sheet model: 1. The climate model

1991 ◽  
Vol 96 (D7) ◽  
pp. 13139 ◽  
Author(s):  
H. Gallée ◽  
J. P. van Ypersele ◽  
Th. Fichefet ◽  
Ch. Tricot ◽  
A. Berger
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.


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.


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.


1997 ◽  
Vol 25 ◽  
pp. 58-65 ◽  
Author(s):  
L. Tarasov ◽  
W. R. Peltier

Significant improvements to the representation of climate forcing and mass-balance response in a coupled two-dimensional global energy balance climate model (EBM) and vertically integrated ice-sheet model (ISM) have led to the prediction of an ice-volume chronology for the most recent ice-age cycle of the Northern Hemisphere that is close to that inferred from the geological record. Most significant is that full glacial termination is delivered by the model without the need for new physical ingredients. In addition, a relatively close match is achieved between the Last Glacial Maximum (LGM) model ice topography and that of the recently-described ICE-4G reconstruction. These results suggest that large-scale climate system reorganization is not required to explain the main variations of the North American (NA) ice sheets over the last glacial cycle. Lack of sea-ice and marine-ice dynamics in the model leaves the situation over the Eurasian (EA) sector much more uncertain.The incorporation of a gravitationally self-consistent description of the glacial isostatic adjustment process demonstrates that the NA and EA bedrock responses can be adequately represented by simpler damped-relaxation models with characteristic time-scales of 3–5ka and 5 ka, respectively. These relaxation times agree with those independently inferred on the basis of postglacial relative sea-level histories.


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.


2015 ◽  
Vol 11 (9) ◽  
pp. 1165-1180 ◽  
Author(s):  
M. Willeit ◽  
A. Ganopolski

Abstract. Permafrost influences a number of processes which are relevant for local and global climate. For example, it is well known that permafrost plays an important role in global carbon and methane cycles. Less is known about the interaction between permafrost and ice sheets. In this study a permafrost module is included in the Earth system model CLIMBER-2, and the coupled Northern Hemisphere (NH) permafrost–ice-sheet evolution over the last glacial cycle is explored. The model performs generally well at reproducing present-day permafrost extent and thickness. Modeled permafrost thickness is sensitive to the values of ground porosity, thermal conductivity and geothermal heat flux. Permafrost extent at the Last Glacial Maximum (LGM) agrees well with reconstructions and previous modeling estimates. Present-day permafrost thickness is far from equilibrium over deep permafrost regions. Over central Siberia and the Arctic Archipelago permafrost is presently up to 200–500 m thicker than it would be at equilibrium. In these areas, present-day permafrost depth strongly depends on the past climate history and simulations indicate that deep permafrost has a memory of surface temperature variations going back to at least 800 ka. Over the last glacial cycle permafrost has a relatively modest impact on simulated NH ice sheet volume except at LGM, when including permafrost increases ice volume by about 15 m sea level equivalent in our model. This is explained by a delayed melting of the ice base from below by the geothermal heat flux when the ice sheet sits on a porous sediment layer and permafrost has to be melted first. Permafrost affects ice sheet dynamics only when ice extends over areas covered by thick sediments, which is the case at LGM.


Sign in / Sign up

Export Citation Format

Share Document