the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 22 Jan 2026
Sensitivity of Andean Glaciers to ice-flow parameters in the Parallel Ice Sheet Model
Abstract. Mountain glaciers are losing mass rapidly due to anthropogenic climate change. Projections of glacier evolution across the Andes under different warming scenarios have primarily been as part of global scale modelling frameworks, rather than dedicated, regionally optimised, simulations. These global-scale models use simplifications of ice flow physics that may be unsuitable for steep topography, such as that which occurs at mountain valley glaciers. More complex models are available, but with that complexity comes further sources of uncertainty. Here, we assess the sensitivity of the Parallel Ice Sheet Model to ice-flow parameters influencing the ice rheology and subglacial sliding characteristics. We find that the resistance of subglacial material has the most impact on modelled ice outputs (e.g., ice volume), followed by the exponent which relates basal shear stress to sliding, and the threshold velocity at which sliding occurs. The ice-flow rheology enhancement factors, the rate of subglacial water decay, and the maximum water thickness within a presumed subglacial drainage network, can either cause minor variations, or no effect at all, on ice outputs. Our study informs what parameters can potentially be negated in future parameter ensemble tests and provides direction on where further investigation is needed.
- Preprint
(3033 KB) - Metadata XML
-
Supplement
(19625 KB) - BibTeX
- EndNote
Status: final response (author comments only)
- RC1: 'Comment on egusphere-2025-6169', Cristina I. Balaban, 19 Feb 2026
-
RC2: 'Comment on egusphere-2025-6169', Adem Candas, 12 May 2026
General comments
The manuscript presents a parameter analysis for users of the Parallel Ice Sheet Model (PISM). Parameter selection and tuning in PISM can be challenging, particularly for researchers applying the model to complex glacier settings. Therefore, I would like to emphasize that a systematic parameter analysis of this kind is highly useful for the community. This study will be beneficial for researchers using PISM for glacier modelling. A manuscript that focuses specifically on parameter analysis is especially useful for researchers working on mountain glaciers. As the authors also point out, PISM has been more commonly used for ice-sheet applications, while the number of mountain-glacier applications remains relatively limited. When mountain glaciers are considered, additional difficulties arise due to steep topography. While reading the manuscript, I was pleased to find answers to some questions that I had also considered in my own work.
Together with the Supplementary Information, this manuscript has the potential to serve as a useful reference for future researchers working on glacier modelling. I offer the following comments with the intention of improving the clarity of the manuscript and making it more useful for other users. My recommendation is that the manuscript would be suitable for publication after these revisions are addressed.
Specific comments
I think that a short explanation summarizing the main outline of the study at the end of the Introduction would be helpful for the reader. This is, of course, a matter of author preference.
The selected study area is highly suitable for this type of parameter analysis. I think that the topographic characteristics of all five study areas should also be provided. For example, it would be useful to see the minimum, maximum, and mean elevation values for each area. It would also be helpful to better understand the slopes at the valley scale. I think Table 2 may be an appropriate place for this information.
The parameters described in Sections 3.1.1–3.1.3 may need to be explained in more detail in order to reach a broader readership. For example, it would be useful to indicate how an increase or decrease in each parameter is expected to affect glacier thickness, volume, area, and velocity.
I suggest that the authors reconsider the flowchart in Figure 2. Its general structure could be revised, and the workflow should be made more intuitive. For example, Table 1 could be included in the first step.
Technical corrections
Some additional technical suggestions are listed below.
L18 – “The ice-flow rheology enhancement factors, the rate of subglacial water decay, and the maximum water thickness within a presumed subglacial drainage network, can either cause minor variations, or no effect at all, on ice outputs.” It may be useful to clarify that this conclusion applies specifically to mountain glaciers or to the model configurations tested in this study.
L53 – The parameters on which the study focuses could be added here.
L109 – In the description column of Table 1, the direction of the parameter effects could be indicated. For example, the sentence “ESIA and ESSA controls how easily the ice deforms” is useful, but it would be clearer if the authors specified whether increasing or decreasing these parameters leads to this effect.
L121 – Several previous studies used E values of 1–6 and 0–1.5. Could the authors briefly explain why a much wider range of 0.2–20 was selected in this study?
L135 – “till friction angle till friction angle” appears to be repeated.
L149 – I recommend making this explanation slightly clearer for a general readership.
L151 – “The till water decay rate is varied between 0.1 and 12 mm/a, while the maximum till water thickness is varied between 0.1 and 10 m.” Why were these specific values selected?
L175 – “Within the PDD scheme, there is stochastic ‘white noise’ to simulate additional undetermined daily variability, as well as a daily temperature standard deviation that is set by default at 5℃ (Winkelmann et al., 2011). These can cause minor fluctuations in the climate, and thereby in the ice extent even under steady state conditions.” This is actually a debatable issue. There are also studies by Seguinot on this topic, for example doi.org/10.3189/2013JoG13J081. Since this is not the direct focus of the present study, keeping the default setting seems reasonable. However, a short clarification or acknowledgement of this issue may be useful.
L223 – I am not raising this point as a criticism of the approach “bed elevation = topography – present-day thickness”. If present-day thickness values are available, the models could also have been initialized from present-day glaciers. In fact, if present-day glacier outlines are known, this could be a useful calibration constraint. I assume that this was not preferred because the present study is designed as a parameter analysis. Could the authors add a short explanation?
L232 – “we applied a lapse-rate correction of 6.5°C/km”. Was this applied to the entire terrain, or only to the mountain peaks? Could this introduce additional uncertainty?
L234 – “Ultimately, we are not concerned about the size and shape of the glaciers produced.” I could not fully understand this sentence. I suggest rephrasing it.
L269 – “Ice peripheral differences in ice thickness arise from internal variability in the PDD model, despite a temperature standard deviation of 5℃.” This seems to be a rather strong statement. If the authors are confident about this interpretation, it may need to be discussed in the text.
L313 – In Table 4 and the following similar tables, I think the “Average” values may not have a clear statistical meaning because the domains and environmental conditions are different.
L347 – “inferred to be related to the domain geometry or resolution effects within this region.” This seems somewhat speculative and needs further support. Is the resolution not the same in all regions?
L357 – “However, C has been shown to have an increased influence over ice sheets settings.” The phrase “increased influence” is not fully clear. Which parameter or output is affected, and in which direction?
L366 – Thank you for including this information.
L367 – “were saw” should be corrected.
L368 – “…between –40.5% to +23.4% across the domains, respectively” Which domains do these values refer to?
L375 – “Increased values of the ϕ parameter see a reduction in basal velocities.” However, there also seem to be regions where velocities increase. This statement may need to be corrected or clarified.
L406 – These are good results. Thank you!
L433 – “As can be seen in Fig. 11, with decreased 𝑈threshold values overall mean velocities decreased,” Does decreasing the threshold value reduce velocities because it facilitates basal motion? Could the authors clarify the mechanism?
L503 – I think it would be useful to repeat the parameters used here.
L505 – “others” is not clear. Which parameters or studies are being referred to?
L526 – Should this be “Acknowledgements”? Or is this referring to a funded project cited in the text? Please clarify.
Citation: https://doi.org/10.5194/egusphere-2025-6169-RC2
Viewed
| HTML | XML | Total | Supplement | BibTeX | EndNote | |
|---|---|---|---|---|---|---|
| 1,042 | 575 | 123 | 1,740 | 364 | 125 | 249 |
- HTML: 1,042
- PDF: 575
- XML: 123
- Total: 1,740
- Supplement: 364
- BibTeX: 125
- EndNote: 249
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
General comments
The present manuscript aims to test the impact of various ice-flow parameters of the Parallel Ice Sheet Model (PISM) over numerically modelled glacier outputs in the Andes. This study is motivated by a research gap on the appropriate parameterization of glaciological parameters in PISM over mountain glacier settings, the simplistic nature of global-scale models in replicating accurate ice flow, and the importance of the study area for water resources under future climate change projections. The manuscript addresses the aim through two modelling experiment ensembles, iterating enhancement factors, subglacial properties- and basal-sliding-related PISM parameters through default, minimal and maximal values over five Andean hydrological catchments, assessing their effect on ice volume and area. The study finds that some parameters, such as enhancement factors (E), the rate of subglacial water decay (C), and maximum water thickness (Tm) have little-to-no impact on the modelled outputs, whereas some of the subglacial properties-related ones, namely the resistance of subglacial material (ϕ), the exponent relating basal shear stress to sliding (q), and the sliding velocity threshold (𝑈𝑡ℎ𝑟𝑒𝑠ℎ𝑜𝑙𝑑) have a significant or higher impact. The authors conclude that such findings can help in the selection of the appropriate ice-flow parameters in future work, improving both research design and conservation of computational resources.
Overall, this manuscript represents a significant contribution to our understanding of how several PISM ice flow parameters influence ice volume and area outputs, which have not been intensively experimented with in mountain glaciation settings before. Thus, I foresee this work being a future reference point for many PISM users, who often resort to default parameter values or guesswork in their choices when modelling the (palaeo)climate conditions favourable for empirically constrained glacier extents. Without accounting for uncertainty/sensitivity of model parameters, results are less robust and convincing, which this work emphasizes through the variability exhibited in several ice volume outputs. Additionally, as many users encounter inefficient use or estimation of required computing resources in High Performance Computing access applications, this work will help streamline the workflow and prioritise the most relevant sensitivity tests. Furthermore, I commend the authors for packaging this work into a concise, well-structured and -illustrated manuscript that is mostly easy to follow. The Zenodo repository files further provide useful PISM scripts that can be checked by interested readers and model users. The Supplementary Information includes clear tables and mapped model outputs that support the manuscript interpretations and conclusions.
Specific comments
My main comment pertains to embedding more literature in the explanations of the chosen parameters and their behaviour. While this is done well at times, the manuscript could engage with more empirical and/or modelling studies that would support their choices and arguments, or when it does not exist, to explain how such choices were made (see in-line comments). Furthermore, I would encourage the authors to emphasize the potential impact of size and climatic zones of the model domains on the ice-flow parameter sensitivity tests – while the former is acknowledged throughout the results and interpretations, there is no unifying conclusion emerging out of it. In the case of the latter, the climatic zones are described earlier in the manuscript (lines 72-79) but their potential implications (or lack thereof?) are not addressed. Therefore, I think a small subsection/paragraph(s) on these aspects could enrich the manuscript.
Technical corrections
Below are some further in-line and Figure/Table comments that include minor suggestions, changes and considerations for the authors to include in the manuscript or consider in their upcoming or future work.
Lines 54-55 – perhaps it is better to be specific about the exact kind of statistical analyses undertaken (these are revealed later, but might be good to explicitly state them here)
Lines 65-67 – this long phrase could benefit from clarity – perhaps write it as two sentences
Figure 1 – clear, nicely labelled maps. One minor suggestion about the DEM colour scheme– at first sight, it is difficult to tell which is higher/lower ground (in the case of the zoomed-in maps). In my opinion, this could be improved by underlaying them with a hillshade model or using a more traditional topographic colour scheme for the DEM (brown-high – green-low). Otherwise, including the elevation label for each zoomed-in map might also work well.
Table 1 – Table is organized very clearly. It might also be useful to include a column with references for each parameter. I am also wondering where the minimum and maximum value for the parameters have been determined from – this could go either in-text or in the table.
Lines 123-124 – it might be useful to explain the variation between the min and max value, both in terms of values and intervals, justifying the rationale
Line 126 – is it called T component because of Till? For the sliding (S) and enhancement (E) parameters, this seems logical – perhaps a mention either in the heading or text would help.
Lines 128-133 – this paragraph is explained very well, but I feel some references might be needed here, e.g., on what basis do you know about the thick sediment in the glacier forefields? Are the glaciological explanations supported by other studies?
Lines 184-185 – could perhaps explain how the choice of the number of ice and bedrock layers was determined – was it by trial or error or suggested by previous work? Are there any effects of this choice on the model outputs and what would these potentially look like?
Line 186 – what is not clear at this point for the reader is whether the mountain glaciers were grown from no ice conditions – perhaps this can be explicitly stated.
Lines 200-205 – this paragraph seems vague – what do you mean by ‘commonly used and extended ranges’? Although the ensemble design becomes a little clearer as the reader progresses through the manuscript, a table exemplifying the combination could be useful – there are some examples in the Supplementary Information that could be cited or brought forward as a table here.
Lines 224-225 - ‘We used Davies (2013) which uses…’ - try to avoid repetition
Lines 230-232 - very good acknowledgement of the WorldClim dataset limitation regarding underestimation of temperature in mountain peaks. Perhaps in future work, the authors could consider the CHELSA dataset (Karger et al., 2017; 2023). It is a downscaling of ERA-interim climatic reanalysis. Although proven to resolve temperature similar to other climate datasets, it yielded more accurate results for precipitation (Karger et al., 2017).
Lines 244-245 – Although somewhat understandable, this sentence could benefit from some clarity/simplification
Line 251- ‘varying’ and ‘varied’ – consider avoiding repetition in the same sentence
Lines 335-337 – this phrase is slightly long and convoluted – it might need simplifying.
Line 362 – remove second ‘.’ After ‘difficult’
Lines 363-365 – good explanation as to why some of the parameters do not have a sufficient effect on the modelled ice outputs. Regarding the PISM hydrology, it might be good to briefly add something here about the default hydrological model and its alternatives (I am aware experimenting with different hydrological models is a suggestion for further work later in the manuscript)
Lines 378-380 – excellent explanation the different ϕ values but might be useful to include some references here
Line 392 – ‘simulations’ instead of typo ‘simulation s’
Line 393 – do you mean ‘centered’ instead of ‘cantered’?
Lines 393 – 395 – good interpretation of the strong influence of ϕ on ice volume outputs – I am wondering whether highlighting this in envelopes on Figure 7 could make them more obvious to the reader? Furthermore, you use catchment numbers in the text (e.g., #4, #5), but these are only labelled with letters on Figure 7 (A-E). Consider adjusting either on the Figure on in the text.
Lines 397-399 – keep font consistent in the final submission
Lines 481-482 – typo ‘parameters mention’ instead of ‘mentioned’. However, it might be better to recap the names of the specific parameters instead
Lines 497-500- good to see appropriate recommendations for sensitivity tests. I would just add that not only is the PISM PDD model sensitive to climatic parameters, but so is the diurnal energy balance model simple (dEBM-simple). This model was argued to be an improved alternative to the PDD model, by accounting for the melt-albedo feedback, without significantly increasing computational time (Zeitz et al., 2021; Garbe et al., 2023). However, to the best of my knowledge, this model has only been applied in ice sheet settings. Therefore, it could represent a research gap for mountain glaciation applications, and it might be worth mentioning as a recommendation for future work.
Line 535 – zenodo link is not easily clickable/searchable – please adjust the link (this one worked: https://zenodo.org/records/17878115 )
Suggested references:
Garbe, J., Zeitz, M., Krebs-Kanzow, U., & Winkelmann, R. (2023). The evolution of future Antarctic surface melt using PISM-dEBM-simple. The Cryosphere, 17(11), 4571-4599. https://doi.org/10.5194/tc-17-4571-2023
Karger, D. N., Lange, S., Hari, C., Reyer, C. P., Conrad, O., Zimmermann, N. E., & Frieler, K. (2023). CHELSA-W5E5: daily 1 km meteorological forcing data for climate impact studies. Earth System Science Data, 15(6), 2445-2464. https://doi.org/10.5194/essd-15-2445-2023
Karger, D.N., Conrad, O., Böhner, J., Kawohl, T., Kreft, H., Soria-Auza, R.W., Zimmermann, N.E., Linder, H.P. and Kessler, M., 2017. Climatologies at high resolution for the earth’s land surface areas. Scientific data, 4(1), pp.1-20. https://doi.org/10.1038/sdata.2017.122
Zeitz, M., Reese, R., Beckmann, J., Krebs-Kanzow, U., & Winkelmann, R. (2021). Impact of the melt–albedo feedback on the future evolution of the Greenland Ice Sheet with PISM-dEBM-simple. The Cryosphere, 15(12), 5739-5764.https://doi.org/10.5194/tc-15-5739-2021