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.
Sergi Molins, Benjamin Andre, Jeffrey Johnson, Glenn Hammond, Benjamin Sulman, Konstantin Lipnikov, Marcus Day, James Beisman, Daniil Svyatsky, Hang Deng, Peter Lichtner, Carl Steefel, and David Moulton
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-108, https://doi.org/10.5194/gmd-2024-108, 2024
Preprint under review for GMD
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
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
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
Quantifying the Buttressing Contribution of Sea Ice to Crane Glacier
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
Biases in ice sheet models from missing noise-induced drift
A new glacier thickness and bed map for Svalbard
Multi-physics ensemble modelling of Arctic tundra snowpack properties
A 3D glacier dynamics–line plume model to estimate the frontal ablation of Hansbreen, Svalbard
Sensitivity of Future Projections of the Wilkes Subglacial Basin Ice Sheet to Grounding Line Melt Parameterizations
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
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
Smoothed particle hydrodynamics implementation of the standard viscous–plastic sea-ice model and validation in simple idealized experiments
Coupled thermo–geophysical inversion for permafrost monitoring
Using specularity content to evaluate eight geothermal heat flow maps of Totten Glacier
Surging of a Hudson Strait-scale ice stream: subglacial hydrology matters but the process details mostly do not
Impact of the Nares Strait sea ice arches on the long-term stability of the Petermann Glacier ice shelf
Coupling between ice flow and subglacial hydrology enhances marine ice-sheet retreat
Regularization and L-curves in ice sheet inverse models: a case study in the Filchner–Ronne catchment
Quantifying the uncertainty in the Eurasian ice-sheet geometry at the Penultimate Glacial Maximum (Marine Isotope Stage 6)
Simulating ice segregation and thaw consolidation in permafrost environments with the CryoGrid community model
Reconciling ice dynamics and bed topography with a versatile and fast ice thickness inversion
The stability of present-day Antarctic grounding lines – Part 2: Onset of irreversible retreat of Amundsen Sea glaciers under current climate on centennial timescales cannot be excluded
The stability of present-day Antarctic grounding lines – Part 1: No indication of marine ice sheet instability in the current geometry
Phase-field models of floe fracture in sea ice
Exploring the ability of the variable-resolution Community Earth System Model to simulate cryospheric–hydrological variables in High Mountain Asia
Investigating the thermal state of permafrost with Bayesian inverse modeling of heat transfer
Modelling the development and decay of cryoconite holes in northwestern Greenland
The effect of partial dissolution on sea-ice chemical transport: a combined model–observational study using poly- and perfluoroalkylated substances (PFASs)
Deep learning subgrid-scale parametrisations for short-term forecasting of sea-ice dynamics with a Maxwell elasto-brittle rheology
Modelling ice mélange based on the viscous-plastic sea-ice rheology
Representation of soil hydrology in permafrost regions may explain large part of inter-model spread in simulated Arctic and subarctic climate
Impact of atmospheric forcing uncertainties on Arctic and Antarctic sea ice simulations in CMIP6 OMIP models
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.
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.
Richard Parsons, Sainan Sun, G. Hilmar Gudmundsson, Jan Wuite, and Thomas Nagler
EGUsphere, https://doi.org/10.5194/egusphere-2024-1499, https://doi.org/10.5194/egusphere-2024-1499, 2024
Short summary
Short summary
In 2022, 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 sea ice was attached to the terminus of the glacier, it could provide a resistive stress against the glacier’s ice-flow, slowing down the rate of ice discharge. We used numerical modelling to quantify this resistive stress and found that the sea ice provided significant support to Crane prior to its disintegration.
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.
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.
Ward van Pelt and Thomas Frank
EGUsphere, https://doi.org/10.5194/egusphere-2024-1525, https://doi.org/10.5194/egusphere-2024-1525, 2024
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 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 analyze the new bed and thickness maps, quantify the ice volume (6,800 km3), and compare against radar data and previous studies.
Georgina Jean Woolley, Nick Rutter, Leanne Wake, Vincent Vionnet, Chris Derksen, Richard Essery, Philip Marsh, Rosamund Tutton, Branden Walker, Matthieu Lafaysse, and David Pritchard
EGUsphere, https://doi.org/10.5194/egusphere-2024-1237, https://doi.org/10.5194/egusphere-2024-1237, 2024
Short summary
Short summary
Parameterisations of Arctic snow processes were implemented into the multi-physics ensemble version of SVS2-Crocus and evaluated using density and SSA measurements at an Arctic tundra site. Optimal combinations of parameterisations that improved the simulation of density and SSA were identified. Top performing ensemble members featured modifications that raise wind speeds to increase compaction in surface layers, prevent snowdrift and increase viscosity in basal layers.
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.
Yu Wang, Chen Zhao, Rupert Gladstone, Thomas Zwinger, Ben Galton-Fenzi, and Poul Christoffersen
EGUsphere, https://doi.org/10.5194/egusphere-2024-1005, https://doi.org/10.5194/egusphere-2024-1005, 2024
Short summary
Short summary
Our research delves into the future ice loss in Antarctica’s Wilkes Subglacial Basin (WSB) and its impact on sea level rise, focusing on how basal melt is implemented at the grounding line in ice flow models. According to our best model results, under high-emission scenarios, the WSB ice sheet could undergo massive and rapid retreat between 2200 and 2300, potentially raising global sea levels by up to 0.34 m by 2500.
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.
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.
Oreste Marquis, Bruno Tremblay, Jean-François Lemieux, and Mohammed Islam
The Cryosphere, 18, 1013–1032, https://doi.org/10.5194/tc-18-1013-2024, https://doi.org/10.5194/tc-18-1013-2024, 2024
Short summary
Short summary
We developed a standard viscous–plastic sea-ice model based on the numerical framework called smoothed particle hydrodynamics. The model conforms to the theory within an error of 1 % in an idealized ridging experiment, and it is able to simulate stable ice arches. However, the method creates a dispersive plastic wave speed. The framework is efficient to simulate fractures and can take full advantage of parallelization, making it a good candidate to investigate sea-ice material properties.
Soňa Tomaškovičová and Thomas Ingeman-Nielsen
The Cryosphere, 18, 321–340, https://doi.org/10.5194/tc-18-321-2024, https://doi.org/10.5194/tc-18-321-2024, 2024
Short summary
Short summary
We present the results of a fully coupled modeling framework for simulating the ground thermal regime using only surface measurements to calibrate the thermal model. The heat conduction model is forced by surface ground temperature measurements and calibrated using the field measurements of time lapse apparent electrical resistivity. The resistivity-calibrated thermal model achieves a performance comparable to the traditional calibration of borehole temperature measurements.
Yan Huang, Liyun Zhao, Michael Wolovick, Yiliang Ma, and John C. Moore
The Cryosphere, 18, 103–119, https://doi.org/10.5194/tc-18-103-2024, https://doi.org/10.5194/tc-18-103-2024, 2024
Short summary
Short summary
Geothermal heat flux (GHF) is an important factor affecting the basal thermal environment of an ice sheet and crucial for its dynamics. But it is poorly defined for the Antarctic ice sheet. We simulate the basal temperature and basal melting rate with eight different GHF datasets. We use specularity content as a two-sided constraint to discriminate between local wet or dry basal conditions. Two medium-magnitude GHF distribution maps rank well, showing that most of the inland bed area is frozen.
Matthew Drew and Lev Tarasov
The Cryosphere, 17, 5391–5415, https://doi.org/10.5194/tc-17-5391-2023, https://doi.org/10.5194/tc-17-5391-2023, 2023
Short summary
Short summary
The interaction of fast-flowing regions of continental ice sheets with their beds governs how quickly they slide and therefore flow. The coupling of fast ice to its bed is controlled by the pressure of meltwater at its base. It is currently poorly understood how the physical details of these hydrologic systems affect ice speedup. Using numerical models we find, surprisingly, that they largely do not, except for the duration of the surge. This suggests that cheap models are sufficient.
Abhay Prakash, Qin Zhou, Tore Hattermann, and Nina Kirchner
The Cryosphere, 17, 5255–5281, https://doi.org/10.5194/tc-17-5255-2023, https://doi.org/10.5194/tc-17-5255-2023, 2023
Short summary
Short summary
Sea ice arch formation in the Nares Strait has shielded the Petermann Glacier ice shelf from enhanced basal melting. However, with the sustained decline of the Arctic sea ice predicted to continue, the ice shelf is likely to be exposed to a year-round mobile and thin sea ice cover. In such a scenario, our modelled results show that elevated temperatures, and more importantly, a stronger ocean circulation in the ice shelf cavity, could result in up to two-thirds increase in basal melt.
George Lu and Jonathan Kingslake
EGUsphere, https://doi.org/10.5194/egusphere-2023-2794, https://doi.org/10.5194/egusphere-2023-2794, 2023
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.
Michael Wolovick, Angelika Humbert, Thomas Kleiner, and Martin Rückamp
The Cryosphere, 17, 5027–5060, https://doi.org/10.5194/tc-17-5027-2023, https://doi.org/10.5194/tc-17-5027-2023, 2023
Short summary
Short summary
The friction underneath ice sheets can be inferred from observed velocity at the top, but this inference requires smoothing. The selection of smoothing has been highly variable in the literature. Here we show how to rigorously select the best smoothing, and we show that the inferred friction converges towards the best knowable field as model resolution improves. We use this to learn about the best description of basal friction and to formulate recommended best practices for other modelers.
Oliver G. Pollard, Natasha L. M. Barlow, Lauren J. Gregoire, Natalya Gomez, Víctor Cartelle, Jeremy C. Ely, and Lachlan C. Astfalck
The Cryosphere, 17, 4751–4777, https://doi.org/10.5194/tc-17-4751-2023, https://doi.org/10.5194/tc-17-4751-2023, 2023
Short summary
Short summary
We use advanced statistical techniques and a simple ice-sheet model to produce an ensemble of plausible 3D shapes of the ice sheet that once stretched across northern Europe during the previous glacial maximum (140,000 years ago). This new reconstruction, equivalent in volume to 48 ± 8 m of global mean sea-level rise, will improve the interpretation of high sea levels recorded from the Last Interglacial period (120 000 years ago) that provide a useful perspective on the future.
Juditha Aga, Julia Boike, Moritz Langer, Thomas Ingeman-Nielsen, and Sebastian Westermann
The Cryosphere, 17, 4179–4206, https://doi.org/10.5194/tc-17-4179-2023, https://doi.org/10.5194/tc-17-4179-2023, 2023
Short summary
Short summary
This study presents a new model scheme for simulating ice segregation and thaw consolidation in permafrost environments, depending on ground properties and climatic forcing. It is embedded in the CryoGrid community model, a land surface model for the terrestrial cryosphere. We describe the model physics and functionalities, followed by a model validation and a sensitivity study of controlling factors.
Thomas Frank, Ward J. J. van Pelt, and Jack Kohler
The Cryosphere, 17, 4021–4045, https://doi.org/10.5194/tc-17-4021-2023, https://doi.org/10.5194/tc-17-4021-2023, 2023
Short summary
Short summary
Since the ice thickness of most glaciers worldwide is unknown, and since it is not feasible to visit every glacier and observe their thickness directly, inverse modelling techniques are needed that can calculate ice thickness from abundant surface observations. Here, we present a new method for doing that. Our methodology relies on modelling the rate of surface elevation change for a given glacier, compare this with observations of the same quantity and change the bed until the two are in line.
Ronja Reese, Julius Garbe, Emily A. Hill, Benoît Urruty, Kaitlin A. Naughten, Olivier Gagliardini, Gaël Durand, Fabien Gillet-Chaulet, G. Hilmar Gudmundsson, David Chandler, Petra M. Langebroek, and Ricarda Winkelmann
The Cryosphere, 17, 3761–3783, https://doi.org/10.5194/tc-17-3761-2023, https://doi.org/10.5194/tc-17-3761-2023, 2023
Short summary
Short summary
We use an ice sheet model to test where current climate conditions in Antarctica might lead. We find that present-day ocean and atmosphere conditions might commit an irreversible collapse of parts of West Antarctica which evolves over centuries to millennia. Importantly, this collapse is not irreversible yet.
Emily A. Hill, Benoît Urruty, Ronja Reese, Julius Garbe, Olivier Gagliardini, Gaël Durand, Fabien Gillet-Chaulet, G. Hilmar Gudmundsson, Ricarda Winkelmann, Mondher Chekki, David Chandler, and Petra M. Langebroek
The Cryosphere, 17, 3739–3759, https://doi.org/10.5194/tc-17-3739-2023, https://doi.org/10.5194/tc-17-3739-2023, 2023
Short summary
Short summary
The grounding lines of the Antarctic Ice Sheet could enter phases of irreversible retreat or advance. We use three ice sheet models to show that the present-day locations of Antarctic grounding lines are reversible with respect to a small perturbation away from their current position. This indicates that present-day retreat of the grounding lines is not yet irreversible or self-enhancing.
Huy Dinh, Dimitrios Giannakis, Joanna Slawinska, and Georg Stadler
The Cryosphere, 17, 3883–3893, https://doi.org/10.5194/tc-17-3883-2023, https://doi.org/10.5194/tc-17-3883-2023, 2023
Short summary
Short summary
We develop a numerical method to simulate the fracture in kilometer-sized chunks of floating ice in the ocean. Our approach uses a mathematical model that balances deformation energy against the energy required for fracture. We study the strength of ice chunks that contain random impurities due to prior damage or refreezing and what types of fractures are likely to occur. Our model shows that crack direction critically depends on the orientation of impurities relative to surrounding forces.
René R. Wijngaard, Adam R. Herrington, William H. Lipscomb, Gunter R. Leguy, and Soon-Il An
The Cryosphere, 17, 3803–3828, https://doi.org/10.5194/tc-17-3803-2023, https://doi.org/10.5194/tc-17-3803-2023, 2023
Short summary
Short summary
We evaluate the ability of the Community Earth System Model (CESM2) to simulate cryospheric–hydrological variables, such as glacier surface mass balance (SMB), over High Mountain Asia (HMA) by using a global grid (~111 km) with regional refinement (~7 km) over HMA. Evaluations of two different simulations show that climatological biases are reduced, and glacier SMB is improved (but still too negative) by modifying the snow and glacier model and using an updated glacier cover dataset.
Brian Groenke, Moritz Langer, Jan Nitzbon, Sebastian Westermann, Guillermo Gallego, and Julia Boike
The Cryosphere, 17, 3505–3533, https://doi.org/10.5194/tc-17-3505-2023, https://doi.org/10.5194/tc-17-3505-2023, 2023
Short summary
Short summary
It is now well known from long-term temperature measurements that Arctic permafrost, i.e., ground that remains continuously frozen for at least 2 years, is warming in response to climate change. Temperature, however, only tells half of the story. In this study, we use computer modeling to better understand how the thawing and freezing of water in the ground affects the way permafrost responds to climate change and what temperature trends can and cannot tell us about how permafrost is changing.
Yukihiko Onuma, Koji Fujita, Nozomu Takeuchi, Masashi Niwano, and Teruo Aoki
The Cryosphere, 17, 3309–3328, https://doi.org/10.5194/tc-17-3309-2023, https://doi.org/10.5194/tc-17-3309-2023, 2023
Short summary
Short summary
We established a novel model that simulates the temporal changes in cryoconite hole (CH) depth using heat budgets calculated independently at the ice surface and CH bottom based on hole shape geometry. The simulations suggest that CH depth is governed by the balance between the intensity of the diffuse component of downward shortwave radiation and the wind speed. The meteorological conditions may be important factors contributing to the recent ice surface darkening via the redistribution of CHs.
Max Thomas, Briana Cate, Jack Garnett, Inga J. Smith, Martin Vancoppenolle, and Crispin Halsall
The Cryosphere, 17, 3193–3201, https://doi.org/10.5194/tc-17-3193-2023, https://doi.org/10.5194/tc-17-3193-2023, 2023
Short summary
Short summary
A recent study showed that pollutants can be enriched in growing sea ice beyond what we would expect from a perfectly dissolved chemical. We hypothesise that this effect is caused by the specific properties of the pollutants working in combination with fluid moving through the sea ice. To test our hypothesis, we replicate this behaviour in a sea-ice model and show that this type of modelling can be applied to predicting the transport of chemicals with complex behaviour in sea ice.
Tobias Sebastian Finn, Charlotte Durand, Alban Farchi, Marc Bocquet, Yumeng Chen, Alberto Carrassi, and Véronique Dansereau
The Cryosphere, 17, 2965–2991, https://doi.org/10.5194/tc-17-2965-2023, https://doi.org/10.5194/tc-17-2965-2023, 2023
Short summary
Short summary
We combine deep learning with a regional sea-ice model to correct model errors in the sea-ice dynamics of low-resolution forecasts towards high-resolution simulations. The combined model improves the forecast by up to 75 % and thereby surpasses the performance of persistence. As the error connection can additionally be used to analyse the shortcomings of the forecasts, this study highlights the potential of combined modelling for short-term sea-ice forecasting.
Saskia Kahl, Carolin Mehlmann, and Dirk Notz
EGUsphere, https://doi.org/10.5194/egusphere-2023-982, https://doi.org/10.5194/egusphere-2023-982, 2023
Short summary
Short summary
Ice mélange is a mixture of sea ice and icebergs, which can have a strong influence on the sea-ice-ocean interaction. So far, ice mélange is not represented in climate models. We include icebergs into the most used sea-ice model by modifying the mathematical equations that describe the material law of sea ice. We show with three test cases that the modification is necessary to represent icebergs. Furthermore we suggest a numerical method to solve the ice mélange equations computational efficient.
Philipp de Vrese, Goran Georgievski, Jesus Fidel Gonzalez Rouco, Dirk Notz, Tobias Stacke, Norman Julius Steinert, Stiig Wilkenskjeld, and Victor Brovkin
The Cryosphere, 17, 2095–2118, https://doi.org/10.5194/tc-17-2095-2023, https://doi.org/10.5194/tc-17-2095-2023, 2023
Short summary
Short summary
The current generation of Earth system models exhibits large inter-model differences in the simulated climate of the Arctic and subarctic zone. We used an adapted version of the Max Planck Institute (MPI) Earth System Model to show that differences in the representation of the soil hydrology in permafrost-affected regions could help explain a large part of this inter-model spread and have pronounced impacts on important elements of Earth systems as far to the south as the tropics.
Xia Lin, François Massonnet, Thierry Fichefet, and Martin Vancoppenolle
The Cryosphere, 17, 1935–1965, https://doi.org/10.5194/tc-17-1935-2023, https://doi.org/10.5194/tc-17-1935-2023, 2023
Short summary
Short summary
This study provides clues on how improved atmospheric reanalysis products influence sea ice simulations in ocean–sea ice models. The summer ice concentration simulation in both hemispheres can be improved with changed surface heat fluxes. The winter Antarctic ice concentration and the Arctic drift speed near the ice edge and the ice velocity direction simulations are improved with changed wind stress. The radiation fluxes and winds in atmospheric reanalyses are crucial for sea ice simulations.
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.
