Home

Awesome

OMS_FreeThaw1D

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.

FreeThaw1D 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.

Alt text

Website

OMS

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 source code

GEOframe programs run on Java 8. The OMS project can be run either within the OMS v3 Console or by using Docker.

OMS v3 Console

OMS from Docker

Please refers to the guidelines at https://hub.docker.com/r/omslab/oms/. Please use the version 0.3_py3.

Project documentation

In the folder Jupyter_Notebook` you can find

Acknowledgements