Modelling paleo-glacier networks in mountain ranges on the millennial timescales requires ice flow approximations. Hybrid models calculating ice flow by combining vertical shearing (shallow ice approximation) and longitudinal stretching (shallow shelf approximation) have been applied to model paleo-glacier networks on steep terrain, yet their validity has not yet been assessed quantitatively. Moreover, hybrid models consistently yield higher ice thicknesses than Last Glacial Maximum geomorphological reconstructions in the European Alps. Here, we compare results based on the hybrid Parallel Ice Sheet Model (PISM) and the Stokes model Elmer/Ice on the Rhine Glacier, a catchment of the former European Alpine Icefield. For PISM, we also test two magnitudes of flux limitation in a scheme that reduces shearing velocities. We find that the flux limitation typically used in PISM yields significantly reduced shearing speeds and increases ice thicknesses by up to 500 m, partly explaining previous overestimations. However, reducing the ice flux limitation allows the hybrid model to minimize this mismatch and captures sliding speeds, ice thicknesses, ice extent and basal temperatures in close agreement with those obtained with the Stokes model.

The evolution of Northern Hemisphere ice sheets through the last glacial cycle is simulated with the glacial index method by using the climate forcing from one General Circulation Model, COSMOS. By comparing the simulated results to geological reconstructions, we first show that the modelled climate is capable of capturing the main features of the ice-sheet evolution. However, large deviations exist, likely due to the absence of nonlinear interactions between ice sheet and other climate components. The model uncertainties of the climate forcing are examined using the output from nine climate models from the Paleoclimate Modelling Intercomparison Project Phase III. The results show a large variability in simulated ice sheets between the different models. We find that the ice-sheet extent pattern resembles summer surface air temperature pattern at the Last Glacial Maximum, confirming the dominant role of surface ablation process for high-latitude Northern Hemisphere ice sheets. This study shows the importance of the upper boundary condition for ice-sheet modelling, and implies that careful constraints on climate output is essential for simulating realistic glacial Northern Hemisphere ice sheets.

The dynamics of a marine ice sheet's grounding lines determine the rate of ice discharge from the grounded part of ice sheet into surrounding oceans. In many locations in West Antarctica ice flows into ice shelves through ice streams experiencing low driving stress. However, existing simple theories of marine ice sheets are developed under the assumption of high basal and driving stress. Here we analyze the grounding line behavior of marine ice streams experiencing low basal shear and driving stress. We find that in this regime, the ice flux at the grounding line is a complex function of the geometry of the ice-stream bed, net accumulation rate and gradient of the net accumulation rate. Our analysis shows that the stability of distinct steady states is determined by the same parameters, suggesting a more complex (in)stability criterion than what is commonly referred to within the context of the ‘marine ice-sheet instability hypothesis’. We also determine characteristic timescales (e-folding time) of ice-sheet configurations perturbed from their steady states. These timescales can be used to determine whether particular configurations can be considered in isolation from other components of the climate system or whether their effects and feedbacks between the ice sheet and the rest of the climate system have to be taken into account.

Ice-core drilling to depths of 200–300 m is an important part of research studies concerned with paleoclimate reconstruction and anthropogenic climate change. However, conventional drilling methods face difficulties due to firn permeability. We have developed an electromechanical ice-core drill with air reverse circulation at the hole bottom. We believe that the new drilling system will recover ice cores faster than shallow auger drills, with high efficiency and low energy consumption. The theoretically estimated up-hole speed of the airflow should be not <7.7 m s−1 to allow proper removal of ice cuttings from the borehole bottom. The computer simulation and test results showed that the design of the new ice-coring drill is feasible. The maximum allowed penetration rate depends by square law on airflow.

The British Antarctic Survey Rapid Access Isotope Drill is an innovative new class of electromechanical ice drill, which has recently been used to drill the deepest dry hole drilled by an electromechanical auger drill. The record-breaking depth of 461.58 m was drilled in just over 104 hours at Little Dome C. The drill collects ice chippings, for water stable isotope analysis, rather than an ice core. By not collecting a core the winch can be geared for speed rather than core breaking and is lightweight. Furthermore, emptying of the chippings is performed by simply reversing the drill motor on the surface reducing the overall drilling time significantly. The borehole is then available for instrumentation. We describe the drill in its current state including modifications carried out since it was last deployed. Test seasons and the lessons learned from each are outlined. Finally, future developments for this class of drill are discussed.

Karakoram glaciers experienced balanced or slightly positive mass budgets since at least the 1970s. Here, we provide an update on the state of balance of Central and Eastern Karakoram glaciers (12 000 km2) between 2008 and 2016 by differencing DEMs derived from satellite optical images. The mass budget of Central Karakoram glaciers was slightly positive (0.12 ± 0.14 m w.e. a−1) while eastern Karakoram glaciers lost mass (−0.24 ± 0.12 m w.e. a−1). The glacier-wide mass balances of surge-type and nonsurge-type glaciers were not statistically different. Our elevation change data also depict the effect of a > 100 Mm3 rock avalanche on Siachen Glacier ablation area in September 2010. It covered a 4 km2 area with a thick debris layer that unexpectedly, led to locally enhanced glacier mass loss during the following years. Enhanced melt opened a > 100 m deep 2 km2 depression and contributed to 6% of the mass loss of Siachen Glacier from 2010 to 2016 (−0.39 m w.e. a−1). We hypothesize that sub- or englacial melt may be responsible for this intriguing behaviour. This study contributes to a better knowledge of the regional pattern of the Karakoram anomaly and of the influence of rock avalanches on glacier mass changes.

Increasing ocean and air temperatures have contributed to the removal of floating ice shelves from several Greenland outlet glaciers; however, the specific contribution of these external forcings remains poorly understood. Here we use atmospheric, oceanographic and glaciological time series data from the ice shelf of Petermann Gletscher, NW Greenland to quantify the forcing of the ocean and atmosphere on the ice shelf at a site ~16 km from the grounding line within a large sub-ice-shelf channel. Basal melt rates here indicate a strong seasonality, rising from a winter mean of 2 m a−1 to a maximum of 80 m a−1 during the summer melt season. This increase in basal melt rates confirms the direct link between summer atmospheric warming around Greenland and enhanced ocean-forced melting of its remaining ice shelves. We attribute this enhanced melting to increased discharge of subglacial runoff into the ocean at the grounding line, which strengthens under-ice currents and drives a greater ocean heat flux toward the ice base.

From field observations, dry snow slab avalanche initiation is associated with fracture within relatively thin weak layers under stronger, cohesive slabs. For risk-based avalanche prediction, it is important to understand the fracture properties of alpine snow. Alpine snow is a quasi-brittle material with a fracture mechanical size effect on nominal shear strength meaning that strength decreases with increasing specimen size. A related size effect is the critical length required for rapid propagation of a shear fracture. In that case, the probability of fracture increases with increasing crack length. In this paper, 45 sets of field-measured critical lengths are presented based on 591 individual tests. From analysis, a probabilistic size effect law based on critical lengths is derived analogous to the deterministic size effect law for nominal shear strength related to fracture mechanics. The new size effect law may be useful in applications, particularly since the critical length is easily measured in the field and it is a principal component of weak layer fracture toughness.

A linear model was used to investigate the behaviour of the Italian Western Tauri glaciers as a result of temperature changes projected for the 21st century. The model estimates the temperature variations once the glacier snout length variations are known and vice versa: it estimates the glacier snout length variations once the air temperature variations are known or predicted. The 46 glaciers of the Italian Western Tauri, 35 (76%) of which have areas smaller than 0.5 km2 and only 7 (15%) larger than 1 km2, are mostly mountain type glaciers and only three are presently valley glaciers. The model has been forced by the air temperature projections of the A1B emission scenario, which indicates an increase in temperature of 2.7°C from 2015 to 2100. The results show a shortening of more than 35% for mountain glaciers by 2100 with a surface loss of more than 60% and smaller reductions for valley glaciers. The consequent fragmentation into smaller units would lead to the extinction of 95% of the existing glaciers by the end of the century, possibly leaving only the valley glaciers surviving.

A new type of ice core drill bit, designed with a vane swirler, was developed for ice core drilling with air reverse circulation. An orthogonal experimental design method was employed to investigate the effects of the swirler structure parameters on the reverse circulation performance of the drill bit including helical angle, number of blades, blade length and blade central angle, etc. The entrainment ratio was used to evaluate the reverse circulation effectiveness of the drill bit. The results show that the helical angle is the dominant factor regardless of whether or not the flushing nozzles are part of the design of the drill bit. The number of blades is the least important factor for the drill bit designed with the flushing nozzles (referred to as drill bit I), while the outlet area of the swirling slot is the least influential factor for the drill bit without flushing nozzles (referred to as drill bit П). In addition, the appearance of the ice core has a certain effect on the air reverse circulation for both drill bits. Within the ranges of this study, the optimal structure of the drill bit was determined based on the range analysis of the orthogonal design.

Changes in glacier length and extent are indicators of contemporary and archives of past climate changes, but this common climate proxy presents a challenge for inferring a climate signal. Modeling studies suggest that length fluctuations can occur due to interannual climate variability within an unchanging mean climate and that changes in interannual climate variability can also drive changes in average length. This paper quantifies the impacts of interannual climate variability on average glacier length and mass balance, using a flowline model coupled to a simplified mass-balance model. Results illustrate that changes in the magnitude of interannual temperature variability can non-linearly affect the mean glacier length through a mass-balance asymmetry between warm and cold years. This asymmetry is present in models where melt only initiates after a temperature threshold is crossed. Glaciers susceptible to this asymmetry can be identified based on the shape of their mass-balance profiles. The presence of mass-balance asymmetries in glaciological databases is evaluated, but current records are too short for high statistical resolving power. While the asymmetry in this study can affect the average length and mass-balance, its impacts are small, and paleoclimate interpretations from glacier-length changes are likely not notably influenced by this process.

Mountain glaciers integrate climate processes to provide an unmatched signal of regional climate forcing. However, extracting the climate signal via intercomparison of regional glacier mass-balance records can be problematic when methods for extrapolating and calibrating direct glaciological measurements are mixed or inconsistent. To address this problem, we reanalyzed and compared long-term mass-balance records from the US Geological Survey Benchmark Glaciers. These five glaciers span maritime and continental climate regimes of the western United States and Alaska. Each glacier exhibits cumulative mass loss since the mid-20th century, with average rates ranging from −0.58 to −0.30 m w.e. a−1. We produced a set of solutions using different extrapolation and calibration methods to inform uncertainty estimates, which range from 0.22 to 0.44 m w.e. a−1. Mass losses are primarily driven by increasing summer warming. Continentality exerts a stronger control on mass loss than latitude. Similar to elevation, topographic shading, snow redistribution and glacier surface features often exert important mass-balance controls. The reanalysis underscores the value of geodetic calibration to resolve mass-balance magnitude, as well as the irreplaceable value of direct measurements in contributing to the process-based understanding of glacier mass balance.

Multiple studies on regional glacier mass balance in the Pamir Mountains have been conducted using the geodetic method, but they were rarely extended to the period before 2000. In this study, we used KH-9 imagery acquired in 1975 to generate the historical DEM for the central Pamir, and then obtained the glacier elevation change by comparing this with the SRTM C-band DEM. The penetration depth of the C-band radar was corrected for different glacier surfaces, i.e. 2.96, 1.68 and 0 m for firn/snow cover, bare ice and debris-covered areas, respectively. The final results suggest that the central Pamir glaciers, overall, experienced a near-zero mass balance of −0.03 ± 0.24 m w.e. a−1 for 1975–99. Due to glacier surge activity, the elevation change patterns of individual glaciers were highly variable, and their mass balances varied from −0.12 ± 0.26 to 0.63 ± 0.20 m w.e. a−1. The mean mass budgets of the surge-type glaciers and non-surge-type glaciers were 0.03 ± 0.14 and −0.05 ± 0.28 m w.e. a−1, respectively. Concurring with previous studies, we conclude that the central Pamir glaciers may have been in a state of approximately balanced mass budget or slight mass deficit from the mid-1970s to the mid-2010s.

The spatial distribution of basal water critically impacts the evolution of ice sheets. Current estimates of basal water distribution beneath the Greenland Ice Sheet (GrIS) contain large uncertainties due to poorly constrained boundary conditions, primarily from geothermal heat flux (GHF). The existing GHF models often contradict each other and implementing them in numerical ice-sheet models cannot reproduce the measured temperatures at ice core locations. Here we utilize two datasets of radar-detected basal water in Greenland to constrain the GHF at regions with a thawed bed. Using the three-dimensional ice-sheet model SICOPOLIS, we iteratively adjust the GHF to find the minimum GHF required to reach the bed to the pressure melting point, GHFpmp, at locations of radar-detected basal water. We identify parts of the central-east, south and northwest Greenland with significantly high GHFpmp. Conversely, we find that the majority of low-elevation regions of west Greenland and parts of northeast have very low GHFpmp. We compare the estimated constraints with the available GHF models for Greenland and show that GHF models often do not honor the estimated constraints. Our results highlight the need for community effort to reconcile the discrepancies between radar data, GHF models, and ice core information.

We derive the surface and basal radar reflectance and backscatter coefficients of the southern McMurdo Ice Shelf (SMIS) and part of the nearby Ross Ice Shelf (RIS), Antarctica, from radar statistical reconnaissance using a 60-MHZ airborne survey. The surface coefficients are further inverted in terms of snow density and roughness, providing a spatial distribution of the processes contributing to the surface boundary conditions. We disentangle the basal coefficients from surface transmission losses, and we provide the basal coherent content, an indicator of the boundary geometric disorder that is also self-corrected from englacial attenuation. The basal radar properties exhibit sharp gradients along specific iso-depths, suggesting an abrupt modification of the ice composition and geometric structure. We interpret this behavior as locations where the pressure-melting point is reached, outlining fields of freezing and melting ice. Basal steps are observed at both SMIS and RIS, suggesting a common geometric expression of widespread basal processes. This technique offers a simultaneous view of both the surface and basal boundary conditions to help investigate the ice-shelf stability, while its application to airborne data significantly improves coverage of the difficult-to-observe ice–ocean boundary. It also provides constraints on thermohaline circulation in ice shelves cavities, which are analogs for ice-covered ocean worlds.

Glacier sliding over small obstacles relies on melting on their upstream sides and refreezing downstream. Previous treatments have appealed to ‘pressure melting’ as the cause of the spatial variations in melting temperature that drive this regelation process. However, we show that typical liquid pressure variations across small obstacles are negligible and therefore variations in ice pressure closely approximate variations in effective stress. For a given change in effective stress, the equilibrium melting temperature changes by an order of magnitude more than when the pressure of ice and liquid both change by an equal amount. In consequence, the temperature gradients that drive heat flow across small obstacles are larger than previously recognized and the rate of regelation is faster. Under typical conditions, the transition wavelength at which ice deformation and regelation contribute equally is of m-scale, several times longer than previous predictions, which have been reported to underestimate field inferences.

To study subglacial hydrological condition and its influence on the glacier dynamics, we drilled Johnsons Glacier on Livingston Island in the Antarctic Peninsula region. Subglacial water pressure was recorded in boreholes at two locations over 2 years, accompanied by high-frequency ice-speed measurements during two summer melt seasons. Water pressure showed two different regimes, namely high frequency and large amplitude variations during the melt season (January–April) and small fluctuations near the overburden pressure the rest of the year. Speed-up events were observed several times in each summer measurement period. Ice motion during these events substantially contributed to total glacier motion, for example, fast ice flow over 1 week accounted for ~70% of the total displacement over a 25-day long measurement period. We did not find a clear relationship between subglacial water pressure and ice speed. This was probably because subglacial hydraulic conditions were spatially inhomogeneous and thus our borehole data did not always represent a large-scale subglacial condition. Ice temperature measurements in the boreholes confirmed the existence of a cold ice layer near the glacier surface. Our data provide a basis to better understand the dynamic and hydrological conditions of relatively unstudied glaciers in the Antarctic Peninsula region.