The accumulation region of Fedchenko Glacier represents an extensive snow reservoir in the Pamir Mountains feeding the longest glacier in Central Asia. Observed elevation changes indicate a continuous ice loss in the ablation region of Fedchenko Glacier since 1928, while the mass balance of the accumulation region is largely unknown. In this study, we show that accumulation varies considerably in the main accumulation basin, with accumulation rates up to 2400 mm w.e. a−1 in the West, decreasing to <1000 mm w.e. a−1 in the center, although the elevation difference is <200 m. The combination of snow/firn samples and ground-penetrating radar profiles suggests that this accumulation pattern is persistent during the recent past. The recent accumulation history is reconstructed from internal radar reflectors using a firn densification model and shows strong interannual variations, but near constant mean values since 2002. Modeling of trajectories, based on accumulation and glacier geometry, results in an estimate of the depth/age relation close to the main divide. This region provides one of the most suitable locations for retrieving climate information with temporal high resolution for the last millennium, with a potential to cover most of the Holocene in less detail.

Negribreen, a tidewater glacier located in central eastern Svalbard, began actively surging after it experienced an initial collapse in summer 2016. The surge resulted in horizontal surface velocities of more than 25 m d−1, making it one of the fastest-flowing glaciers in the archipelago. The last surge of Negribreen likely occurred in the 1930s, but due to a long quiescent phase, investigations of this glacier have been limited. As Negribreen is part of the Negribreen Glacier System, one of the largest glacier systems in Svalbard, investigating its current surge event provides important information on surge behaviour among tidewater glaciers within the region. Here, we demonstrate the surge development and discuss triggering mechanisms using time series of digital elevation models (1969–2018), surface velocities (1995–2018), crevasse patterns and glacier extents from various data sources. We find that the active surge results from a four-stage process. Stage 1 (quiescent phase) involves a long-term, gradual geometry change due to high subglacial friction towards the terminus. These changes allow the onset of Stage 2, an accelerating frontal destabilization, which ultimately results in the collapse (Stage 3) and active surge (Stage 4).

The high-Alpine ice-core drilling site Colle Gnifetti (CG), Monte Rosa, Swiss/Italian Alps, provides climate records over the last millennium and beyond. However, the full exploitation of the oldest part of the existing ice cores requires complementary knowledge of the intricate glacio-meteorological settings, including glacier dynamics. Here, we present new ice-flow modeling studies of CG, focused on characterizing the flow at two neighboring drill sites in the eastern part of the glacier. The3-D full Stokes ice-flow model is thermo-mechanically coupled and includes firn rheology, firn densification and enthalpy transport, and is implemented using the finite element software Elmer/Ice. Measurements of surface velocities, accumulation, borehole inclination, density and englacial temperatures are used to validate the model output. We calculate backward trajectories and map the catchment areas. This constrains, for the first time at this site, the so-called upstream effects for the stable water isotope time series of the two ice cores drilled in 2005 and 2013. The model also provides a 3-D age field of the glacier and independent ice-core chronologies for five ice-core sites. Model results are a valuable addition to the existing glaciological and ice-core datasets. This especially concerns the quantitative estimate of upstream conditions affecting the interpretation of the deep ice-core layers.

We examine the spatial patterns of near-surface air temperature (Ta) over a melting glacier using a multi-annual dataset from McCall Glacier, Alaska. The dataset consists of a 10-year (2005–2014) meteorological record along the glacier centreline up to an upper glacier cirque, spanning an elevation difference of 900 m. We test the validity of on-glacier linear lapse rates, and a model that calculates Ta based on the influence of katabatic winds and other heat sources along the glacier flow line. During the coldest hours of each summer (10% of time), average lapse rates across the entire glacier range from −4.7 to −6.7°C km−1, with a strong relationship between Ta and elevation (R2 > 0.7). During warm conditions, Ta shows more complex, non-linear patterns that are better explained by the flow line-dependent model, reducing errors by up to 0.5°C compared with linear lapse rates, although more uncertainty might be associated with these observations due to occasionally poor sensor ventilation. We conclude that Ta spatial distribution can vary significantly from year to year, and from one glacier section to another. Importantly, extrapolations using linear lapse rates from the ablation zone might lead to large underestimations of Ta on the upper glacier areas.

We present changes of the ice margin in Northwest Greenland at the Eastern part of the Nunatarssuaq Ice Cap (NIC) over six decades. The ice margin in this area terminates as a near-vertical ice cliff of between 9 and 33 m thickness. During the years 1954–1957 and in 1965 multi-disciplinary studies were performed. We digitise and orthorectify material, that is often difficult to access, in order to use the historical data as an absolute starting point of our change assessment. We compare the cliff morphology of the mid-1950s and the mid-1960s with various time-steps between 1985 and 2017. The studied ice margin remained remarkably constant with very subtle changes of changing sign: rather slow advance rates are reported from the 1950s and 1960s that accelerated until 1985 and were followed by a general retreat until 2012 and a subsequent advance until 2017. Thickness changes are negative throughout the entire time-period, however, different rates of thinning are shown and there is a positive relationship with air temperature anomalies. Compared to similar elevations on the adjacent Greenland ice sheet, we find significantly weaker thinning rates at the NIC.

Glaciers and the periglacial environment in Argentina have been protected by the Law since 2010. This legislation required the development of the first National Glacier Inventory (NGI), which was officially presented in May 2018 and based on satellite images spanning between 2004 and 2016. Here, we present the methods and results of the NGI, summarize the glaciers’ morphological and spatial characteristics, and compare our results to previous regional and global inventories. The NGI reveals an impressive variety of ice masses including rock glaciers, permanent snowfields, mountain and valley glaciers with varying amounts of debris-cover and large outlet glaciers. The Argentinean Andes contain 16 078 ice masses covering an area of 5769 km2 between 200 and 6900 m a.s.l. Comparison of the combined national inventories of Argentina and Chile (~30 000 glaciers and 28 400 km2) with the Randolph Glacier Inventory 6.0 for the Southern Andes (~16 000 glaciers and 29 400 km2), shows that there are large differences in extent and number of glaciers in some sub-regions. The NGI represents an improvement for a better understanding of Argentina's freshwater reservoirs and provides detailed information for the preservation and study of ice masses along 4000 km of the Southern Andes.

Mount Achernar moraine is a terrestrial sediment archive that preserves a record of ice-sheet dynamics and climate over multiple glacial cycles. Similar records exist in other blue ice moraines elsewhere on the continent, but an understanding of how these moraines form is limited. We propose a model to explain the formation of extensive, coherent blue ice moraine sequences based on the integration of ground-penetrating radar (GPR) data with ice velocity and surface exposure ages. GPR transects (100 and 25 MHz) both perpendicular and parallel to moraine ridges at Mount Achernar reveal an internal structure defined by alternating relatively clean ice and steeply dipping debris bands extending to depth, and where visible, to the underlying bedrock surface. Sediment is carried to the surface from depth along these debris bands, and sublimates out of the ice, accumulating over time (>300 ka). The internal pattern of dipping reflectors, combined with increasing surface exposure ages, suggest sequential exposure of the sediment where ice and debris accretes laterally to form the moraine. Subsurface structure varies across the moraine and can be linked to changes in basal entrainment conditions. We speculate that higher concentrations of debris may have been entrained in the ice during colder glacial periods or entrained more proximal to the moraine sequence.

This work quantifies the increased temperature sensitivity of the constitutive behavior of ice with proximity to the melting point in terms of dislocation mechanics. An analysis of quasistatic and dynamic cyclic loading data for several ice types leads to the conclusion that high temperature (e.g. T ≥ −8°C) behavior is the result of a thermally induced increase in the number of mobile dislocations rather than an increase in the activation energy of dislocation glide or the introduction of a new deformation mechanism. The relationship between dislocation density and temperature is quantified and the model is shown to adequately represent the published minimum creep rate vs stress data for isotropic granular freshwater ice for −48 ≤ T ≤ −0.01°C.

Streaming ice accounts for a major fraction of global ice flux, yet we cannot yet fully explain the dominant controls on its kinematics. In this contribution, we use an anisotropic full-Stokes thermomechanical flow solver to characterize how mechanical anisotropy and temperature distribution affect ice flux. For the ice stream and glacier geometries we explored, we found that the ice flux increases 1–3% per °C temperature increase in the margin. Glaciers and ice streams with crystallographic fabric oriented approximately normal to the shear plane increase by comparable amounts: an otherwise isotropic ice stream containing a concentrated transverse single maximum fabric in the margin flows 15% faster than the reference case. Fabric and temperature variations independently impact ice flux, with slightly nonlinear interactions. We find that realistic variations in temperature and crystallographic fabric both affect ice flux to similar degrees, with the exact effect a function of the local fabric and temperature distributions. Given this sensitivity, direct field-based measurements and models incorporating additional factors, such as water content and temporal evolution, are essential for explaining and predicting streaming ice dynamics.

Diurnal depth cycles of decimeter scale are observed in a supraglacial lake on the McMurdo Ice Shelf, Antarctica. We evaluate two possible causes: (1) tidal tilt of the ice shelf in response to the underlying ocean tide, and (2) meltwater input variation. We find the latter to be the most likely explanation of our observations. However, we do not rule out tidal tilt as a source of centimeter scale variations, and point to the possibility that other, larger supraglacial lake systems, particularly those on ice shelves that experience higher amplitude tidal tilts, such as in the Weddell Sea, may have depth cycles driven by ocean tide. The broader significance of diurnal cycles in meltwater depth is that, under circumstances where the ice shelf is thin, tidal-tilt amplitudes are high, and meltwater runoff rates are large, there may be associated flexure stresses that can contribute to ice-shelf fracture and destabilization. For the McMurdo Ice Shelf (~20–50 m thickness, ~ 1 m tidal amplitude and ~10 cm water-depth variations), these stresses amount to several 10's of kPa.

Ice mélange has been postulated to impact glacier and fjord dynamics through a variety of mechanical and thermodynamic couplings. However, observations of these interactions are very limited. Here, we report on glaciological and oceanographic data that were collected from 2016 to 2017 at LeConte Glacier and Bay, Alaska, and serendipitously captured the formation, flow and break-up of ephemeral ice mélange. Sea ice formed overnight in mid-February. Over the subsequent week, the sea ice and icebergs were compacted by the advancing glacier terminus, after which the ice mélange flowed quasi-statically. The presence of ice mélange coincided with the lowest glacier velocities and frontal ablation rates in our record. In early April, increasing glacier runoff and the formation of a sub-ice-mélange plume began to melt and pull apart the ice mélange. The plume, outgoing tides and large calving events contributed to its break-up, which took place over a week and occurred in pulses. Unlike observations from elsewhere, the loss of ice mélange integrity did not coincide with the onset of seasonal glacier retreat. Our observations provide a challenge to ice mélange models aimed at quantifying the mechanical and thermodynamic couplings between ice mélange, glaciers and fjords.

Sources and implications of black carbon (BC) and mineral dust (MD) on two glaciers on the central Tibetan Plateau were estimated based on in situ measurements and modeling. The results indicated that BC and MD accounted for ~11 ± 1% and 4 ± 0% of the albedo reduction relative to clean snow, while the radiative forcing varied between 11 and 196 and 1–89 W m−2, respectively. Assessment of BC and MD contributions to the glacier melt can reach up 88 to 434 and 35 to 187 mm w.e., respectively, contributing ~9–23 and 4–10% of the total glacier melt. A footprint analysis indicated that BC and MD deposited on the glaciers originated mainly from the Middle East, Central Asia, North China and South Asia during the study period. Moreover, a potentially large fraction of BC may have originated from local and regional fossil fuel combustion. This study suggests that BC and MD will enhance glacier melt and provides a scientific basis for regional mitigation efforts.

The retention of meltwater in the accumulation area of the Greenland ice sheet and other Arctic ice masses buffers their contribution to sea level change. However, sustained warming also results in impermeable ice layers or ‘ice slabs’ that seal the underlying pore space. Here, we use a 1-D, physically based, high-resolution model to simulate the surface mass balance (SMB), percolation, refreezing, ice layer formation and runoff from across the high-elevation area of Devon Ice Cap, Canada, from 2001 to 2016. We vary the thickness of the ‘impermeable’ ice layer at which underlying firn becomes inaccessible to meltwater. Thick near-surface ice layers are established by an initial deep percolation, the formation of decimetre ice layers and the infilling of interleaving pore space. The cumulative SMB increases by 48% by varying impermeable layer thickness between 0.01 and 5 m. Within this range we identify narrower range (0.25–1 m) that can simulate both the temporal variability in SMB and the observed near-surface density structure. Across this range, cumulative SMB variation is limited to 6% and 45–49% of mass retention takes place within the annually replenished snowpack. Our results indicate cooler summers after intense mid-2000s warming have led to a partial replenishment of pore space.

Turbulent fluxes make a substantial and growing contribution to the energy balance of ice surfaces globally, but are poorly constrained owing to challenges in estimating the aerodynamic roughness length (z0). Here, we used structure from motion (SfM) photogrammetry and terrestrial laser scanning (TLS) surveys to make plot-scale 2-D and 3-D microtopographic estimations of z0 and upscale these to map z0 across an ablating mountain glacier. At plot scales, we found spatial variability in z0 estimates of over two orders of magnitude with unpredictable z0 trajectories, even when classified into ice surface types. TLS-derived surface roughness exhibited strong relationships with plot-scale SfM z0 estimates. At the glacier scale, a consistent increase in z0 of ~0.1 mm d−1 was observed. Space-for-time substitution based on time since surface ice was exposed by snow melt confirmed this gradual increase in z0 over 60 d. These measurements permit us to propose a scale-dependent temporal z0 evolution model where unpredictable variability at the plot scale gives way to more predictable changes of z0 at the glacier scale. This model provides a critical step towards deriving spatially and temporally distributed representations of z0 that are currently lacking in the parameterisation of distributed glacier surface energy balance models.

Computing predictions of future sea level that include well-defined uncertainty bounds requires models that are capable of robustly simulating the evolution of ice sheets and glaciers. Ice flow behaviour is known to be sensitive to the location and geometry of dynamic ice boundaries such as the grounding line (GRL), terminus position and ice surface elevation, so that any such model should track these interfaces with a high degree of accuracy. To address this challenge, we implement a numerical approach that uses the level-set method (LSM) that accurately models the evolution of the ice–air and ice–water interface as well as capturing topological changes in ice-sheet geometry. This approach is evaluated by comparing simulations of grounded and marine-terminating ice sheets to various analytical and numerical benchmark solutions. A particular advantage of the LSM is its ability to explicitly track the moving margin and GRL while using a fixed grid finite-difference scheme. Our results demonstrate that the LSM is an accurate and robust approach for tracking the ice surface interface and terminus for advancing and retreating ice sheets, including the transient marine ice-sheet interface and GRL positions.

Hydrofracturing can enhance the depth to which crevasses propagate and, in some cases, allow full depth crevasse penetration and iceberg detachment. However, many existing crevasse models either do not fully account for the stress field driving the hydrofracture process and/or treat glacier ice as elastic, neglecting the non-linear viscous rheology. Here, we present a non-local continuum poro-damage mechanics (CPDM) model for hydrofracturing and implement it within a full Stokes finite element formulation. We use the CPDM model to simulate the propagation of water-filled crevasses in idealized grounded glaciers, and compare crevasse depths predicted by this model with those from linear elastic fracture mechanics (LEFM) and zero stress models. We find that the CPDM model is in good agreement with the LEFM model for isolated crevasses and with the zero stress model for closely-spaced crevasses, until the glacier approaches buoyancy. When the glacier approaches buoyancy, we find that the CPDM model does not allow the propagation of water-filled crevasses due to the much smaller size of the tensile stress region concentrated near the crevasse tip. Our study suggests that the combination of non-linear viscous and damage processes in ice near the tip of a water-filled crevasse can alter calving outcomes.

Current sea-level rise partly stems from increased surface melting and meltwater runoff from the Greenland ice sheet. Multi-year snow, also known as firn, covers about 80% of the ice sheet and retains part of the surface meltwater. Since the firn cold content integrates its physical and thermal characteristics, it is a valuable tool for determining the meltwater-retention potential of firn. We use gap-filled climatological data from nine automatic weather stations in the ice-sheet accumulation area to drive a surface-energy-budget and firn model, validated against firn density and temperature observations, over the 1998–2017 period. Our results show a stable top 20 m firn cold content (CC20) at most sites. Only at the lower-elevation Dye-2 site did CC20 decrease, by 24% in 2012, before recovering to its original value by 2017. Heat conduction towards the surface is the main process feeding CC20 at all nine sites, while CC20 reduction occurs through low-cold-content fresh-snow addition at the surface during snowfall and latent-heat release when meltwater refreezes. Our simulations suggest that firn densification, while reducing pore space for meltwater retention, increases the firn cold content, enhances near-surface meltwater refreezing and potentially sets favourable conditions for ice-slab formation.

Deglaciation of high mountain rockwalls alters slope stability as rockwalls become more sensitive to modifications in environmental factors (e.g. seasonal temperature variations). In the past decades, increasing efforts focused on studying deglaciated Alpine rockwalls. Yet, currently deglaciating rockfaces remain unstudied. Here, we quantify surface area variations of massive ice bodies lying on high mountain rockwalls (ice aprons) in the French sector of the Mont Blanc massif between the end of the Little Ice Age (LIA) and 2018. Surface area estimates are computed from terrestrial and aerial oblique photographs via photogrammetry. This technique allows using photographs taken without scientific intent, and to tap into diverse historical or recent photographic catalogs. We derive an ice apron surface area model from precipitation records and the positive degree-days. The studied ice aprons shrank from 1854 to the 1950s, before expanding until the end of the 1990s. The beginning of the 21st century shows a decrease in surface area, leading to the complete melt of one of the studied ice aprons in 2017. Observed variations correlate with modeled surface area, suggesting strong sensitivity of ice aprons to changes in climatic variables. By studying site-specific correlations, we explore the importance of local drivers over the balance of ice aprons.

Glaciers depicted on old maps reveal their historical extents, before the advent of aerial and satellite remote sensing. Digital glacier inventories produced from these maps can be employed in assessments of centennial-scale glacier change. This study reconstructs the ~1899 (covering the period 1882–1916) glacier extent in Nordland, northern Norway, from historical gradteigskart maps, with an emphasis on examining the accuracy of the mapped glaciers. Glacier outlines were digitised from georectified scans of the analogue maps in a raster graphics editor and were subsequently inventoried in a GIS. The accuracy of the historical glacier extent was established from written descriptions and landscape photographs created during the original field surveys, and further validated against independent glacier outlines of (1) the maximum Little Ice Age extent derived from geomorphological evidence, and (2) the 1945 extent derived from vertical aerial photographs. An overall uncertainty of ±17% is associated with our inventory. Nordland's glaciers covered an area of 1712 ± 291 km2 in 1899. By 2000, total ice cover had decreased by 47% (807 ± 137 km2) at a rate of 6% 10 a−1 (80 ± 14 km2 10 a−1). The approach presented here may serve as a blueprint for future studies intending to derive glacier inventories from historical maps.

Seismic studies of glaciers yield insights into spatio-temporal processes within and beneath glaciers on scales relevant to flow and deformation of the ice. These methods enable direct monitoring of the bed in ways that complement other geophysical techniques, such as geodetic or ground penetrating radar observations. In this work, we report on the analysis of passive seismic data collected from the interior of the North East Greenland Ice Stream, the Greenland ice sheet's largest outlet glacier. We record thousands of basal earthquakes, many of which repeat with nearly identical waveforms. We also record many long-duration glacial tremor episodes that migrate across the seismic network with slow velocities (e.g. ~4–12 m s−1). Analysis of the basal earthquakes indicates a transition between times of individual event activity and times of tremor activity. We suggest that both processes are produced by shear slip at localized asperities along the bed. The transition between discrete and quasi-continuous slipping modes may be driven by pore-water pressure transients or heterogeneous strain accumulation in the ice due to strength contrasts of the underlying till.