Dustin M. Schroeder

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.

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.

Rapid submarine ice melting in the grounding zones of ice shelves in West Antarctica

Enhanced submarine ice-shelf melting strongly controls ice loss in the Amundsen Sea embayment (ASE) of West Antarctica, but its magnitude is not well known in the critical grounding zones of the ASEs major glaciers. Here we directly quantify bottom ice losses along tens of kilometres with airborne radar sounding of the Dotson and Crosson ice shelves, which buttress the rapidly changing Smith, Pope and Kohler glaciers. Melting in the grounding zones is found to be much higher than steady-state levels, removing 300–490 m of solid ice between 2002 and 2009 beneath the retreating Smith Glacier. The vigorous, unbalanced melting supports the hypothesis that a significant increase in ocean heat influx into ASE sub-ice-shelf cavities took place in the mid-2000s. The synchronous but diverse evolutions of these glaciers illustrate how combinations of oceanography and topography modulate rapid submarine melting to hasten mass loss and glacier retreat from West Antarctica.

Surface slope control on firn density at Thwaites Glacier, West Antarctica: Results from airborne radar sounding

We derive the surface density variations over Thwaites Glacier from a statistical analysis of airborne High-Capability Radar Sounder surface returns acquired in the 2004–2005 summer. We produce a 5 km gridded map with an estimated ±12.5 kg m−3 accuracy. The background pattern of densities decreases inland from ~ 480 kg m−3 to ~ 350 kg m−3. A remarkable ~ 30 km wide vein-shaped anomaly of up to 570 kg m−3 is located ~ 100 km from the coastline. Density anomalies correspond with surface slope breaks but are not necessarily coincident with smaller slopes. They could result from complex wind-driven snow redistribution and/or refreezing of former snowmelt. This inversion technique can significantly improve surface mass balance calculations to understanding of glacier dynamics at regional scales and is valuable to verify and improve Antarctic climate models. It is also a promising approach for future surface analysis of icy moons by planetary radars.

Airborne radar sounding evidence for deformable sediments and outcropping bedrock beneath Thwaites Glacier, West Antarctica

The geologic and morphologic records of prior ice sheet configurations show evidence of rapid, back-stepping, meltwater intensive retreats. However, the potential for such a retreat in a contemporary glacier depends on the lithology of the current ice sheet bed, which lies beneath kilometers of ice, making its physical properties difficult to constrain. We use radar sounding and marine bathymetry data to compare the bed configuration of Thwaites Glacier to the bed of paleo-Pine Island Glacier. Using observed and modeled radar scattering, we show that the tributaries and upper trunk of Thwaites Glacier are underlain by ice flow-aligned bedforms consistent with deformable sediment and that the lower trunk is grounded on a region of high bed roughness consistent with outcropping bedrock. This is the same configuration as paleo-Pine Island Glacier during its retreat across the inner continental shelf.

The distribution of basal water between Antarctic subglacial lakes from radar sounding

Antarcticas subglacial lakes have two end member geophysical expressions: as hydraulically flat, radar reflective regions highlighted in ice surface topography and radar sounding profiles (‘definite lakes’), and as localized sites of elevation change identified from repeat elevation observations (‘active lakes’) that are often found in fast flowing ice streams or enhanced ice flow tributaries. While ‘definite lakes’ can be identified readily by high bed reflectivity in radar sounding, the identification and characterization of less distinct subglacial lakes and water systems with radar sounding are complicated by variable radio-wave attenuation in the overlying ice. When relying on repeat elevation observations, the relatively short times series and biased distribution of elevation observations, along with the episodic nature of ‘active lake’ outflow and replenishment, limit our understanding of how water flows under the ice sheet. Using recently developed methods for quantifying the radar scattering behaviour of the basal interface of the ice, we can avoid the problem of attenuation, and observe the plumbing of the subglacial landscape. In West Antarcticas Ross Sea Embayment, we confirm that extensive distributed water systems underlie these ice streams. Distributed water sheets are upstream in the onset regions of fast flow, while canal systems underly downstream regions of fast flow. In East Antarctica, we use specularity analysis to recover substantial hydraulic connectivity extending beyond previous knowledge, connecting the lakes already delineated by traditional radar sounding or surface elevation transients.

Estimating Subglacial Water Geometry Using Radar Bed Echo Specularity: Application to Thwaites Glacier, West Antarctica

Airborne radar sounding is an established tool for observing the bed conditions and subglacial hydrology of ice sheets and glaciers. The specularity content of radar bed echoes has also been used to detect the hydrologic transition of a subglacial water system from a network of distributed canals to a network of concentrated channels beneath the Thwaites Glacier. However, the physical dimensions of the distributed water bodies in these networks have not been constrained by observations. In this letter, we use a variety of simple radar scattering, attenuation, and cross-sectional models to provide a first estimate of the subglacial water body geometries capable of producing the observed anisotropic specularity of the Thwaites Glacier catchment. This approach leads to estimates of ice/water interface root mean square roughnesses less than about 15 cm, thicknesses of more than about 5 cm, lengths of more than about 15 m, and widths between about 0.5 and 5 m.

Extensive winter subglacial water storage beneath the Greenland Ice Sheet

W.C. is a recipient of the NASA Earth and Space Science Fellowship. D.M.S. is supported by a grant from the NASA Cryospheric Sciences Program. H.S. is supported by grants from the NASA Cryospheric Sciences and Sea Level Rise Programs. T.T.C and R.E.B are supported by grants from National Science Foundation (NSF) and NASA Cryospheric Sciences. S.P. is supported by the Natural Environment Research Council’s Centre for Polar Observation