Volcanic glass properties from 1459 C.E. volcanic event in South Pole ice core dismiss Kuwae caldera as a potential source.

A large volcanic sulfate increase observed in ice core records around 1450 C.E. has been attributed in previous studies to a volcanic eruption from the submarine Kuwae caldera in Vanuatu. Both EPMA-WDS (electron microprobe analysis using a wavelength dispersive spectrometer) and SEM-EDS (scanning electron microscopy analysis using an energy dispersive spectrometer) analyses of five microscopic volcanic ash (cryptotephra) particles extracted from the ice interval associated with a rise in sulfate ca. 1458 C.E. in the South Pole ice core (SPICEcore) indicate that the tephra deposits are chemically distinct from those erupted from the Kuwae caldera. Recognizing that the sulfate peak is not associated with the Kuwae volcano, and likely not a large stratospheric tropical eruption, requires revision of the stratospheric sulfate injection mass that is used for parameterization of paleoclimate models. Future work is needed to confirm that a volcanic eruption from Mt. Reclus is one of the possible sources of the 1458 C.E. sulfate anomaly in Antarctic ice cores.

Introduction

The second largest sulfate loading event in the last 2500 years has long been thought to originate from the Kuwae caldera in Vanuatu. This event would be second only to the 1259 C.E. Samalas eruption in Indonesia[1] in sulfate emission, and would have had a similar to or possibly larger impact on the climate system than the 1815 C.E. Tambora eruption[2-4]. Ice core studies first suggested the occurrence of a large volcanic event around 1450 C.E. based on glaciochemical evidence in a South Pole ice core[5], a large electrical conductivity spike at 1460 C.E. in East Antarctica[6], and a sulfate increase between 1380–1460 C.E. in a South Pole ice core[7]. The Greenland ice sheet acidity record also shows a modest increase at 1450 C.E.[8], which was not thought to be a tropical volcanic event at the time. Later, Pang[9] tied the southern hemisphere ice core anomaly to a large volcanic event from Kuwae volcano in Vanuatu primarily utilizing historical reports and tree ring data. The Kuwae caldera collapse was first suggested to have occurred in 1425 C.E., based on historical reconstructions and radiocarbon dating[10]. The timing was shortly revised by Robin et al.[11] to 1452–1453 C.E. by combining 14C and ice core chronologies. Published tree ring frost damage[12] and density data[13] were interpreted as being the result of a strong cooling in the northern hemisphere in 1453 C.E., seemingly confirming 1452 C.E. as the start of the major volcanic eruption.

A number of Greenland and Antarctic ice core studies[2,3,14-16] accepted the 1452 C.E. age and Kuwae as the source caldera. At odds with this are proximal volcanological data collected at Tongoa and Tongariki islands[17], reported as a sequence of subaqueous deposits with no evidence of widespread fall deposits that would be expected for a large caldera-forming eruption. The same report also disputed early historical evidences, and pointed to early radiocarbon dates[17,18]. This information did not influence the attribution of the largest sulfate signal observed in the developed 2000-yr Law Dome ice core record around 1458 C.E. to the large caldera forming event from Kuwae[19].

Recently, the eVolv2k compilation of stratospheric volcanic injections listed the 1458 C.E. event as the third largest since 500 B.C.E., but Kuwae was not identified as a potential source[20]. Updated tree-ring based paleoclimate reconstructions strongly support a major cooling event in 1453 C.E. in the Northern Hemisphere[21], implying that a large tropical volcanic event occurred in 1452 C.E. There are no similar reconstructions yet available for the Southern Hemisphere.

Recent developments in methodology and instrumentation in the last decade have made it possible to detect and measure physicochemical properties of ultra fine tephra particles (cryptotephra) in various depositional media[22]. The new methodology has made it possible to correlate tephra deposits located thousands of kilometers from source volcanoes[23,24]. Here we test the proposed source of this large 1450s enigmatic sulfate signal by investigating the geochemical composition of cryptotephra (micron sized volcanic glass particles) extracted from ice sampled between 72.53–72.63 meters depth of the South Pole Ice Core (SPICEcore) from Antarctica (89.99°S, 98.16°W; Fig. 1). According to the SP19 timescale[25], developed using annual layer counting and volcanic signal synchronization with the WDC-06A ice core, this depth interval corresponds to 1458 C.E (Supplementary Fig. S1). To date, no ice-core-based studies have analyzed cryptotephra extracted from the mid-1450s interval. Cryptotephra geochemistry[23] is instrumental for independent verification of developed time scales and source confirmation of identified volcanic events, which is what this work was originally designed to do.

Figure 1

Location map with ice core sites and volcanic centers discussed in this paper. Ice core sites are represented as purple and orange dots. Different colors and sizes are associated with the sulfate flux ratios observed in the ice core records. Possible sources of volcanic products discussed in the paper are marked in red. Mapped using the Generic Mapping Tools software package version 5.4.5.[42]. Grounded ice boundary data from[43]. Sulfate flux values from[3,38,44] are provided in the Supplementary Table S1.

Results

The particles found in the SPICEcore are rhyolitic with ~75 wt.% SiO2, ~3 wt.% Na2O and K2O, and ~2 wt.% FeO. Figure 2 and Table S2 show the major element concentrations (expressed as oxide weight percent) for two particles measured by EMPA-WDS and for five particles measured by SEM-EDS. All individual point measurements were averaged after normalization of data from both analyses. In previous ice core cryptotephra work studies[26], a 2% uncertainty for SEM-EDS analyses was established and observed standard deviation of a secondary standard rhyolite glass fit within this parameterization (see Supplemental section for details).

Figure 2

Geochemical composition of individual glass shards extracted from the SPICEcore compared with possible volcanic sources using Harker variation diagrams. Open red triangles represent ice core volcanic glass particle data analyzed using SEM-EDS. The same volcanic glass particles analyzed using EPMA-WDS method are shown with filled red triangles. Symbol size is larger than analytical precision. Tabulated data with SPICEcore cryptotephra and referenced geochemical data are provided in the Supplementary Data section. A subset of this figure is shown in Supplementary Fig. S3.

We observed (Fig. 2) slight variations between the two measurement techniques and within populations measured using the same instrument. The SEM-EDS population variability, notably in K2O and CaO, can be explained by the particles being unpolished. Rough surface geometry is known to cause some geochemical variations[27]. Iverson and others[26] established that for some major element oxides, unpolished SEM-EDS measurements are robust enough to be useful for geochemical comparison despite a high standard deviation. In addition, the microprobe results show generally higher Na2O wt% compared to the SEM-EDS measurements. This is expected because microprobes are configured to measure Na2O in the first counting period during EPMA tephra analyses, whereas all elements are measured simultaneously during SEM-EDS analyses (longer counting time could result in Na loss on the SEM/EDS), and the SEM-EDS measurements are on unpolished particles, which increases variability of the results.

Discussion

The geochemical composition of volcanic glass particles extracted from 72.53–72.63 meters depth in the South Pole ice core is very distinct from Kuwae eruptive products. The Kuwae glass analyses published by Robin et al.[11] define a basic to intermediate geochemical composition (Fig. 2). The high SiO2 observed in ice core tephra particles is inconsistent with any known Kuwae eruptive products. An argument could be made that these particles represent a more evolved composition of the Kuwae glass; however, low K2O values do not support this (Fig. 2). The expected trend for the measured ice core cryptotephra SiO2/K2O ratio is in disagreement with values expected for Kuwae eruptive products, which suggests that the ice core cryptotephra originated from a different volcanic source. The composition of ice core cryptotephra is chemically distinct from Kuwae with a similarity coefficient (SC)[22,28,29] around 0.5 (see Table S3). In ideal cases, an SC exceeding 0.92 is considered to be a reasonable match[29]. Therefore, other volcanic sources for these particles must be considered.

Relatively large particle size (Fig. 3), implies rapid (within weeks after injection into the atmosphere) particle transport, so a range of regional volcanic sources are possible. Typical regional sources for tephra in Antarctic ice cores are Antarctic, sub-Antarctic, South American, and New Zealand volcanoes[30-32].

Figure 3

Selected cryptotephra particle images captured from 72.53–72.63 m SPICEcore sample. Each BSE image was taken at 15 kV and 0.5 nA on the Tescan Vega-II XMU before polishing and analysis by EPMA-WDS. Numbers in upper left corners represent the particle number with composition reported in Table S2.

The last 1000 years of the New Zealand volcanism record is fairly comprehensive[33], therefore missing a large eruption from a New Zealand source would be unlikely. One large New Zealand rhyolitic event (Fig. 2) was the Kaharoa eruption which occurred in 1314 C.E. from Taupo Volcanic Zone[34]. Because the tephra was extracted from a relatively modern part of the ice core record, the uncertainty in the layer counting-based SP19 timescale[25] is relatively low (±3 years). Therefore, the age alone makes Kaharoa tephra an unlikely candidate for the SPICEcore tephra deposits match. In a very comprehensive volcanism record from New Zealand there are no other identified eruptions from New Zealand around the 1450’s, which makes New Zealand volcanism an unlikely source for the tephra extracted from the SPICEcore, despite a very similar geochemical signature (SC range 0.7–0.9, see Table S3) of Kaharoa tephra from the eruption of Mt. Tarawera. The affiliation of tephra with the large sulfate increase also rules out potential transport of wind blown material from New Zealand.

The present characterization of South American and Antarctic volcanism in the 1450’s is less comprehensive than other localities (i.e. New Zealand), and a volcanic eruption from these regions could be a source for SPICEcore tephra deposits. Antarctic volcanism during the last 1000 years typically has a basic composition[35] and no known rhyolitic eruptions occurred in the 15 century in Antarctica. However, the Malta Plateau volcanic area is known for compositionally evolved lavas[35]. A comparison of a general geochemical suite of eruptive products from this area is shown in Fig. 2. It is not a perfect match because FeO values from Malta Plateau volcanic products are within a range of 3–10 wt.%, which is higher than the 1–2.5 wt.% observed in SPICEcore tephra. According to Thouret et al.[36], an eruption was observed in the 1450s C.E. at El Misti volcano in Peru. Geochemical data from this volcanic center are included in Fig. 2 for comparison. The SPICEcore tephra SiO2 and K2O geochemical fingerprint is inconsistent with the El Misti signature. A geochemical match (SC ranging from 0.7 to 0.8, Table S3) with Reclus tephra and evidence of recent activity of that volcano point to the possibility that a relatively large, unknown volcanic eruption from this volcanic center is responsible for this largest increase in sulfate in Antarctic ice cores in the last two millennia[37].

The highly variable magnitude and duration of the sulfate peak observed in many Antarctic ice core records (Figs 1, S1) favors a South American source volcanic eruption around 1458 C.E. The 1453 C.E. temperature anomalies in the Northern Hemisphere and the sulfate increase in Greenland ice cores around 1454 and 1460 C.E.[38] can possibly be attributed to different volcanic eruptions occurring in the northern hemisphere. Several moderate eruptions occurring in both hemispheres over a short time frame could have complex global effects on atmospheric circulation patterns and temperatures[39]. Location sources and potentially different timing for signals attributed to a large volcanic eruption from the Kuwae caldera call for revision of parameterization values for ice core based volcanic forcing time series used in climate models since the pioneering work of[40].

In addition, the 1458 C.E. Antarctic sulfate signal timing should not be used to adjust the Greenland sulfate signal originally observed at 1452 C.E.[3,12-14,21]. Therefore, the proposed use of a 1458 C.E. event for bipolar synchronization[1] should be reevaluated.

Conclusion

This study presents the first quantitative EPMA-WDS and EDS-based geochemical data set on volcanic glass particles extracted from the 72.53–72.63 m depth of the SPICEcore. These new data challenge the majority of previous studies that have suggested that the large sulfate peak present during the ~1450’s C.E. time interval is associated with a large tropical eruption from the Kuwae volcano in Vanuatu, while providing further evidence toward the emerging volcanological and historical data that question such attribution. An accurate determination of the volcanic eruption source and magnitude is important and data presented here suggest that Kuwae has been incorrectly used in past climate model simulations as a source of the ca. 1450s C.E. large sulfate signal in polar ice cores. Although no volcanic source is unambiguously identified in this research, particle size and sulfate appearance in many Antarctic ice core records suggest the eruption is likely affiliated with a tropospherically transported aerosol cloud. A reasonably good geochemical match of SPICEcore tephra with the Reclus volcano and relatively recent, compositionally similar reworked tephra in Laguna Potrok Aike, provide a basis for a preliminary conclusion that possibly an unknown Reclus eruption is responsible for the volcanic signal in Antarctic ice cores. These data contribute to the Antarctic tephra framework needed for future tephrochronological studies, and suggest that paleoclimate models using Kuwae as a climate-forcing event for the ca. 1450 C.E. climate anomaly need to be revised.

Methods

The sampled interval targeted to locate and extract tephra associated with the 1450s sulfate anomaly was initially selected based on the preliminary version of the SP19 timescale. We sampled two depth intervals: 72.43–72.525 m (sample AntT-335) and 72.525–72.63 m (AntT-336). The sampling was guided by the age of an electrical conductivity (ECM) anomaly in the SPICEcore[41]. In the final version of the SP19 time scale, developed using ECM, visual layer counting, sodium and magnesium ions measurements, interpolation between 251 stratigraphic volcanic tie points to the WD2014 chronology, and previously published South Pole ice core records, the increase in electrical conductivity starts around 1459 (±2) C.E. and lasts for two years[25]. The remarkable manifestation of the sulfate-based volcanic record at the South Pole is attributable to low contributions of sulfate from dust, biogenic, and anthropogenic sources as well as stratospheric air masses[5]. Supplementary Figs S1 and S2 highlight the studied SPICEcore interval and correlated volcanic signals in other Antarctic ice core records.

Using improved ice core tephra mounting methodologies (see online Supplementary section for more details), five tephra particles (ca. 5–50 μm in size) were captured and analyzed using a Tescan Vega-II XMU scanning electron microscope (SEM) equipped with an EDAX Apollo SSD40 energy dispersive spectrometer (EDS). The particles were analyzed on a single spot with a 140 pA beam for 100 seconds of live time at 15 kV of accelerating voltage. Major element oxides were calculated and adjusted using the semi-quantitative EDAX Genesis PhiRhoZ internal quantification procedure and the USNM 72854 VG-568 rhyolitic glass standard. The two largest particles were then polished and subsequently reanalyzed using a Cameca SX-100 electron microprobe equipped with five wavelength dispersive spectrometers (WDS). All measurements were made using instrumentation at the University of Maine. Analytical settings are discussed in Supplementary methods.

Data Availability

All analytical cryptotephra geochemical data shown in Figs 2 and S3 are available in the Supplement Table S2. The data that support findings of this study are available within the paper and its supplementary files. All tephra related microbeam measurements will be archived with the publication of this paper at the U.S. Antarctic Program Data Center (USAP-DC) http://www.usap-dc.org.

Supplementary materials

Table S1

Table S2

Table S3

Dataset Figure 1

Dataset Figure 2 and S3