Wunderlich, F., & Mitchell, D. M. (217). Revisiting the observed surface climate response to large volcanic eruptions. Atmospheric Chemistry and Physics, 17(1), 485-499. DOI: 1.5194/acp-17-485-217 Publisher's PDF, also known as Version of record License (if available): CC BY Link to published version (if available): 1.5194/acp-17-485-217 Link to publication record in Explore Bristol Research PDF-document This is the final published version of the article (version of record). It first appeared online via EGU at http://www.atmos-chem-phys.net/17/485/217. Please refer to any applicable terms of use of the publisher. University of Bristol - Explore Bristol Research General rights This document is made available in accordance with publisher policies. Please cite only the published version using the reference above. Full terms of use are available: http://www.bristol.ac.uk/pure/about/ebr-terms
Supplement of Atmos. Chem. Phys., 17, 485 499, 217 http://www.atmos-chem-phys.net/17/485/217/ doi:1.5194/acp-17-485-217-supplement Author(s) 217. CC Attribution 3. License. Supplement of Revisiting the observed surface climate response to large volcanic eruptions Fabian Wunderlich and Daniel M. Mitchell Correspondence to: Fabian Wunderlich (fabian.wunderlich@met.fu-berlin.de) The copyright of individual parts of the supplement might differ from the CC-BY 3. licence.
ERA-2C JRA-55 NOAA-CIRES 2CR ERA-4 MERRA NCEP CFSR ERA-Interim NCEP/NCAR I JRA-25 NCEP/DOE II -5. -4. -3. -2. -1.. 1. 2. 3. 4. 5. -4.5-3.5-2.5-1.5 -.5.5 1.5 2.5 3.5 4.5 Figure S1. SLP anomalies (hpa) in the Northern Hemisphere in the first post volcanic winters (DJF) averaged over the two volcanic eruptions of Mt. Pinatubo (1991) and El Chichón (1982) for all 1 reanalysis datasets. Anomalies are calculated with respect to the mean for the years 1979-212, excluding the following two years after the eruptions. Single diagonal lines correspond to the 9% and double diagonal lines to the 95% confidence level obtained with a Monte Carlo test of two independent samples. 2
ERA-2C JRA-55 NOAA-CIRES 2CR ERA-4 MERRA NCEP CFSR ERA-Interim NCEP/NCAR I JRA-25 NCEP/DOE II -5. -4. -3. -2. -1.. 1. 2. 3. 4. 5. -4.5-3.5-2.5-1.5 -.5.5 1.5 2.5 3.5 4.5 Figure S2. As Figure 1, but for the second post volcanic winter. 3
a) December b) January 65 65 52 52 39 Santa Maria, 192 39 Santa Maria, 192 26 26 13 13-4 -3-2 -1 1 2 3 4 c) February -4-3 -2-1 1 2 3 4 d) Winter (DJF) 65 52 39 Santa Maria, 192 65 52 39 NOAA 2CR HadSLP All Reanal. Santa Maria, 192 26 26 13 13-4 -3-2 -1 1 2 3 4-4 -3-2 -1 1 2 3 4 Figure S3. Monthly North Atlantic Oscillation (NAO) index of the second winter after 5 volcanic eruptions (a-c) and the winter mean (d), calculated with HadSLP2 observation data and reanalysis data. All data is calculated with respect to the mean for the years 1979-212, excluding the following two years after the two eruption. The EOF is calculated over the period 1979-212 for every product separately. The histogram shows the NAO index of the 163 years of observation data (185-212). The blue lines show the reanalysis data spread of the NAO index for the winter after Mt. Pinatubo and El Chichón. The orange lines show the 9% ensemble spread of the NAO after all 5 volcanic eruptions calculated with the NCEP-2CR dataset. 3 2 1 Krakatau Tarawera Bandaisan Colima Pelee/Soufr. Ksudach Novarupta Colima Quizapu Kharimkotan Fuego Mt. Pinatubo St. Helens El Chichon -1-2 Bezymianny Mt. Agung Askja Santa Maria Fernandina -3 1871 1881 1891 191 1911 1921 1931 1941 1951 1961 1971 1981 1991 21 Figure S4. Winter mean NAO index from 1871 until present, calculated with HadSLP2 observation data and NCEP-2CR ensemble data. All data is calculated with respect to the mean for the years 1979-212, excluding the following two years after the two eruption. The EOF is calculated over the period 1979-212. Every single thin coloured line represents one ensemble of the NCEP-2CR reanalysis. Thick coloured line shows the NCEP-2CR ensemble mean and black line shows HadSLP2 data. All volcanoes from Figure 2 are indicated by a grey line at the first winter after the eruption or at the second winter after the eruptions in the case of Fuego and Santa María. All volcano names which are above the zero line were followed by a positive NAO in the first (or second) winter after the eruption according to the HadSLP2 data and all volcano names which are below the zero line were followed by a negative NAO. 4
a) HadCRUT4, 1st year b) HadCRUT4, 2nd year c) Reanalysis Mean, 1st year d) Reanalysis Mean, 2nd year -3. -2.4-1.8-1.2 -.6..6 1.2 1.8 2.4 3. -2.7-2.1-1.5 -.9 -.3.3.9 1.5 2.1 2.7 Figure S5. Observed and multi-reanalysis mean TAS anomalies (K) averaged over the first year (a and c) and second year (b and d) after the eruption and over the three volcanic eruptions Mt. Pinatubo (1991), El Chichón (1982) and Mt. Agung (1963). a) and b) show the mean anomaly for the HadCRUT4 observations and c) and d) show respectively the mean anomaly of the five reanalysis datasets, containing this period (Table 1). Anomalies are calculated with respect to the average over the years 1961-199, consistent to the HadCRUT4 dataset. Single diagonal lines correspond to the 9% and double diagonal lines to the 95% confidence level obtained by a Monte Carlo test of three independent samples. The Figure is similar to Figure 6 with the difference, that no ENSO removal is applied to the data. 5
ERA-2C JRA-55 NOAA-CIRES 2CR ERA-4 MERRA NCEP CFSR ERA-Interim NCEP/NCAR I JRA-25 NCEP/DOE II -2. -1.6-1.8-1.2-1.4 -.8-1. -.4 -.6. -.2.4.2.8.6 1.2 1. 1.6 1.4 2. 1.8 Figure S6. TAS anomalies (K) averaged over the first year after the eruption and over the two volcanic eruptions Mt. Pinatubo (1991) and El Chichón (1982) for all 1 reanalysis products. Anomalies are calculated with respect to the average over the years 1979-212. Single diagonal lines correspond to the 9% and double diagonal lines to the 95% confidence level obtained with a Monte Carlo test of two independent samples. 6
a) First year after eruption b) Second year after eruption.3 HadCRUT4 NOAA 2CR 5 Reanalyses CMIP5 his CMIP5 vol.3 probability density.2.1 Santa Maria, 192 Santa Maria, 192.2.1 probability density -.6 -.3..3.6 TAS anomalies -.6 -.3..3.6 TAS anomalies Figure S7. Global TAS anomalies (K) averaged and over the first year (a) and second year (b) after five volcanic eruptions, calculated with HadCRUT4 observation data and reanalysis data. Anomalies are calculated with respect to the average over the years 1961-199, consistent to the HadCRUT4 dataset. The histogram shows the anomalies of 153 years of observation data (1855-27). The first and last five years are leaved out because of the trend removal. The blue lines show the reanalysis data spread of the TAS anomalies for the years after all eruptions except Krakatau and Santa María. The orange lines show the 9% ensemble spread of TAS after all five volcanic eruptions calculated with the NOAA-2CR dataset. The mean response is indicated by the dot between the whiskers. 7