Heterogeneous melting near the Thwaites Glacier grounding line

Thwaites Glacier represents 15% of the ice discharge from the West Antarctic Ice Sheet and influences a wider catchment^1-3. Because it is grounded below sea level^4,5, Thwaites Glacier is thought to be susceptible to runaway retreat triggered at the grounding line (GL) at which the glacier reaches the ocean^6,7. Recent ice-flow acceleration^2,8 and retreat of the ice front^8-10 and GL^11,12 indicate that ice loss will continue. The relative impacts of mechanisms underlying recent retreat are however uncertain. Here we show sustained GL retreat from at least 2011 to 2020 and resolve mechanisms of ice-shelf melt at the submetre scale. Our conclusions are based on observations of the Thwaites Eastern Ice Shelf (TEIS) from an underwater vehicle, extending from the GL to 3 km oceanward and from the ice-ocean interface to the sea floor. These observations show a rough ice base above a sea floor sloping upward towards the GL and an ocean cavity in which the warmest water exceeds 2 °C above freezing. Data closest to the ice base show that enhanced melting occurs along sloped surfaces that initiate near the GL and evolve into steep-sided terraces. This pronounced melting along steep ice faces, including in crevasses, produces stratification that suppresses melt along flat interfaces. These data imply that slope-dependent melting sculpts the ice base and acts as an important response to ocean warming.

Heterogeneous retreat and ice melt of Thwaites Glacier, West Antarctica

The glaciers flowing into the Amundsen Sea Embayment, West Antarctica, have undergone acceleration and grounding line retreat over the past few decades that may yield an irreversible mass loss. Using a constellation of satellites, we detect the evolution of ice velocity, ice thinning, and grounding line retreat of Thwaites Glacier from 1992 to 2017. The results reveal a complex pattern of retreat and ice melt, with sectors retreating at 0.8 km/year and floating ice melting at 200 m/year, while others retreat at 0.3 km/year with ice melting 10 times slower. We interpret the results in terms of buoyancy/slope-driven seawater intrusion along preferential channels at tidal frequencies leading to more efficient melt in newly formed cavities. Such complexities in ice-ocean interaction are not currently represented in coupled ice sheet/ocean models.

A Century of Stability of Avannarleq and Kujalleq Glaciers, West Greenland, Explained Using High-Resolution Airborne Gravity and Other Data

The evolution of Greenland glaciers in a warming climate depends on their depth below sea level, flow speed, surface melt, and ocean-induced undercutting at the calving front. We present an innovative mapping of bed topography in the frontal regions of Sermeq Avannarleq and Kujalleq, two major glaciers flowing into the ice-choked Torssukatak Fjord, central west Greenland. The mapping combines a mass conservation algorithm inland, multibeam echo sounding data in the fjord, and high-resolution airborne gravity data at the ice-ocean transition where other approaches have traditionally failed. We obtain a reliable, precision (±40 m) solution for bed topography across the ice-ocean boundary. The results reveal a 700 m deep fjord that abruptly ends on a 100-300 m deep sill along the calving fronts. The shallow sills explain the presence of stranded icebergs, the resilience of the glaciers to ocean-induced undercutting by warm Atlantic water, and their remarkable stability over the past century.

Vulnerability of Southeast Greenland Glaciers to Warm Atlantic Water From Operation IceBridge and Ocean Melting Greenland Data

We employ National Aeronautics and Space Administration (NASA)'s Operation IceBridge high-resolution airborne gravity from 2016, NASA's Ocean Melting Greenland bathymetry from 2015, ice thickness from Operation IceBridge from 2010 to 2015, and BedMachine v3 to analyze 20 major southeast Greenland glaciers. The results reveal glacial fjords several hundreds of meters deeper than previously thought; the full extent of the marine-based portions of the glaciers; deep troughs enabling warm, salty Atlantic Water (AW) to reach the glacier fronts and melt them from below; and few shallow sills that limit the access of AW. The new oceanographic and topographic data help to fully resolve the complex pattern of historical ice front positions from the 1930s to 2017: glaciers exposed to AW and resting on retrograde beds have retreated rapidly, while glaciers perched on shallow sills or standing in colder waters or with major sills in the fjords have remained stable.

BedMachine v3: Complete Bed Topography and Ocean Bathymetry Mapping of Greenland From Multibeam Echo Sounding Combined With Mass Conservation

Greenland's bed topography is a primary control on ice flow, grounding line migration, calving dynamics, and subglacial drainage. Moreover, fjord bathymetry regulates the penetration of warm Atlantic water (AW) that rapidly melts and undercuts Greenland's marine-terminating glaciers. Here we present a new compilation of Greenland bed topography that assimilates seafloor bathymetry and ice thickness data through a mass conservation approach. A new 150 m horizontal resolution bed topography/bathymetric map of Greenland is constructed with seamless transitions at the ice/ocean interface, yielding major improvements over previous data sets, particularly in the marine-terminating sectors of northwest and southeast Greenland. Our map reveals that the total sea level potential of the Greenland ice sheet is 7.42 ± 0.05 m, which is 7 cm greater than previous estimates. Furthermore, it explains recent calving front response of numerous outlet glaciers and reveals new pathways by which AW can access glaciers with marine-based basins, thereby highlighting sectors of Greenland that are most vulnerable to future oceanic forcing.

Fast retreat of Zachariæ Isstrøm, northeast Greenland

After 8 years of decay of its ice shelf, Zachariæ Isstrøm, a major glacier of northeast Greenland that holds a 0.5-meter sea-level rise equivalent, entered a phase of accelerated retreat in fall 2012. The acceleration rate of its ice velocity tripled, melting of its residual ice shelf and thinning of its grounded portion doubled, and calving is now occurring at its grounding line. Warmer air and ocean temperatures have caused the glacier to detach from a stabilizing sill and retreat rapidly along a downward-sloping, marine-based bed. Its equal-ice-volume neighbor, Nioghalvfjerdsfjorden, is also melting rapidly but retreating slowly along an upward-sloping bed. The destabilization of this marine-based sector will increase sea-level rise from the Greenland Ice Sheet for decades to come.

Observed thinning of Totten Glacier is linked to coastal polynya variability

Analysis of ICESat-1 data (2003-2008) shows significant surface lowering of Totten Glacier, the glacier discharging the largest volume of ice in East Antarctica, and less change on nearby Moscow University Glacier. After accounting for firn compaction anomalies, the thinning appears to coincide with fast-flowing ice indicating a dynamical origin. Here, to elucidate these observations, we apply high-resolution ice-ocean modelling. Totten Ice Shelf is simulated to have higher, more variable basal melting rates. We link this variability to the volume of cold water, originating in polynyas upon sea ice formation, reaching the sub-ice-shelf cavity. Hence, we propose that the observed increased thinning of Totten Glacier is due to enhanced basal melting caused by a decrease in cold polynya water reaching its cavity. We support this hypothesis with passive microwave data of polynya extent variability. Considering the widespread changes in sea ice conditions, this mechanism could be contributing extensively to ice-shelf instability.

Ice-shelf melting around Antarctica

We compare the volume flux divergence of Antarctic ice shelves in 2007 and 2008 with 1979 to 2010 surface accumulation and 2003 to 2008 thinning to determine their rates of melting and mass balance. Basal melt of 1325 ± 235 gigatons per year (Gt/year) exceeds a calving flux of 1089 ± 139 Gt/year, making ice-shelf melting the largest ablation process in Antarctica. The giant cold-cavity Ross, Filchner, and Ronne ice shelves covering two-thirds of the total ice-shelf area account for only 15% of net melting. Half of the meltwater comes from 10 small, warm-cavity Southeast Pacific ice shelves occupying 8% of the area. A similar high melt/area ratio is found for six East Antarctic ice shelves, implying undocumented strong ocean thermal forcing on their deep grounding lines.

Accelerated sea-level rise from West Antarctica

Recent aircraft and satellite laser altimeter surveys of the Amundsen Sea sector of West Antarctica show that local glaciers are discharging about 250 cubic kilometers of ice per year to the ocean, almost 60% more than is accumulated within their catchment basins. This discharge is sufficient to raise sea level by more than 0.2 millimeters per year. Glacier thinning rates near the coast during 2002-2003 are much larger than those observed during the 1990s. Most of these glaciers flow into floating ice shelves over bedrock up to hundreds of meters deeper than previous estimates, providing exit routes for ice from further inland if ice-sheet collapse is under way.

Segmentation of polarimetric synthetic aperture radar data

A statistical image model is proposed for segmenting polarimetric synthetic aperture radar (SAR) data into regions of homogeneous and similar polarimetric backscatter characteristics. A model for the conditional distribution of the polarimetric complex data is combined with a Markov random field representation for the distribution of the region labels to obtain the posterior distribution. Optimal region labeling of the data is then defined as maximizing the posterior distribution of the region labels given the polarimetric SAR complex data (maximum a posteriori (MAP) estimate). Two procedures for selecting the characteristics of the regions are then discussed. Results using real multilook polarimetric SAR complex data are given to illustrate the potential of the two selection procedures and evaluate the performance of the MAP segmentation technique. It is also shown that dual polarization SAR data can yield segmentation resultS similar to those obtained with fully polarimetric SAR data.