REVIEW 3 major objections 5 minor 20 references
Integrated Experimental and Numerical Investigations on the Thermo-Hydro-Mechanical Behavior of Clays and Argillaceous Rocks: A Perspective
T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read This perspective claims that one family of critical-state constitutive models can unify the coupled thermal, hydraulic, and mechanical behavior of clays and claystones, from desiccation cracks in surface soils to multi-year heating around…
desk verdict A clear but heavily self-referential roadmap to one group's THM modeling program; the flagship validation claim is asserted, not shown. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is a family of critical-state based constitutive models, meaning models in which the soil's stress-strain response is organized around the critical-state line and the yield surface is coupled to temperature, suction, and structure. The family grew from a non-isothermal extension of a basic unsaturated-soil model into a unified elasto-viscoplastic framework for stiff claystones, then into versions for cyclic slip surfaces and hard-soils–weak-rocks. This machinery does the work of converting laboratory observations and field sensor data into predictions: it provides the stress-strain relations that the fully coupled THM finite-element code solves, and it supplies the calibrated parameters that allow the same model to be scaled from a small triaxial specimen to a repository gallery.
What would settle it
A direct check would take one of the calibrated models and run it cold against a second in-situ heating or excavation dataset that was not used in calibration, then compare computed pore-pressure and displacement transients with sensor measurements. A systematic mismatch in the timing or magnitude of delayed swelling, or in pore-pressure dissipation, would show that the claimed unified validation does not transfer to new conditions.
Extended reading notes
Core claim
The central claim is that the full range of thermally, hydraulically, and mechanically coupled behaviors seen in clays and argillaceous rocks can be captured by one constitutive family built on critical-state soil mechanics. The models unify four mechanisms—temperature-dependent plasticity, time-dependent creep, suction hysteresis in unsaturated behavior, and progressive loss of bonding or structure—under a single thermodynamic potential, and they are implemented in an implicit, fully coupled THM code. The paper reports that this implementation reproduced key field observations: delayed anisotropic convergence around tunnels in claystone, pore-pressure transients during multi-year in-situ heating, permeability and saturation changes in excavation-damaged zones, irreversible cracking and rebound in desiccated expansive strata, and degradation of shaft resistance in energy foundations under cyclic heating. Validation rests on comparisons with laboratory triaxial tests, wetting-drying cycles, and sensors installed in underground and surface field sites.
Load-bearing premise
The claim rests on the accuracy of the author's previously published calibration and validation work: this paper asserts that simulations reproduced key trends in pore pressure and displacement, but it shows no data or quantitative comparisons itself.
Editorial extensions
If this is right
- Repository design can rely on coupled THM simulations that reproduce anisotropic convergence and delayed swelling around disposal tunnels rather than on simplified uncoupled estimates.
- Permeability and saturation changes in excavation-damaged zones under sustained heating can be quantified, informing sealing and safety assessments for high-level waste cells.
- Critical suction thresholds for desiccation cracking can be identified, allowing slope and foundation designs to anticipate irreversible fissuring before it appears.
- Cyclic thermal loads on energy foundations can be simulated, so shaft-capacity degradation can be built into the design of geothermal piles and retaining structures.
- The Bayesian-calibrated framework can supply design charts and digital twins for site-specific repository layouts and climate-sensitive slope stability assessments.
Reading between the lines
- The paper does not say this, but the same constitutive family could be transplanted to other low-permeability geomaterials, such as shales around carbon-dioxide or hydrogen storage sites, because the coupling structure is not specific to clay mineralogy.
- The paper reports fitted comparisons rather than blind forecasts; a stronger public test would be to freeze the calibrated parameters and predict a later, unseen heating experiment.
- The Bayesian calibration workflow implies that probabilistic, not just deterministic, performance assessments are within reach, which would change how safety cases report uncertainty.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript is a perspective-style synthesis of the author's roughly decade-long research programme on the thermo-hydro-mechanical (THM) behaviour of clays and argillaceous rocks. It opens with a high-level abstract claiming that a family of critical-state based constitutive models unifies temperature-dependent elasto-plasticity, viscoplastic creep, suction hysteresis and structure degradation within a single thermodynamically consistent framework, and that these models have been validated against laboratory and in-situ data. Section 2 recounts the evolution of the constitutive models, Section 3 describes applications to Callovo-Oxfordian claystone in-situ heating tests and micro-tunnel simulations, Section 4 addresses desiccation cracking and slope processes at the Počáry test site, and Section 5 outlines future directions including chemo-THM coupling, upscaling, digital twins and machine-learning-based parameter inversion. The text contains no equations, no data tables, no error bars and no quantitative comparison with measurements; all validation claims are supported only by citations to the author's own prior publications, several of which are peer-reviewed journal articles.
Significance. If the underlying claims are accurate, the described framework could be of considerable practical significance for nuclear waste repository design, energy geostructures and climate-resilient infrastructure, and the consolidated overview may help practitioners identify relevant models and datasets. The paper's main strength is its breadth and the clear organization of a decade of related work, and the author is honest in presenting it as a perspective rather than as a new technical contribution. However, as a standalone scholarly document it is not verifiable: the central assertion of a validated, thermodynamically consistent unified model is presented without any of the mathematical or empirical evidence that would allow a reader to check it. The paper therefore currently functions more as an annotated bibliography of the author's own work than as a critical or independent synthesis.
major comments (3)
- [Abstract and §3] The abstract and Section 3 state that the THM framework 'successfully reproduced key trends' in thermal profiles, pore pressure and radial displacements, and the abstract claims the models have been 'validated against laboratory and in-situ data', but the manuscript contains no data, no quantitative comparisons, no error metrics and no figures. Since this validation claim is the central claim of the paper, the author should include at least one representative quantitative comparison, for example measured and computed temperature and pore-pressure time series from the heating test described in [13] or the convergence data in [11,12], or alternatively should explicitly state in the abstract and introduction that this is a perspective and give precise pointers to the specific figures and tables in [13], [14] and [20] that support each claimed outcome. Without such support, a reader of this manuscript alone cannot assess whether the flagship claim is warranted.
- [§2] Section 2 asserts that the family of critical-state models forms 'a single thermodynamically consistent framework', but no governing equations are provided: there is no free-energy function, dissipation potential, yield function, flow rule, hardening law or coupling term. Thermodynamic consistency is a mathematical property that must be demonstrated, not asserted. The author should either append a minimal set of constitutive equations (or state explicitly which specific equations in [4] and [7] establish the property) or temper the wording to say that thermodynamic consistency was shown in the cited works. As written, the label cannot be checked by the reader.
- [§3 and §4] The evaluation of the framework is entirely self-referential: all 'validations' involve comparisons of the author's own simulations with data that the author's group also interpreted, and no independent alternative models or external benchmarks are considered. For a perspective article this is a significant gap in scholarly balance. The author should add a comparative discussion of at least one independent THM modelling exercise for Callovo-Oxfordian claystone (for example, the international DECOVALEX-type benchmarks or the ALC1605 exercise) and should explicitly discuss any cases where the proposed model required recalibration, failed to match observations, or exhibited significant discrepancies. Including such limitations would increase the credibility of the synthesis.
minor comments (5)
- [§3] The manuscript twice writes 'HL W' with a space; this should be 'HLW'.
- [References] Reference [20] has a leading space before 'Journal of Rock Mechanics' and the author name 'Počáry' contains a rendering artifact; the reference should be cleaned up.
- [Overall structure] The paper has no conclusions or limitations section; it ends abruptly with the Outlook section. A short 'Concluding remarks' paragraph that explicitly states what the framework can and cannot currently do would improve the perspective.
- [Abstract vs. body] The abstract mentions a 'Bayesian workflow using sensor data', but the body text does not describe any Bayesian calibration or inversion procedure; this methodological component is introduced and then never discussed. Either add a brief description in Section 2 or remove the claim from the abstract.
- [§5] The terms 'digital twins', 'machine-learning pipelines' and 'reduced-order surrogates' appear in the future-work section without technical context or references; this is acceptable for a perspective, but the author should be careful to frame them as envisioned concepts rather than as existing capabilities.
Circularity Check
The central validation claim reduces to a calibration consistency check: the model is Bayesian-calibrated on sensor data, then the same sensor-measured trends are reported as 'successfully reproduced', with no holdout or quantitative comparison shown.
-
fitted input called prediction
[Abstract; Section 3, second paragraph ('In [13, 14]...')]
"The models are implemented in an implicit, fully coupled THM code and calibrated with a Bayesian workflow using sensor data from instrumented tunnels, boreholes and surface observatories. ... In [13, 14], we used a fully coupled THM framework to simulate the long-duration heating phase and its effect on pore pressure buildup, thermal expansion, and stress redistribution in COx. ... Simulations successfully reproduced key trends in thermal profiles, pore pressure evolution, and radial displacements measured by sensors embedded in the host rock."
The calibration input and the reported validation output are the same genre of data. The abstract says the models were Bayesian-calibrated using sensor data from instrumented tunnels, boreholes, and surface observatories; Section 3 then presents as its main evidence that simulations 'successfully reproduced key trends' in thermal profiles, pore pressures, and radial displacements measured by sensors in the host rock. No holdout split, withheld borehole, or out-of-sample comparison is described, and no quantitative misfit is reported.
full rationale
This paper is a perspective, not a new derivation, so most of its claims are summaries of the author's prior work. The one clear circular step is in the validation narrative: the abstract states that the models were calibrated with a Bayesian workflow using sensor data from instrumented tunnels, boreholes, and observatories, while Section 3 reports as the key evidence that simulations reproduced trends in thermal profiles, pore pressures, and displacements measured by embedded sensors. Without an explicit statement that the validation sensors were withheld from calibration, the successful reproduction is a fit-based consistency check rather than an out-of-sample prediction. I do not count the many self-citations themselves as circular, since the cited papers are peer-reviewed and are based on external field observations, but this manuscript does not show any of those comparisons, so it cannot convert them into independent evidence here. No free-energy function or dissipation inequality is displayed, so the 'thermodynamically consistent' assertion is unsupported in the text, though that is a completeness issue rather than a circularity. Overall, the paper's central validation claim partially reduces by construction to its calibration input, giving a partial circularity score of 6.
Assumptions & free parameters
free parameters (4)
- temperature-dependent permeability parameters
- anisotropic thermal conductivity coefficients
- viscoplastic creep parameters
- critical suction thresholds for cracking
assumptions (5)
- domain assumption Critical state soil mechanics is applicable to clays and claystones under thermal and hydraulic loading.
- domain assumption Thermo-hydro-mechanical coupling is necessary to reproduce observed behavior.
- domain assumption The cited numerical implementations are correct and their results are accurate.
- domain assumption Bayesian calibration using sensor data yields reliable parameter estimates.
- domain assumption The in situ sensor data are reliable and the simulations captured key trends.
Cite this review
Pith. "Pith review of Integrated Experimental and Numerical Investigations on the Thermo-Hydro-Mechanical Behavior of Clays and Argillaceous Rocks: A Perspective." pith.science (2026). https://pith.science/paper/JWCERL5O
@misc{pith2026250605818,
author = {Pith},
title = {Pith review of: Integrated Experimental and Numerical Investigations on the Thermo-Hydro-Mechanical Behavior of Clays and Argillaceous Rocks: A Perspective},
year = {2026},
howpublished = {\url{https://pith.science/paper/JWCERL5O}},
note = {Machine review of arXiv:2506.05818}
}
read the original abstract
This paper synthesizes nearly a decade of research on the coupled thermo-hydro-mechanical (THM) behavior of clays and argillaceous rocks. Drawing from experimental observations, numerical model development, and field-scale simulations, it presents a consolidated view of soil-structure interaction under thermal loading, desiccation cracking, and long-term excavation impacts. Key findings are drawn from constitutive modeling, in situ tests, and energy geostructure applications, offering a practical THM framework for nuclear waste repositories and climate-resilient infrastructure.
Reference graph
Works this paper leans on
-
[13]
S. Tourchi, J. Vaunat, A. Gens, F. Bumbieler, M. Vu, and G. Armand, “A full- scale in situ heating test in callovo-oxfordian claystone: observations, analysis and interpretation,” Computers and Geotechnics , vol. 133, p. 104045, 2021
work page 2021
-
[14]
F. Bumbieler, C. Pl´ ua, S. Tourchi, M. Vu, J. Vaunat, A. Gens, and G. Armand, “Feasibility of constructing a full-scale radioactive high-level waste disposal cell and characterization of its thermo-hydro-mechanical behavior,” International Journal of Rock Mechanics and Mining Sciences , vol. 137, p. 104555, 2021
work page 2021
-
[20]
S. Tourchi, M. Jabbarzadeh, A. Alimardani Lavasan, H. Sadeghi, and O. Racek, “Thermo-hydro-mechanical dynamics of a rock slope: Integrated field and numerical analysis at the poˇ z´ ary test site in the czech republic,”Journal of Rock Mechanics and Geotechnical Engineering, 2024. 7
work page 2024
-
[4]
Thm analysis of argillaceous rocks with application to nuclear waste un- derground storage,
S. Tourchi, “Thm analysis of argillaceous rocks with application to nuclear waste un- derground storage,” Ph.D. dissertation, UPC [Department of Civil and Environmen- tal Engineering Section of Geotechnical Engineering and Geosciences], Barcelona, Spain, 2020
work page 2020
-
[7]
S. Tourchi, M. M´ anica, A. Gens, J. Vaunat, M. Vu, and G. Armand, “A thermome- chanical model for argillaceous hard soils–weak rocks: application to thm simulation of deep excavations in claystone,” G´ eotechnique, vol. 75, no. 2, pp. 139–152, 2025
work page 2025
-
[1]
Thermo-mechanical constitutive modeling of unsatu- rated clays based on the critical state concepts,
S. Tourchi and A. Hamidi, “Thermo-mechanical constitutive modeling of unsatu- rated clays based on the critical state concepts,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 7, pp. 193–198, 2015
work page 2015
-
[2]
A thermomechanical constitutive model for unsaturated clays,
A. Hamidi and S. Tourchi, “A thermomechanical constitutive model for unsaturated clays,” International Journal of Geotechnical Engineering , 2016
work page 2016
-
[3]
A critical state based thermo-elasto- plastic constitutive model for structured clays,
A. Hamidi, S. Tourchi, and F. Kardooni, “A critical state based thermo-elasto- plastic constitutive model for structured clays,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 9, pp. 1094–1103, 2017
work page 2017
Show all 20 references
-
[5]
A thermomechanical model for argillaceous rocks,
S. Tourchi, A. Gens, J. Vaunat, M. M´ anica, and G. Scaringi, “A thermomechanical model for argillaceous rocks,” E3S Web of Conferences , vol. 205, p. 13014, 2020
2020
-
[6]
A non-isothermal viscoplastic model for clay slip surface: formulation and validation,
S. Tourchi, “A non-isothermal viscoplastic model for clay slip surface: formulation and validation,” Submitted to Computers and Geotechnics , 2025
2025
-
[8]
Numerical modeling of volumetric behavior of unsaturated expansive soil under wetting and drying cycles,
M. Ghandilou, S. Tourchi, A. G. Darzi, and H. Sadeghi, “Numerical modeling of volumetric behavior of unsaturated expansive soil under wetting and drying cycles,” in 5th Iranian Conference on Geotechnical Engineering , 2023. 6
2023
-
[9]
Numerical investigation of cyclic wetting and drying of boom clay based on the barcelona expansive model,
M. J. Ghandilou, S. Tourchi, and H. Sadeghi, “Numerical investigation of cyclic wetting and drying of boom clay based on the barcelona expansive model,” in 84th EAGE Annual Conference & Exhibition , 2023
2023
-
[10]
Thermomechanical behaviour of silty sandy clays: An experimental and numerical investigation,
H. Hoseinimighani, S. Tourchi, and J. Szendefy, “Thermomechanical behaviour of silty sandy clays: An experimental and numerical investigation,” in 10th European Conference on Numerical Methods in Geotechnical Engineering , 2023
2023
-
[11]
Coupled thm analysis of long-term anisotropic convergence in the full-scale micro tunnel excavated in the callovo-oxfordian argillite,
S. Tourchi, J. Vaunat, A. Gens, M. Vu, and F. Bumbieler, “Coupled thm analysis of long-term anisotropic convergence in the full-scale micro tunnel excavated in the callovo-oxfordian argillite,” COUPLED PROBLEMS 2019 , pp. 291–299, 2019
2019
-
[12]
Thermo-hydro-mechanical simulation of a full-scale steel-lined micro-tunnel excavated in the callovo-oxfordian argillite,
——, “Thermo-hydro-mechanical simulation of a full-scale steel-lined micro-tunnel excavated in the callovo-oxfordian argillite,” in COMPLAS 2019, 2019, pp. 544–552
2019
-
[15]
Thermo-hydro-mechanical simu- lation of deep excavations in claystone,
S. Tourchi, M. M´ anica, A. Gens, and J. Vaunat, “Thermo-hydro-mechanical simu- lation of deep excavations in claystone,” in 10th European Conference on Numerical Methods in Geotechnical Engineering, 2023
2023
-
[16]
Non-isothermal behavior of excavation dam- aged zone around deep radioactive waste disposal,
S. Tourchi, A. Lavasan, and A. Gens, “Non-isothermal behavior of excavation dam- aged zone around deep radioactive waste disposal,” in Clay Conference 2024 , 2024
2024
-
[17]
Thermo-hydro-mechanical analysis of soil strata suffering from desiccation cracking,
S. Tourchi, M. Jabbarzadeh, and H. Sadeghi, “Thermo-hydro-mechanical analysis of soil strata suffering from desiccation cracking,” in XVIII European Conference on Soil Mechanics and Geotechnical Engineering , 2024, pp. 1362–1367
2024
-
[18]
Thermo-hydraulic analysis of desiccation cracked soil strata considering ground temperature and mois- ture dynamics under the influence of soil-atmosphere interactions,
M. Jabbarzadeh, H. Sadeghi, S. Tourchi, and A. Golaghaei Darzi, “Thermo-hydraulic analysis of desiccation cracked soil strata considering ground temperature and mois- ture dynamics under the influence of soil-atmosphere interactions,” Geomechanics for Energy and the Environmen...
2024
-
[19]
Thermo-hydro-mechanical mod- elling of the heterogeneous subsidence and swelling in the desiccation cracked clayey strata,
H. Sadeghi, M. Jabbarzadeh, and S. Tourchi, “Thermo-hydro-mechanical mod- elling of the heterogeneous subsidence and swelling in the desiccation cracked clayey strata,” Engineering Geology, vol. 343, p. 107798, 2024
2024
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.