Toward Forecasting Volcanic Eruptions using Seismic Noise

During inter-eruption periods, magma pressurization yields subtle changes of the elastic properties of volcanic edifices. We use the reproducibility properties of the ambient seismic noise recorded on the Piton de la Fournaise volcano to measure relative seismic velocity variations of less than 0.1 % with a temporal resolution of one day. Our results show that five studied volcanic eruptions were preceded by clearly detectable seismic velocity decreases within the zone of magma injection. These precursors reflect the edifice dilatation induced by magma pressurization and can be useful indicators to improve the forecasting of volcanic eruptions.Comment: Supplementary information: http://www-lgit.obs.ujf-grenoble.fr/~fbrengui/brenguier_SI.pdf Supplementary video: http://www-lgit.obs.ujf-grenoble.fr/~fbrengui/brenguierMovieVolcano.av

Source Mechanism of Small Long-Period Events at Mount St. Helens in July 2005 Using Template Matching, Phase-Weighted Stacking, and Full-Waveform Inversion

Long-period (LP, 0.5-5 Hz) seismicity, observed at volcanoes worldwide, is a recognized signature of unrest and eruption. Cyclic LP “drumbeating” was the characteristic seismicity accompanying the sustained dome-building phase of the 2004–2008 eruption of Mount St. Helens (MSH), WA. However, together with the LP drumbeating was a near-continuous, randomly occurring series of tiny LP seismic events (LP “subevents”), which may hold important additional information on the mechanism of seismogenesis at restless volcanoes. We employ template matching, phase-weighted stacking, and full-waveform inversion to image the source mechanism of one multiplet of these LP subevents at MSH in July 2005. The signal-to-noise ratios of the individual events are too low to produce reliable waveform inversion results, but the events are repetitive and can be stacked. We apply network-based template matching to 8 days of continuous velocity waveform data from 29 June to 7 July 2005 using a master event to detect 822 network triggers. We stack waveforms for 359 high-quality triggers at each station and component, using a combination of linear and phase-weighted stacking to produce clean stacks for use in waveform inversion. The derived source mechanism points to the volumetric oscillation (∼10 m3) of a subhorizontal crack located at shallow depth (∼30 m) in an area to the south of Crater Glacier in the southern portion of the breached MSH crater. A possible excitation mechanism is the sudden condensation of metastable steam from a shallow pressurized hydrothermal system as it encounters cool meteoric water in the outer parts of the edifice, perhaps supplied from snow melt

Rapid ascent of rhyolitic magma at Chaitén volcano, Chile

International audienceAlthough rhyolite magma has fuelled some of the Earth's largest explosive volcanic eruptions, our understanding of these events is incomplete due to the previous lack of direct observation of these eruptions. On 1 May 2008, Chaitén volcano in Chile erupted rhyolite magma unexpectedly and explosively. Here, petrological and experimental data are presented that indicate that the hydrous rhyolite magma at Chaitén ascended very rapidly from storage depth to near-surface, with velocities of the order of one metre per second

Earthquake source parameters and scaling relationships in Hungary (central Pannonian basin)

Abstract Fifty earthquakes that occurred in Hungary (central part of the Pannonian basin) with local magnitude ML ranging from 0.8 to 4.5 have been analyzed. The digital seismograms used in this study were recorded by six permanent broad-band stations and twenty short-period ones at hypocentral distances between 10 and 327 km. The displacement spectra for P- and SH-waves were analyzed according to Brune’s source model. Observed spectra were corrected for path-dependent attenuation effects using an independent regional estimate of the quality factor QS. To correct spectra for near-surface attenuation, the k parameterwas calculated, obtaining it fromwaveforms recorded at short epicentral distances. The values of the k parameter vary between 0.01 to 0.06 s with a mean of 0.03 s for P-waves and between 0.01 to 0.09 s with a mean of 0.04 s for SH-waves. After correction for attenuation effects, spectral parameters (corner frequency and low-frequency spectral level) were estimated by a grid search algorithm. The obtained seismic moments range from4.21×1011 to 3.41×1015 Nm (1.7≤Mw ≤4.3). The source radii are between 125 and 1343 m. Stress drop values vary between 0.14 and 32.4 bars with a logarithmic mean of 2.59 bars (1 bar = 105 Pa). From the results, a linear relationship between local andmomentmagnitudes has been established. The obtained scaling relations show slight evidence of self-similarity violation. However, due to the high scatter of our data, the existence of self-similarity cannot be excluded

Crystal and melt inclusion timescales reveal the evolution of magma migration before eruption

Volatile element concentrations measured in melt inclusions are a key tool used to understand magma migration and degassing, although their original values may be affected by different re-equilibration processes. Additionally, the inclusion-bearing crystals can have a wide range of origins and ages, further complicating the interpretation of magmatic processes. To clarify some of these issues, here we combined olivine diffusion chronometry and melt inclusion data from the 2008 eruption of Llaima volcano (Chile). We found that magma intrusion occurred about 4 years before the eruption at a minimum depth of approximately 8 km. Magma migration and reaction became shallower with time, and about 6 months before the eruption magma reached 3–4 km depth. This can be linked to reported seismicity and ash emissions. Although some ambiguities of interpretation still remain, crystal zoning and melt inclusion studies allow a more complete understanding of magma ascent, degassing, and volcano monitoring data.NRF (Natl Research Foundation, S’pore)MOE (Min. of Education, S’pore)Published versio

The propagation and seismicity of dyke injection, new experimental evidence

To reach the surface, dykes must overcome the inherent tensile strength of the country rock. As they do, they generate swarms of seismic signals, frequently used for forecasting. In this study we pressurize and inject molten acrylic into an encapsulating host rocks of (1) Etna basalt and (2) Comiso limestone, at 30 MPa of confining pressure. Fracture was achieved at 12 MPa for Etna basalt and 7.2 MPa for Comiso limestone. The generation of radial fractures was accompanied by acoustic emissions (AE) at a dominant frequency of 600 kHz. During “magma” movement in the dykes, AE events of approximately 150 kHz dominant frequency were recorded. We interpret our data using AE location and dominant frequency analysis, concluding that the seismicity associated with magma transport in dykes peaks during initial dyke creation but remains significant as long as magma movement continues. These results have important implications for seismic monitoring of active volcanoes

Dynamic Evolution of Microscopic Wet Cracking Noises

Characterizing the interaction between water and microscopic defects is one of the long-standing challenges in understanding a broad range of cracking processes. Different physical aspects of microscopic events, driven or influenced by water, have been extensively discussed in atomistic calculations but have not been accessible in microscale experiments. Through the analysis of the emitted noises during the evolution of individual, dynamic microcracking events, we show that the onset of a secondary instability known as hybrid events occurs during the fast healing phase of microcracking, which leads to (local) sudden increase of pore water pressure in the process zone, inducing a secondary instability, which is followed by a fast-locking phase on the microscopic faults (pulse-like rupture)