Search

A large earthquake (Mw 7.7) occurred on 16 April 2016 within the source region of the 1906 earthquake in the Ecuador-Colombia subduction zone. The 1906 event has been interpreted as a megathrust earthquake (Mw 8.8) that ruptured the source regions of smaller earthquakes in 1942, 1958, and 1979 in this subduction. Our seismic analysis indicated that the spatial distribution of the 2016 earthquake and its aftershocks correlated with patches of high interplate coupling strength and was similar to those of the 1942 earthquake and its aftershocks, suggesting that the 2016 and 1942 earthquakes ruptured the same asperity. Our analysis of tsunami waveforms of the 1906 event indicated Mw around 8.4 and showed that large slip occurred near the trench off the source regions of the above three historical and the 2016 earthquakes, suggesting that a depth-dependent complex rupture mode exists along this subduction zone.

Last August, stations on a newly deployed permanent ocean floor observation network recorded rarely seen, very low frequency signals from shallow earthquakes.

[1] We systematically used two approaches to analyze broadband seismic signals for monitoring active volcanoes: one is waveform inversion of very-long-period (VLP) signals assuming possible source mechanisms; the other is a source location method of long-period (LP) events and tremor using their amplitudes. The deterministic approach of the waveform inversion is useful to constrain the source mechanism and location but is basically only applicable to VLP signals with periods longer than a few seconds. The source location method assumes isotropic radiation of S waves and uses seismic amplitudes corrected for site amplifications. This simple approach provides reasonable source locations for various seismic signals such as a VLP event accompanying LP signals, an explosion event, and tremor associated with lahars and pyroclastic flows observed at five or fewer stations. Our results indicate that a frequency band of about 5–12 Hz and a Q factor of about 60 are appropriate for the determination of the source locations. In this frequency band the assumption of isotropic radiation may become valid because of the path effect caused by the scattering of seismic waves. The source location method may be categorized as a stochastic approach based on the nature of scattering waves. Systematic use of these two approaches provides a way to better utilize broadband seismic signals observed at a limited number of stations for improved monitoring of active volcanoes.

[1] We propose a method to investigate the dimensions and oscillation characteristics of the source of volcano-seismic signals based on waveform inversion for an extended source. An extended source is realized by a set of point sources distributed on a grid surrounding the centroid of the source in accordance with the source geometry and orientation. The source-time functions for all point sources are estimated simultaneously by waveform inversion carried out in the frequency domain. We apply a smoothing constraint to suppress short-scale noisy fluctuations of source-time functions between adjacent sources. The strength of the smoothing constraint we select is that which minimizes the Akaike Bayesian Information Criterion (ABIC). We perform a series of numerical tests to investigate the capability of our method to recover the dimensions of the source and reconstruct its oscillation characteristics. First, we use synthesized waveforms radiated by a kinematic source model that mimics the radiation from an oscillating crack. Our results demonstrate almost complete recovery of the input source dimensions and source-time function of each point source, but also point to a weaker resolution of the higher modes of crack oscillation. Second, we use synthetic waveforms generated by the acoustic resonance of a fluid-filled crack, and consider two sets of waveforms dominated by the modes with wavelengths 2L/3 and 2W/3, or L and 2L/5, where W and L are the crack width and length, respectively. Results from these tests indicate that the oscillating signature of the 2L/3 and 2W/3 modes are successfully reconstructed. The oscillating signature of the L mode is also well recovered, in contrast to results obtained for a point source for which the moment tensor description is inadequate. However, the oscillating signature of the 2L/5 mode is poorly recovered owing to weaker resolution of short-scale crack wall motions. The triggering excitations of the oscillating cracks are successfully reconstructed.

[1] We present a model of gas exsolution and bubble expansion in a melt supersaturated in response to a sudden pressure drop. In our model, the melt contains a suspension of gas bubbles of identical sizes and is encased in a penny-shaped crack embedded in an elastic solid. The suspension is modeled as a three-dimensional lattice of spherical cells with slight overlap, where each elementary cell consists of a gas bubble surrounded by a shell of volatile-rich melt. The melt is then subjected to a step drop in pressure, which induces gas exsolution and bubble expansion, resulting in the compression of the melt and volumetric expansion of the crack. The dynamics of diffusion-driven bubble growth and volumetric crack expansion span 9 decades in time. The model demonstrates that the speed of the crack response depends strongly on volatile diffusivity in the melt and bubble number density and is markedly sensitive to the ratio of crack thickness to crack radius and initial bubble radius but is relatively insensitive to melt viscosity. The net drop in gas concentration in the melt after pressure recovery represents only a small fraction of the initial concentration prior to the drop, suggesting the melt may undergo numerous pressure transients before becoming significantly depleted of gases. The magnitude of pressure and volume recovery in the crack depends sensitively on the size of the input-pressure transient, becoming relatively larger for smaller-size transients in a melt containing bubbles with initial radii less than 10−5 m. Amplification of the input transient may be large enough to disrupt the crack wall and induce brittle failure in the rock matrix surrounding the crack. Our results provide additional basis for the interpretation of volume changes in the magma conduit under Popocatepetl Volcano during Vulcanian degassing bursts in its eruptive activity in April–May 2000.

[1] We propose an approach to quantify the source of volcano-seismic signals assuming possible source geometries. Such an assumption reduces the number of free parameters in a waveform inversion, so we can quantify the source of these signals observed by a small number of seismic stations. We test this method by using the synthesized waveforms radiated by a fluid-filled crack and with seismic data observed at Kusatsu-Shirane Volcano. We find that waveform data from 2 to 3 three-component seismic stations surrounding the source enable us to estimate the source mechanism and location in the present approach.

In the early morning of 27 May 2006, local time, central Java was jolted by strong seismic ground motion. In spite of the moderate size of the earthquake (Mw = 6.4), it caused severe damage nearYogyakarta city. According to a report from the Social Affairs Ministry of Indonesia, more than 5700 people were killed, 38,000 injured, and 423,000 evacuated. As a result of the shaking, more than 126,000 buildings collapsed, and more than 392,000—including those of the famous Prambanan temple complex, located about 17 kilometers east of Yogyakarta—were severely damaged (Figure 1).

[1] The source mechanism of long-period (LP) events is examined using synthetic waveforms generated by the acoustic resonance of a fluid-filled crack. We perform a series of numerical tests in which the oscillatory signatures of synthetic LP waveforms are used to determine the source time functions of the six moment tensor components from waveform inversions assuming a point source. The results indicate that the moment tensor representation is valid for the odd modes of crack resonance with wavelengths 2L/n, 2W/n, n = 3, 5, 7, …, where L and W are the crack length and width, respectively. For the even modes with wavelengths 2L/n, 2W/n, n = 2, 4, 6, …, a generalized source representation using higher-order tensors is required, although the efficiency of seismic waves radiated by the even modes is expected to be small. We apply the moment tensor inversion to the oscillatory signatures of an LP event observed at Kusatsu-Shirane Volcano, central Japan. Our results point to the resonance of a subhorizontal crack located a few hundred meters beneath the summit crater lakes. The present approach may be useful to quantify the source location, geometry, and force system of LP events, and opens the way for moment tensor inversions of tremor.