Decadal-to-centennial increases of volcanic aerosols from Iceland challenge the concept of a Medieval Quiet Period

Existing global volcanic radiative aerosol forcing estimates portray the period 700 to 1000 as volcanically quiescent, void of major volcanic eruptions. However, this disagrees with proximal Icelandic geological records and regional Greenland ice-core records of sulfate. Here, we use cryptotephra analyses, high-resolution sulfur isotope analyses, and glaciochemical volcanic tracers on an array of Greenland ice cores to characterise volcanic activity and climatically important sulfuric aerosols across the period 700 to 1000. We identify a prolonged episode of volcanic sulfur dioxide emissions (751 – 940) dominated by Icelandic volcanism, that we term the Icelandic Active Period. This period commences with the Hrafnkatla episode (751 – 763), which coincided with strong winter cooling anomalies across Europe. This study reveals an important contribution of prolonged volcanic sulfate emissions to the pre-industrial atmospheric aerosol burden, currently not considered in existing forcing estimates, and highlights the need for further research to disentangle their associated climate feedbacks.


Check for updates
Imogen Gabriel 1 , Gill Plunkett 2 , Peter M. Abbott 1 , Melanie Behrens 3 , Andrea Burke 4 , Nathan Chellman 5 , Eliza Cook 6 , Dominik Fleitmann 7 , Maria Hörhold 3 , William Hutchison 4 , Joseph R. McConnell 5 , Bergrún A. Óladóttir 8 , Johannes Preiser-Kapeller 9 , Jakub T. Sliwinski 4,11 , Patrick Sugden 4 , Birthe Twarloh 3 & Michael Sigl 1,10 Existing global volcanic radiative aerosol forcing estimates portray the period 700 to 1000 as volcanically quiescent, void of major volcanic eruptions.However, this disagrees with proximal Icelandic geological records and regional Greenland ice-core records of sulfate.Here, we use cryptotephra analyses, high-resolution sulfur isotope analyses, and glaciochemical volcanic tracers on an array of Greenland ice cores to characterise volcanic activity and climatically important sulfuric aerosols across the period 700 to 1000.We identify a prolonged episode of volcanic sulfur dioxide emissions (751-940) dominated by Icelandic volcanism, that we term the Icelandic Active Period.This period commences with the Hrafnkatla episode (751-763), which coincided with strong winter cooling anomalies across Europe.This study reveals an important contribution of prolonged volcanic sulfate emissions to the pre-industrial atmospheric aerosol burden, currently not considered in existing forcing estimates, and highlights the need for further research to disentangle their associated climate feedbacks.
Volcanic aerosols are considered major drivers of natural climate variability, modifying the climate system through the reflection of incoming solar radiation, and subsequent cooling of the Earth's surface 1,2 .As a result, reconstructions of volcanic aerosol radiative forcing underpin climate model simulations, and the generation of accurate forcing estimates are therefore essential 3,4 .
The primary archives used for such reconstructions are the polar ice-core records.These highly resolved records provide insights into past volcanic activity over longer time periods (e.g., the last 11,500 years), through the preservation of sulfuric acid (H 2 SO 4 ) peaks.The synchronisation of these ice-core records on a global scale (Greenland and Antarctica), facilitates the generation of global volcanic radiative aerosol forcing estimates (Fig. 1a [5][6][7] ), which reveal periods of increased and reduced volcanic activity, for example, the Medieval Quiet Period 8 .The Medieval Quiet Period (700-1000 Common Era (CE)) was defined based on the scarcity of major tropical volcanic eruptions and minimal perturbations in solar forcing 8 .However, across this 300-year window, the Greenland ice-core records show an increase in background sulfate 9 (Fig. 1b).At the same time, proximal geological records of Icelandic volcanism, a source which has dominated the Greenland ice-core records through time 6,[10][11][12][13] , suggest an increase in eruption frequency between 750 and 1000 CE [14][15][16][17][18] (Fig. 1c).Despite this, few volcanic eruptions have been detected within the Greenland ice-core records across the period 19,20 (Table 1).
Existing volcanic detection methods, which facilitate the distinction between the non-volcanic (i.e., marine, biogenic, anthropogenic) and volcanic sulfate deposited within the ice, as sulfur has a number of different sources 21 , use the implicit assumption that volcanic sulfate depositional signals are sharp and short-lived (up to 3 years), due to the short atmospheric lifetime of sulfate 1,10 .However, geological records show that this assumption may not hold true, with evidence of eruptions lasting for several decades to centuries 22,23 , including eruptions from Iceland 6,18,24,25 .Due to the proximity of Greenland to Iceland, such long-lasting eruptions may be imprinted in the ice-core sulfate records, resulting in long-term changes which introduce ambiguity for secure detection and quantification of volcanic contributions.At present, clear diagnostic tracers unique to volcanic eruptions are missing within the polar ice which would facilitate in both their detection and source attribution and thus allow the distinction between both short-lived eruptions and longer lasting volcanic episodes.
Here, we present a set of four distinct diagnostic tracers (cryptotephra, sulfur isotopes, halogens, and heavy metals) within the ice-core records, which when combined, can characterize important source parameters i.e., exact location and plume height of past eruptions.We apply this multiproxy toolkit across the time period 700-1000 CE, resolving existing mismatches between proximal Icelandic records and global volcanic radiative aerosol forcing estimates through the attribution of these events within the Greenland ice cores.Finally, we explore the climate impacts at the onset of this period of elevated background sulfate, using proxies and documentary evidence.

Results
The seven Greenland ice cores (NEEM-2011-S1, NGRIP1, NGRIP2, EGRIP, RECAP, B19, and TUNU2013) used within this study are located across the Greenland ice sheet (see Fig. 2, Supplementary Table 1) between a b c NEEM-2011-S1 Greenland Fig. 1 | Global, Regional, and Proximal records of volcanic activity between 500 and 1500 CE. a Radiative forcing 5,7 between 500 and 1500 CE.A reduction in forcing can be observed between 700 and 1000 CE; b 21-year, 31-year, and 61-year median non-sea-salt Sulfur (nss-S) concentrations for NEEM-2011-S1.Concentrations are elevated compared to the periods before and after; c estimated eruptive frequency of Katla, Grímsvötn and Bárðarbunga-Veiðivötn volcanic systems in 250-year time windows, determined using the Catalog of Icelandic Volcanoes 154 and published literature 16,32,[155][156][157][158] .An increase in eruptive frequency from all volcanic systems is observed during the period 751-1000 CE.The Medieval Quiet Period (700-1000 CE) is indicated by the gray shading.
elevations of 2100 m above sea level (a.s.l) and 3000 m a.s.l.These ice cores are all dated by annual-layer counting (NEEM-2011-S1, NGRIP1, EGRIP, RECAP, B19, and TUNU2013 records are all volcanically synchronised to the NS1-2011 chronology 26 , whilst the NGRIP2 record uses the DRI_N-GRIP2 chronology 27 ).Together, they fully cover the period 700-1000 CE, with a dating uncertainty during this 300-year period better than ±2 years, constrained by two calendar dated cosmogenic age markers in 774 CE and 993 CE (see Methods, Supplementary Table 2).For consistency, we will use the mean of these two age models which may differ by ±1 year as will be shown later in individual figures.Details on the analytical methods employed for each ice core are provided (Methods, Supplementary Table 1).
We present stacked annual non-sea-salt (nss) sulfur concentrations (N = 6), which show an abrupt increase in 751 CE; prior to this, sulfur concentrations were lower overall, with less frequent variability (Fig. 3a).Median sulfur concentrations between 751 and 940 CE (the latter date marking the Eldgjá fissure eruption from the Katla volcano 28 ) were elevated by a factor of 1.6 compared to the period 700-750 CE (Fig. 3a).Superimposed on these increased background values are frequent sulfate peaks lasting between 1-12 years, with a sulfur concentration increase in 751-763 CE marking this transition.
Stacked nss chlorine (N = 6), fluorine (N = 2), and non-crustal (nc) bismuth and thallium (N = 4) concentration records strongly resemble the overall structure of the nss sulfur record, but with a more pronounced step-change at 751 CE (Fig. 3b, c).Throughout the 751-763 CE episode, element concentrations are elevated (8-fold increase for F -, 4-fold increase for Cl, 12-fold increase for Tl, and 5-fold increase for Bi) and remain elevated above pre-event levels for ~12 years.Though concentrations decrease again after 765 CE, the overall concentrations remain elevated (3-fold increase for F -and 2-fold increase for Cl, Tl, and Bi) and highly variable in contrast to the first half of the 8th century, where concentrations of nss chlorine (2.9 ng/g), fluorine (0.2 ng/g) and the heavy metals thallium (0.04 pg/g) and bismuth (0.06 pg/g) are low.Correlations between volatile metals (Tl, Bi) and volcanic volatiles (Cl, S) are strong throughout the period 751-940 CE, particularly for the 751-763 CE episode, suggesting a common source, transport, and deposition of these species.This is in agreement with ref. 29 who used Tl to assess past volcanic emissions of heavy metals (Bi, lead, cadmium) in Antarctic ice during the Common Era and ref. 30 who propose that volatile metals from volcanic sources can be transported to nearby icecore sites (i.e., 400 km distance) as water soluble aerosols, forming metalhalogen complexes.
To constrain the source of the volcanic events we searched for cryptotephra in the B19, TUNU2013, and RECAP ice cores.Guided by exactly co-registered analyses of size resolved insoluble particle concentrations and sulfur (see Fig. 3a, Supplementary Fig. 1), we successfully identified cryptotephra in the following years: 753 CE, 757 CE, 763 CE, 766 CE, 781 CE, 877 CE in B19 and TUNU2013.Notably, ice samples with cross sectional areas of up to 15 cm 2 and correlated to the years 763 CE and 877 CE, had high abundances of volcanic glass shards (up to 4000).Through geochemical analyses (see Methods), correlations were made with the Katla, Grímsvötn, Bárðarbunga-Veiðivötn, and Torfajökull Icelandic volcanic systems, as outlined below, with sample identification based on ice-core site, depth, and estimated age.A detailed description of shard characteristics and concentrations can be found in the Supplementary Data.
Major (see Fig. 4 and Supplementary Fig. 2) and trace element (see Supplementary Fig. 3) analysis of individual glass shards from six samples obtained from the three ice cores between 751 and 763 CE can be correlated to the Grímsvötn (B19_149.88m_753CE)and Katla (B19_149.60m_757CE;B19_148.88m_763CE;TUNU2013_158.20m_762CE)volcanic systems.The 753 CE glass shards were yellowish-green in color with thin-platy morphologies (see Supplementary Fig. 4a), whilst the 757 CE and 763 ± 1 CE glass shards (up to 50 μm in size within all ice samples) were yellowishbrown in color, with platy morphologies, and a TiO 2 -rich basaltic composition characteristic of the Katla volcanic system (see Supplementary Fig. 4b-f).We geochemically correlated the 763 CE ice-core tephra (B19_148.88m_763CEand TUNU2013_158.20m_762CE,RECAP_ 297.44m_762CE) with proximal material from the Icelandic Hrafnkatla eruption, including glass from identified tephra horizons in soil sections ~20 km from the main Katla vent (see Methods).The Hrafnkatla eruption 14,16,31 , also referred to as AT-8 or E2 within the literature 15,32 , was a major tephra-producing eruption from the Katla volcanic system, with an estimated total bulk volume of >4.5 km 3,15, .This is the earliest tephraproducing Katla volcanic eruption during the Common Era to be geochemically identified within the Greenland ice-core records and marks the end of a ~12-year period of continuous volcanic activity from Icelandic volcanic systems which here we term the Hrafnkatla episode.Limited geochemical analyses were obtained for the 766 ± 1 CE sample (B19_148.67m_766CE)and therefore, a definitive source could not be identified.However, due to the abundance of other material observed in the sample, including minerals and black particles of a similar size to the tephra shards, it is possible that the glass shards were deposited at the site with these other particles via the same secondary processes 33 .Additional geochemical analyses of this sample (on both the tephra glass and other particles) could help test this hypothesis.
Analysis of the 781 ± 2 CE sample (B19_147.29m_781CE)provided a correlation to a Grímsvötn tephra identified in Lake Lögurinn, Iceland 16 (see Fig. 4 and Supplementary Fig. 5).The B19_147.29m_781CE glass shards (up to 50 μm in size within all ice samples) were greenish in color, with platy morphologies (see Supplementary Fig. 4g).Within the lake, this horizon has been dated to 788 CE using the lake chronology (constructed using radiocarbon dating and historical tephra).This horizon is thus considered to be the most likely correlative for B19_147.29m_781CE as it is the only Grímsvötn event reported in the proximal records within ±50 years.
A correlation was also made between the 877 CE event in the B19 ice core to the Settlement Layer (877 ± 2 CE) in Iceland using proximal material which was analyzed simultaneously (see Fig. 4 and Supplementary Fig. 6).Both the silicic and basaltic components of the event, attributed to the Torfajökull and Bárðarbunga-Veiðivötn systems respectively, were identified in the B19 ice core, with glass shards having clear, greenish platy morphologies (see Supplementary Fig. 4h).The Settlement Layer is widely dispersed across Iceland and is considered an important horizon for separating the pre-historic (i.e., before 877 CE) and historic (i.e., after 877 CE) periods, defining the approximate date of the Viking settlement of Iceland 34 .The silicic and basaltic components of this horizon have previously been identified in the GRIP, GISP2, and TUNU2013 ice cores 19,35,36 , providing an important tie-point for the synchronisation of records.
The multiple correlations made within the Greenland ice-core records to Icelandic volcanic systems during the 751-940 CE period align well with reconstructions of volcanic activity from proximal Icelandic geological records.Proximal geological records show an increase in activity (i.e., eruption frequency and cumulative tephra volume) of the Katla and Bárðarbunga-Veiðivötn systems during the period 751-940 CE compared to earlier and later time-periods 14,15 (see Fig. 1c).Centennial mean tephra volume (minimum volume dense rock equivalent (DRE)) from Katla volcanic eruptions was 0.1 km 3 between 1-750 CE, increasing to 2.0 km 3 between 751-940 CE and decreasing to 0.1 km 3 between 941-1755 CE 15 .This fluctuation has been linked by volcanologists to the presence of a relatively shallow magma chamber that was active prior to the 939 CE Eldgjá eruption 15 .Owing to the strong agreement between the two records which suggest Icelandic volcanic systems were highly active and productive over  68 .The red rectangle highlights the location of the correlated volcanic systems within this study, with a detailed map found in Supplementary Fig. 10.
the period 751-940 CE, we subsequently term this period the Icelandic Active Period (IAP).
To further investigate this abrupt shift in background sulfate, and thus the onset of the IAP which commenced with the Hrafnkatla episode (751-763 CE), we measured S isotopes in sulfates at high-temporal resolution across this 12-year period on the TUNU2013 ice core (see full details in Methods).Mean δ 34 S values (see Methods for details) across the Hrafnkatla episode are +6.4‰(see Fig. 5b).These values are lower compared to previously measured mean δ 34 S TUNU2013 background samples (+9.5‰ ± 0.93‰) 37 , and thus suggest a largely volcanic origin.Background sulfate consists of both marine and volcanic emissions, and an increase in volcanic sulfate results in lower δ 34 S values within the ice (i.e., values closer to 0‰), as observed across the Hrafnkatla episode.Background corrected δ 34 S values over this period are +2.6‰(±2.5‰) (see Supplementary Fig. 7).The anti-correlation between the local S peaks and the δ 34 S further suggests large contributions of sulfur with a volcanic origin.Cryptotephra was identified in samples with low (at or below −0.6‰) δ 34 S values (see Fig. 5a).
Previous polar ice-core studies [37][38][39][40][41][42] have successfully utilised Δ³³S (see Methods for details) isotopic values to determine whether sulfate aerosol formation occurred in the stratosphere, or more precisely above the ozone layer, which above Iceland is located between 11 and 15 km 43 , thus providing an indication of relative plume height.Our analysis of Δ 33 S across the Hrafnkatla episode has revealed that this volcanic sulfate mainly has values varying between -0.2 and +0.2‰ for ~12 years (see Fig. 5c).These values are indistinguishable from the 2. s.d. of a non-mass-independentfractionation secondary standard (see Methods) and indicate a ~0‰ Δ 33 S value, and hence a dominant tropospheric origin for the S fallout.However, in 762/3 CE, values fall to -0.5‰, suggesting that this prolonged elevated tropospheric phase culminated in a short-lived stratospheric phase, which through cryptotephra geochemical analysis we have correlated to the Katla Hrafnkatla volcanic eruption (see Fig. 5a).This is the first detection within the ice-core records, to the best of our knowledge, of anomalous Δ 33 S values for an Icelandic eruption, and unlike tropical stratospheric eruptions (i.e., Samalas 1257 CE and Tambora 1815 CE 37 ), the Δ 33 S anomaly is much smaller for this high-latitude stratospheric Icelandic eruption.This difference can be attributed to the simultaneous deposition of tropospheric and stratospheric S, resulting in a muted anomalous Δ 33 S value 44 .
No sulfur isotope analyses have been performed for the remaining part of the IAP, but sulfur concentrations in the Greenland ice cores frequently exceed typical levels which were prevalent for the first half of the 20th century during periods of intensive industrial activity and related sulfur emissions across the North Atlantic region 45 (Fig. 6a).Considering the high activity of these Icelandic volcanic systems across the study period, as revealed through the application of the above proxies, traditional approaches to isolate and quantify volcanic eruptions and their sulfur emissions are inevitably biased.This has resulted in an underestimation in the number, total duration, and SO 2 emission strength of volcanic events.For example, both the reconstructions of ref. 46, the recommended volcanic forcing for the Paleoclimate modeling intercomparison project (PMIP) 3 model experiments, and ref. 47, did not detect any volcanic events between 751 and 763 CE (see Supplementary Fig. 8).This is because existing reconstructions of volcanic sulfate deposition used in these forcing estimates (e.g., NGRIP1 47 and NEEM-2011-S1 9 ) across the 700-1000 CE period have attributed most of the background variation to non-volcanic sources.Annual mean volcanic sulfate deposition across this period was only 1.0 kg km −2 yr −1 of sulfate and 2.1 kg km −2 yr −1 of sulfate for NGRIP1 and NEEM-2011-S1 respectively.We therefore revised existing reconstructions of volcanic sulfate deposition, now attributing long-term background variations of sulfur as volcanic in origin (see Methods).As a result, we were able to reconstruct overall higher volcanic sulfate deposition rates for NGRIP1 (2.4 kg km −2 yr −1 of sulfate) and NEEM-2011-S1 (4.1 kg km −2 yr −1 of sulfate) (Fig. 6b).Additionally, we reconstruct new volcanic sulfate depositional records from NGRIP2, B19, TUNU2013, and EGRIP (see Fig. 6c), thus doubling the number of ice-core sulfate records from Greenland currently used in reconstructions of volcanic forcing (i.e., by PMIP4) during the Medieval period 3,5 .

Discussion
Previous suggestions that the period 700-1000 CE was a time of quiescence are not supported by the proximal Icelandic geological records or the Greenland ice-core records.Our application of a multi-proxy geochemical toolkit has revealed heightened volcanic activity from several Icelandic volcanic systems and includes both long-lasting, multi-centennial episodes such as the Icelandic Active Period (IAP) (751-940 CE) and short-lived episodes such as the Hrafnkatla episode (751-763 CE).Both are characterized by prolonged and elevated enrichments of volcanic metals (Tl, Bi) and volatiles (F, Cl, S).
Within Antarctic ice cores, a pronounced elevated halogen and metal enrichment lasting 192 years (~17.7 ka) was attributed to regional West Antarctic volcanism 23,30,48 .Mason et al. 30 suggested that the long-lasting excessive deposition of these volatile species is indicative of a proximal eruption, with magma source chemistry, eruptive processes, and plume transport also contributing to the signature within the ice.Whilst a great deal is still unknown about the behavior of many volatiles during transport, prolonged stratospheric transport (i.e., from a distal eruption) is likely to result in greater loss due to scavenging, and thus reduced or negligible deposition on polar ice sheets.Modern plume and ground-based measurements of the Hawaiian Kīlauea 2018 volcanic eruption support this, with ref. 49 observing a decrease in atmospheric burdens of HF, HCl, and volatile metals with distance from the source due to their solubility.Owing to this decay, and the proximity of Iceland to our Greenland ice-core sites (~1000-1500 km), we suggest that in the Greenland ice cores this metal signature could be used as an indicator for Icelandic eruptions in the absence of cryptotephra.These trace metals are now routinely analyzed over the entire length of deep ice cores and thus are readily available for improved source attribution of volcanic signatures within ice cores 23,27 .
Across the IAP, we propose the prolonged elevated concentrations of volcanic volatiles and volatile metals that are observed can best be explained by continuous volcanic activity dominated by several Icelandic sources (Katla, Bárðarbunga-Veiðivötn, and Grímsvötn).This is supported by the general increased eruption frequency from these systems, with Katla producing 10 times more tephra volume between 750 and 950 CE compared to time-periods before and after 15 .Based on proximal geological and regional glaciochemical records, we postulate that eruptions across the IAP likely encompassed a wide range of volcanic eruption types, involving explosive stratospheric eruptions as well as prolonged fissure and flood basalt eruptions and quiescent degassing.This includes, for example, the long-lasting Hallmundarhraun lava flow, which occurred between 880 and 910 CE, and is considered the third largest effusive eruption since the settlement of Iceland (~877 CE) with an estimated erupted volume of up to 4 km 3[ 50 .Overall, as supported by S isotopes, a large fraction of these aerosol emissions remained predominantly in the troposphere and lowermost stratosphere.This strong and long-lasting tropospheric S contribution led to the estimate of volcanic sulfate deposition in Greenland across this 300-year period being revised and increased by a factor of 2 (see Fig. 6b, c), thereby challenging the notion of this being a Medieval Quiet Period.Previous underestimations in volcanic sulfate may largely stem from assumptions and conservative approaches used in the detection of volcanic events within the ice-core records.Prior calculations of volcanic sulfate deposition rates across this 300-year period assumed that the background levels of sulfate in the ice were not volcanic in origin or time variable 9,46,47 .This is inherently linked to the tailoring of detection methods 2,10,51 towards the isolation of short-lived stratospheric perturbing volcanic signals, as these eruptions have typically been associated with stronger climatic impacts on account of the longer residence time of the aerosols within the atmosphere 1,2 .Consequently, the 939 CE Icelandic Eldgjá eruption 28 was identified as the only large eruption across this period, which when accompanied by the relative stability of solar activity across the period 700-1000 CE, subsequently resulted in the period being considered one of quiescence.Therefore, our revised estimates across this 300-year period, with the inclusion of a long-lasting tropospheric S contribution, have important implications for climate model simulations, such as those from PMIP3 and PMIP4, which use the year 850 CE as the start of the past1000 simulations covering the pre-industrial millennium (850-1849 CE) 3,4 .These models use only stratospheric volcanic forcing where the tropospheric aerosols are not included, in contrast to historic climate simulations (starting in 1850 CE), which also include tropospheric aerosol forcing.
The role of tropospheric volcanic aerosols in Earth's system should not be neglected.These aerosols not only have direct impacts (i.e., reflection of incoming radiation and surface cooling 52 ), but can act as ice-nucleating particles and cloud condensation nuclei, thereby influencing cloud properties and the hydrological cycle 53 .Despite the importance of tropospheric aerosols, little is known about pre-industrial concentrations, a time when the atmosphere showed greater sensitivity towards tropospheric volcanic aerosol inputs (e.g., from quiescent degassing) due to the absence of anthropogenic aerosols 54,55 .Developing a baseline of pre-industrial sulfuric aerosol concentrations would therefore help constrain uncertainties in radiative forcing estimates which are used in climate model simulations.
These ongoing uncertainties and different aerosol baseline conditions during the pre-industrial era have presented challenges when assessing the climate impacts of prolonged volcanic episodes.Here, we further investigate the Hrafnkatla episode, as this marks the onset of the IAP and this abrupt shift in background sulfate, by comparing forcing reconstructions to palaeoenvironmental records from the North Atlantic region where they are abundant and well-dated.
Due to atmospheric circulation patterns, which confine aerosols from high-latitude volcanic eruptions to their respective hemispheres, the climatic impacts of large Icelandic eruptions are more spatially limited than those of tropical eruptions of a similar magnitude.Here, we consider climatic impacts to be those which perturb the system on a seasonal to multidecadal time scale.The climatic impacts associated with the Hrafnkatla event (the stratospheric eruption which marks the end of Hrafnkatla episode in 762 CE) are generally comparable to other large Icelandic volcanic eruptions such as Eldgjá 939 CE and Laki 1783 CE, with strong winter cooling anomalies across Europe 28,56,57 and perturbations to the hydrological cycle 56,[58][59][60] .At the onset of the Hrafnkatla episode in 751 CE, a cooling signature is recorded in two speleothem records from central Europe (Milandre cave, Switzerland 61 and Spannagel cave, Austria 62,63 ), and is most pronounced (up to 1 °C) at around 762-763 CE (see Fig. 7d).These records have a high-resolution (~2 years per sample) and age error of around ±16 years (1.s.d) across the study period.Although several processes may govern the isotopic values in speleothems (i.e., climatic and post-depositional), these archives have been shown to reliably capture past climate variability with a demonstrated sensitivity to the post-climatic effects of volcanic eruptions 61 .
Historical records cite the winters of the same and subsequent years of this pronounced cooling anomaly (762-764 CE) among the harshest of the Common Era, with widespread frost and ice, crop failures, and famines documented from Ireland to the Black Sea [64][65][66] (see Fig. 2).Indeed, the firsthand eyewitness account of Theophanes the Confessor documents an extremely severe winter in the area of the Black Sea and the Bosporus from October 763 CE to February 764 CE and describes that "in early October [763 CE], there was very bitter cold, not only in our land, but even more so to the east, the north, and the west, so that on the north coast of the Pontos [the Black Sea] to a distance of 100 miles the sea froze […] In the month of February [764 CE] this ice was […] split up into many different mountainlike sections.[…] One of the icebergs struck the jetty of the Acropolis [of Constantinople] and crushed it" [67][68][69] .These harsh conditions documented across Europe coincide with a period of extreme annual low stand values of the river Nile (20 of the 52 lowest 5 th percentiles occurred between 750 to 790 CE) 70 (Fig. 7c).The flow of the river Nile is highly sensitive to shifts in the Intertropical Convergence Zone (ITCZ), with the headwaters of the Nile located in the equatorial regions.The ITCZ is displaced southward in response to the high atmospheric burden of sulfate aerosols and subsequent cooling within the Northern Hemisphere atmosphere following highlatitude eruptions, resulting in varied precipitation patterns in the equatorial regions and reduced stream flow of the Nile 56,71,72 .The low flow following these short-lived high-latitude eruptions is of short duration (i.e., the year after the eruption), and therefore, the clustering of extreme annual low stand values occurring during and following the Hrafnkatla episode may suggest a different, more persistent tropical hydroclimate response following prolonged volcanic emissions as was the case during the 20th century which was dominated by anthropogenic sulfur emissions [73][74][75][76][77] .Climate models may provide further insight, as they can better explore the mechanisms and feedbacks following prolonged volcanic emissions.
In contrast to many historic large eruptions, no major summer cooling anomaly is recorded in Northern Hemisphere tree-ring reconstructions in the year following the explosive stratospheric end-phase of the Hrafnkatla episode (i.e., 764 CE) [78][79][80] .The cause for this absence of a summer cooling anomaly is unknown; however, given that the explosive Hrafnkatla event followed a 12-year period of sustained aerosol emissions, increased background aerosol concentrations may have modified the climatic response.Arctic sea-ice plays an important role in amplifying and sustaining climatic responses to large explosive volcanic eruptions, and typically, following an explosive volcanic eruption there is an expansion in sea-ice [81][82][83] .In contrast, methanesulphonic acid (MSA) concentrations (a proxy for sea-ice) from the TUNU2013 Greenland ice core 84 (Fig. 7b) show no clear indication of a strong sea-ice expansion around the 760s CE.They do however, indicate a general reduction in Arctic sea-ice across the IAP, coincidental with an increase in surface air temperature, peaking in the late 10th century (Fig. 7a).This transition into a North Atlantic Medieval Warm Period (MWP) coincided with the Viking expansion across the Arctic region to settle first in Iceland (~877 CE) and then in southern Greenland (~985 CE) 85,86 .Arguably, these warmer temperatures and ice-free conditions may have facilitated such an expansion.Ultimately, this is an area that requires further research, with continued uncertainties regarding the climatic response and feedbacks associated with prolonged volcanic episodes, particularly during the preindustrial era.
In conclusion, our application of a multi-proxy geochemical toolkit has revealed that the period 751-940 CE was a time of heightened Icelandic volcanic activity, which commenced with the Hrafnkatla episode in 751 CE.This is in contrast to previous suggestions of the period 700-1000 CE being a time of quiescence, as indicated by global radiative aerosol forcing estimates.This period was characterized by heavy metal and halogen enrichments, with these species most pronounced at the time of local Icelandic eruptions.We have therefore postulated that both heavy metals and halogens can be used as distinct tracers of Icelandic volcanic eruptions within the Greenland ice core records and aid in discerning the source in the absence of cryptotephra.Through our attribution of several volcanic signals within the Greenland ice-core records to Icelandic volcanic systems, we have been able to resolve existing mismatches between these regional records and proximal geological records.Our revision of volcanic sulfur deposition across this period identified a greater tropospheric contribution of sulfate aerosols, which is currently not considered in existing volcanic forcing estimates.We found that the explosive stratospheric end-phase of the Hrafnkatla episode coincided with a strong winter cooling anomaly throughout Europe and could be considered as a likely cause, owing to the absence of major tropical eruptions during this time period.However, unlike other large explosive volcanic eruptions, no major summer temperature anomaly has been observed thus far.Finally, we suggest that more research into prolonged volcanic episodes, such as the IAP and Hrafnkatla episode, is required due to the uncertainties in the associated climate feedbacks to sustained tropospheric aerosol emissions, especially during the pre-industrial period.Ultimately, resolving these uncertainties will be of particular importance when investigating other time intervals such as the last deglaciation to early Holocene, where evidence suggests long-lasting volcanic episodes occurred frequently, with some lasting up to 70 years 6,24 .

Methods
Ice-core analysis Information regarding the seven Greenland ice-core records used within this study NEEM-2011-S1, NGRIP1, NGRIP2, EGRIP, RECAP, B19, and TUNU2013) can be found in Supplementary Table 1.Elemental data (including sulfur, chlorine, bismuth, and thallium) for NGRIP2, NEEM-2011-S1, and TUNU2013 9,26,27 ice cores were utilised alongside new measurements of these elements in the B19 ice core (78.00°N36.40°E,150.4 m depth, 2234 m a.s.l.) 87 .All data from these four ice cores were obtained through continuous flow analysis (CFA) coupled to two inductively coupled plasma mass spectrometers (ICP-MS, Element2, Thermo Scientific) at the Desert Research Institute, USA (DRI), providing a simultaneous determination of a broad spectrum of elements at a sub-annual resolution 23 .Additional sulfate (SO 4 2− ), chloride (Cl -), and fluoride (F -) measurements for EGRIP and NGRIP1 were undertaken at the Alfred Wegner Institute, Germany (AWI), and Physics of Ice, Climate and Earth, Denmark (PICE), respectively, using ion chromatography.NGRIP1 had been sampled discretely at 5 cm depth resolution; EGRIP samples were taken during the CFA analyses at the University of Bern using a fraction collector with a nominal depth resolution of ~6 cm 88 .For the determination of anionic species in EGRIP, a two-channel Dionex ICS 5000+ reagent-free HPIC system (Thermo Fisher) with a Dionex IonPac AS18-Fast-4 μm analytical column was used.Details of the entire sampling procedure and analytical conditions are described in detail in ref. 89.Between 700 and 1000 CE the time resolution is about 0.5 years in EGRIP.For quantification of ionic species in NGRIP1 a two-channel Dionex 500 microbore ion chromatograph was used 47,90,91 resulting in a time resolution of about four samples per year during the first millennium CE.In addition, methanesulphonic acid (MSA) was measured at DRI continuously over the TUNU2013 ice core using an electrospray triple-quad mass spectrometer (ESI-MS-MS) 84 alongside stable water isotope ratios (δ 18 O), determined using a cavity ring-down water isotope analyzer (L2130-i, Picarro Inc) 92 .All parameters measured for each ice core and used within this study can be found in Supplementary Table 1.

Ice-core data processing
We estimated non-sea-salt sulfur (nss-S) concentrations by correcting total S concentrations for the sulfur contribution from sea salt using coregistered Na concentrations as a sea salt tracer and assuming a sulfur to sodium mass ratio in bulk sea water of 0.084 93 .Similarly, we estimated non-sea-salt chlorine (nss-Cl) concentrations by correcting total Cl concentrations for the chlorine contribution from sea salt using coregistered Na concentrations as a sea salt tracer and assuming a chlorine to sodium mass ratio in bulk sea water of 1.796 94 .To calculate the noncrustal fractions (nc-Tl and nc-Bi) and enrichment factors (Tl-EFCe and Bi-EFCe, see Supplementary Fig. 9), we used cerium (Ce) as a conservative tracer for mineral dust and the "mean continental crust" Tl/Ce ratio of 0.01172 and Bi/Ce of 0.00198, respectively 95 .Ce concentrations were scaled by a factor of 1.7 to account for an under-recovery of Ce (of 60%) during the continuous-flow analysis at the Desert Research Institute, whereas recovery was ~100% for both Tl and Bi 96 .
Targeted sampling CFA measurements of insoluble particles at two size fractions (2.6-4.5 μm and 4.5-9.5 μm) were used alongside peaks in sulfur to assist with isolating events and sampling for tephra.Particle measurements were undertaken using an inline Abakus® laser particle counter 97 during the CFA analyses at DRI.Although the particle concentrations alone cannot be considered as a direct proxy for tephra, the application of this approach, based on the exact stratigraphic co-registration of insoluble particulates and volcanic volatiles in the polar ice cores, has proven highly successful for detecting cryptotephra layers 19,33,[98][99][100][101][102][103] .One 11 cmlong sample was obtained for tephra analysis from the RECAP ice core between 297.44-297.55 m for which a non-targeted continuous sampling approach was followed as described by ref. 104.

Ice-core chronologies
The NS1-2011 26 and DRI_NGRIP2 27 chronologies are used within this study.Both chronologies were constructed through annual layer counting of seasonally varying impurity tracers constrained by well-dated stratigraphic markers such as historic volcanic eruptions, volcanic dust veils, and solar energetic particle events (see Supplementary Tables 2, 3).Dating uncertainty in both chronologies during the period 700-1000 CE is ±2 years.This is supported by the occurrence of short-lived cosmogenic nuclide events in 774 CE and 993 CE, which are observed in both the ice core and tree-ring records 26,105,106 and by exact matches of dendrochronologically dated volcanic eruptions and corresponding ice-core tephra identifications in 822/ 823 CE 33,105 and 946/947 CE 28 .Both chronologies are within errors consistent with the GICC21 chronology 107 , but are more precise during the period 700-1000 CE because of the additional age constraints.To transfer the age of the RECAP ice core tephra sample (available on the GICC05 chronology 108 ) to the NS1-2011 chronology, we used an age transfer of +7 years following ref.26.

Cryptotephra analysis
Targeted tephra samples were obtained for 11 volcanic events from TUNU2013, B19, and RECAP ice cores and analyzed at Queen's University Belfast, UK (QUB); University of Bern, Switzerland (UniBe); and the University of Edinburgh, UK (Edin), respectively (see Supplementary Data).The targeted sampling approach allowed us to sample comparably large cross sections of archived ice ranging between 16.5 cm 2 (RECAP) and 15 cm 2 (B19) thus increasing the chance of retrieving greater abundances of tephra shards.In addition, samples that were analyzed for S isotopes were also centrifuged prior to analysis, facilitating the screening for cryptotephra.The method for the preparation for tephra identification and geochemical analysis of samples analyzed at QUB and Edin can be found in ref. 98 and ref. 104 respectively, and secondary standards for analyses obtained in this study can be found in the Supplementary Data.The preparation of samples for tephra identification and geochemical analysis at UniBe followed a method adapted from several recent papers [109][110][111][112] and the full methodology used within this study can be found in Supplementary Methods A. Samples where glass tephra shards were identified were analyzed using Electron probe microanalysis (EPMA) at the Institute of Geological Sciences, University of Bern, Switzerland.The samples were carbon coated prior to major elemental analysis using a JEOL JXA-8200 probe, with operating conditions of an accelerating voltage of 15 kV, beam current of 2 nA, and beam diameter of 5 μm.Secondary standards of Cannetto Lami Lava, Lipari and BCR-2Ga were measured during the analytical session to ensure stability and allow session comparisons.Full details about the operating conditions and secondary standards at UniBe can be found in the Supplementary Data.
Following the filtering of all obtained data to totals >90 wt% and the removal of Cl and P 2 O 5 (to facilitate correlation with other data sets), the data were normalized to an anhydrous basis following standardised practices 113 .The data were further inspected to remove any outliers, both those occurring from analytical issues (e.g., analysis of the resin) and those which appeared outside of the main population.All raw data can be found in the Supplementary Data.Data were then correlated to published tephra geochemical data, with potential correlations explored graphically.In order to make robust correlations, samples were analyzed alongside proximal material when available, which was obtained in summer 2021/ 2022 fieldwork.The full preparation method for proximal material and sampling strategy can be found in Supplementary Methods B.
Trace element analysis was conducted on individual glass shards from four ice core and three proximal samples covering the period 751-763 CE using laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) at St Andrews Isotope Geochemistry lab, UK (StAiG).Laser ablation used an Applied Spectra RESOlution system, featuring a Laurin Technic S155 2-volume cell under a pure He atmosphere (370 mL min -1 ) with supplemental N 2 (5 mL min −1 ) to improve sensitivity in the ICP-MS.Analyses were undertaken using a 20 μm beam diameter and a laser energy of 3.5 J cm −2 with a 4 Hz repetition rate, 20-s background measurement, and 30-s ablation time.The laser cell gases and ablated sample were mixed with Ar (1 L min −1 ) before being homogenized in a signal smoothing device prior to introduction to an Agilent 8900-QQQ for mass spectrometry.Measurements of glass standards of Old Crow, Sheep Track, and ML3B, and synthetic glass standards of NIST610 and T1G which were used to correct for drift, bracketing measurements, can be found in the Supplementary Data.Samples were normalized using Al 27 , which was determined from the same samples during previous EPMA at the UniBe.A detailed methodology can be found in ref. 114.All raw data can be found in the Supplementary Data.

Sulfur isotope analysis
The ratio of isotopic 34 S, 32 S, can be used to determine the origin of the sulfate (i.e., marine or volcanic) [115][116][117] .This distinction can be made through the measurement of δ 34 S (defined as, δ x S ¼ ð x S= 32 SÞ sample =ð x S= 32 where x is 33 or 34).Marine sulfate typically has values between 16-20‰, whilst volcanic isotopic values are closer to zero [117][118][119][120][121] .By measuring the minor isotope, 33 S, one can assess Mass-independent Fractionation (MIF) by calculating Δ³³S (defined as, Δ³³S isotopic values can be used to determine whether sulfate aerosol formation occurred in the stratosphere (i.e., above the ozone layer) since UV photochemical reactions generate MIF and non-zero Δ³³S 40,[122][123][124] .As a result, sulfur deposited over the ice sheet that has reached the stratosphere has Δ³³S ≠ 0. Due to the ozone layer, which filters out radiation of the specific wavelength that causes MIF, the sulfate which remains in the troposphere/lowermost stratosphere is not exposed to UV radiation.Therefore, Mass-Dependent Fractionation (MDF) is generated and tropospheric/ lowermost stratospheric sulfate has a Δ³³S isotopic value that is equal to zero.Through these different fractionation processes, it is possible to distinguish stratospheric sulfate from tropospheric/lowermost stratospheric sulfate and provide an insight into the explosivity and plausible plume heights of past eruptions.
We undertook S isotopic measurements on samples from TUNU2013 across an elevated period of sulfate concentrations.TUNU2013 was continuously sampled at a 3 cm depth resolution (cross sections of 7 cm 2 ) between 740 and 765 CE for sulfur isotope analysis at the University of St Andrews Isotope Geochemistry Laboratory (StAiG), UK.The detailed methodology can be found in the methodology of ref. 37. Briefly, the sulfate concentration of each sample was first determined by ion chromatography using a Metrohm 930 Compact IC Flex.Sulfate was then extracted from the sample using column chemistry and triple sulfur isotopes ( 32 S, 33 S, and 34 S) were subsequently measured using a Neptune Plus Multi-Collector-ICPMS 125 .These measurements were conducted alongside an established secondary standard (Switzer Falls 126 ) and procedural blanks.The procedural blank values (0.18 ± 0.08 nmol S, δ 34 S = 6.8 ± 4.3‰, 1 s.d., n = 11) were used to correct all measured ratios.The δ 33 S and δ 34 S values were used to calculate Δ 33 S. The external reproducibility of the secondary standards run alongside these samples (0.19‰, 0.08‰, and 0.09‰ for δ 34 S, δ 33 S, and Δ 33 S respectively, 1 s.d, n = 9) is used as the reported uncertainty on isotope ratios.For some samples there was only enough sulfate for a single isotope measurement, and thus the uncertainty on these samples is slightly elevated and was calculated from the external reproducibility of the secondary standard when only one replicate was considered (0.22‰, 0.1‰, and 0.1‰ for δ 34 S, δ 33 S and Δ 33 S respectively, 1 s.d, n = 9).The isotope data can be found in the Supplementary Data.Following the approach of ref. 38.Gautier et al. 127 , we used isotope mass balance to solve for the sulfur isotopic composition of volcanic sulfate (see Supplementary Fig. 7) using the following equation: where δ x S meas is the measured isotopic value, f background is the fraction of background sulfate (f background = sulfate concentration of background/total sulfate), f volc is the fraction of volcanic sulfate (f volc = 1-f background ), and δ x S background is the isotopic composition of the background sulfate.We use the mean and standard deviation of sulfate concentration and sulfur isotope composition from 13 TUNU2013 non-volcanic background samples analyzed throughout the Common Era to calculate f background and δ x S background .Uncertainties on this calculation are propagated using Monte Carlo simulations, and only samples with more than 57% volcanic sulfate are plotted, similar to that suggested by ref. 127.

Volcano detection
Several statistical methods have been applied in previous studies to the icecore sulfate records to isolate volcanic events from the background signal 128,129 .These methods often use running medians of multi-decadal lengths for the separation of time-variant volcanic and background sulfate (a combination of time invariant marine and volcanic contributions).Their application has proven successful in the detection of climate impacting (medium and large sized) volcanic events, but ultimately has resulted in an underestimation in the total amount of volcanic sulfate.Within this study, we therefore adjusted the previous methods 9,128,129 , across the NEEM-2011-S1, NGRIP1, NGRIP2, EGRIP, B19, and TUNU2013 ice core records, as follows 1 : a 91-year Running Median (RM) was used on the annually averaged nss-S record on periods unaffected by strong changes in volcanic background emissions (i.e., 300-700 CE for all cores but B19, for which 950-1250 CE was used) to estimate the natural background sulfate levels 2 ; a Median of Absolute Deviation (MAD), calculated from the RM, was used in peak detection.As with the calculation of the RM, the MAD for the period before (i.e., 300-700 CE and 950-1250 CE for B19) was used.A peak was considered volcanic if it passed the upper threshold (K = 3, estimated as RM plus 3* MAD).The duration of the event was determined when it passed the lower threshold (K = 1, estimated as the RM plus 1* MAD).These upper and lower thresholds were selected by validation on well-known historic eruptions 3 ; volcanic peaks were then removed to calculate the non-volcanic background (S RRMi).To further calculate the amount of sulfate deposited, S RRMi (from 300-700 CE and 950-1250 CE for B19) was subtracted from the average annual nss-S and then multiplied by the accumulation rate of the drill site 4 ; finally, cumulative volcanic sulfate deposition was calculated for each event by summing the sulfate deposited across the total duration of the event.
Temperature reconstruction for Milandre cave (Switzerland) 600-1200 CE The Milandre cave (see Fig. 2) temperature reconstruction is based on oxygen isotope measurements (δ 18 O) of speleothem calcite and the data can be found in the Supplementary Data.The statistically significant correlation (r 2 = 0.56, p = 0.0001) between stalagmite calcite δ 18 O 61 and historical and instrumental temperature records from western Switzerland 130 allows the development of a transfer function to calculate mean annual air temperatures (MAAT) for the last two millennia for western Switzerland.
Temperature reconstructions from Greenland δ 18 O stack 600-1120 CE As a proxy for annual mean surface air temperatures (SAT) over Greenland, we collated published stable oxygen isotope (δ 18 O) records across Greenland fulfilling the three selection criteria: continuous coverage between 600 and 1120 CE, with annual time resolution, and dating precision and accuracy better than ±3 years.A total of six records met these requirements, from the Dye-3, NGRIP, NEEM, GRIP, EGRIP, and TUNU2013 ice cores 107,131,132 .
We then compiled our dataset of these individual δ 18 O anomaly records (relative to a 1901-1960 reference interval) into a single stack by calculating the arithmetic mean δ 18 O anomaly for each year.We further apply a 31-year running mean filter to our time series and estimated temperature anomalies, by applying the spatial calibration slope for Greenland of 1/0.67 °C per ‰ 133,134 .
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.© The Author(s) 2024

Fig. 2 |
Fig. 2 | Cooling anomalies in the North Atlantic region.Map showing the location of the Greenland ice core sites (B19, TUNU2013, NEEM-2011-S1, NGRIP, EGRIP, and RECAP) used within this study.The locations of the Milandre (MIL) and Spannagel (SPA) speleothem records (red crosses) and River Nile low stand record (blue cross) utilised within this study are indicated.The black diamonds denote harsh winters in the years following the 763 ± 1 CE Hrafnkatla volcanic eruption across Europe, as cited in the historical written records68 .The red rectangle highlights the location of the correlated volcanic systems within this study, with a detailed map found in Supplementary Fig.10.

Fig. 3 |
Fig. 3 | Greenland ice core volcanic tracers across 700-1000 CE. a Non-sea-salt Sulfur (nss-S) stack (N = 6) of Greenland ice cores and associated insoluble particle record (3-10 μm).The locations of existing sampled tephra (gray) across this time period are denoted, where correlations have been made to the 822 CE eruption of Katla, Iceland 33 , 853 CE eruption of Mt Churchill, USA 20,100 , and 946 CE Changbaishan, China/DPRK 20,103 .New samples which form part of this study are shown by the red circles; b Non-sea-salt chlorine (Nss-Cl) stack (N = 6) and Fluoride (F -) stack (N = 2) of Greenland ice cores; c Non-crustal thallium (Nc-Tl), non-crustal bismuth (Nc-Bi) stack (N = 4) of Greenland ice cores.The blue shading denotes the Hrafnkatla episode (751-763 CE), and the beige shading shows the duration of the Icelandic Active Period (751-940 CE).

Fig. 5 |
Fig.5| Hrafnkatla episode.a Non-sea-salt Sulfur (nss-S) from TUNU2013 ice core and associated insoluble particle record (5-10 μm) across the period 742-768 CE.Targeted tephra sample (gray shading) and sulfur isotope (green shading) locations (which were screened for tephra) are shown, with the geochemically confirmed occurrences of cryptotephra indicated by red circles; b Measured δ 34 S (‰) and sulfur concentrations versus time recorded in the TUNU2013 ice core; c Measured Δ 33 S (‰) and sulfur concentrations versus time recorded in the TUNU2013 ice core.Blue shading indicates stratospheric Δ 33 S range (0.0 ± 0.2‰); d Non-crustal thallium (Nc-Tl), non-crustal bismuth (Nc-Bi) and non-sea-salt chlorine (Nss-Cl) record for TUNU2013 ice core.A pronounced enrichment can be observed at the beginning of the Hrafnkatla episode.

Fig. 6 |
Fig. 6 | Volcanic sulfate deposition in Greenland ice cores.a Nss-S (stack N = 6) across the period 700-900 CE and 1800-2000 CE (stack N = 4).Mean nss-S for the periods 1901-1950 (red dashed line) and 1951-1996 (brown dashed line) are shown for comparison; b Annual volcanic sulfate deposition from NGRIP147 (purple line) and NEEM-2011-S1 9 (red line) across the period 700-1000 CE; c Annual volcanic sulfate deposition (this study) using a stack (N = 6) of Greenland ice cores, where calculated cumulative volcanic sulfate for the period 750-1000 CE is greater (658 kg km −2 compared to 305 kg km −2 and 608 kg km -2 for NGRIP1 and NEEM-S1 respectively).Tephra occurrences with confidently attributed sources throughout this period are indicated by the red circles.A comprehensive overview can be found in ref.20.

Fig. 7 |
Fig. 7 | Palaeoenvironmental and proxy records 600-1120 CE. a Surface air temperature (SAT [°C]) 131 inferred from a δ 18 O stack (N = 6) of Greenland ice cores with 31-year mean (yellow line).The timing of the Hrafnkatla episode (blue) is shown, alongside the onset of the Medieval Warm Period (red) and the settlement of Iceland and Greenland (green dashed line); b TUNU2013 MSA 84 with 31-year mean (yellow line).Outliers have been detected and removed with a threshold of 1.5 times the Interquartile range; c Nile annual low stand 70 with 31-year mean (dark blue line), with the lowest 5th percentiles marked by brown circles; d Spannagel cave 62,63 (gray line) and Milandre cave (blue line)[this study] speleothem temperature anomalies [°C] with reference to mean annual air temperature (MAAT) across the 600-1120 CE period are shown.Icelandic (red triangles) and NHET (gray triangles) volcanic eruptions during the Icelandic Active Period (yellow) are shown.