Articles | Volume 15, issue 12
https://doi.org/10.5194/tc-15-5423-2021
© Author(s) 2021. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/tc-15-5423-2021
© Author(s) 2021. 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 2: A modular and extendable Eulerian–Lagrangian numerical scheme for coupled transport, phase changes and settling processes
Anna Simson
CORRESPONDING AUTHOR
AICES Graduate School, RWTH Aachen University, Schinkelstr. 2a, 52062 Aachen, Germany
Methods for Model-based Development in Computational Engineering, RWTH Aachen University, Eilfschornsteinstraße 18, 52062 Aachen, Germany
Henning Löwe
WSL Institute for Snow and Avalanche Research SLF, Flüelastr. 11, 7260 Davos, Switzerland
Julia Kowalski
Computational Geoscience, University of Göttingen, Goldschmidtstr. 1, 37077 Göttingen, Germany
Methods for Model-based Development in Computational Engineering, RWTH Aachen University, Eilfschornsteinstraße 18, 52062 Aachen, Germany
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.
Konstantin Schürholt, Julia Kowalski, and Henning Löwe
The Cryosphere, 16, 903–923, https://doi.org/10.5194/tc-16-903-2022, https://doi.org/10.5194/tc-16-903-2022, 2022
Short summary
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.
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
Audet, D. and Fowler, A.:
A mathematical model for compaction in sedimentary basins,
Geophys. J. Int.,
110, 577–590, https://doi.org/10.1111/j.1365-246x.1992.tb02093.x, 1992. a, b
Bartelt, P. and Christen, M.:
A computational procedure for instationary temperature-dependent snow creep,
Springer Berlin Heidelberg, https://doi.org/10.1007/BFb0104195, 2007. a, b
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, b, c
Brun, E., David, P., Sudul, M., and Brunot, G.:
A numerical model to simulate snow-cover stratigraphy for operational avalanche forecasting,
J. Glaciol.,
38, 13–22, https://doi.org/10.1017/s0022143000009552, 1992. 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
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
Domine, F., Picard, G., Morin, S., Barrere, M., Madore, J.-B., and Langlois, A.:
Major Issues in Simulating Some Arctic Snowpack Properties Using Current Detailed Snow Physics Models: Consequences for the Thermal Regime and Water Budget of Permafrost,
J. Adv. Model. Earth Sy.,
11, 34–44, https://doi.org/10.1029/2018ms001445, 2019. a, b
Farlow, S. J.:
Partial Differential Equations for Scientists and Engineers,
Dover Publications, Mineola, New York, USA,
reprint of the John Wiley & Sons, New York, USA, 1982 edition, 1993. a
Helfricht, K., Hartl, L., Koch, R., Marty, C., and Olefs, M.: Obtaining sub-daily new snow density from automated measurements in high mountain regions, Hydrol. Earth Syst. Sci., 22, 2655–2668, https://doi.org/10.5194/hess-22-2655-2018, 2018. a
Huang, K., Zhang, R., and van Genuchten, M. T.:
An Eulerian-Lagrangian approach with an adaptively corrected method of characteristics to simulate variably saturated water flow,
Water Resour. Res.,
30, 499–507, https://doi.org/10.1029/93WR02881, 1994. a
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, e, f
Johnson, J. B.: Snow Deformation, in: Encyclopedia of Snow, Ice and Glaciers, 2011
edition, edited by: Singh, V. P., Singh, P., and Haritashya, UK, Springer,
Dordrecht, the Netherlands, 1041–1045, https://doi.org/10.1007/978-90-481-2642-2_501, 2011. a, b
Judson, A. and Doesken, N.:
Density of Freshly Fallen Snow in the Central Rocky Mountains,
B. Am. Meteorol. Soc.,
81, 1577–1587, https://doi.org/10.1175/1520-0477(2000)081<1577:DOFFSI>2.3.CO;2, 2000. a
Kirchner, H. K., Michot, G., Narita, H., and Suzuki, T.:
Snow as a foam of ice: Plasticity, fracture and the brittle-to-ductile transition,
Philos. Mag. A,
81, 2161–2181, https://doi.org/10.1080/01418610108217141, 2001. a
Kowalski, J. and Torrilhon, M.:
Moment Approximations and Model Cascades for Shallow Flow,
Commun. Comput. Phys.,
25, 669–702, https://doi.org/10.4208/cicp.oa-2017-0263, 2019. a
Krinner, G., Derksen, C., Essery, R., Flanner, M., Hagemann, S., Clark, M., Hall, A., Rott, H., Brutel-Vuilmet, C., Kim, H., Ménard, C. B., Mudryk, L., Thackeray, C., Wang, L., Arduini, G., Balsamo, G., Bartlett, P., Boike, J., Boone, A., Chéruy, F., Colin, J., Cuntz, M., Dai, Y., Decharme, B., Derry, J., Ducharne, A., Dutra, E., Fang, X., Fierz, C., Ghattas, J., Gusev, Y., Haverd, V., Kontu, A., Lafaysse, M., Law, R., Lawrence, D., Li, W., Marke, T., Marks, D., Ménégoz, M., Nasonova, O., Nitta, T., Niwano, M., Pomeroy, J., Raleigh, M. S., Schaedler, G., Semenov, V., Smirnova, T. G., Stacke, T., Strasser, U., Svenson, S., Turkov, D., Wang, T., Wever, N., Yuan, H., Zhou, W., and Zhu, D.: ESM-SnowMIP: assessing snow models and quantifying snow-related climate feedbacks, Geosci. Model Dev., 11, 5027–5049, https://doi.org/10.5194/gmd-11-5027-2018, 2018. a, b
Krol, Q. and Löwe, H.:
Analysis of local ice crystal growth in snow,
J. Glaciol.,
62, 378–390, https://doi.org/10.1017/jog.2016.32, 2016. a
Lacroix, M. and Garon, A.:
Numerical Solution of phase change problems: an Eulerian-Lagrangian approach,
Numer. Heat Tr. B-Fund.,
19, 57–78, https://doi.org/10.1080/10407799208944922, 1992. a
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, b
LeVeque, R. J.: Finite Volume Methods for Hyperbolic Problems,
Cambridge Texts in Applied Mathematics, Cambridge University Press, Cambridge, https://doi.org/10.1017/CBO9780511791253, 2002. a
Lundin, J. M., Stevens, C. M., Arthern, R., Buizert, C., Orsi, A., Ligtenberg, S. R., Simonsen, S. B., Cummings, E., Essery, R., Leahy, W., Harris, P., Helsen, M. M., and Waddington, E. D.:
Firn Model Intercomparison Experiment (FirnMICE),
J. Glaciol.,
63, 401–422, https://doi.org/10.1017/jog.2016.114, 2017. a, b, c
Mellor, G. L. and Blumberg, A. F.:
Modeling Vertical and Horizontal Diffusivities with the Sigma Coordinate System,
Mon. Weather Rev.,
113, 1379–1383, https://doi.org/10.1175/1520-0493(1985)113<1379:mvahdw>2.0.co;2, 1985. a
Meyer, C. R., Keegan, K. M., Baker, I., and Hawley, R. L.: A model for French-press experiments of dry snow compaction, The Cryosphere, 14, 1449–1458, https://doi.org/10.5194/tc-14-1449-2020, 2020. a
Morland, L., Kelly, R., and Morris, E.:
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
Morland, L. W.:
A fixed domain method for diffusion with a moving boundary,
J. Eng. Math.,
16, 259–269, https://doi.org/10.1007/bf00042720, 1982. a
Paterson, W.: The Physics of Glaciers, third edition, Elsevier, Oxford, England,
https://doi.org/10.1016/C2009-0-14802-X, 1994. a
Schleef, S., Löwe, H., and Schneebeli, M.: Influence of stress, temperature and crystal morphology on isothermal densification and specific surface area decrease of new snow, The Cryosphere, 8, 1825–1838, https://doi.org/10.5194/tc-8-1825-2014, 2014. a
Schweizer, J., Jamieson, J. B., and Schneebeli, M.:
Snow avalanche formation,
Rev. Geophys.,
41, 1016, https://doi.org/10.1029/2002rg000123, 2003. a
Simson, A. and Kowalski, J.: Eulerian_Lagrangian_snow_solver: final paper submission TC, Zenodo [code], https://doi.org/10.5281/zenodo.5588308, 2021. a
Viallon-Galinier, L., Hagenmuller, P., and Lafaysse, M.:
Forcing and evaluating detailed snow cover models with stratigraphy observations,
Cold Reg. Sci. Technol.,
180, 103163, https://doi.org/10.1016/j.coldregions.2020.103163, 2020. 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
Wever, N., Valero, C. V., and Fierz, C.:
Assessing wet snow avalanche activity using detailed physics based snowpack simulations,
Geophys. Res. Lett.,
43, 5732–5740, https://doi.org/10.1002/2016GL068428, 2016. a
Wiese, M. and Schneebeli, M.:
Early-stage interaction between settlement and temperature-gradient metamorphism,
J. Glaciol.,
63, 652–662, https://doi.org/10.1017/jog.2017.31, 2017. a
Wingham, D. J.:
Small fluctuations in the density and thickness of a dry firn column,
J. Glaciol.,
46, 399–411, https://doi.org/10.3189/172756500781833089, 2000. a, b
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.
This companion paper deals with numerical particularities of partial differential equations...