Your search

During the Quaternary, ice sheets experienced several retreat–advance cycles, strongly influencing climate patterns. In order to properly simulate these phenomena, it is preferable to use physics-based models instead of parameterizations to estimate the surface mass balance (SMB), which strongly influences the evolution of the ice sheet. To further investigate the potential of these SMB models, this work evaluates the BErgen Snow SImulator (BESSI), a multi-layer snow model with high computational efficiency, as an alternative to providing the SMB for the Earth system model iLOVECLIM for multi-millennial simulations as in paleostudies. We compare the behaviors of BESSI and insolation temperature melt (ITM), an existing SMB scheme of iLOVECLIM during the Last Interglacial (LIG). Firstly, we validate the two SMB models using the regional climate model Mod- èle Atmosphérique Régional (MAR) as forcing and reference for the present-day climate over the Greenland and Antarctic ice sheets. The evolution of the SMB over the LIG (130–116 ka) is computed by forcing BESSI and ITM with transient climate forcing obtained from iLOVECLIM for both ice sheets. For present-day climate conditions, both BESSI and ITM exhibit good performance compared to MAR despite a much simpler model setup. While BESSI performs well for both Antarctica and Greenland for the same set of parame- ters, the ITM parameters need to be adapted specifically for each ice sheet. This suggests that the physics embedded in BESSI allows better capture of SMB changes across varying climate conditions, while ITM displays a much stronger sen- sitivity to its tunable parameters. The findings suggest that BESSI can provide more reliable SMB estimations for the iLOVECLIM framework to improve the model simulations of the ice sheet evolution and interactions with climate for multi-millennial simulations.

Water stable isotope records in polar ice cores have been largely used to reconstruct past local temperatures and other climatic information such as evaporative source region conditions of the precipitation reaching the ice core sites. However, recent studies have identified post-depositional processes taking place at the ice sheet's surface, modifying the original precipitation signal and challenging the traditional interpretation of ice core isotopic records. In this study, we use a combination of existing and new datasets of precipitation, snow surface, and subsurface isotopic compositions (δ18O and deuterium excess (d-excess)); meteorological parameters; ERA5 reanalyses; outputs from the isotope-enabled climate model ECHAM6-wiso; and a simple modelling approach to investigate the transfer function of water stable isotopes from precipitation to the snow surface and subsurface at Dome C in East Antarctica. We first show that water vapour fluxes at the surface of the ice sheet result in a net annual sublimation of snow, from 3.1 to 3.7 mm w.e. yr−1 (water equivalent) between 2018 and 2020, corresponding to 12 % to 15 % of the annual surface mass balance. We find that the precipitation isotopic signal cannot fully explain the mean, nor the variability in the isotopic composition observed in the snow, from annual to intra-monthly timescales. We observe that the mean effect of post-depositional processes over the study period enriches the snow surface in δ18O by 3.0 ‰ to 3.3 ‰ and lowers the snow surface d-excess by 3.4 ‰ to 3.5 ‰ compared to the incoming precipitation isotopic signal. We also show that the mean isotopic composition of the snow subsurface is not statistically different from that of the snow surface, indicating the preservation of the mean isotopic composition of the snow surface in the top centimetres of the snowpack. This study confirms previous findings about the complex interpretation of the water stable isotopic signal in the snow and provides the first quantitative estimation of the impact of post-depositional processes on the snow isotopic composition at Dome C, a crucial step for the accurate interpretation of isotopic records from ice cores.

Future climate and sea level projections depend sensitively on the response of the Antarctic Ice Sheet to ocean-driven melting and the resulting freshwater fluxes into the Southern Ocean. Circumpolar Deep Water (CDW) transport across the Antarctic continental shelf and into cavities beneath ice shelves is increasingly recognised as a crucial heat source for ice shelf melt. Quantifying past changes in the temperature of CDW is therefore of great benefit for modelling ice sheet response to past warm climates, for validating paleoclimate models, and for putting recent and projected changes in CDW temperature into context. Here we compile the available bottom water temperature reconstructions representative of CDW over the past 800 kyr. Estimated interglacial warming reached anomalies of +0.6 +/- 0.4 degrees C (MIS 11) and +0.5 +/- 0.5 degrees C (MIS 5) relative to present. Glacial cooling typically reached anomalies of ca. -1.5 to -2 degrees C, therefore maintaining positive thermal forcing for ice shelf melt even during glacials in the Amundsen Sea region of West Antarctica. Despite high variance amongst a small number of records and poor (4 kyr) temporal resolution, we find persistent and close relationships between our estimated CDW temperature and Southern Ocean sea surface temperature, Antarctic surface air temperature, and global deep-water temperature reconstructions at glacial-cycle timescales. Given the important role that CDW plays in connecting the world's three main ocean basins and in driving Antarctic Ice Sheet mass loss, additional temperature reconstructions targeting CDW are urgently needed to increase temporal and spatial resolution and to decrease uncertainty in past CDW temperatures - whether for use as a boundary condition, for model validation, or for understanding past oceanographic changes.

Earth’s obliquity and eccentricity cycles are strongly imprinted on Earth’s climate and widely used to measure geological time. However, the record of these imprints on the oxygen isotope record in deep-sea benthic foraminifera (δ18Ob) shows contradictory signals that violate isotopic principles and cause controversy over climate-ice sheet interactions. Here, we present a δ18Ob record of high fidelity from International Ocean Drilling Program (IODP) Site U1406 in the northwest Atlantic Ocean. We compare our record to other records for the time interval between 28 and 20 million years ago, when Earth was warmer than today, and only Antarctic ice sheets existed. The imprint of eccentricity on δ18Ob is remarkably consistent globally whereas the obliquity signal is inconsistent between sites, indicating that eccentricity was the primary pacemaker of land ice volume. The larger eccentricity-paced early Antarctic ice ages were vulnerable to rapid termination. These findings imply that the self-stabilizing hysteresis effects of large land-based early Antarctic ice sheets were strong enough to maintain ice growth despite consecutive insolation-induced polar warming episodes. However, rapid ice age terminations indicate that resistance to melting was weaker than simulated by numerical models and regularly overpowered, sometimes abruptly.

During the mid-Pliocene warm period (mPWP; 3.264–3.025 Ma), atmospheric CO2 concentrations were approximately 400 ppm, and the Antarctic Ice Sheet was substantially reduced compared to today. Antarctica is surrounded by the Southern Ocean, which plays a crucial role in the global oceanic circulation and climate regulation. Using results from the Pliocene Model Intercomparison Project (PlioMIP2), we investigate Southern Ocean conditions during the mPWP with respect to the pre-industrial period. We find that the mean sea surface temperature (SST) warming in the Southern Ocean is 2.8 °C, while global mean SST warming is 2.4 °C. The enhanced warming is strongly tied to a dramatic decrease in sea ice cover over the mPWP Southern Ocean. We also see a freshening of the ocean (sub)surface, driven by an increase in precipitation over the Southern Ocean and Antarctica. The warmer and fresher surface leads to a highly stratified Southern Ocean that can be related to weakening of the deep abyssal overturning circulation. Sensitivity simulations show that the decrease in sea ice cover and enhanced warming is largely a consequence of the reduction in the Antarctic Ice Sheet. In addition, the mPWP geographic boundary conditions are responsible for approximately half of the increase in mPWP SST warming, sea ice loss, precipitation, and stratification increase over the Southern Ocean. From these results, we conclude that a strongly reduced Antarctic Ice Sheet during the mPWP has a substantial influence on the state of the Southern Ocean and exacerbates the changes that are induced by a higher CO2 concentration alone. This is relevant for the long-term future of the Southern Ocean, as we expect melting of the western Antarctic Ice Sheet in the future, an effect that is not currently taken into account in future projections by Coupled Model Intercomparison Project (CMIP) ensembles.

Ice streams regulate most ice mass loss in Antarctica. Determining ice stream response to warmer conditions during the Pliocene could provide insights into their future behaviour, but this is hindered by a poor representation of subglacial topography in ice-sheet models. We address this limitation using a high-resolution model for Dronning Maud Land (East Antarctica). We show that contrary to dynamic thinning of the region’s ice streams following ice-shelf collapse, the largest ice stream, Jutulstraumen, thickens by 700 m despite lying on a retrograde bed slope. We attribute this counterintuitive thickening to a shallower Pliocene subglacial topography and inherent high lateral stresses at its flux gate. These conditions constrict ice drainage and, combined with increased snowfall, allow ice accumulation upstream. Similar stress balances and increased precipitation projections occur across 27% of present-day East Antarctica, and understanding how lateral stresses regulate ice-stream discharge is necessary for accurately assessing Antarctica’s future sea-level rise contribution.

New Zealand was among the last habitable places on earth to be colonized by humans. Charcoal records indicate that wildfires were rare prior to colonization and widespread following the 13th- to 14th-century Māori settlement, but the precise timing and magnitude of associated biomass-burning emissions are unknown, as are effects on light-absorbing black carbon aerosol concentrations over the pristine Southern Ocean and Antarctica. Here we used an array of well-dated Antarctic ice-core records to show that while black carbon deposition rates were stable over continental Antarctica during the past two millennia, they were approximately threefold higher over the northern Antarctic Peninsula during the past 700 years. Aerosol modelling demonstrates that the observed deposition could result only from increased emissions poleward of 40° S—implicating fires in Tasmania, New Zealand and Patagonia—but only New Zealand palaeofire records indicate coincident increases. Rapid deposition increases started in 1297 (±30 s.d.) in the northern Antarctic Peninsula, consistent with the late 13th-century Māori settlement and New Zealand black carbon emissions of 36 (±21 2 s.d.) Gg y−1 during peak deposition in the 16th century. While charcoal and pollen records suggest earlier, climate-modulated burning in Tasmania and southern Patagonia, deposition in Antarctica shows that black carbon emissions from burning in New Zealand dwarfed other preindustrial emissions in these regions during the past 2,000 years, providing clear evidence of large-scale environmental effects associated with early human activities across the remote Southern Hemisphere.

The dominant feature of large-scale mass transfer in the modern ocean is the Atlantic meridional overturning circulation (AMOC). The geometry and vigour of this circulation influences global climate on various timescales. Palaeoceanographic evidence suggests that during glacial periods of the past 1.5 million years the AMOC had markedly different features from today1; in the Atlantic basin, deep waters of Southern Ocean origin increased in volume while above them the core of the North Atlantic Deep Water (NADW) shoaled2. An absence of evidence on the origin of this phenomenon means that the sequence of events leading to global glacial conditions remains unclear. Here we present multi-proxy evidence showing that northward shifts in Antarctic iceberg melt in the Indian–Atlantic Southern Ocean (0–50° E) systematically preceded deep-water mass reorganizations by one to two thousand years during Pleistocene-era glaciations. With the aid of iceberg-trajectory model experiments, we demonstrate that such a shift in iceberg trajectories during glacial periods can result in a considerable redistribution of freshwater in the Southern Ocean. We suggest that this, in concert with increased sea-ice cover, enabled positive buoyancy anomalies to ‘escape’ into the upper limb of the AMOC, providing a teleconnection between surface Southern Ocean conditions and the formation of NADW. The magnitude and pacing of this mechanism evolved substantially across the mid-Pleistocene transition, and the coeval increase in magnitude of the ‘southern escape’ and deep circulation perturbations implicate this mechanism as a key feedback in the transition to the ‘100-kyr world’, in which glacial–interglacial cycles occur at roughly 100,000-year periods.

Understanding climate proxy records that preserve physical characteristics of past climate is a prerequisite to reconstruct long-term climatic conditions. Water stable isotope ratios (δ18O) constitute a widely used proxy in ice cores to reconstruct temperature and climate. However, the original climate signal is altered between the formation of precipitation and the ice, especially in low-accumulation areas such as the East Antarctic Plateau. Atmospheric conditions under which the isotopic signal is acquired at Aurora Basin North (ABN), East Antarctica, are characterized with the regional atmospheric model Modèle Atmosphérique Régional (MAR). The model shows that 50% of the snow is accumulated in less than 24 days year−1. Snowfall occurs throughout the year and intensifies during winter, with 64% of total accumulation between April and September, leading to a cold bias of −0.86°C in temperatures above inversion compared to the annual mean of −29.7°C. Large snowfall events are associated with high-pressure systems forcing warm oceanic air masses toward the Antarctic interior, which causes a warm bias of +2.83°C. The temperature-δ18O relationship, assessed with the global atmospheric model ECHAM5-wiso, is primarily constrained by the winter variability, but the observed slope is valid year-round. Three snow δ18O records covering 2004–2014 indicate that the anomalies recorded in the ice core are attributable to the occurrence of warm winter storms bringing precipitation to ABN and support the interpretation of δ18O in this region as a marker of temperature changes related to large-scale atmospheric conditions, particularly blocking events and variations in the Southern Annular Mode.