Ice Front Research Articles

Topographic modulation of outlet glaciers in Greenland: a review

AbstractBed 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.

The Influence of Pine Island Ice Shelf Calving on Basal Melting

AbstractThe combination of the Pine Island Ice Shelf (PIIS) draft and a seabed ridge beneath it form a topographic barrier, restricting access of warm Circumpolar Deep Water to a cavity inshore of the ridge, thus exerting an important control on PIIS basal ablation. In addition, PIIS has recently experienced several large calving events and further calving could significantly alter the cavity geometry. Changes in the ice front location, together with changes in ice thickness, might relax the topographic barrier and thus significantly change basal melt rates. Here, we consider the impact of past, and possible future, calving events on PIIS melt rates. We use a high‐resolution ocean model to simulate melt rates in both an idealized domain whose geometry captures the salient features of Pine Island Glacier, and a realistic geometry accurately resembling it, to explore how calving affects melt rates. The idealized simulations reveal that the melt response to calving has a sensitive dependence on the thickness of the gap between PIIS and the seabed ridge and inform our interpretation of the realistic simulations, which suggest that PIIS melt rates did not respond significantly to recent calving. However, the mean melt rate increases approximately linearly with further calving, and is amplified by approximately 10% relative to present day once the ice front reaches the ridge‐crest, taking less than one decade if calving maintains its present rate. This provides strong evidence that calving may represent an important, but as yet unexplored, contribution to the ice‐ocean sensitivity of the West Antarctic Ice Sheet.

Steep ice – progress and future challenges in research on ice cliffs

AbstractIce 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.

Progress toward globally complete frontal ablation estimates of marine-terminating glaciers

AbstractKnowledge of frontal ablation from marine-terminating glaciers (i.e., mass lost at the calving face) is critical for constraining glacier mass balance, improving projections of mass change, and identifying the processes that govern frontal mass loss. Here, we discuss the challenges involved in computing frontal ablation and the unique issues pertaining to both glaciers and ice sheets. Frontal ablation estimates require numerous datasets, including glacier terminus area change, thickness, surface velocity, density, and climatic mass balance. Observations and models of these variables have improved over the past decade, but significant gaps and regional discrepancies remain, and better quantification of temporal variability in frontal ablation is needed. Despite major advances in satellite-derived large-scale datasets, large uncertainties remain with respect to ice thickness, depth-averaged velocities, and the bulk density of glacier ice close to calving termini or grounding lines. We suggest ways in which we can move toward globally complete frontal ablation estimates, highlighting areas where we need improved datasets and increased collaboration.

Rapid retreat of Thwaites Glacier in the pre-satellite era

Understanding the recent history of Thwaites Glacier, and the processes controlling its ongoing retreat, is key to projecting Antarctic contributions to future sea-level rise. Of particular concern is how the glacier grounding zone might evolve over coming decades where it is stabilized by sea-floor bathymetric highs. Here we use geophysical data from an autonomous underwater vehicle deployed at the Thwaites Glacier ice front, to document the ocean-floor imprint of past retreat from a sea-bed promontory. We show patterns of back-stepping sedimentary ridges formed daily by a mechanism of tidal lifting and settling at the grounding line at a time when Thwaites Glacier was more advanced than it is today. Over a duration of 5.5 months, Thwaites grounding zone retreated at a rate of >2.1 km per year—twice the rate observed by satellite at the fastest retreating part of the grounding zone between 2011 and 2019. Our results suggest that sustained pulses of rapid retreat have occurred at Thwaites Glacier in the past two centuries. Similar rapid retreat pulses are likely to occur in the near future when the grounding zone migrates back off stabilizing high points on the sea floor.

Spaceborne InSAR mapping of landslides and subsidence in rapidly deglaciating terrain, Glacier Bay National Park and Preserve and vicinity, Alaska and British Columbia

The Glacier Bay area in southeastern Alaska and British Columbia, encompassing Glacier Bay National Park and Preserve, has experienced rapid glacier retreat since the end of the Little Ice Age in the mid-1800s. The impact that rapid deglaciation has had on the slope stability of valley walls and on the sedimentation of fans and deltas adjacent to fjords and inlets is an ongoing research topic. Using 3-year (2018–2020) Sentinel-1 datasets, and an automated time-series persistent scatterer interferometric synthetic aperture radar (PSInSAR) processing method, we detected landslides or delta subsidence at 27 sites within a vast 180 × 180 km remote region encompassing Glacier Bay proper. Most of the sites that we identified had not been previously identified. We categorized the hazard source areas that we identified into three general types:1) slow-moving landslides on steep rocky slopes not near (> 2 km away from) present-day glacier termini (e.g., Tidal Inlet landslide), 2) slow-moving landslides directly adjacent to (< 2 km away from), and associated with glacier thinning and retreat, and 3) subsidence of glacial outwash fan deltas. In categories 1 and 2, we observed 22 landslides moving at velocities ranging from 0.5 to 4 cm/yr. In category 3, we detected five fan deltas subsiding at velocities ranging from 0.5 to 6 cm/yr. Within our measurement error, these velocities were consistent during the monitoring period. Because acceleration was not observed, the issuance of warnings of imminent rapid failure is not warranted, however, continued remote monitoring is warranted. Our interferometric synthetic aperture radar (InSAR) results could be combined with other data sets including field observations, subaerial and submarine landslide inventories, bedrock fabric mapping from newly available light detection and ranging (lidar) data, and geologic maps to produce an inherent susceptibility map for landslides in bedrock and fan deltas. This map could be used to forecast susceptibility for both earthquake and climatically induced landslides.

Influence of nanoparticles on freezing inside container equipped with fins

With loading of different shapes of nanoparticles, the solidification speed can be changed which was scrutinized in current work. Although the nanoparticles dispersion can decline the heat capacity, the conduction mode can be improved with such technique and changing the styles of nano-powders can alter the strength of conduction. The velocity terms were neglected in freezing, thus, the main equations include two equations with unsteady form for scalars of solid fraction and temperature. Grid adaption with position of ice front has been considered in simulations utilizing FEM. The upper sinusoidal and inner rectangular walls maintain cold temperature and freezing starts from these regions. Adding nanomaterial can expedite the process around 15.75% (for m = 4.8) and 29.8% (for m = 8.6). Also, utilizing particles with shapes of blade form can augment the freezing rate around 16.69%. The efficacy of m on freezing process rises around 4% with elevate of concentration of nanoparticles.

Modelled dynamic retreat of Kangerlussuaq Glacier, East Greenland, strongly influenced by the consecutive absence of an ice mélange in Kangerlussuaq Fjord

AbstractMass loss at the Greenland Ice Sheet is influenced by atmospheric processes controlling its surface mass balance, and by submarine melt and calving where glaciers terminate in fjords. There, an ice mélange – a composite matrix of calved ice bergs and sea ice – may provide a buttressing force on a glacier terminus and control terminus dynamics. Kangerlussuaq Glacier is a major outlet of the Greenland Ice Sheet, for which recent major retreat events in 2004/2005 and 2016–2018 coincided with the absence of an ice mélange in Kangerlussuaq Fjord. To better understand the response of Kangerlussuaq Glacier to climatic and oceanic drivers, a 2D flowline model is employed. Results indicate that an ice mélange buttressing force exerts a major control on calving frequency and rapid retreat. When an ice mélange forms in Kangerlussuaq Fjord, it provides stabilising forces and conditions favourable for winter terminus re-advance. When it fails to form during consecutive years, model results indicate that Kangerlussuaq Glacier is primed to retreat into the large overdeepenings in Kangerlussuaq Fjord, and to terminus positions more than 30 km farther inland, implying that excessive mass loss from Kangerlussuaq Glacier by the year 2065 cannot be excluded.

Halving of Swiss glacier volume since 1931 observed from terrestrial image photogrammetry

Abstract. 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.

Interactive movements of outlet glacier tongue and landfast sea ice in Lützow-Holm Bay, East Antarctica, detected by ALOS-2/PALSAR-2 imagery

The Shirase Glacier, one of the fastest-flowing outlet glaciers in East Antarctica, flows into a bay that is usually covered by landfast sea ice. Although the presence of landfast ice is considered to stabilize glacier tongue, quantitative assessments are insufficient. We investigated the spatiotemporal variations in the flow velocities of both the glacier tongue and the surrounding landfast ice using a correlation method with ALOS-2/PALSAR-2 images, before and after a landfast-ice collapse in April 2017. After the collapse, the glacier tongue flow accelerated by 8–15% (0.20–0.38 km a−1). After reestablishment of landfast ice, the flow velocity of the landfast ice in contact with the glacier terminus is 20%–90% of the flow velocity of the glacier terminus The estimates of interactive changes in the motions of both glacier ice and sea ice can be used to quantify the buttress effect, which is a key factor affecting the mass budget of Antarctic outlet glaciers.

TermPicks: a century of Greenland glacier terminus data for use in scientific and machine learning applications

Abstract. Marine-terminating outlet glacier terminus traces, mapped from satellite and aerial imagery, have been used extensively in understanding how outlet glaciers adjust to climate change variability over a range of timescales. Numerous studies have digitized termini manually, but this process is labor intensive, and no consistent approach exists. A lack of coordination leads to duplication of efforts, particularly for Greenland, which is a major scientific research focus. At the same time, machine learning techniques are rapidly making progress in their ability to automate accurate extraction of glacier termini, with promising developments across a number of optical and synthetic aperture radar (SAR) satellite sensors. These techniques rely on high-quality, manually digitized terminus traces to be used as training data for robust automatic traces. Here we present a database of manually digitized terminus traces for machine learning and scientific applications. These data have been collected, cleaned, assigned with appropriate metadata including image scenes, and compiled so they can be easily accessed by scientists. The TermPicks data set includes 39 060 individual terminus traces for 278 glaciers with a mean of 136 ± 190 and median of 93 of traces per glacier. Across all glaciers, 32 567 dates have been digitized, of which 4467 have traces from more than one author, and there is a duplication rate of 17 %. We find a median error of ∼ 100 m among manually traced termini. Most traces are obtained after 1999, when Landsat 7 was launched. We also provide an overview of an updated version of the Google Earth Engine Digitization Tool (GEEDiT), which has been developed specifically for future manual picking of the Greenland Ice Sheet.

Evaluating the variability of glacial lakes in the Kyrgyz and Teskey ranges, Tien Shan

The northern Tien Shan Mountains are one of the most hazardous areas in Kyrgyzstan due to frequent glacial lake outburst floods (GLOFs). The current status of glacial lakes in this region has not been reported in recent years. Therefore, this study investigates the variations of glacial lakes in the Kyrgyz and Teskey ranges using Sentinel-2 and PlanetScope optical satellite images from 2017 to 2019. Results confirm that there are 800 glacial lakes, of which 351 were in the Kyrgyz and 449 were in the Teskey Range. We found seasonality in the character of glacial lakes, most of which have been found to be unstable during our observation period. 242 glacial lakes showed significant variations in area. We found 46 newly emerging glacial lakes, five of which have rapidly expanded in the year since they formed, while 18 were found to have vanished. Thirteen lakes were found to be short-lived or unstable. In this study, 14 lakes showed a significant expansion of area. Such large variations in the glacial lakes are closely linked to local geomorphological conditions, local climate and glacier melt. We also identified a significant shift in the number of glacial lakes at altitudes over 3,500 m asl. It is associated with the mountain permafrost zone and variations of the glacier terminus. The current study recommends detailed investigation of glacial lakes in each season with the use of remote sensing and field surveys.

Maritime glacier retreat and terminus area change in Kenai Fjords National Park, Alaska, between 1984 and 2021

AbstractGlacier change in Kenai Fjords National Park in southcentral Alaska affects local terrestrial, fresh water and marine ecosystems and will likely impact ecotourism. We used Landsat 4–8 imagery from 1984 through 2021 to manually map lower glacier ice margins for 19 maritime glaciers in Kenai Fjords National Park. Of these glaciers, six are tidewater, three are lake-terminating, six are land-terminating and four terminated in more than one environment throughout the study period. We used the mapped ice margins to quantify seasonal terminus position and areal change, including distinguishing between ice loss at glacier termini and along glacier margins. Overall, 13 glaciers substantially retreated (more than 2σ), 14 lost substantial area and only two underwent both net advance and area gain. The glaciers that had insubstantial length and area changes were predominantly tidewater. Cumulatively, the lower reaches of these 19 glaciers lost 42 km2 of ice, which was nearly evenly distributed between the terminus and the lateral margins. The rapid rate of glacier change and subsequent land cover changes are highly visible to visitors and locals at Kenai Fjords National Park, and this study quantifies those changes in terms of glacier length and area.

Monitoring glacier terminus and surface velocity changes over different time scales using massive imagery analysis and offset tracking at the Hoh Xil World Heritage Site, Qinghai-Tibet Plateau

Due to climate change, the Qinghai Hoh Xil World Heritage Site is under a significant threat from glacier melting and related glaciological hazards. A comprehensive understanding of glacier change trends will be important for launching the policy of sustainable development and for taking precautionary actions when abnormal glacier movements are recognized. Current research neglects the phenomenon in glacier movement monitoring, where the glacier variation pattern responds abnormally over different time scales. In our work, the most advanced techniques, Google Earth Engine (GEE) and offset tracking (OT), were utilized to construct and compare long-term series of glacier changes over decadal and annual scales, respectively. In total, 264 glaciers were extracted, and similar trends, namely, general degradation except one glacier expressing a fast increase in XinQing Feng (XQF), slow recession in Malan (ML), and complex and seasonal degradation in Galadandong (GLDD), were confirmed by GEE and OT. By comparing the performance of these two methods, GEE shows its merits when applied to detect the large magnitude of geomorphic and glaciological change with an average uncertainty error of approximately 8 m/yr, while OT is more suitable for detecting regular and continuous displacement with high accuracy up to 1 m/yr. Our results highlight the asynchrony of glacier variation within different time scales and the significant complementary role of integrating GEE and OT approaches for world heritage site conservation and protection.

Simulation of physical behavior of nanomaterial during freezing within a container

Galerkin method for scrutinizing the treatment of system during solidification of water has been implemented in existence of alumina nanoparticle. Cylinder and blade shapes of nanoparticles have been dispersed in to water in considering concentration less than 0.04. With neglecting the gravity force impact during freezing process, there exist only energy equations which should involve phase change term as mathematical model. Such equation has been coupled with equation for fraction of solid phase and final equations have been solved numerically. The configuration of grid has been considered variable according to position of ice front. As amount of [Formula: see text] increases, the freezing period declines around 10.7% owing to greater conductivity of Nanoparticle Enhanced Phase Change Material (NEPCM). With dispersing nanoparticles, the required time for solidification declines around 47.7%.

Structural evolution during a surge in the Paulabreen glacier system, Svalbard

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

A probabilistic framework for quantifying the role of anthropogenic climate change in marine-terminating glacier retreats

Abstract. Many marine-terminating outlet glaciers have retreated rapidly in recent decades, but these changes have not been formally attributed to anthropogenic climate change. A key challenge for such an attribution assessment is that if glacier termini are sufficiently perturbed from bathymetric highs, ice-dynamic feedbacks can cause rapid retreat even without further climate forcing. In the presence of internal climate variability, attribution thus depends on understanding whether (or how frequently) these rapid retreats could be triggered by climatic noise alone. Our simulations with idealized glaciers show that in a noisy climate, rapid retreat is a stochastic phenomenon. We therefore propose a probabilistic approach to attribution and present a framework for analysis that uses ensembles of many simulations with independent realizations of random climate variability. Synthetic experiments show that century-scale climate trends substantially increase the likelihood of rapid glacier retreat. This effect depends on the timescales over which ice dynamics integrate forcing. For a population of synthetic glaciers with different topographies, we find that external trends increase the number of large retreats triggered within the population, offering a metric for regional attribution. Our analyses suggest that formal attribution studies are tractable and should be further pursued to clarify the human role in recent ice-sheet change. We emphasize that early-industrial-era constraints on glacier and climate state are likely to be crucial for such studies.

Reconstruction of The Ice Front Within an Icing Droplet Using High Speed Laser Induced Fluorescence Imaging

An experimental setup is implemented in order to study the solidification of a drop impinging a subcooled substrate. One of the major problems relying the study of droplet solidification consist of the visualization of the solidification front within the droplet. Indeed, shadowgraphy measurement only allows the observation of the solidification front along the tri-junction liquid/solid/air. The use of the Laser Induced Fluorescence provides information on the evolution of the horizontal solidification front in a droplet over t ime. The images reconstruction of the solidification front is made by using two high-speed cameras (side view and top view). The measurements allow for the first time to observe the evolution of the solidification front geometry over t ime. The measurements, carried out on a duralumin substrate, were then compared with a 2D numerical model taking into account the heat transfers conjugated with the substrate.

Chlorine-36 Surface Exposure Dating of Late Holocene Moraines and Glacial Mass Balance Modeling, Monte Sierra Nevada, South-Central Chilean Andes (38°S)

The development of robust chronologies of Neoglaciation from individual glaciers throughout the high-altitude Andes can provide fundamental knowledge of influences such as regional temperature and precipitation variability, and aid in predicting future changes in the Andean climate system. However, records of Late Holocene glaciation from the Central Chilean Andes are sparse, and often poorly constrained. Here, we present 36Cl surface exposure ages, dendrochronologic constraints, and glacial mass balance modeling simulations of Late Holocene glacier fluctuations in the Central-South Chilean Andes. A series of concentric moraine ridges were identified on Monte Sierra Nevada (38°S), where exposure dating of basaltic boulders was used to establish a chronology of ice recession. We infer that moraine abandonment of the most distal ridge in the valley commenced by ∼4.2 ka, and was followed by glacier margin retreat to an up-valley position. Exposure ages of the oldest Late Holocene boulders (∼2.5–0.8 ka) along the marginal extents of the moraine complex indicate fluctuations of the glacier terminus prior to ∼0.65 ka. A final expansion of the ice margin reoccupied the position of the 4.2 ka moraine, with abatement from the outermost composite moraine occurring by ∼0.70 ka, as constrained by tree-ring data from live Araucaria araucana trees. Finally, a series of nested moraines dating to ∼0.45–0.30 ka, formed from a pulsed ice recession during the latest Holocene when the lower reaches of the glacial snout was most likely debris mantled. A distributed temperature index model combined with a glacier flow model was used to quantify an envelope of possible climatic conditions of Late Holocene glaciation. The glacial modeling results suggest conditions were ∼1.5°C colder and 20% wetter during peak Neoglaciation relative to modern conditions. These records also suggest a near-coeval record of Late Holocene climate variability between the middle and high southern latitudes. Furthermore, this study presents some of the youngest 36Cl exposure ages reported for moraines in the Andes, further supporting this method as a valuable geochronologic tool for assessing Late Holocene landscape development.

Effect of topography on the changes of Urumqi Glacier No. 1 in the Chinese Tianshan Mountains

Topography plays an important role in determining the glacier changes. However, topography has often been oversimplified in the studies of the glacier changes. No systematic studies have been conducted to evaluate the relationship between the glacier changes and topographic features. The present study provided a detailed insight into the changes in the two branches (east branch and west branch) of Urumqi Glacier No. 1 in the Chinese Tianshan Mountains since 1993 and systematically discussed the effect of topography on the glacier parameters. This study analyzed comprehensive recently observed data (from 1992/1993 to 2018/2019), including mass balance, ice thickness, surface elevation, ice velocity, terminus, and area, and then determined the differences in the changes of the two branches and explored the effect of topography on the glacier changes. We also applied a topographic solar radiation model to analyze the influence of topography on the incoming shortwave radiation (SWin) across the entire glacier, focusing on the difference in the SWin between the two branches. The glacier mass balance of the east branch was more negative than that of the west branch from 1992/1993 to 2018/2019, and this was mainly attributed to the lower average altitude of the east branch. Compared with the west branch, the decrease rate of the ice velocity was lower in the east branch owing to its relatively increased slope. The narrow shape of the west branch and its southeast aspect in the earlier period resulted in a larger glacier terminus retreat of the west branch. The spatial variability of the SWin across the glacier surface became much larger as altitude increased. The SWin received by the east branch was slightly larger than that received by the west branch, and the northern aspect could receive more SWin, leading to glacier melting. In the future, the difference of the glacier changes between the two branches will continue to exist due to their topographic differences. This work is fundamental to understanding how topographic features affect the glacier changes, and provides information for building different types of relationship between the glacier area and ice volume to promote further studies on the basin-scale glacier classification.