the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 27 Nov 2025
Multi-component reactive transport in near-saturated deformable porous media
Abstract. This study develops a hydro–mechanical–chemical (HMC) framework for simulating reactive solute migration in near-saturated, deformable porous media. The model couples the pore-water mass balance, force equilibrium, and advection–dispersion equations, and further incorporates a flexible geochemical reaction module to address both single-reaction and multi-component, multi-mineral systems. Numerical results indicate that deformation, mechanical loading, saturation, and mineral reactions jointly control the distribution and evolution of the solute. Compression and stronger mechanical loads accelerate solute transport in the early stage but later hinder migration as the pore structure tightens. Moreover, reduced saturation promotes concentration build-up by enhancing advective transport and limiting the ability of the aqueous phase to dilute accumulated solutes. The framework improves the predictive capability for long-term plume behaviour and mineral alteration in reactive porous systems where mechanical, hydraulic, and geochemical processes interact.
- Preprint
(2103 KB) - Metadata XML
- BibTeX
- EndNote
Status: open (until 10 Feb 2026)
- CC1: 'Comment on egusphere-2025-5497', Giacomo Medici, 19 Dec 2025 reply
-
CC2: 'Scope and validity of the near-saturated formulation need to be pinned down', Hannah Menke, 28 Dec 2025
reply
Right now, near-saturated risks sounding like a convenience assumption rather than a controlled approximation. You explore S_r 0.80–1.0, but (as written) saturation appears prescribed rather than evolving, which is a big modeling decision for landfill liners and many barrier systems. Please add a short model scope paragraph early (end of Introduction or start of Methods) stating explicitly: (1) saturation is constant in time (if that’s the case), (2) when this is reasonable (e.g., quasi-saturated low gas mobility, slow infiltration variability relative to consolidation time, etc.), (3) when it is not (drying fronts, strong unsaturated flow, gas phase transport). Consider adding a brief comparison to Richards/unsaturated THMC literature (even if you don’t implement it). And be explicit whether near-saturated here means “two-phase but simplified compressibility” vs “single-phase with effective compressibility.”
Citation: https://doi.org/10.5194/egusphere-2025-5497-CC2 -
CC3: 'Validation is currently insufficient for a framework paper', Hannah Menke, 28 Dec 2025
reply
You have two verification benchmarks, but the paper itself acknowledges no lab/field quantitative validation. While this isnt a fatal flaw, you do need to strengthen credibility with additional numerical evidence. Consider adding (1) Grid/time-step sensitivity for one representative case (even a coarse/medium/fine study + one plot at the monitoring point). (2) Splitting / sub-stepping sensitivity: show that results don’t materially change when the chemistry sub-step size changes within a window (or quantify the trade-off). Your own text highlights synchronisation issues, so you should demonstrate control. (3) For benchmark #2, quantify mismatch: e.g., L2 error vs time, and show it decreases with refinement or explain the irreducible discrepancy (boundary formulation differences are mentioned but not demonstrated).
Citation: https://doi.org/10.5194/egusphere-2025-5497-CC3 -
CC4: 'More information needed for reproducibility', Hannah Menke, 28 Dec 2025
reply
At minimum, include: (1) The PHREEQC database used, the key input blocks (KINETICS/EQUILIBRIUM_PHASES/etc.), and where they are provided (supplement/Git repo). (2) OpenFOAM discretisation choices for transport (schemes, limiters), linear solvers/preconditioners, coupling tolerances, and convergence criteria (you mention a prescribed tolerance but not what it is). (3) Use a public repository plus archived case files. Without this, the ChemWindow contribution especially will be hard to credit.
Citation: https://doi.org/10.5194/egusphere-2025-5497-CC4 -
CC5: 'More information needed for reproducibility', Hannah Menke, 28 Dec 2025
reply
At minimum, include: (1) The PHREEQC database used, the key input blocks (KINETICS/EQUILIBRIUM_PHASES/etc.), and where they are provided (supplement/Git repo). (2) OpenFOAM discretisation choices for transport (schemes, limiters), linear solvers/preconditioners, coupling tolerances, and convergence criteria (you mention a prescribed tolerance but not what it is). (3) Use a public repository plus archived case files. Without this, the ChemWindow contribution especially will be hard to credit.
Citation: https://doi.org/10.5194/egusphere-2025-5497-CC5 -
CC6: 'Internal consistency of the transport and chemistry modelling choices', Hannah Menke, 28 Dec 2025
reply
I have some concerns about hybrid decisions. (1) You use element-based ADE transport, with a weighted molecular diffusion coefficient built from prescribed species fractions, while noting that fully coupled PhreeqcRM could update speciation and thereby D_m dynamically, but you don't do that in this study. (2) You also include a linear adsorption K_d term in the ADE, while stating that more detailed sorption could be handled by PHREEQC reaction terms. Please add a why this choice paragraph: is the goal robustness and speed? If so, say that explicitly and acknowledge the trade-off. Specify whether K_d is used in the case studies and, if yes, whether PHREEQC includes any overlapping sorption reactions.
Citation: https://doi.org/10.5194/egusphere-2025-5497-CC6 -
CC7: 'The case studies are plausible, but are they physics appropriate?', Hannah Menke, 28 Dec 2025
reply
You model landfill-liner-like consolidation with chemistry. There are some questions with (1) the mechanical model (elastic vs elasto-plastic/creep typical in clays), (2) the lack of reaction-induced porosity/permeability evolution (I know you already acknowledge this which is good, but it reduces interpretability of long-term plume changes), (3) boundary conditions (drainage, injection as fluxes at multiple points) and whether they represent a realistic exposure. Consider adding a short justification for the constitutive law used and explicitly call the case studies “demonstration problems” rather than predictive landfill designs. Add one sentence in the Conclusions about what would change when property feedbacks are included (directional expectation is fine, but label it as expectation).
Also some minor comments below:
(1) In benchmark #2 you say early-time deviations come from transport formulation and boundary representations; please add one sentence indicating which boundary treatment differs (Dirichlet/Neumann/mixed; inlet formulation), and whether refinement improves it.
(2) Consistently distinguish “hydraulic conductivity K” vs “permeability k” (the symbol list emphasises K).
(3) Injection strengths are given as fluxes at J1–J3; ensure it’s always explicit whether that is imposed as boundary flux (mol/(m²·s)) and over what area.
(4) Computational performance section: helpful, but reviewers will ask for scaling beyond 4 MPI ranks and/or mesh size dependence. Even a small strong-scaling plot (1,2,4,8,16) on one case would help.
(5) Consider moving the limitations to earlier in the discussion.
Citation: https://doi.org/10.5194/egusphere-2025-5497-CC7 -
RC1: 'Reproducibility', Anonymous Referee #1, 06 Jan 2026
reply
E.g. validation : Line 280 - information regarding mesh, parameters, boundary conditions and similar are not reported. This circmstance makes the reproducibility of the results a major concern.
Citation: https://doi.org/10.5194/egusphere-2025-5497-RC1 -
RC2: 'Use of saturation', Anonymous Referee #1, 06 Jan 2026
reply
please better clarify how saturation can remain constant upon loading and compaction.
why wasn't the water content considered instead?
Citation: https://doi.org/10.5194/egusphere-2025-5497-RC2
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 194 | 102 | 39 | 335 | 13 | 13 |
- HTML: 194
- PDF: 102
- XML: 39
- Total: 335
- BibTeX: 13
- EndNote: 13
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
General comments
Good research on contaminant transport in porous media that needs some improvement and further detail. See my specific comments, they will improve the potential new version of the manuscript.
Specific comments
Line 23….“has been the most widely applied framework in the field”. Last words not backed up by references. Add these papers that incorporate discussion on evidence of the Fick’s law in the field:
- Agbotui, P.Y., Firouzbehi, F., Medici G. 2025. Review of effective porosity in sandstone aquifers: insights for representation of contaminant transport. Sustainability 17, no. 14 (2025): 6469.
- Parker, B.L., Cherry, J.A., Wanner, P. 2022. Determining effective diffusion coefficients of chlorohydrocarbons in natural clays: unique results from highly resolved controlled release field experiments. Journal of Contaminant Hydrology 250, 104075.
Line 30. This sentence does not work as it is. Before “purely physical transport” and then you involve geochemistry. Please, revise the structure.
Line 76. You need to disclose the general aim of the research at the end of the introduction.
Line 76. You need to describe the specific objectives of your research by using numbers (e.g., i, ii, and iii) at the end of the introduction.
Line 199. Insert reference to Parkhurst for PHREEQC.
Line 250. This one looks to me a figure, not a table or an equation. Am I correct?
Line 700. More detail on the acknowledgement. Funding bodies are unclear.
Figures and tables
Figure 3. Make the image larger.
Figure 5. Same here, make the figure larger.
Two figures 11. This is not ok.
Figure 9. Make the figure larger also here.