Estimating Antarctic Ice Sheet Contributions to Future Sea Level Rise Using a Coupled Climate-Ice Sheet Model

Author(s):  
Jun-Young Park ◽  
Fabian Schloesser ◽  
Axel Timmermann ◽  
Dipayan Choudhury ◽  
June-Yi Lee ◽  
...  

<p>One of the largest uncertainties in projecting future global mean sea level (GSML) rise in response to anthropogenic global warming originates from the Antarctic ice sheet (AIS) contribution. Previous studies suggested that a potential AIS collapse due to the Marine Ice Sheet Instability (MISI) and Marine Ice Cliff Instability (MICI) may contribute up to 1m GMSL rise by the year 2100. However, these estimates were based on uncoupled ice sheet models that do not capture interactions between the AIS and the ocean and atmosphere. Here, we explore future GMSL projections using a three-dimensional coupled climate-ice sheet model (LOVECLIP) that simulates ice sheet dynamics in both hemispheres. The model was forced by increasing CO<sub>2</sub> concentrations following the Shared Socioeconomic Pathway (SSP) 1-1.9, 2-4.5 and 5-8.5 scenarios. Over the next 80 years, the corresponding GMSL contribution from AIS amounts to about 2cm, 8cm and 11cm, respectively. Additional sensitivity experiments show that AIS meltwater flux in response to the SSP 5-8.5 CO<sub>2</sub> concentrations causes subsurface Southern Ocean warming which leads to an additional 20% AIS melting and a reduction in Southern Hemispheric future warming.</p>

2020 ◽  
Author(s):  
Jim Jordan ◽  
Hilmar Gudmundsson ◽  
Adrian Jenkins ◽  
Chris Stokes ◽  
Stewart Jamieson ◽  
...  

<p>The East Antarctic Ice Sheet (EAIS) is the single largest potential contributor to future global mean sea level rise, containing a water mass equivalent of 53 m. Recent work has found the overall mass balance of the EAIS to be approximately in equilibrium, albeit with large uncertainties. However, changes in oceanic conditions have the potential to upset this balance. This could happen by both a general warming of the ocean and also by shifts in oceanic conditions allowing warmer water masses to intrude into ice shelf cavities.</p><p>We use the Úa numerical ice-flow model, combined with ocean-melt rates parameterized by the PICO box mode, to predict the future contribution to global-mean sea level of the EAIS. Results are shown for the next 100 years under a range of emission scenarios and oceanic conditions on a region by region basis, as well as for the whole of the EAIS. </p>


2014 ◽  
Vol 7 (2) ◽  
pp. 1791-1827
Author(s):  
G. Shaffer

Abstract. The Dcess Antarctic Ice Sheet (DAIS) model is presented. Model hindcasts of Antarctic Ice Sheet (AIS) sea level equivalent are forced by reconstructed Antarctic temperatures, global mean sea level and high-latitude, subsurface ocean temperatures, the latter calculated using the Danish Center for Earth System Science (DCESS) Earth System Model forced by reconstructed global mean atmospheric temperatures. The model is calibrated by comparing such hindcasts for different model configurations with paleoreconstructions of AIS sea level equivalent from the last interglacial, the last glacial maximum and the mid-Holocene. The calibrated model is then validated against present estimates of the rate of AIS ice loss. It is found that a high-order dependency of ice flow at the grounding line on water depth there is needed to capture the observed response of the AIS at ice age terminations. Furthermore it is found that a dependency of this ice flow on ocean subsurface temperature by way of ice shelf demise and a resulting buttressing decrease is needed to explain the contribution of the AIS to global mean sea level rise at the last interglacial. When forced and calibrated in this way, model hindcasts of the rate of present day AIS ice loss agree with recent, data-based estimates of this ice loss rate.


2014 ◽  
Vol 7 (4) ◽  
pp. 1803-1818 ◽  
Author(s):  
G. Shaffer

Abstract. The DCESS (Danish Center for Earth System Science) Antarctic Ice Sheet (DAIS) model is presented. Model hindcasts of Antarctic ice sheet (AIS) sea level equivalent are forced by reconstructed Antarctic temperatures, global mean sea level and high-latitude, ocean subsurface temperatures, the latter calculated using the DCESS model forced by reconstructed global mean atmospheric temperatures. The model is calibrated by comparing such hindcasts for different model configurations with paleoreconstructions of AIS sea level equivalent from the last interglacial, the last glacial maximum and the mid-Holocene. The calibrated model is then validated against present estimates of the rate of AIS ice loss. It is found that a high-order dependency of ice flow at the grounding line on water depth there is needed to capture the observed response of the AIS at ice age terminations. Furthermore, it is found that a dependency of this ice flow on ocean subsurface temperature by way of ice shelf demise and a resulting buttressing decrease is needed to explain the contribution of the AIS to global mean sea level rise at the last interglacial. When forced and calibrated in this way, model hindcasts of the rate of present-day AIS ice loss agree with recent, data-based estimates of this ice loss rate.


2018 ◽  
Author(s):  
Javier Blasco ◽  
Ilaria Tabone ◽  
Jorge Alvarez-Solas ◽  
Alexander Robinson ◽  
Marisa Montoya

Abstract. The Antarctic Ice Sheet (AIS) is the largest ice sheet on Earth and hence a major potential contributor to future global sea-level rise. A wealth of studies suggest that increasing oceanic temperatures could cause a collapse of its marine-based western sector, the West Antarctic Ice Sheet, through the mechanism of marine ice-sheet instability, leading to a sea-level increase of 3–5 m. Thus, it is crucial to constrain the sensitivity of the AIS to rapid climate changes. The Last Glacial Period is an ideal benchmark period for this purpose as it was punctuated by abrupt Dansgaard-Oeschger events at millennial timescales. Because their centre of action was in the North Atlantic, where their climate impacts were largest, modelling studies have mainly focused on the millennial-scale evolution of Northern Hemisphere (NH) paleo ice sheets. Sea-level reconstructions attribute the origin of millennial-scale sea-level variations mainly to NH paleo ice sheets, with a minor but not negligible role to the AIS. Here we investigate the AIS response to millennial-scale climate variability for the first time. To this end we use a three-dimensional, thermomechanical hybrid, ice-sheet-shelf model. Different oceanic sensitivities are tested and the sea-level equivalent (SLE) contributions computed. We find that whereas atmospheric variability has no appreciable effect on the AIS, changes in submarine melting rates can have a strong impact on it. We show that in contrast to the widespread assumption that the AIS is a slow reactive and static ice sheet that responds at orbital timescales only, it can lead to ice discharges of almost 15 m of SLE involving substantial grounding line migrations at millennial timescales.


2019 ◽  
Author(s):  
Constantijn J. Berends ◽  
Bas de Boer ◽  
Aisling M. Dolan ◽  
Daniel J. Hill ◽  
Roderik S. W. van de Wal

Abstract. 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.


2021 ◽  
Author(s):  
Emily A. Hill ◽  
Sebastian H. R. Rosier ◽  
G. Hilmar Gudmundsson ◽  
Matthew Collins

Abstract. The future of the Antarctic Ice Sheet in response to climate warming is one of the largest sources of uncertainty in estimates of future changes in global mean sea level (∆GMSL). Mass loss is currently concentrated in regions of warm circumpolar deep water, but it is unclear how ice shelves currently surrounded by relatively cold ocean waters will respond to climatic changes in the future. Studies suggest that warm water could flush the Filchner-Ronne (FR) ice shelf cavity during the 21st century, but the inland ice sheet response to a drastic increase in ice shelf melt rates, is poorly known. Here, we use an ice flow model and uncertainty quantification approach to project the GMSL contribution of the FR basin under RCP emissions scenarios, and assess the forward propagation and proportional contribution of uncertainties in model parameters (related to ice dynamics, and atmospheric/oceanic forcing) on these projections. Our probabilistic projections, derived from an extensive sample of the parameter space using a surrogate model, reveal that the FR basin is unlikely to contribute positively to sea level rise by the 23rd century. This is primarily due to the mitigating effect of increased accumulation with warming, which is capable of suppressing ice loss associated with ocean–driven increases in sub-shelf melt. Mass gain (negative ∆GMSL) from the FR basin increases with warming, but uncertainties in these projections also become larger. In the highest emission scenario RCP 8.5, ∆GMSL is likely to range from −103 to 26 mm, and this large spread can be apportioned predominantly to uncertainties in parameters driving increases in precipitation (30 %) and sub-shelf melting (44 %). There is potential, within the bounds of our input parameter space, for major collapse and retreat of ice streams feeding the FR ice shelf, and a substantial positive contribution to GMSL (up to approx. 300 mm), but we consider such a scenario to be very unlikely. Adopting uncertainty quantification techniques in future studies will help to provide robust estimates of potential sea level rise and further identify target areas for constraining projections.


2019 ◽  
Author(s):  
Anders Levermann ◽  
Ricarda Winkelmann ◽  
Torsten Albrecht ◽  
Heiko Goelzer ◽  
Nicholas R. Golledge ◽  
...  

Abstract. The sea level contribution of the Antarctic ice sheet constitutes a large uncertainty in future sea level projections. Here we apply a linear response theory approach to 16 state-of-the-art ice sheet models to estimate the Antarctic ice sheet contribution from basal ice shelf melting within the 21st century. The purpose of this computation is to estimate the uncertainty that arises from large uncertainty in the external forcing that future warming may exert onto the ice sheet. While ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean, the approach is neglecting a number of processes such as surface mass balance related contributions and mechanisms. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any dampening or self-amplifying processes. This is particularly relevant in situations where an instability is dominating the ice loss. Results obtained here are thus relevant in particular wherever the ice loss is dominated by the forcing as opposed to an internal instability, for example in strong warming scenarios. In order to allow for comparison the methodology was chosen to be exactly the same as in an earlier study (Levermann et al., 2014), but with 16 instead of 5 ice sheet models. We include uncertainty in the atmospheric warming response to carbon emissions (full range of CMIP-5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP-5 models in relation to the global mean surface warming) and the observed range of responses of basal ice shelf melting to oceanic warming outside the ice shelf cavity. This uncertainty in basal ice shelf melting is then convoluted with the linear response functions of each of the 16 ice sheet models to obtain the ice flow response to the individual global warming path. The model median for the observational period from 1992 to 2017 is 9.6 mm with a likely range between 5.2 mm and 20.3 mm compared to the observed sea-level contribution from Antarctica of 7.4 mm with a standard deviation of 3.7 mm (Shepherd et al., 2018). For the so-called business-as-usual warming path, RCP-8.5, we obtain a median contribution of the Antarctic ice sheet to global mean sea-level rise within the 21st century of 17 cm with a likely range (66-percentile around the mean) between 9 cm and 36 cm and a very likely range (90-percentile around the mean) between 6 cm and 59 cm. For the RCP-2.6 warming path which will keep the global mean temperature below two degrees of global warming and is thus consistent with the Paris Climate Agreement yields a median of 13 cm of global mean sea-level contribution. The likely range for the RCP-2.6 scenario is between 7 cm and 25 cm and the very likely range is between 5 cm and 39 cm. The structural uncertainties in the method do not allow an interpretation of any higher uncertainty percentiles. We provide projections for the five Antarctic regions and for each model and each scenario, separately. The rate of sea level contribution is highest under the RCP-8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade with a likely range between 2 cm/dec and 8 cm/dec and a very likely range between 1 cm/dec and 13 cm/dec.


2020 ◽  
Vol 11 (1) ◽  
pp. 35-76 ◽  
Author(s):  
Anders Levermann ◽  
Ricarda Winkelmann ◽  
Torsten Albrecht ◽  
Heiko Goelzer ◽  
Nicholas R. Golledge ◽  
...  

Abstract. The sea level contribution of the Antarctic ice sheet constitutes a large uncertainty in future sea level projections. Here we apply a linear response theory approach to 16 state-of-the-art ice sheet models to estimate the Antarctic ice sheet contribution from basal ice shelf melting within the 21st century. The purpose of this computation is to estimate the uncertainty of Antarctica's future contribution to global sea level rise that arises from large uncertainty in the oceanic forcing and the associated ice shelf melting. Ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean. However, by computing only the sea level contribution in response to ice shelf melting, our study is neglecting a number of processes such as surface-mass-balance-related contributions. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any self-dampening or self-amplifying processes. This is particularly relevant in situations in which an instability is dominating the ice loss. The results obtained here are thus relevant, in particular wherever the ice loss is dominated by the forcing as opposed to an internal instability, for example in strong ocean warming scenarios. In order to allow for comparison the methodology was chosen to be exactly the same as in an earlier study (Levermann et al., 2014) but with 16 instead of 5 ice sheet models. We include uncertainty in the atmospheric warming response to carbon emissions (full range of CMIP5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP5 models in relation to the global mean surface warming), and the observed range of responses of basal ice shelf melting to oceanic warming outside the ice shelf cavity. This uncertainty in basal ice shelf melting is then convoluted with the linear response functions of each of the 16 ice sheet models to obtain the ice flow response to the individual global warming path. The model median for the observational period from 1992 to 2017 of the ice loss due to basal ice shelf melting is 10.2 mm, with a likely range between 5.2 and 21.3 mm. For the same period the Antarctic ice sheet lost mass equivalent to 7.4 mm of global sea level rise, with a standard deviation of 3.7 mm (Shepherd et al., 2018) including all processes, especially surface-mass-balance changes. For the unabated warming path, Representative Concentration Pathway 8.5 (RCP8.5), we obtain a median contribution of the Antarctic ice sheet to global mean sea level rise from basal ice shelf melting within the 21st century of 17 cm, with a likely range (66th percentile around the mean) between 9 and 36 cm and a very likely range (90th percentile around the mean) between 6 and 58 cm. For the RCP2.6 warming path, which will keep the global mean temperature below 2 ∘C of global warming and is thus consistent with the Paris Climate Agreement, the procedure yields a median of 13 cm of global mean sea level contribution. The likely range for the RCP2.6 scenario is between 7 and 24 cm, and the very likely range is between 4 and 37 cm. The structural uncertainties in the method do not allow for an interpretation of any higher uncertainty percentiles. We provide projections for the five Antarctic regions and for each model and each scenario separately. The rate of sea level contribution is highest under the RCP8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade, with a likely range between 2 and 9 cm per decade and a very likely range between 1 and 14 cm per decade.


Sign in / Sign up

Export Citation Format

Share Document