129 Long Period (LP) events, divided into two families of similar events, were recorded by the 50 stations deployed on Mount Etna in the second half of June 2008. During this period lava was flowing from a lateral fracture after a summit strombolian
eruption. In order to understand the mechanisms of these events, we perform moment tensor
During the second half of June, 2008, 50 broadband seismic stations were deployed on Mt Etna volcano in close proximity to the summit, allowing us to observe seismic activity with exceptionally high resolution. 129 long period events (LP) with domina
nt frequencies ranging between 0.3 and 1.2 Hz, were extracted from this dataset. These events form two families of similar waveforms with different temporal distributions. Event locations are performed by cross-correlating signals for all pairs of stations in a two-step scheme. In the first step, the absolute location of the centre of the clusters was found. In the second step, all events are located using this position. The hypocentres are found at shallow depths (20 to 700 m deep) below the summit craters. The very high location resolution allows us to detect the temporal migration of the events along a dike-like structure and 2 pipe shaped bodies, yielding an unprecedented view of some elements of the shallow plumbing system at Mount Etna. These events do not seem to be a direct indicator of the ongoing lava flow or magma upwelling.
We investigate the possibility to extract information contained in seismic waveforms propagating in fluid-filled porous media by developing and using a full waveform inversion procedure valid for layered structures. To reach this objective, we first
solve the forward problem by implementing the Biot theory in a reflectivity-type simulation program. We then study the sensitivity of the seismic response of stratified media to the poroelastic parameters. Our numerical tests indicate that the porosity and consolidation parameter are the most sensitive parameters in forward and inverse modeling, whereas the permeability has only a very limited influence on the seismic response. Next, the analytical expressions of the sensitivity operators are introduced in a generalized least-square inversion algorithm based on an iterative modeling of the seismic waveforms. The application of this inversion procedure to synthetic data shows that the porosity as well as the fluid and solid parameters can be correctly reconstructed as long as the other parameters are well known. However, the strong seismic coupling between some of the model parameters makes it difficult to fully characterize the medium by a multi-parameter inversion scheme. One solution to circumvent this difficulty is to combine several model parameters according to rock physics laws to invert for composite parameters. Another possibility is to invert the seismic data for the perturbations of the medium properties, such as those resulting from a gas injection.
Analytical formulas are derived to compute the first-order effects produced by plane inhomogeneities on the point source seismic response of a fluid-filled stratified porous medium. The derivation is achieved by a perturbation analysis of the poro-el
astic wave equations in the plane-wave domain using the Born approximation. This approach yields the Frechet derivatives of the P -- SV - and SH-wave responses in terms of the Greens functions of the unperturbed medium. The accuracy and stability of the derived operators are checked by comparing, in the time-distance domain, differential seismograms computed from these analytical expressions with complete solutions obtained by introducing discrete perturbations into the model properties. For vertical and horizontal point forces, it is found that the Frechet derivative approach is remarkably accurate for small and localized perturbations of the medium properties which are consistent with the Born approximation requirements. Furthermore, the first-order formulation appears to be stable at all source-receiver offsets. The porosity, consolidation parameter, solid density and mineral shear modulus emerge as the most sensitive parameters in forward and inverse modeling problems. Finally, the Amplitude-Versus-Angle response of a thin layer shows strong coupling effects between several model parameters.
An array of ten broadband stations was installed on the Popocatepetl volcano (Mexico) for five months between October 2002 and February 2003. 26 regional and teleseismic earthquakes were selected and filtered in the frequency time domain to extract t
he fundamental mode of the Rayleigh wave. The average dispersion curve was obtained in two steps. Firstly, phase velocities were measured in the period range [2-50] s from the phase difference between pairs of stations, using Wiener filtering. Secondly, the average dispersion curve was calculated by combining observations from all events in order to reduce diffraction effects. The inversion of the mean phase velocity yielded a crustal model for the volcano which is consistent with previous models of the Mexican Volcanic Belt. The overall crustal structure beneath Popocatepetl is therefore not different from the surrounding area, and the velocities in the lower crust are confirmed to be relatively low. Lateral variations of the structure were also investigated by dividing the network into four parts and by applying the same procedure to each sub-array. No well-defined anomalies appeared for the two sub-arrays for which it was possible to measure a dispersion curve. However, dispersion curves associated with individual events reveal important diffraction for 6 s to 12 s periods which could correspond to strong lateral variations at 5 to 10 km depth.
