Introduction
Recent studies show that contemporary ice sheets, as well as ice sheets throughout the Quaternary, are susceptible to changes on decadal to millennial timescales (Reference Joughin, Tulaczyk, Bindschadler and PriceJoughin and others, 2002; Reference HemmingHemming, 2004; Reference Stearns, Jezek and van der VeenStearns and others, 2005). Changes in the geometry of ice sheets have a direct effect on global sea level as more mass is transferred from the frozen reservoirs to the world’s oceans. Sea level is currently rising and grounded ice is likely to be a major contributor to sea- level rise in the future (Reference Church and HoughtonChurch and others, 2001; Reference Huybrechts, Gregory, Janssens and WildHuybrechts and others, 2004). Most of the ice discharge from grounded ice sheets is by means of ice streams and outlet glaciers. In West Antarctica, drainage is achieved primarily through meandering fast-moving ice streams embedded in the main body of more sluggishly moving ice. These ice streams, in particular those on the Siple Coast draining into the Ross Ice Shelf, have been the subject of extensive studies and field campaigns since the early 1980s (e.g. Reference Alley and BindschadlerAlley and Bindschadler, 2001). From these studies, it has become evident that ice streams are one of the most critical components of the West Antarctic ice sheet system, yet they are also the least understood.
The causes for ice streaming, and the mechanical constraints on its dimensions and high speeds, must be understood if observed imbalances and future behavior of the West Antarctic ice sheet are to be comprehended. Ice streams are characterized by high speeds achieved under exceptionally low gravitational driving stress. On Whillans Ice Stream (formerly known as Ice Stream B), the driving stress ranges between 10 and 20 kPa, while speeds up to 800 ma–1 have been measured (Reference Whillans and van der VeenWhillans and others, 2001). This observation indicates that resistance to flow must be exceptionally small, most likely due to a weakening of subglacial sediments (Reference Alley, Blankenship, Bentley and RooneyAlley and others, 1986; Reference Tulaczyk, Kamb and EngelhardtTulaczyk and others, 2001). Reference Whillans and van der VeenWhillans and others (2001) presented an estimate of the budget of forces along the length of Whillans Ice Stream and concluded that most or all of the driving stress is balanced by drag originating at the lateral margins where faster-moving ice shears past the almost stagnant interstream ridges. Subsequent studies (Reference Joughin, Tulaczyk, Bindschadler and PriceJoughin and others, 2002; Reference Stearns, Jezek and van der VeenStearns and others, 2005) have confirmed this result.
Perhaps the most striking characteristic of the active Siple Coast ice streams is the 100-fold contrast in speed between the fast-moving ice streams and adjacent near-stagnant interstream ridges. The lateral shear margins represent the narrow along-flow zone over which the speed contrast is developed. Strong shearing across vertical planes paralleling the flow results in bands of severely crevassed ice that can be traced from deep in the interior (Reference Hodge and DoppelhammerHodge and Doppelhammer, 1996; Reference JezekJezek, 1999) to beyond the grounding line on the Ross Ice Shelf (Reference Fahnestock, Bamber, Alley and BindschadlerFahnestock and Bamber, 2001). One unresolved question is what determines the location of shear margins, and how these margins migrate laterally. In earlier models of ice-stream flow, it was assumed that the motion of ice streams was controlled from the bed, by the extent of contact with the underlying rock, or by the strength of possible mobile sediments under the ice. However, several studies have indicated that the active ice streams are largely controlled from their sides (Reference Bindschadler, Stephenson, MacAyeal and ShabtaieBindschadler and others, 1987; Reference Echelmeyer, Harrison, Larsen and MitchellEchelmeyer and others, 1994; Reference MacAyeal, Bindschadler and ScambosMacAyeal and others, 1995; Whillans and Reference Van der VeenVan der Veen, 1997; Reference Whillans and van der VeenWhillans and others, 2001; Reference Joughin, Tulaczyk, Bindschadler and PriceJoughin and others, 2002). The observation that shear margins are rather narrow indicates that the transition between well-lubricated flow on the ice stream to a flow regime dominated by basal friction on the interstream ridge must be rather sharp. On the ice stream, vertical shear is negligible and the dominant strain rate is generated by lateral shear. On the ridges, lateral shear is insignificant and flow is dominated by vertical shear. The shear margins represent the region where the shear stress is transformed from acting on vertical planes parallel to the sides to shear stress acting on planes parallel to the bed (Reference Whillans and van der VeenWhillans and Van der Veen, 1997, Reference Whillans and van der Veen2001). The nature of this transfer is not clear, however, due to the complexity of this stress transfer.
In the studies referred to above, the role of lateral drag is estimated from observations made on the ice stream itself. A different approach was adopted by Reference Whillans and van der VeenWhillans and Van der Veen (2001) who considered the stagnant ridges adjacent to the ice stream where measured speeds are typically of the order of a few meters per year (Whillans and Van der Veen, 1993). They supposed that if lateral drag supports all or part of the driving stress on the ice stream, then there must be a region of large basal drag elsewhere. It then follows that lateral drag supporting the ice stream must be transferred to basal drag under the interstream ridges, such that, integrated over the full width of the ice stream and adjacent ridges, the integrated basal drag equals the total driving force. Another possibility would be that the reduced basal drag under one part of the ice stream is accompanied by increased basal drag elsewhere under the ice stream. However, this requires along-flow transfer of flow resistance through gradients in longitudinal stress. Whillans and Van der Veen (1993) and Reference Whillans and van der VeenWhillans and others (2001) used measurements of surface velocity along the length of Whillans Ice Stream to show that longitudinal stress gradients are unimportant to the large- scale dynamics of this ice stream, indicating no along-flow transfer of resistance over larger distances. On a local scale, such as along the strain grid near the UpB Camp (located near the center of the ice stream between the ‘snake’ and ‘Dragon’), gradients in longitudinal stress can become important (Reference Hulbe and WhillansHulbe and Whillans, 1997; see also Reference Price, B¡ndschadler, Hulbe and BlankenshipPrice and others, 2002).
The objective of the present contribution is to further investigate the transfer of lateral drag to basal drag under the interstream ridge using previously unpublished measurements of surface speed derived from strain grids in the shear margins of three West Antarctic ice streams. Here, we focus on results for Whillans Ice Stream; results for the strain grids on Bindschadler and MacAyeal Ice Streams will be presented in a future paper.
Measurements
Four lines of poles were installed at the northern margin of Whillans Ice Stream, extending ~-6km from the interstream ridge into the shear margin (the ‘snake’). Poles were placed ~250 m apart, approximately perpendicular to the direction of ice-stream flow. Due to intense crevassing at shear margins, all surveys were conducted on foot by a four- or five-person rope team. The global positioning system (GPS) equipment and a Trimble 4000SSI receiver were pulled on a sled with the compact L1/L2 antenna mounted on a pole. The transects were initially surveyed in November 1998 and again in January 2000. GPS surveying was carried out using the ‘stop-and-go’ kinematic method, following Hulbe and Whillans (Reference Hulbe and Whillans1993, Reference Hulbe and Whillans1994). A continuously recording GPS receiver was located at the base camp (SNKE), and operated for the duration of each survey (typically a week or so). A fixed control in each transect, usually the pole closest to the base station, acted as a tie point between the line of poles and the base station. Kinematic GPS data (between poles) and the static GPS data at each subsequent pole were linked to this fixed point. The base station collected data at 30 s intervals. All other static (at each pole and the fixed control station) and kinematic surveys collected data at 5 s intervals.
The base station and fixed control stations were post- processed using the GIPSY-OASIS II software package, which improves the precision of the point position (Reference Zumberge, Heflin, Jefferson, Watkins and WebbZumberge and others, 1997). The program incorporates the Earth’s orientation and precise satellite orbit, as well as clock solutions derived from NASA Jet Propulsion Laboratory’s independent analysis of global fiducial stations. Formal 1σ (rms) uncertainties for the daily solutions are typically better than 0.01 m for the x, y and z positions (Reference Spikes, Csatho and WhillansSpikes and others, 2003). The daily solutions are referenced to the International Terrestrial Reference Frame 1996 (ITRF96) using the World Geodetic System 1984 (WGS84) ellipsoidal model. The stop-and-go kinematic survey data for each transect were processed relative to the fixed control station using Trimble’s proprietary software and precise satellite orbits. This processing procedure has been successfully used to calculate kinematic solutions over baselines of several hundred kilometers during airborne survey work and extended ground traverses in West Antarctica, and should be appropriate for the short baselines used in this study (Reference Spikes, Csatho and WhillansSpikes and others, 2003).
Primary products from the GPS surveys are positions in geodetic coordinates. For further analysis, these coordinates were transformed to Universal Transverse Mercator (UTM) coordinates using the coordinate transformation routine available in ENVI. These coordinates were subsequently transformed to a local Cartesian system with the origin defined by the position of the base station during the first survey, and the x axis approximately in the direction of ice- stream discharge. For both surveys, these positions are given in Table 1.
Whillans Ice Stream survey positions
continued.
Results
Elevation and surface speed
Three of the four pole lines are considered in this study (Fig. 1). The short transect formed by poles B69–B76 is not discussed here because it is very close to transect B58–B68 and results for both lines are essentially the same. Relative elevations and derived velocities are shown in Figure 2. The origin of the local coordinate system is at the base station (SNKE) and the x axis is chosen to be the direction of speed nearest to the center of the ice stream (i.e. at station B18). Crevasse patterns indicated in this and following figures are discussed in a later section.
Location of survey lines across the northern shear margin of Whillans Ice Stream overlain on a portion of the RADARSAT-1 Antarctic Mapping Project (RAMP) mosaic.
Surface elevation contours obtained from gridding GPS measurements at survey stations (black dots; estimated elevation accuracy 0.01 m). Arrows represent ice velocity at survey stations. The dark gray band corresponds to the zone of arcuate crevasses, while the lighter gray band represents a narrow zone of relatively undisturbed surface topography, separating the arcuate crevasses from the chaotic zone.
In contrast to the southern shear margin where elevation contours are mostly parallel to the direction of ice-stream flow (Reference Whillans and van der VeenWhillans and Van der Veen, 2001, fig. 4), the contour map in Figure 2 suggests a significant local component of driving stress acting in the flow direction. At the most inboard station, B18, the local along-flow driving stress is ~22 kPa, or about twice the width-averaged driving stress on the ice stream proper. At the eastern edge of the study area (near transect B58–B68), surface slopes are even steeper, resulting in exceptionally large driving stress (~90 kPa). (Note that elevations along the B69–B76 transect are incorporated into the contour map shown in Figure 2, so the steep gradients in elevation near x = 8 km are not an edge effect of the gridding routine. Of course, due to the lack of stations between the three lines, the exact shape of the contours in these regions is poorly constrained.) The significance of the comparatively large driving stress is not entirely clear, as this survey was conducted mostly in the chaotic zone where the surface topography is highly irregular. Moreover, inspection of elevations from the OSU (Ohio State University) digital elevation model (Reference JezekJezek and others, 1999; Reference Liu, Jezek and LiLiu and others, 1999) shows these slopes to be very localized. Considered over greater distances (~10km), the ice-stream parallel component of surface slope on the ridge is small, resulting in a lower driving stress of ~12 kPa.
Comparison between GPS-derived velocities along the B30–B42 transect and velocities obtained from interferometry.
Velocities measured along the three transects are shown in Figure 3. The speed rapidly increases from near zero on the interstream ridge to several hundred meters per year at the inboard edge of the chaotic zone. Near the center line of Whillans Ice Stream, ice speed is 450ma–1 (Whillans and Van der Veen, 1993), i.e. over a 3–4km narrow zone, ice speed increases from essentially zero to about 75% of the center line speed.
Measured surface speed along the three transects. The dark gray band corresponds to the zone of arcuate crevasses, while the light gray band represents the chaotic zone.
Figure 4 shows a comparison between the velocities along the B30–B42 gridline and those derived from interferometric synthetic aperture radar (InSAR) (Reference Joughin, Tulaczyk, Bindschadler and PriceJoughin and others, 2002). The latter dataset covers much of the Siple Coast ice streams, revealing a complex network of feeder tributaries (Reference JoughinJoughin and others, 1999). The velocity field was derived from 6 days of repeat coverage from the RADARSAT-1 mission in October 1997, which achieved the first complete high-resolution (25 m) radar imaging of the entire Antarctic continent (Reference JezekJezek, 1999, Reference Jezek2003). Generally, there is a good agreement between the two unrelated velocity datasets; the different shapes of the profiles likely reflect that the two profiles are not exactly overlapping.
Shear strain rate and shear stress
The lateral shearing rate and associated shear stress are necessary to determine the importance of lateral drag on the dynamics of the ice stream. At the surface, the dominant velocity gradient is ∂u/∂y, where u represents the along- flow component of velocity, and y the horizontal direction perpendicular to the ice stream. The shear strain rate may then be calculated from
with the corresponding surface shear stress
where B is the temperature-dependent rate factor and n = 3 is the flow-law exponent. For the rate factor, the value B = 700 kPa aш is used, appropriate for the colder upper layers (Reference Whillans and van der VeenWhillans and Van der Veen, 2001). Results of these calculations are shown in Figures 5 and 6.
Derived shear strain rate along the three transects. The dark gray band corresponds to the zone of arcuate crevasses, while the light gray band represents the chaotic zone.
Shear stress along the three transects, calculated from the strain rates shown in Figure 5 using Glen’s flow law with a rate factor B = 700 kPa a1/3. The dark gray band corresponds to the zone of arcuate crevasses, while the light gray band represents the chaotic zone.
The shear strain rate increases rapidly from zero on the interstream ridge to a maximum of 0.06 a–1 at the edge of the chaotic zone. The corresponding near-surface shear stress is ~270 kPa.
Comparison with radar sounding
In 1998, a surface-based ice-penetrating radar system was used to measure surface and bed elevations as well as internal layers across the margins of Whillans and MacAyeal Ice Streams (Reference Raymond, Catania, Nereson and van der VeenRaymond and others, 2006). In the margin of Whillans Ice Stream, radar profiles extend from stations B7 and B68 onto the interstream ridge. For both transects, the position of the shear margin does not appear to be associated with any significant bed topography.
Across the northern shear margin of Whillans Ice Stream, internal layers beneath the heavily crevassed chaotic zone are warped but remain continuous. None of the radar profiles show a change in bed reflectivity or energy returned from internal layers, as would be expected with a transition from a well-lubricated bed under the ice stream to the high- basal-drag ridge (Reference Raymond, Catania, Nereson and van der VeenRaymond and others, 2006). The implication is that the region of maximum shear does not directly overlie a transition in bed morphology. A similar result was obtained by Reference Catania, Conway, Gades, Raymond and EngelhardtCatania and others (2003) who considered bed reflectivity across several relict ice-stream margins. Their measurements show a jump in bed reflectivity from low to high values several kilometers to the ice-stream side of the outermost buried crevasses. They interpret this jump in reflectivity as a transition from frozen to thawed basal conditions. While it is possible that the transition from frozen to thawed conditions has migrated inward since the shutdown of the ice streams studied, as noted by Reference Catania, Conway, Gades, Raymond and EngelhardtCatania and others (2003), the results of Reference Raymond, Catania, Nereson and van der VeenRaymond and others (2006) suggest that the transition from streaming flow to internal deformation on the ridges may be more complex and involve a broader region than suggested by measurements of surface speed only.
Crevasse patterns
Based on aerial photography, Reference Vornberger and WhillansVornberger and Whillans (1986) provided a description of visible features across the margins of Whillans Ice Stream. The outermost part of each margin consists of a band of large, widely spaced (50–100m) arcuate crevasses. Bordering this band on the ice-stream side is the chaotic zone, with intersecting crevasses dissecting the surface. Reference Vornberger and WhillansVornberger and Whillans (1986) presume that major shearing between the ice stream and adjacent ridge occurs within this chaotic zone. The measurements of Echelmeyer and others (Reference Echelmeyer, Harrison, Larsen and Mitchell1994, fig. 2), confirmed this inference. Overlaying the survey stations on the RADARSAT map (Fig. 1) allows the arcuate crevasses and chaotic zone to be mapped relative to these stations. This comparison shows that maximum shearing occurs inside the chaotic zone (Fig. 5).
Vornberger and Whillans (1990) developed a simple model for the development of crevasses by calculating particle trajectories in a specified two-dimensional velocity field. The results indicate that the outermost curved crevasses are initiated by lateral shear and subsequently rotated under the combined action of flow in the down-glacier direction and transverse flow from the ridge into the ice stream. The narrow band of arcuate crevasses starts where the lateral shear stress at the surface is ~Ί30kPa (Fig. 6). Neglecting other components of stress acting in the horizontal plane, the corresponding principal tensile stress is ~90 kPa, oriented at 45°. Such a stress is sufficient to initiate fracturing and crevasse formation on glacial ice (Reference Van der VeenVan der Veen, 1998, Reference Van der Veen1999b).
The RADARSAT image (Fig. 1) indicates that the chaotic zone does not directly abut the band of arcuate crevasses. Instead, there appears to be a narrow zone of relatively undisturbed surface. Such a band has not been reported by previous workers, although the SPOT (Systéme Probatoire pour l’Observation de la Terre) image of the ‘Dragon’ shear margin in Reference Whillans and van der VeenWhillans and Van der Veen (2001) also suggests the presence of a band with a less chaotic surface morphology separating the true chaotic zone from the outermost crevasses. A possible explanation for the crevasse-free zone is that the arcuate crevasses develop at isolated locations along the flank of the shear margin, perhaps in regions where there are transverse increases in shear stress that push the margin of the chaotic zone further towards the interstream ridge. The crevasse-free zone develops as the motion of the arcuate crevasses causes the now relict features to diverge from the edge of the chaotic zone.
Interpretation
Force balance on the ice stream
The role of lateral drag in supporting the driving stress on the ice stream can be readily estimated from the following expression, derived by integrating the depth-averaged force- balance equation over the half-width of the ice stream (Reference Van der VeenVan der Veen, 1999a, section 5.6):
where H represents the ice thickness at the margin, τs the maximum value of Rxy averaged over the ice column and W the half-width of the ice stream. For Whillans Ice Stream, H = 1000 m, τs = 200 kPa and W = 17 km, giving F lat = 12kPa, compared to a width-averaged driving stress of 13 kPa (Reference Stearns, Jezek and van der VeenStearns and others, 2005). Thus, basal drag under the ice stream accounts for only ~1 kPa of resistance to flow, confirming earlier estimates of the relative roles of basal and lateral drags. There is some uncertainty due to possible weakening of ice in the shear margin (Reference Jackson and KambJackson and Kamb, 1997), but this does not affect the general conclusion that lateral drag provides most of the resistance to flow of the ice stream. It should be noted that the value for τs is estimated from the lateral shear stress at the surface, using a depth- averaged value for the rate factor. While this is a sufficiently good first approximation for estimating the role of lateral drag on the motion of the ice stream, this procedure is not fully correct, as suggested by the results discussed below.
Lateral shear concentrated at the ice-stream margins must be transferred to basal drag under the shear margin or interstream ridge. The unsupported driving stress (12 kPa) acts over a half-width of 17 km of the ice stream;thus, outboard of the ice stream must be a region of increased basal drag that supports a total of 2 × 108 Pa m. To investigate this zone of high basal drag and the nature of transfer of stress from shear in horizontal planes to vertical shear, the model developed by Reference Whillans and van der VeenWhillans and Van der Veen (2001) is applied to strain measurements along the central line of survey poles.
Basal drag across the margin
The vertical shearing rate is linked to the vertical shear stress, Rxz , by invoking Glen’s flow law (Reference Van der VeenVan der Veen, 1999a, p. 38):
where
represents the effective shear stress. The assumption is made here that other components of stress do not contribute significantly to the effective stress. Note that all parameters in these expressions represent quantities that vary over the ice thickness and are thus functions of the vertical coordinate z.
Integrating Equation (4) over the full ice thickness gives
where z = h represents the upper ice surface and z = b the ice/substrate interface. To evaluate the integral, the depth variation of the integrand needs to be known. Following Reference Whillans and van der VeenWhillans and Van der Veen (2001), the expression
is adopted, in which τb = Rxz(b) represents basal drag in the direction of ice-stream flow and B b represents the rate factor corresponding to the warmer, basal ice layers, where most of the vertical shearing is concentrated. Equation (7) effectively prescribes how the vertical shear strain rate, and hence the deformational component of ice velocity, varies with depth. For isothermal laminar flow, m = 1 and the vertical shear stress increases linearly with depth from zero at the surface to the maximum value (equal to basal drag) at the bed (Reference Van der VeenVan der Veen, 1999a, section 5.1). In the following analysis, the value m = 2 is used to account for the depth variation in the rate factor. As discussed below in the sensitivity analysis, this parameter choice has little effect on the results. For the deep-ice rate factor, the value B b = 120 kPa aш is adopted, following Reference Whillans and van der VeenWhillans and Van der Veen (2001).
Evaluating the integral in Equation (6) yields an expression from which basal drag can be estimated:
(Reference Whillans and van der VeenWhillans and Van der Veen, 2001). In this expression, S = u(b)/u(h) represents the ratio of sliding speed and surface velocity.
To estimate basal drag from Equation (8), the sliding ratio, S, must be known. Because it is not a priori obvious where the transition from no-sliding to sliding is located, we start with the assumption that S = 0 (no basal motion) along the gridline to arrive at a first guess for basal drag and integrated basal resistance under the ridge. Based on that result, the sliding parameter is estimated via forward modeling to vary in such a way that integrated basal resistance under the ridge equals the width-integrated driving stress on the ice stream.
Figure 7 shows calculated basal drag for the case of no basal slip (S = 0). For ease of discussion, measured surface speeds and derived strain rates are shown in the upper two panels. Along the survey line, basal drag increases towards the ice stream, as was to be expected because basal sliding is not included in this calculation. The bottom panel in Figure 7 shows the integrated excess basal resistance defined as
where τd(y) represents the local driving stress (taken to be 12kPa). The small contribution to the integral from the region not covered by the poles is neglected. F bed represents resistance available to support the nearly frictionless ice stream. The total resistance actually required to support the adjacent ice stream (2 × 108 Pa m) is indicated by the dashed line in the bottom panel. Thus, outboard of station B14 at y = 4.1 km, excess basal resistance is sufficient to oppose the driving stress on the ice stream. Consequently, inboard of this station, basal drag must be significantly smaller: equal to the local driving stress or less. To achieve this, sliding must be included. One possibility is to set the sliding ratio, S, equal to unity inboard of B14. More realistic, perhaps, is to assume this ratio increases more gradually over a distance of one or two ice thicknesses. It may be noted that if the width of this transition region is too great, basal drag remains too small for Ғbed to reach the value required for balancing the driving stress on the ice stream. In other words, the transition from all sliding on the ice-stream side to no basal slip on the interstream ridge must occur over a narrow zone. Furthermore, the general shape of the functional change of the sliding ratio across the shear margin must be similar to that adopted here. For example, allowing a small component of deformational velocity on the ice stream (that is, S < 1) results in too much integrated basal resistance, while shifting the transition from no-sliding to sliding more towards the interstream ridge leads to insufficient basal resistance to support the main body of the ice stream. The prescribed sliding ratio is shown in the top panel of Figure 8; corresponding basal drag and integrated basal resistance are shown in the second and third panels. After reaching a maximum slightly less than 100 kPa, basal drag decreases rapidly towards the ice stream, due to the transition from nosliding to sliding. On the ice-stream side of the survey line, basal drag is negligible. With this result, the transfer of stress can be investigated further.
Results of the force-balance calculation for the central line of survey poles for the case of no basal sliding. The upper panels show measured surface speed and shear strain rate. The third panel shows basal drag calculated from Equation (8). The bottom panel shows integrated basal resistance as defined in Equation (9); the horizontal dashed line corresponds to net basal resistance required to balance the driving stress on the ice stream.
Results of the force-balance calculation for the central line of poles with inclusion of basal sliding. The top panel shows the prescribed sliding ratio with the transition from no sliding (S = 0) to all sliding (S = 1) occurring over a zone of ~2 km wide. The second panel shows calculated basal drag and the third panel the integrated basal resistance. The bottom panel shows the stress guide calculated from Equation (13).
Transmission of stress
The balance of forces integrated over the ice thickness reduces to (Reference Whillans and van der VeenWhillans and Van der Veen, 2001)
in which the y axis is perpendicular to the shear margin. Following Reference Whillans and van der VeenWhillans and Van der Veen (2001), a stress guide, φ(y), is introduced:
This stress guide accounts for the depth variation in lateral shear stress and represents the fractional depth over which the surface stress acts. For the case that this stress is constant over the full ice thickness, φ(y) = 1. With this parameterization, the balance equation becomes
which can be integrated to find an expression for the stress guide:
Here the assumption is made that the stress-guide factor is zero at the outboard end of the survey line (y = 0), thus neglecting any contribution to the integral from excess basal drag outboard of the survey line. In Equation (13), Rxy represents the lateral shear stress at the surface shown in the middle panel in Figure 6.
The calculated stress guide is shown in the bottom panel of Figure 8. On the ridge side, the stress guide is small, as was to be expected since the integrated basal resistance is small there. Going towards the ice stream, the stress guide increases to 0.94, indicating that the surface lateral shear acts over 94% of the ice thickness. As basal drag drops to zero and the integrated basal resistance starts to decrease, the stress guide begins to decrease also. The survey line does not extend far enough into the ice stream to evaluate the stress guide on the ice stream proper, but its value there should be equal to the ratio of the depth-averaged rate factor to the surface value of the rate factor (~ 0.77).
Variation of stresses across the shear margin. The top panel shows calculated basal resistance; the middle panel shows the lateral shear stress at the surface evaluated from measured surface velocities, while the third panel shows the inferred depth-averaged lateral shear stress obtained by multiplying the surface value of the lateral shear stress and the stress guide.
The variation of the stress guide along the survey line is essentially different from that presented in Whillans and Van der Veen (Reference Whillans and van der Veen2001, fig. 5e), which shows a decrease in stress guide from unity on the interstream ridge to about 0.4 at the ice-stream end of the Dragon grid. That result is wrong and our recalculation of the stress guide for the Dragon grid shows a pattern similar to that shown in Figure 8. A programming error caused the erroneous result presented by Reference Whillans and van der VeenWhillans and Van der Veen (2001). Apart from the discussion in Reference Whillans and van der VeenWhillans and Van der Veen (2001) pertaining to the stress guide and its variation across the shear margin, the conclusions reached in that paper are still supported by their analysis, and, in particular, calculated basal drag (Reference Whillans and van der VeenWhillans and Van der Veen, 2001, fig. 5f) is correct inasmuch as the adopted model is correct, of course.
The maximum in basal drag is ~1 km inboard of where lateral shearing at the surface reaches its maximum (Fig. 9). While this might appear to be violating the horizontal balance equation (10), it should be kept in mind that this balance equation calls for the depth-integrated lateral shear stress. This average is obtained by multiplying the shear stress at the surface with the inferred stress guide and is shown in the third panel in Figure 9. As should be the case, the depth-averaged shear stress is maximum on the ice- stream side of the shear margin, where basal drag vanishes. Thus, if the margin of the ice stream is defined as the location where basal drag becomes important, that position does not coincide with where shearing at the surface is greatest. Additionally, our results indicate that the depth- averaged shear stress cannot be simply estimated from velocity gradients measured at the surface and applying an appropriate depth-averaged rate factor. This results from the contribution of both lateral and vertical shear to the effective stress in the marginal zone.
Sensitivity of model results
Two sources of uncertainty contribute to uncertainty in calculated basal drag, namely measurement errors and unknown parameter values. In Equation (8), from which basal drag is calculated, the measured quantities are the surface velocity, u(h), and ice thickness, H, while the deep rate factor, Bb, and the exponent, m, represent parameters whose values are uncertain.
How measurement errors affect calculated basal drag can be evaluated from the standard formula for error propagation which yields for the relative error
The estimated error in GPS-determined positions is 0.01 m or less (Reference Spikes, Csatho and WhillansSpikes and others, 2003), giving an error in velocity of about 0.015 ma–1. On the ridge side of the survey line, the surface speed is 1.5m a–1 and the relative error in velocity is 10%; on the ice-stream side, this relative error is significantly smaller. The corresponding relative error in calculated basal drag is thus ~3% on the ridge side of the survey line and much smaller on the other end of the line (using n = 3). For the ice thickness, we estimate a relative error of 10%, giving a relative error in basal drag of 3%. Combining the two error estimates, the relative error in calculated basal drag is 6% or less. Compared to the effects of parameter uncertainties, this error is inconsequential.
Figure 10 illustrates the sensitivity of basal drag, τ b, and depth-averaged lateral shear stress, τ S, to the values of the deep rate factor, Bb, and exponent, m. Of these two parameters, uncertainty in the deep rate factor has the greatest impact on calculated basal drag and lateral shear stress. A small value (curves labeled ‘1’) produces smaller values of basal drag and hence a wider region where drag at the glacier base exceeds the local driving stress. Increasing the rate factor has the opposite effect and leads to larger values in basal drag but over a narrower region (curves ‘2’). As noted previously by Reference Whillans and van der VeenWhillans and Van der Veen (2001), the choice for the value of the exponent m in Equation (7) has only a minor effect on the model results (curves ‘3’–’5’).
Sensitivity of calculated basal drag (upper panel) and depth-averaged lateral shear stress (lower panel) to parameter values. Labeled curves correspond to: 1: Bb = 80kPaa1/3, m = 2; 2: B b = 200 kPa a1/3, m = 2; 3: B b = 120kPaa1/3, m = 1; 4: B b = 120kPaa1/3, m = 3; 5: B b = 120kPaa1/3, m = 4.
Meltwater production
The narrow zone of elevated basal drag suggests that meltwater production may be important under the ridge. Assuming the basal ice is at the pressure-melting point, the rate of melting or freezing is determined by the balance between the geothermal heat flux, frictional heating associated with sliding over the bed and upward heat conduction. That is (Reference Lingle, Brown, Van der Veen and Oerle mansLingle and Brown, 1987),
where L i = 333.5 kJ kg–1 represents the latent heat of fusion, ρ = 917 kg m–3 is the density of ice and k = 2.1 W(m K)–1 is the thermal conductivity of ice. For the geothermal heat flux, the value G = 0.06Wm–2 (Reference Whillans and BindschadlerWhillans and Bindschadler, 1988) is adopted. No temperature profiles to the bed are available for the interstream ridge, and the high value measured on the adjacent ice stream, ∂T/дz(b) = 0.04 Km–1 (Reference Engelhardt, Humphrey, Kamb and FahnestockEngelhardt and others, 1990) is used here. This high value may overestimate upward heat conduction and thus underestimate calculated melt rates, but the effect is minor given that frictional heating (the third term in parentheses in Equation (15)) is dominant in the region considered.
Frictional heating across the shear margin can be estimated from the results shown in Figure 8 recognizing that the basal velocity equals the surface speed multiplied by the sliding parameter. The result of this calculation is shown in Figure 11. Under the ridge where sliding is zero, the basal freeze-on rate is ~2mma–1 (this could be somewhat less if a smaller value for the basal temperature gradient is used). In the region of onset of basal sliding, the combination of sliding and high basal drag results in melting rates reaching 40 mm a–1. Going further towards the ice stream, the melt rate decreases as basal drag decreases rapidly. Near the center of the ice stream, the estimated melt rate is ~2mma–1 (Reference Stearns, Jezek and van der VeenStearns and others, 2005). Thus, melting under the chaotic zone is an order of magnitude greater than melt rates under the main body of the ice stream. We propose that this concentrated melting is the process by which shear margins migrate outward.
Estimated basal melt rate under the middle transect. Negative melt rates correspond to basal freeze-on.
The estimated void ratio of till beneath Whillans Ice Stream is 0.6–0.7 (Reference KambKamb, 1991; Reference Tulaczyk, Kamb and EngelhardtTulaczyk and others, 2000a). Assuming a void ratio of 0.65, a 5 m thick sediment layer contains the equivalent of 2 m of water, corresponding to 2000 years of basal melting at a rate of 1 mm a–1 (Tulaczyk and others, 2000b). Under the chaotic zone, the larger basal melt rate may allow basal sediments to become water- saturated more rapidly, resulting in a widening of the ice stream. With a melt rate of 40 mm a–1, a 5 m thick sediment layer would be fully water-saturated after ~50 years, resulting in lowering of basal resistance and outward migration of the region of large basal drag.
Discussion
While the model developed by Reference Whillans and van der VeenWhillans and Van der Veen (2001) is based on a number of simplifying assumptions allowing a solution to be derived, it does incorporate the essential physical processes believed to be important in the transfer from shearing in horizontal planes to vertical shear. Perhaps the most contentious model assumption is the prescribed vertical variation in vertical shear (Equation (7)). In essence, this corresponds to prescribing the depth variation of the effective viscosity of ice. Other functional forms could be adopted, such as constant effective ice viscosity (Reference SchoofSchoof, 2004), or exponentially increasing with height above the bed (Reference Guðmundsson, Raymond and BindschadlerGuðmundsson and others, 1998). However, qualitatively, the results presented in Figures 7–9 are not importantly affected by the choice for the exponent, m, suggesting that the precise nature of the depth variation of effective viscosity is of secondary importance for calculated basal drag and lateral shear stress. The value of the deep rate factor, B b, affects the magnitude of calculated basal drag and the width of the region of large drag, but qualitatively, results are unaffected (Fig. 10). Thus, the results obtained here appear to be robust and can be used to qualitatively describe the stress transition in this shear margin. It should be pointed out that the model as formulated applied here is not a true inversion problem as implied by Schoof (Reference Schoof2004, p.217). Rather, measurements of surface speed are used to infer basal drag by integrating the vertical shear strain rate over the entire ice thickness, adopting Glen’s flow law to link the strain rate to the vertical shear stress. By introducing the sliding ratio, the difficulty of selecting an appropriate sliding relation is circumvented here.
On the ice stream, resistance from the bed is near zero and the driving stress is balanced by gradients in lateral shearing. As the margins are approached, basal drag increases rapidly and basal conditions transition from well lubricated to no basal slip. Basal drag decreases further onto the ridge. For the northern shear margin, the zone of elevated basal drag is ~4km wide. This width is not significantly affected by the prescribed transition from sliding to no-sliding. Over this region, lateral drag decreases towards the ice ridge and becomes increasingly limited to acting over the upper layers only. The maximum in basal drag lies inboard of where lateral shearing at the surface is maximum.
The present model does not include enough constraints to yield a unique solution. In particular, the transition from sliding to no-sliding may be different from that assumed here (Fig. 8, top panel). Nevertheless, some constraints are imposed by the requirement that basal drag across the margin must be large enough to support the adjacent ice stream. In the present model, this requires a rather narrow zone over which the sliding fraction increases from zero on the ridge side, to unity on the ice-stream side. In this narrow zone, a significant amount of meltwater is produced.
Clearly, given the inferred large value of basal resistance under the chaotic zone, this region cannot be underlain by a weak sedimentary layer as is the case under the ice stream. At the same time, some lubricating water layer must be present to allow for basal sliding. A thin water layer (several mm or less) separating the subglacial sediments from the ice is consistent with the low-reflectivity returns observed by Reference Raymond, Catania, Nereson and van der VeenRaymond and others (2006), but it appears that most of the remaining meltwater produced in the region of elevated basal drag is diverted elsewhere. Some of this meltwater may be frozen on outboard of this region, but the bulk of the water produced must either percolate downward into the sediments, or be diverted under the ice stream. Wholesale storage in subglacial sediments appears to be incompatible with the results of the radar profiling as well as the requirement of a sufficiently strong bed to support large basal drag, as this would rapidly produce a weak water- saturated sediment layer of substantial thickness. More likely is the possibility that much of the meltwater produced under the shear margin is rerouted towards the center of the ice stream, providing additional lubrication and aiding in maintaining water-saturated conditions under the ice stream required for fast flow with near-zero basal resistance.
The observation that only a fraction of the locally produced meltwater is used to mobilize subglacial sediments indicates that the available void space in these subglacial deposits must be small. As noted above, the void ratio of till under Whillans Ice Stream ranges from 0.6 to 0.7 (Reference KambKamb, 1991; Reference Tulaczyk, Kamb and EngelhardtTulaczyk and others, 2000a), which would accommodate a substantial influx of meltwater. Since this does not happen, much of the void space must be filled with frozen water, and mobilizing this sediment is achieved primarily by melting of the ice occupying these voids. Given that this requires transferal of heat generated at the base of the ice to within the sediments, mobilization may be expected to be less rapid than for the case of saturating basal sediments with a significant void ratio. A further inference is that the shear margins may contribute significantly to the water balance under the ice stream. If all meltwater produced under the chaotic zone is distributed evenly under the ice stream, this input would correspond to ~2 mm a–1 of additional water (per unit area), or about twice the estimated melt rate under the ice stream.
Concluding Remarks
The present study extends the previous work of Reference Whillans and van der VeenWhillans and Van der Veen (2001) in evaluating the transmission of stress from a fast-moving ice stream to the adjacent near-stagnant interstream ridge. The survey line considered in this study covers the region from vanishing basal drag under the ice stream to small basal drag on the interstream ridge. Under the chaotic zone is a narrow region (~-4km wide) where drag at the glacier base is large; this excess basal drag essentially supports the driving stress on the ice stream. Further, the transition from ice motion dominated by basal sliding to internal deformation only is restricted to a narrow band, 1–2 km wide.
As concluded by Reference Raymond, Catania, Nereson and van der VeenRaymond and others (2006), lubrication under the shear margin must have different characteristics from that under the main body of the ice stream. Some thin layer of dilated till may be present under the margin, but most of the locally produced meltwater must be diverted under the ice stream. Our results suggest that the amount of meltwater produced under the chaotic zone is about twice the amount produced under the ice stream. This suggests that the shear margins may be controlling the motion of the ice stream not only by providing support for the gravitational driving stress, but also by facilitating weakening of sediment under the ice stream. If this view is correct, the implication is that studies of ice-stream stability that consider only the coupling between glacier speed, local basal drag and meltwater production under the ice stream may not fully incorporate all of the controlling processes.
Acknowledgements
This study would not have been initiated without the enthusiastic efforts of the late I. Whillans. We are indebted to G. Hamilton, B. Spikes, P. Braddock and S. Norman for conducting the surveys, and to C. Raymond for encouraging completion of the manuscript. Deployment of strain grids and initial data processing was supported by the US National Science Foundation (NSF) through grant OPP- 9615127. Subsequent data analysis and interpretation was made possible by support from NASA through NAG5-11383 and by NSF through the Center for Remote Sensing of Ice Sheets. This is Byrd Polar Research Center contribution No. C-1346.
