1. Introduction
Changes in the resistance to motion of ice streams over their bases is a primary concern with regard to accelerated discharge of ice from the West Antarctic ice sheet (WAIS) and the resulting sea-level rise. The conditions at the base of ice streams are best examined directly with access by boreholes (Reference Kamb, Alley and BindschadlerKamb, 2001). Densely spaced information along profile paths can also be found from remote sensing using seismic (Reference Blankenship, Bentley, Rooney and AlleyBlankenship and others,1987; Reference Anandakrishnan, Blankenship, Alley and StoffaAnandakrishnan and others, 1998) and radar (Reference Bentley, Lord and LiuBentley and others, 1998; Reference Gades, Raymond, Conway and JacobelGades and others, 2000) data. Inversion of the disturbance of surface topography and velocity field arising from bumps and sticky spots on the bed provides another line of attack that has the potential to provide broad map-plane spatial coverage utilizing satellite remote sensing (Reference MacAyeal, Bindschadler and ScambosMacAyeal and others,1995). The spatial pattern is key to examining the relationship between basal topography and local lubrication, which is still unknown for any ice stream.
Small-scale undulations on the surface of ice streams and flow variations associated with them are pervasive and characteristic of streaming flow. An example of an ice-stream surface is shown in Figure 1. In this paper we use the theory of glacier flow perturbations caused by small variations in bed topography and/or basal lubrication (Reference GudmundssonGudmundsson, 2003) to investigate the sources of these surface features and their implication for basal conditions retarding ice-stream motion. The theory has been used to show that flow stripes on ice streams can form from flow over localized bedrock bumps or sticky spots (Reference Gudmundsson, Raymond and BindschadlerGudmundsson and others,1998).
Landsat image of a 64 × 32 km2 area on Ice Stream E, West Antarctica. Surface undulations of various wavelengths are made visible by oblique solar illumination from the top of the image. Direction of motion is westward, approximately from top to bottom of the figure. Coordinates of the corners of the outlined box are: lower left: 79.9890° S, 144.3307° W; upper left: 79.6543° S, 140.6127° W; upper right: 79.9958° S, 139.5263° W; lower right: 80.3421° S, 143.3374° W.
The goal of this paper is to use maps of the horizontal variation in x and y of surface elevation s(x, y) and components of horizontal velocity u(x, y) and v(x, y) to infer bed elevation b(x, y) and basal lubrication c(x, y) under an ice stream where c(x, y) relates motion tangential to the bed to the local shear stress. In principle, bed topography b(x, y) could be measured with radio-echo sounding, which would allow a more constrained computation of lubrication c(x, y) and a more secure examination of any correlation between b and c. In practice, map-plane measurements of b are not likely to be realized in the near future at a very fine spatial scale over large areas. For this reason we attempt to infer both b and c.
The primary assumption of the theory is that b(x, y) and c(x, y) disturbances to flow and associated topography are small compared to mean conditions, thus allowing analysis using perturbation theory (Reference Hutter, Legerer and SpringHutter and others, 1981). In addition, the ice is assumed to be linearly viscous with an exponential vertical variation in viscosity, which can partly account for effects from vertical variation of temperature and non-linearity of the ice flow law.
The theory can be formulated in terms of transfer functions that relate bed topography, b, and basal lubrication, c, to surface topography, s, and velocity components, u, v and w (Reference GudmundssonGudmundsson, 2003). It is an extension of earlier work by Reference BuddBudd (1970), Reference Hutter, Legerer and SpringHutter and others (1981), Reference HutterHutter (1983),Reference Whillans and JohnsenWhillans and Johnsen (1983), Reference Balise and RaymondBalise and Raymond (1985), Reference ReehReeh (1987) and Reference JóhannessonJóhannesson (1992) to include full three-dimensionality in the flow and time dependence in the basal boundary condition. The theory shows that effects on the ice flow from localized perturbations in the basal conditions (b or c) are strongly attenuated toward the surface. This attenuation complicates the inverse calculation of bed conditions from surface measurements (Reference Bahr, Pfeffer and MeierBahr and others, 1994). The bed-to-surface attenuation is smaller when the mean basal motion is large compared to motion from internal deformation. Thus, ice streams that move dominantly by motion over their beds are amongst locations where this kind of inversion has some chance of success.
To test inversion procedures, we specify various patterns of b(x, y) and c(x, y) and generate corresponding synthetic data for s(x, y), u(x, y) and v(x, y) with added noise to simulate errors.We find that successful inversion should be possible for a useful range of spatial scales in the fast-flowing ice streams.
Based on the experience gained from inversion of synthetic data, we investigate the covariation of b and c in a local area of Ice Stream E,West Antarctica, where measurements of s, u and v are available. This area lies within a larger area examined earlier by Reference MacAyeal, Bindschadler and ScambosMacAyeal and others (1995). Our analysis differs from theirs in several ways. The fully three-dimensional treatment of the flow enables us to probe for small-spatial-scale features that cannot be taken into account by the vertically integrated treatment used by Reference MacAyeal, Bindschadler and ScambosMacAyeal and others (1995). We input detailed information on the surface topography derived from recent altimetry and photoclinometry (Reference Bindschadler, Vornberger, Blankenship, Scambos and JacobelBindschadler and others, 1996), whereas Reference MacAyeal, Bindschadler and ScambosMacAyeal and others (1995) used only velocity pattern as input and inverted for the surface topography. We also use updated surface velocity measurements (Reference JoughinJoughin and others, 1999). The main disadvantage of our perturbation-based method is that it is restricted to small variations from mean bed conditions, and boundary conditions on the edges of the domain cannot be arbitrarily specified. Thus, we cannot examine questions such as the relative role of side drag and longitudinal forces associated with large amplitude and spatial-scale variations in bed conditions; these questions are better attacked by the Reference MacAyeal, Bindschadler and ScambosMacAyeal and others (1995) approach.
We have found that spatially correlated variations in basal topography and lubrication are required to explain the combined features of surface elevation and velocity in Ice Stream E.
2. Description of the Forward Theory
The forward theory is described in detail by Reference GudmundssonGudmundsson (2003). It solves equations for creeping flow of an incompressible Newtonian fluid for deviations in velocity field and associated surface evolution from mean steady conditions as caused by small, possibly time-dependent changes in the basal topography and lubrication. In this paper, we assume time-independent basal conditions on a bed z = b(x, y), where the vertical coordinate is z. The upper surface elevation is taken as z = s(x, y). The fluid viscosity is allowed to vary exponentially with depth s – z beneath the surface given as
where η s is the viscosity at the upper surface and μ scales the vertical variation of viscosity.
Boundary conditions on the bed are specified by the equation:
where n is the local normal to the bed, T(n) is the stress vector acting on the bed, and t is the local tangent to the bed in the plane common to T and n. This proscribes no motion normal to the bed (u n = 0) and bed-tangential motion proportional to and in the direction of the local basal shear stress. The proportionality factor c(x, y) describes the local lubrication of the bed. Corresponding boundary conditions on the upper surface s(x, y) are no motion normal to the local surface and no shear stress tangential to it. Thus, accumulation, ablation and time evolution of the upper surface are not taken into account.
The bed height and its lubrication are partitioned as b(x, y) = b 0(x, y) + b 1(x, y) and by c(x, y) = c 0 + c 1(x, y). b 0(x, y) is a mean bed represented by a plane of slope α0.We take the x axis to be down this slope, so b 0(x, y) = b 0(x). c 0 is a constant mean lubrication.When b 1 and c 1 are zero and the ice thickness is uniform, the problem reduces to a uni-formly sloped slab that moves by shearing over its thickness at a rate
and moves uniformly over its bed at a rate
where h 0 = s 0(x) – b 0(x) is the contant thickness of the slab. This solution is regarded as a mean reference condition defined by h 0, α0, c 0, ηs and μ.When b 1 ≪ h 0 and c 1 ≪ c 0 are specified, perturbation analysis of the full three-dimensional field equations for incompressible, creeping flow and the boundary conditions gives solutions for the full distributions of all components of velocity (u 1(x, y, z), v 1(x, y, z), w 1(x, y, z)) and strain rate, stress and the height of the surface (s 1, x, y)). Setting the perturbation equations in non-dimensional form (Reference GudmundssonGudmundsson, 2003) shows that the characteristics of the solutions are determined by three dimensionless parameters: the slope 0, a measure of the vertical variation of viscosity
and the slip ratio
Definition of a mean reference condition in an actual ice stream is straightforward in principle. h 0 and α0 can be found from sparse radio-echo sounding measurements and elevation models from satellite altimetry. In manyWAIS ice streams, u d ≪ u b so c 0 can be estimated directly from the surface velocity.η s and μ can be chosen to best match the vertical variation of effective viscosity associated with the non-linear ice-flow law for the mean flow Ξ and Ω the temperature distribution. and can be determined in a similar manner.
We focus on the connection of perturbations on the upper surface s 1(x, y), u 1(x, y; s) and v 1(x, y; s), which are accessible to measurement, with those on the bed b 1(x, y) and c 1(x, y), which are not accessible to direct measurement. For notational compactness we drop the subscript 1, and s(x, y), u(x, y), v(x, y), b(x, y) and c(x, y) will henceforth represent the perturbations.
Description of the general results and much of the application of the theory is done in wavenumber space. Capital letters S(k; l), B(k; l), C(k; l), U(k; l) and V (k; l) denote the Fourier transform 
of surface elevation s(x, y), bed elevation b(x, y), basal lubrication c(x, y), longitudinal u(x, y) and lateral v(x, y) velocity, respectively. For example 
, which we will write as S.
The solutions for surface topography S and velocity components U and V are then represented in wavenumber space as
where the transfer functions 
depend on α
0, Ξ and Ω Figure 2 shows transfer functions as a function of longitudinal (k) and lateral (l) wavenumber for the case where α = 0.11°,Ω = 100 and Ξ = 2.5, which are values appropriate for the case of Ice Stream E considered below. A number of narrow peaks can be seen for k
2 + l
2 ≈ 1 in the velocity transfer functions (Fig. 2). These are presumably related to the introduction of numerical errors, as a 6 × 6 matrix must be inverted for each wavenumber to solve for integration constants entering the transfer function. A singular-value decomposition of this matrix shows that for certain combinations of wavenumbers, mean surface slope and vertical variation of viscosity, its singular values become small and the matrix almost singular. These spurious peaks are narrow and isolated, and although it is not clear why they appear, they do not affect the results of the inversion.
The amplitude of the transfer functions 
as a function of longitudinal (k) and lateral (l) wavenumber for Ω = 102 Ξ = 2.5, and α0 = 0:11°. Left column: transfer function for bed topography B to surface topography S 
, longitudinalvelocity 
and lateral velocity 
. Right column: lubrication C to surface topography 
, longitudinal velocity 
and lateral velocity 
. Units in plots are: 
, 
, where [C] = u0/(ρgα0h0). Note that the grayscales are calibrated using log10 for 
and 
and that l = 0 and K → 0, 
, although not seen on the grid here.
3. The Inverse Problem
We use a least-squares inversion of Equations (7–9) to find B and C that give the best consistency with specified S, U and V. Error-weighted residuals are defined as
where Σ
S
, Σ
U
and Σ
V
are the rms observation errors for S, U and V, which are assumed to be independent. The error weightings non-dimensionalize the expected residual to about 1 in each case. Following Reference MenkeMenke (1989) minimization of 
gives
where
E is an error matrix, here assumed to be diagonal with 
, along the diagonal, and superscript H denotes the Hermetian transpose,
For compactness, we define
Then the solutions for the bed B and lubrication C are
where
The solution for B and C involves dividing by D, which is pictured in Figure 3. The inversion is problematic for those (k, l) where D is zero or small. This is the case for all large wavenumbers (wavelengths shorter than about 2h 0) and k = 0 (some combinations of cross-flow variations in b(.,y) and c(.,y) which we will call longitudinal ridges and stripes). With regard to short-wavelength bed features, such fluctuations in the basal boundary condition will be accom-modated by flow disturbances near the bed that do not influence the surface in a measurable way. It is interesting that longitudinal ridges and stripes are not transferred to the surface at all scales, and there will be a problem even for those inversion procedures that are concerned only with longer scale (e.g. Reference MacAyeal, Bindschadler and ScambosMacAyeal and others,1995). An easily understood example is indeterminancy of the combination of bed topography and slip condition across a valley glacier based on measurement of down-glacier velocity across its transversely flat upper surface. A shallower bed can be compensated by a more lubricated bed (or visa versa) to produce an equivalent fit to the down-valley speed on the surface. More generally, the problem is that lateral variations in b or c do not generate any change in surface topography or lateral velocity, and by arranging the phase and amplitude of b and c perturbations, the total surface velocity of the combination cancels. Then, there is no signal generated by b and c at the surface at all.These kinds of bed features that do not transfer to the surface are in the null space of the forward problem. Such features can be added to any inversion for b and c without affecting the fit to surface data about s, u and v. They display the non-uniqueness of the inversion. Correspondingly, it is not possible to invert for them on the basis of surface data.
D(k,l) defined in Equation (19) for Ω = 102, Ξ = 2:5 and α0 = 0:11° with Σs = 2 m, Σu = Σv = 5 m a−1. The grayscale gives D in units of log10(h0C)2 = (ρgα0/h0)2.
To avoid this ill conditioning of the inverse, we use a truncated version of D as a filter F to remove variations in the problematic wavenumber ranges from the input data S, U and V. The filter is
It allows all wavelengths where D is larger than P to go through, but filters out other wavelengths. P is therefore an adjustable trade-off parameter. Smaller values of P give more detail, but may unjustifiably over-fit the surface data in view of errors. Larger values of P under-fit the data and may leave out features actually demanded by the data. We describe P = max(|D(k, l)|)ΩP, where the power p ≤ 0 is varied to adjust P.
Filtering is essential to avoid blow-up of errors in the data or the model, that could be interpreted erroneously as being generated by huge features on the bed. In practice, measurements of S, U and V may contain other noise that physically could have been generated from the bed and would not be removed by the filter F. In addition to measurement errors, such noise could arise from interference outside the model region transmitted or advected into it, for example flow stripes (Reference Gudmundsson, Raymond and BindschadlerGudmundsson and others, 1998} or non-steady surface topography.
The capabilities and limitations of the inversion procedure can be explored using synthetic data. Our procedure is to use the forward theory (Equations (7–9)), to calculate s, u and v due to prescribed distributions of b and c. We combine the patterns generated by b and c and add random noise to each of the resulting total s, u and v signals. The noise amplitude is set at 50\% of the maximum amplitude of each (combined) signal. This is similar to the noise level in data from the surface of Ice Stream E that we examine Ice Stream E (Reference Bindschadler, Vornberger, Blankenship, Scambos and JacobelBindschadler and others, 1996) at locations below. We then invert for b and c using the combined s, u and v distributions with the added noise. A number of these synthetic inversions have been done. Only two examples are given here that illustrate the limitations and successes.
The checkerboard pattern shown in Figure 4 combines both lateral and longitudinal variations in b and c and provides an example where some information about the bed is transferred to the surface and some is not. The substantial lateral variation at the bed in both b and c appears on the surface only as a small lateral variation in u. Correspondingly, the inverted bi and ci (far right in Fig. 4) captures only the longitudinal variations. The filtering procedure (Equation (20)) suppresses production of large unrealistic lateral variations in b and c that would otherwise be generated by the noise. The small high-frequency lateral variations are a result of noise that has not been fully taken into account by the filtering. This example illustrates the unfortunate reality that the inversion is blind to certain features of b and c, but also demonstrates that the filtering stabilizes the inversion against the corresponding indeterminancy.
Surface topography (s) and velocity (u and v) generated by superposition of longitudinal and lateral variations in b and c. The bed and lubrication patterns returned by the inversion of the surface “data”, added together in equal proportions with added noise in each signal (s, u and v) equal to 50% of its maximum amplitude, are shown on the far right as b i and c i.
Gaussian shape perturbations in b and c provide an example of localized bed features with both longitudinal and lateral spatial variation. Figure 5 shows input and inverted b and c. Figure 6 shows a longitudinal section along y = 0. Remember that, in this and the previous example, the signals (s, u, v) generated by b and c are added together and degraded by noise before the inversion. This example shows that this inversion method can separate effects from b and c perturbations and accurately recover both longitudinal and lateral features of their spatial patterns including their relative positions.
Surface topography (s) and velocity (u and v) generated by Gaussian peak distributions for b and c. Note that the b and c patterns are slightly offset in the along-flow direction, and normalized to have a zero mean amplitude.The far right panels show the inferred bed and lubrication patterns, b i and c i, using the surface data in equal proportions, and adding noise that is 50% of the maximum amplitude of the signal (s, u and v).
Along-flow cross-section through y = 0 of the surface s, bed topography b and lubrication c, generated by the Gaussian bed and lubrication shown in Figure 5.The thin linesare the input. The dashed line in the top plot is the surface with added noise.The thick lines are the surface, bed elevation and lubrication inverted using the noise-degraded surface signals.
4. Application to Ice Stream E
Figure 7 shows surface elevation s(x, y) and velocity u(x, y), v(x, y) that have been extracted from satellite imagery Ice Stream E (Reference Bindschadler, Vornberger, Blankenship, Scambos and JacobelBindschadler and others, 1996) at locations covering the area shown by the box in Figure 1. We invert for the basal conditions within this box based on the data in Figure 7. All of s, u and v have a linear feature running from the lower left corner to the upper right corner. Most of our discussion will focus on the features at the bed that are responsible for these most prominent surface features.
Ice Stream E topography $s$, with the mean elevation and gradient removed, and velocity along flow (u) and lateral to flow(v).
At this location h 0 = 1.1 km and α0 = 0.11°. The mean longitudinal speed is 370 m a−1. There is a mean lateral shear strain rate of 0.04 a−1. The deformational velocity u d is estimated to be 5.4 m a−1 using the standard temperature-dependent, non-linear flow law for ice (Reference PatersonPaterson,1994, p. 97). Temperature was estimated following the theory of Reference RobinRobin (1955), assuming steady state with a surface temperature of 247 K (Engelhardt and others, :1990}, a bed temperature at the pressure-melting point, a surface accumulation rate of 0.1 m a−1 and a linear vertical variation of vertical velocity. This estimate yields a temperature profile in reasonable agreement with temperature measured through Whillans Ice Stream (Reference Engelhardt, Humphrey, Kamb and FahnestockEngelhardt and others, 1990). The effect of the mean lateral shear strain rate on viscosity and an 18% reduction in the surface-parallel shear stress by side drag (Reference RaymondRaymond, 2000) were taken into account. The implied dimensionless value for Ω is then 370/5.4 = 68 (i.e. of order 10−2 which is the value used in the inversion described below). The effective viscosity decreases with depth beneath the surface in an approximately, although not exactly, exponential manner corresponding to dimensionless $\Xi = \mu h_0$ = 2.5. The value of U d calculated as above is reproduced by Equation (5) with this value of and η s = 3:7 × 1014 Pa s.
Since the transfer-function theory analysis in the wavenumber domain treats the boundaries as periodic, discontinuities are introduced at the boundaries, and the inversion will not be reliable near those boundaries. To alleviate this problem, we first linearly taper the data to zero over 10 km at up- and downstream edges and over 5 km on the sides. Then we remove the mean elevation and gradient (slope α0) from the surface topography data. We also remove mean longitudinal velocity and the lateral velocity gradient from u.
Error levels (Equations (10–12)) are estimated to be Σ
s
= 2.2 m, Σ
u
= 10 m a−1 and Σ
v
= 20 m a−1. Velocity data are from feature tracking, and errors are expected to be uncorrelated from one location to another independent of distance. Thus, we expect a flat error spectrum for u and v. On the other hand, s is determined by integration of slopes, which would produce a red error spectrum. We nevertheless use a flat error spectrum for s, since we are suppressing short-scale features on the surface in any case because of the the transfer characteristics. The error-weighted root-mean-square (rms) residual is 
Here 
, where N is the number of gridpoints and x stands for s, u or v.
Figure 8 shows the scaled residuals that result from the inversion for different cut-off for the filtering as described by P (Equation (20)). A value of p = 2 gives residuals at a level consistent with estimated errors and we use this value for our preferred inversion.
Residual r2 for inversion as a function of the power p in the filter defined by Equation (20) and following text.
Figure 9 shows the inferred bed topography (b) and basal lubrication (c) using these data. The overall agreement is consistent with expected errors (Table 1, first row). The error-weighted rms residual is r 2 = 0.404. Table 1 shows the rx values for each of the surface variables, using the error levels suggested above,Σ s = 2.2 m, Σ u = 10 m a−1 and Σv = 20 m a−1.
Preferred inversion for b and c using s, u and v (Fig. 7) with h0 = 1100 m, α 0 = 0.11°,Ω = 100 and Ξ = 2.5.
The rx and r2 values for inversions for both b and c, only b or only c (see text)
The main bed feature identified in the inversion is a prominent linear zone involving both b(x, y) and c(x, y) that runs diagonally across the region from x/h 0 ∼ −21 on the left side (y/h 0 ∼ −15) to x/h 0 ∼ 15 onthe right side (y/h 0 ∼ + 15). It is a deep trough in the topography and a low in the lubrication that approximately coincide. It lies slightly upflow from the linear surface feature in s, u and v apparent in Figure 7.
Outside of this linear zone, there also appears to be a correlation between lows in the topography and lows in lubrication. However, that correlation is less compelling, especially recognizing that the inversion is less reliable near the boundaries of the solution region.
Figure 10 compares the observed (top row) and recovered (bottom row) surface topography (s) and velocity components (u, v). The spatial patterns are consistent. The amplitudes are recovered less well, especially for u and v. Figure 11 (top panel) gives a more detailed comparison of observed and recovered surface topography (s) along flow through the center of the solution region (y/h 0 = 0). The bottom panels of Figure 11 show the inferred bed topography (b) and lubrication (c) along the same line. The trough feature (Fig. 9, left panel) is displayed in diagonal section.
Input data (Fig. 7 with gradients removed and filtered (Equation (20)). Bottom row shows s, u and v recovered by applying forward theory to b and c inferred from inversion (Fig. 9). The surface topography s is in m, and the velocities are in m a−1.
Observed (thin line) and recovered (thick line) surface topography (s), bed topography (b), and lubrication (c/C0) along y = 0.
Figure 12 shows a transverse profile of the bed at x/h 0 ∼ 2:5. This is aligned with a radar profile giving a measured bed topography (Reference Bindschadler, Vornberger, Blankenship, Scambos and JacobelBindschadler and others, 1996, fig. 10b) plotted here as a thin solid curve. The dashed curve in Figure 12 also shows the bed deduced by Reference MacAyeal, Bindschadler and ScambosMacAyeal and others (1995, Fig. 6) for comparison. Neither inversion agrees completely with the radar measurements.The MacAyeal inversion gives a smoothed picture of the bed as it was designed to do. This inversion (thick solid line), following its design, captures more detail. It succeeds in identifying the deep, narrow trough evident in radar data, but there is an offset. Edge effects in the inversion may explain some of the offset. Co-registration errors in aligning the radar profile relative to the data for s, u and v may also contribute. However, the disagreement could arise because of nonuniqueness in the inversions discussed in section 3. Indeed, if the bed trough were aligned with the flow direction, it would produce a very small (or possibly zero) surface effect and would not be detectable with measurements of s, u and v.
Comparison of bed measured by radar as given by (Reference Bindschadler, Vornberger, Blankenship, Scambos and JacobelBindschadler and others (1996) (thin line), the inversion by (Reference MacAyeal, Bindschadler and ScambosMacAyeal and others (1995) (dashed line) and the inferred bed (thick line) near x/h 0 = 2:5.
Now we investigate how trustworthy this inversion is.We first examine sensitivity of the inversionto changes in Ω and Ξ (Table 2). It is encouraging that the best fits to the data are possible for combinations of Ω and Ξ that are close to the a priori values of 2.5 and 102 estimated as described above. That suggests the underlying physics of the model is reasonable. It also suggests that one might actually invert for Ω and Ξ. However, it is clear that the minimum in residuals is quite broad, and difficult to identify precisely. Inversions for different pairs of Ω and Ξ predict similar spatial features for b and c, but different amplitudes (higher amplitude on the bed as Ξ increases) and short scale variations (higher frequency on the bed if $\Omega$ increases). Therefore, if Ξ and Ω are not known, this inversion technique can retrieve the pattern of b and c but not the amplitudes.
The r2 values for various combinations of Ω and Ξ but fixed Σ‘s (see text)
Inversion for only the bed topography bor only basal lubrication c does not yield results as good as the combination (Table 1). The r 2 values for b and c alone are nevertheless still reasonable, and we cannot reject the corresponding inversions on this basis alone. The inversion for b only yields larger peaks on the bed to produce the high in the surface topography (near x/h 0 = 20), and downstream the topography looks more like the c pattern from the b and c inversion. The c-only pattern stays similar to the c pattern from the inversion for both b and c, although not with as strong a ridge-like feature, and added slippery and sticky spots to create the topography high near x/h 0 = 20. However, our a priori expectation is that there are both b and c variations under the ice stream, and on this basis we prefer the combined b and c inversion result.
We now examine whether the correlation between the low in the topography and low in lubrication is a robust feature of the inversion under the supposition that spatial variation of both b and c is present. Figure 13 shows that the residual r 2 increases substantially when the c(x, y) pattern is shifted by small Δx/h 0 and Δy/h 0 relative to the b(x, y) pattern, especially for shifts that are perpendicular to the lineations. Reversing the sign of the lubrication or topography also produces clear disagreement with the surface data. If perturbations in b and c are both present, then the lows in b and c are evidently nearly coincident.
Sensitivity of the weighted residual r2 to shifts in the c(x,y) pattern relative to b(x,y) by Δx and Δy.
The existence of low lubrication (stickyness) in the bed trough can have two basic sources: (1) a distinctive bed morphologically such as a rough, hard rock or consolidated sediments, or (2) low water pressure resulting in relatively high effective normal stress. There are several factors that argue against (2). A trough should be a location of relatively low upward conduction of heat through the thick, overlying ice. There should be relatively high basal heat generation from frictionwhere the bed is sticky.These factors would combine to produce higher basal melt rate in a sticky trough relative to its surroundings. The hydraulic gradient required to evacuate the extra meltwater would tend to produce high water pressure and low effective normal stress. The density difference between ice and water would also promote low effective normal stress in the trough. Finally, the bed low would tend to be a preferred flow path for water coming from upstream. These hydraulic effects associated with the low bed would tend to enhance lubrication. Therefore, we suggest that the first of the above explanations for low lubrication in the trough is more likely. We leave open the interesting question of why this sticky trough is there.
5. Conclusions
It is possible to invert measurements of surface topography and velocity for some features of bed topography and lubrication. Short-wavelength features (λ ≃ 2h 0) and longitudinal ridges/troughs and stripes in lubrication at all lateral scales cannot be recovered, leading to intrinsic non-uniqueness in the inversion even with precise input data. Application of the inversion method to available observation on parts of the surface of Ice Stream E identifies a lineation on the bed of the ice stream that runs diagonally across the flow direction.While the available surface data allow the possibility that the basal lineation is a topographic trough with no variation of lubrication or that it is a low-lubrication (sticky) zone with no topography, our preferred interpretation is a topographic trough with low lubrication in its base.
Acknowledgements
We wish to thank T. Jóhannesson who provided valuable guidance and suggestions as Scientific Editor, and reviewers R. Greve and T. Payne for their constructive comments. This work was supported by U.S. National Science Foundation grant No. OPP-9526707.
