Krakatoa lives: The effect of volcanic eruptions on ocean heat content and thermal expansion

A suite of climate model experiments indicates that 20th Century increases in ocean heat content and sea-level ( via thermal expansion) were substantially reduced by the 1883 eruption of Krakatoa. The volcanically-induced cooling of the ocean surface is subducted into deeper ocean layers, where it persists for decades. Temporary reductions in ocean heat content associated with the comparable eruptions of El Chichon ( 1982) and Pinatubo ( 1991) were much shorter lived because they occurred relative to a non-stationary background of large, anthropogenically-forced ocean warming. Our results suggest that inclusion of the effects of Krakatoa ( and perhaps even earlier eruptions) is important for reliable simulation of 20th century ocean heat uptake and thermal expansion. Inter-model differences in the oceanic thermal response to Krakatoa are large and arise from differences in external forcing, model physics, and experimental design. Systematic experimentation is required to quantify the relative importance of these factors. The next generation of historical forcing experiments may require more careful treatment of pre-industrial volcanic aerosol loadings.

The significance of volcanic eruption strength and frequency for climate

A simple physical model of the atmospheric effects of large explosive volcanic eruptions is developed. Using only one input parameter - the initial amount of sulphur dioxide injected into the stratosphere - the global-average stratospheric optical-depth perturbation and surface temperature response are modelled. The simplicity of this model avoids issues of incomplete data (applicable to more comprehensive models), making it a powerful and useful tool for atmospheric diagnostics of this climate forcing mechanism. It may also provide a computationally inexpensive and accurate way of introducing volcanic activity into larger climate models. The modelled surface temperature response for an initial sulphur-dioxide injection, coupled with emission-history statistics, is used to demonstrate that the most climatically significant volcanic eruptions are those of sufficient explosivity to just reach into the stratosphere (and achieve longevity). This study also highlights the fact that this measure of significance is highly sensitive to the representation of the climatic response and the frequency data used, and that we are far from producing a definitive history of explosive volcanism for at least the past 1000 years. Given this high degree of uncertainty, these results suggest that eruptions that release around and above 0.1 Mt SO2 into the stratosphere have the maximum climatic impact.

The climatic effects of the direct injection of water vapour into the stratosphere by large volcanic eruptions

We describe a novel mechanism that can significantly lower the amplitude of the climatic response to certain large volcanic eruptions and examine its impact with a coupled ocean-atmosphere climate model. If sufficiently large amounts of water vapour enter the stratosphere, a climatically significant amount of water vapour can be left over in the lower stratosphere after the eruption, even after sulphate aerosol formation. This excess stratospheric humidity warms the tropospheric climate, and acts to balance the climatic cooling induced by the volcanic aerosol, especially because the humidity anomaly lasts for a period that is longer than the residence time of aerosol in the stratosphere. In particular, northern hemisphere high latitude cooling is reduced in magnitude. We discuss this mechanism in the context of the discrepancy between the observed and modelled cooling following the Krakatau eruption in 1883. We hypothesize that moist coignimbrite plumes caused by pyroclastic flows travelling over ocean rather than land, resulting from an eruption close enough to the ocean, might provide the additional source of stratospheric water vapour.

Long-term effect of volcanic forcing on ocean heat content

Explosive volcanic eruptions cause episodic negative radiative forcing of the climate system. Using coupled atmosphere-ocean general circulation models (AOGCMs) subjected to historical forcing since the late nineteenth century, previous authors have shown that each large volcanic eruption is associated with a sudden drop in ocean heat content and sea-level from which the subsequent recovery is slow. Here we show that this effect may be an artefact of experimental design, caused by the AOGCMs not having been spun up to a steady state with volcanic forcing before the historical integrations begin. Because volcanic forcing has a long-term negative average, a cooling tendency is thus imposed on the ocean in the historical simulation. We recommend that an extra experiment be carried out in parallel to the historical simulation, with constant time-mean historical volcanic forcing, in order to correct for this effect and avoid misinterpretation of ocean heat content changes

Tephra fallout hazard assessment for a Plinian eruption scenario at Volcán de Colima (Mexico)

Volcanic ash fallout associated with renewal of explosive activity at Colima, represents a serious threat to the surrounding urbanized area. Here we assess the tephra fallout hazard associated with a Plinian eruption scenario. The eruptive history of Volcán de Colima shows that Plinian eruptions occur approximately every 100 years and the last eruption, the 1913, represents the largest historic eruption of this volcano. We used the last eruption as a reference to discuss volcanic hazard and risk scenarios connected with ash fallout. Tephra fallout deposits are modeled using HAZMAP, a model based on a semi-analytical solution of the advection– diffusion–sedimentation equation for volcanic particles. Based on a statistical study of wind profiles at Colima region, we first reconstructed ash loading maps and then computed ground load probability maps for different seasons. The obtained results show that a Plinian eruptive scenario at Volcán de Colima, could seriously damage more than 10 small towns and ranches, and potentially affect big cities located at tens of kilometers from the eruptive center. The probability maps obtained are aimed to give support to the risk mitigation strategies

Ultra-distal tephra deposits from super-eruptions: examples from Toba, Indonesia and Taupo Volcanic Zone, New Zealand

Voluminous rhyolitic eruptions from Toba, Indonesia, and Taupo Volcanic Zone (TVZ), New Zealand, have dispersed volcanic ash over vast areas in the late Quaternary. The ~74 ka Youngest Toba Tuff (YTT) eruption deposited ash over the Bay of Bengal and the Indian subcontinent to the west. The ~340 ka Whakamaru eruption (TVZ) deposited the widespread Rangitawa Tephra, dominantly to the southeast (in addition to occurrences northwest of vent), extending across the landmass of New Zealand, and the South Pacific Ocean and Tasman Sea, with distal terrestrial exposures on the Chatham Islands. These super-eruptions involved ~2500 km^3 and ~1500 km3 of magma (dense-rock equivalent; DRE), respectively. Ultra-distal terrestrial exposures of YTT at two localities in India, Middle Son Valley, Madhya Pradesh, and Jurreru River Valley, Andhra Pradesh, at distances of >2000 km from the source caldera, show a basal ‘primary’ ashfall unit ~4 cm thick, although deposits containing reworked ash are up to ~3 m in total thickness. Exposures of Rangitawa Tephra on the Chatham Islands, >900 km from the source caldera, are ~15-30 cm thick. At more proximal localities (~200 km from source), Rangitawa Tephra is ~55-70 cm thick and characterized by a crystal-rich basal layer and normal grading. Both distal tephra deposits are characterized by very-fine ash (with high PM10 fractions) and are crystal-poor. Glass chemistry, stratigraphy and grain-size data for these distal tephra deposits are presented with comparisons of their correlation, dispersal and preservation. Using field observations, ash transport and deposition were modeled for both eruptions using a semi-analytical model (HAZMAP), with assumptions concerning average wind direction and strength during eruption, column shape and vent size. Model outputs provide new insights into eruption dynamics and better estimates of eruption volumes associ- ated with tephra fallout. Modeling based on observed YTT distal tephra thicknesses indicate a relatively low (<40 km high), very turbulent eruption column, consistent with deposition from a co-ignimbrite cloud extending over a broad region. Similarly, the Whakamaru eruption was modeled as producing a predominantly Plinian column (~45 km high), with dispersal to the southeast by strong prevailing winds. Significant ash fallout of the main dispersal direction, to the northwest of source, cannot be replicated in this modeling. The widespread dispersal of large volumes of fine ash from both eruptions may have had global environmental consequences, acutely affecting areas up to thousands of kilometers from vent.

Rossby wave-breaking analysis of explosive cyclones in the Euro-Atlantic sector

The two-way relationship between Rossby Wave-Breaking (RWB) and intensification of extra tropical cyclones is analysed over the Euro-Atlantic sector. In particular, the timing, intensity and location of cyclone development are related to RWB occurrences. For this purpose, two potential-temperature based indices are used to detect and classify anticyclonic and cyclonic RWB episodes from ERA-40 Re-Analysis data. Results show that explosive cyclogenesis over the North Atlantic (NA) is fostered by enhanced occurrence of RWB on days prior to the cyclone’s maximum intensification. Under such conditions, the eddy-driven jet stream is accelerated over the NA, thus enhancing conditions for cyclogenesis. For explosive cyclogenesis over the eastern NA, enhanced cyclonic RWB over eastern Greenland and anticyclonic RWB over the sub-tropical NA are observed. Typically only one of these is present in any given case, with the RWB over eastern Greenland being more frequent than its southern counterpart. This leads to an intensification of the jet over the eastern NA and enhanced probability of windstorms reaching Western Europe. Explosive cyclones evolving under simultaneous RWB on both sides of the jet feature a higher mean intensity and deepening rates than cyclones preceded by a single RWB event. Explosive developments over the western NA are typically linked to a single area of enhanced cyclonic RWB over western Greenland. Here, the eddy-driven jet is accelerated over the western NA. Enhanced occurrence of cyclonic RWB over southern Greenland and anticyclonic RWB over Europe is also observed after explosive cyclogenesis, potentially leading to the onset of Scandinavian Blocking. However, only very intense developments have a considerable influence on the large-scale atmospheric flow. Non-explosive cyclones depict no sign of enhanced RWB over the whole NA area. We conclude that the links between RWB and cyclogenesis over the Euro-Atlantic sector are sensitive to the cyclone’s maximum intensity, deepening rate and location.

Diagnosing the influence of diabatic processes on the explosive deepening of extratropical cyclones

A novel version of the classical surface pressure tendency equation (PTE) is applied to ERA-Interim reanalysis data to quantitatively assess the contribution of diabatic processes to the deepening of extratropical cyclones relative to effects of temperature advection and vertical motions. The five cyclone cases selected, Lothar and Martin in December 1999, Kyrill in January 2007, Klaus in January 2009, and Xynthia in February 2010, all showed explosive deepening and brought considerable damage to parts of Europe. For Xynthia, Klaus and Lothar diabatic processes contribute more to the observed surface pressure fall than horizontal temperature advection during their respective explosive deepening phases, while Kyrill and Martin appear to be more baroclinically driven storms. The powerful new diagnostic tool presented here can easily be applied to large numbers of cyclones and will help to better understand the role of diabatic processes in future changes in extratropical storminess.

Explosive development of winter storm Xynthia over the subtropical North Atlantic Ocean

In winter of 2009–2010 south-western Europe was hit by several destructive windstorms. The most important was Xynthia (26–28 February 2010), which caused 64 reported casualties and was classified as the 2nd most expensive natural hazard event for 2010 in terms of economic losses. In this work we assess the synoptic evolution, dynamical characteristics and the main impacts of storm Xynthia, whose genesis, development and path were very uncommon. Wind speed gusts observed at more than 500 stations across Europe are evaluated as well as the wind gust field obtained with a regional climate model simulation for the entire North Atlantic and European area. Storm Xynthia was first identified on 25 February around 30° N, 50° W over the subtropical North Atlantic Ocean. Its genesis occurred on a region characterized by warm and moist air under the influence of a strong upper level wave embedded in the westerlies. Xynthia followed an unusual SW–NE path towards Iberia, France and central Europe. The role of moist air masses on the explosive development of Xynthia is analysed by considering the evaporative sources. A lagrangian model is used to identify the moisture sources, sinks and moisture transport associated with the cyclone during its development phase. The main supply of moisture is located over an elongated region of the subtropical North Atlantic Ocean with anomalously high SST, confirming that the explosive development of storm Xynthia had a significant contribution from the subtropics.

The role of anomalous SST and surface fluxes over the Southeastern North Atlantic in the explosive development of windstorm Xynthia

In late February 2010 the extraordinary windstorm Xynthia crossed over Southwestern and Central Europe and caused severe damage, affecting particularly the Spanish and French Atlantic coasts. The storm was embedded in uncommon large-scale atmospheric and boundary conditions prior to and during its development, namely enhanced sea surface temperatures (SST) within the low-level entrainment zone of air masses, an unusual southerly position of the polar jet stream, and a remarkable split jet structure in the upper troposphere. To analyse the processes that led to the rapid intensification of this exceptional storm originating close to the subtropics (30°N), the sensitivity of the cyclone intensification to latent heat release is determined using the regional climate model COSMO-CLM forced with ERA-Interim data. A control simulation with observed SST shows that moist and warm air masses originating from the subtropical North Atlantic were involved in the cyclogenesis process and led to the formation of a vertical tower with high values of potential vorticity (PV). Sensitivity studies with reduced SST or increased laminar boundary roughness for heat led to reduced surface latent heat fluxes. This induced both a weaker and partly retarded development of the cyclone and a weakening of the PV-tower together with reduced diabatic heating rates, particularly at lower and mid levels. We infer that diabatic processes played a crucial role during the phase of rapid deepening of Xynthia and thus to its intensity over the Southeastern North Atlantic. We suggest that windstorms like Xynthia may occur more frequently under future climate conditions due to the warming SSTs and potentially enhanced latent heat release, thus increasing the windstorm risk for Southwestern Europe.

Fireballs, flares, and flickering: A semianalytic, LTE, explosive model from accretion disks to supernovae

We derive simple analytic expressions for the continuum light curves and spectra of flaring and flickering events that occur over a wide range of astrophysical systems. We compare these results to data taken from the cataclysmic variable SS Cygni and also from SN 1987A, deriving physical parameters for the material involved. Fits to the data indicate a nearly time-independent photospheric temperature arising from the strong temperature dependence of opacity when hydrogen is partially ionized.

On the relation between explosive cyclones affecting Europe and the North Atlantic oscillation

Intense winter cyclones often lead to hazardous weather over Europe. Previous studies have pointed to a link between the North Atlantic Oscillation (NAO) and strong European windstorms. However, the robustness of this relation for cyclones of varying intensities remains largely unexplored. In this paper, the bi-directional relation between the NAO and cyclones impacting Europe is analyzed for the period 1950–2010 focusing on the sensitivity to storm intensity. Evidence is given that explosive (EC) and non-explosive cyclones (NoEC) predominantly develop under different large-scale circulation conditions over the North Atlantic. Whereas NoEC evolve more frequently under negative and neutral NAO phases, the number of EC is larger under a positive NAO phase, typically characterized by an intensified jet toward Western Europe. Important differences are also found on the dynamics of NAO evolution after peak intensity for both cyclone populations.

Climate models without pre-industrial volcanic forcing underestimate historical ocean thermal expansion

Episodic explosive volcanic eruptions are a natural part of the climate system but are often omitted from atmosphere-ocean general circulation model (AOGCM) preindustrial spin-up and control experiments. This omission imposes a negative bias on ocean heat uptake in simulations of the historical period. In models of a range of complexity, we find that global-mean sea level rise due to thermal expansion during the last ∼ 150 years is consequently underestimated by 5–30 mm, which is a substantial proportion of the model mean of 50 mm in Coupled Model Intercomparison Project Phase 3 AOGCMs with anthropogenic forcing only, and is therefore important in accounting for 20th century sea level rise. We test and recommend a procedure for removing the bias.

Counteracting the climate effects of volcanic eruptions using short-lived greenhouse gases

A large volcanic eruption might constitute a climate emergency, significantly altering global temperature and precipitation for several years. Major future eruptions will occur, but their size or timing cannot be predicted. We show, for the first time, that it may be possible to counteract these climate effects through deliberate emissions of short-lived greenhouse gases, dampening the abrupt impact of an eruption. We estimate an emission pathway countering a hypothetical eruption 3 times the size of Mount Pinatubo in 1991. We use a global climate model to evaluate global and regional responses to the eruption, with and without counteremissions. We then raise practical, financial, and ethical questions related to such a strategy. Unlike the more commonly discussed geoengineering to mitigate warming from long-lived greenhouse gases, designed emissions to counter temporary cooling would not have the disadvantage of needing to be sustained over long periods. Nevertheless, implementation would still face significant challenges.