SeismoStats: A Python Package for Statistical Seismology

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.