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

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.

2017 ◽  
Author(s):  
Rubén Banderas ◽  
Jorge Alvarez-Solas ◽  
Alexander Robinson ◽  
Marisa Montoya

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.


2019 ◽  
Vol 15 (3) ◽  
pp. 957-979 ◽  
Author(s):  
Jorge Alvarez-Solas ◽  
Rubén Banderas ◽  
Alexander Robinson ◽  
Marisa Montoya

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.


2016 ◽  
Author(s):  
David A. Hodell ◽  
James E.T. Channell

Abstract. We present a 3.2-Myr record of stable isotopes and physical properties at IODP Site U1308 (re-occupation of DSDP Site 609) located within the ice-rafted detritus (IRD) belt of the North Atlantic. We compare the isotope and lithological proxies at Site U1308 with other North Atlantic records (e.g., Sites 982, 607/U1313 and U1304) to reconstruct the history of orbital and millennial-scale climate variability during the Quaternary. The Site U1308 record documents a progressive increase in the intensity of Northern Hemisphere glacial-interglacial cycles during the late Pliocene and Quaternary with mode transitions at ~ 2.7, 1.5, 0.9 and 0.65 Ma. These transitions mark times of change in the growth and stability of Northern Hemisphere ice sheets. They also coincide with increases in vertical carbon isotope gradients between the intermediate and deep ocean, suggesting changes in deep carbon storage and atmospheric CO2. Orbital and millennial climate variability co-evolved during the Quaternary such that the trend towards larger ice sheets was accompanied by changes in the style, frequency and intensity of millennial-scale variability. This co-evolution may be important for explaining the observed patterns of Quaternary climate change.


2011 ◽  
Vol 57 (206) ◽  
pp. 1017-1026 ◽  
Author(s):  
Tsutomu Uchida ◽  
Atsushi Miyamoto ◽  
Atsushi Shin’yama ◽  
Takeo Hondoh

AbstractAir-hydrate crystals store most of the ancient air contained in deep ice sheets. We carried out microscopic observations of air-hydrate crystals below 2000 m depth within the ice core from Dome Fuji, Antarctica, to obtain their number and size distributions. We found that the number density continuously decreased with depth, whereas the average size increased, in contrast to findings from shallower depths. In addition, the characteristic perturbations in both number density and average size distribution with climatic changes almost disappeared, although they are clearly observed in shallow cores. These results indicate that the air-hydrate crystals grow considerably in deeper parts of the ice sheet, and this growth is accompanied by the diffusion of air molecules in the ice. The permeation coefficient of the air molecules in the ice sheet was estimated from the geometric parameters of the air-hydrate distributions. This is the first practical evidence comparable to the previous model estimations. It allows us to evaluate the impacts of the air-molecule migration in the ice sheet on the paleoclimatic information recorded in the deep ice cores.


Author(s):  
David J. A. Evans

To reconstruct the former extent and dynamics of ice sheets and glaciers requires a knowledge of process-form relationships that goes beyond individual landform types. Instead, glacial geomorphologists need to analyse large areas of glaciated terrain in a more holistic way, combining the whole range of glacial landforms and sediments to reconstruct glacier systems of the past, a subject now known as palaeoglaciology. ‘Glaciers of the past’ explains how the combination of aerial imagery and landform analysis is used in palaeoglaciological reconstruction. Increasingly powerful computers are making it possible to compile sophisticated numerical models that use our knowledge of glaciological processes and ice-core-derived palaeoclimate data to create three-dimensional glacier and ice sheet reconstructions.


2006 ◽  
Vol 2 (1) ◽  
pp. 43-55 ◽  
Author(s):  
W. F. Ruddiman

Abstract. The origin of the major ice-sheet variations during the last 2.7 million years is a long-standing mystery. Neither the dominant 41 000-year cycles in δ18O/ice-volume during the late Pliocene and early Pleistocene nor the late-Pleistocene oscillations near 100 000 years is a linear ("Milankovitch") response to summer insolation forcing. Both responses must result from non-linear behavior within the climate system. Greenhouse gases (primarily CO2) are a plausible source of the required non-linearity, but confusion has persisted over whether the gases force ice volume or are a positive feedback. During the last several hundred thousand years, CO2 and ice volume (marine δ18O) have varied in phase at the 41 000-year obliquity cycle and nearly in phase within the ~100 000-year band. This timing rules out greenhouse-gas forcing of a very slow ice response and instead favors ice control of a fast CO2 response. In the schematic model proposed here, ice sheets responded linearly to insolation forcing at the precession and obliquity cycles prior to 0.9 million years ago, but CO2 feedback amplified the ice response at the 41 000-year period by a factor of approximately two. After 0.9 million years ago, with slow polar cooling, ablation weakened. CO2 feedback continued to amplify ice-sheet growth every 41 000 years, but weaker ablation permitted some ice to survive insolation maxima of low intensity. Step-wise growth of these longer-lived ice sheets continued until peaks in northern summer insolation produced abrupt deglaciations every ~85 000 to ~115 000 years. Most of the deglacial ice melting resulted from the same CO2/temperature feedback that had built the ice sheets. Several processes have the northern geographic origin, as well as the requisite orbital tempo and phasing, to be candidate mechanisms for ice-sheet control of CO2 and their own feedback.


2020 ◽  
Vol 14 (3) ◽  
pp. 833-840 ◽  
Author(s):  
Heiko Goelzer ◽  
Violaine Coulon ◽  
Frank Pattyn ◽  
Bas de Boer ◽  
Roderik van de Wal

Abstract. Estimating the contribution of marine ice sheets to sea-level rise is complicated by ice grounded below sea level that is replaced by ocean water when melted. The common approach is to only consider the ice volume above floatation, defined as the volume of ice to be removed from an ice column to become afloat. With isostatic adjustment of the bedrock and external sea-level forcing that is not a result of mass changes of the ice sheet under consideration, this approach breaks down, because ice volume above floatation can be modified without actual changes in the sea-level contribution. We discuss a consistent and generalised approach for estimating the sea-level contribution from marine ice sheets.


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.


Author(s):  
S. Xiong ◽  
J.-P. Muller

Accumulation of snow and ice over time result in ice sheet layers. These can be remotely sensed where there is a contrast in electromagnetic properties, which reflect variations of the ice density, acidity and fabric orientation. Internal ice layers are assumed to be isochronous, deep beneath the ice surface, and parallel to the direction of ice flow. The distribution of internal layers is related to ice sheet dynamics, such as the basal melt rate, basal elevation variation and changes in ice flow mode, which are important parameters to model the ice sheet. Radar echo sounder is an effective instrument used to study the sedimentology of the Earth and planets. Ice Penetrating Radar (IPR) is specific kind of radar echo sounder, which extends studies of ice sheets from surface to subsurface to deep internal ice sheets depending on the frequency utilised. In this study, we examine a study site where folded ice occurs in the internal ice sheet south of the North Greenland Eemian ice drilling (NEEM) station, where two intersected radar echograms acquired by the Multi-channel Coherent Radar Depth Sounder (MCoRDS) employed in the NASA’s Operation IceBridge (OIB) mission imaged this folded ice. We propose a slice processing flow based on a Radon Transform to trace and extract these two sets of curved ice sheet layers, which can then be viewed in 3-D, demonstrating the 3-D structure of the ice folds.


2018 ◽  
Author(s):  
Jorge Alvarez-Solas ◽  
Rubén Banderas ◽  
Alexander Robinson ◽  
Marisa Montoya

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.


Sign in / Sign up

Export Citation Format

Share Document