P-wave attenuation tomography of Mount St. Helens: preliminary results from coda-normalized spectra

The P-wave velocity structure of Mount St. Helens has been recently imaged with local earthquake tomography, using the data recorded since its eruption in 1980. Part of this dataset has been processed to obtain a preliminary frequency dependent Qp image of the crust below the volcanic cone. We extended the so-called coda-normalization method, usually applied to S-waves, to the measurement of path-dependent P-wave attenuation. A scattering model has been developed in order to select the best time-window to measure the P-wave spectral amplitude on each trace. The objective is to average the effect of the source radiation pattern with the properties of early coda. The weighting matrix in the final inversion is dependent on the source radiation pattern and the spectral amplitude of noise. The Discrete Picard Condition and the Discrepancy Principle have been applied to investigate the maximum resolution available in each part of the medium. Truncated Singular Value Decomposition as well as Zeroth-, first- and second-order Tikhonov regularization techniques have been investigated by using the multi-resolution inversion code (MuReATA). The interpretation of the preliminary results is carried out by using cluster analysis on velocity and attenuation measurements.