Search

The European Beyond EPICA project aims to extract a continuous ice core of up to 1.5 Ma, with a maximum age density of 20 kyr m−1 at Little Dome C (LDC). We present a 1D numerical model which calculates the age of the ice around Dome C. The model inverts for basal conditions and accounts either for melting or for a layer of stagnant ice above the bedrock. It is constrained by internal reflecting horizons traced in radargrams and dated using the EPICA Dome C (EDC) ice core age profile. We used three different radar datasets ranging from a 10 000 km2 airborne survey down to 5 km long ground-based radar transects over LDC. We find that stagnant ice exists in many places, including above the LDC relief where the new Beyond EPICA drill site (BELDC) is located. The modelled thickness of this layer of stagnant ice roughly corresponds to the thickness of the basal unit observed in one of the radar surveys and in the autonomous phase-sensitive radio-echo sounder (ApRES) dataset. At BELDC, the modelled stagnant ice thickness is 198±44 m and the modelled oldest age of ice is 1.45±0.16 Ma at a depth of 2494±30 m. This is very similar to all sites situated on the LDC relief, including that of the Million Year Ice Core project being conducted by the Australian Antarctic Division. The model was also applied to radar data in the area 10–15 km north of EDC (North Patch), where we find either a thin layer of stagnant ice (generally m) or a negligible melt rate ( mm yr−1). The modelled maximum age at North Patch is over 2 Ma in most places, with ice at 1.5 Ma having a resolution of 9–12 kyr m−1, making it an exciting prospect for a future Oldest Ice drill site.

Megadunes are peculiar features formed by the interaction between the atmosphere and cryosphere and are known to be present only on the East Antarctic Plateau and other planets (Mars and Pluto). In this study, we have analysed the glaciological dynamic of megadunes, their spectral properties and morphology on two sample areas of the East Antarctic Plateau where in the past international field activities were carried out (EAIIST, East Antarctic International Ice Sheet Traverse; It-ITASE, Italian International Trans-Antarctic Scientific Expedition). Using satellite images spanning 7 years, we analysed the spatial and temporal variability in megadune surface characteristics, i.e. near-infrared (NIR) albedo, thermal brightness temperature (BT) and slope along the prevailing wind direction (SPWD), useful for mapping them. These parameters allowed us to characterize and perform an automated detection of the glazed surfaces, and we determined the influence of the SPWD by evaluating different combinations of these parameters. The inclusion of the SPWD significantly increased the accuracy of the method, doubling it in certain analysed scenes. Using remote and field observations, for the first time we surveyed all the components of upwind migration (absolute, sedimentological and ice flow), finding an absolute value of about 10 m a−1. The analysis shows that the migration is driven by the snow accumulation on the crest and trough prograding upwind on the previous windward flanks characterized by glazed surface. Our results present significant implications for the surface mass balance estimation, paleo-climate reconstruction using ice cores, and the measurements using optical and radar images/data in the megadune areas.

An airborne microwave wide-band radiometer (500–2000 MHz) was operated for the first time in Antarctica to better understand the emission properties of sea ice, outlet glaciers and the interior ice sheet from Terra Nova Bay to Dome C. The different glaciological regimes were revealed to exhibit unique spectral signatures in this portion of the microwave spectrum. Generally, the brightness temperatures over a vertically homogeneous ice sheet are warmest at the lowest frequencies, consistent with models that predict that those channels sensed the deeper, warmer parts of the ice sheet. Vertical heterogeneities in the ice property profiles can alter this basic interpretation of the signal. Spectra along the lengths of outlet glaciers were modulated by the deposition and erosion of snow, driven by strong katabatic winds. Similar to previous experiments in Greenland, the brightness temperatures across the frequency band were low in crevasse areas. Variations in brightness temperature were consistent with spatial changes in sea ice type identified in satellite imagery and in situ ground-penetrating radar data. The results contribute to a better understanding of the utility of microwave wide-band radiometry for cryospheric studies and also advance knowledge of the important physics underlying existing L-band radiometers operating in space.

Thanks to its insolubility, mineral dust is considered a stable proxy in polar ice cores. With this study we show that the Talos Dome ice core (TALDICE, Ross Sea sector of East Antarctica) displays evident and progressive signs of post-depositional processes affecting the mineral dust record below 1000 m deep. We apply a suite of established and cutting-edge techniques to investigate the properties of dust in TALDICE, ranging from concentration and grain size to elemental composition and Fe mineralogy. Results show that through acidic/oxidative weathering, the conditions of deep ice at Talos Dome promote the dissolution of specific minerals and the englacial formation of others, affecting primitive dust features. The expulsion of acidic atmospheric species from ice grains and their concentration in localized environments is likely the main process responsible for englacial reactions. Deep ice can be seen as a “geochemical reactor” capable of fostering complex reactions which involve both soluble and insoluble impurities. Fe-bearing minerals can efficiently help in exploring such transformations.

The area near Dome C, East Antarctica, is thought to be one of the most promising targets for recovering a continuous ice-core record spanning more than a million years. The European Beyond EPICA consortium has selected Little Dome C (LDC), an area ∼ 35 km southeast of Concordia Station, to attempt to recover such a record. Here, we present the results of the final ice-penetrating radar survey used to refine the exact drill site. These data were acquired during the 2019–2020 austral summer using a new, multi-channel high-resolution very high frequency (VHF) radar operating in the frequency range of 170–230 MHz. This new instrument is able to detect reflectors in the near-basal region, where previous surveys were largely unable to detect horizons. The radar stratigraphy is used to transfer the timescale of the EPICA Dome C ice core (EDC) to the area of Little Dome C, using radar isochrones dating back past 600 ka. We use these data to derive the expected depth–age relationship through the ice column at the now-chosen drill site, termed BELDC (Beyond EPICA LDC). These new data indicate that the ice at BELDC is considerably older than that at EDC at the same depth and that there is about 375 m of ice older than 600 kyr at BELDC. Stratigraphy is well preserved to 2565 m, ∼ 93 % of the ice thickness, below which there is a basal unit with unknown properties. An ice-flow model tuned to the isochrones suggests ages likely reach 1.5 Myr near 2500 m, ∼ 65 m above the basal unit and ∼ 265 m above the bed, with sufficient resolution (19 ± 2 kyr m−1) to resolve 41 kyr glacial cycles.

Results from ground-penetrating radar (GPR) measurements and shallow ice cores carried out during a scientific traverse between Dome Concordia (DC) and Vostok stations are presented in order to infer both spatial and temporal characteristics of snow accumulation over the East Antarctic Plateau. Spatially continuous accumulation rates along the traverse are computed from the identification of three equally spaced radar reflections spanning about the last 600 years. Accurate dating of these internal reflection horizons (IRHs) is obtained from a depth–age relationship derived from volcanic horizons and bomb testing fallouts on a DC ice core and shows a very good consistency when tested against extra ice cores drilled along the radar profile. Accumulation rates are then inferred by accounting for density profiles down to each IRH. For the latter purpose, a careful error analysis showed that using a single and more accurate density profile along a DC core provided more reliable results than trying to include the potential spatial variability in density from extra (but less accurate) ice cores distributed along the profile. The most striking feature is an accumulation pattern that remains constant through time with persistent gradients such as a marked decrease from 26 mm w.e. yr−1 at DC to 20 mm w.e. yr−1 at the south-west end of the profile over the last 234 years on average (with a similar decrease from 25 to 19 mm w.e. yr−1 over the last 592 years). As for the time dependency, despite an overall consistency with similar measurements carried out along the main East Antarctic divides, interpreting possible trends remains difficult. Indeed, error bars in our measurements are still too large to unambiguously infer an apparent time increase in accumulation rate. For the proposed absolute values, maximum margins of error are in the range 4 mm w.e. yr−1 (last 234 years) to 2 mm w.e. yr−1 (last 592 years), a decrease with depth mainly resulting from the time-averaging when computing accumulation rates.

We reconstruct the pattern of surface accumulation in the region around Dome C, East Antarctica, since the last glacial. We use a set of 18 isochrones spanning all observable depths of the ice column, interpreted from various ice-penetrating radar surveys and a 1-D ice flow model to invert for accumulation rates in the region. The shallowest four isochrones are then used to calculate paleoaccumulation rates between isochrone pairs using a 1-D assumption where horizontal advection is negligible in the time interval of each layer. We observe that the large-scale (100s km) surface accumulation gradient is spatially stable through the last 73 kyr, which reflects current modeled and observed precipitation gradients in the region. We also observe small-scale (10 s km) accumulation variations linked to snow redistribution at the surface, due to changes in its slope and curvature in the prevailing wind direction that remain spatially stationary since the last glacial.

Basal melt rate is the most important physical quantity to be evaluated when looking for an old-ice drilling site, and it depends to a great extent on the geothermal flux (GF), which is poorly known under the East Antarctic ice sheet. Given that wet bedrock has higher reflectivity than dry bedrock, the wetness of the ice–bed interface can be assessed using radar echoes from the bedrock. But, since basal conditions depend on heat transfer forced by climate but lagged by the thick ice, the basal ice may currently be frozen whereas in the past it was generally melting. For that reason, the risk of bias between present and past conditions has to be evaluated. The objective of this study is to assess which locations in the Dome C area could have been protected from basal melting at any time in the past, which requires evaluating GF. We used an inverse approach to retrieve GF from radar-inferred distribution of wet and dry beds. A 1-D heat model is run over the last 800 ka to constrain the value of GF by assessing a critical ice thickness, i.e. the minimum ice thickness that would allow the present local distribution of basal melting. A regional map of the GF was then inferred over a 80 km × 130 km area, with a N–S gradient and with values ranging from 48 to 60 mW m−2. The forward model was then emulated by a polynomial function to compute a time-averaged value of the spatially variable basal melt rate over the region. Three main subregions appear to be free of basal melting, two because of a thin overlying ice and one, north of Dome C, because of a low GF.

A high-resolution (1 km line spacing) aerogeophysical survey was conducted over a region near the East Antarctic Ice Sheet's Dome C that may hold a 1.5 Myr climate record. We combined new ice thickness data derived from an airborne coherent radar sounder with unpublished data that was in part unavailable for earlier compilations, and we were able to remove older data with high positional uncertainties. We generated a revised high-resolution digital elevation model (DEM) to investigate the potential for an old ice record in this region, and used laser altimetry to confirm a Cryosat-2 derived DEM for inferring the glaciological state of the candidate area. By measuring the specularity content of the bed, we were able to find an additional 50 subglacial lakes near the candidate site, and by Doppler focusing the radar data, we were able to map out the roughness of the bed at length scales of hundreds of meters. We find that the primary candidate region contains elevated rough topography interspersed with scattered subglacial lakes and some regions of smoother bed. Free subglacial water appears to be restricted from bed overlain by ice thicknesses of less than 3000 m. A site near the ice divide was selected for further investigation. The high resolution of this ice thickness data set also allows us to explore the nature of ice thickness uncertainties in the context of radar geometry and processing.

Halogen chemistry in the polar regions occurs through the release of halogen elements from different sources. Bromine is primarily emitted from sea salt aerosols and other saline condensed phases associated with sea ice surfaces, while iodine is affected by the release of organic compounds from algae colonies living within the sea ice environment. Measurements of halogen species in polar snow samples are limited to a few sites although there is some evidence that they are related to sea ice extent. We examine here total bromine, iodine and sodium concentrations in a series of 2 m cores collected during a traverse from Talos Dome (72°48' S, 159°06' E) to GV7 (70°41' S, 158°51' E) analyzed by inductively coupled plasma-sector field mass spectrometry (ICP-SFMS) at a resolution of 5 cm. We find a distinct seasonality of the bromine enrichment signal in most of the cores, with maxima during the austral spring. Iodine shows average concentrations of 0.04 ppb with little variability. No distinct seasonality is found for iodine and sodium. The transect reveals homogeneous air-to-snow fluxes for the three chemical species along the transect due to competing effects of air masses originating from the Ross Sea and the Southern Ocean.

In 2011 four ice cores were extracted from the summit of Alto dell'Ortles (3859 m), the highest glacier of South Tyrol in the Italian Alps. This drilling site is located only 37 km southwest from where the Tyrolean Iceman, ∼ 5.3 kyrs old, was discovered emerging from the ablating ice field of Tisenjoch (3210 m, near the Italian–Austrian border) in 1991. The excellent preservation of this mummy suggested that the Tyrolean Iceman was continuously embedded in prehistoric ice and that additional ancient ice was likely preserved elsewhere in South Tyrol. Dating of the ice cores from Alto dell'Ortles based on 210Pb, tritium, beta activity and 14C determinations, combined with an empirical model (COPRA), provides evidence for a chronologically ordered ice stratigraphy from the modern glacier surface down to the bottom ice layers with an age of ∼ 7 kyrs, which confirms the hypothesis. Our results indicate that the drilling site has continuously been glaciated on frozen bedrock since ∼ 7 kyrs BP. Absence of older ice on the highest glacier of South Tyrol is consistent with the removal of basal ice from bedrock during the Northern Hemisphere Climatic Optimum (6–9 kyrs BP), the warmest interval in the European Alps during the Holocene. Borehole inclinometric measurements of the current glacier flow combined with surface ground penetration radar (GPR) measurements indicate that, due to the sustained atmospheric warming since the 1980s, an acceleration of the glacier Alto dell'Ortles flow has just recently begun. Given the stratigraphic–chronological continuity of the Mt. Ortles cores over millennia, it can be argued that this behaviour has been unprecedented at this location since the Northern Hemisphere Climatic Optimum.

Global climate models suggest that Antarctic snowfall should increase in a warming climate and mitigate rises in the sea level. Several processes affect surface mass balance (SMB), introducing large uncertainties in past, present and future ice sheet mass balance. To provide an extended perspective on the past SMB of Antarctica, we used 67 firn/ice core records to reconstruct the temporal variability in the SMB over the past 800 yr and, in greater detail, over the last 200 yr. Our SMB reconstructions indicate that the SMB changes over most of Antarctica are statistically negligible and that the current SMB is not exceptionally high compared to the last 800 yr. High-accumulation periods have occurred in the past, specifically during the 1370s and 1610s. However, a clear increase in accumulation of more than 10% has occurred in high SMB coastal regions and over the highest part of the East Antarctic ice divide since the 1960s. To explain the differences in behaviour between the coastal/ice divide sites and the rest of Antarctica, we suggest that a higher frequency of blocking anticyclones increases the precipitation at coastal sites, leading to the advection of moist air in the highest areas, whereas blowing snow and/or erosion have significant negative impacts on the SMB at windy sites. Eight hundred years of stacked records of the SMB mimic the total solar irradiance during the 13th and 18th centuries. The link between those two variables is probably indirect and linked to a teleconnection in atmospheric circulation that forces complex feedback between the tropical Pacific and Antarctica via the generation and propagation of a large-scale atmospheric wave train.

Ice thickness data over much of East Antarctica are sparse and irregularly distributed. This poses difficulties for reconstructing the homogeneous coverage needed to properly assess underlying sub-glacial morphology and fundamental geometric constraints on sea level rise. Here we introduce a new physically-based ice thickness interpolation scheme and apply this to existing ice thickness data in the Aurora Subglacial Basin region. The skill and robustness of the new reconstruction is demonstrated by comparison with new data from the ICECAP project. The interpolated morphology shows an extensive marine-based ice sheet, with considerably more area below sea-level than shown by prior studies. It also shows deep features connecting the coastal grounding zone with the deepest regions in the interior. This has implications for ice sheet response to a warming ocean and underscores the importance of obtaining additional high resolution data in these marginal zones for modelling ice sheet evolution.