[1] Here we present the results of local source tomographic inversion beneath central Java. The data set was collected by a temporary seismic network. More than 100 stations were operated for almost half a year. About 13,000 P and S arrival times from 292 events were used to obtain three-dimensional (3-D) Vp, Vs, and Vp/Vs models of the crust and the mantle wedge beneath central Java. Source location and determination of the 3-D velocity models were performed simultaneously based on a new iterative tomographic algorithm, LOTOS-06. Final event locations clearly image the shape of the subduction zone beneath central Java. The dipping angle of the slab increases gradually from almost horizontal to about 70°. A double seismic zone is observed in the slab between 80 and 150 km depth. The most striking feature of the resulting P and S models is a pronounced low-velocity anomaly in the crust, just north of the volcanic arc (Merapi-Lawu anomaly (MLA)). An algorithm for estimation of the amplitude value, which is presented in the paper, shows that the difference between the fore arc and MLA velocities at a depth of 10 km reaches 30% and 36% in P and S models, respectively. The value of the Vp/Vs ratio inside the MLA is more than 1.9. This shows a probable high content of fluids and partial melts within the crust. In the upper mantle we observe an inclined low-velocity anomaly which links the cluster of seismicity at 100 km depth with MLA. This anomaly might reflect ascending paths of fluids released from the slab. The reliability of all these patterns was tested thoroughly.Citation: Koulakov, I., et al. (2007), P and S velocity structure of the crust and the upper mantle beneath central Java from local tomography inversion,

SUMMARY The seismic structure of the stratovolcano Merapi (Java, Indonesia) was studied using an active seismic experiment. Three 3 km long seismic profiles each consisting of up to 30 three‐component seismometers with an interstation distance of 100 m were built up in an altitude range between 1000 and 2000 m above sea level. The detailed study of the seismic properties of the propagation media in active volcanic regions is important to understand the natural seismic signals used for eruption forecasting. The seismic experiment at Merapi therefore concentrates on the heterogeneous structure within a radius of 5 km from the active dome, where the sources of most of the natural volcanic seismic events are located. The cone of Merapi volcano consists of different materials changing on a small scale due to the layering of eruptive material. Additionally, the topography of the erosion valleys leads to an irregular deposition, which cannot be described by a simple 1‐D layering. These inhomogeneities have a strong influence on seismic signals. The direct P and S waves are attenuated quickly and show only small amplitudes on seismograms. The energy lost from the direct waves, however, is not changed into heat but scattered and can be observed as seismic coda following the direct waves. The observed seismograms show a spindle‐like amplitude increase after the direct P phase. This shape of the envelope can be explained by the diffusion model. According to this model there are so many strong inhomogeneities that the direct wave can be neglected and all energy is concentrated in multiple scattered waves. Besides the envelope, the coherence and polarization properties of the wavefield also indicate strong scattering. Only the first onset shows coherence over a station spacing of 100 m, whereas the late phases carrying the major part of the energy are mainly incoherent. The horizontal components of the seismograms have larger amplitudes than the vertical component, but within the horizontal plane the polarization is almost arbitrary, corresponding to waves arriving from scatterers located arbitrarily in space. As a result of the inversion using the diffusion model we obtain values of the S‐wave scattering attenuation coefficient, ηs, and the S‐wave intrinsic absorption coefficient, ηi. In the frequency range of 4–20 Hz used in this study the scattering attenuation is at least one order of magnitude larger than the intrinsic absorption (ηs≫ηi). The mean free path of S waves is as low as 100 m (ηs−1≈100 m). The scattering coefficient is independent of frequency (ηs∼f0.0), whereas the coefficient of intrinsic attenuation increases with increasing frequency (ηi∼f1.6). The natural seismic signals at Merapi volcano show similar characteristics to the artificial shots. The first onsets have only small amplitudes and the energy maximum arrives delayed compared to the direct waves. Therefore, these signals appear to be strongly affected by multiple scattering also.

Mt. Merapi is one of the most dangerous volcanoes in Indonesia, located within the tectonically active region of south‐central Java. This study investigates how Mt. Merapi affected ‐ and was affected by ‐ nearby tectonic earthquakes. In 2001, a Mw6.3 earthquake occurred in conjunction with an increase in fumarole temperature at Mt. Merapi. In 2006, another Mw6.3 earthquake took place, concomitant with an increase of magma extrusion and pyroclastic flows. Here, we develop theoretical models to study the amount of stress transfer between the earthquakes and the volcano, showing that dynamic, rather than static, stress changes are likely responsible for the temporal and spatial proximity of these events. Our examination of the 2001 and 2006 events implies that volcanic activity at Mt. Merapi is influenced by stress changes related to remote tectonic earthquakes, a finding that is important for volcano hazard assessment in this densely inhabited area.

We infer temporal changes in the elastic properties of the edifice of Merapi volcano (Java, Indonesia) before its eruption in 1998 by analyzing multiply scattered elastic waves excited by a repeatable controlled seismic source. A pre‐eruptive increase of shear wave velocity, which correlates well with pre‐eruptive seismicity and dome‐growth is revealed. The method can be used as a “pressure‐gauge” for pressure changes inside of volcanoes, because increasing pressures in rocks are known to cause proportionally increasing elastic wave velocities.

No abstract