Articles | Volume 16, issue 3
https://doi.org/10.5194/tc-16-903-2022
© Author(s) 2022. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/tc-16-903-2022
© Author(s) 2022. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Elements of future snowpack modeling – Part 1: A physical instability arising from the nonlinear coupling of transport and phase changes
Konstantin Schürholt
WSL Institute for Snow and Avalanche Research SLF, Flüelastr. 11, 7260 Davos, Switzerland
Julia Kowalski
AICES Graduate School, RWTH Aachen University, Schinkelstr. 2a, 52062 Aachen, Germany
Computational Geoscience, University of Göttingen, Goldschmidtstr. 1, 37077 Göttingen, Germany
WSL Institute for Snow and Avalanche Research SLF, Flüelastr. 11, 7260 Davos, Switzerland
Related authors
No articles found.
V Mithlesh Kumar, Anil Yildiz, and Julia Kowalski
EGUsphere, https://doi.org/10.5194/egusphere-2025-4531, https://doi.org/10.5194/egusphere-2025-4531, 2025
This preprint is open for discussion and under review for Nonlinear Processes in Geophysics (NPG).
Short summary
Short summary
The reliability of Bayesian calibration depends on the quality and availability of observational data. But are we choosing the right data? We address this question by measuring the information gained during calibration to quantify how data selection influences the Bayesian calibration of physics-based landslide runout models. We find that more data does not always yield better results – observations that capture the dynamics governed by a parameter are more effective for its calibration.
Alexander Hermanns, Anne Caroline Lange, Julia Kowalski, Hendrik Fuchs, and Philipp Franke
EGUsphere, https://doi.org/10.5194/egusphere-2025-450, https://doi.org/10.5194/egusphere-2025-450, 2025
Short summary
Short summary
For air quality analyses, data assimilation models split available data into assimilation and validation data sets. The former is used to generate the analysis, the latter to verify the simulations. A preprocessor classifying the observations by the data characteristics is developed based on clustering algorithms. The assimilation and validation data sets are compiled by equally allocating data of each cluster. The resulting improvement of the analysis is evaluated with EURAD-IM.
Matthias Rauter and Julia Kowalski
Geosci. Model Dev., 17, 6545–6569, https://doi.org/10.5194/gmd-17-6545-2024, https://doi.org/10.5194/gmd-17-6545-2024, 2024
Short summary
Short summary
Snow avalanches can form large powder clouds that substantially exceed the velocity and reach of the dense core. Only a few complex models exist to simulate this phenomenon, and the respective hazard is hard to predict. This work provides a novel flow model that focuses on simple relations while still encapsulating the significant behaviour. The model is applied to reconstruct two catastrophic powder snow avalanche events in Austria.
Benjamin Walter, Hagen Weigel, Sonja Wahl, and Henning Löwe
The Cryosphere, 18, 3633–3652, https://doi.org/10.5194/tc-18-3633-2024, https://doi.org/10.5194/tc-18-3633-2024, 2024
Short summary
Short summary
The topmost layer of a snowpack forms the interface to the atmosphere and is critical for the reflectance of solar radiation and avalanche formation. The effect of wind on the surface snow microstructure during precipitation events is poorly understood and quantified. We performed controlled lab experiments in a ring wind tunnel to systematically quantify the snow microstructure for different wind speeds, temperatures and precipitation intensities and to identify the relevant processes.
Kévin Fourteau, Johannes Freitag, Mika Malinen, and Henning Löwe
The Cryosphere, 18, 2831–2846, https://doi.org/10.5194/tc-18-2831-2024, https://doi.org/10.5194/tc-18-2831-2024, 2024
Short summary
Short summary
Understanding the settling of snow under its own weight has applications from avalanche forecasts to ice core interpretations. We study how this settling can be modeled using 3D images of the internal structure of snow and ice deformation mechanics. We found that classical ice mechanics, as used, for instance, in glacier flow, explain the compaction of dense polar snow but not that of lighter seasonal snow. How, exactly, the ice deforms during light snow compaction thus remains an open question.
Julien Meloche, Melody Sandells, Henning Löwe, Nick Rutter, Richard Essery, Ghislain Picard, Randall K. Scharien, Alexandre Langlois, Matthias Jaggi, Josh King, Peter Toose, Jérôme Bouffard, Alessandro Di Bella, and Michele Scagliola
EGUsphere, https://doi.org/10.5194/egusphere-2024-1583, https://doi.org/10.5194/egusphere-2024-1583, 2024
Preprint archived
Short summary
Short summary
Sea ice thickness is essential for climate studies. Radar altimetry has provided sea ice thickness measurement, but uncertainty arises from interaction of the signal with the snow cover. Therefore, modelling the signal interaction with the snow is necessary to improve retrieval. A radar model was used to simulate the radar signal from the snow-covered sea ice. This work paved the way to improved physical algorithm to retrieve snow depth and sea ice thickness for radar altimeter missions.
Anna Braun, Kévin Fourteau, and Henning Löwe
The Cryosphere, 18, 1653–1668, https://doi.org/10.5194/tc-18-1653-2024, https://doi.org/10.5194/tc-18-1653-2024, 2024
Short summary
Short summary
The specific surface of snow dictates key physical properties and continuously evolves in natural snowpacks. This is referred to as metamorphism. This work develops a rigorous physical model for this evolution, which is able to reproduce X-ray tomography measurements without using unphysical tuning parameters. Our results emphasize that snow crystal growth at the micrometer scale ultimately controls the pace of metamorphism.
Kavitha Sundu, Johannes Freitag, Kévin Fourteau, and Henning Löwe
The Cryosphere, 18, 1579–1596, https://doi.org/10.5194/tc-18-1579-2024, https://doi.org/10.5194/tc-18-1579-2024, 2024
Short summary
Short summary
Ice crystals often show a rod-like, vertical orientation in snow and firn; they are said to be anisotropic. The stiffness in the vertical direction therefore differs from the horizontal, which, for example, impacts the propagation of seismic waves. To quantify this anisotropy, we conducted finite-element simulations of 391 snow, firn, and ice core microstructures obtained from X-ray tomography. We then derived a parameterization that may be employed for advanced seismic studies in polar regions.
Amy R. Macfarlane, Henning Löwe, Lucille Gimenes, David N. Wagner, Ruzica Dadic, Rafael Ottersberg, Stefan Hämmerle, and Martin Schneebeli
The Cryosphere, 17, 5417–5434, https://doi.org/10.5194/tc-17-5417-2023, https://doi.org/10.5194/tc-17-5417-2023, 2023
Short summary
Short summary
Snow acts as an insulating blanket on Arctic sea ice, keeping the underlying ice "warm", relative to the atmosphere. Knowing the snow's thermal conductivity is essential for understanding winter ice growth. During the MOSAiC expedition, we measured the thermal conductivity of snow. We found spatial and vertical variability to overpower any temporal variability or dependency on underlying ice type and the thermal resistance to be directly influenced by snow height.
Julien Brondex, Kévin Fourteau, Marie Dumont, Pascal Hagenmuller, Neige Calonne, François Tuzet, and Henning Löwe
Geosci. Model Dev., 16, 7075–7106, https://doi.org/10.5194/gmd-16-7075-2023, https://doi.org/10.5194/gmd-16-7075-2023, 2023
Short summary
Short summary
Vapor diffusion is one of the main processes governing snowpack evolution, and it must be accounted for in models. Recent attempts to represent vapor diffusion in numerical models have faced several difficulties regarding computational cost and mass and energy conservation. Here, we develop our own finite-element software to explore numerical approaches and enable us to overcome these difficulties. We illustrate the capability of these approaches on established numerical benchmarks.
Ghislain Picard, Henning Löwe, and Christian Mätzler
The Cryosphere, 16, 3861–3866, https://doi.org/10.5194/tc-16-3861-2022, https://doi.org/10.5194/tc-16-3861-2022, 2022
Short summary
Short summary
Microwave satellite observations used to monitor the cryosphere require radiative transfer models for their interpretation. These models represent how microwaves are scattered by snow and ice. However no existing theory is suitable for all types of snow and ice found on Earth. We adapted a recently published generic scattering theory to snow and show how it may improve the representation of snows with intermediate densities (~500 kg/m3) and/or with coarse grains at high microwave frequencies.
Clara Betancourt, Timo T. Stomberg, Ann-Kathrin Edrich, Ankit Patnala, Martin G. Schultz, Ribana Roscher, Julia Kowalski, and Scarlet Stadtler
Geosci. Model Dev., 15, 4331–4354, https://doi.org/10.5194/gmd-15-4331-2022, https://doi.org/10.5194/gmd-15-4331-2022, 2022
Short summary
Short summary
Ozone is a toxic greenhouse gas with high spatial variability. We present a machine-learning-based ozone-mapping workflow generating a transparent and reliable product. Going beyond standard mapping methods, this work combines explainable machine learning with uncertainty assessment to increase the integrity of the produced map.
Anna Simson, Henning Löwe, and Julia Kowalski
The Cryosphere, 15, 5423–5445, https://doi.org/10.5194/tc-15-5423-2021, https://doi.org/10.5194/tc-15-5423-2021, 2021
Short summary
Short summary
This companion paper deals with numerical particularities of partial differential equations underlying one-dimensional snow models. In this second part we include mechanical settling and develop a new hybrid (Eulerian–Lagrangian) method for solving the advection-dominated ice mass conservation on a moving mesh alongside Eulerian diffusion (heat and vapor) and phase changes. The scheme facilitates a modular and extendable solver strategy while retaining controls on numerical accuracy.
Sebastian Hellmann, Melchior Grab, Johanna Kerch, Henning Löwe, Andreas Bauder, Ilka Weikusat, and Hansruedi Maurer
The Cryosphere, 15, 3507–3521, https://doi.org/10.5194/tc-15-3507-2021, https://doi.org/10.5194/tc-15-3507-2021, 2021
Short summary
Short summary
In this study, we analyse whether ultrasonic measurements on ice core samples could be employed to derive information about the particular ice crystal orientation in these samples. We discuss if such ultrasonic scans of ice core samples could provide similarly detailed results as the established methods, which usually destroy the ice samples. Our geophysical approach is minimally invasive and could support the existing methods with additional and (semi-)continuous data points along the ice core.
Cited articles
Alnæs, M. S., Blechta, J., Hake, J., Johansson, A., Kehlet, B., Logg, A., Richardson, C., Ring, J., Rognes, M. E., and Wells, G. N.: The FEniCS Project Version 1.5, Archive of Numerical Software, 3, 9–23, https://doi.org/10.11588/ans.2015.100.20553, 2015. a, b, c
Barrere, M., Domine, F., Decharme, B., Morin, S., Vionnet, V., and Lafaysse, M.: Evaluating the performance of coupled snow–soil models in SURFEXv8 to simulate the permafrost thermal regime at a high Arctic site, Geosci. Model Dev., 10, 3461–3479, https://doi.org/10.5194/gmd-10-3461-2017, 2017. a
Brun, E., Martin, E., Simon, V., Gendre, C., and Coleou, C.: An Energy and Mass Model of Snow Cover Suitable for Operational Avalanche Forecasting, J. Glaciol., 35, 333–342, https://doi.org/10.3189/S0022143000009254, 1989. a
Calonne, N., Flin, F., Morin, S., Lesaffre, B., du Roscoat, S. R., and Geindreau, C.: Numerical and experimental investigations of the effective thermal conductivity of snow, Geophys. Res. Lett., 38, L23501, https://doi.org/10.1029/2011GL049234, 2011. a
Colbeck, S. C. and Jamieson, J.: The formation of faceted layers above crusts, Cold Reg. Sci. Technol., 33, 247–252, https://doi.org/10.1016/S0165-232X(01)00045-3, 2001.
a
Cross, M. C. and Hohenberg, P. C.: Pattern formation outside of equilibrium, Rev. Mod. Phys., 65, 851–1112, https://doi.org/10.1103/RevModPhys.65.851, 1993. a
Cummings, E. M.: Modeling the cryosphere with FEniCS, arXiv [preprint], arXiv:1609.02190, 9 September 2016. a
Domine, F., Barrere, M., and Sarrazin, D.: Seasonal evolution of the effective thermal conductivity of the snow and the soil in high Arctic herb tundra at Bylot Island, Canada, The Cryosphere, 10, 2573–2588, https://doi.org/10.5194/tc-10-2573-2016, 2016. a, b
Du, Q., Kamm, J. R., Lehoucq, R. B., and Parks, M. L.: A new approach for a nonlocal, nonlinear conservation law, SIAM J. Appl. Math., 72, 464–487, https://doi.org/10.1137/110833233, 2012. a
Fourteau, K., Domine, F., and Hagenmuller, P.: Macroscopic water vapor diffusion is not enhanced in snow, The Cryosphere, 15, 389–406, https://doi.org/10.5194/tc-15-389-2021, 2021. a
Hammonds, K., Lieb-Lappen, R., Baker, I., and Wang, X.: Investigating the thermophysical properties of the ice–snow interface under a controlled temperature gradient: Part I: Experiments & Observations, Cold Reg. Sci. Technol., 120, 157–167, https://doi.org/10.1016/j.coldregions.2015.09.006 2015. a, b
Jafari, M., Gouttevin, I., Couttet, M., Wever, N., Michel, A., Sharma, V., Rossmann, L., Maass, N., Nicolaus, M., and Lehning, M.: The Impact of Diffusive Water Vapor Transport on Snow Profiles in Deep and Shallow Snow Covers and on Sea Ice, Front. Earth Sci., 8, 249, https://doi.org/10.3389/feart.2020.00249, 2020. a, b, c, d
Krol, Q. and Löwe, H.: Upscaling ice crystal growth dynamics in snow: Rigorous modeling and comparison to 4D X-ray tomography data, Acta Mater., 151, 478–487, https://doi.org/10.1016/j.actamat.2018.03.010, 2018. a, b
Lehning, M., Bartelt, P., Brown, B., Fierz, C., and Satyawali, P.: A physical snowpack model for the Swiss avalanche warning: Part II. Snow microstructure, Cold Reg. Sci. Technol., 35, 147–167, https://doi.org/10.1016/S0165-232X(02)00073-3, 2002. a
Li, J. and Zwally, H. J.: Modeling the density variation in the shallow firn layer, Ann. Glaciol., 38, 309–313, https://doi.org/10.3189/172756404781814988, 2004. a, b
Libbrecht, K. G.: Physical properties of ice, SnowCrystals.com, http://www.cco.caltech.edu/~atomic/snowcrystals/ice/ice.htm (last access: 29 July 2021), 1999. a
Morland, L. W., Kelly, R. J., and Morris, E. M.: A mixture theory for a phase-changing snowpack, Cold Reg. Sci. Technol., 17, 271–285,
https://doi.org/10.1016/S0165-232X(05)80006-0, 1990. a
Schürholt, K., Kowalski, J., and Löwe, H.: A numerical solver for heat and mass transport in snow based on FEniCS, EnviDat [code], https://doi.org/10.16904/envidat.298, 2022. a
Shu, C.-W. and Osher, S.: Efficient implementation of essentially non-oscillatory shock-capturing schemes, J. Comput. Phys., 77, 439–471, https://doi.org/10.1016/0021-9991(88)90177-5, 1988. a
Simson, A., Löwe, H., and Kowalski, J.: Elements of future snowpack modeling – Part 2: A modular and extendable Eulerian–Lagrangian numerical scheme for coupled transport, phase changes and settling processes, The Cryosphere, 15, 5423–5445, https://doi.org/10.5194/tc-15-5423-2021, 2021. a, b
Sturm, M. and Benson, C. S.: Vapor transport, grain growth and depth-hoar development in the subarctic snow, J. Glaciol., 43, 42–59, https://doi.org/10.3189/S0022143000002793, 1997. a
Touzeau, A., Landais, A., Morin, S., Arnaud, L., and Picard, G.: Numerical experiments on vapor diffusion in polar snow and firn and its impact on isotopes using the multi-layer energy balance model Crocus in SURFEX v8.0, Geosci. Model Dev., 11, 2393–2418, https://doi.org/10.5194/gmd-11-2393-2018, 2018.
a, b
Vidal-Henriquez, E., Zykov, V., Bodenschatz, E., and Gholami, A.: Convective instability and boundary driven oscillations in a reaction-diffusion-advection model, Chaos, 27, 103110, https://doi.org/10.1063/1.4986153, 2017. a
Wever, N., Fierz, C., Mitterer, C., Hirashima, H., and Lehning, M.: Solving Richards Equation for snow improves snowpack meltwater runoff estimations in detailed multi-layer snowpack model, The Cryosphere, 8, 257–274, https://doi.org/10.5194/tc-8-257-2014, 2014. a
Short summary
This companion paper deals with numerical particularities of partial differential equations underlying 1D snow models. In this first part we neglect mechanical settling and demonstrate that the nonlinear coupling between diffusive transport (heat and vapor), phase changes and ice mass conservation contains a wave instability that may be relevant for weak layer formation. Numerical requirements are discussed in view of the underlying homogenization scheme.
This companion paper deals with numerical particularities of partial differential equations...