Search

Glaciers and ice sheets are an integral part of Earth's system, advancing and retreating in response to changes in climate. Clues about the past, present and future behaviour of these ice masses are found throughout current and former glaciated landscapes. In this commentary, we outline recent scientific advances from a collection of articles in which geomorphological evidence is used to inform us about the behaviour of glaciers and ice sheets across a range of spatial (landform to continent) and temporal (seasons to millennia) scales. Through a diversity of approaches including field measurements, remote sensing and numerical modelling, these studies build on an extensive background literature to deepen our understanding of how ice flows, how glaciers and ice sheets respond to climate change, and of the processes of ice advance and retreat and the stability of the system. Further integration of knowledge across the fields of geomorphology and glaciology will have tangible benefits for managing the societal and environmental impacts of glacier change and for improved projections of sea‐level rise over the coming decades to centuries. [ABSTRACT FROM AUTHOR]

This review of the dynamics of the flow of glaciers and ice sheets focusses on the mathematical models which have been developed to explain a number of observations in the behaviour of these large ice masses. It takes a personal view, starting with an account of the historical development of the subject, followed by a tour of some of the observational puzzles and theoretical solutions which I have found most intriguing. Most of the discourse involves basal processes: glacier sliding, subglacial hydrology, bedform generation. I describe in some detail wave ogives, surging glaciers, subglacial floods and subglacial bedforms. I variously explore issues in sliding theory: sub-temperate sliding, slip or till deformation; the rheology of till; the status of the instability theory of drumlin formation. And in an effort to broaden the field of view, I include discussion of Dansgaard-Oeschger and Heinrich events, grounding line stability, and the formation of eskers and, particularly, flutes. [ABSTRACT FROM AUTHOR]

Glacier flow modulates sea level and is governed largely by the viscous deformation of ice. Multiple molecular-scale mechanisms facilitate viscous deformation, but it remains unclear how each contributes to glacier-scale deformation. Here, we present a model of ice deformation that bridges laboratory and glacier scales, unifies existing estimates of the viscous parameters, and provides a framework for estimating the parameters from observations and incorporating flow laws derived from laboratory observations into glacier-flow models. Our results yield a map of the dominant deformation mechanisms in the Antarctic Ice Sheet, showing that, contrary to long-standing assumptions, dislocation creep, characterized by a value of the stress exponent n = 4, likely dominates in all fast-flowing areas. This increase from the canonical value of n = 3 dramatically alters the climate conditions under which marine ice sheets may become unstable and drive rapid rates of sea-level rise. [ABSTRACT FROM AUTHOR]

The objectives of this paper are to investigate the trade-offs between a physically constrained neural network and a deep, convolutional neural network and to design a combined ML approach ("VarioCNN"). Our solution is provided in the framework of a cyberinfrastructure that includes a newly designed ML software, GEOCLASS-image (v1.0), modern high-resolution satellite image data sets (Maxar WorldView data), and instructions/descriptions that may facilitate solving similar spatial classification problems. Combining the advantages of the physically-driven connectionist-geostatistical classification method with those of an efficient CNN, VarioCNN provides a means for rapid and efficient extraction of complex geophysical information from submeter resolution satellite imagery. A retraining loop overcomes the difficulties of creating a labeled training data set. Computational analyses and developments are centered on a specific, but generalizable, geophysical problem: The classification of crevasse types that form during the surge of a glacier system. A surge is a glacial catastrophe, an acceleration of a glacier to typically 100–200 times its normal velocity. GEOCLASS-image is applied to study the current (2016-2024) surge in the Negribreen Glacier System, Svalbard. The geophysical result is a description of the structural evolution and expansion of the surge, based on crevasse types that capture ice deformation in six simplified classes. [ABSTRACT FROM AUTHOR]

Stratovolcanoes are common globally, with high‐altitude summit regions that are often glacier‐clad and intersect the seasonal and perennial snow line. During an eruption, interaction between snow/ice and hot, pyroclastic deposits will potentially lead to extensive melt and steam production. This is particularly pertinent when pyroclastic density currents (PDCs) are emplaced onto and propagate over glacierised substrates. Generated melt and steam are incorporated into the flow, which can cause a transformation from a hot, dry granular flow, to a water‐saturated, sediment‐laden flow, termed a lahar. Both PDCs and ice‐melt lahars are highly hazardous due to their high energy during flow and long runout distances. Knowledge of the physics that underpin these interactions and the transformation to ice‐melt lahar is extremely limited, preventing accurate descriptions within hazard models. To physically constrain the thermal interactions we conduct static melting experiments, where a hot granular layer was emplaced onto an ice substrate. The rate of heat transfer through the particle layer, melt and steam generation were quantified. Experiments revealed systematic increases in melt and steam with increasing particle layer thicknesses and temperatures. We also present a one‐dimensional numerical model for heat transfer, calibrated against experimental data, capable of accurately predicting temperature and associated melting. Furthermore, similarity solutions are presented for early‐time melting which are used to benchmark our numerical scheme, and to provide rapid estimates for meltwater flux hydrographs. These data are vital for predicting melt volume and incorporation into PDCs required to facilitate the transformation to and evolution of ice‐melt lahars. Plain Language Summary: When volcanoes explosively erupt they may produce avalanches of hot, dry volcanic ash. When these volcanic avalanches occur on snow and glacier‐covered volcanoes, they produce steam and melt, that can mix with the volcanic avalanche, transforming it to a cool, wet volcanic mudflow. Both volcanic avalanches and mudflows are extremely destructive and dangerous due to their high speeds and long flow paths. Historically, these flows have resulted in many fatalities and extensive building and infrastructure damage. We investigate the conditions under which transformation from volcanic avalanches to mudflows can occur. We use small‐scale laboratory experiments to measure the transfer of heat, steam, and melt generation when a hot ash layer is emplaced onto an ice layer. A numerical model is presented to describe this heat transfer at large scales, like in natural volcanic settings. This can be used to estimate the amount of melt required to cause this transformation from volcanic avalanche to mudflow. This can help us predict the destructiveness of these interactive events, and help us convey the hazard to stakeholders, and populations living in regions affected by volcano‐ice interactions. Key Points: Pyroclast‐ice interactions are investigated using laboratory experiments and numerical modelingThe parameter space that governs melting and steam generation is estimatedMeltwater source flux hydrographs show similarities with rainfall‐driven lahar source hydrographs [ABSTRACT FROM AUTHOR]

Wind‐driven redistribution of surface snow is one of the key factors leading to heterogeneous accumulation of snow at small scales. Understanding the processes that lead to this heterogeneous accumulation is, therefore, of great importance to many glaciological and hydrological questions. High‐quality information on the wind field is necessary to realistically represent drifting snow. Here, we introduce a novel, intermediate‐complexity drifting and blowing snow module for the Weather Research and Forecasting (WRF) model that integrates seamlessly into the standard WRF infrastructure. The module also accounts for snow particle sublimation and considers the thermodynamic feedback on the atmospheric fields. In an idealized model environment the module was tested for physical consistency. Sensitivity experiments showed that drifting snow sublimation has on the one hand local effects on the erosion and deposition patterns and on the other hand non‐local effects on the larger‐scale atmosphere. The simple and computationally efficient implementation allows this module to be used for small‐scale and large‐scale simulations in polar and glaciated regions. Plain Language Summary: Transport of snow by the wind can have high impact on local glacier mass changes as it leads to non‐uniform amounts of snow on the ground. In order to simulate and better understand this process we introduce a new modeling framework that is included into a widely used atmospheric model. Test simulations and sensitivity experiments show the physical consistency of the model. Complex interactions between different processes like snow erosion, drifting snow sublimation and the wind field show the necessity of coupling the snow and atmospheric models. Key Points: We introduce a new module to simulate drifting snow within the framework of the Weather Research and Forecasting (WRF) modelThe module is tested in an idealized setting for physical consistencyWe show that local sublimational cooling has larger‐scale non‐local effects on the wind field [ABSTRACT FROM AUTHOR]

Herein, fast fracture initiation in glacier ice is modeled using a Material Point Method and a simplified constitutive law describing tensile strain softening. Relying on a simple configuration where ice flows over a vertical step, crevasse patterns emerge and are consistent with previous observations reported in the literature. The model's few parameters allows identification of a single dimensionless number controlling fracture spacing and depth. This scaling law delineates two regimes. In the first one, ice thickness does not play a role and only ice tensile strength controls the spacing, giving rise to numerous surface crevasses, as observed in crevasse fields. In this regime, scaling can recover classical values for ice tensile strength from macroscopic field observations. The second regime, governed by ice bending, produces large‐scale, deep fractures resembling serac falls or calving events. Plain Language Summary: In ice sheets and alpine glaciers, fast‐flowing sections are often characterized by crevasse fields that play a significant role in the cryo‐hydrologic system by facilitating meltwater flow, enhancing basal sliding, weakening the ice, and impacting glacier thermodynamics. Modeling these fractures at the glacier scale remains challenging and often necessitates integrating diverse models which hinders the straightforward consideration of physical issues associated with crevasse fields on a large scale. Here, a new numerical framework allows us to conduct field‐scale experiments and paves the way for a scaling law to elucidate the macroscopic factors influencing fracture fields and to easily incorporate crevasse depth and spacing into large‐scale models. A newly discovered scaling law highlights the transition between a mechanical behavior where the regular crevasse spacing is unaffected by geometry to a regime where geometry plays a significant role, particularly in large‐scale fracture processes like glacier calving. While the numerical experiments in this paper focus on glaciers, the model and conceptual framework is versatile and can address the mechanical behavior of fractures in broader geophysical contexts such as snow, rock or ice avalanches, tectonics and landslides. Key Points: Fractures in glacier flow are modeled using material point method with elastoplasticity and tensile strain softeningA dimensional analysis reveals a key dimensionless number characterizing two different regimes of fast fractureOne regime predicts acknowledged ice tensile strength from field observations and characterizes the regular crevasse spacing [ABSTRACT FROM AUTHOR]

Reconstructing Pleistocene glaciation timing and extent is vital for understanding paleoclimate. Whereas late Pleistocene glaciation has been studied extensively in western North American mountain ranges, the glacial history of the western Elk Range in Colorado remains understudied, particularly in the East River watershed, a site of intense scientific focus. Here we use cosmogenic nuclide exposure and depth–profile dating methods to determine the timing of glaciation in the East River watershed. We use glacier modeling to reconstruct paleoglacier extents and quantify past climate conditions. Our findings indicate that the East River glacier retreated from its maximum position approximately 17–18 ka, moving to recessional positions between 13 and 15 ka, before experiencing more substantial retreat to high-elevation cirques around 13 ka. Glacier modeling suggests that the maximum ice extents at 17–18 ka could have been sustained by temperature depressions of approximately −6.5°C compared to modern conditions, assuming consistent precipitation. Additionally, the ice position at 13–15 ka could have been supported by temperature depressions of around −4.0°C. These results offer insights into the deglaciation timeline in the East River watershed and broader western Elk Range as well as paleoclimate conditions during the late Pleistocene, which may aid future research on critical zone evolution in the East River watershed. [ABSTRACT FROM AUTHOR]

Today, relatively warm Circumpolar Deep Water is melting Thwaites Glacier at the base of its ice shelf and at the grounding zone, contributing to significant ice retreat. Accelerating ice loss has been observed since the 1970s; however, it is unclear when this phase of significant melting initiated. We analyzed the marine sedimentary record to reconstruct Thwaites Glacier's history from the early Holocene to present. Marine geophysical surveys were carried out along the floating ice-shelf margin to identify core locations from various geomorphic settings. We use sedimentological data and physical properties to define sedimentary facies at seven core sites. Glaciomarine sediment deposits reveal that the grounded ice in the Amundsen Sea Embayment had already retreated to within ~45 km of the modern grounding zone prior to ca. 9,400 y ago. Sediments deposited within the past 100+ y record abrupt changes in environmental conditions. On seafloor highs, these shifts document ice-shelf thinning initiating at least as early as the 1940s. Sediments recovered from deep basins reflect a transition from ice proximal to slightly more distal conditions, suggesting ongoing grounding-zone retreat since the 1950s. The timing of ice-shelf unpinning from the seafloor for Thwaites Glacier coincides with similar records from neighboring Pine Island Glacier. Our work provides robust new evidence that glacier retreat in the Amundsen Sea was initiated in the mid-twentieth century, likely associated with climate variability. [ABSTRACT FROM AUTHOR]

Copyright of Social Anthropology / Anthropologie Sociale is the property of Berghahn Books and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

As the main contributor to global sea-level rise, the Greenland Ice Sheet (GrIS) has undergone significant mass change over the last two decades. The satellite mission of GRACE (Gravity Recovery And Climate Experiment) and its follow-on mission (GRACE-FO) provide accurate observations but low-spatial resolution. In contrast, satellite altimetry provides observations at a high-spatial resolution but with large uncertainties, limiting the understanding of glacier-scale mass change. To derive accurate and high-spatial resolution mass change estimates from GRACE/GRACE-FO observations, we present a novel constraint mascon method in which the regularization matrix is constructed with the signal variances from satellite altimetry. Based on the proposed method, we derive a series of high-resolution (25 km × 25 km) monthly mascon solutions from 2002 April to September. The glacier-scale estimates from the input–output method agree better with those from our mascon solutions than those from the global mascons of CSR (Center for Space Research, the University of Texas), JPL (Jet Propulsion Laboratory), and GSFC (Goddard Space Flight Center), with a higher linear regression coefficient of 0.71. Benefitting from the greatly improved spatial resolution, our estimates provide the first accurate monthly glacier-scale mass change estimates from GRACE/GRACE-FO observations over the GrIS, to our knowledge. The results show that 20 of the 260 glaciers contributed to more than 42 per cent of the ice loss in the GrIS from 2002 to 2022. Most strikingly, the mass loss of Jakobshavn Isbrae was the most significant at –18.7 ± 0.05 Gt yr−1, accounting for 7.4 per cent of the total in the GrIS during the study period. Furthermore, we find that the SMB (surface mass balance) and ice-dynamics-related mass changes contribute nearly equally to the observed mass changes, but the corresponding spatiotemporal characteristics differ. SMB contributed the most to the mass change in 2019, while ice dynamics played the most dominant role in 2018. Moreover, the SMB effect is widespread, and ice-dynamics-related mass loss is mainly concentrated in five small areas. [ABSTRACT FROM AUTHOR]

The increased melting and rapid retreat of marine-terminating glaciers is a key contributor to sea-level rise. In glacial fjords with shallow sills common in Patagonia, Alaska, and other systems, these bathymetric features can act as a first-order control on the dynamics. However, our understanding of how this shallow bathymetry interacts with the subglacial discharge from the glacier and impacts the fjord circulation, water properties, and rates of submarine melting is limited. To address this gap, we conduct idealized numerical simulations using a coupled plume–ocean fjord model spanning a wide range of initial ocean conditions, sill depths, and subglacial discharge. A previously documented circulation regime leads to strong mixing and vertical transport over the sill, where up to ∼ 70 % of the colder water from the upper-layer outflow is refluxed into the deeper layer, cooling the incoming warm oceanic water by as much as 1 ∘C and reducing the stratification near the glacier front. When the initial stratification is relatively strong or the subglacial discharge is relatively weak, an additional unsteady circulation regime arises where the freshwater flow can become trapped below the sill depth for weeks to months, creating an effective cooling mechanism for the deep water. We also find that submarine melting often increases when a shallow sill is added to a glacial fjord due to the reduction of stratification – which increases submarine melting – dominating over the cooling effect as the oceanic inflow is modified by the presence of the sill. These results underscore that shallow-silled fjords can have distinct dynamics that strongly modulate oceanic properties and the melting rates of marine-terminating glaciers. [ABSTRACT FROM AUTHOR]

Climate change is driving the loss of alpine glaciers globally, yet investigations about the water quality of rivers stemming from them are few. Here we provide an overview assessment of a biogeochemical data set containing 200+ parameters that we collected between 2019 and 2021 from the headwaters of three such rivers (Sunwapta‐Athabasca, North Saskatchewan, and Bow) which originate from the glacierized eastern slopes of the Canadian Rocky Mountains. We used regional hydrometric data sets to accurately model discharge at our 14 sampling sites. We created a Local Meteoric Water Line (LMWL) using riverine water isotope signatures and compared it to collected regional rain, snow, and glacial ice signatures. Principal component analyses of river physicochemical measures revealed distance from glacier explained more data variability than other spatiotemporal factors (i.e., season, year, or river). Discharge, chemical concentrations, and watershed areas were then used to model site‐specific open water season yields for 25 parameters. Chemical yields followed what would generally be expected along river continuums from glacierized to montane altitudinal life zones, with landscape characteristics driving chemical sources and sinks. For instance, particulate chemical yields were generally highest near source glaciers with proglacial lakes acting as settling ponds, whereas most dissolved yields varied by parameter and site. As these headwaters continue to evolve with glacier mass loss, the data set and analyses presented here can be used as a contemporary baseline to mark future change against. Further, following this initial assessment of our data set, we encourage others to mine it for additional biogeochemical studies. Plain Language Summary: Alpine glaciers are vulnerable to climate change, with their numbers and sizes expected to decline dramatically within the current century. It is known that a decrease in glacier mass will have direct consequences on glacial meltwater quantity, but how this may impact the quality of receiving freshwaters is less understood. Our goal was to record a water quality data set of 200+ physical and chemical parameters for the headwaters of three major rivers draining glaciers on the eastern slopes of the Canadian Rocky Mountains. We first used statistics to model how fast the rivers were flowing at our 14 sampling locations over our 2‐year sampling period. We then looked at how and why water isotopes collected from river, rain, snow, and glacier ice samples differed. We discovered that distance downstream from the glacier explained more statistical variation in our physical and chemical measures than season, year, or river. Finally, at each of our sampling sites, we calculated chemical yields, or the amount of chemical constituent being exported downriver per unit watershed area, and found landscape features like lakes, outwash plains, and forests influence them. We ultimately hope that others use our comprehensive data set to investigate this changing region further. Key Points: We assessed a 200+ parameter data set collected from glacial river headwaters on the eastern slopes of the Canadian Rocky MountainsOur findings can act as a contemporary reference for future investigation of glacial headwaters in rapidly evolving alpine regionsThe possibilities for data exploration of our 260,000+ measure data set are numerous and we encourage others to mine it for new studies [ABSTRACT FROM AUTHOR]

The Bering‐Bagley Glacier System (BBGS), Alaska, Earth's largest temperate surging glacier, surged in 2008–2013. We use numerical modeling and satellite observations to investigate how surging in a large and complex glacier system differs from surging in smaller glaciers for which our current understanding of the surge phenomenon is based. With numerical simulations of a long quiescent phase and a short surge phase in the BBGS, we show that surging is more spatiotemporally complex in larger glaciers with multiple reservoir areas forming during quiescence which interact in a cascading manner when ice accelerates during the surge phase. For each phase, we analyze the simulated elevation‐change and ice‐velocity pattern, infer information on the evolving basal drainage system through hydropotential analysis, and supplement these findings with observational data such as CryoSat‐2 digital elevation maps. During the quiescent simulation, water drainage paths become increasingly lateral and hydropotential wells form indicating an expanding storage capacity of subglacial water. These results are attributed to local bedrock topography characterized by large subglacial ridges that dam the down‐glacier flow of ice and water. In the surge simulation, we model surge evolution through Bering Glacier's trunk by imposing a basal friction representation that mimics a propagating surge wave. As the surge progresses, drainage efficiency further degrades in the active surging‐zone from its already inefficient, end‐of‐quiescence state. Results from this study improve our knowledge of surging in large and complex systems which generalizes to glacial accelerations observed in outlet glaciers of Greenland, thus reducing uncertainty in modeling sea‐level rise. Plain Language Summary: The Bering‐Bagley Glacier System (BBGS), Alaska, Earth's largest temperate surging glacier, recently surged in 2008–2013. A surge glacier cycles between a long period of normal flow and a short period of accelerated flow where large‐scale deformations, such as crevasses, occur. This paper focuses on investigating a surge in a large and complex system rather than a small glacier where most studies on surges have been conducted. We use a numerical model to simulate glacier evolution for both the quiescent phase and the initial surge phase of the BBGS. For each phase, we analyze the simulated elevation‐change and ice‐velocity, and infer information on the evolving hydrologic drainage system. During the quiescent phase, ice‐mass builds up at locations consistent with those observed and water drainage paths become longer with expanding capacity to store subglacial water. These results are attributed to local bedrock topography characterized by large subglacial ridges that act to dam the down‐glacier flow of ice and water. In the surge simulation, we model surge evolution through Bering Glacier by implementing a new friction representation that mimics a propagating wave. As the surge progresses through the glacier, drainage efficiency further degrades in the areas of fast‐moving ice. Key Points: Using a full‐Stokes approach informed by satellite observations, a quiescent and surge phase of the Bering‐Bagley Glacier System (BBGS) are simulatedLocal bed topography controls the formation of several reservoir areas, which lengthen down‐glacier drainage paths during quiescenceA friction representation for the surge phase is implemented based on observed properties of kinematic surge waves in the BBGS [ABSTRACT FROM AUTHOR]

In this study, we develop a model that allows to combine gravimetric and geometric data. By the combination, we improve the spatial resolution of the resulting mass balance estimate compared to a purely gravimetric one. The equivalent ice or firn density of the changing ice volume is estimated within a mathematical inversion model, which includes geometric information about the volumetric change of the ice sheet and the resulting gravity change. This gravity change is computed from monthly GRACE gravity fields. They have a limited spatial resolution of a few 100 km, but allow direct conclusions about the true mass changes over Greenland. The ice-volume changes are described by a product by the Climate Change Initiative of European Space Agency, which is based on altimetry data. They have a very fine spatial resolution (down to a few km), but are not directly sensitive to mass changes. By combining both data sets in a common mathematical model, the advantages of both data types (direct sensitivity to mass versus high spatial resolution) are made use of. In this way, we improve the spatial resolution of mass balance estimates over Greenland. This leads to a map of mass trends, which has the same spatial resolution as the input map of geometric changes, but which is consistent with the input gravity fields. It will enable improving the localization of mass change signals of ice sheets and glaciers, which are usually rather small-scale. We compare our estimates to the results of complementary studies regarding the total mass loss of the Greenland ice sheet and its surrounding land surface. Our study leads to a value of |$-213\pm 37\, \text{Gt}\,\text{a}^{-1}$| in the time span from 2011 to 2015. We also discuss the problem of separating the mass contribution of the Greenland ice sheet itself and its surrounding region. [ABSTRACT FROM AUTHOR]

In addition to having far‐reaching impacts on human health, agriculture, wildfires, ecosystems, and infrastructure, heatwaves control streamflow through the melting of seasonal snow and glacier ice. Despite their importance, there is limited understanding of how heatwaves modify streamflow at regional scales, how these impacts vary by heatwave timing and duration, and how glaciers control the streamflow response. Here, we use a deep learning hydrological model, which has previously been trained, evaluated, and interpreted in southwestern Canada, to simulate the streamflow response to heatwaves at 111 basins in the region. The model, driven by gridded ERA5 reanalysis temperature and precipitation data from 1979 to 2015, is forced by synthetic heatwave conditions that vary in their duration and onset throughout the year. We consider how the streamflow response to heatwaves is sensitive to annual temperatures by adding spatially and temporally uniform warming of 2°C across the study region, under the assumption that the underlying hydrological system behavior remains unchanged. We find that heatwaves, particularly in spring and summer, induce an initial streamflow surplus followed by a streamflow deficit, relative to the non‐heatwave case. In summer, glacier contributions to streamflow partially compensate for streamflow deficits that arise from heatwaves earlier in the melt season. In the scenario with 2°C warmer annual temperatures, heatwaves induce a lesser streamflow response in spring when the seasonal streamflow is most increased due to the advancing freshet. Our findings demonstrate how glaciers buffer the impacts of heatwaves on streamflow, but this buffering effect is expected to diminish as glaciers retreat. Plain Language Summary: Heatwaves are periods of extremely warm temperatures that can have profound consequences on not only human health, agriculture, wildfires, ecosystems, and infrastructure, but also on river flows. We use a model to simulate river flows throughout southwestern Canada under heatwave conditions. We find that 1 week of extreme heat can influence river flows for several months by changing the timing and amount of snow and glacier melt. Excess melt during heatwaves can lead to river flows that initially are higher than they otherwise would be; however, river flows later in the season are then lower than they otherwise would be. Glacierized rivers are unique in that glaciers melt during summer, which can compensate for heatwaves that melt snow earlier in spring. If average temperatures were warmer by 2°C, heatwaves could melt more snow and ice during winter and early spring, but less in late spring and summer. Key Points: We use a deep learning hydrological model to quantify streamflow responses to simulated heatwaves that vary in duration and date of onsetHeatwaves drive initial streamflow surpluses and delayed streamflow deficits that are linked to heatwave timing and glacier coverageWith warmer annual temperatures, spring heatwaves drive lesser streamflow responses when streamflow is most increased by the earlier freshet [ABSTRACT FROM AUTHOR]

Small glaciers are one of the best indicators of climatic variations and their short‐term effects. Located in the Spanish Pyrenees, the Maladeta is one of these glaciers. Its systematic observation began in the 1980s, being one of the few Pyrenean glaciers with a tongue‐shaped front. This study presents the evolution of the Maladeta glacial tongue over a decade (2010–2020) through multiple geomatic techniques. Surveys have ranged from Total Stations and Global Navigation Satellite System devices to massive data capture techniques such as Terrestrial Laser Scanners or Unmanned Aerial Vehicles photogrammetry. The aim is to analyze in detail the loss of surface area and thickness of the glacier and its transition from being a glacier with a tongue partially determined by climate to a topoclimatically determined cirque glacier. The results reveal a tongue retreat of over 5 m/yr and area losses of over 0.2 ha/yr, along with ice thickness and volume losses of −1.7 m/yr and over −21 × 103 m3/yr, respectively. If this trend continues, the tongue, and possibly the Maladeta glacier, could disappear by the end of the 2030s. [ABSTRACT FROM AUTHOR]

A large ensemble of ice sheet projections to the end of the 21st century have been compiled within community-based initiatives. These ensembles allow for assessment of uncertainties in projections associated with climate forcing and a wide range of parameters governing ice sheet and shelf dynamics, including ice-ocean interactions. Herein, we compute geographically variable sea level 'fingerprints' associated with ∼320 simulations of polar ice sheet projections included in the Ice Sheet Model Intercomparison Project for CMIP6 and ∼180 projections of glacier mass changes from the Glacier Model Intercomparison Project. We find a strong correlation (coefficient > 0.97) between all fingerprints of Greenland Ice Sheet projections when considering a global region outside the near field of the ice sheet. Consistency in the fingerprints for the Antarctic Ice Sheet (AIS) projections is much weaker, though correlation coefficients > 0.80 were found for all projections with global mean sea level (GMSL) greater than 10 cm. The far-field variability in the fingerprints associated with the AIS is due in large part to the sea level change driven by Earth rotation changes. The size and position of the AIS on the south pole makes the rotational component of the sea level fingerprint highly sensitive to the geometry of the ice mass flux, a geometry that becomes more consistent as the GMSL associated with the ice sheet projection increases. Finally, the fingerprints of glacier mass flux show an intermediate level of consistency, with contributions from Antarctic glaciers being the primary driver of decorrelation. [ABSTRACT FROM AUTHOR]

In Greenland, 87% of the glacierized area terminates in the ocean, but mass lost at the ice‐ocean interface, or frontal ablation, has not yet been fully quantified. Using measurements and models we calculate frontal ablation of Greenland's 213 outlet and 537 peripheral glaciers and find a total frontal ablation of 481.8 ± 24.0 for 2000–2010 and 510.2 ± 18.6 Gt a−1 for 2010–2020. Ice discharge accounted for ∼90% of frontal ablation during both periods, while mass loss due to terminus retreat comprised the remainder. Only 16 glaciers were responsible for the majority (>50%) of frontal ablation from 2010 to 2020. These estimates, along with the climatic‐basal balance, allow for a more complete accounting of Greenland Ice Sheet and peripheral glacier mass balance. In total, Greenland accounted for ∼90% of Northern Hemisphere frontal ablation for 2000–2010 and 2010–2020. Plain Language Summary: We estimate the mass of ice lost from all Greenland glaciers that entered the ocean during each of the last two decades. This ice loss at the front of these marine‐terminating glaciers is called frontal ablation and is approximately equal to the mass of icebergs entering the ocean. Frontal ablation is important because 87% of glacier area in Greenland ends in the ocean, through 750 outlets, and previous work has only approximated frontal ablation. This study quantifies it for the first time, helping to close the mass budget for the Greenland Ice Sheet and better partition its mass balance into components. We find that Greenland accounts for ∼90% of all Northern Hemisphere frontal ablation and, of that contribution, just 17 glaciers for 2000–2010 and 16 glaciers for 2010–2020 account for more than half of total Greenland frontal ablation. Key Points: Frontal ablation of the Greenland Ice Sheet averaged 510.2 ± 18.6 Gt a−1 for 2010–2020, ∼90% of which came from ice dischargeThe frontal ablation we measured is larger than the total mass loss from the ice sheet, indicating a positive climatic‐basal balanceOnly 16 glaciers account for 50% of the total frontal ablation from the Greenland Ice Sheet [ABSTRACT FROM AUTHOR]

What is the relationship between vision and sound in more‐than‐human environmental sensing? This article traces an ethnography of glaciologists' experiences with technological sensing systems that surpass human sensing capabilities, producing an expansion of sensory knowledge that enmeshes both imagery and acoustics. Sound and vision emerge no longer as separate modalities, but in a united vocabulary of sensing in which the human and the machine are collaborators, producing a multi‐eared and multi‐eyed system by which icebergs and glaciers are observed and perceived as a quickly morphing process, rather than as a static object. It examines the manifold aesthetic variations and conversions of data from acoustic sensing systems into sounds and images that work together to reshape the scientific imagination of the cryosphere. This process ultimately reveals how technology changes perceptions and how glaciers, human bodies, and machines become intertwined with each other in a more‐than‐human system that holds the promise to diversify knowledge. [ABSTRACT FROM AUTHOR]

Bed topography is a critical parameter for determining the modern-day and future dynamics of ice sheets and their outlet glaciers. This is because the topography controls the state of stress for glaciers. At glacier termini, topography can influence the timing of terminus retreat by controlling access to warm ocean waters and/or by influencing the ability of a glacier terminus to retreat over bed bumps (moraines). Inland from the terminus, the topography can also influence where glacier retreat and thinning can stabilize. In part, this is because of knickpoints in bed topography created through glacial erosion that may influence the extent to which thinning can diffuse inland for an individual glacier and thus, the timing and magnitude of long-term mass loss. Here we provide a review of the current literature on these topics. While much of the reviewed literature assumes that topography is stable on relevant timescales to humans, new research suggests that topography may change much faster than previously thought and this further complicates our ability to project future outlet glacier change. [ABSTRACT FROM AUTHOR]

Ice cliffs are features along ice sheet margins, along tropical mountain glaciers, at termini of mountain glaciers and on debris-covered glacier tongues, that have received scattered attention in literature. They cover small relative areas of glacier or margin surface respectively, but have been involved in two apparent anomalies. On the one hand, they have been identified as potential hotspots of extreme melt rates on debris-covered tongues contributing to their relatively rapid ablation, compared to the surrounding glacier surface. On the other hand, they appear where the ice margin is stable (or temporarily advancing) even under conditions of negative mass balance. In this manuscript, we recapitulate why ice cliffs remain interesting features to investigate and what we know about them so far. We conclude by suggesting to further investigate their genesis and variable morphology and their potential as windows into past climates and processes. [ABSTRACT FROM AUTHOR]

Greenland's glacial fjords modulate the exchange between the ice sheet and ocean. Subglacial‐discharge‐driven plumes adjacent to glaciers may exert an important influence on fjord water properties, submarine glacier melting and the export of glacially‐modified waters to the shelf. Here we use a numerical plume model in conjunction with observations from proximal to 14 glaciers in northwest Greenland to assess the impact of these plumes on near‐glacier water properties. We find that in late summer, waters emanating from glacial plumes often make up >50% of the fjord water composition at intermediate depths. These plume waters are comprised largely of upwelled Atlantic Water, warming the near‐glacier water profile and likely increasing submarine melting. Our findings demonstrate the key role played by plumes in driving water modification in Greenland's fjords, and the potential for simple models to capture these impacts across a range of settings. Plain Language Summary: Fjords form an important point of contact between the Greenland Ice Sheet and surrounding ocean. Warm ocean waters in these fjords melt the submerged sections of marine‐terminating glaciers, whilst freshwater draining from the ice sheet modifies coastal water properties, impacting ocean circulation. In this study, we use observations of fjord temperature and salinity in combination with numerical modeling to evaluate the mechanisms through which glaciers modify fjord circulation and water properties. Across 14 glaciers in northwest Greenland, we find that the input of glacial meltwater drives the upwelling of warm, deep Atlantic‐origin waters, which mix with the cooler Arctic‐origin waters typically found above ∼250 m depth. These glacially modified waters thus cause a substantial warming of fjord waters relative to those off the coast, a feedback that will amplify the oceanic melting of these glaciers. Key Points: We assess the impact of subglacial‐discharge plumes on water properties proximal to 14 tidewater glaciers in northwest GreenlandIn late summer, waters emanating from glacial plumes often make up >50% of the fjord water composition at intermediate depthsThese plume waters are comprised largely of upwelled Atlantic Water, warming the near‐glacier water profile and likely increasing submarine melting [ABSTRACT FROM AUTHOR]

Icequakes associated with crevassing provide key information on stress fields and water flow pathways inside glaciers, conditions which are sensitive to temperature and precipitation. In this study, we analyse seismic data from 321 days recorded by five broad-band seismometers on Urumqi Glacier No. 1, Tianshan, China, to investigate the source characteristics of icequakes. We first apply a multidimensional autoregressive maximum-likelihood algorithm to detect 12 template events from 62 days of continuous waveform data. The centroid locations of the 12 template events are located using a grid-search approach. We then obtain a catalogue of 41   895 icequakes by compiling all seismic data through a template-matching method. The local magnitudes of all icequakes range between −3.3 and 0.5, with a completeness magnitude of −2.95 and Gutenberg–Richter b -value of 1.08. The particle motion and frequency contents of the 12 template events show that the dominant phases on the vertical component are Rayleigh waves at 5–40 Hz, suggesting that these icequakes were generated by near-surface crevasses. The average daily icequake number for each month and diurnal icequake variations show clear seasonal features. The icequakes may be related to growth of the surface crevasses due to additional melting water and precipitation in summer and to the thermal contraction of the ice body in winter. [ABSTRACT FROM AUTHOR]

Greenland glaciers have three primary seasonal ice flow patterns, or "types": terminus‐driven, runoff‐driven, and runoff‐adapting. To date, glacier types have been identified by analyzing flow at a single location near the terminus; information at all other locations is discarded. Here, we use principal component (PC)/empirical orthogonal function (EOF) analysis to decompose multi‐year time series of glacier speed, combined from three satellite‐derived products at four glaciers feeding Sermilik Fjord, Greenland. This improves on single‐point methods by yielding spatial patterns (EOFs), which ensure the result reflects data at all locations on the glacier, and associated temporal patterns (PCs), which allow identification of glacier type. We find that the leading EOF is uniformly signed over the entire glacier domain, that this mode explains the majority of the variance in speed, and therefore that glacier type can be inferred from the leading PC. We find that Helheim Glacier was terminus‐driven, Fenris Glacier and Midgard Glacier were runoff‐adapting, and Pourquoi Pas Glacier was runoff‐driven over 2016–2021. Our classification agrees with previous work for Helheim and Midgard Glaciers, but differs at the other two. At all but Fenris Glacier, the leading PC correlates significantly with the speed pattern observed at the single point used in previous analyses. Thus, Fenris Glacier has more complex flow patterns than single‐point analysis can capture, and wider spatial analysis techniques such as EOF/PC are required. We suggest that, due to its low computational cost and inclusion in standard analysis packages, EOF/PC analysis should be used for assessing glacier type. Plain Language Summary: Glaciers change their flow speed throughout the year. Most glaciers move slowest in winter and more quickly during summer, but subtle differences in the timing give us clues as to what controls the speed of the glacier. Some glaciers respond to ocean conditions ("terminus‐driven"), others respond to ice melting at the top surface of the glacier ("runoff‐driven"), and still others adjust their water systems as the ice melts, responding to the same melt in a different way ("runoff‐adapting"). We do not know the types of all glaciers, nor what causes a glacier to be a certain type. Previous work identified the types of a few dozen glaciers around Greenland by manually examining ice flow at a single point on each glacier. Here, we take a wider approach by mathematically extracting patterns in ice flow across the entire glacier, not just a single point. We analyze four glaciers that flow into the same fjord in East Greenland. Our results broadly agree with the previous simpler analyses, but differ at two glaciers and identify possible "multi‐type" glaciers. Thus, our method holds promise in the quest to discover what controls the seasonal flow patterns of Greenland glaciers. Key Points: Glacier type varies across Sermilik Fjord: Helheim is terminus‐driven, Fenris and Midgard are runoff‐adapting, and Pourquoi Pas is runoff‐drivenDecomposition by empirical orthogonal functions and principal components reveals glacier type using data at all points in the glacier domainInferred glacier types differ slightly from previous work because we weight all locations equally and isolate the seasonal component [ABSTRACT FROM AUTHOR]

The largest regional drivers of current surface elevation increases in the Antarctic Ice Sheet are associated with ice flow reconfiguration in previously active ice streams, highlighting the important role of ice dynamics in responding to climate change. Here, we investigate controls on the evolution of the flow configuration of the Vanderford and Totten Glaciers – key outlet glaciers of the Aurora Subglacial Basin, the most rapidly thinning region of the East Antarctic Ice Sheet. We review factors that influence the ice flow in this region, and use an ice sheet model to investigate the sensitivity of the catchment divide location to thinning at Vanderford Glacier associated with ongoing retreat, and thickening at Totten Glacier associated with an intensification of the east-west snowfall gradient. The present-day catchment divide between the Totten and Vanderford Glaciers is not constrained by the geology or topography, but is determined by the large-scale ice sheet geometry and its long-term evolution in response to climate forcing. Furthermore, the catchment divide is subject to migration under relatively small changes in surface elevation, leading to ice flow and basal water piracy from Totten to Vanderford Glacier. Our findings show that ice flow reconfigurations do not only occur in regions of West Antarctica like the Siple Coast, but also in the east, motivating further investigations of past, and potential for future, ice flow reconfigurations around the whole Antarctic coastline. Such modelling of ice flow and basal water piracy may require coupled ice sheet thermomechanical and subglacial hydrology models, constrained by field observations of subglacial conditions. Our results also have implications for ice sheet mass budget studies that integrate over catchments, and the validity of the zero flow assumption when selecting sites for ice core records of past climate. [ABSTRACT FROM AUTHOR]

Warming of the ocean waters surrounding Greenland plays a major role in driving glacier retreat and the contribution of glaciers to sea level rise. The melt rate at the junction of the ocean with grounded ice—or grounding line—is, however, not well known. Here, we employ a time series of satellite radar interferometry data from the German TanDEM-X mission, the Italian COSMO-SkyMed constellation, and the Finnish ICEYE constellation to document the grounding line migration and basal melt rates of Petermann Glacier, a major marine-based glacier of Northwest Greenland. We find that the grounding line migrates at tidal frequencies over a kilometer-wide (2 to 6 km) grounding zone, which is one order of magnitude larger than expected for grounding lines on a rigid bed. The highest ice shelf melt rates are recorded within the grounding zone with values from 60 ± 13 to 80 ± 15 m/y along laterally confined channels. As the grounding line retreated by 3.8 km in 2016 to 2022, it carved a cavity about 204 m in height where melt rates increased from 40 ± 11 m/y in 2016 to 2019 to 60 ± 15 m/y in 2020 to 2021. In 2022, the cavity remained open during the entire tidal cycle. Such high melt rates concentrated in kilometer-wide grounding zones contrast with the traditional plume model of grounding line melt which predicts zero melt. High rates of simulated basal melting in grounded glacier ice in numerical models will increase the glacier sensitivity to ocean warming and potentially double projections of sea level rise. [ABSTRACT FROM AUTHOR]

Glacial landscapes preserve records of past climate change. Investigating the glacier–climate system over the Late Quaternary provides information about past climate change and context for present‐day glacier response to climate warming. Using 28 beryllium‐10 (10Be) surface exposure dates and snowline reconstructions, we present glacier fluctuations and climate changes for the Antarctic Cold Reversal in the Ahuriri River catchment, Southern Alps of New Zealand (44°7′50″S, 169°38′29″E). Prominent terminal and lateral moraine features from the upper right tributary of the Ahuriri River valley have exposure ages of 14.5 ± 0.3, 13.6 ± 0.3 and 12.6 ± 0.2 ka, suggesting retreat of the glacier during the Antarctic Cold Reversal. Maximum elevation of lateral moraines (MELM) and accumulation area ratio (AAR) suggest snowline elevations at these ages were ≤700, ≤630 and ~360 m lower than today, respectively. This equates to air temperatures ≤3.9, ≤3.5 and 2.3 ± 0.7 °C lower than today (1981–2010), assuming no changes in past precipitation. Ice‐sculpted bedrock surfaces bound by a lateral moraine at nearby Canyon Creek have an age of 13.1 ± 0.3 ka, indicating the moraine correlates with those in the Ahuriri upper right tributary. MELM and AAR reconstructions from the Canyon Creek suggest that snowline elevations at 14.5–13.6 ka were ≤500 or ~380 m lower than today, corresponding to air temperatures ≤2.8 or 2.4 ± 0.7 °C lower than the present‐day (1981–2010). Our results provide insight into the structure of the Antarctic Cold Reversal in the Southern Alps, showing that the largest glacier advance occurred at the start of this interval at c. 14.5 ± 0.3 ka and was followed by gradual retreat. We hypothesize that the early cooling and glacier readvance in New Zealand at the onset of the Antarctic Cold Reversal were triggered by a latitudinal shift of the Southern Hemisphere westerly wind belt. [ABSTRACT FROM AUTHOR]

Seismicity at glaciers, ice sheets and ice shelves provides observational constraint of a number of glaciological processes. Detecting and locating this seismicity, specifically icequakes, is a necessary first step in studying processes such as basal slip, crevassing, and imaging ice fabric, for example. Most glacier deployments to date use conventional seismic networks, comprised of seismometers distributed over the entire area of interest. However, smaller aperture seismic arrays can also be used, which are typically sensitive to seismicity distal from the array footprint and require a smaller number of instruments. Here, we investigate the potential of arrays and array-processing methods to detect and locate seismicity in the cryosphere, benchmarking performance against conventional seismic network-based methods. We also provide an array-processing recipe for cryosphere applications. Results from an array and network deployed at Rutford Ice Stream, Antarctica, show that arrays and networks both have strengths and weaknesses. Arrays can detect icequakes from further distances whereas networks outperform arrays for more comprehensive studies of a process within the network extent, due to greater hypocentral constraint and a smaller magnitude of completeness. We also gain new insights into seismic behaviour at Rutford Ice Stream. The array detects basal icequakes in what was previously interpreted to be an aseismic region of the bed, as well as new icequake observations at the ice stream shear-margins, where it would be challenging to deploy instruments. Finally, we make some practical recommendations for future array deployments at glaciers. [ABSTRACT FROM AUTHOR]

In 2023, the 150th anniversary of the death of Jean-Louis-Rodolphe Agassiz, a Swiss naturalist recognized for his analysis of glacier movement and his proposal of the existence of glaciations in the Pleistocene, is commemorated. Although initially rejected, his theory was eventually accepted and contributed to the systematic study of evidence of glaciations. Agassiz was recognized for his work on glaciers and his influence on the scientific community of his time. [Extracted from the article]

We reconstruct the mass balance of the Greenland Ice Sheet using a comprehensive survey of thickness, surface elevation, velocity, and surface mass balance (SMB) of 260 glaciers from 1972 to 2018. We calculate mass discharge, D, into the ocean directly for 107 glaciers (85% of D) and indirectly for 110 glaciers (15%) using velocity-scaled reference fluxes. The decadal mass balance switched from a mass gain of +47 ± 21 Gt/y in 1972-1980 to a loss of 51 ± 17 Gt/y in 1980-1990. The mass loss increased from 41 ± 17 Gt/y in 1990-2000, to 187 ± 17 Gt/y in 2000-2010, to 286 ± 20 Gt/y in 2010-2018, or sixfold since the 1980s, or 80 ± 6 Gt/y per decade, on average. The acceleration in mass loss switched from positive in 2000-2010 to negative in 2010-2018 due to a series of cold summers, which illustrates the difficulty of extrapolating short records into longer-term trends. Cumulated since 1972, the largest contributions to global sea level rise are from northwest (4.4 ± 0.2 mm), southeast (3.0 ± 0.3 mm), and central west (2.0 ± 0.2 mm) Greenland, with a total 13.7 ± 1.1 mm for the ice sheet. The mass loss is controlled at 66 ± 8% by glacier dynamics (9.1 mm) and 34 ± 8% by SMB (4.6 mm). Even in years of high SMB, enhanced glacier discharge has remained sufficiently high above equilibrium to maintain an annual mass loss every year since 1998.

Glacial detachment is a newly emerged glacial hazard type, and more than 20 cases have been identified worldwide. The Xiao Ma Glacier in the Amney Machen Mountains, northeastern Tibetan Plateau, experienced an intensive detachment of the terminal low-slope glacier in 2004, forming a large dammed lake downstream and burying 2.27 km2 of grassland. An advanced multi-phase r.avaflow model was used to simulate this glacier detachment event to reconstruct its slide mechanism and dynamic process. In total, 9.77 × 106 m3 mass detached were categorized into three different phases according to their properties: the old ice remaining in the middle layer following the 1988–1997 glacier surge (47% of total volume), new ice in the upper layer originating from the 2003–2004 surge, which directly caused the 2004 detachment (48%), and subglacial water trapped in the bedrock (5%). The modeling results showed that the glacier detachment–induced debris flow reached a peak discharge of 0.25 × 106 m3/s at 40 s at the mouth of the Xiao Ma Valley, and the entire propagation process was completed within 140 s. The initially trapped meltwater dominated the debris flow state, with a maximum velocity of 103 m/s. The underlying old glacial ice was pushed and deposited to the front of the fan, and the new ice in the upper layer was primarily deposited over the middle fan. Additionally, several important parameters in the r.avaflow model were recommended for the optimal interval of values via optimality analysis, and the potential of the model to be extensively applied for large-scale risk assessment was discussed under relatively simple settings. The modeling of the 2004 Amney Machen glacier detachment has resulted in a comprehensive understanding of the dynamic processes of such events, which is significant for future glacier disaster simulation and assessment thereby improving regional strategies in hazard early warning, prevention, and mitigation. [ABSTRACT FROM AUTHOR]

Arctic ice shelves have declined over the past several decades, one of many indications of a rapidly changing cryosphere. Here we use a collection of off-nadir Landsat 8 images, a 1978 digital orthophotograph and photogrammetrically derived DEM, satellite altimetry and other data to examine the causes of an Arctic ice-shelf retreat in northernmost Greenland, the Hunt Fjord Ice Shelf (HFIS). HFIS has several distinct provenance regions comprised of glacier-derived ice and corrugated multi-decadal fast ice, with varying ice thicknesses (5–64 m). Available imagery shows little change in HFIS between 1978 and 2012, after which several midsummer calving events occurred (2012, 2016 and 2019) that reduced the HFIS by 42.5 km2 (~56%). Shelf area losses began as the number of surface melt days on the adjacent ice sheet more than doubled relative to the 1980s. Recent calving events also occurred during open-water periods at the ice-shelf front. Prior to mid-2012, there were no calving events during similar open-water periods. HFIS tributary glaciers have thinned by 3–20 m near their grounding zones, and may have accelerated since the 1980s, likely due to increased basal melting from contact with warm Atlantic Water. [ABSTRACT FROM AUTHOR]

High-resolution, downscaled climate model data are used in a wide variety of applications across environmental sciences. Here we introduce a new, high-resolution dataset, CHELSA-TraCE21k. It is obtained by downscaling TraCE-21k data, using the "Climatologies at high resolution for the earth's land surface areas" (CHELSA) V1.2 algorithm with the objective to create global monthly climatologies for temperature and precipitation at 30 arcsec spatial resolution in 100-year time steps for the last 21 000 years. Paleo-orography at high spatial resolution and for each time step is created by combining high-resolution information on glacial cover from current and Last Glacial Maximum (LGM) glacier databases and interpolations using data from a global model of glacial isostasy (ICE-6G_C) and a coupling to mean annual temperatures from TraCE21k (Transient Climate Evolution of the last 21 000 years) based on the Community Climate System Model version 3 (CCSM3). Based on the reconstructed paleo-orography, mean annual temperature and precipitation were downscaled using the CHELSA V1.2 algorithm. The data were validated by comparisons with the glacial extent of the Laurentide ice sheet based on expert delineations, proxy data from Greenland ice cores, historical climate data from meteorological stations, and a dynamic simulation of species distributions throughout the Holocene. Validations show that the CHELSA-TraCE21k V1.0 dataset reasonably represents the distribution of temperature and precipitation through time at an unprecedented 1 km spatial resolution, and simulations based on the data are capable of detecting known LGM refugia of species. [ABSTRACT FROM AUTHOR]

When Ötzi, the Iceman, was found in a gully in the Tisenjoch pass in the Tyrolean Alps in 1991, he was a huge surprise for the archaeological community. The lead initial investigator of the find argued that it was unique, preserved by serendipitous circumstances. It was hypothesised that the mummy with associated artefacts had been quickly covered by glacier ice and stayed buried until the melt-out in 1991. It is now more than 30 years since Ötzi appeared. In this paper, we take a closer look at how the find can be understood today, benefitting from increased knowledge gained from more than two decades of investigations of other glacial archaeological sites, and from previous palaeo-biological investigations of the find assemblage. In the light of radiocarbon dates from the gully and new glaciological evidence regarding mass balance, it is likely that Ötzi was not permanently buried in ice immediately after his death, but that the gully where he lay was repeatedly exposed over the next 1500 years. We discuss the nature of the ice covering the site, which is commonly described as a basally sliding glacier. Based on the available evidence, this ice is better understood as a non-moving, stationary field of snow and ice, frozen to the bedrock. The damaged artefacts found with Ötzi were probably broken by typical postdepositional processes on glacial archaeological sites, and not, as previously claimed, during conflict prior to Ötzi's flight from the valley below. [ABSTRACT FROM AUTHOR]

Subglacial channels are morphologically and morphometrically distinct in comparison to fluvial channels, yet their identification from remote sensing data is still problematic. To contribute to the current set of criteria used to identify such channels, we performed a detailed analysis of two subglacial channel networks on Devon Island, Nunavut, Canada. In planform, these channels are isolated, finger-like networks that drain into a main stem and have distinct cross-sectional and longitudinal profiles. Cross-sections are flat-bottomed with steep walls and longitudinal profiles are convex and exhibit undulations, typical of pressurized water flow (i.e., subglacial flow). To facilitate remote sensing identification, we interrogated how well-known scaling relationships capturing hydraulics and mass balance dynamics of fluvial systems differ in subglacial channels. Scaling relationships typically used to discern connections between discharge and channel and catchment size in fluvial systems were applied to both networks, yielding trends distinct from the fluvial literature. We suggest that the weakly correlated relationship we found between channel discharge and the size of the drainage area indicates a discrete point or line source of water, such as a moulin or crevasse. • Subglacial channels related to polythermal glaciers are morphologically distinct. • Channels form networks display longitudinal undulations and U-shape cross-section. • Fluvial scaling relationships applied for the first time to subglacial channels • Scaling relationship reveals metric indicative of lack of drainage basin. • Metric may be useful as additional diagnostic criteria for subglacial meltwater channels. [ABSTRACT FROM AUTHOR]

Links between hydrology and sliding of the Greenland Ice Sheet (GrIS) are poorly understood. Here, we monitored meltwater's propagation through the glacial hydrologic system for catchments at different elevations by quantifying the lag cascade as daily meltwater pulses traveled through the supraglacial, englacial, and subglacial drainage systems. We found that meltwater's residence time within supraglacial catchments—depending upon area, snow cover, and degree of channelization—controls the timing of peak moulin head, resulting in the 2 hr later peak observed at higher elevations. Unlike at lower elevations where peak moulin head and peak sliding coincided, at higher elevations peak sliding lagged peak moulin head by ∼2.8 hr. This delay was likely caused by the area's lower moulin density, which required diurnal pressure oscillations to migrate further into the distributed drainage system to elicit the observed velocity response. These observations highlight the supraglacial drainage system's control on coupling GrIS subglacial hydrology and sliding. Plain Language Summary: Each summer, melting snow and ice collects within streams and rivers on the Greenland Ice Sheet's surface until reaching the bed through crevasses or moulins—near‐vertical conduits that penetrate the entire ice thickness—where this meltwater can lubricate the bed, causing the overlying ice to slide more rapidly. Despite the important role of meltwater in modulating sliding speeds, little is known about how relationships between melting and sliding vary spatially or through time. Here, we take the novel approach of monitoring meltwater's propagation through the entire glacial hydraulic system at two elevations. We find that longer delays in the timing of meltwater delivery to moulins draining larger, higher‐elevation catchments, caused peak moulin water level (i.e., peak pressurization) to occur 2 hr later in the day than at smaller, lower‐elevation catchments. Unlike at lower elevations where peak moulin water level and sliding coincided, at higher elevations sliding lagged peak moulin water level by 2.8 hr. This delay was likely caused by the fewer number of moulins which require a single moulin to pressurize a larger area. This work reveals the importance of the supraglacial drainage system in controlling the timing of meltwater reaching the bed and its relationship with sliding. Key Points: Larger catchments within the Greenland Ice Sheet's ablation area impart significant delays on the timing of meltwater delivery to moulinsLonger delays in meltwater delivery to moulins caused peak moulin head to occur 1–3.25 hr later in the day at higher elevationsPeak moulin head and peak sliding speeds are not coincident at higher elevations where moulin density is low [ABSTRACT FROM AUTHOR]

Contemporary and significant spatiotemporal changes in glaciers are a result of rapidly evolving regional and global climate, and continuous monitoring is imperative for understanding the complexities of glacio-climatic interactions. Another glacier parameter that has often been correlated with mass changes is glacier snow line altitude, prompting for its seasonal and annual monitoring. The effectiveness of digital elevation models (DEMs) in estimating geodetic glacier mass balances becomes impacted by the presence of voids in them, and Seehaus et al. [[9]] targeted this issue of void filling using several approaches for generating more reliable glacier elevation change datasets. [Extracted from the article]

Bed topography is a critical parameter for determining the modern-day and future dynamics of ice sheets and their outlet glaciers. This is because the topography controls the state of stress for glaciers. At glacier termini, topography can influence the timing of terminus retreat by controlling access to warm ocean waters and/or by influencing the ability of a glacier terminus to retreat over bed bumps (moraines). Inland from the terminus, the topography can also influence where glacier retreat and thinning can stabilize. In part, this is because of knickpoints in bed topography created through glacial erosion that may influence the extent to which thinning can diffuse inland for an individual glacier and thus, the timing and magnitude of long-term mass loss. Here we provide a review of the current literature on these topics. While much of the reviewed literature assumes that topography is stable on relevant timescales to humans, new research suggests that topography may change much faster than previously thought and this further complicates our ability to project future outlet glacier change. [ABSTRACT FROM AUTHOR]

Ice cliffs are features along ice sheet margins, along tropical mountain glaciers, at termini of mountain glaciers and on debris-covered glacier tongues, that have received scattered attention in literature. They cover small relative areas of glacier or margin surface respectively, but have been involved in two apparent anomalies. On the one hand, they have been identified as potential hotspots of extreme melt rates on debris-covered tongues contributing to their relatively rapid ablation, compared to the surrounding glacier surface. On the other hand, they appear where the ice margin is stable (or temporarily advancing) even under conditions of negative mass balance. In this manuscript, we recapitulate why ice cliffs remain interesting features to investigate and what we know about them so far. We conclude by suggesting to further investigate their genesis and variable morphology and their potential as windows into past climates and processes. [ABSTRACT FROM AUTHOR]

Seismic signals generated by iceberg calving can be used to monitor ice loss at tidewater glaciers with high temporal resolution and independent of visibility. We combine the empirical matched field (EMF) method and machine learning using convolutional neural networks (CNNs) for calving event detection at the Spitsbergen (SPITS) seismic array and the single broad-band station KBS on the Arctic Archipelago of Svalbard. EMF detection with seismic arrays seeks to identify all signals generated by events in a confined target region similar to single P and/or S phase templates by assessing the beam power obtained using empirical phase delays between the array stations. The false detection rate depends on threshold settings and therefore needs appropriate tuning or, alternatively, post-processing. We combine the EMF detector at the SPITS array, as well as an STA/LTA (short term average/long term average) detector at the KBS station, with a post-detection classification step using CNNs. The CNN classifier uses waveforms of the three-component record at KBS as input. We apply the methodology to detect and classify calving events at tidewater glaciers close to the KBS station in the Kongsfjord region in Northwestern Svalbard. In a previous study, a simpler method was implemented to find these calving events in KBS data, and we use it as the baseline in our attempt to improve the detection and classification performance. The CNN classifier is trained using classes of confirmed calving signals from four different glaciers in the Kongsfjord region, seismic noise examples and regional tectonic seismic events. Subsequently, we process continuous data of six months in 2016. We test different CNN architectures and data augmentations to deal with the limited training data set available. Targeting Kronebreen, one of the most active glaciers in the Kongsfjord region, we show that the best performing models significantly improve the baseline classifier. This result is achieved for both the STA/LTA detection at KBS followed by CNN classification, as well as EMF detection at SPITS combined with a CNN classifier at KBS, despite of SPITS being located at 100 km distance from the target glacier in contrast to KBS at 15 km distance. Our results will further increase confidence in estimates of ice loss at Kronebreen derived from seismic observations which in turn can help to better understand the impact of climate change in Svalbard. [ABSTRACT FROM AUTHOR]

The monitoring of glaciers in Switzerland has a long tradition, yet glacier changes during the 20th century are only known through sparse observations. Here, we estimate a halving of Swiss glacier volumes between 1931 and 2016 by mapping historical glacier elevation changes at high resolution. Our analysis relies on a terrestrial image archive known as TerrA, which covers about 86 % of the Swiss glacierised area with 21 703 images acquired during the period 1916–1947 (with a median date of 1931). We developed a semi-automated workflow to generate digital elevation models (DEMs) from these images, resulting in a 45 % total glacier coverage. Using the geodetic method, we estimate a Swiss-wide glacier mass balance of -0.52 ± 0.09 m w.e. a -1 between 1931 and 2016. This equates to a 51.5 ± 8.0 % loss in glacier volume. We find that low-elevation, high-debris-cover, and gently sloping glacier termini are conducive to particularly high mass losses. In addition to these glacier-specific, quasi-centennial elevation changes, we present a new inventory of glacier outlines with known timestamps and complete attributes from around 1931. The fragmented spatial coverage and temporal heterogeneity of the TerrA archive are the largest sources of uncertainty in our glacier-specific estimates, reaching up to 0.50 m w.e. a -1. We suggest that the high-resolution mapping of historical surface elevations could also unlock great potential for research fields other than glaciology. [ABSTRACT FROM AUTHOR]

Artificial glacier melt reduction is gaining increasing attention because of rapid glacier retreats and the projected acceleration of future mass losses. However, quantifying the effect of artificial melt reduction on glaciers in China has not been currently reported. Therefore, the case of Urumqi Glacier No.1 (eastern Tien Shan, China) is used to conduct a scientific evaluation of glacier cover efficiency for melt reduction between 24 June and 28 August 2021. By combining two high-resolution digital elevation models derived from terrestrial laser scanning and unmanned aerial vehicles, albedo, and meteorological data, glacier ablation mitigation under three different cover materials was assessed. The results revealed that up to 32% of mass loss was preserved in the protected areas compared with that of the unprotected areas. In contrast to the unprotected glacier surface, the nanofiber material reduced the glacier melt by up to 56%, which was significantly higher than that achieved by geotextiles (29%). This outcome could be attributed to the albedo of the materials and local climate factors. The nanofiber material showed higher albedo than the two geotextiles, dirty snow, clean ice, and dirty ice. Although clean snow had a higher albedo than the other materials, its impact on slowing glacier melt was minor due to the lower snowfall and relatively high air temperature after snowfall in the study area. This indicates that the efficiencies of nanofiber material and geotextiles can be beneficial in high-mountain areas. In general, the results of our study demonstrate that the high potential of glacier cover can help mitigate issues related to regions of higher glacier melt or lacking water resources, as well as tourist attractions. [ABSTRACT FROM AUTHOR]

Climate warming has become indisputable, and it is now crucial to increase our understanding of both the mechanisms and consequences of climate change. The Antarctic region is subjected to substantial changes, the trends of which have been recognized for several decades. In the South Shetland Islands, the most visible effect of climate change is progressive deglaciation. The following review focuses on past glaciological studies conducted on King George Island (KGI). The results of collected cryosphere element observations are discussed herein in a comprehensive manner. Our analysis showed that there is a lack of temporal as well as spatial continuity for studies on the basic mass balance parameters on the entire KGI ice dome and only Bellingshausen Dome has a relatively long history of data collection. The methodologies of past work, which have improved over time, are also discussed. When studying the glacier front fluctuations, the authors most frequently use a 1956 aerial photography as reference ice coverage. This was the case for seven papers, while other sources are seldomly mentioned. In other papers as many as 41 other sources were used, and therefore comparison to photos taken up to 60 years later can give misleading trends, as small glaciers may have both advanced and retreated in that time. In the case of glacial velocities there is also an apparent lack of consistency, as different glaciers were indicated as the fastest on KGI. Only Lange, Anna, Crystal, Eldred, and eastern part of Usher glaciers were determined by more than one author as the fastest. Additionally, there are gaps in the KGI Ground Penetrating Radar (GPR) survey area, which includes three ice domes: the Warszawa Icefield, the Krakow Icefield, and eastern part of King George Island. Ideas for further work on the topic are also suggested, allowing for easier access to data and thus contributing to a better understanding of glacier development mechanisms. [ABSTRACT FROM AUTHOR]

Numerical models have become an indispensable tool for understanding and predicting the flow of ice sheets and glaciers. Here we present the full-Stokes software package Underworld to the glaciological community. The code is already well established in simulating complex geodynamic systems. Advantages for glaciology are that it provides a full-Stokes solution for elasto-visco-plastic materials and includes mechanical anisotropy. Underworld uses a material point method to track the full history information of Lagrangian material points, of stratigraphic layers and of free surfaces. We show that Underworld successfully reproduces the results of other full-Stokes models for the benchmark experiments of the ISMIPHOM project. Furthermore,we test FE meshes with different geometries, and highlight the need to be able to adapt the finiteelement grid to discontinuous interfaces between materials with strongly different properties, such as the ice-bedrock boundary. [ABSTRACT FROM AUTHOR]

Our ability to improve prognostic modeling of the Greenland Ice Sheet relies on understanding the long-term relationships between climate and mass flux (via iceberg calving) from marine-terminating tidewater glaciers (TWGs). Observations of recent TWG behavior are widely available, but long-term records of TWG advance are currently lacking. We present glacial geomorphological, sedimentological, archaeological, and modeling data to reconstruct the -20 km advance of Kangiata Nunaata Sermia (KNS; the largest tidewater glacier in southwest Greenland) during the first half of the past millennium. The data show that KNS advanced -15 km during the 12th and 13th centuries CE at a rate of -115 m a-1, contemporaneous with regional climate cooling toward the Little Ice Age and comparable to rates of TWG retreat witnessed over the past -200 years. Presence of Norse farmsteads proximal to KNS demonstrates their resilience to climate change, manifest as a rapidly advancing TWG in a cooling climate. The results place limits on the magnitude of ice-margin advance and demonstrate TWG sensitivity to climate cooling as well as warming. These data combined with our grounding-line stability analysis provide a long-term record that validates approaches to numerical modeling aiming to link calving to climate. [ABSTRACT FROM AUTHOR]

We present the first fully coupled 3D full-Stokes model of a tidewater glacier, incorporating ice flow, subglacial hydrology, plume-induced frontal melting and calving. We apply the model to Store Glacier (Sermeq Kujalleq) in west Greenland to simulate a year of high melt (2012) and one of low melt (2017). In terms of modelled hydrology, we find perennial channels extending 5 km inland from the terminus and up to 41 and 29 km inland in summer 2012 and 2017, respectively. We also report a hydrodynamic feedback that suppresses channel growth under thicker ice inland and allows water to be stored in the distributed system. At the terminus, we find hydrodynamic feedbacks exert a major control on calving through their impact on velocity. We show that 2012 marked a year in which Store Glacier developed a fully channelised drainage system, unlike 2017, where it remained only partially developed. This contrast in modelled behaviour indicates that tidewater glaciers can experience a strong hydrological, as well as oceanic, control, which is consistent with observations showing glaciers switching between types of behaviour. The fully coupled nature of the model allows us to demonstrate the likely lack of any hydrological or ice-dynamic memory at Store Glacier. [ABSTRACT FROM AUTHOR]

Since the discovery of thick ice deposits on the polar caps of Mars, scientists have proposed that they flow as glaciers or ice sheets, and simultaneously, others have suggested that erosional activity, primarily by wind, has been the main driver shaping them. This debate began in 1971, and following those original interpretations, numerous studies have suggested that the PLDs flow, up to 1 m per year at steep marginal scarps, and numerous subsequent studies have negated those interpretations. Further, each time a testable prediction has been made in support of flow, is has been falsified by observations. To date, no flow has been observed at the PLDs. This creates an enduring problem; according to conventional wisdom, the polar caps "should" be flowing, but they have not in any appreciable, measurable way. Here, I review the evidence for and against flow and offer and test four hypotheses to retard flow of the NPLD based. The first three hypotheses overpredict strain and fail to account for other observations. The fourth hypothesis, based on stratified layers of varying rheology, successfully explains why the polar ice has no measurable flow. The polar ice caps of Mars together comprise as much ice as the Greenland Ice Sheet on Earth. They reach 2–3 km thick and span >1000 km. On Earth, where ice is as many as 100° warmer, places with ice ∼1 km thick flow at a measurable rate. However, at the poles of Mars, all predictions of a flowing ice sheet have been negated by observations – any flow would have to be so miniscule as to have no effect on the surface or be detectable in the stratigraphy. Mars does have abundant evidence of past flowing ice in the mid-latitudes, mostly as glaciers, and the poles have experienced temperatures warm enough for flow in the recent past, so explaining this lack of evidence for flow has eluded scientists. This manuscript proposes four ideas to inhibit flow and finds one based on stacking layers of alternating viscosity to successfully match observations. The more resistant layers act to reduce overall flow to an undetectable level. • The polar ice on Mars "should" flow; however, all past interpretations of flowing polar ice on Mars have been invalidated • I provide and test four hypotheses to explain the lack of observations supporting flow over 50+ years • Stratified materials that impede bulk flow explains the observations better than all past interpretations or models [ABSTRACT FROM AUTHOR]