<figure><img src="image_1.png" /></figure>

HAL
open science

# Brief communication: Non-linear sensitivity of glacier mass balance to climate attested by temperature-index models

Christian Vincent, Emmanuel Thibert

► To cite this version:

Christian Vincent, Emmanuel Thibert. Brief communication: Non-linear sensitivity of glacier mass balance to climate attested by temperature-index models. The Cryosphere, 2023, 17 (5), pp.1989-1995. 10.5194/tc-17-1989-2023. hal-04305821

HAL Id: hal-04305821

https://hal.inrae.fr/hal-04305821

Submitted on 25 Nov 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

# Brief communication: Non-linear sensitivity of glacier mass balance to climate attested by temperature-index models

Christian Vincent and Emmanuel Thibert

Université Grenoble Alpes, CNRS, IRD, Grenoble-INP, INRAE, Institut des Géosciences de l'Environnement (IGE, UMR 5001), 38000 Grenoble, France

**Correspondence:** Christian Vincent (christian.vincent@univ-grenoble-alpes.fr) and Emmanuel Thibert (emmanuel.thibert@inrae.fr)

Received: 24 October 2022 – Discussion started: 28 October 2022

Revised: 28 February 2023 – Accepted: 2 April 2023 – Published: 12 May 2023

**Abstract.** Temperature-index models have been widely used for glacier-mass projections spanning the 21st century. The ability of temperature-index models to capture non-linear responses of glacier surface mass balance (SMB) to high deviations in air temperature and solid precipitation was recently discussed in the context of mass-balance simulations employing advanced machine-learning techniques. Here, we performed numerical experiments with a classic temperature-index model and confirmed that such models are capable of detecting non-linear responses of glacier SMB to temperature and precipitation changes. Non-linearities derive from the change in the degree-day factor over the ablation season and from the lengthening of the ablation season.

## 1 Introduction

Glacier surface-mass-balance (SMB) projections in response to climate change up to the end of the 21st century can be analysed via physical approaches using energy-balance calculations and empirical approaches linking simple meteorological variables to SMB such as temperature-index models. Most glacier-mass projections in response to climate change in large-scale studies spanning the 21st century have been based on temperature-index models (Huss and Hock, 2015; Fox-Kemper et al., 2021), given the lack of available or reliable information on detailed future meteorological variables (Réveillet et al., 2018). The deep artificial-neural-network (ANN) approach is a promising new empirical approach to simulate SMB in the future (Bolibar et al., 2020). A neural network is a collection of interconnected simple processing

elements called neurons. These processing elements are assigned coefficients or weights, which constitute the neural-network structure. Each weight is generated by the training process for the ANN (Agatonovic-Kustrin and Beresford, 2000).

Recently, Bolibar et al. (2022) analysed the sensitivity of glacier SMB to future climate change using a deep ANN. They write that their deep-learning approach captures non-linear responses of glacier SMB to high deviations in air temperature and solid precipitation, improving the representation of extreme SMBs. Bolibar et al. (2022) argue that temperature-index models, widely used to simulate the large-scale evolution of glaciers, can be suitable for steep mountain glaciers but may be less suitable for some scenarios and flatter glaciers and ice caps due to linear sensitivities in such mass-balance models. In this paper we perform numerical experiments with a classic and simple temperature-index model. Our unique purpose is to demonstrate that temperature-index models are able to capture non-linear responses of glacier mass balance (MB) to high deviations in air temperature and solid precipitation.

## 2 Data

For our numerical experiments, we selected two very different glaciers in the French Alps. The first, the Argentière Glacier, is located in the Mont Blanc range (45°55' N, 6°57' E). Its surface area was approximately 10.9 km² in 2018. The glacier extends from an altitude of approximately 3400 m a.s.l. at the upper bergschrund down to 1600 m a.s.l. at the snout. It faces north-west, except for a large part of the
accumulation area (south-west-facing tributaries). The second, the Sarennes Glacier, is a small south-facing glacier (0.51 km²) with a limited altitude range between 2820 and 3160 m (mean values over the period used for the present study) located in the Grande Rousses range (45°07′N, 6°07′E). The field SMB observations of the Argentière and Sarennes glaciers come from the French glacier-monitoring programme GLACIOCLIM (Les GLACIers, un Observatoire du CLIMat; https://glacioclim.osug.fr/, last access: 25 April 2023). Annual SMBs were monitored in the ablation area of the Argentière Glacier between 1975 and 1993 using 20 to 30 ablation stakes. Since 1993, systematic winter and summer mass-balance measurements (May and September respectively) have been carried out over the entire surface of the glacier. Approximately 40 sites were selected at various elevations representative of the whole surface. Moreover, geodetic mass balances have been calculated using digital elevation models on the basis of an old map from 1905 and photogrammetric measurements carried out in 1949, 1980, 1993, 1998, 2003, 2008, and 2019 (Vincent et al., 2009). Since 1949, systematic winter and summer mass-balance measurements have been carried out on the Sarennes Glacier, from which annual balances are calculated (Thibert et al., 2013).

We used the atmospheric temperature and precipitation data from the SAFRAN (Système d'Analyse Fournissant des Renseignements Adaptés à la Nivologie analysis system for the provision of information for snow research) reanalysis process that are available from 1958 to date (Durand et al., 2009; Verfaillie et al., 2018). SAFRAN disaggregates large-scale meteorological analyses and observations in the French Alps. The analyses provide hourly meteorological data as a function of seven slope exposures (N, S, E, W, SE, SW, and flat) and altitude (at 300 m intervals up to 3600 m a.s.l.) that differ for each mountain range (e.g. the Mont Blanc, Vanoise, and Grandes Rousses ranges).

### 3 Method

We ran numerical experiments with a classic simple temperature-index model (Hock, 2003; Reveillet et al., 2017) and using SAFRAN reanalysis data (Durand et al., 2009; Verfaillie et al., 2018). These numerical experiments were run on the two very different French glaciers, Argentière and Sarennes, observed over several decades (Thibert et al., 2013; Vincent et al., 2009). The SMB model was run for each day using the equation

$$SMB = DDF_{\text{snow/ice}} \cdot T + k \cdot P,$$

where T is the difference between the mean daily air temperature and the melting point, DDF_snow/ice is the degree-day factor for snow and ice, DDF = 0 if T < 0 °C, P is the precipitation (m w.e.), k is a ratio between snow accumulation and precipitation, and k = 0 if T > 0 °C.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Glacier-wide mass balance of the Argentière Glacier (1990–2015) and the Sarennes Glacier (1949–2015). Observations and simulations from the simple degree-day model used in our experiments.</figcaption></figure>

The degree-day factors for snow and ice were 0.0035 and 0.0055 m w.e. K⁻¹ d⁻¹ for the Argentière Glacier (Reveillet et al., 2017) and 0.0041 and 0.0068 m w.e. K⁻¹ d⁻¹ for the Sarennes Glacier (Thibert et al., 2013). The point mass balances were calculated for each elevation for the Argentière and Sarennes glaciers. In addition, we calculated the glacier-wide mass balance of the Argentière Glacier using the point mass balances for the elevation range and the geodetic mass balances (Vincent et al., 2009). Parameter k depends on the site elevation in accounting for the precipitation gradient and is determined from the winter-balance measurements and precipitation data.

Other enhanced temperature-index models including potential direct solar radiation could be used for our study, but here the purpose is to show that responses in SMB are not linear to temperature or precipitation changes even when using a simple degree-day model.

### 4 Results

The reconstruction of the glacier-wide MBs of these glaciers from our simple temperature-index model shows good agreement with data (Fig. 1). Using these reconstructed MBs, we calculated the SMB sensitivities to temperature and winter precipitation at 2750 and 3100 m on the Argentière and Sarennes glaciers respectively (Fig. 2). These altitudes were selected because they correspond to the approximate centre of the glaciers. For each day of each series, we calculated an annual SMB anomaly by adding a temperature anomaly or a precipitation anomaly. The anomaly was generated as a shift (increment/decrement) in the mean of the distribution of the original data in temperatures and winter balances. The distribution around the means was unchanged (same year-to-year variability as found in the original data).
<figure><img src="image_6.png" /></figure>

**Figure 2.** Response of mass balance to climate forcing using a temperature-index model (green line) at 2750 and 3100 m on the Argentière (a, c) and Sarennes (b, d) glaciers respectively. The red dashed lines are the best linear fit. Note that, in such graphs, the sensitivity of the mass balance to temperature and winter accumulation changes is the slope of the curves.

<figure><img src="image_7.png" /></figure>

**Figure 3.** Response of annual mass balance to air temperature (left panel) and to winter accumulation (right panel) using a temperature-index model (green line) on the Argentière Glacier. The red dashed lines are the best linear fit forced through the origin.
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Positive degree-day model running on synthetic data (response to air temperature). Evolution of air temperatures (a), ablation rates (b), and mass balance (c) over the year, according to different temperature scenarios, calculated at 2800 m. Note the jump in ablation rates when ablation shifts from snow to ice. This occurs earlier with temperature forcing. Note also the lengthening of the ablation season with a rise in temperature.</figcaption></figure>

We report the results in Fig. 2 to mirror Fig. 3 of Bolibar et al. (2022) and make the comparison easier. We also ran these numerical experiments at different altitudes and over the entire glacier surface of the Argentière Glacier (Fig. 3).

From our experiments, we found first that the response of SMB to temperature, using a temperature-index model, is not linear. As expected, the sensitivity of annual SMB (i.e. the slope of the green curves in the graphs of Fig. 2) increases with the positive-degree-day (PDD) anomaly. To explain the physical processes involved in non-linearity, we again used our PDD model but with synthetic data for atmospheric temperature changes over a year (Fig. 4a). The reference scenario (unforced temperature and winter-balance reference conditions) of synthetic data is typical for a location in the upper ablation area of an Alpine glacier (cumulative PDD of 800 degree-days from early May to early October; 1700 mm of winter balance). We use increments of ±1 K (−5 K; +5 K) to analyse the response of SMB. PDD factors for snow and ice come from Thibert et al. (2013). As shown in Fig. 4, the non-linearity with respect to temperature forcing (the spread between SMB plots in Fig. 4c) comes from (i) the lengthening

of the ablation season (Fig. 4a) and (ii) the earlier disappearance of the winter snow cover which increases the ablation rate due to the change in the degree-day factor from snow to ice (Fig. 4b).

Concerning the winter balance, runs of our PDD model on synthetic data under different conditions of winter balance (Fig. 5) used a reference scenario of 1700 mm of winter balance changed by increments of ±300 mm in precipitation. We found a non-linear response of SMBs to winter precipitation with our PDD model. For instance, with winter accumulation decreased by −1500 mm, ice ablation starts very early (by the end of May), and the annual MB is −5.55 m w.e. a⁻¹ in October. With winter accumulation increased by +1500 mm, ice ablation starts in mid-September and the annual MB is −0.21 m w.e. a⁻¹ in October. This asymmetry clearly shows that the response to winter accumulation is not linear. Results show that the increase in sensitivity can be physically explained by the earlier disappearance of the winter snow cover. The earlier and abrupt increase in the ablation rate under lower conditions of winter balance (Fig. 5a) results in non-linearity attested by the spread be-
<figure><img src="image_9.png" /><figcaption>**Figure 5.** Positive degree-day model running on synthetic data (response to winter balance). Change in ablation rates (a) and mass balance (b) over the year according to different winter-balance scenarios calculated at 2800 m. Note the jump in ablation rates when ablation shifts from snow to ice. This occurs earlier under lower winter-balance conditions. Note that the duration of the ablation season is unchanged under variable winter-balance conditions.</figcaption></figure>

tween SMB plots in Fig. 5b. We detect sensitivity to winter accumulation, contrary to the Bolibar et al. (2022) findings using their ANN (Figs. 2 and 3). Indeed, MB sensitivity increases with low winter-accumulation anomalies using our model but decreases in the deep-learning model of Bolibar et al. (2022). Our results are consistent with direct in situ observations (Six and Vincent, 2014) and are also consistent with the results reported by Réveillet et al. (2018) from observations and energy-balance modelling. The opposite results obtained from the deep-learning model may be due to an issue in the calibration of the model.

Summing up, the ability of PDD models to provide non-linear sensitivity to air temperature and solid precipitation is due to the different ablation rates and the associated change in the degree-day factor that can be involved depending on snow or ice conditions at the glacier surface. An additional non-linearity to temperature forcing is caused by changes in the ablation duration.

Another question arises in the discussion section of Bolibar et al. (2022) concerning the comparison between their results and those from other studies. According to this paper, all glacier models in the Glacier Model Intercomparison Project (GlacierMIP) (Hock et al., 2019) rely on

SMB models with linear relationships between PDDs, melt, and precipitation. The authors argue that these PDD models present behaviour very similar to the linear-build statistical Least Absolute Shrinkage and Selection Operator (LASSO) model. However, most of the temperature-index models used in GlacierMIP include two degree-day factors. Consequently, they cannot provide a linear response to climate forcing as shown above. In the Bolibar et al. (2022) paper, the MB anomalies in response to climate forcing were obtained using a linear LASSO SMB model, which is based on a regularized multi-linear regression, although the authors also used the Global Glacier Evolution Model (GloGEMflow) model in their paper (see their discussion section), which is a temperature-index model widely used for glacier projections (Zekollari et al., 2019).

## 5 Conclusions

From numerical experiments with a classic and simple temperature-index model, our results highlight that temperature-index models are able to capture non-linear responses of glacier mass balance (MB) to high deviations in
air temperature and solid precipitation, unlike the Bolibar et al. (2022) study.

Bolibar et al. (2022) compared the response of SMB to climate forcing (air temperature, winter and summer snowfall) using a deep-learning approach and a LASSO model. From this comparison, they conclude that deep learning provides a non-linear response, contrary to the LASSO model. The conclusions of Bolibar et al. (2022) may be due to the use of a linear LASSO SMB model instead of a temperature-index model. We would suggest testing the capability of an ANN to capture non-linearity by comparing its results with that of the GloGEM PDD model that they used in their paper.

Regarding specifically SMB changes due to solid pre-cipitation, the deep-learning model used by Bolibar et al. (2022) foresees decreasing sensitivity under low winter-accumulation conditions. We point out that this result directly contradicts PDD model outcomes. We explain in physical terms why a PDD model projects higher sensitivity to low winter accumulation but do not yet understand why the approach of Bolibar et al. (2022) does not.

Given that detailed meteorological variables are highly unpredictable in the future, most glacier-mass projections in response to climate change in large-scale studies spanning the 21st century are still today based on temperature-index models with simple temperature and precipitation variables. It follows that the questions raised here relative to the non-linear responses of surface SMB to meteorological variables are crucial.

***Data availability.*** This commentary does not include original data. All data referred to in the text have been published elsewhere. Field data are accessible through the project website at https://glacioclim.osug.fr, last access: 11 May 2023 (Sarennes: https://nextcloud.osug.fr/index.php/s/ToNQAYbQdZyqEkK, Six, 2023a; Argentière: https://nextcloud.osug.fr/index.php/s/65R2kJP8DMkFGzC, Six, 2023b).

Results from the PDD simulations on synthetic data are accessible from the Open Data Repository: https://doi.org/10.5281/zenodo.7603415 (Vincent and Emmanuel, 2023).

***Author contributions.*** ET and CV ran the numerical modelling calculations and produced the analysis. CV supervised the study and wrote the paper. Both authors contributed to discussion of the results.

***Competing interests.*** The contact author has declared that neither of the authors has any competing interests.

***Disclaimer.*** Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

***Acknowledgements.*** This study was funded by the Observatoire des Sciences de l'Univers de Grenoble (OSUG) and Institut des Sciences de l'Univers (INSU-CNRS) in the framework of the French GLACIOCLIM (Les GLACIers, un Observatoire du CLIMat) programme. We thank all those who conducted the field measurements. We are grateful to Cary Bartsch for reviewing the English.

***Financial support.*** This study was funded by Observatoire des Sciences de l'Univers de Grenoble (OSUG) and Institut des Sciences de l'Univers (INSU) in the framework of the French "GLACIOCLIM (Les GLACIers comme Observatoire du CLIMat)".

***Review statement.*** This paper was edited by Daniel Farinotti and reviewed by Jordi Bolibar and three anonymous referees.

## References

Agatonovic-Kustrin, S. and Beresford, R.: Basic concepts of artificial neural network (ANN) modeling and its application in pharmaceutical research, J. Pharmaceut. Biomed., 22, 5, https://doi.org/10.1016/s0731-7085(99)00272-1, 2000.

Bolibar, J., Rabatel, A., Gouttevin, I., Galiez, C., Condom, T., and Sauquet, E.: Deep learning applied to glacier evolution modelling, The Cryosphere, 14, 565–584, https://doi.org/10.5194/tc-14-565-2020, 2020.

Bolibar, J., Rabatel, A., Gouttevin, I., Zekollari, H., and Galiez, C.: Nonlinear sensitivity of glacier mass balance to future climate change unveiled by deep learning, Nat. Commun., 13, 409, https://doi.org/10.1038/s41467-022-28033-0, 2022.

Durand, Y., Laternser, M., Giraud, G., Etchevers, P., Lesafre, B., and Mérindol, L.: Reanalysis of 44 yr of climate in the French Alps (1958–2002): Methodology, model validation, climatology, and trends for air temperature and precipitation, J. Appl. Meteorol. Clim., 48, 429–449, https://doi.org/10.1175/2008JAMC1808.1, 2009.

Fox-Kemper, B., Hewitt, H. T., Xiao, C., Aðalgeirsdóttir, G., Drijfhout, S. S., Edwards, T. L., Golledge, N. R., Hemer, M., Kopp, R. E., Krinner, G., Mix, A., Notz, D., Nowicki, S., Nurhati, I. S., Ruiz, L., Sallée, J.-B., Slangen, A. B. A., and Yu, Y.: Ocean, Cryosphere and Sea Level Change, in: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Masson-Delmotte, V., Zhai, P., Pirotani, A., Connors, S. L., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M. I., Huang, M., Leitzell, K., Lonnoy, E., Matthews, J. B. R., Maycock, T. K., Waterfield, T., Yelekçi, O., Yu, R., and Zhou, B., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 1211–1362, https://doi.org/10.1017/9781009157896.011, 2021.

Hock, R.: Temperature index melt modelling in mountain areas, J. Hydrol., 282, 104–115, https://doi.org/10.1016/S0022-1694(03)00257-9, 2003.

Hock, R., Bliss, A., Marzeion, B., Giesen, R. H., Hirabayashi, Y., Huss, M., Radic, V., and Slangen, A. B. A.: A model intercomparison of global-scale glacier mass-
balance models and projections, J. Glaciol., 65, 453-467,
https://doi.org/10.1017/jog.2019.22, 2019.

Huss, M. and Hock, R.: New model for global glacier
change and sea-level rise, Front. Earth Sc., 3,
https://doi.org/10.3389/feart.2015.00054, 2015.

Reveillet, M., Vincent, C., Six, D., and Rabatel, A.: Which empir-
ical model is best suited to simulate glacier mass balances?, J.
Glaciol., 63, 39–54, https://doi.org/10.1017/jog.2016.110, 2017.

Réveillet, M., Six, D., Vincent, C., Rabatel, A., Dumont, M., Lafaysse, M., Morin, S., Vionnet, V., and Litt, M.: Relative performance of empirical and physical models in assessing the seasonal and annual glacier surface mass balance of Saint-Sorlin Glacier (French Alps), The Cryosphere, 12, 1367–1386, https://doi.org/10.5194/tc-12-1367-2018, 2018. Six, D.: Mass balance of Sarennes glacier (French Alps), https://nextcloud.osug.fr/index.php/s/ToNQAYbQdZyqEkK, last access: 9 May 2023a. Six, D.: Mass balance of Argentiere glacier (French Alps), https://nextcloud.osug.fr/index.php/s/65R2kJP8DMkFGzC, last access: 9 May 2023b.

Six, D. and Vincent, C.: Sensitivity of mass balance and
equilibrium-line altitude to climate change in the French Alps, J.
Glaciol., 60, 867–878, https://doi.org/10.3189/2014JoG14J014,
2014.

Thibert, E., Eckert, N., and Vincent, C.: Climatic drivers of seasonal glacier mass balances: an analysis of 6 decades at Glacier de Sarennes (French Alps), The Cryosphere, 7, 47–66, https://doi.org/10.5194/tc-7-47-2013, 2013.

Verfaillie, D., Lafaysse, M., Déqué, M., Eckert, N., Lejeune, Y., and Morin, S.: Multi-component ensembles of future meteorological and natural snow conditions for 1500 m altitude in the Chartreuse mountain range, Northern French Alps, The Cryosphere, 12, 1249–1271, https://doi.org/10.5194/tc-12-1249-2018, 2018.

Vincent, C. and Thibert, E.: Nonlinear sensitivity of
glacier-mass balance to climate attested by temperature-
index models; synthetic data, Zenodo [data set],
https://doi.org/10.5281/zenodo.7603416, 2023.

Vincent, C., Soruco, A., Six, D., and Le Meur, E.: Glacier thickening and decay analysis from 50 years of glaciological observations performed on Glacier d'Argentière, Mont Blanc area, France, Ann. Glaciol., 50, 73-79, https://doi.org/10.3189/172756409787769500, 2009.

Zekollari, H., Huss, M., and Farinotti, D.: Modelling the future evolution of glaciers in the European Alps under the EURO-CORDEX RCM ensemble, The Cryosphere, 13, 1125–1146, https://doi.org/10.5194/tc-13-1125-2019, 2019.