Refine
Has Fulltext
- no (14)
Year of publication
Document Type
- Article (13)
- Doctoral Thesis (1)
Language
- English (14)
Is part of the Bibliography
- yes (14)
Keywords
- Antarctica (1)
- Larsen Ice Shelf (1)
- glaciology (1)
- iceberg calving (1)
- numerical ice modeling (1)
- sea level (1)
Institute
- Institut für Physik und Astronomie (14) (remove)
Recent observations and modeling studies emphasize the crucial role of fracture mechanics for the stability of ice shelves and thereby the evolution of ice sheets. Here we introduce a macroscopic fracture-density field into a prognostic continuum ice-flow model and compute its evolution incorporating the initiation and growth of fractures as well as their advection with two-dimensional ice flow. To a first approximation, fracture growth is assumed to depend on the spreading rate only, while fracture initiation is defined in terms of principal stresses. The inferred fracture-density fields compare well with observed elongate surface structures. Since crevasses and other deep-reaching fracture structures have been shown to influence the overall ice-shelf dynamics, we propose the fracture-density field introduced here be used as a measure for ice softening and decoupling of the ice flow in fracture-weakened zones. This may yield more accurate and realistic velocity patterns in prognostic simulations. Additionally, the memory of past fracture events links the calving front to the upstream dynamics. Thus the fracture-density field proposed here may be employed in fracture-based calving parameterizations. The aim of this study is to introduce the field and investigate which of the observed surface structures can be reproduced by the simplest physically motivated fracture source terms.
Antarctic ice-discharge constitutes the largest uncertainty in future sea-level projections. Floating ice shelves, fringing most of Antarctica, exert retentive forces onto the ice flow. While abrupt ice-shelf retreat has been observed, it is generally considered a localized phenomenon. Here we show that the disintegration of an ice shelf may induce the spontaneous retreat of its neighbor. As an example, we reproduce the spontaneous but gradual retreat of the Larsen B ice front as observed after the disintegration of the adjacent Larsen A ice shelf. We show that the Larsen A collapse yields a change in spreading rate in Larsen B via their connecting ice channels and thereby causes a retreat of the ice front to its observed position of the year 2000, prior to its collapse. This mechanism might be particularly relevant for the role of East Antarctica and the Antarctic Peninsula in future sea level.
Floating ice shelves can exert a retentive and hence stabilizing force onto the inland ice sheet of Antarctica. However, this effect has been observed to diminish by the dynamic effects of fracture processes within the protective ice shelves, leading to accelerated ice flow and hence to a sea-level contribution. In order to account for the macroscopic effect of fracture processes on large-scale viscous ice dynamics (i.e., ice-shelf scale) we apply a continuum representation of fractures and related fracture growth into the prognostic Parallel Ice Sheet Model (PISM) and compare the results to observations. To this end we introduce a higher order accuracy advection scheme for the transport of the two-dimensional fracture density across the regular computational grid. Dynamic coupling of fractures and ice flow is attained by a reduction of effective ice viscosity proportional to the inferred fracture density. This formulation implies the possibility of non-linear threshold behavior due to self-amplified fracturing in shear regions triggered by small variations in the fracture-initiation threshold. As a result of prognostic flow simulations, sharp across-flow velocity gradients appear in fracture-weakened regions. These modeled gradients compare well in magnitude and location with those in observed flow patterns. This model framework is in principle expandable to grounded ice streams and provides simple means of investigating climate-induced effects on fracturing (e. g., hydro fracturing) and hence on the ice flow. It further constitutes a physically sound basis for an enhanced fracture-based calving parameterization.
Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM)
(2020)
The Parallel Ice Sheet Model (PISM) is applied to the Antarctic Ice Sheet over the last two glacial cycles (approximate to 210 000 years) with a resolution of 16 km. An ensemble of 256 model runs is analyzed in which four relevant model parameters have been systematically varied using full-factorial parameter sampling. Parameters and plausible parameter ranges have been identified in a companion paper (Albrecht et al., 2020) and are associated with ice dynamics, climatic forcing, basal sliding and bed deformation and represent distinct classes of model uncertainties. The model is scored against both modern and geologic data, including reconstructed grounding-line locations, elevation-age data, ice thickness, surface velocities and uplift rates. An aggregated score is computed for each ensemble member that measures the overall model-data misfit, including measurement uncertainty in terms of a Gaussian error model (Briggs and Tarasov, 2013). The statistical method used to analyze the ensemble simulation results follows closely the simple averaging method described in Pollard et al. (2016).
This analysis reveals clusters of best-fit parameter combinations, and hence a likely range of relevant model and boundary parameters, rather than individual best-fit parameters. The ensemble of reconstructed histories of Antarctic Ice Sheet volumes provides a score-weighted likely range of sea-level contributions since the Last Glacial Maximum (LGM) of 9.4 +/- 4.1m (or 6.5 +/- 2.0 x 10(6) km(3)), which is at the upper range of most previous studies. The last deglaciation occurs in all ensemble simulations after around 12 000 years before present and hence after the meltwater pulse 1A (MWP1a). Our ensemble analysis also provides an estimate of parametric uncertainty bounds for the present-day state that can be used for PISM projections of future sea-level contributions from the Antarctic Ice Sheet.
More than half of Earth's freshwater resources are held by the Antarctic Ice Sheet, which thus represents by far the largest potential source for global sea-level rise under future warming conditions(1). Its long-term stability determines the fate of our coastal cities and cultural heritage. Feedbacks between ice, atmosphere, ocean, and the solid Earth give rise to potential nonlinearities in its response to temperature changes. So far, we are lacking a comprehensive stability analysis of the Antarctic Ice Sheet for different amounts of global warming. Here we show that the Antarctic Ice Sheet exhibits a multitude of temperature thresholds beyond which ice loss is irreversible. Consistent with palaeodata(2)we find, using the Parallel Ice Sheet Model(3-5), that at global warming levels around 2 degrees Celsius above pre-industrial levels, West Antarctica is committed to long-term partial collapse owing to the marine ice-sheet instability. Between 6 and 9 degrees of warming above pre-industrial levels, the loss of more than 70 per cent of the present-day ice volume is triggered, mainly caused by the surface elevation feedback. At more than 10 degrees of warming above pre-industrial levels, Antarctica is committed to become virtually ice-free. The ice sheet's temperature sensitivity is 1.3 metres of sea-level equivalent per degree of warming up to 2 degrees above pre-industrial levels, almost doubling to 2.4 metres per degree of warming between 2 and 6 degrees and increasing to about 10 metres per degree of warming between 6 and 9 degrees. Each of these thresholds gives rise to hysteresis behaviour: that is, the currently observed ice-sheet configuration is not regained even if temperatures are reversed to present-day levels. In particular, the West Antarctic Ice Sheet does not regrow to its modern extent until temperatures are at least one degree Celsius lower than pre-industrial levels. Our results show that if the Paris Agreement is not met, Antarctica's long-term sea-level contribution will dramatically increase and exceed that of all other sources. <br /> Modelling shows that the Antarctic Ice Sheet exhibits multiple temperature thresholds beyond which ice loss would become irreversible, and once melted, the ice sheet can regain its previous mass only if the climate cools well below pre-industrial temperatures.
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.
Predictions of marine ice-sheet behaviour require models able to simulate grounding-line migration. We present results of an intercomparison experiment for plan-view marine ice-sheet models. Verification is effected by comparison with approximate analytical solutions for flux across the grounding line using simplified geometrical configurations (no lateral variations, no buttressing effects from lateral drag). Perturbation experiments specifying spatial variation in basal sliding parameters permitted the evolution of curved grounding lines, generating buttressing effects. The experiments showed regions of compression and extensional flow across the grounding line, thereby invalidating the boundary layer theory. Steady-state grounding-line positions were found to be dependent on the level of physical model approximation. Resolving grounding lines requires inclusion of membrane stresses, a sufficiently small grid size (<500 m), or subgrid interpolation of the grounding line. The latter still requires nominal grid sizes of <5 km. For larger grid spacings, appropriate parameterizations for ice flux may be imposed at the grounding line, but the short-time transient behaviour is then incorrect and different from models that do not incorporate grounding-line parameterizations. The numerical error associated with predicting grounding-line motion can be reduced significantly below the errors associated with parameter ignorance and uncertainties in future scenarios.
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.
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.