Search

The Greenland Ice Sheet (GIS) is losing mass at a high rate1. Given the short-term nature of the observational record, it is difficult to assess the historical importance of this mass-loss trend. Unlike records of greenhouse gas concentrations and global temperature, in which observations have been merged with palaeoclimate datasets, there are no comparably long records for rates of GIS mass change. Here we reveal unprecedented mass loss from the GIS this century, by placing contemporary and future rates of GIS mass loss within the context of the natural variability over the past 12,000 years. We force a high-resolution ice-sheet model with an ensemble of climate histories constrained by ice-core data2. Our simulation domain covers southwestern Greenland, the mass change of which is dominated by surface mass balance. The results agree favourably with an independent chronology of the history of the GIS margin3,4. The largest pre-industrial rates of mass loss (up to 6,000 billion tonnes per century) occurred in the early Holocene, and were similar to the contemporary (ad 2000–2018) rate of around 6,100 billion tonnes per century5. Simulations of future mass loss from southwestern GIS, based on Representative Concentration Pathway (RCP) scenarios corresponding to low (RCP2.6) and high (RCP8.5) greenhouse gas concentration trajectories6, predict mass loss of between 8,800 and 35,900 billion tonnes over the twenty-first century. These rates of GIS mass loss exceed the maximum rates over the past 12,000 years. Because rates of mass loss from the southwestern GIS scale linearly5 with the GIS as a whole, our results indicate, with high confidence, that the rate of mass loss from the GIS will exceed Holocene rates this century. Rates of ice-mass loss from southwestern Greenland this century will exceed the maximum rate over the past 12,000 years, and would not be the result of natural variation.

The land ice contribution to global mean sea level rise has not yet been predicted 1 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 9 and climate models 10, 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.

Projections of the sea level contribution from the Greenland and Antarctic ice sheets (GrIS and AIS) rely on atmospheric and oceanic drivers obtained from climate models. The Earth System Models participating in the Coupled Model Intercomparison Project phase 6 (CMIP6) generally project greater future warming compared with the previous Coupled Model Intercomparison Project phase 5 (CMIP5) effort. Here we use four CMIP6 models and a selection of CMIP5 models to force multiple ice sheet models as part of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). We find that the projected sea level contribution at 2100 from the ice sheet model ensemble under the CMIP6 scenarios falls within the CMIP5 range for the Antarctic ice sheet but is significantly increased for Greenland. Warmer atmosphere in CMIP6 models results in higher Greenland mass loss due to surface melt. For Antarctica, CMIP6 forcing is similar to CMIP5 and mass gain from increased snowfall counteracts increased loss due to ocean warming.

Paleoclimate proxies are being used in conjunction with ice sheet modeling experiments to determine how the Greenland ice sheet responded to past changes, particularly during the last deglaciation. Although these comparisons have been a critical component in our understanding of the Greenland ice sheet sensitivity to past warming, they often rely on modeling experiments that favor minimizing computational expense over increased model physics. Over Paleoclimate timescales, simulating the thermal structure of the ice sheet has large implications on the modeled ice viscosity, which can feedback onto the basal sliding and ice flow. To accurately capture the thermal field, models often require a high number of vertical layers. This is not the case for the stress balance computation, however, where a high vertical resolution is not necessary. Consequently, since stress balance and thermal equations are generally performed on the same mesh, more time is spent on the stress balance computation than is otherwise necessary. For these reasons, running a higher-order ice sheet model (e.g., Blatter-Pattyn) over timescales equivalent to the paleoclimate record has not been possible without incurring a large computational expense. To mitigate this issue, we propose a method that can be implemented within ice sheet models, whereby the vertical interpolation along the z axis relies on higher-order polynomials, rather than the traditional linear interpolation. This method is tested within the Ice Sheet System Model (ISSM) using quadratic and cubic finite elements for the vertical interpolation on an idealized case and a realistic Greenland configuration. A transient experiment for the ice thickness evolution of a single-dome ice sheet demonstrates improved accuracy using the higher-order vertical interpolation compared to models using the linear vertical interpolation, despite having fewer degrees of freedom. This method is also shown to improve a model's ability to capture sharp thermal gradients in an ice sheet particularly close to the bed, when compared to models using a linear vertical interpolation. This is corroborated in a thermal steady-state simulation of the Greenland ice sheet using a higher-order model. In general, we find that using a higher-order vertical interpolation decreases the need for a high number of vertical layers, while dramatically reducing model runtime for transient simulations. Results indicate that when using a higher-order vertical interpolation, runtimes for a transient ice sheet relaxation are upwards of 5 to 7 times faster than using a model which has a linear vertical interpolation, and this thus requires a higher number of vertical layers to achieve a similar result in simulated ice volume, basal temperature, and ice divide thickness. The findings suggest that this method will allow higher-order models to be used in studies investigating ice sheet behavior over paleoclimate timescales at a fraction of the computational cost than would otherwise be needed for a model using a linear vertical interpolation.

We report 80 10Be ages on 14 moraines from Irish cirques that show a previously unrecognized signal of at least eight millennial-scale fluctuations between 24.5±0.7 ka and 11.0±0.3 ka. Several moraine ages may be correlative with abrupt warming at the onset of the Bolling (14.7 ka) and the end of the Younger Dryas (11.7 ka), suggesting a forced response. Our ages also identify glacier fluctuations that occurred when regional temperatures were relatively stable. This finding is consistent with modeling results showing several hundred-meter-scale glacier fluctuations in response to interannual variability. At the same time, our composite record of cirque-glacier average equilibrium line altitudes (ELAs) shows a response to warming due to increasing greenhouse gases and summer insolation modulated by abrupt climate changes. Our new 10Be chronology thus records both forced and unforced millennial-scale glacier fluctuations superimposed on a lower frequency ELA signal of forced response to climate change.

Geologic archives constraining the variability of the Greenland ice sheet (GrIS) during the Holocene provide targets for ice sheet models to test sensitivities to variations in past climate and model formulation. Even as data–model comparisons are becoming more common, many models simulating the behavior of the GrIS during the past rely on meshes with coarse horizontal resolutions (≥10 km). In this study, we explore the impact of model resolution on the simulated nature of retreat across southwestern Greenland during the Holocene. Four simulations are performed using the Ice Sheet System Model (ISSM): three that use a uniform mesh and horizontal mesh resolutions of 20, 10, and 5 km, and one that uses a nonuniform mesh with a resolution ranging from 2 to 15 km. We find that the simulated retreat can vary significantly between models with different horizontal resolutions based on how well the bed topography is resolved. In areas of low topographic relief, the horizontal resolution plays a negligible role in simulated differences in retreat, with each model instead responding similarly to retreat driven by surface mass balance (SMB). Conversely, in areas where the bed topography is complex and high in relief, such as fjords, the lower-resolution models (10 and 20 km) simulate unrealistic retreat that occurs as ice surface lowering intersects bumps in the bed topography that would otherwise be resolved as troughs using the higher-resolution grids. Our results highlight the important role that high-resolution grids play in simulating retreat in areas of complex bed topography, but also suggest that models using nonuniform grids can save computational resources through coarsening the mesh in areas of noncomplex bed topography where the SMB predominantly drives retreat. Additionally, these results emphasize that care must be taken with ice sheet models when tuning model parameters to match reconstructed margins, particularly for lower-resolution models in regions where complex bed topography is poorly resolved.

Despite elevated summer insolation forcing during the early Holocene, global ice sheets retained nearly half of their volume from the Last Glacial Maximum, as indicated by deglacial records of global mean sea level (GMSL). Partitioning the GMSL rise among potential sources requires accurate dating of ice-sheet extent to estimate ice-sheet volume. Here, we date the final retreat of the Laurentide Ice Sheet with 10 Be surface exposure ages for the Labrador Dome, the largest of the remnant Laurentide ice domes during the Holocene. We show that the Labrador Dome deposited moraines during North Atlantic cold events at ∼10.3 ka, 9.3 ka and 8.2 ka, suggesting that these regional climate events helped stabilize the retreating Labrador Dome in the early Holocene. After Hudson Bay became seasonally ice free at ∼8.2 ka, the majority of Laurentide ice-sheet melted abruptly within a few centuries. We demonstrate through high-resolution regional climate model simulations that the thermal properties of a seasonally ice-free Hudson Bay would have increased Laurentide ice-sheet ablation and thus contributed to the subsequent rapid Labrador Dome retreat. Finally, our new 10 Be chronology indicates full Laurentide ice-sheet had completely deglaciated by 6.7 ± 0.4 ka, which re quires that Antarctic ice sheets contributed 3.6–6.5 m to GMSL rise since 6.3–7.1 ka.

Understanding patterns of ice-sheet deglaciation is key for predicting the rate of future ice-sheet melt, yet the processes underlying deglaciation remain elusive. The early Holocene (11.7 ka to 8.2 ka; Greenlandian Stage) represents the most recent period when the Laurentide and Greenland ice sheets underwent large-scale recession. Moreover, this ice-sheet recession occurred under the backdrop of regional temperatures that were similar to or warmer than today, and comparable to those projected for the upcoming centuries. Reconstructing Laurentide and Greenland ice sheet behavior during the early Holocene, and elucidating the mechanisms dictating this behavior may serve as a partial analog for future Greenland ice-sheet change in a warming world. Here, we present 123 new 10Be surface exposure ages from two sites on Baffin Island and southwestern Greenland that constrain the behavior of the Laurentide and Greenland ice sheets, and an independent alpine glacier during the early Holocene. On Baffin Island, we focus on a unique area where moraines deposited by the Laurentide Ice Sheet rest directly adjacent to moraines deposited by an independent alpine glacier. Sixty-one 10Be ages reveal that advances and/or stillstands of the Laurentide Ice Sheet and an alpine glacier occurred in unison around 11.8 ka, 10.3 ka, and 9.2 ka. Sixty-two 10Be ages from southwestern Greenland indicate that the Greenland Ice Sheet margin experienced re-advances or stillstands around 11.6 ka, 10.4 ka, 9.1 ka, 8.1 ka, and 7.3 ka. Our results reveal that ice sheets respond to climate perturbations on the same centennial timescale as small alpine glaciers. We hypothesize that during the warming climate of the early Holocene, freshening of the North Atlantic Ocean induced by ice-sheet melt resulted in regional cooling and brief periods of ice-sheet stabilization superimposed on net glacier recession. These observations point to a negative feedback mechanism inherent to melting ice sheets in the Baffin Bay region that slows ice-sheet recession during intervals of otherwise rapid deglaciation.

The last deglaciation of the Scandinavian Ice Sheet (SIS) from ∼ 21 , 000 to 13,000 yr ago is well-constrained by several hundred 10Be and 14C ages. The subsequent retreat history, however, is established primarily from minimum-limiting 14C ages and incomplete Baltic-Sea varve records, leaving a substantial fraction of final SIS retreat history poorly constrained. Here we develop a high-resolution chronology for the final deglaciation of the SIS based on 79 10Be cosmogenic exposure dates sampled along three transects spanning southern to northern Sweden and Finland. Combining this new chronology with existing 10Be ages on deglaciation since the Last Glacial Maximum shows that rates of SIS margin retreat were strongly influenced by deglacial millennial-scale climate variability and its effect on surface mass balance, with regional modulation of retreat associated with dynamical controls. Ice-volume estimates constrained by our new chronology suggest that the SIS contributed ∼ 8 m sea-level equivalent to global sea-level rise between ∼14.5 ka and 10 ka. Final deglaciation was largely complete by ∼10.5 ka, with highest rates of sea-level rise occurring during the Bolling–Allerod, a 50% decrease during the Younger Dryas, and a rapid increase during the early Holocene. Combining our SIS volume estimates with estimated contributions from other remaining Northern Hemisphere ice sheets suggests that the Antarctic Ice Sheet (AIS) contributed 14.4 ± 5.9 m to global sea-level rise since ∼13 ka. This new constraint supports those studies that indicate that an ice volume of 15 m or more of equivalent sea-level rise was lost from the AIS during the last deglaciation.

This study uses reanalysis data from ECMWF ERA-Interim and GCM output from the CCSM3 to investigate how sea ice and clouds interact locally (within individual grid boxes) and whether similar variability between the two datasets is captured. During autumn (October), the vertically integrated low cloud amount increases over increased sea ice in the reanalysis, but decreases in the GCM output. Closer inspection, however, reveals that both datasets have more low cloud cover over increased sea ice within the lower boundary layer (1000–925 hPa for the reanalysis and 1000–975 hPa for the GCM output), but they differ in their integrated response within the lower troposphere. These results highlight the differences between the datasets and show the importance of understanding where cloud changes occur, because clouds vary in their effect on the radiation budget as a function of height.