Evolution and drivers of the large-scale surface freeze-back onset in Siberian permafrost regions with the ERA5-Land reanalysis
Abstract. Permafrost is a defining feature of the Arctic and sub-Arctic environments. The extent of the permafrost region accounts for about a quarter of the Northern Hemisphere’s terrestrial surface. While most research on permafrost has focused on the summer season—when the active layer is thawed and carbon emissions peak—the late shoulder season, marking the transition season between summer and winter and ending by surface freeze-back at the large scale, has received less attention. Yet, about 14 % of the annual mean methane emissions from the permafrost occur during the refreezing period of the active layer. Understanding the seasonality, interannual variability and long-term trends of the surface freeze-back is therefore crucial to better constrain the high-latitude atmospheric carbon budget and improve Earth System Model projections. In this study, we analyze the evolution of surface freeze-back onset from 1950 to 2020 using the ERA5-Land reanalysis (0.1° spatial resolution) over a large region of Siberia encompassing the four main permafrost types. We find that surface freeze-back onset has been delayed by five days on average over that 70-yr period. Through spatial regression modeling, we show that, while several climatic and geographic factors influence freeze-back timing, the driving factor at the large scale (~kilometers) is the date when the 2-meter air temperature first falls below 0 °C, followed by the snow cover depth. These findings complement previous research that focused on the small scale (~meters), which emphasized the importance of the vegetation type and the snow cover characteristics at these spatial scales. Our results provide new insights into changes during the late shoulder season in one of the world’s fastest-warming regions and identify key variables to monitor for improving sub-seasonal forecasts that could become relevant to infrastructure upgrade and logistics planning in permafrost-affected areas.
Authors presented evolution and drivers of the surface freeze-back onset in Siberian permafrost regions based on the ERA5-Land reanalysis. The results show a delay in the surface freeze-back onset of approximately 5 days since 1950. The most significant controlling factors were identified as near-surface air temperature and snow depth. Furthermore, the positive trend in the freeze-back onset was found to be weakest in regions with continuous permafrost.
However, in its current form, the manuscript does not appear to contribute substantial new insights to the scientific understanding of how environmental changes interact with soil freeze-back processes. This limitation stems partly from the somewhat subjective and potentially problematic approach taken in the selection and testing of candidate predictors. The paper could be significantly strengthened through more in-depth spatial analysis. General comments and specific suggestions are provided below for the authors' consideration.
General comments
1) Candidate predictors
The authors dedicate a substantial portion of the manuscript to discussing the importance of the selected candidate predictors. However, the results indicate that the influence of most predictors (3 out of 5) is not statistically significant. While the coarse scale and inherent biases of the ERA5-Land dataset may explain why such influences are not captured, I strongly recommend that the authors condense this section, particularly the discussions concerning vegetation type, groundwater content, and latitude. Additionally, the authors mention the potential importance of snowfall timing but do not provide further analysis on this point.
2) Manuscript Structure
Abstract: Please consider making the abstract more concise by removing information that is either well-known (e.g., the spatial resolution of ERA5-Land and general statements on permafrost extent) or not directly relevant to the study's core findings (e.g., methane emissions). The abstract should function as a summary and contain only the most critical points of the article.
Introduction: A significant part of the introduction (2 out of 5 paragraphs) focuses on the carbon cycle feedback. While the permafrost-carbon cycle is undoubtedly important, the introduction would benefit from a better balance. Please consider reformulating this section to also adequately highlight other crucial processes significantly affected by the freeze-back onset, such as hydrological cycles and phenology.
Furthermore, the extensive body of research on surface freeze/thaw status derived from remote sensing (e.g., Kim et al., 2012, 2017) is overlooked and should be incorporated into the literature review.
Results: The readability of this section is compromised because it frequently incorporates methodological explanations (e.g., L233-235, 251-255, 273-274). These details should be relocated to the Methods section to allow the Results to focus solely on the findings and their implications.
3) Nomenclature
The nomenclature for some variables is unconventional. For instance, in Equation 4, "Qno snow" could be changed to the more standard notation "Qsf" (for snow-free conditions) or "Qsc" (for snow-covered conditions). Additionally, the variable "dayair temp" in Equation 6 is not commonly used and should be revised for clarity.
4) Methology
Daily mean ground temperature: is defined as the arithmetic average of the temperature at 00:00 and 12:00 UTC. Given that ERA5-Land provides hourly soil temperature data, could the authors clarify why only two specific timesteps were selected for this calculation? It would be more representative to use the full diurnal cycle
Snow Cover: While snow insulation plays an important role in the soil thermal regime at seasonal and annual scales, its relevance to the freeze-back onset is less clear, as this transition typically occurs before a snow depth exceeding 14 cm is established. In my opinion, the freeze-back onset may be more closely linked to early snowfall events, which can cool the ground via the latent heat during snowmelt.
The use of a 14 cm snow depth as an insulating threshold (based on Eq. 5) depends on several factors, including snow thermal conductivity, the chosen insulation threshold (e.g., 5%), and the values of parameters k and u. Although the use of typical values is a reasonable starting point, additional sensitivity experiments are needed to assess the influence of these parameters and improve the robustness of the algorithm. Furthermore, please provide a reference for the snow thermal conductivity value of 0.25 W m-1 K-1—this value appears relatively low for a snow density of 350 kg m-3.
5) Spatial analysis
This study could be significantly enhanced by additional spatial analysis. For instance, mapping the onset trend with significant patterns across the entire study area would be particularly informative.
Specific Comments
L28: Cohen et al., 2020? This is not a permafrost mapping research
L84: Consider using the updated referece of Obu et al., 2019.
L199: Typo here?
L224: 1950-200?
References:
Kim, Y., Kimball, J. S., Zhang, K., and McDonald, K. C.: Satellite detection of increasing Northern Hemisphere non-frozen seasons from 1979 to 2008: Implications for regional vegetation growth, Remote Sensing of Environment, 121, 472–487, 2012.
Kim, Y., Kimball, J. S., Glassy, J., and Du, J.: An extended global Earth system data record on daily landscape freeze–thaw status determined from satellite passive microwave remote sensing, 2017.
Obu, J., Westermann, S., Bartsch, A., Berdnikov, N., Christiansen, H. H., Dashtseren, A., Delaloye, R., Elberling, B., Etzelmüller, B., Kholodov, A., Khomutov, A., Kääb, A., Leibman, M. O., Lewkowicz, A. G., Panda, S. K., Romanovsky, V., Way, R. G., Westergaard-Nielsen, A., Wu, T., Yamkhin, J., and Zou, D.: Northern Hemisphere permafrost map based on TTOP modelling for 2000–2016 at 1 km2 scale, Earth-Science Reviews, 193, 299–316, 2019.