SeismoStats: A Python Package for Statistical Seismology

TL;DR

SeismoStats is a comprehensive, open-source Python package designed for statistical seismology, facilitating earthquake data analysis, visualization, and estimation of key seismic parameters, thus supporting community-driven development.

Contribution

It is the first well-tested, well-documented Python package offering a wide range of statistical seismology tools for the community.

Findings

Provides tools for downloading and manipulating earthquake catalogs.

Includes functionalities for estimating Gutenberg-Richter law parameters.

Enables estimation of earthquake catalog completeness.

Abstract

We introduce SeismoStats, a Python package that enables essential statistical seismology analyses, with a focus on well-established methods. The package provides user-friendly tools to download and manipulate earthquake catalogs, but also plotting functionalities to visualize them, as well as means to perform analyses such as estimating the a- and b-value of the Gutenberg-Richter law, or estimating the magnitude of completeness of any earthquake catalog. This is the first well-tested, well-documented, and openly accessible Python package with all these features. It is intended to serve as the nucleus of a long-term community effort, continually expanding in functionality through shared contributions. We invite seismologists and developers to contribute ideas and code to support and shape its future development.