Snowpack depth hoar impacts on the permafrost thermal regime in a land surface model
Abstract. In cold, dry regions, snowpacks can build a basal depth hoar layer that has a lower density that effectively increases the insulative ability and thus warms permafrost. However, the precise magnitude of the depth hoar insulation and the subsequent impact on permafrost thermal regimes remains highly uncertain. Current Earth and land surface models do not consider the depth hoar effect in their simulations. This study coupled a depth hoar scheme into the snowpack processes within a land surface model. The depth hoar effect on the permafrost thermal regimes in present (1980–2016) and future (2015–2100) climates was assessed by model simulations excluding and including depth hoar. Snow and soil thermal properties observed at the field sites were used to evaluate the model results. The modeling assessment revealed that the depth hoar produced permafrost warming of 2.5–4.0 °C in February under the present climate and deepened the active layer thickness by up to 20 cm. The insulative ability of depth hoar was further significant for colder air temperatures and a thicker snowpack. The effect was demonstrated by increased temperature differences between the air and soil depth at 20 cm with increasing snow depth under colder air temperature conditions. Under a future warming climate, the depth hoar functioned as an insulator enhancing permafrost degradation. Our analysis reveals a new perspective on the underlying influence of depth hoar as a driver of permafrost degradation in present and future climates. We recommend improvements to model representation to include the depth hoar process to improve simulations of the permafrost-related changes under warming climate.
This manuscript presents a valuable and highly relevant investigation of how the representation of depth hoar affects simulated snow properties and the thermal regime of permafrost under present and future climatic conditions. The study is scientifically sound, the experimental design and methods are appropriate, and the analyses clearly support the conclusions drawn by the authors. In particular, the evaluation against observations demonstrates that the implemented depth-hoar formulation improves not only the simulation of snow depth and snowpack properties, but also the resulting soil temperatures. The manuscript therefore makes an important contribution to land surface modelling by addressing a consequential shortcoming in current representations of Arctic snow and by demonstrating the relevance of snow stratigraphy for simulating snow–permafrost interactions. I consider the study fully worthy of publication once the presentation has been substantially improved.
My principal concern is the quality and clarity of the English throughout the manuscript. In many places, I was able to infer the intended meaning because I am familiar with the subject matter. However, the wording itself is frequently unclear, imprecise, or potentially misleading, and readers who are less familiar with the topic may interpret some statements differently from what the authors appear to intend. This is therefore more than a matter of style or occasional grammatical errors: the language sometimes obscures the scientific reasoning and, in some cases, changes the apparent meaning of a statement. The manuscript requires thorough language revision, ideally by a fluent English speaker with sufficient familiarity with the subject to preserve the intended scientific meaning. My detailed comments below identify a number of examples, but they should not be regarded as an exhaustive language correction of the manuscript.
Please also find a number of detailed comments unrelated to the language issue below:
Abstract
Page 1, Line 13: I would suggest to mention which model here.
Page 1, Line 18: I find this sentence unclear. Do you mean that depth hoar had more impact in zour simulations under colder air temperatures and thicker snow packs than under warmer temperatures and more shallow snow packs? If so, please reformulate to make the sentence clearer.
Page 1, line 19: “soil at 20cm depth” instead of “soil depth at 20 cm"
Introduction
Page 2, line 35: “increase of snow depth in eastern” instead of “increase of snow depth in the eastern”
Page 2, line 47 “snow accumulation and usually found at the” instead of “snow accumulation and is usually found at the”
Page 2, line 48-49: This sentence is unclear to me. Do you mean to say that depth hoar increases the role of the snow pack as thermal insulator? Should be rephrased.
Page 2, line 53: Maybe consider adding Domine et al., 2019, who showed the same effects as Dutch at al with two other models, https://doi.org/10.1029/2018MS001445
Page 2, line 53: The term anomaly is ambiguous here. Maybe use biases, or even errors.
Page 2, line 60-63: I would suggest to add that the analysis in Dutch et al. is for one observational site in Canada.
Page 3, line 66: “depth hoar schemes within land surface models” instead of “depth hoar scheme within land surface models”
Page 3, line 68-69: Damseaux et al did not study how depth hoar forms, they only discussed physical depth hoar properties taken from other studies. This citation should be removed here.
Page 3, line 90-91: I would suggest to acknowledge here that most CMIP6 models don't represent permafrost carbon processes well (eg Schädel et al., 2024, https://doi.org/10.1038/s41558-023-01909-9), so projections of the permafrost carbon feedback have a high uncertainty (IPCC), which is intensified by biases in snow
2.2 Snow model
2.2.1 General description
Page 4, line 126-127: For a better flow of your argument, consider changing this sentence into something like "Since the density and liquid water content of the snowpack can in reality have vertically heterogeneous ..."
2.3 Study site and observation
Page 9, line 244: How is ALT monitored at these sites?
2.4 Model simulation
Page 9, line 249: Is there a reason why ERA-Interim was used instead of the more recent ERA5 reanalysis?
Page 9, line 270: From this sentence, it remains unclear whether you used the IPSL forcing downscaled by ISIMIP or forcing from another model.
3 Result
3.1.1 Snow depth and snow water equivalent
Page 10, line 285: Please make the mechanisms resulting in different snow depths and swe explicit here. I assume the differences between the snow depths in the two simulations arise from the difference in snow density, causing similar amounts of SWE to result in greater snow depths. If there is identical forcing, what causes the differences in swe during the build up of the snow pack?
3.1.3 Active layer thickness
Page 12, line 315: How exactly is active layer thickness derived from the model?
3.2 Seasonal variation of snow grain form
Page 13, line 335 (figure 5): Is this also an 11-year average? How is snow type averaged?
3.4 Correlation between air-soil temperature difference and snow depth
Page 16, line 402 (figure 9): In figure 8, and also in the text, you state that depth hoar has the strongest impact on insulation in the regime with the coldest air temperatures. Figure 9 seems to contradict that, the shifts in the pdfs of the moderate and warm air temperature regime are much more pronounced than those in the cold air temperature regime. Why is that?
Discussion
4.1 Impact of depth hoar on snow properties
Page 17, lines 425–426: The model was configured so that depth hoar has a greater insulating capacity; the simulation therefore does not independently confirm this. In addition, no observations from this study are presented in which depth hoar itself was measured. Please reformulate this statement.
Page 17, line 436: I find this formulation misleading. Parameterizations of snow thermal conductivity based on snow density are vulnerable to errors in simulated snow density, but this applies to any scheme that uses snow density as its basis. As you point out later, using different methods to estimate snow thermal conductivity can improve the results.
4.2 Impact of depth hoar on soil temperature
Page 18, lines 446–447: The differences between the two simulations show the impact of implementing depth hoar; they do not demonstrate that including depth hoar improved the simulation.
Page 19, lines 465–466: Does this statement contradict Figure 9, in which the probability density functions for the depth-hoar and no-depth-hoar simulations appear most similar at the coldest air temperatures?
Page 19, lines 486–488: In addition to reduced depth-hoar formation in a warmer climate, I would expect the shortening of the snow season to play a role. Could the reduced influence of the depth-hoar formulation during the shoulder seasons in the future projections also result from a decrease in snow cover itself?
4.3 Uncertainty and future work
Page 21, line 516: Did you mean ‘colder and warmer’?
Page 21, line 518: In the description of the model simulation, you write that ‘the soil organic carbon profile was simulated by CHANGE at each time step’. Does this mean that CHANGE calculates soil organic carbon, but that it does not affect soil heat capacity or soil thermal conductivity?