A proxy-based method of estimating glacier heat sources and transport mechanisms to map thermal processes at a regional scale

Glacier thermal structure exerts a strong control on ice dynamics, yet it remains poorly characterized at the global scale. Despite its influence on glacier flow and, consequently, on melt and sea-level projections, thermal structure is typically neglected in large-scale glacier models. This omission stems from the fact that current methods of determining thermal structure require extensive field observations or computationally expensive modelling, resulting in only dozens of glaciers having well-defined thermal structures. In this work, we developed and applied a proxy-based approach to estimate the major heat sources (meltwater refreezing and strain heating) and transport mechanisms (advection and diffusion) in the heat transfer equation, which shapes glacier thermal structure. Proxies derived from publicly available observations and model output are calculated for glaciers in RGI-01 (Alaska) and RGI-07 (Svalbard and Jan Mayen). This framework enables a rapid regional assessment of thermal characteristics without requiring coupled thermomechanical modelling.Across both Alaska and Svalbard, meltwater refreezing is the dominant heat source for 73% of glaciers, while diffusion dominates heat transport for 75% of glaciers. Note that this model only considers refreezing in the accumulation area, and 21% of these glaciers do not have an accumulation area. The majority (76%) of glaciers exhibit warmer accumulation-area ice transported toward cooler ablation areas. When examined by proxy rank, the most common pattern (47% of glaciers) is characterized by refreezing as the dominant heat source over strain heating, and vertical diffusion over horizontal advection as the primary transport mechanism. This pattern represents 25% of glacierized area, corresponding primarily to smaller glaciers (