We developed automated ablation stakes to measure colocated in situ changes in relative glacier-surface elevation and climatological drivers of ablation. The designs, refined over 10 years of development and deployments, implement open-source hardware and common building materials. The ablation stakes record distance to the snow/ice surface, air temperature and relative humidity every 1–15 min. Using these high-frequency data, we demonstrate that melt factors calculated using standard melt-rate vs temperature regressions converge over averaging windows of approximately 12 h or greater. Furthermore, we evaluate an integral approach to estimating temperature-index melt factors for ablation. In a test case on Glaciar Perito Moreno, Argentina, this integral approach reveals an overall positive-degree-day melt factor of 7.5 mm w.e. $^\circ$C−1 d−1. We describe four deployments with iteratively improved designs and provide a list of materials required to construct an automated ablation stake.

Observations of glacier melt and runoff are of fundamental interest in the study of glaciers and their interactions with their environment. Considerable recent interest has developed around distributed acoustic sensing (DAS), a sensing technique which utilizes Rayleigh backscatter in fiber optic cables to measure the seismo-acoustic wavefield in high spatial and temporal resolution. Here, we present data from a month-long, 9 km DAS deployment extending through the ablation and accumulation zones on Rhonegletscher, Switzerland, during the 2020 melt season. While testing several types of machine learning (ML) models, we establish a regression problem, using the DAS data as the dependent variable, to infer the glacier discharge observed at a proglacial stream gauge. We also compare two predictive models that only depend on meteorological station data. We find that the seismo-acoustic wavefield recorded by DAS can be utilized to infer proglacial discharge. Models using DAS data outperform the two models trained on meteorological data with mean absolute errors of 0.64, 2.25 and 2.72 m3 s−1, respectively. This study demonstrates the ability of in situ glacier DAS to be used for quantifying proglacial discharge and points the way to a new approach to measuring glacier runoff.

Glaciers in the Southern Patagonia Icefield (SPI) have been shrinking in recent decades, but due to a lack of field observations, understanding of the drivers of ablation is limited. We present a distributed surface energy balance model, forced with meteorological observations from a west–east transect located in the north of the SPI. Between October 2015 and June 2016, humid and warm on-glacier conditions prevailed on the western side compared to dry and cold conditions on the eastern side. Controls of ablation differ along the transect, although at glacier-wide scale sensible heat (mean of 72 W m−2 to the west and 51 W m−2 to the east) and net shortwave radiation (mean of 54 W m−2 to the west and 52 W m−2 to the east) provided the main energy sources. Net longwave radiation was an energy sink, while latent heat was the most spatially variable flux, being an energy sink in the east (−4 W m−2) and a source in the west (20 W m−2). Ablation was high, but at comparable elevations, it was greater to the west. These results provide new insights into the spatial variability of energy-balance fluxes and their control over the ablation of Patagonian glaciers.

This 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 daily qsp was 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.

Processes governing meltwater penetration into cold firn remain poorly constrained. Here, in situ experiments are used to develop a grain-scale model to investigate physical limitations on meltwater infiltration in firn. At two sites in Greenland, drilling pumped water into cold firn to >75 m depth, and the thermo-hydrologic evolution of the firn column was measured. Rather than filling all available pore space, the water formed perched aquifers with downward penetration halted by thermal and density conditions. The two sites formed deep aquifers at ~40 m depth and at densities considerably less than the air pore close-off density (~725 kg m−3 at −18°C, and ~750 kg m−3 at −14°C), demonstrating that some pore space at depth remains inaccessible. A geometric grain-scale model of firn is constructed to quantify the limits of a descending fully saturated wetting front in cold firn. Agreement between the model and field data implies the model includes the first-order effects of water and heat flow in a firn lattice. The model constrains the relative importance of firn density, temperature and grain/pore size in inhibiting wetting front migration. Results imply that deep infiltration, including that which leads to firn aquifer formation, does not have access to all available firn pore space.

Relatively little is known about the physical mechanisms that drive the dynamics of the East Antarctic outlet glaciers. Here we conduct a remote-sensing investigation of the Polar Record Glacier (PRG), East Antarctica to analyze its ice flow acceleration, ice front variations and ice surface melting. Ice flow speeds at PRG increased by up to 15% from 2005 to 2015, with substantial interannual fluctuations. The ice velocities also showed seasonal variations, accelerating by up to 9% between September and January. Multiple mechanisms contribute to the observed seasonal variations: the initial acceleration may result from the lost back-stress provided by the sea ice in the austral spring and the later speedup relate to the surface meltwater that leads to weakened ice shelf and shear margins. The sensitivity of the PRG to oceanic forcing is confirmed by comparing the secular ice velocity increases with ocean temperatures. These measurements suggest that the dynamics of East Antarctic ice shelves are sensitive to melt at both the surface and base, at a range of timescales.

Changing Arctic sea-ice extent and melt season duration, and increasing economic interest in the Arctic have prompted the need for enhanced marine ecosystem studies and improvements to dynamical and forecast models. Sea-ice melt pond fraction fp has been shown to be correlated with the September minimum ice extent due to its impact on ice albedo and heat uptake. Ice forecasts should benefit from knowledge of fp as melt ponds form several months in advance of ice retreat. This study goes further back by examining the potential to predict fp during winter using backscatter data from the commonly available Sentinel-1 synthetic aperture radar. An object-based image analysis links the winter and spring thermodynamic states of first-year and multiyear sea-ice types. Strong correlations between winter backscatter and spring fp, detected from high-resolution visible to near infrared imagery, are observed, and models for the retrieval of fp from Sentinel-1 data are provided (r2 ≥ 0.72). The models utilize HH polarization channel backscatter that is routinely acquired over the Arctic from the two-satellite Sentinel-1 constellation mission, as well as other past, current and future SAR missions operating in the same C-band frequency. Predicted fp is generally representative of major ice types first-year ice and multiyear ice during the stage in seasonal melt pond evolution where fp is closely related to spatial variations in ice topography.