Search

The Hadley Centre Global Environment Model version 3 (HadGEM3) is the first coupled climate model to simulate an ice-free Arctic during the Last Interglacial (LIG), 127 000 years ago. This simulation appears to yield accurate Arctic surface temperatures during the summer season. Here, we investigate the causes and impacts of this extreme simulated ice loss. We find that the summer ice melt was predominantly driven by thermodynamic processes: atmospheric and ocean circulation changes did not significantly contribute to the ice loss. We demonstrate these thermodynamic processes were significantly impacted by melt ponds, which formed on average 8 d earlier during the LIG than during the pre-industrial control (PI) simulation. This relatively small difference significantly changed the LIG surface energy balance and impacted the albedo feedback. Compared to the PI simulation: in mid-June, of the absorbed flux at the surface over ice-covered cells (sea-ice concentration > 0.15), ponds accounted for 45 %–50 %, open water 35 %–45 %, and bare ice and snow 5 %–10 %. We show that the simulated ice loss led to large Arctic sea surface salinity and temperature changes. The sea surface temperature and salinity signals we identify here provide a means to verify, in marine observations, if and when an ice-free Arctic occurred during the LIG. Strong LIG correlations between spring melt pond and summer ice area indicate that, as Arctic ice continues to thin in future, the spring melt pond area will likely become an increasingly reliable predictor of the September sea-ice area. Finally, we note that models with explicitly modelled melt ponds seem to simulate particularly low LIG sea-ice area. These results show that models with explicit (as opposed to parameterised) melt ponds can simulate very different sea-ice behaviour under forcings other than the present day. This is of concern for future projections of sea-ice loss.

Changes in the Greenland ice sheet (GIS) affect global sea level. Greenland stable water isotope (δ18O) records from ice cores offer information on past changes in the surface of the GIS. Here, we use the isotope-enabled Hadley Centre Coupled Model version 3 (HadCM3) climate model to simulate a set of last interglacial (LIG) idealised GIS surface elevation change scenarios focusing on GIS ice core sites. We investigate how δ18O depends on the magnitude and sign of GIS elevation change and evaluate how the response is altered by sea ice changes. We find that modifying GIS elevation induces changes in Northern Hemisphere atmospheric circulation, sea ice and precipitation patterns. These climate feedbacks lead to ice-core-averaged isotopic lapse rates of 0.49 ‰ (100 m)−1 for the lowered GIS states and 0.29 ‰ (100 m)−1 for the enlarged GIS states. This is lower than the spatially derived Greenland lapse rates of 0.62–0.72 ‰ (100 m)−1. These results thus suggest non-linearities in the isotope–elevation relationship and have consequences for the interpretation of past elevation and climate changes across Greenland. In particular, our results suggest that winter sea ice changes may significantly influence isotope–elevation gradients: winter sea ice effect can decrease (increase) modelled core-averaged isotopic lapse rate values by about −19 % (and +28 %) for the lowered (enlarged) GIS states, respectively. The largest influence of sea ice on δ18O changes is found in coastal regions like the Camp Century site.

The Last Interglacial (LIG), a warmer period 130,000–116,000 years before present, is a potential analogue for future climate change. Stronger LIG summertime insolation at high northern latitudes drove Arctic land summer temperatures 4–5 °C higher than in the pre-industrial era. Climate model simulations have previously failed to capture these elevated temperatures, possibly because they were unable to correctly capture LIG sea-ice changes. Here, we show that the latest version of the fully coupled UK Hadley Center climate model (HadGEM3) simulates a more accurate Arctic LIG climate, including elevated temperatures. Improved model physics, including a sophisticated sea-ice melt-pond scheme, result in a complete simulated loss of Arctic sea ice in summer during the LIG, which has yet to be simulated in past generations of models. This ice-free Arctic yields a compelling solution to the long-standing puzzle of what drove LIG Arctic warmth and supports a fast retreat of future Arctic summer sea ice. Arctic climate in the Last Interglacial (LIG)—a warm period 130,000–116,000 years ago—is poorly simulated by modern climate models. A model with improved sea-ice melt-pond physics reproduces LIG Arctic temperatures, suggests an ice-free Arctic during this period and predicts the same by 2035.

When the same weather or climate simulation is run on different high-performance computing (HPC) platforms, model outputs may not be identical for a given initial condition. While the role of HPC platforms in delivering better climate projections is to some extent discussed in the literature, attention is mainly focused on scalability and performance rather than on the impact of machine-dependent processes on the numerical solution. Here we investigate the behaviour of the Preindustrial (PI) simulation prepared by the UK Met Office for the forthcoming CMIP6 (Coupled Model Intercomparison Project Phase 6) under different computing environments. Discrepancies between the means of key climate variables were analysed at different timescales, from decadal to centennial. We found that for the two simulations to be statistically indistinguishable, a 200-year averaging period must be used for the analysis of the results. Thus, constant-forcing climate simulations using the HadGEM3-GC3.1 model are reproducible on different HPC platforms provided that a sufficiently long duration of simulation is used. In regions where El Niño–Southern Oscillation (ENSO) teleconnection patterns were detected, we found large sea surface temperature and sea ice concentration differences on centennial timescales. This indicates that a 100-year constant-forcing climate simulation may not be long enough to adequately capture the internal variability of the HadGEM3-GC3.1 model, despite this being the minimum simulation length recommended by CMIP6 protocols for many MIP (Model Intercomparison Project) experiments. On the basis of our findings, we recommend a minimum simulation length of 200 years whenever possible.

Marine Isotope Stage (MIS) 5e (130–116 ka) represents an important ‘process analogue’ for understanding the climatic feedbacks and responses likely active under future anthropogenic warming. Reconstructing the Southern Ocean (SO) palaeoenvironment during MIS 5e and comparing it to the present day provides insights into the different responses of the SO sectors to a warmer climate. This study presents new records from seven marine sediment cores for MIS 5e together with their surface sediment records; all cores are located south of 55 oS. We investigate changes in diatom species assemblage and the accompanying variations in sea surface temperatures, winter sea-ice extent (WSIE) and glacial meltwater flux. All records show warmer conditions and a reduced WSIE during MIS 5e relative to the surface sediments. While the Pacific and Indian Sector records present very stable conditions throughout MIS 5e, the Atlantic Sector records display much more changeable conditions, particularly with respect to the WSIE. These variable conditions are attributed to higher iceberg and glacial meltwater flux in the Weddell Sea. This evidence for increased iceberg and glacial meltwater flux in the Weddell Sea during MIS 5e may have significant implications for understanding the stability of the West Antarctic Ice Sheet, both during MIS 5e and under future warming.

Environmental conditions during Marine Isotope Stage (MIS) 5e (130–116 ka) represent an important “process analogue” for understanding the climatic responses to present and future anthropogenic warming. The response of Antarctic sea ice to global warming is particularly uncertain due to the short length of the observational record. Reconstructing Antarctic winter sea-ice extent during MIS 5e therefore provides insights into the temporal and spatial patterns of sea-ice change under a warmer-than-present climate. This study presents new MIS 5e records from nine marine sediment cores located south of the Antarctic Polar Front between 55 and 70∘ S. Winter sea-ice extent and sea-surface temperatures are reconstructed using marine diatom assemblages and a modern analogue technique transfer function, and changes in these environmental variables between the three Southern Ocean sectors are investigated. The Atlantic and East Indian sector records show much more variable MIS 5e winter sea-ice extent and sea-surface temperatures than the Pacific sector records. High variability in the Atlantic sector winter sea-ice extent is attributed to high glacial meltwater flux in the Weddell Sea, indicated by increased abundances of the diatom species Eucampia antarctica and Fragilariopsis cylindrus. The high variability in the East Indian sector winter sea-ice extent is conversely believed to result from large latitudinal migrations of the flow bands of the Antarctic Circumpolar Current, inferred from latitudinal shifts in the sea-surface temperature isotherms. Overall, these findings suggest that Pacific sector winter sea ice displays a low sensitivity to warmer climates. The different variability and sensitivity of Antarctic winter sea-ice extent in the three Southern Ocean sectors during MIS 5e may have significant implications for the Southern Hemisphere climatic system under future warming.

Improving the understanding of changes in the mean and variability of climate variables as well as their interrelation is crucial for reliable climate change projections. Comparisons between general circulation models and paleoclimate archives using indirect proxies for temperature or precipitation have been used to test and validate the capability of climate models to represent climate changes. The oxygen isotopic ratio δ18O, a proxy for many different climate variables, is routinely measured in speleothem samples at decadal or higher resolution, and single specimens can cover full glacial–interglacial cycles. The calcium carbonate cave deposits are precisely dateable and provide well preserved (semi-)continuous albeit multivariate climate signals in the lower and mid-latitudes, where the measured δ18O in the mineral does not directly represent temperature or precipitation. Therefore, speleothems represent suitable archives to assess climate model abilities to simulate climate variability beyond the timescales covered by meteorological observations (101–102 years). Here, we present three transient isotope-enabled simulations from the Hadley Center Climate Model version 3 (iHadCM3) covering the last millennium (850–1850 CE) and compare them to a large global dataset of speleothem δ18O records from the Speleothem Isotopes Synthesis and AnaLysis (SISAL) database version 2 (Comas-Bru et al., 2020b). We systematically evaluate offsets in mean and variance of simulated δ18O and test for the main climate drivers recorded in δ18O for individual records or regions. The time-mean spatial offsets between the simulated δ18O and the speleothem data are fairly small. However, using robust filters and spectral analysis, we show that the observed archive-based variability of δ18O is lower than simulated by iHadCM3 on decadal and higher on centennial timescales. Most of this difference can likely be attributed to the records' lower temporal resolution and averaging or smoothing processes affecting the δ18O signal, e.g., through soil water residence times. Using cross-correlation analyses at site level and modeled grid-box level, we find evidence for highly variable but generally low signal-to-noise ratios in the proxy data. This points to a high influence of cave-internal processes and regional climate particularities and could suggest low regional representativity of individual sites. Long-range strong positive correlations dominate the speleothem correlation network but are much weaker in the simulation. One reason for this could lie in a lack of long-term internal climate variability in these model simulations, which could be tested by repeating similar comparisons with other isotope-enabled climate models and paleoclimate databases.

The 8.2 ka event is a well-known cooling event in the Northern Hemisphere, but is poorly understood in Madagascar. Here, we compare paleoclimate data and outputs from paleoclimate simulations to better understand it. Records from Madagascar suggest two distinct sub-events (8.3 ka and 8.2 ka), that seem to correlate with records from northern high latitude. This could indicate causal relationships via changes in the Atlantic Meridional Overturning Circulation (AMOC) with changes in moisture source's δ18O, and changes in the mean position of the Inter-Tropical Convergence Zone (ITCZ), as climate modelling suggests. These two sub-events are also apparent in other terrestrial records, but the climatic signals are different. The prominent 8.2 ka sub-event records a clear antiphase relationship between the northern and southern hemisphere monsoons, whereas such relationship is less evident during the first 8.3 ka sub-event. Data–model comparison have also shown a mismatch between the paleoclimate data and the model outputs, the causes of which are more or less understood and may lie in the proxies, in the model, or in both data and model. Knowing that paleoclimate proxies and climate models produce different sets of variables, further research is needed to improve the data–model comparison approach, so that both paleoclimate data and paleoclimate models will better predict the likely climate status of a region during a specified time in the past with minimal uncertainties.

The modeling of paleoclimate, using physically based tools, is increasingly seen as a strong out-of-sample test of the models that are used for the projection of future climate changes. New to the Coupled Model Intercomparison Project (CMIP6) is the Tier 1 Last Interglacial experiment for 127 000 years ago (lig127k), designed to address the climate responses to stronger orbital forcing than the midHolocene experiment, using the same state-of-the-art models as for the future and following a common experimental protocol. Here we present a first analysis of a multi-model ensemble of 17 climate models, all of which have completed the CMIP6 DECK (Diagnostic, Evaluation and Characterization of Klima) experiments. The equilibrium climate sensitivity (ECS) of these models varies from 1.8 to 5.6 ∘C. The seasonal character of the insolation anomalies results in strong summer warming over the Northern Hemisphere continents in the lig127k ensemble as compared to the CMIP6 piControl and much-reduced minimum sea ice in the Arctic. The multi-model results indicate enhanced summer monsoonal precipitation in the Northern Hemisphere and reductions in the Southern Hemisphere. These responses are greater in the lig127k than the CMIP6 midHolocene simulations as expected from the larger insolation anomalies at 127 than 6 ka. New synthesis for surface temperature and precipitation, targeted for 127 ka, have been developed for comparison to the multi-model ensemble. The lig127k model ensemble and data reconstructions are in good agreement for summer temperature anomalies over Canada, Scandinavia, and the North Atlantic and for precipitation over the Northern Hemisphere continents. The model–data comparisons and mismatches point to further study of the sensitivity of the simulations to uncertainties in the boundary conditions and of the uncertainties and sparse coverage in current proxy reconstructions. The CMIP6–Paleoclimate Modeling Intercomparison Project (PMIP4) lig127k simulations, in combination with the proxy record, improve our confidence in future projections of monsoons, surface temperature, and Arctic sea ice, thus providing a key target for model evaluation and optimization.

The modeling of paleoclimate, using physically based tools, is increasingly seen as a strong out-of-sample test of the models that are used for the projection of future climate changes. New to CMIP6 is the Tier 1 lig127k experiment, designed to address the climate responses to stronger orbital forcing than the midHolocene experiment, using the same state-of-the-art models and following a common experimental protocol. We present a multi-model ensemble of 17 climate models, all of which (except for two) have also completed the CMIP6 DECK experiments. The Equilibrium Climate Sensitivity (ECS) of these models varies from 2.1 to 5.6 °C. The seasonal character of the insolation anomalies results in strong warming over the Northern Hemisphere (NH) continents in the lig127k ensemble as compared to the piControl in June–July–August and a much-reduced minimum (August–September) summer sea ice extent in the Arctic. The multi-model results indicate enhanced summer monsoonal precipitation and areal extent in the Northern Hemisphere and reductions in the Southern Hemisphere. These responses are greater in the lig127k than midHolocene simulations as expected from the larger insolation anomalies at 127 ka than 6 ka. New syntheses for surface temperature and precipitation, targeted for 127 ka, have been developed for comparison to the multi-model ensemble. The lig127k model ensemble and data reconstructions are in good agreement for summer temperature anomalies over Canada, Scandinavia, and the North Atlantic and precipitation over the Northern Hemisphere continents. The model-data comparisons and mismatches point to further study of the sensitivity of the simulations to uncertainties in the specified boundary conditions and of the uncertainties and sparse coverage in current proxy reconstructions. The CMIP6-PMIP4 lig127k simulations, in combination with the proxy record, have potential implications for confidence in future projections of monsoons, surface temperature, Arctic sea ice, and the stability of the Greenland ice sheet.

The shape and extent of the Greenland Ice Sheet (GIS) during the Last Interglacial (LIG) is a matter of controversy, with different studies proposing a wide range of reconstructions. Here, for the first time, we combine stable water isotopic information from ice cores with isotope-enabled climate model outputs to investigate the problem. Exploring the space of possible ice sheet geometries by simulation is prohibitively expensive. We address this problem by using a Gaussian process emulator as a statistical surrogate of the full climate model. The emulator is calibrated using the results of a small number of carefully chosen simulations and then permits fast, probabilistic predictions of the simulator outputs at untried inputs. The inputs are GIS morphologies, parameterized through a dimension-reduction technique adapted to the spherical geometry of the setting. Based on the emulator predictions, the characteristics of morphologies compatible with the available ice core measurements are explored, leading to a reduction in uncertainty on the LIG GIS morphology. Moreover, a scenario-based approach allows to assess the gains in uncertainty reduction which would result from the availability of better dated LIG measurements in Greenland ice cores.

Palaeoclimate model simulations are an important tool to improve our understanding of the mechanisms of climate change. These simulations also provide tests of the ability of models to simulate climates very different to today. Here we present the results from two simulations using the latest version of the UK’s physical climate model, HadGEM3-GC3.1; the mid-Holocene (~ 6 ka) and Last Interglacial (~ 127 ka) simulations, both conducted under the auspices of CMIP6/PMIP4. These periods are of particular interest to PMIP4 because they represent the two most recent warm periods in Earth history, where atmospheric concentration of greenhouse gases and continental configuration is similar to the pre-industrial period but where there were significant changes to the Earth’s orbital configuration, resulting in a very different seasonal cycle of radiative forcing. Results for these simulations are assessed against proxy data, previous versions of the UK model, and models from the previous CMIP5 exercise. When the current version is compared to the previous generation of the UK model, the most recent version suggests limited improvement. In common with these previous model versions, the simulations reproduce global land and ocean temperatures (both surface and at 1.5 m) and a West African monsoon that is consistent with the latitudinal and seasonal distribution of insolation. The Last Interglacial simulation appears to accurately capture Northern Hemisphere temperature changes, but without the addition of Last Interglacial meltwater forcing cannot capture the magnitude of Southern Hemisphere changes. Model-data comparisons indicate that some geographical regions, and some seasons, produce better matches to the palaeodata (relative to pre-industrial) than others. Model-model comparisons, relative to previous generations same model and other models, indicate similarity between generations in terms of both the intensity and northward enhancement of the mid-Holocene West African monsoon, both of which are underestimated. On the Saharan greening which occurred the mid-Holocene African Humid Period, simulation results are likewise consistent with other models. The most recent version of the UK model appears to still be unable to reproduce the amount of rainfall necessary to support grassland across the Sahara.

The peak of the Last Interglacial, Marine Isotope Stage (MIS) 5e (130–116 ka), provides a valuable ‘process analogue’ for validating the climatic feedbacks and forcings likely active under future anthropogenic warming. Reconstructing exact timings of MIS 5e peak warming and minimum winter sea-ice extent (WSIE) throughout the Southern Ocean (SO) will help to identify the interactions and feedbacks within the ice-ocean system. Here we present a new MIS 5e marine sediment record from the SW Atlantic sector together with 28 published core records (chronologies standardised to the LR04 δ18O benthic stack; Lisiecki and Raymo, 2005) to investigate the timing and sequence of minimum WSIE and peak warming across the SO. Sea-surface temperatures (SSTs) peaked earliest in the Indian (20oE–150oE) and Atlantic (70oW–20oE) sectors, at 128.7 ± 0.8 ka and 127.4 ± 1.1 ka respectively, followed by the Pacific sector (150oE–70oW) at 124.9 ± 3.6 ka. The interval of minimum WSIE for all three sectors occurred within the period from 129–125 ka, consistent with the ∼128 ka sea salt flux minimum in Antarctic ice cores. Minimum WSIE appears to have coincided with peak July insolation at 55 oS, suggesting it could be linked with the mildest winters. The reduced WSIE during MIS 5e would have likely reduced the production of deep- and bottom water masses, inhibiting storage of CO2 in the abyssal ocean and lowering nutrient availability in SO surface waters. Examining a wide spatial range of proxy records for MIS 5e is a critical step forward in understanding climatic interactions and processes that will be active under warmer global temperatures.

Improving our knowledge of the temporal and spatial variability of the Antarctic Ice Sheet (AIS) surface mass balance (SMB) is crucial to reduce the uncertainties of past, present, and future Antarctic contributions to sea level rise. An examination of the surface air temperature–SMB relationship in model simulations demonstrates a strong link between the two. Reconstructions based on ice cores display a weaker relationship, indicating a model–data discrepancy that may be due to model biases or to the non-climatic noise present in the records. We find that, on the regional scale, the modeled relationship between surface air temperature and SMB is often stronger than between temperature and δ18O. This suggests that SMB data can be used to reconstruct past surface air temperature. Using this finding, we assimilate isotope-enabled SMB and δ18O model output with ice core observations to generate a new surface air temperature reconstruction. Although an independent evaluation of the skill is difficult because of the short observational time series, this new reconstruction outperforms the previous reconstructions for the continental-mean temperature that were based on δ18O alone. The improvement is most significant for the East Antarctic region, where the uncertainties are particularly large. Finally, using the same data assimilation method as for the surface air temperature reconstruction, we provide a spatial SMB reconstruction for the AIS over the last 2 centuries, showing large variability in SMB trends at a regional scale, with an increase (0.82 Gt yr−2) in West Antarctica over 1957–2000 and a decrease in East Antarctica during the same period (−0.13 Gt yr−2). As expected, this is consistent with the recent reconstruction used as a constraint in the data assimilation.

El Niño–Southern Oscillation (ENSO) is the strongest mode of interannual climate variability in the current climate, influencing ecosystems, agriculture, and weather systems across the globe, but future projections of ENSO frequency and amplitude remain highly uncertain. A comparison of changes in ENSO in a range of past and future climate simulations can provide insights into the sensitivity of ENSO to changes in the mean state, including changes in the seasonality of incoming solar radiation, global average temperatures, and spatial patterns of sea surface temperatures. As a comprehensive set of coupled model simulations is now available for both palaeoclimate time slices (the Last Glacial Maximum, mid-Holocene, and last interglacial) and idealised future warming scenarios (1 % per year CO2 increase, abrupt four-time CO2 increase), this allows a detailed evaluation of ENSO changes in this wide range of climates. Such a comparison can assist in constraining uncertainty in future projections, providing insights into model agreement and the sensitivity of ENSO to a range of factors. The majority of models simulate a consistent weakening of ENSO activity in the last interglacial and mid-Holocene experiments, and there is an ensemble mean reduction of variability in the western equatorial Pacific in the Last Glacial Maximum experiments. Changes in global temperature produce a weaker precipitation response to ENSO in the cold Last Glacial Maximum experiments and an enhanced precipitation response to ENSO in the warm increased CO2 experiments. No consistent relationship between changes in ENSO amplitude and annual cycle was identified across experiments.

Improving our knowledge of the temporal and spatial variability of the Antarctic Ice Sheet (AIS) Surface Mass Balance (SMB) is crucial to reduce the uncertainties of past, present and future Antarctic contributions to sea level rise. Here, we show that Global Climate Models (GCMs) can reproduce the present-day (1979–2005) AIS SMB and the temporal variations over the last two centuries. An examination of the surface temperature–SMB relationship in model simulations demonstrates a strong link between the two. Reconstructions based on ice cores display a weaker relationship, indicating a model-data discrepancy that may be due to model biases or to the non-climatic noise present in the records. We find that, on the regional scale, the modelled temperature-SMB relationship is stronger than the relationship between δ18O-temperature. This suggests that SMB data can be used to reconstruct past surface temperatures. Using this finding, we assimilate isotope-enabled model SMB and δ18O output with ice-core observations, to generate a new surface temperature reconstruction. Although an independent evaluation of the skill is difficult because of the short observational time series, this new reconstruction outperforms the previous reconstructions for the continental-mean temperature that were based on δ18O alone with a linear correlation coefficient with the observed surface temperatures (1958–2010 CE) of 0.73. The improvement is largest for the East Antarctic region, where the uncertainties are particularly large. Finally, we provide a spatial SMB reconstruction of the AIS over the last two centuries showing 1) large variability in SMB trends at regional scale; and 2) a large SMB increase (0.82 Gt year−2) in West Antarctica over 1957–2000 while at the same time, East Antarctica has experienced a large SMB decrease (−3.3 Gt year−2), which is consistent with a recent reconstruction.

When the same weather or climate simulation is run on different High Performance Computing (HPC) platforms, model outputs may not be identical for a given initial condition. While the role of HPC platforms in delivering better climate projections is often discussed in literature, attention is mainly focused on scalability and performance rather than on the impact of machine-dependent processes on the numerical solution. At the same time, machine dependence is an overlooked source of uncertainty when it comes to discussing the model spread observed within the Coupled Model Intercomparison Projects (CMIP). Here we investigate the impact of machine dependence on model results and quantify, for a selected case study, the magnitude of the uncertainty. We consider the Preindustrial (PI) simulation prepared by the UK Met Office for the forthcoming CMIP6. We compare key climate variables between PI control simulations run on the UK Met Office supercomputer and the ARCHER HPC platform. Discrepancies strongly depend on the timescale. Decadal means show substantial differences of up to 0.2 °C for global mean air temperature, 1 W/m2 for TOA outgoing longwave flux and 1.2 million km2 for Southern Hemisphere sea ice area. However, on multi-centennial timescales the differences are not significant and the long-term statistics of the two runs are similar. Differences between the two simulations can be linked to variations in the strongest modes of climate variability. In the Southern Hemisphere, this results in large SST anomalies where ENSO teleconnection patterns are expected that can reach 0.6 °C (and SNR > 1) even on centennial timescales.

Significance The Dansgaard–Oeschger events contained in Greenland ice cores constitute the archetypal record of abrupt climate change. An accurate understanding of these events hinges on interpretation of Greenland records of oxygen and nitrogen isotopes. We present here the important results from a suite of modeled Dansgaard–Oeschger events. These simulations show that the change in oxygen isotope per degree of warming becomes smaller during larger events. Abrupt reductions in sea ice also emerge as a strong control on ice core oxygen isotopes because of the influence on both the moisture source and the regional temperature increase. This work confirms the significance of sea ice for past abrupt warming events., Greenland ice cores provide excellent evidence of past abrupt climate changes. However, there is no universally accepted theory of how and why these Dansgaard–Oeschger (DO) events occur. Several mechanisms have been proposed to explain DO events, including sea ice, ice shelf buildup, ice sheets, atmospheric circulation, and meltwater changes. DO event temperature reconstructions depend on the stable water isotope (δ18O) and nitrogen isotope measurements from Greenland ice cores: interpretation of these measurements holds the key to understanding the nature of DO events. Here, we demonstrate the primary importance of sea ice as a control on Greenland ice core δ18O: 95% of the variability in δ18O in southern Greenland is explained by DO event sea ice changes. Our suite of DO events, simulated using a general circulation model, accurately captures the amplitude of δ18O enrichment during the abrupt DO event onsets. Simulated geographical variability is broadly consistent with available ice core evidence. We find an hitherto unknown sensitivity of the δ18O paleothermometer to the magnitude of DO event temperature increase: the change in δ18O per Kelvin temperature increase reduces with DO event amplitude. We show that this effect is controlled by precipitation seasonality.

Creating more accurate reconstructions of past Antarctic ice sheet shapes allows us to better predict how they will vary in the changing climate and contribute to future sea level changes. In this research, we use expert elicitation to create a subjective prior distribution of the Antarctic ice sheets at the Last Glacial Maximum (LGM), 21Ka. A design of shapes from this distribution will be run through the global climate model HadCM3, providing us with output that we can compare with proxy data to find a better estimate of the ice sheet shape at the LGM.

Last Interglacial (LIG), stable water isotope values (δ18O) measured in Greenland deep ice cores are at least 2.5‰ higher compared to the present day. Previous isotopic climate simulations of the LIG do not capture the observed Greenland δ18O increases. Here, we use the isotope-enabled HadCM3 (UK Met Office coupled atmosphere-ocean general circulation model) to investigate whether a retreat of Northern Hemisphere sea ice was responsible for this model-data disagreement. Our results highlight the potential significance of sea ice changes on the LIG Greenland isotopic maximum. Sea ice loss in combination with increased sea surface temperatures, over the Arctic, affect δ18O: water vapour enriched in heavy isotopes and a shorter distillation path may both increase δ18O values over Greenland. We show, for the first time, that simulations of the response to Arctic sea ice reduction are capable of producing the likely magnitude of LIG δ18O increases at NEEM, NGRIP, GIPS2 and Camp Century ice core sites. However, we may underestimate δ18O changes at the Renland, DYE3 and GRIP ice core locations. Accounting for possible ice sheet changes is likely to be required to produce a better fit to the LIG ice core δ18O values.

The Southern Hemisphere westerly winds (SHW) play an important role in regulating the capacity of the Southern Ocean carbon sink. They modulate upwelling of carbon-rich deep water and, with sea ice, determine the ocean surface area available for air–sea gas exchange. Some models indicate that the current strengthening and poleward shift of these winds will weaken the carbon sink. If correct, centennial- to millennial-scale reconstructions of the SHW intensity should be linked with past changes in atmospheric CO2, temperature and sea ice. Here we present a 12,300-year reconstruction of wind strength based on three independent proxies that track inputs of sea-salt aerosols and minerogenic particles accumulating in lake sediments on sub-Antarctic Macquarie Island. Between about 12.1 thousand years ago (ka) and 11.2 ka, and since about 7 ka, the wind intensities were above their long-term mean and corresponded with increasing atmospheric CO2. Conversely, from about 11.2 to 7.2 ka, the wind intensities were below their long-term mean and corresponded with decreasing atmospheric CO2. These observations are consistent with model inferences of enhanced SHW contributing to the long-term outgassing of CO2 from the Southern Ocean. The strength of the Southern Hemisphere westerly winds varied throughout the Holocene, according to a reconstruction from lake sediments, with periods of stronger winds coincident with higher atmospheric CO2 levels.

We compare multi-ice core data with δ18O model output for the early last interglacial Antarctic sea ice minimum. The spatial pattern of δ18O across Antarctica is sensitive to the spatial pattern of sea ice retreat. Local sea ice retreat increases the proportion of winter precipitation, depleting δ18O at ice core sites. However, retreat also enriches δ18O because of the reduced source-to-site distance for atmospheric vapor. The joint overall effect is for δ18O to increase as sea ice is reduced. Our data-model comparison indicates a winter sea ice retreat of 67, 59, and 43% relative to preindustrial in the Atlantic, Indian, and Pacific sectors of the Southern Ocean. A compilation of Southern Ocean sea ice proxy data provides weak support for this reconstruction. However, most published marine core sites are located too far north of the 128,000 years B.P. sea ice edge, preventing independent corroboration for this sea ice reconstruction. Plain Language Summary The Antarctic isotope and temperature maximum, which occurred approximately 128,000 years before present (B.P.) during the warmer than present last interglacial period, is associated with a major retreat of Antarctic sea ice. Understanding the details of this major sea ice retreat is crucial in order to understand the sensitivity of the Southern Hemisphere sea ice system and to evaluate the performance of climate model simulations in response to future warming. This work uses a multi-ice and ocean core data-model evaluation to assess the magnitude and spatial pattern of this sea ice retreat. Our results suggest that sea ice retreat was greatest in the Atlantic and Indian sectors of the Southern Ocean and less in the Pacific sector. These results may have had serious implications for the stability of marine terminating glaciers around the Antarctic Ice Sheet and their contribution to the last interglacial sea level rise. These results also support a hypothesized slowdown in northward ocean heat transport during the early last interglacial.

Warmer than present Antarctic and Southern Ocean temperatures during the last interglacial, approximately 128,000 years ago, have been attributed to changes in north-south ocean heat transport, causing opposing hemispheric temperature anomalies. We investigate the magnitude of Antarctic warming and Antarctic ice core isotopic enrichment in response to Northern Hemisphere meltwater input during the early last interglacial. A 1,600 year HadCM3 simulation driven by 0.25 Sv of meltwater input reproduces 50-60 % of the peak Southern Ocean summer sea surface temperature anomaly, sea ice retreat and ice core isotope enrichment. We also find a robust increase in the proportion of cold season precipitation during the last interglacial, leading to lower isotopic values at the Antarctic ice core sites. These results suggest that a HadCM3 simulation including 0.25 Sv for 3,000-4,000 years would reconcile the last interglacial observations, providing a potential solution for the last interglacial 'missing heat' problem.

Changes in the strength and position of Southern Hemisphere westerly winds during the Last Glacial cycle have been invoked to explain both millennial and glacial–interglacial climate fluctuations. However, neither paleo models nor paleodata agree on the magnitude, or even the sign, of the change in wind strength and latitude during the most studied glacial period, the Last Glacial Maximum (LGM), compared to the recent past. This paper synthesizes paleo-environmental data that have been used to infer changes in LGM winds. Data compilations are provided for changes in terrestrial moisture, dust deposition, sea surface temperatures and ocean fronts, and ocean productivity, and existing data on Southern Hemisphere ocean circulation changes during the LGM are summarized. We find that any hypothesis of LGM wind and climate change needs to provide a plausible explanation for increased moisture on the west coast of continents, cooler temperatures and higher productivity in the Subantarctic Zone, and reductions in Agulhas leakage around southern Africa. Our comparison suggests that an overall strengthening, an equatorward displacement, or no change at all in winds could all be interpreted as consistent with observations. If a single cause related to the southern westerlies is sought for all the evidence presented, then an equatorward displacement or strengthening of the winds would be consistent with the largest proportion of the observations. However, other processes, such as weakening or poleward shifts in winds, a weakened hydrological cycle, extended sea-ice cover, and changed buoyancy fluxes, cannot be ruled out as potential explanations of observed changes in moisture, surface temperature, and productivity. We contend that resolving the position and strength of westerly winds during the LGM remains elusive based on data reconstructions alone. However, we believe that these data reconstructions of environmental conditions can be used in conjunction with model simulations to identify which processes best represent westerly wind conditions during the LGM.

[1] Oxygen isotope palaeoclimate records, preserved in moss tissue cellulose, are complicated by environmental influences on the relationships between source water inputs and evaporative conditions. We carried out stable isotope analyses of precipitation collected from the maritime Antarctic and cellulose extracted from co-located Chorisodontium aciphyllum dominated moss peat bank deposits accumulated since 1870 A.D. Analyses of stable oxygen and hydrogen isotope composition of summer precipitation on Signy Island (60.7°S, 45.6°W) established a local meteoric water line (LMWL) similar to both the global MWL and other LMWLs, and almost identical to the HadAM3 isotope-enabled global circulation model output. The oxygen isotopic composition of cellulose (δ18OC) revealed little temporal variation between four moss peat banks on Signy Island since 1870. However, δ18OC followed two patterns with Sites A and D consistently 3‰ enriched relative to δ18OC values from Sites B and C. The growing moss surfaces at Sites A and D are likely to have been hydrated by isotopically heavier summer precipitation, whilst at Sites B and C, the moss banks are regularly saturated by the isotopically depleted snow melt streams. Laboratory experiments revealed that evaporative enrichment of C. aciphyllum moss leaf water by 5‰ occurred rapidly following saturation (ecologically equivalent to post-rainfall or snow melt periods). In addition to the recognized source water-cellulose fractionation extent of 27 ± 3‰, such a shift would account for the 32‰ difference measured between δ18O of Signy Island precipitation and cellulose.

Latitudinal shifts in the Southern Ocean westerly wind jet could drive changes in the glacial to interglacial ocean CO2 inventory. However, whilst CMIP5 model results feature consistent future-warming jet shifts, there is considerable disagreement in deglacial-warming jet shifts. We find here that the dependence of pre-industrial (PI) to Last Glacial Maximum (LGM) jet shifts on PI jet position, or state dependency, explains less of the shifts in jet simulated by the models for the LGM compared with future-warming scenarios. State dependence is also weaker for intensity changes, compared to latitudinal shifts in the jet. Winter sea ice was considerably more extensive during the LGM, along with state dependence, changes in surface heat fluxes due this sea ice change probably had a large impact on the jet. Models which both simulate realistically large expansions in sea ice and feature PI jets which are south of 50 °S, show an increase in wind speed around 55 °S, and can show a poleward shift in the jet between the PI and the LGM. However models with the PI jet positioned equatorwards of around 47 °S do not show this response: the sea ice edge is too far from the jet for it to respond. In models with accurately positioned PI jets; a +1° difference in the latitude of the sea ice edge tends to be associated with a −0.85° shift in the 850 hPa jet. However, it seems that around 5° of expansion of LGM sea ice is necessary to hold the jet in its PI position. Since observational data supports an expansion of more than 5°, this result suggests that a slight poleward shift and intensification, was the most likely jet change between the PI and the LGM. Without the effect of sea ice, models simulated polewards shifted westerlies in warming climates and equatorward shifted westerlies in colder climates. However, the feedback of sea ice counters and reverses the equatorward trend in cooler climates so that the LGM winds were more likely to have also been shifted slightly poleward.

Reconstructions of salinity are used to diagnose changes in the hydrological cycle and ocean circulation. A widely used method of determining past salinity uses oxygen isotope (δOw) residuals after the extraction of the global ice volume and temperature components. This method relies on a constant relationship between δOw and salinity throughout time. Here we use the isotope-enabled fully coupled General Circulation Model (GCM) HadCM3 to test the application of spatially and time-independent relationships in the reconstruction of past ocean salinity. Simulations of the Late Holocene (LH), Last Glacial Maximum (LGM), and Last Interglacial (LIG) climates are performed and benchmarked against existing compilations of stable oxygen isotopes in carbonates (δOc), which primarily reflect δOw and temperature. We find that HadCM3 produces an accurate representation of the surface ocean δOc distribution for the LH and LGM. Our simulations show considerable variability in spatial and temporal δOw-salinity relationships. Spatial gradients are generally shallower but within ~50% of the actual simulated LH to LGM and LH to LIG temporal gradients and temporal gradients calculated from multi-decadal variability are generally shallower than both spatial and actual simulated gradients. The largest sources of uncertainty in salinity reconstructions are found to be caused by changes in regional freshwater budgets, ocean circulation, and sea ice regimes. These can cause errors in salinity estimates exceeding 4 psu. Our results suggest that paleosalinity reconstructions in the South Atlantic, Indian and Tropical Pacific Oceans should be most robust, since these regions exhibit relatively constant δOw-salinity relationships across spatial and temporal scales. Largest uncertainties will affect North Atlantic and high latitude paleosalinity reconstructions. Finally, the results show that it is difficult to generate reliable salinity estimates for regions of dynamic oceanography, such as the North Atlantic, without additional constraints.

The sensitivity of interglacial Greenland temperature to orbital and CO2 forcing is investigated using the NorthGRIP ice core data and coupled ocean-atmosphere IPSL-CM4 model simulations. These simulations were conducted in response to different interglacial orbital configurations, and to increased CO2 concentrations. These different forcings cause very distinct simulated seasonal and latitudinal temperature and water cycle changes, limiting the analogies between the last interglacial and future climate. However, the IPSL-CM4 model shows similar magnitudes of Arctic summer warming and climate feedbacks in response to 2 × CO2 and orbital forcing of the last interglacial period (126 000 years ago). The IPSL-CM4 model produces a remarkably linear relationship between TOA incoming summer solar radiation and simulated changes in summer and annual mean central Greenland temperature. This contrasts with the stable isotope record from the Greenland ice cores, showing a multi-millennial lagged response to summer insolation. During the early part of interglacials, the observed lags may be explained by ice sheet-ocean feedbacks linked with changes in ice sheet elevation and the impact of meltwater on ocean circulation, as investigated with sensitivity studies. A quantitative comparison between ice core data and climate simulations requires stability of the stable isotope – temperature relationship to be explored. Atmospheric simulations including water stable isotopes have been conducted with the LMDZiso model under different boundary conditions. This set of simulations allows calculation of a temporal Greenland isotope-temperature slope (0.3–0.4‰ per °C) during warmer-than-present Arctic climates, in response to increased CO2, increased ocean temperature and orbital forcing. This temporal slope appears half as large as the modern spatial gradient and is consistent with other ice core estimates. It may, however, be model-dependent, as indicated by preliminary comparison with other models. This suggests that further simulations and detailed inter-model comparisons are also likely to be of benefit. Comparisons with Greenland ice core stable isotope data reveals that IPSL-CM4/LMDZiso simulations strongly underestimate the amplitude of the ice core signal during the last interglacial, which could reach +8–10 °C at fixed-elevation. While the model-data mismatch may result from missing positive feedbacks (e.g. vegetation), it could also be explained by a reduced elevation of the central Greenland ice sheet surface by 300–400 m.

We present a novel automated processing method for obtaining layer dip from radio-echo sounding (RES) data. The method is robust, easily applicable and can be used to process large (several terabytes) ground and airborne RES datasets using modest computing resources. We give test results from the application of the method to two Antarctic datasets: the Fletcher Promontory ground-based radar dataset and the Wilkes Subglacial Basin airborne radar dataset. The automated RES processing (ARESP) method comprises the basic steps: (1) RES noise reduction; (2) radar layer identification; (3) isolation of individual ‘layer objects’; (4) measurement of orientation and other object properties; (5) elimination of noise in the orientation data; and (6) collation of the valid dip information. The apparent dip datasets produced by the method will aid glaciologists seeking to understand ice-flow dynamics in Greenland and Antarctica: ARESP could enable a shift from selective regional case studies to ice-sheet-scale studies.

Reconstructions of temperature variations from Antarctic ice cores rely on the assumption that the relationship between hydrogen and oxygen isotopic ratios and temperature are stable in space and time. Sime et al. analyse three 340,000-year-old ice cores from East Antarctica and use an isotope-enabled general circulation model to show that instead, the relationship is nonlinear. During warm periods, the ratios are less sensitive to temperature, so previous estimates of interglacial temperatures are likely to be about 3 °C too low. This is consistent with peak Antarctic interglacial temperatures at least 6 °C higher than today. This work suggests that there are serious deficiencies in our understanding of climates that are warmer than today's. Reconstructions of temperature variations from Antarctic ice cores rely on the assumption that the relationship between hydrogen and oxygen isotope ratios and temperature are stable in space and time. Three East Antarctic ice core records are now analysed alongside input from general circulation models to reveal that during warmer interglacial periods the isotope ratios are less sensitive to temperature than during colder interglacials; consequently, previous estimates of interglacial temperatures are probably too cold. Stable isotope ratios of oxygen and hydrogen in the Antarctic ice core record have revolutionized our understanding of Pleistocene climate variations and have allowed reconstructions of Antarctic temperature over the past 800,000 years (800 kyr; refs 1, 2). The relationship between the D/H ratio of mean annual precipitation and mean annual surface air temperature is said to be uniform ±10% over East Antarctica3 and constant with time ±20% (refs 3–5). In the absence of strong independent temperature proxy evidence allowing us to calibrate individual ice cores, prior general circulation model (GCM) studies have supported the assumption of constant uniform conversion for climates cooler than that of the present day3,5. Here we analyse the three available 340 kyr East Antarctic ice core records alongside input from GCM modelling. We show that for warmer interglacial periods the relationship between temperature and the isotopic signature varies among ice core sites, and that therefore the conversions must be nonlinear for at least some sites. Model results indicate that the isotopic composition of East Antarctic ice is less sensitive to temperature changes during warmer climates. We conclude that previous temperature estimates from interglacial climates are likely to be too low. The available evidence is consistent with a peak Antarctic interglacial temperature that was at least 6 K higher than that of the present day —approximately double the widely quoted 3 ± 1.5 K (refs 5, 6).

Several studies have suggested that sea-level rise during the last interglacial implies retreat of the West Antarctic Ice Sheet (WAIS). The prevalent hypothesis is that the retreat coincided with the peak Antarctic temperature and stable water isotope values from 128,000 years ago (128 ka); very early in the last interglacial. Here, by analysing climate model simulations of last interglacial WAIS loss featuring water isotopes, we show instead that the isotopic response to WAIS loss is in opposition to the isotopic evidence at 128 ka. Instead, a reduction in winter sea ice area of 65±7% fully explains the 128 ka ice core evidence. Our finding of a marked retreat of the sea ice at 128 ka demonstrates the sensitivity of Antarctic sea ice extent to climate warming., The peak in Antarctic ice core isotope values, 128,000 years before present, was concurrent with a significantly warmer-than-present Antarctic climate. Here, the authors show that this isotope maximum was associated with a major retreat of sea ice and not a collapse of the West Antarctic Ice Sheet.

International audience; Measurements of Last Interglacial stable water isotopes in ice cores show that central Greenland d18O increased by at least 3‰ compared to present day. Attempting to quantify the Greenland interglacial temperature change from these ice core measurements rests on our ability to interpret the stable water isotope content of Greenland snow. Current orbitally driven interglacial simulations do not show d18O or temperature rises of the correct magnitude, leading to difficulty in using only these experiments to inform our understanding of higher interglacial d18O. Here, analysis of greenhouse gas warmed simulations from two isotope-enabled general circulation models, in conjunction with a set of Last Interglacial sea surface observations, indicates a possible explanation for the interglacial d18O rise. A reduction in the winter time sea ice concentration around the northern half of Greenland, together with an increase in sea surface temperatures over the same region, is found to be sufficient to drive a >3‰ interglacial enrichment in central Greenland snow. Warm climate d18O and dD in precipitation falling on Greenland are shown to be strongly influenced by local sea surface condition changes: local sea surface warming and a shrunken sea ice extent increase the proportion of water vapour from local (isotopically enriched) sources, compared to that from distal (isotopically depleted) sources. Precipitation intermittency changes, under warmer conditions, leads to geographical variability in the d18O against temperature gradients across Greenland. Little sea surface warming around the northern areas of Greenland leads to low d18O against temperature gradients (0.1-0.3‰ per °C), whilst large sea surface warmings in these regions leads to higher gradients (0.3-0.7‰ per °C). These gradients imply a wide possible range of present day to interglacial temperature increases (4 to >10 °C). Thus, we find that uncertainty about local interglacial sea surface conditions, rather than precipitation intermittency changes, may lead to the largest uncertainties in interpreting temperature from Greenland ice cores. We find that interglacial sea surface change observational records are currently insufficient to enable discrimination between these different d18O against temperature gradients. In conclusion, further information on interglacial sea surface temperatures and sea ice changes around northern Greenland should indicate whether +5 °C during the Last Interglacial is sufficient to drive the observed ice core d18O increase, or whether a larger temperature increases or ice sheet changes are also required to explain the ice core observations.

The Southern Hemisphere (SH) westerly winds are thought to be critical to global ocean circulation, productivity, and carbon storage. For example, an equatorward shift in the winds, though its affect on the Southern Ocean circulation, has been suggested as the leading cause for the reduction in atmospheric CO2 during the Last Glacial period. Despite the importance of the winds, it is currently not clear, from observations or model results, how they behave during the Last Glacial. Here, an atmospheric modelling study is performed to help determine likely changes in the SH westerly winds during the Last Glacial Maximum (LGM). Using LGM boundary conditions, the maximum in SH westerlies is strengthened by ∼+1 m s−1 and moved southward by ∼2° at the 850 hPa pressure level. Boundary layer stabilisation effects over equatorward extended LGM sea-ice can lead to a small apparent equatorward shift in the wind band at the surface. Further sensitivity analysis with individual boundary condition changes indicate that changes in sea surface temperatures are the strongest factor behind the wind change. The HadAM3 atmospheric simulations, along with published PMIP2 coupled climate model simulations, are then assessed against the newly synthesised database of moisture observations for the LGM. Although the moisture data is the most commonly cited evidence in support of a large equatorward shift in the SH winds during the LGM, none of the models that produce realistic LGM precipitation changes show such a large equatorward shift. In fact, the model which best simulates the moisture proxy data is the HadAM3 LGM simulation which shows a small poleward wind shift. While we cannot prove here that a large equatorward shift would not be able to reproduce the moisture data as well, we show that the moisture proxies do not provide an observational evidence base for it.

We know from the ice record that the concentration of atmospheric methane, [CH4], at the Last Glacial Maximum (LGM) was roughly half that in the pre-industrial era (PI), buthow much of the difference was source-driven, and how much was sink-driven, remains uncertain. Recent developments include: a higher estimate of the LGM-PI change in methane emissions from wetlands―the dominant, natural methane source; and the possible recycling of OH consumed in isoprene oxidation―the principal methane sink. Here, in view of these developments, we use an atmospheric chemistry-transport model to re-examine the main factors affecting OH during this period: changes in air temperature and emissions of non-methane volatile organic compounds from vegetation. We find that their net effect was negligible(with and without an OH recycling mechanism), implyingthe change in [CH4] was almost entirely source driven―a conclusion that, though subject to significant uncertainties,can be reconciled with recent methane source estimates.

The resemblance of the orbitally filtered isotope signal from the past 340 kyr in Antarctic ice cores to Northern Hemisphere summer insolation intensity has been used to suggest that the northern hemisphere may drive orbital-scale global climate changes. A recent Letter by Laepple et al. suggests that, contrary to this interpretation, this semblance may instead be explained by weighting the orbitally controlled Antarctic seasonal insolation cycle with a static (present-day) estimate of the seasonal cycle of accumulation. We suggest, however, that both time variability in accumulation seasonality and alternative stable seasonality can markedly alter the weighted insolation signal. This indicates that, if the last 340 kyr of Antarctic accumulation has not always looked like the estimate of precipitation and accumulation seasonality made by Laepple et al., this particular accumulation weighting explanation of the Antarctic orbital-scale isotopic signal might not be robust.

Recent developments in ice melter systems and continuous flow analysis (CFA) techniques now allow higher-resolution ice core analysis. Here, we present a new method to aid interpretation of high-resolution ice core stable water isotope records. Using a set of simple isotopic recording and postdepositional assumptions, the European Centre for Medium-Range Weather Forecasts' 40 year reanalysis time series of temperature and precipitation are converted to “virtual core” depth series across the Antarctic Peninsula, helping us to understand what information can be gleaned from the CFA high-resolution observations. Virtual core temperatures are transferred onto time using three different depth-age transfer assumptions: (1) a perfect depth-age model, (2) a depth-age model constructed from single or dual annual photochemical tie points, and (3) a cross-dated depth-age model. Comparing the sampled temperatures on the various depth-age models with the original time series allows quantification of the effect of ice core sample resolution and dating. We show that accurate annual layer count depth-age models should allow some subseasonal temperature anomalies to be recovered using a sample resolution of around 40 mm, or 10–20 samples per year. Seasonal temperature anomalies may be recovered using sample lengths closer to 60 mm, or about 7–14 samples per year. These results tend to confirm the value of current CFA ice core sampling strategies and indicate that it should be possible to recover about a third of subannual (but not synoptic) temperature anomaly information from annually “layer-counted” peninsula ice cores.

We compare the present and last interglacial periods as recorded in Antarctic water stable isotope records now available at various temporal resolutions from six East Antarctic ice cores: Vostok, Taylor Dome, EPICA Dome C (EDC), EPICA Dronning Maud Land (EDML), Dome Fuji and the recent TALDICE ice core from Talos Dome. We first review the different modern site characteristics in terms of ice flow, meteorological conditions, precipitation intermittency and moisture origin, as depicted by meteorological data, atmospheric reanalyses and Lagrangian moisture source diagnostics. These different factors can indeed alter the relationships between temperature and water stable isotopes. Using five records with sufficient resolution on the EDC3 age scale, common features are quantified through principal component analyses. Consistent with instrumental records and atmospheric model results, the ice core data depict rather coherent and homogenous patterns in East Antarctica during the last two interglacials. Across the East Antarctic plateau, regional differences, with respect to the common East Antarctic signal, appear to have similar patterns during the current and last interglacials. We identify two abrupt shifts in isotopic records during glacial inception at TALDICE and EDML, likely caused by regional sea ice expansion. These regional differences are discussed in terms of moisture origin and in terms of past changes in local elevation histories which are compared to ice sheet model results. Our results suggest that, for coastal sites, elevation changes may contribute significantly to inter-site differences. These elevation changes may be underestimated by current ice sheet models.

The ice-core record of the carbon-13 content of atmospheric methane (δ13CH4) has largely been used to constrain past changes in methane sources. The aim of this paper is to explore, for the first time, the contribution that changes in the strength of a minor methane sink―oxidation by atomic chlorine in the marine boundary layer (ClMBL)―could make to changes in δ13CH4 on glacial-interglacial timescales. Combining wind and temperature data from a variety of general circulation models with a simple formulation for the concentration of ClMBL, we find that changes in the strength of this sink, driven solely by changes in the atmospheric circulation, could have been responsible for changes in δ13CH4 of the order of 10% of the glacial-interglacial difference observed. We thus highlight the need to quantify past changes in the strength of this sink, including those relating to changes in the sea-ice source of sea salt aerosol.

[1] This study uses ice core methanesulphonic acid (MSA) records from the Antarctic Peninsula, where temperatures have been warming faster than anywhere else in the Southern Hemisphere, to reconstruct the 20th century history of sea ice change in the adjacent Bellingshausen Sea. Using satellite‐derived sea ice and meteorological data, we show that ice core MSA records from this region are a reliable proxy for regional sea ice change, with years of increased winter sea ice extent recorded by increased ice core MSA concentrations. Our reconstruction suggests that the satellite‐observed sea ice decline in the Bellingshausen Sea during recent decades is part of a long‐term regional trend that has occurred throughout the 20th century. The long‐term perspective on sea ice in the Bellingshausen Sea is consistent with evidence of 20th century warming on the Antarctic Peninsula and may reflect a progressive deepening of the Amundsen Sea Low due to increasing greenhouse gas concentrations and, more recently, stratospheric ozone depletion. As a first‐order estimate, our MSA‐based reconstruction suggests that sea ice in the Bellingshausen Sea has retreated southward by ∼0.7° during the 20th century. Comparison with other 20th century sea ice observations, reconstructions, and model simulations provides a coherent picture of Antarctic sea ice decline during the 20th century, although with regional‐scale differences evident in the timing and magnitude of this sea ice decline. This longer‐term perspective contrasts with the small overall increase in Antarctic sea ice that is observed in post‐1979 satellite data.

International audience; The EPICA ice core from Dome C extends 3259 m in depth, and encompasses 800 ka of datable and sequential ice. Numerous chemical species have been measured along the length of the cores. Here we concentrate on interpreting the main low-resolution patterns of major ions. We extend the published record for non-sea-salt calcium, sea-salt sodium and non-sea-salt sulfate flux to 800 ka. The non-sea-salt calcium record confirms that terrestrial dust originating from South America closely mirrored Antarctic climate, both at orbital and millennial timescales. A major cause of the main trends is most likely climate in southern South America, which could be sensitive to subtle changes in atmospheric circulation. Sea-salt sodium also follows temperature, but with a threshold at low temperature. We re-examine the use of sodium as a sea ice proxy, concluding that it is probably reflecting extent, with high salt concentrations reflecting larger ice extents. With this interpretation, the sodium flux record indicates low ice extent operating as an amplifier in warm interglacials. Non-sea-salt sulfate flux is almost constant along the core, confirming the lack of change in marine productivity (for sulfur-producing organisms) in the areas of the Southern Ocean contributing to the flux at Dome C. For the first time we also present long records of reversible species such as nitrate and chloride, and show that the pattern of post-depositional losses described for shallower ice is maintained in older ice. It appears possible to use these concentrations to constrain snow accumulation rates in interglacial ice at this site, and the results suggest a possible correction to accumulation rates in one early interglacial. Taken together the chemistry records offer a number of constraints on the way the Earth system combined to give the major climate fluctuations of the late Quaternary period.

Stable water isotopes have been added to the full hydrological cycle of the Hadley Centre Climate model (HadCM3) coupled atmosphere-ocean GCM. Simulations of delta O-18 in precipitation and at the ocean surface compare well with observations for the present-day climate. The model has been used to investigate the isotopic anomalies associated with ENSO; it is found that the anomalous delta O-18 in precipitation is correlated with the anomalous precipitation amount in accordance with the "amount effect.'' The El Nino delta O-18 anomaly at the ocean surface is largest in coastal regions because of the mixing of ocean water and the more depleted runoff from the land surface. Coral delta O-18 anomalies were estimated, using an established empirical relationship, and generally reflect ocean surface delta O-18 anomalies in coastal regions and sea surface temperatures away from the coast. The spatial relationship between tropical precipitation and delta O-18 was investigated for the El Nino anomaly simulated by HadCM3. Weighting the El Nino precipitation anomaly by the precipitation amount at each grid box gave a large increase in the spatial correlation between tropical precipitation and delta O-18. This improvement was most apparent over land points and between 10 and 20 degrees of latitude.

Analysis of ERA40 temperature and accumulation data suggests that annual mean isotopic fluctuations due to temperature change will be geographically very variable across the Peninsula: isotopic variations of 0.4 parts per thousand at James Ross Island; 0.9 parts per thousand at Dyer; and 1.3 parts per thousand at Gomez are all likely to indicate an identical magnitude of temperature change. The reduction in the magnitude of the isotopic signal in the north and east is due to climatically dependent synoptic covariance between temperature and accumulation; whilst in the west and south seasonal covariance amplifies the isotopic temperature signal. Additionally we show that the relationship between accumulation and temperature is rather weak in the north-east regions but is stronger in the central and southerly regions. Therefore isotopes may record 11% to 30% of the variance in annual mean temperatures in the north east; 75% in central regions; and 70% in the south. This study enables physically based reconstructions of Peninsula climate based on multi-core analysis. Citation: Sime, L. C., G. J. Marshall, R. Mulvaney, and E. R. Thomas ( 2009), Interpreting temperature information from ice cores along the Antarctic Peninsula: ERA40 analysis, Geophys. Res. Lett., 36, L18801, doi:10.1029/2009GL038982.

[1] Results from an isotope-enabled general circulation model are presented in order to determine the isotopic signal of a warmer climate on Antarctica. The warming is forced using CO2 forecasts for the next century. For unforced interannual climate variability the temporal gradient and correlation between stable water isotopes and surface temperature is small. The relationship is much stronger for the CO2 forced event. There is little regional coherence between temporal gradients for the forced and unforced climates, implying that correlations between stable water isotopes and temperature from instrumental records of a couple of decades cannot be applied to larger warming events. Additionally, there are strong discrepancies between the forced warming temporal gradients and present-day spatial gradients of isotopes against temperature. We show that it is difficult to obtain a local spatial gradient since it is systematically affected by the geographical size of the spatial sample. For the forced warming, the temporal gradient derived for the warming event over Dome C is less than half the value generally applied. We determine, through means of a new frequency decomposition, that a large portion of this decrease from the expected value is due to changes in the seasonal precipitation temperature covariance. This low isotopic sensitivity to a CO2 driven warming implies that current and future warming trends may have rather small isotopic signals in Antarctica.

[1] Moving boat acoustic Doppler current profiling (ADCP) is increasingly used to measure discharge in large rivers. We investigate whether useful information about bed shear stress can be recovered from such data. Alternative ways to estimate local bed shear stress using the logarithmic law of the wall and spatial averaging are tested using ADCP transects across lower Fraser River, Canada. Repeatability is assessed by comparing estimates from outward and return boat tracks. The most precise method uses the vertically averaged mean velocity and a zero-velocity height based on bed grain size information. The accuracy of the assumed zero-velocity height can be judged by consistency between estimates using mean velocity and near-bed velocity. Shear stress estimates from unconstrained log-law fits are less repeatable and tend to overpredict, and mean shear stress estimates using the depth-slope product are unreliable in this river because of nonuniform flow.

Plain Language Summary\ud The Last Interglacial period (LIG, 116,000 to 130,000 years ago) was globally ∼ 0.8 °C warmer than today at its peak, with substantially more warming at the poles. It is a valuable analogue for future global temperature rise, especially for understanding rates and sources of polar ice melt and subsequent global sea level rise. Records of water stable isotopes from Antarctic ice cores have been crucial for understanding past polar temperature during the LIG. However we currently lack a framework for estimating how changes in the ice sheet elevation, alongside sea‐ice feedbacks, affect these water stable isotopes. To address this, we examine the effect of the Antarctic Ice Sheet (AIS) elevation on water stable isotopes, using an ensemble of climate simulations where we vary the AIS elevation. We observe that (i) water stable isotope values lower with increasing AIS elevation following linear relationships, (ii) the effect of sea‐ice induced by AIS elevation is small so the effect of AIS elevation can be isolated. Finally, this study provides appropriate elevation‐water stable isotope gradients for the reconstruction of the AIS topography using ice cores.\ud Abstract\ud Changes of the topography of the Antarctic ice sheet (AIS) can complicate the interpretation of ice core water stable isotope measurements in terms of temperature. Here, we use a set of idealised AIS elevation change scenarios to investigate this for the warm Last Interglacial (LIG). We show that LIG δ 18 O against elevation relationships are not uniform across Antarctica, and that the LIG response to elevation is lower than the preindustrial response. The effect of LIG elevation‐induced sea ice changes on δ 18 O is small, allowing us to isolate the effect of elevation change alone. Our results help to define the effect of AIS changes on the LIG δ 18 O signals, and should be invaluable to those seeking to use AIS ice core measurements for these purposes. Especially, our simulations strengthen the conclusion that ice core measurements from the Talos Dome core exclude the loss of the Wilkes Basin at around 128 ky.

The many different proxy records from the European Project for Ice Coring in Antarctica (EPICA) Dome C ice core allow for the first time a comparison of nine glacial terminations in great detail. Despite the fact that all terminations cover the transition from a glacial maximum into an interglacial, there are large differences between single terminations. For some terminations, Antarctic temperature increased only moderately, while for others, the amplitude of change at the termination was much larger. For the different terminations, the rate of change in temperature is more similar than the magnitude or duration of change. These temperature changes were accompanied by vast changes in dust and sea salt deposition all over Antarctica. Here we investigate the phasing between a South American dust proxy (non-sea-salt calcium flux, nssCa2+), a sea ice proxy (sea salt sodium flux, ssNa+) and a proxy for Antarctic temperature (deuterium, δD). In particular, we look into whether a similar sequence of events applies to all terminations, despite their different characteristics. All proxies are derived from the EPICA Dome C ice core, resulting in a relative dating uncertainty between the proxies of less than 20 years. At the start of the terminations, the temperature (δD) increase and dust (nssCa2+ flux) decrease start synchronously. The sea ice proxy (ssNa+ flux), however, only changes once the temperature has reached a particular threshold, approximately 5°C below present day temperatures (corresponding to a δD value of −420‰). This reflects to a large extent the limited sensitivity of the sea ice proxy during very cold periods with large sea ice extent. At terminations where this threshold is not reached (TVI, TVIII), ssNa+ flux shows no changes. Above this threshold, the sea ice proxy is closely coupled to the Antarctic temperature, and interglacial levels are reached at the same time for both ssNa+ and δD. On the other hand, once another threshold at approximately 2°C below present day temperature is passed (corresponding to a δD value of −402‰), nssCa2+ flux has reached interglacial levels and does not change any more, despite further warming. This threshold behaviour most likely results from a combination of changes to the threshold friction velocity for dust entrainment and to the distribution of surface wind speeds in the dust source region.