Abstract:The tomographic inversion of shear wave splitting data for upper mantle anisotropy has been a longstanding challenge. This is due to the ray-based approximation of classical approaches and the near-vertical incidence of the core-mantle converted phases such as SKS that are often used. Recent developments include the calculation of finite-frequency sensitivity kernels for SKS splitting intensity observations, which allows us to accurately take into account the sensitivity to anisotropic structure with depth. A requirement of this tomographic technique is a dense station spacing, which results in overlapping sensitivity kernels at depth and allows for the localization of anisotropic structure. This is satisfied by a growing number of temporary seismic deployments, which motivates the desire to image anisotropic complexities with depth. Here, we introduce and make available a toolbox for the MATLAB environment that facilitates the application of finite-frequency splitting intensity tomography to dense seismic arrays. Our implementation includes several key features, including: 1) A forward calculation of splitting intensities and sensitivity kernels for a complex anisotropic model space. 2) Consideration of the dominant period of the wave, allowing for multiple-frequency analysis, as well as the incoming wave's non-vertical incidence. 3) The inversion can be based on a classical gradient descent, on a form of the conjugate gradient method known as the BFGS algorithm, or on a gradient-informed stochastic reversible jump algorithm, allowing for a data-driven parametrization of the model space. 4) Importing splitting intensity measurements from waveforms processed in SplitRacer allows for fast pre-processing of large data sets due to its fully automatic design. To illustrate our method, we present both synthetic tests and an application to real data. We apply our inversion procedure to data from the Swath-D network, which densely covers the transition of the Central to the Eastern Alps. Previous studies showed evidence for an abrupt lateral change of layered seismic anisotropy that had been attributed to an opening for channeled asthenospheric flow. Using an SKS splitting intensity tomography approach, we can confirm previous inferences while providing additional constraints on the distribution of anisotropy laterally and with depth.

PyRefra – Refraction seismic data treatment and inversion

Refrapy: A Python program for seismic refraction data analysis

DiffraPy: An open-source Python software for seismic diffraction imaging

PMTec: A New MATLAB Toolbox for Absolute Plate Motion Reconstructions from Paleomagnetism

Reflected wave least squares reverse time migration with angle illumination compensation

Pyseistr: A Python Package for Structural Denoising and Interpolation of Multichannel Seismic Data

SEISMIC , a Python-based code of the Quantas package to calculate the phase and group acoustic velocities in crystals

surfQuake: A new Python toolbox for the workflow process of seismic sources

Downward continuation of marine seismic reflection data: An undervalued tool to improve velocity models

EMagPy: Open-source standalone software for processing, forward modeling and inversion of electromagnetic induction data

TomograPy: A Fast, Instrument-Independent, Solar Tomography Software

Refraction waves full waveform inversion of deep reflection seismic profiles in the central part of Lhasa Terrane

A Compact Program for 3D Passive Seismic Source-Location Imaging

ObsPy: a bridge for seismology into the scientific Python ecosystem

Rapid Development of Interferometric Software Using MIRIAD and Python

SItomo – A toolbox for Splitting Intensity Tomography and Application in the Eastern Alps

Refractive COLMAP: Refractive Structure-from-Motion Revisited

PyIRI: Whole‐Globe Approach to the International Reference Ionosphere Modeling Implemented in Python

SAIPy: A Python Package for single station Earthquake Monitoring using Deep Learning

Open source software for simulations and inversions of airborne electromagnetic data

A multipurpose numerical method for imaging studies and tomographic reconstruction