Preparatory and precursory processes leading up to the 2014 phreatic eruption of Mount Ontake, Japan

Abstract

We analyzed seismicity linked to the 2014 phreatic eruption of Mount Ontake, Japan, on 27 September 2014. We first relocated shallow volcano tectonic (VT) earthquakes and long-period (LP) events from August to September 2014. By applying a matched-filter technique to continuous waveforms using these relocated earthquakes, we detected numerous additional micro-earthquakes beneath the craters. The relocated VT earthquakes aligned on a near-vertical plane oriented NNW-SSE, suggesting they occurred around a conduit related to the intrusion of magmatic-hydrothermal fluids into the craters. The frequency of VT earthquakes gradually increased from 6 September 2014 and reached a peak on 11 September 2014. After the peak, seismicity levels remained elevated until the eruption. b-values gradually increased from 1.2 to 1.7 from 11 to 16 September 2014 then declined gradually and dropped to 0.8 just before the eruption. During the 10-min period immediately preceding the phreatic eruption, VT earthquakes migrated in the up-dip direction as well as laterally along the NNW-SSE feature. The migrating seismicity coincided with an accelerated increase of pre-eruptive tremor amplitude and with an anomalous tiltmeter signal that indicated summit upheaval. Therefore, the migrating seismicity suggests that the vertical conduit was filled with pressurized fluids, which rapidly propagated to the surface during the final 10 min before the eruption.

Kato et al. Earth, Planets and Space (2015) 67:111 DOI 10.1186/s40623-015-0288-x

FRONTIER LETTER Open Access

http://crossmark.crossref.org/dialog/?doi=10.1186/s40623-015-0288-x&domain=pdf

Web End = Preparatory and precursory processes leading up to the 2014 phreatic eruption of Mount Ontake, Japan

Aitaro Kato1,2*, Toshiko Terakawa1, Yoshiko Yamanaka1, Yuta Maeda1, Shinichiro Horikawa1, Kenjiro Matsuhiro1 and Takashi Okuda1

Abstract

We analyzed seismicity linked to the 2014 phreatic eruption of Mount Ontake, Japan, on 27 September 2014. We first relocated shallow volcano tectonic (VT) earthquakes and long-period (LP) events from August to September 2014. By applying a matched-filter technique to continuous waveforms using these relocated earthquakes, we detected numerous additional micro-earthquakes beneath the craters. The relocated VT earthquakes aligned on a near-vertical plane oriented NNWSSE, suggesting they occurred around a conduit related to the intrusion of magmatichydrothermal fluids into the craters. The frequency of VT earthquakes gradually increased from 6 September 2014 and reached a peak on 11 September 2014. After the peak, seismicity levels remained elevated until the eruption. b-values gradually increased from 1.2 to 1.7 from 11 to 16 September 2014 then declined gradually and dropped to 0.8 just before the eruption. During the 10-min period immediately preceding the phreatic eruption, VT earthquakes migrated in the up-dip direction as well as laterally along the NNWSSE feature. The migrating seismicity coincided with an accelerated increase of pre-eruptive tremor amplitude and with an anomalous tiltmeter signal that indicated summit upheaval. Therefore, the migrating seismicity suggests that the vertical conduit was filled with pressurized fluids, which rapidly propagated to the surface during the final 10 min before the eruption.

Keywords: Mount Ontake; 2014 phreatic eruption; Precursor; Earthquake; Relocation; Matched-filter technique; b-value

Findings

Introduction

Mount Ontake is located at the southern end of the Norikura volcanic zone, central Japan (Fig. 1a). In 1979, Mount Ontake produced its first phreatic eruption in recorded history, ejecting ash, gas, and debris into the sky (volcanic explosivity index (VEI) = 2, Smithsonian 2015). The 1979 eruption was followed by minor phreatic explosions in 1991 and 2007 (e.g., Nakamichi et al. 2009; Japan Meteorological 2014). On 27 September 2014, at 11:52 a.m. JST (UTC + 9), Mount Ontake produced another phreatic (steam type) eruption, spraying volcanic ash, gas, and debris on surrounding areas. Based on ejecta volume and the maximum ash column height of

~7000 m, the VEI for this eruption was 2 (Earthquake Research Institute in the University of Tokyo 2014; Japan Meteorological 2014). The 2014 eruption resulted in the deaths of 57 climbers, with a further 6 missing and presumed dead (e.g., Cyranoski 2014). The 2014 phreatic eruption had similar features to the 1979 eruption, including eruption style and rough positions of vents.

It is generally difficult to predict the timing and likely size (VEI) of phreatic eruptions because their measurable precursors are weak and highly localized (e.g., Barberi et al. 1992) or absent (e.g., Maeda et al. 2015). Phreatic eruptions often occur during periods of elevated seismicity and high heat flow through the volcanic system, which can be potential precursors (e.g., Barberi et al. 1992; Aoyama and Oshima 2008; Mordret et al. 2010). However, the fundamental processes related to preparatory and precursory stages of phreatic eruptions are poorly known. Because seismicity is one of the most powerful tools for identifying volcanic processes that could precede phreatic

* Correspondence: mailto:[email protected]

Web End [email protected]

1Earthquake and Volcano Research Center, Graduate School of Environmental Studies, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8601, Japan

2Earthquake Research Institute, University of Tokyo, 1-1-1 Yayoi, Bunkyo-ku, Tokyo 113-0032, Japan

2016 Kato et al. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0

Web End =http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.

Kato et al. Earth, Planets and Space (2015) 67:111 Page 2 of 11

Fig. 1 Tectonic setting of Mount Ontake. a Bathymetry of the central Japan region. Red triangles denote active volcanoes on the Japanese islands. Red contours are isodepth contours of the upper surface of the subducting Philippine Sea plate at 10-km intervals (Hirose et al. 2008), and blue dashed contours denote the upper surface of the Pacific plate at 50-km intervals (Nakajima and Hasegawa 2006). The black square indicates the location of the study area (Fig. 1b). b Topographic map of Mount Ontake. Black circles denote the epicenters of relocated events (Fig. 4). Blue and red squares represent short-period (1-s natural resonance) and broadband seismic stations whose waveforms were analyzed in the present study. Stations with prefix NU, V, and N have been operated by Nagoya University, JMA, and NIED, respectively. ROPW and GNDT have been operated by Nagano Prefecture and Gifu Prefecture. Seismometer of V.ONTN, N.ASSH and N.KADH are emplaced in a borehole at a depth of ~100 m. Green arrow shows downward tilt changes recorded during the final 6 min before the 2014 phreatic eruption

eruptions (e.g., Chouet 1996; Roman and Cashman 2006), it is essential to investigate the spatialtemporal evolution of earthquakes leading up to the 2014 phreatic eruption of Mount Ontake.

A manually constructed earthquake catalog for the region under the summit of Mount Ontake shows that an average monthly rate of volcano tectonic (VT) and long-period (LP) earthquakes (magnitude < 1.0) during 1 year before the eruption was about 3.3/month, and0.17/month, respectively (Japan Meteorological 2014). Then, the micro-seismicity sharply increased beginning around 2 weeks before the eruption. To constrain the fundamental processes associated with the phreatic eruption, we relocated VT and LP earthquakes beneath the summit, using double differential travel times extracted from waveform cross-correlation method. We then searched for new earthquakes by applying a matched-filter detection technique (e.g., Shelly et al. 2007; Kato et al. 2013) to continuous seismograms recorded near the Mount Ontake summit (Fig. 1b) from 23 August to 30 September 2014, using the waveforms of well-located earthquakes as template events. We discuss the spatialtemporal evolution of the elevated seismicity, including temporal changes in b-values from the beginning of September to the time of eruption. Furthermore, we focus on the immediate precursors that began less than 10 min before the phreatic eruption, including the up-dip and lateral migrations of

earthquake hypocenters, an accelerated increase in the amplitudes of pre-eruptive tremors, and tiltmeter changes observed near the summit.

Data and methods

We used three-component velocity seismograms retrieved at 13 seismic stations around the summit of Mount Ontake, operated by Nagoya University, the Japan Meteorological Agency (JMA), the National Research Institute for Earth Science and Disaster Prevention (NIED), Nagano Prefecture, and Gifu Prefecture (Fig. 1b). Data were continuously recorded by each station at a sampling rate of 100 Hz (the station ROPW restarted the recording from 19 September 2014). Most earthquakes before the eruption were classified as VT earthquakes, having clear high-frequency P- and S-wave arrivals with broadband frequencies greater than 5 Hz (Fig. 2a). In contrast to VT earthquakes, LP events were characterized by a weak high-frequency onset, quasi-monochromatic oscillations with a peak frequency near 1.5 Hz, and durations of more than 20 s (Fig. 2b). Comparisons of the spectrograms show a marked difference between VT and LP earthquakes in terms of the spectral amplitude of the seismogram versus frequency. These classifications are similar to those proposed previously for other volcanoes (e.g., Lahr et al. 1994; Chouet 1996).

Kato et al. Earth, Planets and Space (2015) 67:111 Page 3 of 11

VT : 2014/09/30 05:13:38 (M 0.3)V.ONTN (NS)

a

-6

Velocity (m/s) Frequency (Hz)

2.5

0 0 2.5 5 7.5 10 12.5 15 17.5 20 22.5 25 27.5 30

0

-2.5

20

x 10

2.5 5 7.5 10 12.5 15 17.5 20 22.5 25 27.5

15

10

5

b

-6

LP : 2014/09/24 17:00:11 (M 0.2)

Velocity (m/s) Frequency (Hz)

5

0 0 2.5 5 7.5 10 12.5 15 17.5 20 22.5 25 27.5 30

V.ONTN (NS)

0

-5

20

x 10

2.5 5 7.5 10 12.5 15 17.5 20 22.5 25 27.5

15

10

5

Time (s)

Fig. 2 Velocity seismograms and associated spectrograms of a VT earthquake (a) and an LP event (b) recorded at the V.ONTN station (NS component). Each spectrogram is calculated from the spectral amplitudes of velocity versus frequency in a moving 1-s window running through the entire seismogram (98 % overlap), as shown above each spectrogram

Relocation of template earthquakes

We visually picked the arrival times of the P- and S-waves of 96 VT earthquakes between 31 August and 30 September 2014 and of 2 LP earthquakes on 16 and 24 September 2014. The number of P- and S-wave arrival times for each event exceeded 8 and 4, respectively. The picking accuracy of P- and S-wave arrivals for VT earthquakes was estimated to be 0.04 s and 0.1 s, respectively. We relocated these 98 template earthquakes by applying a relocation algorithm (Zhang and Thurber 2003) to both arrival times and differential arrival times fixing the velocity model. Our P-wave velocity model was constructed from previous seismic surveys (Fig. 3), assuming a Poisson ratio of 0.25. For depths above sea level, we used a P-wave velocity structure determined by seismic surveys at other active volcanoes in Japan because the velocity structure of Mount Ontake has not been previously investigated (e.g., Tanaka et al. 2002; Aoki et al. 2009). For the deeper velocity structure, we used a regional velocity model for central Japan based on Ikami et al. (1986) and Kato et al. (2007).

The numbers of absolute arrival times used in the relocation are 979 for P-waves and 722 for S-waves. The numbers of differential arrival times for the manually picked P- and S-waves are 6697 and 4003, respectively. We also used differential arrival times obtained by waveform cross-correlation method. The correlation measurements were conducted using a 0.8-s window length beginning 0.4 s before each visually picked arrival time for waveform data bandpass filtered between 8 and 16 Hz. This produced a data set of differential arrival times that contained 2477 P-wave and 367 S-wave observation pairs, all of which had a normalized cross-correlation coefficient of 0.90. The root mean square of the travel time residual was reduced from 0.118 to 0.06 s after 21 iterations.

To assess the reliability of the relocated hypocenters, we applied a statistical resampling approach termed as bootstrap method to the events. We calculated synthetic arrival time data by adding random Gaussian noises with a standard deviation of 0.05 s (average of 0 s) to the observed arrival time data. We then relocated events

Kato et al. Earth, Planets and Space (2015) 67:111 Page 4 of 11

Vp (km/s)

0.0 2.0 4.0 6.0

Depth (km)

2

0

2

4

6

8

2

0

2

4

6

8

Fig. 3 One-dimensional P-wave velocity structure employed for the earthquake relocation. We assumed a Poissons ratio of 0.25

to determine the shift in location based on the corrupted arrival times. The process was repeated 500 times with new Gaussian noise generated for each iteration. Based on the cumulative results, the average horizontal and depth errors are 160 and 340 m, respectively. The error ellipsoids of each hypocenter on the horizontal surface are slightly elongated along the NESW directions. The hypocenter accuracy of the LP events are not so well compared with that of the VT events, due to a shortage of high-frequency components at onsets of P- and S-wave arrivals. Furthermore, absolute locations were sensitive to the assumed velocity structure, but the relative distribution of earthquake was less sensitive. For an example, relocated hypocenters tend to be deeper and more clustered for an assumed velocity model with a small velocity gradient near the surface.

Detection of earthquakes by a matched-filter technique

To more precisely characterize pre- and post-eruption temporal changes in seismicity, we searched for additional events using a matched-filter technique (e.g., Shelly et al. 2007; Kato et al. 2012, 2013). Waveform data from the relocatable earthquakes were used to construct a dictionary

of template events. We then matched these waveforms against continuous data recorded between 23 August 2014 and 30 September 2014.

To search for VT earthquakes, both continuous data and template waveforms were bandpass filtered from 4 to 12 Hz and decimated to 50 Hz. Similarly, for LP earthquakes, data were bandpass filtered from 0.5 to 1 Hz and decimated to 10 Hz. As template events, we used the relocated earthquakes in Fig. 4. A synthetic S-wave arrival time for each template event was calculated using the one-dimensional seismic velocity structure used for relocation. To obtain template event waveforms and target waveforms, we used a 4.0-s window for VT earthquakes and a 20.0-s window for LP earthquakes, beginning 2.0 s prior to the synthetic S-wave arrival time. We computed correlation coefficients between a template event waveform and a target waveform using incremental detection-window shifts of 0.02 s for VT earthquakes and0.1 s for LP events. We set a threshold for event detection equal to 10 times the median absolute deviation (MAD) of the average correlation coefficients calculated throughout the day of interest for VT earthquakes and 13 times the MAD for LP events. To remove multiple detections, we allocated the location of the detected event to that of the template event with the highest correlation coefficient within the 4-s window for VT earthquakes and 20-s window LP events. Fig. 5 shows an example of waveforms at the detection time for a template event. We identified 3300 VT and 28 LP earthquakes (Fig. 6). For more details of the matched-filter technique, see Kato et al. (2013).

Calculation of b-values

Using the newly constructed earthquake catalog, we calculated b-values using the maximum likelihood estimate (Utsu 1965) in overlapping 400-event windows shifted by 20 events per calculation. We simultaneously determined the temporal change in the catalog completeness magnitude using the goodness-of-fit method (Wiemer and Wyss 2000). We used the 90 % goodness-of-fit level as the completeness magnitude, which implies that 90 % of the variability in the observed frequencymagnitude distribution can be modeled with a synthetic power law.

Results

Figure 4a shows that the relocated hypocenters are tightly clustered to the southwest of the Mount Ontake summit, with ~500-m absolute offset. The VT and LP events occurred in the same area at shallow depths close to sea level. The relocated hypocenters (98 events) roughly coincide with alignments of eruptive vents identified by remote sensing images obtained by Geospatial Information Authority of Japan (Geospatial Information Authority of Japan 2014) (white masked area in Fig. 4a). The epicentral distribution was oriented roughly NNWSSE. Depth

Kato et al. Earth, Planets and Space (2015) 67:111 Page 5 of 11

Fig. 4 Earthquake relocations. a Epicentral distribution of the relocated earthquakes superimposed on topography. Yellow and blue circles denote VT earthquakes before and after the 2014 eruption, respectively. Circles are scaled to magnitude. Yellow stars represent LP events prior to the eruption. White masked areas and a red triangle are eruptive craters identified by remote sensing (Geospatial Information Authority of Japan 2014) and the summit of Mount Ontake, respectively. b Vertical depth section of the relocated hypocenters perpendicular to the NNWSSE orientation. The positive vertical axis is below sea level

section for relocated hypocenters shows that the earthquakes are aligned along a vertically dipping plane (Fig. 4b) approximately 1.5 km high (dip direction) and 0.7 km wide. Additionally, while VT earthquakes before the eruption were located at greater depths, post-eruption VT earthquakes occurred at depths shallower than the preceding seismicity by approximately 1.0 km.

Figure 6a, b shows spacetime diagrams of the newly detected seismicity. From 6 September 2014, the source region of VT earthquakes slightly expanded vertically as well as laterally along the NNWSSE feature, accompanied by a gradual increase in VT and LP earthquakes. The VT earthquake rate peaked on 11 September and remained relatively high until the eruption, with a slight decay (Fig. 6c). The number of LP events showed a sharp increase between 15 and 16 September, a delay of about 5 days from the peak of VT seismicity.

After the eruption, the depths of VT earthquakes became systematically shallower by ~1.0 km relative to pre-eruption hypocenters (Fig. 6b), even after adding numerous events with the matched-filter technique. Furthermore, we recognized a slight lateral expansion of the earthquake alignment in the NNWSSE direction just after the eruption (Fig. 6a).

Based on the newly detected seismicity, we identified several immediate precursors to the phreatic eruption of27 September 2014. About 10 min prior to the eruption, hypocenters of VT earthquakes migrated in the up-dip direction as well as laterally along the NNWSSE

feature in a near-vertical plane (Fig. 7a, b). The rate of lateral earthquake migration ranged from 0.35 to 0.7 m/s (Fig. 7a). Along the depth section, the VT earthquakes become gradually shallower from 11:40 to 11:45 (JST) on 27 September then rapidly rose toward the craters (Fig. 7b). Simultaneously, the amplitudes of the pre-eruptive tremors increased at an accelerated rate (Fig. 7c). Each spike in the pre-eruptive tremor corresponds to a VT earthquake detected by the matched-filter technique, resulting in a gradual increase in large-magnitude events (Fig. 7a, b). In addition to both earthquake migrations and accelerating tremor amplitude, tiltmeter signals were also characterized by downward acceleration in the ESE direction, beginning ~6 min before the phreatic eruption (Fig. 1b and 7d).

After the peak in VT seismicity on 11 September 2014, b-values gradually increased from ~1.2 to ~1.7 on 16 September 2014 (Fig. 8). Shortly after attaining this maximum, b-values declined until the eruption, with only minor fluctuations. This reduction in b-values is well correlated with an increased frequency of relatively large-magnitude events after 21 September 2014 (Fig. 6d). The final drop in b-values was due to the relative increase in large-magnitude events during the 10 min immediately preceding the eruption. The representative plots of frequencymagnitude distributions shown in Fig. 8b indicate that the temporal changes in b-values described above are statistically significant.

Kato et al. Earth, Planets and Space (2015) 67:111 Page 6 of 11

Fig. 5 An example of waveforms at the detection time of M-1.1 event using a template event (2014/09/17. 09:34:47.75). The gray and red lines are the continuous and template event waveforms for seismic stations aligned according to epicentral distance, respectively. The window of the template event is 4 s, starting from 2.0 s prior to the synthetic S-wave arrival time. The station names and components are given to the left of each trace. Waveforms were bandpass filtered between 4 and 12 Hz, and the amplitudes of waveforms were normalized by the maximum value at each station

Discussion and conclusions

By relocating shallow VT and LP earthquakes and applying a matched-filter detection technique to seismic data recorded before the 2014 phreatic eruption of Mount Ontake, we reconstructed the earthquake catalog beneath the craters, which provides new insights into

the preparatory and precursory stages of the eruption. The relocated VT earthquakes are aligned on a near-vertical plane oriented NNWSSE, defining a near-vertical conduit related to the intrusion of hydrothermal fluids (steam, gases, and liquid water) into the craters (Fig. 4).

Kato et al. Earth, Planets and Space (2015) 67:111 Page 7 of 11

Eruption

N15W

0.4

2.0

1.5

1.0

0.5

0.0

0.5

1.0

0

a

Long-period earthquake (LP)

0.2

Depth (km BSL) VTnumber Distance (km)

S15E

0.0

0.2

M1 M0 M-1

0.4

b

1.0 0.5 0.0 Depth (km BSL)

3000

c

50

LPnumber

40

2000

30

1000

20

LP

10

0

2

2

d

Magnitude

1

1

0

0

1

1

2

2

26 31 06 11 16 21 26Date 2014 Sep.

2014 Aug.

Fig. 6 Seismicity detected by matched-filter processing. a Spacetime diagram of all detected events before and after the 2014 phreatic eruption, from 23 August to 30 September. Blue circles and red stars denote VT earthquakes and LP events, respectively. The diagram shows earthquake origin times and locations projected onto the strike of the hypocenter alignments. b Depth variations in detected VT earthquakes over time. Circles are scaled to magnitude and color-coded to depth (below sea level (BSL), with positive values being below sea level). c Cumulative number of VT earthquakes (blue curve) and LP events (red curve) over time. d Magnitude versus time

The initial acceleration of VT seismicity, as well as LP events during the middle of September, suggests that infiltration of hot fluids triggered the abundant seismicity beneath the craters (Fig. 9a). Following the peak in the occurrence rate of VT earthquakes, shallow LP events sharply increased. It is widely believed that the source mechanism of shallow LP events is resonant excitation of a fluid-filled crack (e.g., Chouet 1996; Kumagai and Chouet 2000; Waite et al. 2008). In the present case, the excitation could be driven by hydrothermal fluids. We thus infer that hydrothermal fluids infiltrated into preexisting faults/cracks during the period of LP activity.

Before the minor phreatic eruption in 2007, a very-long-period (VLP) earthquake with equivalent moment magnitude 3.4, as well as a burst of LP events, was recorded by broadband seismic stations (Nakamichi et al. 2009). Waveform inversion yielded a location 0.6 km above sea level beneath the craters of Mount Ontake and suggests that the source mechanism was inflation followed by deflation of the crack associated with vaporization of water and the subsequent discharge of steam into the hydrothermal system. Interestingly, the azimuth of the crack exciting the VLP (N20W) was sub-parallel to that of the vertical conduit indicated by the relocated hypocenters

Kato et al. Earth, Planets and Space (2015) 67:111 Page 8 of 11

in the present study (N15W). Whether reactivation of an old conduit or coincidental similarity, the similar geometries suggest the preferred orientation of pre-existing cracks along NNWSSE beneath the summit of Mount Ontake. However, a VLP event was not recorded during the elevated VT seismicity before the 2014 eruption.

During the period of increased LP seismicity, b-values increased to nearly 1.7 (Fig. 8). Based on spatial mapping of b-values at many volcanoes (e.g., Wyss et al. 1997; McNutt 2005; Bridges and Gao 2006), high b-values have been widely observed in areas with high thermal gradients due to the interaction of hot fluids (vapors, gases, and liquid water) with adjacent magma bodies. Thus, we interpret that the high b-values reflect an increase in the density of smaller faults/cracks due to the infiltration of hot fluids (Fig. 9a).

Between 11 and 16 September 2014, the combination of an increase in LP events and an increase in b-values suggests that hot fluids infiltrated into pre-existing faults/ cracks above a magma chamber. Fluid pressures within the faults/cracks further increased, resulting in a reduction of the effective normal stress on faults/cracks (e.g., Hubbert and Rubey 1959; Terakawa et al. 2013; Terakawa 2014). This reduction in fault strength facilitated the growth of small faults/cracks into larger ones, producing higher-magnitude earthquakes (Fig. 6d), which can also explain the gradual decrease in b-values after 16 September 2014 (Fig. 8).

Based on these observations, we postulate that hot fluids were pervasive to pre-existing faults/cracks and pressurized by heat supplied from the underlying magma chamber from the middle of September 2014 (Fig. 9a). The shallow portion beneath the summit perhaps withstood stress concentration derived from the increase in internal pressure of the underlying hydrothermal system.

During the 10 min prior to the phreatic eruption, the up-dip and lateral migrations of VT earthquakes coincided with an increase in pre-eruptive tremor amplitude and an anomalous tiltmeter signal near the craters. Because the tiltmeter signal was consistent with summit upheaval (Figs. 1 and 9b), the migrating seismicity indicates the rapid propagation of pressurized fluids (vapor, gases, and liquid water) through a vertical conduit. The lateral propagating speed of the conduit tip ranged from 0.35 to 0.7 m/s (Fig. 7a). Interestingly, these speeds are roughly comparable to those of migrating VT earthquakes triggered by dike intrusion (0.060.3 m/s) (e.g., Hayashi and Morita 2003; Shelly and Hill 2011). Finally, failure of the sealed surface of the vertical conduit led to the 2014 phreatic eruption. The sealed surface may be produced by the precipitation of altered clay minerals or hydrothermal secondary deposits near the surface (e.g., Bertrand et al. 2012; Maeda et al. 2015).

Eruption

0.4

1.0

0.5

0.0

1.0

0.5

0.0

0.70 m/s

0.35 m/s

M1

M0

M-1

0.2

Distance (km)

Depth (km BSL)

0.0

1.0

0.5

0.0

Depth (km BSL)

0.2

a

0.4

5.0 m/s

0.7 m/s

0.3 m/s

b

Velocity (4-12 Hz)

Envelope (1-4 Hz)

c

0.1mm/s

0.02mm/s

V.ONTN

d

N, E up

Tilt: NS

Tilt: EW

1nrad

V.ONTN

11:30 11:40 11:50 12:002014/09/27

Time

Fig. 7 Precursors of the 2014 phreatic eruption. a Horizontal distance of detected VT earthquakes from 11:30 to 12:00 JST on 27 September 2014, projected onto the strike of the hypocenter alignments (N15W-S15E). Circles are scaled to magnitude and color-coded to depth. Black dashed lines approximate the locations of the fronts of earthquake migrations. b Depth variations in VT earthquakes, colored to depth. Circles are scaled to magnitude. Black dashed lines approximate the locations of the fronts of earthquake migrations. c Bandpass filtered waveform between 4 and 12 Hz (black curve) and envelope between 1 and 4 Hz (red curve) recorded at the V.ONTN station. d Time series of tiltmeter records observed at the V.ONTN station. NS and EW components denote the northward and eastward ground-up components, respectively

Kato et al. Earth, Planets and Space (2015) 67:111 Page 9 of 11

Eruption

2.0

a

b-value

1.5

1.0

0.5

26 31 06 11 16 21 26Date 2014 Sep.

2014 Aug.

b

b = 0.81 (2014/09/22 to Eruption)

100

1.5 1.0 0.5 0.0 0.5 1.0 1.5 Magnitude

100

1.5 1.0 0.5 0.0 0.5 1.0 1.5

100

1.5 1.0 0.5 0.0 0.5 1.0 1.5

Frequency

100

1.5 1.0 0.5 0.0 0.5 1.0 1.5

100

1.5 1.0 0.5 0.0 0.5 1.0 1.5

100

1.5 1.0 0.5 0.0 0.5 1.0 1.5

10

10

10

10

10

10

b = 1.73 (2014/09/15 to 09/18)

b = 1.19 (2014/09/18 to 09/23)

1

1

1

1

1

1

Fig. 8 Temporal changes in b-values. a Calculated b-values are represented by yellow squares. Uncertainties are indicated by vertical and horizontal error bars. b-value uncertainties (vertical error bars) were estimated following Shi and Bolt (1982). The horizontal error bars indicate the time window used for each b-value calculation (containing 400 events). Colored squares represent the selected magnitudefrequency distributions given by the triangles in (b). b Selected cumulative magnitudefrequency distributions

a

Middle of September ~10 minutes before the eruption

b

WSW ENE WSW ENE

Water table

Propagation of crack

Tilt meter

Tilt meter

LP events

Low b-value Infiltration of

fluids along large faults

Hydrothermal

system

Up-dip migration of VT events

VT events

Water table

High b-value

Infiltration of hot

fluids along small faults

Heat supply

Heat supply

Hydrothermal

system

Magma

Magma

Fig. 9 Schematic figure showing the preparation and precursory stages of the 2014 phreatic eruption of Mount Ontake. a Around the middle of September 2014 (1 to 2 weeks prior to the eruption). b Approximately 10 min before the eruption. Tilt change showing down in the ESE direction is consistent with opening of the shallow vertical conduit

Kato et al. Earth, Planets and Space (2015) 67:111 Page 10 of 11

To our knowledge, the present study was the first to describe the rapid propagation of a vertical conduit during the accelerating stage prior to a phreatic eruption. It is important to capture short-term precursors as early as possible and to raise the alarm for climbers to evacuate from potential eruption sites. Further understanding of the excitation mechanism of phreatic eruptions will require continuous monitoring of seismicity and crustal movement in the vicinity of the summit because the geophysical signals related to the pre-eruption propagation of a vertical crack are too localized to be observed in the far field.

Competing interestsThe authors declare that they have no competing interests.

Authors contributionsAK carried out the data processing, phase picking, relocation, and analysis of the continuous data and drafted the manuscript. TT carried out the seismic observation and the data processing and helped draft the manuscript. YY carried out the seismic observation and the data processing and helped draft the manuscript. YM helped draft the manuscript. SH, KM, and TO operated and maintained the seismic stations. All authors read and approved the final manuscript.

AcknowledgementsWe thank JMA, NIED, Nagano Prefecture, and Gifu Prefecture for allowing us to use waveform data collected by their permanent stations. JMA provided the tiltmeter data recorded near the summit (V.ONTN). We are grateful to Y. Shibayama for the dedicated help with the data processing. We thank D. Roman and one anonymous reviewer for useful comments and suggestions. Figures were created using GMT (Wessel and Smith 1995) and spectrogram analysis tool (Miyakawa and Sakai 2008). This study was supported by the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan, under its Earthquake and Volcano Hazards Observation and Research Program and under KAKEN (26900002). This study was also supported by the Earthquake Research Institute cooperative research program.

Received: 15 April 2015 Accepted: 2 July 2015

ReferencesAoki Y, Takeo M, Aoyama H, Fujimatsu J, Matsumoto S, Miyamachi H, Nakamichi H,

Ohkura T, Ohminato T, Oikawa J, Tanada R, Tsutsui T, Yamamoto K, Yamamoto M, Yamasato H, Yamawaki T (2009) P-wave velocity structure beneath Asama Volcano, Japan, inferred from active source seismic experiment. J Volcanol Geotherm Res 187:2727Aoyama H, Oshima H (2008) Tilt change recorded by broadband seismometer prior to small phreatic explosion of Meakan-dake volcano, Hokkaido, Japan. Geophys Res Lett 35:L06307. doi:http://dx.doi.org/10.1029/2007GL032988

Web End =10.1029/2007GL032988 Barberi F, Bertagnini A, Landi P, Principe C (1992) A review on phreatic eruptions and their precursors. J Volcanol Geotherm Res 52(4):23146. doi:http://dx.doi.org/10.1016/0377-0273(92)90046-G

Web End =10.1016/ http://dx.doi.org/10.1016/0377-0273(92)90046-G

Web End =0377-0273(92)90046-G Bertrand EA, Caldwell TG, Hill GJ, Wallin EL, Bennie SL, Cozens N, Onacha SA,

Ryan GA, Walter C, Zaino A, Wameyo P (2012) Magnetotelluric imaging of upper-crustal convection plumes beneath the Taupo Volcanic Zone, New Zealand. Geophys Res Lett 39:L02304. doi:http://dx.doi.org/10.1029/2011GL050177

Web End =10.1029/2011GL050177 Bridges DL, Gao S (2006) Spatial variation of seismic b-values beneath Makushin Volcano, Unalaska Island, Alaska. Earth Planet Sci Lett, 245, 408415 doi:http://dx.doi.org/10.1016/j.epsl.2006.03.010

Web End =10.1016/j.epsl.2006.03.010.

Chouet B (1996) Long-period volcano seismicity: its source and use in eruption forecasting. Nature 380:30916Cyranoski D (2014) Why Japan missed volcanos warning signs. Nature. doi:http://dx.doi.org/10.1038/nature.2014.16022

Web End =10.1038/nature.2014.16022 Earthquake Research Institute in the University of Tokyo (2014), http://www.data.jma.go.jp/svd/vois/data/tokyo/STOCK/kaisetsu/CCPVE/shiryo/130/130_no01.pdf

Web End =http:// http://www.data.jma.go.jp/svd/vois/data/tokyo/STOCK/kaisetsu/CCPVE/shiryo/130/130_no01.pdf

Web End =www.data.jma.go.jp/svd/vois/data/tokyo/STOCK/kaisetsu/CCPVE/shiryo/130/ http://www.data.jma.go.jp/svd/vois/data/tokyo/STOCK/kaisetsu/CCPVE/shiryo/130/130_no01.pdf

Web End =130_no01.pdf . Accessed 1 April 2015.

Geospatial Information Authority of Japan (2014), http://www.gsi.go.jp/BOUSAI/h26-ontake-index.html

Web End =http://www.gsi.go.jp/BOUSAI/ http://www.gsi.go.jp/BOUSAI/h26-ontake-index.html

Web End =h26-ontake-index.html . Accessed 1 April 2015.

Hayashi Y, Morita Y (2003) An image of a magma intrusion process inferred from precise hypocentral migrations of the earthquake swarm east of the Izu Peninsula. Geophys J Int 153:15974

Hirose F, Nakajima J, Hasegawa A (2008) Three-dimensional seismic velocity structure and configuration of the Philippine Sea slab in southwestern Japan estimated by double-difference tomography. J Geophys Res 113:B09315. doi:http://dx.doi.org/10.1029/2007JB005274

Web End =10.1029/2007JB005274

Hubbert MK, Rubey WW (1959) Role of fluid pressure in mechanics of overthrust faulting 1. Mechanics of fluid-filled porous solids and its application to overthrust faulting. Geol Soc Am Bull 70(2):11566

Ikami A, Yoshii T, Kubota S, Sasaki Y, Hasemi A, Moriya T, Miyamachi H, Matsuura

R, Wada K (1986) A seismic refraction profile in and around Nagano Prefecture, Central Japan. J Phys Earth 34:45774Japan Meteorological Agency (2014), Rep. Coordin. Committee on Prediction of

Volcanic Eruption, http://www.data.jma.go.jp/svd/vois/data/tokyo/STOCK/kaisetsu/CCPVE/shiryo/130/130_no01.pdf

Web End =http://www.data.jma.go.jp/svd/vois/data/tokyo/STOCK/ http://www.data.jma.go.jp/svd/vois/data/tokyo/STOCK/kaisetsu/CCPVE/shiryo/130/130_no01.pdf

Web End =kaisetsu/CCPVE/shiryo/130/130_no01.pdf . Accessed 1 April 2015.

Kato A, Iidaka T, Kurashimo E, Nakagawa S, Hirata N, Iwasaki T (2007)

Delineation of probable asperities on the Atotsugawa fault, central Japan, using a dense temporary seismic network. Geophys Res Lett 34:L09318. doi:http://dx.doi.org/10.1029/2007GL029604

Web End =10.1029/2007GL029604 Kato A, Obara K, Igarashi T, Tsuruoka H, Nakagawa S, Hirata N (2012) Propagation of slow slip leading up to the 2011 Mw 9.0 Tohoku-Oki earthquake. Science 335:7058. doi:http://dx.doi.org/10.1126/science.1215141

Web End =10.1126/science.1215141 Kato A, Fukuda J, Obara K (2013) Response of seismicity to static and dynamic stress changes induced by the 2011 M9.0 Tohoku-Oki earthquake. Geophys Res Lett 40, doi:http://dx.doi.org/10.1002/grl.50699

Web End =10.1002/grl.50699 .

Kumagai H, Chouet BA (2000) Acoustic properties of a crack containing magmatic or hydrothermal fluids. J Geophys Res 105(B11):25,49325,512 Lahr JC, Chouet BA, Stephens CD, Power JA, Page RA (1994) Earthquake classification, location and error analysis in a volcanic environment: implications for the magmatic system of the 19891990 eruptions at Redoubt Volcano, Alaska. J Volcanol Geotherm Res 62:13751Maeda Y, Kumagai H, Lacson RJ, Figueroa MS, Yamashina T, Ohkura T, Baloloy AV

(2015) A phreatic explosion model inferred from a very long period seismic event at Mayon Volcano, Philippines. J Geophys Res Solid Earth 120:22642. doi:http://dx.doi.org/10.1002/2014JB011440

Web End =10.1002/2014JB011440 McNutt SR (2005) Volcanic seismology. Annu Rev Earth Planet Sci 32:46191 Miyakawa K, Sakai S (2008) Development of a spectrogram analysis tool for seismic waveform data and its application to MeSO-net for noise survey, Technical Research Report, Earthquake Research Institute, the University of Tokyo, 14, 1322.

Mordret A, Jolly AD, Duputel Z, Fournier N (2010) Monitoring of phreatic eruptions using Interferometry on Retrieved Cross-Correlation Function from Ambient Seismic Noise: results from Mt. Ruapehu, New Zealand. J Volcanol Geotherm Res 191:4659Nakajima J, Hasegawa A (2006) Anomalous low-velocity zone and linear alignment of seismicity along it in the subducted Pacific slab beneath Kanto, Japan: reactivation of subducted fracture zone? Geophys Res Lett 33:L16309. doi:http://dx.doi.org/10.1029/2006GL026773

Web End =10.1029/2006GL026773 Nakamichi H, Kumagai H, Nakano M, Okubo M, Kimata F, Ito Y, Obara K (2009) Source mechanism of very-long-period event at Mt. Ontake, central Japan: response of a hydrothermal system to magma intrusion beneath the summit. J Volcanol Geotherm Res 187:16777Roman DC, Cashman KV (2006) The origin of volcano-tectonic earthquake swarms. Geology 34:45760Shelly DR, Hill DP (2011) Migrating swarms of brittlefailure earthquakes in the lower crust beneath Mammoth Mountain, California. Geophys Res Lett 38:L20307. doi:http://dx.doi.org/10.1029/2011GL049336

Web End =10.1029/2011GL049336 Shelly DR, Beroza GC, Ide S (2007) Non-volcanic tremor and low-frequency earthquake swarms. Nature 446:3057. doi:http://dx.doi.org/10.1038/nature05666

Web End =10.1038/nature05666 Shi Y, Bolt BA (1982) The standard error of the magnitude-frequency b-value. Bull

Seism Soc Am 72:167787Smithsonian Institution (2015) Bulletin of the Global Volcanism Network, http://www.volcano.si.edu/search_eruption.cfm

Web End =http://www.volcano.si.edu/search_eruption.cfm . Accessed 1 April 2015. Tanaka S, Hamaguchi H, Nishimura T, Yamawaki T, Ueki S, Nakamichi H, Tsutsui T,

Miyamachi H, Matsuwo N, Oikawa J, Ohminato T, Miyaoka K, Onizawa S, Mori T, Aizawa K (2002) Three-dimensional P-wave velocity structure of Iwate volcano, Japan from active seismic survey. Geophys Res Lett 29, 10.1029/2002GL014983. Terakawa T (2014) Evolution of pore fluid pressures in a stimulated geothermal reservoir inferred from earthquake focal mechanisms. Geophys Res Lett 41:746876. doi:http://dx.doi.org/10.1002/2014GL061908

Web End =10.1002/2014GL061908

Kato et al. Earth, Planets and Space (2015) 67:111 Page 11 of 11

Terakawa T, Yamanaka Y, Nakamichi H, Watanabe T, Yamazaki F, Horikawa S,

Okuda T (2013) Effects of pore fluid pressure and tectonic stress on diverse seismic activities around the Mt. Ontake volcano, central Japan. Tectonophysics 608:13848Utsu T (1965) A method for determining the value of b in a formula logn = a - bM showing the magnitude frequency for earthquakes. Geophys Bull Hokkaido Univ 13:99103Waite GP, Chouet BA, Dawson PB (2008) Eruption dynamics at Mount St. Helens imaged from broadband seismic waveforms: interaction of the shallow magmatic and hydrothermal systems. J Geophys Res 113, B02305, doi:http://dx.doi.org/10.1029/2007JB005259

Web End =10.1029/2007JB005259 .

Wessel P, Smith WHF (1995) New version of the generic mapping tools released.

Eos Trans AGU 76:329Wiemer S, Wyss M (2000) Minimum magnitude of completeness in earthquake catalogs: examples from Alaska, the western United States, and Japan. Bull Seismol Soc Am 90:85969Wyss M, Shimazaki K, Wiemer S (1997) Mapping active magma chambers by b-value beneath the off-Ito volcano, Japan. J Geophys Res 102:2041322 Zhang H, Thurber CH (2003) Double-difference tomography: the method and its application to the Hayward fault, California. Bull Seismol Soc Am 93:187589. doi:http://dx.doi.org/10.1785/0120020190

Web End =10.1785/0120020190

Submit your manuscript to a journal and benet from:

7 Convenient online submission7 Rigorous peer review7 Immediate publication on acceptance7 Open access: articles freely available online 7 High visibility within the eld7 Retaining the copyright to your article

Submit your next manuscript at 7 springeropen.com