ALBERT

Description: A new 21.3m firn core was drilled in 2015 at a coastal Antarctic high accumulation site in Adélie Land (66.78°S; 139.56°E, 602ma.s.l.). The core was dated by annual layers counting based on non-sea-salt sulfate and methanesulfonate summer peaks, refined by a comparison between the reconstructed surface mass balance (hereafter, SMB) and the closest available stake data. The mean reconstructed SMB of 75.2 ± 15.0cmw.e. y−1 is consistent with local stake data, and remarkably high for coastal East Antarctica. The resulting inter-annual and sub-annual variations in isotopic records (δ18O and deuterium excess, hereafter d-excess) are explored for 1998–2014 and are systematically compared with a couple of climatic time series: an updated database of Antarctic surface snow isotopic composition, SMB stake data, meteorological observations from Dumont d'Urville station, sea-ice concentration based on passive microwave satellite data, precipitation outputs of atmospheric reanalyses, climate and water stable isotope outputs from the atmospheric general circulation model ECHAM5-wiso, as well as air mass origins diagnosed using 5-days back-trajectories. The mean isotopic values (−19.3 ± 3.1‰ for δ18O and 5.4 ± 2.2‰ for d-excess) are consistent with other coastal Antarctic values. No significant isotope-temperature relationship can be evidenced at any timescale, ruling out a simple interpretation of in terms of local temperature. An observed asymmetry in the δ18O seasonal cycle may be explained by the precipitation of air masses coming from Indian and Pacific/West Antarctic Ice Sheet sectors in autumn and winter times, recorded in the d-excess signal showing outstanding values in austral spring versus autumn. Significant positive trends are observed in the annual d-excess record and local sea-ice extent (135°E–145°E) over the period 1998–2014. However, processes studies focusing on resulting isotopic compositions and particularly the d-excess-δ18O relationship, evidenced as a potential fingerprint of moisture origins, as well as the collection of more isotopic measurements in Adélie Land are needed for an accurate interpretation of our signals.

Description: The present data consists in a database of isotope (δO18, δD and deuterium excess) data from precipitations, snow and firn/ice cores, gathering the following data: - the isotopic surface snow data from Masson et al. (2008, doi:10.1175/2007JCLI2139.1) - the Antarctica2k database from Stenni et al. (2017, doi:10.5194/cp-13-1609-2017), available on https://www.ncdc.noaa.gov/paleo-search/study/22589 - the data from Fernandoy et al. (2012, doi:10.5194/tc-6-313-2012) - the precipitation data from Rozanski et al. (1993, doi:10.1029/GM078p0001) and available on the IAEA/GNIP platform -data personnally communicated See below for full references of articles and datasets. The "averages" xls file give necessary informations to retrieve the data a its original temperoral scale, as well as time-averaged, standard deviations and extremum values. They are completed isotope time-averages and standard deviations from the ECHAM5-wiso model forced to the ERA-interim reanalysis and run at the daily scale over the period 1979-2013 (Werner et al., 2011; doi:10.1029/2011JD015681). The "seasonal_snowfall.xls" file give the seasonal cycles of precipitation, temperature, δO18 and deuterium excecss of snowfall data, as used in the associated manuscript.

Description: Glaciers in the Andes are among the fastest shrinking and largest contributors to sea level rise in the world. They also represent crucial water resources in many tropical and semi-arid portions of this extensive mountain range. Yet, the magnitude of the recent ice loss over the entire Andes cordillera is still debated. Here we present Andean glacier mass changes (from 10°N to 56°S) between 2000 and 2018 using time series of digital elevation models derived from ASTER stereo images. The total mass change is -22.9 ± 5.9 Gt yr-1 (-0.72 ± 0.22 m w.e. yr-1), with the most negative mass balances in the Patagonic Andes (-0.78 ± 0.25 m w.e. yr-1) and the Tropical Andes (-0.42 ± 0.24 m w.e. yr-1), and relatively moderate losses (-0.28 ± 0.18 m w.e. yr-1) in the Dry Andes. Sub-period analysis (2000-2009 vs. 2009-2018) reveals steady mass loss in the Tropics and south of 45°S. Conversely, a shift from slightly positive to strongly negative mass balance is measured between 26° and 45°S. In this latter region, the rapid glacier loss in recent years coincides with extremely dry conditions since 2010 and partially helped to mitigate the negative hydrological impacts of this severe drought. These results provide a comprehensive, high resolution and multi-decadal dataset of recent Andes-wide glacier mass changes that constitutes a relevant basis for calibration and validation of hydrological and glaciological models intended to project future glacier changes and their hydrological impacts.

Description: Fluctuations in the length of 72 glaciers in the Northern and Southern Patagonia Icefield (NPI and SPI, respectively) and the Cordillera Darwin Icefield (CDI) were estimated between 1945 and 2005. The information obtained from historical maps based on 1945 aerial photographs was compared to ASTER and Landsat satellite images and to information found in the literature. The majority of glaciers have retreated considerably, with maximum values of 12.2 km for Marinelli Glacier in the CDI, 11.6 km for O'Higgins Glacier in the SPI and 5.7 km for San Rafael Glacier in the NPI. Among the 20 glaciers that have retreated the most relative to their size, small (less than 50 km**2) and medium (between 50 and 200 km**2) glaciers are the most affected. However, no direct relation between glacier retreat and size was found for the 72 glaciers studied. The highest percentage retreat in the CDI was by the CDI-03 Glacier (37.9%) and Marinelli Glacier (37.6%). In the SPI, relative retreats were heterogeneous and fluctuated between 27.2% (Amelia Glacier) and 0.4% (Viedma Glacier). In the NPI, relative retreat was very high for Strindberg and Cachet glaciers (35.9% and 27.6%, respectively) but for the remaining glaciers in this icefield it ranged between 11.8% (Piscis Glacier) and 3.6% (San Quintin Glacier). In addition to surface area, the surface slope (calculated on the basis of the DEM SRTM) was also related to the relative retreat and no straightforward relation was found. From a global point of view, we suggest that glacier retreat in the region is controlled firstly by atmospheric warming, as it has been reported in this area. Besides the general increase in temperature observed, no signal of a geographical pattern for the fluctuations in glacier length was found. Consequently, glaciers appear to initially react to local conditions most probably induced by their exposition, geometry and hypsometry. The heterogeneity of rates of retreat suggests that differences in basin geometry, glacier dynamics and response time are key features to explain fluctuations of each glacier.

Description: A new 21.3m firn core was drilled in 2015 at a coastal Antarctic high-accumulation site in Adélie Land (66.78◦ S; 139.56◦ E, 602 m a.s.l.), named Terre Adélie 192A (TA192A). The mean isotopic values (−19.3 ‰ ± 3.1 ‰ for δ18O and 5.4 ‰±2.2 ‰ for deuterium excess) are consistent with other coastal Antarctic values. No significant isotope–temperature relationship can be evidenced at any timescale. This rules out a simple interpretation in terms of local temperature. An observed asymmetry in the δ18O seasonal cycle may be explained by the precipitation of air masses coming from the eastern and western sectors in autumn and winter, recorded in the d-excess signal showing outstanding values in austral spring versus autumn. Significant positive trends are observed in the annual d-excess record and local sea ice extent (135–145◦ E) over the period 1998–2014. However, process studies focusing on resulting isotopic compositions and particularly the deuterium excess–δ18O relationship, evidenced as a potential fingerprint of moisture origins, as well as the collection of more isotopic measurements in Adélie Land are needed for an accurate interpretation of our signals.

Description: The Sentinel Application Platform (SNAP) architecture facilitates Earth Observation data processing. In this work, we present results from a new Snow Processor for SNAP. We also describe physical principles behind the developed snow property retrieval technique based on the analysis of Ocean and Land Colour Instrument (OLCI) onboard Sentinel-3A/B measurements over clean and polluted snow fields. Using OLCI spectral reflectance measurements in the range 400–1020 nm, we derived important snow properties such as spectral and broadband albedo, snow specific surface area, snow extent and grain size on a spatial grid of 300 m. The algorithm also incorporated cloud screening and atmospheric correction procedures over snow surfaces. We present validation results using ground measurements from Antarctica, the Greenland ice sheet and the French Alps. We find the spectral albedo retrieved with accuracy of better than 3% on average, making our retrievals sufficient for a variety of applications. Broadband albedo is retrieved with the average accuracy of about 5% over snow. Therefore, the uncertainties of satellite retrievals are close to experimental errors of ground measurements. The retrieved surface grain size shows good agreement with ground observations. Snow specific surface area observations are also consistent with our OLCI retrievals. We present snow albedo and grain size mapping over the inland ice sheet of Greenland for areas including dry snow, melted/melting snow and impurity rich bare ice. The algorithm can be applied to OLCI Sentinel-3 measurements providing an opportunity for creation of long-term snow property records essential for climate monitoring and data assimilation studies—especially in the Arctic region, where we face rapid environmental changes including reduction of snow/ice extent and, therefore, planetary albedo.

Description: Polar ice core water isotope records are commonly used to infer past changes in Antarctic temperature, motivating an improved understanding and quantification of the temporal relationship between δ18O and temperature. This can be achieved using simulations performed by atmospheric general circulation models equipped with water stable isotopes. Here, we evaluate the skills of the high-resolution water-isotope-enabled atmospheric general circulation model ECHAM5-wiso (the European Centre Hamburg Model) nudged to European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis using simulations covering the period 1960–2013 over the Antarctic continent. We compare model outputs with field data, first with a focus on regional climate variables and second on water stable isotopes, using our updated dataset of water stable isotope measurements from precipitation, snow, and firn–ice core samples. ECHAM5-wiso simulates a large increase in temperature from 1978 to 1979, possibly caused by a discontinuity in the European Reanalyses (ERA) linked to the assimilation of remote sensing data starting in 1979. Although some model–data mismatches are observed, the (precipitation minus evaporation) outputs are found to be realistic products for surface mass balance. A warm model bias over central East Antarctica and a cold model bias over coastal regions explain first-order δ18O model biases by too strong isotopic depletion on coastal areas and underestimated depletion inland. At the second order, despite these biases, ECHAM5-wiso correctly captures the observed spatial patterns of deuterium excess. The results of model–data comparisons for the inter-annual δ18O standard deviation difer when using precipitation or ice core data. Further studies should explore the importance of deposition and post-deposition processes affecting ice core signals and not resolved in the model. These results build trust in the use of ECHAM5-wiso outputs to investigate the spatial, seasonal, and inter-annual δ18O–temperature relationships. We thus make the first Antarctica-wide synthesis of prior results. First, we show that local spatial or seasonal slopes are not a correct surrogate for inter-annual temporal slopes, leading to the conclusion that the same isotope–temperature slope cannot be applied for the climatic interpretation of Antarctic ice core for all timescales. Finally, we explore the phasing between the seasonal cycles of deuterium excess and δ18O as a source of information on changes in moisture sources affecting the δ18O–temperature relationship. The few available records and ECHAM5-wiso show different phase relationships in coastal, intermediate, and central regions. This work evaluates the use of the ECHAM5-wiso model as a tool for the investigation of water stable isotopes in Antarctic precipitation and calls for extended studies to improve our understanding of such proxies.

Description: A 22.4 m-long shallow firn core was extracted during the 2006/2007 field season from coastal Adélie Land. Annual layer counting based on subannual analyses of δ18O and major chemical components was combined with 5 reference years associated with nuclear tests and non-retreat of summer sea ice to build the initial ice-core chronology (1946– 2006), stressing uncertain counting for 8 years. We focus here on the resulting δ18O and accumulation records. With an average value of 21.8 ± 6.9 cm w.e. yr−1 , local accumulation shows multi-decadal variations peaking in the 1980s, but no long-term trend. Similar results are obtained for δ18O, also characterised by a remarkably low and variable amplitude of the seasonal cycle. The ice-core records are compared with regional records of temperature, stake area accumulation measurements and variations in sea-ice extent, and outputs from two models nudged to ERA (European Reanalysis) atmospheric reanalyses: the high-resolution atmospheric general circulation model (AGCM), including stable water isotopes ECHAM5-wiso (European Centre Hamburg model), and the regional atmospheric model Modèle Atmosphérique Régional ( AR). A significant linear correlation is identified between decadal variations in δ18O and regional temperature. No significant relationship appears with regional sea-ice extent. A weak and significant correlation appears with Dumont d’Urville wind speed, increasing after 1979. The model-data comparison highlights the inadequacy of ECHAM5-wiso simulations prior to 1979, possibly due to the lack of data assimilation to constrain atmospheric reanalyses. Systematic biases are identified in the ECHAM5-wiso simulation, such as an overestimation of the mean accumulation rate and its interannual variability, a strong cold bias and an underestimation of the mean δ18O value and its interannual variability. As a result, relationships between simulated δ18O and temperature are weaker than observed. Such systematic precipitation and temperature biases are not displayed by MAR, suggesting that the model resolution plays a key role along the Antarctic ice sheet coastal topography. Interannual variations in ECHAM5-wiso temperature and precipitation accurately capture signals from meteorological data and stake observations and are used to refine the initial ice-core chronology within 2 years. After this adjustment, remarkable positive (negative) δ18O anomalies are identified in the ice-core record and the ECHAM5-wiso simulation in 1986 and 2002 (1998–1999), respectively. Despite uncertainties associated with post-deposition processes and signal-to-noise issues, in one single coastal ice-core record, we conclude that the S1C1 core can correctly capture major annual anomalies in δ18O as well as multi-decadal variations. These findings highlight the importance of improving the network of coastal high-resolution ice-core records, and stress the skills and limitations of atmospheric models for accumulation and δ18O in coastal Antarctic areas. This is particularly important for the overall East Antarctic ice sheet mass balance.

Description: The Weddell Sea sector is one of the main formation sites for Antarctic Bottom Water and an outlet for about one fifth of Antarctica’s continental ice volume. Over the last few decades, studies on glacialegeological records in this sector have provided conflicting reconstructions of changes in ice-sheet extent and ice-sheet thickness since the Last Glacial Maximum (LGM at ca 23e19 calibrated kiloyears before present, cal ka BP). Terrestrial geomorphological records and exposure ages obtained from rocks in the hinterland of the Weddell Sea, ice-sheet thickness constraints from ice cores and some radiocarbon dates on offshore sediments were interpreted to indicate no significant ice thickening and locally restricted grounding-line advance at the LGM. Other marine geological and geophysical studies concluded that subglacial bedforms mapped on theWeddell Sea continental shelf, subglacial deposits and sediments over-compacted by overriding ice recovered in cores, and the few available radiocarbon ages from marine sediments are consistent with major ice-sheet advance at the LGM. Reflecting the geological interpretations, different icesheet models have reconstructed conflicting LGM ice-sheet configurations for the Weddell Sea sector. Consequently, the estimated contributions of ice-sheet build-up in the Weddell Sea sector to the LGM sealevel low-stand of w130 m vary considerably. In this paper, we summarise and review the geological records of past ice-sheet margins and past icesheet elevations in the Weddell Sea sector. We compile marine and terrestrial chronological data constraining former ice-sheet size, thereby highlighting different levels of certainty, and present two alternative scenarios of the LGM ice-sheet configuration, including time-slice reconstructions for post- LGM grounding-line retreat. Moreover, we discuss consistencies and possible reasons for inconsistencies between the various reconstructions and propose objectives for future research. The aim of our study is to provide two alternative interpretations of glacialegeological datasets on Antarctic Ice- Sheet History for the Weddell Sea sector, which can be utilised to test and improve numerical icesheet models