Initialization of coupled thermo-hydro-mechanical models of permafrost terrain using the frozen-ground-fem package
Bibliographic record
Abstract
Permafrost degradation under climate change poses growing risks to infrastructure and environmental stability in cold regions. Processes such as frost heave, thaw settlement, and thermokarst are intensifying due to rising temperatures, changing precipitation, and evolving snow and surface hydrology. While many thermal models capture heat transfer in permafrost, they often overlook the coupled thermal, hydraulic, and mechanical (THM) interactions that are crucial for predicting thaw-induced deformation and consolidation. A critical yet underreported aspect of permafrost modelling is the “spin-up” phase: an initialization procedure that applies cyclic boundary forcing to stabilize subsurface conditions before transient simulations. Prior studies, including Ross et al. (2022) , demonstrated that thermal spin-up alone can require over 10,000 cycles for permafrost ground, but neglected hydromechanical effects such as void ratio evolution and stress redistribution. In this study, we extend Ross et al.'s approach by implementing a coupled THM spin-up strategy using the open-source Python 3 package frozen-ground-fem. This one-dimensional, large-strain finite element framework solves for both temperature and void ratio, incorporating phase change, cryosuction, nonlinear constitutive laws, and residual stress initialization. Examples using Ross et al.'s input data illustrate how coupled initialization improves stability and accuracy of permafrost predictions. Our results demonstrate that joint thermal and hydromechanical spin-up significantly alters equilibrium profiles and enhances simulation reliability. The frozen-ground-fem model is freely available on GitHub, supporting open and reproducible development of climate-resilient permafrost infrastructure models. • Thermo-hydro-mechanical spin-up using the open-source frozen-ground-fem package. • Coupled temperature and void ratio profiles establish initial ground conditions. • Sub-Arctic and Arctic sites show distinct thermal and mechanical ground conditions. • Reproducible initialization method for robust THM modelling of permafrost terrain. • Coupled spin-up enhances model stability for long-term permafrost simulations.
Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.
How this classification was reachedexpand
Full frame distilled prediction
Teacher imitationNot calibrated prevalence, not ground truth. Human validation pending. Learned from the 10,348 direct Codex labels and 10,348 direct Gemma labels. Candidate is the union of thresholded teacher heads; consensus is their intersection. These outputs are machine_predicted_unvalidated and are not human labels or direct frontier model labels.
Codex and Gemma teacher scores by category
| Category | Codex | Gemma |
|---|---|---|
| Metaresearch | 0.000 | 0.000 |
| Meta-epidemiology (narrow) | 0.000 | 0.000 |
| Meta-epidemiology (broad) | 0.000 | 0.000 |
| Bibliometrics | 0.000 | 0.002 |
| Science and technology studies | 0.000 | 0.002 |
| Scholarly communication | 0.000 | 0.000 |
| Open science | 0.000 | 0.000 |
| Research integrity | 0.000 | 0.000 |
| Insufficient payload (model declined to judge) | 0.000 | 0.000 |
Machine scores (provisional)
The two teacher heads of the student model, read on this work. A score orders the frame for review; it never asserts a category, and the validation status ships verbatim with every row.
Baseline scores from an immature model (maturity gate not passed, 7 training rounds). Scores rank; they never assert a category.
score_only:v0-immature-baseline · verbatim from the scoring run: score_only means the number may rank works, and no category label ships from itClassification
machine, unvalidatedMachine predicted; a candidate call from one teacher head, not a consensus.
How this classification was reached, model by model and score by score, is at the end of the page under "How this classification was reached".