the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Modelling snowpack on ice surfaces with the ORCHIDEE land surface model: Application to the Greenland ice sheet
Abstract. Current climate warming is accelerating mass loss from glaciers and ice sheets. In Greenland, the rates of mass changes are now dominated by changes in surface mass balance (SMB) due to increased surface melting. To improve the future sea-level rise projections, it is therefore critical to have an accurate estimate of the SMB, which depends on the representation of the processes occurring within the snowpack. The snow scheme (ES) implemented in the land surface model ORCHIDEE has not yet been adapted to ice-covered areas. Here, we present the preliminary developments we made to apply the ES model to glaciers and ice sheets. Our analysis mainly concerns the model’s ability to represent ablation-related processes. At the regional scale, our results are compared to the MAR regional atmospheric model outputs and to MODIS albedo retrievals.
Using different albedo parameterizations, we performed offline ES simulations forced by the MAR model over the 2000–2019 period. Our results reveal a strong sensitivity of the modeled SMB components to the albedo parameterization. Results inferred with albedo parameters obtained with a manual tuning approach present a very good agreement with the MAR outputs. Conversely, with the albedo parameterization used in the standard ORCHIDEE version, runoff and sublimation were underestimated. We also tested parameters found from a previous data assimilation experiment calibrating the ablation processes using MODIS snow albedo. While these parameters greatly improve the modelled albedo over the entire ice sheet, they degrade the other model outputs compared to those obtained with the manually-tuned approach. This is likely due to the model overfitting to the calibration albedo dataset without any constraint applied to the other processes controlling the state of the snowpack. This underlines the need for performing a “multi-objective” optimisation using auxiliary observations related to snowpack internal processes. Although there is still room for further improvements, the developments reported in the present study constitute an important advance in assessing the Greenland SMB with possible extension to mountain glaciers or the Antarctic ice sheet.
-
Notice on discussion status
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
-
Preprint
(2640 KB)
-
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(2640 KB) - Metadata XML
- BibTeX
- EndNote
- Final revised paper
Journal article(s) based on this preprint
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2024-285', Anonymous Referee #1, 05 Apr 2024
My comments are given in the attached .pdf file.
-
AC1: 'Reply on RC1', Sylvie Charbit, 16 Jun 2024
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2024/egusphere-2024-285/egusphere-2024-285-AC1-supplement.pdf
-
AC1: 'Reply on RC1', Sylvie Charbit, 16 Jun 2024
-
RC2: 'Comment on egusphere-2024-285', Anonymous Referee #2, 17 Apr 2024
This paper describes new modeling of Greenland snow and ice in the ORCHIDEE surface model. It is well written and pleasant to read. The work is substantial and certainly deserves to be published. The idea of using MAR as a forcing and comparing it to its outputs is very pertinent. Nevertheless, it seems to me that certain modeling choices need to be discussed and model validation a little improved.
Major comments :
1 - For Figures 1 to 4, you give the raw distributions of MAR and each simulation. But to show differences would be also interesting. In addition, for the others (Figure 5 to 8), you give the differences, but not the raw distributions. As it is possible for you to have additional information, it would be nice to provide all the details, i.e. both the raw distributions and the differences for each variable analyzed (and thus for each figure 1 to 8).
2 - By the way, regarding this validation, why not provide the PDF of each quantity, i.e. compare the MAR (or MODIS) PDF with the ORCHIDEE PDFs ? When you look at Figures 1 to 8, it's hard to see whether one version is much better than another. An objective method of comparison is missing. Comparing PDFs could be a solution but there is perhaps another ways.
3 - It seems that spatial statistics (correlation, etc.) are missing for each field analyzed. In other words, two or three objective criteria to determine whether in Figure 2 (for example) the OPT-12L spatial distribution (e) is better than the others compared to MAR (a). That is, each panel (b, c , d, e) should have its spatial correlation (spatial rmse, etc.) with MAR. The fact that, for example, the spatial distribution of OPT-3L refreezing (Figure 3e) is closer to MAR (Figure 3a) is not trivial to see with the naked eye. Anyway, I hope this comment is understandable. You lack objective statistical criteria in your assessment of all the figures showing comparisons of spatial distribution. The simple Table 2 obtained via a spatial average is not enough.
4 - The modeling choices made could have been discussed. In particular, the parametrization of snow albedo. Moreover, I don't understand why this new parametrization compared to Wang et al. (2013) is in section 2.1 (existing parm) and not rather in section 2.2 (new param). This new parametrization is a bit outdated today when there are more robust parametrizations in land surface models accounting for spectral albedo and solar absorption calculation as in CLM with SNICAR (Flanner and Zender 2006) or ISBA-ES (Decharme et al. 2016). What's more, this more robust representation already exists for ES (Decharme et al. 2016). Why not use it ? This choice is debatable in view of the importance of albedo on the SMB. Please discuss about that.
Other comments:
1 - Independent MAR observations of the Greenland SMB based on GRACE data could have been used in section 5.3 (Schlegel et al. 2016, Wang et al. 2024).
2 - From what I understand, some parameterizations that are in section 2.1 are new compared to Wang et al. (2013), and should therefore be in section 2.2. : snow fraction, albedo
3 - On the implementation of the ice layer (section 2.2.2), why didn't you use ES directly to model this. On line 293, you say that snow density is limited to 750kg/m³. But if you had raised this limit to 900 or 950kg/m³, it seems to me that all the "snow" equations converge to "ice" equations, at least that's what comes out when we compare equation 26 to 20, 27 to 21, etc. In theory, if ES is done right, snow that has reached a certain density should be able to become ice. It would then be sufficient to initialize the height and density of the snowpack accordingly (e.g. the last 6 layers with an ice density and a total height of 100m for example). I don’t know if it’s possible but this could be discussed.
4 - Line 629 - 632: I understand here that the improved runoff modeling in OPT-12L would not be due to bias compensation. Well, I'm really not sure. What I understand from looking at your results is that to improve runoff compared to MAR, you need to set an albedo lower than MAR (Figure 6e), which inevitably induces a surface temperature (and surely an internal temperature of the snowpack) that is too high (Figure 7e) compared to MAR. I have the impression that this is also what Figure (8f) reveals. So, to obtain the same runoff than MAR, ORCHIDEE-ICE have to simulate a lower albedo than MAR to capture more energy. If it is true, it is perhaps due to the non representation of solar absorption by snow or a poor simulation of snowpack density.
5 - This last remark also underlines the fact that other important variables concerning the internal properties of the snowpack could be shown/analyzed, such as the temperature and density of the simulated snowpack compared with MAR, etc. This would enable a better understanding of the processes involved in the improvements related to one or another process claimed by the authors.
References :
Flanner, M. G., and C. S. Zender (2006), Linking snowpack microphysics and albedo evolution, J. Geophys. Res., 111, D12208, doi:10.1029/2005JD006834.
Schlegel, N.-J., Wiese, D. N., Larour, E. Y., Watkins, M. M., Box, J. E., Fettweis, X., and van den Broeke, M. R.: Application of GRACE to the assessment of model-based estimates of monthly Greenland Ice Sheet mass balance (2003–2012), The Cryosphere, 10, 1965–1989, https://doi.org/10.5194/tc-10-1965-2016, 2016.
Wang W., Yunzhong Shen, Qiujie Chen, Fengwei Wang, High-resolution mascon solutions reveal glacier-scale mass changes over the Greenland Ice Sheet from 2002 to 2022, Geophysical Journal International, Volume 236, Issue 1, January 2024, Pages 494–515, https://doi.org/10.1093/gji/ggad439
Citation: https://doi.org/10.5194/egusphere-2024-285-RC2 -
AC2: 'Reply on RC2', Sylvie Charbit, 16 Jun 2024
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2024/egusphere-2024-285/egusphere-2024-285-AC2-supplement.pdf
-
AC2: 'Reply on RC2', Sylvie Charbit, 16 Jun 2024
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2024-285', Anonymous Referee #1, 05 Apr 2024
My comments are given in the attached .pdf file.
-
AC1: 'Reply on RC1', Sylvie Charbit, 16 Jun 2024
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2024/egusphere-2024-285/egusphere-2024-285-AC1-supplement.pdf
-
AC1: 'Reply on RC1', Sylvie Charbit, 16 Jun 2024
-
RC2: 'Comment on egusphere-2024-285', Anonymous Referee #2, 17 Apr 2024
This paper describes new modeling of Greenland snow and ice in the ORCHIDEE surface model. It is well written and pleasant to read. The work is substantial and certainly deserves to be published. The idea of using MAR as a forcing and comparing it to its outputs is very pertinent. Nevertheless, it seems to me that certain modeling choices need to be discussed and model validation a little improved.
Major comments :
1 - For Figures 1 to 4, you give the raw distributions of MAR and each simulation. But to show differences would be also interesting. In addition, for the others (Figure 5 to 8), you give the differences, but not the raw distributions. As it is possible for you to have additional information, it would be nice to provide all the details, i.e. both the raw distributions and the differences for each variable analyzed (and thus for each figure 1 to 8).
2 - By the way, regarding this validation, why not provide the PDF of each quantity, i.e. compare the MAR (or MODIS) PDF with the ORCHIDEE PDFs ? When you look at Figures 1 to 8, it's hard to see whether one version is much better than another. An objective method of comparison is missing. Comparing PDFs could be a solution but there is perhaps another ways.
3 - It seems that spatial statistics (correlation, etc.) are missing for each field analyzed. In other words, two or three objective criteria to determine whether in Figure 2 (for example) the OPT-12L spatial distribution (e) is better than the others compared to MAR (a). That is, each panel (b, c , d, e) should have its spatial correlation (spatial rmse, etc.) with MAR. The fact that, for example, the spatial distribution of OPT-3L refreezing (Figure 3e) is closer to MAR (Figure 3a) is not trivial to see with the naked eye. Anyway, I hope this comment is understandable. You lack objective statistical criteria in your assessment of all the figures showing comparisons of spatial distribution. The simple Table 2 obtained via a spatial average is not enough.
4 - The modeling choices made could have been discussed. In particular, the parametrization of snow albedo. Moreover, I don't understand why this new parametrization compared to Wang et al. (2013) is in section 2.1 (existing parm) and not rather in section 2.2 (new param). This new parametrization is a bit outdated today when there are more robust parametrizations in land surface models accounting for spectral albedo and solar absorption calculation as in CLM with SNICAR (Flanner and Zender 2006) or ISBA-ES (Decharme et al. 2016). What's more, this more robust representation already exists for ES (Decharme et al. 2016). Why not use it ? This choice is debatable in view of the importance of albedo on the SMB. Please discuss about that.
Other comments:
1 - Independent MAR observations of the Greenland SMB based on GRACE data could have been used in section 5.3 (Schlegel et al. 2016, Wang et al. 2024).
2 - From what I understand, some parameterizations that are in section 2.1 are new compared to Wang et al. (2013), and should therefore be in section 2.2. : snow fraction, albedo
3 - On the implementation of the ice layer (section 2.2.2), why didn't you use ES directly to model this. On line 293, you say that snow density is limited to 750kg/m³. But if you had raised this limit to 900 or 950kg/m³, it seems to me that all the "snow" equations converge to "ice" equations, at least that's what comes out when we compare equation 26 to 20, 27 to 21, etc. In theory, if ES is done right, snow that has reached a certain density should be able to become ice. It would then be sufficient to initialize the height and density of the snowpack accordingly (e.g. the last 6 layers with an ice density and a total height of 100m for example). I don’t know if it’s possible but this could be discussed.
4 - Line 629 - 632: I understand here that the improved runoff modeling in OPT-12L would not be due to bias compensation. Well, I'm really not sure. What I understand from looking at your results is that to improve runoff compared to MAR, you need to set an albedo lower than MAR (Figure 6e), which inevitably induces a surface temperature (and surely an internal temperature of the snowpack) that is too high (Figure 7e) compared to MAR. I have the impression that this is also what Figure (8f) reveals. So, to obtain the same runoff than MAR, ORCHIDEE-ICE have to simulate a lower albedo than MAR to capture more energy. If it is true, it is perhaps due to the non representation of solar absorption by snow or a poor simulation of snowpack density.
5 - This last remark also underlines the fact that other important variables concerning the internal properties of the snowpack could be shown/analyzed, such as the temperature and density of the simulated snowpack compared with MAR, etc. This would enable a better understanding of the processes involved in the improvements related to one or another process claimed by the authors.
References :
Flanner, M. G., and C. S. Zender (2006), Linking snowpack microphysics and albedo evolution, J. Geophys. Res., 111, D12208, doi:10.1029/2005JD006834.
Schlegel, N.-J., Wiese, D. N., Larour, E. Y., Watkins, M. M., Box, J. E., Fettweis, X., and van den Broeke, M. R.: Application of GRACE to the assessment of model-based estimates of monthly Greenland Ice Sheet mass balance (2003–2012), The Cryosphere, 10, 1965–1989, https://doi.org/10.5194/tc-10-1965-2016, 2016.
Wang W., Yunzhong Shen, Qiujie Chen, Fengwei Wang, High-resolution mascon solutions reveal glacier-scale mass changes over the Greenland Ice Sheet from 2002 to 2022, Geophysical Journal International, Volume 236, Issue 1, January 2024, Pages 494–515, https://doi.org/10.1093/gji/ggad439
Citation: https://doi.org/10.5194/egusphere-2024-285-RC2 -
AC2: 'Reply on RC2', Sylvie Charbit, 16 Jun 2024
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2024/egusphere-2024-285/egusphere-2024-285-AC2-supplement.pdf
-
AC2: 'Reply on RC2', Sylvie Charbit, 16 Jun 2024
Peer review completion
Journal article(s) based on this preprint
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
338 | 126 | 34 | 498 | 24 | 20 |
- HTML: 338
- PDF: 126
- XML: 34
- Total: 498
- BibTeX: 24
- EndNote: 20
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1
Cited
Sylvie Charbit
Christophe Dumas
Fabienne Maignan
Catherine Ottlé
Nina Raoult
Xavier Fettweis
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(2640 KB) - Metadata XML