SEITCOM-1D: An Interactive 1D Code for temperature and composition modelling of the crust and mantle from seismological data

Abstract. We present SEITCOM-1D, a software to model the Earth’s thermochemical and geophysical structure from the surface down to the core-mantle boundary (CMB). The code is designed to estimate geophysical parameters of the Earth’s crust and mantle from petrological and thermal information within a thermodynamically consistent framework and to perform forward 1D coupled geophysical-petrological modelling of the structure of the Earth. Developed in Julia Language, the open-source code is intended to be an easy-to-use, flexible, and fast. SEITCOM-1D includes tools to exploit the large repertoire of 1D seismological data available, namely: surface wave dispersion curves (of fundamental and higher modes of Rayleigh and Love waves) and receiver functions (of P, S, and SKS waves). Surface heat flow and isostatic topography can also be modelled. Four simple examples that illustrate the capabilities of the code are presented to show the sensitivity of Rayleigh wave phase velocity curves and P-to-S receiver functions to compositional and temperature variations.