Climate response to the Toba super-eruption: Regional changes
Introduction
The Young Toba Tuff (YTT), eruption in northern Sumatra is considered to represent the most recent “super eruption” on Earth based on the amount of solid ejecta making this eruption a magnitude 81 event (2500–3000 km3, Rose and Chesner, 1987, Chesner et al., 1991). The timing of the event near the stadial climate prevailing about 73 to 72 ka (Ninkovich et al., 1978, Chesner et al., 1991) initially led to the suggestion that the massive eruption accelerated the long-term cooling trend to the stadial (Rampino and Self, 1992). Zielinski et al. (1996) reported that cooling already had begun before the sulphate peak in the Greenland GISP2 ice core at 71 ± 5 ka, which they attributed to the YTT eruption. Higher resolution oxygen isotope stratigraphy by Stuiver and Grootes (2000) places the sulphate peak near the end of the cooling phase. Schulz et al. (2002) dated the YTT ash layer in their high-resolution benthic and planktonic δ18O records slightly younger than the Marine Isotope Stage 5–4 boundary, leading to a YTT eruption age of the Toba volcano of about 70 ka, matching Greenland (GISP2 and GRIP) ice-core chronologies. The exact timing of the YTT eruption suffers from large uncertainties in ice core chronologies and other proxy data for the era. No traces of YTT ash have been found in ice cores that would allow unambiguous identification of the origin of isotope peaks. Sediment cores also have small sampling rates and may suffer from reworking. Hence, at best, it can currently be stated that the YTT eruption occurred at some time during the cooling towards the Marine Isotope Stage 5-4 boundary. The exact timing of the eruption is important, as different responses of the climate system would be expected, especially of vegetation cover, to the strong volcanic forcing under different climate background conditions.
Besides the timing problems, there are also other uncertainties of the effects of the YTT eruption. Although the total amount of ejected ash and its sediment pattern is reasonably well documented (e.g. Rose and Chesner, 1987, Chesner et al., 1991, Acharyya and Basu, 1993, Westgate et al., 1998, Jones, 2010), there are considerable uncertainties regarding the sulphur yield of the eruption, which is the parameter mainly determining the global climatic consequences of the eruption. Until the publication of sulphate concentrations in the Antarctic EPICA-Dome C ice core (Traversi et al., 2009), there was no evidence of any sulphate peak around the time of the YTT eruption over Antarctica. Although the data of Traversi et al. (2009) show sulphate concentrations in reasonable agreement with the GISP2 core (∼1200 ppb SO42+ vs. ∼2000 ppb SO42+), uncertainty in the dating accuracy (ranging between 1.5 ka and 3 ka), next to the missing direct identification of the ash, leaves it an open question whether or not both spikes result from the same source. If they were indeed a result of the YTT eruption, the estimated sulphur emission would be 60–100 times (510–850 Mt S) that from Mount Pinatubo (June 1991). This estimate is in contrast with results from Scaillet et al. (1998) and Chesner and Luhr (2010) suggesting a sulphur yield of not more than 2–5 times the one from Pinatubo for the YTT eruption, i.e. similar to the 1815 eruption of Tambora. Hence, the main parameter determining the climatic consequences of the YTT eruption is highly uncertain. Its estimates vary by more than one order of magnitude. From the publications currently available, the commonly assumed stratospheric sulphur injection of 100 times Pinatubo is certainly an upper limit. Recent simulations with coupled atmosphere-ocean models (Jones et al., 2005, Robock et al., 2009) yielded a decade of severe cooling by up to 10 K (global mean) for a YTT eruption comparable to 100 times Pinatubo: however, there was no further sustained cooling. Earth System Model (ESM) simulations (Timmreck et al., 2010), also using 100 times Pinatubo sulphur emission but, for the first time, also considering microphysical processes as well, indicate that global mean temperature anomalies are three times smaller than suggested by Jones et al. (2005) and Robock et al. (2009). The reduced response mainly arises from a more complete treatment of stratospheric aerosol formation and growth, leading to much weaker radiative forcing due to larger particle sizes and a faster removal rate.
Climatic consequences of the YTT eruption are a crucial argument in the current discussion about the fate of modern humans especially in India and for African refugees around 70 ka (e.g. Ambrose, 1998, Williams et al., 2009, Haslam and Petraglia, 2010). Of particular interest is temperature, precipitation, runoff, and food supply (e.g. vegetation) changes in potential regions of human refugees, which directly impact the evolutionary aspects (see Fig. 1 for an overview of relevant YTT eruption climate effects). Following Timmreck et al. (2010) who have already analyzed the simulated temperature and precipitation anomalies after the YTT eruption over the Indian subcontinent, the focus of ESM analysis in this study is extended to more regions relevant to human evolutionary issues during that time (see Fig. 2): most importantly, Southern and Eastern Africa (Scholz et al., 2007, Jacobs et al., 2008, Drake et al., 2011 and others), but also South East Asia where evidence was found of a high survival rate of mammalian mega fauna during the time of the YTT eruption (Louys, 2007, Louys, 2012 this volume). The sensitivity of the ESM model results are tested with respect to the assumed sulphur amount of a Tambora-like eruption which is considered to be the lower limit for the YTT eruption. Robock et al. (2009) already investigated the uncertainty in the sulphur emission for the YTT eruption, considering values in the range of 33–900 times the Pinatubo injection, and obtained global maximum cooling between 8 and 17 K. However they did not take into account the temporal evolution of the aerosol size distribution, which turned out to be an important diminishing factor for the volcanic climate signal (Timmreck et al., 2010). Furthermore, the hydrological cycle is investigated in detail by analyzing large river catchments within or close to the regions of interest: Ganges/Brahmaputra for the northern and Krishna for the southern part of India, Orange for Southern Africa, Nile for North East Africa, and Mekong for South East Asia. Vegetation plays the dominant role in the food chain, and therefore the climatic disturbances simulated by the ESM are used as drivers of a dynamical global vegetation model (DGVM) to estimate the lasting effects of the transient climate anomalies on vegetation. In order to account for the uncertainties in the timing of the YTT eruption, the DGVM was forced with climate anomalies under glacial and interglacial background conditions. This allows for an estimate of the upper and lower bounds of vegetation anomalies.
Section snippets
Model description
For the simulation of the climatic impacts of the YTT eruption, the two-step modelling approach described in Timmreck et al. (2010) is applied. In the first step, a simulation of the volcanic aerosol forcing is performed with the global aerosol atmospheric model MAECHAM5/HAM (Niemeier et al., 2009). In the second step, an ensemble of climate simulations using the forcing parameters derived from step one is performed with the ESM developed at the Max Planck Institute for Meteorology (MPI-ESM) (
Global and regional cooling
As demonstrated in Timmreck et al. (2010), the consideration of microphysical processes leads to strongly reduced radiative flux anomalies with maximum values of −18 W/m2 at the top of the atmosphere and −49 W/m2 at the surface lasting over four years. The temperature response in the MPI-ESM simulations of the YTT eruption is therefore shorter and weaker than estimated by previous model studies (Jones et al., 2005, Robock et al., 2009). Ensemble-mean global temperature anomaly patterns averaged
Uncertainty in the input parameters
To estimate the climatic effects of very large volcanic eruptions, the stratospheric sulphur (S) emission of the erupting magma is a key aspect, which poses an important uncertainty in the model simulations. This analysis has assumed a sulphur injection rate of 850 Mt S, which is at the upper end of recent estimates. According to the study by Chesner and Luhr (2010), the amount of sulphate mass released by the YTT eruption is more similar to the Tambora eruption. Thus, the sensitivity of the
Summary and conclusion
This paper has presented results of an ensemble simulation of the YTT eruption of the Toba volcano. The main simulation is based on a stratospheric sulphur emission in the order of 100 times that of the Mt. Pinatubo (1991) eruption and under interglacial background climate conditions, which also leads to a maximum estimate of the resulting climate anomalies.
Global and hemispheric scale maximum cooling is less severe than in previous simulations, as the model provides a more realistic
Acknowledgements
This work benefits from stimulating discussions within the MPI-M Super Volcano and Millennium projects D.Z. acknowledges funding by the ENIGMA project from the Max Planck Society. This publication is contribution no. 233 of the Sonderforschungsbereich 574 “Volatiles and Fluids in Subduction Zones” at Kiel University. Computations were done at the German Climate Computer Center (DKRZ).
References (61)
- et al.
Toba ash on the Indian subcontinent and its implications for correlation of Late Pleistocene alluvium
Quaternary Research
(1993) Late Pleistocene human population bottlenecks, volcanic winter, and differentiation of modern humans
Journal of Human Evolution
(1998)- et al.
Melt inclusion study of the Toba Tuffs, Sumatra, Indonesia
Journal of Volcanology and Geothermal Research
(2010) - et al.
Comment on "Environmental impact of the 73 ka Toba super-eruption in South Asia
Palaeogeography, Palaeoclimatology, Palaeoecology
(2010) Palaeoenvironmental response to the ∼74 ka Toba ash-fall in the Jurreru and Middle Son valleys in southern and north-central India
Quaternary Research
(2010)Limited effect of the Quaternary’s largest super-eruption (Toba) on land mammals from Southeast Asia
Quaternary Science Reviews
(2007)- et al.
The Max Planck Institute global ocean/sea ice model with orthogonal curvilinear coordinates
Ocean Modelling
(2003) - et al.
The Toba volcanic event and interstadial/stadial climates at the Marine Isotopic Stage 5 to 4 transition in the northern Indian Ocean
Quaternary Research
(2002) - et al.
GISP2 oxygen isotope Ratios
Quaternary Research
(2000) - et al.
All Toba tephra occurrences across Peninsular India belong to the 75,000 yr. B.P. Eruption
Quaternary Research
(1998)
Environmental impact of the 73 ka Toba super-eruption in South Asia
Palaeogeography, Palaeoclimatology, Palaeoecology
Modeling the "Year without summer" 1816 with the CCM SOCOL
Geophysical Research Abstracts
Impact of volcanic super-eruptions on the global carbon cycle
Geophysical Research Abstracts
Oxidation of volcanic SO2: a sink for stratospheric OH and H2O
Geophysical Research Letters
The role of microphysical and chemical processes in prolonging the climate forcing of the Toba Eruption
Geophysical Research Letters
Impact of an improved shortwave radiation scheme in the MAECHAM5 General Circulation Model
Atmospheric Chemistry and Physics
A. Eruptive history of Earth’s largest Quaternary caldera (Toba, Indonesia) clarified
Geology
Volcanism and the little ice age
PAGES Newsletter
Agricultural net primary production in relation to that liberated by the extinction of Pleistocene mega-herbivores: an estimate of agricultural carrying capacity?
Environmental Research Letters
Ancient watercourses and biogeography of the Sahara explain the peopling of the desert
Proceedings of the National Academy of Science
Computations of solar heating of the earth’s atmosphere: a new parameterization
Contributions to Atmospheric Physics
Climatology and forcing of the quasi-biennial oscillation in the MAECHAM5 model
Journal of Climate
Evaluation of the hydrological cycle in the ECHAM5 model
Journal of Climate
The Energetic Metabolism of Societies. Part II: empirical examples
Journal of Industrial Ecology
Applying the three-dimensional model ATHAM to volcanic plumes: dynamic of large co-ignimbrite eruptions and associated injection heights for volcanic gases
Geophysical Research Letters
A New Sea surface temperature and Sea ice boundary Dataset for the community atmosphere model
Journal of Climate
Ages for the middle Stone age of southern Africa: implications for human Behavior and dispersal
Science
An AOGCM simulation of the climate response to a volcanic super-eruption
Climate Dynamics
Climate and carbon-cycle variability over the last millennium
Climate of the Past
Mammal community structure of Sundanese fossil assemblages from the Late Pleistocene, and a discussion on the ecological effects of the Toba eruption
Quaternary International
Cited by (59)
Bipolar ice-core records constrain possible dates and global radiative forcing following the ∼74 ka Toba eruption
2023, Quaternary Science ReviewsA review study on the artifacts, vertebrate remains, and volcanic ash of peninsular Malaysia and their significance to reconstruct paleoclimate
2023, The Role of Tropics in Climate Change: Global Case StudiesHad Youngest Toba Tuff (YTT, ca. 75 ka) eruption really destroyed living media explicitly in entire Southeast Asia or just a theoretical debate? An extensive review of its catastrophic event
2022, Journal of Asian Earth Sciences: XCitation Excerpt :Based on the analytical work and climatic simulations, several studies made in the past predicted the concentration of considerable numbers of the humans, animals and vegetation in the southern, western and central parts of the India subcontinent long before the YTT eruption and also survived its primary ash fall event (Petraglia et al., 2007, 2012; Haslam et al., 2010; Timmreck et al., 2012; Clarkson et al., 2020). Using isotopic concentrations of carbonates associated with the YTT ash in Jurreru valley of southern India, Haslam et al. (2010) predicted that the cooling trend in the Indian Peninsular began considerably before the YTT eruption and warmer conditions introduced just before the said eruption (Timmreck et al., 2012). According to Petraglia et al. (2012), the humans and animals were present in the continent over 140 ka and had survived minimal cooling created by the YTT event.
Volcanic emissions, aerosol processes, and climatic effects
2022, Aerosols and ClimateUnderstanding the overestimated impact of the Toba volcanic super-eruption on global environments and ancient hominins
2020, Quaternary International