Search Results

Now showing 1 - 10 of 37
  • 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
    Exploring the impact of atmospheric forcing and basal drag on the Antarctic Ice Sheet under Last Glacial Maximum conditions
    (Katlenburg-Lindau : Copernicus, 2021-1-18) Blasco, Javier; Alvarez-Solas, Jorge; Robinson, Alexander; Montoya, Marisa
    Little is known about the distribution of ice in the Antarctic Ice Sheet (AIS) during the Last Glacial Maximum (LGM). Whereas marine and terrestrial geological data indicate that the grounded ice advanced to a position close to the continental-shelf break, the total ice volume is unclear. Glacial boundary conditions are potentially important sources of uncertainty, in particular basal friction and climatic boundary conditions. Basal friction exerts a strong control on the large-scale dynamics of the ice sheet and thus affects its size and is not well constrained. Glacial climatic boundary conditions determine the net accumulation and ice temperature and are also poorly known. Here we explore the effect of the uncertainty in both features on the total simulated ice storage of the AIS at the LGM. For this purpose we use a hybrid ice sheet shelf model that is forced with different basal drag choices and glacial background climatic conditions obtained from the LGM ensemble climate simulations of the third phase of the Paleoclimate Modelling Intercomparison Project (PMIP3). Overall, we find that the spread in the simulated ice volume for the tested basal drag parameterizations is about the same range as for the different general circulation model (GCM) forcings (4 to 6 m sea level equivalent). For a wide range of plausible basal friction configurations, the simulated ice dynamics vary widely but all simulations produce fully extended ice sheets towards the continental-shelf break. More dynamically active ice sheets correspond to lower ice volumes, while they remain consistent with the available constraints on ice extent. Thus, this work points to the possibility of an AIS with very active ice streams during the LGM. In addition, we find that the surface boundary temperature field plays a crucial role in determining the ice extent through its effect on viscosity. For ice sheets of a similar extent and comparable dynamics, we find that the precipitation field determines the total AIS volume. However, precipitation is highly uncertain. Climatic fields simulated by climate models show more precipitation in coastal regions than a spatially uniform anomaly, which can lead to larger ice volumes. Our results strongly support using these paleoclimatic fields to simulate and study the LGM and potentially other time periods like the last interglacial. However, their accuracy must be assessed as well, as differences between climate model forcing lead to a large spread in the simulated ice volume and extension.
  • Item
    Drivers of Pine Island Glacier speed-up between 1996 and 2016
    (Katlenburg-Lindau : Copernicus, 2021-1-7) De Rydt, Jan; Reese, Ronja; Paolo, Fernando S.; Gudmundsson, G. Hilmar
    Pine Island Glacier in West Antarctica is among the fastest changing glaciers worldwide. Over the last 2 decades, the glacier has lost in excess of a trillion tons of ice, or the equivalent of 3 mm of sea level rise. The ongoing changes are thought to have been triggered by ocean-induced thinning of its floating ice shelf, grounding line retreat, and the associated reduction in buttressing forces. However, other drivers of change, such as large-scale calving and changes in ice rheology and basal slipperiness, could play a vital, yet unquantified, role in controlling the ongoing and future evolution of the glacier. In addition, recent studies have shown that mechanical properties of the bed are key to explaining the observed speed-up. Here we used a combination of the latest remote sensing datasets between 1996 and 2016, data assimilation tools, and numerical perturbation experiments to quantify the relative importance of all processes in driving the recent changes in Pine Island Glacier dynamics. We show that (1) calving and ice shelf thinning have caused a comparable reduction in ice shelf buttressing over the past 2 decades; that (2) simulated changes in ice flow over a viscously deforming bed are only compatible with observations if large and widespread changes in ice viscosity and/or basal slipperiness are taken into account; and that (3) a spatially varying, predominantly plastic bed rheology can closely reproduce observed changes in flow without marked variations in ice-internal and basal properties. Our results demonstrate that, in addition to its evolving ice thickness, calving processes and a heterogeneous bed rheology play a key role in the contemporary evolution of Pine Island Glacier.
  • Item
    Revisiting temperature sensitivity: how does Antarctic precipitation change with temperature?
    (Katlenburg-Lindau : EGU, 2023) Nicola, Lena; Notz, Dirk; Winkelmann, Ricarda
    With progressing global warming, snowfall in Antarctica is expected to increase, which could counteract or even temporarily overcompensate increased ice-sheet mass losses caused by increased ice discharge and melting. For sea-level projections it is therefore vital to understand the processes determining snowfall changes in Antarctica. Here we revisit the relationship between Antarctic temperature changes and precipitation changes, identifying and explaining regional differences and deviations from the theoretical approach based on the Clausius-Clapeyron relationship. Analysing the latest estimates from global (CMIP6, Coupled Model Intercomparison Project Phase 6) and regional (RACMO2.3) model projections, we find an average increase of 5.5 % in annual precipitation over Antarctica per degree of warming, with a minimum sensitivity of 2 % K-1 near Siple Coast and a maximum sensitivity of > 10 % K-1 at the East Antarctic plateau region. This large range can be explained by the prevailing climatic conditions, with local temperatures determining the Clausius-Clapeyron sensitivity that is counteracted in some regions by the prevalence of the coastal wind regime. We compare different approaches of deriving the sensitivity factor, which in some cases can lead to sensitivity changes of up to 7 percentage points for the same model. Importantly, local sensitivity factors are found to be strongly dependent on the warming level, suggesting that some ice-sheet models which base their precipitation estimates on parameterisations derived from these sensitivity factors might overestimate warming-induced snowfall changes, particularly in high-emission scenarios. This would have consequences for Antarctic sea-level projections for this century and beyond.
  • Item
    Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM) – Part 1: Boundary conditions and climatic forcing
    (Katlenburg-Lindau : Copernicus, 2020) Albrecht, Torsten; Winkelmann, Ricarda; Levermann, Anders
    Simulations of the glacial–interglacial history of the Antarctic Ice Sheet provide insights into dynamic threshold behavior and estimates of the ice sheet's contributions to global sea-level changes for the past, present and future. However, boundary conditions are weakly constrained, in particular at the interface of the ice sheet and the bedrock. Also climatic forcing covering the last glacial cycles is uncertain, as it is based on sparse proxy data. We use the Parallel Ice Sheet Model (PISM) to investigate the dynamic effects of different choices of input data, e.g., for modern basal heat flux or reconstructions of past changes of sea level and surface temperature. As computational resources are limited, glacial-cycle simulations are performed using a comparably coarse model grid of 16 km and various parameterizations, e.g., for basal sliding, iceberg calving, or for past variations in precipitation and ocean temperatures. In this study we evaluate the model's transient sensitivity to corresponding parameter choices and to different boundary conditions over the last two glacial cycles and provide estimates of involved uncertainties. We also discuss isolated and combined effects of climate and sea-level forcing. Hence, this study serves as a “cookbook” for the growing community of PISM users and paleo-ice sheet modelers in general. For each of the different model uncertainties with regard to climatic forcing, ice and Earth dynamics, and basal processes, we select one representative model parameter that captures relevant uncertainties and motivates corresponding parameter ranges that bound the observed ice volume at present. The four selected parameters are systematically varied in a parameter ensemble analysis, which is described in a companion paper.
  • Item
    The role of history and strength of the oceanic forcing in sea level projections from Antarctica with the Parallel Ice Sheet Model
    (Katlenburg-Lindau : Copernicus, 2020) Reese, Ronja; Levermann, Anders; Albrecht, Torsten; Seroussi, Hélène; Winkelmann, Ricarda
    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.
  • Item
    Timescales of outlet-glacier flow with negligible basal friction: Theory, observations and modeling
    (Katlenburg-Lindau : Copernicus, 2023) Feldmann, Johannes; Levermann, Anders
    The timescales of the flow and retreat of Greenland's and Antarctica's outlet glaciers and their potential instabilities are arguably the largest uncertainty in future sea-level projections. Here we derive a scaling relation that allows the comparison of the timescales of observed complex ice flow fields with geometric similarity. The scaling relation is derived under the assumption of fast, laterally confined, geometrically similar outlet-glacier flow over a slippery bed, i.e., with negligible basal friction. According to the relation, the time scaling of the outlet flow is determined by the product of the inverse of (1) the fourth power of the width-To-length ratio of its confinement, (2) the third power of the confinement depth and (3) the temperature-dependent ice softness. For the outflow at the grounding line of streams with negligible basal friction, this means that the volume flux is proportional to the ice softness and the bed depth, but goes with the fourth power of the gradient of the bed and with the fifth power of the width of the stream. We show that the theoretically derived scaling relation is supported by the observed velocity scaling of outlet glaciers across Greenland as well as by idealized numerical simulations of marine ice-sheet instabilities (MISIs) as found in Antarctica. Assuming that changes in the ice-flow velocity due to ice-dynamic imbalance are proportional to the equilibrium velocity, we combine the scaling relation with a statistical analysis of the topography of 13 MISI-prone Antarctic outlets. Under these assumptions, the timescales in response to a potential destabilization are fastest for Thwaites Glacier in West Antarctica and Mellor, Ninnis and Cook Glaciers in East Antarctica; between 16 and 67 times faster than for Pine Island Glacier. While the applicability of our results is limited by several strong assumptions, the utilization and potential further development of the presented scaling approach may help to constrain timescale estimates of outlet-glacier flow, augmenting the commonly exploited and comparatively computationally expensive approach of numerical modeling.
  • Item
    Stabilizing effect of mélange buttressing on the marine ice-cliff instability of the West Antarctic Ice Sheet
    (Katlenburg-Lindau : Copernicus, 2022) Schlemm, Tanja; Feldmann, Johannes; Winkelmann, Ricarda; Levermann, Anders
    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 mélange-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.6m within 100a. 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 mélange 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-mélange 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.
  • Item
    The tipping points and early warning indicators for Pine Island Glacier, West Antarctica
    (Katlenburg-Lindau : Copernicus, 2021-3-25) Rosier, Sebastian H. R.; Reese, Ronja; Donges, Jonathan F.; De Rydt, Jan; Gudmundsson, G. Hilmar; Winkelmann, Ricarda
    Mass loss from the Antarctic Ice Sheet is the main source of uncertainty in projections of future sea-level rise, with important implications for coastal regions worldwide. Central to ongoing and future changes is the marine ice sheet instability: once a critical threshold, or tipping point, is crossed, ice internal dynamics can drive a self-sustaining retreat committing a glacier to irreversible, rapid and substantial ice loss. This process might have already been triggered in the Amundsen Sea region, where Pine Island and Thwaites glaciers dominate the current mass loss from Antarctica, but modelling and observational techniques have not been able to establish this rigorously, leading to divergent views on the future mass loss of the West Antarctic Ice Sheet. Here, we aim at closing this knowledge gap by conducting a systematic investigation of the stability regime of Pine Island Glacier. To this end we show that early warning indicators in model simulations robustly detect the onset of the marine ice sheet instability. We are thereby able to identify three distinct tipping points in response to increases in ocean-induced melt. The third and final event, triggered by an ocean warming of approximately 1.2 ∘C from the steady-state model configuration, leads to a retreat of the entire glacier that could initiate a collapse of the West Antarctic Ice Sheet.
  • 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.