Search

During the Mid-Pleistocene transition (MPT) the dominant periodicity of glacial cycles increased from 41 thousand years (kyr) to an average of 100 kyr, without any appreciable change in the orbital pacing. As the MPT is not a linear response to orbital forcing, it must have resulted from feedback processes in the Earth system. However, the precise mechanisms underlying the transition are still under debate. In this study, we investigate the MPT by simulating the Northern Hemisphere ice sheet evolution over the past 1.5 million years. The transient climate forcing of the ice-sheet model was obtained using a matrix method, by interpolating between two snapshots of global climate model simulations. Changes in climate forcing are caused by variations in CO2, insolation, as well as implicit climate–ice sheet feedbacks. Using this method, we were able to capture glacial-interglacial variability during the past 1.5 million years and reproduce the shift from 41 kyr to 100 kyr cycles without any additional drivers. Instead, the modelled frequency change results from the prescribed CO2 combined with orbital forcing, and ice sheet feedbacks. Early Pleistocene terminations are initiated by insolation maxima. After the MPT, low CO2 levels can compensate insolation maxima which favour deglaciation, leading to an increasing glacial cycle periodicity. These deglaciations are also prevented by a relatively small North American ice sheet, which, through its location and feedback processes, can generate a relatively stable climate. Larger North American ice sheets become more sensitive to small temperature increases. Therefore, Late Pleistocene terminations are facilitated by the large ice-sheet volume, were small changes in temperature lead to self-sustained melt instead. This concept is confirmed by experiments using constant insolation or CO2. The constant CO2 experiments generally capture only the Early Pleistocene cycles, while those with constant insolation only capture Late Pleistocene cycles. Additionally, we find that a lowering of CO2concentrations leads to an increasing number of insolation maxima that fail to initiate terminations. These results therefore suggest a regime shift, where during the Early Pleistocene, glacial cycles are dominated by orbital oscillations, while Late Pleistocene cycles tend to be more dominated by CO2. This implies that the MPT can be explained by a decrease in glacial CO2 concentration superimposed on orbital forcing. [ABSTRACT FROM AUTHOR]

During the glacial cycles of the past 800 000 years, Eurasia and North America were periodically covered by large ice sheets, causing up to 100 m of sea-level change. While Late Pleistocene glacial cycles typically lasted 80 000–120 000 years, the termination phases were completed in only 10 000 years. During these glacial terminations, the North American and Eurasian ice sheets retreated, which created large proglacial lakes in front of the ice-sheet margin. Proglacial lakes accelerate deglaciation as they facilitate the formation of ice shelves at the southern margins of the North American and Eurasian ice sheets. These ice shelves are characterized by basal melting, low surface elevations, and negligible friction at the base. Here, we use an ice-sheet model to quantify the (combined) effects of proglacial lakes on Late Pleistocene glacial terminations by examining their interplay with glacial isostatic adjustment (GIA) and basal sliding. We find that proglacial lakes accelerate the deglaciation of ice sheets mainly because there is an absence of basal friction underneath ice shelves. If friction underneath grounded ice is applied to floating ice, full deglaciation is postponed by a few millennia, resulting in more ice remaining during interglacial periods and no extensive ice shelves forming. Additionally, the large uncertainty in melt rates underneath lacustrine ice shelves translates to an uncertainty in the timing of the termination of up to a millennium. Proglacial lakes are created by depressions in the landscape that remain after an ice sheet has retreated. The depth, size, and timing of proglacial lakes depend on the rate of bedrock rebound. We find that if bedrock rebounds within a few centuries (rather than a few millennia), the mass loss rate of the ice sheet is substantially reduced. This is because fast bedrock rebound prevents the formation of extensive proglacial lakes. Additionally, a decrease in ice thickness is partly compensated for by faster bedrock rebound, resulting in a higher surface elevation; lower temperatures; and a higher surface mass balance, which delays deglaciation. We find that a very long bedrock relaxation time does not substantially affect terminations, but it may lead to a delayed onset of the next glacial period. This is because inception regions, such as northwestern Canada, remain below sea level throughout the preceding interglacial period. [ABSTRACT FROM AUTHOR]

Deltas and coastal plains worldwide developed under the influence of relative sea level rise (RSLR) during the Holocene. In the Netherlands, Holocene RSLR results from both regional sea-level rise and regional subsidence patterns, mainly caused by glacial isostatic adjustment (GIA: Scandinavian forebulge collapse) and longer-term North Sea Basin tectono-sedimentary subsidence. Past coastal and inland water levels are preserved in geological indicators marking the gradual drowning of an area, for example basal peats. Such geological water-level indicators have been used in the Netherlands for varying types of research. However, uniform overviews of these data exist only for smaller local subsets and not for the entire Netherlands. In this paper we present a data set of 712 Holocene water-level indicators from the Dutch coastal plain that are relevant for studying RSLR and regional subsidence, compiled in HOLSEA workbook format. This format was expanded to allow for registering basal-peat type geological indicators, documenting Dutch-setting specific parameters and accompanying uncertainties, to assess indicative meaning, and to appropriately correct the raw vertical positions of the indicators. Overall, our new, internally consistent, expanded documentation provided for the water-level indicators encourages users to choose the information relevant for their research and report RSLR uncertainties transparently. From the indicators, 59 % was collected in 1950–2000, mainly in academic studies and survey mapping campaigns; 37 % was collected in 2000–2020 in academic studies and archaeological surveying projects, 4 % was newly collected (this study), the latter mainly in previously under sampled central and northern Netherlands regions. Prominent regional differences exist in the vertical position and abundance of the indicators. Older indicators in our data set are mostly located in the deeper seaward area of the Netherlands. These indicators correspond well with previous transgression reconstructions, that are partly based on the same data. The younger, landwards set of indicators in the Rhine-Meuse central and Flevoland regions corresponds with the transgression phase reaching further inland, from 8000 cal. BP onwards. Northern indicators of Middle Holocene age (8–5 ka cal. BP), in general lie 2–3 meters lower compared to those in the south. For younger data this difference is less, showing spatial and temporal variation in RSLR throughout the Netherlands. [ABSTRACT FROM AUTHOR]

When different flooding drivers co‐occur, they can cause compound floods. Despite the potential impact of compound flooding, few studies have projected how the joint probability of flooding drivers may change. Furthermore, existing projections may not be very robust, as they are based on only 5 to 6 climate model simulations. Here, we use a large ensemble of simulations from the Coupled Model Intercomparison Project 6 (CMIP6) to project changes in the joint probability of extreme storm surges and precipitation at European tide gauges under a medium and high emissions scenario, enabled by data‐proximate cloud computing and statistical storm surge modeling. We find that the joint probability will increase in the northwest and decrease in most of the southwest of Europe. Averaged over Europe, the absolute magnitude of these changes is 36%–49% by 2080, depending on the scenario. The large‐scale changes in the joint probability of extreme storm surges and precipitation are similar to those in the joint probability of extreme wind speeds and precipitation, but locally, differences can exceed the changes themselves. Due to internal climate variability and inter‐model differences, projections based on simulations of only 5 to 6 randomly chosen CMIP6 models have a probability of higher than 10% to differ qualitatively from projections based on all CMIP6 simulations in multiple regions, especially under the medium emissions scenario and earlier in the twenty‐first century. Therefore, our results provide a more robust and less uncertain representation of changes in the potential for compound flooding in Europe than previous projections. Plain Language Summary: Extreme storm surges, rainfall or river discharge can cause flooding. When these events happen at the same time, even more severe flooding may follow. Climate change could affect the odds that drivers of flooding coincide, potentially leading to larger flood risk. However, few scientists have tried to compute such changes, using only a few different computer models of our climate. Here, we use a much larger set of climate models to compute how the odds that an extreme storm surge coincides with extreme precipitation could change in the future. We find that at the coasts of northwestern Europe, those odds will increase, whereas in southwestern Europe, they will mostly decrease. On average, the changes will be as large as 36%–49% of the current odds, depending on whether the concentration of greenhouse gases in the atmosphere will increase by a medium or a large amount. When we use smaller sets of climate models for our calculations, we get substantially different results in some cases. In conclusion, by using a larger set of climate models than previous studies, we have made more robust computations of how the odds that extreme storm surges and precipitation coincide will change in Europe. Key Points: We project changes in the joint probability of storm surge and precipitation extremes based on a large ensemble of model simulations from the Coupled Model Intercomparison Project 6The joint probability will increase in the northwest and decrease in the southwest of Europe, with an average absolute magnitude of 36%–49%Especially under lower emissions, often more than 5 or 6 climate model simulations are needed to draw robust conclusions on these changes [ABSTRACT FROM AUTHOR]

The strongly varying benthic δ18 O levels of the early and mid-Miocene (23 to 14 Myr ago) are primarily caused by a combination of changes in Antarctic Ice Sheet (AIS) volume and deep-ocean temperatures. These factors are coupled since AIS changes affect deep-ocean temperatures. It has recently been argued that this is due to changes in ice sheet area rather than volume because area changes affect the surface albedo. This finding would be important when the transient AIS grows relatively faster in extent than in thickness, which we test here. We analyse simulations of Miocene AIS variability carried out using the three-dimensional ice sheet model IMAU-ICE forced by warm (high CO 2 , no ice) and cold (low CO 2 , large East AIS) climate snapshots. These simulations comprise equilibrium and idealized quasi-orbital transient runs with strongly varying CO 2 levels (280 to 840 ppm). Our simulations show a limited direct effect of East AIS changes on Miocene orbital-timescale benthic δ18 O variability because of the slow build-up of volume. However, we find that relative to the equilibrium ice sheet size, the AIS area adapts significantly faster and more strongly than volume to the applied forcing variability. Consequently, during certain intervals the ice sheet is receding at the margins, while ice is still building up in the interior. That means the AIS does not adapt to a changing equilibrium size at the same rate or with the same sign everywhere. Our results indicate that the Miocene Antarctic Ice Sheet affects deep-ocean temperatures more than its volume suggests. [ABSTRACT FROM AUTHOR]

Future projections of sea-level rise under strong warming scenarios are dominated by mass loss in the marine-grounded sectors of West Antarctica, where thinning shelves as a result of warming oceans can lead to reduced buttressing. This consequently leads to accelerated flow from the upstream grounded ice. However, the relation between warming oceans and increased melt rates under the shelves is very uncertain, especially when interactions with the changing shelf geometry are considered. Here, we compare six widely used, highly parameterised formulations relating sub-shelf melt to thermal forcing. We implemented them in an ice-sheet model, and applied the resulting set-up to an idealised-geometry setting, as well as to the Antarctic ice sheet. In our simulations, the differences in modelled ice-sheet evolution resulting from the choice of parameterisation, as well as the choice of numerical scheme used to apply sub-shelf melt near the grounding line, generally are larger than differences from ice-dynamical processes such as basal sliding, as well as uncertainties from the forcing scenario of the model providing the ocean forcing. This holds for the idealised-geometry experiments as well as for the experiments using a realistic Antarctic topography. [ABSTRACT FROM AUTHOR]

Glacial isostatic adjustment (GIA) has a stabilizing effect on the evolution of the Antarctic ice sheet by reducing the grounding line migration following ice melt. The timescale and strength of this feedback depends on the spatially varying viscosity of the Earth's mantle. Most studies assume a relatively long and laterally homogenous response time of the bedrock. However, the mantle viscosity is spatially variable, with a high mantle viscosity beneath East Antarctica and a low mantle viscosity beneath West Antarctica. For this study, we have developed a new method to couple a 3D GIA model and an ice sheet model to study the interaction between the solid Earth and the Antarctic ice sheet during the last glacial cycle. With this method, the ice sheet model and GIA model exchange ice thickness and bedrock elevation during a fully coupled transient experiment. The feedback effect is taken into account with a high temporal resolution, where the coupling time steps between the ice sheet and GIA model are 5000 years over the glaciation phase and vary between 500 and 1000 years over the deglaciation phase of the last glacial cycle. During each coupling time step, the bedrock elevation is adjusted at every ice sheet model time step, and the deformation is computed for a linearly changing ice load. We applied the method using the ice sheet model ANICE and a 3D GIA finite element model. We used results from a regional seismic model for Antarctica embedded in the global seismic model SMEAN2 to determine the patterns in the mantle viscosity. The results of simulations over the last glacial cycle show that differences in mantle viscosity of an order of magnitude can lead to differences in the grounding line position up to 700 km and to differences in ice thickness of the order of 2 km for the present day near the Ross Embayment. These results underline and quantify the importance of including local GIA feedback effects in ice sheet models when simulating the Antarctic ice sheet evolution over the last glacial cycle. [ABSTRACT FROM AUTHOR]

A major source of uncertainty in future sea level projections is the ocean-driven basal melt of Antarctic ice shelves. While ice sheet models require a kilometre-scale resolution to realistically resolve ice shelf stability and grounding line migration, global or regional 3D ocean models are computationally too expensive to produce basal melt forcing fields at this resolution on long timescales. To bridge this resolution gap, we introduce the 2D numerical model LADDIE (one-layer Antarctic model for dynamical downscaling of ice–ocean exchanges), which allows for the computationally efficient modelling of detailed basal melt fields. The model is open source and can be applied easily to different geometries or different ocean forcings. The aim of this study is threefold: to introduce the model to the community, to demonstrate its application and performance in two use cases, and to describe and interpret new basal melt patterns simulated by this model. The two use cases are the small Crosson–Dotson Ice Shelf in the warm Amundsen Sea region and the large Filchner–Ronne Ice Shelf in the cold Weddell Sea. At ice-shelf-wide scales, LADDIE reproduces observed patterns of basal melting and freezing in warm and cold environments without the need to re-tune parameters for individual ice shelves. At scales of 0.5–5 km, which are typically unresolved by 3D ocean models and poorly constrained by observations, LADDIE produces plausible basal melt patterns. Most significantly, the simulated basal melt patterns are physically consistent with the applied ice shelf topography. These patterns are governed by the topographic steering and Coriolis deflection of meltwater flows, two processes that are poorly represented in basal melt parameterisations. The kilometre-scale melt patterns simulated by LADDIE include enhanced melt rates in grounding zones and basal channels and enhanced melt or freezing in shear margins. As these regions are critical for ice shelf stability, we conclude that LADDIE can provide detailed basal melt patterns at the essential resolution that ice sheet models require. The physical consistency between the applied geometry and the simulated basal melt fields indicates that LADDIE can play a valuable role in the development of coupled ice–ocean modelling. [ABSTRACT FROM AUTHOR]

During the glacial cycles of the past 800 thousand years, Eurasia and North America were periodically covered by large ice sheets. While the Late Pleistocene glacial cycles typically lasted 80 - 120 thousand years, the termination phases only took 10 thousand years to complete. During these glacial terminations, the North American and Eurasian ice sheets retreated which created large proglacial lakes in front of the ice sheet margin. Proglacial lakes accelerate the deglaciation as they can facilitate ice shelves in the southern margins of the North American and the Eurasian ice sheets. Ice shelves are characterized by basal melting, low surface elevations and negligible friction at the base. Here we quantify the effect of proglacial lakes, and the combined effect with glacial isostatic adjustment (GIA) on Late Pleistocene glacial terminations. We find that proglacial lakes accelerate the deglaciation of the ice sheets mainly because of the absence of basal friction underneath ice shelves. If the friction underneath grounded ice is applied to floating ice, we find that full deglaciation is postponed by a few millennia, the Barents-Kara Sea region does not fully deglaciate, and there are no extensive ice shelves. Additionally, the large uncertainty in melt rates underneath lacustrine ice shelves translates to an uncertainty in the timing of the termination of only a few centuries at most. Proglacial lakes are created by the depression in the landscape that linger after the ice sheet has retreated. The depth, size and timing of proglacial lakes depend on the bedrock rebound. We find that if the bedrock rebounds within a few centuries, instead of a few millennia, the mass loss rate of the ice sheet is substantially reduced. This is because fast bedrock rebound prevents the formation of extensive proglacial lakes. Additionally, a decrease in thickness is partly compensated by the faster bedrock rebound, resulting in a higher surface elevation with lower temperatures and higher surface mass balance delaying deglaciation. We find that a very long bedrock relaxation time does not affect terminations substantially, but will lead to a later inception of the next glacial period. This is because initial inception regions, such as North-Western Canada, remain below sea level throughout the preceding interglacial period. [ABSTRACT FROM AUTHOR]

The strongly varying benthic δ18O levels of the early and mid-Miocene (23 to 14 Myr ago) are primarily caused by a combination of changes in Antarctic ice sheet (AIS) volume and deep ocean temperatures. These factors are coupled since AIS changes affect deep ocean temperatures. It has recently been argued that this is due to changes in ice sheet area rather than volume, because area changes affect the surface albedo. This would be important when the transient AIS grows relatively faster in extent than in thickness, which we test here. We analyse simulations of Miocene AIS variability carried out using the three-dimensional ice-sheet model IMAU-ICE forced by warm (high CO2, no ice) and cold (low CO2, large East-AIS) climate snapshots. These simulations comprise equilibrium and idealised quasi-orbital transient runs with strongly varying CO2 levels (280 to 840 ppm). Our simulations show limited direct effect of East-AIS changes on Miocene orbital timescale benthic δ18O variability, because of the slow build-up of volume. However, we find that AIS area responds significantly faster and more strongly than volume to the applied forcing variability. Consequently, during certain intervals the ice sheet is receding at the margins, while ice is still building up in the interior. That means the AIS does not adapt to a changing equilibrium size at the same rate or with the same sign everywhere. Our results indicate that the Miocene Antarctic ice sheet affects deep ocean temperatures more than its volume suggests. [ABSTRACT FROM AUTHOR]

Subglacial bed roughness is one of the main factors controlling the rate of future Antarctic ice-sheet retreat and also one of the most uncertain. A common technique to constrain the bed roughness using ice-sheet models is basal inversion, tuning the roughness to reproduce the observed present-day ice-sheet geometry and/or surface velocity. However, many other factors affecting ice-sheet evolution, such as the englacial temperature and viscosity, the surface and basal mass balance, and the subglacial topography, also contain substantial uncertainties. Using a basal inversion technique intrinsically causes any errors in these other quantities to lead to compensating errors in the inverted bed roughness. Using a set of idealised-geometry experiments, we quantify these compensating errors and investigate their effect on the dynamic response of the ice sheet to a prescribed forcing. We find that relatively small errors in ice viscosity and subglacial topography require substantial compensating errors in the bed roughness in order to produce the same steady-state ice sheet, obscuring the realistic spatial variability in the bed roughness. When subjected to a retreat-inducing forcing, we find that these different parameter combinations, which per definition of the inversion procedure result in the same steady-state geometry, lead to a rate of ice volume loss that can differ by as much as a factor of 2. This implies that ice-sheet models that use basal inversion to initialise their model state can still display a substantial model bias despite having an initial state which is close to the observations. [ABSTRACT FROM AUTHOR]

During the last glacial cycle (LGC), ice sheets covered large parts of Eurasia and North America, which resulted in ∼120 m of sea level change. Ice sheet–climate interactions have considerable influence on temperature and precipitation patterns and therefore need to be included when simulating this time period. Ideally, ice sheet–climate interactions are simulated by a high-resolution Earth system model. While these models are capable of simulating climates at a certain point in time, such as the pre-industrial (PI) or the Last Glacial Maximum (LGM; 21 000 years ago), a full transient glacial cycle is currently computationally unfeasible as it requires a too-large amount of computation time. Nevertheless, ice sheet models require forcing that captures the gradual change in climate over time to calculate the accumulation and melt of ice and its effect on ice sheet extent and volume changes. Here we simulate the LGC using an ice sheet model forced by LGM and PI climates. The gradual change in climate is modelled by transiently interpolating between pre-calculated results from a climate model for the LGM and the PI. To assess the influence of ice sheet–climate interactions, we use two different interpolation methods: the climate matrix method, which includes a temperature–albedo and precipitation–topography feedback, and the glacial index method, which does not. To investigate the sensitivity of the results to the prescribed climate forcing, we use the output of several models that are part of the Paleoclimate Modelling Intercomparison Project Phase III (PMIP3). In these simulations, ice volume is prescribed, and the climate is reconstructed with a general circulation model (GCM). Here we test those models by using their climate to drive an ice sheet model over the LGC. We find that the ice volume differences caused by the climate forcing exceed the differences caused by the interpolation method. Some GCMs produced unrealistic LGM volumes, and only four resulted in reasonable ice sheets, with LGM Northern Hemisphere sea level contribution ranging between 74–113 m with respect to the present day. The glacial index and climate matrix methods result in similar ice volumes at the LGM but yield a different ice evolution with different ice domes during the inception phase of the glacial cycle and different sea level rates during the deglaciation phase. The temperature–albedo feedback is the main cause of differences between the glacial index and climate matrix methods. [ABSTRACT FROM AUTHOR]

Sea-level science has seen many recent developments in observations and modelling of the different contributions and the total mean sea-level change. In this overview, we discuss (1) the evolution of the Intergovernmental Panel on Climate Change (IPCC) projections, (2) how the projections compare to observations and (3) the outlook for further improving projections. We start by discussing how the model projections of 21st century sea-level change have changed from the IPCC AR5 report (2013) to SROCC (2019) and AR6 (2021), highlighting similarities and differences in the methodologies and comparing the global mean and regional projections. This shows that there is good agreement in the median values, but also highlights some differences. In addition, we discuss how the different reports included high-end projections. We then show how the AR5 projections (from 2007 onwards) compare against the observations and find that they are highly consistent with each other. Finally, we discuss how to further improve sea-level projections using high-resolution ocean modelling and recent vertical land motion estimates. [ABSTRACT FROM AUTHOR]

Ice-dynamical processes constitute a large uncertainty in future projections of sea-level rise caused by anthropogenic climate change. Improving our understanding of these processes requires ice-sheet models that perform well at simulating both past and future ice-sheet evolution. Here, we present version 2.0 of the ice-sheet model IMAU-ICE, which uses the depth-integrated viscosity approximation (DIVA) to solve the stress balance. We evaluate its performance in a range of benchmark experiments, including simple analytical solutions and both schematic and realistic model intercomparison exercises. IMAU-ICE has adopted recent developments in the numerical treatment of englacial stress and sub-shelf melt near the grounding line, which result in good performance in experiments concerning grounding-line migration (MISMIP, MISMIP +) and buttressing (ABUMIP). This makes it a model that is robust, versatile, and user-friendly, which will provide a firm basis for (palaeo-)glaciological research in the coming years. [ABSTRACT FROM AUTHOR]

A major source of uncertainty in future sea-level projections is the ocean-driven basal melt of Antarctic ice shelves. Whereas ice sheet models require a kilometer-scale resolution to realistically resolve ice shelf stability and grounding line migration, global or regional 3D ocean models are computationally too expensive to produce basal melt forcing fields at this resolution. To bridge this resolution gap, we introduce the 2D numerical model LADDIE (one-Layer Antarctic model for Dynamical Downscaling of Ice-ocean Exchanges) which allows for the computationally efficient modeling of basal melt rates. The model is flexible, and can be forced with output from coarse 3D ocean models or with vertical profiles of offshore temperature and salinity. In this study, we describe the model equations and numerics. To illustrate and validate the model performance, we apply the model to two test cases: the small Crosson-Dotson Ice Shelf in the warm Amundsen Sea region, and the large Filchner-Ronne Ice Shelf in the cold Weddell Sea. At ice-shelf wide scales, LADDIE reproduces observed patterns of basal melt and freezing that are also well reproduced by 3D ocean models. At scales of 0.5-5 km, which are unresolved by 3D ocean models and poorly constrained by observations, LADDIE produces plausible basal melt patterns. Most significantly, the simulated basal melt patterns are physically consistent with the applied ice shelf topography. These patterns are governed by the topographic steering and Coriolis deflection of meltwater flows, two processes that are poorly represented in basal melt parameterisations. The kilometer-scale melt patterns simulated by LADDIE include enhanced melt rates in basal channels, in some shear margins, and nearby grounding lines. As these regions are critical for ice shelf stability, we conclude that LADDIE can provide detailed basal melt patterns at the essential resolution that ice sheet models require. The physical consistency between the applied geometry and the simulated basal melt fields indicates that LADDIE can play a valuable role in the development of coupled ice-ocean modeling. [ABSTRACT FROM AUTHOR]

Glacial Isostatic Adjustment (GIA) has a stabilizing effect on the evolution of the Antarctic Ice Sheet by reducing the grounding line migration that follows ice melt. The timescale and strength of this feedback depend on the spatially varying viscosity of the Earth's mantle. Most studies assume a relatively high laterally homogenous response time of the bedrock. However, viscosity is spatially variable with a high viscosity beneath East Antarctica, and a low viscosity beneath West Antarctica. For this study, we have developed a new method to couple a 3D GIA model and an ice-sheet model to study the interaction between the Solid Earth and the Antarctic Ice Sheet during the last glacial cycle. The feedback effect into account on a high temporal resolution by using coupling time steps of 500 years. We applied the method using the ice-sheet model ANICE, a 3D GIA FE model, and results from a seismic model to determine the patterns in the viscosity. The results of simulations over the Last Glacial Cycle show that differences in viscosity of an order of magnitude can lead to differences in grounding line position up to 500 km, to differences in ice thickness in the order of 1.5 km. These results underline and quantify the importance of including local GIA feedback effects in ice-sheet models when simulating the Antarctic Ice Sheet evolution over the Last Glacial Cycle. [ABSTRACT FROM AUTHOR]

During the Last Glacial Cycle (LGC), ice sheets covered large parts of Eurasia and North America which resulted in ~120 meters of sea level change. Ice sheet - climate interactions have considerable influence on temperature and precipitation patterns, and therefore need to be included when simulating this time period. Ideally, ice sheet - climate interactions are simulated by a high-resolution earth system model. While these models are capable of simulating climates at a certain point in time, such as the Pre-Industrial (PI) or the Last Glacial Maximum (LGM; 21,000 years ago), a full glacial cycle is currently unfeasible as it requires a too large amount of computation time. Nevertheless, ice-sheet models require forcing that captures the gradual change in climate over time to calculate the accumulation and melt of ice and its effect on ice sheet extent and volume changes. Here we simulate the LGC using an ice sheet model forced by LGM and PI climates. The gradual change in climate is modelled by transiently interpolating between pre-calculated results from a climate model for the LGM and the PI. To assess the influence of ice sheet - climate interactions, we use two different interpolation methods: The climate matrix method, which includes these interactions, and the glacial index method, which does not. To investigate the sensitivity of the results to the prescribed climate forcing, we use the output of several models that are part of the Paleoclimate Modelling Intercomparison Project Phase III (PMIP3). In these simulations, ice volume is prescribed and the climate is reconstructed. Here we test those models by using their climate to drive an ice sheet model over the LGC. We find that the differences caused by the climate forcing exceeds the differences caused by the interpolation method. Some General Circulation Models (GCMs) produced unrealistic LGM volumes and only four resulted in reasonable ice sheets with LGM Northern Hemisphere sea level contribution ranging between 74 - 113 meters with respect to the present day. The glacial index and climate matrix methods result in similar ice volumes at LGM but yield a different ice evolution with different ice domes during the inception phase of the glacial cycle, and different sea-level rates during the deglaciation phase. The temperature-albedo feedback is the main cause of differences between the glacial index and climate matrix methods. [ABSTRACT FROM AUTHOR]

Subglacial bed roughness is one of the main factors controlling the rate of future Antarctic ice-sheet retreat, and also one of the most uncertain. A common technique to constrain the bed roughness using ice-sheet models is basal inversion, tuning the roughness to reproduce the observed present-day ice-sheet geometry and/or surface velocity. However, many other factors affecting ice-sheet evolution, such as the englacial temperature and viscosity, the surface and basal mass balance, and the subglacial topography, also contain substantial uncertainties. Using a basal inversion technique intrinsically causes any errors in these other quantities, to lead to compensating errors in the inverted bed roughness. Using a set of idealised-geometry experiments, we quantify these compensating errors and investigate their effect on the dynamic response of the ice-sheet to a prescribed forcing. We find that relatively small errors in ice viscosity and subglacial topography require substantial compensating errors in the bed roughness in order to produce the same steady-state ice sheet, obscuring the realistic spatial variability in the bed roughness. When subjected to a retreat-inducing forcing, we find that these different parameter combinations, which per definition of the inversion procedure result in the same steady-state geometry, lead to a rate of ice volume loss that can differ by as much as a factor of two. This implies that ice-sheet models that use basal inversion to initialise their model state can still display a substantial model bias despite having an initial state which is close to the observations. [ABSTRACT FROM AUTHOR]

Benthic δ18 O levels vary strongly during the warmer-than-modern early and mid-Miocene (23 to 14 Myr ago), suggesting a dynamic Antarctic ice sheet (AIS). So far, however, realistic simulations of the Miocene AIS have been limited to equilibrium states under different CO 2 levels and orbital settings. Earlier transient simulations lacked ice-sheet–atmosphere interactions and used a present-day rather than Miocene Antarctic bedrock topography. Here, we quantify the effect of ice-sheet–atmosphere interactions, running the ice-sheet model IMAU-ICE using climate forcing from Miocene simulations by the general circulation model GENESIS. Utilising a recently developed matrix interpolation method enables us to interpolate the climate forcing based on CO 2 levels (between 280 and 840 ppm), as well as varying ice-sheet configurations (between no ice and a large East Antarctic Ice Sheet). We furthermore implement recent reconstructions of Miocene Antarctic bedrock topography. We find that the positive albedo–temperature feedback, partly compensated for by a negative feedback between ice volume and precipitation, increases hysteresis in the relation between CO 2 and ice volume. Together, these ice-sheet–atmosphere interactions decrease the amplitude of Miocene AIS variability in idealised transient simulations. Forced by quasi-orbital 40 kyr forcing CO 2 cycles, the ice volume variability reduces by 21 % when ice-sheet–atmosphere interactions are included compared to when forcing variability is only based on CO 2 changes. Thereby, these interactions also diminish the contribution of AIS variability to benthic δ18 O fluctuations. Evolving bedrock topography during the early and mid-Miocene also reduces ice volume variability by 10 % under equal 40 kyr cycles of atmosphere and ocean forcing. [ABSTRACT FROM AUTHOR]

Accurate estimates of surface mass balance over the Greenland ice sheet (GrIS) would contribute to understanding the cause of recent changes and would help to better estimate the future contribution of the GrIS to sea‐level rise. Given the limitations of in‐situ measurements, modeling, and remote sensing, it is critical to explore the opportunity to merge the available data to better characterize the spatial and temporal variation of the GrIS surface mass balance (SMB). This work utilizes a particle batch smoother data assimilation technique that yields SMB estimates that benefit from the snow model Crocus and a 16‐day albedo product derived from satellite remote sensing data. Comparison of the results against in‐situ SMB measurements shows that the assimilation of the albedo product reduces the root mean square error of the posterior estimates of SMB by 51% and reduces bias by 95%. Plain Language Summary: Greenland ice sheet (GrIS) is losing mass through ice discharge from outlet glaciers and surface processes (e.g., meltwater runoff, sublimation, and evaporation). Recent studies suggest that meltwater runoff will be the dominant mass loss process over the GrIS in the future as it will increase under climate warming. Accurate estimates of the GrIS surface mass balance (SMB) are a critical objective, which, despite its importance, continues to contain large uncertainties from significant errors in forcing data as well as model errors. This work uses a data assimilation framework (which has not been used in estimation of the GrIS SMB) and a satellite‐derived 16‐day albedo product to produce a reanalysis estimates of SMB along the Kangerlussuaq transect (K‐transect) stations in west Greenland. We used the K‐transect in‐situ SMB measurements to validate our results over the 2000–2014 hydrological year. The data assimilation technique (i.e., particle batch smoother) reduces the spatial root‐mean‐square error of SMB over the K‐transect stations by 51% from 858 millimeter water equivalent (mmWE) to 423 mmWE and the bias in the estimates by 95%, from −70 to 3.5 mmWE. It was shown that this methodology has the potential to resolve the spatial variability of the surface processes along the K‐transect stations and in particular of the bare ice surface albedo that is not resolved by the model at a resolution of 25 km (i.e., the model uses a constant bare ice albedo). The results suggest that the methodology can be applied over the entire GrIS using MODIS albedo observations to generate an improved reanalysis of SMB estimates. Key Points: A data assimilation method was used to generate a reanalysis estimates of the surface mass balance of the Greenland ice sheet along the K‐transect stationsA particle batch smoother technique was used to condition the prior estimates of surface mass balance on 16‐day MODIS albedoResults show that the assimilation of albedo reduces the root mean square error of the surface mass balance estimates by 51% [ABSTRACT FROM AUTHOR]

River deltas will likely experience significant land loss because of relative sea‐level rise (RSLR), but predictions have not been tested against observations. Here, we use global data of RSLR and river sediment supply to build a model of delta response to RSLR for 6,402 deltas, representing 86% of global delta land. We validate this model against delta land area change observations from 1985–2015, and project future land area change for IPCC RSLR scenarios. For 2100, we find widely ranging delta scenarios, from +94 ± 125 (2 s.d.) km2 yr−1 for representative concentration pathway (RCP) 2.6 to −1,026 ± 281 km2 yr−1 for RCP8.5. River dams, subsidence, and sea‐level rise have had a comparable influence on reduced delta growth over the past decades, but if we follow RCP8.5 to 2100, more than 85% of delta land loss will be caused by climate‐change driven sea‐level rise, resulting in a loss of ∼5% of global delta land. Plain Language Summary: River deltas can erode and lose land from sea‐level rise. Here we make model predictions of the effects of sea‐level rise for global delta land area change up to 2100. Our model is validated against observations of delta land area change from 1985–2015. For 2100, we find that most climate change scenarios lead to net delta land loss. Worst‐case scenarios for 2100 lead to a global river delta land loss of ∼5% of delta land, at a rate of 1,000 km2 per year. Key Points: Model of river delta response to sea level rise and river sediment applied to 6,402 deltas and tested against land area change observationsProjections suggest net delta land loss for GMSLR >5.5 mm/yr, and for RSLR under IPCC RCP4.5 and RCP8.5 by 2100Climate‐change driven sea‐level rise will likely exceed the effects of dams and subsidence on global delta land loss by 2100 [ABSTRACT FROM AUTHOR]

Improving our confidence in future projections of sea-level rise requires models that can simulate ice-sheet evolution both in the future and in the geological past. A physically accurate treatment of large changes in ice-sheet geometry requires a proper treatment of processes near the margin, like grounding line dynamics, which in turn requires a high spatial resolution in that specific region, so that small-scale topographical features are resolved. This leads to a demand for computationally efficient models, where such a high resolution can be feasibly applied in simulations of 105 – 107 years in duration. Here, we present and evaluate a new ice-sheet model that solves the hybrid SIA–SSA approximation of the stress balance, including a heuristic rule for the grounding-line flux. This is done on a dynamic adaptive mesh which is adapted to the modelled ice-sheet geometry during a simulation. Mesh resolution can be configured to be fine only at specified areas, such as the calving front or the grounding line, as well as specified point locations such as ice-core drill sites. This strongly reduces the number of grid points where the equations need to be solved, increasing the computational efficiency. A high resolution allows the model to resolve small geometrical features, such as outlet glaciers and sub-shelf pinning points, which can significantly affect large-scale ice-sheet dynamics. We show that the model reproduces the analytical solutions or model intercomparison benchmarks for a number of schematic ice-sheet configurations, indicating that the numerical approach is valid. Because of the unstructured triangular mesh, the number of vertices increases less rapidly with resolution than in a square-grid model, greatly reducing the required computation time for high resolutions. A simulation of all four continental ice sheets during an entire 120 kyr glacial cycle, with a 4 km resolution near the grounding line, is expected to take 100–200 wall clock hours on a 16-core system (1600–3200 core hours), implying that this model can be feasibly used for high-resolution palaeo-ice-sheet simulations. [ABSTRACT FROM AUTHOR]

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 δ18 O 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 δ18 O 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. [ABSTRACT FROM AUTHOR]

Improving our confidence in future projections of sea-level rise requires models that can simulate ice-sheet evolution both in the future and in the geological past. A physically accurate treatment of large changes in ice-sheet geometry requires a proper treatment of processes near the margin, like grounding line dynamics, which in turn requires a high spatial resolution in that specific region. This leads to a demand for computationally efficient models, where such a high resolution can be feasibly applied in simulations of 105-107 yr in duration. Here, we present and evaluate a new ice-sheet model that solves the SIA and SSA approximations of the stress balance on a fully adaptive, unstructured triangular mesh. This strongly reduces the number of grid points where the equations need to be solved, increasing the computational efficiency. We show that the model reproduces the analytical solutions or model intercomparison benchmarks for a number of schematic ice-sheet configurations, indicating that the numerical approach is valid. Because of the unstructured triangular mesh, the number of vertices increases less rapidly with resolution than in a square-grid model, greatly reducing the required computation time for high resolutions. A simulation of all four continental ice sheets during an entire 120 kyr glacial cycle, with a 4 km resolution near the grounding line, is expected to take 100-200 wall clock hours on a 16-core system (1,600-3,200 core hours), implying that this model can be feasibly used for high-resolution paleo-ice-sheet simulations. [ABSTRACT FROM AUTHOR]

Future sea-level change projections with process-based stand-alone ice sheet models are typically driven with surface mass balance (SMB) forcing derived from climate models. In this work we address the problems arising from a mismatch of the modelled ice sheet geometry with the geometry used by the climate model. We present a method for applying SMB forcing from climate models to a wide range of Greenland ice sheet models with varying and temporally evolving geometries. In order to achieve that, we translate a given SMB anomaly field as a function of absolute location to a function of surface elevation for 25 regional drainage basins, which can then be applied to different modelled ice sheet geometries. The key feature of the approach is the non-locality of this remapping process. The method reproduces the original forcing data closely when remapped to the original geometry. When remapped to different modelled geometries it produces a physically meaningful forcing with smooth and continuous SMB anomalies across basin divides. The method considerably reduces non-physical biases that would arise by applying the SMB anomaly derived for the climate model geometry directly to a large range of modelled ice sheet model geometries. [ABSTRACT FROM AUTHOR]

Understanding the evolution of, and the interactions between, ice sheets and the global climate over geological time is important for being able to constrain earth system sensitivity. However, direct observational evidence of past CO2 concentrations 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 such as Marine Isotope Stage (MIS) KM5c, to 235-250 ppmv for the MIS M2 glacial excursion. Sea level is relatively stable during this period, with a high stand of 6-14 m, and a drop of 12-26 m during MIS M2. 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. During 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 ppmv to 205-225 ppmv during this period. After the MPT, when the glacial cycles change from 40 kyr 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. [ABSTRACT FROM AUTHOR]

Changes in sea level lead to some of the most severe impacts of anthropogenic climate change. Consequently, they are a subject of great interest in both scientific research and public policy. This paper defines concepts and terminology associated with sea level and sea-level changes in order to facilitate progress in sea-level science, in which communication is sometimes hindered by inconsistent and unclear language. We identify key terms and clarify their physical and mathematical meanings, make links between concepts and across disciplines, draw distinctions where there is ambiguity, and propose new terminology where it is lacking or where existing terminology is confusing. We include formulae and diagrams to support the definitions. [ABSTRACT FROM AUTHOR]

The interaction between ice shelves and the ocean is an important process for the development of marine ice sheets. However, it is difficult to model in full detail due to the high computational cost of coupled ice–ocean simulations, so that simplified basal-melt parameterizations are required. In this work, a new analytical expression for basal melt is derived from the theory of buoyant meltwater plumes moving upward under the ice shelf and driving the overturning circulation within the ice-shelf cavity. The governing equations are nondimensionalized in the case of an ice shelf with constant basal slope and uniform ambient ocean conditions. An asymptotic analysis of these equations in terms of small slopes and small thermal driving, assumed typical for Antarctic ice shelves, leads to an equation that can be solved analytically for the dimensionless melt rate. This analytical expression describes a universal melt-rate curve onto which the scaled results of the original plume model collapse. Its key features are a positive melt peak close to the grounding line and a transition to refreezing further away. Comparing the analytical expression with numerical solutions of the plume model generally shows a close agreement between the two, even for more general cases than the idealized geometry considered in the derivation. The results show how the melt rates adapt naturally to changes in the geometry and ambient ocean temperature. The new expression can readily be used for improving ice-sheet models that currently still lack a sufficiently realistic description of basal melt. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Ice flow forced by gravity is governed by the full Stokes (FS) equations, which are computationally expensive to solve due to the nonlinearity introduced by the rheology. Therefore, approximations to the FS equations are commonly used, especially when modeling a marine ice sheet (ice sheet, ice shelf, and/or ice stream) for 103 years or longer. The shallow ice approximation (SIA) and shallow shelf approximation (SSA) are commonly used but are accurate only for certain parts of an ice sheet. Here, we report a novel way of iteratively coupling FS and SSA that has been implemented in Elmer/Ice and applied to conceptual marine ice sheets. The FS–SSA coupling appears to be very accurate; the relative error in velocity compared to FS is below 0.5 % for diagnostic runs and below 5 % for prognostic runs. Results for grounding line dynamics obtained with the FS–SSA coupling are similar to those obtained from an FS model in an experiment with a periodical temperature forcing over 3000 years that induces grounding line advance and retreat. The rapid convergence of the FS–SSA coupling shows a large potential for reducing computation time, such that modeling a marine ice sheet for thousands of years should become feasible in the near future. Despite inefficient matrix assembly in the current implementation, computation time is reduced by 32 %, when the coupling is applied to a 3-D ice shelf. [ABSTRACT FROM AUTHOR]

Abstract: We reanalyze existing paleodata of global mean surface temperature ΔTg and radiative forcing ΔR of CO2 and land ice albedo for the last 800,000 years to show that a state‐dependency in paleoclimate sensitivity S, as previously suggested, is only found if ΔTg is based on reconstructions, and not when ΔTg is based on model simulations. Furthermore, during times of decreasing obliquity (periods of land ice sheet growth and sea level fall) the multimillennial component of reconstructed ΔTg diverges from CO2, while in simulations both variables vary more synchronously, suggesting that the differences during these times are due to relatively low rates of simulated land ice growth and associated cooling. To produce a reconstruction‐based extrapolation of S for the future, we exclude intervals with strong ΔTg‐CO2 divergence and find that S is less state‐dependent, or even constant state‐independent), yielding a mean equilibrium warming of 2–4 K for a doubling of CO2. [ABSTRACT FROM AUTHOR]

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 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 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. [ABSTRACT FROM AUTHOR]

Observational evidence, including offshore moraines and sediment cores, confirm that at the Last Glacial Maximum (LGM) the Greenland ice sheet (GrIS) expanded to a significantly larger spatial extent than seen at present, grounding into Baffin Bay and out onto the continental shelf break. Given this larger spatial extent and its close proximity to the neighbouring Laurentide Ice Sheet (LIS) and Innuitian Ice Sheet (IIS), it is likely these ice sheets will have had a strong non-local influence on the spatial and temporal behaviour of the GrIS. Most previous paleo ice-sheet modelling simulations recreated an ice sheet that either did not extend out onto the continental shelf or utilized a simplified marine ice parameterization which did not fully include the effect of ice shelves or neglected the sensitivity of the GrIS to this non-local bedrock signal from the surrounding ice sheets. In this paper, we investigated the evolution of the GrIS over the two most recent glacial-interglacial cycles (240 ka BP to the present day) using the ice-sheet-ice-shelf model IMAU-ICE. We investigated the solid earth influence of the LIS and IIS via an offline relative sea level (RSL) forcing Observational evidence, including offshore moraines and sediment cores, confirm that at the Last Glacial Maximum (LGM) the Greenland ice sheet (GrIS) expanded to a significantly larger spatial extent than seen at present, grounding into Baffin Bay and out onto the continental shelf break. Given this larger spatial extent and its close proximity to the neighbouring Laurentide Ice Sheet (LIS) and Innuitian Ice Sheet (IIS), it is likely these ice sheets will have had a strong non-local influence on the spatial and temporal behaviour of the GrIS. Most previous paleo ice-sheet modelling simulations recreated an ice sheet that either did not extend out onto the continental shelf or utilized a simplified marine ice parameterization which did not fully include the effect of ice shelves or neglected the sensitivity of the GrIS to this non-local bedrock signal from the surrounding ice sheets. In this paper, we investigated the evolution of the GrIS over the two most recent glacial-interglacial cycles (240 ka BP to the present day) using the ice-sheet-ice-shelf model IMAU-ICE. We investigated the solid earth influence of the LIS and IIS via an offline relative sea level (RSL) forcing generated by a glacial isostatic adjustment (GIA) model. The RSL forcing governed the spatial and temporal pattern of sub-ice-shelf melting via changes in the water depth below the ice shelves. In the ensemble of simulations, at the glacial maximums, the GrIS coalesced with the IIS to the north and expanded to the continental shelf break to the southwest but remained too restricted to the northeast. In terms of the global mean sea level contribution, at the Last Interglacial (LIG) and LGM the ice sheet added 1.46 and -2.59 m, respectively. This LGM contribution by the GrIS is considerably higher (~1.26 m) than most previous studies whereas the contribution to the LIG highstand is lower (~0.7 m). The spatial and temporal behaviour of the northern margin was highly variable in all simulations, controlled by the sub-ice-shelf melting which was dictated by the RSL forcing and the glacial history of the IIS and LIS. In contrast, the southwestern part of the ice sheet was insensitive to these forcings, with a uniform response in all simulations controlled by the surface air temperature, derived from ice cores. [ABSTRACT FROM AUTHOR]

We evaluate modelled Greenland ice sheet (GrIS) near-surface climate, surface energy balance (SEB) and surface mass balance (SMB) from the updated regional climate model RACMO2 (1958-2016). The new model version, referred to as RACMO2.3p2, incorporates updated glacier outlines, topography and ice albedo fields. Parameters in the cloud scheme governing the conversion of cloud condensate into precipitation have been tuned to correct inland snowfall underestimation: snow properties are modified to reduce drifting snow and melt production in the ice sheet percolation zone. The ice albedo prescribed in the updated model is lower at the ice sheet margins, increasing ice melt locally. RACMO2.3p2 shows good agreement compared to in situ meteorological data and point SEB/SMB measurements, and better resolves the spatial patterns and temporal variability of SMB compared with the previous model version, notably in the north-east, south-east and along the K-transect in south-western Greenland. This new model version provides updated, high-resolution gridded fields of the GrIS presentday climate and SMB, and will be used for projections of the GrIS climate and SMB in response to a future climate scenario in a forthcoming study. [ABSTRACT FROM AUTHOR]

Basal melting below ice shelves is a major factor in mass loss from the Antarctic Ice Sheet, which can contribute significantly to possible future sea-level rise. Therefore, it is important to have an adequate description of the basal melt rates for use in ice-dynamical models. Most current ice models use rather simple parametrizations based on the local balance of heat between ice and ocean. In this work, however, we use a recently derived parametrization of the melt rates based on a buoyant meltwater plume travelling upward beneath an ice shelf. This plume parametrization combines a non-linear ocean temperature sensitivity with an inherent geometry dependence, which is mainly described by the grounding-line depth and the local slope of the ice-shelf base. For the first time, this type of parametrization is evaluated on a two-dimensional grid covering the entire Antarctic continent. In order to apply the essentially one-dimensional parametrization to realistic ice-shelf geometries, we present an algorithm that determines effective values for the grounding-line depth and basal slope in any point beneath an ice shelf. Furthermore, since detailed knowledge of temperatures and circulation patterns in the ice-shelf cavities is sparse or absent, we construct an effective ocean temperature field from observational data with the purpose of matching (area-averaged) melt rates from the model with observed present-day melt rates. Our results qualitatively replicate large-scale observed features in basal melt rates around Antarctica, not only in terms of average values, but also in terms of the spatial pattern, with high melt rates typically occurring near the grounding line. The plume parametrization and the effective temperature field presented here are therefore promising tools for future simulations of the Antarctic Ice Sheet requiring a more realistic oceanic forcing. [ABSTRACT FROM AUTHOR]

Currently a paradigm shift is made from global averaged to spatially variable sea level change (SLC) projections. Traditionally, the contribution from ice sheet mass loss to SLC is considered to be symmetrically distributed. However, several assessments suggest that the probability distribution of dynamical ice sheet mass loss is asymmetrically distributed towards higher SLC values. Here we show how asymmetric probability distributions of dynamical ice sheet mass loss impact the high-end uncertainties of regional SLC projections across the globe. For this purpose we use distributions of dynamical ice sheet mass loss presented by Church et al. (2013), De Vries and Van de Wal (2015) and Ritz et al. (2015). The global average median can be 0.18m higher compared to symmetric distributions based on IPCCAR5, but the change in the global average 95th percentile SLC is considerably larger with a shift of 0.32 m. Locally the 90th, 95th and 97.5th SLC percentiles exceed +1.4, +1.6 and C1.8 m. The high-end percentiles of SLC projections are highly sensitive to the precise shape of the probability distributions of dynamical ice sheet mass loss. The shift towards higher values is of importance for coastal safety strategies as they are based on the high-end percentiles of projections. [ABSTRACT FROM AUTHOR]