@article{KreuzerReeseHuiskampetal.2021, author = {Kreuzer, Moritz and Reese, Ronja and Huiskamp, Willem Nicholas and Petri, Stefan and Albrecht, Torsten and Feulner, Georg and Winkelmann, Ricarda}, title = {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}, series = {Geoscientific model development : an interactive open access journal of the European Geosciences Union}, volume = {14}, journal = {Geoscientific model development : an interactive open access journal of the European Geosciences Union}, number = {6}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1991-959X}, doi = {10.5194/gmd-14-3697-2021}, pages = {3697 -- 3714}, year = {2021}, abstract = {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.}, language = {en} } @article{SeroussiNowickiPayneetal.2020, author = {Seroussi, Helene and Nowicki, Sophie and Payne, Antony J. and Goelzer, Heiko and Lipscomb, William H. and Abe-Ouchi, Ayako and Agosta, Cecile and Albrecht, Torsten and Asay-Davis, Xylar and Barthel, Alice and Calov, Reinhard and Cullather, Richard and Dumas, Christophe and Galton-Fenzi, Benjamin K. and Gladstone, Rupert and Golledge, Nicholas R. and Gregory, Jonathan M. and Greve, Ralf and Hattermann, Tore and Hoffman, Matthew J. and Humbert, Angelika and Huybrechts, Philippe and Jourdain, Nicolas C. and Kleiner, Thomas and Larour, Eric and Leguy, Gunter R. and Lowry, Daniel P. and Little, Chistopher M. and Morlighem, Mathieu and Pattyn, Frank and Pelle, Tyler and Price, Stephen F. and Quiquet, Aurelien and Reese, Ronja and Schlegel, Nicole-Jeanne and Shepherd, Andrew and Simon, Erika and Smith, Robin S. and Straneo, Fiammetta and Sun, Sainan and Trusel, Luke D. and Van Breedam, Jonas and van de Wal, Roderik S. W. and Winkelmann, Ricarda and Zhao, Chen and Zhang, Tong and Zwinger, Thomas}, title = {ISMIP6 Antarctica}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {14}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {9}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-14-3033-2020}, pages = {3033 -- 3070}, year = {2020}, abstract = {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.}, language = {en} } @article{FeldmannReeseWinkelmannetal.2022, author = {Feldmann, Johannes and Reese, Ronja and Winkelmann, Ricarda and Levermann, Anders}, title = {Shear-margin melting causes stronger transient ice discharge than ice-stream melting in idealized simulations}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {16}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {5}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-16-1927-2022}, pages = {1927 -- 1940}, year = {2022}, abstract = {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.}, language = {en} } @article{ZeitzReeseBeckmannetal.2021, author = {Zeitz, Maria and Reese, Ronja and Beckmann, Johanna and Krebs-Kanzow, Uta and Winkelmann, Ricarda}, title = {Impact of the melt-albedo feedback on the future evolution of the Greenland Ice Sheet with PISM-dEBM-simple}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {15}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {12}, publisher = {Copernicus}, address = {Katlenburg-Lindau}, issn = {1994-0416}, doi = {10.5194/tc-15-5739-2021}, pages = {5739 -- 5764}, year = {2021}, abstract = {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.}, language = {en} } @article{LevermannWinkelmann2016, author = {Levermann, Anders and Winkelmann, Ricarda}, title = {A simple equation for the melt elevation feedback of ice sheets}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {10}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-10-1799-2016}, pages = {1799 -- 1807}, year = {2016}, abstract = {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.}, language = {en} } @article{SeroussiNowickiSimonetal.2019, author = {Seroussi, Helene and Nowicki, Sophie and Simon, Erika and Abe-Ouchi, Ayako and Albrecht, Torsten and Brondex, Julien and Cornford, Stephen and Dumas, Christophe and Gillet-Chaulet, Fabien and Goelzer, Heiko and Golledge, Nicholas R. and Gregory, Jonathan M. and Greve, Ralf and Hoffman, Matthew J. and Humbert, Angelika and Huybrechts, Philippe and Kleiner, Thomas and Larourl, Eric and Leguy, Gunter and Lipscomb, William H. and Lowry, Daniel and Mengel, Matthias and Morlighem, Mathieu and Pattyn, Frank and Payne, Anthony J. and Pollard, David and Price, Stephen F. and Quiquet, Aurelien and Reerink, Thomas J. and Reese, Ronja and Rodehacke, Christian B. and Schlegel, Nicole-Jeanne and Shepherd, Andrew and Sun, Sainan and Sutter, Johannes and Van Breedam, Jonas and van de Wal, Roderik S. W. and Winkelmann, Ricarda and Zhang, Tong}, title = {initMIP-Antarctica}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {13}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {5}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-13-1441-2019}, pages = {1441 -- 1471}, year = {2019}, abstract = {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.}, language = {en} } @article{CiemerBoersHirotaetal.2019, author = {Ciemer, Catrin and Boers, Niklas and Hirota, Marina and Kurths, J{\"u}rgen and M{\"u}ller-Hansen, Finn and Oliveira, Rafael S. and Winkelmann, Ricarda}, title = {Higher resilience to climatic disturbances in tropical vegetation exposed to more variable rainfall}, series = {Nature geoscience}, volume = {12}, journal = {Nature geoscience}, number = {3}, publisher = {Nature Publ. Group}, address = {New York}, issn = {1752-0894}, doi = {10.1038/s41561-019-0312-z}, pages = {174 -- 179}, year = {2019}, abstract = {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.}, language = {en} } @misc{CiemerRehmKurthsetal.2020, author = {Ciemer, Catrin and Rehm, Lars and Kurths, J{\"u}rgen and Donner, Reik Volker and Winkelmann, Ricarda and Boers, Niklas}, title = {An early-warning indicator for Amazon droughts exclusively based on tropical Atlantic sea surface temperatures}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {9}, issn = {1866-8372}, doi = {10.25932/publishup-52586}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-525863}, pages = {12}, year = {2020}, abstract = {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.}, language = {en} } @article{ReeseWinkelmannGudmundsson2018, author = {Reese, Ronja and Winkelmann, Ricarda and Gudmundsson, Gudmundur Hilmar}, title = {Grounding-line flux formula applied as a flux condition in numerical simulations fails for buttressed Antarctic ice streams}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {12}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {10}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-12-3229-2018}, pages = {3229 -- 3242}, year = {2018}, abstract = {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.}, language = {en} } @article{CiemerRehmKurthsetal.2020, author = {Ciemer, Catrin and Rehm, Lars and Kurths, J{\"u}rgen and Donner, Reik Volker and Winkelmann, Ricarda and Boers, Niklas}, title = {An early-warning indicator for Amazon droughts exclusively based on tropical Atlantic sea surface temperatures}, series = {Environmental Research Letters}, volume = {15}, journal = {Environmental Research Letters}, number = {9}, publisher = {IOP - Institute of Physics Publishing}, address = {Bristol}, pages = {10}, year = {2020}, abstract = {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.}, language = {en} }