Volcanic eruptions and climate

Understanding natural causes of climate change is vital to evaluate the relative impacts of human pollution and land surface modification on climate. We have investigated one of the most important natural causes of climate change, volcanic eruptions, by using 54 ice core records from both the Arctic and Antarctica. Our recently collected suite of ice core data, more than double the number of cores ever used before, reduces errors inherent in reconstructions based on a single or small number of cores, which enables us to obtain much higher accuracy in both detection of events and quantification of the radiative effects. We extracted volcanic deposition signals from each ice core record by applying a high‐pass loess filter to the time series and examining peaks that exceed twice the 31‐year running median absolute deviation. We then studied the spatial pattern of volcanic sulfate deposition on Greenland and Antarctica and combined this knowledge with a new understanding of stratospheric transport of volcanic aerosols to produce a forcing data set as a function of month, latitude, and altitude for the past 1500 years. We estimated the uncertainties associated with the choice of volcanic signal extraction criteria, ice core sulfate deposition to stratospheric loading calibration factor, and the season for the eruptions without a recorded month. We forced an energy balance climate model with this new volcanic forcing data set, together with solar and anthropogenic forcing, to simulate the large‐scale temperature response. The results agree well with instrumental observations for the past 150 years and with proxy records for the entire period. Through better characterization of the natural causes of climate change, this new data set will lead to improved prediction of anthropogenic impacts on climate. The new data set of stratospheric sulfate injections from volcanic eruptions for the past 1500 years, as a function of latitude, altitude, and month, is available for download in a format suitable for forcing general circulation models of the climate system.

Long-term changes of the diurnal temperature cycle: implications about mechanisms of global climate change

<strong class="journal-contentHeaderColor">Abstract.</strong> 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 11500 year BP to 1900 CE), constituting an extension of the previous record by 7000 years. The database incorporates new-generation ice-core aerosol records with sub-annual temporal resolution and 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 previous 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 TgS, of which 329 (39 %) are located in the low latitudes with bipolar sulfate deposition, 426 (50 %) are located in the Northern Hemisphere (NH) extratropics and 88 (10 %) are located in the Southern Hemisphere (SH) extratropics. The spatial distribution of reconstructed eruption locations is in agreement with prior reconstructions for the past 2,500 years, and follows the global distribution of landmasses. In total, these eruptions injected 7410 TgS in the stratosphere, for which tropical eruptions accounted for 70 % and NH extratropics for 25 %. 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&ndash;7000 BCE) experienced a higher number of volcanic eruptions (+16 %) and cumulative VSSI (+86 %) compared to the past 2,500 years. This increase coincides with the rapid retreat of ice sheets during deglaciation, providing context for potential future increases of volcanic activity in regions under projected glacier melting in the 21st century. The reconstructed VSSI and SAOD data are available at <a href="https://doi.pangaea.de/10.1594/PANGAEA.928646" target="_blank" rel="noopener">https://doi.pangaea.de/10.1594/PANGAEA.928646</a> (Sigl et al., 2021).

Abstract. Large tropical (TROP) volcanic eruptions can influence North Atlantic climate by inducing a positive shift of the North Atlantic Oscillation (NAO), typically resulting in winter warming across northern Eurasia. However, these changes remain highly uncertain, as they may coincide with strong internal variability in Northern Hemisphere wintertime climate. In contrast, Northern Hemisphere extratropical (NHET) eruptions are proposed to have opposite impacts, but they have been comparatively less studied, and large uncertainties remain regarding the ability of climate models to capture volcanic responses. This study examines winter North Atlantic climate responses to TROP and NHET eruptions by comparing temperature and atmospheric circulation patterns from last millennium simulations with multiple proxy-based reconstructions. We find distinct differences in NAO-related climate changes in reconstructions, with TROP eruptions followed by a shift towards positive NAO and NHET eruptions associated with a negative NAO. In comparison, modelled responses exhibit a wide spread with strong dependence on the choice of volcanic forcing dataset. Notably, simulations using the latest volcanic forcing data show improved agreement with reconstructions, particularly for TROP eruptions. This model-proxy agreement provides a useful basis for investigating the mechanisms that drive positive NAO responses after TROP eruptions. However, the simulated impacts of NHET eruptions are less consistent and remain unclear. These results highlight the importance of improved volcanic forcing datasets, refined paleoclimate reconstructions, and robust statistical approaches to better constrain uncertainties in assessing volcanic impacts on North Atlantic climate.

&amp;lt;p&amp;gt;The cumulative impacts of frequent moderate-magnitude eruptions on stratospheric aerosols were identified among the factors in recent decadal climate trends. Moderate volcanic eruptions are a recurrent source of sulfur dioxide (SO2) in the Upper Troposphere and Lower Stratosphere (UTLS) region and the resulting formation of sulfuric acid aerosol particles from the SO2 emitted provides sites for chemical reactions leading to enhancement of stratospheric optical depth (SAOD) and ozone depletion. Modelling properly the volcanic aerosol content and its evolution in this region is important for radiative impact issues. In this work, we explore the variability of the tropical UTLS aerosol content between 2013 and 2019, a period which was particularly impacted by moderate tropical and mid-latitude volcanic eruptions. For that purpose, space-borne observations from OMPS (version 2, datasets from GES DISC), and IASI, together with simulations by the Whole Atmosphere Community Climate Model (WACCM) coupled with the Community Aerosol and Radiation Model for Atmospheres (CARMA), are used. Different model sensitive experiments, particularly for the injection altitude and timing, have been conducted to evaluate how the model captures the aerosol plume in terms of content, optical and microphysical properties, transport and residence time. We find that the decay of the Calbuco and Kelud plumes observed by OMPS version 2 is well reproduced by the model. Comparisons with unique datasets in the tropical southern hemisphere from the NDACC Ma&amp;amp;#239;do observatory (Reunion Island, France, 20.5&amp;amp;#176;S, 55.5&amp;amp;#176;E) show good agreement between the lidar SAOD observations and WACCM-CARMA SAOD simulations although we observe a difference in the altitude of the maximum aerosol concentration between the model and the in situ profile after Calbuco eruption in April 2015. A particular focus is also made on recent eruptions like Raikoke, Ambae and Ulawun. The plume of the Ambae volcano (15&amp;amp;#176;S, 167&amp;amp;#176;E) which erupted in July 2018 is shown to propagate to the northern hemisphere with some influence until summer 2019 in the Asian monsoon region. For the year 2019, we investigate how the Ulawun (5&amp;amp;#176;S, 151&amp;amp;#176;E; ~0.14 Tg of SO2) tropical eruption and the Raikoke mid-latitude eruption (48&amp;amp;#176;N, 153&amp;amp;#176;E; ~1.5Tg of SO2), have influenced the aerosol burden in the tropics.&amp;lt;/p&amp;gt;

Previous studies in the literature show that the annual cycle of surface air temperature (SAT) is changing in both amplitude and phase, and the SAT departures from the annual cycle are long-term correlated. However, the classical definition of temperature anomalies is based on the assumption that the annual cycle is constant, which contradicts the fact of changing annual cycle. How to quantify the impact of the changing annual cycle on the long-term correlation of temperature anomaly variability still remains open. In this paper, a recently developed data adaptive analysis tool, the nonlinear mode decomposition (NMD), is used to extract and remove time-varying annual cycle to reach the new defined temperature anomalies in which time-dependent amplitude of annual cycle has been considered. By means of detrended fluctuation analysis, the impact induced by inter-annual variability from the time-dependent amplitude of annual cycle has been quantified on the estimation of long-term correlation of long historical temperature anomalies in Europe. The results show that the classical climatology annual cycle is supposed to lack inter-annual fluctuation which will lead to a maximum artificial deviation centering around 600 days. This maximum artificial deviation is crucial to defining the scaling range and estimating the long-term persistence exponent accurately. Selecting different scaling range could lead to an overestimation or underestimation of the long-term persistence exponent. By using NMD method to extract the inter-annual fluctuations of annual cycle, this artificial crossover can be weakened to extend a wider scaling range with fewer uncertainties.

The annual energy budget of the Arctic Ocean is characterized by a net heat loss at the air-sea interface that is balanced by oceanic heat transport into the Arctic. The energy loss at the air-sea interface is due to the combined effects of radiative, sensible, and latent heat fluxes. The inflow of heat by the ocean can be divided into two components: the transport of water masses of different temperatures between the Arctic and the Atlantic and Pacific Oceans and the export of sea ice, primarily through Fram Strait. Two 150-year simulations (1950-2099) of a global climate model are used to examine how this balance might change if atmospheric greenhouse gases (GHGs) increase. One is a control simulation for the present climate with constant 1950 atmospheric composition, and the other is a transient experiment with observed GHGs from 1950 to 1990 and 0.5% annual compounded increases of CO2 after 1990. For the present climate the model agrees well with observations of radiative fluxes at the top of the atmosphere, atmospheric advective energy transport into the Arctic, and surface air temperature. It also simulates the seasonal cycle and summer increase of cloud cover and the seasonal cycle of sea-ice cover. In addition, the changes in high-latitude surface air temperature and sea-ice cover in the GHG experiment are consistent with observed changes during the last 40 and 20 years, respectively. Relative to the control, the last 50-year period of the GHG experiment indicates that even though the net annual incident solar radiation at the surface decreases by 4.6 W(per square meters) (because of greater cloud cover and increased cloud optical depth), the absorbed solar radiation increases by 2.8 W(per square meters) (because of less sea ice). Increased cloud cover and warmer air also cause increased downward thermal radiation at the surface so that the net radiation into the ocean increases by 5.0 Wm-2. The annual increase in radiation into the ocean, however, is compensated by larger increases in sensible and latent heat fluxes out of the ocean. Although the net energy loss from the ocean surface increases by 0.8 W (per square meters), this is less than the interannual variability, and the increase may not indicate a long-term trend. The seasonal cycle of heat fluxes is significantly enhanced. The downward surface heat flux increases in summer (maximum 2 of 19 W per square meters or 23% in June) while the upward heat flux increases in winter (maximum of 16 W per square meters or 28% in November). The increased downward flux in summer is due to a combination of increases in absorbed solar and thermal radiation and smaller losses of sensible and latent heat. The increased heat loss in winter is due to increased sensible and latent heat fluxes, which in turn are due to reduced sea-ice cover. On the other hand, the seasonal cycle of surface air temperature is damped, as there is a large increase in winter temperature but little change in summer.

We examine whether significant changes in ocean temperatures can be detected in recent decades and if so whether they can be attributed to anthropogenic or natural factors. We compare ocean temperature changes for 1960–2005 in four observational datasets and in historical simulations by atmosphere-ocean general circulation models (AOGCMs) from the Coupled Model Intercomparison Project phase 5 (CMIP5). Observations and CMIP5 models show that the upper 2000 m has warmed with a signal that has a well-defined geographical pattern that gradually propagates to deeper layers over time. Greenhouse gas forcing has contributed most to increasing the temperature of the ocean, a warming which has been offset by other anthropogenic forcing (mainly aerosols), and volcanic eruptions which cause episodic cooling. By characterizing the ocean temperature change response to these forcings we construct multi-model mean fingerprints of time-depth changes in temperature and carry out two detection and attribution analysis. We consider first a two-signal separation into anthropogenic and natural forcings. Then, for the first time, we consider a three signal separation into greenhouse gas, anthropogenic aerosols and natural forcings. We show that all three signals are simultaneously detectable. Using multiple depth levels decreases the uncertainty of the results. Limiting the observations and model fields to locations where there are observations increases the detectability of the signal.

The pH of polar ice is important for the stability and mobility of impurities in ice cores and can be strongly influenced by volcanic eruptions or anthropogenic emissions. We present a simple optical method for continuous determination of acidity in ice cores based on spectroscopically determined color changes of two common pH-indicator dyes, bromophenol blue, and chlorophenol red. The sealed-system method described here is not equilibrated with CO2, making it simpler than existing methods for pH determination in ice cores and offering a 10-90% peak response time of 45 s and a combined uncertainty of 9%. The method is applied to Holocene ice core sections from Greenland and Antarctica and compared to standard techniques such as electrical conductivity measurement (ECM) conducted on the solid ice, and electrolytic meltwater conductivity, EMWC. Acidity measured in the Greenland NGRIP ice core shows good agreement with acidity calculated from ion chromatography. Conductivity and dye-based acidity Hdye+ are found to be highly correlated in the Greenland NEGIS firn core (75.38° N, 35.56° W), with all signals greater than 3σ variability coinciding with either volcanic eruptions or possible wild fire activity. In contrast, the Antarctic Roosevelt Island ice core (79.36° S, 161.71° W) features an anticorrelation between conductivity and Hdye+, likely due to strong influence of marine salts.

Fe2+ in ice cores as a new potential proxy to detect past volcanic eruptions

Large volcanic eruptions occurring in the last glacial period can be detected in terms of their deposited sulfuric acid 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 ka years before AD 2000 (b2k). The ice cores are synchronized over most of the investigated interval 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 to a large variability in the sulfate background signal, particularly in the Greenland glacial climate, we only detect 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 740 eruptions in Antarctica within the 51 ka period – where the sulfate deposition of 85 eruptions is defined at both poles (bipolar eruptions). Based on the relative Greenland and Antarctic sulfate deposition, we estimate the latitudinal band of the bipolar eruptions and assess their approximate climatic forcing based on established methods. The climate forcing of the five largest eruptions is estimated to be higher than −70 W m−2. Twenty-seven 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 VEI-7 Tambora eruption that occurred in Indonesia in 1815 AD. The frequency of eruptions larger than the typical VEI-7 (VEI-8) eruption by the comparison of sulfur emission strength is found to be 5.3 (7) times higher than estimated from geological evidence. 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. These volcanoes documented in ice cores provide atmospheric sulfate burden and climate forcing for further research on climate impact and understanding the mechanism of the Earth system.

Abstract. 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.

Abstract. 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.

Abstract. The major volcanic eruption of Mount Pinatubo in 1991 has been shown to have significant effects on stratospheric chemistry and ozone depletion even at midlatitudes. Since then, only moderate but recurrent volcanic eruptions have modulated the stratospheric aerosol loading and are assumed to be one cause for the reported increase in the global aerosol content over the past 15 years. This particularly enhanced aerosol context raises questions about the effects on stratospheric chemistry which depend on the latitude, altitude and season of injection. In this study, we focus on the midlatitude Sarychev volcano eruption in June 2009, which injected 0.9 Tg of sulfur dioxide (about 20 times less than Pinatubo) into a lower stratosphere mainly governed by high-stratospheric temperatures. Together with in situ measurements of aerosol amounts, we analyse high-resolution in situ and/or remote-sensing observations of NO2, HNO3 and BrO from balloon-borne infrared and UV–visible spectrometers launched in Sweden in August–September 2009. It is shown that differences between observations and three-dimensional (3-D) chemistry-transport model (CTM) outputs are not due to transport calculation issues but rather reflect the chemical impact of the volcanic plume below 19 km altitude. Good measurement–model agreement is obtained when the CTM is driven by volcanic aerosol loadings derived from in situ or space-borne data. As a result of enhanced N2O5 hydrolysis in the Sarychev volcanic aerosol conditions, the model calculates reductions of ∼ 45 % and increases of ∼ 11 % in NO2 and HNO3 amounts respectively over the August–September 2009 period. The decrease in NOx abundances is limited due to the expected saturation effect for high aerosol loadings. The links between the various chemical catalytic cycles involving chlorine, bromine, nitrogen and HOx compounds in the lower stratosphere are discussed. The increased BrO amounts (∼ 22 %) compare rather well with the balloon-borne observations when volcanic aerosol levels are accounted for in the CTM and appear to be mainly controlled by the coupling with nitrogen chemistry rather than by enhanced BrONO2 hydrolysis. We show that the chlorine partitioning is significantly controlled by enhanced BrONO2 hydrolysis. However, simulated effects of the Sarychev eruption on chlorine activation are very limited in the high-temperature conditions in the stratosphere in the period considered, inhibiting the effect of ClONO2 hydrolysis. As a consequence, the simulated chemical ozone loss due to the Sarychev aerosols is low with a reduction of −22 ppbv (−1.5 %) of the ozone budget around 16 km. This is at least 10 times lower than the maximum ozone depletion from chemical processes (up to −20 %) reported in the Northern Hemisphere lower stratosphere over the first year following the Pinatubo eruption. This study suggests that moderate volcanic eruptions have limited chemical effects when occurring at midlatitudes (restricted residence times) and outside winter periods (high-temperature conditions). However, it would be of interest to investigate longer-lasting tropical volcanic plumes or sulfur injections in the wintertime low-temperature conditions.

The June 15, 1991 Mount Pinatubo eruption was a large but relatively shortlived shock to the Earth’s atmosphere. It thus provided an excellent opportunity to study the workings of the climate system, to test climate models, and to examine the impacts of climate change on life. The largest eruption of the 20th Century inspired a large amount of research on the connection between volcanic eruptions and the Earth’s atmosphere in the 12 years since that eruption, as exemplified by the chapters in this book. Here several additional examples of our new understanding of these connections are presented. While the global cooling after Pinatubo was not surprising, the observed winter warming over Northern Hemisphere continents in the two winters following the eruption is now understood as a dynamic response to volcanically produced temperature gradients in the lower stratosphere from aerosol heating and ozone depletion, and to reduced tropospheric storminess. Interactions of the Quasi-Biennial Oscillation of tropical stratospheric winds with the climate system are also now better understood by examining their role in the Pinatubo response. We have more confidence in the sensitivity of climate models used for attribution and projection of anthropogenic effects on climate because the strength of the water vapor feedback has been validated with Pinatubo simulations. The response of the biosphere to the Pinatubo eruption also illustrates its sensitivity to climate change and clarifies portions of the carbon cycle. Death of coral in the Red Sea in the winter of 1991–1992 and an unusually large number of polar bear cubs born in the summer of 1992 were two responses to the characteristic winter and summer temperature responses of the climate system. This strengthens our concern about negative impacts of global warming on polar bears and other wildlife. Enhanced vegetation growth from more diffuse and less direct solar radiation took more carbon dioxide out of the atmosphere than normal, temporarily reducing the observed long-term increase in carbon dioxide. Continued research on the Pinatubo eruption and its aftermath will undoubtedly enhance our understanding of the climate system.