Multi-scale reasonable attenuation tomography analysis (MuRAT): An imaging algorithm designed for volcanic regions

L. De Siena*, C. Thomas, R. Aster

*Corresponding author for this work

Research output: Contribution to journalArticle

11 Citations (Scopus)

Abstract

The attenuation of body-wave amplitudes with propagation distance can be used to provide detailed tomographic images of seismic interfaces, fluid reservoirs, and melt batches in the crust The high sensitivity of body-wave energies to high-scattering structures becomes an obstacle when we try to apply attenuation tomography to small-scale volcanic media, where we must take into account the complexities induced by strong heterogeneous scattering, topography, and uncertain source modeling in the recorded wave-fields. The MuRAT code uses a source- and site-independent coda-normalization method to obtain frequency-dependent measurements of P-to-coda and S-to-coda energy ratios. The code inverts these data for both the geometrical spreading factor and the spatially-dependent quality factors (Q), providing additional attenuation information in the regions where velocity tomography is available. The high sensitivity of coda-waves to highly heterogeneous structures highlights zones of anomalous scattering, which may corrupt amplitude-dependent attenuation measurements, and where basal assumptions of linear optics may go unfulfilled. A multi-step tomographic inversion increases the stability of the results obtained in regions of high heterogeneity (e.g., the volcanic edifice) by the inclusion of data corresponding to either sources or stations located in regions of lower heterogeneity. On the other hand, a mere increase in the number of rays entirely contained in the heterogeneous structures affects both the stability and the effective resolution of the results. We apply the code to two small waveform datasets recorded at an active (Mount St. Helens) and at a quiescent (Mount Vesuvius) volcano. The results show that the seismicity located inside or under the volcanic edifice produces an increase of the low-frequency energy ratios with travel time in both areas. In our interpretation, the anomalous concentration of energy which affects any waveform recorded on the cone, produced inside the volcanic edifice or in the feeding system of the volcano is due to seismic source- or medium-dependent resonance. The results also provide spatial and frequency limits on the feasibility of attenuation tomography in these two regions with larger datasets. (C) 2014 Elsevier B.V. All rights reserved.

Original languageEnglish
Pages (from-to)22-35
Number of pages14
JournalJournal of Volcanology and Geothermal Research
Volume277
Early online date8 Apr 2014
DOIs
Publication statusPublished - 1 May 2014

Keywords

  • attenuation tomography
  • multi-step code
  • scattering
  • volcano imaging
  • local earthquake tomography
  • Mount-ST-Helens
  • active seismic experiment
  • velocity structure
  • P-wave
  • multiple-scattering
  • 3-dimensional P
  • Mt. Vesuvius
  • inversion
  • topography

Cite this

Multi-scale reasonable attenuation tomography analysis (MuRAT) : An imaging algorithm designed for volcanic regions. / De Siena, L.; Thomas, C.; Aster, R.

In: Journal of Volcanology and Geothermal Research, Vol. 277, 01.05.2014, p. 22-35.

Research output: Contribution to journalArticle

@article{c7b8aa98fc524d6c9539f9808ef32d8c,
title = "Multi-scale reasonable attenuation tomography analysis (MuRAT): An imaging algorithm designed for volcanic regions",
abstract = "The attenuation of body-wave amplitudes with propagation distance can be used to provide detailed tomographic images of seismic interfaces, fluid reservoirs, and melt batches in the crust The high sensitivity of body-wave energies to high-scattering structures becomes an obstacle when we try to apply attenuation tomography to small-scale volcanic media, where we must take into account the complexities induced by strong heterogeneous scattering, topography, and uncertain source modeling in the recorded wave-fields. The MuRAT code uses a source- and site-independent coda-normalization method to obtain frequency-dependent measurements of P-to-coda and S-to-coda energy ratios. The code inverts these data for both the geometrical spreading factor and the spatially-dependent quality factors (Q), providing additional attenuation information in the regions where velocity tomography is available. The high sensitivity of coda-waves to highly heterogeneous structures highlights zones of anomalous scattering, which may corrupt amplitude-dependent attenuation measurements, and where basal assumptions of linear optics may go unfulfilled. A multi-step tomographic inversion increases the stability of the results obtained in regions of high heterogeneity (e.g., the volcanic edifice) by the inclusion of data corresponding to either sources or stations located in regions of lower heterogeneity. On the other hand, a mere increase in the number of rays entirely contained in the heterogeneous structures affects both the stability and the effective resolution of the results. We apply the code to two small waveform datasets recorded at an active (Mount St. Helens) and at a quiescent (Mount Vesuvius) volcano. The results show that the seismicity located inside or under the volcanic edifice produces an increase of the low-frequency energy ratios with travel time in both areas. In our interpretation, the anomalous concentration of energy which affects any waveform recorded on the cone, produced inside the volcanic edifice or in the feeding system of the volcano is due to seismic source- or medium-dependent resonance. The results also provide spatial and frequency limits on the feasibility of attenuation tomography in these two regions with larger datasets. (C) 2014 Elsevier B.V. All rights reserved.",
keywords = "attenuation tomography, multi-step code, scattering, volcano imaging, local earthquake tomography, Mount-ST-Helens, active seismic experiment, velocity structure, P-wave, multiple-scattering, 3-dimensional P, Mt. Vesuvius, inversion, topography",
author = "{De Siena}, L. and C. Thomas and R. Aster",
year = "2014",
month = "5",
day = "1",
doi = "10.1016/j.jvolgeores.2014.03.009",
language = "English",
volume = "277",
pages = "22--35",
journal = "Journal of Volcanology and Geothermal Research",
issn = "0377-0273",
publisher = "Elsevier",

}

TY - JOUR

T1 - Multi-scale reasonable attenuation tomography analysis (MuRAT)

T2 - An imaging algorithm designed for volcanic regions

AU - De Siena, L.

AU - Thomas, C.

AU - Aster, R.

PY - 2014/5/1

Y1 - 2014/5/1

N2 - The attenuation of body-wave amplitudes with propagation distance can be used to provide detailed tomographic images of seismic interfaces, fluid reservoirs, and melt batches in the crust The high sensitivity of body-wave energies to high-scattering structures becomes an obstacle when we try to apply attenuation tomography to small-scale volcanic media, where we must take into account the complexities induced by strong heterogeneous scattering, topography, and uncertain source modeling in the recorded wave-fields. The MuRAT code uses a source- and site-independent coda-normalization method to obtain frequency-dependent measurements of P-to-coda and S-to-coda energy ratios. The code inverts these data for both the geometrical spreading factor and the spatially-dependent quality factors (Q), providing additional attenuation information in the regions where velocity tomography is available. The high sensitivity of coda-waves to highly heterogeneous structures highlights zones of anomalous scattering, which may corrupt amplitude-dependent attenuation measurements, and where basal assumptions of linear optics may go unfulfilled. A multi-step tomographic inversion increases the stability of the results obtained in regions of high heterogeneity (e.g., the volcanic edifice) by the inclusion of data corresponding to either sources or stations located in regions of lower heterogeneity. On the other hand, a mere increase in the number of rays entirely contained in the heterogeneous structures affects both the stability and the effective resolution of the results. We apply the code to two small waveform datasets recorded at an active (Mount St. Helens) and at a quiescent (Mount Vesuvius) volcano. The results show that the seismicity located inside or under the volcanic edifice produces an increase of the low-frequency energy ratios with travel time in both areas. In our interpretation, the anomalous concentration of energy which affects any waveform recorded on the cone, produced inside the volcanic edifice or in the feeding system of the volcano is due to seismic source- or medium-dependent resonance. The results also provide spatial and frequency limits on the feasibility of attenuation tomography in these two regions with larger datasets. (C) 2014 Elsevier B.V. All rights reserved.

AB - The attenuation of body-wave amplitudes with propagation distance can be used to provide detailed tomographic images of seismic interfaces, fluid reservoirs, and melt batches in the crust The high sensitivity of body-wave energies to high-scattering structures becomes an obstacle when we try to apply attenuation tomography to small-scale volcanic media, where we must take into account the complexities induced by strong heterogeneous scattering, topography, and uncertain source modeling in the recorded wave-fields. The MuRAT code uses a source- and site-independent coda-normalization method to obtain frequency-dependent measurements of P-to-coda and S-to-coda energy ratios. The code inverts these data for both the geometrical spreading factor and the spatially-dependent quality factors (Q), providing additional attenuation information in the regions where velocity tomography is available. The high sensitivity of coda-waves to highly heterogeneous structures highlights zones of anomalous scattering, which may corrupt amplitude-dependent attenuation measurements, and where basal assumptions of linear optics may go unfulfilled. A multi-step tomographic inversion increases the stability of the results obtained in regions of high heterogeneity (e.g., the volcanic edifice) by the inclusion of data corresponding to either sources or stations located in regions of lower heterogeneity. On the other hand, a mere increase in the number of rays entirely contained in the heterogeneous structures affects both the stability and the effective resolution of the results. We apply the code to two small waveform datasets recorded at an active (Mount St. Helens) and at a quiescent (Mount Vesuvius) volcano. The results show that the seismicity located inside or under the volcanic edifice produces an increase of the low-frequency energy ratios with travel time in both areas. In our interpretation, the anomalous concentration of energy which affects any waveform recorded on the cone, produced inside the volcanic edifice or in the feeding system of the volcano is due to seismic source- or medium-dependent resonance. The results also provide spatial and frequency limits on the feasibility of attenuation tomography in these two regions with larger datasets. (C) 2014 Elsevier B.V. All rights reserved.

KW - attenuation tomography

KW - multi-step code

KW - scattering

KW - volcano imaging

KW - local earthquake tomography

KW - Mount-ST-Helens

KW - active seismic experiment

KW - velocity structure

KW - P-wave

KW - multiple-scattering

KW - 3-dimensional P

KW - Mt. Vesuvius

KW - inversion

KW - topography

U2 - 10.1016/j.jvolgeores.2014.03.009

DO - 10.1016/j.jvolgeores.2014.03.009

M3 - Article

VL - 277

SP - 22

EP - 35

JO - Journal of Volcanology and Geothermal Research

JF - Journal of Volcanology and Geothermal Research

SN - 0377-0273

ER -