We present a quantitative ground-penetrating radar (GPR) analysis of the Lower Calderone Glacieret to highlight the recent evolution of one of the southernmost glacial systems in Europe. The Upper and Lower Calderone Glacierets are the last two perennial ice bodies in the Apennine Mountains (Italy), and their continuous monitoring is important for glaciology, hydrology and climate science. We applied a proprietary auto-picking algorithm to track reflections accurately and objectively within three pseudo-3-D GPR data sets that were acquired over the Lower Glacieret in different years. After the time-to-depth conversion, the basal reflections were projected onto the normal versors of a 3-D topographic model of the survey area, at the different GPR trace positions. We then applied an Adjusted Inverse-Distance Weighted Spatial Interpolation method to extrapolate the ice-bedrock interface within the areas not directly covered by the GPR profiles and compare it with the topographic surface to recover the glacieret volume. In this paper, we critically examine the accuracy of the reconstructed models, including possible issues related to the challenging survey areas, such as local artifacts in the interpolated interface caused by irregular GPR coverage. We further discuss the various advantages of the implemented procedure with respect to more traditional glacier monitoring techniques.

Surface height changes above three previously undetected subglacial lakes in northeastern Greenland are documented using CryoSat, DEMs and ICESat-2. Between 7 February and 6 March 2012, the central ice region (22.6 km2) above the largest lake dropped by ~37 m followed by a further drop of 12 m in the following 29 days. This implies a subglacial water outflow, or jökulhlaup, of at least 1 km3 at rates of hundreds of cubic meters per second. A comparable outflow occurred again between 23 July and 15 September 2019, with smaller outflows in the fall of 2014 and 2016. In contrast, a second smaller subglacial lake at a higher elevation had two subglacial outbursts of ~0.3 km3 in 2012 and 2019 but the lake filling was gradual and not strongly seasonal or episodic. Water remained in both lakes after the outflows but this may not be the case for the third smallest and lowest subglacial lake. While there appears to be some hydrological link between the three lakes, the flux of water moving under the ice in this area appears to be larger than would be expected from local summer melt. However, the source of the excess water remains uncertain.

Recent studies have reported a shift in the Arctic sea ice to a younger state after around 2007. This study reveals that this shift can be primarily attributed to a stepwise-type reduction in the extent of 4 years or older (4+ year) ice and its linked survivability. After this shift, the fraction of 4+ year ice extent relative to the total ice changed from 30.5 to 10.0%. Sea-ice survivability can serve as a key indicator of sea-ice persistence in response to other factors. We demonstrate that the decrease of 4+ year ice is controlled by the decrease of its linked survivability in a non-linear manner, signifying small alterations in the survivability can result in relatively large changes in the extent of 4+ year ice. The decrease in survivability is affected by both winter and summer processes. Summer melting contributed the most, while the contribution of the export through Fram Strait was minor. However, the significant rise in residual loss during the growth season suggests that other winter processes may also have played an important role.

We demonstrate that a ~20 km long valley glacier in the St. Elias Mountains, Yukon, can experience both partial and full surges, likely controlled by the presence of a topographic constriction and the formation and drainage of supraglacial lakes. Based on analysis of air photos, satellite images and field observations since the 1940s, we identify a full surge of ‘Little Kluane Glacier’ from 2013 to 2018, and a partial surge of just the upper north arm between 1963 and 1972. Repeat digital elevation models and velocity profiles indicate that the recent surge initiated from the upper north arm in 2013, which developed into a full surge of the main trunk from 2017 to 2018 with peak velocities of ~3600 m a−1 and frontal advance of ~1.7 km from May to September 2018. In 2016, a mass movement from the north arm to the main trunk generated a surface depression in a region immediately downstream of a topographic constriction, which promoted the formation and rapid drainage of supraglacial lakes to the glacier bed, and likely established the conditions to propel the initial partial surge into a full surge. Our results underscore the complex interplay between glacier geometry, surface hydrology and topography required to drive full surges of this glacier.

Ice-contact proglacial lakes affect ice dynamics and the transition of glacier termini from land- to lake-terminating has been shown to cause ice flow acceleration. In recent decades, the number and size of Greenlandic ice-marginal lakes has increased, highlighting the need to further understand these lake-terminating ice-margins as their influence on ice sheet mass balance increases. Here, time series of satellite-derived observations of ice velocity, surface elevation, and terminus position were generated at a lake-terminating outlet glacier, Isortuarsuup Sermia, and the nearby land-terminating Kangaasarsuup Sermia in south-west Greenland. At Isortuarsuup Sermia, annual surface velocity at the terminus increased by a factor of 2.5 to 214 ± 4 m yr−1 (2013–2021), with the magnitude of this acceleration declining with distance up-glacier. Meanwhile, near-terminus surface elevation changed at a rate of −2.3 ± 1.1 m yr−1 (2012–2021). Conversely, velocity change at Kangaasarsuup Sermia was minimal, while surface elevation change was approximately half at comparable elevations (−1.2 ± 0.3 m yr−1). We attribute these dynamic differences to thinning at Isortuarsuup Sermia and subsequent retreat from a stabilising sublacustrine moraine, and emphasise the potential of proglacial lakes to enhance future rates of mass loss from the Greenland Ice Sheet.

The objective of this research is to model snow's response to dynamic, impact loading. Two constitutive relationships are considered: elastic and Maxwell-viscoelastic. These material models are applied to laboratory experiments consisting of 1000 individual impacts across 22 snow column configurations. The columns are 60 cm tall with a 30 cm by 30 cm cross-section. The snow ranges in density from 135 to 428 kg m−3 and is loaded with both short-duration (~1 ms) and long-duration (~10 ms) impacts. The Maxwell-viscoelastic model more accurately describes snow's response because it contains a mechanism for energy dissipation, which the elastic model does not. Furthermore, the ascertained model parameters show a clear dependence on impact duration; shorter duration impacts resulted in higher wave speeds and greater damping coefficients. The stress wave's magnitude is amplified when it hits a stiffer material because of the positive interference between incident and reflected waves. This phenomenon is observed in the laboratory and modeled with the governing equations.

Glacier and ice sheet research is frequently justified on the basis of potential benefits to those communities that are most vulnerable to glacier change. In this glaciology research, funding priorities and communication to the broader public are strongly affected by the experiences and values of glaciology researchers. Using population data and newly available survey data from research organizations including glaciologists, we show that there is a substantial misalignment between the demographics of those who stand to benefit from glaciological research and those who produce glaciological knowledge. We discuss the potential negative consequences of this misalignment, which causes scientific research to be less effective, valuable and usable for communities. We conclude by outlining twenty evidence-based strategies that individuals and organizations can adopt to improve the recruitment and retention of a more representative group of scientists in glaciological research and encourage co-production with communities.

This study investigates the geodetic mass balance of nearly all glaciers in the Ladakh region, which are crucial for local water security. Utilizing multiple digital elevation models from 2000 and 2021, we estimate glacier mass balances. Climatic drivers of glacier mass balances are explored using ERA5-Land reanalysis data, evaluated by in situ climate data. The study also examines the role of nonclimatic (morphological) variables on glacier mass balances. Results indicate Ladakh glaciers experienced negative mass balances during 2000–2021, with significant spatial variability. Western Ladakh glaciers lost slightly higher mass (−0.35 ± 0.07 to −0.37 ± 0.07 m w.e. a−1) than eastern Ladakh glaciers (−0.21 ± 0.07 to −0.33 ± 0.05 m w.e. a−1). While warming is the main driver of widespread mass loss in Ladakh, the spatial variability in mass loss is attributed to changes in regional precipitation and glacier morphological settings. Eastern Ladakh glaciers, being smaller and at higher elevations, experience lower mass loss, whereas western Ladakh glaciers, larger and at lower elevations, are more susceptible to the impact of temperature, resulting in higher mass loss. The study underscores the potentially greater vulnerability of western Ladakh glaciers to a warming climate compared to their eastern counterparts.

We present a new map of bed topography and ice thickness together with a corresponding ice volume estimate representative of the years ~2010 for all Scandinavian ice caps and glaciers. Starting from surface observations, we invert for ice thickness by iteratively running an innovative ice dynamics model on a distributed grid and updating bed topography until modelled and observed glacier dynamics as represented by their rate of surface elevation change (dh/dt) fields align. The ice flow model used is the instructed glacier model (Jouvet and Cordonnier, 2023, Journal of Glaciology 1–15), a generic physics-informed deep-learning emulator that models higher-order ice flow with high-computational efficiency. We calibrate the modelled thicknesses against >11 000 ice thickness observations, resulting in a final ice volume estimate of 302.7 km3 for Norway, 18.4 km3 for Sweden and 321.1 km3 for the whole of Scandinavia with an error estimate of ~$\pm 11\%$. The validation statistics computed indicate good agreement between modelled and observed thicknesses (RMSE = 55 m, Pearson's r = 0.87, bias = 0.8 m), outperforming all other ice thickness maps available for the region. The modelled bed shapes thus provide unprecedented detail in the subglacial topography, especially for ice caps where we produce the first maps that show ice-dynamically realistic flow features.

Basal channels, which are troughs carved into the undersides of ice shelves by buoyant plumes of water, are modulators of ice-shelf basal melt and structural stability. In this study, we track the evolution of 12 large basal channels beneath ice shelves of the Amundsen and Bellingshausen seas region in West Antarctica using the Landsat record since its start in the 1970s through 2020. We observe examples of channel growth, interactions with ice-shelf features, and systematic changes in sinuosity that give insight into the life cycles of basal channels. We use the last two decades of the record, combined with contemporary ice-flow velocity datasets, to separate channel-path evolution into components related to advection by ice flow and those controlled by other forcings, such as ocean melt or surface accumulation. Our results show that ice-flow-independent lateral channel migration is overwhelmingly to the left when viewed down-flow, suggesting that it is dominated by Coriolis-influenced ocean melt. By applying a model of channel-path evolution dominantly controlled by ice flow and ocean melt, we show that the majority of channels surveyed exhibit non-steady behavior that serves as a novel proxy for increased ocean forcing in West Antarctica starting at least in the early 1970s.

In the current study, cryoconite samples were collected from six glaciers on the Tibetan Plateau to analyze n-alkanes and n-alkenes. The findings revealed that the concentrations of n-alkanes and n-alkenes varied from 40.1 to 496.1 μg g−1 and from 4.6 to 13.8 μg g−1, respectively. The carbon preference index of the long-chain n-alkanes ranged from 3.3 to 8.4, and the average chain length ranged from 28.7 to 29.3. Moreover, the δ13C of the n-alkanes in cryoconites were within the range of C3 plants, demonstrating that the n-alkanes in cryoconites were only derived from vascular plants. However, the δDmean were more negative than that of C3 plants, which could be caused by dry and humid conditions of glaciers. Unlike n-alkanes, n-alkenes ranged from C17:1 to C30:1 and showed a weak even-over-odd carbon number preference in the Dongkemadi, Yuzhufeng, Laohugou and Tianshan glacier, but a weak odd carbon preference in the Qiyi glacier. The n-alkenes in the YL Snow Mountains showed an obvious odd-over-even carbon number predominance from C17:1 to C22:1 with Cmax at C19:1, and the even-over-odd carbon number preference from C23:1 to C30:1 with Cmax at C28:1. This demonstrated that the n-alkenes of cryoconites may be mainly derived from in situ production in glaciers.

This study reports and discusses the first case of glaciotectonic landforms in the Shyok valley of the Trans-Himalayan Karakoram Range, Ladakh, where a large decomposed granite megablock (8.2 km2) along with underlying diamicton is thrust over the unconsolidated Quaternary glaciofluvial sediments along a fault gouge zone near the village of Khalsar. The absence of deformation signatures below the fault gouge indicates that the brittle fault acted as a décollement surface under frozen conditions along which the glaciotectonic megablock was translated. The other deformation features include slickensides, ductile shear, thrust propagation fold noses, clastic dykes and rafts of granite and slate within the diamicton sediments. These features indicate a subglacial glaciotectonic nappe origin of the landform. The presence of juxtaposed brittle to ductile deformation fabric, clastic dykes and the superimposition of deformed decomposed granite and diamicton over the undisturbed fluvial sediments indicates a permafrost glacial margin and proglacial environment under sufficient subglacial hydrodynamic conditions for the entrapment and transportation of the glaciotectonic megablock. The deformation fabric consistently shows a southeast orientation, indicating an advancing glacier motion from northwest to southeast. The Siachen Glacier which formerly flowed down the Nubra valley is the most likely cause of the Khalsar glaciotectonic landform.

Ice tongues are unconfined by land on their lateral margins and are sensitive to external forcing from the ocean. They are found sporadically around the Antarctic coast but are common in the western Ross Sea. Lateral flexure creates bending stresses within these ice tongues which is likely to contribute to their fragility and may restrict their spatial distribution. A three-year time series (2017–2020) of synthetic aperture radar interferometry (InSAR) and differential interferometry (DInSAR) is used to observe the lateral flexure of the 10 km long Erebus Ice Tongue as a result of ocean currents. The fast ice area around the ice tongue was mapped during the same period. It was found that when fast ice was absent (34.7% of the time), the average lateral movement of the ice tongue was twice as high (0.44 m) as when it was embedded in fast ice (0.19 m). A significant correlation (0.45) between flexure and tidal currents was found when fast ice was absent. An analytical model tuned to observations suggests that even without sea ice for stabilisation, the lateral bending stresses induced by the ocean are insufficient to cause calving without additional amplifying factors.

One fifth of Earth's volcanoes are covered by snow or ice and many have active geothermal systems that interact with the overlying ice. These glaciovolcanic interactions can melt voids into glaciers, and are subject to controls exerted by ice dynamics and geothermal heat output. Glaciovolcanic voids have been observed to form prior to volcanic eruptions, which raised concerns when such features were discovered within Job Glacier on Qw̓elqw̓elústen (Mount Meager Volcanic Complex), British Columbia, Canada. In this study we model the formation, evolution, and steady-state morphology of glaciovolcanic voids using analytical and numerical models. Analytical steady-state void geometries show cave height limited to one quarter of the ice thickness, while numerical model results suggest the void height h scales with ice thickness H and geothermal heat flux $\dot {Q}$ as $h/H = a H^b \dot {Q}^c$, with exponents b = −n/2 and c = 1/2 where n is the creep exponent. Applying this scaling to the glaciovolcanic voids within Job Glacier suggests the potential for total geothermal heat flux in excess of 10 MW. Our results show that relative changes in ice thickness are more influential in glaciovolcanic void formation and evolution than relative changes in geothermal heat flux.

The number of glacial lakes has grown globally concurrently with the retreat of glaciers in the last few decades, increasing the risk of potentially hazardous glacial lake outburst floods (GLOFs) and posing a threat to downstream communities. Norway has several known ice-dammed lakes that produce repeated GLOFs but as we show here, the existing GLOF database is incomplete and needs to be improved through continuous monitoring of glaciers and glacial lakes. This study examines the case of an ice cap in central Norway hosting at least four drainage-prone lakes. We reconstruct the sequence of lake drainage events through a combination of remote sensing, ground-truthing and citizen science while scrutinising the applicability of the PlanetScope imagery vs Sentinel-2 and Landsat-8 OLI products. As opposed to the Landsat imagery that often fails to resolve even the largest glacial lakes of Tystigbreen, both PlanetScope and Sentinel-2 are helpful in identifying previously unrecognised glacial lakes and undocumented drainage events. Our analysis suggests that a fusion of the two satellite products may be beneficial for automated tracking of glacial lake changes. We also demonstrate that local knowledge and systematic involvement of citizens in data collection have a potential to enrich GLOF databases.

New experiments have revealed that the addition of a small amount of wood ash doubles the angle of repose of granulated ice at −10°C, but has little effect at −30°C. It is hypothesized that this behavior may be explained in terms of the freezing of water that is formed at the higher temperature through the reaction of salt within the ash.

Recently released Moderate-Resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) collection 6.1 (C6.1) products are useful for understanding ice–atmosphere interactions over East Antarctica, but their accuracy should be known prior to application. This study assessed Level 2 and Level 3 MODIS C6.1 LST products (MxD11_L2 and MxD11C1) in comparison with the radiance-derived in situ LSTs from 12 weather stations. Significant cloud-related issues were identified in both LST products. By utilizing a stricter filter based on automatic weather station cloud data, despite losing 29.4% of the data, accuracy of MODIS LST was greatly improved. The cloud-screened MODIS LST exhibited cold biases (−5.18 to −0.07°C, and root mean square errors from 2.37 to 6.28°C) than in situ LSTs at most stations, with smaller cold biases at inland stations, but larger ones at coastal regions and the edge of plateau. The accuracy was notably higher during warm periods (October–March) than during cold periods (April–September). The cloud-screened MODIS C6.1 LST did not show significant improvements over C5 (Collection 5) version across East Antarctica. Ice-crystal precipitation occurring during temperature inversions at the surface (Tair-Tsurface) played a crucial role in MODIS LST accuracy on inland plateau. In coastal regions, larger MODIS LST biases were observed when the original measurements were lower.