@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{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{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{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} } @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} } @article{WinkelmannLevermannMartinetal.2012, author = {Winkelmann, Ricarda and Levermann, Anders and Martin, Maria A. and Frieler, Katja}, title = {Increased future ice discharge from Antarctica owing to higher snowfall}, series = {Nature : the international weekly journal of science}, volume = {492}, journal = {Nature : the international weekly journal of science}, number = {7428}, publisher = {Nature Publ. Group}, address = {London}, issn = {0028-0836}, doi = {10.1038/nature11616}, pages = {239 -- +}, year = {2012}, abstract = {Anthropogenic climate change is likely to cause continuing global sea level rise(1), but some processes within the Earth system may mitigate the magnitude of the projected effect. Regional and global climate models simulate enhanced snowfall over Antarctica, which would provide a direct offset of the future contribution to global sea level rise from cryospheric mass loss(2,3) and ocean expansion(4). Uncertainties exist in modelled snowfall(5), but even larger uncertainties exist in the potential changes of dynamic ice discharge from Antarctica(1,6) and thus in the ultimate fate of the precipitation-deposited ice mass. Here we show that snowfall and discharge are not independent, but that future ice discharge will increase by up to three times as a result of additional snowfall under global warming. Our results, based on an ice-sheet model(7) forced by climate simulations through to the end of 2500 (ref. 8), show that the enhanced discharge effect exceeds the effect of surface warming as well as that of basal ice-shelf melting, and is due to the difference in surface elevation change caused by snowfall on grounded versus floating ice. Although different underlying forcings drive ice loss from basal melting versus increased snowfall, similar ice dynamical processes are nonetheless at work in both; therefore results are relatively independent of the specific representation of the transition zone. In an ensemble of simulations designed to capture ice-physics uncertainty, the additional dynamic ice loss along the coastline compensates between 30 and 65 per cent of the ice gain due to enhanced snowfall over the entire continent. This results in a dynamic ice loss of up to 1.25 metres in the year 2500 for the strongest warming scenario. The reported effect thus strongly counters a potential negative contribution to global sea level by the Antarctic Ice Sheet.}, language = {en} } @article{LevermannAlbrechtWinkelmannetal.2012, author = {Levermann, Anders and Albrecht, Tanja and Winkelmann, Ricarda and Martin, Maria A. and Haseloff, Monika and Joughin, I.}, title = {Kinematic first-order calving law implies potential for abrupt ice-shelf retreat}, series = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, volume = {6}, journal = {The Cryosphere : TC ; an interactive open access journal of the European Geosciences Union}, number = {2}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-6-273-2012}, pages = {273 -- 286}, year = {2012}, abstract = {Recently observed large-scale disintegration of Antarctic ice shelves has moved their fronts closer towards grounded ice. In response, ice-sheet discharge into the ocean has accelerated, contributing to global sea-level rise and emphasizing the importance of calving-front dynamics. The position of the ice front strongly influences the stress field within the entire sheet-shelf-system and thereby the mass flow across the grounding line. While theories for an advance of the ice-front are readily available, no general rule exists for its retreat, making it difficult to incorporate the retreat in predictive models. Here we extract the first-order large-scale kinematic contribution to calving which is consistent with large-scale observation. We emphasize that the proposed equation does not constitute a comprehensive calving law but represents the first-order kinematic contribution which can and should be complemented by higher order contributions as well as the influence of potentially heterogeneous material properties of the ice. When applied as a calving law, the equation naturally incorporates the stabilizing effect of pinning points and inhibits ice shelf growth outside of embayments. It depends only on local ice properties which are, however, determined by the full topography of the ice shelf. In numerical simulations the parameterization reproduces multiple stable fronts as observed for the Larsen A and B Ice Shelves including abrupt transitions between them which may be caused by localized ice weaknesses. We also find multiple stable states of the Ross Ice Shelf at the gateway of the West Antarctic Ice Sheet with back stresses onto the sheet reduced by up to 90 \% compared to the present state.}, 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} } @article{AlbrechtMartinHaseloffetal.2011, author = {Albrecht, Tanja and Martin, M. and Haseloff, M. and Winkelmann, Ricarda and Levermann, Anders}, title = {Parameterization for subgrid-scale motion of ice-shelf calving fronts}, 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 = {1}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1994-0416}, doi = {10.5194/tc-5-35-2011}, pages = {35 -- 44}, year = {2011}, abstract = {A parameterization for the motion of ice-shelf fronts on a Cartesian grid in finite-difference land-ice models is presented. The scheme prevents artificial thinning of the ice shelf at its edge, which occurs due to the finite resolution of the model. The intuitive numerical implementation diminishes numerical dispersion at the ice front and enables the application of physical boundary conditions to improve the calculation of stress and velocity fields throughout the ice-sheet-shelf system. Numerical properties of this subgrid modification are assessed in the Potsdam Parallel Ice Sheet Model (PISM-PIK) for different geometries in one and two horizontal dimensions and are verified against an analytical solution in a flow-line setup.}, language = {en} }