Search

Changes in sea ice conditions and atmospheric circulation over the Southern Ocean play an important role in modulating Antarctic climate. However, observations of both sea ice and wind conditions are limited in Antarctica and the Southern Ocean, both temporally and spatially, prior to the satellite era (1970 onwards). Ice core chemistry data can be used to reconstruct changes over annual, decadal, and millennial timescales. To facilitate sea ice and wind reconstructions, the CLIVASH2k (CLimate Variability in Antarctica and the Southern Hemisphere over the past 2000 years) working group has compiled a database of two species, sodium [Na+] and sulfate [SO42-], commonly measured ionic species. The database (https://doi.org/10.5285/9E0ED16E-F2AB-4372-8DF3-FDE7E388C9A7; Thomas et al., 2022) comprises records from 105 Antarctic ice cores, containing records with a maximum age duration of 2000 years. An initial filter has been applied, based on evaluation against sea ice concentration, geopotential height (500 hPa), and surface wind fields to identify sites suitable for reconstructing past sea ice conditions, wind strength, or atmospheric circulation.

The injection of sulfur into the stratosphere by volcanic eruptions is the dominant driver of natural climate variability on interannual to multidecadal timescales. Based on a set of continuous sulfate and sulfur records from a suite of ice cores from Greenland and Antarctica, the HolVol v.1.0 database includes estimates of the magnitudes and approximate source latitudes of major volcanic stratospheric sulfur injection (VSSI) events for the Holocene (from 9500 BCE or 11 500 years BP to 1900 CE), constituting an extension of the previous record by 7000 years. The database incorporates new-generation ice-core aerosol records with a sub-annual temporal resolution and a demonstrated sub-decadal dating accuracy and precision. By tightly aligning and stacking the ice-core records on the WD2014 chronology from Antarctica, we resolve long-standing inconsistencies in the dating of ancient volcanic eruptions that arise from biased (i.e., dated too old) ice-core chronologies over the Holocene for Greenland. We reconstruct a total of 850 volcanic eruptions with injections in excess of 1 teragram of sulfur (Tg S); of these eruptions, 329 (39 %) are located in the low latitudes with bipolar sulfate deposition, 426 (50 %) are located in the Northern Hemisphere extratropics (NHET) and 88 (10 %) are located in the Southern Hemisphere extratropics (SHET). The spatial distribution of the reconstructed eruption locations is in agreement with prior reconstructions for the past 2500 years. In total, these eruptions injected 7410 Tg S into the stratosphere: 70 % from tropical eruptions and 25 % from NH extratropical eruptions. A long-term latitudinally and monthly resolved stratospheric aerosol optical depth (SAOD) time series is reconstructed from the HolVol VSSI estimates, representing the first Holocene-scale reconstruction constrained by Greenland and Antarctica ice cores. These new long-term reconstructions of past VSSI and SAOD variability confirm evidence from regional volcanic eruption chronologies (e.g., from Iceland) in showing that the Early Holocene (9500–7000 BCE) experienced a higher number of volcanic eruptions (+16 %) and cumulative VSSI (+86 %) compared with the past 2500 years. This increase coincides with the rapid retreat of ice sheets during deglaciation, providing context for potential future increases in volcanic activity in regions under projected glacier melting in the 21st century. The reconstructed VSSI and SAOD data are available at https://doi.org/10.1594/PANGAEA.928646 (Sigl et al., 2021).

The 852/3 CE eruption of Mount Churchill, Alaska, was one of the largest first-millennium volcanic events, with a magnitude of 6.7 (VEI 6) and a tephra volume of 39.4–61.9 km3 (95 % confidence). The spatial extent of the ash fallout from this event is considerable and the cryptotephra (White River Ash east; WRAe) extends as far as Finland and Poland. Proximal ecosystem and societal disturbances have been linked with this eruption; however, wider eruption impacts on climate and society are unknown. Greenland ice core records show that the eruption occurred in winter 852/3 ± 1 CE and that the eruption is associated with a relatively moderate sulfate aerosol loading but large abundances of volcanic ash and chlorine. Here we assess the potential broader impact of this eruption using palaeoenvironmental reconstructions, historical records and climate model simulations. We also use the fortuitous timing of the 852/3 CE Churchill eruption and its extensively widespread tephra deposition of the White River Ash (east) (WRAe) to examine the climatic expression of the warm Medieval Climate Anomaly period (MCA; ca. 950–1250 CE) from precisely linked peatlands in the North Atlantic region. The reconstructed climate forcing potential of the 852/3 CE Churchill eruption is moderate compared with the eruption magnitude, but tree-ring-inferred temperatures report a significant atmospheric cooling of 0.8 ∘C in summer 853 CE. Modelled climate scenarios also show a cooling in 853 CE, although the average magnitude of cooling is smaller (0.3 ∘C). The simulated spatial patterns of cooling are generally similar to those generated using the tree-ring-inferred temperature reconstructions. Tree-ring-inferred cooling begins prior to the date of the eruption suggesting that natural internal climate variability may have increased the climate system's susceptibility to further cooling. The magnitude of the reconstructed cooling could also suggest that the climate forcing potential of this eruption may be underestimated, thereby highlighting the need for greater insight into, and consideration of, the role of halogens and volcanic ash when estimating eruption climate forcing potential. Precise comparisons of palaeoenvironmental records from peatlands across North America and Europe, facilitated by the presence of the WRAe isochron, reveal no consistent MCA signal. These findings contribute to the growing body of evidence that characterises the MCA hydroclimate as time-transgressive and heterogeneous rather than a well-defined climatic period. The presence of the WRAe isochron also demonstrates that no long-term (multidecadal) climatic or societal impacts from the 852/3 CE Churchill eruption were identified beyond areas proximal to the eruption. Historical evidence in Europe for subsistence crises demonstrate a degree of temporal correspondence on interannual timescales, but similar events were reported outside of the eruption period and were common in the 9th century. The 852/3 CE Churchill eruption exemplifies the difficulties of identifying and confirming volcanic impacts for a single eruption, even when the eruption has a small age uncertainty.

The mid-17th century is characterized by a cluster of explosive volcanic eruptions in the 1630s and 1640s, climatic conditions culminating in the Maunder Minimum, and political instability and famine in regions of western and northern Europe as well as China and Japan. This contribution investigates the sources of the eruptions of the 1630s and 1640s and their possible impact on contemporary climate using ice core, tree-ring, and historical evidence but will also look into the socio-political context in which they occurred and the human responses they may have triggered. Three distinct sulfur peaks are found in the Greenland ice core record in 1637, 1641–1642, and 1646. In Antarctica, only one unambiguous sulfate spike is recorded, peaking in 1642. The resulting bipolar sulfur peak in 1641–1642 can likely be ascribed to the eruption of Mount Parker (6∘ N, Philippines) on 26 December 1640, but sulfate emitted from Komaga-take (42∘ N, Japan) volcano on 31 July 1641 has potentially also contributed to the sulfate concentrations observed in Greenland at this time. The smaller peaks in 1637 and 1646 can be potentially attributed to the eruptions of Hekla (63∘ N, Iceland) and Shiveluch (56∘ N, Russia), respectively. To date, however, none of the candidate volcanoes for the mid-17th century sulfate peaks have been confirmed with tephra preserved in ice cores. Tree-ring and written sources point to cold conditions in the late 1630s and early 1640s in various parts of Europe and to poor harvests. Yet the early 17th century was also characterized by widespread warfare across Europe – and in particular the Thirty Years' War (1618–1648) – rendering any attribution of socio-economic crisis to volcanism challenging. In China and Japan, historical sources point to extreme droughts and famines starting in 1638 (China) and 1640 (Japan), thereby preceding the eruptions of Komaga-take (31 July 1640) and Mount Parker (4 January 1641). The case of the eruption cluster between 1637 and 1646 and the climatic and societal conditions recorded in its aftermath thus offer a textbook example of difficulties in (i) unambiguously distinguishing volcanically induced cooling, wetting, or drying from natural climate variability and (ii) attributing political instability, harvest failure, and famines solely to volcanic climatic impacts. This example shows that while the impacts of past volcanism must always be studied within the contemporary socio-economic contexts, it is also time to move past reductive framings and sometimes reactionary oppositional stances in which climate (and environment more broadly) either is or is not deemed an important contributor to major historical events.

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.

Volcanic fallout in polar ice sheets provides important opportunities to date and correlate ice-core records as well as to investigate the environmental impacts of eruptions. Only the geochemical characterization of volcanic ash (tephra) embedded in the ice strata can confirm the source of the eruption, however, and is a requisite if historical eruption ages are to be used as valid chronological checks on annual ice layer counting. Here we report the investigation of ash particles in a Greenland ice core that are associated with a volcanic sulfuric acid layer previously attributed to the 79 CE eruption of Vesuvius. Major and trace element composition of the particles indicates that the tephra does not derive from Vesuvius but most likely originates from an unidentified eruption in the Aleutian arc. Using ash dispersal modeling, we find that only an eruption large enough to include stratospheric injection is likely to account for the sizable (24–85 µm) ash particles observed in the Greenland ice at this time. Despite its likely explosivity, this event does not appear to have triggered significant climate perturbations, unlike some other large extratropical eruptions. In light of a recent re-evaluation of the Greenland ice-core chronologies, our findings further challenge the previous assignation of this volcanic event to 79 CE. We highlight the need for the revised Common Era ice-core chronology to be formally accepted by the wider ice-core and climate modeling communities in order to ensure robust age linkages to precisely dated historical and paleoclimate proxy records.

We present a transient simulation of global vegetation and climate patterns of the mid- and late Holocene using the MPI-ESM (Max Planck Institute for Meteorology Earth System Model) at T63 resolution. The simulated vegetation trend is discussed in the context of the simulated Holocene climate change. Our model captures the main trends found in reconstructions. Most prominent are the southward retreat of the northern treeline that is combined with the strong decrease of forest in the high northern latitudes during the Holocene and the vast increase of the Saharan desert, embedded in a general decrease in precipitation and vegetation in the Northern Hemisphere monsoon margin regions. The Southern Hemisphere experiences weaker changes in total vegetation cover during the last 8000 years. However, the monsoon-related increase in precipitation and the insolation-induced cooling of the winter climate lead to shifts in the vegetation composition, mainly between the woody plant functional types (PFTs). The large-scale global patterns of vegetation almost linearly follow the subtle, approximately linear, orbital forcing. In some regions, however, non-linear, more rapid changes in vegetation are found in the simulation. The most striking region is the Sahel–Sahara domain with rapid vegetation transitions to a rather desertic state, despite a gradual insolation forcing. Rapid shifts in the simulated vegetation also occur in the high northern latitudes, in South Asia and in the monsoon margins of the Southern Hemisphere. These rapid changes are mainly triggered by changes in the winter temperatures, which go into, or move out of, the bioclimatic tolerance range of individual PFTs. The dynamics of the transitions are determined by dynamics of the net primary production (NPP) and the competition between PFTs. These changes mainly occur on timescales of centuries. More rapid changes in PFTs that occur within a few decades are mainly associated with the timescales of mortality and the bioclimatic thresholds implicit in the dynamic vegetation model, which have to be interpreted with caution. Most of the simulated Holocene vegetation changes outside the high northern latitudes are associated with modifications in the intensity of the global summer monsoon dynamics that also affect the circulation in the extra tropics via teleconnections. Based on our simulations, we thus identify the global monsoons as the key player in Holocene climate and vegetation change.

In this study, we analyze extreme daily precipitation during the pre-industrial period from 1501 BCE to 1849 CE in simulations from the Community Earth System Model version 1.2.2. A peak-over-threshold (POT) extreme value analysis is employed to examine characteristics of extreme precipitation and to identify connections of extreme precipitation with the external forcing and with modes of internal variability. The POT analysis shows that extreme precipitation with similar statistical characteristics, i.e., the probability density distributions, tends to cluster spatially. There are differences in the distribution of extreme precipitation between the Pacific and Atlantic sectors and between the northern high and southern low latitudes. Extreme precipitation during the pre-industrial period is largely influenced by modes of internal variability, such as El Niño–Southern Oscillation (ENSO), the Pacific North American, and Pacific South American patterns, among others, and regional surface temperatures. In general, the modes of variability exhibit a statistically significant connection to extreme precipitation in the vicinity to their regions of action. The exception is ENSO, which shows more widespread influence on extreme precipitation across the Earth. In addition, the regions with which extreme precipitation is more associated, either by a mode of variability or by the regional surface temperature, are distinguished. Regional surface temperatures are associated with extreme precipitation over lands at the extratropical latitudes and over the tropical oceans. In other regions, the influence of modes of variability is still dominant. Effects of the changes in the orbital parameters on extreme precipitation are rather weak compared to those of the modes of internal variability and of the regional surface temperatures. Still, some regions in central Africa, southern Asia, and the tropical Atlantic ocean show statistically significant connections between extreme precipitation and orbital forcing, implying that in these regions, extreme precipitation has increased linearly during the 3351-year pre-industrial period. Tropical volcanic eruptions affect extreme precipitation more clearly in the short term up to a few years, altering both the intensity and frequency of extreme precipitation. However, more apparent changes are found in the frequency than the intensity of extreme precipitation. After eruptions, the return periods of extreme precipitation increase over the extratropical regions and the tropical Pacific, while a decrease is found in other regions. The post-eruption changes in the frequency of extreme precipitation are associated with ENSO, which itself is influenced by tropical eruptions. Overall, the results show that climate simulations are useful to complement the information on pre-industrial extreme precipitation, as they elucidate statistical characteristics and long-term connections of extreme events with natural variability.

Volcanic eruptions are a key source of climatic variability, and reconstructing their past impact can improve our understanding of the operation of the climate system and increase the accuracy of future climate projections. Two annually resolved and independently dated palaeoarchives – tree rings and polar ice cores – can be used in tandem to assess the timing, strength and climatic impact of volcanic eruptions over the past ∼ 2500 years. The quantification of post-volcanic climate responses, however, has at times been hampered by differences between simulated and observed temperature responses that raised questions regarding the robustness of the chronologies of both archives. While many chronological mismatches have been resolved, the precise timing and climatic impact of two major sulfate-emitting volcanic eruptions during the 1450s CE, including the largest atmospheric sulfate-loading event in the last 700 years, have not been constrained. Here we explore this issue through a combination of tephrochronological evidence and high-resolution ice-core chemistry measurements from a Greenland ice core, the TUNU2013 record. We identify tephra from the historically dated 1477 CE eruption of the Icelandic Veiðivötn–Bárðarbunga volcanic system in direct association with a notable sulfate peak in TUNU2013 attributed to this event, confirming that this peak can be used as a reliable and precise time marker. Using seasonal cycles in several chemical elements and 1477 CE as a fixed chronological point shows that ages of 1453 CE and 1458 CE can be attributed, with high precision, to the start of two other notable sulfate peaks. This confirms the accuracy of a recent Greenland ice-core chronology over the middle to late 15th century and corroborates the findings of recent volcanic reconstructions from Greenland and Antarctica. Overall, this implies that large-scale Northern Hemisphere climatic cooling affecting tree-ring growth in 1453 CE was caused by a Northern Hemisphere volcanic eruption in 1452 or early 1453 CE, and then a Southern Hemisphere eruption, previously assumed to have triggered the cooling, occurred later in 1457 or 1458 CE. The direct attribution of the 1477 CE sulfate peak to the eruption of Veiðivötn, one of the most explosive from Iceland in the last 1200 years, also provides the opportunity to assess the eruption's climatic impact. A tree-ring-based reconstruction of Northern Hemisphere summer temperatures shows a cooling in the aftermath of the eruption of −0.35 ∘C relative to a 1961–1990 CE reference period and −0.1 ∘C relative to the 30-year period around the event, as well as a relatively weak and spatially incoherent climatic response in comparison to the less explosive but longer-lasting Icelandic Eldgjá 939 CE and Laki 1783 CE eruptions. In addition, the Veiðivötn 1477 CE eruption occurred around the inception of the Little Ice Age and could be used as a chronostratigraphic marker to constrain the phasing and spatial variability of climate changes over this transition if it can be traced in more regional palaeoclimatic archives.

Recent large wildfires, such as those in Portugal in 2017, have devastating impacts on societies, economy, ecosystems and environments. However, wildfires are a natural phenomenon, which has been exacerbated by land use during the past millennia. Ice cores are one of the archives preserving information on fire occurrences over these timescales. A difficulty is that emission sensitivity of ice cores is often unknown, which constitutes a source of uncertainty in the interpretation of such archives. Information from specific and well-documented case studies is therefore useful to better understand the spatial representation of ice-core burning records. The wildfires near Pedrógão Grande in central Portugal in 2017 provided a test bed to link a fire event to its footprint left in a high-alpine snowpack considered a surrogate for high-alpine ice-core sites. Here, we (1) analysed black carbon (BC) and microscopic charcoal particles deposited in the snowpack close to the high-alpine research station Jungfraujoch in the Swiss Alps; (2) calculated backward trajectories based on ERA-Interim reanalysis data and simulated the transport of these carbonaceous particles using a global aerosol-climate model; and (3) analysed the fire spread, its spatial and temporal extent, and its intensity with remote-sensing (e.g. MODIS) Active Fire and Burned Area products. According to modelled emissions of the FINN v1.6 database, the fire emitted a total amount of 203.5 t BC from a total burned area of 501 km2 as observed on the basis of satellite fire products. Backward trajectories unambiguously linked a peak of atmospheric-equivalent BC observed at the Jungfraujoch research station on 22 June – with elevated levels until 25 June – with the highly intensive fires in Portugal. The atmospheric signal is in correspondence with an outstanding peak in microscopic charcoal observed in the snow layer, depositing nearly as many charcoal particles as during an average year in other ice archives. In contrast to charcoal, the amount of atmospheric BC deposited during the fire episode was minor due to a lack of precipitation. Simulations with a global aerosol-climate model suggest that the observed microscopic charcoal particles originated from the fires in Portugal and that their contribution to the BC signal in snow was negligible. Our study revealed that microscopic charcoal can be transported over long distances (1500 km) and that snow and ice archives are much more sensitive to distant events than sedimentary archives, for which the signal is dominated by local fires. The findings are important for future ice-core studies as they document that, for BC as a fire tracer, the signal preservation depends on precipitation. Single events, like this example, might not be preserved due to unfavourable meteorological conditions.

The last glacial period is characterized by a number of millennial climate events that have been identified in both Greenland and Antarctic ice cores and that are abrupt in Greenland climate records. The mechanisms governing this climate variability remain a puzzle that requires a precise synchronization of ice cores from the two hemispheres to be resolved. Previously, Greenland and Antarctic ice cores have been synchronized primarily via their common records of gas concentrations or isotopes from the trapped air and via cosmogenic isotopes measured on the ice. In this work, we apply ice core volcanic proxies and annual layer counting to identify large volcanic eruptions that have left a signature in both Greenland and Antarctica. Generally, no tephra is associated with those eruptions in the ice cores, so the source of the eruptions cannot be identified. Instead, we identify and match sequences of volcanic eruptions with bipolar distribution of sulfate, i.e. unique patterns of volcanic events separated by the same number of years at the two poles. Using this approach, we pinpoint 82 large bipolar volcanic eruptions throughout the second half of the last glacial period (12–60 ka). This improved ice core synchronization is applied to determine the bipolar phasing of abrupt climate change events at decadal-scale precision. In response to Greenland abrupt climatic transitions, we find a response in the Antarctic water isotope signals (δ18O and deuterium excess) that is both more immediate and more abrupt than that found with previous gas-based interpolar synchronizations, providing additional support for our volcanic framework. On average, the Antarctic bipolar seesaw climate response lags the midpoint of Greenland abrupt δ18O transitions by 122±24 years. The time difference between Antarctic signals in deuterium excess and δ18O, which likewise informs the time needed to propagate the signal as described by the theory of the bipolar seesaw but is less sensitive to synchronization errors, suggests an Antarctic δ18O lag behind Greenland of 152±37 years. These estimates are shorter than the 200 years suggested by earlier gas-based synchronizations. As before, we find variations in the timing and duration between the response at different sites and for different events suggesting an interaction of oceanic and atmospheric teleconnection patterns as well as internal climate variability.

The South Pole Ice Core (SPICEcore) was drilled in 2014–2016 to provide a detailed multi-proxy archive of paleoclimate conditions in East Antarctica during the Holocene and late Pleistocene. Interpretation of these records requires an accurate depth–age relationship. Here, we present the SPICEcore (SP19) timescale for the age of the ice of SPICEcore. SP19 is synchronized to the WD2014 chronology from the West Antarctic Ice Sheet Divide (WAIS Divide) ice core using stratigraphic matching of 251 volcanic events. These events indicate an age of 54 302±519 BP (years before 1950) at the bottom of SPICEcore. Annual layers identified in sodium and magnesium ions to 11 341 BP were used to interpolate between stratigraphic volcanic tie points, yielding an annually resolved chronology through the Holocene. Estimated timescale uncertainty during the Holocene is less than 18 years relative to WD2014, with the exception of the interval between 1800 to 3100 BP when uncertainty estimates reach ±25 years due to widely spaced volcanic tie points. Prior to the Holocene, uncertainties remain within 124 years relative to WD2014. Results show an average Holocene accumulation rate of 7.4 cm yr−1 (water equivalent). The time variability of accumulation rate is consistent with expectations for steady-state ice flow through the modern spatial pattern of accumulation rate. Time variations in nitrate concentration, nitrate seasonal amplitude and δ15N of N2 in turn are as expected for the accumulation rate variations. The highly variable yet well-constrained Holocene accumulation history at the site can help improve scientific understanding of deposition-sensitive climate proxies such as δ15N of N2 and photolyzed chemical compounds.

The Amazon Basin is one of the major contributors to global biomass burning emissions. However, regional paleofire trends remain particularly unknown. Due to their proximity to the Amazon Basin, Andean ice cores are suitable to reconstruct paleofire trends in South America and improve our understanding of the complex linkages between fires, climate and humans. Here we present the first refractory black carbon (rBC) ice-core record from the Andes as a proxy for biomass burning emissions in the Amazon Basin, derived from an ice core drilled at 6300 m a.s.l. from the Illimani glacier in the Bolivian Andes and spanning the entire Holocene back to the last deglaciation 13 000 years ago. The Illimani rBC record displays a strong seasonality with low values during the wet season and high values during the dry season due to the combination of enhanced biomass burning emissions in the Amazon Basin and less precipitation at the Illimani site. Significant positive (negative) correlations were found with reanalyzed temperature (precipitation) data for regions in eastern Bolivia and western Brazil characterized by substantial fire activity. rBC long-term trends indirectly reflect regional climatic variations through changing biomass burning emissions as they show higher (lower) concentrations during warm–dry (cold–wet) periods, in line with climate variations such as the Younger Dryas, the 8.2 ka event, the Holocene Climatic Optimum, the Medieval Warm Period and the Little Ice Age. The highest rBC concentrations of the entire record occurred during the Holocene Climatic Optimum between 7000 and 3000 BCE, suggesting that this exceptionally warm and dry period caused high levels of biomass burning activity, unprecedented in the context of the past 13 000 years. Recent rBC levels, rising since 1730 CE in the context of increasing temperatures and deforestation, are similar to those of the Medieval Warm Period. No decrease in fire activity was observed in the 20th century, in contradiction to global biomass burning reconstructions based on charcoal data.

Light absorbing aerosols in the atmosphere and cryosphere play an important role in the climate system. Their presence in ambient air and snow changes the radiative properties of these systems, thus contributing to increased atmospheric warming and snowmelt. High spatio-temporal variability of aerosol concentrations and a shortage of long-term observations contribute to large uncertainties in properly assigning the climate effects of aerosols through time.Starting around AD 1860, many glaciers in the European Alps began to retreat from their maximum mid-19th century terminus positions, thereby visualizing the end of the Little Ice Age in Europe. Radiative forcing by increasing deposition of industrial black carbon to snow has been suggested as the main driver of the abrupt glacier retreats in the Alps. The basis for this hypothesis was model simulations using elemental carbon concentrations at low temporal resolution from two ice cores in the Alps.Here we present sub-annually resolved concentration records of refractory black carbon (rBC; using soot photometry) as well as distinctive tracers for mineral dust, biomass burning and industrial pollution from the Colle Gnifetti ice core in the Alps from AD 1741 to 2015. These records allow precise assessment of a potential relation between the timing of observed acceleration of glacier melt in the mid-19th century with an increase of rBC deposition on the glacier caused by the industrialization of Western Europe. Our study reveals that in AD 1875, the time when rBC ice-core concentrations started to significantly increase, the majority of Alpine glaciers had already experienced more than 80 % of their total 19th century length reduction, casting doubt on a leading role for soot in terminating of the Little Ice Age. Attribution of glacial retreat requires expansion of the spatial network and sampling density of high alpine ice cores to balance potential biasing effects arising from transport, deposition, and snow conservation in individual ice-core records.

Produced by the incomplete combustion of fossil fuel and biomass, black carbon (BC) contributes to Arctic warming by reducing snow albedo and thus triggering a snow-albedo feedback leading to increased snowmelt. Therefore, it is of high importance to assess past BC emissions to better understand and constrain their role. However, only a few long-term BC records are available from the Arctic, mainly originating from Greenland ice cores. Here, we present the first long-term and high-resolution refractory black carbon (rBC) record from Svalbard, derived from the analysis of two ice cores drilled at the Lomonosovfonna ice field in 2009 (LF-09) and 2011 (LF-11) and covering 800 years of atmospheric emissions. Our results show that rBC concentrations strongly increased from 1860 on due to anthropogenic emissions and reached two maxima, at the end of the 19th century and in the middle of the 20th century. No increase in rBC concentrations during the last decades was observed, which is corroborated by atmospheric measurements elsewhere in the Arctic but contradicts a previous study from another ice core from Svalbard. While melting may affect BC concentrations during periods of high temperatures, rBC concentrations remain well preserved prior to the 20th century due to lower temperatures inducing little melt. Therefore, the preindustrial rBC record (before 1800), along with ammonium (NH4+), formate (HCOO−) and specific organic markers (vanillic acid, VA, and p-hydroxybenzoic acid, p-HBA), was used as a proxy for biomass burning. Despite numerous single events, no long-term trend was observed over the time period 1222–1800 for rBC and NH4+. In contrast, formate, VA, and p-HBA experience multi-decadal peaks reflecting periods of enhanced biomass burning. Most of the background variations and single peak events are corroborated by other ice core records from Greenland and Siberia. We suggest that the paleofire record from the LF ice core primarily reflects biomass burning episodes from northern Eurasia, induced by decadal-scale climatic variations.

Microfossil records from ice archives allow vegetation, fire and land-use activity reconstructions on broad spatial scales. Samples typically contain low microfossil concentrations. Therefore, large ice volumes are often needed for palynology. Hence, it is crucial to extract maximum microfossil numbers through appropriate physical-chemical treatments. We compare six methods covering the main water reduction procedures: evaporation, filtration and centrifugation with snow samples. Adding a known number of Lycopodium marker spores prior to sample treatment and a second marker (Eucalyptus) after laboratory processing allows a quantitative microfossil loss assessment during pollen extraction. We applied the best-performing method (average loss of 22%) to high-alpine firn cores from Colle Gnifetti glacier for validation with a natural archive containing extremely low microfossil concentrations. We conclude that samples processed with different microfossil extraction protocols may give different results for pollen concentrations, percentages and ratios between different pollen types, especially if vesiculate conifer pollen is an important pollen assemblage component. We recommend a new evaporation-based method which delivers the smallest and least variable losses among the tested approaches. Since microfossil losses are inevitable during laboratory procedure, adding markers prior to sample processing is mandatory to achieve reliable microfossil concentration and influx estimates.

The eruption of Mt. Tambora in 1815 was the largest volcanic eruption of the past 500 years. The eruption had significant climatic impacts, leading to the 1816 year without a summer, and remains a valuable event from which to understand the climatic effects of large stratospheric volcanic sulfur dioxide injections. The eruption also resulted in one of the strongest and most easily identifiable volcanic sulfate signals in polar ice cores, which are widely used to reconstruct the timing and atmospheric sulfate loading of past eruptions. As part of the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP), five state-of-the-art global aerosol models simulated this eruption. We analyse both simulated background (no Tambora) and volcanic (with Tambora) sulfate deposition to polar regions and compare to ice core records. The models simulate overall similar patterns of background sulfate deposition, although there are differences in regional details and magnitude. However, the volcanic sulfate deposition varies considerably between the models with differences in timing, spatial pattern and magnitude. Mean simulated deposited sulfate on Antarctica ranges from 19 to 264 kg km−2 and on Greenland from 31 to 194 kg km−2, as compared to the mean ice-core-derived estimates of roughly 50 kg km−2 for both Greenland and Antarctica. The ratio of the hemispheric atmospheric sulfate aerosol burden after the eruption to the average ice sheet deposited sulfate varies between models by up to a factor of 15. Sources of this inter-model variability include differences in both the formation and the transport of sulfate aerosol. Our results suggest that deriving relationships between sulfate deposited on ice sheets and atmospheric sulfate burdens from model simulations may be associated with greater uncertainties than previously thought.

The injection of sulfur into the stratosphere by explosive volcanic eruptions is the cause of significant climate variability. Based on sulfate records from a suite of ice cores from Greenland and Antarctica, the eVolv2k database includes estimates of the magnitudes and approximate source latitudes of major volcanic stratospheric sulfur injection (VSSI) events from 500 BCE to 1900 CE, constituting an update of prior reconstructions and an extension of the record by 1000 years. The database incorporates improvements to the ice core records (in terms of synchronisation and dating) and refinements to the methods used to estimate VSSI from ice core records, and it includes first estimates of the random uncertainties in VSSI values. VSSI estimates for many of the largest eruptions, including Samalas (1257), Tambora (1815), and Laki (1783), are within 10 % of prior estimates. A number of strong events are included in eVolv2k which are largely underestimated or not included in earlier VSSI reconstructions, including events in 540, 574, 682, and 1108 CE. The long-term annual mean VSSI from major volcanic eruptions is estimated to be ∼ 0.5 Tg [S] yr−1, ∼ 50 % greater than a prior reconstruction due to the identification of more events and an increase in the magnitude of many intermediate events. A long-term latitudinally and monthly resolved stratospheric aerosol optical depth (SAOD) time series is reconstructed from the eVolv2k VSSI estimates, and the resulting global mean SAOD is found to be similar (within 33 %) to a prior reconstruction for most of the largest eruptions. The long-term (500 BCE–1900 CE) average global mean SAOD estimated from the eVolv2k VSSI estimates including a constant background injection of stratospheric sulfur is ∼ 0.014, 30 % greater than a prior reconstruction. These new long-term reconstructions of past VSSI and SAOD variability give context to recent volcanic forcing, suggesting that the 20th century was a period of somewhat weaker than average volcanic forcing, with current best estimates of 20th century mean VSSI and SAOD values being 25 and 14 % less, respectively, than the mean of the 500 BCE to 1900 CE period. The reconstructed VSSI and SAOD data are available at https://doi.org/10.1594/WDCC/eVolv2k_v2.

The pre-industrial millennium is among the periods selected by the Paleoclimate Model Intercomparison Project (PMIP) for experiments contributing to the sixth phase of the Coupled Model Intercomparison Project (CMIP6) and the fourth phase of the PMIP (PMIP4). The past1000 transient simulations serve to investigate the response to (mainly) natural forcing under background conditions not too different from today, and to discriminate between forced and internally generated variability on interannual to centennial timescales. This paper describes the motivation and the experimental set-ups for the PMIP4-CMIP6 past1000 simulations, and discusses the forcing agents orbital, solar, volcanic, and land use/land cover changes, and variations in greenhouse gas concentrations. The past1000 simulations covering the pre-industrial millennium from 850 Common Era (CE) to 1849 CE have to be complemented by historical simulations (1850 to 2014 CE) following the CMIP6 protocol. The external forcings for the past1000 experiments have been adapted to provide a seamless transition across these time periods. Protocols for the past1000 simulations have been divided into three tiers. A default forcing data set has been defined for the Tier 1 (the CMIP6 past1000) experiment. However, the PMIP community has maintained the flexibility to conduct coordinated sensitivity experiments to explore uncertainty in forcing reconstructions as well as parameter uncertainty in dedicated Tier 2 simulations. Additional experiments (Tier 3) are defined to foster collaborative model experiments focusing on the early instrumental period and to extend the temporal range and the scope of the simulations. This paper outlines current and future research foci and common analyses for collaborative work between the PMIP and the observational communities (reconstructions, instrumental data).

Englacial ice contains a significant reservoir of organic material (OM), preserving a chronological record of materials from Earth's past. Here, we investigate if OM composition surveys in ice core research can provide paleoecological information on the dynamic nature of our Earth through time. Temporal trends in OM composition from the early Holocene extending back to the Last Glacial Maximum (LGM) of the West Antarctic Ice Sheet Divide (WD) ice core were measured by fluorescence spectroscopy. Multivariate parallel factor (PARAFAC) analysis is widely used to isolate the chemical components that best describe the observed variation across three-dimensional fluorescence spectroscopy (excitation–emission matrices; EEMs) assays. Fluorescent OM markers identified by PARAFAC modeling of the EEMs from the LGM (27.0–18.0 kyr BP; before present 1950) through the last deglaciation (LD; 18.0–11.5 kyr BP), to the mid-Holocene (11.5–6.0 kyr BP) provided evidence of different types of fluorescent OM composition and origin in the WD ice core over 21.0 kyr. Low excitation–emission wavelength fluorescent PARAFAC component one (C1), associated with chemical species similar to simple lignin phenols was the greatest contributor throughout the ice core, suggesting a strong signature of terrestrial OM in all climate periods. The component two (C2) OM marker, encompassed distinct variability in the ice core describing chemical species similar to tannin- and phenylalanine-like material. Component three (C3), associated with humic-like terrestrial material further resistant to biodegradation, was only characteristic of the Holocene, suggesting that more complex organic polymers such as lignins or tannins may be an ecological marker of warmer climates. We suggest that fluorescent OM markers observed during the LGM were the result of greater continental dust loading of lignin precursor (monolignol) material in a drier climate, with lower marine influences when sea ice extent was higher and continents had more expansive tundra cover. As the climate warmed, the record of OM markers in the WD ice core changed, reflecting shifts in carbon productivity as a result of global ecosystem response.

Wildfires and their emissions have significant impacts on ecosystems, climate, atmospheric chemistry, and carbon cycling. Well-dated proxy records are needed to study the long-term climatic controls on biomass burning and the associated climate feedbacks. There is a particular lack of information about long-term biomass burning variations in Siberia, the largest forested area in the Northern Hemisphere. In this study we report analyses of aromatic acids (vanillic and para-hydroxybenzoic acids) over the past 2600 years in the Eurasian Arctic Akademii Nauk ice core. These compounds are aerosol-borne, semi-volatile organic compounds derived from lignin combustion. The analyses were made using ion chromatography with electrospray mass spectrometric detection. The levels of these aromatic acids ranged from below the detection limit (0.01 to 0.05 ppb; 1 ppb = 1000 ng L−1) to about 1 ppb, with roughly 30 % of the samples above the detection limit. In the preindustrial late Holocene, highly elevated aromatic acid levels are observed during three distinct periods (650–300 BCE, 340–660 CE, and 1460–1660 CE). The timing of the two most recent periods coincides with the episodic pulsing of ice-rafted debris in the North Atlantic known as Bond events and a weakened Asian monsoon, suggesting a link between fires and large-scale climate variability on millennial timescales. Aromatic acid levels also are elevated during the onset of the industrial period from 1780 to 1860 CE, but with a different ratio of vanillic and para-hydroxybenzoic acid than is observed during the preindustrial period. This study provides the first millennial-scale record of aromatic acids. This study clearly demonstrates that coherent aromatic acid signals are recorded in polar ice cores that can be used as proxies for past trends in biomass burning.

Reconstruction of past changes in Arctic sea ice extent may be critical for understanding its future evolution. Methanesulfonate (MSA) and bromine concentrations preserved in ice cores have both been proposed as indicators of past sea ice conditions. In this study, two ice cores from central and north-eastern Greenland were analysed at sub-annual resolution for MSA (CH3SO3H) and bromine, covering the time period 1750–2010. We examine correlations between ice core MSA and the HadISST1 ICE sea ice dataset and consult back trajectories to infer the likely source regions. A strong correlation between the low-frequency MSA and bromine records during pre-industrial times indicates that both chemical species are likely linked to processes occurring on or near sea ice in the same source regions. The positive correlation between ice core MSA and bromine persists until the mid-20th century, when the acidity of Greenland ice begins to increase markedly due to increased fossil fuel emissions. After that time, MSA levels decrease as a result of declining sea ice extent but bromine levels increase. We consider several possible explanations and ultimately suggest that increased acidity, specifically nitric acid, of snow on sea ice stimulates the release of reactive Br from sea ice, resulting in increased transport and deposition on the Greenland ice sheet.

High-altitude glaciers and ice caps from midlatitudes and tropical regions contain valuable signals of past climatic and environmental conditions as well as human activities, but for a meaningful interpretation this information needs to be placed in a precise chronological context. For dating the upper part of ice cores from such sites, several relatively precise methods exist, but they fail in the older and deeper parts, where plastic deformation of the ice results in strong annual layer thinning and a non-linear age–depth relationship. If sufficient organic matter such as plant, wood or insect fragments were found, radiocarbon (14C) analysis would have thus been the only option for a direct and absolute dating of deeper ice core sections. However such fragments are rarely found and, even then, they would not be very likely to occur at the desired depth and resolution. About 10 years ago, a new, complementary dating tool was therefore introduced by our group. It is based on extracting the µg-amounts of the water-insoluble organic carbon (WIOC) fraction of carbonaceous aerosols embedded in the ice matrix for subsequent 14C dating. Since then this new approach has been improved considerably by reducing the measurement time and improving the overall precision. Samples with ∼ 10 µg WIOC mass can now be dated with reasonable uncertainty of around 10–20 % (variable depending on sample age). This requires about 300 to 800 g of ice for WIOC concentrations typically found in midlatitude and low-latitude glacier ice. Dating polar ice with satisfactory age precision is still not possible since WIOC concentrations are around 1 order of magnitude lower. The accuracy of the WIOC 14C method was validated by applying it to independently dated ice. With this method, the deepest parts of the ice cores from Colle Gnifetti and the Mt Ortles glacier in the European Alps, Illimani glacier in the Bolivian Andes, Tsambagarav ice cap in the Mongolian Altai, and Belukha glacier in the Siberian Altai have been dated. In all cases a strong annual layer thinning towards the bedrock was observed and the oldest ages obtained were in the range of 10 000 years. WIOC 14C dating was not only crucial for interpretation of the embedded environmental and climatic histories, but additionally gave a better insight into glacier flow dynamics close to the bedrock and past glacier coverage. For this the availability of multiple dating points in the deepest parts was essential, which is the strength of the presented WIOC 14C dating method, allowing determination of absolute ages from principally every piece of ice.

Here, we review different attempts made since the early 1990s to reconstruct past forest fire activity using chemical signals recorded in ice cores extracted from the Greenland ice sheet and a few mid-northern latitude, high-elevation glaciers. We first examined the quality of various inorganic (ammonium, nitrate, potassium) and organic (black carbon, various organic carbon compounds including levoglucosan and numerous carboxylic acids) species proposed as fire proxies in ice, particularly in Greenland. We discuss limitations in their use during recent vs. pre-industrial times, atmospheric lifetimes, and the relative importance of other non-biomass-burning sources. Different high-resolution records from several Greenland drill sites and covering various timescales, including the last century and Holocene, are discussed. We explore the extent to which atmospheric transport can modulate the record of boreal fires from Canada as recorded in Greenland ice. Ammonium, organic fractions (black and organic carbon), and specific organic compounds such as formate and vanillic acid are found to be good proxies for tracing past boreal fires in Greenland ice. We show that use of other species – potassium, nitrate, and carboxylates (except formate) – is complicated by either post-depositional effects or existence of large non-biomass-burning sources. The quality of levoglucosan with respect to other proxies is not addressed here because of a lack of high-resolution profiles for this species, preventing a fair comparison. Several Greenland ice records of ammonium consistently indicate changing fire activity in Canada in response to past climatic conditions that occurred during the last millennium and since the last large climatic transition. Based on this review, we make recommendations for further study to increase reliability of the reconstructed history of forest fires occurring in a given region.

Advances in trace gas analysis allow localised, non-atmospheric features to be resolved in ice cores, superimposed on the coherent atmospheric signal. These high-frequency signals could not have survived the low-pass filter effect that gas diffusion in the firn exerts on the atmospheric history and therefore do not result from changes in the atmospheric composition at the ice sheet surface. Using continuous methane (CH4) records obtained from five polar ice cores, we characterise these non-atmospheric signals and explore their origin. Isolated samples, enriched in CH4 in the Tunu13 (Greenland) record are linked to the presence of melt layers. Melting can enrich the methane concentration due to a solubility effect, but we find that an additional in situ process is required to generate the full magnitude of these anomalies. Furthermore, in all the ice cores studied there is evidence of reproducible, decimetre-scale CH4 variability. Through a series of tests, we demonstrate that this is an artifact of layered bubble trapping in a heterogeneous-density firn column; we use the term “trapping signal” for this phenomenon. The peak-to-peak amplitude of the trapping signal is typically 5 ppb, but may exceed 40 ppb. Signal magnitude increases with atmospheric CH4 growth rate and seasonal density contrast, and decreases with accumulation rate. Significant annual periodicity is present in the CH4 variability of two Greenland ice cores, suggesting that layered gas trapping at these sites is controlled by regular, seasonal variations in the physical properties of the firn. Future analytical campaigns should anticipate high-frequency artifacts at high-melt ice core sites or during time periods with high atmospheric CH4 growth rate in order to avoid misinterpretation of such features as past changes in atmospheric composition.

We present the WD2014 chronology for the upper part (0–2850 m; 31.2 ka BP) of the West Antarctic Ice Sheet (WAIS) Divide (WD) ice core. The chronology is based on counting of annual layers observed in the chemical, dust and electrical conductivity records. These layers are caused by seasonal changes in the source, transport, and deposition of aerosols. The measurements were interpreted manually and with the aid of two automated methods. We validated the chronology by comparing to two high-accuracy, absolutely dated chronologies. For the Holocene, the cosmogenic isotope records of 10Be from WAIS Divide and 14C for IntCal13 demonstrated that WD2014 was consistently accurate to better than 0.5 % of the age. For the glacial period, comparisons to the Hulu Cave chronology demonstrated that WD2014 had an accuracy of better than 1 % of the age at three abrupt climate change events between 27 and 31 ka. WD2014 has consistently younger ages than Greenland ice core chronologies during most of the Holocene. For the Younger Dryas–Preboreal transition (11.595 ka; 24 years younger) and the Bølling–Allerød Warming (14.621 ka; 7 years younger), WD2014 ages are within the combined uncertainties of the timescales. Given its high accuracy, WD2014 can become a reference chronology for the Southern Hemisphere, with synchronization to other chronologies feasible using high-quality proxies of volcanism, solar activity, atmospheric mineral dust, and atmospheric methane concentrations.

The West Antarctic Ice Sheet Divide (WAIS Divide, WD) ice core is a newly drilled, high-accumulation deep ice core that provides Antarctic climate records of the past ∼68 ka at unprecedented temporal resolution. The upper 2850 m (back to 31.2 ka BP) have been dated using annual-layer counting. Here we present a chronology for the deep part of the core (67.8–31.2 ka BP), which is based on stratigraphic matching to annual-layer-counted Greenland ice cores using globally well-mixed atmospheric methane. We calculate the WD gas age–ice age difference (Δage) using a combination of firn densification modeling, ice-flow modeling, and a data set of δ15N-N2, a proxy for past firn column thickness. The largest Δage at WD occurs during the Last Glacial Maximum, and is 525 ± 120 years. Internally consistent solutions can be found only when assuming little to no influence of impurity content on densification rates, contrary to a recently proposed hypothesis. We synchronize the WD chronology to a linearly scaled version of the layer-counted Greenland Ice Core Chronology (GICC05), which brings the age of Dansgaard–Oeschger (DO) events into agreement with the U/Th absolutely dated Hulu Cave speleothem record. The small Δage at WD provides valuable opportunities to investigate the timing of atmospheric greenhouse gas variations relative to Antarctic climate, as well as the interhemispheric phasing of the "bipolar seesaw".