Articles | Volume 8, issue 5
https://doi.org/10.5194/tc-8-1935-2014
© Author(s) 2014. This work is distributed under
the Creative Commons Attribution 3.0 License.
the Creative Commons Attribution 3.0 License.
https://doi.org/10.5194/tc-8-1935-2014
© Author(s) 2014. This work is distributed under
the Creative Commons Attribution 3.0 License.
the Creative Commons Attribution 3.0 License.
Research article
| 
23 Oct 2014
Research article |
| 23 Oct 2014
Three-phase numerical model for subsurface hydrology in permafrost-affected regions (PFLOTRAN-ICE v1.0)
S. Karra
Computational Earth Science Group, Earth and Environmental Sciences Division, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
S. L. Painter
Computational Earth Science Group, Earth and Environmental Sciences Division, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
current address: Climate Change Science Institute, Environmental Sciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA
P. C. Lichtner
OFM Research, 28430 NE 47th Pl, Redmond, WA 98053, USA
Related authors
Scott K. Hansen, Sachin Pandey, Satish Karra, and Velimir V. Vesselinov
Geosci. Model Dev., 10, 4525–4538, https://doi.org/10.5194/gmd-10-4525-2017, https://doi.org/10.5194/gmd-10-4525-2017, 2017
Short summary
Short summary
Groundwater contamination by heavy metals is a critical environmental problem for which in situ remediation is frequently the only viable treatment option. We developed a computer model for simulating such interventions, which includes custom reaction kinetics equations for the coupled interactions between the metal, the existing microbial population, and a number of possible chemical amendments. We describe installation and usage and present two instructive example simulations.
Michael Nole, Katherine Muller, Glenn Hammond, Xiaoliang He, and Peter Lichtner
EGUsphere, https://doi.org/10.5194/egusphere-2025-1343, https://doi.org/10.5194/egusphere-2025-1343, 2025
This preprint is open for discussion and under review for Geoscientific Model Development (GMD).
Short summary
Short summary
Subsurface injection of carbon dioxide (CO2) can be used for a variety of purposes including geologic carbon storage and enhanced oil recovery. Recently, CO2 injection into reactive host rocks has been explored as a way to transform CO2 into dense solid minerals. We present a simulation framework for modeling flow of CO2 due to injection and subsequent reactions that take place to mineralize CO2.
Sergi Molins, Benjamin J. Andre, Jeffrey N. Johnson, Glenn E. Hammond, Benjamin N. Sulman, Konstantin Lipnikov, Marcus S. Day, James J. Beisman, Daniil Svyatsky, Hang Deng, Peter C. Lichtner, Carl I. Steefel, and J. David Moulton
Geosci. Model Dev., 18, 3241–3263, https://doi.org/10.5194/gmd-18-3241-2025, https://doi.org/10.5194/gmd-18-3241-2025, 2025
Short summary
Short summary
Developing scientific software and making sure it functions properly requires a significant effort. As we advance our understanding of natural systems, however, there is the need to develop yet more complex models and codes. In this work, we present a piece of software that facilitates this work, specifically with regard to reactive processes. Existing tried-and-true codes are made available via this new interface, freeing up resources to focus on the new aspects of the problems at hand.
Scott K. Hansen, Sachin Pandey, Satish Karra, and Velimir V. Vesselinov
Geosci. Model Dev., 10, 4525–4538, https://doi.org/10.5194/gmd-10-4525-2017, https://doi.org/10.5194/gmd-10-4525-2017, 2017
Short summary
Short summary
Groundwater contamination by heavy metals is a critical environmental problem for which in situ remediation is frequently the only viable treatment option. We developed a computer model for simulating such interventions, which includes custom reaction kinetics equations for the coupled interactions between the metal, the existing microbial population, and a number of possible chemical amendments. We describe installation and usage and present two instructive example simulations.
Guoping Tang, Fengming Yuan, Gautam Bisht, Glenn E. Hammond, Peter C. Lichtner, Jitendra Kumar, Richard T. Mills, Xiaofeng Xu, Ben Andre, Forrest M. Hoffman, Scott L. Painter, and Peter E. Thornton
Geosci. Model Dev., 9, 927–946, https://doi.org/10.5194/gmd-9-927-2016, https://doi.org/10.5194/gmd-9-927-2016, 2016
Short summary
Short summary
We demonstrate that CLM-PFLOTRAN predictions are consistent with CLM4.5 for Arctic, temperate, and tropical sites. A tight relative tolerance may be needed to avoid false convergence when scaling, clipping, or log transformation is used to avoid negative concentration in implicit time stepping and Newton-Raphson methods. The log transformation method is accurate and robust while relaxing relative tolerance or using the clipping or scaling method can result in efficient solutions.
Related subject area
Numerical Modelling
Calibrated sea level contribution from the Amundsen Sea sector, West Antarctica, under RCP8.5 and Paris 2C scenarios
Calving driven by horizontal forces in a revised crevasse-depth framework
The demise of the world's largest piedmont glacier: a probabilistic forecast
Improved basal drag of the West Antarctic Ice Sheet from L-curve analysis of inverse models utilizing subglacial hydrology simulations
Disentangling the oceanic drivers behind the post-2000 retreat of Sermeq Kujalleq, Greenland (Jakobshavn Isbræ)
Glacier damage evolution over ice flow timescales
Inter-model differences in 21st century glacier runoff for the world's major river basins
Quantifying radiative effects of light-absorbing particle deposition on snow at the SnowMIP sites
A minimal machine-learning glacier mass balance model
Improving large-scale snow albedo modeling using a climatology of light-absorbing particle deposition
Physically based modelling of glacier evolution under climate change in the tropical Andes
The thermal state of permafrost in under climate change on the Qinghai-Tibet Plateau from 1980 to 2022: A case study of the West Kunlun
Present-day mass loss rates are a precursor for West Antarctic Ice Sheet collapse
A hybrid ice-mélange model based on particle and continuum methods
New glacier thickness and bed topography maps for Svalbard
Sea-level rise contribution from Ryder Glacier in Northern Greenland varies by an order of magnitude by 2300 depending on future emissions
Quantifying the buttressing contribution of landfast sea ice and melange to Crane Glacier, Antarctic Peninsula
Multi-physics ensemble modelling of Arctic tundra snowpack properties
Combining observational data and numerical models to obtain a seamless high temporal resolution seasonal cycle of snow and ice mass balance at the MOSAiC Central Observatory
Two-way coupling between ice flow and channelized subglacial drainage enhances modeled marine-ice-sheet retreat
Sensitivity of the future evolution of the Wilkes Subglacial Basin ice sheet to grounding-line melt parameterizations
Modelling snowpack on ice surfaces with the ORCHIDEE land surface model: application to the Greenland ice sheet
Application of a regularised Coulomb sliding law to Jakobshavn Isbræ, western Greenland
Brief communication: Stalagmite damage by cave ice flow quantitatively assessed by fluid–structure interaction simulations
Exploring the decision-making process in model development: focus on the Arctic snowpack
Exploring the potential of forest snow modeling at the tree and snowpack layer scale
Simulating lake ice phenology using a coupled atmosphere–lake model at Nam Co, a typical deep alpine lake on the Tibetan Plateau
Modelling the effect of free convection on permafrost melting rates in frozen rock clefts
How many parameters are needed to represent polar sea ice surface patterns and heterogeneity?
Antarctic sensitivity to oceanic melting parameterizations
Analytical solutions for the advective–diffusive ice column in the presence of strain heating
Ice viscosity governs hydraulic fracture that causes rapid drainage of supraglacial lakes
Investigating the impact of reanalysis snow input on an observationally calibrated snow-on-sea-ice reconstruction
Microstructure-based modelling of snow mechanics: experimental evaluation of the cone penetration test
Snow redistribution in an intermediate-complexity snow hydrology modelling framework
Increasing numerical stability of mountain valley glacier simulations: implementation and testing of free-surface stabilization in Elmer/Ice
Analyzing the sensitivity of a blowing snow model (SnowPappus) to precipitation forcing, blowing snow, and spatial resolution
Exploring non-Gaussian sea ice characteristics via observing system simulation experiments
Past and future of the Arctic sea ice in High-Resolution Model Intercomparison Project (HighResMIP) climate models
Age, thinning and spatial origin of the Beyond EPICA ice from a 2.5D ice flow model
Biases in ice sheet models from missing noise-induced drift
A Novel Transformation of the Ice Sheet Stokes Equations and Some of its Properties and Applications
A 3D glacier dynamics–line plume model to estimate the frontal ablation of Hansbreen, Svalbard
Data-driven surrogate modeling of high-resolution sea-ice thickness in the Arctic
Using Icepack to reproduce ice mass balance buoy observations in landfast ice: improvements from the mushy-layer thermodynamics
Modeling the timing of Patagonian Ice Sheet retreat in the Chilean Lake District from 22–10 ka
Decadal re-forecasts of glacier climatic mass balance
Understanding the influence of ocean waves on Arctic sea ice simulation: a modeling study with an atmosphere–ocean–wave–sea ice coupled model
Sea ice cover in the Copernicus Arctic Regional Reanalysis
Regime shifts in Arctic terrestrial hydrology manifested from impacts of climate warming
Sebastian H. R. Rosier, G. Hilmar Gudmundsson, Adrian Jenkins, and Kaitlin A. Naughten
The Cryosphere, 19, 2527–2557, https://doi.org/10.5194/tc-19-2527-2025, https://doi.org/10.5194/tc-19-2527-2025, 2025
Short summary
Short summary
Glaciers in the Amundsen Sea region of Antarctica have been retreating and losing mass, but their future contribution to global sea level rise remains highly uncertain. We use an ice sheet model and uncertainty quantification methods to evaluate the probable range of mass loss from this region for two future climate scenarios. We find that the rate of ice loss until 2100 will likely remain similar to present-day observations, with little sensitivity to climate scenario over this short time frame.
Donald A. Slater and Till J. W. Wagner
The Cryosphere, 19, 2475–2493, https://doi.org/10.5194/tc-19-2475-2025, https://doi.org/10.5194/tc-19-2475-2025, 2025
Short summary
Short summary
Calving is when icebergs break off glaciers and fall into the ocean. It is an important process determining how ice sheets will respond to changes in climate, but it is currently poorly understood and hard to include in numerical models that are used for sea-level projections. We adapted and extended an existing theory for how this process works, better explaining observations showing that calving style depends on how thick the ice is.
Douglas J. Brinkerhoff, Brandon S. Tober, Michael Daniel, Victor Devaux-Chupin, Michael S. Christoffersen, John W. Holt, Christopher F. Larsen, Mark Fahnestock, Michael G. Loso, Kristin M. F. Timm, Russell C. Mitchell, and Martin Truffer
The Cryosphere, 19, 2321–2353, https://doi.org/10.5194/tc-19-2321-2025, https://doi.org/10.5194/tc-19-2321-2025, 2025
Short summary
Short summary
Sít' Tlein is one of the largest glaciers in the world outside of the polar regions, and we know that it has been rapidly thinning. To forecast how this glacier will change in the future, we combine a computer model of ice flow with measurements from many different sources. Our model tells us that with high probability, Sít' Tlein's lower reaches are going to disappear in the next century and a half, creating a new bay or lake along Alaska's coastline.
Lea-Sophie Höyns, Thomas Kleiner, Andreas Rademacher, Martin Rückamp, Michael Wolovick, and Angelika Humbert
The Cryosphere, 19, 2133–2158, https://doi.org/10.5194/tc-19-2133-2025, https://doi.org/10.5194/tc-19-2133-2025, 2025
Short summary
Short summary
The sliding of glaciers over bedrock is influenced by water pressure in the underlying hydrological system and the roughness of the land underneath the glacier. We estimate this roughness through a modeling approach that optimizes this unknown parameter. Additionally, we simulate water pressure, enhancing the reliability of the computed drag at the ice sheet base. The resulting data are provided to other modelers and scientists conducting geophysical field observations.
Ziad Rashed, Alexander A. Robel, and Hélène Seroussi
The Cryosphere, 19, 1775–1788, https://doi.org/10.5194/tc-19-1775-2025, https://doi.org/10.5194/tc-19-1775-2025, 2025
Short summary
Short summary
Sermeq Kujalleq, Greenland's largest glacier, has significantly retreated since the late 1990s in response to warming ocean temperatures. Using a large-ensemble approach, our simulations show that the retreat is mainly initiated by the arrival of warm water but sustained and accelerated by the glacier's position over deeper bed troughs and vigorous calving. We highlight the need for models of ice mélange to project glacier behavior under rapid calving regimes.
Meghana Ranganathan, Alexander A. Robel, Alexander Huth, and Ravindra Duddu
The Cryosphere, 19, 1599–1619, https://doi.org/10.5194/tc-19-1599-2025, https://doi.org/10.5194/tc-19-1599-2025, 2025
Short summary
Short summary
The rate of ice loss from ice sheets is controlled by the flow of ice from the center of the ice sheet and by the internal fracturing of the ice. These processes are coupled; fractures reduce the viscosity of ice and enable more rapid flow, and rapid flow causes the fracturing of ice. We present a simplified way of representing damage that is applicable to long-timescale flow estimates. Using this model, we find that including fracturing in an ice sheet simulation can increase the loss of ice by 13–29 %.
Finn Wimberly, Lizz Ultee, Lilian Schuster, Matthias Huss, David R. Rounce, Fabien Maussion, Sloan Coats, Jonathan Mackay, and Erik Holmgren
The Cryosphere, 19, 1491–1511, https://doi.org/10.5194/tc-19-1491-2025, https://doi.org/10.5194/tc-19-1491-2025, 2025
Short summary
Short summary
Glacier models have historically been used to understand glacier melt’s contribution to sea level rise. The capacity to project seasonal glacier runoff is a relatively recent development for these models. In this study we provide the first model intercomparison of runoff projections for the glacier evolution models capable of simulating future runoff globally. We compare model projections from 2000 to 2100 for all major river basins larger than 3000 km2 with over 30 km2 of initial glacier cover.
Enrico Zorzetto, Paul Ginoux, Sergey Malyshev, and Elena Shevliakova
The Cryosphere, 19, 1313–1334, https://doi.org/10.5194/tc-19-1313-2025, https://doi.org/10.5194/tc-19-1313-2025, 2025
Short summary
Short summary
Light-absorbing particle (LAP) deposition on snow leads to a darkening of the snow surface and can thus accelerate snow melt. Understanding the extent to which different types of LAPs contribute to snow melt is important to both predict changes in water availability and improve global climate model predictions. Here, we extend a recently developed snow model to account for the deposition of LAPs in the snowpack and evaluate the effect of snow darkening on accelerating snow melt.
Marijn van der Meer, Harry Zekollari, Matthias Huss, Jordi Bolibar, Kamilla Hauknes Sjursen, and Daniel Farinotti
The Cryosphere, 19, 805–826, https://doi.org/10.5194/tc-19-805-2025, https://doi.org/10.5194/tc-19-805-2025, 2025
Short summary
Short summary
Glacier retreat poses big challenges, making understanding how climate affects glaciers vital. But glacier measurements worldwide are limited. We created a simple machine-learning model called miniML-MB, which estimates annual changes in glacier mass in the Swiss Alps. As input, miniML-MB uses two climate variables: average temperature (May–Aug) and total precipitation (Oct–Feb). Our model can accurately predict glacier mass from 1961 to 2021 but struggles for extreme years (2022 and 2023).
Manon Gaillard, Vincent Vionnet, Matthieu Lafaysse, Marie Dumont, and Paul Ginoux
The Cryosphere, 19, 769–792, https://doi.org/10.5194/tc-19-769-2025, https://doi.org/10.5194/tc-19-769-2025, 2025
Short summary
Short summary
This study presents an efficient method to improve large-scale snow albedo simulations by considering the spatial variability in light-absorbing particles (LAPs) like black carbon and dust. A global climatology of LAP deposition was created and used to optimize a parameter in the Crocus snow model. Testing at 10 global sites improved albedo predictions by 10 % on average and over 25 % in the Arctic. This method can enhance other snow models' predictions without complex simulations.
Jonathan D. Mackay, Nicholas E. Barrand, David M. Hannah, Emily Potter, Nilton Montoya, and Wouter Buytaert
The Cryosphere, 19, 685–712, https://doi.org/10.5194/tc-19-685-2025, https://doi.org/10.5194/tc-19-685-2025, 2025
Short summary
Short summary
We combine two globally capable glacier evolution models to include processes that are typically neglected but thought to control tropical glacier retreat (e.g. sublimation). We apply the model to Peru's Vilcanota-Urubamba Basin. The model captures observed glacier mass changes,but struggles with surface albedo dynamics. Projections show glacier mass shrinking to 17 % or 6 % of 2000 levels by 2100 under moderate- and high-emission scenarios, respectively.
Jianting Zhao, Lin Zhao, Ze Sun, Guojie Hu, Defu Zou, Minxuan Xiao, Guangyue Liu, Qiangqiang Pang, Erji Du, Zhibin Li, Xiaodong Wu, Yao Xiao, Lingxiao Wang, and Wenxin Zhang
EGUsphere, https://doi.org/10.5194/egusphere-2024-3956, https://doi.org/10.5194/egusphere-2024-3956, 2025
Short summary
Short summary
The thermal regime is a key indicator of permafrost evolution. We quantitatively analyzed the spatiotemporal dynamics of the permafrost status in western Tibet since the 1980s, based on numerical simulations using the enhanced, model-forcing-driven Moving-Grid Permafrost Model. Our simulated results indicated that slow and lagged response of permafrost to climate warming, which closely linked to historical thermal conditions.
Tim van den Akker, William H. Lipscomb, Gunter R. Leguy, Jorjo Bernales, Constantijn J. Berends, Willem Jan van de Berg, and Roderik S. W. van de Wal
The Cryosphere, 19, 283–301, https://doi.org/10.5194/tc-19-283-2025, https://doi.org/10.5194/tc-19-283-2025, 2025
Short summary
Short summary
In this study, we present an improved way of representing ice thickness change rates in an ice sheet model. We apply this method using two ice sheet models of the Antarctic Ice Sheet. We found that the two largest outlet glaciers on the Antarctic Ice Sheet, Thwaites Glacier and Pine Island Glacier, will collapse without further warming on a timescale of centuries. This would cause a sea level rise of about 1.2 m globally.
Saskia Kahl, Carolin Mehlmann, and Dirk Notz
The Cryosphere, 19, 129–141, https://doi.org/10.5194/tc-19-129-2025, https://doi.org/10.5194/tc-19-129-2025, 2025
Short summary
Short summary
Ice mélange, a mixture of sea ice and icebergs, can impact sea-ice–ocean interactions. But climate models do not yet represent it due to computational limits. To address this shortcoming and include ice mélange into climate models, we suggest representing icebergs as particles. We integrate their feedback into mathematical equations used to model the sea-ice motion in climate models. The setup is computationally efficient due to the iceberg particle usage and enables a realistic representation.
Ward van Pelt and Thomas Frank
The Cryosphere, 19, 1–17, https://doi.org/10.5194/tc-19-1-2025, https://doi.org/10.5194/tc-19-1-2025, 2025
Short summary
Short summary
Accurate information on the ice thickness of Svalbard's glaciers is important for assessing the contribution to sea level rise in a present and a future climate. However, direct observations of the glacier bed are scarce. Here, we use an inverse approach and high-resolution surface observations to infer basal conditions. We present and analyse the new bed and thickness maps, quantify the ice volume (6800 km3), and compare these against radar data and previous studies.
Felicity Alice Holmes, Jamie Barnett, Henning Åkesson, Mathieu Morlighem, Johan Nilsson, Nina Kirchner, and Martin Jakobsson
EGUsphere, https://doi.org/10.5194/egusphere-2024-3839, https://doi.org/10.5194/egusphere-2024-3839, 2024
Short summary
Short summary
Northern Greenland contains some of the ice sheet's last remaining glaciers with floating ice tongues. One of these is Ryder Glacier, which has been relatively stable in recent decades in contrast to nearby glaciers. Here, we use a computer model to simulate Ryder Glacier until 2300 under both a low and high emissions scenario. Very high levels of surface melt under a high emissions future leads to a sea-level rise contribution an order of magnitude higher than under a low emissions future.
Richard Parsons, Sainan Sun, G. Hilmar Gudmundsson, Jan Wuite, and Thomas Nagler
The Cryosphere, 18, 5789–5801, https://doi.org/10.5194/tc-18-5789-2024, https://doi.org/10.5194/tc-18-5789-2024, 2024
Short summary
Short summary
In 2022, multi-year landfast sea ice in Antarctica's Larsen B embayment disintegrated, after which time an increase in the rate at which Crane Glacier discharged ice into the ocean was observed. As the fast ice was joined to the glacier terminus, it could provide resistance against the glacier's flow, slowing down the rate of ice discharge. We used numerical modelling to quantify this resistive stress and found that the fast ice provided significant support to Crane prior to its disintegration.
Georgina J. Woolley, Nick Rutter, Leanne Wake, Vincent Vionnet, Chris Derksen, Richard Essery, Philip Marsh, Rosamond Tutton, Branden Walker, Matthieu Lafaysse, and David Pritchard
The Cryosphere, 18, 5685–5711, https://doi.org/10.5194/tc-18-5685-2024, https://doi.org/10.5194/tc-18-5685-2024, 2024
Short summary
Short summary
Parameterisations of Arctic snow processes were implemented into the multi-physics ensemble version of the snow model Crocus (embedded within the Soil, Vegetation, and Snow version 2 land surface model) and evaluated at an Arctic tundra site. Optimal combinations of parameterisations that improved the simulation of density and specific surface area featured modifications that raise wind speeds to increase compaction in surface layers, prevent snowdrift, and increase viscosity in basal layers.
Polona Itkin and Glen E. Liston
EGUsphere, https://doi.org/10.5194/egusphere-2024-3402, https://doi.org/10.5194/egusphere-2024-3402, 2024
Short summary
Short summary
The MOSAiC project provided a year of observations of Arctic snow and sea ice, though some data were interrupted, especially during summer melt onset. We developed a data-model fusion system to produce continuous, high-resolution time series of snow and sea ice parameters. On all three analyzed three ice types snow redistribution correlated with sea ice deformation and level ice thickness was governed by the thinnest fraction of snow cover.
George Lu and Jonathan Kingslake
The Cryosphere, 18, 5301–5321, https://doi.org/10.5194/tc-18-5301-2024, https://doi.org/10.5194/tc-18-5301-2024, 2024
Short summary
Short summary
Water below ice sheets affects ice-sheet motion, while the evolution of ice sheets likewise affects the water below. We create a model that allows for water and ice to affect each other and use it to see how this coupling or lack thereof may impact ice-sheet retreat. We find that coupling an evolving water system with the ice sheet results in more retreat than if we assume unchanging conditions under the ice, which indicates a need to better represent the effects of water in ice-sheet models.
Yu Wang, Chen Zhao, Rupert Gladstone, Thomas Zwinger, Benjamin K. Galton-Fenzi, and Poul Christoffersen
The Cryosphere, 18, 5117–5137, https://doi.org/10.5194/tc-18-5117-2024, https://doi.org/10.5194/tc-18-5117-2024, 2024
Short summary
Short summary
Our research delves into the future evolution of Antarctica's Wilkes Subglacial Basin (WSB) and its potential contribution to sea level rise, focusing on how basal melt is implemented at the grounding line in ice flow models. Our findings suggest that these implementation methods can significantly impact the magnitude of future ice loss projections. Under a high-emission scenario, the WSB ice sheet could undergo massive and rapid retreat between 2200 and 2300.
Sylvie Charbit, Christophe Dumas, Fabienne Maignan, Catherine Ottlé, Nina Raoult, Xavier Fettweis, and Philippe Conesa
The Cryosphere, 18, 5067–5099, https://doi.org/10.5194/tc-18-5067-2024, https://doi.org/10.5194/tc-18-5067-2024, 2024
Short summary
Short summary
The evolution of the Greenland ice sheet is highly dependent on surface melting and therefore on the processes operating at the snow–atmosphere interface and within the snow cover. Here we present new developments to apply a snow model to the Greenland ice sheet. The performance of this model is analysed in terms of its ability to simulate ablation processes. Our analysis shows that the model performs well when compared with the MAR regional polar atmospheric model.
Matt Trevers, Antony J. Payne, and Stephen L. Cornford
The Cryosphere, 18, 5101–5115, https://doi.org/10.5194/tc-18-5101-2024, https://doi.org/10.5194/tc-18-5101-2024, 2024
Short summary
Short summary
The form of the friction law which determines the speed of ice sliding over the bedrock remains a major source of uncertainty in ice sheet model projections of future sea level rise. Jakobshavn Isbræ, the fastest-flowing glacier in Greenland, which has undergone significant changes in the last few decades, is an ideal case for testing sliding laws. We find that a regularised Coulomb friction law reproduces the large seasonal and inter-annual flow speed variations most accurately.
Alexander H. Jarosch, Paul Hofer, and Christoph Spötl
The Cryosphere, 18, 4811–4816, https://doi.org/10.5194/tc-18-4811-2024, https://doi.org/10.5194/tc-18-4811-2024, 2024
Short summary
Short summary
Mechanical damage to stalagmites is commonly observed in mid-latitude caves. In this study we investigate ice flow along the cave bed as a possible mechanism for stalagmite damage. Utilizing models which simulate forces created by ice flow, we study the structural integrity of different stalagmite geometries. Our results suggest that structural failure of stalagmites caused by ice flow is possible, albeit unlikely.
Cecile B. Menard, Sirpa Rasmus, Ioanna Merkouriadi, Gianpaolo Balsamo, Annett Bartsch, Chris Derksen, Florent Domine, Marie Dumont, Dorothee Ehrich, Richard Essery, Bruce C. Forbes, Gerhard Krinner, David Lawrence, Glen Liston, Heidrun Matthes, Nick Rutter, Melody Sandells, Martin Schneebeli, and Sari Stark
The Cryosphere, 18, 4671–4686, https://doi.org/10.5194/tc-18-4671-2024, https://doi.org/10.5194/tc-18-4671-2024, 2024
Short summary
Short summary
Computer models, like those used in climate change studies, are written by modellers who have to decide how best to construct the models in order to satisfy the purpose they serve. Using snow modelling as an example, we examine the process behind the decisions to understand what motivates or limits modellers in their decision-making. We find that the context in which research is undertaken is often more crucial than scientific limitations. We argue for more transparency in our research practice.
Giulia Mazzotti, Jari-Pekka Nousu, Vincent Vionnet, Tobias Jonas, Rafife Nheili, and Matthieu Lafaysse
The Cryosphere, 18, 4607–4632, https://doi.org/10.5194/tc-18-4607-2024, https://doi.org/10.5194/tc-18-4607-2024, 2024
Short summary
Short summary
As many boreal and alpine forests have seasonal snow, models are needed to predict forest snow under future environmental conditions. We have created a new forest snow model by combining existing, very detailed model components for the canopy and the snowpack. We applied it to forests in Switzerland and Finland and showed how complex forest cover leads to a snowpack layering that is very variable in space and time because different processes prevail at different locations in the forest.
Xu Zhou, Binbin Wang, Xiaogang Ma, Zhu La, and Kun Yang
The Cryosphere, 18, 4589–4605, https://doi.org/10.5194/tc-18-4589-2024, https://doi.org/10.5194/tc-18-4589-2024, 2024
Short summary
Short summary
The simulation of the ice phenology of Nam Co by WRF is investigated. Compared with the default model, improving the key lake schemes, such as water surface roughness length for heat fluxes and the shortwave radiation transfer for lake ice, can better simulate the lake ice phenology. The still existing errors in the spatial patterns of lake ice phenology imply that challenges still exist in modelling key lake and non-lake physics such as grid-scale water circulation and snow-related processes.
Amir Sedaghatkish, Frédéric Doumenc, Pierre-Yves Jeannin, and Marc Luetscher
The Cryosphere, 18, 4531–4546, https://doi.org/10.5194/tc-18-4531-2024, https://doi.org/10.5194/tc-18-4531-2024, 2024
Short summary
Short summary
We developed a model to simulate the natural convection of water within frozen rock crevices subject to daily warming in mountain permafrost regions. Traditional models relying on conduction and latent heat flux typically overlook free convection. The results reveal that free convection can significantly accelerate the melting rate by an order of magnitude compared to conduction-based models. Our results are important for assessing the impact of climate change on mountain infrastructure.
Joseph Fogarty, Elie Bou-Zeid, Mitchell Bushuk, and Linette Boisvert
The Cryosphere, 18, 4335–4354, https://doi.org/10.5194/tc-18-4335-2024, https://doi.org/10.5194/tc-18-4335-2024, 2024
Short summary
Short summary
We hypothesize that using a broad set of surface characterization metrics for polar sea ice surfaces will lead to more accurate representations in general circulation models. However, the first step is to identify the minimum set of metrics required. We show via numerical simulations that sea ice surface patterns can play a crucial role in determining boundary layer structures. We then statistically analyze a set of high-resolution sea ice surface images to obtain this minimal set of parameters.
Antonio Juarez-Martinez, Javier Blasco, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
The Cryosphere, 18, 4257–4283, https://doi.org/10.5194/tc-18-4257-2024, https://doi.org/10.5194/tc-18-4257-2024, 2024
Short summary
Short summary
We present sea level projections for Antarctica in the context of ISMIP6-2300 with several forcings but extend the simulations to 2500, showing that more than 3 m of sea level contribution could be reached. We also test the sensitivity on a basal melting parameter and determine the timing of the loss of ice in the west region. All the simulations were carried out with the ice sheet model Yelmo.
Daniel Moreno-Parada, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
The Cryosphere, 18, 4215–4232, https://doi.org/10.5194/tc-18-4215-2024, https://doi.org/10.5194/tc-18-4215-2024, 2024
Short summary
Short summary
Our study tries to understand how the ice temperature evolves in a large mass as in the case of Antarctica. We found a relation that tells us the ice temperature at any point. These results are important because they also determine how the ice moves. In general, ice moves due to slow deformation (as if pouring honey from a jar). Nevertheless, in some regions the ice base warms enough and melts. The liquid water then serves as lubricant and the ice slides and its velocity increases rapidly.
Tim Hageman, Jessica Mejía, Ravindra Duddu, and Emilio Martínez-Pañeda
The Cryosphere, 18, 3991–4009, https://doi.org/10.5194/tc-18-3991-2024, https://doi.org/10.5194/tc-18-3991-2024, 2024
Short summary
Short summary
Due to surface melting, meltwater lakes seasonally form on the surface of glaciers. These lakes drive hydrofractures that rapidly transfer water to the base of ice sheets. This paper presents a computational method to capture the complicated hydrofracturing process. Our work reveals that viscous ice rheology has a great influence on the short-term propagation of fractures, enabling fast lake drainage, whereas thermal effects (frictional heating, conduction, and freezing) have little influence.
Alex Cabaj, Paul J. Kushner, and Alek A. Petty
EGUsphere, https://doi.org/10.5194/egusphere-2024-2562, https://doi.org/10.5194/egusphere-2024-2562, 2024
Short summary
Short summary
The output of snow-on-sea-ice models is influenced by the choice of snowfall input used. We ran such a model with different snowfall inputs and calibrated it to observations, produced a new calibrated snow product, and regionally compared the model outputs to another snow-on-sea-ice model. The two models agree best on the seasonal cycle of snow in the central Arctic Ocean. However, estimated snow trends in some regions can depend more on the snowfall input than on the choice of model.
Clémence Herny, Pascal Hagenmuller, Guillaume Chambon, Isabel Peinke, and Jacques Roulle
The Cryosphere, 18, 3787–3805, https://doi.org/10.5194/tc-18-3787-2024, https://doi.org/10.5194/tc-18-3787-2024, 2024
Short summary
Short summary
This paper presents the evaluation of a numerical discrete element method (DEM) by simulating cone penetration tests in different snow samples. The DEM model demonstrated a good ability to reproduce the measured mechanical behaviour of the snow, namely the force evolution on the cone and the grain displacement field. Systematic sensitivity tests showed that the mechanical response depends not only on the microstructure of the sample but also on the mechanical parameters of grain contacts.
Louis Quéno, Rebecca Mott, Paul Morin, Bertrand Cluzet, Giulia Mazzotti, and Tobias Jonas
The Cryosphere, 18, 3533–3557, https://doi.org/10.5194/tc-18-3533-2024, https://doi.org/10.5194/tc-18-3533-2024, 2024
Short summary
Short summary
Snow redistribution by wind and avalanches strongly influences snow hydrology in mountains. This study presents a novel modelling approach to best represent these processes in an operational context. The evaluation of the simulations against airborne snow depth measurements showed remarkable improvement in the snow distribution in mountains of the eastern Swiss Alps, with a representation of snow accumulation and erosion areas, suggesting promising benefits for operational snow melt forecasts.
André Löfgren, Thomas Zwinger, Peter Råback, Christian Helanow, and Josefin Ahlkrona
The Cryosphere, 18, 3453–3470, https://doi.org/10.5194/tc-18-3453-2024, https://doi.org/10.5194/tc-18-3453-2024, 2024
Short summary
Short summary
This paper investigates a stabilization method for free-surface flows in the context of glacier simulations. Previous applications of the stabilization on ice flows have only considered simple ice-sheet benchmark problems; in particular the method had not been tested on real-world glacier domains. This work addresses this shortcoming by demonstrating that the stabilization works well also in this case and increases stability and robustness without negatively impacting computation times.
Ange Haddjeri, Matthieu Baron, Matthieu Lafaysse, Louis Le Toumelin, César Deschamps-Berger, Vincent Vionnet, Simon Gascoin, Matthieu Vernay, and Marie Dumont
The Cryosphere, 18, 3081–3116, https://doi.org/10.5194/tc-18-3081-2024, https://doi.org/10.5194/tc-18-3081-2024, 2024
Short summary
Short summary
Our study addresses the complex challenge of evaluating distributed alpine snow simulations with snow transport against snow depths from Pléiades stereo imagery and snow melt-out dates from Sentinel-2 and Landsat-8 satellites. Additionally, we disentangle error contributions between blowing snow, precipitation heterogeneity, and unresolved subgrid variability. Snow transport enhances the snow simulations at high elevations, while precipitation biases are the main error source in other areas.
Christopher Riedel and Jeffrey Anderson
The Cryosphere, 18, 2875–2896, https://doi.org/10.5194/tc-18-2875-2024, https://doi.org/10.5194/tc-18-2875-2024, 2024
Short summary
Short summary
Accurate sea ice conditions are crucial for quality sea ice projections, which have been connected to rapid warming over the Arctic. Knowing which observations to assimilate into models will help produce more accurate sea ice conditions. We found that not assimilating sea ice concentration led to more accurate sea ice states. The methods typically used to assimilate observations in our models apply assumptions to variables that are not well suited for sea ice because they are bounded variables.
Julia Selivanova, Doroteaciro Iovino, and Francesco Cocetta
The Cryosphere, 18, 2739–2763, https://doi.org/10.5194/tc-18-2739-2024, https://doi.org/10.5194/tc-18-2739-2024, 2024
Short summary
Short summary
Climate models show differences in sea ice representation in comparison to observations. Increasing the model resolution is a recognized way to improve model realism and obtain more reliable future projections. We find no strong impact of resolution on sea ice representation; it rather depends on the analysed variable and the model used. By 2050, the marginal ice zone (MIZ) becomes a dominant feature of the Arctic ice cover, suggesting a shift to a new regime similar to that in Antarctica.
Ailsa Chung, Frédéric Parrenin, Robert Mulvaney, Luca Vittuari, Massimo Frezzotti, Antonio Zanutta, David A. Lilien, Marie G. P. Cavitte, and Olaf Eisen
EGUsphere, https://doi.org/10.5194/egusphere-2024-1650, https://doi.org/10.5194/egusphere-2024-1650, 2024
Short summary
Short summary
We applied an ice flow model to a flow line from the summit of Dome C to the Beyond EPICA ice core drill site on Little Dome C in Antarctica. Results show that the oldest ice at the drill site may be 1.12 Ma (at age density of 20 kyr/m) and originate from around 15 km upstream. We also discuss the nature of the 200–250 m thick basal layer which could be composed of accreted ice, stagnant ice, or even disturbed ice containing debris.
Alexander A. Robel, Vincent Verjans, and Aminat A. Ambelorun
The Cryosphere, 18, 2613–2623, https://doi.org/10.5194/tc-18-2613-2024, https://doi.org/10.5194/tc-18-2613-2024, 2024
Short summary
Short summary
The average size of many glaciers and ice sheets changes when noise is added to the system. The reasons for this drift in glacier state is intrinsic to the dynamics of how ice flows and the bumpiness of the Earth's surface. We argue that not including noise in projections of ice sheet evolution over coming decades and centuries is a pervasive source of bias in these computer models, and so realistic variability in glacier and climate processes must be included in models.
John K. Dukowicz
EGUsphere, https://doi.org/10.5194/egusphere-2024-1052, https://doi.org/10.5194/egusphere-2024-1052, 2024
Short summary
Short summary
A novel transformation of the Stokes ice sheet equations is presented that expands the scope of traditional methods. The new formulation is closely related to a widely used Stokes approximation, the Blatter-Pattyn model, such that an ice sheet model may be easily switched between the two formulations, allowing for adaptive applications. The new formulation also facilitates new approximations that improve on the Blatter-Pattyn model, heretofore the best approximate ice sheet model.
José M. Muñoz-Hermosilla, Jaime Otero, Eva De Andrés, Kaian Shahateet, Francisco Navarro, and Iván Pérez-Doña
The Cryosphere, 18, 1911–1924, https://doi.org/10.5194/tc-18-1911-2024, https://doi.org/10.5194/tc-18-1911-2024, 2024
Short summary
Short summary
A large fraction of the mass loss from marine-terminating glaciers is attributed to frontal ablation. In this study, we used a 3D ice flow model of a real glacier that includes the effects of calving and submarine melting. Over a 30-month simulation, we found that the model reproduced the seasonal cycle for this glacier. Besides, the front positions were in good agreement with observations in the central part of the front, with longitudinal differences, on average, below 15 m.
Charlotte Durand, Tobias Sebastian Finn, Alban Farchi, Marc Bocquet, Guillaume Boutin, and Einar Ólason
The Cryosphere, 18, 1791–1815, https://doi.org/10.5194/tc-18-1791-2024, https://doi.org/10.5194/tc-18-1791-2024, 2024
Short summary
Short summary
This paper focuses on predicting Arctic-wide sea-ice thickness using surrogate modeling with deep learning. The model has a predictive power of 12 h up to 6 months. For this forecast horizon, persistence and daily climatology are systematically outperformed, a result of learned thermodynamics and advection. Consequently, surrogate modeling with deep learning proves to be effective at capturing the complex behavior of sea ice.
Mathieu Plante, Jean-François Lemieux, L. Bruno Tremblay, Adrienne Tivy, Joey Angnatok, François Roy, Gregory Smith, Frédéric Dupont, and Adrian K. Turner
The Cryosphere, 18, 1685–1708, https://doi.org/10.5194/tc-18-1685-2024, https://doi.org/10.5194/tc-18-1685-2024, 2024
Short summary
Short summary
We use a sea ice model to reproduce ice growth observations from two buoys deployed on coastal sea ice and analyze the improvements brought by new physics that represent the presence of saline liquid water in the ice interior. We find that the new physics with default parameters degrade the model performance, with overly rapid ice growth and overly early snow flooding on top of the ice. The performance is largely improved by simple modifications to the ice growth and snow-flooding algorithms.
Joshua Cuzzone, Matias Romero, and Shaun A. Marcott
The Cryosphere, 18, 1381–1398, https://doi.org/10.5194/tc-18-1381-2024, https://doi.org/10.5194/tc-18-1381-2024, 2024
Short summary
Short summary
We simulate the retreat history of the Patagonian Ice Sheet (PIS) across the Chilean Lake District from 22–10 ka. These results improve our understanding of the response of the PIS to deglacial warming and the patterns of deglacial ice margin retreat where gaps in the geologic record still exist, and they indicate that changes in large-scale precipitation during the last deglaciation played an important role in modulating the response of ice margin change across the PIS to deglacial warming.
Larissa van der Laan, Anouk Vlug, Adam A. Scaife, Fabien Maussion, and Kristian Förster
EGUsphere, https://doi.org/10.5194/egusphere-2024-387, https://doi.org/10.5194/egusphere-2024-387, 2024
Short summary
Short summary
Usually, glacier models are supplied with climate information from long (e.g. 100 year) simulations by global climate models. In this paper, we test the feasibility of supplying glacier models with shorter simulations, to get more accurate information on 5–10 year time scales. Reliable information on these time scales is very important, especially for water management experts to know how much meltwater to expect, for rivers, agriculture and drinking water.
Chao-Yuan Yang, Jiping Liu, and Dake Chen
The Cryosphere, 18, 1215–1239, https://doi.org/10.5194/tc-18-1215-2024, https://doi.org/10.5194/tc-18-1215-2024, 2024
Short summary
Short summary
We present a new atmosphere–ocean–wave–sea ice coupled model to study the influences of ocean waves on Arctic sea ice simulation. Our results show (1) smaller ice-floe size with wave breaking increases ice melt, (2) the responses in the atmosphere and ocean to smaller floe size partially reduce the effect of the enhanced ice melt, (3) the limited oceanic energy is a strong constraint for ice melt enhancement, and (4) ocean waves can indirectly affect sea ice through the atmosphere and the ocean.
Yurii Batrak, Bin Cheng, and Viivi Kallio-Myers
The Cryosphere, 18, 1157–1183, https://doi.org/10.5194/tc-18-1157-2024, https://doi.org/10.5194/tc-18-1157-2024, 2024
Short summary
Short summary
Atmospheric reanalyses provide consistent series of atmospheric and surface parameters in a convenient gridded form. In this paper, we study the quality of sea ice in a recently released regional reanalysis and assess its added value compared to a global reanalysis. We show that the regional reanalysis, having a more complex sea ice model, gives an improved representation of sea ice, although there are limitations indicating potential benefits in using more advanced approaches in the future.
Michael A. Rawlins and Ambarish V. Karmalkar
The Cryosphere, 18, 1033–1052, https://doi.org/10.5194/tc-18-1033-2024, https://doi.org/10.5194/tc-18-1033-2024, 2024
Short summary
Short summary
Flows of water, carbon, and materials by Arctic rivers are being altered by climate warming. We used simulations from a permafrost hydrology model to investigate future changes in quantities influencing river exports. By 2100 Arctic rivers will receive more runoff from the far north where abundant soil carbon can leach in. More water will enter them via subsurface pathways particularly in summer and autumn. An enhanced water cycle and permafrost thaw are changing river flows to coastal areas.
Cited articles
Akbari, G., Basirat Tabrizi, H., and Damangir, E.: Numerical and experimental investigation of variable phase transformation number effect in porous media during freezing process, Heat Mass Trans., 45, 407–416, 2009.
Bear, J.: Dynamics of Fluids in Porous Media, Dover Publications, Inc., New York, 2013.
Bense, V., Ferguson, G., and Kooi, H.: Evolution of shallow groundwater flow systems in areas of degrading permafrost, Geophys. Res. Lett., 36, L22401, https://doi.org/10.1029/2009GL039225, 2009.
Bosson, E., Selroos, J. O., Stigsson, M., Gustafsson, L. G., and Destouni, G.: Exchange and pathways of deep and shallow groundwater in different climate and permafrost conditions using the Forsmark site, Sweden, as an example catchment, Hydrogeol. J., 21, 225–237, 2013.
Dall'Amico, M., Endrizzi, S., Gruber, S., and Rigon, R.: A robust and energy-conserving model of freezing variably-saturated soil, The Cryosphere, 5, 469–484, https://doi.org/10.5194/tc-5-469-2011, 2011.
Edlefsen, N. and Anderson, A.: Thermodynamics of soil moisture, Hilgardia, 15, 31–298, 1943.
Frampton, A., Painter, S., Lyon, S., and Destouni, G.: Non-isothermal, three-phase simulations of near-surface flows in a model permafrost system under seasonal variability and climate change, J. Hydrol., 403, 352–359, https://doi.org/10.1016/j.jhydrol.2011.04.010, 2011.
Grenier, C., Régnier, D., Mouche, E., Benabderrahmane, H., Costard, F., and Davy, P.: Impact of permafrost development on groundwater flow patterns: a numerical study considering freezing cycles on a two-dimensional vertical cut through a generic river-plain system, Hydrogeol. J., 21, 257–270, 2013.
Grimm, R. and Painter, S.: On the secular evolution of groundwater on Mars, Geophys. Res. Lett., 36, L24803, https://doi.org/10.1029/2009GL041018, 2009.
Guymon, G. and Luthin, J.: A coupled heat and moisture transport model for arctic soils, Water Resour. Res., 10, 995–1001, 1974.
Hammond, G. E., Lichtner, P. C., and Mills, R. T.: Evaluating the performance of parallel subsurface simulators: An illustrative example with PFLOTRAN, Water Resour. Res., 50, 208–228, 2014.
Hansen, J., Ruedy, R., Glascoe, J., and Sato, M.: GISS analysis of surface temperature change, J. Geophys. Res., 104, 30997–31022, 1999.
Hansson, K., Simunek, J., Mizoguchi, M., Lundin, L, C., and van Genuchten, M.: Water flow and heat transport in frozen soil: numerical solution and freeze-thaw applications, Vadose Zone J., 3, 693–704, 2004.
Harlan, R.: Analysis of coupled heat-fluid transport in partially frozen soil, Water Resour. Res., 9, 1314–1323, 1973.
Jame, Y. and Norum, D.: Heat and mass transfer in a freezing unsaturated porous medium, Water Resour. Res., 16, 811–819, 1980.
Kane, D., Hinzman, L., and Zarling, J.: Thermal response of the active layer to climatic warming in a permafrost environment, Cold Reg. Sci. Technol., 19, 111–122, 1991.
Lawrence, D. M., Oleson, K. W., Flanner, M. G., Thornton, P. E., Swenson, S. C., Lawrence, P. J., Zeng, X., Yang, Z.-L., Levis, S., Sakaguchi, K., Bonan, G. B., and Slater, A. G.: Parameterization improvements and functional and structural advances in Version 4 of the Community Land Model, J. Adv. Model Earth Sy., 3, M03001, https://doi.org/10.1029/2011MS00045, 2011.
Lichtner, P. C., Hammond, G. E., Lu, C., Karra, S., Bisht, G., Andre, B., Mills, R. T., and Kumar, J.: PFLOTRAN User Manual, Tech. rep., http://www.pflotran.org/docs/user_manual.pdf (last access: October 2014), 2013.
Ling, F. and Zhang, T.: A numerical model for surface energy balance and thermal regime of the active layer and permafrost containing unfrozen water, Cold Reg. Sci. Technol., 38, 1–15, 2004.
Lu, T., Du, J., Lei, S., and Wang, B.: Heat and mass transfer in unsaturated porous media with solid–liquid change, Heat Mass Trans., 37, 237–242, 2001.
Lunardini, V.: Climatic warming and the degradation of warm permafrost, Permafrost Periglac., 7, 311–320, 1996.
McKenzie, J. M. and Voss, C. I.: Permafrost thaw in a nested groundwater-flow system, Hydrogeol. J., 21, 299–316, 2013.
McKenzie, J. M., Voss, C. I., and Siegel, D. I.: Groundwater flow with energy transport and water–ice phase change: numerical simulations, benchmarks, and application to freezing in peat bogs, Adv. Water Resour., 30, 966–983, 2007.
Mualem, Y.: A new model for predicting the hydraulic conductivity of unsaturated porous media, Water Resour. Res., 12, 513–322, 1976.
Nakano, Y. and Brown, J.: Effect of a freezing zone of finite width on the thermal regime of soils, Water Resour. Res., 7, 1226–1233, 1971.
Oleson, K., Lawrence, D., Bonan, G., Flanner, M., Kluzek, E., Lawrence, P., Levis, S., Swenson, S., Thornton, P., Dai, A., Decker, M., Dickinson, R., Feddema, J., Heald, C., Hoffman, F., Lamarque, J.-F., Mahowald, N., Niu, G.-Y., Qian, T., Randerson, J., Running, S., Sakaguchi, K., Slater, A., Stockli, R., Wang, A., Yang, Z.-L., Zeng, X., and Zeng, X.: Technical Description of version 4.0 of the Community Land Model (CLM), NCAR Technical Note NCAR/TN−478+STR, National Center for Atmospheric Research, Boulder, CO, 2010.
Osterkamp, T. and Romanovsky, V.: Evidence for warming and thawing of discontinuous permafrost in Alaska, Permafrost Periglac., 10, 17–37, 1999.
Painter, S.: Three-phase numerical model of water migration in partially frozen geological media: model formulation, validation, and applications, Comput. Geosci., 15, 69–85, 2011.
Painter, S. and Karra, S.: Constitutive model for unfrozen water content in subfreezing unsaturated soils, Vadose Zone J., 13, https://doi.org/10.2136/vzj2013.04.0071, 2014.
Painter, S., Moulton, J., and Wilson, C.: Modeling challenges for predicting hydrologic response to degrading permafrost, Hydrogeol. J., 21, 221–224, 2013.
Schuur, E., Bockheim, J., Canadell, J., Euskirchen, E., Field, C., Goryachkin, S., Hagemann, S., Kuhry, P., Lafleur, P., Lee, H., Mazhitova, G., Nelson, F. E., Rinke, A., Romanovsky, V. E., Shiklomanov, N., Tarnocai, C., Venevsky, S., Vogel, J. G., and Zimov, S. A.: Vulnerability of permafrost carbon to climate change: implications for the global carbon cycle, BioScience, 58, 701–714, 2008.
Sheshukov, A. and Nieber, J.: One-dimensional freezing of nonheaving unsaturated soils: model formulation and similarity solution, Water Resour. Res., 47, W11519, https://doi.org/10.1029/2011WR010512, 2011.
Tarnocai, C., Canadell, J., Schuur, E., Kuhry, P., Mazhitova, G., and Zimov, S.: Soil organic carbon pools in the northern circumpolar permafrost region, Global Biogeochem. Cy., 23, GB2023, https://doi.org/10.1029/2008GB003327, 2009.
Turner, J., Overland, J., and Walsh, J.: An Arctic and Antarctic perspective on recent climate change, Int. J. Climatol., 27, 277–293, 2007.
van Genuchten, M.: A closed-form equation for predicting the hydraulic conductivity of unsaturated soils, Soil Sci. Soc. Am. J., 44, 892–898, 1980.
Vidstrand, P., Follin, S., Selroos, J.-O., Näslund, J.-O., and Rhén, I.: Modeling of groundwater flow at depth in crystalline rock beneath a moving ice-sheet margin, exemplified by the Fennoscandian Shield, Sweden, Hydrogeol. J., 21, 239–255, 2013.
Watanabe, K. and Wake, T.: Measurement of unfrozen water content and relative permittivity of frozen unsaturated soil using NMR and TDR, Cold Reg. Sci. Technol., 59, 34–41, 2009.
Webb, S. W. and Ho, C. K.: Review of enhanced vapor diffusion in porous media, in: Proceedings of the TOUGH Workshop'98, Lawrence Berkeley National Laboratory Report LBNL-41995, LBNL technical report, 257–262, 1998.
Wen, Z., Ma, W., Feng, W., Deng, Y., Wang, D., Fan, Z., and Zhou, C.: Experimental study on unfrozen water content and soil matric potential of Qinghai-Tibetan silty clay, Environ. Earth Sci., 66, 1467–1476, 2012.
White, M.: Theory and numerical application of subsurface flow and transport for transient freezing conditions, in: Proceedings of the Fifteenth Annual Americal Geophysical Hydrology Days, 339–352, 1995.
White, M. and Oostrom, M.: STOMP Subsurface Transport Over Multiple Phases: Users Guide PNNL-15782, Pacific Northwest National Laboratory, Richland, 2006.
Zhang, Y., Carey, S., and Quinton, W.: Evaluation of the algorithms and parameterizations for ground thawing and freezing simulation in permafrost regions, J. Geophys. Res., 113, D17116, https://doi.org/10.1029/2007JD009343, 2008.
Zhao, L., Gray, D., and Male, D.: Numerical analysis of simultaneous heat and mass transfer during infiltration into frozen ground, J. Hydrol., 200, 345–363, 1997.
Zhou, Y. and Zhou, G.: Numerical simulation of coupled heat-fluid transport in freezing soils using finite volume method, Heat Mass Trans., 46, 989–998, 2010.
