We collected ~1300 km of ground-penetrating radar profiles over McMurdo Ice Shelf, Antarctica, using frequencies between 40 and 400 MHz to determine extent, continuity and depth to the brine. We also used profiles to determine meteoric ice thickness and locate englacial features, which may suggest ice shelf instability. The brine extends 9–13 km inland from the ice shelf terminus and covers the entire region between Ross, White and Black Islands. Jump unconformities and basal fractures exist in the brine and ice shelf, respectively, suggesting prior fracturing and re-suturing. One 100 MHz profile, the most distal from the ice shelf edge while still being situated over the brine, simultaneously imaged the brine and bottom of meteoric ice. This suggests a negative brine salinity gradient moving away from the terminus. The meteoric ice bottom was also imaged in a few select locations through blue ice in the ablation zone near Black Island. We suggest that brine, sediment-rich ice and poor antenna coupling on rough ice attenuates the signal in this area. When combined with other recent mass-balance and structural glaciology studies of MIS, our results could contribute to one of the most high-resolution physical models of an ice shelf in Antarctica.

An improved understanding of fresh water stored in the Himalaya is crucial for water resource management in South Asia and can be inferred from glacier mass-balance estimates. However, field investigations in the rugged Himalaya are limited to a few individual glaciers and short duration. Therefore, we have recently developed an approach that combines satellite-derived snowlines, a temperature-index melt model and the accumulation-area ratio method to estimate annual mass balance of glaciers at basin scale and for a long period. In this investigation, the mass balance of 146 glaciers in the Chandra basin, western Himalaya, is estimated from 1984 to 2012. We estimate the trend in equilibrium line altitude of the basin as +113 m decade−1 and the mean mass balance as −0.61 ± 0.46 m w.e. a−1. Our basin-wide mass-balance estimates are in agreement with the geodetic method during 1999–2012. Sensitivity analysis suggests that a 20% increase in precipitation can offset changes in mass balance for a 1 °C temperature rise. A water loss of 18% of the total basin volume is estimated, and 67% for small and low-altitude glaciers during 1984–2012, indicating a looming water scarcity crisis for villages in this valley.

This sensitivity study applies the offline Canadian Land Surface Scheme (CLASS) version 3.6 to simulate snowpack evolution in idealized topography using observations at Likely, British Columbia, Canada over 1 July 2008 to 30 June 2009. A strategy for a subgrid-scale snow (SSS) parameterization is developed to incorporate two key features: ten elevation bands at 100 m intervals to capture air temperature lapse rates, and five slope angles on four aspects to resolve solar radiation impacts on the evolution of snow depth and SWE. Simulations reveal strong elevational dependencies of snow depth and SWE when adjusting temperatures using a moist adiabatic lapse rate with elevation, with 26% peak SWE differences between that at the average elevation versus the mean of the remainder of the elevation bands. Differences in peak SWE on north- and south-facing slopes increase from 3.0 mm at 10° slope to 17.9 mm at 50° slope. When applied to elevation, slope and aspect combinations derived from a high-resolution digital elevation model, elevation dominates the control of peak SWE values. Inclusion of the range of SSS effects into a regional climate model will improve snowpack and hydrological simulations of western North America's snow-dominated, mountainous watersheds.

Avalanching is a prominent source of accumulation on glaciers that have high and steep valley-walls surrounding their accumulation zones. These glaciers are typically characterised by an extensive supraglacial debris cover and a low accumulation area ratio. Despite an abundance of such glaciers in the rugged landscapes of the High Himalaya, attempts to quantify the net avalanche contribution to mass balance and its long-term variation are almost missing. We first discuss diagnostic criteria to identify strongly avalanche-fed glaciers. Second, we develop an approximate method to quantify the magnitude of the avalanche accumulation exploiting its expected control on the dynamics of these glaciers. The procedure is based on a simplified flowline model description of the glacier concerned and utilises the known glaciological mass-balance, velocity and surface-elevation profiles of the glacier. We apply the method to three Himalayan glaciers and show that the data on the recent dynamics of these glaciers are consistent with a dominant contribution of avalanches to the total accumulation. As a control experiment, we also simulate another Himalayan glacier where no significant avalanche contribution is expected, and reproduce the recent changes in that glacier without any additional avalanche contribution.

In this study, seasonality and interannual variability of snow depth at two stations (Hornsund and Barentsburg) located in western Spitsbergen are investigated. For this purpose, the novel Moving Average over Shifting Horizon method combined with trend estimation was used. The Hornsund and Barentsburg stations exhibit similar snow depth trends during early autumn and late spring when statistically significant decreases were estimated at both stations (for August 1984–July 2016). In the remaining period, there are differences in outcomes between stations. The results for Barentsburg from October till the end of May are characterised by the lack of a trend while at Hornsund decreases of snow depth were estimated. The largest changes occur in May when the snow depth was at its maximum. Differences in the estimated tendencies were explained with the help of a trend analysis for air temperature and precipitation. An analysis of maximum snow depth, snow onset date, snow disappearance date and snow-cover duration is included. The results of the assessment depend on the location, with a lack of statistically significant changes in Barentsburg, and later snow onset date, shorter duration and decrease of maximum depth in Hornsund.

Ocean tides influence the flow of marine-terminating glaciers. Observations indicate that the large fortnightly variations in ice flow at Rutford Ice Stream in West Antarctica originate in the floating ice shelf. We show that nonlinear variations in ice shelf buttressing driven by tides can produce such fortnightly variations in ice flow. These nonlinearities in the tidal modulation of buttressing stresses can be caused by asymmetries in the contact stress from migration of the grounding line and bathymetric pinning points beneath the ice shelf. Using a simple viscoelastic model, we demonstrate that a combination of buttressing and hydrostatic stress variations can explain a diverse range of tidal variations in ice shelf flow, including the period, phase and amplitude of flow variations observed at Rutford and Bindschadler Ice Streams.

Albedo is an important parameter in the energy balance of bare-ice surfaces and modulates glacier melt rates. The prolongation of the ablation period enforces the albedo feedback and highlights the need for profound knowledge on impacts of bare-ice albedo on glacier mass balance. In this study, we assess the mass balance sensitivity of 12 Swiss glaciers with abundant long-term in-situ data on changes in bare-ice albedo. We use pixel-based bare-ice albedo derived from Landsat 8. A distributed mass-balance model is applied to the period 1997–2016 and experiments are performed to assess the impact of albedo changes on glacier mass balance. Our results indicate that glacier-wide mass-balance sensitivities to changes in bare-ice albedo correlate strongly with mean annual mass balances (r2 = 0.81). Large alpine glaciers react more sensitively to bare-ice albedo changes due to their ablation areas being situated at lower elevations. We find average sensitivities of glacier-wide mass balance of −0.14 m w.e. a−1 per 0.1 albedo decrease. Although this value is considerably smaller than sensitivity to air temperature change, we stress the importance of the enhanced albedo feedback that will be amplified due to atmospheric warming and a suspected darkening of glacier surface in the near future.

Seasonal snow-cover modulates water and energy budgets across large areas of the Northern Hemisphere. Previous research, based on satellite imagery interpreted and curated by the Rutgers University Snow Laboratory, has identified significant negative and positive trends in annual snow-covered duration and area at hemispheric and continental scales between 1971 and 2014. This study uses the same dataset to generate more detailed descriptions of spatial variations in these trends, maps intraannual variations in sign, statistical significance and strength, and quantifies associations with latitude and elevation. It also considers the limitations and uncertainties associated with a binary classification of this type, and the implications for trend magnitudes of adopting alternatives to the conventional assumption of 100% (0%) actual fractional snow-covered area in ‘snow-covered’ (‘snow-free’) spatial units at different stages of the snow-season. This prompts adoption of alternative terminology, referring to ‘snow-dominated’ area and duration. In response to questions about the dataset's veracity raised by some prior studies, it discusses climatological factors of potential relevance in explaining spatio-temporal trend patterns, and considers how biases might possibly have been introduced as a result of extraneous influences.

Ross Ice Shelf (RIS) is known to experience transient thickness change due to changes in the flow of its tributary ice streams and glaciers and this may complicate identification of external, climate-forced signals in contemporary observations of ice shelf thinning and thickening. Flux changes at the lateral boundaries produce both instantaneous and longer timescale adjustments in the coupled velocity and thickness fields. Here, we adapt a statistical approach to output from a numerical model of ice shelf flow to identify characteristic patterns (spatial response surfaces) associated with stepped and cyclic perturbations to boundary fluxes. Once known, characteristic patterns identified in observational data may be attributed to specific sources. An example involving discharge of Byrd Glacier into RIS is described. We find that spatial response surfaces for thickness and velocity generated in individual flow model experiments appear to be independent of flux perturbation shape and magnitude. Additionally, recent acceleration of Byrd Glacier is apparent in ICESat-detected change in RIS thickness.

A third-order linear glacier length model is used to analyze if the retreat of Hintereisferner in the Austrian Alps over the past 160 years is exceptional, or lies within the range of the natural variability inherent to a stationary climate. A detailed uncertainty analysis takes into account glacier geometry, model parameters and initial conditions. A Monte Carlo ensemble strengthens the result that the observed retreat cannot be explained by natural variability and therefore affirms regional climate change. Finally the observed temperature trend at Hintereisferner lies outside the range of natural variability from an ensemble of climate models, but is consistent with the modeled range of responses to anthropogenic forcing.

Information on the occurrence, spatial distribution and morphometric characteristics of Antarctic icebergs is needed in a large number of applications including navigation, heat and freshwater balance calculations, biochemistry of the ocean and climatology. Using over 60 000 ship observations of icebergs in the Southern Ocean collected since the end of the 1940s we have produced a detailed map of the distribution of Antarctic icebergs as well as maps of related statistics including the standard deviation, minimum and maximum values of the iceberg concentration and the probability of iceberg-free observations. The study incorporated small and medium-sized icebergs with a length of <10 nautical miles. Most observations were taken during the warm period of the year, from December to April. It is shown that the iceberg distribution across the Southern Ocean is determined by the location of calving regions and peculiarities of the atmospheric circulation and ocean currents. Iceberg concentration data combined with information on the iceberg size and shape distribution have been used to evaluate the area-integrated characteristics of Antarctic icebergs. The instantaneous number of icebergs in the Southern Ocean was estimated as 132 269 with an uncertainty of 7%. The area and volume of icebergs were equal correspondingly to 55 805 km2 and 16 893 km3 with uncertainties of 32–33%.

The goal of this study is to test applicability of the conceptual catchment run-off HBV model to simulate discharge in small non-glaciated Arctic catchment. Within two ablation seasons, 2014 and 2015, in the Fuglebekken catchment (Spitsbergen, Svalbard), selected hydro-meteorological measurements were conducted, including discharge measurements in 10 min interval by Nivus PCM-F device with active Doppler sensor. The model parameters were calibrated on discharge measurements from both years separately and verified independently. As the transformation from rainfall to runoff includes a number of processes with different dynamics and timescales, the proper description of the processes and their simulation of discharge depend on the temporal resolution of the data. For that purpose, the relationships between the calibration and validation results, and optimal model parameters with different time steps were analyzed. It was found that calculated fit of simulated to observed discharge, depends on the year, time step and data averaging. The best results were obtained for the model from year 2015 for 3 and 6 h using averaged input data.

The history of glacier length fluctuations serves as a reliable indicator of the past climate. In this paper, a numerical flowline model has been used to study the relationship between length variations of Chhota Shigri glacier and local climate since 1876. The simulated front positions of Chhota Shigri glacier are in agreement with those observed. After a successful simulation of the past retreat, the model was also used to predict future evolution of the glacier for the next 100 years under different climatic scenarios. These simulations indicate that the Chhota Shigri glacier may lose ~90% of its present volume by 2100 if the local temperature increases by 2.4 K, and for a temperature rise of 5.5 K, the glacier loses almost all its volume.

Measurements made by an underwater glider deployed near the Ross Ice Shelf were used to identify the presence of Ice Shelf Water (ISW), which is defined as seawater with its potential temperature lower than its surface freezing point temperature. Properties logged by the glider included in situ temperature, electrical conductivity, pressure, GPS location at surfacings and time. For most of the first 30 recorded dives of its deployment, evidence suggests the glider was prevented from surfacing due to being under the ice shelf. For dives under the ice shelf, farthest from the ice shelf front, ISW layers of varying thicknesses and depth locations were observed; between 2 m thick (centred at 231 m depth) to >93 m thick (centred at >360 m). For dives under the ice shelf, close to the ice shelf front, either no ISW was observed or ISW layers were centred at shallower depths (116–127 m). Thicker ISW layers (e.g. up to 250 m thickness centred at 421 m) were observed for some glider dives in open water in front of the Ross Ice Shelf. No in situ supercooling (water colder than the pressure-dependent freezing point temperature) was observed.

Meltwater runoff in the catchment area containing Chhota Shigri glacier (Western Himalaya) is simulated for the period 1951–2099. The applied mass-balance model is forced by downscaled products from four regional climate models with different horizontal resolution. For the future climate scenarios we use high resolution time series of 5 km grid spacing, generated using the newly developed Intermediate Complexity Atmospheric Research Model. The meteorological input is downscaled to 300 m horizontal resolution. The use of an ice flow model provides annually updated glacier area for the mass-balance calculations. The mass-balance model calculates daily snow accumulation, melt, runoff, as well as the individual runoff components (glacial melt, snowmelt and rain). The resulting glacier area decreases by 35% (representative concentration pathway (RCP) 4.5 scenario) to 70% (RCP 8.5 scenario) by 2099 relative to 2000. The average annual mass balance over the whole model period (1951–2099) was –0.4 (±0.3) m w.e. a–1. Average annual runoff does not differ substantially between the two climate scenarios. However, for the years after 2040 our results show a shift towards earlier snowmelt onset that increases runoff in May and June, and reduced glacier melt that decreases runoff in August and September. This shift is much stronger pronounced in the RCP 8.5 scenario.

The seasonal cycle of fast ice thickness in Prydz Bay, East Antarctica, was observed between March and December 2012. In March, we observed a 0.16 m thickness gain of 0.22 m-thick first-year ice (FYI), while 1.16 m-thick second-year ice (SYI) nearby simultaneously ablated by 0.59 m. A 1-D thermodynamic sea-ice model was applied to identify the factors that led to the simultaneous growth of FYI and melt of SYI. The different evolutions were explained by the difference in the conductive heat flux between the FYI and SYI. As the FYI was thin, there was a large temperature gradient between the ice base and the colder ice surface. This generated an upward conductive heat flux, which was larger than the heat flux from the ocean to the ice base, yielding basal growth of ice. In the case of the thicker SYI the temperature gradient and, hence, the conductive heat flux were smaller, and not sufficient to balance the oceanic heat flux at the ice base, yielding basal ablation. Penetration of solar radiation affected the conductive heat flux in both cases, and the model results were sensitive to the initial ice temperature profile and the uncertainty of the oceanic heat flux.

The Drygalski Ice Tongue presents an ~80 km long floating obstacle to alongshore flows in the Victoria Land coastal ocean region of the Western Ross Sea. Here we use oceanographic data from near to the tongue to explore the interplay between the floating glacier and the local currents and stratification. A vessel-based circuit of the glacier, recording ocean temperature and salinity profiles, reveals the southwest corner to be the coldest and most complex in terms of vertical structure. The southwest corner structure beneath the surface warm, salty layer sustains a block of very cold water extending to 200 m depth. In this same location there was a distinct layer at 370 m not seen anywhere else of water at ~−1.93°C. The new observations broadly, but not directly, support the presence of a coherent Victoria Land Coastal Current. The data suggest the northward moving coastal current turns against the Coriolis force and works its way anticlockwise around the glacier, but with leakage beneath the glacier through the highly ‘rippled’ underside, resulting in a spatially heterogeneous supply to the Terra Nova Bay Polynya region – an important location for the formation of high-salinity shelf water.

As the focus of intensive glaciological studies in the 1960–70s, White Glacier on Axel Heiberg Island, Canada, has played an important role in understanding the dynamics of mostly-cold polythermal glaciers in the high Arctic. In this study, we examine the magnitude, duration and timing of peak velocity events in the summers of 2013–15 using continuous dual-frequency GPS observations, and compare them with similar measurements made in 1968. Summer speed-up events in 1968 and 2014, in which ice velocities reached 200% above winter values, were found to occur in conjunction with formation and drainage of an ice-marginal lake. Despite thinning of the glacier by >20 m and a decrease in annual surface velocities of 15–35% since the 1960s, the relative magnitude and duration of these peak events has increased, particularly at lower elevations, in comparison with the observations at the same locations many decades ago. Given the long-term slowdown of the glacier, the relative contribution of summer displacement to the net annual motion has therefore increased significantly, with summer motion over the span of <2 months now accounting for nearly half of the total annual displacement.

Relationships between stable isotopes (δD–δ18O), ice facies and glacier structures have hitherto gone untested in the mid-latitude maritime glaciers of the Southern Hemisphere. Here, we present δD–δ18O values as part of a broader study of the structural glaciology of Fox Glacier, New Zealand. We analyzed 94 samples of δD–δ18O from a range of ice facies to investigate whether isotopes have potential for structural glaciological studies of a rapidly deforming glacier. The δD–δ18O measurements were aided by structural mapping and imagery from terminus time-lapse cameras. The current retreat phase was preceded by an advance of 1 km between 1984 and 2009, with the isotopic sampling and analysis undertaken at the end of that advance (2010/11). Stable isotopes from debris-bearing shear planes near the terminus, interpreted as thrust faults, are isotopically enriched compared with the surrounding ice. When plotted on co-isotopic diagrams (δD–δ18O), ice sampled from the shear planes appears to show a subtle, but distinctive isotopic signal compared with the surrounding clean ice on the lower glacier. Hence, stable isotopes (δD–δ18O) have potential within the structural glaciology field, but larger sample numbers than reported here may be required to establish isotopic contrasts between a broad range of ice facies and glacier structures.

Ice algae are a key component in polar marine food webs and have an active role in large-scale biogeochemical cycles. They remain extremely under-sampled due to the coarse nature of traditional point sampling methods compounded by the general logistical limitations of surveying in polar regions. This study provides a first assessment of hyperspectral imaging as an under-ice remote-sensing method to capture sea-ice algae biomass spatial variability at the ice/water interface. Ice-algal cultures were inoculated in a unique inverted sea-ice simulation tank at increasing concentrations over designated cylinder enclosures and sparsely across the ice/water interface. Hyperspectral images of the sea ice were acquired with a pushbroom sensor attaining 0.9 mm square pixel spatial resolution for three different spectral resolutions (1.7, 3.4, 6.7 nm). Image analysis revealed biomass distribution matching the inoculated chlorophyll a concentrations within each cylinder. While spectral resolutions >6 nm hindered biomass differentiation, 1.7 and 3.4 nm were able to resolve spatial variation in ice algal biomass implying a coherent sensor selection. The inverted ice tank provided a suitable sea-ice analogue platform for testing key parameters of the methodology. The results highlight the potential of hyperspectral imaging to capture sea-ice algal biomass variability at unprecedented scales in a non-invasive way.