Search

Assessment of landslide hazards often depends on the ability to track possible changes in natural slopes. To that end, historical air photos can be useful, particularly when slope stability is compromised by visible cracking. Undocumented landsliding rejuvenates a paleo-landslide on a busy motorway connecting Jerusalem to a small Jewish settlement. Recently, a plan for broadening the motorway was approved, and we were asked to study the hazards of the road by Israeli NGOs and Palestinian residents of the area. We captured high-resolution topography around the unstable slope using drone-borne photogrammetry and LiDAR surveys. The modern data allow us to analyze historic air photos and topo maps to assess the level of sliding prior to and during modern landscaping. Our results indicate horizontal offsets of ~0.9–1.8 m and vertical offsets of 1.54–2.95 m at selected sites. We next assess the possible role of anthropogenic versus natural factors in compromising slope stability. We analyze monthly rain records together with seismic catalogs spanning several decades. Shortly after the motorway construction in 1995, a January 1996 rainstorm triggered a massive rockfall. The rockfall blocked traffic with up to 4 m-diameter boulders. We found that while a certain level of rain is a necessary condition for mobilizing the rock mass, it is the anthropogenic intervention that caused the rockfall in this site. We conclude that the recent plan for broadening the motorway jeopardizes the lives of vehicle passengers and the lives of future residents should the development materialize. [ABSTRACT FROM AUTHOR]

Assessment of landslide hazards often depends on the ability to track possible changes in natural slopes. To that end, historical air photos can be useful, particularly when slope stability is compromised by visible cracking. Undocumented landsliding rejuvenates a paleo-landslide on a busy motorway connecting Jerusalem to a small Jewish settlement. Recently, a plan for broadening the motorway was approved, and we were asked to study the hazards of the road by Israeli NGOs and Palestinian residents of the area. We captured high-resolution topography around the unstable slope using drone-borne photogrammetry and LiDAR surveys. The modern data allow us to analyze historic air photos and topo maps to assess the level of sliding prior to and during modern landscaping. Our results indicate horizontal offsets of ~0.9–1.8 m and vertical offsets of 1.54–2.95 m at selected sites. We next assess the possible role of anthropogenic versus natural factors in compromising slope stability. We analyze monthly rain records together with seismic catalogs spanning several decades. Shortly after the motorway construction in 1995, a January 1996 rainstorm triggered a massive rockfall. The rockfall blocked traffic with up to 4 m-diameter boulders. We found that while a certain level of rain is a necessary condition for mobilizing the rock mass, it is the anthropogenic intervention that caused the rockfall in this site. We conclude that the recent plan for broadening the motorway jeopardizes the lives of vehicle passengers and the lives of future residents should the development materialize.

Understanding glaciers structural heterogeneity is crucial for assessing their fate. Yet, places where structure changes are strong, such as crevasses fields, are often inaccessible for direct instrumentation. To overcome this limitation, we introduce an innovative technique that transforms seismic sources, here generated by crevasses, into virtual receivers using source‐to‐receiver spatial reciprocity. We demonstrate that phase interference patterns between well‐localized seismic sources can be leveraged to retrieve phase velocity maps using Seismic Michelson Interferometry. The obtained phase velocity exhibits sensitivity to changes in glacier structure, offering insights into the origins of mechanical property changes, with spatial resolution surpassing traditional methods by a factor of five. In particular, we observe sharp variations in phase velocity related to strongly damaged subsurface areas indicating a complex 3‐D medium. Applying this method more systematically and in other contexts will enhance our understanding of the structure of glaciers and other seismogenic environments. Key Points: We transform seismic sources from crevasses into virtual receivers using source‐to‐receiver spatial reciprocityWe derive phase velocity maps in previously inaccessible areas with a resolution five times larger than traditional approachesWe retrieve the influence of glacier geometry and structural heterogeneity on the glacier mechanical properties [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]

The evolution of crevasses on ice shelves plays a pivotal role in maintaining their stability. Investigating the developmental pattern of crevasses on ice shelf surfaces necessitates long-term dynamic monitoring, high-precision automatic extraction algorithms, and quantitative analysis tools. Accordingly, based on characteristics such as the long span of crevasses and their meandering distribution, we proposed Strip_Unet for the identification of crevasses. Within this, the Strip_block module utilizes strip convolution in conjunction with traditional convolution to capture the contextual information from the surrounding area, and the ASPP module is used to enhance the receptive field. Leveraging this model, we successfully identified crevasses within the damaged area of the Amery Ice Shelf from 2003 through 2023, utilizing 15-meter resolution Landsat imagery. Subsequently, we retrieved the depth details of these crevasses by applying ATL06 data, allowing us to quantitatively assess the evolution of crevasses. Our study represents the first investigation into the development of the damage area on the Amery Ice Shelf through long-term time-series analysis. Our findings indicate that crevasses originate from specific areas and follow recurring propagation cycles, exhibiting a consistent annual increase in total length. Additionally, changes in the width of these crevasses coincide with alterations in their depth. Notably, the downstream regions of the damaged area demonstrate higher crevasses density and the fastest movement speed. The proposed automatic identification of crevasses and three-dimensional dynamic monitoring hold significant potential for application throughout the Antarctic region, providing valuable insights into the long-term stability of ice shelves.

Abstract Understanding glaciers structural heterogeneity is crucial for assessing their fate. Yet, places where structure changes are strong, such as crevasses fields, are often inaccessible for direct instrumentation. To overcome this limitation, we introduce an innovative technique that transforms seismic sources, here generated by crevasses, into virtual receivers using source‐to‐receiver spatial reciprocity. We demonstrate that phase interference patterns between well‐localized seismic sources can be leveraged to retrieve phase velocity maps using Seismic Michelson Interferometry. The obtained phase velocity exhibits sensitivity to changes in glacier structure, offering insights into the origins of mechanical property changes, with spatial resolution surpassing traditional methods by a factor of five. In particular, we observe sharp variations in phase velocity related to strongly damaged subsurface areas indicating a complex 3‐D medium. Applying this method more systematically and in other contexts will enhance our understanding of the structure of glaciers and other seismogenic environments.

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

Crevasses pose severe risks for mountaineers and field glaciologists. Smaller cracks between 0.5 and 2 m are still dangerous, but often not visible in medium resolution satellite imagery. If they are snow covered, they are completely undetectable by optical sensors. We set out to develop an approach to detect potentially crevassed areas by a minimum of geometric data, and to make the method generally applicable to glacier regions. On Vernagtferner, we compared a reference dataset of crevasses observed in high-resolution optical imagery with the curvature of the ice surface and the spatial gradients in driving stress. Both parameters can be derived from a digital surface model and a bedrock model, derived from ice thickness measurements. The correlation patterns show that crevasses preferably form in convex areas and in areas where the driving stress rapidly increases. This corresponds with the theory of crevasse formation. Although the method still misclassifies larger parts, the approach has the potential to define probable non-crevassed areas as well as to aid the planning of safe routes.

In Southeast Greenland, summer melt and high winter snowfall rates give rise to firn aquifers: vast stores of meltwater buried beneath the ice-sheet surface. Previous detailed studies of a single Greenland firn aquifer site suggest that the water drains into crevasses, but this is not known at a regional scale. We develop and use a tool in Ghub, an online gateway of shared datasets, tools and supercomputing resources for glaciology, to identify crevasses from elevation data collected by NASA's Airborne Topographic Mapper across 29000 km2 of Southeast Greenland. We find crevasses within 3 km of the previously mapped downglacier boundary of the firn aquifer at 20 of 25 flightline crossings. Our data suggest that crevasses widen until they reach the downglacier boundary of the firn aquifer, implying that crevasses collect firn-aquifer water, but we did not find this trend with statistical significance. The median crevasse width, 27 meters, implies an aspect ratio consistent with the crevasses reaching the bed. Our results support the idea that most water in Southeast Greenland firn aquifers drains through crevasses. Less common fates are discharge at the ice-sheet surface (3 of 25 sites) and refreezing at the aquifer bottom (1 of 25 sites).

Mass loss from iceberg calving at marine-terminating glaciers is one of the largest and most poorly constrained contributors to sea-level rise. However, our understanding of the processes controlling ice fracturing and crevasse evolution is incomplete. Here, we use Gabor filter banks to automatically map crevasse density and orientation through time on a ~150 km2 terminus region of Narsap Sermia, an outlet glacier of the southwest Greenland ice sheet. We find that Narsap Sermia is dominated by transverse (flow-perpendicular) crevasses near the ice front and longitudinal (flow-aligned) crevasses across its central region. Measured crevasse orientation varies on sub-annual timescales by more than 45$^\circ$ in response to seasonal velocity changes, and also on multi-annual timescales in response to broader dynamic changes and glacier retreat. Our results show a gradual up-glacier propagation of the zone of flow-transverse crevassing coincident with frontal retreat and acceleration occurring in 2020/21, in addition to sub-annual crevasse changes primarily in transition zones between longitudinal to transverse crevasse orientation. This provides new insight into the dynamics of crevassing at large marine-terminating glaciers and a potential approach for the rapid identification of glacier dynamic change from a single pair of satellite images.

Basal crevasses are macroscopic structural discontinuities at the base of ice sheets and glaciers that arise by fracture. Motivated by observations and by the mechanics of elastic fracture, we hypothesise that spatial variations in basal stress (in the presence of basal water pressure) can promote and localise basal crevassing. We quantify this process in the theoretical context of linear elastic fracture mechanics. We develop a model evaluating the effect of shear-stress variation on the growth of basal crevasses. Our results indicate that sticky patches promote the propagation of basal crevasses, increase their length of propagation into the ice and, under some conditions, give them curved trajectories that incline upstream. A detailed exploration of the parameter space is conducted to gain a better understanding of the conditions under which sticky-patch-induced basal crevassing is expected beneath ice sheets and glaciers.

In Southeast Greenland, summer melt and high winter snowfall rates give rise to firn aquifers: vast stores of meltwater buried beneath the ice-sheet surface. Previous detailed studies of a single Greenland firn aquifer site suggest that the water drains into crevasses, but this is not known at a regional scale. We develop and use a tool in Ghub, an online gateway of shared datasets, tools and supercomputing resources for glaciology, to identify crevasses from elevation data collected by NASA's Airborne Topographic Mapper across 29000 km2 of Southeast Greenland. We find crevasses within 3 km of the previously mapped downglacier boundary of the firn aquifer at 20 of 25 flightline crossings. Our data suggest that crevasses widen until they reach the downglacier boundary of the firn aquifer, implying that crevasses collect firn-aquifer water, but we did not find this trend with statistical significance. The median crevasse width, 27 meters, implies an aspect ratio consistent with the crevasses reaching the bed. Our results support the idea that most water in Southeast Greenland firn aquifers drains through crevasses. Less common fates are discharge at the ice-sheet surface (3 of 25 sites) and refreezing at the aquifer bottom (1 of 25 sites). [ABSTRACT FROM AUTHOR]

Crevasses pose severe risks for mountaineers and field glaciologists. Smaller cracks between 0.5 and 2 m are still dangerous, but often not visible in medium resolution satellite imagery. If they are snow covered, they are completely undetectable by optical sensors. We set out to develop an approach to detect potentially crevassed areas by a minimum of geometric data, and to make the method generally applicable to glacier regions. On Vernagtferner, we compared a reference dataset of crevasses observed in high-resolution optical imagery with the curvature of the ice surface and the spatial gradients in driving stress. Both parameters can be derived from a digital surface model and a bedrock model, derived from ice thickness measurements. The correlation patterns show that crevasses preferably form in convex areas and in areas where the driving stress rapidly increases. This corresponds with the theory of crevasse formation. Although the method still misclassifies larger parts, the approach has the potential to define probable non-crevassed areas as well as to aid the planning of safe routes. [ABSTRACT FROM AUTHOR]

Mass loss from iceberg calving at marine-terminating glaciers is one of the largest and most poorly constrained contributors to sea-level rise. However, our understanding of the processes controlling ice fracturing and crevasse evolution is incomplete. Here, we use Gabor filter banks to automatically map crevasse density and orientation through time on a ~150 km2 terminus region of Narsap Sermia, an outlet glacier of the southwest Greenland ice sheet. We find that Narsap Sermia is dominated by transverse (flow-perpendicular) crevasses near the ice front and longitudinal (flow-aligned) crevasses across its central region. Measured crevasse orientation varies on sub-annual timescales by more than 45 $^\circ$ in response to seasonal velocity changes, and also on multi-annual timescales in response to broader dynamic changes and glacier retreat. Our results show a gradual up-glacier propagation of the zone of flow-transverse crevassing coincident with frontal retreat and acceleration occurring in 2020/21, in addition to sub-annual crevasse changes primarily in transition zones between longitudinal to transverse crevasse orientation. This provides new insight into the dynamics of crevassing at large marine-terminating glaciers and a potential approach for the rapid identification of glacier dynamic change from a single pair of satellite images. [ABSTRACT FROM AUTHOR]

Tasman Glacier, a temperate maritime glacier in the New Zealand Southern Alps, is rapidly receding. Climate warming is resulting in lengthening of the ablation season, meaning crevasses in the accumulation area are becoming exposed at the surface for longer. We combine measurements of air temperature and wind speed from inside crevasses with surface meteorological data, finding that during summer, in-crevasse air temperature is frequently positive, and can at times exceed surface air temperature. Greatest warming occurred in the widest crevasses during clear-sky conditions, but full depth warming of crevasses also occurred at night. Net shortwave radiation contributes to heating of air in the upper regions of crevasses, but turbulent sensible heat transfer was responsible for driving warm air deeper into crevasses. Crevasses orientated to maximise radiation retrieval, and running parallel to wind flow, have the greatest potential for warming and heat storage. We hypothesise a positive feedback loop in the surface energy-balance system, where crevasses entrain and trap heat, which enhances melting, that in turn enlarges the crevasses, enabling greater heat storage and further melting. Energy-balance models that treat accumulation areas of alpine glaciers as homogeneous surfaces will therefore underestimate snow melt and overestimate mass balance.

Basal crevasses are macroscopic structural discontinuities at the base of ice sheets and glaciers that arise by fracture. Motivated by observations and by the mechanics of elastic fracture, we hypothesise that spatial variations in basal stress (in the presence of basal water pressure) can promote and localise basal crevassing. We quantify this process in the theoretical context of linear elastic fracture mechanics. We develop a model evaluating the effect of shear-stress variation on the growth of basal crevasses. Our results indicate that sticky patches promote the propagation of basal crevasses, increase their length of propagation into the ice and, under some conditions, give them curved trajectories that incline upstream. A detailed exploration of the parameter space is conducted to gain a better understanding of the conditions under which sticky-patch-induced basal crevassing is expected beneath ice sheets and glaciers. [ABSTRACT FROM AUTHOR]

Iceberg calving, the process where icebergs detach from glaciers, remains poorly understood. Moreover, few parameterizations of the calving process can easily be integrated into numerical models to accurately capture observations, resulting in large uncertainties in projected sea level rise. Recent efforts have focused on estimating crevasse depths assuming tensile failure occurs when crevasses fully penetrate the glacier thickness. However, these approaches often ignore the role of advecting crevasses on calculations of crevasse depth. Here, we examine a more general crevasse depth calving model that includes crevasse advection. We apply this model to idealized prograde and retrograde bed geometries as well as a prograde geometry with a sill. Neglecting crevasse advection results in steady glacier advance and ice tongue formation for all ice temperatures, sliding law coefficients and constant slope bed geometries considered. In contrast, crevasse advection suppresses ice tongue formation and increases calving rates, leading to glacier retreat. Furthermore, crevasse advection allows a grounded calving front to stabilize on top of sills. These results suggest that crevasse advection can radically alter calving rates and hint that future parameterizations of fracture and failure need to more carefully consider the lifecycle of crevasses and the role this plays in the calving process.

Tasman Glacier, a temperate maritime glacier in the New Zealand Southern Alps, is rapidly receding. Climate warming is resulting in lengthening of the ablation season, meaning crevasses in the accumulation area are becoming exposed at the surface for longer. We combine measurements of air temperature and wind speed from inside crevasses with surface meteorological data, finding that during summer, in-crevasse air temperature is frequently positive, and can at times exceed surface air temperature. Greatest warming occurred in the widest crevasses during clear-sky conditions, but full depth warming of crevasses also occurred at night. Net shortwave radiation contributes to heating of air in the upper regions of crevasses, but turbulent sensible heat transfer was responsible for driving warm air deeper into crevasses. Crevasses orientated to maximise radiation retrieval, and running parallel to wind flow, have the greatest potential for warming and heat storage. We hypothesise a positive feedback loop in the surface energy-balance system, where crevasses entrain and trap heat, which enhances melting, that in turn enlarges the crevasses, enabling greater heat storage and further melting. Energy-balance models that treat accumulation areas of alpine glaciers as homogeneous surfaces will therefore underestimate snow melt and overestimate mass balance. [ABSTRACT FROM AUTHOR]

Snow covered crevasses in the Siachen glacier (Karakoram Range) cause great danger during glacier movement, and the knowledge of their spatial distribution is important for the safe travel. In the present study, crevasses in the Siachen glacier have been detected and further categorized as permanent open/hidden and seasonal open/hidden using a combination of optical (Landsat-8 and Sentinel-2) and microwave (ALOS-2) satellite data. The study is carried out for the year 2018. Initially, the locations of crevasses are manually marked using ALOS-2 data and further, their categorization in open and hidden is done using Sentinel-2 data. Apart from manual marking, the band ratio method is applied on Landsat-8 data to detect the permanent open crevasses in an automatic manner. Gray Level Co-occurrence Matrix (GLCM) technique has also been successfully attempted for the automatic classification of the crevasse and non-crevasse zones. The open crevasse locations using band ratioing are compared with manually marked and a good correlation is found having an accuracy of ~ 93%. A total of 140 crevasse zones have been found within the study area, with 32% as permanent open, 29% as permanent hidden, and 39% as seasonal open/hidden. Manual digitization is done for estimating crevasse dimension (length and width) as this is important input required before the start of any movement. The study reveals that the highly varying nature of crevasses in terms of changes from hidden to open in lower altitude regions of the glacier is mainly observed between July and September. The present methodology of crevasse detection and categorization leads to crevasse information system, which will be used in future to monitor the opening and closing of crevasses in the Siachen glacier on regular basis. [ABSTRACT FROM AUTHOR]

Ice cliff distribution plays a major role in determining the melt of debris‐covered glaciers but its controls are largely unknown. We assembled a data set of 37,537 ice cliffs and determined their characteristics across 86 debris‐covered glaciers within High Mountain Asia (HMA). We find that 38.9% of the cliffs are stream‐influenced, 19.5% pond‐influenced and 19.7% are crevasse‐originated. Surface velocity is the main predictor of cliff distribution at both local and glacier scale, indicating its dependence on the dynamic state and hence evolution stage of debris‐covered glacier tongues. Supraglacial ponds contribute to maintaining cliffs in areas of thicker debris, but this is only possible if water accumulates at the surface. Overall, total cliff density decreases exponentially with debris thickness as soon as the debris layer reaches a thickness of over 10 cm. Plain Language Summary: Debris‐covered glaciers are common throughout the world's mountain ranges and are characterized by the presence of steep ice cliffs among the debris‐covered ice. It is well‐known that the cliffs are responsible for a large portion of the melt of these glaciers but the controls on their formation, development and distribution across glaciers remains poorly understood. Novel mapping approaches combined with high‐resolution satellite and drone products enabled us to disentangle some of these controls and to show that the ice cliffs are generally formed and maintained by the surface hydrology (ponds or streams) or by the opening of crevasses. As a result, they depend both at the local and glacier scale on the dynamic state of the glaciers as well as the evolution stage of their debris cover. This provides a pathway to better represent their contribution to glacier melt in predictive glacier models. Key Points: We derived an unprecedented data set of 37,537 ice cliffs and their characteristics across 86 debris‐covered glaciers in High Mountain AsiaWe find that 38.9% of the cliffs are stream‐influenced, 19.5% pond‐influenced and 19.7% are crevasse‐originatedIce cliff distribution can be predicted by velocity, as an indicator of the dynamics and state of evolution of debris‐covered glaciers [ABSTRACT FROM AUTHOR]

Abstract Ice cliff distribution plays a major role in determining the melt of debris‐covered glaciers but its controls are largely unknown. We assembled a data set of 37,537 ice cliffs and determined their characteristics across 86 debris‐covered glaciers within High Mountain Asia (HMA). We find that 38.9% of the cliffs are stream‐influenced, 19.5% pond‐influenced and 19.7% are crevasse‐originated. Surface velocity is the main predictor of cliff distribution at both local and glacier scale, indicating its dependence on the dynamic state and hence evolution stage of debris‐covered glacier tongues. Supraglacial ponds contribute to maintaining cliffs in areas of thicker debris, but this is only possible if water accumulates at the surface. Overall, total cliff density decreases exponentially with debris thickness as soon as the debris layer reaches a thickness of over 10 cm.

Iceberg calving, the process where icebergs detach from glaciers, remains poorly understood. Moreover, few parameterizations of the calving process can easily be integrated into numerical models to accurately capture observations, resulting in large uncertainties in projected sea level rise. Recent efforts have focused on estimating crevasse depths assuming tensile failure occurs when crevasses fully penetrate the glacier thickness. However, these approaches often ignore the role of advecting crevasses on calculations of crevasse depth. Here, we examine a more general crevasse depth calving model that includes crevasse advection. We apply this model to idealized prograde and retrograde bed geometries as well as a prograde geometry with a sill. Neglecting crevasse advection results in steady glacier advance and ice tongue formation for all ice temperatures, sliding law coefficients and constant slope bed geometries considered. In contrast, crevasse advection suppresses ice tongue formation and increases calving rates, leading to glacier retreat. Furthermore, crevasse advection allows a grounded calving front to stabilize on top of sills. These results suggest that crevasse advection can radically alter calving rates and hint that future parameterizations of fracture and failure need to more carefully consider the lifecycle of crevasses and the role this plays in the calving process. [ABSTRACT FROM AUTHOR]

Crary Ice Rise formed after the Ross Ice Shelf re-grounded ~1 kyr BP. We present new ice-penetrating radar data from two systems operating at center frequencies of 7 and 750 MHz that confirm the ice rise is composed of a former ice shelf buried by subsequent accumulation. Stacks of englacial diffraction hyperbolas are present almost everywhere across the central ice rise and extend up to ~350 m above the bed. In many cases, bed reflections beneath the diffraction hyperbolas are obscured for distances up to 1 km. Waveform modeling indicates that the diffraction hyperbolas are likely caused by marine ice deposits in former basal crevasses and rifts. The in-filling of rifts and basal crevasses may have strengthened the connection between the ice rise and the surrounding ice shelf, which could have influenced local and regional ice dynamics. Three internal reflection horizons mark the upper limit of disturbed ice and diffraction hyperbolas in different sections of the ice rise, indicating at least three stages of flow stabilization across the ice rise. A surface lineation visible in MODIS imagery corresponds spatially to deepening and strong deformation of these layers, consistent with the characteristics of former grounding lines observed elsewhere in Antarctica.

A subglacial groundwater system beneath Taylor Glacier, Antarctica, discharges hypersaline, iron‐rich brine episodically at the glacier surface to create Blood Falls. However, the triggering mechanism for these brine release events is not yet understood. Identifying which fracture processes are observed seismically can help us better characterize the hydrological system at Taylor Glacier, and more generally, provide us with a broader understanding of englacial hydrologic activity in cold glaciers. We document wintertime brine discharge using time‐lapse photography. Subfreezing air temperatures during the brine discharge indicate that surface melt‐induced hydrofracture is an unlikely trigger for brine release. Further, we analyze local seismic data to test a hypothesis that fracturing generates elevated surface wave energy preceding and/or coinciding with brine release events. Our results show no discernible elevated Rayleigh wave activity prior to or during Blood Falls brine release. Instead, we find a pattern of seismic events dominated by a seasonal signal, with more Rayleigh events occurring in the summer than the winter from the Blood Falls source area. We calculate that the volumetric opening of cracks that would generate Rayleigh waves at our detection limits are of similar size to myriad cracks in glacier ice, lake ice, and frozen sediment in the terminus area. We therefore propose that any fracturing coincident with brine release activity likely consists of a series of smaller opening events that are masked by other seismicity in the local environment. Plain Language Summary: Blood Falls is a reddish feature that forms at the terminus of Taylor Glacier in Antarctica when hypersaline, iron‐rich brine flows from sediment under the glacier up through the ice to emerge from cracks at the surface. We used data from a time‐lapse camera and nearby seismic sensors to document a brine release. Our images show a brine release event starting in May 2014 (austral winter), and fractures in the glacier surface observed following the event encouraged us to hypothesize that we would detect an increase in Rayleigh waves (a type of seismic wave that can be generated by surface crack opening). However, we do not observe an increase in Rayleigh wave activity prior to or during brine release. When we estimate the average size of the fractures that we detect, we find the size is similar to many types of cracks in the nearby environment (for instance, cracks in the lake ice). We conclude that any Rayleigh wave seismicity that occurred during the Blood Falls brine release must be from a series of crack openings smaller than our detection limit, and that other cracks opening in the nearby environment may mask any signal specific to the Blood Falls release. Key Points: Time‐lapse photos capture wintertime Blood Falls brine releaseBrine release occurs without evidence for enhanced Rayleigh wave seismicity near the release pointThe Blood Falls crack may open as a series of small fracture events, masked by local seismicity [ABSTRACT FROM AUTHOR]

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

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

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

Russian (former Soviet) systematic studies of Antarctica by radio-echo sounding (RES) and ground-penetrating radar technique (GPR) were commenced in 1964. Since that time airborne RES surveys have covered about 5.5 × 106 km2 of the icy continent discovering remarkable geographic objects such as Subglacial Gamburtsev Mountains, and allowed studies of Filchner-Ronne Ice Shelf, Amery Ice Shelf and Lambert Glacier. Ground-based investigations during the 1990s and 2000s revealed the structure of the Lake Vostok area and surveyed along the Mirny to Vostok and Progress to Vostok traverse routes. GPR studies during the 2010s were to select the site for a new snow-runway at Mirny Station, with the resumption of the aviation after a 25 year hiatus.

Retreat of calving glaciers worldwide has contributed substantially to sea-level rise in recent decades. Mass loss by calving contributes significantly to the uncertainty of sea-level rise projections. At Bowdoin Glacier, Northwest Greenland, most calving occurs by a few large events resulting from kilometre-scale fractures forming parallel to the calving front. High-resolution terrestrial radar interferometry data of such an event reveal that crevasse opening is fastest at low tide and accelerates during the final 36 h before calving. Using the ice flow model Elmer/Ice, we identify the crevasse water level as a key driver of modelled opening rates. Sea water-level variations in the range of local tidal amplitude (1 m) can reproduce observed opening rate fluctuations, provided crevasse water level is at least 4 m above the low-tide sea level. The accelerated opening rates within the final 36 h before calving can be modelled by additional meltwater input into the crevasse, enhanced ice cliff undercutting by submarine melt, ice damage increase due to tidal cyclic fatigue, crevasse deepening or a combination of these processes. Our results highlight the influence of surface meltwater and tides on crevasse opening leading to major calving events at grounded tidewater glaciers such as Bowdoin.

Finely resolved geodetic data provide an opportunity to assess the extent and morphology of crevasses and their change over time. Crevasses have the potential to bias geodetic measurements of elevation and mass change unless they are properly accounted for. We developed a framework that automatically maps and extracts crevasse geometry and masks them where they interfere with surface mass-balance assessment. Our study examines airborne light detection and ranging digital elevation models (LiDAR DEMs) from Haig Glacier, which is experiencing a transient response in its crevassed upper regions as the glacier thins, using a self-organizing map algorithm. This method successfully extracts and characterizes ~1000 crevasses, with an overall accuracy of 94%. The resulting map provides insight into stress and flow conditions. The crevasse mask also enables refined geodetic estimates of summer mass balance. From differencing of September and April LiDAR DEMs, the raw LiDAR DEM gives a 9% overestimate in the magnitude of glacier thinning over the summer: −5.48 m compared with a mean elevation change of −5.02 m when crevasses are masked out. Without identification and removal of crevasses, the LiDAR-derived summer mass balance therefore has a negative bias relative to the glaciological surface mass balance.

While large-scale observations of intensified fracture and rifting can be observed through remote-sensing observations, understanding crevasse initiation may best be achieved with small-scale observations in which crevasses can be directly observed. Here we investigate the kinematic drivers of crevasse initiation in the McMurdo Shear Zone (MSZ), Antarctica. We delineated 420 crevasses from ~95 km of 400 MHz frequency ground-penetrating radar data and compared these data with kinematic outputs derived from remotely-sensed ice surface velocities to develop a statistical method to estimate crevasse initiation threshold strain rate values. We found the MSZ to be dominated by simple shear and that surface shear strain rates proved best for predicting crevasse features, with regions of higher shear strain rate more likely to have a greater number of crevasses. In the surveyed portion of our study region, values of shear strain rate and vorticity rate derived from the MEaSUREs2 velocity dataset range between 0.005–0.020 and 0.006–0.022 a−1, respectively, with crevasses located at ≥0.011 and ≥0.013 a−1. While threshold values from this study cannot be directly applied to other glacial environments, the method described here should allow for the study of shear margin evolution and assessment of localized damage and weakening processes in other locations where in situ data are available.

Surface crevasses on the Greenland Ice Sheet (GrIS) capture nearly half of the seasonal runoff, yet their role in transferring meltwater to the bed has received little attention relative to that of supraglacial lakes and moulins. Here, we present observations of crevasse ponding and investigate controls on their hydrological behavior at a fast‐moving, marine‐terminating sector of the GrIS. We map surface meltwater, crevasses, and surface‐parallel stress across a ∼2,700 km2 region using satellite data and contemporaneous uncrewed aerial vehicle (UAV) surveys. From 2017 to 2019 an average of 26% of the crevassed area exhibited ponding at locations that remained persistent between years despite rapid advection. We find that the spatial distribution of ponded crevasses does not relate to previously proposed controls on the distribution of supraglacial lakes (elevation and topography) or crevasses (von Mises stress thresholds), suggesting the operation of some other physical control(s). Ponded crevasse fields were preferentially located in regions of compressive surface‐parallel mean stress, which we interpret to result from the hydraulic isolation of these systems. This contrasts with unponded crevasse fields, which we suggest are readily able to transport meltwater into the wider supraglacial and englacial network. UAV observations show that ponded crevasses can drain episodically and rapidly, likely through hydrofracture. We therefore propose that the surface stress regime influences a spatially heterogeneous transfer of meltwater through crevasses to the bed of ice sheets, with consequences for processes, such as subglacial drainage and the heating of ice via latent heat release by refreezing meltwater. Plain Language Summary: Surface crevasses on the Greenland Ice Sheet (GrIS) transfer nearly half of all the ice sheet's meltwater to the bed, yet when, where, and how this occurs is poorly understood. We use large‐scale satellite analysis and uncrewed aerial vehicle (UAV) surveys to assess the spatial variability of meltwater ponding in crevasses across a fast‐flowing sector of the western GrIS. Between 2017 and 2019 an average of a quarter of crevasse fields ponded, in roughly the same area every year. However, the locations of such ponding cannot be explained in the same way as supraglacial lakes, which collect in surface basins. Instead, we find that ponded crevasses exist in regions of compressive surface stress. We suggest that this is because compressive regimes close pathways that elsewhere allow crevasses to drain into the wider surface hydrological system. Using UAV surveys, we show that these ponded crevasses instead drain rapidly to the bed by hydrofracture. Differing drainage processes in regions of compressive and extensional regimes may have distinct consequences for subglacial drainage and the heating of the ice sheet due to energy release during meltwater refreezing. Key Points: At a marine‐terminating sector of the Greenland Ice Sheet, 26% of all crevasses ponded between 2017 and 2019; the rest remained unfilledPonded crevasse locations persist between years but, unlike supraglacial lakes, this distribution is not topographically controlledCrevasses pond in compressive stress regimes, likely due to hydraulic isolation by creep closure, and drain predominantly via hydrofracture [ABSTRACT FROM AUTHOR]

Crary Ice Rise formed after the Ross Ice Shelf re-grounded ~1 kyr BP. We present new ice-penetrating radar data from two systems operating at center frequencies of 7 and 750 MHz that confirm the ice rise is composed of a former ice shelf buried by subsequent accumulation. Stacks of englacial diffraction hyperbolas are present almost everywhere across the central ice rise and extend up to ~350 m above the bed. In many cases, bed reflections beneath the diffraction hyperbolas are obscured for distances up to 1 km. Waveform modeling indicates that the diffraction hyperbolas are likely caused by marine ice deposits in former basal crevasses and rifts. The in-filling of rifts and basal crevasses may have strengthened the connection between the ice rise and the surrounding ice shelf, which could have influenced local and regional ice dynamics. Three internal reflection horizons mark the upper limit of disturbed ice and diffraction hyperbolas in different sections of the ice rise, indicating at least three stages of flow stabilization across the ice rise. A surface lineation visible in MODIS imagery corresponds spatially to deepening and strong deformation of these layers, consistent with the characteristics of former grounding lines observed elsewhere in Antarctica. [ABSTRACT FROM AUTHOR]

Seasonal lake Gornersee forms at the confluence of Gornergletscher and Grenzgletscher, Switzerland, and experiences outburst floods annually in midsummer. To study the interplay between lake drainage, glacier movement and crevasse activity, high-frequency seismometers and GPS receivers were deployed in networks near Gornersee during the summer ablation seasons of 2004, 2006 and 2007. We use a Rayleigh wave coherence method to locate 3289, 7939 and 4087 icequakes, respectively, primarily along well-defined surface crevasses. We calculate two-dimensional strains from triads of GPS stations and find mean differential strain rates of ~300 × 10−6 d−1 with diurnal variations up to 800 × 10−6 d−1. Crevasse icequake activity and glacial velocity are highest during early season, then decrease as meltwater channels erode and subglacial water pressure decreases. Glacial response to Gornersee drainage varied year-to-year, with icequake activity promoted at some crevasses and inhibited at others, suggesting syn-drainage icequakes may be indicative of local drainage patterns and small-scale features of the stress field. Diurnal pulses in icequake activity exhibit peak activity at different times of day in different locations, coincident with a southeast-to-northwest trending concentrated shear zone near the Gornergletscher–Grenzgletscher confluence, likely due to differences in the timing of peak strain rate in these regions.

Changing rates of water input can affect both the flow of glaciers and ice sheets and their propensity to crevasse. Here we examine geodetic and seismic observations during two substantial (10–18-times background velocity) rain-induced glacier accelerations at Haupapa/Tasman Glacier, New Zealand. Changes in rain rate result in glacier acceleration and associated uplift, which propagate down-glacier. This pattern of acceleration results in a change to the strain rate field, which correlates with an order of magnitude increase in the apparent seismicity rate and an overall down-glacier migration in located seismicity. After each acceleration event the apparent seismicity rate decreases to below the pre-acceleration rate for 3 days. This suggests that seismic events associated with surface crevasse growth occur early during phases of glacier acceleration due to elevated extensional stresses, and then do not occur again until stresses recover.

Icequakes at or near the bed of a glacier have the potential to allow us to investigate the interaction of ice with the underlying till or bedrock. Understanding this interaction is important for studying basal sliding of glaciers and ice streams, a critical process in ice dynamics models used to constrain future sea-level rise projections. However, seismic observations on glaciers can be dominated by seismic energy from surface crevassing. We present a method of automatically detecting basal icequakes and discriminating them from surface crevassing, comparing this method to a commonly used spectrum-based method of detecting icequakes. We use data from Skeidararjökull, an outlet glacier of the Vatnajökull Ice Cap, South-East Iceland, to demonstrate that our method outperforms the commonly used spectrum-based method. Our method detects a higher number of basal icequakes, has a lower rate of incorrectly identifying crevassing as basal icequakes and detects an additional, spatially independent basal icequake cluster. We also show independently that the icequakes do not originate from near the glacier surface. We conclude that the method described here is more effective than currently implemented methods for detecting and discriminating basal icequakes from surface crevassing.

Crevasses and englacial fracture networks route meltwater from a glacier's surface to the subglacial drainage system and thus influence glacial hydraulics. However, rapid fracture growth may also lead to sudden and potentially hazardous structural failure of unstable glaciers and ice dams, rifting of ice shelves, or iceberg calving. Here, we use passive seismic recordings to investigate the englacial fracture network on Glacier de la Plaine Morte, Switzerland. Glacier dynamics and the drainage of an ice-marginal lake give rise to numerous icequakes, the majority of which generate dispersed, high-frequency Rayleigh waves. A wide distribution of events allows us to study azimuthal anisotropy between 10 and 30 Hz in order to extract englacial seismic velocities in regions of preferentially oriented crevasses. Beamforming applied to a 100-m-aperture array reveals azimuthal anisotropy of Rayleigh-wave phase velocities reaching a strength of 8% at high frequencies. In addition, we find that the fast direction of wave propagation coincides with the observed surface strike of the narrow crevasses. Forward modeling and inversion of dispersion curves suggest that the azimuthal anisotropy is induced by a 40-m-thick crevassed layer at the surface of the glacier with 8% anisotropy in shear-wave velocity.

Passive seismology allows measurement of the structure of glaciers and ice sheets. However, most techniques used so far in this context are based on horizontally homogeneous media where parameters vary only with depth (1-D approximations), which are appropriate only for a subset of glaciers. Here, we analyze seismic noise records from three different types of glaciers (plateau, valley and avalanching glacier) to characterize the influence of the glacier geometry on the seismic wavefield. Using horizontal-to-vertical spectral ratios, polarization analysis and modal analysis, we show that the plateau glacier and the valley glacier can be seen as 1-D, whereas the relatively small avalanching glacier shows 3-D effects due to its bed topography and the deep crevasses. In principle, the techniques proposed here might allow monitoring such crevasses and their depth, and thus to constrain a key parameter of avalanching and calving glacier fronts.

GPS measurements of tidal modulation of ice flow and seismicity within the grounding zone of Beardmore Glacier show that tidally induced fluctuations of horizontal flow are largest near the grounding line and decrease downstream. Seismic activity is continuous, but peaks occur on falling and rising tides. Beamforming methods reveal that most seismic events originate from two distinct locations, one on the grid-north side of the grounding zone, and one on the grid-south side. The broad pattern of deformation generated as Beardmore Glacier merges with the Ross Ice Shelf results in net extension along the grid-north side of the grounding zone and net compression along the grid-south side. During falling tides, seismic activity peaks on both sides because of increased vertical flexure across the grounding line. During rising tides, seismic activity in the region of extension on the grid-north side is relatively low because the tidal influence on both horizontal strain rate and vertical flexure is small. On the grid-south side during rising tides, however, tidally induced horizontal strain rates promote increased seismicity in regions of long-term compressional flow paths. Our study highlights how concurrent geodetic and seismic measurements provide insight into grounding-zone mechanics and their influence on ice-shelf buttressing.

Glacial seismicity provides important insights into glacier dynamic processes. We study the temporal distribution of cryogenic seismic signals (icequakes) at Holtedahlfonna, Svalbard, between April and August 2016 using a single three-component sensor. We investigate sources of observed icequakes using polarization analysis and waveform modeling. Processes responsible for five icequake categories are suggested, incorporating observations of previous studies into our interpretation. We infer that the most dominant icequake type is generated by surface crevasse opening through hydrofracturing. Secondly, bursts of high-frequency signals are presumably caused by repeated near-surface crevassing due to high strain rates during glacier fast-flow episodes. Furthermore, signals related to resonance in water-filled cracks, fracturing or settling events in dry firn or snow before the melt season, and processes at the glacier bed are observed. Amplitude of seismic background noise is clearly related to glacier runoff. We process ambient seismic noise to invert horizontal-to-vertical spectral ratios for a sub-surface seismic velocity model used to model icequake signals. Our study shows that a single seismic sensor provides useful information about seasonal ice dynamics in case deployment of a network is not feasible.

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

Natural levees are common features in river, delta and tidal landscapes. They are elevated near‐channel morphological features that determine the connection between channel and floodbasin, and consequently affect long‐term evolution up to delta‐scales. Despite their relevance in shaping fluvial‐tidal systems, research on levees is sparse and often limited to fluvial or non‐tidal case studies. There is also a general lack of understanding of the role of vegetation in shaping these geomorphic units, and how levee morphology and dimensions vary in the transition from fluvial to coastal environments, where tides are increasingly important. Our goal is to unravel the effects of fluvial‐tidal boundary conditions, sediment supply and vegetation on levee characteristics and floodbasin evolution. These conditions were systematically explored by 60 large‐scale idealized morphodynamic simulations in Delft3D which self‐developed levees over the course of one century. We compared our results to a global levee dataset compilation of natural levee dimensions. We found that levee height is determined by the maximum water level, provided sufficient levee building sediments are available. Discharge fluctuations increased levee width and triggered more levee breaches, i.e. crevasses, that effectively filled the fluvio‐tidal floodbasin. The presence of wood‐type (sparse) vegetation further increased the number of crevasses in comparison with the non‐vegetated scenarios. Conversely, reed‐type (dense) vegetation strongly dampened tidal amplitude and reduced the accommodation space and sedimentation further into the floodbasin, resulting in narrower levees, no crevasses and limited floodbasin accretion. However, dense vegetation reduced tidal forces which allowed levee growth further downstream. Ultimately, the levees merged with the coastal barrier, eliminating the floodbasin tides entirely. Our results elucidate the mechanisms by which levee and crevasse formation, and vegetation may fill fluvio‐tidal wetlands and affect estuary evolution. This brings new insights for geological reconstructions as well as for the future management of deltas and estuaries under sea‐level rise. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd [ABSTRACT FROM AUTHOR]

Projections of future ice sheet mass loss and thus sea level rise rely on the parametrization of iceberg calving in ice sheet models. The interconnection between submarine melt-induced undercutting and calving is still poorly understood, which makes predicted contributions of tidewater glaciers to sea level rise uncertain. Here, we compare detailed 3-D simulations of fracture initiation obtained with the Helsinki Discrete Element Model (HiDEM) to observations, prior to a major calving event at Bowdoin Glacier, Northwest Greenland. Observations of a plume surfacing at the calving location suggest that local melt-undercutting influenced the size of the major calving event. Therefore, several experiments are conducted with various local and distributed (front-wide) undercut geometries. Although the number of undercut experiments is limited by computational requirements, one of the conjectured undercut geometries reproduces the crevasse leading to the observed major calving event in great detail. Our simulations show that undercutting leads to initiation of wider fractures more than 100 m upstream of the terminus, well-beyond the directly undercut region. When combining a moderate distributed undercut with local amplified undercuts at the two observed plumes, fracture initiation also increases in between the local undercuts. Thus, our results agree with previous studies suggesting the existence of a “calving amplifier” effect by submarine melt, both upglacier and across-glacier. Consequently, the simulations show the potentially large impact of submarine melt-induced undercutting on iceberg size.

Crevasses are formed by glacier movement and the stresses within glacier ice. Knowledge of the crevasses’ distribution is critical for understanding the glacier and ice shelf stability. In this study, we propose an automated crevasse extraction framework based on Sentinel-1 SAR imagery and an improved U-Net network. The spatial distribution of crevasses on Antarctic ice shelves in 2020 was mapped with a spatial resolution of ~40 m, and the characteristics of crevasses on the Nickerson Ice Shelf, Jelbart Ice Shelf, Amery Ice Shelf, Thwaites Glacier, and Shackleton Ice Shelf were analyzed. The results indicated the extraction accuracy of our method was 84.2% and the F1 score was 72.5%. Compared with previous published studies, the identification of the crevasse areas had good visual consistency. However, in some scenes, the recall rate was relatively lower due to the quality of the SAR image, terrain surrounding the crevasses, and observation geometry. The crevasses on different ice shelves had different characteristics in terms of length, density, type, and spatial pattern, implying the different stress structures of ice shelves. The Thwaites Glacier and the Nickerson Ice Shelf in the West Antarctica Ice Sheet (WAIS) had shorter ice crevasses, whereas the lengths of ice crevasses on the Jelbart Ice Shelf and the Amery Ice Shelf in the East Antarctica Ice Sheet (EAIS) were relatively long. Nevertheless, there are more closely spaced crevasses on the ice shelf in WAIS compared to that in the EAIS. For the distribution of crevasse types, the Nickerson Ice Shelf and the Shackleton Ice Shelf had various forms of crevasses. There were mainly transverse crevasses developed on the Jelbart Ice Shelf and the Amery Ice Shelf. This study provides a helpful reference and guidance for automated crevasse extraction. The method proposed by this study manifests great application potential and the efficacy of producing a time-series crevasse data set with higher spatial resolution and larger coverage. In the future, more Sentinel-1 SAR imagery will be applied and the effect of temporal and spatial variations in crevasses on the stability of ice shelves will be investigated, which will contribute to project the ice shelf stability and explore the sea level rise implications of recent and future cryosphere changes.

We present a 1:42,000 scale map of Isunguata Sermia, a land-terminating outlet glacier draining the western-sector of the Greenland Ice Sheet. Structure-from-Motion software applied to ∼3,600 aerial images collected by a fixed-wing unmanned aerial vehicle in July 2015 allowed us to produce a high resolution (0.3 m ground sampling distance (GSD)) orthomosaic and digital elevation model (DEM; 1.5 m GSD).These products were used to map and describe the structural, geomorphological and hydrological features of the lower 16 km terminus of Isunguata Sermia and include many thousands of crevasses, crevasse traces and supraglacial channels. Additionally, several geomorphological features and pro-glacial hydrological features were identified, including debris-covered ice, lateral moraines and ice-marginal lakes. The map has potential for informing and reconstructing the long-term dynamic history of the glacier, including its response to variable environmental forcing.

Aoraki Mount Cook is New Zealand's highest mountain and a popular destination for climbers. This study combined historical accounts, photos, and geophysical surveys with modern-day spatial analysis, popular narratives, and interviews to explore how conditions on the most popular climbing route up the mountain, the Linda Glacier route, have changed over time. Results highlight significant change on the lower section of the route due to ongoing downwasting of the Tasman Glacier; but higher on the mountain, changes in the route are more strongly associated with year-to-year variability in snow conditions. Even so, recent observations of new rock exposures and a shortening of the climbing season due to earlier crevasse exposure may be indications of longer-term glacier change.

We investigate the appropriateness of calving or crevasse models from the literature using linear elastic fracture mechanics (LEFM). To this end, we compare LEFM model-predicted stress intensity factors (SIFs) against numerically computed SIFs using the displacement correlation method in conjunction with the finite element method. We present several benchmark simulations wherein we calculate the SIF at the tips of water-filled surface and basal crevasses penetrating through rectangular ice slabs under different boundary conditions, including grounded and floating conditions. Our simulation results indicate that the basal boundary condition significantly influences the SIF at the crevasse tips. We find that the existing calving models using LEFM are not generally accurate for evaluating SIFs in grounded glaciers or floating ice shelves. We also illustrate that using the ‘single edge crack’ weight function in the LEFM formulations may be appropriate for predicting calving from floating ice shelves, owing to the low fracture toughness of ice; whereas, using the ‘double edge crack’ or ‘central through crack’ weight functions is more appropriate for predicting calving from grounded glaciers. To conclude, we recommend using the displacement correlation method for SIF evaluation in real glaciers and ice shelves with complex geometries and boundary conditions.

Over recent decades, Greenland ice sheet surface melt has shown an increase both in intensity and spatial extent. Part of this water probably reaches the bed and can enhance glacier speed, advecting a larger volume of ice into the ablation area. In the context of a warming climate, this mechanism could contribute to the future rate of thinning and retreat of land-terminating glaciers of Greenland. These changes in ice flow conditions will in turn influence surface crevassing and thus the ability of water to reach the bed at higher elevations. Here, using a coupled basal hydrology and prognostic ice flow model, the evolution of a Greenland-type glacier subject to increasing surface melt is studied over a few decades. For different scenarios of surface melt increase over the next decades, the evolution of crevassed areas and the ability of water to reach the bed is inferred. Our results indicate that the currently observed crevasse distribution is likely to extend further upstream which will allow water to reach the bed at higher elevations. This will lead to an increase in ice flux into the ablation area which, in turn, accelerates the mass loss of land-terminating glaciers.

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