Introduction
Debris-covered glaciers are a common glacier type in alpine environments, especially in the high-relief Hindu Kush–Himalaya (HKH) range (e.g. Reference Moribayashi and HiguchiMoribayashi and Higuchi, 1977; Reference Benn, Kirkbride, Owen, Brazier and EvansBenn and others, 2003), where most of the ablation areas are partly or even totally covered with supraglacial debris. Given that the HKH region is the largest ice mass outside the polar regions (Reference Dyurgerov and MeierDyurgerov and Meier, 2005) and the most populated on Earth, it is potentially one of the most critical parts of the world in terms of the social and economic impacts of glacier shrinkage (Reference Barnett, Adam and LettenmaierBarnett and others, 2005; Reference Immerzeel, Van Beeke and BierkensImmerzeel and others, 2010; Reference BolchBolch and others, 2012). Predicting the evolution of HKH glaciers is therefore a key issue as their melting may (1) negatively affect regional water supply in the next decades (Reference Barnett, Adam and LettenmaierBarnett and others, 2005; Reference Immerzeel, Van Beeke and BierkensImmerzeel and others, 2010), (2) significantly contribute to ongoing sea-level rise (Reference Kaser, Cogley, Dyurgerov, Meier and OhmuraKaser and others, 2006) and (3) increase natural hazards linked to glaciers, especially glacial lake outburst floods (e.g. Reference Mool, Bajracharya and JoshiMool and others, 2001). Recent research has not reached a consensus, in terms of the role played by debris at the surface, about glacier mass balance in the HKH and the response of debris-covered glaciers to climate change (Reference Scherler, Bookhagen and StreckerScherler and others, 2011; Reference Kääb, Berthier, Nuth, Gardelle and ArnaudKääb and others, 2012). The role of debris has hitherto been neglected in predictions of future water availability (Reference Immerzeel, Van Beeke and BierkensImmerzeel and others, 2010) or sea-level rise at the global scale.
A supraglacial debris layer significantly influences the melting of the underlying ice and has a specific surface energy balance very different to that of debris-free glaciers. The main physical characteristics of a debris layer (in contrast to snow or ice) are its thermal conductivity, its albedo and the fact that its temperature can rise above 273.15 K (Nicholson and Benn, in press). The debris has a conflicting effect on the underlying ice through a reduced albedo, dominant for thin debris layers, and increased thermal shielding, dominant for thicker debris cover. Reference ØstremØstrem (1959) established empirical relationships between supra-glacial debris thickness and ice-melt rates, and introduced the concept of critical debris thickness, below which ice ablation is enhanced and above which it is reduced, compared to a debris-free surface. This pattern has been confirmed in numerous subsequent studies (e.g. Reference LoomisLoomis, 1970; Reference Mattson, Gardner and YoungMattson and others, 1993; Reference Conway and RasmussenConway and Rasmussen, 2000; Reference Kayastha, Takeuchi, Nakawo and AgetaKayastha and others, 2000; Nicholson and Benn, in press), which further show that the critical thickness varies under the influence of debris lithology and local climate.
In order to estimate the mass balance of debris-covered glaciers, empirical methods, such as the degree-day approach, have been applied to a few glaciers (Reference Kayastha, Takeuchi, Nakawo and AgetaKayastha and others, 2000), but they are too site-specific and too dependent on the conditions prevailing during the measurement period to allow any reliable spatial or temporal extrapolation (e.g. Reference HockHock, 2003; Reference Nicholson and BennNicholson and Benn, 2006). To predict both short-term melt rates in response to meteorological conditions and long-term glacier ablation regimes, which influence the dynamic response of debris-covered glaciers to climate forcing, a physically based ablation model is needed (Reference Nicholson and BennNicholson and Benn, 2006). Consequently, an effort has recently been made by glaciologists to develop physically based energy-balance models to assess how the debris layer affects glacier melt rates (Reference Nicholson and BennNicholson and Benn, 2006; Reference Reid and BrockReid and Brock, 2010; Reference Reid, Carenzo, Pellicciotti and BrockReid and others, 2012). The numerical models introduced by Reference Nicholson and BennNicholson and Benn (2006) and Reference Reid and BrockReid and Brock (2010) were developed only for the ablation season, assuming that the ice/debris interface always remains at the melting point. This assumption is sometimes not true during clear nights where refreezing is efficient, and is obviously not fulfilled outside the ablation season. Contrary to the model of Reference Nicholson and BennNicholson and Benn (2006), which operates at a daily time resolution, assuming a linear temperature profile between the upper and lower surfaces of the debris layer, the model of Reference Reid and BrockReid and Brock (2010) fully solves the heat conduction equation through the debris using a time-step of 1 hour. Working at such a short timescale is necessary to properly simulate the processes responsible for the diurnal melting of the underlying ice, adequately accounting for the competition between the surface energy budget and thermal diffusion through the debris. The model of Reference Reid and BrockReid and Brock (2010) is thus able to reproduce the existence of the critical debris thickness initially described by Reference ØstremØstrem (1959). Here, we present a physically based model able to simulate mass and energy fluxes of a debris-covered glacier all year round, i.e. able to calculate the temperature of the ice/debris interface, to simulate the ice melt and to manage the presence of a transient snowpack on top of the debris. This model is called Crocus-DEB, because it stems from the detailed snowpack model (Reference Brun, Martin, Simon, Gendre and ColéouBrun and others, 1989, Reference Brun, David, Sudul and Brunot1992) which has been used in the past for distributed simulations of glacier mass balance (Reference Gerbaux, Genthon, Etchevers, Vincent and DedieuGerbaux and others, 2005; Reference LejeuneLejeune, 2009; Reference Dumont, Durand, Arnaud and SixDumont and others, 2012). The suffix ‘-DEB’ used to discriminate between the standard and modified versions of Crocus refers to the first syllable of the word ‘debris’, consistent with Reference Reid and BrockReid and Brock (2010). The detailed snowpack model Crocus explicitly evaluates mass and energy exchange between the snowpack and the low-level atmosphere as a function of meteorological conditions. In Crocus-DEB, the debris layer has been inserted as stacks of snow layers, featuring physical properties corresponding to the debris.
In this study, we first give a detailed description of the model, focusing on the adaptations made to include debris layers in the snowpack. Second, we present the results of simulations performed at the instrumented experimental site Col de Porte (Chartreuse range, French Alps, 1325 m a.s.l.; 45°17′ N, 05°45′ E), where we artificially poured a 4 cm and a 15 cm thick debris layer over seasonal snow cover in the middle of winter 2011, to quantify the effect of the debris over the snowmelt, thereby providing appropriate driving and evaluation data to inform and test the Crocus-DEB model. Third, the performance of the model is also evaluated using a comprehensive meteorological and glaciological dataset obtained on a Nepalese debris-covered glacier, Changri Nup Glacier (5360 m a.s.l.; 27°59′ N, 86°47′ E), between November and December 2010. In the final sections we critically discuss the performance of the model and show how it can be used to explore the relationships between the critical debris thickness and the meteorological and environmental conditions of a given glacier.
Material and methods
Model description
In simple terms, the Crocus-DEB model was developed from the multilayer Crocus snowpack model, building on the idea that a debris layer can be viewed as a snow layer featuring different physical properties. Below, we briefly recall the main characteristics of the Crocus snowpack model which are relevant to the description of Crocus-DEB.
Brief overview of the Crocus snowpack model
Crocus is a one-dimensional, multilayer, physically based snowpack model that simulates in detail the time evolution of the physical properties of layers making up the snowpack, solely driven by meteorological conditions at the surface (air temperature, wind speed, relative humidity, incoming long-wave and shortwave radiation, snowfall and precipitation rate). For each numerical snow layer, i, within the snow-pack, the model prognostic variables are temperature, T(i) (K), dry density, ρ(i) (kg m−3), i.e. the density of the ice matrix, thickness, dz(i) (m), liquid water content, w(i) (kg m−3), snow age (days since snowfall) and grain-type properties. The latter are described by four empirical variables (dendricity, sphericity, size and metamorphic state) representing the snow microstructure. More details on the variables used to describe the snow layers in Crocus are given by Reference Bouilloud and MartinBouilloud and Martin (2006), Reference WillemetWillemet (2010) and Reference VionnetVionnet and others (2012). The energy budget of each snow layer, i, is written as follows:
where Cp (i) is the specific heat capacity of ice, depending on temperature, and L f is the latent heat for ice fusion (3.33 × 105 J kg−1). On the right-hand side of the equations, W p and W represent the liquid water input through rainfall or percolation, respectively. S abs (i) represents the fraction of incoming shortwave radiation absorbed by layer i, L net is the net longwave radiation flux, H and LE are the turbulent fluxes for sensible and latent heat, respectively, P is the sensible heat flux due to precipitation, Q c (i) represents the divergence of the conduction flux within layer i and Q g is the basal heat flux. Figure 1 provides an overview of the different processes accounted for by the Crocus snowpack model.
The computation of energy fluxes within and at the boundaries of the snowpack by Crocus has been described extensively in previous publications (e.g. Reference Brun, Martin, Simon, Gendre and ColéouBrun and others, 1989, Reference Brun, David, Sudul and Brunot1992; Reference Martin and LejeuneMartin and Lejeune, 1998; Reference WillemetWillemet, 2010; Reference VionnetVionnet and others, 2012) and is not repeated here.
From Crocus to Crocus-DEB
Debris-covered glaciers are modeled as a vertical stack of ice, debris and a potentially developing snowpack at the surface. In this regard, the basic principles of Crocus, i.e. the time evolution of the physical processes of a snowpack consisting of several numerical layers, remain unchanged. We describe below the handling of debris layers by Crocus-DEB, including the simplifying hypotheses made.
Variables. To differentiate between numerical layers consisting of debris and those of snow and ice, we use an arbitrarily out-of-range value for one of the prognostic variables used by Crocus to handle snow metamorphism (namely, the variable representing either sphericity or grain size). The only true prognostic variable considered for a debris layer is its temperature. In contrast to a snow or ice layer, the temperature of a debris layer can take values >273.15 K. Debris layers are attributed a given thickness, density, thermal conductivity and specific heat capacity. These properties remain constant throughout a given simulation. In this work, all the numerical layers making up the debris cover were attributed the same physical properties, but this is not imposed by the model structure. We further hypothesize that the liquid water content of debris layers is always zero, implying that liquid water percolation through the debris is instantaneous and that phase changes do not occur in the interstitial volume in the debris layers.
Processes. When there is a debris layer at the surface, the computation of the surface energy budget is carried out in a different way to the case where snow is present at the surface. When debris is present at the surface, different values than those for snow are used for the thermal emissivity, ε, roughness length, z 0, and spectral albedo. The spectral absorption coefficient of the debris is considered infinite, i.e. all of the net shortwave radiation budget is affected to the uppermost debris layer. The thermal emissivity of the debris, ε d, is kept constant. The spectral albedo of the debris is fixed throughout a given simulation. The sensible heat flux, H, is computed similarly to snow; however the debris can be attributed a different roughness length, z 0. The latent heat flux, LE, is zero because the debris is considered to be always dry. Liquid water from rainfall is supplied to the uppermost debris layer. The corresponding heat flux is computed similarly to the rainon-snow case (i.e. rainfall temperature is taken equal to air temperature) (Reference Brun, Martin, Simon, Gendre and ColéouBrun and others, 1989, Reference Brun, David, Sudul and Brunot1992; Reference VionnetVionnet and others, 2012). For the case of snowfall on top of the debris (in a similar way to how Crocus handles snowfall on the ground) a snow cover may or may not develop, depending on the thermal state of the debris. When the debris temperature is below the melting point for ice, a snow cover develops on top of the debris. For higher debris temperatures, development of a snow cover depends on the balance between the heat content of the debris and the falling snow.
When debris is covered by snow, all surface fluxes are computed similarly to the standard Crocus simulation. We use the parameterization of snow albedo introduced by Reference LejeuneLejeune and others (2007), which serves to lower albedos, α sd, of the snowpack under a given critical thickness:
where α d (j) and α s (j) are the albedo of debris and snow, respectively, in a given spectral band j. Fα is a weighting function defined as
where h sd is the total snow depth over the debris, is the critical snowpack thickness under which the debris has an impact on the snowpack albedo () and x is an empirical coefficient (x = 0.33) (Reference LejeuneLejeune and others, 2007). This approach allows us to account for the fact that a thin snowpack on top of the debris cover or on top of the ground is influenced by the underlying absorbing surface, effectively lowering its albedo, and also represents the fact that thin snowpacks are often uneven and patchy, which tends to accelerate their melting.
Handling of numerical layers. Crocus features an evolved scheme to dynamically adapt the vertical grid mesh to the snowpack stratification (Reference Brun, David, Sudul and BrunotBrun and others, 1992; Reference WillemetWillemet, 2010; Reference VionnetVionnet and others, 2012), leading to automated splitting and aggregating of the numerical snow layers during a model run, depending on the similarity of their physical properties and on the vertical discretization of the snowpack. This feature is switched off for debris layers, i.e. the latter do not undergo any splitting/aggregation throughout a given model run, and keep a constant thickness. The thickness of the numerical debris layers has to be chosen sufficiently small to accurately solve the heat diffusion equation, with a minimum imposed by the numerical scheme and the time-step of the model. A good compromise was found to be a numerical debris-layer thickness of 1 cm and an internal model time-step of 15 min. In this work, the debris cover was split into as many 1 cm thick layers as necessary. Note that the thickness of each debris-cover layer can be set individually, and can be adapted depending on the application.
Below-debris temperature. The reference implementation of Crocus-DEB solves the heat diffusion equation (including the phase change in snow and ice) through the whole snow/ debris/ice continuum. However, for the purpose of comparing Crocus-DEB with previous models of debris-covered glaciers assuming that ice below the debris is at the melting point (e.g. Reference Nicholson and BennNicholson and Benn, 2006; Reference Reid and BrockReid and Brock, 2010), a modified version of Crocus-DEB was developed in which the temperature of sub-debris layers is forced to the melting point. At each time-step, the temperature of below-debris layers is raised to the melting point, and the corresponding amount of energy artificially added to the system is recorded. This model version is referred to as Crocus-DEBm below, for the sake of brevity.
Field and numerical experiments at Col de Porte, France
Field experiment
The research station Col de Porte (CDP) was used to test and evaluate Crocus-DEB. The CDP experimental snow and meteorological station has been recording the physical properties of the seasonal snowpack since 1959, along with the meteorological conditions. In addition to a detailed description of the site, Reference MorinMorin and others (2012) provide a quality-controlled snow and meteorological dataset spanning the period 1 August 1993 to 31 July 2011 at the site, at hourly time-steps, which includes the four components of the surface radiative balance, air temperature and relative humidity, wind speed and direction, precipitation (amount and phase), snow depth and snow water equivalent (SWE), and surface temperature. In addition, at weekly time resolution, in situ snow-pit observations are manually performed at three locations within the field site to measure the SWE, as well as the vertical profile of the physical properties of the snowpack (grain type, density, temperature, liquid water content). Due to differential shading and occasional snowdrift, the snowpack can be spatially heterogeneous (Reference MorinMorin and others, 2012). Nevertheless, the whole area is considered here as the reference debris-free test area (TA), TA0, for our experiment.
Between February and April 2011, inside the CDP experimental field, a dedicated field experiment was carried out to investigate the impact of a well-characterized debris cover on the underlying snowpack, and to provide data to drive and evaluate the model. The seasonal snowpack was used as a surrogate for the ice mass underlying the debris in the case of a real debris-covered glacier. The gravels were deposited on top of the snowpack on 1 February 2011 at 12:00. One of the two test areas of 1.5 m × 3 m was covered with a 4 cm thick layer, the other with a 15 cm layer; these ‘debris-like’ layers consisted of rounded gravel (mixed rock types including igneous, sedimentary and metamorphic rocks extracted from the Isère river in a quarry near Grenoble, France) with a density of 1460 kg m−3 and a bead size ranging from 10 to 20 mm. These two test areas are referred to hereafter as TA4 and TA15, respectively, and their results are compared with those of the reference debris-free area, TA0. Figure 2 provides an overview of the experimental set-up. Sensors were placed to continuously measure the temperature at the bottom, in the middle and at the top of artificial debris layers (Pt100 temperature probes), and the conductive heat flux at the bottom of the debris (Hukseflux HFP01). In addition, incoming and reflected shortwave radiation were measured on top of TA15 with a pyranometer (Hukseflux NR01). From 8 February 2011 at 16:00 to 7 April 2011 at 08:00, approximately weekly visits were made to the study site to carry out manual measurements of snow depth and SWE above and below the artificial debris, as well as on TA0. Manual measurements were carried out on one side of the TA, while automated measurements were carried out on the opposite side, to minimize disturbance to the continuous measurements. Manual snow corings providing the sum of SWE below and above the debris were performed on 8, 11 and 24 February and 21, 25 and 29 March 2011. Measurements limited to the upper snow layers above the debris were additionally carried out on 24 February and 11 and 15 March 2011.
Numerical simulations
The Crocus-DEB model is driven using the same meteorological forcing as Crocus (Reference Brun, David, Sudul and BrunotBrun and others, 1992; Reference WillemetWillemet, 2010), i.e. hourly records of air temperature, T a, relative humidity, H a, incoming longwave and shortwave radiation, L in and S in, wind speed and snow- and rainfall. Model simulations were carried out on TA0, TA4 and TA15, to numerically reproduce the conditions of the field experiment held at CDP. The values for roughness length, z 0, s, and thermal emissivity, εs, for snow were the usual values of 0.004 m and 1.0, respectively. The thermal emissivity, εd, thermal conductivity, k d, and specific heat, C p,d, for the debris were assigned the values 0.95, 0.7 W m−1 K−1 and 950 J kg−1 K−1, consistent with Reference Reid and BrockReid and Brock (2010) for their study on Miage glacier, Italy. The values of density, ρ d, and albedo, α d, of the debris were determined experimentally and the obtained values of 1460 kg m−3 and 0.2, respectively, were used for the numerical simulations. The value of the roughness length when debris is at the surface, z 0,d, was set to a value of 0.012, which is higher than that of snow but takes into account the fact that the test areas were surrounded by snow. Given that turbulent heat fluxes at the surface/atmosphere interface operate over spatial scales of a few metres at least, it was necessary to account for spatial heterogeneity of the surface in setting the z 0,d value. Nevertheless, results were not significantly modified when the z 0,d value was varied within a reasonable range. The main adjustable parameters used for the simulation at CDP are summarized in Table 1.
The model runs were set up and initialized as follows. First, a full model run without any debris was carried out from 1 September 2010 at 00:00 to 12 April 2011 at 0:00, corresponding to the reference simulation on TA0. In order to carry out the simulation including the presence of debris, on 1 February 2011 at 12:00, i.e. on the date when artificial debris was poured on top of the snowpack, the simulated vertical profile of the physical properties of the snowpack was altered by adding the required number of debris-cover layers to simulate a 4 and a 15 cm thick cover. The initial temperature of all debris layers was 273.15 K. Crocus-DEB model runs were then performed for TA4 and TA15 conditions starting on 1 February 2011 at 12:00, until 8 February 2011 at 16:00, the date of the first field measurement of the total snow depth and SWE below the artificial debris cover. The simulated profiles on 8 February at 16:00 were adjusted for TA4 and TA15, to match the observed snow depth and SWE below the artificial debris cover on each TA. Model runs for TA4 and TA15 were then initialized using the distinct starting conditions on 8 February at 16:00 and run until 12 April 2011. The same procedure was carried out when using Crocus-DEBm.
Field and numerical experiments at Changri Nup Glacier
In addition to the dedicated field experiment carried out at CDP over a temporary snow cover, a full in situ dataset obtained on Changri Nup Glacier, Nepal, was used to further evaluate Crocus-DEB. This glacier is located a few kilometres west of Khumbu Glacier (Everest area) and was equipped with an automatic weather station (AWS) on 31 October 2010, providing all the driving data needed for Crocus-DEB, i.e. half-hourly records of incoming shortwave and longwave radiation (Kipp & Zonen CNR4 sensor), air temperature and relative humidity (artificially aspirated Vaisala HMP45C sensor), wind speed (Young 05103 sensor) and snowfall occurrence (Ultrasonic SR50 sensor). Additionally, Pt100 temperature probes, referred to as T2.5, T5, T7.5 and T10, were inserted within the 10 cm thick debris layer, at 2.5, 5, 7.5 and 10 cm below the surface, respectively. Probes were inserted to be as much as possible in contact with both clasts and air pockets, in order to be as representative as possible of the physical state of the debris. Figure 3 provides a picture of this AWS just after its installation, showing the overall aspect of the debris layer in the vicinity. The debris particles are mainly of granite type and their size varies from <1 cm to a few tens of centimetres, making a continuous heterogeneous debris cover at this elevation (5360 m a.s.l.). Below the AWS, where the temperature probes were inserted, the debris is made up of homogeneous centimetric broken gravels, roughly similar to those artificially poured at CDP. The adjustable parameters of Crocus-DEB for the Changri Nup simulation are summarized in Table 1. The main difference between this and the CDP simulation is the value of the roughness length, which is higher than at CDP, because the measurements were carried out in the presence of debris of various sizes near the measurement site, while at CDP the test areas were surrounded by snow. The simulation was initialized using 30 m of ice, overlaid with 10 cm of debris with a flat initial temperature profile of 273.15 K in the debris and in the ice underneath. This approach is similar to the method used successfully in previous investigations (Reference Gerbaux, Genthon, Etchevers, Vincent and DedieuGerbaux and others, 2005; Reference LejeuneLejeune, 2009; Reference Dumont, Durand, Arnaud and SixDumont and others, 2012) using Crocus to simulate the mass balance of mountain glaciers.
Results – col de porte
Meteorological conditions and general overview of the field experiment
Figure 4 provides an overview of the meteorological and snow conditions (on TA0) encountered at CDP from 8 February to 11 April 2011. Figure 5 shows an overview of the temperature at the top and bottom of the debris layers on TA4 and TA15, as well as the albedo and snow depth on top of the debris on TA15. The winter of 2010/11 was exceptionally dry, with little and light snowfall. From December 2010 to April 2011, cumulative precipitation was only 382 kgm−2, the lowest value recorded at CDP since 1960 (for comparison, the mean value for the period 1971–2000 is 817 kg m−2). On 8 February 2011 the snow depth was only 0.28 m (Fig. 4d), which is only 30% of the mean value for the 1960–2010 time period at this date (1.10 m). The SWE of the snowpack was 140 kg m−2 (Fig. 4e). These peculiar conditions during winter 2010/11 were appropriate for carrying out our experiment. Indeed, due to the absence of snowfall between 13 January and 1 February, the snowpack was compact enough (ρ ≃ 360 kgm−3) to receive gravels without being physically or thermally disturbed. Consequently, at the beginning of the experiment, the snowpack below the debris on TA4 or on TA15 can be considered physically and thermally similar to the snowpack on TA0. During the 50 days of the experiment, cumulative precipitation was only 60 kg m−2, mostly occurring during two main snowfall events (Fig. 4c), the first on 20 February between 02.00 and 16.00 UTC (10 kg m−2) and the second between 26 February, 19.00 UTC and 28 February, 0.00 UTC (48 kg m−2). After the snowfalls, the snow depth increased by a few centimetres (Fig. 4d), and concurrently the albedo of the debris-covered surface recorded a significant rise (Fig. 5c). After these snowfalls, debris was temporarily covered by snow. Such periods with snow over debris were characterized by a high albedo (Fig. 5c) and a strong attenuation with depth of the daily amplitude of the surface temperature, due to the insulating effect of snow, as observed between 20 and 25 February for TA4 and TA15 or from 27 February to 16 March for TA4 or even 18 March for TA15 (Fig. 5a and b).
Besides its interest for evaluating the Crocus-DEB model, this experiment provides interesting information, per se, in terms of the impact of debris on the melt rate of the underlying snow. According to the values reported by the temperature sensors placed at the bottom of the debris cover (in particular the first occurrence of a value above the melting point), the total melt below the debris occurred on 29 March (Fig. 5a) and 4 April (Fig. 5b) for TA4 and TA15, respectively. The complete melting of the snowpack was thus quicker on TA4 than on TA15, which is in line with the fact that a thicker debris layer tends to suppress melting (Reference ØstremØstrem, 1959; Reference Kayastha, Takeuchi, Nakawo and AgetaKayastha and others, 2000). The snow cover on TA0 disappeared between 29 March and 3 April, depending on the measurement chosen (e.g. automatic snow depth, manual snow depth, visual inspection) (Fig. 4d and e). This indicates that the critical thickness relevant to our experimental conditions lies between 4 and 15 cm. However, the variability of the melt date on TA0 stresses the need for a well-defined debris-free reference against which melt rates below the debris can be compared.
Delineation of four typical periods
Based on the meteorological conditions and the presence or not of a transient snowpack on top of the debris, the measurement campaign was split into four periods (Figs 4–7; Table 2), and summarized along with their typical meteorological characteristics. Due to the peculiar conditions that prevailed during the 2010/11 winter, these periods were exceptionally homogeneous in terms of their meteorological conditions.
Overview of the model performance throughout the experimental campaign
For both TA4 and TA15, the performance of Crocus-DEB was analysed by comparing observed and simulated snow depth, SWE, and internal and basal temperatures of the debris, as well as heat fluxes at the bottom of the debris layers at hourly or daily timescales, until there was no snow left on the ground. Daily heat flux values at the bottom of the debris layers were cumulated over the whole experimental period, as an additional proxy for the total energy available to melt the underlying snow. Figures 6 and 7 provide an overview of these measured and simulated variables for TA4 and TA15, respectively. Below we summarize the main results of these comparisons.
Snow on top of debris
Table 3 provides measured and simulated SWE and snow-depth data above the debris on three dates during the experimental campaign, after the two significant snowfalls on 20 February and 26–27 February (Figs 4c and 5c). It shows that these two integrated snow properties are well simulated by Crocus-DEB. These snowfalls are characterized by a SWE increase visible in Figures 6a and 7a.
Total snow water equivalent
In terms of total SWE (i.e. taking into account snow below and above the debris layers), there is overall good agreement between the model results and the observations (Figs 6a and 7a). On TA4, observed and calculated SWE agree well throughout the whole period of the experiment, except on 24 February, when the measured and simulated SWE are 145.5 and 114.7 kg m−2, respectively. The agreement is very good in March, leading to a simulated complete melting of the snow cover on 30 March that is within 1 day of the total melt date identified using temperature measurements. On TA15, in March, the melting rate is well reproduced by the model but there is systematic overestimation of the SWE (of 32, 19 and 23 kg m−2 on 21, 25 and 29 March, respectively), leading to complete disappearance of the snow on 10 April, 6 days later than the field observations.
Albedo
Figure 7a shows the observed and simulated albedo on TA15. It shows overall agreement between the two, with lower values, ∼0.20, in the absence of snow, and higher values, up to 0.8, in the presence of snow at the surface. Consistent with other evaluations of Crocus, in the presence of snow the simulated albedo is generally higher than the observations (despite the implementation of the formula reducing the snow albedo for shallow snowpacks; see above). This may explain why the complete melting of the snow on top of the debris is delayed by a few days in the simulation (cf. the observations).
Heat flux and temperature within the debris
Heat fluxes and temperatures at the top, in the middle and at the bottom of the debris-cover layers are reasonably well simulated by the model for both debris thicknesses (Figs 6b–f and 7b–f). Indeed, the model satisfactorily captures the dynamics and range of variation of the temperature and heat flux within and at the boundary of the debris cover. The main discrepancies occurred on 25 and 26 February. For these days, on both TA4 and TA15, the model overestimates the melt rate of the thin snow layer occurring on top of the debris for just a few hours (total snowfall precipitation of only 10 kg m−2), explaining why temperatures and heat fluxes then show a large daily variability, in contrast to the observations. Consequently, the calculated melting during these two days is 20 and 10 kg m−2 over TA4 and TA15, respectively, although it was insignificant from the measurements. Between 2 and 13 March, simulated temperatures were below the melting point, leading to slightly negative daily basal fluxes, although the debris always remained at or slightly below freezing. Consequently, the model erroneously simulated a cooling of the snow below the debris during period P2-2. Finally, the relative differences of the observed and calculated diurnal amplitudes of the debris temperatures and basal fluxes are higher over TA15 than over TA4 during period P1. This is also the case during period P3, but exclusively for debris temperatures and not for heat flux. Indeed, between 20 and 27 March anomalously low temperatures were simulated at the bottom of the debris over TA15, albeit daily basal heat fluxes were well simulated. Although the date of the complete disappearance of the snow is well simulated over TA4, the calculated cumulative heat flux on 29 March is 120 W d m−2 higher than the observed cumulative flux, i.e. a discrepancy of ∼20%. In contrast, over TA15, the simulated and measured cumulative fluxes are in good agreement, but the date of snow disappearance is simulated to be too late. These discrepancies clearly show that there are cases of error compensation in the model.
Mean diurnal cycle of thermal conditions during periods of homogeneous meteorological conditions
To follow the thermal evolution of the debris layers during the day, the mean diurnal cycle of the observed and simulated temperatures at the top, in the middle and at the bottom of the debris layer were calculated, as well as the heat fluxes at the bottom of the debris for every period given in Table 2. Since diurnal cycles are very similar between periods P1 and P3 (no snow on debris) and between periods P2-1 and P2-2 (debris covered by snow), only the mean diurnal cycles of periods P2-2 (17.5 days) and P3 (11 days) are reported in Figures 8 and 9, respectively. Table 4 summarizes the statistics of temperatures and fluxes observed and calculated on both debris-covered test areas during the four periods considered.
Snow covering the debris (period P2-2)
During period P2-2, i.e. when snow was covering debris on TA4 and TA15, field data confirm that the presence of snow had a particularly strong attenuating effect on heat transfer, as seen in Figure 8, where the snow layer was 10–30 cm thick. At TA4 (Fig. 8a), temperatures at the three measurement levels (top, middle and bottom) always remain slightly below the melting point.T bottom–4–meas ranged between 272.9 and 273.0 K with a daily mean of 273.0 K. Daily means of T mid–4–meas and T top–4–meas were 273.1 and 273.0 K, respectively, and their respective daily amplitudes of 0.1 and 0.2 K were barely higher than that of T bottom–4–meas. At TA15, debris always remained at the melting temperature (Fig. 8b) and the basal flux, Q c–bottom–15–meas, was 0 (Fig. 8f). At TA4 the daily basal flux was only −0.2 W m−2 (Fig. 8e).
The results of the simulation during period P2-2 are given in Figure 8c, d, g and h. They show that the numerical simulations overestimate the diurnal amplitude of debris temperature. At both TA4 and TA15, simulated temperatures are always below the melting point, and minimal at midday. Temperature increases from the top of the debris layer toward the bottom. Although measured basal fluxes were always 0, or very close to 0, simulated fluxes are always negative and range between −1 and −3.5 W m−2 at TA4 and between 0 and −2.5 W m−2 at TA15. The model tends to simulate the downward propagation of a ‘cold wave’, although measurements indicate that the snow below the debris is almost perfectly insulated from the surface.
Melting period without snow on top of the debris (period P3)
During period P3, corresponding to the main melt period without snow on top of the debris, experimental results (Fig. 9a and b) indicate that the surface temperatures measured on the debris over both test areas (T top–4–meas and T top–15–meas) minimize and maximize at the same time of the day, at around 5.00 and 13.00 UTC. The diurnal amplitude of T top–4–meas is 22 K, slightly lower than that of T top–15–meas (26 K). At the surface, the debris follows a thermal cycle, which appears to be similar regardless of the debris thickness, even if the maximum temperature is 5 K higher over TA15 than over TA4. In the middle of the debris layer, the diurnal amplitude of the temperature measured at -2 cm on TA4 (T mid–4–meas) is 15 K, three times higher than at −7.5 cm on TA15 (T mid–15–meas) (Fig. 9a and b). Also visible is a clear delay for T mid–15–meas compared to T mid–4–meas, the minimum and maximum of T mid–15–meas being reached at 9.00 and 16.00 UTC, respectively, ∼4 and 3 hours after T mid–4–meas, respectively. The heat flux propagating within the debris is then all the more attenuated and delayed as the debris is deeper in the layer. At the bottom of the debris layer, sensors were installed at the snow/debris interface, preventing the temperature from exceeding the melting point. At the bottom of TA15, the temperature, T bottom–15–meas, always remains at 273.15 K (melting), although on TA4, underlying snow melts for only 6 hours a day, between 10.00 and 16.00 UTC. The rest of the time, T bottom–4–meas is below the melting point, reaching its minimum value at 6.00 UTC (271.3 K). The underlying snow thus alternates between thawing and refreezing along the diurnal cycle.
The measured conductive heat fluxes at the snow/debris interface on TA4 and TA15, Q c–bottom–4–meas and Q c–bottom–15–meas, respectively, are plotted in Figure 9e and f. On TA4, Q c–bottom–4–meas is positive during daytime between 8:00 and 18:00, with a maximum value at 14:00 (100 W m−2), and remains slightly negative (∼−7 W m−2) the rest of the time. On TA15, Q c–bottom–15–meas always remains positive, with significant values only between 10:00 and 24:00, a maximum of 50 W m−2 being reached at around 16:00. Even if Q c–bottom–4–meas is negative at night, which means that refreezing occurs under TA4 debris, its daily mean is 19 W m−2 (cf. 14 W m−2 for Q c–bottom–15–meas), indicating that daily melting is higher below TA4 than below TA15.
Figure 9c, d, g and h show the results of the simulations. Comparison of Figure 9a and c forTA4, or Figure 9b and d for TA15, shows that simulated temperatures agree fairly well with observations, as also confirmed by correlation coefficients in Table 4 (r 2 values 0.74–0.95, except for temperature at the bottom of TA15 where r 2 values are only 0.41 but the temperature variations are small). Nevertheless, the model fails to consistently reproduce the maximal heating at the surface of TA4 and TA15, with underestimations as high as 7 and 9 K, respectively. At the surface and in the middle of TA4, mean daily values of T top–4–sim and T mid–4–sim are 1.3 and 0.6 K, respectively, lower than the measurements (Table 4). At TA15, daily T top–15–sim is 1.3 K lower than T top–15–meas but daily T mid–15–sim is 0.3 K higher than T mid–15–meas. The diurnal amplitude of T mid–4–sim is also weaker than that of T mid–4–meas (Fig. 9a and c), although that of T mid–15–sim is higher than that of T mid–15–meas (Fig 9b and d). Finally, at the bottom of both debris layers, daily mean simulated temperatures, T bottom–4–sim and T bottom–15–sim, are slightly lower (∼0.2 K) than observations. There is good agreement between T bottom–4–sim and T bottom–4–meas with observed and simulated melting conditions encountered between 10.00 and 16.00 UTC, but T bottom–15–sim remains slightly below the melting point between 4.00 and 12.00 UTC, although observations show that underlying snow is still melting.
The comparisons between Figure 9e and g and between Figure 9f and h show that the model is able to consistently reproduce the timing of the daily cycle of the heat flux at the bottom of TA4 or TA15 (with r 2 between simulated and measured fluxes as high as 0.97 and 0.88, respectively; Table 4), but fails to properly simulate its daily amplitude. Indeed, the daily amplitude of Q c–bottom–4–sim (220 W m−2) is twice as high as that of Q c–bottom–4–meas and the daily amplitude of Q c–bottom–15–sim (75 W m−2) is one-third higher than that of Q c–bottom–15–meas. Looking at daily means, the agreement is much better, with daily values of Q c–bottom–4–sim and Q c–bottom–4–meas equal to 17 and 19 W m−2, respectively, and daily values of Q c–bottom–15–sim and Q c–bottom–15–meas equal to 3 and 14 W m−2, respectively.
Impact of assumptions about under-debris temperature
Model runs using the reference version of Crocus-DEB and the modified version setting the temperature of the snow below the debris to the melting point (Crocus-DEBm) were compared in terms of net heat flux at the debris/belowdebris interface, also considering the additional (artificial) amount of energy used to raise the temperature below the debris to the melt temperature. The latter was converted to a flux, allowing comparison with the net interface flux, although, strictly speaking, it does not correspond to a flux through an interface, because in Crocus-DEBm energy may be artificially added to any layer below the debris. (Note that as the Changri Nup Glacier simulation does not correspond to the melt season, it is unreasonable to compare Crocus-DEB model runs with Crocus-DEBm there.) At CDP, total melt simulated by Crocus-DEBm occurred 5.7 and 4.3 days before the standard Crocus-DEB simulation at TA4 and TA15, respectively. Table 5 shows the results in terms of heat flux at the debris/underlying-snow interface and basal runoff, respectively, at CDP. During periods P2-1 and P2-2, i.e. in the presence of snow on top of the debris, heat fluxes at the debris/underlying-snow interface are negligible (<1 W m−2) for both models, although the energy added to the underlying snow corresponds to 2 W m−2 or more. During the cold, dry period P1, discrepancies are higher between the two versions of the model and between TA4 and TA15. Using Crocus-DEBm decreases the net heat flux through the debris (as much as halving it in the case of TA4), while the additional equivalent heat flux is of the order of twice the heat flux through the debris simulated by Crocus-DEB. However, it should be noted that this period does not correspond to a melt period, so we anticipate that melt-period debris-cover models will not be used under such conditions. During the main melt season (P3), the heat flux at the bottom of the debris on TA4 using Crocus-DEBm is less than half that of Crocus-DEB, although Crocus-DEBm leads to a total melt occurring several days earlier than Crocus-DEB. Using Crocus-DEBm thus leads to overestimation of the melt rate below the debris, even if this cannot be concluded when looking only at the net heat flux through the interface below the debris. The heat flux at the bottom of the debris on TA15 is rather similar in the case of Crocus-DEBm and Crocus-DEB, consistent with the fact that the temperature at the base of a thicker debris cover exhibits less variation and, in particular during the melt period, the melt operates more continuously (Figs 6 and 7). Over the full measurement period, the same situation as during P3 is found, consistent with the fact that this period is the longest and also that much of the net flux through the debris occurs during the melt period. In terms of basal runoff, the use of Crocus-DEBm leads to a simulated basal runoff, which is directly related to the melt flux (but accounts for the impact of rain), that is -10% higher during the melt period than the results of Crocus-DEB. During periods not dominated by melt, the discrepancy is even higher, probably because, in addition to increasing the melting below the debris, Crocus-DEBm does not allow any meltwater to refreeze within the debris layers. When liqud water is present in the ice/snow layer located just above the debris layers, it is instantaneously drained through the debris to the ice/snow layer just below the debris layers. Note that period P1, which was not simulated to be a melt period using Crocus-DEB, owing to cold and dry meteorological conditions in February, becomes a period with significant melt using Crocus-DEBm, with on average 8 kg m-2 d−1 runoff flux below TA4.
Results – Changri Nup
Meteorological conditions
An overview of the meteorological conditions between 16 November and 31 December 2010 at Changri Nup Glacier is given in Figure 10. This shows that during this period the glacier was exposed to typical winter dry conditions prevailing in this part of the Himalaya (Reference Bookhagen and BurbankBookhagen and Burbank, 2010), i.e. no precipitation, almost no clouds, moderate wind speed (2.6 m s−1), low air temperature (265.95 K) and humidity (19.1%), high incoming solar radiation (194 W m−2) and low incoming longwave radiation (201 W m−2). This period was selected because the specific and extremely homogeneous conditions encountered make it challenging for the model to reproduce the thermal behaviour of the debris layer as well as the ice/ debris interface. Indeed, there was no snow cover over the debris (albedo = 0:27), favouring the absorption of solar radiation by the debris surface during daytime, and causing much energy loss during the night through outgoing long-wave radiation. Consequently, energy exchanges across the ice/debris/atmosphere continuum are high and very variable over a day, as evidenced by the large daily amplitude of air temperature (sometimes >15 K; Fig. 10).
Mean diurnal cycle of heat transfer within the natural debris layers
Figure 11 shows the hourly evolution of measured and calculated temperatures of the debris over the entire study. As seen before, the model yields good results for simulating the thermal evolution of the debris, especially near the surface (at −2.5 cm) and at the ice/debris interface (at −10 cm). To study the ability of the model to simulate the day-to-day variability of the debris temperatures, we selected two contrasting days: the warmest of the period (25 November 2010; daily mean air temperature 272.75 K) and the coldest (19 December 2010; daily air temperature 259.15 K). Table 6 gives the mean measured and simulated temperatures of the debris for both days. On 25 November, temperature decreases from the surface of the debris toward the bottom, and the profile is reversed during the cold day (19 December). The fact that in both cases vertical temperature profiles (Fig. 11) and daily temperatures (Table 6) are well simulated gives confidence in the ability of the model to properly reproduce the day-to-day variations in the thermal conditions in the debris. Nevertheless, there are some discrepancies between observations and simulations. Mainly at the beginning of the study period, at −2.5 and −5.0 cm simulated maximal temperatures often exceed the melting point, sometimes by more than 10 K, although measured temperatures do not exceed 273.15 K at −5 cm and are limited to 280 K at −2.5 cm. At −5 cm, the simulated daily amplitude is also somewhat overestimated compared with measurements.
The mean diurnal cycles over the study period (5 November to 30 December) of measured and simulated temperatures at various depths within the debris (−2.5, −5.0, −7.5 and 10 cm) on Changri Nup Glacier are shown in Figure 12. A statistical summary is also presented in Table 7. The warming of the debris follows the diurnal cycle of incident solar radiation. The temperature at −2.5 cm increases from 7:00 to 13:00 and decreases quickly until the sunset and then slowly but continuously until sunrise. In daytime, a warm front progressively propagates within the debris layer from the surface, with a delay of temperature maximum and attenuation of the diurnal thermal amplitude, more pronounced for deeper debris (Fig. 12a). At night, the opposite pattern is observed, as a cold front propagates into the debris. A comparison between Figure 12a and b, and the statistical analysis in Table 7 shows that Crocus-DEB is able to consistently reproduce the thermal evolution of the debris, but tends to slightly overestimate the diurnal amplitude at −7.5 cm and at the bottom (4.1 and 2.9 K, respectively).
Discussion
Ability of the model to simulate the specific processes driving the melting of debris-covered glaciers
In this study, the snow model Crocus was adapted to deal with a debris layer as observed on debris-covered glaciers. Crocus-DEB is therefore a physically based model able to simulate the energy and mass fluxes across the ice/debris/ snow/atmosphere continuum.
To thoroughly test the performance of the model, Crocus-DEB was applied to and calibrated on 4 and 15 cm thick debris layers, artificially poured over a seasonal alpine snow cover where meteorological forcing and snow conditions are continuously extensively measured, before being directly applied (without extra calibration) over a natural Nepalese debris-covered glacier. Over this glacier, we selected a dry period (no precipitation and debris exposed at the surface) in November and December 2011, with highly variable meteorological conditions between nights and days, and between the beginning (melting conditions encountered for a few hours in daytime) and the end (cold conditions and no melting) of the period (Fig. 12). On the whole, there is good agreement between observations and simulations, while considering the temporal evolution of the SWE above the ground or debris temperatures or heat fluxes at the bottom of the debris (Figs 8, 9 and 12; Tables 3, 4 and 6). This means that the main processes that control the melting of the underlying snow/ice layers below the debris are well simulated by the model, i.e. that the conflicting effect between the enhanced absorption of incoming solar radiation, due to the low albedo of the debris, and the insulating effect is well reproduced by the model. However, several discrepancies between measurements and simulations have also been shown, in particular regarding the magnitude of the hourly variations of temperature and heat flux within or at the bottom of the debris layer (e.g. Figs 9 and 12). Daily integrations of the heat flux at the bottom of the debris cover are well simulated at Col de Porte; we further note that, as soon as the temperature remains below the melting point throughout the day at the Changri Nup site, the agreement between measurements and simulations becomes almost perfect (Fig. 11d). It appears very likely that the discrepancies reported here stem from the inadequate handling of phase change and hydraulics within the debris. As indicated above, one key assumption of the model is that liquid water entering the debris layer is immediately transferred below it. In reality, even if it is very coarsely grained, the debris cover is able to retain liquid or ice water in interstitial spaces, often as a coating covering debris surfaces. Diurnal phase changes of the water film covering debris are likely to induce significant disturbances in the time variations of the energy budget of the debris within one day; however, as long as the liquid water flux is small, freeze/thaw cycles should compensate each other, leading to an insignificant net effect on the total energy transferred through the debris. This explains why, even if there are some discrepancies between hourly variations of the temperature or heat flux within or at the base of the debris layers, the net daily integrated simulated flux remains consistent with observations. Note that such an issue is not specific to Crocus-DEB and is likely to also have a significant impact on previously existing physically based debris-cover models (Reference Nicholson and BennNicholson and Benn, 2006; Reference Reid and BrockReid and Brock, 2010). Various other reasons can be invoked to explain the observed discrepancies, including, for example, non-conductive processes within the debris cover, which were not accounted for in the present study (Nicholson and Benn, in press). More work is needed to incorporate a physically consistent phase-change and water percolation scheme within the debris in Crocus-DEB.
Comparison with alternative models
Besides the ability to simulate the inception, development and melt of a transient snowpack on top of a debris layer, which directly derives from the use of the detailed snowpack model Crocus, the main difference between Crocus-DEB and existing numerical models applied to the melt of debris-covered glaciers (Reference Nicholson and BennNicholson and Benn, 2006; Reference Reid and BrockReid and Brock, 2010; Reference Reid, Carenzo, Pellicciotti and BrockReid and others, 2012) is the explicit simulation of the temperature profile below the debris. Previous models are generally only applied during the melt season, and assume that the temperature at the base of the debris is the melt temperature. In Crocus-DEB, this hypothesis is relaxed, which allows us to simulate the surface energy and mass balance year-round. Numerical experiments were carried out to investigate the impact of this assumption on the surface energy balance. In the case of a thick debris cover (15 cm), the impact of this assumption during the melt season is limited, although it is detectable. When the debris thickness is lower (4 cm in the tested example), even during the melt season the net heat flux through the debris is underestimated by a factor of two when the temperature profile below the debris is imposed to remain at the melting temperature. In this case, the artificial addition of energy to the below-debris snow layers exceeds the heat flux through the debris. In such a case, the adjustment in surface energy and mass-balance models of some of the physical properties not easily measured in the field, such as the effective thermal conductivity of the debris, may be prone to errors due to the compensation of an inadequate assumption on the basal temperature of the debris. We also note that this assumption, if implemented in algorithms used to map the physical properties of debris cover from remotely sensed thermal data, would lead to systematic errors in the retrieved properties of the debris, especially when the debris cover is thin (i.e. less than −15 cm).
An example of model application: a study of the critical thickness at CDP
The ability of the model to represent the opposing effects of the debris cover, iconically summarized through the concept of the critical thickness (Reference ØstremØstrem, 1959), was tested using model runs at CDP. Crocus-DEB runs were performed using the same method as described above for TA4 and TA15, but for a series of debris thickness values ranging from 1 to 30 cm in steps of 1 cm. Moreover, three additional runs were done without debris. For these latter runs, the broadband albedo of all the initial snow layers (i.e. snow layers older than 8 February 2012) was forced to 0.4, 0.3 and 0.2, respectively. These runs mimic the temporal evolution of the surface energy and mass balance of debris-free glaciers with varying bare-ice albedo. The model was run from 8 February 2011 with identical snowpack properties below the debris and for the debris-free runs. The initial total SWE of the underlying snowpack or the debris-free snowpack was taken equal to 150.6 kg m−2. From these model runs, daily ablation rates integrated over various periods were computed, and the date of the total melt was extracted. These results are shown in Figure 13, which can be viewed as various flavours of so-called Østrem curves (Reference ØstremØstrem, 1959; Reference Reid and BrockReid and Brock, 2010).
Figure 13 confirms that, as expected, melting is accelerated for thinner debris layers, and slowed for thicker layers, in full agreement with the Reference ØstremØstrem (1959) curve. Nevertheless, this artificial example illustrates several facts relevant to the concept of the critical thickness. In particular, significantvariations in critical debris thickness may, in some cases, simply stem from a change in reference debris-free glacier albedo. In addition, the critical debris thickness computed can differ widely, depending on the period chosen, and in some cases it cannot even be defined (e.g. when the melt rate under the thinnest debris layer is lower than the melt rate in the debris-free area). This occurs in our case when the albedo of the debris-free area is of the same order of magnitude as the albedo of the debris. The dependence of critical debris thickness on meteorological conditions that we observed experimentally is also revealed by the numerical experiments, and confirms that this concept, while extremely powerful in explaining concisely the physical impact of debris covers on glacier energy and mass balance, should not be used further than for pedagogic purposes.
Conclusion and outlook
In this study, we have introduced the physically based Crocus-DEB model able to simulate year-round the energy budget of the surface and the mass balance of debris-covered glaciers, taking into account not only snow layers that can deposit over the debris, but also the thermal evolution of the snow/debris/ice continuum without any assumption about the temperature at the base of the debris layer. The model was evaluated against data from a series of field experiments using an artificial debris-covered snow-pack, as well as in situ data from the debris-covered glacier Changri Nup, in the Himalaya. The model performs well in terms of surface energy balance and net daily heat flux through the debris, and can be used at the spatial scale of the whole glacier to study the impact of the debris-cover layer on the glacier energy and mass balance. Improving the handling of vertical water transfer and phase changes within the debris cover is the next step of model development, which will improve the time evolution of the temperature and heat flux within and at the boundaries of the debris cover. The model was also used to briefly explore the concept of a critical debris thickness. The critical debris thickness depends strongly on the reference debris-free glacier, in particular its ice broadband albedo, and on the meteorological conditions, implying that defining a sole value of the critical debris thickness for a given glacier may be misleading.
Acknowledgements
We thank J.-M. Panel, B. Lesaffre and D. Poncet for support during the field experiment at Col de Porte, and J.-M. Willemet for help in developing the model. This work was supported by the French Service d’Observation GLACIOCLIM and the Italian Consortium Ev-K2-CNR, which provided all the logistical support on Changri Nup Glacier. Pyramid and Ev-K2-CNR staff members are especially acknowledged for their help and their participation in the fieldwork in Nepal. Funding from INSU-LEFE and Agence Nationale de la Recherche (ANR-09-CEP-005-01/ PAPRIKA) is gratefully acknowledged. We thank the Scientific Editor, T. Moelg, T. Reid and an anonymous reviewer for very useful comments and suggestions.