Whillans Ice Stream Research Articles

Monitoring Activity of the Subglacial Lake Mercer Based on CryoSat-2 and ICESat-2 Altimetry Data from 2011 to 2023

Abstract. Active subglacial lakes are capable of exchanging and transporting water through water flow paths, thus having an important impact on ice sheet movement and even on the regional mass balance. The use of satellite altimetry can indirectly reveal internal subglacial hydrological activity through direct surface elevation observation, which plays an important role in the study of the evolution of active subglacial lakes. This paper monitors and analyzes the activity of the subglacial lake Mercer (SLM) in the Whillans and Mercer Ice Stream (WIS/MIS), using data from the CryoSat-2 satellite radar altimetry and ICESat-2 satellite laser altimetry. We adopt the differential DEM model and the repeat orbit model based on their respective data characteristics. Finally, we reveal the temporal and spatial evolution patterns of SLM during 2011-2023, and construct the elevation anomaly time series. The results show that the central and western parts of SLM were the areas with more significant elevation anomalies, and there were three significant fill-drain cycles during the time periods studied in this research, with the elevation anomaly range reaching about 8 m and 4 m for the first and second cycles, respectively, and the lake is currently in the third significant draining phase. Furthermore, the elevation anomalies obtained during the mission overlap time of CryoSat-2 and ICESat-2 data are in good agreement, which confirms the accuracy of our method.

Entrained Water in Basal Ice Suppresses Radar Bed‐Echo Power at Active Subglacial Lakes

AbstractSubglacial lakes have been mapped across Antarctica with two methods, radio‐echo sounding (RES) and ice‐surface deformation. At sites where both are coincident, these methods typically provide conflicting interpretations about the ice‐bed interface. With a single exception, active subglacial lakes identified by surface deformation do not display the expected flat, bright, and specular bed reflection in RES data, characteristic of non‐active lakes. This observational conundrum suggests that our understanding of Antarctic subglacial hydrology, especially beneath important fast‐moving ice streams, remains incomplete. Here, we use an airborne RES campaign that surveyed a well‐characterized group of active subglacial lakes on lower Mercer and Whillans ice streams, West Antarctica, to explore inconsistency between the two observational techniques. We test hypotheses of increased scattering and attenuation due to the presence of an active subglacial lake system that could suppress reflected bed‐echo power for RES observations in these locations, finding that entrained water is most plausible.

Towards the systematic reconnaissance of seismic signals from glaciers and ice sheets – Part 1: Event detection for cryoseismology

Abstract. Cryoseismology is a powerful toolset for progressing the understanding of the structure and dynamics of glaciers and ice sheets. It can enable the detection of hidden processes such as brittle fracture, basal sliding, transient hydrological processes, and calving. Addressing the challenge of detecting signals from many different processes, we present a novel approach for the semi-automated detection of events and event-like noise, which is well-suited for use as Part 1 of a workflow where unsupervised machine learning will be used as Part 2 (Latto et al., 2024) to facilitate the main reconnaissance of diverse detected event types. Implemented in the open-source and widely used ObsPy Python package, the multi-STA/LTA algorithm constructs a hybrid characteristic function from a set of short-term average (sta)–long-term average (lta) pairs (refer to Sect. 2 in the main text for an explanation of how uppercase and lowercase STA/sta and LTA/lta abbreviations are differentiated). We apply the algorithm to data from a seismic array deployed on the Whillans Ice Stream (WIS) in West Antarctica (austral summer 2010–2011) to form a “catch-all” catalogue of events and event-like noise. The new algorithm compares favorably with standard approaches, yielding a diversity of seismic events, including all previously identified stick-slip events (Pratt et al., 2014), teleseisms, and other noise-type signals. In terms of a catalogue overview, we investigate a partial association of seismicity with the tidal cycle and a slight association with ice temperature changes of the Antarctic summer. The new algorithm and workflow will assist in the comparison of different glacier environments using seismology, the identification of process change over time, and the targeting of possible subsequent high-resolution studies.

Towards the systematic reconnaissance of seismic signals from glaciers and ice sheets – Part 2: Unsupervised learning for source process characterization

Abstract. Given the high number and diversity of events in a typical cryoseismic dataset, in particular those recorded on ice sheet margins, it is desirable to use a semi-automated method of grouping similar events for reconnaissance and ongoing analysis. We present a workflow for employing semi-unsupervised cluster analysis to inform investigations of the processes occurring in glaciers and ice sheets. In this demonstration study, we make use of a seismic event catalogue previously compiled for the Whillans Ice Stream, for the 2010–2011 austral summer (outlined in Part 1, Latto et al., 2024). We address the challenges of seismic event analysis for a complex wave field by clustering similar seismic events into groups using characteristic temporal, spectral, and polarization attributes of seismic time series with the k-means++ algorithm. This provides the basis for a reconnaissance analysis of a seismic wave field that contains local events (from the ice stream) set in an ambient wave field that itself contains a diversity of signals (mostly from the Ross Ice Shelf). As one result, we find that two clusters include stick-slip events that diverge in terms of length and initiation locality (i.e., central sticky spot and/or the grounding line). We also identify a swarm of high-frequency signals on 16–17 January 2011 that are potentially associated with a surface melt event from the Ross Ice Shelf. Used together with the event detection presented in Part 1, the semi-automated workflow could readily be generalized to other locations and, as a possible benchmark procedure, could enable the monitoring of remote glaciers over time and comparisons between locations.

Ross Ice Shelf Displacement and Elastic Plate Waves Induced by Whillans Ice Stream Slip Events

AbstractIce shelves are assumed to flow steadily from their grounding lines to the ice front. We report the detection of ice‐propagating extensional Lamb (plate) waves accompanied by pulses of permanent ice shelf displacement observed by co‐located Global Navigation Satellite System receivers and seismographs on the Ross Ice Shelf. The extensional waves and associated ice shelf displacement are produced by tidally triggered basal slip events of the Whillans Ice Stream, which flows into the ice shelf. The propagation velocity of 2,800 m/s is intermediate between shear and compressional ice velocities, with velocity and particle motions consistent with predictions for extensional Lamb waves. During the passage of the Lamb waves the entire ice shelf is displaced about 60 mm with a velocity more than an order of magnitude above its long‐term flow rate. Observed displacements indicate a peak dynamic strain of 10−7, comparable to that of earthquake surface waves that trigger ice quakes.

Seismic attenuation in Antarctic firn

Abstract. We estimate the seismic attenuation of P and S waves in the polar firn and underlying ice by spectral analysis of diving, refracted, and reflected waves from active-source three-component seismic signals obtained in 2010 on the Whillans Ice Stream (WIS), a fast-flowing ice stream in West Antarctica. The resulting quality factors are then successfully modeled using a rock-physics theory of wave propagation that combines White's mesoscopic attenuation theory of interlayer flow with that of Biot/squirt flow. The first theory describes an equivalent viscoelastic medium consisting of a stack of two alternating thin porous layers, both of which have thicknesses that are much greater than the pore size but smaller than the wavelength. On the other hand, in the so-called Biot/squirt-flow model, there are two loss mechanisms, namely the global Biot flow and the local flow from fluid-filled microcracks (or grain contacts) to the pore space and back, where the former is dominant over the latter. The fluid saturating the pores is assumed to be fluidized snow, defined as a mixture of snow particles and air, such as powder, with a rigidity modulus of zero.

Insights into the vulnerability of Antarctic glaciers from the ISMIP6 ice sheet model ensemble and associated uncertainty

Abstract. The Antarctic Ice Sheet represents the largest source of uncertainty in future sea level rise projections, with a contribution to sea level by 2100 ranging from −5 to 43 cm of sea level equivalent under high carbon emission scenarios estimated by the recent Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). ISMIP6 highlighted the different behaviors of the East and West Antarctic ice sheets, as well as the possible role of increased surface mass balance in offsetting the dynamic ice loss in response to changing oceanic conditions in ice shelf cavities. However, the detailed contribution of individual glaciers, as well as the partitioning of uncertainty associated with this ensemble, have not yet been investigated. Here, we analyze the ISMIP6 results for high carbon emission scenarios, focusing on key glaciers around the Antarctic Ice Sheet, and we quantify their projected dynamic mass loss, defined here as mass loss through increased ice discharge into the ocean in response to changing oceanic conditions. We highlight glaciers contributing the most to sea level rise, as well as their vulnerability to changes in oceanic conditions. We then investigate the different sources of uncertainty and their relative role in projections, for the entire continent and for key individual glaciers. We show that, in addition to Thwaites and Pine Island glaciers in West Antarctica, Totten and Moscow University glaciers in East Antarctica present comparable future dynamic mass loss and high sensitivity to ice shelf basal melt. The overall uncertainty in additional dynamic mass loss in response to changing oceanic conditions, compared to a scenario with constant oceanic conditions, is dominated by the choice of ice sheet model, accounting for 52 % of the total uncertainty of the Antarctic dynamic mass loss in 2100. Its relative role for the most dynamic glaciers varies between 14 % for MacAyeal and Whillans ice streams and 56 % for Pine Island Glacier at the end of the century. The uncertainty associated with the choice of climate model increases over time and reaches 13 % of the uncertainty by 2100 for the Antarctic Ice Sheet but varies between 4 % for Thwaites Glacier and 53 % for Whillans Ice Stream. The uncertainty associated with the ice–climate interaction, which captures different treatments of oceanic forcings such as the choice of melt parameterization, its calibration, and simulated ice shelf geometries, accounts for 22 % of the uncertainty at the ice sheet scale but reaches 36 % and 39 % for Institute Ice Stream and Thwaites Glacier, respectively, by 2100. Overall, this study helps inform future research by highlighting the sectors of the ice sheet most vulnerable to oceanic warming over the 21st century and by quantifying the main sources of uncertainty.

Trace element, rare earth element and trace carbon compounds in Subglacial Lake Whillans, West Antarctica

Whillans Subglacial Lake (SLW) lies beneath 801 m of ice in the lower portion of the Whillans Ice Stream (WIS) in West Antarctica and is part of an extensive and active subglacial drainage network. Here, the geochemical characterization of SLW rare earth elements (REE), trace elements (TE), free amino acids (FAA), and phenolic compounds (PC) measured in lakewater and sediment porewater are reported. The results show, on average, higher values of REEs in the lakewater than in the porewater, and clear changes in all REE concentrations and select redox sensitive trace element concentrations in porewaters at a depth of ~15 cm in the 38 cm lake sediment core. This is consistent with prior results on the lake sediment redox conditions based on gas chemistry and microbiological data. Low concentrations of vanillyl phenols were measured in the SLW water column with higher concentrations in porewater samples and their concentration profiles in the sediments may also reflect changing redox conditions in the sediments. Vanillin concentrations increased with depth in the sediments as oxygenation decreases, while the concentrations of vanillic acid, the more oxidized component, were higher in the more oxygenated surface sediments. Collectively these results indicate redox changes occurring with the upper 38 cm of sediment in SLW and provide support for the existence of a seawater source, already hypothesized, in the sediments below the lowest measured depth, and of a complex and dynamic geochemical system beneath the West Antarctic Ice Sheet.Our results are the first to detail geochemical properties from an Antarctic subglacial environment using direct sampling technology. Due to their isolation from the wider environment, subglacial lakes represent one of our planets last pristine environments that provide habitats for microbial life and natural biogeochemical cycles but also impact the basal hydrology and can cause ice flow variations.

Identification and Activity of Subglacial Lakes Beneath the Mercer and Whillans Ice Streams and Slessor Glacier

Identification and Activity of Subglacial Lakes Beneath the Mercer and Whillans Ice Streams and Slessor Glacier

Radar attenuation demonstrates advective cooling in the Siple Coast ice streams

AbstractIce streams are warmed by shear strain, both vertical shear near the bed and lateral shear at the margins. Warm ice deforms more easily, establishing a positive feedback loop in an ice stream where fast flow leads to warm ice and then to even faster flow. Here, we use radar attenuation measurements to show that the Siple Coast ice streams are colder than previously thought, which we hypothesize is due to along-flow advection of cold ice from upstream. We interpret the attenuation results within the context of previous ice-temperature measurements from nearby sites where hot-water boreholes were drilled. These in-situ temperatures are notably colder than model predictions, both in the ice streams and in an ice-stream shear margin. We then model ice temperature using a 1.5-dimensional numerical model which includes a parameterization for along-flow advection. Compared to analytical solutions, we find depth-averaged temperatures that are colder by 0.7°C in the Bindschadler Ice Stream, 2.7°C in the Kamb Ice Stream and 6.2–8.2°C in the Dragon Shear Margin of Whillans Ice Stream, closer to the borehole measurements at all locations. Modelled cooling corresponds to shear-margin thermal strengthening by 3–3.5 times compared to the warm-ice case, which must be compensated by some other weakening mechanism such as material damage or ice-crystal fabric anisotropy.

A dynamic saline groundwater system mapped beneath an Antarctic ice stream

Antarctica's fast-flowing ice streams drain the ice sheet, with their velocity modulated by subglacial water systems. Current knowledge of these water systems is limited to the shallow portions near the ice-bed interface, but hypothesized deeper groundwater could also influence ice streaming. Here, we use magnetotelluric and passive seismic data from Whillans Ice Stream, West Antarctica, to provide the first observations of deep sub-ice stream groundwater. Our data reveal a volume of groundwater within a >1-kilometer-thick sedimentary basin that is more than an order of magnitude larger than the known subglacial system. A vertical salinity gradient indicates exchange between paleo seawater at depth and contemporary basal meltwater above. Our results provide new constraints for subglacial water systems that affect ice streaming and subglacial biogeochemical processes.

Post-Processing Synchronized Bistatic Radar for Long Offset Glacier Sounding

Radar tomography of glaciers promises to improve imaging and estimates of subsurface ice-sheet structures and properties, including temperature distributions, basal materials, ice fabric, and englacial water content. However, bistatic radar data with long (i.e., larger than the ice thickness) walk-away surveys are required to constrain high-fidelity tomographic inversions. These long-offset data have proven difficult to collect due to the hardware complexity of existing synchronization techniques. Therefore, we remove the hardware complexity required for real-time synchronization by synchronizing in postprocessing. Our technique transforms an Autonomous phase-sensitive Radio Echo Sounder (ApRES) system and a software-defined radio receiver into a coherent bistatic radar capable of recovering basal echoes at long offsets. We validated our system at Whillans Ice Stream, West Antarctica, with a walk-away survey up to 1300 m (797 m thick) and at Store Glacier, Greenland, up to 1450 m (1028 m thick). At both field sites, we measured the basal echo at angles beyond the point of total internal reflection (TIR), whose previous literature had set as a hard physical limit. We support our experimental results with high-frequency structure simulation, which shows that ground-based radar systems capture evanescent waves and are not hindered by TIR. Our analysis and experiments demonstrate a system capable of executing wide-angle bistatic radar surveys for improved geometric and radiometric resolution of inversions for englacial and subglacial properties.

Did Holocene climate changes drive West Antarctic grounding line retreat and readvance?

Abstract. Knowledge of past ice sheet configurations is useful for informing projections of future ice sheet dynamics and for calibrating ice sheet models. The topology of grounding line retreat in the Ross Sea sector of Antarctica has been much debated, but it has generally been assumed that the modern ice sheet is as small as it has been for more than 100 000 years (Conway et al., 1999; Lee et al., 2017; Lowry et al., 2019; McKay et al., 2016; Scherer et al., 1998). Recent findings suggest that the West Antarctic Ice Sheet (WAIS) grounding line retreated beyond its current location earlier in the Holocene and subsequently readvanced to reach its modern position (Bradley et al., 2015; Kingslake et al., 2018). Here, we further constrain the post-LGM (Last Glacial Maximum) grounding line retreat and readvance in the Ross Sea sector using a two-phase model of radiocarbon input and decay in subglacial sediments from six sub-ice sampling locations. In addition, we reinterpret high basal temperature gradients, measured previously at three sites in this region (Engelhardt, 2004), which we explain as resulting from recent ice shelf re-grounding accompanying grounding line readvance. At one location – Whillans Subglacial Lake (SLW) – for which a sediment porewater chemistry profile is known, we estimate the grounding line readvance by simulating ionic diffusion. Collectively, our analyses indicate that the grounding line retreated over SLW 4300-2500+1500 years ago, and over sites on Whillans Ice Stream (WIS), Kamb Ice Stream (KIS), and Bindschadler Ice Stream (BIS) 4700-2300+1500, 1800-700+2700, and 1700-600+2800 years ago, respectively. The grounding line only recently readvanced back over those sites 1100-100+200, 1500-200+500, 1000-300+200, and 800±100 years ago for SLW, WIS, KIS, and BIS, respectively. The timing of grounding line retreat coincided with a warm period in the mid-Holocene to late Holocene. Conversely, grounding line readvance is coincident with cooling climate in the last 1000–2000 years. Our estimates for the timing of grounding line retreat and readvance are also consistent with relatively low carbon-to-nitrogen ratios measured in our subglacial sediment samples (suggesting a marine source of organic matter) and with the lack of grounding zone wedges in front of modern grounding lines. Based on these results, we propose that the Siple Coast grounding line motions in the mid-Holocene to late Holocene were primarily driven by relatively modest changes in regional climate, rather than by ice sheet dynamics and glacioisostatic rebound, as hypothesized previously (Kingslake et al., 2018).

Frictional Origin of Slip Events of the Whillans Ice Stream, Antarctica

AbstractIce sheet evolution depends on subglacial conditions, with the ice‐bed interface's strength exerting an outsized role on the ice dynamics. Along fast‐flowing glaciers, this strength is often controlled by the deformation of subglacial till, making quantification of spatial variations of till strength essential for understanding ice‐sheet contribution to sea‐level. This task remains challenging due to a lack of in situ observations. We analyze continuous seismic data from the Whillans Ice Plain (WIP), West Antarctica, to uncover spatio‐temporal patterns in subglacial conditions. We exploit tidally modulated stick‐slip events as a natural source of sliding variability. We observe a significant reduction of the till seismic wave‐speed between the WIP sticky‐spots. These observations are consistent with a poroelastic model where the bed experiences relative porosity and effective pressure increases of >11% during stick‐slips. We conclude that dilatant strengthening appears to be an essential mechanism in stabilizing the rapid motion of fast‐flowing ice streams.

Grounding zone subglacial properties from calibrated active-source seismic methods

Abstract. The grounding zone of Whillans Ice Stream, West Antarctica, exhibits an abrupt transition in basal properties from the grounded ice to the ocean cavity over distances of less than 0.5–1 km. Active-source seismic methods reveal the downglacier-most grounded portion of the ice stream is underlain by a relatively stiff substrate (relatively high shear wave velocities of 1100±430 m s−1) compared to the deformable till found elsewhere beneath the ice stream. Changes in basal reflectivity in our study area cannot be explained by the stage of the tide. Several kilometres upstream of the grounding zone, layers of subglacial water are detected, as are regions that appear to be water layers but are less than the thickness resolvable by our technique. The presence of stiff subglacial sediment and thin water layers upstream of the grounding zone supports previous studies that have proposed the dewatering of sediment within the grounding zone and the trapping of subglacial water upstream of the ocean cavity. The setting enables calibration of our methodology using returns from the floating ice shelf. This allows a comparison of different techniques used to estimate the sizes of the seismic sources, a constraint essential for the accurate recovery of subglacial properties. We find a strong correlation (coefficient of determination=0.46) between our calibrated method and a commonly used multiple-bounce method, but our results also highlight the incomplete knowledge of other factors affecting the amplitude of seismic sources and reflections in the cryosphere.

Mid‐Holocene Grounding Line Retreat and Readvance at Whillans Ice Stream, West Antarctica

AbstractUnderstanding ice sheet evolution through the geologic past can help constrain ice sheet models that predict future ice dynamics. Existing geological records of grounding line retreat in the Ross Sea, Antarctica, have been confined to ice‐free and terrestrial archives, which reflect dynamics from periods of more extensive ice cover. Therefore, our perspective of grounding line retreat since the Last Glacial Maximum remains incomplete. Sediments beneath Ross Ice Shelf and grounded ice offer complementary insight into the southernmost extent of grounding line retreat, yielding a more complete view of ice dynamics during deglaciation. Here we thermochemically separate the youngest organic carbon to estimate ages from sediments extracted near the Whillans Ice Stream grounding line to provide direct evidence for grounding line retreat in that region as recent as the mid‐Holocene (7.2 kyr B.P.). Our study demonstrates the utility of accurately dated, grounding‐line‐proximal sediment deposits for reconstructing past interactions between marine and subglacial environments.

Environmentally clean access to Antarctic subglacial aquatic environments

AbstractSubglacial Antarctic aquatic environments are important targets for scientific exploration due to the unique ecosystems they support and their sediments containing palaeoenvironmental records. Directly accessing these environments while preventing forward contamination and demonstrating that it has not been introduced is logistically challenging. The Whillans Ice Stream Subglacial Access Research Drilling (WISSARD) project designed, tested and implemented a microbiologically and chemically clean method of hot-water drilling that was subsequently used to access subglacial aquatic environments. We report microbiological and biogeochemical data collected from the drilling system and underlying water columns during sub-ice explorations beneath the McMurdo and Ross ice shelves and Whillans Ice Stream. Our method reduced microbial concentrations in the drill water to values three orders of magnitude lower than those observed in Whillans Subglacial Lake. Furthermore, the water chemistry and composition of microorganisms in the drill water were distinct from those in the subglacial water cavities. The submicron filtration and ultraviolet irradiation of the water provided drilling conditions that satisfied environmental recommendations made for such activities by national and international committees. Our approach to minimizing forward chemical and microbiological contamination serves as a prototype for future efforts to access subglacial aquatic environments beneath glaciers and ice sheets.

Biogeochemical Connectivity Between Freshwater Ecosystems beneath the West Antarctic Ice Sheet and the Sub‐Ice Marine Environment

AbstractAlthough subglacial aquatic environments are widespread beneath the Antarctic ice sheet, subglacial biogeochemistry is not well understood, and the contribution of subglacial water to coastal ocean carbon and nutrient cycling remains poorly constrained. The Whillans Subglacial Lake (SLW) ecosystem is upstream from West Antarctica's Gould‐Siple Coast ~800 m beneath the surface of the Whillans Ice Stream. SLW hosts an active microbial ecosystem and is part of an active hydrological system that drains into the marine cavity beneath the adjacent Ross Ice Shelf. Here we examine sources and sinks for organic matter in the lake and estimate the freshwater carbon and nutrient delivery from discharges into the coastal embayment. Fluorescence‐based characterization of dissolved organic matter revealed microbially driven differences between sediment pore waters and lake water, with an increasing contribution from relict humic‐like dissolved organic matter with sediment depth. Mass balance calculations indicated that the pool of dissolved organic carbon in the SLW water column could be produced in 4.8 to 11.9 yr, which is a time frame similar to that of the lakes’ fill‐drain cycle. Based on these estimates, subglacial lake water discharged at the Siple Coast could supply an average of 5,400% more than the heterotrophic carbon demand within Siple Coast embayments (6.5% for the entire Ross Ice Shelf cavity). Our results suggest that subglacial discharge represents a heretofore unappreciated source of microbially processed dissolved organic carbon and other nutrients to the Southern Ocean.

Estimation of ice fabric within Whillans Ice Stream using polarimetric phase-sensitive radar sounding

AbstractHere we use polarimetric measurements from an Autonomous phase-sensitive Radio-Echo Sounder (ApRES) to investigate ice fabric within Whillans Ice Stream, West Antarctica. The survey traverse is bounded at one end by the suture zone with the Mercer Ice Stream and at the other end by a basal ‘sticky spot’. Our data analysis employs a phase-based polarimetric coherence method to estimate horizontal ice fabric properties: the fabric orientation and the magnitude of the horizontal fabric asymmetry. We infer an azimuthal rotation in the prevailing horizontalc-axis between the near-surface (z≈ 10–50 m) and deeper ice (z≈ 170–360 m), with the near-surface orientated closer to perpendicular to flow and deeper ice closer to parallel. In the near-surface, the fabric asymmetry increases toward the center of Whillans Ice Stream which is consistent with the surface compression direction. By contrast, the fabric orientation in deeper ice is not aligned with the surface compression direction but is consistent with englacial ice reacting to longitudinal compression associated with basal resistance from the nearby sticky spot.

Interpreting englacial layer deformation in the presence of complex ice flow history with synthetic radargrams

AbstractFast ice flow on the Antarctic continent constitutes much of the mass loss from the ice sheet. However, geophysical methods struggle to constrain ice flow history at depth, or separate the signatures of topography, ice dynamics and basal conditions on layer structure. We develop and demonstrate a methodology to compare layer signatures in multiple airborne radar transects in order to characterize ice flow at depth, or improve coverage of existing radar surveys. We apply this technique to generate synthetic, along-flow radargrams and compare different deformation regimes to observed radargram structure. Specifically, we investigate flow around the central sticky spot of Whillans Ice Stream, West Antarctica. Our study suggests that present-day velocity flowlines are insufficient to characterize flow at depth as expressed in layer geometry, and streaklines provide a better characterization of flow around a basal sticky spot. For Whillans Ice Stream, this suggests that ice flow wraps around the central sticky spot, supported by idealized flow simulations. While tracking isochrone translation and rotation across survey lines is complex, we demonstrate that our approach to combine radargram interpretation and modeling can reveal critical details of past ice flow.