Grounding-line basal melt rates determined using radar-derived internal stratigraphy

We use ice-penetrating radar data across grounding lines of Siple Dome and Roosevelt Island, Antarctica, to measure the spatial pattern, magnitude and duration of sub-ice-shelf melting at these locations. Stratigraphic layers across the grounding line show, in places, a large-amplitude downwarp at, or slightly downstream of, the grounding line due to sub-ice-shelf basal melting. Localized downwarping indicates that melting is transient; melt rates, or the grounding line position, have changed within a few hundred years in order to produce the observed stratigraphy. Elsewhere, no meltrelated stratigraphic signature is preserved. In part, heterogeneity in the amount of sub-ice-shelf melt is due to regional circulation patterns in the sub-shelf cavity, but local (on the order of tens of kilometers) heterogeneity in the melt pattern may reflect small differences in the shape of the ice-shelf base at the grounding line. We find that all of the grounding lines crossed have been in place for at most ~400 years

The Link Between Climate Warming and Break-Up of Ice Shelves in the Antarctic Peninsula

A review of in situ and remote-sensing data covering the ice shelves of the Antarctic Peninsula provides a series of characteristics closely associated with rapid shelf retreat: deeply embayed ice fronts; calving of myriad small elongate bergs in punctuated events; increasing flow speed; and the presence of melt ponds on the ice-shelf surface in the vicinity of the break-ups. As climate has warmed in the Antarctic Peninsula region, melt-season duration and the extent of ponding have increased. Most break-up events have occurred during longer melt seasons, suggesting that meltwater itself, not just warming, is responsible. Regions that show melting without pond formation are relatively unchanged. Melt ponds thus appear to be a robust harbinger of ice-shelf retreat. We use these observations to guide a model of ice-shelf flow and the effects of meltwater. Crevasses present in a region of surface ponding will likely fill to the brim with water. We hypothesize (building on Weertman (1973), Hughes (1983) and Van der Veen (1998)) that crevasse propagation by meltwater is the main mechanism by which ice shelves weaken and retreat. A thermodynamic finite-element model is used to evaluate ice flow and the strain field, and simple extensions of this model are used to investigate crack propagation by meltwater. The model results support the hypothesis

Post-Stagnation Behavior in the Upstream Regions of Ice Stream C, West Antarctica

The region where two active tributaries feed into the now stagnant Ice Stream C (ISC), West Antarctica, is thickening. In this region, we observe a correlation between faster ice flow (the tributaries) and elevated topography. We conclude that stagnation of ISC resulted in compression and thickening along the tributaries, eventually forming a bulge on the ice-sheet surface. Modern hydraulic potential gradients would divert basal meltwater from ISC to Ice Stream B (ISB). These gradients are primarily controlled by the bulge topography, and so likely formed subsequent to trunk stagnation. As such, we argue against water piracy as being the cause for ISC\u27s stagnation. Kinematic-wave theory suggests that thickness perturbations propagate downstream over time, but that kinematic-wave speed decreases near the stagnant trunk. This and modest diffusion rates combine to trap most of the tributary-fed ice in the bulge region. Using interferometric synthetic aperture radar velocity measurements, we observe that half of the ice within ISC\u27s southern tributary flows into ISB. That flow pattern and other observations of non-steady flow in the region likely result from stagnation-induced thickening along upper ISC combined with a longer period of thinning on upper ISB. If current trends in thickness change continue, more ice from upper ISC will be diverted to ISB

The shallow shelf approximation as a "sliding law" in a thermomechanically coupled ice sheet model

The shallow shelf approximation is a better ``sliding law'' for ice sheet modeling than those sliding laws in which basal velocity is a function of driving stress. The shallow shelf approximation as formulated by \emph{Schoof} [2006a] is well-suited to this use. Our new thermomechanically coupled sliding scheme is based on a plasticity assumption about the strength of the saturated till underlying the ice sheet in which the till yield stress is given by a Mohr-Coulomb formula using a modeled pore water pressure. Using this scheme, our prognostic whole ice sheet model has convincing ice streams. Driving stress is balanced in part by membrane stresses, the model is computable at high spatial resolution in parallel, it is stable with respect to parameter changes, and it produces surface velocities seen in actual ice streams.Comment: 12 pages of text; 4 tables; 27 figures; submitted to JGR Earth Surfac

Consistent generation of ice-streams via thermo-viscous instabilities modulated by membrane stresses

Accurate computation of ice-stream location and dynamics is a key aspiration for theoretical glaciology. Ice-sheet models with thermo-viscous coupling have been shown to exhibit stream-like instabilities using shallow-ice approximation mechanics, but the location and width of these streams depends on the numerical implementation and are not unique. We present results from thermo-viscously coupled ice-sheet models incorporating membrane stresses. Spontaneous generation of fast-flowing linear features still occurs under certain parameter regimes, with computed stream widths between 20 km to 100 km, comparable with observations. These features are maintained as the grid-size is decreased. The thermo-viscous feedback mechanism that generates ice-streams under the shallow ice approximation still operates, now selecting a unique stream size. Computations of thermo-viscous ice flows should include membrane stresses when the bed is approximately flat, e. g. parts of Antarctica and former ice-sheets of the Northern hemisphere. Previous calculations of spontaneous ice-stream generation using the shallow ice approximation should be reassessed. Citation: Hindmarsh, R. C. A. (2009), Consistent generation of ice-streams via thermo-viscous instabilities modulated by membrane stresses, Geophys. Res. Lett., 36, L06502, doi:10.1029/2008GL036877

Multidecadal Basal Melt Rates and Structure of the Ross Ice Shelf, Antarctica, Using Airborne Ice Penetrating Radar

Basal melting of ice shelves is a major source of mass loss from the Antarctic Ice Sheet. In situ measurements of ice shelf basal melt rates are sparse, while the more extensive estimates from satellite altimetry require precise information about firn density and characteristics of near‐surface layers. We describe a novel method for estimating multidecadal basal melt rates using airborne ice penetrating radar data acquired during a 3‐year survey of the Ross Ice Shelf. These data revealed an ice column with distinct upper and lower units whose thicknesses change as ice flows from the grounding line toward the ice front. We interpret the lower unit as continental meteoric ice that has flowed across the grounding line and the upper unit as ice formed from snowfall onto the relatively flat ice shelf. We used the ice thickness difference and strain‐induced thickness change of the lower unit between the survey lines, combined with ice velocities, to derive basal melt rates averaged over one to six decades. Our results are similar to satellite laser altimetry estimates for the period 2003–2009, suggesting that the Ross Ice Shelf melt rates have been fairly stable for several decades. We identify five sites of elevated basal melt rates, in the range 0.5–2 m a⁻¹, near the ice shelf front. These hot spots indicate pathways into the sub‐ice‐shelf ocean cavity for warm seawater, likely a combination of summer‐warmed Antarctic Surface Water and modified Circumpolar Deep Water, and are potential areas of ice shelf weakening if the ocean warms

Seismic and geodetic evidence for grounding-line control of Whillans Ice Stream stick-slip events

The tidally modulated, stick‐slip events of Whillans Ice Stream in West Antarctica produce seismic energy from three locations near the grounding line. Using ice velocity records obtained by combining time series from colocated broadband seismometers and GPS receivers installed on the ice stream during the 2010–2011 and 2011–2012 austral summers, along with far‐field seismic recordings of elastic waves, we locate regions of high rupture velocity and stress drop. These regions, which are analogous to “asperities” in traditional seismic fault studies, are areas of elevated friction at the base of the ice stream. Slip events consistently initiate at one of two locations: near the center of the ice stream, where events associated with the Ross Sea high tide originate, or a grounding‐line spot, where events associated with the Ross Sea low tide initiate, as well as occasional high‐tide events following a skipped low‐tide event. The grounding‐line site, but not the central site, produces Rayleigh waves observable up to 1000 km away, through fast expansion of the slip area. Grounding‐line initiation events also show strong directivity in the downstream direction, indicating initial rupture propagation at 1.5 km/s, compared to an average of 0.150 km/s for the entire slip event. Following slip initiation, additional seismic energy is produced from two sources located near the grounding line: first at the downstream end of Subglacial Lake Engelhardt and second toward the farthest downstream extent of the ice stream. This evidence suggests that the stronger, higher‐friction material along the grounding line controls motion throughout the stick‐slip region