Antarctic Ice Cores Research Articles

Sulfur-isotope anomalies recorded in Antarctic ice cores as a potential proxy for tracing past ozone layer depletion events.

Changes in the cosmic-ray background of the Earth can impact the ozone layer. High-energy cosmic events [e.g. supernova (SN)] or rapid changes in the Earth's magnetic field [e.g. geomagnetic Excursion (GE)] can lead to a cascade of cosmic rays. Ensuing chemical reactions can then cause thinning/destruction of the ozone layer-leading to enhanced penetration of harmful ultraviolet (UV) radiation toward the Earth's surface. However, observational evidence for such UV "windows" is still lacking. Here, we conduct a pilot study and investigate this notion during two well-known events: the multiple SN event (≈10 kBP) and the Laschamp GE event (≈41 kBP). We hypothesize that ice-core-Δ33S records-originally used as volcanic fingerprints-can reveal UV-induced background-tropospheric-photochemical imprints during such events. Indeed, we find nonvolcanic S-isotopic anomalies (Δ33S ≠ 0‰) in background Antarctic ice-core sulfate during GE/SN periods, thereby confirming our hypothesis. This suggests that ice-core-Δ33S records can serve as a proxy for past ozone-layer-depletion events.

Compilation of Southern Ocean sea-ice records covering the last glacial-interglacial cycle (12–130 ka)

Abstract. Antarctic sea ice forms a critical part of the Southern Ocean and global climate system. The behaviour of Antarctic sea ice throughout the last glacial-interglacial (G-IG) cycle (12 000–130 000 years) allows us to investigate the interactions between sea ice and climate under a large range of mean climate states. Understanding both temporal and spatial variations in Antarctic sea ice across a G-IG cycle is crucial to a better understanding of the G-IG regulation of atmospheric CO2, ocean circulation, nutrient cycling and productivity. This study presents 28 published qualitative and quantitative estimates of G-IG sea ice from 24 marine sediment cores and an Antarctic ice core. Sea ice is reconstructed from the sediment core records using diatom assemblages and from the ice core record using sea-salt sodium flux. Whilst all regions of the Southern Ocean display the same overall pattern in G-IG sea-ice variations, the magnitudes and timings vary between regions. Sea-ice cover is most sensitive to changing climate in the regions of high sea-ice outflow from the Weddell Sea and Ross Sea gyres, as indicated by the greatest magnitude changes in sea ice in these areas. In contrast the Scotia Sea sea-ice cover is much more resilient to moderate climatic warming, likely due to the meltwater stratification from high iceberg flux through “iceberg alley” helping to sustain high sea-ice cover outside of full glacial intervals. The differing sensitivities of sea ice to climatic shifts between different regions of the Southern Ocean has important implications for the spatial pattern of nutrient supply and primary productivity, which subsequently impact carbon uptake and atmospheric CO2 concentrations changes across a G-IG cycle.

Rapid Warming Over East Antarctica Since the 1940s Caused by Increasing Influence of El Niño Southern Oscillation and Southern Annular Mode

El Niño Southern Oscillation (ENSO), Interdecadal Pacific Oscillations (IPO), and their phase relation with the Southern Annular Mode (SAM) largely control Antarctic climate variability. The relative roles of these climate modes remain elusive, particularly in the backdrop of global warming. In this study, we present a seasonally resolved new ice core (IND33) record of oxygen isotope (δ18O) for the past two centuries (1809–2013 CE) from coastal Dronning Maud Land (DML) to investigate the role of these climate modes in the Antarctic temperature variability and trend. Our investigation based on this record combined with available records from the DML region reveals that ∼32% variability in δ18O records is related to late spring to summer (Nov–Dec–Jan) temperature rather than the mean annual temperature. This indicates that reconstructed annual temperature based on Antarctic ice core δ18O records could be biased toward the temperature of the months/seasons of higher precipitation with low-moderate wind speed, which are suitable for better preservation of the ice core signal. We have reconstructed the DML temperature record of the past two centuries (1809–2019 CE) at an annual resolution based on the δ18O ice core record (1809–1993 CE) combined with the recent ERA5 surface air temperature record (1994–2019 CE). The reconstructed temperature anomaly record reveals a significant cooling trend in the 19th century during 1809–1907 CE with a rate of −0.164 ± 0.045°C decade−1 followed by a warming trend from the mid-20th to early 21st centuries (1942–2019 CE) with a rate of +0.452 ± 0.056°C decade−1. This long-term warming trend since the 1940s coincides with the increase in ENSO events and its strong antiphase relation with SAM, suggesting an increasing influence of SAM–ENSO coupling in modulating the DML temperature in recent decades.

Antarctic sea ice over the past 130 000 years – Part 1: a review of what proxy records tell us

Abstract. Antarctic sea ice plays a critical role in the Earth system, influencing energy, heat and freshwater fluxes, air–sea gas exchange, ice shelf dynamics, ocean circulation, nutrient cycling, marine productivity and global carbon cycling. However, accurate simulation of recent sea-ice changes remains challenging and, therefore, projecting future sea-ice changes and their influence on the global climate system is uncertain. Reconstructing past changes in sea-ice cover can provide additional insights into climate feedbacks within the Earth system at different timescales. This paper is the first of two review papers from the Cycles of Sea Ice Dynamics in the Earth system (C-SIDE) working group. In this first paper, we review marine- and ice core-based sea-ice proxies and reconstructions of sea-ice changes throughout the last glacial–interglacial cycle. Antarctic sea-ice reconstructions rely mainly on diatom fossil assemblages and highly branched isoprenoid (HBI) alkenes in marine sediments, supported by chemical proxies in Antarctic ice cores. Most reconstructions for the Last Glacial Maximum (LGM) suggest that winter sea ice expanded all around Antarctica and covered almost twice its modern surface extent. In contrast, LGM summer sea ice expanded mainly in the regions off the Weddell and Ross seas. The difference between winter and summer sea ice during the LGM led to a larger seasonal cycle than today. More recent efforts have focused on reconstructing Antarctic sea ice during warm periods, such as the Holocene and the Last Interglacial (LIG), which may serve as an analogue for the future. Notwithstanding regional heterogeneities, existing reconstructions suggest that sea-ice cover increased from the warm mid-Holocene to the colder Late Holocene with pervasive decadal- to millennial-scale variability throughout the Holocene. Studies, supported by proxy modelling experiments, suggest that sea-ice cover was halved during the warmer LIG when global average temperatures were ∼2 ∘C above the pre-industrial (PI). There are limited marine (14) and ice core (4) sea-ice proxy records covering the complete 130 000 year (130 ka) last glacial cycle. The glacial–interglacial pattern of sea-ice advance and retreat appears relatively similar in each basin of the Southern Ocean. Rapid retreat of sea ice occurred during Terminations II and I while the expansion of sea ice during the last glaciation appears more gradual especially in ice core data sets. Marine records suggest that the first prominent expansion occurred during Marine Isotope Stage (MIS) 4 and that sea ice reached maximum extent during MIS 2. We, however, note that additional sea-ice records and transient model simulations are required to better identify the underlying drivers and feedbacks of Antarctic sea-ice changes over the last 130 ka. This understanding is critical to improve future predictions.

Drivers of Last Millennium Antarctic Climate Evolution in an Ensemble of Community Earth System Model Simulations

Improved understanding of the drivers of climate variability, particularly over the last millennium, and its influence on Antarctic ice melt have important implications for projecting ice sheet resilience in a changing climate. Here, we investigated the variability in Antarctic climate and sea ice extent during the last millennium (850–1850 CE) by comparing paleoenvironmental reconstructions with simulations from the Community Earth System Model Last Millennium Ensemble (CESM-LME). Atmospheric and oceanic response to external forcing in CESM-LME simulations typically take the form of an Antarctic dipole: cooling over most of Antarctica and warming east of the Antarctic Peninsula. This configuration is also observed in ice core records. Unforced variability and a dipole response to large volcanic eruptions contribute to weaker cooling in the Antarctic than the Arctic, consistent with the absence of a strong volcanic signal in Antarctic ice core records. The ensemble does not support a clear link between the dipole pattern and baseline shifts in the Southern Annular Mode and El Niño-Southern Oscillation proposed by some paleoclimate reconstructions. Our analysis provides a point of comparison for paleoclimate reconstructions and highlights the role of internal climate variability in driving modeled last millennium climate evolution in the Antarctic.

Sunlight-driven nitrate loss records Antarctic surface mass balance

Standard proxies for reconstructing surface mass balance (SMB) in Antarctic ice cores are often inaccurate or coarsely resolved when applied to more complicated environments away from dome summits. Here, we propose an alternative SMB proxy based on photolytic fractionation of nitrogen isotopes in nitrate observed at 114 sites throughout East Antarctica. Applying this proxy approach to nitrate in a shallow core drilled at a moderate SMB site (Aurora Basin North), we reconstruct 700 years of SMB changes that agree well with changes estimated from ice core density and upstream surface topography. For the under-sampled transition zones between dome summits and the coast, we show that this proxy can provide past and present SMB values that reflect the immediate local environment and are derived independently from existing techniques.

Gas isotope thermometry in the South Pole and Dome Fuji ice cores provides evidence for seasonal rectification of ice core gas records

Abstract. Gas isotope thermometry using the isotopes of molecular nitrogen and argon has been used extensively to reconstruct past surface temperature change from Greenland ice cores. The gas isotope ratios δ15N and δ40Ar in the ice core are each set by the amount of gravitational and thermal fractionation in the firn. The gravitational component of fractionation is proportional to the firn thickness, and the thermal component is proportional to the temperature difference between the top and bottom of the firn column, which can be related to surface temperature change. Compared to Greenland, Antarctic climate change is typically more gradual and smaller in magnitude, which results in smaller thermal fractionation signals that are harder to detect. This has hampered application of gas isotope thermometry to Antarctic ice cores. Here, we present an analytical method for measuring δ15N and δ40Ar with a precision of 0.002 ‰ per atomic mass unit, a two-fold improvement on previous work. This allows us to reconstruct changes in firn thickness and temperature difference at the South Pole between 30 and 5 kyr BP. We find that variability in firn thickness is controlled in part by changes in snow accumulation rate, which is, in turn, influenced strongly by the along-flowline topography upstream of the ice core site. Variability in our firn temperature difference record cannot be explained by annual-mean processes. We therefore propose that the ice core gas isotopes contain a seasonal bias due to rectification of seasonal signals in the upper firn. The strength of the rectification also appears to be linked to fluctuations in the upstream topography. As further evidence for the existence of rectification, we present new data from the Dome Fuji ice core that are also consistent with a seasonal bias throughout the Holocene. Our findings have important implications for the interpretation of ice core gas records. For example, we show that the effects of upstream topography on ice core records can be significant at flank sites like the South Pole – they are responsible for some of the largest signals in our record. Presumably upstream signals impact other flank-flow ice cores such as EDML, Vostok, and EGRIP similarly. Additionally, future work is required to confirm the existence of seasonal rectification in polar firn, to determine how spatially and temporally widespread rectifier effects are, and to incorporate the relevant physics into firn air models.

Regional validation of the use of diatoms in ice cores from the Antarctic Peninsula as a Southern Hemisphere westerly wind proxy

Abstract. The Southern Hemisphere westerly winds are among the most important drivers of recently observed environmental changes in West Antarctica. However, the lack of long-term wind records in this region hinders our ability to assess the long-term context of these variations. Ice core proxy records yield valuable information about past environmental changes, although current proxies present limitations when aiming to reconstruct past winds. Here we present the first regional wind study based on the novel use of diatoms preserved in Antarctic ice cores. We assess the temporal variability in diatom abundance and its relation to regional environmental parameters spanning a 20-year period across three sites in the southern Antarctic Peninsula and Ellsworth Land, Antarctica. Correlation analyses reveal that the temporal variability of diatom abundance from high-elevation ice core sites is driven by changes in wind strength over the core of the Southern Hemisphere westerly wind belt, validating the use of diatoms preserved in ice cores from high-elevation inland sites in the southern Antarctic Peninsula and Ellsworth Land as a proxy for reconstructing past variations in wind strength over the Pacific sector of the Southern Hemisphere westerly wind belt.

The Holocene hypsithermal in the Australian region

Close examination of key and well-dated Holocene sites, both on land and at sea in the Australian region indicates that at the very beginning of the Holocene, as a result of strong westerlies, there must have been a continuous positive Southern Annular Mode (SAM). Following from that, the entire region switched to a negative SAM scenario and, during that time, the westerlies must have retreated further south. Afterwards, a period of time spaning ∼8200 to ∼5500 years ago temperatures were higher than today. We refer to it as the Holocene Hyspithermal. Coincident to this period, lake levels and postulated rainfall were extraordinarily high and vegetation spectra in places very different compared to today. The extent of this period varies by a few centuries between sites, but this may result from the level of resolution and also appears to be controlled by latitude. There is also clear indication that the influence of the westerlies was reduced over Australia during those two and a half millennia.Nevertheless, air temperatures recognised in Antarctic ice cores are at the opposite to those recognised in Australia. In addition, during the Australian Holocene Hypsithermal, CO2 levels were at their lowest in Antarctic ice cores.Climatic conditions then progressively deteriorated everywhere a bit after ∼6000 years BP until recent times as ENSO signals with alternating El Niño and La Niña conditions across the entire Pacific region as already described by Perner et al. (2018) based on the same cores studied here.Brief mention is also made to the presence of humans in SE Australia during the Holocene. It seems that human activities changed well after the period of high temperatures and rainfall, with more sedentary activities along the major rivers, with an enhancement of food production in organized settings suggestive of villages.

Can Ammonium Records in Antarctic Ice Cores be a Proxy for Sea Ice Fluctuations?

Abstract Ammonium (NH4+) accumulation and variability in Antarctic ice cores possess signals of various biological, anthropogenic, and atmospheric processes. However, the use of NH4+ variability in Antarctic ice cores for paleoclimate reconstructions is not well established. Addressing this, we have utilized, published NH4+ ice core records from the Dronning Maud Land (DML), Antarctica to assess the sources of NH4+, its transport, major forcing factors, and their mechanism in controlling the variability of ammonium signals in coastal ice cores. The first principle component (PC1) derived from PC analysis of NH4+ data comprises ~43% of total variability and also depicts a significant spatial correlation with the NSIDC Sea ice concentration (SIC) variability in the Weddell Sea region of the Southern Ocean. In addition, the PC1 shows strong 2-9 year periodicities related to El Niño Southern Oscillation along with two major decadal to centennial periodicities of ~40 and 100 years. The control of sea ice fluctuations and related oceanic processes (ENSO) provide evidence of the marine sources and influence of ocean-related processes as a major controlling factor for the NH4+ variability in DML ice cores. The cumulative accumulation of the NH4+ produced from enhanced biological decomposition, and unutilized ammonium during nitrification under high SIC conditions are possible mechanisms for the increased supply. The correlation of ERA5 surface air temperatures with Nino3.4 suggests that during El Niño conditions, the air temperatures significantly dropped in the Weddell Sea sector, fostering the high SIC conditions. The findings of this study regarding the role of the SIC in the ammonium variability in Antarctic ice cores would help us better interpret the paleo NH4+ records and their applications to deduce the paleoclimate conditions in coastal Antarctic regions.

Dust correlation and oxygen isotope stratigraphy in the Southern Ocean over the last 450 kyrs: An Indian sector perspective

The chronology of Southern Ocean (SO) marine sediment cores forms the basis of understanding the SO paleoceanography, with significant implications for global climate. Because tuning of the oxygen isotope (δ18O) record of a marine sediment core to a δ18O stack (δ18O stratigraphy) is difficult in the SO because of a general paucity of calcareous foraminifera, tuning of the dust proxy signal of a marine sediment core to the dust record of an ice core (dust correlation) is a promising way to construct an age–depth model. However, the reliability of dust correlation has not been established, especially beyond the last ∼300 kyrs, and such work has been performed more in the Atlantic and Pacific sectors of the SO than in the Indian sector. Here we present a new dust correlation using the rock magnetic record, together with continuous δ18O stratigraphy, for a marine sediment core in the Indian sector of the SO over the last ∼410 kyrs. The δ18O stratigraphy is consistent with the dust correlation within their chronological uncertainties, supporting the reliability of the latter chronology for glacial–interglacial timescales. However, the dust correlation often produces older ages than the δ18O stratigraphy by up to 3 kyrs (410–126 ka and 84–50 ka). We additionally compiled, based on dust correlation, available dust proxy records in the circum-Antarctic SO for the last 450 kyrs to discuss variations in the dust proxy record within the Indian sector and between all three sectors. For the Indian sector, there is a marked difference in dust proxy signals for Marine Isotope Stage (MIS) 6 between two cores that have markedly different latitudinal positions relative to the oceanic fronts. Increased biogenic magnetite production by iron fertilization in the Subantarctic Zone during MIS 6 may partly explain the difference between the two cores. For the circum-Antarctic SO, moving correlation coefficients were computed between the marine sediment-core dust proxy signals and the dust flux of an Antarctic ice core. Strong correlation was recognized for intervals of high dust flux in the ice core, except for MIS 6. Although minor, the slightly lower correlation in the Indian sector than in the other sectors indicates a contribution from local dust source and volcanic materials and thus suggests the necessity for caution when performing dust correlation in the Indian sector of the SO.

Polarimetric radar reveals the spatial distribution of ice fabric at domes and divides in East Antarctica

Abstract. Ice crystals are mechanically and dielectrically anisotropic. They progressively align under cumulative deformation, forming an ice-crystal-orientation fabric that, in turn, impacts ice deformation. However, almost all the observations of ice fabric are from ice core analysis, and its influence on the ice flow is unclear. Here, we present a non-linear inverse approach to process co- and cross-polarized phase-sensitive radar data. We estimate the continuous depth profile of georeferenced ice fabric orientation along with the reflection ratio and horizontal anisotropy of the ice column. Our method approximates the complete second-order orientation tensor and all the ice fabric eigenvalues. As a result, we infer the vertical ice fabric anisotropy, which is an essential factor to better understand ice deformation using anisotropic ice flow models. The approach is validated at two Antarctic ice core sites (EPICA (European Project for Ice Coring in Antarctica) Dome C and EPICA Dronning Maud Land) in contrasting flow regimes. Spatial variability in ice fabric characteristics in the dome-to-flank transition near Dome C is quantified with 20 more sites located along with a 36 km long cross-section. Local horizontal anisotropy increases under the dome summit and decreases away from the dome summit. We suggest that this is a consequence of the non-linear rheology of ice, also known as the Raymond effect. On larger spatial scales, horizontal anisotropy increases with increasing distance from the dome. At most of the sites, the main driver of ice fabric evolution is vertical compression, yet our data show that the horizontal distribution of the ice fabric is consistent with the present horizontal flow. This method uses polarimetric-radar data, which are suitable for profiling radar applications and are able to constrain ice fabric distribution on a spatial scale comparable to ice flow observations and models.

Variation in magnetic susceptibility in the Bellingshausen Sea continental rise since the last glacial period and implications for terrigenous material input mechanisms

Magnetic susceptibility (MS) is routinely used as a proxy for terrigenous material input to Antarctic continental margin sediments. Variations in terrigenous input on glacial-interglacial timescales in this setting are related to a range of environmental changes; including to the cryosphere, atmosphere, and hydrosphere. However, the exact environmental factors and depositional processes that control MS values in sedimentary cores from the Antarctic continental margin remains unresolved. Here, we explore in detail sedimentary physical characteristics and MS values of core BS17-GC01 collected from Bellingshausen Sea continental rise, to understand how depositional process have varied during the late Quaternary. We found that input of terrigenous material increased markedly during Marine Isotope Stages (MIS) 2 and 4, compared to deposition during the adjacent interstadial and Holocene times. Variation in the bulk MS values were associated with changes in grain size. Size specific mass-normalized MS measurements show the fine fraction (<16 μm) has the highest MS values, followed by the coarse silt fraction (16–63 μm). However, because coarse silt is the predominant grain size, the increased bulk MS during MIS 2 and 4 are mainly driven by changes in this fraction. We propose that the coarse silt with high bulk MS values observed in glacial sediments were transported to the site as iceberg rafted debris, originating from the Bellingshausen Sea shelf region. In contrast, grains >1 mm – a commonly applied proxy for glacial transport in this setting – were observed during MIS 1 and 3 only. We conclude these contrasting signatures of iceberg rafted debris in this Bellingshausen Sea core arise from changes in the to the erosive dynamics of the adjacent ice shelf and continental shelf on glacial/interglacial timescales. • The MS in the Bellingshausen Sea co-varied with Antarctic ice core dust record. • The correlations of MS and dust record can be used for obtaining chronological points. • The MS is controlled by source/transport mechanism changes resulting in grain size change. • MS increases during MIS 2 and 4 are ascribed to increased coarse silt having high MS. • The coarse silt fraction during glacial periods was transported by icebergs.

Annually Resolved Profiles of δ34S and Sulfate in Shallow Ice Core DF01 (Dome Fuji, Antarctica) Spanning the Nineteenth Century and Their Geochemical Implications

AbstractWe report the results of a high‐sensitivity sulfur‐isotope (δ34S) analysis of samples from shallow ice core DF01 from Dome Fuji, East Antarctica, and present chronological profiles of the data with an annual resolution covering the entire nineteenth century. In the δ34S and non‐sea‐salt (nss)‐sulfate ([SO42−]nss) profiles, the signals of five well‐known nineteenth century volcanic eruptions were confirmed: Krakatau 1883 (and Tarawera 1886), Makian 1861, Cosiguina 1835, Tambora 1815, and Unknown 1809. Due to large influence of sulfate caused by volcanic eruptions in the first half of the nineteenth century, the averages of nss‐sulfur‐isotope composition (δ34Snss) and [SO42−]nss for the second half (+14.6‰ and 67.6 μg/L), calculated after excluding these prominent signals, were found to differ from those of the first half century (+9.1‰ and 84.1 μg/L). We also examined anomalous variations other than the five major volcanic signals in the δ34S and [SO42−]nss profiles based on not only variations in δ34S and [SO42−]nss but also variations in the concentrations of other ions. We then identified the presence of three additional volcanic signals: Babuyan 1831, Galunggung 1822, and an unknown volcanic eruption around 1807. The latter two signals confirmed for the first time in Antarctic ice cores. Other δ34S and [SO42−]nss anomalies may be related to an influx of coastal air masses, a nonstratospheric volcanic eruption near Antarctica or an increase in the relative contribution of a stratospheric component to the precipitation in inland Antarctica. These demonstrate the importance of obtaining the annually resolved, continuous δ34S profile over a certain period.

Antiphased dust deposition and productivity in the Antarctic Zone over 1.5 million years

The Southern Ocean paleoceanography provides key insights into how iron fertilization and oceanic productivity developed through Pleistocene ice-ages and their role in influencing the carbon cycle. We report a high-resolution record of dust deposition and ocean productivity for the Antarctic Zone, close to the main dust source, Patagonia. Our deep-ocean records cover the last 1.5 Ma, thus doubling that from Antarctic ice-cores. We find a 5 to 15-fold increase in dust deposition during glacials and a 2 to 5-fold increase in biogenic silica deposition, reflecting higher ocean productivity during interglacials. This antiphasing persisted throughout the last 25 glacial cycles. Dust deposition became more pronounced across the Mid-Pleistocene Transition (MPT) in the Southern Hemisphere, with an abrupt shift suggesting more severe glaciations since ~0.9 Ma. Productivity was intermediate pre-MPT, lowest during the MPT and highest since 0.4 Ma. Generally, glacials experienced extended sea-ice cover, reduced bottom-water export and Weddell Gyre dynamics, which helped lower atmospheric CO2 levels.

Magnitude, frequency and climate forcing of global volcanism during the last glacial period as seen in Greenland and Antarctic ice cores (60–9 ka)

Abstract. Large volcanic eruptions occurring in the last glacial period can be detected by their accompanying sulfuric acid deposition in continuous ice cores. Here we employ continuous sulfate and sulfur records from three Greenland and three Antarctic ice cores to estimate the emission strength, the frequency and the climatic forcing of large volcanic eruptions that occurred during the second half of the last glacial period and the early Holocene, 60–9 kyr before 2000 CE (b2k). Over most of the investigated interval the ice cores are synchronized, making it possible to distinguish large eruptions with a global sulfate distribution from eruptions detectable in one hemisphere only. Due to limited data resolution and large variability in the sulfate background signal, particularly in the Greenland glacial climate, we only list Greenland sulfate depositions larger than 20 kg km−2 and Antarctic sulfate depositions larger than 10 kg km−2. With those restrictions, we identify 1113 volcanic eruptions in Greenland and 737 eruptions in Antarctica within the 51 kyr period – for which the sulfate deposition of 85 eruptions is found at both poles (bipolar eruptions). Based on the ratio of Greenland and Antarctic sulfate deposition, we estimate the latitudinal band of the bipolar eruptions and assess their approximate climatic forcing based on established methods. A total of 25 of the identified bipolar eruptions are larger than any volcanic eruption occurring in the last 2500 years, and 69 eruptions are estimated to have larger sulfur emission strengths than the Tambora, Indonesia, eruption (1815 CE). Throughout the investigated period, the frequency of volcanic eruptions is rather constant and comparable to that of recent times. During the deglacial period (16–9 ka b2k), however, there is a notable increase in the frequency of volcanic events recorded in Greenland and an obvious increase in the fraction of very large eruptions. For Antarctica, the deglacial period cannot be distinguished from other periods. This confirms the suggestion that the isostatic unloading of the Northern Hemisphere (NH) ice sheets may be related to the enhanced NH volcanic activity. Our ice-core-based volcanic sulfate records provide the atmospheric sulfate burden and estimates of climate forcing for further research on climate impact and understanding the mechanism of the Earth system.

Regional variability of diatoms in ice cores from the Antarctic Peninsula and Ellsworth Land, Antarctica

Abstract. The presence of marine microfossils (diatoms) in glacier ice and ice cores has been documented from numerous sites in Antarctica, Greenland, as well as from sites in the Andes and the Altai mountains, and attributed to entrainment and transport by winds. However, their presence and diversity in snow and ice, especially in polar regions, are not well documented and still poorly understood. Here we present the first data to resolve the regional and temporal distribution of diatoms in ice cores, spanning a 20-year period across four sites in the Antarctic Peninsula and Ellsworth Land, Antarctica. We assess the regional variability in diatom composition and abundance at annual and sub-annual resolution across all four sites. These data corroborate the prevalence of contemporary marine diatoms in Antarctic Peninsula ice cores, reveal that the timing and amount of diatoms deposited vary between low- and high-elevation sites, and support existing evidence that marine diatoms have the potential to yield a novel palaeoenvironmental proxy for ice cores in Antarctica.

The dynamics of warming during the last deglaciation in high-elevation regions of Eastern Equatorial Africa

Tropical mountain environments, such as the Rwenzori Mountains in equatorial Africa, are thought to be particularly sensitive to climate change. Ongoing warming in the Rwenzori is impacting local environments and communities through glacial retreat, fires, and flooding. Paleoclimate reconstructions from elsewhere in Africa suggest considerable warming accompanied glacier retreat during the last glacial termination, from ∼21 thousand years before present (ka) through the early to mid-Holocene. Quantifying these changes has been difficult but could help to assess future impacts in the Rwenzori. Here, we present a ∼21 thousand-year (kyr) temperature reconstruction based on the relative abundance of branched glycerol dialkyl glycerol tetraethers (brGDGTs) from Lake Mahoma (2,990 m above sea level; m asl) in the Rwenzori Mountains, Uganda. Our record, paired with existing Rwenzori glacial moraine 10Be exposure ages, suggests that deglacial warming and glacial retreat began by ∼20 ka and accelerated at ∼18–18.5 ka. The timing of the onset of rapid warming matches the timing of the post-glacial rise in radiative forcing from atmospheric greenhouse gases (GHGs) from Antarctic ice cores (Brook et al., 1996; Marcott et al., 2014; Monnin et al., 2004; Schilt et al., 2010). Our temperature reconstruction registers ∼4.9 °C warming from the Last Glacial Maximum (LGM) to the late Holocene. This increase is larger than the average ∼2-4 °C warming observed in records from lower elevation sites in tropical East Africa, but similar to that observed at other high-elevation sites in this region. The increased warming at higher elevations thus confirms that the temperature lapse rate steepened during the LGM over this region. Our results also indicate ∼3 °C of warming during the mid-Holocene relative to the late Holocene. This suggests that the freezing-level height rose above Rwenzori summit elevations at that time, likely causing complete deglaciation of the Rwenzori Mountains from ∼5 to 7 ka. The mid-Holocene is thus a potential analog for the glacial and environmental changes that these mountains are likely to experience in the coming decades. Overall, the timing and magnitude of temperature change observed in our record has important implications for climate model projections of future warming in tropical Africa.

Pacific decadal variability over the last 2000 years and implications for climatic risk

The Interdecadal Pacific Oscillation, an index which defines decadal climate variability throughout the Pacific, is generally assumed to have positive and negative phases that each last 20-30 years. Here we present a 2000-year reconstruction of the Interdecadal Pacific Oscillation, obtained using information preserved in Antarctic ice cores, that shows negative phases are short (7 ± 5 years) and infrequent (occurring 10% of the time) departures from a predominantly neutral-positive state that lasts decades (61 ± 56 years). These findings suggest that Pacific Basin climate risk is poorly characterised due to over-representation of negative phases in post-1900 observations. We demonstrate the implications of this for eastern Australia, where drought risk is elevated during neutral-positive phases, and highlight the need for a re-evaluation of climate risk for all locations affected by the Interdecadal Pacific Oscillation. The initiation and future frequency of negative phases should also be a research priority given their prevalence in more recent centuries.

New constraints on the postglacial shallow-water carbonate accumulation in the Great Barrier Reef

More accurate global volumetric estimations of shallow-water reef deposits are needed to better inform climate and carbon cycle models. Using recently acquired datasets and International Ocean Discovery Program (IODP) Expedition 325 cores, we calculated shallow-water CaCO3 volumetrics and mass for the Great Barrier Reef region and extrapolated these results globally. In our estimates, we include deposits that have been neglected in global carbonate budgets: Holocene Halimeda bioherms located on the shelf, and postglacial pre-Holocene (now) drowned coral reefs located on the shelf edge. Our results show that in the Great Barrier Reef alone, these drowned reef deposits represent ca. 135 Gt CaCO3, comparatively representing 16–20% of the younger Holocene reef deposits. Globally, under plausible assumptions, we estimate the presence of ca. 8100 Gt CaCO3 of Holocene reef deposits, ca. 1500 Gt CaCO3 of drowned reef deposits and ca. 590 Gt CaCO3 of Halimeda shelf bioherms. Significantly, we found that in our scenarios the periods of pronounced reefal mass accumulation broadly encompass the occurrence of the Younger Dryas and periods of CO2 surge (14.9–14.4 ka, 13.0–11.5 ka) observed in Antarctic ice cores. Our estimations are consistent with reef accretion episodes inferred from previous global carbon cycle models and with the chronology from reef cores from the shelf edge of the Great Barrier Reef.