Fine resolution topographic data derived from methods such as Structure from Motion (SfM) and Multi-View Stereo (MVS) have the potential to provide detailed observations of geomorphological change, but have thus far been limited by the logistical constraints of conducting repeat surveys in the field. Here, we present the results from an automated time-lapse camera array, deployed around an ice-marginal lake on the western margin of the Greenland ice sheet. Fifteen cameras acquired imagery three-times per day over a 426 day period, yielding a dataset of ~19 000 images. From these data we derived 18 point clouds of the ice-margin across a range of seasons and successfully identified calving events (ranging from 234 to 1475 m2 in area and 815–8725 m3 in volume) induced by ice cliff undercutting at the waterline and the collapse of spalling flakes. Low ambient light levels, locally reflective surfaces and the large survey range hindered analysis of smaller scale ice-margin dynamics. Nevertheless, this study demonstrates that an integrated SfM-MVS and time-lapse approach can be employed to generate long-term 3-D topographic datasets and thus quantify ice-margin dynamics at a fine spatio-temporal scale. This approach provides a template for future studies of geomorphological change.

Wind-packed snow in the form of slabs or crusts is an important part of alpine and polar snow covers. Yet, the formation process of such layers is poorly understood. For example, it remains unclear whether drifting snow is necessary for wind-packing or not. A better understanding of wind-packing could improve snow-cover models and avalanche danger forecasts and contribute to the assessment of mass balances in polar regions. We designed a closed-circuit, obround wind tunnel to study the process of wind crust formation. A SnowMicroPen was used to measure how the hardness of the snow evolved. The results show that no crust forms without saltation. Drifting snow is a necessary but not sufficient condition for wind-packing. The dynamics of erosion and deposition appear to be equally important.

The Darwin–Hatherton Glacial system (DHGS) connects the East Antarctic Ice Sheet (EAIS) with the Ross Ice Shelf and is a key area for understanding past variations in ice thickness of surrounding ice masses. Here we present the first detailed measurements of ice thickness and grounding zone characteristics of the DHGS as well as new measurements of ice velocity. The results illustrate the changes that occur in glacier geometry and ice flux as ice flows from the polar plateau and into the Ross Ice Shelf. The ice discharge and the mean basal ice shelf melt for the first 8.5 km downstream of the grounding line amount to 0.24 ± 0.05 km3 a−1 and 0.3 ± 0.1 m a−1, respectively. As the ice begins to float, ice thickness decreases rapidly and basal terraces develop. Constructed maps of glacier geometry suggest that ice drainage from the EAIS into the Darwin Glacier occurs primarily through a deep subglacial canyon. By contrast, ice thins to <200 m at the head of the much slower flowing Hatherton Glacier. The glaciological field study establishes an improved basis for the interpretation of glacial drift sheets at the link between the EAIS and the Ross Ice Sheet.

The spatio-temporal distribution of air temperature over mountain glaciers can demonstrate complex patterns, yet it is often represented simplistically using linear vertical temperature gradients (VTGs) extrapolated from off-glacier locations. We analyse a network of centreline and lateral air temperature observations at Tsanteleina Glacier, Italy, during summer 2015. On average, VTGs are steep (<−0.0065 °C m−1), but they are shallow under warm ambient conditions when the correlation between air temperature and elevation becomes weaker. Published along-flowline temperature distribution methods explain centreline observations well, including warming on the lower glacier tongue, but cannot estimate lateral temperature variability. Application of temperature distribution methods improves simulation of melt rates (RMSE) in an energy-balance model by up to 36% compared to the environmental lapse rate extrapolated from an off-glacier station. However, results suggest that model parameters are not easily transferable to glaciers with a small fetch without recalibration. Such methods have potential to improve estimates of temperature across a glacier, but their parameter transferability should be further linked to the glacier and atmospheric characteristics. Furthermore, ‘cold spots’, which can be >2°C cooler than expected for their elevation, whose occurrence is not predicted by the temperature distribution models, are identified at one-quarter of the measurement sites.

The thickness of supraglacial debris cover controls how it impacts the ablation rate of underlying glacier ice, yet this quantity remains challenging to measure, particularly at glacier scales. We present a relatively straightforward, and cost-effective method to estimate debris thickness exposed above ice cliffs using simplified geometrical measurements from a high-resolution digital surface model (DSM), derived from a terrestrial photographic survey and a Structure from Motion with Multi-View Stereo workflow (SfM-MVS). As the ice surface relief beneath the debris cover is unknown, we assume it to be horizontal and provide error bounds based on characteristic ice-surface slope at the visible debris/ice interface. Debris thickness around the three sampled ice cliffs was highly variable (interquartile range of 0.80–2.85 m) and negatively skewed with a mean thickness of 2.08 ± 0.68 m. Manual, and high-frequency radar, determinations of debris thickness in the same area show similar thickness distributions, but statistically different mean debris thickness, due to local heterogeneity. Debris thickness values derived in this study all exceed estimates from satellite surface temperature inversions. Wider application of the method presented here would provide useful data for improving debris thickness approximations from satellite imagery.

The basal interface of glaciers is generally not directly observable. Geophysical inverse methods are therefore used to infer basal parameters from surface observations. Such methods can also provide information about potential inadequacies of the forward model. Ideally an inverse problem can be regularized so that the differences between modeled and observed surface velocities reflect observational errors. However, deficiencies in the forward model usually result in additional errors. Here we use the spatial pattern of velocity residuals to discuss the main error sources for basal stress inversions for Jakobshavn Isbræ, Greenland. Synthetic tests with prescribed patterns of basal yield stress with varying length scales are then used to investigate different weighting functions for the data-model misfit and for the ability of the inversion to resolve details in basal yield stress. We also test real-data inversions for their sensitivities to prior estimate, forward model parameters, data gaps, and temperature fields. We find that velocity errors are not sufficient to explain the residual patterns of real-data inversions. Conversely, ice-geometry errors and especially simulated errors in model simplifications are capable of reproducing similar error patterns and magnitudes. We suggest that residual patterns can provide useful guidance for forward model improvements.

We evaluate the variability in basal friction for Kronebreen, Svalbard, a fast-flowing tidewater glacier. We invert 3 years (2013–15) of surface velocities at high temporal resolution (generally 11 days), to estimate the changing basal properties of the glacier. Our results suggest that sliding behaviour of Kronebreen within a year is primarily influenced by changes in water input patterns during the meltwater season and basal friction is highly variable from a year to another. At present, models usually employ parameterisations to encompass the complex physics of glacier sliding by mathematically simulate their net effect. For such ice masses with strong seasonal variations of surface melt, the spatio-temporal patterns of basal friction imply that it is neither possible nor appropriate to use a parameterisation for bed friction that is fixed in space and/or time, at least in a timescale of a few years. Basal sliding may not only be governed by local processes such as basal topography or summer melt, but also be mediated by factors that vary over a larger distance and over a longer time period such as subglacial hydrology organisation, ice-thickness changes or calving front geometry.

Global climate change is causing widespread glacier retreat, with many small glaciers disappearing from the world's mountain ranges. We obtained the annual mass balance of a small glacier (the Ningchan No.1 Glacier) located on the northeastern Tibetan Plateau, from the years 2010 to 2015 using glaciological and geodetic methods. We also measured the glacier's thickness in 2014 using ground-penetrating radar. Employing topographical maps and ZY-3 images, we obtained Digital Elevation Models for 1972 and 2014. Our results showed that the mean annual mass balance from 2010 to 2015 was ~−0.9 ± 0.5 m w.e. The mean equilibrium line altitude was ~4680 m in the period 2010–15, which exceeds the maximum elevation of the glacier. The glacier has lost area and mass across its elevation range. The mean ice thickness was 24.0 ± 2.5 m in 2014. From 1972 to 2014, the glacier's area shrank from 0.77 ± 0.05 to 0.39 ± 0.04 km2, and the ice volume decreased by (14.96 ± 0.97) × 106 m3, equivalent to (12.72 ± 0.82) × 106 t w.e. over the same period.

Flow-frictional resistance at the base of glaciers and ice sheets is strongly linked to subglacial water pressure. Understanding the physical mechanisms that govern meltwater fluxes in subglacial channels is hence critical for constraining variations in ice flow. Previous mathematical descriptions of soft-bed subglacial channels assume a viscous till rheology, which is inconsistent with laboratory data and the majority of field studies. Here, we use a grain-scale numerical formulation coupled to pore-water dynamics to analyze the structural stability of channels carved into soft beds. Contrary to the soft-bed channel models assuming viscous till rheology, we show that the flanks of till channels can support substantial ice loads without creep closure of the channel, because the sediment has finite frictional strength. Increased normal stress on the channel flanks causes plastic failure of the sediment, and the channel rapidly shrinks to increase the ice-bed contact area. We derive a new parameterization for subglacial channelized flow on soft beds and show that channel dynamics are dominated by fluvial erosion and deposition processes with thresholds linked to the plastic rheology of subglacial tills. We infer that the described limits to channel size may cause subglacial drainage to arrange in networks of multiple closely spaced channels.

We present the first dedicated study into the phenomenon of ice sails. These are clean ice structures that protrude from the surface of a small number of debris-covered glaciers and can grow to heights of over 25 m. We draw together what is known about them from the academic/exploration literature and then analyse imagery. We show here that ice sails can develop by one of two mechanisms, both of which require clean ice to become surrounded by debris-covered ice, where the debris layer is shallow enough for the ice beneath it to melt faster than the clean ice. Once formed, ice sails can persist for decades, in an apparently steady state, before debris layer thickening eventually causes a reversal in the relative melt rates and the ice sails decay to merge back with the surrounding glacier surface. We support our image-based analysis with a surface energy-balance model and show that it compares well with available observations from Baltoro Glacier in the Karakoram. A sensitivity analysis of the model is performed and confirms the results from our empirical study that ice sails require a relatively high evaporative heat flux and/or a relatively low sensible heat flux in order to exist.

Surface melt in the ablation zone is dominated by atmospheric temperature and surface albedo. We developed a surface mass-balance model with a dynamic component of glacier ice albedo which includes surface properties, clouds and the angle of the sun. The ice albedo reduction is mainly caused by impurity accumulation of non-biological origin such as dust and black carbon (BC), which is currently not included in other surface mass-balance models. Simulations show that dust from meltout is the main source of impurity mass at the melting glacier ice surface, and current rates of atmospheric deposition of dust play only a minor role. However, for BC the atmospheric deposition is the main source where ice melt rates are below 1 m, and atmospheric deposition is most likely from intercontinental transport due to the scarce population and lack of forests in Greenland.

Glacial earthquakes are caused by large iceberg calving events, which are an important mechanism for mass loss from the Greenland ice sheet. The number of glacial earthquakes in Greenland has increased sixfold over the past two decades. We use teleseismic surface waves to analyze the 145 glacial earthquakes that occurred in Greenland from 2011 through 2013, and successfully determine source parameters for 139 events at 13 marine-terminating glaciers. Our analysis increases the number of events in the glacial-earthquake catalog by nearly 50% and extends it to 21 years. The period 2011–13 was the most prolific 3-year period of glacial earthquakes on record, with most of the increase over earlier years occurring at glaciers on Greenland's west coast. We investigate changes in earthquake productivity and geometry at several individual glaciers and link patterns in glacial-earthquake production and cessation to the absence or presence of a floating ice tongue. We attribute changes in earthquake force orientations to changes in calving-front geometry, some of which occur on timescales of days to months. Our results illustrate the utility of glacial earthquakes as a remote-sensing tool to identify the type of calving event, the grounded state of a glacier, and the orientation of an active calving front.

Knowledge of debris-free and debris-covered glaciers is important for understanding the varying response of glaciers to climate change. Measurements at the debris-free Parlung No. 4 Glacier and the debris-covered 24 K Glacier in the southeastern Tibetan Plateau were carried out to compare the meteorology and surface energy fluxes and to understand the factors controlling the melting process. The meteorological comparisons displayed temporally synchronous fluctuations in air temperature, relative humidity and incoming longwave radiation (Lin), but notable differences in precipitation, incoming shortwave radiation (Sin) and wind speed. Under the prevailing regional precipitation and debris conditions, more Lin (42 W m−2) was supplied from warmer and more humid air and more Sin (58 W m−2) was absorbed at the 24 K Glacier. The relatively high energy supply led mainly to an increased energy output via turbulent heat fluxes and outgoing longwave radiation, rather than glacier melting beneath the thick debris. The sensitivity experiment showed that melt rates were sensitive to energy supply at debris thicknesses <~10 cm. In contrast, energy supply to the Parlung No. 4 Glacier mainly resulted in snow/ice melting, the magnitude of which was significantly influenced by energy supplied by Sin and the sensible heat flux.

Quantifying historic changes in glacier size and mass balance is important for understanding how the cryosphere responds to climate variability and change. Airborne photogrammetry enables glacier extent and equilibrium line altitudes (ELAs) to be monitored for more glaciers at lower cost than traditional mass-balance programs and other remote-sensing techniques. Since 1977, end-of-summer-snowlines, which are a proxy for annual ELAs, have been recorded for 50 glaciers in the Southern Alps of New Zealand using oblique aerial photographs. In this study, we use structure from motion photogrammetry to estimate the camera parameters, including position, for historic photographs, which we then use to measure glacier change. We apply this method to a small maritime New Zealand glacier (Brewster Glacier, 1670–2400 m a.s.l.) to derive annual ELA and length records between 1981 and 2017, and quantify the uncertainties associated with the method. Our length reconstruction shows largely continuous terminus retreat of 365 ± 12 m for Brewster Glacier since 1981. The ELA record, which compares well with glaciological mass-balance data measured between 2005 and 2015, shows pronounced interannual variability. Mean ELAs range from 1707 ± 6 to 2303 ± 5 m a.s.l., with the highest ELAs occurring in the last decade.