Characterising the structures within glaciers can give unique insight into ice motion processes. On debris-covered glaciers, traditional structural glaciological mapping is challenging because the lower glacier is hidden by the supraglacial debris layer. Here, we use high-resolution optical televiewer (OPTV) image logs from four boreholes drilled into Khumbu Glacier, Nepal, to overcome this limitation and investigate englacial structural features within a Himalayan debris-covered glacier. The OPTV logs show structural features that are up to an order of magnitude thinner than those observed at the glacier surface and reveal five structural units: (I) primary stratification of ice; (II) debris-rich planes that conform with the primary stratification; (III) water-healed crevasse traces; (IV) healed crevasse traces; and (V) steeply dipping planes of basally derived fine sediment near the glacier terminus. The OPTV logs also reveal that the primary stratification both decreases in dip with depth (by up to 56° over 20 m) and rotates with depth (by up to 100° over 20 m) towards parallelism with the proximal lateral moraine. This transformation and the presence of relict layers of basally derived sediment raised into an englacial position – possibly involving thrusting – near the glacier's now stagnant terminus reveal a previously more dynamic glacier regime.

We assess the evolution of glaciological structures during the 2003–05 surge in the Paulabreen glacier system, Svalbard. Glaciological structures on the glacier surface were mapped using aerial photographs captured in the early stages of the surge (2003) and 5 years after surge termination (2011). Three-dimensional measurements of glaciological structures were collected at the tidewater front in 2013. These datasets document the physical changes during (1) the late quiescent phase; (2) the early phase of the surge as the surge front propagated down Skobreen and advanced into Paulabreen and (3) the final stages of the surge following the surge front reaching the glacier terminus. Crevasse patterns and clusters of arcuate shear planes record zones of compressive and extensional flow associated with the downglacier progression of the surge front. The transfer of surging ice from Skobreen into Paulabreen caused lateral displacement of the medial moraines to the northeast. At the ice front, this movement tilted glaciological structures in the same direction. Structures at the southwest margin record strike–slip faulting and the elevation of debris into the ice in a zone of compression and transpression. We summarise these observations in a schematic reconstruction of structural evolution during the surge.

Since the mid-20th century, ice shelves around the Antarctic Peninsula have declined in extent and thickness, and some have shown signs of structural instability. Here, using satellite imagery from 1999/2000 to 2019/20 (Landsat 7 and 8, Sentinel-2 and ASTER), we measure areal changes, calculate surface flow speeds, and quantify structural changes of Bach, Stange and George VI ice shelves, located in the southwest Antarctic Peninsula. We recorded a total area loss of 797.5 km2 from 2009/10 to 2019/20, though spatial and temporal patterns varied at individual ice fronts. Flow speeds remained largely stable over the observation periods, but notable acceleration was calculated for Bach Ice Shelf, and at the northern and southern extents of George VI Ice Shelf. Open fractures widened and lengthened between 2009/10 and 2019/20 on all three ice shelves. We conclude that Stange Ice Shelf is stable, and not under any immediate threat of enhanced recession. Continued ice-mass loss and consequential speed up of George VI South may cause further fracturing and destabilisation in the coming decades. Of more immediate concern are the glaciological changes noted for Bach Ice Shelf and George VI North; substantial areas of stabilising ice have already, or will soon be removed, that may lead to enhanced recession within the next decade.

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.

Microstructures, including crystallographic fabric, within the margin of streaming ice can exert strong control on flow dynamics. To characterize a natural setting, we retrieved three cores, two of which reached bed, from the flank of Jarvis Glacier, eastern Alaska Range, Alaska. The core sites lie ~1 km downstream of the source, with abundant water present in the extracted cores and at the base of the glacier. All cores exhibit dipping layers, a combination of debris bands and bubble-free domains. Grain sizes coarsen on average approaching the lateral margin. Crystallographic orientations are more clustered and with c-axes closer to horizontal nearer the lateral margin. The measured fabric is sufficiently weak to induce little mechanical anisotropy, but the data suggest that despite the challenging conditions of warm ice, abundant water and a short flow distance, many aspects of the microstructure, including measurable crystallographic fabric, evolved in systematic ways.

In this study, we explore small-scale (~1 to 20 km) thermal-refractive effects on basal geothermal heat flux (BGHF) at subglacial boundaries resulting from lateral thermal conductivity contrasts associated with subglacial topography and geologic contacts. We construct a series of two-dimensional, conductive, steady-state models that exclude many of the complexities of ice sheets in order to demonstrate the effect of thermal refraction. We show that heat can preferentially flow into or around a subglacial valley depending on the thermal conductivity contrast with underlying bedrock, with anomalies of local BGHF at the ice–bedrock interface between 80 and 120% of regional BGHF and temperature anomalies on the order of ±15% for the typical range of bedrock conductivities. In the absence of bed topography, subglacial contacts can produce significant heat flux and temperature anomalies that are locally extensive (>10 km). Thermal refraction can result in either an increase or decrease in the likelihood of melting and ice-sheet stability depending on the conductivity contrast and bed topography. While our models exclude many of the physical complexities of ice behavior, they illustrate the need to include refractive effects created by realistic geology into future glacial models to improve the prediction of subglacial melting and ice viscosity.

Ever since the first deep ice cores were drilled, it has been a challenge to determine their original, in-situ orientation. In general, the orientation of an ice core is lost as the drill is free to rotate during transport to the surface. For shallow ice cores, it is usually possible to match the adjacent core breaks, which preserves the orientation of the ice column. However, this method fails for deep ice cores, such as the EastGRIP ice core in Northeast Greenland. We provide a method to reconstruct ice core orientation using visual stratigraphy and borehole geometry. As the EastGRIP ice core is drilled through the Northeast Greenland Ice Stream, we use information about the directional structures to perform a full geographical re-orientation. We compared the core orientation with logging data from core break matching and the pattern of the stereographic projections of the crystals’ c-axis orientations. Both comparisons agree very well with the proposed orientation method. The method works well for 441 out of 451 samples from a depth of 1375–2120 m in the EastGRIP ice core. It can also be applied to other ice cores, providing a better foundation for interpreting physical properties and understanding the flow of ice.

We investigated the effects of microparticles and grain size on the microstructural evolutions and mechanical properties of polycrystalline ice. Uniaxial compression tests were conducted using fine-grained pure ice and silica-dispersed ice under various conditions. Deformation behavior of fine-grained ice was found to be characterized by stress exponent n ≈ 2 and activation energy Q ≈ 60 kJ mol−1. The derived strain rates of fine-grained ice were ≈ 1 order of magnitude larger than those of coarse-grained ice obtained in previous studies, and they were found to be independent of particle dispersion and dependent on the mean grain size of ice, with grain size exponent p ≈ 1.4. Work hardening was observed in dislocation creep, while the strain rate continued to decrease. These results indicate that the deformation mechanism of fine-grained ice is different from typical dislocation creep, often associated with n = 3. Although microparticles restricted grain growth, there was little direct effect on the deformation of fine-grained ice. Microstructural observations of the ice samples indicated that the grain boundaries were straight and that the subgrain boundary densities increased after deformation. Our experiments suggest that grain size and boundaries play important roles in the deformation processes of polycrystalline ice.

Over recent decades, the number of outlet glaciers terminating in lakes in Iceland has increased in line with climate warming. The mass-balance changes of these lake-terminating outlet glaciers are sensitive to rising air temperatures, due to altered glacier dynamics and increased surface melt. This study aims to better understand the relationship between proglacial lake development, climate, glacier dynamics and glacier structure at Fjallsjökull, a large, lake-terminating outlet glacier in south-east Iceland. We used satellite imagery to map glacier terminus position and lake extent between 1973 and 2016, and a combination of aerial and satellite imagery to map the structural architecture of the glacier's terminus in 1982, 1994 and 2011. The temporal evolution of ice surface velocities between 1990 and 2018 was calculated using feature tracking. Statistically significant increases in the rate of terminus retreat and lake expansion were identified in 2001, 2009 and 2011. Our surface velocity and structural datasets revealed the development of localised flow ‘corridors’ over time, which conveyed relatively faster flow towards the glacier's terminus. We attribute the overall changes in dynamics and structural architecture at Fjallsjökull to rising air temperatures, but argue that the spatial complexities are driven by glacier specific factors, such as basal topography.

This contribution discusses results obtained from 3-D neutron diffraction and 2-D fabric analyser in situ deformation experiments on laboratory-prepared polycrystalline deuterated ice and ice containing a second phase. The two-phase samples used in the experiments are composed of an ice matrix with (1) air bubbles, (2) rigid, rhombohedral-shaped calcite and (3) rheologically soft, platy graphite. Samples were tested at 10°C below the melting point of deuterated ice at ambient pressures, and two strain rates of 1 × 10−5 s−1 (fast) and 2.5 × 10−6 s−1 (medium). Nature and distribution of the second phase controlled the rheological behaviour of the ice by pinning grain boundary migration. Peak stresses increased with the presence of second-phase particles and during fast strain rate cycles. Ice-only samples exhibit well-developed crystallographic preferred orientations (CPOs) and dynamically recrystallized microstructures, typifying deformation via dislocation creep, where the CPO intensity is influenced in part by the strain rate. CPOs are accompanied by a concentration of [c]-axes in cones about the compression axis, coinciding with increasing activity of prismatic-<a> slip activity. Ice with second phases, deformed in a relatively slower strain rate regime, exhibit greater grain boundary migration and stronger CPO intensities than samples deformed at higher strain rates or strain rate cycles.

We analyse ice flow and structural glaciology of Shackleton Glacier, a dendritic glacier with multiple icefalls in the Canadian Rockies. A major tributary-trunk junction allows us to investigate the potential of tributaries to alter trunk flow and structure, and the formation of bedrock steps at confluences. Multi-year velocity-stake data and structural glaciology up-glacier from the junction were assimilated with glacier-wide velocity derived from Radarsat-2 speckle tracking. Maximum flow speeds are 65 m a−1 in the trunk and 175 m a−1 in icefalls. Field and remote-sensing velocities are in good agreement, except where velocity gradients are high. Although compression occurs in the trunk up-glacier of the tributary entrance, glacier flux is steady state because flow speed increases at the junction due to the funnelling of trunk ice towards an icefall related to a bedrock step. Drawing on a published erosion model, we relate the heights of the step and the hanging valley to the relative fluxes of the tributary and trunk. It is the first time that an extant glacier is used to test and support such model. Our study elucidates the inherent complexity of tributary/trunk interactions and provides a conceptual model for trunk flow restriction by a tributary in surge-type glaciers.

The surface of the Ross Ice Shelf (RIS) is textured by flow stripes, crevasses and other features related to ice flow and deformation. Here, moderate resolution optical satellite images are used to map and classify regions of the RIS characterized by different surface textures. Because the textures arise from ice deformation, the map is used to identify structural provinces with common deformation history. We classify four province types: regions associated with large outlet glaciers, shear zones, extension downstream of obstacles and suture zones between provinces with different upstream sources. Adjacent provinces with contrasting histories are in some locations deforming at different rates, suggesting that our province map is also an ice fabric map. Structural provinces have more complicated shapes in the part of the ice shelf fed by West Antarctic ice streams than in the part fed by outlet glaciers from the Transantarctic Mountains. The map may be used to infer past variations in stress conditions and flow events that cannot be inferred from flow traces alone.

Large numbers of small valley glaciers on Svalbard were thicker and more extensive during the Little Ice Age (LIA), demonstrated by prominent ice-cored moraines up to several kilometres beyond present-day margins. The majority of these glaciers have since experienced a long period of strongly negative mass balance during the 20th century and are now largely frozen to their beds, indicating they are likely to have undergone a thermal transition from a polythermal to a cold-based regime. We present evidence for such a switch by reconstructing the former flow dynamics and thermal regime of Tellbreen, a small cold-based valley glacier in central Spitsbergen, based on its basal sequence and glaciological structures. Within the basal sequence, the underlying matrix-supported diamict is interpreted as saturated subglacial traction till which has frozen at the bed, indicating that the thermal switch has resulted in a cessation of subglacial sediment deformation due to freezing of the former deforming layer. This is overlain by debris-poor dispersed facies ice, interpreted to have formed through strain-induced metamorphism of englacial ice. The sequential development of structures includes arcuate fracture traces, interpreted as shear planes formed in a compressive/transpressive stress regime; and fracture traces, interpreted as healed extensional crevasses. The formation of these sediment/ice facies and structures is indicative of dynamic, warm-based flow, most likely during the LIA when the glacier was significantly thicker.

Our objective is to map dynamic provinces and investigate dynamic changes in Jakobshavn Isbræ, Greenland. We use an approach that combines structural glaciology and remote-sensing data analysis, facilitated by mathematical characterization of generalized spatial surface roughness that provides parameters related to ice dynamics, deformation and interaction of the ice with bed topography. The approach is applied to derive time series of elevation and roughness changes and to attribute changes during rapid retreat. Different dynamic types of fast- and slow-moving ice can be mapped from ICESat Geoscience Laser Altimeter System data (2003–09) and Airborne Topographic Mapper data, using spatial roughness characterization, validated with ASTER and bed-topographic data. Results of comparative analysis of elevation changes and roughness changes of Jakobshavn south ice stream indicate (1) surface lowering of 10–15 m a-1 between 2004 and 2009 and (2) no change in surface roughness and dynamic types. These findings are consistent with a front retreat as part of a fjord-glacier cycle or following warming of fjord water and with climatic warming, but not with an internal dynamic acceleration as a cause of the observed changes during rapid retreat. Relationships to changes in basal water pressure are discussed. All glaciodynamic changes appear to have initiated near the front and propagated up-glacier.

To understand the dynamics of ice shelves, a knowledge of their internal and basal structure is very important. As the capacity to perform local surveys is limited, remote sensing provides an opportunity to obtain the relevant information. We must prove, however, that the relevant information can be obtained from remote sensing of the surface. That is the aim of this study. The Jelbart Ice Shelf, Antarctica, exhibits a variety of surface structures appearing as stripe-like features in radar imagery. We performed an airborne geophysical survey across these features and compared the results to TerraSAR-X imagery. We find that the stripe-like structures indicate surface troughs coinciding with the location of basal channels and crevasse-like features, revealed by radio-echo sounding. HH and VV polarizations do not show different magnitude. In surface troughs, the local accumulation rate is larger than at the flat surface. Viscoelastic modelling is used to gain an understanding of the surface undulations and their origin. The surface displacement, computed with a Maxwell model, matches the observed surface reasonably well. Our simulations show that the surface troughs develop over decadal to centennial timescales.

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.