“…However, due to the complex boundary layer meteorology of mountainous areas and the general lack of detailed on-glacier measurements (Hanna et al, 2017), constant and linear LRs are commonly used for glacier ablation estimations, rather than distributed air temperature fields for glacier ablation estimations (Ayala et al, 2015). This is a major simplification, as it has been widely recognized that air temperature LRs are spatially and temporally variable in mountainous regions (Petersen & Pellicciotti, 2011), both on-glacier (Ayala et al, 2015;Hanna et al, 2017;Shaw et al, 2017) and off-glacier (Heynen et al, 2016;Shen et al, 2016). Many studies use off-glacier data that do not account for the variability of the air temperature associated with katabatic boundary layer flows and the damping and ice surface cooling effect observed over glacier surfaces (Ayala et al, 2015;Carturan et al, 2015;Petersen & Pellicciotti, 2011;Petersen et al, 2013;Shaw et al, 2016).…”