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.

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