@article{ReeseAlbrechtMengeletal.2018, author = {Reese, Ronja and Albrecht, Torsten and Mengel, Matthias and Asay-Davis, Xylar and Winkelmann, Ricarda}, title = {Antarctic sub-shelf melt rates via PICO}, 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 = {6}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-12-1969-2018}, pages = {1969 -- 1985}, year = {2018}, abstract = {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.}, language = {en} } @article{LevermannWinkelmannNowickietal.2014, author = {Levermann, Anders and Winkelmann, Ricarda and Nowicki, S. and Fastook, J. L. and Frieler, Katja and Greve, R. and Hellmer, H. H. and Martin, M. A. and Meinshausen, Malte and Mengel, Matthias and Payne, A. J. and Pollard, D. and Sato, T. and Timmermann, R. and Wang, Wei Li and Bindschadler, Robert A.}, title = {Projecting antarctic ice discharge using response functions from SeaRISE ice-sheet models}, series = {Earth system dynamics}, volume = {5}, journal = {Earth system dynamics}, number = {2}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {2190-4979}, doi = {10.5194/esd-5-271-2014}, pages = {271 -- 293}, year = {2014}, abstract = {The largest uncertainty in projections of future sea-level change results from the potentially changing dynamical ice discharge from Antarctica. Basal ice-shelf melting induced by a warming ocean has been identified as a major cause for additional ice flow across the grounding line. Here we attempt to estimate the uncertainty range of future ice discharge from Antarctica by combining uncertainty in the climatic forcing, the oceanic response and the ice-sheet model response. The uncertainty in the global mean temperature increase is obtained from historically constrained emulations with the MAGICC-6.0 (Model for the Assessment of Greenhouse gas Induced Climate Change) model. The oceanic forcing is derived from scaling of the subsurface with the atmospheric warming from 19 comprehensive climate models of the Coupled Model Intercomparison Project (CMIP-5) and two ocean models from the EU-project Ice2Sea. The dynamic ice-sheet response is derived from linear response functions for basal ice-shelf melting for four different Antarctic drainage regions using experiments from the Sea-level Response to Ice Sheet Evolution (SeaRISE) intercomparison project with five different Antarctic ice-sheet models. The resulting uncertainty range for the historic Antarctic contribution to global sea-level rise from 1992 to 2011 agrees with the observed contribution for this period if we use the three ice-sheet models with an explicit representation of ice-shelf dynamics and account for the time-delayed warming of the oceanic subsurface compared to the surface air temperature. The median of the additional ice loss for the 21st century is computed to 0.07 m (66\% range: 0.02-0.14 m; 90\% range: 0.0-0.23 m) of global sea-level equivalent for the low-emission RCP-2.6 (Representative Concentration Pathway) scenario and 0.09 m (66\% range: 0.04-0.21 m; 90\% range: 0.01-0.37 m) for the strongest RCP-8.5. Assuming no time delay between the atmospheric warming and the oceanic subsurface, these values increase to 0.09 m (66\% range: 0.04-0.17 m; 90\% range: 0.02-0.25 m) for RCP-2.6 and 0.15 m (66\% range: 0.07-0.28 m; 90\% range: 0.04-0.43 m) for RCP-8.5. All probability distributions are highly skewed towards high values. The applied ice-sheet models are coarse resolution with limitations in the representation of grounding-line motion. Within the constraints of the applied methods, the uncertainty induced from different ice-sheet models is smaller than that induced by the external forcing to the ice sheets.}, language = {en} } @article{SteffenRoeckstromRichardsonetal.2018, author = {Steffen, Will and R{\"o}ckstrom, Johan and Richardson, Katherine and Lenton, Timothy M. and Folke, Carl and Liverman, Diana and Summerhayes, Colin P. and Barnosky, Anthony D. and Cornell, Sarah E. and Crucifix, Michel and Donges, Jonathan and Fetzer, Ingo and Lade, Steven J. and Scheffer, Marten and Winkelmann, Ricarda and Schellnhuber, Hans Joachim}, title = {Trajectories of the Earth System in the Anthropocene}, series = {Proceedings of the National Academy of Sciences of the United States of America}, volume = {115}, journal = {Proceedings of the National Academy of Sciences of the United States of America}, number = {33}, publisher = {National Acad. of Sciences}, address = {Washington}, issn = {0027-8424}, doi = {10.1073/pnas.1810141115}, pages = {8252 -- 8259}, year = {2018}, abstract = {We explore the risk that self-reinforcing feedbacks could push the Earth System toward a planetary threshold that, if crossed, could prevent stabilization of the climate at intermediate temperature rises and cause continued warming on a "Hothouse Earth" pathway even as human emissions are reduced. Crossing the threshold would lead to a much higher global average temperature than any interglacial in the past 1.2 million years and to sea levels significantly higher than at any time in the Holocene. We examine the evidence that such a threshold might exist and where it might be. If the threshold is crossed, the resulting trajectory would likely cause serious disruptions to ecosystems, society, and economies. Collective human action is required to steer the Earth System away from a potential threshold and stabilize it in a habitable interglacial-like state. Such action entails stewardship of the entire Earth System-biosphere, climate, and societies-and could include decarbonization of the global economy, enhancement of biosphere carbon sinks, behavioral changes, technological innovations, new governance arrangements, and transformed social values.}, language = {en} } @article{ClarkShakunMarcottetal.2016, author = {Clark, Peter U. and Shakun, Jeremy D. and Marcott, Shaun A. and Mix, Alan C. and Eby, Michael and Kulp, Scott and Levermann, Anders and Milne, Glenn A. and Pfister, Patrik L. and Santer, Benjamin D. and Schrag, Daniel P. and Solomon, Susan and Stocker, Thomas F. and Strauss, Benjamin H. and Weaver, Andrew J. and Winkelmann, Ricarda and Archer, David and Bard, Edouard and Goldner, Aaron and Lambeck, Kurt and Pierrehumbert, Raymond T. and Plattner, Gian-Kasper}, title = {Consequences of twenty-first-century policy for multi-millennial climate and sea-level change}, series = {Nature climate change}, volume = {6}, journal = {Nature climate change}, publisher = {Nature Publ. Group}, address = {London}, issn = {1758-678X}, doi = {10.1038/NCLIMATE2923}, pages = {360 -- 369}, year = {2016}, abstract = {Most of the policy debate surrounding the actions needed to mitigate and adapt to anthropogenic climate change has been framed by observations of the past 150 years as well as climate and sea-level projections for the twenty-first century. The focus on this 250-year window, however, obscures some of the most profound problems associated with climate change. Here, we argue that the twentieth and twenty-first centuries, a period during which the overwhelming majority of human-caused carbon emissions are likely to occur, need to be placed into a long-term context that includes the past 20 millennia, when the last Ice Age ended and human civilization developed, and the next ten millennia, over which time the projected impacts of anthropogenic climate change will grow and persist. This long-term perspective illustrates that policy decisions made in the next few years to decades will have profound impacts on global climate, ecosystems and human societies - not just for this century, but for the next ten millennia and beyond.}, 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} } @article{GanopolskiWinkelmannSchellnhuber2016, author = {Ganopolski, A. and Winkelmann, Ricarda and Schellnhuber, Hans Joachim}, title = {Critical insolation-CO2 relation for diagnosing past and future glacial inception}, series = {Nature : the international weekly journal of science}, volume = {529}, journal = {Nature : the international weekly journal of science}, publisher = {Nature Publ. Group}, address = {London}, issn = {0028-0836}, doi = {10.1038/nature16494}, pages = {200 -- U159}, year = {2016}, abstract = {The past rapid growth of Northern Hemisphere continental ice sheets, which terminated warm and stable climate periods, is generally attributed to reduced summer insolation in boreal latitudes(1-3). Yet such summer insolation is near to its minimum at present(4), and there are no signs of a new ice age(5). This challenges our understanding of the mechanisms driving glacial cycles and our ability to predict the next glacial inception(6). Here we propose a critical functional relationship between boreal summer insolation and global carbon dioxide (CO2) concentration, which explains the beginning of the past eight glacial cycles and might anticipate future periods of glacial inception. Using an ensemble of simulations generated by an Earth system model of intermediate complexity constrained by palaeoclimatic data, we suggest that glacial inception was narrowly missed before the beginning of the Industrial Revolution. The missed inception can be accounted for by the combined effect of relatively high late-Holocene CO2 concentrations and the low orbital eccentricity of the Earth(7). Additionally, our analysis suggests that even in the absence of human perturbations no substantial build-up of ice sheets would occur within the next several thousand years and that the current interglacial would probably last for another 50,000 years. However, moderate anthropogenic cumulative CO2 emissions of 1,000 to 1,500 gigatonnes of carbon will postpone the next glacial inception by at least 100,000 years(8,9). Our simulations demonstrate that under natural conditions alone the Earth system would be expected to remain in the present delicately balanced interglacial climate state, steering clear of both large-scale glaciation of the Northern Hemisphere and its complete deglaciation, for an unusually long time.}, 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} } @article{FrielerClarkHeetal.2015, author = {Frieler, Katja and Clark, Peter U. and He, Feng and Buizert, Christo and Reese, Ronja and Ligtenberg, Stefan R. M. and van den Broeke, Michiel R. and Winkelmann, Ricarda and Levermann, Anders}, title = {Consistent evidence of increasing Antarctic accumulation with warming}, series = {Nature climate change}, volume = {5}, journal = {Nature climate change}, number = {4}, publisher = {Nature Publ. Group}, address = {London}, issn = {1758-678X}, doi = {10.1038/nclimate2574}, pages = {348 -- 352}, year = {2015}, abstract = {Projections of changes in Antarctic Ice Sheet (AIS) surface mass balance indicate a negative contribution to sea level because of the expected increase in precipitation due to the higher moisture holding capacity of warmer air(1). Observations over the past decades, however, are unable to constrain the relation between temperature and accumulation changes because both are dominated by strong natural variability(2-5). Here we derive a consistent continental-scale increase in accumulation of approximately 5 +/- 1\% K-1, through the assessment of ice-core data (spanning the large temperature change during the last deglaciation, 21,000 to 10,000 years ago), in combination with palaeo-simulations, future projections by 35 general circulation models (GCMs), and one high-resolution future simulation. The ice-core data and modelling results for the last deglaciation agree, showing uniform local sensitivities of similar to 6\% K-1. The palaeo-simulation allows for a continental-scale aggregation of accumulation changes reaching 4.3\% K-1. Despite the different timescales, these sensitivities agree with the multi-model mean of 6.1 +/- 2.6\% K-1 (GCMprojections) and the continental-scale sensitivity of 4.9\% K-1 (high-resolution future simulation). Because some of the mass gain of the AIS is offset by dynamical losses induced by accumulation(6,7), we provide a response function allowing projections of sea-level fall in terms of continental-scale accumulation changes that compete with surface melting and dynamical losses induced by other mechanisms(6,8,9).}, 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{WinkelmannMartinHaseloffetal.2011, author = {Winkelmann, Ricarda and Martin, Maria A. and Haseloff, Monika and Albrecht, Torsten and Bueler, Ed and Khroulev, C. and Levermann, Anders}, title = {The Potsdam parallel ice sheet model (PISM-PIK) - Part 1: Model description}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {5}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {3}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-5-715-2011}, pages = {715 -- 726}, year = {2011}, abstract = {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).}, language = {en} } @article{MartinWinkelmannHaseloffetal.2011, author = {Martin, Maria A. and Winkelmann, Ricarda and Haseloff, M. and Albrecht, Tanja and Bueler, Ed and Khroulev, C. and Levermann, Anders}, title = {The Potsdam parallel ice sheet model (PISM-PIK) - Part 2: Dynamic equilibrium simulation of the Antarctic ice sheet}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {5}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {3}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-5-727-2011}, pages = {727 -- 740}, year = {2011}, abstract = {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.}, language = {en} } @phdthesis{Reese2018, author = {Reese, Ronja}, title = {The far reach of ice-shelf thinning in Antarctica}, school = {Universit{\"a}t Potsdam}, pages = {227}, year = {2018}, language = {en} } @phdthesis{Winkelmann2012, author = {Winkelmann, Ricarda}, title = {The future sea-level contribution from antartica: projections of solid ice discharge}, address = {Potsdam}, pages = {140 S.}, year = {2012}, language = {en} } @phdthesis{Wunderling2021, author = {Wunderling, Nico}, title = {Nichtlineare Dynamiken und Interaktionen von Kippelementen im Erdsystem}, doi = {10.25932/publishup-52514}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-525140}, school = {Universit{\"a}t Potsdam}, pages = {ix, 303}, year = {2021}, abstract = {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{\"o}s-R{\´e}nyi, Barab{\´a}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.}, 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{ReeseLevermannAlbrechtetal.2020, author = {Reese, Ronja and Levermann, Anders and Albrecht, Torsten and Seroussi, Helene and Winkelmann, Ricarda}, title = {The role of history and strength of the oceanic forcing in sea level projections from Antarctica with the Parallel Ice Sheet Model}, 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-3097-2020}, pages = {3097 -- 3110}, year = {2020}, abstract = {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.}, language = {en} } @article{ZeitzLevermannWinkelmann2020, author = {Zeitz, Maria and Levermann, Anders and Winkelmann, Ricarda}, title = {Sensitivity of ice loss to uncertainty in flow law parameters in an idealized one-dimensional geometry}, 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 = {10}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-14-3537-2020}, pages = {3537 -- 3550}, year = {2020}, abstract = {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.}, language = {en} } @article{ReeseGudmundssonLevermannetal.2017, author = {Reese, Ronja and Gudmundsson, Gudmundur Hilmar and Levermann, Anders and Winkelmann, Ricarda}, title = {The far reach of ice-shelf thinning in Antarctica}, series = {Nature climate change}, volume = {8}, journal = {Nature climate change}, number = {1}, publisher = {Nature Publ. Group}, address = {London}, issn = {1758-678X}, doi = {10.1038/s41558-017-0020-x}, pages = {53 -- 57}, year = {2017}, abstract = {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.}, language = {en} } @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} }