Integrating a global glacier model into local hydrological modeling: Impacts on melt contributions

Berg, Justine; Horton, Pascal; Kauzlaric, Martina; von der Esch, Alexandra; Schaefli, Bettina

doi:10.5194/egusphere-2026-439

Preprints

https://doi.org/10.5194/egusphere-2026-439

Preprints

25 Feb 2026

| 25 Feb 2026

Integrating a global glacier model into local hydrological modeling: Impacts on melt contributions

Justine Berg, Pascal Horton, Martina Kauzlaric, Alexandra von der Esch, and Bettina Schaefli

Abstract. Snow and glacier meltwater are critical water resources for mountain regions, yet their accurate representation in hydrological models remains challenging. As climate change alters the timing and magnitude of melt contributions, accurate partitioning between snow and glacier sources becomes increasingly important. To address this challenge, this study couples the hydrological model HBV (as implemented in the Raven modeling framework) with the Global Glacier Evolution Model GloGEM and integrates a snow redistribution scheme for application across 14 glaciated headwater catchments in Switzerland. At this local catchment scale, we investigate whether glacier constraints and the addition of snow redistribution reduce parameter equifinality and increase the reliability of melt contribution estimates. Simulations are evaluated against observed streamflow, gridded snow water equivalent data, and glacier melt data based on observed mass balance. The gravitational snow redistribution algorithm successfully prevents unrealistic high-elevation snow accumulation and improves catchment average SWE simulation performance. Uncoupled HBV configurations outperform coupled HBV-GloGEM setups according to streamflow metrics. However, this superior performance is achieved by simulating glacier melt rates exceeding GloGEM estimates and glacier storage change data by factors of 2–3 in some catchments, effectively using glacier ice to offset precipitation biases in forcing data, which can be critical for climate change impact studies. Glacier melt contributions can vary by nearly an order of magnitude among best-performing parameter sets, highlighting severe parameter equifinality. Coupling with GloGEM, calibrated using glacier-specific geodetic mass balance, produces glacier melt consistent with observations and substantially improves identifiability of both glacier and snowmelt contributions. Despite lower streamflow performance metrics, glacier melt constraints prevent compensatory errors that would compromise projections as melt dynamics shift under climate change. Applying this framework to future climate scenarios and integrating additional constraints such as snow observations may further improve the reliability of melt contribution projections.

Received: 26 Jan 2026 – Discussion started: 25 Feb 2026

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.

Download & links

Preprint (PDF, 10118 KB)

Supplement (2494 KB)

Download & links

Preprint (10118 KB)
Metadata XML
Supplement (2494 KB)
BibTeX
EndNote

Justine Berg, Pascal Horton, Martina Kauzlaric, Alexandra von der Esch, and Bettina Schaefli

Status: closed

RC1:
'Comment on egusphere-2026-439', Anonymous Referee #1, 07 Mar 2026
This manuscript presents a solid study that integrates the global glacier model GloGEM with the HBV model implemented in the Raven framework, together with a snow redistribution scheme, across 14 glacierized headwater catchments in Switzerland. The study offers a systematic assessment of how glacier constraints, snow redistribution, and precipitation correction affect melt partitioning and parameter behavior at the catchment scale. Overall, the work is carefully conducted and valuable. That said, I have several concerns regarding the experimental design and the novelty of the main findings. My main comments are as follows:
The main contribution of the paper seems stronger as a systematic regional assessment than as a genuinely new methodological advance. The study is solid and useful, but the novelty of the main findings appears somewhat limited. The main messages of the paper are that (i) good streamflow performance does not necessarily imply realistic cryospheric process representation, (ii) uncoupled HBV can compensate forcing biases through excessive glacier melt, and (iii) introducing a precipitation correction factor can aggravate equifinality and uncertainty. These are important points, but they are not new in cold-region hydrology and glacier-hydrology modeling. In this sense, the main value of the manuscript seems to lie in the systematic multi-catchment comparison across 14 headwater basins, rather than in delivering new methodological or process-level insight.

The role of the precipitation correction factor (PCF) needs to be explained more clearly. The manuscript lists PRECIP_CORR as a precipitation correction factor in Table 2 and treats it as a separate configuration dimension in Table 3 (e.g., HBV-PCF and HBV-GloGEM-PCF). However, its hydrological meaning and implementation remain insufficiently explained. At present, the manuscript provides a label, but does not clearly explain what exactly is corrected, how the factor operates within the HBV/Raven setup, or what physical bias it is intended to represent. It is also unclear why calibration of this parameter is treated as a dedicated scenario rather than simply as one additional calibration parameter.

The contrast between the two model setups may not primarily reflect “whether GloGEM is coupled,” but rather whether glacier information is used to constrain the glacier component. In the current design, the coupled setup benefits from glacier-specific calibration against geodetic mass balance, whereas the original HBV setup is calibrated only against streamflow. This makes it difficult to attribute the differences solely to coupling with GloGEM. In practical terms, both the GloGEM-based runoff and the melt generated by the simplified HBV glacier routine ultimately act as inputs to the hydrological simulation. It is therefore plausible that, if the native HBV glacier routine were also constrained using glacier data, it could produce more comparable glacier-melt estimates. The manuscript should discuss this point more explicitly and clarify whether the key distinction is really coupling itself, or the use of glacier observations to constrain the glacier contribution.

The snow redistribution module should be introduced and visualized more clearly. The manuscript explains why snow redistribution is needed and briefly describes how it is implemented. However, its role is not clearly shown in the workflow figure, which makes the overall methodology harder to follow. I suggest either adding this component explicitly to Figure 2 or clarifying its position and function more clearly in the text.

The treatment of snow-covered area (SCA) remains unclear, and the lack of SCA-based evaluation weakens confidence in the snow and melt partitioning results. I may have overlooked it, but I could not find a clear explanation of how snow-covered area is determined in the model. In addition, the model does not appear to be evaluated against SCA data. This limits confidence in the snowmelt simulations. Even if catchment-average SWE is reproduced reasonably well, the spatial distribution of snow and the timing and location of melt may still remain uncertain without an independent constraint on snow extent. I therefore recommend that the authors clarify whether SCA is represented or diagnosed in the current workflow, and discuss the implications of not evaluating SCA for the interpretation of snow storage and melt contributions.

Minor issues
The naming of the model configurations in Table 3 is somewhat confusing. In particular, the HBV-PCF and HBV-GloGEM-PCF configurations also include snow redistribution, but this is not reflected in their names. A more consistent naming convention, or a clearer explanation when these configurations are first introduced, would improve readability.

In Table 2, BASEFLOW_COEFF is listed twice for the fast and slow reservoirs. Although this may be correct in the model implementation, the presentation could be clearer for readers who are not familiar with the Raven/HBV setup.

Please clarify explicitly whether SWE was used only for validation or also during calibration. From the current text, SWE appears to have been used for validation rather than calibration. Since SWE results are shown prominently, this point should be stated unambiguously.

There are several minor language and typographical issues that should be corrected during revision. For example, “largly” should be “largely” (line 85), “relativity similar” should be “relatively similar” (line 230). Figure 2 appears to contain a typo in the label “Precipiation per HRU”, which should be corrected.
Citation: https://doi.org/10.5194/egusphere-2026-439-RC1
- AC1: 'Reply on RC1', Justine Berg, 16 Mar 2026
  
  Dear Reviewer, thank you for your thorough review of our manuscript. We appreciate your valuable comments. Please find our point-by-point response in the attached document.
  
  Citation: https://doi.org/10.5194/egusphere-2026-439-AC1
RC2:
'Comment on egusphere-2026-439', Anonymous Referee #2, 07 Jun 2026

The manuscript presents a systematic modelling experiment in which the authors couple the global glacier evolution model GloGEM with the HBV model implemented in the Raven framework, and assess the added value of glacier constraints, snow redistribution, and precipitation correction for simulating streamflow and partitioning snowmelt and glacier-melt contributions across 14 glacierized Swiss headwater catchments. The topic is relevant, and the paper addresses the fact that good streamflow performance does not necessarily imply a realistic representation of snow and glacier melt processes. Overall, I find the study valuable, but I think that some methodological choices need to be clarified before publication. The major issues in my opinion are:
Role and validation of snow redistribution: snow redistrbution appears to have an important influence on the partitioning between snowmelt and glacier melt. In particular, it may affect how long glacierized areas remain snow-covered, and therefore how much ice melt is simulated. Given this importance, the snow redistribution scheme should be described more clearly, including how it operates within the semi-distributed HRU structure. I also think that an independent evaluation would be useful, for example using maps snow-covered area from satellite observations. Catchment-average SWE alone may not be sufficient to verify whether the spatial distribution and timing of snowmelt are realistic.
Routing of GloGEM meltwater: the treatment of glacier runoff in the coupled configuration is not fully clear to me. The manuscript states that GloGEM runoff bypasses intermediate storage components and is directly added to the runoff–streamflow transfer function. However, the native HBV glacier module seems to include a glacier storage coefficient controlling meltwater release. Are glacier melt contributions routed in a comparable way in the HBV and HBV-GloGEM configurations, or does the coupling remove part of the glacier storage/routing representation? Since this may affect the timing of simulated discharge, especially during the melt season, this point should be clarified and justified.
Precipitation correction factor: the hydrological meaning and calibration of the precipitation correction factor should be explained more clearly (i.e., is this factor intended to correct precipitation undercatch, orographic effects, winter precipitation bias,etc.). Is this PCF different from the fc,prec (line 205) calibrated for the GloGEM? How do these two parameters are linked? If they are different, does this means that there is a precipitation correction in all the HBV-GloGEM coupled configuration? Would this change the volumes reaching the catcment between configurations?
Initialization of model states: the initialization and spin-up of the different configurations should be clarified. Are all configurations initialized with the same SWE and storage conditions? Is SPASS used only for validation or also for initialization? Why not for calibration? In Figure 5b, the different configurations already seem to show different SWE values at the beginning of the validation period. This initial offset should be explained, as it may affect the interpretation of SWE evolution and “snow towers” development.
“Model parameters” section: this section is very short and the entire paper could benefit from a more extended description (and discussion in the “discussion” section) of the changes in model parameters (not only the melt factor) among the different configurations and how such differences are related to model overcompensation to fit the streamflow calibration objective. I see the authors have included figures in the supplement but I would appreciate a more extended presentation in the text.

Citation: https://doi.org/10.5194/egusphere-2026-439-RC2
- AC2: 'Reply on RC2', Justine Berg, 23 Jun 2026
  
  Dear Reviewer, thank you for your thorough and constructive review of our manuscript. We appreciate your comments and respond to each of them in the attached point-by-point document.
  
  Citation: https://doi.org/10.5194/egusphere-2026-439-AC2

Status: closed

RC1:
'Comment on egusphere-2026-439', Anonymous Referee #1, 07 Mar 2026
This manuscript presents a solid study that integrates the global glacier model GloGEM with the HBV model implemented in the Raven framework, together with a snow redistribution scheme, across 14 glacierized headwater catchments in Switzerland. The study offers a systematic assessment of how glacier constraints, snow redistribution, and precipitation correction affect melt partitioning and parameter behavior at the catchment scale. Overall, the work is carefully conducted and valuable. That said, I have several concerns regarding the experimental design and the novelty of the main findings. My main comments are as follows:
The main contribution of the paper seems stronger as a systematic regional assessment than as a genuinely new methodological advance. The study is solid and useful, but the novelty of the main findings appears somewhat limited. The main messages of the paper are that (i) good streamflow performance does not necessarily imply realistic cryospheric process representation, (ii) uncoupled HBV can compensate forcing biases through excessive glacier melt, and (iii) introducing a precipitation correction factor can aggravate equifinality and uncertainty. These are important points, but they are not new in cold-region hydrology and glacier-hydrology modeling. In this sense, the main value of the manuscript seems to lie in the systematic multi-catchment comparison across 14 headwater basins, rather than in delivering new methodological or process-level insight.

The role of the precipitation correction factor (PCF) needs to be explained more clearly. The manuscript lists PRECIP_CORR as a precipitation correction factor in Table 2 and treats it as a separate configuration dimension in Table 3 (e.g., HBV-PCF and HBV-GloGEM-PCF). However, its hydrological meaning and implementation remain insufficiently explained. At present, the manuscript provides a label, but does not clearly explain what exactly is corrected, how the factor operates within the HBV/Raven setup, or what physical bias it is intended to represent. It is also unclear why calibration of this parameter is treated as a dedicated scenario rather than simply as one additional calibration parameter.

The contrast between the two model setups may not primarily reflect “whether GloGEM is coupled,” but rather whether glacier information is used to constrain the glacier component. In the current design, the coupled setup benefits from glacier-specific calibration against geodetic mass balance, whereas the original HBV setup is calibrated only against streamflow. This makes it difficult to attribute the differences solely to coupling with GloGEM. In practical terms, both the GloGEM-based runoff and the melt generated by the simplified HBV glacier routine ultimately act as inputs to the hydrological simulation. It is therefore plausible that, if the native HBV glacier routine were also constrained using glacier data, it could produce more comparable glacier-melt estimates. The manuscript should discuss this point more explicitly and clarify whether the key distinction is really coupling itself, or the use of glacier observations to constrain the glacier contribution.

The snow redistribution module should be introduced and visualized more clearly. The manuscript explains why snow redistribution is needed and briefly describes how it is implemented. However, its role is not clearly shown in the workflow figure, which makes the overall methodology harder to follow. I suggest either adding this component explicitly to Figure 2 or clarifying its position and function more clearly in the text.

The treatment of snow-covered area (SCA) remains unclear, and the lack of SCA-based evaluation weakens confidence in the snow and melt partitioning results. I may have overlooked it, but I could not find a clear explanation of how snow-covered area is determined in the model. In addition, the model does not appear to be evaluated against SCA data. This limits confidence in the snowmelt simulations. Even if catchment-average SWE is reproduced reasonably well, the spatial distribution of snow and the timing and location of melt may still remain uncertain without an independent constraint on snow extent. I therefore recommend that the authors clarify whether SCA is represented or diagnosed in the current workflow, and discuss the implications of not evaluating SCA for the interpretation of snow storage and melt contributions.

Minor issues
The naming of the model configurations in Table 3 is somewhat confusing. In particular, the HBV-PCF and HBV-GloGEM-PCF configurations also include snow redistribution, but this is not reflected in their names. A more consistent naming convention, or a clearer explanation when these configurations are first introduced, would improve readability.

In Table 2, BASEFLOW_COEFF is listed twice for the fast and slow reservoirs. Although this may be correct in the model implementation, the presentation could be clearer for readers who are not familiar with the Raven/HBV setup.

Please clarify explicitly whether SWE was used only for validation or also during calibration. From the current text, SWE appears to have been used for validation rather than calibration. Since SWE results are shown prominently, this point should be stated unambiguously.

There are several minor language and typographical issues that should be corrected during revision. For example, “largly” should be “largely” (line 85), “relativity similar” should be “relatively similar” (line 230). Figure 2 appears to contain a typo in the label “Precipiation per HRU”, which should be corrected.
Citation: https://doi.org/10.5194/egusphere-2026-439-RC1
- AC1: 'Reply on RC1', Justine Berg, 16 Mar 2026
  
  Dear Reviewer, thank you for your thorough review of our manuscript. We appreciate your valuable comments. Please find our point-by-point response in the attached document.
  
  Citation: https://doi.org/10.5194/egusphere-2026-439-AC1
RC2:
'Comment on egusphere-2026-439', Anonymous Referee #2, 07 Jun 2026

The manuscript presents a systematic modelling experiment in which the authors couple the global glacier evolution model GloGEM with the HBV model implemented in the Raven framework, and assess the added value of glacier constraints, snow redistribution, and precipitation correction for simulating streamflow and partitioning snowmelt and glacier-melt contributions across 14 glacierized Swiss headwater catchments. The topic is relevant, and the paper addresses the fact that good streamflow performance does not necessarily imply a realistic representation of snow and glacier melt processes. Overall, I find the study valuable, but I think that some methodological choices need to be clarified before publication. The major issues in my opinion are:
Role and validation of snow redistribution: snow redistrbution appears to have an important influence on the partitioning between snowmelt and glacier melt. In particular, it may affect how long glacierized areas remain snow-covered, and therefore how much ice melt is simulated. Given this importance, the snow redistribution scheme should be described more clearly, including how it operates within the semi-distributed HRU structure. I also think that an independent evaluation would be useful, for example using maps snow-covered area from satellite observations. Catchment-average SWE alone may not be sufficient to verify whether the spatial distribution and timing of snowmelt are realistic.
Routing of GloGEM meltwater: the treatment of glacier runoff in the coupled configuration is not fully clear to me. The manuscript states that GloGEM runoff bypasses intermediate storage components and is directly added to the runoff–streamflow transfer function. However, the native HBV glacier module seems to include a glacier storage coefficient controlling meltwater release. Are glacier melt contributions routed in a comparable way in the HBV and HBV-GloGEM configurations, or does the coupling remove part of the glacier storage/routing representation? Since this may affect the timing of simulated discharge, especially during the melt season, this point should be clarified and justified.
Precipitation correction factor: the hydrological meaning and calibration of the precipitation correction factor should be explained more clearly (i.e., is this factor intended to correct precipitation undercatch, orographic effects, winter precipitation bias,etc.). Is this PCF different from the fc,prec (line 205) calibrated for the GloGEM? How do these two parameters are linked? If they are different, does this means that there is a precipitation correction in all the HBV-GloGEM coupled configuration? Would this change the volumes reaching the catcment between configurations?
Initialization of model states: the initialization and spin-up of the different configurations should be clarified. Are all configurations initialized with the same SWE and storage conditions? Is SPASS used only for validation or also for initialization? Why not for calibration? In Figure 5b, the different configurations already seem to show different SWE values at the beginning of the validation period. This initial offset should be explained, as it may affect the interpretation of SWE evolution and “snow towers” development.
“Model parameters” section: this section is very short and the entire paper could benefit from a more extended description (and discussion in the “discussion” section) of the changes in model parameters (not only the melt factor) among the different configurations and how such differences are related to model overcompensation to fit the streamflow calibration objective. I see the authors have included figures in the supplement but I would appreciate a more extended presentation in the text.

Citation: https://doi.org/10.5194/egusphere-2026-439-RC2
- AC2: 'Reply on RC2', Justine Berg, 23 Jun 2026
  
  Dear Reviewer, thank you for your thorough and constructive review of our manuscript. We appreciate your comments and respond to each of them in the attached point-by-point document.
  
  Citation: https://doi.org/10.5194/egusphere-2026-439-AC2

Justine Berg, Pascal Horton, Martina Kauzlaric, Alexandra von der Esch, and Bettina Schaefli

Supplement

https://doi.org/10.5194/egusphere-2026-439-supplement

Justine Berg, Pascal Horton, Martina Kauzlaric, Alexandra von der Esch, and Bettina Schaefli

Viewed

Total article views: 2,227 (including HTML, PDF, and XML)

HTML	PDF	XML	Total	Supplement	BibTeX	EndNote
1,268	848	111	2,227	226	67	117

HTML: 1,268
PDF: 848
XML: 111
Total: 2,227
Supplement: 226
BibTeX: 67
EndNote: 117

Views and downloads (calculated since 25 Feb 2026)

Month	HTML	PDF	XML	Total
Feb 2026	404	188	12	604
Mar 2026	615	421	76	1,112
Apr 2026	88	95	1	184
May 2026	26	65	0	91
Jun 2026	17	12	4	33
Jul 2026	25	9	3	37
Aug 2026	15	12	0	27
Sep 2026	54	36	9	99
Oct 2026	24	10	6	40

Cumulative views and downloads (calculated since 25 Feb 2026)

Month	HTML	PDF	XML	Total
Feb 2026	404	188	12	604
Mar 2026	615	421	76	1,112
Apr 2026	88	95	1	184
May 2026	26	65	0	91
Jun 2026	17	12	4	33
Jul 2026	25	9	3	37
Aug 2026	15	12	0	27
Sep 2026	54	36	9	99
Oct 2026	24	10	6	40

Viewed (geographical distribution)

Total article views: 2,232 (including HTML, PDF, and XML) Thereof 2,232 with geography defined and 0 with unknown origin.

Country	#	Views	%

Latest update: 11 Oct 2026

Short summary

Snow and glacier melt are essential for water supply in mountain regions, but it is challenging to determine how much each source contributes to streamflow. To improve these estimates, we simulated 14 headwater catchments in Switzerland using a hydrological model modified to include snow redistribution and integrating glacier melt from a glacier evolution model. This allowed more accurate estimates of snow and glacier melt contributions, supporting future predictions of water availability.


Total:	0
HTML:	0
PDF:	0
XML:	0