Refine
Year of publication
Document Type
- Article (31)
- Postprint (2)
- Doctoral Thesis (1)
Language
- English (34)
Is part of the Bibliography
- yes (34)
Keywords
- Amazon rainforest (2)
- complex networks (2)
- droughts (2)
- prediction (2)
- 21st-Century (1)
- Anthropocene (1)
- Earth System trajectories (1)
- Greenland (1)
- anthropocene (1)
- biosphere feedbacks (1)
In recent decades, the Greenland Ice Sheet has been losing mass and has thereby contributed to global sea-level rise. The rate of ice loss is highly relevant for coastal protection worldwide. The ice loss is likely to increase under future warming. Beyond a critical temperature threshold, a meltdown of the Greenland Ice Sheet is induced by the self-enforcing feedback between its lowering surface elevation and its increasing surface mass loss: the more ice that is lost, the lower the ice surface and the warmer the surface air temperature, which fosters further melting and ice loss. The computation of this rate so far relies on complex numerical models which are the appropriate tools for capturing the complexity of the problem. By contrast we aim here at gaining a conceptual understanding by deriving a purposefully simple equation for the self-enforcing feedback which is then used to estimate the melt time for different levels of warming using three observable characteristics of the ice sheet itself and its surroundings. The analysis is purely conceptual in nature. It is missing important processes like ice dynamics for it to be useful for applications to sea-level rise on centennial timescales, but if the volume loss is dominated by the feedback, the resulting logarithmic equation unifies existing numerical simulations and shows that the melt time depends strongly on the level of warming with a critical slow-down near the threshold: the median time to lose 10% of the present-day ice volume varies between about 3500 years for a temperature level of 0.5 degrees C above the threshold and 500 years for 5 degrees C. Unless future observations show a significantly higher melting sensitivity than currently observed, a complete meltdown is unlikely within the next 2000 years without significant ice-dynamical contributions.
In recent decades, the Greenland Ice Sheet has been losing mass and has thereby contributed to global sea-level rise. The rate of ice loss is highly relevant for coastal protection worldwide. The ice loss is likely to increase under future warming. Beyond a critical temperature threshold, a meltdown of the Greenland Ice Sheet is induced by the self-enforcing feedback between its lowering surface elevation and its increasing surface mass loss: the more ice that is lost, the lower the ice surface and the warmer the surface air temperature, which fosters further melting and ice loss. The computation of this rate so far relies on complex numerical models which are the appropriate tools for capturing the complexity of the problem. By contrast we aim here at gaining a conceptual understanding by deriving a purposefully simple equation for the self-enforcing feedback which is then used to estimate the melt time for different levels of warming using three observable characteristics of the ice sheet itself and its surroundings. The analysis is purely conceptual in nature. It is missing important processes like ice dynamics for it to be useful for applications to sea-level rise on centennial timescales, but if the volume loss is dominated by the feedback, the resulting logarithmic equation unifies existing numerical simulations and shows that the melt time depends strongly on the level of warming with a critical slow-down near the threshold: the median time to lose 10% of the present-day ice volume varies between about 3500 years for a temperature level of 0.5 degrees C above the threshold and 500 years for 5 degrees C. Unless future observations show a significantly higher melting sensitivity than currently observed, a complete meltdown is unlikely within the next 2000 years without significant ice-dynamical contributions.
Droughts in tropical South America have an imminent and severe impact on the Amazon rainforest and affect the livelihoods of millions of people. Extremely dry conditions in Amazonia have been previously linked to sea surface temperature (SST) anomalies in the adjacent tropical oceans. Although the sources and impacts of such droughts have been widely studied, establishing reliable multi-year lead statistical forecasts of their occurrence is still an ongoing challenge. Here, we further investigate the relationship between SST and rainfall anomalies using a complex network approach. We identify four ocean regions which exhibit the strongest overall SST correlations with central Amazon rainfall, including two particularly prominent regions in the northern and southern tropical Atlantic. Based on the time-dependent correlation between SST anomalies in these two regions alone, we establish a new early-warning method for droughts in the central Amazon basin and demonstrate its robustness in hindcasting past major drought events with lead-times up to 18 months.
Droughts in tropical South America have an imminent and severe impact on the Amazon rainforest and affect the livelihoods of millions of people. Extremely dry conditions in Amazonia have been previously linked to sea surface temperature (SST) anomalies in the adjacent tropical oceans. Although the sources and impacts of such droughts have been widely studied, establishing reliable multi-year lead statistical forecasts of their occurrence is still an ongoing challenge. Here, we further investigate the relationship between SST and rainfall anomalies using a complex network approach. We identify four ocean regions which exhibit the strongest overall SST correlations with central Amazon rainfall, including two particularly prominent regions in the northern and southern tropical Atlantic. Based on the time-dependent correlation between SST anomalies in these two regions alone, we establish a new early-warning method for droughts in the central Amazon basin and demonstrate its robustness in hindcasting past major drought events with lead-times up to 18 months.
Ocean-induced melting below ice shelves is one of the dominant drivers for mass loss from the Antarctic Ice Sheet at present. An appropriate representation of sub-shelf melt rates is therefore essential for model simulations of marine-based ice sheet evolution. Continental-scale ice sheet models often rely on simple melt-parameterizations, in particular for long-term simulations, when fully coupled ice-ocean interaction becomes computationally too expensive. Such parameterizations can account for the influence of the local depth of the ice-shelf draft or its slope on melting. However, they do not capture the effect of ocean circulation underneath the ice shelf. Here we present the Potsdam Ice-shelf Cavity mOdel (PICO), which simulates the vertical overturning circulation in ice-shelf cavities and thus enables the computation of sub-shelf melt rates consistent with this circulation. PICO is based on an ocean box model that coarsely resolves ice shelf cavities and uses a boundary layer melt formulation. We implement it as a module of the Parallel Ice Sheet Model (PISM) and evaluate its performance under present-day conditions of the Southern Ocean. We identify a set of parameters that yield two-dimensional melt rate fields that qualitatively reproduce the typical pattern of comparably high melting near the grounding line and lower melting or refreezing towards the calving front. PICO captures the wide range of melt rates observed for Antarctic ice shelves, with an average of about 0.1 ma(-1) for cold sub-shelf cavities, for example, underneath Ross or Ronne ice shelves, to 16 ma(-1) for warm cavities such as in the Amundsen Sea region. This makes PICO a computationally feasible and more physical alternative to melt parameterizations purely based on ice draft geometry.
Consequences of twenty-first-century policy for multi-millennial climate
and sea-level change
(2016)
Most of the policy debate surrounding the actions needed to mitigate and adapt to anthropogenic climate change has been framed by observations of the past 150 years as well as climate and sea-level projections for the twenty-first century. The focus on this 250-year window, however, obscures some of the most profound problems associated with climate change. Here, we argue that the twentieth and twenty-first centuries, a period during which the overwhelming majority of human-caused carbon emissions are likely to occur, need to be placed into a long-term context that includes the past 20 millennia, when the last Ice Age ended and human civilization developed, and the next ten millennia, over which time the projected impacts of anthropogenic climate change will grow and persist. This long-term perspective illustrates that policy decisions made in the next few years to decades will have profound impacts on global climate, ecosystems and human societies - not just for this century, but for the next ten millennia and beyond.
Projections of changes in Antarctic Ice Sheet (AIS) surface mass balance indicate a negative contribution to sea level because of the expected increase in precipitation due to the higher moisture holding capacity of warmer air(1). Observations over the past decades, however, are unable to constrain the relation between temperature and accumulation changes because both are dominated by strong natural variability(2-5). Here we derive a consistent continental-scale increase in accumulation of approximately 5 +/- 1% K-1, through the assessment of ice-core data (spanning the large temperature change during the last deglaciation, 21,000 to 10,000 years ago), in combination with palaeo-simulations, future projections by 35 general circulation models (GCMs), and one high-resolution future simulation. The ice-core data and modelling results for the last deglaciation agree, showing uniform local sensitivities of similar to 6% K-1. The palaeo-simulation allows for a continental-scale aggregation of accumulation changes reaching 4.3% K-1. Despite the different timescales, these sensitivities agree with the multi-model mean of 6.1 +/- 2.6% K-1 (GCMprojections) and the continental-scale sensitivity of 4.9% K-1 (high-resolution future simulation). Because some of the mass gain of the AIS is offset by dynamical losses induced by accumulation(6,7), we provide a response function allowing projections of sea-level fall in terms of continental-scale accumulation changes that compete with surface melting and dynamical losses induced by other mechanisms(6,8,9).
The past and future evolution of the Antarctic Ice Sheet is largely controlled by interactions between the ocean and floating ice shelves. To investigate these interactions, coupled ocean and ice sheet model configurations are required. Previous modelling studies have mostly relied on high-resolution configurations, limiting these studies to individual glaciers or regions over short timescales of decades to a few centuries. We present a framework to couple the dynamic ice sheet model PISM (Parallel Ice Sheet Model) with the global ocean general circulation model MOM5 (Modular Ocean Model) via the ice shelf cavity model PICO (Pots-dam Ice-shelf Cavity mOdel). As ice shelf cavities are not resolved by MOM5 but are parameterized with the PICO box model, the framework allows the ice sheet and ocean components to be run at resolutions of 16 km and 3 degrees respectively. This approach makes the coupled configuration a useful tool for the analysis of interactions between the Antarctic Ice Sheet and the global ocean over time spans of the order of centuries to millennia. In this study, we describe the technical implementation of this coupling framework: sub-shelf melting in the ice sheet component is calculated by PICO from modelled ocean temperatures and salinities at the depth of the continental shelf, and, vice versa, the resulting mass and energy fluxes from melting at the ice-ocean interface are transferred to the ocean component. Mass and energy fluxes are shown to be conserved to machine precision across the considered component domains. The implementation is computationally efficient as it introduces only minimal overhead. Furthermore, the coupled model is evaluated in a 4000 year simulation under constant present-day climate forcing and is found to be stable with respect to the ocean and ice sheet spin-up states. The framework deals with heterogeneous spatial grid geometries, varying grid resolutions, and timescales between the ice and ocean component in a generic way; thus, it can be adopted to a wide range of model set-ups.
The past rapid growth of Northern Hemisphere continental ice sheets, which terminated warm and stable climate periods, is generally attributed to reduced summer insolation in boreal latitudes(1-3). Yet such summer insolation is near to its minimum at present(4), and there are no signs of a new ice age(5). This challenges our understanding of the mechanisms driving glacial cycles and our ability to predict the next glacial inception(6). Here we propose a critical functional relationship between boreal summer insolation and global carbon dioxide (CO2) concentration, which explains the beginning of the past eight glacial cycles and might anticipate future periods of glacial inception. Using an ensemble of simulations generated by an Earth system model of intermediate complexity constrained by palaeoclimatic data, we suggest that glacial inception was narrowly missed before the beginning of the Industrial Revolution. The missed inception can be accounted for by the combined effect of relatively high late-Holocene CO2 concentrations and the low orbital eccentricity of the Earth(7). Additionally, our analysis suggests that even in the absence of human perturbations no substantial build-up of ice sheets would occur within the next several thousand years and that the current interglacial would probably last for another 50,000 years. However, moderate anthropogenic cumulative CO2 emissions of 1,000 to 1,500 gigatonnes of carbon will postpone the next glacial inception by at least 100,000 years(8,9). Our simulations demonstrate that under natural conditions alone the Earth system would be expected to remain in the present delicately balanced interglacial climate state, steering clear of both large-scale glaciation of the Northern Hemisphere and its complete deglaciation, for an unusually long time.
The stability of the Greenland Ice Sheet under global warming is governed by a number of dynamic processes and interacting feedback mechanisms in the ice sheet, atmosphere and solid Earth.
Here we study the long-term effects due to the interplay of the competing melt-elevation and glacial isostatic adjustment (GIA) feedbacks for different temperature step forcing experiments with a coupled ice-sheet and solid-Earth model.
Our model results show that for warming levels above 2 degrees C, Greenland could become essentially ice-free within several millennia, mainly as a result of surface melting and acceleration of ice flow. These ice losses are mitigated, however, in some cases with strong GIA feedback even promoting an incomplete recovery of the Greenland ice volume. We further explore the full-factorial parameter space determining the relative strengths of the two feedbacks: our findings suggest distinct dynamic regimes of the Greenland Ice Sheets on the route to destabilization under global warming - from incomplete recovery, via quasi-periodic oscillations in ice volume to ice-sheet collapse.
In the incomplete recovery regime, the initial ice loss due to warming is essentially reversed within 50 000 years, and the ice volume stabilizes at 61 %-93 % of the present-day volume. For certain combinations of temperature increase, atmospheric lapse rate and mantle viscosity, the interaction of the GIA feedback and the melt-elevation feedback leads to self-sustained, long-term oscillations in ice-sheet volume with oscillation periods between 74 000 and over 300 000 years and oscillation amplitudes between 15 %-70 % of present-day ice volume.
This oscillatory regime reveals a possible mode of internal climatic variability in the Earth system on timescales on the order of 100 000 years that may be excited by or synchronized with orbital forcing or interact with glacial cycles and other slow modes of variability. Our findings are not meant as scenario-based near-term projections of ice losses but rather providing insight into of the feedback loops governing the "deep future" and, thus, long-term resilience of the Greenland Ice Sheet.
Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM)
(2020)
The Parallel Ice Sheet Model (PISM) is applied to the Antarctic Ice Sheet over the last two glacial cycles (approximate to 210 000 years) with a resolution of 16 km. An ensemble of 256 model runs is analyzed in which four relevant model parameters have been systematically varied using full-factorial parameter sampling. Parameters and plausible parameter ranges have been identified in a companion paper (Albrecht et al., 2020) and are associated with ice dynamics, climatic forcing, basal sliding and bed deformation and represent distinct classes of model uncertainties. The model is scored against both modern and geologic data, including reconstructed grounding-line locations, elevation-age data, ice thickness, surface velocities and uplift rates. An aggregated score is computed for each ensemble member that measures the overall model-data misfit, including measurement uncertainty in terms of a Gaussian error model (Briggs and Tarasov, 2013). The statistical method used to analyze the ensemble simulation results follows closely the simple averaging method described in Pollard et al. (2016).
This analysis reveals clusters of best-fit parameter combinations, and hence a likely range of relevant model and boundary parameters, rather than individual best-fit parameters. The ensemble of reconstructed histories of Antarctic Ice Sheet volumes provides a score-weighted likely range of sea-level contributions since the Last Glacial Maximum (LGM) of 9.4 +/- 4.1m (or 6.5 +/- 2.0 x 10(6) km(3)), which is at the upper range of most previous studies. The last deglaciation occurs in all ensemble simulations after around 12 000 years before present and hence after the meltwater pulse 1A (MWP1a). Our ensemble analysis also provides an estimate of parametric uncertainty bounds for the present-day state that can be used for PISM projections of future sea-level contributions from the Antarctic Ice Sheet.
Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM)
(2020)
Simulations of the glacial-interglacial history of the Antarctic Ice Sheet provide insights into dynamic threshold behavior and estimates of the ice sheet's contributions to global sea-level changes for the past, present and future. However, boundary conditions are weakly constrained, in particular at the interface of the ice sheet and the bedrock. Also climatic forcing covering the last glacial cycles is uncertain, as it is based on sparse proxy data. <br /> We use the Parallel Ice Sheet Model (PISM) to investigate the dynamic effects of different choices of input data, e.g., for modern basal heat flux or reconstructions of past changes of sea level and surface temperature. As computational resources are limited, glacial-cycle simulations are performed using a comparably coarse model grid of 16 km and various parameterizations, e.g., for basal sliding, iceberg calving, or for past variations in precipitation and ocean temperatures. In this study we evaluate the model's transient sensitivity to corresponding parameter choices and to different boundary conditions over the last two glacial cycles and provide estimates of involved uncertainties. We also discuss isolated and combined effects of climate and sea-level forcing. Hence, this study serves as a "cookbook" for the growing community of PISM users and paleo-ice sheet modelers in general. <br /> For each of the different model uncertainties with regard to climatic forcing, ice and Earth dynamics, and basal processes, we select one representative model parameter that captures relevant uncertainties and motivates corresponding parameter ranges that bound the observed ice volume at present. The four selected parameters are systematically varied in a parameter ensemble analysis, which is described in a companion paper.
Several large-scale cryosphere elements such as the Arctic summer sea ice, the mountain glaciers, the Greenland and West Antarctic Ice Sheet have changed substantially during the last century due to anthropogenic global warming. However, the impacts of their possible future disintegration on global mean temperature (GMT) and climate feedbacks have not yet been comprehensively evaluated. Here, we quantify this response using an Earth system model of intermediate complexity. Overall, we find a median additional global warming of 0.43 degrees C (interquartile range: 0.39-0.46 degrees C) at a CO2 concentration of 400 ppm. Most of this response (55%) is caused by albedo changes, but lapse rate together with water vapour (30%) and cloud feedbacks (15%) also contribute significantly. While a decay of the ice sheets would occur on centennial to millennial time scales, the Arctic might become ice-free during summer within the 21st century. Our findings imply an additional increase of the GMT on intermediate to long time scales. The disintegration of cryosphere elements such as the Arctic summer sea ice, mountain glaciers, Greenland and West Antarctica is associated with temperature and radiative feedbacks. In this work, the authors quantify these feedbacks and find an additional global warming of 0.43 degrees C.
Currently, several large-scale ice-flow models impose a condition on ice flux across grounding lines using an analytically motivated parameterisation of grounding-line flux. It has been suggested that employing this analytical expression alleviates the need for highly resolved computational domains around grounding lines of marine ice sheets. While the analytical flux formula is expected to be accurate in an unbuttressed flow-line setting, its validity has hitherto not been assessed for complex and realistic geometries such as those of the Antarctic Ice Sheet. Here the accuracy of this analytical flux formula is tested against an optimised ice flow model that uses a highly resolved computational mesh around the Antarctic grounding lines. We find that when applied to the Antarctic Ice Sheet the analytical expression provides inaccurate estimates of ice fluxes for almost all grounding lines. Furthermore, in many instances direct application of the analytical formula gives rise to unphysical complex-valued ice fluxes. We conclude that grounding lines of the Antarctic Ice Sheet are, in general, too highly buttressed for the analytical parameterisation to be of practical value for the calculation of grounding-line fluxes.
Higher resilience to climatic disturbances in tropical vegetation exposed to more variable rainfall
(2019)
With ongoing global warming, the amount and frequency of precipitation in the tropics is projected to change substantially. While it has been shown that tropical forests and savannahs are sustained within the same intermediate mean annual precipitation range, the mechanisms that lead to the resilience of these ecosystems are still not fully understood. In particular, the long-term impact of rainfall variability on resilience is as yet unclear. Here we present observational evidence that both tropical forest and savannah exposed to a higher rainfall variability-in particular on interannual scales-during their long-term past are overall more resilient against climatic disturbances. Based on precipitation and tree cover data in the Brazilian Amazon basin, we constructed potential landscapes that enable us to systematically measure the resilience of the different ecosystems. Additionally, we infer that shifts from forest to savannah due to decreasing precipitation in the future are more likely to occur in regions with a precursory lower rainfall variability. Long-term rainfall variability thus needs to be taken into account in resilience analyses and projections of vegetation response to climate change.
Surface melting of the Greenland Ice Sheet contributes a large amount to current and future sea level rise. Increased surface melt may lower the reflectivity of the ice sheet surface and thereby increase melt rates: the so-called melt-albedo feedback describes this self-sustaining increase in surface melting. In order to test the effect of the melt-albedo feedback in a prognostic ice sheet model, we implement dEBM-simple, a simplified version of the diurnal Energy Balance Model dEBM, in the Parallel Ice Sheet Model (PISM). The implementation includes a simple representation of the melt-albedo feedback and can thereby replace the positive-degree-day melt scheme. Using PISM-dEBM-simple, we find that this feedback increases ice loss through surface warming by 60 % until 2300 for the high-emission scenario RCP8.5 when compared to a scenario in which the albedo remains constant at its present-day values. With an increase of 90 % compared to a fixed-albedo scenario, the effect is more pronounced for lower surface warming under RCP2.6. Furthermore, assuming an immediate darkening of the ice surface over all summer months, we estimate an upper bound for this effect to be 70 % in the RCP8.5 scenario and a more than 4-fold increase under RCP2.6. With dEBM-simple implemented in PISM, we find that the melt-albedo feedback is an essential contributor to mass loss in dynamic simulations of the Greenland Ice Sheet under future warming.
Anthropogenic climate change is likely to cause continuing global sea level rise(1), but some processes within the Earth system may mitigate the magnitude of the projected effect. Regional and global climate models simulate enhanced snowfall over Antarctica, which would provide a direct offset of the future contribution to global sea level rise from cryospheric mass loss(2,3) and ocean expansion(4). Uncertainties exist in modelled snowfall(5), but even larger uncertainties exist in the potential changes of dynamic ice discharge from Antarctica(1,6) and thus in the ultimate fate of the precipitation-deposited ice mass. Here we show that snowfall and discharge are not independent, but that future ice discharge will increase by up to three times as a result of additional snowfall under global warming. Our results, based on an ice-sheet model(7) forced by climate simulations through to the end of 2500 (ref. 8), show that the enhanced discharge effect exceeds the effect of surface warming as well as that of basal ice-shelf melting, and is due to the difference in surface elevation change caused by snowfall on grounded versus floating ice. Although different underlying forcings drive ice loss from basal melting versus increased snowfall, similar ice dynamical processes are nonetheless at work in both; therefore results are relatively independent of the specific representation of the transition zone. In an ensemble of simulations designed to capture ice-physics uncertainty, the additional dynamic ice loss along the coastline compensates between 30 and 65 per cent of the ice gain due to enhanced snowfall over the entire continent. This results in a dynamic ice loss of up to 1.25 metres in the year 2500 for the strongest warming scenario. The reported effect thus strongly counters a potential negative contribution to global sea level by the Antarctic Ice Sheet.
initMIP-Antarctica
(2019)
Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMlP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMlP-Greenland, initMlP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMlP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.
ISMIP6 Antarctica
(2020)
Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and assess the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimates of the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes, forcings employed and initial states of ice sheet models. This study presents results from ice flow model simulations from 13 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015-2100 as part of the Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). They are forced with outputs from a subset of models from the Coupled Model Intercomparison Project Phase 5 (CMIP5), representative of the spread in climate model results. Simulations of the Antarctic ice sheet contribution to sea level rise in response to increased warming during this period varies between 7:8 and 30.0 cm of sea level equivalent (SLE) under Representative Concentration Pathway (RCP) 8.5 scenario forcing. These numbers are relative to a control experiment with constant climate conditions and should therefore be added to the mass loss contribution under climate conditions similar to present-day conditions over the same period. The simulated evolution of the West Antarctic ice sheet varies widely among models, with an overall mass loss, up to 18.0 cm SLE, in response to changes in oceanic conditions. East Antarctica mass change varies between 6 :1 and 8.3 cm SLE in the simulations, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional simulated mass loss of 28mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the climate forcing, the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 climate models show an additional mass loss of 0 and 3 cm of SLE on average compared to simulations done under present-day conditions for the two CMIP5 forcings used and display limited mass gain in East Antarctica.
Recently observed large-scale disintegration of Antarctic ice shelves has moved their fronts closer towards grounded ice. In response, ice-sheet discharge into the ocean has accelerated, contributing to global sea-level rise and emphasizing the importance of calving-front dynamics. The position of the ice front strongly influences the stress field within the entire sheet-shelf-system and thereby the mass flow across the grounding line. While theories for an advance of the ice-front are readily available, no general rule exists for its retreat, making it difficult to incorporate the retreat in predictive models. Here we extract the first-order large-scale kinematic contribution to calving which is consistent with large-scale observation. We emphasize that the proposed equation does not constitute a comprehensive calving law but represents the first-order kinematic contribution which can and should be complemented by higher order contributions as well as the influence of potentially heterogeneous material properties of the ice. When applied as a calving law, the equation naturally incorporates the stabilizing effect of pinning points and inhibits ice shelf growth outside of embayments. It depends only on local ice properties which are, however, determined by the full topography of the ice shelf. In numerical simulations the parameterization reproduces multiple stable fronts as observed for the Larsen A and B Ice Shelves including abrupt transitions between them which may be caused by localized ice weaknesses. We also find multiple stable states of the Ross Ice Shelf at the gateway of the West Antarctic Ice Sheet with back stresses onto the sheet reduced by up to 90 % compared to the present state.