Options
Surface-wave tomography using SeisLib: a Python package for multiscale seismic imaging
Author(s)
Language
English
Obiettivo Specifico
1T. Struttura della Terra
Status
Published
JCR Journal
JCR Journal
Peer review journal
Yes
Title of the book
Issue/vol(year)
/231(2022)
ISSN
0956-540X
Publisher
Oxford University Press - The Royal Astronomical Society
Pages (printed)
1011-1030
Issued date
June 2022
Abstract
To improve our understanding of the Earth’s interior, seismologists often have to deal with enormous amounts of data, requiring automatic tools for their analyses. It is the purpose of this study to present SeisLib, an open-source Python package for multiscale seismic imaging. At present, SeisLib includes routines for carrying out surface-wave tomography tasks based on seismic ambient noise and teleseismic earthquakes. We illustrate here these functionalities, both from the theoretical and algorithmic point of view and by application of our library to seismic data from North America. We first show how SeisLib retrieves surface-wave phase velocities from the ambient noise recorded at pairs of receivers, based on the zero crossings of their normalized cross-spectrum. We then present our implementation of the two-station method, to measure phase velocities from pairs of receivers approximately lying on the same great-circle path as the epicentre of distant earthquakes. We apply these methods to calcu- late dispersion curves across the conterminous United States, using continuous seismograms from the transportable component of USArray and earthquake recordings from the permanent networks. Overall, we measure 144 272 ambient-noise and 2055 earthquake-based dispersion curves, that we invert for Rayleigh-wave phase-velocity maps. To map the lateral variations in surface-wave velocity, SeisLib exploits a least-squares inversion algorithm based on ray theory. Our implementation supports both equal-area and adaptive parametrizations, with the latter al- lowing for a finer resolution in the areas characterized by high density of measurements. In the broad period range 4–100 s, the retrieved velocity maps of North America are highly correlated (on average, 96 per cent) and present very small average differences (0.14 ± 0.1 per cent) with those reported in the literature. This points to the robustness of our algorithms. We also produce a global phase-velocity map at the period of 40 s, combining our dispersion measurements with those collected at global scale in previous studies. This allows us to demonstrate the reliability and optimized computational speed of SeisLib, even in presence of very large seismic inverse problems and strong variability in the data coverage. The last part of the manuscript deals with the attenuation of Rayleigh waves, which can be estimated through SeisLib based on the seismic ambient noise recorded at dense arrays of receivers. We apply our algorithm to produce an at- tenuation map of the United States at the period of 4 s, which we find consistent with the relevant literature.
Sponsors
Funded by the Deutsche Forschungsgemeinschaft (DFG - German Research Foundation) under the Individual Research Project: SI 1748/4- 1.
German Science Foundation: Deutsche Forschungsgemeinschaft (www.dfg.de; SPP-2017, Project Ha 2403/21-1).
German Science Foundation: Deutsche Forschungsgemeinschaft (www.dfg.de; SPP-2017, Project Ha 2403/21-1).
Description
This article has been accepted for publication in Geophysical Journal International ©:The Author(s) 2022. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved.Uploaded in accordance with the publisher's self-archiving policy.
All rights reserved.
All rights reserved.
Type
article
File(s)
Loading...
Name
magrini_seislib_python.pdf
Description
Article
Size
11 MB
Format
Adobe PDF
Checksum (MD5)
e0902888804060c4322a1d0e34192dd7