TY - JOUR A1 - Kreuzer, Moritz A1 - Reese, Ronja A1 - Huiskamp, Willem Nicholas A1 - Petri, Stefan A1 - Albrecht, Torsten A1 - Feulner, Georg A1 - Winkelmann, Ricarda T1 - Coupling framework (1.0) for the PISM (1.1.4) ice sheet model and the MOMS (5.1.0) ocean model via the PICO ice shelf cavity model in an Antarctic domain JF - Geoscientific model development : an interactive open access journal of the European Geosciences Union N2 - 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. Y1 - 2021 U6 - https://doi.org/10.5194/gmd-14-3697-2021 SN - 1991-959X SN - 1991-9603 VL - 14 IS - 6 SP - 3697 EP - 3714 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Seroussi, Helene A1 - Nowicki, Sophie A1 - Payne, Antony J. A1 - Goelzer, Heiko A1 - Lipscomb, William H. A1 - Abe-Ouchi, Ayako A1 - Agosta, Cecile A1 - Albrecht, Torsten A1 - Asay-Davis, Xylar A1 - Barthel, Alice A1 - Calov, Reinhard A1 - Cullather, Richard A1 - Dumas, Christophe A1 - Galton-Fenzi, Benjamin K. A1 - Gladstone, Rupert A1 - Golledge, Nicholas R. A1 - Gregory, Jonathan M. A1 - Greve, Ralf A1 - Hattermann, Tore A1 - Hoffman, Matthew J. A1 - Humbert, Angelika A1 - Huybrechts, Philippe A1 - Jourdain, Nicolas C. A1 - Kleiner, Thomas A1 - Larour, Eric A1 - Leguy, Gunter R. A1 - Lowry, Daniel P. A1 - Little, Chistopher M. A1 - Morlighem, Mathieu A1 - Pattyn, Frank A1 - Pelle, Tyler A1 - Price, Stephen F. A1 - Quiquet, Aurelien A1 - Reese, Ronja A1 - Schlegel, Nicole-Jeanne A1 - Shepherd, Andrew A1 - Simon, Erika A1 - Smith, Robin S. A1 - Straneo, Fiammetta A1 - Sun, Sainan A1 - Trusel, Luke D. A1 - Van Breedam, Jonas A1 - van de Wal, Roderik S. W. A1 - Winkelmann, Ricarda A1 - Zhao, Chen A1 - Zhang, Tong A1 - Zwinger, Thomas T1 - ISMIP6 Antarctica BT - a multi-model ensemble of the Antarctic ice sheet evolution over the 21st century JF - The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union N2 - 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. Y1 - 2020 U6 - https://doi.org/10.5194/tc-14-3033-2020 SN - 1994-0416 SN - 1994-0424 VL - 14 IS - 9 SP - 3033 EP - 3070 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Feldmann, Johannes A1 - Reese, Ronja A1 - Winkelmann, Ricarda A1 - Levermann, Anders T1 - Shear-margin melting causes stronger transient ice discharge than ice-stream melting in idealized simulations JF - The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union N2 - Basal ice-shelf melting is the key driver of Antarctica's increasing sea-level contribution. In diminishing the buttressing force of the ice shelves that fringe the ice sheet, the melting increases the ice discharge into the ocean. Here we contrast the influence of basal melting in two different ice-shelf regions on the time-dependent response of an isothermal, inherently buttressed ice-sheet-shelf system. In the idealized numerical simulations, the basal-melt perturbations are applied close to the grounding line in the ice-shelf's (1) ice-stream region, where the ice shelf is fed by the fastest ice masses that stream through the upstream bed trough and (2) shear margins, where the ice flow is slower. The results show that melting below one or both of the shear margins can cause a decadal to centennial increase in ice discharge that is more than twice as large compared to a similar perturbation in the ice-stream region. We attribute this to the fact that melt-induced ice-shelf thinning in the central grounding-line region is attenuated very effectively by the fast flow of the central ice stream. In contrast, the much slower ice dynamics in the lateral shear margins of the ice shelf facilitate sustained ice-shelf thinning and thereby foster buttressing reduction. Regardless of the melt location, a higher melt concentration toward the grounding line generally goes along with a stronger response. Our results highlight the vulnerability of outlet glaciers to basal melting in stagnant, buttressing-relevant ice-shelf regions, a mechanism that may gain importance under future global warming. Y1 - 2022 U6 - https://doi.org/10.5194/tc-16-1927-2022 SN - 1994-0416 SN - 1994-0424 VL - 16 IS - 5 SP - 1927 EP - 1940 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Zeitz, Maria A1 - Reese, Ronja A1 - Beckmann, Johanna A1 - Krebs-Kanzow, Uta A1 - Winkelmann, Ricarda T1 - Impact of the melt-albedo feedback on the future evolution of the Greenland Ice Sheet with PISM-dEBM-simple JF - The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union N2 - 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. Y1 - 2021 U6 - https://doi.org/10.5194/tc-15-5739-2021 SN - 1994-0416 SN - 1994-0424 VL - 15 IS - 12 SP - 5739 EP - 5764 PB - Copernicus CY - Katlenburg-Lindau ER - TY - JOUR A1 - Seroussi, Helene A1 - Nowicki, Sophie A1 - Simon, Erika A1 - Abe-Ouchi, Ayako A1 - Albrecht, Torsten A1 - Brondex, Julien A1 - Cornford, Stephen A1 - Dumas, Christophe A1 - Gillet-Chaulet, Fabien A1 - Goelzer, Heiko A1 - Golledge, Nicholas R. A1 - Gregory, Jonathan M. A1 - Greve, Ralf A1 - Hoffman, Matthew J. A1 - Humbert, Angelika A1 - Huybrechts, Philippe A1 - Kleiner, Thomas A1 - Larourl, Eric A1 - Leguy, Gunter A1 - Lipscomb, William H. A1 - Lowry, Daniel A1 - Mengel, Matthias A1 - Morlighem, Mathieu A1 - Pattyn, Frank A1 - Payne, Anthony J. A1 - Pollard, David A1 - Price, Stephen F. A1 - Quiquet, Aurelien A1 - Reerink, Thomas J. A1 - Reese, Ronja A1 - Rodehacke, Christian B. A1 - Schlegel, Nicole-Jeanne A1 - Shepherd, Andrew A1 - Sun, Sainan A1 - Sutter, Johannes A1 - Van Breedam, Jonas A1 - van de Wal, Roderik S. W. A1 - Winkelmann, Ricarda A1 - Zhang, Tong T1 - initMIP-Antarctica BT - an ice sheet model initialization experiment of ISMIP6 JF - The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union N2 - 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. Y1 - 2019 U6 - https://doi.org/10.5194/tc-13-1441-2019 SN - 1994-0416 SN - 1994-0424 VL - 13 IS - 5 SP - 1441 EP - 1471 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Reese, Ronja A1 - Winkelmann, Ricarda A1 - Gudmundsson, Gudmundur Hilmar T1 - Grounding-line flux formula applied as a flux condition in numerical simulations fails for buttressed Antarctic ice streams JF - The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union N2 - 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. Y1 - 2018 U6 - https://doi.org/10.5194/tc-12-3229-2018 SN - 1994-0416 SN - 1994-0424 VL - 12 IS - 10 SP - 3229 EP - 3242 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Reese, Ronja A1 - Albrecht, Torsten A1 - Mengel, Matthias A1 - Asay-Davis, Xylar A1 - Winkelmann, Ricarda T1 - Antarctic sub-shelf melt rates via PICO JF - The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union N2 - 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. Y1 - 2018 U6 - https://doi.org/10.5194/tc-12-1969-2018 SN - 1994-0416 SN - 1994-0424 VL - 12 IS - 6 SP - 1969 EP - 1985 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Reese, Ronja A1 - Levermann, Anders A1 - Albrecht, Torsten A1 - Seroussi, Helene A1 - Winkelmann, Ricarda T1 - The role of history and strength of the oceanic forcing in sea level projections from Antarctica with the Parallel Ice Sheet Model JF - The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union N2 - Mass loss from the Antarctic Ice Sheet constitutes the largest uncertainty in projections of future sea level rise. Ocean-driven melting underneath the floating ice shelves and subsequent acceleration of the inland ice streams are the major reasons for currently observed mass loss from Antarctica and are expected to become more important in the future. Here we show that for projections of future mass loss from the Antarctic Ice Sheet, it is essential (1) to better constrain the sensitivity of sub-shelf melt rates to ocean warming and (2) to include the historic trajectory of the ice sheet. In particular, we find that while the ice sheet response in simulations using the Parallel Ice Sheet Model is comparable to the median response of models in three Antarctic Ice Sheet Intercomparison projects - initMIP, LARMIP-2 and ISMIP6 - conducted with a range of ice sheet models, the projected 21st century sea level contribution differs significantly depending on these two factors. For the highest emission scenario RCP8.5, this leads to projected ice loss ranging from 1:4 to 4:0 cm of sea level equivalent in simulations in which ISMIP6 ocean forcing drives the PICO ocean box model where parameter tuning leads to a comparably low sub-shelf melt sensitivity and in which no surface forcing is applied. This is opposed to a likely range of 9:1 to 35:8 cm using the exact same initial setup, but emulated from the LARMIP-2 experiments with a higher melt sensitivity, even though both projects use forcing from climate models and melt rates are calibrated with previous oceanographic studies. Furthermore, using two initial states, one with a previous historic simulation from 1850 to 2014 and one starting from a steady state, we show that while differences between the ice sheet configurations in 2015 seem marginal at first sight, the historic simulation increases the susceptibility of the ice sheet to ocean warming, thereby increasing mass loss from 2015 to 2100 by 5% to 50 %. Hindcasting past ice sheet changes with numerical models would thus provide valuable tools to better constrain projections. Our results emphasize that the uncertainty that arises from the forcing is of the same order of magnitude as the ice dynamic response for future sea level projections. Y1 - 2020 U6 - https://doi.org/10.5194/tc-14-3097-2020 SN - 1994-0416 SN - 1994-0424 VL - 14 IS - 9 SP - 3097 EP - 3110 PB - Copernicus CY - Göttingen ER -