Stratigraphic templates for ice core records of the past 1.5 million years

The international ice core community has a target to obtain continuous ice cores stretching back as far as 1.5 million years. This would provide vital data (including a CO2 profile) allowing us to assess ideas about the cause of the Mid-Pleistocene Transition (MPT). The European Beyond EPICA project and the Australian Million Year Ice Core project each plan to drill such a core in the region known as Little Dome C. Dating the cores will be challenging, and one approach will be to match some of the records obtained with existing marine sediment datasets, informed by similarities in the existing 800 kyr period. Water isotopes in Antarctica have been shown to closely mirror deepwater temperature, estimated from Mg / Ca ratios of benthic foraminifera, in a marine core on the Chatham Rise near to New Zealand. The dust record in ice cores resembles very closely a South Atlantic marine record of iron accumulation rate. By assuming these relationships continue beyond 800 ka, our ice core record could be synchronised to dated marine sediments. This could be supplemented, and allow synchronisation at higher resolution, by the identification of rapid millennial scale-events that are observed both in Antarctic methane records and in emerging records of planktic oxygen isotopes and alkenone sea surface temperature (SST) from the Portuguese Margin. Although published data remain quite sparse, it should also be possible to match 10Be from ice cores to records of geomagnetic palaeointensity and authigenic 10Be/9Be in marine sediments. However, there are a number of issues that have to be resolved before the ice core 10Be record can be used. The approach of matching records to a template will be most successful if the new core is in stratigraphic order, but should also provide constraints on disordered records, if used in combination with absolute radiogenic ages.

New insights into the ∼ 74 ka Toba eruption from sulfur isotopes of polar ice cores

The ∼74 ka Toba eruption was one of the largest volcanic events of the Quaternary. There is much interest in determining the impact of such a large event, particularly on the climate and hominid populations at the time. Although the Toba eruption has been identified in both land and marine archives as the Youngest Toba Tuff, its precise place in the ice core record is ambiguous. Several volcanic sulfate signals have been identified in both Antarctic and Greenland ice cores and span the Toba eruption 40Ar/39Ar age uncertainty. Here, we measure sulfur isotope compositions in Antarctic ice samples from the Dome C (EDC) and Dronning Maud Land (EDML) ice cores at high temporal resolution across 11 of these potential Toba sulfate peaks to identify candidates with sulfur mass-independent fractionation (S-MIF), indicative of an eruption whose plume reached altitudes at or above the stratospheric ozone layer. Using this method, we identify several candidate sulfate peaks that contain stratospheric sulfur. We further narrow down potential candidates based on the isotope signatures by identifying sulfate peaks that are due to a volcanic event at tropical latitudes. In one of these sulfate peaks at 73.67 ka, we find the largest ever reported magnitude of S-MIF in volcanic sulfate in polar ice, with a Δ33S value of −4.75 ‰. As there is a positive correlation between the magnitude of the S-MIF signal recorded in ice cores and eruptive plume height, this could be a likely candidate for the Toba super-eruption, with a plume top height in excess of 45 km. These results support the 73.7±0.3 ka (1σ) 40Ar/39Ar age estimate for the eruption, with ice core ages of our candidates with the largest magnitude S-MIF at 73.67 and 73.74 ka. Finally, since these candidate eruptions occurred on the transition into Greenland Stadial 20, the relative timing suggests that Toba was not the trigger for the large Northern Hemisphere cooling at this time although we cannot rule out an amplifying effect.

Thermal legacy of a large paleolake in Taylor Valley, East Antarctica, as evidenced by an airborne electromagnetic survey

Previous studies of the lakes of the McMurdo Dry Valleys have attempted to constrain lake level history, and results suggest the lakes have undergone hundreds of meters of lake level change within the last 20 000 years. Past studies have utilized the interpretation of geologic deposits, lake chemistry, and ice sheet history to deduce lake level history; however a substantial amount of disagreement remains between the findings, indicating a need for further investigation using new techniques. This study utilizes a regional airborne resistivity survey to provide novel insight into the paleohydrology of the region. Mean resistivity maps revealed an extensive brine beneath the Lake Fryxell basin, which is interpreted as a legacy groundwater signal from higher lake levels in the past. Resistivity data suggest that active permafrost formation has been ongoing since the onset of lake drainage and that as recently as 1500–4000 years BP, lake levels were over 60 m higher than present. This coincides with a warmer-than-modern paleoclimate throughout the Holocene inferred by the nearby Taylor Dome ice core record. Our results indicate Mid to Late Holocene lake level high stands, which runs counter to previous research finding a colder and drier era with little hydrologic activity throughout the last 5000 years.

Brief communication: New radar constraints support presence of ice older than 1.5 Myr at Little Dome C

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.

New insights into the ~74 ka Toba eruption from sulfur isotopes of polar ice cores

The ~74 ka Toba eruption was one of the largest volcanic events of the Quaternary. There is much interest in determining the impact of such a huge event, particularly on the climate and hominid populations at the time. Although the Toba eruption has been identified in both land and marine archives as the Youngest Toba Tuff, its precise place in the ice core record is ambiguous. Multiple volcanic sulfate signals have been identified in both Antarctic and Greenland ice cores within the uncertainty of age estimates as possible events for the Toba eruption. We measure sulfur isotope compositions in Antarctic ice samples at high temporal resolution across 11 of these potential Toba sulfate peaks in two cores to identify candidates with sulfur mass-independent fractionation (S-MIF), indicative of an eruption whose plume reached altitudes at or above the ozone layer in the stratosphere. Using this method, we identify several candidate sulfate peaks that contain stratospheric sulfur. We further narrow down potential candidates based on the isotope signatures by identifying sulfate peaks that are due to a volcanic event at tropical latitudes. In one of these sulfate peaks at 73.67 ka, we find the largest ever reported magnitude of S-MIF in volcanic sulfate in polar ice, with a Δ33S value of −4.75 ‰. As there is a positive correlation between the magnitude of the S-MIF signal recorded in ice cores and eruptive plume height, this could be a likely candidate for the Toba supereruption, with a plume height in excess of 45 km. These results support the 73.7 ± 0.3 ka (1σ) ka Ar/Ar age estimate for the eruption, with ice core ages of our candidates with the largest magnitude S-MIF at 73.67 and 73.74 ka. Finally, since these candidate eruptions occurred on the transition into Greenland Stadial 20, the relative timing suggests that Toba was not the trigger for the large Northern Hemisphere cooling at this time although we cannot rule out an amplifying effect.

Drought-induced biomass burning as a source of black carbon to the central Himalaya since 1781 CE as reconstructed from the Dasuopu ice core

Himalayan glaciers are melting due to atmospheric warming, with the potential to limit access to water for more than 25 % of the global population that resides in these glacier meltwater catchments. Black carbon has been implicated as a factor that is contributing to Himalayan glacier melt, but its sources and mechanisms of delivery to the Himalayas remain controversial. Here, we provide a 211-year ice core record spanning 1781–1992 CE for refractory black carbon (rBC) deposition from the Dasuopu glacier ice core that has to date provided the highest-elevation ice core record (7200 m). We report an average rBC concentration of 1.5 µg L−1 (SD=5.0, n=1628) over the 211-year period. An increase in the frequency and magnitude of rBC deposition occurs after 1877 CE, accompanied by decreased snow accumulation associated with a shift in the North Atlantic Oscillation Index to a positive phase. Typically, rBC is deposited onto Dasuopu glacier during the non-monsoon season, and short-lived increases in rBC concentration are associated with periods of drought within neighboring regions in northwestern India, Afghanistan, and Pakistan. Using a combination of spectral and back-trajectory analyses, as well as a comparison with a concurrent analysis of trace metals at equivalent depths in the same ice core, we show that biomass burning resulting from dry conditions is a source of rBC to the central Himalaya and is responsible for deposition that is up to 60 times higher than the average rBC concentration over the time period analyzed. We suggest that biomass burning is a significant source of rBC to the central Himalaya and that the rBC record can be used to identify periods of drought in nearby regions that are upwind of Dasuopu glacier.

On the relationship between stable isotopes and major impurity species as inferred from a two-dimensional firn sampling approach at EDML, East Antarctica

<p>Ice cores constitute a major palaeoclimate archive by recording, among many others, the atmospheric variations of stable oxygen and hydrogen isotopic composition of water and of soluble ionic impurities. While impurities are used as proxies for, e.g., variations in sea ice, marine biological activity and volcanism, stable isotope records are the main source of information for the reconstruction of polar temperature changes.</p><p>However, such reconstruction efforts are complicated by the fact that temperature is by far not the only driver of isotopic composition changes. A single isotopic ice-core record will comprise variations caused by a multitude of processes, from variable atmospheric circulation and moisture pathways to the intermittency of precipitation and finally to the mixing and re-location of surface snow by wind drift (stratigraphic noise). Under the assumption that specific trace components are originally deposited with the precipitated snow and its isotopic composition, the retrieved impurity records should display a similar spatial and seasonal to interannual variability as the isotope records, caused by local stratigraphic noise as well as the time-variable and intermittent precipitation patterns, respectively.</p><p>In this contribution, we investigate the possible relationship between isotope and impurity data at the East Antarctic low-accumulation site EDML. We sampled and analysed isotopic composition and major impurity species on a four metre deep and 50 metre long trench. This enables us (1) to study the spatial (horizontal times vertical) relationship in the data, and (2) to analyse and compare the seasonal and interannual variability after removing the strong contribution of local stratigraphic noise. By this, the study improves our understanding of the depositional mechanisms that play an important role for the formation of ice-core records, and it offers to investigate the potential of using impurities to correct isotopic variability in order to improve temperature reconstructions.</p>

Influence of Enso in Perú's Cordillera Blanca Glaciers

<p>In areas located over 2000 m.a.s.l., the warm phase of ENSO (El Niño) is characterized by a decrease in precipitation and an increase in temperature which can reach values above the annual average, while in the cold phase of ENSO (La Niña), precipitation increases and temperature decreases compared to the annual average. In both cases ENSO has an influence on the glacier evolution of the Andes.</p><p>The objective of the present investigation is to determine the influence of ENSO in the Cordillera Blanca through satellite images (glacier coverage delimitation) and climatic proxy (ice core) in the Shallap and Artesonraju glaciers respectively for the hydrological years between 2009/2010 to 2018/2019.</p><p>The climate analysis in both glaciers showed higher annual temperatures and lower precipitation, revealing the influence of the 2015/2016 El Niño on the studied glaciers. There was a prominent reduction in glacier coverage in Shallap, which is supported by the ice core record extracted from Artesonraju, presenting an equivalent accumulated water decrease and an <sup>18</sup>O enrichment for this period. These findings point out the influence of the 2015/2016 El Niño that significantly reduced the glacier coverage in both studied areas. On the other hand, the 2011/2012 La Niña event displayed the opposite effect, that is, colder temperatures, less glacier coverage reduction, an increase in the volume of accumulated water and an impoverishment of <sup>18</sup>O.</p><p>Given the results, it can be affirmed that during an El Niño year the loss of glacier coverage is greater, causing less equivalent water accumulation and an enrichment of <sup>18</sup>O; inversely for a La Niña year. These results support previous findings shown in research about glaciers in Peru.</p>

A smaller glacial ocean carbon inventory?

<p>Previous studies attempting to explain Pleistocene atmospheric CO<sub>2</sub> variations have focused on mechanisms that transfer carbon (C) between the oceanic, atmospheric and terrestrial reservoirs, with the underlying assumption that the total C inventory in these three Earth’s surface reservoirs remained constant during glacial-interglacial cycles. Under this framework, ocean C inventory would have been marginally increased by 500-1000 GtC (1-2%) during the glacial period. Here, we show that past ocean C inventory can be revealed by reconstructed bulk ocean <sup>14</sup>C/<sup>12</sup>C (denoted as ∆<sup>14</sup>C) and atmospheric <sup>14</sup>C production rates with an Earth system model - cGENIE. First, we develop a bulk ocean ∆<sup>14</sup>C record that spans the last 40 ka from thousands of benthic foraminifera and deep sea coral ∆<sup>14</sup>C data with a fairly good coverage of the global seafloor. We then run cGENIE under constant pre-industrial boundary conditions, with the only forcing being atmospheric <sup>14</sup>C production rates reconstructed by geomagnetic field intensity records and ice core record of <sup>10</sup>Be fluxes. Under most of the <sup>14</sup>C production scenarios, the simulated bulk ocean ∆<sup>14</sup>C are significantly lower than our composite during the Last Glacial Maximum as well as the early deglaciation. Bulk ocean ∆<sup>14</sup>C is a metric controlled by <sup>14</sup>C production rates and ocean C inventory, with the state of ocean circulation playing a minor role.  Our finding suggests either glacial <sup>14</sup>C production was much higher and/or glacial C inventory was much lower than previously thought. Implications of both possibilities are discussed. In particular, the second possibility highlight the exchange of C and ALK between Earth’s surface and geological reservoirs as a critical missing piece in searching for a complete theory of glacial-interglacial atmospheric CO<sub>2 </sub>variability.</p>

Brief Communication: New radar constraints support presence of ice older than 1.5 Ma at Little Dome C

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, an area ~35 km south-east 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 VHF radar operating in the frequency range of 170–230 MHz. This new instrument is able to detect reflections in the near-basal region, where previous surveys were unable to trace continuous 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. 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 ka at BELDC. Stratigraphy is well preserved to 2565 m, below which there is a basal unit with unknown properties. A simple ice flow model tuned to the isochrones suggests ages likely reach 1.5 Ma near 2525 m, ~40 m above the basal unit and ~240 m above the bed, with sufficient resolution (14±1 ka m−1) to resolve 41 ka glacial cycles.