133 The aim of this study is to develop a numerical method to evaluate the optimum seismic station distribution for imaging the source rupture process of an earthquake. Based on the beam-forming technique, the source rupture distribution of an earthquake can be reconstructed through theoretical travel time correction and waveform stacking. Numerical tests show that this method successfully reco...

Using array analysis, the direction and distance to a seismic P wave source can be determined. However, individual arrays are limited in their geographical coverage and by their resolving capability, which is determined by the array aperture, configuration, and number of stations. We demonstrate these limitations on three large seismic arrays located in Japan, Europe, and California, and find t...

With seismic interferometry one can retrieve the response to a virtual source inside an unknown medium, if there is a receiver at the position of the virtual source. Using inverse scattering theory, we demonstrate that, for a 1D medium, the requirement of having an actual receiver inside the medium can be circumvented, going beyond seismic interferometry. In this case, the wavefield can be focu...

2-D elementspace and beamspaceextensions of Friedlander’s wideband interpolated root-MUSIC technique are applied to azimuthvelocity source location using real seismic data from the GERESS array (Germany). We demonstrate that 2-D interpolated rootMUSIC is able to estimate the parameters of a typical seismic source with a good accuracy. The use of interpolated root-MUSIC and its beamspace modific...

[1] The volcanic eruption on Miyakejima, Japan, in 2000 was marked by the largest earthquake swarm ever recorded in Japan, a seismicity migration accompanying a dike intrusion as the dike propagated from Miyakejima to the northwest, and formation of a caldera on Mount Oyama on Miyakejima. In this study, we propose a seismic source model which can be used to model both seismic and geodetic displ...

This report presents an interim appraisal of capabilities of a single Large Aperture Seismic Array system to perform the following functions: (i) preprocess arriving seismic signals to increase their detectability, (ii) use such preprocessed signals to perform on-line automatic detection and location, (iii) process recordings of LASA data off-line, and (iv) use the results of the off-line proce...

Accurate interpretation of seismic traveltimes and amplitudes in the exploration and global scales is complicated by the band-limited nature of seismic data. We discovered a stochastic method to reduce a seismic waveform into a most probable constituent spike train. Model waveforms were constructed from a set of candidate spike trains convolved with a source wavelet estimate. For each model wav...

Seismic inversion is a method that extracts acoustic impedance data from the seismic traces. Source wavelets are band-limited, and thus seismic traces do not contain low and high frequency information. Therefore, there is a serious problem when the deterministic seismic inversion is applied to real data and the result of deterministic inversion is smooth. Low frequency component is obtained fro...

Meng (1997) used the McConnell crosswell data to show that transmitted arrival traveltimes can be extracted from common receiver gathers (CRG) where the effective source is a tube wave generator. In this report I use the McElroy crosswell data to show that reflection events can be extracted from CRG’s where the effective source is a tube wave generator. This result opens the possibility that mo...

We investigate the effects of noise on controlled-source interferometry by multi-dimensional deconvolution (MDD) for diffusive electromagnetic (EM) fields and for seismic wavefields. An electromagnetic and a seismic dataset are modeled and decomposed into upand downwards decaying and upand downgoing fields, respectively. From the decomposed fields a reflection response can be retrieved using le...