Search Results
Simple models for the simulation of submarine melt for a Greenland glacial system model
Two hundred marine-terminating Greenland outlet glaciers deliver more than half of the annually accumulated ice into the ocean and have played an important role in the Greenland ice sheet mass loss observed since the mid-1990s. Submarine melt may play a crucial role in the mass balance and position of the grounding line of these outlet glaciers. As the ocean warms, it is expected that submarine melt will increase, potentially driving outlet glaciers retreat and contributing to sea level rise. Projections of the future contribution of outlet glaciers to sea level rise are hampered by the necessity to use models with extremely high resolution of the order of a few hundred meters. That requirement in not only demanded when modeling outlet glaciers as a stand alone model but also when coupling them with high-resolution 3-D ocean models. In addition, fjord bathymetry data are mostly missing or inaccurate (errors of several hundreds of meters), which questions the benefit of using computationally expensive 3-D models for future predictions. Here we propose an alternative approach built on the use of a computationally efficient simple model of submarine melt based on turbulent plume theory. We show that such a simple model is in reasonable agreement with several available modeling studies. We performed a suite of experiments to analyze sensitivity of these simple models to model parameters and climate characteristics. We found that the computationally cheap plume model demonstrates qualitatively similar behavior as 3-D general circulation models. To match results of the 3-D models in a quantitative manner, a scaling factor of the order of 1 is needed for the plume models. We applied this approach to model submarine melt for six representative Greenland glaciers and found that the application of a line plume can produce submarine melt compatible with observational data. Our results show that the line plume model is more appropriate than the cone plume model for simulating the average submarine melting of real glaciers in Greenland.
An efficient regional energy-moisture balance model for simulation of the Greenland Ice Sheet response to climate change
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.
Simulating the Greenland ice sheet under present-day and palaeo constraints including a new discharge parameterization
In this paper, we propose a new sub-grid scale parameterization for the ice discharge into the ocean through outlet glaciers and inspect the role of different observational and palaeo constraints for the choice of an optimal set of model parameters. This parameterization was introduced into the polythermal ice-sheet model SICOPOLIS, which is coupled to the regional climate model of intermediate complexity REMBO. Using the coupled model, we performed large ensemble simulations over the last two glacial cycles by varying two major parameters: a melt parameter in the surface melt scheme of REMBO and a discharge scaling parameter in our parameterization of ice discharge. Our empirical constraints are the present-day Greenland ice sheet surface elevation, the surface mass balance partition (ratio between total ice discharge and total precipitation) and the Eemian interglacial elevation drop relative to present day in the vicinity of the NEEM ice core. We show that the ice discharge parameterization enables us to simulate both the correct ice-sheet shape and mass balance partition at the same time without explicitly resolving the Greenland outlet glaciers. For model verification, we compare the simulated total and sectoral ice discharge with other estimates. For the model versions that are consistent with the range of observational and palaeo constraints, our simulated Greenland ice sheet contribution to Eemian sea-level rise relative to present-day amounts to 1.4 m on average (in the range of 0.6 and 2.5 m).
The Potsdam Parallel Ice Sheet Model (PISM-PIK) - Part 2: Dynamic equilibrium simulation of the Antarctic ice sheet
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.
SEMIC: An efficient surface energy and mass balance model applied to the Greenland ice sheet
We present SEMIC, a Surface Energy and Mass balance model of Intermediate Complexity for snow- and ice-covered surfaces such as the Greenland ice sheet. SEMIC is fast enough for glacial cycle applications, making it a suitable replacement for simpler methods such as the positive degree day (PDD) method often used in ice sheet modelling. Our model explicitly calculates the main processes involved in the surface energy and mass balance, while maintaining a simple interface and requiring minimal data input to drive it. In this novel approach, we parameterise diurnal temperature variations in order to more realistically capture the daily thaw–freeze cycles that characterise the ice sheet mass balance. We show how to derive optimal model parameters for SEMIC specifically to reproduce surface characteristics and day-to-day variations similar to the regional climate model MAR (Modèle Atmosphérique Régional, version 2) and its incorporated multilayer snowpack model SISVAT (Soil Ice Snow Vegetation Atmosphere Transfer). A validation test shows that SEMIC simulates future changes in surface temperature and surface mass balance in good agreement with the more sophisticated multilayer snowpack model SISVAT included in MAR. With this paper, we present a physically based surface model to the ice sheet modelling community that is general enough to be used with in situ observations, climate model, or reanalysis data, and that is at the same time computationally fast enough for long-term integrations, such as glacial cycles or future climate change scenarios.
Fracture-induced softening for large-scale ice dynamics
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.
Antarctic sub-shelf melt rates via PICO
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.
Parameterization for subgrid-scale motion of ice-shelf calving fronts
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.
From cyclic ice streaming to Heinrich-like events: the grow-and-surge instability in the parallel ice sheet model
Here we report on a cyclic, physical ice-discharge instability in the Parallel Ice Sheet Model, simulating the flow of a three-dimensional, inherently buttressed ice-sheet-shelf system which periodically surges on a millennial timescale. The thermomechanically coupled model on 1 km horizontal resolution includes an enthalpy-based formulation of the thermodynamics, a nonlinear stress-balance-based sliding law and a very simple subglacial hydrology. The simulated unforced surging is characterized by rapid ice streaming through a bed trough, resulting in abrupt discharge of ice across the grounding line which is eventually calved into the ocean. We visualize the central feedbacks that dominate the subsequent phases of ice buildup, surge and stabilization which emerge from the interaction between ice dynamics, thermodynamics and the subglacial till layer. Results from the variation of surface mass balance and basal roughness suggest that ice sheets of medium thickness may be more susceptible to surging than relatively thin or thick ones for which the surge feedback loop is damped. We also investigate the influence of different basal sliding laws (ranging from purely plastic to nonlinear to linear) on possible surging. The presented mechanisms underlying our simulations of self-maintained, periodic ice growth and destabilization may play a role in large-scale ice-sheet surging, such as the surging of the Laurentide Ice Sheet, which is associated with Heinrich events, and ice-stream shutdown and reactivation, such as observed in the Siple Coast region of West Antarctica.
Design and results of the ice sheet model initialisation experiments initMIP-Greenland: an ISMIP6 intercomparison
Earlier large-scale Greenland ice sheet sea-level projections (e.g. those run during the ice2sea and SeaRISE initiatives) have shown that ice sheet initial conditions have a large effect on the projections and give rise to important uncertainties. The goal of this initMIP-Greenland intercomparison exercise is to compare, evaluate, and improve the initialisation techniques used in the ice sheet modelling community and to estimate the associated uncertainties in modelled mass changes. initMIP-Greenland is the first in a series of ice sheet model intercomparison activities within ISMIP6 (the Ice Sheet Model Intercomparison Project for CMIP6), which is the primary activity within the Coupled Model Intercomparison Project Phase 6 (CMIP6) focusing on the ice sheets. Two experiments for the large-scale Greenland ice sheet have been designed to allow intercomparison between participating models of (1) the initial present-day state of the ice sheet and (2) the response in two idealised forward experiments. The forward experiments serve to evaluate the initialisation in terms of model drift (forward run without additional forcing) and in response to a large perturbation (prescribed surface mass balance anomaly); they should not be interpreted as sea-level projections. We present and discuss results that highlight the diversity of data sets, boundary conditions, and initialisation techniques used in the community to generate initial states of the Greenland ice sheet. We find good agreement across the ensemble for the dynamic response to surface mass balance changes in areas where the simulated ice sheets overlap but differences arising from the initial size of the ice sheet. The model drift in the control experiment is reduced for models that participated in earlier intercomparison exercises.