Recognising the scarcity of glacier mass-balance data in the Southern Hemisphere, a mass-balance measurement programme was started at Brewster Glacier in the Southern Alps of New Zealand in 2004. Evolution of the measurement regime over the 11 years of data recorded means there are differences in the spatial density of data obtained. To ensure the temporal integrity of the dataset a new geostatistical approach is developed to calculate mass balance. Spatial co-variance between elevation and snow depth allows a digital elevation model to be used in a co-kriging approach to develop a snow depth index (SDI). By capturing the observed spatial variability in snow depth, the SDI is a more reliable predictor than elevation and is used to adjust each year of measurements consistently despite variability in sampling spatial density. The SDI also resolves the spatial structure of summer balance better than elevation. Co-kriging is used again to spatially interpolate a derived mean summer balance index using SDI as a co-variate, which yields a spatial predictor for summer balance. The average glacier-wide surface winter, summer and annual balances over the period 2005–15 are 2484, −2586 and −102 mm w.e., respectively, with changes in summer balance explaining most of the variability in annual balance.

We present climate data, direct surface mass balance (SMB) observations and model results for Mocho Glacier in the Chilean Lake District. Mean annual temperature on a nunatak of Mocho Glacier at an elevation of ~2000 m was +2.6°C in 2006–15 and mean annual precipitation in Puerto Fuy (13 km from the glacier, at an elevation of 600 m) was 4000 mm for the same period. High interannual variations in the SMB of Mocho Glacier were observed. A simple SMB model is able to reproduce the observed annual variations in SMB, but fails to predict the steep observed mass-balance gradient. The average of the measured annual glacier mass balances in the four hydrological years 2009/10–2012/13 was −0.90 m w.e. a−1 and the average modelled annual glacier mass balance 2006/07–2014/15 was −1.05 m w.e. a−1. The observed distributed ablation shows a clear altitudinal dependency, while accumulation is determined by patterns of snow drift as well. These patterns are only poorly represented in the model and have to be included in order to be able to reproduce a realistic SMB map of the glacier.

Microtremor measurements and the horizontal-to-vertical spectral ratio (HVSR) technique, generally used for site effect studies as well as to determine the thickness of soft sedimentary layers, can effectively be applied to map the thickness of glaciers. In this work the radio-echo sounding, geoelectric and active seismic methods, widely employed to image the earth interior, are applied to verify the reliability of the HVSR technique in Alpine and Antarctic glacial environments. The technique has been used to analyze passive seismic data from glaciers of the Adamello and Ortles-Cevedale massifs (Italy), the Bernese Oberland Alps (Switzerland) and from the Whillans Ice Stream (West Antarctica). Comparing with the results obtained from the different geophysical imaging methods, we show that the resonance frequency in the HVSR spectra correlates well with the ice thickness at the site, in a wide range from a few tens of meters to more than 800 m. The reliability of the method mainly depends on the coupling of sensors at the glacier surface and on the basal impedance contrast. This passive seismic technique offers a logistically efficient and cost effective method to map glacier and ice-sheet thicknesses. Moreover, under certain conditions, it allows reliable estimations of the basal seismic properties.

Particle size, particle speed and airflow speed have been measured in the powder snow clouds of avalanches to investigate the suspension and transportation processes of snow particles. The avalanches were artificially triggered at the Lautaret full-scale avalanche test-site (French Alps) where an ultrasonic anemometer and a snow particle counter were setup in an avalanche track for measurements. Relatively large particles were observed during passage of the avalanche head and then the size of the particles slightly decreased as the core of the avalanche passed the measurement station. The particle size distribution was well fitted by a gamma distribution function. A condition for suspension of particles within the cloud based on the ratio of vertical velocity fluctuation to particle settling velocity suggests that the large particles near the avalanche head are not lifted up by turbulent diffusion, but rather ejected by a process involving collisions between the avalanche flow and the rough snow surface. Particle speeds were lower than the airflow speed when large particles were present in the powder cloud.

Despite renewed efforts to better understand glacier change and recognize glacier change trends in the Andes, relatively large areas in the Andes of Argentina and Chile are still not investigated. In this study, we report on glacier elevation and mass changes in the outer region of the Northern and Southern Patagonian Icefields in the Southern Patagonian Andes. A newly-compiled Landsat ETM+ derived glacier inventory (consisting of 2253 glaciers and ~1314 ± 66 km2 of ice area) and differencing of the SRTM and SPOT5 DEMs were used to derive glacier-specific elevation changes over the 2000–12 period. The investigated glaciers showed a volume change of −0.71 ± 0.55 km3 a−1, yielding a surface lowering of 0.52 ± 0.35 m a−1 on average and an overall mass loss of 0.46 ± 0.37 m w.e. a−1. Highly variable individual glacier responses were observed and interestingly, they were less negative than previously reported for the neighboring Patagonian Icefields.

Glacier area changes on the Tibetan Plateau were studied in different drainage basins based on Landsat satellite images from three epochs: 263 in the mid-1970s, 150 in 1999–2002 and 148 in 2013/14. Three mosaics (M1976, M2001 and M2013) with minimal cloud and snow cover were constructed, and the uncertainty due to each epoch having a finite span was accounted for. Glacier outlines (TPG1976, TPG2001 and TPG2013) were digitized manually with guidance from the SRTM DEM v4.1 and Google Earth imagery. To achieve complete multi-temporal coverage in a reasonable time, only debris-free ice was delineated. Area mapping uncertainty was evaluated at three study sites, Mount Qomolangma (Everest), Mount Naimona'Nyi, Mount Geladandong, where the largest differences between present and earlier measurements were within ~±4%. Area differences with previous inventories ranged from −19.6% (TPG1976 minus the first Chinese Glacier Inventory) to −3.6% and −1.1% (TPG2013 and TPG2001, respectively minus the second Chinese Glacier Inventory), while the difference TPG2001 minus the GAMDAM Glacier Inventory was +10.4%. Glacier area on the plateau decreased from 44 366 ± 2827 km2 (1.7% of the study area) in the 1970s to 42 210 ± 1621 km2 in 2001 and 41 137 ± 1616 km2 in 2013. Shrinkage was faster in external drainage basins of the southeast than in the interior basins of the northwest, from a maximum of −0.43% a−1 (−1.60% a−1 during 1994–2013) in the Mekong catchment down to a minimum of −0.12% a−1 in the Tarim interior drainage.

We assess ice loss and velocity changes between 1985 and 2014 of three tidewater and five-land terminating glaciers in Godthåbsfjord (Nuup Kangerlua), Greenland. Glacier thinning accounted for 43.8 ± 0.2 km3 of ice loss, equivalent to 0.10 mm eustatic sea-level rise. An additional 3.5 ± 0.3 km3 was lost to the calving retreats of Kangiata Nunaata Sermia (KNS) and Narsap Sermia (NS), two tidewater glaciers that exhibited asynchronous behavior over the study period. KNS has retreated 22 km from its Little Ice Age (LIA) maximum (1761 AD), of which 0.8 km since 1985. KNS has stabilized in shallow water, but seasonally advects a 2 km long floating tongue. In contrast, NS began retreating from its LIA moraine in 2004–06 (0.6 km), re-stabilized, then retreated 3.3 km during 2010–14 into an over-deepened basin. Velocities at KNS ranged 5–6 km a−1, while at NS they increased from 1.5 to 5.5 km a−1 between 2004 and 2014. We present comprehensive analyses of glacier thinning, runoff, surface mass balance, ocean conditions, submarine melting, bed topography, ice mélange and conclude that the 2010–14 NS retreat was triggered by a combination of factors but primarily by an increase in submarine melting.

Understanding the drivers of recent change at Greenlandic tidewater glaciers is of great importance if we are to predict how these glaciers will respond to climatic warming. A poorly constrained component of tidewater glacier processes is the near-terminus subglacial hydrology. Here we present a novel method for constraining near-terminus subglacial hydrology with application to marine-terminating Kangiata Nunata Sermia in South-west Greenland. By simulating proglacial plume dynamics using buoyant plume theory and a general circulation model, we assess the critical subglacial discharge, if delivered through a single compact channel, required to generate a plume that reaches the fjord surface. We then compare catchment runoff to a time series of plume visibility acquired from a time-lapse camera. We identify extended periods throughout the 2009 melt season where catchment runoff significantly exceeds the discharge required for a plume to reach the fjord surface, yet we observe no plume. We attribute these observations to spatial spreading of runoff across the grounding line. Persistent distributed drainage near the terminus would lead to more spatially homogeneous submarine melting and may promote more rapid basal sliding during warmer summers, potentially providing a mechanism independent of ocean forcing for increases in atmospheric temperature to drive tidewater glacier acceleration.

Petermann Gletscher drains ~4% of the Greenland ice sheet (GrIS) area, with ~80% of its mass loss occurring by basal melting of its ice shelf. We use a high-resolution coupled ocean and sea-ice model with a thermodynamic glacial ice shelf to diagnose ocean-controlled seasonality in basal melting of the Petermann ice shelf. Basal melt rates increase by ~20% in summer due to a seasonal shift in ocean circulation within Nares Strait that is associated with the transition from landfast sea ice to mobile sea ice. Under landfast ice, cold near-surface waters are maintained on the eastern side of the strait and within Petermann Fjord, reducing basal melt and insulating the ice shelf. Under mobile sea ice, warm waters are upwelled on the eastern side of the strait and, mediated by local instabilities and eddies, enter Petermann Fjord, enhancing basal melt down to depths of 200 m. The transition between these states occurs rapidly, and seasonal changes within Nares Strait are conveyed into the fjord within the same season. These results suggest that long-term changes in the length of the landfast sea-ice season will substantially alter the structure of Petermann ice shelf and its contribution to GrIS mass loss.

An anomalously slight glacier mass gain during 2000 to the 2010s has recently been reported in the Karakoram region. However, to date, no investigations of the region-wide glacier mass balance in the Karakoram prior to 2000 have been reported, leaving a knowledge gap for assessing glacier responses to climate change. We calculated elevation and mass change using DEMs generated from KH-9 images acquired during 1973–1980 and the 1 arc-second SRTM DEM. We find a slight mass loss of −0.09 ± 0.03 m w.e. a−1 (12 366 km2) for 1973–2000, which is less negative than the global average rate for 1971–2009 (−0.31 ± 0.19 m w.e. a−1). Within the Karakoram, the glacier change patterns are spatially and temporally heterogeneous. In particular, a nearly stable state in the central Karakoram (−0.04 ± 0.05 m w.e. a−1 during the period 1974–2000) implies that the Karakoram anomaly dates back to the 1970s. Combined with the previous studies, we conclude that the Karakoram glaciers as a whole were in a nearly balanced state during the 1970s to the 2010s.

We analyzed glacier surface elevations (1957, 2010 and 2015) and surface mass-balance measurements (2008–2015) on the 30 km2 Eklutna Glacier, in the Chugach Mountains of southcentral Alaska. The geodetic mass balances from 1957 to 2010 and 2010 to 2015 are −0.52 ± 0.46 and −0.74 ± 0.10 m w.e. a−1, respectively. The glaciological mass balance of −0.73 m w.e. a−1 from 2010 to 2015 is indistinguishable from the geodetic value. Even after accounting for loss of firn in the accumulation zone, we found most of the mass loss over both time periods was from a broad, low-slope basin that includes much of the accumulation zone of the main branch. Ice-equivalent surface elevation changes in the basin were −1.0 ± 0.8 m a−1 from 1957 to 2010, and −0.6 ± 0.1 m a−1 from 2010 to 2015, shifting the glacier hypsometry downward and resulting in more negative mass balances: an altitude-mass-balance feedback. Net mass loss from Eklutna Glacier accounts for 7 ± 1% of the average inflow to Eklutna Reservoir, which is entirely used for water and power by Anchorage, Alaska's largest city. If the altitude-mass-balance feedback continues, this ‘deglaciation discharge dividend’ is likely to increase over the short-term before it eventually decreases due to diminishing glacier area.

The instrumented sample holder Snowbreeder 5 is used to investigate the simultaneous influence of settlement on temperature-gradient snow metamorphism in time-lapse micro-computed tomography experiments. So far, experiments have only been done on temperature-gradient snow metamorphism without settlement or settlement under isothermal conditions. With the new device we can impose a constant temperature gradient on a snow sample and induce settlement by placing a passive load on top of the snow sample. The weight of the load can be varied, simulating various snow heights on top of the snow sample. Snow-temperature measurements on the passive load are possible due to wireless data transfer via Bluetooth. The temperature gradient is set by controlling the air temperature inside the computer tomograph and by a Peltier element at the bottom of the snow sample. First experiments under isothermal conditions and a constant temperature gradient of 43 K m−1 showed that the settlement was reduced to almost half as soon as a temperature gradient was applied under otherwise almost equal snow conditions. The compactive viscosity in the isothermal experiment was in the range of literature values.

Improvements in understanding glacial extents and chronologies for the southeastern slope of the western Nyaiqentanglha Shan on the Tibetan Plateau are required to understand regional climate changes during the Last Glacial cycle. A two-dimensional numerical model of mass balance, based on snow–ice melting factors, and of ice flow for mountain glaciers is used to assess the glacier sensitivity to climatic change in a catchment of the region. The model can reproduce valley glaciers, wide-tongued glaciers and a coalescing glacier within step temperature lowering and precipitation increasing experiments. The model sensitivity experiments also indicate that the dependence of glacier growth on temperature and/or precipitation is nonlinear. The model results suggest that the valley glaciers respond more sensitively to an imposed climate change than wide-tongued and coalescing glaciers. Guided by field geological evidence of former glacier extent and other independent paleoclimate reconstructions, the model is also used to constrain the most realistic multi-year mean temperatures to be 2.9–4.6°C and 1.8–2.5°C lower than present in the glacial stages of the Last Glacial Maximum and middle marine oxygen isotope stage 3, respectively.

Versatile instruments assembled from off-the-shelf sensors and open-source electronics are used to record wave propagation and damping measured by Inertial Motion Units (IMUs) in a grease ice slick near the shore in Adventfjorden, Svalbard. Viscous attenuation of waves due to the grease ice slick is clearly visible by comparing the IMU data recorded by the different instruments. The frequency dependent spatial damping of the waves is computed by comparing the power spectral density obtained from the different IMUs. We model wave attenuation using the one-layer model of Weber from 1987. The best-fit value for the effective viscosity is ν = (0.95 ± 0.05 × 10−2)m2 s−1, and the coefficient of determination is R2 = 0.89. The mean absolute error and RMSE of the damping coefficient are 0.037 and 0.044m−1, respectively. These results provide continued support for improving instrument design for recording wave propagation in ice-covered regions, which is necessary to this area of research as many authors have underlined the need for more field data.