Search USGS⌕ Search

Geology topics

Hiroyuki Kumagai

Publications and source records attributed to Hiroyuki Kumagai.

13 recordsLinked to original sources
USGS2007

Waveform inversion of volcano-seismic signals for an extended source

[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 2 L /3 and 2 W /3, or L and 2 L /5, where W and L are the crack width and length, respectively. Results from these tests indicate that the oscillating signature of the 2 L /3 and 2 W /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 2 L /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.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Waveform inversion of volcano-seismic signals assuming possible source geometries

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. Copyright 2005 by the American Geophysical Union.

Geophysical Research Letters↗
USGS2005

Response of a hydrothermal system to magmatic heat inferred from temporal variations in the complex frequencies of long-period events at Kusatsu-Shirane Volcano, Japan

We investigate temporal variations in the complex frequencies (frequency and quality factor Q) of long-period (LP) events that occurred at Kusatsu-Shirane Volcano, central Japan. We analyze LP waveforms observed at this volcano in the period between 1988 and 1995, which covers a seismically active period between 1989 and 1993. Systematic temporal variations in the complex frequencies are observed in October-November 1989, July-October 1991, and September 1992-January 1993. We use acoustic properties of a crack filled with hydrothermal fluids to interpret the observed temporal variations in the complex frequencies. The temporal variations in October-November 1989 can be divided into two periods, which are explained by a gradual decrease and increase of a gas-volume fraction in a water-steam mixture in a crack, respectively. The temporal variations in July-October 1991 can be also divided into two periods. These variations in the first and second periods are similar to those observed in November 1989 and in September-November 1992, respectively, and are interpreted as drying of a water-steam mixture and misty gas in a crack, respectively. The repeated nature of the temporal variations observed in similar seasons between July and November suggests the existence of seasonality in the occurrence of LP events. This may be caused by a seasonally variable meteoritic water supply to a hydrothermal system, which may have been heated by the flux of volcanic gases from magma beneath this volcano. ?? 2005 Elsevier B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗
USGS2005

Source process of a long-period event at Kilauea volcano, Hawaii

We analyse a long-period (LP) event observed by a dense seismic network temporarily operated at Kilauea volcano, Hawaii, in 1996. We systematically perform spectral analyses, waveform inversions and forward modeling of the LP event to quantify its source process. Spectral analyses identify two dominant spectral frequencies at 0.6 and 1.3 Hz with associated Q values in the range 10-20. Results from waveform inversions assuming six moment-tensor and three single-force components point to the resonance of a horizontal crack located at a depth of approximately 150 m near the northeastern rim of the Halemaumau pit crater. Waveform simulations based on a fluid-filled crack model suggest that the observed frequencies and Q values can be explained by a crack filled with a hydrothermal fluid in the form of either bubbly water or steam. The shallow hydrothermal crack located directly above the magma conduit may have been heated by volcanic gases leaking from the conduit. The enhanced flux of heat raised the overall pressure of the hydrothermal fluid in the crack and induced a rapid discharge of fluid from the crack, which triggered the acoustic vibrations of the resonator generating the LP waveform. The present study provides further support to the idea that LP events originate in the resonance of a crack. 2005 RAS.

Hawaii↗
USGS2003

Source mechanism of long-period events at Kusatsu-Shirane Volcano, Japan, inferred from waveform inversion of the effective excitation functions

We investigate the source mechanism of long-period (LP) events observed at Kusatsu-Shirane Volcano, Japan, based on waveform inversions of their effective excitation functions. The effective excitation function, which represents the apparent excitation observed at individual receivers, is estimated by applying an autoregressive filter to the LP waveform. Assuming a point source, we apply this method to seven LP events the waveforms of which are characterized by simple decaying and nearly monochromatic oscillations with frequency in the range 1-3 Hz. The results of the waveform inversions show dominant volumetric change components accompanied by single force components, common to all the events analyzed, and suggesting a repeated activation of a sub-horizontal crack located 300 m beneath the summit crater lakes. Based on these results, we propose a model of the source process of LP seismicity, in which a gradual buildup of steam pressure in a hydrothermal crack in response to magmatic heat causes repeated discharges of steam from the crack. The rapid discharge of fluid causes the collapse of the fluid-filled crack and excites acoustic oscillations of the crack, which produce the characteristic waveforms observed in the LP events. The presence of a single force synchronous with the collapse of the crack is interpreted as the release of gravitational energy that occurs as the slug of steam ejected from the crack ascends toward the surface and is replaced by cooler water flowing downward in a fluid-filled conduit linking the crack and the base of the crater lake. ?? 2003 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗
USGS2002

Waveform inversion of oscillatory signatures in long-period events beneath volcanoes

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 2 L / n , 2 W / n , n = 3, 5, 7, …, where L and W are the crack length and width, respectively. For the even modes with wavelengths 2 L / n , 2 W / 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.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Temporal evolution of a hydrothermal system in Kusatsu-Shirane Volcano, Japan, inferred from the complex frequencies of long-period events

We present a detailed description of temporal variations in the complex frequencies of long-period (LP) events observed at Kusatsu-Shirane Volcano. Using the Sompi method, we analyze 35 LP events that occurred during the period from August 1992 through January 1993. The observed temporal variations in the complex frequencies can be divided into three periods. During the first period the dominant frequency rapidly decreases from 5 to 1 Hz, and Q of the dominant spectral peak remains roughly constant with an average value near 100. During the second period the dominant frequency gradually increases up to 3 Hz, and Q gradually decreases from 160 to 30. During the third period the dominant frequency increases more rapidly from 3 to 5 Hz, and Q shows an abrupt increase at the beginning of this period and then remains roughly constant with an average value near 100. Such temporal variations can be consistently explained by the dynamic response of a hydrothermal crack to a magmatic heat pulse. During the first period, crack growth occurs in response to the overall pressure increase in the hydrothermal system caused by the heat pulse. Once crack formation is complete, heat gradually changes the fluid in the crack from a wet misty gas to a dry gas during the second period. As heating of the hydrothermal system gradually subsides, the overall pressure in this system starts to decrease, causing the collapse of the crack during the third period.

Kusatsu-Shirane Volcano↗
USGS2001

The dependence of acoustic properties of a crack on the resonance mode and geometry

We examine the dependence of the acoustic properties of a crack containing magmatic or hydrothermal fluids on the resonance mode and geometry to quantify the source properties of long-period (LP) events observed in volcanic areas. Our results, based on spectral analyses of synthetic waveforms generated with a fluid-driven crack model, indicate that the basic features of the dimensionless frequency (??) and quality factor (Qr) for a crack containing various types of fluids are not strongly affected by the choice of mode, although the actual ranges of Q?? and ?? both depend on the mode. The dimensionless complex frequency systematically varies with changes in the crack geometry, showing increases in both Qr and ?? as the crack length to aperture ratio decreases. The present results may be useful for the interpretation of spatial and temporal variations in the observed complex frequencies of LP events.

Geophysical Research Letters↗
USGS2001

Tomographic inversion of P-wave velocity and Q structures beneath the Kirishima volcanic complex, Southern Japan, based on finite difference calculations of complex traveltimes

We estimate the P-wave velocity and attenuation structures beneath the Kirishima volcanic complex, southern Japan, by inverting the complex traveltimes (arrival times and pulse widths) of waveform data obtained during an active seismic experiment conducted in 1994. In this experiment, six 200-250 kg shots were recorded at 163 temporary seismic stations deployed on the volcanic complex. We use first-arrival times for the shots, which were hand-measured interactively. The waveform data are Fourier transformed into the frequency domain and analysed using a new method based on autoregressive modelling of complex decaying oscillations in the frequency domain to determine pulse widths for the first-arrival phases. A non-linear inversion method is used to invert 893 first-arrival times and 325 pulse widths to estimate the velocity and attenuation structures of the volcanic complex. Wavefronts for the inversion are calculated with a finite difference method based on the Eikonal equation, which is well suited to estimating the complex traveltimes for the structures of the Kirishima volcano complex, where large structural heterogeneities are expected. The attenuation structure is derived using ray paths derived from the velocity structure. We obtain 3-D velocity and attenuation structures down to 1.5 and 0.5 km below sea level, respectively. High-velocity pipe-like structures with correspondingly low attenuation are found under the summit craters. These pipe-like structures are interpreted as remnant conduits of solidified magma. No evidence of a shallow magma chamber is visible in the tomographic images.

Geophysical Journal International↗
USGS2000

Acoustic properties of a crack containing magmatic or hydrothermal fluids

We estimate the acoustic properties of a crack containing magmatic or hydrothermal fluids to quantify the source properties of long-period (LP) events observed in volcanic areas assuming that a crack-like structure is the source of LP events. The tails of synthetic waveforms obtained from a model of a fluid-driven crack are analyzed by the Sompi method to determine the complex frequencies of one of the modes of crack resonance over a wide range of the model parameters α/ a and ρ f /ρ s , where αis the P wave velocity of the rock matrix, a is the sound speed of the fluid, and ρ f and p s are the densities of the fluid and rock matrix, respectively. The quality factor due to radiation loss ( Q r ) for the selected mode almost monotonically increases with increasing α/ a , while the dimensionless frequency ( v ) of the mode decreases with increasing α/ a and ρ f /ρ s . These results are used to estimate Q and v for a crack containing various types of fluids (gas-gas mixtures, liquid-gas mixtures, and dusty and misty gases) for values of a , ρ f , and quality factor due to intrinsic losses ( Q i ) appropriate for these types of fluids, in which Q is given by Q −l = Q r −l + Q i −1 . For a crack containing such fluids, we obtain Q ranging from almost unity to several hundred, which consistently explains the wide variety of quality factors measured in LP events observed at various volcanoes. We underscore the importance of dusty and misty gases containing small-size particles with radii around 1 μm to explain long-lasting oscillations with Q significantly larger than 100. Our results may provide a basis for the interpretation of spatial and temporal variations in the observed complex frequencies of LP events in terms of fluid compositions beneath volcanoes.

Journal of Geophysical Research B: Solid Earth↗
USGS1999

The complex frequencies of long-period seismic events as probes of fluid composition beneath volcanoes

Long-period (LP) events have been widely observed in relation to magmatic and hydrothermal activities in volcanic areas. LP waveforms characterized by their harmonic signature have been interpreted as oscillations of a fluid-filled resonator, and mixtures of liquid and gas in the form of bubbly liquids have been mainly assumed for the fluid. To investigate the characteristic properties of the resonator system, we analyse waveforms of LP events observed at four different volcanoes in Hawaii, Alaska, Colombia and Japan using a newly developed spectral method. This method allows an estimation of the complex frequencies of decaying sinusoids based on an autoregressive model. The results of our analysis show a wide variety of Q factors, ranging from tens to several hundred. We compare these complex frequencies with those predicted by the fluid-filled crack model for various mixtures of liquid, gas and ash. Although the oscillations of LP events with Q smaller than 50 can be explained by various combinations of liquids and gases, we find that ash-laden gases are required to explain long-lasting oscillations with Q larger than 100. The complex frequencies of LP events yield useful information on the types of fluids. Temporal and spatial variations of the complex frequencies can be used as probes of fluid composition beneath volcanoes.

Geophysical Journal International↗
USGS1998

The excitation and characteristic frequency of the long-period volcanic event: An approach based on an inhomogeneous autoregressive model of a linear dynamic system

We present a method to quantify the source excitation function and characteristic frequencies of long-period volcanic events. The method is based on an inhomogeneous autoregressive (AR) model of a linear dynamic system, in which the excitation is assumed to be a time-localized function applied at the beginning of the event. The tail of an exponentially decaying harmonic waveform is used to determine the characteristic complex frequencies of the event by the Sompi method. The excitation function is then derived by operating an AR filter constructed from the characteristic frequencies to the entire seismogram of the event, including the inhomogeneous part of the signal. We apply this method to three long-period events at Kusatsu-Shirane Volcano, central Japan, whose waveforms display simple decaying monochromatic oscillations except for the beginning of the events. We recover time-localized excitation functions lasting roughly 1 s at the start of each event and find that the estimated functions are very similar to each other at all the stations of the seismic network for each event. The phases of the characteristic oscillations referred to the estimated excitation function fall within a narrow range for almost all the stations. These results strongly suggest that the excitation and mode of oscillation are both dominated by volumetric change components. Each excitation function starts with a pronounced dilatation consistent with a sudden deflation of the volumetric source which may be interpreted in terms of a choked-flow transport mechanism. The frequency and Q of the characteristic oscillation both display a temporal evolution from event to event. Assuming a crack filled with bubbly water as seismic source for these events, we apply the Van Wijngaarden-Papanicolaou model to estimate the acoustic properties of the bubbly liquid and find that the observed changes in the frequencies and Q are consistently explained by a temporal change in the radii of the bubbles characterizing the bubbly water in the crack.

Journal of Geophysical Research B: Solid Earth↗