Search Results

Now showing 1 - 9 of 9
  • Item
    A simple stress-based cliff-calving law
    (Göttingen : Copernicus GmbH, 2019) Schlemm, T.; Levermann, A.
    Over large coastal regions in Greenland and Antarctica the ice sheet calves directly into the ocean. In contrast to ice-shelf calving, an increase in calving from grounded glaciers contributes directly to sea-level rise. Ice cliffs with a glacier freeboard larger than ≈100 m are currently not observed, but it has been shown that such ice cliffs are increasingly unstable with increasing ice thickness. This cliff calving can constitute a self-amplifying ice loss mechanism that may significantly alter sea-level projections both of Greenland and Antarctica. Here we seek to derive a minimalist stress-based parametrization for cliff calving from grounded glaciers whose freeboards exceed the 100 m stability limit derived in previous studies. This will be an extension of existing calving laws for tidewater glaciers to higher ice cliffs.

    To this end we compute the stress field for a glacier with a simplified two-dimensional geometry from the two-dimensional Stokes equation. First we assume a constant yield stress to derive the failure region at the glacier front from the stress field within the glacier. Secondly, we assume a constant response time of ice failure due to exceedance of the yield stress. With this strongly constraining but very simple set of assumptions we propose a cliff-calving law where the calving rate follows a power-law dependence on the freeboard of the ice with exponents between 2 and 3, depending on the relative water depth at the calving front. The critical freeboard below which the ice front is stable decreases with increasing relative water depth of the calving front. For a dry water front it is, for example, 75 m. The purpose of this study is not to provide a comprehensive calving law but to derive a particularly simple equation with a transparent and minimalist set of assumptions.

  • Item
    Scaling of instability timescales of Antarctic outlet glaciers based on one-dimensional similitude analysis
    (Göttingen : Copernicus GmbH, 2019) Levermann, A.; Feldmann, J.
    Recent observations and ice-dynamic modeling suggest that a marine ice-sheet instability (MISI) might have been triggered in West Antarctica. The corresponding outlet glaciers, Pine Island Glacier (PIG) and Thwaites Glacier (TG), showed significant retreat during at least the last 2 decades. While other regions in Antarctica have the topographic predisposition for the same kind of instability, it is so far unclear how fast these instabilities would unfold if they were initiated. Here we employ the concept of similitude to estimate the characteristic timescales of several potentially MISI-prone outlet glaciers around the Antarctic coast. Our results suggest that TG and PIG have the fastest response time of all investigated outlets, with TG responding about 1.25 to 2 times as fast as PIG, while other outlets around Antarctica would be up to 10 times slower if destabilized. These results have to be viewed in light of the strong assumptions made in their derivation. These include the absence of ice-shelf buttressing, the one-dimensionality of the approach and the uncertainty of the available data. We argue however that the current topographic situation and the physical conditions of the MISI-prone outlet glaciers carry the information of their respective timescale and that this information can be partially extracted through a similitude analysis.
  • Item
    Modeling Antarctic tides in response to ice shelf thinning and retreat
    (Hoboken, NJ : Blackwell Publishing Ltd, 2014) Rosier, S.H.R.; Green, J.A.M.; Scourse, J.D.; Winkelmann, R.
    Tides play an important role in ice sheet dynamics by modulating ice stream velocity, fracturing, and moving ice shelves and mixing water beneath them. Any changes in ice shelf extent or thickness will alter the tidal dynamics through modification of water column thickness and coastal topography but these will in turn feed back onto the overall ice shelf stability. Here, we show that removal or reduction in extent and/or thickness of the Ross and Ronne-Filchner ice shelves would have a significant impact on the tides around Antarctica. The Ronne-Filchner appears particularly vulnerable, with an increase in M2 amplitude of over 0.5 m beneath much of the ice shelf potentially leading to tidally induced feedbacks on ice shelf/sheet dynamics. These results highlight the importance of understanding tidal feedbacks on ice shelves/streams due to their influence on ice sheet dynamics.
  • Item
    Experimental design for three interrelated marine ice sheet and ocean model intercomparison projects: MISMIP v. 3 (MISMIP +), ISOMIP v. 2 (ISOMIP +) and MISOMIP v. 1 (MISOMIP1)
    (München : European Geopyhsical Union, 2016) Asay-Davis, Xylar S.; Cornford, Stephen L.; Durand, Gaël; Galton-Fenzi, Benjamin K.; Gladstone, Rupert M.; Gudmundsson, G. Hilmar; Hattermann, Tore; Holland, David M.; Holland, Denise; Holland, Paul R.; Martin, Daniel F.; Mathiot, Pierre; Pattyn, Frank; Seroussi, Hélène
    Coupled ice sheet–ocean models capable of simulating moving grounding lines are just becoming available. Such models have a broad range of potential applications in studying the dynamics of marine ice sheets and tidewater glaciers, from process studies to future projections of ice mass loss and sea level rise. The Marine Ice Sheet–Ocean Model Intercomparison Project (MISOMIP) is a community effort aimed at designing and coordinating a series of model intercomparison projects (MIPs) for model evaluation in idealized setups, model verification based on observations, and future projections for key regions of the West Antarctic Ice Sheet (WAIS). Here we describe computational experiments constituting three interrelated MIPs for marine ice sheet models and regional ocean circulation models incorporating ice shelf cavities. These consist of ice sheet experiments under the Marine Ice Sheet MIP third phase (MISMIP+), ocean experiments under the Ice Shelf-Ocean MIP second phase (ISOMIP+) and coupled ice sheet–ocean experiments under the MISOMIP first phase (MISOMIP1). All three MIPs use a shared domain with idealized bedrock topography and forcing, allowing the coupled simulations (MISOMIP1) to be compared directly to the individual component simulations (MISMIP+ and ISOMIP+). The experiments, which have qualitative similarities to Pine Island Glacier Ice Shelf and the adjacent region of the Amundsen Sea, are designed to explore the effects of changes in ocean conditions, specifically the temperature at depth, on basal melting and ice dynamics. In future work, differences between model results will form the basis for the evaluation of the participating models.
  • Item
    Parameterization for subgrid-scale motion of ice-shelf calving fronts
    (München : European Geopyhsical Union, 2011) Albrecht, T.; Martin, M.; Haseloff, M.; Winkelmann, R.; Levermann, A.
    In order to explore the response of the Greenland ice sheet (GIS) to climate change on long (centennial to multi-millennial) time scales, a regional energy-moisture balance model has been developed. This model simulates seasonal variations of temperature and precipitation over Greenland and explicitly accounts for elevation and albedo feedbacks. From these fields, the annual mean surface temperature and surface mass balance can be determined and used to force an ice sheet model. The melt component of the surface mass balance is computed here using both a positive degree day approach and a more physically-based alternative that includes insolation and albedo explicitly. As a validation of the climate model, we first simulated temperature and precipitation over Greenland for the prescribed, present-day topography. Our simulated climatology compares well to observations and does not differ significantly from that of a simple parameterization used in many previous simulations. Furthermore, the calculated surface mass balance using both melt schemes falls within the range of recent regional climate model results. For a prescribed, ice-free state, the differences in simulated climatology between the regional energy-moisture balance model and the simple parameterization become significant, with our model showing much stronger summer warming. When coupled to a three-dimensional ice sheet model and initialized with present-day conditions, the two melt schemes both allow realistic simulations of the present-day GIS.
  • Item
    Antarctic sub-shelf melt rates via PICO
    (München : European Geopyhsical Union, 2018) Reese, Ronja; Albrecht, Torsten; Mengel, Matthias; Asay-Davis, Xylar; Winkelmann, Ricarda
    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 m a−1 for cold sub-shelf cavities, for example, underneath Ross or Ronne ice shelves, to 16 m a−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.
  • Item
    Interaction of marine ice-sheet instabilities in two drainage basins: Simple scaling of geometry and transition time
    (München : European Geopyhsical Union, 2015) Feldmann, J.; Levermann, A.
    The initiation of a marine ice-sheet instability (MISI) is generally discussed from the ocean side of the ice sheet. It has been shown that the reduction in ice-shelf buttressing and softening of the coastal ice can destabilize a marine ice sheet if the bedrock is sloping upward towards the ocean. Using a conceptional flow-line geometry, we investigate the possibility of whether a MISI can be triggered from the direction of the ice divide as opposed to coastal forcing and explore the interaction between connected basins. We find that the initiation of a MISI in one basin can induce a destabilization in the other. The underlying mechanism of basin interaction is based on dynamic thinning and a consecutive motion of the ice divide which induces a thinning in the adjacent basin and a successive initiation of the instability. Our simplified and symmetric topographic setup allows scaling both the geometry and the transition time between both instabilities. We find that the ice profile follows a universal shape that is scaled with the horizontal extent of the ice sheet and that the same exponent of 1/2 applies for the scaling relation between central surface elevation and horizontal extent as in the pure shallow ice approximation (Vialov profile). Altering the central bed elevation, we find that the extent of grounding-line retreat in one basin determines the degree of interaction with the other. Different scenarios of basin interaction are discussed based on our modeling results as well as on a conceptual flux-balance analysis. We conclude that for the three-dimensional case, the possibility of drainage basin interaction on timescales on the order of 1 kyr or larger cannot be excluded and hence needs further investigation.
  • Item
    Fracture-induced softening for large-scale ice dynamics
    (München : European Geopyhsical Union, 2014) Albrecht, T.; Levermann, A.
    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.
  • Item
    The Potsdam Parallel Ice Sheet Model (PISM-PIK) - Part 2: Dynamic equilibrium simulation of the Antarctic ice sheet
    (München : European Geopyhsical Union, 2011) Martin, M.A.; Winkelmann, R.; Haseloff, M.; Albrecht, T.; Bueler, E.; Khroulev, C.; Levermann, A.
    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.