Triggering of volcanic activity by large earthquakes

Statistical analyses of temporal relationships between large earthquakes and volcanic eruptions suggest seismic waves may trigger eruptions even over great (>1000 km) distances, although the causative mechanism is not well constrained. In this study the relationship between large earthquakes and subtle changes in volcanic activity was investigated in order to gain greater insight into the relationship between dynamic stresses propagated by surface waves and volcanic response. Daily measurements from the Ozone Monitoring Instrument (OMI), onboard the Aura satellite, provide constraints on volcanic sulfur-dioxide (SO2) emission rates as a measure of subtle changes in activity. Time series of SO2 emission rates were produced from OMI data for thirteen persistently active volcanoes from 1 October 2004 to 30 September 2010. In order to quantify the affect of earthquakes at teleseismic distances, we modeled surface-wave amplitudes from the source mechanisms of moment magnitude (Mw) ≥7 earthquakes, and calculated the Peak Dynamic Stress (PDS). We assessed the influence of earthquakes on volcanic activity in two ways: 1) by identifying increases in the SO2 time series data and looking for causative earthquakes and 2) by examining the average emission rate before and after each earthquake. In the first, the SO2 time series for each volcano was used to calculate a baseline threshold for comparison with post-earthquake emission. Next, we generated a catalog of responses based on sustained SO2 emission increases above this baseline. Delay times between each SO2 response and each prior earthquake were analyzed using both the actual earthquake catalog, and a randomly generated catalog of earthquakes. This process was repeated for each volcano. Despite varying multiple parameters, this analysis did not demonstrate a clear relationship between earthquake-generated PDS and SO2 emission. However, the second analysis, which was based on the occurrence of large earthquakes indicated a response at most volcanoes. Using the PDS calculations as a filtering criterion for the earthquake catalog, the SO2 mass for each volcano was analyzed in 28-day windows centered on the earthquake origin time. If the average SO2 mass after the earthquake was greater than an arbitrary percentage of pre-earthquake mass, we identified the volcano as having a response to the event. This window analysis provided insight on what type of volcanic activity is more susceptible to triggering by dynamic stress. The volcanoes with very open systems included in this study, Ambrym, Gaua, Villarrica, Erta Ale and, Turrialba, showed a clear response to dynamic stress while the volcanoes with more closed systems, Merapi, Semeru, Fuego, Pacaya, and Bagana, showed no response.

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.

STATIC AND DYNAMIC STRESS CHANGE AT 27 VOLCANOES OF THE CENTRAL AMERICAN VOLCANIC ARC AFTER THE MW. 7.6 COSTA RICA EARTHQUAKE OF 5 SEPTEMBER 2012

Large earthquakes may strongly influence the activity of volcanoes through static and dynamic processes. In this study, we quantify the static and dynamic stress change on 27 volcanoes in Central America, after the Mw 7.6 Costa Rica earthquake of 5 September 2012. Following this event, 8 volcanoes showed signs of activity. We calculated the static stress change due to the earthquake on hypothetical faults under these volcanoes with Coulomb 3.3. For the dynamic stress change, we computed synthetic seismograms to simulate the waveforms at these volcanoes. We then calculated the Peak Dynamic Stress (PDS) from the modeled peak ground velocities. The resulting values are from moderate to minor changes in stress (10-1-10-2 MPa) with the PDS values generally an order of magnitude larger than the static stress change. Although these values are small, they may be enough to trigger a response by the volcanoes, and are on the order of stress changes implicated in many other studies of volcano and earthquake triggering by large earthquakes. This study provides insight into the poorly-constrained mechanism for remote triggering.

A Microscopic Study of some Volcanic Rocks in the Vicinity of Butte

The rhyolite rocks in the vicinity of Big Butte present a very interesting, though complicated study. It seems that no extensive or conclusive work has ever been done upon these extrusives. I have found that the rhyolites present a very interesting problem, especially microscopically, and that no doubt an intensive study along these lines will solve the problems concerning age relations of different flows and origin of the magma that produced the material for emanations.

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.

The Geology of Central Chile's Long-Lived Volcanic Arc

Discussion of subduction zones that are associated with volcanic arcs and major chemical exchanges between the Earth's surface and underlying mantle.

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.

Petrogenesis of rift-related tephrites, phonolites and trachytes (Central European Volcanic Province, Rhon, FRG): Constraints from Sr, Nd, Pb and O isotopes

The volcanic rocks of the Rhön area (Central European Volcanic Province, Germany) belong to a moderately alkali basaltic suite that is associated with minor tephriphonolites, phonotephrites, tephrites, phonolites and trachytes. Based on isotope sytematics (87Sr/86Sr: 0.7033–0.7042; 143Nd/144Nd: 0.51279–0.51287; 206Pb/204Pb: 19.1–19.5), the inferred parental magmas formed by variable degrees of partial melting of a common asthenospheric mantle source (EAR: European Asthenospheric Reservoir of Cebriá and Wilson, 1995). Tephrites, tephriphonolites, phonotephrites, phonolites and trachytes show depletions and enrichments in some trace elements (Sr, Ba, Nb, Zr, Y) indicating that they were generated by broadly similar differentiation processes that were dominated by fractionation of olivine, clinopyroxene, amphibole, apatite and titaniferous magnetite ± plagioclase ± alkalifeldspar. The fractionated samples seem to have evolved by two distinct processes. One is characterized by pure fractional crystallization indicated by increasing Nb (and other incompatible trace element) concentrations at virtually constant 143Nd/144Nd ~ 0.51280 and 87Sr/86Sr ~ 0.7035. The other process involved an assimilation–fractional crystallization (AFC) process where moderate assimilation to crystallization rates produced evolved magmas characterized by higher Nb concentrations at slightly lower 143Nd/144Nd down to 0.51275. Literature data for some of the evolved rocks show more variable 87Sr/86Sr ranging from 0.7037 to 0.7089 at constant 143Nd/144Nd ~ 0.51280. These features may result from assimilation of upper crustal rocks by highly differentiated low-Sr (< 100 ppm Sr) lavas. However, based on the displacement of the differentiated rocks from this study towards lower 143Nd/144Nd ratios and modeled AFC processes in 143Nd/144Nd vs. 87Sr/86Sr and 207Pb/204Pb vs. 143Nd/144Nd space assimilation of lower crustal rocks seems more likely. The view that assimilation of lower crustal rocks played a role is confirmed by high-precision double-spike Pb isotope data that reveal higher 207Pb/204Pb ratios (15.62–15.63) in the differentiated rocks than in the primitive basanites (15.58–15.61). This is compatible with incorporation of radiogenic Pb from lower crustal xenoliths (207Pb/204Pb: 15.63–15.69) into the melt. However, 206Pb/204Pb ratios are similar for the differentiated rocks (19.13–19.35) and the primitive basanites (19.12–19.55) implying that assimilation involved an ancient crustal end member with a higher U/Pb ratio than the mantle source of the basanites. In addition, alteration-corrected δ18O values of the differentiated rocks range from c. 5 to 7‰ which is the same range as observed in the primitive alkaline rocks. This study confirms previous interpretations that highlighted the role of AFC processes in the evolution of alkaline volcanic rocks in the Rhön area of the Central European Volcanic Province.

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.

The history of subaquatic volcanism recorded in the sediments of Lake Kivu

Subaquatic volcanic activity has been ongoing in Lake Kivu since the early Holocene and has a dynamic effect on the biological productivity in the surface water, and the preservation of carbonate in the deep anoxic water. Groundwater discharge into the lake’s deepwater propels the upward advection of the water column that ultimately supplies nutrients to the surface water for biological production. The amount of nutrients supplied from the deepwater can be increased suddenly by (1) a cold meteorological event that drives deep seasonal mixing resulting in increased nutrients from below and oxygen from above, or (2) subaquatic volcanic activity that induces a buoyant hydrothermal plume, which entrains nutrients from the deepwater and results in anoxia or suboxic conditions in the surface water. Previous sedimentological studies in Lake Kivu have hypothesized that regional climatic changes are responsible for sudden changes in the preservation of carbonates in the Main Basin. Here we reveal that sublacustrine volcanic events most likely induce the abrupt changes to the geochemistry in the sediment in Lake Kivu. An unprecedented look into the sediment stratigraphy and geochemistry from high-resolution seismic-reflection, and 15N-isotope analyses was conducted in the Main Basin. The results reveal that buoyant hydrothermal plumes caused by subaquatic volcanic activity are a possible trigger for increased biological productivity and organic matter preservation, and that ongoing hydrothermal activity increases the alkalinity in the deepwater, leading to carbonate preservation. The onset of carbonate preservation since the 1970s that is currently observed in the sediment could indicate that hydrothermal discharge has recently increased in the lake.