The <i>Multitaper</i> Spectrum Analysis Package in Python
G. A. Prieto
Abstract
Abstract Spectral analysis has been a fundamental tool in analyzing seismic signals for studying the earthquake source, propagation of seismic waveforms through the Earth, and even monitoring changes in Earth’s structure. I present an open-source Python package, multitaper, for spectral analysis using the multitaper algorithm. The package not only includes power spectral density estimation (with confidence intervals) but also includes bivariate problems such as coherence, dual-frequency correlations, and deconvolution estimation. Implementation of the sine and quadratic multitaper algorithms is also available. For the reader to quickly learn how to use the package, I briefly present several examples using earthquake records from the 2019 Mw 6.0 Mesetas, Colombia, earthquake and its aftershocks recorded at regional distances for estimating time–frequency spectrograms, spectral ratio and source time functions, and correlations between neighboring frequencies. Jupyter Notebooks are shared to reproduce the figures.