Design and results of the ice sheet model initialisation experiments initMIP-Greenland: an ISMIP6 intercomparison

Holocene thinning in central Greenland controlled by the Northeast Greenland Ice Stream

Ice-core records from the interior of the Greenland ice sheet suggest widespread thinning during the Holocene. However, the recurring underestimation of this thinning in numerical models raises concerns about both the veracity of such reconstructions and the reliability of glaciological models. Recent work suggests the 8000-year-old Northeast Greenland Ice Stream (NEGIS), including a now-extinct northern tributary, may have been an early influence on Greenland ice-sheet dynamics. Yet, the inaccurate reproduction of NEGIS-like dynamics in most models hampers investigation of whether this feature played a role in Holocene ice-sheet thinning. Here we show that grounding-line retreat in northeast Greenland triggers elevation changes at the northern summit via ice-dynamic effects modulated by the paleo NEGIS system. In our simulations, fast ice-stream flow caused by transiently imposed reduced basal shear stress following the northeast retreat explains 55% ( ± 18%) of the estimated ice thinning, showing that ice-stream dynamics is one of the main drivers of the NGRIP Holocene surface elevation drop. Our findings show that the ice-flow in northeast Greenland plays a large role in ice-surface elevation changes in central Greenland.

Rapid disintegration and weakening of ice shelves in North Greenland

The glaciers of North Greenland are hosting enough ice to raise sea level by 2.1 m, and have long considered to be stable. This part of Greenland is buttressed by the last remaining ice shelves of the ice sheet. Here, we show that since 1978, ice shelves in North Greenland have lost more than 35% of their total volume, three of them collapsing completely. For the floating ice shelves that remain we observe a widespread increase in ice shelf mass losses, that are dominated by enhanced basal melting rates. Between 2000 and 2020, there was a widespread increase in basal melt rates that closely follows a rise in the ocean temperature. These glaciers are showing a direct dynamical response to ice shelf changes with retreating grounding lines and increased ice discharge. These results suggest that, under future projections of ocean thermal forcing, basal melting rates will continue to rise or remain at high level, which may have dramatic consequences for the stability of Greenlandic glaciers.

Climate model differences contribute deep uncertainty in future Antarctic ice loss

Future projections of ice sheets in response to different climate scenarios and their associated contributions to sea level changes are subject to deep uncertainty due to ice sheet instability processes, hampering a proper risk assessment of sea level rise and enaction of mitigation/adaptation strategies. For a systematic evaluation of the uncertainty due to climate model fields used as input to the ice sheet models, we drive a three-dimensional model of the Antarctic Ice Sheet (AIS) with the output from 36 climate models to simulate past and future changes in the AIS. Simulations show that a few climate models result in partial collapse of the West AIS under modeled preindustrial climates, and the spread in future changes in the AIS's volume is comparable to the structural uncertainty originating from differing ice sheet models. These results highlight the need for improved representations of physical processes important for polar climate in climate models.

The Transient Sea Level Response to External Forcing in CMIP6 Models

Earth is warming and sea levels are rising as land-based ice is lost to melt, and oceans expand due to accumulation of heat. The pace of ice loss and steric expansion is linked to the intensity of warming. How much faster sea level will rise as climate warms is, however, highly uncertain and difficult to model. Here, we quantify the transient sea level sensitivity of the sea level budget in both models and observations. Models show little change in sensitivity to warming between the first and second half of the twenty-first century for most contributors. The exception is glaciers and ice caps (GIC) that have a greater sensitivity pre-2050 (2.8 ± 0.4 mm/yr/K) compared to later (0.7 ± 0.1 mm/yr/K). We attribute this change to the short response time of glaciers and their changing area over time. Model sensitivities of steric expansion (1.5 ± 0.2 mm/yr/K), and Greenland Ice Sheet mass loss (0.8 ± 0.2 mm/yr/K) are greater than, but still compatible with, corresponding estimates from historical data (1.4 ± 0.5 and 0.4 ± 0.2 mm/yr/K). Antarctic Ice Sheet (AIS) models tends to show lower rates of sea level rise (SLR) with warming (-0.0 ± 0.3 mm/yr/K) in contrast to historical estimates (0.4 ± 0.2 mm/yr/K). This apparent low bias in AIS sensitivity is only partly able to account for a similar low bias identified in the sensitivity of global mean sea level excluding GIC (3.1 ± 0.4 vs. 2.3 ± 0.4 mm/yr/K). The balance temperature, where SLR is zero, lies close to the pre-industrial value, implying that SLR can only be mitigated by substantial global cooling.

Impact of paleoclimate on present and future evolution of the Greenland Ice Sheet

Using transient climate forcing based on simulations from the Alfred Wegener Institute Earth System Model (AWI-ESM), we simulate the evolution of the Greenland Ice Sheet (GrIS) from the last interglacial (125 ka, kiloyear before present) to 2100 AD with the Parallel Ice Sheet Model (PISM). The impact of paleoclimate, especially Holocene climate, on the present and future evolution of the GrIS is explored. Our simulations of the past show close agreement with reconstructions with respect to the recent timing of the peaks in ice volume and the climate of Greenland. The maximum and minimum ice volume at around 18-17 ka and 6-5 ka lag the respective extremes in climate by several thousand years, implying that the ice volume response of the GrIS strongly lags climatic changes. Given that Greenland's climate was getting colder from the Holocene Thermal Maximum (i.e., 8 ka) to the Pre-Industrial era, our simulation implies that the GrIS experienced growth from the mid-Holocene to the industrial era. Due to this background trend, the GrIS still gains mass until the second half of the 20th century, even though anthropogenic warming begins around 1850 AD. This is also in agreement with observational evidence showing mass loss of the GrIS does not begin earlier than the late 20th century. Our results highlight that the present evolution of the GrIS is not only controlled by the recent climate changes, but is also affected by paleoclimate, especially the relatively warm Holocene climate. We propose that the GrIS was not in equilibrium throughout the entire Holocene and that the slow response to Holocene climate needs to be represented in ice sheet simulations in order to predict ice mass loss, and therefore sea level rise, accurately.

Projected land ice contributions to twenty-first-century sea level rise

The land ice contribution to global mean sea level rise has not yet been predicted¹ using ice sheet and glacier models for the latest set of socio-economic scenarios, nor using coordinated exploration of uncertainties arising from the various computer models involved. Two recent international projects generated a large suite of projections using multiple models^2-8, but primarily used previous-generation scenarios⁹ and climate models¹⁰, and could not fully explore known uncertainties. Here we estimate probability distributions for these projections under the new scenarios^11,12 using statistical emulation of the ice sheet and glacier models. We find that limiting global warming to 1.5 degrees Celsius would halve the land ice contribution to twenty-first-century sea level rise, relative to current emissions pledges. The median decreases from 25 to 13 centimetres sea level equivalent (SLE) by 2100, with glaciers responsible for half the sea level contribution. The projected Antarctic contribution does not show a clear response to the emissions scenario, owing to uncertainties in the competing processes of increasing ice loss and snowfall accumulation in a warming climate. However, under risk-averse (pessimistic) assumptions, Antarctic ice loss could be five times higher, increasing the median land ice contribution to 42 centimetres SLE under current policies and pledges, with the 95th percentile projection exceeding half a metre even under 1.5 degrees Celsius warming. This would severely limit the possibility of mitigating future coastal flooding. Given this large range (between 13 centimetres SLE using the main projections under 1.5 degrees Celsius warming and 42 centimetres SLE using risk-averse projections under current pledges), adaptation planning for twenty-first-century sea level rise must account for a factor-of-three uncertainty in the land ice contribution until climate policies and the Antarctic response are further constrained.

Retreat of Humboldt Gletscher, North Greenland, Driven by Undercutting From a Warmer Ocean

Humboldt Gletscher is a 100-km wide, slow-moving glacier in north Greenland which holds a 19-cm global sea level equivalent. Humboldt has been the fourth largest contributor to sea level rise since 1972 but the cause of its mass loss has not been elucidated. Multi-beam echo sounding data collected in 2019 indicate a seabed 200 m deeper than previously known. Conductivity temperature depth data reveal the presence of warm water of Atlantic origin at 0°C at the glacier front and a warming of the ocean waters by 0.9 ± 0.1°C since 1962. Using an ocean model, we reconstruct grounded ice undercutting by the ocean, combine it with calculated retreat caused by ice thinning to floatation, and are able to fully explain the observed retreat. Two thirds of the retreat are caused by undercutting of grounded ice, which is a physical process not included in most ice sheet models.

Understanding of Contemporary Regional Sea-Level Change and the Implications for the Future

Global sea level provides an important indicator of the state of the warming climate, but changes in regional sea level are most relevant for coastal communities around the world. With improvements to the sea-level observing system, the knowledge of regional sea-level change has advanced dramatically in recent years. Satellite measurements coupled with in situ observations have allowed for comprehensive study and improved understanding of the diverse set of drivers that lead to variations in sea level in space and time. Despite the advances, gaps in the understanding of contemporary sea-level change remain and inhibit the ability to predict how the relevant processes may lead to future change. These gaps arise in part due to the complexity of the linkages between the drivers of sea-level change. Here we review the individual processes which lead to sea-level change and then describe how they combine and vary regionally. The intent of the paper is to provide an overview of the current state of understanding of the processes that cause regional sea-level change and to identify and discuss limitations and uncertainty in our understanding of these processes. Areas where the lack of understanding or gaps in knowledge inhibit the ability to provide the needed information for comprehensive planning efforts are of particular focus. Finally, a goal of this paper is to highlight the role of the expanded sea-level observation network-particularly as related to satellite observations-in the improved scientific understanding of the contributors to regional sea-level change.

An Efficient Ice Sheet/Earth System Model Spin-up Procedure for CESM2-CISM2: Description, Evaluation, and Broader Applicability

Spinning up a highly complex, coupled Earth system model (ESM) is a time consuming and computationally demanding exercise. For models with interactive ice sheet components, this becomes a major challenge, as ice sheets are sensitive to bidirectional feedback processes and equilibrate over glacial timescales of up to many millennia. This work describes and demonstrates a computationally tractable, iterative procedure for spinning up a contemporary, highly complex ESM that includes an interactive ice sheet component. The procedure alternates between a computationally expensive coupled configuration and a computationally cheaper configuration where the atmospheric component is replaced by a data model. By periodically regenerating atmospheric forcing consistent with the coupled system, the data atmosphere remains adequately constrained to ensure that the broader model state evolves realistically. The applicability of the method is demonstrated by spinning up the preindustrial climate in the Community Earth System Model Version 2 (CESM2), coupled to the Community Ice Sheet Model Version 2 (CISM2) over Greenland. The equilibrium climate state is similar to the control climate from a coupled simulation with a prescribed Greenland ice sheet, indicating that the iterative procedure is consistent with a traditional spin-up approach without interactive ice sheets. These results suggest that the iterative method presented here provides a faster and computationally cheaper method for spinning up a highly complex ESM, with or without interactive ice sheet components. The method described here has been used to develop the climate/ice sheet initial conditions for transient, ice sheet-enabled simulations with CESM2-CISM2 in the Coupled Model Intercomparison Project Phase 6 (CMIP6).

Data Reduction Using Statistical and Regression Approaches for Ice Velocity Derived by Landsat-8, Sentinel-1 and Sentinel-2

During the last decade, the number of available satellite observations has increased significantly, allowing for far more frequent measurements of the glacier speed. Appropriate methods of post-processing need to be developed to efficiently deal with the large volumes of data generated and relatively large intrinsic errors associated with the measurements. Here, we process and combine together measurements of ice velocity of Russell Gletscher in Greenland from three satellites—Sentinel-1, Sentinel-2, and Landsat-8, creating a multi-year velocity database with high temporal and spatial resolution. We then investigate post-processing methodologies with the aim of generating corrected, ordered, and simplified time series. We tested rolling mean and median, cubic spline regression, and linear non-parametric local regression (LOWESS) smoothing algorithms to reduce data noise, evaluated the results against ground-based GPS in one location, and compared the results between two locations with different characteristics. We found that LOWESS provides the best solution for noisy measurements that are unevenly distributed in time. Using this methodology with these sensors, we can robustly derive time series with temporal resolution of 2–3 weeks and improve the accuracy on the ice velocity to about 10 m/yr, or a factor of three compared to the initial measurements. The presented methodology could be applied to the entire Greenland ice sheet with an aim of reconstructing comprehensive sub-seasonal ice flow dynamics and mass balance.

An Integrated View of Greenland Ice Sheet Mass Changes Based on Models and Satellite Observations

The Greenland ice sheet is a major contributor to sea level rise, adding on average 0.47 ± 0.23 mm year − 1 to global mean sea level between 1991 and 2015. The cryosphere as a whole has contributed around 45% of observed global sea level rise since 1993. Understanding the present-day state of the Greenland ice sheet is therefore vital for understanding the processes controlling the modern-day rates of sea level change and for making projections of sea level rise into the future. Here, we provide an overview of the current state of the mass budget of Greenland based on a diverse range of remote sensing observations to produce the essential climate variables (ECVs) of ice velocity, surface elevation change, grounding line location, calving front location, and gravimetric mass balance as well as numerical modelling that together build a consistent picture of a shrinking ice sheet. We also combine these observations with output from a regional climate model and from an ice sheet model to gain insight into existing biases in ice sheet dynamics and surface mass balance processes. Observations show surface lowering across virtually all regions of the ice sheet and at some locations up to −2.65 m year − 1 between 1995 and 2017 based on radar altimetry analysis. In addition, calving fronts at 28 study sites, representing a sample of typical glaciers, have retreated all around Greenland since the 1990s and in only two out of 28 study locations have they remained stable. During the same period, two of five floating ice shelves have collapsed while the locations of grounding lines at the remaining three floating ice shelves have remained stable over the observation period. In a detailed case study with a fracture model at Petermann glacier, we demonstrate the potential sensitivity of these floating ice shelves to future warming. GRACE gravimetrically-derived mass balance (GMB) data shows that overall Greenland has lost 255 ± 15 Gt year − 1 of ice over the period 2003 to 2016, consistent with that shown by IMBIE and a marked increase compared to a rate of loss of 83 ± 63 Gt year − 1 in the 1993–2003 period. Regional climate model and ice sheet model simulations show that surface mass processes dominate the Greenland ice sheet mass budget over most of the interior. However, in areas of high ice velocity there is a significant contribution to mass loss by ice dynamical processes. Marked differences between models and observations indicate that not all processes are captured accurately within models, indicating areas for future research.

The paradigm shift in Antarctic ice sheet modelling

The Antarctic ice sheet is one of the largest potential contributors to future sea level rise. Predicting its future behaviour using physically-based ice sheet models has been a bottleneck for the past decades, but major advances are ongoing.

Rising Oceans Guaranteed: Arctic Land Ice Loss and Sea Level Rise

PURPOSE OF REVIEW

This paper reviews sea level contributions from land ice across the Arctic, including Greenland. We summarize ice loss measurement methods, ice loss mechanisms, and recent observations and projections, and highlight research advances over the last 3-5 years and remaining scientific challenges.

RECENT FINDINGS

Mass loss across the Arctic began to accelerate during the late twentieth century, with projections of continued loss across all future greenhouse gas emission scenarios. Recent research has improved knowledge of ice hydrology and surface processes, influences of atmospheric and oceanic changes on land ice, and boundary conditions such as subglacial topography. New computer models can also more accurately simulate glacier and ice sheet evolution.

SUMMARY

Rapid Arctic ice loss is underway, and future ice loss and sea level rise are guaranteed. Research continues to better understand and model physical processes and to improve projections of ice loss rates, especially after 2050.