We thank the reviewer for the careful reading of the manuscript, the appreciations towards our work and the very constructive comments. Below we provide a point-by-point response.

General comment

For the evaluation of the land surface datasets against IT-SNOW (for SWE and snow cover duration) or SCIA (for 2-m temperature and precipitation), the authors look into spatial biases between the different products, but do not evaluate the statistical significance of these biases, which make me wonder whether the biases presented are statistically significant or not. As the authors work with a relatively short period (2010-2024) and a large number of grid points, I would suggest assessing the statistical significance of SWE, snow cover duration, temperature and precipitation biases using a two-tailed students T-test combined with a False Discovery Rate (FDR) approach following Wilks (2016).

 Reply: We would like to thank the Reviewer for this highly pertinent and constructive observation, which significantly enhances the statistical robustness of our findings. We fully agree that evaluating the statistical significance of the spatial biases is crucial, particularly given the relatively short evaluation period (2010–2024) and the large number of grid points involved in the analysis. Following your valuable suggestion, we have already carried out the analysis assessing the statistical significance of the biases for all the variables using the proposed methodology. In the revised version of the manuscript, we will update the results section accordingly. Specifically, rather than overlaying stippling directly onto the main manuscript maps—which would severely compromise their readability given the spatial resolution (9 km), the small pixel size, and the compact format of the figures—we plan to provide dedicated additional materials. These will likely be included in a new Supplement and will probably consist of (i) separate maps explicitly showing the spatial distribution of the statistical significance, and (ii) additional graphical outputs illustrating the relationship between the magnitude of the bias and the significance of the estimate. However, the results of this new analysis confirm our initial findings. The statistical significance is generally consistent with the interpretations provided in the original draft; therefore, the core message and the overall conclusions of the Technical Note will remain unchanged, but they will now be supported by a much more rigorous statistical framework. We sincerely thank the Reviewer again for pointing this out and for helping us improve the quality of our manuscript.

Specific comments (SC)

SC1: L124 and Figure 1: Could the authors add more geographical information to the map of Italy (other than “Alps” and “Apennines”) to help readers locate places more easily? For example, could they indicate the location of the Ligurian sector, Sardinia, Sicily, the Strait of Messina, etc.? I know where these places are, but I imagine that not all readers do, so it would be helpful to include a little more geographical background information.

Reply to SC1: We agree that adding more geographical references to Figure 1 will help the international readership better orient themselves and follow the regional differences discussed in the text. In the revised version of the manuscript, we will update Figure 1 to include additional geographical background information.

SC2: L146-147: Could the authors be a little bit more specific on how the precipitation-phase partitioning is implemented in IT-SNOW?

Reply to SC2: We thank the Reviewer for pointing this out. Phase partitioning in S3M relies on both air temperature and relative humidity, which generally provides more reliable results than using air temperature alone (Zhang et al., 2017). Rain and snow proportions are derived based on the simple method by Froidurot et al. (2014). The separation is performed at an hourly time scale. In the revised manuscript, we will include these additional details.

SC3: L148: Is the hybrid temperature-index and radiation-driven melt approach in IT-SNOW similar to the enhanced temperature-index approach of Pellicciotti et al., 2005 (https://doi.org/10.3189/172756505781829124)?

Reply to SC3: Yes, the approach in S3M was inspired by the approach of Pellicciotti et al. (2005). The only, significant deviation is the addition of two modulation parameters that reduce snowmelt based on antecedent temperature. This was done to include in the approach some forms of cold content, that is, thermal inertia of snow. While a comprehensive description is available in Avanzi et al. (2023), we agree that briefly clarifying this connection will benefit the readers. Accordingly, we will include these details and the reference to Pellicciotti et al. (2005) in the revised text.

SC4: Appendix Figures B1-B3: For the middle row figures (representing the temperature biases) I recommend reversing the color bar. In the current state, the colors representing the Tbias are somewhat counterintuitive in my opinion, i.e. warmer is purple and colder is red whereas a positive SWE bias is blue (usually associated with colder temperature conditions) and a negative SWE bias is red (usually associated with warmer conditions).

Reply to SC4: We completely agree that the current color combination for the temperature bias might be counterintuitive for the readers. In the revised manuscript, we will update the color scale. Specifically, we will reverse the color bar, using red to indicate temperature overestimation (positive bias).

SC5: L358-359: Another potential explanation for the overestimation of snow depth in the higher parts of the Alps could be that (according to Orsolini et al., 2019; https://doi.org/10.5194/tc-13-2221-2019) data assimilation of snow cover (using IMS – Interactive Multisensor Snow and Ice Mapping System) was discontinued in ERA5 above 1500 m a.s.l. This could explain the snow biases in ERA5 itself, but as snow biases also affect ERA5 temperature (via snow-albedo feedbacks) I can imagine that the negative temperature biases propagated from ERA5 into ERA5-Land also have affected the snow conditions in ERA5-Land, particularly at high Alpine elevations. And what about snow-albedo feedback in ERA5-Land itself? Could these feedbacks also have contributed to the SWE bias in ERA5-Land?

Reply to SC5: We thank the Reviewer for this highly insightful comment and for pointing out the relevance of Orsolini et al. (2019). We agree that these represent very plausible hypotheses that could help explain the observed biases. First, it is reasonable to hypothesize that the discontinue IMS snow cover assimilation above 1500 m a.s.l. in ERA5 introduces snow and temperature biases in the atmospheric reanalysis. These biases may then propagate into ERA5-Land through the atmospheric forcing. Specifically, the temperature at the “lowest model level” (10 m), which drives the offline land surface model after being adjusted for elevation differences between the two grids via a standard lapse rate. Second, regarding the snow-albedo feedback, we agree that this could also play a role. Specifically, while ERA5-Land runs offline (meaning there is no two-way land-atmosphere coupling), such a mechanism could be driven by the internal thermodynamics of its underlying land surface model CHTESSEL. An initial overestimation of SWE could lead to an anomalously high surface albedo, which would reduce the absorbed shortwave radiation. This process might lower the computed skin temperature and reduce the sensible heat flux, thereby potentially delaying snowmelt. Consequently, this internal thermodynamic loop may act to exacerbate and prolong the positive SWE biases, particularly at higher Alpine elevations. In the revised manuscript, we will expand the discussion (around L358-359) to explicitly cite Orsolini et al. (2019) and to discuss these two intertwined, potential mechanisms—the propagation of biases from the ERA5 forcing and the possible internal thermodynamic loop within CHTESSEL—as plausible drivers of the overestimation in the Alpine sector.

SC6: L371-374 and beyond: I guess the negative SWE biases in the Apennines or some regions of the Alps, like the Trentino-Alto Adige region  could also be related to differences in how ERA5-Land (via precipitation forcing from ERA5) or CERRA-Land/VR-REA-IT and the S3M model (the basis for IT-SNOW) deal with the precipitation phase, particularly for temperature regimes around 0°C. I don’t know how the precipitation-phase partitioning in S3M works, but could it be that precipitation falls as snow in S3M, whereas in ERA5-Land or the other 2 models precipitation falls as rain? To address this question, I recommend adding two more rows to Figures B1-B3 that visualize the seasonal biases in rainfall and snowfall.

Reply to SC6: We fully agree that differences in the precipitation-phase partitioning (solid vs. liquid) around 0°C between the evaluated products and the S3M model are highly likely to be a key driver of the SWE biases observed in such regions. We agree with the rationale behind your suggestion to visualize the seasonal biases in rainfall and snowfall. However, we are unfortunately unable to perform this specific evaluation due to data availability constraints. To ensure a consistent and "symmetrical" evaluation, our analysis relies strictly on the variables that are publicly and openly distributed for all datasets. While ERA5-Land does provide partitioned snowfall and rainfall outputs, these specific fluxes are not available in the public data releases for CERRA-Land and VHR-REA_IT. Also, the reference dataset (IT-SNOW) includes only the final snow variables (SWE, snow depth, bulk snow density and liquid water yield), but not the snowfall and rainfall fields internally computed by the S3M model. Consequently, computing a direct, comparative spatial bias for snowfall against the reference dataset is not feasible. While a quantitative evaluation is not feasible, during the revision process we will evaluate how best to acknowledge this aspect in the text. We plan to briefly mention the possibility that differences in the precipitation-phase partitioning approaches could potentially influence the SWE discrepancies, particularly in regions where temperatures frequently fluctuate around the freezing point.

SC7: Also, I am wondering whether there are differences in the density of ultrasonic sensors (used for assimilation of IT-SNOW) among the different regions? Could it be that the density of sensors is higher in the central Apennines and some regions of the Alps, such as Trentino Alto-Adige? I am asking this because the negative SWE biases presented in Figure 3 seem to appear consistently over the same regions, independent from the comparison with the different land surface datasets. Even the SWE biases between ERA5 Land and IT-SNOW show patches of negative SWE biases in the Trentino Alto-Adige region, which seem to appear in the lower-elevated areas of the respective region.

Reply to SC7: Yes, there are indeed regional differences: the bulk of the snow data are in the northern regions, particularly within the Po river basin in the northwest. While a lower sensor density might indeed affect the performance of the model, previous comprehensive assessments of IT-SNOW against a spatially distributed Sentinel-1 snow depth dataset did not reveal any clear spatial bias patterns linked to the sensor network density (Avanzi et al., 2023). So we are inclined to attribute most of these patterns to the evaluated datasets, rather than to IT-SNOW. We appreciate the Reviewer for prompting this clarification. In the revised manuscript, we will briefly expand upon this aspect in the text.

References

Avanzi, F., Gabellani, S., Delogu, F., Silvestro, F., Pignone, F., Bruno, G., Pulvirenti, L., Squicciarino, G., Fiori, E., Rossi, L., Puca, S., Toniazzo, A., Giordano, P., Falzacappa, M., Ratto, S., Stevenin, H., Cardillo, A., Fioletti, M., Cazzuli, O., Cremonese, E., Morra di Cella, U., Ferraris, L., 2023. IT-SNOW: a snow reanalysis for Italy blending modeling, in situ data, and satellite observations (2010–2021). Earth Syst. Sci. Data 15, 639–660. https://doi.org/10.5194/essd-15-639-2023

Froidurot, S., Zin, I., Hingray, B., Gautheron, A., 2014. Sensitivity of Precipitation Phase over the Swiss Alps to Different Meteorological Variables. J. Hydrometeorol. 15, 685–696. https://doi.org/10.1175/JHM-D-13-073.1

Orsolini, Y., Wegmann, M., Dutra, E., Liu, B., Balsamo, G., Yang, K., De Rosnay, P., Zhu, C., Wang, W., Senan, R., Arduini, G., 2019. Evaluation of snow depth and snow cover over the Tibetan Plateau in global reanalyses using in situ and satellite remote sensing observations. The Cryosphere 13, 2221–2239. https://doi.org/10.5194/tc-13-2221-2019

Pellicciotti, F., Brock, B., Strasser, U., Burlando, P., Funk, M., Corripio, J., 2005. An enhanced temperature-index glacier melt model including the shortwave radiation balance: development and testing for Haut Glacier d’Arolla, Switzerland. J. Glaciol. 51, 573–587. https://doi.org/10.3189/172756505781829124

Zhang, Z., Glaser, S., Bales, R., Conklin, M., Rice, R., Marks, D., 2017. Insights into mountain precipitation and snowpack from a basin‐scale wireless‐sensor network. Water Resour. Res. 53, 6626–6641. https://doi.org/10.1002/2016WR018825

We thank Prof. Margulis for the comments and feedback on our manuscript. We provide our reply in the attached pdf file.