Mountain glaciers provide an important service in sustaining river flows for large populations downstream of High Mountain Asia (HMA) but these glaciers are retreating and the future of this water resource is highly uncertain. Glacier thickness measurements are vital for accurate mapping of the remaining ice reserve and for predicting where and how fast it will decline under climate change, but such measurements are severely lacking in this region due to the difficulties of surveying in remote, high-altitude settings. We report on a uniquely extensive new thickness dataset for eleven glaciers in the Khumbu Himal around Mount Everest that we collected in late 2019 using a novel, low-frequency helicopter-borne radar. To aid in interpreting the survey radargrams we developed a terrain clutter model, and we succeeded in mapping ice thickness with a precision of around +/- 7 % and horizontal spacing of around 40 m, for thicknesses of up to 445 m and spanning a total of 119 line-km, approximately doubling the length of previous thickness surveys in HMA. To demonstrate the utility of our new measurements, we compare them to existing modelled thickness products and find that the models struggle to reproduce the distribution of ice in these complex, steep, rapidly slowing, thinning and stagnating glaciers, with widespread systematic thin and thick biases equivalent to around half of the measured ice thickness or more. This new dataset (10.5285/e39647f5-fb72-4d16-acbd-9784ed2167b8 (Pritchard al., 2025) permits for the first time a detailed analysis of model performance on Himalayan glaciers, a key step in improving model skill and hence the accuracy of modelled thickness distributions and future ice loss on the mountain-range scale.

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.

The accuracy of sea-level rise predictions is highly dependent on reliably understanding the subglacial environment beneath ice streams. Bedforms result from the interaction between ice and its substrate, and therefore have the potential to shed light on poorly understood basal conditions influencing ice dynamics. However, theoretical models of bedform formation are often based on observations from deglaciated areas or from sparsely-sampled geophysical surveys over glaciated regions. Here, we use high-resolution three- and two-dimensional radar and seismic data to reveal details of the initiation and evolution of a subglacial bedform beneath Rutford Ice Stream, West Antarctica. Radar surveys at 20 m and 50 m line spacing allow detailed imaging of bed topography, including a moat up to 55 m deep, surrounding the upstream end of a 50-m high and >18-km long bedform. Many models rely on either a topographical or a locally resistant seed point to initiate bedform formation. The bedform described here is mostly composed of soft sediment (porosity >0.3) and lacks a boulder or outcrop, suggesting the bedform initiated without a topographic seed point. Sediment at the upstream end of bedforms appears stiffer over a distance of 2.3 km. We suggest sediment inhomogeneities in the initially flat bed cause the deposition of sediment, which, assuming resistant enough, acts as a seed point for bedform extension and moat erosion. The moat's geometry and its truncation of other bedforms suggest that it was eroded after the deposition of surrounding bedforms. These observations from a modern ice stream deliver information of subglacial processes involved in the initiation as well as in situ high-resolution topography and properties of bedforms and moats. Using these observations numerical models can be tested and developed accordingly.

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.

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.

Basal properties beneath ice streams and glaciers are known to be a control for ice flow dynamics, hence knowledge of them is crucial for predicting sea level due to changes in glacial dynamics. Basal properties, processes and topography also drive the formation of subglacial landforms. Bed properties beneath Rutford Ice Stream (West Antarctica) have previously been described using seismic acoustic impedance measurements at a sparse spatial coverage. Here, we derive bed properties in a 15 × 17 km grid of surface radar data with coverage and sampling much higher than previous seismic studies. Bed reflection amplitudes in surface radar data were calibrated using sediment porosities (ranging from 0.4–0.5) derived from seismic acoustic impedance. We find the bed properties are spatially variable, consisting of low porosity material in some areas and soft sediment in other areas. Comparison of seismic and surface radar data imply the low porosity material to be a consolidated sediment or sedimentary rock. Mega‐scale glacial lineations (MSGLs) are ubiquitous on the bed and consist of soft, high porosity, probably deforming sediment, consistent with previous interpretations of MSGLs. We find some MSGLs have high reflectivity on their crest, interpreted as water bodies overlying high porosity sediment, whereas the trough around and the upstream end of some landforms consist of low porosity material. Integrating these different observations, we place constraints on possible explanations for the occurrence of water on the crest of landforms.

Subglacial bedforms such as mega-scale glacial lineations and drumlins are commonly thought to form during active ice flow. They are often present in deglaciated areas with various elongation ratios, consisting of different materials , information which led to the development of different formation theories. However, these exposed examples were formed not only by subglacial processes during glaciation but also altered by processes during and after deglaciation. Here, we analyse in-situ properties and topography beneath Rutford Ice Steam, a fast flowing ice stream in West Antarctica to evaluate current theoretical models of bedform formation. We present a combination of seismic and radar data, including high-resolution 3D radar topography covering the upstream end of a bedform. Data acquisition and processing of the high-resolution 3D radar dataset result in a horizontal resolution of 24 m along- and across-track and a vertical resolution of 12 m. Using seismic acoustic impedance and calibrated radar reflectivity subglacial properties of the bedforms as well as the surrounding area are identified. A depression around the upstream end of a 360 m wide, 50 m high and more than 13 km long bedform was observed for the first time analysing the high-resolution 3D radar topography. The depression consists of a deepening up to 45 m deep and 360 m wide and is seen to extend around 10.5 km downstream. Radar reflectivity reveals that the material the depression is excavated into at least partly consists of low porosity material. Radar reflectivity and seismic acoustic impedance along the bedform imply a stiffer upstream end which softens along flow. The subglacial topography and properties give evidence that the bedform and the depression are formed by a combination of erosional and depositional processes. Both processes are likely interlinked, as implied by the comparable volume of the moat and the bedform at the upstream end of the bedform. Based on these observations we support or reject common bedform formation theories beneath Rutford Ice Stream.

Earth and Space Science Open Archive This preprint has been submitted to and is under consideration at Journal of Geophysical Research - Earth Surface. ESSOAr is a venue for early communication or feedback before peer review. Data may be preliminary.Learn more about preprints preprintOpen AccessYou are viewing an older version [v1]Go to new versionSpatial Variation in Radar-Derived Basal Properties of Rutford Ice Stream, West AntarcticaAuthorsRebeccaSchlegeliDTaviMurrayiDAndrew MarkSmithiDAlex MarkBrisbourneiDAdam DBoothiDEdward CharlesKingRoger ArthurClarkiDSee all authors Rebecca SchlegeliDCorresponding Author• Submitting AuthorSwansea UniversityiDhttps://orcid.org/0000-0003-1149-2816view email addressThe email was not providedcopy email addressTavi MurrayiDSwansea UniversityiDhttps://orcid.org/0000-0001-6714-6512view email addressThe email was not providedcopy email addressAndrew Mark SmithiDBritish Antarctic SurveyiDhttps://orcid.org/0000-0001-8577-482Xview email addressThe email was not providedcopy email addressAlex Mark BrisbourneiDBritish Antarctic SurveyiDhttps://orcid.org/0000-0002-9887-7120view email addressThe email was not providedcopy email addressAdam D BoothiDUniversity of LeedsiDhttps://orcid.org/0000-0002-8166-9608view email addressThe email was not providedcopy email addressEdward Charles KingBritish Antarctic Surveyview email addressThe email was not providedcopy email addressRoger Arthur ClarkiDUniversity of LeedsiDhttps://orcid.org/0000-0002-3508-8183view email addressThe email was not providedcopy email address

Recent advances in the measurement of bedforms beneath active ice streams have been made using ground-based grid profiling using impulse radar systems operating with centre frequencies in the 3-5 MHz range. Surveys of Rutford Ice Stream and Pine Island Glacier have shown that features such as mega-scale glacial lineations with topographic relief of as little as 3 m can be traced for many kilometres downstream under more than 2 km of fast-moving ice. In the discussion of these data, it is often asked 'How is it possible to map such fine-scale topography with such a low-frequency radar'. In answering that question, the key point is the distinction between the precision of a radar range measurement to a single, isolated reflective interface and the ability to resolve the presence of two closely-spaced interfaces of similar reflectivity (commonly referred to as the vertical resolution). This paper will discuss and illustrate this distinction and use the case study of data acquired over Pine Island Glacier to examine the limits of precision of the radar range measurement.

The thickness of glaciers in High-Mountain Asia (HMA) is critical in determining when the ice reserve will be lost as these glaciers thin but is remarkably poorly known because very few measurements have been made. Through a series of ground-based and airborne field tests, we have adapted a low-frequency ice-penetrating radar developed originally for Antarctic over-snow surveys, for deployment as a helicopter-borne system to increase the number of measurements. The manoeuvrability provided by helicopters and the ability of our system to detect glacier beds through thick, dirty, temperate ice makes it well suited to increase greatly the sample of measurements available for calibrating ice thickness models on the regional and global scale. The Bedmap Himalayas radar-survey system can reduce the uncertainty in present-day ice volumes and therefore in projections of when HMA's river catchments will lose this hydrological buffer against drought.

Suture zones are abundant on Antarctic ice shelves and widely observed to impede fracture propagation, greatly enhancing ice-shelf stability. Using seismic and radar observations on the Larsen C Ice Shelf of the Antarctic Peninsula, we confirm that such zones are highly heterogeneous, consisting of multiple meteoric and marine ice bodies of diverse provenance fused together. Here we demonstrate that fracture detainment is predominantly controlled by enhanced seawater content in suture zones, rather than by enhanced temperature as previously thought. We show that interstitial seawater can reduce fracture-driving stress by orders of magnitude, promoting both viscous relaxation and the development of micro cracks, the incidence of which scales inversely with stress intensity. We show how simple analysis of viscous buckles in ice-penetrating radar data can quantify the seawater content of suture zones and their modification of the ice-shelf's stress regime. By limiting fracture, enhancing stability and restraining continental ice discharge into the ocean, suture zones act as vital regulators of Antarctic mass balance.

Projections of sea-level rise contributions from West Antarctica's dynamically thinning ice streams contain high uncertainty because some of the key processes involved are extremely challenging to observe. An especially poorly observed parameter is sub-decadal stability of ice-stream beds, which may be important for subglacial traction, till continuity and landform development. Only two previous studies have made repeated geophysical measurements of ice-stream beds at the same locations in different years, but both studies were limited in spatial extent. Here, we present the results from repeat radar measurements of the bed of Pine Island Glacier, West Antarctica, conducted 3–6 years apart, along a cumulative ∼ 60 km of profiles. Analysis of the correlation of bed picks between repeat surveys shows that 90 % of the bed displays no significant change despite the glacier increasing in speed by up to 40 % over the last decade. We attribute the negligible detection of morphological change at the bed of Pine Island Glacier to the ubiquitous presence of a deforming till layer, wherein sediment transport is in steady state, such that sediment is transported along the basal interface without inducing morphological change to the radar-sounded basal interface. Given the precision of our measurements, the upper limit of subglacial erosion observed here is 500 mm a−1, far exceeding erosion rates reported for glacial settings from proglacial sediment yields, but substantially below subglacial erosion rates of 1.0 m a−1 previously reported from repeat geophysical surveys in West Antarctica.

Flowlines on ice sheets and glaciers form complex patterns. To explore their role in ice routing and extend the language for studying such patterns, we develop a theory of flow convergence and curvature in plan view. These geometric quantities respectively equal the negative divergence of the vector field of ice-flow direction and the curl of this field. From the first of these two fundamental results, we show that flow in individual catchments of an ice sheet can converge (despite its overall spreading) because ice divides are loci of strong divergence, and that a sign bifurcation in convergence occurs during ice-sheet "symmetry breaking" (the transition from near-radial spreading to spreading with substantial azimuthal velocities) and during the formation of ice-stream tributary networks. We also uncover the topological control behind balance-flux distributions across ice masses. Notably, convergence participates in mass conservation along flowlines to amplify ice flux via a positive feedback; thus the convergence field governs the form of ice-stream networks simulated by balance-velocity models. The theory provides a roadmap for understanding the tower-shaped plot of flow speed versus convergence for the Antarctic Ice Sheet.

We report on the recent reactivation of a large rift in the Brunt Ice Shelf, East Antarctica, in December 2012 and the formation of a 50 km long new rift in October 2016. Observations from a suite of ground-based and remote sensing instruments between January 2000 and July 2017 were used to track progress of both rifts in unprecedented detail. Results reveal a steady accelerating trend in their width, in combination with alternating episodes of fast ( > 600 m day−1) and slow propagation of the rift tip, controlled by the heterogeneous structure of the ice shelf. A numerical ice flow model and a simple propagation algorithm based on the stress distribution in the ice shelf were successfully used to hindcast the observed trajectories and to simulate future rift progression under different assumptions. Results show a high likelihood of ice loss at the McDonald Ice Rumples, the only pinning point of the ice shelf. The nascent iceberg calving and associated reduction in pinning of the Brunt Ice Shelf may provide a uniquely monitored natural experiment of ice shelf variability and provoke a deeper understanding of similar processes elsewhere in Antarctica.

Lakes beneath the Antarctic Ice Sheet are of fundamental scientific interest for their ability to contain unique records of ice sheet history and microbial life in their sediments. However, no records of subglacial lake sedimentation have yet been acquired from beneath the interior of the ice sheet, and understanding of sediment pathways, processes and structure in subglacial lake environments remains uncertain. Here we present an analysis of seismic data from Subglacial Lake Ellsworth, showing that the lake bed comprises very fine-grained sediments deposited in a low energy environment, with low water- and sediment-fluxes. Minimum sediment thickness is 6 m, the result of prolonged low sedimentation rates. Based on the few available analogues, we speculate this sediment age range is a minimum of 150 ka, and possibly >1 Ma. Sediment mass movements have occurred, but they are rare and have been buried by subsequent sedimentation. We present a new conceptual model of subglacial lake sedimentation, allowing a framework for evaluating processes in subglacial lake environments, and for determining future lake access locations and interpreting subglacial lake samples.

The aminostratigraphy of up to 1000m of glacial and interglacial sediments in the North Sea Basin (NSB) is compiled from multiple boreholes sites and dated using strontium isotope (Sr) analysis to provide a chronological framework extending throughout the Quaternary. The Quaternary aminostratigraphy within the NSB can be divided into five amino zones (AZ), AZ1 (youngest) to AZ5 based on alle/Ile ratios which are used to correlate depositional events recorded in sedimentary facies from multiple borehole sites. The new relative and absolute chronology also ties to seismostratigraphy and multi-generational tunnel valley fill. Early Pleistocene sediments produce AZ5 ratios with Sr ages around 1.4–1.5Ma deposited in marine to glacimarine conditions. Sediments with AZ4 ratios are constrained by Sr ages and magnetic polarity to 1.3–0.9Ma, spanning the eccentricity to axial tilt insolation periodicity change of the Mid Pleistocene Transition (MPT). Some coarse-grained sedimentary units with AZ4 ratios may be linked to ice proximal deposition during the MPT. Sediments with AZ3 ratios are poorly chronologically constrained to between 1.1 and 0.26Ma. At least one episode of diminished coarse clasts and increased foraminiferal diversity during this Mid Pleistocene time suggests a progression from ice proximal to warmer, possibly interglacial conditions. Furthermore, the stratigraphically deepest generation of tunnel valleys cut strata with AZ3 ratios indicating subglacial conditions at least once during the Mid Pleistocene. Sediments with AZ2 ratios corresponding to MIS7-5 are tied to seismic units including a thick (up to 40m) deformation till, suggesting glacial incursion into the NSB during MIS6 followed by tunnel valley cutting of the this till, thought to be an MIS6 event also. The geometry and cross-cutting complexities of such tunnel valleys precludes further correlation of this ice extent. The youngest aminostratigraphic zone, AZ1, is recognised in all boreholes and well constrained to MIS4-2 and the youngest set of tunnel valleys, attributed to a single or multiple glacial advances at this time.