1. Introduction
Understanding the dynamical processes that govern the location of iceberg-calving fronts of ice shelves in the Ross and Weddell Seas is a long-standing problem in Antarctic ice-sheet stability and response to climate change. These processes influence directly the areal extent of the ice sheet and the net heat and salt exchanges between air and sea. It is clear that the mechanical connection between an ice shelf and the bay walls that confine it and the ice rises and islands that penetrate into it is important in maintaining a stable shelf-front position. The nature ofthat connection is not well understood, because extensive crcvassing and iceberg calving make field observations of ice-shelf margins both difficult and dangerous.
Historical data indicate that calving from large ice shelves is episodic, with long periods of ice-front advance punctuated by short periods of relatively rapid retreat (see, e.g., Reference Jacobs, MacAyeal and ArdaiJacobs and others, 1986; Rott and others, 1996). This cyclic behavior suggests that the mean front position depends on glaciological factors such as the rate of ice-spreading required to maintain contact with confining bay walls (Reference SandersonSanderson, 1979) and the formation of large-scale rifts that provide the planes of weakness along which tabular icebergs calve. However, correlations between retreat of small ice shelves and southward migration of the -5°C mean annual atmospheric surface isotherm around the Antarctic Peninsula (Reference Vaughan and DoakeVaughan and Doake, 1996) suggest that environmental factors may determine the mean ice-front position. A successful description of ice-shelf front stability must accommodate those two seemingly independent controls on ice-front location, one an “internal” glaciological control and the other an “externa]” climate control. Here, a link between internal and external controls is identified, using interfero-metric observations and finite-element modelling of the separation of the Filchner-Ronne Ice Shelf (FRIS) from Hemmen Ice Rise (HIR) and from the Lassiter Coast (LC).
The HIR and LC corners of the FRIS calving front (Fig. 1) illustrate two typical shelf-front geometries: at HIR the ice-shelf front appears to be anchored by an intruding ice rise (and by Berkner Island), and along the LC the shelf front diverges from the wall of a widening bay. Synthetic aperture radar (SAR) images (figs 2 and 3; from Rignot and MacAyeal, in press) reveal that at both locations the rifts and coastal separation cavities are filled with what appears to be a melange of multi-year sea ice, ice-shelf fragments and snow. Interferograms constructed by Rignot and MacAyeal (in press) from these SAR images show that the melange deforms, presumably in response to ice-shelf flow. Finite-element models that are tuned to match the interferograms suggest that at HIR the melange must possess sufficient mechanical strength to influence ice-shelf flow and maintain some mechanical interaction between the ice shelf and coast after the ice shelf has separated from the coast. in contrast, the melange along the LC does not appear to influence icc-shclf flow.
The link between glaciological and climate controls on the position of an ice-shelf front, suggested by the present work, is as follows. Large-scale rifting and coastal separation occur at a number of locations within an ice shelf (e.g. around ice rises and along coastal promontories). The seaward front of the ice shelf forms where weaknesses introduced by those processes result in iceberg calving. Rifting and separation in cold environmental conditions fail to produce an ice-shelf from, because the ice melange that fills the voids reconnects the ice shelf mechanically across rifts and with the coast. Where local environmental conditions do not permit the melange to bind the ice shelf together, rifting and separation produce an ice-shelf front. Moreover, we suggest that because the thin melange will melt more rapidly than the thicker icc-shclf ice, it may facilitate rapid shelf-front retreat in future-warming scenarios.
2. Data and Model Description
Icc-shclf flow and tidal motion in the LC area are observed using SAR interferograms derived from eight passes of the European remote-sensing satellites (ERS-1/2) during their tandem mission in early 1996. A sequence of seven SAR images from the 1992 ERS-1 “ice phase” mission is used to observe the area around HIR. Both sets of observations are described in Rignot and MacAyeal (in press). One ascending- and one descending-satellite-pass interferogram, filtered for tidal motion (using the interferogram differencing methods described by Rignot and MacAyeal), are used for each location.
Interpreting the pattern of ice motion depicted in the interferograms is slightly complicated because the motion is in the line-of-sight of the orbiting satellite, rather than in the plane of the ice shelf or along an ice flowline. Also, the dynamical implications of ice displacement may not be obvious. We address these limitations by visual comparison between the observations and artificial interferograms that correspond to the satellite look direction of the observations. The artificial interferograms are created using a finite-element numerical model of ice-shelf flow in the same area as that covered by the interferograms. The model (described in MacAyeal and others, in press; see also MacAyeal and others, 1996) computes ice velocity for given ice thicknesses and boundary conditions, using the typical ice-shelf stress-balance equations and assumptions about ice rheology. The finite-element mesh resolution is fine enough to trace rift boundaries and accommodate narrow boundary layers.
3. Model-Sensitivity Experiments
Two series of sensitivity experiments, in which ice softness and melange ice thickness are varied, are used to investigate the importance of coastal boundary-layer softening and of the ice melange that fills shelf-front rifts at HIR and along the LC. The cause of ice softening is not specified, but it is reasonable to assume that some softening due to strain heating and ice-crystal alignment will occur with in the rapidly shearing margin between fast-flowing ice-shelf ice and the adjacent coast (see, e.g., Reference Echelmeyer, Harrison, Larsen and MitchellEchelmeyer and others, 1994; Reference Fujita and MaeFujita and Mae, 1994). The melange is idealized simply as thin ice, with the same rhéologie properties and dynamical behavior as the surrounding ice shelf. Melange thickness will vary, due to ice-flow divergence, the age of the melange material, and environmental conditions, but because we do not intend to reproduce the interferograms exactly (as might be possible with a formal inverse method), that simplification is acceptable. Our goal is to determine the combination of model parameters which best accounts for the observed dynamical processes.
A qualitative assessment of model performance for the assumptions in each experiment is made by comparing artificial and observed interferograms. We compare fringe spacing and orientation but not absolute velocities, because the processing required to derive velocity from the interferograms cannot be performed everywhere in the images and because velocity errors are large (at best, about 15%; discussed by Rignot and MacAyeal, in press). Error in the comparisons made here should be limited to error in the observed interferometric fringes, typically about 1/28 of a fringe, and to error due to unfiltered tidal displacement.
Four experiments are conducted with the LC model. Experiment 1 assumes no boundary-layer ice softening and has open-water rifts. Experiment 2 assumes a boundary-layer ice softening of 50% (a flow-law-rate constant of 1 x 108 Pa s −1/3, reduced from 2 x 108 Pas−1/3) in a 5 km wide band along the unnamed promontory east of Hansen Inlet, extending downstream in the crevassed wake of the promontory, and has open-water rifts. Experiment 3 assumes no boundary-layer ice softening and fills the rifts with a 25 m thick ice melange. Experiment 4 assumes a boundary-layer ice softening of 50% and has a melange thickness of 25 m.
The individual HIR experiments are similar to the LC experiments and so are not described here; their details can be found in MacAyeal and others (in press).
Results of the LC experiments are depicted in Figure 4. The best match between artificial and observed interferograms is achieved in experiment 2. Boundary-layer softening is required to reproduce the narrow fringe spacing along and downstream of the promontory east of Hansen Inlet. Fringe spacing and orientation near the front of the ice shelf are best reproduced when the rifts are simulated as being ice-freC. Thus, the rift-filling ice melange observed in the LC images and interferograms does not influence flow in the main ice shelf and must not be a mechanically competent unit.
The most favorable comparisons for the HIR experiments (Fig. 5; from MacAyeal and others, in press, figs 6 and 7) are obtained with 50% coastal boundary-layer softening and a melange thickness of 10 m. Boundary-layer ice softening is required to reproduce the narrow fringe spacing around HIR and along the coast of Berkner Island, and the thin, contiguous layer of rift-filling ice must be present to reproduce the observed fringe spacing and orientation with in the main ice shelf. The melange ice must therefore be mechanically competent enough to transmit stress across rifts and between ice-shelf fragments. Its effect is to bind together the rifled shelf front.
4. Discussion
The comparison between model-derived and observed interferograms presented here demonstrates the dynamical importance of two features of ice-shelf front corners. First, coastal boundary-layer ice must be substantially softer than adjacent icc-shclf ice. At both HIR and the LC, the soft ice is advected downstream into areas where rifting begins ( in the wake of HIR in Figure 2 and downstream of the promontory east of Hansen Inlet in Figure 3). Weakness introduced into the ice shelf by this softening may determine where rifting begins and shelf fronts form. That weakness may be opposed by the second feature studied here, the ice melange that fills shelf-front rifts and coastal separation cavities. At HIR, the melange provides mechanical coupling across rift walls atid between the ice shelf and the adjacent coasts of HIR and Berkner Island. That coupling binds together the rift walls, along which large icebergs eventually calve, and it maintains a connection to Berkner Island farther seaward than horizontal spreading of ice-shell ice would allow. The ice melange probably allows the shelf front to extend beyond the stable position that would be determined by internal controls alone, by preventing or delaying iceberg calving.
The ice melange may also be the ice-shelf front's Achilles' heel, because it would be more vulnerable to melting than thicker ice-shelf ice. Indeed, melange ice, though present along the LC, oilers no mechanical coupling between the shelf front and the adjacent coast. The reason for that difference in melange competence at HIR and the LC may reflect a difference in melange thickness or temperature. Τ he LC ice melange could be warmed and melted as high-salinity shelf water Hows below the ice shelf, along the western boundary of the Weddell Sea (.Jacobs and others, 1985). Another tantalizing possibility is that the melange provides the link between shelf-front position, ice dynamics and climate. The southward progression of atmospheric warming observed around the Antarctic Peninsula (Vaughan and Doake, 1996) may be weakening the thin layer of ice. The LC is 3° north of HIR, and so the ice melange there may be responding to regional warming sooner than at HIR.
An example of the vulnerability of rift-filling ice to regional warming may be the sudden disaggregation of the Larsen-A Ice Shelf in January 1995 (Rott and others, 1996). The break-up, during a storm, followed several consecutive warm summer seasons. It is possible that water-filled crevasses penetrated the entire thickness of the shelf, creating a network of wounds held together by remaining bridges between crevasses and by frozen scawater with in the crevasses, similar to the melange observed here with in the LC and HIR rifts. The crevasse-filling melange, weakened by a particularly warm summer (suggested by extensive melt ponds on the surface of the ice shelf), would have been more vulnerable to failure during a large storm event than the thicker ice shelf. The atmospheric and sea-surface warming currently under way around the Antarctic Peninsula may progress southward, toward the Antarctic interior. That warming could weaken rift-filling sca-icc melange, making events like those on Larsen-A more common around Antarctica.
Acknowledgements
Support for this research was provided by the European Space Agency (ERS-1/2 project code: AO2.USA.160), the U.S. National Aeronautics and Space Administration (NAGW-5005) and a NASA Graduate Student Fellowship (awarded to C. L. Hulbe). We thank A. Gabriel for help in the project inception. D. Vaughan, C. Doake, R. Frölich and other colleagues at the British Antarctic Survey provided ice-thickness data for the FRIS and many suggestions during the course of the research reported here. The manuscript was improved by comments from J. A. van Zanten, G. H. Gudmundsson and an anonymous reviewer.