Articles | Volume 15, issue 6
https://doi.org/10.5194/tc-15-2541-2021
© Author(s) 2021. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/tc-15-2541-2021
© Author(s) 2021. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Research article
| 
04 Jun 2021
Research article |
| 04 Jun 2021
| 04 Jun 2021
A method for solving heat transfer with phase change in ice or soil that allows for large time steps while guaranteeing energy conservation
Department of Civil, Environmental and Mechanical Engineering, University of Trento, Trento, Italy
Stephan Gruber
Department of Geography and Environmental Studies, Carleton University, Ottawa, ON, K1S 5B6, Canada
Riccardo Rigon
Department of Civil, Environmental and Mechanical Engineering, University of Trento, Trento, Italy
Center Agriculture Food Environment, University of Trento, Trento, Italy
Related authors
Riccardo Rigon, Giuseppe Formetta, Marialaura Bancheri, Niccolò Tubini, Concetta D'Amato, Olaf David, and Christian Massari
Hydrol. Earth Syst. Sci., 26, 4773–4800, https://doi.org/10.5194/hess-26-4773-2022, https://doi.org/10.5194/hess-26-4773-2022, 2022
Short summary
Short summary
The
Digital Earth(DE) metaphor is very useful for both end users and hydrological modelers. We analyse different categories of models, with the view of making them part of a Digital eARth Twin Hydrology system (called DARTH). We also stress the idea that DARTHs are not models in and of themselves, rather they need to be built on an appropriate information technology infrastructure. It is remarked that DARTHs have to, by construction, support the open-science movement and its ideas.
Niccolò Tubini and Riccardo Rigon
Geosci. Model Dev., 15, 75–104, https://doi.org/10.5194/gmd-15-75-2022, https://doi.org/10.5194/gmd-15-75-2022, 2022
Short summary
Short summary
This paper presents WHETGEO and its 1D deployment: a new physically based model simulating the water and energy budgets in a soil column. WHETGEO-1D is intended to be the first building block of a new customisable land-surface model that is integrated with process-based hydrology. WHETGEO is developed as an open-source code and is fully integrated into the GEOframe/OMS3 system, allowing the use of the many ancillary tools it provides.
Niccolò Tubini and Stephan Gruber
EGUsphere, https://doi.org/10.5194/egusphere-2025-2649, https://doi.org/10.5194/egusphere-2025-2649, 2025
This preprint is open for discussion and under review for The Cryosphere (TC).
Short summary
Short summary
This research introduces a new model for simulating how melting ground ice in permafrost reshapes the land surface over time. It shows that small differences in soil and the depth where ice is found can cause large differences in how the ground sinks or rises. This helps improves our ability to predict future impacts on terrain, ecosystems, and infrastructure as the climate warms.
Nicholas Brown and Stephan Gruber
EGUsphere, https://doi.org/10.5194/egusphere-2025-2658, https://doi.org/10.5194/egusphere-2025-2658, 2025
This preprint is open for discussion and under review for The Cryosphere (TC).
Short summary
Short summary
This study improves how we track changes in permafrost by testing new ways to use ground temperature data. A set of five simple but powerful metrics was found to give a clearer picture of thawing than current methods. The results also show that the depth where sensors are placed can strongly affect measured warming rates. These findings help make permafrost monitoring more accurate and support better planning for a changing climate.
Hosein Fereydooni, Stephan Gruber, David Stillman, and Derek Cronmiller
EGUsphere, https://doi.org/10.5194/egusphere-2025-1801, https://doi.org/10.5194/egusphere-2025-1801, 2025
Short summary
Short summary
Detecting ground ice in permafrost is crucial for climate research and infrastructure, but traditional methods often struggle to distinguish it. This study examines the dielectric properties of ground ice as a unique fingerprint. Field measurements were taken at two Yukon permafrost sites: a retrogressive thaw slump and a pingo. Comparing these with electrical resistivity and impedance results, we found relaxation time is a more reliable indicator for ground ice detection.
Bin Cao and Stephan Gruber
EGUsphere, https://doi.org/10.5194/egusphere-2025-575, https://doi.org/10.5194/egusphere-2025-575, 2025
Short summary
Short summary
The climate-driven changes in cold regions have an outsized importance for local resilient communities and for global climate through teleconnections. We shows that reanalyses are less accurate in cold regions compared to other more populated regions, coincident with the low density of observations. Our findings likely point to similar gaps in our knowledge and capabilities for climate research and services in cold regions.
Concetta D'Amato, Niccolò Tubini, and Riccardo Rigon
EGUsphere, https://doi.org/10.5194/egusphere-2024-4128, https://doi.org/10.5194/egusphere-2024-4128, 2025
Short summary
Short summary
This paper presents GEOSPACE and its 1D implementation: an open-source tool for simulating soil-plant-atmosphere continuum (SPAC) interactions. Using object-oriented programming, GEOSPACE modularizes SPAC processes, focusing on infiltration, evapotranspiration, and root water uptake. The 1D deployment integrates plant transpiration, soil evaporation, and root growth, providing a flexible and validated framework for ecohydrological modeling and applications.
Shima Azimi, Christian Massari, Giuseppe Formetta, Silvia Barbetta, Alberto Tazioli, Davide Fronzi, Sara Modanesi, Angelica Tarpanelli, and Riccardo Rigon
Hydrol. Earth Syst. Sci., 27, 4485–4503, https://doi.org/10.5194/hess-27-4485-2023, https://doi.org/10.5194/hess-27-4485-2023, 2023
Short summary
Short summary
We analyzed the water budget of nested karst catchments using simple methods and modeling. By utilizing the available data on precipitation and discharge, we were able to determine the response lag-time by adopting new techniques. Additionally, we modeled snow cover dynamics and evapotranspiration with the use of Earth observations, providing a concise overview of the water budget for the basin and its subbasins. We have made the data, models, and workflows accessible for further study.
Alessandro Cicoira, Samuel Weber, Andreas Biri, Ben Buchli, Reynald Delaloye, Reto Da Forno, Isabelle Gärtner-Roer, Stephan Gruber, Tonio Gsell, Andreas Hasler, Roman Lim, Philippe Limpach, Raphael Mayoraz, Matthias Meyer, Jeannette Noetzli, Marcia Phillips, Eric Pointner, Hugo Raetzo, Cristian Scapozza, Tazio Strozzi, Lothar Thiele, Andreas Vieli, Daniel Vonder Mühll, Vanessa Wirz, and Jan Beutel
Earth Syst. Sci. Data, 14, 5061–5091, https://doi.org/10.5194/essd-14-5061-2022, https://doi.org/10.5194/essd-14-5061-2022, 2022
Short summary
Short summary
This paper documents a monitoring network of 54 positions, located on different periglacial landforms in the Swiss Alps: rock glaciers, landslides, and steep rock walls. The data serve basic research but also decision-making and mitigation of natural hazards. It is the largest dataset of its kind, comprising over 209 000 daily positions and additional weather data.
Francisco José Cuesta-Valero, Hugo Beltrami, Stephan Gruber, Almudena García-García, and J. Fidel González-Rouco
Geosci. Model Dev., 15, 7913–7932, https://doi.org/10.5194/gmd-15-7913-2022, https://doi.org/10.5194/gmd-15-7913-2022, 2022
Short summary
Short summary
Inversions of subsurface temperature profiles provide past long-term estimates of ground surface temperature histories and ground heat flux histories at timescales of decades to millennia. Theses estimates complement high-frequency proxy temperature reconstructions and are the basis for studying continental heat storage. We develop and release a new bootstrap method to derive meaningful confidence intervals for the average surface temperature and heat flux histories from any number of profiles.
Riccardo Rigon, Giuseppe Formetta, Marialaura Bancheri, Niccolò Tubini, Concetta D'Amato, Olaf David, and Christian Massari
Hydrol. Earth Syst. Sci., 26, 4773–4800, https://doi.org/10.5194/hess-26-4773-2022, https://doi.org/10.5194/hess-26-4773-2022, 2022
Short summary
Short summary
The
Digital Earth(DE) metaphor is very useful for both end users and hydrological modelers. We analyse different categories of models, with the view of making them part of a Digital eARth Twin Hydrology system (called DARTH). We also stress the idea that DARTHs are not models in and of themselves, rather they need to be built on an appropriate information technology infrastructure. It is remarked that DARTHs have to, by construction, support the open-science movement and its ideas.
Élise G. Devoie, Stephan Gruber, and Jeffrey M. McKenzie
Earth Syst. Sci. Data, 14, 3365–3377, https://doi.org/10.5194/essd-14-3365-2022, https://doi.org/10.5194/essd-14-3365-2022, 2022
Short summary
Short summary
Soil freezing characteristic curves (SFCCs) relate the temperature of a soil to its ice content. SFCCs are needed in all physically based numerical models representing freezing and thawing soils, and they affect the movement of water in the subsurface, biogeochemical processes, soil mechanics, and ecology. Over a century of SFCC data exist, showing high variability in SFCCs based on soil texture, water content, and other factors. This repository summarizes all available SFCC data and metadata.
Niccolò Tubini and Riccardo Rigon
Geosci. Model Dev., 15, 75–104, https://doi.org/10.5194/gmd-15-75-2022, https://doi.org/10.5194/gmd-15-75-2022, 2022
Short summary
Short summary
This paper presents WHETGEO and its 1D deployment: a new physically based model simulating the water and energy budgets in a soil column. WHETGEO-1D is intended to be the first building block of a new customisable land-surface model that is integrated with process-based hydrology. WHETGEO is developed as an open-source code and is fully integrated into the GEOframe/OMS3 system, allowing the use of the many ancillary tools it provides.
John Mohd Wani, Renoj J. Thayyen, Chandra Shekhar Prasad Ojha, and Stephan Gruber
The Cryosphere, 15, 2273–2293, https://doi.org/10.5194/tc-15-2273-2021, https://doi.org/10.5194/tc-15-2273-2021, 2021
Short summary
Short summary
We study the surface energy balance from a cold-arid permafrost environment in the Indian Himalayan region. The GEOtop model was used for the modelling of surface energy balance. Our results show that the variability in the turbulent heat fluxes is similar to that reported from the seasonally frozen ground and permafrost regions of the Tibetan Plateau. Further, the low relative humidity could be playing a critical role in the surface energy balance and the permafrost processes.
Rupesh Subedi, Steven V. Kokelj, and Stephan Gruber
The Cryosphere, 14, 4341–4364, https://doi.org/10.5194/tc-14-4341-2020, https://doi.org/10.5194/tc-14-4341-2020, 2020
Short summary
Short summary
Permafrost beneath tundra near Lac de Gras (Northwest Territories, Canada) contains more ice and less organic carbon than shown in global compilations. Excess-ice content of 20–60 %, likely remnant Laurentide basal ice, is found in upland till. This study is based on 24 boreholes up to 10 m deep. Findings highlight geology and glacial legacy as determinants of a mosaic of permafrost characteristics with potential for thaw subsidence up to several metres in some locations.
Cited articles
Anderson, D. M. and Tice, A. R.: Predicting unfrozen water contents in frozen
soils from surface area measurements, Highway research record, 393, 12–18,
1972. a
Andreas, E. L.: Handbook of physical constants and functions for use in
atmospheric boundary layer studies, Cold Regions Research and Engineering Laboratory, US Army Engineer Research and Development Center, 2005. a
Aschwanden, A. and Blatter, H.: Meltwater production due to strain heating in
Storglaciären, Sweden, J. Geophys. Res.-Earth Surf.,
110, F04024, https://doi.org/10.1029/2005JF000328, 2005. a
Bancheri, M.: A flexible approach to the estimation of water budgets and its
connection to the travel time theory, PhD thesis, University of Trento, Trento,
2017. a
Bartelt, P. and Lehning, M.: A physical SNOWPACK model for the Swiss avalanche
warning: Part I: numerical model, Cold Reg. Sci. Technol., 35,
123–145, 2002. a
Bouyoucos, G.: Degree of temperature to which soils can be cooled without
freezing, Mon. Weather Rev., 48, 718–718, 1920. a
Bouyoucos, G. and McCool, M.: The freezing point method as a new means of
measuring the concentration of the soil solution directly in the soil,
Mich, Agr. Exp. Sta. Tech. Bull., 24, 44 pp., 1915. a
Bouyoucos, G. J.: An investigation of soil temperature and some of the most
important facters influencing it, Technical Bulletin of Michigan Agriculture
Experimental Station, 17, 1–196, 1913. a
Bouyoucos, G. J.: Movement of soil moisture from small capillaries to the large
capillaries of the soil upon freezing, J. Agric. Res., 24,
427–432, 1923. a
Brugnano, L. and Casulli, V.: Iterative solution of piecewise linear systems,
SIAM J. Sci. Comput., 30, 463–472, 2008. a
Brugnano, L. and Casulli, V.: Iterative solution of piecewise linear systems
and applications to flows in porous media, SIAM J. Sci.
Comput., 31, 1858–1873, 2009. a
Brun, E., Martin, E., Simon, V., Gendre, C., and Coleou, C.: An energy and mass
model of snow cover suitable for operational avalanche forecasting, J. Glaciol., 35, 333–342, 1989. a
Callen, H. B.: Thermodynamics and an Introduction to Thermostatistics, John
Wiley and Sons, Inc., Republic of Singapore, 1985. a
Casulli, V. and Walters, R. A.: An unstructured grid, three-dimensional model
based on the shallow water equations, Int. J. Numer.
Meth. Fluids, 32, 331–348, 2000. a
Celia, M. A., Bouloutas, E. T., and Zarba, R. L.: A general mass-conservative
numerical solution for the unsaturated flow equation, Water Resour.
Res., 26, 1483–1496, 1990. a
Chistyakov, V.: On mappings of bounded variation, J. Dyn.
Control Syst., 3, 261–289, 1997. a
Clow, G. D.: CVPM 1.1: a flexible heat-transfer modeling system for permafrost, Geosci. Model Dev., 11, 4889–4908, https://doi.org/10.5194/gmd-11-4889-2018, 2018. a, b
Dai, Y., Zeng, X., Dickinson, R. E., Baker, I., Bonan, G. B., Bosilovich,
M. G., Denning, A. S., Dirmeyer, P. A., Houser, P. R., Niu, G., Dai, Y., Zeng, X., Dickinson, R. E., Baker, I., Bonan, G. B., Bosilovich, M. G., Denning, A. S., Dirmeyer, P. A., Houser, P. R., Niu, G., Oleson, K. W., Schlosser, C. A., and Yang, Z.-L.: The
common land model, B. Am. Meteorol. Soc., 84,
1013–1024, 2003. a
D'Amboise, C. J. L., Müller, K., Oxarango, L., Morin, S., and Schuler, T. V.: Implementation of a physically based water percolation routine in the Crocus/SURFEX (V7.3) snowpack model, Geosci. Model Dev., 10, 3547–3566, https://doi.org/10.5194/gmd-10-3547-2017, 2017. a
David, O., Ascough II, J. C., Lloyd, W., Green, T. R., Rojas, K., Leavesley,
G. H., and Ahuja, L. R.: A software engineering perspective on environmental
modeling framework design: The Object Modeling System, Environ.
Model. Softw., 39, 201–213, 2013. a
De Lorenzo, S., Di Renzo, V., Civetta, L., D'antonio, M., and Gasparini, P.:
Thermal model of the Vesuvius magma chamber, Geophys. Res. Lett.,
33, L17302, https://doi.org/10.1029/2006GL026587, 2006. a
Ekici, A., Beer, C., Hagemann, S., Boike, J., Langer, M., and Hauck, C.: Simulating high-latitude permafrost regions by the JSBACH terrestrial ecosystem model, Geosci. Model Dev., 7, 631–647, https://doi.org/10.5194/gmd-7-631-2014, 2014. a, b
Endrizzi, S., Gruber, S., Dall'Amico, M., and Rigon, R.: GEOtop 2.0: simulating the combined energy and water balance at and below the land surface accounting for soil freezing, snow cover and terrain effects, Geosci. Model Dev., 7, 2831–2857, https://doi.org/10.5194/gmd-7-2831-2014, 2014. a
Foley, J. A., Prentice, I. C., Ramankutty, N., Levis, S., Pollard, D., Sitch,
S., and Haxeltine, A.: An integrated biosphere model of land surface
processes, terrestrial carbon balance, and vegetation dynamics, Global
Biogeochem. Cycles, 10, 603–628, 1996. a
Formetta, G., Antonello, A., Franceschi, S., David, O., and Rigon, R.:
Hydrological modelling with components: A GIS-based open-source framework,
Environ. Model. Softw., 55, 190–200, 2014. a
Goodrich, L.: Some results of a numerical study of ground thermal regimes, in:
Proceedings of the Third International Conference on Permafrost,
National Research Council of Canada, Ottawa, Edmonton, Canada,
1978. a
Gubler, S., Endrizzi, S., Gruber, S., and Purves, R. S.: Sensitivities and uncertainties of modeled ground temperatures in mountain environments, Geosci. Model Dev., 6, 1319–1336, https://doi.org/10.5194/gmd-6-1319-2013, 2013. a, b
Harris, C., Arenson, L. U., Christiansen, H. H., Etzelmüller, B., Frauenfelder, R., Gruber, S., Haeberli, W., Hauck, C., Hölzle, M., Humlum, O., Isaksen, K., Kääb, A., Kern-Lütschg, M. A., Lehning, M., Matsuoka, N., Murton, J. B., Noetzli, J., Phillips, M., Ross, N., Seppälä, M., Springman, S. M., and Vonder Mühll, D.: Permafrost and climate in Europe: Monitoring and modelling
thermal, geomorphological and geotechnical responses, Earth-Science Rev.,
92, 117–171, 2009. a
Hollesen, J., Elberling, B., and Jansson, P.-E.: Future active layer dynamics
and carbon dioxide production from thawing permafrost layers in Northeast
Greenland, Global Change Biol., 17, 911–926, 2011. a
Hu, H. and Argyropoulos, S. A.: Mathematical modelling of solidification and
melting: a review, Model. Simul. Mater. Sc., 4, 371–396, 1996. a
InterFrost Project: InterFrost Project,
available at: https://wiki.lsce.ipsl.fr/interfrost/doku.php?id=test_cases:one,
last access: 20 August 2020. a
Jansson, P. and Karlberg, L.: Coupled heat and mass transfer model for
soil-plant-atmosphere systems, Royal Institute of Technology, Dept. of Civil
and Environmental Engineering, Stockholm, 2011. a
Kozlowski, T.: A semi-empirical model for phase composition of water in
clay–water systems, Cold Reg. Sci. Technol., 49, 226–236, 2007. a
Krinner, G., Viovy, N., de Noblet-Ducoudré, N., Ogée, J., Polcher, J.,
Friedlingstein, P., Ciais, P., Sitch, S., and Prentice, I. C.: A dynamic
global vegetation model for studies of the coupled atmosphere-biosphere
system, Global Biogeochem. Cycles, 19, GB1015, https://doi.org/10.1029/2003GB002199, 2005. a
Kurylyk, B. L., MacQuarrie, K. T., and McKenzie, J. M.: Climate change impacts
on groundwater and soil temperatures in cold and temperate regions:
Implications, mathematical theory, and emerging simulation tools,
Earth-Sci. Rev., 138, 313–334, 2014a. a
Lafaysse, M., Cluzet, B., Dumont, M., Lejeune, Y., Vionnet, V., and Morin, S.: A multiphysical ensemble system of numerical snow modelling, The Cryosphere, 11, 1173–1198, https://doi.org/10.5194/tc-11-1173-2017, 2017. a
Lawrence, D., Fisher, R., Koven, C., Oleson, K., Swenson, S., and Vertenstein,
M.: CLM5 documentation, Tech. rep., Tech. rep., National Center
for Atmospheric Research, Boulder, CO, 2019. a
Lewis, R. and Ravindran, K.: Finite element simulation of metal casting,
Int. J. Numer. Meth. Eng., 47, 29–59, 2000. a
Lliboutry, L. and Duval, P.: Various isotropic and anisotropic ices found in
glaciers and polar ice caps and their corresponding rheologies, Ann.
Geophys., 3, 207–224, 1985. a
Marchenko, S., Romanovsky, V., and Tipenko, G.: Numerical modeling of spatial
permafrost dynamics in Alaska, in: Proceedings of the ninth international
conference on permafrost, 29, 1125–1130, Institute of Northern
Engineering, University of Alaska, Fairbanks, 2008. a
McKenzie, J. M. and Voss, C. I.: Permafrost thaw in a nested groundwater-flow
system, Hydrogeol. J., 21, 299–316, 2013. a
Mongibello, L., Bianco, N., Caliano, M., and Graditi, G.: Numerical simulation of an aluminum container including a phase change material for cooling energy storage, Applied System Innovation, 1, 34, 2018. a
Morin, S., Lejeune, Y., Lesaffre, B., Panel, J.-M., Poncet, D., David, P., and Sudul, M.: An 18-yr long (1993–2011) snow and meteorological dataset from a mid-altitude mountain site (Col de Porte, France, 1325 m alt.) for driving and evaluating snowpack models, Earth Syst. Sci. Data, 4, 13–21, https://doi.org/10.5194/essd-4-13-2012, 2012. a
Nazzi Ehms, J. H., De Césaro Oliveski, R., Oliveira Rocha, L. A., Biserni, C., and Garai, M.: Fixed grid numerical models for solidification and melting of phase change materials (PCMs), Appl. Sci., 9, 4334, 2019. a
Nedjar, B.: An enthalpy-based finite element method for nonlinear heat problems
involving phase change, Comput. Struct., 80, 9–21, 2002. a
Nicolsky, D., Romanovsky, V., Alexeev, V., and Lawrence, D.: Improved modeling
of permafrost dynamics in a GCM land-surface scheme, Geophys. Res.
Lett., 34, L08501, https://doi.org/10.1029/2007GL029525, 2007a. a, b
Nicolsky, D. J., Romanovsky, V. E., and Tipenko, G. S.: Using in-situ temperature measurements to estimate saturated soil thermal properties by solving a sequence of optimization problems, The Cryosphere, 1, 41–58, https://doi.org/10.5194/tc-1-41-2007, 2007.b. a, b, c
Nye, J. F. and Frank, F. C.: Hydrology of the intergranular veins in a temperate glacier.InInternational Association of Scientific Hydrology Publication, Symp. Cambridge, 1969, Hydrol. Glaciers, 95, 157–161,
1973. a
Rigon, R., Bertoldi, G., and Over, T. M.: GEOtop: A distributed hydrological
model with coupled water and energy budgets, J. Hydrometeorol., 7,
371–388, 2006. a
Ruhaak, W., Anbergen, H., Grenier, C., McKenzie, J., Kurylyk, B., Molson, J.,
Roux, N., and Sass, I.: Benchmarking numerical freeze/thaw models, Energy
Procedia, 76, 301–310, 2015. a
Schuur, E. A. G., David McGuire, A., Schädel, C., Grosse, G., Harden, J. W., Hayes, D. J., Hugelius, G., Koven, C. D., Kuhry, P., Lawrence, D. M., Natali, S. M., Olefeldt, D., Romanovsky, V. E., Schaefer, K., Turetsky, M. R., Treat, C. C., and Vonk, J. E.:
Climate change and the permafrost carbon feedback, Nature, 520, 171–179,
2015. a
Sheshukov, A. Y. and Nieber, J. L.: One-dimensional freezing of nonheaving
unsaturated soils: Model formulation and similarity solution, Water Resour.
Res., 47, W11519, https://doi.org/10.1029/2011WR010512, 2011. a
Shewchuk, J. R.: An Introduction to the Conjugate Gradient Method Without the
Agonizing Pain, Tech. rep., USA, 1994. a
Streletskiy, D. A., Suter, L. J., Shiklomanov, N. I., Porfiriev, B. N., and
Eliseev, D. O.: Assessment of climate change impacts on buildings, structures
and infrastructure in the Russian regions on permafrost, Environ.
Res. Lett., 14, 025003, https://doi.org/10.1088/1748-9326/aaf5e6, 2019. a
Tan, X., Chen, W., Tian, H., and Cao, J.: Water flow and heat transport
including ice/water phase change in porous media: numerical simulation and
application, Cold Reg. Sci. Technol., 68, 74–84, 2011. a
Teng, J., Kou, J., Yan, X., Zhang, S., and Sheng, D.: Parameterization of soil freezing characteristic curve for unsaturated soils, Cold Reg. Sci. Technol., 170, 102928, https://doi.org/10.1016/j.neubiorev.2019.07.019, 2020. a
Tubini, N.: OMS project,
https://github.com/GEOframeOMSProjects/OMS_FreeThaw1D [Code], last
access: 16 October 2020a. a
Tubini, N.: source code,
https://github.com/geoframecomponents/FreeThaw1D [Code], last
access: 16 October 2020b. a
Tubini, N., Gruber, S., and Rigon, R.: FreeThaw1D (Version v-0.98) [Code], Zenodo, https://doi.org/10.5281/zenodo.4017668, 2020. a
Verseghy, D. L.: CLASS – A Canadian land surface scheme for GCMs. I. Soil
model, Int. J. Climatol., 11, 111–133, 1991. a
Vionnet, V., Brun, E., Morin, S., Boone, A., Faroux, S., Le Moigne, P., Martin, E., and Willemet, J.-M.: The detailed snowpack scheme Crocus and its implementation in SURFEX v7.2, Geosci. Model Dev., 5, 773–791, https://doi.org/10.5194/gmd-5-773-2012, 2012. a
Voller, V. and Cross, M.: Accurate solutions of moving boundary problems using
the enthalpy method, Int. J. Heat Mass Transf., 24,
545–556, 1981. a
Voller, V. R.: Fast implicit finite-difference method for the analysis of phase
change problems, Numer. Heat Transf., 17, 155–169, 1990. a
Voss, C. I. and Provost, A.: SUTRA: A model for 2D or 3D saturated-unsaturated,
variable-density ground-water flow with solute or energy transport, Tech.
rep., U.S. Geological Survey Water-Resources Investigations Report, Reston, Virginia, 2002. a
Vuik, C.: Some historical notes about the Stefan problem, Nieuw Arch. Wiskd., 11, 157–167, 1993. a
Watanabe, K., Kito, T., Wake, T., and Sakai, M.: Freezing experiments on
unsaturated sand, loam and silt loam, Ann. Glaciol., 52, 37–43, 2011. a
Westermann, S., Schuler, T. V., Gisnås, K., and Etzelmüller, B.: Transient thermal modeling of permafrost conditions in Southern Norway, The Cryosphere, 7, 719–739, https://doi.org/10.5194/tc-7-719-2013, 2013. a, b
Westermann, S., Langer, M., Boike, J., Heikenfeld, M., Peter, M., Etzelmüller, B., and Krinner, G.: Simulating the thermal regime and thaw processes of ice-rich permafrost ground with the land-surface model CryoGrid 3, Geosci. Model Dev., 9, 523–546, https://doi.org/10.5194/gmd-9-523-2016, 2016. a
Wever, N., Fierz, C., Mitterer, C., Hirashima, H., and Lehning, M.: Solving Richards Equation for snow improves snowpack meltwater runoff estimations in detailed multi-layer snowpack model, The Cryosphere, 8, 257–274, https://doi.org/10.5194/tc-8-257-2014, 2014. a
Zhang, Y., Chen, W., and Cihlar, J.: A process-based model for quantifying the
impact of climate change on permafrost thermal regimes, J. Geophys. Res.-Atmos., 108, 4695, 2003. a
Zhang, Y., Carey, S. K., and Quinton, W. L.: Evaluation of the algorithms and
parameterizations for ground thawing and freezing simulation in permafrost
regions, J. Geophys. Res.-Atmos., 113, D17116, https://doi.org/10.1029/2007JD009343, 2008. a, b, c
Short summary
We present a new method to compute temperature changes with melting and freezing – a fundamental challenge in cryosphere research – extremely efficiently and with guaranteed correctness of the energy balance for any time step size. This is a key feature since the integration time step can then be chosen according to the timescale of the processes to be studied, from seconds to days.
We present a new method to compute temperature changes with melting and freezing – a fundamental...
