1. Introduction

Ice from the interior of Antarctica is delivered to ice shelves and oceans through fast-flowing ice streams and glaciers. The ice streams flow up to two orders of magnitude faster than the surrounding ice and they are the most dynamic components of the ice-sheet system. These fast-flowing ice streams provide a connection between the floating ice shelves, which have relatively short-period responses, and the slower-moving inland ice sheet with a longer response time.

Global positioning system (GPS) receivers, which use data from satellites and do not need a visible reference point, are now commonly used in glaciological fieldwork. This method, with high temporal resolution and rigorous processing methods (Reference KingKing, 2004), has provided the breakthrough required to observe short-term changes in ice flow and to study ice dynamics on timescales that have not previously been accessible (e.g. Reference O’Neel, Echelmeyer and MotykaO’Neel and others, 2001; Reference DoakeDoake and others, 2002; Reference Anandakrishnan, Voigt, Alley and KingAnandakrishnan and others, 2003; Reference Bindschadler, Vornberger, King and PadmanBindschadler and others, 2003b; Reference Legrésy, Wendt, Tabacco, Rémy and DietrichLegrésy and others, 2004; Reference GudmundssonGudmundsson, 2006). In addition, the deployment of passive seismic recording stations that monitor microearthquakes has proven to be valuable in mapping ice-stream basal conditions (Reference Blankenship, Anandakrishnan, Kempf and BentleyBlankenship and others, 1987; Reference Anandakrishnan and BentleyAnandakrishnan and Bentley, 1993; Reference SmithSmith, 2006). Combining the two observation methods is a novel way to enhance the understanding of ice-stream dynamics and to relate what happens at the bed to the resulting flow regime.

It is now possible to observe highly variable flow dynamics and basal conditions on the ice streams in Antarctica and there appear to be some similarities, but also differences, between the ice streams that have been measured. On the Siple Coast, several ice-stream flow dynamics have been characterized. The speed of Bindschadler Ice Stream (former Ice Stream D) has been measured to vary by a factor of three over the course of a day (Reference Anandakrishnan, Voigt, Alley and KingAnandakrishnan and others, 2003). At the grounding line the ice stream speeds up during falling tide and slows during rising tide. This effect was observed to propagate upstream at 5.6 m s⁻¹ and can be detected at least 80 km upstream. On the ice plain at the mouth of Whillans Ice Stream (former Ice Stream B), Reference Bindschadler, Vornberger, King and PadmanBindschadler and others (2003b) showed a strong stick–slip motion behaviour, where the ice stream accelerates to speeds greater than 1 m h⁻¹ (30 times faster than normal) in as little as 30 s. This flow behaviour has been modelled with varying ocean back pressure, as the adjacent ice shelf flexes with the tides (Reference Bindschadler, King, Alley, Anandakrishnan and PadmanBindschadler and others, 2003a). On Kamb Ice Stream (former Ice Stream C), there is an anticorrelation between mean basal seismicity (averaged over many days) and mean flow speed (Reference Anandakrishnan and AlleyAnandakrishnan and Alley, 1997a). On shorter timescales, periods of basal seismic activity, at a site 10 km from the grounding line, correlate with low tide on the ice shelf. In this case, the tidal forcing travels up the ice stream at ∼1.6 m s⁻¹ (Reference Anandakrishnan and AlleyAnandakrishnan and Alley, 1997b). On the opposite side of West Antarctica, where ice flows into the Filchner–Ronne Ice Shelf, the surface velocity of Rutford Ice Stream has been measured to vary with a fortnightly period up to 40 km upstream from the grounding line; close to the grounding line the variation is about 20% (Reference GudmundssonGudmundsson, 2006). The longest available velocity record, of ∼2 years, reveals that this variability also occurs at annual and semi-annual periods with higher and more variable velocity around equinoxes (21 March and 22 September) compared to variations around solstices (21 June and 21 December), when most measurements are made (Reference Murray, Smith, King and WeedonMurray and others, 2007). These observations indicate that ice streams are sensitive to sea level with a short response time and could respond rapidly to changes in the mean sea level or any consequent changes in tidal range.

In this paper, surface velocities derived from GPS position measurements and passive seismic data from Rutford Ice Stream, 40 km upstream from the grounding line, are presented. The GPS measurements were recorded during a 6 week period in 2004/05, and the passive seismic data 7 years previously during a 4 week period in 1997. These data show a complex interaction between ocean tides, basal condition and ice-stream flow speed.

2. Location

Rutford Ice Stream is part of the West Antarctic ice sheet (WAIS). It flows in a south-southeast direction until it discharges into the Ronne Ice Shelf (Fig. 1). It is estimated that the area that feeds the Rutford Ice Stream is ∼49 000 km², which includes the northern part of the Ellsworth Mountains and borders the drainage basins of Pine Island Glacier, Carlson Inlet and Evans Ice Stream (Reference Doake, Alley and BindschadlerDoake and others, 2001). The thickness of the ice stream is ∼2 km, and the ice is grounded below sea level along its entire length. The surface elevation of the grounded areas is >30 m above that needed for hydrostatic equilibrium (Reference Doake, Frolich, Mantripp, Smith and VaughanDoake and others, 1987).

Fig. 1. Site map showing Rutford Ice Stream flowing into the Ronne Ice Shelf and the neighbouring Carlson Inlet. The five GPS antennas were located ∼40 km upstream from the grounding line, and the reference station is 31 km from the GPS array. The left inset shows the geometry of the measurement array. Diamonds indicate the geographic location of the GPS receivers, filled circles the passive seismic receivers that ran throughout the period shown in Figure 4 and open circles the ones that ran during only the first part of the experiment. The seismic and GPS receivers were not run simultaneously (see text for details).

Ground and airborne radio-echo sounding measurements on the ice stream have revealed variable ice thickness; the bed across the ice stream is W-shaped, with deep trenches running parallel to the flow (Reference Stephenson and DoakeStephenson and Doake, 1982). The location of the grounding line has been found by use of tiltmeters and flexure analysis (Reference Stephenson and DoakeStephenson and Doake, 1982; Reference StephensonStephenson, 1984; Reference SmithSmith, 1991) and synthetic aperture radar interferometry (Reference RignotRignot, 1998). The ice-stream bed at the grounding line is ∼1500 m below sea level and the basal topography is complicated. There is a pinning point on which the ice rests in the centre of the ice stream, which coincides with a topographic high in the bed (Reference Doake, Frolich, Mantripp, Smith and VaughanDoake and others, 1987), and interferograms reveal areas of partial grounding on the ice shelf (Reference RignotRignot, 1998).

Repeated seismic surveys on the ice stream, undertaken in the same area as the measurements presented here, have shown that there are both spatial (Reference SmithSmith, 1997a, b) and temporal (Reference SmithSmith and others, 2007) variations in the basal conditions. Despite a relatively steady ice-flow rate, rapid erosion, drumlin formation and changing hydrology have been observed over a short period of time (<6 years) (Reference SmithSmith and others, 2007).

3. GPS Data Processing

The GPS reference station, which was located in the Ellsworth Mountains (Fig. 1), 31 km away and 894 m above the field site, allowed a relative processing procedure to be used. The position of the reference station was determined using the Precise Point Positioning (PPP) approach implemented in the GIPSY/OASIS II software (Reference Zumberge, Heflin, Jefferson, Watkins and WebbZumberge and others, 1997). To ensure that systematic errors were not introduced during the data processing, rigorous processing strategies were adopted (Reference KingKing, 2004). The differential-phase kinematic positioning program, TRACK (part of the GAMIT software package (R.W. King and Y. Bock, unpublished information)), was used to process the continuous GPS data on the ice stream relative to the base station. The motion of each antenna was modelled as a random walk, with a constraint on the epoch-to-epoch motion set to 5–10 mm. GPS formal errors are generally over-optimistic, and hence we compute the uncertainties by examining the detrended coordinate time series computed between the base station and another linearly moving site ∼80 km away on the other side of Rutford Ice Stream. The root-mean-square of this time series suggests a precision of better than 1 cm (horizontal) and 3 cm (vertical). Given that this baseline is longer than the baselines to Rutford Ice Stream, we therefore expect the coordinate precisions of the velocity stakes to be better than 1 and 3 cm in the horizontal and vertical components, respectively. Visual inspection of the horizontal-position time series confirms this. An example of horizontal movement of stake D over 1 hour is shown in Figure 2.

Fig. 2. An example of the computed path of antenna D during 1 hour. Ice flow is from top to bottom. The estimated absolute accuracy in the horizontal direction is shown with error bars. A solution is computed every 10 s; larger crosses are drawn on the path every 15 min.

Inspection of the horizontal movement of all the stakes during the whole period (at 10 s measurement resolution) reveals that the motion is smooth, with no sign of stick–slip behaviour.

The average horizontal velocity of each antenna was computed by linear least squares, fitting a line through the antenna path. The resulting values are presented in Table 1. The flow velocity increases in the downstream direction and also southwest across the ice stream.

Table 1. From top to bottom: the along–flowline average velocity (m a⁻¹;estimated error 0.06 m a⁻¹),the direction of the average velocity (east from north; estimated error 0.01 °), the standard deviation from mean velocity (m a⁻¹), the percentage of the variance in the velocity record that can be explained by tidal components and the amplitude (m a⁻¹) of the six largest tidal constituents resulting from harmonic analysis of the velocity records. These values are plotted against the corresponding amplitudes of the tidal signal on the ice shelf in Figure 6

The coordinate system for each antenna was then rotated to have the x axis along the mean flow direction, and the time series of position was then detrended by subtracting the mean linear motion. Neither the surface elevation nor the lateral horizontal movement (i.e. perpendicular to the mean flow direction) shows anything other than smooth motion during the recording period. In contrast, the detrended, along-flow horizontal motion (Fig. 3) shows clear, periodic flow variations: a dominant fortnightly variation, with diurnal and semi-diurnal frequencies superimposed upon it. This flow variation is analogous to similar GPS data of horizontal along-flow positions from Rutford Ice Stream that were collected at the grounding line, 10, 20 and 40 km upstream from the grounding line (Reference GudmundssonGudmundsson, 2006). The new data therefore confirm that tides have a persistent influence on the flow of the ice stream at least as far as 40 km upstream from the grounding line. Below we analyze a high-resolution velocity time series derived from these along-flow horizontal motion records.

Fig. 3. The detrended inline position of the five antennas showing the fortnightly flow variability of the ice stream. Deviation from mean longitudinal motion is shown: negative values indicate that the antenna has been moving slower than the mean velocity. Top to bottom: U, NE, C, SW, D. For clarity, curves for NE and SW are shifted +0.2 and −0.2 m, respectively, and lines for U and D are shifted +0.4 and −0.4 m, respectively. The inset in the top left corner shows 24 hours of data for SW during 10 January 2005.

4. Passive Seismic Data

The characteristics and interpretation of the passive seismic data were presented in detail by Reference SmithSmith (2006), including examples of typical waveforms (see fig. 2 therein). In summary, events from the bed were identified from waveform characteristics (Reference Blankenship, Anandakrishnan, Kempf and BentleyBlankenship and others, 1987; Reference Anandakrishnan and BentleyAnandakrishnan and Bentley, 1993) and from the difference in their P- and S-wave travel times, knowing the ice thickness from seismic reflection surveys and typical seismic velocities in ice (∼3600–3900 m s⁻¹ for P waves and ∼1700–1950 m s⁻¹ for S waves). For each bed event, this difference in travel times was consistent with P and S waves originating from a single event at the ice-stream bed. Events from all other origins were discarded. No bed events could be identified at more than one recording station, and the maximum observed difference in P- and S-wave travel times indicated a source with an epicentral distance of no greater than 1 km from the recording station. Hence, the bed events are low-energy and each station was sensitive only to those occurring within an epicentral radius of ∼1 km. This means that all the basal events detected at a given recording station originate within a relatively small area of the bed (radius 1 km) centred directly beneath that station. The hourly mean number of basal events was calculated across the whole array (Fig. 4). The seismic events from the ice-stream bed all have similar amplitudes (Reference SmithSmith, 2006), and hence, although the coarse sample rate (2 ms) precludes a rigorous quantitative analysis of wavelet energies, the number of events detected is still a reasonable indication of the relative levels of elastic energy being released at the bed.

Fig. 4. Mean basal seismicity across the array, showing the number of events per hour recorded at all active stations (ten stations during the first period (30 November 1997 to 12 December 1997), five stations in the second period (15–23 December 1997)).

Mean basal seismicity across the array (Fig. 4) shows active periods of a few hours, separated by inactive periods of between 1 hour and 1 day. Basal seismicity varies across the array (Reference SmithSmith, 2006), so the reduced number of stations from 15 December 1997 onwards (five rather than ten) means the basal seismicity magnitude cannot necessarily be compared directly between the two recording periods in Figure 4. However, although the level of seismicity varies, the timing of any relatively active and inactive periods of basal seismicity should be largely independent of the number of stations or their locations. Hence, we concentrate primarily on the temporal aspects of the passive seismic record, i.e. whether a particular time is relatively active or not.

5. Tidal Record at the Grounding Line

To allow an assessment of the influence of ocean tides on the observed ice-stream velocity and basal seismicity variations, we generated tidal predictions of the surface elevation a short distance downstream of the grounding line, where the ice-shelf vertical displacement is driven by tides in the ocean beneath. The highest tidal amplitudes around Antarctica are found in the southern Filchner–Ronne Ice Shelf area, south of Henry and Korff ice rises. The spring-tide peak-to-peak range exceeds 7 m in the mouth of Rutford Ice Stream, and the major semi-diurnal constituents M2 and S2 dominate tidal records (Reference StephensonStephenson, 1984; Reference DoakeDoake, 1992; Reference Padman, Fricker, Coleman, Howard and ErofeevaPadman and others, 2002). A 56 day GPS data series of the vertical motion close to the grounding line (personal communication from G.H. Gudmundsson, 2005) has been analyzed with standard tidal analysis methods using the T-TIDE software package (Reference Pawlowicz, Beardsley and LentzPawlowicz and others, 2002). The four main constituents are the O1 and K1 (diurnal) and M2 and S2 (semi-diurnal) constituents. Two other constituents are significant, Q1 and N2, and two constituents, K2 and P1, were inferred from S2 and K1, respectively (Reference PughPugh, 1987), using amplitude and phase relationships from the CATS02.01 tide model for Antarctic ice shelves and seas (Reference Padman, Fricker, Coleman, Howard and ErofeevaPadman and others, 2002). The amplitude and phase of fortnightly and monthly tidal constituents could not be resolved from the original GPS record, so these were also obtained from the CATS02.01 tide model. Tidal predictions for the periods of the GPS and passive seismic experiments were constructed from these ten constituents. The resulting tidal prediction for the period of the GPS experiment is shown in Figure 5a.

Fig. 5. (a) Tidal signal constructed with eight diurnal and semi-diurnal constituents and the lunar fortnightly and monthly (Mf and Mm) constituents for the period when the GPS measurements on Rutford Ice Stream were obtained. (b) The velocity variation from the mean velocity computed from the positions shown in Figure 3 as explained in the text. Top to bottom: U, NE, C, SW, D. For clarity, curves for NE and SW are shifted +150 and −150 m a⁻¹,respectively, and curves for U and D are shifted +300 and −300 m a⁻¹.The inset in the bottom left corner shows the computed velocity of SW during 10 January 2005, the same period as shown in the inset of Figure 3.

6. GPS-Data Velocity Analysis

6.3. Harmonic analysis for tidal species

The general frequency analysis described above indicates that the GPS data are dominated by signals with typical ocean-tide frequency values. Hence, to allow a more rigorous comparison between the GPS data and the tidal record from the ice shelf, the GPS velocity records were analyzed specifically for the semi-diurnal, diurnal and fortnightly constituents present in the ice-shelf record using the T-TIDE software package (Reference Pawlowicz, Beardsley and LentzPawlowicz and others, 2002). The data-sets are too short to include the longer-period components such as monthly, semi-annual and annual, although analysis of a longer time series at this location demonstrates that these components are present (Reference Murray, Smith, King and WeedonMurray and others, 2007). For consistency, the tidal prediction was also analyzed in the same way. Although the constituents used to generate the prediction are known, this approach means that uncertainties due to the limited record lengths, particularly in the longer-period constituents, should be similar for both the GPS velocity analysis and that of the tidal prediction. This allows a rigorous comparison to be made between the two datasets. It also gives some assessment of the proportion of the GPS record that can be explained by the anticipated ocean forcing, and how much must be attributed to other mechanisms, or to unresolved tidal components.

Two fortnightly constituents were specified in the harmonic analysis: MSf (14.76 day period), the lunisolar synodic fortnightly constituent; and Mf, the lunar fortnightly constituent (period 13.66 days). The main frequency in the observed time series is closer to 15 days. Globally, Mf is typically the dominant fortnightly constituent, due to direct gravitational forcing. However, due to non-linear interaction between M2 and S2 in shallow water (Reference PughPugh, 1987), MSf may be the dominant fortnightly term around much of Antarctica, particularly in the sub-ice-shelf regions. There are, at present, too few GPS or tilt data on the Ronne Ice Shelf near the grounding zone to separate Mf and MSf reliably and hence determine their respective magnitudes.

Three separate harmonic analyses were done using different fortnightly components: with MSf only, with Mf only and finally with both. The best result (giving the highest signal-to-noise ratio) was obtained when both constituents were used in the analysis. In this case, MSf was about two times larger in amplitude than Mf. The length of our velocity record is, however, only two fortnightly cycles, which is rather short to determine fortnightly characteristics accurately. Similar analysis of a 2 year velocity record nearby (Reference Murray, Smith, King and WeedonMurray and others, 2007) showed an amplitude of MSf about five times larger than Mf. Adding the two fortnightly components in our analysis gives the total contribution at the period, which is ∼10–12% of the total velocity. This is in good agreement with similar data from Reference GudmundssonGudmundsson (2006), who analyzed the in-line position.

The harmonic analysis shows that the semi-diurnal, diurnal and the two fortnightly tidal components can account for up to 86% of the variance in the observed velocity record (Table 1). The remaining variance is due to limitations of this analysis due to the short length of the record, to other tides (including higher and lower frequencies) and noise and/or potential non-tidal contributions. The amplitudes of the constituents at each station are plotted against the corresponding constituent in the ice-shelf tidal signal in Figure 6, and the amplitudes of the six largest constituents are given in Table 1. The amplitude at fortnightly periods decreases upstream, with the amplitude at D and C of similar value but lower at U. The semi-diurnal constituents (M2, S2 and K2) have larger amplitudes than the diurnal ones, ∼4–12% of the total velocity. The diurnal components (O1, P1 and K1) have ∼3–5% amplitude compared to the total velocity. The estimates for the phase of each constituent have too high an error for conclusions to be drawn about the change in phase across or along the ice stream over the stake array.

Fig. 6. The amplitude of each tidal component of the velocity records on the ice stream plotted against the corresponding amplitude on the ice shelf. Formal errors are shown as error bars at the 95% confidence level.

Figure 6 shows that for constituents with small amplitude on the ice shelf (Q1, P1, N2 and K2), the ice-stream response is smaller than for those of larger amplitude (O1, K1, S2 and M2). Three further details of Figure 6 are also worth noting:

1. 1. The response of the ice stream at the fortnightly periods is dominant, despite the relatively small vertical tidal forcing of that period on the ice shelf.

2. 2. For the two largest constituents in the ice-shelf signal (semi-diurnal S2 and M2) the response of the ice stream is of similar magnitude, even though the amplitude of M2 is more than 1.5 times that of S2 on the ice shelf.

3. 3. For most tidal frequencies, the errors in the harmonic analysis of the velocity data are such that the ice-stream response at each station cannot be claimed to differ significantly from that at any other site.

These observations have a number of implications. In the case of the fortnightly periods, which have a small amplitude on the ice shelf, the ice stream cannot solely be responding to that forcing but rather to the tidal range, which also has a period of ∼15 days (two per lunar month of 29.5 days, or 14.75 days). Compare Figure 5a and b and observe that the velocity is high when the tidal range is large and low when the tidal range is smaller. Both autocorrelation of each of the velocity series, which indicates a periodicity of 15 days, and power spectra of FFT analysis, which shows by far the most energy in the 15 day period, support this. The second observation indicates that the amplitude of the forcing at the ice shelf is not the only factor controlling the ice-stream response. This has also been observed for the longer than fortnightly periods; semi-annual response of the ice stream does not scale with the forcing on the ice shelf (Reference Murray, Smith, King and WeedonMurray and others, 2007). The third observation indicates that no significant difference in the ice-stream response to ocean tide forcing can be detected across an array of these dimensions (6 km × 6 km) using the data acquisition, processing and harmonic analysis approaches we have employed.

7. Tidal Correlation with Surface Velocity and Basal Seismicity

To compare the ice-flow variability and basal seismicity with the tidal record at the ice shelf, the time series were plotted together (Fig. 7). The ocean tide record is presented as the magnitude of the deviation from mean tide height (negative values of the ocean tide record (low tides in Fig. 5a) are drawn as positive values with a grey curve). It must be kept in mind that the tidal record is based on analysis of data from a different period and is not error-free. An example of the velocity record from one station during a 12 day period from spring tide to neap tide, as well as the first part of the seismic record on the same timescale are shown in Figure 7. It should also be kept in mind that the velocity perturbation propagates upstream with a velocity that can be variable along the ice stream. The cross-correlation of the position records indicates a propagation velocity of 10 ± 4 m s⁻¹, which would separate the velocity and the tidal signal by 0.8–1.8 hours over the 40 km distance between the grounding line (assumed to be the location of the forcing) and the location of the measurements. Similar measurements over a larger area of the ice stream, from the grounding line to the measurement site, by Reference GudmundssonGudmundsson (2006) indicate a propagation velocity in the range 1–2 m s⁻¹,which separates the two records by 5–11 hours. The unknown separation of the records in Figure 7 makes it impossible to be certain whether a velocity or seismic peak occurs at high, falling, low or rising tide. However, comparing the two records without separation gives information about the occurrence of velocity or seismic peaks at equivalent times within the tidal cycle, as the separation can be assumed to be the same throughout the records.

Fig. 7. (a) An example 12 day velocity record from C, drawn together with the modelled tidal record on the ice shelf. (b) The mean basal seismicity (events per hour) plotted on top of the modelled tidal record on the ice shelf. Only the first 12 days, while ten receivers were running, are shown. The predicted tidal record in both panels is shown as positive deviation from the mean. The high tide is plotted with a thick black curve, the low tide with a grey curve.

Inspection of the entire record for all five stations, as well as the seismic record, reveals that there is no simple relationship between the ocean tide and the velocity and basal seismicity of the ice stream. In the velocity record there is, however, a distinct pattern visible, which is characterized by two velocity peaks during one tidal cycle (61–72% of the tidal cycles have two velocity peaks). During about a third of the observed tidal cycles, there were three velocity peaks. For the basal seismicity record there are similarly two peaks in seismicity during 67% of the tidal cycles. The variability in seismic activity is larger than that in the velocity record.

Dividing the semi-diurnal tidal cycle into four equally long parts and counting the occurrence of peaks within each part gives information about relative timing of velocity or seismicity peaks within the tidal cycle. For the velocity record, simple peak counting reveals that about three-quarters (68–83%) of the first part of the tidal cycle coincides with velocity peaks, there are velocity peaks in 29–49% of the next part, 55–64% of the third part and 9–17% of the last part (the percentage ranges represent the values from the five stations). The five stations do not have the same response (e.g. the occurrence of three velocity peaks or a missing double peak does not happen concurrently at all sites). The exceptions to the general behaviour (missing a double peak or three velocity peaks) occur randomly over the array, with more exceptions during neap tide than spring tide. The pattern between the velocity and the tidal record at the ice shelf is stronger during spring tide than neap tide. Similar counting of peaks in seismic activity within the four equally long parts of the tidal cycle shows that during 35 tidal cycles there are 24 peaks in the first part, 11 in the second, 24 in the third and 4 in the fourth part. The record of basal seismicity is not long enough, nor sufficiently complete, to determine with certainty whether the level of basal seismicity is greater during periods of spring tides. However, we note that the correlation with tides is still seen on some occasions during neap tide periods, as shown in the velocity record. Longer seismic records would be required to investigate this further.

The comparison of the velocity and seismic records with the tidal record at the ice shelf shows that both velocity and seismicity increase at the same time within the tidal cycle. The first and third parts of the semi-diurnal tidal cycle have highest velocity and seismicity, and the fourth part the lowest. The timing is, however, not precise, as can be seen in Figure 7.

8. Discussion

The observed flow of Rutford Ice Stream is strongly affected by the tides downstream of the measurement site. The fourth row in Table 1 shows that the semi-diurnal, diurnal and fortnightly ocean tides can account for up to 86% of the velocity variations evident in the GPS data. Comparison of rows three and four in Table 1 shows that the largest variation in velocity is at D and this decreases upstream, but the tidal components can explain the largest percentage of the variance in the velocity at U. This indicates that contributions of noise, non-tidal disturbance and higher or lower tidal periods to the total velocity variance are not constant between the sites. The unexplained variance may well be the result of a combination of factors. Longer-period (e.g. monthly, semi-annual and annual) tidal components have been observed in a similar velocity record from the same location (Reference Murray, Smith, King and WeedonMurray and others, 2007), but the records presented here are too short to analyze for them, and nonlinear forcing at grounding lines can lead to significant higher-frequency signals (Reference Pedley, Paren and PotterPedley and others, 1986), for which we have also not analyzed.

The analysis of the 6 week high-resolution and continuous-velocity records during the 2004/05 austral summer confirms the conclusion of Reference Murray, Smith, King and WeedonMurray and others (2007) that the magnitude of the downstream velocity response at each frequency does not scale linearly with the amplitude of the corresponding vertical tidal forcing. Even though the measurements are made during solstice when the ice flow is slowest and the semi-diurnal variation is smallest (Reference Murray, Smith, King and WeedonMurray and others, 2007), there are two or three velocity peaks observed during each tidal cycle.

The analysis of the velocity and seismic records from Rutford Ice Stream reveals some similarities and some differences between this ice stream and the Siple Coast ice streams. The position records, which have high temporal resolution (10 s), exclude any stick–slip motion on Rutford Ice Stream. Increased velocity is observed during all parts of the tidal cycle, and therefore also during some of the falling tides, indicating that the same effect as on Bindschadler Ice Stream is acting (Reference Anandakrishnan, Voigt, Alley and KingAnandakrishnan and others, 2003). In contrast to observations on Kamb Ice Stream, where there is anticorrelation between basal seismicity and mean flow speed and, on shorter timescales, correlation of basal seismicity and low tide on the ice shelf (Reference Anandakrishnan and AlleyAnandakrishnan and Alley, 1997a, b), we observe basal seismicity correlating with velocity and highest seismic activity during spring tide when the velocity is highest.

Rather than attempting to model all the different mechanisms that appear to be interacting simultaneously on Rutford Ice Stream, we make the following list which can contribute to the understanding of the transfer of the tidal force at the ice shelf to the flow of the ice stream, in particular, the increase in both velocity and seismic activity two (and sometimes three) times during each tidal cycle.

9. Conclusion

During the austral summer, Rutford Ice Stream flows at a speed of ∼370 m a⁻¹ (367–380 m a⁻¹),40 km upstream from the grounding line. Comparison between five velocity time series and the passive seismic record and the tidal record indicates a complicated response to the tidal forcing. Both the surface velocity and basal seismicity show two peaks, and sometimes three, during a tidal cycle, indicating that ocean tides at the grounding line modulate the ice-stream flow. The transfer of the ocean forcing up the ice stream is complicated and the effects of irregular grounding zone with areas of partial grounding (Reference RignotRignot, 1998) and both spatially and temporally variable basal conditions (Reference SmithSmith, 1997a, Reference Smithb; Reference SmithSmith and others, 2007) contribute to the complexity. The variation in the velocity with fortnightly period has an amplitude of ∼10–12% of the total velocity, but at diurnal and semi-diurnal period the variation amplitude is 3–5% and 4–12%, respectively. We summarize our observations as follows:

These observations suggest that numerous mechanisms are acting simultaneously. It is clear that the ice stream is sensitive to sea-level variations, at least at tidal frequencies. This sensitivity could have global implications, as the observed sea-level rise could be accelerated by increased ice-stream discharge.