Search

This paper provides the first description of the open-source Glacier Energy and Mass Balance model. GEMB models the ice sheet and glacier surface–atmospheric energy and mass exchange, as well as the firn state. It is a column model (no horizontal communication) of intermediate complexity that includes those processes deemed most relevant to glacier studies. GEMB prioritizes computational efficiency to accommodate the very long (thousands of years) spin-ups necessary for initializing deep firn columns and sensitivity experiments needed to characterize model uncertainty on continental scales. The model is one-way coupled with the atmosphere, which allows the model to be run offline with a diversity of climate forcing but neglects feedback to the atmosphere. GEMB provides numerous parameterization choices for various key processes (e.g., albedo, subsurface shortwave absorption, and compaction), making it well suited for uncertainty quantification and model exploration. The model is evaluated against the current state of the art and in situ observations and is shown to perform well.

International audience; Ice losses from the Greenland and Antarctic ice sheets have accelerated since the 1990s, accounting for a significant increase in the global mean sea level. Here, we present a new 29-year record of ice sheet mass balance from 1992 to 2020 from the Ice Sheet Mass Balance Inter-comparison Exercise (IMBIE). We compare and combine 50 independent estimates of ice sheet mass balance derived from satellite observations of temporal changes in ice sheet flow, in ice sheet volume, and in Earth's gravity field. Between 1992 and 2020, the ice sheets contributed 21.0±1.9 mm to global mean sea level, with the rate of mass loss rising from 105 Gt yr−1 between 1992 and 1996 to 372 Gt yr−1 between 2016 and 2020. In Greenland, the rate of mass loss is 169±9 Gt yr−1 between 1992 and 2020, but there are large inter-annual variations in mass balance, with mass loss ranging from 86 Gt yr−1 in 2017 to 444 Gt yr−1 in 2019 due to large variability in surface mass balance. In Antarctica, ice losses continue to be dominated by mass loss from West Antarctica (82±9 Gt yr−1) and, to a lesser extent, from the Antarctic Peninsula (13±5 Gt yr−1). East Antarctica remains close to a state of balance, with a small gain of 3±15 Gt yr−1, but is the most uncertain component of Antarctica's mass balance. The dataset is publicly available at https://doi.org/10.5285/77B64C55-7166-4A06-9DEF-2E400398E452 (IMBIE Team, 2021).

Numerical simulations of the Greenland Ice Sheet (GrIS) over geologic timescales can greatly improve our knowledge of the critical factors driving GrIS demise during climatically warm periods, which has clear relevance for better predicting GrIS behavior over the upcoming centuries. To assess the fidelity of these modeling efforts, however, observational constraints of past ice sheet change are needed. Across southwestern Greenland, geologic records detail Holocene ice retreat across both terrestrial-based and marine-terminating environments, providing an ideal opportunity to rigorously benchmark model simulations against geologic reconstructions of ice sheet change. Here, we present regional ice sheet modeling results using the Ice-sheet and Sea-level System Model (ISSM) of Holocene ice sheet history across an extensive fjord region in southwestern Greenland covering the landscape around the Kangiata Nunaata Sermia (KNS) glacier and extending outward along the 200 km Nuup Kangerula (Godthåbsfjord). Our simulations, forced by reconstructions of Holocene climate and recently implemented calving laws, assess the sensitivity of ice retreat across the KNS region to atmospheric and oceanic forcing. Our simulations reveal that the geologically reconstructed ice retreat across the terrestrial landscape in the study area was likely driven by fluctuations in surface mass balance in response to Early Holocene warming – and was likely not influenced significantly by the response of adjacent outlet glaciers to calving and ocean-induced melting. The impact of ice calving within fjords, however, plays a significant role by enhancing ice discharge at the terminus, leading to interior thinning up to the ice divide that is consistent with reconstructed magnitudes of Early Holocene ice thinning. Our results, benchmarked against geologic constraints of past ice-margin change, suggest that while calving did not strongly influence Holocene ice-margin migration across terrestrial portions of the KNS forefield, it strongly impacted regional mass loss. While these results imply that the implementation and resolution of ice calving in paleo-ice-flow models is important towards making more robust estimations of past ice mass change, they also illustrate the importance these processes have on contemporary and future long-term ice mass change across similar fjord-dominated regions of the GrIS.

Ice sheet models are often initialized with data assimilation of present-day conditions, in which unknown model parameters are estimated using the inverse method. While assimilation of snapshot observations has been widely used for regional and large scale ice sheet simulations, data assimilation based on time dependent data has recently started to emerge to constrain model parameters while capturing the transient evolution of the system. However, this method has been applied only to a few glaciers with fixed ice front positions, using spatially and temporally limited observations, and has not been applied to marine terminating glaciers of the Greenland ice sheet that have been retreating over the last 30 years. In this study, we assimilate time series of surface velocity into a model of Kjer glacier in West Greenland to better capture the observed acceleration over the past three decades. We compare snapshot and transient inverse methods and investigate the impact of initialization procedures on the parameters inferred, as well as model projections. We find that transient-calibrated simulations better capture past trends and better reproduce changes after the calibration period, even when a short period of observations is used. The results show the feasibility and clear benefits of a time-dependent data assimilation for initializing ice sheet models. This approach is now possible with the development of longer observational records, though it remains computationally challenging.

Determining the future evolution of the Antarctic Ice Sheet is critical for understanding and narrowing the large existing uncertainties in century-scale global mean sea-level-rise (SLR) projections. One of the most significant glaciers and ice streams in Antarctica, Thwaites Glacier, is at risk of destabilization and, if destabilized, has the potential to be the largest regional-scale contributor of SLR on Earth. This is because Thwaites Glacier is vulnerable to the marine ice-sheet instability as its grounding line is significantly influenced by ocean-driven basal melting rates, and its bedrock topography retrogrades into kilometer-deep troughs. In this study, we investigate how bedrock topography features influence the grounding line migration beneath Thwaites Glacier when extreme ocean-driven basal melt rates are applied. Specifically, we design experiments using the Ice-sheet and Sea-level System Model (ISSM) to quantify the SLR projection uncertainty due to reported errors in the current bedrock topography maps that are often used by ice-sheet models. We find that spread in model estimates of sea-level-rise contribution from Thwaites Glacier due to the reported bedrock topography error could be as large as 21.9 cm after 200 years of extreme ocean warming. Next, we perturb the bedrock topography beneath Thwaites Glacier using wavelet decomposition techniques to introduce realistic noise (within error). We explore the model space with multiple realizations of noise to quantify what spatial and vertical resolutions in bedrock topography are required to minimize the uncertainty in our 200-year experiment. We conclude that at least a 2 km spatial and 8 m vertical resolution would independently constrain possible SLR to ±2 cm over 200 years, fulfilling requirements outlined by the 2017 Decadal Survey for Earth Science. Lastly, we perform an ensemble of simulations to determine in which regions our model of Thwaites Glacier is most sensitive to perturbations in bedrock topography. Our results suggest that the retreat of the grounding line is most sensitive to bedrock topography in proximity to the grounding line's initial position. Additionally, we find that the location and amplitude of the bedrock perturbation is more significant than its sharpness and shape. Overall, these findings inform and benchmark observational requirements for future missions that will measure ice-sheet bedrock topography, not only in the case of Thwaites Glacier but for Antarctica on the continental scale.

The evolution of the Antarctic Ice Sheet (AIS) represents one of the most important and uncertain contributions to sea level rise in the upcoming centuries. Thwaites Glacier and the Amundsen Sea sector of the West Antarctic Ice Sheet (WAIS) have been identified as the continent's most critical areas. The retreat of Thwaites Glacier's grounding line - the transition area where ice is no longer grounded and becomes afloat - is the subject of considerable study for modelers as it governs the collapse of the glacier. Recent advances towards the coupling of dynamical ice models with Glacial Isostatic Adjustment (GIA) models have provided the means to improve grounding line projections by considering solid-Earth processes and their interactions with the cryosphere and hydrosphere. However, the spatial and temporal model resolutions necessary to fully capture these interactions, and the sensitivity to model parametrization, remain elusive. We investigate the grounding line retreat of Thwaites Glacier through 2300 using the parallelized coupled physics capabilities of the Ice-sheet and Sea-level System Model (ISSM) which capture the complex interactions between solid-Earth, ice-sheets, and ocean. We incorporate realistic climatology, ocean melt rates, and GIA models and we discuss the impact of spatial and temporal model resolution, and solid-Earth parametrization, on the grounding line retreat and sea level change.

Ice losses from the Greenland and Antarctic Ice Sheets have accelerated since the 1990s, accounting for a significant increase in global mean sea level. Here, we present a new 29-year record of ice sheet mass balance from 1992 to 2020 from the Ice Sheet Mass Balance Inter-comparison Exercise (IMBIE). We compare and combine 50 independent estimates of ice sheet mass balance derived from satellite observations of temporal changes in ice sheet flow, in ice sheet volume and in Earth’s gravity field. Between 1992 and 2020, the ice sheets contributed 21.0 ± 1.9 mm to global mean sea-level, with the rate of mass loss rising from 105 Gt yr-1 between 1992 and 1996 to 372 Gt yr-1 between 2016 and 2020. In Greenland, the rate of mass loss is 169 ± 9 Gt yr-1 between 1992 and 2020 but there are large inter-annual variations in mass balance with mass loss ranging from 86 Gt yr-1 in 2017 to 444 Gt yr-1 in 2019 due to large variability in surface mass balance. In Antarctica, ice losses continue to be dominated by mass loss from West Antarctica (-82 ± 9 Gt yr-1) and to a lesser extent from the Antarctic Peninsula (-13 ± 5 Gt yr-1). East Antarctica remains close to a state of balance (3 ± 15 Gt yr-1), but is the most uncertain component of Antarctica’s mass balance.

This paper provides the first description of version 1.0 of the open-source Glacier Energy and Mass Balance model. GEMB models the ice sheet and glacier surface-atmospheric energy and mass exchange, and firn state. It is a column model (no horizontal communication) of intermediate complexity that includes those processes deemed most relevant to glacier studies with the goal of retaining computational efficiency that can accommodate the very long (thousands of years) spin-ups necessary for initializing deep firn columns and for running sensitivity experiments to characterize model uncertainty on continental scales. The model is one-way coupled with the atmosphere which allows the model to be run off-line with a diversity of climate forcing but neglects feedback to the atmosphere. GEMB provides numerous parameterization choices for various key processes (e.g. albedo, subsurface shortwave absorption, and compaction), making it well suited for uncertainty quantification and model exploration. The model is evaluated against the current state-of-the-art and in situ observations and is shown to perform well.

Ice shelves tend to grow through a steady influx of glacial ice and retreat in discrete calving events that occur on subannual to multidecadal timescales. The impacts of ice shelf calving and retreat are far-reaching, but the evolution of Antarctica’s coastline has not been well characterized, owing to the difficulty of delineating ice fronts in limited satellite data. To create an annual coastline dataset that spans the past quarter century, we combine data from multiple satellite sensors, and we use the known physics of ice flow to constrain ice front positions and fill gaps in the data record. We find that since 1997, Antarctica’s coastlines have retreated by 37,000 km2, led by major calving events from the Ross and Ronne ice shelves in the early 2000s, and sustained by countless loss events from smaller ice shelves ever since. Calving losses total nearly 6000 Gt, which is roughly equivalent to the total mass that has been lost to ice shelf thinning over the same period. Using an ice sheet model, we examine the impacts of observed coastal changes on the buttressing strength of Antarctica’s ice shelves.

All LIDAR instruments are not the same, and advancement of LIDAR technology requires an ongoing interest and demand from the community to foster further development of the required components. The ...

The Greenland Ice Sheet (GIS) is losing mass at a high rate

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 that arises from large uncertainty in the external forcing that future warming may exert onto the ice sheet. While ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean, the approach is neglecting a number of processes such as surface mass balance related contributions and mechanisms. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any dampening or self-amplifying processes. This is particularly relevant in situations where an instability is dominating the ice loss. 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 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 CMIP-5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP-5 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 is 9.6 mm with a likely range between 5.2 mm and 20.3 mm compared to the observed sea-level contribution from Antarctica of 7.4 mm with a standard deviation of 3.7 mm (Shepherd et al., 2018). For the so-called business-as-usual warming path, RCP-8.5, we obtain a median contribution of the Antarctic ice sheet to global mean sea-level rise within the 21st century of 17 cm with a likely range (66-percentile around the mean) between 9 cm and 36 cm and a very likely range (90-percentile around the mean) between 6 cm and 59 cm. For the RCP-2.6 warming path which will keep the global mean temperature below two degrees of global warming and is thus consistent with the Paris Climate Agreement yields a median of 13 cm of global mean sea-level contribution. The likely range for the RCP-2.6 scenario is between 7 cm and 25 cm and the very likely range is between 5 cm and 39 cm. The structural uncertainties in the method do not allow 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 RCP-8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade with a likely range between 2 cm/dec and 8 cm/dec and a very likely range between 1 cm/dec and 13 cm/dec.

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., SCOPUS: ar.j, info:eu-repo/semantics/published