Mapping pressurized volcanic fluids from induced crustal seismic velocity drops

See allHide authors and affiliations

Science  04 Jul 2014:
Vol. 345, Issue 6192, pp. 80-82
DOI: 10.1126/science.1254073

Seismic noise reveals volcanic plumbing

Monitoring the way in which seismic noise passes through Earth's crust after a large earthquake can clarify how volcanoes erupt. Japan has the highest-density seismic network in the world. Brenguier et al. observed reductions in seismic velocity below volcanic regions of Japan from before, to the weeks and months after the 2011 Tohoku-Oki earthquake (see the Perspective by Prejean and Haney). This indicates that pressurized fluids below volcanoes can weaken in response to dynamic stress perturbations.

Science, this issue p. 80; see also p. 39


Volcanic eruptions are caused by the release of pressure that has accumulated due to hot volcanic fluids at depth. Here, we show that the extent of the regions affected by pressurized fluids can be imaged through the measurement of their response to transient stress perturbations. We used records of seismic noise from the Japanese Hi-net seismic network to measure the crustal seismic velocity changes below volcanic regions caused by the 2011 moment magnitude (Mw) 9.0 Tohoku-Oki earthquake. We interpret coseismic crustal seismic velocity reductions as related to the mechanical weakening of the pressurized crust by the dynamic stress associated with the seismic waves. We suggest, therefore, that mapping seismic velocity susceptibility to dynamic stress perturbations can be used for the imaging and characterization of volcanic systems.

Large volcanic eruptions are preceded by long-term pressure buildup in volcano magmatic and hydrothermal systems. Therefore, knowledge of the extent and state of these pressurized volcanic fluids at depth will help in the better anticipation of future eruptions. In particular, seismic tomography is often used to delineate volcano-feeding systems at different scales (1, 2). However, a major difficulty of traditional seismic imaging of volcanoes is that the geological contrasts of the host rock might dominate the final tomographic images, which are only partially sensitive to the content and state of volcanic fluids (3).

Recent geodetic observations have shown that volcanic areas are characterized by anomalous responses to crustal deformation induced by large earthquakes, as demonstrated by the subsidence of volcanoes in Chile and Japan after the 2010 Maule and 2011 Tohoku-Oki earthquakes (4, 5). This sensitivity to strong coseismic deformation and shaking is probably associated with the presence of pressurized hydrothermal and magmatic fluids at depth in a fractured medium. We explored the responses of volcanoes to transient stress perturbations by investigating the temporal evolution of crustal seismic velocities in Japan in response to the seismic shaking and deformation caused by the March 2011 moment magnitude (Mw) 9.0 Tohoku-Oki earthquake.

The Hi-net, Japanese high-sensitivity seismograph network, is among the densest in the world; thus, the 2011 Tohoku-Oki earthquake remains the best-recorded large earthquake to date. It was associated with large, widespread static ground deformation and ground shaking (Fig. 1). In this study, we used seismic noise-based monitoring (6) to characterize the response of the upper crust to the coseismic shaking and deformation caused by the earthquake. We analyzed 1 year of continuous seismic records from a portion of the dense Hi-net seismic network (600 stations, as shown in the inset to Fig. 1A), spanning from 6 months before to 6 months after the earthquake occurrence.

Fig. 1 Static strain and ground shaking caused by the Tohoku-Oki earthquake.

(A) Modeled coseismic dilatation static strain at 5 km in depth (7). The red star shows the position of the epicenter of the Tohoku-Oki earthquake. (Inset) Positions of the Hi-net seismic stations (red points). (B) Averaged peak ground velocity measured using the KiK-net strong-motion network (7).

We computed the daily vertical-vertical noise cross-correlation functions using a processing scheme that minimized the effects of the strong aftershock activity that followed the Tohoku-Oki earthquake (7). To avoid the choice of an arbitrary reference cross-correlation function and to improve the precision of the measurements, we separately computed velocity changes for all of the possible daily cross-correlation functions for each station pair. Using a Bayesian least-squares inversion, we retrieved accurate daily continuous velocity change time series for every station pair (7).

We computed the seismic velocity changes averaged over the day of the Tohoku-Oki earthquake and 4 days after, relative to the seismic velocity changes time series averaged over 6 months before the Tohoku-Oki earthquake (Fig. 2A) (7). These changes mainly correspond to the response of the upper crust to the coseismic shaking and deformation. Similar to previous studies of coseismic velocity variations (6, 8), a reduction in velocity was widespread over Honshu Island. Furthermore, the strongest velocity drops were not observed in the area closest to the epicenter or within large sedimentary basins, as would be expected. The patterns of the observed velocity reductions did not correlate with the intensity of the ground shaking or with the coseismic deformation (Fig. 1); instead, the strongest coseismic velocity reductions occurred under volcanic regions. In particular, a large part of the Tohoku volcanic front and the Mt. Fuji volcanic region are well delineated.

Fig. 2 Crustal seismic velocity perturbations caused by the Tohoku-Oki earthquake.

(A) Coseismic crustal seismic velocity changes induced by the 2011 Tohoku-Oki earthquake. (Inset) Velocity changes averaged over 5 days preceding the day of the Tohoku-Oki earthquake. (B) Seismic velocity susceptibility computed using the seismic velocity changes shown in (A). Black triangles denote Quaternary period volcanoes, and the red line depicts the main volcanic fronts.

The mechanism by which seismic velocities decrease in response to stress perturbations is commonly described as related to the opening of cracks (9, 10), which might also induce an increase in permeability and a transfer of fluids at depth and may lead to further triggering of earthquakes. Over long distances, large earthquakes are known to trigger anomalous hydrothermal activity (11), aftershocks on a global scale (12), tectonic tremor activity (13), and slow-slip events (14). The origin for this remote triggering of activity is believed to be the associated dynamic stress that is caused by the passing of the seismic waves.

We used the approach of Gomberg and Agnew (15) to estimate the level of dynamic strain Δξ and dynamic stress Δσ from the observed peak ground velocity (PGV), such that Δξ ≈ ν/c and Δσ ≈ μν/c, where μ is the mean crustal shear modulus (~30 × 109 Pa), ν is the PGV (measured by the KiK-net, strong-motion seismograph network installed in boreholes together with the Hi-net sensors), and c is the mean wave phase velocity of the Rayleigh waves that propagate within the upper crust (~3 km/s). The dynamic strain caused by the passing of the seismic waves was one to two orders of magnitude higher than the static coseismic strain for Honshu Island. We thus conclude that the dynamic stress associated with the seismic waves emitted by the Tohoku-Oki earthquake was the main cause of the large seismic velocity reductions under the volcanic regions—in particular, the Mt. Fuji area, where the static stress change can be considered negligible. We then defined the seismic velocity susceptibility as the ratio between the observed reductions in the seismic velocity and the estimated dynamic stress. The distribution of these seismic velocity susceptibilities correlates with the main volcanic areas (Fig. 2B).

The sensitivity of the seismic velocity to stress changes in the rock increases with decreasing effective pressure (16, 17). Under volcanic areas, the effective pressure in the crust can be reduced because of the presence of highly pressurized hydrothermal and magmatic volcanic fluids at depth. We thus argue that the observed strong coseismic velocity reductions delineate the regions where such pressurized volcanic fluids are present in the upper crust. An important implication of our observation is that the seismic velocity susceptibility to stress can be used as a proxy to the level of pressurization of the hydrothermal and/or magmatic fluids in volcanic areas. So far, this susceptibility is greatest in the Mt. Fuji area and along the Tohoku volcanic arc, where it reached 15 × 10−4 MPa−1, whereas it is more than 10 times smaller for the Cretaceous stiff plutonic regions of eastern Tohoku (Fig. 2B).

Fluids are also known to have important roles in earthquake nucleation (3). The volcanic areas where large seismic velocity susceptibility was observed were also characterized by large triggered seismic activity after the Tohoku-Oki earthquake (18, 19), including a particularly strong (magnitude 6.4) earthquake that occurred 4 days after the main shock, near Mt. Fuji. This confirms that the crust in these areas is quite sensitive to strong transient stress perturbations. We argue that mapping the susceptibility of seismic velocities to dynamic stress changes can be used to image and characterize regions with low effective pressure, such as volcanic systems.

Supplementary Materials

Materials and Methods

Figs. S1 and S2

References (2032)

References and Notes

  1. See supplementary materials on Science Online.
  2. Acknowledgments: All of the seismological data used in this study are archived at the Japanese National Research Institute for Earth Science and Disaster Prevention (NIED) ( and The seismological data used in this study are available upon request at the NIED data center. We thank Geospatial Information Authority of Japan for access to Global Navigation Satellite System data. All of the computations in this study were performed using the High-Performance Computing (CIMENT) infrastructure (, which is supported by the Rhône-Alpes region (GRANT CPER07_13 CIRA:, France Grilles (, and the CNRS MASTODONS program. We acknowledge the French National Research Agency and the Japan Science and Technology Agency for funding support (project 2011-JAPN-006-NAMAZU). We received support from the European Union through the projects European Research Council advanced grant 227507 “Whisper” and FP7 RI283542 “VERCE.” We thank T. Nishimura, P. Johnson, G. Olivier, G. Hillers, C. Jaupart, and T. Lecocq for useful discussions.
View Abstract

Navigate This Article