How warm was Greenland during the last interglacial period?

The last interglacial period (LIG, ~ 129–116 thousand years ago) provides the most recent case study for multi-millennial polar warming above pre-industrial level and a respective response of the Greenland and Antarctic ice sheets to this warming, as well as a test bed for climate and ice sheet models. Past changes in Greenland ice sheet thickness and surface temperature during this period were recently derived from the NEEM ice core records, North-West Greenland. The NEEM paradox has emerged from an estimated large local warming above pre-industrial level (7.5 ± 1.8 °C at the deposition site 126 ka ago without correction for any overall ice sheet altitude changes between the LIG and pre-industrial) based on water isotopes, together with limited local ice thinning, suggesting more resilience of the real Greenland ice sheet than shown in some ice sheet models. Here, we provide an independent assessment of the average LIG Greenland surface warming using ice core air isotopic composition (δ15N) and relationships between accumulation rate and temperature. The LIG surface temperature at the upstream NEEM deposition site without ice sheet altitude correction is estimated to be warmer by +7 to +11 °C (+8 °C being the most likely estimate according to constraints on past accumulation rate) compared to the pre-industrial period. This temperature estimate is consistent with the 7.5 ± 1.8 °C warming initially determined from NEEM water isotopes. Moreover, we show that under such warm temperatures, melting of snow probably led to a significant firn shrinking by ~ 15 m. Climate simulations performed with present day ice sheet topography lead to much smaller warming but larger amplitudes (up to 5 °C) can be obtained from changes in sea ice extent and ice sheet topography. Still, ice sheet simulations forced by 5 °C surface warming lead to large ice sheet decay that are not compatible with existing data. Our new, independent temperature constrain therefore reinforces the NEEM paradox

Assessing the robustness of Antarctic temperature reconstructions over the past 2 millennia using pseudoproxy and data assimilation experiments

The Antarctic temperature changes over the past millennia remain more uncertain than in many other continental regions. This has several origins: (1) the number of high-resolution ice cores is small, in particular on the East Antarctic plateau and in some coastal areas in East Antarctica; (2) the short and spatially sparse instrumental records limit the calibration period for reconstructions and the assessment of the methodologies; (3) the link between isotope records from ice cores and local climate is usually complex and dependent on the spatial scales and timescales investigated. Here, we use climate model results, pseudoproxy experiments and data assimilation experiments to assess the potential for reconstructing the Antarctic temperature over the last 2 millennia based on a new database of stable oxygen isotopes in ice cores compiled in the framework of Antarctica2k (Stenni et al.,). The well-known covariance between δ 18 O and temperature is reproduced in the two isotope-enabled models used (ECHAM5/MPI-OM and ECHAM5-wiso), but is generally weak over the different Antarctic regions, limiting the skill of the reconstructions. Furthermore, the strength of the link displays large variations over the past millennium, further affecting the potential skill of temperature reconstructions based on statistical methods which rely on the assumption that the last decades are a good estimate for longer temperature reconstructions. Using a data assimilation technique allows, in theory, for changes in the δ 18 O-temperature link through time and space to be taken into account. Pseudoproxy experiments confirm the benefits of using data assimilation methods instead of statistical methods that provide reconstructions with unrealistic variances in some Antarctic subregions. They also confirm that the relatively weak link between both variables leads to a limited potential for reconstructing temperature based on δ 18 O. However, the reconstruction skill is higher and more uniform among reconstruction methods when the reconstruction target is the Antarctic as a whole rather than smaller Antarctic subregions. This consistency between the methods at the large scale is also observed when reconstructing temperature based on the real δ 18 O regional composites of Stenni et al. (2017). In this case, temperature reconstructions based on data assimilation confirm the long-term cooling over Antarctica during the last millennium, and the later onset of anthropogenic warming compared with the simulations without data assimilation, which is especially visible in West Antarctica. Data assimilation also allows for models and direct observations to be reconciled by reproducing the east-west contrast in the recent temperature trends. This recent warming pattern is likely mostly driven by internal variability given the large spread of individual Paleoclimate Modelling Intercomparison Project (PMIP)/Coupled Model Intercomparison Project (CMIP) model realizations in simulating it. As in the pseudoproxy framework, the reconstruction methods perform differently at the subregional scale, especially in terms of the variance of the time series produced. While the potential benefits of using a data assimilation method instead of a statistical method have been highlighted in a pseudoproxy framework, the instrumental series are too short to confirm this in a realistic setup

Surface formation, preservation, and history of low-porosity crusts at the WAIS Divide site, West Antarctica

International audienceObservations at the West Antarctic Ice Sheet (WAIS) Divide site show that near-surface snow is strongly altered by weather-related processes such as strong winds and temperature fluctuations, producing features that are recognizable in the deep ice core. Prominent glazed surface crusts develop frequently at the site during summer seasons. Surface, snow pit, and ice core observations made in this study during summer field seasons from 2008–2009 to 2012–2013, supplemented by automated weather station (AWS) data with short- and longwave radiation sensors, revealed that such crusts formed during relatively low-wind, low-humidity, clear-sky periods with intense daytime sunshine. After formation, such glazed surfaces typically developed cracks in a polygonal pattern likely from thermal contraction at night. Cracking was commonest when several clear days occurred in succession and was generally followed by surface hoar growth; vapor escaping through the cracks during sunny days may have contributed to the high humidity that favored nighttime formation of surface hoar. Temperature and radiation observations show that daytime solar heating often warmed the near-surface snow above the air temperature, contributing to upward mass transfer, favoring crust formation from below, and then surface hoar formation. A simple surface energy calculation supports this observation. Subsequent examination of the WDC06A deep ice core revealed that crusts are preserved through the bubbly ice, and some occur in snow accumulated during winters, although not as commonly as in summertime deposits. Although no one has been on site to observe crust formation during winter, it may be favored by greater wintertime wind packing from stronger peak winds, high temperatures and steep temperature gradients from rapid midwinter warmings reaching as high as −15 °C, and perhaps longer intervals of surface stability. Time variations in crust occurrence in the core may provide paleoclimatic information, although additional studies are required. Discontinuity and cracking of crusts likely explain why crusts do not produce significant anomalies in other paleoclimatic records

Understanding the drivers of near-surface winds in Adélie Land, East Antarctica

International audienceAbstract. Near-surface winds play a crucial role in the climate of Antarctica, but accurately quantifying and understanding their drivers is complex. They result from the contribution of two distinct families of drivers: the large-scale pressure gradient and surface-induced pressure gradients known as katabatic and thermal wind. The extrapolation of vertical potential temperature above the boundary layer down to the surface enables us to separate and quantify the contribution of these different pressure gradients in the momentum budget equations. Using this method applied to outputs of the regional atmospheric model MAR at a 3-hourly resolution, we find that the seasonal and spatial variability in near-surface winds in Adélie Land is dominated by surface processes. On the other hand, high-frequency temporal variability (3-hourly) is mainly controlled by large-scale variability everywhere in Antarctica, except on the coast. In coastal regions, although the katabatic acceleration surpasses all other accelerations in magnitude, none of the katabatic or large-scale accelerations can be identified as the single primary driver of near-surface wind variability. The angle between the large-scale acceleration and the surface slope is a key factor in explaining strong wind speed events: the highest-wind-speed events happen when the katabatic and large-scale forcing are aligned, although each acceleration, when acting alone, can also cause strong wind speed. This study underlines the complexity of the drivers of Antarctic surface winds and the value of the momentum budget decomposition to identify drivers at different spatial and temporal scales

Control of the temperature signal in Antarctic proxies by snowfall dynamics

International audienceAbstract. Antarctica, the coldest and driest continent, is home to the largest ice sheet, whose mass is predominantly recharged by snowfall. A common feature of polar regions is the warming associated with snowfall, as moist oceanic air and cloud cover increase the surface temperature. Consequently, snow that accumulates on the ice sheet is deposited under unusually warm conditions. Here we use a polar-oriented regional atmospheric model to study the statistical difference between average and snowfall-weighted temperatures. During snowfall, the warm anomaly scales with snowfall amount, with the strongest sensitivity occurring at low-accumulation sites. Heavier snowfall in winter helps to decrease the annual snowfall-weighted temperature, but this effect is overwritten by the event-scale warming associated with precipitating atmospheric systems, which particularly contrast with the extremely cold conditions that occur in winter. Consequently, the seasonal range of snowfall-weighted temperature is reduced by 20 %. On the other hand, the annual snowfall-weighted temperature shows 80 % more interannual variability than the annual temperature due to the irregularity of snowfall occurrence and its associated temperature anomaly. Disturbances of the apparent annual temperature cycle and interannual variability have important consequences for the interpretation of water isotopes in precipitation, which are deposited with snowfall and commonly used for paleotemperature reconstructions from ice cores

Comparison of different methods to retrieve optical-equivalent snow grain size in central Antarctica

International audienceThe optical-equivalent snow grain size affects the reflectivity of snow surfaces and, thus, the local surface energy budget in particular in polar regions. Therefore, the specific surface area (SSA), from which the optical snow grain size is derived, was observed for a 2-month period in central Antarctica (Kohnen research station) during austral summer 2013/14. The data were retrieved on the basis of ground-based spectral surface albedo measurements collected by the COmpact RAdiation measurement System (CORAS) and airborne observations with the Spectral Modular Airborne Radiation measurement sysTem (SMART). The snow grain size and pollution amount (SGSP) algorithm, originally developed to analyze spaceborne reflectance measurements by the MODerate Resolution Imaging Spectroradiometer (MODIS), was modified in order to reduce the impact of the solar zenith angle on the retrieval results and to cover measurements in overcast conditions. Spectral ratios of surface albedo at 1280 and 1100 nm wavelength were used to reduce the retrieval uncertainty. The retrieval was applied to the ground-based and airborne observations and validated against optical in situ observations of SSA utilizing an IceCube device. The SSA retrieved from CORAS observations varied between 27 and 89 m2 kg−1. Snowfall events caused distinct relative maxima of the SSA which were followed by a gradual decrease in SSA due to snow metamorphism and wind-induced transport of freshly fallen ice crystals. The ability of the modified algorithm to include measurements in overcast conditions improved the data coverage, in particular at times when precipitation events occurred and the SSA changed quickly. SSA retrieved from measurements with CORAS and MODIS agree with the in situ observations within the ranges given by the measurement uncertainties. However, SSA retrieved from the airborne SMART data slightly underestimated the ground-based results

A 2000-year temperature reconstruction on the East Antarctic plateau from argon–nitrogen and water stable isotopes in the Aurora Basin North ice core

International audienceThe temperature of the Earth is one of the most important climate parameters. Proxy records of past climate changes, in particular temperature, represent a fundamental tool for exploring internal climate processes and natural climate forcings. Despite the excellent information provided by ice core records in Antarctica, the temperature variability of the past 2000 years is difficult to evaluate from the low-accumulation sites in the Antarctic continent interior. Here we present the results from the Aurora Basin North (ABN) ice core (71∘ S, 111∘ E, 2690 m a.s.l.) in the lower part of the East Antarctic plateau, where accumulation is substantially higher than other ice core drilling sites on the plateau, and provide unprecedented insight into East Antarctic past temperature variability. We reconstructed the temperature of the last 2000 years using two independent methods: the widely used water stable isotopes (δ18O) and by inverse modelling of borehole temperature and past temperature gradients estimated from the inert gas stable isotopes (δ40Ar and δ15N). This second reconstruction is based on three independent measurement types: borehole temperature, firn thickness, and firn temperature gradient. The δ18O temperature reconstruction supports stable temperature conditions within 1 ∘C over the past 2000 years, in agreement with other ice core δ18O records in the region. However, the gas and borehole temperature reconstruction suggests that surface conditions 2 ∘C cooler than average prevailed in the 1000–1400 CE period and supports a 20th century warming of 1 ∘C. A precipitation hiatus during cold periods could explain why water isotope temperature reconstruction underestimates the temperature changes. Both reconstructions arguably record climate in their own way, with a focus on atmospheric and hydrologic cycles for water isotopes, as opposed to surface temperature for gas isotopes and boreholes. This study demonstrates the importance of using a variety of sources for comprehensive paleoclimate reconstructions

Experimental determination and theoretical framework of kinetic fractionation at the water vapour–ice interface at low temperature

International audienc