the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 27 Aug 2025
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.
- Preprint
(4283 KB) - Metadata XML
- BibTeX
- EndNote
Status: open (until 27 Nov 2025)
- RC1: 'Comment on egusphere-2025-2802', Anonymous Referee #1, 21 Oct 2025 reply
-
RC2: 'Comment on egusphere-2025-2802', Roberto Cabieces, 13 Nov 2025
reply
Dear Editor,
I have carefully reviewed the manuscript entitled “SEITCOM-1D: An Interactive 1D Code for Temperature and Composition Modelling of the Crust and Mantle from Seismological Data” by Arnaiz-Rodríguez and Fullea. The authors present a valuable open-source Julia implementation (SEITCOM-1D and SWD.jl) that effectively integrates petrological phase equilibria with forward seismological modelling of surface waves and receiver functions. The work is both timely and significant, as it provides a coherent framework that links seismological observables with thermochemical interpretation. The manuscript also includes appropriate validation tests and well-chosen illustrative examples.
Recommendation: Minor revisions. The manuscript is generally sound and well structured; however, I request several clarifications and minor additions (detailed below) prior to publication.
In summary, I consider this manuscript a meaningful contribution to the field and recommend its acceptance following the requested minor revisions.
Sincerely, Roberto Cabieces
QUESTIONS FOR AUTHORS:
1. The SWD.jl dispersion engine is presented as a fast, approximate alternative to Mineos.jl, with error behaviour described (e.g. larger errors for strong attenuation and high overtones). Could the authors provide: (a) runtime benchmarks (typical CPU/time for representative models) and (b) recommended decision rules for users when to prefer SWD.jl vs Mineos.jl?
2. The model assumes isotropy and a 1-D thermal parametrization. Could the authors discuss more explicitly the implications and limits of the isotropy assumption on (a) synthetic RF and (b) surface-wave sensitivity — and quantify (if possible) how large radial anisotropy would have to be to bias their inferred thermochemical interpretations?
3. The chemical parameterization in the mantle uses Al₂O₃ and FeO as main free variables and derives other oxides from statistical correlations. Please justify the sensitivity of the main results to these chosen parameterizations (e.g., would different correlations or a different independent oxide change the key conclusions?) and comment on applicability to very non-peridotitic domains.
4. Porosity corrections are applied for sediments but not for igneous crustal rocks (discussion and choices are stated in the text). Could the authors (a) justify this decision with references/measurements for typical secondary porosity ranges, and (b) discuss possible biases for fractured upper crustal settings (e.g. faulted regions)?
5. Validation: the manuscript shows comparisons with other codes and reports small discrepancies for selected tests. Could the authors (a) provide the exact model files / scripts used for those comparisons (so readers can reproduce Fig.5 and Fig.6), and (b) please state whether the code repository includes automated tests.
6. Attenuation (Qs) and anelastic formulations are important for SWD performance. Could the authors clarify which anelastic parameters are user-exposed vs fixed (e.g., grain-size, activation energies), and give guidance on how to choose them for realistic mantle scenarios? If possible, include an example sensitivity plot (Qs → phase shifts).
7. Code availability and reproducibility: the Data Availability states the package is on GitHub. Please (a) add a DOI / release tag (e.g., Zenodo) for the version used in this manuscript, (b) include a short README excerpt in the manuscript or supplement showing a minimal run example (input + expected output), and (c) indicate required dependencies and approximate memory/CPU requirements.
8. It would be advisable to include illustrative pseudocode for the principal core methods implemented in the software. Additionally, I recommend providing a block diagram that outlines the overall architecture and workflow of the full code.
MINOR FIXES Typos / small writing edits & places to improve.
Below I list specific phrases I found that should be corrected or improved.
These are quick fixes that would tidy the manuscript. “compressional and sheave waves velocity anomalies” — sheave → shear.
“between 0 and 4 GaP every 0.05 GPa” — likely 4 GPa (unit typo: capital P and a stray ‘a’). Please correct to GPa. “K-feldespars” — should be K-feldspars (or “K-feldspar” depending on intended form). “(NO MELT, dotter lines).” — dotter → dotted. “(màpé)” and similar accented forms appear many times — ensure consistency and define the term (you use “màpé” as shorthand; prefer to define once and then use “mapé” or “mape” consistently). “ferrocerase” (or similar) — this looks like a misspelling in the description of D’’ / lower-mantle phases (I think the intended mineral is ferropericlase or ferroperovskite — check and correct).
Reference / author name inconsistency: the text cites “Haker et al., 2003” in one place and “Hacker et al., 2004” elsewhere — unify/correct the spelling and years for that source. A few sentences are long and could be split for readability — e.g., the paragraph introducing the code’s aims would benefit from bulleting or shorter sentences (improves clarity for broad readership).
Figure captions: some captions could be more explicit about the parameter choices (e.g., state the exact composition and temperature fields used in each panel) — this helps readers reproduce the figures. Units & notation: ensure spacing is consistent around symbols (e.g., “T = 50 s” vs “T=50 s”), and check notation for subscripts/superscripts (some glyphs in the PDF look misaligned).
Replace the Data Availability Statement by:
“The full SEITCOM-1D v1.0 package is openly available for download at the project’s GitHub repository (https://github.com/marianoarnaiz/SeitComp.jl). The software requires a functional installation of the Julia Language (version 1.7 or higher) to run.”
Citation: https://doi.org/10.5194/egusphere-2025-2802-RC2
Model code and software
SeitComp Mariano Arnaiz and Javier Fullea https://github.com/marianoarnaiz/SeitComp.jl
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 1,267 | 46 | 18 | 1,331 | 74 | 65 |
- HTML: 1,267
- PDF: 46
- XML: 18
- Total: 1,331
- BibTeX: 74
- EndNote: 65
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
The Authors provided an excellent manuscript in terms of clarity, organization and quality of the content. I admit I've enjoyed the reading and have been impressed by the work done to assembly the SEITCOM-1D code. It is exactly the type of comprehensive, flexible, fast and reliable code that was missing for the forward computation of surface waves dispersion curves and RF. I'd like to congratulate the authors for this.
I have hardly any issue to mention, except for some minor suggestions:
- The name SEITCOM-1D does not seem to be an acronym, or at least this is not reported in the text. Consider to find a more appealing name, which sounds more related to "geoscience". "SEITCOM" sounds more related to the "telecomunication" field.
- L140: 1300°C --> convert to K for consistency with the rest of the text.
- L274: explain what mantle fertility means, as not all interested readers might be familiar.
- L428: I assume that dc refers to phase velocities and dU to group velocities. Please state.
- Figure 5: missing x-label for periods
- ALL FIGURES: please make all fonts/size the same across labels. Now they seem changing from one another (e.g. in Fig8-9 the "Temperature" labels are in bold...