the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 12 May 2025
Simulating snow properties and Ku-band backscatter across the forest-tundra ecotone
Abstract. Sophisticated snowpack models are required to provide accurate estimation of snowpack properties across the forest-tundra ecotone where in situ measurements are sparse. As snowpack properties strongly influence radar scattering signals, accurate simulation is crucial for the success of spaceborne radar missions to retrieve snow water equivalent (SWE). In this study, we evaluate the ability of default and Arctic versions of Crocus embedded within the Soil, Vegetation and Snow version 2 (SVS2-Crocus) land surface model to simulate snowpack properties (e.g. depth, density, SWE, specific surface area; SSA) across a 40-km transect of the Northwest Territories, Canada, using two winter seasons (2021–22 & 2022–23) of in situ measurements. An ensemble of simulated snowpack properties (120 members from default and Arctic SVS2-Crocus) was used in the Snow Microwave Radiative Transfer (SMRT) model to simulate Ku-band (13.5 GHz) backscatter. Modelled backscatter using multi-layer SVS2-Crocus snowpack simulations were compared to backscatter using a simplified 3-layer radar-equivalent snowpack. Results highlight that Arctic SVS2-Crocus wind-induced compaction modifications were spatially transferable across the forest-tundra ecotone, reducing the RMSE of surface density by an average of 29 %. Basal vegetation modifications were less effective in simulating low-density basal snow layers at all sites (2022 & 2023; default RMSE: 67 kg m-3; Arctic RMSE: 69 kg m-3) but were necessary to simulate a physically representative Arctic density profile. SVS2-Crocus underestimated SSA leading to high errors in the simulation of snow backscatter (2022 & 2023; default RMSE 3.5 dB; Arctic RMSE: 4.8 dB). RMSE of backscatter was reduced by implementing a minimum SSA value (8.7 m2 kg-1; 2022 & 2023; default RMSE: 1.5 dB; Arctic RMSE: 1.5 dB). A radar-equivalent snowpack was effective in retaining the scattering behaviour of the multi-layer snowpack (RMSE < 1 dB) providing a means to estimate SWE with increased computational efficiency.
Competing interests: Some authors are members of the editorial board of The Cryosphere.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.- Preprint
(2779 KB) - Metadata XML
- BibTeX
- EndNote
Status: final response (author comments only)
- RC1: 'Comment on egusphere-2025-1498', Anonymous Referee #1, 13 Jun 2025
-
RC2: 'Comment on egusphere-2025-1498', Anonymous Referee #2, 22 Aug 2025
This manuscript presents a comprehensive assessment of Arctic-modified SVS2-Crocus snow model performance across the forest-tundra ecotone and its implications for microwave remote sensing applications. Although the work addresses important questions regarding snow model transferability and microwave retrieval preparation, several methodological and interpretational concerns limit its impact.
Major comments:
- The meteorological forcing data are from HRDPS, which is distributed in space with 25 km spatial resolution. You simulated snow cover based on a point scale, then you validated the simulated results using ground observations. I am confused about how the 2.5 km HRDPS data is applied to the model (point extraction vs. spatial interpolation) and whether simulations are run as single points or on a spatial grid. Did you downscale the forcing data before inputting them into the snow model? The spatial resolution of 25 km is too coarse, and it’s hard to validate simulated snow cover from HRDPS using point-based ground observations, especially in the complex terrain. That’s also one of the reasons for the uncertainties of simulation associated with wind speed within the HRDPS. Therefore, it could be more reasonable and decrease uncertainties to downscale the forcing data first if you did not do that.
- It is not clear how the 120 ensemble members are statistically processed. When reporting metrics like "default RMSE and Arctic RMSE" (Lines 288-291), it is unclear whether these represent ensemble mean performance, median values, or some other aggregation method.
- The MS repeatedly claims that improved backscatter simulation will advance SWE retrieval capabilities, but this logic is confused. If SVS2-Crocus provides snow density and depth, SWE calculation is trivial and does not require backscatter simulation. The authors have not explained clearly why simulating backscatter (forward modeling) helps retrieve SWE from measured backscatter (inverse problem). The MS discovers substantial backscatter simulation errors (Lines 364-372) that would severely compromise any retrieval algorithm. The authors resort to ad hoc corrections (minimum SSA values, scaling factors) that lack physical justification and would not be transferable to operational scenarios and more study regions. In other words, the MS combines snow model evaluation with microwave retrieval algorithm development without clearly articulating which problem it aims to solve. If the goal is snow model improvement, the microwave component adds unnecessary complexity. If the final goal is advancing SWE retrieval, the methodology does not address the fundamental challenges of operational retrieval algorithms. Therefore, please clarify the study aims or objectives specifically.
Minor comments:
Lines 107-108, regarding the two forest sites, can you give a more detailed description of their location? They are under the canopy or canopy gaps?
Lines 157-158, how and when (summer or winter) did you measure the polar vegetation heights? These heights (0.1-0.35 m) seem static, but shrub bending under snow load is dynamic. Did you consider that?
Lines 177-178, you mentioned the range of polar vegetation height (0.1-0.35 m) before. When you used polar vegetation height in Arctic SVS2-Crocus parameterization, did you use a fixed value or changing values? It’s so simple to create a binary threshold at all sites (tundra, shrub, and forest): below this height = vegetation effects active, above = normal snow physics, especially several vegetation types in your study region. For example, why would 0.35 m shrub effects apply in 10 m tall forests?
Line 238, I know “TVC” represents Trail Valley Creek, but this is the first time you've used the abbreviation. You should also indicate its full name; similar problems also exist for other sites. In addition, I'm confused about how you named the seven sites. Either all of them are named after places, or all of them are named after vegetation types.
Figs. 5 and 6, It could be better to show different pit measurements in the legend.
Lines 459-461, "some simulated profiles can be shallower than measured profiles as a function of the precipitation inputs meaning some polar vegetation heights encompass much of the simulated profile", is the snow thermal conductivity changes influenced by shrub considered during the parameterization processes? Except for the decreased wind-induced snow compaction, the changes in snow thermal conductivity are also important to snow energy and mass balance as well as soil thermal regime.
The Hedstrom and Pomeroy (1998) interception model requires some vegetation information, such as LAI and canopy coverage. Where did you get them?
Citation: https://doi.org/10.5194/egusphere-2025-1498-RC2
Data sets
Meteorological forcing data, SVS2-Crocus and SMRT simulated output and measurements of snowpack properties (2022-2023) Georgina Woolley https://doi.org/10.5281/zenodo.15091024
Model code and software
SVS2-Crocus Model Code Georgina Woolley et al. https://doi.org/10.5281/zenodo.15091095
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 739 | 111 | 20 | 870 | 17 | 34 |
- HTML: 739
- PDF: 111
- XML: 20
- Total: 870
- BibTeX: 17
- EndNote: 34
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
General
The study compares simulated snow parameters using SVS2-Crocus to in situ measurements at a study area in the NWT, Canada. Two versions of the model are applied; the default model and an Arctic-specific modification. Several locations with differing vegetation and snow conditions are analyzed. Furthermore, a forward model is used to simulate microwave backscatter from SVS2-Crocus outputs, comparing these to backscatter simulations using the in situ data directly. A microwave-effective snowpack concept, which aggregates the data to three representative layers, is used. The paper is of interest to the scientific community as it has become clear that some type of fusion of modelled and remote sensing information is required for successful retrieval of SWE from spaceborne radar. This concerns, in particular, approaches using Ku-band SAR backscatter. This paper takes some steps to attempt to quantify how an advanced coupled land-surface -snow process model reproduces natural snowpacks in a challenging environment, and what are the implications on (simulated) backscatter. As such, the study is worthy to consider for publication.
The study is also generally well written although some specifics in the methodology are hard to grasp and require several readings. For example, from the abstract it was not at first obvious that actual observations of backscatter are not used, only simulations. Furthermore, it is hard to discern where exactly the three layer aggregates (radar equivalent snowpack) where used: in simulations based on SVS2-Crocus, in simulations from snowpit data, or perhaps both? Terminology referring to the sites also changes occasionally, sometimes referring to the biome (tundra, forest), sometimes to the names defined in Figure 1. These are the main examples, but the complexity of the diverse model settings makes the paper somewhat hard to follow and requires particular care in describing the experimental setups.
My main concern, however, is related to the usefulness of the backscatter simulation setup itself. The SMRT model is treated as a black box, testing which kind of numbers come out with each version of the input data. It seems that the pit data forcing is used as the “truth”, with SVS2-Crocus -based simulations representing deviations (“errors”) from this. No real effort is placed on which parameters actually induce these differences, beyond testing different approaches for tuning the optical grain size/SSA. Since the study is based on only simulated backscatter, one could expect a thorough sensitivity study on different parameters, or something similar. Perhaps using the 120 ensemble members in SVS2-Crocus makes this a challenge; however, you could then consider dropping the ensemble approach, and use individual, controlled simulations?
Further, some confusion is created by first choosing to fix a scaling parameter in the forward model (the polydispersity factor K) at unity, only to re-introduce another scaling parameter with basically the same end result, as adjusting K would have had. I realized only after some time that the idea is to try to scale the SVS2-Crocus SSA to observations; however, the choice of the scaling factors tested seems arbitrary. Of some interest could be to try to derive your own optimal scaling for 1) SSA 2) density 3) possibly another parameter, such as snow depth, and with the (average) optical scaling reconduct the simulation exercise.
Please see major comments below for specifics.
Major comments
Throughout the paper, please choose which name you use when referring to sites. Now, sometimes “forested sites” and “tundra” are used, sometimes “Havikpak” etc. You could also use always both to be explicit e.g. “Upper plateau (tundra)”
I also fail to see why the scaling factor of 0.63 derived in a paper from 2011 could be relevant. I would suggest another approach: 1) derive your optimal scaling factors required to match SVS2-Crocus to field data 2) do this at least for SSA, density and snow depth 3) analyze which is the most important factor to scale (which has the largest impact) for each surface type. As an optimal case, you could test scaling all three+ parameters to match the field data, and see how much RMSD remains from variability in the ensembles.
Again, it is not clear to which data the n-layer -> 3-layer conversion is applied. Apparently to SVS2-Crocus at least, but is it applied to field data as well? Or are filed data always simulated “as is”?
Whole section 4.0 on results. It is quite tedious to read endless RMSE values in the text, and the point is quickly lost. This really gets out of hand e.g. on p. 16. Please tabulate the results, refer in the text to the tables, highlighting only the most important results.
It is also not clear against which SD data the RMSEs are calculated against. Magnaprobe, pit, or both?
Same question as above for SWE. are the RMSEs against Pit or SWE Tube values?
For density, can you not calculate a reference value from SWE tubes? Do you have the snow depth from the SWE core site recorded?
Minor comments