Description: We present continuous records of 1 cm resolved sulfate concentrations, insoluble particle number and mass concentration, and liquid conductivity using a Continuous Flow Analysis (CFA) system (Bigler et al., 2002; Bigler et al., 2011; Erhardt et al., 2022) from the North Greenland Ice Core Project (NGRIP; 75.10°N, 42.33°W; 2941 m a.s.l.) and estimated volcanic sulfate mass depositions for the time period 79.14 and 80.48 ka BP on the GICC05modelext chronology, transferred to the AICC2012 chronology (Veres et al., 2013; Seierstad et al., 2014) using a common volcanic marker in EPICA Dome C ice core dated 79.51 ka BP. The reconstruction is based on sulfate measurements employing high-resolution continuous flow analysis. Volcanic eruptions are detected when annual sulfate concentrations exceeded the background concentrations + 4 times the median of the absolute deviation. Background concentrations are estimated using a 181-point running median. Volcanic sulfate deposition rates are calculated by subtracting the background concentrations from total sulfate concentrations using thinning corrected estimates of mean ice accumulation rates at the ice-core site.

Description: The West Antarctic Ice Sheet (WAIS) Divide deep ice core WD2014 chronology, consisting of ice age, gas age, delta-age and uncertainties therein. 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; Sigl et al., 2015, Sigl et al., 2016) have been dated using annual-layer counting based on counting of annual layers observed in the chemical, dust and electrical conductivity records. The measurements were interpreted manually and with the aid of two automated methods. We validated the chronology by comparing of the cosmogenic isotope records of 10Be from WAIS Divide and 14C for IntCal13. We demonstrated that over the Holocene WD2014 was consistently accurate to better than 0.5% of the age. The chronology for the deep part of the core (below 2850m; 67.8-31.2 ka BP; Buizert et al., 2015) 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 (Delta age) using a combination of firn densification modeling, ice-flow modeling, and a data set of d15N-N2, a proxy for past firn column thickness. The largest Delta age at WD occurs during the Last Glacial Maximum, and is 525 +/- 120 years. We synchronized 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.

Description: We present continuous records of ca. 2 cm resolved sulfate concentrations using a Continuous Flow Analysis (CFA) system (Röthlisberger et al., 2000; Severi et al., 2007; Wegner et al., 2015) from the EPICA Dronning Maud Land (EDML; 75.00°S, 00.07°E; 2892 m a.s.l.) ice core and estimated volcanic sulfate mass depositions for the time period 78.97 to 79.97 ka BP on the AICC2012 chronology (Veres et al., 2013) using a common volcanic marker in EDC dated 79.51 ka BP. The reconstruction is based on sulfate measurements employing high-resolution continuous flow analysis coupled to Fast Ion Chromatography. Volcanic eruptions are detected when annual sulfur concentrations exceeded the background concentrations + 4 times the median of the absolute deviation. Background concentrations are estimated using a 101-point running median. Volcanic sulfate deposition rates are calculated by subtracting the background concentrations from total sulfate concentrations using thinning corrected estimates of mean ice accumulation rates at the ice-core site.

Description: We used four polar ice cores - two from Greenland and two from Antarctica - to investigate volcanic sulphate deposition on polar ice sheets between 13,200 and 12,800 years BP. From Greenland, we employed discrete sulphate measurements from the Greenland Ice Sheet Project Two (GISP2; 72.97°N, 38.80°W) ice core (Mayewski et al., 1997) as well as sulphate measurements using Continuous Flow Analysis (CFA) (Bigler et al., 2002; Bigler et al., 2011) from the North Greenland Ice Core Project (NGRIP; 75.10°N, 42.33°W; 2941 m a.s.l.). From Antarctica, we used Fast Ion Chromatography (FIC) measurements from the EPICA Dronning Maud Land (EDML; 75.00°S, 00.07°E; 2892 m a.s.l.) ice core (Severi et al., 2007) as well as sulphur measurements (converted to sulphate) using Continuous Flow Analysis (CFA) coupled to an inductively coupled plasma mass spectrometer (ICPMS) (McConnell et al., 2017) from the WAIS (West Antarctic Ice Sheet) Divide ice core project (WD; 79.48°S, 112.11°W, 1766 m a.s.l.) (Sigl et al., 2016). The time resolution of the data range between multi-annual (i.e. 4 years) for GISP2 to sub-annual for NGRIP and WD. Based on these ice cores and by applying the age synchronization lattice previously developed for these cores (Seierstad et al., 2014; Buizert et al., 2018; Svensson et al., 2020) with linear interpolation between the volcanic markers, we reconstructed volcanic sulphate deposition over Greenland and Antarctica using established methods (Sigl et al., 2014; Sigl et al., 2015). Using the methodology developed and applied to a similar network of ice cores over the past 2,500 years (Toohey & Sigl 2017) we further estimated stratospheric sulphur injection (SSI) from 30 eruptions with SSI in excess of 1 Tg S.

Description: We present continuous records of 1 cm resolved sodium, sulfur, non-sea-salt sulfur, size-resolved insoluble particle concentrations and liquid conductivity using a Continuous Flow Analysis (CFA) system (McConnell et al., 2017) from the Greenland Ice Sheet Project Two (GISP2; 72.97°N, 38.80°W) ice core and estimated volcanic sulfate mass depositions for the depth interval 2635-2638 m (equivalent to an age of c. 80 ka BP on GICC05modelext (Seierstad et al., 2014) and 79.5 ka BP on AICC2012 (Veres et al., 2013). The reconstruction is based on sulfur measurements employing high-resolution continuous flow analysis coupled to inductively coupled plasma mass spectrometry performed at the Desert Research Institute (Reno, NV, USA) under (class 100) cleanroom conditions. Volcanic eruptions are detected when annual sulfur concentrations exceeded the background concentrations + 4 times the median of the absolute deviation. Background concentrations are estimated using a 101-point running median. Volcanic sulfate deposition rates are calculated by subtracting the background concentrations from total sulfate equivalent (i.e. sulfur x 3) concentrations using thinning corrected estimates of mean ice accumulation rates at the ice-core site.

Description: Sulfate concentration and triple sulfur isotope (32S, 33S, 34S) data measured in TUNU2013 ice-core samples and calculated background-corrected isotope composition (Gabriel et al., in review). S-isotope analyses were made on 65 discrete samples (cross-sections of 6 cm2) cut from archived ice-core sections from TUNU2013, including prior background and full acid deposition lasting over 12 years. The sample resolution over the entire acid deposition event is 3 cm corresponding to a nominal 3-to-4-month age resolution. Sulfate concentration was measured by ion chromatography on the discrete samples, and the sulfate was purified from the melted ice using anion exchange columns. Triple sulfur isotopes (32S, 33S, and 34S) were measured on the samples using a Neptune Plus multi-collector inductively-coupled mass spectrometer (MC-ICP-MS) at the St Andrews Isotope Geochemistry Lab (STAiG lab) (Burke et al., 2019; Burke et al., 2023) and are reported as δ34S and D33S relative to Vienna-Canyon Diablo Troilite (V-CDT), where δxS =(xS/32S)sample/(xS/32S)V-CDT−1. A sample is considered to have a mass-independent fractionation signature if it has a nonzero value of D33S = δ33S -((δ34S + 1)^0.515 − 1), outside of 2σ uncertainty. Full procedural blanks and an in-house secondary standard (Switzer Falls) were processed alongside samples. All data was blank corrected for the process blanks, and uncertainties were propagated with Monte Carlo simulations. The isotopic composition of the volcanic sulfate (δ34Svolc and D33Svolc) was calculated using isotope mass balance and the concentration and isotopic composition of background ice from TUNU2013 taken prior and after the volcanic peaks in sulfate. Data is shown on the NS1-2011 chronology (Sigl et al., 2015). Isotopic composition of the volcanic sulfate is only provided for samples with more than 65% volcanic sulfate following Burke et al. (2019).

Description: Annual volcanic sulfate mass deposition rates from the West Antarctic Ice Sheet (WAIS) Divide ice core project for the time period 12.68-13.20 ka BP on the annual-layer counted WD2014 chronology and transfered to the Greenland GICC05 chronology using common volcanic markers. The reconstruction is based on sulfur measurements employing high-resolution continuous flow analysis coupled to inductively coupled plasma mass spectrometry performed at the Desert Research Institute (Reno, NV, USA). Volcanic eruptions are detected when annual sulfur concentrations exceeded the background concentrations + 2 times the median of the absolute deviation. Background concentrations are estimated using a 101-year running median. Volcanic sulfate deposition rates are calculated by subtracting the background concentrations from total sulfate concentrations using thinning corrected estimates of mean accumulation rates at the ice-core site.

Description: Annual-resolved non-sea-salt sulfur concentrations from a four ice-core stack, NEEM-2011-S1 (Sigl et al., 2013), NGRIP1 (Plummer et al., 2012), TUNU2013 (Sigl et al., 2015) and B19, between 1731 and 1996 CE including volcanic samples and with volcanic samples replaced by a 11-year running median.

Description: Annual-resolved sulfur and non-sea-salt sulfur concentrations and inferred volcanic sulfate deposition rates from the six ice cores NEEM-2011-S1 (Sigl et al., 2013), NGRIP1 (Plummer et al., 2012), NGRIP2 (McConnell et al., 2018), TUNU2013 (Sigl et al., 2015) and B19 between 699 and 1001 CE.

Description: Annual-resolved previously published volcanic sulfate deposition rates from two ice cores NEEM-2011-S1 (Sigl et al., 2013) and NGRIP1 (Plummer et al., 2012) between 699 and 1001 CE.