Duncan A. Young

Radar Sounding Evidence for Buried Glaciers in the Southern Mid-Latitudes of Mars

Lobate features abutting massifs and escarpments in the middle latitudes of Mars have been recognized in images for decades, but their true nature has been controversial, with hypotheses of origin such as ice-lubricated debris flows or glaciers covered by a layer of surface debris. These models imply an ice content ranging from minor and interstitial to massive and relatively pure. Soundings of these deposits in the eastern Hellas region by the Shallow Radar on the Mars Reconnaissance Orbiter reveal radar properties entirely consistent with massive water ice, supporting the debris-covered glacier hypothesis. The results imply that these glaciers formed in a previous climate conducive to glaciation at middle latitudes. Such features may collectively represent the most extensive nonpolar ice yet recognized on Mars.

New Boundary Conditions for the West Antarctic Ice Sheet: Subglacial Topography of the Thwaites and Smith Glacier Catchments

Airborne radar sounding over the Thwaites Glacier (TG) catchment and its surroundings provides the first comprehensive view of subglacial topography in this dynamic part of the West Antarctic Ice Sheet (WAIS) and reveals that TG is underlain by a single, broad basin fed by a dendritic pattern of valleys, while Smith Glacier lies within an extremely deep, narrow trench. Subglacial topography in the TG catchment slopes inland from a broad, low-relief coastal sill to the thickest ice of the WAIS and makes deep connections to both Pine Island Glacier and the Ross Sea Embayment enabling dynamic interactions across the WAIS during deglaciation. Simple isostatic rebound modeling shows that most of this landscape would be submarine after deglaciation, aside from an island chain near the present-day Ross-Amundsen ice divide. The lack of topographic confinement along TGs eastern margin implies that it may continue to widen in response to grounding line retreat.

A dynamic early East Antarctic Ice Sheet suggested by ice-covered fjord landscapes

The first Cenozoic ice sheets initiated in Antarctica from the Gamburtsev Subglacial Mountains and other highlands as a result of rapid global cooling ∼34 million years ago. In the subsequent 20 million years, at a time of declining atmospheric carbon dioxide concentrations and an evolving Antarctic circumpolar current, sedimentary sequence interpretation and numerical modelling suggest that cyclical periods of ice-sheet expansion to the continental margin, followed by retreat to the subglacial highlands, occurred up to thirty times. These fluctuations were paced by orbital changes and were a major influence on global sea levels. Ice-sheet models show that the nature of such oscillations is critically dependent on the pattern and extent of Antarctic topographic lowlands. Here we show that the basal topography of the Aurora Subglacial Basin of East Antarctica, at present overlain by 2–4.5 km of ice, is characterized by a series of well-defined topographic channels within a mountain block landscape. The identification of this fjord landscape, based on new data from ice-penetrating radar, provides an improved understanding of the topography of the Aurora Subglacial Basin and its surroundings, and reveals a complex surface sculpted by a succession of ice-sheet configurations substantially different from today’s. At different stages during its fluctuations, the edge of the East Antarctic Ice Sheet lay pinned along the margins of the Aurora Subglacial Basin, the upland boundaries of which are currently above sea level and the deepest parts of which are more than 1 km below sea level. Although the timing of the channel incision remains uncertain, our results suggest that the fjord landscape was carved by at least two iceflow regimes of different scales and directions, each of which would have over-deepened existing topographic depressions, reversing valley floor slopes.

Along-Track Focusing of Airborne Radar Sounding Data From West Antarctica for Improving Basal Reflection Analysis and Layer Detection

This paper presents focused synthetic aperture radar (SAR) processing of airborne radar sounding data acquired with the High-Capability Radar Sounder system at 60 MHz. The motivation is to improve basal reflection analysis for water detection and to improve layer detection and tracking. The processing and reflection analyses are applied to data from Kamb Ice Stream, West Antarctica. The SAR processor correlates the radar data with reference echoes from subsurface point targets. The references are 1-D responses limited by the pulse nadir footprint or 2-D responses that include echo tails. Unfocused SAR and incoherent integration are included for comparison. Echoes are accurately preserved from along-track slopes up to about 0.5deg for unfocused SAR, 3deg for 1-D correlations, and 10deg for 2-D correlations. The noise/clutter levels increase from unfocused SAR to 1-D and 2-D correlations, but additional gain compensates at the basal interface. The basal echo signal-to-noise ratio improvement is typically about 5 dB, and up to 10 dB for 2-D correlations in rough regions. The increased noise degrades the clarity of internal layers in the 2-D correlations, but detection of layers with slopes greater than 3deg is improved. Reflection coefficients are computed for basal water detection, and the results are compared for the different processing methods. There is a significant increase in the detected water from unfocused SAR to 1-D correlations, indicating that substantial basal water exists on moderately sloped interfaces. Very little additional water is detected from the 2-D correlations. The results from incoherent integration are close to the focused SAR results, but the noise/clutter levels are much greater.

Evidence for elevated and spatially variable geothermal flux beneath the West Antarctic Ice Sheet

Significance Thwaites Glacier is one of the West Antarcticas most prominent, rapidly evolving, and potentially unstable contributors to global sea level rise. Uncertainty in the amount and spatial pattern of geothermal flux and melting beneath this glacier is a major limitation in predicting its future behavior and sea level contribution. In this paper, a combination of radar sounding and subglacial water routing is used to show that large areas at the base of Thwaites Glacier are actively melting in response to geothermal flux consistent with rift-associated magma migration and volcanism. This supports the hypothesis that heterogeneous geothermal flux and local magmatic processes could be critical factors in determining the future behavior of the West Antarctic Ice Sheet. Heterogeneous hydrologic, lithologic, and geologic basal boundary conditions can exert strong control on the evolution, stability, and sea level contribution of marine ice sheets. Geothermal flux is one of the most dynamically critical ice sheet boundary conditions but is extremely difficult to constrain at the scale required to understand and predict the behavior of rapidly changing glaciers. This lack of observational constraint on geothermal flux is particularly problematic for the glacier catchments of the West Antarctic Ice Sheet within the low topography of the West Antarctic Rift System where geothermal fluxes are expected to be high, heterogeneous, and possibly transient. We use airborne radar sounding data with a subglacial water routing model to estimate the distribution of basal melting and geothermal flux beneath Thwaites Glacier, West Antarctica. We show that the Thwaites Glacier catchment has a minimum average geothermal flux of ∼114 ± 10 mW/m2 with areas of high flux exceeding 200 mW/m2 consistent with hypothesized rift-associated magmatic migration and volcanism. These areas of highest geothermal flux include the westernmost tributary of Thwaites Glacier adjacent to the subaerial Mount Takahe volcano and the upper reaches of the central tributary near the West Antarctic Ice Sheet Divide ice core drilling site.

Evidence for a water system transition beneath Thwaites Glacier, West Antarctica

Thwaites Glacier is one of the largest, most rapidly changing glaciers on Earth, and its landward-sloping bed reaches the interior of the marine West Antarctic Ice Sheet, which impounds enough ice to yield meters of sea-level rise. Marine ice sheets with landward-sloping beds have a potentially unstable configuration in which acceleration can initiate or modulate grounding-line retreat and ice loss. Subglacial water has been observed and theorized to accelerate the flow of overlying ice dependent on whether it is hydrologically distributed or concentrated. However, the subglacial water systems of Thwaites Glacier and their control on ice flow have not been characterized by geophysical analysis. The only practical means of observing these water systems is airborne ice-penetrating radar, but existing radar analysis approaches cannot discriminate between their dynamically critical states. We use the angular distribution of energy in radar bed echoes to characterize both the extent and hydrologic state of subglacial water systems across Thwaites Glacier. We validate this approach with radar imaging, showing that substantial water volumes are ponding in a system of distributed canals upstream of a bedrock ridge that is breached and bordered by a system of concentrated channels. The transition between these systems occurs with increasing surface slope, melt-water flux, and basal shear stress. This indicates a feedback between the subglacial water system and overlying ice dynamics, which raises the possibility that subglacial water could trigger or facilitate a grounding-line retreat in Thwaites Glacier capable of spreading into the interior of the West Antarctic Ice Sheet.

Evidence of a hydrological connection between the ice divide and ice sheet margin in the Aurora Subglacial Basin, East Antarctica

Subglacial hydrology in East Antarctica is poorly understood, yet may be critical to the manner in which ice flows. Data from a new regional airborne geophysical survey (ICECAP) have transformed our understanding of the topography and glaciology associated with the 287,000 km2 Aurora Subglacial Basin in East Antarctica. Using these data, in conjunction with numerical ice sheet modeling, we present a suite of analyses that demonstrate the potential of the 1000 km-long basin as a route for subglacial water drainage from the ice sheet interior to the ice sheet margin. We present results from our analysis of basal topography, bed roughness and radar power reflectance and from our modeling of ice sheet flow and basal ice temperatures. Although no clear-cut subglacial lakes are found within the Aurora Basin itself, dozens of lake-like reflectors are observed that, in conjunction with other results reported here, support the hypothesis that the basin acts as a pathway allowing discharge from subglacial lakes near the Dome C ice divide to reach the coast via the Totten Glacier.

Using radar-sounding data to identify the distribution and sources of subglacial water: application to Dome C, East Antarctica

Basal radar reflectivity is the most important measurement for the detection of subglacial water. However, dielectric loss in the overlying ice column complicates the determination of basal reflectivity. Dielectric attenuation is a function of ice temperature and impurity concentration. Temperature distribution is a function of climate history, basal heat flow and vertical strain rate, all of which can be partially inferred from the structure of dated internal layers. Using 11 dated layers, isotope records from Dome C, East Antarctica, and a model of the spatial variation of geothermal flux, we calculate the vertical strain rate and accumulation-rate history, allowing identification of areas where the basal melt rate exceeds 1.5 mm a -1 . The accumulation-rate history and vertical strain rates are then used as inputs for a transient temperature model. The model outputs for the present-day temperature distribution are then combined with depth-dependent ionic concentrations to model dielectric loss and infer basal reflectivity. The resulting reflection coefficients are consistent (� -5 dB) across a variety of subglacial water bodies. We also identify a high reflectivity >-15 dB in Concordia Trench and along suspected subglacial water-flow routes in Vincennes Basin. Highland areas tend to have highly variable reflection coefficients near -30 dB, consistent with an ice-bedrock interface. This combined model also identifies three areas of enhanced basal melting along Concordia Ridge, Concordia Subglacial Lake and Vincennes Basin. Melt at Concordia Subglacial Lake exceeds 5 mm a -1 . The inferred basal melt at these locations is not possible without enhanced geothermal flux. We demonstrate how radar-sounding data can provide both input and verification for a self-consistent model of vertical strain, vertical temperature distribution and meltwater distribution.

Repeated large-scale retreat and advance of Totten Glacier indicated by inland bed erosion

Climate variations cause ice sheets to retreat and advance, raising or lowering sea level by metres to decametres. The basic relationship is unambiguous, but the timing, magnitude and sources of sea-level change remain unclear; in particular, the contribution of the East Antarctic Ice Sheet (EAIS) is ill defined, restricting our appreciation of potential future change. Several lines of evidence suggest possible collapse of the Totten Glacier into interior basins during past warm periods, most notably the Pliocene epoch, causing several metres of sea-level rise. However, the structure and long-term evolution of the ice sheet in this region have been understood insufficiently to constrain past ice-sheet extents. Here we show that deep ice-sheet erosion—enough to expose basement rocks—has occurred in two regions: the head of the Totten Glacier, within 150 kilometres of today’s grounding line; and deep within the Sabrina Subglacial Basin, 350–550 kilometres from this grounding line. Our results, based on ICECAP aerogeophysical data, demarcate the marginal zones of two distinct quasi-stable EAIS configurations, corresponding to the ‘modern-scale’ ice sheet (with a marginal zone near the present ice-sheet margin) and the retreated ice sheet (with the marginal zone located far inland). The transitional region of 200–250 kilometres in width is less eroded, suggesting shorter-lived exposure to eroding conditions during repeated retreat–advance events, which are probably driven by ocean-forced instabilities. Representative ice-sheet models indicate that the global sea-level increase resulting from retreat in this sector can be up to 0.9 metres in the modern-scale configuration, and exceeds 2 metres in the retreated configuration.

Modeling 5 years of subglacial lake activity in the MacAyeal Ice Stream (Antarctica) catchment through assimilation of ICESat laser altimetry

Subglacial lakes beneath Antarcticas fast-moving ice streams are known to undergo ∼1 km 3 volume changes on annual timescales. Focusing on the MacAyeal Ice Stream (MacIS) lake system, we create a simple model for the response of subglacial water distribution to lake discharge events through assimilation of lake volume changes estimated from Ice, Cloud and land Elevation Satellite (ICESat) laser altimetry. We construct a steady-state water transport model in which known subglacial lakes are treated as either sinks or sources depending on the ICESat-derived filling or draining rates. The modeled volume change rates of five large subglacial lakes in the downstream portion of MacIS are shown to be consistent with observed filling rates if the dynamics of all upstream lakes are considered. However, the variable filling rate of the northernmost lake suggests the presence of an undetected lake of similar size upstream. Overall, we show that, for this fast-flowing ice stream, most subglacial lakes receive >90% of their water from distant distributed sources throughout the catchment, and we confirm that water is transported from regions of net basal melt to regions of net basal freezing. Our study provides a geophysically based means of validating subglacial water models in Antarctica and is a potential way to parameterize subglacial lake discharge events in large-scale ice-sheet models where adequate data are available.