Ian R. Joughin

Interferometric estimation of three-dimensional ice-flow using ascending and descending passes

Satellite radar interferometry (SRI) provides an important new tool for determining ice-flow velocity. Interferometric measurements made from a single-track direction are sensitive only to a single component of the three-component velocity vector. Observations from along three different track directions would allow the full velocity vector to be determined. A north/south-looking synthetic aperture radar (SAR) could provide these observations over large portions of the globe, but not over large areas of the polar ice sheets. The authors develop and demonstrate a technique that allows the three-component velocity vector to be estimated from data acquired along two track directions (ascending and descending) under a surface-parallel flow assumption. This technique requires that there are accurate estimates of the surface slope, which are also determined interferometrically. To demonstrate the technique, the authors estimate the three-component velocity field for the Ryder Glacier, Greenland. Their results are promising, although they do not have yet ground-truth data with which to determine the accuracy of their estimates.

Changes in west Antarctic ice stream velocities: Observation and analysis

[1]xa0We have produced a map of velocity covering much of the Siple Coast ice streams. The map confirms earlier estimates of deceleration on Whillans Ice Stream. Comparison with bed elevation data indicates that subglacial topography and the location of consolidated sediment play a strong role in determining the location of the tributaries feeding the ice streams. Force balance estimates based on these data indicate that the tributaries have beds nearly an order of magnitude stronger than those beneath many of the ice streams. We have used a theoretical analysis to examine the controls on fast flow. This analysis suggests that ice plains (very wide ice streams) are inherently unstable. This instability may be responsible for the current deceleration on the Ice Plain of Whillans Ice Stream and the shutdown of Ice Stream C 150 years ago. Thinning-induced reductions in driving stress may also explain some of the observed deceleration, particularly in upstream areas. The active portions of Ice Stream C coincide well with the areas where we estimate that melt should be taking place. Current topography and inferences of large thickening following a shutdown suggest the upstream migration of a stagnation front that initiated at the ice plain. Uncertainty remains about the basal conditions on Ice Stream D, while the basal resistance on Ice Stream E is large enough to ensure basal melting.

A Mini-Surge on the Ryder Glacier, Greenland, Observed by Satellite Radar Interferometry

Satellite radar interferometry reveals that the speed of the Ryder Glacier increased roughly threefold and then returned to normal (100 to 500 meters/year) over a 7-week period near the end of the 1995 melt season. The accelerated flow represents a substantial, though short-lived, change in ice discharge. During the period of rapid motion, meltwater-filled supraglacial lakes may have drained, which could have increased basal water pressure and caused the mini-surge. There are too few velocity measurements on other large outlet glaciers to determine whether this type of event is a widespread phenomenon in Greenland, but because most other outlet glaciers are at lower latitudes, they should experience more extensive melting, making them more susceptible to meltwater-induced surges.

Timing of Recent Accelerations of Pine Island Glacier, Antarctica

[1]xa0We have used Interferometric Synthetic Aperture Radar (InSAR) data and sequential Landsat imagery to identify and temporally constrain two acceleration events on Pine Island Glacier (PIG). These two events are separated by a period of at least seven years (1987–1994). The change in discharge between two flux gates indicates that the majority of the increase in discharge associated with the second acceleration originates well inland (>80 km) from the grounding line. An analysis indicates that changes in driving stress consistent with observed thinning rates are sufficient in magnitude to explain much of the acceleration.

Melting and freezing beneath Filchner‐Ronne Ice Shelf, Antarctica

[1]xa0We use remote-sensing data sets to evaluate the spatial distribution of melt beneath the Filchner-Ronne Ice Shelf (FRIS). The net melt rate of 83.4 ± 24.8 Gtons/yr is 2.5–5 times lower than previous glaciological estimates, but is similar to existing oceanographic estimates. The spatial distribution, however, differs significantly from standard conceptual and numerical models in which most melt occurs along the grounding lines. Our results suggest most grounding-line melt is refrozen, while the dominant Ice Shelf Water (ISW) source is melting near the ice shelf front, probably associated with tidal action. This suggests that changes in ice shelf extent can impact ISW production rates in the Weddell Sea.

Response of subglacial sediments to basal freeze‐on 2. Application in numerical modeling of the recent stoppage of Ice Stream C, West Antarctica

[1]xa0Ross ice streams supply over 90% of the ice volume flowing out of the Ross sector of the West Antarctic ice sheet (WAIS). Stoppage of Ice Stream C (ISC) ca. 150 years ago appears to have pushed this sector of WAIS from negative into positive mass balance [Joughin and Tulaczyk, 2002]. We propose an explanation for the unsteady behavior of ISC using a new numerical ice-stream model, which includes an explicit treatment of a subglacial till layer. When constrained by initial conditions emulating prestoppage geometry, dynamics, and mass balance of ISC, the model yields a rapid (∼100 years) stoppage of the main ice-stream trunk. The stoppage is triggered by basal freeze-on, which consolidates and strengthens the subglacial till. Our numerical simulations produce results consistent with a number of existing observations, for example, continuing activity of the two tributaries of ISC. The model always yields rapid stoppage unless we specify ice-stream width that is smaller than its prestoppage values (maximum of ∼80 km). We conjecture that if ISC was active for at least a few thousand years before slowdown, its width was significantly smaller than today to sustain the long active phase. Ice-stream width is a key control that helps determine whether ice-stream flow is sustainable over a long term. Our work indicates that the recent stoppage of Ice Stream C could have been part of inherent ice-stream cyclicity, and it leaves open the possibility that other active ice streams may evolve in the future toward rapid shutdowns.

Balance velocities of the Greenland ice sheet

We present a map of balance velocities for the Greenland ice sheet. The resolution of the underlying DEM, which was derived primarily from radar altimetery data, yields far greater detail than earlier balance velocity estimates for Greenland. The velocity contours reveal in striking detail the location of an ice stream in northeastern Greenland, which was only recently discovered using satellite imagery. Enhanced flow associated with all of the major outlets is clearly visible, although small errors in the source data result in less accurate estimates of the absolute flow speeds. Nevertheless, the balance map is useful for ice-sheet modelling, mass balance studies, and field planning.

Ice-stream-related patterns of ice flow in the interior of northeast Greenland

Ice flow in the interior of the NE quadrant of the Greenland ice sheet is focused on the large ice stream draining the north side of the summit dome. The rapid ice flow in the stream is apparent in the surface features in the stream and at the margins, in the broad scale topography that drives ice flow, and in satellite-derived ice motion information. The patterns of ice flow in the upper half of the stream are remarkable for their level of organization, simple geometry, and effects on local surface topography. The stream begins less than 100 km from the ice divide as a current 15 km wide and then broadens symmetrically downstream by the addition of ice from the sides to a width of more than 60 km. Elevation data and visible-band imagery show that the stream has marginal troughs tens of meters deep in its upper reach which are coincident with regions of high shear strain rate. The topography of the margins and undulating surface of the stream is generated by the ice flow; the surface undulations in the stream are fixed in location and shape over the 5 year period from 1994 to 1999. The enhanced flow presents a challenge for researchers trying to understand the history of ice discharge from a significant area in the interior of the ice sheet.

Antarctic Firn Compaction Rates from Repeat-Track Airborne Radar Data: I. Methods

Abstract While measurements of ice-sheet surface elevation change are increasingly used to assess mass change, the processes that control the elevation fluctuations not related to ice-flow dynamics (e.g. firn compaction and accumulation) remain difficult to measure. Here we use radar data from the Thwaites Glacier (West Antarctica) catchment to measure the rate of thickness change between horizons of constant age over different time intervals: 2009–10, 2010–11 and 2009–11. The average compaction rate to ~25 m depth is 0.33 m a−1, with largest compaction rates near the surface. Our measurements indicate that the accumulation rate controls much of the spatio-temporal variations in the compaction rate while the role of temperature is unclear due to a lack of measurements. Based on a semi-empirical, steady-state densification model, we find that surveying older firn horizons minimizes the potential bias resulting from the variable depth of the constant age horizon. Our results suggest that the spatio-temporal variations in the firn compaction rate are an important consideration when converting surface elevation change to ice mass change. Compaction rates varied by up to 0.12 m a−1 over distances <6 km and were on average >20% larger during the 2010–11 interval than during 2009–10.

A comparison of balance velocities, measured velocities and thermomechanically modelled velocities for the Greenland Ice Sheet

Abstract Balance velocities for the Greenland ice sheet have been calculated from a new digital-elevation model, accumulation-rates compilation and an existing ice-thickness grid, using a two-dimensional finite-difference scheme. The pattern of velocities over the ice sheet is presented and compared with velocities derived from synthetic-aperture-radar interferometry for part of northern Greenland and a limited number of global positioning system data. This comparison indicated that the balance-velocity scheme and boundary conditions used here provide a remarkably good representation of the dynamics of the ice sheet inland from the margins. It is suggested, therefore, that these balance-velocity data could provide a valuable method of constraining a numerical ice-sheet model. The balance velocities were compared with the diagnostic velocity field calculated from several different configurations of a numerical ice-sheet model. The general pattern of flow agrees well. The detail, however, is quite different. For example, the large (>300km) ice stream in the northeast is not generated by the numerical model and much of the detailed flow pattern is completely lost due to the limited model resolution and limitations in the model physics.