Search

We review successes and challenges from five recent subglacial bedrock drilling campaigns intended to find evidence for Antarctic Ice Sheet retreat during warm periods in the geologic past. Insights into times when the polar ice sheets were smaller than present serve as guiding information for modeling efforts that aim to predict the rate and magnitude of future sea level rise that would accompany major retreat of the Antarctic Ice Sheet. One method to provide direct evidence for the timing of deglaciations and minimum extent of prior ice sheets is to extract subglacial bedrock cores for cosmogenic nuclide analysis from beneath the modern ice sheet surface. Here we summarize the lessons learned from five field seasons tasked with obtaining bedrock cores from shallow depths (

Subglacial lake water-volume changes produce ice-elevation anomalies that provide clues about water flow beneath glaciers and ice sheets. Significant challenges remain in the quantitative interpretation of these elevation-change anomalies because the surface expression of subglacial lake activity depends on basal conditions, rate of water-volume change, and ice rheology. To address these challenges, we introduce an inverse method that reconstructs subglacial lake activity from altimetry data while accounting for the effects of viscous ice flow. We use a linearized approximation of a Stokes ice-flow model under the assumption that subglacial lake activity only induces small perturbations relative to a reference ice-flow state. We validate this assumption by accurately reconstructing lake activity from synthetic data that are produced with a fully nonlinear model. We then apply the method to estimate the water-volume changes of several active subglacial lakes in Antarctica by inverting data from NASA's Ice, Cloud, and land Elevation Satellite 2 (ICESat-2) laser altimetry mission. The results show that there can be substantial discrepancies (20% or more) between the inversion and traditional estimation methods due to the effects of viscous ice flow. The inverse method will help refine estimates of subglacial water transport and further constrain the role of subglacial hydrology in ice-sheet evolution.

Stable water isotope records of six firn cores retrieved from two adjacent plateaus on the northern Antarctic Peninsula between 2014 and 2016 are presented and investigated for their connections with firn-core glacio-chemical data, meteorological records and modelling results. Average annual accumulation rates of 2500 kg m−2 a−1 largely reduce the modification of isotopic signals in the snowpack by post-depositional processes, allowing excellent signal preservation in space and time. Comparison of firn-core and ECHAM6-wiso modelled δ18O and d-excess records reveals a large agreement on annual and sub-annual scales, suggesting firn-core stable water isotopes to be representative of specific synoptic situations. The six firn cores exhibit highly similar isotopic patterns in the overlapping period (2013), which seem to be related to temporal changes in moisture sources rather than local near-surface air temperatures. Backward trajectories calculated with the HYSPLIT model suggest that prominent δ18O minima in 2013 associated with elevated sea salt concentrations are related to long-range moisture transport dominated by westerly winds during positive SAM phases. In contrast, a broad δ18O maximum in the same year accompanied by increased concentrations of black carbon and mineral dust corresponds to the advection of more locally derived moisture with northerly flow components (South America) when the SAM is negative.

Converting measurements of ice-sheet surface elevation change to mass change requires measurements of accumulation and knowledge of the evolution of the density profile in the firn. Most firn-densification models are tuned using measured depth–density profiles, a method which is based on an assumption that the density profile in the firn is invariant through time. Here we present continuous measurements of firn-compaction rates in 12 boreholes near the South Pole over a 2 year period. To our knowledge, these are the first continuous measurements of firn compaction on the Antarctic plateau. We use the data to derive a new firn-densification algorithm framed as a constitutive relationship. We also compare our measurements to compaction rates predicted by several existing firn-densification models. Results indicate that an activation energy of 60 kJ mol−1, a value within the range used by current models, best predicts the seasonal cycle in compaction rates on the Antarctic plateau. Our results suggest models can predict firn-compaction rates with at best 7% uncertainty and cumulative firn compaction on a 2 year timescale with at best 8% uncertainty.

A major consideration for maritime activity in the Southern Hemisphere is the northern limit of icebergs, or the Southern Ocean Limit Of Known Ice (SOLOKI). This analysis of historical reports of icebergs during Southern Hemisphere voyages from 1687 to 1933 provides a basis for examination of their geographical and chronological occurrence during ~250 years. The analyses use tabulated data from 742 voyages and other reports from many sources, some including first-person descriptions. While these data are dependent on icebergs being reported by mariners, as well as the variable frequency of voyages, they demonstrate distinct periods of exceptional frequency of icebergs occurring in certain localities, particularly the far South Atlantic. Based upon historical records the evidence suggests unprecedented numbers of icebergs were present in southern shipping channels in the 1890s. When these historical observations are combined with modern iceberg drift trajectories, their possible origin can be elucidated. Owing to the numbers of icebergs seen and their geographical spread, our results suggest that this was possibly the largest near-synchronous calvings in the last 300 years, and the northernmost extent of the SOLOKI.

From 2017 to 2020, three significant calving events took place on Pine Island Glacier, West Antarctica. Ice-shelf velocities changed over this period and the calving events have been suggested as possible drivers. However, satellite observations also show significant changes in the areal extent of fracture zones, especially in the marginal areas responsible for providing lateral support to the ice shelf. Here, we conduct a model study to identify and quantify drivers of recent ice-flow changes of the Pine Island Ice Shelf. In agreement with recent studies, we find that the calving events caused significant velocity changes over the ice shelf. However, calving alone cannot explain observed velocity changes. Changes in the structural rigidity, i.e. ice damage, further significantly impacted ice flow. We suggest that ice damage evolution of the ice-shelf margins may have influenced recent calving events, and these two processes are linked.

Hydrologic pathways beneath ice sheets and glaciers play an important role in regulating ice flow. Antarctica has experienced, and will continue to experience, changes in ice dynamics and geometry, but the associated changes in subglacial hydrology have received less attention. Here, we use the GlaDS subglacial hydrology model to examine drainage evolution beneath an idealised Antarctic glacier in response to steepening ice surface slopes, accelerating ice velocities and subglacial lake drainages. Ice surface slope changes exerted a dominant influence, redirecting basal water to different outlet locations and substantially increasing channelised discharge crossing the grounding line. Faster ice velocities had comparatively negligible effects. Subglacial lake drainage results indicated that lake refilling times play a key role in drainage system evolution, with lake flux more readily accommodated following shorter refilling times. Our findings are significant for vulnerable Antarctic regions currently experiencing dynamic thinning since subglacial water re-routing could destabilise ice shelves through enhanced sub-shelf melting, potentially hastening irreversible retreat. These changes could also affect subglacial lake activity. We, therefore, emphasise that including a nuanced and complex representation of subglacial hydrology in ice-sheet models could provide critical information on the timing and magnitude of sea-level change contributions from Antarctica.

Rifts are full-thickness fractures that propagate laterally across an ice shelf. They cause ice-shelf weakening and calving of tabular icebergs, and control the initial size of calved icebergs. Here, we present a joint inverse and forward computational modeling framework to capture rifting by combining the vertically integrated momentum balance and anisotropic continuum damage mechanics formulations. We incorporate rift–flank boundary processes to investigate how the rift path is influenced by the pressure on rift–flank walls from seawater, contact between flanks, and ice mélange that may also transmit stress between flanks. To illustrate the viability of the framework, we simulate the final 2 years of rift propagation associated with the calving of tabular iceberg A68 in 2017. We find that the rift path can change with varying ice mélange conditions and the extent of contact between rift flanks. Combinations of parameters associated with slower rift widening rates yield simulated rift paths that best match observations. Our modeling framework lays the foundation for robust simulation of rifting and tabular calving processes, which can enable future studies on ice-sheet–climate interactions, and the effects of ice-shelf buttressing on land ice flow.

We developed a novel laser melting sampler (LMS) for ice cores to measure the stable water isotope ratios (δ18O and δD) as temperature proxies at sub-centimeter depth resolutions. In this LMS system, a 2 mm diameter movable evacuation nozzle holds an optical fiber through which a laser beam irradiates the ice core. The movable nozzle intrudes into the ice core, the laser radiation meanwhile melts the ice cylindrically, and the meltwater is pumped away simultaneously through the same nozzle and transferred to a vial for analysis. To avoid isotopic fractionation of the ice through vaporization, the laser power is adjusted to ensure that the temperature of the meltwater is always kept well below its boiling point. A segment of a Dome Fuji shallow ice core (Antarctica), using the LMS, was then demonstrated to have been discretely sampled with a depth resolution as small as 3 mm: subsequent analysis of δ18O, δD, and deuterium excess (d) was consistent with results obtained by hand segmentation within measurement uncertainties. With system software to control sampling resolution, the LMS will enable us to identify temperature variations that may be detectable only at sub-centimeter resolutions in ice cores.

To investigate the mechanisms driving recent changes in outlet glaciers in Antarctica, we measured the glacier front position, flow velocity and surface elevation of five outlet glaciers flowing into Lützow-Holm Bay in East Antarctica. After a steady advance from 2008 to 2015, all the glaciers synchronously retreated by 0.4–6.0 km between 2016 and 2018. The initiation of the retreat coincided with the breakup of land-fast sea ice in Lützow-Holm Bay in 2016, which resulted in the largest sea-ice loss in the region since 1998. Similar flow variations and surface elevation changes were observed near the grounding line of Shirase, Skallen and Telen glaciers. The slowdown in 2011–15 (by 13%) and the speedup in 2016–18 (by 7%) coincided with the respective increase and decrease in surface elevation. Simultaneous retreat and acceleration after the land-fast sea-ice breakup implies that sea ice has a significant influence on glacier dynamics. Thickening/thinning observed near the grounding line was attributed to a reduced/enhanced stretching flow regime during the deceleration/acceleration period. Our results demonstrate that land-fast sea ice affects not only terminus positions, but also the flow speed and ice thickness of the Antarctic glaciers.

Changes in iceberg calving fluxes and oceanographic conditions around Antarctica have likely influenced the spatial and temporal distribution of iceberg fresh water fluxes to the surrounding ocean basins. However, Antarctic iceberg melt rate estimates have been limited to very large icebergs in the open ocean. Here we use a remote-sensing approach to estimate iceberg melt rates from 2011 to 2022 for 15 study sites around Antarctica. Melt rates generally increase with iceberg draft and follow large-scale variations in ocean temperature: maximum melt rates for the western peninsula, western ice sheet, eastern ice sheet and eastern peninsula are ~50, ~40, ~5 and ~5 m a−1, respectively. Iceberg melt sensitivity to thermal forcing varies widely, with a best-estimate increase in melting of ~24 m a−1°C−1 and range from near-zero to ~100 m a−1°C−1. Variations in water shear likely contribute to the apparent spread in thermal forcing sensitivity across sites. Although the sensitivity of iceberg melt rates to water shear prevents the use of melt rates as a proxy to infer coastal water mass temperature variability, additional coastal iceberg melt observations will likely improve models of Southern Ocean fresh water fluxes and have potential for subglacial discharge plume mapping.

Mass loss from the Amundsen Sea Embayment of the West Antarctic Ice Sheet is a major contributor to global sea-level rise (SLR) and has been increasing over recent decades. Predictions of future SLR are increasingly modelled using ensembles of simulations within which model parameters and external forcings are varied within credible ranges. Accurately reporting the uncertainty associated with these predictions is crucial in enabling effective planning for, and construction of defences against, rising sea levels. Calibrating model simulations against current observations of ice-sheet behaviour enables the uncertainty to be reduced. Here we calibrate an ensemble of BISICLES ice-sheet model simulations of ice loss from the Amundsen Sea Embayment using remotely sensed observations of surface elevation and ice speed. Each calibration type is shown to be capable of reducing the 90% credibility bounds of predicted contributions to SLR by 34 and 43% respectively.

One of the largest contributors to uncertainty in predictions of sea-level rise from ice-sheet models is a lack of knowledge about the bed topography beneath ice sheets. Bed topography maps are normally made by interpolating between linear radar surveys using methods that include kriging, mass conservation and flowline diffusion, all of which may miss influential mesoscale (2–30 km) bedforms. Previous works have explored an Ice-Flow Perturbation Analysis (IFPA) approach for estimating bed topography using the surface expression of these mesoscale bedforms. Using regions of Pine Island Glacier that have been intensively surveyed by ice-penetrating radar as test sites, and a refined IFPA methodology, we find that IFPA detects bedforms capable of influencing ice flow which are not represented in Bedmachine Antarctica and other interpolated bed products. We further explore the ability of IFPA to estimate relative bed slipperiness, finding higher slipperiness in the main trunk and tributaries. Alongside other methods which estimate ice thickness, bed topography maps from IFPA have the potential to constrain projections of future sea-level rise, especially where radar data are sparse.

Estimates of ice shelf mass loss are typically based on surface height measurements, assuming hydrostatic equilibrium and estimated firn thickness. Recent investigations, however, challenge the assumption that ice shelves are freely floating, particularly in proximity to narrow structures such as basal channels and shear margins. We compare contemporaneous measurements of Antarctic ice shelf thickness, from ice-penetrating radar, to freeboard height, from laser altimetry, acquired during multiple airborne surveys. On average, the hydrostatic thickness differs from observed thickness by at least ~17 ± 98 m, but this difference varies well beyond the propagated error within and among ice shelves, and depends on the corrections applied. We find that uncertainty in firn thickness can account for most, but not all, of the imbalance. Overall, errors in hydrostatic thickness do not significantly impact estimated basal melt rates. Our results indicate that localized approaches to estimating ice shelf thickness and rates of change are not applicable at large scales, and vice versa, and point to the need for more abundant and accurate firn and ice thickness measurements to improve estimates and predictions of ice shelf mass loss.

As part of the Coupled Model Intercomparison Project Phase 6 (CMIP6), the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) was devised to assess the likely sea-level-rise contribution from the Earth's ice sheets. Here, we construct an ensemble of climate forcings for Antarctica until the year 2300 based on original ISMIP6 forcings until 2100, combined with climate indices from simulations with the MIROC4m climate model until 2300. We then use these forcings to run simulations for the Antarctic ice sheet with the SICOPOLIS model. For the unabated warming pathway RCP8.5/SSP5-8.5, the ice sheet suffers a severe mass loss, amounting to ~ 1.5 m SLE (sea-level equivalent) for the fourteen-experiment mean, and ~ 3.3 m SLE for the most sensitive experiment. Most of this loss originates from West Antarctica. For the reduced emissions pathway RCP2.6/SSP1-2.6, the loss is limited to a three-experiment mean of ~ 0.16 m SLE. The means are approximately two times larger than what was found in a previous study (Chambers and others, 2022, doi:10.1017/jog.2021.124) that assumed a sustained late-21st-century climate beyond 2100, demonstrating the importance of post-2100 climate trends on Antarctic mass changes in the 22nd and 23rd centuries.

Although over 600 Antarctic subglacial lakes have been identified using radar and satellite observations, the bathymetry and bed properties, which are key to understanding conditions within the lake, have been determined in very few localities. We present measurements of water column thickness and lakebed properties from Lago Subglacial CECs (SLC), located beneath 2653 m of ice at the Rutford-Institute-Minnesota divide in Antarctica. Seismic profiles indicate a maximum water column thickness of 301.3 ± 1.5 m, at the widest part of the lake, with an estimated lake volume of 2.5 ± 0.3 km3. Seismic imaging and measurements of the reflection strength at the ice base and lakebed indicate >15 m of high-porosity fine-grained sediment in the central section of the lakebed, consistent with a depositional sequence with an age of up to 0.5 Ma. These observations, along with previous radar measurements and modelling, indicate a low-energy sedimentary environment with a long water-residence time. As such, SLC is a suitable target for exploration via direct access to recover sediment records of ice sheet and climate history and investigate microbial life with long periods of isolation.

The Antarctic Peninsula's widespread glacier retreat and ice shelf collapse have been attributed to atmospheric and oceanic warming. Following the initial post-collapse period of retreat, several former tributary glaciers of the Larsen A and B ice shelves have been slowly re-advancing for more than a decade. Here, we use a flowline model of Crane Glacier to gauge the sensitivity of former tributary glaciers to future climate change following this period of long-term dynamic adjustment. The glacier's long-term geometry and speed changes are similar to those of other former Larsen A and B tributaries, suggesting that Crane Glacier is a reasonable representation of regional dynamics. For the unperturbed climate simulations, discharge remains nearly unchanged in 2018–2100, indicating that dynamic readjustment to shelf collapse took ~15 years. Despite large uncertainties in Crane Glacier's past and future climate forcing, a wide range of future climate scenarios leads to a relatively modest range in grounding line discharge (0.8–1.5 Gt a−1) by 2100. Based on the model results for Crane, we infer that although former ice shelf tributaries may readvance following collapse, similar to the tidewater glacier cycle, their dynamic response to future climate perturbations should be less than their response to ice shelf collapse.

Antarctica's ice shelves stabilize the ice sheet and, therefore, understanding processes affecting the mass budgets of ice shelves is important for estimating grounded ice loss. To study the ice shelf dynamics, we analyzed seismological and GNSS data from the Ekström Ice Shelf in Dronning Maud Land. We extracted probabilistic power spectral densities (PPSD) in the frequency band 3.4–6.8 Hz, typical of icequakes, from seismological data and observed pronounced signals in the PPSD with near 3 and 4 cycles per day (cpd) corresponding to tidal overharmonics, in addition to the main tidal constituents near 1 and 2 cpd. GNSS data reveal the same components in ice flow speed but not in vertical displacements. Generally, tide-induced grounding line migration modulates the flow velocity of an entire ice shelf. We find that this velocity modulation causes the increased icequake activity in the tidal overharmonics with 3 and 4 cpd in an ice shear zone where the flow velocity drops to nearly zero.

Over the last two decades, anomalous warming events have been observed in coastal Antarctic regions. While these events have been documented in the Ross Sea sector, the Antarctic interior is believed to have been buffered from warming. In this work, we present data from lakes located near Mt. Heekin and Thanksgiving Valley (~85° S) along the Shackleton Glacier, which are believed to be the southern-most Antarctic dry valley lakes. In 2018, the lakes were characterized, repeat satellite images were examined, and lake water chemistry was measured. Our analysis shows that lake areas recently increased, and the water-soluble ion chemistry indicates a flushing of salts from periglacial soils, likely from increased glacial melt as illustrated by water isotope data. Our results show that high southern latitude ice-free areas have likely been affected by warm pulses over the past 60 years and these pulses may be quasi-synchronous throughout the Transantarctic Mountains.

At any one time 130 000 icebergs are afloat in the Southern Ocean; 97% of these are too small to be registered in current satellite-based databases, yet the melting of these small icebergs provides a major input to the Southern Ocean. We use a unique set of visual size observations of 53 000 icebergs in the South Atlantic Ocean, the SCAR International Iceberg Database, to derive average iceberg dissolution rates. Fracture into two parts is the dominant dissolution process for tabular icebergs, with an average half-life of 30 days for icebergs

During Cook's 1772–75 Antarctic circumnavigation on the HMS Resolution, he recorded the positions of hundreds of icebergs. This paper compares Cook's observations and those of Halley in 1700, Bouvet in 1739 and Riou in 1789, with the Brigham Young University/National Ice Center (BYU/NIC) and the Alfred Wegener Institute datasets. Cook's description of the iceberg plume east of the Amery Ice Shelf and the iceberg distributions in the Weddell, Ross and Amundsen Seas agree with modern data. In January 1774, Cook reached his farthest south on the shelf of the Amundsen Sea Embayment, the site of the current International Thwaites Glacier Collaboration field study. Cook's largest iceberg had a 2.5 km diameter, where power-law models show that icebergs of this size or smaller comprise 92% of their total number. In the eastern Weddell, Cook's observation of a sea-ice tongue with a much greater extent than in satellite imagery remains unexplained. Although Riou's icebergs lie 1000 km east of the BYU/NIC trajectories, application of the England and others (2020) fracture and drift model to the trajectories removes the discrepancy and means that all the ship observations are consistent with modern observations and theory.

The analysis of Martian meteorites is a key research to understand the mineralogical composition of Mars. However, they suffer different types of alteration due to the environment where they fall on Earth. These differences should be identified in order to characterize correctly the original Martian compounds. Most of the meteorites found on Earth are collected in Antarctica where the environmental conditions are such that, in general and based on many previous investigations, produce fewer alterations compared with other terrestrial environments such as hot deserts. In this study, the weathering alterations of minerals from the MIL 090030 Martian Nakhlite found in Antarctica were analyzed to determine which minerals are formed by Antarctic environmental conditions. It was confirmed that the south-polar region environment can contaminate Martian meteorites and, in this study, several minerals like halite, nitratine and niter, were detected as weathering alteration products that have not been referenced in the literature before.

The Antarctic Peninsula Ice Sheet (APIS) has become a significant contributor to sea-level rise over recent decades. Accurately estimating the ice discharge from the outlet glaciers of the APIS is crucial to quantify the mass balance of the Antarctic Peninsula. We here compute the ice discharge from the outlet glaciers of the APIS north of 70${^\circ }$S for the five most widely used ice-thickness reconstructions, using a common surface velocity field and a common set of flux gates, so the differences in ice discharge can be solely attributed to the differences in ice thickness at the flux gates. The total volumetric ice discharge for 2015–2017 ranges within 45–141 km3 a−1, depending on the ice-thickness model, with a mean of 87 ± 44 km3 a−1. The substantial differences between the ice-discharge results, and a multi-model normalized root-mean-squared deviation of 0.91 for the whole data set, reveal large differences and inconsistencies between the ice-thickness models, giving an indication of the large uncertainty in the current ice-discharge estimates for the APIS. This manifests a fundamental problem of the region: the scarcity of appropriate ice-thickness measurements and the difficulty of the current models to reconstruct the ice-thickness distribution in this complex region.

Hercules Dome, Antarctica, has long been identified as a prospective deep ice core site due to the undisturbed internal layering, climatic setting and potential to obtain proxy records from the Last Interglacial (LIG) period when the West Antarctic ice sheet may have collapsed. We performed a geophysical survey using multiple ice-penetrating radar systems to identify potential locations for a deep ice core at Hercules Dome. The surface topography, as revealed with recent satellite observations, is more complex than previously recognized. The most prominent dome, which we term ‘West Dome’, is the most promising region for a deep ice core for the following reasons: (1) bed-conformal radar reflections indicate minimal layer disturbance and extend to within tens of meters of the ice bottom; (2) the bed is likely frozen, as evidenced by both the shape of the measured vertical ice velocity profiles beneath the divide and modeled ice temperature using three remotely sensed estimates of geothermal flux and (3) models of layer thinning have 132 ka old ice at 45–90 m above the bed with an annual layer thickness of ~1 mm, satisfying the resolution and preservation needed for detailed analysis of the LIG period.

The British Antarctic Survey's (BAS) Rapid Access Isotope Drill (RAID), designed for rapid drilling to survey prospective ice core sites, has been deployed at multiple Antarctic locations over 6 years. This drilling method creates ice chippings that can be discretely sampled and analysed for their chemical and water isotopic composition. Ice sampling methods have evolved since the first uses of the BAS RAID, enabling a more quantifiable sample resolution. Here, we show that water isotope records obtained from RAID ice are comparable to those of equivalent depth resolution from proximal ice cores. Records of chemical impurities also show good agreement with nearby cores. Our findings suggest that the RAID is suitable for both chemical and isotopic reconnaissance of drilling sites. Residual contamination of certain ions is discussed, with proposed design changes to avoid this issue with future use.

Accurate ice flow velocity data are essential for studying the mass balance of the Antarctic ice sheet. However, there is a lack of ice velocity maps of 1960s–80s in basin-wide regions or the entire ice sheet. In this study, an enhanced hierarchical network densification approach is developed for basin-wide Antarctic velocity mapping using historical ARGON and Landsat images. The produced multiple historical velocity maps from 1963 to 1989 in the region of the Fimbul and Jelbart ice shelves, East Antarctica, achieved an accuracy better than 29 m a−1. They revealed that the ice flow velocity had no significant changes over the period. Combining the surface mass balance estimate with the ice discharge estimated from our historical velocity maps and recently published velocity maps, we estimated a positive mass balance of 8.6 ± 3.9 Gt a−1 in the study area from 1963 and 2015. Our results indicate that the region's positive mass balance, as estimated in recently published studies, has been maintained since the 1960s. It is also in concordance with the low level of mass balance from 1992 to 2017 in East Antarctica. This suggests that the study area has been stable since the 1960s.

The retreat rates of Triangular Glacier since 1979 and its mass changes during the period 2014/15–2019/20 indicate the sensitive response of small ice masses on the eastern side of the Antarctic Peninsula to air temperature evolution. This cirque glacier in the northern part of James Ross Island receded rapidly during the period of regional warming in the late 20th century, losing 30.8% of its surface area between 1979 and 2006 (−1.7% a−1). The retreat rate then dropped to −0.3% a−1 following the regional cooling trend, but started to accelerate again (−0.8 to −2.3% a−1) with increasing air temperature since the summer 2014/15. Since the glaciological year 2015/16, Triangular Glacier has experienced enhanced snow melt, wind scour and permanent mass loss with annual mass balance ranging from −0.08 ± 0.35 to −0.56 ± 0.25 m w.e. The largest mass loss was observed in the glaciological year 2019/20, which included the warmest summer of the observation period. The cumulative mass balance of −1.66 ± 0.83 m w.e. over the years 2014/15–2019/20 is consistent with the termination of the positive mass-balance period that occurred in the north-eastern Antarctic Peninsula from 2009/10 to 2014/15.

The future contributions of the Antarctic Ice Sheet to sea level rise will depend on the evolution of its surface mass balance (SMB), which could amplify/dampen mass losses increasingly observed at the ice sheet's edge. In situ constraints of SMB over annual-to-decadal timescales consist mostly of firn/ice cores that have a surface footprint $\sim$cm$^{2}$. SMB constraints also come from climate models, which have a higher temporal resolution but a larger surface footprint of several km$^{2}$. We use ice-penetrating radar data to obtain an intermediate spatial and temporal resolution SMB record over three ice rises along the Princess Ragnhild Coast. The co-located ice cores allow us to obtain absolute radar-derived SMB rates at a multi-annual-to-decadal temporal resolution. By comparing the ice core SMB measurements and the radar-derived SMB records, we determine that pointwise measurements of SMB are representative of a small surface area, $\sim 200-500$ m radius extending from the ice core drill site for the ice rises studied here, and that the pointwise measurements are systematically 7–15 cm w.e. a$^{-1}$ lower than the mean SMB value calculated for the whole ice rises. However, ice core records are representative of an entire ice rise's temporal variability at the temporal resolution examined.

Since the mid-20th century, ice shelves around the Antarctic Peninsula have declined in extent and thickness, and some have shown signs of structural instability. Here, using satellite imagery from 1999/2000 to 2019/20 (Landsat 7 and 8, Sentinel-2 and ASTER), we measure areal changes, calculate surface flow speeds, and quantify structural changes of Bach, Stange and George VI ice shelves, located in the southwest Antarctic Peninsula. We recorded a total area loss of 797.5 km2 from 2009/10 to 2019/20, though spatial and temporal patterns varied at individual ice fronts. Flow speeds remained largely stable over the observation periods, but notable acceleration was calculated for Bach Ice Shelf, and at the northern and southern extents of George VI Ice Shelf. Open fractures widened and lengthened between 2009/10 and 2019/20 on all three ice shelves. We conclude that Stange Ice Shelf is stable, and not under any immediate threat of enhanced recession. Continued ice-mass loss and consequential speed up of George VI South may cause further fracturing and destabilisation in the coming decades. Of more immediate concern are the glaciological changes noted for Bach Ice Shelf and George VI North; substantial areas of stabilising ice have already, or will soon be removed, that may lead to enhanced recession within the next decade.

Sea-ice drift in the Antarctic marginal ice zone (MIZ) is discussed using data from a 4-month-long drift of a buoy deployed on a pancake ice floe during the winter sea-ice expansion. We demonstrate increased meandering and drift speeds, and changes in the dynamical regimes of the absolute dispersion during cyclone activity, together with high correlations between drift velocities and wind from atmospheric reanalyses. This indicates a dominant physical control of wind forcing on ice drift and the persistence of free-drift conditions. These conditions occurred despite the buoy remaining largely in >80% ice concentrations and at distances >200 km from the estimated ice edge. The drift is additionally characterised by a strong inertial signature at 13.47 h, which appears initiated by passing cyclones. A wavelet analysis of the buoy's velocity confirms that the momentum transfer from winds at the multi-day frequencies is due to atmospheric forcing, while the initiation of inertial oscillations of sea ice has been identified as the secondary effect. Propagating storm-generated waves may initiate inertial oscillations by increasing the mobility of floes and enhance the drag of the inertial current. This analysis indicates that the Antarctic MIZ in the Indian Ocean sector remains much wider and mobile, during austral winter-to-spring, than defined by sea-ice concentration.

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.

Ice cliffs are features along ice sheet margins, along tropical mountain glaciers, at termini of mountain glaciers and on debris-covered glacier tongues, that have received scattered attention in literature. They cover small relative areas of glacier or margin surface respectively, but have been involved in two apparent anomalies. On the one hand, they have been identified as potential hotspots of extreme melt rates on debris-covered tongues contributing to their relatively rapid ablation, compared to the surrounding glacier surface. On the other hand, they appear where the ice margin is stable (or temporarily advancing) even under conditions of negative mass balance. In this manuscript, we recapitulate why ice cliffs remain interesting features to investigate and what we know about them so far. We conclude by suggesting to further investigate their genesis and variable morphology and their potential as windows into past climates and processes.

During the second half of the 20th century, the Antarctic Peninsula region has undergone a long and sustained warming period, followed by a shorter but also sustained cooling period, and then a very recent return to warming conditions. All of these have profoundly impacted the glaciers peripheral to the Antarctic Peninsula. This paper focuses on the analysis of the surface mass balance monitoring of such glaciers by the glaciological method, complemented by the analysis of mass-balance estimates by geodetic methods, as well as frontal ablation estimates. We aim to summarize the current knowledge and outline the main challenges faced by investigating the mass balance of such peripheral glaciers and their current contribution to sea-level rise.

In the evening of 24 December 1789, 2100 km southeast of Cape Town and after encountering three icebergs, the HMS Guardian under Captain Edward Riou collided with the submerged foot of a large iceberg. Despite severe damage to the ship and its abandonment by many of its crew and passengers, Riou sailed the hulk back to Cape Town, arriving on 22 February 1790. From present-day research and field studies, the formation of the foot in the collision is consistent with the above-freezing seawater temperatures inferred from Riou's commentary. Further, the observed 60 m iceberg height suggests that it calved from the Filchner Ice Shelf in the Weddell Sea. Comparison of the positions of Riou's icebergs with historic sightings, satellite observations and iceberg drift and fracture models also shows that they originated in the Weddell Sea and that their likelihood of occurrence in the collision region is small.

Subglacial hydrology is an important component of the ice dynamic system in Antarctica but is challenging to investigate due to the large spatial scales of the catchment systems, the ice thickness, and remote location. Here I discuss key discoveries about Antarctic subglacial hydrology from the Glacier Drainage System (GlaDS) model, including the presence of long, often high-pressure, subglacial channels. These channels pump tens of cubic metres per second of freshwater into ice-shelf cavities and directly affect melt rates at the critical grounding zone regions. Future ice dynamics and ice-shelf cavity models should take subglacial hydrology into account if they are to accurately predict future behaviour of the Antarctic Ice Sheet.

Firn is the pervasive surface material across Antarctica, and its structures reflect its formation and history in response to environmental perturbations. In addition to the role of firn in thermally isolating underlying glacial ice, it defines near-surface elastic and density structure and strongly influences high-frequency (> 5 Hz) seismic phenomena observed near the surface. We investigate high-frequency seismic data collected with an array of seismographs deployed on the West Antarctic Ice Sheet (WAIS) near WAIS Divide camp in January 2019. Cross-correlations of anthropogenic noise originating from the approximately 5 km-distant camp were constructed using a 1 km-diameter circular array of 22 seismographs. We distinguish three Rayleigh (elastic surface) wave modes at frequencies up to 50 Hz that exhibit systematic spatially varying particle motion characteristics. The horizontal-to-vertical ratio for the second mode shows a spatial pattern of peak frequencies that matches particle motion transitions for both the fundamental and second Rayleigh modes. This pattern is further evident in the appearance of narrow band spectral peaks. We find that shallow lateral structural variations are consistent with these observations, and model spectral peaks as Rayleigh wave amplifications within similarly scaled shallow basin-like structures delineated by the strong velocity and density gradients typical of Antarctic firn.

Ice-sheet simulations of Antarctica extending to the year 3000 are analysed to investigate the long-term impacts of 21st-century warming. Climate projections are used as forcing until 2100 and afterwards no climate trend is applied. Fourteen experiments are for the ‘unabated warming’ pathway, and three are for the ‘reduced emissions’ pathway. For the unabated warming path simulations, West Antarctica suffers a much more severe ice loss than East Antarctica. In these cases, the mass loss amounts to an ensemble average of ~3.5 m sea-level equivalent (SLE) by the year 3000 and ~5.3 m for the most sensitive experiment. Four phases of mass loss occur during the collapse of the West Antarctic ice sheet. For the reduced emissions pathway, the mean mass loss is ~0.24 m SLE. By demonstrating that the consequences of the 21st century unabated warming path forcing are large and long term, the results present a different perspective to ISMIP6 (Ice Sheet Model Intercomparison Project for CMIP6). Extended ABUMIP (Antarctic BUttressing Model Intercomparison Project) simulations, assuming sudden and sustained ice-shelf collapse, with and without bedrock rebound, corroborate a negative feedback for ice loss found in previous studies, where bedrock rebound acts to slow the rate of ice loss.

The Wilkes Subglacial Basin in East Antarctica contains ice equivalent to 3–4 m of global mean sea level rise and is primarily drained by Cook Glacier. Of concern is that recent observations (since the 1970s) show an acceleration in ice speed over the grounding line of both the Eastern and Western portions of Cook Glacier. Here, we use a numerical ice-flow model (Úa) to simulate the instantaneous effects of observed changes at the terminus of Cook Glacier in order to understand the link between these changes and recently observed ice acceleration. Simulations suggest that the acceleration of Cook West was caused by a retreat in calving-front position in the 1970s, potentially enhanced by grounding-line retreat, while acceleration of Cook East was likely caused by ice-shelf thinning and grounding-line retreat in the mid-1990s. Moreover, we show that the instantaneous ice discharge at Cook East would increase by up to 85% if the whole ice shelf is removed and it ungrounds from a pinning point; and that the discharge at Cook West could increase by ~300% if its grounding line retreated by 10 km.

Surface melting on Amery Ice Shelf (AIS), East Antarctica, produces an extensive supraglacial drainage system consisting of hundreds of lakes connected by surface channels. This drainage system forms most summers on the southern portion of AIS, transporting meltwater large distances northward, toward the ice front and terminating in lakes. Here we use satellite imagery, Landsat (1, 4 and 8), MODIS multispectral and Sentinel-1 synthetic aperture radar to examine the seasonal and interannual evolution of the drainage system over nearly five decades (1972–2019). We estimate seasonal meltwater input to one lake by integrating output from the regional climate model [Regional Atmospheric Climate Model (RACMO 2.3p2)] over its catchment defined using the Reference Elevation Model of Antarctica. We find only weak positive relationships between modeled seasonal meltwater input and lake area and between meltwater input and lake volume. Consecutive years of extensive melting lead to year-on-year expansion of the drainage system, potentially through a link between melt production, refreezing in firn and the maximum extent of the lakes at the downstream termini of drainage. These mechanisms are important when evaluating the potential of drainage systems to grow in response to increased melting, delivering meltwater to areas of ice shelves vulnerable to hydrofracture.

In the afternoon of 30 January 2017, a catastrophic outburst flood occurred in the Larsemann Hills (Princess Elizabeth Land, East Antarctica). The rapid drainage of both a thin supraglacial layer of water (near Boulder Lake) and Lake Ledyanoe into the englacial Lake Dålk provoked its overfill and outburst. As a result, a depression of 183 m × 220 m was formed in the place where Lake Dålk was located. This study summarises and clarifies the current state of knowledge on the flood that occurred in 2017. We present a phenomenological model of depression formation. We specify the reasons for the outburst of the system of lakes Boulder, Ledyanoe and Dålk. In addition, we carry out mathematical modelling of the outburst of each of the three lakes and estimate the flood severity. Outburst hydrographs, channel diameters, volume and duration of floods were calculated. Particular simulation results were validated with field data. In conclusion, we give an overview of the new outburst cycle of the lake system, which began in 2020 with the drainage of the lakes Boulder and Ledyanoe, and the new formation of Lake Dålk. Further research is required to improve our understanding of the lake system responses to changing external factors.

Earlier gas measurements of firn air (atmosphere in open pore channels) at polar sites have revealed the occurrence of gas fractionation phenomena during bubble close-off, in addition to well-known thermal and gravitational gas separation. Nevertheless, because of difficulties posed by measurement, little is known about the distribution of air constituents in already closed pores (bubbles) in firn. Herein, we describe the application of high-sensitivity pure rotational Raman spectroscopy, combined with sample immersion in the fluorocarbon-based inert fluid for removing the optical disturbance by diffused reflection. That application efficiently elicits information about nitrogen and oxygen composition ratios (N2/O2 or O2/N2) for each air bubble in firn. The developed methodology presents important implications for elucidating how gas records are formed and modified in the course of pore close-off in polar firn.

The Japanese second deep ice coring project was carried out at Dome Fuji, Antarctica. Following the drilling of the pilot hole in 2001, deep ice core drilling led by the Japanese Antarctic Research Expedition (JARE) was conducted over four austral summer seasons, beginning with the 2003/04 season and reached a depth of 3035.22 m near the bedrock in January 2007. The new drill was designed and developed with the goals of (1) solving the problems encountered during the first JARE deep coring drill and (2) achieving more efficient drilling. In particular, the maximum core length that can be drilled at one time was increased from 2.30 m to 3.84 m and the chip storage efficiency was enhanced by a special pipe with many small holes. This paper gives an outline of the improved drilling system, the progress of drilling and various drilling data.

Rapid Access Ice Drill is a new drilling technology capable of quickly accessing the glacial bed of Antarctic ice sheets, retrieving ice core and rock core samples, and providing boreholes for downhole logging of physical properties. Scientific goals include searching for old ice near the glacial bed and sampling subglacial bedrock. During field trials near McMurdo Station on a piedmont glacier at Minna Bluff in the 2019–20 austral summer, we successfully completed a ‘top-to-bottom’ operational sequence in three boreholes by (1) augering through firn, (2) creating a borehole packer seal in non-porous ice, (3) establishing fluid circulation, (4) quickly drilling a borehole in ice at penetration rates up to 1.2 m min−1, (5) acquiring a short ice core at depth, (6) penetrating the glacial bed at a depth of ~677 m, (7) recovering a 3.2 m core of ice, basal till and subglacial bedrock, (8) optically logging the borehole on wireline, (9) testing hydrofracture potential by overpressuring the borehole fluid and (10) operating in an environmentally benign yet rapid field mode. Minna Bluff testing, therefore, demonstrates the effectiveness of this integrated system to drill rapidly through thick ice and penetrate across the glacial bed to take cores of bedrock.

Three holes were drilled to the bed of Rutford Ice Stream, through ice up to 2154 m thick, to investigate the basal processes and conditions associated with fast ice flow and the glacial history of the West Antarctic Ice Sheet. A narrative of the drilling, measuring and sampling activities, as well as some preliminary results and initial interpretations of subglacial conditions, is given. These were the deepest subglacial access holes ever drilled using the hot-water drilling method. Samples of bed and englacial sediments were recovered, and a number of instruments were installed in the ice column and the bed. The ice–bed interface was found to be unfrozen, with an existing, well-developed subglacial hydrological system at high pressure, within ~1% of the ice overburden. The bed itself comprises soft, water-saturated sediments, consistent with previous geophysical interpretations. Englacial sediment quantity varies significantly between two locations ~2 km apart, and possibly over even shorter (~20 m) distances. Difficulties and unusual observations encountered while connecting to the subglacial hydrological system in one hole possibly resulted from the presence of a large clast embedded in the bottom of the ice.

The Subglacial Antarctic Lakes Scientific Access (SALSA) Project accessed Mercer Subglacial Lake using environmentally clean hot-water drilling to examine interactions among ice, water, sediment, rock, microbes and carbon reservoirs within the lake water column and underlying sediments. A ~0.4 m diameter borehole was melted through 1087 m of ice and maintained over ~10 days, allowing observation of ice properties and collection of water and sediment with various tools. Over this period, SALSA collected: 60 L of lake water and 10 L of deep borehole water; microbes >0.2 μm in diameter from in situ filtration of ~100 L of lake water; 10 multicores 0.32–0.49 m long; 1.0 and 1.76 m long gravity cores; three conductivity–temperature–depth profiles of borehole and lake water; five discrete depth current meter measurements in the lake and images of ice, the lake water–ice interface and lake sediments. Temperature and conductivity data showed the hydrodynamic character of water mixing between the borehole and lake after entry. Models simulating melting of the ~6 m thick basal accreted ice layer imply that debris fall-out through the ~15 m water column to the lake sediments from borehole melting had little effect on the stratigraphy of surficial sediment cores.

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.

During the last few decades, bed-elevation profiles from radar sounders have been used to quantify bed roughness. Various methods have been employed, such as the ‘two-parameter’ technique that considers vertical and slope irregularities in topography, but they struggle to incorporate roughness at multiple spatial scales leading to a breakdown in their depiction of bed roughness where the relief is most complex. In this article, we describe a new algorithm, analogous to wavelet transformations, to quantify the bed roughness at multiple scales. The ‘Self-Adaptive Two-Parameter’ system calculates the roughness of a bed profile using a frequency-domain method, allowing the extraction of three characteristic factors: (1) slope, (2) skewness and (3) coefficient of variation. The multi-scale roughness is derived by weighted-summing of these frequency-related factors. We use idealized bed elevations to initially validate the algorithm, and then actual bed-elevation data are used to compare the new roughness index with other methods. We show the new technique is an effective tool for quantifying bed roughness from radar data, paving the way for improved continental-wide depictions of bed roughness and incorporation of this information into ice flow models.

Ice stream discharge responds to a balance between gravity, basal friction and lateral drag. Appreciable viscous heating occurs in shear margins between ice streams and adjacent slow-moving ice ridges, altering the temperature-dependent viscosity distribution that connects lateral drag to marginal strain rates and ice stream velocity. Warmer ice deforms more easily and accommodates faster flow, whereas cold ice supplied from ice ridges drives advective cooling that counteracts viscous heating. Here, we present a two-dimensional (three velocity component), steady-state model designed to explore the thermal controls on ice stream shear margins. We validate our treatment through comparison with observed velocities for Bindschadler Ice Stream and verify that calculated temperatures are consistent with results from previous studies. Sweeping through a parameter range that encompasses conditions representative of ice streams in Antarctica, we show that modeled steady-state velocity has a modest response to different choices in forcing up until temperate zones develop in the shear margins. When temperate zones are present, velocity is much more sensitive to changes in forcing. We identify key scalings for the emergence of temperate conditions in our idealized treatment that can be used to identify where thermo-mechanical feedbacks influence the evolution of the ice sheet.

GRACE and ICESat Antarctic mass-balance differences are resolved utilizing their dependencies on corrections for changes in mass and volume of the same underlying mantle material forced by ice-loading changes. Modeled gravimetry corrections are 5.22 times altimetry corrections over East Antarctica (EA) and 4.51 times over West Antarctica (WA), with inferred mantle densities 4.75 and 4.11 g cm−3. Derived sensitivities (Sg, Sa) to bedrock motion enable calculation of motion (δB0) needed to equalize GRACE and ICESat mass changes during 2003–08. For EA, δB0 is −2.2 mm a−1 subsidence with mass matching at 150 Gt a−1, inland WA is −3.5 mm a−1 at 66 Gt a−1, and coastal WA is only −0.35 mm a−1 at −95 Gt a−1. WA subsidence is attributed to low mantle viscosity with faster responses to post-LGM deglaciation and to ice growth during Holocene grounding-line readvance. EA subsidence is attributed to Holocene dynamic thickening. With Antarctic Peninsula loss of −26 Gt a−1, the Antarctic total gain is 95 ± 25 Gt a−1 during 2003–08, compared to 144 ± 61 Gt a−1 from ERS1/2 during 1992–2001. Beginning in 2009, large increases in coastal WA dynamic losses overcame long-term EA and inland WA gains bringing Antarctica close to balance at −12 ± 64 Gt a−1 by 2012–16.