the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 31 Mar 2026
SeidarT: An open-source full-waveform seismic and electromagnetic wave propagation modeling toolbox demonstrated in snow and ice
Abstract. We present SeidarT, an open-source, community-driven software package for full-waveform finite-difference time-domain (FDTD) modeling of elastic and electromagnetic wave propagation in heterogeneous, anisotropic, and attenuating media. SeidarT natively incorporates full-tensor anisotropy (all 21 stiffness coefficients), frequency-independent attenuation through a generalized Q formulation, and unified treatment of both seismic and electromagnetic wave physics on a simple Cartesian grid. The software prioritizes accessibility and extensibility by combining the computational efficiency of FORTRAN with the user-friendly scripting capabilities of Python. Model construction leverages an intuitive image-based geometry workflow and JSON project files, eliminating the need for complex mesh generation while allowing flexible specification of arbitrary stiffness or permittivity tensors, material distributions, and boundary conditions. We implement the Convolutional Perfectly Matched Layer (CPML) with explicit tuning strategies adapted for anisotropic media, and provide automatic stability checking via the Courant-Friedrichs-Lévy and wavenumber-bandlimit criteria. The software is validated against analytical solutions for elastic wave propagation and empirically constrained through comparison with ground-penetrating radar and seismic field observations in snow and ice. We document the physical property parameterizations for ice and snow as functions of temperature, pressure, and liquid water content, and provide multiple material homogenization schemes (Hill average, Gassmann substitution, Self-Consistent Approximation) to accommodate variable porosity and fluid saturation regimes. SeidarT is designed to lower economic and technical barriers for scientists, engineers, and students integrating sophisticated wave-physics simulations into workflows spanning cryospheric research, environmental monitoring, subsurface characterization, and civil infrastructure assessment. The open-source development model on GitHub and PyPI encourages community contributions and iterative improvements, positioning SeidarT as a versatile platform for advancing both fundamental understanding and applied geophysical imaging.
- Preprint
(7476 KB) - Metadata XML
- BibTeX
- EndNote
Status: final response (author comments only)
-
RC1: 'Comment on egusphere-2026-1357', Henry Moore, 26 Apr 2026
-
AC1: 'Reply on RC1', Steven Bernsen, 30 Apr 2026
Hi Henry,
Thank you for the constructive feedback and input. I do not have any questions at this time, and I will address your comments.
Steve
Citation: https://doi.org/10.5194/egusphere-2026-1357-AC1
-
AC1: 'Reply on RC1', Steven Bernsen, 30 Apr 2026
-
RC2: 'Comment on egusphere-2026-1357', Matthias Steiner, 29 Apr 2026
This manuscript presents SeidarT, a novel toolbox for modeling seismic and electromagnetic wave propagation in cryospheric settings. Given its focus on ice and snow scenarios, it is of particular interest to the readership of The Cryosphere. The open-source nature of the project and the sufficiently documented Git repository facilitate community contribution and ease of use.
SeidarT’s workflow, specifically the use of JSON files for parameterization and PNG files for domain definition, is a significant step toward lowering the barrier to entry for numerical modeling. This approach holds great potential for research, but also opens new perspectives in teaching and for citizen science projects.
The manuscript provides a thorough description of SeidarT’s computational framework. This allows potential users to understand exactly how the toolbox handles forward modeling, which (petro-)physical relationships are considered for both seismic and electromagnetic waves, and how the authors have balanced computational efficiency with numerical accuracy.
While SeidarT and the accompanying manuscript are of utmost relevance, the presentation of the toolbox needs to be enhanced prior to publication in The Cryosphere:
Manuscript Preparation
- The manuscript contains numerous citation errors, typos, unreferenced figures, and incomplete sentences that hinder readability.
- The manuscript is currently unbalanced. While the computational framework is described at length, the Validation section is disproportionately brief. An in-depth description of the numerical models and a more robust discussion of results versus theoretical expectations are necessary.
Figures
- Figures lack subplot labels, which makes description provided in the main text of the manuscript hard to follow.
- Some figures seem to be unfinished due to the lack of axis labels, legends etc.
- Several figures are not sufficiently described or interpreted within the main body of the text.
Numerical Studies
- The settings for the synthetic experiments are not sufficiently detailed. These should be supported by graphical illustrations of the experimental design (e.g., layer geometry and source/receiver configurations).
- Information regarding model parameters is not always provided with the required detail or scattered throughout the text. This should be consolidated at the beginning of the relevant sections.
- The authors should explicitly describe the expected physical outcomes of their experiments rather than assuming the reader is familiar with the underlying theoretical considerations or cited literature.
Detailed comments, additional questions and suggestions are provided in the attached commented PDF version of the manuscript.
I hope the review and comments serve as constructive input that is useful to the authors for preparing the manuscript for publication in The Cryosphere.
-
AC2: 'Reply on RC2', Steven Bernsen, 30 Apr 2026
Hi Matthias,
Thank you for the constructive feedback and input. I do not have any questions at this time, and I will address your comments.
Steve
Citation: https://doi.org/10.5194/egusphere-2026-1357-AC2
-
RC3: 'Comment on egusphere-2026-1357', Anonymous Referee #3, 06 May 2026
The manuscript presents SeidarT, an open-source software package for full-waveform finite-difference time-domain (FDTD) modeling of seismic and electromagnetic wave propagation in heterogeneous, anisotropic, and attenuating media. The code integrates a unified treatment of elastic and electromagnetic physics within a single framework.
The topic is relevant to the Solid Earth community, particularly given the increasing demand for accessible, flexible, and reproducible numerical tools. The manuscript is generally well-structured and includes validation against analytical solutions as well as field observations.
However, in its current form, the manuscript reads more as a technical/software description with illustrative examples than as a rigorous validation and benchmarking study. The manuscript claims “computational efficiency,” but no benchmarks are provided. There is no discussion of code performance, scalability, or resource usage. The authors should include benchmark (e.g. runtime vs model size, CPU/GPU usage, memory requirements) and ideally scaling test to support this claim.
One valuable contribution of the software is the image-based geometry workflow (PNG + JSON) for model construction. This approach significantly simplifies model setup compared to traditional mesh-based or script-defined geometries, especially for geometrically complex or highly heterogeneous media. However, the authors should also discuss the potential limitations of raster-based geometries, such as resolution dependence and staircasing effects, and their impact on numerical accuracy. In particular, the relationship between pixel size and model resolution could introduce aliasing effects. It would be very useful to provide a practical guideline, for example, in terms of pixels per wavelength.
Figure 4: In the caption, profile ZOP2 is mentioned, but this profile is not introduced or described in the text, leaving the reader without context. As ZOP2 is not defined, the figure is difficult to interpret.
Figure 8: The left panel lacks a label on the horizontal axis.
Language is generally clear, but occasionally verbose and repetitive. And some sentences are incomplete (see, for instance, lines 483-484) and should be revised for clarity.
Citation: https://doi.org/10.5194/egusphere-2026-1357-RC3
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 288 | 179 | 37 | 504 | 20 | 21 |
- HTML: 288
- PDF: 179
- XML: 37
- Total: 504
- BibTeX: 20
- EndNote: 21
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
This manuscript presents the new modeling toolbox, SeidarT, for modeling seismic and electromagnetic wave propagation in snow/ice settings. Modeling presented in this manuscript is of particular interest to readers of The Cryosphere in the combined forward modeling approach of seismic and GPR data commonly collected in snow/ice settings. A comprehensive, computationally efficient, and open-source framework is presented with the manuscript, encouraging community use of the forward model. Parameterization and model limitations are well established throughout the manuscript, while also suggesting the possibilities of future research. Results from the forward model contrasted with field data reveal reliable agreement and show promise for interested researchers. Two main comments regarding the manuscript are highlighted below for the consideration of the authors with revisions.
General Comment 1:
In the Introduction, SeidarT is presented as a forward model capable of handling anisotropy in any environment in which seismic/GPR data are collected. Upon reading the rest of the manuscript, many parameters and governing equations are tuned only to the application in snow/ice settings. The Introduction should be tailored to applications in snow/ice environments and not more generally to ensure capabilities are not “over-promised”. While it is acknowledged that such a model could be extended to other applications, the model is not parameterized for such applications, and the manuscript does not present enough case studies for broader applications. Rather than targeting these broader applications, it is suggested that authors simply make this note in the Summary/Conclusions that future work could be done to examine these applications in other fields.
General Comment 2:
Within the governing equations section, further parallels between the seismic and GPR methods could be touched on. Specifically, the attenuation of the EM wave and seismic wave could use additional explanation as to the same amplitude decay governing the propagation in the subsurface. Clarification on the use of the term “permittivity” should also be considered in the GPR section for reader understanding of either the “relative dielectric permittivity” or the “relative permittivity”.
Final Comment:
As a reviewer I see the value of this model and the fit of the manuscript to The Cryosphere journal. Additional minor grammatical and specific queries are provided in the attached “mark-up” version of the manuscript. The suggestions presented above are intended to be constructive, with the aim that they will be useful to the authors.