Climate change has led to a significant shrinkage of glaciers in the Tropical Andes during the last decades. Recent multi-temporal quantifications of ice mass loss at mountain range to regional scale are missing. However, this is fundamental information for future water resource planning and glacier change projections. In this study, we measure temporally consistent glacier area changes and geodetic mass balances throughout the Bolivian Cordillera Real and Tres Cruces based on multi-sensor remote-sensing data. By analyzing multi-spectral satellite images and interferometric SAR data, a glacier recession of 81 ± 18 km2 (29%; 5.1 ± 1.1 km2 a−1), a geodetic mass balance of −403 ± 78 kg m−2 a−1 and a total ice mass loss of 1.8 ± 0.5 Gt is derived for 2000–2016. In the period 2013–2016, ice mass loss was 21% above the average rate. A retreat rate of 15 ± 5 km2 a−1 and a mass budget of −487 ± 349 kg m−2 a−1 are found in this more recent period. These higher change rates can be attributed to the strong El Niño event in 2015/16. The analyses of individual glacier changes and topographic variables confirmed the dependency of the mass budget and glacier recession on glacier aspect and median elevation.

Evolving conditions at the terminus of Thwaites Glacier will be important in determining the rate of its future sea-level contribution over the coming decades. Here, we use remote-sensing observations to investigate recent changes (2000–2018) in the structure and velocity of Thwaites Glacier and its floating tongue. We show that the main trunk of Thwaites Glacier has accelerated by 38% over this period, while its previously intact floating tongue has transitioned to a weaker mélange of fractured icebergs bounded by sea ice. However, the rate of structural weakening and acceleration was not uniform across the observational period and we identify two periods of rapid acceleration and structural weakening (2006–2012; 2016–2018), separated by a period of deceleration and re-advance of the structurally-intact shear margin boundary (2012–2015). The timing of these accelerations/decelerations strongly suggests a link to variable ocean forcing. The weakened tongue now has some dependency on landfast sea ice for structural integrity and is vulnerable to changes in landfast ice persistency. Future reductions in landfast sea ice could manifest from changes in climate and/or the imminent removal of the B-22A iceberg from the Thwaites embayment. Such changes could have important implications for the integrity of the ice tongue and future glacier discharge.

The representation of iceberg calving in numerical models is a key source of uncertainty in century-scale sea-level rise projections. Parameters central to model representations of calving, including the tensile strength of glacier ice, remain poorly constrained. Grain-size and sample-size dependence make it difficult to reconcile laboratory and in situ estimates of ice tensile strength. Further, assumptions of various numerical models obscure comparison of the ‘strength’ parameter with a physically observable value. Here, we address the problem of fracture during calving using an analogous natural laboratory: a viscoelastic analysis of observed surface deformation and associated stresses in the 2015 collapse of eastern Skaftá cauldron, Vatnajökull ice cap, Iceland. We find that the ice within the cauldron could have experienced instantaneous elastic stress on the order of several MPa. We observe surface crevasses at the cauldron edges and center, but find that large areas of ice remain intact despite high stress. Our findings suggest a tensile strength of glacier ice on the order of 1 MPa, consistent with laboratory estimates but exceeding previous glacier-specific estimates.

In this paper, we give an overview of changes in area, length, surface elevation and mass balance of glaciers in mainland Norway since the 1960s. Frontal advances have been recorded in all regions except the northernmost glaciers in Troms and Finnmark (Storsteinsfjellbreen, Lyngen and Langfjordjøkelen). More than half of the observed glaciers, 27 of 49, had marked advances in the 1990s. The glaciological mass-balance values for the period 1962–2018, where 43 glaciers have been measured, show great inter-annual variability. The results reveal accelerated deficit since 2000, the most negative decade being 2001–2010. Some years with a positive mass balance (or less negative) after 2010s can be attributed to variations in large-scale atmospheric circulation. A surface elevation change and geodetic mass balance were calculated for a sample of 131 glaciers covering 817 km2 in the ‘1960s’ and 734 km2 in the ‘2010s’, giving an area reduction of 84 km2, or 10%. The sample covers many of the largest glaciers in Norway, and they had an overall change in surface elevation of −15.5 m for the ~50 year period. Converted to a geodetic mass balance this gives a mean mass balance of −0.27 ± 0.05 m w.e. a−1.

The Antarctic Ice Sheet loses mass via its ice shelves predominantly through two processes: basal melting and iceberg calving. Iceberg calving is episodic and infrequent, and not well parameterized in ice-sheet models. Here, we investigate the impact of hydrostatic forces on calving. We develop two-dimensional elastic and viscous numerical frameworks to model the ‘footloose’ calving mechanism. This mechanism is triggered by submerged ice protrusions at the ice front, which induce unbalanced buoyancy forces that can lead to fracturing. We compare the results to identify the different roles that viscous and elastic deformations play in setting the rate and magnitude of calving events. Our results show that, although the bending stresses in both frameworks share some characteristics, their differences have important implications for modeling the calving process. In particular, the elastic model predicts that maximum stresses arise farther from the ice front than in the viscous model, leading to larger calving events. We also find that the elastic model would likely lead to more frequent events than the viscous one. Our work provides a theoretical framework for the development of a better understanding of the physical processes that govern glacier and ice-shelf calving cycles.

Climate-induced cryospheric changes can have a significant impact on the downstream water availability. In this study, the Open Global Glacier Model (OGGM) and the Glacio-hydrological Degree-day Model (GDM) are integrated to project the response of cryospheric and hydrological systems to climate change until 2100. The study area comprises six sub-basins of glacierized Koshi River basin covering Nepalese and Chinese territories. The output from OGGM is provided as input to GDM along with the spatial and hydro-meteorological data. The average glacier area change in all the sub-basins from 2021 to 2100 is estimated as 65 and 85% decrease and the average glacier volume change is estimated as 76 and 86% decrease for RCP 4.5 and 8.5 scenarios, respectively. The future simulated discharge shows an increasing trend in pre-monsoon and monsoon seasons and a decreasing trend in post-monsoon and winter seasons after 2060 in all the sub-basins, which can lead to wetter wet seasons and drier dry seasons in the far future. A shift in peak flow is observed from August to July in most of the sub-basins. The coupled modelling technique used in this study can largely improve our understanding of glacio-hydrological dynamics in the Himalayan region.

Using an ensemble of close- and long-range remote sensing, lake bathymetry and regional meteorological data, we present a detailed assessment of the geometric changes of El Morado Glacier in the Central Andes of Chile and its adjacent proglacial lake between 1932 and 2019. Overall, the results revealed a period of marked glacier down wasting, with a mean geodetic glacier mass balance of −0.39 ± 0.15 m w.e.a−1 observed for the entire glacier between 1955 and 2015 with an area loss of 40% between 1955 and 2019. We estimate an ice elevation change of −1.00 ± 0.17 m a−1 for the glacier tongue between 1932 and 2019. The increase in the ice thinning rates and area loss during the last decade is coincident with the severe drought in this region (2010–present), which our minimal surface mass-balance model is able to reproduce. As a result of the glacier changes observed, the proglacial lake increased in area substantially between 1955 and 2019, with bathymetry data suggesting a water volume of 3.6 million m3 in 2017. This study highlights the need for further monitoring of glacierised areas in the Central Andes. Such efforts would facilitate a better understanding of the downstream impacts of glacier downwasting.

On 30 September 2017, an Air France Airbus A380-800 suffered a failure of its fourth engine while over Greenland. This failure resulted in the loss of the engine fan hub, fan blades and surrounding structure. An initial search recovered 30 pieces of light debris, but the primary part of interest, a ~220 kg titanium fan hub, was not recovered because it had a different fall trajectory than the light debris, impacted into the ice-sheet's snow surface, and was quickly covered by drifting snow. Here we describe the methods used for the detection of the fan hub and details of the field campaigns. The search area included two crevasse fields of at least 50 snow-covered crevasses 1 to ~30 m wide with similar snow bridge thicknesses. After 21 months and six campaigns, using airborne synthetic aperture radar, ground-penetrating radar, transient electromagnetics and an autonomous vehicle to survey the crevasse fields, the fan hub was found within ~1 m of a crevasse at a depth of ~3.3 to 4 m and was excavated with shovels, chain saws, an electric winch, sleds and a gasoline heater, by workers using fall-arrest systems.

Glaciers in the Alps and several other regions in the world have experienced strong negative mass balances over the past few decades. Some of them are disappearing, undergoing exceptionally negative mass balances that impact the mean regional value, and require replacement. In this study, we analyse the geomorphometric characteristics of 46 mass-balance glaciers in the Alps and the long-term mass-balance time series for a subset of nine reference glaciers. We identify regime shifts in the mass-balance time series (when non-climatic controls started impacting) and develop a glacier vulnerability index (GVI) as a proxy for their possible future development, based on criteria such as hypsometric index, breaks in slope, thickness distribution and elevation change pattern. We found that the subset of 46 mass-balance glaciers reflects the characteristics of the total glacier sample very well and identified a region-specific variability of the mass balance. As the GVI is strongly related to cumulative glacier mass balances, it can be used as a pre-selector of future mass-balance glaciers. We conclude that measurements on rapidly shrinking glaciers should be continued as long as possible to identify regime shifts in hind-cast and better understand the impacts of climatic variability on such glaciers.

Accurate estimations of ice thickness and volume are indispensable for ice flow modelling, hydrological forecasts and sea-level rise projections. We present a new ice thickness estimation model based on a mass-conserving forward model and a Bayesian inversion scheme. The forward model calculates flux in an elevation-band flow-line model, and translates this into ice thickness and surface ice speed using a shallow ice formulation. Both ice thickness and speed are then extrapolated to the map plane. The model assimilates observations of ice thickness and speed using a Bayesian scheme implemented with a Markov chain Monte Carlo method, which calculates estimates of ice thickness and their error. We illustrate the model's capabilities by applying it to a mountain glacier, validate the model using 733 glaciers from four regions with ice thickness measurements, and demonstrate that the model can be used for large-scale studies by fitting it to over 30 000 glaciers from five regions. The results show that the model performs best when a few thickness observations are available; that the proposed scheme by which parameter-knowledge from a set of glaciers is transferred to others works but has room for improvements; and that the inferred regional ice volumes are consistent with recent estimates.

Rapidly-rising jökulhlaups, or glacial outburst floods, are a phenomenon with a high potential for damage. The initiation and propagation processes of a rapidly-rising jökulhlaup are still not fully understood. Seismic monitoring can contribute to an improved process understanding, but comprehensive long-term seismic monitoring campaigns capturing the dynamics of a rapidly-rising jökulhlaup have not been reported so far. To fill this gap, we installed a seismic network at the marginal, ice-dammed lake of the A.P. Olsen Ice Cap (APO) in NE-Greenland. Episodic outbursts from the lake cause flood waves in the Zackenberg river, characterized by a rapid discharge increase within a few hours. Our 6 months long seismic dataset comprises the whole fill-and-drain cycle of the ice-dammed lake in 2012 and includes one of the most destructive floods recorded so far for the Zackenberg river. Seismic event detection and localization reveals abundant surface crevassing and correlates with changes of the river discharge. Seismic interferometry suggests the existence of a thin basal sedimentary layer. We show that the ballistic part of the first surface waves can potentially be used to infer medium changes in both the ice body and the basal layer. Interpretation of time-lapse interferograms is challenged by a varying ambient noise source distribution.

Snow appears as a granular material in most engineering applications. We examined the role of grain shape and cohesion in angle of repose experiments, which are a common means for the characterization of granular materials. The role of shape was examined by investigating diverse snow types with discernable shape and spherical ice beads. Two geometrical shape parameters were calculated from X-ray micro-computed-tomography images after a particle segmentation was performed with a watershed algorithm. Cohesion was examined by conducting experiments at six different temperatures between −40 and −2°C, assuming sintering as its cause, which accelerates with increasing temperature. As a cohesionless reference, experiments with glass beads were performed. We found that both shape and cohesion exerted about equally strong influence on the angle of repose. We utilized our results for an empirical model that describes the influence of shape and cohesion as additive corrections of the angle of repose of cohesionless spheres and explains all experiments with a correlation coefficient r2 = 0.95. With temperature and the chosen shape parameter as fitting variables, previous experiments with another snow type could be consistently included. The experiments highlight the relevance of these parameters in granular snow mechanics and can be used for model calibration.

Geothermal heat flux (GHF) is an important control on the dynamics of Antarctica's ice sheet because it controls basal melt and internal deformation. However, it is hard to estimate because of a lack of in-situ measurements. Estimating GHF from ice-borehole temperature profiles is possible by combining a heat-transfer equation and the physical properties of the ice sheet in a numerical model. In this study, we truncate ice-borehole temperature profiles to determine the minimum ratio of temperature profile depth to ice-sheet thickness required to produce acceptable GHF estimations. For Law Dome, a temperature profile that is within 60% of the local ice thickness is sufficient for an estimation that is within approximately one median absolute deviation of the whole-profile GHF estimation. This result is compared with the temperature profiles at Dome Fuji and the West Antarctic Ice Sheet divide which require a temperature profile that is 80% and more than 91% of the ice thickness, respectively, for comparable accuracy. In deriving GHF median estimations from truncated temperature profiles, it is possible to discriminate between available GHF models. This is valuable for assessing and constraining future GHF models.

Several global datasets of glacier thickness exist, but the number of observations from western Canada are sparse and spatially biased. To supplement these limited observations, we measured ice thickness with ice penetrating radar on five glaciers in the Columbia Mountains, Canada. Our radar surveys, when combined with previous surveys for two glaciers in the Rocky Mountains, total 182 km of transects that represent 34 672 point measurements of ice thickness. Our measurements are, on average, 38% thicker than previous surface inversion model estimates of glacier thickness. Using our measurements within a cross-validation scheme, we model ice thickness with the Open Global Glacier Model (OGGM) driven with recent observations of surface mass balance and glacier elevation. We calibrated OGGM ice thickness by optimizing the ice creep parameter in the model. The optimized OGGM yields an ice volume for Columbia Basin of 122.5 ± 22.4 km3 for the year 2000, which is 23% greater than the range of previous estimates. At current rates of glacier mass loss for this region, glaciers would disappear from the basin in about 65–80 years. Disappearance of these glaciers will negatively affect the basin's surface hydrology, freshwater availability and aquatic ecosystems.

Passive seismology is becoming increasingly popular for glacier/ice-sheet structure investigations in Polar regions. Single-station passive seismic methods including P-wave receiver functions (PRFs), horizontal-to-vertical spectral ratio (HVSR) and a recently proposed autocorrelation method have been used to retrieve glacier/ice-sheet structures. Despite their successful applications, analysis regarding their detection abilities in different glaciological environments has not been reported. In this study, we compare ice thicknesses and vp/vs ratios obtained from the three methods using data collected at GAMSEIS and POLENET/ANET seismic arrays in Antarctica. Ice thickness estimates derived from the three methods are found to be consistent. Comparisons conducted under various model setups, including those involving tiled layers and sedimentary layers, show that the effectiveness of the autocorrelation method is not superior to the PRF method for retrieving ice-sheet structures. The autocorrelation method however can complement other methods as it only requires a single component seismic record.

Each summer, surface melting of the margin of the Greenland Ice Sheet exposes a distinctive visible stratigraphy that is related to past variability in subaerial dust deposition across the accumulation zone and subsequent ice flow toward the margin. Here we map this surface stratigraphy along the northern margin of the ice sheet using mosaicked Sentinel-2 multispectral satellite imagery from the end of the 2019 melt season and finer-resolution WorldView-2/3 imagery for smaller regions of interest. We trace three distinct transitions in apparent dust concentration and the top of a darker basal layer. The three dust transitions have been identified previously as representing late-Pleistocene climatic transitions, allowing us to develop a coarse margin chronostratigraphy for northern Greenland. Substantial folding of late-Pleistocene stratigraphy is observed but uncommon. The oldest conformal surface-exposed ice in northern Greenland is likely located adjacent to Warming Land and may be up to ~55 thousand years old. Basal ice is commonly exposed hundreds of metres from the ice margin and may indicate a widespread frozen basal thermal state. We conclude that the ice margin across northern Greenland offers multiple opportunities to recover paleoclimatically distinct ice relative to previously studied regions in southwestern Greenland.

Ice shelves play a critical role in modulating dynamic loss of ice from the grounded portion of the Antarctic Ice Sheet and its contribution to sea-level rise. Measurements of ice-shelf motion provide insights into processes modifying buttressing. Here we investigate the effect of seasonal variability of basal melting on ice flow of Ross Ice Shelf. Velocities were measured from November 2015 to December 2016 at 12 GPS stations deployed from the ice front to 430 km upstream. The flow-parallel velocity anomaly at each station, relative to the annual mean, was small during early austral summer (November–January), negative during February–April, and positive during austral winter (May–September). The maximum velocity anomaly reached several metres per year at most stations. We used a 2-D ice-sheet model of the RIS and its grounded tributaries to explore the seasonal response of the ice sheet to time-varying basal melt rates. We find that melt-rate response to changes in summer upper-ocean heating near the ice front will affect the future flow of RIS and its tributary glaciers. However, modelled seasonal flow variations from increased summer basal melting near the ice front are much smaller than observed, suggesting that other as-yet-unidentified seasonal processes are currently dominant.

The cross-sectional profile of a glacial valley can be obtained with a variational principle in which the friction against the valley walls and the glacier bed is extremized, subject to a Lagrangian constraint. We show that the actual valley profile maximizes the friction, thus settling an old debate.