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 <60 m) or a negligible melt rate (<0.1 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.

The future contributions of the Antarctic Ice Sheet to sea level rise will depend on the evolution of its surface mass balance (SMB), which could amplify/dampen mass losses increasingly observed at the ice sheet's edge. In situ constraints of SMB over annual-to-decadal timescales consist mostly of firn/ice cores that have a surface footprint cm. SMB constraints also come from climate models, which have a higher temporal resolution but a larger surface footprint of several km. We use ice-penetrating radar data to obtain an intermediate spatial and temporal resolution SMB record over three ice rises along the Princess Ragnhild Coast. The co-located ice cores allow us to obtain absolute radar-derived SMB rates at a multi-annual-to-decadal temporal resolution. By comparing the ice core SMB measurements and the radar-derived SMB records, we determine that pointwise measurements of SMB are representative of a small surface area, m radius extending from the ice core drill site for the ice rises studied here, and that the pointwise measurements are systematically 7–15 cm w.e. a lower than the mean SMB value calculated for the whole ice rises. However, ice core records are representative of an entire ice rise's temporal variability at the temporal resolution examined.

The future contributions of the Antarctic Ice Sheet to sea level rise will depend on the evolution of its surface mass balance (SMB), which could amplify/dampen mass losses increasingly observed at the ice sheet's edge. In situ constraints of SMB over annual-to-decadal timescales consist mostly of firn/ice cores that have a surface footprint $sim$ cm $^{2}$ .SMB constraints also come from climate models, which have a higher temporal resolution but a larger surface footprint of several km $^{2}$ .We use ice-penetrating radar data to obtain an intermediate spatial and temporal resolution SMB record over three ice rises along the Princess Ragnhild Coast. The co-located ice cores allow us to obtain absolute radar-derived SMB rates at a multi-annual-to-decadal temporal resolution. By comparing the ice core SMB measurements and the radar-derived SMB records, we determine that pointwise measurements of SMB are representative of a small surface area, $sim 200-500$ m radius extending from the ice core drill site for the ice rises studied here, and that the pointwise measurements are systematically 7–15 cm w.e. a $^{-1}$ lower than the mean SMB value calculated for the whole ice rises. However, ice core records are representative of an entire ice rise's temporal variability at the temporal resolution examined., SCOPUS: ar.j, info:eu-repo/semantics/published

The future contributions of the Antarctic Ice Sheet to sea level rise will depend on the evolution of its surface mass balance (SMB), which could amplify/dampen mass losses increasingly observed at the ice sheet's edge. In situ constraints of SMB over annual-to-decadal timescales consist mostly of firn/ice cores that have a surface footprint cm. SMB constraints also come from climate models, which have a higher temporal resolution but a larger surface footprint of several km. We use ice-penetrating radar data to obtain an intermediate spatial and temporal resolution SMB record over three ice rises along the Princess Ragnhild Coast. The co-located ice cores allow us to obtain absolute radar-derived SMB rates at a multi-annual-to-decadal temporal resolution. By comparing the ice core SMB measurements and the radar-derived SMB records, we determine that pointwise measurements of SMB are representative of a small surface area, m radius extending from the ice core drill site for the ice rises studied here, and that the pointwise measurements are systematically 7–15 cm w.e. a lower than the mean SMB value calculated for the whole ice rises. However, ice core records are representative of an entire ice rise's temporal variability at the temporal resolution examined., Mathematical Geodesy and Positioning

We present an ice-penetrating radar data set which consists of 26 internal reflecting horizons (IRHs) that cover the entire Dome C area of the East Antarctic plateau, the most extensive to date in the region. This data set uses radar surveys collected over the space of 10 years, starting with an airborne international collaboration in 2008 to explore the region, up to the detailed ground-based surveys in support of the Beyond EPICA – Oldest Ice (BE-OI) European Consortium. Through direct correlation with the EPICA-DC ice core, we date 19 IRHs that span the past four glacial cycles, from 10 ka, beginning of the Holocene, to over 350 ka, ranging from 10 % to 83 % of the ice thickness at the EPICA-DC ice core site. We indirectly date and provide stratigraphic information for seven older IRHs using a 1D ice flow inverse model, going back to an estimated 700 ka. Depth and age uncertainties are quantified for all IRHs and provided as part of the data set. The IRH data set presented in this study is available at the US Antarctic Program Data Center (USAP-DC) (https://doi.org/10.15784/601411, Cavitte et al., 2020) and represents a contribution to the SCAR AntArchitecture action group (AntArchitecture, 2017).

Based on sparse data, Titan Dome has been identified as having a higher probability of containing ice that would capture the middle Pleistocene transition (1.25 to 0.7 Ma). New aerogeophysical observations (radar and laser altimetry) collected over Titan Dome, located about 200 km from the South Pole within the East Antarctic Ice Sheet, were used to characterize the region (e.g., geometry, internal structure, bed reflectivity, and flow history) and assess its suitability as a paleoclimate ice core site. The radar coupled with an available ice core chronology enabled the tracing of dated internal reflecting horizons throughout the region, which also served as constraints on basal ice age modeling. The results of the survey revealed new basal topographic detail and better constrain the ice topographical location of Titan Dome, which differs between community datasets. Titan Dome is not expected to be relevant to the study of the middle Pleistocene transition due to a combination of past fast flow dynamics, the basal ice likely being too young, and the temporal resolution likely being too coarse if 1 Ma ice were to exist.

We present an ice-penetrating radar data set which consists of 26 internal reflecting horizons (IRHs) that cover the entire Dome C area of the East Antarctic plateau, the most extensive to date in the region. This data set uses radar surveys collected over the space of 10 years, starting with an airborne international collaboration in 2008 to explore the region, up to the detailed ground-based surveys in support of the Beyond EPICA – Oldest Ice (BE-OI) European Consortium. Through direct correlation with the EPICA-DC ice core, we date 19 IRHs that span the past four glacial cycles, from 10 ka, beginning of the Holocene, to over 350 ka, ranging from 10 % to 83 % of the ice thickness at the EPICA-DC ice core site. We indirectly date and provide stratigraphic information for seven older IRHs using a 1D ice flow inverse model, going back to an estimated 700 ka. Depth and age uncertainties are quantified for all IRHs and provided as part of the data set. The IRH data set presented in this study is available at the US Antarctic Program Data Center (USAP-DC) (https://doi.org/10.15784/601411, Cavitte et al., 2020) and represents a contribution to the SCAR AntArchitecture action group (AntArchitecture, 2017).

Large open-ocean polynyas, defined as ice-free areas within the sea ice pack, have only rarely been observed in the Southern Ocean over the past decades. In addition to smaller recent events, an impressive sequence occurred in the Weddell Sea in 1974, 1975 and 1976 with openings of more than 300 000 km2 that lasted the full winter. These big events have a huge impact on the sea ice cover, deep-water formation, and, more generally, on the Southern Ocean and the Antarctic climate. However, we have no estimate of the frequency of the occurrence of such large open-ocean polynyas before the 1970s. Our goal here is to test if polynya activity could be reconstructed using continental records and, specifically, observations derived from ice cores. The fingerprint of big open-ocean polynyas is first described in reconstructions based on data from weather stations, in ice cores for the 1970s and in climate models. It shows a signal characterized by a surface air warming and increased precipitation in coastal regions adjacent to the eastern part of the Weddell Sea, where several high-resolution ice cores have been collected. The signal of the isotopic composition of precipitation is more ambiguous; thus, we base our reconstructions on surface mass balance records alone. A first reconstruction is obtained by performing a simple average of standardized records. Given the similarity between the observed signal and the one simulated in models, we also use data assimilation to reconstruct past polynya activity. The impact of open-ocean polynyas on the continent is not large enough, compared with the changes due to factors such as atmospheric variability, to detect the polynya signal without ambiguity, and additional observations would be required to clearly discriminate the years with and without open-ocean polynya. Thus, it is reasonable to consider that, in these preliminary reconstructions, some high snow accumulation events may be wrongly interpreted as the consequence of polyn

Ice cores are an important record of the past surface mass balance (SMB) of ice sheets, with SMB mitigating the ice sheets' sea level impact over the recent decades. For the Antarctic Ice Sheet (AIS), SMB is dominated by large-scale atmospheric circulation, which collects warm moist air from further north and releases it in the form of snow as widespread accumulation or focused atmospheric rivers on the continent. This suggests that the snow deposited at the surface of the AIS should record strongly coupled SMB and surface air temperature (SAT) variations. Ice cores use δ18O as a proxy for SAT as they do not record SAT directly. Here, using isotope-enabled global climate models and the RACMO2.3 regional climate model, we calculate positive SMB–SAT and SMB–δ18O annual correlations over ∼90 % of the AIS. The high spatial resolution of the RACMO2.3 model allows us to highlight a number of areas where SMB and SAT are not correlated, and we show that wind-driven processes acting locally, such as foehn and katabatic effects, can overwhelm the large-scale atmospheric contribution in SMB and SAT responsible for the positive SMB–SAT annual correlations. We focus in particular on Dronning Maud Land, East Antarctica, where the ice promontories clearly show these wind-induced effects. However, using the PAGES2k ice core compilations of SMB and δ18O of Thomas et al. (2017) and Stenni et al. (2017), we obtain a weak annual correlation, on the order of 0.1, between SMB and δ18O over the past ∼150 years. We obtain an equivalently weak annual correlation between ice core SMB and the SAT reconstruction of Nicolas and Bromwich (2014) over the past ∼50 years, although the ice core sites are not spatially co-located with the areas displaying a low SMB–SAT annual correlation in the models. To resolve the discrepancy between the measured and modeled signals, we show that averaging the ice core records in close spatial proximity increases their SMB–SAT annual correlation. This increase sho

Over the last century, the increase in snow accumulation has partly mitigated the total dynamic Antarctic Ice Sheet mass loss. However, the mechanisms behind this increase are poorly understood. Here we analyze the Antarctic Ice Sheet atmospheric moisture budget based on climate reanalysis and model simulations to reveal that the interannual variability of regional snow accumulation is controlled by both the large-scale atmospheric circulation and short-lived synoptic-scale events (i.e. storm systems). Yet, when considering the entire continent at the multi-decadal scale, only the synoptic-scale events can explain the recent and expected future snow accumulation increase. In a warmer climate induced by climate change, these synoptic-scale events transport air that can contain more humidity due to the increasing temperatures leading to more precipitation on the continent. Our findings highlight that the multi-decadal and interannual snow accumulation variability is governed by different processes, and that we thus cannot rely directly on the mechanisms driving interannual variations to predict long-term changes in snow accumulation in the future.

To resolve the mechanisms behind the major climate reorganisation which occurred between 0.9 and 1.2 Ma, the recovery of a suitable 1.5 million-year-old ice core is fundamental. The quest for such an Oldest Ice core requires a number of key boundary conditions, of which the poorly known basal geothermal heat flux (GHF) is lacking. We use a transient thermodynamical 1D vertical model that solves for the rate of change of temperature in the vertical, with surface temperature and modelled GHF as boundary conditions. For each point on the ice sheet, the model is forced with variations in atmospheric conditions over the last 2 Ma, and modelled ice-thickness variations. The process is repeated for a range of GHF values to determine the value of GHF that marks the limit between frozen and melting conditions over the whole ice sheet, taking into account 2 Ma of climate history. These threshold values of GHF are statistically compared to existing GHF data sets. The new probabilistic GHF fields obtained for the ice sheet thus provide the missing boundary conditions in the search for Oldest Ice. High spatial resolution radar data are examined locally in the Dome Fuji and Dome C regions, as these represent the ice core community’s primary drilling sites. GHF, bedrock variability, ice thickness and other essential criteria combined highlight a dozen major potential Oldest Ice sites in the vicinity of Dome Fuji and Dome C, where GHF allows for Oldest Ice.

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-averagi

The search for ice as old as 1.5 Myr requires the identification of places that maximize our chances to retrieve old, well-resolved, undisturbed and datable ice. One of these locations is very likely southwest of the Dome C summit, where elevated bedrock makes the ice thin enough to limit basal melting. A 3-D ice flow simulation is used to calculate five selection criteria, which together delineate the areas with the most appropriate glaciological properties. These selected areas (a few square kilometers) lie on the flanks of a bedrock high, where a balance is found between risks of basal melting, stratigraphic disturbances and sufficient age resolution. Within these areas, several sites of potential 1.5 Myr old ice are proposed, situated on local bedrock summits or ridges. The trajectories of the ice particles towards these locations are short, and the ice flows over a smoothly undulating bedrock. These sites will help to choose where new high-resolution ground radar surveys should be conducted in upcoming field seasons.

To resolve the mechanisms behind the major climate reorganisation, which occurred between 0.9 and 1.2 Ma, the recovery of a suitable 1.5 million-year-old ice core is fundamental. The quest for an Oldest Ice core requires a number of key boundary conditions, of which the poorly known basal geothermal heat flux (GHF) is lacking. We use a transient thermodynamical 1-D vertical model that solves for the rate of change of temperature in the vertical, with surface temperature and modelled GHF as boundary conditions. For each point on the ice sheet, the model is forced with variations in atmospheric conditions over the last 2 Ma and modelled ice-thickness variations. The process is repeated for a range of GHF values to determine the value of GHF that marks the limit between frozen and melting conditions over the whole ice sheet, taking into account 2 Ma of climate history. These threshold values of GHF are statistically compared to existing GHF data sets. The new probabilistic GHF fields obtained for the ice sheet thus provide the missing boundary conditions in the search for Oldest Ice. High spatial resolution radar data are examined locally in the Dome Fuji and Dome C regions, as these represent the ice core community's primary drilling sites. GHF, bedrock variability, ice thickness and other essential criteria combined highlight a dozen major potential Oldest Ice sites in the vicinity of Dome Fuji and Dome C, where GHF could allow for Oldest Ice.

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 co

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., SCOPUS: ar.j, info:eu-repo/semantics/published