Exploring parameter uncertainty in a model of the Antarctic Ice Sheet

Author(s):  
Steven Phipps ◽  
Jason Roberts ◽  
Matt King

<p>Physical processes within ice sheet models are sometimes described by simplified schemes known as parameterisations. The values of the parameters within these schemes can be poorly constrained by theory or observation. Uncertainty in the parameter values translates into uncertainty in the outputs of the models. Proper quantification of the uncertainty in model predictions therefore requires a systematic approach for sampling parameter space. We demonstrate a simple and efficient approach to identify regions of multi-dimensional parameter space that are consistent with observations. Using the Parallel Ice Sheet Model to simulate the present-day state of the Antarctic Ice Sheet, we find that co-dependencies between parameters preclude the identification of a single optimal set of parameter values. Approaches such as large ensemble modelling are therefore required in order to generate model predictions, such as projections of future global sea level rise, that incorporate proper quantification of the uncertainty arising from the parameterisation of physical processes.</p>

2020 ◽  
Author(s):  
Steven J. Phipps ◽  
Jason L. Roberts ◽  
Matt A. King

Abstract. Physical processes within geoscientific models are sometimes described by simplified schemes known as parameterisations. The values of the parameters within these schemes can be poorly constrained by theory or observation. Uncertainty in the parameter values translates into uncertainty in the outputs of the models. Proper quantification of the uncertainty in model predictions therefore requires a systematic approach for sampling parameter space. In this study, we develop a simple and efficient approach to identify regions of multi-dimensional parameter space that are consistent with observations. Using the Parallel Ice Sheet Model to simulate the present-day state of the Antarctic Ice Sheet, we find that co-dependencies between parameters preclude the identification of a single optimal set of parameter values. Approaches such as large ensemble modelling are therefore required in order to generate model predictions that incorporate proper quantification of the uncertainty arising from the parameterisation of physical processes.


2021 ◽  
Vol 14 (8) ◽  
pp. 5107-5124
Author(s):  
Steven J. Phipps ◽  
Jason L. Roberts ◽  
Matt A. King

Abstract. Physical processes within geoscientific models are sometimes described by simplified schemes known as parameterisations. The values of the parameters within these schemes can be poorly constrained by theory or observation. Uncertainty in the parameter values translates into uncertainty in the outputs of the models. Proper quantification of the uncertainty in model predictions therefore requires a systematic approach for sampling parameter space. In this study, we develop a simple and efficient approach to identify regions of multi-dimensional parameter space that are consistent with observations. Using the Parallel Ice Sheet Model to simulate the present-day state of the Antarctic Ice Sheet, we find that co-dependencies between parameters preclude any simple identification of a single optimal set of parameter values. Approaches such as large ensemble modelling are therefore required in order to generate model predictions that incorporate proper quantification of the uncertainty arising from the parameterisation of physical processes.


2019 ◽  
Vol 13 (10) ◽  
pp. 2615-2631 ◽  
Author(s):  
Michelle Tigchelaar ◽  
Axel Timmermann ◽  
Tobias Friedrich ◽  
Malte Heinemann ◽  
David Pollard

Abstract. Antarctic ice volume has varied substantially during the late Quaternary, with reconstructions suggesting a glacial ice sheet extending to the continental shelf break and interglacial sea level highstands of several meters. Throughout this period, changes in the Antarctic Ice Sheet were driven by changes in atmospheric and oceanic conditions and global sea level; yet, so far modeling studies have not addressed which of these environmental forcings dominate and how they interact in the dynamical ice sheet response. Here, we force an Antarctic Ice Sheet model with global sea level reconstructions and transient, spatially explicit boundary conditions from a 408 ka climate model simulation, not only in concert with each other but, for the first time, also separately. We find that together these forcings drive glacial–interglacial ice volume changes of 12–14 ms.l.e., in line with reconstructions and previous modeling studies. None of the individual drivers – atmospheric temperature and precipitation, ocean temperatures, or sea level – single-handedly explains the full ice sheet response. In fact, the sum of the individual ice volume changes amounts to less than half of the full ice volume response, indicating the existence of strong nonlinearities and forcing synergy. Both sea level and atmospheric forcing are necessary to create full glacial ice sheet growth, whereas the contribution of ocean melt changes is found to be more a function of ice sheet geometry than climatic change. Our results highlight the importance of accurately representing the relative timing of forcings of past ice sheet simulations and underscore the need for developing coupled climate–ice sheet modeling frameworks that properly capture key feedbacks.


2015 ◽  
Vol 1 (8) ◽  
pp. e1500589 ◽  
Author(s):  
Ricarda Winkelmann ◽  
Anders Levermann ◽  
Andy Ridgwell ◽  
Ken Caldeira

The Antarctic Ice Sheet stores water equivalent to 58 m in global sea-level rise. We show in simulations using the Parallel Ice Sheet Model that burning the currently attainable fossil fuel resources is sufficient to eliminate the ice sheet. With cumulative fossil fuel emissions of 10,000 gigatonnes of carbon (GtC), Antarctica is projected to become almost ice-free with an average contribution to sea-level rise exceeding 3 m per century during the first millennium. Consistent with recent observations and simulations, the West Antarctic Ice Sheet becomes unstable with 600 to 800 GtC of additional carbon emissions. Beyond this additional carbon release, the destabilization of ice basins in both West and East Antarctica results in a threshold increase in global sea level. Unabated carbon emissions thus threaten the Antarctic Ice Sheet in its entirety with associated sea-level rise that far exceeds that of all other possible sources.


2020 ◽  
Author(s):  
Torsten Albrecht ◽  
Ricarda Winkelmann ◽  
Anders Levermann

<p>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. We use the Parallel Ice Sheet Model (PISM) to investigate the dynamic effects of different choices of input data and of various parameterizations on the sea-level relevant ice volume. 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 present isolated and combined effects of climate and sea-level forcing on glacial time scales. </p>


2020 ◽  
Author(s):  
James O'Neill ◽  
Tamsin Edwards ◽  
Lauren Gregoire ◽  
Niall Gandy ◽  
Aisling Dolan ◽  
...  

<p>The Antarctic ice sheet is a deeply uncertain component of future sea level under anthropogenic climate change. To shed light on the ice sheets response to warmer climates in the past and its’ response to future warming, periods in Earth’s geological record can serve as instructive modelling targets. The mid-Pliocene warm period (3.3 – 3.0 Ma) is characterised by global mean surface temperatures ~2.7-4<sup>o</sup>C above pre-industrial, atmospheric CO<sub>2</sub> concentrations of ~400ppm and eustatic sea level rise on the order of ~10-30m above modern. The mid-Pliocene sea level record is subject to large uncertainties. The upper end of this record implies a significant contribution from Antarctica and possible collapse of regions of the ice sheet, driven by marine ice sheet instabilities.</p><p>We present a suite of BISICLES ice sheet model simulations, forced with a subset of Pliocene Modelling Intercomparison Project (PlioMIP phase 1) coupled atmosphere-ocean climate models, that represent the Pliocene Antarctic ice sheet. This ensemble captures a range of possible ice sheet model responses to a warm Pliocene-like climate under different parameter choices, sampled in a Latin hypercube design. Modelled Antarctic sea level contribution is compared to reconstructions of Pliocene sea level, to explore the extent to which available data with large uncertainties can constrain the model parameter values.</p><p>Our aim with this work is to provide insights on Antarctic contribution to sea level in the warm mid-Pliocene. We seek to characterise the role of ice-ocean, ice-atmosphere and ice-bedrock parameter uncertainty in BISICLES on the ice sheet sea level contribution range, and whether cliff instability processes are necessary in reproduce high Pliocene sea levels in this ice sheet model.</p>


2020 ◽  
Vol 11 (4) ◽  
pp. 953-976
Author(s):  
Jonas Van Breedam ◽  
Heiko Goelzer ◽  
Philippe Huybrechts

Abstract. The emphasis for informing policy makers on future sea-level rise has been on projections by the end of the 21st century. However, due to the long lifetime of atmospheric CO2, the thermal inertia of the climate system and the slow equilibration of the ice sheets, global sea level will continue to rise on a multi-millennial timescale even when anthropogenic CO2 emissions cease completely during the coming decades to centuries. Here we present global sea-level change projections due to the melting of land ice combined with steric sea effects during the next 10 000 years calculated in a fully interactive way with the Earth system model of intermediate complexity LOVECLIMv1.3. The greenhouse forcing is based on the Extended Concentration Pathways defined until 2300 CE with no carbon dioxide emissions thereafter, equivalent to a cumulative CO2 release of between 460 and 5300 GtC. We performed one additional experiment for the highest-forcing scenario with the inclusion of a methane emission feedback where methane is slowly released due to a strong increase in surface and oceanic temperatures. After 10 000 years, the sea-level change rate drops below 0.05 m per century and a semi-equilibrated state is reached. The Greenland ice sheet is found to nearly disappear for all forcing scenarios. The Antarctic ice sheet contributes only about 1.6 m to sea level for the lowest forcing scenario with a limited retreat of the grounding line in West Antarctica. For the higher-forcing scenarios, the marine basins of the East Antarctic Ice Sheet also become ice free, resulting in a sea-level rise of up to 27 m. The global mean sea-level change after 10 000 years ranges from 9.2 to more than 37 m. For the highest-forcing scenario, the model uncertainty does not exclude the complete melting of the Antarctic ice sheet during the next 10 000 years.


2014 ◽  
Vol 8 (1) ◽  
pp. 85-120 ◽  
Author(s):  
M. N. A. Maris ◽  
S. R. M. Ligtenberg ◽  
M. Crucifix ◽  
B. de Boer ◽  
J. Oerlemans

Abstract. We present the effects of changing two sliding parameters, a deformational velocity parameter and two bedrock deflection parameters on the evolution of the Antarctic Ice Sheet over the period from the last interglacial until the present. These sensitivity experiments have been conducted by running the ice-dynamical model ANICE forward in time. The climatological forcing over time is established by interpolating between two climate states from a regional climate model over time. The interpolation is done in such a way that both temperature and surface mass balance follow the Epica Dome C ice-core proxy record for temperature. We have determined an optimal set of parameter values, for which a realistic grounding line retreat history and present-day ice sheet can be simulated, the simulation with this set of parameter values is defined as the reference simulation. An increase of sliding with respect to this reference simulation leads to a decrease of the Antarctic ice volume due to enhanced ice velocities on mainly the West Antarctic Ice Sheet. The effect of changing the deformational velocity parameter mainly yields a change in East-Antarctic ice volume. Furthermore, we have found a minimum in the Antarctic ice volume during the mid-Holocene. This is a robust feature in our model results, where the strength and the timing of this minimum are both dependent on the investigated parameters. More sliding and a slower responding bedrock lead to a stronger minimum which emerges at an earlier time. From the model results we conclude that the Antarctic Ice Sheet has contributed 10.7 ± 1.3 m of eustatic sea level to the global ocean from the Last Glacial Maximum (about 16 kyr ago for the Antarctic Ice Sheet) until the present.


2021 ◽  
Author(s):  
Clara Burgard ◽  
Nicolas Jourdain

<p>Ocean-induced melting at the base of ice shelves is one of the main drivers of the currently observed mass loss of the Antarctic Ice Sheet. A good understanding of the interaction between ice and ocean at the base of the ice shelves is therefore crucial to understand and project the Antarctic contribution to global sea-level rise. </p><p>Due to the high difficulty to monitor these regions, our understanding of the processes at work beneath ice shelves is limited. Still, several parameterisations of varying complexity have been developed in past decades to describe the ocean-induced sub-shelf melting. These parameterisations can be implemented into standalone ice-sheet models, for example when conducting long-term projections forced with climate model output.</p><p>An assessment of the performance of these parameterisations was conducted in an idealised setup (Favier et al, 2019). However, the application of the better-performing parameterisations in a more realistic setup (e.g. Jourdain et al., 2020) has shown that individual adjustments and corrections are needed for each ice shelf.</p><p>In this study, we revisit the assessment of the parameterisations, this time in a more realistic setup than previous studies. To do so, we apply the different parameterisations on several ice shelves around Antarctica and compare the resulting melt rates to satellite and oceanographic estimates. Based on this comparison, we will refine the parameters and propose an approach to reduce uncertainties in long-term sub-shelf melting projections.</p><p><em>References</em><br><em>- Favier, L., Jourdain, N. C., Jenkins, A., Merino, N., Durand, G., Gagliardini, O., Gillet-Chaulet, F., and Mathiot, P.: Assessment of sub-shelf melting parameterisations using the ocean–ice-sheet coupled model NEMO(v3.6)–Elmer/Ice(v8.3) , Geosci. Model Dev., 12, 2255–2283, https://doi.org/10.5194/gmd-12-2255-2019, 2019. </em><br><em>- Jourdain, N. C., Asay-Davis, X., Hattermann, T., Straneo, F., Seroussi, H., Little, C. M., and Nowicki, S.: A protocol for calculating basal melt rates in the ISMIP6 Antarctic ice sheet projections, The Cryosphere, 14, 3111–3134, https://doi.org/10.5194/tc-14-3111-2020, 2020. </em></p>


2014 ◽  
Vol 8 (4) ◽  
pp. 1347-1360 ◽  
Author(s):  
M. N. A. Maris ◽  
B. de Boer ◽  
S. R. M. Ligtenberg ◽  
M. Crucifix ◽  
W. J. van de Berg ◽  
...  

Abstract. We present the effects of changing two sliding parameters, a deformational velocity parameter and two bedrock deflection parameters on the evolution of the Antarctic ice sheet over the period from the last interglacial until the present. These sensitivity experiments have been conducted by running the dynamic ice model ANICE forward in time. The temporal climatological forcing is established by interpolating between two temporal climate states created with a regional climate model. The interpolation is done in such a way that both temperature and surface mass balance follow the European Project for Ice Coring in Antarctica (EPICA) Dome C ice-core proxy record for temperature. We have determined an optimal set of parameter values, for which a realistic grounding-line retreat history and present-day ice sheet can be simulated; the simulation with this set of parameter values is defined as the reference simulation. An increase of sliding with respect to this reference simulation leads to a decrease of the Antarctic ice volume due to enhanced ice velocities on mainly the West Antarctic ice sheet. The effect of changing the deformational velocity parameter mainly yields a change in east Antarctic ice volume. Furthermore, we have found a minimum in the Antarctic ice volume during the mid-Holocene, in accordance with observations. This is a robust feature in our model results, where the strength and the timing of this minimum are both dependent on the investigated parameters. More sliding and a slower responding bedrock lead to a stronger minimum which emerges at an earlier time. From the model results, we conclude that the Antarctic ice sheet has contributed 10.7 ± 1.3 m of eustatic sea level to the global ocean from the last glacial maximum (about 16 ka for the Antarctic ice sheet) until the present.


Sign in / Sign up

Export Citation Format

Share Document