Monitoring of Volcanic Activity by Sub-mm Geodetic Analyses

Volcanic earthquakes have been occurring beneath Zao volcano in northern Honshu, Japan since 2013, following the increase of deep low frequency earthquakes from 2012. On account of a burst of seismicity initiated in April 2015, the JMA announced a warning of eruption, however, the seismicity gradually decreased for the next two months and the warning was canceled in June. In the same time period, minor expansive deformation was observed by GNSS. Small earthquakes are still occurring, and low-freq. earthquakes (LPE) occur sometimes accompanied by static tilt changes. In this study, we try to extract the sub-mm displacements from the LPE waveforms observed by broadband seismometers (BBS) and utilize them for geodetic inversion to monitor volcanic activities. Thun et al. (2015, 2016) devised an efficient method using a running median filter (RMF) to remove LP noises, which contaminate displacement waveforms. They demonstrated the reproducibility of the waveforms corresponding to the experimentally given sub-mm displacements in the laboratory. They also apply the method to the field LPE data obtained from several volcanoes to show static displacements. The procedure is outlined as follows: (1) Unfiltered removal of the instrument response, (2) LP noise estimate by LPF with a corner frequency of 5/M, where M (seconds) is the time window of the RMF and should be at least three times the length of the rise time. (3) Subtract the noise estimated from step (2). (4) Integrate to obtain displacement waveforms. We apply the method to the BBS waveform at a distance of about 1.5 km ESE from the summit crater of Zao Volcano associated with a LPE on April 1, 2017. Assuming the time window M as 300 seconds, we successfully obtained the displacement history: taking the rise time of about 2 minutes, the site was gradually uplifted with the amount of about 50-60 µm and then subsided with HF displacements in the next 2 minutes resulting about 20-30 µm static upheaval. Comparing the waveforms obtained by the RMF and the conventional methods, the onset of the upheaval is delayed a little in the former one, while the latter one is consistent with the inclination change observed by the tiltmeter installed at the same site. It may be caused by the characteristics of the RMF and attention should be paid to detailed discussions of the source time histories.