Search Results

Now showing 1 - 10 of 29
  • Item
    Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM) – Part 2: Parameter ensemble analysis
    (Katlenburg-Lindau : Copernicus, 2020) Albrecht, Torsten; Winkelmann, Ricarda; Levermann, Anders
    The Parallel Ice Sheet Model (PISM) is applied to the Antarctic Ice Sheet over the last two glacial cycles (≈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.1 m (or 6.5±2.0×106km3 ), 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.
  • Item
    Grounding-line flux formula applied as a flux condition in numerical simulations fails for buttressed Antarctic ice streams
    (Katlenburg-Lindau : Copernicus, 2018) Reese, Ronja; Winkelmann, Ricarda; Gudmundsson, G. Hilmar
    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.
  • Item
    Shear-margin melting causes stronger transient ice discharge than ice-stream melting in idealized simulations
    (Katlenburg-Lindau : Copernicus, 2022) Feldmann, Johannes; Reese, Ronja; Winkelmann, Ricarda; Levermann, Anders
    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.
  • Item
    Projecting Antarctica's contribution to future sea level rise from basal ice shelf melt using linear response functions of 16 ice sheet models (LARMIP-2)
    (Göttingen : Copernicus Publ., 2020) Levermann, Anders; Winkelmann, Ricarda; Albrecht, Torsten; Goelzer, Heiko; Golledge, Nicholas R.; Greve, Ralf; Huybrechts, Philippe; Jordan, Jim; Leguy, Gunter; Martin, Daniel; Morlighem, Mathieu; Pattyn, Frank; Pollard, David; Quiquet, Aurelien; Rodehacke, Christian; Seroussi, Helene; Sutter, Johannes; Zhang, Tong; Van Breedam, Jonas; Calov, Reinhard; DeConto, Robert; Dumas, Christophe; Garbe, Julius; Gudmundsson, G. Hilmar; Hoffman, Matthew J.; Humbert, Angelika; Kleiner, Thomas; Lipscomb, William H.; Meinshausen, Malte; Ng, Esmond; Nowicki, Sophie M.J.; Perego, Mauro; Price, Stephen F.; Saito, Fuyuki; Schlegel, Nicole-Jeanne; Sun, Sainan; van de Wal, Roderik S.W.
    The sea level contribution of the Antarctic ice sheet constitutes a large uncertainty in future sea level projections. Here we apply a linear response theory approach to 16 state-of-the-art ice sheet models to estimate the Antarctic ice sheet contribution from basal ice shelf melting within the 21st century. The purpose of this computation is to estimate the uncertainty of Antarctica's future contribution to global sea level rise that arises from large uncertainty in the oceanic forcing and the associated ice shelf melting. Ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean. However, by computing only the sea level contribution in response to ice shelf melting, our study is neglecting a number of processes such as surface-mass-balance-related contributions. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any self-dampening or self-amplifying processes. This is particularly relevant in situations in which an instability is dominating the ice loss. The results obtained here are thus relevant, in particular wherever the ice loss is dominated by the forcing as opposed to an internal instability, for example in strong ocean warming scenarios. In order to allow for comparison the methodology was chosen to be exactly the same as in an earlier study (Levermann et al., 2014) but with 16 instead of 5 ice sheet models. We include uncertainty in the atmospheric warming response to carbon emissions (full range of CMIP5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP5 models in relation to the global mean surface warming), and the observed range of responses of basal ice shelf melting to oceanic warming outside the ice shelf cavity. This uncertainty in basal ice shelf melting is then convoluted with the linear response functions of each of the 16 ice sheet models to obtain the ice flow response to the individual global warming path. The model median for the observational period from 1992 to 2017 of the ice loss due to basal ice shelf melting is 10.2 mm, with a likely range between 5.2 and 21.3 mm. For the same period the Antarctic ice sheet lost mass equivalent to 7.4mm of global sea level rise, with a standard deviation of 3.7mm (Shepherd et al., 2018) including all processes, especially surface-mass-balance changes. For the unabated warming path, Representative Concentration Pathway 8.5 (RCP8.5), we obtain a median contribution of the Antarctic ice sheet to global mean sea level rise from basal ice shelf melting within the 21st century of 17 cm, with a likely range (66th percentile around the mean) between 9 and 36 cm and a very likely range (90th percentile around the mean) between 6 and 58 cm. For the RCP2.6 warming path, which will keep the global mean temperature below 2 °C of global warming and is thus consistent with the Paris Climate Agreement, the procedure yields a median of 13 cm of global mean sea level contribution. The likely range for the RCP2.6 scenario is between 7 and 24 cm, and the very likely range is between 4 and 37 cm. The structural uncertainties in the method do not allow for an interpretation of any higher uncertainty percentiles.We provide projections for the five Antarctic regions and for each model and each scenario separately. The rate of sea level contribution is highest under the RCP8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade, with a likely range between 2 and 9 cm per decade and a very likely range between 1 and 14 cm per decade. © Author(s) 2020.
  • Item
    Interacting tipping elements increase risk of climate domino effects under global warming
    (Göttingen : Copernicus, 2021) Wunderling, Nico; Donges, Jonathan F.; Kurths, Jürgen; Winkelmann, Ricarda
    With progressing global warming, there is an increased risk that one or several tipping elements in the climate system might cross a critical threshold, resulting in severe consequences for the global climate, ecosystems and human societies. While the underlying processes are fairly well-understood, it is unclear how their interactions might impact the overall stability of the Earth's climate system. As of yet, this cannot be fully analysed with state-of-the-art Earth system models due to computational constraints as well as some missing and uncertain process representations of certain tipping elements. Here, we explicitly study the effects of known physical interactions among the Greenland and West Antarctic ice sheets, the Atlantic Meridional Overturning Circulation (AMOC) and the Amazon rainforest using a conceptual network approach. We analyse the risk of domino effects being triggered by each of the individual tipping elements under global warming in equilibrium experiments. In these experiments, we propagate the uncertainties in critical temperature thresholds, interaction strengths and interaction structure via large ensembles of simulations in a Monte Carlo approach. Overall, we find that the interactions tend to destabilise the network of tipping elements. Furthermore, our analysis reveals the qualitative role of each of the four tipping elements within the network, showing that the polar ice sheets on Greenland and West Antarctica are oftentimes the initiators of tipping cascades, while the AMOC acts as a mediator transmitting cascades. This indicates that the ice sheets, which are already at risk of transgressing their temperature thresholds within the Paris range of 1.5 to 2 ∘C, are of particular importance for the stability of the climate system as a whole.
  • Item
    Dynamic regimes of the Greenland Ice Sheet emerging from interacting melt–elevation and glacial isostatic adjustment feedbacks
    (Göttingen : Copernicus Publ., 2022) Zeitz, Maria; Haacker, Jan M.; Donges, Jonathan F.; Albrecht, Torsten; Winkelmann, Ricarda
    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 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 50000years, 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 74000 and over 300000 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 100000years 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.
  • Item
    Sensitivity of ice loss to uncertainty in flow law parameters in an idealized one-dimensional geometry
    (Katlenburg-Lindau : Copernicus, 2021) Zeitz, Maria; Levermann, Anders; Winkelmann, Ricarda
    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.
  • Item
    initMIP-Antarctica: an ice sheet model initialization experiment of ISMIP6
    (Katlenburg-Lindau : Copernicus, 2019) Seroussi, Hélène; Nowicki, Sophie; Simon, Erika; Abe-Ouchi, Ayako; Albrecht, Torsten; Brondex, Julien; Cornford, Stephen; Dumas, Christophe; Gillet-Chaulet, Fabien; Goelzer, Heiko; Golledge, Nicholas R.; Gregory, Jonathan M.; Greve, Ralf; Hoffman, Matthew J.; Humbert, Angelika; Huybrechts, Philippe; Kleiner, Thomas; Larour, Eric; Leguy, Gunter; Lipscomb, William H.; Lowry, Daniel; Mengel, Matthias; Morlighem, Mathieu; Pattyn, Frank; Payne, Anthony J.; Pollard, David; Price, Stephen F.; Quiquet, Aurélien; Reerink, Thomas J.; Reese, Ronja; Rodehacke, Christian B.; Schlegel, Nicole-Jeanne; Shepherd, Andrew; Sun, Sainan; Sutter, Johannes; Van Breedam, Jonas; van de Wal, Roderik S. W.; Winkelmann, Ricarda; Zhang, Tong
    Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.
  • Item
    Effects of extreme melt events on ice flow and sea level rise of the Greenland Ice Sheet
    (Katlenburg-Lindau : Copernicus, 2023) Beckmann, Johanna; Winkelmann, Ricarda
    Over the past decade, Greenland has experienced several extreme melt events, the most pronounced ones in the years 2010, 2012 and 2019. With progressing climate change, such extreme melt events can be expected to occur more frequently and potentially become more severe and persistent. So far, however, projections of ice loss and sea level change from Greenland typically rely on scenarios which only take gradual changes in the climate into account. Using the Parallel Ice Sheet Model (PISM), we investigate the effect of extreme melt events on the overall mass balance of the Greenland Ice Sheet and the changes in ice flow, invoked by the altered surface topography. As a first constraint, this study estimates the overall effect of extreme melt events on the cumulative mass loss of the Greenland Ice Sheet. We find that the sea level contribution from Greenland might increase by 2 to 45 cm (0.2 % to 14 %) by the year 2300 if extreme events occur more frequently in the future under a Representative Concentration Pathway 8.5 (RCP8.5) scenario, and the ice sheet area might be reduced by an additional 6000 to 26 000 km2 by 2300 in comparison to future warming scenarios without extremes. In conclusion, projecting the future sea level contribution from the Greenland Ice Sheet requires consideration of the changes in both the frequency and intensity of extreme events. It is crucial to individually address these extremes at a monthly resolution as temperature forcing with the same excess temperature but evenly distributed over longer timescales (e.g., seasonal) leads to less sea level rise than for the simulations of the resolved extremes. Extremes lead to additional mass loss and thinning. This, in turn, reduces the driving stress and surface velocities, ultimately dampening the ice loss attributed to ice flow and discharge. Overall, we find that the surface elevation feedback largely amplifies melting for scenarios with and without extremes, with additional mass loss attributed to this feedback having the greatest impact on projected sea level.
  • Item
    ISMIP6 Antarctica: A multi-model ensemble of the Antarctic ice sheet evolution over the 21st century
    (Katlenburg-Lindau : Copernicus, 2020) Seroussi, Hélène; Nowicki, Sophie; Payne, Antony J.; Goelzer, Heiko; Lipscomb, William H.; Abe-Ouchi, Ayako; Agosta, Cécile; Albrecht, Torsten; Asay-Davis, Xylar; Barthel, Alice; Calov, Reinhard; Cullather, Richard; Dumas, Christophe; Galton-Fenzi, Benjamin K.; Gladstone, Rupert; Golledge, Nicholas R.; Gregory, Jonathan M.; Greve, Ralf; Hattermann, Tore; Hoffman, Matthew J.; Humbert, Angelika; Huybrechts, Philippe; Jourdain, Nicolas C.; Kleiner, Thomas; Larour, Eric; Leguy, Gunter R.; Lowry, Daniel P.; Little, Chistopher M.; Morlighem, Mathieu; Pattyn, Frank; Pelle, Tyler; Price, Stephen F.; Quiquet, Aurélien; Reese, Ronja; Schlegel, Nicole-Jeanne; Shepherd, Andrew; Simon, Erika; Smith, Robin S.; Straneo, Fiammetta; Sun, Sainan; Trusel, Luke D.; Van Breedam, Jonas; van de Wal, Roderik S. W.; Winkelmann, Ricarda; Zhao, Chen; Zhang, Tong; Zwinger, Thomas
    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 presentday 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. © Author(s) 2020.