Supraglacial debris cover regulates the melt rates of many glaciers in mountainous regions around the world, thereby modifying the availability and quality of downstream water resources. However, the influence of supraglacial debris is often poorly represented within glaciological models, due to the absence of a technique to provide high-precision, spatially continuous measurements of debris thickness. Here, we use high-resolution UAV-derived thermal imagery, in conjunction with local meteorological data, visible UAV imagery and vertically profiled debris temperature time series, to model the spatially distributed debris thickness across a portion of Llaca Glacier in the Cordillera Blanca of Peru. Based on our results, we simulate daily sub-debris melt rates over a 3-month period during 2019. We demonstrate that, by effectively calibrating the radiometric thermal imagery and accounting for temporal and spatial variations in meteorological variables during UAV surveys, thermal UAV data can be used to more precisely represent the highly heterogeneous patterns of debris thickness and sub-debris melt on debris-covered glaciers. Additionally, our results indicate a mean sub-debris melt rate nearly three times greater than the mean melt rate simulated from satellite-derived debris thicknesses, emphasising the importance of acquiring further high-precision debris thickness data for the purposes of investigating glacier-scale melt processes, calibrating regional melt models and improving the accuracy of runoff predictions.

Subglacial hydrological systems require innovative technological solutions to access and observe. Wireless sensor platforms can be used to collect and return data, but their performance in deep and fast-moving ice requires quantification. We report experimental results from Cryoegg: a spherical probe that can be deployed into a borehole or moulin and transit through the subglacial hydrological system. The probe measures temperature, pressure and electrical conductivity in situ and returns all data wirelessly via a radio link. We demonstrate Cryoegg's utility in studying englacial channels and moulins, including in situ salt dilution gauging. Cryoegg uses VHF radio to transmit data to a surface receiving array. We demonstrate transmission through up to 1.3 km of cold ice – a significant improvement on the previous design. The wireless transmission uses Wireless M-Bus on 169 MHz; we present a simple radio link budget model for its performance in cold ice and experimentally confirm its validity. Cryoegg has also been tested successfully in temperate ice. The battery capacity should allow measurements to be made every 2 h for more than a year. Future iterations of the radio system will enable Cryoegg to transmit data through up to 2.5 km of ice.

The detection and monitoring of meltwater within firn presents a significant monitoring challenge. We explore the potential of small wireless sensors (ETracer+, ET+) to measure temperature, pressure, electrical conductivity and thus the presence or absence of meltwater within firn, through tests in the dry snow zone at the East Greenland Ice Core Project site. The tested sensor platforms are small, robust and low cost, and communicate data via a VHF radio link to surface receivers. The sensors were deployed in low-temperature firn at the centre and shear margins of an ice stream for 4 weeks, and a ‘bucket experiment’ was used to test the detection of water within otherwise dry firn. The tests showed the ET+ could log subsurface temperatures and transmit the recorded data through up to 150 m dry firn. Two VHF receivers were tested: an autonomous phase-sensitive radio-echo sounder (ApRES) and a WinRadio. The ApRES can combine high-resolution imaging of the firn layers (by radio-echo sounding) with in situ measurements from the sensors, to build up a high spatial and temporal resolution picture of the subsurface. These results indicate that wireless sensors have great potential for long-term monitoring of firn processes.

Aquatic dissolved organic matter (DOM) is a major reservoir of reduced organic carbon and has a significant influence on heterotrophic biological productivity and water quality in marine and freshwater environments. Although the forms and transformations of DOM in temperate aquatic and soil environments have been studied extensively, this is not the case for glacial environments. In this study, fluorescent excitation–emission matrices (EEMs), parallel factor analysis (PARAFAC) and cluster analysis were used to characterize the fluorescing components of DOM in ice and water samples from supraglacial, englacial, subglacial and proglacial environments of seven glaciers in the Canadian Arctic, Norway and Antarctica. At least five significant fluorescent DOM fractions were identified, which accounted for 98.2% of the variance in the dataset. These included four protein-like components and one humic-like component. The predominantly proteinaceous character of DOM from these glaciers is very different from the more humic character of DOM described previously from lacustrine, fluvial, estuarine and marine environments. DOM from the sampled glaciers is broadly similar in character despite their geographically distinct locations, different thermal regimes and inter- and intra-site differences in potential organic matter sources. Glacier ice samples had a relatively low ratio of humic-like :protein-like fluorescence while meltwater samples had a higher ratio.

Glacier mass balance and hydrology are strongly influenced by the distribution of snow accumulation at the start of the melt season. Two successive end-of-winter snow-cover surveys at Finsterwalderbreen, Svalbard, are here used to investigate the interannual variability in the spatial distribution of accumulation, and its relationship with topography. 40–62% of the variance in snow depth was not determined by elevation (assessed by linear regression of snow depth on surface elevation), which could not therefore necessarily be used as a sole predictor of the spatial distribution of accumulation here. Principal components (PC) analysis of the topographic variables elevation, slope, north–south and east–west aspects shows that only two of six PCs, determined for 2years’ sampling locations, had maximum loadings on altitude; aspect was more important, with maximum loadings on four PCs. Hierarchical cluster analysis was then applied to these PCs: significant correlations with accumulation in each of two terrain clusters were given by (1) elevation and slope, (2) east–west aspect only (1999); (1) elevation only, (2) no significant correlations (2000). There is strong interannual variability not only in the magnitude of winter accumulation (0.41 mw.e. in 1999, 0.58 mw.e. in 2000), but also in its spatial distribution, and its relationship with topography.

Quantifying the biogeochemical cycling of carbon in glacial ecosystems is of great significance for regional, and potentially global, carbon flow estimations. The concentration and quality of organic carbon (OC) is an important indicator of biogeochemical and physical processes that prevail in an ice-sheet ecosystem. Here we determine the content and quality of OC in debris from the surface of the Greenland ice sheet (GrIS) using microscopic, chromatographic, spectrophotometric and high-temperature combustion techniques. The total OC content in the debris increased with distance from the edge of the ice sheet, from virtually zero to >6% dry weight at 50 km inland, and there was a peak in the carbohydrate proportion and the microbial abundance at ∼6km inland. The highest (galactose + mannose)/(arabinose + xylose) ratios, indicating maximum autochthonous microbial production, were found at >10km inland. We propose that three key processes influence the carbon cycling on the GrIS: aeolian input of microbial inoculum and nutrients, in situ biological C transformation and the wash-away of supraglacial debris by meltwaters. We show that all these processes have significant spatial variability. While the total OC content of the debris on the ice sheet is probably controlled by the physical processes of wind transport and wash-away by meltwater, the microbial abundance and the quantity of the labile cell-contained OC within the debris is likely to be driven by the balance between the wash-away and the microbial productivity.

A surface-energy/mass-balance model with an explicit calculation of meltwater refreezing and superimposed ice formation is applied to midre Lovénbreen, Spitsbergen, Svalbard. The model is run with meteorological measurements to represent the present climate, and run with scenarios taken from global climate model predictions based on the IS92a emissions scenario to represent future climates. Model results indicate that superimposed ice accounts for on average 37% of the total net accumulation under present conditions. The model is found to be highly sensitive to changes in the mean annual air temperature and much less sensitive to changes in the total annual precipitation. A 0.5˚C decade–1 temperature increase is predicted to cause an average mass-balance change of –0.43 ma–1, while a 2% decade–1 increase in precipitation will result in only a +0.02 ma–1 change in mass balance. An increase in temperature results in a significant decrease in the size of the accumulation area at midre Lovénbreen and hence a similar decrease in the net volume of superimposed ice. The model predicts, however, that the relative importance of superimposed ice will increase to account for >50% of the total accumulation by 2050. The results show that the refreezing of meltwater and in particular the formation of superimposed ice make an important positive contribution to the mass balance of midre Lovénbreen under present conditions and will play a vital future role in slowing down the response of glacier mass balance to climate change.

Superimposed ice forms when meltwater refreezes onto a sub-freezing glacier surface. The accumulation zones of many Arctic glaciers include large areas of superimposed ice, which for mass-balance purposes have to be distinguished from the ablation zone consisting of glacier ice. We examine the ability of synthetic aperture radar (SAR) satellite sensors to detect superimposed ice on the glaciers Kongsvegen and midre Lovénbreen on Svalbard. Structural analysis of ice cores as well as surface observations from these glaciers in 1999 and 2000 provide a spatial record of superimposed ice. Winter SAR images show three distinct zones, which correspond closely to areas of glacier ice, superimposed ice and firn. This is seen very clearly on Kongsvegen, but not as clearly on the much smaller midre Lovénbreen. One possible explanation for the contrasting SAR signal may relate to the differing air-bubble content of firn, superimposed ice and glacier ice. Thin layers of winter-formed superimposed ice (510 cm) in some places are not seen on the SAR images, indicating that a certain thickness is needed for detection. The equilibrium-line altitude cannot be detected since the SAR cannot differentiate old superimposed ice, superimposed ice formed currently in the accumulation area in summer and superimposed ice formed currently in the ablation zone in autumn and winter.

Interpretation of ice mass elevation changes observed by satellite altimetry demands quantification of the proportion of elevation change which is attributable to variations in firn densification. Detailed stratigraphic logging of snowpack structure and density was carried out at ~1km intervals along a 47 km transect on Devon Ice Cap, Canada, in spring (pre-melt) and autumn (during/ after melt) 2004 and 2006 to characterize seasonal snowpack variability across the full range of snow facies. Simultaneous meteorological measurements were gathered. Spring (pre-melt) snowpacks show low variability over large spatial scales, with low-magnitude changes in density. The end-of-summer/ autumn density profiles show high variability in both 2004 and 2006, with vastly different melt regimes generating dissimilar patterns of ice-layer formation over the two melt seasons. Dye-tracing experiments from spring to autumn 2006 reveal that vertical and horizontal distribution of meltwater flow within and below the annual snowpack is strongly affected by the pre-existing, often subtle stratigraphic interfaces in the snowpack, rather than its bulk properties. Strong interannual variability suggests that using a simple relationship between air temperature, elevation and snowpack densification to derive mass change from measurements of elevation change across High Arctic ice caps may be misguided. Melt timing and duration are important extrinsic factors governing snowpack densification and ice-layer formation in summer, rather than averaged air temperatures.

Meteorological and snow-temperature data from midtre Lovénbreen, Spitsbergen, Svalbard, indicate two distinct annual phases of rapid snowpack warming and superimposed-ice formation in 1998/99. Short periods of positive air temperatures in early winter, lasting up to 36 hours and often coinciding with rainfall, caused rapid glacier-wide melting. Percolating water froze to form superimposed ice on the lower half of the glacier, and wetted–refrozen snow and ice lenses at higher altitudes. The second period of superimposed-ice formation commenced in May/June 1999 and continued for 5 weeks at low altitudes and throughout the summer at high altitudes. These observations at midtre Lovénbreen are typical of Spitsbergen glaciers and reflect the unique climatology of the region. They contrast with those from glaciers in more continental climatic settings where superimposed-ice formation is confined to a single period during summer. There are significant implications for glacier mass balance, with superimposed ice locally comprising up to 20% of winter balances and accounting for ∼16–25% of the annual accumulation. Since projected climatic warming is greatest during the winter months in Arctic regions, superimposed ice may become an increasingly important component of the winter, and potentially the net, balance of Spitsbergen glaciers.

The environments underneath ice sheets are of high scientific interest. Wireless sensors offer the prospect of sustained, distributed remote sensing in the subglacial environment. Typically, wireless sensor networks use radio-frequency (RF) electromagnetic communications, but these are highly attenuated in wet environments. In such environments, acoustic communications may be more power-efficient. Here we review the literature on acoustic and RF attenuation through ice and other relevant media, and present the results of new experiments on acoustic attenuation in glacial ice. Link budgets for communications from a range of subglacial environments show that acoustic communications are a viable strategy for transmission through water and ice where RF is too highly attenuated to be detected. Acoustic communication at 30 kHz is predicted to be possible through 1 km of glacial ice, using a 1 W transmitter. Such a strategy may be appropriate for shallow ice-stream environments around the Antarctic and Greenland ice sheet margins.

The ability to perform microbial detection and characterization in-field at extreme environments, rather than on returned samples, has the potential to improve the efficiency, relevance and quantity of data from field campaigns. To date, few examples of this approach have been reported. Therefore, we demonstrate that the approach is feasible in subglacial environments by deploying four techniques for microbial detection: real-time polymerase chain reaction; microscopic fluorescence cell counts, adenosine triphosphate bioluminescence assay and recombinant Factor C assay (to detect lipopolysaccharide). Each technique was applied to 12 subglacial ice samples, 12 meltwater samples and two snow samples from Engabreen, Northern Norway. Using this multi-technique approach, the detected biomarker levels were as expected, being highest in debris-rich subglacial ice, moderate in glacial meltwater and low in clean ice (debris-poor) and snow. Principal component analysis was applied to the resulting dataset and could be performed in-field to rapidly aid the allocation of resources for further sample analysis. We anticipate that in-field data collection will allow for multiple rounds of sampling, analysis, interpretation and refinement within a single field campaign, resulting in the collection of larger and more appropriate datasets, ultimately with more efficient science return.