MOIST: a MATLAB-based fully coupled one-dimensional isotope and soil water transport model
- 1Soil Science Department, University of Saskatchewan, Saskatoon, S7N 5A8, Canada
- 2College of Resources and Environmental Engineering, Ludong University, Yantai, Shandong Province, 264025, China
- 1Soil Science Department, University of Saskatchewan, Saskatoon, S7N 5A8, Canada
- 2College of Resources and Environmental Engineering, Ludong University, Yantai, Shandong Province, 264025, China
Abstract. Modeling water stable isotope transport in soil is crucial to sharpen our understanding of water cycles in terrestrial ecosystems. However, isotope and soil water transport are not fully coupled in current models. In this study, we developed MOIST: a MATLAB-based one-dimensional isotope and soil water transport model, a program that solves one-dimensional water, heat, and isotope transport equations simultaneously. Results showed that the MOIST model has good agreements to the theoretical tests and semi-analytical solutions of isotope transport under fixed boundary conditions. Furthermore, we validated the program with short- and long-term measurements from lysimeters studies. The overall Nash-Sutcliff efficiency coefficient (NSE) of soil water and deuterium (2H) transport for the short-term measurements are 0.66 and 0.69, respectively, with respective determine coefficient (R2) of 0.82 and 0.70, mean absolute error (MAE) of 0.02 m3 m-3 and 11.84 ‰. For the long-term lysimeter study, the overall NSE, R2, and MAE of simulated δ18O are 0.47, 0.49, and 0.92 ‰, respectively. These indices indicated the excellent performance of the MOIST model in simulating water flow and isotope transport in simplified ecosystems, suggesting a great potential of our program in promoting understandings of ecohydrological processes in terrestrial ecosystems.
Han Fu et al.
Status: open (until 13 Mar 2023)
-
RC1: 'Comment on hess-2022-422', Anonymous Referee #1, 03 Feb 2023
reply
I suggest rejection for the following reasons:
Equation (1) is the water content-based Richards equation. It is well known that this kind of formulation cannot handle saturated problems and is not well posed at the interface between two layers, because water content is discontinuous. The mixed form of Richards equation should be used.
The heat transfer equation (2) is not correct. Csoil depends on the water content and should be embed in the time derivative.
More information should be provided concerning the solvers (ode113, ode23tb).
The tests 1 to 6 are very qualitatively discussed. Only different types of processes are checked. Physical processes can be verified but it does not mean that the computed variables and the process kinetics are correct. Moreover, these tests are development tests. They do not provide any new information on processes and therefore, should not be part of the manuscript. It is expected that models overcome these kinds of tests before publication.
L479-480: the reason for poor MAE value is unclear to me.
L521-540: The analysis of the difference between fully coupled or sequential approach (segregation) is convincing but it applies for an explicit time scheme discretization whereas HYDRUS and SiSPat use an implicit scheme. Moreover, the flow equation is written in terms of water content for MOIST, the other codes are using pressure based or a mixed form of Richards equation.
L608-610: The discussion about boundary conditions and intermodal conductivity is very popular. There are key papers not cited in the manuscript that review some of the techniques (see for example Belfort et al.,. On equivalent hydraulic conductivity for oscillation–free solutions of Richards equation. Journal of Hydrology, 2013, 505, pp.202-217).
MOIST was used to simulate two types of experiments and the authors concluded that MOIST is more accurate and reliable. This is not supported by the provided results. These results only show that MOIST might be better calibrated not that the numerical scheme – fully coupled- is better than other schemes. Parameters used by MOIST and the other models should be given.
The comparisons do not provide any information on the code accuracy and efficiency. To demonstrate the ‘excellent performance of the MOIST,’ the authors should compare their code with other existing codes (for example looking at breakthrough curves at different locations) and check detailed mass balances, time and space discretization sensitivity and computer time.
Han Fu et al.
Model code and software
MOIST Han Fu, Bingcheng Si https://github.com/HAN-2/MOIST
Han Fu et al.
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
261 | 85 | 6 | 352 | 2 | 4 |
- HTML: 261
- PDF: 85
- XML: 6
- Total: 352
- BibTeX: 2
- EndNote: 4
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1