Abstract Understanding the coastal zone of the Antarctic Ice Sheet (AIS), where it interacts with the Southern Ocean and warmer air masses, is crucial for predicting Antarctica's influence on the global climate and sea level. This region has multiple tipping mechanisms that could trigger large, rapid, and potentially irreversible changes in the AIS, the Southern Ocean and their global connections in the coming centuries. The AIS remains the largest source of uncertainty in future sea‐level projections. Bed topography beneath the ice shelves and the coastal ice sheet is not yet well documented, and is a major source of this uncertainty. This review assesses current knowledge of the coastal zone and highlights methods to investigate it, including aerogeophysical surveys, ground‐ and ship‐based measurements, satellite observations, and computer modeling. An ensemble analysis of published bed topography data sets identifies significant data gaps and their regional distribution, framed in the context of current ice‐sheet behavior and potential instability. We propose scientific priorities and guidelines for future aerogeophysical surveys, advocating for a comprehensive, coordinated international effort to build a next‐generation data set of Antarctic bed properties. Such an initiative would significantly advance understanding of the role of coastal processes in ice‐sheet dynamics, reducing uncertainties in sea‐level rise projections and improving predictions of future ocean and climate changes.

The process of laser light reflecting from surfaces made of scattering materials that do not strongly absorb at the wavelength of the laser can involve reflections from hundreds or thousands of individual grains, which can introduce delays in the time between light entering and leaving the surface. These time-of-flight biases depend on the grain size and density of the medium, and thus they can result in spatially and temporally varying surface height biases estimated from laser altimeters, such as NASA's ICESat-2 (Ice Cloud, and land Elevation Satellite-2) mission. Modeling suggests that ICESat-2 might experience a bias difference as large as 0.1–0.2 m between coarse-grained melting snow and fine-grained wintertime snow (Smith et al., 2018), which exceeds the mission's requirement to measure seasonal height differences to an accuracy better than 0.1 m (Markus et al., 2017). In this study, we investigate these biases using a model of subsurface scattering, laser altimetry measurements from NASA's ATM (Airborne Topographic Mapper) system, and grain size estimates based on optical imagery of the ice sheet. We demonstrate that distortions in the shapes of waveforms measured using ATM are related to the optical grain size of the surface estimated using optical reflectance measurements and show that they can be used to estimate an effective grain radius for the surface. Using this effective grain radius as a proxy for the severity of subsurface scattering, we use our model with grain size estimates from optical imagery to simulate corrections for biases in ICESat-2 data due to subsurface scattering and demonstrate that, on the basis of large-scale averages, the corrections calculated based on the satellite optical imagery match the biases in the data. This work demonstrates that waveform-based altimetry data can measure the optical properties of granular surfaces and that corrections based on optical grain size estimates can correct for subsurface-scattering biases in ICESat-2 data.

We present Bedmap3, the latest suite of gridded products describing surface elevation, ice-thickness and the seafloor and subglacial bed elevation of the Antarctic south of 60 °S. Bedmap3 incorporates and adds to all post-1950s datasets previously used for Bedmap2, including 84 new aero-geophysical surveys by 15 data providers, an additional 52 million data points and 1.9 million line-kilometres of measurement. These efforts have filled notable gaps including in major mountain ranges and the deep interior of East Antarctica, along West Antarctic coastlines and on the Antarctic Peninsula. Our new Bedmap3/RINGS grounding line similarly consolidates multiple recent mappings into a single, spatially coherent feature. Combined with updated maps of surface topography, ice shelf thickness, rock outcrops and bathymetry, Bedmap3 reveals in much greater detail the subglacial landscape and distribution of Antarctica’s ice, providing new opportunities to interpret continental-scale landscape evolution and to model the past and future evolution of the Antarctic ice sheets.

Differential penetration of green laser light into snow and ice has long been considered a possible cause of range and thus elevation bias in laser altimeters. Over snow, ice, and water, green photons can penetrate the surface and experience multiple scattering events in the subsurface volume before being scattered back to the surface and subsequently the instrument's detector, therefore biasing the range of the measurement. Newly formed sea ice adjacent to open-water leads provides an opportunity to identify differential penetration without the need for an absolute reference surface or dual-color lidar data. We use co-located, coincident high-resolution natural-color imagery and airborne lidar data to identify surface and ice types and evaluate elevation differences between those surfaces. The lidar data reveals that apparent elevations of thin ice and finger-rafted thin ice can be several tens of centimeters below the water surface of surrounding leads, but not over dry snow. These lower elevations coincide with broadening of the laser pulse, suggesting that subsurface volume scattering is causing the pulse broadening and elevation shift. To complement our analysis of pulse shapes and help interpret the physical mechanism behind the observed elevation biases, we match the waveform shapes with a model of scattering of light in snow and ice that predicts the shape of lidar waveforms reflecting from snow and ice surfaces based on the shape of the transmitted pulse, the surface roughness, and the optical scattering properties of the medium. We parameterize the scattering in our model based on the scattering length Lscat, the mean distance a photon travels between isotropic scattering events. The largest scattering lengths are found for thin ice that exhibits the largest negative elevation biases, where scattering lengths of several centimeters allow photons to build up considerable range biases over multiple scattering events, indicating that biased elevations exist in lower-level Airborne Topographic Mapper (ATM) data products. Preliminary analysis of ICESat-2 ATL10 data shows that a similar relationship between subsurface elevations (restored negative freeboard) and “pulse width” is present in ICESat-2 data over sea ice, suggesting that biased elevations caused by differential penetration likely also exist in lower-level ICESat-2 data products. The spatial correlation of observed differential penetration in ATM data with surface and ice type suggests that elevation biases could also have a seasonal component, increasing the challenge of applying a simple bias correction.

One of the key components of this research has been the mapping of Antarctic bed topography and ice thickness parameters that are crucial for modelling ice flow and hence for predicting future ice loss and the ensuing sea level rise. Supported by the Scientific Committee on Antarctic Research (SCAR), the Bedmap3 Action Group aims not only to produce new gridded maps of ice thickness and bed topography for the international scientific community, but also to standardize and make available all the geophysical survey data points used in producing the Bedmap gridded products. Here, we document the survey data used in the latest iteration, Bedmap3, incorporating and adding to all of the datasets previously used for Bedmap1 and Bedmap2, including ice bed, surface and thickness point data from all Antarctic geophysical campaigns since the 1950s. More specifically, we describe the processes used to standardize and make these and future surveys and gridded datasets accessible under the Findable, Accessible, Interoperable, and Reusable (FAIR) data principles. With the goals of making the gridding process reproducible and allowing scientists to re-use the data freely for their own analysis, we introduce the new SCAR Bedmap Data Portal (https://bedmap.scar.org, last access: 1 March 2023) created to provide unprecedented open access to these important datasets through a web-map interface. We believe that this data release will be a valuable asset to Antarctic research and will greatly extend the life cycle of the data held within it. Data are available from the UK Polar Data Centre: https://data.bas.ac.uk (last access: 5 May 2023​​​​​​​). See the Data availability section for the complete list of datasets.

Over the past 60 years, scientists have strived to understand the past, present and future of the Antarctic Ice Sheet. One of the key components of this research has been the mapping of Antarctic bed topography and ice thickness parameters that are crucial for modelling ice flow and hence for predicting future ice loss and ensuing sea level rise. Supported by the Scientific Committee on Antarctic Research (SCAR), the Bedmap3 Action Group aims not only to produce new gridded maps of ice thickness and bed topography for the international scientific community, but also to standardize and make available all the geophysical survey data points used in producing the Bedmap gridded products. Here, we document the survey data used in the latest iteration, Bedmap3, incorporating and adding to all of the datasets previously used for Bedmap1 and Bedmap2, including ice-bed, surface and thickness point data from all Antarctic geophysical campaigns since the 1950s. More specifically, we describe the processes used to standardize and make these and future survey and gridded datasets accessible under the ‘Findable, Accessible, Interoperable and Reusable’ (FAIR) data principles. With the goals to make the gridding process reproducible and to allow scientists to re-use the data freely for their own analysis, we introduce the new SCAR Bedmap Data Portal (bedmap.scar.org, last access: 18 October 2022) created to provide unprecedented open access to these important datasets, through a user-friendly webmap interface. We believe that this data release will be a valuable asset to Antarctic research and will greatly extend the life cycle of the data held within it. Data are available from the UK Polar Data Centre: https://data.bas.ac.uk.

Seasonal meltwater pools on the surface of the Greenland Ice Sheet (GrIS) during late spring and summer in lakes on the surface and transforms the ice sheet's surface into a wet environment in the ablation zone below the equilibrium line. These supraglacial lakes in topographic lows on the ice surface are connected by a dendritic pattern of meandering streams and channels that together form a hydrological system consisting of supra-, en-, and subglacial components. Here, we use lidar data from NASA's Airborne Topographic Mapper (ATM) instrument suite and high-resolution optical imagery collected as part of Operation IceBridge (OIB) in spring 2019 over the GrIS to develop methods for the study of supraglacial hydrological features. While airborne surveys have a limited temporal and spatial coverage compared to imaging spaceborne sensors, their high footprint density and high-resolution imagery reveal a level of detail that is currently not obtainable from spaceborne measurements. The accuracy and resolution of airborne measurements complement spaceborne measurements, can support calibration and validation of spaceborne methods, and provide information necessary for high-resolution process studies of the supraglacial hydrological system on the GrIS that currently cannot be achieved from spaceborne observations alone.

The basal thermal state (frozen or thawed) of the Greenland Ice Sheet is under-constrained due to few direct measurements, yet knowledge of this state is becoming increasingly important to interpret modern changes in ice flow. The first synthesis of this state relied on inferences from widespread airborne and satellite observations and numerical models, for which most of the underlying datasets have since been updated. Further, new and independent constraints on the basal thermal state have been developed from analysis of basal and englacial reflections observed by airborne radar sounding. Here we synthesize constraints on the Greenland Ice Sheet's basal thermal state from boreholes, thermomechanical ice-flow models that participated in the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6; Coupled Model Intercomparison Project Phase 6), IceBridge BedMachine Greenland v4 bed topography, Making Earth Science Data Records for Use in Research Environments (MEaSUREs) Multi-Year Greenland Ice Sheet Velocity Mosaic v1 and multiple inferences of a thawed bed from airborne radar sounding. Most constraints can only identify where the bed is likely thawed rather than where it is frozen. This revised synthesis of the Greenland likely Basal Thermal State version 2 (GBaTSv2) indicates that 33 % of the ice sheet's bed is likely thawed, 40 % is likely frozen and the remainder (28 %) is too uncertain to specify. The spatial pattern of GBaTSv2 is broadly similar to the previous synthesis, including a scalloped frozen core and thawed outlet-glacier systems. Although the likely basal thermal state of nearly half (46 %) of the ice sheet changed designation, the assigned state changed from likely frozen to likely thawed (or vice versa) for less than 6 % of the ice sheet. This revised synthesis suggests that more of northern Greenland is likely thawed at its bed and conversely that more of southern Greenland is likely frozen, both of which influence interpretation of the ice sheet's present subglacial hydrology and models of its future evolution. The GBaTSv2 dataset, including both code that performed the analysis and the resulting datasets, is freely available at https://doi.org/10.5281/zenodo.6759384 (MacGregor, 2022).

Abstract. Over the past 60 years, scientists have strived to understand the past, present and future of the Antarctic Ice Sheet. One of the key components of this research has been the mapping of Antarctic bed topography and ice thickness parameters that are crucial for modelling ice flow and hence for predicting future ice loss and ensuing sea level rise. Supported by the Scientific Committee on Antarctic Research (SCAR), the Bedmap3 Action Group aims not only to produce new gridded maps of ice thickness and bed topography for the international scientific community, but also to standardize and make available all the geophysical survey data points used in producing the Bedmap gridded products. Here, we document the survey data used in the latest iteration, Bedmap3, incorporating and adding to all of the datasets previously used for Bedmap1 and Bedmap2, including ice-bed, surface and thickness point data from all Antarctic geophysical campaigns since the 1950s. More specifically, we describe the processes used to standardize and make these and future survey and gridded datasets accessible under the ‘Findable, Accessible, Interoperable and Reusable’ (FAIR) data principles. With the goals to make the gridding process reproducible and to allow scientists to re-use the data freely for their own analysis, we introduce the new SCAR Bedmap Data Portal (bedmap.scar.org, last access: 18 October 2022) created to provide unprecedented open access to these important datasets, through a user-friendly webmap interface. We believe that this data release will be a valuable asset to Antarctic research and will greatly extend the life cycle of the data held within it. Data are available from the UK Polar Data Centre: https://data.bas.ac.uk.

The National Aeronautics and Space Administration (NASA)’s Operation IceBridge (OIB) was a 13‐year (2009–2021) airborne mission to survey land and sea ice across the Arctic, Antarctic, and Alaska. Here, we review OIB’s goals, instruments, campaigns, key scientific results, and implications for future investigations of the cryosphere. OIB’s primary goal was to use airborne laser altimetry to bridge the gap in fine‐resolution elevation measurements of ice from space between the conclusion of NASA’s Ice, Cloud, and land Elevation Satellite (ICESat; 2003–2009) and its follow‐on, ICESat‐2 (launched 2018). Additional scientific requirements were intended to contextualize observed elevation changes using a multisensor suite of radar sounders, gravimeters, magnetometers, and cameras. Using 15 different aircraft, OIB conducted 968 science flights, of which 42% were repeat surveys of land ice, 42% were surveys of previously unmapped terrain across the Greenland and Antarctic ice sheets, Arctic ice caps, and Alaskan glaciers, and 16% were surveys of sea ice. The combination of an expansive instrument suite and breadth of surveys enabled numerous fundamental advances in our understanding of the Earth’s cryosphere. For land ice, OIB dramatically improved knowledge of interannual outlet‐glacier variability, ice‐sheet, and outlet‐glacier thicknesses, snowfall rates on ice sheets, fjord and sub‐ice‐shelf bathymetry, and ice‐sheet hydrology. Unanticipated discoveries included a reliable method for constraining the thickness within difficult‐to‐sound incised troughs beneath ice sheets, the extent of the firn aquifer within the Greenland Ice Sheet, the vulnerability of many Greenland and Antarctic outlet glaciers to ocean‐driven melting at their grounding zones, and the dominance of surface‐melt‐driven mass loss of Alaskan glaciers. For sea ice, OIB significantly advanced our understanding of spatiotemporal variability in sea ice freeboard and its snow cover, especially through combined analysis of fine‐resolution altimetry, visible imagery, and snow radar measurements of the overlying snow thickness. Such analyses led to the unanticipated discovery of an interdecadal decrease in snow thickness on Arctic sea ice and numerous opportunities to validate sea ice freeboards from satellite radar altimetry. While many of its data sets have yet to be fully explored, OIB’s scientific legacy has already demonstrated the value of sustained investment in reliable airborne platforms, airborne instrument development, interagency and international collaboration, and open and rapid data access to advance our understanding of Earth’s remote polar regions and their role in the Earth system.

Abstract. Radar sounding of the thickness of temperate glaciers is more challenging than for polar ice sheets, due to the former's greater volume scattering (englacial water), surface scattering (crevasses and debris) and dielectric attenuation rate (warmer ice). Lower frequency (~1–100 MHz) radar sounders are commonly deployed to mitigate these effects, but the lack of a synthesis of existing radar-sounding surveys of temperate glaciers limits progress in system and survey design. Here we use a recent global synthesis of measured glacier thickness to evaluate the relation between the radar center frequency and maximum thickness. From a maximum reported thickness of ~1500 m near 1 MHz, the maximum thickness sounded decreases with increasing frequency by ~500 m per frequency decade. Newer airborne radar sounders generally outperform older, ground-based ones at comparable frequencies, so radar-sounder success is also influenced by system design and processing methods. Based on globally modeled glacier thicknesses, we conclude that a multi-element airborne radar sounder with a center frequency of ≤ 30 MHz could survey most temperate glaciers more efficiently than presently available systems.

We use repeat high-resolution airborne geophysical data consisting of laser altimetry, snow, and Ku-band radar and optical imagery acquired in 2014, 2016, and 2017 to analyze the spatial and temporal variability in surface roughness, slope, wind deposition, and snow accumulation at 88∘ S, an elevation bias validation site for ICESat-2 and potential validation site for CryoSat-2. We find significant small-scale variability (<10 km) in snow accumulation based on the snow radar subsurface stratigraphy, indicating areas of strong wind redistribution are prevalent at 88∘ S. In general, highs in snow accumulation rate correspond with topographic lows, resulting in a negative correlation coefficient of r2=-0.32 between accumulation rate and MSWD (mean slope in the mean wind direction). This relationship is strongest in areas where the dominant wind direction is parallel to the survey profile, which is expected as the geophysical surveys only capture a two-dimensional cross section of snow redistribution. Variability in snow accumulation appears to correlate with variability in MSWD. The correlation coefficient between the standard deviations of accumulation rate and MSWD is r2=0.48, indicating a stronger link between the standard deviations than the actual parameters. Our analysis shows that there is no simple relationship between surface slope, wind direction, and snow accumulation rates for the overall survey area. We find high variability in surface roughness derived from laser altimetry measurements on length scales smaller than 10 km, sometimes with very distinct and sharp transitions. Some areas also show significant temporal variability over the course of the 3 survey years. Ultimately, there is no statistically significant slope-independent relationship between surface roughness and accumulation rates within our survey area. The observed correspondence between the small-scale temporal and spatial variability in surface roughness and backscatter, as evidenced by Ku-band radar signal strength retrievals, will make it difficult to develop elevation bias corrections for radar altimeter retrieval algorithms.

Surface height and total freeboard from the Ice, Cloud, and Land Elevation Satellite‐2 (ICESat‐2, IS‐2) sea ice data products (ATL07/ATL10) are assessed with near‐coincident retrievals from the Airborne Topographic Mapper (ATM) lidar in four dedicated underflights during the 2019 Operation IceBridge Arctic deployment. Over a mix of seasonal and older ice, we find remarkable correlations between the ATM and IS‐2 height profiles and roughness (in ninety‐nine 10‐km segments) that averages to >0.95 and > 0.97, respectively. Regression slopes near unity, between 0.93 and 0.99, indicate close agreement of the height estimates. Larger differences between the surface heights are seen in rougher areas where it is more difficult for the photon heights (used in IS‐2 surface finding) to capture the surface distributions at short length scales. Total freeboard in 10‐km segments, calculated using three different approaches, show variability of 0.02 to 0.04 m. Sources of residual variance, attributable to differences between the two instruments, are discussed.

Following the discovery of the Hiawatha impact crater beneath the northwest margin of the Greenland Ice Sheet, we explored satellite and aerogeophysical data in search of additional such craters. Here we report the discovery of a possible second subglacial impact crater that is 36.5‐km wide and 183 km southeast of the Hiawatha impact crater. Although buried by 2 km of ice, the structure's rim induces a conspicuously circular surface expression, it possesses a central uplift, and it causes a negative gravity anomaly. The existence of two closely spaced and similarly sized complex craters raises the possibility that they formed during related impact events. However, the second structure's morphology is shallower, its overlying ice is conformal and older, and such an event can be explained by chance. We conclude that the identified structure is very likely an impact crater, but it is unlikely to be a twin of the Hiawatha impact crater.

Using field observations followed by petrological, geochemical, geochronological, and geophysical data we infer the presence of a previously unknown Miocene subglacial volcanic center ~230 km from the South Pole. Evidence of volcanism is from boulders of olivine-bearing amygdaloidal/vesicular basalt and hyaloclastite deposited in a moraine in the southern Transantarctic Mountains. 40Ar/39Ar ages from five specimens plus U-Pb ages of detrital zircon from glacial till indicate igneous activity 25-17 Ma. The likely source of the volcanism is a circular -735 nT magnetic anomaly 60 km upflow from the sampling site. Subaqueous textures of the volcanics indicate eruption beneath ice or into water at the margin of an ice mass during the early Miocene. These rocks record the southernmost Cenozoic volcanism in Antarctica and expand the known extent of the oldest lavas associated with West Antarctic rift system. They may be an expression of lithospheric foundering beneath the southern Transantarctic Mountains.