Description: table water isotopes in polar ice cores are widely used to reconstruct past temperature variations over several orbital climatic cycles. One way to calibrate the isotope– temperature relationship is to apply the present-day spatial relationship as a surrogate for the temporal one. However, this method leads to large uncertainties because several fac- tors like the sea surface conditions or the origin and transport of water vapor influence the isotope–temperature temporal slope. In this study, we investigate how the sea surface tem- perature (SST), the sea ice extent, and the strength of the Atlantic Meridional Overturning Circulation (AMOC) affect these temporal slopes in Greenland and Antarctica for Last Glacial Maximum (LGM, ∼ 21 000 years ago) to preindus- trial climate change. For that, we use the isotope-enabled at- mosphere climate model ECHAM6-wiso, forced with a set of sea surface boundary condition datasets based on recon- structions (e.g., GLOMAP) or MIROC 4m simulation out- puts. We found that the isotope–temperature temporal slopes in East Antarctic coastal areas are mainly controlled by the sea ice extent, while the sea surface temperature cooling af- fects the temporal slope values inland more. On the other hand, ECHAM6-wiso simulates the impact of sea ice extent on the EPICA Dome C (EDC) and Vostok sites through the contribution of water vapor from lower latitudes. Effects of sea surface boundary condition changes on modeled isotope-emperature temporal slopes are variable in West Antarctica. This is partly due to the transport of water vapor from the Southern Ocean to this area that can dampen the influence of local temperature on the changes in the isotopic compo- sition of precipitation and snow. In the Greenland area, the isotope–temperature temporal slopes are influenced by the sea surface temperatures near the coasts of the continent. The greater the LGM cooling off the coast of southeastern Green- land, the greater the transport of water vapor from the North Atlantic, and the larger the temporal slopes. The presence or absence of sea ice very near the coast has a large influ- ence in Baffin Bay and the Greenland Sea and influences the slopes at some inland ice core stations. The extent of the sea ice far south slightly influences the temporal slopes in Greenland through the transport of more depleted water vapor from lower latitudes to this area. The seasonal vari- ations of sea ice distribution, especially its retreat in sum- mer, influence the isotopic composition of the water vapor in this region and the modeled isotope–temperature tempo- ral slopes in the eastern part of Greenland. A stronger LGM AMOC decreases LGM-to-preindustrial isotopic anomalies in precipitation in Greenland, degrading the isotopic model– data agreement. The AMOC strength modifies the temporal slopes over inner Greenland slightly and by a little on the coasts along the Greenland Sea where the changes in surface temperature and sea ice distribution due to the AMOC strength mainly occur.

Description: The Antarctic Ice Sheet represents the largest source of uncertainty in future sea level rise projections, with a contribution to sea level by 2100 ranging from −5 to 43 cm of sea level equivalent under high carbon emission scenarios estimated by the recent Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). ISMIP6 highlighted the different behaviors of the East and West Antarctic ice sheets, as well as the possible role of increased surface mass balance in offsetting the dynamic ice loss in response to changing oceanic conditions in ice shelf cavities. However, the detailed contribution of individual glaciers, as well as the partitioning of uncertainty associated with this ensemble, have not yet been investigated. Here, we analyze the ISMIP6 results for high carbon emission scenarios, focusing on key glaciers around the Antarctic Ice Sheet, and we quantify their projected dynamic mass loss, defined here as mass loss through increased ice discharge into the ocean in response to changing oceanic conditions. We highlight glaciers contributing the most to sea level rise, as well as their vulnerability to changes in oceanic conditions. We then investigate the different sources of uncertainty and their relative role in projections, for the entire continent and for key individual glaciers. We show that, in addition to Thwaites and Pine Island glaciers in West Antarctica, Totten and Moscow University glaciers in East Antarctica present comparable future dynamic mass loss and high sensitivity to ice shelf basal melt. The overall uncertainty in additional dynamic mass loss in response to changing oceanic conditions, compared to a scenario with constant oceanic conditions, is dominated by the choice of ice sheet model, accounting for 52 % of the total uncertainty of the Antarctic dynamic mass loss in 2100. Its relative role for the most dynamic glaciers varies between 14 % for MacAyeal and Whillans ice streams and 56 % for Pine Island Glacier at the end of the century. The uncertainty associated with the choice of climate model increases over time and reaches 13 % of the uncertainty by 2100 for the Antarctic Ice Sheet but varies between 4 % for Thwaites Glacier and 53 % for Whillans Ice Stream. The uncertainty associated with the ice–climate interaction, which captures different treatments of oceanic forcings such as the choice of melt parameterization, its calibration, and simulated ice shelf geometries, accounts for 22 % of the uncertainty at the ice sheet scale but reaches 36 % and 39 % for Institute Ice Stream and Thwaites Glacier, respectively, by 2100. Overall, this study helps inform future research by highlighting the sectors of the ice sheet most vulnerable to oceanic warming over the 21st century and by quantifying the main sources of uncertainty.