Rift propagation, rather than basal melt, drives the destabilization and disintegration of the Thwaites Eastern Ice Shelf. Since 2016, rifts have episodically advanced throughout the central ice-shelf area, with rapid propagation events occurring during austral spring. The ice shelf's speed has increased by ~70% during this period, transitioning from a rate of 1.65 m d−1 in 2019 to 2.85 m d−1 by early 2023 in the central area. The increase in longitudinal strain rates near the grounding zone has led to full-thickness rifts and melange-filled gaps since 2020. A recent sea-ice break out has accelerated retreat at the western calving front, effectively separating the ice shelf from what remained of its northwestern pinning point. Meanwhile, a distributed set of phase-sensitive radar measurements indicates that the basal melting rate is generally small, likely due to a widespread robust ocean stratification beneath the ice–ocean interface that suppresses basal melt despite the presence of substantial oceanic heat at depth. These observations in combination with damage modeling show that, while ocean forcing is responsible for triggering the current West Antarctic ice retreat, the Thwaites Eastern Ice Shelf is experiencing dynamic feedbacks over decadal timescales that are driving ice-shelf disintegration, now independent of basal melt.

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

To increase inclusivity, diversity, equity and accessibility in Antarctic science, we must build more positive and inclusive Antarctic field work environments. The International Thwaites Glacier Collaboration (ITGC) has engaged in efforts to contribute to that goal through a variety of activities since 2018, including creating an open-access ‘Field and Ship Best Practices’ guide, engaging in pre-field season team dynamics meetings, and surveying post-field season reflections and experiences. We report specific actions taken by ITGC and their outcomes. We found that strong and supported early career researchers brought new and important perspectives regarding strategies for transforming culture. We discovered that engaged and involved senior leadership was also critical for expanding participation and securing funding to support efforts. Pre-field discussions involving all field team members were particularly helpful for setting expectations, improving sense of belonging, describing field work best practices, and co-creating a positive work culture.

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

Basal channels, which form where buoyant plumes of ocean water and meltwater carve troughs upwards into ice-shelf bases, are widespread on Antarctic ice shelves. The formation of these features modulates ice-shelf basal melt by influencing the flow of buoyant plumes, and influences structural stability through concentration of strain and interactions with fractures. Because of these effects, and because basal channels can change rapidly, on timescales similar to those of ice-shelf evolution, constraining the impacts of basal channels on ice shelves is necessary for predicting future ice-shelf destabilization and retreat. We suggest that future research priorities should include constraining patterns and rates of basal channel change, determining mechanisms and detailed patterns of basal melt, and quantifying the influence that channel-related fractures have on ice-shelf stability.

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

Strain rates are fundamental measures of ice flow and are used in a wide variety of glaciological applications including investigations of bed properties, calculations of basal mass balance on ice shelves, and constraints on ice rheological models. However, despite their extensive application, strain rates are calculated using a variety of methods and length scales and the details are often not specified. In this study, we compare the results of nominal and logarithmic strain-rate calculations based on a satellite-derived velocity field of the Antarctic ice sheet generated from Landsat 8 satellite data. Our comparison highlights the differences between the two common approaches in the glaciological literature. We evaluate the errors introduced by each approach and their impacts on the results. We also demonstrate the importance of choosing and specifying a length scale over which strain-rate calculations are made, which can strongly influence other derived quantities such as basal mass balance on ice shelves. Finally, we present strain-rate data products calculated using an approximate viscous length-scale with satellite observations of ice velocity for the Antarctic continent.

The surface of the Ross Ice Shelf (RIS) is textured by flow stripes, crevasses and other features related to ice flow and deformation. Here, moderate resolution optical satellite images are used to map and classify regions of the RIS characterized by different surface textures. Because the textures arise from ice deformation, the map is used to identify structural provinces with common deformation history. We classify four province types: regions associated with large outlet glaciers, shear zones, extension downstream of obstacles and suture zones between provinces with different upstream sources. Adjacent provinces with contrasting histories are in some locations deforming at different rates, suggesting that our province map is also an ice fabric map. Structural provinces have more complicated shapes in the part of the ice shelf fed by West Antarctic ice streams than in the part fed by outlet glaciers from the Transantarctic Mountains. The map may be used to infer past variations in stress conditions and flow events that cannot be inferred from flow traces alone.

We have assembled an elevation grid for the Greenland ice sheet using a combination of the best current digital elevation model (DEM) (Bamber and others, 2000a, 2001) and 44 Advanced Very High Resolution Radiometer satellite images acquired in spring 1997. The images are used to quantitatively enhance the representation of surface undulations through photoclinometry. Gridcell spacing of the new DEM is 625 m. To validate the new DEM, we compared profiles extracted from it and the Bamber and others DEM with airborne laser altimetry profiles collected in the 1990s by the Airborne Topographic Mapper (Krabill and others, 1995). The image-enhanced DEM has a greatly improved representation of decameter-relief surface features <15 km in lateral extent, and reduces the mean elevation error in regions having these features by 20–50%. Root-mean-squared errors are typically 7–15m in the Bamber DEM, and 4–10m after image enhancement. However, the photoclinometry process adds some noise. In very smooth portions of the ice sheet where decameter undulationsare absent, the photoclinometry process caused a slight increase in the rms error, from ~1 min the Bamber and others DEM to ∼2.5 min the image-enhanced DEM. The image-enhanced DEM will be useful for inferring accumulation-rate variations over the undulation field, or for improving maps of bedrock elevation through inversion of surface elevation, for example. We briefly explore the preliminary steps of this latter application.

Ship-borne and airborne infrared radiometric measurements during the Arise cruise of September–October 2003 permitted in Situ validation Studies of two Satellite-based ice Surface Skin temperature algorithms: the AVHRR Polar Pathfinder Ice Surface Temperature and the MODIS Sea Ice Surface Temperature. Observations of Sea ice from the Aurora Australis Ship’s rail using a KT-19.82 radiometer were conducted between 25 September and 21 October during clear-sky overflights by AVHRR (41 passes) and MODIS (17 passes) on their respective Satellite platforms. Data from both Sensors Show highly linear fits to 1 min integrated radiometer Spot measurements, Spanning the range 245–270 K with a ±1.4˚C, 1σ (AVHRR) and ±1.0˚C (MODIS) variation relative to a 1: 1 relationship. There was no Significant offset. Helicopter observations made with a KT-19.85 radiometer on three dates (8, 19 and 20 October) provided more data (236 gridcell Sites total), but over a more limited Sea-ice Skin temperature range (252–268 K), with higher variation (±1.7˚C, 1σ) due to mixed-pixel issues. Comparison of MODIS and AVHRR algorithms directly, with both images acquired during a helicopter flight, indicates very high correlation and near-unity Slope for the two Satellite-based algorithms. Ship air-temperature data during the validation indicated moderate to Strong inversions over Sea ice under clear Skies. These formed and decayed rapidly (tens of minutes) as clouds moved out of and into the zenith area.

We use laser altimetry from the Ice, Cloud, and land Elevation Satellite (ICESat) to map the grounding zone (GZ) of the Ross Ice Shelf, Antarctica, at 491 locations where ICESat tracks cross the grounding line (GL). Ice flexure in the GZ occurs as the ice shelf responds to short-term sea-level changes due primarily to tides. ICESat repeat-track analysis can be used to detect this region of flexure since each repeated pass is acquired at a different tidal phase; the technique provides estimates for both the landward limit of flexure and the point where the ice becomes hydrostatically balanced. We find that the ICESat-derived landward limits of tidal flexure are, in many places, offset by several km (and up to ∼60km) from the GL mapped previously using other satellite methods. We discuss the reasons why different mapping methods lead to different GL estimates, including: instrument limitations; variability in the surface topographic structure of the GZ; and the presence of ice plains. We conclude that reliable and accurate mapping of the GL is most likely to be achieved when based on synthesis of several satellite datasets.

Ice dolines on the Larsen Ice Shelf, Antarctica, are observed to be elongated depressions a few hundred meters across and up to 19 m deep. One-meter resolution imagery is used to quantify these dimensions. Elevation profiles across five dolines are derived by photoclinometry. Landsat and radar imagery is also used to show that dolines can form in a single melt season and persist for years. Dolines occur in clusters and in direct proximity to surface meltwater lakes. Field observations suggest dolines form by collapse into a subsurface cavity. A direct hydraulic connection with the underlying ocean is believed necessary to drain water that would otherwise collect in dolines. A formation hypothesis is discussed consistent with these observations and with energy-and hydrostatic-imbalance considerations.

We identify a series of basal crevasses along a 31 km transect across the northern sector of the Larsen C ice shelf, Antarctica, using in situ ground-penetrating radar. The basal crevasses propagate from a region of multiple, shallow basal fractures to form widely spaced (0.5–2.0 km) but deeply incised (70–134 m) features. Surface troughs, observed in visible imagery, exist above the basal crevasses as the ice vertically shears to reach hydrostatic equilibrium, while widespread surface crevassing occurs along the crests and on the flanks of the undulations, primarily aligned with the topography. We suggest, based on the location of the surface crevasses and the along-flow evolution of the basal crevasses, that the former are induced by a bending stress created by gradients in hydrostatic forces. Using a linear elastic fracture mechanics model, we investigate the sensitivity of basal crevasse propagation to observed trends of ice-shelf thinning and acceleration. Basal crevasses are large-scale structural weaknesses that can both control meltwater ponding and induce surface crevassing. Together, these features may represent an important mechanism in both past and future ice-shelf disintegration events on the Antarctic Peninsula.

Ice surface altimetry from ICESat-1 and NASA aircraft altimeter overflights spanning 2002–09 indicate that a region of lower Crane Glacier, Antarctic Peninsula, shows an unusual temporal pattern of elevation loss: a period of very rapid drawdown (~91ma–1 between September 2004 and September 2005) bounded by periods of large but more moderate rates (23ma–1 until September 2004; 12ma–1 after September 2005). The region of increased drawdown is ~4.5 km ×2.2 km based on satellite (ASTER and SPOT-5) stereo-image digital elevation model (DEM) differencing spanning the event. In a later differential DEM the anomalous drawdown feature is not seen. Bathymetry in Crane Glacier fjord reveals a series of flat-lying, formerly subglacial deeps interpreted as lake sediment basins. We conclude that the elevation-change feature resulted from drainage of a small, deep subglacial lake. We infer that the drainage event was induced by hydraulic forcing of subglacial water past a downstream obstruction. However, only a fraction of Crane Glacier’s increase in flow speed that occurred near the time of lake drainage (derived from image feature tracking) appears to be directly attributable to the event; instead, retreat of the ice front off a subglacial ridge 6 km downstream of the lake is likely the dominant cause of renewed fast flow and more negative mass balance in the subsequent 4 years.

Snow and firn properties control the transport of vapor, gases and water between the atmosphere and the underlying strata. An understanding of this transport and the properties that control it is important for predicting air–snow transfer of chemical species and for interpreting ice cores. Remote-sensing images of East Antarctica show large areas of alternating light and dark bands. These low-amplitude, long-wavelength features have glazed downwind faces and rough upwind faces and are called megadunes. The first linked measurements of the permeability and the associated microstructure for a glazed area within a well-defined megadune area are reported in this paper. Permeability and density were measured, along with grain-scale properties derived from digital image processing of preserved thick sections, at this cold, low-accumulation glazed site. A clear layering pattern exists. In the top meter the firn density ranges from 0.24 to 0.50 g cm–3. Permeability measurements range from 50 x 10–10 to 200 x 10–10μ2, several times greater than corresponding profiles from warmer, higher-accumulation sites like Siple Dome, Antarctica. It is shown that buoyancy-driven natural convection may be important in post-depositional processes in very cold, low-accumulation sites like this.

We use satellite images to track seasonal and interannual variations in blue-ice extent over the past 30 years near Byrd Glacier on the East Antarctic plateau. The study areas have low slope and few nearby nunataks, which may increase their climate sensitivity. A threshold-based algorithm sensitive to snow grain-size is used to analyze 56 Moderate Resolution Imaging Spectroradiometer (MODIS) images over three recent summer seasons. Seasonal blue-ice exposure grows rapidly in late spring, and peaks by late December. Exposure is relatively constant between late December and mid-January, then declines in February. We interpret this cycle as due to removal and re-accumulation of patchy snow. Interannual changes in blue-ice area may be estimated by tracking the near-constant summer maximum extent period. Fifteen mid-summer Landsat images, spanning 1974–2002, were analyzed to determine long-term variations. Interannual area changes are 10–30%; however, the MODIS data revealed that the exposed blue-ice area can be sharply reduced for up to 2 weeks after a snowfall event; and in the 2001/02 season, patchy snow cover persisted for the entire summer. The combination of MODIS seasonal and Landsat interannual data indicates that blue-ice areas can be climate-sensitive. The strong feedback between snow cover and surface energy balance implies that blue-ice areas could rapidly decrease due to climate-related increases in snowfall or reduced ablation.

Automatic feature tracking on two Landsat images (acquired inJanuary 2000 and December 2001) generates a complete and accurate velocity field of Mertz Glacier, East Antarctica. This velocity field shows two main tributaries to the ice stream. Between the tributaries, a likely obstruction feature in the bedrock results in a slow-down of the flow. A third Landsat image, acquired in 1989 and combined with the 2000 image, permits the determination of the glacier mean velocity during the 1990s. Although some parts of the Mertz Glacier system show evidence of slight speed increase, we conclude that the Mertz flow speed is constant within our uncertainty (35 m a−1). Using this complete velocity field, new estimates of the ice discharge flux, 17.8 km3 a−1 (16.4 Gt a−1), and of the basal melting of the tongue, 11 m a−1 of ice, are given. Our results lead to an apparent imbalance of the drainage basin (ice discharge 3.5 km3 a−1 lower than the accumulation). Considering previous studies in the Mertz Glacier area, we then discuss the uncertainty of this imbalance and the problems with accumulation mapping for this region.

We examine patterns of localized surface elevation change in lower Mercer and Whillans Ice Streams, West Antarctica, which we interpret as subglacial water movement through a system of lakes and channels. We detect and measure the lake activity using repeat-track laser altimetry from ICESat and image differencing from MODIS image pairs. A hydrostatic-potential map for the region shows that the lakes are distributed across three distinct hydrologic regimes. Our analysis shows that, within these regimes, some of the subglacial lakes appear to be linked, with drainage events in one reservoir causing filling and follow-on drainage in adjacent lakes. We also observe changes near ice raft ‘a’ in lower Whillans Ice Stream, and interpret them as evidence of subglacial water and other changes at the bed. The study provides quantitative information about the properties of this complex subglacial hydrologic system, and a relatively unstudied component of ice-sheet mass balance: subglacial drainage across the grounding line.

We use models constrained by remotely sensed data from Pine Island and Thwaites Glaciers, West Antarctica, to infer basal properties that are difficult to observe directly. The results indicate strong basal melting in areas upstream of the grounding lines of both glaciers, where the ice flow is fast and the basal shear stress is large. Farther inland, we find that both glaciers have ‘mixed’ bed conditions, with extensive areas of both bedrock and weak till. In particular, there are weak areas along much of Pine Island Glacier’s main trunk that could prove unstable if it retreats past the band of strong bed just above its current grounding line. In agreement with earlier studies, our forward ice-stream model shows a strong sensitivity to small perturbations in the grounding line position. These results also reveal a large sensitivity to the assumed bed (sliding or deforming) model, with non-linear sliding laws producing substantially greater dynamic response than earlier simulations that assume a linear-viscous till rheology. Finally, comparison indicates that our results using a plastic bed are compatible with the limited observational constraints and theoretical work that suggests an upper bound exists on maximum basal shear stress.

Using RADARSAT synthetic aperture radar data, we have mapped the flow velocity over much of the Greenland ice sheet for the winters of 2000/01 and 2005/06. These maps provide a detailed view of the ice-sheet flow, including that of the hundreds of glaciers draining the interior. The focused patterns of flow at the coast suggest a strong influence of bedrock topography. Differences between our two maps confirm numerous early observations of accelerated outlet glacier flow as well as revealing previously unrecognized changes. The overall pattern is one of speed-up accompanied by terminus retreat, but there are also several instances of surge behavior and a few cases of glacier slowdown. Comprehensive mappings such as these, at regular intervals, provide an important new observational capability for understanding ice-sheet variability.