Explore the vast range of titles available.

Radar sounding is a powerful geophysical approach for characterizing the subsurface conditions of terrestrial and planetary ice masses at local to global scales. As a result, a wide array of orbital, airborne, ground-based, and in situ instruments, platforms and data analysis approaches for radioglaciology have been developed, applied or proposed. Terrestrially, airborne radar sounding has been used in glaciology to observe ice thickness, basal topography and englacial layers for five decades. More recently, radar sounding data have also been exploited to estimate the extent and configuration of subglacial water, the geometry of subglacial bedforms and the subglacial and englacial thermal states of ice sheets. Planetary radar sounders have observed, or are planned to observe, the subsurfaces and near-surfaces of Mars, Earth's Moon, comets and the icy moons of Jupiter. In this review paper, and the thematic issue of the Annals of Glaciology on ‘Five decades of radioglaciology’ to which it belongs, we present recent advances in the fields of radar systems, missions, signal processing, data analysis, modeling and scientific interpretation. Our review presents progress in these fields since the last radio-glaciological Annals of Glaciology issue of 2014, the context of their history and future prospects.

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 subglacial lake, Lake Snow Eagle (LSE), was identified in East Antarctica by airborne geophysical surveys. LSE, contained within a subglacial canyon, likely hosts a valuable sediment record of the geological and glaciological changes of interior East Antarctica. Understanding past lake activity is crucial for interpreting this record. Here, we present the englacial radiostratigraphy in the LSE area mapped by airborne ice-penetrating radar, which reveals a localized high-amplitude variation in ice unit thickness that is estimated to be ∼12 ka old. Using an ice-flow model that simulates englacial stratigraphy, we investigate the origin of this feature and its relationship to changes in ice dynamical boundary conditions. Our results reveal that local snowfall redistribution initiated around the early Holocene is likely the primary cause, resulting from a short-wavelength (∼10 km) high-amplitude (∼20 m) ice surface slope variation caused by basal lubrication over a large subglacial lake. This finding indicates an increase in LSE water volume during the Holocene, illustrating the sensitivity in volume of a major topographically constrained subglacial lake across a single glacial cycle. This study demonstrates how englacial stratigraphy can provide valuable insight into subglacial hydrological changes before modern satellite observations, both for LSE and potentially at other locations.

Viscous flow features known as glacier-like forms on Mars have been observed emerging from alcoves that resemble cirques on Earth. However, many alcoves exist without associated glacier-like forms, and these features have never been studied or categorized at a regional population scale. On Earth, cirques form when depressions on mountain slopes accumulate snow, which gradually compacts into glacial ice. As the glacier flows downhill, it deepens the depression through erosion. Most of this erosion is driven by wet-based glaciers, although cold-based glaciers can also contribute to some headward and sidewall retreat. Here, we present evidence that cirque-like alcoves on Mars, similar to terrestrial cirques, are shaped by glacial erosion. To assess which alcoves on Mars are most “cirque-like”, we mapped a regional population of ∼ 2000 alcoves in Deuteronilus Mensae (40–48° N, 16–35° E), a region in the mid-latitudes of Mars characterized by mesas surrounded by glacial remnants. Based on visual characteristics and morphometrics, we refined our dataset to 434 “cirque-like alcoves” – nearly six times the amount of glacier-like forms in the region – and used this to better understand the potential contribution of glaciation to landscape evolution in Deuteronilus Mensae, Mars. High-resolution imagery reveals geomorphic evidence for past glacial occupation of these cirque-like alcoves, including flow features, linear terrain, moraine-like ridges, mound-and-tail terrain, polygonal terrain, rectilinear-ridge terrain, and washboard terrain, as well as an ice-rich mantling unit. Most cirque-like alcoves face south to southeast, similar to gullies poleward of 40°. One possibility to explain this trend is that southward facing cirque-like alcoves in the northern mid-latitudes formed when conditions were more favorable for ice accumulation during periods of high obliquity. Using wet-based glacial erosion rates, assuming a mean annual temperature of 0 °C (compared to present-day temperatures of ∼-63 °C), the timescales for Martian cirque-like alcove formation align with both the estimated ages of glacier-like forms (millions to tens of millions of years) and other viscous flow features such as lobate debris aprons (hundreds of millions of years). In contrast, using a temperature assumption of −50 to −68 °C, cold-based erosion rates are only consistent with the older ages of lobate debris aprons. By mapping cirque-like alcoves at a large scale for the first time, we expand the catalog of features attributed to glacial erosion on Mars. Future work is needed to specify the timing of the formation of cirque-like alcoves and whether their formation required warm-based erosion.

Radio-echo sounding (RES) has revealed an internal architecture within both the West and East Antarctic ice sheets that records their depositional, deformational and melting histories. Crucially, RES-imaged internal-reflecting horizons, tied to ice-core age–depth profiles, can be treated as isochrones that record the age–depth structure across the Antarctic ice sheets. These enable the reconstruction of past climate and ice dynamical processes on large scales, which are complementary to but more spatially extensive than commonly used proxy records (e.g. former ice limits constrained by cosmogenic dating or offshore sediment sequences) around Antarctica. We review the progress towards building a pan-Antarctic age–depth model from these data by first introducing the relevant RES datasets that have been acquired across Antarctica over the last 6 decades (focussing specifically on those that detected internal-reflecting horizons) and outlining the processing steps typically undertaken to visualise, trace and date (by intersection with ice cores or modelling) the RES-imaged isochrones. We summarise the scientific applications for which Antarctica's internal architecture has been used to date and present a pathway to expanding Antarctic radiostratigraphy across the continent to provide a benchmark for a wider range of investigations: (1) identification of optimal sites for retrieving new ice-core palaeoclimate records targeting different periods; (2) reconstruction of surface mass balance on millennial or historical timescales; (3) estimation of basal melting and geothermal heat flux from radiostratigraphy and comprehensive mapping of basal-ice units to complement inferences from other geophysical and geological methods; (4) advancement of the knowledge of volcanic activity and fallout across Antarctica; and (5) refinement of numerical models that leverage radiostratigraphy to tune time-varying accumulation, basal melting and ice flow, firstly to reconstruct past behaviour and then to reduce uncertainties in projecting future ice-sheet behaviour.

Spatial and temporal variations in past accumulation, ice thickness, and ice flow of polar ice sheets are weakly constrained on Earth, and are fundamental unknowns on Mars. On Earth, the spatial and temporal histories of accumulation and ice-sheet flow are necessary to recreate ice-volume and sea-level histories, and are important to properly interpret ice-core chemistry. On Mars, accumulation and ice-flow histories are necessary to decipher the connection between climate and ice-mass formation, evolution, and observable structure. Internal layers in ice sheets on Earth and on Mars have been observed with ice-penetrating radar. These layers preserve information about how the ice sheet responded to past spatial and/or temporal changes in accumulation rate and ice flow, and present-day internal-layer shapes observed by radar are the most accessible remaining record of this past information. Deeper layers contain information from further in the past, making them highly valuable, but they are more difficult to decipher. In this work, an inverse problem is solved to infer transients in accumulation rate, ice-sheet thickness, and ice flow from the shapes of deep internal layers. While some details of these histories can be recovered from ice cores, ice cores represent conditions at only a single point. However, the approach presented here is more robust in combination with ice-core data. If internal layers are dated, for example by an intersecting ice core, then radar-observed internal layers provide both spatial and temporal information. Each layer represents a past surface of a particular age that has been subsequently buried by accumulation and also modified by ice flow. In this work, the goal of solving this inverse problem is to find a set of model parameters (e.g. accumulation-rate history) that have the minimum variation required to explain the data (e.g. internal-layer shapes). The process of internal-layer formation is described with a 2.5-D thermomechanical ice-flow flowband model. Estimates of the data are matched to measured values within their uncertainties, and to an expected tolerance. We seek an accumulation pattern that is spatially smooth, and a parameter set that is consistent with characteristic values of the parameters. This dissertation presents this inverse approach, and discusses applications to data from Antarctica and from Mars.

Chronologies of glacier deposits in the Transantarctic Mountains provide important constraints on grounding-line retreat during the last deglaciation in the Ross Sea. However, between Beardmore Glacier and Ross Island – a distance of some 600 km – the existing chronologies are generally sparse and far from the modern grounding line, leaving the past dynamics of this vast region largely unconstrained. We present exposure ages of glacial deposits at three locations alongside the Darwin–Hatherton Glacier System – including within 10 km of the modern grounding line – that record several hundred meters of Late Pleistocene to Early Holocene thickening relative to present. As the ice sheet grounding line in the Ross Sea retreated, Hatherton Glacier thinned steadily from about 9 until about 3 ka. Our data are equivocal about the maximum thickness and Mid-Holocene to Early Holocene history at the mouth of Darwin Glacier, allowing for two conflicting deglaciation scenarios: (1) ∼500 m of thinning from 9 to 3 ka, similar to Hatherton Glacier, or (2) ∼950 m of thinning, with a rapid pulse of ∼600 m thinning at around 5 ka. We test these two scenarios using a 1.5-dimensional flowband model, forced by ice thickness changes at the mouth of Darwin Glacier and evaluated by fit to the chronology of deposits at Hatherton Glacier. The constraints from Hatherton Glacier are consistent with the interpretation that the mouth of Darwin Glacier thinned steadily by ∼500 m from 9 to 3 ka. Rapid pulses of thinning at the mouth of Darwin Glacier are ruled out by the data at Hatherton Glacier. This contrasts with some of the available records from the mouths of other outlet glaciers in the Transantarctic Mountains, many of which thinned by hundreds of meters over roughly a 1000-year period in the Early Holocene. The deglaciation histories of Darwin and Hatherton glaciers are best matched by a steady decrease in catchment area through the Holocene, suggesting that Byrd and/or Mulock glaciers may have captured roughly half of the catchment area of Darwin and Hatherton glaciers during the last deglaciation. An ensemble of three-dimensional ice sheet model simulations suggest that Darwin and Hatherton glaciers are strongly buttressed by convergent flow with ice from neighboring Byrd and Mulock glaciers, and by lateral drag past Minna Bluff, which could have led to a pattern of retreat distinct from other glaciers throughout the Transantarctic Mountains.

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.