Snow height on sea ice and sea ice drift from autonomous measurements from buoy 2014S24, deployed during Antarctic Fast Ice Network 2014 (AFIN 2014)

Snow height was measured by the Snow Depth Buoy 2014S24, an autonomous platform, installed close to Neumayer III Base, Antarctic during Antarctic Fast Ice Network 2014 (AFIN 2014). The resulting time series describes the evolution of snow depth as a function of place and time between 2014-03-07 and 2014-05-16 in sample intervals of 1 hour. The Snow Depth Buoy consists of four independent sonar measurements representing the area (approx. 10 m**2) around the buoy. The buoy was installed on the ice shelf. In addition to snow depth, geographic position (GPS), barometric pressure, air temperature, and ice surface temperature were measured. Negative values of snow depth occur if surface ablation continues into the sea ice. Thus, these measurements describe the position of the sea ice surface relative to the original snow-ice interface. Differences between single sensors indicate small-scale variability of the snow pack around the buoy. The data set has been processed, including the removal of obvious inconsistencies (missing values). Records without any snow depth may still be used for sea ice drift analyses. Note: This data set contains only relative changes in snow depth, because no initial readings of absolute snow depth are available.

Ice and snow measurements during Nathaniel B. Palmer cruise NBP0709 (SIMBA experiment)

The Sea Ice Mass Balance in the Antarctic (SIMBA) experiment was conducted from the RVIB N.B. Palmer in September and October 2007 in the Bellingshausen Sea in an area recently experiencing considerable changes in both climate and sea ice cover. Snow and ice properties were observed at 3 short-term stations and a 27-day drift station (Ice Station Belgica, ISB) during the winter-spring transition. Repeat measurements were performed on sea ice and snow cover at 5 ISB sites, each having different physical characteristics, with mean ice (snow) thicknesses varying from 0.6 m (0.1 m) to 2.3 m (0.7 m). Ice cores retrieved every five days from 2 sites and measured for physical, biological, and chemical properties. Three ice mass-balance buoys (IMBs) provided continuous records of snow and ice thickness and temperature. Meteorological conditions changed from warm fronts with high winds and precipitation followed by cold and calm periods through four cycles during ISB. The snow cover regulated temperature flux and controlled the physical regime in which sea ice morphology changed. Level thin ice areas had little snow accumulation and experienced greater thermal fluctuations resulting in brine salinity and volume changes, and winter maximum thermodynamic growth of ~0.6 m in this region. Flooding and snow-ice formation occurred during cold spells in ice and snow of intermediate thickness. In contrast, little snow-ice formed in flooded areas with thicker ice and snow cover, instead nearly isothermal, highly permeable ice persisted. In spring, short-lived cold air episodes did not effectively penetrate the sea ice nor overcome the effect of ocean heat flux, thus favoring net ice thinning from bottom melt over ice thickening from snow-ice growth, in all cases. These warm ice conditions were consistent with regional remote sensing observations of earlier ice breakup and a shorter sea ice season, more recently observed in the Bellingshausen Sea.

(Appendix A) Geochemical composition of snow and eolian mineral dust from Berkner Island ice sheet, Antarctica

The Sr and Nd isotopic composition of dust extracted from recent snow layers at the top of Berkner Island ice sheet (located within the Filchner-Ronne Ice Shelf at the southern end of the Weddell Sea) enables us, for the first time, to document dust provenance in Antarctica outside the East Antarctic Plateau (EAP) where all previous studies based on isotopic fingerprinting were carried out. Berkner dust displays an overall crust-like isotopic signature, characterized by more radiogenic 87Sr/86Sr and much less radiogenic 143Nd/144Nd compared to dust deposited on the EAP during glacial periods. Differences with EAP interglacial dust are not as marked but still significant, indicating that present-day Berkner dust provenance is distinct, at least to some extent, from that of the dust reaching the EAP. The fourteen snow-pit sub-seasonal samples that were obtained span a two-year period (2002-2003) and their dust Sr and Nd isotopic composition reveals that multiple sources are at play over a yearly time period. Southern South America, Patagonia in particular, likely accounts for part of the observed spring/summer dust deposition maxima, when isotopic composition is shifted towards 'younger' isotopic signatures. In the spring, possible additional inputs from Australian sources would also be supported by the data. Most of the year, however, the measured isotopic signatures would be best explained by a sustained background supply from putative local sources in East Antarctica, which carry old-crust-like isotopic fingerprints. Whether the restricted East Antarctic ice-free areas produce sufficient eolian material has yet to be substantiated however. The fact that large (> 5 µm) particles represent a significant fraction of the samples throughout the entire time-series supports scenarios that involve contributions from proximal sources, either in Patagonia and/or Antarctica (possibly including snow-free areas in the Antarctic Peninsula and other areas as well). This also indicates that additional dust transport, which does not reach the EAP, must occur at low-tropospheric levels to this coastal sector of Antarctica.

Rare earth elements determined in Antarctic ice by inductively coupled plasma

Inductively coupled plasma mass spectrometry (ICP-MS) is a suitable tool for multi-element analysis at low concentration levels. Rare earth element (REE) determinations in standard reference materials and small volumes of molten ice core samples from Antarctica have been performed with an ICP-time of flight-MS (ICP-TOF-MS) system. Recovery rates for REE in e.g. SPS-SW1 amounted to not, vert, similar ~103%, and the relative standard deviations were 3.4% for replicate analysis at REE concentrations in the lower ng/l range. Analyses of REE concentrations in Antarctic ice core samples showed that the ICP-TOF-MS technique meets the demands of restricted sample mass. The data obtained are in good agreement with ICP-Quadrupole-MS (ICP-Q-MS) and ICP-Sector Field-MS (ICP-SF-MS) results. The ICP-TOF-MS system determines accurately and precisely REE concentrations exceeding 5 ng/l while between 0.5 and 5 ng/l accuracy and precision are element dependent.

(Table 1) Sea ice and snow characteristics and heat fluxes observed during R/V Aurora Australis cruise to East Antarctica in September/October 2007

The properties of snow on East Antarctic sea ice off Wilkes Land were examined during the Sea Ice Physics and Ecosystem Experiment (SIPEX) in late winter of 2007, focusing on the interaction with sea ice. This observation includes 11 transect lines for the measurement of ice thickness, freeboard, and snow depth, 50 snow pits on 13 ice floes, and diurnal variation of surface heat flux on three ice floes. The detailed profiling of topography along the transects and the d18O, salinity, and density datasets of snow made it possible to examine the snow-sea-ice interaction quantitatively for the first time in this area. In general, the snow displayed significant heterogeneity in types, thickness (mean: 0.14 +- 0.13 m), and density (325 +- 38 kg/m**3), as reported in other East Antarctic regions. High salinity was confined to the lowest 0.1 m. Salinity and d18O data within this layer revealed that saline water originated from the surface brine of sea ice in 20% of the total sites and from seawater in 80%. From the vertical profiles of snow density, bulk thermal conductivity of snow was estimated as 0.15 W/K/m on average, only half of the value used for numerical sea-ice models. Although the upward heat flux within snow estimated with this value was significantly lower than that within ice, it turned out that a higher value of thermal conductivity (0.3 to 0.4 W/K/m) is preferable for estimating ice growth amount in current numerical models. Diurnal measurements showed that upward conductive heat flux within the snow and net long-wave radiation at the surface seem to play important roles in the formation of snow ice from slush. The detailed surface topography allowed us to compare the air-ice drag coefficients of ice and snow surfaces under neutral conditions, and to examine the possibility of the retrieval of ice thickness distribution from satellite remote sensing. It was found that overall snow cover works to enhance the surface roughness of sea ice rather than moderate it, and increases the drag coefficient by about 10%. As for thickness retrieval, mean ice thickness had a higher correlation with ice surface roughness than mean freeboard or surface elevation, which indicates the potential usefulness of satellite L-band SAR in estimating the ice thickness distribution in the seasonal sea-ice zone.

Snow-pit measurements on Atka Bay landfast sea ice during ANT-Land 2012/2013 season

Up to now, snow cover on Antarctic sea ice and its impact on radar backscatter, particularly after the onset of freeze/thaw processes, are not well understood. Here we present a combined analysis of in situ observations of snow properties from the landfast sea ice in Atka Bay, Antarctica, and high-resolution TerraSAR-X backscatter data, for the transition from austral spring (November 2012) to summer (January 2013). The physical changes in the seasonal snow cover during that time are reflected in the evolution of TerraSAR-X backscatter. We are able to explain 76-93% of the spatio-temporal variability of the TerraSAR-X backscatter signal with up to four snowpack parameters with a root-mean-squared error of 0.87-1.62 dB, using a simple multiple linear model. Over the complete study, and especially after the onset of early-melt processes and freeze/thaw cycles, the majority of variability in the backscatter is influenced by changes in snow/ice interface temperature, snow depth and top-layer grain size. This suggests it may be possible to retrieve snow physical properties over Antarctic sea ice from X-band SAR backscatter.

(Table 3) Details of field-based ice and snow profile measurements from the SIMBA 2007 experiment

Although sea-ice extent in the Bellingshausen-Amundsen (BA) seas sector of the Antarctic has shown significant decline over several decades, there is not enough data to draw any conclusion on sea-ice thickness and its change for the BA sector, or for the entire Southern Ocean. This paper presents our results of snow and ice thickness distributions from the SIMBA 2007 experiment in the Bellingshausen Sea, using four different methods (ASPeCt ship observations, downward-looking camera imaging, ship-based electromagnetic induction (EM) sounding, and in situ measurements using ice drills). A snow freeboard and ice thickness model generated from in situ measurements was then applied to contemporaneous ICESat (satellite laser altimetry) measured freeboard to derive ice thickness at the ICESat footprint scale. Errors from in situ measurements and from ICESat freeboard estimations were incorporated into the model, so a thorough evaluation of the model and uncertainty of the ice thickness estimation from ICESat are possible. Our results indicate that ICESat derived snow freeboard and ice thickness distributions (asymmetrical unimodal tailing to right) for first-year ice (0.29 ± 0.14 m for mean snow freeboard and 1.06 ± 0.40 m for mean ice thickness), multi-year ice (0.48 ± 0.26 and 1.59 ± 0.75 m, respectively), and all ice together (0.42 ± 0.24 and 1.38 ± 0.70 m, respectively) for the study area seem reasonable compared with those values from the in situ measurements, ASPeCt observations, and EM measurements. The EM measurements can act as an appropriate supplement for ASPeCt observations taken hourly from the ship's bridge and provide reasonable ice and snow distributions under homogeneous ice conditions. Our proposed approaches: (1) of using empirical equations relating snow freeboard to ice thickness based on in situ measurements and (2) of using isostatic equations that replace snow depth with snow freeboard (or empirical equations that convert freeboard to snow depth), are efficient and important ways to derive ice thickness from ICESat altimetry at the footprint scale for Antarctic sea ice. Spatial and temporal snow and ice thickness from satellite altimetry for the BA sector and for the entire Southern Ocean is therefore possible.

Lagrangian back-tracing results for Kohnen station; output of time dependent numerical experiment with a 3D nested Antarctic ice sheet model

A nested ice flow model was developed for eastern Dronning Maud Land to assist with the dating and interpretation of the EDML deep ice core. The model consists of a high-resolution higher-order ice dynamic flow model that was nested into a comprehensive 3-D thermomechanical model of the whole Antarctic ice sheet. As the drill site is on a flank position the calculations specifically take into account the effects of horizontal advection as deeper ice in the core originated from higher inland. First the regional velocity field and ice sheet geometry is obtained from a forward experiment over the last 8 glacial cycles. The result is subsequently employed in a Lagrangian backtracing algorithm to provide particle paths back to their time and place of deposition. The procedure directly yields the depth-age distribution, surface conditions at particle origin, and a suite of relevant parameters such as initial annual layer thickness. This paper discusses the method and the main results of the experiment, including the ice core chronology, the non-climatic corrections needed to extract the climatic part of the signal, and the thinning function. The focus is on the upper 89% of the ice core (appr. 170 kyears) as the dating below that is increasingly less robust owing to the unknown value of the geothermal heat flux. It is found that the temperature biases resulting from variations of surface elevation are up to half of the magnitude of the climatic changes themselves.

Snow depths of volcanic events and accumulation rates of firn cores from the IPY traverse, Antarctica

Volcanic signatures in ice-core records provide an excellent means to date the cores and obtain information about accumulation rates. From several ice cores it is thus possible to extract a spatio-temporal accumulation pattern. We show records of electrical conductivity and sulfur from 13 firn cores from the Norwegian-USA scientific traverse during the International Polar Year 2007-2009 (IPY) through East Antarctica. Major volcanic eruptions are identified and used to assess century-scale accumulation changes. The largest changes seem to occur in the most recent decades with accumulation over the period 1963-2007/08 being up to 25% different from the long-term record. There is no clear overall trend, some sites show an increase in accumulation over the period 1963 to present while others show a decrease. Almost all of the sites above 3200 m above sea level (asl) suggest a decrease. These sites also show a significantly lower accumulation value than large-scale assessments both for the period 1963 to present and for the long-term mean at the respective drill sites. The spatial accumulation distribution is influenced mainly by elevation and distance to the ocean (continentality), as expected. Ground-penetrating radar data around the drill sites show a spatial variability within 10-20% over several tens of kilometers, indicating that our drill sites are well representative for the area around them. Our results are important for large-scale assessments of Antarctic mass balance and model validation.

Atmospheric carbon dioxide, methane, deuterium, and calculated Antarctic temperature of EPICA Dome C ice core

Understanding the role of atmospheric CO2 during past climate changes requires clear knowledge of how it varies in time relative to temperature. Antarctic ice cores preserve highly resolved records of atmospheric CO2 and Antarctic temperature for the past 800,000 years. Here we propose a revised relative age scale for the concentration of atmospheric CO2 and Antarctic temperature for the last deglacial warming, using data from five Antarctic ice cores. We infer the phasing between CO2 concentration and Antarctic temperature at four times when their trends change abruptly. We find no significant asynchrony between them, indicating that Antarctic temperature did not begin to rise hundreds of years before the concentration of atmospheric CO2, as has been suggested by earlier studies.

Measurements of stable carbon isotope ratios of CO2 over the last 24000 years and CO2 concentration measurements on Antarctic ice cores using three different methods

The stable carbon isotope ratio of atmospheric CO2 (d13Catm) is a key parameter in deciphering past carbon cycle changes. Here we present d13Catm data for the past 24,000 years derived from three independent records from two Antarctic ice cores. We conclude that a pronounced 0.3 per mil decrease in d13Catm during the early deglaciation can be best explained by upwelling of old, carbon-enriched waters in the Southern Ocean. Later in the deglaciation, regrowth of the terrestrial biosphere, changes in sea surface temperature, and ocean circulation governed the d13Catm evolution. During the Last Glacial Maximum, d13Catm and atmospheric CO2 concentration were essentially constant, which suggests that the carbon cycle was in dynamic equilibrium and that the net transfer of carbon to the deep ocean had occurred before then.