the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
On the performance of the snow model Crocus driven by in situ and reanalysis data at Villum Research Station in northeast Greenland
Abstract. Reliable and detailed snow data are limited in the Arctic. We aim at overcoming this issue by addressing two questions: (1) Can the reanalysis ERA5 replace limited in situ measurements in high latitudes to drive snow models? (2) Can the Alpine model Crocus simulate reliably Arctic snow depth and stratigraphy? We compare atmospheric in situ measurements and ERA5 reanalysis and evaluate simulated and measured snow depth, density and specific surface area (SSA) in northeast Greenland (October 2014–October 2018). To account for differences between Alpine and Arctic region, we introduce a new parametrisation for the density of new snow.
Our results show a good agreement between in situ and ERA5 atmospheric variables except for precipitation, wind speed and direction. ERA5’s resolution is too coarse to resolve the topography in the study area adequately, leading presumably to the detected biases. Nevertheless, measured snow depth agrees better with ERA5 forced simulations than forced with in situ measurements.
Crocus can simulate satisfactory the evolution of snow depth, but simulations of SSA and density profiles for both forcings are biased compared to field measurements. Adjusting the new snow density parametrisation leads to improvements in the simulated snow stratigraphy. In conclusion, ERA5 can be used instead of in situ measurements to force snow models but the use of Crocus in the Arctic is affected by limitations likely due to the missing vertical water vapour transport and snow redistribution during strong winds. These limitations strongly affect the accuracy of the vertical profiles of physical snow properties.
- Preprint
(2588 KB) - Metadata XML
- BibTeX
- EndNote
Status: closed
-
RC1: 'Comment on tc-2021-100', Anonymous Referee #1, 07 May 2021
This paper presents first a validation of ERA5 over the Villum Research Station (NE Greenland) in the aim afterwards to force the CROCUS snowmodel with ERA5 and to evaluate its results with observations made at this station. While the paper is well written and uncertainties well discussed, the main scientific message, relevance and originality of this paper is poor for a journal such TC. ERA5 has already been evaluated over Greenland, CROCUS has already been applied in Arctic/Greenland and the fresh snow density improvements proposed in this paper seems to not improve a lot the comparison.
Moreover, this paper shows notably that using observed precipitation as forcing of CROCUS overestimates the snow height and that, although improvements in the prescribed fresh snow density, CROCUS still underestimates snow density. This suggests for me that:
=> as the modelled snow density is biased, both modelled and observed snowpack should be compared in water equivalent as observed density is available. A large part of discrepancies shown in Fig6 could be due to problems in the modelled snow density and not to overestimation/underestimation of prescribed precipitation as forcing. This could change the discussion in this paper if it is well the case. In Fig6, the modelled snow height can at least be corrected with respect to the snow density bias with observation and in this case, we could likely see a systematic underestimation of modelled snowpack using ERA5 as forcing.
=> as mentioned in the discussion, it is also very likely that this area is highly impacted by the wind erosion (increasing snow density and eroding the snow pack) as well as by the wind transport of the snow (biasing the snowfall measurements and potentially bringing snow from elsewhere) decreasing a lot the relevance of this study knowing that the blowing snow processes are not resolved here. Therefore, it is very likely that a large part of discrepancies between CROCUS and observations could just be due to the fact that this process is missing here and this does not necesserly mean that CROCUS is not working well in Arctic. Comparison with RCM (RACMO, MAR) taking into account these processes could help to evaluate the importance of these processes in this area.
Some more minor remarks:
- line78: I assume that this data is not assimilated into ERA5. Where is the closest weather station for which its observation is assimilated into ERA5?
- line 124: As the altitude of the ERA5 pixel is too high with respect to the station altitude, ERA5 is too cold and its surface pressure is too low. Converting all the precipitation to snow below 1°C will overestimate a lot the amount of snowfall with respect to rainfall knowing that ERA5 is already too cold.
- line 169, section 2.3: According to Fausto et al. (2018, A Snow Density Dataset for Improving Surface Boundary Conditions in Greenland Ice Sheet Firn Modeling, https://doi.org/10.3389/feart.2018.00051) using a fixed fresh snow density of 315 km/m3 is better than any parametrisation based on temperature and wind. In MARv3.11 model, a fixed value of 300 kg/m3 is for example used according to Amory et al. (2021, GMD). Have you tried to use for example as minimum value of 250kg/m3 for fresh snow density?
- line 141: what is the amount of measured precipitation which is missing by comparing the amount of ERA5 precipitation during the measured period vs the ERA5 precipitation during the missing period?
Citation: https://doi.org/10.5194/tc-2021-100-RC1 -
RC2: 'Comment on tc-2021-100', Nander Wever, 18 May 2021
The manuscript discusses the application of the physics based snow cover model Crocus at a field site in Greenland that is seasonally snow covered. First, the manuscript verifies representativeness of ERA5 using available weather station data. A new snow density parameterization scheme is proposed which should better account for the influence of wind on new snow density. Model output is compared to in-situ observations of snow density and specific surface area. Agreement is found to vary, which the authors attribute to missing snow transport by wind and water vapour transport in the Crocus model.
Principally, the manuscript is well written (it was a very pleasant read) with mostly clear and informative figures. Interesting field data is described, using, among others, the SnowMicroPen and IceCube instruments. I think the manuscript is the result of substantial work, for which I commend the authors.However, I do feel that the manuscript falls short of providing an adequate approach to address the research questions posed in the introduction. Firstly, the stated aims to validate the use of ERA5 and Crocus for Greenland cannot be achieved when considering only one field site. But more importantly, it is clear that the site is susceptible to snow transport by wind, which apparently is insufficiently described in the Crocus model. I also feel that it is an inadequate approach to focus on new snow density parameterizations in such cases, as I will explain below.
As shown in Fig. 9a/10a, it is clear that when Crocus is driven by ERA5 forcing data, there is only a small improvement from the new new snow density scheme (ERA5-CTRL vs ERA5-NSD). The differences between both schemes seem smaller than the discrepancy with observed density. It is also clear that the precipitation data from the in-situ measurements capture a lot of drifting snow (authors discuss this in L515-517). This likely also constitutes snow that has been eroded from the surface and not only precipitation. When Crocus is driven with in-situ data, the presence of drifting snow in the precipitation data leads to an overestimation of snow depth, even with 35% missing precipitation data which was set to 0 by the authors. More importantly, we can see that using in-situ precipitation gives a good agreement in snow density (Figs. 9b, 10b). Fig 10b illustrates that the default new snow scheme in Crocus actually performs more or less equally well as the proposed new scheme (even though based on Fig. 9b, the new scheme does work better). So this to me clearly indicates that it is the timing of when snow is actually incorporated into the snowpack that is important, rather than the new snow density formulation. When using in-situ precipitation data, it captures the drifting snow episodes that obviously occur during high wind speeds. Then, when the drifting snow is considered snow fall and is added to the Crocus model domain, it is assigned a higher density due to the high wind speed. This cannot be achieved when relying on ERA5 precipitation timings. To conclude: I think the results demonstrate that a better agreement in snow density is primarily achieved from the difference in timing between in-situ and ERA5 precipitation, rather than the choice of new snow density parameterization. Note that Fig. 5 suggests that wind speed was adequately simulated by ERA5 during the events before the field measurements, so I don't think that there is a problem with ERA5 wind speeds here.
The importance of timing when snow is incorporated into the snowpack has been addressed previously. In Groot Zwaaftink et al. (2013), a method is proposed to only add snow to the firn layer when there are episodes with high wind speed. This is based on the notion that in windy environments, low density snow will always be eroded at some point and deposited with higher density. It was also the motivation for the implementation of an erosion/deposition mechanism in the SNOWPACK model (Keenan et al. 2021).
So my judgement is that in spite of the large data acquisition and modelling efforts by the authors, the manuscript unfortunately is not addressing the problems at hand in a useful way, I think. For example, I think it would have been better to focus on the near surface densification, mentioned in L97-100 (see also Brun et al. 1997, Groot Zwaaftink 2013). I consider that a more important tuning knob than the new snow density scheme, since the near surface compaction by wind is implemented in Crocus precisely to describe the effect of drifting snow on the near surface density.
I wonder if the manuscript would be better suited as a brief communication, highlighting the issue of using reanalysis to drive snow models in windy environments. It is a clear demonstration that the snow stratigraphy in windy environments is impacted more by drifting snow episodes than by precipitation. Relying only on reanalysis precipitation without considering drifting snow will always fall short to correctly describe the stratigraphy. A new snow density parameterization applied when reanalysis provides precipitation will never be able to cover cases when the snow is actually incorporated into the snowpack during high wind speed conditions (which can occur even in the absence of precipitation). It probably would require some additional analysis of timing differences between precipitation recorded at the weather station and in ERA5.
Other comments:
Introduction: L50: sentences like "while other studies found deficiencies" need more explanation what deficiencies were found. Similarly L55 "Adjustments of the original model can help..." What adjustments are meant here?Fig 1 should be improved for readability. First of all, the caption claims that ERA5 grid cells are coloured green, whereas I'm under the impression that the crosses are coloured according to elevation (presumably the elevation in the ERA5 model?) Second, it's somewhat hard to see which crosses have a black outline, so why not use plus signs for those grid points, and replace the plus signs for the station and study site to a square and a circle for example? I think it's always best when different markers differ as much as possible (both shape and colour for example).
L95: what is the difference between "compaction" and "near surface densification"? I guess the latter is the compaction due to wind?
L138-140: following the discussion by the authors, the main source of discrepancy seems to be that the precipitation measurements are heavily impacted by drifting snow, rather than different spatial scales.
L610-615: I definitely think that when a newer version of the SMP was used in the field, the newer published calibrations should be favoured over Proksch et al. (2015). I think it is known that the Proksch et al. (2015) calibration does not hold for newer versions of the instrument. L155 should specify the version of the instrument used.
References:
- Brun, E., Martin, E., & Spiridonov, V. (1997). Coupling a multi-layered snow model with a GCM. Annals of Glaciology, 25, 66-72. doi:10.3189/S0260305500013811- Groot Zwaaftink, C. D., Cagnati, A., Crepaz, A., Fierz, C., Macelloni, G., Valt, M., and Lehning, M.: Event-driven deposition of snow on the Antarctic Plateau: analyzing field measurements with SNOWPACK, The Cryosphere, 7, 333–347, https://doi.org/10.5194/tc-7-333-2013, 2013.
- Keenan, E., Wever, N., Dattler, M., Lenaerts, J. T. M., Medley, B., Kuipers Munneke, P., and Reijmer, C.: Physics-based SNOWPACK model improves representation of near-surface Antarctic snow and firn density, The Cryosphere, 15, 1065–1085, https://doi.org/10.5194/tc-15-1065-2021, 2021.
Citation: https://doi.org/10.5194/tc-2021-100-RC2
Status: closed
-
RC1: 'Comment on tc-2021-100', Anonymous Referee #1, 07 May 2021
This paper presents first a validation of ERA5 over the Villum Research Station (NE Greenland) in the aim afterwards to force the CROCUS snowmodel with ERA5 and to evaluate its results with observations made at this station. While the paper is well written and uncertainties well discussed, the main scientific message, relevance and originality of this paper is poor for a journal such TC. ERA5 has already been evaluated over Greenland, CROCUS has already been applied in Arctic/Greenland and the fresh snow density improvements proposed in this paper seems to not improve a lot the comparison.
Moreover, this paper shows notably that using observed precipitation as forcing of CROCUS overestimates the snow height and that, although improvements in the prescribed fresh snow density, CROCUS still underestimates snow density. This suggests for me that:
=> as the modelled snow density is biased, both modelled and observed snowpack should be compared in water equivalent as observed density is available. A large part of discrepancies shown in Fig6 could be due to problems in the modelled snow density and not to overestimation/underestimation of prescribed precipitation as forcing. This could change the discussion in this paper if it is well the case. In Fig6, the modelled snow height can at least be corrected with respect to the snow density bias with observation and in this case, we could likely see a systematic underestimation of modelled snowpack using ERA5 as forcing.
=> as mentioned in the discussion, it is also very likely that this area is highly impacted by the wind erosion (increasing snow density and eroding the snow pack) as well as by the wind transport of the snow (biasing the snowfall measurements and potentially bringing snow from elsewhere) decreasing a lot the relevance of this study knowing that the blowing snow processes are not resolved here. Therefore, it is very likely that a large part of discrepancies between CROCUS and observations could just be due to the fact that this process is missing here and this does not necesserly mean that CROCUS is not working well in Arctic. Comparison with RCM (RACMO, MAR) taking into account these processes could help to evaluate the importance of these processes in this area.
Some more minor remarks:
- line78: I assume that this data is not assimilated into ERA5. Where is the closest weather station for which its observation is assimilated into ERA5?
- line 124: As the altitude of the ERA5 pixel is too high with respect to the station altitude, ERA5 is too cold and its surface pressure is too low. Converting all the precipitation to snow below 1°C will overestimate a lot the amount of snowfall with respect to rainfall knowing that ERA5 is already too cold.
- line 169, section 2.3: According to Fausto et al. (2018, A Snow Density Dataset for Improving Surface Boundary Conditions in Greenland Ice Sheet Firn Modeling, https://doi.org/10.3389/feart.2018.00051) using a fixed fresh snow density of 315 km/m3 is better than any parametrisation based on temperature and wind. In MARv3.11 model, a fixed value of 300 kg/m3 is for example used according to Amory et al. (2021, GMD). Have you tried to use for example as minimum value of 250kg/m3 for fresh snow density?
- line 141: what is the amount of measured precipitation which is missing by comparing the amount of ERA5 precipitation during the measured period vs the ERA5 precipitation during the missing period?
Citation: https://doi.org/10.5194/tc-2021-100-RC1 -
RC2: 'Comment on tc-2021-100', Nander Wever, 18 May 2021
The manuscript discusses the application of the physics based snow cover model Crocus at a field site in Greenland that is seasonally snow covered. First, the manuscript verifies representativeness of ERA5 using available weather station data. A new snow density parameterization scheme is proposed which should better account for the influence of wind on new snow density. Model output is compared to in-situ observations of snow density and specific surface area. Agreement is found to vary, which the authors attribute to missing snow transport by wind and water vapour transport in the Crocus model.
Principally, the manuscript is well written (it was a very pleasant read) with mostly clear and informative figures. Interesting field data is described, using, among others, the SnowMicroPen and IceCube instruments. I think the manuscript is the result of substantial work, for which I commend the authors.However, I do feel that the manuscript falls short of providing an adequate approach to address the research questions posed in the introduction. Firstly, the stated aims to validate the use of ERA5 and Crocus for Greenland cannot be achieved when considering only one field site. But more importantly, it is clear that the site is susceptible to snow transport by wind, which apparently is insufficiently described in the Crocus model. I also feel that it is an inadequate approach to focus on new snow density parameterizations in such cases, as I will explain below.
As shown in Fig. 9a/10a, it is clear that when Crocus is driven by ERA5 forcing data, there is only a small improvement from the new new snow density scheme (ERA5-CTRL vs ERA5-NSD). The differences between both schemes seem smaller than the discrepancy with observed density. It is also clear that the precipitation data from the in-situ measurements capture a lot of drifting snow (authors discuss this in L515-517). This likely also constitutes snow that has been eroded from the surface and not only precipitation. When Crocus is driven with in-situ data, the presence of drifting snow in the precipitation data leads to an overestimation of snow depth, even with 35% missing precipitation data which was set to 0 by the authors. More importantly, we can see that using in-situ precipitation gives a good agreement in snow density (Figs. 9b, 10b). Fig 10b illustrates that the default new snow scheme in Crocus actually performs more or less equally well as the proposed new scheme (even though based on Fig. 9b, the new scheme does work better). So this to me clearly indicates that it is the timing of when snow is actually incorporated into the snowpack that is important, rather than the new snow density formulation. When using in-situ precipitation data, it captures the drifting snow episodes that obviously occur during high wind speeds. Then, when the drifting snow is considered snow fall and is added to the Crocus model domain, it is assigned a higher density due to the high wind speed. This cannot be achieved when relying on ERA5 precipitation timings. To conclude: I think the results demonstrate that a better agreement in snow density is primarily achieved from the difference in timing between in-situ and ERA5 precipitation, rather than the choice of new snow density parameterization. Note that Fig. 5 suggests that wind speed was adequately simulated by ERA5 during the events before the field measurements, so I don't think that there is a problem with ERA5 wind speeds here.
The importance of timing when snow is incorporated into the snowpack has been addressed previously. In Groot Zwaaftink et al. (2013), a method is proposed to only add snow to the firn layer when there are episodes with high wind speed. This is based on the notion that in windy environments, low density snow will always be eroded at some point and deposited with higher density. It was also the motivation for the implementation of an erosion/deposition mechanism in the SNOWPACK model (Keenan et al. 2021).
So my judgement is that in spite of the large data acquisition and modelling efforts by the authors, the manuscript unfortunately is not addressing the problems at hand in a useful way, I think. For example, I think it would have been better to focus on the near surface densification, mentioned in L97-100 (see also Brun et al. 1997, Groot Zwaaftink 2013). I consider that a more important tuning knob than the new snow density scheme, since the near surface compaction by wind is implemented in Crocus precisely to describe the effect of drifting snow on the near surface density.
I wonder if the manuscript would be better suited as a brief communication, highlighting the issue of using reanalysis to drive snow models in windy environments. It is a clear demonstration that the snow stratigraphy in windy environments is impacted more by drifting snow episodes than by precipitation. Relying only on reanalysis precipitation without considering drifting snow will always fall short to correctly describe the stratigraphy. A new snow density parameterization applied when reanalysis provides precipitation will never be able to cover cases when the snow is actually incorporated into the snowpack during high wind speed conditions (which can occur even in the absence of precipitation). It probably would require some additional analysis of timing differences between precipitation recorded at the weather station and in ERA5.
Other comments:
Introduction: L50: sentences like "while other studies found deficiencies" need more explanation what deficiencies were found. Similarly L55 "Adjustments of the original model can help..." What adjustments are meant here?Fig 1 should be improved for readability. First of all, the caption claims that ERA5 grid cells are coloured green, whereas I'm under the impression that the crosses are coloured according to elevation (presumably the elevation in the ERA5 model?) Second, it's somewhat hard to see which crosses have a black outline, so why not use plus signs for those grid points, and replace the plus signs for the station and study site to a square and a circle for example? I think it's always best when different markers differ as much as possible (both shape and colour for example).
L95: what is the difference between "compaction" and "near surface densification"? I guess the latter is the compaction due to wind?
L138-140: following the discussion by the authors, the main source of discrepancy seems to be that the precipitation measurements are heavily impacted by drifting snow, rather than different spatial scales.
L610-615: I definitely think that when a newer version of the SMP was used in the field, the newer published calibrations should be favoured over Proksch et al. (2015). I think it is known that the Proksch et al. (2015) calibration does not hold for newer versions of the instrument. L155 should specify the version of the instrument used.
References:
- Brun, E., Martin, E., & Spiridonov, V. (1997). Coupling a multi-layered snow model with a GCM. Annals of Glaciology, 25, 66-72. doi:10.3189/S0260305500013811- Groot Zwaaftink, C. D., Cagnati, A., Crepaz, A., Fierz, C., Macelloni, G., Valt, M., and Lehning, M.: Event-driven deposition of snow on the Antarctic Plateau: analyzing field measurements with SNOWPACK, The Cryosphere, 7, 333–347, https://doi.org/10.5194/tc-7-333-2013, 2013.
- Keenan, E., Wever, N., Dattler, M., Lenaerts, J. T. M., Medley, B., Kuipers Munneke, P., and Reijmer, C.: Physics-based SNOWPACK model improves representation of near-surface Antarctic snow and firn density, The Cryosphere, 15, 1065–1085, https://doi.org/10.5194/tc-15-1065-2021, 2021.
Citation: https://doi.org/10.5194/tc-2021-100-RC2
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
1,056 | 496 | 56 | 1,608 | 60 | 53 |
- HTML: 1,056
- PDF: 496
- XML: 56
- Total: 1,608
- BibTeX: 60
- EndNote: 53
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1
Cited
2 citations as recorded by crossref.
- Snow and meteorological conditions at Villum Research Station, Northeast Greenland: on the adequacy of using atmospheric reanalysis for detailed snow simulations D. Krampe et al. 10.3389/feart.2023.1053918
- Multi-physics ensemble modelling of Arctic tundra snowpack properties G. Woolley et al. 10.5194/tc-18-5685-2024