Atmospheric CO₂ response to volcanic eruptions: The role of ENSO, season, and variability

Tropical explosive volcanism is one of the most important natural factors that significantly impact the climate system and the carbon cycle on annual to multi-decadal time scales. The three largest explosive eruptions in the last 50�years�Agung, El Chichón, and Pinatubo�occurred in spring/summer in conjunction with El Niño events and left distinct negative signals in the observational temperature and CO2 records. However, confounding factors such as seasonal variability and El Niño-Southern Oscillation (ENSO) may obscure the forcing-response relationship. We determine for the first time the extent to which initial conditions, i.e., season and phase of the ENSO, and internal variability influence the coupled climate and carbon cycle response to volcanic forcing and how this affects estimates of the terrestrial and oceanic carbon sinks. Ensemble simulations with the Earth System Model (Climate System Model 1.4-carbon) predict that the atmospheric CO2 response is �60 larger when a volcanic eruption occurs during El Niño and in winter than during La Niña conditions. Our simulations suggest that the Pinatubo eruption contributed 11�±�6 to the 25�Pg terrestrial carbon sink inferred over the decade 1990�1999 and �2�±�1 to the 22�Pg oceanic carbon sink. In contrast to recent claims, trends in the airborne fraction of anthropogenic carbon cannot be detected when accounting for the decadal-scale influence of explosive volcanism and related uncertainties. Our results highlight the importance of considering the role of natural variability in the carbon cycle for interpretation of observations and for data-model intercomparison.

European summer temperature response to annually dated volcanic eruptions over the past nine centuries

The drop in temperature following large volcanic eruptions has been identified as an important component of natural climate variability. However, due to the limited number of large eruptions that occurred during the period of instrumental observations, the precise amplitude of post-volcanic cooling is not well constrained. Here we present new evidence on summer temperature cooling over Europe in years following volcanic eruptions. We compile and analyze an updated network of tree-ring maximum latewood density chronologies, spanning the past nine centuries, and compare cooling signatures in this network with exceptionally long instrumental station records and state-of-the-art general circulation models. Results indicate post-volcanic June–August cooling is strongest in Northern Europe 2 years after an eruption (−0.52 ± 0.05 °C), whereas in Central Europe the temperature response is smaller and occurs 1 year after an eruption (−0.18 ± 0.07 °C). We validate these estimates by comparison with the shorter instrumental network and evaluate the statistical significance of post-volcanic summer temperature cooling in the context of natural climate variability over the past nine centuries. Finding no significant post-volcanic temperature cooling lasting longer than 2 years, our results question the ability of large eruptions to initiate long-term temperature changes through feedback mechanisms in the climate system. We discuss the implications of these findings with respect to the response seen in general circulation models and emphasize the importance of considering well-documented, annually dated eruptions when assessing the significance of volcanic forcing on continental-scale temperature variations.

Volcanic forcing for climate modeling: a new microphysics-based data set covering years 1600–present

As the understanding and representation of the impacts of volcanic eruptions on climate have improved in the last decades, uncertainties in the stratospheric aerosol forcing from large eruptions are now linked not only to visible optical depth estimates on a global scale but also to details on the size, latitude and altitude distributions of the stratospheric aerosols. Based on our understanding of these uncertainties, we propose a new model-based approach to generating a volcanic forcing for general circulation model (GCM) and chemistry–climate model (CCM) simulations. This new volcanic forcing, covering the 1600–present period, uses an aerosol microphysical model to provide a realistic, physically consistent treatment of the stratospheric sulfate aerosols. Twenty-six eruptions were modeled individually using the latest available ice cores aerosol mass estimates and historical data on the latitude and date of eruptions. The evolution of aerosol spatial and size distribution after the sulfur dioxide discharge are hence characterized for each volcanic eruption. Large variations are seen in hemispheric partitioning and size distributions in relation to location/date of eruptions and injected SO2 masses. Results for recent eruptions show reasonable agreement with observations. By providing these new estimates of spatial distributions of shortwave and long-wave radiative perturbations, this volcanic forcing may help to better constrain the climate model responses to volcanic eruptions in the 1600–present period. The final data set consists of 3-D values (with constant longitude) of spectrally resolved extinction coefficients, single scattering albedos and asymmetry factors calculated for different wavelength bands upon request. Surface area densities for heterogeneous chemistry are also provided.

Antarctic Volcanic Flux Ratios from Law Dome Ice Cores

Explosive volcanic eruptions can inject large quantities of sulphur dioxide into the stratosphere. The aerosols that result from oxidation of the sulphur dioxide can produce significant cooling of the troposphere by reflecting or absorbing solar radiation. It is possible to obtain an estimate of the relative stratospheric sulphur aerosol concentration produced by different volcanoes by comparing sulphuric acid fluxes determined by analysis of polar ice cores. Here, we use a non-sea-salt sulphate time series derived from three well-dated Law Dome ice cores to investigate sulphuric acid flux ratios for major eruptions over the period AD 1301-1995. We use additional data from other cores to investigate systematic spatial variability in the ratios. Only for the Kuwae eruption (Law Dome ice date AD 1459.5) was the H2SO4 flux larger than that deposited by Tambora (Law Dome ice date AD 1816.7).

Timing of Tertiary Extension in the Railroad Valley Pioche Transect, Nevada: Constraints from Ar-40/Ar-39 Ages of Volcanic Rocks

Time-space relations of extension and volcanism place critical constraints on models of Basin and Range extensional processes. This paper addresses such relations in a 130-km-wide transect in the eastern Great Basin, bounded on the east by the Ely Springs Range and on the west by the Grant and Quinn Canyon ranges. Stratigraphic and structural data, combined with 40Ar/39Ar isotopic ages of volcanic rocks, document a protracted but distinctly episodic extensional history. Field relations indicate four periods of faulting. Only one of these periods was synchronous with nearby volcanic activity, which implies that volcanism and faulting need not be associated closely in space and time. Based on published dates and the analyses reported here, the periods of extension were (1) prevolcanic (pre-32 Ma), (2) early synvolcanic (30 to 27 Ma), (3) immediately postvolcanic (about 16 to 14 Ma), and (4) Pliocene to Quaternary. The break between the second and third periods is distinct. The minimum gap between the first two periods is 2 Ma, but the separation may be much larger. Temporal separation of the last two periods is only suggested by the stratigraphic record and cannot be rigorously demonstrated with present data. The three younger periods of faulting apparently occurred across the entire transect. The oldest period is recognized only at the eastern end of the transect, but appears to correlate about 150 km northward along strike with extension in the Northern Snake Range-Kern Mountains area. Therefore the oldest period also is regional in extent, but affected a different area than that affected by younger periods. This relation suggests that distinct extensional structures and master detachment faults were active at different times. The correlation of deformation periods of a few million years duration across the Railroad Valley-Pioche transect suggests that the scale of active extensional domains in the Great Basin may be greater than 100 km across strike.

The 1452 or 1453 A.D. Kuwae Eruption Signal Derived from Multiple Ice Core Records: Greatest Volcanic Sulfate Event of the Past 700 Years

We combined 33 ice core records, 13 from the Northern Hemisphere and 20 from the Southern Hemisphere, to determine the timing and magnitude of the great Kuwae eruption in the mid-15th century. We extracted volcanic deposition signals by applying a high-pass loess filter to the time series and examining peaks that exceed twice the 31 year running median absolute deviation. By accounting for the dating uncertainties associated with each record, these ice core records together reveal a large volcanogenic acid deposition event during 1453 - 1457 A. D. The results suggest only one major stratospheric injection from the Kuwae eruption and confirm previous findings that the Kuwae eruption took place in late 1452 or early 1453, which may serve as a reference to evaluate and improve the dating of ice core records. The average total sulfate deposition from the Kuwae eruption was 93 kg SO4/km(2) in Antarctica and 25 kg SO4/km(2) in Greenland. The deposition in Greenland was probably underestimated since it was the average value of only two northern Greenland sites with very low accumulation rates. After taking the spatial variation into consideration, the average Kuwae deposition in Greenland was estimated to be 45 kg SO4/km(2). By applying the same technique to the other major eruptions of the past 700 years our result suggests that the Kuwae eruption was the largest stratospheric sulfate event of that period, probably surpassing the total sulfate deposition of the Tambora eruption of 1815, which produced 59 kg SO4/km(2) in Antarctica and 50 kg SO4/km(2) in Greenland.

Impact of solar versus volcanic activity variations on tropospheric temperatures and precipitation during the Dalton Minimum

The aim of this work is to elucidate the impact of changes in solar irradiance and energetic particles versus volcanic eruptions on tropospheric global climate during the Dalton Minimum (DM, AD 1780–1840). Separate variations in the (i) solar irradiance in the UV-C with wavelengths λ < 250 nm, (ii) irradiance at wavelengths λ > 250 nm, (iii) in energetic particle spectrum, and (iv) volcanic aerosol forcing were analyzed separately, and (v) in combination, by means of small ensemble calculations using a coupled atmosphere–ocean chemistry–climate model. Global and hemispheric mean surface temperatures show a significant dependence on solar irradiance at λ > 250 nm. Also, powerful volcanic eruptions in 1809, 1815, 1831 and 1835 significantly decreased global mean temperature by up to 0.5 K for 2–3 years after the eruption. However, while the volcanic effect is clearly discernible in the Southern Hemispheric mean temperature, it is less significant in the Northern Hemisphere, partly because the two largest volcanic eruptions occurred in the SH tropics and during seasons when the aerosols were mainly transported southward, partly because of the higher northern internal variability. In the simulation including all forcings, temperatures are in reasonable agreement with the tree ring-based temperature anomalies of the Northern Hemisphere. Interestingly, the model suggests that solar irradiance changes at λ < 250 nm and in energetic particle spectra have only an insignificant impact on the climate during the Dalton Minimum. This downscales the importance of top–down processes (stemming from changes at λ < 250 nm) relative to bottom–up processes (from λ > 250 nm). Reduction of irradiance at λ > 250 nm leads to a significant (up to 2%) decrease in the ocean heat content (OHC) between 0 and 300 m in depth, whereas the changes in irradiance at λ < 250 nm or in energetic particles have virtually no effect. Also, volcanic aerosol yields a very strong response, reducing the OHC of the upper ocean by up to 1.5%. In the simulation with all forcings, the OHC of the uppermost levels recovers after 8–15 years after volcanic eruption, while the solar signal and the different volcanic eruptions dominate the OHC changes in the deeper ocean and prevent its recovery during the DM. Finally, the simulations suggest that the volcanic eruptions during the DM had a significant impact on the precipitation patterns caused by a widening of the Hadley cell and a shift in the intertropical convergence zone.

Volcanic Influence on European Summer Precipitation through Monsoons: Possible Cause for “Years without Summer”*

Strong tropical volcanic eruptions have significant effects on global and regional temperatures. Their effects on precipitation, however, are less well understood. Analyzing hydroclimatic anomalies after 14 strong eruptions during the last 400 years in climate reconstructions and model simulations, a reduction of the Asian and African summer monsoons and an increase of south-central European summer precipitation in the year following the eruption was found. The simulations provide evidence for a dynamical link between these phenomena. The weaker monsoon circulations weaken the northern branch of the Hadley circulation, alter the atmospheric circulation over the Atlantic–European sector, and increase precipitation over Europe. This mechanism is able to explain, for instance, the wet summer in parts of Europe during the “year without a summer” of 1816, which up to now has not been explained. This study underlines the importance of atmospheric teleconnections between the tropics and midlatitudes to better understand the regional climate response to stratospheric volcanic aerosols.

Ice Core Evidence for a Second Volcanic Eruption Around 1809 in the Northern Hemisphere

A volcanic signal observed in ice cores from both polar regions six years prior to Tambora is attributed to an unknown tropical eruption in 1809. Recovery of dacitic tephra from the 1809 horizon in a Yukon ice core ( Eclipse) that is chemically distinct from andesitic 1809 tephra found in Antarctic ice cores indicates a second eruption in the Northern Hemisphere at this time. Together with the similar magnitude and timing of the 1809 volcanic signal in the Arctic and Antarctic, this could suggest a large tropical eruption produced the sulfate and Antarctic tephra and a minor Northern Hemisphere eruption produced the Eclipse tephra. Nonetheless, the possibility that there were coincidental eruptions of similar magnitude in both hemispheres, rather than a single tropical eruption, should not be discounted. Correctly attributing the source of the 1809 volcanic signal has important implications for modeling the magnitude and latitudinal distribution of volcanic radiative forcing.

Late Quaternary history of the Vakinankaratra volcanic field (central Madagascar): insights from luminescence dating of phreatomagmatic eruption deposits

The Quaternary Vakinankaratra volcanic field in the central Madagascar highlands consists of scoria cones, lava flows, tuff rings, and maars. These volcanic landforms are the result of processes triggered by intracontinental rifting and overlie Precambrian basement or Neogene volcanic rocks. Infrared-stimulated luminescence (IRSL) dating was applied to 13 samples taken from phreatomagmatic eruption deposits in the Antsirabe–Betafo region with the aim of constraining the chronology of the volcanic activity. Establishing such a chronology is important for evaluating volcanic hazards in this densely populated area. Stratigraphic correlations of eruption deposits and IRSL ages suggest at least five phreatomagmatic eruption events in Late Pleistocene times. In the Lake Andraikiba region, two such eruption layers can be clearly distinguished. The older one yields ages between 109 ± 15 and 90 ± 11 ka and is possibly related to an eruption at the Amboniloha volcanic complex to the north. The younger one gives ages between 58 ± 4 and 47 ± 7 ka and is clearly related to the phreatomagmatic eruption that formed Lake Andraikiba. IRSL ages of a similar eruption deposit directly overlying basement laterite in the vicinity of the Fizinana and Ampasamihaiky volcanic complexes yield coherent ages of 68 ± 7 and 65 ± 8 ka. These ages provide the upper age limit for the subsequently developed Iavoko, Antsifotra, and Fizinana scoria cones and their associated lava flows. Two phreatomagmatic deposits, identified near Lake Tritrivakely, yield the youngest IRSL ages in the region, with respective ages of 32 ± 3 and 19 ± 2 ka. The reported K-feldspar IRSL ages are the first recorded numerical ages of phreatomagmatic eruption deposits in Madagascar, and our results confirm the huge potential of this dating approach for reconstructing the volcanic activity of Late Pleistocene to Holocene volcanic provinces.

A tephra lattice for Greenland and a reconstruction of volcanic events spanning 25–45 ka b2k

Tephra layers preserved within the Greenland ice-cores are crucial for the independent synchronisation of these high-resolution records to other palaeoclimatic archives. Here we present a new and detailed tephrochronological framework for the time period 25,000 e 45,000 a b2k that brings together results from 4 deep Greenland ice-cores. In total, 99 tephra deposits, the majority of which are preserved as cryptotephra, are described from the NGRIP, NEEM, GRIP and DYE-3 records. The major element signatures of single glass shards within these deposits indicate that 93 are basaltic in composition all originating from Iceland. Specifically, 43 originate from Grimsv € otn, 20 are thought to be sourced from the Katla volcanic system and 17 show affinity to the Kverkfj € oll system. Robust geochemical characterisations, independent ages derived from the GICC05 ice-core chronology, and the stratigraphic positions of these deposits relative to the Dansgaard-Oeschger climate events represent a key framework that provides new information on the frequency and nature of volcanic events in the North Atlantic region between GS-3 and GI-12. Of particular importance are 19 tephra deposits that lie on the rapid climatic transitions that punctuate the last glacial period. This framework of well-constrained, time-synchronous tie-lines represents an important step towards the independent synchronisation of marine, terrestrial and ice-core records from the North Atlantic region, in order to assess the phasing of rapid climatic changes during the last glacial period.

Timing and climate forcing of volcanic eruptions for the past 2,500 years

Volcanic eruptions contribute to climate variability, but quantifying these contributions has been limited by inconsistencies in the timing of atmospheric volcanic aerosol loading determined from ice cores and subsequent cooling from climate proxies such as tree rings. Here we resolve these inconsistencies and show that large eruptions in the tropics and high latitudes were primary drivers of interannual-to-decadal temperature variability in the Northern Hemisphere during the past 2,500 years. Our results are based on new records of atmospheric aerosol loading developed from high-resolution, multi-parameter measurements from an array of Greenland and Antarctic ice cores as well as distinctive age markers to constrain chronologies. Overall, cooling was proportional to the magnitude of volcanic forcing and persisted for up to ten years after some of the largest eruptive episodes. Our revised timescale more firmly implicates volcanic eruptions as catalysts in the major sixth-century pandemics, famines, and socioeconomic disruptions in Eurasia and Mesoamerica while allowing multi-millennium quantification of climate response to volcanic forcing.

Time-lapse pressure tomography for characterizing CO2 plume evolution in a deep saline aquifer

A time-lapse pressure tomography inversion approach is applied to characterize the CO2 plume development in a virtual deep saline aquifer. Deep CO2 injection leads to flow properties of the mixed-phase, which vary depending on the CO2 saturation. Analogous to the crossed ray paths of a seismic tomographic experiment, pressure tomography creates streamline patterns by injecting brine prior to CO2 injection or by injecting small amounts of CO2 into the two-phase (brine and CO2) system at different depths. In a first step, the introduced pressure responses at observation locations are utilized for a computationally rapid and efficient eikonal equation based inversion to reconstruct the heterogeneity of the subsurface with diffusivity (D) tomograms. Information about the plume shape can be derived by comparing D-tomograms of the aquifer at different times. In a second step, the aquifer is subdivided into two zones of constant values of hydraulic conductivity (K) and specific storage (Ss) through a clustering approach. For the CO2 plume, mixed-phase K and Ss values are estimated by minimizing the difference between calculated and “true” pressure responses using a single-phase flow simulator to reduce the computing complexity. Finally, the estimated flow property is converted to gas saturation by a single-phase proxy, which represents an integrated value of the plume. This novel approach is tested first with a doublet well configuration, and it reveals a great potential of pressure tomography based concepts for characterizing and monitoring deep aquifers, as well as the evolution of a CO2 plume. Still, field-testing will be required for better assessing the applicability of this approach.