We present a theoretical framework that integrates the dynamics of glaciers with and without the topographic confinement. This Part 1 paper concerns the former, which may exhibit surge cycles when subjected to thermal switches associated with the bed condition. With the topographic trough setting the glacier width and curbing the lateral drainage of the meltwater, the problem falls under the purview of the undrained plastic bed (UPB) formalism. Employing the UPB, we shall examine the external controls of the glacial behavior and test them against observations. Through our non-dimensionalization scheme, we construct a 2-D regime diagram, which allows a ready prognosis of the glacial properties over the full range of the external conditions, both climate- and size-related. We first discern the boundaries separating the glacial regimes of steady-creep, cyclic-surging and steady-sliding. We then apply the regime diagram to observed glaciers for quantitative comparisons. These include the Svalbard glaciers of both normal and surge types, Northeast Greenland Ice Stream characterized by steady-sliding, and Hudson Strait Ice Stream exhibiting cyclic surges. The quantitative validation of our model containing no free parameters suggests that the thermal switch may unify the dynamics of these diverse glaciers.

In Part 1, we have considered the dynamics of topographically confined glaciers, which may undergo surge cycles when the bed becomes temperate. In this Part 2, we consider the ice discharge over a flatbed, which would self-organize into alternating stream/ridge pairs of wet/frozen beds. The meltwater drainage, no longer curbed by the bed trough, would counter the conductive cooling to render a minimum bed strength at some intermediate width, toward which the stream would evolve over centennial timescale. At this stationary state, the stream width is roughly twice the geometric mean of the stream height and length, which is commensurate with its observed width. Over a flatbed, streams invariably interact, and we deduce that the neighboring ones would exhibit compensating cycles of maximum velocity and stagnation over the centennial timescale. This deduction is consistent with observed time variation of Ross ice streams B and C (ISB/C), which is thus a manifestation of the natural cycle. Moreover, the model uncovers an overlooked mechanism of the ISC stagnation: as ISB widens following its reactivation, it narrows ISC to augment the loss of the meltwater, leading to its stagnation. This stagnation is preceded by ice thickening hence opposite to the thinning-induced surge termination.

This paper presents a feasibility study of in situ field measurements of unsaturated meltwater percolation flux within the vertical profile of a snowpack, using the self-potential (SP) method. On-site snowmelt column tests calibrated the SP measurements. The SP data measured electrical field strength with an electrode spacing of 20 cm, and coincident water saturation (Sw) measurements using time domain reflectometry allowed calculation of SP-modeled vertical percolation flux (qsp), expressed as Darcy velocity. The results reflected transient diurnal snowmelt dynamics, with peak flux lagging arrival of a saturation wetting front. Peak daily qsp was 60 to >300 mm d−1, whereas daily snowmelt was 20–50 mm w.e. Surface refreezing events appeared to cause upward flow, possibly representing water redistribution toward the freezing boundary. Calculated fluxes were comparable to actual fluxes, although average errors ranged from −15 to +46% compared to average of melt expected from surface energy-balance and ablation stake measurements. By advancing method development to measure unsaturated meltwater percolation flux in snowpacks this study creates opportunities to study fundamental snowmelt processes, may improve mathematical modeling and may supplement glacier mass-balance studies and studies of snowmelt interactions with avalanches, groundwater and surface water.

This study deploys RTK-GNSS in 2012, TLS in 2015 and UAV in 2018 to monitor the changes of Urumqi Glacier No. 1 (UG1), eastern Tien Shan, and analyzes the feasibility of three technologies in monitoring the mountain glaciers. DEM differencing shows that UG1 has experienced a pronounced thinning and mass loss for the period of 2012–18. The glacier surface elevation change of −0.83 ± 0.57 m w.e. a−1 has been recorded for 2012–15, whereas the changes of glacier tongue surface elevation in 2015–18 and 2012–18 were −2.03 ± 0.95 and −1.34 ± 0.88 m w.e. a−1, respectively. The glacier area shrunk by 0.07 ± 0.07 × 10−3 km2 and the terminus retreat rate was 6.28 ± 0.83 m a−1 during 2012–18. The good agreement between the glaciological and geodetic specific mass-balances is promising, showing the combination of the three technologies is suitable to monitor glacier mass change. We recommend application of the three technologies to assess each other in different locations of the glacier, e.g. RTK-GNSS base stations, ground control points, glacier tongue and terminus, in order to avoid the inherent limitations of each technology and to provide reliable data for the future studies of mountain glacier changes in western China.

Recent regional cooling has impacted the natural systems of the Antarctic Peninsula (AP); however, little is known concerning the changes in the high parts of the glacial systems. Dry-snow line (DSL), situated in the high parts of glaciers, is the uppermost limit of frequent or occasional surface melt. We analyse dry-snow line altitude (DSLA) changes on the AP during 2004–2020 using C-band synthetic aperture radar time series data. We demonstrate that the DSLA in the eastern part of the AP is usually higher than that of the western part. Moreover, using simulated climatic variables from regional climate models, the lowering in altitude of DSL of glaciers in most areas is identified as a response to a decrease in snowmelt and an increase in precipitation. Furthermore, correlation analyses between simulated climatic variables and the DSLA are conducted. These results present the sensitive response of variations in DSLA to meteorological conditions, and the capability of DSLA being a proxy of polar local climate in high-altitude areas with no in situ meteorological observations.

Polycrystalline ice is known to exhibit macroscopic anisotropy in relative permittivity (ɛ) depending on the crystal orientation fabric (COF). Using a new system designed to measure the tensorial components of ɛ, we investigated the dielectric anisotropy (Δɛ) of a deep ice core sample obtained from Dome Fuji, East Antarctica. This technique permits the continuous nondestructive assessment of the COF in thick ice sections. Measurements of vertical prism sections along the core showed that the Δɛ values in the vertical direction increased with increasing depth, supporting previous findings of c-axis clustering around the vertical direction. Analyses of horizontal disk sections demonstrated that the magnitude of Δɛ in the horizontal plane was 10–15% of that in the vertical plane. In addition, the directions of the principal axes of tensorial ɛ in the horizontal plane corresponded to the long or short axis of the elliptically elongated single-pole maximum COF. The data confirmed that Δɛ in the vertical and horizontal planes adequately indicated the preferred orientations of the c-axes, and that Δɛ can be considered to represent a direct substitute for the normalized COF eigenvalues. This new method could be extremely useful as a means of investigating continuous and depth-dependent variations in COF.

The influence of surface melt on the flow of Greenland's largest outlet glaciers remains poorly known and in situ observations are few. We use field observations to link surface meltwater forcing to glacier-wide diurnal velocity variations on East Greenland's Helheim Glacier over two summer melt seasons. We observe diurnal variations in glacier speed that peak ~6.5 h after daily maximum insolation and extend from the terminus region to the equilibrium line. Both the amplitude of the diurnal speed variation and its sensitivity to daily melt are largest at the glacier terminus and decrease up-glacier, suggesting that the magnitude of the response is controlled not only by melt input volume and temporal variability, but also by background effective pressure, which approaches zero at the terminus. Our results provide evidence that basal lubrication by meltwater drives diurnal velocity variations at Greenland's marine-terminating glaciers in a similar manner to alpine glaciers and Greenland's land-terminating outlet glaciers.

Dissolved and particulate sodium, magnesium and calcium are analyzed in ice cores to determine past changes in sea ice extent, terrestrial dust variability and atmospheric aerosol transport efficiency. They are also used to date ice cores if annual layers are visible. Multiple methods have been developed to analyze these important compounds in ice cores. Continuous flow analysis (CFA) is implemented with instruments that sample the meltstream continuously. In this study, CFA with ICP-MS (inductively coupled-plasma mass spectrometry) and fast ion chromatography (FIC) methods are compared for analysis of sodium and magnesium. ICP-MS, FIC and fluorescence methods are compared for analysis of calcium. Respective analysis of a 10 m section of the Antarctic WACSWAIN Skytrain Ice Rise ice core shows that all of the methods result in similar levels of the compounds. The ICP-MS method is the most suitable for analysis of the Skytrain ice core due to its superior precision (relative standard deviation: 1.6% for Na, 1.3% for Mg and 1.2% for Ca) and sampling frequency compared to the FIC method. The fluorescence detection method may be preferred for calcium analysis due to its higher depth resolution (1.4 cm) relative to the ICP-MS and FIC methods (~4 cm).

Supraglacial ponds and ice cliffs can dramatically enhance ablation rates on debris-covered glaciers. Supraglacial ponds can also coalesce, forming moraine-dammed lakes at risk of glacial lake outburst flood (GLOF). Given Bhutanese glaciers have some of the highest ice loss rates in the Himalaya and GLOF vulnerability is high, we seek to advance our understanding of the spatial distribution and evolution of supraglacial ponds and ice cliffs. Here, we use high-resolution (3 m) Planet Labs satellite imagery to provide the first short-term, high-resolution dataset of supraglacial pond and ice cliff evolution for three glaciers along the Bhutan–Tibet border from 2016 to 2018. A total of 5754 ponds and 2088 ice cliffs were identified. Large intra-annual changes were observed, with ponded area changes and drainage events coinciding with the seasonality of the Indian Summer Monsoon. On average, ~19% of the total number of ponds had a coincident ice cliff. Pond spatial distribution was driven by ice-surface velocities, with higher numbers of ponds found in areas of low velocity (<8 m a−1). Our study provides the first detailed, quantitative investigation of supraglacial ponds and ice cliffs in Bhutan, providing a framework for further monitoring in this understudied, yet important, region of the Himalaya.

The Vadret da Tschierva (Vd Tschierva) is a 4 km long glacier in the Swiss Alps spanning an altitude range of 2400–4049 m a.s.l. Length observations since 1855 show steady retreat interrupted by a period of advance from 1965 until 1985. The total retreat is ~2200 m (period 1855–2018). We have studied the Vd Tschierva with a flowline model, combined with ‘buckets’ that represent steep hanging glaciers and ice-free rock faces delivering mass to the main stream. The model is calibrated by a control method, in which an ELA history is objectively determined by finding the best match between observed and simulated glacier length. There is a modest correlation between the reconstructed ELA and an ELA record based on meteorological observations at Segl-Maria (only 8 km away from the glacier). It is difficult to reproduce the observed length record when the glacier model is driven by climate model output (Coupled Model Intercomparison Project 5). We have calculated the future evolution of the Vd Tschierva for different rates of ELA rise. For a constant rise of 4 ${\rm m\;}{\rm a}^{ \hbox{-} 1}$, we predict that the glacier length will change from the current 3.2 km to ~1.7 km in the year 2100.

This paper examines the effect of basal topography and strength on the grounding-line position, flux and stability of rapidly-sliding ice streams. It does so by supposing that the buoyancy of the ice stream is small, and of the same order as the longitudinal stress gradient. Making this scaling assumption makes the role of the basal gradient and accumulation rate explicit in the lowest order expression for the ice flux at the grounding line and also provides the transcendental equation for the grounding-line position. It also introduces into the stability condition terms in the basal curvature and accumulation-rate gradient. These expressions revert to well-established expressions in circumstances in which the thickness gradient is large at the grounding line, a result which is shown to be the consequence of the non-linearity of the flow. The behaviour of the grounding-line flux is illustrated for a range of bed topographies and strengths. We show that, when bed topography at a horizontal scale of several tens of ice thicknesses is present, the grounding-line flux and stability have more complex dependencies on bed gradient than that associated with the ‘marine ice-sheet instability hypothesis’, and that unstable grounding-line positions can occur on prograde beds as well as stable positions on retrograde beds.

Firn compaction models inform mass-balance estimates and paleo-climate reconstructions, but current models introduce key uncertainties. For example, models disagree on the dependence of density and compaction on accumulation rate. Observations of compaction to test these models are rare, partly because in situ methods for measuring englacial strain are time-consuming and expensive. Moreover, shallow measurements may confound strain due to compaction with strain due to ice-sheet flow. We show that phase-sensitive radio-echo sounder (pRES) systems, typically deployed to measure sub-shelf melting or ice-sheet deformation, can be used to measure firn compaction and test firn models. We present two complementary methods for extracting compaction information from pRES data, along with a method for comparing compaction models to pRES observations. The methods make different assumptions about the density structure and vary in their need for independent density measurements. Compaction profiles computed from pRES data collected on three ice rises in West Antarctica are largely consistent with measured densities and a physics-based model. With their minimal logistic requirements, new pRES systems, such as autonomous pRES, could be inexpensively deployed to monitor firn compaction more widely. Existing phase-sensitive radar data may contain compaction information even when surveys targeted other processes.

Arctic landfast sea ice (LFSI) represents an important quasi-stationary coastal zone. Its evolution is determined by the regional climate and bathymetry. This study investigated the seasonal cycle and interannual variations of LFSI along the northwest coast of Kotelny Island. Initial freezing, rapid ice formation, stable and decay stages were identified in the seasonal cycle based on application of the visual inspection approach (VIA) to MODIS/Envisat imagery and results from a thermodynamic snow/ice model. The modeled annual maximum ice thickness in 1995–2014 was 2.02 ± 0.12 m showing a trend of −0.13 m decade−1. Shortened ice season length (−22 d decade−1) from model results associated with substantial spring (2.3°C decade−1) and fall (1.9°C decade−1) warming. LFSI break-up resulted from combined fracturing and melting, and the local spatiotemporal patterns of break-up were associated with the irregular bathymetry. Melting dominated the LFSI break-up in the nearshore sheltered area, and the ice thickness decreased to an average of 0.50 m before the LFSI disappeared. For the LFSI adjacent to drift ice, fracturing was the dominant process and the average ice thickness was 1.56 m at the occurrence of the fracturing. The LFSI stages detected by VIA were supported by the model results.

The long-term refrigerated storage of melted snow and/or ice samples for analyses of insoluble microparticles (hereafter, microparticles) may be limited by increases in the biological particle concentration caused by microbial growth after ~1–2 weeks. In this study, we examined an ultraviolet (UV) disinfection method for the storage of melted snow and/or ice samples and determined the effects of this method on microparticles. Surface snow obtained from Glacier No. 31 in the Suntar-Khayata Range, eastern Siberia, Russia was divided into two portions for UV treatment and untreated controls. Microparticle concentrations and size distributions (in the range of 0.52–12.0 μm) in the samples were measured using a Coulter counter. Whereas the microparticle concentration in untreated samples increased, no obvious increase was observed over 53 d in the samples subjected to UV treatment. Microbial growth was detected in only untreated samples using a viable particle counter. In addition, the original microparticle concentrations and size distributions were unaffected by UV treatment. Our results demonstrated that the microparticle size distribution in untreated melted water samples reflects the growth, decomposition and succession of microorganisms over time and further indicate that UV irradiation is effective for long-term storage for microparticle analysis.

The coastal Droning Maud Land in East Antarctica is characterized by small ice shelves with numbers of promontories and locally grounded isles, both called ice rises. These ice rises are typically dome-shaped and surface elevations are hundreds of meters above the surrounding ice shelves, which cause strong orographic effects on surface mass balance (SMB). We conducted shallow ice-penetrating radar sounding to visualize firn stratigraphy in the top 35 m over ~400 km of profiles across the Nivlisen Ice Shelf, and in a grid pattern over two adjacent ice rises (Djupranen and Leningradkollen). We tracked six reflectors (isochrones) and dated them using two ice cores taken at the ice rise summits, from which SMB over six periods in the past three decades was retrieved. The overall SMB pattern across the ice shelf remained similar for all periods; however, the eastwest contrast in SMB varies by a factor of 1.5–2 between the Leningradkollen and Djupranen grounding lines. The SMB patterns over the ice rises are more varied owing to complex interactions between topography, snowfall and wind. We use our results to evaluate the regional climate model RACMO2.3p2 in terms of the spatial SMB distribution and temporal changes over the ice shelf and ice rises at regional scale.

Surface speeds in Greenland's ablation zone undergo substantial variability on an annual basis which are presumed to mainly be driven by changes in sliding. Yet, meltwater-forced changes in ice–bed coupling can also produce variable deformation motion, which impacts the magnitude of sliding changes inferred from surface measurements and provides important context to flow dynamics. We examine spatiotemporal changes in deformation, sliding and surface velocities over a 2-year period using GPS and a dense network of inclinometers installed in borehole grid drilled in western Greenland's ablation zone. We find time variations in deformation motion track sliding changes through the summer and entire measurement period. A distinct spatial deformation and sliding pattern is also observed within the borehole grid which remains similar during winter and summer flow. We suggest that positively covarying sliding and deformation across seasonal timescales is characteristic of passive areas that are coupled to regions undergoing transient forcing, and the spatial patterns are consistent with variations in the local bed topography. The covarying deformation and sliding result in a 1.5–17% overestimate of sliding changes during summer compared to that inferred from surface velocity changes alone. This suggests that summer sliding increases are likely overestimated in many locations across Greenland.

The importance of glacier sliding has motivated a rich literature describing the thermomechanical interactions between ice, liquid water and bed materials. Early recognition of the gradient in melting temperature across small bed obstacles led to focused studies of regelation. An appreciation for the limits on ice deformation rates downstream of larger obstacles highlighted a role for cavitation, which has subsequently gained prominence in descriptions of subglacial drainage. Here, we show that the changes in melting temperature that accompany changes in normal stress along a sliding ice interface near cavities and other macroscopic drainage elements cause appreciable supercooling and basal mass exchange. This provides the basis of a novel formation mechanism for widely observed laminated debris-rich basal ice layers.

In the ablation zone of glacier habitats, cryoconite holes are known to harbor diverse microbial communities, including unique diatom floras distinct from those of surrounding aquatic and terrestrial systems. Besides descriptive studies, little is known about the diversity of cryoconite diatoms and their response to environmental stressors, particularly in low-latitude glaciers. This paper documents an extremely diversified diatom community in Antisana Glacier (Ecuador), reporting 278 taxa found in 54 surface holes, although with low individual abundances. Contrary to our expectations, assemblage structure did not respond to water physical or chemical characteristics, nor to cryoconite hole morphology, but to elevation. We demonstrate that elevation is a driver of diatom assemblages. Both alpha diversity (measured as Fisher's index) and species richness (corrected for unequal sample sizes) correlated negatively with elevation, suggesting a replacement toward simplified, poorer communities along this gradient. The taxonomic composition also changed significantly, as revealed by multivariate statistics. In summary, cryoconite holes are sites of high taxonomic diversity composed of taxa that are allochthonous in origin.