TY - JOUR A1 - Winkelmann, Ricarda A1 - Martin, Maria A. A1 - Haseloff, Monika A1 - Albrecht, Torsten A1 - Bueler, Ed A1 - Khroulev, C. A1 - Levermann, Anders T1 - The Potsdam parallel ice sheet model (PISM-PIK) - Part 1: Model description JF - The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union N2 - We present the Potsdam Parallel Ice Sheet Model (PISM-PIK), developed at the Potsdam Institute for Climate Impact Research to be used for simulations of large-scale ice sheet-shelf systems. It is derived from the Parallel Ice Sheet Model (Bueler and Brown, 2009). Velocities are calculated by superposition of two shallow stress balance approximations within the entire ice covered region: the shallow ice approximation (SIA) is dominant in grounded regions and accounts for shear deformation parallel to the geoid. The plug-flow type shallow shelf approximation (SSA) dominates the velocity field in ice shelf regions and serves as a basal sliding velocity in grounded regions. Ice streams can be identified diagnostically as regions with a significant contribution of membrane stresses to the local momentum balance. All lateral boundaries in PISM-PIK are free to evolve, including the grounding line and ice fronts. Ice shelf margins in particular are modeled using Neumann boundary conditions for the SSA equations, reflecting a hydrostatic stress imbalance along the vertical calving face. The ice front position is modeled using a subgrid-scale representation of calving front motion (Albrecht et al., 2011) and a physically-motivated calving law based on horizontal spreading rates. The model is tested in experiments from the Marine Ice Sheet Model Intercomparison Project (MISMIP). A dynamic equilibrium simulation of Antarctica under present-day conditions is presented in Martin et al. (2011). Y1 - 2011 U6 - https://doi.org/10.5194/tc-5-715-2011 SN - 1994-0416 VL - 5 IS - 3 SP - 715 EP - 726 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Martin, Maria A. A1 - Winkelmann, Ricarda A1 - Haseloff, M. A1 - Albrecht, Tanja A1 - Bueler, Ed A1 - Khroulev, C. A1 - Levermann, Anders T1 - The Potsdam parallel ice sheet model (PISM-PIK) - Part 2: Dynamic equilibrium simulation of the Antarctic ice sheet JF - The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union N2 - We present a dynamic equilibrium simulation of the ice sheet-shelf system on Antarctica with the Potsdam Parallel Ice Sheet Model (PISM-PIK). The simulation is initialized with present-day conditions for bed topography and ice thickness and then run to steady state with constant present-day surface mass balance. Surface temperature and sub-shelf basal melt distribution are parameterized. Grounding lines and calving fronts are free to evolve, and their modeled equilibrium state is compared to observational data. A physically-motivated calving law based on horizontal spreading rates allows for realistic calving fronts for various types of shelves. Steady-state dynamics including surface velocity and ice flux are analyzed for whole Antarctica and the Ronne-Filchner and Ross ice shelf areas in particular. The results show that the different flow regimes in sheet and shelves, and the transition zone between them, are captured reasonably well, supporting the approach of superposition of SIA and SSA for the representation of fast motion of grounded ice. This approach also leads to a natural emergence of sliding-dominated flow in stream-like features in this new 3-D marine ice sheet model. Y1 - 2011 U6 - https://doi.org/10.5194/tc-5-727-2011 SN - 1994-0416 VL - 5 IS - 3 SP - 727 EP - 740 PB - Copernicus CY - Göttingen ER - TY - THES A1 - Reese, Ronja T1 - The far reach of ice-shelf thinning in Antarctica Y1 - 2018 ER - TY - THES A1 - Winkelmann, Ricarda T1 - The future sea-level contribution from antartica: projections of solid ice discharge Y1 - 2012 CY - Potsdam ER - TY - THES A1 - Wunderling, Nico T1 - Nichtlineare Dynamiken und Interaktionen von Kippelementen im Erdsystem T1 - Nonlinear dynamics and interactions of tipping elements in the Earth system N2 - With ongoing anthropogenic global warming, some of the most vulnerable components of the Earth system might become unstable and undergo a critical transition. These subsystems are the so-called tipping elements. They are believed to exhibit threshold behaviour and would, if triggered, result in severe consequences for the biosphere and human societies. Furthermore, it has been shown that climate tipping elements are not isolated entities, but interact across the entire Earth system. Therefore, this thesis aims at mapping out the potential for tipping events and feedbacks in the Earth system mainly by the use of complex dynamical systems and network science approaches, but partially also by more detailed process-based models of the Earth system. In the first part of this thesis, the theoretical foundations are laid by the investigation of networks of interacting tipping elements. For this purpose, the conditions for the emergence of global cascades are analysed against the structure of paradigmatic network types such as Erdös-Rényi, Barabási-Albert, Watts-Strogatz and explicitly spatially embedded networks. Furthermore, micro-scale structures are detected that are decisive for the transition of local to global cascades. These so-called motifs link the micro- to the macro-scale in the network of tipping elements. Alongside a model description paper, all these results are entered into the Python software package PyCascades, which is publicly available on github. In the second part of this dissertation, the tipping element framework is first applied to components of the Earth system such as the cryosphere and to parts of the biosphere. Afterwards it is applied to a set of interacting climate tipping elements on a global scale. Using the Earth system Model of Intermediate Complexity (EMIC) CLIMBER-2, the temperature feedbacks are quantified, which would arise if some of the large cryosphere elements disintegrate over a long span of time. The cryosphere components that are investigated are the Arctic summer sea ice, the mountain glaciers, the Greenland and the West Antarctic Ice Sheets. The committed temperature increase, in case the ice masses disintegrate, is on the order of an additional half a degree on a global average (0.39-0.46 °C), while local to regional additional temperature increases can exceed 5 °C. This means that, once tipping has begun, additional reinforcing feedbacks are able to increase global warming and with that the risk of further tipping events. This is also the case in the Amazon rainforest, whose parts are dependent on each other via the so-called moisture-recycling feedback. In this thesis, the importance of drought-induced tipping events in the Amazon rainforest is investigated in detail. Despite the Amazon rainforest is assumed to be adapted to past environmental conditions, it is found that tipping events sharply increase if the drought conditions become too intense in a too short amount of time, outpacing the adaptive capacity of the Amazon rainforest. In these cases, the frequency of tipping cascades also increases to 50% (or above) of all tipping events. In the model that was developed in this study, the southeastern region of the Amazon basin is hit hardest by the simulated drought patterns. This is also the region that already nowadays suffers a lot from extensive human-induced changes due to large-scale deforestation, cattle ranching or infrastructure projects. Moreover, on the larger Earth system wide scale, a network of conceptualised climate tipping elements is constructed in this dissertation making use of a large literature review, expert knowledge and topological properties of the tipping elements. In global warming scenarios, tipping cascades are detected even under modest scenarios of climate change, limiting global warming to 2 °C above pre-industrial levels. In addition, the structural roles of the climate tipping elements in the network are revealed. While the large ice sheets on Greenland and Antarctica are the initiators of tipping cascades, the Atlantic Meridional Overturning Circulation (AMOC) acts as the transmitter of cascades. Furthermore, in our conceptual climate tipping element model, it is found that the ice sheets are of particular importance for the stability of the entire system of investigated climate tipping elements. In the last part of this thesis, the results from the temperature feedback study with the EMIC CLIMBER-2 are combined with the conceptual model of climate tipping elements. There, it is observed that the likelihood of further tipping events slightly increases due to the temperature feedbacks even if no further CO$_2$ would be added to the atmosphere. Although the developed network model is of conceptual nature, it is possible with this work for the first time to quantify the risk of tipping events between interacting components of the Earth system under global warming scenarios, by allowing for dynamic temperature feedbacks at the same time. N2 - Bei fortdauerndem anthropogenem Klimawandel, könnten einige der vulnerabelsten Komponenten des Erdsystem instabil werden und in einen anderen Zustand übergehen. Diese Komponenten des Erdsystems sind die sogenannten Kippelemente. Bei ihnen wird angenommen, dass sie einen Kipppunkt besitzen ab dem sie in einen qualitativ anderen Zustand übergehen können. Sollte das passieren, hätte das schwerwiegende Konsequenzen für die Biosphäre und menschliche Gesellschaften. Des Weiteren ist gezeigt worden, dass Kippelemente keine isolierte Reigionen oder Prozesse sind, sondern über das gesamte Erdsystem hinweg interagieren. Das Ziel dieser Arbeit ist es daher, die Wahrscheinlichkeit für Kippereignisse sowie deren Feedbacks im Erdsystem zu quantifizieren. Zu diesem Zweck kommen vor allem Frameworks aus der Wissenschaft komplexer Systeme und Netzwerke zum Einsatz. Für einige Teilaspekte dieser Arbeit wird aber auch ein detaillierteres und prozessbasierteres Erdsystemmodell verwendet. Im ersten Teil dieser Arbeit werden die theoretischen Grundlagen gelegt, indem komplexe Netzwerke bestehend aus interagierenden Kippelementen untersucht werden. Hier werden Voraussetzungen für das Auftreten globaler Kippkaskaden anhand der Struktur paradigmatischer Netzwerktypen analysiert. Diese Typen sind Netzwerke wie Erdös-Rényi, Barabási-Albert, Watts-Strogatz Netzwerke oder auch explizit räumlich eingebettete Netzwerke. Darüber hinaus sind bestimmte Mikrostrukturen in Netzwerken dafür entscheidend, ob sich eine lokale Kaskaden auf das globale Netzwerk ausbreiten kann. Diese Strukturen sind das Bindeglied zwischen der Mikro- und der Makroebene des Netzwerks und werden Motive genannt. Zusammen mit einer Publikation zur Modellbeschreibung, werden alle diese Ergebnisse im Python-Softwarepaket PyCascades veröffentlicht, das auf github öffentlich verfügbar ist. Im zweiten Teil dieser Dissertation wird das Kippelementframework zunächst auf Kompenenten des Erdsystems angewendet wie der Kryosphäre und Teilen der Biosphäre, und danach auf globaler Skala für interagierende Klimakippelemente. In einem ersten Schritt werden mit dem Erdsystemmodell mittlerer Komplexität CLIMBER-2 die Temperaturfeedbacks ermittelt, die entstehen würden, wenn große Gebiete der Kryosphäre auf lange Sicht eisfrei werden. In dieser Berechnung werden das arktische Sommermeereis, die Gebirgsgletscher, der grönländische und der westantarktische Eisschild berücksichtigt. Die quantifizierte Temperaturerhöhung liegt in der Größenordnung von einem halben Grad zusätzlicher globaler Erwärmung (0.39--0.46°C). Lokale bis regionale Temperaturerhöhungen können allerdings 5°C übersteigen. Wenn also das Kippen einiger Elemente begonnen hat, bedeutet dieses Ergebnis, dass Temperaturfeedbacks in der Lage sind, das Risiko weiterer Kippereignisse zu erhöhen. Dies ist auch der Fall im Amazonasregenwald, dessen Unterregionen über den sogenannten Feuchtig-keits-Recycling-Feedback miteinander in Beziehung stehen und voneinander abhängen. In dieser Dissertation wird die Bedeutung von Kippereignissen im Detail untersucht, die aufgrund von Dürreperioden zustande kommen. Obwohl man davon ausgehen kann, dass der Regenwald sich an zurückliegende und gegenwärtige Klimabedingungen angepasst hat, kann festgestellt werden, dass die Häu-figkeit von Kippereignissen stark zunimmt, wenn die jeweilige Trockenperiode eine gewisse Intensität übersteigt und damit die Anpassungsfähigkeit des Amazonasregenwalds überschritten wird. In solchen Fällen steigt auch die Häufigkeit von Kippkaskaden unter allen Kippereignissen auf 50% (und mehr) an. In dem Modell, das in dieser Studie entwickelt wurde, zeigt sich, dass der Südosten des Amazonasbeckens am stärksten von den simulierten Trockenheitsmustern betroffen ist. Das ist gleichzeitig die Region, die bereits heute stark unter anthropogener Veränderung leidet, unter anderem aufgrund von großflächiger Abholzung, Viehzucht oder Infrastrukturprojekten. Zudem wird in dieser Dissertation auf der größeren, erdsystemweiten Skala ein Netzwerk konzeptionalisierter Klimakippelemente aufgebaut. Zu diesem Zweck wird eine umfangreiche Literaturrecherche durchgeführt, die zusammen mit Expertenwissen und den topologischen Eigenschaften der Kippelemente in die Studien mit einfließt. In Klimawandelszenarien können dann Kippkaskaden beobachtet werden, selbst wenn die globale Erderwärmung auf 2°C über dem vorindustriellen Niveau begrenzt werden kann. Außerdem werden die strukturellen Rollen der Klimakippelemente im Netzwerk ermittelt. Während die großen Eisschilde auf Grönland und der Westantarktis viele Kippkaskaden initiieren, ist die Atlantische Umwälzzirkulation für die Weitergabe vieler dieser Kaskaden verantwortlich. In unserem konzeptionellen Modell für Klimakippelemente wird darüber hinaus festgestellt, dass die Eisschilde von besonderer Bedeutung für die Stabilität des Gesamtsystems sind. Im letzen Teil dieser Dissertation werden die Ergebnisse der Feedbackstudie (CLIMBER-2-Studie) zusammengebracht mit dem konzeptionellen Klimakippelementmodell. Dabei zeigt sich, dass die Wahrscheinlichkeit zusätzlicher Kippereignisse aufgrund der berücksichtigten Temperaturfeedbacks auch ohne das Zuführen eines zusätzlichen CO2-Eintrags in die Atmosphäre leicht ansteigt. Trotz der konzeptionellen Natur des entwickelten Netzwerkmodells, ist es mit dieser Arbeit erstmals möglich eine Risikoabschätzung über das Auftreten von Kippkaskaden im Erdsystem vorzunehmen. Darüber hinaus können, unter der Annahme globaler Erwärmungsszenarien, auch dynamische Temperaturfeedbacks berücksichtigt werden. KW - tipping element KW - nonlinear dynamics KW - tipping cascade KW - climate change KW - complex networks KW - Klimawandel KW - komplexe Netzwerke KW - nichtlineare Dynamiken KW - Kippkaskade KW - Kippelement Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-525140 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 - 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 - TY - JOUR A1 - Zeitz, Maria A1 - Levermann, Anders A1 - Winkelmann, Ricarda T1 - Sensitivity of ice loss to uncertainty in flow law parameters in an idealized one-dimensional geometry JF - The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union N2 - Acceleration of the flow of ice drives mass losses in both the Antarctic and the Greenland Ice Sheet. The projections of possible future sea-level rise rely on numerical ice-sheet models, which solve the physics of ice flow, melt, and calving. While major advancements have been made by the ice-sheet modeling community in addressing several of the related uncertainties, the flow law, which is at the center of most process-based ice-sheet models, is not in the focus of the current scientific debate. However, recent studies show that the flow law parameters are highly uncertain and might be different from the widely accepted standard values. Here, we use an idealized flow-line setup to investigate how these uncertainties in the flow law translate into uncertainties in flow-driven mass loss. In order to disentangle the effect of future warming on the ice flow from other effects, we perform a suite of experiments with the Parallel Ice Sheet Model (PISM), deliberately excluding changes in the surface mass balance. We find that changes in the flow parameters within the observed range can lead up to a doubling of the flow-driven mass loss within the first centuries of warming, compared to standard parameters. The spread of ice loss due to the uncertainty in flow parameters is on the same order of magnitude as the increase in mass loss due to surface warming. While this study focuses on an idealized flow-line geometry, it is likely that this uncertainty carries over to realistic three-dimensional simulations of Greenland and Antarctica. Y1 - 2020 U6 - https://doi.org/10.5194/tc-14-3537-2020 SN - 1994-0416 SN - 1994-0424 VL - 14 IS - 10 SP - 3537 EP - 3550 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Reese, Ronja A1 - Gudmundsson, Gudmundur Hilmar A1 - Levermann, Anders A1 - Winkelmann, Ricarda T1 - The far reach of ice-shelf thinning in Antarctica JF - Nature climate change N2 - Floating ice shelves, which fringe most of Antarctica’s coastline, regulate ice flow into the Southern Ocean1,2,3. Their thinning4,5,6,7 or disintegration8,9 can cause upstream acceleration of grounded ice and raise global sea levels. So far the effect has not been quantified in a comprehensive and spatially explicit manner. Here, using a finite-element model, we diagnose the immediate, continent-wide flux response to different spatial patterns of ice-shelf mass loss. We show that highly localized ice-shelf thinning can reach across the entire shelf and accelerate ice flow in regions far from the initial perturbation. As an example, this ‘tele-buttressing’ enhances outflow from Bindschadler Ice Stream in response to thinning near Ross Island more than 900 km away. We further find that the integrated flux response across all grounding lines is highly dependent on the location of imposed changes: the strongest response is caused not only near ice streams and ice rises, but also by thinning, for instance, well-within the Filchner–Ronne and Ross Ice Shelves. The most critical regions in all major ice shelves are often located in regions easily accessible to the intrusion of warm ocean waters10,11,12, stressing Antarctica’s vulnerability to changes in its surrounding ocean. Y1 - 2017 U6 - https://doi.org/10.1038/s41558-017-0020-x SN - 1758-678X SN - 1758-6798 VL - 8 IS - 1 SP - 53 EP - 57 PB - Nature Publ. Group CY - London ER - 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 - Zeitz, Maria A1 - Haacker, Jan M. A1 - Donges, Jonathan A1 - Albrecht, Torsten A1 - Winkelmann, Ricarda T1 - Dynamic regimes of the Greenland Ice Sheet emerging from interacting melt-elevation and glacial isostatic adjustment feedbacks JF - Earth system dynamics N2 - 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. Y1 - 2022 U6 - https://doi.org/10.5194/esd-13-1077-2022 SN - 2190-4979 SN - 2190-4987 VL - 13 IS - 3 SP - 1077 EP - 1096 PB - Copernicus Publ. CY - Göttingen ER - TY - JOUR A1 - Albrecht, Torsten A1 - Winkelmann, Ricarda A1 - Levermann, Anders T1 - Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM) BT - part 2: parameter ensemble analysis JF - The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union N2 - 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. Y1 - 2020 U6 - https://doi.org/10.5194/tc-14-633-2020 SN - 1994-0416 SN - 1994-0424 VL - 14 IS - 2 SP - 633 EP - 656 PB - Copernicus Publ. CY - Göttingen ER - TY - JOUR A1 - Garbe, Julius A1 - Albrecht, Torsten A1 - Levermann, Anders A1 - Donges, Jonathan A1 - Winkelmann, Ricarda T1 - The hysteresis of the Antarctic Ice Sheet JF - Nature : the international weekly journal of science N2 - More than half of Earth's freshwater resources are held by the Antarctic Ice Sheet, which thus represents by far the largest potential source for global sea-level rise under future warming conditions(1). Its long-term stability determines the fate of our coastal cities and cultural heritage. Feedbacks between ice, atmosphere, ocean, and the solid Earth give rise to potential nonlinearities in its response to temperature changes. So far, we are lacking a comprehensive stability analysis of the Antarctic Ice Sheet for different amounts of global warming. Here we show that the Antarctic Ice Sheet exhibits a multitude of temperature thresholds beyond which ice loss is irreversible. Consistent with palaeodata(2)we find, using the Parallel Ice Sheet Model(3-5), that at global warming levels around 2 degrees Celsius above pre-industrial levels, West Antarctica is committed to long-term partial collapse owing to the marine ice-sheet instability. Between 6 and 9 degrees of warming above pre-industrial levels, the loss of more than 70 per cent of the present-day ice volume is triggered, mainly caused by the surface elevation feedback. At more than 10 degrees of warming above pre-industrial levels, Antarctica is committed to become virtually ice-free. The ice sheet's temperature sensitivity is 1.3 metres of sea-level equivalent per degree of warming up to 2 degrees above pre-industrial levels, almost doubling to 2.4 metres per degree of warming between 2 and 6 degrees and increasing to about 10 metres per degree of warming between 6 and 9 degrees. Each of these thresholds gives rise to hysteresis behaviour: that is, the currently observed ice-sheet configuration is not regained even if temperatures are reversed to present-day levels. In particular, the West Antarctic Ice Sheet does not regrow to its modern extent until temperatures are at least one degree Celsius lower than pre-industrial levels. Our results show that if the Paris Agreement is not met, Antarctica's long-term sea-level contribution will dramatically increase and exceed that of all other sources.
Modelling shows that the Antarctic Ice Sheet exhibits multiple temperature thresholds beyond which ice loss would become irreversible, and once melted, the ice sheet can regain its previous mass only if the climate cools well below pre-industrial temperatures. Y1 - 2020 U6 - https://doi.org/10.1038/s41586-020-2727-5 SN - 0028-0836 SN - 1476-4687 VL - 585 IS - 7826 SP - 538 EP - 544 PB - Macmillan Publishers Limited CY - Berlin ER - TY - JOUR A1 - Schlemm, Tanja A1 - Feldmann, Johannes A1 - Winkelmann, Ricarda A1 - Levermann, Anders T1 - Stabilizing effect of melange buttressing on the marine ice-cliff instability of the West Antarctic Ice Sheet JF - The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union N2 - Owing to global warming and particularly high regional ocean warming, both Thwaites and Pine Island Glaciers in the Amundsen region of the Antarctic Ice Sheet could lose their buttressing ice shelves over time. We analyse the possible consequences using the parallel ice sheet model (PISM), applying a simple cliff-calving parameterization and an ice melange-buttressing model. We find that the instantaneous loss of ice-shelf buttressing, due to enforced ice-shelf melting, initiates grounding-line retreat and triggers marine ice sheet instability (MISI). As a consequence, the grounding line progresses into the interior of the West Antarctic Ice Sheet and leads to a sea level contribution of 0.6 m within 100 a. By subjecting the exposed ice cliffs to cliff calving using our simplified parameterization, we also analyse marine ice cliff instability (MICI). In our simulations it can double or even triple the sea level contribution depending on the only loosely constrained parameter that determines the maximum cliff-calving rate. The speed of MICI depends on this upper bound of the calving rate, which is given by the ice melange buttressing the glacier. However, stabilization of MICI may occur for geometric reasons. Because the embayment geometry changes as MICI advances into the interior of the ice sheet, the upper bound on calving rates is reduced and the progress of MICI is slowed down. Although we cannot claim that our simulations bear relevant quantitative estimates of the effect of ice-melange buttressing on MICI, the mechanism has the potential to stop the instability. Further research is needed to evaluate its role for the past and future evolution of the Antarctic Ice Sheet. Y1 - 2022 U6 - https://doi.org/10.5194/tc-16-1979-2022 SN - 1994-0416 SN - 1994-0424 VL - 16 IS - 5 SP - 1979 EP - 1996 PB - Copernicus CY - Göttingen ER - TY - THES A1 - Bastian, Martin T1 - An emergent machine learning approach for seasonal cyclone activity forecasts N2 - Seasonal forecasts are of great interest in many areas. Knowing the amount of precipitation for the upcoming season in regions of water scarcity would facilitate a better water management. If farmers knew the weather conditions of the upcoming summer at sowing time, they could select those cereal species that are best adapted to these conditions. This would allow farmers to improve the harvest and potentially even reduce the amount of pesticides used. However, the undoubted advantages of seasonal forecasts are often opposed by their high degree of uncertainty. The great challenge of generating seasonal forecasts with lead times of several months mainly originates from the chaotic nature of the earth system. In a chaotic system, even tiny differences in the initial conditions can lead to strong deviations in the system’s state in the long run. In this dissertation we propose an emergent machine learning approach for seasonal forecasting, called the AnlgModel. The AnlgModel combines the analogue method with myopic feature selection and bootstrapping. To benchmark the abilities of the AnlgModel we apply it to seasonal cyclone activity forecasts in the North Atlantic and Northwest Pacific. The AnlgModel demonstrates competitive hindcast skills with two operational forecasts and even outperforms these for long lead times. In the second chapter we comprehend the forecasting strategy of the Anlg-Model. We thereby analyse the analogue selection process for the 2017 North Atlantic and the 2018 Northwest Pacific seasonal cyclone activity. The analysis shows that those climate indices which are known to influence the seasonal cyclone activity, such as the Niño 3.4 SST, are correctly represented among the selected analogues. Furthermore the selected analogues reflect large-scale climate patterns that were identified by expert reports as being determinative for these particular seasons. In the third chapter we analyse the features that are used by the AnlgModel for its predictions. We therefore inspect the feature relevance (FR). The FR patterns learned by the AnlgModel show a high congruence with the predictor regions used by the operational forecasts. However, the AnlgModel also discovered new features, such as the SST anomaly in the Gulf of Guinea during November. This SST pattern exhibits a remarkably high predictive potential for the upcoming Atlantic hurricane activity. In the final chapter we investigate potential mechanisms, that link two of these regions with high feature relevance to the Atlantic hurricane activity. We mainly focus on ocean surface transport. The ocean surface flow paths are calculated using Lagrangian particle analysis. We demonstrate that the FR patterns in the region of the Canary islands do not correspond with ocean surface transport. It is instead likely that these FR patterns fingerprint a wind transport of latent heat. The second region to be studied is situated in the Gulf of Guinea. Our analysis shows that the FR patterns seen there do fingerprint ocean surface transport. However, our simulations also show that at least one other mechanism is involved in linking the Gulf of Guinea SST anomaly in November to the hurricane activity of the upcoming season. In this work the AnlgModel does not only demonstrate its outstanding forecast skills but also shows its capabilities as research tool for detecting oceanic and atmospheric mechanisms. KW - seasonal cyclone activity forecasts Y1 - 2023 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 - Wunderling, Nico A1 - Willeit, Matteo A1 - Donges, Jonathan A1 - Winkelmann, Ricarda T1 - Global warming due to loss of large ice masses and Arctic summer sea ice JF - Nature Communications N2 - 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. Y1 - 2020 U6 - https://doi.org/10.1038/s41467-020-18934-3 SN - 2041-1723 VL - 11 IS - 1 PB - Nature Publishing Group CY - Berlin 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 - Klose, Ann Kristin A1 - Wunderling, Nico A1 - Winkelmann, Ricarda A1 - Donges, Jonathan T1 - What do we mean, 'tipping cascade'? JF - Environmental research letters : ERL N2 - Based on suggested interactions of potential tipping elements in the Earth's climate and in ecological systems, tipping cascades as possible dynamics are increasingly discussed and studied. The activation of such tipping cascades would impose a considerable risk for human societies and biosphere integrity. However, there are ambiguities in the description of tipping cascades within the literature so far. Here we illustrate how different patterns of multiple tipping dynamics emerge from a very simple coupling of two previously studied idealized tipping elements. In particular, we distinguish between a two phase cascade, a domino cascade and a joint cascade. A mitigation of an unfolding two phase cascade may be possible and common early warning indicators are sensitive to upcoming critical transitions to a certain degree. In contrast, a domino cascade may hardly be stopped once initiated and critical slowing down-based indicators fail to indicate tipping of the following element. These different potentials for intervention and anticipation across the distinct patterns of multiple tipping dynamics should be seen as a call to be more precise in future analyses of cascading dynamics arising from tipping element interactions in the Earth system. KW - tipping cascade KW - domino effect KW - tipping interactions KW - cascading regime KW - shifts KW - early warning indicators Y1 - 2021 U6 - https://doi.org/10.1088/1748-9326/ac3955 SN - 1748-9326 VL - 16 IS - 12 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Durand, Gael A1 - van den Broeke, Michiel R. A1 - Le Cozannet, Goneri A1 - Edwards, Tamsin L. A1 - Holland, Paul R. A1 - Jourdain, Nicolas C. A1 - Marzeion, Ben A1 - Mottram, Ruth A1 - Nicholls, Robert J. A1 - Pattyn, Frank A1 - Paul, Frank A1 - Slangen, Aimee B. A. A1 - Winkelmann, Ricarda A1 - Burgard, Clara A1 - van Calcar, Caroline J. A1 - Barre, Jean-Baptiste A1 - Bataille, Amelie A1 - Chapuis, Anne T1 - Sea-Level rise: from global perspectives to local services JF - Frontiers in Marine Science N2 - Coastal areas are highly diverse, ecologically rich, regions of key socio-economic activity, and are particularly sensitive to sea-level change. Over most of the 20th century, global mean sea level has risen mainly due to warming and subsequent expansion of the upper ocean layers as well as the melting of glaciers and ice caps. Over the last three decades, increased mass loss of the Greenland and Antarctic ice sheets has also started to contribute significantly to contemporary sea-level rise. The future mass loss of the two ice sheets, which combined represent a sea-level rise potential of similar to 65 m, constitutes the main source of uncertainty in long-term (centennial to millennial) sea-level rise projections. Improved knowledge of the magnitude and rate of future sea-level change is therefore of utmost importance. Moreover, sea level does not change uniformly across the globe and can differ greatly at both regional and local scales. The most appropriate and feasible sea level mitigation and adaptation measures in coastal regions strongly depend on local land use and associated risk aversion. Here, we advocate that addressing the problem of future sea-level rise and its impacts requires (i) bringing together a transdisciplinary scientific community, from climate and cryospheric scientists to coastal impact specialists, and (ii) interacting closely and iteratively with users and local stakeholders to co-design and co-build coastal climate services, including addressing the high-end risks. KW - sea-level rise KW - Antarctic KW - Greenland KW - glaciers KW - local impact Y1 - 2022 U6 - https://doi.org/10.3389/fmars.2021.709595 SN - 2296-7745 VL - 8 PB - Frontiers Media CY - Lausanne ER - TY - JOUR A1 - Shin, Jaeoh A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Self-subdiffusion in solutions of star-shaped crowders: non-monotonic effects of inter-particle interactions JF - New journal of physics : the open-access journal for physics N2 - We examine by extensive computer simulations the self-diffusion of anisotropic star-like particles in crowded two-dimensional solutions. We investigate the implications of the area coverage fraction phi of the crowders and the crowder-crowder adhesion properties on the regime of transient anomalous diffusion. We systematically compute the mean squared displacement (MSD) of the particles, their time averaged MSD, and the effective diffusion coefficient. The diffusion is ergodic in the limit of long traces, such that the mean time averaged MSD converges towards the ensemble averaged MSD, and features a small residual amplitude spread of the time averaged MSD from individual trajectories. At intermediate time scales, we quantify the anomalous diffusion in the system. Also, we show that the translational-but not rotational-diffusivity of the particles Dis a nonmonotonic function of the attraction strength between them. Both diffusion coefficients decrease as the power law D(phi) similar to (1 - phi/phi*)(2 ... 2.4) with the area fraction phi occupied by the crowders and the critical value phi*. Our results might be applicable to rationalising the experimental observations of non-Brownian diffusion for a number of standard macromolecular crowders used in vitro to mimic the cytoplasmic conditions of living cells. KW - anomalous diffusion KW - crowded fluids KW - stochastic processes Y1 - 2015 U6 - https://doi.org/10.1088/1367-2630/17/11/113028 SN - 1367-2630 VL - 17 PB - IOP Publ. Ltd. CY - Bristol ER -