Ablation Stakes Research Articles

Sea ice mass balance during the MOSAiC drift experiment: Results from manual ice and snow thickness gauges

Precise measurements of Arctic sea ice mass balance are necessary to understand the rapidly changing sea ice cover and its representation in climate models. During the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition, we made repeat point measurements of snow and ice thickness on primarily level first- and second-year ice (FYI, SYI) using ablation stakes and ice thickness gauges. This technique enabled us to distinguish surface and bottom (basal) melt and characterize the importance of oceanic versus atmospheric forcing. We also evaluated the time series of ice growth and melt in the context of other MOSAiC observations and historical mass balance observations from the Surface Heat Budget of the Arctic (SHEBA) campaign and the North Pole Environmental Observatory (NPEO). Despite similar freezing degree days, average ice growth at MOSAiC was greater on FYI (1.67 m) and SYI (1.23 m) than at SHEBA (1.45 m, 0.53 m), due in part to initially thinner ice and snow conditions on MOSAiC. Our estimates of effective snow thermal conductivity, which agree with SHEBA results and other MOSAiC observations, are unlikely to explain the difference. On MOSAiC, FYI grew more and faster than SYI, demonstrating a feedback loop that acts to increase ice production after multi-year ice loss. Surface melt on MOSAiC (mean of 0.50 m) was greater than at NPEO (0.18 m), with considerable spatial variability that correlated with surface albedo variability. Basal melt was relatively small (mean of 0.12 m), and higher than NPEO observations (0.07 m). Finally, we present observations showing that false bottoms reduced basal melt rates in some FYI cases, in agreement with other observations at MOSAiC. These detailed mass balance observations will allow further investigation into connections between the carefully observed surface energy budget, ocean heat fluxes, sea ice, and ecosystem at MOSAiC and during other campaigns.

High temporal resolution records of the velocity of Hansbreen, a tidewater glacier in Svalbard

Abstract. Monitoring changes in glacial dynamics is essential for understanding the environmental response to accelerated climate warming in the Arctic. However, geodetic surveys in polar regions continue to present considerable challenges because of the harsh environmental conditions and the polar night. This study records a 14-year-long time series (2006–2019) of global navigation satellite system (GNSS) surveys of the positions of 16 ablation stakes distributed across Hansbreen, a tidewater glacier in southern Svalbard. The measurements were conducted with an exceptionally high temporal resolution, from about 1 week to about 1 month, and covering altitudes ranging from 20 to 490 m a.s.l. The position of one stake was surveyed every day. The primary data products consist of the stake coordinates and velocities. Time series of annual and seasonal velocities are also provided. This dataset may be a subject of further studies of glacier dynamics in relation to the long-term and seasonal impact of climate change on ice flow in the region. It also offers unique material for tuning numerical models of glacier dynamics and for validating satellite-derived products such as velocity and digital elevation models. The dataset described here has been made publicly available through the Zenodo repository: https://doi.org/10.5281/zenodo.8289380 (Błaszczyk et al., 2023).

Thinning and dynamics of the glacier terminus at the Pasterze Glacier (Austria), 2016–2021

As a result of global warming, glaciers worldwide recede and lose mass. Glacier ablation is influenced by surface characteristics, in particular the existence of debris cover. The evolution of glaciers may also be influenced by the existence of proglacial lakes, which are worldwide increasing in number and size. We examine the recent changes of the terminus of Pasterze Glacier by discussing the thinning and dynamics of the glacier terminus along with the changes of the glacial-proglacial transition zone between September 2016 and September 2021. By considering land cover in the interpretation, we find that glacier surface characteristics show different rates of elevation changes and surface movement. We mainly use orthophotos and digital elevation models (DEMs) resulting from airborne surveys of 12 different acquisition dates, but also draw on ablation stake measurements, lake water level data, water temperature data, as well as thermal infrared (TIR) imagery. While debris-covered ice lowered typically −5.0 m a−1 ± 4.6 m a−1 (mean ± standard deviation), debris-free ice was affected by ablation in the range of −7.0 m a−1 ± 2.7 m a−1, thus demonstrating that melt rates in similar elevations are related to surface characteristics. The elevation change rates were, in most cases, in a similar range throughout the observed periods of time. The horizontal displacement close to glacier-marginal crevasses amounted to up to ∼8 m a−1, but was mostly in the range of ∼1–4 m a−1. Considering both horizontal displacement and elevation changes, we find that clean and dirty ice surfaces show larger thickness losses and are moving faster compared to debris-covered glacier surfaces. Diurnal and seasonal fluctuations of lake water levels in the range of decimeters and lake water temperatures mostly clearly exceeding 0 °C are indicative of the potential influence of the proglacial lake on the glacier terminus. Given the decelerating and decaying behavior of the Pasterze Glacier in recent decades, we interpret the change rates of surface lowering and horizontal displacement 2016–2021 to be influenced by other factors such as the recently formed proglacial lake. This is also supported by direct ablation measurements that mostly show lesser ice melt differences compared to DEM-based surface elevation changes, suggesting that other factors such as thermally-induced melting contributed to an increased elevation difference.

Monitoring the annual geodetic mass balance of Bordu and Sary-Tor glaciers using UAV data

Abstract The geodetic mass balance of a glacier corresponds to glacier-wide volume changes, converted to mass changes using density assumptions. It is typically calculated by differencing multi-temporal digital elevation models. In this study, we show how the annual geodetic mass balance of a glacier can be derived from uncrewed aerial vehicle (UAV) data. The presented workflow is applied to two small- to medium-sized glaciers in the Kyrgyz Tien Shan (Central Asia): Bordu glacier and Sary-Tor glacier. The obtained geodetic mass balance is compared with the glaciological mass balance derived from a network of ablation stakes and snow pits. A previously calibrated mass-balance model is used to correct for the difference in acquisition dates. The results show that the determined geodetic mass balance matches closely with the glaciological mass balance. Besides, for both glaciers the geodetic mass balance does not seem to be particularly sensitive to the assumptions regarding volume-to-mass conversion. Therefore, our results demonstrate that UAVs can serve as a valuable instrument to quantify the annual geodetic mass balance and to validate the glaciological mass balance. The conventional glaciological mass-balance estimation often relies on interpolation and extrapolation methods, whereas UAVs offer the potential for direct data acquisition over the entire glacier surface.

Glaciological and meteorological investigations of an Alpine debris-covered glacier: the case study of Amola Glacier (Italy)

Debris-covered glaciers are common in many regions of the world, and accurately modelling their melt is of increasing importance for water resources planning, and biological and ecological research. In this study, we investigate meteorological and glaciological conditions and estimate the melt of Amola Glacier (a small debris-covered glacier in the Adamello-Presanella Massif, Italian Alps) using an empirical approach, based on shortwave radiation, surface temperature, debris thickness and thermal resistance. Meteorological conditions are determined from a supraglacial automatic weather station, while the model is calibrated using i) field data acquired during the ablation season 2020, including a network of ablation stakes and thermistors, ii) modelled solar radiation and iii) thermal imagery from Landsat 8 TIRS.The analysis of glacier meteorological conditions shows a high prevalence of cloud-covered (50.60% of daytime observations) and humid conditions, with a high daily air temperature range (22.24 °C). Analysis of thermistor data suggests that a linear thermal gradient of the debris layer can be assumed when the model is run at daily resolution. Modelled debris thickness, surface temperatures and melt capture patterns observed on the field, including the decrease in debris thickness and increasing melt with elevation and their variability across the glacier surface. The root mean square error between measured and observed melt is 0.16 m, corresponding to 22% of the average observed melt, in line with the performance of empirical models for debris-free and debris-covered ice. The model could thus be used to provide a first estimate of debris-covered ice melt for glaciers in the Adamello region. Improvements to the model would require measuring all energy fluxes on the glacier from a weather station and investigating their spatial distribution on the glacier surface.

Climate and Surface Mass Balance at Glaciar Perito Moreno, Southern Patagonia

Abstract The mass budget of southern Patagonian glaciers is characterized by an extreme amount of surface ablation. To understand the processes controlling surface mass balance, we analyzed in situ data including meteorological variables and ablation stakes for the 25 years between 1996 and 2020 near the terminus of Glaciar Perito Moreno in southern Patagonia in South America. The mean annual temperature has increased over the study period at a rate of 0.2°C decade−1. An energy-balance model was applied to calculate a point surface mass balance, based on meteorological records. The average point surface mass balance is estimated to be −16.3 m water equivalent (w.e.) yr−1 between 1996 and 2020, decreasing at a rate in the range from −0.4 to −0.9 m w.e. yr−1 decade−1. The greatest contribution to the surface energy balance was due to the sensible heat flux, and its variation drove the surface mass balance variation. The meteorological and surface mass balance records were compared with the Southern Annular Mode and El Niño–Southern Oscillation, which change the atmospheric circulation over southern Patagonia and influence surface mass balance near the terminus of the glacier. Our long-term dataset investigates the detailed meteorological conditions and surface mass balance and their connection with the large-scale climate variability over the last 25 years, reported for the first time in Patagonia.

Modelling Debris-Covered Glacier Ablation Using the Simultaneous Heat and Water Transport Model. Part 1: Model Development and Application to North Changri Nup

Modelling ablation of glacier ice under a layer of mineral debris is increasingly important, because the extent of supraglacial debris is expanding worldwide due to glacier recession. Physically based models have been developed, but the uncertainty in predictions is not yet well constrained. A new one-dimensional model of debris-covered ice ablation that is based on the Simultaneous Heat and Water transfer model is introduced here. SHAW-Glacier is a physically based, vertically integrated, fully coupled, water and energy balance model, which includes the advection of heat by rainwater and lateral flow. SHAW-Glacier was applied to North Changri Nup, a high elevation alpine glacier in the monsoon-dominated Central Himalaya. Simulations were compared with observed debris temperature profiles, snow depth, and ablation stake measurements for debris 0.03–0.41 m thick, in a 2500 m2 study area. Prediction uncertainty was estimated in a Monte Carlo analysis. SHAW-Glacier simulated the characteristic pattern of decreasing ablation with increasing debris thickness. However, the observations of ablation did not follow the characteristic pattern; annual ablation was highest where the debris was thickest. Recursive partitioning revealed a substantial, non-linear sensitivity to the snow threshold air temperature, suggesting a sensitivity to the duration of snow cover. Photographs showed patches of snow persisting through the ablation season, and the observational data were consistent with uneven persistence of snow patches. The analyses indicate that patchy snow cover in the ablation season can overwhelm the sensitivity of sub-debris ablation to debris thickness. Patchy snow cover may be an unquantified source of uncertainty in predictions of sub-debris ablation.

The Experimental Research of Ablation of Adishi Glacier Tongue

he study of glacier mass balance is gaining significant interest worldwide against modern global climate change. Using Drone gave good results in glacier mass balance studies. There is a sizeable difference between the values ​​obtained by the modern and previous used methods, which is primarily because of the shortcomings of the old one. We believe the results obtained cause the revision of the old data of the mass balance ablation component. Therefore, it became essential to restore the observation series and determine the compatibility of the data obtained by the old, traditional methods with the data obtained using new, modern technologies. The study aimed to compare the results obtained by standard ablation studies (ablation stake) and modern aerial photogrammetric methods based on actual observations in the Adishi glacier tongue.

The implications of permafrost thaw and land cover change on snow water equivalent accumulation, melt and runoff in discontinuous permafrost peatlands

AbstractIn the discontinuous permafrost zone of the Northwest Territories (NWT), Canada, snow covers the ground surface for half the year. Snowmelt constitutes a primary source of moisture supply for the short growing season and strongly influences stream hydrographs. Permafrost thaw has changed the landscape by increasing the proportional coverage of permafrost‐free wetlands at the expense of permafrost‐cored peat plateau forests. The biophysical characteristics of each feature affect snow water equivalent (SWE) accumulation and melt rates. In headwater streams in the southern Dehcho region of the NWT, snowmelt runoff has significantly increased over the past 50 years, despite no significant change in annual SWE. At the Fort Simpson A climate station, we found that SWE measurements made by Environment and Climate Change Canada using a Nipher precipitation gauge were more accurate than the Adjusted and Homogenized Canadian Climate Dataset which was derived from snow depth measurements. Here, we: (a) provide 13 years of snow survey data to demonstrate differences in end‐of‐season SWE between wetlands and plateau forests; (b) provide ablation stake and radiation measurements to document differences in snow melt patterns among wetlands, plateau forests, and upland forests; and (c) evaluate the potential impact of permafrost‐thaw induced wetland expansion on SWE accumulation, melt, and runoff. We found that plateaus retain significantly (p < 0.01) more SWE than wetlands. However, the differences are too small (123 mm and 111 mm, respectively) to cause any substantial change in basin SWE. During the snowmelt period in 2015, wetlands were the first feature to become snow‐free in mid‐April, followed by plateau forests (7 days after wetlands) and upland forests (18 days after wetlands). A transition to a higher percentage cover of wetlands may lead to more rapid snowmelt and provide a more hydrologically‐connected landscape, a plausible mechanism driving the observed increase in spring freshet runoff.

Quantifying heterogeneous monsoonal melt on a debris-covered glacier in Nepal Himalaya using repeat uncrewed aerial system (UAS) photogrammetry

AbstractThe ablation zones of debris-covered glaciers in Himalaya exhibit heterogeneous processes and melt patterns. Although sub-debris melt is measured at ablation stakes, the high variability of debris thickness necessitates distributed melt measurements at the glacier scale. Focusing on Annapurna III Glacier, we used uncrewed aerial system (UAS) photogrammetry to estimate total volume loss and slope-perpendicular glacier melt between May and November 2019 using flow-corrected point clouds. Results indicated the average elevation change was −1.10 ± 0.19 m, while the mean melt was −0.87 m w.e., equating to a mean melt rate of −0.47 cm w.e. d−1. However, the spatial pattern was highly variable due to complex local processes necessitating future study over short intervals. The evaluation of specific areas showed the interplay of debris thickness variability, subseasonal debris redistribution, supraglacial channel reconfiguration and the imprint of relict ice cliffs in leading to contemporary melt rates. Ice cliffs had higher melt distances (mean −3.9 ± 0.19 m) compared to non-cliff areas (mean −0.75 ± 0.19 m) and were the predominant control on the spatial patterns of seasonal melt rates. Crucially, the definition of ice cliff areas from thinning data has a profound impact on derived melt rates and melt enhancement. Our study demonstrates the possibility and utility of deriving fully-distributed slope-perpendicular melt measurements.

Measurement of unsaturated meltwater percolation flux in seasonal snowpack using self-potential

AbstractThis paper presents a feasibility study of in situ field measurements of unsaturated meltwater percolation flux within the vertical profile of a snowpack, using the self-potential (SP) method. On-site snowmelt column tests calibrated the SP measurements. The SP data measured electrical field strength with an electrode spacing of 20 cm, and coincident water saturation (Sw) measurements using time domain reflectometry allowed calculation of SP-modeled vertical percolation flux (qsp), expressed as Darcy velocity. The results reflected transient diurnal snowmelt dynamics, with peak flux lagging arrival of a saturation wetting front. Peak dailyqspwas 60 to >300 mm d−1, whereas daily snowmelt was 20–50 mm w.e. Surface refreezing events appeared to cause upward flow, possibly representing water redistribution toward the freezing boundary. Calculated fluxes were comparable to actual fluxes, although average errors ranged from −15 to +46% compared to average of melt expected from surface energy-balance and ablation stake measurements. By advancing method development to measure unsaturated meltwater percolation flux in snowpacks this study creates opportunities to study fundamental snowmelt processes, may improve mathematical modeling and may supplement glacier mass-balance studies and studies of snowmelt interactions with avalanches, groundwater and surface water.

Surface Velocity Analysis of Surge Region of Karayaylak Glacier from 2014 to 2020 in the Pamir Plateau

The west branch of Karayaylak Glacier (eastern Pamir Plateau) surged in May 2015, significantly impacting on local socio-economic development. This event was also of great significance for studies of surging glaciers. Using Sentinel-1 imagery analyzed by offset tracking, based on normalized cross-correlation (NCC), and with the support of the Google Earth Engine (GEE) platform, we quantified the ice surface velocity of the west branch and terminus of Karayaylak Glacier from 13 October 2014 to 17 October 2020. Sentinel-1 images were acquired at intervals of 12 or 24 days. We also used a three-dimensional (3-D) laser scanner to measure the velocity of 3 ablation stakes and 56 feature points in the study region from 15 August to 6 October 2015, for the purpose of accuracy assessment. We set up an automatic meteorological station to record the air temperature in the same period and combined this with data from Tashkurgan meteorological station from 1957 to 2015. Analysis of this dataset provided insights into the glacier surge mechanism, with the following conclusions. (1) Surface velocity of the west branch and terminus of Karayaylak Glacier increased sharply after October 2014. The velocity then dropped significantly in the two months after the surge, and stayed at low values for nearly a year. After 2017, the velocity was slightly higher than in the previous period. (2) The surge event occurred from 11 April to 17 May 2015; the average surface velocity in this phase attained 2395 m a−1 with a maximum velocity of 4265 m a−1 at the west branch terminus. (3) From 2017 to 2020, the velocity showed periodic annual changes. (4) Based on the meteorological data analysis, we conclude that this surge resulted from the interaction between thermal and hydrological control mechanisms. Simultaneously, we demonstrate the high potential of the GEE platform and Sentinel-1 data to extract glacier surface velocity.

Estimation of the recession rate of Gangotri glacier, Garhwal Himalaya (India) through kinematic GPS survey and satellite data

Snout monitoring of the Gangotri glacier (Uttarakhand, India) during the ablation season (May to September) in years 2005 and 2015 by using rapid static and kinematic GPS survey reveals that the retreating rate has been comparatively more declined than shown by the earlier studies. Our study is based on the individual measurement by the remote sensing, added by the ground observations by using Differential global positioning system (DGPS) to determine the precise recession rate of the glacier at sub-centimeter level of accuracy. The GPS dataset show that the total average retreat along the snout has been 102.57 ± 0.05 m from 2005 to 2015 with an average rate as 10.26 ± 0.05 m/yr. Additionally, the shift in snout position was also measured through multi-temporal satellite data from 1989 to 2016. The results indicate that the Gangotri glacier snout has retreated by 585.62 ± 38.30 m during this period with an average retreat of 26.75 ± 4.36 m/yr from 1989 to 1999, 21.58 ± 3.77 m/yr from 1999 to 2009 and 14.60 ± 4.81 m/yr from 2009 to 2016. Such a decline in retreat is further confirmed by the satellite data set. A close examination of meltwater discharge and retreating rate (r2 = 0.95) show that both parameters are strongly correlated. Therefore, we suggest that a consistent decrease in meltwater discharge from 1999 to 2015 is in agreement with decreasing trend of retreating rate during the recent years. To determine the possible causes of decreased retreating rate, a relationship between debris thickness and melt rate was also established by ablation stakes. Further, we infer that the declining trend in the glacier retreat is not only controlled by prevailing weather conditions (rainfall and air temperature) but is also governed by increased debris cover on the glacier surface which prevents the ice to melt.

Processes at the margins of supraglacial debris cover: Quantifying dirty ice ablation and debris redistribution

AbstractCurrent glacier ablation models have difficulty simulating the high‐melt transition zone between clean and debris‐covered ice. In this zone, thin debris cover is thought to increase ablation compared to clean ice, but often this cover is patchy rather than continuous. There is a need to understand ablation and debris dynamics in this transition zone to improve the accuracy of ablation models and the predictions of future debris cover extent. To quantify the ablation of partially debris‐covered ice (or ‘dirty ice’), a high‐resolution, spatially continuous ablation map was created from repeat unmanned aerial systems surveys, corrected for glacier flow in a novel way using on‐glacier ablation stakes. Surprisingly, ablation is similar (range ~5 mm w.e. per day) across a wide range of percentage debris covers (~30–80%) due to the opposing effects of a positive correlation between percentage debris cover and clast size, countered by a negative correlation with albedo. Once debris cover becomes continuous, ablation is significantly reduced (by 61.6% compared to a partial debris cover), and there is some evidence that the cleanest ice (<~15% debris cover) has a lower ablation than dirty ice (by 3.7%). High‐resolution feature tracking of clast movement revealed a strong modal clast velocity where debris was continuous, indicating that debris moves by creep down moraine slopes, in turn promoting debris cover growth at the slope toe. However, not all slope margins gain debris due to the removal of clasts by supraglacial streams. Clast velocities in the dirty ice area were twice as fast as clasts within the continuously debris‐covered area, as clasts moved by sliding off their boulder tables. These new quantitative insights into the interplay between debris cover characteristics and ablation can be used to improve the treatment of dirty ice in ablation models, in turn improving estimates of glacial meltwater production. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd

Surface melt and the importance of water flow – an analysis based on high-resolution unmanned aerial vehicle (UAV) data for an Arctic glacier

Abstract. Models of glacier surface melt are commonly used in studies of glacier mass balance and runoff; however, with limited data available, most models are validated based on ablation stakes and data from automatic weather stations (AWSs). The technological advances of unmanned aerial vehicles (UAVs) and structure from motion (SfM) have made it possible to measure glacier surface melt in detail over larger portions of a glacier. In this study, we use melt measured using SfM processing of UAV imagery to assess the performance of an energy balance (EB) and enhanced temperature index (ETI) melt model in two dimensions. Imagery collected over a portion of the ablation zone of Fountain Glacier, Nunavut, on 21, 23, and 24 July 2016 was previously used to determine distributed surface melt. An AWS on the glacier provides some measured inputs for both models as well as an additional check on model performance. Modelled incoming solar radiation and albedo derived from UAV imagery are also used as inputs for both models, which were used to estimate melt from 21 to 24 July 2016. Both models estimate total melt at the AWS within 16 % of observations (4 % for ETI). Across the study area the median model error, calculated as the difference between modelled and measured melt (EB = −0.064 m, ETI = −0.050 m), is within the uncertainty of the measurements. The errors in both models were strongly correlated to the density of water flow features on the glacier surface. The relation between water flow and model error suggests that energy from surface water flow contributes significantly to surface melt on Fountain Glacier. Deep surface streams with highly asymmetrical banks are observed on Fountain Glacier, but the processes leading to their formation are missing in the model assessed here. The failure of the model to capture flow-induced melt would lead to significant underestimation of surface melt should the model be used to project future change.

ОЦЕНКА БАЛАНСА МАССЫ ЛЕДНИКА АЛЬДЕГОНДА (ЗАПАДНЫЙ ШПИЦБЕРГЕН) В 2015-2018 ГГ. НА ОСНОВЕ МОДЕЛИ ARCTICDEM, ГЕОДЕЗИЧЕСКИХ И ГЛЯЦИОЛОГИЧЕСКИХ ДАННЫХ

The Aldegonda (Aldegondabreen) Glacier, located on the Nordenskjold Land, West Spitsbergen, covers the area of about 6 km2 (in 2018) and does constantly retreat since the very first observations of 1936. In August 2018, a topographic survey of the glacier was carried out. By comparing the results with the ArcticDEM model, built from space images of 2015, the difference in heights of the surface over three years had been calculated. Comparison of this difference with in situ data of glaciological measurements by the ablation stakes, made during the same period 2015–2018, demonstrated a high correlation between them. Considering the almost complete absence of snow cover on the glacier at the end of the summer season, the difference was recalculated into the spatial distribution of the specific mass balance by multiplying the ice density (0.88 g cm−3). Using the empirical dependence of the specific mass balance on the altitude above sea level, the obtained values were extrapolated to that part of the glacier which was not surveyed in 2018. The total loss of the Aldegonda Clacier mass for 2015–2018, calculated on the basis of topographic survey and the ArcticDEM, was estimated as 30.3 million tons (about 10.1 million tons per year). This magnitude gives the value of mean annual specific balance of approximately −1.76 m w.e, which is almost 2.5 times larger modulo than the previously published mean for the period 1936–1990, but close to the values of the early 2000s. Despite the small difference in the values obtained by geodetic and glaciological methods, the measurements does not show a systematic shift relative to each other and demonstrate approximately the same intervals of specific balance from the glacier tongue to its upper reaches (−1.08 ÷ −3.01 m w.e). This makes possible to conclude that the ArcticDEM model has the satisfactory vertical accuracy (both relative and absolute) to study on its basis changes in the surface height of an individual glacier.

Automatic measurement of glacier ice ablation using thermistor strings

AbstractIn this work we tested the suitability of thermistor strings as automatic tools for the continuous measurement of glacier ice ablation. Experimental data collected in summer 2017 over an Italian glacier provided ice ablation readings with accuracy similar to manual measurements with ablation stakes and other automatic systems, like the draw-wire method and the Bøggild ablatometer. Thermistor strings have potential for future applications in remote glacier monitoring, thanks to their flexibility, simple construction, and robustness.

Detecting Short-Term Surface Melt on an Arctic Glacier Using UAV Surveys

Current understanding of glacier mass balance changes under changing climate is limited by scarcity of in situ measurements in both time and space, as well as resolution of remote sensing products. Recent innovations in unmanned aerial vehicles (UAVs), as well as structure-from-motion photogrammetry (SfM), have led to increased use of digital imagery to derive topographic data in great detail in many fields, including glaciology. This study tested the capability of UAV surveys to detect surface changes over glacier ice during a three-day period in July 2016. Three UAV imaging missions were conducted during this time over 0.185 km 2 of the ablation area of Fountain Glacier, NU. These were processed with the SfM algorithms in Agisoft Photoscan Professional and overall accuracies of the resulting point clouds ranged from 0.030 to 0.043 m. The high accuracy of point clouds achieved here is primarily a result of a small ground sampling distance (0.018 m), and is also influenced by GPS precision. Glacier surface change was measured through differencing of point clouds and change was compared to ablation stake measurements. Surface change measured with the UAV-SfM method agreed with the coincident ablation stake measurements in most instances, with RMSE values of 0.033, 0.028, and 0.042 m for one-, two-, and three-day periods, respectively. Total specific melt over the study area measured with the UAV was 0.170 m water equivalent (w.e.), while interpolation of ablation measurements resulted in 0.144 m w.e. Using UAVs to measure small changes in glacier surfaces will allow for new investigations of distribution of mass balance measurements.

Assessing glacier melt contribution to streamflow at Universidad Glacier, central Andes of Chile

Abstract. Glacier melt is an important source of water for high Andean rivers in central Chile, especially in dry years, when it can be an important contributor to flows during late summer and autumn. However, few studies have quantified glacier melt contribution to streamflow in this region. To address this shortcoming, we present an analysis of meteorological conditions and ablation for Universidad Glacier, one of the largest valley glaciers in the central Andes of Chile at the head of the Tinguiririca River, for the 2009–2010 ablation season. We used meteorological measurements from two automatic weather stations installed on the glacier to drive a distributed temperature-index and runoff routing model. The temperature-index model was calibrated at the lower weather station site and showed good agreement with melt estimates from an ablation stake and sonic ranger, and with a physically based energy balance model. Total modelled glacier melt is compared with river flow measurements at three sites located between 0.5 and 50 km downstream. Universidad Glacier shows extremely high melt rates over the ablation season which may exceed 10 m water equivalent in the lower ablation area, representing between 10 and 13 % of the mean monthly streamflow at the outlet of the Tinguiririca River Basin between December 2009 and March 2010. This contribution rises to a monthly maximum of almost 20 % in March 2010, demonstrating the importance of glacier runoff to streamflow, particularly in dry years such as 2009–2010. The temperature-index approach benefits from the availability of on-glacier meteorological data, enabling the calculation of the local hourly variable lapse rate, and is suited to high melt regimes, but would not be easily applicable to glaciers further north in Chile where sublimation is more significant.

An enhanced temperature index model for debris-covered glaciers accounting for thickness effect

Debris-covered glaciers are increasingly studied because it is assumed that debris cover extent and thickness could increase in a warming climate, with more regular rockfalls from the surrounding slopes and more englacial melt-out material. Debris energy-balance models have been developed to account for the melt rate enhancement/reduction due to a thin/thick debris layer, respectively. However, such models require a large amount of input data that are not often available, especially in remote mountain areas such as the Himalaya, and can be difficult to extrapolate. Due to their lower data requirements, empirical models have been used extensively in clean glacier melt modelling. For debris-covered glaciers, however, they generally simplify the debris effect by using a single melt-reduction factor which does not account for the influence of varying debris thickness on melt and prescribe a constant reduction for the entire melt across a glacier.In this paper, we present a new temperature-index model that accounts for debris thickness in the computation of melt rates at the debris-ice interface. The model empirical parameters are optimized at the point scale for varying debris thicknesses against melt rates simulated by a physically-based debris energy balance model. The latter is validated against ablation stake readings and surface temperature measurements. Each parameter is then related to a plausible set of debris thickness values to provide a general and transferable parameterization. We develop the model on Miage Glacier, Italy, and then test its transferability on Haut Glacier d’Arolla, Switzerland.The performance of the new debris temperature-index (DETI) model in simulating the glacier melt rate at the point scale is comparable to the one of the physically based approach, and the definition of model parameters as a function of debris thickness allows the simulation of the nonlinear relationship of melt rate to debris thickness, summarised by the Østrem curve. Its large number of parameters might be a limitation, but we show that the model is transferable in time and space to a second glacier with little loss of performance. We thus suggest that the new DETI model can be included in continuous mass balance models of debris-covered glaciers, because of its limited data requirements. As such, we expect its application to lead to an improvement in simulations of the debris-covered glacier response to climate in comparison with models that simply recalibrate empirical parameters to prescribe a constant across glacier reduction in melt.