Ultraviolet digital imaging of volcanic plumes : implementation and application to magmatic processes at basaltic volcanoes

Magmatic volatiles play a crucial role in volcanism, from magma production at depth to generation of seismic phenomena to control of eruption style. Accordingly, many models of volcano dynamics rely heavily on behavior of such volatiles. Yet measurements of emission rates of volcanic gases have historically been limited, which has restricted model verification to processes on the order of days or longer. UV cameras are a recent advancement in the field of remote sensing of volcanic SO2 emissions. They offer enhanced temporal and spatial resolution over previous measurement techniques, but need development before they can be widely adopted and achieve the promise of integration with other geophysical datasets. Large datasets require a means by which to quickly and efficiently use imagery to calculate emission rates. We present a suite of programs designed to semi-automatically determine emission rates of SO2 from series of UV images. Extraction of high temporal resolution SO2 emission rates via this software facilitates comparison of gas data to geophysical data for the purposes of evaluating models of volcanic activity and has already proven useful at several volcanoes. Integrated UV camera and seismic measurements recorded in January 2009 at Fuego volcano, Guatemala, provide new insight into the system’s shallow conduit processes. High temporal resolution SO2 data reveal patterns of SO2 emission rate relative to explosions and seismic tremor that indicate tremor and degassing share a common source process. Progressive decreases in emission rate appear to represent inhibition of gas loss from magma as a result of rheological stiffening in the upper conduit. Measurements of emission rate from two closely-spaced vents, made possible by the high spatial resolution of the camera, help constrain this model. UV camera measurements at Kilauea volcano, Hawaii, in May of 2010 captured two occurrences of lava filling and draining within the summit vent. Accompanying high lava stands were diminished SO2 emission rates, decreased seismic and infrasonic tremor, minor deflation, and slowed lava lake surface velocity. Incorporation of UV camera data into the multi-parameter dataset gives credence to the likelihood of shallow gas accumulation as the cause of such events.

Observations of volcanic emissions using satellite remote sensing

Volcanoes pose a threat to the human population at regional and global scales and so efficient monitoring is essential in order to effectively manage and mitigate the risks that they pose. Volcano monitoring from space has been possible for over thirty years and now, more than ever, a suite of instruments exists with the capability to observe emissions of gas and ash from a unique perspective. The goal of this research is to demonstrate the use of a range of satellite-based sensors in order to detect and quantify volcanic sulphur dioxide, and to assess the relative performances of each sensor against one another. Such comparisons are important in order to standardise retrievals and permit better estimations of the global contribution of sulphur dioxide to the atmosphere from volcanoes for climate modelling. In this work, retrievals of volcanic sulphur dioxide from a number of instruments are compared, and the individual performances at quantifying emissions from large, explosive volcanic eruptions are assessed. Retrievals vary widely from sensor to sensor, and often the use of a number of sensors in synergy can provide the most complete picture, rather than just one instrument alone. Volcanic emissions have the ability to result significant economic loses by grounding aircraft due to the high risk associated with ash encountering aircraft. As sulphur dioxide is often easier to measure than ash, it is often used as a proxy. This work examines whether this is a reasonable assumption, using the Icelandic eruption in early 2010 as a case study. Results indicate that although the two species are for the most part collocated, separation can occur under some conditions, meaning that it is essential to accurately measure both species in order to provide effective hazard mitigation. Finally, the usefulness of satellite remote sensing in quantifying the passive degassing from Turrialba, Costa Rica is demonstrated. The increase in activity from 2005 – 2010 can be observed in satellite data prior to the phreatic phase of early 2010, and can therefore potentially provide a useful indication of changing activity at some volcanoes.

EXPLORATION OF THE MTSAT2 SATELLITE CAPABILITIES FOR REAL TIME DETECTION AND CHARACTERIZATION OF VOLCANIC EMISSIONS

In this report, we attempt to define the capabilities of the infrared satellite remote sensor, Multifunctional Transport Satellite-2 (MTSAT-2) (i.e. a geosynchronous instrument), in characterizing volcanic eruptive behavior in the highly active region of Indonesia. Sulfur dioxide data from NASA's Ozone Monitoring Instrument (OMI) (i.e. a polar orbiting instrument) are presented here for validation of the processes interpreted using the thermal infrared datasets. Data provided from two case studies are analyzed specifically for eruptive products producing large thermal anomalies (i.e. lava flows, lava domes, etc.), volcanic ash and SO2 clouds; three distinctly characteristic and abundant volcanic emissions. Two primary methods used for detection of heat signatures are used and compared in this report including, single-channel thermal radiance (4-µm) and the normalized thermal index (NTI) algorithm. For automated purposes, fixed thresholds must be determined for these methods. A base minimum detection limit (MDL) for single-channel thermal radiance of 2.30E+05 Wm- 2sr-1m-1 and -0.925 for NTI generate false alarm rates of 35.78% and 34.16%, respectively. A spatial comparison method, developed here specifically for use in Indonesia and used as a second parameter for detection, is implemented to address the high false alarm rate. For the single-channel thermal radiance method, the utilization of the spatial comparison method eliminated 100% of the false alarms while maintaining every true anomaly. The NTI algorithm showed similar results with only 2 false alarms remaining. No definitive difference is observed between the two thermal detection methods for automated use; however, the single-channel thermal radiance method coupled with the SO2 mass abundance data can be used to interpret volcanic processes including the identification of lava dome activity at Sinabung as well as the mechanism for the dome emplacement (i.e. endogenous or exogenous). Only one technique, the brightness temperature difference (BTD) method, is used for the detection of ash. Trends of ash area, water/ice area, and their respective concentrations yield interpretations of increased ice formation, aggregation, and sedimentation processes that only a high-temporal resolution instrument like the MTSAT-2 can analyze. A conceptual model of a secondary zone of aggregation occurring in the migrating Kelut ash cloud, which decreases the distal fine-ash component and hazards to flight paths, is presented in this report. Unfortunately, SO2 data was unable to definitively reinforce the concept of a secondary zone of aggregation due to the lack of a sufficient temporal resolution. However, a detailed study of the Kelut SO2 cloud is used to determine that there was no climatic impacts generated from this eruption due to the atmospheric residence times and e-folding rate of ~14 days for the SO2. This report applies the complementary assets offered by utilizing a high-temporal and a high-spatial resolution satellite, and it demonstrates that these two instruments can provide unparalleled observations of dynamic volcanic processes.

EVALUATION OF THE EVOLVING STRESS FIELD OF THE YELLOWSTONE VOLCANIC PLATEAU FROM 1988 TO 2010 FROM EARTHQUAKE FIRST MOTIONS INVERSION

Within the Yellowstone National Park, Wyoming, the silicic Yellowstone volcanic field is one of the most active volcanic systems all over the world. Although the last rhyolite eruption occurred around 70,000 years ago, Yellowstone is still believed to be volcanically active, due to high hydrothermal and seismic activity. The earthquake data used in this study cover the period of time between 1988 and 2010. Earthquake relocations and a set of 369 well-constrained, double-couple, focal mechanism solutions were computed. Events were grouped according to location and time to investigate trends in faulting. The majority of the events has oblique, normal-faulting solutions. The overall direction of extension throughout the 0.64 Ma Yellowstone caldera looks nearly ENE, consistently with the direction of alignments of volcanic vents within the caldera, but detailed study revealed spatial and temporal variations. Stress-field solutions for different areas and time periods were calculated from earthquake focal mechanism inversion. A well-resolved rotation of σ3 was found, from NNE-SSW near the Hebgen Lake fault zone, to ENE-WSW near Norris Junction. In particular, the σ3 direction changed throughout the years in the Norris Junction area, from being ENE-WSW, as calculated in the study by Waite and Smith (2004), to NNE-SSW, while the other σ3 directions are mostly unchanged over time. The Yellowstone caldera was subject to periods of net uplift and subsidence over the past century, explained in previous studies as caused by expanding or contracting sills, at different depths. Based on the models used to explain these deformation periods, we investigated the relationship between variability in aseismic deformation and seismic activity and faulting styles. Focal mechanisms and P and T axes were divided into temporal and depth intervals, in order to identify spatial or temporal trends in deformation. The presence of “chocolate tablet” structures, with composite dilational faults, was identified in many stages of the deformation history both in the Norris Geyser Basin area and inside the caldera. Strike-slip component movement was found in a depth interval below a contracting sill, indicating the movement of magma towards the caldera.

Impact of volcanic stratospheric aerosols on diurnal temperature range in Europe over the past 200 years: Observations versus model simulations

We analyze the impact of stratospheric volcanic aerosols on the diurnal temperature range (DTR) over Europe using long-term subdaily station records. We compare the results with a 28-member ensemble of European Centre/Hamburg version 5.4 (ECHAM5.4) general circulation model simulations. Eight stratospheric volcanic eruptions during the instrumental period are investigated. Seasonal all- and clear-sky DTR anomalies are compared with contemporary (approximately 20 year) reference periods. Clear sky is used to eliminate cloud effects and better estimate the signal from the direct radiative forcing of the volcanic aerosols. We do not find a consistent effect of stratospheric aerosols on all-sky DTR. For clear skies, we find average DTR anomalies of −0.08°C (−0.13°C) in the observations (in the model), with the largest effect in the second winter after the eruption. Although the clear-sky DTR anomalies from different stations, volcanic eruptions, and seasons show heterogeneous signals in terms of order of magnitude and sign, the significantly negative DTR anomalies (e.g., after the Tambora eruption) are qualitatively consistent with other studies. Referencing with clear-sky DTR anomalies to the radiative forcing from stratospheric volcanic eruptions, we find the resulting sensitivity to be of the same order of magnitude as previously published estimates for tropospheric aerosols during the so-called “global dimming” period (i.e., 1950s to 1980s). Analyzing cloud cover changes after volcanic eruptions reveals an increase in clear-sky days in both data sets. Quantifying the impact of stratospheric volcanic eruptions on clear-sky DTR over Europe provides valuable information for the study of the radiative effect of stratospheric aerosols and for geo-engineering purposes.

Distinct lake level lowstand in Lake Prespa (SE Europe) at the time of the 74 (75) ka Toba eruption

Abstract. The 74 (75) ka Toba eruption in Sumatra, Indonesia, is considered to be one of the largest volcanic events during the Quaternary. Tephra from the Toba eruption has been found in many terrestrial and marine sedimentary deposits, and acidity peaks related to the eruption have been used to synchronize ice core records from Greenland and Antarctica. Seismic profiles and sedimentological data from Lake Prespa on the Balkan Peninsula, SE Europe, indicate a lake level lowstand at 73.6±7.7 ka based on ESR dating of shells. Tephrostratigraphy, radiocarbon dating and tuning of the total organic carbon content with the NGRIP isotope record, corroborate that the lake level lowstand was a short-term event superimposed on the general cooling trend at the end of MIS 5, most likely at the onset of the Greenland Stadial GS- 20. Acknowledging that tectonic events or karst processes could have triggered this lake level lowstand, the chronological correspondence between the lowstand and the Toba eruption is intriguing. Therefore a Toba-driven short-term shift to aridity in the Balkan region, leading to lake level changes and triggering spatial expansion events in one of the lake’s most abundant benthic species, the carino mussel Dreissena presbensis, cannot be excluded.

Northern hemispheric winter warming pattern after tropical volcanic eruptions: Sensitivity to the ozone climatology

An important key for the understanding of the dynamic response to large tropical volcanic eruptions is the warming of the tropical lower stratosphere and the concomitant intensification of the polar vortices. Although this mechanism is reproduced by most general circulation models today, most models still fail in producing an appropriate winter warming pattern in the Northern Hemisphere. In this study ensemble sensitivity experiments were carried out with a coupled atmosphere-ocean model to assess the influence of different ozone climatologies on the atmospheric dynamics and in particular on the northern hemispheric winter warming. The ensemble experiments were perturbed by a single Tambora-like eruption. Larger meridional gradients in the lower stratospheric ozone favor the coupling of zonal wind anomalies between the stratosphere and the troposphere after the eruption. The associated sea level pressure, temperature, and precipitation patterns are more pronounced and the northern hemispheric winter warming is highly significant. Conversely, weaker meridional ozone gradients lead to a weaker response of the winter warming and the associated patterns. The differences in the number of stratosphere-troposphere coupling events between the ensembles experiments indicate a nonlinear response behavior of the dynamics with respect to the ozone and the volcanic forcing.

Direct linking of Greenland and Antarctic ice cores at the Toba eruption (74 ka BP)

The Toba eruption that occurred some 74 ka ago in Sumatra, Indonesia, is among the largest volcanic events on Earth over the last 2 million years. Tephra from this eruption has been spread over vast areas in Asia, where it constitutes a major time marker close to the Marine Isotope Stage 4/5 boundary. As yet, no tephra associated with Toba has been identified in Greenland or Antarctic ice cores. Based on new accurate dating of Toba tephra and on accurately dated European stalagmites, the Toba event is known to occur between the onsets of Greenland interstadials (GI) 19 and 20. Furthermore, the existing linking of Greenland and Antarctic ice cores by gas records and by the bipolar seesaw hypothesis suggests that the Antarctic counterpart is situated between Antarctic Isotope Maxima (AIM) 19 and 20. In this work we suggest a direct synchronization of Greenland (NGRIP) and Antarctic (EDML) ice cores at the Toba eruption based on matching of a pattern of bipolar volcanic spikes. Annual layer counting between volcanic spikes in both cores allows for a unique match. We first demonstrate this bipolar matching technique at the already synchronized Laschamp geomagnetic excursion (41 ka BP) before we apply it to the suggested Toba interval. The Toba synchronization pattern covers some 2000 yr in GI-20 and AIM-19/20 and includes nine acidity peaks that are recognized in both ice cores. The suggested bipolar Toba synchronization has decadal precision. It thus allows a determination of the exact phasing of inter-hemispheric climate in a time interval of poorly constrained ice core records, and it allows for a discussion of the climatic impact of the Toba eruption in a global perspective. The bipolar linking gives no support for a long-term global cooling caused by the Toba eruption as Antarctica experiences a major warming shortly after the event. Furthermore, our bipolar match provides a way to place palaeo-environmental records other than ice cores into a precise climatic context.

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.

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

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.

The impacts of volcanic aerosol on stratospheric ozone and the Northern Hemisphere polar vortex: separating radiative-dynamical changes from direct effects due to enhanced aerosol heterogeneous chemistry

After major volcanic eruptions the enhanced aerosol causes ozone changes due to greater heterogeneous chemistry on the particle surfaces (HET-AER) and from dynamical effects related to the radiative heating of the lower stratosphere (RAD-DYN). We carry out a series of experiments with an atmosphere–ocean–chemistry–climate model to assess how these two processes change stratospheric ozone and Northern Hemispheric (NH) polar vortex dynamics. Ensemble simulations are performed under present day and preindustrial conditions, and with aerosol forcings representative of different eruption strength, to investigate changes in the response behaviour. We show that the halogen component of the HET-AER effect dominates under present-day conditions with a global reduction of ozone (−21 DU for the strongest eruption) particularly at high latitudes, whereas the HET-AER effect increases stratospheric ozone due to N2O5 hydrolysis in a preindustrial atmosphere (maximum anomalies +4 DU). The halogen-induced ozone changes in the present-day atmosphere offset part of the strengthening of the NH polar vortex during mid-winter (reduction of up to −16 m s-1 in January) and slightly amplify the dynamical changes in the polar stratosphere in late winter (+11 m s-1 in March). The RAD-DYN mechanism leads to positive column ozone anomalies which are reduced in a present-day atmosphere by amplified polar ozone depletion (maximum anomalies +12 and +18 DU for present day and preindustrial, respectively). For preindustrial conditions, the ozone response is consequently dominated by RAD-DYN processes, while under present-day conditions, HET-AER effects dominate. The dynamical response of the stratosphere is dominated by the RAD-DYN mechanism showing an intensification of the NH polar vortex in winter (up to +10 m s-1 in January). Ozone changes due to the RAD-DYN mechanism slightly reduce the response of the polar vortex after the eruption under present-day conditions.

Major and trace element analyses of ODP Leg 120 volcanic rocks (Table 1)

Widespread Lower Cretaceous magmatism occurred along the Indian-Australian/Antarctic margins, and in the juvenile Indian Ocean, during the rifting of eastern Gondwana. The formation of this magmatic province probably began around 120-130 Ma with the eruption of basalts on the Naturaliste Plateau and at Bunbury, western Australia. On the northeast margin of India, activity began around 117 Ma with the Rajmahal continental basalts and associated lamprophyre intrusions. The formation of the Kerguelen Plateau in the Indian Ocean began no later than 114 Ma. Ultramafic lamprophyres (alnoites) were emplaced in the Prince Charles Mountains near the Antarctic continental margin at ~ 110 Ma. These events are considered to be related to a major mantle plume, the remnant of which is situated beneath the region of Kerguelen and Heard islands at the present day. Geochemical data are presented for each of these volcanic suites and are indicative of complex interactions between asthenosphere-derived magmas and the continental lithosphere. Kerguelen Plateau basalts have Sr and Nd isotopic compositions lying outside the field for Indian Ocean mid-ocean ridge basalts (MORB) but, with the exception of Site 738 at the southern end of the plateau, within the range of more recent hotspot basalts from Kerguelen and Heard Islands. However, a number of the plateau tholeiites are characterized by lower 206Pb/204Pb ratios than are basalts from Kerguelen Island, and many also have anomalously high La/Nb ratios. These features suggest that the source of the Kerguelen Plateau basalts suffered contamination by components derived from the Gondwana continental lithosphere. An extreme expression of this lithospheric signature is shown by a tholeiite from Site 738, suggesting that the southernmost part of the Kerguelen Plateau may be underlain by continental crust. The Rajmahal tholeiites mostly fall into two distinct geochemical groups. Some Group I tholeiites have Sr and Nd isotopic compositions and incompatible element abundances, similar to Kerguelen Plateau tholeiites from Sites 749 and 750, indicating that the Kerguelen-Heard mantle plume may have directly furnished Rajmahal volcanism. However, their elevated 207Pb/204Pb ratios indicate that these magmas did not totally escape contamination by continental lithosphere. In contrast to the Group I tholeiites, significant contamination is suggested for Group II Rajmahal tholeiites, on the basis of incompatible element abundances and isotopic compositions. The Naturaliste Plateau and the Bunbury Basalt samples show varying degrees of enrichment in incompatible elements over normal MORB. The Naturaliste Plateau samples (and Bunbury Basalt) have high La/Nb ratios, a feature not inconsistent with the notion that the plateau may consist of stretched continental lithosphere, near the ocean-continent divide.