Description: We present a revised and extended high Arctic air temperature reconstruction from a single proxy that spans the past ∼12,000 y (up to 2009 CE). Our reconstruction from the Agassiz ice cap (Ellesmere Island, Canada) indicates an earlier and warmer Holocene thermal maximum with early Holocene temperatures that are 4-5 °C warmer compared with a previous reconstruction, and regularly exceed contemporary values for a period of ∼3,000 y. Our results show that air temperatures in this region are now at their warmest in the past 6,800-7,800 y, and that the recent rate of temperature change is unprecedented over the entire Holocene. The warmer early Holocene inferred from the Agassiz ice core leads to an estimated ∼1 km of ice thinning in northwest Greenland during the early Holocene using the Camp Century ice core. Ice modeling results show that this large thinning is consistent with our air temperature reconstruction. The modeling results also demonstrate the broader significance of the enhanced warming, with a retreat of the northern ice margin behind its present position in the mid Holocene and a ∼25% increase in total Greenland ice sheet mass loss (∼1.4 m sea-level equivalent) during the last deglaciation, both of which have implications for interpreting geodetic measurements of land uplift and gravity changes in northern Greenland.

Description: We provide a global 0.5-degree grid of vertical land motion (in mm/a) of the LM17.3 glacial isostatic adjustment (GIA) model. The radially varying earth model part is profile VM5a (Peltier et al. 2015). The ice load is different to any other GIA model and combines regional ice loads without taking care of balancing the global sea-level equivalent of all ice sheets and glaciers with that expected from paleo-sea-level indicators. The regional models are: * GLAC-1D for North America (Tarasov et al. 2012), * HUY3 for Greenland (Lecavalier et al. 2014), * GLAC #71340 for Fennoscandia/Barents Sea (Tarasov et al., 2014), * ANU-ICE for Iceland, High Mountain Areas, Siberian Mountains and Tibet (Lambeck et al. 2014), * IJ04_Patagonia for Patagonia (updated from Ivins & James 2004), * ICE-6G_C for New Zealand (Argus et al. 2014, Peltier et al. 2015), * GLAC-1D for Antarctica (Briggs et al. 2014). Additional models (W12, Whitehouse et al. 2012, and IJ05_R2, Ivins et al. 2013, for Antarctica; ANU-ICE, Lambeck et al. 2017, and NAIce, Gowan et al. 2016, for North America) were tested in the development of the model but not used in the end. Little ice age is not included nor any ice mass change during the last 100 years. The eustatic sea-level equivalent at last glacial maximum amounts to 113.8 m for all ice sheets and glaciers together. Because we use an ice model that has not been tuned to fit global constraints, it may highlight areas which cannot match commonly used GIA observations. However, we note that the earth model used in our calculations is different to the earth model used in the development of some regional ice models, e.g. HUY3, ANU-ICE, IJ04_Patagonia (see respective references), thus some differences can be related to this. The LM17.3 model was introduced in Jäggi et al. (2019), and its DDK5-filtered geoid and water heights can be found in the EGSIEM plotter (http://plot.egsiem.eu/index.php?p=timeseries). The GIA model uses material compressibility and includes time-dependent coastlines and rotational feedback. The vertical land motion can be used/tested in sea-level investigations and projections. Work towards a model that incorporates 3D earth structure, and an updated ice model, is ongoing.

Description: The Antarctic Ice Sheet is an important indicator of climate change and driver of sea-level rise. Here we combine satellite observations of its changing volume, flow and gravitational attraction with modelling of its surface mass balance to show that it lost 2,720 ± 1,390 billion tonnes of ice between 1992 and 2017, which corresponds to an increase in mean sea level of 7.6 ± 3.9 millimetres (errors are one standard deviation). Over this period, ocean-driven melting has caused rates of ice loss from West Antarctica to increase from 53 ± 29 billion to 159 ± 26 billion tonnes per year; ice-shelf collapse has increased the rate of ice loss from the Antarctic Peninsula from 7 ± 13 billion to 33 ± 16 billion tonnes per year. We find large variations in and among model estimates of surface mass balance and glacial isostatic adjustment for East Antarctica, with its average rate of mass gain over the period 1992–2017 (5 ± 46 billion tonnes per year) being the least certain.

Description: The Greenland Ice Sheet has been a major contributor to global sea-level rise in recent decades1,2, and it is expected to continue to be so3. Although increases in glacier flow4–6 and surface melting7–9 have been driven by oceanic10–12 and atmospheric13,14 warming, the magnitude and trajectory of the ice sheet’s mass imbalance remain uncertain. Here we compare and combine 26 individual satellite measurements of changes in the ice sheet’s volume, flow and gravitational potential to produce a reconciled estimate of its mass balance. The ice sheet was close to a state of balance in the 1990s, but annual losses have risen since then, peaking at 345 ± 66 billion tonnes per year in 2011. In all, Greenland lost 3,902 ± 342 billion tonnes of ice between 1992 and 2018, causing the mean sea level to rise by 10.8 ± 0.9 millimetres. Using three regional climate models, we show that the reduced surface mass balance has driven 1,964 ± 565 billion tonnes (50.3 per cent) of the ice loss owing to increased meltwater runoff. The remaining 1,938 ± 541 billion tonnes (49.7 per cent) of ice loss was due to increased glacier dynamical imbalance, which rose from 46 ± 37 billion tonnes per year in the 1990s to 87 ± 25 billion tonnes per year since then. The total rate of ice loss slowed to 222 ± 30 billion tonnes per year between 2013 and 2017, on average, as atmospheric circulation favoured cooler conditions15 and ocean temperatures fell at the terminus of Jakobshavn Isbræ16. Cumulative ice losses from Greenland as a whole have been close to the rates predicted by the Intergovernmental Panel on Climate Change for their high-end climate warming scenario17, which forecast an additional 70 to 130 millimetres of global sea-level rise by 2100 compared with their central estimate.

Description: The Last Glacial Maximum (LGM, 21 000 years ago) is one of the suite of paleoclimate simulations included in the current phase of the Coupled Model Intercomparison Project (CMIP6). It is an interval when insolation was similar to the present, but global ice volume was at a maximum, eustatic sea level was at or close to a minimum, greenhouse gas concentrations were lower, atmospheric aerosol loadings were higher than today, and vegetation and land-surface characteristics were different from today. The LGM has been a focus for the Paleoclimate Modelling Intercomparison Project (PMIP) since its inception, and thus many of the problems that might be associated with simulating such a radically different climate are well documented. The LGM state provides an ideal case study for evaluating climate model performance because the changes in forcing and temperature between the LGM and pre-industrial are of the same order of magnitude as those projected for the end of the 21st century. Thus, the CMIP6 LGM experiment could provide additional information that can be used to constrain estimates of climate sensitivity. The design of the Tier 1 LGM experiment (lgm) includes an assessment of uncertainties in boundary conditions, in particular through the use of different reconstructions of the ice sheets and of the change in dust forcing. Additional (Tier 2) sensitivity experiments have been designed to quantify feedbacks associated with land-surface changes and aerosol loadings, and to isolate the role of individual forcings. Model analysis and evaluation will capitalize on the relative abundance of paleoenvironmental observations and quantitative climate reconstructions already available for the LGM.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Love, R., Andres, H. J., Condron, A., & Tarasov, L. Freshwater routing in eddy-permitting simulations of the last deglacial: the impact of realistic freshwater discharge. Climate of the Past, 17(6), (2021): 2327–2341. https://doi.org/10.5194/cp-17-2327-2021.

Description: Freshwater, in the form of glacial runoff, is hypothesized to play a critical role in centennial- to millennial-scale climate variability, such as the Younger Dryas and Dansgaard–Oeschger events, but this relationship is not straightforward. Large-scale glacial runoff events, such as Meltwater Pulse 1a (MWP1a), are not always temporally proximal to subsequent large-scale cooling. Moreover, the typical design of hosing experiments that support this relationship tends to artificially amplify the climate response. This study explores the impact that limitations in the representation of runoff in conventional “hosing” simulations has on our understanding of this relationship by examining where coastally released freshwater is transported when it reaches the ocean. We particularly focus on the impact of (1) the injection of freshwater directly over sites of deep-water formation (DWF) rather than at runoff locations (i.e. hosing), (2) excessive freshwater injection volumes (often by a factor of 5), and (3) the use of present-day (rather than palaeo) ocean gateways. We track the routing of glaciologically constrained freshwater volumes from four different inferred injection locations in a suite of eddy-permitting glacial ocean simulations using the Massachusetts Institute of Technology General Circulation Model (MITgcm) under both open and closed Bering Strait conditions. Restricting freshwater forcing values to realistic ranges results in less spreading of freshwater across the North Atlantic and indicates that the freshwater anomalies over DWF sites depend strongly on the geographical location of meltwater input. In particular, freshwater released into the Gulf of Mexico generates a very weak freshwater signal over DWF regions as a result of entrainment by the turbulent Gulf Stream. In contrast, freshwater released into the Arctic with an open Bering Strait or from the Eurasian ice sheet is found to generate the largest salinity anomalies over DWF regions in the North Atlantic and GIN (Greenland–Iceland–Norwegian) seas region respectively. Experiments show that when the Bering Strait is open, the Mackenzie River source exhibits more than twice as much freshening of the North Atlantic deep-water formation regions as when the Bering Strait is closed. Our results illustrate that applying freshwater hosing directly into the North Atlantic with even “realistic” freshwater amounts still overestimates the amount of terrestrial runoff reaching DWF regions. Given the simulated salinity anomaly distributions and the lack of reconstructed impact on deep-water formation during the Bølling–Allerød, our results support that the majority of the North American contribution to MWP1a was not routed through the Mackenzie River.

Description: This research has been supported by the Bundesministerium für Bildung und Forschung (Research for Sustainability initiative, FONA, through the PalMod project), the Natural Sciences and Engineering Research Council of Canada Discovery Grant, and by the NSF (grant no. OCE-1903427).

Description: The land ice contribution to global mean sea level rise has not yet been predicted1 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 models2,3,4,5,6,7,8, but primarily used previous-generation scenarios9 and climate models10, and could not fully explore known uncertainties. Here we estimate probability distributions for these projections under the new scenarios11,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.