Signal restoration through deconvolution applied to deep mantle seismic probes

Research output: Contribution to journalArticle

8 Citations (Scopus)

Abstract

We present a method of signal restoration to improve the signal-to-noise ratio, sharpen seismic arrival onset, and act as an empirical source deconvolution of specific seismic arrivals. Observed time-series gi are modelled as a convolution of a simpler time-series fi, and an invariant point spread function (PSF) h that attempts to account for the earthquake source process. The method is used on the shear wave time window containing SKS and S, whereby using a Gaussian PSF produces more impulsive, narrower, signals in the wave train. The resulting restored time-series facilitates more accurate and objective relative traveltime estimation of the individual seismic arrivals. We demonstrate the accuracy of the reconstruction method on synthetic seismograms generated by the reflectivity method. Clean and sharp reconstructions are obtained with real data, even for signals with relatively high noise content. Reconstructed signals are simpler, more impulsive, and narrower, which allows highlighting of some details of arrivals that are not readily apparent in raw waveforms. In particular, phases nearly coincident in time can be separately identified after processing. This is demonstrated for two seismic wave pairs used to probe deep mantle and core-mantle boundary structure: (1) the Sab and Scd arrivals, which travel above and within, respectively, a 200-300-km-thick, higher than average shear wave velocity layer at the base of the mantle, observable in the 88-92 deg epicentral distance range and (2) SKS and SPdiffKS, which are core waves with the latter having short arcs of P-wave diffraction, and are nearly identical in timing near 108-110 deg in distance. A Java/ Matlab algorithm was developed for the signal restoration, which can be downloaded from the authors web page, along with example data and synthetic seismograms.

Original languageEnglish (US)
Pages (from-to)1353-1362
Number of pages10
JournalGeophysical Journal International
Volume167
Issue number3
DOIs
StatePublished - Dec 2006

Fingerprint

Deconvolution
deconvolution
restoration
Restoration
arrivals
Time series
Earth mantle
Shear waves
probe
Optical transfer function
mantle
probes
synthetic seismogram
seismograms
time series
point spread functions
Seismic waves
S waves
S-wave
Convolution

Keywords

  • Core-mantle boundary
  • DDeconvolution
  • Lower mantle
  • Regularization
  • Total variation

ASJC Scopus subject areas

  • Geochemistry and Petrology
  • Geophysics

Cite this

Signal restoration through deconvolution applied to deep mantle seismic probes. / Stefan, W.; Garnero, Edward; Renaut, Rosemary.

In: Geophysical Journal International, Vol. 167, No. 3, 12.2006, p. 1353-1362.

Research output: Contribution to journalArticle

@article{e95a3103c6dd4d1c84589fa6554e8ec8,
title = "Signal restoration through deconvolution applied to deep mantle seismic probes",
abstract = "We present a method of signal restoration to improve the signal-to-noise ratio, sharpen seismic arrival onset, and act as an empirical source deconvolution of specific seismic arrivals. Observed time-series gi are modelled as a convolution of a simpler time-series fi, and an invariant point spread function (PSF) h that attempts to account for the earthquake source process. The method is used on the shear wave time window containing SKS and S, whereby using a Gaussian PSF produces more impulsive, narrower, signals in the wave train. The resulting restored time-series facilitates more accurate and objective relative traveltime estimation of the individual seismic arrivals. We demonstrate the accuracy of the reconstruction method on synthetic seismograms generated by the reflectivity method. Clean and sharp reconstructions are obtained with real data, even for signals with relatively high noise content. Reconstructed signals are simpler, more impulsive, and narrower, which allows highlighting of some details of arrivals that are not readily apparent in raw waveforms. In particular, phases nearly coincident in time can be separately identified after processing. This is demonstrated for two seismic wave pairs used to probe deep mantle and core-mantle boundary structure: (1) the Sab and Scd arrivals, which travel above and within, respectively, a 200-300-km-thick, higher than average shear wave velocity layer at the base of the mantle, observable in the 88-92 deg epicentral distance range and (2) SKS and SPdiffKS, which are core waves with the latter having short arcs of P-wave diffraction, and are nearly identical in timing near 108-110 deg in distance. A Java/ Matlab algorithm was developed for the signal restoration, which can be downloaded from the authors web page, along with example data and synthetic seismograms.",
keywords = "Core-mantle boundary, DDeconvolution, Lower mantle, Regularization, Total variation",
author = "W. Stefan and Edward Garnero and Rosemary Renaut",
year = "2006",
month = "12",
doi = "10.1111/j.1365-246X.2006.03124.x",
language = "English (US)",
volume = "167",
pages = "1353--1362",
journal = "Geophysical Journal International",
issn = "0956-540X",
publisher = "Wiley-Blackwell",
number = "3",

}

TY - JOUR

T1 - Signal restoration through deconvolution applied to deep mantle seismic probes

AU - Stefan, W.

AU - Garnero, Edward

AU - Renaut, Rosemary

PY - 2006/12

Y1 - 2006/12

N2 - We present a method of signal restoration to improve the signal-to-noise ratio, sharpen seismic arrival onset, and act as an empirical source deconvolution of specific seismic arrivals. Observed time-series gi are modelled as a convolution of a simpler time-series fi, and an invariant point spread function (PSF) h that attempts to account for the earthquake source process. The method is used on the shear wave time window containing SKS and S, whereby using a Gaussian PSF produces more impulsive, narrower, signals in the wave train. The resulting restored time-series facilitates more accurate and objective relative traveltime estimation of the individual seismic arrivals. We demonstrate the accuracy of the reconstruction method on synthetic seismograms generated by the reflectivity method. Clean and sharp reconstructions are obtained with real data, even for signals with relatively high noise content. Reconstructed signals are simpler, more impulsive, and narrower, which allows highlighting of some details of arrivals that are not readily apparent in raw waveforms. In particular, phases nearly coincident in time can be separately identified after processing. This is demonstrated for two seismic wave pairs used to probe deep mantle and core-mantle boundary structure: (1) the Sab and Scd arrivals, which travel above and within, respectively, a 200-300-km-thick, higher than average shear wave velocity layer at the base of the mantle, observable in the 88-92 deg epicentral distance range and (2) SKS and SPdiffKS, which are core waves with the latter having short arcs of P-wave diffraction, and are nearly identical in timing near 108-110 deg in distance. A Java/ Matlab algorithm was developed for the signal restoration, which can be downloaded from the authors web page, along with example data and synthetic seismograms.

AB - We present a method of signal restoration to improve the signal-to-noise ratio, sharpen seismic arrival onset, and act as an empirical source deconvolution of specific seismic arrivals. Observed time-series gi are modelled as a convolution of a simpler time-series fi, and an invariant point spread function (PSF) h that attempts to account for the earthquake source process. The method is used on the shear wave time window containing SKS and S, whereby using a Gaussian PSF produces more impulsive, narrower, signals in the wave train. The resulting restored time-series facilitates more accurate and objective relative traveltime estimation of the individual seismic arrivals. We demonstrate the accuracy of the reconstruction method on synthetic seismograms generated by the reflectivity method. Clean and sharp reconstructions are obtained with real data, even for signals with relatively high noise content. Reconstructed signals are simpler, more impulsive, and narrower, which allows highlighting of some details of arrivals that are not readily apparent in raw waveforms. In particular, phases nearly coincident in time can be separately identified after processing. This is demonstrated for two seismic wave pairs used to probe deep mantle and core-mantle boundary structure: (1) the Sab and Scd arrivals, which travel above and within, respectively, a 200-300-km-thick, higher than average shear wave velocity layer at the base of the mantle, observable in the 88-92 deg epicentral distance range and (2) SKS and SPdiffKS, which are core waves with the latter having short arcs of P-wave diffraction, and are nearly identical in timing near 108-110 deg in distance. A Java/ Matlab algorithm was developed for the signal restoration, which can be downloaded from the authors web page, along with example data and synthetic seismograms.

KW - Core-mantle boundary

KW - DDeconvolution

KW - Lower mantle

KW - Regularization

KW - Total variation

UR - http://www.scopus.com/inward/record.url?scp=33845325329&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=33845325329&partnerID=8YFLogxK

U2 - 10.1111/j.1365-246X.2006.03124.x

DO - 10.1111/j.1365-246X.2006.03124.x

M3 - Article

AN - SCOPUS:33845325329

VL - 167

SP - 1353

EP - 1362

JO - Geophysical Journal International

JF - Geophysical Journal International

SN - 0956-540X

IS - 3

ER -