Freezing and thawing of soils affect a wide range of biogeochemical and hydrological (Walvoord and Kurylyk, 2016; Schuur et al., 2015) processes and interact with engineered structures in cold regions. This foundational importance together with the vast extent of cold regions globally makes the simulation of freezing and thawing soil an important and well-researched topic (Streletskiy et al., 2019; Walvoord and Kurylyk, 2016; Harris et al., 2009). The interest in understanding and anticipating the consequences of permafrost thaw driven by anthropogenic climate change has added additional urgency.
Studying frozen soil can be done at different timescales. At short timescale, the latent heat of fusion of water determines a time-lag in the cooling and warming of the soil thus affecting the surface heat fluxes. Studies have shown that proper frozen soil scheme help improve land surface and climate model simulation (Viterbo et al., 1999; Smirnova et al., 2000). On the other hand, permafrost evolution is characterized by a long timescale and its evolution is manly controlled by climate change. Permafrost degradation has significant impact on the ecosystems and on the social and economic life in cold regions (Bao et al., 2016) Therefore, it is desirable to have a numerical model that can efficiently simulate frozen soil at both short timescale and long timescale. At present, a common problem of physical-based frozen soil model concerns the treatment of the nonlinear behavior of the energy and enthalpy as function of temperature, water content and other governing variables. Because of this nonlinearity, the model may fail to compute the correct solution, or a short time step must be chosen leading to an increase of the computational cost.
FreThaw1D model solves the so called enthalpy - or conservative - formulation. The novelty regards the use of the nested Newton-Casulli-Zanolli (NCZ) algorithm, (Casulli and Zanolli) to linearize the nonlinear system resulting from the approximation of the governing equation. Using the enthalpy formulation and the NCZ algorithm the convergence of the solver is guarateed for any time step size. This is a key feature since the integration time step can be choosen accordingly to the time scale of the processes to study, from seconds to days.
FreeThaw1D is written in Java, works under the OMS3(David et al., 2013) framework and is part of the GEOframe system (Formetta et al., 2014, Bancheri, 2017). It was produced as part of the Ph.D. work by Niccolò Tubini.
Here you can find the OMS project
- Niccolò Tubini, Stephan Gruber, Riccardo Rigon developed the theoretical aspects of the model (Authors).
- Niccolò Tubini, Riccardo Rigon designed the first version of the code (Authors)
- Niccolò Tubini implemented and deployed it (Authors)
- Riccardo Rigon and Stephan Gruber provided financial support
- Niccolò Tubini and Riccardo Rigon wrote the documentation in the Notebooks
- Niccolò Tubini was partially supported by a Ph.D. grant by DICAM-UniTrento and by the NSERC PermafrostNet project.
- We thank Professor Vincenzo Casulli and Professor Michael Dumbser for their fruitful discussions on the numerical aspects of the work.