scholarly journals A new approach for simulating the paleo evolution of the Northern Hemisphere ice sheets

2017 ◽  
Author(s):  
Rubén Banderas ◽  
Jorge Alvarez-Solas ◽  
Alexander Robinson ◽  
Marisa Montoya
Keyword(s):  
Climate Variability ◽  
Northern Hemisphere ◽  
Ice Sheets ◽  
Ice Sheet ◽  
Last Glacial Period ◽  
New Approach ◽  
Anomaly Field ◽  
Ice Volume ◽  
The Last Glacial ◽  
Millennial Scale

Abstract. The aim of this study is to assess and improve the methods currently used to force ice sheet models offline. To this end, three different synthetic transient forcing climatologies are developed for the past 120 kyr following a perturbative approach and applied to an ice-sheet model. The results are used to evaluate their consequences for simulating the paleo evolution of the Northern Hemisphere ice sheets. The first method follows the usual approach in which temperature anomalies relative to present are calculated by combining a present-day climatology with a simulated glacial-interglacial climatic anomaly field interpolated through an index derived from ice-core data. In the second approach the representation of millennial-scale climate variability is improved by incorporating a simulated stadial-interstadial anomaly field. The third is a refinement of the second one in which the amplitudes of both orbital and millennial-scale variations are corrected to provide a perfect agreement with a recent absolute temperature reconstruction over Greenland. The comparison of the three climate forcing methods highlights the tendency of the usual approach to overestimate the temperature variability over North America and Eurasia at millennial timescales. This leads to a relatively high Northern Hemisphere (NH) ice-volume variability on these timescales. Through enhanced ablation, this results in too low an ice volume throughout the last glacial period (LGP), below or at the lower end of the uncertainty range of estimations. Improving the representation of millennial-scale variability alone yields an important increase of ice volume in all NH ice sheets, but especially in the Fennoscandian ice sheet (FIS). Optimizing the amplitude of the temperature anomalies to match the Greenland reconstruction results in a further increase of the simulated ice-sheet volume throughout the LGP. Our new method provides a more realistic representation of orbital and millennial scale climate variability and represents an improvement in the transient forcing of ice sheets during the last glacial period. Interestingly, our new approach underestimates ice-volume variations on millennial timescales as indicated by sea-level records. This suggests that either the origin of the latter is not the NH or that processes not represented in our study, notably variations in oceanic conditions, need to be invoked to account for an important role of millennial-scale climate variability on millennial-scale ice- volume fluctuations. We finally provide here both our derived climate evolution of the LGP using the three methods as well as the resulting ice-sheet configurations. These could be of interest for future studies dealing with the atmospheric or/and oceanic consequences of transient ice-sheet evolution throughout the LGP, and as a source of climate input to other ice sheet models.

scholarly journals A new approach for simulating the paleo-evolution of the Northern Hemisphere ice sheets

2018 ◽  
Vol 11 (6) ◽  
pp. 2299-2314 ◽  
Author(s):  
Rubén Banderas ◽  
Jorge Alvarez-Solas ◽  
Alexander Robinson ◽  
Marisa Montoya
Keyword(s):  
Climate Variability ◽  
Ice Core ◽  
Ice Sheets ◽  
Ice Sheet ◽  
Temperature Reconstruction ◽  
Anomaly Field ◽  
Temperature Anomalies ◽  
Ice Volume ◽  
Index Approach ◽  
Millennial Scale

Abstract. Offline forcing methods for ice-sheet models often make use of an index approach in which temperature anomalies relative to the present are calculated by combining a simulated glacial–interglacial climatic anomaly field, interpolated through an index derived from the Greenland ice-core temperature reconstruction, with present-day climatologies. An important drawback of this approach is that it clearly misrepresents climate variability at millennial timescales. The reason for this is that the spatial glacial–interglacial anomaly field used is associated with orbital climatic variations, while it is scaled following the characteristic time evolution of the index, which includes orbital and millennial-scale climate variability. The spatial patterns of orbital and millennial variability are clearly not the same, as indicated by a wealth of models and data. As a result, this method can be expected to lead to a misrepresentation of climate variability and thus of the past evolution of Northern Hemisphere (NH) ice sheets. Here we illustrate the problems derived from this approach and propose a new offline climate forcing method that attempts to better represent the characteristic pattern of millennial-scale climate variability by including an additional spatial anomaly field associated with this timescale. To this end, three different synthetic transient forcing climatologies are developed for the past 120 kyr following a perturbative approach and are applied to an ice-sheet model. The impact of the climatologies on the paleo-evolution of the NH ice sheets is evaluated. The first method follows the usual index approach in which temperature anomalies relative to the present are calculated by combining a simulated glacial–interglacial climatic anomaly field, interpolated through an index derived from ice-core data, with present-day climatologies. In the second approach the representation of millennial-scale climate variability is improved by incorporating a simulated stadial–interstadial anomaly field. The third is a refinement of the second one in which the amplitudes of both orbital and millennial-scale variations are tuned to provide perfect agreement with a recently published absolute temperature reconstruction over Greenland. The comparison of the three climate forcing methods highlights the tendency of the usual index approach to overestimate the temperature variability over North America and Eurasia at millennial timescales. This leads to a relatively high NH ice-volume variability on these timescales. Through enhanced ablation, this results in too low an ice volume throughout the last glacial period (LGP), below or at the lower end of the uncertainty range of estimations. Improving the representation of millennial-scale variability alone yields an important increase in ice volume in all NH ice sheets but especially in the Fennoscandian Ice Sheet (FIS). Optimizing the amplitude of the temperature anomalies to match the Greenland reconstruction results in a further increase in the simulated ice-sheet volume throughout the LGP. Our new method provides a more realistic representation of orbital and millennial-scale climate variability and improves the transient forcing of ice sheets during the LGP. Interestingly, our new approach underestimates ice-volume variations on millennial timescales as indicated by sea-level records. This suggests that either the origin of the latter is not the NH or that processes not represented in our study, notably variations in oceanic conditions, need to be invoked to explain millennial-scale ice-volume fluctuations. We finally provide here both our derived climate evolution of the LGP using the three methods as well as the resulting ice-sheet configurations. These could be of interest for future studies dealing with the atmospheric or/and oceanic consequences of transient ice-sheet evolution throughout the LGP and as a source of climate input to other ice-sheet models.

scholarly journals Ocean-driven millennial-scale variability of the Eurasian ice sheet during the last glacial period simulated with a hybrid ice-sheet–shelf model

2019 ◽  
Vol 15 (3) ◽  
pp. 957-979 ◽  
Author(s):  
Jorge Alvarez-Solas ◽  
Rubén Banderas ◽  
Alexander Robinson ◽  
Marisa Montoya
Keyword(s):  
Climate Variability ◽  
Climate Changes ◽  
Ice Sheets ◽  
Ice Sheet ◽  
Glacial Period ◽  
Last Glacial Period ◽  
Perturbative Approach ◽  
Last Glacial ◽  
The Last Glacial ◽  
Millennial Scale

Abstract. The last glacial period (LGP; ca. 110–10 kyr BP) was marked by the existence of two types of abrupt climatic changes, Dansgaard–Oeschger (DO) and Heinrich (H) events. Although the mechanisms behind these are not fully understood, it is generally accepted that the presence of ice sheets played an important role in their occurrence. While an important effort has been made to investigate the dynamics and evolution of the Laurentide ice sheet (LIS) during this period, the Eurasian ice sheet (EIS) has not received much attention, in particular from a modeling perspective. However, meltwater discharge from this and other ice sheets surrounding the Nordic seas is often implied as a potential cause of ocean instabilities that lead to glacial abrupt climate changes. Thus, a better comprehension of the evolution of the EIS during the LGP is important to understand its role in glacial abrupt climate changes. Here we investigate the response of the EIS to millennial-scale climate variability during the LGP. We use a hybrid, three-dimensional, thermomechanical ice-sheet model that includes ice shelves and ice streams. The model is forced off-line via a novel perturbative approach that, as opposed to conventional methods, clearly differentiates between the spatial patterns of millennial-scale and orbital-scale climate variability. Thus, it provides a more realistic treatment of the forcing at millennial timescales. The effect of both atmospheric and oceanic variations are included. Our results show that the EIS responds with enhanced ice discharge in phase with interstadial warming in the North Atlantic when forced with surface ocean temperatures. Conversely, when subsurface ocean temperatures are used, enhanced ice discharge occurs both during stadials and at the beginning of the interstadials. Separating the atmospheric and oceanic effects demonstrates the major role of the ocean in controlling the dynamics of the EIS on millennial timescales. While the atmospheric forcing alone is only able to produce modest iceberg discharges, warming of the ocean leads to higher rates of iceberg discharges as a result of relatively strong basal melting at the margins of the ice sheet. Our results clearly show the capability of the EIS to react to glacial abrupt climate changes, and highlight the need for stronger constraints on the ice sheet's glacial dynamics and climate–ocean interactions.

scholarly journals Oceanic forcing of the Eurasian Ice Sheet on millennial time scales during the Last Glacial Period

2018 ◽  
Author(s):  
Jorge Alvarez-Solas ◽  
Rubén Banderas ◽  
Alexander Robinson ◽  
Marisa Montoya
Keyword(s):  
Time Scales ◽  
North Atlantic ◽  
Climate Changes ◽  
Ice Sheets ◽  
Ice Sheet ◽  
Glacial Period ◽  
Last Glacial Period ◽  
Last Glacial ◽  
The Last Glacial ◽  
Millennial Scale

Abstract. The last glacial period (LGP; ca.110–10 ka BP) was marked by the existence of two types of abrupt climatic changes, Dansgaard-Oeschger (DO) and Heinrich (H) events. Although the mechanisms behind these are not fully understood, it is generally accepted that the presence of ice sheets played an important role in their occurrence. While an important effort has been made to investigate the dynamics and evolution of the Laurentide Ice Sheet (LIS) during this period, the Eurasian Ice Sheet (EIS) has not received much attention, in particular from a modeling perspective. However, meltwater discharge from this and other ice sheets surrounding the Nordic Seas is often implied as a potential cause of ocean instabilities that lead to glacial abrupt climate changes. Thus, a better understanding of its variations during the LGP is important to understand its role in glacial abrupt climate changes. Here we investigate the response of the EIS to millennial-scale climate variability during the LGP. We use a hybrid, three-dimensional, thermomechanical ice-sheet model that includes ice shelves and ice streams. The model is forced offline through a novel perturbative approach that includes the effect of both atmospheric and oceanic variations and provides a more realistic treatment of millennial-scale climatic variability than conventional methods. Our results show that the EIS responds with enhanced ice discharge in phase with interstadial warming in the North Atlantic when forced with surface ocean temperatures. Conversely, when subsurface ocean temperatures are used, enhanced ice discharge occurs both during stadials and at the beginning of the interstadials. Separating the atmospheric and oceanic effects demonstrates the major role of the ocean in controlling the dynamics of the EIS on millennial time scales. While the atmospheric forcing alone is only able to produce modest iceberg discharges, warming of the ocean leads to higher rates of iceberg discharges as a result of relatively strong basal melting at the margins of the ice sheet. Together with previous work, our results provide a consistent explanation for the response of the LIS and the EIS to glacial abrupt climate changes, and highlight the need for stronger constraints on the local North Atlantic behavior in order to improve our understanding of the ice sheet's glacial dynamics.

scholarly journals Numerical reconstructions of the Northern Hemisphere ice sheets through the last glacial-interglacial cycle

2007 ◽  
Vol 3 (1) ◽  
pp. 15-37 ◽  
Author(s):  
S. Charbit ◽  
C. Ritz ◽  
G. Philippon ◽  
V. Peyaud ◽  
M. Kageyama
Keyword(s):  
Northern Hemisphere ◽  
General Circulation ◽  
Ice Sheets ◽  
Ice Sheet ◽  
Maximum Amplitude ◽  
The Other ◽  
Climate Forcing ◽  
Last Glacial ◽  
Ice Volume ◽  
The Last Glacial

Abstract. A 3-dimensional thermo-mechanical ice-sheet model is used to simulate the evolution of the Northern Hemisphere ice sheets through the last glacial-interglacial cycle. The ice-sheet model is forced by the results from six different atmospheric general circulation models (AGCMs). The climate evolution over the period under study is reconstructed using two climate equilibrium simulations performed for the Last Glacial Maximum (LGM) and for the present-day periods and an interpolation through time between these snapshots using a glacial index calibrated against the GRIP δ18O record. Since it is driven by the timing of the GRIP signal, the temporal evolution of the ice volume and the ice-covered area is approximately the same from one simulation to the other. However, both ice volume curves and spatial distributions of the ice sheets present some major differences from one AGCM forcing to the other. The origin of these differences, which are most visible in the maximum amplitude of the ice volume, is analyzed in terms of differences in climate forcing. This analysis allows for a partial evaluation of the ability of GCMs to simulate climates consistent with the reconstructions of past ice sheets. Although some models properly reproduce the advance or retreat of ice sheets in some specific areas, none of them is able to reproduce both North American or Eurasian ice complexes in full agreement with observed sea-level variations and geological data. These deviations can be attributed to shortcomings in the climate forcing and in the LGM ice-sheet reconstruction used as a boundary condition for GCM runs, but also to missing processes in the ice-sheet model itself.

Water vapour, CO 2 and insolation over the last glacial-interglacial cycles

1993 ◽  
Vol 341 (1297) ◽  
pp. 253-261 ◽  
Keyword(s):  
Water Vapour ◽  
Northern Hemisphere ◽  
Ice Sheets ◽  
Ice Sheet ◽  
Last Glacial ◽  
Ice Volume ◽  
Water Vapour Feedback ◽  
Astronomical Forcing ◽  
The Last Glacial Maximum ◽  
The Last Glacial

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.

scholarly journals Numerical reconstructions of the penultimate glacial maximum Northern Hemisphere ice sheets: sensitivity to climate forcing and model parameters

2016 ◽  
Vol 62 (234) ◽  
pp. 607-622 ◽  
Author(s):  
CLAUDIA WEKERLE ◽  
FLORENCE COLLEONI ◽  
JENS-OVE NÄSLUND ◽  
JENNY BRANDEFELT ◽  
SIMONA MASINA
Keyword(s):  
Northern Hemisphere ◽  
Ice Sheets ◽  
Ice Sheet ◽  
Glacial Maximum ◽  
Laurentide Ice Sheet ◽  
Model Parameters ◽  
Geological Evidence ◽  
Last Glacial ◽  
Ice Volume ◽  
The Last Glacial

ABSTRACTNumerous ice-sheet reconstructions of the last glacial cycle have been proposed, however due to limited geological evidence, reconstructing older Northern Hemisphere ice sheets remains a difficult exercise. Here we focus on the penultimate glacial maximum (PGM; ~140 ka BP) over the Northern Hemisphere. While some evidence of the PGM Eurasian ice-sheet extent were found, this is not the case for the corresponding Laurentide ice sheet. To improve the glaciological reconstructions of the PGM Northern Hemisphere ice sheets, we explore the parameter space of ice-sheet model uncertainties and carry out numerous univariate ice-sheet steady-state sensitivity simulations. We use two PGM climate simulations to force the ice-sheet model, differing in the prescribed Laurentide ice topography (small and large). The simulated Northern Hemisphere ice volume ranges from 124.7 to 152 m SLE when using the climate accounting for a small Laurentide ice sheet, which is compatible with global sea-level reconstructions of this period (−92 to −150 m). Conversely, using the climate simulation with a Laurentide ice sheet comparable in size to that of the last glacial maximum results in too large ice volumes. Changes in basal drag provide the upper bound ice volume of our experiments, whereas changes in the distribution of ice streams provide the lower bound.

scholarly journals Numerical reconstructions of the Northern Hemisphere ice sheets through the last glacial-interglacial cycle

2006 ◽  
Vol 2 (5) ◽  
pp. 879-921 ◽  
Author(s):  
S. Charbit ◽  
C. Ritz ◽  
G. Philippon ◽  
V. Peyaud ◽  
M. Kageyama
Keyword(s):  
Northern Hemisphere ◽  
General Circulation ◽  
Ice Sheets ◽  
Ice Sheet ◽  
Maximum Amplitude ◽  
The Other ◽  
Climate Forcing ◽  
Last Glacial ◽  
Ice Volume ◽  
The Last Glacial

Abstract. A 3-dimensional thermo-mechanical ice-sheet model is used to simulate the evolution of the Northern hemisphere ice sheets through the last glacial-interglacial cycle. The ice-sheet model is forced by the results from six different atmospheric general circulation models (AGCMs). Two climate snapshots simulations performed for the last glacial maximum (LGM) and for the present-day periods are interpolated through time using a glacial index calibrated against the GRIP δ18O record to reconstruct the climate evolution over the period under study. Since it is driven by the timing of the GRIP signal, the temporal evolution of the ice volume and the ice-covered area is approximately the same from one simulation to the other. However, both ice volume curves and spatial distributions of the ice sheets present some major differences from one AGCM to the other. The origin of these differences, which are most visible in the maximum amplitude of the ice volume, is analyzed in terms of differences in climate forcing. The analysis of the results allows an evaluation of the ability of GCMs to simulate climates consistent with the reconstructions of past ice sheets to be evaluated. Although some models properly reproduce the advance or retreat of ice sheets in some specific areas, none of them is able to reproduce both North American or Eurasian ice complexes in full agreement with observed sea-level variations and geological data. These deviations can be attributed to shortcomings in the climate forcing and in the LGM ice-sheet reconstruction used as a boundary condition for GCM runs, but also to missing processes in the ice-sheet model itself.

scholarly journals Oceanic forcing of the Eurasian Ice Sheet on millennial time scales during the Last Glacial Period

2017 ◽  
Author(s):  
Jorge Alvarez-Solas ◽  
Rubén Banderas ◽  
Alexander Robinson ◽  
Marisa Montoya
Keyword(s):  
Time Scales ◽  
Climate Changes ◽  
Ice Sheets ◽  
Ice Sheet ◽  
Glacial Period ◽  
Last Glacial Period ◽  
Perturbative Approach ◽  
Last Glacial ◽  
The Last Glacial ◽  
Millennial Scale

Abstract. The last glacial period (LGP; ca. 110–10 ka BP) was marked by the existence of two types of abrupt climatic changes, Dansgaard-Oeschger (D/O) and Heinrich (H) events. Although the mechanisms behind these are not fully understood, it is generally accepted that the presence of ice sheets played an important role in their occurrence. While an important effort has been made to investigate the dynamics and evolution of the Laurentide Ice Sheet (LIS) during this period, the Eurasian Ice Sheet (EIS) has not received much attention, in particular from a modeling perspective. However, meltwater discharge from this and other ice sheets surrounding the Nordic Seas is often implied as a potential cause of ocean instabilities that lead to glacial abrupt climate changes. Thus, a better understanding of its variations during the LGP is important to understand its role in glacial abrupt climate changes. Here we investigate the response of the EIS to millennial-scale climate variability during the LGP. We use a hybrid, three-dimensional, thermomechanical ice-sheet model that includes ice shelves and ice streams. The model is forced offline through a novel perturbative approach that includes the effect of both atmospheric and oceanic variations and provides a more realistic treatment of millennial-scale climatic variability than conventional methods. Our results show that the EIS responds with enhanced iceberg discharges in phase with interstadial warming in the North Atlantic. Separating the atmospheric and oceanic effects demonstrates the major role of the ocean in controlling the dynamics of the EIS on millennial time scales. While the atmospheric forcing alone is only able to produce modest iceberg discharges, warming of oceanic surface waters leads to much higher rates of iceberg discharges as a result of relatively strong basal melting within the margins of the ice sheet. Together with previous work, our results provide a consistent explanation for the timing of the responses of the LIS and the EIS to glacial abrupt climate changes.

scholarly journals Simulations of the Northern Hemisphere through the last Glacial-interglacial cycle with a vertically integrated and a three-dimensional thermomechanical ice-sheet model coupled to a climate model

1998 ◽  
Vol 27 ◽  
pp. 169-176 ◽  
Author(s):  
R. Calov ◽  
I. Marsiat
Keyword(s):  
Northern Hemisphere ◽  
Enhancement Factor ◽  
Climate Model ◽  
Three Dimensional ◽  
Ice Sheets ◽  
Ice Sheet ◽  
Two Dimensional ◽  
Thermomechanical Model ◽  
Ice Volume ◽  
The Last Glacial

We present simulations of the Northern Hemisphere land ice through the last glacial-interglacial cycle with a vertically integrated ice-sheet model and a three-dimensional thermomechanical ice-sheet model. Both models are coupled asynchronously to the zonally averaged Louvain-la-Neuve climate model, which includes simplified treatments of the atmosphere, ocean and sea ice. The two-dimensional vertically integrated ice-sheet model, which contains no thermomechanical coupling, was developed in spherical coordinates (Marsiat, 1994). The three-dimensional thermomechanical ice-sheet model was developed using the two-dimensional vertically integrated model as source. We compare results of the vertically integrated with those of the thermomechanical ice-sheet model. in the thermomechanical model the deformation properties of ice depend on the temperature within the ice and the enhancement factor; the latter is introduced to model, in a simplified approach, the different flow properties of Pleistocene and Holocene ice due to varying dust content. The computations with the thermomechanical model show that the growth and decay of the Northern Hemisphere ice sheets can be modelled with a common enhancement factor for all ice sheets. It is shown that there are model set-ups for the thermomechanical model yielding temporal developments of the total ice volume comparable to those of the vertically integrated model. Furthermore, we demonstrate that for the coupled climate/cryosphere system the total ice volume depends considerably on the enhancement factor.

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 ◽  
...  
Keyword(s):  
Northern Hemisphere ◽  
Climate Model ◽  
Ice Sheets ◽  
Ice Sheet ◽  
Co2 Concentration ◽  
Atmospheric Co2 ◽  
Glacial Cycle ◽  
Last Glacial ◽  
Atmospheric Co2 Concentration ◽  
The Last Glacial

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.

Sign in / Sign up

Export Citation Format

Share Document