ALBERT

Description: We present results from simulations of the Miocene Antarctic ice sheet, that were performed using the 3D thermodynamical ice-sheet model IMAU-ICE (v1.1.1-MIO). Five steady-state present-day simulations were conducted for reference (PI_ref), and 12 experiments using Miocene settings. Each Miocene experiment comprises 11 steady-state and 4 transient simulations. In the README file, the experiments and simulations are listed. IMAU-ICE was run using a 40x40km grid covering the Antarctic continent. Initial conditions were obtained from reconstructions of the Antarctic bathymetry and bedrock topography pertaining to 23 to 24 million years (Myr) ago (dataset doi:10.1594/PANGAEA.923109). The simulations were forced by climate input data obtained from GENESIS simulations with varying CO2 levels (280 to 840 ppm) and Antarctic ice sheet cover (no ice to a large East-Antarctic ice sheet), and with present-day insolation. We utilized a matrix interpolation method to construct the time-varying climate forcing, based on the prescribed CO2 levels and ice cover simulated by IMAU-ICE. For each simulation, we provide the run script, 1D output variables including CO2 level and the sea level contribution of the Antarctic ice sheet, and 3D output variables including ice thickness, bedrock and surface height, surface mass balance, basal mass balance, ice velocities, and ice temperatures. For more information, please contact L.B. Stap at l.b.stap@uu.nl.

Description: In order to investigate the relation between ice sheets and climate in a warmer-than-present world, recent research has focussed on the Late Pliocene, 3.6 to 2.58 million years ago. It is the most recent period in Earth's history when such a warm climate state existed for a significant duration of time. Marine Isotope Stage (MIS) M2 (∼3.3 Myr ago) is a strong positive excursion in benthic oxygen records in the middle of the otherwise warm and relatively stable Late Pliocene. However, the relative contributions to the benthic δ18O signal from deep ocean cooling and growing ice sheets are still uncertain. Here, we present results from simulations of the Late Pliocene with a hybrid ice-sheet–climate model, showing a reconstruction of ice sheet geometry, sea level and atmospheric CO2. Initial experiments simulating the last four glacial cycles indicate that this model yields results which are in good agreement with proxy records in terms of global mean sea level, benthic oxygen isotope abundance, ice-core-derived surface temperature and atmospheric CO2 concentration. For the Late Pliocene, our results show an atmospheric CO2 concentration during MIS M2 of 233–249 ppmv and a drop in global mean sea level of 10 to 25 m. Uncertainties are larger during the warmer periods leading up to and following MIS M2. CO2 concentrations during the warm intervals in the Pliocene, with sea-level high stands of 8–14 m above the present day, varied between 320 and 400 ppmv, lower than indicated by some proxy records but in line with earlier model reconstructions.

Description: Many processes govern the deglaciation of ice sheets. One of the processes that is usually ignored is the calving of ice in lakes that temporarily surround the ice sheet. In order to capture this process a “flood-fill algorithm” is needed. Here we present and evaluate several optimizations to a standard flood-fill algorithm in terms of computational efficiency. As an example, we determine the land–ocean mask for a 1 km resolution digital elevation model (DEM) of North America and Greenland, a geographical area of roughly 7000 by 5000 km (roughly 35 million elements), about half of which is covered by ocean. Determining the land–ocean mask with our improved flood-fill algorithm reduces computation time by 90 % relative to using a standard stack-based flood-fill algorithm. This implies that it is now feasible to include the calving of ice in lakes as a dynamical process inside an ice-sheet model. We demonstrate this by using bedrock elevation, ice thickness and geoid perturbation fields from the output of a coupled ice-sheet–sea-level equation model at 30 000 years before present and determine the extent of Lake Agassiz, using both the standard and improved versions of the flood-fill algorithm. We show that several optimizations to the flood-fill algorithm used for filling a depression up to a water level, which is not defined beforehand, decrease the computation time by up to 99 %. The resulting reduction in computation time allows determination of the extent and volume of depressions in a DEM over large geographical grids or repeatedly over long periods of time, where computation time might otherwise be a limiting factor. The algorithm can be used for all glaciological and hydrological models, which need to trace the evolution over time of lakes or drainage basins in general.

Description: In order to investigate the relation between ice sheets and climate in a warmer-than-present world, recent research has focussed on the Late Pliocene, 3.6 to 2.58 million years ago. It is the most recent period in Earth history when such a climate state existed for a significant duration of time. Marine Isotope Stage (MIS) M2 (~ 3.3 Myr ago) is a strong positive excursion in benthic oxygen records in the middle of the otherwise warm and relatively stable Late Pliocene. However, the relative contributions to the benthic δ18O signal from deep-ocean cooling and growing ice sheets are still uncertain. Here, we present results from simulations of the late Pliocene with a hybrid ice-sheet–climate model, showing a reconstruction of ice sheet geometry, sea-level and atmospheric CO2. Initial experiments simulating the last four glacial cycles indicate that this model yields results which are in good agreement with proxy records in terms of global mean sea level, benthic oxygen isotope abundance, ice core-derived surface temperature and atmospheric CO2 concentration. For the Late Pliocene, our results show an atmospheric CO2 concentration during MIS M2 of 233–249 ppmv, and a drop in global mean sea level of 10 to 25 m. Uncertainties are larger during the warmer periods leading up to and following MIS M2. CO2 concentrations during the warm intervals in the Pliocene, with sea-level high stands of 8–14 m above present-day, varied between 320 and 400 ppmv, lower than indicated by some proxy records but in line with earlier model reconstructions.

Description: Fully coupled ice-sheet-climate modelling over 10,000–100,000-year time scales on high spatial and temporal resolution remains beyond the capability of current computational systems. Hybrid GCM-ice-sheet modelling offers a middle ground, balancing the need to accurately capture both long-term processes, in particular circulation driven changes in precipitation, and processes requiring a high spatial resolution like ablation. Here, we present and evaluate a model set-up that forces the ANICE 3D thermodynamic ice-sheet-shelf model calculating all ice on Earth, with pre-calculated output from several steady-state simulations with the HadCM3 general circulation model (GCM), using a so-called matrix method of coupling both components, where simulations with various levels of pCO2 and ice-sheet configuration are combined to form a time-continuous transient climate forcing consistent with the modelled ice-sheets. We address the difficulties in downscaling low-resolution GCM output to the higher-resolution grid of an ice-sheet model, and account for differences between GCM and ice-sheet model surface topography ranging from interglacial to glacial conditions. As a benchmark experiment to assess the validity of this model set-up, we perform a simulation of the entire last glacial cycle, from 120kyr ago to present-day. The simulated eustatic sea-level drop at the Last Glacial maximum (LGM) for the combined Antarctic, Greenland, Eurasian and North-American ice-sheets amounts to 100m, in line with many other studies. The simulated ice-sheets at LGM agree well with the ICE-5G reconstruction and the more recent DATED-1 reconstruction in terms of total volume and geographical location of the ice sheets. Moreover, modelled benthic oxygen isotope abundance and the relative contributions from global ice volume and deep-water temperature agree well with available data, as do surface temperature histories for the Greenland and Antarctic ice-sheets. This model strategy can be used to create time-continuous ice-sheet distribution and sea-level reconstructions for geological periods up to several millions of years in duration, capturing climate model driven variations in the mass balance of the ice sheet.

Description: Fully coupled ice-sheet–climate modelling over 10 000–100 000-year timescales at high spatial and temporal resolution remains beyond the capability of current computational systems. Forcing an ice-sheet model with precalculated output from a general circulation model (GCM) offers a middle ground, balancing the need to accurately capture both long-term processes, in particular circulation-driven changes in precipitation, and processes requiring a high spatial resolution like ablation. Here, we present and evaluate a model set-up that forces the ANICE 3-D thermodynamic ice-sheet–shelf model calculating the four large continental ice sheets (Antarctica, Greenland, North America, and Eurasia) with precalculated output from two steady-state simulations with the HadCM3 (GCM) using a so-called matrix method of coupling both components, whereby simulations with various levels of pCO2 and ice-sheet configuration are combined to form a time-continuous transient climate forcing consistent with the modelled ice sheets. We address the difficulties in downscaling low-resolution GCM output to the higher-resolution grid of an ice-sheet model and account for differences between GCM and ice-sheet model surface topography ranging from interglacial to glacial conditions. Although the approach presented here can be applied to a matrix with any number of GCM snapshots, we limited our experiments to a matrix of only two snapshots. As a benchmark experiment to assess the validity of this model set-up, we perform a simulation of the entire last glacial cycle from 120 kyr ago to present day. The simulated eustatic sea-level drop at the Last Glacial Maximum (LGM) for the combined Antarctic, Greenland, Eurasian, and North American ice sheets amounts to 100 m, in line with many other studies. The simulated ice sheets at the LGM agree well with the ICE-5G reconstruction and the more recent DATED-1 reconstruction in terms of total volume and geographical location of the ice sheets. Moreover, modelled benthic oxygen isotope abundance and the relative contributions from global ice volume and deep-water temperature agree well with available data, as do surface temperature histories for the Greenland and Antarctic ice sheets. This model strategy can be used to create time-continuous ice-sheet distribution and sea-level reconstructions for geological periods up to several million years in duration, capturing climate-model-driven variations in the mass balance of the ice sheet.

Description: Observations and models agree that the Greenland Ice Sheet (GrIS) surface mass balance (SMB) has decreased since the end of the 1990s due to an increase in meltwater runoff and that this trend will accelerate in the future. However, large uncertainties remain, partly due to different approaches for modelling the GrIS SMB, which have to weigh physical complexity or low computing time, different spatial and temporal resolutions, different forcing fields, and different ice sheet topographies and extents, which collectively make an inter-comparison difficult. Our GrIS SMB model intercomparison project (GrSMBMIP) aims to refine these uncertainties by intercomparing 13 models of four types which were forced with the same ERA-Interim reanalysis forcing fields, except for two global models. We interpolate all modelled SMB fields onto a common ice sheet mask at 1 km horizontal resolution for the period 1980–2012 and score the outputs against (1) SMB estimates from a combination of gravimetric remote sensing data from GRACE and measured ice discharge; (2) ice cores, snow pits and in situ SMB observations; and (3) remotely sensed bare ice extent from MODerate-resolution Imaging Spectroradiometer (MODIS). Spatially, the largest spread among models can be found around the margins of the ice sheet, highlighting model deficiencies in an accurate representation of the GrIS ablation zone extent and processes related to surface melt and runoff. Overall, polar regional climate models (RCMs) perform the best compared to observations, in particular for simulating precipitation patterns. However, other simpler and faster models have biases of the same order as RCMs compared with observations and therefore remain useful tools for long-term simulations or coupling with ice sheet models. Finally, it is interesting to note that the ensemble mean of the 13 models produces the best estimate of the present-day SMB relative to observations, suggesting that biases are not systematic among models and that this ensemble estimate can be used as a reference for current climate when carrying out future model developments. However, a higher density of in situ SMB observations is required, especially in the south-east accumulation zone, where the model spread can reach 2 m w.e. yr−1 due to large discrepancies in modelled snowfall accumulation.

Description: Understanding the evolution of, and the interactions between, ice sheets and the global climate over geological timescales is important for being able to project their future evolution. However, direct observational evidence of past CO2 concentrations, and the implied radiative forcing, only exists for the past 800 000 years. Records of benthic δ18O date back millions of years but contain signals from both land ice volume and ocean temperature. In recent years, inverse forward modelling has been developed as a method to disentangle these two signals, resulting in mutually consistent reconstructions of ice volume, temperature, and CO2. We use this approach to force a hybrid ice-sheet–climate model with a benthic δ18O stack, reconstructing the evolution of the ice sheets, global mean sea level, and atmospheric CO2 during the late Pliocene and the Pleistocene, from 3.6 million years (Myr) ago to the present day. During the warmer-than-present climates of the late Pliocene, reconstructed CO2 varies widely, from 320–440 ppmv for warm periods to 235–250 ppmv for the early glacial excursion ∼3.3 million years ago. Sea level is relatively stable during this period, with maxima of 6–14 m and minima of 12–26 m during glacial episodes. Both CO2 and sea level are within the wide ranges of values covered by available proxy data for this period. Our results for the Pleistocene agree well with the ice-core CO2 record, as well as with different available sea-level proxy data. For the Early Pleistocene, 2.6–1.2 Myr ago, we simulate 40 kyr glacial cycles, with interglacial CO2 decreasing from 280–300 ppmv at the beginning of the Pleistocene to 250–280 ppmv just before the Mid-Pleistocene Transition (MPT). Peak glacial CO2 decreases from 220–250 to 205–225 ppmv during this period. After the MPT, when the glacial cycles change from 40 to 80 120 kyr cyclicity, the glacial–interglacial contrast increases, with interglacial CO2 varying between 250–320 ppmv and peak glacial values decreasing to 170–210 ppmv.

Description: The Greenland Ice Sheet (GrIS) mass loss has been accelerating at a rate of about 20 ± 10 Gt/yr2 since the end of the 1990's, with around 60 % of this mass loss directly attributed to enhanced surface meltwater runoff. However, in the climate and glaciology communities, different approaches exist on how to model the different surface mass balance (SMB) components using: (1) complex physically-based climate models which are computationally expensive; (2) intermediate complexity energy balance models; (3) simple and fast positive degree day models which base their inferences on statistical principles and are computationally highly efficient. Additionally, many of these models compute the SMB components based on different spatial and temporal resolutions, with different forcing fields as well as different ice sheet topographies and extents, making inter-comparison difficult. In the GrIS SMB model intercomparison project (GrSMBMIP) we address these issues by forcing each model with the same data (i.e., the ERA-Interim reanalysis) except for two global models for which this forcing is limited to the oceanic conditions, and at the same time by interpolating all modelled results onto a common ice sheet mask at 1 km horizontal resolution for the common period 1980–2012. The SMB outputs from 13 models are then compared over the GrIS to (1) SMB estimates using a combination of gravimetric remote sensing data from GRACE and measured ice discharge, (2) ice cores, snow pits, in-situ SMB observations, and (3) remotely sensed bare ice extent from MODerate-resolution Imaging Spectroradiometer (MODIS). Our results reveal that the mean GrIS SMB of all 13 models has been positive between 1980 and 2012 with an average of 340 ± Gt/yr, but has decreased at an average rate of −7.3 Gt/yr2 (with a significance of 96 %), mainly driven by an increase of 8.0 Gt/yr2 (with a significance of 98 %) in meltwater runoff. Spatially, the largest spread among models can be found around the margins of the ice sheet, highlighting the need for accurate representation of the GrIS ablation zone extent and processes driving the surface melt. In addition, a higher density of in-situ SMB observations is required, especially in the south-east accumulation zone, where the model spread can reach 2 mWE/yr due to large discrepancies in modelled snowfall accumulation. Overall, polar regional climate models (RCMs) perform the best compared to observations, in particular for simulating precipitation patterns. However, other simpler and faster models have biases of same order than RCMs with observations and remain then useful tools for long-term simulations. Finally, it is interesting to note that the ensemble mean of the 13 models produces the best estimate of the present day SMB relative to observations, suggesting that biases are not systematic among models.