A novel numerical implementation for the surface energy budget of melting snowpacks and glaciers
-
Published:2024-03-01
Issue:4
Volume:17
Page:1903-1929
-
ISSN:1991-9603
-
Container-title:Geoscientific Model Development
-
language:en
-
Short-container-title:Geosci. Model Dev.
Author:
Fourteau KévinORCID, Brondex JulienORCID, Brun FannyORCID, Dumont MarieORCID
Abstract
Abstract. The surface energy budget drives the melt of the snow cover and glacier ice and its computation is thus of crucial importance in numerical models. This surface energy budget is the result of various surface energy fluxes, which depend on the input meteorological variables and surface temperature; of heat conduction towards the interior of the snow/ice; and potentially of surface melting if the melt temperature is reached. The surface temperature and melt rate of a snowpack or ice are thus driven by coupled processes. In addition, these energy fluxes are non-linear with respect to the surface temperature, making their numerical treatment challenging. To handle this complexity, some of the current numerical models tend to rely on a sequential treatment of the involved physical processes, in which surface fluxes, heat conduction, and melting are treated with some degree of decoupling. Similarly, some models do not explicitly define a surface temperature and rather use the temperature of the internal point closest to the surface instead. While these kinds of approaches simplify the implementation and increase the modularity of models, they can also introduce several problems, such as instabilities and mesh sensitivity. Here, we present a numerical methodology to treat the surface and internal energy budgets of snowpacks and glaciers in a tightly coupled manner, including potential surface melting when the melt temperature is reached. Specific care is provided to ensure that the proposed numerical scheme is as fast and robust as classical numerical treatment of the surface energy budget. Comparisons based on simple test cases show that the proposed methodology yields smaller errors for almost all time steps and mesh sizes considered and does not suffer from numerical instabilities, contrary to some classical treatments.
Funder
H2020 European Research Council
Publisher
Copernicus GmbH
Reference50 articles.
1. Albert, M. R.: Computer models for two-dimensional transient heat conduction, https://apps.dtic.mil/sti/pdfs/ADA134893.pdf (last accessed: 30 November 2023), 1983. a 2. Anderson, E. A.: A point energy and mass balance model of a snow cover, https://repository.library.noaa.gov/view/noaa/6392 (last accessed: 30 November 2023), 1976. a, b, c 3. Barnett, T. P., Adam, J. C., and Lettenmaier, D. P.: Potential impacts of a warming climate on water availability in snow-dominated regions, Nature, 438, 303–309, https://doi.org/10.1038/nature04141, 2005. a 4. Barrett, A. I., Wellmann, C., Seifert, A., Hoose, C., Vogel, B., and Kunz, M.: One Step at a Time: How Model Time Step Significantly Affects Convection-Permitting Simulations, J. Adv. Model. Earth Sy., 11, 641–658, https://doi.org/10.1029/2018MS001418, 2019. a 5. Bartelt, P. and Lehning, M.: A physical SNOWPACK model for the Swiss avalanche warning: Part I: numerical model, Cold Reg. Sci. Tech., 35, 123–145, https://doi.org/10.1016/S0165-232X(02)00074-5, 2002. a, b, c, d, e, f, g, h, i, j
|
|