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.

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.

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.

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.

The Seal Nunataks ice shelf (SNIS, ~743 km2 in 2013) is an unofficial name for a remnant area between the former Larsen A and Larsen B ice shelves off the northeastern Antarctic Peninsula. Analyses using Landsat 7 ETM+ and Terra ASTER images from 2001 to 13 and ICESat altimetry from 2003 to 09 show it has retreated and thinned following the Larsen A (1995) and Larsen B (2002) disintegrations. Despite some regional cooling and more fast ice since 2008, SNIS continues to lose ice along its margins and may be losing contact with some nunataks. Detailed analysis of data from four ICESat tracks indicates that ice shelf thinning rates range between 1.9 and 2.7 m a−1, and generally increase from west to east. An ICESat repeat track crossing the adjacent Robertson Island shows a mean elevation loss of 1.8 m a−1. Two tracks crossing the SNIS's remaining tributary, Rogosh Glacier, show sub-meter elevation losses. Comparing shelf remnant and grounded ice thinning rates implies that basal ocean melting augments SNIS thinning by ~1 m a−1, a rate that is consistent with other estimates of ocean-driven shelf thinning in the region.

Glaciers in Greenland are changing rapidly. To better understand these changes, we have produced a series of seven synthetic aperture radar (SAR) backscatter mosaics for seven winters during the period from 2000 to 2013. Six of the mosaics were created using RADARSAT Fine-Beam data and the seventh used ALOS-PALSAR Fine-Beam Single-Polarization data. The RADARSAT mosaics are radiometrically calibrated and capture changes in the backscatter coefficient related to melt and other events, particularly the strong melting in the summer of 2012. Comparison of features in the ascending-orbit ALOS mosaic and the descending-orbit RADARSAT mosaics indicate that in areas of smooth to moderate topography their locations are consistent to within a few tens of meters. The locations of features identifiable in the RADARSAT mosaics, which were collected with the same imaging parameters, generally agree to within better than the 20 m posting of the data. With such geometric accuracy, these data establish a record of change in Greenland for the early part of the 21st century, thus providing a baseline that can be compared with new radar and optical datasets.

Observations of Southern Hemisphere sea ice from passive microwave satellite measurements show that a new record maximum extent of 19.58 x 106 km2 was reached on 30 September 2013; the extent is just over two standard deviations above the 1979-2012 mean and follows a similar record (19.48x 106km2) in 2012. On the record day in 2013, sea-ice extent was greater than the 30 year average (1981-2010) in nearly all Southern Ocean regions. For the year as a whole, Southern Hemisphere sea-ice area and extent were well above average, and numerous monthly and daily records were broken. Analysis of anomaly patterns and the atmospheric and oceanic events suggests that a sequence of regional wind and cold-freshened surface waters is likely responsible for the record maximum and the generally high 2013 extent. In particular, the Ross Sea sector experienced a combination of cold southerly winds associated with the position and depth of the Amundsen Sea low, and lower than normal sea surface temperatures (up to 2°C below normal). The resulting very high anomaly in ice extent in this region was a major component of the overall record maximum.

Blue-ice areas (BIAs) and their geographical distribution in Antarctica were mapped using Landsat-7 ETM+ images with 15 m spatial resolution obtained during the 1999–2003 austral summers and covering the area north of 82.5° S, and a snow grain-size image of the MODIS-based Mosaic of Antarctica (MOA) dataset with 125 m grid spacing acquired during the 2003/04 austral summer from 82.5°S to the South Pole. A map of BIAs was created with algorithms of thresholds based on band ratio and reflectance for ETM+ data and thresholds based on snow grain size for the MOA dataset. The underlying principle is that blue ice can be separated from snow or rock by their spectral discrepancies and by different grain sizes of snow and ice. We estimate the total area of BIAs in Antarctica during the data acquisition period is 234 549 km2, or 1.67% of the area of the continent. Blue ice is scattered widely over the continent but is generally located in coastal or mountainous regions. The BIA dataset presented in this study is the first map covering the entire Antarctic continent sourced solely from ETM+ and MODIS data. This dataset can potentially benefit other studies in glaciology, meteorology, climatology and paleoclimate, meteorite collection and airstrip site selection.

A connected system of active subglacial lakes was revealed beneath Recovery Ice Stream, East Antarctica, by ICESat laser altimetry. Here we combine repeat-track analysis of ICESat (2003–09), Operation IceBridge laser altimetry and radio-echo sounding (2011 and 2012), and MODIS image differencing (2009–2011) to learn more about the lake activity history, the surface and bedrock topographic setting of the lakes and the constraints on water flow through the system. We extend the lake activity time series to 2012 for the three lower lakes and capture two major lake drainages. One lake underwent a large deflation between 2009 and 2011 while another lake, which had been continuously filling between 2003 and 2010, started to drain after 2011. Most of the active lakes are located in a ~ 1000 km long bedrock trough under the main trunk of Recovery Ice Stream, whose base is ~ 1500– 2000 m below present-day sea level. The hydrologic system beneath Recovery Ice Stream is controlled by this unusually pronounced bedrock topography, in contrast to most Antarctic systems studied to date, which are controlled by the ice surface topography. Hydrologic connections among the lakes appear to be direct and responsive, and we reproduce the lake activity using a simple subglacial water model. We discuss potential causes of non-steady hydrologic behavior in major Antarctic catchments.

The Arctic pole of inaccessibility (API), defined as the point on the Arctic Ocean that is farthest from any land, is commonly asserted to lie at 84° 03′ N, 174° 51′ W. We show that the true position is 85° 48′ N, 176° 09′ E, over 200 km from the traditional location. The reason for this error is unknown.

We present an update of the ‘key points’ from the Antarctic Climate Change and the Environment (ACCE) report that was published by the Scientific Committee on Antarctic Research (SCAR) in 2009. We summarise subsequent advances in knowledge concerning how the climates of the Antarctic and Southern Ocean have changed in the past, how they might change in the future, and examine the associated impacts on the marine and terrestrial biota. We also incorporate relevant material presented by SCAR to the Antarctic Treaty Consultative Meetings, and make use of emerging results that will form part of the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report.

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.

We investigate the elevation and mass-balance response of tributary glaciers following the loss of the Larsen A and B ice shelves, Antarctic Peninsula (in 1995 and 2002 respectively). Our study uses MODIS imagery to track ice extent, and ASTER and SPOT5 digital elevation models (DEMs) plus ATM and ICESat laser altimetry to track elevation changes, spanning the period 2001–09. The measured Larsen B tributary glaciers (Hektoria, Green, Evans, Punchbowl, Jorum and Crane) lost up to 160 m in elevation during 2001–06, and thinning continued into 2009. Elevation changes were small for the more southerly Flask and Leppard Glaciers, which are still constrained by a Larsen B ice shelf remnant. In the northern embayment, continued thinning of >3 m a−1 on Drygalski Glacier, 14 years after the Larsen A ice shelf disintegrated, suggests that mass losses for the exposed Larsen B tributaries will continue for years into the future. Grounded ice volume losses exceed 13 km3 for Crane Glacier and 30 km3 for the Hektoria–Green–Evans glaciers. The combined mean loss rate for 2001–06 is at least 11.2 Gt a−1. Our values differ significantly from published mass-budget-based estimates for these embayments, but are a reasonable fraction of GRACE-derived rates for the region (∼40 Gt a−1).

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.

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.

We present an analysis of the active hydrologic system of MacAyeal Ice Stream (MacIS), West Antarctica, from a synthesis of multiple remote-sensing techniques: satellite laser altimetry; satellite image differencing; and hydrologic potential mapping (using a satellite-derived DEM and a bedrock DEM from airborne radio-echo sounding). Combining these techniques augments the information provided by each one individually, and allows us to develop a protocol for studying subglacial hydrologic systems in a holistic manner. Our study reveals five large active subglacial lakes under MacIS, the largest of which undergoes volume changes of at least 1.0 km3. We discuss the hydrologic properties of this system and present evidence for links between the lakes. At least three of the lakes are co-located with sticky spots, i.e. regions of high local basal shear stress. We also find evidence for surface elevation changes due to ice-dynamic effects (not just water movement) caused by changes in basal resistance. Lastly, we show that satellite radar altimetry is of limited use for monitoring lake activity on fast-flowing ice streams with surfaces that undulate on ∼10 km length scales.

We use a combination of satellite techniques (interferometric synthetic aperture radar (InSAR), visible-band imagery, and repeat-track laser altimetry) to develop a benchmark map for the Amery Ice Shelf (AIS) grounding zone (GZ), including its islands and ice rises. The break-in-slope, as an indirect estimate of grounding line location, was mapped for the entire AIS. We have also mapped ∼55% of the landward edge and ∼30% of the seaward edge of the ice shelf flexure boundary for the AIS perimeter. Vertical ice motion from Global Positioning System receivers confirms the location of the satellite-derived GZ in two regions. Our map redefines the extent of floating ice in the south-western AIS and identifies several previously unmapped grounded regions, improving our understanding of the stresses supporting the current dynamical state of the ice shelf. Finally, we identify three along-flow channels in the ice shelf basal topography, approximately 10 km apart, 1.5 km wide and 300–500 m deep, near the southern GZ. These channels, which form at the suture zones between ice streams, may represent zones of potential weakness in the ice shelf and may influence sub-ice-shelf ocean circulation.